Changjiang Xu 1,2 and Jiabo Yin 1 and Shenglian Guo 1 and Zhangjun Liu 1 and Xingjun Hong 1
Academic Editor:Renata Archetti
1, State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
2, Bureau of Hydrology, Changjiang Water Resources Commission, Wuhan 430010, China
Received 18 October 2015; Revised 19 January 2016; Accepted 14 February 2016
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
A reservoir dam can be designed with the design flood hydrograph (DFH), which is a hydrograph adopted according to design standards to ensure the safety of a structure. DFH generally refers to the computed flood hydrograph at the dam site that may be estimated from precipitation or flow records, depending on engineering judgment after consideration of all the pertinent data, including the extent and reliability of the flow data and rainfall data [1]. The DFH is characterized by the joint behavior of several correlated random variables such as flood peak, volume, and duration. Knowledge of the flood peak or volumes alone is not sufficient to design the dam spillway; the entire flood hydrograph must be utilized. The quantiles of the flood peak and volumes corresponding to a particular design return period are necessary for deriving DFH. However, the conventional flood frequency analysis methods for deriving DFH recommended by many countries are based on the univariate distribution, mainly concentrated on the analysis of annual peak discharge or flood volume series without analyzing the inherent relationship between flood peak and volumes [2]. For example, the Chinese univariate flood frequency analysis method for deriving DFH implicitly assumed that flood peak and flood volumes are independent and amplified with the same frequency, which has shortcomings in practice [3, 4]. Over- or underestimation of risk would be produced by the univariate frequency analysis method because of the unclear relationship between the correlated hydrological variables [5, 6].
To solve this problem, many authors have used statistical methods to derive the DFH by examining the dependence of flood peak, flood volume, and its shape. Xiao et al. [2] proposed a multicharacteristic synthesis index method based on the MSI quantile and derived the DFH. Pramanik et al. [7] used the probability density functions to fit the shape of hydrographs and subsequently to develop design flood hydrographs for various return periods. Mediero et al. [8] estimated flood volume for a given peak discharge by fitting a regional log-log regression equation over the observed pairs, and a Monte Carlo experiment was conducted to generate an ensemble of DHFs that maintain the statistical properties of marginal distributions of the peaks, volumes, and durations. Serinaldi and Grimaldi [9] studied the advantages and shortcomings of using simple distribution functions with finite support (namely, beta and generalized standard two-sided power distributions) to represent and synthesize direct runoff hydrographs. In their study, the dependence among peak discharge, volume, and duration was explored on a few flood events selected by a recursive digital filter algorithm and an over threshold approach. Dominguez and Arganis [10] introduced and used the IINGEN method to estimate the design floods for Malpaso Dam, and the method makes it possible for considering flood peak, volume, and shape simultaneously. Óscar et al. [11] proposed a method to select the time base of the DFH using spectral analysis method and defined the number of days necessary for taking into account the complete dynamics of the hydrological floods.
Recently, the copula functions have been increasingly applied in exploring the inherent correlation between flood peak and volumes for deriving DFH. Copula functions allow for more flexibility in the marginal distributions and the dependence between peaks and volumes. Chowdhary et al. [12] discussed the identification of suitable copulas for bivariate frequency analysis of flood peak and flood volume data, and Gumbel-Hougaard copula of Archimedean family was found to be the most suitable dependence model for flood peak and volume. Bacova and Halmova [13] explored the analysis and statistical evaluation of the joint probability of the occurrence of peak discharge and volume; they also calculated joint return periods and conditional return periods for the hydrological pair. There have been numerous studies on the degree of the dependence between peaks and volumes that is needed to estimate the multivariate quantiles (e.g., [5, 6]) and the choice of the copula function (e.g., [12, 14]). Gaal et al. [15] studied the peak-volume relationships for the region of Austria, which featured a diverse spectrum of hydrological flood processes and discussed whether flood peak-volume relationships can be typified by comparative hydrology. Sraj et al. [16] compared different bivariate copulas from three families and analyzed the bivariate flood frequency analysis in the Litija station in the Sava River. The previous researches [12-16] derived the multivariate quantile curve under one joint return period, and infinite combinations on the isoline or isosurface could be selected to evaluate the effects of different hydrological loads on the structure [17]. However, one has to select a quantiles combination out of events which all share the same return period for deriving DFH [14]. How to solve this problem has attracted the attention of many researchers (e.g., [6, 14]). Li et al. [6] proposed a copula identical frequency method and used it to estimate bivariate joint flood quantiles in Three Gorges reservoir. Graler et al. [14] gave a critical and practical review focusing on synthetic design hydrograph estimation. Their approaches were based on regression analysis, bivariate conditional distributions, bivariate joint distributions, and Kendall distribution function, highlighting theoretic and practical issues of multivariate frequency analysis. Chebana and Ouarda [18] proposed the decomposition of the level curve into a naïve part (tail) and the proper part (central); they assumed that the naïve part was composed of two segments starting at the end of each extremity of the proper part. Salvadori et al. [19] introduced two basic design realizations, that is, componentwise excess design realization and most likely design realization. Volpi and Fiori [20] identified a subset of the critical combinations set that includes a fixed and arbitrarily chosen percentage in probability of the events, on the basis of their probability of occurrence.
In the multivariate domain for deriving DFH, the return period must be defined in terms of the acceptable risk for the reservoir. Different return periods estimated by copula function have been developed for the case of a multivariate flood frequency analysis. Eight types of possible joint events were presented by Yue and Rasmussen [21] and Shiau [22], in which the OR, AND, conditional cases are of greatest interest in hydrological applications. Recently, another definition of the multivariate return period was given by Salvadori and De Michele [23], called Kendall return period. This return period was introduced to identify in a multivariate context a univariate critical threshold. Until now, only a very limited number of studies actually applied this kind of return period (e.g., [24]). Volpi and Fiori [25] proposed a general, structure-based return period for the design and risk assessment of hydrological structures in a bivariate environment. Their work draws attention of practitioners to the importance of considering the structure in hydraulic design and/or risk assessment problems in a multivariate environment, thus advising against the uncritical use of design event-based approaches which neglect the interplay between the structure and the hydrological loads acting on it. More information about the return periods could be found in Shiau [22] and Volpi and Fiori [25]. Until now, a clear assessment of the use of return periods is still a matter of debates [26, 27].
In this paper, we study the peak-volume relationships by copula functions for Danjiangkou reservoir in the Hanjiang basin extensively. Specifically, we investigate the inherent relationship on the basis of their probability of occurrence. Different copula functions from Archimedean family were applied and compared. Then the conditional most likely composition (CMLC) and conditional expectation composition (CEC) methods were proposed, inspired by the idea of Graler et al. [14] and Salvadori et al. [19]. The main aim of this study is to provide a simple tool for the identification of the relationship between flood peak and volumes and derive DFH.
To achieve these objectives, the paper is organized as follows: introduction of the Chinese method of deriving DFH is illustrated in Section 2. The methodologies used in this study are presented in Section 3. In Section 4, results and discussions of case study are shown, followed by the Conclusions.
2. Chinese Method of Deriving DFH
Pearson type III (P3) distribution is recommended by Ministry of Water Resources of China for deriving DFH [2, 3]. For a given sample series [figure omitted; refer to PDF] , the probability density function (PDF) of P3 distribution is expressed as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are the shape, scale, and location parameters, respectively, and [figure omitted; refer to PDF] is the gamma function.
The L-moment (LM) method is preferred for parameter estimation because of its robust properties in the presence of usually small or large values (outliers) [28, 29]. Univariate flood frequency analysis model (P3/LM), that is, P3 distribution and L-moment method, is used to estimate flood quantile for given return periods.
The peak and volume amplitude (PVA) method implicitly assumed that flood peak and flood volumes are amplified with the same frequency (return period). For example, the values of peak discharge and flood volumes are individually estimated by each marginal distribution under given return periods firstly. Then the DFH is constructed by multiplying each discharge ordinate of the typical flood hydrograph (TFH) by an amplifier.
The amplifier is either the ratio of the design peak with a given return period to the TFH peak or the ratio of the design volume with a given return period to the TFH volume. Suppose that the characteristics of annual maximum flood hydrograph consist of flood peak [figure omitted; refer to PDF] , 1-day maximum flood volume [figure omitted; refer to PDF] , 3-day maximum flood volume [figure omitted; refer to PDF] , 7-day maximum flood volume [figure omitted; refer to PDF] , and 15-day maximum flood volume [figure omitted; refer to PDF] , the quantiles with design return period [figure omitted; refer to PDF] are estimated by P3/LM and denoted by [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively. The corresponding characteristics of selected TFH are denoted by [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively. These amplifiers are calculated as follows.
Amplifier [figure omitted; refer to PDF] for flood peak discharge is as follows: [figure omitted; refer to PDF]
Amplifier [figure omitted; refer to PDF] for the 1-day maximum flood volume (excluding flood peak) is as follows: [figure omitted; refer to PDF]
Amplifier [figure omitted; refer to PDF] for the 3-day maximum flood volume except for the 1-day maximum flood volume is as follows: [figure omitted; refer to PDF]
Amplifier [figure omitted; refer to PDF] for the 7-day maximum flood volume except for the 3-day maximum flood volume is as follows: [figure omitted; refer to PDF]
Amplifier [figure omitted; refer to PDF] for the 15-day maximum flood volume except for the 7-day maximum flood volume is as follows: [figure omitted; refer to PDF]
The dividing portion between two contiguous segments of the constructed DFH is not continuous and should be modified by hand smooth treatment in engineering practice [2, 3].
3. Derivation of DFH Based on Copula Function
3.1. Joint Distribution Based on Copula Function
Copula function is a function that describes and models the dependence structure between random variables, independent of the marginal distributions involved (e.g., [26, 28, 30]). Let [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] ) be the cumulative distribution function (CDF) of [figure omitted; refer to PDF] . The objective is to determine the multivariate distribution, denoted as [figure omitted; refer to PDF] or simply [figure omitted; refer to PDF] . Thus, the multivariate probability distribution [figure omitted; refer to PDF] is expressed in terms of its marginal and the associated dependence function as Sklar's theorem: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , called the copula function, is uniquely determined whenever [figure omitted; refer to PDF] are continuous and captures the essential features of the dependence among the random variables.
3.2. Conditional Most Likely Composition (CMLC) Method
The P3 distribution is selected as marginal distributions of peak discharge ( [figure omitted; refer to PDF] ) and flood volumes ( [figure omitted; refer to PDF] ), which are denoted as [figure omitted; refer to PDF] , [figure omitted; refer to PDF] . Then the corresponding density functions are denoted as [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , in which [figure omitted; refer to PDF] represented the number of flood volumes needed for estimating DFH. For a given design standard (return period), it is necessary to find an appropriate combination of flood quantiles since different combinations would result in different DFH [18, 20]. How to select appropriate combinations (i.e., flood quantiles composition method) is extremely important in practice. The conventional univariate hydrological frequency analysis method for deriving DFH implicitly assumed that flood peak and volumes are independent and amplified with the same frequency, which could not explore the inherent correlation of hydrological variables [2]. Under given peak discharge, the combinations differ in terms of their probability of occurrence.
Based on copula function, the joint distribution function [figure omitted; refer to PDF] can be expressed as [30] [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is copula function; [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , are marginal distribution of peak discharge and flood volumes, respectively.
The corresponding conditional probability distribution is [figure omitted; refer to PDF]
The density function of conditional probability distribution is [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is density function of copula.
Maximize [figure omitted; refer to PDF] ; the corresponding values of [figure omitted; refer to PDF] are conditional most likely composition (CMLC).
[figure omitted; refer to PDF] is a single function of [figure omitted; refer to PDF] ; its derivation to [figure omitted; refer to PDF] is given by [figure omitted; refer to PDF] in which [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] is the derivative function of [figure omitted; refer to PDF] .
To maximize [figure omitted; refer to PDF] , the following equation should be satisfied: [figure omitted; refer to PDF]
If [figure omitted; refer to PDF] is following P3 distribution as (1), then the following equation should be satisfied: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are shape, scale, and location parameters of marginal distribution function [figure omitted; refer to PDF] , respectively.
Substitute (13) into (12) and simplify it as follows: [figure omitted; refer to PDF]
Equation (14) is the general formula of the CMLC method. If the peak discharge and flood volumes are independent, then the value of [figure omitted; refer to PDF] would be 0, and the value of [figure omitted; refer to PDF] would be 1. Substituting the values of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] into (14), it could be simplified as follows: [figure omitted; refer to PDF] . However, the independent assumption is only a special case in practice.
The CMLC method is an approach to describe the compositions of flood peak and volumes by using the conditional density function to measure the occurrence likelihood of flood events. It is the mode of the conditional distribution [figure omitted; refer to PDF] . Based on the principle of maximizing conditional joint probability density function, (14) can be solved by numerical computation methods, such as harmonic mean Newton's method [31].
3.3. Conditional Expectation Composition (CEC) Method
The CEC method is also proposed to estimate the volumes under given hydrological variables, which is another way to analyze the inherent dependence between correlated flood peak and volumes. If the variables [figure omitted; refer to PDF] are known, then the conditional expectation [figure omitted; refer to PDF] is used to estimate flood volume; that is, [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is [figure omitted; refer to PDF] reverse function.
Equation (15) is the general formula of conditional expectation composition (CEC) method, which can be solved by Gauss-Legendre quadrature rule [32].
3.4. Confidence Interval of Conditional Probability Distribution
If [figure omitted; refer to PDF] are known, let [figure omitted; refer to PDF]
For a given significance level [figure omitted; refer to PDF] , the lower and upper limits [figure omitted; refer to PDF] and [figure omitted; refer to PDF] identifying the confidence interval of flood volume [figure omitted; refer to PDF] are derived by [figure omitted; refer to PDF]
Specifically, if [figure omitted; refer to PDF] , then the 5% and 95% quantiles compose the 90% confidence intervals of conditional probability distribution.
3.5. Evaluation Indexes
To assess and compare different methods, the root mean square error (RMSE) and relative error (Bias) were selected as evaluation indexes: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are observed and estimated values and [figure omitted; refer to PDF] is the length of sample series. The smaller these indexes are, the better the fitting or estimating methods are.
4. Case Study
4.1. Danjiangkou Reservoir
The Hanjiang basin is the largest tributary of the Yangtze River; it passes through the provinces of Shanxi and Hubei of China and merges into the Yangtze River at Wuhan city. The river's length is 1570 km and the basin area is 159,000 km2 with a subtropical monsoon climate. The regional precipitation and flow are characterized by high seasonal variability, with 75% of rainfall occurring from June to October, where the rare rainstorms in the early summer and long-lasting rains in the autumn often result in great floods.
The Danjiangkou reservoir located in the middle reach of the basin (Figure 1) is the source of water for the middle route of the South-North Water Diversion Project in China. The available water resources in Hanjiang basin and the impact of water diversion have been discussed by many authors [33]. Since it is a multipurpose reservoir, there is serious conflict between flood prevention, water supply, and hydropower generation. The normal pool level is 170 m and corresponding storage capacity is 29.05 billion m3 . The annual maximum peak discharges and flood volumes of the inflow of Danjiangkou reservoir are available with a systematic record of 61 years (1954-2014). The flood data series was provided by the Bureau of Danjiangkou Water Resources Management. This data series has been restored and checked for possible long-term trends and abrupt changes, and the test results demonstrate that the data set is stationary and represented. The sample statistics of the flood peak and volume series were listed in Table 1.
Table 1: Statistics of flood peak and volume series at Danjiangkou reservoir.
| Flood peak (103 m3 /s) | 3-day flood volume (108 m3 ) | 7-day flood volume (108 m3 ) | 15-day flood volume (108 m3 ) |
[figure omitted; refer to PDF] Min, max] | [figure omitted; refer to PDF] 2.3, 34.3] | [figure omitted; refer to PDF] 4.61, 66.43] | [figure omitted; refer to PDF] 10.04, 93.12] | [figure omitted; refer to PDF] 20.96, 134.40] |
Median | 15.50 | 18.45 | 37.67 | 67.95 |
Mean | 15.77 | 28.03 | 46.07 | 71.82 |
Variance | 8.32 | 15.27 | 24.07 | 35.17 |
Coefficient of skewness | 0.384 | 0.628 | 0.551 | 0.481 |
Figure 1: Location of Danjiangkou reservoir in the Hanjiang basin.
[figure omitted; refer to PDF]
4.2. Flood Quantiles Estimated by P3/LM
Univariate hydrological frequency analysis model (P3/LM) was used to estimate flood quantile for given return periods. The estimated parameters of the P3 distribution for peak discharge and flood volumes series were listed in Table 2. A Chi-Square [figure omitted; refer to PDF] goodness-of-fit test was performed to test the assumption, [figure omitted; refer to PDF] , that the flood magnitudes follow the P3 distribution. Table 2 shows that the assumption could not be rejected at the 5% significance level.
Table 2: Estimated parameters of the P3 distribution and hypothesis test.
Variable | Sample statistic values | P3 parameter | [figure omitted; refer to PDF] statistics | [figure omitted; refer to PDF] | ||||
Mean | [figure omitted; refer to PDF] - [figure omitted; refer to PDF] | [figure omitted; refer to PDF] - [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | |||
[figure omitted; refer to PDF] (m3 /s) | 15768 | 0.54 | 0.58 | 11.8906 | 0.0004 | -13593 | 0.28 | 3.84 |
[figure omitted; refer to PDF] (108 m3 ) | 11.94 | 0.55 | 0.64 | 9.7656 | 0.4759 | -8.5819 | 0.76 | 3.84 |
[figure omitted; refer to PDF] (108 m3 ) | 28.03 | 0.56 | 0.82 | 5.9488 | 0.1554 | -10.2549 | 0.88 | 3.84 |
[figure omitted; refer to PDF] (108 m3 ) | 46.07 | 0.54 | 0.84 | 5.6689 | 0.0957 | -13.1629 | 1.03 | 3.84 |
[figure omitted; refer to PDF] (108 m3 ) | 71.82 | 0.5 | 0.71 | 7.9349 | 0.0784 | -29.3349 | 1.29 | 3.84 |
For given different return periods ( [figure omitted; refer to PDF] = 5, 10, 20, 50, 100, 200, 500, and 1000 years), the corresponding quantiles of flood peaks and volumes were estimated based on the univariate P3/LM model, and the results were listed in Table 3.
Table 3: Comparison of design floods estimated by P3/LM and CMLC methods.
Return period (year) | [figure omitted; refer to PDF] (m3 /s) | [figure omitted; refer to PDF] (108 m3 ) | [figure omitted; refer to PDF] (108 m3 ) | ||||
P3/LM | P3/LM | CMLC | Bias (%) | P3/LM | CMLC | Bias (%) | |
5 | 22625 | 65.4 | 64.7 | -1.07 | 100.2 | 98.1 | -0.02 |
10 | 27135 | 79.3 | 78.6 | -0.88 | 119.7 | 116.6 | -0.03 |
20 | 31122 | 92.0 | 91.4 | -0.65 | 137.2 | 134.2 | -0.02 |
50 | 35901 | 107.5 | 107.0 | -0.47 | 158.4 | 154.6 | -0.02 |
100 | 39261 | 118.7 | 118.2 | -0.42 | 174.5 | 170.7 | -0.02 |
200 | 42461 | 129.4 | 128.0 | -1.09 | 189.9 | 184.3 | -0.03 |
500 | 46498 | 143.0 | 140.9 | -1.47 | 208.2 | 201.8 | -0.03 |
1000 | 49437 | 156.1 | 148.3 | -4.50 | 224.6 | 211.7 | -0.06 |
The empirical and P3 distribution frequency curves of flood peaks and flood volumes were drawn in Figure 2, on which the Weibull plotting position formula [34] was used. It is shown that theoretical values can fit the observed values very well.
Figure 2: The fit of P3 probability curves with observed flood peak and flood volumes.
(a) Flood peak
[figure omitted; refer to PDF]
(b) Three-day flood volume
[figure omitted; refer to PDF]
(c) Seven-day flood volume
[figure omitted; refer to PDF]
(d) 15-day flood volume
[figure omitted; refer to PDF]
4.3. Joint Distribution of Flood Peak and Volumes Based on Copulas
Different families of copula functions have been proposed and described by Nelsen [30]. The Archimedean copula family is more desirable for hydrologic analyses, because most of the family can easily be constructed and some of the family can be applied whether the correlation among the hydrological variables is positive or negative [16, 26]. The copula functions (Table 4) were applied to construct the bivariate joint distributions between flood peak ( [figure omitted; refer to PDF] ) and flood volumes ( [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] ). The values of Kendall correlation coefficient ( [figure omitted; refer to PDF] ) are 0.86, 0.80, 0.75, and 0.73 for the four pairs [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , respectively. The parameters of bivariate copulas were estimated by the method of the inversion of Kendall's [figure omitted; refer to PDF] . Based on the results of Kendall's [figure omitted; refer to PDF] , the AMH copula could not be used because it is only suitable when Kendall's coefficient is below 0.34 [30].
Table 4: Bivariate copula functions and the relationship between their parameter and Kendall correlation coefficient.
Copula | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
Gumbel-Hougaard (G-H) | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
Clayton | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
Ali-Mikhail-Haq (AMH) | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
Frank | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
The G-H, Clayton, and Frank copula functions were compared with different statistical tests. The [figure omitted; refer to PDF] values were calculated based on the parametric bootstrap or multipliers [35] procedure with 10 000 runs. The results of the Cramer-von Mises test ( [figure omitted; refer to PDF] ) and Kolmogorov-Smirnov test ( [figure omitted; refer to PDF] ) were listed in Table 5.
Table 5: Estimated parameters of copulas and results of statistical test.
Copula | G-H | Clayton | Frank | ||||||
Pair | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] value) | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] value) | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] value) | [figure omitted; refer to PDF] |
[figure omitted; refer to PDF] - [figure omitted; refer to PDF] | 6.91 | 0.010 (0.39) | 0.023 | 11.81 | 0.022 (0.21) | 0.037 | 25.87 | 0.019 (0.22) | 0.035 |
[figure omitted; refer to PDF] - [figure omitted; refer to PDF] | 5.08 | 0.012 (0.38) | 0.027 | 8.17 | 0.026 (0.17) | 0.040 | 18.53 | 0.022 (0.19) | 0.030 |
[figure omitted; refer to PDF] - [figure omitted; refer to PDF] | 4.15 | 0.011 (0.35) | 0.026 | 6.30 | 0.024 (0.18) | 0.039 | 14.75 | 0.020 (0.21) | 0.032 |
[figure omitted; refer to PDF] - [figure omitted; refer to PDF] | 3.77 | 0.013 (0.36) | 0.028 | 5.55 | 0.021 (0.20) | 0.035 | 13.21 | 0.017 (0.27) | 0.031 |
The comparison between simulated (sample size 10 000) and observed values for three chosen copula functions was plotted in Figure 3. On the basis of statistical tests and graphical goodness-of-fit tests, the G-H copula is much suitable for constructing the bivariate joint distribution of peak discharge and flood volumes. This finding is in accordance with Li et al. [6] and Sraj et al. [16].
Figure 3: Comparison between the observed and simulated values for three chosen copula functions.
(a) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (G-H)
[figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Clayton)
[figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Frank)
[figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (G-H)
[figure omitted; refer to PDF]
(e) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Clayton)
[figure omitted; refer to PDF]
(f) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Frank)
[figure omitted; refer to PDF]
(g) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (G-H)
[figure omitted; refer to PDF]
(h) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Clayton)
[figure omitted; refer to PDF]
(i) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Frank)
[figure omitted; refer to PDF]
(j) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (G-H)
[figure omitted; refer to PDF]
(k) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Clayton)
[figure omitted; refer to PDF]
(l) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] (Frank)
[figure omitted; refer to PDF]
Therefore, the G-H copula was used to model the dependence between the flood peak and flood volumes in this study. The linear correlation between empirical and theoretical frequency values of bivariate G-H copulas for flood peak and volumes was plotted in Figure 4, on which the linear correlation coefficients are 0.9994, 0.9991, 0.9991, and 0.9991, respectively. The comparison of empirical plots and theoretic frequency curves of bivariate G-H copulas was also plotted in Figure 5, which shows that the theoretical frequency line (or curve) fit empirical values very well with high correlation coefficients.
Figure 4: Linear correlation between empirical and theoretical frequency values of bivariate G-H copulas for flood peak and volumes.
(a) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
Figure 5: Comparison of the empirical plots and theoretical frequency curves of bivariate G-H copula for flood peak and volumes.
(a) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
The G-H, Clayton, and Frank copulas from Archimedean family were also used to establish trivariate symmetric and asymmetric joint distributions of [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF] , respectively. Table 6 shows selected trivariate copula functions, where parameter of symmetric copulas is [figure omitted; refer to PDF] and asymmetric copulas are [figure omitted; refer to PDF] and [figure omitted; refer to PDF] ( [figure omitted; refer to PDF] ). Parameters of trivariate copula functions were estimated by the maximum likelihood method [26]. To assess and compare six trivariate copula functions, the Kolmogorov-Smirnov test value [figure omitted; refer to PDF] and Akaike information criterion (AIC) [29] were selected and the results were listed in Table 7. It is shown that AIC and [figure omitted; refer to PDF] values of asymmetrical G-H copulas are less than those of other copulas. Therefore, the asymmetrical G-H copulas are selected to establish 3-dimensional joint distribution of flood peak and volumes in this study. Similar to bivariate joint distribution, the linear correlation relationship and comparison of the empirical plots and theoretical frequency line (or curve) of asymmetrical trivariate G-H copulas were plotted in Figures 6 and 7, which show that the theoretical frequency line (or curve) fit empirical values very well with high correlation coefficients.
Table 6: Trivariate symmetric and asymmetric Archimedean copula functions.
Copula functions | [figure omitted; refer to PDF] or [figure omitted; refer to PDF] |
Symmetric Gumbel-Hougaard (S-G-H) | [figure omitted; refer to PDF] |
| |
Symmetric Clayton (S-Clayton) | [figure omitted; refer to PDF] |
| |
Symmetric Frank (S-Frank) | [figure omitted; refer to PDF] |
| |
Asymmetric Gumbel-Hougaard-M6 (A-G-H) | [figure omitted; refer to PDF] |
| |
Asymmetric Clayton-M4 (A-Clayton) | [figure omitted; refer to PDF] |
| |
Asymmetric Frank-M3 (A-Frank) | [figure omitted; refer to PDF] |
Table 7: Estimated parameters of trivariate copulas and results of statistical test.
Pair | [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF] | [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF] | [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF] | ||||||
Copula | [figure omitted; refer to PDF] | AIC | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | AIC | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | AIC | [figure omitted; refer to PDF] |
S-G-H | 4.79 | 204.3 | 0.042 | 3.82 | 198.4 | 0.045 | 3.63 | 218.4 | 0.046 |
S-Clayton | 7.25 | 210.3 | 0.047 | 5.84 | 206.4 | 0.049 | 5.24 | 225.7 | 0.050 |
S-Frank | 15.34 | 208.9 | 0.044 | 12.54 | 204.8 | 0.047 | 11.24 | 220.4 | 0.049 |
A-G-H | [figure omitted; refer to PDF] 4.56, 6.45] | 200.5 | 0.035 | [figure omitted; refer to PDF] 3.77, 4.03] | 193.7 | 0.037 | [figure omitted; refer to PDF] 3.51, 3.69] | 206.5 | 0.040 |
A-Clayton | [figure omitted; refer to PDF] 7.01, 7.89] | 204.6 | 0.043 | [figure omitted; refer to PDF] 5.72, 6.24] | 197.2 | 0.046 | [figure omitted; refer to PDF] 5.12, 5.65] | 213.9 | 0.047 |
A-Frank | [figure omitted; refer to PDF] 15.03, 16.02] | 202.5 | 0.039 | [figure omitted; refer to PDF] 11.89, 12.67] | 200.4 | 0.041 | [figure omitted; refer to PDF] 11.09, 12.34] | 209.6 | 0.044 |
Figure 6: Linear correlation between the empirical and theoretical frequency values of trivariate asymmetric G-H copulas for flood peaks and volumes.
(a) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
Figure 7: Comparison of the empirical plots and theoretical frequency curves of trivariate asymmetric G-H copulas for flood peaks and volumes.
(a) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF] - [figure omitted; refer to PDF] - [figure omitted; refer to PDF]
[figure omitted; refer to PDF]
4.4. CMLC and CEC Estimation and Assessment
Flood volumes were estimated by the CMLC, CEC methods and compared with P3/LM model. For bivariate distribution, the flood volumes ( [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] ) based on peak discharge were estimated by (14) and the results were denoted as [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively. Similarly, the CEC method was used to estimate [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] by (15) for given peak discharge. As the P3/LM model is based on the assumption that volumes and peak discharge were in the same frequency, flood volumes ( [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] ) were estimated for given return period.
The RMSE and Bias were selected as evaluation indexes to assess the performance of different methods. Table 8 shows that the CMLC method has smallest RMSE and Bias while CEC method has the largest RMSE and Bias values. For a given significance level [figure omitted; refer to PDF] , then 90% confidence intervals of flood volumes were evaluated by (17). Figure 8 shows 7-day and 15-day flood volumes as well as confidence intervals (CI) estimated by CMLC and CEC methods. Figure 8 shows that most of observations are within confidence intervals. It could be concluded that CMLC method is reasonable and performs best in bivariate distribution.
Table 8: Comparison of flood volumes estimated by univariate and bivariate distributions.
Flood volumes | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | ||||
Index | RMSE | Bias (%) | RMSE | Bias (%) | RMSE | Bias (%) | RMSE | Bias (%) |
| ||||||||
CMLC | 0.1580 | 1.14 | 0.1833 | 1.41 | 0.2070 | 1.59 | 0.2299 | 2.53 |
CEC | 0.1680 | 2.78 | 0.2025 | 4.08 | 0.2205 | 4.70 | 0.2503 | 5.70 |
P3/LM | 0.1620 | 2.21 | 0.1878 | 2.49 | 0.2074 | 2.67 | 0.2359 | 4.21 |
Figure 8: Estimated volumes and confidence intervals for given peak discharges.
(a) Seven-day flood volume
[figure omitted; refer to PDF]
(b) 15-day flood volume
[figure omitted; refer to PDF]
To compare the performance of different methods in trivariate distribution, 15-day volumes were estimated by CMLC and CEC method, respectively, based on peak discharges and 7-day volumes. [figure omitted; refer to PDF] ) and [figure omitted; refer to PDF] were estimated by CMLC and CEC methods, respectively. Table 9 shows that CMLC method performs best with the smallest RMSE and Bias values, while CEC method has the largest values of indexes and univariate P3/LM model performs reasonably well.
Table 9: Comparison of flood volumes estimated by univariate and trivariate distributions.
Method | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | |||
Index | RMSE | Bias (%) | RMSE | Bias (%) | RMSE | Bias (%) |
| ||||||
3-dimension CMLC | 0.1525 | 0.87 | 0.1501 | -0.58 | 0.1518 | -0.33 |
3-dimension CEC | 0.1976 | 4.58 | 0.2741 | 4.71 | 0.2807 | 6.69 |
P3/LM | 0.1833 | 1.41 | 0.2070 | 1.59 | 0.2299 | 2.53 |
4.5. Derivation of Design Flood Hydrograph
All of these results demonstrate that the CMLC method performs best for analyzing the inherent relationship between flood peak and volumes. The proposed method was used to derive the DFH in Danjiangkou reservoir. There are so many definitions of return period, such as AND, OR, and Kendall return periods [23, 24]. The conditional return period [15, 27] was selected just as an example to demonstrate the proposed method to derive DFH. In Danjiangkou reservoir, the peak discharge, 7-day and 15-day volumes were always indicated as the control hydrological variables for deriving DFH. For different return periods (5, 10, 20, 50, 100, 200, 500, and 1000 years), the peak discharge was estimated by marginal distribution firstly. The second step is to estimate 7-day volumes under the calculated peak discharge by bivariate CMLC method. The third step is to estimate 15-day volumes under the calculated peak discharge and 7-day volumes by trivariate CMLC method. The design flood volumes estimated by CMLC and Chinese method were listed in Table 3. It is shown that 7-day and 15-day flood volumes estimated by CMLC method are less than those of P3/LM method for different return periods.
After getting the values of flood peak and volumes under each return period, the DFH could be derived by amplifying the typical flood hydrographs (TFH). The observed annual maximum flood hydrographs of 1975 and 1983 that have high peaks and large volumes with posterior-peak shape at the Danjiangkou reservoir were selected as TFH, respectively. The amplifiers of flood peak and volumes were calculated by (2) to (6), respectively. The DFH is derived by multiplying discharge ordinate of the TFH by corresponding amplification factor. As an illustrative example, only 1000-year DFH derived by the proposed CLMC and conventional Chinese methods were plotted in Figure 9 based on 1975 and 1983 TFH. It is shown that the flood quantiles estimated by the Chinese method are larger than that by the proposed CLMC method.
Figure 9: 1000-year DFH derived by the CMLC and Chinese methods based on 1975 and 1983 TFH.
(a) 1975 TFH
[figure omitted; refer to PDF]
(b) 1983 TFH
[figure omitted; refer to PDF]
5. Conclusions
The derivation of design flood hydrograph is very important in multivariate flood frequency analysis frameworks. Unique design floods under given standard are necessary for engineering practice. The simple regression method may be preferred due to its ease for application. To avoid the arbitrariness of design floods, the proposed conditional methods which have statistical basis might be one of the possible options for deriving flood quantiles and DFH. In this study, the general formulae of CMLC and CEC methods were proposed and developed to consider the relationship between the flood peaks and volumes. The proposed methods were applied in the Danjiangkou reservoir and compared with conventional Chinese univariate flood frequency analysis method. The following conclusions were drawn from this study.
(1) The Chinese method for deriving DFH implicitly assumed that flood peak and flood volumes are independent and amplified with the same frequency, which might produce over- or underestimation of risk.
(2) The Gumbel-Hougaard copula is more suitable for describing the correlation between flood peak and volumes. The results of case study in Danjiangkou reservoir demonstrate that the CMLC method performs best while the CEC method performs worst.
(3) The proposed CMLC method not only has strong statistical basis, but also has rational results with unique composition scheme and will provide a new approach for derivation of design flood hydrograph.
Acknowledgments
This study was supported by the Chinese National Nature Research Foundation (Grant nos. 51539009, 51579183).
[1] ICOLD (International Congress on Large Dams), "Selection of design flood," Bulletin Report , no. 82, ICOLD, Paris, France, 1992.
[2] Y. Xiao, S. L. Guo, P. Liu, B. Yan, L. Chen, "Design flood hydrograph based on multicharacteristic synthesis index method," Journal of Hydrologic Engineering , vol. 14, no. 12, pp. 1359-1364, 2009.
[3] MWR (Ministry of Water Resources) 2006 (Chinese) Regulation for Calculating Design Flood of Water Resources and Hydropower Projects , Water Resources & Hydropower Press, Beijing, China
[4] P. Liu, L. Li, S. Guo, L. Xiong, W. Zhang, J. Zhang, C. Xu, "Optimal design of seasonal flood limited water levels and its application for the Three Gorges Reservoir," Journal of Hydrology , vol. 527, pp. 1045-1053, 2015.
[5] L. Chen, V. P. Singh, G. Shenglian, Z. Hao, T. Li, "Flood coincidence risk analysis using multivariate copula functions," Journal of Hydrologic Engineering , vol. 17, no. 6, pp. 742-755, 2012.
[6] T. Y. Li, S. L. Guo, L. Chen, J. Guo, "Bivariate flood frequency analysis with historical information based on copula," Journal of Hydrologic Engineering , vol. 18, no. 8, pp. 1018-1030, 2013.
[7] N. Pramanik, R. K. Panda, D. Sen, "Development of design flood hydrographs using probability density functions," Hydrological Processes , vol. 24, no. 4, pp. 415-428, 2010.
[8] L. Mediero, A. Jimenez-Álvarez, L. Garrote, "Design flood hydrographs from the relationship between flood peak and volume," Hydrology and Earth System Sciences , vol. 14, no. 12, pp. 2495-2505, 2010.
[9] F. Serinaldi, S. Grimaldi, "Synthetic design hydrographs based on distribution functions with finite support," Journal of Hydrologic Engineering , vol. 16, no. 5, pp. 434-446, 2011.
[10] M. R. Dominguez, J. M. L. Arganis, "Validation of methods to estimate design discharge flow rates for dam spillways with large regulating capacity," Hydrological Sciences Journal , vol. 57, no. 3, pp. 460-478, 2012.
[11] A. Óscar, D. M. Ramon, A. J. Maritza, "Estimate of design hydrographs for the Angostura Dam, Sonora, using statistical and spectral methods," Water Resources Management , vol. 29, no. 11, pp. 4021-4043, 2015.
[12] H. Chowdhary, L. A. Escobar, V. P. Singh, "Identification of suitable copulas for bivariate frequency analysis of flood peak and flood volume data," Hydrology Research , vol. 42, no. 2-3, pp. 193-216, 2011.
[13] M. V. Bacova, D. Halmova, "Joint modeling of flood peak discharges, volume and duration: a case study of the Danube River in Bratislava," Journal of Hydrology and Hydromechanics , vol. 62, no. 3, pp. 186-196, 2014.
[14] B. Graler, M. J. VandenBerg, S. Vandenberghe, A. Petroselli, S. Grimaldi, B. De Baets, N. E. C. Verhoest, "Multivariate return periods in hydrology: a critical and practical review focusing on synthetic design hydrograph estimation," Hydrology and Earth System Sciences , vol. 17, no. 4, pp. 1281-1296, 2013.
[15] L. Gaal, J. Szolgay, S. Kohnova, K. Hlavcova, J. Parajka, A. Viglione, R. Merz, G. Blöschl, "Dependence between flood peaks and volumes: a case study on climate and hydrological controls," Hydrological Sciences Journal , vol. 60, no. 6, pp. 968-984, 2015.
[16] M. Sraj, N. Bezak, M. Brilly, "Bivariate flood frequency analysis using the copula function: a case study of the Litija station on the Sava River," Hydrological Processes , vol. 29, no. 2, pp. 225-238, 2015.
[17] D. Nijssen, A. Schumann, M. Pahlow, B. Klein, "Planning of technical flood retention measures in large river basins under consideration of imprecise probabilities of multivariate hydrological loads," Natural Hazards and Earth System Science , vol. 9, no. 4, pp. 1349-1363, 2009.
[18] F. Chebana, T. B. Ouarda, "Multivariate quantiles in hydrological frequency analysis," Environmetrics , vol. 22, no. 1, pp. 63-78, 2011.
[19] G. Salvadori, C. De Michele, F. Durante, "Multivariate design via Copulas," Hydrology and Earth System Sciences Discussions , vol. 8, no. 3, pp. 5523-5558, 2011.
[20] E. Volpi, A. Fiori, "Design event selection in bivariate hydrological frequency analysis," Hydrological Sciences Journal , vol. 57, no. 8, pp. 1506-1515, 2012.
[21] S. Yue, P. Rasmussen, "Bivariate frequency analysis: discussion of some useful concepts in hydrological application," Hydrological Processes , vol. 16, no. 14, pp. 2881-2898, 2002.
[22] J. T. Shiau, "Return period of bivariate distributed extreme hydrological events," Stochastic Environmental Research and Risk Assessment , vol. 17, no. 1-2, pp. 42-57, 2003.
[23] G. Salvadori, C. De Michele, "Frequency analysis via copulas: theoretical aspects and applications to hydrological events," Water Resources Research , vol. 40, no. 12, pp. 1-17, 2004.
[24] S. Vandenberghe, N. E. C. Verhoest, C. Onof, B. De Baets, "A comparative copula-based bivariate frequency analysis of observed and simulated storm events: a case study on Bartlett-Lewis modeled rainfall," Water Resources Research , vol. 47, no. 7, 2011.
[25] E. Volpi, A. Fiori, "Hydraulic structures subject to bivariate hydrological loads: return period, design, and risk assessment," Water Resources Research , vol. 50, no. 2, pp. 885-897, 2014.
[26] G. Salvadori, C. De Michele, "Multivariate real-time assessment of droughts via copula-based multi-site Hazard Trajectories and Fans," Journal of Hydrology , vol. 526, pp. 101-115, 2015.
[27] S. Yue, T. B. M. J. Ouarda, B. Bobee, P. Legendre, P. Bruneau, "The Gumbel mixed model for flood frequency analysis," Journal of Hydrology , vol. 226, no. 1-2, pp. 88-100, 1999.
[28] Q. Zhang, T. Qi, V. P. Singh, Y. D. Chen, M. Xiao, "Regional frequency analysis of droughts in China: a multivariate perspective," Water Resources Management , vol. 29, no. 6, pp. 1767-1787, 2015.
[29] L. Chen, S. L. Guo, B. W. Yan, "A new seasonal design flood method based on bivariate joint distribution of flood magnitude and date of occurrence," Hydrological Sciences Journal , vol. 7, no. 11, pp. 1264-1280, 2015.
[30] R. B. Nelsen An Introduction to Copulas , of Springer Series in Statistics, Springer, New York, NY, USA, 2006., 2nd.
[31] A. Y. Özban, "Some new variants of Newton's method," Applied Mathematics Letters , vol. 17, no. 6, pp. 677-682, 2004.
[32] E. Babolian, M. J. Masjed, M. R. Eslahchi, "On numerical improvement of Gauss-Legendre quadrature rules," Applied Mathematics and Computation , vol. 160, no. 3, pp. 779-789, 2005.
[33] S. L. Guo, J. X. Wang, L. H. Xiong, A. Ying, D. Li, "A macro-scale and semi-distributed monthly water balance model to predict climate change impacts in China," Journal of Hydrology , vol. 268, no. 1-4, pp. 1-15, 2002.
[34] S. L. Guo, C. Cunnane, "Evaluation of the usefulness of historical and palaeological floods in quantile estimation," Journal of Hydrology , vol. 129, no. 1-4, pp. 245-262, 1991.
[35] I. Kojadinovic, J. Yan, M. Holmes, "Fast large-sample goodness-of-fit tests for copulas," Statistica Sinica , vol. 21, no. 2, pp. 841-871, 2011.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2016 Changjiang Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
Design flood hydrograph (DFH) for a dam is the flood of suitable probability and magnitude adopted to ensure safety of the dam in accordance with appropriate design standards. Estimated quantiles of peak discharge and flood volumes are necessary for deriving the DFH, which are mutually correlated and need to be described by multivariate analysis methods. The joint probability distributions of peak discharge and flood volumes were established using copula functions. Then the general formulae of conditional most likely composition (CMLC) and conditional expectation composition (CEC) methods that consider the inherent relationship between flood peak and volumes were derived for estimating DFH. The Danjiangkou reservoir in Hanjiang basin was selected as a case study. The design values of flood volumes and 90% confidence intervals with different peak discharges were estimated by the proposed methods. The performance of CMLC and CEC methods was also compared with conventional flood frequency analysis, and the results show that CMLC method performs best for both bivariate and trivariate distributions which has the smallest relative error and root mean square error. The proposed CMLC method has strong statistical basis with unique design flood composition scheme and provides an alternative way for deriving DFH.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer