Spatial probabilistic calibration of a high-resolution Amundsen Sea Embayment ice sheet model with satellite altimeter data

• Published in: The Cryosphere
• Main Authors: A. Wernecke, T. L. Edwards, I. J. Nias, P. B. Holden, N. R. Edwards
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2020
• Subjects: geo; envir; Amundsen Sea; Antarc*; Antarctica; Ice Sheet; The Cryosphere
• Online Access: https://doi.org/10.5194/tc-14-1459-2020; https://www.the-cryosphere.net/14/1459/2020/tc-14-1459-2020.pdf; https://doaj.org/article/e9892cc212334400bd02b5b77f4def84

Probabilistic predictions of the sea level contribution from Antarctica often have large uncertainty intervals. Calibration of model simulations with observations can reduce uncertainties and improve confidence in projections, particularly if this exploits as much of the available information as possible (such as spatial characteristics), but the necessary statistical treatment is often challenging and can be computationally prohibitive. Ice sheet models with sufficient spatial resolution to resolve grounding line evolution are also computationally expensive. Here we address these challenges by adopting and comparing dimension-reduced calibration approaches based on a principal component decomposition of the adaptive mesh model BISICLES. The effects model parameters have on these principal components are then gathered in statistical emulators to allow for smooth probability density estimates. With the help of a published perturbed parameter ice sheet model ensemble of the Amundsen Sea Embayment (ASE), we show how the use of principal components in combination with spatially resolved observations can improve probabilistic calibrations. In synthetic model experiments (calibrating the model with altered model results) we can identify the correct basal traction and ice viscosity scaling parameters as well as the bedrock map with spatial calibrations. In comparison a simpler calibration against an aggregated observation, the net sea level contribution, imposes only weaker constraints by allowing a wide range of basal traction and viscosity scaling factors. Uncertainties in sea level rise contribution of 50-year simulations from the current state of the ASE can be reduced with satellite observations of recent ice thickness change by nearly 90 %; median and 90 % confidence intervals are 18.9 [13.9, 24.8] mm SLE (sea level equivalent) for the proposed spatial calibration approach, 16.8 [7.7, 25.6] mm SLE for the net sea level calibration and 23.1 [−8.4, 94.5] mm SLE for the uncalibrated ensemble. The spatial model ...