Assessment of Surface Irrigation Potential of the Dhidhessa River Basin, Ethiopia

Assessing available water resources and their potential for irrigation water use is vital for sustainable agricultural development and planning. This is particularly of interest in developing countries like Ethiopia, where a small portion of largely accessible land for surface irrigation applications has been utilized, despite the majority of the population relying on agricultural productivity. This study utilized the Dhidhessa River Basin (Ethiopia) as a case study and analyzed the main challenges to balance the sustainable water resources utilization and enhance agricultural productivity of the basin. The study mainly focused on estimating the available water resources and their potential for surface irrigation water use in the basin. This was achieved by utilizing Geographic Information System (GIS)-based tools, a hydrological Soil and Water Assessment Tool (SWAT) model, and a Crop Water and Irrigation Requirements Program of FAO (CROPWAT) model. While the SWAT estimated the water availability in the basin, GIS-tools such as Model Builder were used to map the irrigation potential of the basin. For irrigation water potential assessment, we selected six crops (cabbage, maize, tomato, pepper, groundnut and sugarcane) and estimated their irrigation water requirements using the CROPWAT model. We developed the SWAT model for the period from 1986 to 2012 using the available hydro-meteorological and geo-spatial data. Due to many parameters used in the model, we first performed a parameter sensitivity analysis and identified the most essential/sensitivity parameters via Sequential Uncertainty Fitting-II (SUFI-2). The identified sensitive parameters were subsequently used for model calibration (1989–2000) and validation (2001–2012) procedures achieved via SUFI-2. SWAT was able to reproduce the observed monthly streamflow values with a coefficient of determination (R2) and Nash-Sutcliffe Coefficient (NSE) of 0.85 and 0.87 for the calibration period and 0.91 and 0.89 for the validation period, respectively. The findings generally indicated a “good” performance of the model in simulating the hydrology. The annual available water of the basin is 9.26 billion cubic meters (BCM) whereas the 70% and 80% dependable flow is 7.56 and 6.97 BCM, respectively. Based on the Model Builder of ArcGIS, the SWAT estimated available water can potentially irrigate an area of 259,028 ha for slope less than 8%, 643,162 ha for slopes less than 15% and 1,023,581 ha for slopes less than 30%. Moreover, the irrigation water requirements were calculated by the CROPWAT model for the six selected crops indicated that although the need for irrigation water varies depending on the season, the potential irrigation area of the Dhidhessa River Basin is greater than its irrigated land. Therefore, it is concluded that the basin’s surface irrigation Hydrology 2020, 7, 68; doi:10.3390/hydrology7030068 www.mdpi.com/journal/hydrology Hydrology 2020, 7, 68 2 of 21 systems need to be expanded to enhance the agricultural productivity and improve the livelihood of the basin’s communities and similar basins elsewhere.


Introduction
Due to population growth and increases in firewater resource demand, it has been become a challenge to sustain the unevenly distributed water resources worldwide [1]. In addition, the geo-morphological conditions and climate change/variability have further exacerbated freshwater availability and thus increases its demand, both temporally and spatially [1][2][3]. Globally, climate change has been putting pressure on the problem, modifying hydrological responses and shifts in time of rainfall-runoff responses [1,4,5]. Consequently, it has been documented that about two-thirds of the world's population have been, or will be, suffering from water shortage or unavailability [1,6,7]. These negative implications have been more pronounced in developing countries that mainly depend on rain-fed agricultural systems. Thus, such countries, including Ethiopia, may need to migrate from the traditional agricultural practices towards irrigation-based systems [8]. However, the amount of water available for irrigation depends on the availability of freshwater resources [9][10][11][12][13]. At the same time, agricultural areas have also been further threatened by the expansion of the residential and development areas, commercial/industrials, and other sectors [14,15].
Some studies have already shown the negative and potential impact of climate change on water resources supply for agriculture, urbanization, and other environment uses [1,[16][17][18]. While both groundwater and surface water are very susceptible to climate change [16,17,[19][20][21], some researchers have reported that there is a limited understanding of the dynamic relationship between groundwater and climate change [19,22]. In addition, the impact of climate change could be intensified by human interference such as irrigation water use and land use change [19,22]. Therefore, the world has been facing the problem of tackling the impact of climate and land use changes, population growth and urbanization on water resources and agricultural production and products [16,23].
The available water from surface water, renewable groundwater, and rainwater can be used for agriculture and irrigation [24][25][26]. However, such water sources have been severely under-developed and under-utilized in Ethiopia, implying the importance of enhancing the irrigation water potential of the available water resources [27,28]. Ethiopia mainly depends on rain-fed agriculture with limited use of irrigation for agricultural production [29]. Approximately >85% of the country's food supply has been generated from rain fed agriculture and low-productive subsistence crops [30]. Hence, rainfall is the most determining factor for food production, supply and the economy of the country [30]. However, increases in temperature and the alteration of rainfall patterns in Ethiopia have direct effects on agricultural crop growth and yields [31][32][33]. The major problem associated with rainfall-dependent agriculture in the country is the high degree of rainfall variability and unreliability [34,35]. Such conditions might cause crop failures due to the frequent occurrence of dry spells and droughts. This in turn often results in famine, particularly affecting the livelihoods of the rural and poor communities. Consequently, food and nutrition insecurity are a frequent and bottleneck problem in Ethiopia. Therefore, the country needs alternative approaches, for example enhancing irrigation water use from its abundant rivers, such as the Dhidhessa River Basin (DRB), and thus improving the agricultural productivity of the rural and poor communities [36].
The Dhidhessa River is one of the biggest tributaries of the upper Blue Nile River (the Abbay River in Ethiopia) [37]. The Dhidhessa River accounts for approximately about one-fourth of the Blue Nile River (BNR) [38][39][40][41][42]. Given its major contributions to the BNR, the Dhidhessa River can be considered as one of the potential water resources to alleviate agricultural productivity and food security problems through the use of advanced irrigation systems [38]. The land and water resource potential assessment using hydrological models has been important to ensure the sustainable development [43][44][45]. On the other hand, the evaluation of basin irrigation potential and irrigation water requirement is crucial in order to enhance food security [2,31,[46][47][48][49]. Previous studies have proven the applicability of the Soil and Water Assessment Tool (SWAT) model for irrigation potential assessment of the BNR Basins even with limited data [50][51][52]. SWAT is an open access and widely used model in poorly gauged basins like the BNR basin [50,52,53]. For example, the irrigation potential of Lake Tana, which is one of the contributors to the BNR, was extensively studied compared to the other sub basins of the BNR [54][55][56][57][58]. Although the Dhidhessa River Basin has comparatively sufficient hydrological and meteorological time series data (compared to the other sub-basins of the BNR), its hydrological and land resource development have not been well investigated yet. More specifically, the basin has not been well studied and understood as compared with the northern sub-basins of the Blue Nile River such as Tana sub-basin. There is limited information regarding the water resources development and surface irrigation potential of the Dhidhessa River Basin. There is a need to secure food security and water resource development in the country. The potential irrigation assessment is the first step towards ensuring sound planning and sustainability of irrigation development and ultimately to secure food and nutrition security [59]. Hence, this study was performed using the Soil and Water Assessment Tool (SWAT) and GIS as a tool for assessing the irrigation potential in the Dhidhessa River sub-basin. Furthermore, the study attempted to estimate the water resource potential of the river catchments in the basin and identify irrigable areas. The objective of this study is to demonstrate the available resources of the basin for sustainable development. As such, irrigation potential assessment is the first step towards ensuring planning and suitability of irrigation development in the area.

Description of the Study Area
The Dhidhessa River Basin (DRB) is located in the western part of Ethiopia between approximately 7 • -9 • N and 35 • -36 • E. The overall elevation in the study area varies between 954 and 3154 m above the mean sea level (Figure 1). The river drains about 14,710 km 2 . Most of the DRB is classified under a humid tropical climate with heavy rainfall with an annual amount that is mainly received during the rainy season [36]. The highest and lowest temperature ranges between 21-31 • C and 10-15 • C, respectively. The mean annual flow of the DRB at the Arjo station ( Figure 2) is about 341 and 318 m 3 /s, respectively. The river reaches its maximum flow occurring during August to September.  The basin is characterized with diversified land use/land cover that includes moderately and intensively cultivated land (50%), dense woodland (7%), wooded grassland (8%), open woodland (7%), natural forest cover (9%), natural forest with coffee (5%), coffee farm with shade trees (6%), riverine forest, bamboo forest, plantation forest, settlement (1%), shrublands and open grassland (7%). The cultivated land consists of mixed cultivation, coffee production, livestock production, subsistence and commercial forest products utilization, non-timber products utilization, beekeeping, wildlife management and utilization, infrastructure development, mining and different investment activities [60,61].

Data Collection and Processing
Daily river discharge and meteorological data (rainfall, temperature, relative humidity, sunshine hours, and wind speed) were used in this study [62,63]. Digital elevation model (DEM) of 30 × 30 m resolution was used as input for the SWAT model. Based on the DEM data, sub-basins were delineated and other essential physiographic parameters were estimated for water resources estimation. The physiographic parameters include, among others, the hydrological response units The basin is characterized with diversified land use/land cover that includes moderately and intensively cultivated land (50%), dense woodland (7%), wooded grassland (8%), open woodland (7%), natural forest cover (9%), natural forest with coffee (5%), coffee farm with shade trees (6%), riverine forest, bamboo forest, plantation forest, settlement (1%), shrublands and open grassland (7%). The cultivated land consists of mixed cultivation, coffee production, livestock production, subsistence and commercial forest products utilization, non-timber products utilization, beekeeping, wildlife management and utilization, infrastructure development, mining and different investment activities [60,61].

Data Collection and Processing
Daily river discharge and meteorological data (rainfall, temperature, relative humidity, sunshine hours, and wind speed) were used in this study [62,63]. Digital elevation model (DEM) of 30 × 30 m resolution was used as input for the SWAT model. Based on the DEM data, sub-basins were delineated and other essential physiographic parameters were estimated for water resources estimation. The physiographic parameters include, among others, the hydrological response units (HRUs) and the slope. Moreover, soil and land use maps were used to analyze the irrigation potential of the basin [64].

Climate
The daily meteorological data including daily rainfall, temperature, wind speed, sunshine hours, and relative humidity of existing stations were obtained from the Ethiopian National Meteorological Service Agency (ENMSA). These data include stations that are within and around the DRB and are classified from class 1 to 4. The definition of each class is given on the ENMSA's official website (www.ethiomet.gov.et/). The time-series data of all stations have got continuous missing values, as indicated in the missing column in Table 1. The missing data can be filled automatically in the SWAT model using the weather generator. However, due to lack of completed datasets to create SWAT's weather generator, in this study the missing climatic data were filled manually using recorded data from the nearby stations. The latter helped to appropriately check the data quality and conduct the objective of the hydrological analysis of this study. For this, eight stations that fall into class 1 were screened based on their location to the study area (Table 1). Then missing values were filled based on the recorded values in the surrounding stations. This was completed by using quality controls description and performing direct replace or rainfall contour lines to fill the rainfall gap and average methods based on data from the surrounding stations. The locations of class 1 stations are shown in Figure 2.

Discharge Data
The streamflow data collected at Dhidhessa near the Arjo station by the Ministry of Water, Irrigation and Energy (MoWIE) features more extended time-series data (1989 to 2012). Visually gross data errors, such as missed erroneous and noises on recorded flows were checked to ensure the data quality.

Physiographical Data
Sub-basin parameters, such as the slope gradient, and slope length, and the stream network characteristics, such as the channel slope, length, and width, were derived from the DEM data ( Figure 2). The land use and soil parameters were derived from Figure 3 and used to identify potentially irrigable land and water resources. All the secondary data, such as the hydrological, meteorological and agricultural data were obtained from the Oromia Water Works Design and Supervision Enterprise (OWWDSE).
Hydrology 2020, 7, x FOR PEER REVIEW 7 of 23 Figure 3. Physiographic data (land use, land cover and major soil type).

SWAT Model
The water availability in the DRB was simulated using the SWAT model, which is a semidistributed, and physically based model [65]. This model has been widely used to evaluate land management practices, sediment transport, and nutrient loading and discharge at the basin scale. SWAT quantifies and characterizes the effect of land use management practices on water, agricultural

SWAT Model
The water availability in the DRB was simulated using the SWAT model, which is a semi-distributed, and physically based model [65]. This model has been widely used to evaluate land management practices, sediment transport, and nutrient loading and discharge at the basin scale. SWAT quantifies and characterizes the effect of land use management practices on water, agricultural nutrients, and sediment yields of the basin at different time spans [53,66,67]. In SWAT, a basin is sub-divided into sub-basins based on DEM [65]. The sub-basins are then further sub-divided into several hydrologic response units (HRUs) that represent a unique and homogeneous combination of slope, land use, and soil type within a sub-basin. SWAT computes the hydrological processes using the water balance concept and equations, applied at the HRU level of the model [65].

SWAT Sensitivity Analysis
The Sequential Uncertainty Fitting-II (SUFI-2) algorithm of the SWAT Calibration and Uncertainty Program (SWAT-CUP) was used for the sensitivity analysis (SA). SUFI-2 uses a One-factor-At-a-Time (OAT) global SA approach [68]. In this study, SA was only performed for discharge-related parameters ( Table 2). Daily observed data were used for SA. The sensitivity of a parameter was determined based on the t-stat and p-value. The smaller the p-value and larger the absolute value of the t-stat are, the more sensitive the parameter will be, whereas a larger p-value and smaller t-stat indicate less sensitivity for a given watershed. p-values close to zero have a greater significance. Accordingly, sensitivity analysis significantly lessens the relative sensitivity of parameter identification, increases the precision of calibration, and decreases the uncertainty and time necessary for calculations. We used the global sensitivity analysis option of the SUFI-2 algorithm, implemented in SWAT-CUP. The sensitivity of the parameters was computed to identify the most essential and sensitive parameters for the model calibration and validation processes, including uncertainty quantification [69]. The parameter sensitivity was ranked based on the t-statistics and level of significance at 5%. The 19 sensitivity parameters are presented in Table 2. All the selected parameters are considered sensitive and associated mainly with soil and groundwater, including EPCO, ESCO, CANMX, SOL_AWC, PLAPS, TLAPS, SLSOIL, SOL_DB, CH_K, CN2, SLSUBBSN, HRU_SLP, OV_N, ALPHA_BF, GWQMN, GW_DELAY, and GW_RECAP. The most sensitive parameters screened by SUFI-2 are enumerated in (Table 2).

Model Calibration and Validation
The SWAT model was set-up for the period from 1986 to 2012. While the first three years were used for model warm up, the periods from 1989-2000 and 2001-2012 were used for model calibration and validation, respectively.
The SUFI-2 algorithm, which is included in the SWAT-CUP interface, was used to perform the SWAT model calibration and validation processes [67]. The performance of the model was evaluated by graphical comparison and using statistical indices [53,70]. The statistical criteria include the Root Mean Square Error (RMSE), Mean Absolute Error (MAE), Coefficient of determination (R 2 ), Nash-Sutcliffe Coefficient (NSE), percent bias (PBIAS) [53,67,70]. The value of R 2 ranges from 0 to 1. As the value of R 2 approaches 1, the better is the performance of the model and values of R 2 that are less than 0.6 indicate poor performance of the model. The figure value of more than 0.6 is found in the satisfactory range for the coefficient of determination and a 0 value means no correlation at all; whereas 1 indicates that the desperation of the prediction is equal to that of the measured.
Additionally, we considered the model prediction uncertainty that is mapped into input model parameters. The parameter values are represented by uniform distributions. In SUFI-2, the model prediction uncertainty is quantified for the 95% prediction uncertainty (95PPU) interval. The 95PPU estimates simulated values at the 2.5% and 97.5% levels from hundreds of simulations. In this study, the 995PPU confidence interval was estimated for 500 simulations. The SUFI-2 algorithm introduces two efficiency criteria, P-factor and r-factor, that provide a measure of the model's ability to capture uncertainties and a measure of the quality of the calibration data [65,71].

Model Performance Evaluation
The results of the calibration and validation were evaluated using graphical and statistical procedures based on the quality criteria [53,70]. Daily data were used for the SWAT model, where the first two years were considered as the warmup period. Then, sequential Uncertainty Fitting (SUFI-2)-based calibrations were done for a period of 12 years starting from 1989 to 2000 whereas, validation was done for the 12 years from 2001 to 2012.
The quantitative statistical parameters employed for measuring the model quality are the Root Mean Square Error (RMSE), Mean Absolute Error (MAE), Coefficient of determination (R 2 ), Nash-Sutcliffe Coefficient (NSE), percent bias (PBIAS), p-factor, and r-factor [65,71]. The value of R 2 ranges from 0 to 1. The more the value of R 2 approaches 1, the better the performance of the model will be; conversely, values of R 2 less than 0.6 indicate poor performance of the model. A figure value of more than 0.6 represents a satisfactory range for the coefficient of determination, while a value of 0 means no correlation at all; conversely, 1 indicates that the determination of the prediction is equal to that of the measured value.
where i is the i th measured or simulated datum, Q s is the simulated value, Q m is the observed value, Q m is the mean measured value and Q s is the mean simulated value. The Nash Sutcliffe Efficiency (NSE) was used to measure the relative magnitude of the residual variance compared to the measured data variance and define the extent of the variation between the simulated and observed data discrepancies, and its value ranges from −∞ to 1.

Assessing Water Availability
To analyze the irrigation potential of the DRB, the soil map, land use, and land cover were used as inputs for the suitability assessment. The DEM data were processed for the extraction of flow direction, flow accumulation, stream network generation and delineation of the watershed and sub-basins using ArcGIS (Figure 4). Information on low flow is required to estimate the amount of water available for surface water irrigation application during the dry season [21]. Low-flow characteristics were estimated using a flow duration curve (FDC) and by determining the 70, 80, 85 and 90-percentile available flow values. The latter is described as the flow exceedance at 70% (Q70), 80% (Q80), 85% (Q85), and 90% (Q90) of the time. A FDC at a particular outlet on a river shows the proportion of time during which the discharge at the location equals or exceeds specific values. The FDCs for long periods of runoff are useful for deciding what proportion of flow should be used for particular purposes since the area under the curve represents the volume.
The land use and soil map along with their respective lookup tables (for reclassification of according to SWAT coding conventions) were prepared as inputs for the SWAT model. The ArcGIS platform of SWAT, which referred as "ArcSWAT", provides the user with a complete set of GIS tools for developing, running, and editing the hydrologic and management inputs of the model. The spatially distributed data required for Arc SWAT includes the Digital Elevation Model (DEM), soil and land use maps, either as shapefiles or as grid format (Figures 4 and 5).  The land use and soil map along with their respective lookup tables (for reclassification of according to SWAT coding conventions) were prepared as inputs for the SWAT model. The ArcGIS platform of SWAT, which referred as "ArcSWAT", provides the user with a complete set of GIS tools for developing, running, and editing the hydrologic and management inputs of the model. The spatially distributed data required for Arc SWAT includes the Digital Elevation Model (DEM), soil and land use maps, either as shapefiles or as grid format (Figures 4 and 5).

Computation of the Irrigation Potential Area
The irrigation potential area of the DRB was determined using the volume of water that could be supplied at varying levels of risk from the major river during the irrigation period of the selected crops (cabbage, maize, tomato, pepper, groundnut and sugarcane). Figure 6 summarizes the overall methodology used to analyze the irrigation water potential and water requirement of the upper DRB.

Computation of the Irrigation Potential Area
The irrigation potential area of the DRB was determined using the volume of water that could be supplied at varying levels of risk from the major river during the irrigation period of the selected crops (cabbage, maize, tomato, pepper, groundnut and sugarcane). Figure 6 summarizes the overall methodology used to analyze the irrigation water potential and water requirement of the upper DRB.
The gross irrigation water requirement (GIWR) is the amount of water to be extracted (by diversion or pumping) and applied to the irrigation command and scheme. The potential irrigable area was calculated as an irrigation water supply (IWS) from the river during the irrigation period divided by the GIWR. Thus, the possible irrigable area of the basin is computed based on the following. PIA = IWS/GIWR (6) where PIA is the Potential Irrigable Area (ha); IWS is the Irrigation Water Supply from the river during the irrigation period (m 3 ), and GIWR is the Gross Irrigation Water Requirement (m 3 /ha). Hydrology 2020, 7, x FOR PEER REVIEW 13 of 23 The gross irrigation water requirement (GIWR) is the amount of water to be extracted (by diversion or pumping) and applied to the irrigation command and scheme. The potential irrigable area was calculated as an irrigation water supply (IWS) from the river during the irrigation period divided by the GIWR. Thus, the possible irrigable area of the basin is computed based on the following.

PIA=IWS/GIWR (6)
Area (ha); IWS is the Irrigation Water Supply from the river during the irrigation period (m 3 ), and GIWR is the Gross Irrigation Water Requirement (m 3 /ha). The irrigation potential area is the minimum area that can be irrigated with the available water supply (m 3 ) divided by the annual irrigation water demand (m 3 /ha). This value indicates the areas that can be irrigated with the available streamflow of the sub-basin. The irrigation potential of a given sub-basin was obtained by comparing the irrigation requirements of the identified land suitable for surface irrigation system and the available mean monthly flows in the sub-basin's river based on the suggested method [46].
To estimate the GIWR of a given crop, we used the climatic data from three different stations (Arjo, Bedelle and Angar), which represent each climatic zone of the basin. Some specific procedures recommended by the FAO for estimating missing climatic data [72,73] were applied in this study. The irrigation potential assessment mainly based on the irrigation water requirement, which was computed using the CROPWAT model depending on the cropping pattern and climate of the area The irrigation potential area is the minimum area that can be irrigated with the available water supply (m 3 ) divided by the annual irrigation water demand (m 3 /ha). This value indicates the areas that can be irrigated with the available streamflow of the sub-basin. The irrigation potential of a given sub-basin was obtained by comparing the irrigation requirements of the identified land suitable for surface irrigation system and the available mean monthly flows in the sub-basin's river based on the suggested method [46].
To estimate the GIWR of a given crop, we used the climatic data from three different stations (Arjo, Bedelle and Angar), which represent each climatic zone of the basin. Some specific procedures recommended by the FAO for estimating missing climatic data [72,73] were applied in this study. The irrigation potential assessment mainly based on the irrigation water requirement, which was computed using the CROPWAT model depending on the cropping pattern and climate of the area [46,72,73]. The CROPWAT model was used to estimate the reference evapotranspiration (ETo) and net irrigation water requirement (NIWR). To facilitate the analysis, CROPWAT software subdivides the net irrigable area into different crop units with different planting dates at intervals of 2 to 10 days (Figure 7). The length of the growing season is the period (in days) during a year when precipitation exceeds half the evapotranspiration (i.e., PCP > 0.5 ETo). The normal growing period (also called a type three season) is characterized by a dry period, moist period (also called the intermediate period), and wet period (or also referred as the humid period). net irrigation water requirement (NIWR). To facilitate the analysis, CROPWAT software subdivides the net irrigable area into different crop units with different planting dates at intervals of 2 to 10 days (Figure 7). The length of the growing season is the period (in days) during a year when precipitation exceeds half the evapotranspiration (i.e., PCP > 0.5 ETo). The normal growing period (also called a type three season) is characterized by a dry period, moist period (also called the intermediate period), and wet period (or also referred as the humid period). Irrigation efficiency is a fraction of the gross irrigation water depth-that is, the net irrigation water depth needed to effectively reach the crop root zone. Irrigation water is required when rainfall is insufficient to compensate for the water lost by evapotranspiration. By calculating the soil water balance of the root zone on a daily basis, the timing and depth of future irrigations can be planned. The net irrigation water requirement was estimated using CROPWAT8. The gross irrigation water requirements at the high land, midland, and low land for different water application methods (surface, sprinkler, and drip irrigation methods) were calculated based on the net irrigation requirements.

Land Suitability Evaluation
This study considers the irrigation potential and offers feasible alternatives for sustainable land use and productive, without negatively impacting natural resources. To determine the irrigation potential, soil factors such as drainage, depth, texture, slope, and nutrient availability (such as pH) were evaluated. The soils in the study area have remarkable potential for a variety of rain-fed and irrigated crops. Hence the texture, depth, drainage, and pH of the soil raster, land use raster, land cover raster, and slope raster maps were generated. Moreover, a suitability model was created using the GIS model builder in the Arc tools box and the tools from the spatial analysis tool sets. Then, after individual suitability was assessed, the irrigation suitability factors considered in this study, such as the slope factor, soil factor, and land cover /use factor, were used as the input for the irrigation Irrigation efficiency is a fraction of the gross irrigation water depth-that is, the net irrigation water depth needed to effectively reach the crop root zone. Irrigation water is required when rainfall is insufficient to compensate for the water lost by evapotranspiration. By calculating the soil water balance of the root zone on a daily basis, the timing and depth of future irrigations can be planned. The net irrigation water requirement was estimated using CROPWAT8. The gross irrigation water requirements at the high land, midland, and low land for different water application methods (surface, sprinkler, and drip irrigation methods) were calculated based on the net irrigation requirements.

Land Suitability Evaluation
This study considers the irrigation potential and offers feasible alternatives for sustainable land use and productive, without negatively impacting natural resources. To determine the irrigation potential, soil factors such as drainage, depth, texture, slope, and nutrient availability (such as pH) were evaluated. The soils in the study area have remarkable potential for a variety of rain-fed and irrigated crops. Hence the texture, depth, drainage, and pH of the soil raster, land use raster, land cover raster, and slope raster maps were generated. Moreover, a suitability model was created using the GIS model builder in the Arc tools box and the tools from the spatial analysis tool sets. Then, after individual suitability was assessed, the irrigation suitability factors considered in this study, such as the slope factor, soil factor, and land cover /use factor, were used as the input for the irrigation suitability model to determine the most suitable land for surface irrigation. The three layers (slope, soil, and land cover raster) were overlaid to find the final and suitable land.

Surface Water Modelling
To assess irrigation water availability, the SWAT model was used. The sensitivity of the SWAT model parameters was determined for the flow-related parameters [69]. Before estimating the available water for irrigation, we calibrated and validated the SWAT model against observed stream flows. Table 3 summarizes the model performance using multiple statistical criteria. The values of the statistical indices of R 2 and NSE were 0.85 and 0.87 for the calibration period ( Figure 8) and 0.91 and 0.89 for the validation period, respectively ( Figure 9). These values indicate the "good" performance of the model in simulating the hydrology of the basin [53,71]. The p-factor and r-factor were 0.77 and 0.75 for the calibration period and 0.87 and 0.81 for the validation period, respectively. SUFI-2 was found to be more convenient and simpler to use compared to other automatic calibration techniques. The SWAT model performance indicates good performance that fulfills the necessary criteria [53,70]. The Dhidhessa catchment satisfies the following criteria ( Table 3): The coefficient of determination is best when R 2 > 0.60, and the Nash-Sutcliffe (1970) coefficient of efficiency is best when NS > 0.5. The average annual surface water potential of the upper Dhidhessa River Basin was estimated from the SWAT calibrated and validated results. It was found that the average annual precipitation and water yield of the basin are 1787 and 1605 mm of the water yield, surface water runoff contributed 146 mm, while the base flow accounts for 320 mm. Considering the basin's area of 14,710.41 km 2 , it was estimated that 2147 million m 3 of water is delivered from the basin to the river, indicating about 9. 26 BCM is annually available for surface irrigation application.   The average annual surface water potential of the upper Dhidhessa River Basin was estimated from the SWAT calibrated and validated results. It was found that the average annual precipitation and water yield of the basin are 1787 and 1605 mm of the water yield, surface water runoff contributed 146 mm, while the base flow accounts for 320 mm. Considering the basin's area of 14,710.41 km 2 , it was estimated that 2147 million m 3 of water is delivered from the basin to the river, indicating about 9. 26 BCM is annually available for surface irrigation application.

Surface Water Availablity Assessment
To determine how much the available water can be reliably abstracted from the rivers for irrigation, flow duration analyses were carried out at a selected outlet on the river and using mean monthly simulated flow time series. The amounts of water for the different dependable flow values that can be reliably abstracted were estimated for each sub-watershed outlet.
The flows were determined by ranking all the mean monthly discharges. The purpose of the flow duration curve (FDC) analysis was to evaluate the dependable flow at various percentages for surface irrigation water use. The FDC gives the duration of occurrences in the whole range of the flow values of the river. The result indicates that the annual dependable flow values at 70% and 80% are 7.56 and 6.97 BCM, respectively (Table 4).

Potential Irrigable Area
The analyses revealed that the monthly irrigation requirements are less than the available mean monthly flow values. In addition, these mean monthly flows are slightly higher than the irrigation water requirements of the crops in their corresponding command areas.
The result of the annually available surface water in the sub-basin is 9.26 BCM. Based on the annual water supply and demand, the available water is capable of irrigating the available land suitable for agriculture. The result indicates that there is an excess amount of water resources available to satisfy the irrigation demand. Moreover, the annual irrigation potential, with respect to the physical resources and irrigation water requirements for different slope classes and dependable flow values, are shown in Table 5. The irrigation potential area was estimated as the quotient of the monthly available flow and the monthly gross irrigation water requirements during the growing season of the crop. The result indicates that the potential irrigation area for the sub-basin was greater than the potentially irrigable land identified using physiographic factors for all irrigation application methods. This study also indicates that the assessed irrigation potential of the basin is higher compared to that of the Tana basin of the BNR [46].

Land Suitability Assessment Results
The land suitability shows that as the maximum slope limit increases, the suitable land area for irrigation also increases ( Table 6). This analysis gives useful insights into the suitable land resources development for irrigation at different slope levels. Overall, there is a large amount of land resources and potential irrigable area (1,023,581.08 ha (69.58%)) available for irrigation development and applications of the DRB.

Conclusions and Recommendations
In this paper, a SWAT model was developed to assess the available water resources of the Dhidhessa River Basin (DRB) for the potential application of surface water irrigation and improvement of agricultural crop productivity. The SWAT simulated flows were further analyzed using flow duration curves whereby 70% and 80% dependable flow values were estimated.
The study also estimated the irrigation water requirement of six selected crops (cabbage, maize, tomato, pepper, groundnut and sugarcane) and successfully identified their irrigation command areas. This study indicates that there is enough water to expand irrigable areas and improve the livelihoods of the local communities through implementing an enhanced irrigation system and improving agricultural crop productivity. However, improvement in irrigation efficiency depends on the choice irrigation type that plays a critical role in the water resource development and management of the Dhidhessa River. This could have an effect on satisfying the downstream environmental flow needs and water users. Therefore, to avoid more water loses in the irrigation field and to improve the irrigation water efficiency, it is highly recommended to adopt sprinkler or drip irrigation methods. This is expected to increase the irrigable area especially for command areas located in slopes ≤ 15%. This study further noticed that the identified surface irrigation potential of the DRB can assist in the decision making process and implementation of irrigation development projects in the basin. Finally, the application of remote sensing and GIS was found to be helpful in assessing the surface irrigation potential of the study area. Therefore, future irrigation development and activity would exploit these updated resources and combined with the developed SWAT model for better assessment of the land and water resources the study area and elsewhere.