Drought indices based on precipitation are commonly used to identify and characterize droughts. Due to the general complexity of droughts, the comparison of index-identified events with droughts at different levels of the complete system, including soil humidity or river discharges, relies typically on model simulations of the latter, entailing potentially significant uncertainties.
The present study explores the potential of using precipitation-based indices to reproduce observed droughts in the lower part of the Jinsha River basin (JRB), proposing an innovative approach for a catchment-wide drought detection and characterization. Two indicators, namely the Overall Drought Extension (ODE) and the Overall Drought Indicator (ODI), have been defined. These indicators aim at identifying and characterizing drought events on the basin scale, using results from four meteorological drought indices (standardized precipitation index, SPI; rainfall anomaly index, RAI; percent of normal precipitation, PN; deciles, DEC) calculated at different locations of the basin and for different timescales. Collected historical information on drought events is used to contrast results obtained with the indicators.
This method has been successfully applied to the lower Jinsha River basin in China, a region prone to frequent and severe droughts. Historical drought events that occurred from 1960 to 2014 have been compiled and cataloged from different sources, in a challenging process. The analysis of the indicators shows a good agreement with the recorded historical drought events on the basin scale. It has been found that the timescale that best reproduces observed events across all the indices is the 6-month timescale.
Drought is a natural phenomenon that results from persistent deficiency of precipitations over an extended period of time compared with some long-term average condition (e.g., precipitation), resulting in a water shortage for some activity, group, or environmental sector (Landsberg, 1982). It generally affects larger areas than other hazards and more people than any other natural catastrophe (Keyantash and Dracup, 2002; Wilhite, 2000).
In China, droughts represent the most severe natural threat for socioeconomic development and ecosystems (Mei and Yang, 2014). Drought events occur in the Jinsha River basin (JRB) and surrounding regions with high frequency. They affect a wide range of areas and cause huge losses to the agriculture sector (He et al., 2013). The clustering of severe and sustained droughts in southwest China during the last decade has resulted in tremendous losses, including crop failure, lack of drinking water, ecosystem degradation, health problems and even deaths (Wang et al., 2015).
To reduce and anticipate such drought impacts, a comprehensive characterization of the phenomenon is essential. Effective and accurate analysis of hydrometeorological data is a key input. Drought indices are useful for tracking droughts and providing a quantitative assessment of the severity, location, timing and duration of such events (World Meteorological Organization and Global Water Partnership, 2016), but also for real-time monitoring (Niemeyer, 2008), risk analysis (Hayes et al., 2004) and drought early warning (Kogan, 2000).
Some organizations and agencies already rely on the use of indices in their decision-making processes, thus enhancing proactive drought management policies (Wilhite, 2000). An example is the U.S. Drought Monitor (USDM, 2017), an index-based drought map that policymakers use in discussions of drought and in allocating drought relief. Other platforms such as the European Drought Observatory (Joint Research Centre, 2017); China's Department of Climate Change, National Development and Reform Commission (Department of Climate Change, National Development and Reform Commission, 2017); or the experimental African Flood and Drought Monitor (Land Surface Hydrology Group – Princeton University, 2017) also use this approach for the assessment, diagnosing and forecasting of droughts.
The choice of index should be based on the type of drought (meteorological, agricultural, hydrological or socioeconomical), the climate regime and the regions affected, as well as the available data. It was found that measured meteorological data were limited in the study region and that precipitation was the single most reliable type of exploitable information. The present study thus focuses on the use of meteorological indices based only on precipitation data. The main advantages are their ease of use, the limited data requirements and the capacity for early detection of drought events, while extensive literature and calculation tools are widely accessible (World Meteorological Organization and Global Water Partnership, 2016).
The use of integrated indices such as the Palmer drought severity index (PDSI; Palmer, 1965) or standardized precipitation–evapotranspiration index (SPEI; Vicente-Serrano et al., 2010), relying also on potential evapotranspiration (PET), could improve the scope and quality of this study. However, no reliable information necessary to calculate PET was accessible for the study region. Although approximations may be applied to estimate this variable, for example, by only considering temperature data, some studies (Jeevananda Reddy, 1995; Shaw and Riha, 2011; Staage et al., 2014) showed a high sensitivity of the PET to the chosen equation. A deeper analysis that helps in selecting and applying such methods should be performed prior to the use of these indices. Therefore, it has been decided to base this study on the standardized precipitation index (SPI, McKee et al., 1993, 1995), the rainfall anomaly index (RAI, Van Rooy, 1965), the percent of normal precipitation (PN, Barua et al., 2011) and the deciles index (DEC, Gibbs and Maher, 1967).
To fill the lack of specific drought-related information, most studies assess the performance of drought indices against results from hydrological soil water models (Halwatura et al., 2016; Hao and AghaKouchak, 2013; Trambauer et al., 2014; Vasiliades et al., 2011; Wanders et al., 2010). However, the performance of these types of studies depends on the accuracy of the models. Their limitations and uncertainties represent an important drawback and should be addressed (Mishra and Singh, 2011). An alternative that often requires more time-consuming work is the compilation of historical records of drought events from different sources. Consequently, their duration, the water scarcity levels, and the drought impacts on population and agriculture can be estimated and then integrated into the analysis. This enables one to identify other types of droughts such as socioeconomical droughts that are hard to assess with hydrological models.
Regarding their spatial resolution, the available drought indices may be based on local measurements (Zhou et al., 2012) and index calculations are usually applied to stations or cells of gridded precipitation datasets; overall spatial patterns on catchment or sub-catchment scales are thus hardly captured. As stated above, droughts affect large areas whose limits are often vaguely demarcated. In addition, water resources are part of a more complex interrelated network which links the source to the point of consumption, where isolated rainfall deficiencies do not necessarily imply a shortage of water availability or even a drought event. Some work (Bhalme and Mooley, 1980; Fleig et al., 2011; Mitchell et al., 1979) suggests the use of drought area indices for the study of droughts that considers areal coverage. The use of overall indicators capable of capturing in a single value the effect of the rainfall deficiency at a regional level is thus convenient and will be applied in this study based on the above-mentioned work.
The objective of this study is to capitalize on the collection of drought events that the authors have registered in the lower part of the JRB since 1960 to evaluate and calibrate two indicators capable of identifying drought occurrence and characterizing their intensity on the catchment scale. These indicators are based on commonly used meteorological drought indices for particular timescales.
The JRB is a sensitive zone in terms of water resources, food security,
ecosystem management and human well-being where glacier and climatic
variability greatly influence the water regimes and availability.
Originating from the southern glacier at Jianggendiru peak, the highest
point of the Geladaindong snowy mountain in the middle of the Tanggula
Mountains, the JRB, constitutes the upper part of the Yangtze River basin. It
is located between
The lower part of JRB is a hot–dry valley region characterized by a
southwest monsoon climate. The hydrologic regime is characterized by a
pronounced seasonal cycle with an annual average precipitation
of 600–800 mm yr
Droughts occurring in the lower JRB and surrounding areas affect a wide
range of areas, causing huge losses in agriculture (Wu et al.,
2011): more than 4 million people and 3 million livestock face a drinking
water shortage, and more than 1 million hm
Figure 1 shows the division of the JRB in three parts (upper, middle and lower) and the locations of the meteorological stations used. This study focuses on the analysis of drought events in the lower JRB. The precipitation data needed in this study have been obtained from the China Meteorological Data Service Center (CMA) and downloaded from its data sharing service system (CMDC, 2017). A preliminary quality check and correction of datasets (including data gap filling) is performed by CMA before uploading them to the system. The monthly precipitation data of 29 meteorological stations within or around JRB, recorded from 1960 to 2014, have been collected and processed. More than 50 years of continuous data are thus available, except for the Batang and Yuanmou stations where only 46 years are available. The spatial distribution of the stations is assumed to be adequate for the purposes of the study: the stations are distributed relatively evenly both in the zonal and meridional directions, with no zones having a significantly denser presence of stations that could overestimate their importance.
In order to obtain a good basis for the evaluation of drought indices' performance, historical drought events have been collected since 1960. The information required for the identification and characterization of major droughts in the lower JRB has been compiled from different sources, including scientific literature, inventories (e.g., international disaster database, Chinese inventories), governmental reports and yearbooks, and newspaper and internet articles.
Detailed information is available for the major drought events over the past 20 years. Before 1980, much less information about droughts in the lower JRB is available. Moreover, detailed information prior to 1960 could not be found.
Compiling and harmonizing the information from these different data sources
was a challenge. The drought event descriptions in the scientific literature
often give an overview about the entire event in a descriptive way without
detailed information about the affected area and damage. The literature often
provides information about the meteorological conditions and the duration of the
event. Available information from government reports and databases, in
contrast, generally contain information with a high level of detail for a specific county
(e.g., affected areas in km
Information at very different levels of detail and various contents were collected. As a first step, all available information was registered in a database. For this purpose, a web-based event registration platform and database (GEOTEST AG, 2017) has been developed to provide a standardized analytical framework with a quantitative description of the drought characteristics. In particular, drought events since 2000 are all mentioned in a different source, which significantly enhances the reliability of their existence and related information. Most of the drought events before 2000 are documented in detail in He (2010), and some of these major dry periods are also mentioned in scientific literature such as He et al. (2016) and Wang et al. (2015). Even if the amount of data and level of detail is lower for these older events, their occurrence and temporal positioning can be assumed to be reliable.
In a second step, the level of detail was harmonized for the most relevant information and summarized in Table 1. This catalog of the most relevant drought characteristics focuses on the affected area, the start date of the drought events, their duration, their spatial and temporal distribution, their severity, and the impacts on the population and agriculture, including damage and financial losses. The time of occurrence and duration is given in seasonal units, the affected area is described on a county scale, and the indicators and impacts are summarized in a descriptive way, as better accuracy was not feasible for all events.
Catalog of historical drought collected for the lower JRB (DJF: December–January–February; MAM: March–April–May; JJA: June–July–August; SON: September–October–November).
In total, 13 major drought episodes have been registered from 1960 until 2014. However, this dataset is probably not complete, as non-documented events likely have occurred. A clustering of severe and sustained droughts in the JRB has been observed from 2009 to 2014. Another period with high drought activity and severity can be detected between 1980 and 1990. Although the droughts identified from 2009 until 2014 were extremely serious, this was not the worst period in the long term because the drought episodes that occurred around 1940 were of similar intensity and duration (Wang and Chen, 2012). The registered drought events are often strongly correlated with low precipitation, but the analyses also reveal that the registered droughts often occur during periods with temperatures above average.
Four different commonly used meteorological drought indices have been applied in this study: the standardized precipitation index, the rainfall anomaly index, the percent of normal precipitation and the deciles index. Their definition basically rests upon the comparison of precipitation values with the normal value (the definition of “normality” may vary from one index to another), resulting in a single number. This allows characterizing drought conditions and thus facilitating the index's interpretation and use in strategic planning and operational applications (Tigkas et al., 2013).
This comparison must be month or season specific. For instance, for the index calculation of January 2000, the precipitation of this month should be compared to the normal precipitation extracted taking into account only the Januaries from a reference period. The same applies when calculating the index for the time window January–February–March 2000: the precipitation for these 3 months will be compared to the sum of precipitation of all the groups of January-February-March registered in the reference period.
The widely used standardized precipitation index was formulated by McKee et al. (1993, 1995) to quantify the precipitation deficit from long-term recording and for multiple timescales.
The long-term record of precipitation values is fitted to a probability distribution which is then transformed into a standard normal distribution, of which mean and variance are 0 and 1, respectively (Edwards and McKee, 1997). The datasets are most commonly adjusted to the gamma function (McKee et al., 1993; Sönmez et al., 2005; Tsakiris et al., 2007) although some studies show better adjustments to other functions (Akbari et al., 2015).
A classification of drought conditions based on the SPI values was established by McKee et al. (1993) to define drought intensities and is presented in Table 2. Positive SPI values indicate greater than normal precipitation, and negative values indicate lower than normal precipitation.
Classification of drought conditions according to the SPI values.
As mentioned earlier, the SPI was designed to quantify the precipitation deficit
for multiple timescales or moving time windows (World
Meteorological Organization, 2012). These timescales reflect the drought
impacts on different water resources which are needed by decision makers:
The 3-month SPI reflects short- and medium-term moisture conditions and
provides a seasonal estimation of precipitation. The 6-month SPI indicates seasonal to medium-term trends in precipitation and
may be very effective in showing the precipitation anomaly over distinct
seasons. Information from a 6-month SPI may also be associated with
anomalous streamflow and reservoir levels, depending on the region and time
of year. The 12-month up to 24-month SPI reflects long-term precipitation patterns and
is usually tied to streamflow, reservoir levels and even groundwater levels
on longer timescales.
The rainfall anomaly index was developed by Van Rooy (1965). The RAI indices are computed by comparing the average precipitation over a given time window with the mean of the 10 highest (for positive anomalies) and the 10 lowest (for negative anomalies) precipitation records. Despite its simplicity, this index requires a series of complete data to be calculated.
The RAI values are classified (Van Rooy, 1965) as shown in Table 3. Olukayode Oladipo (1985) found that differences between the RAI and the more complicated indices of the Palmer drought index (Palmer, 1965) and Bhalme and Mooley drought index (Bhalme and Mooley, 1980) were negligible.
Classification of the period according to the values of the RAI.
The percent of normal precipitation is one of the simplest measurements of precipitation value for a location. It is calculated by dividing precipitation during a given time window by normal precipitation of that same time window over the reference period (typically considered to be a 30-year average). For PN values over 100 %, the precipitation is higher than the average precipitation (and vice versa): the higher the PN value, the wetter the considered month is.
The main advantage of this index is its simplicity and transparency, which makes it practical to communicate drought levels to the public (Keyantash and Dracup, 2002). Analyses using PN are very effective when used for a single region and/or a specific season.
Even if no threshold ranges have been widely established in the technical literature for the PN, some studies (Barua et al., 2011; Morid et al., 2006) propose a classification similar to the SPI. For this study, the classification proposed by Barua et al. (2011) has been adopted (Table 4).
Classification of drought conditions according to the PN values.
Another drought-monitoring technique consists in dividing the monthly precipitation data into deciles. This method, developed by Gibbs and Maher (1967), was selected as the meteorological measurement of drought for the Australian Drought Watch System (Lee, 1979; Sivakumar et al., 2010) because it is relatively simple to calculate and requires less data and fewer assumptions than the Palmer drought severity index (Smith et al., 1993). The procedures have also been adopted by the World Meteorological Organization to monitor drought on a worldwide scale (World Meteorological Organization, 1985).
The threshold ranges of deciles used to classify drought conditions are presented in Table 5 (Gibbs and Maher, 1967).
Classification of drought conditions according to the values of the deciles.
An indicator (or indicators) capable of adequately characterizing historical
droughts must be able to capture the following characteristics:
The beginning and the end of the event, which defines its duration. The drought intensity, derived from the index value. The geographical area affected by the drought.
The following guidelines specify the approach proposed in this study to
characterize drought events on the basin scale based on precipitation data
available at each station and how to contrast these results with the
cataloged historical events.
First, following the previous definitions (Sect. 4), precipitation data are used to calculate the four above-described meteorological drought indices (SPI, PN, RAI and DEC) for each station and for different timescales (1-, 3-, 6-, 12-, 24- and 48-month timescales) using the 1951–2000 reference period. Then, according to the criteria presented below, these values are used to detect potential drought events at a given station and at a given time.
In order to aggregate results from all stations of the basin, two indicators are proposed in this study: the Overall Drought Extension (ODE) and Overall Drought Indicator (ODI). The results of these indicators will then be contrasted with historical recorded events to define the best combination of the index and timescale used for the definition of the ODE and ODI indicators.
According to McKee et al. (1993), a drought event occurs at
the station level any time the SPI is continuously negative and the SPI
reaches a value of
Example of the SPI-6, SPI-12 and SPI-24 series at the Chuxiong station, indicating drought periods in orange (lower threshold in red).
Figure 2 shows an example of the SPI-6, SPI-12 and SPI-24 series calculated at the Chuxiong station. Drought periods are colored in orange and the lower threshold that defines their occurrence in red. The influence of the timescale on the number and duration of detected droughts is clearly apparent. It is worth noting that there are periods in Fig. 2 that identify very short droughts (1 or 2 months long), which is due to the identification criteria based on the index values.
In the present study, the above-mentioned principles used to detect drought
events based on the SPI classification (McKee et al., 1993)
have been standardized to be applicable to the other three indices (PN, RAI
and DEC) as follows:
A drought event occurs any time the index is continuously below its normal
value and reaches the moderately dry condition class. The drought is considered to begin when the index first falls below its
normal value. The drought ends when the index exceeds its normal value.
Table 6 summarizes the thresholds for each index that specify the drought
event's start and end criteria, which correspond, respectively, to the limit
of the moderately dry class and to the index normal value. Although the
“normal value” of DEC would be 50 % (which corresponds to the median of
the precipitation records), in this study the drought end criterion
suggested for this index is 60 %, which is the limit between near- and
above-normal conditions
(Jain et al., 2015;
Tsakiris et al., 2007).
Values of the thresholds defining the start and the end of the drought events for each index.
When analyzing directly these meteorological indices, the results only concern each station's surroundings without capturing the patterns of neighboring areas. However, available historical records refer to regional droughts characterized by larger areas that cover several stations.
In order to consider the basin as a whole in the definition of drought
occurrence, duration and intensity, the resulting indices must be
consistently extended to the entire area and then combined in overall
indicators. For that purpose, a regular grid divides the lower JRB into a
For this study, we followed the approaches taken for the definition of drought area indices (Bhalme and Mooley, 1980; Mitchell et al., 1979) and regional drought area indices (Fleig et al., 2011). It is thus considered that a basin-wide event is ongoing when a substantial part of the basin is under drought conditions. It is therefore necessary to identify the portion of the territory for which the calculated index indicates a drought. An indicator to detect drought occurrence on the basin scale has been set up based on the criteria described above to identify an event considering the index values (Fig. 2).
Based on the interpolation of the index, drought events are detected for
each time step and at each grid cell of the described raster. This allows us
to define an indicator, named here the Overall Drought Extension,
expressed as the percentage of the lower JRB area suffering a drought. It is
calculated as the number of cells indicating a drought at a precise date
(N_drought) divided by the total number of cells of the
raster (in this case, N_TOTAL
Regarding the intensity of the droughts, in this study a complementary
indicator is applied to integrate the intensities computed at every grid
cell. The Overall Drought Indicator is defined as the average index
value across the cells under drought conditions at a precise date, as shown
in Eq. (2):
Only cells under drought conditions have been considered to define this indicator. If the ODI had been calculated as an average value for the entire basin (as adopted for instance in Trambauer et al., 2014) higher (or lower) indicator values in a part of the basin may have compensated lower (or higher) indicator values in the rest of the basin, yielding an overall value close to normal precipitation. Therefore, the ODI must always be used together with the ODE: whenever a drought has been detected with the ODE, its overall intensity may be assessed with the corresponding value of the ODI.
In order to support the choice of an index and timescale combination for the definition of the ODE and ODI, an assessment of the quality of the forecasts performed with the different variants is recommended. The hypothesis is that detected drought events (i.e., the forecasts) correspond to the cases when the ODE value exceeds a given threshold, which indicates that a certain area is affected by an event. The temporal coincidence of these forecasts has to be then contrasted with the occurrence of recorded droughts (i.e., the observations). As stated above, Fig. 2 shows that very short index-based events risk being forecasted. In order to avoid an overestimation of droughts, an additional 3-month criterion for beginning and ending forecasted droughts was established: an event will be effectively detected when the ODE value exceeds the threshold for at least 3 consecutive months.
For the matching between forecasts and observations, two monthly series of
events were created (one for the events detected according to the ODE
values and one for the historical events), where for each month either a
“drought” or a “no drought” condition is assigned. Different scores for
contrasting this type of dichotomous forecasts (occurrence vs. no
occurrence) exist: the Peirce skill score, PSS
(Hanssen and Kuipers, 1965; Murphy and Daan,
1985; Peirce, 1884); the Heidke skill score, HSS (Heidke, 1926);
the Gilbert skill score, GSS (Schaefer, 1990); or the odds
ratio skill score, ORSS (Stephenson, 2000). As recommended by
Candogan Yossef et al. (2012), the PSS is used in this
study. For its calculation, the miss rate ( hit: when one detected drought corresponds with an observed drought; false alarm: when a drought appears during a month where no observed event
has occurred; miss: when, during a month where a drought has been observed, no event has
been detected; correct rejection: when, during a month where no drought has been observed,
no drought is detected.
Contingency table of the comparison between forecasts and observations.
The miss rate (
However, high values of the PSS score may be obtained purely by chance,
especially when using only a small number of forecasts. Such is the case of
the present work, where only 13 independent events have been documented
during the 55 years of record keeping. This could lead to overestimating the
goodness of a combination of the index and timescale. Therefore, sampling
uncertainties must be addressed. A statistical test was applied to check if
the calculated PSS values were significantly different from zero, at least
at a 95 % confidence level. Assuming independence of the miss and false alarm rates, the standard error in the Peirce skill score is simply the square
root of the sum of the squared standard errors in the miss and false alarm rates (Stephenson, 2000), as expressed in Eq. (6):
Standard error in estimated miss or false alarm rate calculated using the 95 % score confidence interval as discussed in Agresti and Coull (1998).
Following the previous approach, the series of the SPI, PN, RAI and DEC indices were calculated for different timescales (1-, 3-, 6-, 12-, 24- and 48-month timescales) within the period of 1960–2014. First computed at the 29 stations, these indices were then extrapolated to the rest of the lower JRB. Figure 3 shows the example of the standardized precipitation index for a 6-month timescale (SPI-6) calculated in October 2012 (corresponding to the drought event XII) and spatially distributed at the lower JRB, where brown colors represent regions under drier conditions.
Extrapolated SPI-6 values in October 2012 for the entire lower JRB.
According to the criteria proposed in Table 6, detected drought events were identified based on the index values. Then, the ODE and ODI indicators were calculated for the lower JRB. An example of the resulting ODE and ODI series is shown in Fig. 4 for the SPI-6 and RAI-6 combinations and in Appendix A for all the timescales and indices analyzed, along with the recorded historical droughts shaded in orange.
ODE and ODI values using the 6-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
The objective is to establish a combination of timescale and index that offers an optimum identification of historical droughts. As stated before, the main criteria used to contrast the performance of the forecasts is that a drought event is supposed to happen when the ODE value exceeds a threshold that is to be defined. The combination finally retained should maximize the number of hits and minimize the misses between the forecasts and the observed events.
The 1-month-scale overall indices show rapid fluctuations that correspond to short periods of precipitation deficiency not captured in the catalog of historical droughts. This is mainly due to punctual, large rainfall events that have an important influence in the indices, which may indicate that the drought had ceased when it is not the case (Barua et al., 2011). The use of this timescale is not recommended for drought monitoring since long drought events are hardly identified. The opposite effect occurs when using the 48-month scale. The inertia of the rainfall shortage tendencies may mask shorter droughts and overestimate their durations. Since most of the episodes last 1 year or less (Table 1), they are hardly detected using the 48-month scale. The droughts which occurred from 2009 to 2014 (droughts IX to XIII) illustrate this phenomenon: even if five different droughts have been cataloged, a unique one is detected using the 48-month scale, according to the ODE time series. Therefore, using the 1- and 48-month scales do not provide any substantial information about the occurrence and duration of the droughts and have been excluded from the performance analysis.
For the rest of the timescales (3-, 6-, 12- and 24-month timescale), the ODE threshold indicating the occurrence of a drought is required for the computing of the PSS that will serve as a support for the selection of the best combination of the index and timescale. Traditionally, cross-validation techniques are used to define optimum thresholds, for when within a training subset the threshold maximizing the PSS is identified and validated in a non-overlapping validation subset. However, the authors consider that the limited number of 13 independent events recorded represents an important limitation to a robust and meaningful cross validation. Instead, another approach was followed. A sensitivity analysis was performed using the same threshold across all of the combinations and exploring the effect of varying it in a reasonable range (in this case, from 0.3 to 1 by 0.1 steps). The resulting PSS values are shown in Fig. B1 of Appendix B and an example of PSS for an ODE threshold of 0.4 is presented in Fig. 5. The black bars indicate the statistical error estimates (confidence intervals) at 95 % confidence, due to sampling uncertainties, assessed with the statistical significance test described in Sect. 5.4, which allows indicating whether the score is significantly different from zero.
Example of graphic representation of the PSS results for an ODE
threshold of 0.4, with the black error bars representing the 95 %
confidence interval (
According to the results, most of the PSS confidence intervals do not
include zero, disproving that skill scores could have identified drought
events by chance sampling fluctuations. For some of the indices on the
24-month scale (e.g., RAI-24 for an ODE threshold
Attending to the PSS values (Fig. B1), results show a consistent tendency across all ODE thresholds of higher PSS at the 3- and the 6-month timescales. Moreover, there is no single index that clearly produces better results. Indeed, based on the PSS values and taking into account their uncertainty, there are no statistically significant differences across the different indices for the 3- and 6-month timescales. This indicates that, for these timescales, all the indices perform similarly well in capturing the events, which is consistent with the fact that they all rely on the same type of data (precipitation). PSS results are independent of the specific threshold and thus they are considered robust. However, it is worth mentioning that, in general, higher PSS values for the 3- and the 6-month timescales are produced using ODE thresholds between 0.4 and 0.6.
Regarding the 6-month ODE series (Figs. A5 and A6 of Appendix A), it is
important to highlight some relevant aspects:
All the observed drought events have their corresponding ODE peaks.
Although event VIII has an estimated duration of 3 months, ODE and ODI
results consistently show a longer drought. The exact period of this drought
is not well defined as indicated in the catalog, leaving room for a longer
duration of the real episode. In general, all the indices are well correlated, identifying most of the
recorded droughts. Several droughts are consistently detected between event I (1962) and II
(1979) even if no drought has been chronicled (false alarms). This may
correspond to the above-mentioned scarcity of reliable information on
droughts prior to 1980. The drought events IX, X, XI, XII and XIII are well captured. As shown in
Figs. A5 and A6, the different events during this period (2009–2014)
match with the consecutive increases in the ODE values for all the indices
(DEC, PN, RAI, SPI). However, the 6-month series of ODE suggest some false positive detections:
more drought events than the observed are calculated.
Regarding the 3-month ODE series (Figs. A3 and A4), results suggest an
overestimation of the number of detected events, as sometimes several
detected events combine into one (longer) observed event. The 6-month timescale appears as more appropriate.
In summary, according to the ODE series presented in Appendix A and to the forecast verification carried out with the Peirce skill score (Appendix B), it seems that the best timescale for the identification of droughts is at 6-months. Results show an equally effective performance of the ODE series for all the indices. However, the risk of false positives must be addressed carefully, as the observation record likely misses events, in particular between 1962 and 1979.
Despite the good performance shown by the overall indicator ODE in drought detection, caution is advised. In particular, the choice of meteorological indices as a basis for the calculation of the ODE and ODI can lead to errors when assessing drought occurrence. Temperature variability, not considered here, can play a significant role in the onset of agricultural drought. Meteorological indices may not be fully capable of capturing the impacts on water scarcity and could be complemented with other types of indices, such as agricultural or hydrological. The same approach proposed in this study is recommended using more comprehensive indices in order to better capture the complex drought processes.
The performance assessment of the ODE indicator to detect droughts relies basically on the comparison with the historical events cataloged in this study. The search and compilation of this information from different data sources, often scarce and ambiguous, represents a challenge. Different information sources often provide only partial information for one episode, and for some events the differences in the available information complicate the harmonization of data. As a result, the accuracy of the collected information may impact on the applicability of the developed methodology.
This study aims at applying overall drought indicators representing the drought status within the entire lower JRB investigation area. This work represents an attempt at building a tool for drought monitoring and risk management purposes on the basin scale. It is based on established meteorological indices for the identification of droughts and a method for a catchment-wide drought assessment and characterization, which is compared to historical drought events of the lower JRB.
The information used for the identification and characterization of major historic droughts was compiled from different sources. A total of 13 major droughts between 1960 and 2014 were identified in the lower JRB and cataloged using a web-based registration platform, allowing for a comparison of the different events.
Drought indices typically assess local water deficits while available historical records usually refer to regional droughts. To overcome this problem, two drought area indicators, the Overall Drought Extension and the Overall Drought Indicator, have been used to characterize the occurrence and intensity of an event within a specific investigation area. These indicators are based on four common meteorological indices on different timescales: the standardized precipitation index, the rainfall anomaly index, the percent of normal precipitation and the deciles index. By relying exclusively on precipitation, the proposed procedure serves as a basis for further studies in other regions where only precipitation data are available.
The performance of the ODE in drought detection has been assessed by contrasting the results of this indicator with historical recorded events, offering promising results. It seems that the best results are independent of the index used and produced using the 6-month timescale. Although results suggest the same patterns for all ODE thresholds, it has been noticed that the highest PSS values are produced for thresholds between 0.4 and 0.6, which can be defined as a trigger to detect the occurrence of a drought in the lower JRB.
Considering the challenge that the compilation of historical drought information represents and the identified limitations, this is a good method for the monitoring of drought episodes within an entire catchment. The use and contrast of drought indicators on the basin scale with historical collected information represent the main innovative aspects of this study. Since meteorological droughts are the first stage in the progression of subsequent agricultural or hydrological droughts, this methodology could be used to activate a management response for a drought event, which starts at a specific threshold value. Additionally, this methodology can be used to complete lacking information on droughts' duration, geographical extension or intensity.
Precipitation data analyzed in the current study are available
from the China Meteorological Data Service Center (CMA) and can be downloaded from its
data sharing service system (
The series of the Overall Drought Extension and the Overall Drought Indicator have been calculated for the 1-, 3-, 6-, 24 and 48-month timescales. Graphic results are presented in Figs. A1 to A12 below.
ODE and ODI values using the 1-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 1-month timescales of DEC and PN indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 3-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 3-month timescales of DEC and PN indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 6-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 6-month timescales of DEC and PN indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 12-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 12-month timescales of DEC and PN indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 24-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 24-month timescales of DEC and PN indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 48-month timescales of SPI and RAI indices, compared with the 13 detected historical droughts (in orange).
ODE and ODI values using the 48-month timescales of DEC and PN indices, compared with the 13 detected historical droughts (in orange).
Graphical PSS results for different index and timescale combinations, for a range of ODE thresholds between 0.3 and 1. The black error bars represent the 95 % confidence interval.
The authors declare that they have no conflict of interest.
This study is based on the project “Jinsha River Basin (JRB): Integrated
Water Resources and Risk Management under a Changing Climate” funded by the
Swiss Agency for Development and Cooperation (SDC) and has been supported
financially by the International Science and Technology Cooperation Program of
China (grant no. 2014DFA71910), the National Natural Science Foundation of
China (grant no. 51709009), and the Fundamental Research Funds for Central
Public Welfare Research Institutes (grant no. CKSF2016024). The authors would
like to acknowledge and thank the Bureau Of Hydrology (BOH) and the
Changjiang River Scientific Research Institute (CRSRI) from the Changjiang
Water Resources Commission for providing the data and the enterprise Ernst
Basler