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.     

Regional bivariate modeling of droughts using L-comoments and copulas

The regional bivariate modeling of drought characteristics using the copulas provides valuable information for water resources management and drought risk assessment. The regional frequency analysis (RFA) can specify the similar sites within a region using L-comoments approach. One of the important steps in the RFA is estimating regional parameters of the copula function. In the present study, an optimization-based method along with the adjusted charged system search are introduced and applied to estimate the regional parameters of the copula models. The capability of the proposed methodology is illustrated by copula functions on drought events. Three commonly used copulas containing Clayton, Frank and Gumbel are employed to derive the joint distribution of drought severity and duration. The result of the new method are compared to the method of moments and after applying several goodness-of-fit tests, the results indicate that the new method provides higher accuracy than the classic one. Furthermore, the results of the upper tail dependence coefficient indicate that the Gumbel copula is the best-fitted copula among the other ones for modeling drought characteristics.     

Assessment of hydrological droughts for the Yellow River,China, using copulas

Droughts are one of the normal and recurrent climatic phenomena on Earth. However, recurring prolonged droughts have caused far‐reaching and diverse impacts because of water deficits. This study aims to investigate the hydrological droughts of the Yellow River in northern China. Since drought duration and drought severity exhibit significant correlation, a bivariate distribution is used to model the drought duration and severity jointly. However, drought duration and drought severity are often modelled by different distributions; the commonly used bivariate distributions cannot be applied. In this study, a copula is employed to construct the bivariate drought distribution. The copula is a function that links the univariate marginal distributions to form the bivariate distribution. The bivariate return periods are also established to explore the drought characteristics of the historically noticeable droughts. The results show that the return period of the drought that occurred in late 1920s to early 1930s is 105 years. The significant 1997 dry‐up phenomenon that occurred in the downstream Yellow River (resulting from the 1997–1998 drought) only has a return period of 4·4 years and is probably induced by two successive droughts and deteriorated by other factors, such as human activities. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of copulas for derivation of drought severity–duration–frequency curves

This study presents copula‐based multivariate probabilistic approach to model severity–duration–frequency (S‐D‐F) relationship of drought events in western Rajasthan, India. Drought occurrences are analysed using standardized precipitation index computed on monthly mean areal precipitation, aggregated at a time scale of 6 months. After testing with a series of probability density functions, the drought variable severity is found to be better represented with log‐normal distribution, whereas duration is well fitted with exponential distribution. Four different classes of bivariate copulas – Archimedean, extreme value, Plackett, and elliptical families are evaluated for modelling joint distribution of drought characteristics. It is observed that the extreme value copula – Gumbel–Hougaard copula – performed better as compared with other classes of copulas, based on results of various statistical tests and upper tail dependence coefficient. The joint distribution obtained from best performing copula is then employed to determine conditional return period and to derive drought severity‐duration‐frequency (S‐D‐F) curves for the study region. The results of the study suggests that the copula method can be used effectively to derive the drought S‐D‐F curves, which can be helpful in planning and adopting suitable drought mitigation strategies in drought‐prone areas. Copyright © 2011 John Wiley & Sons, Ltd.     

Frequency analysis of droughts using the Plackett copula and parameter estimation by genetic algorithm

This study aims to model the joint probability distribution of drought duration, severity and inter-arrival time using a trivariate Plackett copula. The drought duration and inter-arrival time each follow the Weibull distribution and the drought severity follows the gamma distribution. Parameters of these univariate distributions are estimated using the method of moments (MOM), maximum likelihood method (MLM), probability weighted moments (PWM), and a genetic algorithm (GA); whereas parameters of the bivariate and trivariate Plackett copulas are estimated using the log-pseudolikelihood function method (LPLF) and GA. Streamflow data from three gaging stations, Zhuangtou, Taian and Tianyang, located in the Wei River basin, China, are employed to test the trivariate Plackett copula. The results show that the Plackett copula is capable of yielding bivariate and trivariate probability distributions of correlated drought variables.     

Bivariate drought frequency curves and confidence intervals: a case study using monthly rainfall generation

Although water resources management practices recently use bivariate distribution functions to assess drought severity and its frequency, the lack of systematic measurements is the major hindrance in achieving quantitative results. This study aims to suggest a statistical scheme for the bivariate drought frequency analysis to provide comprehensive and consistent drought severities using observed rainfalls and their uncertainty using synthesized rainfalls. First, this study developed a multi-variate regression model to generate synthetic monthly rainfalls using climate variables as causative variables. The causative variables were generated to preserve their correlations using copula functions. This study then focused on constructing bivariate drought frequency curves using bivariate kernel functions and estimating their confidence intervals from 1,000 likely replica sets of drought frequency curves. The confidence intervals achieved in this study may be useful for making a comprehensive drought management plan through providing feasible ranges of drought severity.     

Uncertainty and variability in bivariate modeling of hydrological droughts

There are two kinds of uncertainty factors in modeling the bivariate distribution of hydrological droughts: the alteration of predefined critical ratios for pooling droughts and excluding minor droughts and the temporal variability of univariate and/or bivariate characteristics of droughts due to the impact of human activities. Daily flow data covering a period of 56 hydrological years from two gauging stations from a humid region in South China are used. The influences of alterations of threshold values of flow and critical ratios of pooling droughts and excluding minor droughts on drought properties are analyzed. Six conventional univariate models and three Archimedean copulas are employed to fit the marginal and joint distributions of drought properties, the Kolmogorov–Smirnov and Anderson–Darling methods are used for testing the goodness-of-fit of the univariate model, and the Cramer-von Mises method based on Rosenblatt’s transform is applied for the test of the bivariate model. The change point analysis of the copula parameter of bivariate distribution of droughts is first made. Results demonstrate that both the statistical characteristics of each drought property and their bivariate joint distributions are sensitive to the critical ratio of excluding minor droughts. A model can be selected to fit the marginal distribution for drought deficit volume or maximum deficit, but it is not determined for drought duration with the lower ratios of the pooling and excluding droughts. The statistical uncertainty of drought duration makes the modeling of bivariate joint distribution of drought duration and deficit volume or of drought duration and maximum deficit undermined. Change points significantly occurred in the period from the late 1970s to the middle 1980s for a single drought property and the copula parameter of their joint distribution due to the impact of human activities. The difference between two subseries separated by the change point is remarkable in the magnitudes of drought properties and the joint return periods. A copula function can be selected to optimally fit the bivariate distribution, provided that the critical ratios of pooling and excluding droughts are great enough such as the optimal value of 0.4 in the case study. It is valuable that the modeling and designing of the bivariate joint correlation and distribution of drought properties can be performed on the subseries separated by the change point of the copula parameter.     

Drought hotspot analysis and risk assessment using probabilistic drought monitoring and severity–duration–frequency analysis

Drought hotspot identification requires continuous drought monitoring and spatial risk assessment. The present study analysed drought events in the agriculture‐dominated mid‐Mahanadi River Basin in Odisha, India, using crop water stress as a drought indicator. This drought index incorporated different factors that affect crop water deficit such as the cropping pattern, soil characteristics, and surface soil moisture. The drought monitoring framework utilized a relevance vector machine model‐based classification that provided the uncertainty associated with drought categorization. Using the proposed framework, drought hotspots are identified in the study region and compared with indices based on precipitation and soil moisture. Further, a bivariate copula is employed to model the agricultural drought characteristics and develop the drought severity–duration–frequency (S–D–F) relationships. The drought hotspot maps and S–D–F curves are developed for different locations in the region. These provided useful information on the site‐specific drought patterns and the characteristics of the devastating droughts of 2002 and 2012, characterized by an average drought duration of 7 months at several locations. The site‐specific risk of short‐ and long‐term agricultural droughts are then investigated using the conditional copula. The results suggest that the conditional return periods and the S–D–F curves are valuable tools to assess the spatial variability of drought risk in the region.     

Multivariate modeling of droughts using copulas and meta-heuristic methods

This study investigated the utility of two meta-heuristic algorithms to estimate parameters of copula models and for derivation of drought severity–duration–frequency (S–D–F) curves. Drought is a natural event, which has huge impact on both the society and the natural environment. Drought events are mainly characterized by their severity, duration and intensity. The study adopts standardized precipitation index for drought characterization, and copula method for multivariate risk analysis of droughts. For accurate estimation of copula model parameters, two meta-heuristic methods namely genetic algorithm and particle swarm optimization are applied. The proposed methodology is applied to a case study in Trans Pecos, an arid region in Texas, USA. First, drought severity and duration are separately modeled by various probability distribution functions and then the best fitted models are selected for copula modeling. For modeling the joint dependence of drought variables, different classes of copulas, namely, extreme value copulas, Plackett and Student’s t copulas are employed and their performance is evaluated using standard performance measures. It is found that for the study region, the Gumbel–Hougaard copula is the best fitted copula model as compared to the others and is used for the development of drought S–D–F curves. Results of the study suggest that the meta-heuristic methods have greater utility in copula-based multivariate risk assessment of droughts.     

Spatio-temporal analysis and derivation of copula-based intensity–area–frequency curves for droughts in western Rajasthan (India)

This study presents spatio-temporal analysis of droughts in one of the most drought prone region in India–western Rajasthan and develops drought intensity-area-frequency curves for the region. The meteorological drought conditions are analyzed using 6-month standardized precipitation index (SPI-6) estimated at spatial resolution of 0.5° × 0.5°. Spatio-temporal analysis of SPI-6 indicates increase in frequency of droughts at the central part of the region. The non-parametric Mann–Kendall test for seasonal trend analysis showed increase in number of grids under drought during the study period. Further, bivariate frequency analysis of drought characteristics—intensity and areal extent is carried out using copula methods. For modeling joint dependence between drought variables, three copula families namely Gumbel-Hougaard, Frank and Plackett copulas are evaluated. Based on goodness-of-fit as well as upper tail dependence tests, it is found that the Gumbel-Hougaard copula best represents the drought properties. The copula-based joint distribution is used to compute conditional return periods and drought intensity–area–frequency (I–A–F) curves. The I–A–F curves could be helpful in risk evaluation of droughts in the region.     

A bivariate pareto model for drought

Univariate Pareto distributions have been so widely used in hydrology. It seems however that bivariate or multivariate Pareto distributions have not yet found applications in hydrology, especially with respect to drought. In this note, a drought application is described by assuming a bivariate Pareto model for the joint distribution drought durations and drought severity in the State of Nebraska. Based on this model, exact distributions are derived for the inter arrival time, magnitude and the proportion of droughts. Estimates of 2, 5, 10, 20, 50 and 100 year return periods are derived for the three variables, drought duration, drought severity and the pairwise combinations: (drought duration, drought severity), (inter arrival time of drought, proportion of drought) and (drought duration, drought magnitude). These return period estimates could have an important role in hydrology, for example, with respect to measures of vegetation water stress for plants in water-controlled ecosystems.     

西江水文干旱历时与强度的遭遇概率分析

利用西江下游马口水文站1959 2009年月径流量数据计算径流干旱指数,经游程理论提取了水文干旱特征值.应用Copula函数分析水文干旱强度和历时之间的联合概率分布.对构建的干旱历时和强度联合分布模式进行分析,结果表明:(1)径流干旱历时和强度之间具有高关联性,秩相关系数达0.617;(2)三参数Weibull分布较好地描述了干旱历时和强度的边缘分布特征;(3)经拟合优度检验结果优选的干旱历时和强度之间的较优连接函数为Archimedean类的Gumbel-Hougaard Copula函数;(4)5～10年重现期和20年重现期的水文干旱分别达到了重旱级别和特旱级别;(5)干旱历时和强度之间的遭遇概率可为特定干旱历时与水文干旱级别或特定干旱强度与干旱历时之间的对应关系提供概率意义上的干旱特征诊断与预测.     

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.     

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     

On the generation of drought events using an alternating renewal-reward model

This research study focused on the hypothesis that extreme drought and high streamflow events come from different independent populations with different probability distributions which need to be studied separately, rather than considering the streamflow population as a whole. The inability of traditional streamflow generator models to consistently reproduce the frequency of occurrence of severe droughts observed in the historical record has been questioned by many researchers. Our study focused on the development of astochastic event generator model which would be capable of doing so. This was accomplished in a two-step process by first generating the drought event, and then deriving the streamflows which comprised that event. The model considered for this analysis was an alternating renewal-reward procedure that cycles between eventon andoff times, and is representative of drought or high streamflow event duration. The reward gained while the event ison oroff represents drought severity or high streamflow surplus. Geometric and gamma distributions were considered for drought duration and deficit respectively. Model validation was performed using calculated required capacities from the sequent peak algorithm.     

Assessment of regional drought risk under climate change using bivariate frequency analysis

Since water supply failure is one of the primary impacts of drought, drought risk should be quantified in the context of a lack of available water. To assess the drought risk, water supply system performance indices such as reliability, resiliency, and vulnerability are usually introduced as they correspond to primary drought characteristics, i.e., frequency, duration, and magnitude. In this study, we developed a drought risk index (DRI) through weighted averaging the performance indices derived using bivariate drought frequency analysis. We suggested two types of DRI: observed DRI (DRI_O) and designed DRI (DRI_D). DRI_O was calculated using an observed (or synthesized) time series of water shortages. DRI_D was estimated from the bivariate drought frequency curves, which are the probabilistic magnitudes of water shortages corresponding to a particular duration. The historical maximum drought event that represents the maximum DRI_O has generally been used as the target security level. However, we could establish a practically applicable target security level considering that the future water supply failure risk is represented by DRI_D. We defined regional drought safety criteria in this study by comparing DRI_O and DRI_D. Application of the criteria to the Nakdong river basin in South Korea showed that W1 (Byeongseongcheon) and W4 (Hyeongsangang) had the lowest and highest drought risk, respectively, and the drought safety criteria showed an average range of 5–20 years.     

A drought frequency formula

Abstract

The important elements of a drought phenomenon are the longest duration and the largest severity for a desired return period. These elements form a basis for designing water storage systems to cope with droughts. At times, a third element, drought intensity, is also used and is defined as the ratio of severity to duration. The commonly available statistics for the causative drought variables such as annual rainfall or runoff sequences are the mean, the coefficient of variation and the lag one serial correlation coefficient, and occasionally some indication of the probability distribution function (pdf) of the sequences. The extremal values of the duration and severity are modelled in the present paper using information on the aforesaid parameters at the truncation level equal to the mean of the drought sequence, which is generally taken as the truncation level in the analysis of droughts. The drought severity has been modelled as the product of the duration and intensity with the assumption of independence between them. An estimate of drought intensity has been realized from the concept of the truncated normal distribution of the standardized form of the drought sequences in the normalized domain. A formula in terms of the extremal severity and the T-year return period has been suggested similar to the flood frequency formulae, commonly cited in hydrological texts.     

A multimodel assessment of the climate change effect on the drought severity–duration–frequency relationship

In this study, the patterns of past and future drought occurrences in the Seoul region were analysed using observed historical data from the Seoul weather station located in the Korean Peninsula and four different types of general circulation models (GCMs), namely, GFDL:CM2_1, CONS:ECHO‐G, MRI:CGCM2_3_2 and UKMO:HADGEM1. To analyse statistical properties such as drought frequency duration and return period, the Standardized Precipitation Index was used to derive the severity–duration–frequency (SDF) curve from the drought frequency analysis. In addition, a drought spell analysis was conducted to estimate the frequency and change of drought duration for each drought classification. The results of the analysis suggested a decrease in the frequency of mild droughts and an increase in the frequency of severe and extreme droughts in the future. Furthermore, the average duration of droughts is expected to increase. A comparison of the SDF relationship derived from the observed data with that derived via the GCMs indicated that the drought severity for each return period was reduced as drought duration increased and that the drought severity derived from the GCMs was severer than the severity obtained using the observed data for the same duration and return period. Furthermore, among the four types of GCMs used in this study, the MRI model predicted the most severe future drought for the Seoul region, and the SDF curve derived using the MRI model also resulted in the highest degree of drought severity compared with the other GCMs. Copyright © 2012 John Wiley & Sons, Ltd.     

A new bivariate Gamma distribution generated from functional scale parameter with application to drought data

Univariate and bivariate Gamma distributions are among the most widely used distributions in hydrological statistical modeling and applications. This article presents the construction of a new bivariate Gamma distribution which is generated from the functional scale parameter. The utilization of the proposed bivariate Gamma distribution for drought modeling is described by deriving the exact distribution of the inter-arrival time and the proportion of drought along with their moments, assuming that both the lengths of drought duration (X) and non-drought duration (Y) follow this bivariate Gamma distribution. The model parameters of this distribution are estimated by maximum likelihood method and an objective Bayesian analysis using Jeffreys prior and Markov Chain Monte Carlo method. These methods are applied to a real drought dataset from the State of Colorado, USA.