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.     

Return period of bivariate distributed extreme hydrological events

 Extreme hydrological events are inevitable and stochastic in nature. Characterized by multiple properties, the multivariate distribution is a better approach to represent this complex phenomenon than the univariate frequency analysis. However, it requires considerably more data and more sophisticated mathematical analysis. Therefore, a bivariate distribution is the most common method for modeling these extreme events. The return periods for a bivariate distribution can be defined using either separate single random variables or two joint random variables. In the latter case, the return periods can be defined using one random variable equaling or exceeding a certain magnitude and/or another random variable equaling or exceeding another magnitude or the conditional return periods of one random variable given another random variable equaling or exceeding a certain magnitude. In this study, the bivariate extreme value distribution with the Gumbel marginal distributions is used to model extreme flood events characterized by flood volume and flood peak. The proposed methodology is applied to the recorded daily streamflow from Ichu of the Pachang River located in Southern Taiwan. The results show a good agreement between the theoretical models and observed flood data. The author wishes to thank the two anonymous reviewers for their constructive comments that improving the quality of this work.     

Trivariate generalized extreme value distribution in flood frequency analysis

Abstract

The multivariate extension of the logistic model with generalized extreme value (GEV) marginals is applied to provide a regional at-site flood estimate. The maximum likelihood estimators of the parameters were obtained numerically by using a multivariable constrained optimization algorithm. The asymptotic results were checked by distribution sampling techniques in order to establish whether or not those results can be utilized for small samples. A region in northern Mexico with 21 gauging stations was selected to apply the model. Results were compared with those obtained by the most popular univariate distributions, the bivariate approach of the logistic model and three regional methods: station-year, index flood and L-moments. These show that there is a reduction in the standard error of fit when estimating the parameters of the marginal distribution with the trivariate distribution instead of its univariate and bivariate counterpart, and differences between at-site and regional at-site design events can be significant as return period increases.     

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

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

The bivariate lognormal distribution to model a multivariate flood episode

Complex hydrological events such as floods always appear to be multivariate events that are characterized by a few correlated variables. A complete understanding of these events needs to investigate joint probabilistic behaviours of these correlated variables. The lognormal distribution is one of frequently selected candidates for flood‐frequency analysis. The multivariate lognormal distribution will serve as an important tool for analysing a multivariate flood episode. This article presents a procedure for using the bivariate lognormal distribution to describe the joint distributions of correlated flood peaks and volumes, and correlated flood volumes and durations. Joint distributions, conditional distributions, and the associated return periods of these random variables can be readily derived from their marginal distributions. The approach is verified using observed streamflow data from the Nord river basin, located in the Province of Quebec, Canada. The theoretical distributions show a good fit to observed ones. Copyright © 2000 John Wiley & Sons, Ltd.     

INTERNATIONAL HYDROLOGICAL DECADE / WORKING GROUP ON “THE INFLUENCE OF MAN ON THE HYDROLOGICAL CYCLE”

Abstract

Floods, as extreme hydrological phenomena, can be described by more than one correlated characteristic, such as peak, volume and duration. These characteristics should be jointly considered since they are generally not independent. For an ungauged site, univariate regional flood frequency analysis (FA) provides a limited assessment of flood events. A recent study proposed a procedure for regional FA in a multivariate framework. This procedure represents a multivariate version of the index-flood model and is based on copulas and multivariate quantiles. The performance of the proposed procedure was evaluated by simulation. However, the model was not tested on a real-world case study data. In the present paper, practical aspects are investigated jointly for flood peak (Q) and volume (V) of a dataset from the Côte-Nord region in the province of Quebec, Canada. The application of the proposed procedure requires the identification of the appropriate marginal distribution, the estimation of the index flood and the selection of an appropriate copula. The results of the case study show that the regional bivariate FA procedure performed well. This performance depends strongly on the performance of the two univariate models and, more specifically, the univariate model of Q. The results show also the impact of the homogeneity of the region on the performance of the univariate and bivariate models.

Editor D. Koutsoyiannis     

Copula‐based flood frequency (COFF) analysis at the confluences of river systems

Many civil infrastructures are located near the confluence of two streams, where they may be subject to inundation by high flows from either stream or both. These infrastructures, such as highway bridges, are designed to meet specified performance objectives for floods of a specified return period (e.g. the 100 year flood). Because the flooding of structures on one stream can be affected by high flows on the other stream, it is important to know the relationship between the coincident exceedence probabilities on the confluent stream pair in many hydrological engineering practices. Currently, the National Flood Frequency Program (NFF), which was developed by the US Geological Survey (USGS) and based on regional analysis, is probably the most popular model for ungauged site flood estimation and could be employed to estimate flood probabilities at the confluence points. The need for improved infrastructure design at such sites has motivated a renewed interest in the development of more rigorous joint probability distributions of the coincident flows. To accomplish this, a practical procedure is needed to determine the crucial bivariate distributions of design flows at stream confluences. In the past, the copula method provided a way to construct multivariate distribution functions. This paper aims to develop the Copula‐based Flood Frequency (COFF) method at the confluence points with any type of marginal distributions via the use of Archimedean copulas and dependent parameters. The practical implementation was assessed and tested against the standard NFF approach by a case study in Iowa's Des Moines River. Monte Carlo simulations proved the success of the generalized copula‐based joint distribution algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

Design event selection in bivariate hydrological frequency analysis

Abstract

In the bivariate analysis of hydrological events, such as rainfall storms or flood hydrographs, the choice of an appropriate return period for structure design leads to infinite combinations of values of the related random variables (e.g. peak and volume in the analysis of floods). These combinations are generally not equivalent, from a practical point of view. In this paper, a methodology is proposed to identify 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. Therefore, several combinations can be selected within the subset, taking into account the specific characteristic of the design problem, in order to evaluate the effects of different hydrological loads on a structure. The proposed method is applicable to any type of bivariate distribution, thus providing a simple but effective rule to narrow down the infinite possible choices for the hydrological design variables. In order to illustrate how the proposed methodology can be easily used in practice, it is applied to a study case in the context of bivariate flood frequency analysis.

Editor Z.W. Kundzewicz; Associate editor Sheng Yue

Citation Volpi, E. and Fiori, A., 2012. Design event selection in bivariate hydrological frequency analysis. Hydrological Sciences Journal, 57 (8), 1506–1515.     

Copula models for frequency analysis what can be learned from a Bayesian perspective?

Large spring floods in the Québec region exhibit correlated peakflow, duration and volume. Consequently, traditional univariate hydrological frequency analyses must be complemented by multivariate probabilistic assessment to provide a meaningful design flood level as requested in hydrological engineering (based on return period evaluation of a single quantity of interest). In this paper we study 47 years of a peak/volume dataset for the Romaine River with a parametric copula model. The margins are modeled with a normal or gamma distribution and the dependence is depicted through a parametric family of copulas (Arch 12 or Arch 14). Parameter joint inference and model selection are performed under the Bayesian paradigm. This approach enlightens specific features of interest for hydrological engineering: (i) cross correlation between margin parameters are stronger than expected , (ii) marginal distributions cannot be forgotten in the model selection process and (iii) special attention must be addressed to model validation as far as extreme values are of concern.     

基于二次重现期的多变量洪水风险评估

由于洪水是一种具有多个特征属性的随机事件,频率分析成为洪水风险评估的一种有效手段,多变量重现期与设计值的定义与计算则是洪水频率分析中的重点和难点.本文通过构造洪水历时、洪峰与洪量的联合分布,介绍了一种新的多变量重现期定义——二次重现期,并探讨了"或"重现期、"且"重现期和二次重现期对安全与危险域识别的差异性,以及在洪水风险管理与工程设计中的合理性与可靠性.传统的"或"和"且"多变量重现期对安全与危险域的识别存在局限性,利用Kendall函数定义的二次重现期则提供了更加合理的安全与风险域识别,避免了对安全事件与危险事件的错误判定,更有利于指导洪水风险的管理.在给定的二次重现期条件下,依据出现概率最大原则推算的历时、洪峰与洪量设计值组合可以满足工程设计以较低成本承受较大风险的追求,相比于单变量设计值,考虑了洪水多个属性联合特征的多变量设计值提供了更加全面和可靠的参考信息.     

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.     

Applicability of the Nagao–Kadoya bivariate exponential distribution for modeling two correlated exponentially distributed variates

 The open literature reveals several types of bivariate exponential distributions. Of them only the Nagao–Kadoya distribution (Nagao and Kadoya, 1970, 1971) has a general form with marginals that are standard exponential distributions and the correlation coefficient being 0≤ρ<1. On the basis of the principle that if a theoretical probability distribution can represent statistical properties of sample data, then the computed probabilities from the theoretical model should provide a good fit to observed ones, numerical experiments are executed to investigate the applicability of the Nagao–Kadoya bivariate exponential distribution for modeling the joint distribution of two correlated random variables with exponential marginals. Results indicate that this model is suitable for analyzing the joint distribution of two exponentially distributed variables. The procedure for the use of this model to represent the joint statistical properties of two correlated exponentially distributed variables is also presented.     

Statistical analysis of flash flood events for designing water harvesting systems in an extremely arid environment

A water harvesting system for research purposes has been established in the Lisan Peninsula of the Dead Sea in the middle of the Jordan Rift Valley, where no authorized guideline is available for designing water harvesting systems. Rainfall and runoff, which occurred as flash floods, were observed at the downstream end of a gorge with a 1.12 km² barren catchment area from October 2014 through July 2019. Due to the extremely arid environment, runoff from the catchment is ephemeral, and the flash flood events can be clearly distinguishable from each other. Thirteen flash flood events with a total runoff volume of more than 100 m³ were successfully recorded during the five rainy seasons. Pearson and Spearman correlations between duration, total rainfall depths at two points, total runoff volume, maximum runoff discharge, bulk runoff coefficient, total variation in runoff discharge and maximum variation in runoff discharge of each flash flood event were examined, revealing no straightforward relationship between rainfall and runoff. The performance of the conventional SCS runoff curve number method was also deficient in reproducing any rainfall–runoff relationship. Therefore, probability distribution fitting was performed for each random variable, focusing on the lognormal distribution with three parameters and the generalized extreme value distribution. The maximum goodness-of-fit estimation turns out to be a more rational and efficient method in obtaining the parameter values of those probability distributions rather than the standard maximum likelihood estimation, which has known disadvantages. Results support the design of the water harvesting system and provide quantitative information for designing and operating similar systems in the future.     

A new seasonal design flood method based on bivariate joint distribution of flood magnitude and date of occurrence

Abstract

Seasonal design floods which consider information on seasonal variation are very important for reservoir operation and management. The seasonal design flood method currently used in China is based on seasonal maximum (SM) samples and assumes that the seasonal design frequency is equal to the annual design frequency. Since the return period associated with annual maximum floods is taken as the standard in China, the current seasonal design flood cannot satisfy flood prevention standards. A new seasonal design flood method, which considers dates of flood occurrence and magnitudes of the peaks (runoff), was proposed and established based on copula function. The mixed von Mises distribution was selected as marginal distribution of flood occurrence dates. The Pearson Type III and exponential distributions were selected as the marginal distribution of flood magnitude for annual maximum flood series and peak-over-threshold samples, respectively. The proposed method was applied at the Geheyan Reservoir, China, and then compared with the currently used seasonal design flood methods. The case study results show that the proposed method can satisfy the flood prevention standard, and provide more information about the flood occurrence probabilities in each sub-season. The results of economic analysis show that the proposed design flood method can enhance the floodwater utilization rate and give economic benefits without lowering the annual flood protection standard.

Citation Chen, L., Guo, S. L., Yan, B. W., Liu, P. & Fang, B. (2010) A new seasonal design flood method based on bivariate joint distribution of flood magnitude and date of occurrence. Hydrol. Sci. J. 55(8), 1264–1280.     

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.     

Dam risk assessment based on univariate versus bivariate statistical approaches: a case study for Argentina

Abstract

Considering floods as multivariate events allows a better statistical representation of their complexity. In this work the relevance of multivariate analysis of floods for designing or assessing the safety of hydraulic structures is discussed. A flood event is characterized by its peak flow and volume. The dependence between the variables is modelled with a copula. One thousand random pairs of variables are transformed to hydrographs, applying the Beta distribution function. Synthetic floods are routed through a reservoir to assess the risk of overtopping a dam. The resulting maximum water levels are compared to estimations considering the peak flow and volume separately. The analysis is performed using daily flows observed in the River Agrio in Neuquén Province, Argentina, a catchment area of 7300 km². The bivariate approach results in higher maximum water level values. Therefore the multivariate approach should be preferred for the estimation of design variables.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

基于Copula函数的连河通江湖泊防洪安全设计:以洪泽湖为例

湖泊作为一种蓄水单元,尤其是大型过水性湖泊,是一种典型的平原型水库,在功能上与山谷型水库具有许多相似之处,但由于其特殊的地理地形构造,使得入湖洪水过程与入库洪水过程存在着较大的差异.在防洪安全设计研究中,山谷型水库关注的多是坝址洪水,即总的入库洪水过程,而对于湖泊来说,还需要关注各个分区的入湖洪水过程对湖区洪水演进的影响.针对大型过水性湖泊入湖洪水特征,本文采用Copula函数构造了多个联合分布函数,提出了一套基于总的入湖洪水过程推导各个分区入湖洪水过程置信区间的方法.以洪泽湖为应用实例,结果表明:1)在联合重现期已知的情况下,该方法能够确定总入湖洪量与洪峰的95%置信区间;2)该方法通过径流相关性分析对入湖河道合并聚类,形成分区入湖过程,既考虑了河道间天然的水文、水力联系,又避免了联合分布函数维度过高的问题;3)在总入湖洪量已知的情况下,该方法能够确定各分区入湖洪量分配95%置信区域.该方法具有较强的统计理论基础,拓展了多变量洪水频率分析技术在水利工程实际中的应用范围.     

Delineation of flood generating processes and their hydrological response

Knowledge about flood generating processes can be beneficial for numerous applications. Especially in the context of climate change impact assessment, daily patterns of meteorological and catchment state conditions leading to flood events (i.e., storylines) may be of value. Here, we propose an approach to identify storylines of flood generation using daily weather and snow cover observations. The approach is tested for and applied to a typical pre‐Alpine catchment in the period between 1961 and 2014. Five precipitation parameters were determined that describe temporal and spatial characteristics of the flood associated precipitation events. The catchment's snow coverage was derived using statistical relationships between a satellite‐derived snow cover climatology and station snow measurements. Moreover, (pre‐) event snow melt sums were estimated using a temperature‐index model. Based on the precipitation and catchment state parameters, 5 storylines were identified with a cluster analysis: These are (a) long duration, low intensity precipitation events with high precipitation depths, (b) long duration precipitation events with high precipitation depths and episodes of high intensities, (c) shorter duration events with high or (d) low precipitation intensity, respectively, and (e) rain‐on‐snow events. The event groups have distinct hydrological characteristics that can largely be explained by the storylines' respective properties. The long duration, high intensity storyline leads to the most adverse hydrological response, namely, a combination of high peak magnitudes, high volumes, and long durations of threshold exceedance. The results show that flood generating processes in mesoscale catchments can be distinguished on the basis of daily meteorological and catchment state parameters and that these process types can explain the hydrological flood properties in a qualitative way. Hydrological simulations of daily resolution can thus be analysed with respect to flood generating processes.