Assessing the Water-Resources Potential and Soil Erosion Hotspot Areas for Sustainable Land Management in the Gidabo Watershed, Rift Valley Lake Basin of Ethiopia

For development of a comprehensive sediment management plan, it is crucial to categorize watersheds on the basis of soil erosion hotspot areas to extend the useful life of water bodies (e.g., Gidam reservoir). The goal of this study was to assess the surface water potential and identify erosion hotspot areas of the Gidabo watershed in Ethiopia using the Soil and Water Assessment Tool (SWAT) model. The SUFI-2 (Sequential Uncertainty Fitting Version 2) program was used to calibrate the model, and the model’s performance was evaluated. According to the catchment prioritization analysis, some of the sub-basins with similar land use, land cover, and soil type but with higher slope would generate higher sediment yield. Furthermore, the soil conservation scenarios were developed in SWAT, and the model result showed that average annual sediment yield could be reduced by the application of grassed waterway, filter strips, terracing, and contouring by 49%, 37.53%, 62.32%, and 54.6% respectively. It was concluded that sediment yield reduction by applying terracing was more effective than other conservation measures for affected sub-basins. The surface water potential of the watershed varies spatially from sub-basin to sub-basin, and the mean monthly surface water potential of the watershed is 33 million cubic meters. These findings can help decision-makers to develop appropriate strategies to minimize the erosion rate from erosion hotspot areas and to allocate the watershed water potential for different types of water demands. Strip planting, terracing, or contour farming may be necessary on chosen hotspot erosion sites to reduce the effect of slopes on surface runoff flow velocity and sediment transport capacity.


Introduction
Water is a vital and limited natural resource. Due to the increasing water need throughout the world, the accessibility of freshwater in numerous areas is decreasing due to the effect of the growing population, industrial development, high ET rate due to warming, and land-use land-cover changes [1]. Land-use land-cover change plays a major role in the availability of freshwater and is just one of the challenges currently facing water resources. The rapidly growing population together with the high dependence of the economy on agriculture contributes to increased demand and increased competition for inadequate water resources. Hence, quantifying the water resources and sediment yield of the basin is essential for the planning and management of available land water resources [2].
In a similar manner, estimating the sediment yields of the river basin is critical for evaluating the effectiveness of catchment management strategies in order evaluate the impacts of deposition of sediment on water bodies [3]. The shortage of information on the degree of erosion and sediment yields, as well as the transportation of fine sediment, impedes the design and application of better watershed management plans [4][5][6]. Sediment yield provides an important index of land degradation, severity, and trends, and it also Sustainability 2022, 14, 5262 3 of 19 in halting soil erosion [41,43], the method is hampered by the lack of intervention prioritization strategies that identify locations extremely vulnerable to erosion. Today, Ethiopia is struggling with poverty. The implementation of irrigation by constructing water storage dams is the main approach to reducing poverty by increasing productivity. However, our dams are facing a serious problem with sedimentation. Thus, effective modeling of sediment yield and accurate determination of the water potential of the basin are crucial steps for developing watershed management plans. The Gidabo river basin is a central highland basin in Ethiopia's Rift Valley Region, where soil erosion is widespread [44]. To address this, it is necessary to identify the basin's most erosion-prone locations so that effective conservation measures can be implemented.
To achieve the overall goal of the study, the following objectives were addressed: • To estimate the current water budget of the Gidabo watershed; • To identify and prioritize erosion hotspot sub-catchments on the basis of estimated runoff and sediment yield; • To propose appropriate watershed management possibilities for reducing soil degradation complications. particular and reverse land deterioration in general [40,43]. Although significant progress has been made in halting soil erosion [41,43], the method is hampered by the lack of intervention prioritization strategies that identify locations extremely vulnerable to erosion. Today, Ethiopia is struggling with poverty. The implementation of irrigation by constructing water storage dams is the main approach to reducing poverty by increasing productivity. However, our dams are facing a serious problem with sedimentation. Thus, effective modeling of sediment yield and accurate determination of the water potential of the basin are crucial steps for developing watershed management plans. The Gidabo river basin is a central highland basin in Ethiopia's Rift Valley Region, where soil erosion is widespread [44]. To address this, it is necessary to identify the basin's most erosion-prone locations so that effective conservation measures can be implemented. To achieve the overall goal of the study, the following objectives were addressed:

Materials and Methods
• To estimate the current water budget of the Gidabo watershed; • To identify and prioritize erosion hotspot sub-catchments on the basis of estimated runoff and sediment yield; • To propose appropriate watershed management possibilities for reducing soil degradation complications.

Location and Topography
The Gidabo watershed is situated in the Rift Valley Basin of Ethiopia between 6°20′ and 7° N and between 38°05′ and 38°38′ E ( Figure 1) with a coverage area of 3386 square kilometers.   The altitude of the catchment area ranges from 1183 to 3182 m above mean sea level (a.m.s.l). The major physiographic units in this area are undulating plains, valleys, steep stream banks, and major escarpments in the eastern region of the watershed.

Climate
There are rainy and dry seasons in the Gidabo watershed. The main rainy season is from April to October with a peak rainy season from April to May and a second peak rainy season from September to October. The mean monthly rainfall of the watershed varies between 37 mm and 188 mm, and the mean annual rainfall varies from 1219 mm to 1593 mm ( Figure 2). The monthly average maximum temperature of the watershed varies between 24.2 • C and 33.06 • C; similarly, the average monthly minimum temperature of the watershed varies between 10.7 • C and 17.43 • C.

Land Use and Land Cover
Agriculture is the dominant land-use type inside the basin (Figure 3a). The entire basin is home to comprehensive farming practices, and various crops are grown during the two known rainy seasons, locally named Kiremt (June to September) and Belgi (February to May in Ethiopia). An assessment conducted by Ethiopia [44] indicated dynamic landuse and land-cover changes inside the watershed. The assessment reported the presence of comprehensive farming systems extending to peripheral lands, suffering from weak agricultural technology together with poor catchment management. This leads to the occurrence of high-degree soil erosion, a decrease in soil fertility, and land degradation as a whole.

Soil and Geology
The dominant soil in the watershed includes Eutric Vertisols with sandy loam in the middle sub-basin, Orthic Luvisols and Chromic Luvisols with clay to sandy loam in the upper sub-basin, Pellic Vertisols, Dystric Nitosols, and Eutric Cambisols in the upper and middle sub-basins [44] (Figure 3b). Geologically, the southern region predominantly features volcanic quaternary rhyolites and trachyte, while the western region features Oligocene to Miocene basalts. The lowlands of the watershed are covered with Holocene alluvial and Eolian deposits of Rift Valley [44].

Hydrology
The Gidabo basin with a total area of 3384 square kilometers is located in the Rift Valley Lake basin, particularly in the Abaya Lake sub-basin. The river is composed of various streams initiating from the highland regions of the basin. The river flows from the northeast and crosses the highway of Addis Abeba to Dilla at Aposto Town where it is gauged. The average flow of the river at the Aposto station is 17.35 m 3 /s [44].

River Basin Water Balance
Water balance is the driving force behind all processes in the Soil and Water Assessment Tool (SWAT) model [45]. SWAT simulates the hydrologic cycle as a function of the principle of water balance [46].
where SW t is the soil water content (mm), SW o is the initial water content, t is the time in (days), R day . is the rainfall amount on day I (mm), Q surf . is the surface runoff quantity on day I (mm), E a is the evapotranspiration amount on day I (mm), W seep . is the quantity of water entering the groundwater from the soil on day I (mm), and Q gw is the quantity of return flow on day I (mm). from April to October with a peak rainy season from April to May and a second peak rainy season from September to October. The mean monthly rainfall of the watershed varies between 37 mm and 188 mm, and the mean annual rainfall varies from 1219 mm to 1593 mm ( Figure 2). The monthly average maximum temperature of the watershed varies between 24.2 °C and 33.06 °C; similarly, the average monthly minimum temperature of the watershed varies between 10.7 °C and 17.43 °C. the two known rainy seasons, locally named Kiremt (June to September) and Belgi (February to May in Ethiopia). An assessment conducted by Ethiopia [44] indicated dynamic land-use and land-cover changes inside the watershed. The assessment reported the presence of comprehensive farming systems extending to peripheral lands, suffering from weak agricultural technology together with poor catchment management. This leads to the occurrence of high-degree soil erosion, a decrease in soil fertility, and land degradation as a whole.

Soil and Geology
The dominant soil in the watershed includes Eutric Vertisols with sandy loam in the middle sub-basin, Orthic Luvisols and Chromic Luvisols with clay to sandy loam in the upper sub-basin, Pellic Vertisols, Dystric Nitosols, and Eutric Cambisols in the upper and middle sub-basins [44] (Figure 3b). Geologically, the southern region predominantly features volcanic quaternary rhyolites and trachyte, while the western region features Oligocene to Miocene basalts. The lowlands of the watershed are covered with Holocene alluvial and Eolian deposits of Rift Valley [44].

Hydrology
The Gidabo basin with a total area of 3384 square kilometers is located in the Rift Valley Lake basin, particularly in the Abaya Lake sub-basin. The river is composed of various streams initiating from the highland regions of the basin. The river flows from the northeast and crosses the highway of Addis Abeba to Dilla at Aposto Town where it is gauged. The average flow of the river at the Aposto station is 17.35 m 3 /s [44]. In the watershed, the SWAT model simulates runoff using the curve number method of the Soil Conservation Service (SCS) [47]. It estimates the surface runoff using the following equation: where Q surf is the accumulated runoff or rainfall excess (mm), R day is the daily rainfall depth (mm water), I a is the initial abstraction before the runoff event (mm water), and S is the retention parameter (mm water). The variation in land surface features such as land use, slope, and management measures causes spatial variation in retention parameters. This can be expressed mathematically as where CN is the daily curve number, which is governed mainly by land use, soil permeability, and a hydraulic group of soils. The initial abstraction I a is mostly taken as 0.2 S, and Equation (2) can be derived as The SWAT model predicts the retention parameter with the help of two methods; the first method of calculation uses the soil moisture content in its profile, in which runoff is overestimated in shallow soil. CN is mainly governed by the antecedent climate, and the value is rarely governed by the soil storage. The CN estimation is based on evapo- transpiration. The other method predicts the retention parameter using the accumulated plant evapotranspiration.
where S is the retention parameter for a specified day (mm), S max is the maximum retention parameter on any specified day (mm), SW is the moisture content of the soil excluding the water held at wilting point (mm), and w 1 and w 2 are shape coefficients. When evapotranspiration governs the retention parameter, the result of the retention parameter at the end of the day can be updated as where S is the retention parameter for a given day (mm), S prev is the previous day's retention parameter (mm), E o is the daily potential evapotranspiration (mm per day), cncoef is the average coefficient to estimate the retention coefficient for the daily curve number predictions, S max is the daily higher retention parameter value (mm), R day is the daily rainfall depth (mm), and Q surf is the surface runoff (mm). Evapotranspiration (PET) is one of the parameters for studying the water balance. In the SWAT model, three evapotranspiration prediction methods (the Penman-Monteith method [48], the Priestley-Taylor method [49], and the Hargreaves method [50]) are integrated.
The Penman-Monteith method requires solar radiation, air temperature, relative humidity, and wind sped; the Priestley-Taylor method requires solar radiation, air temperature, and relative humidity; the Hargreaves method requires air temperature only. For this study, the Penman-Monteith method was used to estimate evapotranspiration.
The Penman-Monteith equation that estimates evapotranspiration is as follows: where ET is the daily reference crop evapotranspiration (mm·day −1 ), R net is the net radiation flux (MJ·m −2 ·day −1 ), G is the heat flux density in the soil (which is very small and can be neglected) (MJ·m −2 ·day −1 ), T is the mean daily air temperature ( • C), γ is the psychometric constant (kPa· • C −1 ), U. is the wind speed measured at 2 m height (m·s −1 ), e s is the saturation vapor pressure (e a = e s × RH/100 (kPa)), RH is the relative humidity (%), and ∆. is the slope of the saturation vapor pressure curve (a· • C −1 ). In a similar manner, the model estimates the ground water flow using Equation (8).
where aq sh,i is the water accumulated in the shallow aquifer on day I (mm), aq sh,i−1 . is the water accumulated in the shallow aquifer on the previous day (mm), W rchrg is the recharge percolating to the aquifer on day I (mm), Q revap is the water entering into the soil zone on day I (mm), W deep is the water percolating from the shallow to the deep aquifer on day I (mm), and W pump,sh is the water removed from the shallow aquifer by pumping on day I (mm).

Sediment Balance
SWAT predicts soil erosion and sedimentation by applying the Modified Universal Soil Loss Equation (MUSLE) [51]. In MUSLE, the rainfall energy factor used for USLE is replaced with a runoff factor to improve the sediment yield prediction, which eliminates the need for delivery ratios and allows the equation to be applied to individual storm events. Delivery ratios are not needed with MUSLE because the runoff factor represents energy used in detaching and transporting sediment. (9) where Sed is the sediment yield on a given day (metric tons), Q sur f is the surface runoff volume (mm /ha), Q peak is the peak runoff rate (m 3 /s), area hru is the area of the HRU (ha), K USLE is the soil erodibility factor, C USLE is the cover and management factor, P USLE is the support practice factor, LS USLE is the topographic factor, CFRG is the coarse fragment factor, and 0.56 is the delivery ratio.
SWAT simulates the sediment flow in channel networks using the following formula [52]: where Sed ch is the suspended sediment available in the reach, Sed ch,i is the suspended sediment provided in the reach at the beginning of the time period, Sed dep is the deposited sediment in the reach segment, and Sed deg is sediment amount reentering the reach segment. According to [52], SWAT simulates the quantity of sediment moved out of the channel as follows: where Sed out is the amount of sediment moved out of the channel, Sed ch is the suspended sediment available in the reach, V out is the outflow volume, and V ch is the volume in the reach segment.

Model Development and Input Description
Two datasets describing watershed features and meteorological data are required to develop the SWAT model.
Land-use/cover data were obtained from Landsat-8 OLI (Operational Land Imager); for this study, Landsat-8 images from 2015 were used for land-use/cover classification. The satellite images of the watershed with high resolution and zero cloud cover were collected from the US Geological Survey (USGS) Center for Earth Resources Observation and Science (EROS) via https://earthexplorer.usgs.gov (accessed on 5 January 2015). ER-DAS IMAGINE was used for image classification purposes, and Arc GIS10.3 was used for mapping purposes.
In the processes of image classification, the main task is assigning pixels of a constant raster images to the predefined land-cover classes. In this study, supervised classification was applied. This is the most common type of classification method in which all pixels with comparable spectral values are automatically categorized into land-cover classes. Supervised classification, which depends on the prior knowledge of pattern recognition of the research area, was used. For this study, the land-cover map was produced depending on the pixel-based supervised classification. For the land-use/cover classification, a general correctness assessment and kappa coefficient of 95.7% and 0.9402, respectively, were realized, in agreement with the minimum accuracy level (85%) recommended by [53].
The essential soil properties needed to set up the Soil and Water Assessment Tool Model are soil texture, size percentage, soil saturated hydraulic conductivity, grain, bulk density, texture class, and soil available water. These soil characteristics were obtained from laboratory analysis presented in the Rift Valley Lake Basin Master Plan study document [44]. During the study of the Ethiopian Rift Valley Lake Basin Master Plan, soil samples were collected from all soil units of the basin, and physical and chemical laboratory analyses were conducted in the Ethiopia Water Works Design and Supervision Enterprise (WWDSE) laboratory. From 13 soil units in the basin, 203 soil samples were collected, and their physical and chemical properties were analyzed. Hence, the soil database of the SWAT model was set up for the basin using the analyzed soil properties. The basin's soil erodibility (K) factor was calculated using the equation shown in EPIC [54] from the analyzed soil parameters.
The meteorological data included daily precipitation, maximum and minimum temperature, daily wind speed, daily sunshine hours, and daily relative humidity, which were obtained from meteorological stations (Table 1) available within and nearby the study area. Daily data over 27 years (1991-2017) were collected for the study. To fill the missing values of climate elements, a weather generator model was used. The required statistical parameters were computed using the computer program developed by [55]. As shown in Table 1, in the river basin, only one meteorological station (Dilla) was used to establish the weather generator database.

The Hydrological Data
In the Gidabo river basin (total area of 3386 square kilometers), there are four stream gauging stations, namely, Aposto (area of 703 square kilometers), Kola (area of 145 square kilometers), Badessa (area of 76 square kilometers), and Measso (area of 2462 square kilometers) with data periods of 1997-2015 for Aposto, Kolla, and Badessa, whereas Measso has stream flow data from 1997-2006. As shown in Figure 4, the gauging station Measso is located near to the outlet of the watershed, and it was selected for model calibration and validation. The daily flow data of all stations were collected from the Ministry of Water and Energy of Ethiopia. For all stations (except Measso), daily flow data were available up the year 2015, and a regression model (r 2 = 0.85) was used to generate the data of Measso from its upper station Aposto.
As shown in the Figure 4, 27.3% of the catchment was ungauged. To predict the water and sediment budget of the ungauged regions, an empirical regression model relating sediment yield and flow of the gauged stations to several catchment characteristics, namely, drainage area, slope, and average annual rainfall, was used. Three explanatory factors were calculated for gauged and ungauged river catchments, i.e., areas of catchments and slopes were processed from the digital elevation model (DEM), and their mean annual area rainfall was calculated on the basis of the inverse distance-weighted interpolation (IDW) of the nearby meteorological stations (Figure 2).

Sub-Catchment Delineation
Using the 30 × 30 m resolution DEM data, the contributing upstream area was delineated using ArcGIS 10.3 (Esri, Redlands, CA, USA). Accordingly, the entire study area was divided into 23 sub-watersheds. According to the model setup, the lowest area threshold values for land use, soil, and slope were set as 5%, 10%, and 10%, respectively, and 105 hydrological response units (HRUs) were identified, denoting unique combinations of land use, soil type, and slope.

Model Calibration and Validation
For the purposes of calibration and validation of the SWAT model, the daily discharge data of the river were obtained from the MoWE of Ethiopia, while the sediment data were generated using the sediment discharge rating curve from the Rift Valley Lake Basin Master Plan study [44] with an R 2 value of 0.8. For calibration and validation of the model, the bed load component of the sediment was considered suspended as the model simulates total sediment load plus bed load, which was commonly overlooked in most studies in the country due to measurement limitations. In most streams, bed load to suspended load varies from 10% to 30% [56]. For instance, a study in the Ziway Lake Basin in Ethiopia [3] used 10% of the suspended sediment as bed load. Since Lake Ziway and Gidabo River are located in the Ethiopian Rift Valley, we assumed the bed load as 10% of the suspended sediment load gained from the rating curve, which was included for model calibration.
(Dilla) was used to establish the weather generator database.

The Hydrological Data
In the Gidabo river basin (total area of 3386 square kilometers), there are four stream gauging stations, namely, Aposto (area of 703 square kilometers), Kola (area of 145 square kilometers), Badessa (area of 76 square kilometers), and Measso (area of 2462 square kilometers) with data periods of 1997-2015 for Aposto, Kolla, and Badessa, whereas Measso has stream flow data from 1997-2006. As shown in Figure 4, the gauging station Measso is located near to the outlet of the watershed, and it was selected for model calibration and validation. The daily flow data of all stations were collected from the Ministry of Water and Energy of Ethiopia. For all stations (except Measso), daily flow data were available up the year 2015, and a regression model (r 2 = 0.85) was used to generate the data of Measso from its upper station Aposto. As shown in the Figure 4, 27.3% of the catchment was ungauged. To predict the water and sediment budget of the ungauged regions, an empirical regression model relating sediment yield and flow of the gauged stations to several catchment characteristics, namely, drainage area, slope, and average annual rainfall, was used. Three explanatory factors were calculated for gauged and ungauged river catchments, i.e., areas of catchments and slopes were processed from the digital elevation model (DEM), and their mean annual area rainfall was calculated on the basis of the inverse distance-weighted interpolation (IDW) of the nearby meteorological stations (Figure 2). The SWAT model was calibrated and validated to ensure its ability to prioritize the sub-watershed according to sediment contribution. The model was run for the simulation period of 1 January 1998 through December 2015 for the Measso gauging station. River discharge and sediment data over 10 years from 2001 to 2010 were used for calibration, and data from the subsequent 5 years (2011-2015) were then used for the validation period. The first 3 years in the calibration run were used for model preparation. For this study, the duration lengths used for calibration and validation were fixed according to the observed data records. The sensitivity analyses were conducted automatically using the SUFI-2 program in SWATCUP software during calibration.

Quality Control and Statistical Analysis of Data
Checking for continuity and consistency of both metrological and hydrological data is required before using them for further analysis. Quality assessment and filling of data gaps are needed for the analyses. Continuous daily series data are exposed to outliers caused by instrumental and/or human error. To perceive these outliers, concentrations less than half or more than twice the expected amount were removed from further analysis according to Grubb's test method [57]. No more than 1% of the dataset was rejected in all cases. Four broadly used statistical indices, namely, the Nash-Sutcliffe coefficient (NSE), percentage bias (PBIAS), coefficient of determination (R 2 ), and RSR (ratio of the root-mean-square error (RMSE)) were used to evaluate model performance by fitting the measured constituent data to the standard deviation of the measured data [58].

Analysis of Sensitivity
To find the most sensitive model parameter in the Gidabo watershed, a stream flow sensitivity analysis was carried out after simulating the model using daily data. A total of 10 parameters were used for sensitivity analysis, and the parameters governing runoff in the Gidabo watershed were ranked as high, medium, and low according to the associated p-value and equivalent t-stat value; the results are shown in Table 2.

Stream Flow Calibration and Validation
The six most influential flow parameters with high and medium sensitivity ( Table 2) were used for further iterations in the calibration period (Table 3). The model performance evaluation result is presented in Table 4 according to [58]. The result of flow calibration was very good. In a similar manner, for the validation period, the model performance statistics were determined, revealing good results (Table 4). Moreover, the graphical representations of measured and simulated stream flows correlated well for both calibration and validation periods ( Figure 5). As shown in Figure 6, the calibration of the SWAT model was successfully achieved using the observed daily sediment flow. Similarly, during the validation stage, the SWAT model adequately mimicked the sediment inflow of the river.  According to the outcome of calibration and validation, the model performance was examined according to statistical indicators (Table 7).

Water Balance Components
The Gidabo watershed annual water balance components were simulated by the SWAT model using the principle of conservation of mass. The simulated annual water balance components showed that an average annual precipitation of 1400.4 mm and evaporation loss of 652.7 mm in the basin, accounting for 49.95% of the annual water budget. Surface runoff (SURQ), lateral flow (LATQ), and ground water flow (GWQ) were 342.24 mm (19.58%), 275.6 mm (8.17%), and 21.35 mm (6.82%), respectively. The average stream flow (WYLD) is a combination of surface runoff, lateral, and ground water flow, and it accounted for 421.63 mm (34.57%) of the annual water budget. The deep percolation (PERCOLATE) and initial soil water content (SW o ) accounted for 107.10 mm (8.78%) and 167.14 mm (13.71%), respectively, of the annual water budget. The final soil water content (SW t ) determined using the water balance equation was 81.72 mm.
The water yield of ungauged regions was estimated using an empirical regression model relating the flow of gauged stations (Aposto, Kolla, Badessa and Measso) to several catchment characteristics, namely, drainage area, slope, and average annual rainfall. As a result, the average annual flow rate of ungauged regions was determined as 2201.714 CMC.

Sensitivity Analysis
Twelve sediment parameters were checked using SUFI-2 in SWAT-CUP to identify the most sensitive parameters (Table 5).

Sediment Yield Calibration and Validation
After stream flow calibration, sediment flow calibration was conducted using sensitive parameters related to soil loss from each HRU. To determine the magnitude of catchment sediment yield, the initial sensitive sediment parameters were calibrated using the global Sustainability 2022, 14, 5262 13 of 19 sensitivity analysis procedure. As a result, nine parameters with high to medium sensitivity (Table 6) were identified for sediment calibration process. As shown in Figure 6, the calibration of the SWAT model was successfully achieved using the observed daily sediment flow. Similarly, during the validation stage, the SWAT model adequately mimicked the sediment inflow of the river.  According to the outcome of calibration and validation, the model performance was examined according to statistical indicators (Table 7).  According to the outcome of calibration and validation, the model performance was examined according to statistical indicators (Table 7). As can be seen in Table 7, all numerical model performance measures were in an acceptable range, indicating that the SWAT model replicated the observed sediment yields. Similarly, the graphical representation of measured and simulated flows correlated well for both calibration and validation periods ( Figure 6). Hence, the result can be used to identify major sediment source areas within the sub-watersheds.

Spatial Distribution of Sediment Generation Hotspot Areas
The calibrated and validated SWAT model was used to simulate the effect of management/conservation measures on water and sediment yield in the watershed. The spatial variability of erosion rate was identified as shown in Figure 7.
To quantify the effects of catchment characteristics on the sediment yield rates of the sub-basin, the correlation of hotspot areas and catchment characteristics was determined.
As indicated in Table 8, all of the selected catchments in areas with high erosion risk (sub-basins 22, 2, 1, 19, and 3) were characterized by a single land-use class (agricultural) and one soil type (Chromic). Furthermore, the catchment prioritization study showed that soil loss is modest in low-lying locations with the same land use and soil type. This indicates that the variation of sediment yield is more sensitive to terrain slopes due to poor land-use practices in steep areas.
acceptable range, indicating that the SWAT model replicated the observed sediment yields. Similarly, the graphical representation of measured and simulated flows correlated well for both calibration and validation periods ( Figure 6). Hence, the result can be used to identify major sediment source areas within the sub-watersheds.

Spatial Distribution of Sediment Generation Hotspot Areas
The calibrated and validated SWAT model was used to simulate the effect of management/conservation measures on water and sediment yield in the watershed. The spatial variability of erosion rate was identified as shown in Figure 7.  On the basis of this finding, potential areas of intervention were prioritized, and four management scenarios were designed and simulated using the SWAT model in order to assess the basin's most appropriate management/conservation actions: scenario I (grassed waterway), scenario II (filter strip), scenario III (terracing), and scenario IV (contouring). The baseline scenario was used as a reference for comparisons of the effectiveness of the developed sediment reduction scenarios (Table 9).
In conclusion, after the application of terracing, the greatest reduction in sediment yield (62.32%) was observed. As the Gidabo irrigation dam is located at the outlet of the Measso gauging station (where the river basin was modeled), the application of the best selected scenario would decrease the sedimentation rate of the dam by 62.32%.  Table 9. Estimated sediment reduction due to conservation structures and best management practices as compared to the baseline scenario. The reductions in sediment yield for grassed waterway, filter strip, terracing, and contouring were 49%, 37.53%, 62.32%, and 54.6%, respectively.

Conclusions
The goal of this research was to use the SWAT model to analyze surface runoff generation and soil erosion rates in the Gidabo watershed. Both runoff and sediment fluxes were calibrated and validated using field-measured soil characteristics. To assess the degree of connection between measured and simulated monthly datasets, the sensitive flow and sediment parameters were considered in the process of calibration. During both calibration and validation periods, the model's performance was determined to be very good for both stream and sediment flow.
For each HRU, the sediment yield validated using the Soil and Water Assessment Tool model was shown to be correlated with land use, soil, and topography. Possible sediment basin regions for watershed prioritization and erosion regulation were distinguished on the basis of the sediment yield of the sub-basins. All of the sub-basins identified as sedimentation source locations shared a typical soil type and a typical land-use class, with varying terrain slopes. This shows that the difference in sediment yield is more sensitive to terrain slopes because of poor catchment management in steep areas. Moreover, the study findings indicated that the best choice to lessen the sediment flow of the watershed is to apply soil protection practices that decrease the catchment slope length. Accordingly, the application of terracing would result in the greatest reduction in sediment yield (62.32%).
In order to decrease the land degradation of the watershed and lengthen the beneficial lifespan of the dam located below the watershed, an organized strategy might be required.
As a result, the produced map ( Figure 7) indicating high-erosion areas can be used by decision-makers to undertake appropriate soil conservation measures in particular areas. In these erosion-prone areas, mandatory actions may also include strip planting, terracing, or contour farming to decrease the effect of the slope on surface runoff flow velocity and sediment transport capacity.