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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Ensemble prediction of regional droughts using climate inputs and the SVM–copula approach

In this study, the climate teleconnections with meteorological droughts are analysed and used to develop ensemble drought prediction models using a support vector machine (SVM)–copula approach over Western Rajasthan (India). The meteorological droughts are identified using the Standardized Precipitation Index (SPI). In the analysis of large‐scale climate forcing represented by climate indices such as El Niño Southern Oscillation, Indian Ocean Dipole Mode and Atlantic Multidecadal Oscillation on regional droughts, it is found that regional droughts exhibits interannual as well as interdecadal variability. On the basis of potential teleconnections between regional droughts and climate indices, SPI‐based drought forecasting models are developed with up to 3 months' lead time. As traditional statistical forecast models are unable to capture nonlinearity and nonstationarity associated with drought forecasts, a machine learning technique, namely, support vector regression (SVR), is adopted to forecast the drought index, and the copula method is used to model the joint distribution of observed and predicted drought index. The copula‐based conditional distribution of an observed drought index conditioned on predicted drought index is utilized to simulate ensembles of drought forecasts. Two variants of drought forecast models are developed, namely a single model for all the periods in a year and separate models for each of the four seasons in a year. The performance of developed models is validated for predicting drought time series for 10 years' data. Improvement in ensemble prediction of drought indices is observed for combined seasonal model over the single model without seasonal partitions. The results show that the proposed SVM–copula approach improves the drought prediction capability and provides estimation of uncertainty associated with drought predictions. Copyright © 2013 John Wiley & Sons, Ltd.     

Socio-economic drought assessment in Lake Mead,USA, based on a multivariate standardized water-scarcity index

Drought events are increasing worldwide. Socio-economic drought is the least investigated type of drought, and is the only type for which water demand is taken into consideration. In this research, socio-economic drought was studied in Lake Mead, USA, using a multivariate standardized water-scarcity index (MSWSI) over the period 1990–2014, combining two water-scarcity indices based on time series of inflow, outflow, reservoir storage, and water demand. The inflow and outflow were determined from streamgage data, and reservoir storage from lake level data; demand was based on water pumped by the Southern Nevada Water Authority. Missing observations in input streamgage data were filled through regression modeling. The results indicate that Lake Mead has been in socio-economic drought since 2000, with the most severe drought occurring between 2006 and 2012, and the highest intensity drought in April–July 2014. The Lake Mead droughts revealed through the MSWSI are consistent with those reported in US Drought Monitor (USDM) products. The temporal behavior of MSWSI provides an insight into the socio-economic effects of droughts not captured by USDM products.     

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.     

Multivariate Standardized Drought Index: A parametric multi-index model

Defining droughts based on a single variable/index (e.g., precipitation, soil moisture, or runoff) may not be sufficient for reliable risk assessment and decision-making. In this paper, a multivariate, multi-index drought-modeling approach is proposed using the concept of copulas. The proposed model, named Multivariate Standardized Drought Index (MSDI), probabilistically combines the Standardized Precipitation Index (SPI) and the Standardized Soil Moisture Index (SSI) for drought characterization. In other words, MSDI incorporates the meteorological and agricultural drought conditions for overall characterization of drought. In this study, the proposed MSDI is utilized to characterize the drought conditions over several Climate Divisions in California and North Carolina. The MSDI-based drought analyses are then compared with SPI and SSI. The results reveal that MSDI indicates the drought onset and termination based on the combination of SPI and SSI, with onset being dominated by SPI and drought persistence being more similar to SSI behavior. Overall, the proposed MSDI is shown to be a reasonable model for combining multiple indices probabilistically.     

Space–time modelling of catchment scale drought characteristics

Drought may affect all components of the water cycle and covers commonly a large part of the catchment area. This paper examines drought propagation at the catchment scale using spatially aggregated drought characteristics and illustrates the importance of catchment processes in modifying the drought signal in both time and space. Analysis is conducted using monthly time series covering the period 1961–1997 for the Pang catchment, UK. The time series include observed rainfall and groundwater recharge, head and discharge simulated by physically-based soil water and groundwater models. Drought events derived separately for each unit area and variable are combined to yield catchment scale drought characteristics. The study reveals relatively large differences in the spatial and temporal characteristics of drought for the different variables. Meteorological droughts cover frequently the whole catchment; and they are more numerous and last for a short time (1–2 months). In comparison, droughts in recharge and hydraulic head cover typically a smaller area and last longer (4–5 months). Hydraulic head and groundwater discharge exhibit similar drought characteristics, which can be expected in a groundwater fed catchment. Deficit volume is considered a robust measure of the severity of a drought event over the catchment area for all variables; whereas, duration is less sensitive, particular for rainfall. Spatial variability in drought characteristics for groundwater recharge, head and discharge are primarily controlled by catchment properties. It is recommended not to use drought area separately as a measure of drought severity at the catchment scale, rather it should be used in combination with other drought characteristics like duration and deficit volume.     

Analysis of streamflow droughts using fixed and variable thresholds

Hydrological drought analysis is very important in the design of hydrotechnical projects and water resources management and planning. In this study, a methodology is proposed for the analysis of streamflow droughts using the threshold level approach. The method has been applied to Yermasoyia semiarid basin in Cyprus based on 30‐year daily discharge data. Severity was defined as the accumulated water deficit volume occurring during a drought event, in respect with a target threshold. Fixed and variable thresholds (seasonal, monthly, and daily) were employed to derive the drought characteristics. The threshold levels were determined based on the Q₅₀ percentiles of flow extracted from the corresponding flow duration curves for each threshold. The aim is to investigate the sensitivity of these thresholds in the estimation of maximum drought severities for various return periods and the derivation of severity–duration–frequency curves. The block maxima and the peaks over threshold approaches were used to perform the extreme value analysis. Three pooling procedures (moving average, interevent time criterion, and interevent time and volume criterion) were employed to remove the dependent and minor droughts. The application showed that the interevent time and volume criterion is the most unbiased pooling method. Therefore, it was selected to estimate the drought characteristics. The results of this study indicate that monthly and daily variable thresholds are able to capture abnormal drought events that occur during the whole hydrological year whereas the other two, only the severe ones. They are also more sensitive in the estimation of maximum drought severities and the derivation of the curves because they incorporate better the effect of drought durations.     

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.     

Probabilistic spatial prediction of categorical data using elliptical copulas

This study uses elliptical copulas and transition probabilities for uncertainty modeling of categorical spatial data. It begins by discussing the expressions of the cumulative distribution function and probability density function of two major elliptical copulas: Gaussian copula and t copula. The basic form of spatial copula discriminant function is then derived based on Bayes’ theorem, which consists of three parts: the prior probability, the conditional marginal densities, and the conditional copula density. Finally, three kinds of parameter estimation methods are discussed, including maximum likelihood estimation, inference functions for margins and canonical maximum likelihood (CML). To avoid making assumptions on the form of marginal distributions, the CML approach is adopted in the real-world case study. Results show that the occurrence probability maps generated by these two elliptical copulas are similar to each other. However, the prediction map interpolated by Gaussian copula has a relatively higher classification accuracy than t copula.     

Evaluation, modification, and application of the Effective Drought Index to 200-Year drought climatology of Seoul, Korea

In this paper, progressive methods for assessing drought severity from diverse points of view were conceived. To select a fundamental drought index, the performances of the Effective Drought Index (EDI) and 1-, 3-, 6-, 9-, 12-, and 24-month Standardized Precipitation Indices (SPIs) were compared for drought monitoring data accumulated over 200-year period from 1807 to 2006 for Seoul, Korea. The results confirmed that the EDI was more efficient than the SPIs in assessing both short and long-term droughts.We then proposed the following methods for modifying and supplementing the EDI: (1) CEDI, a corrected EDI that considers the rapid runoff of water resources after heavy rainfall; (2) AEDI, an accumulated EDI that considers the drought severity and duration of individual drought events; and (3) YAEDI, a year-accumulated negative EDI representing annual drought severity. In addition to these indices, to more accurately measure and diagnose droughts, we proposed the utilization of (4) the Available Water Resources Index (AWRI), an existing index that expresses the actual amount of available water.Using the improved methods above, we assessed and summarized important droughts that have occurred in Seoul over the 200 years from 1807 to 2006.