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 OSTI2021

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Exploring ice sheet model sensitivity to ocean thermal forcing and basal sliding using the Community Ice Sheet Model (CISM)

Abstract. Multi-meter sea level rise (SLR) is thought to be possible within the next few centuries, with most of the uncertainty originating from the Antarctic land ice contribution. One source of uncertainty relates to the ice sheet model initialization. Since ice sheets have a long response time (compared to other Earth system components such as the atmosphere), ice sheet model initialization methods can have significant impacts on how the ice sheet responds to future forcings. To assess this, we generated 25 different ice sheet spin-ups, using the Community Ice Sheet Model (CISM) at a 4 km resolution. During each spin-up, we varied two key parameters known to impact the sensitivity of the ice sheet to future forcing: one related to the sensitivity of the ice shelf melt rate to ocean thermal forcing (TF) and the other related to the basal friction. The spin-ups all nudge toward observed thickness and enforce a no-advance calving criterion, such that all final spin-up states resemble observations but differ in their melt and friction parameter settings. Each spin-up was then forced with future ocean thermal forcings from 13 different CMIP6 models under the Shared Socioeconomic Pathway (SSP)5-8.5 emissions scenario and modern climatological surface mass balance data. Our results show that the effects of the ice sheet and ocean parameter settings used during the spin-up are capable of impacting multi-century future SLR predictions by as much as 2 m. By the end of this century, the effects of these choices are more modest, but still significant, with differences of up to 0.2 m of SLR. We have identified a combined ocean and ice parameter space that leads to widespread mass loss within 500 years (low friction and high melt rate sensitivity). To explore temperature thresholds, we also ran a synthetically forced CISM ensemble that is focused on the Amundsen region only. Given certain ocean and ice parameter choices, Amundsen mass loss can be triggered with thermal forcing anomalies between 1.5 and 2 ∘C relative to the spin-up. Our results emphasize the critical importance of considering ice sheet and ocean parameter choices during spin-up for SLR predictions and suggest the importance of including glacial isostatic adjustment in ice sheet simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2‐CISM2 Coupled Climate–Ice Sheet Model

Abstract During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uncertainty in Reconstructing Paleo‐Elevation of the Antarctic Ice Sheet From Temperature‐Sensitive Ice Core Records

Paleotemperature reconstructions from ice cores are mixed signals of changes in climate and ice‐surface elevation. A common, temperature‐based paleoaltimetry method suggests these signals can be disentangled by comparing two proxy locations with similar climates. The difference between the records is assumed to be due to elevation, which is estimated by scaling the temperature difference by a lapse rate. We investigate the uncertainty associated with this approach using a case study of the Antarctic Ice Sheet during the Last Glacial Maximum. From an ensemble of climate simulations, we extract modeled temperatures at locations of real ice cores. We find uncertainty on the order of hundreds of meters that results from spatial heterogeneity in non‐adiabatic temperature change, which itself stems in part from elevation‐induced atmospheric circulation change. Our findings suggest that caution is needed when interpreting temperature‐based paleoaltimetry results for ice sheets.

58 GEOSCIENCES↗
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 OSTI2020

Experimental protocol for sea level projections from ISMIP6 stand-alone ice sheet models

Projection of the contribution of ice sheets to sea level change as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the form of simulations from coupled ice sheet–climate models and stand-alone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea level change projections to be performed with stand-alone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice–ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 stand-alone ice sheet simulations, document the experimental framework and implementation, and present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century

Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

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 OSTI2019

Deglaciation of the Greenland and Laurentide ice sheets interrupted by glacier advance during abrupt coolings

Understanding patterns of ice-sheet deglaciation is key for predicting the rate of future ice-sheet melt, yet the processes underlying deglaciation remain elusive. The early Holocene (11.7 ka to 8.2 ka; Greenlandian Stage) represents the most recent period when the Laurentide and Greenland ice sheets underwent large-scale recession. Moreover, this ice-sheet recession occurred under the backdrop of regional temperatures that were similar to or warmer than today, and comparable to those projected for the upcoming centuries. Reconstructing Laurentide and Greenland ice sheet behavior during the early Holocene, and elucidating the mechanisms dictating this behavior may serve as a partial analog for future Greenland ice-sheet change in a warming world. Here, we present 123 new 10 Be surface exposure ages from two sites on Baffin Island and southwestern Greenland that constrain the behavior of the Laurentide and Greenland ice sheets, and an independent alpine glacier during the early Holocene. On Baffin Island, we focus on a unique area where moraines deposited by the Laurentide Ice Sheet rest directly adjacent to moraines deposited by an independent alpine glacier. Sixty-one 10 Be ages reveal that advances and/or stillstands of the Laurentide Ice Sheet and an alpine glacier occurred in unison around 11.8 ka, 10.3 ka, and 9.2 ka. Sixty-two 10 Be ages from southwestern Greenland indicate that the Greenland Ice Sheet margin experienced re-advances or stillstands around 11.6 ka, 10.4 ka, 9.1 ka, 8.1 ka, and 7.3 ka. Our results reveal that ice sheets respond to climate perturbations on the same centennial timescale as small alpine glaciers. Further, we hypothesize that during the warming climate of the early Holocene, freshening of the North Atlantic Ocean induced by ice-sheet melt resulted in regional cooling and brief periods of ice-sheet stabilization superimposed on net glacier recession. These observations point to a negative feedback mechanism inherent to melting ice sheets in the Baffin Bay region that slows ice-sheet recession during intervals of otherwise rapid deglaciation.

58 GEOSCIENCES↗
DOE OSTI2022

A detection of the sea level fingerprint of Greenland Ice Sheet melt

Rapid melting of ice sheets and glaciers drives a unique geometry, or fingerprint, of sea level change. However, the detection of individual fingerprints has been challenging because of sparse observations at high latitudes and the difficulty of disentangling ocean dynamic variability from the signal. Here we predict the fingerprint of Greenland Ice Sheet (GrIS) melt using recent ice mass loss estimates from radar altimetry data and model reconstructions of nearby glaciers and compare this prediction to an independent, altimetry-derived sea surface height trend corrected for ocean dynamic variability in the region adjacent to the ice sheet. A statistically significant correlation between the two fields (P < 0.001) provides an unambiguous observational detection of the near-field sea level fingerprint of recent GrIS melting in our warming world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica

Abstract. The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications.

54 ENVIRONMENTAL SCIENCES↗
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↗