Next Article in Journal
Sand Dune Albedo Feedback
Next Article in Special Issue
Drought Prediction System for Central Europe and Its Validation
Previous Article in Journal
Regional Scale Sea Cliff Hazard Assessment at Sintra and Cascais Counties, Western Coast of Portugal
Previous Article in Special Issue
ENSO Index-Based Insurance for Agricultural Protection in Southern Peru
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Performance of Six Drought Indices in Characterizing Historical Drought for the Upper Blue Nile Basin, Ethiopia

1
IHE Delft Institute for Water Education, Westvest 7, 2611 AX Delft, The Netherlands
2
National Drought Mitigation Center, University of Nebraska-Lincoln, 815 Hardin Hall, 3310 Holdrege St., Lincoln, NE 68583, USA
3
Department of Civil Engineering, Addis Ababa Institute of Technology, King George VI St., Addis Ababa 1000, Ethiopia
4
Department of Water Science and Engineering, Vrije Universiteit Brussel, Boulevard de la Plaine 2, 1050 Ixelles, Belgium
5
Water Resources Section, Delft University of Technology, 2600 GA Delft, The Netherlands
6
Water Problems Institute of RAS, Gubkina st., 3, Moscow 119333, Russia
*
Author to whom correspondence should be addressed.
Submission received: 17 October 2017 / Revised: 22 February 2018 / Accepted: 22 February 2018 / Published: 27 February 2018
(This article belongs to the Special Issue Drought Monitoring and Prediction)

Abstract

:
The Upper Blue Nile (UBN) basin is less-explored in terms of drought studies as compared to other parts of Ethiopia and lacks a basin-specific drought monitoring system. This study compares six drought indices: Standardized Precipitation Index (SPI), Standardized Precipitation Evaporation Index (SPEI), Evapotranspiration Deficit Index (ETDI), Soil Moisture Deficit Index (SMDI), Aggregate Drought Index (ADI), and Standardized Runoff-discharge Index (SRI), and evaluates their performance with respect to identifying historic drought events in the UBN basin. The indices were calculated using monthly time series of observed precipitation, average temperature, river discharge, and modeled evapotranspiration and soil moisture from 1970 to 2010. The Pearson’s correlation coefficients between the six drought indices were analyzed. SPI and SPEI at 3-month aggregate period showed high correlation with ETDI and SMDI (r > 0.62), while SPI and SPEI at 12-month aggregate period correlate better with SRI. The performance of the six drought indices in identifying historic droughts: 1973–1974, 1983–1984, 1994–1995, and 2003–2004 was analyzed using data obtained from Emergency Events Database (EM-DAT) and previous studies. When drought onset dates indicated by the six drought indices are compared with that in the EM-DAT. SPI, and SPEI showed early onsets of drought events, except 2003–2004 drought for which the onset date was unavailable in EM-DAT. Similarly, ETDI, SMDI and SRI-3 showed early onset for two drought events and late onsets in one-drought event. In contrast, ADI showed late onsets for two drought events and early onset for one drought event. None of the six drought indices could individually identify the onsets of all the selected historic drought events; however, they may identify the onsets when combined by considering several input variables at different aggregate periods.

1. Introduction

Drought is abnormally dry condition sustained for a prolonged period. Unlike other natural disasters, such as floods and tornadoes, drought’s impacts do not manifest immediately. However, in its aftermath, crops wilt, livestock dies, rivers and reservoirs dry up, and socioeconomic sectors are affected [1,2]. Drought can be characterized by its severity, duration, and areal extent. In the past few decades, a number of devastating drought events have been occurred. For example, in 2009, severe drought affected different parts of the globe, but more people were affected in eastern Africa than anywhere else. About 20 million people were estimated to have been suffering from food insecurity in eastern Africa [3]. The year 2009 was also recorded as a severe drought year in Ethiopia; six million people were affected and needed food aid from international emergency services [3]. Other historic drought events have also occurred in Ethiopia in the past few decades, most notably 2003/2004 and 1984/1985 [4]. These drought years were recorded as the worst in the history of the country, in which millions of people lost their lives and many more people were forced into further destitution [5]. Natural drought events cannot be prevented, but through monitoring, forecasting, and early warning, drought mitigation can be better managed, and impacts can be reduced.
Drought is mainly categorized into four major types: meteorological, agricultural, hydrological, and socioeconomic [1,2]. Meteorological drought is characterized by a prolonged deficit of precipitation from its long-term average [1]. Agricultural drought is characterized by deficits of total soil moisture and results mainly from the shortage of precipitation. Hydrological drought is related to the impacts of persistent shortage of precipitation on surface and/or subsurface water supply (i.e., streamflow, reservoir and lake levels, and groundwater) [1,6,7,8,9,10]. Socioeconomic drought is associated with the impacts of meteorological, agricultural, and hydrological droughts on the socioeconomic sectors. Drought monitoring employs widely used drought indices, such as the Palmer Drought Severity Index (PDSI) [11,12], Standardized Precipitation Index (SPI) [12,13,14,15], Standardized Precipitation Evapotranspiration Index (SPEI) [16], Standardized Runoff Index (SRI) [17], Surface Water Supply Index [18], Decile Index [19], and Vegetation Drought Response Index [20]. Detailed reviews of drought indices can be found in Zargar et al. [21], and regional specific drought studies can be found in Keyantash and Dracup [22], Ntale and Gan [23], Morid et al. [24], Barua et al. [25], Dogan et al. [26], and Trambauer et al. [27].
Because of the inherent complexity of the drought phenomenon and variation in the available good quality data, the performance of drought indices in characterizing historic drought events may vary from place to place [28]. Identifying an appropriate drought index for a specific region is crucial for mitigating and preparing for drought-related disasters. Comparative studies of drought indices have been carried out in many regions or river basins. Jain et al. [28] compared the suitability of six drought indices in the Ken River basin, India, and indicated that the Effective Drought Index (EDI) is a more suitable drought index for the study basin. Morid et al. [24] compared the performances of seven drought indices and recommended using the EDI and SPI for drought monitoring in the basin. Their study also indicated that the EDI was more responsive to drought and performed better than the SPI. Okpara and Tarhule [29] evaluated the performance of three drought indices in the Upper Niger basin based on six decision criteria. They reported that the SPI ranked first among meteorological drought indices. Barua et al. [25] evaluated the performance of five drought indices for the Yarra river catchment in Australia using five subjective decision criteria (robustness, tractability, sophistication, transparency, and extendibility). They found that the Aggregate Drought Index (ADI) was superior and preferable to other indices for the study catchment. Many similar studies have been conducted in other basins [30,31,32,33]. In addition, drought related studies in Ethiopia have been reported in the literature, especially focusing on areas in the northern and eastern parts of the country [4,34,35], including the Upper Blue Nile (UBN) basin [5,34,36]. However, the UBN lacks a basin-specific drought monitoring system.
The main objective of this study is to evaluate and compare the performance of multiple drought indices for indicating the onset, severity, and duration of historical drought events in the Upper Blue Nile (UBN) basin. The six drought indices used are the Standardized Precipitation Index (SPI), Standardized Precipitation Evaporation Index (SPEI), Evapotranspiration Deficit Index (ETDI), Soil Moisture Deficit Index (SMDI), Aggregate Drought Index (ADI), and Standardized Runoff-discharge Index (SRI). These drought indices were calculated using observation data and hydrological model results. This study is expected to help in establishing basin-specific drought monitoring by advancing the understanding of how different drought indices perform in identifying historical drought events.

2. Study Area and Data

2.1. Study Area: Upper Blue Nile Basin (UBN)

The UBN basin is located in the north-western region of Ethiopia between 7°40′ N and 12°51′ N latitudes and 34°25′ E and 39°49′ E longitudes (Figure 1), and covers a total area of 176,000 km2 upstream of the Ethiopia-Sudan border [37]. The topography of the basin signifies two distinct features: the highlands with rugged mountainous areas in the central and eastern part of the basin, and the lowlands in the western part of the basin. The altitude in the basin ranges from 490 to 4260 m. While the highlands are the main source of water, the lowlands have expanses of flat lands through which the river flow travels from the highlands to the lower riparian countries (Sudan and Egypt). Tekleab et al. [38] reported that the mean annual temperature (data years 1995–2004) ranges from 13 °C in south-eastern parts to 26 °C in the south-western part near the Ethiopia-Sudan border. The annual rainfall ranges from 780 to 2200 mm, with the highlands having the highest rainfall (ranging from 1500 to 2200 mm) and the lowlands receiving less than 1500 mm [37,39,40]. Following the seasonality of the rainfall, 82% of the annual flow occurs during the main rainy season in June to September [41].
The land cover in the basin is dominated by dry-land crops, pastures, savannah, grassland, woodland, water bodies, and sparsely vegetated plants [42]. Volcanic rock and Precambrian basement rock are the most widely available geological formations in the basin, and small areas are covered by sedimentary rock [36]. The dominant soil types are Leptosols and Alisols (21%), Nitosols (16%), Vertisols (15%), and Cambisols (9%) [43].
Administratively, 47% of the area of the UBN basin is in the region of Amhara, 31% in the region of Oromia, and 23% in the region of Benishangul Gumuz. In the earlier division of the country (before 1992), the basin covers substantial parts of Gondar, Wollo, Gojjam, Welega, and Shoa provinces of Ethiopia. The drought onset and number of affected population data in the basin prior to 1992 are associated with the former political administrative provinces (Note: Ethiopia was divided into provinces until they were replaced by regions (“Kililoch”) and chartered cities in 1992).

2.2. Data Used

2.2.1. Historical Drought Events

The history of drought in the basin is poorly documented. However, major historical drought events within the study period (1970–2010) were identified from previous studies and Emergency Events Database (EM-DAT), the international disaster database (http://www.emdat.be/database) [5,33,44]. According to previous studies, 1973–1974, 1983–1984, 1994–1995, and 2003–2004 were reported as the major drought years in Ethiopia [5,33]. EM-DAT (Table 1) includes these drought events as well, but reports the 1973 and 1994 events as part of longer drought periods from 1973 to 1978 and from 1989 to 1994, respectively. This possibly reflects the occurrence of multiple drought events during these years in different parts of the country. For evaluating the six drought indices (Section 4.2), the EM-DAT drought periods are taken as a reference, while keeping in mind their reporting of prolonged droughts. The EM-DAT database is compiled from various sources, including UN agencies, non-governmental organizations, insurance companies, research institutes, and press agencies. The available information for the selected drought events is presented in Table 1. A previous report shows that the quality and reliability of the EM-DAT is complementary with other disasters data reporting organizations (e.g., NatCat and Sigma) since the EM-DAT data are derived mainly from humanitarian agencies and development organizations [45]. Its main limitation is a lack of data at a finer spatial resolutions (e.g., at county, or district level) so that the particular drought start date (Table 1) is assumed to be the same for all of the provinces. The year 1983–1984 stands out as the drought year with the highest percentage (22%) of the population affected. For other drought events, this percentage is 16% or less. The provinces indicated in bold (Table 1) are situated (partly) within the UBN basin. For example, according to the previous division of provinces (before 1992), Gondar includes the northern part of the basin, whereas Wollo and Shoa include eastern and northeast, and southeast parts of the basin, respectively. According to the current regional classification, Amhara covers the central, northern, northeastern, and northwest parts, and Oromia covers the southern, southeast, and southwest parts of the UBN basin. The eastern and some of the central parts of the basin have not been reported drought affected area in EM-DAT. Both the start and end dates of the historic drought events are shown in Table 1, but the months of the end dates are not indicated. Hence, the evaluation of the drought indices in identifying drought events focuses on the onset of the drought events.

2.2.2. Actual Evapotranspiration (ET) and Soil Moisture Data

We used a process based semi-distributed hydrological model (Soil and Water Assessment Tool; SWAT) to simulate the actual evapotranspiration and soil moisture time series data. The SWAT model was setup using mainly daily rainfall, and temperature; and calibrated and validated using river discharge data from 1970 to 2010. The monthly calibration results in terms of coefficient of determination (r2) and Nash-Sutcliffe efficiency (NSE) between observed and simulated river discharge in five stations vary from 0.83 to 0.93 (for r2) and 0.84 to 0.91 (for NSE). The Hargreaves method [46] was used for potential ET in SWAT. Weather data of relative humidity, solar radiation, wind speed, and sunshine hour available from 10 stations are used for updating the SWAT database. Other input data (Digital Elevation Model (DEM), land use, and soil maps) were collected from difference sources. The 90 m spatial resolution DEM data were obtained from the Shuttle Radar Topography Mission (SRTM) (http://www.cgiar-csi.org/data/srtm-90m-digital-elevation-database-v4-1). Land use classification data of the year 2000 at a spatial resolution of 5 km by 5 km were acquired from the Ministry of Agriculture and Rural Development (MARD) [47] of Ethiopia. The basin is classified into 32 land use classes. Cultivated land is the dominant land cover in the basin. The soil map and physical and chemical properties at different soil layers were obtained from the Ministry of Water, Irrigation, and Electricity of Ethiopia [48]. Actual evapotranspiration and soil moisture from the model output were used to calculate evapotranspiration deficit index (ETDI) and soil moisture deficit index (SMDI), respectively. Previous studies have shown the capability of SWAT to model the hydrological processes in the UBN basin (e.g., [42,49,50,51,52,53]).

2.2.3. Rainfall and Temperature Data

Weather data, such as precipitation, and maximum and minimum temperature, were collected from the National Meteorological Agency (NMA) of Ethiopia. Data from 34 stations with records from 1970 to 2010 were used in this study. According to the NMA’s classification, the majority of the stations are class three, which only record precipitation and temperature data. In this study, ten “class I” stations that have additional observations of relative humidity, wind speed, and sunshine hour, were also included for updating the SWAT model database. The data were on daily time step and the location of the stations is shown in Figure 1. The quality of the rainfall data has been checked in the previous study [5] and the same stations were used in this study. The rainfall and temperature data were used to calculate the SPI and SPEI indices as well as to run the SWAT model.

2.2.4. River Discharge Data

The river discharge data measured at five selected gauging stations (Gilgel Abay, Ribb, Gumara, Abbay at Kessie, and Ethiopia-Sudan border) were obtained from the Ministry of Water, Irrigation, and Electricity of Ethiopia (Figure 1). The data are on a daily time step (except at the Ethiopia-Sudan border, where only monthly data were available). The river discharge data at these selected stations were used to calculate the time series values of SRI, and to calibrate and validate the SWAT model. The data length for the period from 1989 to 1996, and 1997 to 2002 were used to calibrate and validate the model respectively.
Table 2 summarizes the input data used in this study. The details about the data can be referred in Section 2.2.1, Section 2.2.2, Section 2.2.3 and Section 2.2.4.

3. Drought Indicators

The drought indices used in this study collectively characterizes meteorological, agricultural, and hydrological drought types. The SPI (based on precipitation) and SPEI based on precipitation and temperature) characterize primarily meteorological drought. The ETDI (based on evapotranspiration) and SMDI (based on soil moisture) characterize agricultural drought. Similarly, the SRI, which is based on discharge, characterizes the occurrence of hydrological drought. ADI aggregates multiple input data sets that represent meteorological, agricultural and hydrological variables. These drought indices have been selected because of their less data intense, potential for characterizing drought at different aggregate period, and describe the dry and wet periods. These drought indices are described below.

3.1. Meteorological Drought Indicator

Standardized Precipitation Index (SPI): The SPI utilizes monthly time series of rainfall data preferably (30 years or more) to compute its value. SPI measures the standardized departure that the observed precipitation on a given time scale deviates from the long-term average precipitation [22]. The computation fits a probability density function to the frequency distribution of precipitation summed over the period of interest then transfers to a normal distribution [53,54]. The details about the calculation procedure of the SPI can be found in Bayissa et al. [5]. Based on SPI values, meteorological droughts are classified into different levels of severity, as shown in Table 3.
Standardized Precipitation Evaporation Index (SPEI): The SPEI differ from SPI in that it incorporates the effect of potential evaporation (calculated based on temperature) in addition to precipitation [16,27,55]. The calculation procedure of the SPEI is also similar to that of the SPI, except that it uses the difference between precipitation and potential evaporation instead of precipitation alone. Like the SPI, the SPEI is calculated at different time scales (e.g., 1-, 2-, and 3-month). In this study, a log-logistic distribution is applied because it fitted the best in the majority of rainfall stations that were used [5]. Negative SPEI values indicate dry condition due to less precipitation and/or higher potential evaporation (dry conditions) as compared to the historical mean. Unlike the SPI, the SPEI has not been widely applied and tested in the UBN basin [5].

3.2. Agricultural Drought Indicators

Evapotranspiration Deficit Index (ETDI): The ETDI is based on the anomaly of water stress to its long-term average [9], in which monthly water stress is defined using potential and actual evapotranspiration (Equation (1)).
WS = PET AET PET ,
The WS ranges from 1 (no evapotranspiration) to 0 (evapotranspiration occurring at the same rate as PET). Next, monthly water stress anomaly (WSA) is calculated as:
WSA i , j = MWS j WS i , j MWS j minWS j × 100 , if   WS i , j MWS j ,
WSA i , j = MWS j WS i , j maxWS j MWS j × 100 , if   WS i , j > MWS j ,
where MWSj is the long-term median of water stress of month j, maxMWSj is the long-term maximum water stress of month j, minWSj is the long-term minimum water stress of month j, and WS is the monthly water stress. The subscripts i and j are used for years and months, respectively. The water stress anomaly ranges from −100 to +100 for each month, indicating very dry to very wet conditions [9].
Finally, ETDI is calculated using Equation (4). In the original formula developed by Narasimhan and Srinivasan [9], WSAj is divided by 50 to scale the ETDI between −4 and 4 to compare with PDSI values. Here, as suggested by Trambauer et al. [27], the WSAj values are divided by 100 to scale the ETDI between −2 and 2 to compare with SPI, SPEI, and SRI values.
ETDI j = 0.5 ETDI j 1 + WSA j 100 ,
Soil Moisture Deficit Index (SMDI): The SMDI is calculated in the same way as ETDI, but with the available soil water content in the soil profile [9]. First, the median, maximum, and minimum values for each month were extracted using soil moisture time series. The median was chosen instead of the mean because it is less affected by the outliers. The SMDI values (deficit or excess) for the 39 years (1970–2010) were calculated using Equations (5)–(7).
SD i , j = SW i , j MSW j MSW j minSW j × 100 , if   SW i , j MSW j ,
SD i , j = SW i , j MSW j maxSW j MSW j × 100 , if   SW i , j > MSW j ,
where SDi,j is the soil water deficit (%) ranging from −100 (very dry condition) to +100 (very wet condition); SWi,j is monthly soil water available in the soil profile (mm); and, MSWj, maxSWj, and minSWj are long-term median, maximum, and minimum available soil water in the soil profile (mm), respectively.
Thus, the SMDI in any given month is determined by:
SMDI j = 0.5 SMDI j 1 + SD j 100 ,
SMDI ranges from −2 to +2, with negative values indicating drought.

3.3. Hydrological Drought Indicator

Standardized Runoff-discharge Index (SRI): The SRI is based on river discharge and its computation procedure is similar to that of SPI [17]. The gamma distribution is also used here to fit the river discharge data. The SRI and SPI drought categories are similarly calculated. We used river discharge data at the five selected gauging stations from the study area.

3.4. Aggregate Drought Index (ADI)

The ADI [22] is a multivariate drought index that aggregates the three types of drought (meteorological, agricultural, and hydrological) through considering several input variables representing these three drought types (e.g., rainfall, evapotranspiration, soil moisture, and stream flow). In this study, rainfall, evapotranspiration, and soil moisture are used to derive a single (aggregated) index, ADI, using Principal Components Analysis (PCA). The z-score values (1-month time scale) of these three input variables were calculated to avoid inconsistency as a result of different units of measurement. PCA constructs a symmetric correlation matrix (p × p, where p is the number of variables) between the z-score values of the input variables. PCA transform the original p-variables data set to a number of uncorrelated (principal) components zj (1 < j ≤ p) [25] using Equation (8).
Z = XE ,
where Z is the (n × p) matrix of PCs (i.e., uncorrelated components), in which n is the number of observations, X is the (n × p) matrix of the observation data, and E is the (p × p) matrix of eigenvectors of the correlation matrix.
The ADI was calculated for the first PC (PC1) normalized by the standard deviation. PC1 described more than 65% of the variation in the input data.
ADI i , k = Z i , 1 , k / σ k ,
where ADIi,k is the ADI value for month k in year i, Zi,1,k is the first PC for month k in year i, and σk is the sample standard deviation of Zi,1,k for all years I and months k. The ADI was calculated for each of the station locations and each month. The threshold of ADI was determined based on McKee classifications of SPI drought categories (Table 3), which were developed for stations across Colorado in the United States. According to McKee, mild drought occurs 24% of the time, moderate drought occurs 9.2% of the time, severe drought occurs 4.4% of the time, and extreme drought occurs 2.3% of the time. In this study, we determine the ADI drought category through constructing cumulative distribution of the time series of ADI across the UBN basin. The drought categories of ADI at the corresponding percentile are shown in Table 4.ADI values of −0.96 or lower indicate a drought at a different severity level (Table 4).
The summary of the drought indices used in this study is shown in Table 5.

4. Methods

4.1. Correlation between Drought Indices

First, we derived monthly time series of drought indices that are defined in the previous section. The drought indices were calculated at the corresponding locations of each meteorological station. The spatial extent of SWAT based evapotranspiration and soil moisture were in hydrologic response units (HRUs). The HRUs are the smallest spatial units of the model defined based on lamping similar land uses, soils, and slopes [56]. In SWAT modelling, the basin was divided into two major parts: upstream and downstream of the Lake outlet. This division has been made to exclude the Lake Tana from the modelling processes of the basin due to the lack of appropriate data to characterize the lake. The upstream and downstream parts are further divided into 14 and 139 sub-basins and 104 and 1027 HRUs, respectively. The HRUs were defined based on percentage combinations of land use, soil, and slope. In this study, 10% land use, 20% soil, and 10% slope threshold combinations were adopted, based on the findings of Setegn et al. [50]. The HRUs values of evapotranspiration, and soil moisture at the locations of each meteorological station were extracted and used to calculate ETDI and SMDI drought indices. For all indices, data from 1970 to 2010 are used. For SPI, SPEI, and SRI indices were derived for five different aggregation periods, namely 1-, 3-, 6-, 9-, and 12-months. Thus we have a total of 18 indices; five each from SPI, SPEI, and SRI and one each from ETDI, SMDI, and ADI. Next, Pearson correlation coefficients were derived for paired time series values of drought indices. Each drought index is paired with every other drought indices, resulting in 18 by 18 correlation coefficient matrix.

4.2. Comparison of Drought Indices Based on Drought Onset, Duration, and Severity

The comparison of the drought indices based on the drought characteristics, such as percentage of drought months, maximum drought intensity, and drought duration was analyzed for the selected 16 stations that are representing majority of the study area. The percentage of drought months was calculated by taking the percentage of the ratio of the total number of months that show drought condition (including mild, moderate, severe, and extreme drought) with the total events in the study period. The maximum drought intensity represents the smallest value of the drought index within the study period (1970–2010). The average value of the maximum intensity of the different aggregate periods (1-, 3-, 6-, 9-, and 12-month) was considered for SPI, SPEI, and SRI, and was compared with the self-defined single time scale indices (i.e., ETDI, SMDI, and ADI). The drought duration is defined as the consecutive months that show the drought condition (below normal conditions). Drought duration is also considered as one of the evaluation criteria to compare the drought indices for the selected stations.
Other drought indices comparison criteria (i.e., drought onset, duration, and severity) have been used by dividing the UBN basin into upper, middle, and lower parts, representing the areas upstream of the river gauging stations Abbay at Bahirdar, Kessie, and Ethiopia-Sudan border, respectively. The drought indices were calculated using areal averages of the variables, namely rainfall, temperature, soil moisture, and actual ET based on the simple arithmetic mean technique. The drought severity, duration, and onset values were extracted for the known historic drought year (1973–1978, 1983–1984, 1989–1994, and 2003–2004). The time series values of the drought indices for these selected drought events were analyzed to calculate the drought severity. The consecutive negative values of each drought index were considered to quantify the different drought characteristics. For example, mean intensity (M) and maximum intensity (Mmax) of drought are the average and maximum values within the consecutive negative drought index values, respectively. Drought duration (D) represents the time span of the consecutive negative index values. Drought severity is defined as the product of drought duration (D) and the mean intensity (M) [1,8,9]. The severity is interpreted as near normal, moderate, severe, and extreme drought based on the drought category classification of each drought index, as described in Section 3. Based on drought severity and duration, each drought index was compared whether they characterize and indicate the severity and persistency of drought years in the basin. Moreover, the drought onset (starting month) that was estimated by each index was compared with similar data obtained from EM-DAT (Section 2.2.1).

5. Results and Discussion

5.1. Time Series of the Drought Indices

The time series (1970–2010) of the drought indices were produced for all the stations considered. Only the result obtained at Gondar, and Debremarkos stations were shown in Figure 2 for further discussion in this section. These two stations are representative of the upstream and downstream parts of the study basin. Figure 2 depicts the performance of the drought indicators on indicating the historic drought events. In general, the majority of the drought indices indicated the historic drought years (1973–1974, 1983–1984, 1994–1995, 2003–2004, 2008–2009), except observing a frequent jump between the drought and normal condition at a lower time scale (1- and 3-month). The persistency of the historic drought years was indicated at a larger time scales (e.g., 6-, 9-, 12-month). The meteorological drought indices showed severe drought conditions in the year 1991–1992 in the majority of the stations; however, the severity of this drought condition has not been reflected by other drought indices. The agricultural drought indices showed a similar patter with the meteorological drought indices of smaller time scales (1- and 3-month). Moreover, the agricultural drought indices showed smaller magnitude of the drought severity in some of the stations (e.g., Gondar) as compared to other drought indices. The hydrological drought indices indicated the persistency of the drought condition even at lower time scale as compared to the other indices. The possible reason is perhaps the fluctuation of river flow is very gradual as compared to the other input variables, such as rainfall, evapotranspiration, etc. The majority of the drought indices indicated the severity of the historic drought condition at different severity level that ranges from mild to extreme drought condition. Thus, each drought index could potentially be used to assess some of the historic drought conditions in the UBN basin.

5.2. Correlation between Drought Indices

The Pearson’s correlation coefficient matrix was developed for the six drought indices for each meteorological station. The results obtained at Gondar (Table 6a) and Debremarkos (Table 6b) stations are presented as an example. A correlation matrix of 18 rows by 18 columns was created to investigate the relationship between the drought indices (Table 6). SPI-3 showed a relatively high correlation coefficient with a majority of other indices, except SPI among itself at other aggregation periods. For example, SPI-6 correlates better with SPI-9 (0.77) and SPI-12 (0.72) than SPI-3 does with SPI-9 (0.73) and SPI-12 (0.68) at Gondar station. However, SPI-3 correlated better with SPI-1 (0.65) than SPI-6 (0.53), SPI-9 (0.48), and SPI-12 (0.45). Table 6 shows that the correlation coefficients between SPI and SRI or SPEI and SRI consistently increase as the aggregation period for SPI and SPEI increases. For example, SPI-12 and SPEI-12 correlated better with the SRI than the SPI and SPEI at shorter aggregation periods. In general, longer aggregation periods of SPI and SPEI better correlate with SRI. The highest correlation for SRI was reached when both indices are aggregated over 12 months (correlation of about 0.55). Compared to correlations between other pairs of indices, correlations involving SRI are generally low. SRI is based on river flow, which is a combination of surface flow, interflow and base flow (or ground water flow), and is expected to have certain lag time with meteorological and agricultural droughts. This might be one of the possible reasons why the SRI is not correlated well with other drought indices. Table 6a,b, first five columns, show that there is a strong correlation between drought indices of the same aggregation period. SMDI and ETDI (agricultural drought indices based on monthly data) have a higher correlation with SPI-3 and SPEI-3 (meteorological drought indices based on 3-month aggregated data) than with the other indices or aggregate periods. This indicates that in the UBN basin the meteorological drought indices at 3-month aggregate period explains best the agricultural drought indices and perhaps the occurrence of an agricultural drought. Similar results were obtained for Debremarkos (Table 6b), and other meteorological stations (can be referred from the Supplementary Material Table S1).
The statistical significance of the correlation coefficient values of each station was tested using t-test distribution. The result of the test is shown in Table 7 for each drought index. The result shows that, for the majority of the drought indices, the p-values are less than or equal to 0.01 (shown in green in Table 7), meaning they are statistically significant at 99% or higher confident level. Few other correlation values (e.g., SPEI-1 and SPEI-12 with SRI at different aggregation periods, and ADI with all SPEIs) are significant at 90% or higher confidence level (i.e., p ≤ 0.1) (shown in yellow in Table 7). However, the correlation values of ADI with the majority of the other drought indices (e.g., with SPI-6, SRI-9, ETDI, and SMDI) are statistical insignificant at 90% confident level (i.e., p > 0.10), and needs further investigation for future studies.

5.3. Comparison of Drought Indices Based on Drought Characteristics

The comparison of the drought indices based on drought characteristics, such as the percentage of drought months, maximum drought intensity, and drought duration was analyzed and used as additional comparison criteria for each index. The results obtained at the selected 16 stations that are representing the majority of the study area were discussed in this section. The percentage of the drought months represents the proportion of the total number of drought months (including mild, moderate, severe, and extreme droughts) within the study periods (1970–2010) of each station. The resulting graph is shown in Figure 3 and each spider web represents the value of each drought index across the selected stations. In general, the percentages of the drought months of the hydrological (SRI), agricultural (ETDI and SMDI), and ADI depict relatively larger values as compared to the meteorological drought indices (i.e., SPI, and SPEI) in majority of the stations for 1-month, 3-month, 6-month, and 12-month temporal scales. However, meteorological indices showed a smaller improvement (increase) of the percentage of drought months as the time scale increases in majority of the stations. The comparative analysis of the percentage of the drought months was extended separately for each drought severity categories, such as mild, moderate, severe, and extreme drought. Only the result obtained for severe and extreme drought categories are shown in Figure 4. Interestingly, a similar result (as above) was observed only for mild and moderate drought conditions (not shown). However, an opposite result that shows a relatively larger percentage of drought months was observed by meteorological drought indices (Figure 4) for severe and extreme drought condition as compared to other drought indices in the majority of the stations. This perhaps shows meteorological drought indices are influenced by the variability of the amount of rainfall during drought event. ADI also showed large percentage of the drought months, which may indicate the dominance of the rainfall in ADI index than the other two input variables (evapotranspiration and soil moisture).
The maximum duration of drought was also considered as another comparison criteria of drought indices for all time scales (Figure 5). The number of consecutive drought months that shows drought severity values of -1 and less (representing moderate, severe, and extreme) were counted for the study period. Evaluating drought indices based on the maximum duration of drought helps to evaluate their performance on indicating the persistency of historic drought events. In general, the result shows an increase in the maximum drought duration as the time scale increases for SPI, SPEI, and SRI in the majority of the stations. This shows persistent drought is indicated at larger time scales (e.g., 12-month) than small time scales (e.g., 1-month). ETDI, SMDI, and ADI shows less persistent drought in the majority of the stations and the result corresponds with the persistency indicated by SPI and SPEI at lower time scales (1-, and 3-month). This shows the comparability of the agricultural drought with the meteorological drought indices at lower temporal scales. The comparison of maximum duration of drought of SPI, SPEI, and SRI at the same time scales (e.g., 3-month) shows the performance of SRI on indicated relatively a large drought duration month than SPI, and SPEI. The time series of SRI is derived based on river flow (resulted from the catchment process) so that SRI is less affected by the extreme wet event in between a particular historic drought years. Interestingly, SPEI also showed relatively larger drought duration months as compared to SPI at the same time scale. The possible reason might be the use of additional input variables (e.g., potential evaporation, in the case of SPEI) could help to show the persistency of drought. The maximum duration of drought indicated by ETDI, SMDI, and ADI is relatively less in the majority of the stations. For several aggregate periods and for several locations, SPI, SPEI, and SRI showed maximum consecutive drought duration of more than 12-months. It may, therefore, be of interest to analyze drought indices at higher time scales (>12-months) in future studies.
Similar approach/method that was followed in this study was also tested in other watershed to compare drought indices [28,30]. In general, the results of these studies showed the capability of the method to identify the drought indices that can better characterize the drought condition in a specific location.

5.4. Comparison of Drought Indices through Characterizing the Hhistoric Drought Events

Comparison of the drought indices was carried out at three river gauging stations: Abbay at Bahirdar, Kessie, and Ethiopia-Sudan border (their locations are shown in Figure 1). These stations represent the upper, middle, and lower parts of the UBN basin, respectively. The time series values of the SRI-3 were calculated using the river flow data measured at these three gauging stations, whereas the areal average values were considered for the SPI-3, SPEI-3, ETDI, SMDI, and ADI. The result that was obtained for the upper part of the basin is presented for further discussion. Table 8 shows for the historic drought events, the characteristics resulting from each of the six indices. Based on the severity of the drought, the majority of the drought indices indicated 2003–2004 and 1983–1984 as the most severe drought years in the basin, except the SPEI, which showed 1989–1994 as the most severe drought year. The result further reveals that persistent droughts were observed in 2003–2004 and 1983–1984. Based on the mean (M) and maximum (Mmax) intensities, the SRI indicated the severity of 2003–2004 and 1984–1985 to have been higher than the other drought indices, except the SPEI (Mmax = −2.63). The SRI indicated the persistence of the 1983–1984 drought, which extended beyond the start and end year as provided by EM-DAT. The SPI ranks 1973–1974 as the most severe drought year based on mean intensity (M = −1.21) and maximum intensity (Mmax = −2.57). However, the drought durations of the SPI and SPEI for the same year are relatively small, most likely resulting from the earlier start of the drought in 1972. The onset month of each drought year is used as the other comparison criteria, and the result is shown in Table 8. According to EM-DAT, the 1973–1978 drought started in December 1973 and lasted until 1978. The ETDI showed the closest onset (November 1973) for this year. The SMDI showed a slightly earlier onset (July 1973) of the 1973–1978 drought, whereas the SPI and SPEI (January 1973) and SRI (March 1973) indicated an early start of this drought. The ADI showed the latest onset (July 1976) of this particular drought event. According to EM-DAT, the 1983–1984 drought started in May 1983, but all indices except the ETDI and ADI showed an early start (January 1983) of this drought event. The ETDI and ADI show a late start (December 1983 and September 1983, respectively) of the 1983–1984 drought. The 1989–1994 drought started in October 1989; the SPI and SPEI showed an August 1989 start, the SRI showed a December 1989, and the ETDI and SMDI indicated drought earlier than all the other indices. The drought start month indicated by ADI is December 1993. While this does not correspond with EM-DAT, it corresponds to the historic drought 1993–1994 reported in other studies (see Section 2.2.1). Note that it is likely that the 1993–1994 has been included in the EM-DAT database as part of the prolonged drought period from 1989 to 1994. For 2003–2004, there is no reference onset from EM-DAT to compare with the six drought indices. However, all of the indices indicated the onset of the drought in January and February 2003, except the SRI, which indicated the onset of drought in July 2003. Overall, the meteorological drought indices (SPI and SPEI) indicate the start date of four or more months before the reported start month, whereas agriculture drought indices onset often lagged by three or more months than the reported starting date of a historic drought. The aggregate drought index (ADI) most often lagged by some months. It can be noted from this result that no single drought index consistently indicates the exact onset of the drought. Similar results are obtained for the middle, and lower parts of the UBN basin.
Figure 6 compares the time series patterns of the six drought indices for the historic drought years 1973–1978, 1983–1984, 1989–1994, and 2003–2004. In the 1973–1978 drought (Figure 6a), the majority of the drought indices showed the occurrence of moderate to severe drought condition in several months between 1973 and 1975 (index value <−1.0 and lower), whereas mild drought and wet conditions were observed from 1976 to 1978 except for the ADI, which showed severe drought for ten months. Relatively persistent drought was indicated by the SRI and ETDI in 1973–1974, whereas other indices show a frequent jump between drought and non-drought conditions. In this study, persistent drought is defined if consecutive months (>10 months) showed below average values (zero). The majority of drought indices indicated mild to severe drought condition for several months (>6) in the year 1983/1984 (Figure 6b). The SRI showed the severity and persistency of the 1983–1984 drought and indicated drought conditions throughout 1983/1984. The SPI also showed mild to severe drought conditions for several months, except for non-drought conditions in mid-1984. The 1989–1994 drought is not well defined except for 1990, 1991, and 1994; the majority of the indices showed the occurrence of moderate to severe drought condition during these years. Like the 1983/1984 drought, the majority of the drought indices showed mild to extreme in 2003–2004, although the SRI, SPEI, and ADI showed non-drought conditions for some months. SRI showed severe to extreme drought condition, whereas SMDI showed mild to moderate drought conditions for several months in 1983–1984. SPI-3, SPEI-3, ETDI, and ADI showed out of drought condition for the period of 4 to 10 months. However, SPI-3, SMDI, and ETDI indicated the drought condition in the year 2003–2004, whereas SRI showed an opposite pattern than other indices. In general, SPI and SPEI showed similar patterns in most cases while ADI showed a frequent jump between drought and non-drought conditions. Moreover, ADI showed relatively higher positive and negative values perhaps because of a direct consideration of the eigenvector as weights/coefficients in combining the input variables. Other approach, such as the percentage of contribution, could be considered for future studies. SRI showed the persistency of the drought and non-drought condition for several months possibly because of less fluctuation of the river flow with time as compared to the rainfall.

6. Conclusions

In this study, we compared six drought indices—SPI and SPEI (meteorological indices), ETDI and SMDI (agricultural indices), SRI (hydrological index), and ADI (aggregate index)—to investigate how well characteristics of historic drought events in the Upper Blue Nile basin are identified by these indicators. Observed data were used for the precipitation, temperature, and streamflow inputs. Soil moisture and actual evapotranspiration data were estimated using the SWAT hydrological model.
In general, meteorological drought indices SPI-3 and SPEI-3 show a higher correlation with the agricultural drought indices SMDI and ETDI than with the hydrological drought index (SRI), whereas hydrological drought index (SRI) correlates better with meteorological drought indices at a 12-month aggregate period (SPI-12 and SPEI-12). This indicates that there exist some interconnection between drought indices in the sense that one drought index can explain more than one specific drought category to a certain degree. Moreover, each index has the potential to characterize and explain at least one recorded (historical) drought condition. It seems that a combination of more than one drought indices, appropriately selected for the specific region, is usually required for drought monitoring.
When comparing the drought onset dates indicated by the six indices for the historical (recorded) droughts, the meteorological drought indices (SPI and SPEI) showed early onsets when compared to the other drought indices except the 2003–2004 drought where SPEI showed a late onset. The agricultural (ETDI, SMDI, and ADI) and hydrological (SRI) drought indices showed late onsets, particularly the ADI, which lagged by several months for all events, except 2003–2004. When the onset dates indicated by the six drought indices are compared with the EM-DAT onset dates, meteorological drought indices (SPI-3, and SPEI-3) showed earlier onsets except for 2003–2004 drought where EM-DAT onset was unavailable. Similarly, the agricultural (ETDI and SMDI) and hydrological (SRI-3) drought indices showed earlier onsets of drought for two drought events and a late onset for one drought event. In contrast, ADI showed late onsets for two drought events and an early onset for 1983–1984. The comparison showed that none of the six drought indices could individually identify the onset of the four selected historic events. They could identify all events if all six indices (including different aggregation periods) were combined. Note that the ADI index tested in this study considers only three input variables at a 1-month aggregate period. The years 2003/2004 and 1983/1984 were indicated as the most severe drought years in the basin. An effective drought monitoring and planning is essential to mitigate drought impacts in the basin. For future drought monitoring for the Upper Blue Nile basin, developing a method that makes optimal use of multiple drought indices is recommended.

Supplementary Materials

The following are available online at www.mdpi.com/2076-3263/8/3/81/s1, Table S1: The Pearson’s correlation coefficient matrix computed between the paired values of the drought indices at multiple time scales for the selected stations.

Acknowledgments

The lead author was financially supported through the PhD fellowship of the Netherlands Fellowship Program (NFP). We are grateful to the National Meteorological Agency (NMA) and the Ministry of Water, Irrigation and Electricity of Ethiopia for providing the long-term climate and river flow data.

Author Contributions

All the authors were involved on framing the research. The lead author carried out analysis and prepared the first draft of the manuscript. All the authors contributed to interpretation of the results and preparation of the manuscript. Further all authors reviewed and provided comments/inputs to the manuscript at several stages.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wilhite, D.A. Drought as a natural hazard: Concepts and definitions. In Drought: A Global Assessment. Wilhite; Routledge: London, UK, 2000; Volume 1, pp. 3–18. [Google Scholar]
  2. Wilhite, A.; Svoboda, D.; Hayes, J. Understanding the complex impacts of drought: A key to enhancing drought mitigation and preparedness. Water Resour. Manag. 2007, 21, 763–774. [Google Scholar] [CrossRef]
  3. Sheffield, J.; Wood, E.F. Drought: Past Problems and Future Scenarios; Earthscan: London, UK, 2012. [Google Scholar]
  4. Tagel, G.; Van Der Veen, A.; Maathuis, B. Spatial and temporal assessment of drought in the northern highlands of Ethiopia. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 309–321. [Google Scholar]
  5. Bayissa, Y.; Semu, A.; Yunqing, X.; Schalk, A.; Shreedhar, M.; Dimitri, S.; Griensven, A.; Tadesse, T. Spatio-temporal assessment of meteorological drought under the influence of varying record length: The case of Upper Blue Nile Basin, Ethiopia. Hydrol. Sci. J. 2015, 60, 1927–1942. [Google Scholar] [CrossRef]
  6. Wilhite, D.; Buchanan-Smith, M. Drought as hazard: Understanding the natural and social context. In Drought and Water Crises: Science, Technology, and Management Issues; Wilhite, D.A., Ed.; CRC Press: Boca Raton, FL, USA, 2005; pp. 3–29. [Google Scholar]
  7. Valliyodan, B.T.; Nguyen, H.T. Understanding regulatory networks and engineering for enhanced drought tolerance in plants. Curr. Opin. Plant Biol. 2006, 9, 189–195. [Google Scholar] [CrossRef] [PubMed]
  8. Mishra, A.K.; Singh, V.P. A review of drought concepts. J. Hydrol. 2010, 391, 202–216. [Google Scholar] [CrossRef]
  9. Narasimhan, B.; Srinivasan, R. Development and evaluation of Soil Moisture Deficit Index (SMDI) and Evapotranspiration Deficit Index (ETDI) for agricultural drought monitoring. Agric. For. Meteorol. 2005, 133, 69–88. [Google Scholar] [CrossRef]
  10. Kogan, F.N. Global drought watch from space. Bull. Am. Meteorol. Soc. 1997, 78, 621–636. [Google Scholar] [CrossRef]
  11. Palmer, W.C. Meteorological Drought, 1st ed.; U.S. Weather Bureau: Washington, DC, USA, 1965.
  12. Guttman, N.B. Comparing the Palmer Drought Index and the Standardized Precipitation Index. J. Am. Water Resour. Assoc. 1998, 34, 113–121. [Google Scholar] [CrossRef]
  13. McKee, T.B.; Doesken, N.J.; Kleist, J. The relationship of drought frequency and duration to time scales. In Proceedings of the 8th Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993; pp. 179–184. [Google Scholar]
  14. McKee, T.B.; Doesken, N.J.; Kleist, J. Drought monitoring with multiple time scales. In Proceedings of the Ninth Conference on Applied Climatology, Dallas, TX, USA, 15–20 January 1995; pp. 233–236. [Google Scholar]
  15. Wu, H.; Hayes, M.J.; Weiss, A.; Hu, Q. An evaluation the standardized precipitation index, the China-Z index and the statistical Z-score. Int. J. Climatol. 2001, 21, 745–758. [Google Scholar] [CrossRef]
  16. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A multiscalar drought index sensitive to global warming: The standardized precipitation evapotranspiration index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef]
  17. Shukla, S.; Wood, A.W. Use of a standardized runoff index for characterizing hydrologic drought. Geophys. Res. Lett. 2008, 35, 1–7. [Google Scholar] [CrossRef]
  18. Shafer, B.A.; Dezman, L.E. Development of a surface water supply index (SWSI) to assess the severity of drought conditions in snowpack runoff areas. In Proceedings of the Western Snow Conference, Colorado State University, Fort Collins, CO, USA, 19–23 April 1982; pp. 164–175. [Google Scholar]
  19. Gibbs, W.J.; Maher, J.V. Rainfall Deciles as Drought Indicators; Bureau of Meteorology: Melbourne, Australia, 1967; p. 29.
  20. Brown, J.F.; Wardlow, B.D.; Tadesse, T.; Hayes, M.J.; Reed, B.C. The vegetation drought response index (VegDRI): A new integrated approach for monitoring drought stress in vegetation. GISci. Remote Sens. 2008, 45, 16–46. [Google Scholar] [CrossRef]
  21. Zargar, A.; Sadiq, R.; Naser, B.; Khan, F.I. A review of drought indices. Environ. Rev. 2011, 19, 333–349. [Google Scholar] [CrossRef]
  22. Keyantash, J.A.; Dracup, J.A. An aggregate drought index: Assessing drought severity based on fluctuations in the hydrologic cycle and surface water storage. Water Resour. Res. 2004, 40. [Google Scholar] [CrossRef]
  23. Ntale, H.K.; Gan, T. Drought Indices and their Application to East Africa. Int. J. Climatol. 2003, 23, 1335–1357. [Google Scholar] [CrossRef]
  24. Morid, S.; Smakhtin, V.; Moghaddasi, M. Comparison of seven meteorological indices for drought monitoring in Iran. Int. J. Climatol. 2006, 26, 971–985. [Google Scholar] [CrossRef]
  25. Barua, S.; Perera, B.J.C.; Ng, A.W.M. A comparative drought assessment of Yarra River Catchment in Victoria, Australia. In Proceedings of the 18th World IMACS/MODSIM Congress, Cairns, Australia, 13–17 July 2009; pp. 13–17. [Google Scholar]
  26. Dogan, S.; Berktay, A.; Singh, V.P. Comparison of multi-monthly rainfall based drought severity indices, with application to semi-arid Konya closed basin Turkey. J. Hydrol. 2012, 470, 255–268. [Google Scholar] [CrossRef]
  27. Trambauer, P.; Maskey, S.; Werner, M.; Pappenberger, F.; Van Beek, L.; Uhlenbrook, S. Identification and simulation of space-time variability of past hydrological drought events in the Limpopo river basin, Southern Africa. Hydrol. Earth Syst. Sci. Discuss. 2014, 11, 2639–2677. [Google Scholar] [CrossRef]
  28. Jain, V.K.; Pandey, R.P.; Jain, M.K.; Byun, H.R. Comparison of drought indices for appraisal of drought characteristics in the Ken River Basin. Weather Clim. Extrem. 2015, 8, 1–11. [Google Scholar] [CrossRef]
  29. Okpara, J.N.; Tarhule, A. Evaluation of drought indices in the Niger Basin, West Africa. J. Geogr. Earth Sci. 2015, 3, 1–32. [Google Scholar]
  30. Naumann, G.; Dutra, E.; Barbosa, P.; Pappenberger, F.; Wetterhall, F.; Vogt, J.V. Comparison of drought indicators derived from multiple data sets over Africa. Hydrol. Earth Syst. Sci. 2013, 18, 1625–1640. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, J.; Li, Y.; Ren, Y.; Liu, Y. Comparison among several drought indices in the Yellow River Valley. J. Nat. Resour. 2013, 28, 1337–1349. [Google Scholar]
  32. Houcine, A.; Bargaoui, Z. Comparison of rainfall based SPI drought indices with SMDI and ETDI indices derived from a soil water budget model. J. Geophys. Res. 2012, 14, EGU2012-2666. [Google Scholar]
  33. Heim, R.R. A review of Twentieth-Century drought indices used in the United States. Bull. Am. Meteorol. Soc. 2002, 83, 1149–1165. [Google Scholar] [CrossRef]
  34. Viste, E.; Korecha, D.; Sorteberg, A. Recent drought and precipitation tendencies in Ethiopia. Theor. Appl. Climatol. 2013, 112, 535–551. [Google Scholar] [CrossRef] [Green Version]
  35. Edossa, D.; Babel, S.; Gupta, D. Drought analysis in the Awash River Basin, Ethiopia. Water Resour. Manag. 2009, 24, 1441–1460. [Google Scholar] [CrossRef]
  36. Dessie, M.D.; Verhoest, N.; Pauwels, V.; Negatu, T.A.; Poesen, J.; Adgo, E.; Deckers, J.; Nyssen, J. Analyzing runoff processes through conceptual hydrological modeling in the Upper Blue Nile Basin, Ethiopia. Hydrol. Earth Syst. Sci. 2014, 18, 5149–5167. [Google Scholar] [CrossRef] [Green Version]
  37. Conway, D. The climate and hydrology of the Upper Blue Nile River. Geograph. J. 2000, 166, 49–62. [Google Scholar] [CrossRef] [Green Version]
  38. Tekleab, S.; Mohamed, Y.; Uhlenbrook, S. Hydro-Climatic trends in the Abay/Upper Blue Nile basin, Ethiopia. Phys. Chem. Earth 2013, 61, 32–42. [Google Scholar] [CrossRef]
  39. Kebede, S.; Travi, Y.; Alemayehu, T.; Marc, V. Water balance of Lake Tana and its sensitivity to fluctuations in rainfall, Blue Nile basin, Ethiopia. J. Hydrol. 2006, 316, 233–247. [Google Scholar] [CrossRef]
  40. Characterization and Atlas of the Blue Nile Basin and Its Sub Basins. Available online: http://publications.iwmi.org/pdf/H042502.pdf (accessed on 24 February 2018).
  41. Mellander, P.; Gebrehiwot, S.; Gärdenäs, A.; Bewket, W.; Bishop, K. Summer rains and dry seasons in the Upper Blue Nile Basin: the predictability of half a century of past and future spatiotemporal patterns. PLoS ONE 2013, 8, e68461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Gebremicael, T.G.; Mohamed, Y.A.; Betrie, G.D.; van der Zaag, P.; Teferi, E. Trend analysis of runoff and sediment fluxes in the Upper Blue Nile basin: A combined analysis of statistical tests, physically-based models and landuse maps. J. Hydrol. 2013, 482, 57–68. [Google Scholar] [CrossRef]
  43. Betrie, G.D.; Mohamed, Y.A.; van Griensven, A.; Srinivasan, R. Sediment management modelling in the Blue Nile Basin using SWAT model. Hydrol. Earth Syst. Sci. 2011, 15, 807–818. [Google Scholar] [CrossRef]
  44. Annual Disaster Statistical Review 2010: The Numbers and Trends. Available online: http://www.cred.be/downloadFile.php?file=sites/default/files/ADSR_2010.pdf (accessed on 24 February 2018).
  45. Guha-Sapir, D.; Below, R. The Quality and Accuracy of Disaster Data: A Comparative Analyse of 3 Global Data Sets; Working Paper ID 191; World Bank: Washington, DC, USA, 2002. [Google Scholar]
  46. Hargreaves, G.H.; Samani, Z.A. Estimating potential evapotranspiration. J. Irrig. Drain. Div. 1982, 108, 225–230. [Google Scholar]
  47. MARD. Land Use/Cover Classification of Ethiopia; Ministry of Agriculture and Rural Development: Addis Ababa, Ethiopia, 2004.
  48. Bceom French Engineer Consultant. Abay River Basin Integrated Master Plan, Main Report; Ministry of Water Resources (MOWR): Addis Ababa, Ethiopia, 1999.
  49. Mengistu, D.; Sorteberg, A. Validation of SWAT simulated streamflow in the Eastern Nile and sensitivity to climate change. Hydrol. Earth Syst. Sci. Discuss. 2011, 8, 9005–9062. [Google Scholar] [CrossRef]
  50. Setegn, S.G.; Srinivasan, R.; Dargahi, B. Hydrological modelling in the Lake Tana Basin, Ethiopia using SWAT model. Open Hydrol. J. 2008, 2, 49–62. [Google Scholar] [CrossRef]
  51. Easton, Z.M.; Fuka, D.R.; White, E.D.; Collick, A.S.; Biruk, A.B.; McCartney, M.; Awulachew, S.B.; Ahmed, A.A.; Steenhuis, T.S. A multi basin SWAT model analysis of runoff and sedimentation in the Blue Nile, Ethiopia. Hydrol. Earth Syst. Sci. 2010, 14, 1827–1841. [Google Scholar] [CrossRef]
  52. Van Griensven, A.; Ndomba, P.; Yalew, S.; Kilonzo, F. Critical review of SWAT applications in the upper Nile basin countries. Hydrol. Earth Syst. Sci. 2012, 16, 371–3381. [Google Scholar] [CrossRef]
  53. Steenhuis, T.S.; Collick, A.S.; Easton, Z.M.; Leggesse, E.S.; Bayabil, H.K.; White, E.D.; Awulachew, S.B.; Adgo, E.; Ahmed, A.A. Predicting discharge and sediment for the Abay (Blue Nile) with a simple model. Hydrol. Process. 2009, 23, 3728–3737. [Google Scholar]
  54. Khan, S.; Gabriel, H.; Rana, T. Standard precipitation index to track drought and assess impact of rainfall on watertables in irrigation areas. Irrig. Drain. Syst. 2008, 22, 159–177. [Google Scholar] [CrossRef]
  55. Moreira, E.E.; Coelho, C.A.; Paulo, A.A.; Pereira, L.S.; Mexia, J.T. SPI-based drought category prediction using loglinear models. J. Hydrol. 2008, 354, 116–130. [Google Scholar] [CrossRef]
  56. Kalcic, M.M.; Chaubey, I.; Frankenberger, J. Defining soil and water assessment tool (SWAT) hydrologic response units (HRUs) by field boundaries. Int. J. Agric. Biol. Eng. 2015, 8, 69–80. [Google Scholar]
Figure 1. Upper Blue Nile basin with river network, gauging stations, and the sub-basin discretisation used by the Soil and Water Assessment Tool (SWAT) model to estimate soil moisture and actual evapotranspiration (Section 2.2.2), upstream and downstream of the Lake Tana outlet.
Figure 1. Upper Blue Nile basin with river network, gauging stations, and the sub-basin discretisation used by the Soil and Water Assessment Tool (SWAT) model to estimate soil moisture and actual evapotranspiration (Section 2.2.2), upstream and downstream of the Lake Tana outlet.
Geosciences 08 00081 g001
Figure 2. The time series plots of the drought indices for Gondar (a) and Debremarkos (b) stations. The plots for Standardized Precipitation Index (SPI), Standardized Precipitation Evapotranspiration Index (SPEI), and Standardized Runoff Index (SRI) at 6-, 9-, and 12-month time scales are shown in the single panel since it is easier to visualize the trend. Similar approach was considered for Evapotranspiration Deficit Index (ETDI), Soil Moisture Deficit Index (SMDI), and Aggregate Drought Index (ADI).
Figure 2. The time series plots of the drought indices for Gondar (a) and Debremarkos (b) stations. The plots for Standardized Precipitation Index (SPI), Standardized Precipitation Evapotranspiration Index (SPEI), and Standardized Runoff Index (SRI) at 6-, 9-, and 12-month time scales are shown in the single panel since it is easier to visualize the trend. Similar approach was considered for Evapotranspiration Deficit Index (ETDI), Soil Moisture Deficit Index (SMDI), and Aggregate Drought Index (ADI).
Geosciences 08 00081 g002aGeosciences 08 00081 g002b
Figure 3. Comparison of SPI, SPEI, and SRI with ETDI, SMDI, and ADI at 1-month (a), 3-month (b), 6-month (c), 9-month (d), and 12-month (e) time scales based on percentage of drought months at the locations of the selected meteorological stations. A line is used to connect values of the same variable at distinct locations to increase the visibility of the points in this figure. .
Figure 3. Comparison of SPI, SPEI, and SRI with ETDI, SMDI, and ADI at 1-month (a), 3-month (b), 6-month (c), 9-month (d), and 12-month (e) time scales based on percentage of drought months at the locations of the selected meteorological stations. A line is used to connect values of the same variable at distinct locations to increase the visibility of the points in this figure. .
Geosciences 08 00081 g003aGeosciences 08 00081 g003b
Figure 4. Percentage of severe (a), and extreme (b) droughts months observed for the drought indices within the study period (1970–2010) for the selected meteorological stations.
Figure 4. Percentage of severe (a), and extreme (b) droughts months observed for the drought indices within the study period (1970–2010) for the selected meteorological stations.
Geosciences 08 00081 g004
Figure 5. The maximum duration of consecutive drought indicated by SPI (a), SPEI (b), SRI (c) and other indices (d) such as ETDI, SMDI, and ADI. Five temporal scales (i.e., 1-, 3-, 6-, 9-, and 12-months) were considered for SPI, SPEI, and SRI.
Figure 5. The maximum duration of consecutive drought indicated by SPI (a), SPEI (b), SRI (c) and other indices (d) such as ETDI, SMDI, and ADI. Five temporal scales (i.e., 1-, 3-, 6-, 9-, and 12-months) were considered for SPI, SPEI, and SRI.
Geosciences 08 00081 g005
Figure 6. The time series plot of the six drought indices for 1973/1978 (a), 1983/1984 (b), 1989/1994 (c), and 2003/2004 (d) historic drought years in the UBN basin.
Figure 6. The time series plot of the six drought indices for 1973/1978 (a), 1983/1984 (b), 1989/1994 (c), and 2003/2004 (d) historic drought years in the UBN basin.
Geosciences 08 00081 g006aGeosciences 08 00081 g006b
Table 1. The start and end date of the historic drought events and the total number of people affected during the years 1970–2010. The provinces indicated in bold are situated (partly) within the UBN basin.
Table 1. The start and end date of the historic drought events and the total number of people affected during the years 1970–2010. The provinces indicated in bold are situated (partly) within the UBN basin.
Start DateEnd DateLocationTotal Affected Population (106)Total Number of Population in Ethiopia (106)Ratio b/n Column 4 and 5 (%)
December 19731978Wollo, North Shoa, Tigray, Kangra province332.579
May 19831984Wollo, Gondar, Gore, Tigray, Shoa, Harerge, Sidamo7.7535.2422
June 19871987Ogaden, Tigray, Wollo, Shewa, Gamo Gofa, Sidama, Gondar, Bale748.0615
October 19891994Northern Ethiopia, Tigray, Wollo, Gondar, Harerge6.548.0614
20032004Tigray, Oromia, Amhara, Somali, Afar province12.676.6116
May 2008October 2009Oromia, Somali, Amhara, Afar, Tigray, SNNPR province6.487.567
Source: EM-DAT International Disaster Database. Centre for Research on the Epidemiology of Disasters-CRED; http://www.emdat.be/database (last access: 17 July 2016); http://www.worldometers.info/world-population/ethiopia-population/. The provinces indicated in bold are situated (partly) within the UBN basin.
Table 2. Summary of the input data used in this study.
Table 2. Summary of the input data used in this study.
VariablesData LengthResolutionsSource
TemporalSpatial
Actual Evapotranspiration1970–2010MonthlyHRUSWAT based
soil moisture data1970–2010MonthlyHRUSWAT based
Rainfall1970–2010DailyPointNMA Ethiopia
Temperature1970–2010DailyPointNMA Ethiopia
Relative Humidity1970–2010DailyPointNMA Ethiopia
Solar Radiation1970–2010DailyPointNMA Ethiopia
Wind Speed1970–2010DailyPointNMA Ethiopia
River discharge1970–2010Daily/monthlyPointNMA Ethiopia
DEM 90 mSRTM
Land Use2000 5 kmMARD
Table 3. SPI values that show different categories of drought severity (McKee et al. [13]).
Table 3. SPI values that show different categories of drought severity (McKee et al. [13]).
SPI ValueDrought Category
−2.00 and lessExtreme
−1.50 to −1.99Severe
−1.00 to −1.49Moderate
0 to −0.99Near normal or mild
Above 0No drought
Table 4. Aggregate drought index (ADI) drought category classification.
Table 4. Aggregate drought index (ADI) drought category classification.
above 0.92Wet
−0.95 to 0.92Near normal
−1.40 to −0.96Moderate drought
−1.69 to −1.41Severe drought
−1.70 or lessExtreme drought
Table 5. Summary of the drought indices used in this study.
Table 5. Summary of the drought indices used in this study.
Drought IndexVariableDistribution Fitted
SPIPrecipitationGamma
SPEIPrecipitation, PETLog-logistic
SRIRunoffCDF standardized to Gaussian values
ETDIActual evapotranspiration
SMDISoil moisture
ADIPrecipitation, actual evapotranspiration, and soil moisture
Table 6. The Pearson’s correlation coefficient matrix computed between the paired values of the drought indices at multiple time scales for (a) Gondar station and (b) Debremarkos station. The SRI used in the cases are calculated using flow records at Bahirdar station and Debremarkos station respectively.
Table 6. The Pearson’s correlation coefficient matrix computed between the paired values of the drought indices at multiple time scales for (a) Gondar station and (b) Debremarkos station. The SRI used in the cases are calculated using flow records at Bahirdar station and Debremarkos station respectively.
(a)
SPI-11
SPI-30.651
SPI-60.530.771
SPI-90.480.730.771
SPI-120.450.680.720.961
SPEI-10.850.550.360.280.261
SPEI-30.490.890.610.480.450.631
SPEI-60.350.590.930.720.650.420.691
SPEI-90.320.530.750.940.830.330.550.791
SPEI-120.300.510.690.830.940.310.520.720.891
SRI-10.200.280.330.400.410.160.230.160.210.251
SRI-30.220.310.310.440.480.110.260.250.260.300.811
SRI-60.200.240.370.490.520.070.190.310.310.330.700.881
SRI-90.240.320.430.510.530.080.160.260.330.340.630.800.931
SRI-120.210.310.410.520.550.040.150.240.300.360.570.740.870.951
ETDI0.530.670.590.430.450.590.490.340.280.280.230.210.220.280.271
SMDI0.440.650.450.380.290.430.550.370.320.320.250.280.210.180.170.681
ADI0.750.590.400.320.300.780.590.390.380.360.330.230.150.120.130.680.631
(b)
SPI-11
SPI-30.551
SPI-60.350.651
SPI-90.270.480.741
SPI-120.200.390.610.811
SPEI-10.870.540.330.260.201
SPEI-30.480.890.590.430.340.591
SPEI-60.300.580.910.660.530.360.661
SPEI-90.240.440.690.920.730.290.490.741
SPEI-120.180.350.560.750.910.220.390.600.811
SRI-10.100.180.170.140.160.070.150.150.080.071
SRI-30.100.190.230.170.180.050.140.200.120.080.861
SRI-60.080.200.280.260.230.050.130.210.200.130.630.831
SRI-90.060.190.290.320.300.050.140.220.240.210.490.670.871
SRI-120.050.160.270.320.360.030.120.210.230.240.480.560.710.841
ETDI0.270.510.390.250.220.290.490.370.270.260.010.000.030.040.091
SMDI0.240.420.350.300.280.230.370.280.230.220.140.130.120.130.080.381
ADI0.070.100.080.060.030.090.060.050.010.030.160.140.090.020.080.040.001
Table 7. The summary of the significant test for the selected 16 stations.
Table 7. The summary of the significant test for the selected 16 stations.
SPI-1SPI-3SPI-6SPI-9SPI-12SPEI-1SPEI-3SPEI-6SPEI-9SPEI-12SRI-1SRI-3SRI-6SRI-9SRI-12ETDISMDI
SPI-3
SPI-6
SPI-9
SPI-12
SPEI-1
SPEI-3
SPEI-6
SPEI-9
SPEI-12
SRI-1
SRI-3
SRI-6
SRI-9
SRI-12
ETDI
SMDI
ADI
Geosciences 08 00081 i001 p < 0.01, Geosciences 08 00081 i002 p < 0.1, Geosciences 08 00081 i003 p > 0.1.
Table 8. Characteristics of the historic drought events as identified by the six drought indices for the Upper Blue Nile basin upstream of Bahirdar. Note that if the index indicated the onset to be before January of the reported calendar year of the historic drought, this is presented in the table with “before” and if, according to the index, a drought extended beyond December of the reported calendar year, this is presented with “after”. The drought duration in these cases is based on a start in January and an end in December of the reported drought year.
Table 8. Characteristics of the historic drought events as identified by the six drought indices for the Upper Blue Nile basin upstream of Bahirdar. Note that if the index indicated the onset to be before January of the reported calendar year of the historic drought, this is presented in the table with “before” and if, according to the index, a drought extended beyond December of the reported calendar year, this is presented with “after”. The drought duration in these cases is based on a start in January and an end in December of the reported drought year.
Drought IndexSPI-3SPEI-3SRI-3ETDISMDIADI
1973–1978 drought (Onset on December according to EM-DAT)
Starting dateBefore 01/1973Before 01/197303/197311/197307/197307/1976
Ending date04/197304/197306/197406/197405/197404/1977
Mean Intensity, M−1.21−1.23−0.75−0.48−0.42−0.94
Maximum intensity, Mmax−2.57−2.33−1.36−1.47−1.19−2.43
Duration, D (months)441681110
Severity, S−4.84−4.91−11.94−3.85−4.57−9.38
1983–1984 drought (Onset on May according to EM-DAT)
Starting dateBefore 01/1983Before 01/1983Before 01/198312/1983Before 01/198309/1983
Ending date05/198407/1983After 12/1984After 12/1984After 12/198404/1984
Mean Intensity, M−0.99−0.86−1.36−-0.66−0.53−0.46
Maximum intensity, Mmax−2.45−1.95−2.53−1.53−1.78−1.98
Duration, D (months)1772413248
Severity, S−16.82−6.02−32.59−8.62−12.79−3.68
1989–1994 drought (Onset on October according to EM-DAT)
Starting date06/198907/198912/198904/199002/199012/1993
Ending date08/199008/199012/199002/199101/1991After 12/1994
Mean Intensity, M−0.69−0.65−0.70−0.85−0.46−0.57
Maximum intensity, Mmax−2.49−2.09−1.38−1.92−1.18−1.57
Duration, D (months)15141391113
Severity, S−10.30−9.13−-9.06−7.62−5.08−7.38
2003–2004 drought (Onset not defined according to EM-DAT)
Starting dateBefore 01/200302/200407/2004Before 01/2003Before 01/2003Before 01/2003
Ending date10/200412/200412/2004After 12/2004After 12/200412/2003
Mean Intensity, M−0.59−0.71−1.61−0.51−0.45−0.92
Maximum intensity, Mmax−1.95−2.63−2.13−1.23−1.36−2.27
Duration, D (months)22116242412
Severity, S−12.95−7.79−9.64−12.26−10.85−11.09

Share and Cite

MDPI and ACS Style

Bayissa, Y.; Maskey, S.; Tadesse, T.; Van Andel, S.J.; Moges, S.; Van Griensven, A.; Solomatine, D. Comparison of the Performance of Six Drought Indices in Characterizing Historical Drought for the Upper Blue Nile Basin, Ethiopia. Geosciences 2018, 8, 81. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences8030081

AMA Style

Bayissa Y, Maskey S, Tadesse T, Van Andel SJ, Moges S, Van Griensven A, Solomatine D. Comparison of the Performance of Six Drought Indices in Characterizing Historical Drought for the Upper Blue Nile Basin, Ethiopia. Geosciences. 2018; 8(3):81. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences8030081

Chicago/Turabian Style

Bayissa, Yared, Shreedhar Maskey, Tsegaye Tadesse, Schalk Jan Van Andel, Semu Moges, Ann Van Griensven, and Dimitri Solomatine. 2018. "Comparison of the Performance of Six Drought Indices in Characterizing Historical Drought for the Upper Blue Nile Basin, Ethiopia" Geosciences 8, no. 3: 81. https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences8030081

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop