Characteristics of 2022–2023 annual cycle. We start our analysis by characterizing the development of SIE from March 2022 to February 2023 and by placing it in the broader historical context. Figure 1a displays the evolution of observed SIE from March 2022 to February 2023 and matches it with the corresponding observations for 2016–2017 and 2021–2022 (i.e., the two lowest years on record before the event under study), 2017–2018 (i.e., the year following the 2016–2017 summer minimum), as well as the climatological average (1981–2010). We find that the 2023 record results from the combination of two characteristics. First, a consistently lower-than-average SIE prevailed during the whole ice growth season, which found its origin as early as February 2022 (i.e., the previous record), suggesting a strong role for preconditioning. Second, melt rates (measured as day-to-day changes in sea ice extent; Fig. 1b) were unusually high in December 2022, reaching the absolute maximum in melting rate for all the investigated years. These exceptional melt rates further accelerated sea ice loss towards the summer minimum. It is also worth noting that the 2022–2023 sea ice started to retreat at a normal timing in September (Figs. 1a, 1b), unlike the earlier retreat during former two record lows12,44.
The SIE is a limited indicator of changes as it neglects spatial aspects. A look at the spatial expression of the anomalies discussed above can further direct the search for the mechanisms underlying these changes. On 21 February, when the SIE reached its minimum value, SIC anomalies displayed a nearly circumpolar pattern except in the western Indian Ocean (Fig. 1c), while sea ice anomalies were more regional in the former minimum events15. The negative SIC anomalies in 2023 are particularly evident in the Amundsen, Bellingshausen and Ross Seas regions. A look at the SIC melt rate anomalies in December (Fig. 1d) reveals that the intense melt rates in SIE in December are dominated by strong SIC changes in the north-eastern Ross Sea.
Based on these observations, the 2023 summer SIE minimum can be characterized as (1) having inherited low SIE conditions from the previous summer (2022) and persisted throughout the year, (2) having resulted from exceptional late-spring/early-summer melt rates, which mostly operated in the north-eastern Ross Sea where sea ice normally survives into the summer, and (3) having displayed consistent reductions in other sectors where summer sea ice is climatologically extensive. One exception is the western Indian sector, but the contribution of this sector to the total SIE is negligible (5% for the SIE climatological average)45.
Recent studies pointed out that subsurface warming played a decisive role in causing more frequent sea ice extremes13,21. A potential regime shift to persistent low-extent sea ice state was also suggested. To explore if the subsurface warming has had an effect on the 2023 minimum, we inspected the upper-ocean temperature anomalies in different regions derived from the PAROCE simulation during 2022–2023 (Supplementary Fig. 4). We find that warmer subsurface ocean conditions have prevailed in all basins of the Southern Ocean in the first 100 m of the ocean, except in the Amundsen and Bellingshausen Seas in summer. The warming was particularly prominent in the Indian and Pacific sectors. We note that the oceanic warming is not specifically high in the eastern Ross Sea, which indicates the potential impacts of other processes on the extreme sea ice melt there. From a long-term perspective (Supplementary Fig. 5), it appears that a remarkable warming has also been initiated in these sectors since 2019. A warmer ocean state can favor the basal and lateral melting or hinder the basal growth in such event, as basal melt is a dominant source for sea ice melting in the Southern Ocean46,47, but the detailed mechanism of ice-ocean interactions is beyond the scope of this study. In the following, we will focus on the impacts of local atmospheric processes on the sea ice anomalies and especially investigate the cause of the high melting rates in the eastern Ross Sea.
Local impacts of atmospheric anomalies. Now that we have described the temporal and spatial characteristics of SIE throughout 2022–2023, we can turn our attention to the candidate mechanisms that drove these changes. We examine here the atmospheric variables and SIC anomalies and their potential mutual connections during the four seasons of 2022–2023. Note that the division of the Southern Ocean here (Fig. 2m) is made as to highlight the special feature of the eastern Ross Sea.
In autumn, a slightly deepened ASL established in the eastern Ross Sea (Fig. 2a) and air temperature was anomalously warm in the eastern Pacific region (ABS + ERS; Fig. 2e). The intraseasonal air temperature changes mainly come from strong heat transport from lower latitudes through warm air advection (Fig. 2i), hindering offshore sea ice growth in the eastern Pacific region. In the coastal region of the Ross Sea, less sea ice and warmer air temperature can be seen while cold advection started to prevail and led to decreased temperature anomalies at the end of autumn.
During winter, the negative SLP anomaly moved eastward to the ABS region (Fig. 2b), causing positive temperature anomalies in the Weddell Sea, Bellingshausen Sea, eastern Indian Ocean (Fig. 2f) with warm advection and negative SIC anomalies there (Fig. 2j, 2n). In contrast, more sea ice was produced due to cold advection in the Ross Sea, as expected with the typical ASL pattern.
During spring, the SLP anomaly deepened in the Bellingshausen Sea (Fig. 2c), accompanied by strong warm advection in the Bellingshausen Sea and Eastern Antarctic and cold advection in the Ross Sea (Fig. 2k). Consequently, a SIC dipole (positive anomalies in the Ross Sea and negative ones in the Bellingshausen Sea) appeared and SIC decreased in the Eastern Antarctic (Fig. 2o). Due to anomalous northerly winds, more ice was transported to the northern Ross Sea and the ice edge was pushed further north by about 2 degrees of latitude till late November (Supplementary Fig. 2a). After that, sea ice edge anomalies retreated abruptly by about 5 degrees of latitude in December. Little ice survived in the Bellingshausen Sea in October after the consistent warm advection pattern that had prevailed.
In summer, SLP presented an annular negative anomaly over the whole Antarctic and the regional temperature anomalies were limited. Sea ice experienced stronger melting than usual in most regions, consistent with warm advections in the ABS and Eastern Antarctic (Fig. 2l). Nevertheless, a large polynya appeared along the coast of Amundsen Sea in December (see the red domain highlighted in Fig. 2p). Interestingly, air temperature anomalies and advection anomalies were only weakly positive over that domain in summer (Fig. 2h). This apparent inconsistency between the sea ice state and the atmospheric forcing pushed us to elucidate the origins of this polynya from the aspect of preconditions since the SIE reduction in this region (70–80S, 120–180W) accounts for 37% of the total reduction in February.
Impact of spring preconditions in the eastern Ross Sea. Guided by the analysis above, we now focus on the preceding atmospheric impacts on sea ice anomalies in the eastern Ross Sea. In this region, sea ice edge anomalies are exclusively positive in spring and dropped abruptly afterwards (Supplementary Fig. 2a), while there was an extreme melting rate and the development of a coastal polynya in December.
Since the sea ice anomalies happened without significant atmospheric anomalies, we investigate the impacts of preceding sea ice volume (SIV) anomalies by showing the modeled SIC and SIV changes during October-December 2022 and climatology from the PAROCE configuration in the polynya region (Fig. 3a, 3b). The performance of the model output has been assessed (see Methods, Supplementary Figs. 1–3) and the model is considered to be sufficiently realistic to reproduce sea ice anomalies in this region. Although SIC showed few anomalies and remained high before December, the SIV has already presented large negative anomalies since the beginning of sea ice growth in 2022 based on the model outputs, indicating that the generally thinner ice along the coast in spring is a potential precondition for the formation of polynya. Different from the conditions inside (onshore) the polynya region, sea ice displayed higher SIC than average outside (offshore) the region (Fig. 3c) and positive anomalies in SIV (Fig. 3d) before December. After that, both SIC and SIV started to decrease. This dipole between negative SIV anomalies close to the coast and positive ones offshore suggests an impact of northward transport from coastal region to northeastern Ross Sea. A SIT tendency budget analysis is further conducted below with the model outputs to confirm this hypothesis (note that the modeled SIT refers to the grid-cell average thickness and is equal to the SIV divided by the grid cell area).
Figure 4 shows the modeled SIT budgets inside and outside the Amundsen polynya region during October–December 2022 and climatology and each term has been quantified in Supplementary Table 1–2. The daily SIT changes have been divided into growth terms (basal growth, lateral growth, growth due to snow-ice formation, growth due to dynamic transportation) and melt terms (surface melt, lateral melt, basal melt, melt due to dynamic transportation). Inside the polynya region, growth terms in 2022 are generally similar to the climatology, except some transportation-induced growth in December due to anomalous northerly winds (Fig. 2l). However, extreme transportation-induced melt dominated the total SIT tendency from October to late November 2022 (more than twice the climatology), indicating a strong outwards SIV transport. From mid-November to end-December 2022, anomalous melting developed rapidly. This melt was first primarily driven by surface melt (three-fold of climatology), accompanied by growing negative SIC anomalies (Fig. 3a) and thus more open water fraction. From late November onward, oceanic melt (basal melt plus lateral melt) took over and accounted for the largest contribution of the total melt (as also observed for the climatology). In late December, the total melt decreased due to little SIV left (Fig. 3b). Correspondingly, outside the polynya region, transportation-induced growth, which originates from sea ice transportation, is significantly larger during October to November 2022 than climatology. Surface melt is also anomalously large from late-November to mid-December 2022 while oceanic melt is comparable to climatology. Checking the modeled snow depth anomalies over the polynya region (Supplementary Fig. 6), we find that the snow depth presented year-round negative anomalies and a marked drop from September to December, suggesting that surface melt was accompanied by a thinner late-winter/spring snow pack on sea ice.
We also analyze the net surface heat flux anomalies from ERA5 reanalysis, and further separate it into long-wave radiation, short-wave radiation, latent heat flux and sensible heat flux (Supplementary Fig. 7). The net surface heat flux shows increased absorption in the polynya region, dominated by net short-wave radiation flux anomalies. By splitting the net short-wave flux into a downward and an upward component, we can attribute the net surface heat flux anomalies to a less reflective surface, thus implying that the positive ice-albedo feedback48 controls those net short-wave anomalies.
The deepened ASL in spring has been demonstrated as the potential trigger of positive ice-albedo feedback accompanied with a large sea ice area flux outwards the region27. Here, we use the model outputs to calculate the monthly SIV transport outwards the three gates of Amundsen polynya region, and to connect them with the changing ASL index with the ice-albedo feedback initiation. Specifically, supplementary Fig. 8 displays the anomalies of the ASL depth, the ASL longitudinal position, the ERA5 meridional winds, the SIV transport, and the modeled albedo during August–December 2022. The ASL remained deeper than in the climatology throughout the entire winter, spring and summer seasons. In October, though the depth anomalies were slightly smaller, the ASL was located more eastward than usual, bringing strong southerly winds in the polynya region. Similar ASL conditions lasted until November, resulting in the largest SIV transport anomalies outwards. Consequently, large fractions of open water appeared, surface albedo decreased, and upper-ocean was heated by solar energy, triggering the positive ice-albedo feedback.
Concluding remarks. The new Antarctic sea ice minimum in summer 2023 happened following the last minimum in 2022. The 2022–2023 sea ice evolution featured year-round lower-than-average SIE since March, extremely high melting rate in December in the northeastern Ross Sea and a circumpolar melting in summer. In general, negative sea ice anomalies are consistent with advection-induced air temperature anomalies, especially in the Weddell Sea, Bellingshausen Sea and Eastern Antarctic. In the Eastern Ross Sea, strong southerly winds driven by ASL transported sea ice along the coast northward, forming a large polynya and increasing the solar absorption of upper ocean. As a result, severe surface melting happened and the ice-albedo feedback amplified the initial change.
The recurrent low summer sea ice minimums in recent years as well as the exceptionally low winter sea ice maximum in 2023 pose a direct threat to coastal areas. The potential implications encompass coastal erosion, reduced ice-shelf stability, and disruptions to ecosystems. Consecutive observations and improved model performances are imperative to understand the multiscale variabilities inherent in the evolving state of Antarctic sea ice.