Hydro-meteorological drought in Addis Ababa: A characterization study

.


Introduction
Drought is a common natural disaster linked with scarcity of water due to reduction in the amount of rainfall over an extended period of time.Drought occurrences are expected to become more severe and frequent in future years [1].Rainfall is the primary driver of any drought.Drought leads to a reduction of surface water (reservoir storage level).Surface and subsurface water resources are commonly affected by an extended period of drought and this results from rainfall variability [2,3].Meteorological droughts are impermanent, recurring natural disasters, which occurs from the shortage of rainfall and can bring significant economic losses.There are numerous indicators and hypothesis that are being employed for drought monitoring [4,5].Based on drought indices, various categories can be identified-permanent, seasonal, contingent, and invisible drought.Drought maybe meteorological, agricultural or hydrological especially in its manifestation and impact [6].
As a country, Ethiopia is under rapid development both in rural and urban settings; there are enormous constructions in all sectors, rapid urbanization, industrialization, floriculture and big farming are being practiced.These anthropogenic activities along with climate change are affecting the water resources of the country and it also exacerbate the urban water supply shortage.The economic impact of severe drought is that people spend more money to access water for household and industrial purpose, since reservoir level and even ground water levels have decreased.Drought also affects plants, animals, and human health.Hence, many drought indices are proposed to examine the historical drought trends based on rainfall and temperature data which supports water security in the country.Feyissa et al. [7] studied downscaling of future temperature and precipitation extremes in Addis Ababa under climate change.The study result shows that the maximum temperature increases were in the range of 0.9 °C (RCP4.5) in 2020 to 2.1 °C (CGCM3A2) in 2080 at Addis Ababa Observatory.In Ethiopia, rainfall is very variable both in amount and distribution through regions and seasons [8].The frequency and ardency of extremes are likely impacted by climate changes and therefore the changes in climate extremes can affect other forms of biological and physical systems and their impacts can even lead to global warming [9].In Addis Ababa-the capital of Ethiopia, the temperature is usually affected by human activities together with global climate change [10].The rise in temperature will aggravate the urban heat highland effects in warm seasons and an increase in rainfall is expected along with a possible risk of flooding [7].The rising trend of temperature owing to the adverse effects of climate change and other factors can cause climate extremes within the city of Addis Ababa [11].Hence, the drought in Addis Ababa can be attributed to the fluctuations in the climate pattern of the city.
The prediction of drought can play a vital role in the mitigation of its effects.In practice, drought indices are usually employed in the assessment of drought events.In recent times, drought indices have supported meteorological and hydrological stations easily to predict drought events using the relevant climatic data on the basis of practical applications.Drought indices has become the foremost option for drought monitoring and characterization [12].In the literature, many studies have used the two common drought indices-the Standardized Precipitation Index (SPI) and the Reconnaissance Drought Index (RDI) for monitoring and forecasting of droughts in numerous areas.The RDI was used to determine the temporal and spatial analysis of meteorological drought in the Awash River Basin (Ethiopia) [13].Likewise, many research studies has also proposed the RDI as one climatic index for the detection of potential climatic changes [14].However, drought indices are typically continuous functions of temperature and/or rainfall, river discharge, or additional assessable hydro-meteorological variables.Consequently, the SPI is also used as one of the large numbers of drought indices [15].SPIs were utilized in drought analysis based on monthly rainfall data to examine the drought characteristics of Meric-Ergene, Gediz, Seyhan, and Ceyhan hydrological basins in Turkey [16].Pashiardis and Michaelides [17] examined regional drought assessment by calculating the SPI and RDI indices which are based on past monthly rainfall and temperature data.Both indices are similar and can be used efficiently to analyze and monitor drought situations and events.Similarly, in Iran, another study applied the SPI and RDI indices for different time scales in forty meteorological stations.The study found that drought severity was obtainable in the drought year from 1999 to 2000 for the two indices.It also showed that the two indices were more significant in the 3-months, 6-months, and 9-months than longer time scales [18].
There are other numerous examples of SPI and RDI applications such as Ansarifard and Shamsnia [19] who used the SPI and RDI for the analysis of rainfall for drought monitoring.The SPI result showed near-normal condition, while the RDI categorized as moderately dry condition for the year 1999-2000.The SPI has several features that give it an edge over earlier indices, including its simplicity and temporal flexibility [20].Similarly, the SPI was used for the examination of drought intensity assessment in northeastern Nigeria and the decadal SPI result revealed a mild drought [6].Zarei et al. [21] employed the modified RDI for the assessment of drought severity in the south of Iran.The study result showed that the trends for annual drought were significantly increasing at 95% confidence level for the areas with normal, dry and extremely dry conditions.Likewise, the SPI was used to examine the rainfall deficit for multiple timescales in Tegal City, Central Java (Indonesia).From all scarcity periods, the result of the study revealed that the foremost severe drought within the study temporal horizon fell in 2015 [22].
We recognize that there are other studies on this subject matter, and thus, we have just presented a snapshot of a few.The foregoing indicates that drought assessment of a given jurisdiction involves the computation of the two (2) well-known drought indices namely; the RDI and SPI.The RDI is computed based on the historic monthly precipitation and temperature data.Whereas, the SPI is calculated only with rainfall data.Both indices can be employed effectively to analyze the occurrence of drought history and used to predict drought trends.However, RDI has more advantage to SPI since it includes temperature data.Accordingly, this study is focused on the precipitation and/or temperature variation in the study area over three stations-Sendafa, Bole, and Observatory in Addis Ababa, Ethiopia.The historic precipitation and temperature data used for the three stations are from 1983 to 2015.
The main objective of this paper is to examine the occurrence of historical drought trend in Addis Ababa city.It seeks to provide scientific information for the water sector professionals and support for better management of urban water resources.With respect to scholarly contributions, this paper adds to the literature in this domain by introducing one of the earliest case studies to examine historical drought trends in Ethiopia.In a useful manner, the study applies the two well-known drought indices for a jurisdiction that have scanty information on its drought trend.Moreover, beyond Ethiopia, the analytical frameworks presented here can be used by other researchers to examine historical drought trends in other locations of the world.Therefore, while this paper is purposed within a national perspective, we believe that it will be relevant for international audience as it also supports the realization of the UN Sustainable Development Goals (SDG) of clean water & sanitation (SDG 6), and climate action (SDG 13).
Section 2 of the paper outlines the study methodology employed in the paper/investigation which contains the study area description, quality control of the data, SPI and RDI etc. Section 3 presents the study results and in-depth discussion.Next is section 4 which states the study conclusions.

Materials and method
We used the Drought Indices Calculator (DrinC) for our study.DrinC is a software package which was established for providing easy, though flexible interface for the computation of drought indices [29].The Reconnaissance Drought Indices (RDI) and Standard Precipitation Index (SPI) were employed to examine the occurrence of drought history at a longer time scale.Figure 1 shows the simplified study methodological framework.

Location of the study area
Addis Ababa is geographically located at 38° 44' East and 9° 1' North longitude and latitude respectively, and has a mild climate.It have three layers of Government hierarchy namely; City Government, Sub City Administrations, and Kebele Administrations at the top, middle, and bottom levels respectively.The average elevation of the city is around 2,500 meters above sea level.The period from March to May, is the warmest of the year.Ethiopia has approximately 1,200 mm of rain fall throughout the year, with the maximum being from June through September.Figure 2 shows the study area.

Quality control of data
The data quality was checked visually and statistically.Visually, data cleaning was performed to check the data entry error.Whereas, statistically the quality control were carried out by using boxplot and SPSS software package.The tool is capable of identifying outliers and missed data.An outlier is defined as the value above the upper fence or below the lower fence.The upper and lower fence represent values more than 75 th percentile (3 rd quartiles) and 25 th percentile (1 st quartiles) respectively.The rainfall and temperature data were verified for outlier and missing data to avoid mistakes that can cause changes in the final result.Finally, the DrinC was employed for the analysis of drought indices and used to investigate the RDI and SPI indices [29].

Standard Precipitation Index (SPI)
The SPI is a widely used meteorological drought index and it was first introduced by [15].It is basically the transformation of the precipitation (rainfall) time series into a standardized normal distribution.The SPI is calculated as per Eq 1.

SPI =
Where xi is monthly rainfall,  ̅ is mean, and σ is the standard deviation of the data.SPI can be computed for different time scales.SPI is computed for short time (up to 3 month, medium term (3-12 months), and long term accumulation periods (12-24 months).The present study computed the SPI values for medium-term cumulated (3 to 12 months, that is SPI 3, and SPI 12) since it is more suitable for measuring the impact on stream flow and reservoir levels [15,28].But, the SPI values were also computed for short-term time scales (1-month SPI) to highlight the agricultural drought.Table 1 shows the SPI value used for drought classification.

Reconnaissance Drought Index (RDI)
The Reconnaissance Drought Index (RDI) is one among the meteorological drought indices used to assess drought severity.The RDI is expressed in three forms which are commonly called the initial value of Reconnaissance Drought Index (αk), normalized Reconnaissance Drought Index (RDIn), and standardized Reconnaissance Drought Index (RDIst) [34,29].The standardized expression of RDI is not significantly affected by the method of potential evapotranspiration (PET) calculation, whereas the initial and the normalized RDI are directly influenced by other methods.However, if precipitation and air temperature data are the only available data, then the simplified version of FAO Penman-Monteith (or Hargreaves) method, could be used [35].In cases where adequate long time data of several meteorological variables are obtainable, the researchers should use the widely accepted FAO Penman-Monteith method [36].
In this study, the potential evapotranspiration (PET) is computed based on rainfall and temperature data since solar radiation, saturation vapor pressure, relative humidity, soil heat flux density, wind speed, vapor pressure deficit data are not available.Thus, the potential evapotranspiration is calculated by the following Eq 2 [37].PET = 0.0023 ( Where T max is the maximum temperature, T min is the minimum temperature, Ra is the extraterrestrial solar radiation, which depends on the year as well as latitude of the observation site.The initial value of the Reconnaissance Drought Index (αk) is presented in an aggregated form employing a monthly time step and will be calculated for every month of the hydrological year (or an entire year).The αk will be computed by Eq 3.
Where: Pij and PETij are the precipitation (rainfall) and potential evapotranspiration (PET) of the j th month of the i th year and i usually starts from September as per Ethiopian hydrological time.Hence, September k = 1 and N is the overall entire number of years for the available data.The Normalized values of the Reconnaissance Drought Index (RDIn) is calculated using Eq 4.

RDI 𝑛
Where, α k () ̅̅̅̅̅ is the arithmetic mean of α  () values and calculated for the N years of given data.The values of αk follows satisfactorily both the lognormal as well as the gamma distributions in which they were tested [32].
The initial formulation of RDIst used the thought that α  values follow the log-normal (LN) distribution.So, the following equations are often used for the computation of RDIst: In which, yk is ln(αk (i) ), yk ̅̅̅ is the arithmetic mean of yk and std αyk is its standard deviation (SD).
However, the αk values follow both the ln and the gamma distribution values at the majority of several locations and time scales.But in most of the cases, the gamma distribution was proved to be more successful.Therefore, the calculation of RDIst might be performed better by fitting the gamma probability density function (pdf) to the given frequency distribution of ak.The positive values of RDIst indicate the wet periods, whereas the negative values indicate dry periods compared to the conditions of the studied area [38,34].The RDI value used for drought classification is stated in Table 2.

Assessment of rainfall variability
The given input parameters are derived via descriptive statistical information.The procedure for the descriptive statistical summary analysis can be displayed in the form of a table.The statistical summary of parameters are the sample size (the number of observations), minimum, maximum, mean, standard deviation (SD), variance, and coefficient of variance (CV).The Minimum value shows the lowest category while the Maximum value shows the highest category of the given variable.Hence, the mean is calculated as stated below (Eq.6).

Mean (𝑋) ̅̅̅
Here, Xi is the data value, and n is the sample size.The sample variance indicates the classical measure of spread.Similar to the mean, it is strongly affected by outliers.In a population, the standard deviation and variance can both be used to measure the variability of the given population.Variance simply refers to the average squared deviations from the arithmetic mean.Whereas, the SD refers to the square root of the variance.Consequently, the variance is simply the square of the SD (Eq.7).Similarly, assessment of rainfall variability were determine by the CV as stated below (Eq.9).

Data collection and processing
The daily climatic data for temperature and rainfall (precipitation) were taken from the National Meteorological Agency (NMA).The past (historic) precipitation and temperature data were collected from 1983 to 2015 (NY=33) for all stations (i.e.Sendafa, Bole, and Observatory station).The whole research data for this study were collected based on the secondary data sources to address the objectives of the study.The data was applied to characterize the hydro-meteorological drought and its trend in the city of Addis Ababa.

Data analysis tools
The SPI and RDI indices were employed to examine the occurrence of drought history at a longer time scale.Both indices were computed using the Drought Indices Calculator (DrinC) software.Meanwhile, the descriptive statistical techniques such as mean, standard deviation, minimum (min), maximum (max), and coefficient of variance (CV), and also annual and seasonal graph were computed with the aid of Microsoft excel and SPSS software package.The analyzed data were used to detect the historic rainfall and drought indices classification.

SPI and RDI drought indices
The RDI and SPI were used to find the long-term drought trends as well as to examine the occurrence of drought history at a longer time scale.Figure 3 shows the graph for mean precipitation of Sendafa, Bole, and Observatory stations from 1983 to 2015, whereas, Table 3 shows annual and seasonal SPI values for Bole, Sendafa, and Observatory stations from 1983 to 2015. Figure 4a, b, and  c show annual and seasonal precipitation series from 1983 to 2015 for Bole, Sendafa, and Observatory stations respectively.Additionally, Figure 8 shows the Standardized 12-month RDI value (RDIst_12) for Bole and Observatory stations.
From climatological or hydrological perspective, rainfall during the summer season is very common and it is important for deciding the total amount of rainfall received annually.In Ethiopia, the rainy season is usually during the summer ("Kiremt" season) and as Figure 3 shows, the average summer precipitation was very high for Sendafa, Bole, and Observatory stations in comparison with other time scales (i.e Autumn, Winter, and the Spring season).The result found in this study is in line with the results of a previous study whose outcomes show that in the study area, the seasonal cycle shows more rain in May and June [42].The SPI values for 12 months for Bole Station in Addis Ababa was high in 1995 and the result found in this SPI value plot agrees with the outcomes of standardized RDI indices (RDI st) for the year of 1995 for Bole station in Addis Ababa (Figures 5a and Figure 8).SPI can be computed for different time scales.The current study computed SPI related to mediumterm cumulated values (3 to 12 months, that is SPI 3 and SPI 12) and the SPI indices is also recommended by other authors as it is more appropriate for measuring the impact on stream flow and reservoir storage [5,28].Figure 5a, b, and c show the plot of SPI values for 12 months for Bole, Sendafa, and Observatory Station in Addis Ababa city.Likewise, the SPI values for 3 months range (seasonal) and 12-month (annual) SPI values for Bole, Sendafa, and Observatory Station in Addis Ababa city is computed as shown in Table 3.
Furthermore, Figure 6 shows the plot of the initial value of RDI_12 (α12) for a series of 33 years for Bole and Observatory station.The dashed line represents the average of α12 for each station.The research finding suggests the Reconnaissance Drought Index (RDI) as a climatic index for detecting possible drought classification for a given data.Moreover, the RDI is more advantageous than SPI since the RDI incorporates both temperature and precipitation data in a single index whereas the SPI includes only Precipitation (rainfall) data (Figures 5-7, Tables 1 and 2).Similar to this study, for drought monitoring in given stations, Jamshidi et al. [43] showed that the RDI is more sensitive than the SPI to climatic conditions and therefore the authors recommended that the role of evapotranspiration is very important in drought assessment and should not be ignored.However, many scholars from different countries such as Algeria, India, USA, Nepal, Tunisia, Kuwait, etc., also recommend the importance and suitability of SPI and RDI indices for drought monitoring, assessing and comparing for meteorological and hydrological droughts [44][45][46][47][48][49].
Figure 7 revealed that both SPI and RDI show an identical plot/ graph.Both indices did not show large difference and thus in the present study the primary driver of the severe drought is both rainfall and temperature even if RDI is more sensitive to drought severity.The meteorological variable of both rainfall and temperature are interrelated (Figure.

Drought severity classification
The RDI value is used for drought severity classification (Table 2).Therefore, Figure 9 shows the pie chart of drought conditions for Addis Ababa (Ethiopia) expressed by annual RDIst for 33 hydrological years and each drought class represented by a different color.The drought severity classification for the normal condition for Bole and Observatory station were 85% and 83% respectively.Whereas, the moderate drought classifications were 12% for Bole and 15% for Observatory station.Similarly, the severe drought classification for both stations (i.e.Bole and Observation) were 3% (Figure 9).The foregoing result concurs with the findings of an earlier study on this subject matter, which showed that the temperature trend was significantly increased due to climate change in Bole station (Addis Ababa) [10].Accordingly, the severe drought and moderate drought indices due to reduction of rainfall, climate change and other variables can lead to a further indication of a shortage of urban water supply.The calculated SPI and RDI values were classified based on the drought classification range (Tables 1 and 2).However, Table 4 shows the drought classes classified by SPI and RDI values which are analyzed from the monthly rainfall and temperature data and therefore can be used for drought monitoring.The results indicate that severe drought conditions were observed for SPI and RDI in the year 2013 for Bole station, but for the same year, the drought condition for Observatory station was normal for SPI and RDI indices.Similarly, the SPI and RDI indices of drought classification for all stations were medium droughts for the years 1991 and 2002.On the other hand, for the year 1996,  4 and Figure 8).Likewise, the SPI and RDI indices of drought classification for all stations were at normal conditions for the year 2015.The results found in this study agrees with the discoveries of a previous study [42] on this subject matter, which showed that for the year 2015 in Addis Ababa city (no significant trend), the seasonal cycle shows more rain in May and June than in the same months in Combolcha (risk ratio of 2) and this occurs from the difference in location.The SPI and RDI of drought classification for both stations (i.e.Bole and Observation) in Addis Ababa were normal drought conditions for the year 2004.In contrast to this result, NMSA [50] reported that in Addis Ababa, Robe, Adigrat etc., for the year 2004, the extreme minimum temperature below 5 °C were displayed and therefore this minimum temperature caused crop damage over same areas of Ethiopia (i.e.central, eastern, and northeastern Ethiopia).Also, as shown in Table 4, the SPI and RDI indices for the year 1984 and 2015 indicated that normal drought conditions were observed for Bole and Observatory stations (Addis Ababa).In contrast, our result disagrees with the findings of an earlier study in this domain which opined that the year 1984 was extreme drought whereas in 2015 the extremity and severity of drought is very high in the upper part of Awash River basin, Ethiopia [13].The differences in the findings of these studies can be partly attributed to the climatic and geographical variations.

Rainfall variability
Table 5 shows the statistical summary of annual and seasonal rainfall at the three meteorological stations.The average annual minimum precipitation for Sendafa, Bole, and Observatory stations are computed as 836.5 mm, 817.5 mm, and 999.2 mm respectively.Similarly, the average annual maximum precipitation for Sendafa, Bole, and Observatory station was 1481.4 mm, 1466.8 mm, and 1624.7 mm respectively.In Ethiopia, hydrological time starts from September and ends in August.Hence, the rainy season usually occurs in the summer ("Kiremt" season).Table 5 shows that the average summer precipitation was 778.58 mm, 399.12 mm, and 709.52 mm for Sendafa, Bole, and Observatory station respectively.Whereas, the time scale of spring ("Bega" season) are low raining season and hence the average spring precipitation for Sendafa, Bole, and Observatory station were 41.81 mm, 46.86 mm, and 53.95 mm respectively (Figure 4 and Table 5).These are in line with the findings of the previous study of Temam et al. [51], whose results show that the "Bega" season (from October to January) has the lowest rainfall amounts compared to other seasons in the given year (~2% of the annual precipitation of the country).Similar assumptions is obtainable in the work by Cermak et al. [52].Table 3 and 5 show SPI values and seasonal rainfall at the three Meteorological Stations calculated for the 3-month period (i.e.Sep. to Nov. for autumn, Dec. to Feb. for winter, Mar. to May for spring, and Jun. to Aug. for summer).Similarly, Table 5 revealed that Observatory station, Sendafa Station, and Bole station have highest annual mean rainfall and highest standard deviations.In comparison, similar results were obtained by Yahaya et al. [6], who found that the highest monthly mean rainfall and highest standard deviation were recorded in Bauchi station, Ibi station, and Gombe station (Nigeria)-indicating high variability of annual rainfall for the study year.
Larger coefficient of variation (higher variability) of rainfall during winter seasons (scarcity of rainfall observed) were 95.82%, 95.93%, and 77.91% for Sendafa, Bole, and Observatory station respectively.Whereas, the lower coefficient of variation during annual rainfall (this is most commonly affected with the summer season) was 15.59% for Sendafa station, 14.38% for Bole station, and 13.98% for Observatory station.As a result, this statistical summary gives a good snapshot of precipitation variability for the region and this has important implications for urban water management while showing an intensification of the adverse impact of climate change.Table 5 also indicates a larger coefficient of variation (higher variability) of rainfall during winter seasons over all stations.

Conclusion
The study indicated that severe drought conditions were observed for SPI and RDI indices in the year 2013 for Bole station, but for the same year, the drought condition for Observatory station was normal for SPI and RDI indices and this happens due to the temperature and latitude variation of the station.In a similar bid, the SPI and RDI of drought classification for all stations show medium droughts for the years 1991 and 2002.On the other hand, for the year 1996, the SPI indices showed moderate drought conditions for the Observatory station.Whereas for RDI indices for 1996, the drought classification indicated severe drought for the Observatory station.The drought classification difference in the year 1996 occurred due to the indices method difference.It can therefore be concluded that the long-term drought trend analysis of annual precipitation shows that almost three percent of severe drought and twelve percent of moderate drought were recorded in the study area.The severe and moderate drought indices due to the reduction of rainfall, climate change, and other factors can lead to an extra indication of urban water supply shortage.Lack of analyzed long term drought trends were observed in many African countries and this remains the limitation of water sector professionals in managing the water supply sector.Thus, the results of this study will be of importance to water sector professionals in forecasting weather variations and better management of urban water resources.
Further research study is recommended to conduct drought analysis using the application of soil moisture which is used to estimate the severity of agricultural droughts particularly for sub-Saharan African countries.

Figure 3 .
Figure 3. Mean precipitation of Sendafa, Bole, and Observatory stations from 1983 to 2015.(Time scale classification indicated for autumn: Sep. to Nov., winter: Dec. to Feb., spring: Mar. to May, summer: Jun. to Aug., and time scale for annual is the sum of all 12 months).

Figure. 4 (Figure. 5 (
Figure. 4 (a) Annual and seasonal precipitation series for Bole station from 1983 to 2015, (b) Annual and seasonal precipitation series for Sendafa station from 1983 to 2015, (c) Annual and seasonal precipitation series for Observatory station from 1983 to 2015.

Figure 6 .
Figure 6.The initial value of RDI_12 (α12) for a time series of 33 years for Bole and Observatory station (dashed line represents the average of α12 for each station).

Figure. 7 (
Figure.7 (a) Comparison of SPI and RDI values for Bole station, (b) Comparison of SPI and RDI values for Observatory station.

Figure 9 .
Figure 9. Drought condition of a region (Addis Ababa, Ethiopia) expressed by annual RDIst for 33 hydrological year (each drought class represents by different color).
displayed moderate drought conditions for the Observatory station.Whereas for RDI indices for the year (1996), the drought classification indicated severe drought for the Observatory station (Table

Table 3 .
Annual and seasonal SPI values for Bole, Sendafa, and Observatory stations from 1983 to 2015.

Table 4 .
Drought Classification based on SPI_12 and RDI_12 indices.

Table 5 .
Statistical Summary of annual and seasonal Rainfall at three Meteorological Stations.