Estimation of confidence regions for random excursion sets with application to large-scale ice-sheet simulations

In many applications, including in evaluations of failure regions and in geophysical hazard assessment, we are interested in evaluating excursion sets, that is, regions in the spatial domain where the response function exceeds some critical value. Determining such excursion sets in the presence of u...

Full description

Bibliographic Details
Main Authors: Bulthuis, Kevin, Pattyn, Frank, Arnst, Maarten
Format: Conference Object
Language:English
Published: 2019
Subjects:
Online Access:https://orbi.uliege.be/handle/2268/238637
https://orbi.uliege.be/bitstream/2268/238637/1/abstract_uncecomp2019.pdf
Description
Summary:In many applications, including in evaluations of failure regions and in geophysical hazard assessment, we are interested in evaluating excursion sets, that is, regions in the spatial domain where the response function exceeds some critical value. Determining such excursion sets in the presence of uncertainties in a model is an interesting problem connected with the theory of random sets. A first issue connected with random excursion sets is in defining confidence regions that can properly represent the uncertainty in the excursion sets . Here, we adopt a definition based on a generalization of the concept of confidence regions based on previous works by [Bolin, 2015] and [French, 2015]. An outer or inner confidence region is defined as a region that contains or is contained in the excursion set with a given level of probability, respectively. Such confidence regions are approximated numerically as optimal subsets within a parametric family of subsets with the appropriate coverage probability, which provides nested approximations for the confidence regions. A second issue, related to this numerical approximation of confidence regions, stems from the numerical approximation of the coverage probability, which may prove challenging for computationally intensive models and small probability levels. Here, we explore methods based on a hybrid surrogate-based approach [Li, 2010] and subset simulation [Au, 2001] to evaluate the coverage probability. We apply this methodology to the evaluation of confidence regions for the retreat of grounded ice in largescale simulation of the Antarctic ice sheet subject to parametric uncertainties. [1] [Au, 2001] Au, S.-K. and Beck, J. L. (2001). Estimation of small failure probabilities in high dimensions by subset simulation. Probabilistic Engineering Mechanics, 16(4):263–277. [2] [Bolin, 2015] Bolin, D. and Lindgren, F. (2015). Excursion and contour uncertainty regions for latent gaussian models. Journal of the Royal Statistical Society: Series B (Statistical Methodology), ...