Bivariate frequency analysis: discussion of some useful concepts in hydrological application

Basic concepts such as conditional probability distributions, conditional return periods, and joint return periods are important to understand and to interpret multivariate hydrological events such as floods and storms. However, these concepts are not well documented in the open literature. This paper assembles and clarifies these concepts, and illustrates their practical utility. Relationships between joint return periods and univariate return periods are also derived. These concepts and relationships are demonstrated by applying a bivariate extreme value distribution to represent the joint distribution of flood peak and volume from an actual basin. Copyright © 2002 John Wiley & Sons, Ltd.     

A bivariate gamma distribution for use in multivariate flood frequency analysis

A gamma distribution is one of the most frequently selected distribution types for hydrological frequency analysis. The bivariate gamma distribution with gamma marginals may be useful for analysing multivariate hydrological events. This study investigates the applicability of a bivariate gamma model with five parameters for describing the joint probability behavior of multivariate flood events. The parameters are proposed to be estimated from the marginal distributions by the method of moments. The joint distribution, the conditional distribution, and the associated return periods are derived from marginals. The usefulness of the model is demonstrated by representing the joint probabilistic behaviour between correlated flood peak and flood volume and between correlated flood volume and flood duration in the Madawask River basin in the province of Quebec, Canada. Copyright © 2001 John Wiley & Sons, Ltd.     

An advanced data collection procedure in bivariate drought frequency analysis

In recent years, the bivariate frequency analysis of drought duration and severity using independent drought events and copula functions has been under extensive application. Meanwhile, emphasis on the procedure of independent drought data collection leads to the omission of the actual potential of short-term extreme droughts within a long-term independent drought. However, a long-term individual continuous drought as an Unconnected Drought Runs can be considered as a combination of short-term Connected Drought Runs. Thus, an advanced and new procedure of data collection in the bivariate drought characteristics analysis has been developed in this study. The results indicated a high relative advantage of the new proposed procedure in analysing bivariate drought characteristics (i.e., drought duration and severity frequency analysis). This advantage has been reflected in the more appropriate determination of the best copula and significant reduction in the uncertainty of bivariate drought frequency analysis.     

Spatial mapping of drought in Zambia using regional frequency analysis

Regional frequency analysis based on L-moments was applied to assess the spatial extent of meteorological droughts in tandem with their return periods in Zambia. Weather station monthly rainfall data were screened to form homogeneous sub-regions-, validated by a homogeneity criterion and fitted by a generalized extreme value distribution using goodness-of-fit test statistics. Predictor equations at regional scale for L-moment ratios and mean annual precipitation were developed to generate spatial maps of meteorological drought recurrences. The 80% of normal rainfall level and two thresholds of 60% and 70% were synonymous with moderate and severe droughts, respectively. Droughts were more severe in the south than in the north of Zambia. The return periods for severe and moderate droughts showed an overlapping pattern in their occurrence at many locations, indicating that in certain years droughts can affect the entire country. The extreme south of Zambia is the most prone to drought.     

Data‐based analysis of bivariate copula tail dependence for drought duration and severity

In recent decades, copula functions have been applied in bivariate drought duration and severity frequency analysis. Among several potential copulas, Clayton has been mostly used in drought analysis. In this research, we studied the influence of the tail shape of various copula functions (i.e. Gumbel, Frank, Clayton and Gaussian) on drought bivariate frequency analysis. The appropriateness of Clayton copula for the characterization of drought characteristics is also investigated. Drought data are extracted from standardized precipitation index time series for four stations in Canada (La Tuque and Grande Prairie) and Iran (Anzali and Zahedan). Both duration and severity data sets are positively skewed. Different marginal distributions were first fitted to drought duration and severity data. The gamma and exponential distributions were selected for drought duration and severity, respectively, according to the positive skewness and Kolmogorov–Smirnov test. The results of copula modelling show that the Clayton copula function is not an appropriate choice for the used data sets in the current study and does not give more drought risk information than an independent model for which the duration and severity dependence is not significant. The reason is that the dependence of two variables in the upper tail of Clayton copula is very weak and similar to the independent case, whereas the observed data in the transformed domain of cumulative density function show high association in the upper tail. Instead, the Frank and Gumbel copula functions show better performance than Clayton function for drought bivariate frequency analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

Performance of flood frequency pooling analysis in a low CV context

Abstract

Pooling of flood data is widely used to provide a framework to estimate design floods by the Index Flood method. Design flood estimation with this approach involves derivation of a growth curve which shows the relationship between X_T and the return period T, where X_T ?=?Q_T /Q_I and Q_I is the index flood at the site of interest. An implicit assumption with the Index Flood procedure of pooling analysis is that the X_T –T relationship is the same at all sites in a homogeneous pooling group, although this assumption would generally be violated to some extent in practical cases, i.e. some degree of heterogeneity exists. In fact, in only some cases is the homogeneity criterion effectively satisfied for Irish conditions. In this paper, the performance of the index-flood pooling analysis is assessed in the Irish low CV (coefficient of variation) hydrology context considering that heterogeneity is taken into account. It is found that the performance of the pooling method is satisfactory provided there are at least 350 station years of data included. Also it is found that, in a highly heterogeneous group, it is more desirable to have many sites with short record lengths than a smaller number of sites with long record lengths. Increased heterogeneity decreases the advantage of pooling group-based estimation over at-site estimation. Only a heterogeneity measure (H1) less than 4.0 can render the pooled estimation preferable to that obtained for at-site estimation for the estimation of 100-year flood. In moderately to highly heterogeneous regions it is preferable to conduct at-site analysis for the estimation of 100-year flood if the record length at the site concerned exceeds 50.

Editor Z.W. Kundzewicz; Associate editor A. Carsteanu

Citation Das, S. and Cunnane, C., 2012. Performance of flood frequency pooling analysis in a low CV context. Hydrological Sciences Journal, 57 (3), 433–444.     

Flood frequency analysis of Turkey using L‐moments method

The index flood procedure coupled with the L‐moments method is applied to the annual flood peaks data taken at all stream‐gauging stations in Turkey having at least 15‐year‐long records. First, screening of the data is done based on the discordancy measure (D_i) in terms of the L‐moments. Homogeneity of the total geographical area of Turkey is tested using the L‐moments based heterogeneity measure, H, computed on 500 simulations generated using the four parameter Kappa distribution. The L‐moments analysis of the recorded annual flood peaks data at 543 gauged sites indicates that Turkey as a whole is hydrologically heterogeneous, and 45 of 543 gauged sites are discordant which are discarded from further analyses. The catchment areas of these 543 sites vary from 9·9 to 75121 km² and their mean annual peak floods vary from 1·72 to 3739·5 m³ s^?1. The probability distributions used in the analyses, whose parameters are computed by the L‐moments method are the general extreme values (GEV), generalized logistic (GLO), generalized normal (GNO), Pearson type III (PE3), generalized Pareto (GPA), and five‐parameter Wakeby (WAK). Based on the L‐moment ratio diagrams and the |Z_dist|‐statistic criteria, the GEV distribution is identified as the robust distribution for the study area (498 gauged sites). Hence, for estimation of flood magnitudes of various return periods in Turkey, a regional flood frequency relationship is developed using the GEV distribution. Next, the quantiles computed at all of 543 gauged sites by the GEV and the Wakeby distributions are compared with the observed values of the same probability based on two criteria, mean absolute relative error and determination coefficient. Results of these comparisons indicate that both distributions of GEV and Wakeby, whose parameters are computed by the L‐moments method, are adequate in predicting quantile estimates. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulation of probability distributions commonly used in hydrological frequency analysis

Random variable simulation has been applied to many applications in hydrological modelling, flood risk analysis, environmental impact assessment, etc. However, computer codes for simulation of distributions commonly used in hydrological frequency analysis are not available in most software libraries. This paper presents a frequency‐factor‐based method for random number generation of five distributions (normal, log–normal, extreme‐value type I, Pearson type III and log‐Pearson type III) commonly used in hydrological frequency analysis. The proposed method is shown to produce random numbers of desired distributions through three means of validation: (1) graphical comparison of cumulative distribution functions (CDFs) and empirical CDFs derived from generated data; (2) properties of estimated parameters; (3) type I error of goodness‐of‐fit test. An advantage of the method is that it does not require CDF inversion, and frequency factors of the five commonly used distributions involves only the standard normal deviate. Copyright © 2006 John Wiley & Sons, Ltd.     

A comparison of three methods for regional flood frequency analysis in Saudi Arabia

The objective of the study was to compare the relative accuracy of three methodologies of regional flood frequency analysis in areas of limited flood records. Thirty two drainage basins of different characteristics, located mainly in the southwest region of Saudi Arabia, were selected for the study. In the first methodology, region curves were developed and used together with the mean annual flood, estimated from the characteristics of drainage basin, to estimate flood flows at a location in the basin. The second methodology was to fit probability distribution functions to annual maximum rainfall intensity in a drainage basin. The best fitted probability function was used together with common peak flow models to estimate the annual maximum flood flows in the basin. In the third methodology, duration reduction curves were developed and used together with the average flood flow in a basin to estimate the peak flood flows in the basin. The results obtained from each methodology were compared to the flood records of the selected stations using three statistical measures of goodness-of-fit. The first methodology was found best in a case of having short length of record at a drainage basin. The second methodology produced satisfactory results. Thus, it is recommended in areas where data are not sufficient and/or reliable to utilise the first methodology.     

Regional flood frequency analysis using Bayesian generalized least squares: a comparison between quantile and parameter regression techniques

Regression‐based regional flood frequency analysis (RFFA) methods are widely adopted in hydrology. This paper compares two regression‐based RFFA methods using a Bayesian generalized least squares (GLS) modelling framework; the two are quantile regression technique (QRT) and parameter regression technique (PRT). In this study, the QRT focuses on the development of prediction equations for a flood quantile in the range of 2 to 100 years average recurrence intervals (ARI), while the PRT develops prediction equations for the first three moments of the log Pearson Type 3 (LP3) distribution, which are the mean, standard deviation and skew of the logarithms of the annual maximum flows; these regional parameters are then used to fit the LP3 distribution to estimate the desired flood quantiles at a given site. It has been shown that using a method similar to stepwise regression and by employing a number of statistics such as the model error variance, average variance of prediction, Bayesian information criterion and Akaike information criterion, the best set of explanatory variables in the GLS regression can be identified. In this study, a range of statistics and diagnostic plots have been adopted to evaluate the regression models. The method has been applied to 53 catchments in Tasmania, Australia. It has been found that catchment area and design rainfall intensity are the most important explanatory variables in predicting flood quantiles using the QRT. For the PRT, a total of four explanatory variables were adopted for predicting the mean, standard deviation and skew. The developed regression models satisfy the underlying model assumptions quite well; of importance, no outlier sites are detected in the plots of the regression diagnostics of the adopted regression equations. Based on ‘one‐at‐a‐time cross validation’ and a number of evaluation statistics, it has been found that for Tasmania the QRT provides more accurate flood quantile estimates for the higher ARIs while the PRT provides relatively better estimates for the smaller ARIs. The RFFA techniques presented here can easily be adapted to other Australian states and countries to derive more accurate regional flood predictions. Copyright © 2011 John Wiley & Sons, Ltd.     

Applicability of Wakeby distribution in flood frequency analysis: a case study for eastern Australia

Parametric method of flood frequency analysis (FFA) involves fitting of a probability distribution to the observed flood data at the site of interest. When record length at a given site is relatively longer and flood data exhibits skewness, a distribution having more than three parameters is often used in FFA such as log‐Pearson type 3 distribution. This paper examines the suitability of a five‐parameter Wakeby distribution for the annual maximum flood data in eastern Australia. We adopt a Monte Carlo simulation technique to select an appropriate plotting position formula and to derive a probability plot correlation coefficient (PPCC) test statistic for Wakeby distribution. The Weibull plotting position formula has been found to be the most appropriate for the Wakeby distribution. Regression equations for the PPCC tests statistics associated with the Wakeby distribution for different levels of significance have been derived. Furthermore, a power study to estimate the rejection rate associated with the derived PPCC test statistics has been undertaken. Finally, an application using annual maximum flood series data from 91 catchments in eastern Australia has been presented. Results show that the developed regression equations can be used with a high degree of confidence to test whether the Wakeby distribution fits the annual maximum flood series data at a given station. The methodology developed in this paper can be adapted to other probability distributions and to other study areas. Copyright © 2014 John Wiley & Sons, Ltd.     

Implications of heterogeneous flood‐frequency distributions on traditional stream‐discharge prediction techniques

Traditional flood‐frequency analysis involves the assumption of homogeneity of the flood distribution. However, floods are often generated by heterogeneous distributions composed of a mixture of two or more populations. Differences between the populations may be the result of a number of factors, including seasonal variations in the flood‐producing mechanisms, changes in weather patterns resulting from low‐frequency climate shifts and/or El Niño/La Nina oscillations, changes in channel routing owing to the dominance of within‐channel or floodplain flow, and basin variability resulting from changes in antecedent soil moisture. Not recognizing these physical processes in conventional flood‐frequency analysis probably is the main reason that many frequency distributions do not provide an acceptable fit to flood data. In this paper, we use long‐term hydroclimatic records from the Gila River basin of south‐east and central Arizona in the USA to explore the extent and significance of mixed populations. First, we discuss the probable causes of heterogeneity in the frequency distribution of annual flood and present evidence of its occurrence. Second, we investigate the implications of using various popular homogeneous distributions for predicting peak flows for basins that exhibit mixed population characteristics. Third, we demonstrate how alternative frequency models that explicitly account for floods generated by a mixture of two or more populations are both hydrologically and statistically more appropriate. We illustrate how the selection of the most plausible distribution for flood‐frequency analysis also should be based on hydrological reasoning as opposed to the sole application of the traditional statistical goodness‐of‐fit tests. Copyright © 2002 John Wiley & Sons, Ltd.     

Validation of methods to estimate design discharge flow rates for dam spillways with large regulating capacity

Abstract

The estimation and review of discharge flow rates in hydraulic works is a fundamental problem in water management. In the case of dams with large regulating capacity, in order to estimate return periods of discharge flow rates from the spillways, it becomes necessary to consider both peak flow and volume of the incoming floods. In this paper, the results of the validation for several methods of assessing design floods for spillways of dams with a large flood control capacity are presented; the validation is performed by comparing the maximum outflows (or the maximum levels reached in the reservoir) obtained from the routing of the design floods with those obtained from the routing of the historical annual maximum floods. The basin of Malpaso Dam, Mexico, is used as the case study.

Editor D. Koutsoyiannis

Citation Domínguez, M.R. and Arganis, J.M.L., 2012. Validation of methods to estimate design discharge flow rates for dam spillways with large regulating capacity. Hydrological Sciences Journal, 57 (3), 460–478.     

Asymmetric copula in multivariate flood frequency analysis

The univariate flood frequency analysis is widely used in hydrological studies. Often only flood peak or flood volume is statistically analyzed. For a more complete analysis the three main characteristics of a flood event i.e. peak, volume and duration are required. To fully understand these variables and their relationships, a multivariate statistical approach is necessary. The main aim of this paper is to define the trivariate probability density and cumulative distribution functions. When the joint distribution is known, it is possible to define the bivariate distribution of volume and duration conditioned on the peak discharge. Consequently volume–duration pairs, statistically linked to peak values, become available. The authors build trivariate joint distribution of flood event variables using the fully nested or asymmetric Archimedean copula functions. They describe properties of this copula class and perform extensive simulations to highlight differences with the well-known symmetric Archimedean copulas. They apply asymmetric distributions to observed flood data and compare the results those obtained using distributions built with symmetric copula and the standard Gumbel Logistic model.     

基于IDL的地磁数据频谱分析系统设计

根据地磁台站数据在地震预报研究中的应用特点,研究了地磁数据获取、处理和存储的规范和流程,并根据规范和流程,利用可视化交互数据语言(IDL)和DataMiner接口实现了一个地磁台站数据频谱分析处理系统.在系统中,采用先处理存储后读写显示的方案,降低了系统的耦合度,改善了显示效果,提高了用户的体验度.     

水文干旱多变量联合设计及水库影响评估

基于东江流域博罗站月径流数据,采用游程理论提取水文干旱事件.选用Meta-Gaussian Copula函数,统计模拟水文干旱指标的多变量联合分布.采用Kendall联合重现期和最大可能权函数,设计给定联合超越重现期的水文干旱指标组合值,并定量评估水库径流调节作用对水文干旱多变量联合特征的影响.结果表明:东江流域水文干旱历时、强度和峰值的统计优选分布均为韦布尔分布.干旱指标之间具有较高的正相关性,Meta-Gaussian Copula能够很好地模拟水文干旱指标两变量和三变量联合分布.基于任意两个变量联合设计和三变量联合设计,干旱指标设计组合值位于同频位置附近,且同一个干旱指标设计值在不同变量组合之间差别较小.水库径流调节作用对于缓解东江流域水文干旱效果明显,同一组干旱指标的多变量联合超越重现期在水库影响下明显变大.联合超越重现期越小,水库对联合设计值的影响程度越大.根据目前水库运行模式,若要满足河道内最小管理流量目标,联合超越重现期10 a一遇的干旱历时、强度和峰值依然达到了约3.89~4.04月、7.20~7.97亿m3和2.99~3.12亿m3.     

Storm‐coverage effect on dynamic flood‐frequency analysis: empirical data analysis

In this study, a dynamic flood‐frequency analysis model considering the storm coverage effect is proposed and applied to six sub‐basins in the Pyungchang River basin, Korea. The model proposed is composed of the rectangular pulse Poisson process model for rainfall, the Soil Conservation Service curve number method for infiltration and the geomorphoclimatic instantaneous unit hydrograph for runoff estimation. Also, the model developed by Marco and Valdes is adopted for quantifying the storm‐coverage characteristics. By comparing the results from the same model with and without the storm‐coverage effect consideration, we could quantify the storm‐coverage effect on the flood‐frequency analysis. As a result of that, we found the storm‐coverage effect was so significant that overestimation of the design flood was unavoidable without its consideration. This also becomes more serious for larger basins where the probability of complete storm coverage is quite low. However, for smaller basins, the limited number of rain gauges is found to hamper the proper quantification of the storm‐coverage characteristics. Provided with a relationship curve between the basin size and the storm coverage (as in this study), this problem could be overcome with an acceptable accuracy level. Copyright © 2003 John Wiley & Sons, Ltd.     

Configuration of hydrological processes in a steep hillslope through causality analysis of soil moisture measurements

The vertical and lateral profiles of temporal variations in soil moisture are important for understanding the hydrological process along hillside transects. In this study, relationships among measured soil moistures were explored to configure the hydrological contributions of different flowpaths. All the measured soil moistures included a common stochastic structure because rainfall, the hydrometeological driver, is the main factor that determines the soil moisture response feature, and the infiltration process through the topsoil at a shallow depth is also common in all measured soil moisture histories. Therefore, the relationships between the measured series are also affected by both rainfall and topsoil infiltration. The common stochastic structure of the soil moisture series was removed via a prewhitening procedure. A systematic analysis procedure is presented to delineate the exclusive causal relationships among multiple soil moisture measurements. A monitoring system based on multiplexed time domain reflectometry was used to obtain soil moisture time series along two transects on a steep hillslope during the rainy season. The application of the proposed method for monitoring points in two adjacent locations provided 8, 12, 14, and 13, 16, 22 causal relationships for vertical, lateral in parallel, and diagonal directions, respectively, along the two transects. The point‐based contributions of the internal flowpath can be evaluated as the correlation is normalized in the context of inflow and outflow. The hydrological processes in the soil layer, vertical flow, lateral flow, downslope recharge, and return flow were quantified, and the relative importance of each hydrological component was determined to improve our understanding of the hydrological processes along the two transects of the study area. Copyright © 2011 John Wiley & Sons, Ltd.