Impact of Grazing Intensity on Soil Properties in Teltele Rangeland, Ethiopia

Grazing intensity (GI) is a major determining factor that controls the functioning of rangelands and the overall nutrient cycle. The Teltele rangeland is used for communal grazing area by the local pastorals; however, to date, there is no documented study data about the impact of GI. The objective of this study was to evaluate the impacts of grazing intensity on selected soil properties in the Teltele rangeland, Ethiopia. Soil samples were collected from different GI sites using different elevation gradient and soil depth from both open grazing and bush-encroached grazing land sand-assessed soil properties. Grazing intensity, elevation, and soil depth significantly (p < 0.05) affected both soils’ physical and chemical properties but rangeland types had no significant effect. The correlation analysis of soil characteristics with the principal component analysis axes showed significant variation. The highly weighted and correlated properties under principal component 1 (PC1) were electrical conductivity, organic carbon, total nitrogen, available phosphorus, and potassium, and under principal component 2, sand and bulk density with equal loaded value (r = −0.998), clay and silt, with silt (0.962) a more loaded one. Soil pH (0.743) demonstrated a significant (p < 0.05) positive correlation with sodium (−0.960) at PC1 (r = 0.610). Based on our results, we recommend further model-based studies on spatial–temporal change of soil properties due to impact of grazing intensity, combined with GIS and remote sensing data to be developed for sustainable rangeland management.

Introduction

Rangelands are lands on which the indigenous vegetation is predominantly grasses, grass-like plants, forbs, or shrubs and is managed as a natural ecosystem (Raj, 2005). Arid and semi-arid rangelands are heterogeneous in space and time because of variation in biotic and abiotic factors related to vegetation and soil properties and provide multiple ecosystem functions and services (Wang et al., 2016; Yang et al., 2016). Rangeland heterogeneity shapes vegetation structure and productivity (IPCC, 2013; Yigini and Panagos, 2016; Ademe et al., 2017). Variability in soil properties is a major main cause of rangeland heterogeneity (Ayalew, 2011). The major properties include soil textural, electrical conductivity (EC), organic matter (OM), and soil pH (Liu et al., 2011b; Abdalla et al., 2018). The primary use of the Teltele rangelands of Ethiopia is for livestock grazing (Derner et al., 2006). These rangelands are almost entirely occupied by a pastoral population using a system of communal resources for livestock production (Solomon et al., 2007). Grazing intensity is a major determining factor controlling rangeland functioning and the overall nutrient cycle (Hafner et al., 2012). Intensive livestock production and grazing gradually modify the soil characteristics, in particular organic carbon (OC), OM, EC, total nitrogen (TN), available phosphorus(P), exchangeable potassium (K), sodium (Na), texture, bulk density (BD), and pH (Pellant et al., 2000; Dessalegn et al., 2015; Zhou et al., 2017). Overgrazing can also cause soil erosion by reducing rangeland productivity and vegetation cover and in the long term, results in loss of environmental services, siltation of dams and river beds, reduction of groundwater, and social-community losses due to malnutrition and poverty (Adimassu et al., 2017). Furthermore, the removal of palatable species due to overgrazing suppresses their growth and facilitates the rapid encroachment of less desirable invasive species, mostly bush and shrubs plants species (Lin et al., 2010; Chen et al., 2015; Hailu et al., 2020).

Previous studies indicated that overgrazing increases soil heterogeneity (Su et al., 2006), while others reported that soil heterogeneity and vegetation diversity decrease with an increasing grazing intensity (Zhou et al., 2010; Zhao et al., 2011; Zhou et al., 2017). However, the majority of studies indicated that continuous and significant grazing intensities are generally accepted as having negative effects on OC (Piñeiro et al., 2010). Similarly, the effects of grazing on the spatial heterogeneity of grassland ecosystems related to soil properties have been inconsistent and need to be clarified in the Borana rangelands of Ethiopia. Evaluating dynamics of soil properties through grazing intensity (GI) requires clear and measurable data using comparable spatial methods at the study site. Therefore, understanding of soil properties is essential for rangeland management because such properties are among the primary factors that determine the forage production potential of an area in a particular climate (Hardy and Mentis, 1986; Aynekulu et al., 2017; Zhang et al., 2018). For centuries, the Teltele rangeland was used for communal grazing area by the local pastoralists, however, to date, there is no documented study data about the impact of GI on the soil properties in the rangeland area. This becomes one of the major gaps for sustainable rangeland management through balancing grazing capacity and maintaining rangeland productivity and livestock performance. Therefore, the objective of this study was to evaluate the soil properties in relation to the GI across different levels of altitude and grazing areas in Teltele rangeland, Ethiopia. So, this study aimed to address the following basic questions that can be used for effective implementation of management strategies and fill the knowledge gap mentioned above: 1) Is the significant difference observed in the soil properties due to variation of GI? 2) Does variation of grazing land type (GLT) had an impact on soil structure? 3) What is the interaction impact of GI with elevation (E) and soil depth (SD) on the soil properties? We hypothesized that 1) GI strongly affected soil properties, 2) GI had a similar impact on the soil properties both at the open grazing site and bush-encroached grazing site, and 3) interaction impact of GI, E, and SD is significant on the Teltele rangeland.

Materials and Methods

Study Area

This study was conducted from January–December 2019 in Teltele district, Borana zone, Ethiopia (Figure 1), which covers an area of 15,430 km² of which 68% (10,492 km²) is rangeland (Billi et al., 2015). The Teltele rangeland is 666 km south of Addis Ababa. The area is situated approximately between 4° 56′ 23″ and 5° 49′ 21″ N latitude and 37° 41′ 51″ and 38° 39′ 37″ E longitude. Mean elevation is 500–1,500 m with a maximum of 2059 m. The mean annual temperature ranges from 28–33°C with little seasonal variation. The rainfall is distributed with 60% from March–May and 27% from September–November and with high temporal and spatial fluctuations (Dalle et al., 2015). Potential evapotranspiration is 700–3,000 mm (Billi et al., 2015). The soil in the study area includes 53% red sandy loam soil, 30% black clay, and volcanic light-colored silt clay and 17% silt, and the vegetation mainly dominated by encroaching woody species, and those that frequently thinned out, including Senegalia mellifera, Vachellia reficiens, and Vachellia oerfota (Coppock and scarnecchia, 1994; Gemdeo et al., 2005). The 2017 national census data reported a total population of 100,501 in this district, 51,670 men and 48,831 women. Cattle, goats, sheep, camels, mules, donkeys, and horses are the main livestock species.

FIGURE 1

FIGURE 1. Location map of the study area.

Grazing Site Selection

A reconnaissance survey and discussion were conducted with local pastoralists and district Pastoral Development Officers on grazing intensity issues. The sampling site was selected both from open (free from any bush encroachment) and from bush-encroached grazing site since both grazing land types were available in the study area with different GI. In each grazing land type, three grazing treatments were categorized based on GI. Grazing intensity data were collected using the same procedures described by Fenetahun et al. (2020,2021), the same authors at the same study site. Based on the discussion and survey data, grazing sites were chosen based on similar, uniform, and same soil series. All sites were on the northeast part of Teltele woreda and have laid in similar slope and elevation range, and all sampling sites had been grazed by livestock for several decades up to date and almost have the same seasonal and environmental features. Cattle, goats, and sheep are among the dominant livestock. The grazing sites were used for yearly round, seasonal, and some are already fenced by the government for conservation and rehabilitation purposes in order to use it when the harsh environment like drought will happen. The status of pasture and rangeland condition of grazing site was used to estimate the level of GI (Morteza et al., 2012). We selected a site with two treatments: a non-grazing (NG) (as a control) and a grazing site (moderately grazing and overgrazing) that was considered to see the effect of grazing intensity based on grazing intensity gradient (Fenetahun et al., 2021). The rate of GI was described as follows: non-grazed (NG) (livestock have been excluded from the pasture by fence and the ground was almost completely covered by vegetation) and moderately grazed (MG) (pasture has been used for grazing in regular rotational basis, that is, used during non-dry seasons but not in the rainy season and vegetation covers almost 50–55% and overgrazed (OG) (pasture is used for grazing constantly throughout the year and totally grazed and undergoes degradation, and vegetation cover was in most cases less than 15%) for the last 1.5 years, and also, GI was divided into MG and OG based on the current carrying capacity potential (Fenetahun et al., 2020; Fenetahun et al., 2021). The treatments of sample collection involved at NG (∼0 ha AU⁻¹Y⁻¹), MG (6 ha AU⁻¹ Y⁻¹), and OG (12 ha AU⁻¹ Y⁻¹ and above) grazing area based on the current carrying capacity of rangeland was calculated by Fenetahun et al. (2020) and physical field observation. Also, we have selected the sampling site which has almost similar rainfall pattern and temperature in order to reduce the climate difference effect on soil composition.

Soil Sampling and Analyses

For soil sample collection, we applied the judgment sampling method (USEPA, 2002) to locate sampling sites for both open-grazing land (OGL) and bush encouraged (BE), landscape E lower (LE) (700–1,000 m), medium (ME) (1,001–1,300 m), and higher (HE) (1,301–1,600 m), and at soil depths (SD) of 0–10 cm and 10–20 cm using an auger at all GI because most grass roots are found within this layer (Mekuria et al., 2018; Zhu et al., 2015). Then, established a 5-km transect at each site, three main plots (50 m × 50 m), the western plot having a GI of NG, the middle plot a GI of MG, and the eastern plot a GI of OG, were marked and had a 1 km interval between each GI plots, and in each marked plot, five quadrats of (1 m × 1 m) were placed at 5 m buffer zone for each sampling quadrat that was used for soil sample collection (Figure 2). Soil sampling was done both during the dry season (end of February 2019) and the rainy season (end of May 2019) along each of the grazing sites in order to overcome the season variation effect and we took the mean value. Thus, a total of 108 quadrat samples were collected (2 grazing land types × 3 landscape E × 3 GI plots × 2 SD) × 3 replications. The samples were mixed at the point of sampling and 0.5 kg sub-samples from each sampling point were taken in the laboratory in a plastic bag and were oven-dried at 105°C for 24 h in order to avoid delay. Samples for the BD estimation were collected two days after a rainfall event, a 5 cm × 5 cm soil corer of volume 98.13 cm³ was used, making sure not to disturb soil aggregation (Masebo et al., 2014; Hishe et al., 2017). Samples were analyzed in the soil laboratory at Yabello Pastoral and Dryland Agricultural Soil Research Center. After drying, soil samples were crushed to pass a 2-mm stainless steel sieve to remove foreign bodies (Hishe et al., 2017). Analysis was performed for OC, TN, P, K, Na, pH, EC contents, and particle size distribution (clay, silt, and sand) following standard procedures is described (Table 1) below.

FIGURE 2

FIGURE 2. Sampling lay out.

TABLE 1

TABLE 1. Standard procedures and methods used to analyses soil properties.

Data Analysis

The effect of GLT (OGL and BE), GI, and E impact on selected soil properties was analyzed by using analysis of variance (ANOVA) effect in SAS version 9.1.3 (Statistical Analysis System) with different statistical packages. Some of the properties possessed extreme outliers and did not meet the normality assumptions. Then, an assessment of significant differences evaluation at p < 0.05 was used to analyze the impact of treatment using the LSMEANS procedure (Yang and Luo, 2011). A Spearman rank correlation analysis and matrix were carried out to investigate the impact of soil factors resulting from the different GI, and a full set of soil properties data across GI, E, SD, and GLT were subjected to principal component analysis (PCA) to evaluate the impact of GI on each soil properties and also in order to evaluate the bonding characteristics within each soil properties. The criterion used for selecting the optimal subset of the main component (PC) is to select a subset with eigenvalues greater than 1.

Results

Physical Properties

From our result, we can understand that both levels of GI, E, and SD, had significant (p < 0.05) effects whereas GLT (OGL and BE) had no significant (p > 0.05) effect on all of the soil physical properties. Sand soil content was highest on the OG level of grazing than the MG and NG, at both SD and E, particularly at HE grazing position and at the 0–10 cm of SD and lowest on the NG level of grazing, at LE grazing position and at the 10–20 cm of SD in both GLT. Clay and silt soil content were highest on NG level of grazing than MG and OG at both SD and E, particularly at LE grazing position and at 10–20 cm of SD and lowest on OG level of grazing, at HE grazing position and at 0–10 cm of SD (Table 2). Also, the interaction (X) effect of both GI, E, and SD, across GLT had significant effects on all of the physical properties of the soil, both within and across the different grazing sites (Table 3 and Figure 3). The highest sand soil particle distribution was recorded at OG X BE X HE X 0–10 cm depth and the lowest was recorded at NG X BE X LE X 10–20 cm depth. The highest clay and silt soil particle distribution were recorded at NG X BE X LE X 10–20 cm depth and the lowest clay and silt were recorded at OG X BE X HE X 0–10 cm depth. The highest concentration of BD was observed at 0–10 cm depth in all GI and E. The highest BD was recorded at OG, in the BE site, at HE of 0–10 cm depth and lowest was recorded at NG, in the BE site at LE of 10–20 cm depth. This indicated that the interaction effect of GLT had significant effects on all of the physical properties of the soil.s

TABLE 2

TABLE 2. Mean (±SE) value of soil physical particle size distribution under different grazing land type, elevation, grazing intensity, and soil depths.

TABLE 3

TABLE 3. Mean (±SE) value of soil physical particle size distribution under the interaction effects of grazing intensity with grazing land type, elevation, and soil depths.

FIGURE 3

FIGURE 3. Mean ± SE particle size distribution across grazing intensity. Bars with different letters are significantly different for each grazing intensity site, PSD, particle size distribution; NG, non-grazing; MG, moderately grazing; OG, overgrazing.

Chemical Properties

Results indicated that both levels of GI, E, and SD, exhibited a significant (p < 0.05) effect but GLT had no significant (p > 0.05) effect on all of the soil chemical properties. The highest EC, OC, TN, P, and K contents were recorded at NG than the MG and OG level of grazing at both SD and E, particularly at HE grazing position and at the 10–20 cm of SD and lowest at OG level of grazing at both SD and E, particularly at LE grazing position and at the 0–10 cm of SD in both GLT. The highest pH and Na values were recorded at OG than NG and MG level of grazing at both SD and E, particularly at LE and at the 0–10 cm of SD and lowest at NG level of grazing, especially at HE and at the 10–20 cm of SD (Table 4). Interaction (X) effects of both GI, E, and SD, across GLT had significant effects on all of the chemical properties of the soil, both within and across the different grazing sites (Table 5 and Figure 4). The highest EC, OC, TN, P, and K contents were recorded at NG X OGL X HE X 10–20 cm depth and the lowest was recorded at OG X BE X LE X 0–10 cm depth. The highest pH and Na values were recorded at OG X BE X LE X 0–10 cm depth and the lowest was recorded at NG X OGL X HE X 10–20 cm depth. Further, the interaction effect of bush encroachment also had affected the soil chemical properties distribution in the grazing rangeland.

TABLE 4

TABLE 4. Mean (±SE) value of soil chemical properties distribution under different grazing land type, elevation, grazing intensity, and soil depths.

TABLE 5

TABLE 5. Mean (±SE) value of soil chemical properties distribution under the interaction effects of grazing intensity with grazing land type, elevation and soil depths.

FIGURE 4

FIGURE 4. Mean ± SE soil chemical properties abundance across grazing intensity. Bars with different letters are significantly different for each grazing intensity site. SCP, soil chemical properties; GI, grazing intensity; NG, non-grazing; MG, moderately grazing; OG, overgrazing.

Correlation Analysis of Soil Properties With Grazing Intensity

The correlation analysis of the regression lines describes the relationship between the soil properties and the GI. In relation to soil physical properties, there is a positive correlation between sand and BD contents with the GI but a negative correlation with clay and silt contents with the rate of GI (Figure 5). In the case of soil chemical properties, GI showed a negative correlation with the soil EC, OC, TN, P, and K contents and a positive correlation with soil pH and Na contents (Figure 6).

FIGURE 5

FIGURE 5. Correlations between soil physical properties and grazing intensity. GI, grazing intensity; NG, non-grazing; MG, moderately grazing; OG, overgrazing.

FIGURE 6

FIGURE 6. Correlations between soil chemical properties and grazing intensity (GI, grazing intensity; NG, non-grazing; MG, moderately grazing; OG, overgrazing.

Correlation Matrix and Principal Component Analysis of Soil Particle

Observing the correlation matrix and PCA that were computed for each pair of soil properties at each GLT along with the interaction effect of GI, E, and SD, and EC, OC, TN, P, and K showed positive correlations with each other and with clay and silt soil contents, whereas, a negative correlation with Na, pH, and sand soil contents across both GLT, GI, E, and SD (Table 6 and Figure 7B). The location of soil properties under different regions of the PCA axes is based on the correlation coefficient between each variable. Since the main components are orthogonal, this defines a projection of the data on vector space spanned by the first two principal components. Thus, we used two PCs with eigenvalues >1 for our study (Table 6 and Figure 7A). Therefore, the location of each soil property in the PCA diagram is very significant and important. The highly weighted and positively correlated properties under PC1 were EC, OC, TN, P, and K and highly impacted or negatively correlated with Na and pH (Figure 7B). Under PC2, the highly weighted and positively correlated properties were sand and BD with a high negative correlation of the other weighted and positively correlated properties of clay and silt. Soil pH showed a strong positive correlation with Na at PC1 with a more loaded value of Na, and this result strongly agrees with the data reported by Kane (2015).

TABLE 6

TABLE 6. Correlation matrix within each pair of soil properties across the interaction effect of impacting factors.

FIGURE 7

FIGURE 7. (A) Scree plot for the different components considered for the principal component analysis with eigenvalues greater than one and (B) principal component analysis (PCA) of the overall data set. OC, organic carbon; TN, total nitrogen; Av. p, available phosphorus; Av. K, available potassium; Ex. Na, exchangble sodium; BD, bulk density; EC, electrical conductivity.

Discussion

Impact of Grazing Intensity on Soil Properties

Grazing Intensity Impact on Soil Physical Properties

The Teltele rangeland site was selected because it is one of the driest parts of the Borana region and therefore the pastoral communities in this region are the most vulnerable to rangeland degradation due to overgrazing. The GI significantly affected all the soil properties that we measured. The clay and silt soil contents distribution were decreased with increasing GI and decreased SD and decreased with increasing E, while the sand soil content distribution was increased with decreasing GI and increasing E and decreased with increasing SD (Stavi et al., 2008; Larreguy et al., 2014; Yun and Wesche, 2016; Lu et al., 2017; Zhang et al., 2018). In general, in this study, grazing was associated with higher BD, increase in the sand, decreased clay and silt, decreased soil moisture, decreased diversity and coverage of grass species, and results increasing in unpalatable plant species. The result showed a positive relationship between sand and BD soil contents with GI, that is, increased while GI increasing and vice versa. This indicates the direct linkage of BD with sand content and an inverse linkage with clay and silt content and results in line with data reported by Gashaw et al. (2017), Tufa et al. (2019). In addition, the major cause for increasing of BD on the grazing site is the presence of high sand soil content distribution and higher porous spaces compared with other soil particles (Adesodun et al., 2007; Wolka et al., 2011; Chaudhari et al., 2013; Adimassu et al., 2017; Zhou et al., 2017; Paulo et al., 2018). The soil structure, moisture content, OM composition, and GI determine the soil content distribution in the grazing site (Liu et al., 2011b; Ademe et al., 2017; Tufa et al., 2019). For instance, the NG site has a better OM and moisture content and in the NG grazing level, clay and silt soil contents were observed dominantly and mainly at LE position of 10–20 cm of SD. This is due to the upper part of 0–10 cm SD, easily exposed to any livestock trampling during grazing and facilitates wind and flood erosion. This speedup water runoff and loss of soil moisture and changes to bare land, resulting in clay and silt soil particle easily eroded from the site and dominated with sand soil content (Fei et al., 2010; Zhou et al., 2010; Morteza et al., 2012; Descalzi et al., 2018). The data reported by Alemayehu and Fisseha (2018), Demelash and Karl (2010) indicated that higher clay and lower sand content was recorded at non-conserving or OG sites without any influence of the difference in E in contrast to our current result. Consequently, the linear regression analysis showed that the GI had a positive relationship with sand and BD soil contents and a negative relationship with clay and silt contents (Ademe et al., 2017).

Grazing Intensity Impact on Soil Chemical Properties

The soil chemical properties are soil management properties dependent on the soil structure, air and water conductivity, and highly influenced by grazing management. As a result, grazing significantly affected the soil chemical properties concentration found at the grazing rangeland site. At the grazing site of LE, the density of livestock was higher than that of ME and HE and had a great effect on the soil chemical properties distribution of the grazing site, even under similar GI and GLT. The soil EC was high at the NG grazing level than that in the MG and OG. The reason for the highest EC observed at the NG site mainly at HE in the soil surface of 10–20 cm is due to the high availability of OC, TN, P, and K, resulting in high cation and anion concentration, and EC is the sum of two ions (Kidane, 2006; Yan et al., 2013; Zhang et al., 2018; Guo et al., 2019). The major reason for the distribution of OC, TN, P, and K was higher in the managed area or NG area is due to the availability of higher grass biomass which results in increased availability of soil nutrients during decomposition (Hamilton and Frank, 2001; Wang et al., 2018),but when the rangeland was exposed to continuous grazing and transformed into OG degraded area, the aboveground grass biomass declined and the formation of OM was affected. The rangeland site exposed to different erosion agents like wind and water resulted in a lower distribution of OC, TN, P, and K and lead to the distribution of Na and pH become higher (Zhang et al., 2009; Chaudhari et al., 2013; Guo et al., 2019). The other possible reason for the highest EC, OC, TN, P, and K values was recorded at soil surface of 10–20 cm depth compared to the soil surface of 0–10 cm was due to high distribution of clay and silt soil particle at soil surface of 10–20 cm than 0–10 cm (Demelash and Karl, 2010; Aytenew and Kibret, 2016; Ademe et al., 2017). Those soil chemical properties had an abundance of direct linkage between them. The pH and Na showed higher values at 0–10 cm of SD across all grazing level, especially at the OG site of HE, and this was due to the high distribution of sand soil particles and a positive relationship with pH and Na (Tufa et al., 2019; Yang et al., 2016; Wolka et al., 2011). As pH increases, the availability of certain major basic cations like Na, clay, and CEC are positively correlated because clay minerals provide the negative charge to attract the cations. The variation of GLT does not show a significant effect on the soil chemical properties except with slight variation. This means the values of EC, OC, TN, P, and K showed slightly higher value at OGL compared to BE grazing site, whereas pH and Na values showed higher at BE at OGL grazing site and our result is in agreement with the data reported by Mulder et al. (2015). Further, the linear regression analysis showed that a significant negative relationship with the soil EC, OC, TN, P, and K and a positive relationship with Na and pH as GI increased, and this result supported with the data reported by Kate (2019), Hao and He (2019). Overall, grazing has shown a decreasing effect on the abundance of soil organic matter and the water-holding capacity, which leads to an increase in soil BD, sandy soil particles, and Na and pH and decrease in EC, OC, TN, P, and K; therefore, our hypothesis was accepted. Managing rangeland grazing used to protect not only the soil contents of the rangeland but also the water availability, and this research will be used as a reference and initiative for further research and aware the pastoralist community through showing how over stoking currently impacts both the rangeland soil and productivity in the Teltele rangeland, and this was the implication of our current work.

Conclusion

In the Teltele rangeland, grazing intensity strongly influenced the soil properties of the grazing rangeland. The increase in the distribution of the apparent density of the soil was mainly due to the increase in the distribution of sand soil particles and decline of sand and silt soil particle. This caused speedup of water infiltration, changes in chemical properties, and fertility of the soil and is among the major impacts of overgrazing. It can be concluded that managing the level of grazing is an essential technique used for improving arid and semi-arid rangeland areas including Teltele. Managing the grazing period and balancing the number of livestock grazed on the grazing land help reduce the grazing density and restore the soil properties, through improved vegetation cover and reduced runoff and erosion. Variation of grazing land type had less impact on the soil properties as compared with grazing intensity, elevation, and soil depth difference effect. Further, studies are needed for better understanding of seasonal management of grazing intensity for a better improvement of the grazing land soil properties and enhanced general ecosystem function of rangeland. Timely reform and balancing of carrying capacity of livestock at a certain grazing area are needed for proper rehabilitation of rangeland through providing a recovery period and for the proper implementation of the limited number of livestock and reduce grazing intensity, introduction and application of appropriate laws that is formulated by both the local communities and government official that govern the way how to use of communal grazing areas. Based on our result, we also recommend that the influence of grazing intensity should be further studied by combining GIS, remote sensing, and NDVI data to see temporal and spatial changes due to the effect of grazing and model-based data that showed a change of soil properties is more reliable in the Teltele rangeland site.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

YF: data collection, writing up, gap assessment and design of the experiments. XX, WD, TF, and VN: editing, proofing, laboratory assistance, provide important advice as well as supervision of the whole work.

Funding

This study received financial support from CAS-TWAS fellowship program and African Great Green Wall Adaptation Technical Cooperation Research and Demonstration (Grant No. 2018YFE0106000), Science and Technology Partnership Program, Ministry of Science and Technology of China (Grant No. KY 201702010), Integration and application of appropriate technologies for desertification control in Africa (Grant No. SAJC2021080, and International cooperation and Exchanges NSFC (Grant No. 41861144020).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank the University of Chinese Academy of Science and CAS-TWAS fellowship program that provides funding and the PhD Scholarship for the first author. They also acknowledge the African Great Green Wall Adaptation Technical Cooperation Research and Demonstration (2018YFE0106000), Science and Technology Partnership Program, Ministry of Science and Technology of China (Grant Nos. KY 201702010), and International cooperation and Exchanges NSFC (Grant Nos. 41861144020) financial support to do this paper, also our great thanks go to the local community and stakeholder of the Teltele district for giving us the basic information that is still the challenge for them for our next research step. Once again, Yabello Pastoral and Dryland Agricultural Soil Research Center deserve many thanks for providing us the laboratory and other facilitates while conducting our experiment.

References

Abdalla, M., Hastings, A., Chadwick, D. R., Jones, D. L., Evans, C. D., Jones, M. B., et al. (2018). Critical Review of the Impacts of Grazing Intensity on Soil Organic Carbon Storage and Other Soil Quality Indicators in Extensively Managed Grasslands. Agric. Ecosyst. Environ. 253, 62–81. doi:10.1016/j.agee.2017.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Ademe, Y., Temesgen, K., Alemayehu, M., and Toyiba, S. (2017). Evaluation of the Effectiveness of Soil and Water Conservation Practices on Improving Selected Soil Properties in Wonago District, Southern Ethiopia. J. Soil Sci. Environ. Manage. 8, 70–79. doi:10.5897/jssem2016.0601

CrossRef Full Text | Google Scholar

Adesodun, J., Adeyemi, E., and Oyegoke, C. (2007). Distribution of Nutrient Elements within Water-Stable Aggregates of Two Tropical Agro-Ecological Soils under Different Land Uses. Soil Tillage Res. 92, 190–197. doi:10.1016/j.still.2006.03.003

CrossRef Full Text | Google Scholar

Adimassu, Z., Langan, S., Johnston, R., Mekuria, W., and Amede, T. (2017). Impacts of Soil and Water Conservation Practices on Crop Yield, Run-Off, Soil Loss and Nutrient Loss in Ethiopia: Review and Synthesis. Environ. Manage. 59, 87–101. doi:10.1007/s00267-016-0776-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Alemayehu, T., and Fisseha, G. (2018). Effects of Soil and Water Conservation Practices on Selected Soil Physico-Chemical Properties in Debre-Yakob Micro-watershed, Northwest Ethiopia. Eth J. Sci. Technol. 11, 29–38. doi:10.4314/ejst.v11i1.3

CrossRef Full Text | Google Scholar

Ayalew, A. (2011). Construction of Soil Conservation Structures for Improvement of Crops and Soil Productivity in the Southern Ethiopia. J. Environ. Earth Sci. 1, 21–29.

Google Scholar

Aynekulu, E., Mekuria, W., Tsegaye, D., Feyissa, K., Angassa, A., De leeuw, J., et al. (2017). Long-term Livestock Exclosure Did Not Affect Soil Carbon in Southern Ethiopian Rangelands. Geoderma 307, 1–7. doi:10.1016/j.geoderma.2017.07.030

CrossRef Full Text | Google Scholar

Aytenew, M., and Kibret, K. (2016). Assessment of Soil Fertility Status at Dawja Watershed in Enebse Sar Midir District, Northwestern Ethiopia. Ijpss 11, 1–13. doi:10.9734/ijpss/2016/21646

CrossRef Full Text | Google Scholar

Billi, P., Alemu, Y. T., and Ciampalini, R. (2015). Increased Frequency of Flash Floods in Dire Dawa, Ethiopia: Change in Rainfall Intensity or Human Impact? Nat. Hazards 76 (2), 1373–1394. doi:10.1007/s11069-014-1554-0

CrossRef Full Text | Google Scholar

Chaudhari, P. R., Ahire, D. V., Ahire, V. D., Chkravarty, M., and Maity, S. (2013). Soil Bulk Density as Related to Soil Texture, Organic Matter Content and Available Total Nutrients of Coimbatore Soil. Int. J. Sci. Res. Publ. 3, 1–8.

Google Scholar

Chen, W., Huang, D., Liu, N., Zhang, Y., Badgery, W. B., Wang, X., et al. (2015). Improved Grazing Management May Increase Soil Carbon Sequestration in Temperate Steppe. Scientific Rep. 5, 10892. doi:10.1038/srep10892

PubMed Abstract | CrossRef Full Text | Google Scholar

Coppock, D. L., and Scarnecchia, D.L., (1994). The Borana Plateau of Southern Ethiopia: Synthesis of Pastoral Research, Development and Change. J. Range Management. 48(6):566. doi:10.2307/4003073

CrossRef Full Text | Google Scholar

Dalle, G., Maass, B. L., and Isselstein, J. (2015). Rangeland Condition and Trend in the Semi-arid Borana Lowlands, Southern Oromia, Ethiopia. Afr. J. Range Forage Sci. 23 (1), 49–58. doi:10.2989/10220110609485886

CrossRef Full Text | Google Scholar

Demelash, M., and Karl, S. (2010). Assessment of Integrated Soil and Water Conservation Measures on Key Soil Properties in South Gonder, North-Western Highlands of Ethiopia. J. Soil Sci. Environ. Manage. 1, 164–176.

Google Scholar

Derner, J. D., Boutton, T. W., and Briske, D. D. (2006). Grazing and Ecosystem Carbon Storage in the North American Great Plains. Plant and Soil 280 (1-2), 77–90. doi:10.1007/s11104-005-2554-3

CrossRef Full Text | Google Scholar

Descalzi, C., Balocchi, O., López, I., Kemp, P., and Dörner, L. (2018). Different Soil Structure and Water Conditions Affect the Growing Response of Lolium Perenne L. And Bromus Valdivianus Phil. Growing Alone or in Mixture. J. Soil Sci. Plant Nutr. 18 (3), 617–635. doi:10.4067/s0718-95162018005001901

CrossRef Full Text | Google Scholar

Dessalegn, C. D., Christian, D. G., Assefa, D. Z., Tigist, Y. T., Menelik, G., Solomon, A., et al. (2015). Impact of Conservation Practices on Runoff and Soil Loss in the Sub-humid Ethiopian Highlands: the Debre Mawi Watershed. J. Hydrol. Hydromechanics. 63, 210–219. doi:10.1515/johh-2015-0021

CrossRef Full Text | Google Scholar

Fei, P., Tao, W., Xian, X., and Fang, Z. (2010). Soil and Plant Responses to Degradation of Alpine Grassland in Source Region of the Yellow River. Sci. Cold Arid Regions 2 (4), 0364–0370. doi:10.1016/j.proenv.2010.01.090

CrossRef Full Text | Google Scholar

Fenetahun, Y., Wang, Y. D., You, Y., and Xu, X. (2020). Dynamics of Forage and Land Cover Changes in Teltele District of Borana Rangelands, Southern Ethiopia: Using Geospatial and Field Survey Data. BMC Ecol. 20, 55. doi:10.1186/s12898-020-00320-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenetahun, Y., Yuan, Y., Xinwen, X., and Yongdong, W. (2021). Effects of Grazing Enclosures on Species Diversity, Phenology, Biomass, and Carrying Capacity in Borana Rangeland, Southern Ethiopia. Front. Ecol. Evol. 8, 623627. doi:10.3389/fevo.2020.623627

CrossRef Full Text | Google Scholar

Gashaw, T., Taffa, T., and Mekuria, A. (2017). Erosion Risk Assessment for Prioritization of Conservation Measures in Geleda Watershed, Blue Nile Basin, Ethiopia. Environ. Syst. Res. 6, 1–14. doi:10.1186/s40068-016-0078-x

CrossRef Full Text | Google Scholar

Gemedo, D., Maass, B. L., and Isselstein, J. (2005). Plant Communities and Their Species Diversity in the Semi-arid Rangelands of Borana Lowlands, Southern Oromia, Ethiopia. Community Ecol. 6 (2), 167–176. doi:10.1556/comec.6.2005.2.5

CrossRef Full Text | Google Scholar

Guo, X., Chen, L., Zheng, R., Zhang, K., Qiu, Y., and Yue, H. (2019). Differences in Soil Nitrogen Availability and Transformation in Relation to Land Use in the Napahai Wetland, Southwest China. J. Soil Sci. Plant Nutr. 19, 92–97. doi:10.1007/s42729-019-0013-0

CrossRef Full Text | Google Scholar

Hafner, S., Unteregelsbacher, S., Seeber, E., Lena, B., Xu, X., Li, X., et al. (2012). Effect of Grazing on Carbon Stocks and Assimilate Partitioning in a Tibetan Montane Pasture Revealed by 13CO2 Pulse Labeling. Glob. Change Biol. 18 (2), 528–538. doi:10.1111/j.1365-2486.2011.02557.x

CrossRef Full Text | Google Scholar

Hailu, T., Mekuria, A., Lulseged, T., Kindu, M., John, R., and Dawit, S. (2020). Effects of Sustainable Land Management Interventions on Selected Soil Properties in Geda Watershed, Central Highlands of Ethiopia. Ecol. Process. 9, 14. doi:10.1186/s13717-020-0216-2

CrossRef Full Text | Google Scholar

Haldar, A., and Sakar, D. (2005). Physical and Chemical Method in Soil Analysis: Fundamental Concepts of Analytical Chemistry and Instrumental Techniques. New Delhi: New Age International (P) Ltd. doi:10.1145/1065907.1066039

CrossRef Full Text

Hamilton, E. W., and Frank, D. A. (2001). Can Plants Stimulate Soil Microbes and Their Own Nutrient Supply? Evidence from a Grazing Tolerant Grass. Ecology 82, 2397–2402. doi:10.1890/0012-9658(2001)082[2397:cpssma]2.0.co;2

CrossRef Full Text | Google Scholar

Hao, Y., and He, Z. (2019). Effects of Grazing Patterns on Grassland Biomass and Soil Environments in China: A Meta-Analysis. PLoS ONE 14 (4), e0215223. doi:10.1371/journal.pone.0215223

PubMed Abstract | CrossRef Full Text | Google Scholar

Hardy, M. B., and Mentis, M. T. (1986). Grazing Dynamics in Sour Grassveld. S. Afr. J. Sci. 82, 566–572.

Google Scholar

Hishe, S., Lyimo, J., and Bewket, W. (2017). Soil and Water Conservation Effects on Soil Properties in the Middle Silluh Valley, Northern Ethiopia. Int. Soil Water Conservation Res. 5, 231–240. doi:10.1016/j.iswcr.2017.06.005

CrossRef Full Text | Google Scholar

IPCC (2013). Climate Change 2013-The Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press. doi:10.1017/CBO9781107415324

Kane, D. (2015). Carbon Sequestration Potential on Agricultural Lands: A Review of Current Science and Available practicesAssociation with: National Sustainable Agriculture Coalition Breakthrough Strategies and Solutions. Washington, DC, United States: LLC.

Kate, K. (2019). in Sustainable International Agriculture Specialization in Tropical Agriculture Faculty of Agricultural Sciences Institute of Grassland Science (Georg-August-Universität- Göttingen, Germany). Model-Based Assessment of Grazing Impact on Soil Organic Carbon Stocks and Dynamics of a Kenyan. Rangeland. M.Sc. Thesis. Nairobi-Kenya.

Kidane, G. (2006). Rangeland Potential, Quality and Restoration Strategies in North-Eastern Ethiopia: a Case Study Conducted in the Southern Afar Region. Dissertation Submitted in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy at the University of Stellenbosch.

Google Scholar

Larreguy, C., Carrera, A. L., and Bertiller, M. B. (2014). Effects of Long-Term Grazing Disturbance on the Belowground Storage of Organic Carbon in the Patagonian Monte, Argentina. J. Environ. Manage. 134, 47–55. doi:10.1016/j.jenvman.2013.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y., Hong, M., Han, G., Zhao, M., Bai, Y., and Chang, S. X. (2010). Grazing Intensity Affected Spatial Patterns of Vegetation and Soil Fertility in a Desert Steppe. Agric. Ecosyst. Environ. 138, 282–292. doi:10.1016/j.agee.2010.05.013

CrossRef Full Text | Google Scholar

Liu, T., Nan, Z., and Hou, F. (2011b). Grazing Intensity Effects on Soil Nitrogen Mineralization in Semi-arid Grassland on the Loess Plateau of Northern China. Nutr. Cycl Agroecosyst 91, 67–75. doi:10.1007/s10705-011-9445-1

CrossRef Full Text | Google Scholar

Lu, X., Kelsey, K. C., Yan, Y., Sun, J., Wang, X., Cheng, G., et al. (2017). Effects of Grazing on Ecosystem Structure and Function of Alpine Grasslands in Qinghai-Tibetan Plateau: a Synthesis. Ecosphere 8. doi:10.1002/ecs2.1656

CrossRef Full Text | Google Scholar

Masebo, N., Abdellkadir, A., and Mohammed, A. (2014). Evaluating the Effect of Agroforestry-Based Soil and Water Conservation Measures on Selected Soil Properties at Tembaro District, SNNPR, Ethiopia. Direct Res. J. Agric. Food Sci. 2, 141–146.

Google Scholar

Mekuria, W., Menale, W., Tadele, A., Asmare, W., Tesfaye, F., and Birru, Y. (2018). Restoration of Degraded Landscapes for Ecosystem Services in North-Western Ethiopia. Heliyon 4, e00764. doi:10.1016/j.heliyon.2018.e00764

PubMed Abstract | CrossRef Full Text | Google Scholar

Miles, N., and Farina, M. (2013). Soil Acidity and its Management in Crop Production, Available from http://www.grainsa.co. doi:10.4324/9781315024653

CrossRef Full Text | Google Scholar

Morteza, M., Mehdi, R., Hassan, A., Ali, E., Esmaeil, S., and Javad, M. (2012). Impact of Grazing on Chemical, Physical and Biological Properties of Soils in the Mountain Rangelands of Sahand, Iran. Rangeland J. 34, 297–303. doi:10.1071/RJ11087

CrossRef Full Text | Google Scholar

Motsara, N. R., and Roy, R. N. (2008). Guide to Laboratory Establishment for Plant Nutrient analysisFAO Fertilizer and Plant Nutrition Bulletin. United Nations: Food and agriculture organization.

Mulder, V. L., Lacoste, M., Martin, M. P., Richer‐de‐Forges, A., and Arrouays, D. (2015). Understanding Large‐extent Controls of Soil Organic Carbon Storage in Relation to Soil Depth and Soil‐landscape Systems. Glob. Biogeochem. Cycles 29 (8), 1210–1229. doi:10.1002/2015gb005178

CrossRef Full Text | Google Scholar

Paulo, R., Paulo, H. R., Luiz, F. M., Felipe, V. A., and Eduardo, S. M. (2018). Distribution of C and Inorganic Phosphorus Fractions in Different Aggregate Sizes under Forestry, Agroforestry System and Pasture. J. Soil Sci. Plant Nutr. 18 (2), 361–375. doi:10.4067/S0718-95162018005001202

CrossRef Full Text | Google Scholar

Pellant, M., Shaver, P., Pyke, D., and Herrick, J. (2000). Interpreting Indicators of Rangeland Health. available at: ftp://ftp.ftw.nrcs.usda.gov/pub/glti/IntIndRangeHealth.pdf.Version 3. Technical Reference 1734–1736,

Piñeiro, G., Paruelo, J. M., Oesterheld, M., and Jobbágy, E. G. (2010). Pathways of Grazing Effects on Soil Organic Carbon and Nitrogen. Rangeland Ecol. Management 63 (1), 109–119. doi:10.2111/08-255.1

CrossRef Full Text | Google Scholar

Raj, K. (2005). The Impact of Rangeland Condition and Trend to the Grazing Resources of a Semi-arid Environment in Kenya. J. Hum. Ecol. 17 (2), 143–147. doi:10.1080/09709274.2005.11905769

CrossRef Full Text | Google Scholar

Solomon, T., Snyman, H. A., and Smit, G. N. (2007). Rangeland Dynamics in Southern Ethiopia: (1) Botanical Composition of Grasses and Soil Characteristics in Relation to Land-Use and Distance from Water in Semi-arid Borana Rangelands. J. Environ. Manage. 85, 429–442. doi:10.1016/j.jenvman.2006.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Stavi, I., Ungar, E. D., Lavee, H., and Sarah, P. (2008). Grazing-induced Spatial Variability of Soil Bulk Density and Content of Moisture, Organic Carbon and Calcium Carbonate in a Semi-arid Rangeland. Catena 75, 288–296. doi:10.1016/j.catena.2008.07.007

CrossRef Full Text | Google Scholar

Su, Y. Z., Li, Y. L., and Zhao, H. L. (2006). Soil Properties and Their Spatial Pattern in a Degraded Sandy Grassland under Post-grazing Restoration, Inner Mongolia, Northern China. Biogeochemistry 79, 297–314. doi:10.1007/s10533-005-5273-1

CrossRef Full Text | Google Scholar

Tufa, M., Melese, A., and Tena, W. (2019). Effects of Land Use Types on Selected Soil Physical and Chemical Properties: The Case of Kuyu District, Ethiopia. Eur. J. Soil Sci. 8, 94–109. doi:10.18393/ejss.510744

CrossRef Full Text | Google Scholar

United States Environmental Protection Agency (2002). Guidance on Choosing a Sampling Design for Environmental Data Collection. Washington, DC: Office of Environmental Information, 20460.

Van der Waal, B. (2009). The Influence of Acacia Mearnsii Invasion on Soil Properties in the Kouga Mountains, Eastern Cape, South Africa. Master’s Thesis. Grahamstown: Rhodes University. Online.

Google Scholar

Van Reeuwijk, L. P. (2002). Procedures for Soil Analysis. International Soil Reference and Information Center. Technical paper 96th edn. doi:10.1145/583810.583831

CrossRef Full Text | Google Scholar

Walkeley, A., and Black, C. A. (1934). An Examination of Degtjareff Method for Determining Soil Organic Matter and Proposed Modification of the Chromic Acid Titration Method. Soil Sci. 37, 29–38. doi:10.1097/00010694-193401000-00003

CrossRef Full Text | Google Scholar

Wang, L., Gan, Y., and Wiesmeier, M. (2018). Grazing Exclusion. An Effective Approach for Naturally Restoring Degraded Grasslands in Northern China. L China. Land Degrad. Dev. 1-18. doi:10.1002/ldr.3191

Google Scholar

Wang, Z., Johnson, D. A., Rong, Y., and Wang, K. (2016). Grazing Effects on Soil Characteristics and Vegetation of Grassland in Northern China. Solid Earth 7, 55–65. doi:10.5194/se-7-55-2016

CrossRef Full Text | Google Scholar

Wilke, B. M. (2005). “Determination of Physical and Chemical Soil Properties,” in Manual of Soil Analysis: Monitoring and Assessing Soil Bioremediation (3). Editors R. Margesin, and F. Schinner (Berlin: Springer-Verlag).

Google Scholar

Wolka, K., Moges, A., and Yimer, F. (2011). Effects of Level Soil Bunds and Stone Bunds on Soil Properties and its Implications for Crop Production: the Case of Bokole Watershed, Dawuro Zone, Southern Ethiopia. As 02, 357–363. doi:10.4236/as.2011.23047

CrossRef Full Text | Google Scholar

Yan, L., Zhou, G., and Zhang, F. (2013). Effects of Different Grazing Intensities on Grassland Production in China: a Meta-Analysis. PLoS ONE 8 (12), 1466–1475. doi:10.1371/journal.pone.0081466

CrossRef Full Text | Google Scholar

Yang, Y. H., and Luo, Y. Q. (2011). Isometric Biomass Partitioning Pattern in Forest Ecosystems: Evidence from Temporal Observations during Stand Development. J. Ecol. 99, 431–437. doi:10.1111/j.1365-2745.2010.01774.x

CrossRef Full Text | Google Scholar

Yang, Z., Xiong, W., Xu, Y., Jiang, L., Zhu, E., Zhan, W., et al. (2016). Soil Properties and Species Composition under Different Grazing Intensity in an Alpine Meadow on the Eastern Tibetan Plateau, China. Environ. Monit. Assess. 188, 678–688. doi:10.1007/s10661-016-5663-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yigini, Y., and Panagos, P. (2016). Assessment of Soil Organic Carbon Stocks under Future Climate and Land Cover Changes in Europe. Sci. Total Environ. 557-558, 838–850. doi:10.1016/j.scitotenv.2016.03.085

PubMed Abstract | CrossRef Full Text | Google Scholar

Yun, W., and Wesche, K. (2016). Vegetation and Soil Responses to Livestock Grazing in Central Asian Grasslands: a Review of Chinese Literature. Biodivers Conserv 25 (12), 2401–2240. doi:10.1007/s10531-015-1034-1

CrossRef Full Text | Google Scholar

Zhang, H., Xu, M., Zhang, W., and He, X. (2009). Factors Affecting Potassium Fixation in Seven Soils under 15-year Long-Term Fertilization. Sci. Bull. 54, 1773–1780. doi:10.1007/s11434-009-0164-9

CrossRef Full Text | Google Scholar

Zhang, X., Liu, M., Zhao, X., Li, Y., Zhao, W., Li, A., et al. (2018). Topography and Grazing Effects on Storage of Soil Organic Carbon and Nitrogen in the Northern China Grasslands. Ecol. Indicators 93, 45–53. doi:10.1016/j.ecolind.2018.04.068

CrossRef Full Text | Google Scholar

Zhao, Y., Peth, S., Hallett, P., Wang, X., Giese, M., Gao, Y., et al. (2011). Factors Controlling the Spatial Patterns of Soil Moisture in a Grazed Semi-arid Steppe Investigated by Multivariate Geostatistics. Ecohydrol. 4, 36–48. doi:10.1002/eco.121

CrossRef Full Text | Google Scholar

Zhou, G., Zhou, X., He, Y., Shao, J., Hu, Z., Liu, R., et al. (2017). Grazing Intensity Significantly Affects Belowground Carbon and Nitrogen Cycling in Grassland Ecosystems: a Meta-Analysis. Glob. Change Biol. 23, 1167–1179. doi:10.1111/gcb.13431

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z. C., Gan, Z. T., Shangguan, Z. P., and Dong, Z. B. (2010). Effects of Grazing on Soil Physical Properties and Soil Erodibility in Semiarid Grassland of the Northern Loess Plateau (China). Catena 82, 87–91. doi:10.1016/j.catena.2010.05.005

CrossRef Full Text | Google Scholar

Zhu, L., Johnson, D. A., Wang, W., Ma, L., and Rong, Y. (2015). Grazing Effects on Carbon Fluxes in a Northern China Grassland. J. Arid Environments 114, 41–48. doi:10.1016/j.jaridenv.2014.11.004

CrossRef Full Text | Google Scholar