Understanding how watersheds respond to ecological changes and how LULC alteration affects watershed hydrology is crucial for water and soil resource management. LULC changes in the Ajora-Woybo watershed, Ethiopia, have momentously affected the water and soil resources. The researchers aimed to see how LULC changes affect hydrological components (HCs) and sediment yield (SED) in the watershed, both historically and in the future. The Soil and Water Assessment Tool (SWAT2012) and Partial Least Squares Regression (PLSR) models were used to investigate the contribution of each of the LULC classes to achieve the goal. The findings revealed a continual growth of cultivated land, built-up areas, and bare land, and a retreat of shrub land and forest land during the 2000 to 2020 periods, which is expected to continue in the 2035 and 2050 periods. Changes in LULC that happened over the historical era increased yearly surface runoff (23.5%), water yield (5.7%), and sediment yield (23.5%). On the other hand, the observed modifications have reduced lateral flow (12.8%) and groundwater flow (10.9%). Except for the 2020 LULC period, evapotranspiration decreased during the studied years. The future impacts of LULC states are predicted to increase in line with the historical trend. The PLSR results showed that cultivated land and built-up areas had a positive association with surface runoff and sediment yield, but shrub land and forest land had a negative correlation. This highlights the importance of controlling the LULC change as soon as possible to maintain long-term watershed stability.
Introduction
One of the critical concerns in global change studies has been the impact of natural and man-made variables on land cover change, particularly the transition from vegetation to another LULC (Assfaw, 2020). LULC refers to the bio-physical cover on the Earth’s surface and close sub-surfaces that affects water, sediment, and energy flows (Chilagane et al., 2021; Choto & Fetene, 2019). LULC changes are considered one of the most important elements influencing water and soil resources, which are the fundamental hydrological parameters for basin water resource planning and management (Leta et al., 2021). Because of their unavoidable consequences, human activities demand additional consideration when analyzing water and soil losses (Huang & Lo, 2015). The impact of LULC change on HCs such as surface runoff (SURQ), lateral flow (LATQ), groundwater flow (GWQ), and evapotranspiration (ET) in a river basin is a well-studied issue around the world (Abe et al., 2018; Chilagane et al., 2021; Kumar et al., 2018). LULC variations have been studied at various geographical and temporal ranges for their impact on hydrology and soil loss (Assfaw, 2020; Baig et al., 2022; Chaemiso et al., 2021; Choto & Fetene, 2019; Gebremichael, 2021; Wang et al., 2022). Furthermore, evaluations of multiple research papers indicate that LULC modification might modify streamflow, hence altering a basin’s hydrological components (HCs) and sediment yield (SED) (Afonso de Oliveira Serrão et al., 2022; Kumar et al., 2018; Negese, 2021).
Anthropogenic causes hastened LULC changes, which affect natural resources in general, and water and soil resources in particular (Negese, 2021). All of these causes contribute to a decline in living circumstances (Munoth & Goyal, 2020). Increased SURQ and decreased GWQ contribute to more extreme weather events, runoff, sedimentation, siltation, and water shortages during the dry season, all of which impede socioeconomic development (Leta et al., 2021). Droughts and floods, which have an impact on hydropower and irrigation production, may become more frequent and stay longer as a result of LULC change (Leta et al., 2021). Deprived LULC management and a lack of proper soil conservation methods are the main threats to watershed management and agricultural productivity (Pandey et al., 2021). As a result, understanding how LULC changes affect water and soil resources is critical for watershed management (Nie et al., 2011).
Although the hydrologic cycle at the watershed level is a complex process impacted by geology, geography, soil conditions, and vegetation, it is also influenced by human actions (e.g. land-use practices, reservoir restrictions, and water transfer projects) (Samie et al., 2019). Also, the hydrology of local watersheds can vary greatly, and water flow patterns are frequently determined by a unique combination of soil, LULC, and elevation variables (Tadesse et al., 2015). River basins are essential for hydrological and environmental improvements because they create consistent streamflow from baseflow and runoff (Teshome et al., 2020). As a result, knowing the impact of LULC change on water and soil resources is critical, particularly in nations whose civilizations rely heavily on rain-fed agriculture (Daramola et al., 2022).
According to Liu, Schmalz et al. (2022), understanding historical LULC changes and their interactions with hydrological processes is critical for better supporting future land-use strategies. Conversely, most scholars do not quantify the future repercussions and contributions of certain LULC classes to diverse hydrological responses (Nie et al., 2011). Understanding in what way different LULC types affect watershed hydrology will extensively advance the liability of LULC dynamics on HCs for water resource management (Woldesenbet et al., 2017). Because external factors that influenced the previous shift in LULC change continue to have an impact on the future direction of LULC change (Schilling et al., 2008). For that reason, a better knowledge of LULC change and watershed hydrology at the watershed scale can help to prevent negative consequences.
In Ethiopia, particularly in the Ajora-Woybo watershed, the impact of the LULC change study on hydrological components, soil erosion, and sedimentation is still not addressed, mainly on the scales that are most important for local actors. According to Stefanidis and Stathis (2018), developing cost-effective watershed management and identifying critical soil erosion-prone areas are necessary for the implementation of appropriate mitigation measures. Soil and water conservation have thus been implemented in the watershed for approximately a decade. Nevertheless, its implementation has been hindered by the absence of site-specific, scientifically calculated soil erosion data and priority criteria. It is also unclear how changes in each LULC class affect the HCs and SED of the watershed. The response to this question will aid in estimating the hydrological implications of each LULC modification, which will be useful in future water resource planning and management.
Diverse models are now utilized to estimate the performance of watersheds under LULC changes on hydrological processes and sediment yield with modest modifications and direct application. The combination of LULC change and HCs modeling provides a valuable framework for assessing the environmental impacts of a watershed LULC change (Shao et al., 2013). Hydrological modeling (HM) is a good way to look into the impact of land use, climate, and soil interactions on natural ecosystems (Scott-Shaw et al., 2020). An HM approach was also employed to evaluate the potential implications of past as well as future LULC on HCs and SED (Abe et al., 2018). SWAT is a popular model among hydrologists because it places a heavy emphasis on the environment and hydrological interactions, making it appropriate for LULC-based water and soil losses (Aragaw et al., 2021; Assfaw, 2020; Kumar et al., 2018). SWAT is a semi-distributed and continuous time-step model based on GIS that allows for speedier implementation of input information across diverse river basins and gives an efficient and cost-effective understanding of the effects of LULC change on water and soil resources (Arnold et al., 2012; Samal & Gedam, 2021). Some of the studies that demonstrated the applicability of SWAT in the study basin (Chaemiso et al., 2021; Choto & Fetene, 2019). Therefore, in this study, an integrated technique of SWAT and the PLSR model was applied to evaluate the situation.
In general, this research assessed the HCs and SED in the Ajora-Woybo watershed, utilizing historical and future LULC change scenarios. Understanding the effects of LULC change is necessary for long-term water resource sustainability, adaptation options, and effective water resource management in the future. Additionally, the study evaluated the localization of erosion-prone sites. Accurate quantitative estimates of soil loss rates are crucial for developing and putting into practice appropriate soil and water conservation policies. This research could also aid policymakers in making more long-term judgments when executing watershed development plans in a particular area.
Material and Methods
Description of the study area
The Omo-Gibe basin in Ethiopia is 79,000 km2 in area, measuring 550 km in length and 140 km in width. The Ajora-Woybo watershed, located on the eastern edge of the middle Omo-Gibe basin, is one of the basin’s principal watersheds. It was chosen as the case study for this investigation. The Ajora-Woybo watershed contains four perennial rivers (Ajancho, Soke, Shapa, and Woybo) that flow downstream to the lower valley and finally form two rivers, Ajora and Woybo, which eventually join the Gibe III dam reservoir (Figure 1). In addition, the rivers are fed by several tiny intermittent tributaries. The watershed is located between latitudes 7°2′0″N and 7°50′0″N and longitudes 37°30′’0″E and 38°0′’0″E. The climate of the Ajora-Woybo watershed ranges from a hot arid environment in the floodplain to a tropical humid climate in the highlands, which encompass the extreme north and northeastern parts of the watershed. The climate is tropical and sub-humid. Weather patterns and agricultural production practices vary by zone. Annual precipitation varies substantially from north to south, ranging from approximately 1,900 mm in the north to less than 300 mm in the south. Furthermore, the watershed rainfall regime is bimodal. The mean annual temperature ranges from 16°C in the northern highlands to more than 29°C in the southern lowlands.
Data
The SWAT model requires the following input data, as stated in Table 1.
Table 1.
Data for the Model.
DEM: It was entered into the SWAT model to determine the watershed, extract topography or slope data, and evaluate terrain drainage patterns. A 20 m × 20 m DEM was used for this study.
LULC: The historical LULC images were obtained from Landsat images (4–5 TM, 7 ETM+ and 8 OLI) and classified using supervised classification in Earth Resource Data Analysis System (ERDAS) imagine 2015 model. Based on classified historical satellite photographs, the future LULC was anticipated. The Land Change Modeler included in TerrSet model version 18.31 was used for this. As a result, the classified historical LULC data for three time periods (2000, 2010, and 2020) as well as the two forecasted time periods (2035 and 2050) were displayed in detail by Toma et al. (2022) and utilized in this study (Figure 2). The classified LULC classes of the Ajora-Woybo watershed are cultivated land (CL), built-up areas (BUA), shrub land (SHL), forest land (FL), bare land (BL), and water body (WB) (Toma et al., 2022).
Soil: Physicochemical and hydrological parameters, were determined for each soil polygon. Dystric nitisols (45.8%), Pellic vertisols (38.5%), Chromic luvisols (10.5%), Eutric cambisols (4.6%) and Leptosols (0.4%) are the major soils in the watershed (Figure 3).
Description and set up of SWAT model
The model was developed over a 30-year period by the US Section of Agricultural Research Service. It is used to simulate the impact of LULC difference on HCs and SED over lengthy time periods with changing LULC, climate, soils, and management (Arnold et al., 2012). For this investigation, a SWAT2012 interface compatible with ArcGIS version 10.3 was used.
The delineation of the watershed is the first stage in creating SWAT model input. Prior to dealing with raster input data such as the DEM, LULC, and soil maps, which were all projected into Ethiopia’s projection parameters (UTM Zone 37N). All of these physical parameters were overlaid for hydrological response units (HRUs) definition after reclassifying LULC, soil, and slope in the SWAT database. The lowest area threshold values for land use, soil, and slope were defined as 5%, 5%, and 5%, respectively, and 25 sub-basins and 192 HRUs were identified, representing inimitable amalgamations of land use, soil type, and slope.
The simulation process is allocated into two parts: (1) the land phase and (2) the routing phase, which are used to simulate the HCs and SED of the Ajora-Woybo watershed. The first phase takes into account HCs such as SURQ, LATQ, GWQ, water yield (WTYD), and ET. The watershed channel network is used to track the flow of water and SED in the second phase. The simulation of the model is established on the water balance equation (Arnold et al., 2012; Neitsch et al., 2005).
Where, SWT is the last soil water satisfied, SWO is the early water satisfied, n is the time, Rday is the amount of rainfall, Qsurf is the volume of SURQ, Ea is the extent of ET, Wseep is the extent of water arriving the vadose zone, and Qgw is the degree of return flow. Except for n, which is measured in days, other physiognomies are measured in millimeters.
The Soil Conservation Service curve number (SCS-CN) approach is used to estimate SURQ (Neitsch et al., 2005). The SCS-CN equation is:
Where, Qsurf is the total SURQ, Rday is the daily rainfall depth, S is the retention. The equation that defines the S limit was:
Groundwater flow: It is calculated using the hydraulic conductivity, the aloofness from the sub basin to the main channel, and the water table height.
Where, aqsh,i is the water kept in the aquifer on day i, aqsh,i−1 is the water kept in the aquifer on day i−1, Wrchrg is the recharge entering the aquifer, Qgw is the GWQ, Wdeep is the deep percolating water, and Wpump, sh is the water extracted from the aquifer. All parameters are in mm.
ET: Evaporation and transpiration are calculated separately: evaporation is assessed by means of a linear function of potential evapotranspiration and leaf area index, while transpiration is assessed by means of an exponential function of soil depth and water content. In this study, the Penman–Monteith technique was used to compute ET.
WTYD: The residual water yield to the stream channel was calculated by adding SURQ, LATQ, GWQ, and transmission loss (Tloss).
SED: The SED for individual sub-basin is computed, and the total SED is directed to the watershed outlet. The Modified Universal Soil Loss Equation (MUSLE) was applied to calculate SED at an HRU level (MUSLE). Since SURQ is a result of precursor moisture conditions as well as rainfall drive, this strategy enhanced sediment output prediction (Chimdessa et al., 2018; Neitsch et al., 2005).
Where, Sed is the SED (metric tons/day), Qsurf is the SURQ (mm/ha), qpeak is the peak runoff (m3/s), areahru is the area (ha), K is the soil erodibility factor, C is the cover factor (dimensionless), P is the USLE conservation practice factor (dimensionless), LS is the topographic factor (dimensionless) and CFRG is the course fragment factor (dimensionless).
Weather Generator (WGEN): All-climate variables must have daily values derived from measured data, as well as values derived from monthly average data over several years. Nonetheless, there is a scarcity of complete and representative long-term climatic data in the study area. As a result, the weather generator answers this issue by producing data from observed data (Sodo station). The weather generator produces precipitation (as PCP), maximum and minimum temperatures (as TMP), solar radiation (as SLR), relative humidity (as HMD) and wind speed (as WND) individually.
Sensitivity analysis
During model simulation, a model sensitive analysis is required to find the most sensitive parameters and decrease parameter joblessness (Pandey et al., 2021). The sensitivity of each parameter was determined using a global sensitivity analysis. The most sensitive parameters that greatly influence the flow and sedimentation processes were found individually using the sensitivity analysis. The parameters were ranked using the most recent iterations of t-stat and p-value. As a result, a higher absolute value t-test implies greater sensitivity, while a p-value of zero shows greater significance. Based on the literature and the weight of parameters, sixteen and nine parameters for flow and sediment, respectively, were examined and tested for sensitivity analysis (Tables 2 and 3).
Table 2.
Sensitive Parameters for Flow.
Table 3.
Sensitive Parameters for Sediment Yield.
Calibration and validation
Monthly flow (1990–2015) and suspended sediment load (1990–2015) at Woybo and Soke River gaging stations were utilized independently to calibrate and validate the model. To do this, simulations from 1990 to 1992 were utilized for warm-up, 1993 to 1997, 2001 to 2006, and 2011 to 2013 for calibration, and 1998 to 2000, 2007 to 2010, and 2014 to 2015 for validation. The method in this algorithm attempts to capture the majority of the measured data inside the 95% prediction uncertainty (95PPU) band. Validation was performed after calibration with no further changes in the values of the sensitive calibration parameters. The Sequential Uncertainty Fitting (SUFI-2) approach, which typifies uncertainty, was used in SWAT-CUP2019.
Model competence
It was evaluated by comparing the model simulation results in relation to the observed data (Moriasi et al., 2007). Most HM are evaluated by comparing simulated and measured hydrographs using coefficients of determination (R2), Nash-Sutcliffe efficiency (NSE) and percent bias (PBIAS) (Table 4).
Table 4.
Model Efficiency Scores by Moriasi et al. (2007).
Where, Xi = measured flow (m3/s), Yi = simulated flow (m3/s), Xav = average measured flow (m3/s) and Yav is average simulated flow (m3/s).
Sediment rating curve (SRC)
It is a frequently used method for assessing the suspended sediment load (SSL) carried by a river by establishing a link between river discharge and sediment load (SL) (Assfaw, 2020). Because SSL records lack a continuous time step, the SRC for this investigation was generated from empirical relationships between river discharges and related SSL (Arnold et al., 2012). The overall correlation of the SRC includes SL (ton/day), discharge (cumecs) and constants (a and b).
The initial task was to convert the measured SSL (mg/l) data at Woybo and Soke stations into SL (ton/day) using:
Where, S is the SL (ton/day), Q is the discharge (m3/s), C is the SSL (mg/l), and 0.0864 is the conversion factor. Following the calculation of the SL, the next step was to establish a relationship between the continuous observed discharge and SL in a daily time step (Assfaw, 2020) (Figure 4).
Scenario development
After determining the most sensitive values for the specified parameters, the SURQ and SED were simulated by altering the LULC map under constant DEM, soil, and climate conditions. LULC maps from the years 2000, 2010, and 2020 were utilized for historical scenarios, while LULC maps from the years 2035 and 2050 were used for future scenarios.
Partial least square regression
The PLSR approach is a multivariate statistical technique that combines multiple regressions and principal component analysis (Abdi, 2010). It uses a set of predictors to foresee a set of dependent variables. The PLSR model has recently been employed for a variety of hydrological applications, including investigating the effects of changes in specific LULC classes on watershed hydrology at the sub-basin level at a 95% significant level (Gashaw et al., 2018; Nie et al., 2011; Twisa et al., 2020; Woldesenbet et al., 2017). As a result, the PLSR model is beneficial in determining the impact of particular factors (Gashaw et al., 2018; Woldesenbet et al., 2017). The model explains the correlation between a response and a set of predictor variables in such a way that:
Where, a and x are the regression coefficients generated from the data, that is, the intercept and variables (1 through n). As a result, it is easy to figure out which LULC types interact the most with the HCs and SED. Annual HCs and SED were employed as dependent variables in PLSR in this study, while individual LULC classes were used as independent variables. The variable importance of projection (VIP) was also utilized to determine the degree of correlation between independent and dependent variables (Aragaw et al., 2021). The importance of the LULC classes for the dependent variables is better explained by predictors with higher VIP values. In general, the scientific community considers a VIP value of 0.8 to be the minimum acceptable value (Abdi, 2010; Gashaw et al., 2018). The changes in dependent and independent variables were considered using LULC maps between 2020 and 2000 for this study.
Identifying places that are prone to erosion
The simulated result of annual average SED mapping was performed after utilizing the calibrated and validated SWAT model. Following that, significant sub-basins prone to soil erosion were identified and selected based on the simulated results for the current (2020LULC) timeframe. The Ajora-Woybo watershed was divided into low, moderate, high, very high, severe, and very severe risk areas, with 0–10 t ha−1 y−1, 10–20 t ha−1 y−1, 20–30 t ha−1 y−1, 30–50 t ha−1 y−1, 50–75 t ha−1 y−1, and >75 t ha−1 y−1, respectively, to better understand the prioritization and classification of SED severity. The soil loss ranks were established based on Ethiopian study findings (Assfaw, 2020; Chaemiso et al., 2021; Chimdessa et al., 2018; Dibaba et al., 2021).
Results and Discussions
Land use land cover change
Figure 5 depicts the LULC stages of the Ajora-Woybo watershed between the years 2000 and 2050. The most substantial changes occurred in six LULC classes, according to a comparison of LULC maps for the years 2000, 2010, 2020, 2035, and 2050. The assessment of LULC change shows that SHL and FL have been decreasing. However, high increasing rate of CL, BUA, and BL. Similarly, an increase in WB occurred from 2020 to 2050. It is clear from the results that modifications occurred in all LULC classes (Toma et al., 2022). Details on the historical and prospective LULC change analyses, as well as the key causes of LULC change in the Ajora-Woybo watershed, were conducted according to (Toma et al., 2022).
Sensitivity analysis
The greatest influential sensitive parameters for flow and SED were discovered and prioritized according to their sensitivity using the t-stat and the p-value at Woybo and Soke gaging stations in this study (Tables 5 and 6). The sensitivity analysis revealed that the parameters played a significant role in affecting stream flow and SED in the Ajora-Woybo watershed. These variables were also mentioned as being important in the Omo-Gibe basin (Chaemiso et al., 2021; Choto & Fetene, 2019), though the rank and fitted value differ. As a consequence, the first 13 flow parameters and the first six SED parameters were used to calibrate and validate the model.
Table 5.
Flow Sensitive Parameters.
Table 6.
Sediment Sensitive Parameters.
Calibration and validation
For both the calibration and validation periods, monthly scales for the Woybo and Soke gage stations were verified (Table 7 and Figures 6 and 7). Based on the recommendation standards delivered by (Moriasi et al., 2007), the overall SWAT model performance for the research region was judged to be good to very good throughout both the calibration and validation periods. The R2, ENS, and PBIAS values for monthly flow rate and sediment yield imply the model predictive capacity, and it is appropriate to estimate the HCs and SED owing to LULC changes in the Ajora-Woybo watershed. Similar model efficiency perfections were also recorded during the calibration and validation periods by Chaemiso et al. (2021) and Choto and Fetene (2019). Following that, the information can be used to identify and designate sub-basin HCs and important sediment source locations.
Table 7.
SWAT Model Performance Evaluation.
Response of HCs and SED to LULC change
Based on the simulation results, HCs and SED responses to LULC changes were examined on a seasonal basis, namely the main rainy season, also known as kiremt (June–September), the small rainy season, also known as Belg (February–May), and the dry season, also known as Bega (October–January).
Seasonal analysis shows (Figures 8 and 9), during the main rainy season (June-September), SURQ for the 2000, 2010, 2020, 2035 and 2050 LULC was 142.2, 171.2, 204.2, 208.9, and 218.9 mm, while the WTYD was 264.2, 275.8, 288.3, 295.3, and 306.3 mm, respectively. During the same season, GWQ was 142.9, 138.9, 133.7, 130.1, and 124.1 mm, while the ET was 91.4, 95.2, 96.2, 96.8, and 97.7 mm. Additionally, in the same season, SED was 17.2, 19.4, 23.5, 25.9, and 28.2 t/ha, respectively. During the small rainy season (February-May), SURQ for the 2000, 2010, 2020, 2035 and 2050 LULC was 104.5, 116.6, 128.4, 130.7, and 131.6 mm, while the WTYD was 132.1, 137.7, 142.4, 144.5, and 145.8 mm, respectively. During the same season, GWQ was 94.8, 78.3, 73.1, 69.5, and 63.4 mm, while the ET was 84.9, 86.4, 88.5, 89.1, and 89.9 mm. Additionally, in the same season, SED was 9.3, 10.1, 11.7, 12.5, and 13.5 t/ha. During the dry season (October-January), SURQ for the 2000, 2010, 2020, 2035 and 2050 LULC was 105.4, 103.4, 102, 100.8, and 98.1 mm, while the WTYD was 156, 155.2, 151.8, 147.8, and 142.8 mm, respectively. During the same season, GWQ was 86.1, 85.1, 81.7, 78.8, and 75 mm, while the ET was 67.2, 65.8, 64.1, 63.3, 62.8, and 61.2 mm. Additionally, in the same season, SED was 6.3, 5.8, 5.1, 4.3, and 3.75 t/ha.
According to the results, the kiremt season accounts for more than 60% and 45% of SURQ and GWQ for all LULC time periods, respectively, whereas the belg and bega rainy seasons account for less than 40% of SURQ (Figure 8). In addition, the SED rate is the highest in the kiremt season (60%), followed by the belg season (29%) and the bega season (11%). The seasonal results show that due to LULC changes, the SURQ in the watershed increased during the kiremt and belg season and declined during the bega. Due to the increase in agricultural land and other changes in land use during the rainy season, the SURQ of the Ajora-Woybo watershed rises more during the rainy season than during the dry. Similarly, the groundwater component decreased overall across all seasons. This suggests that the LULC alterations have reduced groundwater yield in the watershed. Similarly, changes in lateral flow and evapotranspiration due to LULC variations are inferred on a monthly and seasonal basis (Figure 8). The seasonal results show that ET in the watershed increased during the kiremt and belg seasons and reduced during the bega season. This could be due to soil moisture and the condition of the canopy cover during the rainy season.
Heavy rainfall during kiremt season causes large seasonal floods and SURQ over Ethiopia. The reason for this could be that SURQ was more vulnerable in the kiremt than in other seasons. The shift in HCs caused by LULC change lowered runoff during the dry season, which is largely supplied by base flow, and increased runoff during the rainy season, which is mostly supplied by SURQ (Leta et al., 2021; Zhang et al., 2016). According to Woldesenbet et al. (2017), during the dry season, decreased actual ET and GWQ have an opposing effect to streamflow availability. Furthermore, Benegas et al. (2014) proved that improved infiltration is responsible for a decrease in SURQ. High surface runoff caused by LULC changes, as observed in this study, combined with heavy rainfall could result in significant SED during the rainy season. This was also mentioned by Woldesenbet et al. (2017), who stated that increased SURQ during the rainy season has the potential to produce flooding and extreme soil loss. Changes in LULC also have an effect, particularly during the dry season, which has a negative influence on biotic components of ecosystems both inside and beyond the watershed, as well as hydropower generation downstream of the watershed.
The annual variations in the HCs and SED caused by LULC changes revealed effects on the components. Table 8 summarizes the findings, showing that the HCs and SED in the watershed were altered in diverse ways. SURQ, WTYD and SED all increased from the first scenario (2000LULC) to the fifth scenario (2050LULC), whereas LATQ and GWQ dropped. Except for the 2020 LULC period, ET decreased during the studied years. This could be related to the occurrence of a significant body of water from the Gibe III reservoir backflow in the outlet area of the Ajora-Woybo watershed.
Table 8.
Annual Average HCs and SED in the Ajora-Woybo Watershed.
Furthermore, SURQ increased by 25.5%, WTYD climbed by 5.7%, and SED increased by 23.5% between 2000 and 2020, whereas LATQ and GWQ declined by 12.8% and 10.9%, respectively. SURQ (by 3.1%), WTYD (by 2.3%), and SED (by 12.8%) increased during the research period (2020–2050), while LATQ (by 14.6%), GWQ (by 8.9%), and ET decreased (by 2.9%). The comparison of SURQ increment and changes in LULCs reveals that the rise in average annual SURQ might be attributable mostly to the growth of cultivated land, built-up areas, and bare land. Besides, the study found that the decrease in LATQ, GWQ and ET are associated with changes in shrub land and forest cover. This was also noted by Nie et al. (2011). They contend that loss of canopy cover reduces percolation and base flow.
For HCs and SED indicators, individual values in 25 sub-basins showed a decreasing and growing pattern between 2000 and 2020. In contrast to HCs, sub-basins with greater growth in agricultural and built-up regions generate more runoff and net water yield. The influence of LULC change on watershed hydrology and sediment output was greater at the sub-basin level than at the watershed level (Figure 10). According to Aragaw et al. (2021), depending on the extent of changes in LULC classes, the influence of LULC class alterations on HCs varies among sub-basins. Choto and Fetene (2019), also observed that sub-basins with a higher share of agriculturally used areas have higher runoff and sediment outputs. As a result, fluctuations in LULC can have a variety of effects on the river flow and SED generating features of any sub-basin. This shows at a smaller scale, the special effects of LULC variations were more pronounced.
The conclusions of this study are reliable with other investigations (Aragaw et al., 2021; Assfaw, 2020; Chaemiso et al., 2021; Chimdessa et al., 2018; Choto & Fetene, 2019; Woldesenbet et al., 2017). They determined that the amount of cultivated land, built-up areas, and bare land is increasing at the expense of vegetation cover (shrub land and forest land), resulting in increased surface runoff and sediment accumulation, as well as decreased GWQ. Many recent articles from around the world have discussed how LULC change influences stream flow and sediment output, as well as the effects of LULC change on water and soil resources (Abe et al., 2018; Chilagane et al., 2021; Daramola et al., 2022).
Individual LULC change effects on HCs and SED
Table 9 shows the regression correlations of the six LULC class with HCs and SED. According to the findings, almost every LULC type had a substantial link with the specific HCs and SED. For instance, cultivated land, built-up areas, and bare land showed positive (significant) correlations with SURQ (0.98), WTYD (0.97), ET (0.91) and SED (0.98), and negative correlations with LATQ (−0.98) and GWQ (−0.94). According to the findings, cultivated land and bare land are the largest contributors to both SURQ and SED, followed by built-up regions, which contribute a lot to SURQ but not much to SED. Shrub land and forest land were shown to contribute to SURQ and SED in opposite directions, with some variance in their contributions to SURQ. Table 10 summarizes the four PLSR models that were built separately.
Table 9.
PLSR Correlation for LULC Changes (2000–2020) in HCs and SED.
Table 10.
VIP and Weights (W*) of HCs and SED in the Ajora-Woybo Watershed.
To differentiate the kin status of predictors, cultivated land, shrub land and forest land (VIP > 1.0) has additional significance in this model. Changes in SURQ and WTYD were shown to be dependent on cultivated land, built-up areas, bare land, and water bodies (on the positive regression coefficient), as well as shrub and forest land (on the negative regression coefficient). The highest VIP values for SURQ and WTYD are obtained for shrub land (VIP = 1.046), followed by cultivated land (VIP = 1.034), and forest land (VIP = 1.024). VIP values for LATQ and GWQ are highest for shrub land (VIP = 1.038), followed by forest land (VIP = 1.036), cultivated land (VIP = 1.023), and bare land (VIP = 1.004). The water body has the greatest VIP value for ET (VIP = 1.306), followed by cultivated land (VIP = 1.206) and shrub land (VIP = 1.172). Cultivated land (VIP = 1.092) has the greatest VIP rating for SED, followed by shrub land (VIP = 1.091), and water body (VIP = 1.051). The weight value in Table 10 demonstrated the relatively low and high impact of LULC classes in influencing HCs and SED of the Ajora-Woybo watershed.
Considering separable LULC classes for each HCs and SED, conversion in cultivation land exaggerated SURQ and WTYD positively, but prejudiced LATQ and GWQ negatively (Tables 9 and 10). This suggests that as the Ajora-Woybo watershed’s cultivated area expanded, SURQ climbed while LATQ and GWQ decreased. The increased SURQ, coupled with a decrease in actual ET due to the expansion of cultivation land and a decline in shrub land and forest land, would result in increased streamflow mainly in the rainy season and decreased in the dry season. Deforestation in the watershed reduced actual ET while increasing surface run-off SURQ and WTYD. Similar to the findings of this study, an escalation in surface runoff, as well as a decrease in groundwater flow, and ET, was seen in the Gidabo basin by Aragaw et al. (2021) and the Andassa watershed by Gashaw et al. (2018), in two watersheds of central America by Benegas et al. (2014), in the Wami basins of Tanzania by Twisa et al. (2020) and the Xunwu watershed of China by Liu, Schmalz et al. (2022).
Identification of hotspot areas in SED
The sub-basin was classified in to different categories of soil erosion categorization based on study experiences in Ethiopia (Assfaw, 2020; Chimdessa et al., 2018; Choto & Fetene, 2019; Megersa et al., 2019). It was assumed to allocate the priority stages of I–VI (Table 11). For present scenarios (2020 LULC), the largest possible SED in each sub-basin was analyzed and plotted, as shown in Figure 11. Sub basins 1 and 2 are falling under very severe erosion classes, while sub-basin 3, 7, 8, 12, 16, and 17 are falling under severe erosion classes. Furthermore, sub-basins 5, 9, 11, 14, 19, 21, and 22 fall under very high erosion classes, while the remaining sub-basins fall under high, moderate and low erosion classes. This graphic depicts how much the watershed is currently vulnerable to erosion as a result of more intensive agricultural activity and urbanization. It is necessary to implement proper best management scenarios in places that are prone to erosion.
Table 11.
Sediment yield (t/ha) and treatment priority order based on the 2020 LULC.
Conclusions
This study predicted the potential impacts of LULC change on annual and seasonal HCs and SED distribution, in the Ajora-Woybo watershed of Omo-Gibe River basin, Ethiopia. According to the model efficiency norms used in this work, the SWAT model was able to predict the HCs and SED processes in the Ajora-Woybo watershed adequately. As a consequence, it was able to simulate changes in watershed HCs and SED in response to various LULC situations. The PLSR model was utilized to determine the relationship between LULC type, HCs, and SED.
The CN2, ESCO, ALPHA_BF, GW_SPYLD, and SOL_BD model parameters have the greatest sensitivity on streamflow, whereas the SPEXP, USLE C, USLE P, and CH_EQN model parameters have the greatest sensitivity on SED. The model runs for different LULC periods show that land use change can induce significant HCs variance and an increase in SED. This is primarily due to land degradation (conversion of shrub land and forest land to cultivated land, built-up areas and bare land). Thus, cultivated land development, urbanization, and exposed land are the most significant contributors to SURQ, WTYD, and SED, and can be regarded as the main environmental stressor affecting HCs and SED in the Ajora-Woybo watershed. A decrease in LATQ and GWQ would immediately reduce recharge to shallow aquifers, and would thus be considered a negative influence on watersheds. Individual changes in cultivated land, built-up areas, shrub land, and forested regions were the most important factors of SURQ and SED at any particular time. Changes in cultivated land, built-up areas, and bare land were shown to correlate favorably with changes in SURQ and SED. However, shrub land and forest land were found to correlate negatively with changes in SURQ and SED. The PLSR model was beneficial in examining the effects of LULC classes on HCs and SED.
These findings are predicted to be relevant in the future for sustainable watershed management in the Ajora-Woybo watershed in particular and the Omo-Gibe basin in general. This research also shows how SWAT-based watershed models can be used to analyze LULC change in a watershed and how they can be simply implemented in similar settings around the world. Still, the outcomes of this study could assist environmental managers in developing long-term watershed based LULC adjustments.
Acknowledgements
We would like to thank, Ministry of Water and Energy, and National Meteorological Agency of Ethiopia, for providing hydrological and meteorological data for the study.
Author Contributions
Conceptualization, M.B.T., M.D.B. and M.D.U.; methodology, M.B.T., M.D.B. and M.D.U; software, M.B.T.; validation, M.B.T., M.D.B. and M.D.U; formal analysis, M.B.T., M.D.B and M.D.U; investigation, M.B.T.; writing—original draft preparation, M.B.T.; writing—review and editing, M.B.T., M.D.B. and M.D.U.; supervision, M.D.B. and M.D.U. All authors have read and agreed to the published version of the manuscript.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data Availability Statement
Data are available from the corresponding authors upon reasonable request.