Hydro-Climatic Variability: A Characterisation and Trend Study of the Awash River Basin, Ethiopia

: The objective of this study was to characterize, quantify and validate the variability and trends of hydro-climatic variables in the Awash River Basin (ARB) in Ethiopia using graphical and statistical methods. The rainfall and streamﬂow trends and their relationships were evaluated using the regression method, Mann–Kendall (MK) test and correlation analysis. The analysis focused on rainfall and streamﬂow collected from 28 and 18 stations, respectively. About 85.7% and 75.3% of the rainfall stations exhibited normal to moderate variability in annual and June to September rainfall, respectively, whereas 96.43% of rainfall stations showed high variability in March to May. The MK test showed that most of the signiﬁcant trends in annual rainfall were decreasing except in two stations. These research ﬁndings provide valuable information on the characteristics, variability, and trend of rainfall and streamﬂow necessary for the design of sustainable water management strategies and to reduce the impact of droughts and ﬂoods in the ARB.


Introduction
In recent decades, concerns around water resources and climate change have gained global significance. Water resources are exposed to high spatial and temporal hydro-climatic variability, leading to major challenges affecting the economic development [1]. Climate change and an altered water supply regime have necessitated the need to investigate the relationships between climate and water resources. The availability and accessibility of water resources are mainly reliant on climate and environment of the area. The rainfall variability, alteration in temperature and other climatic variables significantly affect the variation in streamflow level and time of hydrologic events [2]. The level of streamflow will be affected by several factors which include rainfall characteristics, soil physical properties, watershed factors and human activities [3].
Water resources in Ethiopia are highly vulnerable to climate and hydrological variability due to its topography and human-induced factors, such as land degradation, high population density and water management practices [4][5][6]. Ethiopia has twelve river basins, out of which the Awash River Basin (ARB) is intensively utilized and environmentally vulnerable. Freshwater biodiversity, basin habitat and human livelihood are exposed to a risk due to over-abstraction of river water, population growth and poor water management [7]. The Intergovernmental Panel on Climate Change (IPCC) Africa report indicated that in the highlands of Ethiopia, where there is a high and heterogeneous topography, rainfall and floods are expected to increase during the final period of the 21st century [1].

Study Area
The study was conducted in the ARB of Ethiopia which is located between latitudes 7 • 53 ′ N-12 • N and longitudes 37 • 57 ′ E-43 • 25 ′ E ( Figure 1). The basin has an area of 112,696 km 2 , with an annual flow of 4.9 billion cubic metres (BM 3 ) and 0.14 BM 3 estimated groundwater potential [20]. The total length of Awash River is 1200 km and starts flowing at an altitude of 3000 m in the Central Highlands of Ethiopia and flows through different locations, and finally joins Lake Abe at a 250 m altitude above mean sea level (AMSL). The basin covers diverse parts of the country, namely: Amhara, Oromia, Afar, Dire Dawa and Addis Ababa. The ARB has a potentially irrigable land of about 206,000 ha [16].
There are three main dams in Awash River Basin-Koka, Kessem and Tendaho. Koka Dam is located in the upper Awash and is important for hydropower generation and irrigation development. The Koka reservoir is a main source of water supply for users of downstream dams from Amibara to Gewane. Kessem and Tendaho dams are located in the Middle and Lower Awash, respectively, for irrigation and household water supply. Major crops grown in the area are sugarcane, vegetables, fruit and cotton. The Awash River gives regular irrigation water to productive and larger agricultural areas. Moreover, it generates hydropower energy and provides water supply for different towns and cities along its course, which includes Legadadi, Gefersa, Dire and Aba-Samuel reservoirs. A number of large and small irrigation schemes are found in ARB, which includes Wonji sugar factory and Wonchiti irrigation project in Upper Awash (around 33,900 ha) and Tendaho irrigation project (Dubti, Assayita, Afambo vegetable and sugarcane corporation) in Lower Awash (around 62,500 ha). Large and small irrigation schemes in the Middle Awash (around 55,000 ha) include Abadir and Merti canal, Amibara irrigation, Tebli irrigation project, Nura hera, Gewane and Kessem [21,22]. = Figure 1. Location of the study area.

Data Pre-Processing
Digital Elevation Model (DEM) (30 m resolution) of ARB and the hydrological shape files obtained from the Ministry of Water, Irrigation, Electric and Energy (MWIEE) in Ethiopia were used in ArcGIS to delineate the study area. The long-term hydro-climatic data were collected from the National Meteorology Agency (NMA), MWIEE and Ethiopian Institute of Agricultural Research (EIAR). The current analysis focused on the two hydro-climatic variables: rainfall and streamflow collected from 28 and 18 stations, respectively ( Figure 1). The rainfall and streamflow stations were selected based on three criteria: long-term data availability, percentage of missing data and representativeness of the basin. The duration selected was more than 30 years, from 1954 to 2016 for rainfall except for two stations, and more than 25 years (1970-2016) for streamflow except for one station. In this study, sites with less than 10% of the missing data were selected for both rainfall and streamflow stations, but three rainfall stations and two streamflow stations were accepted with missing values of less than 15% in order to have a better representative station distribution within the basin. Missing gaps greater than 5% were completed using normal ratio method and simple linear regression. The normal ratio method was used to complete missing data in rainfall stations that have a closer nearby rainfall station (Equation (1) [23]).
where P x and P i are annual average rainfall at a station with missing value and neighbouring gauges respectively, P i is the rainfall in the neighbouring gauge and n is the number of neighbouring stations. Primary data visualization of rainfall and streamflow was conducted with time series plots, and then normality, homogeneity and outlier tests were performed using Excel, Sigma XL and XLSTAT software (SigmaXL Version 8. In the present study, annual, monthly and seasonal data were arranged from daily data of rainfall and streamflow. The basin was subdivided into three based on their location and physical factors: Upper Awash covering the upstream parts from Koka Dam; Middle Awash including areas between Koka and Awash station; and lastly, Lower Awash covering stations between Awash, Tendaho and the eastern stations [24,25] (Figure 1). The stations included in each sub-basin are shown in Appendix A: Table A1. A similar time series was required to evaluate the relation of inter-annual and temporal variability over different locations. One of the major challenges in ARB is less data availability with varying duration. The rainfall and streamflow data have a different period with a different start and end date. In order to address this challenge and to get representative stations within ARB, stations with the similar period were clustered into one group. Therefore, depending on the duration of the time series, streamflow and rainfall stations were clustered into three and four groups respectively (Appendix A: Tables A2 and A3). The period of data used in each rainfall station groups was 60 years (1954-2013) for Group 1, 52 years (1965Group 1, 52 years ( -2016 for Group 2, 30 years (1975Group 2, 30 years ( -2004 for Group 3 and28 years (1985-2012) for Group 4. The period of data used in each streamflow station groups was 35 years  for Group 1, 26 years (1983Group 1, 26 years ( -2008 for Group 2 and21 years (1981-2001) for Group 3. Furthermore, the rainfall and streamflow gauge stations, having both rainfall and streamflow data, were selected to analyse the relationship between rainfall and streamflow. Moreover, Koka and Kessem streamflow stations are regulated at Koka and Kessem reservoirs, respectively, while Tendaho and Mille stations are monitored at Tendaho reservoir.

Statistical Analysis
The hydro-climatic data analysis was undertaken to characterise, quantify and validate the variability and trend in the basin. The variability and trend can be tested using different statistical and graphical methods. In this study, the coefficient of variation (CV), standardised anomaly index (SAI) and graphical methods were used. The rainfall and streamflow trends and their relationships were evaluated using the regression method, Mann-Kendall (MK) test, Sen's slope estimator and correlation analysis (see the subsection below for further details).
The CV is a statistical measure that shows the variation of data about the mean; it is the ratio of standard deviation to the mean (Equation (2)). In this study, a CV of <20, ≥20≤30, and >30 was considered as normal, moderate and highly variable, respectively. In previous studies, the CV was used to characterise the variability of rainfall [8,15,26].
where σ and X denote standard deviation and mean of rainfall, respectively. The SAI indicates a measure of distance between the data and its mean (Equation (3)). It shows the probability of existence of an observed rainfall amount related with the long-term mean rainfall. Examples of previous studies that have used SAI to analyse rainfall variability include [12,26,27].
where x is rainfall, X and σ denote the mean standard deviation of the rainfall data.

MK Test and Sen's Slope
There are different parametric and non-parametric trend-detecting techniques. A non-parametric trend test, the MK test [28,29], is the widely preferred method, due to its independence in the data distribution and not being affected by missing values and outliers [30]. The probability of MK trend test to give a significant trend increases when serial autocorrelation exists; therefore, lag1-autocorrelation was tested for rainfall and streamflow on annual and seasonal periods before Mann-Kendall trend test [30,31]. A non-significant autocorrelation was obtained in annual rainfall and streamflow for twenty-two and fourteen stations, respectively. For stations that showed significant autocorrelation, the trend was tested and verified using linear regression. Once the hydro-climatic trends of stations were analysed, the common duration of record in the same stations were selected to investigate the relationship of their trends using correlation. The test statistic S was calculated by: where x j denote the consecutive data values, n denotes the number of the data set, and The variance, Var(s) was calculated by where n denote the number of observations, m is the number of tied groups and t i is the number of ties of the sample at time i. The standard test statistic Zs was calculated using: In a Zs test, the null hypothesis implies no trend, whereas the alternative hypothesis means a significant change has occurred over time [30,31]. The values of Zs > 0 and Zs < 0 indicate an increasing and decreasing trend, respectively. Moreover, a significant trend at 0.05 and 0.01 significance levels exist when the /Z/ > 1.96 and /Z/ > 2.58, respectively [32]. In this study, 0.05 and 0.01 significance levels were considered to analyse the hypothesis.
The Sen's slope estimator is a non-parametric analysis to estimate the magnitude and direction of the trend in the sample of N pairs of data [32]. The slope (Q i ) was calculated by where x j is the data values at j and x k presents the data values at k (j > k). When there is one data point in each time period, then N = n(n−1) 2 , where n is the number of time periods. But, if there are more data points then, N < n(n−1) 2 , where n is the total number of observations. The values of N are arranged from smallest to biggest. Then, the median of slope or Sen's slope estimator (Q mid ) was caluclated by , i f N is even (9) When the value Q mid > 0, it indicates an increasing trend while Q mid < 0 represents a decreasing trend.

Monthly and Seasonal Characteristics of Rainfall
The mean monthly rainfall of the stations in the Upper Awash basin varied from 6 to 223 mm in the period 1975-2004 ( Figure 2). Comparatively, the monthly rainfall was low from October to February, but started to increase in June ( Figure 2). Moreover, relatively intensive rainfall was received between June and August, with the maximum mean monthly rainfall received in July at the Tulubolo station. The minimum monthly rainfall was recorded at Tulubolo in December and in all stations the lowest rainfall occurred in November and December. In the Middle Awash, the mean monthly rainfall varied from 13 to 222 mm ( Figure 2). The mean monthly rainfall in the Lower Awash varied from 16 to 186 mm ( Figure 2). The minimum rainfall was recorded in December and the maximum in August.    The variability of rainfall in Groups 3 and 4 is shown in Figure 5. In Group 3, the annual rainfall varied from 44 mm at Mille to 1590 mm at Aleltu, while the annual temporal rainfall variability indicated a similar pattern of rise and decline in a majority of stations. For example, in 1980 and 1997, decreasing annual rainfall was recorded in all stations except Wuchale and Majete; however, there was an increase in 1988 except at the Wuchale station. In Group 4, the annual rainfall varied from 107 at Adiaytu to 1361 mm at Shola Gebeya. Stations in this group showed temporal variation with a different pattern; an increase in some and decrease in others. For instance, in 1995 six stations presented a decreasing pattern while there was an increase in three others.

Inter-Annual and Temporal Rainfall Variability
The summary of statistical analysis for the 28 rainfall stations located in ARB is shown in the Table 1. The annual skewness showed that 75% of the stations were approximately symmetrical while 17.86% were moderately symmetrical. Moreover, all stations indicated a Z-score of Skewness and Kurtosis closer to zero, indicating no significant difference from normality. In about half the stations, the annual mean rainfall had a CV of < 20. On the other hand, 35.71% and 14.29% of the stations had a CV of 20-30% (moderate variability) and >30% (high variability), respectively (Table 1). In this study, high annual rainfall variability was detected in four stations: Adaitu, Debresina, Mille and Koka dams.
The seasonal variability was also tested in order to evaluate the detailed rainfall characteristics that occurred in the period. The two rainfall seasons selected were June to September (JJAS) and March to May (MAM). The JJAS rainfall showed slightly similar result as that of the annual rainfall variability. Out of twenty-eight rainfall stations, 21.43% of the stations showed a CV which was less than 20%, while 53.57% of stations indicated a CV between 20-30% and 25% of stations showed a CV which was greater than 30%. Seven stations showed high rainfall variability in the JJAS seasonal period, and, out of these, four stations indicated the same trend in the annual time series. However, the period MAM rainfall showed 96.43% of the stations had a CV > 30 (Table 1). One station (Abomassa) showed a moderate CV, but other stations indicated high rainfall variability in the period of MAM.   Figure 6 shows the SAI of the annual rainfall for selected rainfall stations in the ARB. The figure shows that the negative anomaly varied from 40% to 64% of the total number of observations at Holetta and Aleltu, respectively. On seasonal time series, MAM and JJAS showed a negative anomaly between 47% to 63% and 43% to 59% of the observations, respectively. The MAM rainfall displayed a higher negative anomaly at Shola Gebeya, whereas JJAS showed a higher negative anomaly at Aliyu Amba. In some stations a negative anomaly was observed in consecutive years, for example, Koka dam, Aleltu and Melkasa.

Monthly and Seasonal Streamflow Characteristics
The mean monthly streamflow of clustered river gauged stations in the ARB is shown in Figure 7. The Awash River discharge was low from November to mid-June, but started to increase in July (Figure 7). The highest flow was recorded in August at 63.6% of the stations, and in September for the rest of the stations. The maximum streamflow was 206.78 cubic metres per second (m 3 /s) in August at Hombole, while the minimum was 0.28 m 3 /s in February at Holetta. Group 2 hydrograph indicated that the minimum flow occurred during November to June before starting to increase in July. The maximum flow was 121.32 m 3 /s at Melkasedi and the minimum record was 4.07 m 3 /s at Akaki. Group 3 hydrograph indicated a slightly different pattern-the minimum flow was in December to February and then slightly increased in March and April; again it declined in May and June but increased in July. The results from most of these stations showed a low flow during November to June but the increase started in July.      Table 2 gives the summary of descriptive statistics of the 18 streamflow stations. The annual streamflow in the majority of the stations indicated a Z-score of Skewness and Kurtosis closer to zero, indicating no significant difference from normality, while four stations slightly deviated from the range.  (Table 2). Seasonal streamflow variability was tested for similar periods as rainfall seasonality (JJAS and MAM) analysis. The streamflow in MAM and JJAS showed a high CV at 94.4% of the stations. Only one station for each season showed a moderate variability, with a CV of 28.7% and 29.1% at Wonji in MAM and at Hombole in JJAS, respectively.

Rainfall
The summary of rainfall trend analysis using MK test and Sen's slope estimator is shown in Table 3. The MK analysis presented a significant (0.01 and 0.05) annual rainfall trend at seven stations, while five of the stations showed a negative trend. The positive annual significant (0.01) trend was observed at two stations, Koka Dam and Melkasa. The JJAS rainfall indicated a significant (0.05 and 0.01) trend at eight stations, while three of the stations showed a positive significant trend. The MAM rainfall showed a significant (0.05 and 0.01) trend at four stations. In MAM rainfall, only one station showed a positive significant trend at 5% significance level. Furthermore, a negative Zs value was detected in about 71.43%, 64.28% and 53.57% of the stations in annual, MAM and JJAS seasons, respectively.

Streamflow
The summary of streamflow trend analysis using MK test and Sen's slope estimator is shown in Table 4. The MK analysis presented a positive significant (0.01 and 0.05) annual streamflow trend at five stations. The highest positive annual significant (0.01) trend was observed at Kessem stations. The MAM streamflow showed a positive significant (0.05) trend at two stations and a negative significance trend at one station. In JJAS streamflow, only one station showed a positive significance trend at 5% significance level. Furthermore, a positive Zs value was detected in about 77.22%, 72.22% and 88.89% of the stations in annual, MAM and JJAS seasons, respectively. Moreover, in MAM, a negative Zs value was indicated in three stations located in the Lower Awash Basin.

Relationship between Rainfall and Streamflow
The annual rainfall and streamflow at six stations indicated a positive correlation coefficient which varied from 0.24 to 0.68 ( Figure 11). The correlation coefficient was 0.5, 0.55 and 0.68 at Holetta, Addis Alem and Mille, respectively ( Figure 11). By comparison, the correlation was low at two stations: Metehara (0.24) and Koka Dam (0.36). On the contrary, two stations (Adaiytu and Mojo) showed both a positive and negative relationship in a different time series. The current finding revealed the existence of a positive correlation in a majority of the tested stations.

Discussion
The mean monthly rainfall showed a bimodal pattern with two rainy seasons in the Middle Awash and Lower Awash, while unimodal rainfall was indicated in the Upper Awash. The two rainfall regimes in the Middle and Lower Awash occurred in JJAS and MAM, which presented the main and shorter rainy periods, respectively. Moreover, in all the rainfall stations, the minimum rainfall was recorded in December while the maximum rainfall occurred in August and September. These rainfall regimes coincide with the main agricultural seasons in Ethiopia, whereby agricultural production is predominantly reliant on rainfall. The study by [16,24] indicated a similar rainfall pattern; which was unimodal rainfall in Upper Awash and a bimodal rainfall in the Middle Awash and Lower Awash.
The finding of this study indicated the occurrence of temporal variability with an increasing pattern in some stations, and a decreasing pattern in other stations. The JJAS rainfall showed slightly similar result as that of the annual rainfall variability. About 85.7% and 75.3% of the rainfall stations exhibited normal to moderate variability in annual and JJAS rainfall, respectively, whereas, 96.43% of rainfall stations showed a high variability in MAM. The rainfall variability has a direct influence on agricultural water availability and irrigation, especially for countries like Ethiopia, where 85% of the livelihood depends on agriculture. Additionally, the high rainfall variability aggravates the frequency of flood occurrence in this basin. Similarly, the study by [24] indicated the occurrence of flood and drought due to annual rainfall variability in the ARB. Additionally, the study by [17] indicated a variability of precipitation and streamflow with 8.49% and 34% coefficients of variation, respectively. The study of [17] also showed an increasing rainfall trend in Kiremt (JJA) and a decreasing rainfall trend in Meher (September to November), Bega (December to February) and Belg (MAM). The SAI in annual and seasonal time series indicated a negative anomaly in a closely similar proportion of the observations. Furthermore, the SAI result showed that most of the years exhibited negative anomaly in MAM rainfall, indicating drier period. Moreover, a lot of large-and small-scale irrigation schemes are found in ARB, especially in the middle Awash, which could be directly affected by low water availability during MAM. The occurrence of high rainfall and streamflow variability combined with the expansion of irrigation in ARB has an influence on agriculture and crop productivity, which aggravates the prevalence of food insecurity and drought in the basin.
The results in most of the streamflow stations revealed that river discharge was low during the months of November to June, but started to increase in July. Furthermore, in all stations, the peak flow occurred in August and September. The monthly river flow was low during the dry period and it started to increase in the rainy season. This meant that the available water supply was limited during the dry season when irrigation is required, and more in the rainy season when hardly any irrigation is needed. Similar to the monthly streamflow, the annual and MAM river discharge indicated a similar pattern except in a few years with extreme high flow records. Moreover, the results showed that the annual temporal variability of streamflow exhibited a positive trend in fourteen stations and a negative trend in four stations.
The Awash River is the major water supply in the ARB that provides regular irrigation water to productive and larger agricultural areas. The high variability during the period MAM and JJAS affects the streamflow level, frequency and the occurrence of extreme hydrological events. As the streamflow variability becomes higher and the water demand increases, water shortage becomes a major threat unless there are adequate water conservation structures and proper water management approaches in this basin. In recent years, there has been significant expansion of irrigation activities in this basin, especially in Upper and Middle Awash region. Furthermore, flood and drought have become major challenges in the ARB [5,16]. The MK trend test indicated a significant (0.01 and 0.05) rainfall trend in seven, four and eight stations in annual, MAM and JJAS seasons, respectively. The MK test showed that most of the significant trends in annual rainfall were decreasing except in two stations. Moreover, the negative Zs and Sen's slope estimator values were detected in about 71.43%, 64.28% and 53.57% of the stations in annual, MAM and JJAS seasons, respectively, indicating a decreasing rainfall trend. Both the graphical and MK test showed a decreasing trend in most of the rainfall stations, indicating a decline in water availability in ARB.
In addition, the MK test showed a significant (0.01 and 0.05) streamflow trend in five, one and three stations in annual, JJAS and MAM seasons, respectively. Furthermore, a positive Zs and Sen's slope estimator values were detected in the majority of the stations in annual, MAM and JJAS seasons, indicating an increasing streamflow trend. However, in MAM, a decreasing streamflow trend was detected in Lower Awash at three stations. The streamflow trend in ARB relates to the level of Awash River that contributes to the water availability in this basin. The increasing trend leads to frequent occurrence of peak flows, which aggravates flooding in the area. Furthermore, the results indicated a decreasing streamflow trend in few stations in the Middle Awash and Lower Awash. This low river flow indicated a decline in water availability, which affects the irrigation water supply and also other water use in the ARB, especially in the lower part. The Awash River flow reduction might be related to irrigation expansion in the Upper and Middle Awash, improper water-abstraction, poor water management, climatic factors, population density and environmental change.
The relationship between rainfall and streamflow trends in six stations showed a positive correlation coefficient with relatively good correlation (>0.5) in three stations and poor correlation (<0.4) in two stations. The correlation observed at the two stations, below Koka dam and Mille, might be affected due to the monitoring effect of Koka and Tendaho reservoirs, respectively. Moreover, the rainfall-streamflow trends showed that rainfall was not the only factor that influenced the variability and availability of streamflow in this basin. Further research is recommended to investigate other factors that influence variability in Awash River flow and the rainfall-streamflow relationships, for instance: the density of physical properties, watershed characteristics and environmental factors like vegetation cover.

Conclusions
The objective of this study was to characterise, quantify and validate the variability and trends of hydro-climatic variables in the ARB in Ethiopia using graphical and statistical methods. The rainfall and streamflow trends and their relationships were analysed using regression method, MK test and correlation analysis. The main findings of this research are summarised as follows:

1.
A bimodal pattern with two rainy seasons was displayed in the mean monthly rainfall in the Middle and Lower Awash, while the pattern in the Upper Awash was unimodal. In the majority of the stations, the rainfall exhibited from normal to moderate variability in the mean annual and the JJAS period, whereas it was highly variable in almost all stations during the MAM rainfall season.

2.
High streamflow variability was observed in MAM and JJAS periods. Besides, the monthly river flow was low during the dry season of the year, but started to increase in the rainy season. Shortage of rainfall in MAM and low monthly streamflow during the dry period leads to a decline in water supply in ARB.

3.
The MK trend test showed that most of the significant trends in annual rainfall were decreasing. The negative Zs and Sen's slope estimator values were detected in about 71.43%, 64.28% and 53.57% of the stations in annual, MAM and JJAS seasons respectively, indicating a decreasing rainfall trend. Both the graphic and MK test indicated a decreasing trend in most of the stations, which may result in a decrease in water availability in ARB.

4.
The MK trend test revealed an increasing annual streamflow trend in 27.8% of the stations, whereas, in the MAM period, a decreasing trend was displayed in the Lower Awash at three stations. Moreover, a positive Zs and Sen's slope estimator values were observed in majority of the stations in annual, MAM and JJAS seasons, indicating an increasing streamflow trend. An increasing trend was shown in most of stations in Upper and Middle Awash, whereas a decreasing trend was revealed in Lower Awash.

5.
The trend relationship between rainfall and streamflow in six stations indicated a positive correlation coefficient, with slightly good correlation in four stations and poor correlation in two stations. The correlation observed downstream of Koka and Tendaho dams might be affected due to the monitoring effect of Koka and Tendaho reservoirs, respectively. Moreover, the rainfall-streamflow trend relationship showed that rainfall was not the only factor that influenced the variability and availability of streamflow in this basin. The other factors that influence variability in Awash River flow and the rainfall-streamflow relationship includes physical property, watershed characteristics and environmental factors, like vegetation cover, which still needs further assessment.
The findings of this study provide valuable information on the characteristics and variability of rainfall and streamflow necessary for planning and designing of sustainable water management strategies and minimizing the impact of droughts and floods in the ARB. The cause of an increasing trend in streamflow might be urbanisation, land degradation, population density and climatic factors. However, the declining trend might be related to irrigation expansion in Upper and Middle Awash, improper water-abstraction and poor water management, climatic factors and environmental changes in the basin. Therefore, further assessment on factors influencing the variability and availability of Awash River flow is needed in order to implement a sustainable water management strategy and also to mitigate the adverse effects of water scarcity and extreme events on the food security, livelihoods, environment and overall development of the country. Additionally, further research on river water management and water saving techniques is recommended to cope with the high variability of river flow and rainfall in this basin.