An official website of the United States government
Ice shelf basal melting is the primary mechanism driving mass loss from the Antarctic Ice Sheet, yet it is unknown how the localized melt enhancement from subglacial discharge will affect future Antarctic glacial retreat. We develop a parameterization of ice shelf basal melt that accounts for both ocean and subglacial discharge forcing and apply it in future projections of Denman and Scott Glaciers, East Antarctica, through 2300. In forward simulations, subglacial discharge accelerates the onset of retreat of these systems into the deepest continental trench on Earth by 25 years. During this retreat, Denman Glacier alone contributes 0.33 millimeters per year to global sea level rise, comparable to half of the contemporary sea level contribution of the entire Antarctic Ice Sheet. Our results stress the importance of resolving complex interactions between the ice, ocean, and subglacial environments in future Antarctic Ice Sheet projections.
more »
« less
Statistical emulation of a perturbed basal melt ensemble of an ice sheet model to better quantify Antarctic sea level rise uncertainties
Abstract. Antarctic ice shelves are vulnerable to warming ocean temperatures, and some have already begun thinning in response to increased basal melt rates.Sea level is therefore expected to rise due to Antarctic contributions, but uncertainties in its amount and timing remain largely unquantified. In particular, there is substantial uncertainty in future basal melt rates arising from multi-model differences in thermal forcing and how melt rates depend on that thermal forcing. To facilitate uncertainty quantification in sea level rise projections, we build, validate, and demonstrate projections from a computationally efficient statistical emulator of a high-resolution (4 km) Antarctic ice sheet model, the Community Ice Sheet Model version 2.1. The emulator is trained to a large (500-member) ensemble of 200-year-long 4 km resolution transient ice sheet simulations, whereby regional basal melt rates are perturbed by idealized (yet physically informed) trajectories. The main advantage of our emulation approach is that by sampling a wide range of possible basal melt trajectories, the emulator can be used to (1) produce probabilistic sea level rise projections over much larger Monte Carlo ensembles than are possible by direct numerical simulation alone, thereby providing better statistical characterization of uncertainties, and (2) predict the simulated ice sheet response under differing assumptions about basal melt characteristics as new oceanographic studies are published, without having to run additional numerical ice sheet simulations. As a proof of concept, we propagate uncertainties about future basal melt rate trajectories, derived from regional ocean models, to generate probabilistic sea level rise estimates for 100 and 200 years into the future.
more »
« less
- Award ID(s):
- 2045075
- PAR ID:
- 10418764
- Date Published:
- Journal Name:
- The Cryosphere
- Volume:
- 15
- Issue:
- 6
- ISSN:
- 1994-0424
- Page Range / eLocation ID:
- 2683 to 2699
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Previous studies have interpreted Last Interglacial (LIG;∼129–116 ka) sea‐level estimates in multiple different ways to calibrate projections of future Antarctic ice‐sheet (AIS) mass loss and associated sea‐level rise. This study systematically explores the extent to which LIG constraints could inform future Antarctic contributions to sea‐level rise. We develop a Gaussian process emulator of an ice‐sheet model to produce continuous probabilistic projections of Antarctic sea‐level contributions over the LIG and a future high‐emissions scenario. We use a Bayesian approach conditioning emulator projections on a set of LIG constraints to find associated likelihoods of model parameterizations. LIG estimates inform both the probability of past and future ice‐sheet instabilities and projections of future sea‐level rise through 2150. Although best‐available LIG estimates do not meaningfully constrain Antarctic mass loss projections or physical processes until 2060, they become increasingly informative over the next 130 years. Uncertainties of up to 50 cm remain in future projections even if LIG Antarctic mass loss is precisely known (±5 cm), indicating that there is a limit to how informative the LIG could be for ice‐sheet model future projections. The efficacy of LIG constraints on Antarctic mass loss also depends on assumptions about the Greenland ice sheet and LIG sea‐level chronology. However, improved field measurements and understanding of LIG sea levels still have potential to improve future sea‐level projections, highlighting the importance of continued observational efforts.more » « less
-
Abstract Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica.more » « less
-
Abstract Mass loss from the Antarctic ice sheet is projected to continue over the coming century. The resultant sea level change will have a regional pattern that evolves over time as the ocean adjusts. Accurate estimation of this evolution is crucial for local communities. Current state-of-the-art climate models typically do not couple ice sheets to the atmosphere–ocean system, and the impact of ice sheet melt has often been studied by injecting meltwater at the model ocean surface. However, observational evidence suggests that most Antarctic meltwater enters the ocean at depth through ice shelf basal melt. A previous study has demonstrated that the regional sea level pattern at a given time depends on meltwater injection depth. Here, we introduce a 2.5-layer model to investigate this dependence and develop a theory for the associated adjustment mechanisms. We find mechanisms consistent with previous literature on the ocean adjustment to changes in forcing, whereby a slower Rossby wave response off the eastern boundary follows a fast response from the western boundary current and Kelvin waves. We demonstrate that faster baroclinic Rossby waves near the surface than at depth explain the injection depth dependence of the adjustment in the 2.5-layer model. The identified Rossby wave mechanism may contribute to the dependence of the ocean’s transient adjustment on meltwater injection depth in more complex models. This work highlights processes that could cause errors in the projection of the time-varying pattern of sea level rise using surface meltwater input to represent Antarctica’s freshwater forcing. Significance StatementSea level rise is expected to be larger in some locations than in others. Accurate projections of the pattern of sea level change, which changes in time as the ocean adjusts, are essential information for local communities. One of the factors that leads to uncertainty in the local sea level change due to Antarctic melt is the depth at which this meltwater is input into an ocean model. We propose a mechanism for a faster response of sea level around the basin when meltwater is injected at the ocean surface compared to when it is injected well below the surface. This mechanism has implications for projections of the regional sea level response to Antarctic melt.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5194/tc-15-2683-2021
