Conrad Wasko

Accurate estimation of annual exceedance probabilities (AEPs) of extreme rainfalls through rainfall frequency analysis is a critical step in engineering design for flood mitigation and disaster response. Here we show how the estimation of rainfall frequency curves can be improved by fitting a four‐parameter Kappa distribution to a peaks‐over‐threshold (POT) series. To fit the Kappa distribution to POT data we present a two‐step fitting approach based on maximum likelihood estimation which separately models storm intensity and the arrival frequency. First, a Generalized Pareto distribution (GPA) describing storm intensity is fitted, followed by a Binomial distribution for storm arrivals. We compare the performance of this two‐step Kappa approach to an analogous two‐step Generalized Extreme Value (GEV) approach and to Kappa and GEV distributions fitted to annual maxima series (AMS), using both synthetic and real‐world data representative of Australian climatic conditions. The experiments show that leveraging additional information from the POT series in the two‐step Kappa approach dramatically improves quantile estimation and reduces uncertainty compared to fitting either the GEV or the Kappa distributions to AMS, particularly for rare quantiles. When skew–kurtosis properties of extreme rainfalls are well represented by the Kappa but not the GEV, the two‐step Kappa approach yields unbiased quantiles under extrapolation, while the use of GEV can lead to highly biased estimates. From these results, we believe the two‐step Kappa approach is suitable for both at‐site and regional rainfall frequency analyses as it can accommodate distributions ranging from GPA to GEV without encountering over‐fitting problems generally associated with four‐parameter distributions.