Received 25 December 2015; accepted 21 February 2016; published 24 February 2016
The atmospheric circulation in the Southern Hemisphere (SH) has a large impact on SH climate and sea ice cover through complex air-ocean-sea ice interaction. Previous studies have identified the hemispheric-scale atmospheric variability in SH comprising of a dominant mode characterized by zonally symmetric see-saw pattern between the mid- and high-latitudes. This mode appears as the leading mode of many atmospheric variables, e.g., mean sea level pressure (MSLP), geopotential heights, zonal winds  - , which has been referred to as Antarctic Oscillation (AAO) or Southern Annular Mode (SAM) and dominates the region poleward of 20˚S     . Second and third modes are identified as the El Niño-Southern Oscillation (ENSO)-driven Pacific- South America (PSA) teleconnection patterns  . These modes are part of a stationary wave train generated by tropical convection that travels towards Polar regions.
Aforementioned studies have focused on the role of global climate indices and circulation cells on Antarctic climate and sea ice. However understanding of Southern Polar Cell’s impact on sea ice system has been fragmentary. Though the impact of Ferrel Cell variability is well documented, as Marshall  provided station- based SAM time-series as a see-saw pattern in Ferrel Cell which circulated in the region bounded by mid-lati- tudes (~30˚S) in the north and Antarctic coast (~65˚S) in the south, there is a significant knowledge gap in the understanding on how the local circulation cell―the Polar Cell―modulates seasonal and interannual Antarctic sea-ice variability. In this study, we focus on the impact of the local atmospheric circulation pattern due to its immediate interaction with and in vicinity of sea ice. Moreover, several studies have suggested that local climate and circulation might have a larger role in modulating the sea ice variability  -  .
This study attempts to identify the modes of variability in the Polar Cell (region polewards of 65˚S). We perform empirical orthogonal function (EOF) analysis on geopotential height anomalies south of 65˚S, in order to highlight the dominant circulation patterns over Antarctic continent and establish its relationship with the observed sea ice changes (e.g.,   ). The Polar Cell is characterized by intense surface divergence at the South Pole (SP) and a convective zone at the Antarctic coast. Hence, a strong pressure gradient forms between high-altitude Antarctic plateau and Antarctic coastal region that drives air circulation in the Polar Cell. However, previous studies have investigated air mass movement in and out of Polar Cell, which is known as the semiannual oscillation (SAO;   ) that is characterised by bi-annual transfer of air masses between the Antarctic and mid-latitudes  . It is manifested by the changes in the position and strength of the Antarctic circumpolar trough (CPT), lying between 60˚ - 65˚S, which has major impact on Antarctic climate and sea ice.
In SAO studies, the choice of latitude bands (50˚ - 65˚S) selected for its computation does not capture the full variability of forces which drives intra-cell circulation in the Polar Cell. Hence, we propose a new polar-coastal pressure gradient index called the Polar Coastal Index (PCI) which is defined as normalized MSLP difference between 65˚S and 90˚S. It is similar to the SAO, but the important difference between PCI and SAO is that the former accounts for strong pressure gradient between SP and Antarctic coast which drives most of the Polar Cell circulation and the CPT, whereas the SAO measures air masses movement between 50˚S and 65˚S (peripheral ocean).
We employ reanalysis products (description of all datasets in Section 2) and use EOF techniques to study the atmospheric variability. As would be shown in results (Section 3), the PCI could explain more than 15% variability observed in the CPT and about 30% of the Amundsen-Bellingshausen Sea Low (ABSL) variability  . Essentially PCI is a measure of the CPT variability and hence coastal winds. Studies have shown that these surface wind fields   and near-shore winds  play a major role in the convective processes and exchange of heat and momentum with ocean. We attempt to investigate modes of variability in the Polar Cell and its relationship with the above mentioned global climate mode indices. We will characterise dominant modes of Polar Cell in relationship with Antarctic sea ice concentration. We employ lead-lag correlation analysis for each climate indices with sea ice, in order to explore causation impact on each other. As dominant modes are extracted from geopotential height anomalies and PCI is generated from MSLP field this procedure will provide us with more information on variability of the Polar Cell across atmospheric levels. We discuss the findings and offer conclusions in section 4.
2. Data and Methods
The sea ice concentration (SIC) with a spatial resolution of 25 × 25 km was obtained from the National Snow and Ice Data Center. The monthly SIC spanning 1979-2013 were computed by using the NASA Team algorithm  which were obtained from SMMR, SSMI, and SSMIS observations. The data used to examine circulation features are monthly mean geopotential heights derived from National Centers for Environmental Prediction- Department of Energy (NCEP-DOE) Reanalysis 2 (NCEP2) dataset  . It is an updated and improved version of the widely used NCEP-NCAR data, assimilated with global satellite soundings network data. The horizontal resolution is 2.5˚ × 2.5˚ and the study period spans 1979 to 2013.
To capture the essential features of Polar Cell, we apply EOF or principal component (PC) analysis to separate observed variability into orthogonal modes, each consists of a spatial pattern and a time series. Before carrying out the EOF analysis the data was processed as follows. First, annual cycle was removed from the datasets by subtracting climatological monthly means from individual months. Second, to ensure equal area weighting for the covariance matrix, the gridded data was weighted by the square root of the cosine of latitude to obtain the EOF loading patterns. The modes are arranged in the descending order of variance explained by them, i.e., the first PC explains more variance than second and so on. The primary modes of variability in the Polar Cell is defined as the first three leading modes of monthly height anomalies at 500 hPa (Z500), 800 hPa (Z800), MSLP in region poleward of 65˚S. Monthly time series of PCs are constructed by projecting the monthly mean anomalies onto the modes. The normalized PC time series have the same units as the anomaly field hence dimensionless.
The variability in the Polar Cell can be extracted by considering two extremes of its circulation, thus the PCI pattern is most profound feature in MSLP field, e.g. average MSLP values over continental Antarctica and coastal region differ by 37.5 hPa over the study period. Monthly PCI was computed from NCEP2 by subtracting normalised MSLP at 65˚S (Antarctic coast) from 90˚S (SP) (Figure 1(d)). The ENSO/Southern Oscillation Index (SOI) was obtained from NOAA (http://www.cpc.ncep.noaa.gov/data/indices/soi) and the SAM time series was obtained from British Antarctic Survey (http://www.nerc-bas.ac.uk/icd/gjma/sam.html).
The time series of these EOF modes and PCI were standardized and lines of least squares fit were fitted to quantify trends. The trend over 34 years was determined as the difference between the values at the two endpoints of the least-squares-fit line divided by the number of years. The long-term trend was removed at each SIC grid point and from each of the indices before computing the correlation coefficients which were judged at 99% significance level by using the standard Student t-test. To focus on interannual and longer timescales, we filtered the monthly time series of climate indices and sea ice anomalies with a low pass filter, an infinite impulse response (IIR) filter with cut-off threshold at 13-month. All time-series are filtered unless otherwise stated.
3.1. Modes of Variability in SH Polar Cell
The spatial pattern of first three leading modes of MSLP, Z800 and Z500 are very similar (Figure S1), and also their time series show high correlation (Table 1). Total variance explained by first three modes of each variable MSLP, Z800 and Z500 are 70%, 78% and 79%, respectively. It is apparent that Polar Cell is barotropic, hence we will use modes of Z500 field to define Polar Cell variability as it explains more variability than any other
Table 1. Cross-correlation between PC’s of Z500, MSLP and Z800, all are significant at 99% level (shared variance in percentages).
field. Spatial pattern of first three PCs of Z500 covers important features contained in Polar Cell (Figure 1), such as a high-pressure over Antarctic continent and a low pressure zone over Amundsen Sea and Antarctic Peninsula (AP).
Leading PC (PC1500) pattern which explains more than 50% of the variability shows anomalous negative values all over Antarctica with three main lobes near coast of Amundsen Sea and high-altitude East Antarctic plateau region. High pressure zone is expected in this region due to permanent cold temperatures, high-altitude and Polar divergence. This quasi-permanent anomaly affects circulation over entire continental Antarctica and divergence zone results in lowest recorded precipitation in these regions (Figure S2(a)).
The second dominant mode (PC2500) explains ~14% of variability and shows a pattern of low pressure zone, as positive anomalous heights in the Amundsen and Bellingshaushen Seas (ABS) (similar to ABS Low) and negative anomalies in the East Antarctica. Reference  pointed out that the low in ABS exists because of the asymmetrical topography of AP and Antarctic continent which actively influence the air circulation in the region. It’s also been suggested that ABSL strength is influenced by large-scale patterns of climate variability, namely SAM  and ENSO   . However, the interaction between these climate modes and topographic conditions, which dictates strength and position of the ABSL, are not yet clear.
The third dominant mode (PC3500), explains ~10% variability, with anomalous positive values moving to AP region. Reference  investigated large number of cyclogenesis events which take place within this low pressure zone (around AP) between the latitudes of 60˚ and 70˚S. They pointed out that the main mechanism behind
Figure 1. Spatial pattern of first three leading modes of Z500 (a) PC1500 (b) PC2500 (c) PC3500 and (d) latitude bands used to calculate PCI (red lines), mean sea ice edge (solid blue) and ±1σ (dashed blue).
meso-cyclones developments were the lows, which are formed inside pre-existing low pressure zone in the AP region, from which synoptic-scale cyclonic weather systems emerges leading to large cyclone density. This cyclonic activity due to deeper lows has significant impact on local climate and sea-ice conditions around AP.
Lastly, the PCI, which describes essentially zonal mode in the Antarctic coastal region, is characterized by the bi-annual heightened meridional pressure gradient similar to SAO  -  which is due to difference in heating of the Antarctic continent and the coastal region. This influence is also visible in the spatial pattern of PCI regression on MSLP field, in which the PCI accounts for more than 15% variability in much of the Antarctic coastal region (e.g., CPT), except for AP. Since, more than 30% variability can be accounted in Amundsen Sea and Ross Sea regions, we suggest that PCI can be used as a reliable proxy to measure variability in the CPT and especially for ABSL phenomena (Figure 2(a)). The low percentage of the variance explained by PCI over the Antarctic continent implies that PC1500 or combination of higher order modes is likely dominant there. Therefore, PCI can vary across the latitudes around the region that it affects, both in terms of strength and extent. Seasonal variations in these modes and PCI are discussed below.
3.2. Seasonality of Modes
From January till onset of winter season (Jun-Jul-Aug), PC1500 weakens (Figure 2(b)), as Antarctic continent is heated up due to large incoming solar radiation, pressures over continent are lower than average; it only peaks in
Figure 2. (a) Spatial distribution of MSLP field variance explained by PCI in percentages (%), and (b) Seasonal variations of first three leading modes and PCI.
winter season when extremely low temperatures return. PC2500 follows half-yearly cycle, which is mainly due to fall and rise in the strength of ABSL. As ABSL strengthens in autumn (Mar-Apr-May) and spring (Sep-Oct-Nov) PC2500 follows it closely. However, PC3500 shows an asymmetric cycle―it peaks in March and in winter months (Jul-Aug), and dips in May and September between these two peaks. It can be explained as follows. After period of large solar radiance at the end of summer season, AP region is heated hence PC3500 index is positive. Again when winter sets in, the CPT contracts and a low pressure zone is created which brings AP region under the influence of westerlies, which is mild maritime air from the north or north-west as found by  .
The PCI is overwhelmingly dominated by half-year cycle like SAO depicting two peaks and two troughs. The PCI represents a bi-annual spread and tightening of CPT, which exhibits a minimum pressure and extends furthest south in the autumn and spring seasons (Figure 2(b)). This pattern occurs due to difference in the heat capacity of the Antarctic and the Southern Ocean; the shift in air-masses described by PCI strives to equalise the intra-cell imbalance. The PCI is associated with the sub-polar convergence zone which exerts influence over the stages of sea ice advance and retreat cycles   . The changes in PCI will have a direct impact on the chronological phases of sea ice formation and melting.
3.3. Trends in Climate Mode Indices
Significant trends exist only in two climate mode indices (PCI and PC2500). Apparently, the trend in PCI (0.21 per decade, p < 0.01) is the largest, followed by PC2500 (−0.16 per decade, p < 0.01) (Figure 3(a) and Figure 3(b)). The positive trend in PCI could mean a deeper CPT and/or stronger high pressure at SP. Our results concurs with the finding of reference  who pointed out that MSLP has decreased all around the CPT. Moreover, the cooling of East Antarctica due to loss of stratospheric ozone  has also been attributed to changes in SAM. Reference  pointed out that cooling could be due to alterations to the circumpolar vortex. All these studies point to favourable conditions in which PCI will have positive trend.
A strong trend in PC2500 also indicates an increasingly lower pressure over ABS (e.g., ABSL) and higher pressure over East Antarctic elevated plateau. Proposed mechanisms for Antarctic cooling has been mentioned above and also Reference  pointed out that over 1979-2008 period, the ABSL has deepened in spring and autumn. We conclude that simultaneous deepening of ABSL and cooling over Antarctic plateau has resulted in the observed PC2500 trend. Overall, simultaneous warming in the AP and cooling over Antarctic plateau is related to broader changes in atmospheric circulation (e.g., SAM) and changes in sea ice   .
Both PCI and PC2500 trends could result from ever more coherent ABSL response to ENSO cold events   . This is expected as both modes draw their existence from CPT in case of former and ABSL for later. It can be argued that deepening of CPT (or ABSL) will provide reinforcement for these two modes (PC2500 and PCI).
3.4. Is Variability in Polar Cell Related to Other Climatic Modes?
To investigate inter-relationship between known climate modes of SH and variability in Polar Cell, we use correlation analysis. Figure 3(a) shows that the correlation between SAM and PC1500 is quite strong (R = 0.69, p < 0.01), indeed PC1500 is a manifestation of SAM in the Polar Cell. This result extends observation made by previous studies that, SAM is a dominant mode of variability in the SH including Polar Cell. It is interesting to note that PC1500 and PCI (R = 0.32, p < 0.01) are not strongly correlated suggesting that both represent independent variability in the Polar Cell, which can be attributed to the difference in spatial structure of PC1500 and PCI (Figure 1(a), Figure 1(d)). Also, SAM is highly correlated with PCI (R = 0.62, p < 0.01). The PCI and SAM describe similar zonal modes of pressure gradients in high latitudes. PCI is a measure of the pressure gradient between 90˚S and 65˚S, while SAM represents pressure difference between 40˚S and 65˚S; so both indices share one region in common. Both indices owe their existence to movement of air masses, back and forth from mid-latitude to coastal Antarctica in the case of SAM, and from coastal Antarctica to elevated Antarctic region in the case of PCI. Since, the fluctuations in pressure differences in Antarctic coastal region and mid-latitudes (SAM’s variability) also modify the pressure distribution that is used to define PCI, so high correlation is expected.
There are differences between the two indices as well. PCI is a semiannual phenomenon whereas SAM has an annual cycle. The SAM influences the strength of the westerlies, while PCI drives katabatic winds in the coastal region, although it also contains zonal variability. Since, SAM and PCI share 38% of observed variance, it further proves our point that PC1500 and PCI are both driven by and the variability in Polar Cell is influenced by
Figure 3. (a) Time-series of SAM, PC1500 and PCI, with PCI trend line (in olive green); (b) Time-series of SOI and PC2500, with PC2500 trend line (in pink).
SAM. Other PCs do not show any significant relationship with SAM.
We also seek to highlight the relationship between variability in Polar Cell and ENSO forcing. The PC2500 is negatively correlated (R = −0.66, p < 0.01) with standardised SOI. Previous studies have also detected positive height anomalies at 500 hPa level over the ABS in response to ENSO events   . Overall, PC2500 and SOI complement each other closely (Figure 3(b)). This is consistent with findings of Mo  that second PC (e.g. PC2500) is the SH (Polar Cell) is response to warm ENSO events and correlates well with the negative SOI. They pointed out that the ENSO events are primary mechanisms for excitation of second mode of variability in the SH atmosphere
These modes seldom act independently despite the fact that modes are driven by diverse processes and computed from different climate variables. For example, reference  found that when an ENSO cold event takes place with the positive SAM, ABSL is intensified. Similarly, the ABS region experiences anomalous high pressure when ENSO warm event is concurrent with the negative phase of the SAM. And these changes in ABSL are shown to affect West Antarctic climate and sea ice around Antarctica  .
3.5. Impact of Polar Cell Modes on Sea Ice
We evaluate the relevance of each index by the strength of correlation of these indices with the SIC pixels. We quantify the influences from these climate indices by calculating number significant correlations with indices lag and lead sea ice by 6 months, in order to explore causative impact on each other. Monthly time series were low pass filtered to focus on longer than sub-annual timescales.
Figure 4 shows number of grid points for correlation coefficient significant at p < 0.01. The total number of grids was 104,912 (332 × 316) out of which sea ice grids were 82,907. With modes lagging SIC by 1 to 6 months, SAM is less influential than PC3500 and PCI, suggesting that the former cannot affect sea ice more than a month in advance (Figure 4). PC1500 and PC2500 modes exhibit more significant correlations than SAM for all lead-lags. SAM is slightly behind PC1500 at all times, while the difference between SAM and PC2500 is as much as 5 times at 6-month mode lag then narrows at 2-month mode lead but still consistently higher than expectation. Above observations clearly show that SAM is relatively less influential on sea ice than other climate modes. Reference  also observed similar results and concluded that local factors have “consistently higher” influence on sea ice than SAM, which confirms that local circulation has major impact on SIC.
Figure 5 shows spatial pattern of peak correlation between these modes and SIC when the former leads the latter by 2 months, except for PCI and PC3500 where they lead SIC by 4 and 12 months, respectively. The PCI and PC1500 have impact on SIC similar to that of SAM (Figure S2(b)), which is expected due to high cross- correlation with SAM. The impact of PCI, PC1500 and PC3500 on SIC is similar with positive correlation around Antarctic coast, except Weddell Sea. In ABS and Ross Sea, we find significant positive correlation, while negative correlation occurs in Weddell Sea (Figure 5(a), Figure 5(b), Figure 5(d)). PC2500 strongly projects the Antarctic Dipole (ADP) SIC anomalies in the western Antarctica (Figure 5(c)). This correlation pattern results from the observed low similar to ABSL (Figure 1(b)), which generates opposite heat fluxes in the Ross Sea and Weddell Sea  .
As ABS is the region with significant variability in Southern Ocean on the interannual time scale, PC2500 represents the largest contribution in ADP regions. The second largest impact comes from the PC1500 pattern. As its spatial loading occurs around East Antarctica, Ross Sea and ABS, we point out that its influence on SIC is evenly distributed among the three lobes (Figure 1(a)). The correlation pattern from the PCI is quite similar to PC2500 in both hemispheres, only the polarity is reversed. The PC3500 mode is formed by anomalous low pressure system which creates milder conditions over AP (Figure 5(d)). In both hemispheres, PC1500 and PC2500 patterns generally have higher correlations with SIC than PCI and PC3500. Though the influences from all modes
Figure 4. Number of significant correlations of each climate mode (PC1500, PC2500, PC3500, PCI and SAM) with 6 months lag and 6 months lead with respect to SIC (on left-hand side modes are lagging while on right-hand side leading the sea-ice).
Figure 5. Spatial correlation of (a) PCI, (b) PC1500, (c) PC2500, and (d) PC3500 with a lead of 4, 2, 2 and 12 months, respectively, with sea ice concentration (SIC); regions with values significant at 99% level are enclosed in green contours.
are present in South Western Pacific Ocean, only PC1500 has significant impact especially off the coast of Wilkes Land.
The highest number of significant correlations occurs at lead of 2 months for PC1500 and PC2500 modes, 12 months for PC3500, and 4 months for PCI, suggesting that SIC reacts to circulation changes after a pause. This analysis delineates cause and effect relationship between the climate indices and SIC. It is noted that these assessments reflect the influence of the each climate mode on SIC, which may not adequately explain SIC trends observed in individual regions. For example, the largest increase in sea ice is in the Ross Sea with minor contributions from the Weddell Sea and Indian Ocean sector, contrasting with decline in the ABS  . These trends cannot be attributed to individual modes as such trends emerge from complex interaction between these modes at different time-scales.
Moreover, our investigation does not point to decadal variability owing to the limited time series. Apart from interannual variability, these climate indices have significant sub-annual variabilities. As aforementioned, PC2500 and PCI have the largest signal in autumn and spring, whereas PC1500 pattern is dominant in winter (Figure 2(b)). To include seasonal variability  , we repeated the analysis with unfiltered time series. The seasonal correlation patterns were similar to Figure 5, but with largely paler magnitude for all modes except for PCI (not shown), suggesting that PCI yields substantial control over sub-annual SIC variability and merit separate study to understand synoptic scale processes.
4. Discussion and Conclusions
In this study, we examined the main characteristics of leading modes of Southern Polar Cell and their relationship with global climate indices along with its impact on Antarctic sea ice. Four climate modes were generated from different atmospheric variables MSLP (PCI) and Z500 (PC1500, PC2500 and PC3500). It was observed that the Polar Cell is barotropic and Z500 field was used for further investigation as it accounted for more observed variability than other variables. First three dominant modes combined explain nearly 80% of the observed variability. Maximum loading of PC1500 occurs at high-plateau regions in the East Antarctica and Marie Byrd Land in the West Antarctica coinciding with high MSLP regions. PC2500 pattern is a well documented feature representing ABSL system forming off the coast of West Antarctica and PC3500 mode is formed by a low pressure centre over AP.
As expected PC1500 is closely related to SAM, which confirms that SAM indeed is dominant atmospheric feature of the SH. Also, PCI is highly correlated with SAM suggesting that latter’s signal is present in all atmospheric variables. The PC2500 was found to share large amount of variance with ENSO confirming Mo’s  and related studies. It is observed that PSA pattern of high pressure in the Amundsen Sea, low pressure in the South Atlantic Ocean and east of New Zealand is formed in response to ENSO warm events and vice versa for cold events. They referred PSA pattern as sum of two distinct modes PSA1 and PSA2 which are high- and low-frequency response to the ENSO signal, respectively. However, a recent study by  has questioned the existence of PSA1 and PSA2 modes as separate modes due to non-compliance to separation criteria of  . If same test is applied to PC2500 and PC3500, the difference in amount of variance explained by these two is not large enough to satisfy such criterion. However, the existence and uniqueness of a climate mode can be judged by amount of impact it has on other fields. Here, we find that PC2500 and PC3500 not only separately exist but are driven by different physical processes, testified by unique impact on sea ice.
One of the interesting results is that only two (PCI and PC2500) climate indices have significant linear trend over the study period. The observed trend results from simultaneous ABSL deepening and cooling over Antarctica. ABSL deepening is in response to La Nina events   . Some studies indicate that Antarctic cooling may have been caused by a combination of positive SAM  , stratospheric ozone depletion  and stratospheric cooling due to CO2  . Seasons of largest changes  in ozone (spring) and wind circulation (summer and autumn) are not synchronised with deepening of the CPT (spring and autumn) suggesting a non-linear relationship. Also PCI and PC2500 show clear seasonality of half-yearly cycle strictly following movement of the CPT (Figure 3).
Based on the correlation patterns produced by these modes (Figure 5), we would anticipate sea ice trends different than current sea ice trends. The correlation maps are very dissimilar to one another and often contradictory as mentioned before. The reported sea ice trends are resultant of complex amalgamation of known and unknown climate modes on different time scales. However, a few points can be noted. All climate modes except PC2500 favour sea ice growth in Indian Ocean with strongest impact from PC1500. As for the West Pacific Ocean sector only PC1500 have significant positive impact, while others modes have mixed influences. Most significant correlations are found in the ADP regions which differentiate one mode from another. All modes except PC2500 intensify sea ice dipole in the Weddell Sea and ABS, similar to the response of ENSO cold events (La Nina). However, PC2500 response is similar to ENSO warm event (El Nino) as revealed by a strong negative correlation between PC2500 and SOI (Figure 2).
Current changes in sea ice convey opposite message compared to what can be expected from correlation pattern for SAM (Figure S2(a)) or PCI/PC1500 modes (Figure 5(a), Figure 5(b)). For example, SAM and PCI would favour positive trend in ABS while negative trend in Weddell Sea, on the contrary the current trend is just reverse of that. We argue that large impact from PC2500 nullifies the effect of other modes and reverses the polarity of the sea ice trends in the above mentioned regions, resulting in observed sea ice trends  . However, the trend is positive for PCI and SAM; and negative for PC2500, hence we suggest that current SIC trend in these regions would reverse, provided polarity of trends of these modes remains same. In addition to the spatial patterns of correlation, Figure 4 reveals that the major air-ice interaction occurs in the forms of proposed climate modes, PC2500, PC1500, PCI and PC3500 in the decreasing order of influence in the interannual time scale.
The interplay of atmosphere and sea ice ensue continuous and delayed response. The continuous response could be due to thermodynamic impact of the surface circulation on the ice, e.g., meridional heat transport, kinetic energy exchanges. For example, PCI could dictate the ocean heat-flux in the open water areas due to its influence on coastal winds. The lagged reaction would occur through prior conditioning of the ocean (due to e.g. ENSO warm/cold events), which will affect ice growth/decay with finite delay. Reference  showed that ENSO signal persists after 3 - 4 seasons in sea ice fields in the ADP regions. Moreover, this work shows that sea ice responds to atmospheric anomalies with a gap, in contrast to marginal ice response to coastal wind forcing at much shorter time scales.
Reference  synthesised the evolution and persistence of ADP anomalies in terms of the location and structure of the jets in south-eastern Pacific during ENSO warm and cold events. This study takes this understanding further by suggesting that during El Niño, the subtropical jet is robust and polar front jet weakens in the Pacific basin and vice-versa for La Niña. Therefore, more (less) heat is transported poleward in the South Pacific (South Atlantic) resulting in less (more) sea ice.
The jet stream changes reflect in the strength of Ferrel Cell, which becomes strong in the South Pacific and frail in the South Atlantic during El Niño. On a similar note, we can expect an enhanced Polar Cell in the South Pacific and weakened Polar Cell in the South Atlantic due anomalous heat flux polewards of 70˚S (Fig. 7 in  ). The changes in polar frontal jet manifest itself as the suppressor of the CPT northward expansion subsequently forming a strong polar vortex  , which will create stronger pressure gradient across Antarctica (e.g. stronger PCI, PC2500) and will have much more impact on sea ice variability.
The new results suggest that ENSO and SAM driven atmospheric modes, especially PC2500, have largest impact on Antarctic sea ice changes in the Polar Cell domain. Despite the fact that all these modes are generated distinctly and computed from different atmospheric variables, they mutually interact within the climate system of the SH. The shared impact from these regional climate modes on sea ice and the feedback generated by the tropical ENSO events are not fully comprehended and need additional inquiry. The investigation on PCI should be continued to understand atmospheric impact on seasonal sea ice variability.
Sea ice data were obtained from http://nsidc.org/data/nsidc-0051.html and NCEP2 data were downloaded from http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis2.html. The NCAR Command Language (NCL, Version 6.2.0) was used to perform EOF analysis and can be found at http://www.ncl.ucar.edu/. This is NCAOR contribution No.: 02/2016.
Figure S1. Spatial pattern of first three leading modes of (a) MSLP, (b) Z800 (800 hPa level), and (c) Z500 (500 hPa level) in the Polar Cell (poleward of 65˚S) with explained variances in percentage.
Figure S2. (a) 1979-2013 mean precipitable water content in the air column over Antarctica (source: NOAA); (b) Correlation between SAM with sea ice concentration (SIC); regions where the values are significant at 99% level are enclosed in green contours.