A Bayesian semi-parametric hybrid model for spatial extremes with unknown dependence structure

Abstract: The max-stable process is an asymptotically justified model for spatial extremes. In particular, we focus on the hierarchical extreme-value process (HEVP), which is a particular max-stable process that is conducive to Bayesian computing. The HEVP and all max-stable process models are parametric and impose strong assumptions including that all marginal distributions belong to the generalized extreme value family and that nearby sites are asymptotically dependent. We generalize the HEVP by relaxing these assumptions to provide a wider class of marginal distributions via a Dirichlet process prior for the spatial random effects distribution. In addition, we present a hybrid max-mixture model that combines the strengths of the parametric and semi-parametric models. We show that this versatile max-mixture model accommodates both asymptotic independence and dependence and can be fit using standard Markov chain Monte Carlo algorithms. The utility of our model is evaluated in Monte Carlo simulation studies and application to Netherlands wind gust data.