## 1. Introduction

Accurate forecasting of hydrologic extreme events plays a significant role in developing appropriate policies to plan for available water resources. Although several studies have proposed promising methods to improve hydrologic forecasts, the observed effects of climate change on floods and droughts across different regions of the globe in recent decades highlights the need for more sophisticated methods in predicting extreme events (Mishra and Singh 2010; Moradkhani et al. 2010; Halmstad et al. 2012; Risley et al. 2011; Madadgar and Moradkhani 2011; Najafi et al. 2012). Motivated by the increasing frequency of droughts throughout the globe, this study proposes an advanced methodology for forecasting droughts within a probabilistic framework.

Drought development is rooted in shortages of precipitation over an extended period of time, resulting in water scarcity (Jaeger et al. 2013). As droughts expand across different spatial and temporal scales, the natural environment and human lives are strongly affected. Almost all the continents throughout the globe have been affected by various drought phenomena during recent decades (Mishra and Singh 2010). According to the U.S. Drought Monitor, more than 70% of the United States currently experiences dry spells of severities ranging from abnormally dry to exceptional droughts (Showstack 2012). Several recent efforts have attempted to enhance forecast accuracy, mitigation policies, and damage estimate of drought events in the globe, specifically in the United States. The Drought Impact Reporter (DIR), launched by National Drought Mitigation Center (NDMC), is a comprehensive database reporting damages caused by recent droughts within the United States. Reported by the North America Drought Monitor of National Oceanic and Atmospheric Administration (NOAA) National Climate Data Center (NCDC; http://www.ncdc.noaa.gov/temp-and-precip/drought/nadm/), droughts with an estimated damage of over 100 billion dollars (Lott and Ross 2000) have been among the costliest natural disasters in the United States since 1980. Lott and Ross (2006) estimated drought and heat wave induced damages to the U.S. economy at $174 billion, between 1980 and 2005. However, drought is a slowly developed phenomenon, and hence there might be a chance to mitigate drought impacts if the events are forecasted within an appropriate time.

There are various studies in the literature examining different aspects of drought phenomena. Drought identification and monitoring (Andreadis and Lettenmaier 2006; Sheffield and Wood 2007; Shukla et al. 2011; Van Loon and Van Lanen 2012; van Huijgevoort et al. 2012), developing drought indices (Niemeyer 2008; Mishra and Singh 2010), characterization of drought events (Shiau 2006; Dupuis 2007; Wong et al. 2010), evaluation of global warming and climate change impacts on future droughts (Blenkinsop and Fowler 2007; Ghosh and Mujumdar 2007; Sheffield and Wood 2008; Burke et al. 2010; Risley et al. 2011; Madadgar and Moradkhani 2011), early warning systems in reservoir operation under drought stress (Huang and Chou 2005, 2008), and many other drought-related analyses have been studied during last years. There are also several studies that have focused on qualified drought forecast for different regions. In an early attempt, Karl et al. (1987) analyzed the probability of receiving the amount of precipitation required to terminate an ongoing dry condition and return to the normal conditions over various time periods. The limitation, however, was using an unconditional gamma distribution—ignoring the dependency and correlation of precipitation in temporal scale—to obtain the probabilities. Lohani and Loganathan (1997) used a nonhomogeneous Markov chain model to generate the transition probability matrix of drought states. In another study, the Markov chain model was employed to evaluate drought transition probabilities, persistence, duration, and frequency within six categories of different severities (Steinemann 2003). Some other studies used the stochastic renewal models, stochastic autoregressive models, and artificial neural networks (ANN) to estimate different characteristics of future droughts and low-flow periods (Kendall and Dracup 1992; Loaiciga and Leipnik 1996; Mishra and Desai 2005, 2006; Barros and Bowden 2008). However, the autoregressive and neural network models were later questioned by Hwang and Carbone (2009) because of limiting the forecasts into the deterministic estimate of the mean drought status. Recently, Özger et al. (2012) developed a wavelet and fuzzy logic combination model for long-lead drought forecasting. The technique was found to outperform fuzzy logic, ANN, or coupled wavelet and fuzzy logic models, yet prior to an application it needs a significant work to find the appropriate independent predictors, which strongly affect the forecast. Without using any frequency-analysis methods, Cancelliere et al. (2007) derived the transition probabilities matrix by revising the statistics underlying the standardized precipitation index (SPI) series. They also questioned the validity of a Markov chain model in forming the transition probability matrix for forecasting SPI values. However, two major limitations of their study are 1) there are promising approaches like copula functions, as discussed later, for frequency analyses of drought status that look promising and avoid overwhelming procedures to analytically derive the transition probability matrix from the index formulas; 2) to reduce the computational burden, they assumed that aggregated monthly precipitations are uncorrelated and normally distributed variables, while this is not a valid assumption specifically when the method is expanded to other hydrologic variables like streamflow. Some other studies took advantage of seasonal climate forecasts to predict the future droughts. Carbone and Dow (2005) and Hwang and Carbone (2009) incorporated the seasonal forecast products of NOAA Climate Prediction Center (CPC) with historical climate records to address the uncertainties of future droughts. Despite the potential of CPC products in drought forecasting, Steinemann (2006) argued that the water managers are not usually able to properly interpret the forecast probabilities and uncertainty information provided by CPC seasonal precipitation outlooks.

This study approaches probabilistic seasonal drought forecast using copula functions (Joe 1997; Nelsen 1999). Despite the fact that copulas are admired in probability theory and statistics, they have yet to be effectively applied in probabilistic forecasts of drought events, except in a few studies in determining drought characteristics. Copulas are distribution functions that join several variables of diverse correlation and dependence structures. Since hydrologic variables are usually correlated to each other, the copula functions have been applied in probabilistic modeling of various hydrologic phenomena over the past few years, including flood analyses (e.g., Favre et al. 2004; Zhang and Singh 2007a; Salvadori and De Michele 2010), rainfall analyses (Zhang and Singh 2007b; Salvadori and De Michele 2006; Kao and Govindaraju 2008), spatial analysis of groundwater quality parameters (Bárdossy 2006; Bárdossy and Li 2008), synthesizing and downscaling of monthly river flow (Li et al. 2013), and low-flow/drought analyses (Shiau 2006; Dupuis 2007; Kao and Govindaraju 2010; Wong et al. 2010; Madadgar and Moradkhani 2011). Recently, Li et al. (2013) employed copula functions to synthesize and downscale the monthly river flows conditional on the past observations. In another recent study, a copula was employed in the postprocessing of hydrologic forecasts (Madadgar et al. 2013), and its performance was tested on 2500 hypothetical case studies as well as a real case study of streamflow forecast ensembles for Sprague River basin in southern Oregon. By reducing the uncertainty and increasing the reliability of forecasts, the developed technique could efficiently improve the streamflow predictions, and it was found to be a robust alternative for postprocessing or bias correction of hydrologic forecasts. This study aims at modeling and forecasting the hydrologic droughts regarding seasonal flows. Since the seasonal flow and corresponding drought conditions are dependent variables in time, copulas as capable tools are incorporated to model the Bayesian network of sequential events. Therefore, a future drought status, given the earlier drought conditions, is estimated using copula functions within a Bayesian network.

Aside from the probabilistic approach presented in this study to develop the seasonal drought forecasting, we also use the forecast tools incorporated by the operational streamflow forecast centers to evaluate how various forecast techniques are in agreement in seasonal drought forecasting. Ensemble Streamflow Prediction (ESP; Twedt et al. 1977; Day 1985) produced by the National Weather Service (NWS) River Forecast Centers (RFC) is a way to address the uncertainties of the future flows. Promising forecasts by ESPs have encouraged the researchers to extensively use them in estimating the uncertainties of various hydrometeorological forecasts, and several techniques have been recently developed to build qualified ESPs (e.g., Schaake et al. 2007; Wood and Lettenmaier 2008; Brown and Seo 2010; Dechant and Moradkhani 2011; Madadgar et al. 2013; Najafi et al. 2012). In seasonal flow forecasting, ESPs are generated using the watershed models driven by the resampled historical meteorology. Since the water supply forecasts are developed in advance when the future climate has not yet being recorded, the historical meteorology with an uncertainty of not knowing the exact future conditions is assumed to reasonably treat the future climatology. Therefore, a watershed model forced by the resampled historical meteorology produces an ensemble of streamflow forecasts that reflect the uncertainties of future unseen climatology.

The paper is organized as follows. Section 2 explains the Bayesian networks and adopts copula modeling to define the conditional probabilities of the Bayesian networks. Section 3 presents drought forecasting techniques and section 4 outlines the study basin. The drought index used in this study is described in section 5 followed by the analysis of historical droughts and developing probabilistic drought forecasts. The forecast products are presented for the first-order and second-order conditional forecasts; last, section 8 provides the final remarks and conclusion.

## 2. Copula-based modeling of Bayesian networks

*i*th variable, and

## 3. Drought forecasting methods

This study aims at forecasting seasonal droughts conditional on the drought status of the earlier seasons with highest correlations. As shown later for the study basin, the target season is correlated the most to its previous season; however, its correlation to the next earlier season (two prior seasons) is not insignificant enough to ignore. Thus, Eq. (9) and (10) are applied in the remainder of this paper to practice the probabilistic drought analysis of the target season given the drought status of either one or two earlier seasons.

To provide an operational insight toward the developed forecast methodology, the presented technique is compared with ESP (Twedt et al. 1977; Day 1985), having been used by various operational hydrological forecast centers including the NWS. ESPs are generated to characterize the uncertainty of hydrological forecasts. They use the observed meteorology in a historical time period to reflect on the unseen future climate. To generate the more reliable ESPs, the climate data might be weighted such that it reflects the similarity of the current climate condition to those of the historical data (Najafi et al. 2012). ESP can also be generated by driving the hydrologic model with the climate forecast. To predict the seasonal flow and generate the ESPs, this study uses the Precipitation-Runoff Modeling System (PRMS; Leavesley et al. 1983), which is a distributed-parameter watershed model developed by the U.S. Geological Survey (USGS) to simulate the effects of various combinations of climate, land use, soil type, etc. on the hydrologic response of the watersheds. The watershed model is driven by the historical climatology during the spinup period before the forecast date. Beginning from the forecast date, the model is forced by the resampled historical meteorology to produce an ensemble of hydrologic forecasts. Since the resampled climate data reasonably reflects the uncertainty of the unseen future meteorology, the generated ESP is assumed to properly model the uncertainty of future hydrology caused by unknown climatology. To implement the retrospective forecast of drought events in the Gunnison River basin, the duration of 1980–90 is taken as the spinup period and the ESPs are generated starting from 1 January of each year during 1990–2011 with a lead time of 6 months. Therefore, each ESP is built for the period of January–June of each year, whereas the PRMS is driven with the resampled meteorology of January–June from each year in the entire period (1980–2011)—except the year that the ESP is generated for.

## 4. Case study and data

The Gunnison River basin is one of the headwater subbasins of the Colorado River basin, located in the southwestern United States (Fig. 1). The Colorado River basin, with an approximate drainage area of 640 000 km^{2}, is divided into upper and lower portions, and encompasses parts of seven states: Wyoming, Colorado, Utah, Nevada, California, New Mexico, and Arizona. The Gunnison River basin includes seven subbasins with a total drainage area of 5400 km^{2} at the conjunction of two upstream reaches: the Tomichi Creek and the Gunnison. Streamflow observations of the upstream reaches immediately before the basin outlet at USGS 09119000 (Tomichi Creek River) and USGS 09114500 (Gunnison River) are accumulated to use as the basin’s total outflow at any time. Streamflow record during 1950–90 is contributed to establish and train the forecast model and the rest (1990–2011) is used for verification purposes.

The Gunnison River basin is a snowmelt dominated watershed. According to the drought summary by Western Water Assessment (WWA) and National Integrated Drought Information System (NIDIS) released in July 2012 (WWA and NIDIS 2012), depletion of the snowpack and the early meltout in the spring of 2012 caused the below-average flow in April–July of 2012. During the past 118 yr, 2012 was the second warmest year on record in the state of Colorado. Regarding the inflows to Lake Powell, which reflects the runoff of the entire upper Colorado River basin, the water year 2012 was the fourth driest year in the past century. The intense drought of 2012 (comparable with the drought of 2002 across the region) in the states of Colorado, Utah, and Wyoming caused insufficient water supply, poor pasture and crop conditions, and regionwide wildfires. However, in general, the region has been undergoing various droughts since 2000, with the most intense drought occurring in 2002. Thus accurate forecast of future droughts is significant for reliable planning and management of available water resources across this area.

## 5. Probabilistic drought forecasting

### a. Standardized streamflow index *(SSI)*

Droughts may occur in different phases of hydrologic cycle. Water movement through the hydrologic cycle is generally slow phenomenon, except for quick mass-transfer events like sudden storms. It happens that in a specific time window some hydrologic variables (e.g., soil moisture) experience a level of drought while some others (e.g., streamflow or water availability) do not undergo any identifiable drought categories. Therefore, drought status would appear differently upon the target hydrologic variable. To assess the drought status of a region, many drought indicators have been developed each using different hydrologic variables. The Palmer drought severity index (PDSI; Palmer 1965), SPI (McKee et al. 1993), crop moisture index (CMI; Palmer 1968), surface water supply index (SWSI; Shafer and Dezman 1982), and vegetation condition index (VCI; Liu and Kogan 1996) are among the most applied indices to characterize different drought types. Though all of these indices are widely used, each one focuses on particular hydrologic variables and has its own specific strengths and weaknesses.

*m*of the year

*y*for the time window of length

*k*(

Drought classification used by the U.S. Drought Monitor.

According to Eq. (11), separate distributions fit the aggregated flows with different starting months. Therefore, this definition of SSI preserves the seasonality effect; otherwise, the seasonal flow pattern would be disregarded if the marginal distribution fitted to the entire series of the aggregated flows.

### b. Analyses of historical droughts

According to Eq. (11), probabilistic drought analysis starts by fitting the marginal distributions to aggregated flows. For the seasonal drought forecast, monthly flow volumes are aggregated starting from each month over a sequence of 3 months (*p* value associated with the K–S test for different distributions fitted to the seasonal flow volumes in the training period of 1950–90. The best distributions, with the smallest AIC and the *p* values greater than the significance level

The AIC and the *p* value of the K–S test found for different distributions fitted to the seasonal flow volumes starting from different months over a sequence of 3 months. The numbers in the parentheses are the *p* values. The statistics of the best distributions are presented in bold.

Histogram against the best-fitted distribution of the aggregated flow volumes over a yearly sequence of 3 months during the training period of 1950–90.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Histogram against the best-fitted distribution of the aggregated flow volumes over a yearly sequence of 3 months during the training period of 1950–90.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Histogram against the best-fitted distribution of the aggregated flow volumes over a yearly sequence of 3 months during the training period of 1950–90.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

After fitting the marginal distributions, a backward calculation in the SSI formula [Eq. (11)] would give us the range of seasonal flow within each drought category used by the U.S. Drought Monitor (Table 1). The bar chart in Fig. 3 shows the range of flow volume within each drought category D0–D4. From this figure, the seasonal pattern is evidently captured as expected from the SSI definition in Eq. (11). The high-flow season shows more variability in the required amount of water to transition to a different category. In other words, the same class of drought is likely to persist within a high-flow season (e.g., the spring) if the seasonal flow does not change significantly. Moreover, while a particular seasonal flow (e.g., 90 KAF) is defined as an absolutely wet condition in winter (January–March or February–April periods), the same amount of flow might lead to an exceptional drought condition (D4) in spring (e.g., April–June and May–July periods). This is the seasonality issue reflected in drought definitions that identifies the high dependency of drought status on the time of the year.

Three-month seasonal flow volumes classified into different drought categories.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Three-month seasonal flow volumes classified into different drought categories.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Three-month seasonal flow volumes classified into different drought categories.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Using the marginal distributions found upon the aggregated flow during the training period (1950–90), the SSI with *k* = 3 [Eq. (11)] is calculated for the entire analysis period (1950–2011) as plotted in Fig. 4. Using the U.S. Drought Monitor categories shown in Table 1,

The SSI3 [SSI with *k* = 3 in Eq. (11)] time series during the analysis period. The line

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The SSI3 [SSI with *k* = 3 in Eq. (11)] time series during the analysis period. The line

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The SSI3 [SSI with *k* = 3 in Eq. (11)] time series during the analysis period. The line

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Matrix plot of seasonal droughts indicating the sequence of droughts during the analysis period.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Matrix plot of seasonal droughts indicating the sequence of droughts during the analysis period.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Matrix plot of seasonal droughts indicating the sequence of droughts during the analysis period.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

### c. Products of probabilistic drought forecast

For the forecast phase, this study selects the spring (April–June) as the target season. To determine the predictor seasons in conditional drought forecasting, the correlation between the margin of spring flow and the margins of prior seasons should be obtained (Fig. 6). Given the correlations in Fig. 6, the transformed spring flow (found from the marginal distribution of spring flow) has the highest correlation with the transformed winter flow (January–March) among the other seasons of the year. Yet the correlation between the seasonal flow of spring and fall or even summer is not insignificant, and the analysis of spring flow should be established upon either one season (winter) or two seasons (winter and fall) earlier. Therefore, two different Bayesian networks are applied to study the conditional probabilities of spring flow. In the first network, the spring flow is assumed to be only dependent on the winter flow and hence the drought forecasting is conducted within Eq. (9). The second network adds the impact of fall status on spring drought and it considers the seasonal flow of both past winter and fall in drought forecasting of spring season [Eq. (10)]. However, given the greater correlation between the spring and winter flow, it is expected that the winter influences the spring flow more than the fall.

The autocorrelation of the transformed spring flow (April–June) to the prior 3-month flows. The moving window locates no further than the previous year spring season (April–June). The circles are associated with the transformed seasonal flows before the spring (winter, fall, and summer).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The autocorrelation of the transformed spring flow (April–June) to the prior 3-month flows. The moving window locates no further than the previous year spring season (April–June). The circles are associated with the transformed seasonal flows before the spring (winter, fall, and summer).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The autocorrelation of the transformed spring flow (April–June) to the prior 3-month flows. The moving window locates no further than the previous year spring season (April–June). The circles are associated with the transformed seasonal flows before the spring (winter, fall, and summer).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

#### 1) First-order conditional forecasts

*t*copulas, from the elliptical family, are tested with parameters estimated by the inference function for margins (IFM) method (Joe 1997). Copula functions, and associated parameter estimation methods, are described in detail in the above-mentioned references. To select the appropriate copula among several alternatives, a visual inspection between the empirical and parametric copulas might be helpful (Madadgar and Moradkhani 2011); however, there are several goodness-of-fit (GOF) tests (Genest et al. 2009) that are more analytically conclusive in copula selection. The parametric bootstrapping GOF test, introduced by Genest and Rémillard (2008), is employed in this study to select the best copula fitted to the dataset. This GOF test calculates the Cramér–von Mises (Anderson 1962) statistic as a measure of distance between the empirical and parametric copulas:

*p*value is greater than a predefined significance level

*p*value (greater than the significance level as well) is selected as the best copula. In this study, the significance level is set to

*p*values are obtained by the parametric bootstrapping procedure with 1000 replications. As summarized in Table 3, the Gaussian copula with the smallest

*p*value is the best choice among others.

Statistics of the GOF test to find the best copula joining the marginal variables of seasonal flows in Eqs. (9) and (10)—*p* value are presented [Eq. (12)] in each case. The results for the superior copulas are shown in bold.

After picking the best-fitted copula for the data, several probabilistic analyses of drought status in spring season can be conducted using Eq. (9). In such an analysis, one might be interested in spring-flow distribution conditional upon a given winter drought status in winter. In this case, the winter drought conditions might be fixed at the very transition point of drought categories, where a drought status turns into another. Figure 7 shows the distribution of spring flow (April–June) conditional on the drought status in winter (January–March). Each curve represents the probability distribution function (PDF) associated with a particular drought status in winter (D0–D4). The analysis is narrowed down to the fixed drought conditions in winter. According to the drought PDFs, as the winter flow increases, which is equivalent to less intense drought in winter, the spring drought is expected to be less intense as well. For example, if a D4 drought occurred in winter, the spring drought status is likely to be more intense than if a D0 winter drought occurred. Moreover, when an intense drought is experienced in winter, the distribution of spring flow is rather narrow around its mode. For instance, the PDF associated with D0 winter droughts is wider than the PDF associated with the D4 winter drought. This leads to a larger range of spring flow given the D0 drought status in winter, as compared to the D4 winter drought.

Distribution of seasonal flow in the spring (Apr–Jun) given the drought status in winter. Each PDF is in accordance with a particular winter drought status.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Distribution of seasonal flow in the spring (Apr–Jun) given the drought status in winter. Each PDF is in accordance with a particular winter drought status.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Distribution of seasonal flow in the spring (Apr–Jun) given the drought status in winter. Each PDF is in accordance with a particular winter drought status.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

To expand the results and demonstrate the usefulness of this approach, the conditional probability of drought in spring season, given the drought condition in winter, is shown in a shaded scheme by Fig. 8. To show all the PDFs in a same plot, the probability distributions are scaled between 0 and 1; where 1 (the darkest area) represents the most probable spring flow (mode of the PDF), and hence drought condition, given the flow magnitude (i.e., hydrologic drought condition) in winter. Similarly, the areas with lighter shade represent the tails of the spring flow PDF given the winter flow. Therefore, the darker areas closely surround the mode of the PDFs; hence the associated range of spring flow is more likely to happen if the given winter flow is observed. The scatterplot of the spring flow against the winter flow during the entire analysis period, 1950–2011, is also shown in the *x–y* range of Fig. 8. The scatterplot is showing both training (1950–90) and validation (1990–2011) periods. As seen, the dark area (most likely situations given the winter flow) captures almost the entire scatterplot, which approves the reliable performance of the forecast methodology in both training and validation phases. The range of seasonal flows is split by the dashed lines, for both winter and spring, to illustrate the range of seasonal flows corresponding to various drought conditions, in each season. Figure 8 helps to find out the most likely drought status in spring, given the winter flow, by simply looking at the dark regions. Given D4 drought in winter, for instance, the spring drought status is most likely to be D4, D3, or D2 depending on the exact value of the winter flow. As another example, if a D2 drought occurs in winter, the spring drought status will likely be located in either D1 or D0 spans. Visual inspection of Fig. 7 also verifies that the PDF of the spring flow is narrower within the D1 and D0 spans, given the winter drought of type D2. Hence, Fig. 7 is a limited version of Fig. 8 where only a few conditional PDFs are illustrated.

Conditional PDF of spring (Apr–Jun) flow given the winter (Jan–Mar) flow. The conditional PDFs are scaled between 0 and 1 for visualization purposes. The shade level in each pixel indicates its conditional probability density. The circles show the spring-flow volume against the winter-flow volume during the analysis period of 1950–2011. The faint dashed vertical (horizontal) lines bounded by diamonds identify the range of winter (spring) seasonal flow within each drought category.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Conditional PDF of spring (Apr–Jun) flow given the winter (Jan–Mar) flow. The conditional PDFs are scaled between 0 and 1 for visualization purposes. The shade level in each pixel indicates its conditional probability density. The circles show the spring-flow volume against the winter-flow volume during the analysis period of 1950–2011. The faint dashed vertical (horizontal) lines bounded by diamonds identify the range of winter (spring) seasonal flow within each drought category.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Conditional PDF of spring (Apr–Jun) flow given the winter (Jan–Mar) flow. The conditional PDFs are scaled between 0 and 1 for visualization purposes. The shade level in each pixel indicates its conditional probability density. The circles show the spring-flow volume against the winter-flow volume during the analysis period of 1950–2011. The faint dashed vertical (horizontal) lines bounded by diamonds identify the range of winter (spring) seasonal flow within each drought category.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

To compare the forecast results of the presented method with ESP, Fig. 9 is developed showing the 90% predictive uncertainty bound of the retrospective forecasts within the validation period (1990–2011). As seen, there are only few spring seasons with drought conditions during 1990–2011. The uncertainty bound for the copula method is the limited representation of Fig. 8 where, given any winter flow, only the 5% and 95% bounds of the PDFs are shown. Unlike the copula-based forecast, the uncertainty bound of the ESP forecasts do not change smoothly as the winter flow increases. The uncertainty bound gradually expands in the copula-based forecasts, while it does not follow a particular trend for the ESP approach. To ensure the validity of the generated ESPs, the PRMS simulations are plotted versus the observations in Fig. 10. As seen, PRMS performs quite reasonably in simulating the 3-monthly flows implying that the generated ESPs (shown in Fig. 9) are reliable to compare with copula-based forecasts. According to Fig. 9, the significance of the copula-based forecast model over the ESP approach is that the generated uncertainty bound is reasonably large to encompass the observations showing a discernible trend against the increase in winter flow. Moreover, according to the expected values of forecasts, the performance of copula-based method is seen to be slightly better than the ESP forecasts.

The spring (Apr–Jun) and winter (Jan –Mar) seasonal flow forecast using the developed method vs the ESP forecasts generated by PRMS. Flow forecasts are corresponding to the validation period (1990–2011).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The spring (Apr–Jun) and winter (Jan –Mar) seasonal flow forecast using the developed method vs the ESP forecasts generated by PRMS. Flow forecasts are corresponding to the validation period (1990–2011).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The spring (Apr–Jun) and winter (Jan –Mar) seasonal flow forecast using the developed method vs the ESP forecasts generated by PRMS. Flow forecasts are corresponding to the validation period (1990–2011).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Observations vs PRMS simulations of 3-monthly flows beginning from each calendar month during 1980–2011.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Observations vs PRMS simulations of 3-monthly flows beginning from each calendar month during 1980–2011.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Observations vs PRMS simulations of 3-monthly flows beginning from each calendar month during 1980–2011.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The conditional probability of spring (Apr–Jun) flow [given the winter flow (Jan –Mar)] exceeding particular threshold associated with the drought status of *D*_{i}.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The conditional probability of spring (Apr–Jun) flow [given the winter flow (Jan –Mar)] exceeding particular threshold associated with the drought status of *D*_{i}.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The conditional probability of spring (Apr–Jun) flow [given the winter flow (Jan –Mar)] exceeding particular threshold associated with the drought status of *D*_{i}.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

#### 2) Second-order conditional forecasts

Given the correlation of spring flow to the past winter and fall seasons (Fig. 6), the spring droughts might be analyzed using the second-order conditional probabilities [Eq. (10)]. In Eq. (10), the random variables

Equation (10) returns the second-order conditional PDF of spring flow given the winter and fall observations. Figure 12 shows the conditional PDFs of spring flow given various droughts in the past seasons. In Fig. 12a, the fall drought is fixed at D4 and the winter drought varies from D0 to D4. Similar to Fig. 7, the modes of PDFs move toward smaller spring flows as the winter drought becomes more intense. The very little difference between Fig. 12a and Fig. 7 indicates that the fall status does not have a noticeable impact on the spring drought. This fact is approved with Fig. 12b as well where, for a given winter drought (D4), even a big change in fall status (D0 to D4) does not make a significant change in spring status. As seen, the PDFs in Fig. 12b are clustered together and are quite similar to the PDF associated with D4 drought in Figs. 7 and 12a. This is evidence of a high influence of winter status on the spring drought. In Figs. 12a,b, all the PDFs are conditional on a particular drought (D4) fixed for either fall or winter. Figure 12c displays the effect of any possible variation in winter and fall status on the spring drought. The PDFs are associated with two situations (D0 and D4) for either fall or winter. While the fall status is fixed at D0 (circle markers), the spring drought would change rather significantly upon the magnitude of change in winter status. In contrast, if the winter status is fixed at a particular drought, for example, D0 (solid markers), the spring status does not considerably change with even a big change in fall status. Hence, the spring drought is found to be more sensitive to the winter status than the fall status, which approves practicing the reduced-dimension form of the forecast model (denoted as the first-order conditional forecast).

Distribution of spring (Apr–Jun) flow given the drought status of (a) *D*_{i} in winter and D4 in fall, (b) D4 in winter and *D*_{i} in fall, (c) D0 to D4 in winter and/or in fall.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Distribution of spring (Apr–Jun) flow given the drought status of (a) *D*_{i} in winter and D4 in fall, (b) D4 in winter and *D*_{i} in fall, (c) D0 to D4 in winter and/or in fall.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Distribution of spring (Apr–Jun) flow given the drought status of (a) *D*_{i} in winter and D4 in fall, (b) D4 in winter and *D*_{i} in fall, (c) D0 to D4 in winter and/or in fall.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Another outlook to the spring drought would be the exceedance probability of spring flow given the winter and fall observations. This is equivalent to the probabilities developed in Fig. 11, with the exception that both fall and winter flows are used and therefore the 2D plot of Fig. 11 turns into a 3D plot (Fig. 13). Given the seasonal flow of winter and fall, Fig. 13 illustrates the probability of spring drought being equally wet or wetter than a particular drought status. For example, the probability layer associated with D0 (the lowest one) gives the probability of spring flow exceeding the threshold for D0 drought. Similar to Fig. 11, as the winter or fall flow increases, the probability layers approach 1 and get close together. This is a valid observation since the high-flow winters and falls are usually followed by wet rather than dry springs. Furthermore, for a particular drought status, the probabilities change more quickly with winter-flow variations than fall-flow variations. Given a particular winter flow and variable fall flow, the probabilities do not change as significantly as if the fall flow was fixed and the winter flow varied. This is consistent with the outcomes of Fig. 12 where the winter flow was found to have stronger impact on the spring status than the fall flow.

The conditional probability of spring (Apr–Jun) flow [given the winter (Jan–Mar) and fall (Oct–Dec) flow] exceeding particular threshold associated with the drought status of *D*_{i}. Each layer represents the exceeding probability of a particular drought in spring.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The conditional probability of spring (Apr–Jun) flow [given the winter (Jan–Mar) and fall (Oct–Dec) flow] exceeding particular threshold associated with the drought status of *D*_{i}. Each layer represents the exceeding probability of a particular drought in spring.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The conditional probability of spring (Apr–Jun) flow [given the winter (Jan–Mar) and fall (Oct–Dec) flow] exceeding particular threshold associated with the drought status of *D*_{i}. Each layer represents the exceeding probability of a particular drought in spring.

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

The seasonal flow hindcast during the validation period (1990–2011) is shown in Fig. 14. This plot is the 3D version of Fig. 9, where the conditional forecast of spring flow is upon the past winter and fall seasonal flows. The mesh grids show the uncertainty bound of spring flow generated by copula model. Similar to Fig. 9, the uncertainty bound associated with the copula model gradually expands as the winter and/or fall flow increases and it also captures the observed spring flows during the validation period. Furthermore, the variation of the uncertainty bound is more dependent on the variation of winter flow than the fall flow, as earlier approved by Figs. 12 and 13. This is the reason that the hindcasts of Figs. 9 and 14 are very similar.

Performance of the copula model vs the ESP approach in spring flow hindcast given the observed flow of the past winter and fall. The mesh grids show the 5%–95% uncertainty bound of spring flow generated by copula model. The observations correspond to the validation period (1990–2011).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Performance of the copula model vs the ESP approach in spring flow hindcast given the observed flow of the past winter and fall. The mesh grids show the 5%–95% uncertainty bound of spring flow generated by copula model. The observations correspond to the validation period (1990–2011).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

Performance of the copula model vs the ESP approach in spring flow hindcast given the observed flow of the past winter and fall. The mesh grids show the 5%–95% uncertainty bound of spring flow generated by copula model. The observations correspond to the validation period (1990–2011).

Citation: Journal of Hydrometeorology 14, 6; 10.1175/JHM-D-13-010.1

## 6. Summary and conclusions

Hydrologic droughts across the Gunnison River basin, a subbasin of the upper Colorado River basin, were analyzed, and the probabilistic forecasts of seasonal droughts were obtained. The standardized streamflow index (SSI) was used to characterize drought events and was also adopted to address the seasonality effects in drought analysis. More fall and winter droughts occurred in the 1950s and 1960s, while the spring and summer droughts were distributed more evenly during the historical period (1950–2011). In the forecast phase, the Bayesian networks contributed to reflect the autocorrelation in the streamflow record and develop the conditional probabilities. An autocorrelation analysis of seasonal flows indicated that the spring flow (forecast variable) has the highest autocorrelation with the previous winter and the second highest correlation with the previous fall. A forecast model reflecting the dependency of seasonal flows was then developed by a family of multivariate distribution functions, the so-called copulas. Copulas are distribution functions that join several variables with different levels of correlation and dependence structures. Owing to the autocorrelation results, two different forecast models were established to obtain the first-order and second-order conditional probabilities of the spring flow. In the former (latter) case, the spring flow was estimated upon the observation of previous winter (winter and fall). The forecast products were presented as the conditional PDF and CDF of spring droughts.

In the first-order conditional model, it was seen that for a given drought status in winter, a less-severe drought is anticipated in spring. In addition, the probabilistic forecasts provided a tool to estimate the most probable status in future (here, spring), while the past status is known (e.g., winter). The probability of terminating the winter drought in the upcoming spring was also estimated for different drought states. A dry spell is less likely in spring if the winter experiences either normal or wet condition. In the second-order conditional model, the fall status displayed very slight impact on spring drought, which was in agreement with the results of autocorrelation analysis. For a given winter drought, even a big change in fall status could not make a significant change in spring status. As seen from the conditional CDFs, the 3D probability layers of spring flow approached 1 and got close together as the winter or fall flow were increasing. This fact indicated that the high-flow winters and falls are usually followed by wet rather than dry springs. According to forecasts, the spring drought was found to be more sensitive to the winter status than the fall status, which approves practicing the reduced-dimension form of the forecast model (denoted as the first-order conditional forecast).

The drought forecasts using the presented methodology were also compared with the conventional ESP forecasts. The ESP retrospective forecasts, starting from 1 January of each year from 1990 to 2011, were produced using the PRMS watershed model. The 90% predictive uncertainty bound of the spring drought forecasts generated by the copula-based method changed smoothly as the given winter and/or fall flow varied gradually. The uncertainty bound made by ESP forecasts did not follow any discernible trend, whereas copula-based forecasts showed a particular trend. The uncertainty bound slowly expands as the winter and or fall flow increases; however, the change is again much smaller in respect with fall flow than the winter flow.

The proposed probabilistic approach helps the decision makers to develop drought mitigation plans and policies with an appropriate insight toward the future drought status. Given the basics of copula modeling, the proposed methodology shows promising results for adequately correlated and dependent variables. However, correlation and dependency are the essential components of the *conditional* forecast methodologies, whereas without dependent variables, the conditional forecast would not be meaningful. According to the potential of the proposed methodology in probabilistic forecast of the future drought status, the next attempt will study the spatial variability of future drought across the region. Reported results in this study are based upon the streamflow record at the basin outlet; however, the spatial variability of dry spells throughout a region using this method may be of interest and this is subject to future studies. An alternative probabilistic approach to the presented method is the use of data assimilation where the uncertainties from observation, initial condition, and model can be integrated into the forecasting framework (Moradkhani et al. 2012; DeChant and Moradkhani 2012; Parrish et al. 2012).

## Acknowledgments

Partial financial support for this study was provided by NOAA-MAPP Grant NA110AR4310140.

## REFERENCES

Akaike, H., 1974: A new look at the statistical model identification.

,*IEEE Trans. Autom. Control***19**, 716–723, doi:10.1109/TAC.1974.1100705.Anderson, T. W., 1962: On the distribution of the two-sample Cramer–von Mises criterion.

,*Ann. Math. Stat.***33**, 1148–1159, doi:10.1214/aoms/1177704477.Andreadis, K. M., and Lettenmaier D. P. , 2006: Trends in 20th century over the continental United States.

,*Geophys. Res. Lett.***33**, L10403, doi:10.1029/2006GL025711.Bárdossy, A., 2006: Copula-based geostatistical models for groundwater quality parameters.

,*Water Resour. Res.***42**, W11416, doi:10.1029/2005WR004754.Bárdossy, A., and Li J. , 2008: Geostatistical interpolation using copulas.

,*Water Resour. Res.***44**, W07412, doi:10.1029/2007WR006115.Barros, A. P., and Bowden G. J. , 2008: Toward long-lead operational forecasts of drought: An experimental study in the Murray-Darling River Basin.

,*J. Hydrol.***357**, 349–367, doi:10.1016/j.jhydrol.2008.05.026.Blenkinsop, S., and Fowler H. J. , 2007: Changes in drought frequency, severity and duration for the British Isles projected by the PRUDENCE regional climate models.

,*J. Hydrol.***342**(1–2), 50–71, doi:10.1016/j.jhydrol.2007.05.003.Brown, J. D., and Seo D. J. , 2010: A nonparametric postprocessor for bias correction of hydrometeorological and hydrologic ensemble forecasts.

,*J. Hydrometeor.***11**, 642–665.Burke, E. J., Perry R. H. J. , and Brown S. J. , 2010: An extreme value analysis of UK drought and projections of change in the future.

,*J. Hydrol.***388**, 131–143, doi:10.1016/j.jhydrol.2010.04.035.Cancelliere, A., Mauro G. D. , Bonaccorso B. , and Rossi G. , 2007: Drought forecasting using the standardized precipitation index.

,*Water Resour. Manag.***21**, 801–819.Carbone, G. J., and Dow K. , 2005: Water resource management and drought forecasts in South Carolina.

,*J. Amer. Water Resour. Assoc.***41**, 145–155, doi:10.1111/j.1752-1688.2005.tb03724.x.Day, G. N., 1985: Extended streamflow forecasting using NWSRFS.

,*J. Water Resour. Plan. Manage.***111**(2), 157–170, doi:10.1061/(ASCE)0733-9496(1985)111:2(157).DeChant, C., and Moradkhani H. , 2011: Improving the characterization of initial condition for ensemble streamflow prediction using data assimilation.

,*Hydrol. Earth Syst. Sci.***15**, 3399–3410, doi:10.5194/hess-15-3399-2011.DeChant, C., and Moradkhani H. , 2012: Examining the effectiveness and robustness of data assimilation methods for calibration and quantification of uncertainty in hydrologic forecasting.

,*Water Resour. Res.***48**, W04518, doi:10.1029/2011WR011011.Dupuis, D. J., 2007: Using copulas in hydrology: Benefits, cautions, and issues.

,*J. Hydrol. Eng.***12**, 381–393, doi:10.1061/(ASCE)1084-0699(2007)12:4(381).Embrechts, P., Lindskog F. , and McNeil A. J. , 2003: Modelling dependence with copulas and applications to risk management.

*Handbook of Heavy Tailed Distributions in Finance,*S. T. Rachev, Ed., Elsevier Science, 329–384.Favre, A. C., Adlouni S. E. , Perreault L. , Thiéonge N. , and Bobée B. , 2004: Multivariate hydrological frequency analysis using copulas.

,*Water Resour. Res.***40**, W01101, doi:10.1029/2003WR002456.Genest, C., and Rémillard B. , 2008: Validity of the parametric bootstrap for goodness-of-fit testing in semiparametric models.

,*Annal. Inst. Henri-Poincaré***44**, 1096–1127.Genest, C., Rémillard B. , and Beaudoin D. , 2009: Goodness-of-fit tests for copulas: A review and a power study.

,*Insur. Math. Econ.***44**, 199–213, doi:10.1016/j.insmatheco.2007.10.005.Ghosh, S., and Mujumdar P. P. , 2007: Nonparametric methods for modeling GCM and scenario uncertainty in drought assessment.

,*Water Resour. Res.***43**, W07405, doi:10.1029/2006WR005351.Halmstad, A., Najafi M. R. , and Moradkhani H. , 2012: Analysis of precipitation extremes with the assessment of regional climate models over the Willamette River Basin, USA.

*Hydrol. Processes,***27,**2579–2590, doi:10.1002/hyp.9376.Huang, W. C., and Chou C. C. , 2005: Drought early warning system in reservoir operation: Theory and practice.

,*Water Resour. Res.***41**, W11406, doi:10.1029/2004WR003830.Huang, W. C., and Chou C. C. , 2008: Risk-based drought early warning system in reservoir operation.

,*Adv. Water Resour.***31**, 649–660, doi:10.1016/j.advwatres.2007.12.004.Hwang, Y., and Carbone G. J. , 2009: Ensemble forecasts of drought indices using a conditional residual resampling technique.

,*J. Appl. Meteor. Climatol.***48**, 1289–1301, doi:10.1175/2009JAMC2071.1.Jaeger, W. K., and Coauthors, 2013: Toward a formal definition of water scarcity in natural-human systems.

,*Water Resour. Res.***49,**4506–4517, doi:10.1002/wrcr.20249.Joe, H., 1997:

*Multivariate Models and Dependence Concepts.*Chapman & Hall, 399 pp.Kao, S., and Govindaraju R. S. , 2008: Trivariate statistical analysis of extreme rainfall events via Plackett family of copulas.

,*Water Resour. Res.***44**, W02415, doi:10.1029/2007WR006261.Kao, S., and Govindaraju R. S. , 2010: A copula-based joint deficit index for droughts.

,*J. Hydrol.***380**, 121–134, doi:10.1016/j.jhydrol.2009.10.029.Karl, T., Quinlan F. , and Ezell D. S. , 1987: Drought termination and amelioration: Its climatological probability.

,*J. Climate Appl. Meteor.***26**, 1198–1209, doi:10.1175/1520-0450(1987)026<1198:DTAAIC>2.0.CO;2.Kendall, J. R., and Dracup D. A. , 1992: On the generation of drought events using an alternating renewal-reward model.

,*Stochastic Hydrol. Hydraul.***6**, 55–68, doi:10.1007/BF01581675.Kolmogorov, A. N., 1933: Sulla Determinazione Empirca di una Legge di Distribuzione, Giornale dell’.

,*Istituto Italiano degli Attuari***4**, 83–91.Leavesley, G. H., Lichty R. W. , Troutman B. M. , and Saindon L. G. , 1983: Precipitation-runoff modeling system: User’s manual. U.S. Geological Survey Water-Resources Investigations Rep. 83-4238, 207 pp.

Li., C., Singh V. P. , and Mishra A. K. , 2013: Monthly river flow simulation with a joint conditional density estimation network.

*Water Resour. Res.,***49,**3229–3242, doi:10.1002/wrcr.20146.Liu, W. T., and Kogan F. N. , 1996: Monitoring regional drought using the vegetation condition index.

,*Int. J. Remote Sens.***17**, 2761–2782, doi:10.1080/01431169608949106.Loaiciga, H. A., and Leipnik R. B. , 1996: Stochastic renewal model of low-flow streamflow sequences.

,*Stochastic Hydrol. Hydraul.***10**, 65–85, doi:10.1007/BF01581794.Lohani, V. K., and Loganathan G. V. , 1997: An early warning system for drought management using the palmer drought index.

,*J. Amer. Water Resour. Assoc.***33**, 1375–1386, doi:10.1111/j.1752-1688.1997.tb03560.x.Lott, N., and Ross T. , 2000: A climatology of recent extreme weather and climate events. NCDC Tech. Rep. 2000-02, 17 pp.

Lott, N., and Ross T. , 2006: Tracking and evaluating U.S. billion dollar weather disasters, 1980–2005.

*Extended Abstracts, AMS Forum: Environmental Risk and Impacts on Society: Successes and Challenges,*Atlanta, GA, Amer. Meteor. Soc., 1.2. [Available online at https://ams.confex.com/ams/Annual2006/techprogram/paper_100686.htm.]Madadgar, S., and Moradkhani H. , 2011: Drought analysis under climate change using copula.

,*J. Hydrol. Eng.***18**(7), 746–759, doi:10.1061/(ASCE)HE.1943-5584.0000532.Madadgar, S., Moradkhani, H., and Garen D. , 2013: Towards improved post-processing of hydrologic forecast ensembles.

*Hydrol. Processes,*doi:10.1002/hyp.9562, in press.Massey, F. J., 1951: The Kolmogorov-Smirnov test for goodness of fit.

,*Amer. Stat. Assoc.***46**(253), 68–78, doi:10.1080/01621459.1951.10500769.McKee, T. B., Doesken N. J. , and Kleist J. , 1993: The relationship of drought frequency and duration to time scales. Preprints,

*Eighth Conf. on Applied Climatology,*Anaheim, CA, Amer. Meteor. Soc., 179–184.Mishra, A. K., and Desai V. R. , 2005: Drought forecasting using stochastic models.

,*Stochastic Environ. Res. Risk Assess.***19**(5), 326–339, doi:10.1007/s00477-005-0238-4.Mishra, A. K., and Desai V. R. , 2006: Drought forecasting using feed-forward recursive neural network.

,*Ecol. Modell.***198**, 127–138, doi:10.1016/j.ecolmodel.2006.04.017.Mishra, A. K., and Singh V. P. , 2010: A review of drought concepts.

,*J. Hydrol.***391**, 202–216, doi:10.1016/j.jhydrol.2010.07.012.Moradkhani, H., and Meier M. , 2010: Long-lead water supply forecast using large-scale climate predictors and independent component Analysis.

,*J. Hydrol. Eng.***15**, 744–762.Moradkhani, H., Baird R. G. , and Wherry S. , 2010: Assessment of climate change impact on floodplain and hydrologic ecotones.

,*J. Hydrol.***395**, 264–278, doi:10.1016/j.jhydrol.2010.10.038.Moradkhani, H., DeChant C. M. , and Sorooshian S. , 2012: Evolution of ensemble data assimilation for uncertainty quantification using the particle filter–Markov Chain Monte Carlo method.

,*Water Resour. Res.***48**, W12520, doi:10.1029/2012WR012144.Najafi, M. R., Moradkhani H. , and Piechota T. , 2012: Ensemble streamflow prediction: Climate signal weighting vs. climate forecast system reanalysis.

,*J. Hydrol.***442–443**, 105–116, doi:10.1016/j.jhydrol.2012.04.003.Nelsen, R. B., 1999:

*An Introduction to Copulas*. Springer, 216 pp.Niemeyer, S., 2008: New drought indices.

*Proc. First Int. Conf. on Drought Management: Scientific and Technological Innovations,*Zaragoza, Spain, CIHEAM, 267–274.Özger, M., Mishra A. K. , and Singh V. P. , 2012: Long lead time drought forecasting using a wavelet and fuzzy logic combination model: A case study in Texas.

,*J. Hydrometeor.***13**, 284–297, doi:10.1175/JHM-D-10-05007.1.Palmer, W. C., 1965: Meteorologic drought. U.S. Department of Commerce, Weather Bureau, Research Paper 45, 58 pp.

Palmer, W. C., 1968: Keeping track of crop moisture conditions, nationwide: The new crop moisture index.

,*Weatherwise***21**, 156–161, doi:10.1080/00431672.1968.9932814.Parrish, M., Moradkhani H. , and DeChant C. M. , 2012: Towards reduction of model uncertainty: Integration of Bayesian model averaging and data assimilation.

,*Water Resour. Res.***48**, W03519, doi:10.1029/2011WR011116.Risley, J., Moradkhani H. , Hay L. , and Markstrom S. , 2011: Statistical comparisons of watershed-scale response to climate change in selected basins across the United States.

,*Earth Interact.***15**, 1–26, doi:10.1175/2010EI364.1.Russell, S. J., and Peter N. , 2009:

*Artificial Intelligence: A Modern Approach.*3rd ed. Prentice Hall, 1152 pp.Salvadori, G., and De Michele C. , 2006: Statistical characterization of temporal structure of storms.

,*Adv. Water Resour.***29**, 827–842, doi:10.1016/j.advwatres.2005.07.013.Salvadori, G., and De Michele C. , 2010: Multivariate multiparameter extreme value models and return periods: A copula approach.

,*Water Resour. Res.***46**, W10501, doi:10.1029/2009WR009040.Schaake, J. C., and Coauthors, 2007: Precipitation and temperature ensemble forecasts from single-value forecasts.

,*Hydrol. Earth Syst. Sci. Discuss.***4**, 655–717.Shafer, B. A., and Dezman L. E. , 1982: Development of a Surface Water Supply Index (SWSI) to assess the severity of drought conditions in snowpack runoff areas. Preprints,

*50th Annual Western Snow Conf.,*Reno, NV, Colorado State University, 164–175.Sheffield, J., and Wood E. F. , 2007: Characteristics of global and regional drought, 1950–2000: Analysis of soil moisture data from off-line simulation of the terrestrial hydrologic cycle.

,*J. Geophys. Res.***112**, D17115, doi:10.1029/2006JD008288.Sheffield, J., and Wood E. F. , 2008: Projected changes in drought occurrence under future global warming from multi-model, multi-scenario, IPCC AR4 simulations.

,*Climate Dyn.***31**, 79–105, doi:10.1007/s00382-007-0340-z.Shiau, J. T., 2006: Fitting drought duration and severity with two-dimensional copulas.

,*Water Resour. Manage.***20**, 795–815, doi:10.1007/s11269-005-9008-9.Showstack, R., 2012: Drought research and monitoring program is focus of congressional hearing.

,*Eos, Trans. Amer. Geophys. Union***93**(32), 310, doi:10.1029/2012EO320002.Shukla, S., Steinemann A. C. , and Lettenmaier D. P. , 2011: Drought monitoring for Washington State: Indicators and applications.

,*J. Hydrometeor.***12**, 66–83, doi:10.1175/2010JHM1307.1.Sklar, K., 1959: Fonctions de repartition à n Dimensions et Leura Marges.

*Publ. Inst. Stat. Univ. Paris,***8,**229–231.Steinemann, A., 2003: Drought indicators and triggers: A stochastic approach to evaluation.

,*J. Amer. Water Resour. Assoc.***39**, 1217–1233, doi:10.1111/j.1752-1688.2003.tb03704.x.Steinemann, A., 2006: Using climate forecasts for drought management.

,*J. Appl. Meteor. Climatol.***45**, 1353–1361, doi:10.1175/JAM2401.1.Thulasiraman, K., and Swamy M. N. S. , 1992:

*Graphs: Theory and Algorithms.*John Wiley and Son, 460 pp.Twedt, T. M., Schaake J. C. , and Peck E. L. , 1977: National Weather Service extended streamflow prediction.

*Proc. 45th Western Snow Conf.,*Albuquerque, NM, 52–57.van Huijgevoort, M. H. J., Hazenberg P. , van Lanen H. A. J. , and Uijlenhoet R. , 2012: A generic method for hydrological drought identification across different climate regions.

,*Hydrol. Earth Syst. Sci.***16**, 2437–2451, doi:10.5194/hess-16-2437-2012.Van Loon, A. F., and Van Lanen H. A. J. , 2012: A process-based typology of hydrological drought.

,*Hydrol. Earth Syst. Sci.***16**, 1915–1946, doi:10.5194/hess-16-1915-2012.Wong, G., Lambert M. F. , Leonard M. , and Metcalfe A. V. , 2010: Drought analysis using trivariate copulas conditional on climatic states.

,*J. Hydrol. Eng.***15**(2), 129–141, doi:10.1061/(ASCE)HE.1943-5584.0000169.Wood, A. W., and Lettenmaier D. P. , 2008: An ensemble approach for attribution of hydrologic prediction uncertainty.

,*Geophys. Res. Lett.***35**, L14401, doi:10.1029/2008GL034648.WWA and NIDIS, 2012: The 2012 Drought in Colorado, Utah and Wyoming: A July 2012 update from the Western Water Assessment and the National Integrated Drought Information System, WWA and NIDIS, 2 pp. [Available online at http://www.drought.gov/imageserver/NIDIS/DEWS/UCRB/docs/WWA-NIDIS_July_2012_Drought_update.pdf.]

Zhang, L., and Singh V. P. , 2007a: Trivariate flood frequency analysis using the Gumbel–Hougaard copula.

,*J. Hydrol. Eng.***12**, 431–439, doi:10.1061/(ASCE)1084-0699(2007)12:4(431).Zhang, L., and Singh V. P. , 2007b: Gumbel–Hougaard copula for trivariate rainfall frequency analysis.

,*J. Hydrol. Eng.***12**, 409–419, doi:10.1061/(ASCE)1084-0699(2007)12:4(409).