IDF curve

An intensity-duration-frequency (IDF) curve is a mathematical function that relates the rainfall intensity with its duration and frequency of occurrence. These curves are commonly used in hydrology for flood forecasting and civil engineering for urban drainage design. However, the IDF curves are also analysed in hydrometeorology because of the interest in the time concentration or time-structure of the rainfall.

Mathematical approaches

The IDF curves can take different mathematical expressions, theoretical or empirically fitted to observed rainfall data. For each duration (e.g. 5, 10, 60, 120, 180 ... minutes), the ECDF or empirical distribution function, and a determined frequency or return period is set. Therefore, the empirical IDF curve is given by the union of the points of equal frequency of occurrence and different duration and intensity Likewise, a theoretical or semi-empirical IDF curve is one whose mathematical expression is physically justified, but presents parameters that must be estimated by empirical fits.

Empirical approaches

There is a large number of empirical approaches that relate the intensity (I), the duration (t) and the return period (p), from fits to power laws such as:

• Sherman's formula, with three parameters (a, c and n), which are a function of the return period, p:

[math]\displaystyle{ I(t)=\frac a {(t+c)^n} }[/math]

• Chow's formula, also with three parameters (a, c and n), for a particular return period p:

[math]\displaystyle{ I(t)= \frac a {t^n+c} }[/math]

• Power law according to Aparicio (1997), with four parameters (a, c, m and n), already adjusted for all return periods of interest:

[math]\displaystyle{ I(t,p)=a \cdot \frac{p^m}{(t+c)^n} }[/math]

In hydrometeorology, the simple power law (taking [math]\displaystyle{ \ c = 0 }[/math]) is used according to Monjo (2016) as a measure of the time-structure of the rainfall:

[math]\displaystyle{ I(t)=\frac a {t^n} = I_o\left( \frac{t_o} t \right)^n }[/math]

where [math]\displaystyle{ \ I_o }[/math] is defined as a intensity of reference for a fixed time [math]\displaystyle{ \ t_o }[/math], i.e. [math]\displaystyle{ \ a=I_o t_o^n }[/math], and [math]\displaystyle{ \ n }[/math] is a non-dimensional parameter known as n-index. In a rainfall event, the equivalent to the IDF curve is called Maximum Averaged Intensity (MAI) curve.

Theoretical approaches

To get an IDF curves from a probability distribution, [math]\displaystyle{ \ F(x) }[/math] it is necessary to mathematically isolate the precipitation[math]\displaystyle{ \ x }[/math], which is directly related to the average intensity [math]\displaystyle{ \ I }[/math] and the duration [math]\displaystyle{ \ t }[/math], by the equation [math]\displaystyle{ \ x = It }[/math], and since the return period [math]\displaystyle{ p }[/math] is defined as the inverse of [math]\displaystyle{ \ 1 - F(x) }[/math], the function [math]\displaystyle{ \ f(p) }[/math] is found as the inverse of [math]\displaystyle{ \ F(x) }[/math], according to:

[math]\displaystyle{ I t = f(p) \quad \Leftarrow \quad p = \frac{1}{1-F(I t)} }[/math]

• Power law with the return period, derived from the Pareto distribution, for a fixed duration [math]\displaystyle{ \ t }[/math]:

[math]\displaystyle{ \ I(p) = kp^m \quad \Leftarrow \quad F(It) = 1 - \left( \frac{kt}{It} \right)^{1/m} = 1 - \frac{1}{p} }[/math]

where the Pareto distribution constant has been redefined as[math]\displaystyle{ \ k' = k t }[/math], since it is a valid distribution for a specific duration of precipitation, [math]\displaystyle{ \ x }[/math], it has been taken as[math]\displaystyle{ \ x = It }[/math].

• Function derived from the generalized Pareto Distribution, for a given duration [math]\displaystyle{ \ t }[/math]:

[math]\displaystyle{ I(p) = \begin{cases} \mu + \frac \sigma m \cdot (p^m-1) \quad \Leftarrow \quad F(I) = 1 - \left(1+ \frac{m(I-\mu)}{\sigma}\right)^{-1/m} = 1 - \frac{1}{p} & \text{if } m \gt 0, \\ \quad \mu + \sigma\ln(p) \quad \quad \Leftarrow \quad F(I) = 1 - \exp \left( - \frac{I-\mu}{\sigma}\right) = 1 - \frac{1}{p} & \text{if } m = 0. \end{cases} }[/math]

Note that for [math]\displaystyle{ \ m \gt 0 }[/math] y [math]\displaystyle{ \ \mu = \frac \sigma m }[/math], the Generalized Pareto Distribution retrieves the simple form of the Pareto distribution, with [math]\displaystyle{ \ k' = \frac \sigma m }[/math]. However, with [math]\displaystyle{ \ m = 0 }[/math] the exponential distribution is retrieved.

• Function deduced from the Gumbel distribution and the opposite Gumbel distribution, for a given duration [math]\displaystyle{ \ t }[/math]:

[math]\displaystyle{ I(p) = \mu + \sigma\ln \left( \ln \left( 1 - \frac{1}{p} \right) \right) \quad \Leftarrow \quad \quad F(I) = \exp \left( - \exp \left( - \frac{I-\mu} \sigma \right) \right) = 1 - \frac{1}{p} }[/math]

[math]\displaystyle{ I(p) = \mu + \sigma\ln(\ln p) \quad \quad \quad \quad \quad \Leftarrow \quad \quad F(I) = 1 - \exp \left( - \exp \left( \frac{I-\mu}{\sigma} \right) \right) = 1 - \frac{1}{p} }[/math]