The Early 1990s Change in ENSO–PSA–SAM Relationships and Its Impact on Southern Hemisphere Climate

This study uncovers an early 1990s change in the relationships between El Ni ñ o–Southern Oscillation (ENSO) and two leading modes of the Southern Hemisphere (SH) atmospheric variability: the southern annular mode (SAM) and the Paciﬁc–South American (PSA) pattern. During austral spring, while the PSA maintained a strong correlation with ENSO throughout the period 1948–2014, the SAM–ENSO correlation changed from being weak before the early 1990s to being strong afterward. Through the ENSO connection, PSA and SAM became more in-phase correlated after the early 1990s. The early 1990s is also the time when ENSO changed from being dominated by the eastern Paciﬁc (EP) type to being dominated by the central Paciﬁc (CP) type. Analyses show that, while the EP ENSO can excite only the PSA, the CP ENSO can excite both the SAM and PSA through tropospheric and stratospheric pathway mechanisms. The more in-phase relationship between SAM and PSA impacted the post-1990s Antarctic climate in at least two aspects: 1) a stronger Antarctic sea ice dipole structure around the Amundsen–Bellingshausen Seas due to intensiﬁed geopotential height anomalies over the region and 2) a shift in the phase relationships of surface air tem- perature anomalies among East Antarctica, West Antarctica, and the Antarctic Peninsula. These ﬁndings imply that ENSO–Antarctic climate relations depend on the dominant ENSO type and that ENSO forcing has become more important to the Antarctic sea ice and surface air temperature variability in the past two decades and will in the coming decades if the dominance of CP ENSO persists.


Introduction
The linkages between the Southern Hemisphere (SH) climate and the tropical Pacific sea surface temperatures (SSTs), particularly those associated with El Niño-Southern Oscillation (ENSO), have been extensively studied (e.g., Karoly 1989;Mo 2000;Yuan 2004;Fogt and Bromwich 2006;L'Heureux and Thompson 2006;Lee et al. 2010). Two leading modes of SH atmospheric variability have been often emphasized: the southern annular mode (SAM; Hartmann and Lo 1998;Thompson and Wallace 2000) and the Pacific-South American (PSA; Mo and Ghil 1987) pattern. The SAM is characterized by outof-phase sea level pressure (SLP) variations between midlatitudes (;408S) and high latitudes (;658S). Previous studies have suggested that the SAM is generated primarily by the internal dynamics of the SH atmosphere, which involves eddy-mean flow interactions (e.g., Yu and Hartmann 1993;Feldstein and Lee 1998;Hartmann and Lo 1998;Limpasuvan and Hartmann 2000;Lorenz and Hartmann 2001). Forcing from the tropics, such as that associated with ENSO, is considered less important, although correlations can be found between SAM and ENSO in austral summer (e.g., Karoly 1989;Seager et al. 2003;Silvestri and Vera 2003;Zhou and Yu 2004;Codron 2005; L'Heureux and Thompson 2006). The PSA is characterized by a stationary Rossby wave train emanating from the tropical central Pacific with major anomaly centers to the east of New Zealand, in the Amundsen-Bellingshausen Seas, and near the southern tip of South America. The positive phase of the PSA is known to be related to La Niña and the negative phase of the PSA is related to El Niño (e.g., Mo 2000). Previous studies showed that ENSO, together with various phases of the SAM and the PSA, can influence the Southern Ocean SSTs (e.g., Ciasto and Thompson 2008;Lee et al. 2010;Yeo and Kim 2015), Antarctic sea ice concentrations (e.g., Liu et al. 2004;Stammerjohn et al. 2008;Yuan and Li 2008), and Antarctic surface air temperatures (e.g., Kwok and Comiso 2002;Ding et al. 2011). Fogt and Bromwich (2006) found that the ENSO influences on the SH strengthened from the 1980s to the 1990s. They noted that the ENSO-excited PSA became more in phase with the SAM in the 1990s to give rise to a stronger ENSO influence. Their study indicates that the relationships among ENSO, PSA, and SAM can change from decade to decade. Interestingly, the early 1990s is the time when ENSO changed from being mostly of the eastern Pacific (EP) type to the central Pacific (CP) type (Yu et al. 2012a(Yu et al. , 2015. This early 1990s change in ENSO properties is a possible reason for the change in ENSO-SH climate relationships. The terms CP ENSO and EP ENSO were introduced by Yu and Kao (2007) and Kao and Yu (2009) to emphasize the different longitudinal positions of these types of ENSO. The EP ENSO is the conventional type of ENSO similar to that portrayed by Rasmusson and Carpenter (1982), which typically onsets and develops near the South American coast and has its SST anomalies centered in the equatorial eastern Pacific. The CP ENSO tends to onset, develop, and decay locally in the equatorial central Pacific. There are indications that the shift in the longitudinal position of ENSO SST variability in recent decades may change the ENSO impacts on the SH climate. For example, Lee et al. (2010) argued that the 2009/10 CP El Niño event caused a record-breaking warming in the South Pacific and western Antarctica associated with a strong anticyclone in the region. Ding et al. (2011) also suggested that the substantial warming in West Antarctica in the recent decades is related to the warming trend of SST in the tropical central Pacific through an enhanced anticyclone over the Amundsen Sea. Lee and McPhaden (2010) suggested that the more frequent occurrence of CP El Niño has contributed to the SST warming trend in the tropical central Pacific. Several recent studies have begun to examine the possible different impacts of the two types of ENSO on the SH climate using observations (Hurwitz et al. 2011;Song et al. 2011;Sun et al. 2013;Ciasto et al. 2015) and model simulations (Zubiaurre and Calvo 2012;Hurwitz et al. 2013;Lim et al. 2013;Wilson et al. 2014).
In this study, we use observational and reanalysis datasets to examine if and specifically how the early 1990s change of ENSO type affects the relationships among ENSO, PSA, and SAM and how these changing relationships impact SH climate variability. We also examine how the two types of ENSO affect the relative importance of two mechanisms suggested to transmit ENSO forcing to high-latitude Southern Hemisphere: 1) an eddy-mean flow interaction mechanism in which an El Niño-induced intensification of the Hadley circulation and equatorward shift in the subtropical jet influence the propagation of midlatitude transient eddies and their associated eddy-mean flow interactions to give rise to a negative phase of the SAM (e.g., Seager et al. 2003) and 2) a stratospheric pathway mechanism in which El Niño events affect the propagation of planetary waves into the stratosphere and induce polar temperature and vortex anomalies that subsequently descend into the troposphere to excite the SAM (e.g., Mechoso et al. 1985;Hurwitz et al. 2011;Son et al. 2013;Evtushevsky et al. 2015). Here, we focus on the austral spring season [September-November (SON)] when ENSO-SH teleconnections are strongest (e.g., Jin and Kirtman 2009), the warming in West Antarctica is most significant (e.g., Schneider et al. 2012;Bromwich et al. 2013), and the maximum stratospheretroposphere coupling occurs (Son et al. 2013).
This paper is organized as follows. Section 2 describes the datasets used in the study. Section 3 examines decadal variations in the relationships among ENSO, PSA, and SAM. Section 4 investigates the linkages between the two SH circulation modes and the two ENSO types. Section 5 discusses the mechanisms enabling the different types of ENSO to influence the SH climate, followed by a discussion of the changing ENSO impacts on SH climate in section 6. Conclusions are presented in section 7.

Datasets
In this study, the following datasets are used: 1) 500-hPa geopotential heights (Z500), zonal and meridional winds, and air temperatures from the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) reanalysis dataset (Kalnay et al. 1996), and 2) the SST and sea ice concentrations from the Hadley Centre Sea Ice and Sea Surface Temperature dataset (HadISST; Rayner et al. 2003). The analyses focused primarily on the period 1948-2014. For sea ice concentrations, however, the data were available only during 1979-2014. Monthly anomalies were obtained by removing the mean seasonal cycles and linear trends.

Decadal variations in the relationships among ENSO, PSA, and SAM
To identify the SAM and PSA, we applied an empirical orthogonal function (EOF) analysis to the covariance matrix of the area-weighted monthly Z500 anomalies over 08-908S following Mo (2000). The two leading EOF modes, which account for 22% and 12% of the total variance, are the SAM and PSA, respectively. The corresponding standardized principal components (PCs) are referred to as the SAM and PSA indices, and their austral spring (SON) values averaged from the monthly PCs are shown in Figs. 1c and 1d. Similar EOFs and PCs were obtained by repeating the EOF analyses using subperiods of the entire data record of 1948-2014 (not shown). In the rest of the study, only the spring SAM and PSA indices are used in the analyses unless stated otherwise. The spatial patterns of these two modes are shown in Figs. 1a and 1b by regressing Z500 anomalies onto the SAM and PSA indices. The positive phase of the SAM (Fig. 1a) is characterized by a zonally symmetric structure that has negative height anomalies over the polar region (;608-908S) and positive anomalies at the midlatitudes (;308-458S) in both the Western and Eastern Hemispheres. In contrast, the positive phase of the PSA (Fig. 1b) is characterized by a zonally asymmetric structure that has a positive height anomaly center to the east of New Zealand (358S, 152.58W), a negative center in the Amundsen-Bellingshausen Seas (608S, 126.58W), and another positive center in the South American-Atlantic sector (658S, 458W). The PSA anomalies are located mostly in the Western Hemisphere and are very small in the Eastern Hemisphere. The SAM and PSA share three anomaly centers (at 47.58S, 1758W; 67.58S, 1258W; and 508S, 608W) in the South Pacific and Atlantic sectors. Various phase relationships between SAM and PSA can result in different impacts on SH climate. For example, an in-phase relationship can result in strong height anomalies over the Amundsen-Bellingshausen Seas by a superposition of the SAM and PSA anomalies.
To investigate the decadal variations in the PSA-SAM relationship, we calculated the running correlation coefficient between the SAM and PSA indices using a 15-yr moving window. As shown in Fig. 1e, the correlation coefficient was small during the 1960s-70s and weakened further in the 1980s, but dramatically increased to large positive values after the early 1990s. The correlation after the early 1990s is statistically significant at the 90% interval using a two-tailed Student's t test. To make sure that the early 1990s increase in correlation coefficient is not a result of the specific method used to construct the SAM and PSA indices, we repeated the running correlation analysis using alternative SAM and PSA indices (hereinafter SAM* and PSA*, respectively) that were defined by other authors. Following Marshall (2003), the SAM* index is defined as the difference of the zonal-mean SLP between 408 and 658S based on observations from 12 weather stations. Following Yuan and Li (2008), the PSA* index is defined as the weighted average of Z500 anomalies at three specific locations (referred to as Z1, Z2, and Z3) in the following way: PSA* 5 (Z1 1 Z2 2 Z3)/3. Here, Z1 is located at 508S, 458W, Z2 is at 458S, 1708W, and Z3 is at 67.58S, 1208W. Figure 1f shows the running correlation coefficient between the SAM* and PSA* indices. The correlation was relatively small in the 1960s-80s but increased to large and significant values after the early 1990s. Both sets of indices show that the PSA and SAM became more in phase in SON after the early 1990s. Such an early 1990s shift in the PSA-SAM correlation was not found in the other seasons (not shown).
We next examined the decadal variations in the relationships between these two SH variability modes and ENSO. We used the cold tongue index (CTI; SST anomalies averaged over 68S-68N and 1808-908W) to represent ENSO. The CTI domain covers both the tropical eastern and central Pacific and, therefore, can reflect the activity of both the EP and CP ENSOs. By superimposing the inverted CTI (green bars) onto the SAM and PSA indices (red lines) in Figs. 1c and 1d, we notice that the SAM index and CTI ( Fig. 1c) do not show an obvious correlation before the early 1990s, but are correlated afterward. As for the PSA index ( Fig. 1d), its variation seems to be more in phase with the variation of the CTI throughout the analysis period. These impressions are confirmed by Fig. 2, which shows the 15-yr running correlation coefficients between the CTI and the SAM and PSA indices. The CTI correlation with SAM ( Fig. 2a) was weak during the 1960s-80s but significant (correlation coefficients around 20.5) after the early 1990s. The CTI correlation with PSA ( Fig. 2b) was significant throughout the analysis period (correlation coefficients around 20.7), consistent with the suggestion that the PSA is excited primarily by ENSO forcing (e.g., Mo 2000). These analyses indicate that the ENSO forcing excited the PSA throughout the past seven decades, but did not become important for SAM until the early 1990s. The in-phase relationship between SAM and PSA in recent decades is likely a result of the increasing contribution of the ENSO forcing to the SAM after the early 1990s. It should be pointed out that the third EOF mode of the Z500 anomalies (not shown) is characterized by a wave train pattern similar to the PSA except for a quadrature phase shift. This EOF mode is often referred to as the PSA2 (Mo 2000;Mo and Paegle 2001) and was suggested to be also excited by tropical convection (e.g., Mo and Higgins 1998;Mo and Paegle 2001) or as a mode of midlatitude intrinsic variability that resembles the PSA pattern (Robertson and Mechoso 2003;Zamboni et al. 2012). We repeated the correlation analyses presented in Figs. 1e and 2b with the PSA2 (not shown) and found the correlation coefficients failed to exceed the 90% confidence intervals. Therefore, the PSA2 is less associated with the SAM and ENSO and is not considered further in this study. ENSO changed from the EP type to the CP type around the early 1990s (Yu et al. 2012a(Yu et al. , 2015. To explore the influence of this change on the ENSO-PSA-SAM relationships, we examined the correlation between tropical SST anomalies with the spring SAM or PSA index for the period before the early 1990s (i.e., 1973-93; pre-1993 period hereinafter) and for the period after the early 1990s (1994-2014; post-1993 period hereinafter). During the pre-1993 period (Fig. 3a), the spring SAM index was not correlated with the tropical SST anomalies, which is consistent with the weak CTI-SAM correlations for that period (Fig. 2a). However, during the post-1993 period (Fig. 3b), negative correlations (90% confidence interval) occur in the tropical central Pacific from May of the El Niño developing year to the following February. The correlations appear first in the northeastern subtropical Pacific in late fall (i.e., May) and later spread toward the tropical central Pacific in winter (i.e., June-August). The strong correlations in the central Pacific grow and persist in the following spring and summer seasons (i.e., September-February). These SST correlations resemble the spatial structure and temporal evolution of the CP ENSO described in Kao and Yu (2009). For the PSA (Figs. 3c and 3d), significant correlations are found in the tropical eastern and central Pacific in both the pre-and post-1993 periods. However, there are noticeable differences between these two periods in the temporal evolution of the correlations. During the pre-1993 period (Fig. 3c), negative SST correlations appear near the South American coast in fall, extend westward to the tropical central Pacific in winter, and develop and peak in the central-toeastern Pacific in the following spring and summer. These SST correlations resemble the spatial structure and temporal evolution of SST anomalies during the EP ENSO described in Kao and Yu (2009). During the post-1993 period (Fig. 3d), negative SST correlations appear off Baja California in fall, spread toward the equatorial central Pacific in winter, and grow and peak in the centralto-eastern Pacific in spring and summer. This SST correlation pattern resembles the evolution of SST anomalies during the CP ENSO. Figure 3 indicates that the SAM is correlated only with the CP ENSO, whereas the PSA is correlated with both the EP ENSO and the CP ENSO. The more frequently occurring CP ENSO after the early 1990s excites both the PSA and SAM to give rise to the in-phase PSA-SAM correlation during the most recent two decades. Based on this result, we hypothesize that the change of ENSO from the EP type to the CP type caused the early 1990s change in the ENSO-PSA-SAM relationships.
To further identify the SAM and PSA signals in the composites, we separated the composite Z500 anomalies into zonally symmetric (i.e., zonal-mean) and zonally asymmetric (i.e., zonal-mean removed) components. Figure 5a compares the zonally symmetric components from the CP and EP El Niño composites with the zonalmean component from the negative phase of the SAM shown in Fig. 1a. The negative SAM (black line) shows a meridional dipole structure with a positive height anomaly centered at high latitudes (;758S) and a negative height anomaly centered at midlatitudes (;458S). The CP El Niño composite (blue line) exhibits a dipole structure with the centers located at similar latitudes of the SAM. Although the EP El Niño composite (red line) also exhibits a dipolar structure, its two anomaly centers do not match those of the SAM. In addition, the zonally symmetric component is stronger in the CP El Niño composite than in the EP El Niño composite, despite the fact that the EP El Niño is typically stronger than the CP El Niño in terms of their tropical Pacific SST anomalies (see Figs. 4a and 4b). As for the zonally asymmetric components, although the CP El Niño composite shows a strong zonally symmetric structure (see Fig. 4c), its zonally asymmetric component is similar to the PSA pattern (Fig. 5b). The EP El Niño composite (Fig. 5c) exhibits an anomaly pattern in the South Pacific and Atlantic sectors similar to the PSA anomaly pattern shown in Fig

Mechanisms linking the two types of ENSO to SAM
Why does the CP ENSO induce the SAM while the EP ENSO does not? Two mechanisms have been suggested in previous studies to link ENSO to the SAM: an eddy-mean flow interaction mechanism in the troposphere and a stratospheric pathway mechanism. The eddy-mean flow interaction mechanism requires that, during El Niño events, the increased meridional temperature gradient associated with a tropical tropospheric warming (Chiang and Sobel 2002) produce a strengthened and equatorward-displaced subtropical jet. The strengthened subtropical jet then affects the propagation of transient eddies to change the location of eddy momentum flux convergence. The changed eddy momentum flux convergence then accelerates westerly anomalies on the equatorward side of the midlatitude jet and accelerates easterly anomalies on the poleward side to give rise to a negative phase of the SAM (e.g., Seager et al. 2003). To examine these processes, we show the vertical structures of the climatology of the zonal-mean zonal wind in SON (Fig. 6a) and the corresponding anomalies composited for the CP and EP El Niño events (Figs. 6b and 6c). Both composites show significant westerly anomalies over the subtropical regions (i.e., 208-308S), which reflect a strengthened and equatorward-displaced subtropical jet (centered at 308S; see Fig. 6a). However, over middle-to-high latitudes (i.e., 308-708S), these two composites are very different. While the EP El Niño composite (Fig. 6c) shows weak and nonsignificant variations, the CP El Niño composite (Fig. 6b) produces a negative phase-like SAM response that is characterized by positive wind anomalies (over 308-458S) equatorward of the midlatitude jet (centered at 508S; see Fig. 6a) and negative wind anomalies (over 508-708S) poleward of the jet. The midlatitude jet is known to be eddy driven, with a central location that coincides with the center of the climatological transient eddy momentum flux convergence (shown in Fig. 6d). We then examine the anomalous transient eddy momentum flux convergence to see if the eddy-mean flow interaction mechanism is at work. Both the CP El Niño (Fig. 6e) and EP El Niño (Fig. 6f)  coincide with the locations of the westerly and easterly wind anomaly bands for the CP El Niño (Fig. 6b) but not for the EP El Niño (Fig. 6c). The results of this analysis indicate that the transient eddy momentum flux convergence and the associated eddy-mean flow interaction can induce the SAM pattern during the CP ENSO but not during the EP ENSO.
To understand why the transient eddy momentum flux convergence cannot induce the SAM for the EP ENSO, we examine the 300-hPa zonal wind in Fig. 7. From this figure (left panels), we notice that the CP El Niño and EP El Niño composites are similar in producing a stationary wave train emanating from the tropical Pacific that represents the PSA. However, the two composites differ over the Indian Ocean sector, where a stationary wave train pattern emanating from the tropics can be found in the EP El Niño composite but not in the CP El Niño composite. This wave train pattern is associated with the Indian Ocean SST anomalies induced by the EP El Niño (see Fig. 4b), which are characterized by an east-west SST anomaly dipole resembling the Indian Ocean dipole (IOD; Saji et al. 1999;Webster et al. 1999). The IOD is known to peak in SON and be capable of exciting a wave train that propagates toward high latitudes (Lau and Nath 2004;Yu and Lau 2005;Cai et al. 2011). We also analyze the momentum flux convergence of the stationary waves in the EP El Niño composite and display its zonal-mean values in Fig. 7b (right panel). It is interesting to note the zonal means (blue line) tend to be of opposite sign to those of the momentum flux convergence produced by transient eddies (red line). As a result of this cancellation, the eddy-mean flow interaction mechanism is not able to produce a SAM pattern during EP El Niño events. In contrast, during the CP El Niño (Fig. 7a), the zonal means of the stationary eddy momentum flux convergence are relatively weak compared to the transient eddy momentum flux convergence. The anomalous transient eddy momentum flux associated with the CP El Niño produces positive and negative wind anomalies, respectively, on the equatorward (shaded in light gray) and poleward (shaded in dark gray) sides of the midlatitude jet over most of the SH, except over the Pacific, to give rise to a SAM pattern. Figure 7 indicates that the different influences of the two types of ENSO on the Indian Ocean SST anomalies enable these two ENSO types to produce difference influences on the SAM through eddy-mean flow interactions. It is noted that the value of the stationary eddy momentum flux convergence is weakly positive at the equator in our Fig. 7 (right panels), which is different from the weakly negative value reported in Seager et al. (2003). We were able to determine that the reason for this difference was the different analysis methods used in these two studies (not shown). We next examine the stratospheric pathway mechanism. This mechanism suggests that ENSO can induce planetary waves leading to a stratospheric polar warming (Hurwitz et al. 2011(Hurwitz et al. , 2013Zubiaurre and Calvo 2012). The associated longitudinal temperature variations weaken the polar vortex, whose associated circulation anomalies can later descend into the troposphere to give rise to the SAM (Son et al. 2013;Evtushevsky et al. 2015). The evolution of the zonal wind and temperature anomalies composited for the two types of El Niño during the spring season are shown in Fig. 8. The CP El Niño composite exhibits a weakening of the polar vortex (i.e., the negative wind anomalies over the region above 16 km and poleward of 608S; Fig. 8a) in September, which is accompanied by a stratospheric warming centered at 20 km in the polar region (Fig. 8b). The wind and temperature anomalies in the stratosphere then descend to the troposphere in the following months (October-November) to give rise to a negative phase of SAM with easterly anomalies around 608S and westerly anomalies around 408S (see Fig. 8a). In contrast, the EP El Niño composite does not show a weakening in the westerlies in the polar stratosphere nor a warming there (Figs. 8c and  8d). To confirm that the composite result is not overly influenced by any particular El Niño event, we repeated the analysis with the 2002 El Niño event excluded (this particular event is known to produce a massive Antarctic stratospheric sudden warming; Newman and Nash 2005) and found similar results (not shown). The analyses we presented in this section suggest that both the eddy-mean flow interaction and the stratospheric pathway mechanisms enable the CP ENSO to influence SAM variability in the SON season.

The changing impacts of ENSO on SH climate
In this section, we investigate how the post-1993 ENSO-PSA-SAM relationships impacted the Antarctic climate. We first examine the intensities of SAM and PSA using the 15-yr running variance of these two indices. The SAM variance increased after the early 1990s (red line in Fig. 9a), which is likely a manifestation of the additional contribution from the CP ENSO forcing to the SAM variability. The PSA variance does not show a substantial change (blue line in Fig. 9a), implying that it is not affected by the change of ENSO type in the early 1990s.
As mentioned, the SAM and PSA share a common anomaly center over the Amundsen-Bellingshausen Seas (ABS). This region is one of the best locations to detect the impact of the changing ENSO-PSA-SAM relationships (or the changing ENSO type) on the SH climate. We defined an ABS index as the Z500 anomalies averaged over 558-758S, 1008-1408W (Figs. 9c and 9d). A 15-yr running variance of this index indicates an increase in the ABS intensity after the early 1990s (Fig. 9b). Correlating the ABS index with the SAM and PSA indices shows that 92% of the total variance of the ABS index is explained by the SAM and PSA indices. Comparing the Z500 anomalies regressed onto the spring ABS index in the pre-1993 period and in the post-1993 period shows that the positive ABS height anomalies intensified during the post-1993 period (Figs. 9c and 9d). Therefore, the strengthened ABS height anomalies in recent decades can be explained as a result of both the intensified SAM and the increasingly in-phase relationship between the SAM and PSA.
The ABS variability has profound impacts on sea ice concentrations (SICs) in the Southern Ocean (e.g., Turner et al. 2013). Geopotential height variations over the ABS region can induce anomalous eddy heat and momentum fluxes to produce the so-called Antarctic sea ice dipole structure (ADP; Yuan and Martinson 2001), which is characterized by out-of-phase SIC anomalies between the Pacific and Atlantic sectors. We examine the variation of the spring ABS index from 1979 to 2014 and note the amplitude of the variation increased after the early 1990s (Fig. 10a, top panel). The ADP can be identified from the longitudinal distribution of the SIC anomalies over the 608-708S belt in early summer [November-January (NDJ)] as the out-of-phase anomalies between the Pacific sector over 1308-1508W and the Atlantic sector over 408-608W (Fig. 10a, bottom panel). Consistent with the increasing ABS index, the amplitude of the ADP also increased after the early 1990s. Following Yuan (2004), we compute the ADP index as the difference between the SIC anomalies averaged in a South Pacific box (608-708S, 1308-1508W) and an Atlantic box (608-708S, 408-608W). Yuan and Li (2008) suggested that the correlations between atmospheric circulation modes and sea ice anomalies are the strongest when the ice lags the circulation modes by 2 months. We regressed SIC anomalies onto the ADP index in the NDJ season separately for the pre-1993pre- and post-1993pre- periods. During the pre-1993pre- period (1979 Fig. 10b), positive SIC anomalies occur in the Pacific sector over the Ross Sea and negative anomalies occur in the Atlantic sector over west of Weddell Sea. During the post-1993 period (1994-2014; Fig. 10c), the dipole structure is similar to that in the previous period. However, the ADP has stronger amplitudes in the post-1993 period than in the pre-1993 period. Our analysis indicates that the sea ice dipolar structure induced by ENSO becomes stronger after ENSO changed from the EP to the CP type in the early 1990s.
We next examined the possible impact of the changing ENSO-PSA-SAM relationships on Antarctic surface air temperatures. We applied an EOF analysis to find the leading variability modes in spring surface air temperature (SAT) anomalies over 308-908S. The structures of the first two leading EOF modes (EOF1 and EOF2, respectively) are shown in Figs. 11a and 11b. Both modes exhibit a dipolar structure with a temperature anomaly center over the Weddell Sea that is out of phase with an anomaly center over the Ross Sea. The two EOFs are different in that surface temperature anomalies over East Antarctica are in phase with the temperature anomaly center over the Ross Sea in EOF1 but are in phase with the anomaly center over the Weddell Sea in the EOF2. EOF1 and EOF2 together suggest that surface air temperature variations over Antarctica are largest in three regions: East Antarctica, West Antarctica (related to the anomaly center over the Ross Sea), and the Antarctica Peninsula (related to the anomaly center over the Weddell Sea). Associated with the anomaly dipole in EOF1 and EOF2, surface temperature variations over the Antarctic Peninsula and West Antarctica tend to be out of phase when either one of the EOF modes dominates. However, surface temperature variations over East Antarctica can be in phase with those over West Antarctica when EOF1 dominates, or in phase with those over the Antarctic Peninsula when EOF2 dominates. When the amplitudes of EOF1 and EOF2 are comparable, their temperature anomalies over East Antarctica tend to cancel each other out. In that case, temperature variations in East Antarctica would be not correlated with either those in West Antarctica or those in the Antarctic Peninsula.
To examine the linkages between these two EOF modes and the SAM and PSA, we regressed surface air temperature anomalies onto the inverted SAM and PSA indices. It is interesting that the SAM regression pattern resembles EOF1 whereas the PSA regression pattern resembles EOF2 (Figs. 11c and 11d). The correlation coefficient between the SAM index and principal component of EOF1 is 20.79, indicating that EOF1 is more associated with the SAM. The correlation coefficient between the PSA and principal component of EOF2 is 20.47, indicating that EOF2 is more associated with the PSA. Based on these linkages, we expect the EP El Niño influence on Antarctic surface air temperature to be close to that of EOF2 because the EP El Niño can excite a negative phase of PSA. Figure 11f shows the surface air temperature anomalies composited from the EP El Niño events. As expected, the EP El Niño composite is similar to EOF2, with the temperature anomalies over the Antarctic Peninsula out of phase with the anomalies over West Antarctica and in phase with the anomalies over East Antarctica. The CP El Niño composite for surface temperature anomalies is shown in Fig. 11e. In this composite, significant temperature anomalies over the Antarctica are found only over West Antarctica and the Antarctic Peninsula, with their anomalies out of phase. Temperature anomalies over East Antarctica are weak and are mostly not statistically significant. These two features can be explained by our finding that the CP El Niño can excite both a negative phase of PSA and a negative phase of SAM. As a result, the CP El Niño composite should resemble the combined pattern of EOF1 and EOF2, which tends to cancel out surface temperature anomalies over East Antarctica and at the same time intensify the out-ofphase relationships between the temperature anomalies over the Antarctic Peninsula and West Antarctica.

Conclusions
In this study, we examined the decadal variations in the relationships among ENSO, the SAM, and the PSA to understand how the ENSO influence on the high-latitude Southern Hemisphere has varied during the past seven decades. It is our conclusion that the change in ENSO type in the early 1990s altered the mechanisms and spatial patterns of ENSO influence on Southern Hemisphere climate. Before the early 1990s, ENSO was dominated by the EP type, which influenced the Southern Hemisphere climate via an excitation of the PSA wave train pattern. The ENSO-Southern Hemisphere influence was strong in some events but weak in others, depending on the random phase relationship between the PSA and the internally generated SAM. After the early 1990s, the dominant CP ENSO modulated not only the PSA but also the SAM. As a result, the SAM and PSA events became more in phase after the early 1990s. The increasing influence of the CP ENSO on the SAM is accomplished through both an eddy-mean flow interaction mechanism and a stratospheric pathway mechanism.
The changing ENSO-PSA-SAM relationships have rendered ENSO more effective in influencing the Southern Hemisphere after the early 1990s. Two such impacts were identified in this study. One of them is related to the sea ice dipole pattern around the Antarctic Peninsula, and the other is related to the phase correlations of the surface air temperature variations in East Antarctica, West Antarctica, and the Antarctic Peninsula. While the sea ice dipole has been known to be related to the ENSO forcing (Yuan 2004), we find the dipole to be stronger after ENSO changed to the CP type. The in-phase relationship between the PSA and SAM enhanced the geopotential height anomalies over the Amundsen-Bellingshausen Seas. The different ENSO types also have different footprints on Antarctic surface air temperature patterns. The EP ENSO tends to link surface temperature variations across the entire Antarctic continent, while the CP ENSO tends to link the temperature variations between West Antarctica and the Antarctic Peninsula only.
The findings in this study can help explain the changing role of ENSO in decadal and interannual variability in high-latitude Southern Hemisphere climate. Such changes have implications for sea ice, ocean circulation, biogeochemistry, and ecology. The relationships found here among ENSO, SAM, and PSA could potentially be used in conjunction with Antarctic paleo-temperature records to examine changes in ENSO type on longer time scales. This study does not examine the reasons why the early 1990s change in the PSA-SAM relationship is observed only during the SON season. Pursuit of this question will

9406
require further analysis of the other three seasons and is beyond the scope of this study. Nevertheless, it is possibly related to the seasonal differences in the SH atmospheric circulation, such as the strength and location of the subtropical and eddy-driven jet streams, which support different wave train and transient eddy responses. It should also be noted that this study is an empirical study and that modeling studies are needed to reach more definite conclusions on the causality of the early 1990s change.