Using short hydrological records to characterise extreme drought for water studies could lead to uncertain results. At first, this article introduces a generic model that extends the length of short hydrological records to improve the analysis of drought. Using simulated data, the extension model performance was tested and found to work well and later the model was used to extend the annual precipitation for Reba region back to 1924. The analysis of the extended precipitation shows the occurrence of 23 drought events that is about 50% more than the number of events detected from the analysis of the short historical precipitation. The analysis shows that events of 1 year duration are common while events of 3 years or more are rare. Furthermore, this article proposes a theoretical model that can be used to estimate the number of drought events that may emerge during a fixed period. The model is shown to be reliable, and the estimated drought numbers can be used in planning for the water supply from the traditional water sources. Moreover, the article suggests a few measures to improve the water supply in Jordan during dry conditions.

Drought is a natural hazard that causes a significant reduction in the water supply below the demand during a specific time (Bonaccorso et al. 2003). Each drought event is characterised by its space coverage, water deficit magnitude, duration upon which the deficit continues, and return period that tells how frequent the drought is (Salas et al. 2005). In recent decades, drought events have caused economic damage to several sectors including rain fed and irrigated agriculture (Rossi & Cancelliere 2013). It should be noted that the prediction of drought as a natural hazard remains challenging due to uncertainties (Van Huijgevoot et al. 2014). One source of the uncertainty could be attributed to the shortness of historical hydrological records used to capture the occurrence of severe drought events.

Yevjevich (1967) proposed the Theory of Runs to recognise drought events. Consequently, several research studies have investigated the magnitude, duration, and frequency of drought events based on probabilistic models (Shiau & Shen 2001; Bonaccorso et al. 2003; Salas et al. 2005; Mishra et al. 2009; Cancelliere & Salas 2010). Furthermore, long proxy data that contain natural climatic information, i.e. tree ring indices, have been used extensively in the literature to reconstruct precipitation or stream flow records back in time and later such records were used to track the past occurrence of severe drought episodes not being detected from the short historical records (Touchan et al. 1999; Gonzales & Valdes 2003; Akkemik et al. 2008; Woodhouse et al. 2009; Fang et al. 2010; Yang et al. 2014). In the literature, the duration as an important feature of the drought has been modelled using the geometric distribution (Shiau & Shen 2001; Salas et al. 2005). The probability distribution of the drought of duration l is given as:
formula
1
where pdw is the conditional probability of observing wet year given that the past year is a dry year.

While most of the water resources in the arid region of Jordan are limited and subjected to intensive stresses due to increasing water demand for the population and economic activities, the occurrence of prolonged drought events rapidly enlarges the gap between the limited supply and the water demand. In Jordan, a few drought studies were conducted employing either the drought indices (Hammouri & El-Naqa 2007) or probabilistic modelling (Tarawneh et al. 2008, 2009; Tarawneh 2011) to characterise the severity of drought. In addition, tree ring indices were used to augment short precipitation records in Jordan and later were used to capture the past occurrence of extreme meteorological droughts (Touchan et al. 1999; Tarawneh & Hadadin 2009). In summary, previous studies conducted in Jordan have characterised the severity of drought in Jordan focusing on the largest drought magnitude, the longest duration upon which the deficit extends, and the frequency of the drought. It should be noted that none of the above mentioned studies have addressed the prediction of drought events in Jordan. Besides the negative impact caused by the drought persistence on the water supply from surface and ground sources, climate change could also induce extra pressure on water resources by limiting the feed water. In the semi-arid region of Jordan, there are few studies that indicated a change in the climate. Smadi & Zghoul (2006) indicated the existence of a downward trend in the annual rainfall amounts for many regions in Jordan. Also, a warming trend in the summer season at a rate of 0.038 °C/year has been detected which would increase the demand on water for irrigating summer crops. In a recent study, Rahman et al. (2015) also indicated the decrease in the annual rainfall amounts in Jordan at an average rate of 1.2 mm/yr. Al-Houri (2014) analysed the rainfall in Jordan and indicated a decreasing trend in the duration of the wet season associated with a decrease in the number of rainy days. In summary, due to the nature of the arid region of Jordan, a slight decrease in rainfall due to climate change will ultimately create extra stress on the surface and groundwater storages.

In Jordan, most of the precipitation records available for water studies, including drought analysis, are short to witness the occurrence of probable extreme drought, i.e. the longest drought. While precipitation records exist for most regions in Jordan, the earliest started in 1924 and the majority of records started after 1950. In the literature, the simple linear regression has been used to extend the length of short hydrologic records benefiting longer records, for instance tree ring indices, when strong cross correlation exists. The simple linear regression model is written as:
formula
2
where is the extension of the short record, xt is the long record, and are the mean values of the short record and long record over the calibration period, b is the model parameter and ne is the extension period. The quality of extensions produced using the linear regression (Equation (2)) relies only on the strength of correlation between the short and the long record, however if the short record exhibits a temporal (serial) correlation, then extensions produced using Equation (2) may not acquire the temporal correlation of the short record, especially when the long record is time independent. Such a drawback could appear in hydrologic studies when river flows that usually show significant serial correlation are reconstructed from a serially independent variable like tree ring data.

One objective of this article is to present a generic theoretical model that extends the length of short hydrological records while accounting for the temporal correlation of the short variable. It is believed that extending the length of short hydrological records will increase the opportunity to capture more extreme drought events, ultimately, water management studies will be improved. Other objectives are to introduce a simple theoretical model that enables the computation of the number of drought events of specific duration, and to analyze the occurrence of drought in the Raba region, Jordan.

All agricultural activities in the region of Raba depend on the rainfall during the wet season that extends over the months October to May. Therefore, the rainy season precipitation for Raba will be used to characterize drought. Monthly precipitation data for Raba were obtained from the Jordan Meteorological Department for the period 1953–2012. The rainy season precipitation was obtained after accumulating the monthly precipitation over the wet season. The 60 years of precipitation in Raba could be short to capture extreme drought events. Therefore, precipitation data for Amman airport gauging station (1924–2012) will be used to extend the precipitation in Raba back to 1924 resulting in 89 years of precipitation, i.e. around 130% larger. Figure 1 shows the geographic location of Amman airport and Raba gauging stations. The simple linear regression model (Equation (2)) can be used to extend the historical precipitation for Raba, however, the precipitation will be extended using the following record extension model. The proposed model will account for the serial autocorrelation of the short hydrologic variable if it exists, i.e. the streamflows that are usually time correlated. The modified record extension model has the following form:
formula
3
where is the one step lagged observation of the targeted variable, a and b are the model parameters. Referring to the simple traditional linear regression (Equation (2)), the proposed model (Equation (3)) with the additional term, , considers the short variable lagged information contained in the past value (time t−1) to predict the current value (time t). Using the method of moments, the model parameters a and b are estimated as (Appendix A, available with the online version of this paper):
formula
4
and
formula
5
where C(XY) is the lag-0 covariance between yt and xt, C(XX) is the lag-0 covariance of xt, C(XY) is the lag-1 covariance between xt and the one time step lagged yt, C(YY) and C(YY) are the lag-0 and lag-1 covariance of yt, respectively.
Figure 1

Location of Raba and Amman airport gauging stations.

Figure 1

Location of Raba and Amman airport gauging stations.

Close modal
In order to improve the estimation of the multiyear drought characteristics for water resources management studies, the extended precipitation record will be used to analyse the drought occurrence in the Raba region. Furthermore, knowledge about the expected number of drought events of particular duration that could arise during the water resources operation time could be useful. Given a hydrologic record of sufficient length, if the number of all drought events detected is N and the number of drought events of specific length is Nl, then the probability that a drought of length l years occurs is estimated empirically as follows:
formula
6
Equating Equations (1) and (6), simplifying and taking expectation on both sides, then the expected number of drought events of the desired length l is:
formula
7
If the hydrologic variable has fixed length (M years), and let pd be the unconditional probability of observing a dry year, then the total number of dry years is pd × M, which is the sum of lengths for all drought events, i.e. N × L = pd × M. Taking the expectation and rearranging, the expected number of all drought events, E[N] is:
formula
8
where E[L] is the expected length of drought.
Substituting Equations (8) in (7), and noting that the expected run length E[L] = 1/pdw (geometric distribution of the run length L), then Equation (7) becomes:
formula
9
Although using simple expressions, Equation (9) enables the computation of the expected number of the l years drought events that could emerge during a fixed time period M. Such a fixed period could be any planning time for a water resource system. Furthermore, Equation (9) is useful to develop an expression that estimates the recurrence time (R) of drought events of specific type. Based on the renewal theory of events, Loaiciga & Leipnik (1996) expressed the recurrence time (R) of hydrologic events emerging during a fixed time t as follows:
formula
10
where Ne is the number of events arising during the period t. If the number of drought events Ne is estimated by its expected number (Equation (9)) and noting that t = M, then the recurrence time (R) is:
formula
11

Equation (11) is the same expression developed by Salas et al. (2005) to evaluate the recurrence time of hydrologic events.

Extension model performance

Although the proposed record extension model has been used here to extend the time independent annual rainy season precipitation for Raba given the annual precipitation for Amman airport, the target is to introduce a general model that can be used under different correlation conditions. For example, it can be used to extend the time dependent streamflows from precipitation or tree ring data (time independent) while preserving the serial correlation of the streamflows that cannot be achieved using the simple linear regression (Equation (2)). In order to check the suitability of the proposed model, the performance of the proposed extension model (Equation (3)) has been tested versus the widely used simple linear regression model using simulation experiments. The performance of both models will be tested at different cross correlation values assuming the variable to be predicted is time dependent. At cross correlation between yt and xt of 0.5, a normally distributed pair (yt and xt) was generated to a length of 150 data points. The mean and standard deviation of both, yt and xt, were selected arbitrarily as 100 and 30, respectively. The serial correlation coefficient of yt was selected as 0.3, i.e. shows significant temporal dependence like stream flows, while the serial correlation coefficient of the predictor variable xt was kept at low value (time independent) similar to most of tree-ring indices used as predictors to extend short precipitation or stream flow records. Out of the yt 150 data points, the last 50 (t = 101, 102, …, 150) were kept aside for model evaluation purposes. Therefore the variable yt was considered the short variable of 100 data points length (n = 100). The proposed and the simple linear regression models were used to estimate the extensions for t = 101, 102, …, 150. The extended part (t = 101, 102, …, 150) was combined with the short part yt (t = 1, 2, …,100) forming the combined record and the serial correlation was computed and checked versus the serial correlation of the short (yt). The performance statistic was computed. Since the proposed record extension model targets the production of extensions that will preserve the serial correlation of the short record, then the model performs well if the statistic is close to 1. The same procedure was repeated 1,000 times and finally the expected value of the performance statistic was computed. The same simulation experiments and analysis procedure were repeated at cross correlation values 0.65, 0.8 and 0.9. Figure 2 shows the expected value of the performance statistic computed using extensions predicted by the proposed and the traditional linear regression models at different cross correlation values.
Figure 2

The expected value of the performance statistic computed from extensions generated using both models versus the cross correlation.

Figure 2

The expected value of the performance statistic computed from extensions generated using both models versus the cross correlation.

Close modal

The expected value of the performance statistic computed from extensions achieved using the traditional regression model (Figure 2) is always less than 1, i.e. the traditional regression model generated extensions that failed to preserve the serial correlation of the short record. In that case, extended values were affected somehow by the independence feature of the long record (xt). The effect of the independence feature of the long record (xt) on the generated values increases as the cross correlation increases, i.e. as the power of the information transfer increases. In conclusion, the performance statistic shown in Figure 2 indicates that extensions produced using the traditional regression model did not maintain the serial correlation of the original short record (yt). As the cross correlation between the short (predicted) and the long record (predictor) increases, the characteristics of the predicted short variable are highly affected by the features of the long recorded variable rather than preserving the features of the original variable. Regarding the performance of the proposed record extension model, Figure 2 shows a consistent performance statistic that is very close to 1 regardless of the variation in the cross correlation values. Although the quality of the information transferred is highly affected by the cross correlation between the xt and yt, the additional term in Equation (3) maintained the serial correlation of the short variable (yt). In general, the proposed extension model succeeded to produce extensions that maintain the original short record serial correlation. Other statistics to measure the performance of both models to generate extensions that preserve features like the mean, variance, cross correlation and sum of squared error were also computed (results are not shown). In general, the analysis of results indicates the ability of both models to generate extensions that reproduce the mean, variance, and sum of squared error of the short variable.

Analysis of multiyear drought events in Raba

Due to the lack of streamflow data in Raba or nearby regions, drought will be analysed using precipitation. The annual rainy season precipitation for Raba covering the period 1953–2012 has been extended back to the year 1924 using Equation (3), benefiting from the annual precipitation for Amman airport that began in 1924. The cross correlation between the precipitation of both sites over the concurrent records period is 0.77 which may indicate a similar climatic pattern. The record extension model parameters, a and b, were estimated over the model calibration period (1953–2012) using Equations (4) and (5). The parameter a= 0.132 and b = 0.973. The combined precipitation series for Raba (1924–2012) consists of the historical part (1953–2012) and the extended part (1924–1952). Figure 3 shows the combined annual precipitation series for Raba region truncated at the level of the long-term mean.
Figure 3

The combined annual precipitation series for Raba (1924–2012). The dashed line represents the truncation level.

Figure 3

The combined annual precipitation series for Raba (1924–2012). The dashed line represents the truncation level.

Close modal

The combined annual precipitation series shown in Figure 3 was truncated at the level of the long-term mean to identify drought events. Similarly, the historical precipitation for Raba was also truncated at the level of the long-term mean for drought characterisation. In the arid region of Jordan, the selection of the long-term mean as a truncation level could be an appropriate choice due to water scarcity, i.e. any small deviation below the average precipitation reduces the water amounts supplied to the rain fed agricultural activities. The analysis of the truncated historical precipitation (1953–2012) indicated the occurrence of 15 drought events, among them the 5-year drought and the 6-year drought that occurred during the periods 1975–1980 and 1958–1963 respectively (Figure 3). On the other hand, the analysis of the truncated combined annual precipitation series (1924–2012) indicated the occurrence of 23 drought events, which is around 1.5 times larger than the number of drought events detected from the analysis of the historical precipitation. Among the 23 drought events, there are two severe drought events each lasting 5 years during the periods 1975–1979 and 1930–1934, and one severe 6-year drought event during the period 1958–1963 (Figure 3). In order to investigate the distribution of the drought events in Raba region, the number of events and their occurrence probability were obtained from the analysis of both truncated precipitation series. Table 1 shows the number of detected drought events and their occurrence probability against the drought length. In general, drought events of short duration like the 1-year and probably the 2-year events, are the common (frequent) types of drought in Raba region with occurrence probability of 0.57 and 0.17 for the 1- and 2-year events respectively. Events of the duration of 3 years or more, although they do occur, are a less common type of drought with low occurrence probability, i.e. P[L ≥3 years] = 0.26. The distribution of drought events of 2 years or more detected from the analysis of the historical precipitation appears to be inconsistent (Table 1) compared to events detected from the combined precipitation (relatively long record), which attributed to the limited number of historical drought events observed.

Table 1

Number and occurrence probability of drought events detected from combined and historical precipitation record for Raba region

 Number of events (probability)
Drought length (year)Combined precipitation 1924–2012Historical precipitation 1953–2012
13 (0.57) 8 (0.53) 
4 (0.17) 2 (0.13) 
2 (0.09) 2 (0.13) 
1 (0.04) 1 (0.07) 
2 (0.09) 1 (0.07) 
1 (0.04) 1 (0.07) 
 Number of events (probability)
Drought length (year)Combined precipitation 1924–2012Historical precipitation 1953–2012
13 (0.57) 8 (0.53) 
4 (0.17) 2 (0.13) 
2 (0.09) 2 (0.13) 
1 (0.04) 1 (0.07) 
2 (0.09) 1 (0.07) 
1 (0.04) 1 (0.07) 

To validate how well Equation (9) enables the computation of the number of drought events of a desired length that could occur during a specified time, it has been used to compute the number of observed drought events that had already occurred during the precipitation period 1953–2012 (historical precipitation, M = 60 years). At first, the unconditional and conditional probabilities pd and pdw were estimated from the truncated historical precipitation. The probabilities pd and pdw are 0.528 and 0.472, respectively. Using Equation (9) with M = 60, the expected number of drought events were computed against the drought length. Figure 4 shows the number of observed (actual) drought events and their expected number computed using Equation (9) versus the length of drought during the historical precipitation period for Raba region. Figure 4 shows the reliability of Equation (9) to closely estimate the number of the observed drought events of specific length during a fixed period. The total number of estimated drought events is 14.6 events, which is very close to the total number of observed events (15 events). For the most common type of events, the 1-year drought, the expected number is seven events compared to eight observed events (Figure 4). The discrepancy between the expected and the observed 2-year drought events shown in Figure 4 is attributed to the limited number (two events) of drought detected from the analysis of the short historical precipitation.
Figure 4

Observed and estimated number of drought events versus their length.

Figure 4

Observed and estimated number of drought events versus their length.

Close modal

The merit of using Equation (9) to estimate the number of drought events during a specified time arises from its simplicity and the usefulness of its product to water resource management studies. It shows flexibility to accommodate any planning period, for example M could be 5, 10, 15 years. Therefore, Equation (9) can be used to construct charts that display the expected number of drought events that could occur in any region given simple statistics like pd and pdw.

In order to estimate the recurrence time of the common drought events in Raba region, Equation (11) has been used. Given the probabilities pd and pdw, the recurrence times for events with a duration of 1 and 2 years are 8.5 and 16.1 years, respectively, whereas the less common events (3 years or more) have recurrence times of more than 30 years.

Measures to reduce the impact of drought and climate change on the water supply

Due to shortage in the water supply from the limited natural water resources, the Ministry of Water and Irrigation in Jordan has constructed several macro-scale rainwater harvesting projects during the past few decades. Such projects were designed for different tasks including groundwater basin recharge and fresh water supply for municipal, agricultural, and industrial sectors. Table 2 shows the storage capacity and the purpose of rainwater harvesting for selected projects in Jordan.

Table 2

Storage capacity in million cubic meters (MCM) and purpose of selected rainwater harvesting projects in Jordan

Project nameCapacity (MCM)PurposeNotes
Mujib dam reservoir 32 Supply fresh water for agriculture and industry Collecting the winter flood and base flow of the Mujib valley 
Wala dam reservoir Supply fresh water for agriculture and recharging groundwater Collecting winter flood 
Qutrana dam reservoir 4.2 Supply fresh water for agriculture Collecting winter flood 
Project nameCapacity (MCM)PurposeNotes
Mujib dam reservoir 32 Supply fresh water for agriculture and industry Collecting the winter flood and base flow of the Mujib valley 
Wala dam reservoir Supply fresh water for agriculture and recharging groundwater Collecting winter flood 
Qutrana dam reservoir 4.2 Supply fresh water for agriculture Collecting winter flood 

Regardless of the intensive efforts made by the administration to harvest rainwater during the short rainy season, such resources are always subjected to stress due to the increasing water demand for drinking, agricultural, and industrial purposes. Furthermore, besides the arid climate of Jordan that usually brings frequent dry periods, the decreased trend in rainfall detected lately due to the potential change in the climate increases the pressure on the limited water resources. For example, considering an overall reduction in the total rainfall of 0.41 mm/yr (Rahman et al. 2015) detected after 1995, Figure 5 shows the potential reduction in the total annual rainfall amount over Jordan after 1995. Figure 5 shows that the feed water has decreased from the long-term average value of 7,300 MCM in 1995 to nearly 6,550 MCM in 2013, which is approximately a 10% reduction due to the change in the climate. In a country where water is scarce, such a reduction is remarkable.
Figure 5

The expected rainfall (bars) after 1995 due to climate change versus the long-term average of the total rainfall in Jordan (solid line).

Figure 5

The expected rainfall (bars) after 1995 due to climate change versus the long-term average of the total rainfall in Jordan (solid line).

Close modal

The following measures, if implemented, could contribute to reducing the devastating impact of drought persistence, the effect of the change in climate, and the increasing demand on the water sector in Jordan. Therefore the following is recommended:

  1. Embed drought persistence in water release policies from surface water resources. In a previous study (Tarawneh 2011), it was concluded that unplanned water releases combined with less attention to drought persistence possibly contributed to the remarkable decline in the water storage in the Mujib dam reservoir during the years 2008 and 2009, upon which the water supply reached a critical level.

  2. Develop new water sources. It could be time in Jordan to develop new water sources to compensate the lowering trend in the annual rainfall amount as a result of climate change and the increasing demand on water. One promising potential source that could provide unlimited amounts of fresh water is sea water desalination. Current desalination practices in Jordan are running at a modest scale in the Jordan Valley after desalting the brackish in small scale private units. The Red Sea-Dead Sea canal project has been implemented, which targets the restoration of the Dead Sea water level and hydropower energy production through conveying seawater to the Dead Sea. After desalination, it is expected that the project will provide Jordan with around 500 MCM of fresh water per year.

  3. Adapt economic incentive plan to encourage water savings and the efficient use of water in irrigation where huge amounts of fresh water are usually consumed. Introducing new irrigation techniques and properly rehabilitating the existing aged water distribution systems are essential measures for saving water. The current water loss percentage reaches 42%, which is relatively high. Further reduction in the loss rate is needed which can be achieved through incentives, probably lowering the price of water when less water is used.

  4. Adapt proper crop changing policy during the dry periods. While around 60% of the water released from surface and groundwater sources in Jordan is used mainly to irrigate crops, irrigating less water demanding crops would possibly lessen the water shortage problem during dry periods. Drought studies would help in advising farmers to adapt less water demanding crops that could survive under dry conditions.

Long dry periods threaten the water supply of traditional water resources in Jordan. To reduce the uncertainty in drought characterisation due to using short hydrological records, this article introduces a generic record extension model that augments short hydrologic records while preserving the features of the hydrologic variable. Using simulation experiments, the model performance was tested at different cross correlation values and extensions produced were found to preserve the characteristics of the original variable. Therefore such a model can be used to extend the length of short hydrologic records targeting the improvement of drought characterisation. Furthermore, the article presents a simple theoretical model that facilitates the estimation of the number of drought events that could emerge during a fixed time. The proposed model was tested to foretell the number of drought events that occurred in Raba region. Such a model can be used to estimate the number of drought events needed for water resources management studies.

After extending the length of the Raba precipitation record, the analysis of the extended precipitation indicates the occurrence of 23 drought events, the longest is a 6-year event. However, events of short duration, one and 2 years, are the common type of events in Raba with recurrence time of 8.5 and 16 years, respectively. While part of the water shortage in Jordan is attributed to a drought persistence problem, it is recommended that the water sector institutions in Jordan consider drought persistence when planning for water releases. Moreover, the reduction in precipitation due to climate change increases the gap between water supply and demand. Therefore, it is recommended to develop new water sources in Jordan through increasing desalination practices to provide fresh water to compensate for the increasing water demand and the lowering trend in the rainfall amounts due to climate change. Furthermore, more efforts are needed to save water through rehabilitating existing water distribution systems, adapting economic incentives to encourage the efficient use of water, and adapting proper crop changing policy during the dry periods.

Akkemik
U.
D'Arrigo
R.
Cherubini
P.
Kose
N.
Jacoby
G.
2008
Tree-ring reconstructions of precipitation and streamflow for North-western Turkey
.
Int. J. Climatol.
28
,
173
183
.
Al-Houri
Z.
2014
Detecting variability and trends in daily rainfall characteristics in Amman-Zarqa basin, Jordan
.
Int. J. Appl. Sci. Technol.
4
(
6
),
11
23
.
Bonaccorso
B.
Bordi
I.
Cancelliere
A.
Rossi
G.
Sutera
A.
2003
Spatial variability of drought: an analysis of the SPI in Sicily
.
Water Resour. Manage.
17
,
273
296
.
Fang
K.
Davi
N.
Gou
X.
Chen
F.
Cook
E.
Li
J.
D'Arrigo
R.
2010
Spatial drought reconstructions for central high Asia based on tree rings
.
Clim. Dyn.
35
,
941
951
.
Gonzales
J.
Valdes
J.
2003
Bivariate drought analysis using tree ring reconstruction
.
J. Hydrol. Eng.
8
(
4
),
247
257
.
Hammouri
N.
El-Naqa
A.
2007
Drought assessment using GIS and remote sensing in Amman-Zarqa basin, Jordan
.
Jordan J. Civil Eng.
1
(
2
),
142
152
.
Loaiciga
H.
Leipnik
R.
1996
Stochastic renewal model of low-flow streamflow sequence
.
Stoch. Hydrol. Hydraul.
10
(
1
),
65
85
.
Mishra
A.
Singh
V.
Desai
V.
2009
Drought characterization: a probabilistic approach
.
Stoch. Environ. Res. Risk Assess.
23
,
41
55
.
Rahman
K.
Gorelick
S.
Dennedy-Frank
P.
Yoon
J.
Rajaratnam
B.
2015
Declining rainfall and regional variability changes in Jordan
.
Water Resour. Res.
51
(
5
),
3828
3835
.
Rossi
G.
Cancelliere
A.
2013
Managing drought risk in water supply systems in Europe: a review
.
Int. J. Water Resour. Dev.
29
(
2
),
272
289
.
Salas
J.
Fu
C.
Cancelliere
A.
Dustin
D.
Bode
D.
Pineda
A.
Vincent
E.
2005
Characterizing the severity and risk of drought in the Poudre River, Colorado
.
J. Water Resour. Plann. Manage.
131
(
5
),
383
393
.
Shiau
J.
Shen
H.
2001
Recurrence analysis of hydrologic droughts of different severity
.
J. Water Resour. Plann. Manage.
127
,
30
40
.
Tarawneh
Z.
Hadadin
N.
2009
Reconstruction of the rainy season precipitation in central Jordan
.
Hydrol. Sci. J.
54
(
1
),
189
198
.
Tarawneh
Z.
Elgaali
E.
Moshrik
H.
2008
Severity of droughts in arid regions
.
Jordan J. Civil Eng.
2
,
228
237
.
Touchan
R.
Meko
D.
Hughes
M.
1999
A 396-year reconstruction of precipitation in Southern Jordan
.
J. Am. Water Resour. Assoc.
35
(
1
),
49
59
.
Yang
B.
Qin
C.
Wang
J.
He
M.
Melvin
T.
Osborn
T.
Briffa
K.
2014
A 3,500-year tree-ring record of annual precipitation on the Northeastern Tibetan Plateau
.
Proc. Natl. Acad. Sci. USA
111
,
2903
2908
.
Yevjevich
V.
1967
An Objective Approach to Definitions and Investigations of Continental Hydrologic Droughts
.
Hydrology paper 23
,
Colorado State University
,
Fort Collins, Colorado
,
USA
.

Supplementary data