Abstract

Determining the future evolution of the Antarctic Ice Sheet is critical for understanding and narrowing the large existing uncertainties in century-scale global mean sea level rise (SLR) projections. One of the most significant glaciers/ice streams in Antarctica, Thwaites Glacier, is at risk of destabilization and, if destabilized, has the potential to be the largest regional-scale contributor of SLR on Earth. This is because Thwaites Glacier is vulnerable to the marine ice sheet instability, as its grounding line is significantly influenced by ocean-driven basal melting rates, and its bedrock topography retrogrades into kilometer deep troughs. In this study, we investigate how bedrock topography features influence the grounding line migration beneath Thwaites Glacier when extreme ocean-driven basal melt rates are applied. Specifically, we design experiments using the Ice-Sheet and Sea-level System Model (ISSM) to quantify the SLR projection uncertainty due to reported errors in the current bedrock topography maps that are often used by ice sheet models. We find that spread in model estimates of sea level rise contribution from Thwaites glacier due to the reported bedrock topography error could be as large as 21.9 cm after 200 years of extreme ocean warming. Next, we perturb the bedrock topography beneath Thwaites Glacier using wavelet decomposition techniques to introduce realistic noise (within error). We explore the model space with multiple realizations of noise to quantify what spatial and vertical resolutions in bedrock topography are required to minimize the uncertainty in our 200-year experiment. We conclude that at least a 2 km spatial and 8 m vertical resolution would independently constrain possible SLR to ±2 cm over 200 years, fulfilling requirements outlined by the 2017 Decadal Survey for Earth Science. Lastly, we perform an ensemble of simulations to determine in which regions our model of Thwaites Glacier is most sensitive to perturbations in bedrock topography. Our results suggest that the retreat of the grounding line is most sensitive to bedrock topography in proximity to the grounding line's initial position. Additionally, we find that the location and amplitude of the bedrock perturbation is more significant than its sharpness and shape. Overall, these findings inform and benchmark observational requirements for future missions that will measure ice sheet bedrock topography, not only in the case of Thwaites Glacier but for Antarctica on the continental scale.

Highlights

  • The future of the West Antarctic Ice Sheet (WAIS) is known to be one of the largest sources of uncertainties in global mean sea level rise (SLR) on a century time scale (Schlegel et al, 2018; Yu et al, 2019)

  • We find that spread in model estimates of sea level rise contribution from Thwaites glacier due to the reported bedrock topography error could be as large as 21.9 cm after 200 years of extreme ocean warming

  • It is important to note that, in comparison to variations in spatial resolution, we find that vertical perturbation in bedrock topography contributes more significantly to uncertainty in ice sheet model estimates of SLR contribution from Thwaites Glacier, especially since even a meter of perturbed vertical bedrock elevation has the potential to alter model results 410 significantly

Read more

Summary

Introduction

The future of the West Antarctic Ice Sheet (WAIS) is known to be one of the largest sources of uncertainties in global mean sea level rise (SLR) on a century time scale (Schlegel et al, 2018; Yu et al, 2019). To characterize how the behavior of various ice shelves throughout WAIS may evolve in the future, the cryosphere community relies heavily on numerical modeling and simulations Such tools allow physically-based prediction and quantification of the sensitivity of grounding line evolution. Recent model-based studies suggest that century-scale uncertainty in SLR potential under different basal melting rate scenarios has a significant spread (Yu et al, 2018; Schlegel et al, 2018) Together, this all suggests that uncertainty in estimates of ocean forcing may remain a significant source (and 65 perhaps the largest) of uncertainty in ice sheet model simulations of Thwaites Glacier into the near future. We design a suite of ice sheet model experiments to investigate how known errors in bedrock topography affect 85 200-year simulations of the response of Thwaites Glacier to an extreme increase in ocean-driven basal melt rates. We perform an uncertainty quantification (UQ) sampling experiment to simulate probabilistic model outcomes followed by a sensitivity test to locate regions where the bedrock topography of Thwaites plays the most significant role in determining grounding line evolution (Sect. 6)

Model Description We use the Ice-Sheet and Sea-level System
Mesh and Boundary Conditions Initialization
Stress Balance Approximation
Grounding Line Migration and Forcing Basal melting rates are calculated using the
Wavelet Decomposition Setup
Experiment 2
Experiment 3
Experiment 4
Experiment 5
Discussion
Conclusion
Initialization Algorithm
Figure Code Combination
Findings
10 Notes
Code and Data Availability
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call