Evidences of relationships between statistics of rainfall extremes and mean annual precipitation: an application for design-storm estimation in northern central Italy

International audienceSeveral hydrological analyses need to be founded on a reliable estimate of the design storm, which is the expected rainfall depth corresponding to a given duration and probability of occurrence, usually expressed in terms of return period. The annual series of precipitation maxima for storm duration ranging from 15 min to 1 day are observed at a dense network of raingauges sited in northern central Italy are statistically analyzed using an approach based on L-moments. The study investigates the statistical properties of rainfall extremes and identifies important relationships between these properties and the mean annual precipitation (MAP). On the basis of these relationships, we develop a regional model for estimating the rainfall depth for a given storm duration and recurrence interval in any location of the study region. The reliability of the regional model is assessed through Monte Carlo simulations. The results are relevant given that the proposed model is able to reproduce the statistical properties of rainfall extremes observed for the study region

Relationships between statistics of rainfall extremes and mean annual precipitation: an application for design-storm estimation in northern central Italy

Several hydrological analyses need to be founded on a reliable estimate of the design storm, which is the expected rainfall depth corresponding to a given duration and probability of occurrence, usually expressed in terms of return period. The annual series of precipitation maxima for storm duration ranging from 15 min to 1 day, observed at a dense network of raingauges sited in northern central Italy, are analyzed using an approach based on L-moments. The analysis investigates the statistical properties of rainfall extremes and detects significant relationships between these properties and the mean annual precipitation (MAP). On the basis of these relationships, we developed a regional model for estimating the rainfall depth for a given storm duration and recurrence interval in any location of the study region. The applicability of the regional model was assessed through Monte Carlo simulations. The uncertainty of the model for ungauged sites was quantified through an extensive cross-validation

A look at the links between drainage density and flood statistics

Abstract. We investigate the links between the drainage density of a river basin and selected flood statistics, namely, mean, standard deviation, coefficient of variation and coefficient of skewness of annual maximum series of peak flows. The investigation is carried out through a three-stage analysis. First, a numerical simulation is performed by using a spatially distributed hydrological model in order to highlight how flood statistics change with varying drainage density. Second, a conceptual hydrological model is used in order to analytically derive the dependence of flood statistics on drainage density. Third, real world data from 44 watersheds located in northern Italy were analysed. The three-level analysis seems to suggest that a critical value of the drainage density exists for which a minimum is attained in both the coefficient of variation and the absolute value of the skewness coefficient. Such minima in the flood statistics correspond to a minimum of the flood quantile for a given exceedance probability (i.e., recurrence interval). Therefore, the results of this study may provide useful indications for flood risk assessment in ungauged basins

Development and assessment of uni- and multivariable flood loss models for Emilia-Romagna (Italy)

Flood loss models are one important source of uncertainty in flood risk assessments. Many countries experience sparseness or absence of comprehensive high-quality flood loss data, which is often rooted in a lack of protocols and reference procedures for compiling loss datasets after flood events. Such data are an important reference for developing and validating flood loss models. We consider the Secchia River flood event of January 2014, when a sudden levee breach caused the inundation of nearly 52km2 in northern Italy. After this event local authorities collected a comprehensive flood loss dataset of affected private households including building footprints and structures and damages to buildings and contents. The dataset was enriched with further information compiled by us, including economic building values, maximum water depths, velocities and flood durations for each building. By analyzing this dataset we tackle the problem of flood damage estimation in Emilia-Romagna (Italy) by identifying empirical uni- and multivariable loss models for residential buildings and contents. The accuracy of the proposed models is compared with that of several flood damage models reported in the literature, providing additional insights into the transferability of the models among different contexts. Our results show that (1) even simple univariable damage models based on local data are significantly more accurate than literature models derived for different contexts; (2) multivariable models that consider several explanatory variables outperform univariable models, which use only water depth. However, multivariable models can only be effectively developed and applied if sufficient and detailed information is available

Estimating the flood frequency distribution at seasonal and annual time scales

Abstract. We propose an original approach to infer the flood frequency distribution at seasonal and annual time scale. Our purpose is to estimate the peak flow that is expected for an assigned return period T, independently of the season in which it occurs (i.e. annual flood frequency regime), as well as in different selected sub-yearly periods (i.e. seasonal flood frequency regime). While a huge literature exists on annual flood frequency analysis, few studies have focused on the estimation of seasonal flood frequencies despite the relevance of the issue, for instance when scheduling along the months of the year the construction phases of river engineering works directly interacting with the active river bed, like for instance dams. An approximate method for joint frequency analysis is presented here that guarantees consistency between fitted annual and seasonal distributions, i.e. the annual cumulative distribution is the product of the seasonal cumulative distribution functions, under the assumption of independence among floods in different seasons. In our method the parameters of the seasonal frequency distributions are fitted by maximising an objective function that accounts for the likelihoods of both seasonal and annual peaks. In contrast to previous studies, our procedure is conceived to allow the users to introduce subjective weights to the components of the objective function in order to emphasize the fitting of specific seasons or of the annual peak flow distribution. An application to the time series of the Blue Nile daily flows at the Sudan–Ethiopia border is presented

Prediction of streamflow regimes over large geographical areas: interpolated flow–duration curves for the Danube region

ABSTRACTFlow–duration curves (FDCs) are essential to support decisions on water resources management, and their regionalization is fundamental for the assessment of ungauged basins. In comparison with calibrated rainfall–runoff models, statistical methods provide data-driven estimates representing a useful benchmark. The objective of this work is the interpolation of FDCs from ~500 discharge gauging stations in the Danube. To this aim we use total negative deviation top-kriging (TNDTK), as multi-regression models are shown to be unsuitable for representing FDCs across all durations and sites. TNDTK shows a high accuracy for the entire Danube region, with overall Nash-Sutcliffe efficiency values computed in a leave-p-out cross-validation scheme (p equal to one site, one-third and half of the sites), all above 0.88. A reliability measure based on kriging variance is attached to each interpolated FDC at ~4000 prediction nodes. The GIS layer of regionalized FDCs is made available for broader use in the region

Streamflow data availability in Europe: a detailed dataset of interpolated flow-duration curves

For about 24 000 river basins across Europe, we provide a continuous representation of the stream-flow regime in terms of empirical flow-duration curves (FDCs), which are key signatures of the hydrological behaviour of a catchment and are widely used for supporting decisions on water resource management as well as for assessing hydrologic change. In this study, FDCs are estimated by means of the geostatistical procedure termed total negative deviation top-kriging (TNDTK), starting from the empirical FDCs made available by the Joint Research Centre of the European Commission (DG-JRC) for about 3000 discharge measurement stations across Europe. Consistent with previous studies, TNDTK is shown to provide high accuracy for the entire study area, even with different degrees of reliability, which varies significantly over the study area. In order to provide this kind of information site by site, together with the estimated FDCs, for each catchment we provide indicators of the accuracy and reliability of the performed large-scale geostatistical prediction. The dataset is freely available at the PANGAEA open-access library (Data Publisher for Earth & Environmental Science) at https://doi.org/10.1594/PANGAEA.938975 (Persiano et al., 2021b)

Machine-learning blends of geomorphic descriptors: value and limitations for flood hazard assessment across large floodplains

Recent literature shows several examples of simplified approaches that perform flood hazard (FH) assessment and mapping across large geographical areas on the basis of fast-computing geomorphic descriptors. These approaches may consider a single index (univariate) or use a set of indices simultaneously (multivariate). What is the potential and accuracy of multivariate approaches relative to univariate ones? Can we effectively use these methods for extrapolation purposes, i.e., FH assessment outside the region used for setting up the model? Our study addresses these open problems by considering two separate issues: (1) mapping flood-prone areas and (2) predicting the expected water depth for a given inundation scenario. We blend seven geomorphic descriptors through decision tree models trained on target FH maps, referring to a large study area (∼ 105 km2). We discuss the potential of multivariate approaches relative to the performance of a selected univariate model and on the basis of multiple extrapolation experiments, where models are tested outside their training region. Our results show that multivariate approaches may (a) significantly enhance flood-prone area delineation (accuracy: 92%) relative to univariate ones (accuracy: 84%), (b) provide accurate predictions of expected inundation depths (determination coefficient ∼0.7), and (c) produce encouraging results in extrapolation

Regional parent flood frequency distributions in Europe – Part 1: Is the GEV model suitable as a pan-European parent?

Abstract. This study addresses the question of the existence of a parent flood frequency distribution on a European scale. A new database of L-moment ratios of flood annual maximum series (AMS) from 4105 catchments was compiled by joining 13 national data sets. Simple exploration of the database presents the generalized extreme value (GEV) distribution as a potential pan-European flood frequency distribution, being the three-parameter statistical model that with the closest resemblance to the estimated average of the sample L-moment ratios. Additional Monte Carlo simulations show that the variability in terms of sample skewness and kurtosis present in the data is larger than in a hypothetical scenario where all the samples were drawn from a GEV model. Overall, the generalized extreme value distribution fails to represent the kurtosis dispersion, especially for the longer sample lengths and medium to high skewness values, and therefore may be rejected in a statistical hypothesis testing framework as a single pan-European parent distribution for annual flood maxima. The results presented in this paper suggest that one single statistical model may not be able to fit the entire variety of flood processes present at a European scale, and presents an opportunity to further investigate the catchment and climatic factors controlling European flood regimes and their effects on the underlying flood frequency distributions

Smooth regional estimation of low-flow indices: physiographical space based interpolation and top-kriging

Recent studies highlight that spatial interpolation techniques of point data can be effectively applied to the problem of regionalization of hydrometric information. This study compares two innovative interpolation techniques for the prediction of low-flows in ungauged basins. The first one, named Physiographical-Space Based Interpolation (PSBI), performs the spatial interpolation of the desired streamflow index (e.g., annual streamflow, low-flow index, flood quantile, etc.) in the space of catchment descriptors. The second technique, named Topological kriging or Top-kriging, predicts the variable of interest along river networks taking both the area and nested nature of catchments into account. PSBI and Top-kriging are applied for the regionalization of <i>Q</i><sub>355</sub> (i.e., a low-flow index that indicates the streamflow that is equalled or exceeded 355 days in a year, on average) over a broad geographical region in central Italy, which contains 51 gauged catchments. The two techniques are cross-validated through a leave-one-out procedure at all available gauges and applied to a subregion to produce a continuous estimation of <i>Q</i><sub>355</sub> along the river network extracted from a 90m elevation model. The results of the study show that Top-kriging and PSBI present complementary features. Top-kriging outperforms PSBI at larger river branches while PSBI outperforms Top-kriging for headwater catchments. Overall, they have comparable performances (Nash-Sutcliffe efficiencies in cross-validation of 0.89 and 0.83, respectively). Both techniques provide plausible and accurate predictions of <i>Q</i><sub>355</sub> in ungauged basins and represent promising opportunities for regionalization of low-flows