Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19
Assessing the impact of uncertainties in ice-sheet models is a major and challenging issue that needs to be faced by the ice-sheet community to provide more robust and reliable model-based projections of ice-sheet mass balance. In recent years, uncertainty quantification (UQ) has been increasingly u...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , |
Format: | Text |
Language: | English |
Published: |
2022
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-15-1195-2022 https://gmd.copernicus.org/articles/15/1195/2022/ |
id |
ftcopernicus:oai:publications.copernicus.org:gmd97879 |
---|---|
record_format |
openpolar |
spelling |
ftcopernicus:oai:publications.copernicus.org:gmd97879 2023-05-15T14:02:17+02:00 Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 Bulthuis, Kevin Larour, Eric 2022-02-10 application/pdf https://doi.org/10.5194/gmd-15-1195-2022 https://gmd.copernicus.org/articles/15/1195/2022/ eng eng doi:10.5194/gmd-15-1195-2022 https://gmd.copernicus.org/articles/15/1195/2022/ eISSN: 1991-9603 Text 2022 ftcopernicus https://doi.org/10.5194/gmd-15-1195-2022 2022-02-14T17:22:14Z Assessing the impact of uncertainties in ice-sheet models is a major and challenging issue that needs to be faced by the ice-sheet community to provide more robust and reliable model-based projections of ice-sheet mass balance. In recent years, uncertainty quantification (UQ) has been increasingly used to characterize and explore uncertainty in ice-sheet models and improve the robustness of their projections. A typical UQ analysis first involves the (probabilistic) characterization of the sources of uncertainty, followed by the propagation and sensitivity analysis of these sources of uncertainty. Previous studies concerned with UQ in ice-sheet models have generally focused on the last two steps but have paid relatively little attention to the preliminary and critical step of the characterization of uncertainty. Sources of uncertainty in ice-sheet models, like uncertainties in ice-sheet geometry or surface mass balance, typically vary in space and potentially in time. For that reason, they are more adequately described as spatio-(temporal) random fields, which account naturally for spatial (and temporal) correlation. As a means of improving the characterization of the sources of uncertainties for forward UQ analysis within the Ice-sheet and Sea-level System Model (ISSM), we present in this paper a stochastic sampler for Gaussian random fields with Matérn covariance function. The class of Matérn covariance functions provides a flexible model able to capture statistical dependence between locations with different degrees of spatial correlation or smoothness properties. The implementation of this stochastic sampler is based on a notable explicit link between Gaussian random fields with Matérn covariance function and a certain stochastic partial differential equation. Discretization of this stochastic partial differential equation by the finite-element method results in a sparse, scalable and computationally efficient representation known as a Gaussian Markov random field. In addition, spatio-temporal samples can be generated by combining an autoregressive temporal model and the Matérn field. The implementation is tested on a set of synthetic experiments to verify that it captures the desired spatial and temporal correlations well. Finally, we illustrate the interest of this stochastic sampler for forward UQ analysis in an application concerned with assessing the impact of various sources of uncertainties on the Pine Island Glacier, West Antarctica. We find that larger spatial and temporal correlations lengths will both likely result in increased uncertainty in the projections. Text Antarc* Antarctica Ice Sheet Pine Island Pine Island Glacier West Antarctica Copernicus Publications: E-Journals Pine Island Glacier ENVELOPE(-101.000,-101.000,-75.000,-75.000) West Antarctica Geoscientific Model Development 15 3 1195 1217 |
institution |
Open Polar |
collection |
Copernicus Publications: E-Journals |
op_collection_id |
ftcopernicus |
language |
English |
description |
Assessing the impact of uncertainties in ice-sheet models is a major and challenging issue that needs to be faced by the ice-sheet community to provide more robust and reliable model-based projections of ice-sheet mass balance. In recent years, uncertainty quantification (UQ) has been increasingly used to characterize and explore uncertainty in ice-sheet models and improve the robustness of their projections. A typical UQ analysis first involves the (probabilistic) characterization of the sources of uncertainty, followed by the propagation and sensitivity analysis of these sources of uncertainty. Previous studies concerned with UQ in ice-sheet models have generally focused on the last two steps but have paid relatively little attention to the preliminary and critical step of the characterization of uncertainty. Sources of uncertainty in ice-sheet models, like uncertainties in ice-sheet geometry or surface mass balance, typically vary in space and potentially in time. For that reason, they are more adequately described as spatio-(temporal) random fields, which account naturally for spatial (and temporal) correlation. As a means of improving the characterization of the sources of uncertainties for forward UQ analysis within the Ice-sheet and Sea-level System Model (ISSM), we present in this paper a stochastic sampler for Gaussian random fields with Matérn covariance function. The class of Matérn covariance functions provides a flexible model able to capture statistical dependence between locations with different degrees of spatial correlation or smoothness properties. The implementation of this stochastic sampler is based on a notable explicit link between Gaussian random fields with Matérn covariance function and a certain stochastic partial differential equation. Discretization of this stochastic partial differential equation by the finite-element method results in a sparse, scalable and computationally efficient representation known as a Gaussian Markov random field. In addition, spatio-temporal samples can be generated by combining an autoregressive temporal model and the Matérn field. The implementation is tested on a set of synthetic experiments to verify that it captures the desired spatial and temporal correlations well. Finally, we illustrate the interest of this stochastic sampler for forward UQ analysis in an application concerned with assessing the impact of various sources of uncertainties on the Pine Island Glacier, West Antarctica. We find that larger spatial and temporal correlations lengths will both likely result in increased uncertainty in the projections. |
format |
Text |
author |
Bulthuis, Kevin Larour, Eric |
spellingShingle |
Bulthuis, Kevin Larour, Eric Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 |
author_facet |
Bulthuis, Kevin Larour, Eric |
author_sort |
Bulthuis, Kevin |
title |
Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 |
title_short |
Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 |
title_full |
Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 |
title_fullStr |
Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 |
title_full_unstemmed |
Implementation of a Gaussian Markov random field sampler for forward uncertainty quantification in the Ice-sheet and Sea-level System Model v4.19 |
title_sort |
implementation of a gaussian markov random field sampler for forward uncertainty quantification in the ice-sheet and sea-level system model v4.19 |
publishDate |
2022 |
url |
https://doi.org/10.5194/gmd-15-1195-2022 https://gmd.copernicus.org/articles/15/1195/2022/ |
long_lat |
ENVELOPE(-101.000,-101.000,-75.000,-75.000) |
geographic |
Pine Island Glacier West Antarctica |
geographic_facet |
Pine Island Glacier West Antarctica |
genre |
Antarc* Antarctica Ice Sheet Pine Island Pine Island Glacier West Antarctica |
genre_facet |
Antarc* Antarctica Ice Sheet Pine Island Pine Island Glacier West Antarctica |
op_source |
eISSN: 1991-9603 |
op_relation |
doi:10.5194/gmd-15-1195-2022 https://gmd.copernicus.org/articles/15/1195/2022/ |
op_doi |
https://doi.org/10.5194/gmd-15-1195-2022 |
container_title |
Geoscientific Model Development |
container_volume |
15 |
container_issue |
3 |
container_start_page |
1195 |
op_container_end_page |
1217 |
_version_ |
1766272485190270976 |