• Ashok, K., , S. K. Behera, , S. A. Rao, , H. Weng, , and T. Yamagata, 2007: El Niño Modoki and its possible teleconnection. J. Geophys. Res., 112, C11007, doi:10.1029/2006JC003798.

    • Search Google Scholar
    • Export Citation
  • Bonaventura, L., , R. Redler, , and R. Budich, 2012: Algorithms, Code Infrastructure and Optimisation. Vol. 2, Earth System Modelling, Springer, 96 pp.

  • Boschat, G., , I. Simmonds, , A. Purich, , T. Cowan, , and A. B. Pezza, 2016: On the use of composite analyses to form physical hypotheses: An example from heat wave–SST associations. Sci. Rep., 6, 29599, doi:10.1038/srep29599.

    • Search Google Scholar
    • Export Citation
  • Bracegirdle, T. J., , and G. J. Marshall, 2012: The reliability of Antarctic tropospheric pressure and temperature in the latest global reanalyses. J. Climate, 25, 71387146, doi:10.1175/JCLI-D-11-00685.1.

    • Search Google Scholar
    • Export Citation
  • Bromwich, D. H., , J. P. Nicolas, , and A. J. Monaghan, 2011: An assessment of precipitation changes over Antarctica and the Southern Ocean since 1989 in contemporary global reanalyses. J. Climate, 24, 41894209, doi:10.1175/2011JCLI4074.1.

    • Search Google Scholar
    • Export Citation
  • Christoph, M., , T. P. Barnett, , and E. Roeckner, 1998: The Antarctic Circumpolar Wave in a coupled ocean–atmosphere GCM. J. Climate, 11, 16591672, doi:10.1175/1520-0442(1998)011<1659:TACWIA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ciasto, L. M., , G. R. Simpkins, , and M. H. England, 2015: Teleconnections between tropical Pacific SST anomalies and extratropical Southern Hemisphere climate. J. Climate, 28, 5665, doi:10.1175/JCLI-D-14-00438.1.

    • Search Google Scholar
    • Export Citation
  • Clem, K. R., , and R. L. Fogt, 2015: South Pacific circulation changes and their connection to the tropics and regional Antarctic warming in austral spring, 1979–2012. J. Geophys. Res. Atmos., 120, 27732792, doi:10.1002/2014JD022940.

    • Search Google Scholar
    • Export Citation
  • Clem, K. R., , and J. A. Renwick, 2015: Austral spring Southern Hemisphere circulation and temperature changes and links to the SPCZ. J. Climate, 28, 73717384, doi:10.1175/JCLI-D-15-0125.1.

    • Search Google Scholar
    • Export Citation
  • Dawson, A., 2016a: eofs: A library for EOF analysis of meteorological, oceanographic and climate data. J. Open Res. Software, 4, e14, doi:10.5334/jors.122.

    • Search Google Scholar
    • Export Citation
  • Dawson, A., 2016b: Windspharm: A High-Level Library for Global Wind Field Computations Using Spherical Harmonics. J. Open Res. Software, 4, e31, doi:10.5334/jors.129.

    • Search Google Scholar
    • Export Citation
  • Dee, D. P., and et al. , 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553597, doi:10.1002/qj.828.

    • Search Google Scholar
    • Export Citation
  • Dee, D. P., , M. Balmaseda, , G. Balsamo, , R. Engelen, , A. J. Simmons, , and J.-N. Thépaut, 2014: Toward a consistent reanalysis of the climate system. Bull. Amer. Meteor. Soc., 95, 12351248, doi:10.1175/BAMS-D-13-00043.1.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , and E. J. Steig, 2013: Temperature change on the Antarctic Peninsula linked to the tropical Pacific. J. Climate, 26, 75707585, doi:10.1175/JCLI-D-12-00729.1.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , E. J. Steig, , D. S. Battisti, , and M. Küttel, 2011: Winter warming in West Antarctica caused by central tropical Pacific warming. Nat. Geosci., 4, 398403, doi:10.1038/ngeo1129.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , E. J. Steig, , D. S. Battisti, , and J. M. Wallace, 2012: Influence of the tropics on the southern annular mode. J. Climate, 25, 63306348, doi:10.1175/JCLI-D-11-00523.1.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , J. M. Wallace, , D. S. Battisti, , E. J. Steig, , A. J. E. Gallant, , H.-J. Kim, , and L. Geng, 2014: Tropical forcing of the recent rapid Arctic warming in northeastern Canada and Greenland. Nature, 509, 209212, doi:10.1038/nature13260.

    • Search Google Scholar
    • Export Citation
  • Fogt, R. L., , and D. H. Bromwich, 2006: Decadal variability of the ENSO teleconnection to the high-latitude South Pacific governed by coupling with the southern annular mode. J. Climate, 19, 979997, doi:10.1175/JCLI3671.1.

    • Search Google Scholar
    • Export Citation
  • Fogt, R. L., , and A. J. Wovrosh, 2015: The relative influence of tropical sea surface temperatures and radiative forcing on the Amundsen Sea low. J. Climate, 28, 85408555, doi:10.1175/JCLI-D-15-0091.1.

    • Search Google Scholar
    • Export Citation
  • Fogt, R. L., , J. M. Jones, , and J. Renwick, 2012: Seasonal zonal asymmetries in the southern annular mode and their impact on regional temperature anomalies. J. Climate, 25, 62536270, doi:10.1175/JCLI-D-11-00474.1.

    • Search Google Scholar
    • Export Citation
  • Ghil, M., , and K. Mo, 1991: Intraseasonal oscillations in the global atmosphere. Part II: Southern Hemisphere. J. Atmos. Sci., 48, 780790, doi:10.1175/1520-0469(1991)048<0780:IOITGA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Glatt, I., , and V. Wirth, 2014: Identifying Rossby wave trains and quantifying their properties. Quart. J. Roy. Meteor. Soc., 140, 384396, doi:10.1002/qj.2139.

    • Search Google Scholar
    • Export Citation
  • Gong, D., , and S. Wang, 1999: Definition of Antarctic Oscillation index. Geophys. Res. Lett., 26, 459462, doi:10.1029/1999GL900003.

  • Grimm, A. M., , and T. Ambrizzi, 2009: Teleconnections into South America from the tropics and extratropics on interannual and intraseasonal timescales. Past Climate Variability in South America and Surrounding Regions, F. Vimeux, F. Sylvestre, and M. Khodri, Eds., Springer, 159–191, doi:10.1007/978-90-481-2672-9_7.

  • Harangozo, S. A., 2004: The relationship of Pacific deep tropical convection to the winter and springtime extratropical atmospheric circulation of the South Pacific in El Niño events. Geophys. Res. Lett., 31, L05206, doi:10.1029/2003GL018667.

    • Search Google Scholar
    • Export Citation
  • Hoskins, B. J., , and D. J. Karoly, 1981: The steady linear response of a spherical atmosphere to thermal and orographic forcing. J. Atmos. Sci., 38, 11791196, doi:10.1175/1520-0469(1981)038<1179:TSLROA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Hunter, J. D., 2007: Matplotlib: A 2D graphics environment. Comput. Sci. Eng., 9, 9095, doi:10.1109/MCSE.2007.55.

  • Irving, D. B., 2016a: A minimum standard for publishing computational results in the weather and climate sciences. Bull. Amer. Meteor. Soc., doi:10.1175/BAMS-D-15-00010.1, in press.

    • Search Google Scholar
    • Export Citation
  • Irving, D. B., 2016b: A new method for identifying the Pacific–South American pattern and its influence on regional climate variability: Supplementary metadata. Figshare, digital media, doi:10.6084/m9.figshare.1612460.

  • Irving, D. B., , and I. Simmonds, 2015: A novel approach to diagnosing Southern Hemisphere planetary wave activity and its influence on regional climate variability. J. Climate, 28, 90419057, doi:10.1175/JCLI-D-15-0287.1.

    • Search Google Scholar
    • Export Citation
  • Karoly, D. J., 1989: Southern Hemisphere circulation features associated with El Niño–Southern Oscillation events. J. Climate, 2, 12391252, doi:10.1175/1520-0442(1989)002<1239:SHCFAW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Kidson, J. W., 1988: Interannual variations in the Southern Hemisphere circulation. J. Climate, 1, 11771198, doi:10.1175/1520-0442(1988)001<1177:IVITSH>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Lachlan-Cope, T., , and W. Connolley, 2006: Teleconnections between the tropical Pacific and the Amundsen-Bellinghausens Sea: Role of the El Niño/Southern Oscillation. J. Geophys. Res., 111, D23101, doi:10.1029/2005JD006386.

    • Search Google Scholar
    • Export Citation
  • Lau, K.-M., , P.-J. Sheu, , and I.-S. Kang, 1994: Multiscale low-frequency circulation modes in the global atmosphere. J. Atmos. Sci., 51, 11691193, doi:10.1175/1520-0469(1994)051<1169:MLFCMI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Li, X., , D. M. Holland, , E. P. Gerber, , and C. Yoo, 2014: Impacts of the north and tropical Atlantic Ocean on the Antarctic Peninsula and sea ice. Nature, 505, 538542, doi:10.1038/nature12945.

    • Search Google Scholar
    • Export Citation
  • Li, X., , D. M. Holland, , E. P. Gerber, , and C. Yoo, 2015: Rossby waves mediate impacts of tropical oceans on West Antarctic atmospheric circulation in austral winter. J. Climate, 28, 81518164, doi:10.1175/JCLI-D-15-0113.1.

    • Search Google Scholar
    • Export Citation
  • Li, Y., , J. Li, , F. F. Jin, , and S. Zhao, 2015: Interhemispheric propagation of stationary Rossby waves in a horizontally nonuniform background flow. J. Atmos. Sci., 72, 32333256, doi:10.1175/JAS-D-14-0239.1.

    • Search Google Scholar
    • Export Citation
  • Liu, Z., , and M. Alexander, 2007: Atmospheric bridge, oceanic tunnel, and global climatic teleconnections. Rev. Geophys., 45, RG2005, doi:10.1029/2005RG000172.

    • Search Google Scholar
    • Export Citation
  • Liu, Z., , Z. Jian, , K. Yoshimura, , N. H. Buenning, , C. J. Poulsen, , and G. J. Bowen, 2015: Recent contrasting winter temperature changes over North America linked to enhanced positive Pacific–North American pattern. Geophys. Res. Lett., 42, 77507757, doi:10.1002/2015GL065656.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., 2000: Relationships between low-frequency variability in the Southern Hemisphere and sea surface temperature anomalies. J. Climate, 13, 35993610, doi:10.1175/1520-0442(2000)013<3599:RBLFVI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., , and M. Ghil, 1987: Statistics and dynamics of persistent anomalies. J. Atmos. Sci., 44, 877902, doi:10.1175/1520-0469(1987)044<0877:SADOPA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., , and R. W. Higgins, 1998: The Pacific–South American modes and tropical convection during the Southern Hemisphere winter. Mon. Wea. Rev., 126, 15811596, doi:10.1175/1520-0493(1998)126<1581:TPSAMA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., , and J. N. Paegle, 2001: The Pacific–South American modes and their downstream effects. Int. J. Climatol., 21, 12111229, doi:10.1002/joc.685.

    • Search Google Scholar
    • Export Citation
  • Nicolas, J. P., , and D. H. Bromwich, 2011: Precipitation changes in high southern latitudes from global reanalyses: A cautionary tale. Surv. Geophys., 32, 475494, doi:10.1007/s10712-011-9114-6.

    • Search Google Scholar
    • Export Citation
  • Nicolas, J. P., , and D. H. Bromwich, 2014: New reconstruction of Antarctic near-surface temperatures: Multidecadal trends and reliability of global reanalyses. J. Climate, 27, 80708093, doi:10.1175/JCLI-D-13-00733.1.

    • Search Google Scholar
    • Export Citation
  • North, G. R., , T. L. Bell, , and R. F. Cahalan, 1982: Sampling errors in the estimation of empirical orthogonal functions. Mon. Wea. Rev., 110, 699706, doi:10.1175/1520-0493(1982)110<0699:SEITEO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Notaro, M., , W.-C. Wang, , and W. Gong, 2006: Model and observational analysis of the northeast U.S. regional climate and its relationship to the PNA and NAO patterns during early winter. Mon. Wea. Rev., 134, 34793505, doi:10.1175/MWR3234.1.

    • Search Google Scholar
    • Export Citation
  • Raphael, M. N., 2004: A zonal wave 3 index for the Southern Hemisphere. Geophys. Res. Lett., 31, L23212, doi:10.1029/2004GL020365.

  • Raphael, M. N., , and W. R. Hobbs, 2014: The influence of the large-scale atmospheric circulation on Antarctic sea ice during ice advance and retreat seasons. Geophys. Res. Lett., 41, 50375045, doi:10.1002/2014GL060365.

    • Search Google Scholar
    • Export Citation
  • Renwick, J. A., , and M. J. Revell, 1999: Blocking over the South Pacific and Rossby wave propagation. Mon. Wea. Rev., 127, 22332247, doi:10.1175/1520-0493(1999)127<2233:BOTSPA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Robertson, A. W., , and C. R. Mechoso, 2003: Circulation regimes and low-frequency oscillations in the South Pacific sector. Mon. Wea. Rev., 131, 15661576, doi:10.1175//2548.1.

    • Search Google Scholar
    • Export Citation
  • Sardeshmukh, P. D., , and B. J. Hoskins, 1988: The generation of global rotational flow by steady idealized tropical divergence. J. Atmos. Sci., 45, 12281251, doi:10.1175/1520-0469(1988)045<1228:TGOGRF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Schneider, D. P., , C. Deser, , and Y. Okumura, 2012: An assessment and interpretation of the observed warming of West Antarctica in the austral spring. Climate Dyn., 38, 323347, doi:10.1007/s00382-010-0985-x.

    • Search Google Scholar
    • Export Citation
  • Screen, J. A., , and I. Simmonds, 2012: Half-century air temperature change above Antarctica: Observed trends and spatial reconstructions. J. Geophys. Res., 117, D16108, doi:10.1029/2012JD017885.

    • Search Google Scholar
    • Export Citation
  • Simmonds, I., 2015: Comparing and contrasting the behaviour of Arctic and Antarctic sea ice over the 35 year period 1979–2013. Ann. Glaciol., 56, 1828, doi:10.3189/2015AoG69A909.

    • Search Google Scholar
    • Export Citation
  • Simpkins, G. R., , S. McGregor, , A. S. Taschetto, , L. M. Ciasto, , and M. H. England, 2014: Tropical connections to climatic change in the extratropical Southern Hemisphere: The role of Atlantic SST trends. J. Climate, 27, 49234936, doi:10.1175/JCLI-D-13-00615.1.

    • Search Google Scholar
    • Export Citation
  • Sinclair, M. R., , J. A. Renwick, , and J. W. Kidson, 1997: Low-frequency variability of Southern Hemisphere sea level pressure and weather system activity. Mon. Wea. Rev., 125, 25312543, doi:10.1175/1520-0493(1997)125<2531:LFVOSH>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Souders, M. B., , B. A. Colle, , and E. K. M. Chang, 2014: A description and evaluation of an automated approach for feature-based tracking of Rossby wave packets. Mon. Wea. Rev., 142, 35053527, doi:10.1175/MWR-D-13-00317.1.

    • Search Google Scholar
    • Export Citation
  • Steig, E., , Q. Ding, , D. Battisti, , and A. Jenkins, 2012: Tropical forcing of circumpolar deep water inflow and outlet glacier thinning in the Amundsen Sea Embayment, West Antarctica. Ann. Glaciol., 53, 1928, doi:10.3189/2012AoG60A110.

    • Search Google Scholar
    • Export Citation
  • Sun, D., , F. Xue, , and T. Zhou, 2013: Impacts of two types of El Niño on atmospheric circulation in the Southern Hemisphere. Adv. Atmos. Sci., 30, 17321742, doi:10.1007/s00376-013-2287-9.

    • Search Google Scholar
    • Export Citation
  • Trenberth, K. E., , and D. P. Stepaniak, 2001: Indices of El Niño evolution. J. Climate, 14, 16971701, doi:10.1175/1520-0442(2001)014<1697:LIOENO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Turner, J., , T. Phillips, , J. S. Hosking, , G. J. Marshall, , and A. Orr, 2013: The Amundsen Sea low. Int. J. Climatol., 33, 18181829, doi:10.1002/joc.3558.

    • Search Google Scholar
    • Export Citation
  • van der Walt, S., , S. C. Colbert, , and G. Varoquaux, 2011: The NumPy array: A structure for efficient numerical computation. Comput. Sci. Eng., 13, 2230, doi:10.1109/MCSE.2011.37.

    • Search Google Scholar
    • Export Citation
  • Wallace, J. M., , and D. S. Gutzler, 1981: Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Wea. Rev., 109, 784812, doi:10.1175/1520-0493(1981)109<0784:TITGHF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Wilks, D. S., 2011: Statistical Methods in the Atmospheric Sciences. 3rd ed. Academic Press, 676 pp.

  • Wilson, A. B., , D. H. Bromwich, , K. M. Hines, , and S.-H. Wang, 2014: El Niño flavors and their simulated impacts on atmospheric circulation in the high southern latitudes. J. Climate, 27, 89348955, doi:10.1175/JCLI-D-14-00296.1.

    • Search Google Scholar
    • Export Citation
  • Yu, L., , Z. Zhang, , M. Zhou, , S. Zhong, , D. Lenschow, , H. Hsu, , H. Wu, , and B. Sun, 2012: Influence of the Antarctic Oscillation, the Pacific–South American modes and the El Niño–Southern Oscillation on the Antarctic surface temperature and pressure variations. Antarct. Sci., 24, 5976, doi:10.1017/S095410201100054X.

    • Search Google Scholar
    • Export Citation
  • Yuan, X., , and D. G. Martinson, 2001: The Antarctic dipole and its predictability. Geophys. Res. Lett., 28, 36093612, doi:10.1029/2001GL012969.

    • Search Google Scholar
    • Export Citation
  • Ziȩba, A., 2010: Effective number of observations and unbiased estimators of variance for autocorrelated data—An overview. Metrol. Meas. Syst., 17, 316, doi:10.2478/v10178-010-0001-0.

    • Search Google Scholar
    • Export Citation
  • Zimin, A. V., , I. Szunyogh, , D. J. Patil, , B. R. Hunt, , and E. Ott, 2003: Extracting envelopes of Rossby wave packets. Mon. Wea. Rev., 131, 10111017, doi:10.1175/1520-0493(2003)131<1011:EEORWP>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Zimin, A. V., , I. Szunyogh, , B. R. Hunt, , and E. Ott, 2006: Extracting envelopes of nonzonally propagating Rossby wave packets. Mon. Wea. Rev., 134, 13291333, doi:10.1175/MWR3122.1.

    • Search Google Scholar
    • Export Citation
  • View in gallery

    Empirical orthogonal function (EOF) analysis of the monthly 500-hPa zonal streamfunction anomaly from the ERA-Interim reanalysis over the period 1979–2014. This is the most common method, variable, and time scale used to investigate the PSA pattern and the data are presented as the correlation of the corresponding principal component with the original field. The second and third EOF modes are degenerate according to the North et al. (1982) rule of thumb, which means the sample eignenvectors (i.e., EOF-2 and EOF-3) represent a random mixture of the true eigenvectors. Different filtering, datasets, time periods, and EOF methodologies influence the location and magnitude of the anomaly centers slightly, but the overall structure of the EOF-1 mode and the degenerate nature of EOF-2 and EOF-3 are a consistent feature. Green lines indicate the search region of interest (the “PSA sector”) and the percentage of variance explained is indicated for each EOF mode.

  • View in gallery

    Atmospheric circulation at 500 hPa for a 30-day mean centered on 18 May 2006. The top panel shows the streamfunction anomaly plotted on a regular global grid (dashed contours indicate negative values and the contour interval is ), while the bottom panel shows the corresponding meridional wind anomaly on a rotated grid where the north pole is located at 20°N, 260°E. The solid green box is bounded by 10°S–10°N and 115°–235°E on the rotated grid and corresponds to the search region of interest (or PSA sector), while the green dashed line corresponds to the equator on the rotated grid.

  • View in gallery

    Fourier analysis of the meridional average (10°S–10°N) 500-hPa rotated meridional wind anomaly for a 30-day mean centered on 18 May 2006 (purple curve). Values outside of the region of interest (115°–235°E) have been set to zero. The individual Fourier components for wavenumbers 1–8 (gray dashed, with wavenumbers 5 and 6 highlighted in green) and the reconstructed signal from an inverse Fourier transform with wavenumbers 4–7 retained (dashed orange) and its corresponding wave envelope (solid orange) are all shown. The inset shows the amplitude of each of the individual Fourier components.

  • View in gallery

    Temporal average (1979–2014) periodograms for the meridional average (10°S–10°N) 500-hPa rotated meridional wind anomaly (wind values outside of 115°–235°E were set to zero). Each curve represents a different running mean window that was applied to the daily time-scale data prior to the analysis. The vertical axis units correspond to Eq. (2).

  • View in gallery

    Life cycle characteristics of PSA-like variability. (a) The duration of all events, and (b) the phase of events lasting more than 10 data times. Events that showed substantial eastward (defined as a linear phase gradient of greater than 0.25°E per data time) or westward (less than −0.25°E per data time) propagation are colored red and blue respectively; otherwise gray shading is used. The intensity of the shading in (b) represents the amplitude of the PSA-like variability.

  • View in gallery

    Phase distribution for all data times displaying PSA-like variability. The bars show the total count for each 0.75°E interval over the period 1979–2014, while the lines represent kernel density estimates for a series of different time periods. Gray shading indicates the phase groupings taken to represent the positive (4.5°–19.5°E) and negative (37.5°–52.5°E) phase of the PSA pattern.

  • View in gallery

    Composite mean 500-hPa streamfunction anomaly for four different phase groupings: positive PSA pattern (4.5°–19.5°E), minimum 1 (22.5°–37.5°E), negative PSA pattern (37.5°–52.5°E), and minimum 2 (50.25°–6.0°E). Dashed contours indicate negative values and the contour interval is .

  • View in gallery

    Variability and trends in the negative phase of the PSA pattern. (a) The total PSA-negative data times for each individual season are shown in panel, and (c) corresponding seasonal linear trends (black represents the annual trend); (b) monthly totals for the entire study period (1979–2014). To account for the fact that not all months have an equal number of days, the counts for each month are presented as a percentage of the total number of days for that month. Years in (a) are defined from December to November (e.g., the “year” 1980 spans December 1979 to November 1980) and trends that are statistically significant at the and level are indicated with a circle and star respectively.

  • View in gallery

    As in Fig. 8, but for the positive phase of the PSA pattern.

  • View in gallery

    SAM vs Niño-3.4 for all data times (shown in black) over the period 1979–2014. Dots corresponding to PSA-positive and PSA-negative data times were recolored red and blue respectively and were arranged in the following order from front to back: blue, red, black (i.e., a blue dot might have red and/or black dots hidden underneath). For clarity, only every seventh data time was plotted. Corresponding histograms and kernel density estimates for the (top) SAM and (right) Niño-3.4 are shown and have been scaled according to density as opposed to frequency (hence the amplitudes are comparable).

  • View in gallery

    Composite mean surface air temperature anomaly, precipitation anomaly, and sea ice fraction anomaly for all data times corresponding to the (top) positive (phase grouping 4.5°–19.5°E) or (bottom) negative (37.5°–52.5°E) phase of the PSA pattern. Black contours show the composite mean 500-hPa streamfunction anomaly (dashed contours indicate negative values, and the contour interval is ), while the hatching shows regions where the difference between the composite mean and climatological mean is significant at the level.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 238 238 48
PDF Downloads 231 231 52

A New Method for Identifying the Pacific–South American Pattern and Its Influence on Regional Climate Variability

View More View Less
  • 1 School of Earth Sciences, University of Melbourne, Parkville, Victoria, Australia
© Get Permissions
Full access

Abstract

The Pacific–South American (PSA) pattern is an important mode of climate variability in the mid-to-high southern latitudes. It is widely recognized as the primary mechanism by which El Niño–Southern Oscillation (ENSO) influences the southeast Pacific and southwest Atlantic and in recent years has also been suggested as a mechanism by which longer-term tropical sea surface temperature trends can influence the Antarctic climate. This study presents a novel methodology for objectively identifying the PSA pattern. By rotating the global coordinate system such that the equator (a great circle) traces the approximate path of the pattern, the identification algorithm utilizes Fourier analysis as opposed to a traditional empirical orthogonal function approach. The climatology arising from the application of this method to ERA-Interim reanalysis data reveals that the PSA pattern has a strong influence on temperature and precipitation variability over West Antarctica and the Antarctic Peninsula and on sea ice variability in the adjacent Amundsen, Bellingshausen, and Weddell Seas. Identified seasonal trends toward the negative phase of the PSA pattern are consistent with warming observed over the Antarctic Peninsula during autumn, but are inconsistent with observed winter warming over West Antarctica. Only a weak relationship is identified between the PSA pattern and ENSO, which suggests that the pattern might be better conceptualized as a preferred regional atmospheric response to various external (and internal) forcings.

Corresponding author address: Damien Irving, School of Earth Sciences, University of Melbourne, Corner Swanston and Elgin Streets, Parkville, Victoria 3010, Australia. E-mail: irving.damien@gmail.com

Abstract

The Pacific–South American (PSA) pattern is an important mode of climate variability in the mid-to-high southern latitudes. It is widely recognized as the primary mechanism by which El Niño–Southern Oscillation (ENSO) influences the southeast Pacific and southwest Atlantic and in recent years has also been suggested as a mechanism by which longer-term tropical sea surface temperature trends can influence the Antarctic climate. This study presents a novel methodology for objectively identifying the PSA pattern. By rotating the global coordinate system such that the equator (a great circle) traces the approximate path of the pattern, the identification algorithm utilizes Fourier analysis as opposed to a traditional empirical orthogonal function approach. The climatology arising from the application of this method to ERA-Interim reanalysis data reveals that the PSA pattern has a strong influence on temperature and precipitation variability over West Antarctica and the Antarctic Peninsula and on sea ice variability in the adjacent Amundsen, Bellingshausen, and Weddell Seas. Identified seasonal trends toward the negative phase of the PSA pattern are consistent with warming observed over the Antarctic Peninsula during autumn, but are inconsistent with observed winter warming over West Antarctica. Only a weak relationship is identified between the PSA pattern and ENSO, which suggests that the pattern might be better conceptualized as a preferred regional atmospheric response to various external (and internal) forcings.

Corresponding author address: Damien Irving, School of Earth Sciences, University of Melbourne, Corner Swanston and Elgin Streets, Parkville, Victoria 3010, Australia. E-mail: irving.damien@gmail.com

1. Introduction

The Pacific–South American (PSA) pattern has long been recognized as an important mode of regional climate variability. First named by Mo and Ghil (1987), the pattern was identified in a number of studies of the large-scale Southern Hemisphere (SH) circulation during the late 1980s and early 1990s (e.g., Kidson 1988; Ghil and Mo 1991; Lau et al. 1994). A link between the pattern and Rossby wave dispersion associated with El Niño–Southern Oscillation (ENSO) was soon found (e.g., Karoly 1989), and this work was followed by a number of detailed analyses of the characteristics of the pattern and its downstream impacts (e.g., Mo and Higgins 1998; Mo 2000; Mo and Paegle 2001). In the period since these initial climatological accounts, substantial advances have been made in the methods and datasets used to identify quasi-stationary Rossby wave patterns. Given that the PSA pattern has been implicated in recent Antarctic temperature and sea ice trends, these advances could be employed to better understand the role of the pattern in high-latitude climate variability and its climatological characteristics more generally (e.g., spatial pattern, propagation, seasonal and interannual variability).

The PSA pattern is most commonly analyzed in terms of a pair of empirical orthogonal function (EOF) modes (e.g., Fig. 1). Known as PSA-1 and PSA-2, these modes are in quadrature and depict a wave train extending along an approximate great circle path from the central Pacific Ocean to the Amundsen and Weddell Seas. Some authors interpret these patterns as a single eastward propagating wave (Mo and Higgins 1998), while others argue that variability in the PSA sector is better described as a set of geographically fixed regimes (Robertson and Mechoso 2003). On a decadal time scale, PSA-1 has been related to sea surface temperature (SST) anomalies over the central and eastern Pacific, whereas on an interannual time scale it appears as a response to ENSO (Mo and Paegle 2001). The association of PSA-2 with tropical variability is less clear, with some authors relating it to the quasi-biennial component of ENSO variability (Mo 2000) and others to the Madden–Julian oscillation (Renwick and Revell 1999). While most of the features of the PSA pattern are consistent with theory and/or modeling of Rossby wave dispersion from anomalous tropical heat sources (e.g., Liu and Alexander 2007; Y. Li et al. 2015), it is recognized that the pattern can also result from internal atmospheric fluctuations caused by instabilities of the basic state (and that both mechanisms likely act in concert; e.g., Grimm and Ambrizzi 2009).

Fig. 1.
Fig. 1.

Empirical orthogonal function (EOF) analysis of the monthly 500-hPa zonal streamfunction anomaly from the ERA-Interim reanalysis over the period 1979–2014. This is the most common method, variable, and time scale used to investigate the PSA pattern and the data are presented as the correlation of the corresponding principal component with the original field. The second and third EOF modes are degenerate according to the North et al. (1982) rule of thumb, which means the sample eignenvectors (i.e., EOF-2 and EOF-3) represent a random mixture of the true eigenvectors. Different filtering, datasets, time periods, and EOF methodologies influence the location and magnitude of the anomaly centers slightly, but the overall structure of the EOF-1 mode and the degenerate nature of EOF-2 and EOF-3 are a consistent feature. Green lines indicate the search region of interest (the “PSA sector”) and the percentage of variance explained is indicated for each EOF mode.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

It has been shown that the PSA pattern plays a role in blocking events (Sinclair et al. 1997; Renwick and Revell 1999) and South American rainfall variability (Mo and Paegle 2001) and is also closely related to prominent regional features such as the Amundsen Sea low (Turner et al. 2013), Antarctic dipole (Yuan and Martinson 2001), Antarctic Circumpolar Wave (Christoph et al. 1998), and southern annular mode (SAM; e.g., Ding et al. 2012; Fogt et al. 2012). While these are all important mid-to-high-latitude impacts and relationships, in recent years the PSA pattern has been mentioned most frequently in the literature in relation to the rapid warming observed over West Antarctica and the Antarctic Peninsula (Nicolas and Bromwich 2014). In particular, it has been suggested that seasonal trends in tropical Pacific SSTs may be responsible, via circulation trends resembling the PSA pattern, for winter (and to a lesser extent spring) surface warming in West Antarctica (Ding et al. 2011), spring surface warming over the western Antarctic Peninsula (Clem and Fogt 2015), and autumn surface warming across the entire Antarctic Peninsula (Ding and Steig 2013). The pattern has also been associated with declines in sea ice in the Amundsen and Bellingshausen Seas (Schneider et al. 2012) and glacier retreat in the Amundsen Sea Embayment (Steig et al. 2012).

In identifying the PSA pattern as a possible contributor to these trends, the aforementioned studies looked through the lens of the variable(s) of interest. For instance, Ding et al. (2011) performed a maximum covariance analysis to examine the relationship between central Pacific SSTs and the broader SH circulation (the 200-hPa geopotential height). The second mode of that analysis revealed a circulation resembling the PSA pattern (and that brings warm air over West Antarctica), and atmospheric model runs forced with the associated central Pacific SSTs produced a PSA-like wave train. While this is certainly a valid research methodology, the result would be more robust if a climatology of PSA pattern activity also displayed trends consistent with warming in West Antarctica. This concept of teleconnection reversibility was recently invoked to question the relationship between Indian Ocean SSTs and heat waves in southwestern Australia (Boschat et al. 2016).

A climatology that addresses issues such as recent trends is currently lacking in the literature, so this study will present an update on our somewhat dated climatological understanding of the PSA pattern (Mo and Higgins 1998; Mo and Paegle 2001). Not only will it utilize a longer, higher-quality reanalysis dataset than previous studies, it will also develop and apply a methodology that fully exploits the capabilities of Fourier analysis, as opposed to relying on a traditional EOF-based approach. This alternative methodology was adapted from a recent climatology of SH zonal wave activity (Irving and Simmonds 2015) and seeks to avoid the issues associated with the stationary nature of spatial EOF modes, which can be problematic when trying to capture phase variations in a wave pattern of interest. This updated climatology will provide new insights into the variability, propagation, and downstream impacts of the PSA pattern, including its role in recent high-latitude trends.

2. Data

Data from the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim; Dee et al. 2011) were used in this study. In particular, the 6-hourly 500-hPa zonal and meridional wind, surface air temperature, sea ice fraction, sea surface temperature, and mean sea level pressure analysis fields were used, from which daily mean time series were calculated for the 36-yr period 1 January 1979 to 31 December 2014. For precipitation, the “total precipitation” forecast fields were used. Each forecast field represents the accumulated precipitation since initialization, so the daily rainfall total was calculated as the sum of the two 12-h post-initialization accumulation fields for each day. The horizontal resolution of the ERA-Interim data was 0.75° latitude × 0.75° longitude.

Relative to the other latest-generation reanalysis datasets, ERA-Interim is thought to best reproduce the precipitation variability (Bromwich et al. 2011; Nicolas and Bromwich 2011), vertical temperature structure (Screen and Simmonds 2012), and mean sea level pressure and 500-hPa geopotential height at station locations (Bracegirdle and Marshall 2012) around Antarctica. While these are encouraging findings, it is worth noting that the sparsity of observational data in the mid-to-high southern latitudes means that ERA-Interim data (like all reanalysis data) still need to be interpreted with caution. There are also well-known difficulties with the representation of low-frequency variability and trends in reanalysis data, due to factors such as changes in the observing system over time (Dee et al. 2014). These issues are highly relevant to the PSA pattern trends discussed in this study, but are somewhat less critical for the results pertaining to seasonal and interannual variability.

3. Computation

A number of different software packages were used in generating the key results presented in this paper. Simple editing of netCDF file attributes and routine data analysis tasks (e.g., anomalies, running means) were performed using a collection of command line utilities known as the NetCDF Operators (NCO) and Climate Data Operators (CDO) respectively, while a Python distribution called Anaconda was used for more complex analysis and visualization. With respect to specific Python libraries, xarray was used for data analysis and the reading/writing of netCDF files, which is a library that builds upon the Numerical Python (NumPy; van der Walt et al. 2011), Pandas, and Scientific Python (SciPy) libraries that come installed by default with Anaconda. Similarly, Iris, Cartopy, and Seaborn build upon Matplotlib (the default Python plotting library; Hunter 2007) and were used to generate many of the figures. Iris was also used for rotating the global coordinate system and meridional wind (via the PROJ.4 Cartographic Projections Library), and the pyqt_fit, eofs (Dawson 2016a), and windspharm (Dawson 2016b) libraries were used for kernel density estimation, for EOF analysis, and for calculating the streamfunction, respectively.

An accompanying Figshare repository has been created to document the computational methodology in more detail (Irving 2016b). It contains the specifics of the software packages discussed above (i.e., version numbers, release dates, web addresses) as well as a supplementary log file for each figure in the paper. Those log files outline the computational steps performed from initial download of the ERA-Interim data through to the final generation of the plot, and a version-controlled repository of the relevant code can be found at https://github.com/DamienIrving/climate-analysis. The rationale behind this approach to computational reproducibility is explained by Irving (2016a).

4. Methodology

While EOF analysis has been useful in highlighting the existence of the PSA pattern, it may not be the best tool for detailed climatological investigation. A particularly problematic shortcoming is that EOF analysis allows for only a crude representation of variations in wave phase (via the PSA-1 and PSA-2 modes, which are 90° out of phase), which makes it difficult to interpret characteristics such as the propagation of the pattern. This issue is further compounded by the degenerate (North et al. 1982) nature of the PSA-2 mode (e.g., Fig. 1; Mo 2000), which means the true second mode is some random combination of the apparent second and third modes.

To obtain more detailed information on variations in wave phase, existing studies of Rossby wave activity have tended to apply Fourier analysis along lines of constant latitude (e.g., Glatt and Wirth 2014). More sophisticated methods have recently been considered for identifying and tracking nonzonal waves (e.g., Zimin et al. 2006; Souders et al. 2014); however, a key insight of the method developed here is that unlike the generalized case of all possible nonzonal propagation, analysis of the PSA pattern can make use of the fact that the waveform follows an approximate great circle path (Hoskins and Karoly 1981). By rotating the global coordinate system such that the equator (itself a great circle) traces the approximate path of the PSA pattern, we were able to simply apply Fourier analysis along the “equator” in the new zonal direction. Such grid rotation is commonly used in ocean modeling to avoid coordinate singularities caused by the convergence of meridians at the poles (i.e., the grid is rotated to place the north pole over land; e.g., Bonaventura et al. 2012) but has not previously been applied in the context of tropospheric planetary wave activity. This new approach to PSA pattern identification is described below, along with the other more general data analysis techniques used in the study.

a. Identification algorithm

1) Grid rotation

To align the new equator with the approximate path of the PSA pattern, a global 0.75° latitude × 0.75° longitude grid was defined (i.e., the same resolution as the original ERA-Interim data) with the north pole located at 20°N, 260°E. The 500-hPa zonal and meridional wind data were used to calculate the meridional wind relative to the new north, and then the temporal anomaly of this new meridional wind was linearly interpolated to the rotated grid for use in the Fourier analysis (e.g., Fig. 2). It should be noted that existing zonal wave studies (e.g., Irving and Simmonds 2015) tend to skip this final step of calculating the anomaly, because in the case of zonal waves the temporal mean of the meridional wind is typically close to zero (and hence waveforms defined by the meridional wind already oscillate about zero).

Fig. 2.
Fig. 2.

Atmospheric circulation at 500 hPa for a 30-day mean centered on 18 May 2006. The top panel shows the streamfunction anomaly plotted on a regular global grid (dashed contours indicate negative values and the contour interval is ), while the bottom panel shows the corresponding meridional wind anomaly on a rotated grid where the north pole is located at 20°N, 260°E. The solid green box is bounded by 10°S–10°N and 115°–235°E on the rotated grid and corresponds to the search region of interest (or PSA sector), while the green dashed line corresponds to the equator on the rotated grid.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

On this rotated grid, the search region of interest was defined as the area bounded by 10°S–10°N, 115°–235°E (this approximate area is referred to as the PSA sector at times throughout the paper). This region was selected via visual comparison with existing definitions of the PSA pattern (e.g., Fig. 1); however, the final results were not sensitive to small changes in pole location or search region bounds.

2) Fourier analysis

To prepare the meridional wind anomaly for Fourier analysis, the meridional mean was calculated over 10°S–10°N (in order to eliminate the latitudinal dimension) and then values outside of 115°–235°E were set to zero. Zero padding is a commonly used technique in signal processing when the waveform of interest does not complete an integer number of cycles in a given domain, and is equivalent to multiplying the original signal (in this case the meridional mean meridional wind anomaly) by a square window function. This multiplication (or convolution) of two waves has consequences in frequency space, such that even a perfectly sinusoidal signal that would repeat exactly six times (for example) over the zero padded domain would show power at more than one frequency. This phenomenon is known as spectral leakage (into the side lobes of the frequency spectrum) and arises from the fact that a square window function is not square in frequency space. In analyses where excessive leakage is undesirable, a Hanning or Hamming window can be used instead. In the frequency space these windows do not display as much spread into the side lobes, but this comes at the expense of the magnitude of the main lobes. Since our selection process (see below) focuses on identifying the main lobes, a square window function was considered most appropriate.

3) Identification and characterization of PSA-like variability

Given that the PSA pattern completes approximately 1.6 to 2.0 cycles (depending on the specific EOF mode) over the 120° search area (see Fig. 1), our analysis focused on data times where a Fourier transform revealed wavenumber 5 and 6 as dominant frequencies over the zero padded 360° 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) 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 Fig. 3 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 below. To quantify the amplitude of PSA-like variability, the wave envelope construct pioneered (in the atmospheric sciences) by Zimin et al. (2003) and recently applied by Irving and Simmonds (2015) was used. The envelope is obtained by performing a Fourier transform, followed by an inverse Fourier transform for only the wavenumbers of interest (this is known as a Hilbert transform in digital signal processing). 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 Fig. 3 for an example envelope).

Fig. 3.
Fig. 3.

Fourier analysis of the meridional average (10°S–10°N) 500-hPa rotated meridional wind anomaly for a 30-day mean centered on 18 May 2006 (purple curve). Values outside of the region of interest (115°–235°E) have been set to zero. The individual Fourier components for wavenumbers 1–8 (gray dashed, with wavenumbers 5 and 6 highlighted in green) and the reconstructed signal from an inverse Fourier transform with wavenumbers 4–7 retained (dashed orange) and its corresponding wave envelope (solid orange) are all shown. The inset shows the amplitude of each of the individual Fourier components.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

4) Time-scale considerations

In applying the identification algorithm to the ERA-Interim dataset, we focus on monthly time-scale data at 500 hPa. This represents a mid-to-upper tropospheric level that is below the tropopause in all seasons and at all latitudes of interest. Given the equivalent barotropic nature of the PSA pattern (i.e., the wave amplitude increases with height but phase lines tend to be vertical) the results do not differ substantially for other levels of the troposphere. Monthly mean data were used for consistency with most previous studies and were obtained by applying a 30-day running mean to the daily (i.e., diurnally averaged) ERA-Interim data, in order to maximize the available monthly time-scale information. As noted by previous authors (e.g., Kidson 1988), potentially useful information may be lost if only 12 (i.e., calendar month) samples are taken every year. Dates were labeled as the middle day of the 30-day period and this middle day was used to determine which season a given data time belonged to (e.g., the labeled date “1979–02–16” spans the period 1 Feb 1979 to 2 Mar 1979 and belongs to DJF).

To explore the implications of this time-scale selection, the Fourier transform used in the identification process was applied to the 500-hPa rotated meridional wind anomaly data (Fig. 4). That analysis revealed wavenumber 7 as the most dominant frequency for daily time-scale data in the PSA sector, with wavenumber 6 dominating the frequency spectrum for a 10–90-day running window. Given that the PSA pattern is itself characterized by wavenumber 5–6 variability in the PSA sector, this result suggests that the PSA pattern is a dominant regional feature on weekly through to seasonal time scales, and by extension the climatological results obtained from 30-day running mean data may also be relevant at those time scales.

Fig. 4.
Fig. 4.

Temporal average (1979–2014) periodograms for the meridional average (10°S–10°N) 500-hPa rotated meridional wind anomaly (wind values outside of 115°–235°E were set to zero). Each curve represents a different running mean window that was applied to the daily time-scale data prior to the analysis. The vertical axis units correspond to Eq. (2).

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

b. Data analysis techniques

The general data analysis techniques described below are similar to those employed in the zonal wave analysis of Irving and Simmonds (2015). The following text is derived from there with minor modifications.

1) Anomalies

All anomaly data discussed in the paper represent the daily anomaly. For instance, in preparing the 30-day running mean surface air temperature anomaly data series, a 30-day running mean was first applied to the daily surface air temperature data. The mean value for each day in this 30-day running mean data series (over the entire 1979–2014 study period) was then calculated to produce a daily climatology (i.e., the multiyear daily mean). The corresponding climatological daily mean value was then subtracted at each data time to obtain the anomaly.

2) Composites

Composite mean fields are presented throughout the paper for various temporal subsets (e.g., all data times corresponding to the positive or negative phase of the PSA pattern). For the composite mean anomalies of surface temperature, precipitation, and sea ice, two-sided, one-sample t tests were applied at each grid point to examine the null hypothesis that the composite mean anomaly had been drawn from a population centered on zero. To account for autocorrelation in the data (which was substantial due to the 30-day running mean applied to the daily time-scale data), the sample size (i.e., the number of data times used in calculating the composite; denoted n) was reduced to an effective sample size according to
e1
where represents the autocorrelation for a given time lag k (Ziȩba 2010).

3) Periodogram

The characteristics of data series that have been Fourier-transformed are often summarized using a plot known as a periodogram or Fourier line spectrum (Wilks 2011). These plots are also referred to as a power or density spectrum, and most commonly display the squared amplitudes of the Fourier transform coefficients as a function of their corresponding frequencies . As an alternative to the squared amplitude, we have chosen to rescale the vertical axis and instead use the statistic commonly computed in regression analysis. The for the kth harmonic is
e2
where is the sample variance and n the length of the data series. This rescaling is particularly useful as it shows the proportion of variance in the original data series accounted for by each harmonic (Wilks 2011).

4) Climate indices

Two of the major modes of SH climate variability are the SAM and ENSO. To assess their relationship with the PSA pattern, the Antarctic Oscillation Index (Gong and Wang 1999) and the Niño-3.4 index (Trenberth and Stepaniak 2001) were calculated from 30-day running mean data (i.e., the same time scale that was used for the PSA pattern analysis). The former represents the normalized difference of zonal mean sea level pressure between 40° and 65°S, while the latter is the sea surface temperature anomaly (relative to the 1981–2000 base period) for the region in the central tropical Pacific Ocean bounded by 5°S–5°N, 190°–240°E.

5. Results

a. General PSA-like variability

Before attempting to isolate the PSA pattern using the phase information obtained from the identification algorithm, it is worth considering the characteristics of all PSA-like variability. In total, 55% (7163 of 13120) of data times were identified as displaying PSA-like variability (i.e., wavenumbers 5 and 6 were among the top three ranked frequencies), which is consistent with the fact that wavenumber 6 dominates the Fourier spectrum at the monthly time scale (Fig. 4). Grouping consecutive identifications into discrete events revealed a mean event duration of 19.7 data times, with a distribution depicted in Fig. 5a. While interpretation of these duration data is complicated by the 30-day running mean applied to the original data (e.g., an event that lasted 10 data times could be said to span anywhere between 10 and 40 days) and the occurrence of short events immediately before or after a long event (i.e., they could conceivably be considered as a single event), it appears that PSA-like variability often persists for up to a few months at a time. Building on this baseline duration data, the life cycle of events lasting longer than 10 data times was investigated in more detail. As depicted in Fig. 5b, the amplitude of these events tended to peak midevent with some longer-lasting events peaking more than once during their lifetime (perhaps suggesting that some events simply merge into the next). The mean (±1 standard deviation) linear phase trend across all events lasting longer than 10 data times was 0.12° ± 0.38°E per data time, which indicates that while there was a tendency for events to propagate to the east, a substantial proportion moved very little (or even to the west) during their lifetime.

Fig. 5.
Fig. 5.

Life cycle characteristics of PSA-like variability. (a) The duration of all events, and (b) the phase of events lasting more than 10 data times. Events that showed substantial eastward (defined as a linear phase gradient of greater than 0.25°E per data time) or westward (less than −0.25°E per data time) propagation are colored red and blue respectively; otherwise gray shading is used. The intensity of the shading in (b) represents the amplitude of the PSA-like variability.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

Important insights were also obtained by considering the phase distribution across all individual PSA-like data times (Fig. 6). On an annual basis the distribution is clearly bimodal, with the two maxima of the kernel density estimate located at 12.75° and 45.0°E. Since the phase was defined as the location of the first local maximum of the wavenumber 6 component of the Fourier transform, this approximate 30° phase separation indicates a pair of spatial patterns that are exactly out of phase (Fig. 7). Taken together these patterns clearly represent the single most dominant mode of variability in the PSA sector, and closely resemble the PSA-1 mode identified by previous authors. On the basis of this finding, it appears that filtering the PSA-like data times according to the location of the two local maxima represents a simple and valid technique for isolating the PSA pattern from the larger population of PSA-like variability.

Fig. 6.
Fig. 6.

Phase distribution for all data times displaying PSA-like variability. The bars show the total count for each 0.75°E interval over the period 1979–2014, while the lines represent kernel density estimates for a series of different time periods. Gray shading indicates the phase groupings taken to represent the positive (4.5°–19.5°E) and negative (37.5°–52.5°E) phase of the PSA pattern.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

Fig. 7.
Fig. 7.

Composite mean 500-hPa streamfunction anomaly for four different phase groupings: positive PSA pattern (4.5°–19.5°E), minimum 1 (22.5°–37.5°E), negative PSA pattern (37.5°–52.5°E), and minimum 2 (50.25°–6.0°E). Dashed contours indicate negative values and the contour interval is .

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

The spatial patterns corresponding to the local minima of the phase distribution are also shown in Fig. 7, as a way to summarize the characteristics of the remaining PSA-like variability. The three anomaly centers associated with these composite mean circulation patterns have different amplitudes (the middle anomaly has a larger amplitude than the others), which indicates that it was often not a coordinated wave pattern that the identification algorithm was picking up (i.e., not the coordinated PSA-2 waveform discussed by previous authors, despite the similarity in wave phase). Looking at the individual data times corresponding to those minima (not shown), they appear to be a mixture of the hemispheric zonal wave 3 pattern (Raphael 2004; Irving and Simmonds 2015), a more meridionally oriented wave train extending from the tropical Pacific to the Amundsen Sea (e.g., Clem and Fogt 2015; Clem and Renwick 2015), and isolated Amundsen Sea low variability.

b. The PSA pattern

In defining the PSA pattern according to the peaks of the PSA-like phase distribution, it was necessary to account for seasonal variations in the location of those peaks (Fig. 6). A spread of 15° was considered sufficient to capture these variations and hence the 15° interval about each local maximum containing the highest mean values (taken from the annual kernel density estimate) was determined. This approach was used to account for the fact that the phase histograms were not symmetrical about the local maxima and it yielded two intervals corresponding to the positive (4.5°–19.5°E) and negative (37.5°–52.5°E) phase of the PSA pattern. Both intervals represented approximately 15% of all data times (14.8% for the positive phase vs 15.8% for the negative), which suggests that the two phases have a similar frequency of occurrence. With this definition in place, it was possible to investigate variability and trends in the PSA pattern as well as its influence on surface temperature, precipitation, and sea ice.

1) Trends and variability

During autumn and winter in particular, the middle years of the study period (1991–2002) were characterized by a predominance of positive PSA pattern activity, whereas negative phase activity was more common in recent years (Fig. 6). This variability is reflected in the linear trends observed over 1979–2014, with negative phase activity showing a statistically significant increasing trend (at the level) on an annual basis and smaller nonsignificant increasing trends for summer, autumn, and winter (Fig. 8). Positive phase activity showed a nonsignificant decreasing trend on an annual basis and also during autumn and winter, with an increasing trend observed for summer (Fig. 9). Consistent with previous studies, both phases of the PSA pattern were most active during winter and spring (Figs. 8 and 9).

Fig. 8.
Fig. 8.

Variability and trends in the negative phase of the PSA pattern. (a) The total PSA-negative data times for each individual season are shown in panel, and (c) corresponding seasonal linear trends (black represents the annual trend); (b) monthly totals for the entire study period (1979–2014). To account for the fact that not all months have an equal number of days, the counts for each month are presented as a percentage of the total number of days for that month. Years in (a) are defined from December to November (e.g., the “year” 1980 spans December 1979 to November 1980) and trends that are statistically significant at the and level are indicated with a circle and star respectively.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

Fig. 9.
Fig. 9.

As in Fig. 8, but for the positive phase of the PSA pattern.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

In attempting to explain annual and decadal variability in PSA pattern activity, previous authors have suggested that coupling between the SAM and ENSO is important (e.g., Fogt and Bromwich 2006). While some degree of coupling is evident in Fig. 10 (e.g., the positive phase of the PSA pattern was most common when positive/warm ENSO events and negative SAM events coincided), it is clear that the SAM has a much stronger association with PSA pattern activity than ENSO. Recent positive trends in the SAM during summer, autumn, and to a lesser extent winter (the latter being smaller and not statistically significant; e.g., Simmonds 2015) are also broadly consistent with the negative trends observed in the PSA pattern during those seasons.

Fig. 10.
Fig. 10.

SAM vs Niño-3.4 for all data times (shown in black) over the period 1979–2014. Dots corresponding to PSA-positive and PSA-negative data times were recolored red and blue respectively and were arranged in the following order from front to back: blue, red, black (i.e., a blue dot might have red and/or black dots hidden underneath). For clarity, only every seventh data time was plotted. Corresponding histograms and kernel density estimates for the (top) SAM and (right) Niño-3.4 are shown and have been scaled according to density as opposed to frequency (hence the amplitudes are comparable).

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

2) Influence on surface variables

To assess the influence of the PSA pattern on regional climate variability, the composite mean surface air temperature anomaly, precipitation anomaly, and sea ice concentration anomaly was calculated for both the positive and negative phase (Fig. 11). On the western flank of the central composite mean streamfunction anomaly associated with positive phase activity, anomalously warm conditions were evident over the Ross Sea, Amundsen Sea, and interior of West Antarctica, particularly during autumn and winter (seasonal composites not shown). The northerly flow responsible for those warm conditions also induced large precipitation increases along the West Antarctic coastline and reduced sea ice in the Amundsen Sea. On the eastern flank, anomalously cool conditions were evident over the Antarctic Peninsula, Patagonia, and the Weddell Sea during all seasons (winter and spring especially; not shown), with the latter also experiencing large increases in sea ice. Anomalously dry conditions were also seen over the Antarctic Peninsula in association with the weaker westerly flow.

Fig. 11.
Fig. 11.

Composite mean surface air temperature anomaly, precipitation anomaly, and sea ice fraction anomaly for all data times corresponding to the (top) positive (phase grouping 4.5°–19.5°E) or (bottom) negative (37.5°–52.5°E) phase of the PSA pattern. Black contours show the composite mean 500-hPa streamfunction anomaly (dashed contours indicate negative values, and the contour interval is ), while the hatching shows regions where the difference between the composite mean and climatological mean is significant at the level.

Citation: Journal of Climate 29, 17; 10.1175/JCLI-D-15-0843.1

The anomalies associated with the negative phase of the PSA pattern were essentially the reverse of the positive phase (Fig. 11). It is also noteworthy that while the hemispheric composite mean streamfunction anomaly associated with the PSA pattern gives the impression of a hemispheric zonal wavenumber 3 pattern, the phase of that pattern and the unremarkable anomalies either side of the Indian Ocean anomaly are inconsistent with the characteristics of the dominant SH zonal wavenumber 3 mode (e.g., Raphael 2004; Irving and Simmonds 2015).

6. Discussion

A novel methodology has been presented for objectively identifying the PSA pattern. By rotating the global coordinate system such that the equator (a great circle path) traces the approximate path of the PSA pattern, the method was able to utilize Fourier analysis to quantify the phase and amplitude of wavelike variability in the PSA sector. In reconciling the results of this Fourier analysis with existing EOF-based definitions of the PSA pattern, a strong resemblance was found between the existing PSA-1 mode and the spatial pattern corresponding to the bimodal phase peaks of wavenumber 5–6 dominant variability in the PSA sector. The lack of a higher-order, multimodal phase distribution questions the physical reality of the existing PSA-2 mode, and may explain the difficulty that researchers have had in identifying a tropical driver for that mode.

These bimodal phase peaks were used as a means to define the positive and negative phase of the PSA pattern. The climatology arising from this definition revealed that the PSA pattern is most active during winter and spring, often persisting for months at a time. It propagates to the east on average, but a substantial number of events remain relatively stationary or even propagate to the west. The pattern was also shown to have a strong influence on regional temperature, precipitation, and sea ice variability. With respect to the former, our results confirm existing relationships established between pattern and station temperatures over the Antarctic Peninsula (e.g., Schneider et al. 2012; Yu et al. 2012), extending the regional picture to highlight equally strong temperature anomalies (of opposite sign) over West Antarctica. Large precipitation anomalies were also identified along the coast of West Antarctica and the Antarctic Peninsula, as well as over South America. These South American anomalies show a more complex spatial pattern than previous analyses (perhaps due to the higher-resolution data), but are otherwise broadly consistent with the results of Mo and Paegle (2001), who found the positive phase of the PSA pattern to be associated with anomalously wet conditions over southern South America and anomalously dry conditions farther north. Previous studies also indicate that the PSA pattern plays an important role in sea ice variability in the Amundsen and Bellingshausen Seas (Raphael and Hobbs 2014). Our results suggest that this role is not uniform across that region, with composites of the positive phase of the PSA pattern simultaneously displaying positive sea ice anomalies in the Bellingshausen Sea and negative in the Amundsen Sea.

With respect to trends in the PSA pattern over the period 1979–2014, a trend toward the negative phase was identified on an annual basis and also during summer, autumn and winter. This autumn trend (and the high-latitude temperature and sea ice anomalies associated with the negative phase of the PSA pattern) is consistent with the work of Ding and Steig (2013), who found that autumn warming over the Antarctic Peninsula and associated sea ice declines over the Bellingshausen Sea are associated with an atmospheric circulation resembling the negative phase of the PSA pattern. While this explanation makes sense on the eastern flank of the central circulation anomaly associated with that pattern, the negative phase of the PSA pattern is also associated with strong cooling over West Antarctica. Autumn temperature declines have not been observed in that region, and thus our results suggest that the PSA-related cooling must have been offset by other factors.

In contrast to the autumn warming over the Antarctic Peninsula, winter warming over West Antarctica has been associated with an atmospheric circulation resembling the positive phase of the PSA pattern (Ding et al. 2011). Our climatology revealed a (albeit nonsignificant) trend toward the negative phase of the PSA pattern during winter, which raises this question: How is it that winter temperature trends over West Antarctica are associated with an atmospheric circulation resembling the positive phase of the PSA pattern, but a climatology of PSA pattern activity reveals a trend that directly opposes that finding? One possible answer to this question comes from X. Li et al. (2015). They analyzed Rossby wave trains associated with observed SST trends in the tropical Atlantic, tropical Indian, west Pacific, and east Pacific regions and found that all four have a center of action over the Amundsen Sea. While none of these individual wave trains resembled the PSA pattern, a linear combination of the four of them did (with the tropical Atlantic and west Pacific identified as most influential). In other words, the integrated influence of tropical SST trends on the atmospheric circulation resembles the positive phase of the PSA pattern, but the waves underpinning that teleconnection do not. This result is consistent with an earlier study that identified the tropical Atlantic as a driver of recent winter trends in West Antarctica (Li et al. 2014). Another possible answer comes from Fogt and Wovrosh (2015), who suggest that radiative forcing has played a role in Amundsen Sea low trends that are consistent with winter warming in West Antarctica. The absence of any springtime trend in the PSA pattern suggests that it has also not played a role in high-latitude warming during that season. Similar to winter, the Atlantic has been linked to warming in West Antarctica during spring (Simpkins et al. 2014), while others point to a more meridionally oriented wave train associated with the Pacific decadal oscillation (PDO; Clem and Fogt 2015; Clem and Renwick 2015).

This idea that radiatively forced Amundsen Sea low variability and/or wave trains associated with the Atlantic or PDO might be responsible for a teleconnection resembling the PSA pattern (i.e., as opposed to changes in actual PSA pattern activity) goes to the heart of the reversibility argument made at the beginning of this paper. For a proposed teleconnection to be robust, it must be evident when looking through the lens of both the variable and mechanism of interest. However, even if these alternative explanations do reconcile the discrepancy between our climatology and winter warming over West Antarctica, the associated circulation anomaly would bring cooler conditions and wind-driven increases in sea ice along the western Antarctic Peninsula, contrary to the observed warming and sea ice declines (Clem and Fogt 2015). One possible explanation is that the negative autumn sea ice anomalies persist into winter (Ding and Steig 2013); however, it is clear that there is still work to be done to fully understand recent winter temperature and sea ice changes in the region.

One topic not addressed here is variability in the east/west location of the PSA pattern. In response to the emergence of central Pacific ENSO events in recent years (e.g., Ashok et al. 2007), some authors have suggested that the PSA pattern moves east or west depending on the precise location of the associated tropical SST anomalies (e.g., Sun et al. 2013; Wilson et al. 2014; Ciasto et al. 2015). Others suggest that the pattern is relatively stationary (e.g., Liu and Alexander 2007; Ding et al. 2012), but either way the broad region (10°N to 10°S in the rotated coordinate system) used by our identification algorithm renders it insensitive to subtle east–west movements. Given that the PSA pattern did not show a strong association with the Niño-3.4 index (an index that is sensitive to both central and eastern Pacific ENSO events), it would be fair to say that even if the location of tropical SSTs does cause the pattern to move slightly, this would represent only a small fraction of all PSA pattern activity.

This weak association with ENSO challenges our fundamental understanding of the PSA pattern. The most commonly held view to date is that the pattern is primarily a response to ENSO forcing (e.g., Mo and Paegle 2001), whereby anomalous ENSO-related SST anomalies modify tropical convection, leading to atmospheric vorticity gradients conducive to Rossby wave generation (Sardeshmukh and Hoskins 1988). A more comprehensive analysis of the relationship between the pattern and tropical convection would be required to confirm this (e.g., lagged correlations with SSTs and other indicators of tropical convection like the outgoing longwave radiation), but our results suggest that the PSA pattern might actually be better conceptualized as a preferred regional atmospheric response to various internal and external forcings (i.e., with ENSO being just one of many players). Tropical convection is almost certainly an important factor, but its role is likely more complicated than can be captured by a broad-scale ENSO index like the Niño-3.4. This complex relationship between tropical SSTs, convection and the PSA pattern has been the focus of a small number of studies (e.g., Harangozo 2004; Lachlan-Cope and Connolley 2006) and is a topic that warrants much more attention. The relationship between the PSA pattern and ENSO has also traditionally been thought to be moderated by the state of the “atmospheric bridge” (Liu and Alexander 2007). In particular, the pattern is thought to be most active when ENSO and the SAM are in phase (Fogt and Bromwich 2006). Rather than casting the SAM in a facilitating or bridging role, the strong association identified here is more consistent with the idea that the PSA pattern is an integral part of the zonally asymmetric structure of the SAM (e.g., Ding et al. 2012; Fogt et al. 2012).

In addition to a more detailed analysis of the relationship between the PSA pattern, tropical convection, and the SAM, our new methodology could also be adapted for use in future studies of other quasi-stationary waveforms. The most obvious candidate is the Pacific–North American (PNA) pattern (Wallace and Gutzler 1981), which plays an important role in winter climate variability over the North Pacific and North America (e.g., Notaro et al. 2006). Like its SH counterpart, the PNA pattern follows an approximate great circle path, has traditionally been analyzed via EOF analysis, and has been implicated in recent mid-to-high latitude trends (e.g., Ding et al. 2014; Liu et al. 2015). Other nonzonal waveforms that do not follow an approximate great circle path would be more challenging; however, methods have been developed for applying Fourier analysis to synoptic-scale, nonzonal waveforms (Zimin et al. 2006; Souders et al. 2014) and may represent a starting point for future research.

Acknowledgments

Parts of this work were made possible by a grant (DP160101997) from the Australian Research Council.

REFERENCES

  • Ashok, K., , S. K. Behera, , S. A. Rao, , H. Weng, , and T. Yamagata, 2007: El Niño Modoki and its possible teleconnection. J. Geophys. Res., 112, C11007, doi:10.1029/2006JC003798.

    • Search Google Scholar
    • Export Citation
  • Bonaventura, L., , R. Redler, , and R. Budich, 2012: Algorithms, Code Infrastructure and Optimisation. Vol. 2, Earth System Modelling, Springer, 96 pp.

  • Boschat, G., , I. Simmonds, , A. Purich, , T. Cowan, , and A. B. Pezza, 2016: On the use of composite analyses to form physical hypotheses: An example from heat wave–SST associations. Sci. Rep., 6, 29599, doi:10.1038/srep29599.

    • Search Google Scholar
    • Export Citation
  • Bracegirdle, T. J., , and G. J. Marshall, 2012: The reliability of Antarctic tropospheric pressure and temperature in the latest global reanalyses. J. Climate, 25, 71387146, doi:10.1175/JCLI-D-11-00685.1.

    • Search Google Scholar
    • Export Citation
  • Bromwich, D. H., , J. P. Nicolas, , and A. J. Monaghan, 2011: An assessment of precipitation changes over Antarctica and the Southern Ocean since 1989 in contemporary global reanalyses. J. Climate, 24, 41894209, doi:10.1175/2011JCLI4074.1.

    • Search Google Scholar
    • Export Citation
  • Christoph, M., , T. P. Barnett, , and E. Roeckner, 1998: The Antarctic Circumpolar Wave in a coupled ocean–atmosphere GCM. J. Climate, 11, 16591672, doi:10.1175/1520-0442(1998)011<1659:TACWIA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Ciasto, L. M., , G. R. Simpkins, , and M. H. England, 2015: Teleconnections between tropical Pacific SST anomalies and extratropical Southern Hemisphere climate. J. Climate, 28, 5665, doi:10.1175/JCLI-D-14-00438.1.

    • Search Google Scholar
    • Export Citation
  • Clem, K. R., , and R. L. Fogt, 2015: South Pacific circulation changes and their connection to the tropics and regional Antarctic warming in austral spring, 1979–2012. J. Geophys. Res. Atmos., 120, 27732792, doi:10.1002/2014JD022940.

    • Search Google Scholar
    • Export Citation
  • Clem, K. R., , and J. A. Renwick, 2015: Austral spring Southern Hemisphere circulation and temperature changes and links to the SPCZ. J. Climate, 28, 73717384, doi:10.1175/JCLI-D-15-0125.1.

    • Search Google Scholar
    • Export Citation
  • Dawson, A., 2016a: eofs: A library for EOF analysis of meteorological, oceanographic and climate data. J. Open Res. Software, 4, e14, doi:10.5334/jors.122.

    • Search Google Scholar
    • Export Citation
  • Dawson, A., 2016b: Windspharm: A High-Level Library for Global Wind Field Computations Using Spherical Harmonics. J. Open Res. Software, 4, e31, doi:10.5334/jors.129.

    • Search Google Scholar
    • Export Citation
  • Dee, D. P., and et al. , 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553597, doi:10.1002/qj.828.

    • Search Google Scholar
    • Export Citation
  • Dee, D. P., , M. Balmaseda, , G. Balsamo, , R. Engelen, , A. J. Simmons, , and J.-N. Thépaut, 2014: Toward a consistent reanalysis of the climate system. Bull. Amer. Meteor. Soc., 95, 12351248, doi:10.1175/BAMS-D-13-00043.1.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , and E. J. Steig, 2013: Temperature change on the Antarctic Peninsula linked to the tropical Pacific. J. Climate, 26, 75707585, doi:10.1175/JCLI-D-12-00729.1.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , E. J. Steig, , D. S. Battisti, , and M. Küttel, 2011: Winter warming in West Antarctica caused by central tropical Pacific warming. Nat. Geosci., 4, 398403, doi:10.1038/ngeo1129.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , E. J. Steig, , D. S. Battisti, , and J. M. Wallace, 2012: Influence of the tropics on the southern annular mode. J. Climate, 25, 63306348, doi:10.1175/JCLI-D-11-00523.1.

    • Search Google Scholar
    • Export Citation
  • Ding, Q., , J. M. Wallace, , D. S. Battisti, , E. J. Steig, , A. J. E. Gallant, , H.-J. Kim, , and L. Geng, 2014: Tropical forcing of the recent rapid Arctic warming in northeastern Canada and Greenland. Nature, 509, 209212, doi:10.1038/nature13260.

    • Search Google Scholar
    • Export Citation
  • Fogt, R. L., , and D. H. Bromwich, 2006: Decadal variability of the ENSO teleconnection to the high-latitude South Pacific governed by coupling with the southern annular mode. J. Climate, 19, 979997, doi:10.1175/JCLI3671.1.

    • Search Google Scholar
    • Export Citation
  • Fogt, R. L., , and A. J. Wovrosh, 2015: The relative influence of tropical sea surface temperatures and radiative forcing on the Amundsen Sea low. J. Climate, 28, 85408555, doi:10.1175/JCLI-D-15-0091.1.

    • Search Google Scholar
    • Export Citation
  • Fogt, R. L., , J. M. Jones, , and J. Renwick, 2012: Seasonal zonal asymmetries in the southern annular mode and their impact on regional temperature anomalies. J. Climate, 25, 62536270, doi:10.1175/JCLI-D-11-00474.1.

    • Search Google Scholar
    • Export Citation
  • Ghil, M., , and K. Mo, 1991: Intraseasonal oscillations in the global atmosphere. Part II: Southern Hemisphere. J. Atmos. Sci., 48, 780790, doi:10.1175/1520-0469(1991)048<0780:IOITGA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Glatt, I., , and V. Wirth, 2014: Identifying Rossby wave trains and quantifying their properties. Quart. J. Roy. Meteor. Soc., 140, 384396, doi:10.1002/qj.2139.

    • Search Google Scholar
    • Export Citation
  • Gong, D., , and S. Wang, 1999: Definition of Antarctic Oscillation index. Geophys. Res. Lett., 26, 459462, doi:10.1029/1999GL900003.

  • Grimm, A. M., , and T. Ambrizzi, 2009: Teleconnections into South America from the tropics and extratropics on interannual and intraseasonal timescales. Past Climate Variability in South America and Surrounding Regions, F. Vimeux, F. Sylvestre, and M. Khodri, Eds., Springer, 159–191, doi:10.1007/978-90-481-2672-9_7.

  • Harangozo, S. A., 2004: The relationship of Pacific deep tropical convection to the winter and springtime extratropical atmospheric circulation of the South Pacific in El Niño events. Geophys. Res. Lett., 31, L05206, doi:10.1029/2003GL018667.

    • Search Google Scholar
    • Export Citation
  • Hoskins, B. J., , and D. J. Karoly, 1981: The steady linear response of a spherical atmosphere to thermal and orographic forcing. J. Atmos. Sci., 38, 11791196, doi:10.1175/1520-0469(1981)038<1179:TSLROA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Hunter, J. D., 2007: Matplotlib: A 2D graphics environment. Comput. Sci. Eng., 9, 9095, doi:10.1109/MCSE.2007.55.

  • Irving, D. B., 2016a: A minimum standard for publishing computational results in the weather and climate sciences. Bull. Amer. Meteor. Soc., doi:10.1175/BAMS-D-15-00010.1, in press.

    • Search Google Scholar
    • Export Citation
  • Irving, D. B., 2016b: A new method for identifying the Pacific–South American pattern and its influence on regional climate variability: Supplementary metadata. Figshare, digital media, doi:10.6084/m9.figshare.1612460.

  • Irving, D. B., , and I. Simmonds, 2015: A novel approach to diagnosing Southern Hemisphere planetary wave activity and its influence on regional climate variability. J. Climate, 28, 90419057, doi:10.1175/JCLI-D-15-0287.1.

    • Search Google Scholar
    • Export Citation
  • Karoly, D. J., 1989: Southern Hemisphere circulation features associated with El Niño–Southern Oscillation events. J. Climate, 2, 12391252, doi:10.1175/1520-0442(1989)002<1239:SHCFAW>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Kidson, J. W., 1988: Interannual variations in the Southern Hemisphere circulation. J. Climate, 1, 11771198, doi:10.1175/1520-0442(1988)001<1177:IVITSH>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Lachlan-Cope, T., , and W. Connolley, 2006: Teleconnections between the tropical Pacific and the Amundsen-Bellinghausens Sea: Role of the El Niño/Southern Oscillation. J. Geophys. Res., 111, D23101, doi:10.1029/2005JD006386.

    • Search Google Scholar
    • Export Citation
  • Lau, K.-M., , P.-J. Sheu, , and I.-S. Kang, 1994: Multiscale low-frequency circulation modes in the global atmosphere. J. Atmos. Sci., 51, 11691193, doi:10.1175/1520-0469(1994)051<1169:MLFCMI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Li, X., , D. M. Holland, , E. P. Gerber, , and C. Yoo, 2014: Impacts of the north and tropical Atlantic Ocean on the Antarctic Peninsula and sea ice. Nature, 505, 538542, doi:10.1038/nature12945.

    • Search Google Scholar
    • Export Citation
  • Li, X., , D. M. Holland, , E. P. Gerber, , and C. Yoo, 2015: Rossby waves mediate impacts of tropical oceans on West Antarctic atmospheric circulation in austral winter. J. Climate, 28, 81518164, doi:10.1175/JCLI-D-15-0113.1.

    • Search Google Scholar
    • Export Citation
  • Li, Y., , J. Li, , F. F. Jin, , and S. Zhao, 2015: Interhemispheric propagation of stationary Rossby waves in a horizontally nonuniform background flow. J. Atmos. Sci., 72, 32333256, doi:10.1175/JAS-D-14-0239.1.

    • Search Google Scholar
    • Export Citation
  • Liu, Z., , and M. Alexander, 2007: Atmospheric bridge, oceanic tunnel, and global climatic teleconnections. Rev. Geophys., 45, RG2005, doi:10.1029/2005RG000172.

    • Search Google Scholar
    • Export Citation
  • Liu, Z., , Z. Jian, , K. Yoshimura, , N. H. Buenning, , C. J. Poulsen, , and G. J. Bowen, 2015: Recent contrasting winter temperature changes over North America linked to enhanced positive Pacific–North American pattern. Geophys. Res. Lett., 42, 77507757, doi:10.1002/2015GL065656.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., 2000: Relationships between low-frequency variability in the Southern Hemisphere and sea surface temperature anomalies. J. Climate, 13, 35993610, doi:10.1175/1520-0442(2000)013<3599:RBLFVI>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., , and M. Ghil, 1987: Statistics and dynamics of persistent anomalies. J. Atmos. Sci., 44, 877902, doi:10.1175/1520-0469(1987)044<0877:SADOPA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., , and R. W. Higgins, 1998: The Pacific–South American modes and tropical convection during the Southern Hemisphere winter. Mon. Wea. Rev., 126, 15811596, doi:10.1175/1520-0493(1998)126<1581:TPSAMA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Mo, K. C., , and J. N. Paegle, 2001: The Pacific–South American modes and their downstream effects. Int. J. Climatol., 21, 12111229, doi:10.1002/joc.685.

    • Search Google Scholar
    • Export Citation
  • Nicolas, J. P., , and D. H. Bromwich, 2011: Precipitation changes in high southern latitudes from global reanalyses: A cautionary tale. Surv. Geophys., 32, 475494, doi:10.1007/s10712-011-9114-6.

    • Search Google Scholar
    • Export Citation
  • Nicolas, J. P., , and D. H. Bromwich, 2014: New reconstruction of Antarctic near-surface temperatures: Multidecadal trends and reliability of global reanalyses. J. Climate, 27, 80708093, doi:10.1175/JCLI-D-13-00733.1.

    • Search Google Scholar
    • Export Citation
  • North, G. R., , T. L. Bell, , and R. F. Cahalan, 1982: Sampling errors in the estimation of empirical orthogonal functions. Mon. Wea. Rev., 110, 699706, doi:10.1175/1520-0493(1982)110<0699:SEITEO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Notaro, M., , W.-C. Wang, , and W. Gong, 2006: Model and observational analysis of the northeast U.S. regional climate and its relationship to the PNA and NAO patterns during early winter. Mon. Wea. Rev., 134, 34793505, doi:10.1175/MWR3234.1.

    • Search Google Scholar
    • Export Citation
  • Raphael, M. N., 2004: A zonal wave 3 index for the Southern Hemisphere. Geophys. Res. Lett., 31, L23212, doi:10.1029/2004GL020365.

  • Raphael, M. N., , and W. R. Hobbs, 2014: The influence of the large-scale atmospheric circulation on Antarctic sea ice during ice advance and retreat seasons. Geophys. Res. Lett., 41, 50375045, doi:10.1002/2014GL060365.

    • Search Google Scholar
    • Export Citation
  • Renwick, J. A., , and M. J. Revell, 1999: Blocking over the South Pacific and Rossby wave propagation. Mon. Wea. Rev., 127, 22332247, doi:10.1175/1520-0493(1999)127<2233:BOTSPA>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Robertson, A. W., , and C. R. Mechoso, 2003: Circulation regimes and low-frequency oscillations in the South Pacific sector. Mon. Wea. Rev., 131, 15661576, doi:10.1175//2548.1.

    • Search Google Scholar
    • Export Citation
  • Sardeshmukh, P. D., , and B. J. Hoskins, 1988: The generation of global rotational flow by steady idealized tropical divergence. J. Atmos. Sci., 45, 12281251, doi:10.1175/1520-0469(1988)045<1228:TGOGRF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Schneider, D. P., , C. Deser, , and Y. Okumura, 2012: An assessment and interpretation of the observed warming of West Antarctica in the austral spring. Climate Dyn., 38, 323347, doi:10.1007/s00382-010-0985-x.

    • Search Google Scholar
    • Export Citation
  • Screen, J. A., , and I. Simmonds, 2012: Half-century air temperature change above Antarctica: Observed trends and spatial reconstructions. J. Geophys. Res., 117, D16108, doi:10.1029/2012JD017885.

    • Search Google Scholar
    • Export Citation
  • Simmonds, I., 2015: Comparing and contrasting the behaviour of Arctic and Antarctic sea ice over the 35 year period 1979–2013. Ann. Glaciol., 56, 1828, doi:10.3189/2015AoG69A909.

    • Search Google Scholar
    • Export Citation
  • Simpkins, G. R., , S. McGregor, , A. S. Taschetto, , L. M. Ciasto, , and M. H. England, 2014: Tropical connections to climatic change in the extratropical Southern Hemisphere: The role of Atlantic SST trends. J. Climate, 27, 49234936, doi:10.1175/JCLI-D-13-00615.1.

    • Search Google Scholar
    • Export Citation
  • Sinclair, M. R., , J. A. Renwick, , and J. W. Kidson, 1997: Low-frequency variability of Southern Hemisphere sea level pressure and weather system activity. Mon. Wea. Rev., 125, 25312543, doi:10.1175/1520-0493(1997)125<2531:LFVOSH>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Souders, M. B., , B. A. Colle, , and E. K. M. Chang, 2014: A description and evaluation of an automated approach for feature-based tracking of Rossby wave packets. Mon. Wea. Rev., 142, 35053527, doi:10.1175/MWR-D-13-00317.1.

    • Search Google Scholar
    • Export Citation
  • Steig, E., , Q. Ding, , D. Battisti, , and A. Jenkins, 2012: Tropical forcing of circumpolar deep water inflow and outlet glacier thinning in the Amundsen Sea Embayment, West Antarctica. Ann. Glaciol., 53, 1928, doi:10.3189/2012AoG60A110.

    • Search Google Scholar
    • Export Citation
  • Sun, D., , F. Xue, , and T. Zhou, 2013: Impacts of two types of El Niño on atmospheric circulation in the Southern Hemisphere. Adv. Atmos. Sci., 30, 17321742, doi:10.1007/s00376-013-2287-9.

    • Search Google Scholar
    • Export Citation
  • Trenberth, K. E., , and D. P. Stepaniak, 2001: Indices of El Niño evolution. J. Climate, 14, 16971701, doi:10.1175/1520-0442(2001)014<1697:LIOENO>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Turner, J., , T. Phillips, , J. S. Hosking, , G. J. Marshall, , and A. Orr, 2013: The Amundsen Sea low. Int. J. Climatol., 33, 18181829, doi:10.1002/joc.3558.

    • Search Google Scholar
    • Export Citation
  • van der Walt, S., , S. C. Colbert, , and G. Varoquaux, 2011: The NumPy array: A structure for efficient numerical computation. Comput. Sci. Eng., 13, 2230, doi:10.1109/MCSE.2011.37.

    • Search Google Scholar
    • Export Citation
  • Wallace, J. M., , and D. S. Gutzler, 1981: Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Wea. Rev., 109, 784812, doi:10.1175/1520-0493(1981)109<0784:TITGHF>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Wilks, D. S., 2011: Statistical Methods in the Atmospheric Sciences. 3rd ed. Academic Press, 676 pp.

  • Wilson, A. B., , D. H. Bromwich, , K. M. Hines, , and S.-H. Wang, 2014: El Niño flavors and their simulated impacts on atmospheric circulation in the high southern latitudes. J. Climate, 27, 89348955, doi:10.1175/JCLI-D-14-00296.1.

    • Search Google Scholar
    • Export Citation
  • Yu, L., , Z. Zhang, , M. Zhou, , S. Zhong, , D. Lenschow, , H. Hsu, , H. Wu, , and B. Sun, 2012: Influence of the Antarctic Oscillation, the Pacific–South American modes and the El Niño–Southern Oscillation on the Antarctic surface temperature and pressure variations. Antarct. Sci., 24, 5976, doi:10.1017/S095410201100054X.

    • Search Google Scholar
    • Export Citation
  • Yuan, X., , and D. G. Martinson, 2001: The Antarctic dipole and its predictability. Geophys. Res. Lett., 28, 36093612, doi:10.1029/2001GL012969.

    • Search Google Scholar
    • Export Citation
  • Ziȩba, A., 2010: Effective number of observations and unbiased estimators of variance for autocorrelated data—An overview. Metrol. Meas. Syst., 17, 316, doi:10.2478/v10178-010-0001-0.

    • Search Google Scholar
    • Export Citation
  • Zimin, A. V., , I. Szunyogh, , D. J. Patil, , B. R. Hunt, , and E. Ott, 2003: Extracting envelopes of Rossby wave packets. Mon. Wea. Rev., 131, 10111017, doi:10.1175/1520-0493(2003)131<1011:EEORWP>2.0.CO;2.

    • Search Google Scholar
    • Export Citation
  • Zimin, A. V., , I. Szunyogh, , B. R. Hunt, , and E. Ott, 2006: Extracting envelopes of nonzonally propagating Rossby wave packets. Mon. Wea. Rev., 134, 13291333, doi:10.1175/MWR3122.1.

    • Search Google Scholar
    • Export Citation
Save