Analysis of meteorological droughts in the Lake’s Region of Ethiopian Rift Valley using reconnaissance drought index (RDI)

This study investigated the intensity, trend and spatio-temporal variability of meteorological drought in the Lakes’ Region of Ethiopian Rift Valley using monthly rainfall and maximum and minimum temperature records for the period 1986–2019. Reconnaissance Drought Index (RDI) was employed to generate the intensity of drought at 3 and 12-months timescale. Mann-Kendall trend test was used to determine the trend of the changes in the RDI time series. The spatial extent of droughts has been interpolated by inverse distance weighted (IDW) method using the spatial analyst tool of ArcGIS. Results indicated the occurrence of different intensity and trend signals across seasons and over space in the study area. A total of 33 extreme drought months were observed in all stations during summer with varying intensity (− 2.01 at Halaba to − 3.52 at Wolaita) and 168 extreme drought months at annual timescale ranging from − 2.10 at Hawassa to − 4.51 at Wolaita. The intensity of drought events observed in Wolaita in 1986 at all timescales (RDI value of − 3.19, − 3.52 and − 4.51 for spring, summer and annual respectively) were very extraordinary and devastating. Drought magnitude showed increasing signal at 6 out of 10 stations, although statistically significant at only two stations (Arsi Negelle at all timescale and Butajira at spring and annual timescale). However, the spatial patterns of drought events didn’t exhibit clear pattern rather more localized distribution and variability. The frequency of drought incidence became intense in the study area from 2008 onwards at all timescales compared to the 1990s and 2000s.The increasing tendency of drought in recent years might be the manifestation of borderless global warming. The empirical evidences showed that drought events and their negative effects are highly localized in the study area and provide useful information for local-scale planning for drought management and response.


Introduction
The warming of earth's climate system is unequivocal and the intensity and frequency of extremes are likely to increase over many areas (IPCC V 2014). The report highlighted that frequency and intensity of extremes like droughts, floods, and heat waves are expected to change as earth's climate changes. Droughts are one of the natural hazards which have devastating effects on the socioeconomic and environmental conditions. Droughts are resulted from many hydrometeorological processes that limit surface water or groundwater availability, creating conditions that are significantly drier than normal (WMO and GWP 2016).
Droughts have significant deleterious impact on many economic sectors and people at any one time. According to the Emergency Event Database, in the period 1995-2015, even though drought accounted for less than 5% of all natural hazards, more than one billion people were affected by droughts, that is more than a quarter of all people affected by all types of weather related disasters worldwide (CRED and UNISDR 2015). In this regard, drought characterization is vital, since it is the basis to find appropriate risk management measures by monitoring its intensity, frequency of occurrence and spatial distribution.
As it is indicated in Degefu and Bewket (2014) and Zargar et al. (2011), the impact of drought is governed by its intensity, duration, frequency and spatial extent of the precipitation deficit. Intensity refers to the amount of precipitation or water storage deficit at a particular place and specific time which is classified as mild, moderate, severe and extreme. Frequency is the average time between successive drought events and duration refers to the extent of time that a given drought event stays. The spatial coverage refers to the areal extent of a specific area affected by a given drought incidence. Usually severe and extreme drought episodes cover wider areas, while mild and moderate drought episodes tend to affect localized areas (Degefu and Bewket 2014).
Ethiopia is affected by frequent drought events that resulted in low crop and livestock productivity which in turn makes Ethiopia dependent on large humanitarian aid (Teshome and Zhang 2019). Ethiopia faced one of the worst droughts the country has seen in decades with over 10.2 million people in need of food aid (UN OCHA 2015). Triggered by El Niño, the rainy seasons in 2015 and 2016 failed so that the drought brought significant impact by limiting agricultural production, straining livelihoods, and exacerbating food insecurity among poor and vulnerable households. Every year, approximately three million Ethiopians are affected by crop production shortfalls adding to the 7.6 million supported every year by the Productive Safety Net Program, a social safety net supporting some of Ethiopia's poorest and most foodinsecure families (Cochrane and Tamiru 2016) because of extreme drought.
It is, therefore, very important to analyze drought events in Ethiopia which is continuously occurred over northern, northeastern, eastern and central parts of the country. Understanding the occurrence of extreme weather events especially rainfall is important for decisionmaking as it is a key factor in prediction and risk assessment (Teshome and Zhang 2019). Understanding drought conditions, societal vulnerability, and their related effects on one another provide us historical lessons that can aid in dealing with future drought conditions. Drought characterization at regional and local scales has significant implications for drought management. To help decision makers and the public address these issues, systematic assessments of meteorological drought is very crucial. The analysis and forecasting of extreme climatic events has become increasingly relevant to make planning effective. Understanding the features of droughts at the regional to local scale is a key for designing response measures that enhance preparedness and early warning systems, and additionally in the long term, it helps in the formulation of effective adaptation planning measures to reduce the related risks.
Many studies (for example, Edossa et al. 2010;Gebrehiwot et al. 2011;Viste et al. 2013;Degefu and Bewket 2014;Gizachew and Shimelis 2014;Muluneh et al. 2014;andMuluneh et al. 2016, Mohammed et al. 2017;Tesfamariam et al. 2019) were carried out to characterize the recurrent nation-wide and local-scale meteorological drought incidences in different parts of Ethiopia. However a few were conducted in this drought-prone rift valley lakes' region of Ethiopia. Gizachew and Shimelis (2014) studied vulnerability of the community to drought in rift valley areas of Ethiopia. Muluneh et al. (2014Muluneh et al. ( , 2016) studied on searching for evidences of changes in extreme rainfall indices and predicted the future precipitation and maize production in central rift valley (LRERV) of Ethiopia. However a few were conducted in this drought-prone rift valley lakes' region of Ethiopia. However, the study area is experiencing high temperature and potential evapotranspiration which has significant influence on the hydrology of the study area. In contrast this study characterized the spatial and temporal patterns of meteorological drought using Reconnaissance Drought Index (RDI) taking precipitation and potential evapotranspiration as an input for the model which help present reliable information about meteorological drought and design appropriate adaptation strategies. Therefore, the objective of this study was to analyze the intensity, frequency and spatiotemporal patterns of meteorological drought using RDI in the lakes' region of Ethiopian rift valley.

Description of the study area
The study covers the lakes' region of Ethiopian rift valley which is located between 37°25′ E and 39°28′ E and between 5°45′ N and 8°28′N (Fig. 1). It has an area of 41, 549.7 km 2 .This area is one of the most unstable and environmentally vulnerable areas of the country (Muluneh et al. 2016). It is part of the Great East African Rift Valley and covers the major dry land parts of the country. The pattern of rainfall in the rift floor is more of stormy type with relatively high intensities (Makin et al. 1975).
There are eight major lakes, namely Koka, Ziway, Langano, Abiyata, Shala, Hawassa, Abaya and Chamo (Starting from north) (Fig. 1). The region is also characterized by three main seasons. The dry period (locally known as Bega) extends between September and February having occasional rains which account for about 10-20% of the annual total, the small rainy season (locally known as Belg) occurs from March to May contributing 20-30% of the annual precipitation, and the long rainy season (locally known as Kiremt) is from June to August during which 50-70% of the mean annual total precipitation is received (Legesse et al. 2004). The characteristics of rainfall in LRERV is more of stormy type with relatively high intensity.
The soils of the rift valley are largely derived fromrecent volcanic rock. The main parent materials are basalt, ignimbrites, lava, gneiss, volcanic ash, alluvium, and pumice. Some of the soil problems include lowphosphorus levels, micronutrient imbalances, and, insome cases; poor physical structure (Makin et al. 1975).
The distribution of plants in the study area is highly influenced by elevation which also impacts the precipitation pattern (Musein 2006). The rift valley floor is largely dominated by deciduous Acacia woodland and wooded grasslands that are increasingly becoming more open, whereas deciduous woodlands (Legesse et al. 2004).

Data source
Monthly minimum and maximum temperature and precipitation data from ten stations for the period (1986-2019) were obtained from Ethiopian National Meteorological Agency (NMA). The max and min temperature were used to generate the potential evapo-transpiration data of each respective station. Initially, we collected max and min temperature and precipitation data from 18 stations found in the study area. We only used the ten stations (Table 1) which were better in their record length and missing values. The remaining stations were not selected because of their missing value and shorter record period. The selected stations are spatially well distributed across the study area (Fig. 1).

Drought indices
In order to understand the complexity of droughts, meteorologists and climatologists around the world developed and used numerous operational drought indices to characterize and monitor the intensity, duration and frequency of drought (WMO 2016;Dai 2011). Most of these indices use the records of precipitation, temperature and evapotranspiration (Steinemann et al. 2005). Palmer Drought Severity Index (PDSI) (Palmer 1965), the Percent of Normal, and the Deciles approach (Gibbs and Maher 1967)  The most important thing which should be considered in selecting the index is the availability and quality of the input data, simplicity of its calculation, and its ability to consistently detect spatial and temporal variations of a drought event (Morid et al. 2006;Keyantash and Dracup 2002).
Understanding the strength and weakness of each index is crucial before deciding to apply. For example, PDSI is considered as important tool in monitoring drought at larger scale (Steinemann et al. 2005). It uses water balance models which require large amount of climate and soil moisture data (Mpelasoka et al. 2008;Ellis et al. 2010) which in turn leads to difficulty in its computation (Degefu and Bewket 2014;Morid et al. 2006). SPI is also a well known indices used for analysing meteorological droughts using only precipitation as an input (El Kenawy et al. 2016;Degefu and Bewket 2014;Viste et al. 2012). However, in arid areas where temperature is highest, it is vital to consider potential evapotranspiration to calculate drought. This study is conducted in lakes' region of Ethiopian Rift Valley where high temperature has profound implication on the use and consumption of water resource. Therefore, it is equally important to calculate potential evapotranspiration (PET) of the study stations in order to calculate meteorological drought. Hence, Reconnaissance Drought Index (RDI) which requires precipitation and PET as an input is used for drought severity assessments in this study.

RDI computation
The RDI was developed to approach the water deficit in a more accurate way, as a sort of balance between input and output in a water system (Tsakiris and Vangelis 2005;Tsakiris et al. 2007b). It is based both on cumulative precipitation (P) and potential evapotranspiration (PET). Due to the unavailability of PET records, DrinC software is used for estimating the variable using Hargreaves method.
This software package was developed for providing a simple, though adaptable interface for the calculation of drought indices (RDI, SPI, SDI). The application of DrinC in several locations, especially in arid and semiarid regions, is gaining ground as a useful research and operational tool for drought analysis (Asadi Zarch et al. 2011;Tigkas et al. 2012).
In this study, RDI values were computed for each month of the year using DrinC (Drought Indices Calculator) software developed at the Centre for the Assessment of Natural Hazards and Proactive Planning and the Laboratory of Reclamation Works and Water Resources Management of the National Technical University of Athens (available at http://drinc.ewra.net/). RDI was designed to quantify the precipitation deficit for multiple timescales. Although RDI can be calculated from 3 month up to 12 months (Tsakiris et al. 2007b), this study computed the RDI values at two time-scales, i.e.3 months (RDI-3) and 12 months or annual (RDI-12). The RDI-3 was used to assess droughts during spring (locally known as Belg) and summer (locally Kiremt) seasons which represent the shorter and longer rain seasons respectively and RDI-12 was used to assess the annual drought. Positive RDI values indicate greater than mean precipitation and negative values indicate less than mean precipitation. Positive values of RDI st indicate wet periods, while negative values indicate dry periods compared with the normal conditions of the area. According to Tigkas et al. (2015), RDI categorized drought severity in to mild, moderate, severe and extreme classes, with corresponding boundary values of RDI st (− 0.5 to − 1.0), (− 1.0 to − 1.5), (− 1.5 to − 2.0) and (< − 2.0), respectively.
The RDI can be formulated in alpha (RDI αk ), normalized (RDI n ) and standardized (RDI st ) forms. The RDI st is used for drought severity evaluation, whereas RDI αk can be applied for assessing the aridity index. Basically, two expressions need to be evaluated viz., initial αk and standardized RDI, i.e. RDI st(k) . The initial RDI αk for year i and reference period of k (months) is computed using a monthly time step as described in Tigkas et al. (2015).
in which Pij and PETij are the precipitation and potential evapotranspiration of the j-th month of the i-th year and N is the total number of years of the available data. The standardized form of RDIst is computed by fitting a log-normal probability density function to the given frequency distribution of αk (Tigkas, 2008). The expression for the Standardized RDI is given as: in which y(i) is the ln (αk (i)), y is its arithmetic mean and σy is its standard deviation. According to Tsakiris et al. (2008) and Asadi Zarch et al. (2011) it has been proved that the calculation of RDIst could be performed better by fitting the gamma probability density function (pdf) at the given frequency distribution of αk, following the procedure described below.
Where: α > 0, is a shape parameter; β > 0, is a scale parameter; x > 0, is the amount of precipitation and Γ (α) is the gamma function, defined as: Fitting of distribution to data requires the estimation of α and β as follows: For n observations, the resulting parameters were then used to find the cumulative probability of an observed precipitation event for the given month or any other time scale: Since, the gamma function is undefined for x = 0 and a precipitation distribution may contain zeros, where the cumulative probability becomes: Where q is the probability of zero precipitation and G(x) is the cumulative probability of the incomplete gamma function. If m be the number of zeros in a αk time scales, then q could be estimated by m/n.
Finally the cumulative probability H(x) is transformed to the standard normal distribution with mean zero and the variance of one (Abramovitz and Stegun, 1965), which is the value of RDIst (Tsakiris et al. 2008).
Figures 4, 5 and 6 present the intensity of drought in each year of the observation period. However, Table 2 presents the total number of drought months in each timescale during the entire observation period. Once the intensity of drought is known in each timescale in the entire observation period, we multiply it by the number of months each time scale covers (i.e multiply the number of intensity by 3 for spring and summer, by 12 for annual timescale) to obtain the total number of drought months in each timescale.

Mapping spatial distribution of drought incidences
The RDI values produced by DrinC software were used as an input to ArcGIS to generate drought severity maps for the study area at 3 and 12-month time scales. To assess the spatial extent of droughts in the study area, RDI time series values of each meteorological station were interpolated by Inverse Distance Weighted (IDW) method using the Spatial Analyst tool of ArcGIS. The IDW method gives better representation for interpolation of precipitation distribution over heterogeneous topographic terrain (Tagel et al. 2011) (Fig. 2).

Trend analysis
In order to detect the changes in the RDI time series, trend analysis is carried out using nonparametric Mann-Kendall trend test. The Mann-Kendall test (Yue et al. 2002) is appropriate for this analysis, since it is distribution-free, can accommodate outliers, is more suitable for non-normal, incomplete, and missing data, which frequently occur in hydro-meteorological studies. The standard Pvalues which will be obtained from it was based on an assumption of independence between observations The test was performed for three timescales according to the reference periods used in the calculation of drought indices.
The Mann-Kendall test statistic S is calculated as follows: where n is the length of the data set, x j and x k are the annual values in years j and k, j > k, respectively, and For samples greater than 10 the normal approximation test is used (Salmi et al. 2002). First the variance of S Var(S) is computed by the following equation that takes ties into account: Here q is the number of tied groups and tp is the number of data values in the pth group. The standardized test statistic Z was calculated using S and var.(s) values as follows: The value of Z in the Mann-Kendall test statistics follows standard normal distribution with average of 0 and variance of 1. The presence of a statistically significant trend is evaluated using the Z values. A positive value of Z indicates a rising trend and vice versa. In a two-tailed test for trend, the null hypothesis Ho (no trend) is rejected if the absolute value of Z is greater than Z1 − α/ 2, where Z1 − α/2 is obtained from the standard normal cumulative distribution tables at α = 0.001, 0.01, 0.05 levels of significance. The computations were undertaken using the MAKESENS add-in software developed by the Finnish Meteorological Institute (Salmi et al. 2002).

Magnitude and frequency of drought events
Lakes' region of Ethiopian rift valley (LRERV) experiences bi-modal rainfall (spring locally called belg and summer locally called kiremt) occurring between March to May and June to August respectively. In the study area, the average contribution of the belg and kiremt rain to the annual mean was 71% in the study period in which they contributed 28 and 43% to mean annual rainfall respectively. Our meteorological drought assessment therefore was for these two rainy seasons and for the annual timescales. Because, the rain-fed agriculture is strongly dependent on the bi-modal rain during these seasons in the study area and generally in the country (Viste et al. 2012;Funk et al. 2015;Tesfamariam et al. 2019). As a result, failure of rainfall in these growing seasons leads to drought incidences (Fig. 3).
The temporal and spatial characteristics of drought in LRERV were identified from RDI time series of multipletime steps. In this study, RDI for 3-months (belg and kiremt) and 12-months (annual) time steps are computed to examine the characteristics of drought in short and medium time periods. Drought happens every time when RDI is negative and its intensity comes to − 0.5 or lower (Tsakiris et al. 2007a). Drought frequency, in this study, was measured by the number of years which experienced negative DRI values in the total time series. Varied drought magnitude classes were identified in 3and 12-months time scales across all stations in each year of the observation period (Figs. 4, 5 and 6) and total number of drought months in each timescales during the observation period (Table 2).

Spring season
The number of drought years with mild, moderate, severe and extreme severity classes computed at 3 month timescale in the study stations was varied (Fig. 4). The total number of drought years was highest in Arsi Negelle, Wolaita, Wondogenet and Ziway stations (Table 2). However, Butajira experienced two extreme droughts (RDI value of − 2.05 and − 2.14) and Hawassa, Kulumsa, Melkassa Wolaita and Wondogent stations experienced one extreme drought with the highest in Wolaita (− 3.19) (Fig. 4). Although the magnitude was varied (Fig. 4), 1999 and 2009 were drought years across all studied stations. Arsi Negelle was experiencing negative RDI values since 2011 which indicates the consecutive occurrence of meteorological drought.

Summer season
The total number of drought months at 3 month timescale in the entire period of analysis was found between 21 and 36 months ( Table 2). The total number of drought years was highest in Arsi Negelle station during summer season. Except Hawassa, all studied stations experienced extreme drought incidents during summer season. The year 2009, similar to spring season, was the driest year in all studied stations except Melkassa and Kulumsa (Fig. 5). The year 2015 was also a drought year of which seven out of ten stations had a negative RDI values ranging from − 0.65 to − 2.33. Arbaminch, Wolaita and Ziway stations experienced an extreme drought event with RDI value of − 3.04, − 3.52 and − 3.04 respectively (Fig. 5). The percentage of drought occurrences at this season was highest in Arsi Negelle (44%) and Hawassa (40%) (Fig. 5). This implies that the areas experienced a drought episode in every 2 to 3 years.

Annual drought
The total number of drought at annual timescale was highest at Arbaminch (144) and Ziway (168). The percentage of drought events at these two stations was 38% and 41% respectively. This indicates that these two areas experienced a drought episode in every 2 to 3 years. Except Arbaminch, all studied stations experienced extreme drought event at annual time scale. However its magnitude was extremely highest at Wolaita with RDI value of − 4.51. Moreover, three extreme droughts events were occurred for three consecutive years (2015-2017) with RDI value of − 2.12,-2.37 and − 2.86 at Butajira. The year 2015 was the worst drought year in the country of which six out of ten stations in the study area also showed extreme drought intensity (Fig. 6). From the total number of extreme drought events (14) registered in all stations in the entire study period , 11 extreme drought incidents were recorded from 2013 to 2017. This indicates that the magnitude of droughts has been increasing in the recent years. The year 2009 was a registered large-scale drought year which was common to the study area and other parts of the country. All stations except Melkassa and Kulumsa experienced drought events in 2009. Similarly, the year 1999, which was known for its dryness in all parts of the country, was also a drought year in six of the stations (Fig. 6). On average, the study area experienced at least moderate drought magnitude in every 2.45 years in the belg, 2.25 years in the kiremt seasons and 2.7 years interval at annual time scale during 1986-2019.

Spatial distribution of drought occurrences
Time series (1986-2019) of drought indices generated by employing the RDI procedure was used as an input to ArcMap/GIS to show the spatial extent of drought events in the study area at multiple time scales.
The spatial patterns of drought at different time steps (3and 12-months) were depicted in 8, The general patterns of drought events in the study region were interpolated based on RDI values of the ten stations under the study period. The frequencies of drought for different intensity classes and timescales show complex and localized patterns. During spring, the southern part of the study area was prone to mild drought intensity with the average 20% probability of drought occurrence. In the same season, large parts of the study area also experienced severe drought episodes (on average 3 during the study period). Typically, Arbaminch, Wolaita Sodo, Arsi Negelle (being the highest), Ziway and Kulumsa stations had severe drought event. During summer, there was no considerable spatial pattern of drought except the mild and moderate droughts events. Southern part of the study area (Arbaminch and Wolaita Sodo) and the middle stations had greater exposure to mild and moderate droughts respectively. At annual time scale, the central and northern parts of the study area experienced moderate drought (the highest in Arsi Negelle and Kulumsa). Generally, the total frequency of drought was highest in Arsi Negelle, Wolaita, Wondo Genet and Ziway in spring season. Whereas Arsi Negelle and Hawassa had highest number of drought events in summer. At the annual time scale, Ziway experienced the highest number of drought events.

Trends of drought occurrences
Trends of drought occurrences for 3-spring and summer and annual timescales were shown in Table 3. The Mann-Kendall trend test during spring showed decreasing changes in RDI values in 6 out of 10 stations (Table 3). This indicates that drought incidence during spring showed an increasing tendency. However, significant increasing trend of drought occurrence was observed at Arsi Negelle and Butajira stations (Table 3). Similarly, the RDI values during summer showed increasing changes in most of the stations. However, significant increasing trend of drought was observed only at Arsi Negelle. Similar significant increasing trend of drought was also observed (annual) timescale at Arsi Negelle and Butajira stations. Although trends at all time-scales were not statistically significant, most of the stations (Table 3) showed increasing tendencies of drought event. Generally more meteorological droughts were observed during spring and annual timescales than summer (Table 3).

Discussion
In this study, we presented a brief drought analysis at annual and seasonal time steps using RDI. The study monitored the intensity, frequency, trend, spatial distribution of drought events over LRERV of Ethiopia during the experimental years of 1986 to 2019 at 3-and 12-month time scales. As the study area experiences bimodal rainfall (spring and summer) (Tesfamariam et al. 2019) occurring between March to May and June to August respectively, we characterized RDI at three-month timescale for both seasons.
Results indicate that occurrences of drought in LRERV showed temporal variation. The highest number of drought events (168 months) was recorded in Ziway at annual timescale which accounts for 40% of the total number of months over the study period. However, during spring, all stations had comparable number of drought episodes ranging from 24 to 33 months (Table 2). There were more number of total drought events in summer season than spring at Halaba, Arsi Negelle and Hawassa. The total number of drought months during summer in Arsi Negelle, Hawassa and Halaba in the entire study period was 36, 33 and 30 respectively ( Table 2).
The prevalence of drought during summer (the main growing season) has multifaceted implications on crop and livestock production, surface hydrology and biodiversity of the study area (Mohammed et al. 2017). Failing to rain at this season means failure of 85-95% of the food crop which is produced in this season (Degefu 1987). Studies indicate that reduction of Ethiopian summer rainfall and the subsequent drought occurrence were highly associated with fluctuations in pressure gradient and sea surface temperature (SST) (Degefu and Bewket, 2017;Diro et al. 2011;Korecha and Barnston 2007;Haile 1988). The development of more warm areas in the equatorial Pacific Ocean attributed to an increase in local precipitation over that ocean and reduced rainfall over Eastern Africa, including Ethiopia (Funk et al. 2014). When El Nino prevails in the equatorial Pacific Ocean, the summer rainfall in the country in general declines to below average, causing drought events (Degefu and Bewket, 2017;Fekadu 2015). In contrast, La Nina events favor further temporal expansion of the summer rain beyond the normal seasonal period of a given area. The clear association between ENSO and the belg rain in Ethiopia was also reported in many studies. During La Nina, the main recipients of spring rain in the country experience shortage of rainfall while during El Nino situation, most of the spring rainfall recipients receive higher amount of rainfall (Tesfamariam et al. 2019;Fekadu 2015). Thus, the frequent drought incidence in the study area during spring and summer seasons could be primarily linked to La Nina and El Nino situations, respectively (Tesfamariam et al. 2019).
The findings indicate that the frequent drought episodes had varying intensities in the study area at all studied timescales. Extreme drought incidents were recorded in almost all stations at summer and annual Compared to other stations, Butajira and Melkasa experienced more number of extreme drought months (6) in this season during the entire study period (Table 2). Similarly, a total of 168 extreme drought months were observed in all stations except Arbaminch at annual timescale with varying magnitude ranging from − 2.10 at Hawassa to − 4.51 at Wolaita (Fig. 6). The most remarkable condition at this timescale is that 156 out of 168 extreme drought months were recorded in recent years (2009 onwards) with the highest in Butajira (36 months). It is again striking that 60 out of 168 extreme drought months were recorded in 2015. This year was known for its dryness due to the existence of El Nino all over the country which brought significant impact by limiting agricultural production, straining livelihoods and exacerbating food insecurity among poor and vulnerable households (Mohammed et al. 2017;UN OCHA 2015).
The other most striking event was that the intensity of drought events observed in Wolaita in 1986 at all timescales (RDI value of − 3.19, − 3.52 and − 4.51 for spring, summer and annual respectively) were very extraordinary and devastating. These had strong correlation with the occurrence of El Nino in 1986 in the country (FAO 2014). Similarly in 1986 most extreme drought incidence at Wolaita (with SPI < − 4.02) at 3-months timescale was reported by Degefu et al., (2014) and Jufare (2008). The occurrence of extreme intensity droughts at 3-month have important implications on rainfed agriculture (Morid et al. 2006;Ellis et al. 2010). The 1986 rainfall variability and drought have caused significant adverse effects on crop production and subsequent food shortages in the most densely populated areas, particularly in the Wolaita Zone (Jufare 2008). The result implies that drought magnitude can vary largely over space and this exceptional extreme drought shows a local character, where small area was affected, while others had less severe drought magnitude or even no drought. This finding is supported by the study of El Kenawy et al. (2016) who found similar drought occurrences in studying the changes in the frequency and severity of hydrological droughts over Ethiopia from 1960 to 2013. It is well known that the complex terrain regulates the Ethiopian climate by creating numerous microclimates having different amounts of rainfall (Viste et al., 2013;Jury 2010;Dinku et al. 2007;Degefu 1987).
The frequency of drought incidence became intense in the study area from 2008 onwards at all timescales compared to the 1990s and 2000s. Arsi Negelle, for example, experienced all drought incidences (189) which constituted 38% of the total number of months at all timescale from 2009 onwards. The year 2009 was one of the few years with drought conditions in all of Ethiopia, both on seasonal and annual scales. On the national level, 2009 was the second driest year in history after 1984 (Viste et al. 2012). Although drought conditions prevailed in all parts of Ethiopia in 1984 and 2009 on both seasonal and annual time scales, the droughts were of a more localized or regional character, affecting only some parts of the country. The increasing tendency of drought in recent years might be the manifestation of borderless global warming. Generally, droughts became more frequent on shorter time scales but their duration was short.
The opposite is true for longer time scales. Their frequency became less but remained for longer time.
Trends of drought events at 3-and 12-month timescales over the study period were depicted in Table 3. Negative RDI values in the MK test indicates the increasing tendency of drought and vice versa. The MK test results indicate that RDI values were ranging from 0.13 to − 0.69 Z unit per decade at the three-month (spring) timescales suggesting increasing tendency of drought episodes. Similarly, at annual timescale, the RDI values were ranging from 0.17 to − 0.73 unit per decade. The trend of drought incidence was significantly increasing in Arsi Negelle at all timescales especially pronounced in recent years and Butajira in spring and annual timescales. The increasing trends of droughts episodes at the two stations were pronounced in recent 10 years which has strong association with global warming.
The spatial patterns of drought events showed more localized distribution and variability when we considered intensity classes and timescales. For example, more frequent mild intensity drought events were observed in the central parts and south western part of the study area during spring. The northern parts of the study area experienced more extreme drought events during spring season. More number of mild droughts was observed in north eastern and south western parts of the study area during summer. Moderate drought events were more pronounced at east central part of the study area especially around Arsi Negelle and Hawassa at summer time scale. Except few pocket areas around Hawassa, all areas of the study area had substantial number of extreme intensity of drought during summer. However, except the southern part of the study area, the numbers of extreme drought incidents were pronounced all over the study area at annual timescale.

Conclusions
This study focused in evaluating the intensity, frequency, trend and spatial patterns of meteorological drought in the LRERV based on rainfall and calculated PET data from ten stations for the period 1986-2019. The study employed RDI to generate drought indices at 3 and 12months timescale. Mann-Kendall trend test was also used to detect trends of change.
Even though almost all stations in the study area were suffering from drought, the temporal distribution and frequency of droughts varied markedly among each station. Although the intensity of drought varied, all  stations experienced frequent drought incidences in recent years (2009 onwards) at all timescales. All areas in LRERV were stricken at least partially by the historical drought events but with varying frequency and intensity. It was evident that most of the stations were exposed to severe and extreme droughts at all timescales in the most known national driest years of 2009 and 2015.
Generally, more frequent droughts were observed in the study area compared to the nationally documented historical drought events, which implies higher vulnerability of LRERV to climate variability. The frequency and intensity of drought events computed at all timescales show complex and localized spatial patterns. However, the total frequency of droughts for all timescales showed that the southern and north eastern part of LRERV experienced the most recurrent drought risks during the study period. Although drought incidents were less frequent at Wolaita Sodo compared to Arsi Negelle and Ziway, the area has experienced the most extreme drought in all time scales reaching to RDI value of − 4.51 and caused significant adverse effects on crop production and subsequent food shortages in the most densely populated areas.
The Mann-Kendall trend test at all timescales indicates the tendency of increasing drought in most of the stations. However, the drought incidents were significantly increasing at Arsi Negelle at all timescales and Butajira at spring and annual timescale. The existence and occurrence of such localized droughts in LRERV suggest differences in local and regional scale rainforming factors in the region. However, this requires a further study, which should consider how local factors (topography, forest, land cover change, water bodies, etc.) affect rainfall amounts in LRERV. Generally, the implication of this study is that drought events and their negative effects are highly localized in the study area and there is a need for local-scale planning for drought management and response. Local climate studies, such as the present study, should be strengthened since they are useful to provide fine scale climate information for drought risk management and adaptation purposes.