Introduction

Antarctic sea ice is characterized by large variations in areal extent, with substantial implications for the exchange of energy, momentum, and gas between the atmosphere and the ocean1. Antarctic sea ice varies on broad timescales from sub-seasonal to decadal and has a long-term trend as well2,3,4,5,6,7. Satellite observation (1979-now) shows Antarctic Sea ice extent (SIE) reached a record high in 2014 and then rapidly dropped to a record low in late 2016 (Fig. 1a). Prior to 2014, total Antarctic SIE had a small positive linear trend over the satellite era, with large regional variability8,9,10,11. This is in stark contrast to the declining Arctic Sea ice12,13. Moreover, the trend for SIE expansion accelerated after 2000 by almost a factor of five relative to the trend for 1979 to 199914. Several hypotheses have been proposed to explain this observed positive trend, including surface wind anomalies induced by the Southern Annual mode (SAM) or tropical SST anomalies15,16,17,18,19,20,21,22, a freshening of the Southern Ocean (SO)23,24,25,26,27,28, internal variability29,30,31, the role of mesoscale eddies32, and successful simulations of sea ice drift velocity33 and SO SST34. The accelerating rate of sea ice expansion during the post-2000 period is suggested to be related to the negative phase of the interdecadal Pacific oscillation (IPO) that can induce atmosphere teleconnection patterns and thus alter the surface wind anomalies over the SO14. A recent study showed that this increasing trend over satellite period is opposite to the longer-term variability in the recent century-long reconstruction35. A linear trend might not be the best way to describe the long-term variability in Antarctic SIE36.

Fig. 1: Variability in Antarctic Sea ice extent (SIE) derived from observations and SPEAR model simulations.
figure 1

a Monthly mean SIE anomaly since 1990 with respect to the 1990–2020 climatology in SPEAR_ECDA (green line), SPEAR_atm_restore (black line), SPEAR_atm_sst_restore (blue line), and National Snow and Ice Data Center (NSIDC) observation (red line). Magenta lines are the time averages of SIE anomalies over the time period 2012–2015 and 2017–2020 in NSIDC. b Sea ice concentration (SIC) difference between time period 2017–2020 and 2012–2015 in NSIDC. c Same as (b) but in SPEAR_atm_sst_restore. Units are 1012m2 for SIE and % for SIC.

After 2014, Antarctic SIE started to decrease and dramatically dropped to a record low in late 201637,38,39,40,41,42,43, after which the unexpected low sea ice state persisted for several years (Fig. 1a). In February 2022, the SIE reached a new record low44. Observations and model experiments by prescribing observed tropical SST anomalies suggested that the 2016 extreme sea ice low was caused by (1) southward warm air advection associated with the zonal wave number 3 teleconnection pattern in the early austral spring induced by tropical SST anomalies, and (2) the strong negative SAM event in the late austral spring38,39,40,41,42. By analyzing sparse ocean observations, previous studies also pointed out that the decadal-long SO warming associated with the SAM and IPO may contribute to the sudden Antarctic SIE decline and sustain a persistent sea ice decrease43. However, it is not clear to what extent this sudden SIE decline can be attributable to the local ocean, particularly the subsurface SO, or the relative role of the subsurface SO compared to the atmosphere contribution. There is also no answer on whether the subsurface heat is large enough to overwhelm the surface ocean variability. We are also not clear if the subsurface SO could lead to potential predictability of the persistent Antarctic SIE declines given the long memory of subsurface ocean.

In this study we use a suite of Geophysical Fluid Dynamic Laboratory (GFDL) SPEAR model simulations45 partially constrained by observations, along with sensitivity experiments (see “Methods” and Table 1 for model details), to evaluate the relative role of the subsurface SO in the recent Antarctic SIE declines. We also use initialized decadal hindcasts/forecasts46 with the SPEAR model to investigate the role of the ocean in determining sea ice cover. We find that the subsurface Southern Ocean (SO) plays a smaller role in the 2016 SIE extreme event than the atmosphere. However, the subsurface SO is critical to the persistence of the negative SIE anomalies after 2016.

Table 1 Description of models used in the current paper.

Results

Simulated Antarctic SIE variability in the GFDL SPEAR models

We show in Fig. 1a the time series of monthly Antarctic SIE anomalies in observations and various model simulations. Note that the SPEAR models we used here well simulate the sea ice climatology compared to observation (Supplementary Figs. S1, S2). The SIE climatology generally has a low bias in austral summer and a high bias in austral winter, or moderate high biases throughout all seasons, which are very common in coupled models (Supplementary Fig. S1). The simulated Antarctic sea ice concentration (SIC) is broadly like that in observation in all seasons, with high spatial correlations with observation (Supplementary Fig. S2). The SIE anomaly in each model is the deviation from each model’s monthly SIE climatology. The sea ice evolution in SPEAR_atm_sst_restore (see “Methods” and Table 1 for details) that is constrained by atmospheric and SST observations is highly correlated with NSIDC, with a correlation reaching 0.75 (p < 0.05) during the period 1990–2020. The SPEAR_atm_sst_restore captures the observed record low SIE in late 2016, and also successfully reproduces the sea ice shift in the most recent decade, with positive SIE anomalies during 2012–2015 and negative SIE anomalies during 2017–2020, albeit with a higher SIE in that later periods compared to observations. The SPEAR_atm_restore, constrained by atmospheric observations but without ERSST restoring, is still able to capture the observed SIE long-term trend, extreme event, and the recent shift. The correlation reaches 0.69 (p < 0.05) and 0.9 (p < 0.05) when comparing with the observations and SPEAR_atm_sst_restore, respectively. These high correlations suggest that the SST restoring north of 60°S does not affect the Antarctic SIE variability very much for the timescales assessed here. On the other hand, when we constrain our model with the ocean temperature and salinity observations but with a free atmosphere (SPEAR_ECDA), we find the SIE correlation with observation in this ocean data assimilation (cor = 0.6, p < 0.05) is lower than the previous two atmospheric constrained runs. This may partly arise from the insufficient ocean observations over the SO particularly in ice-covered regions that are not enough to fully constrain our coupled model. Surprisingly, the SPEAR_ECDA substantially underestimates the observed record low Antarctic SIE in late 2016. Apart from the insufficient ocean data constraints, this may also indicate that the atmosphere, such as internal atmospheric variability, may play a critical role in generating the observed record-breaking SIE extreme event. However, the SPEAR_ECDA broadly captures the SIE shift over the most recent decade, with relatively lower SIE states after 2016 compared to previous years, like that in SPEAR_atm_sst_restore and SPEAR_atm_restore.

To investigate why SPEAR_ECDA fails to reproduce the record low SIE in late 2016, we plot the SIC and sea level pressure (SLP) anomalies averaged in November and December (ND) 2016 (Supplementary Fig. S3a–d) in the SPEAR_atm_sst_restore, SPEAR_ECDA, and observations. In SPEAR_atm_sst_restore, there is a strong negative SAM event in 2016ND characterized by positive sea level pressure anomalies around Antarctica and this is in agreement with observations38,39,40,41. The negative SAM corresponds to a weakening of westerly winds over the SO, which reduces the turbulent heat loss and favors a southward Ekman transport and ocean heat transport, eventually leading to warming anomalies and SIE decline over the SO40,41. The scatter plot of SAM versus the Antarctic SIE in SPEAR_atm_sst_restore further exhibits that the more negative SAM over the SO, the lower the Antarctic SIE (Supplementary Fig. S3c). In stark contrast, we do not see a strong negative SAM event during this season in SPEAR_ECDA (no positive sea level pressure anomalies around Antarctica), for which the atmosphere is not constrained by observations. The absence of weakened westerly winds over the SO precludes the development in SPEAR_ECDA of the observed SIE extreme event, suggesting that wind-driven sea ice dynamics played a critical role in the 2016 event. Similarly, the SPEAR_ECDA underestimates the wave number 3 pattern during September and October (SO) 2016 characterized by three high and three low-pressure centers around the Southern Hemisphere extratropics because of tropical SST anomalies as mentioned by previous studies41,42 compared to the SPEAR_atm_sst_restore and observations due to its unconstrained atmosphere (Supplementary Fig. S4a–d). In SPEAR_atm_sst_restore, the southward warm air advection and dynamical ice convergence associated with the wave number 3 pattern favors a warming of SO and thus a decrease of Antarctic SIE38,39,40,41,42,43. The year with stronger southward wind anomalies corresponds to a lower SIE and vice versa. By comparing the “perfect” atmosphere reanalysis (SPEAR_atm_sst_restore; SPEAR_atm_restore) with the “perfect” ocean reanalysis (SPEAR_ECDA), we conclude that the SAM and zonal wave number 3 pattern over the SO play dominant roles in the record low Antarctic SIE in late 2016.

To explore the persistent SIE declines after 2016, we look at the spatial pattern of SIC differences between 2017–2020 and 2012–2015 (Fig. 1b, c). The sea ice decline occurs in almost all ocean basins, with a maximum over the Weddell Sea. There is an exception over the Amundsen-Bellingshausen Seas which exhibits a small SIC increase in both observations and SPEAR_atm_sst_restore, with bigger increases farther north in the model. Similar characteristics of sea ice decrease are seen in other simulations (Supplementary Fig. S5). The SIE time series in different ocean basins further confirm the total SIE declines after 2016 and the biggest SIE shift over the Weddell Sea (Supplementary Fig. S6). Again, we can see the Amundsen-Bellingshausen Seas show an opposite sea ice response compared to other basins, with increasing SIE anomalies post-2017. Next, we will explore the role of the ocean in this persistently low Antarctic SIE in recent years. Given the dominant variability over the Weddell Sea, we also show the characteristics of Weddell Sea to make a comparison with the whole SO in the next section.

Role of ocean in the persistent Antarctic SIE declines

To highlight the most recent sea ice shift, we display in Fig. 2a, b the total Antarctic SIE and Weddell SIE anomaly time series since 2011, computed with respect to the 2011–2020 climatology. Interestingly, model simulations overall reproduce the persistent sea ice lows after 2016 and match the observation much better over the 2011–2020 period than the full period shown in Fig. 1a. This may benefit from the increased and improved observations in recent years, especially the ocean part (Supplementary Fig. S7), which can better constrain our model. We then show in Fig. 2c, e the SO area averaged temperature and salinity profiles in SPEAR_atm_sst_restore. Similar characteristics are also seen in other simulations (Supplementary Fig. S8). The persistent SIE negative anomalies after 2016 correspond to an overall warm and salty upper ocean. The SIE declines are preceded by a heat buildup in the subsurface ocean (depths of roughly 100–500 m) and a salinification in the surface ocean. Both the subsurface temperature and salinity anomalies tend to destabilize the upper ocean and trigger vertical instabilities. Note that the heat buildup we have seen here takes place in the upper 500 m, which is different from the convection occurrence in the deep ocean shown previously47,48. The Weddell Sea shows similar subsurface ocean features but with a larger magnitude (Fig. 2d, f), consistent with the maximum sea ice decline there (Fig. 1b, c). We also compare these SPEAR simulations with the EN4 reanalysis49 (Supplementary Fig. S9). The subsurface warming and surface salinification signals before ND2016 are also seen in the EN4 reanalysis, although signals are nosier and intermittent presumably due to the presence of small-scale processes in the real ocean. The vertical mixing seems stronger in the EN4 than our models given the stronger penetration of surface temperature and salinity anomalies to the subsurface ocean. It is worth noting that there are still many uncertainties in the EN4 reanalysis, since the SO subsurface observations are very sparse, particularly over the Weddell Sea.

Fig. 2: Ocean changes associated with the Antarctic Sea ice extent (SIE) shift.
figure 2

a Monthly mean SIE anomaly in National Snow and Ice Data Center (NSIDC) observation (red line), SPEAR_ECDA (green line), SPEAR_atm_restore (black line), and SPEAR_atm_sst_restore (blue line). c Time-depth evolutions of the monthly Southern Ocean (SO) area averaged (60°S–70°S, 0°E–360°E) temperature anomalies in SPEAR_atm_sst_restore. e Same as (c) but for the salinity anomalies. g Scatter plot of the SO ocean stratification (defined as the density difference between 400 m and 0 m) averaged in time period 2014–2016 versus the Antarctic SIE difference (between 2017–2020 and 2012–2015). (b), (d), (f), and (h) are the same as (a), (c), (e), and (g) but for the Weddell Sea (55°S–74°S, 50°W–10°E). All anomalies are with respect to the 2011–2020 climatology. Units for the SIE, ocean temperature, salinity, and density are 1012m2, °C, PSU, and kgm−3, respectively.

Heat budget analysis reveals that the SO subsurface (100–500 m) heat reservoir in Fig. 2c is largely due to horizontal and vertical temperature advection terms (Supplementary Fig. S10). As mentioned by previous studies1,24,43, during the 2000–2014/2016 period, the SAM trend was positive, the IPO was in its negative phase and the Atlantic multidecadal oscillation (AMO) was in its positive phase (Supplementary Fig. S10a). All these large-scale climate variabilities correspond to enhanced westerlies and negative wind stress curl anomalies over the SO (Supplementary Fig. S10b). On decadal timescales, these anomalous wind stress and wind stress curls enhance SO upwelling that can bring subsurface warm water to the surface43. The subpolar gyre also spins up due to negative wind stress curl anomalies, which can advect warm water from relatively low latitude to the high latitude. The wind-driven vertical and horizontal temperature advections mutually contribute to the heat buildup in the subsurface ocean (Supplementary Fig. S10c).

Given the preconditioning of persistent sea ice declines by the subsurface heat buildup and surface salinification, we speculate that the ocean stratification (density difference between 400 and 0 m) before 2016ND may control the magnitude of this sea ice shift. We therefore show in Fig. 2g, h the scatter plot of ocean stratification averaged from 2014 to 2016 versus the SIE difference between the negative states (2017–2020) and the positive states (2012–2015) based on 30 ensemble members in SPEAR_atm_sst_restore. We can see the ensemble members with larger sea ice differences correspond to weaker ocean stratification anomalies during 2014–2016 in both the SO and Weddell Sea (Fig. 2g, h). Physically, the larger surface salt or subsurface heat perturbations and therefore weaker ocean stratification are more likely to destabilize the upper ocean and cause large sea ice anomalies and vice versa.

To further confirm the important role of ocean conditions in the Antarctic SIE shift, we conduct a sensitivity experiment (SPEAR_atm_restore_2015_newIC) based on our SPEAR_atm_restore, in which the ocean condition in Jan 1 2015 is replaced by another year (we choose year 2007 when the subsurface water is cool, see “Methods”). In the original SPEAR_atm_restore, the SIE decline persists for several years after 2016 in both the SO and Weddell Sea. The averaged SIE anomalies over the 2017–2020 period are negative (Fig. 3a, b). In contrast, the positive SIE anomalies frequently occur after 2016 when the subsurface ocean condition was changed to cold states in 2015. The averaged SIE anomalies over the 2017–2020 period become weakly positive. These sharp differences suggest that the subsurface warming before 2016 plays a crucial role in the sustained SIE low after 2016. Interestingly, the record low SIE in 2016ND is reproduced very well even when we initialize with the cold conditions, which accounts for ~90% of the low SIE in SPEAR_atm_restore. This enhances our argument in the last section that the subsurface SO prior to 2016 contributes very little to the sea ice extreme event in late 2016, and that this event was essentially driven by the atmospheric forcing.

Fig. 3: Role of ocean initial condition in the Sea ice extent (SIE).
figure 3

a The total Antarctic SIE anomaly for the time period 2011–2020 (shading) and the time-averaged (2017–2020) SIE anomalies (horizontal black line) in SPEAR_atm_restore. The brown lines denote the SIE time series and time-averaged (2017–2020) SIE anomalies in SPEAR_atm_restore_2015_newIC. b Same as (a) but for the Weddell Sea (50°W–10°E). Unit is 1012m2. The time-averaged SIE anomalies between two runs (horizontal black and brown lines) in (a) and (b) are signinficantly different at a 90% confidence level based on the Student’s t-test when we concatenate 5 ensemble members.

Figure 2g, h and the sensitivity experiment (SPEAR_atm_restore_2015_newIC) also indicate a potential predictability of the recent SIE shift since they relate the current ocean state to the future sea ice evolutions. This motivates us to investigate our decadal hindcasts/forecasts that are initialized from SPEAR_atm_sst_restore on January 1 every year (see “Methods” for more model details). Figure 4a, b displays the decadal prediction initialized in 2014, 2015, and 2016 for the total Antarctic SIE and Weddell SIE, respectively, along with the ice time series in observation and SPEAR_atm_sst_restore. All hindcasts produce small SIE values after 2016. The hindcasts initialized in 2014 and 2015 capture a large part of the recent SIE shift, moving the high sea ice states to average level. We also see that none of the hindcasts produced the record-breaking SIE decline in late 2016 (Fig. 4a). The hindcast initialized on 1 January 2016 only produces a very small SIE low in late 2016, which accounts for ~10% of the record low SIE. The hindcasts initialized in 2015 and 2014 have even less skill in predicting the 2016 extreme event. These phenomena suggest that the initial ocean condition plays a small role in the ND2016 sea ice extreme event, consistent with our arguments above and previous studies50. This is not surprising, since the internal atmosphere variability such as the SAM is hard to predict on multiyear timescales. On the other hand, the successful prediction of part of the observed sea ice difference (2017–2020 minus 2012–2015) largely arises from the correct initialization of subsurface heat reservoir (Fig. 4c, d). The initialization captures the subsurface heat evolutions very well primarily due to the long memory of subsurface ocean. The subsurface heat reduces the ocean stratification, leads to a warm and salty surface water (Supplementary Figs. S11 and 12) and eventually decreases the SIE in the surface ocean. It is worth noting that the subsurface heat can drop the high sea ice states to average level but is not enough to trigger tremendous sea ice declines after 2016. The strong interannual variability of sea ice after 2016 is also missing in the decadal prediction. This means that the subsurface ocean status before 2016 only sustains the overall low states of sea ice after 2016. The tremendous decline of sea ice after 2016 is still largely impacted by the simultaneous atmosphere and ocean variabilities. As mentioned in previous studies43, the negative SAM trend and positive IPO phase after 2016 can induce surface warming anomalies through southward warm air advection, Ekman warm advection, and Ekman downwelling and in turn decrease the Antarctic SIE. The reoccurrence of Maud Rise Polynya51,52 observed in 2016 and 2017, which is very likely driven by the strong atmospheric storm events52 can also enhance the Weddell Sea SIE decline from an ocean perspective. Our SPEAR model underestimates these polynya events probably due to the coarse ocean resolution (~0.5°) in the polar region.

Fig. 4: Decadal hindcasts of Sea ice extent (SIE) and subsurface ocean temperature.
figure 4

a Decadal hindcasts of the total Antarctic SIE initialized in 2014 (black), 2015 (blue), and 2016 (green) along with the SIE time series in SPEAR_atm_sst_restore (magenta line) and National Snow and Ice Data Center (NSIDC) (red line). Lines in the decadal hindcast denote the ensemble mean result, while the shadings denote the ensemble spread represented by the one standard deviation. c Decadal hindcasts of ensemble mean Southern Ocean (60°S–70°S, 0°E–360°E) area averaged subsurface temperature evolutions initialized in 2016, 2015, and 2014. b and d are the same as (a) and (c) but for the Weddell Sea (55°S–74°S, 50°W–10°E). Units for the SIE and ocean temperature are 1012m2, °C, respectively.

Discussion

In the present study, we have investigated the relative role of subsurface ocean and atmosphere in the recently observed Antarctic SIE declines. It is confirmed that the record low SIE in late 2016 is largely the result of both meridional temperature advection associated with the zonal wave number 3 teleconnection pattern in the early austral spring induced by tropical SST anomalies and the strong negative SAM event in the late austral spring38,39,40,41, while the subsurface SO plays a much smaller role that only accounts for ~10% of the total change. This argument is supported by three independent model experiments. The first piece of evidence arises from comparing the “perfect” atmosphere simulation with the “perfect” ocean simulation. The former reproduces the record-breaking SIE low very well, while the latter does not. This suggests that the atmosphere and its driven surface ocean and ice variabilities are key for the SIE extreme event. The “perfect” ocean simulation may be questionable due to sparse ocean observations, but our confidence largely enhances as we move to the next two pieces of evidence. The second piece of evidence is from our sensitivity experiment based on SPEAR_atm_restore in which the ocean initial conditions in 2015 were replaced by a cold subsurface state. Even with this change of ocean condition, the record low SIE in late 2016 is reproduced very well that produces ~90% of the SIE low in SPEAR_atm_restore. The third piece of evidence comes from our decadal hindcasts that are initialized from SPEAR_atm_sst_restore. The hindcast initialized on 1 January 2016 only produces a very small SIE low in late 2016, which accounts for ~10% of the record low SIE in SPEAR_atm_sst_restore. The hindcasts initialized in 2015 and 2014 have even less skill in predicting the 2016 extreme event. The second and third evidence further confirm that the record low SIE in late 2016 is not dependent on prior ocean conditions.

In contrast, the subsurface ocean plays a crucial role in the continued Antarctic SIE lows over the last decade. The persistent SIE lows after 2016 are preceded by subsurface warming. This warming is largely due to an enhanced westerly wind trend and thus negative wind stress curl anomalies over the SO in 2000–2015 that leads to an enhanced upwelling and a spin up of subpolar gyres. Our decadal hindcasts initialized in 2014 and 2015 indicate that approximately half of the SIE difference (low ice states in 2017–2020 and high ice states in 2012–2015) can be attributable to this subsurface warming. The subsurface warming persists, gradually destabilizes the ocean from below and in turn moves the high SIE states to the average level. Our sensitivity experiment also points out that the sea ice difference would disappear if the subsurface ocean before 2016 was replaced by cold conditions. This subsurface heat therefore provides a potential prediction source for the recent SIE shift due to the long memory of subsurface ocean. However, we note the subsurface heat only accounts for approximately half of the SIE difference. The other half still arises from the simultaneous atmosphere, ocean, and/or their coupled variabilities. After 2016, the negative SAM trend, the positive phase of IPO, the occurrence of Maud Rise polynyas and potentially other factors can further amplify the SIE decline. Prediction of these simultaneous climate features requires shorter timescale prediction systems and were likely unpredictable on multiyear timescales. Thus, a combined decadal46, seasonal53, and even sub-seasonal54 prediction system may be needed to better characterize the further evolutions of Antarctic SIE. Observations and our SPEAR simulations also show that the persistent sea ice declines over the recent years mainly occur over the Weddell Sea. This suggests that the Weddell Sea is very critical to control the total Antarctic SIE evolutions in the future. However, the Argo observations over this ice-covered region are very sparse compared to other areas (Supplementary Fig. S7). This calls for improved and sustained measurements of the subsurface Weddell Sea using new technologies, which could yield substantial benefits in terms of understanding and predicting future Antarctic SIE evolutions.

Methods

Observations

Here we use the National Snow and Ice Data Center (NSIDC) dataset55 to calculate the SIC spatial pattern, SIE time series and trend. SIE is defined as the cumulative area of all grid cells having at least 15% SIC. The observed subsurface temperature and salinity profiles are calculated from the Met Office Hadley Centre’s EN449 version 2.4 reanalysis.

Model

We use one model from a new set of coupled ocean-atmosphere models developed at the Geophysical Fluid Dynamics Laboratory (GFDL) (Table 1). The set is collectively called SPEAR (Seamless system for Prediction and Earth System Research)45. In this study, we use the low ocean resolution version from this set of models, called SPEAR_LO. The ocean and sea ice components are based on the MOM6 code56 and have a horizontal resolution of ~1° in the subtropics, with refinements to 0.5° at high latitudes and with meridional refinements to 1/3° in the tropics. The ocean model has 75 layers with layer thickness as fine as 2 m near the surface. The physical configurations in the atmospheric and land components are identical to the newly developed GFDL AM4-LM4 model57,58, except that we use the retuned near-infrared isotropic reflectance parameter for cold snow over glacial surfaces.

The SPEAR_LO model is then partially constrained by atmospheric observations and/or ocean observations targeting for initializing decadal or seasonal prediction system46. We have three different reanalysis datasets based on this SPEAR_LO model. The first dataset is SPEAR_atm_sst_restore, in which the 6-h atmospheric winds and temperature were restored toward atmospheric data from the 55-year Japanese Reanalysis (JRA-55)59 and the SST within 60°S–60°N was restored toward the time-varying observations from NOAA Extended Reconstructed Sea Surface Temperature version 5 (ERSSTv5) data60. The strength of the restoring is tapered linearly from 1.0 at 55°N(S) to 0.0 at 60°N(S) and is equivalent to a 10-day restoring timescale for a 50-m mixed layer depth, corresponding to an effective restoring rate of about 240 W/(m2K). This reanalysis has 30 ensemble members. The second dataset is SPEAR_atm_restore, which is identical to the SPEAR_atm_sst_restore but without SST restoring. We find that the SST restoring does not have too many effects on the Antarctic SIE time series for the timescales assessed here. Thus, we deem these two restoring analyses as the “perfect” atmosphere forcing runs. The third reanalysis we used here is called SPEAR_ECDA, which includes the ocean data assimilation (ODA) system in MOM6 as its backbone, incorporates the ocean tendency adjustment (OTA) bias correction scheme and has a free atmosphere61. MOM6 ODA uses the two‐step Ensemble Adjustment Kalman Filter (EAKF) that is an ensemble‐based sequential assimilation scheme that facilitates parallel implementation of online filter computation. We consider this SPEAR_ECDA as the “perfect” ocean reanalysis. Note that the sea ice component in these three reanalyses is not constrained by observations and primarily responds to atmospheric and ocean variations.

We conduct a sensitivity experiment to examine the role of ocean initial condition in the Antarctic SIE shift in recent years. This run is identical to the SPEAR_atm_restore reanalysis except that the ocean conditions on 1 January 2015 is replaced by that in another year. Here we choose the ocean conditions on 1 January 2007 when the subsurface ocean is relatively cool, which is in stark contrast to the warm subsurface ocean in 2015. This new run is named as SPEAR_atm_restore_2015_newIC.

We also conducted a decadal retrospective forecast using the SPEAR_LO model initialized from SPEAR_atm_sst_restore46. A 20-member ensemble of retrospective forecasts was initialized on 1 January every year from 1961 to 2021 and integrated for 10 years with the temporally varying anthropogenic and natural forcing as in the historical simulation. To effectively remove the climate drift assuming that the climate drift is systematic as a function of forecast lead time, the forecast anomalies for each variable were obtained by subtracting out the lead-time-dependent climatology from forecasts.