skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on August 2, 2025

Title: The influence of realistic 3D mantle viscosity on Antarctica’s contribution to future global sea levels
The response of the Antarctic Ice Sheet (AIS) to climate change is the largest uncertainty in projecting future sea level. The impact of three-dimensional (3D) Earth structure on the AIS and future global sea levels is assessed here by coupling a global glacial isostatic adjustment model incorporating 3D Earth structure to a dynamic ice-sheet model. We show that including 3D viscous effects produces rapid uplift in marine sectors and reduces projected ice loss for low greenhouse gas emission scenarios, lowering Antarctica’s contribution to global sea level in the coming centuries by up to ~40%. Under high-emission scenarios, ice retreat outpaces uplift, and sea-level rise is amplified by water expulsion from Antarctic marine areas.  more » « less
Award ID(s):
1745074 1744901 1744852 2035080
PAR ID:
10529771
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
American Association for the Advancement of Science
Date Published:
Journal Name:
Science Advances
Volume:
10
Issue:
31
ISSN:
2375-2548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Over the last 3.3 million years, the Antarctic Ice Sheet (AIS) has undergone phases of ice sheet growth and decay, impacting sea level and climate globally. Presently, the largely marine‐terminating AIS loses mass primarily by iceberg calving and basal melt of ice shelves. Quantifying past rates and timing of AIS melt is vital to understanding future cryosphere and sea level changes. One proxy for past ice sheet instabilities is iceberg rafted debris (IRD) fluxes. However, traditional methods of IRD quantification are labor‐intensive. Here, we present a new method of identifying IRD grains in sediment core X‐ray images using a convolutional neural network machine learning algorithm. We present a 3.3‐million‐year record of AIS IRD melt events using sediment cores from International Ocean Discovery Program Sites U1536, U1537, and U1538 in the Southern Ocean's “Iceberg Alley.” We identify two increases in the IRD fluxes throughout this period, at ∼1.8 and 0.43 Ma. We propose that after 1.8 Ma, the AIS expanded and transitioned from a primarily terrestrial‐terminating to a primarily marine‐terminating ice sheet. Therefore, after 1.8 Ma, glacial terminations and AIS iceberg discharge are associated with variations in global ice volume, presumably through the mechanism of sea level and, therefore, grounding line change. The second AIS regime change occurs during the Mid‐Brunhes Event (∼0.43 Ma). After this time, there are heightened and continuous IRD fluxes at each glacial termination, indicating increased AIS size and instability after this time. 
    more » « less
  2. Observational estimates of Antarctic ice loss have accelerated in recent decades, and worst-case scenarios of modeling studies have suggested potentially catastrophic sea level rise (~2 meters) by the end of the century. However, modeled contributions to global mean sea level from the Antarctic ice-sheet (AIS) in the 21st century are highly uncertain, in part because ice-sheet model parameters are poorly constrained. Individual ice-sheet model runs are also deterministic and not computationally efficient enough to generate the continuous probability distributions required for incorporation into a holistic framework of probabilistic sea-level projections. To address these shortfalls, we statistically emulate an ice-sheet model using Gaussian Process (GP) regression. GP modeling is a non-parametric machine-learning technique which maps inputs (e.g. forcing or model parameters) to target outputs (e.g. sea-level contributions from the Antarctic ice-sheet) and has the inherent and important advantage that emulator uncertainty is explicitly quantified. We construct emulators for the last interglacial period and an RCP8.5 scenario, and separately for the western, eastern, and total AIS. Separate emulation of western and eastern AIS is important because their evolutions and physical responses to climate forcing are distinct. The emulators are trained on 196 ensemble members for each scenario, composed by varying the parameters of maximum rate of ice-cliff wastage and the coefficient of hydrofracturing. We condition the emulators on last interglacial proxy sea-level records and modern GRACE measurements and exclude poor-fitting ensemble members. The resulting emulators are sampled to produce probability distributions that fill intermediate gaps between discrete ice-sheet model outcomes. We invert emulated high and low probability sea-level contributions in 2100 to explore 21st century evolution pathways; results highlight the deep uncertainty of ice-sheet model physics and the importance of using observations to narrow the range of parameters. Our approach is designed to be flexible such that other ice-sheet models or parameter spaces may be substituted and explored with the emulator. 
    more » « less
  3. Abstract Seismic tomography models indicate highly variable Earth structure beneath Antarctica with anomalously low shallow mantle viscosities below West Antarctica. An improved projection of the contribution of the Antarctic Ice Sheet to sea‐level change requires consideration of this complexity to precisely account for water expelled into the ocean from uplifting marine sectors. Here we build a high‐resolution 3‐D viscoelastic structure model based on recent inferences of seismic velocity heterogeneity below the continent. The model serves as input to a global‐scale sea‐level model that we use to investigate the influence of solid Earth deformation in Antarctica on future global mean sea‐level (GMSL) rise. Our calculations are based on a suite of ice mass projections generated with a range of climate forcings and suggest that water expulsion from the rebounding marine basins contributes 4%–16% and 7%–14% to the projected GMSL change at 2100 and 2500, respectively. 
    more » « less
  4. SUMMARY Mass loss from polar ice sheets is becoming the dominant contributor to current sea level changes, as well as one of the largest sources of uncertainty in sea level projections. The spatial pattern of sea level change is sensitive to the geometry of ice sheet mass changes, and local sea level changes can deviate from the global mean sea level change due to gravitational, Earth rotational and deformational (GRD) effects. The pattern of GRD sea level change associated with the melting of an ice sheet is often considered to remain relatively constant in time outside the vicinity of the ice sheet. For example, in the sea level projections from the most recent IPCC sixth assessment report (AR6), the geometry of ice sheet mass loss was treated as constant during the 21st century. However, ice sheet simulations predict that the geometry of ice mass changes across a given ice sheet and the relative mass loss from each ice sheet will vary during the coming century, producing patters of global sea level changes that are spatiotemporally variable. We adopt a sea level model that includes GRD effects and shoreline migration to calculate time-varying sea level patterns associated with projections of the Greenland and Antarctic Ice Sheets during the coming century. We find that in some cases, sea level changes can be substantially amplified above the global mean early in the century, with this amplification diminishing by 2100. We explain these differences by calculating the contributions of Earth rotation as well as gravitational and deformational effects to the projected sea level changes separately. We find in one case, for example, that ice gain on the Antarctic Peninsula can cause an amplification of up to 2.9 times the global mean sea level equivalent along South American coastlines due to positive interference of GRD effects. To explore the uncertainty introduced by differences in predicted ice mass geometry, we predict the sea level changes following end-member mass loss scenarios for various regions of the Antarctic Ice Sheet from the ISMIP6 model ensemblely, and find that sea level amplification above the global mean sea level equivalent differ by up to 1.9 times between different ice mass projections along global coastlines outside of Greenland and Antarctica. This work suggests that assessments of future sea level hazard should consider not only the integrated mass changes of ice sheets, but also temporal variations in the geometry of the ice mass changes across the ice sheets. As well, this study highlights the importance of constraining the relative timing of ice mass changes between the Greenland and Antarctic Ice Sheets. 
    more » « less
  5. Abstract Recent studies have revealed the presence of a complex freshwater system underlying the Aurora Subglacial Basin (ASB), a region of East Antarctica that contains ∼7 m of global sea level potential in ice mainly grounded below sea level. However, the impact that subglacial freshwater has on driving the evolution of the dynamic outlet glaciers that drain this basin has yet to be tested in a coupled ice sheet‐subglacial hydrology numerical modeling framework. Here, we project the evolution of the primary outlet glaciers draining the ASB (Moscow University Ice Shelf, Totten, Vanderford, and Adams Glaciers) in response to an evolving subglacial hydrology system and to ocean forcing through 2100, following low and high CMIP6 emission scenarios. By 2100, ice‐hydrology feedbacks enhance the ASB's 2100 sea level contribution by ∼30% (7.50–9.80 mm) in high emission scenarios and accelerate the retreat of Totten Glacier's main ice stream by 25 years. Ice‐hydrology feedbacks are particularly influential in the retreat of the Vanderford and Adams Glaciers, driving an additional 10 km of retreat in fully coupled simulations relative to uncoupled simulations. Hydrology‐driven ice shelf melt enhancements are the primary cause of domain‐wide mass loss in low emission scenarios, but are secondary to ice sheet frictional feedbacks under high emission scenarios. The results presented here demonstrate that ice‐subglacial hydrology interactions can significantly accelerate retreat of dynamic Antarctic glaciers and that future Antarctic sea level assessments that do not take these interactions into account might be severely underestimating Antarctic Ice Sheet mass loss. 
    more » « less