**Research article**
30 Jul 2021

**Research article** | 30 Jul 2021

# Impact of detention dams on the probability distribution of floods

Salvatore Manfreda Domenico Miglino and Cinzia Albertini

^{1},

^{1},

^{1,2}

**Salvatore Manfreda et al.**Salvatore Manfreda Domenico Miglino and Cinzia Albertini

^{1},

^{1},

^{1,2}

^{1}Dipartimento di Ingegneria Civile, Edile e Ambientale, Università degli Studi di Napoli Federico II, 80125 Naples, Italy^{2}Dipartimento di Scienze Agro Ambientali e Territoriali, Università degli Studi di Bari Aldo Moro, 70126 Bari, Italy

^{1}Dipartimento di Ingegneria Civile, Edile e Ambientale, Università degli Studi di Napoli Federico II, 80125 Naples, Italy^{2}Dipartimento di Scienze Agro Ambientali e Territoriali, Università degli Studi di Bari Aldo Moro, 70126 Bari, Italy

**Correspondence**: Salvatore Manfreda (salvatore.manfreda@unina.it)

**Correspondence**: Salvatore Manfreda (salvatore.manfreda@unina.it)

Received: 25 Jan 2021 – Discussion started: 22 Feb 2021 – Revised: 07 Jun 2021 – Accepted: 08 Jun 2021 – Published: 30 Jul 2021

Detention dams are one of the most effective practices for flood mitigation. Therefore, the impact of these structures on the basin hydrological response is critical for flood management and the design of flood control structures. With the aim of providing a mathematical framework to interpret the effect of flow control systems on river basin dynamics, the functional relationship between inflows and outflows is investigated and derived in a closed form. This allowed the definition of a theoretically derived probability distribution of the peak outflows from in-line detention basins. The model has been derived assuming a rectangular hydrograph shape with a fixed duration and a random flood peak. In the present study, the undisturbed flood peaks are assumed to be Gumbel distributed, but the proposed mathematical formulation can be extended to any other flood-peak probability distribution. A sensitivity analysis of parameters highlighted the influence of detention basin capacity and rainfall event duration on flood mitigation on the probability distribution of the peak outflows. The mathematical framework has been tested using for comparison a Monte Carlo simulation where most of the simplified assumptions used to describe the dam behaviours are removed. This allowed demonstrating that the proposed formulation is reliable for small river basins characterized by an impulsive response. The new approach for the quantification of flood peaks in river basins characterized by the presence of artificial detention basins can be used to improve existing flood mitigation practices and support the design of flood control systems and flood risk analyses.

During the last decades, the growing number of hydrological extremes have raised economic losses and risk perception at the global scale (Peduzzi, 2005; Di Baldassarre et al., 2010; Winsemius et al., 2016). The impact of natural disasters has been quantified in a recent study by Wallemacq and Below (2015), which is based on the Emergency Events Database (EM-DAT). According to EM-DAT, flooding impacted on nearly 2.5 billion people in the period 1994–2013. A more recent study by Munich Re (NatCatSERVICE, Munich Reinsurance Company (MunichRe), 2020) reported 3798 flash flood events that produced economic losses of about USD 592 billion and killed around 100 000 people worldwide during the last two decades (2000–2018).

In this context, climate change and anthropic activities are probably accelerating the number of extremes (Fischer and Knutti, 2016; Papalexiou and Montanari, 2019). In fact, these two factors are significantly modifying river basin hydrology (Di Baldassarre et al., 2017), which now requires renewed methods and techniques. Therefore, the need to find appropriate measures for flood hazard prevention and mitigation is becoming more and more pressing.

Flood risk management is a complex challenge for hydrologists that need to identify dynamic solutions for flood mitigation. Plate (2002) underlined the need to update flood mitigation plans according to the changing conditions (e.g. climate, populations, land use). Hence, the design of flood mitigation systems requires rapid methodologies to evaluate scenarios and intervention measures and obtain effective flood risk management strategies (Franzi et al., 2016).

One of the most effective actions for the attenuation of peak discharges is represented by a detention basin system or coordinated reservoir operations (e.g. UDFCD, 2016; Jacob et al., 2019; Seibert et al., 2014). The proper design of a detention basin is an extremely complex task given the large number of factors that are involved in the decision process: technical engineering features of the structure, site selection and environmental characteristics. For this reason, Bellu et al. (2016) proposed an innovative method to optimize the dimensioning and site selection of a flood mitigation system. The method follows three steps that include a preliminary sizing, a site location and optimization according to environmental objectives. Nevertheless, the sizing of the detention volume is based, in most of the cases, on a given design flood event neglecting the random nature of the phenomena and the impact of the structure on a large spectrum of events.

The flood attenuation caused by the presence of artificial reservoirs is influenced by hydrologic and hydraulic factors, such as flood wave shape and duration, as well as by the storage capacity and geometric parameters of the detention basin. Therefore, it is crucial to build a mathematical scheme able to interpret the functional relationships, even in a simplified form, among the mentioned variables.

With this aim, the present article introduces a theoretically derived probability distribution (TDD) of detention basin outflows, which is obtained assuming the incoming flood peaks randomly distributed and characterized by rectangular hydrographs of fixed duration (see e.g. Manfreda et al., 2018). In the present case, flood distribution is assumed to be Gumbel for the sake of simplicity, but the proposed approach can be applied to any flood-peak probability distribution. Performances of the proposed method are investigated under different configurations and are tested with a numerical simulation of flood mitigation. This method can be applied to improve existing flood mitigation approaches and support the design of flood control systems, flood risk and damage analyses.

This paper is organized as follows: in Sect. 2, the conceptual scheme of the detention basin adopted to derive the mathematical formulation of outflows is described. Furthermore, the hydraulic concepts and assumptions to analytically compute the derived probability distribution of the peak outflows are introduced; in Sect. 3, the proposed methodology is tested under different reservoir configurations and compared with the results of the numerical simulations, and a brief description of the numerical model is also provided; in Sect. 4 main findings and results are discussed.

In-line detention dams (also known as flow-through dams) are constructed solely with the purpose of flood control and mitigation of flood risks in downstream communities and ecosystems. Unlike reservoir dams, which are primarily built for water storage or power generation, the spillway (opening) is located at the same height as the riverbed level, allowing the river to continue its natural flow under normal conditions. When water levels rise above the spillway, the dam restricts the amount flowing through the opening, decreasing peak flow. Since detention dams minimally affect rivers' natural flows, under normal conditions negative environmental and socioeconomic impacts, such as sediment accumulation, restriction of water flow to downstream communities and ecosystems, and breaching during very extreme flood events, can be minimized or avoided altogether.

The schematization of the detention dam has been simplified with the aim to obtain a mathematical description of the outflows associated with a given hydrograph. In particular, we assumed that the dam body has two openings: a low-level opening at the basement and a crest spillway. The former is assumed to let pass the flow below a given control value, while the latter starts functioning only when the volume of the dam is completely filled up to the crest level. Thereafter, the water starts to flow from both the openings, and the control is mainly exerted by the basin volume above the emergency spillway (crest level). This scheme can be described in closed form that may help the construction of a derived probability distribution of the outflow from a detention dam.

## 2.1 The hydraulic characteristics of the problem

Let us first introduce the key equations controlling the dynamics of a system like the one under study. The first equation to introduce is represented by the stage-storage capacity curve which is able to describe the morphology of the gorge closed by a specific dam. The function is generally represented by a power law:

where *W*(*h*) [m^{3}] is the water storage of detention basin,
*w*_{1} [m^{3−n}] is the parameter of the stage-storage capacity
curve, *h* [m] is the water level in the reservoir and *n* [–] is the exponent
influenced by the shape of the control volume. The exponent ranges between 1
and 4.5, where 1 is associated with a prismatic geometry with vertical
surrounding walls and 4.5 is associated with a more complex morphology
closed by more gentle lateral slopes.

The streamflow of the river system reaching the reservoir is altered by it based on the hydraulic characteristics of the dam. The outflow will be controlled by the amount of water accumulated in the system according to the continuity equation:

where *Q*_{in}(*t*) is the incoming flux and
*Q*_{out}(*t*) is the outflow from the reservoir.

The outflow can be computed using the traditional formulation of hydraulics
based on the variation of the water level, *h*, in the reservoir, which can be
derived from the stage-storage capacity curve (Eq. 1) and continuity
equation (Eq. 2). In particular, we can assume that the outflow for a
simple scheme with two openings (the low-level opening and the emergency
spillway) and the characteristics described in Fig. 1 may vary according
to four different discharge laws expressed as a function of the reservoir
water level stage. The outflow can be described as follows:

where *μ*_{f} [–] is the coefficient of discharge of the low-level
opening at the basement (suggested values may range between 0.6 (assuming a
thin sharp edge) and 0.8 (assuming a wall with thickness 2 times larger
than the opening height)), *A* [m^{2}] is the area of the low-level opening
(i.e. the product of the opening width, *b* [m], and the minimum dimension
between the top of the flow surface at the opening exit and the bottom of
the opening), *μ*_{s} [–] is the coefficient of discharge of the spillway
crest (suggested values range around 0.3–0.4 based on the geometry of the
weir), *l*_{f} [*L*] is the height of the low-level opening at the basement,
*L* [m] is the effected crest length, *h*_{f} [m] is the height of the
barycentre of the low-level opening, *h*_{s} [m] is the height of the
spillway crest, and *g* [m/s^{2}] is the acceleration due to gravity.

The above expression includes the three main configurations that may occur
in the proposed scheme with the increase of the water level stage in the
reservoir. A detailed description of the scheme proposed along with
graphical indication of the parameters meaning is given in Fig. 1. In
particular, the opening at the bottom does not exert a significant control
on the incoming flow as long as the water level does not generate a
submergence of the opening (i.e. as long as *h*≤*h*_{s}). The water flow
starts to be limited when the water stage reaches the value *h*_{f}, after
which the opening is submerged. This allows filling the storage volume of
the reservoir up to the level of the crest spillway. After this stage, the
reservoir tends to operate a mitigation that is influenced by the water
storage capacity of the reservoir above the crest level and the hydraulic
characteristics of the spillway. An example of flood mitigation obtained via
numerical simulation is given in Fig. 2, which provides a comparison
between a synthetic hydrograph and the outflow from a detention dam.

For the scope of the present study, the outflow has been simplified through the following set of equations:

where *Q*_{c} is the control value of discharge that is computed using the
discharge equation of the submerged opening and assuming *h*=*h*_{s}.

Assuming a rectangular hydrograph of the incoming flow, it is possible to derive the peak flow associated with an incoming flood peak. Following the simplifying assumption given in Eq. (4), the outflow is not affected by the presence of the dam for lower streamflow values, while it is modified when the inflow exceeds the control discharge. In particular, the outflow remains almost constant as long as the reservoir is filled and, thereafter, it is controlled by the crest spillway. In this last configuration, it is possible to use the linear reservoir concept for the water volume accumulated above the elevation of the crest spillway.

In order to estimate the peak flow associated with a specific rectangular
hydrograph of constant discharge equal to *Q*_{max}, we should recall the
expression of the peak flow generated by a simple linear reservoir, which
can be described as follows:

where *t*_{p} [s] is the event duration, and *k* [s] is the delay constant
of the conceptual linear reservoir.

The above equation should be modified, in the present case, considering that the flood event should fill the dam water storage capacity (or detention basin) before reaching the crest spillway. Therefore, the crest spillway will be activated only after a time:

where *W*_{max} is the volume of water accumulated in the dam at the crest
level *h*_{s}.

When the volume below the spillway crest is totally filled, the crest
spillway starts functioning for discharge values above the control discharge
that is released at the bottom. Therefore, the peak outflow, *Q*_{p,out},
assumes the following form:

where *k*_{eq} [s] is the equivalent delay constant of the conceptual
reservoir associated with the outflow. This parameter can be derived
exploiting the characteristics of the spillway and the stage-storage
capacity curve. In particular, according to the linearity concept, the two
functions should have the same exponent. Under such a hypothesis, the
parameter *k*_{eq} can be estimated as

where *w*_{2} is the coefficient of the rescaled stage-storage capacity
curve above the crest level. Such a parameter should be computed in order to
get the best approximation of the function describing the volumes above the
mentioned level, *h*_{s} imposing a coefficient *n* equal to 1.5 for the
rescaled stage-storage capacity curve (referred to the stage-storage
capacity curve above the crest level). With this aim, the parameter can be
computed by comparing the two functions and imposing that they are equal in
a point *h*_{m} which is representative of the range of variability of the
water level above the crest level.

Within the present article, we assumed the parameter *h*_{m} was set
equal to *h*_{s}.

These assumptions allowed deriving a functional relationship between the reservoir inflow and outflow, which can be used to invert the function with respect to the incoming flow and associate a probability with each flow discharged value by exploiting the theory of derived distributions (Benjamin and Cornell, 2014).

### Estimation of the event duration

The assumption of a rectangular hydrograph may produce a significant
overestimation of the flood volume. Therefore, the parameter *t*_{p} should
be defined accounting for the real volume associated with a realistic flood
hydrograph. In this context, we can rely on the flow duration frequency
reduction curve (FDF) proposed by the NERC (Natural Environment Research
Council, 1975), which describes the maximum average discharge
*q*(*D*) as a function of the event duration *D*:

where *ω* represents the characterizing basin time response that is
frequently associated with the lag time of the river basin.

Adopting the above formulation, Fiorentino (1985) suggested a simple form of synthetic hydrograph redistributing the volume symmetrically with respect to the time of the peak. This leads to the following form of hydrograph:

Based on the above formulation, it is possible to impose that the duration
*t*_{p} of the equivalent rectangular event has the same volume of the
synthetic hydrograph of Eq. (10) in the temporal window of *ω*
around the peak flow. This led to the following equivalence:

Given the above assumption, the term *t*_{p} will be named equivalent event
duration from now on.

## 2.2 Functional relationship between the incoming discharge and the outflow

The possibility to identify the analytical relationship between two processes where one represents the stochastic forcing allows determining the derived probability distribution of the variable. This approach has been used several times for flood maxima (Eagleson, 1972; De Michele and Salvadori, 2002; Gioia et al., 2008), soil moisture (Rodríguez-Iturbe and Porporato, 2007; Manfreda and Fiorentino, 2008), and scour process (Manfreda et al., 2018). In the present case, the methodology has been applied to the laminated flood peak.

With this aim, it is critical to identify the inverse function describing the laminated flood peaks as a function of the maximum floods reaching the detention dam. Equation (6) cannot be inverted in a closed form. Therefore, the function has been approximated using a Taylor expansion of the first order about the point $\mathrm{2}{Q}_{\mathrm{c}}+\frac{{W}_{\mathrm{max}}}{{t}_{\mathrm{p}}}$, where ${Q}_{\mathrm{c}}+\frac{{W}_{\mathrm{max}}}{{t}_{\mathrm{p}}}$ represents the incoming discharge that leads to overspilling. This approximation allows defining a closed form of the inverse function, which represents a good approximation for theoretical derivation.

Given the above approximations, it is possible to derive the inverse function of the peak discharge function of the peak outflow. This equation can be obtained mathematically by inverting Eq. (6) and exploiting the parametrization introduced above. The inverse function assumes the following form:

This expression can be used to analytically compute the derived probability
distribution of the peak outflow from a detention dam characterized by a
storage capacity *W*_{max}, an equivalent delay constant *k*_{eq}, invested
by flood hydrograph of equivalent event duration *t*_{p}. With this aim, any
probability distribution of the flood peaks can be used given the monotonic
nature of the above expression. The expression of the TDD will be (see
Benjamin and Cornell, 2014)

where the derivative of ${g}^{-\mathrm{1}}\left({Q}_{\mathrm{p},\mathrm{out}}\right)$ assumes the following form:

In order to describe the probability distribution of the outflows, we should divide it according to the three potential configurations of the detention dams: (1) undisturbed flow, (2) accumulation of water in the reservoir, and (3) activation of the crest spillway. Based on these assumptions, the probability distribution of the outflow is subdivided into three components and modelled by the following set of equations:

Equation (15) describes the general form of the probability distribution of
the outflow from a detention dam, where the first component coincides with
the distribution of the incoming flow as long as it is below the control
discharge of the lower opening (${Q}_{\mathrm{p},\mathrm{out}}<{Q}_{\mathrm{c}}$). Assuming that the lower
opening is able to control the outflow around *Q*_{c} after submergence,
there is a mass probability in *Q*_{c} depending on the storage volume of
the reservoir (${Q}_{\mathrm{p},\mathrm{out}}={Q}_{\mathrm{c}}$). After these two phases, the outflow is
affected by the lamination due to the water volume accumulation above the
crest level (${Q}_{\mathrm{p},\mathrm{out}}>{Q}_{\mathrm{c}}$).

## 3.1 Examples of application of the TDD of the detention basin outflows under different configurations

In order to explore the behaviour of the proposed formulation, we investigated the effects of different types of parametrization on the derived distribution starting from a single distribution of floods. In Fig. 3, the influence of the storage capacity and hydrograph duration on the outflow of the reservoir is shown. We depicted the probability density functions (PDFs) associated with increasing storage capacity obtained by raising the crest level from 2 m up to 8 m (moving top–down in the figure) and considering two distinct equivalent event durations of 30 min and 1 h. As expected, hydrographs with longer duration tend to saturate sooner the water storage capacity of the reservoir, reducing also the lamination effects. On the other hand, the increase in water storage capacity leads to a proportional growth of flood peak mitigation. These graphs describe the behaviour of a reservoir providing an output consistent with the dynamics of the process.

In Fig. 4, we modified the coefficient of the stage-storage capacity
curve, *w*_{1}, using the values of 5000 and 10 000 and explored equivalent
event durations ranging from half an hour to 2 h. Graphs display how
the PDFs of outflows are altered by the presence of a dam with these
characteristics. In the present example, the impact of the lower opening can
be better appreciated with a mass probability around *Q*_{c} that is equal
to 52 m^{3}/s. It must be clarified that the cross-section of the opening
has been increased in this second example on purpose to emphasize its impact
on the proposed mathematical scheme.

## 3.2 Testing the reliability of the proposed method

In the current work, a numerical simulation of flood mitigation through a detention basin with the characteristics reported in Fig. 1 was carried out using the same forcing adopted for the proposed TDD. Therefore, we adopted the Gumbel distribution as the reference distribution to generate random values of discharge and numerically simulate the dynamics of the detention dam and its peak outflows. It must be clarified that any probability distribution of floods (e.g. generalized extreme value, three-parameter log-normal, generalized logistic and Gumbel distributions) can be applied. An example of the numerical simulation is given in Fig. 2 with the consequent attenuation of the hydrograph due to the detention dam.

The numerical simulation was carried out with the main scope of testing the theoretically derived probability distribution of laminated peak flows and also quantifying the impact of the approximations adopted to obtain a closed form of the solution. Therefore, the comparison of the theoretically derived distribution and numerical outflows helps understand the reliability of the proposed methodology. Results are given in the following graphs.

Figure 5 provides a comparison of different PDFs obtained, modifying the
maximum water storage capacity of the dam and its height. It can be noted
that the values of the probability distribution replicate fairly well those
obtained with the numerical simulations. The approximation of a fixed
discharge from the submerged opening induces a small dispersion of values
around the control value of discharge, *Q*_{c}, that the theoretical
probability distribution is not able to capture. Moreover, the adopted
approximations slightly overestimate the outflows which can be interpreted
as a safety approximation for flood mitigation planning.

Comparing the different PDFs, it should be clarified that the parameters
have been changed, looking for combinations of dam heights, *h*_{s}, and
coefficient *w*_{1} leading to similar water storage capacity on
each row. This allows demonstrating that it is much more effective to
increase the area flooded by the reservoir (the parameter *w*_{1}
represents the rate of increase of the water storage with the water level)
rather than increase the height of the dam.

The last analysis performed is given in Fig. 6, where we tested the performances of the TDD assuming the lower opening closed. Such a condition may be representative of an ordinary dam with an assigned flood retention volume or used for water supply purposes, where the volume above the crest spillway leads to a lamination of floods. This configuration has been compared with the scheme described in the previous sections that includes the presence of the lower opening. The difference between the two configurations is given in Fig. 6a (closed opening) and b (lower opening active). The two graphs show the ability of the mathematical formulation to properly interpret also the present configuration, offering a wide spectrum of potential applications in hydraulic design.

Finally, in order to test the impact of the assumption of rectangular inflow hydrograph, we also compared the results of a numerical simulation where the hydrograph is assumed to be a symmetric exponential one according to the expression given in Eq. (11). The comparison is given in Fig. 7 for three different equivalent event durations. It can be appreciated how the use of an equivalent event duration produced probability distributions of the outflow that look very similar to those obtained with a symmetric exponential hydrograph. It must also be underlined that with the increase of the duration of the event such an approximation tends to deteriorate the result of the proposed model.

The present article introduces a new formulation useful to quantify the impact of detention dams on the probability distribution of floods. We must acknowledge that the formulation was obtained with several simplifying assumptions that include the shape of the incoming hydrograph, the approximation used to interpret the flow through the lower opening, and the approximation of a linear reservoir for the flow above the crest level. Summing all these, it is really satisfying to see that the obtained formulation can interpret the dynamics of such hydraulic infrastructures fairly well, providing an analytical description of the impact of artificial reservoirs on flood dynamics. This may be extremely useful in properly addressing the effects of water infrastructures on floods. The TDD can be used for detention dams, but the formalism can also be applied to ordinary dams just setting the control discharge to zero and assigning a given value of the water level in the dam. Therefore, the formalism is versatile and can be applied in different contexts. The strongest assumption is represented by the rectangular hydrograph which can be realistic for small river basins characterized by relatively short concentration time. With the aim of minimizing the impact of such a choice, we adopted an equivalent event duration in the formulation that allowed accounting for the flood variability during a specific event. However, this assumption may become limitative in large river basins where a hydrograph evolves over large areas and its shape is also not simple to be predicted. In fact, large river basins may display complex hydrographs with multiple peaks that require a specific approach.

The proposed method may be used in some contexts such as the projects and design of small lamination dams and detention dams in small river basins. For instance, there are several river basins along the coastline that drain a high amount of water in short durations, affecting cities and towns developed along the waterfront. These areas are typically exposed to frequent flood events that may impose the need to properly identify potential solutions for flood mitigation.

The scheme can be used to carry out preliminary dimensioning of these structures and eventually could be coupled with other tools to identify optimal configurations for flood mitigation. Moreover, the scheme can be applied to any probability distribution of floods, including the case of floods that are already subject to laminations allowing the description of a scheme of nested dams.

This topic is still under investigation, and its study will be applied to identify optimal solutions in flood control systems quantifying the impact of structure on the full spectra of floods.

α [–] |
Scale parameter of Gumbel distribution |

β [–] |
Location parameter of Gumbel distribution |

μ_{f} [–] |
Coefficient of discharge of the submerged low-level opening |

μ_{s} [–] |
Coefficient of discharge of the crest spillway |

A [m^{2}] |
Area of the low-level opening |

b [m] |
Width of the low-level opening rectangular section |

d [m] |
Height of the low-level opening |

g [m/s^{2}] |
Acceleration due to gravity |

h [m] |
Variable water level within the detention basin |

h_{m} [m] |
Mean of water levels over the spillway crest |

h_{s} [m] |
Height of the spillway crest |

h_{f} [m] |
Height of the barycentre of the low-level opening |

k [s] |
Storage coefficient of the linear reservoir method |

k_{eq} [s] |
Equivalent delay constant of the conceptual reservoir associated with the outflow |

L [m] |
Effected crest length |

n [–] |
Exponent of the stage-storage capacity curve |

p(Q) [–] |
Probability density function of outflows |

l_{f} [m] |
Height of the low-level opening |

Q_{c} [m^{3}/s] |
Design outflow from the low-level opening |

Q_{in} [m^{3}/s] |
Inflow in the detention basin |

Q_{max} [m^{3}/s] |
Peak flow incoming in the detention basin |

Q_{out} [m^{3}/s] |
Outflow from the detention basin |

Q_{p,out} [m^{3}/s] |
Peak outflow from the detention basin |

t [s] |
Time |

D [s] |
Flood event duration |

t_{filling} [s] |
Time after which the crest spillway starts functioning |

t_{p} [s] |
Equivalent flood duration |

ω [s] |
Lag time of the river basin |

W_{max} [m^{3}] |
Water storage capacity at the crest level |

W [m^{3}] |
Variable storage capacity of the detention basin |

w_{1} [m^{3−n}] |
Parameter of the stage-storage capacity curve |

w_{2} [m${}^{\mathrm{3}/\mathrm{2}}$] |
Parameter of the equivalent stage-storage capacity curve |

Codes used to generate the theoretical distributions introduced in the present article are available on MATLAB Central at the following link https://it.mathworks.com/matlabcentral/fileexchange/95813-peak-outflows-of-a-detention-basin (Manfreda, 2021).

SM conceived and coordinated the work, DM supported the coding and numerical simulations, and CA supported the review and interpretation of the results. All authors contributed to the writing and review of the article.

The authors declare that they have no conflict of interest.

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This research has been supported by the Italian Ministry of Environment, Land and Sea with the research project entitled “Hydraulic risk mitigation in coastal basins with in-line expansion tanks: an integrated sizing approach”.

This paper was edited by Giuliano Di Baldassarre and reviewed by two anonymous referees.

Bellu, A., Fernandes, L. F. S., Cortes, R. M., and Pacheco, F. A.: A framework model for the dimensioning and allocation of a detention basin system: The case of a flood-prone mountainous watershed, J. Hydrol., 533, 567–580, 2016.

Benjamin, J. R. and Cornell, C. A.: Probability, statistics, and decision for civil engineers, Courier Corporation, Mineola, New York, 2014.

De Michele, C. A. R. L. O. and Salvadori, G.: On the derived flood frequency distribution: analytical formulation and the influence of antecedent soil moisture condition, J. Hydrol., 262, 245–258, 2002.

Di Baldassarre, G., Montanari, A., Lins, H., Koutsoyiannis, D., Brandimarte, L., and Bloeschl, G.: Flood fatalities in Africa: from diagnosis to mitigation, Geophys. Res. Lett., 37, L22402, https://doi.org/10.1029/2010GL045467, 2010.

Di Baldassarre, G., Martinez, F., Kalantari, Z., and Viglione, A.: Drought and flood in the Anthropocene: feedback mechanisms in reservoir operation, Earth Syst. Dynam., 8, 225–233, https://doi.org/10.5194/esd-8-225-2017, 2017.

Eagleson P. S.: Dynamics of flood frequency, Water Resour Res., 8, 878–98, 1972.

Fiorentino, M.: La valutazione dei volumi dei volumi di piena nelle reti di drenaggio urbano, Idrotecnica, 3, 141–152, 1985.

Fischer, E. M. and Knutti, R. Observed heavy precipitation increase confirms theory and early models, Nat. Clim. Change, 5, 560–564, 2016.

Franzi, L., Pezzoli, A., and Besana, A.: Flood Lamination Strategies for Risk Reduction, River Basin Management, 113, 315, https://doi.org/10.5772/63553, 2016.

Gioia, A., Iacobellis, V., Manfreda, S., and Fiorentino, M.: Runoff thresholds in derived flood frequency distributions, Hydrol. Earth Syst. Sci., 12, 1295–1307, https://doi.org/10.5194/hess-12-1295-2008, 2008.

Jacob, A. C. P., Rezende, O. M., de Sousa, M. M., de França Ribeiro, L. B., de Oliveira, A. K. B., Arrais, C. M., and Miguez, M. G.: Use of detention basin for flood mitigation and urban requalification in Mesquita, Brazil, Water Sci. Technol., 79, 2135–2144, 2019.

Manfreda, S.: Peak outflows of a detention basin, MATLAB Central File Exchange [model code], available at: https://www.mathworks.com/matlabcentral/fileexchange/95813-peak-outflows-of-a-detention-basin, last access: 27 July 2021.

Manfreda, S. and Fiorentino, M.: A stochastic approach for the description of the water balance dynamics in a river basin, Hydrol. Earth Syst. Sci., 12, 1189–1200, https://doi.org/10.5194/hess-12-1189-2008, 2008.

Manfreda S., Link, O., and Pizarro, A.: The Theoretically Derived Probability Distribution of Scour, Water, 10, 1520, https://doi.org/10.3390/w10111520, 2018.

Munich Reinsurance Company [MunichRe]: Geo Risks Research, NatCatSERVICE. © 2020 Münchener Rückversicherungs-Gesellschaft, NatCatSERVICE, available at: https://www.munichre.com/en/solutions/for-industry-clients/natcatservice.html, last access: 26 May 2020.

Natural Environmental Research Council (NERC): Estimation of Flood Volumes over Different Duration, in: Flood Studies Report; Volume I, NERC: London, UK, pp. 352–373, 1975.

Papalexiou, S. M. and Montanari, A.: Global and regional increase of precipitation extremes under global warming, Water Resour. Res., 55, 4901–4914, 2019.

Peduzzi, P.: Is climate change increasing the frequency of hazardous events?, Environment and Poverty Times, 3, p. 7, 2005.

Plate, E. J.: Flood risk and flood management, J. Hydrol., 267, 2–11, 2002.

Rodríguez-Iturbe, I. and Porporato, A.: Ecohydrology of water-controlled ecosystems: soil moisture and plant dynamics, Cambridge University Press, Cambridge, USA, 2007.

Seibert, S. P., Skublics, D., and Ehret, U.: The potential of coordinated reservoir operation for flood mitigation in large basins – A case study on the Bavarian Danube using coupled hydrological – hydrodynamic models, J. Hydrol., 517, 1128–1144, 2014.

Urban Drainage and Flood Control District (UDFCD): Urban Storm Drainage Criteria Manual, Volume 2, Denver, Colorado, available at: https://udfcd.org/wp-content/uploads/uploads/vol2 criteria manual/USDCM Volume 2.pdf (last access: 15 June 2020), 2016.

Wallemacq, P. and Below, R. : The human cost of natural disasters: A global perspective, Centre for Research on the Epidemiology of Disasters, Brussels, Belgium, 2015.

Winsemius, H. C., Aerts, J. C., Van Beek, L. P., Bierkens, M. F., Bouwman, A., Jongman, B., Kwadijk, J. C., Ligtvoet, W., Lucas, P. L., van Vuuren, D. P. and Ward, P. J.: Global drivers of future river flood risk, Nat. Clim. Change, 6, 381–385, 2016.