A 60-Minute Design Rainstorm for the Urban Area of Yangpu District, Shanghai, China

: Rainfall with varied temporal distribution is an essential input to urban ﬂood models. In this study, a 60-min design rainstorm with different return periods for the urban area of Yangpu District, Shanghai, China was derived. The design of areal rainfall amounts with given return periods was calculated through frequency analysis. The temporal distribution of the hyetograph was derived using the Pilgrim and Cordery method, combined with the fuzzy identiﬁcation of seven mode hyetographs for single-peak and double-peak rainstorms separately. The derived hyetographs using the Pilgrim and Cordery method were compared with the classic Chicago rainstorm method. The results indicated that: (1) separating single-peak and double-peak rainstorms to derive respective hyetographs is more practical and rational; (2) a design rainstorm using the proposed methodology is superior to the Chicago rainstorm method.


Introduction
The inevitable trend of global climate change has brought about considerably frequent occurrences of extreme climate events [1][2][3][4]. Recently, under the influence of typhoons, urban floods caused by rainstorms have exerted great harm on the safety of society. In order to ensure the safety and continued operation of cities under extreme climate conditions, the Chinese government is making every effort to solve the problem of urban flooding by constructing and maintaining drainage pipelines. In order to determine the slope and diameter of the pipeline when designing drainage pipelines, models that simulate urban waterlogging are developed. Rainfall is the input of most kinds of hydrological and hydraulic models, whose accuracy and resolution directly affect the correctness of the model output [5,6]. Usually, a design rainstorm with a given return period is used in such models. The total rainfall amount, with specific duration, can be calculated through frequency analysis. In early studies, a hyetograph was assumed to be uniformly distributed on both spatial and temporal scales [7]. However, this assumption is in conflict with reality because in most rainfall events the rainfall intensity is always variable during the whole process [8]. Therefore, the refinement of the temporal and spatial patterns of hyetographs is more conducive to forecasting urban floods.
The precipitation data obtained by rain gauges only account for observation at one point. In order to study the total rainfall amount in a region, areal rainfall should be calculated. The calculation of areal rainfall is normally carried out through the transformation from "point" to "area" by means of various interpolation methods. Conventional areal rainfall estimation methods are the arithmetic average method, Thiessen polygon techniques, the isohyetal method, and the grid method. Nowadays, the development of remote sensing has brought rainfall data from weather radars and satellites with much higher coverage and finer resolution, which makes the calculation of areal rainfall more accurate and convenient [9][10][11]. Examples of such products are TRMM, 3B42/3B43, CMORPH, PERSIANN, TAMSA, and RFE2.0.
Research into design hyetographs with varying temporal patterns started in the early 20th century. There are several classic and widely used methodologies to derive hyetographs. In 1937, Mo et al. [12] analyzed 800 rainfall events in the former Soviet Union and summarized seven types of hyetographs. They are well-known as "mode hyetographs". Keifer and Chu [13] proposed a method which is based on an intensity-duration-frequency (IDF) curve and a rainstorm intensity formula. This method is known as the Chicago rainstorm method. Huff [14] divided rainfall events into four classes based on which quarter of the duration the peak falls in. The hyetograph is then derived from quantiles of rainfall events within each class in a different period. Pilgrim and Cordery [15] developed an order-average hyetograph based on the principle of statistics (referred to as the P&C method). It places the peak of the rainfall in the period where it is most likely to happen. For a small drainage area that is below 10 mi 2 , Yen and Chow [16] proposed a triangular-shaped hyetograph. The position of the peak is determined by the dimensionless first-order moment of the rainstorm event. The Natural Resources Conservation Service (NRCS) [17] developed a synthetic hyetograph, utilizing 6 h and 24 h rainstorm data from the US, which is called a "Soil Conservation Service (SCS) hyetograph". There are four basic types of SCS hyetographs, each of them accounting for a specific region and climate. Besides the above-mentioned classic methods, there are also newly-developed methods to derive hyetographs; refer to Cheng et al. [18], Lin and Wu [19], Powell et al. [20], Lee and Ho [21], and Kottegoda [22] for detailed descriptions.
The classification of seven mode hyetographs is simple and distinct. Among them, Types I, II, and III are single-peak hyetographs, with peak falls at the beginning, middle, and end of the event; Type IV is uniformly distributed; V, VI, and VII are double-peak hyetographs with peaks falling at the beginning and end, beginning and middle, middle and end of the event, respectively. Visual estimation and fuzzy identification are often used for classification.
As the representative of statistical hyetographs, the P&C method is based on the statistical analysis of considerable rainfall data, and thus the derived hyetograph approximates the actual rainfall process. The simplicity of its principle and implementation makes it a convenient method. In most cases, all the rainfall events are used to derive one single hyetograph. However, the rainfall events may have a single peak, double-peaks, or multiple peaks. These rainfall events with different numbers of peaks will have different influences on the confluence process. If one mixes single-peak and multiple-peak rainfalls, the derived hyetograph is distorted to some extent. Therefore, in this study, a combination of mode hyetographs and the P&C method is proposed to derive single-peak hyetographs and double-peak hyetographs.
The main objectives of this study are: (1) to derive the design hyetographs for a 60-min rainstorm with different return periods in Yangpu District in order to serve as the input to an urban runoff model; (2) to explore the influences of separating single-peak and double-peak rainstorms in deriving hyetographs; (3) to compare the methodology of combining mode hyetographs and the P&C method with the classic design hyetograph method. The methodology and procedure used to derive hyetographs in this study can also be a reference for other small cities and regions.

Study Area
Yangpu District in Shanghai, China, was selected as the study area. Shanghai is located in the Southeast of China. Yangpu District is in the central part of the city, and has an area of about 60 km 2 . The Huangpu River winds along the southeast side of the district. The subtropical monsoon climate results in an average annual precipitation of 1060 mm.

Process of Precipitation Data
Because Yangpu District is in the downtown area, there are no adequate precipitation data inside the district. Therefore, precipitation data from three nearby stations: Gaoqiao, Yangjing, and Jiangwan ( Figure 1) were used for areal rainfall estimation, frequency analysis, and the derivation of hyetographs. Due to the availability of the data, the annual maximum 60-min rainfalls from 1978 to 2015 in Gaoqiao and Yangjing were used for areal rainfall estimation and frequency analysis. The precipitation data from Jiangwan Station with 5 min intervals from 2001 to 2015 were used to derive the hyetographs. The duration for the hyetograph design was 60 min. The individual rainfall events were divided before designing the hyetograph. Different criteria for selecting inter-event times are reported in [23][24][25][26][27][28]. In this study, 120 min was selected as the inter-event time based on the regulations in China [29], by which the successive rainfall processes were divided into independent rainfall events. Rainfall events that exceed 16 mm in 60 min were selected (also based on the regulations in China) to derive the hyetographs, and in total 117 events were selected.
Water 2018, 10, x FOR PEER REVIEW 3 of 12

Process of Precipitation Data
Because Yangpu District is in the downtown area, there are no adequate precipitation data inside the district. Therefore, precipitation data from three nearby stations: Gaoqiao, Yangjing, and Jiangwan ( Figure 1) were used for areal rainfall estimation, frequency analysis, and the derivation of hyetographs. Due to the availability of the data, the annual maximum 60-min rainfalls from 1978 to 2015 in Gaoqiao and Yangjing were used for areal rainfall estimation and frequency analysis. The precipitation data from Jiangwan Station with 5 min intervals from 2001 to 2015 were used to derive the hyetographs. The duration for the hyetograph design was 60 min. The individual rainfall events were divided before designing the hyetograph. Different criteria for selecting inter-event times are reported in [23][24][25][26][27][28]. In this study, 120 min was selected as the inter-event time based on the regulations in China [29], by which the successive rainfall processes were divided into independent rainfall events. Rainfall events that exceed 16 mm in 60 min were selected (also based on the regulations in China) to derive the hyetographs, and in total 117 events were selected.

Overall Procedure
The procedure for deriving the 60-min design hyetographs is as follows: 1. Calculation of the areal rainfall in Yangpu District using the arithmetic average method. 2. Frequency analysis using Pearson-III distribution in order to derive the maximum 60-min areal rainfall with return periods of 2, 5, 10, and 20 years. 3. Classification of rainfall events into seven types of mode hyetographs using the fuzzy identification method. 4. Use the P&C method to derive single-peak hyetographs and double-peak hyetographs.
The duration of the design rainstorm was selected based on the following reasons: (1) According to the standard of rainstorm intensity formula and design rainstorm distribution in Shanghai [30], the duration of the design rainstorm should be 5-180 min for drainage system design; (2) The data in the study area are limited, and 60 min is the duration in which the data has the best quality; (3) 60 min is applicable to a large portion of urban drainage work [31].
The design rainstorm using the proposed methodology was compared to the classic design rainstorm to test its superiority. The standard of rainstorm intensity formula and design rainstorm

Overall Procedure
The procedure for deriving the 60-min design hyetographs is as follows: 1.
Calculation of the areal rainfall in Yangpu District using the arithmetic average method.

2.
Frequency analysis using Pearson-III distribution in order to derive the maximum 60-min areal rainfall with return periods of 2, 5, 10, and 20 years.

3.
Classification of rainfall events into seven types of mode hyetographs using the fuzzy identification method.

4.
Use the P&C method to derive single-peak hyetographs and double-peak hyetographs.
The duration of the design rainstorm was selected based on the following reasons: (1) According to the standard of rainstorm intensity formula and design rainstorm distribution in Shanghai [30], the duration of the design rainstorm should be 5-180 min for drainage system design; (2) The data in the study area are limited, and 60 min is the duration in which the data has the best quality; (3) 60 min is applicable to a large portion of urban drainage work [31].
The design rainstorm using the proposed methodology was compared to the classic design rainstorm to test its superiority. The standard of rainstorm intensity formula and design rainstorm distribution in Shanghai suggest the use of the Chicago rainstorm method for the design rainstorm. Therefore, the design rainstorm using the proposed methodology was compared with the one using Chicago rainstorm method.

Areal Rainfall Estimation
The arithmetic average method was used to calculate the areal rainfall due to the limited rainfall data in the region. The average value of the rainfall in all stations was used as the areal rainfall. It is expressed as: where x is the areal rainfall, x i is the rainfall in each station, and n is the number of stations.

Frequency Analysis
The maximum 60-min rainstorm amount can be represented by Pearson-III distribution [32]. It is expressed as: where α and β are calculated as: In this study, x, C s , and C v were estimated by the L-moment method [33]. They can be expressed as follows: If |τ 3 | < 1 3 , let z = 3πτ 2 3 , α can be expressed by: Thus x, C s , and C v can be expressed by: The concept and computation of l 1 , l 2 , l 3 , and τ 3 (τ 3 = l 2 /l 3 ) can be found in Hosking and Wallis [33]. The exceedance probability was estimated using the Weibull equation: where m is the rank number for the sequence in descending order, P m is the empirical probability for the mth term of the sequence, and n is the length of the sequence.

Fuzzy Identification to Classify Hyetograph
The 117 rainfall events were classified into seven types of mode hyetographs. Figure 2 shows the hyetographs of each type.
Water 2018, 10, x FOR PEER REVIEW 5 of 12

Fuzzy Identification to Classify Hyetograph
The 117 rainfall events were classified into seven types of mode hyetographs. Figure 2 shows the hyetographs of each type. Fuzzy identification was used for classification. The basic concept of fuzzy identification is to compare mode matrices with actual indices of rainfall events and classify the hyetograph into the one with the least differences. The steps for the method are as follows: Each rainfall event is divided into n periods, and the percentage of the rainfall in each period is expressed as: where i is the ith period of the rainfall process, Hi is the rainfall in the ith period, and H is the total rainfall. This group of xi is defined as the actual rainfall indices, denoted as vector X: The mode matrices of the seven hyetographs are defined in a similar way, denoted as Vk: where k is the kth type of hyetograph and vki is computed in the same way as xi.
The nearness degree ( ) of actual rainfall indices and mode matrices is computed as: The larger indicates that the rainfall event is closer to a kth type of hyetograph, and thus can be classified into it.
In this study, rainfall events were divided into six equal periods, and the mode matrices of the seven hyetographs are defined as: Fuzzy identification was used for classification. The basic concept of fuzzy identification is to compare mode matrices with actual indices of rainfall events and classify the hyetograph into the one with the least differences. The steps for the method are as follows: Each rainfall event is divided into n periods, and the percentage of the rainfall in each period is expressed as: where i is the ith period of the rainfall process, H i is the rainfall in the ith period, and H is the total rainfall. This group of x i is defined as the actual rainfall indices, denoted as vector X: The mode matrices of the seven hyetographs are defined in a similar way, denoted as V k : where k is the kth type of hyetograph and v ki is computed in the same way as x i . The nearness degree (σ k ) of actual rainfall indices and mode matrices is computed as: The larger σ k indicates that the rainfall event is closer to a kth type of hyetograph, and thus can be classified into it.
In this study, rainfall events were divided into six equal periods, and the mode matrices of the seven hyetographs are defined as: In the matrix, the seven rows indicate seven types of hyetographs. The six columns indicate the rainfall events divided into six periods. Each term in the matrix indicates that the percentage of the rainfall amount in the ith period of the event should be v ki .

Pilgrim and Cordery Method
The steps of t P&C method [15] are as follows:

1.
Each rainfall event is divided into n periods by certain time intervals, and the percentage of the rainfall in each period is computed.

2.
The periods are sorted in descending order according to the percentage of the rainfall, and for each period, the average rank of all the rainfall events is computed as the rank for this period.

3.
For all the rainfall events, the average percentage of the rainfall in the periods of the same rank is computed.

4.
For each period, the rank and the percentage of the rainfall are paired to obtain the temporal distribution of the rainfall process.

5.
The hyetograph can be obtained by multiplying the percentage of the rainfall in each period by the design rainstorm amount with a given return period.

The Chicago Rainstorm Method
The Chicago rainstorm method [13] is based on an intensity-duration-frequency (IDF) curve and the rainstorm intensity formula. According to the standard of rainstorm intensity formula and design rainstorm distribution in Shanghai, the rainfall intensity is calculated as: i(t a ) = where t a is the duration after the peak, t b is the duration before the peak, i(t a ) is the rainfall intensity after the peak (mm/min), i(t b ) is the rainfall intensity before the peak (mm/min), P is the return period, and r is the position coefficient of rainstorm peak (adopt 0.405 according to the standard of rainstorm intensity formula and design rainstorm distribution in Shanghai [30]).

Results and Discussion
The P-III distribution curve and empirical probability points are shown in Figure 3. The maximum value occurred in 1985 and amounted to 116 mm. It can be seen that the P-III distribution curve fits the empirical probability points well, except for the maximum value of the sequence, which seems to be abnormal. The exceedance probability computed using Equation (10) for the point was 1/39. However, one assumes the actual exceedance probability for it should be smaller than 1/39, although there is not concrete evidence for this assumption due to limited data and records in the region. Considering this aspect, the fitted P-III curve should be rational. The Nash-Sutcliffe efficiency coefficient (NSE) between observation and theoretical values was 0.94. The maximum 60-min areal rainfall values with different return periods computed from the fitted P-III distribution are shown in Table 1.   The classification of the 117 rainfall events is presented in Table 2. It can be seen that single-peak hyetographs make up most of the rainfall events, taking up about 60%. Double-peak hyetographs account for approximately 30% of the rainfall events, the rest (about 10%) are uniformly distributed hyetographs. Among single-peak hyetographs, the ones with peak falling at the end of the process (Type III) slightly outnumber the other two types. For double-peak hyetographs, the number of Type V hyetographs-whose peaks fall at the beginning and end-are slightly larger than the number of Type VI and Type VII. Although the number of single-peak hyetographs is two times that of doublepeak hyetographs, the latter still makes up a considerable portion of all the hyetographs. Therefore, it is worth exploring the influences of separating and combining these two kinds of hyetographs.   The classification of the 117 rainfall events is presented in Table 2. It can be seen that single-peak hyetographs make up most of the rainfall events, taking up about 60%. Double-peak hyetographs account for approximately 30% of the rainfall events, the rest (about 10%) are uniformly distributed hyetographs. Among single-peak hyetographs, the ones with peak falling at the end of the process (Type III) slightly outnumber the other two types. For double-peak hyetographs, the number of Type V hyetographs-whose peaks fall at the beginning and end-are slightly larger than the number of Type VI and Type VII. Although the number of single-peak hyetographs is two times that of double-peak hyetographs, the latter still makes up a considerable portion of all the hyetographs. Therefore, it is worth exploring the influences of separating and combining these two kinds of hyetographs.  (Figure 4), the peak falls in the 5th period, which is in the middle of the event. Most of the rainfall concentrates in the 3rd to 6th periods, which take up 69.69% of the total rainfall amount. The single-peak hyetograph seems to be quite smooth and rational. For the double-peak hyetograph of scenario 2 ( Figure 5), the peaks fall in the 2nd and 8th periods, which are at the beginning and end of the event. The rainfall concentrates in the 8th to 10th period, which account for 45.95% of the total rainfall amount. In general, the double-peak hyetograph is rational and smooth, except after the second peak, the rainfall in the 10th period exceeds that in the 9th period. However, the rainfall in the two periods does not have much difference in quantity, and thus it is acceptable. The general hyetograph of scenario 3 ( Figure 6) seems to have three peaks, in the 3rd, 6th, and 8th periods. The peaks are not as obvious, and the rainfall is not as concentrated when compared with other two scenarios. Therefore, it demonstrates that when combining single-peak and double-peak rainfall events, the derived hyetograph will lose important features, and thus is not rational in some senses. rainfall amount. In general, the double-peak hyetograph is rational and smooth, except after the second peak, the rainfall in the 10th period exceeds that in the 9th period. However, the rainfall in the two periods does not have much difference in quantity, and thus it is acceptable. The general hyetograph of scenario 3 ( Figure 6) seems to have three peaks, in the 3rd, 6th, and 8th periods. The peaks are not as obvious, and the rainfall is not as concentrated when compared with other two scenarios. Therefore, it demonstrates that when combining single-peak and double-peak rainfall events, the derived hyetograph will lose important features, and thus is not rational in some senses.   rainfall amount. In general, the double-peak hyetograph is rational and smooth, except after the second peak, the rainfall in the 10th period exceeds that in the 9th period. However, the rainfall in the two periods does not have much difference in quantity, and thus it is acceptable. The general hyetograph of scenario 3 ( Figure 6) seems to have three peaks, in the 3rd, 6th, and 8th periods. The peaks are not as obvious, and the rainfall is not as concentrated when compared with other two scenarios. Therefore, it demonstrates that when combining single-peak and double-peak rainfall events, the derived hyetograph will lose important features, and thus is not rational in some senses.    Percentage of rainfall(%) Periods Figure 6. Rainfall distribution of the design hyetograph for scenario 3: combining single-peak and double-peak rainstorms.
Generally speaking, single-peak rainstorms have more concentrated rainfall, and especially those peaks that appear in the early part of the event are more liable to cause urban flood issues. Therefore, more attention should be paid to single-peak hyetographs. Double-peak rainstorms are less prone to cause waterlogging issues. Nevertheless, if the rainfall is highly concentrated on one of the peaks, it may still cause flood issues. To explore the influences of different hyetographs on urban floods, one can compare the outputs of the waterlogging model when using different hyetographs as an input for the model.
The final hyetograph can be obtained by multiplying the percentage of the rainfall in each period by the design areal rainfall amount with a given return period. For example, a 20-year 60-min design hyetograph with a single peak is presented in Table 3, whereas the double-peak hyetograph is shown in Table 4. Table 3. Rainfall allocation of a single-peak rainstorm with a 60-min duration and a 20-year return period in Yangpu District.  Because the design rainstorm derived by the Chicago rainstorm method (denoted as Chicago hyetograph) can only be a single-peak rainstorm, it is only compared with the single-peak rainstorm hyetograph derived by the P&C method (denoted as P&C hyetograph). The 60-min design hyetographs with 20-year return periods are shown in Figure 7. The total rainfall amount of the Chicago hyetograph (76.46 mm) is smaller than that of the P&C method (82.51 mm), due to the different data used in both methods. The parameter of the intensity rainstorm formula in the standard of rainstorm intensity formula and design rainstorm distribution in Shanghai was estimated using rainfall data from the rain-gauge stations in the range of the whole city of Shanghai, while the P&C hyetograph only used rainfall data of two rain-gauge stations near Yangpu District. The P&C hyetograph is more typical for the study area in this aspect. Both of the peaks fall in the 5th period (20-25 min). In the rising part of the hyetograph, the P&C hyetograph has a higher rainfall amount (45.82 mm, 55.53%) than the Chicago hyetograph (36.39 mm, 47.59%). In general, the P&C hyetograph is similar to the Chicago hyetograph, which proves the rationality of separating single-peak and double-peak rainstorms while applying the P&C method. Furthermore, the P&C hyetograph has a higher rainfall amount before reaching peak, which is less favorable in terms of runoff model results. In addition, the Chicago rainstorm method can only derive single-peak hyetographs, while the P&C method can derive both single-peak and double-peak hyetographs. Therefore, the design rainstorm using the proposed methodology is superior to the classic Chicago rainstorm method. is similar to the Chicago hyetograph, which proves the rationality of separating single-peak and double-peak rainstorms while applying the P&C method. Furthermore, the P&C hyetograph has a higher rainfall amount before reaching peak, which is less favorable in terms of runoff model results.
In addition, the Chicago rainstorm method can only derive single-peak hyetographs, while the P&C method can derive both single-peak and double-peak hyetographs. Therefore, the design rainstorm using the proposed methodology is superior to the classic Chicago rainstorm method.

Conclusions and Recommendations
With the change of climate, rainstorms caused by extreme weather events are becoming increasingly frequent. The time distribution of hyetographs, as the input to various hydrological models, exerts a critical influence on the output of the model. In this study, 60-min design rainstorms with return periods of 2, 5, 10, 20 years were derived for the urban area of Yangpu District. The design areal rainfall of different return periods was calculated based on frequency analysis using the method of L-moments. The rainfall events were initially classified into seven types of mode hyetographs and then used to derive the 60-min duration hyetograph by using the P&C method with respect to three scenarios: single-peak, double-peak, and general hyetographs. The derived hyetograph using the P&C method was compared with the classic Chicago rainstorm method. Based on the results, the following conclusions can be drawn: • The result of the classification of the 117 rainstorm events indicates that single-peak hyetographs account for about 60% of these events, compared to double-peak hyetographs which account for about 30%. The rest are uniformly distributed hyetographs. • Both the single-peak and double-peak hyetographs are smooth and rational. The general hyetograph is irrational due to its relatively even and disordered distribution of rainfall amount.

Conclusions and Recommendations
With the change of climate, rainstorms caused by extreme weather events are becoming increasingly frequent. The time distribution of hyetographs, as the input to various hydrological models, exerts a critical influence on the output of the model. In this study, 60-min design rainstorms with return periods of 2, 5, 10, 20 years were derived for the urban area of Yangpu District. The design areal rainfall of different return periods was calculated based on frequency analysis using the method of L-moments. The rainfall events were initially classified into seven types of mode hyetographs and then used to derive the 60-min duration hyetograph by using the P&C method with respect to three scenarios: single-peak, double-peak, and general hyetographs. The derived hyetograph using the P&C method was compared with the classic Chicago rainstorm method. Based on the results, the following conclusions can be drawn:

•
The result of the classification of the 117 rainstorm events indicates that single-peak hyetographs account for about 60% of these events, compared to double-peak hyetographs which account for about 30%. The rest are uniformly distributed hyetographs. • Both the single-peak and double-peak hyetographs are smooth and rational. The general hyetograph is irrational due to its relatively even and disordered distribution of rainfall amount.

•
The design rainstorm using the proposed methodology in this study shows superiority to the classic Chicago rainstorm method. On one hand, it can derive both single-peak and double-peak hyetographs. On the other hand, it is more disadvantageous as input to runoff models.
In conclusion, rainstorms of different types have different distributions and thus have different influences on urban runoff models. It is recommended to derive hyetographs for each type of rainstorm separately, and then choose the most disadvantageous one as the input for waterlogging model (or other hydrological models).