Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice sheet”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sensitivities of the West Greenland Current to Greenland Ice Sheet Meltwater in a Mesoscale Ocean/Sea Ice Model

Abstract Meltwater from the Greenland Ice Sheet can alter the continental shelf/slope circulation and cross-shelf freshwater fluxes and limit deep convection in adjacent basins through surface freshening. We explore the impacts on the West Greenland Current and eastern Labrador Sea with different vertical distributions of the meltwater forcing. In this study, we present the results from global coupled ocean/sea ice simulations, forced with atmospheric reanalysis, that are mesoscale eddy-active (∼2–3-km horizontal spacing) and eddy-permitting (∼6–7-km horizontal spacing) in the study region. We compare the West Greenland Current in mesoscale eddy-active and eddy-permitting without meltwater to highlight the role of small-scale features. The mesoscale eddy-active configuration is then used to assess the change in the eastern Labrador Sea when meltwater is added to the surface or vertically distributed to account for mixing within fjords. In both simulations with meltwater, the West Greenland and West Greenland Coastal Currents are faster than in the simulation with no meltwater; their mean surface speeds are the highest in the vertical distribution case. In the latter case, there is enhanced baroclinic conversion at the shelf break compared to the simulation with no meltwater. When meltwater is vertically distributed, there is an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy. In addition, in the eastern Labrador Sea, the salinity is lower and the meltwater volume is greater when meltwater is vertically distributed. Therefore, the West Greenland Current is sensitive to how meltwater is added to the ocean with implications for the freshening of the Labrador Sea. Significance Statement Our goal is to understand how the flux of freshwater across the West Greenland continental slope into the Labrador Sea is modified by meltwater from the Greenland Ice Sheet. We compare the simulations of the ocean that capture key dynamics along the West Greenland continental slope that have no meltwater, meltwater added to the ocean surface, and meltwater distributed vertically to represent the mixing within fjords. When meltwater is added, the currents along the continental slope are faster, with the greatest increase when meltwater is vertically distributed. In that case, there is enhanced freshening of the Labrador Sea because modified density gradients generate more eddies. Proper representation of the vertical structure of meltwater is important for projecting the impact of freshwater on the subpolar North Atlantic.

58 GEOSCIENCES↗
DOE OSTI2025

ISMIP6-based Antarctic projections to 2100: simulations with the BISICLES ice sheet model

The contribution of the Antarctic Ice Sheet is one of the most uncertain components of sea level rise to 2100. Ice sheet models are the primary tool for projecting future sea level contribution from continental ice sheets. The Ice Sheet Model Intercomparison for the Coupled Model Intercomparison Phase 6 (ISMIP6) provided projections of the ice sheet contribution to sea level over the 21st century, quantifying uncertainty due to ice sheet model, climate model, emission scenario, and uncertain parameters. We present simulations following the ISMIP6 framework with the BISICLES ice sheet model and new experiments extending the ISMIP6 protocol to more comprehensively sample uncertainties in future climate, ice shelf sensitivity to ocean melting, and their interactions. These results contributed to the land ice projections of Edwards et al. (2021), which formed the basis of sea level projections for the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (AR6). Our experiments show the important interplay between surface mass balance processes and ocean-driven melt in determining Antarctic sea level contribution. Under higher-warming scenarios, high accumulation offsets more ocean-driven mass loss when sensitivity to ocean-driven melt is low. Conversely, we show that when sensitivity to ocean warming is high, ocean melting drives increased mass loss despite high accumulation. Overall, we simulate a sea level contribution range across our experiments from 2 to 178 mm. Finally, we show that collapse of ice shelves due to surface warming increases sea level contribution by 25 mm relative to the no-collapse experiments, for both moderate and high sensitivity of ice shelf melting to ocean forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Plant, insect, and fungi fossils under the center of Greenland’s ice sheet are evidence of ice-free times

The persistence and size of the Greenland Ice Sheet (GrIS) through the Pleistocene is uncertain. This is important because reconstructing changes in the GrIS determines its contribution to sea level rise during prior warm climate periods and informs future projections. To understand better the history of Greenland’s ice, we analyzed glacial till collected in 1993 from below 3 km of ice at Summit, Greenland. The till contains plant fragments, wood, insect parts, fungi, and cosmogenic nuclides showing that the bed of the GrIS at Summit is a long-lived, stable land surface preserving a record of deposition, exposure, and interglacial ecosystems. Knowing that central Greenland was tundra-covered during the Pleistocene informs the understanding of Arctic biosphere response to deglaciation.

58 GEOSCIENCES↗
DOE OSTI2023

Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

Abstract. The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗
DOE OSTI2024

Experimental design for the Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2)

The Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2) is a natural progression of previous and ongoing model intercomparison exercises that have focused on the simulation of ice-sheet and ocean processes in Antarctica. The previous exercises motivate the move towards realistic configurations, as well as more diverse model parameters and resolutions. The main objective of MISOMIP2 is to investigate the performance of existing ocean and coupled ice-sheet–ocean models in a range of Antarctic environments through comparisons to observational data. We will assess the status of ice-sheet–ocean modelling as a community and identify common characteristics of models that are best able to capture observed features. As models are highly tuned based on present-day data, we will also compare their sensitivity to prescribed abrupt atmospheric perturbations leading to either very warm or slightly warmer ocean conditions compared to the present day. The approach of MISOMIP2 is to welcome contributions of models as they are, including global and regional configurations, but we request standardized variables and common grids for the outputs. We target the analysis at two specific regions, the Amundsen Sea and the Weddell Sea, since they describe two different ocean environments and have been relatively well observed compared to other areas of Antarctica. An observational “MIPkit” synthesizing existing ocean and ice-sheet observations for a common period is provided to evaluate ocean and ice-sheet models in these two regions.

58 GEOSCIENCES↗
DOE OSTI2025

East Antarctic Ice Sheet variability in the central Transantarctic Mountains since the mid Miocene

The response of the East Antarctic Ice Sheet to warmer-than-present climate conditions has direct implications for projections of future sea level, ocean circulation, and global radiative forcing. Nonetheless, it remains uncertain whether the ice sheet is likely to undergo net loss due to amplified melting coupled with dynamic instabilities or whether such losses will be balanced, or even offset, by enhanced accumulation under a higher-precipitation regime. The glacial depositional record from the central Transantarctic Mountains (TAM) provides a robust geologic means to reconstruct the past behaviour of the East Antarctic Ice Sheet, including during periods thought to have been warmer than today, such as the mid-Pliocene Warm Period (~3.3–3.0 Ma). This study describes a new surface-exposure-dated moraine record from Otway Massif in the central TAM spanning the last ~9 Myr and synthesises these data in the context of previously published moraine chronologies constrained with cosmogenic nuclides. The resulting record, although fragmentary, represents the majority of direct and unambiguous terrestrial evidence for the existence and size of the East Antarctic Ice Sheet during the last 14 Myr, and it thus provides new insight into the long-term relationship between the ice sheet and global climate. At face value, the existing TAM moraine record does not exhibit a clear signature of the mid-Pliocene Warm Period, thus precluding a definitive verdict on the East Antarctic Ice Sheet's response to this event. In contrast, an apparent hiatus in moraine deposition both at Otway Massif and the neighbouring Roberts Massif suggests that the ice sheet surface in the central TAM was potentially lower than present during the late Miocene and earliest Pliocene.

58 GEOSCIENCES↗
DOE OSTI2024

The Nonlinear and Distinct Responses of Ocean Heat Content and Anthropogenic Carbon to Ice Sheet Freshwater Discharge in a Warming Climate

Anthropogenic climate change will drive extensive mass loss across both the Antarctic (AIS) and Greenland Ice Sheets (GrIS), with the potential for global climate system feedbacks, especially in polar regions. Historically, the high-latitude North Atlantic and Southern Ocean have been critical regions for anthropogenic heat and carbon uptake, but our understanding of how this uptake will be altered by future freshwater discharge is incomplete. We assess each ice sheet's impact on global ocean anthropogenic heat and carbon storage for a high-emission scenario over the 21 st -century using a coupled Earth system model. We explore the impact of contemporaneous mass loss from both ice sheets on anthropogenic heat and carbon storage and quantify their linear and nonlinear contributions. Notably, added freshwater reduces ocean heat and carbon storage by 2,100, and the sum of individual freshwater effects differ from those induced by simultaneous freshwater discharge from both ice sheets. Combined AIS and GrIS freshwater engenders distinct anthropogenic storage anomalies—particularly in the high-latitude Southern Ocean and North Atlantic. From 2080 to 2100, GrIS freshwater exerts primary control on the temporal evolution of global ocean heat storage, while global ocean carbon storage is modulated by the linear AIS and GrIS freshwater impacts. Nonlinear impacts of simultaneous ice sheet discharge have a non-negligible contribution to the evolution of global ocean heat storage. Further, anthropogenic heat changes are realized more quickly in response to ice sheet discharge than anthropogenic carbon. Our results highlight the need to incorporate both ice sheets actively in climate models to accurately project future global climate.

Gorte, Tessa↗
DOE OSTI2025

Improved Understanding of Multicentury Greenland Ice Sheet Response to Strong Warming in the Coupled CESM2‐CISM2 With Regional Grid Refinement

The simulation of ice sheet‐climate interactions, such as surface mass balance fluxes, is sensitive to model grid resolution. Here we simulate the multi‐century evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO 2 increases by 1% yr -1 until quadrupling the pre‐industrial level and then is held fixed). A variable‐resolution (VR) grid with 1/4° regional refinement over the broader Arctic and 1° resolution elsewhere is applied to the atmosphere and land components, and the results are compared with conventional 1° lat‐lon grid simulations to investigate the impact of grid refinement. Compared with the 1° runs, the VR run features a slower rate of surface melt, especially over the western and northern GrIS, where the ice surface slopes gently toward the periphery. This difference pattern originates primarily from higher snow albedo and, thus, weaker albedo feedback in the VR run. The VR grid better captures the CISM ice sheet topography by reducing elevation discrepancies between CAM and CISM and is, therefore, less reliant on the downscaling algorithm, which is known to underestimate albedo gradients. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% less than that of the 1° runs, respectively.

Earth System Model↗
DOE OSTI2025

Improving Projections of Antarctic Ice Sheet Contribution to Sea‐Level Change Through 2300 by Capturing Gravitational, Rotational, and Deformational Effects

Ice-mass change induces regionally varying patterns of sea-level change due to gravitational, rotational, and deformational (GRD) effects, which in turn influence marine-based ice stability in Antarctica. For improved projection of the Antarctic Ice Sheet (AIS), there is a need for including GRD effects in modeling and improving understanding of basin-by-basin sensitivity of ice evolution to GRD effects under a range of climate scenarios. We couple a high-resolution, higher-order ice-sheet model with a 1D global sea-level model that fully captures GRD effects, and simulate ice evolution in Antarctica under the Ice Sheet Model Intercomparison Project for CMIP6 experiments. We perform two sets of coupled simulations incorporating 1D Maxwell solid Earth structure suitable for West and East Antarctica and show that the Amundsen Sea Embayment (ASE) in West Antarctica has the highest sensitivity to GRD effects—in high-emission scenarios, grounding-line retreat accelerates by hundreds of kilometers by 2300 without GRD effects, but GRD effects delay this retreat on a timescale of decades. However, we find that delay times do not show a clear relationship to the strength of climate forcing alone. Furthermore, GRD effects can influence ice-sheet dynamics more than the choice of climate model for a given emissions scenario. In contrast, East Antarctica exhibits minimal sensitivity to GRD effects throughout the study period. These findings underscore the critical role of GRD effects in shaping future West AIS evolution, highlighting the importance of constraining the regional 3D Earth structure and bed topography in West Antarctica, particularly the ASE.

58 GEOSCIENCES↗
DOE OSTI2023

Use of Shallow Ice Core Measurements to Evaluate and Constrain 1980–1990 Global Reanalyses of Ice Sheet Precipitation Rates

Abstract Sea‐level rise (SLR) projections by Earth System Models (ESMs) depend on ice sheet surface mass balances. Accurate, global atmosphere reanalyses would be ideal for providing equilibrated ice sheet model initial conditions in fully coupled ESM simulations. Here we present the first evaluation of 1980–1990 global reanalysis precipitation over Greenland and Antarctica that uses independent observations of net accumulation rates derived from shallow ice cores. Precipitation distributions from both the European Centre for Medium‐Range Weather Forecast's Reanalysis (ERA5) and the Modern‐Era Retrospective Analysis for Research and Applications (MERRA‐2) are highly correlated with contemporaneous co‐located net accumulation rates from Greenland ( r 2 > 0.95) and West Antarctica ( r 2 > 0.7). Three other commonly used reanalyses (WFDE5, CRUNCEP, and GSWP3) exhibit significantly weaker correlations on one or both ice sheets. Our findings imply that ESMs should use ERA5 or MERRA‐2 in data‐forced simulations to validate ice sheet model dynamics and precondition firn for SLR projections.

Schneider, Adam↗
DOE OSTI2026

Antarctic ice sheet model comparison with uncurated geological constraints shows that higher spatial resolution improves deglacial reconstructions

Accurately reconstructing past changes to the shape and volume of the Antarctic ice sheet relies on the use of physically based and thus internally consistent ice sheet modeling, benchmarked against spatially limited geologic data. The challenge in model benchmarking against geologic data is diagnosing whether model-data misfits are the result of an inadequate model, inherently noisy or biased geologic data, and/or incorrect association between modeled quantities and geologic observations. In this work we address this challenge by (i) the development and use of a new model-data evaluation framework applied to an uncurated data set of geologic constraints, and (ii) nested high-spatial-resolution modeling designed to test the hypothesis that model resolution is an important limitation in matching geologic data. While previous approaches to model benchmarking employed highly curated datasets, our approach applies an automated screening and quality control algorithm to an uncurated public dataset of geochronological observations (specifically, cosmogenic-nuclide exposure-age measurements from glacial deposits in ice-free areas). This optimizes data utilization by including more geological constraints, reduces potential interpretive bias, and allows unsupervised assimilation of new data as they are collected. We also incorporate a nested model framework in which high-resolution domains are downscaled from a continent-wide ice sheet model. We highlight the application of this framework by applying these methods to a small ensemble of deglacial ice-sheet model simulations, and demonstrate that the nested approach improves the ability of model simulations to match exposure age data collected from areas of complex topography and ice flow. We develop a range of diagnostic model-data comparison metrics to provide more insight into model performance than possible from a single-valued misfit statistic, showing that different metrics capture different aspects of ice sheet deflation.

Geosciences↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2023

A hybrid deep neural operator/finite element method for ice-sheet modeling

One of the most challenging and consequential problems in climate modeling is to provide probabilistic projections of sea level rise. A large part of the uncertainty of sea level projections is due to uncertainty in ice sheet dynamics. At the moment, accurate quantification of the uncertainty is hindered by the cost of ice sheet computational models. In this work we develop a hybrid approach to approximate existing ice sheet models at a fraction of their cost. Our approach consists of replacing the finite element model for the momentum equations for the ice velocity, the most expensive part of an ice sheet model, with a Deep Operator Network, while we retain a classic finite element discretization for the evolution of the ice thickness. We show that the resulting hybrid model is very accurate and it is an order of magnitude faster than the traditional finite element model. Further, a distinctive feature of the proposed model, compared to other neural network approaches, is that it can handle high-dimensional parameter spaces (parameter fields) such as the basal friction at the bed of the glacier and can therefore be used for generating samples for uncertainty quantification. Further, we study the impact of hyper-parameters, number of unknowns and correlation length of the parameter distribution on the training and accuracy of the Deep Operator Network on a synthetic ice sheet model. We then target the evolution of the Humboldt glacier in Greenland and show that our hybrid model can provide accurate statistics of the glacier mass loss and can be effectively used to accelerate the quantification of uncertainty.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Wind‐Associated Melt Trends and Contrasts Between the Greenland and Antarctic Ice Sheets

Föhn and katabatic winds (downslope winds) can increase ice sheet surface melt, run-off, and ice-shelf vulnerability to hydrofracture and are poorly constrained on the Greenland and Antarctic ice sheets (GIS and AIS). We use regional climate model simulations of the GIS and AIS to quantify and intercompare trends in downslope winds and associated melt since 1960. Results reveal surface melt associated with downslope wind is significant on both the GIS and AIS representing 27.5 ± 4.5% and 19.7 ± 3.8% of total surface melt respectively. Wind-associated melt has decreased 31.8 ± 5.3% on the AIS while total melt decreased 15.4 ± 2.4% due to decreased föhn-induced melt on the Antarctic Peninsula and increasing stratospheric ozone. Wind-associated melt has increased 10.3 ± 2.5% on the GIS, combining with a more positive North Atlantic Oscillation and warmer surface to increase total melt 34 ± 5.8%.

58 GEOSCIENCES↗
DOE OSTI2023

Amundsen Sea Embayment ice-sheet mass-loss predictions to 2050 calibrated using observations of velocity and elevation change

Mass loss from the Amundsen Sea Embayment of the West Antarctic Ice Sheet is a major contributor to global sea-level rise (SLR) and has been increasing over recent decades. Predictions of future SLR are increasingly modelled using ensembles of simulations within which model parameters and external forcings are varied within credible ranges. Accurately reporting the uncertainty associated with these predictions is crucial in enabling effective planning for, and construction of defences against, rising sea levels. Calibrating model simulations against current observations of ice-sheet behaviour enables the uncertainty to be reduced. Here we calibrate an ensemble of BISICLES ice-sheet model simulations of ice loss from the Amundsen Sea Embayment using remotely sensed observations of surface elevation and ice speed. Each calibration type is shown to be capable of reducing the 90% credibility bounds of predicted contributions to SLR by 34 and 43% respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Emulation With Uncertainty Quantification of Regional Sea‐Level Change Caused by the Antarctic Ice Sheet

Abstract Projecting regional sea‐level change under various climate‐change scenarios typically involves running forward simulations of the Earth's gravitational, rotational and deformational (GRD) response to ice‐mass change, which requires substantial computational cost if applied to probabilistic frameworks requiring thousands to millions of samples. Here we build emulators of regional sea‐level change at 27 coastal locations, due to the GRD effects associated with future Antarctic Ice Sheet mass change over the 21st century. The emulators are evaluated against a numerical sea‐level model applied to an ensemble of ice‐sheet model simulations of the Antarctic Ice Sheet through 2100. We build a physics‐based emulator using a recent sensitivity kernel approach and compare it to machine learning based emulators (neural network and conditional variational autoencoder methods). In order to quantify uncertainty, we derive well‐calibrated prediction intervals for regional sea‐level change via split‐conformal inference and linear regression, and show that Monte Carlo dropout does not yield well‐calibrated uncertainties in this instance. We also demonstrate substantial gains in computational efficiency using both the physics‐based emulator and neural networks in comparison to the numerical model for the complete regional sea‐level solution. Overall, we find the physics‐based emulator modestly outperforms the machine learning emulators for this problem.

58 GEOSCIENCES↗