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.     

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.     

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.     

Copula‐based joint probability function for PGA and CAV: a case study from Taiwan

This study aims to develop a joint probability function of peak ground acceleration (PGA) and cumulative absolute velocity (CAV) for the strong ground motion data from Taiwan. First, a total of 40,385 earthquake time histories are collected from the Taiwan Strong Motion Instrumentation Program. Then, the copula approach is introduced and applied to model the joint probability distribution of PGA and CAV. Finally, the correlation results using the PGA‐CAV empirical data and the normalized residuals are compared. The results indicate that there exists a strong positive correlation between PGA and CAV. For both the PGA and CAV empirical data and the normalized residuals, the multivariate lognormal distribution composed of two lognormal marginal distributions and the Gaussian copula provides adequate characterization of the PGA‐CAV joint distribution observed in Taiwan. This finding demonstrates the validity of the conventional two‐step approach for developing empirical ground motion prediction equations (GMPEs) of multiple ground motion parameters from the copula viewpoint. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Investigation on the properties of the relationship between rare and extreme rainfall and flood volumes, under some distributional restrictions

The fact that rainfall data are usually more abundant and more readily regionalized than streamflow data has motivated hydrologists to conceive methods that incorporate the hydrometeorologial information into flood frequency analyses. Some of them, particularly those derived from the French GRADEX method, involve assumptions concerning the relationship between extreme rainfall and flood volumes, under some distributional restrictions. In particular, for rainfall probability distributions exhibiting exponential-like upper tails, it is possible to derive the shape and scale of the probability distribution of flood volumes by hypothesizing the basic properties of such a relationship, under rare and/or extreme conditions. This paper focuses on a parametric mathematical model for the relationship between rare and extreme rainfall and flood volumes under exponentially-tailed distributions. The model is analyzed and fitted to rare and extreme events derived from hydrological simulation of long stochastically-generated synthetic series of rainfall and evaporation for the Indaiá River basin, located in south-central Brazil. The paper also provides a sensitivity analysis of the model parameters in order to better understand flood events under rare and extreme conditions. By working with hydrologically plausible hypothetical events, the modeling approach proved to be a useful way to explore extraordinary rainfall and flood events. The results from this exploratory analysis provide grounds to derive some conclusions regarding the relative positions of the upper tails of the probability distributions of rainfall and flood volumes.     

Multivariate drought characteristics using trivariate Gaussian and Student t copulas

This study aims to investigate the changing properties of drought events in Weihe River basin, China, by modeling the multivariate joint distribution of drought duration, severity and peak using trivariate Gaussian and Student t copulas. Monthly precipitations of Xi'an gauge are used to illustrate the meta‐elliptical copula‐based methodology for a single‐station application. Gaussian and Student t copulas are found to produce a better fit comparing with other six symmetrical and asymmetrical Archimedean copulas, and, checked by the goodness‐of‐fit tests based on a modified bootstrap version of Rosenblatt's transformation, both of them are acceptable to model the multivariate joint distribution of drought variables. Gaussian copula, the best fitting, is employed to construct the dependence structures of positively associated drought variables so as to obtain the multivariate joint and conditional probabilities of droughts. A Kendall's return period (KRP) introduced by Salvadori and De Michele (2010) is then adopted to assess the multivariate recurrent properties of drought events, and its spatial distributions indicate that prolonged droughts are likely to break out with rather short recurrence intervals in the whole region, while drought status in the southeast seems to be severer than the northwest. The study is of some merits in terms of multivariate drought modeling using a preferable copula‐based method, the results of which could serve as a reference for regional drought defense and water resources management. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

13e SESSION DU COMITÉ CONSULTATIF KARACHI 4–15 NOVEMBRE 1957

Abstract

Event-based methods are used in flood estimation to obtain the entire flood hydrograph. Previously, such methods adopted in the UK have relied on pre-determined values of the input variables (e.g. rainfall and antecedent conditions) to a rainfall–runoff model, which is expected to result in an output flood of a particular return period. In contrast, this paper presents a method that allows all the input variables to take on values across the full range of their individual distributions. These values are then brought together in all possible combinations as input to an event-based rainfall–runoff model in a Monte Carlo simulation approach. Further, this simulation strategy produces a long string of events (on average 10 per year), where dependencies from one event to the next, as well as between different variables within a single event, are accounted for. Frequency analysis is then applied to the annual maximum peak flows and flow volumes.

Citation Svensson, C., Kjeldsen, T.R., and Jones, D.A., 2013. Flood frequency estimation using a joint probability approach within a Monte Carlo framework. Hydrological Sciences Journal, 58 (1), 1–20.     

Multivariate homogeneity testing in a northern case study in the province of Quebec,Canada

In regional frequency analysis, the examination of the regional homogeneity represents an important step of the procedure. Flood events possess multivariate characteristics which can not be handled by classical univariate regional procedures. For instance, classical procedures do not allow to assess regional homogeneity while taking into consideration flood peak, volume and duration. Chebana and Ouarda proposed multivariate discordancy and homogeneity tests. They carried out a simulation study to evaluate the performance of these tests. In the present paper, practical aspects are investigated jointly on flood peak and flood volume of a data set from the Côte‐Nord region in the province of Quebec, Canada. It is shown that, after removing the discordant sites, the remaining ones constitute a homogeneous region for the volumes and heterogeneous region for the peaks. However, if both variables are jointly considered, the obtained region is possibly homogeneous. Furthermore, the results demonstrate the usefulness of the bivariate test to take into account the dependence structure between the variables representing the event, and to take advantage of more information from the hydrograph. Copyright © 2009 John Wiley & Sons, Ltd.     

Joint probability distribution of annual maximum storm peaks and amounts as represented by daily rainfalls

Abstract

A procedure is presented for using the bivariate normal distribution to describe the joint distribution of storm peaks (maximum rainfall intensities) and amounts which are mutually correlated. The Box-Cox transformation method is used to normalize original marginal distributions of storm peaks and amounts regardless of the original forms of these distributions. The transformation parameter is estimated using the maximum likelihood method. The joint cumulative distribution function, the conditional cumulative distribution function, and the associated return periods can be readily obtained based on the bivariate normal distribution. The method is tested and validated using two rainfall data sets from two meteorological stations that are located in different climatic regions of Japan. The theoretical distributions show a good fit to observed ones.     

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.     

Response of the flood peak to the spatial distribution of rainfall in the Yom River basin,Thailand

The response of the flood peak to the spatial distribution of rainfall has been reported in basins with nonuniform characteristics. However, prioritization of the influences of these characteristics is still poorly understood. This study evaluated the variability in the flood peak with the spatial distribution of rainfall at Sukhothai (city) in the Yom River basin, Thailand, and investigated the influence of the basin characteristics on the flood peak. For each of the 2-, 5- and 10-y rainfalls with durations of 24, 48 and 72 h, 1000 simulated rainfall events with various spatial distributions were generated according to the observed data by using a Monte Carlo analysis and Cholesky randomization. The floods from these rainfalls were then simulated, and the peak discharges were evaluated. The flood peaks from 24-h rainfalls were usually small but highly variable and could be extremely large when the rainfalls were concentrated over the mountainous region. The flood peaks from 48 to 72-h rainfalls were consistently large and correlated with the rainfalls over the joint area between the mountainous region and plain area. The basin characteristics that influenced the response of the flood peak to the spatial distribution of the rainfall appeared to depend on the rainfall duration and magnitude. For short-duration rainfalls, the response was mainly influenced by the surface storage when the rainfall was small and by the terrain steepness when the rainfall was large. For long-duration rainfalls, the response was mainly influenced by the soil percolation rate.     

Impulse response of a linear diffusion analogy model as a flood frequency probability density function

Abstract

The impulse response of a linear convective-diffusion analogy (LD) model used for flow routing in open channels is proposed as a probability distribution for flood frequency analysis. The flood frequency model has two parameters, which are derived using the methods of moments and maximum likelihood. Also derived are errors in quantiles for these parameter estimation methods. The distribution shows that the two methods are equivalent in terms of producing mean values—the important property in case of unknown true distribution function. The flood frequency model is tested using annual peak discharges for the gauging sections of 39 Polish rivers where the average value of the ratio of the coefficient of skewness to the coefficient of variation equals about 2.52, a value closer to the ratio of the LD model than to the gamma or the lognormal model. The likelihood ratio indicates the preference of the LD over the lognormal for 27 out of 39 cases. It is found that the proposed flood frequency model represents flood frequency characteristics well (measured by the moment ratio) when the LD flood routing model is likely to be the best of all linear flow routing models.     

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.     

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.     

Estimation of flood inundation probabilities using global hazard indexes based on hydrodynamic variables

In this paper a new procedure to derive flood hazard maps incorporating uncertainty concepts is presented. The layout of the procedure can be resumed as follows: (1) stochastic input of flood hydrograph modelled through a direct Monte-Carlo simulation based on flood recorded data. Generation of flood peaks and flow volumes has been obtained via copulas, which describe and model the correlation between these two variables independently of the marginal laws involved. The shape of hydrograph has been generated on the basis of a historical significant flood events, via cluster analysis; (2) modelling of flood propagation using a hyperbolic finite element model based on the DSV equations; (3) definition of global hazard indexes based on hydro-dynamic variables (i.e., water depth and flow velocities). The GLUE methodology has been applied in order to account for parameter uncertainty. The procedure has been tested on a flood prone area located in the southern part of Sicily, Italy. Three hazard maps have been obtained and then compared.     

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