Modelling of hydrological extremes and drought modelling in particular has received much attention over recent decades. The main aim of this study is to apply a statistical model for drought estimation (in this case deficit volume) using extreme value theory and the index-flood method and to reduce the uncertainties in estimation of drought event return levels. Deficit volumes for 133 catchments in the Czech Republic (1901–2015) were simulated by hydrological model BILAN. The validation of severity, intensity and length of simulated drought events revealed good match with the available observed data. To estimate return levels of the deficit volumes, it is assumed (in accord with the index-flood method), that the deficit volumes within a homogeneous region are identically distributed after scaling with a site-specific factor. The parameters of the scaled regional distribution are estimated using L-moments. The goodness-of-fit of the statistical model is assessed by Anderson–Darling test. For the estimation of critical values, sampling methods allowing for handling of years without drought were used. It is shown, that the index-flood model with a Generalized Pareto distribution performs well and substantially reduces the uncertainty related to the estimation of the shape parameter and of the large deficit volume quantiles
