Rainfall Intensity-Duration-Frequency
modelling Level 4

Rainfall Intensity-Duration-Frequency

How intense will the 100-year storm be? Intensity-Duration-Frequency curves relate rainfall intensity to storm duration and return period via extreme value statistics, enabling hydrologic design of drainage systems and flood protection. This model derives IDF equations, implements extreme value analysis, demonstrates design storm calculations, and maps precipitation risk across return periods.

Prerequisites: extreme value statistics, gev distribution, idf relationships, return periods

Updated 13 min read

1. The Question

What rainfall intensity should this culvert be designed for?

Design storm problem:

Every hydraulic structure requires design precipitation.

Too small: Structure fails, flooding occurs
Too large: Unnecessarily expensive

IDF curves provide answer:

Relate three quantities:

  • Intensity (I): Rainfall rate (mm/hr or in/hr)
  • Duration (D): Storm length (minutes to days)
  • Frequency (F): Return period (years)

Return period:

Average recurrence interval between events of given magnitude.

100-year storm: 1% annual exceedance probability

Not: “Occurs once per century”
Rather: “1% chance each year”

Applications:

  • Storm sewer design (10-25 year typical)
  • Culvert sizing (25-50 year)
  • Dam spillway capacity (100-10,000 year)
  • Urban drainage systems
  • Flood insurance rate maps
  • Building codes (rainfall loads)

Key insight:

Short-duration storms are intense (cloudburst).

Long-duration storms are moderate (multi-day rain).

Inverse relationship captured by IDF curves.

2. The Conceptual Model

Extreme Value Theory

Annual maximum series:

For each year, extract maximum rainfall intensity for given duration.

Example - 1-hour duration:

  • 2010: 45 mm/hr
  • 2011: 32 mm/hr
  • 2012: 67 mm/hr (extreme event)
  • …
  • 2025: 41 mm/hr

Statistical distribution:

Fit probability distribution to these annual maxima.

Generalized Extreme Value (GEV) distribution:

\[F(x) = \exp\left[-\left(1 + \xi\frac{x-\mu}{\sigma}\right)^{-1/\xi}\right]\]

Where:

  • $\mu$ = location parameter (central tendency)
  • $\sigma$ = scale parameter (dispersion)
  • $\xi$ = shape parameter (tail behavior)

Shape parameter interpretation:

$\xi > 0$: Fréchet (heavy tail, no upper bound)
$\xi = 0$: Gumbel (exponential tail)
$\xi < 0$: Weibull (bounded upper tail)

Precipitation: Typically $\xi \approx 0$ (Gumbel) or slightly positive

Return period relationship:

\[T = \frac{1}{1 - F(x)}\]

100-year event: $F(x) = 0.99$

Design quantile:

For given return period $T$, solve for rainfall intensity:

\[x_T = \mu - \frac{\sigma}{\xi}\left[1 - \left(-\ln\left(1-\frac{1}{T}\right)\right)^{-\xi}\right]\]

Gumbel simplification ($\xi = 0$):

\[x_T = \mu - \sigma \ln\left(-\ln\left(1 - \frac{1}{T}\right)\right)\]

Sherman Equation

Empirical IDF relationship:

\[I = \frac{a}{(D + b)^c}\]

Where:

  • $I$ = intensity (mm/hr)
  • $D$ = duration (min)
  • $a, b, c$ = fitted parameters (vary with return period)

Parameter $a$: Scales with return period (increases for rarer events)

Parameter $c$: Controls duration dependence (typically 0.7-0.9)

Parameter $b$: Offset (typically 5-20 min)

Physical interpretation:

Short duration ($D$ small): Denominator small → $I$ large (intense)

Long duration ($D$ large): Denominator large → $I$ small (moderate)

Power law decay: Intensity decreases as duration increases

Typical values:

10-year, 1-hour: 50-80 mm/hr (moderate climate)
100-year, 1-hour: 80-120 mm/hr
10-year, 24-hour: 5-10 mm/hr

3. Building the Mathematical Model

Annual Maximum Series Analysis

Step 1: Data extraction

For each duration of interest (5 min, 10 min, 15 min, 30 min, 1 hr, 2 hr, 6 hr, 24 hr):

Extract annual maximum intensity from rainfall records.

Requires: Long-term precipitation data (30+ years ideal)

Step 2: Fit GEV distribution

Maximum likelihood estimation of parameters $\mu, \sigma, \xi$.

Log-likelihood function:

\[\ell(\mu,\sigma,\xi) = -n\ln\sigma - (1+1/\xi)\sum\ln\left[1+\xi\frac{x_i-\mu}{\sigma}\right] - \sum\left[1+\xi\frac{x_i-\mu}{\sigma}\right]^{-1/\xi}\]

Maximize $\ell$ numerically to find parameter estimates.

Step 3: Calculate quantiles

For desired return periods (2, 5, 10, 25, 50, 100, 500 years):

\[I_T = \mu - \frac{\sigma}{\xi}\left[1 - (-\ln(1-1/T))^{-\xi}\right]\]

Step 4: Repeat for all durations

Results in intensity matrix: $I(D, T)$

IDF Equation Fitting

General form:

\[I = \frac{K T^m}{(D+b)^n}\]

Where:

  • $K$ = regional scaling factor
  • $m$ = return period exponent (0.15-0.25)
  • $n$ = duration exponent (0.7-0.9)
  • $b$ = offset parameter

Fitting procedure:

Logarithmic transformation:

\[\ln I = \ln K + m \ln T - n \ln(D + b)\]

Multiple regression: Fit to observed $I(D,T)$ values

Alternatively: Nonlinear least squares on original form

Result: Single equation valid for all durations and return periods

Example - Denver, Colorado:

\[I = \frac{28.5 T^{0.17}}{(D + 8.5)^{0.78}}\]

Units: $I$ in inches/hr, $D$ in minutes

Confidence Intervals

Uncertainty sources:

  • Sampling variability (limited record length)
  • Distribution choice (GEV vs Gumbel vs Log-Pearson III)
  • Non-stationarity (climate change)

Confidence interval:

\[I_T \pm t_{\alpha/2} \times SE(I_T)\]

Where:

  • $t_{\alpha/2}$ = Student’s t value (95% → 1.96)
  • $SE$ = standard error from GEV fit

Typical: 95% CI width = ±20-40% for 100-year event

Increases with return period (extrapolating beyond data)

4. Worked Example by Hand

Problem: Design storm sewer system.

Location: Urban area, Midwest USA

Regional IDF equation:

\[I = \frac{1200 T^{0.2}}{(D+10)^{0.8}}\]

Units: $I$ (mm/hr), $D$ (min), $T$ (years)

Design criteria:

  • Return period: 10 years (urban standard)
  • Time of concentration: 15 minutes (computed from watershed)

Calculate:

  1. Design rainfall intensity
  2. Total rainfall depth
  3. Peak runoff (rational method)

Watershed properties:

  • Area: 5 hectares
  • Imperviousness: 70% (urban commercial)
  • Runoff coefficient: $C = 0.75$

Solution

Step 1: Design intensity

Apply IDF equation with $T = 10$ years, $D = 15$ min:

\[I = \frac{1200 \times 10^{0.2}}{(15+10)^{0.8}}\]

Calculate return period term:

$$10^{0.2} = 10^{1/5} = \sqrt[5]{10} = 1.585$$

Calculate duration term:

$$25^{0.8} = 25^{4/5} = (25^4)^{1/5} = (390625)^{1/5} = 14.62$$

Compute intensity:

\[I = \frac{1200 \times 1.585}{14.62} = \frac{1902}{14.62} = 130.1 \text{ mm/hr}\]

Design intensity: 130 mm/hr

Step 2: Rainfall depth

Convert intensity to total depth over duration:

\[P = I \times D = 130.1 \text{ mm/hr} \times \frac{15 \text{ min}}{60 \text{ min/hr}}\] \[P = 130.1 \times 0.25 = 32.5 \text{ mm}\]

Total rainfall: 32.5 mm in 15 minutes

Step 3: Peak runoff (Rational Method)

\[Q = C \times I \times A\]

Where:

  • $C = 0.75$ (runoff coefficient)
  • $I = 130.1$ mm/hr
  • $A = 5$ ha = 50,000 m²

Unit conversion:

\[Q = 0.75 \times \frac{130.1 \text{ mm/hr}}{1000 \text{ mm/m}} \times \frac{1 \text{ hr}}{3600 \text{ s}} \times 50000 \text{ m}^2\] \[Q = 0.75 \times 0.0361 \text{ m/s} \times 50000 \text{ m}^2\] \[Q = 1354 \text{ m}^3\text{/hr} = 0.376 \text{ m}^3\text{/s} = 376 \text{ L/s}\]

Peak runoff: 376 L/s

Step 4: Sewer sizing

Design sewer pipe to convey 376 L/s.

Manning equation (open channel flow):

\[Q = \frac{1}{n} A R^{2/3} S^{1/2}\]

Where:

  • $n = 0.013$ (concrete pipe)
  • $S = 0.005$ (0.5% slope, typical)
  • $A$ = cross-sectional area
  • $R$ = hydraulic radius

For circular pipe running 80% full (design standard):

Using Manning charts or iteration: 600 mm diameter pipe required

Step 5: Safety check

Verify capacity exceeds design flow with safety margin.

Actual capacity at 80% depth, 600 mm, $S = 0.005$:

\[Q_{capacity} \approx 420 \text{ L/s}\]

Margin: $420/376 = 1.12$ (12% safety margin)

Acceptable (>10% desired)

5. Computational Implementation

Below is an interactive IDF curve generator.

Intensity: -- mm/hr

Rainfall depth: -- mm

Exceedance prob: --%

Design standard: --

Observations:

  • Log-log plot shows power-law relationships
  • Curves decrease monotonically with duration
  • Higher return periods shift curves upward
  • Short durations show steepest slopes
  • Current design point highlighted in red
  • Regional differences significant (Southeast vs Southwest)

Key insights:

  • 5-minute cloudbursts extremely intense (150+ mm/hr)
  • 24-hour events moderate (5-15 mm/hr)
  • Doubling return period increases intensity ~15-20%
  • Regional climate strongly affects IDF parameters

6. Interpretation

Design Standards by Application

Storm sewers (urban):

Return period: 10-25 years

Rationale: Economic balance (cost vs risk)

Consequences of exceedance: Nuisance flooding, property damage

Example: Chicago uses 10-year for combined sewers

Highway drainage:

Return period: 25-50 years

Rationale: Public safety, traffic disruption

Interstate: 50-year
Secondary roads: 25-year

Culverts under roads: 50-100 year

Dam spillways:

Return period: 100-10,000 years

Rationale: Life safety, catastrophic failure potential

Small dams (low hazard): 100-year
Large dams (high hazard): 1000-10,000 year (PMF)

PMF (Probable Maximum Flood): Design standard for critical dams

Based on PMP (Probable Maximum Precipitation), not statistical

Building codes:

Return period: 50-year (roof drainage)

Green roofs: 25-year

Scupper sizing: 100-year

Regional Variations

Climate influences IDF curves:

Southeast USA:

  • High $K$ values (intense convection)
  • Frequent thunderstorms
  • 100-year, 1-hour: 100-130 mm/hr

Southwest USA (arid):

  • Lower $K$ (less moisture)
  • But extreme events still severe
  • Monsoon-driven
  • 100-year, 1-hour: 60-80 mm/hr

Pacific Northwest:

  • Frontal precipitation dominates
  • Less intense, more persistent
  • 100-year, 1-hour: 40-60 mm/hr
  • 100-year, 24-hour: 150-200 mm (high!)

Great Plains:

  • Severe convection
  • Hail-producing storms
  • Flash flood potential
  • 100-year, 1-hour: 90-120 mm/hr

Climate Change Impacts

Clausius-Clapeyron scaling:

\[\frac{dP}{dT} \approx 7\%/°C\]

Warmer atmosphere → more moisture → more intense precipitation

Observed trends:

  • Extreme events intensifying 5-10%
  • Return periods shifting (100-year → 50-year)
  • IDF curves becoming outdated

Design implications:

NOAA Atlas 14: Updated IDF data (2000s-2010s)

Replaces older TP-40 (1960s) and TP-49 (1970s)

Changes: 20-40% increases in some regions

Looking forward:

Climate-adjusted IDF: Incorporate future projections

Stationary assumption violated: Past ≠ future

Adaptive design: Build in flexibility, oversizing

7. What Could Go Wrong?

Short Record Length

Problem:

Estimating 100-year event from 30 years of data.

Extrapolation risk: Uncertainty increases exponentially

Confidence intervals: ±40% for 100-year from 30-year record

Example:

Estimated 100-year rainfall: 100 mm/hr

95% CI: 60-140 mm/hr

Solution:

  • Regional frequency analysis (pool data from similar stations)
  • Longer records when available (>50 years better)
  • Bayesian methods (incorporate prior information)

Non-Stationarity

Climate change violates stationarity assumption.

Past distribution ≠ future distribution

IDF curves shift upward over time

Example - Houston:

Atlas 14 (2018): 100-year, 24-hour = 16 inches
Previous TP-40 (1961): 100-year, 24-hour = 13 inches

23% increase in design rainfall

Consequences:

Existing infrastructure underdesigned.

Solution:

  • Update IDF curves regularly (10-20 year cycle)
  • Climate-adjusted design (add 20-30% safety margin)
  • Green infrastructure (adaptable)

Spatial Variability

IDF curves represent point estimates (rain gauge).

Storms are spatially variable:

Peak intensity may miss gauge.

Areal reduction factors:

\[I_{areal} = I_{point} \times ARF\]

Where ARF < 1, depends on area and duration.

Example: 50 km² watershed, 1-hour storm

ARF ≈ 0.85

Design intensity reduced 15% from point value.

Solution:

  • Apply ARF for large watersheds (>25 km²)
  • Radar-based IDF (spatial coverage)

Mixed Distributions

Two storm types:

Convective: Short, intense (thunderstorms)

Frontal: Long, moderate (multi-day rain)

Single IDF curve may not capture both.

Bimodal distribution: Two peaks

Example - Coastal California:

Short durations: Convective (summer)
Long durations: Atmospheric rivers (winter)

Solution: Separate IDF curves by season or storm type

8. Extension: NOAA Atlas 14

Comprehensive IDF resource for United States.

Coverage:

Released in volumes by region (2004-2018).

All 50 states now covered.

Features:

Online tool: Point-and-click map interface

Gridded data: 2-4 km resolution

Confidence intervals: Quantified uncertainty

Multiple durations: 5 min to 60 days

Multiple return periods: 1-year to 1000-year

Access: https://hdsc.nws.noaa.gov/hdsc/pfds/

Example query:

Select location → Generates IDF table and graphs

Comparison to older methods:

TP-40 (1960s): Outdated, regional maps, interpolation

Atlas 14: Modern, high-resolution, web-based

Updates: Increased design values 10-40% many areas

9. Math Refresher: Probability and Return Periods

Exceedance Probability

Annual Exceedance Probability (AEP):

\[P(X > x) = 1 - F(x)\]

Return period:

\[T = \frac{1}{AEP} = \frac{1}{1 - F(x)}\]

Relationship:

T (years) AEP Probability
2 0.50 50%
5 0.20 20%
10 0.10 10%
25 0.04 4%
50 0.02 2%
100 0.01 1%

Risk Over Design Life

Probability of exceedance during $n$ years:

\[P_{n} = 1 - (1 - AEP)^n\]

Example: 100-year event, 50-year design life

\[P_{50} = 1 - (1 - 0.01)^{50} = 1 - 0.605 = 0.395\]

39.5% chance of exceeding 100-year event in 50 years!

Often misunderstood: “100-year storm won’t happen again for 100 years”

Wrong! Each year independent, always 1% chance.

Summary

  • IDF curves relate rainfall intensity inversely to duration via empirical equations fitted to extreme value distributions
  • Generalized Extreme Value distribution models annual maximum precipitation with three parameters controlling location scale and shape
  • Sherman equation I = a/(D+b)^c provides power-law relationship between intensity and duration
  • Return period T = 1/AEP where annual exceedance probability defines design risk level
  • Typical design standards range from 10-year (urban drainage) to 10,000-year (high-hazard dams)
  • NOAA Atlas 14 provides gridded IDF data for USA with 20-40% increases over historical values
  • Climate change violating stationarity assumption requiring updated IDF curves and adaptive design
  • Regional variations significant with Southeast showing highest intensities and Pacific Northwest lowest short-duration values
  • Confidence intervals widen dramatically for rare events with ±40% typical for 100-year estimates from 30-year records
  • Critical tool enabling hydraulic infrastructure design from storm sewers to dam spillways

References