Bivariate Drought Analysis Using Streamflow Reconstruction with Tree Ring Indices in the Sacramento Basin , California , USA

Long-term streamflow data are vital for analysis of hydrological droughts. Using an artificial neural network (ANN) model and nine tree-ring indices, this study reconstructed the annual streamflow of the Sacramento River for the period from 1560 to 1871. Using the reconstructed streamflow data, the copula method was used for bivariate drought analysis, deriving a hydrological drought return period plot for the Sacramento River basin. Results showed strong correlation among drought characteristics, and the drought with a 20-year return period (17.2 million acre-feet (MAF) per year) in the Sacramento River basin could be considered a critical level of drought for water shortages.


Introduction
Generally, long-term data are recommended for analyzing floods and droughts.However, although precipitation data are available from the 16th century onwards, their quality and reliability in many countries are questionable because of the methods of observation, different periods of observations, uncertainties associated with gaging sites, and temporal resolution of observations [1,2].Many studies have, therefore, used tree-ring data as a way to acquire data for longer periods of time, up to 500 years [3].
Since Ferguson [4] correlated observed hydro-meteorological data and tree-ring data in California, many studies have utilized tree-ring data to reconstruct time series of the past.Fritts [5] did tree ring analysis and correlated the data with climate, suggesting that it can be used in water resources management [6].Also, some studies have correlated hydro-meteorological variables with tree-ring data to reconstruct climate factors [7][8][9][10].Some studies have reconstructed seasonal series, such as precipitation [11,12], natural hazards [13], and temperature [14], based on tree ring width data.However, most of the studies on tree rings and hydrological phenomena have focused on droughts [15][16][17][18][19][20][21][22][23][24].This may be because annual data are normally reconstructed using tree ring data and hence it is difficult to analyze intra-annual hydrological phenomena.Long time-scale data are Water 2016, 8, 122 2 of 16 mainly used for droughts, because it is difficult to determine the onset or end of a drought and often the drought may last for several months or years.The long time-scale occasionally gives rise to a sample size problem [25] for drought analysis that can be overcome with long-term tree ring reconstruction [26].Some studies have directly correlated tree rings with droughts [15,27], drought patterns or oscillations [28][29][30], drought index and its trend [23,[31][32][33], return periods [34], and spatial drought characteristics [35].
Bivariate (or multivariate) analyses of drought characteristics, such as severity, duration, and arrival time, are being increasingly made [36].These analyses have introduced multivariate drought indices, such as Multivariate Drought Index (MDI) including precipitation, runoff, evapotranspiration, and soil moisture [37]; Multivariate Standardized Drought Index (MDSI) which combines the Standardized Precipitation Index and the Standardized Soil Moisture Index [38]; and Vegetation Drought Response Index (VegDRI) which integrates climatic indicators; and satellite-derived vegetation index [39,40].Some studies have employed conventional multivariate analysis for drought indices with PDSI [23,41] and SWSI [42,43], or bivariate frequency analysis [44,45] assuming that all variables had the same probability distribution.To overcome this restriction, the copula method has been developed [46].For doing bivariate drought analysis by the copula method [46,47], tree ring reconstruction can be employed to our advantage.
The objective of this study, therefore, is to reconstruct the annual streamflow of the Sacramento River in California and Oregon using tree-ring width data, and use the reconstructed data for bivariate drought frequency analysis with the copula method.Selected tree-ring data were used in an artificial neural network for streamflow reconstruction and the reconstructed data were verified by comparing with actual observations.The Archimedean copula function was applied to the reconstructed streamflow data and then the return period plot of the hydrological drought in the Sacramento River basin was derived.The advantages of using the copula method have been discussed in many studies [34,36,[48][49][50][51][52][53][54][55].Since drought may last from months to years, as has happened in California, long-term reconstruction based on tree-ring data, addresses the drawback of short-term data [56].

Study Area and Data
Annual streamflow data for four Sacramento rivers and tree-ring data of the nearby region were employed in this study.The tree-ring data from 17 sites (Figure 1 and Table 1) in California and Oregon, which reflect standard chronologies of ring width [57], were obtained from the International Tree-Ring Databank [58].The tree ring width data were standardized [59].Annual streamflow is the sum of four river flows, which are the Sacramento River above Bend Bridge (SBB), the Feather River at the Lake Oroville (FTO), the Yuba River near Smartville (YRS) and the American River at Folsom (AMF), which was obtained from the California Data Exchange Center of California Department of Water Resources [60] for the period 1872 to 1977.It has a long-term mean of 18.9 MAF (million acre-feet; 1.23 ˆ10 9 m 3 ), median of 17.6 MAF, and maximum flow of 51.6 MAF that occurred in the year 1890.These tree ring and streamflow sites are shown in Figure 1 and Table 1.

Drought Definition Using the Run Theory
This study qualitatively defines the hydrological drought (hereafter referred to as "drought") as a significant decrease in the availability of streamflow in the river.Quantitatively, drought was defined using the run theory, which allows us to calculate drought duration, severity, and arrival time [61,62].Thus, drought can be defined as the time when a hydro-meteorological time series x t falls below the truncation level x 0 and that represents a hydro-meteorological event (Figure 2).Drought events are based on the truncation level, so the selection of the level is one of the important issues for proper drought analysis.Generally, the mean value of streamflow has been widely used as the truncation level [63][64][65][66][67][68].However, the Sacramento River streamflow shows high variability, between 5.13 and 51.65 MAF, so a median value of annual streamflow was regarded as a more reliable truncation level in this study.

Artificial Neural Network
The factors that influence annual hydro-meteorological behavior can be roughly classified into four groups: (i) atmospheric-climatic; (ii) geologic-geomorphic; (iii) soil-vegetation; and (iv) runoff-channel factors [69].Tree-ring widths in a trunk of a tree are also influenced by atmospheric-climatic and soil-vegetation factors, such as precipitation, evapotranspiration, and soil moisture.This indicates that an appropriate modeling technique and tree rings that have a correlation with atmospheric-climatic factors can be used to reconstruct annual streamflow, and this study employed an Artificial Neural Network (ANN) model.ANN mimics the structure and functions of a biological neural system, in which neurons are connected through nodes [70].After the perceptron was proposed to categorize information patterns [71], ANNs have been widely used to recognize nonlinear relationships between different variables.The ANN used in this study was comprised of three layers: the input layer that represents observed streamflow data, the output layer that produces simulated streamflow, and the hidden layer that is constituted by a network of neurons that are trained to recognize patterns from observations (Figure 3).Drought events are based on the truncation level, so the selection of the level is one of the important issues for proper drought analysis.Generally, the mean value of streamflow has been widely used as the truncation level [63][64][65][66][67][68].However, the Sacramento River streamflow shows high variability, between 5.13 and 51.65 MAF, so a median value of annual streamflow was regarded as a more reliable truncation level in this study.

Artificial Neural Network
The factors that influence annual hydro-meteorological behavior can be roughly classified into four groups: (i) atmospheric-climatic; (ii) geologic-geomorphic; (iii) soil-vegetation; and (iv) runoff-channel factors [69].Tree-ring widths in a trunk of a tree are also influenced by atmospheric-climatic and soil-vegetation factors, such as precipitation, evapotranspiration, and soil moisture.This indicates that an appropriate modeling technique and tree rings that have a correlation with atmospheric-climatic factors can be used to reconstruct annual streamflow, and this study employed an Artificial Neural Network (ANN) model.ANN mimics the structure and functions of a biological neural system, in which neurons are connected through nodes [70].After the perceptron was proposed to categorize information patterns [71], ANNs have been widely used to recognize nonlinear relationships between different variables.The ANN used in this study was comprised of three layers: the input layer that represents observed streamflow data, the output layer that produces simulated streamflow, and the hidden layer that is constituted by a network of neurons that are trained to recognize patterns from observations (Figure 3).Drought events are based on the truncation level, so the selection of the level is one of the important issues for proper drought analysis.Generally, the mean value of streamflow has been widely used as the truncation level [63][64][65][66][67][68].However, the Sacramento River streamflow shows high variability, between 5.13 and 51.65 MAF, so a median value of annual streamflow was regarded as a more reliable truncation level in this study.

Artificial Neural Network
The factors that influence annual hydro-meteorological behavior can be roughly classified into four groups: (i) atmospheric-climatic; (ii) geologic-geomorphic; (iii) soil-vegetation; and (iv) runoff-channel factors [69].Tree-ring widths in a trunk of a tree are also influenced by atmospheric-climatic and soil-vegetation factors, such as precipitation, evapotranspiration, and soil moisture.This indicates that an appropriate modeling technique and tree rings that have a correlation with atmospheric-climatic factors can be used to reconstruct annual streamflow, and this study employed an Artificial Neural Network (ANN) model.ANN mimics the structure and functions of a biological neural system, in which neurons are connected through nodes [70].After the perceptron was proposed to categorize information patterns [71], ANNs have been widely used to recognize nonlinear relationships between different variables.The ANN used in this study was comprised of three layers: the input layer that represents observed streamflow data, the output layer that produces simulated streamflow, and the hidden layer that is constituted by a network of neurons that are trained to recognize patterns from observations (Figure 3).The back-propagation algorithm was used to train the network through the adjustment of connection strength to learn about the error and optimize the neurons.It calculates the error function with respect to all the weights in the network and the gradient of error functions is fed into the optimization technique, which attempts to minimize the error of the network.Hence, the selection of back-propagation algorithm is one of the challenges when using a neural network [72].The Levenberg-Marquardt-QNBP algorithm was selected as the back-propagation algorithm, because it is known to work well for non-linear problems, such as those related to meteorological and hydrological data [73].Also, the number of hidden layers, which can be optimized by a trial and error method, is important for a proper network.The ANN model used in this study has six hidden layers that are optimized.One of the advantages of ANN is that it can be used as an alternative modeling technique when the data show non-linearity, which may cause error with a linear technique [74].The tree-ring data in California and Oregon have autocorrelation and lagged-correlation characteristics [75], so ANN was employed as an alternative to reconstruct streamflow using tree-ring data.More details on ANN and the back-propagation algorithm can be found in Basheer and Hajmeer [76].

Drought Frequency Analysis Based on Copula
Unlike precipitation or flood occurrence, drought shows a different statistical behavior for a different duration [62].Considering drought duration and severity as mutually related variables, the copula method has been employed to capture the dependence between them [47,77].For a probability distribution F px 1 , . . ., x n q, which has n-dimensional marginal distributions F 1 px 1 q , . . ., F n p x n q, the copula function C that satisfies the relationship between marginal variables can be expressed as: 1{T where ψ α denotes the generating function; ψ α ´1 is the pseudo-inverse of that function, which differs with the copula family; T denotes the return period; and E(L) is the interval between events.Unlike univariate frequency analysis, bivariate frequency indicates the probability that the phenomenon under study occurs if and only if a prior condition takes place.There are several types of copula functions, but the Archimedean copula family, which allows for greater flexibility and simplicity of use, is more commonly used in hydrology [78].From the Archimedean copula family, the Clayton, Gumbel, and Frank copulas were employed and their functions are given in Table 2.
Table 2. Bivariate Archimedean copula family: C is the copula function, t denotes the drought event, α is the copula parameter, and F 1 and F 2 denote cumulative distribution function of each variable [79].

Tree-Ring Data Screening
Appropriate tree rings, which have correlation with atmospheric-climatic factors, can be used as the predictor for annual streamflow.Therefore, selection of the appropriate input for the ANN model was one of the challenges in this study.Generally, a trial and error method with different input variables is employed, but it can lead to the poor performance of neural networks [80,81].Alternatively, Water 2016, 8, 122 6 of 16 cross correlation coefficients were employed to select appropriate inputs in this study in the same way as Meko et al. [75].
Seven tree-rings indices showed correlation with streamflow (Figure 4), including Antelope Lake (Pinus Jeffreyi and Pinus Ponderosa), Felkner Ridge, Frederick Butte, Lemon Canyon, Piute Mountain, and Sorrel Peak; four tree rings had correlation with a one-year time-lag, including Dalton Reservoir, Hager Basin Reservoir and Antelope Lake (Pinus Jeffreyi and Pinus Ponderosa), as shown in Figures 4 and 5. Hence, nine tree rings were selected as predictors for the ANN model (two tree-rings were overlapped in zero-lagged and one-year lagged).Also, four tree-rings, which had one-year time-lag correlation (Antelope with Pinus Jeffreyi and Pinus Ponderosa, Dalton Reservoir, and Hager Basin Reservoir), are located on the nearby lake or reservoir.Therefore, it seems that the groundwater level or soil moisture influenced tree ring width, but there are no clues to estimate the correlation between them and further studies are thus needed.These nine tree-ring data points composed the input dataset for the ANN model.Alternatively, cross correlation coefficients were employed to select appropriate inputs in this study in the same way as Meko et al. [75].
Seven tree-rings indices showed correlation with streamflow (Figure 4), including Antelope Lake (Pinus Jeffreyi and Pinus Ponderosa), Felkner Ridge, Frederick Butte, Lemon Canyon, Piute Mountain, and Sorrel Peak; four tree rings had correlation with a one-year time-lag, including Dalton Reservoir, Hager Basin Reservoir and Antelope Lake (Pinus Jeffreyi and Pinus Ponderosa), as shown in Figures 4 and 5. Hence, nine tree rings were selected as predictors for the ANN model (two tree-rings were overlapped in zero-lagged and one-year lagged).Also, four tree-rings, which had one-year time-lag correlation (Antelope with Pinus Jeffreyi and Pinus Ponderosa, Dalton Reservoir, and Hager Basin Reservoir), are located on the nearby lake or reservoir.Therefore, it seems that the groundwater level or soil moisture influenced tree ring width, but there are no clues to estimate the correlation between them and further studies are thus needed.These nine tree-ring data points composed the input dataset for the ANN model.

Reconstructed ANN Model Calibration and Validation
The ANN model had six hidden layers that were determined by trial and error and was established with selected predictors.The Sacramento streamflow data were divided into calibration period (1872 to 1957) and validation period (1958 to 1977).To evaluate the results of calibration and validation, R 2 and RMSE [82], and the Nash-Sutcliffe model efficiency coefficient [83] were computed.The RMSE describes a measure of average error in prediction and R 2 and Nash-Sutcliffe model efficiency coefficient have been widely used to assess the predictive performance of models [84].The calibration and validation results with selected predictors are shown in Figure 6.

Reconstructed ANN Model Calibration and Validation
The ANN model had six hidden layers that were determined by trial and error and was established with selected predictors.The Sacramento streamflow data were divided into calibration period (1872 to 1957) and validation period (1958 to 1977).To evaluate the results of calibration and validation, R 2 and RMSE [82], and the Nash-Sutcliffe model efficiency coefficient [83] were computed.The RMSE describes a measure of average error in prediction and R 2 and Nash-Sutcliffe model efficiency coefficient have been widely used to assess the predictive performance of models [84].The calibration and validation results with selected predictors are shown in Figure 6.The ANN model showed relatively high values of evaluation measures, with R 2 , RMSE, and Nash-Sutcliffe efficiency of 0.78, 3.66, and 0.78 in the calibration and 0.70, 3.94, and 0.72 in the validation period, respectively.Thus, the reconstructed streamflow, based on the ANN model and selected predictor, could be used as the reconstruction model.Also, it could be used for hydro-meteorological simulations.The variability of each period was 5.74 to 51.64 MAF in the calibration period and 5.13 to 32.5 MAF in the validation period.
The reconstructed streamflow using the selected predictor (nine tree rings) is shown in Figure 7 and the basic statistics are shown in Table 3, with an average of 18.9 MAF observed and 20.4 MAF of reconstructed streamflow.The ANN model showed relatively high values of evaluation measures, with R 2 , RMSE, and Nash-Sutcliffe efficiency of 0.78, 3.66, and 0.78 in the calibration and 0.70, 3.94, and 0.72 in the validation period, respectively.Thus, the reconstructed streamflow, based on the ANN model and selected predictor, could be used as the reconstruction model.Also, it could be used for hydro-meteorological simulations.The variability of each period was 5.74 to 51.64 MAF in the calibration period and 5.13 to 32.5 MAF in the validation period.
The reconstructed streamflow using the selected predictor (nine tree rings) is shown in Figure 7 and the basic statistics are shown in Table 3, with an average of 18.9 MAF observed and 20.4 MAF of reconstructed streamflow.

Bivariate Drought Analysis and Discussion
Before drought analysis based on reconstructed streamflow, the truncation level that defines the relevant streamflow level was determined to define hydrological drought from the streamflow series.The median value of annual streamflow, which was 19.4 MAF, was employed as the

Bivariate Drought Analysis and Discussion
Before drought analysis based on reconstructed streamflow, the truncation level that defines the relevant streamflow level was determined to define hydrological drought from the streamflow series.The median value of annual streamflow, which was 19.4 MAF, was employed as the truncation level to define hydrological droughts for the Sacramento River basin and results are shown in Figure 8.In total, 96 droughts occurred during the period from 1560 to 1977.Their statistical characteristics were: median drought duration of about two years, average drought severity of about 15.8 MAF, and average drought arrival time of about 2.1 years during the 15th to 20th centuries.The longest drought duration estimated was 10 years and had 75.31MAF during 1927 to 1936 (in the observation period), and the severest drought estimated was 76.17 MAF, which had an eight-year drought duration from 1582 to 1589 (in the reconstruction period).

Bivariate Drought Analysis and Discussion
Before drought analysis based on reconstructed streamflow, the truncation level that defines the relevant streamflow level was determined to define hydrological drought from the streamflow series.The median value of annual streamflow, which was 19.4 MAF, was employed as the truncation level to define hydrological droughts for the Sacramento River basin and results are shown in Figure 8.In total, 96 droughts occurred during the period from 1560 to 1977.Their statistical characteristics were: median drought duration of about two years, average drought severity of about 15.8 MAF, and average drought arrival time of about 2.1 years during the 15th to 20th centuries.The longest drought duration estimated was 10 years and had 75.31MAF during 1927 to 1936 (in the observation period), and the severest drought estimated was 76.17 MAF, which had an eight-year drought duration from 1582 to 1589 (in the reconstruction period).For bivariate drought frequency analysis with the copula function, marginal distributions of drought variables (duration and severity) were derived.The drought duration was found to have an "exponential" distribution, if treated as a continuous random variable [34].Also, the drought severity was found to have a "gamma" distribution with 95% confidential level with the PPCC (probability plot correlation coefficient; [85]) distribution goodness of fit test.Then, parameters of the Archimedean family copulas (Frank, Clayton, and Gumbel) were estimated by the method of For bivariate drought frequency analysis with the copula function, marginal distributions of drought variables (duration and severity) were derived.The drought duration was found to have an "exponential" distribution, if treated as a continuous random variable [34].Also, the drought severity was found to have a "gamma" distribution with 95% confidential level with the PPCC (probability plot correlation coefficient; [85]) distribution goodness of fit test.Then, parameters of the Archimedean family copulas (Frank, Clayton, and Gumbel) were estimated by the method of moments according to their relationship between copula parameter and Kendall's tau [86], which has been found adequate for estimating parameters for small sample sizes [78].
The minimum quadratic distance (L 2 ) between the empirical and theoretical values of the K criterion, which describes the most appropriate copula [78], was calculated for each copula.As shown in Figure 9, the Frank copula, which generally fitted well throughout (L 2 = 0.023), was selected for bivariate drought analysis for the Sacramento River.The Frank copula parameter was estimated as 8.03, and the bivariate joint probability of drought for the Sacramento River basin was described as: where, F d ptq and F s ptq are the cumulative distribution functions of drought duration and severity.
Figure 10 shows the joint CDF of the Sacramento River basin drought.
shown in Figure 9, the Frank copula, which generally fitted well throughout (L 2 = 0.023), was selected for bivariate drought analysis for the Sacramento River.The Frank copula parameter was estimated as 8.03, and the bivariate joint probability of drought for the Sacramento River basin was described as: where, ( ) and ( ) are the cumulative distribution functions of drought duration and severity.
Figure 10 shows the joint CDF of the Sacramento River basin drought.shown in Figure 9, the Frank copula, which generally fitted well throughout (L 2 = 0.023), was selected for bivariate drought analysis for the Sacramento River.The Frank copula parameter was estimated as 8.03, and the bivariate joint probability of drought for the Sacramento River basin was described as: where, ( ) and ( ) are the cumulative distribution functions of drought duration and severity.
Figure 10 shows the joint CDF of the Sacramento River basin drought.The return period is described as the average time of occurrence of events with the threshold intensity.The bivariate return period can be divided into the exceedance probabilities of both the drought duration and severity [62].The copula-based return period with the average inter-arrival of occurrences (E pLq), which was 2.1 years, can be defined as: Therefore, the duration and severity of droughts can be expressed in terms of the same return period, which can be illustrated in each "return period plot", as shown in Figure 11.
occurrences ( ( )), which was 2.1 years, can be defined as: Therefore, the duration and severity of droughts can be expressed in terms of the same return period, which can be illustrated in each "return period plot", as shown in Figure 11.The drought event and return period plot in Figure 11 shows the hydrological drought pattern of the Sacramento River basin.Overall, the drought duration and severity seemed to have a positive correlation with each other.To identify the correlation between drought duration and severity, Pearson's linear correlation coefficient, Kendall's rank correlation coefficient (τ), and Spearman's rank correlation coefficient (ρ) were found to be 0.92, 0.73, and 0.85, respectively.Also, the severity of two-year duration drought showed higher variability than did the three-year duration drought.Strong correlation between drought characteristics and higher variability of the two-year drought was the basis for the copula method, especially higher variability of severity or larger statistical irregularity than other durations for the same return period in the univariate frequency analysis [62].Hence, the copula method would be expected to be more reliable for drought analysis for the Sacramento River basin.Furthermore, California has a 22.4 MAF mean annual runoff and used 5.2 MAF annually to supply the southern area [87].So, a drought that becomes over 17.2 MAF will cause the shortage of water supply and is equivalent to approximately a 20-year return period with a two-month duration (median value of drought duration) based on the return period plot in Figure 11.Therefore, any return period that causes an actual water shortage could be the appropriate critical level of drought.Also, the return period plot in Figure 11 can be used as elementary data for water resources planning.For instance, if a decision maker or agency determined a three-year design drought for a dam or reservoir to be having a 20-year return period, then its deficiency would be The drought event and return period plot in Figure 11 shows the hydrological drought pattern of the Sacramento River basin.Overall, the drought duration and severity seemed to have a positive correlation with each other.To identify the correlation between drought duration and severity, Pearson's linear correlation coefficient, Kendall's rank correlation coefficient (τ), and Spearman's rank correlation coefficient (ρ) were found to be 0.92, 0.73, and 0.85, respectively.Also, the severity of two-year duration drought showed higher variability than did the three-year duration drought.Strong correlation between drought characteristics and higher variability of the two-year drought was the basis for the copula method, especially higher variability of severity or larger statistical irregularity than other durations for the same return period in the univariate frequency analysis [62].Hence, the copula method would be expected to be more reliable for drought analysis for the Sacramento River basin.Furthermore, California has a 22.4 MAF mean annual runoff and used 5.2 MAF annually to supply the southern area [87].So, a drought that becomes over 17.2 MAF will cause the shortage of water supply and is equivalent to approximately a 20-year return period with a two-month duration (median value of drought duration) based on the return period plot in Figure 11.Therefore, any return period that causes an actual water shortage could be the appropriate critical level of drought.Also, the return period plot in Figure 11 can be used as elementary data for water resources planning.For instance, if a decision maker or agency determined a three-year design drought for a dam or reservoir to be having a 20-year return period, then its deficiency would be about 11.1 ˆ10 9 m 3 (27.0MAF), and it could be the target storage volume for water resources planning.
Most of the droughts that occurred during the last five centuries did not have more than a 50-year return period, and just six droughts showed 100-year or longer return periods.These high return period events are one of the limitations of the study; for instance, the drought from 1927 to 1936 years, the longest and severest drought in historical data [88], with a 10-year duration and 75.31MAF severity, which equated to a 9.26 ˆ10 9 m 3 streamflow deficiency per year, had approximately a 7500-year return period.That extreme return period was due to the Frank copula and sample size, which shows some bias in high quantile (high return period) events in Figure 9, so it could have overestimated the return period [89], and the number of drought events with high quantiles is also limited.Droughts that have high quantiles or extreme return periods depend significantly on the fitted copula function.Therefore, further studies are needed for generally well-fitted copulas throughout, and also for carefully considering the use of the return period plot in water resources planning.

Conclusions
This study reconstructed the past streamflow of the Sacramento River based on the ANN and tree-ring data, and bivariate drought frequency was analyzed using the Archimedean copula.Results of this study can be summarized as follows: 1.The past streamflow for the period from 1560 to 1871 is reconstructed with the ANN model and tree-ring data, which was found to be the appropriate predictor.As shown by calibration and validation results from 1872 to 1977, the R 2 and Nash values are 0.7 or higher.It is therefore concluded that the ANN model reconstructs streamflow of the Sacramento River satisfactorily.2. Drought characteristics in the Sacramento River basin have strong correlation with each other.
The Archimedean copula is found to be appropriate for bivariate drought frequency analysis.3. It is shown that a drought with a 20-year return period or longer will cause actual water shortages in the perspective of water supply to the southern California area.Hence, it could be considered an appropriate critical level of droughts for actual water shortages.

Figure 1 .
Figure 1.Study area, tree ring sites, and streamflow sites.

Figure 1 .
Figure 1.Study area, tree ring sites, and streamflow sites.

Figure 3 .
Figure 3. ANN schematization:Input i is the input set, output is the result of network delay, and each circle represents neural network[71]; each line indicates nodes between neurons that have their own connection strength.

Figure 3 .
Figure 3. ANN schematization:Input i is the input set, output is the result of network delay, and each circle represents neural network[71]; each line indicates nodes between neurons that have their own connection strength.

Figure 5 .
Figure 5. Study area, tree ring sites, and streamflow sites; a green circle indicates that the tree ring has correlation with streamflow, while a red circle indicates correlation with time-lag (1 year) and a yellow box indicates that there is no correlation with streamflow.

Figure 5 .
Figure 5. Study area, tree ring sites, and streamflow sites; a green circle indicates that the tree ring has correlation with streamflow, while a red circle indicates correlation with time-lag (1 year) and a yellow box indicates that there is no correlation with streamflow.

Figure 7 .
Figure 7. Reconstructed streamflow using the ANN model.

Figure 7 .
Figure 7. Reconstructed streamflow using the ANN model.

Figure 7 .
Figure 7. Reconstructed streamflow using the ANN model.

Figure 10 .
Figure 10.Joint cumulative distribution function (CDF) of the Sacramento River basin drought.

Figure 10 .
Figure 10.Joint cumulative distribution function (CDF) of the Sacramento River basin drought.Figure 10.Joint cumulative distribution function (CDF) of the Sacramento River basin drought.

Figure 10 .
Figure 10.Joint cumulative distribution function (CDF) of the Sacramento River basin drought.Figure 10.Joint cumulative distribution function (CDF) of the Sacramento River basin drought.

Figure 11 .
Figure 11.Return period plot for Sacramento River.

Figure 11 .
Figure 11.Return period plot for Sacramento River.

Table 3 .
Basic statistics of observed and reconstructed streamflow.

Table 3 .
Basic statistics of observed and reconstructed streamflow.