The generalized Pareto distribution (GPD) is often used in the statistical analysis of climate extremes. For a sample of independent and identically distributed observations, the parameters of the GPD can be estimated by the maximum likelihood (ML) method. In this paper, we drop the assumption of identically distributed random variables. We consider independent observations from GPD distributions having a common shape parameter but possibly an increasing trend in the scale parameter. Such a model, with increasing scale parameter, can be used to describe a trend in the observed extremes as time progresses. Estimating an increasing trend in a distribution parameter is common in the field of isotonic regression. We use ideas and tools from that area to compute ML estimates of the GPD parameters. In a simulation experiment, we show that the iterative convex minorant (ICM) algorithm is much faster than the projected gradient (PG) algorithm. We apply the approach to the daily maxima of the central England temperature (CET) data. A clear positive trend in the GPD scale parameter is found, leading to an increase in the 100-year return level from about 31º in the 1880s to 34º in 2015.@en

An 11 year high-quality radar rainfall data set is used to abstract annual maximum rainfall depths for durations of 15 min to 24 h and area sizes of 6 to 1.7 × 103 km2 for the Netherlands. Generalized extreme value (GEV) distributions are fitted to the annual maxima. A new method is presented to describe the distribution of extreme areal rainfall depths by modeling GEV parameters as a function of both duration and area size. This leads to a semiempirical expression from which quantiles of extreme rainfall depths can be obtained for a chosen duration, area size, and return period. The uncertainties in these quantiles are calculated using the bootstrap method. Radar-based areal reduction factors (ARFs) are derived. These ARFs are comparable to those based on high-density rain gauge networks derived from the literature. It is concluded that radar data, after careful quality control, are suitable to estimate extreme areal rainfall depths.@en

An 11 year high-quality radar rainfall data set is used to abstract annual maximum rainfall depths for durations of 15 min to 24 h and area sizes of 6 to 1.7 × 103 km2 for the Netherlands. Generalized extreme value (GEV) distributions are fitted to the annual maxima. A new method is presented to describe the distribution of extreme areal rainfall depths by modeling GEV parameters as a function of both duration and area size. This leads to a semiempirical expression from which quantiles of extreme rainfall depths can be obtained for a chosen duration, area size, and return period. The uncertainties in these quantiles are calculated using the bootstrap method. Radar-based areal reduction factors (ARFs) are derived. These ARFs are comparable to those based on high-density rain gauge networks derived from the literature. It is concluded that radar data, after careful quality control, are suitable to estimate extreme areal rainfall depths.@en

A hurdle in the peaks-over-threshold approach for analyzing extreme values is the selection of the threshold. A method is developed to reduce this obstacle in the presence of multiple, similar data samples. This is for instance the case in many environmental applications. The idea is to combine threshold selection methods into a regional method. Regionalized versions of the threshold stability and the mean excess plot are presented as graphical tools for threshold selection. Moreover, quantitative approaches based on the bootstrap distribution of the spatially averaged Kolmogorov¿Smirnov and Anderson¿Darling test statistics are introduced. It is demonstrated that the proposed regional method leads to an increased sensitivity for too low thresholds, compared to methods that do not take into account the regional information. The approach can be used for a wide range of univariate threshold selection methods. We test the methods using simulated data and present an application to rainfall data from the Dutch water board Vallei en Veluwe.@en

A hurdle in the peaks-over-threshold approach for analyzing extreme values is the selection of the threshold. A method is developed to reduce this obstacle in the presence of multiple, similar data samples. This is for instance the case in many environmental applications. The idea is to combine threshold selection methods into a regional method. Regionalized versions of the threshold stability and the mean excess plot are presented as graphical tools for threshold selection. Moreover, quantitative approaches based on the bootstrap distribution of the spatially averaged Kolmogorov¿Smirnov and Anderson¿Darling test statistics are introduced. It is demonstrated that the proposed regional method leads to an increased sensitivity for too low thresholds, compared to methods that do not take into account the regional information. The approach can be used for a wide range of univariate threshold selection methods. We test the methods using simulated data and present an application to rainfall data from the Dutch water board Vallei en Veluwe.@en