An official website of the United States government
Abstract. Observations of recent mass loss rates of the West Antarctic Ice Sheet (WAIS) raise concerns about its stability since a collapse would increase global sea levels by several meters. Future projections of these mass loss trends are often estimated using numerical ice sheet models, and recent studies have highlighted the need for models to be benchmarked against present-day observed mass change rates. Here, we present an improved initialization method that optimizes local agreement not only with observations of ice thickness and surface velocity but also with satellite-based estimates of mass change rates. This is achieved by a combination of tuned thermal forcing under the floating ice shelves and friction under the ice sheet. Starting from this improved present-day state, we generate an ensemble of future simulations of Antarctic mass change by varying model physical choices and parameter values while fixing the climate forcing at present-day values. The dynamical response shows slow grounding-line retreat over several centuries, followed by a phase of rapid mass loss over about 200 years with a consistent rate of ∼3 mm GMSL yr−1 (global mean sea level). We find that, for all ensemble members, the Thwaites Glacier and Pine Island Glacier collapse. Our results imply that present-day ocean thermal forcing, if held constant over multiple centuries, may be sufficient to deglaciate large parts of the WAIS, raising global mean sea level by at least a meter.
more »
« less
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 scenarioand 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.
more »
« less
- Award ID(s):
- 2045075
- PAR ID:
- 10418689
- Date Published:
- Journal Name:
- The Cryosphere
- Volume:
- 17
- Issue:
- 4
- ISSN:
- 1994-0424
- Page Range / eLocation ID:
- 1513 to 1543
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st centurynull (Ed.)Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution inresponse to different climate scenarios and assess the mass loss that would contribute tofuture sea level rise. However, there is currently no consensus on estimates of the future massbalance of the ice sheet, primarily because of differences in the representation of physicalprocesses, forcings employed and initial states of ice sheet models. This study presentsresults from ice flow model simulations from 13 international groups focusing on the evolutionof the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet ModelIntercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from theCoupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climatemodel results. Simulations of the Antarctic ice sheet contribution to sea level rise in responseto increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent(SLE) under Representative ConcentrationPathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment withconstant climate conditions and should therefore be added to the mass loss contribution underclimate conditions similar to present-day conditions over the same period. The simulated evolution of theWest 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 and8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighingthe increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelfcollapse, here assumed to be caused by large amounts of liquid water ponding at the surface ofice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without iceshelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, thecalibration of these melt rates based on oceanic conditions taken outside of ice shelf cavitiesand the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario basedon two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared tosimulations done under present-day conditions for the two CMIP5 forcings used and displaylimited mass gain in East Antarctica.more » « less
-
Mass loss from the Amundsen Sea sector of the West Antarctic Ice Sheet has increased in recent decades, suggestive of sus- tained ocean forcing or an ongoing, possibly unstable, response to a past climate anomaly. Lengthening satellite records appear to be incompatible with either process, however, revealing both periodic hiatuses in acceleration and intermittent episodes of thinning. Here we use ocean temperature, salinity, dissolved-oxygen and current measurements taken from 2000 to 2016 near the Dotson Ice Shelf to determine temporal changes in net basal melting. A decadal cycle dominates the ocean record, with melt changing by a factor of about four between cool and warm extremes via a nonlinear relationship with ocean temperature. A warm phase that peaked around 2009 coincided with ice-shelf thinning and retreat of the grounding line, which re-advanced during a post-2011 cool phase. These observations demonstrate how discontinuous ice retreat is linked with ocean variability, and that the strength and timing of decadal extremes is more influential than changes in the longer-term mean state. The non- linear response of melting to temperature change heightens the sensitivity of Amundsen Sea ice shelves to such variability, possibly explaining the vulnerability of the ice sheet in that sector, where subsurface ocean temperatures are relatively high.more » « less
-
Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.more » « less
-
Abstract. Societal adaptation to rising sea levels requires robust projections of the Antarctic Ice Sheet’s retreat, particularly due to ocean-driven basal melting of its fringing ice shelves. Recent advances in ocean models that simulate ice-shelf melting offer an opportunity to reduce uncertainties in ice–ocean interactions. Here, we compare several community-contributed, circum-Antarctic ocean simulations to highlight inter-model differences, evaluate agreement with satellite-derived melt rates, and examine underlying physical processes. All but one simulation use a melting formulation depending on both thermal driving (T ⋆) and friction velocity (u⋆), which together represent the thermal and ocean current forcings at the ice–ocean interface. Simulated melt rates range from 650 to 1277 Gt year−1 (m = 0.45 − 0.91 m year−1), driven by variations in model resolution, parameterisations, and sub-ice shelf circulation. Freeze-to-melt ratios span 0.30 to 30.12 %, indicating large differences in how refreezing is represented. The multi-model mean (MMM) produces an averaged melt rate of 0.60 m year−1 from a net mass loss of 842.99 Gt year−1 (876.03 Gt year−1 melting and 33.05 Gt year−1 refreezing), yielding a freeze-to-melt ratio of 3.92 %. We define a thermo-kinematic melt sensitivity, ζ = m/(T ⋆ u⋆) = 4.82 × 10−5 °C−1 for the MMM, with individual models spanning 2.85 × 10−5 to 19.4 × 10−5 °C−1. Higher melt rates typically occur near grounding zones where both T ⋆ and u⋆ exert roughly equal influence. Because friction velocity is critical for turbulent heat exchange, ice-shelf melting must be characterised by both ocean energetics and thermal forcing. Further work to standardise model setups and evaluation of results against in situ observations and satellite data will be essential for increasing model accuracy, reducing uncertainties, to improve our understanding of ice-shelf–ocean interactions and refine sea-level rise predictions.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5194/tc-17-1513-2023
