Response of Pacific-sector Antarctic ice shelves to the El Niño/Southern Oscillation

Amundsen Sea circulation controls bottom upwelling and Antarctic Pine Island and Thwaites ice shelf melting

The Pine Island and Thwaites Ice Shelves (PIIS/TIS) in the Amundsen Sea are melting rapidly and impacting global sea levels. The thermocline depth (TD) variability, the interface between cold Winter Water and warm modified Circumpolar Deep Water (mCDW), at the PIIS/TIS front strongly correlates with basal melt rates, but the drivers of its interannual variability remain uncertain. Here, using an ocean model, we propose that the strength of the eastern Amundsen Sea on-shelf circulation primarily controls TD variability and consequent PIIS/TIS melt rates. The TD variability occurs because the on-shelf circulation meanders following the submarine glacial trough, creating vertical velocity through bottom Ekman dynamics. We suggest that a strong or weak ocean circulation, possibly linked to remote winds in the Bellingshausen Sea, generates corresponding changes in bottom Ekman convergence, which modulates mCDW upwelling and TD variability. We show that interannual variability of off-shelf zonal winds has a minor effect on ocean heat intrusion into PIIS/TIS cavities, contrary to the widely accepted concept.

Progressive unanchoring of Antarctic ice shelves since 1973

Mass loss of the Antarctic Ice Sheet has been driven primarily by the thinning of the floating ice shelves that fringe the ice sheet1, reducing their buttressing potential and causing land ice to accelerate into the ocean2. Observations of ice-shelf thickness change by satellite altimetry stretch back only to 1992 (refs. 1,3-5) and previous information about thinning remains unquantified. However, extending the record of ice-shelf thickness change is possible by proxy, by measuring the change in area of the surface expression of pinning points-local bathymetric highs on which ice shelves are anchored6. Here we measure pinning-point change over three epochs spanning the periods 1973-1989, 1989-2000 and 2000-2022, and thus by proxy infer changes to ice-shelf thickness back to 1973-1989. We show that only small localized pockets of ice shelves were thinning between 1973 and 1989, located primarily in the Amundsen Sea Embayment and the Wilkes Land coastline. Ice-shelf thinning spreads rapidly into the 1990s and 2000s and is best characterized by the proportion of pinning points reducing in extent. Only 15% of pinning points reduced from 1973 to 1989, before increasing to 25% from 1989 to 2000 and 37% from 2000 to 2022. A continuation of this trend would further reduce the buttressing potential of ice shelves, enhancing ice discharge and accelerating the contribution of Antarctica to sea-level rise.

Annual mass budget of Antarctic ice shelves from 1997 to 2021

Antarctic ice shelves moderate the contribution of the Antarctic Ice Sheet to global sea level rise; however, ice shelf health remains poorly constrained. Here, we present the annual mass budget of all Antarctic ice shelves from 1997 to 2021. Out of 162 ice shelves, 71 lost mass, 29 gained mass, and 62 did not change mass significantly. Of the shelves that lost mass, 68 had statistically significant negative mass trends, 48 lost more than 30% of their initial mass, and basal melting was the dominant contributor to that mass loss at a majority (68%). At many ice shelves, mass losses due to basal melting or iceberg calving were significantly positively correlated with grounding line discharge anomalies; however, the strength and form of this relationship varied substantially between ice shelves. Our results illustrate the utility of partitioning high-resolution ice shelf mass balance observations into its components to quantify the contributors to ice shelf mass change and the response of grounded ice.

Sea level rise from West Antarctic mass loss significantly modified by large snowfall anomalies

Mass loss from the West Antarctic Ice Sheet is dominated by glaciers draining into the Amundsen Sea Embayment (ASE), yet the impact of anomalous precipitation on the mass balance of the ASE is poorly known. Here we present a 25-year (1996-2021) record of ASE input-output mass balance and evaluate how two periods of anomalous precipitation affected its sea level contribution. Since 1996, the ASE has lost 3331 ± 424 Gt ice, contributing 9.2 ± 1.2 mm to global sea level. Overall, surface mass balance anomalies contributed little (7.7%) to total mass loss; however, two anomalous precipitation events had larger, albeit short-lived, impacts on rates of mass change. During 2009-2013, persistently low snowfall led to an additional 51 ± 4 Gt yr^-1 mass loss in those years (contributing positively to the total loss of 195 ± 4 Gt yr^-1). Contrastingly, extreme precipitation in the winters of 2019 and 2020 decreased mass loss by 60 ± 16 Gt yr^-1 during those years (contributing negatively to the total loss of 107 ± 15 Gt yr^-1). These results emphasise the important impact of extreme snowfall variability on the short-term sea level contribution from West Antarctica.

Inter-decadal climate variability induces differential ice response along Pacific-facing West Antarctica

West Antarctica has experienced dramatic ice losses contributing to global sea-level rise in recent decades, particularly from Pine Island and Thwaites glaciers. Although these ice losses manifest an ongoing Marine Ice Sheet Instability, projections of their future rate are confounded by limited observations along West Antarctica's coastal perimeter with respect to how the pace of retreat can be modulated by variations in climate forcing. Here, we derive a comprehensive, 12-year record of glacier retreat around West Antarctica's Pacific-facing margin and compare this dataset to contemporaneous estimates of ice flow, mass loss, the state of the Southern Ocean and the atmosphere. Between 2003 and 2015, rates of glacier retreat and acceleration were extensive along the Bellingshausen Sea coastline, but slowed along the Amundsen Sea. We attribute this to an interdecadal suppression of westerly winds in the Amundsen Sea, which reduced warm water inflow to the Amundsen Sea Embayment. Our results provide direct observations that the pace, magnitude and extent of ice destabilization around West Antarctica vary by location, with the Amundsen Sea response most sensitive to interdecadal atmosphere-ocean variability. Thus, model projections accounting for regionally resolved ice-ocean-atmosphere interactions will be important for predicting accurately the short-term evolution of the Antarctic Ice Sheet.

Sea ice variability and trends in the Indian Ocean sector of Antarctica: Interaction with ENSO and SAM

Antarctic sea ice variability is primarily associated with ocean-atmospheric forcing driven by anomalous conditions over the tropical regions of the Pacific and Indian Oceans. The ice-ocean-atmosphere dynamics in the Indian Ocean Sector (IOS) of Antarctica have been studied using monthly satellite and reanalysis observations over four decades (1979-2019). In this study, we revealed that the annual sea ice extent (SIE) in the IOS increases at a rate of 0.7 ± 0.9% decade^-1, with a maximum increase in austral summer (5.9 ± 3.7% decade^-1). The wavelet approach was used to determine the variability in IOS sea ice caused by the El Niño/Southern Oscillation (ENSO) and southern annular mode (SAM). The SIE has a significant association with both indices during the summer and autumn. In comparison to ENSO, the sea ice variability associated with SAM is typically seasonal in nature and lacks distinct patterns. The wavelet coherence analysis revealed a relatively weak relationship between ENSO and SAM but a highly significant coherence between climatic indices and SIE. We observed that sea ice in the IOS is influenced significantly by climatic oscillations during their negative SAM/El Niño or positive SAM/La Niña phases. Furthermore, the study demonstrated a substantial impact of climatic disturbances in determining the sea ice variability in the IOS.

Antarctic Peninsula warming triggers enhanced basal melt rates throughout West Antarctica

The observed acceleration of ice shelf basal melt rates throughout West Antarctica could destabilize continental ice sheets and markedly increase global sea level. Explanations for decadal-scale melt intensification have focused on processes local to shelf seas surrounding the ice shelves. A suite of process-based model experiments, guided by CMIP6 forcing scenarios, show that freshwater forcing from the Antarctic Peninsula, propagated between marginal seas by a coastal boundary current, causes enhanced melting throughout West Antarctica. The freshwater anomaly stratifies the ocean in front of the ice shelves and modifies vertical and lateral heat fluxes, enhancing heat transport into ice shelf cavities and increasing basal melt. Increased glacial runoff at the Antarctic Peninsula, one of the first signatures of a warming climate in Antarctica, emerges as a key trigger for increased ice shelf melt rates in the Amundsen and Bellingshausen Seas.

High spatial and temporal variability in Antarctic ice discharge linked to ice shelf buttressing and bed geometry

Antarctica's contribution to global mean sea level rise has been driven by an increase in ice discharge into the oceans. The rate of change and the mechanisms that drive variability in ice discharge are therefore important to consider in the context of projected future warming. Here, we report observations of both decadal trends and inter-annual variability in ice discharge across the Antarctic Ice Sheet at a variety of spatial scales that range from large drainage basins to individual outlet glacier catchments. Overall, we find a 37 ± 11 Gt year^-1 increase in discharge between 1999 and 2010, but a much smaller increase of 4 ± 8 Gt year^-1 between 2010 and 2018. Furthermore, comparisons reveal that neighbouring outlet glaciers can behave synchronously, but others show opposing trends, despite their close proximity. We link this spatial and temporal variability to changes in ice shelf buttressing and the modulating effect of local glacier geometry.

Decadal-scale onset and termination of Antarctic ice-mass loss during the last deglaciation

Emerging ice-sheet modeling suggests once initiated, retreat of the Antarctic Ice Sheet (AIS) can continue for centuries. Unfortunately, the short observational record cannot resolve the tipping points, rate of change, and timescale of responses. Iceberg-rafted debris data from Iceberg Alley identify eight retreat phases after the Last Glacial Maximum that each destabilized the AIS within a decade, contributing to global sea-level rise for centuries to a millennium, which subsequently re-stabilized equally rapidly. This dynamic response of the AIS is supported by (i) a West Antarctic blue ice record of ice-elevation drawdown >600 m during three such retreat events related to globally recognized deglacial meltwater pulses, (ii) step-wise retreat up to 400 km across the Ross Sea shelf, (iii) independent ice sheet modeling, and (iv) tipping point analysis. Our findings are consistent with a growing body of evidence suggesting the recent acceleration of AIS mass loss may mark the beginning of a prolonged period of ice sheet retreat and substantial global sea level rise.

Ocean warming threatens southern right whale population recovery

Whales contribute to marine ecosystem functioning, and they may play a role in mitigating climate change and supporting the Antarctic krill (Euphausia superba) population, a keystone prey species that sustains the entire Southern Ocean (SO) ecosystem. By analyzing a five-decade (1971–2017) data series of individual southern right whales (SRWs; Eubalaena australis) photo-identified at Península Valdés, Argentina, we found a marked increase in whale mortality rates following El Niño events. By modeling how the population responds to changes in the frequency and intensity of El Niño events, we found that such events are likely to impede SRW population recovery and could even cause population decline. Such outcomes have the potential to disrupt food-web interactions in the SO, weakening that ecosystem’s contribution to the mitigation of climate change at a global scale.

Complex Principal Component Analysis of Antarctic Ice Sheet Mass Balance

Ice sheet changes of the Antarctic are the result of interactions among the ocean, atmosphere, and ice sheet. Studying the ice sheet mass variations helps us to understand the possible reasons for these changes. We used 164 months of Gravity Recovery and Climate Experiment (GRACE) satellite time-varying solutions to study the principal components (PCs) of the Antarctic ice sheet mass change and their time-frequency variation. This assessment was based on complex principal component analysis (CPCA) and the wavelet amplitude-period spectrum (WAPS) method to study the PCs and their time-frequency information. The CPCA results revealed the PCs that affect the ice sheet balance, and the wavelet analysis exposed the time-frequency variation of the quasi-periodic signal in each component. The results show that the first PC, which has a linear term and low-frequency signals with periods greater than five years, dominates the variation trend of ice sheet in the Antarctic. The ratio of its variance to the total variance shows that the first PC explains 83.73% of the mass change in the ice sheet. Similar low-frequency signals are also found in the meridional wind at 700 hPa in the South Pacific and the sea surface temperature anomaly (SSTA) in the equatorial Pacific, with the correlation between the low-frequency periodic signal of SSTA in the equatorial Pacific and the first PC of the ice sheet mass change in Antarctica found to be 0.73. The phase signals in the mass change of West Antarctica indicate the upstream propagation of mass loss information over time from the ocean–ice interface to the southward upslope, which mainly reflects ocean-driven factors such as enhanced ice–ocean interaction and the intrusion of warm saline water into the cavities under ice shelves associated with ice sheets which sit on retrograde slopes. Meanwhile, the phase signals in the mass change of East Antarctica indicate the downstream propagation of mass increase information from the South Pole toward Dronning Maud Land, which mainly reflects atmospheric factors such as precipitation accumulation.

Antarctic ice mass variations from 1979 to 2017 driven by anomalous precipitation accumulation

Antarctic ice mass balance is determined by precipitation and ice discharge, and understanding their relative contributions to contemporary Antarctic ice mass change is important to project future ice mass loss and resulting sea level rise. There has been evidence that anomalous precipitation affects Antarctic ice mass loss estimates, and thus the precipitation contribution should be understood and considered in future projections. In this study, we revisit changes in Antarctic ice mass over recent decades and examine precipitation contributions over this period. We show that accumulated (time-integrated) precipitation explains most inter-annual anomalies of Antarctic ice mass change during the GRACE period (2003–2017). From 1979 to 2017, accumulated Antarctic precipitation contributes to significant ice mass loss acceleration in the Pacific sector and deceleration in the Atlantic-Indian Sectors, forming a bi-polar spatial pattern. Principal component analysis reveals that such a bi-polar pattern is likely modulated by the Southern Annular Mode (SAM). We also find that recent ice mass loss acceleration in 2007 is related to a variation in precipitation accumulation. Overall ice discharge has accelerated at a steady rate since 1992, but has not seen a recent abrupt increase.

Damage accelerates ice shelf instability and mass loss in Amundsen Sea Embayment

Pine Island Glacier and Thwaites Glacier in the Amundsen Sea Embayment are among the fastest changing outlet glaciers in Antarctica. Yet, projecting the future of these glaciers remains a major uncertainty for sea level rise. Here we use satellite imagery to show the development of damage areas with crevasses and open fractures on Pine Island and Thwaites ice shelves. These damage areas are first signs of their structural weakening as they precondition these ice shelves for disintegration. Model results that include the damage mechanism highlight the importance of damage for ice shelf stability, grounding line retreat, and future sea level contributions from Antarctica. Moreover, they underline the need for incorporating damage processes in models to improve sea level rise projections.

Pine Island Glacier and Thwaites Glacier in the Amundsen Sea Embayment are among the fastest changing outlet glaciers in West Antarctica with large consequences for global sea level. Yet, assessing how much and how fast both glaciers will weaken if these changes continue remains a major uncertainty as many of the processes that control their ice shelf weakening and grounding line retreat are not well understood. Here, we combine multisource satellite imagery with modeling to uncover the rapid development of damage areas in the shear zones of Pine Island and Thwaites ice shelves. These damage areas consist of highly crevassed areas and open fractures and are first signs that the shear zones of both ice shelves have structurally weakened over the past decade. Idealized model results reveal moreover that the damage initiates a feedback process where initial ice shelf weakening triggers the development of damage in their shear zones, which results in further speedup, shearing, and weakening, hence promoting additional damage development. This damage feedback potentially preconditions these ice shelves for disintegration and enhances grounding line retreat. The results of this study suggest that damage feedback processes are key to future ice shelf stability, grounding line retreat, and sea level contributions from Antarctica. Moreover, they underline the need for incorporating these feedback processes, which are currently not accounted for in most ice sheet models, to improve sea level rise projections.

Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves

Ocean-driven basal melting of Antarctica's floating ice shelves accounts for about half of their mass loss in steady-state, where gains in ice shelf mass are balanced by losses. Ice shelf thickness changes driven by varying basal melt rates modulate mass loss from the grounded ice sheet and its contribution to sea level, and the changing meltwater fluxes influence climate processes in the Southern Ocean. Existing continent-wide melt rate datasets have no temporal variability, introducing uncertainties in sea level and climate projections. Here, we combine surface height data from satellite radar altimeters with satellite-derived ice velocities and a new model of firn-layer evolution to generate a high-resolution map of time-averaged (2010-2018) basal melt rates, and time series (1994-2018) of meltwater fluxes for most ice shelves. Total basal meltwater flux in 1994 (1090±150 Gt/yr) was not significantly different from the steady-state value (1100±60 Gt/yr), but increased to 1570±140 Gt/yr in 2009, followed by a decline to 1160±150 Gt/yr in 2018. For the four largest "cold-water" ice shelves we partition meltwater fluxes into deep and shallow sources to reveal distinct signatures of temporal variability, providing insights into climate forcing of basal melting and the impact of this melting on the Southern Ocean.

The Southern Ocean and its interaction with the Antarctic Ice Sheet

The Southern Ocean exerts a major influence on the mass balance of the Antarctic Ice Sheet, either indirectly, by its influence on air temperatures and winds, or directly, mostly through its effects on ice shelves. How much melting the ocean causes depends on the temperature of the water, which in turn is controlled by the combination of the thermal structure of the surrounding ocean and local ocean circulation, which in turn is determined largely by winds and bathymetry. As climate warms and atmospheric circulation changes, there will be follow-on changes in the ocean circulation and temperature. These consequences will affect the pace of mass loss of the Antarctic Ice Sheet.

Rebound of shelf water salinity in the Ross Sea

Antarctic Bottom Water (AABW) supplies the lower limb of the global overturning circulation and ventilates the abyssal ocean. In recent decades, AABW has warmed, freshened and reduced in volume. Ross Sea Bottom Water (RSBW), the second largest source of AABW, has experienced the largest freshening. Here we use 23 years of summer measurements to document temporal variability in the salinity of the Ross Sea High Salinity Shelf Water (HSSW), a precursor to RSBW. HSSW salinity decreased between 1995 and 2014, consistent with freshening observed between 1958 and 2008. However, HSSW salinity rebounded sharply after 2014, with values in 2018 similar to those observed in the mid-late 1990s. Near-synchronous interannual fluctuations in salinity observed at five locations on the continental shelf suggest that upstream preconditioning and large-scale forcing influence HSSW salinity. The rate, magnitude and duration of the recent salinity increase are unusual in the context of the (sparse) observational record.

Ross Sea Bottom Water, a major source of Antarctic Bottom Water, has experienced significant freshening in recent decades. Here the authors use 23 years of summer measurements to document temporal variability in the salinity of the Ross Sea High Salinity Shelf Water (HSSW) and found that HSSW salinity decreased between 1995 and 2014 and rebounded sharply after 2014.

Marine ice sheet instability amplifies and skews uncertainty in projections of future sea-level rise

The potential for collapse of the Antarctic ice sheet remains the largest single source of uncertainty in projections of future sea-level rise. This uncertainty comes from an imperfect understanding of ice sheet processes and the internal variability of climate forcing of ice sheets. Using a mathematical technique from statistical physics and large ensembles of state-of-the-art ice sheet simulations, we show that collapse of ice sheets widens the range of possible scenarios for future sea-level rise. We also find that the collapse of marine ice sheets makes worst-case scenarios of rapid sea-level rise more likely in future projections.

Sea-level rise may accelerate significantly if marine ice sheets become unstable. If such instability occurs, there would be considerable uncertainty in future sea-level rise projections due to imperfectly modeled ice sheet processes and unpredictable climate variability. In this study, we use mathematical and computational approaches to identify the ice sheet processes that drive uncertainty in sea-level projections. Using stochastic perturbation theory from statistical physics as a tool, we show mathematically that the marine ice sheet instability greatly amplifies and skews uncertainty in sea-level projections with worst-case scenarios of rapid sea-level rise being more likely than best-case scenarios of slower sea-level rise. We also perform large ensemble simulations with a state-of-the-art ice sheet model of Thwaites Glacier, a marine-terminating glacier in West Antarctica that is thought to be unstable. These ensemble simulations indicate that the uncertainty solely related to internal climate variability can be a large fraction of the total ice loss expected from Thwaites Glacier. We conclude that internal climate variability alone can be responsible for significant uncertainty in projections of sea-level rise and that large ensembles are a necessary tool for quantifying the upper bounds of this uncertainty.

Intrinsic processes drive variability in basal melting of the Totten Glacier Ice Shelf

Over the period 2003–2008, the Totten Ice Shelf (TIS) was shown to be rapidly thinning, likely due to basal melting. However, a recent study using a longer time series found high interannual variability present in TIS surface elevation without any apparent trend. Here we show that low-frequency intrinsic ocean variability potentially accounts for a large fraction of the variability in the basal melting of TIS. Specifically, numerical ocean model simulations show that up to 44% of the modelled variability in basal melting in the 1–5 year timescale (and up to 21% in the 5–10 year timescale) is intrinsic, with a similar response to the full climate forcing. We identify the important role of intrinsic ocean variability in setting the observed interannual variation in TIS surface thickness and velocity. Our results further demonstrate the need to account for intrinsic ocean processes in the detection and attribution of change.

Low frequency intrinsic ocean variability has an unknown impact on Antarctic ice shelves, yet can arise even in the absence of varying climate forcing. Here, the authors show that this variability significantly affects modelled basal melting under the Totten Ice Shelf, with implications for the attribution of change.

Choosing the future of Antarctica

We present two narratives on the future of Antarctica and the Southern Ocean, from the perspective of an observer looking back from 2070. In the first scenario, greenhouse gas emissions remained unchecked, the climate continued to warm, and the policy response was ineffective; this had large ramifications in Antarctica and the Southern Ocean, with worldwide impacts. In the second scenario, ambitious action was taken to limit greenhouse gas emissions and to establish policies that reduced anthropogenic pressure on the environment, slowing the rate of change in Antarctica. Choices made in the next decade will determine what trajectory is realized.