Skip to main content

Gumbel flood, confidence limits, hydrologic risk

Short answer

From Gumbel fit using given floods, mean=4358.4 m³/s, σ=1516.8 m³/s. 200-yr flood: 10849 m³/s, 1000-yr: 13046 m³/s. 95% confidence limits: 200-yr 8044–13653, 1000-yr 9424–16668 m³/s. Risk for 12300 m³/s over 25 yr: 4.2%. Design flow for 99% reliability over 25 yr: 14289 m³/s.

A river, whose annual flood peak can be represented by Gumbel distribution, has 100-year and 500-year return period flood of magnitude 9900 m3^3/s and 12100 m3^3/s respectively. The sample size is n = 30 (yˉn=0.536\bar{y}_n = 0.536, Sn=1.1124S_n = 1.1124).

i) What is the magnitude of 200 year and 1000 year flood?

ii) What are 95% and 80% confidence limits for 200 year and 1000 year flood if f(95%)=1.96f(95\%) = 1.96 and f(80%)=1.28f(80\%) = 1.28.

iii) A hydraulic structure of 25 year life was designed for 12300 m3^3/s peak flow. What is the hydrologic risk of the structure?

iv) What peak flow should be taken into consideration if you want the structure to be 99% reliable for a structure life of 25 years.

Gumbel's method: xT=xˉ+Kσx_T = \bar x + K\sigma, with K=yT−yˉnSnK=\dfrac{y_T-\bar y_n}{S_n} and yT=−ln⁡ ⁣[ln⁡TT−1]y_T = -\ln\!\left[\ln\dfrac{T}{T-1}\right]. Here yˉn=0.536\bar y_n=0.536, Sn=1.1124S_n=1.1124 (n = 30).

(i) 200-year and 1000-year floods

First find xˉ\bar x and σ\sigma from the two given floods.

T (yr)yTy_TK=(yT−0.536)/1.1124K=(y_T-0.536)/1.1124
1004.6003.653
5006.2145.104
σ=12100−99005.104−3.653=1516.8 m3/s,xˉ=9900−3.653×1516.8=4358.4 m3/s\sigma = \frac{12100-9900}{5.104-3.653} = 1516.8 \text{ m}^3/\text{s}, \qquad \bar x = 9900 - 3.653\times 1516.8 = 4358.4 \text{ m}^3/\text{s}
T (yr)yTy_TKKxT=xˉ+Kσx_T=\bar x+K\sigma (m³/s)
2005.2964.27910849
10006.9075.72713046

Answer: x200x_{200} ≈ 10849 m³/s and x1000x_{1000} ≈ 13046 m³/s.

(ii) Confidence limits

xT′=xT±f(c) Se,Se=b σN,b=1+1.3K+1.1K2x_T' = x_T \pm f(c)\,S_e, \qquad S_e = \frac{b\,\sigma}{\sqrt{N}},\quad b=\sqrt{1+1.3K+1.1K^2}
TKbSeS_e (m³/s)95% limits (f = 1.96)80% limits (f = 1.28)
2004.2795.16714318044 to 136539017 to 12680
10005.7276.67318489424 to 1666810680 to 15411

(iii) Hydrologic risk for 12300 m³/s, n = 25 years

K=12300−4358.41516.8=5.236,yT=0.536+5.236×1.1124=6.360K=\frac{12300-4358.4}{1516.8}=5.236, \quad y_T=0.536+5.236\times 1.1124=6.360 T=11−e−e−yT=579 yearsT=\frac{1}{1-e^{-e^{-y_T}}} = 579 \text{ years} Risk=1−(1−1T)n=1−(1−1579)25=0.0423\text{Risk}=1-\left(1-\frac1T\right)^{n}=1-\left(1-\frac{1}{579}\right)^{25} = 0.0423

Answer: T ≈ 579 years and the risk ≈ 4.2%.

(iv) Flood for 99% reliability over 25 years

Reliability = 0.99 means risk = 0.01.

1−1T=(0.99)1/25⇒T=2488 years1-\frac1T = (0.99)^{1/25} \Rightarrow T = 2488 \text{ years} yT=7.819,K=6.547,xT=4358.4+6.547×1516.8=14289 m3/sy_T=7.819, \quad K=6.547, \quad x_T = 4358.4 + 6.547\times1516.8 = 14289 \text{ m}^3/\text{s}

Answer: design peak flow ≈ 14289 m³/s (return period ≈ 2488 years).

✓ Answer checked by Claude Sonnet 5.5

All Engineering Hydrology chapter 6 questions →