Estimation of spatial max-stable models using threshold exceedances

Parametric inference for spatial max-stable processes is difficult since the related likelihoods are unavailable. A composite likelihood approach based on the bivariate distribution of block maxima has been recently proposed. However modeling block maxima is a wasteful approach provided that other information is available. Moreover an approach based on block maxima, typically annual, is unable to take into account the fact that maxima occur or not simultaneously. If time series of, say, daily data are available, then estimation procedures based on exceedances of a high threshold could mitigate such problems. We focus on two approaches for composing likelihoods based on pairs of exceedances. The first one comes from the tail approximation for bivariate distribution proposed by Ledford and Tawn (Biometrika 83:169-187, 1996) when both pairs of observations exceed the fixed threshold. The second one uses the bivariate extension (Rootz,n and Tajvidi in Bernoulli 12:917-930, 2006) of the generalized Pareto distribution which allows to model exceedances when at least one of the components is over the threshold. The two approaches are compared through a simulation study where both processes in a domain of attraction of a max-stable process and max-stable processes are successively considered as time replications, according to different degrees of spatial dependency. Results put forward how the nature of the time replications influences the bias of estimations and highlight the choice of each approach regarding to the strength of the spatial dependencies and the threshold choice.