Probable Maximum Precipitation Estimation Using the Revised -Value Method in Hong Kong
Publication: Journal of Hydrologic Engineering
Volume 22, Issue 8
Abstract
A brief overview of a statistical method to estimate the probable maximum precipitation (PMP) is presented. This study addresses some issues associated with Hershfield’s -value method to estimate PMP in China, which can be solved by the revised Hershfield’s -value method. This new derivation makes it clear that the frequency factor is dependent on only two variables, the standardized variable, , the maximum deviation from the mean scaled by its standard deviation, and the sample size, . It is found that there is a consistent relationship between and . Therefore, can be used to make a preliminary estimate of PMP under some conditions when sufficient rainfall data are available. The advantages and disadvantages of this revised -value method are also discussed here with a case study for the estimation of 24 h PMP in Hong Kong. The 24 h PMP estimate in Hong Kong based on the local rainfall data is approximately 1,753 mm.
Introduction
Probable maximum precipitation (PMP) has been defined as “the greatest depth of precipitation for a given duration meteorologically possible for a design watershed or a given storm area at a particular location at a particular time of year, with no allowance made for long term climatic trends” (WMO 2009). Hydrologists use PMP to estimate the probable maximum flood (PMF) in the design of a particular project in a given watershed, such as the height and reservoir storage capacity of a dam or dimension of flood-carrying structures (spillway and flood carrying tunnel). In this way, the risks of loss of life and damage could be reduced.
The available methods for PMP estimation are mainly categorized in two ways: hydrometeorological estimation methods and statistical methods in practice. Hydrometeorological estimation methods can be generally classified as follows (WMO 2009; Lin and AECOM 2014): a storm model approach (Collier and Hardaker 1996), moisture maximization (Micovic et al. 2015), a generalized method, storm transposition, and a depth-area-duration (D-A-D) method. The statistical estimation is a possible approach, a kind of modified frequency analysis based on historical precipitation records. Among several statistical procedures, Hershfield’s -value method is the best known (National Research Council 1985; WMO 1986, 2009). It was initially developed by Hershfield (1963), who defined a frequency factor omitting the observed maximum value from a sample series for the estimation of PMP. Hershfield combined samples by either including or excluding the maximum value and produced some nomographs with adjustments. Koutsoyiannis and Papalexiou (2016) set up mathematical expressions of Hershfield’s adjustment factors through these nomographs. Subsequently, Hershfield (1965) improved the method by constructing an empirical graph with varying between 5 and 20 depending on the mean value of annual rainfall maxima and rainfall duration. Preliminary appraisal of this technique in Canada (Bruce and Clark 1966) and in the U.S. (Myers 1967) showed that the PMP estimates obtained by this approach are closely comparable to those obtained by the conventional moisture maximization and storm transposition methods. According to Wiesner (1970), this method has the advantage of taking into account the entire rainfall data set, expressing it in terms of statistical parameters. Using the same data set as Hershfield used, Koutsoyiannis (1999) fitted a generalized extreme value (GEV) to the frequency factors and found that the highest value, 15, of corresponds to a 60,000 year return period event. Papalexiou and Koutsoyiannis (2006) studied the maximum precipitation depths derived by moisture maximization at a few stations in the Netherlands. They concluded that a probabilistic approach to estimating extreme precipitation values is more consistent with natural behavior and provides a better basis for estimation. The World Meteorological Organization (WMO) in its various manuals and technical publications (WMO 1969, 1986, 2009) has also recommended this method for estimation of point PMP for those river basins whose daily rainfall data are available for a long period of time.
However, a standardized variable, , the maximum deviation from the mean of a sample scaled by the standard deviation of the sample, derived directly from the frequency equation, was used to replace in China. Its envelope value was used by some researchers to calculate the PMP (Wang 1999; Hua et al. 2007) in China. They recommended that was more reasonable and consistent with the frequency equation. In contrast, Lin (1981) and Lin and Vogel (1993) recommended that replace under some circumstances, and Lin used this revised method to estimate the PMP of the Zhelin Reservoir on the Xiushui River in Jiangxi Province, China, in 1981 and in Austin, Texas, U.S., in 1993.
Regarding the fact that Hershfield’s -value method has been used in China, the purpose of this study was to determine the differences and relationships between the frequency factor defined by Hershfeild and the standardized variable . In this paper, Hershfield’s -value method will be reviewed with some criteria developed for using the method in a reasonable way. The advantages and disadvantages of this revised -value method will also be discussed in connection with a case study on the estimation of 24 h PMP in Hong Kong.
Revised -Value Method
Hershfield’s procedure of estimating PMP was originally developed based on the general frequency equation (Chow 1951)where = rainfall for return period ; and = mean and standard deviation of a series of annual maximum rainfall values for a given duration; and = frequency factor associated with and the probability distribution selected to fit the data. If the maximum observed rainfall, , is substituted for , for , Eq. (1) could be as follows:
(1)
(2)
Then is the number of standard deviations to be added to the mean to obtain . is computed using the equationwhere = highest value for this series; and and are respectively the mean and standard deviation excluding the highest value from the series. By enveloping from a large number of computed , Hershfield estimated PMP using Eq. (4):
(3)
(4)
The statistical expressions for PMP are rederived to obtain a better understanding of the method. The standardized variable is introduced to replace from Eq. (1), soandwhere = standardized variable or maximum deviation from mean scaled by its standard deviation . It is computed directly from the equation
(5)
(6)
(7)
For the distribution of a series of annual rainfall maxima we havewhere and = expected value and population standard deviation, respectively. is the maximum value of . To estimate , the population parameters in Eq. (8) must be substituted, i.e., by and by . The true value of the population parameter would be unknown before happened. If the estimate for is the same as that presented by Hershfield, then is defined as the number of standard deviations to be added to the mean to obtain of the incomplete sample series (), andwhere and = mean and standard deviation of the incomplete sample series, respectively. Thus, the equation for calculating isThe relationship between and is . The definition of is . Combining the foregoing two equations with Eqs. (7) and (10), Eq. (11) can be derived. The detailed derivations of Eq. (11) are included in the Appendix:
(8)
(9)
(10)
(11)
Using Eq. (12) [or Eq. (13)], along with the obvious inequalities and , the following inequalities, which restrict the range of and , can be derived:
(14)
(15)
The lower limit of is approximately 1, while the upper limit of is approximately . Obviously, is always greater than . However, the difference in value between and gets smaller as increases. This means that the frequency factor approaches the upper limit when the sample size increases. For example, in the case of and , Tables 1 and 2 show the variations of . As increases, approaches . Therefore, statistically, has a consistent relationship with and can be used in the estimation of PMP. For practical use, Eq. (4) is applied to estimate PMP. Because the standard deviation is defined as the coefficient of variation multiplied by the mean , Eq. (4) can be written
(16)
27 | 0.9644 | 0.0385 | 5 | 132.378 |
50 | 0.9804 | 0.0204 | 5 | 7.293 |
150 | 0.9934 | 0.0067 | 5 | 5.503 |
500 | 0.9980 | 0.0020 | 5 | 5.136 |
1000 | 0.9990 | 0.0010 | 5 | 5.066 |
5000 | 0.9998 | 0.0002 | 5 | 5.013 |
Note: ; ; .
102 | 0.9903 | 0.0099 | 10 | 1014.938 |
200 | 0.9950 | 0.0050 | 10 | 14.249 |
200 | 0.9967 | 0.0033 | 10 | 12.288 |
600 | 0.9983 | 0.0017 | 10 | 10.967 |
1000 | 0.9990 | 0.0010 | 10 | 10.547 |
5000 | 0.9998 | 0.0002 | 10 | 10.103 |
Note: ; ; .
Since Eq. (12) requires , should be greater than or equal to . Therefore, the minimum data size is deduced a criterion to use the -value eligibly
(17)
From Eq. (12) [or Eq. (13)] the following relation can also be derived:where is defined as the stable data size needed to give statistically reliable results for a particular . If the ratio of to has a value of 1.1, or a relative error of 10%, the relation in Eq. (18) can be further approximated as
(18)
(19)
Eq. (19) can be used to perform a quick PMP estimate with respect to the length of data required to develop a stable estimate. For instance, it requires 104 years of data records as stable data size if , referring to Fig. 1. Furthermore, when , it could raise the uncertainties of , usually up to 50%. Hence, these stations should be discarded to reduce the computation burden as such stations will contribute nothing to the statistical results except errors.
The distribution of the sample mean, , approaches normal distribution as increases (Lin 1978). Let the correction factor of the population mean to the sample mean be , which can be expressed as
(20)
If there was no sampling error, the would be exactly equal to 1. However, due to sampling error may be greater or less than 1 with equal chance. By taking a confidence level of 99.7%, i.e., the is varying within a confidence interval of or . Then could be expressed as for maximum adjustment, and the adjustment iswhere = adjusted sample mean of . The relations of with and are given in Table 3. It shows that decreases toward 1 while decreases when is fixed; and decreases toward 1 while increases when is certain.
(21)
Parameters | ||||||||
---|---|---|---|---|---|---|---|---|
20 | 30 | 40 | 50 | 100 | 500 | 1,000 | 10,000 | |
1.670 | 1.547 | 1.474 | 1.424 | 1.300 | 1.313 | 1.094 | 1.030 | |
1.268 | 1.219 | 1.189 | 1.169 | 1.120 | 1.053 | 1.038 | 1.012 | |
1.067 | 1.054 | 1.047 | 1.042 | 1.030 | 1.013 | 1.009 | 1.003 |
Precipitation data are usually observed at fixed time intervals. For example, 8 a.m. to 8 a.m. means daily, and such data rarely yield the true maximum rainfall amounts for the indicated duration 24 h. As we know, the annual maximum observational day rainfall amounts are very likely to be appreciably less than the annual maximum 24 h rainfall amounts determined by 1,440 consecutive minutes with unrestricted beginning and ending time. Studies of thousands of station-years of rainfall data indicate that a conversion factor of 1.13 to convert 1-day rainfall data to 24 h is reasonable to be used in practice (Weiss 1964).
Now, if we use this revised -value method to estimate PMP at a single station or over a small watershed, it should be taken by the following steps:
1.
Sort the station data in descending order based on the maximum value, ;
2.
Check the required minimum size for each station. If , then the station is discarded;
3.
Check the required stable size , which should meet the inequalities ;
4.
Calculate the values of , , and for each station;
5.
Make the “maximum adjustment” of in terms of sampling error. Furthermore, if precipitation data are observed at fixed time interval, the rainfall should be multiplied by a conversion factor of 1.13; and
6.
Regionalize or pick up the highest values of , and of Eq. (16) over the eligible stations in the concerned region and estimate the PMP.
It is noted that the essence of the reviewed -value method is storm transposition, but instead of transposing the specific rainfall amount of one storm, an abstracted statistic is transposed (WMO 2009). In the study, other parameters such as and are also regionalized within the study area for the sake of maximization.
Application and Results
In Hong Kong, PMP was first derived by Bell and Chin (1968). They used depth-area-duration (DAD) method together with moisture maximization to estimate the 24 h PMP from 21 major storms between 1955 and 1965. They further updated their PMP estimates through the introducing a seasonal adjustment factor. Hong kong Observatory (1999); Chang and Hui (2001) used the same method but without seasonal adjustment to analyze 53 storms occurring between 1966 and 1999, and got a result which was smaller than the former one. The statistical method was not used in the PMP estimates in Hong Kong before. Hence, the revised -value method is applied to estimate 24 h PMP in Hong Kong for the first time. The 5 min data or hourly data of 64 rain-gauge stations in Hong Kong were selected. The annual maximum precipitation series of 24 h for each rain-gauge station was obtained by using a 24 h sliding window over the 5 min or hourly data. Since Shenzhen is very close to Hong Kong, the annual maximum precipitation series of Shenzhen acquired from the Hydrologic Bureau of Guangdong Province was also used in the study. Table 4 gives us the detailed information of all the rain-gauge stations. Fig. 2 shows the locations of the rain-gauge stations.
Number | Station identifier | Period of time | Type |
---|---|---|---|
1–19 | H01, H02, H03, H04, H05, H06, H07, H08, H09, H10, H12, H14, H15, H16, H17, H18, H19, H20, H21 | 1984–2010 | 5 min |
20–27 | K01, K02, K03, K04, K05, K06, K07, K08 | 1984–2010 | 5 min |
28–42 | N01, N02, N03, N04, N05, N06, N07, N08, N09, N10, N11, N12, N13, N14, N15 | 1984–2010 | 5 min |
43 | HKO | 1987–2010 | 5 min |
1885–1986 (1940–1946 data missing due to the World War II) | 1 h | ||
44 | N16 | 1985–2010 | 5 min |
45–61 | R11, R12, R13, R14, R17, R18, R21, R22, R23, R24, R25, R26, R27, R28, R29, R31, R32 | 1987–2010 | 5 min |
62 | R30 | 1987–2008 | 5 min |
63–64 | N17, N18 | 1991–2010 | 5 min |
65 | Shenzhen | 1981–2010 | Annual maximum 24 h rainfall |
As described above, firstly, each rain-gauge station was displayed in descending order of maximum value, . Secondly, the minimum data size, was calculated for each. And all the rain-gauge stations meet the criterion. Thirdly, by applying the criterion of a stable data size , N09 and N15 were discarded. Following, the values of , , and for each eligible rain-gauge station are calculated. Since the annual maximum precipitation series for 24 h is directly obtained by a 24 h sliding window, which is true maxima for 24 h. There was no need to apply a conversion factor of 1.13 to the mean value of stations in Hong Kong. However, the mean daily value series of Shenzhen was multiplied by 1.13 to convert 1-day data to 24 h. Furthermore, each eligible rain-gauge station need a maximized adjustment of mean in terms of sampling error at a confidence level of 99.7%. Finally, the highest values of , , and were obtained. They were 426.30 (of N14), 5.46 (of N14), 0.57 (of R11), respectively. Therefore, according to the Eq. (16), the 24 h PMP estimate for Hong Kong is approximately 1,753 mm.
Table 5 shows the results of parameters of all those rain-gauge stations in descending order of . It is shown that the highest values of and come from the same station N14. As we know, rainfall is quite affected by topography, and rainfall in mountain areas is mostly larger than in plain areas. Observation value at automatic station may be a little higher than the ordinary observation station. N14 is an automatic station located in Tai Mo Shan, the highest elevation in Hong Kong. Therefore, the rainfall amount at station N14 is always larger than the other rain-gauge stations.
Number | Stations | Years | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
1 | N14 | 27 | 956.0 | 0.53 | 426.30 | 5.46 | 3.66 | 16 | 89 | 94 |
2 | N09 | 27 | 800.0 | 0.52 | 334.41 | 7.20 | 4.09 | 19 | 108 | 94 |
3 | N17 | 20 | 745.0 | 0.53 | 389.95 | 4.32 | 2.99 | 11 | 64 | 70 |
4 | R11 | 24 | 735.5 | 0.57 | 365.11 | 4.01 | 3.01 | 12 | 64 | 84 |
5 | HKO | 119 | 697.1 | 0.42 | 233.36 | 5.26 | 4.71 | 25 | 261 | 416 |
6 | N02 | 27 | 587.5 | 0.42 | 313.82 | 4.19 | 3.18 | 13 | 70 | 94 |
7 | N01 | 27 | 570.0 | 0.41 | 331.02 | 3.38 | 2.77 | 10 | 56 | 94 |
8 | N15 | 27 | 562.0 | 0.37 | 284.04 | 5.82 | 3.76 | 17 | 94 | 94 |
9 | K02 | 27 | 508.0 | 0.41 | 319.38 | 2.72 | 2.36 | 8 | 44 | 94 |
10 | N06 | 27 | 508.0 | 0.39 | 320.70 | 2.78 | 2.40 | 8 | 45 | 94 |
11 | K06 | 27 | 506.5 | 0.40 | 301.83 | 3.21 | 2.67 | 10 | 53 | 94 |
12 | R26 | 24 | 504.5 | 0.44 | 315.02 | 2.78 | 2.36 | 8 | 44 | 84 |
13 | H01 | 27 | 496.0 | 0.43 | 278.82 | 3.51 | 2.84 | 11 | 58 | 94 |
14 | N10 | 27 | 492.0 | 0.42 | 295.04 | 3.04 | 2.57 | 9 | 50 | 94 |
15 | R17 | 24 | 487.5 | 0.46 | 262.94 | 3.96 | 2.99 | 11 | 64 | 84 |
16 | H17 | 27 | 486.0 | 0.36 | 303.23 | 3.11 | 2.61 | 9 | 51 | 94 |
17 | K04 | 27 | 484.5 | 0.35 | 297.90 | 3.28 | 2.71 | 10 | 54 | 94 |
18 | H02 | 27 | 483.5 | 0.40 | 296.01 | 2.97 | 2.52 | 9 | 49 | 94 |
19 | H04 | 27 | 477.0 | 0.37 | 310.63 | 2.70 | 2.34 | 8 | 44 | 94 |
20 | K07 | 27 | 473.5 | 0.39 | 312.31 | 2.47 | 2.19 | 7 | 40 | 94 |
21 | N04 | 27 | 470.5 | 0.37 | 288.06 | 3.18 | 2.65 | 10 | 52 | 94 |
22 | N18 | 20 | 470.0 | 0.40 | 322.00 | 2.53 | 2.14 | 7 | 38 | 70 |
23 | N08 | 27 | 468.5 | 0.34 | 292.98 | 3.26 | 2.70 | 10 | 54 | 94 |
24 | H08 | 27 | 468.0 | 0.37 | 300.31 | 2.80 | 2.41 | 8 | 45 | 94 |
25 | H16 | 27 | 468.0 | 0.35 | 299.25 | 2.98 | 2.52 | 9 | 49 | 94 |
26 | N11 | 27 | 467.5 | 0.39 | 274.52 | 3.42 | 2.79 | 10 | 57 | 94 |
27 | R12 | 24 | 467.5 | 0.41 | 274.30 | 3.49 | 2.77 | 10 | 56 | 84 |
28 | N03 | 27 | 466.0 | 0.36 | 299.10 | 2.83 | 2.43 | 8 | 46 | 94 |
29 | H12 | 27 | 465.0 | 0.35 | 308.29 | 2.69 | 2.34 | 8 | 43 | 94 |
30 | H14 | 27 | 460.0 | 0.32 | 257.23 | 4.92 | 3.48 | 15 | 82 | 94 |
31 | N12 | 27 | 458.5 | 0.43 | 265.00 | 3.29 | 2.71 | 10 | 54 | 94 |
32 | H06 | 27 | 457.5 | 0.38 | 300.48 | 2.60 | 2.27 | 8 | 42 | 94 |
33 | H10 | 27 | 457.5 | 0.35 | 305.18 | 2.61 | 2.28 | 8 | 42 | 94 |
34 | R21 | 24 | 453.0 | 0.44 | 287.06 | 2.65 | 2.27 | 8 | 42 | 84 |
35 | N13 | 27 | 452.0 | 0.31 | 287.52 | 3.35 | 2.75 | 10 | 55 | 94 |
36 | K01 | 27 | 447.5 | 0.40 | 294.32 | 2.46 | 2.18 | 7 | 39 | 94 |
37 | R22 | 24 | 435.5 | 0.45 | 249.38 | 3.44 | 2.74 | 10 | 55 | 84 |
38 | R25 | 24 | 434.0 | 0.36 | 258.01 | 3.85 | 2.94 | 11 | 62 | 84 |
39 | R27 | 24 | 434.0 | 0.42 | 283.00 | 2.56 | 2.21 | 7 | 40 | 84 |
40 | H09 | 27 | 433.0 | 0.34 | 307.60 | 2.23 | 2.00 | 7 | 35 | 94 |
41 | H07 | 27 | 432.0 | 0.36 | 289.31 | 2.53 | 2.23 | 7 | 41 | 94 |
42 | K03 | 27 | 426.0 | 0.39 | 287.66 | 2.34 | 2.09 | 7 | 37 | 94 |
43 | H03 | 27 | 424.0 | 0.36 | 261.71 | 3.20 | 2.66 | 10 | 53 | 94 |
44 | K05 | 27 | 420.0 | 0.35 | 281.69 | 2.61 | 2.28 | 8 | 42 | 94 |
45 | N05 | 27 | 419.5 | 0.37 | 266.42 | 2.86 | 2.45 | 8 | 47 | 94 |
46 | H19 | 27 | 419.0 | 0.33 | 278.32 | 2.77 | 2.39 | 8 | 45 | 94 |
47 | R23 | 24 | 415.5 | 0.37 | 294.30 | 2.22 | 1.97 | 6 | 34 | 84 |
48 | H18 | 27 | 410.5 | 0.31 | 281.78 | 2.63 | 2.30 | 8 | 42 | 94 |
49 | R28 | 24 | 405.0 | 0.43 | 275.31 | 2.25 | 2.00 | 6 | 35 | 84 |
50 | N07 | 27 | 399.5 | 0.47 | 268.74 | 2.09 | 1.90 | 6 | 33 | 94 |
51 | H21 | 27 | 395.5 | 0.35 | 253.94 | 2.95 | 2.51 | 9 | 48 | 94 |
52 | R32 | 24 | 393.5 | 0.32 | 259.81 | 3.10 | 2.55 | 9 | 49 | 84 |
53 | R14 | 24 | 390.5 | 0.33 | 254.26 | 3.14 | 2.57 | 9 | 50 | 84 |
54 | K08 | 27 | 390.0 | 0.36 | 277.45 | 2.16 | 1.95 | 6 | 34 | 94 |
55 | H05 | 27 | 388.5 | 0.33 | 273.22 | 2.35 | 2.09 | 7 | 37 | 94 |
56 | Shenzhen | 29 | 386.2 | 0.35 | 278.28 | 2.91 | 2.51 | 9 | 48 | 101 |
57 | H15 | 27 | 385.5 | 0.29 | 253.10 | 3.21 | 2.67 | 10 | 53 | 94 |
58 | N16 | 26 | 377.5 | 0.35 | 281.33 | 1.93 | 1.76 | 6 | 30 | 91 |
59 | R29 | 24 | 374.0 | 0.40 | 237.26 | 2.84 | 2.39 | 8 | 45 | 84 |
60 | R31 | 24 | 364.5 | 0.40 | 237.71 | 2.66 | 2.28 | 8 | 42 | 84 |
61 | R24 | 24 | 360.5 | 0.30 | 272.35 | 2.11 | 1.89 | 6 | 33 | 84 |
62 | H20 | 27 | 355.5 | 0.33 | 251.20 | 2.31 | 2.07 | 7 | 37 | 94 |
63 | R18 | 24 | 334.5 | 0.34 | 233.09 | 2.48 | 2.16 | 7 | 39 | 84 |
64 | R30 | 22 | 332.5 | 0.37 | 236.94 | 2.27 | 1.99 | 6 | 35 | 77 |
65 | R13 | 24 | 284.5 | 0.31 | 217.60 | 1.97 | 1.79 | 6 | 30 | 84 |
Note: The values in bold are the maximum value.
Discussions and Conclusions
Although the frequency factor is always greater than the standardized variable, , it usually rapidly approaches when increases. In this study, it has already proved that had an obvious statistical relationship with . Therefore, some stations are discarded by applying the criteria of the minimum sample size and the stable sample size to reduce the computation burden. And the highest values of , , and are calculated and selected in parameter regionalization. Finally, the PMP at a single station or over a small watershed could be estimated by applying the Eq. (16).
If the selected data length is less than , the revised -value method could not be applied. For example, if , the minimum data size of is at least 27 years. Although representing a stable data length, it cannot guarantee a true estimation of PMP. The criterion just presents a statistically reliable estimation of how long the data series should be. Obviously, the longer data series is better for the method. The estimate of the revised -value method greatly depends on the data availability, including data quality and quantity. The PMP estimate will be more precise with a longer historical data series. In fact, there are few stations in the real world that are long enough to work out reliable PMP estimates. In general, a station having longer data with outlier occurring in history is very welcome for the statistical approach. So, caution must be taken when using this method to estimate PMP in practice.
Although the data used in the revised -value method comes from the whole study area, the PMP result is still a point PMP and mostly to be representative to the potential storm center without temporal and spatial distribution. The major advantage of this approach over others is to provide a quick and preliminary estimate of PMP at a single station or over a small watershed. The suggested PMP estimates of the revised -value method could only be considered as a reference value. It could not be recommended as final estimate for engineering design study. And it will be comparable with values of storm transposition or DAD method.
The 24 h PMP estimate by the revised -value method in Hong Kong is 1,753 mm. It is substantially greater than 1,250 mm given by Chang and Hui (2001). The difference of the two is about 503 mm. The reason for the difference is that they are estimated by totally different method based on different data. The revised -value method is a statistical way with the series of annual maxima rainfall till 2010. While the DAD method and moisture maximization utilized by Chang is a hydrometeorological method with the storm data and dew point data from 1955 to 1999 applied. Per the criteria presented above, if researchers want to get a stable estimate with 10% error in , the data length should be longer than 89 years and 64 years for N14 and R11, respectively. But in this study, historical data are shorter than 30 years for a great majority of Hong Kong rain-gauge stations except for HKO station (119 years available). Therefore, this method cannot guarantee to yield a stable PMP estimate. Furthermore, the storm transposition and DAD method are recommended to apply to estimate the 24 h PMP in Hong Kong in the future if possible to get comparable results.
Supplemental Data
Appendix S1 is available online in the ASCE Library (www.ascelibrary.org).
Supplemental Materials
File (supplemental_data_he.1943-5584.0001517_lin.pdf)
- Download
- 157.85 KB
Appendix. Derivations of Equation
The Hershfield’s procedure of estimating PMP was originally developed based on the general frequency equation (Chow 1951):where = rainfall for return period ; and = mean and standard deviation of a series of annual maximum rainfall values for a given duration, respectively; is defined as a frequency factor, which is associated with and the probability distribution selected to fit the data. If the maximum observed rainfall , is substituted for , for , Eq. (22) could be written as
(22)
(23)
Then is the number of standard deviation to be added to the mean to obtain . is computed by the following equationwhere = highest value for this series; and are, respectively, the mean and the standard deviation excluding the highest value from the series. By enveloping from a large number of computed , Hershfield estimated Probable Maximum Precipitation (PMP) though Eq. (25)
(24)
(25)
The statistical expressions for PMP are re-derived to obtain a better understanding of the method. The standardized variable is introduced to replace the of Eq. (22), soandwhere = standardized variable or the maximum deviation from the mean, scaled by its standard deviation . It is computed directly from the following equation
(26)
(27)
(28)
As known, the relation between and is
(29)
And according to the definition of
(30)
(31)
Because , so Eq. (31) could be
(32)
Both sides of the Eq. (35) were divided by :
(36)
Because both and are always greater than 0. Making , , Eq. (38) could be rewritten as follows:
(39)
Acknowledgments
This work was supported by the 24 h PMP Updating Study Project of the Civil Engineering and Development Department of the Hong Kong Special Administrative Region (HKSAR). The data used in the paper were provided by Hong Kong Observatory (HKO).
References
Bell, G. J., and Chin, P. C. (1968). “The probable maximum rainfall in Hong Kong.”, Royal Observatory, Hong Kong.
Bruce, J. P., and Clark, R. H. (1966). Introduction to hydrmeteorology, Pergaman Press, New York.
Chang, W. L., and Hui, T. W. (2001). “Probable maximum precipitation for Hong Kong.” Hong Kong Observatory, Hong Kong.
Chow, V. T. (1951). “A general formula for hydrologic frequency analysis.” Trans. Am. Geophys. Union, 32(2), 231–237.
Collier, C. G., and Hardaker, P. J. (1996). “Estimating probable maximum precipitation using a storm model approach.” J. Hydrol. 183, 277–306.
Hershfield, D. M. (1963). “Estimating the probable maximum precipitation.” Trans. Am. Soc. Civ. Eng., 128(1), 534–551.
Hershfield, D. M. (1965). “Method for estimating probable maximum precipitation.” J. Am. Water Works Assoc., 57(8), 965–972.
Hong Kong Observatory. (1999). “The 24 h probable maximum precipitation updating study for Hong Kong.” Hong Kong.
Hua, J., Huang, Y., Yang, H., and Liu, L. (2007). “Statistical estimate magnification method for deducing probable maximum precipitation.” J. Hohai Univ., 35(3), 255–257.
Koutsoyinannis, D. (1999). “A probable view of Hershfield’s method for estimating probable maximum precipitation.” Water Resour. Res., 35(4), 1313–1322.
Koutsoyinannis, D., and Papalexiou, S. M. (2016). “Extreme rainfall: Global perspective.” Chow’s handbook of applied hydrology, 2nd Ed., McGraw-Hill, New York.
Lin, B. (1981). “Application of statistical estimation in study of probable maximum precipitation.” J. Hohai Univ., 1, 52–62.
Lin, B., and AECOM. (2014). “24 h PMP updating study, technical report, study of landslides occurring in kowloon and the new territories in 2012 and 2013—Feasibility study Hong Kong.” Hong Kong.
Lin, B., and Vogel, J. (1993). “A new look at the statistical estimation of PMP.” Engineering hydrology, ASCE, Reston, VA.
Lin, S. G. (1978). Basic probability and mathematical statistics, People’s Education Press, Beijing.
Micovic, Z., Schaefer, M. G., and Taylor, G. H. (2015). “Uncertainty analysis for probable maximum precipitation estimates.” J. Hrdrol., 521, 360–373.
Myers, V. A. (1967). “Meteorological estimation of extreme precipitation for spillway design floods.”, U.S. Dept. of Commerce, Environmental Science Services Administration, Weather Bureau, Boulder, CO.
National Research Council. (1985). Safety of dams, flood and earthquake criteria, National Academy Press, Washington, DC.
Papalexiou, S. M., and Koutsoyiannis, D. (2006). “A probabilistic approach to the concept of probable maximum precipitation.” Adv. Geosci., 7, 51–54.
Wang, W. (1999). Principles and methods of PMP/PMF calculations, China Water & Power Press, Beijing.
Weisner, C. J. (1970). Hydrometeorology, Chapman and Hall, London.
Weiss, L. L. (1964). “Ratio of true to fixed-interval maximum rainfall.” J. Hydraul. Div., 90, 77–82.
WMO (World Meteorological Organization). (1969). “Estimation of maximum floods.” Geneva.
WMO (World Meteorological Organization). (1986). “Manual for estimation of probable maximum precipitation (PMP).” Geneva.
WMO (World Meteorological Organization). (2009). “Manual for estimation of probable maximum precipitation.” Geneva.
Information & Authors
Information
Published In
Copyright
This work is made available under the terms of the Creative Commons Attribution 4.0 International license, http://creativecommons.org/licenses/by/4.0/.
History
Received: Jun 10, 2016
Accepted: Dec 14, 2016
Published online: Mar 29, 2017
Published in print: Aug 1, 2017
Discussion open until: Aug 29, 2017
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.