Calibrating an ice sheet model using high-dimensional binary spatial data

Abstract

Rapid retreat of ice in the Amundsen Sea sector of West Antarctica may cause drastic sea level rise, posing significant risks to populations in low-lying coastal regions. Calibration of computer models representing the behavior of the West Antarctic Ice Sheet is key for informative projections of future sea level rise. However, both the relevant observations and the model output are high-dimensional binary spatial data; existing computer model calibration methods are unable to handle such data. Here we present a novel calibration method for computer models whose output is in the form of binary spatial data. To mitigate the computational and inferential challenges posed by our approach, we apply a generalized principal component based dimension reduction method. To demonstrate the utility of our method, we calibrate the PSU3D-ICE model by comparing the output from a 499-member perturbed-parameter ensemble with observations from the Amundsen Sea sector of the ice sheet. Our methods help rigorously characterize the parameter uncertainty even in the presence of systematic data-model discrepancies and dependence in the errors. Our method also helps inform environmental risk analyses by contributing to improved projections of sea level rise from the ice sheets.

References in corpus (4)

• Markov Chain Monte Carlo: Can We Trust the Third Significant Figure?
• Computer model validation with functional output
• Sparse logistic principal components analysis for binary data
• Fast dimension-reduced climate model calibration and the effect of data aggregation

Cited by in corpus (6)

• Uncertainty quantification for computer models with spatial output using calibration-optimal bases
• Efficient calibration for high-dimensional computer model output using basis methods
• Quantifying spatio-temporal boundary condition uncertainty for the North American deglaciation
• Classification of Computer Models with Labelled Outputs
• Bayesian Model Calibration and Sensitivity Analysis for Oscillating Biological Experiments
• Modelling Numerical Systems with Two Distinct Labelled Output Classes