Abstract

Understanding topography effects and assessing the soil properties in different land use is an essential first step for sustainable soil management. Hence, land use type and altitudinal gradient on selected soil parameters were studied in Ayiba watershed, northern Ethiopia. Thirty composite soil samples were collected from 0 to 30 cm of soil depth under four land use types across three altitudinal gradients and were analyzed for selected soil parameters following the standard procedures. A significant main effect of land use and altitudinal gradient on the content of the soil particles was noted. Results also indicated that the bulk density (BD), total porosity (TP), and Pav of the soil are significantly different ( < 0.05) in the watershed because of land use type and altitudinal gradient. Barren land and higher altitude landscapes have the highest BD, and the lowest TP, and grassland and lower altitude landscapes have the lowest BD and high TP. Intensive cultivation accompanied by natural land conversion and erosion due to the rugged landscape nature caused high BD and low TP. Analysis of variance results also shows the significant interaction effect of land use type and altitudinal gradient on EC, SOM, SOC ( < 0.001), and pH-water (1:2.5), and TN ( < 0.01). The Pearson correlation of SOM with TP, TN, MC, and clay content showed a strong positive relationship. But, SOC, TN, and clay content were negatively correlated with BD. Soils of the study watershed are found in low to optimum rating levels in their selected physicochemical properties. Overall, the results show that land use and topography gradient significantly affected soil physicochemical properties in the study watershed. Therefore, soil management options should focus on scenarios that could improve the soil conditions to enhance crop production on a sustainable basis.

1. Introduction

The optimally integrated coexistence of principal natural assets (viz., soil water and vegetation) is vital for maintaining terrestrial ecosystem functions and productivity [1]. Remarkably, the soil, an essential resource of nature and a significant component of the terrestrial ecosystem, plays a vital role in upholding the trio known as NPK and other plant nutrients (https://www.dpi.nsw.gov.au). However, the base of this vitally important natural capital is becoming vulnerable to declination and deteriorating over time [13] and manifested in land/soil degradation. The deriving agent includes high population pressure, excessive livestock pressure of grassland areas, changes in land use, inappropriate agricultural practices, massive deforestation of natural forest (farm expansion, fuel, and construction), agricultural expansion (cultivation of steep slopes and marginal areas), and rugged nature of the landscape [48].

Soil properties varied significantly among soil types and across locations reflecting differences in parent materials, climate, and land use [4]. The northern region and the most highland parts of Ethiopia are among the most seriously affected areas by land/soil degradation [912]. Factors are aggravating land/soil degradation resulting in low soil fertility manifested in lower yield and higher environmental problems [13]. Bekunda et al. [14] specially narrated this situation as “the soil fertility decline is not out of the smallholder farmers’ violation, but rather the consequence of striving for household well being under difficult circumstances.”

The food insecurity and rural poverty in the African smallholder farming system are due to soil nutrient depletion. Pomeroy et al. [15] observed a remarkable decline in soil nutrients due to the deterioration of chemical, physical, and biological properties in East Africa. The main reasons for the decline, besides soil erosion, are a decline in soil organic matter and soil biological activity, degradation of soil structure and loss of other soil physical qualities, reduction in the availability of significant macronutrients (N, P, K) and micronutrients, and increased toxicity due to acidification and salinization [16]. A previous study also confirmed the fact of this negative soil nutrient balance, mainly soils organic carbon by 0–63% following deforestation in sub-Saharan Africa [17] and macronutrients of (N; −22, −122 kg ha−1), (P; −2.5, −13 kg ha−1), and (K; −15, −82 kg ha−1) in Africa and Ethiopia, respectively (the numbers in parenthesis). The negative soil nutrient balance indicates that the primary soil nutrients are already mined exhaustively [18, 19] because of the limitation on sustainable management practices based on the soil system’s understanding. Above all, soil properties’ overall status in a given landscape is controlled mainly by the type of land use system applied. The demand for arable land, grazing, forestry, wildlife, tourism, and community development is more remarkable than land resources available, so that soil quality may worsen [20]. Hence, the soil needs to be protected sustainably without compromising current or future generations’ ability to meet their needs.

Characterizing the spatial variability of soil nutrients concerning site properties, including climate, land use, landscape position (topography), and other variables, is vital for a comprehensive understanding of how the ecosystems work and assessing the effects of future land use change on soil nutrients [21, 22]. Henceforth, under a small watershed scale, topography and land use types may be the dominant factors determining soil properties’ status [23, 24], particularly in areas of multiple land use and complex topography. Accordingly, quantifying the pattern and spatial distribution of basic soil indicators is fundamental to understanding many ecosystem processes. For example, determining the spatial distribution of soil organic carbon is an essential factor in land management decision making, climate change mitigation, and landscape planning [25]. Therefore, landscape-based management, protection, and soil restoration are essential for several environmental and policy reasons.

Topography and land use types are vital factors affecting soil property variability [26]. The topography is a significant factor controlling both hydrological and soil processes at the landscape scale, and its influence is also apparent in the soil catena concept [27]. Environmental factors such as slope aspect induced by microclimate differences, topography, parent materials, and vegetation communities are significantly influenced by the spatial variation of soil properties [28]. Amuyou and Kotingo [29] reported that slope gradients have a marked influence on soil properties expressed in the soils’ distribution along with slope positions. Soil moisture and soil chemistry could affect plant growth, and their distributions on specific ecosystems are governed by topographic attributes [27, 30]. Concerning topographic positions, significant variation for clay, bulk density, organic matter, and total nitrogen was observed [26] in China’s loess hilly area. According to Ezeaku and Eze [31], higher CEC, OM, and TN values were recorded along with the slope compared to the upslope position. Many other authors also reported the influence of topography on soil texture, organic carbon, and CEC [3234] in Africa. Landscape features in most parts of Ethiopia influence soil resources, including depth; thus, cultivation practices aggravate reduced crop production because soil erosion is a function of steep slopes and land use in mountainous areas. Hence, there is a need to improve these tendencies, which lies in the improvement of the agricultural sector’s long-term development.

Besides, changes in land use and its alterations on the ecosystem can significantly affect soil physical, chemical, and biological properties [35, 36]. Beyene [37] reported that, in southern Ethiopia, soils of cultivated lands have higher available phosphorus (P) content, whereas exchangeable bases, OC and CEC, were higher in grasslands. Reduction of OM, available phosphorus, and porosity was observed due to the change of land use from forest to cultivation, while a reverse trend was examined for BD and pH [35, 38]. In a meta-analysis of soil carbon sequestration in agroforestry systems, SOC stocks were found significantly reduced when a pasture/grassland was converted into an agrisilvicultural system. Significant SOC stocks increases were also reported at various soil horizons and depths in the land use change from agriculture to agrisilviculture, forest plantation to silvopasture, and uncultivated/other to agrisilviculture [39]. The study conducted by Hu et al. [26] revealed that forest land use has higher soil OC, TN, and CEC than a cultivated land use system. Many investigations indicated that SOC and other plant nutrients are the soil properties most affected by cultivation over time. Thus, the decline in yield has been attributed to the loss of plant nutrients by plant removal, erosion, leaching, and deterioration of soil physical condition [18]. Majule [40] further justified that soils in areas with continuous cultivation and without appropriate management practices have low fertility levels due to overutilization. The linkage between land use land cover and soil fertility in Kenya and Tanzania across altitudinal gradients was assessed, and higher phosphorus (P) was reported in the lower zones than the middle and upper zones [15, 41]. Cleaning the natural forest in most parts of the upper zones for cultivation has contributed significantly to the reduced level of SOC, N, P, and K in the soil [40].

Soil nutrient depletion and related low agricultural productivity are severe problems in Ethiopia and has forced farmers to abandon their fields and convert into woodlots dominated with Eucalyptus spp. Agricultural activities such as tillage, fertilization, and biomass alteration lead to fundamental changes in the pools and fluxes of carbon (C), nitrogen (N), and phosphorus (P) that initially existed in native ecosystems [42]. Besides, land use change, mainly through the conversion of natural vegetation to cropland and grazed pastures, may influence many ecological properties [43], such as soil carbon dynamics and soil microbial properties [44]. Moreover, land cover changes and deforestation associated with agricultural expansion are among the most significant contributors to anthropogenic carbon emissions worldwide [45]. Thus, in decision making on sustaining agricultural production and land degradation mitigation, a local inventory of soil properties to understand the spatial distribution is essential. A reasonable way of deriving the evidence is using proxy environmental features that have established relationships with the soil properties (e.g., land use and topography gradient factors used in this study). Spatially explicit data on land use and topography effects on soil properties are limited in the study watershed. There is, therefore, a need to assess the status of soil physicochemical properties between landscape positions of each land use types for sound management of land resources and enhance sustainable crop productivity at the watershed level. Therefore, the main objective of this study was to (i) examine the influence of different land use types on soil physicochemical properties at landscape scale and (ii) evaluate selected soil physicochemical properties along the toposequence (lower-middle-upper sequence).

2. Materials and Methods

2.1. Description of the Study Area

Ayiba watershed is part of the Denakil River basin and geographically located between 12°51′18″–12°54′36″N and 39°29′24″–39°35′24″E (Figure 1). It covers an area of about 4099.14 ha. The elevation of the watershed ranges from 2722 to 3944 meter above sea level. The area’s topography is a mountainous and hilly dissected terrain with steep slopes and undulating topography, which characterizes the middle and upper parts, whereas gentle slopes form the lower parts of the study watershed. The soil parent materials of Emba Alaje District (southern Tigray), where the study watershed is situated, are of igneous rock of primary class (I), basic igneous group (IB) of basaltic type (IB2), which are deeply weathered [46]. The major reference soil groups (FAO guideline based) of the study area are Leptosol, Vertisols, Cambisols, and Fluvisols [46]. The area’s climate is generally characterized as a subtropical continental type with an extended dry period of nine to ten months, and the maximum rainy days last for 50–60 days [46]. The watershed’s 20 years of weather data were estimated from four active nearby weather stations (Bora, Maychew, Wedisemero, and Korem) of the Ethiopia National Meteorological Agency. Accordingly, the study area receives a mean monthly rainfall of 72.88 mm and a total annual rainfall of 875 mm. The rainfall pattern in the area is bimodal, with a high amount of rainfall occurring during the primary rainy season (locally called “Kiremti”) lasting between July and September, and the short rainy season extends from March to June. April is the peak period of the small rain (59.94 mm), and August is the peak for the “Kiremti” season (262.54 mm). Similarly, the mean minimum and maximum monthly temperatures of the area are 7.1 and 25.6°C, respectively, with a mean of 16.8°C (Figure 2).

2.2. Current Land Use Cover

Land use and land cover change data provide essential information for environmental management and planning [48]. The current land uses are generally rainfed agriculture and livestock (grazing) dominated by degraded shrublands [46]. Landsat_8 OLI-TIRS-30 m resolution was derived from open access Landsat imagery service (https://earthexplorer.usgs.gov/), acquired with path 168 and row 51 in the dry season (28 February 2017), when spectral differences between land cover types are most significant, and cloud contamination is minimal. To assure that all pixels were precisely aligned and conducted image-to-image coregistration before the classification in the data preprocessing using ENVI software (Exelis Visual Information Solutions, Boulder, Colorado). Data were projected to the Universal Transverse Mercator (UTM) projection system zone 37N and data of World Geodetic System 84 (WGS84), ensuring consistency between datasets during analysis.

The watershed boundary was delineated using Arc hydro tool, and finally, the land use/land cover map of the study watershed was produced in ArcGIS 10.5 environ. Accordingly, four major land use types (barren land, cultivated land, grassland, and plantation forest) of the watershed were identified (Figure 3, Table 1). To generate a training data set for our supervised classification, we identified pure pixel patches of each land use type by visually interpreting Landsat images and high-resolution imagery from Google Earth.

2.3. Soil Sampling and Laboratory Analysis

Before soil sampling, an interview with residents and field observation was conducted to investigate the watershed land use history. Land use types and elevation variations were surveyed before soil samples were collected. Three composite soil samples prepared following the coning and quartering method [49] (0–30 cm depth, each sample consisting of fifteen subsamples) were collected using a hand auger (100 cm3) for each land use type across topography in the watershed in 2017 offseason from 20 January to 30 February. The three topographic gradient categories are lower (2722–2800), middle (2800–3100), and upper (3100–3944), respectively. A global positioning system (GPS) was used to identify the site’s longitude, latitude, and elevation.

Soil sampling sites were selected by stratifying the watershed. Three random strata were aligned vertically along the toposequence from the bottom to the top edge elevation of the mountain, and major adjacent land use types in each stratum across the altitude were identified (four at upper and middle elevation (BL, CL, GL, and PF) and two at a lower elevation (CL and GL). The soil sample was taken 10 m before and 10 m after the border between adjacent land use and ∼150 m away from the outer ridge to avoid edge effects. At each sampling spot, surface litter is scraped, vegetation cover is removed before collecting samples, and sample spot excluding was performed if the considerable difference is noted to minimize soil variability. For bulk density, the cylindrical core method [50] were used to collect undisturbed samples, and a total of 30 undisturbed soil samples were collected by taking three core samples from each land use type along with each topographic class. The samples were air-dried, ground, sieved to remove gravel fractions (>2 mm), weighed, and prepared as required for final laboratory analysis [51]. Laboratory work was carried out at the Tigray Soil Laboratory Centre, Mekelle, Ethiopia, and at the plant nutrition laboratory, college of environmental science resources, Zhejiang University, China.

The chemical analysis for organic carbon was performed by drying samples in an oven at 105°C and subsequently burning using the loss-on-ignition (LOI) method at 400°C [52] to determine soil organic matter. Then, SOM was converted into soil organic carbon according to equation (2) by multiplying the obtained estimated SOM from equation (1) with the “van Bemmelen factor” of 0.58, as described in the work of Brady and Weil [53] and Guo and Gifford [54].where SOMLOI: soil organic matter of loss on ignition, W105°C: the weight of soil before combustion (preignition) after oven-dried at 105°C and W400°C: the weight of soil after combustion (postignition) at 400°C, and SOC: soil organic carbon.

2.4. Statistical Analysis

Soil data obtained from the laboratory work were analyzed using SPSS (version 26 for Windows). Descriptive statistical analyses were conducted using this software. Two-way ANOVA was used following the general linear model (GLM) procedure, where land use types and landscape positions (altitude) were independent variables (factors) and the selected soil physicochemical properties were dependent variables. Post hoc tests for each variable were performed when the variance analysis showed statistically significant differences at  0.05 level of probability. For statistically different parameters at probability 5%, means were separated using Duncan’s test. Pearson correlation analysis was conducted to determine the correlation between all the measured parameters with R “corrplot package” [55].

3. Results and Discussion

3.1. Statistical Analysis of the Selected Soil Physicochemical Properties

The descriptive analyses revealed a numerical difference among the selected physicochemical soil parameters among land use types across altitudinal gradients. The variance analysis also showed significant statistical variation among the selected measured soil physicochemical parameters due to variation in land use type and altitude factors. As shown in Table 2, there was a significant main effect variation in terms of land use and altitudinal gradient on sand fraction and soil moisture (< 0.001), in terms of land use alone on silt fraction (< 0.011) and clay fraction (< 0.01), whereas BD, TP, and Pav (< 0.05), pH-H2O (1 : 2.5) (< 0.01), EC, SOM, SOC, TN, and C : N (< 0.001) were found strongly dependent on land use, topographic positions, and their interactions. However, Kexch did not show any significant difference with land use (> 0.05), altitude (> 0.05), or their interaction (> 0.05) (Table 2).

3.2. The Effect of Land Use and Altitudinal Gradient on Selected Soil Physical Properties
3.2.1. Soil Particle Size Distribution

The variance analysis revealed that all the land use types’ textural class was clay, indicating the similarity in the watershed parent material. The mean sand content ranged from 17.93 ± 1.63% to 29 ± 6.55% across altitude and from 22.95 ± 1.03% to 37.24 ± 1.26% among land use types. Along with the altitudinal range, high clay content was noted. Across altitude, mean silt and clay contents were found in the range of 25.35 ± 45 to 30.73 ± 2.01 and 44.68 ± 4.81 to 51.34 ± 3.07, respectively (Figure 4), with the highest record at GL and lowest record at BL. The highest clay content at GL may be due to the change in land use from the communal-free intensive grazing system to a partially (periodical) protected grazing system and managed cut-and-carry system (personal communication with development agents, 2018). This practice positively contributes to accumulating the fine clay particles in GL. Barren land showed a higher sand percentage, which may be due to vegetation clearance and erosion. Elias [57] confirms that although the texture is an inherent soil property, land use practices can be altered over a more extended period. Other studies also reported that silt and clay particles were found to be significantly ( ≤ 0.05) affected by land use types in Kuyu district, Ethiopia [58].

Relatively high sand contents were observed at the upper topography compared to the lower topography, and along with the soil catena, the distribution of sand fraction was in the following order: upper altitude > middle altitude > low altitude (Figure 4), indicating that sand content decreases down toposequence. Because erosion is a selective process that can remove clay particles leaving the coarse fragments behind and sedimenting downstream (Figure 5), likewise, others observed a decrease in fine fractions and a corresponding increase of sand content in the steeper slope due to the selective removal of fine particles by water [32, 59]. In Ethiopia and elsewhere, similar results have been reported as significantly higher fine soil textural fractions in the low slope gradient than the upper ones [57, 60, 61]. Many other authors also support this result, revealing that deforestation, farming practices, and intensive grazing systems have changed soil texture by aggravating soil erosion [62, 63]. Charan et al. [64] added another possible justification as sand fraction may increase with elevation due to the abundance of sand-forming minerals in parent material at the higher topographic position.

3.2.2. Soil Bulk Density, Total Porosity, and Moisture Content

The mean value for soil bulk density (BD) varied from 0.98 to 1.40 g cm−3, with the highest mean value recorded in the upper topography and the lowest in the lower topography under barren land (BL) and grassland (GL), respectively (Table 2). High BD values in BL might be attributed to deforestation, compaction by animals trampling, continuous erosion, and low organic matter input. The lowest BD at the lower altitude than the other landscape position might be due to managed use of the GL area (cut-and-carry system) since 2004. This practice decreases intensive free grazing and livestock trampling effect and positively contributes to reducing BD due to its relatively highest soil organic matter input from above-ground and below-ground grass biomass. Besides, erodible mineral deposits are transported principally by water as the erosive agent from upper and middle streams, which reduces BD in the lower altitude. The cut-and-carry system is also clearly attributable to the high addition of soil organic matter, which increased soil organic carbon and decreased BD of the surface soil reported by previous studies [1, 65, 66]. Similarly, the lowest BD in GL due to the higher clay content and less disturbance of the soil under grassland was reported [58]. Contrary to this argument, Lelisa and Abebaw [67] reported higher BD in the GL than the adjacent BL and rehabilitated land in the Jigessa rehabilitation site, Borana zone, Ethiopia.

The higher BD in GL of the middle topography than the upper topography indicates that soil compaction was enhanced due to free overgrazing pressure, and livestock disturbance over trampling, especially during the rainy season. The result of this study is in correspondence with the result of the work of Pei et al. [68] and others who revealed BD often increases in soils under heavy grazing pressure by animal trampling, poaching, and pugging [66, 6972]. Recently, Shiferaw et al. [73] added that the lower BD in enclosure land use type resulted from fewer disturbances, higher litterfall, and organic matter accumulation, which indicated better management that led to better soil organic carbon sequestration. Raiesi and Riahi [74] also reported that heavily grazed rangelands had a more significant (22%) BD than ungrazed rangelands in central Iran. Moreover, Mekuria et al. [75], who studied exclosures’ role to restore degraded soils in Tigray due to overgrazing, stated that exclosures lead to a restoration of natural resources such as soil fertility, vegetation biomass and composition, fauna, and water storage. The topsoil BD of the study watershed was found below the typical value suggested by Hazelton and Murphy [76] for agricultural soil, signifying that there is less soil compaction or restriction for root development. Generally, in the present study, an increasing BD trend was observed as elevation increases and followed the Upper > Middle > Lower order (Table 3). Contrary to this finding, Ahmad and Somaiah [77] have reported a decreasing bulk density trend with increase in altitude.

The mean values of total porosity (TP) under different land use types, namely, BL, GL, cultivated land (CL), and plantation forest (PF), were 48.56, 54.55, 57.08, and 55.45%, respectively, and the mean value of TP across topographic classes, namely, upper, middle, and lower, was 51.79, 53.46, and 60.96%, respectively (Table 4). The mean values for TP across the watershed varied from 47.44 to 63.06%, with the highest TP in the lower landscape position under GL (63.05%) and the lowest TP in the upper landscape position under BL (47.44%) (Table 3). These differences might have occurred due to continuous livestock movement for grazing purposes in grassland soils and due to erosion, which exhibits higher BD and lower porosity. The highest TP at lower topography, in PF at the middle topography, and in GL at upper topography could be accredited by high organic matter contents and low BD. Similarly, Gupta [78] and Brady and Weil [79] revealed that the number of pore spaces and soil organic matter is inversely related to BD, and all factors that affect soil pore spaces will also have effects on BD. Hence, a decrease in porosity may be due to a reduction in organic matter content.

In this study, the soil moisture of GL recorded the highest mean value, which was statistically at par with CL and PF, whereas the soil of BL had the least mean value (Table 3). The mean soil moisture content of the studied land uses of the watershed followed GL > CL > PF > BL. Because of the topographic gradient, soil moisture was higher at lower topography (Table 3), which may be ascribed to higher soil clay fractions with possibly good soil structural aggregates and high organic matter accumulating soil surface. The result agrees with the work of Elias [57] and Chimdi et al. [80], who reported higher water holding capacity for soils with higher clay content in Nitisol profiles of some highlands of Ethiopia and Welega, western Ethiopia, respectively. The finding of this result is substantiated by Yimer et al. [81], who reported higher soil moisture content under exclosures in the grassland area than in open grazing land use in Ethiopia’s central rift valley area. Besides, low soil moisture was recorded in middle topography, which may be due to soil degradation, intensive cultivation and free open grazing. Moreover, the post hoc test (Bonferroni) results showed that BD, TP, and soil moisture were found significantly (< 0.05,  < 0.01, and < 0.001, respectively) varied between the different land use types and topographic positions (Table 5).

3.3. The Effects of Land Use Type and Altitudinal Gradient on Selected Soil Chemical Properties
3.3.1. Soil Reaction (Soil pH) and Electric Conductivity (EC)

The overall mean soil pH-H2O ranged from 6.48 ± 0.08 to 7.31 ± 0.08 (Table 6), which is rated as moderately acidic to neutral (Table 4), according to Karltun et al. [56]. But, there was marked variation among land use types across altitudinal gradients (Table 6). The mean EC is 0.20 ± 0.01 mS cm−1 with the highest mean value of 0.37 ± 0.01 cm−1 at grassland, rated as salt free. Comparatively higher soil pH was measured in grassland across all topographic gradients (Table 6), which might be due to the increase in surface soil organic matter, which traps base cations. Additionally, the accumulation of bases that are supposed to be moved laterally by erosion and sediment on grassland and cattle dung wastes deposited on grass fields. The finding was in agreement with Mekuria et al. [75] who justified organic matter accumulation might reduce soil erosion resulting in higher soluble base cations (Ca2+ and Mg2+), which reduce H+ responsible for acidity, which, in turn, increases soil pH in the grassland soils. Tang et al. [82] reported a similar higher pH result (7.11 ± 0.27) of grassland than other land uses in the Horqin Sandy Land of China. On the other hand, plantation forests measured the lowest soil pH compared to other land use. The low pH scale in plantation forests may be ascribed to the effect of allelopathy associated with Eucalyptus trees and intense leaching of basic cations. Hence, chemicals released from the leaves and roots of Eucalyptus trees slow down the litter’s breakdown by affecting microbial activity, reducing the soil nutrients concentration, mainly base cations. In line with this result, Zerga [83] reported the case with allelopathy chemical released from Eucalyptus; Yerima and Van Ranst [84] also added reduced activity of litter decomposing microbial at low pH. Furthermore, low pH at the cultivated land in lower elevation may be associated to the intensive tillage frequency which stimulates leaching of base cations, the use of chemical fertilizers, and weedicides which contain high amounts of cations which neutralize the negative charges and due to the depletion of basic cations through crop harvesting. Supporting our result, Mohammed et al. [85], Mohammed [86], and Yihenew et al. [3] reported that the lowest value of pH under the cultivated land is due to the depletion of basic cations in crop harvest and leached streams in runoff generated from accelerated erosion. Similarly, Bore and Bedadi [65] also found the lower soil reaction under the cultivated land compared to the adjacent forest and grazing land in southern Ethiopia.

The lower and upper elevations recorded the highest mean soil pH value in the study watershed, which might be attributed to the accumulation of bases due to partial exclosures management practice at both landscape positions (Table 6). However, the lowest soil pH at the middle elevation is due to soil erosion, intensive cultivation, and free open grazing pressure. The overall pH value of the studied watershed is almost within the neutral range category, according to EthioSIS [56], for critical levels used for classifying soil fertility analysis for Ethiopian soils. The criteria used by the CASCAPE (CASCAPE: Capacity building for Scaling up of Evidence-based best Practices for increased Agricultural Production in Ethiopia) project [46], which is very detailed in comparison, show that the overall soil pH range of the study watershed is moderately acidic to neutral.

Soil EC is affected by cropping, irrigation, land use, fertilizer, manure, and compost. However, in this study, soil EC was found salt free according to EthioSIS [56] ratings and ranged from 0.16–0.47, 0.12–0.28, and 0.08–0.36 mS cm−1 in lower, middle, and upper topographies, respectively. The highest EC among different land use types were detected under GL at all topographies, followed by PF of both upstream positions (Table 6), whereas the lowest EC was obtained under BL of the upstream, followed by CL of the downstream watershed (Table 6). Contrary to our result, Kizilkaya and Dengiz [87] reported that changing a forest land to CL increased EC values due to high application rates of chemical fertilizers. Since EC represents the soil soluble salt component, the use of basic chemical fertilizers such as ammonium phosphate and urea under farmland can lead to higher EC values above normal (EC > 0.15 mS cm−1), but this did not happen in our study area may be due to runoff and leaching lose. Thus, farmers must avoid complete reliance on chemical inputs and should shift to rely more on organic fertilization to keep EC < 0.15 in their soils. Since the obtained EC value is less than one, the possibility of forming a salinity problem in the study watershed is implausible.

3.3.2. Soil Organic Matter (SOM), Soil Organic Carbon (SOC), Total N (TN), and C : N Ratio

The SOM, SOC, TN, and C : N ratio differed across altitudes and land use types and decreased as elevation increases in Ayiba watershed soils (Table 7). Based on the ratings of EthioSIS [56], the SOM of the study, watershed soil was found optimum at the lower elevation and low at middle and upper elevations, which might be ascribed to temperature decreases as elevation increases, which disfavor plants to survive at a higher altitude; hence, the accumulation of organic matter decreases. There was less SOM by about 43.98% and 48.19% across landscape positions, assuming the lower altitude as the base for comparison. Similarly, Shazia et al. [88] had reported high organic matter (2.89 ± 0.48 at 1660 m.a.s.l) at a lower elevation than at higher elevation (1.82 ± 0.57 at 2133 m.a.s.l.) in Sra Ghurgai (Takatu mountain range) Quetta, Balochistan, and many other studies also discovered the same trend related to elevation with organic carbon accumulation [89, 90].

Moreover, GL recorded the highest SOC (2.46%) in the lower landscape position, and BL recorded the lowest SOC (0.56%) in the upper landscape position (Table 7). Based on the ratings of Hazelton and Murphy [76] (<1.0, 1.0–1.71, 1.72–3.0, 3.1–4.29, and >4.3 leveled as very low, low, medium, high, and very high, respectively), the SOC of the study watershed soil was found from low to medium level across altitudinal gradients and very low to low level among land use types (Table 7). The higher SOM and SOC under GL (at lower and upper elevation) in this study are attributable to the partial exclosures practice, reduced grazing pressure and allowed regeneration of native vegetation, and increased SOC in the topsoil. Particularly, deforestation of native forests for crop production to feed the increasing population in mountainous areas dramatically increases the soil’s vulnerability for rainfall-driven erosion. Following the escalating erosion process, the organic humus and fine particles are transported down the slope from the upper streams with runoff water and accumulated at the lower stream. Therefore, this destructive ecosystem degradation increases SOM and SOC at a lower landscape position.

Recently, the reduction of soil quality and crop productivity has become severe in some Ethiopian highland regions, indicating the degradation of SOM and SOC, including the study area. Moreover, land use changes and subsequent tillage practices resulted in significant decreases in organic matter. Thus, cultivated soils generally have low organic matter content than native ecosystems since cultivation increases the soil’s aeration, enhancing SOM decomposition. Besides, complete crop harvest reduces the SOM in CL. The finding of this research work is in accordance with the work of Elias [57] and Gebreselassie et al. [61], who reported a higher mean value of SOC on a lower slope than the upper. Many other studies conducted in a mountainous area also reported that the cultivation of steep and fragile soils, limited recycling of dung and crop residues, deforestation, and overgrazing [91] are driving agents. Elias [4] and van Beek et al. [92] also reported poor soil management, soil erosion, leaching of nutrients, and organic matter depletion as driving factors for poor soil productivity in the Ethiopian highlands. Moreover, other than Ethiopia, Nsalambi [93] also reported similar results revealing total carbon concentrations, and its fractions were all higher at waterway (depression) and lowest at the summit position at Busby forest in central Missouri. Jendoubi et al. [94] also reported a marked decline in SOC with increased slopes under different land uses and topography in the Mediterranean landscape. Majule [40] also reported that cleaning the natural forest in most upper zones for cultivation has contributed significantly to SOC and N’s reduced soil level.

Given the independent effect of land use, we found the overall means SOM and SOC content was decreasing as GL > PF > CL > BL (Table 7) in this study. Explicitly, the estimated mean SOM and SOC were found higher in GL in all landscape positions except in the middle ones, where the PF was found high. The barren land at upper and middle and CL at lower topographies recorded low SOM and SOC in this study. Our finding agrees with the work of Girmay and Singh [10], who reported higher SOC in exclosures, GL, and plantation areas than CL in Mileba and Gum Slassie areas, northern Ethiopia. Delelegn et al. [95] also reported higher SOC in GL than cropland in Ambo Ber district, Ethiopia. SOC was found low compared to croplands, forestlands, and grazing lands [96]. According to Karltun et al. [56], SOC’s optimum baseline limit is 2%. SOC content in all land use types across topography except GL (2.46%) at a lower elevation was found below the baseline limit. Besides, Gurmu [97] highlighted that lower SOC could result in poor aggregate stability and, thus, aggravate soil degradation. Therefore, a sustained level of SOC is critically important for ecosystem health.

The highest TN content was observed in the GL (0.23%) at low altitude, followed by GL (0.20%) at upper and PF (0.20%) at middle altitudes. The average values of TN in the study watershed ranged from the highest 0.2% in the PF lands to the lowest 0.10% in BL among land use types and from the highest 0.20% at a low altitude to the lowest 0.15% at middle altitude across topography (Table 7). The mean TN content was inconsistent across altitude. The soils are found low to optimum in soil TN as per the EthioSIS ratings [56] with average mean values of 0.17%. The very low value of soil TN observed under BL and CL may be attributed to deforestation, soil erosion, continuous cultivation, plant uptake, and volatilization. Correspondingly, the low carbon input obtained from the subsistence agricultural crop production system could not compensate N losses by OM mineralization, leaching, and denitrification. Other studies similarly reported that deforestation and landscape fragmentation are the two most important factors affecting natural landscape settings, reducing the ecosystem services such as SOC and TN regulation accumulation [5, 98, 99]. On the other hand, TN’s relatively high mean value under PF lands could be attributed to the good microclimate conditions, which might have moderated the soil temperature and decreased TN loss by volatilization because the vegetation cover enriched the SOM. According to Walworth [100], the OM in the soil contributes roughly 95% for soil TN contents. The result is inconsistent with Ethiopia’s observations [67, 101]. There was an increase in soil TN contents of rehabilitated degraded bare land soil.

Carbon-to-Nitrogen (C : N) ratio is important because it directly impacts residue decomposition and nitrogen cycling in our soils with an ideal microbial diet of 24 : 1 (https://advancecovercrops.com/resources-advanced-cover-crops/carbon-nitrogen-ratio/).The C : N ratio indicates the dynamic balance of SOC and TN contents in a substance [79]. It is a good indicator of the degree of decomposition and mineralization rate of SOC and the cycle of carbon and nitrogen in the ecological system and good indicator of the quality of the organic matter held in the soil [102, 103]. The C : N ratios of the study watershed soils were significantly (< 0.001) affected by land use types and altitude (Table 2). The mean C : N ratio of the soils also varied between 6.38 and 8.61 among land use types and between 6.16 and 9.86 across altitude (Table 7). There are six to eight carbon units among land use types and six to nine units of carbon across altitude for each nitrogen unit in the substance. In this study, the C : N ratio showed a decreasing trend as altitude increases. GL at low altitude recorded the highest C : N ratio, and CL at upper altitude recorded the lowest C : N ratio. The present finding is in agreement with the work of Gebreslassie and Ayanna [104], who reported the highest values of the C : N ratio under grazing land in northwestern Ethiopian soils. Researchers reported that the C : N ratio is affected by land use due to competing interests of space for different purposes such as crop production, livestock, and settlement [99, 102]. Severe depletion of organic matter and continuous mining of soil N and S in the Ethiopian highlands is accredited to the practice of complete removal of crop residues as cattle feed and dung burning for household energy, as well as low levels of N fertilizer application (46 kg N ha−1) and no application of S fertilizer until recently [105]. Thus, the protection and management of forest and grassland fragments within the human-modified landscapes and supplying energy choice for the rural community are among the alternative options to keep the balance.

3.3.3. Available Phosphorus (Pav) and Exchangeable Potassium (Kexch)

ANOVA revealed that the Pav values significantly differed among land use ( < 0.05), altitudes ( < 0.01), and their interactions ( < 0.05) (Table 2), suggesting that land use types and the topographic factor of altitude exerted influence on Pav concentration. The mean Pav content of the study watershed soils ranges from 19.17 mg kg−1 to 40.20 mg kg−1 (Figure 6). The contents of Pav in the BL and CL at upper altitude appeared to be significantly higher than other land use types in the watershed. Accordingly, comparatively highest (40.2 mg kg−1) and the lowest (19.7 mg kg−1) Pav contents were recorded under the BL and PF, respectively. The data also revealed that the overall mean Pav was higher (35.98 mg kg−1) in the upper altitude soils than in the middle and lower ones. Based on Pav’s overall mean values, a decline of 28.08%, 16.54%, and 14.62% in PF, GL, and CL is noted compared to BL, respectively, and 25.74% and 18.37% in middle and lower topography compared to the upper topography, respectively. The low P fixation could explain the higher Pav content in the upper topography dominated by Leptosols of young soil sources that are historically less cultivated [106] and may be due to the application of the diammonium phosphate (DAP) fertilizer on cultivated lands. Similarly, Guo et al. [107] reported lower mean contents of Pav in paddy fields (lower and flat locations) than in vegetable lands (higher and steep places) (< 0.05) in southwestern China. This result is also in agreement with the findings of Moges et al. [4] and Elias [108], who specified that higher mean Pav in farmland compared to other land uses may be related to the application history of DAP fertilizer application. Despite substantial DAP application, the low to moderate phosphate levels, particularly in the intensive cereal system, suggest phosphate fixation [105]. Significant spatial variation of soil total phosphorus concentrations among land use types was also reported by Cheng et al. [109] in the Yingwugou watershed of the Dan River, China. According to the critical level adopted by EthioSIS [56], soil Pav of the Ayiba watershed is found in low to optimum status. The post hoc multiple comparison test of Pav among land use types also confirmed that BL had higher Pav content than GL and PF (Table 5).

Exchangeable K (Kexch) content was not significantly affected by land use, altitude, and the interaction of both (Tables 2). The mean Kexch obtained in all sampled land use types across the topography was classified as low for the study watershed according to the rating of EthioSIS [56].

3.4. Correlation Analysis

The correlation values showed an apparent association of the parameters with each other. The correlation analyses between SOM with other soil parameters show that the existence of positive and strong association with TP (r = 0.89, < 0.001), with soil MC and EC (r = 0.83, < 0.001), with TN (r = 0.75, < 0.001), and with clay content (r = 0.68, < 0.001). Pearson coefficient correlation analysis also shows that BD is inversely correlated with SOM (r = −0.89,  < 0.001), TN (r = −0.71,  < 0.001), and clay content (r = −0.68,  < 0.001) (Figure 7) but had a direct relationship with sand content and Pav (r = 0.84, < 0.01, and r = 0.56, < 0.05, respectively). The significant negative associations of sand and BD with SOC and SOM showed that conditions that favor increasing sand fractions and soil compaction (high BD) reduce SOC and SOM status. In turn, these processes reduce soil fertility levels that have a direct relationship with food security, indicating that a change in BD is associated with changes in SOM, clay, and sand content.

4. Conclusions

Understanding soil properties associated with land use and topographical attributes is vital for modeling soil-landscape relationships and establishing sustainable on-field management practices. Increased nutrient mining, soil erosion, and limited nutrient management have led to declining soil quality and reduced productivity in many parts of Ethiopia, which is more pronounced in highly rugged agrarian landscapes where land use influences and topographic variation dominate. In this regard, knowledge of soil properties associated with land use and altitudinal attributes is vital for establishing sustainable site-specific management practices. Therefore, this paper presents a study on the land use type and altitudinal gradient on selected soil physicochemical parameters in the Ayiba watershed, northern Ethiopia.

The variance results analysis shows that the combined effects of land use and variation in the altitudinal gradient caused a significant change in BD, TP, pH-water (1 : 2.5), EC, SOM, SOC, TN, and Pav in the soil of the Ayiba watershed. The result also revealed that BL and CL’s soils had the lowest SOM, SOC, and TN contents. On the other hand, soils of GL and PF had higher values of SOM, SOC, and TN contents. Thus, converting the degraded barren lands into protected lands/exclosures can increase these important soil quality parameters, but the reverse has a negative impact. Especially, expansion of cultivation into steep slopes and marginal areas causes severe loss of soil quality. Moreover, most measured soil parameters in the watershed are higher in the lower landscape position and relatively low in the upper landscape position where soil erosion is maximum due to its rugged nature. Compared to the critical values set, the different land use types and topographic positions’ measured soil physicochemical levels were low to optimum status. Given this, the conversion of natural forest or grassland to cultivation land and degraded barren land subjected to degradation of soil physicochemical properties, thereby ruining soil quality. Therefore, identifying appropriate land use and management practices is of utmost importance to keep soil health sustainably in the study watershed.

In general, the variation in soil physicochemical properties in the Ayiba watershed is not only related to the conversion of land use and topographic variation but also related to other factors including vegetation removal, complete removal of crop residue, and overgrazing, and the decline of fallowing is obvious, leading to soil degradation and loss of productivity. Additionally, livestock dung is baked into a dung cake and burned as household fuel, denying the soil an important nutrient and organic matter source. The most dangerous cropping practices that should change into other appropriate land use systems are the marginal and steep slope cultivated areas’. Therefore, to reverse soil degradation, further research on soil quality indicators on landscape-based should be evaluated and integrated into land use planning for sustainable land resource management practice in the Ayiba watershed.

Data Availability

Data used in this manuscript are available upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest or personal relationships that could have influenced the work reported in this paper.

Authors’ Contributions

WS was involved in conceptualization, methodology, formal analysis, investigation, data curation, writing the original draft, and visualization. EE was involved in conceptualization, methodology, supervision, fund acquisition, writing the review, and editing. GG was involved in conceptualization, methodology, supervision, fund acquisition, data curation, writing the review, and editing.

Acknowledgments

The authors thank the Mekelle University CASCAPE Project (funding for Ph.D. study) for funding this study. The authors are also very grateful to all farmers in the study area who allowed them to conduct their research on their vicinity.