Once the data series is identified and ranked and the plotting position is calculated, a graph of magnitude x versus probability [P(X > x), P(X < x), or T] can be plotted and a distribution fitted graphically. To facilitate this procedure, it is common to use some specially designed probability graph paper rather than linear graph paper. The probability scale in those special papers is chosen such that the resulting probability plot is a straight line. By plotting the data using a particular probability scale and constructing a best-fit straight line through the data, a graphic fit is made to the distribution used in constructing the probability scale. This is a graphic approach to estimate the statistical parameters of the distribution.
Example 3.1 illustrates the graphic approach to the analysis of flood data. The general procedure is as follows: [2] [3] 3 [4]
5. Extend the line to the highest return-period value needed, and read all required return-period values off the line.
Example 3.1 The Boneyard Creek stream gauging station was located near the fire station on the campus of the University of Illinois at Urbana-Champaign. From the USGS Water Supply Papers, the partial duration data of peak discharges above 400 ft3/s between the water years 1961 and 1975 were obtained and listed below. In addition, for the years when there was no flow in a year exceeding 400 ft3/s, the peak flow for that year is given in parenthesis (e. g., 1961).
Year |
Discharge, ft3/s |
Year |
Discharge, ft3 /s |
1961 |
(390) |
1969 |
549, 454 |
1962 |
(374) |
1970 |
414, 410 |
1963 |
(342) |
1971 |
434, 524 |
1964 |
507 |
1972 |
505, 415, 406 |
1965 |
579, 406, 596 |
1973 |
428, 447, 407 |
1966 |
416 |
1974 |
468, 543, 441 |
1967 |
533 |
1975 |
591, 497 |
1968 |
505 |
(a) List the ranked annual maximum series. Also compute and list the corresponding plotting positions (return period) and exceedance probability P (X > x).
(b) Plot the annual maximum series on (i) Gumbel paper and (ii) lognormal paper.
(c) Construct a best-fit line through the nonlinear plots, and estimate the flows for return periods of 2, 10, 25, and 50 years.
Solution n = 15 (a)
Annual Maximum Discharge (ft3/s) |
Rank (m) |
‘T _ n+1 Tm = m (years) |
P(X > X(m)) = 1/Tm |
P (X < x(m)) = 1 — 1/Tm |
596 |
1 |
16.00 |
0.0625 |
0.9375 |
591 |
2 |
8.00 |
0.1250 |
0.8750 |
549 |
3 |
5.33 |
0.1875 |
0.8125 |
543 |
4 |
4.00 |
0.2500 |
0.7500 |
533 |
5 |
3.20 |
0.3125 |
0.6875 |
524 |
6 |
2.67 |
0.3750 |
0.6250 |
507 |
7 |
2.29 |
0.4375 |
0.5625 |
505 |
8 |
2.00 |
0.5000 |
0.5000 |
505 |
9 |
1.78 |
0.5625 |
0.4375 |
447 |
10 |
1.60 |
0.6250 |
0.3750 |
416 |
11 |
1.46 |
0.6875 |
0.3125 |
414 |
12 |
1.33 |
0.7500 |
0.2500 |
390 |
13 |
1.23 |
0.8125 |
0.1875 |
374 |
14 |
1.14 |
0.8750 |
0.1250 |
342 |
15 |
1.06 |
0.9375 |
0.0625 |
(b) Plots of the annual maximum flow series on the Gumbel and lognormal probability papers are shown in Fig. 3.2.
|
||||||||||||||||
|
||||||||||||||||
|