Watershed Hydrological Response to Combined Land Use/Land Cover and Climate Change in Highland Ethiopia: Finchaa Catchment

Land use/land cover (LULC) and climate change affect the availability of water resources by altering the magnitude of surface runoff, aquifer recharge, and river flows. The evaluation helps to identify the level of water resources exposure to the changes that could help to plan for potential adaptive capacity. In this research, Cellular Automata (CA)-Markov in IDRISI software was used to predict the future LULC scenarios and the ensemble mean of four regional climate models (RCMs) in the coordinated regional climate downscaling experiment (CORDEX)-Africa was used for the future climate scenarios. Distribution mapping was used to bias correct the RCMs outputs, with respect to the observed precipitation and temperature. Then, the Soil and Water Assessment Tool (SWAT) model was used to evaluate the watershed hydrological responses of the catchment under separate, and combined, LULC and climate change. The result shows the ensemble mean of the four RCMs reported precipitation decline and increase in future temperature under both representative concentration pathways (RCP4.5 and RCP8.5). The increases in both maximum and minimum temperatures are higher for higher emission scenarios showing that RCP8.5 projection is warmer than RCP4.5. The changes in LULC brings an increase in surface runoff and water yield and a decline in groundwater, while the projected climate change shows a decrease in surface runoff, groundwater and water yield. The combined study of LULC and climate change shows that the effect of the combined scenario is similar to that of climate change only scenario. The overall decline of annual flow is due to the decline in the seasonal flows under combined scenarios. This could bring the reduced availability of water for crop production, which will be a chronic issue of subsistence agriculture. The possibility of surface water and groundwater reduction could also affect the availability of water resources in the catchment and further aggravate water stress in the downstream. The highly rising demands of water, owing to socio-economic progress, population growth and high demand for irrigation water downstream, in addition to the variability temperature and evaporation demands, amplify prolonged water scarcity. Consequently, strong land-use planning and climate-resilient water management policies will be indispensable to manage the risks.


Introduction
Our world is functioning at a mixture of complex systems interacting at a variety of spatial and temporal scales in which both the natural and human factors are thoroughly intertwined. The strong relationship between land use, water and climate from these systems has contributed to prominent concern to the basic planet characteristics and process. The productivity of land, biodiversity, aridity

Study Area
The study was conducted in the upper Blue Nile Basin, Finchaa sub-basin in Oromia Regional State, Ethiopia. Finchaa sub-basin lies in 9°10' to 10°00' North latitude and 37°00' and 37°40' East longitude covering 3781 km 2 areas. The catchment is characterized by high topographic relief, with elevation ranging from 851 to 3213 m above sea level. The area has large upstream water potential sites, intensive irrigable downstream lands, and high hydropower potential. The sub-basin contains three watersheds; Fincha, Amerti and Neshe [2]. The detail study area description is shown in Figure 1. The annual rainfall of Finchaa catchment ranges from 1367 mm to 1842 mm, with the lower rainfall occurring in the northern lowlands, and the higher rainfall is greater than 1500 mm occurring in the southern and western highlands of the sub-basin. June to September is the main rainy season of the catchment, with an average annual rainfall of 1604 mm and a peak occurring between July to August [54]. The mean monthly temperature of the catchment varies between 15.50 °C to 18.62 °C.
Agro ecologically, Finchaa sub-basin is characterized by tepid to cool sub-humid mid highlands in the north-western, moist mid highlands in the south-eastern, and hot to warm moist lowlands in the north-eastern parts of the catchment.

Data
A 30 × 30 m resolution of the Digital Elevation Model (DEM) obtained from the United States Geological Survey (USGS) at https://earthexplorer.usgs.gov/ was used to delineate the watershed of the study area. Moreover, it is used for slope classification, which was the basis for the Hydrological Response Unit (HRU) generation. The areas of different land slope classes as shown in Figure 2 shows that the majority of the land is 31.57%, and is found with a slope greater than 15%, 24.13%, are present between 0% and 4%, 21.77% between 4% and 8% and the remaining 22.53% between 8% and 15%, showing a variation of higher topographic in Finchaa catchment. The annual rainfall of Finchaa catchment ranges from 1367 mm to 1842 mm, with the lower rainfall occurring in the northern lowlands, and the higher rainfall is greater than 1500 mm occurring in the southern and western highlands of the sub-basin. June to September is the main rainy season of the catchment, with an average annual rainfall of 1604 mm and a peak occurring between July to August [54]. The mean monthly temperature of the catchment varies between 15.50 • C to 18.62 • C.
Agro ecologically, Finchaa sub-basin is characterized by tepid to cool sub-humid mid highlands in the north-western, moist mid highlands in the south-eastern, and hot to warm moist lowlands in the north-eastern parts of the catchment.

Data
A 30 × 30 m resolution of the Digital Elevation Model (DEM) obtained from the United States Geological Survey (USGS) at https://earthexplorer.usgs.gov/ was used to delineate the watershed of the study area. Moreover, it is used for slope classification, which was the basis for the Hydrological Response Unit (HRU) generation. The areas of different land slope classes as shown in Figure 2 shows that the majority of the land is 31.57%, and is found with a slope greater than 15%, 24.13%, are present between 0% and 4%, 21.77% between 4% and 8% and the remaining 22.53% between 8% and 15%, showing a variation of higher topographic in Finchaa catchment. The SWAT model requires a basic physic-chemical property of the soil types. In this study, the major soil types were collected from the Ministry of Water, Irrigation and Electricity (MoWIE), a 250 m resolution of soil grids accessed at https://soilgrids.org/#!/?layer=TAXNWRB_250 m&vector=1, the world digital soil map at http://www.fao.org/geonetwork/srv/en/metadata.show?id=14116 and FAO soil [55]. The soil classification used in this study is based on the FAO classification system and customized in the way the SWAT model requires, as shown in Figure 3b. Ten different soil types were identified, as shown by Figure 3a. The most dominant soil type of the catchment is Haplic Alisols covering 42.8% areas of the land. The higher areas of the valley are dominated by haplic Arenosols (7.8%), valley areas are dominantly covered with Eutric Cambisols (13.63%), the upstream reach of Fincha reservoir is covered with marsh (9.1%) and Rhodic Nitosols (13.18%) are found around Nashe reservoir. The daily climate data required by the SWAT model was obtained from the National Meteorological Service Agency of Ethiopia. Based on the record length and the quality of data (less missed data values), only five stations located inside the study catchment were used in this study. The stations that only record rainfall, and stations with many missing values, were not used. Further, for selected stations, the missing data values have been filled with the Inverse Distance Weighting (IDW) method. The consistency of the data was also checked by a double mass curve and found to be consistent. Daily streamflow data for the years 1988-2008 required for SWAT model evaluation was obtained from the MoWIE, hydrology department. The SWAT model requires a basic physic-chemical property of the soil types. In this study, the major soil types were collected from the Ministry of Water, Irrigation and Electricity (MoWIE), a 250 m resolution of soil grids accessed at https://soilgrids.org/#!/?layer=TAXNWRB_250m&vector=1, the world digital soil map at http://www.fao.org/geonetwork/srv/en/metadata.show?id=14116 and FAO soil [55]. The soil classification used in this study is based on the FAO classification system and customized in the way the SWAT model requires, as shown in Figure 3b. Ten different soil types were identified, as shown by Figure 3a. The most dominant soil type of the catchment is Haplic Alisols covering 42.8% areas of the land. The higher areas of the valley are dominated by haplic Arenosols (7.8%), valley areas are dominantly covered with Eutric Cambisols (13.63%), the upstream reach of Fincha reservoir is covered with marsh (9.1%) and Rhodic Nitosols (13.18%) are found around Nashe reservoir. The SWAT model requires a basic physic-chemical property of the soil types. In this study, the major soil types were collected from the Ministry of Water, Irrigation and Electricity (MoWIE), a 250 m resolution of soil grids accessed at https://soilgrids.org/#!/?layer=TAXNWRB_250 m&vector=1, the world digital soil map at http://www.fao.org/geonetwork/srv/en/metadata.show?id=14116 and FAO soil [55]. The soil classification used in this study is based on the FAO classification system and customized in the way the SWAT model requires, as shown in Figure 3b. Ten different soil types were identified, as shown by Figure 3a. The most dominant soil type of the catchment is Haplic Alisols covering 42.8% areas of the land. The higher areas of the valley are dominated by haplic Arenosols (7.8%), valley areas are dominantly covered with Eutric Cambisols (13.63%), the upstream reach of Fincha reservoir is covered with marsh (9.1%) and Rhodic Nitosols (13.18%) are found around Nashe reservoir. The daily climate data required by the SWAT model was obtained from the National Meteorological Service Agency of Ethiopia. Based on the record length and the quality of data (less missed data values), only five stations located inside the study catchment were used in this study. The stations that only record rainfall, and stations with many missing values, were not used. Further, for selected stations, the missing data values have been filled with the Inverse Distance Weighting (IDW) method. The consistency of the data was also checked by a double mass curve and found to be consistent. Daily streamflow data for the years 1988-2008 required for SWAT model evaluation was obtained from the MoWIE, hydrology department. The daily climate data required by the SWAT model was obtained from the National Meteorological Service Agency of Ethiopia. Based on the record length and the quality of data (less missed data values), only five stations located inside the study catchment were used in this study. The stations that only record rainfall, and stations with many missing values, were not used. Further, for selected stations, the missing data values have been filled with the Inverse Distance Weighting (IDW) method. The consistency of the data was also checked by a double mass curve and found to be consistent. Daily streamflow data for the years 1988-2008 required for SWAT model evaluation was obtained from the MoWIE, hydrology department.

Impact Assessment Framework
The hydrological responses of Finchaa catchment to the impacts of LULCC and climate change were investigated for two future periods of 2021-2050 and 2051-2080. The baseline period consists of the years 1986-2015. Three simulation scenarios were developed for the independent and combined effects of the LULCC and climate change. The first scenario (S1) considers only LULCC, the second scenario (S2) considers only climate change and the third scenario (S3) considers both LULCC and climate change (the combination of the two). Details of the scenarios are depicted in Table 1.  [20,56]. The approach of CA-Markov in IDRISI involves two techniques: Markov chain analysis and CA. The procedures of projecting LULCCs with CA-Markov outlined by Pan et al. [22] in the framework of IDRISI software was used in this study. The procedure involves three steps: calculating the transition matrix, creating suitability maps and prediction of the LULC map. The LULC transition matrix was developed by the Markov chain model in IDRISI. Then, the suitability maps of the potential driving factors of the LULC translations are developed for each LULC classes with the multi-criteria evaluation (MCE). Factors such as road proximity, water body proximity, elevation, slope, and urban and built areas were used as criteria to generate suitability maps. Finally, the prediction of the LULC maps are done based on the transition matrix and suitability maps.
For the LULCC analysis, three maps of LULC map for 1987, 2002 and 2017 developed by Dibaba et al. [2] were used. Dibaba et al. [2] used Landsat images from the thematic mapper (TM) for 1987 and 2002, and Landsat 8 operational land imager (OLI) for 2017 to evaluate the spatio-temporal LULCC in Finchaa catchment. The land use maps of 2002 and 2017 were taken as a baseline map and CA-Markov in IDRISI software was used to predict future land use maps of LULC2036 and LULC2055.
Any model prediction requires model calibration and validation. Then, the usefulness of the model output is dependent on the validation results of the model. The performance of the CA-Markov model in predicting the land use maps was assessed using Kappa coefficients calculated by Equation (1) [22,57]. In this study, the LULC maps of 1987 and 2002 were used to simulate the 2017 LULC map. Then, the simulated LULC map of 2017 was compared with the observed 2017 LULC map using the Kappa index: In the above equation, Po is the proportion of correctly simulated cells; Pc is the expected proportion correction by chance between the observed and simulated map. If Kappa ≤ 0.5 shows rare agreement, 0.5 ≤ Kappa ≤ 0.75 shows a medium level of agreement, 0.75 ≤ Kappa ≤ 1 shows a high level of agreement and Kappa = 1 for perfect agreement.

Climate Change Scenarios
After evaluating the performance of six Regional Climate Models (RCMs) in Coordinated Regional Climate Downscaling Experiment (CORDEX) Africa, four RCMs showing a better performance were selected. Detail analysis of the CORDEX evaluation is presented in Dibaba et al. [54]. The study by Dibaba et al. [54] on how well the individual RCMs simulate precipitation and temperature of two Ethiopian highland catchments, including Finchaa, proved that ensemble mean outperforms most of the individual RCMs. A similar study by Gadissa et al. [40] in the central Rift valley basin showed that the ensemble mean of five RCMs have shown better performance than individual RCMs. Using multi-RCMs for climate modelling helps to minimize uncertainties, compared to the use of a single RCMs. However, despite the reasonable performance of the RCMs including the ensemble mean, all RCMs have exhibited noticeable biases. Therefore, corrections to the biases should be done before the climate change impact study.
The RCMs used in this study were: CCLM4-8, HIRHAM5, RACMO22T and RCA4, fed by data of the two representative concentration pathways (RCP) scenarios high emission scenario (RCP8.5) and mid-range mitigation emission (RCP4.5). The years 2021-2050 and 2051-2080 served as future scenario periods and the period from 1986 to 2015 as a historical baseline to evaluate the climate changes. Other weather variables, including solar radiation, relative humidity and wind speed in baseline/historical period were considered in the future scenarios without making any change as the changes in these variables may not have a significant impact in modelling the climate change scenarios on local hydrology [40].

Bias Correction
The climate model data for hydrological modelling (CMhyd) [58], obtained from https://swat.tamu. edu/software/, was used to process the precipitation and temperature bias correction. Teutschbein and Seibert [59] have provided a full review of the bias correction techniques. According to Teutschbein and Seibert [59], all the bias correction techniques have improved the simulation of precipitation and temperature. However, there are differences between the correction methods in the daily precipitation series, standard deviations, and percentiles. Based on mean absolute error ranking, distribution mapping was ranked both, for temperature and precipitation corrections.
The distribution mapping (DM) uses a transfer function to adjust the cumulative distribution of estimated data to the cumulative distribution of rain gauges [59]. Luo et al. [60] applied seven precipitation bias correction methods and five methods for temperature. The results showed that distribution mapping reproduces precipitation and temperature very well. A similar study by Zhang et al. [61] compared five bias correction methods, using CMhyd, and found distribution mapping performed best for climate change impact study on streamflow dynamics of two rivers in Northern lake Erie basin, Canada. The study by Geleta and Gobosho [52] used CMhyd for extraction of CORDEX-NetCDF, and bias correction of minimum and maximum temperature to predict climate change-induced temperature changes in Finchaa catchment. The study showed that distribution mapping was better in improving the simulation. Owing to all the above findings, distribution mapping served as a basis for the precipitation and temperature corrections before using CORDEX-RCMs outputs for a climate impact study.

Hydrological Model
Soil and Water Assessment Tool (SWAT) is a physically-based semi-distributed model that operates on a continues time scale [62]. SWAT model operates on a daily time step and predicts the impact of land management in large complex watersheds with varying soils, land use and management conditions. Major model components include DEM, weather, hydrology, soil and properties and land management [63]. In SWAT, a watershed is divided into multiple sub-watersheds, which are then further subdivided into Hydrologic Response Units (HRUs) that comprise homogeneous land use, slope and soil characteristics. The hydrological components in the model are based on the water balance equation [63] given in Equation (1): where, SWt is the final soil water content(mm), SWo is the initial water content (mm), t is the time (days), Rday is the amount of precipitation on day i (mm), Qsurf is the amount of surface runoff on day i (mm), Ea is the amount of evapotranspiration on day i (mm), Wseep is the amount of water entering the vadose zone from the soil profile on day i (mm) and Qgw is the amount of return flow on day i (mm). In the SWAT model, the simulation of the hydrological process begins with watershed delineation and generating streamflow networks. With a 30 m resolution DEM, Finchaa watershed is delineated into 25 sub-basins and a multiple HRU was defined with classes of 10% land use, 20% soil and 10% slope. Accordingly, 231 HRUs were created. Then, with the input of weather data from five stations (Fincha, Hareto, Kombolcha, Nashe and Shambu), the SWAT model setup was ready for the first simulation, which was used for model evaluation. Shambu station, consisting of all the climate variables, was considered a weather generator.
Before using the output of the SWAT simulation for analysis, the performance of the model was evaluated for the catchment. Although distributed hydrological models can relate spatial changes of LULCC to the hydrological process simulation, it is difficult to transform many calibrated parameters into time-variant conditions of the future, especially if empirical parameters are used [22]. Therefore, parameters of the SWAT model are calibrated and validated by historical data and were assumed to have the extrapolative ability under future scenarios.

Sensitivity Analysis, Calibration and Validation
Owing to a large number of flow parameters in SWAT, identifying the most sensitive parameters is necessary to improve the calibration of the hydrological model. Through the sensitivity analysis, the most sensitive parameters that strongly influence the flow process will be identified. The Sequential Uncertainty Fitting (SUFI-2) embedded in the SWAT-CUP (Calibration and Uncertainty Program) was used to achieve the sensitivity analysis, calibration and validation [64].
The fitness of the model simulation with the observed streamflow was expressed by statistics like coefficients of determination (R 2 ), Nash-Sutcliffe efficiency (NSE), percent bias (PBIAS) and the ratio of the root-mean-square error to the standard deviation of measured data (RSR). The model performance ratings were based on the statistics recommended by Moriasi et al. [65] and Ayele et al. [67] R 2 varies between 0 and 1, where higher value shows less error. NSE ranges from negative infinity to 1(best). PBIAS close to 0 shows best the simulation, a negative value indicates overestimation and a positive value indicates under simulation of the model. RSR varies from zero to a large positive number. The lower RSR shows a better simulation of the model. These statistics are calculated using Equations (2) to (5): Water 2020, 12, 1801 9 of 25 In the above equations, Qm is the measured discharge, Qs is the simulated discharge, Qm is the average measured discharge, and Qs is the average simulated discharge.
The propagation of uncertainties in model outputs in SUFI-2, expressed as the 95% probability distribution, calculated by the 2.5% and 97.5% levels of the cumulative distributions of output variables, is considered as 95PPU [64]. P-factor and R-factor statistics are considered quantifying, the fit between the result expressed as 95PPU and observation. P-factor, the percentage of observations covered by the 95PPU varies from 0 to 1 with the ideal value of 1 while for R-factor, the thickness of the 95PPU optimal value is around 1.

CA-Markov Model Performance and Land Use/Land Cover Change Projections
The observed LULC1987 and LULC2002 were used to facilitate the simulation of LULC2017 using CA-Markov model. Then, the observed LULC2017 were cross-compared with the simulated LULC2017 to evaluate the performance of the model using the kappa index. Accordingly, the calculated Kappa index is calculated as 0.87. The result reveals a high level of agreement between the simulated and observed LULC2017, showing the CA-Markov model is strong to simulate the future LULC in the study area.
The LULC classes of Finchaa catchment were grouped into eight classes: agriculture, commercial farm, forestland, grazing land, rangeland, urban and built up, swampy and water bodies. From the land use maps of 2002 and 2017, Finchaa catchment was characterized by expansion of intensive agriculture and urbanization resulting in a natural vegetation decline. Commercial farm, urban and built-up and agricultural lands have increased as shown by Table 2. The highest decline in forest land from 2002 to 2017 was dominantly due to the expansion of commercial farms. The expansion in commercial lands was carried out, in order to double the sugar production by increasing the area of sugarcane cultivation. The future prediction of LULC 2036 and LULC 2055 were executed using the 2017 land use map as a base map. The transition matrix and transfer probability matrix were developed according to the LULC data of 2002 and 2017. The predicted results of 2017-2036, as shown in Table 2 and Figure 4, show agricultural land, commercial farm and urban and built-up areas are increasing, while forestland, grazing lands, rangelands, swampy and water bodies are decreasing. In 2036 to 2055, commercial farms and urban and built-up was predicted to increase continuously. However, the rate of agricultural land and commercial farm expansion is lower compared to the result of the previous prediction. This might be possible due to the limited area of land for the commercial farm expansion, as the farms are based on irrigation. Furthermore, the increase in range land through continuous future conservation works (soil and water conservation works are already started) might contribute to the decline of agricultural expansions. The high rate of urbanization could also contribute to this. The increase of the shrublands due to the planned afforestation could also be another possible reason.
The SWAT model requires the LULC classes to be reclassified in the code that the model recognizes. Accordingly, the LULC classes are defined with their equivalent representative, as shown in Figure 4.   Figure 4.

Sensitivity Analysis, Calibration and Validation
Sensitivity analysis for the simulated streamflow was performed using a daily observed flow to identify the most sensitive parameter with strong influence on model outputs. Initially, parameters related to surface runoff, groundwater, geomorphology, evaporation and soil water were considered and 11 parameters were identified as the most sensitive parameters for calibration. The given ranks are shown in Table 3. The calibration was carried outfitting the model to the streamflow record of 1990-2000. Afterwards, the calibrated values of these 11 parameters were used for further simulations. The parameter ranking was taken from the last iterations of SUFI-2 based on t-stat and p-stat. The larger in the absolute value of t-stat and the smaller the P-value, the more sensitive is the parameter [64]. As it can be seen from the table, the top three most sensitive parameters are CN2 (moisture condition II curve number), SOL_AWC (Available water capacity of the soil layer) and RCHRG_DP (Deep aquifer percolation fraction).

Sensitivity Analysis, Calibration and Validation
Sensitivity analysis for the simulated streamflow was performed using a daily observed flow to identify the most sensitive parameter with strong influence on model outputs. Initially, parameters related to surface runoff, groundwater, geomorphology, evaporation and soil water were considered and 11 parameters were identified as the most sensitive parameters for calibration. The given ranks are shown in Table 3. The calibration was carried outfitting the model to the streamflow record of 1990-2000. Afterwards, the calibrated values of these 11 parameters were used for further simulations. The parameter ranking was taken from the last iterations of SUFI-2 based on t-stat and p-stat. The larger in the absolute value of t-stat and the smaller the p-value, the more sensitive is the parameter [64]. As it can be seen from the table, the top three most sensitive parameters are CN2 (moisture condition II curve number), SOL_AWC (Available water capacity of the soil layer) and RCHRG_DP (Deep aquifer percolation fraction).
The calibration results on mean monthly flow show that SWAT model is able to capture the observed streamflow with R 2 , NSE, PBIAS and RSR of 0.71, 0.71, −1.1, and 0.54 respectively. The average monthly flow validation indicates R 2 , NSE, PBIAS and RSR of 0.81, 0.76, 14.3 and 0.49 respectively. This indicates that the performance of the SWAT model in the validation is good enough to simulate the stream flow in the Finchaa catchment. Furthermore, p-factor and R-factor statistics showed good agreement with 0.86 and 0.77 during calibration and 0.73, and 0.7 during validation, respectively. Overall, the statistics for goodness fit shows good agreement between the observed and simulated average monthly flow. However, the model is not able to simulate all peak flows well, as shown in Figure 5. respectively. This indicates that the performance of the SWAT model in the validation is good enough to simulate the stream flow in the Finchaa catchment. Furthermore, p-factor and R-factor statistics showed good agreement with 0.86 and 0.77 during calibration and 0.73, and 0.7 during validation, respectively. Overall, the statistics for goodness fit shows good agreement between the observed and simulated average monthly flow. However, the model is not able to simulate all peak flows well, as shown in Figure 5.

Climate Change Projections of the RCMs
To illustrate how each RCMs predicted precipitation and temperature, the projection of individual RCMs and ensemble mean of the RCMs were analyzed. The bias-corrected precipitation, and temperature within distribution mapping for future scenarios, are compared to the baseline data sets . Further, the hydrological responses of the catchment are analyzed by quantifying and comparing the future water balance components (surface runoff, groundwater flow, total water yield and evapotranspiration).

Rainfall
The individual and ensemble mean of the four RCMs were used to assess climate change in two future scenarios; near future (2021-2050) and mid future (2051-2080) with respect to the baseline period (1986-2015) under RCP4.5 and RCP8.5. The results of the projected precipitation by the

Climate Change Projections of the RCMs
To illustrate how each RCMs predicted precipitation and temperature, the projection of individual RCMs and ensemble mean of the RCMs were analyzed. The bias-corrected precipitation, and temperature within distribution mapping for future scenarios, are compared to the baseline data sets . Further, the hydrological responses of the catchment are analyzed by quantifying and comparing the future water balance components (surface runoff, groundwater flow, total water yield and evapotranspiration).

Rainfall
The individual and ensemble mean of the four RCMs were used to assess climate change in two future scenarios; near future (2021-2050) and mid future (2051-2080) with respect to the baseline period (1986-2015) under RCP4.5 and RCP8.5. The results of the projected precipitation by the individual RCMs show a different degree of precipitation changes. Except for HIRHAM5 in mid-term, all the RCMs show precipitation decline under both scenarios. According to CCLM4-8, precipitation decreases by −12.67% under RCP4.5 in near future and −27.47% under RCP8.5 scenarios in a midterm. The highest decline of the precipitation was simulated under RCP8.5 in a mid-term with CCLM4-8 and the lowest decline by HIRHAM5 under RCP4.5 in a midterm. HIRHAM5 shows that precipitation will increase under RCP8.5 scenarios. Summary of RCMs precipitation projections was presented in Table 4. The ensemble mean of the RCMs suggested decreasing precipitation under the two scenarios, like most of the RCMs. The change shows precipitation decline by −8.24% under RCP4.5 and −7.87% under RCP8.5 in the near-future scenario. Under the mid future, the decline was by −11.32% under RCP4.5 and −9.67% under RCP8.5. The use of ensemble mean over the individual RCMs helps to minimize the highest and lowest projections and enables to minimize the uncertainties by working with the average of the RCMs. Furthermore, the analysis of the inter-annular and inter-model variability of the climate variables in Finchaa catchment, presented in Dibaba et al. [54], underlines the need to include several climate models in order to cover the range of possible developments.
Changes in projected precipitation are more profound on seasonal bases compared to the annual bases in the catchment. The projected precipitation shows a declining trend in all seasons, except in April by RCP8.5, for both future time horizons and August by RCP4.5 in a near-future scenario. However, the changes are projected higher in the dry season than a wet season as shown by Figure 6. In particular, the seasons with smaller rainfall (March-April-May) are expected to face high precipitation changes while seasons with high rainfall (June-July-August-September) are expected to face lower precipitation change. For dry seasons, all projections give a higher decreasing signal. precipitation decreases by −12.67% under RCP4.5 in near future and −27.47% under RCP8.5 scenarios in a midterm. The highest decline of the precipitation was simulated under RCP8.5 in a mid-term with CCLM4-8 and the lowest decline by HIRHAM5 under RCP4.5 in a midterm. HIRHAM5 shows that precipitation will increase under RCP8.5 scenarios. Summary of RCMs precipitation projections was presented in Table 4. The ensemble mean of the RCMs suggested decreasing precipitation under the two scenarios, like most of the RCMs. The change shows precipitation decline by −8.24% under RCP4.5 and −7.87% under RCP8.5 in the near-future scenario. Under the mid future, the decline was by −11.32% under RCP4.5 and −9.67% under RCP8.5. The use of ensemble mean over the individual RCMs helps to minimize the highest and lowest projections and enables to minimize the uncertainties by working with the average of the RCMs. Furthermore, the analysis of the inter-annular and inter-model variability of the climate variables in Finchaa catchment, presented in Dibaba et al. [54], underlines the need to include several climate models in order to cover the range of possible developments.
Changes in projected precipitation are more profound on seasonal bases compared to the annual bases in the catchment. The projected precipitation shows a declining trend in all seasons, except in April by RCP8.5, for both future time horizons and August by RCP4.5 in a near-future scenario. However, the changes are projected higher in the dry season than a wet season as shown by Figure  6. In particular, the seasons with smaller rainfall (March-April-May) are expected to face high precipitation changes while seasons with high rainfall (June-July-August-September) are expected to face lower precipitation change. For dry seasons, all projections give a higher decreasing signal. In Ethiopia, a decreasing trend of precipitation intensity was reported in eastern, southwestern and southern regions using data ranges from 1965 to 2002 [68]. However, different studies in the Blue Nile basin reported varying changes of precipitation in the future horizon. Elshamy et al. [41] show In Ethiopia, a decreasing trend of precipitation intensity was reported in eastern, southwestern and southern regions using data ranges from 1965 to 2002 [68]. However, different studies in the Blue Nile basin reported varying changes of precipitation in the future horizon. Elshamy et al. [41] show 10 GCMs out of 17 GCMs reported a reduction of precipitation while the ensemble of the 17 GCMs reports almost no expected precipitation change. On the other hand, Beyene et al. [36] reported an increase in projected precipitation in the late 21st century using 11 GCMs. Large uncertainties among the RCMs-GCMs, regarding the signal of future precipitation, have also been reported in a tropical catchment in a Burkina Faso by de Hipt et al. [69]. The complexity in the nature of precipitation and its dependency on topographic and physical factors could contribute to this.
The results of precipitation projections in this study are consistent with the results presented by climate change studies [16,40]. The study by Daba and Rao [48] on climate change impacts on the hydro-meteorological variables of Finchaa sub-basin using HadCM GCM also reported decreasing precipitation over the study area. Similarly, the study by Girma [70], using CCLM downscaling for the periods 2041-2070 and 2071-2100, have shown declining precipitation by 6.6% and 6.4% over the upper Blue Nile basin.

Temperature
The variation of projected future maximum and minimum temperature is presented in Table 5 for the individual RCMs and their ensemble mean. The result shows that maximum and minimum temperatures increase under both RCPs throughout the study years considered showing the warming trends in the catchment. However, the magnitude of changes by individual RCMs and their mean ensemble is higher for higher emission scenarios (RCP8.5) compared to RCP4.5. Likewise, the change in temperature is higher in the mid future than the near future, under both emission scenarios. Furthermore, there are variations among all the RCMs and the ensemble mean for the degree of the temperature changes. RCA4 have shown the highest increase in maximum temperature under both RCPs during the nearby and mid future. Regarding the ensemble mean of the RCMs, the maximum temperature will increase on average by 1.34 • C and 2.15 • C under medium emission scenario (RCP4.5) and minimum temperature increases by 1.57 • C and 2.67 • C for the nearby and mid future, respectively. Likewise, the maximum temperature will increase by 1.49 • C and 3.21 • C and minimum temperature increases by 1.92 • C and 4.22 • C for the near future and mid future respectively under high emission scenarios RCP8.5. Overall, the study revealed the projection of RCP8.5 is warmer than RCP4.5. The highest temperature change in RCP8.5 is also confirmed by Shiferaw et al. [39] using five different GCMs in Ilala watershed, Northern Ethiopia.
The changes in projected temperature were not only expected to vary annually, but they also vary in all seasons. However, the seasonal variations are higher for the minimum temperature than the maximum temperature, as shown in Figure 7. The seasonal changes of maximum temperature is higher in dry seasons (December, January and February), while the rise in minimum temperature is higher in the wet season (June, July, and August). The highest seasonal rise of maximum temperature is reported by RCP8.5 during 2051-2080 while RCP 4.5 reported the lowest rise during 2021-2050.
Like the annual changes, the highest seasonal rise in temperature is expected by the high emission scenario than the medium emission scenario. Likewise, the highest rise in maximum and minimum temperature is expected in the midterm scenario (2051-2080) than the near future (2021-2050).
higher in the wet season (June, July, and August). The highest seasonal rise of maximum temperature is reported by RCP 8.5 during 2051-2080 while RCP 4.5 reported the lowest rise during 2021-2050.
Like the annual changes, the highest seasonal rise in temperature is expected by the high emission scenario than the medium emission scenario. Likewise, the highest rise in maximum and minimum temperature is expected in the midterm scenario (2051-2080) than the near future (2021-2050). In general, the projection of maximum temperature and minimum temperature in both future time horizons is within the range projected by IPCC and agrees with a range produced by other researchers. Many climate change studies in Ethiopia, in different regions, have shown the temperature is expected to increase. However, the strength of change varies with the techniques of downscaling and climate model types. The study reports by GFDRR [71]on climate risk and adaptation country profile, shows mean annual temperature is projected to increase by 1.1 °C to 3.1 °C by the 2060s and 1.5 °C to 5.1 °C by 2090s. The study by Girma [70] also showed the average annual temperature will increase by 1.5 °C, 2.6 °C and 4.5 °C for the periods 2011-2040, 2041-2070 and 2071-2100 respectively over the whole upper Blue Nile basin. The study on the whole Nile Basin by Beyene et al. [36] using the average of 11 GCMs further reported the annual average temperature change of 0.91 °C to 1.9°C during 2010-2039 for B1 global emission scenarios with respect to 1950-1999 historical periods. Further, the study by Elshamy et al. [41] using 17 GCMs on the Blue Nile projected temperature to increase between 2 °C and 5 °C.

Impacts of Land Use/Land Cover Change
The hydrological responses of the catchment to LULCC was investigated with two LULCC scenarios; first scenario for LULCC of 2017 to 2036 and second scenario for LULCC of 2017 to 2055 with the base land use map of 2017. The hydrological responses of the watershed were considered in terms of the process contributs to the annual surface runoff, annual groundwater, total water yield and potential evapotranspiration as shown in Table 6. Owing to the massive expected changes in land use from 2017 to 2036 and 2017 to 2055, mean annual surface water runoff and total water yield will increase whereas groundwater recharge and potential evapotranspiration will decrease. However, the rate of increase and decline was related to the rate of LULCCs. For example, the rate of agricultural expansion and forest decline was predicted to be higher in the first scenario and lower in the second. Consequently, the rate of increase in surface flow was lower, increasing from 2036 to 2055, compared to the baseline scenario. In general, the impacts of the LULCC are related to the increase of agricultural lands coupled with the decrease of forestlands. In general, the projection of maximum temperature and minimum temperature in both future time horizons is within the range projected by IPCC and agrees with a range produced by other researchers. Many climate change studies in Ethiopia, in different regions, have shown the temperature is expected to increase. However, the strength of change varies with the techniques of downscaling and climate model types. The study reports by GFDRR [71] on climate risk and adaptation country profile, shows mean annual temperature is projected to increase by 1.1 • C to 3.1 • C by the 2060s and 1.5 • C to 5.1 • C by 2090s. The study by Girma [70] also showed the average annual temperature will increase by 1.5 • C, 2.6 • C and 4.5 • C for the periods 2011-2040, 2041-2070 and 2071-2100 respectively over the whole upper Blue Nile basin. The study on the whole Nile Basin by Beyene et al. [36] using the average of 11 GCMs further reported the annual average temperature change of 0.91 • C to 1.9 • C during 2010-2039 for B1 global emission scenarios with respect to 1950-1999 historical periods. Further, the study by Elshamy et al. [41] using 17 GCMs on the Blue Nile projected temperature to increase between 2 • C and 5 • C.

Impacts of Land Use/Land Cover Change
The hydrological responses of the catchment to LULCC was investigated with two LULCC scenarios; first scenario for LULCC of 2017 to 2036 and second scenario for LULCC of 2017 to 2055 with the base land use map of 2017. The hydrological responses of the watershed were considered in terms of the process contributs to the annual surface runoff, annual groundwater, total water yield and potential evapotranspiration as shown in Table 6. Owing to the massive expected changes in land use from 2017 to 2036 and 2017 to 2055, mean annual surface water runoff and total water yield will increase whereas groundwater recharge and potential evapotranspiration will decrease. However, the rate of increase and decline was related to the rate of LULCCs. For example, the rate of agricultural expansion and forest decline was predicted to be higher in the first scenario and lower in the second. Consequently, the rate of increase in surface flow was lower, increasing from 2036 to 2055, compared to the baseline scenario. In general, the impacts of the LULCC are related to the increase of agricultural lands coupled with the decrease of forestlands. The increase in surface runoff and water yield is 5.97% and 1.9% in the first scenario and 7.1% and 2.26% in the second scenario, respectively. The groundwater and potential evapotranspiration declines by 1.09% and 1.15% in the first and 0.91%, and 1.36% in the second scenario, respectively. As it can be seen from Figure 8, the decline of groundwater in the second scenario was lower than the first scenario owing to the increase of the predicted shrublands in 2055, compared to 2036 LULC. Hence, the decreasing rate is higher for the changes in 2017 to 2036, where the rate of forest and shrubland decline is higher. This is a good indication to see the effect of works with changing the land cover in groundwater management i.e., if shrubland increases, then groundwater recharge will increase. The decrease in evapotranspiration was due to the decline of evapotranspiration owing to the decline of forestland. Further, the decline in evapotranspiration could also be related to the shortage of soil moisture as a result of increased runoff. A similar study showed that PET decreases with the shortage of soil moisture as a result of reduced rainfall [72].  The increase in surface runoff and water yield is 5.97% and 1.9% in the first scenario and 7.1% and 2.26% in the second scenario, respectively. The groundwater and potential evapotranspiration declines by 1.09% and 1.15% in the first and 0.91%, and 1.36% in the second scenario, respectively. As it can be seen from Figure 8, the decline of groundwater in the second scenario was lower than the first scenario owing to the increase of the predicted shrublands in 2055, compared to 2036 LULC. Hence, the decreasing rate is higher for the changes in 2017 to 2036, where the rate of forest and shrubland decline is higher. This is a good indication to see the effect of works with changing the land cover in groundwater management i.e., if shrubland increases, then groundwater recharge will increase. The decrease in evapotranspiration was due to the decline of evapotranspiration owing to the decline of forestland. Further, the decline in evapotranspiration could also be related to the shortage of soil moisture as a result of increased runoff. A similar study showed that PET decreases with the shortage of soil moisture as a result of reduced rainfall [72]. The spatial distribution of the hydrological process under the projected LULCC for 2017-2036 is presented by Figure 9a,b,c,d and for 2017-2055 presented by Figure 9e,f,g,h. Each sub-basin has different characteristics of a hydrological process. The highest increase in surface runoff is around the downstream while the upstream area shows a decline of the surface runoff under both scenarios with the highest increase and decrease expected in the second scenario. Areas with a high decline of surface runoff show an increase in groundwater. This demonstrates that surface runoff is inversely correlated with groundwater under both scenarios. The high decline in PET is expected around the central parts of the catchment, while the upstream areas are characterized by low decline. Water yield shows high decline around the northeast part in the first scenario and southwestern part in the second scenario. However, the majority of the sub-basins show an increase in water yield with the highest increase shown by the second scenario. The highest increase in surface flow and decline of groundwater recharge is associated with the highest expansion of agricultural lands and urban expansion and decline of forest and shrubland. This reveals a massive conversion of LULC to intensive agriculture and settlements will reduce the soil infiltrations capacity, causing a large portion of rainfall directly changed to surface runoff. The reduction of soil water infiltration, in turn, causes groundwater flow decline. Furthermore, the increasing demand for the growing population to use groundwater could additionally contribute to the decline of groundwater by increasing groundwater withdrawal. The recent evidence by Dibaba et al. [2] in Finchaa catchment reported that, due to the high changes in LULC, springs in the catchment are dried out, and the level of water in hand-bug wells decline. Further, the finding agrees with the study by Shi et al. [73] who indicated the main effects of LULCC on the water cycle are changes in evapotranspiration (ET), the soil's ability to hold water and in the abilities of vegetation to intercept precipitation. Water 2020, 12, x FOR PEER REVIEW 16 of 26 The highest increase in surface flow and decline of groundwater recharge is associated with the highest expansion of agricultural lands and urban expansion and decline of forest and shrubland. This reveals a massive conversion of LULC to intensive agriculture and settlements will reduce the soil infiltrations capacity, causing a large portion of rainfall directly changed to surface runoff. The reduction of soil water infiltration, in turn, causes groundwater flow decline. Furthermore, the increasing demand for the growing population to use groundwater could additionally contribute to the decline of groundwater by increasing groundwater withdrawal. The recent evidence by Dibaba et al. [2] in Finchaa catchment reported that, due to the high changes in LULC, springs in the catchment are dried out, and the level of water in hand-bug wells decline. Further, the finding agrees with the study by Shi et al. [73] who indicated the main effects of LULCC on the water cycle are changes in evapotranspiration (ET), the soil's ability to hold water and in the abilities of vegetation to intercept precipitation.
The increase in surface water has increased the total water yield of the catchment. This indicates that degraded watersheds can accelerate surface runoff by reducing both the retention as soil water and groundwater recharge. Consequently, the catchment will face severe soil erosion, landslides and flash floods. This could become a great concern where mega hydropower projects exist. Overall, surface runoff, water yield and potential evapotranspiration were found to be more sensitive than the groundwater flow in Finchaa catchment under LULCC.
The findings of the study are consistent with various research studies conducted in other catchments. For instance, Gashaw et al. [74] in Andessa catchment, Ethiopia concluded that the changes in land use have reduced groundwater flow and increased surface runoff. Shawul et al. [18] also showed an increased surface runoff and decline of groundwater between 1974-2014 in upper Awash basin was attributed by the LULCCs of the catchment. However, the study reported that the magnitude of the groundwater fluctuation is much less compared to the surface runoff. A previous study in the Finchaa catchment, in LULCC from 1987 to 2017, shows that construction of Fincha, The increase in surface water has increased the total water yield of the catchment. This indicates that degraded watersheds can accelerate surface runoff by reducing both the retention as soil water and groundwater recharge. Consequently, the catchment will face severe soil erosion, landslides and flash floods. This could become a great concern where mega hydropower projects exist. Overall, surface runoff, water yield and potential evapotranspiration were found to be more sensitive than the groundwater flow in Finchaa catchment under LULCC.
The findings of the study are consistent with various research studies conducted in other catchments. For instance, Gashaw et al. [74] in Andessa catchment, Ethiopia concluded that the changes in land use have reduced groundwater flow and increased surface runoff. Shawul et al. [18] also showed an increased surface runoff and decline of groundwater between 1974-2014 in upper Awash basin was attributed by the LULCCs of the catchment. However, the study reported that the magnitude of the groundwater fluctuation is much less compared to the surface runoff. A previous study in the Finchaa catchment, in LULCC from 1987 to 2017, shows that construction of Fincha, Neshe and Amerti hydropower projects along with irrigation projects have brought a dramatic LULCC [2]. Consequently, the scarcity of land coupled with the need to farm more lands has led the community to cultivate steep slopes. The cultivation of steep slopes, in turn, increased the vulnerability of land to soil erosion, which reduced the water resources of the catchment.

Impacts of Climate Change
Changes in rainfall and temperature were used to predict the impacts of climate change on water balance components of the catchment. Consequently, the projected changes in mean annual precipitation and temperature, under the two RCPs, caused a significant variation in the projected water balance components of Finchaa catchment.
SWAT simulations for the near future and mid future showed that the decline of precipitation and increase of temperature will lead to reduced surface flow, groundwater and overall water yield (Table 7). Changes in temperature and PET are correlated positively. Consequently, the increase in temperature resulted in increased potential evapotranspiration and evaporation, which could be a critical factor for the decline of water yield. As the forms of water are subjected to losses, owing to the changes in temperature, evaporation is also a factor for the future decline of the surface runoff. This is expected considering the warming trends of temperature. In general, the impacts of the climate changes on the watershed hydrology can be explained with the combinations of precipitation, temperature and evapotranspiration. Decreased rainfall reduces surface runoff and increasing temperatures amplified the increase in evapotranspiration, which may have resulted in the decreased runoff. In this regard, the correlation between precipitation and the simulated surface runoff, groundwater and total water yield is positive, while the relationship between the temperature and hydrological process are positive only for PET. However, the relationship between the precipitation and temperature with the hydrological process varies with the RCPs and period of the scenario. For example, changes in precipitation and the predicted decline of groundwater and total water yields are higher during 2051-2080 under both RCPs. However, the decline of surface runoff is only higher under RCP4.5 during the 2051-2080 (14.48%). Under RCP8.5, surface water decline is higher during 2021-2050 (−12.32%) as shown in Figure 10. This shows that change in precipitation alone is not the only variable that affects surface water. This finding was consistent with other research reports in the Nile River [41,75]. According to Coffel et al. [75], dry and hot years driven by increasing regional temperatures are becoming more frequent in the upper Nile.
vulnerability of land to soil erosion, which reduced the water resources of the catchment.

Impacts of Climate Change
Changes in rainfall and temperature were used to predict the impacts of climate change on water balance components of the catchment. Consequently, the projected changes in mean annual precipitation and temperature, under the two RCPs, caused a significant variation in the projected water balance components of Finchaa catchment.
SWAT simulations for the near future and mid future showed that the decline of precipitation and increase of temperature will lead to reduced surface flow, groundwater and overall water yield (Table 7). Changes in temperature and PET are correlated positively. Consequently, the increase in temperature resulted in increased potential evapotranspiration and evaporation, which could be a critical factor for the decline of water yield. As the forms of water are subjected to losses, owing to the changes in temperature, evaporation is also a factor for the future decline of the surface runoff. This is expected considering the warming trends of temperature. In general, the impacts of the climate changes on the watershed hydrology can be explained with the combinations of precipitation, temperature and evapotranspiration. Decreased rainfall reduces surface runoff and increasing temperatures amplified the increase in evapotranspiration, which may have resulted in the decreased runoff. In this regard, the correlation between precipitation and the simulated surface runoff, groundwater and total water yield is positive, while the relationship between the temperature and hydrological process are positive only for PET. However, the relationship between the precipitation and temperature with the hydrological process varies with the RCPs and period of the scenario. For example, changes in precipitation and the predicted decline of groundwater and total water yields are higher during 2051-2080 under both RCPs. However, the decline of surface runoff is only higher under RCP4.5 during the 2051-2080 (14.48%). Under RCP8.5, surface water decline is higher during 2021-2050 (−12.32%) as shown in Figure 10. This shows that change in precipitation alone is not the only variable that affects surface water. This finding was consistent with other research reports in the Nile River [41,75]. According to Coffel et al. [75], dry and hot years driven by increasing regional temperatures are becoming more frequent in the upper Nile. The spatial distribution of the hydrological process under the projected climate change for RCP4.5 and RCP8.5 is presented in Figure 11a-d for surface runoff, Figure 11e-h for groundwater, Figure 11i-l for water yield and Figure 11m-p for PET. Although the impacts of climate change reported a reduction of surface runoff, groundwater and water yield in the catchment, each sub-basin is characterized by different yields. The changing spatiotemporal patterns and magnitudes of precipitation are expected to directly disturb the hydrological process, while indirect impacts are expected from changes in temperature. Commonly, the highest decline of the surface runoff, groundwater, and water yield is projected around the upstream of the catchment, under both scenarios throughout the study projections considered. The projection of surface runoff and water yield decline is higher under RCP4.5 and RCP8.5 scenarios in the mid future than the near future. The lowest decline of surface runoff under both scenario is projected around the southwestern part in the near future and northeastern in the mid future. Groundwater and water yield decline is higher around the southwestern part. The highest rise in PET is shown by RCP8.5 around the outlet of the catchment in the mid future projection.
The spatial distribution of the hydrological process under the projected climate change for RCP4.5 and RCP 8.5 is presented in Figure 11a-11d for surface runoff, Figure 11e-11h for groundwater, Figure 11i-11l for water yield and Figure 11m-11p for PET. Although the impacts of climate change reported a reduction of surface runoff, groundwater and water yield in the catchment, each sub-basin is characterized by different yields. The changing spatiotemporal patterns and magnitudes of precipitation are expected to directly disturb the hydrological process, while indirect impacts are expected from changes in temperature. Commonly, the highest decline of the surface runoff, groundwater, and water yield is projected around the upstream of the catchment, under both scenarios throughout the study projections considered. The projection of surface runoff and water yield decline is higher under RCP4.5 and RCP8.5 scenarios in the mid future than the near future. The lowest decline of surface runoff under both scenario is projected around the southwestern part in the near future and northeastern in the mid future. Groundwater and water yield decline is higher around the southwestern part. The highest rise in PET is shown by RCP8.5 around the outlet of the catchment in the mid future projection. Generally, the decrease in average annual flow over the catchment could be due to the annual and seasonal decrease of precipitation and increase in temperature under both scenarios. The possibility of surface water, groundwater and total water yield reduction in the catchment could affect the availability of water resources in the catchment and further aggravate water stress in the downstream. The acceleration of the increasing water abstractions could also contribute to the further decline in water resources besides the climate change pressure.
The increase in PET due to high-temperature increase coupled with decreasing precipitation could also lead to reduced soil moisture, which is required for plant growth and groundwater storage. This could imply the reduced availability of water for crop production, which will become chronic to the farmers in the catchment whose livelihood is based on agriculture alone. Further, an increase in temperature led to increased PET, this in turn, resulted in an increased need for irrigation.
Climate change studies on Africa, sub-Sahara, in particular, confirmed similar findings as climate projections show a range of warming trends in inland subtropics, a frequent occurrence of Figure 11. Spatial patterns of the changes in surface runoff (a-d), groundwater (e-h), water yield (i-l) and PET (m-p) due to CC.
Generally, the decrease in average annual flow over the catchment could be due to the annual and seasonal decrease of precipitation and increase in temperature under both scenarios. The possibility of surface water, groundwater and total water yield reduction in the catchment could affect the availability of water resources in the catchment and further aggravate water stress in the downstream. The acceleration of the increasing water abstractions could also contribute to the further decline in water resources besides the climate change pressure.
The increase in PET due to high-temperature increase coupled with decreasing precipitation could also lead to reduced soil moisture, which is required for plant growth and groundwater storage. This could imply the reduced availability of water for crop production, which will become chronic to the farmers in the catchment whose livelihood is based on agriculture alone. Further, an increase in temperature led to increased PET, this in turn, resulted in an increased need for irrigation.
Climate change studies on Africa, sub-Sahara, in particular, confirmed similar findings as climate projections show a range of warming trends in inland subtropics, a frequent occurrence of extreme heat events, increasing aridity and changes in rainfall [76]. According to Serdeczny et al. [76] the impacts due to climate change add significant impact on the existing undernutrition, infectious disease, the vulnerability of rain-fed agriculture and flash flooding crises in the region.
The finding also confirms the findings of Beyene et al. [36] which showed the Nile river is expected to decline during the mid. 2040-2069 and late 2070-2099 due to the decline in precipitation and increased evaporation demand. The study by Gebremeskel and Kebede [77] on climate change impacts on water sources of the Tekeze river basin also predicted the decline of surface runoff by 13% and 14%.
The study by Shiferaw et al. [39] using the ensemble of five GCMs also confirmed the decline of the projected surface runoff under both RCPs due to climate change.

Impacts of Land Use/Land Cover and Climate Change
Although the direct effect of a decline in precipitation and increased temperature disturbs the hydrology of the catchment as presented under the impacts of climate change, changes in LULC with soil and vegetation properties, either enhance or decrease the impacts of climate change. In this regard, the effects of combined LULC and climate change on the hydrological process was done to distinguish the effects of LULC and climate change under the isolated and combined scenarios.
The increase of temperature and decreasing precipitation under a varying LULC in Finchaa catchment shows a decline of surface runoff, groundwater and total water yield, whereas the projected evapotranspiration is increasing. The combined effects of LULC and climate change on the watershed hydrological process are presented in Table 8. The highest annual surface runoff decline of 9.09% is projected under the RCP4.5 in a midterm scenario, and the highest groundwater decline of 15.30% is projected under RCP8.5 in the midterm period. As presented in Figure 12, both RCPs projected a higher decline of the water yield and groundwater during the mid-future scenario. Similarly, the highest increase in PET is projected during the mid-future for both RCPs.
the impacts due to climate change add significant impact on the existing undernutrition, infectious disease, the vulnerability of rain-fed agriculture and flash flooding crises in the region.
The finding also confirms the findings of Beyene et al. [36] which showed the Nile river is expected to decline during the mid. 2040-2069 and late 2070-2099 due to the decline in precipitation and increased evaporation demand. The study by Gebremeskel and Kebede [77] on climate change impacts on water sources of the Tekeze river basin also predicted the decline of surface runoff by 13% and 14%. The study by Shiferaw et al. [39] using the ensemble of five GCMs also confirmed the decline of the projected surface runoff under both RCPs due to climate change.

Impacts of Land Use/Land Cover and Climate Change
Although the direct effect of a decline in precipitation and increased temperature disturbs the hydrology of the catchment as presented under the impacts of climate change, changes in LULC with soil and vegetation properties, either enhance or decrease the impacts of climate change. In this regard, the effects of combined LULC and climate change on the hydrological process was done to distinguish the effects of LULC and climate change under the isolated and combined scenarios.
The increase of temperature and decreasing precipitation under a varying LULC in Finchaa catchment shows a decline of surface runoff, groundwater and total water yield, whereas the projected evapotranspiration is increasing. The combined effects of LULC and climate change on the watershed hydrological process are presented in Table 8. The highest annual surface runoff decline of 9.09% is projected under the RCP4.5 in a midterm scenario, and the highest groundwater decline of 15.30% is projected under RCP8.5 in the midterm period. As presented in Figure 12, both RCPs projected a higher decline of the water yield and groundwater during the mid-future scenario. Similarly, the highest increase in PET is projected during the mid-future for both RCPs.  The spatial distribution of the hydrological process under the projected LULC and climate change under the RCP4.5 and RCP8.5 are presented in Figure 13a-d for surface runoff, Figure 13e-h for groundwater, Figure 13i-l for water yield and Figure 13m-p for PET. The reduction in surface runoff is lower around the downstream of the catchment under all scenarios, while the high decline is around upstream of the catchment. Higher groundwater decline is expected around the upstream and outlet of the catchment, with the highest decline under RCP8.5 in the late scenario. Water yield is projected to decline higher around the upstream while a high increase in PET is projected around the downstream of the catchment. the upstream and outlet of the catchment, with the highest decline under RCP8.5 in the late scenario. Water yield is projected to decline higher around the upstream while a high increase in PET is projected around the downstream of the catchment. The increasing temperature is the most dominant factor affecting surface runoff in climate change scenarios while increasing urbanization is the main contributor affecting the surface runoff in LULC scenarios. Overall, the increase in surface runoff and total water yield due to LULCC was offset by the decline of the surface runoff and total water yield, due to climate change impacts. Consequently, the effect of climate change is decisive compared to the LULCC in the Finchaa catchment. The high sensitivity of water balance components and streamflow to the climate change is also reported by other studies in different parts of the world [22,72,73,78,79]. For example, the study by Shi et al. [73] in the upstream of Huai River, China reported that the combined effect of LULC and climate change has increased surface water, evapotranspiration and streamflow. However, the differentiated impacts of LULC and climatic variability on hydrological process shows a change in streamflow because of the effects of climate variability on the hydrological process were offset by the effect of LULCC. A combined study of LULC and climate change in Hoeya River Basin of Korea by Kim et al. [79] also shows that the combined scenario was similar to that of climate change only scenario.
In summary, the seasonal and annual variation of the future temperature shows the increased hot and dry years that will lead to serious water scarcity. Consequently, the increase in warming and evaporation in the catchment highly affects the three reservoirs in the catchment increasing water stress. The recent study by Coffel et al. [75] also revealed that the increase in regional temperature indicates the water balance of the Blue Nile basin may become moisture constrained in the future. The study report also shows the trend continues despite the increasing precipitation projected by climate models. The decline of the runoff amplifies the effects of water security in the catchment. With the LULC and climate changes, the soil moisture required for crop growth is reduced, and land degradation problems of the Finchaa catchment, occur as the landslide and gully become severe. The research studies by Vlek et al. [13] further reported that land degradations are related to the deterioration in climate conditions and human intervention. The increasing temperature is the most dominant factor affecting surface runoff in climate change scenarios while increasing urbanization is the main contributor affecting the surface runoff in LULC scenarios. Overall, the increase in surface runoff and total water yield due to LULCC was offset by the decline of the surface runoff and total water yield, due to climate change impacts. Consequently, the effect of climate change is decisive compared to the LULCC in the Finchaa catchment. The high sensitivity of water balance components and streamflow to the climate change is also reported by other studies in different parts of the world [22,72,73,78,79]. For example, the study by Shi et al. [73] in the upstream of Huai River, China reported that the combined effect of LULC and climate change has increased surface water, evapotranspiration and streamflow. However, the differentiated impacts of LULC and climatic variability on hydrological process shows a change in streamflow because of the effects of climate variability on the hydrological process were offset by the effect of LULCC. A combined study of LULC and climate change in Hoeya River Basin of Korea by Kim et al. [79] also shows that the combined scenario was similar to that of climate change only scenario.
In summary, the seasonal and annual variation of the future temperature shows the increased hot and dry years that will lead to serious water scarcity. Consequently, the increase in warming and evaporation in the catchment highly affects the three reservoirs in the catchment increasing water stress. The recent study by Coffel et al. [75] also revealed that the increase in regional temperature indicates the water balance of the Blue Nile basin may become moisture constrained in the future. The study report also shows the trend continues despite the increasing precipitation projected by climate models. The decline of the runoff amplifies the effects of water security in the catchment. With the LULC and climate changes, the soil moisture required for crop growth is reduced, and land degradation problems of the Finchaa catchment, occur as the landslide and gully become severe. The research studies by Vlek et al. [13] further reported that land degradations are related to the deterioration in climate conditions and human intervention.
Although the study did not cover the quantification of sediment yield, the SWAT simulation, under the LULCC, shows that the sediment yield increases with the changes. However, sediment yield is decreasing with the projected scenarios of climate change. The study by Gadissa et al. [40] on the effects of climate change on sedimentation in the Rift valley basin of Ethiopia reported a decrease of sediment yield by 38% under RCP4.5 and by 23% under RCP8.5. Therefore, the study further suggests the need to study the impacts of the combined LULC and climate change on the sedimentation in Finchaa catchment for better sediment management.

Conclusions
Four RCMs and their ensemble mean, under CORDEX-Africa, were applied for the climate change study. The simulations of all RCMs, except HIRHAM5, show decreasing precipitation under both high and medium-low emission scenarios. With respect to temperature, all RCMs projections show an increasing temperature with varying degree of changes. To avoid the considerable variations of the individual RCMs in projecting precipitation and temperature, the ensemble mean of the RCMs were used for the hydrological impact studies of climate changes, and the combined LULC and climate change. Although uncertainties of climate prediction still exist, the use of the ensemble RCMs was found to be a suited strategy to evaluate the uncertainties of individual RCMs.
The findings from this work were applied to distinguish the effects of LULC and climate change on the water balance components of the Finchaa catchment. Consequently, the metro-hydrological process could change in probability and intensity due to climate change, but also due to local land use/land cover change. If climate change is considered alone, the increase in temperature and decrease of predicted precipitation will decrease the surface runoff, groundwater and total water yield, whereas potential evapotranspiration and evaporation increase. However, the surface runoff increases under land use/land cover changes, due to the highly expanding urbanization and intensive agriculture. Overall, the sensitivity of water resources for land use/land cover and climate change over Finchaa catchment shows the effect of land use/land cover change is stunned by the effects of climate change. Consequently, climate change is found to be predominant over the effects of LULCC.
Further, the seasonal and annual variation of the future precipitation and temperature in the catchment shows the increasingly hot and dry years that will lead to serious water scarcity. Consequently, the catchment is highly vulnerable to climate change and its impacts could range from warming to crop failures as a result of prolonged dry seasons. The Upper Blue Nile basin will face regional water scarcity regardless of whether the precipitation is increased, constant or decreased [75]. The decline of the runoff amplifies the effects of water security in the catchment.
The study analysis suggests that the impact assessments of the combined land use/land cover and climate change in the Finchaa catchments are well-represented by the SWAT simulations. However, the limited availability and quality of hydro-climatic data in the region need urgent attention to improve our understanding of the change in existing and future climate and LULCCs.
The findings of this study provided important information on the relative influences of how the watershed hydrological process in Finchaa catchment respond to the changes in land use/land cover and climate change. This could help to plan proper water resources management interventions. If the degraded sloppy lands are rehabilitated, the ground recharge increase and the surface runoff, which washes the topsoil into the lakes, is reduced. Furthermore, the three lakes in the catchment should be buffered with proper management strategies. In general, the result highlights the need for regional developments and cooperation to urge strong climate-resilient management strategies and to counteract the rapid climate changes in the catchment.