Impact of Climate Change on the Hydrology of the Upper Awash River Basin, Ethiopia

: This study investigated the impacts of climate change on the hydrology of the Upper Awash Basin, Ethiopia. A soil and water assessment tool (SWAT) model was calibrated and validated against observed streamﬂow using SWAT CUP. The Mann–Kendall trend test (MK) was used to assess climate trends. Meteorological drought (SPEI) and hydrological drought (SDI) were also investigated. Based on the ensemble mean of ﬁve global climate models (GCMs), projected increases in mean annual maximum temperature over the period 2015–2100 (compared with a 1983–2014 baseline) range from 1.16 to 1.73 ◦ C, while increases in minimum temperature range between 0.79 and 2.53 ◦ C. Increases in mean annual precipitation range from 1.8% at Addis Ababa to 45.5% over the Hombole area. High streamﬂow (Q 5 ) declines at all stations except Ginchi. Low ﬂows (Q 90 ) also decline with Q 90 equaling 0 m 3 s − 1 (i.e., 100% reduction) at some gauging stations (Akaki and Hombole) for individual GCMs. The SPEI conﬁrmed a signiﬁcant drought trend in the past, while the frequency and severity of drought will increase in the future. The basin experienced conditions that varied from modest dry periods to a very severe hydrological drought between 1986 and 2005. The projected SDI ranges from modestly dry to modestly wet conditions. Climate change in the basin would enhance seasonal variations in hydrological conditions. Both precipitation and streamﬂow will decline in the wet seasons and increase in the dry seasons. These changes are likely to have an impact on agricultural activities and other human demands for water resources throughout the basin and will require the implementation of appropriate mitigation measures.


Introduction
Scientific evidence is unequivocal that the earth's climate is changing. Due to increasing atmospheric concentrations of greenhouse gases, the average surface temperature of the earth has increased by 1 • C with a likely range of 0.8 to 1.2 • C since the preindustrial The current study further deepens knowledge of climate change impacts on the hydrology of the upstream part of the Awash Basin by investigating both the recent historical period  and the remainder of the current century (2015-2100). It includes trend analysis of observed hydrometeorological data, hydrological modelling, and the evaluation of meteorological and hydrological drought indices. The widely used hydrological model Soil and Water Assessment Tool (SWAT) developed by the US Department of Agriculture Research Service [23] is used to simulate the hydrology of the basin. The model is calibrated and validated using the SUFI2 algorithm in the SWAT CUP program [24]. Projections from five GCMs for the representative concentration pathway (RCP) 4.5 emission scenarios were selected based on data availability and coverage. These projections have been downscaled and bias-corrected [25] using historical climate data for the period 1983-2014. Uncertainties associated with GCMs and the downscaling techniques could result in significant differences in simulated hydrological conditions [26]. These uncertainties, particularly those related to inherent errors of the models, can be reduced by using a multimodel ensemble mean of the GCMs, and several studies have suggested using a multimodel ensemble to minimize the biases and uncertainties of future climate simulations [26][27][28]. Accordingly, this study uses both individual downscaled GCMs and the ensemble mean from all of the GCMs. The downscaled climate data were then used as input to the SWAT model to simulate impacts on river discharge. The novelty of the study lies in the combination of trend assessment for both past decades and future projections, evaluation of climate change impacts on hydrological extremes, and the implications for both meteorological and hydrological drought.

Study Area
The Upper Awash River Basin is located in the North West Rift Valley of Ethiopia and has an area of around 11,430 km 2 . The river rises on the high plateau near Ginchi town, west of Addis Ababa, at an altitude of more than 3000 m. The total length of the main course of the river is 336 km (Figure 1).  The climate of the basin varies from humid subtropical to semiarid. Based on data from meteorological stations for the period 1983 to 2014 obtained from the National Meteorological Agency of Ethiopia, the mean maximum annual temperature ranges from 23.3 to 28.0 • C, while the mean annual minimum temperature ranges from 9.5 to 13.3 • C. The basin receives most of its precipitation in two main seasons, spring and summer. Spring, which is locally named Bulg, is a short rainy season from March to May. Summer, known locally as Kiremt, is the main rain and crop-growing season and extends from June to September. Hombole and Tulu Bolo receive over 60% of their annual total rainfall in summer and 20% in spring. The other four stations receive more than 40% in summer and 30% in spring. The mean annual rainfall ranges from 600 mm at the Hombole meteorological station to 1094 mm at the Ginchi station.

Downscaling Climate Variables and Simulation Performance
This study used downscaled outputs from five GCMs (Table 1) for the RCP4.5 emission scenario. This is an intermediate emission scenario of long-term global emissions of greenhouse gases, short-lived species, and land use/land cover, which stabilizes radiative forcing at 4.5 W m −2 in the year 2100 [29]. RCP4.5 represents generally good agreement with the latest policies for lower greenhouse gas emissions by the global community [30]. The GCMs provide projections of future climate at resolutions of hundreds of kilometers and must be resampled at a higher spatial resolution to be used in studies of the impact of change, including those that focus on hydrological impacts using catchment modelling [31]. Therefore, the GCM projections were downscaled using a statistical downscaling model (SDSM) to project future temperature and precipitation across the Upper Awash River Basin.
Before downscaling the climate variables for each station, the simulation performance of each GCM was assessed to validate their accuracy in representing observed climate. For this purpose, we used observed and historical GCM daily data for the period 1983-2005. The mean relative error (MRE) (Equation (1)) was used to define the model simulation error, the correlation coefficient (CORR) (Equation (2)) identified the similarity or otherwise of observed and simulated climate data, and the Nash-Sutcliffe efficiency coefficient (NSE) (Equation (3)) was further used to evaluate the GCMs' capability in simulating climate variables. The smaller the value of the MRE is, the higher the GCM simulation capability is [32], while the closer the NSE value is to 1, the better is the model simulation performance [32].
where P i and O i are the values at the ith time step (daily) in the downscaled GCM simulation and observational time series, respectively, P and O are the mean value of the downscaled GCM simulation and observational sequence, and N is the total sample number [32]. An SDSM is an established decision support tool for assessing local climate change impacts using a robust statistical downscaling technique (see [33] for a detailed review of its conceptual basis). It is computationally inexpensive and thus can be easily applied to outputs from multiple GCM experiments to provide specific local climate information [34]. The tool is designed to identify local or regional climate information by determining a statistical model that relates large-scale climate variables (the predictors) to regional or local variables (the predictand) [31].
An SDSM was used to define statistical models using daily observed data (i.e., local climate data for a specific location for the predictand and larger-scale National Centers for Environmental Prediction (NCEP) data for the predictors). These models are then used with GCM-derived predictors to obtain daily weather data at a specific site for a future period [34]. This study used data from the six meteorological stations located within the basin (Ginchi, Alem Tena, Addis Ababa, Hombole, Mojo, and Tulu Bolo). The observed data covered the period 1983-2014, while the GCM simulated historical dataset covered the period 1983-2005 and projected data for the period 2006-2100. The GCM data (precipitation and maximum and minimum temperature) were downscaled to the station resolution. During downscaling, the linear scaling bias correction method [24], refs. [34][35][36] within an SDSM was employed to correct the GCM simulation systematic or random errors.

Climate Change and Impact Assessment
The Mann-Kendall trend test (MK) was employed to assess the climate change trends for both the historical period and future projections. The trends for precipitation, as well as the minimum and maximum temperature, were investigated for all stations. Furthermore, the Standardized Precipitation Evapotranspiration Index (SPEI) [37] was calculated at a 12-month temporal resolution, while the Streamflow Drought Index (SDI) [38] was evaluated to assess the impact of climate change on meteorological and hydrological drought. The SPEI is a new drought index used to detect the influence of global warming on meteorological drought, whereas the SDI is used to describe the severity of hydrological drought [39]. Hydrological drought refers to a significant reduction of water in the hydrological system, low streamflow in rivers, and abnormally low levels in lakes, reservoirs, and groundwater [38]. The trends for the SPEI and SDI were investigated for both the historical period and future projections at all stations. The procedure and description for the three indices are described as follows:

Mann-Kendall Trend Test (MK)
The MK is a nonparametric test that is commonly used to detect monotonically (increasing or decreasing) trends in climate and hydrological time series data [40]. The nonparametric MK test is suitable for trend detection because it is less affected by outliers in the dataset [41]. In addition, the MK test enables the identification of whether an established trend is statistically significant. The MK test statistic "S" is computed as follows: where x i and x j denote sequential data values of the time series of j and i (j > i), and n is the number of days in the data series. The sign function can be calculated as: When the sample size is greater than or equal to 10 (n ≥ 10), the S statistic is approximately a standard normal distribution with the mean equal to zero; then the variance is calculated with the following equation: where n is the number of data points, m is the number of tied groups (sample data with the same value, where there is zero variance between the compared values), and t i is the number of data values in the mth group. The Z-test is calculated with the following equation: where s denotes the variance and σ denotes the square root of the variance. Positive Z values indicate an increasing trend, while a negative value indicates a decreasing trend. The trend magnitude was calculated using Sen's slope (S i ): where x j and x i denote the data value at time j and i (j > i), respectively. The median of these N values of S i is denoted as Sen's estimator of the slope, which is expressed as: The sign of Q i shows whether the trend is increasing or decreasing. In this study, two specific significance levels, α = 0.05 and α = 0.01, are used (the corresponding threshold values of the MK value are ±1.96 and ±2.58, respectively). When the MK value is greater than 1.96 or less than −1.96, the changing trend is significant. When the MK value is greater than 2.58 or less than −2.58, the changing trend is extremely significant.

Standardized Precipitation Evapotranspiration Index (SPEI)
The SPEI uses the monthly difference between precipitation and potential evapotranspiration (PET), computed using the Penman-Monteith equation [42]. The water balance is computed as: where P and PET denote precipitation and potential evapotranspiration for the ith months, respectively. D i is a simple measure of the water surplus or deficit for the analyzed month. As regards the Standard Precipitation Index (SPI) [43], the probability distribution of the cumulative D i series was aggregated at different time scales. The log-logistic distribution was selected for standardizing the D series. The probability density function of a threeparameter log-logistic distributed variable is expressed as: where α, β, and γ are scale, shape, and origin parameters, respectively, for D values in the range (γ > D < ∞). The parameters of the log-logistic distribution can be obtained following [44]. The probability distribution function for standardizing the D series for all time scales, according to the log-logistic distribution, is determined as: With F(x), the SPEI can be obtained as the standardized values of F(x): where W = − ln(P) for P ≤ 0.5. Where P is the probability of exceeding a determined D value, P = 1 − F(D). If P > 0.5, then P is replaced by 1 − P, and the sign of the resultant SPEI is reversed. The constants are C 0 = 2.515517, C 1 = 0.802853, C 2 = 0.010328, d 1 = 1.432788, d 2 = 0.189269, and d 3 = 0.001308 [37]. The average value of the SPEI is 0, and the standard deviation is 1. The SPEI is a standardized variable, and it can therefore be compared with other SPEI values over time and space [37]. The SPEI classification is described in Table 2 below.  The SDI was calculated for the six gauging stations for both the historical period and the future projections. The historical SDI was calculated based on observed streamflow data and for future projections simulated by the SWAT model (see below). The SDI is based on the time series of monthly streamflow volumes (Q i,j ), where i denotes the hydrological year and j is the month within that hydrological year (from October to September). Streamflow volume is calculated using the following equation: Q i,j i = 1, 2, 3 . . . , 12 j = 1, 2 . . . 12 k = 1, 2, 3, 4 where V i,k is the cumulative streamflow volume for the i-th hydrological year and the k-th reference period, k = 1 for October-December, k = 2 for October-March, k = 3 for October-June, and k = 4 for October-September. The SDI is defined based on cumulative streamflow volumes V i,k for each reference period k of the i-th hydrological year as follows: where V k and S k are, respectively, the mean and the standard deviation of cumulative streamflow volumes of the reference period k as these are estimated over a long period. SWAT is a basin-scale, continuous-time model that operates at a daily time step and is designed to simulate the impact of management on water, sediment, and agricultural chemical yields in river catchments [45]. The hydrological component of the model simulates surface runoff, potential evapotranspiration, percolation, lateral subsurface flow, groundwater flow to streams from shallow aquifers, snowmelt, transmission losses from streams, and water storage and losses from ponds [23]. The model includes climate, nutrient cycling, soil temperature, sediment movement, crop growth, agricultural management, and pesticide dynamics [46], although most of these elements were not used in the current study given its focus on streamflow. The model is widely used (e.g., [47][48][49][50][51]), and previous applications include simulating the impacts of climate change and land management practices on water, sediment, and agricultural yields in large and complex watersheds with varying soils, land uses, and management conditions over long periods [52].
The input data for the SWAT model include climate time series and the location of meteorological stations, a digital elevation model (DEM), land use and soil distributions, and hydraulic characteristics. The data type used in the development of the SWAT model of the Upper Awash Basin and their sources are summarized in Table 3. SWAT divides a river basin into multiple subcatchments, which are then further subdivided into hydrologic response units (HRUs), which consist of homogeneous land use, management, and soil characteristics. In the case of Upper Awash, the model comprised a total of 7 sub-catchments and 93 HRUs. The model uses daily precipitation, maximum and minimum temperature, solar radiation, relative humidity, and wind speed from historical (observed) and downscaled GCM data from the station nearest to the centroid of each subcatchment. The surface runoff, potential evapotranspiration, percolation, lateral subsurface flow, groundwater flow to streams from shallow aquifers, and soil water were simulated for each HRU and then aggregated for the subcatchment using a weighted average [53]. These simulated variables for each subcatchment were then routed through the river system. Channel routing was simulated using the variable storage method. Surface runoff was simulated using a modified SCS curve number method, which is a function of local land use, soil type, and antecedent soil moisture condition [23,53].
SWAT calculates evaporation from soil and plants separately. The model simulates potential soil water evaporation as a function of PET and the leaf area index, which is the area of plant leaves relative to the soil surface area. Actual soil evaporation is estimated by using exponential functions of soil depth and water content. Plant water evaporation is simulated as a linear function of PET, leaf area index, and root depth. We employed Hargreaves methods to estimate PET [47].

Model Calibration and Validation
Calibration and validation of the SWAT hydrological model were undertaken using the sequential uncertainty fitting ver-2 (SUFI-2) algorithm in SWAT CUP [54]. In SUFI-2, all sources of uncertainties, such as uncertainty in driving variables, conceptual model, parameters, and measured data, are mapped to a set of parameter ranges [54]. The uncertainty is quantified at the 2.5% and 97.5% levels of the cumulative distribution of an output variable obtained through Latin hypercube sampling, and it is referred to as the 95% prediction uncertainty (95PPU) [55,56]. In SUFI-2, two indices are used to quantify the strength of the calibration and validation [47]. These are referred to as the "P-factor", which is the fraction of measured data bracketed by the 95PPU, and the "R-factor", defined as the ratio of the average width of the 95PPU band and the standard deviation of the measured variable. The P-factor varies from 0% to 100%, while the R-factor ranges between 0 and infinity. A P-factor of 1 and an R-factor of 0 indicate perfect model performance considering the uncertainty [57]. Furthermore, the algorithm allows the use of 10 different objective functions, including r 2 , Nash-Sutcliffe (NS), and mean square error (MSE). This study employed the P-factor and R-factor as well as the r 2 and NS objective functions.
A sensitivity analysis was undertaken using 16 SWAT input parameters with initial values collected from the literature [5,55,57,58]. The main purpose of the sensitivity analysis was to identify the key parameters that affect model performance during calibration and validation. We used similar initial parameter ranges for the six gauging stations as SWAT CUP allows for the calibration and validation of all stations simultaneously. SUFI-2 allows two types of sensitivity analysis: global sensitivity and one-at-a-time sensitivity. In a global sensitivity analysis, parameter sensitivities are determined using a multiple regression approach, which regresses the Latin hypercube-generated parameters against the objective function values. The t-stat values were used to assess the sensitivity of parameters [57] with larger absolute t-stat values indicating greater sensitivity. P-values are used to determine the significance of the sensitivity where the parameter becomes significant if the P-values are close to 0 [57]. The local sensitivity analysis determines the sensitivity of a variable to the changes in a single parameter if all other parameters are kept constant at their original value [54].
The sensitivity analysis results showed that runoff curve number was the most sensitive parameter, followed by average slope length, groundwater delay, and Manning's (roughness) value for the main channel, respectively (Table 4). After sensitivity analysis, calibration and validation were undertaken for each gauging station based on recorded monthly streamflow data. Accordingly, we used recorded streamflow data from 1993 to 2010 for Ginchi, 1985 to 2009 for Holeta, 1985 to 2004 for Asigori, 1990 to 2004 for Akaki, 1985 to 2008 for Melka Kuntrie, and 1983 to 2014 for Hombole. Two-thirds of the observed data were used for calibration, and the remaining one-third was for validation. During the calibration and validation process, the relative change method was used to parameterize, or regionalize, parameter values. This is more suitable than the alternative replacement method for those parameters that exhibit spatial variations, such as hydraulic conductivity and curve number.

GCM Simulation Performance
The simulation performance assessment showed that all of the selected GCMs have good simulation capacity over the Upper Awash River Basin. The MREs for GFDLE-CM3, GISS-ET-R, and MPI-ESM-LR are all less than 5% for precipitation and less than 0.55 for temperature. The other GCMs have values of MRE in the range of 7.8% to 42%. The CORR and NSE values for precipitation are all greater than 0.89 and 0.71, respectively. The values of these statistics for air temperatures are similarly high and indicate superior GCM performance than for precipitation [32]; they range between 0.94 and 0.99, and between 0.83 and 0.99, respectively (Appendix A).

Hydrological Model Calibration and Validation
Calibration and validation of the SWAT model were performed against recorded streamflow data at the six gauging stations ( Figure 2). Simulation performance for each gauging station and both the calibration and validation periods is described in Table 5. Performance was satisfactory at the Asigori, Melka Kuntrie, Hombole, and Holeta gauging stations. Performance at Ginchi was relatively poor for the calibration period but good for the validation period. In contrast, the model performed poorly for the Akaki gauging station in both the calibration and validation periods as indicated by the lower R 2 and NS values relative to the other gauging stations.

Changes in Air Temperature
The mean monthly minimum temperature (Tmin) and maximum temperature (Tmax) for the baseline and each GCM, as well as the ensemble mean based on all of the GCMs, are presented in Figures 3 and 4. Across the basin, Tmin for the baseline period ranges between 5.67 and 15.24 • C, while Tmax ranges from 16.70 to 31.06 • C. The mean monthly Tmin for the ensemble means ranges between 7.25 and 16.48 • C, while the corresponding ensemble mean Tmax ranges from 19.10 to 32.65 • C.
The mean annual change in Tmax ( • C) and Tmin ( • C) for the different GCMs relative to the baseline are presented in Table 6. The baseline (1983-2014) and projected (2015-2100) rate of change of annual Tmax and Tmin, which is defined by Sen's slope in the MK test, are also presented in Table 7. Both mean annual Tmax and Tmin are projected to increase relative to the baseline at all stations for all five GCMs. Based on the ensemble mean of the GCMs, the change in Tmax across the basin ranges from 1.16 to 1.73 • C, while the Tmin changes range from 0.79 to 2.53 • C. However, individual GCMs project increases in Tmax and Tmin of up to 1.9 and 3.47 • C, respectively (Table 6).   The mean monthly minimum temperature (Tmin) and maximum temperature (Tmax) for the baseline and each GCM, as well as the ensemble mean based on all of the GCMs, are presented in Figures 3 and 4. Across the basin, Tmin for the baseline period ranges between 5.67 and 15.24 °C, while Tmax ranges from 16.70 to 31.06 °C. The mean monthly Tmin for the ensemble means ranges between 7.25 and 16.48 °C, while the corresponding ensemble mean Tmax ranges from 19.10 to 32.65 °C.     The mean annual change in Tmax (°C) and Tmin (°C) for the different GCMs relative to the baseline are presented in Table 6. The baseline (1983-2014) and projected (2015-2100) rate of change of annual Tmax and Tmin, which is defined by Sen's slope in the MK test, are also presented in Table 7. Both mean annual Tmax and Tmin are projected to increase relative to the baseline at all stations for all five GCMs. Based on the ensemble mean of the GCMs, the change in Tmax across the basin ranges from 1.16 to 1.73 °C, while the Tmin changes range from 0.79 to 2.53 °C. However, individual GCMs project increases in Tmax and Tmin of up to 1.9 and 3.47 °C, respectively (Table 6).    Based on the MK trend tests, an extremely significant (α ≤ 0.001) increasing trend is indicated for both baseline and projected annual Tmax and Tmin ( Table 7). The rate of change of temperature for the baseline period ranges from 0.2 to 0.8 • C per decade for both Tmax and Tmin. Projected annual rates of temperature change, however, vary between GCMs. The rate of change of mean annual Tmax is 0.

Changes in Precipitation
Mean annual baseline precipitation for each station and changes (percent) relative to this baseline for the different GCMs and the ensemble mean are presented in Table 8. Increases in mean annual precipitation are indicated in most cases, although GISS E2 R projects decline at five of the six stations (the exception being Hombole). The largest overall increase in mean annual precipitation is projected for Hombole (a 60.5% increase for MPI ESM LR), while the largest decline is at Addis Ababa (14.9%, GISS E2 R).  trend, again, not significant, in autumn and winter at all stations. Monthly precipitation trends: In the majority of cases, there were nonsignificant (α ≤ 0.05) trends in most of the months at all stations ( Figure 5). The five GCMs' ensemble mean also projected a nonsignificant trend in most months and at most stations. The exception was April for Ginchi (significant decrease), October, and November for Hombole (significant increase).

Changes in Streamflow
Change in streamflow extremes: Table 10 presents mean discharge (Qmean) (m 3 s −1 ), discharges equaled or exceeded for 5% of the time (Q5) (m 3 s −1 ), and discharges equaled or exceeded for 90% for the time (Q90) (m 3 s −1 ) for each of the five gauging stations for the baseline as well as the changes (percent) from these baseline figures for each GCM and

Changes in Streamflow
Change in streamflow extremes: Table 10 presents mean discharge (Q mean ) (m 3 s −1 ), discharges equaled or exceeded for 5% of the time (Q 5 ) (m 3 s −1 ), and discharges equaled or exceeded for 90% for the time (Q 90 ) (m 3 s −1 ) for each of the five gauging stations for the baseline as well as the changes (percent) from these baseline figures for each GCM and the ensemble mean. Following the approach of [59], Q mean , Q 5 , and Q 90 were derived from monthly baseline and projected discharge time series data. The annual and seasonal streamflow trends for each gauging station for each of the GCMs and the ensemble mean are described in Table 11. Table 10. Q mean , Q 5 , and Q 90 discharge in m 3 s −1 for the baseline and change (percent) from the baseline for each GCM and ensemble mean (Minus sign (−) indicates declines relative to the baseline).  The projected Q mean at Ginchi and Holeta show an increase for all five GCMs and the ensemble mean. In contrast, Akaki Q mean shows a decrease for all GCMs except CNRM CM5. Asigori and Kuntrie, on the other hand, show a decrease for CMCC CMS and GISS E2 R and an increase for CNRM CM5, GFDL CM3, and MPI ESM LR. Q mean at Hombole shows an increase for all GCMs except CNRM CM5. The largest increase in Q mean is projected in the upper catchment at Ginchi (198.3% for MPI ESM LR), while the largest decrease (55.2% for CNRM CM5) is projected in the lower catchment at Hombole. The ensemble mean projects an increase at all stations except Akaki.
All GCMs project increases in streamflow during periods of high flow (i.e., Q 5 ) at Ginchi and Holeta and decreases at Asigori and Akaki. Q 5 declines at Kuntrie for CMCC CMS and increases for the remaining GCMs. Results for Hombole, on the other hand, show declines in high flows for CMCC CMS, CNRM CM5, and GFDL CM3 and increases for GISS E2 R and MPI ESM LR. Q 5 for the ensemble mean indicates declines during the high flow period at all stations except Ginchi. Low streamflow (i.e., Q 90 ) increases at Holeta and decreases at Akaki and Hombole. Q 90 at Kuntrie, on the other hand, increases for CMCC CMS, CNRM CM5, and MPI ESM LR and decreases for GFDL CM3 and GISS E2 R. Q 90 for the ensemble mean declines at Akaki and increases at Holeta, Kuntrie, and Hombole catchments.
Annual streamflow trends: The annual baseline streamflow trends (Table 11) reveal increasing discharge at Melka Kuntrie, Hombole, and Akaki, while decreasing trends are evident at Ginchi, Asigori, and Holeta. However, only Melka Kuntrie showed a statistically significant increasing trend (6.9 m 3 s −1 per decade). In common with projected precipitation and temperature, there is no overall agreement on the projected trends in both the magnitude and direction of change for streamflow across the GCMs (Table 11). Based on individual GCM projections, significant increasing trends in annual mean discharge were identified at Ginchi tively). Significant increasing trends (0.1 m 3 s −1 per decade) are projected at Asigori from September to January, while significant decreasing trends are projected for April and July (0.1 m 3 s −1 and 0.3 m 3 s −1 per decade, respectively). At both Akaki and Melka Kuntrie, the ensemble mean projects significant increasing trends for October and July (0.2 m 3 s −1 and 1.2 m 3 s −1 per decade), while at Hombole significant decreasing trends are projected for April (1.2 m 3 s −1 per decade) and July (3.3 m 3 s −1 per decade).

Climate Drought Index
Using the historical climate data (1984-2020) and the five GCM ensemble mean projections (2021-2100), the historical and potential future climatic droughts were analyzed for each meteorological station using the SPEI at a time scale of 12 months (Figure 7). Accordingly, the SPEI for the historical period shows that the years 1996,1997,2007,2008,2013,2014,2016, and 2018 were very dry, while 1999 was a very wet year. Figure 8 clearly shows an increasing trend towards meteorological drought in recent decades (2006-2020) compared with the beginning of the period (1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994). The SPEI trend analysis demonstrates a decrease towards a negative value at all stations. This trend was statistically significant at α ≤ 0.01, confirming the increasing tendency towards meteorological drought.
The projected SPEI for the period 2021-2100 indicates an increase in the frequency and severity of drought especially for the period 2045-2072 compared with 2021-2047. The increasing drought trend is significant (α ≤ 0.01) at all stations.

Climate Drought Index
Using the historical climate data (1984-2020) and the five GCM ensemble mean projections (2021-2100), the historical and potential future climatic droughts were analyzed for each meteorological station using the SPEI at a time scale of 12 months (Figure 7). Accordingly, the SPEI for the historical period shows that the years 1996,1997,2007,2008,2013,2014,2016, and 2018 were very dry, while 1999 was a very wet year. Figure 8 clearly shows an increasing trend towards meteorological drought in recent decades (2006-2020) compared with the beginning of the period (1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994). The SPEI trend analysis demonstrates a decrease towards a negative value at all stations. This trend was statistically significant at α ≤ 0.01, confirming the increasing tendency towards meteorological drought.
The projected SPEI for the period 2021-2100 indicates an increase in the frequency and severity of drought especially for the period 2045-2072 compared with 2021-2047. The increasing drought trend is significant (α ≤ 0.01) at all stations.

Hydrological Drought Index
The historical SDI is dominated by negative values, indicating dry years at most of the stations, the exception being Akaki (Figure 8). The basin experienced a range of conditions from modestly dry to very severe hydrological drought between 1986 and 2005. The MK trend test confirmed a significant increasing trend towards negative values of SDI at Ginchi, Asigori, Holeta, and Melka Kuntrie, supporting the assertation of increasing hydrological drought over the past three decades. However, a significant increasing trend towards positive values (indicating wetter conditions) was identified for Akaki, while no significant trend was established for Hombole.
The projected SDI, on the other hand, indicates conditions that range from modestly dry to modestly wet at all stations. The MK trend analysis confirmed an increasing trend towards positive values (i.e., wetter years) at all stations. Trends were statistically significant at Asigori, Ginchi, and Akaki.

Hydrological Drought Index
The historical SDI is dominated by negative values, indicating dry years at most of the stations, the exception being Akaki (Figure 8). The basin experienced a range of conditions from modestly dry to very severe hydrological drought between 1986 and 2005. The MK trend test confirmed a significant increasing trend towards negative values of SDI at Ginchi, Asigori, Holeta, and Melka Kuntrie, supporting the assertation of increasing hydrological drought over the past three decades. However, a significant increasing trend towards positive values (indicating wetter conditions) was identified for Akaki, while no significant trend was established for Hombole.
The projected SDI, on the other hand, indicates conditions that range from modestly dry to modestly wet at all stations. The MK trend analysis confirmed an increasing trend towards positive values (i.e., wetter years) at all stations. Trends were statistically significant at Asigori, Ginchi, and Akaki. trend towards positive values (indicating wetter conditions) was identified for Akaki, while no significant trend was established for Hombole.
The projected SDI, on the other hand, indicates conditions that range from modestly dry to modestly wet at all stations. The MK trend analysis confirmed an increasing trend towards positive values (i.e., wetter years) at all stations. Trends were statistically significant at Asigori, Ginchi, and Akaki.

The Impacts of Climate Change on Precipitation
This study suggests that some changes in the climate of the Upper Awash River Basin have occurred over the period beginning in 1984 and extending to 2014. Declines in precipitation in spring and summer and increases in autumn and winter were identified over this baseline period, although these trends were not statistically significant. Results show that at some stations precipitation decreased by up to 22.4 and 8.8 mm per decade in spring and summer, respectively, between 1983 and 2014. This study, therefore, supports earlier investigations regarding the impacts of climate change on precipitation and their implication for Ethiopia [60][61][62][63][64]. For instance, a study in the southern part of the country revealed that precipitation in spring and summer declined between 1972 and 2011 [65]. In another study in the northeastern part of the Blue Nile basin, a nonsignificant decreasing trend was observed in spring between 1981 and 2010 [66]. A further study also showed a 4. Discussion

The Impacts of Climate Change on Precipitation
This study suggests that some changes in the climate of the Upper Awash River Basin have occurred over the period beginning in 1984 and extending to 2014. Declines in precipitation in spring and summer and increases in autumn and winter were identified over this baseline period, although these trends were not statistically significant. Results show that at some stations precipitation decreased by up to 22.4 and 8.8 mm per decade in spring and summer, respectively, between 1983 and 2014. This study, therefore, supports earlier investigations regarding the impacts of climate change on precipitation and their implication for Ethiopia [60][61][62][63][64]. For instance, a study in the southern part of the country revealed that precipitation in spring and summer declined between 1972 and 2011 [65]. In another study in the northeastern part of the Blue Nile basin, a nonsignificant decreasing trend was observed in spring between 1981 and 2010 [66]. A further study also showed a decreasing trend in precipitation in spring throughout the country between 1980 and 2010 [67].
The IPCC third assessment report [2] stated that because of the diversity of climatic conditions and limited meteorological stations, projected precipitation in Africa, particularly at regional and local scales, is not well defined. Under intermediate warming scenarios, one study showed an increase in projected rainfall during the dry season (December to February) and a decrease in the rainy season (June to August) over parts of East Africa, including the current study area [2]. This study showed that for the ensemble mean, projected future precipitation, especially in spring and summer and as a result annually, showed a decreasing trend at most stations. The magnitude and directions of change do, however, vary between different GCMs, stations, and seasons, echoing the variability identified in earlier studies employing a number (in most cases, fewer than the current study) of climate models [20,21,67].
Declining precipitation will have a significant implication for agricultural productivities. As one of the major agricultural regions of Ethiopia, spring and summer are overwhelmingly important in the Upper Awash River Basin because more than 95% of crop production is produced at this time of year [68]. Decreases in precipitation together with the extremely variable nature of rainfall, and increasing temperature, are likely to have a significant impact on agricultural production as well as the water sector of the basin. Furthermore, the economic dependency of both rural and urban communities on agriculture, combined with the dominant rain-fed nature of the farming system, makes the area particularly vulnerable to climate variability and change. It has been suggested that climate change would affect Ethiopia's GDP growth by 0.5-2.5% per year due to declining agricultural productivity [69]. Another study showed that, due to rainfall variability, about USD 2 billion will be lost in the agricultural sector over future decades [70].

The Implication of Climate Change on Streamflow
Streamflow is strongly influenced by changes in precipitation, temperature, land use, and other factors, such as the withdrawal of water for human use. Change in streamflow can affect the amount of water available for irrigation, water supply, and hydropower generation. Seasonal and spatial variability of streamflow and their changes into the future, therefore, have the potential to exert significant impacts on the economy and local development. The area downstream of the Upper Awash River Basin is relatively developed for irrigation. These irrigation developments include the Wonji and Metehara sugar plantations; upper, middle, and lower Awash state farms developed for fruit and vegetable production; and cotton plantations. Such schemes have flourished following the construction of the Koka Dam, which defines the outlet of the current study area [71]. Streamflow changes in the upper catchment will have an impact on reservoir inflows, and they, therefore, have the potential for large impacts on downstream projects.
This study suggests that discharges during the high flow period (as indicated by Q 5 ) decline by up to 61.9% compared with the baseline at some gauging stations with the largest changes in the opposite direction indicating increases of 84.3%. The ensemble mean, however, projects a decline in these flows at five of the stations (mean of those declining = −40.8%). Conversely, mean flows for the ensemble mean are projected to increase at five stations (declines at the one remaining station), although there is a notable decline in the magnitude of the increases in a downstream direction (just 1.4% at the lowest station). Across the different GCMs, there is considerable variability in the magnitude of changes in mean flow, including both increases and decreases, which concur with those results of other multi-GCM studies of change over the basin [22]. The low flow period is projected to experience both declines and increases in discharge with Q 90 reducing by 100% (i.e., flows cease) at Akaki and Hombole for individual GCM projections (Table 10). The magnitude of changes in percentage terms is much higher for low flows than both mean and high flows. These larger percentage changes in low flows agree with an earlier study undertaken within the Upper Awash Basin using a number of GCMs, including some used in the current study [22]. It also supports a recent global analysis showing greater risks of change in low flow conditions compared with high flows [72]. However, for the ensemble mean there are more declines in high flows than in low flows across the different gauging stations. Reductions in streamflow during the high flow period results from precipitation reductions in spring and summer, the main rainy seasons in the basin. A similar study in the Blue Nile basin identified declines in low streamflow (Q 95 ) of 18.1% between 1982 and 2000, while high streamflow (Q5) increased by 7.6% over the same period [73].
The trend analysis also confirmed that streamflow in the upper catchment (Ginchi, Asigori, and Holeta gauging stations) experienced a decline over the period 1985-2010, while it increased in the lower catchment (Akaki, Melka Kuntrie, and Hombole) in most months and seasons. The projected streamflow, on the other hand, showed no overall agreement in terms of both magnitude and direction of change between different GCMs. As an example, GISS E2 R showed a dominant decreasing trend, while GFDL CM3 projected an increasing trend in most months at Ginchi, Asigori, Hombole, and Kuntrie. CMCC CMS projected a nonsignificant increasing trend, while for MPI ESMLR there was a significant increasing trend in most of the months at all stations. The ensemble mean showed an increasing monthly trend from September to February and a decreasing trend in April and July. The seasonal results showed a decreasing trend in spring and summer, while an increasing trend was identified for autumn and winter. At an annual level, results indicated a decreasing trend at Holeta and Hombole and an increasing trend at the remaining stations.

The Implication of Climate Change on Meteorological and Hydrological Drought
Changes in temperature and precipitation can exacerbate extreme events, such as meteorological and hydrological droughts. Drought is a common phenomenon that occurs frequently in many parts of Ethiopia, including the study area. The frequency of drought in the country has increased, and devastating droughts have occurred in the recent decade [74,75]. The historical meteorological drought index used in the current study indicated an extremely significant increasing trend across the basin (Figure 8). The projected drought index also suggested that the basin will experience more severe droughts over the current century. Many studies have stated that droughts in Ethiopia are El Niño induced and that this climatic phenomenon is responsible for extreme dryness and major rainfall failures in many parts of the country [21,64,71,72]. Historical droughts have been associated with severe water and food shortages, widespread hunger, and diseases that destroyed the livelihood of millions of people across a large part of Ethiopia [74]. With the projected increase in drought incidence and severity, these extreme conditions are likely to remain significant concerns in the country. The hydrological drought index, on the other hand, showed significant dry years at Asigori and Holeta and very wet years at Ginchi, Akaki, and Melka Kuntrie between 1985 and 2020. In contrast, the projected SDI indicated that more wet years are expected than dry years.

Conclusions
This study demonstrated that over the baseline period (1983-2014), the Upper Awash River Basin experienced increases in mean annual maximum temperature (Tmax) in the ranges of 0.60 to 1.90 • C, while mean annual minimum temperature (Tmin) increased between 0.63 and 3.47 • C. The rate of change of temperature for the baseline period ranges from 0.2 to 0.8 • C per decade for both Tmax and Tmin. The projected annual rate of change in temperature over the remainder of the current century varies in magnitude between the five GCMs (between 0.1 and 0.4 • C per decade).
The trend analysis for precipitation shows a decreasing trend over the baseline period for most of the stations. Mean annual precipitation is projected to increase at all stations for all GCMs except GISS ET R. The overall range of change across all GCMs and stations is −14.9% to 60.0%. There is notable variation in the magnitude and direction of change in monthly, seasonal, and annual precipitation across the different GCMs. In the majority of cases, there were nonsignificant (α ≤ 0.05) trends for most months at all stations. The ensemble mean based on the five GCMs also projected a nonsignificant trend in most months and at most stations. However, both baseline and projected precipitation showed a decrease in spring and summer and increases in autumn and winter.
Based on projections using the ensemble mean, streamflow declines during the high flow period at all stations except Ginchi. On the other hand, low flows decrease at Akaki and increase at Holeta, Kuntrie, and Hombole. The mean annual streamflow trend shows an increase at all stations except Akaki. Individual GCMs, however, project changes that differ in terms of both magnitude and direction. Changes range from a decrease of 55% (Hombole for CNRM CM5) and an increase of 198% (Ginchi for MPI ESM LR). Streamflow during wet and dry periods also showed an increase at Ginchi, Holeta, and Kuntrie and a decrease at Asigori, Holeta, and Hombole. The baseline and projected mean monthly, seasonal, and annual streamflow trends indicated a nonsignificant decrease at Ginchi, Asigori, and Holeta and an increase at Akaki, Hombole, and Melka Kuntrie. There was no agreement over the projected streamflow trend in terms of both the magnitude and direction across the different GCMs. The GISS E2 R GCM projected a decreasing trend, while other GCMs projected an increasing trend in most months at all stations. Statistically, the CMCC CMS, CNRM CM5, and GISS E2 R GCMs showed a nonsignificant trend, while the GFDL CM3 and MPI ESMLR GCMs indicated a significant increasing trend at all stations. The seasonal ensemble mean streamflow, however, showed a decreasing trend in spring and summer at all stations except Asigori. The ensemble mean also showed a decrease at Holeta and Hombole and an increase at the remaining stations.
The SPEI confirmed that the basin experienced significant droughts in the past, and the frequency and severity of drought will increase in the future. The historical SDI, on the other hand, showed hydrological drought at Asigori and Holeta and very wet years at the remaining stations. The trend of the projected SDI showed more wet years than dry years.
Climate change in the basin would create seasonal variation in precipitation and streamflow. Both precipitation and streamflow will decline in the wet seasons and increase in the dry seasons. These changes are likely to have an impact on the agricultural system and water resource use within the basin. As such, there will be an urgent need to develop mitigation approaches in response to changing hydrometeorological conditions. Agriculture. The funds are administered by the International Fund for Agricultural Development (IFAD), Rome, Italy, while the project is implemented by Alliance Bioversity-CIAT. The initial stages of the research were financially supported by the UK Department for International Development (DfID) and the Climate Impact Research Capacity and Leadership Enhancement (CIRCLE) program (grant number 201871). The final part of the study was financially supported through a grant from the World Bank "Accelerating Impacts of CGIAR Climate Research for Africa (AICCRA)" ESA regional project (Grant number D7540) signed between CIAT and the International Development Association (IDA).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to confidentiality requirements of the funding organisations for the project.

Acknowledgments:
The authors would like to thank the Ethiopian National Meteorological Agency; Ministry of Irrigation, Energy, and Water Resource of Ethiopia; and Global Runoff Data Centre for providing the required climate and streamflow data.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.