TY - JOUR
T1 - Modelling Sub-daily Precipitation Extremes with the Blended Generalised Extreme Value Distribution
AU - Vandeskog, Silius M.
AU - Martino, Sara
AU - Castro-Camilo, Daniela
AU - Rue, Haavard
N1 - KAUST Repository Item: Exported on 2022-06-17
Acknowledgements: Open access funding provided by orgName. This publication is part of the World of Wild Waters (WoWW) project, which falls under the umbrella of Norwegian University of Science and Technology (NTNU)’s Digital Transformation initiative. We thank Thordis L. Thorarinsdottir and Geir-Arne Fuglstad for helpful discussions.
PY - 2022/6/4
Y1 - 2022/6/4
N2 - A new method is proposed for modelling the yearly maxima of sub-daily precipitation, with the aim of producing spatial maps of return level estimates. Yearly precipitation maxima are modelled using a Bayesian hierarchical model with a latent Gaussian field, with the blended generalised extreme value (bGEV) distribution used as a substitute for the more standard generalised extreme value (GEV) distribution. Inference is made less wasteful with a novel two-step procedure that performs separate modelling of the scale parameter of the bGEV distribution using peaks over threshold data. Fast inference is performed using integrated nested Laplace approximations (INLA) together with the stochastic partial differential equation approach, both implemented in R-INLA. Heuristics for improving the numerical stability of R-INLA with the GEV and bGEV distributions are also presented. The model is fitted to yearly maxima of sub-daily precipitation from the south of Norway and is able to quickly produce high-resolution return level maps with uncertainty. The proposed two-step procedure provides an improved model fit over standard inference techniques when modelling the yearly maxima of sub-daily precipitation with the bGEV distribution. Supplementary materials accompanying this paper appear on-line.
AB - A new method is proposed for modelling the yearly maxima of sub-daily precipitation, with the aim of producing spatial maps of return level estimates. Yearly precipitation maxima are modelled using a Bayesian hierarchical model with a latent Gaussian field, with the blended generalised extreme value (bGEV) distribution used as a substitute for the more standard generalised extreme value (GEV) distribution. Inference is made less wasteful with a novel two-step procedure that performs separate modelling of the scale parameter of the bGEV distribution using peaks over threshold data. Fast inference is performed using integrated nested Laplace approximations (INLA) together with the stochastic partial differential equation approach, both implemented in R-INLA. Heuristics for improving the numerical stability of R-INLA with the GEV and bGEV distributions are also presented. The model is fitted to yearly maxima of sub-daily precipitation from the south of Norway and is able to quickly produce high-resolution return level maps with uncertainty. The proposed two-step procedure provides an improved model fit over standard inference techniques when modelling the yearly maxima of sub-daily precipitation with the bGEV distribution. Supplementary materials accompanying this paper appear on-line.
UR - http://hdl.handle.net/10754/669211
UR - https://link.springer.com/10.1007/s13253-022-00500-7
U2 - 10.1007/s13253-022-00500-7
DO - 10.1007/s13253-022-00500-7
M3 - Article
SN - 1537-2693
JO - JOURNAL OF AGRICULTURAL BIOLOGICAL AND ENVIRONMENTAL STATISTICS
JF - JOURNAL OF AGRICULTURAL BIOLOGICAL AND ENVIRONMENTAL STATISTICS
ER -