Damien Irving edited methods_id.tex  over 8 years ago

Commit id: 90d457334a77cf5e31d270ed41afdb152364d2aa

deletions | additions      

       

Given that the PSA pattern completes approximately 1.6 to 2.0 cycles (depending on the specific EOF mode) over the 120$^{\circ}$ search area (see Figure \ref{fig:eof}), our analysis focused data times where a Fourier transform revealed wavenumber 5 and 6 as dominant frequencies over the zero padded 360$^{\circ}$ domain. In particular, a data time was said to display PSA-like variability (and hence was selected for further analysis) if the amplitude of the wavenumber 5 and 6 components of the Fourier transform were ranked in the top three of all frequencies. The vague `PSA-like' descriptor is used because a number of features besides the PSA pattern (e.g. Antarctic Dipole, Amundsen Sea Low, zonal wave three pattern) can exhibit wavenumber 5-6 variability in the PSA sector.  Once these data times were selected, additional information from the Fourier transform was used to characterize the phase and amplitude of the PSA-like variability. With respect to the former, it can be seen from Figure \ref{fig:transform} that within the search area the phase of the wavenumber 5 and 6 components of the transform (and usually also adjacent frequencies like wavenumber 4 and 7) tend to align both with each other and also with the phase of the actual signal. The phase of the wavenumber 6 component of the Fourier transform was therefore used as a proxy for the phase of the signal as a whole, and this information was used to separate data times displaying the actual PSA pattern from the larger population of PSA-like variability (note that similar results were obtained using wavenumber 5). The details of this separation process (e.g. the phase ranges used to define the PSA pattern) are discussed in Section \ref{s:psa_results}. below.  In order to quantify the amplitude of PSA-like variability, the wave envelope construct pioneered (in the atmospheric sciences) by \citet{Zimin2003} and recently applied by \citet{IrvingSimmonds2015} was used. The envelope is obtained by performing a Fourier transform, followed by an inverse Fourier transform for only the wavenumbers of interest. The complex number amplitude of the resulting waveform represents the envelope. Since the envelope of the complete signal (i.e. with all wavenumbers retained) can be quite noisy, the amplitude of PSA-like variability was defined as the maximum value of the envelope when only wavenumbers 4 to 7 are retained (see Figure \ref{fig:transform} for an example envelope). \subsubsection{Timescale considerations}