Predictive Inference for Spatio-temporal Precipitation Data and Its Extremes

TL;DR

This paper introduces a Bayesian hierarchical model for spatio-temporal precipitation data that effectively captures the entire distribution, including extremes, by accounting for skewness, clustering, and variability, with applications to Australian rainfall.

Contribution

The novel model uses all available data efficiently, modeling the full rainfall distribution and its extremes, unlike traditional extreme value theory approaches.

Findings

Model accurately captures precipitation distribution and extremes.

Provides spatial and temporal predictions with uncertainty estimates.

Effectively models skewness and clustering in environmental data.

Abstract

Modelling of precipitation and its extremes is important for urban and agriculture planning purposes. We present a method for producing spatial predictions and measures of uncertainty for spatio-temporal data that is heavy-tailed and subject to substaintial skewness which often arise in measurements of many environmental processes, and we apply the method to precipitation data in south-west Western Australia. A generalised hyperbolic Bayesian hierarchical model is constructed for the intensity, frequency and duration of daily precipitation, including the extremes. Unlike models based on extreme value theory, which only model maxima of finite-sized blocks or exceedances above a large threshold, the proposed model uses all the data available efficiently, and hence not only fits the extremes but also models the entire rainfall distribution. It captures spatial and temporal clustering, as well as spatially and temporally varying volatility and skewness. The model assumes that the regional precipitation is driven by a latent process characterised by geographical and climatological covariates. Effects not fully described by the covariates are captured by spatial and temporal structure in the hierarchies. Inference is provided by MCMC using a Metropolis-Hastings algorithm and spatial interpolation method, which provide a natural approach for estimating uncertainty. Similarly both spatial and temporal predictions with uncertainty can be produced with the model.