TY - JOUR
T1 - Probabilistic calibration of a Greenland ice sheet model using spatially resolved synthetic observations
T2 - Toward projections of ice mass loss with uncertainties
AU - Chang, W.
AU - Applegate, P. J.
AU - Haran, M.
AU - Keller, K.
N1 - Publisher Copyright:
© Author(s) 2014. CC Attribution 3.0 License.
PY - 2014/9/5
Y1 - 2014/9/5
N2 - Computer models of ice sheet behavior are important tools for projecting future sea level rise. The simulated modern ice sheets generated by these models differ markedly as input parameters are varied. To ensure accurate ice sheet mass loss projections, these parameters must be constrained using observational data. Which model parameter combinations make sense, given observations? Our method assigns probabilities to parameter combinations based on how well the model reproduces the Greenland Ice Sheet profile. We improve on the previous state of the art by accounting for spatial information and by carefully sampling the full range of realistic parameter combinations, using statistically rigorous methods. Specifically, we estimate the joint posterior probability density function of model parameters using Gaussian process-based emulation and calibration. This method is an important step toward calibrated probabilistic projections of ice sheet contributions to sea level rise, in that it uses data-model fusion to learn about parameter values. This information can, in turn, be used to make projections while taking into account various sources of uncertainty, including parametric uncertainty, data-model discrepancy, and spatial correlation in the error structure. We demonstrate the utility of our method using a perfect model experiment, which shows that many different parameter combinations can generate similar modern ice sheet profiles. This result suggests that the large divergence of projections from different ice sheet models is partly due to parametric uncertainty. Moreover, our method enables insight into ice sheet processes represented by parameter interactions in the model.
AB - Computer models of ice sheet behavior are important tools for projecting future sea level rise. The simulated modern ice sheets generated by these models differ markedly as input parameters are varied. To ensure accurate ice sheet mass loss projections, these parameters must be constrained using observational data. Which model parameter combinations make sense, given observations? Our method assigns probabilities to parameter combinations based on how well the model reproduces the Greenland Ice Sheet profile. We improve on the previous state of the art by accounting for spatial information and by carefully sampling the full range of realistic parameter combinations, using statistically rigorous methods. Specifically, we estimate the joint posterior probability density function of model parameters using Gaussian process-based emulation and calibration. This method is an important step toward calibrated probabilistic projections of ice sheet contributions to sea level rise, in that it uses data-model fusion to learn about parameter values. This information can, in turn, be used to make projections while taking into account various sources of uncertainty, including parametric uncertainty, data-model discrepancy, and spatial correlation in the error structure. We demonstrate the utility of our method using a perfect model experiment, which shows that many different parameter combinations can generate similar modern ice sheet profiles. This result suggests that the large divergence of projections from different ice sheet models is partly due to parametric uncertainty. Moreover, our method enables insight into ice sheet processes represented by parameter interactions in the model.
UR - http://www.scopus.com/inward/record.url?scp=84907008575&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84907008575&partnerID=8YFLogxK
U2 - 10.5194/gmd-7-1933-2014
DO - 10.5194/gmd-7-1933-2014
M3 - Article
AN - SCOPUS:84907008575
SN - 1991-959X
VL - 7
SP - 1933
EP - 1943
JO - Geoscientific Model Development
JF - Geoscientific Model Development
IS - 5
ER -