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.

Non-Gaussian Matérn fields with an application to precipitation modeling

High-resolution global precipitation downscaling with latent Gaussian models and non-stationary stochastic partial differential equation structure

Efficient Modeling of Spatial Extremes over Large Geographical Domains

Probabilistic quantitative precipitation field forecasting using a two-stage spatial model

Nonstationary Spatial Gaussian Markov Random Fields

Modeling of Multisite Precipitation Occurrences Using Latent Gaussian-based Multivariate Binary Response Time Series

Deformed SPDE models with an application to spatial modeling of significant wave height

An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach

Projecting Flood-Inducing Precipitation with a Bayesian Analogue Model

Sparse Variational Contaminated Noise Gaussian Process Regression with Applications in Geomagnetic Perturbations Forecasting

A Matern based multivariate Gaussian random process for a consistent model of the horizontal wind components and related variables

Computationally efficient spatial modeling of annual maximum 24 hour precipitation. An application to data from Iceland

Spatial Bayesian hierarchical modeling of precipitation extremes over a large domain

Stochastic Precipitation Generation for the Chesapeake Bay Watershed using Hidden Markov Models with Variational Bayes Parameter Estimation

Latent Gaussian Models for High-Dimensional Spatial Extremes

Nearest-Neighbor Mixture Models for Non-Gaussian Spatial Processes

Exploring the Distributional Properties of the Non-Gaussian Random Field Models

Spatial Modeling of Heavy Precipitation by Coupling Weather Station Recordings and Ensemble Forecasts with Max-Stable Processes

Spatial Modelling of Temperature and Humidity using Systems of Stochastic Partial Differential Equations

Multivariate Matérn Models -- A Spectral Approach

Bayesian hierarchical modeling of extreme hourly precipitation in Norway