Monitoring soil quality of different land use systems: a case study in Suha watershed, northwestern highlands of Ethiopia

The problem of soil quality degradation has been becoming more severe in the highlands of Ethiopia due to soil erosion; land use and land cover change, and poor land management. The level of soil quality degradation was not well known and documented in the study area and the results of this study could provide new information to improve soil conditions. The present study was conducted to evaluate soil quality in terms of its physical and chemical fertility under different land use types in the Suha watershed, northwestern highlands of Ethiopia. A total of 27 composite surface soil samples (0–30 cm) were collected from adjacently located land-uses in three replications from two elevation gradients. Standard procedures were followed to analyze selected soil physical and chemical quality indicators. The differences in the mean values of the parameters were tested using a two-way analysis of variance. In addition, Soil Quality Degradation Index was evaluated to see the direction and magnitude of change in soil quality indicators. The analysis of variance results revealed that soil quality indicators such as index of soil aggregate stability, organic carbon (OC), total nitrogen (TN), and C:N ratio were significantly decreased in the cultivated land use sys-tem compared to other land use systems. On the other hand, the content of available Phosphorus was significantly higher in the cultivated land. Soil quality deterioration index values were highly negative for SOC (− 71.3%) and TN (− 67.7%) in the cultivated land, followed by grazing land (SOM = − 35.5% and TN = − 27.7%). Aggregated Soil Quality Index values also indicated that the status of soil quality under cultivated fields is rated as low, grazing land as optimal, and forest land as high. Generally, results indicated that land use and cover changes had adverse effects on soil quality indicators. Hence, soil management strategies, mainly Integrated Soil Fertility Management which integrates soil and water conservation strategies, are required to alleviate the problem of soil quality deterioration and improve agricultural productivity.


Introduction
Soil degradation, manifested by the deterioration of soil's physical, chemical, and biological properties, is a global environmental problem that affects agricultural productivity, food security, climate change, and biodiversity (Young et al. 2015;Dagnachew et al. 2019).According to Nkonya et al. (2016), the negative consequences of soil degradation have resulted in a 60% reduction in the ecosystem services provided by soil, which now cost the globe between 6.3 and 10.6 trillion dollars annually.In Sub Sahara Africa (SSA), soil degradation has threatened the livelihood of the population since 67% of the land resource is degraded in different severity classes (Sileshi et al. 2019).In these countries, the major soil constraints include nutrient depletion, soil erosion, Al toxicity, low soil depth, high P-fixation, low CEC, salinity, and vertic properties (Zingore et al. 2015).Soil degradation is more severe in Ethiopia, and nutrient mining is the highest compared to SSA (Teferi et al. 2016).Soil degradation in the form of soil erosion and loss of soil fertility in the highlands of Ethiopia causes a reduction in agricultural productivity (Gashaw et al. 2014;Meseret 2016), which in turn impacts the livelihoods of the rural community and the economy of the country (Kassa et al. 2013;Adugna and Abegaz 2015).The study of Zingore et al. (2015) revealed that negative nutrient balances of macronutrients (N = 40-70 kg/ha, P = 7-10 kg/ha, and K = 33-50 kg/ha) were the highest in Ethiopia compared to sub-Saharan Africa (SSA).
Soil quality can be defined as the capacity of soil to provide ecosystem services in a natural or managed environment (Andrews et al. 2004).This parameter cannot be measured directly; rather, it is determined based on soil quality indicators (physical, chemical, and biological properties) (Zornoza et al. 2015).Land use and land cover changes are the major drivers of soil quality deterioration and other environmental changes like soil erosion, climate change, loss of biodiversity, the carbon cycle, and impairment of ecosystem services (Wang et al. 2016;Singh and Benbi 2018;Lepcha and Devi 2020).In Africa, soil quality deterioration has been aggravated due to population pressure and subsistence farming associated with poor land use systems, which in turn has resulted in declining agricultural productivity and food insecurity (Obalum et al. 2012;Alarima et al. 2020).Previous studies conducted in Ethiopia and elsewhere indicated that land use and land cover changes have detrimental effects on soil quality indicators (Tesfahunegn 2016;Delelegn et al. 2017;Pham et al. 2018).Changes from natural forest to agricultural land are the major causes of extensive soil erosion, soil quality deterioration, depletion of soil nutrients, and declining agricultural production in the highlands of Ethiopia (Aredehey et al. 2019;Molla et al. 2022).Dagnachew et al. (2019) and Teferi et al. (2016) also reported that key soil quality indicators are significantly affected due to human-induced activities in the Upper Blue Nile basin.The highlands of Ethiopia faced critical soil-related problems from the perspective of soil quality and fertility.These include extensive soil erosion, soil acidity that covers about 40% of cultivated fields, a high rate of soil organic matter loss, and depletion of macronutrients (high rate of negative nutrient balances).Addressing the problems of soil degradation, particularly the soil fertility component, has become a major policy issue in Ethiopia since the agricultural sector is the major pillar of economic and social development (Neglo et al. 2021;Haregeweyn et al. 2015).
The Suha watershed has a high potential for crop production in the northwestern highlands of Ethiopia.However, the natural resource base, particularly the soil resource, has been under stress due to high population pressure, and the soil condition has been deteriorating over time.The problems of soil fertility remain a major challenge and have not yet been reversed in this watershed (Simane et al. 2013).Maintaining soil quality is a key factor for agricultural productivity and environmental sustainability (Elias 2004).However, soil quality deterioration and nutrient depletion are major problems observed in the north western highlands of Ethiopia including the study watershed.Moreover, the alteration of soil quality indicators is site-specific due to the difference in soil management strategies, agroecologies, soil types and landscape positions.Therefore, regular monitoring of soil quality is vital to know whether this parameter is maintained or degraded and to apply appropriate and site-specific management strategies to improve soil conditions.However, the level of soil quality degradation was not well known and documented in the study area and the results of this study could provide new information to improve soil conditions.Hence, the results of this study could provide information to land managers and experts on how to select sets of site-specific soil management strategies to increase agricultural productivity.Therefore, the general objective of this study was to evaluate the status of soil quality in response to land use and cover change in terms of soil physical and chemical fertility along the topo-sequence of the Suha watershed, in the northwestern highlands of Ethiopia.

Study area description
The study area, the Suha watershed, is located in the north-western highlands of Ethiopia and geographically lies from 37° 56′ 15″ to 38° 18′ 49″ E and from 10° 06′ 46″ to 10° 41′ 56″ N (Fig. 1).Its elevation ranges from 1040 to 3986 m above sea level and covers 80,343 ha.This watershed has seven topographic categories (from flat to very steep) based on FAO's (2006) classification.

Climate
The climate is the most important factor that determines the type of land use system and has a greater impact on agricultural productivity.The study area has a unimodal type of rainfall that extends from June to September, and based on the data obtained from nearby stations (Bichena, Dejen, Kuy, and Robgebeya), it receives annual rainfall that ranges from 1213 mm at the lower part to 1396 mm in the upper part of the watershed.The mean minimum and mean maximum temperatures are 8.3 °C and 23.6 °C respectively.

Site selection and soil sample collection
Before soil sample collection, the watershed was stratified into three strata based on its altitudinal gradient using the digital elevation model (DEM) of the watershed.These strata were classified as the upper (2700-3900 m.a.s.l.), middle (2400-2700 m.a.s.l.), and lower (1080-2400 m.a.s.l.) parts of the watershed.Then, tentative soil sampling sites were fixed on the land use and land cover maps along the top sequence of the watershed.In addition, a reconnaissance field survey and discussion with development agents were carried out to get a general overview of the watershed.A total of 27 composite surface soil samples (0-30 cm) were collected using a soil auger in three replications from adjacently located land use systems (agricultural land, grazing land, and forest land) having similar slope gradients and soil types in each stratum (upper, middle, and lower parts of the watershed) (Fig. 2).Ten to fifteen primary samples were collected to prepare one composite sample.Soil quality indicators including soil particle size distribution, ISS, soil water contents, soil pH-H 2 O, organic carbon (OC), total nitrogen (TN), carbon to nitrogen ratio (C:N), available phosphorus (Av.P), Exchangeable bases (Na + , K + , Ca 2+ , and  Mg 2+ ), PBS, and cation exchange capacity (CEC) were considered in this study.Data on the geographical locations of the sampling sites were also gathered using a global positioning system (GPS).

Sample preparation and laboratory analysis
One kilogram of composite samples was prepared, packed in plastic bags, labeled, and taken to Ethiopia Water Works, Design, and Supervision Corporation for analysis.Soil samples were air-dried, crushed, and sieved through a 2 mm and 0.5 mm mesh sieve to remove coarser particles.Standard laboratory procedures were followed in the analysis of soil physical and chemical quality indicators.

Analysis of soil physical properties
Particle size distribution was determined by the hydrometric method (Bouycous 1962).The index of soil aggregate stability (ISS) was determined using Eq. ( 1) as described by (Pieri 1992).The soil water contents, FC and PWP, were determined at 1/3 bar and 15 bar, respectively, by the pressure membrane suction method (Estefan et al. 2013).Then, plant available water holding capacity (AWHC) was calculated using Eq. ( 2).
where ISS = index of soil aggregate stability; OC = organic carbon; L = silt proportion, and A = clay proportion. (1) (2) where AWHC is the available water holding capacity, FC is the field capacity; and PWP is the permanent wilting point.

Analysis of soil chemical properties
Soil pH was determined in water (pH-H 2 O) using a 1:2.5 soil-to-water solution ratio with a pH meter, as outlined in Van Reeuwijk (2006).The Walkley and Black wet digestion method was used to analyze soil OC content (Van Ranst 1993).Total N was analyzed using the Kjeldahl digestion, distillation, and titration method as described by Bremner and Mulvaney (1982).Available soil P was analyzed according to the standard procedure of Olsen et al. (1954).Exchangeable bases (Ca 2+ , Mg 2+ , K + , and Na + ) were determined by ammonium acetate (1N NH 4 OAc) at pH 7.0.Monovalent (K + and Na + ) and divalent (Ca 2+ and Mg 2+ ) cations were determined by using a flame photometer and an atomic absorption spectrophotometer, respectively (van Ranst 1993).The CEC of the soil was determined from ammonium-saturated samples that were subsequently replaced by Na from a percolated sodium chloride solution.The excess salt was removed by washing with alcohol, and the ammonium that was displaced by sodium was measured using the Kjeldahl method (Chapman 1965).
Percent base saturation (PBS) was calculated by dividing the sum of basic cations (Ca 2+ , Mg 2+ , K + , and Na + ) by the CEC of the soil (Eq.3).
where BS% is base saturation in percent (PBS).
where ESP is the exchangeable sodium percentage.

Soil quality deterioration index (SQDI) for each soil quality indicator
The soil quality deterioration index (SQDI) for each soil quality indicator under different land use systems and elevation gradients was computed using an undisturbed ecosystem (forest land) as a reference and evaluated the soil parameters of other land systems against this reference as described by (Abera and Assen 2019;Gui et al. 2009) The model developed by Bajracharya et al. (2007) and also used by other researchers (Feleke et al. 2019) was adopted for this study.In this approach, weighting factors and ratings of selected soil quality indicators were used to develop soil quality ratings (SQR).The ratings for these parameters were adopted from EthioSIS (2014).
The model has the following form: where a, b, c, and d are weighting values for the four soil properties.RSTC = ranking value for soil texture; RpH = value for soil pH; ROM = ranking value for soil organic matter; and RNPK = value for soil nutrient contents of nitrogen, phosphorus, and potassium.These soil parameters are commonly used to evaluate the soil quality index for a given land use system or ecosystem.

Statistical analysis
A two-way analysis of variance (ANOVA) following the generalized linear model (GLM) procedure was applied to test the difference between the mean values of soil quality indicators in different land use systems, the elevation gradient, and their interaction effects.The mean values of statistically different soil quality indicators were separated by Duncan's Multiple Range Test using SAS software version 9.0.In addition, linear correlation analysis was performed to explore the relationships among soil properties.The study area and soil maps were prepared using ArcGIS software version 9.5.

Results
The results of descriptive statistics illustrated the minimum and maximum values of soil quality indicators which revealed the existence of a high degree of variability in the study watershed (Table 2).The differences ( 5) in the land use systems and elevation gradient might be the factors contributing to this variability.The results of ANOVA, correlation analysis, and PCA are presented and discussed sequentially.

Effects of land use systems and elevation gradients on selected soil physical quality indicators
The analysis of variance results of soil physical quality indicators under different land use types and elevation gradients were depicted in (Table 3).The results revealed that soil particles (sand, silt, and clay) showed significant differences (p < 0.05) between agricultural and forest land use types in the upper part of the watershed.In the other cases, non-significant differences were observed in the main and interaction effects of land use systems and locations.The mean value of sand ranged from 2.55 to 35.12%; for the clay fraction, it ranged from 19.85 to 72.66%; and for the silt fraction, it ranged from 24.78 to 47.36%.The textural class of soil in most locations and land use systems is clay; other textural classes include clay loam, silt clay, and sand clay.When the clay fraction is removed from the surface by soil erosion, the content of sand increases, as evidenced by the negative and significant correlation (r = − 0.755, p < 0.01) between these fractions.Significant differences (p < 0.05) were also observed for the Index of Soil Aggregate Stability (ISS) among treatments (land use systems and elevation gradients) and their interaction effects.The lowest (2.06%) and highest (7.21%) values were recorded in the agricultural land and forest land use systems, respectively.Pieri (1992) classified structural degradation of soils as stable structure when ISS > 9%; low risk of structural degradation when 7% < ISS < 9%; high risk of structural degradation when 5% < ISS < 7%; and structurally degraded soil when ISS < 5%.Based on these classifications, soils under agricultural and grazing lands are rated as structurally degraded soils, whereas soils under forest land are rated as having a low risk of structural degradation.
The mean value of soil water content at field capacity (FC) ranged from 36.94 to 59.08% in the clay loam textural class and from 40.31 to 72.41% for the clay textural class, the values of PWP ranged from 25.52 to 43.91% and from 28.84 to 49.28% for the same textural classes.For the above textural classes, the values of AWC ranged from 11.42 to 15.17%; and from 11.47 to 23.13% respectively.Significant differences were observed between the mean values of PWP in the agricultural land and forest land in the midland of the watershed.Though there is no statistically significant difference among most of the soil physical quality indicators, marked differences were observed as they are altered by land use systems and elevation gradients.

Effects of land use systems and elevation gradients on selected chemical soil quality indicators Soil pH, OC, TN, C:N, and available phosphorus (Av. P)
The analysis of variance results revealed that the mean values of soil pH ranged from 5.56 to 6.13 (Table 4), which is classified as moderately acidic (EthioSIS 2014).The lowest value (5.56) and the highest value (6.13) were found in the mid-altitude and low altitudes of the cultivated fields.Non-significant variations were observed between land use systems, elevation gradients, and their interaction effects in the upper and mid altitudes; whereas, in the lower altitude, significant variations were observed between cultivated land and other land use systems.
The mean value of soil organic carbon (SOC) ranged from 1.01 to 3.89%.In all elevation gradients, the lowest values were recorded in the agricultural land, whereas the highest values were in the forest land (Table 4).Significant differences were observed between the main treatments (land use types and locations) and their interaction effects (LU*L).Based on the ratings of EthioSIS ( 2014), the contents of SOC are very low in the agricultural lands in all locations, low in the grazing lands, and moderate in the forest land use system.This indicates that soil quality, particularly soil biological indicators, is highly degraded in the cultivated and grazing land use systems (Fig. 3).
The results also showed that the mean values of soil TN ranged from 0.08% (agricultural land) to 0.21% (forest land) (Table 4).Significant differences (p < 0.05) were observed in the main and interaction effects of land use systems and elevation gradients.A significant difference was observed between agricultural land and the other two land use systems.Based on the ratings of EthioSIS ( 2014), the status of TN is rated as very low in agricultural land and moderate in forest  A significant difference was observed between cultivated land and forest land.Though non-significant variations were observed, the numerical values were higher on the grazing land as compared to the cultivated land.The same is true for forest land and grazing land (Table 4).
The value of available phosphorus was significantly higher (42.33 mg kg −1 ) in the cultivated land, followed by the forest land (21.31 mg kg −1 ) and the grazing land (17.95 mg kg −1 ).Significant differences (p < 0.05) were also observed in the mean values among locations, where the higher value was observed in the lower part of the watershed and the lowest value in the middle part where there is Vertisol coverage (Table 4).The status of available phosphorus is low in the grazing and forest land use systems and moderate in the agricultural land, whereas the overall mean value (27.2 mg kg −1 ) was rated as low (EthioSIS 2014).

and Mg 2+ ), CEC, and PBS
There is no significant difference in all exchangeable bases in the main and interaction effects of land use systems and elevation gradients (Table 5).However, marked variations were observed among land use systems.In all locations, the lowest values of exchangeable Na + and K + were found in the cultivated land, and the highest values in the forest land.The same is true for exchangeable Ca 2+ and mg 2+ .Since the critical value of exchangeable Na + varied due to the variation in CEC, clay content, and type of crop grown, its  2014), the contents of exchangeable bases are rated as high in all locations and land use systems.
The mean values of CEC didn't show significant differences between the main treatments (land use type and locations) and their interaction effects (LU*L), except in the midland of the watershed, where a significant difference was observed between cultivated land and forest land (Table 5).Even though non-significant differences were found in the other cases, the numerical values substantially varied between land use types and locations.The highest overall mean value (74.07 cmol(+) kg −1 ) was found under the forest land, followed by grazing land (68.37 cmol(+) kg −1 ) and the lowest value (65.69 cmol(+) kg −1 ) under the cultivated land (Table 5).Based on Landon's (2014) classifications, the contents of CEC are rated as very high in all land use systems and locations, indicating that soils are inherently fertile due to the parent materials from which they are developed (Elias 2019;William 2016).The highest mean value (86.12%) of PBS was observed in the forest land, followed by grazing land (78.89%), and the lowest value in the cultivated field (74.13%).These values are rated as high and are indicators of good soil fertility.

Soil quality deterioration index (SQDI) for each quality indicator
Soil quality deterioration index results revealed that most of the soil quality indicators have negative values in the cultivated and grazing lands, indicating that the soil condition is becoming poor due to poor management in these land use systems.In the cultivated field, the highest negative index value was for ISS (− 71.3%), followed by SOC (− 67.7%).On the other hand, the highest positive index value for available P (98.6%) was in the cultivated land.For grazing land, the highest negative values were recorded for ISS (− 35.1%) and SOC (− 27.7%) (Fig. 4).The total SQDI values of cultivated and grazing lands were − 12.34% and − 10.0, respectively.

Soil quality index (SQI) computation from production and soil erosion perspectives
The overall SQI mean value under the cultivated land use system was significantly lower than the forest and grazing lands.The mean values were 0.38, 0.60, and 0.52 for cultivated land, forest land, and grazing land, respectively (Table 6).Based on Li et al. 's (2018) classification, the soil quality was rated as low, moderate, and high for cultivated land, grazing land, and forest land respectively.

Effects of land use types and elevation gradients on selected physical soil quality indicators
The relative proportion of separate soil groups determines the soil textural class, which significantly affects other soil properties (physical, chemical, and biological properties).The silt fraction is the most mobile particle due to water erosion, and the proportions of these fractions are lower in poorly managed and intensively cultivated fields.In the current study, the lowest value of silt fraction was found in the cultivated field, followed by grazing land, which implies the problems of soil erosion in these land use systems.On the other hand, the highest value of clay fraction in the cultivated fields might be attributed to accelerated weathering because of soil disturbance.
Soil texture is an inherent property of soil and is the result of soil-forming processes.Its variation could be due to the illuviation of clay particles, the removal of clay by soil erosion, biological activities, upward movement of the sand fraction by swelling and shrinking phenomena, or combinations of these two or more processes (IUSS Working Group 2006).Even though it is less influenced by anthropogenic activities in a short period (Brejda et al. 2000), in the highland region where there is high rainfall and poor land use systems, smaller particle sizes will be eroded and transported to other locations, resulting in variations in the relative proportions of soil particles and textural classes.This illustrates the non-significant differences in soil particle sizes between land-use systems in the study area.This result agreed with some previous research findings.For instance, Mulat et al. (2021) demonstrated that the distribution of soil in separate groups didn't show significant differences among land use systems in the Kersa watershed, Oromia region.Teferi et al. (2016) also reported similar findings from the Jedeb watershed in northwestern Ethiopia.Soil textural class significantly determines the content of soil water, its retention, and the availability of soil nutrients.For instance, the presence of a higher proportion of sand fraction will negatively affect water-holding capacity and nutrient availability (Abera and Assen 2019;Warra et al. 2015).
Soil aggregate stability is considerably influenced by land use systems and land management strategies.In this study, it was negatively affected in cultivated fields compared to soils under other land use systems.This could be attributed to the depletion of soil organic matter in the cultivated fields and the better accumulation of SOM in the forest and grazing lands.The findings of Delelegn et al. (2017) also explained that the value of soil aggregate stability in cultivated land is significantly lower than in other land use systems.Wei et al. (2014) also emphasized that intensive cultivation significantly decreases soil aggregate stability, whereas, in non-tilled landscapes (forest lands), macro-aggregates of soils are enhanced.Soil organic matter is a key factor in binding soil particles into aggregates.Due to tillage operations, soil organic matter could be further depleted, and hence macro aggregates will be broken down in the cultivated fields.The presence of an adequate amount of SOM in the soil enhances soil microbial and fungal biomass, which will improve soil aggregate stability (Wu et al. 2015).Mulat et al. (2021) reported higher soil aggregate stability in grazing land use as compared to fallow land and cultivated fields.Peng et al. (2016) also explained the role of soil organic matter in maintaining soil aggregate stability.Similarly, Gebreyesus et al. (2014) examined the negative impacts of soil organic carbon depletion on soil aggregate stability and its consequences for soil degradation.Devine et al. (2014) emphasized that undisturbed ecologies like soils under forest land have better aggregate stability than disturbed ecologies.
Exploring the status of soil water content is essential to devising an irrigation schedule and improving water use efficiency depending on the water requirements of plants.In the present study, the highest amounts of soil water contents (FC and PWP) were obtained in the cultivated field, and the lowest values were obtained in the forest land with different land use systems and elevation gradients.This variation could be attributed to the contents of the clay fraction and soil organic matter (SOM).This is also evidenced by the positive and significant correlation of clay with FC (r = 0.441, p < 0.05) and PWP (r = 0.932, p < 0.01) (Table 7).The variation in the contents of available water holding capacity (AWHC) among different land use systems could be due to the difference in the contents of clay fraction and soil organic matter (Alawamy et al. 2021).

Effects of land use systems and elevation gradients on selected chemical soil quality indicators Soil pH, SOC, TN, C:N, and available phosphorus (Av. P)
Even though the values of soil pH didn't show significant differences among land use systems except at location 3 (lower elevation), the soil condition is moderately acidic in all land use systems and locations.Even though the study area is found in the northwestern highlands of Ethiopia, where there is high rainfall and rocky topography, the problem of soil acidity is not critical.This might be attributed to the content of SOM and high CEC, in which basic cations are adsorbed on the surfaces of these colloidal particles.This is also evidenced by positive and significant correlations of pH with CEC (r = 0.456, p < 0.01), PBS (r = 0.425, p < 0.05), exchangeable Ca 2+ (r = 0.486, p < 0.01), exchangeable K + (r = 0.491, p < 0.01), and exchangeable Mg 2+ (r = 0.490, p < 0.01) (Table 7).On the other hand, the lowest values of soil pH in some locations and land uses could be due to the removal of basic cations through soil erosion, leaching, application of acid-forming fertilizers, crop harvest, and residue in the cultivated fields, and leaching and removal through soil erosion in the forest and cultivated lands, respectively.This result is in agreement with the findings of Teferi et al. (2016), who reported a non-significant difference among different land use systems in the Jedeb watershed, in the northwestern highlands of Ethiopia.On the contrary, Mulat et al. (2021) found significant differences among land use systems in the Kersa watershed, Oromia region.These differences could be due to differences in soil types, climatic conditions, and management strategies.Soil pH indicates the toxicity level of aluminum (Al) and the status (deficiency) of soil micronutrients.Soil organic matter (SOC) is an important soil parameter that highly influences soil's physical, chemical, and biological properties.It is the key indicator of soil quality and is significantly influenced by anthropogenic activities, vegetation cover, and climatic conditions (Yu et al. 2020).In this study, the lowest values were found in cultivated lands in all elevation gradients, and the values are far below the critical level.The content is also low in the grazing land but optimal in the forest land.Further decomposition of SOM due to continuous and intensive cultivation in the cultivated fields, depletion through soil erosion, complete removal of crop residues from the field (sources of SOM), and low input of organic fertilizers (crop residues and manure) are the major factors for extensive depletion of SOM in the cultivated fields.This leads to a high loss of soil nutrients from these fields, as the SOM is the most important soil parameter that determines soil fertility.On the other hand, low soil erosion, continuous litter fall to the soil, and the microclimatic conditions under the forest land (low organic matter decomposition) contribute to low soil organic carbon loss and the buildup of SOM in the forest land.These results agreed with the findings of Mulat et al. (2021), who reported a higher content of SOM in the grazing land than in cultivated and fallow lands in the Kersa watershed in eastern and Northern Ethiopia.Feleke et al. (2019) and Ryan et al. (2018) also demonstrated a higher content of SOM in protected grassland than in open grassland due to the frequent turnover of grass and dense root biomass in the protected land use.Similarly, Elias (2019) emphasized the negative impacts of the complete removal of crop residues from cultivated fields and using cow dung as a source of fuel instead of as soil amendments on the loss of soil nutrients and SOM.
Land use and land management strategies significantly affect the content of soil TN in the study area.In poorly managed cultivated fields, the content of TN is significantly lower than in other land use systems.This variation might be attributed to the amount of soil organic matter in the forest and grazing land, the depletion of SOM, and the leaching of N in the form of nitrate ions (NO 3− ) from cultivated fields.The amount of SOM present in the soil is a key factor for the presence of TN in adequate amounts, as evidenced by a significant and positive correlation of TN with SOM (r = 0.908, p < 0.01).Similar findings were reported by Feleke et al. (2019), who demonstrated that the content of TN in protected grassland was higher than in unprotected grassland because of the variation in the contents of SOM in these land use types.
In this study, in all land use systems and locations, the mean values of C:N ratios were less than 20:1, indicating the mineralization process.It also demonstrated the potential impacts of different land use systems and elevation gradients on the C:N ratio.The lowest value was found in the cultivated land in the upper part of the watershed.This could be attributed to further decomposition and removal of SOM as a result of intensive cultivation and soil erosion.On the other hand, the highest mean value was found in the forest land of the same elevation due to its high SOC content.This result is in agreement with Seifu et al. (2020), who explained that the highest C:N ratio is in the grazing land where there is high SOM and the lowest value is in the bare land.
The content of available phosphorus could be influenced by land use systems, soil management strategy (application of external inputs), and the amount of soil organic matter.The highest mean value of available phosphorus was observed in the agricultural land as compared to other land use systems in the study area.This might be attributed to long-term applications of phosphorus-containing fertilizer (DAP), which was confirmed during our field survey work and data gathered from farmers and experts.On the other hand, the lowest value in the grazing land might be due to the low turnover of grass biomass to the soil because of overgrazing.This in turn results in lower soil organic matter and organically bound phosphorus, which could be released into the soil solution through mineralization.In forest land, available phosphorus is higher than in grazing land because of better accumulation of SOM.Parallel with this result, Seifu et al. (2020) reported the highest mean value of available phosphorus in the cultivated land as compared to the forest and grazing land systems in the semiarid watershed in northwestern Ethiopia.Abera and Assen (2019) also reported similar results from the Wanka watershed in northwestern Ethiopia.Likewise, De et al. (2022) and Tellen and Yerima (2018) explained the higher content of available phosphorus and nitrogen in cultivated land because of the continuous application of chemical fertilizers.The high availability of phosphorus is an indication of low p-fixation, which is highly dependent on soil pH (availability decreases in both acidic and alkaline soil conditions).The content of SOM also determines the availability of phosphorus.The findings of other researchers (Stephens et al. 2014;Feleke et al. 2019) indicated that phosphorus is positively and significantly correlated with SOM and that its availability is dependent on the mineralization of SOM.However, in the current study, phosphorus was significantly but negatively correlated with SOM.The reason might be due to the continuous application of phosphorus fertilizer in the cultivated field; conversely, the lowest value of SOM was found in the same field.

Exchangeable bases (Na + , K + , Ca 2+ , and Mg 2+ ), CEC, and PBS
The lowest values of exchangeable bases except Mg 2+ were found in the cultivated land, and this might be due to intensive cultivation, leaching of basic cations, soil erosion, low SOM content, and removal of cations with crop harvest from this land use type.On the other hand, the highest value of CEC was in the forest land, which could be linked with the amount of SOM.Exchangeable K + and exchangeable Ca 2+ were positively and significantly (r = 0.645, p < 0.01 and r = 0.459, p < 0.01) correlated with SOM; exchangeable Na + and Mg 2+ also showed a positive correlation with SOM (r = 0.154 and r = 0.060).From the results of exchangeable bases, it is observed that the exchange sites of the colloidal particle were dominated by divalent cations (Ca 2+ , and Mg 2+ ), According to Bohn et al. (2001), for agriculturally productive soils, the order of exchangeable bases should be Ca 2+ > Mg 2+ > K + > Na + , and any deviation from this could cause a cation imbalance for plants.Another important point is the ratio of Ca to Mg cations (Ca:Mg).For most crops, the optimal Ca:Mg ratio should be between 3:1 and 4:1 (Landon 2014).If the ratio is less than 3:1, it will inhibit the uptake of phosphorus.
In the present study, the higher content of CEC in forest land followed by grazing land might be attributed to the content of SOM.This could be evidenced by a positive (r = 0.31, p > 0.05) correlation of CEC with SOC.In agreement with this result, Abera and Assen (2019) demonstrated that the CEC in the natural forest land was higher than in grazing and cultivated lands because of the high SOM content in the forest land.The content of CE also varied along the top sequence of the watershed.The highest value was found at the lower part; this could be due to the large coverage of Vertisol and Leptosol (soil distribution map), with clay minerals having a high surface area.This could be substantiated by a positive and significant correlation (r = 0.548, p < 0.01) of CEC with clay fraction (Table 7).In line with this result, Elias (2019) reported high CEC in Vertisols and Leptosols, which are dominated by smectite minerals with high surface areas and pH-dependent charges, in the cultivated fields of the Ethiopian highlands.

Soil quality status under different land use systems
The soil quality of the Suha watershed is declining because of anthropogenic activities, as evidenced by the negative values of SQDI for most soil attributes under cultivated and grazing lands.In these land use systems, ISS, SOC, TN, and exchangeable K are the most deteriorating soil quality indicators.This might be due to soil erosion, leaching, depletion of SOC as a result of intensive cultivation, low returns of organic inputs, and poor management.On the other hand, high positive values of available phosphorus and clay fraction in the cultivated field could be because of the regular application of chemical fertilizer and cultivation, which enhances the weathering process and hence increases clay fraction.The total SQDI values also showed that soil conditions under cultivated and grazing land use systems have been experiencing progressive deterioration as a result of poor management strategies.These results agreed with the findings of Abera and Assen (2019), who reported declining soil quality indicators in disturbed ecosystems (agricultural land and grazing land) in the Wanka watershed, northwestern Ethiopia.Feleke et al. (2019) also reported that a higher value of SQI was found in protected grassland and a lower value in unprotected grassland in Farta District, northwestern Ethiopia.Other researchers (Eyayu et al. 2009;Gui et al. 2009) also reported similar findings.

Conclusions
The present study was conducted to evaluate the impacts of land use and cover change on soil quality in terms of its physical and chemical indicators along the top sequence of the Suha watershed in north-western Ethiopia.The results illustrated that land use and land cover change (particularly from forest land to agricultural land) had detrimental impacts on soil physical and chemical quality indicators that resulted in the deterioration of soil conditions and a reduction in agricultural productivity.Soil quality indicators that are sensitive to land use change and management are below their threshold values, indicating that soil condition is deteriorating.Poor management of agricultural lands is the main factor in the decline of soil quality indicators, which could be explained by the depletion of soil organic carbon (SOC) and nitrogen deficiency.In addition, the soil quality degradation index values of most soil quality indicators are highly negative in the cultivated field, indicating the existence of poor management and deterioration of soil conditions.Moreover, aggregated soil quality index values in the agricultural land use systems are below the optimum range, which is evidence of poor management in these land use systems.Therefore, soil and site-specific management strategies are required to reverse degraded soil conditions and maintain soil quality indicators to improve soil ecosystem services, enhance agricultural productivity, and ensure food security.In this regard, the application of integrated soil fertility management (ISFM) strategies and soil conservation technologies are suggested as remedial actions.Land use and land cover change are two of the major factors that impact soil quality; other factors such as soil types, slope classes, and soil management strategies should be explored to design and implement holistic and effective soil management strategies.

Fig. 1
Fig. 1 Map of the study area

Fig. 2
Fig. 2 The general framework of the methodology used

Fig. 4
Fig. 4 SQDI values of soil quality indicators under cultivated and grazing lands based on the referenced undisturbed ecosystem (forest land).SQDI soil quality deterioration index, CL cultivated land, GL grazing land

Table 1
Land use and cover types and their area coverage in the Suha watershed

Table 2
Descriptive statistics of the analyzed soil physicochemical properties in the Suha Watershed

Table 3
Mean values of sand, clay, silt, ISS, FC, PWP, and AWC under different land use systems and elevation gradientsThe mean values for the same letter are not significant

Table 4
Mean values of pH, OC, TN, C:N, and Av.P under different land use systems and elevation gradientsThe means in the rows followed by the same letter are not significant The relationship between OC and ISS under different land use systems.OC organic carbon, ISS index of soil aggregate stability, CL cultivated land, GL grazing land, FL forest land and grazing land use systems.The mean values of the C:N ratio ranged from 11.22 to 18.71, and the overall mean values were 12.98, 14.77, and 17.5 in the cultivated land, grazing land, and forest land, respectively.

Table 5
Mean values of exchangeable bases (Na + , K + , Ca 2+ , and Mg 2+ ), CEC, and PBS under different land use systems and elevation gradients Means in the same row followed by the same letter are not significant (Anon 1954value was more explained by ESP, in which an ESP value of 15% is critical for most crops(Anon 1954).The lowest (1.92%) and highest (2.12%) values of ESP were observed in the grazing and forest land use systems, respectively.The soils of the study area are free of sodicity problems as the values of ESP are far less than 15%.Based on the classifications of EthioSIS (

Table 6
(Li et al. 2018)uality under different land use systems based on defined ratings(Li et al. 2018)CL cultivated land, FL forest land, GL grazing land, SQI soil quality index, STC soil texture, OC organic carbon, NPK nitrogen, phosphorus, and potassium

Table 7
Pearson's correlation matrix of selected physical and chemical soil properties