English

Latent Gaussian Models for High-Dimensional Spatial Extremes

Methodology 2021-10-07 v1 Computation

Abstract

In this chapter, we show how to efficiently model high-dimensional extreme peaks-over-threshold events over space in complex non-stationary settings, using extended latent Gaussian Models (LGMs), and how to exploit the fitted model in practice for the computation of long-term return levels. The extended LGM framework assumes that the data follow a specific parametric distribution, whose unknown parameters are transformed using a multivariate link function and are then further modeled at the latent level in terms of fixed and random effects that have a joint Gaussian distribution. In the extremal context, we here assume that the data level distribution is described in terms of a Poisson point process likelihood, motivated by asymptotic extreme-value theory, and which conveniently exploits information from all threshold exceedances. This contrasts with the more common data-wasteful approach based on block maxima, which are typically modeled with the generalized extreme-value (GEV) distribution. When conditional independence can be assumed at the data level and latent random effects have a sparse probabilistic structure, fast approximate Bayesian inference becomes possible in very high dimensions, and we here present the recently proposed inference approach called "Max-and-Smooth", which provides exceptional speed-up compared to alternative methods. The proposed methodology is illustrated by application to satellite-derived precipitation data over Saudi Arabia, obtained from the Tropical Rainfall Measuring Mission, with 2738 grid cells and about 20 million spatio-temporal observations in total. Our fitted model captures the spatial variability of extreme precipitation satisfactorily and our results show that the most intense precipitation events are expected near the south-western part of Saudi Arabia, along the Red Sea coastline.

Keywords

Cite

@article{arxiv.2110.02680,
  title  = {Latent Gaussian Models for High-Dimensional Spatial Extremes},
  author = {Arnab Hazra and Raphaël Huser and Árni V. Jóhannesson},
  journal= {arXiv preprint arXiv:2110.02680},
  year   = {2021}
}

Comments

This paper (after peer-review) will be a book chapter of the forthcoming book entitled "Statistical modeling using latent Gaussian models - with applications in geophysics and environmental sciences", expected to be published by Springer in 2022

R2 v1 2026-06-24T06:39:58.915Z