Soil anti-scourability characteristics and influencing factors under different land use types in the black soil region of Northeast China

To explore the characteristics of soil anti-scourability under different land use types in the black soil region of Northeast China, surface soil samples (0–10 cm) from three main land use types in a typical black soil region, including natural grassland (uncultivated), farmland (converted from natural grassland), and the artificial forestland (grain plots converted into forests), were selected as the research objects. Soil anti-scourability test was adopted in the study, through which the soil anti-scourability coefficient was calculated and obtained. The study analyzed the differences in soil anti-scourability coefficients among various land use types, clarified the effects of soil physicochemical properties and root characteristics on soil anti-scourability, and thereby identified the key factors influencing soil anti-scourability. The results showed that (i) the soil loss amount under different land use types changes drastically in the first 2 minutes of scouring and gradually stabilizes with the extension of erosion time. Overall, the soil loss of natural grassland is the smallest, while that from farmland is the largest. (ii) The soil anti-scourability under different land use types is as follows: natural grassland > the artificial forestland > farmland. Compared with farmland, the artificial forestland significantly improves soil anti-scourability and notably reduces soil loss amount. (iii) The root geometric characteristic parameters of natural grassland and the artificial forestland are significantly larger than those of farmland, and the root volume density has the greatest impact on the soil anti-scourability coefficient. For the three land use types, the effects of fine roots of different diameter classes on soil anti-scourability are all shown as follows: 1.5~2 mm > 1~1.5 mm > 0.5~1 mm > 0~0.5 mm. The study can provide a theoretical basis for evaluating the soil and water conservation capacity of vegetation in the black soil region of Northeast China.

1 Introduction

Black soils have made tremendous contributions to the improvement of global grain production capacity and the safeguarding of food security. In the black soil region of Northeast China, soil erosion has developed rapidly over the past 50 years, with the area of soil and water loss continuously expanding (1, 2). Concurrently, the quality of black soil has exhibited a trend of gradual decline (3). However, in recent years, to restore black soil fertility, the “Grain for Green” program has become an important measure, which has led to the emergence of different land use types in the black soil region. Notably, land use and management regimes are the most direct factors affecting changes in soil quality (4, 5). Different land use types can influence soil anti-scourability by affecting soil physical and chemical conditions, soil fertility, and soil permeability. Specifically, soil physical, chemical, and biological properties vary significantly under different land use types (6, 7, 8). Consequently, reasonable land use types can improve soil quality, enhance water and fertilizer retention capacities, and enhance soil resilience to external environmental changes, while unreasonable land use types can damage soil structure, cause nutrient loss, and exacerbate erosion and degradation (9, 10).

Soil anti-scourability, defined as the ability of soil to resist runoff-induced mechanical damage, is a key component of soil erosion resistance (11–13), primarily dependent on the cementation force between soil particles and microstructures, as well as the resistance of soil aggregates to runoff scouring and dispersion (14). Soil physicochemical properties (inherent traits including texture structure and looseness (15)) and vegetation roots are critical factors regulating soil anti-scourability and its variation. Different land use types alter soil physicochemical properties (e.g., particle composition, organic matter content, pH) via variations in vegetation cover and management practices (e.g., tillage, fertilization), thereby affecting soil anti-scourability (16–18). Vegetation roots enhance soil anti-scourability by improving soil aggregate stability, physically binding soil mass, and protecting soil nutrients. As the key factor influencing soil anti-scourability (19), root traits (e.g., fine root proportion, root length density, root surface area) are closely related to soil physical and hydrological properties (20), with soil anti-scourability showing a highly significant positive correlation with root biomass and density (21, 22) and soil detachment rate decreasing exponentially with increasing root length or root density (23). Mechanistically, forest and grass roots increase soil shear strength by modifying soil particle internal friction angle and cohesion (24), while root exudate adhesion to soil particles is a primary mechanism enhancing root-mediated soil network protection (25). Additionally, roots promote soil organic matter content and stability, further improving aggregate stability (8, 26, 27).

Therefore, the topsoil (0–10 cm) under three main land use types, including natural grassland (uncultivated), farmland (converted from natural grassland), and the artificial forestland (grain plots converted into forests), was selected as the research object in the typical black soil region. The objectives of the study are to analyze the effects of soil physical and chemical properties and root characteristics on soil anti-scourability, clarify the differences in soil anti-scourability coefficients among different land use types, and identify the key factors influencing soil anti-scourability. Ultimately, the research findings will provide a theoretical basis for studying the soil and water conservation capacity of vegetation in the black soil region of Northeast China.

2 Materials and methods

2.1 Description of the study area

The study area is located in Keshan County, Heilongjiang Province (125°07’40”E -125°37’30”E, 48°11’15”N-48°24’07”N), situated at the western foot of the Xiao Hinggan Mountains and the northeastern part of the Songnen Plain, characterized by a hilly and rolling terrain. The total area is approximately 363,200 hectares. This region has a cold-temperate subhumid monsoon climate, with an average annual temperature of 1.3°C and an annual precipitation of around 500 mm, mostly concentrated in summer and autumn. The main land use type in this area is farmland, accounting for about 80.60% of the total area, followed by forestland (approximately 6.43%) and natural grassland (approximately 9.86%). The dominant soil type in this region is black soil, and the soil texture is silty clay loam.

2.2 Sampling protocol

Three land-use types were selected as the research objects in the study. Specifically, natural grassland represents the earliest documented land use type in the black soil region of Northeast China, with no direct human disturbance, and its main dominant species are Artemisia annua Linn and Carex spp. The farmland was converted from natural grassland and used for food crop cultivation, with maize as the primary crop. The artificial Pinus sylvestris plantation had a prior cultivation history; following this, Pinus sylvestris was artificially planted, and it has now developed into a mature forest with a stand age of 30 years.

For each land-use type, three quadrats (20 m×20 m) were randomly established, with three sampling points arranged along the diagonal of each quadrat. After removing the 0–5 cm thick layer of undecomposed litter, cutting ring and aluminum box samples were collected at each sampling point to determine soil bulk density, soil moisture content, and total porosity. To ensure the representativeness of soil chemical properties and minimize microscale spatial heterogeneity, a composite soil sample was prepared for each sampling point (along the diagonal of each quadrat) prior to chemical analysis. Subsequently, 500 g of undisturbed soil was collected from each sampling point to measure the mass percentage of soil aggregates, soil mechanical composition, soil pH, soil organic matter content, soil total phosphorus content, and soil total nitrogen content. Furthermore, surface-layer undisturbed soil was collected using an undisturbed soil sampler (diameter 16 cm × height 6 cm) at each sampling point for the determination of soil anti-scourability. After intact soil samples were transported to the laboratory, they were manually broken along the natural soil texture into clods with a diameter of less than 1 mm. Following air-drying, the wet-sieving method was adopted to determine the particle size distribution of soil aggregates (28). The clay content is 43.04 ± 0.26%, the silt content is 50.21 ± 0.60%, the sand content is 6.75 ± 0.68%. The basic information of the study plots is presented in Table 1.

Table 1

Table 1. Basic information and soil aggregate mass percentage of study plots.

2.3 Determination of soil physical and chemical properties

Soil bulk density, soil moisture content, total porosity, and soil aggregates were determined by conventional methods. Soil mechanical composition was measured via the pipette method. Soil pH was determined using the potentiometric method. Soil organic carbon content and total nitrogen content were analyzed with an elemental analyzer. Soil total phosphorus content was determined by the sulfuric acid-perchloric acid digestion-molybdenum antimony anti-colorimetric method.

2.4 Determination of root geometric characteristic parameters

After the soil anti-scouribility test, the soil clods were repeatedly rinsed through a fine sieve mesh (pore size 0.25 mm) to obtain clean plant roots, which were subsequently scanned with a root scanner and then analyzed using a root analysis system (EPSON LA 2400; 29) according to four root diameter classes (0-0.5 mm, 0.5-1.0 mm, 1.0-1.5 mm, and 1.5-2.0 mm) to acquire key plant root parameters including root length, root surface area, and root volume. Specifically, root length density (RLD) refers to the root length contained per unit volume of soil and can reflect the degree of root penetration and entanglement in the soil, root surface area density (RSD) is the root surface area per unit volume of soil and reflects the closeness of the root-soil interface, and root volume density (RVD) reflects the space occupied by plant roots (22), all of which (RLD, RSD, and RVD) were computed using Equations 1–3.

$\begin{array}{ll}RLD(cm/c{m}^3)=\frac{L}{D}& (1)\end{array}$

$\begin{array}{ll}RSD(c{m}^2/c{m}^3)=\frac{S}{D}& (2)\end{array}$

$\begin{array}{ll}RVD(cm/c{m}^3)=\frac{V}{D}& (3)\end{array}$

Where L is the root length (cm); S is the root surface area (cm²); V is the root volume (cm³); D is the volume of the sampler (cm³).

2.5 Soil anti-scourability test

The runoff scouring device is composed of a scouring tank, a flow-stabilizing tank, and an upper water supply device (Figure 1). The scouring tank has a dimension of 126 cm (length) × 20 cm (width), and a cylindrical soil tank with a diameter of 16 cm and a height of 6 cm is embedded in the lower half, 10 cm above the bottom, for observing the erosion process. The flow-stabilizing tank measures 20 cm (length) × 10 cm (width) × 16 cm (height), and a flow-retarding porous plate is installed in the middle, which mainly functions to maintain the stability of water flow. The water supply device consists of a water supply pipe, a constant head tank, an overflow pipe, and a connecting water pipe.

Figure 1

Figure 1. Soil scouring device.

Prior to the scouring test, the soil samples were placed in a flat-bottomed water basin with the water level not exceeding the upper edge of the sampler, and soaked for 24 h to achieve capillary water saturation; subsequently, the saturated undisturbed soil was placed on an iron stand for 8 h to remove gravitational water. The test parameters were set as follows: the scouring slope was fixed at 6°; the scouring discharge (unit flow rate) was calculated based on the maximum runoff generated by local standard runoff plots (20 m×5 m), determined to be 30 L/min; the scouring duration was set to 15 min, as preliminary tests confirmed that runoff sediment yield could reach a stable state within this period. During the scouring test, runoff-sediment samples were collected every 1 min with sampling buckets for the first 3 min after runoff initiation, and then every 2 min thereafter, totaling 9 samples. After the test, all sampling buckets were allowed to stand for clarification; once the sediment was completely settled, the supernatant was discarded, and the remaining sediment-water mixture was transferred to iron boxes. These boxes were dried in an oven at 105°C for 8 h to determine the dry mass of sediment (g), which was used to calculate the soil anti-scourability coefficient (AS). Defined as the volume of water required to erode 1 g of soil (unit: L/g), a higher AS value indicates stronger soil anti-scourability (30), and the specific calculation was performed using Equation 4.

$\begin{array}{ll} AS(L/g)=\frac{f\times t}{w}& (4)\end{array}$

Where AS is the soil anti-scouring coefficient (L/g); f is the flow discharge (L/min); t is the scouring duration (min); w is the dry mass of sediment (g).

2.6 Statistical analysis

Excel and Origin 2021 software were used for empirical equation fitting and data processing, while SPSS 21.0 was employed to perform correlation analysis and regression analysis on the data. Both the soil anti-scourability coefficient of different land-use types (Figure 2) and the correlation between soil physical and chemical properties and soil anti-scourability coefficient (Figure 3) were based on a sample size of n = 27, and all data followed a normal distribution.

Figure 2

Figure 2. Soil anti-scourability coefficient of different land-use types. Error bars represent standard deviation (SD).

Figure 3

Figure 3. Correlation between soil physical and chemical properties and soil anti-scourability coefficient.

3 Results

3.1 Soil physical and chemical properties

The soil physical and chemical properties under different land use types are as shown in Table 2. The soil bulk density of natural grassland was the smallest, while that of farmland was the largest, which is 1.781 and 1.199 times higher than that of natural grassland and artificial forest, respectively (P < 0.05). Among the three land-use types, the soil moisture content of natural grassland was the highest, reaching 70.837%, which was 3.308 and 2.914 times that of the artificial forestland and farmland, respectively, with significant differences from the other two land-use types (P < 0.05). The total porosity of natural grassland was significantly higher than that of the other two land-use types (P < 0.05), being 1.215 and 1.240 times that of the artificial forestland and farmland, respectively. Among the three land-use types, the mechanical composition was dominated by silt, with contents ranging from 50.207% to 60.227%, followed by clay and sand. As shown in Table 2, the soil texture of natural grassland and the artificial forestland was silty clay loam, while that of cultivated land was silty clay.

Table 2

Table 2. The soil physical and chemical properties under different land use types.

Overall, the soil pH of the three land-use types was weakly acidic, among which farmland had the highest pH value and the artificial forestland had the lowest. For soil organic matter (SOM) content and total nitrogen (TN) content across the three land-use types, both ranked in the order of natural grassland > the artificial forestland > farmland. Specifically, the SOM content of natural grassland was 2.146 and 2.632 times that of the artificial forestland and farmland, respectively, while its TN content was 2.107 and 2.636 times that of the two land-use types mentioned above, respectively. The total phosphorus (TP) content of the three land-use types was in the order of natural grassland > farmland > the artificial forestland, with the TP content of natural grassland being 1.239 and 1.093 times that of the artificial forestland and farmland, respectively.

3.2 Root geometric characteristics

The variations of soil root geometric characteristic parameters under different land-use types are shown in Table 3. Specifically, the root length density (RLD), root surface area density (RSD), and root volume density (RVD) of natural grassland all followed the order of 0~0.5 mm > 0.5~1 mm > 1~1.5 mm > 1.5~2 mm. Among these, the RLD of the 0~0.5 mm diameter class was significantly different from that of other diameter classes (P < 0.05), while no significant differences were observed in RSD and RVD among the four diameter classes. The variations of root geometric characteristic parameters in the artificial forestland and farmland were consistent with those in natural grassland. For all four diameter classes, the RLD, RSD, and RVD of natural grassland were significantly higher than those of the artificial forestland and farmland (P < 0.05).

Table 3

Table 3. The variations of soil root geometric characteristic parameters under different land-use types.

Furthermore, under three land-use types, for fine roots in the 0~0.5 mm diameter class, their root length density (RLD) accounted for 82.906%~89.945% of the total root RLD, root surface area density (RSD) for 55.414%~69.595% of the total RSD, and root volume density (RVD) for 35.714%~42.857% of the total RVD. For fine roots in the 0~1 mm diameter class, their RLD accounted for 96.764%~98.700% of the total RLD, RSD for 85.987%~93.216% of the total RSD, and RVD for 71.429%~85.714% of the total RVD. Thus, fine roots in the 0~1 mm diameter class constituted the dominant component of the root system.

3.3 Soil loss characteristics

The soil loss rates of different land-use types are shown in Figure 4. Overall, the soil loss rate of natural grassland was 1.947 g/min during the scouring time of 0~1 min, decreased rapidly to a lower level during 1~2 min, and tended to stabilize after 6 min of scouring. For the artificial forestland, the soil loss rate was 4.421 g/min during 0~1 min, dropped sharply during 1~2 min, and became stable after 10 min of scouring. As for forestland, the soil loss rate reached 25.249 g/min during 0~1 min, decreased rapidly to only 2.042 g/min during 1~2 min, and stabilized after 2 min of scouring. In summary, the soil loss rate of all three land-use types peaked during the 0~1 min scouring time. Among them, cultivated land had the highest soil loss rate, which was 12.968 and 5.711 times that of natural grassland and the artificial forestland, respectively. Additionally, with the increase of scouring time, the soil loss rates of the three land-use types gradually tended to stabilize.

Figure 4

Figure 4. The soil loss rates of different land-use types.

3.4 Soil anti-scourability characteristics

As shown in Figure 2, the soil anti-scourability coefficients of natural grassland and the artificial forestland were both higher than that of farmland, with significant differences from farmland (P < 0.05). Specifically, the anti-scourability coefficient of natural grassland was 152.498 L/g, which was 117.487 times that of farmland; the anti-scourability coefficient of the artificial forestland was 136.598 L/g, being 105.237 times that of farmland.

3.5 Correlation between various factors and soil anti-scourability coefficient

As revealed by the correlation analysis between soil physical and chemical properties and soil anti-scourability coefficient under the three land-use types (Figure 3), soil silt content showed an extremely significant positive correlation with soil anti-scourability coefficient (P < 0.01), while soil sand content and bulk density exhibited an extremely significant negative correlation with the soil anti-scourability coefficient (P < 0.01). Soil pH was significantly negatively correlated with soil anti-scourability coefficient (P < 0.05), and soil moisture content, total porosity, organic matter content, and total nitrogen (TN) content were significantly positively correlated with the soil anti-scourability coefficient (P < 0.05). No significant correlation was observed between soil clay content, total phosphorus (TP) content and soil anti-scourability coefficient. The order of correlation strength from strongest to weakest was: silt content > sand content > soil bulk density > total nitrogen content > organic matter content > total porosity > soil moisture content > soil pH.

As revealed by the correlation analysis between root characteristic parameters and soil anti-scourability coefficient under the three land-use types (Figure 5), the root length density (RLD), root surface area density (RSD), and root volume density (RVD) of fine roots in the 0.5~1 mm, 1~1.5 mm, and 1.5~2 mm diameter classes showed an extremely significant positive correlation with the soil anti-scourability coefficient (P<0.01). In contrast, the RLD, RSD, and RVD of fine roots in the 0~0.5 mm diameter class exhibited a significantly positive correlation with the soil anti-scourability coefficient (P<0.05). In terms of correlation strength, the order from strongest to weakest was: RLD (1.5-2) > RSD _(1.5-2) > RVD _(1.5-2) > RVD _(1-1.5) > RSD _(1-1.5) > RLD _(1-1.5) > RVD _(0.5-1) > RSD _(0.5-1) > RLD _(0.5-1) > RVD _(0-0.5) > RSD _(0-0.5) > RLD _(0-0.5).

Figure 5

Figure 5. Correlation between soil root geometric characteristics and soil anti-scourability coefficient.

4 Discussion

4.1 Effects of soil physical and chemical properties on soil anti-scourability under different land-use types

Soil bulk density and porosity are the main indicators for evaluating soil compaction. The higher the soil compaction, the less likely the soil is to be eroded and transported, and the stronger its anti-scourability. Among different land-use types, soil bulk density followed the order of farmland > the artificial forestland > natural grassland, while total porosity showed the opposite trend (Figure 2). One reason is the minimal human disturbance (5). In addition, plant roots exert network consolidation and biochemical effects, which promote the formation of soil aggregates (29, 31). These aggregates further loosen the soil structure, increase soil porosity, and ultimately reduce soil bulk density (32). One reason is the minimal human disturbance; in addition, the network consolidation and biochemical effects of plant roots promote the formation of soil aggregates, loosen the soil, increase soil porosity, and reduce soil bulk density (31, 32). On the other hand, compared with forestland, the vegetation in natural grassland has more sufficient light and water conditions, leading to active soil microorganisms and high organic matter content, which is manifested in a loose soil structure and reduced bulk density (30).

Furthermore, farmland had the highest bulk density, with significant differences from the other two land-use types, which is consistent with the conclusion of Schäffer et al. (33). This is primarily attributed to mechanical tillage. Tillage disturbs the soil matrix, breaking down large soil aggregates into smaller particles (34). The loss of aggregate structure then increases soil bulk density and reduces soil structural stability. The lowest porosity of farmland is partly due to underdeveloped root systems and the inhibition of soil biological activities by chemical fertilizers and pesticides. Additionally, agricultural practices such as irrigation and sowing, coupled with frequent trampling, cause soil compaction and hardening, impairing the connectivity of soil pores (35).

Under different land-use types, soil pH followed the order of farmland > natural grassland > the artificial forestland, while both organic matter content and total nitrogen (TN) content showed the trend of natural grassland > the artificial forestland > farmland. Combined with Table 2, the soil organic matter content and TN content of farmland were significantly lower than those of natural grassland and the artificial forestland, indicating severe soil nutrient loss in farmland. This is because artificial mechanical tillage damages the soil structure, which further reduces the physical protection of soil organic matter, thereby increasing the mineralization and decomposition of soil organic carbon and nitrogen (36). Meanwhile, a decrease in organic matter content impairs aggregate formation and stability, thereby reducing the soil’s ability to resist erosion. This ultimately leads farmland into a vicious circle of gradual degradation (37).

In contrast, the significant increase in soil organic matter content and TN content of the artificial forestland (converted from farmland to forest) indicates that due to the reduction in human activities, litter in the forestland can return to the soil under natural conditions. The litter contains a large amount of tannins, resins, and lignin, whose decomposition produces acidic substances, improving the utilization efficiency of soil nutrients. Additionally, the high humus content in forestland contributes to the acidification of the topsoil (5). Moreover, the type of litter also affects soil pH. This suggests that the large amount of litter returning to the soil in the artificial forestland improves soil texture, thereby enhancing soil anti-scourability to a certain extent (38). Natural grassland has a well-developed root system, and the root network formed by intertwined fine roots facilitates the renewal of the ecosystem and accelerates its rate, enabling better production, storage, and protection of organic matter and nitrogen. However, there are still significant differences in soil organic matter content and TN content between the artificial forestland and the undisturbed natural grassland.

4.2 Effects of root geometric characteristics on soil anti-scourability under different land-use types

The enhancing effect of plant roots on soil anti-scourability is an important aspect of vegetation control of soil and water loss (39). Plant roots can enhance soil anti-scourability through network connection and root-soil bonding effects (40). On the one hand, roots penetrate and entangle in the soil to form a root network, which consolidates the soil mass. This mechanical network effect stabilizes the soil structure by locking soil particles between root strands, making it less likely to be scoured and fragmented by water flow. On the other hand, roots improve soil physicochemical properties by exuding organic compounds (e.g., mucilage), which act as cementing agents to promote aggregate formation. This optimized soil structure further enhances anti-scourability (29). As shown in Table 3, the changes in soil root geometric characteristic parameters under different land-use types indicate that the four diameter class parameters of root length density (RLD), root surface area density (RSD), and root volume density (RVD) all followed the order of natural grassland > the artificial forestland > farmland. It can be seen that the roots of natural grassland have high density in the topsoil and large contact area with the soil mass, forming a good root network. The consolidation, network connection, root-soil bonding, and biochemical effects of vegetation roots weaken the ability of runoff to scour the soil, indicating that vegetation with abundant roots can significantly improve soil anti-scourability (41). The root geometric characteristic parameters of the artificial forestland are 3 to 5 times those of farmland, suggesting that the conversion of farmland to forest can increase soil root parameters and enhance root consolidation capacity.

Correlation analysis between root geometric parameters and soil anti-scourability within the same diameter class shows the following trend: root volume density (RVD) > root surface area density (RSD) > root length density (RLD) for all diameter classes except 1.5~2 mm fine roots, which showed RLD > RSD > RVD. Among these parameters, RVD had the highest average correlation coefficient. This is because RVD represents the proportion of space occupied by roots in the soil—the larger the RVD, the greater the space occupied by roots. Due to the small size of roots themselves, minor changes in RVD can induce significant changes in soil anti-scourability (14). An excessively small RVD has a limited effect on improving soil anti-scourability. When RVD increases, the vegetation roots become dense, which greatly enhances soil resistance to water scouring and erosion (42). Overall, the 1.5~2 mm diameter fine roots have a closer relationship with the soil anti-scourability coefficient, and all root characteristic parameters show a positive correlation with the coefficient. This indicates that the longer the roots, the larger the surface area they occupy, and the greater their volume, the stronger the soil anti-scourability.

5 Conclusions

Soil loss under different land-use types changed drastically within the first scouring time of 2 min. With the extension of scouring time, it tended to stabilize after scouring time of 6 min. Under the same conditions, natural grassland had the smallest soil loss, while farmland had the largest.

Different land-use types result in differences in soil anti-scourability, following the order of natural grassland > the artificial forestland > farmland. The artificial forestland significantly improved soil anti-scourability and exhibited a significant reduction effect on soil loss.

Root geometric characteristic parameters of vegetation in natural grassland and the artificial forestland were significantly higher than those in farmland, and root volume density had the greatest impact on the soil anti-scourability coefficient. For the three land-use types, the impact of fine roots of different diameter classes on soil anti-scourability followed the order: 1.5~2 mm > 1~1.5 mm > 0.5~1 mm > 0~0.5 mm. It is suggested that the characteristic indicators of 1.5~2 mm diameter fine roots be used as an indicator for evaluating the soil and water conservation capacity of vegetation in the black soil region of Northeast China.

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

LB: Resources, Writing – original draft, Visualization, Project administration, Formal analysis, Conceptualization, Methodology, Investigation. JL: Investigation, Writing – review & editing, Data curation, Validation, Software, Supervision, Funding acquisition.

Funding

The author(s) declared that financial support was not received for this work and/or its publication.

Conflict of interest

Author LB was employed by the company Songliao Water Resources & Hydropower Development Co.Ltd,

The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Wang H, Yang S, Wang Y, Gu Z, Xiong S, Huang X, et al. Rates and causes of black soil erosion in Northeast China. Catena. (2022) 214:106250. doi: 10.1016/j.catena.2022.106250

Crossref Full Text | Google Scholar

2. Jiao P, Ou Y, Pang SJ, Yan BX, Zhang Y, Xu WX, et al. Impacts of landscape factors on gully retreat and its morphological characteristics in hilly areas of Northeast China. Soil Till. Res. (2025) 248:106434. doi: 10.1016/j.still.2024.106434

Crossref Full Text | Google Scholar

3. Duan XW, Xie Y, Ou TH, and Lu HM. Effects of soil erosion on long-term soil productivity in the black soil region of northeastern China. Caten. (2011) 87:268–75. doi: 10.1016/j.catena.2011.06.012

Crossref Full Text | Google Scholar

4. Han SJ, Wang EH, Chen XW, and Fu Y. Accumulation of SOM fractions in croplands and plantations converted from cropland with black soil. Land. Degrad. Dev. (2022) 33:638–48. doi: 10.1002/ldr.4187

Crossref Full Text | Google Scholar

5. He YY, Zhang C, Xue TT, Chen XW, and Fu Y. Effects of land-use changes on soil organic nitrogen fractions in the black soil region of Northeast China. Soil Use Manage. (2023) 39:805–16. doi: 10.1016/j.ecolind.2021.108443

Crossref Full Text | Google Scholar

6. Raciti SM, Groffman PM, Jenkins JC, Pouyat RV, Fahey TJ, Pickett STA, et al. Accumulation of carbon and nitrogen in residential soils with different land-use histories. Ecosystems. (2011) 14:287–97. doi: 10.1007/s10021-010-9409-3

Crossref Full Text | Google Scholar

7. Rahimi MR, Ayoubi S, and Abdi MR. Magnetic susceptibility and Cs-137 inventory variability as influenced by land use change and slope positions in a hilly, semiarid region of west-central Iran. J Appl Geophys. (2013) 89:68–75. doi: 10.1016/j.jappgeo.2012.11.009

Crossref Full Text | Google Scholar

8. Koyama T, Enggrob KL, Rasmussen J, Martins JT, and Peixoto L. Substrate quantity and quality affect microbial carbon use efficiency and priming effects of root exudates investigated with micro dialysis. Soil Biol Biochem. (2025) 209:109869. doi: 10.1016/j.soilbio.2025.109869

Crossref Full Text | Google Scholar

9. Kelishadi H, Mosaddeghi MR, Hajabbasi MA, and Ayoubi S. Near-saturated soil hydraulic properties as influenced by land use management systems in Koohrang region of central Zagros, Iran. Geoderma. (2014) 213:426–34. doi: 10.1016/j.geoderma.2013.08.008

Crossref Full Text | Google Scholar

10. Ayoubi S and Moazzeni Dehaghani S. Identifying impacts of land use change on soil redistribution at different slope positions using magnetic susceptibility. Arab. J GeoSci. (2020) 13:1–11. doi: 10.1007/s12517-020-05383-x

Crossref Full Text | Google Scholar

11. Wang B, Zheng FL, Mathias JR, and Darboux F. Soil erodibility for water erosion: A perspective and Chinese experiences. Geomorphology. (2013) 187:1–10. doi: 10.1016/j.geomorph.2013.01.018

Crossref Full Text | Google Scholar

12. Borrelli P, Robinson DA, Fleischer LR, Lugato E, Ballabio C, Alewell C, et al. An assessment of the global impact of 21st century land use change on soil erosion. Nat Commun. (2017) 8:20135. doi: 10.1038/s41467-017-02142-7

PubMed Abstract | Crossref Full Text | Google Scholar

13. Fu Y, Wang DC, Sun WP, and Guo MM. Impacts of grass planting density and components on overland flow hydraulics and soil loss. Land. Degrad. Dev. (2023) 34:234–49. doi: 10.1002/ldr.4456

Crossref Full Text | Google Scholar

14. Ma R, Zheng ZC, Li TX, He SQ, Zhang XZ, Wang YD, et al. Temporal variation of soil erosion resistance on sloping farmland during the growth stages of maize (Zea mays L.). Hydrol. Process. (2021) 35:e14353. doi: 10.1002/hyp.14353

Crossref Full Text | Google Scholar

15. Schlüter S, Wu M, Phalempin M, Philipp L, Blagodatskaya E, Reitz T, et al. Divergence in physical, chemical, and biological soil properties caused by different long-term bare fallow management and natural succession. Geoderma. (2025) 459:117361. doi: 10.1016/j.geoderma.2025.117361

Crossref Full Text | Google Scholar

16. Havaee S, Ayoubi S, Mosaddeghi MR, and Keller T. Impacts of land use on soil organic matter and degree of compactness in calcareous soils of central Iran. Soil Use Manage. (2014) 30:2–9. doi: 10.1111/sum.12092

Crossref Full Text | Google Scholar

17. Zolfaghari Z, Mosaddeghi MR, Ayoubi S, and Kelishadi H. Soil atterberg limits and consistency indices as influenced by land use and slope position in Western Iran. J Mount. Sci. (2015) 12:1471–83. doi: 10.1007/s11629-014-3339-z

Crossref Full Text | Google Scholar

18. Geng TT, Wang ZH, Wang QK, Yang QS, and Wang JG. Effects of land-use type and soil depth on soil physicochemical properties, enzyme activities and bacterial community composition in black soil region of Northeast China. Catena. (2025) 258:109276. doi: 10.1016/j.catena.2025.109276

Crossref Full Text | Google Scholar

19. Reubens B, Achten WMJ, Maes WH, Danjon F, Aerts R, Poesen J, et al. More than biofuel? Jatropha curcas root system symmetry and potential for soil erosion control. J Arid. Environ. (2011) 75:201–5. doi: 10.1016/j.jaridenv.2010.09.011

Crossref Full Text | Google Scholar

20. Baets SD, Poesen J, Gyssels G, and Knapen A. Effects of grass roots on the erodibility of topsoils during concentrated flow. Geomorphology. (2006) 76:54–67. doi: 10.1016/j.geomorph.2005.10.002

Crossref Full Text | Google Scholar

21. Guo MM, Wang WL, Kang HL, and Yang B. Effect of natural vegetation restoration age on slope soil anti-scourability in gully region of Loess Plateau. Trans CSAE. (2018) 34:138–46. doi: 10.11975/j.issn.1002-6819.2018.22.017

Crossref Full Text | Google Scholar

22. Wang JF, Yang YF, Wang B, and Liu GB. Seasonal variations in soil erosion resistance under tap and fibrous root systems grasslands o-n the Chinese Loess Plateau. Geoderma. (2025) 458:117350. doi: 10.1016/j.geoderma.2025.117350

Crossref Full Text | Google Scholar

23. Gyssels G, Poesen J, Bochet E, and Li Y. Impact of plant roots on the resistance of soils to erosion by water: A review. Prog Phys Geography.: Earth Environ. (2005) 29:189–217. doi: 10.1016/B978-0-12-820106-0.00016-6

Crossref Full Text | Google Scholar

24. Guo MM, Wang WL, Wang TC, Wang WX, and Kang HL. Impacts of different vegetation restoration options on gully head soil resistance and soil erosion in loess tablelands. Earth Surf. Process. Landforms. (2020) 45:1049. doi: 10.1002/esp.4798

Crossref Full Text | Google Scholar

25. Yan SB, Yin LM, Dijkstra FA, Wang P, and Cheng WX. Pr-iming effect on soil carbon decomposition by root exudate surrogates: A meta-analysis. Soil Biol Biochem. (2023) 178:108955. doi: 10.1016/j.soilbio.2023.108955

Crossref Full Text | Google Scholar

26. An S, Mentler A, Mayer H, and Blum WEH. Soil aggregation, aggregate stability, organic carbon and nitrogen in different soil aggregate fractions under forest and shrub vegetation on the Loess Plateau, China. Catena. (2010) 81:226–33. doi: 10.1016/j.catena.2010.04.002

Crossref Full Text | Google Scholar

27. Ma J, Ma B, Li ZB, Wang CG, Zhang ZY, and Guo L. Effect of roots on the soil detachment process in grassland and shrubland. Land. Degrad. Dev. (2023) 34:3017–32. doi: 10.1002/ldr.4664

Crossref Full Text | Google Scholar

28. Yoder RE. A direct method of aggregate analysis of soils and a study of the physical nature of erosion losses. J Am Soc Agron. (1936) 28:337–51. doi: 10.2136/SSSAJ1936.036159950B1720010046X

Crossref Full Text | Google Scholar

29. Zhao YK, Fu Y, Xu JZ, Li Y, Zhang YP, Zhang XY, et al. Moisture and root systems affect the turnover of soil aggregates during various crop growth periods. Catena. (2025) 252:108840. doi: 10.1016/j.catena.2025.108840

Crossref Full Text | Google Scholar

30. Coonan EC, Kirkby CA, Kirkegaard JA, Amidy MR, Strong CL, and Richardson AE. Microorganisms and nutrient stoichiometry as mediators of soil organic matter dynamics. Nutr. Cycling. Agroecosyst. (2020) 117:273–98. doi: 10.1007/s10705-020-10076-8

Crossref Full Text | Google Scholar

31. Pinno BD and Wilson SD. Fine root response to soil resource heterogeneity differs between grassland and forest. Plant Ecol. (2013) 214:821–9. doi: 10.1007/s11258-013-0210-1

Crossref Full Text | Google Scholar

32. Li JY, Yuan XL, Ge L, Li Q, Li ZG, Wang L, et al. Rhizosphere effects promote soil aggregate stability and associated organic carbon sequestration in rocky areas of desertification. Agric. Ecosyst. Environ. (2020) 304:107126. doi: 10.1016/j.agee.2020.107126

Crossref Full Text | Google Scholar

33. Schäffer B, Stauber M, Müller R, and Schulin R. Changes in the m-acro-pore structure of restored soil caused by compaction beneath heavy a-gricultural machinery: a morphometric study. Eur J Soil Sci. (2007) 58:1062–73. doi: 10.1111/j.1365-2389.2007.00886.x

Crossref Full Text | Google Scholar

34. Liebhard G, Toth M, Stumpp C, Bodner G, Klik A, Zhang X, et al. Developing topsoil structure through conservation management to protect subsoil from compaction. Soil Till. Res. (2025) 253:10669. doi: 10.1016/j.still.2025.106669

Crossref Full Text | Google Scholar

35. Keller T, Sandin M, Colombi T, Horn R, and Dani O. Historical increase in agricultural machinery weights enhanced soil stress levels and adversely affected soil functioning. Soil Till. Res. (2019) 194:104293. doi: 10.1016/j.still.2019.104293

Crossref Full Text | Google Scholar

36. Fu W, Fan J, Wang S, Wang H, Dai ZJ, Zhao X, et al. Woody peat addition increases soil organic matter but its mineralization is affected by soil clay in the four degenerated erodible soils. Agric Ecosyst Environ. (2021) 318:107495. doi: 10.1016/j.agee.2021.107495

Crossref Full Text | Google Scholar

37. Chen S, Zhang G, and Wang C. Temporal variation in soil erodibility indicators of sloping croplands with different straw-incorporation rates. Soil Till. Res. (2025) 246:106340. doi: 10.1016/j.still.2024.106340

Crossref Full Text | Google Scholar

38. Zhao L, Yebra M, Van Dijk AIJM, and Cary GJ. Representing vapour and capillary rise from the soil improves a leaf litter moisture model. J Hydrol. (2022) 612:128087. doi: 10.1016/j.jhydrol.2022.128087

Crossref Full Text | Google Scholar

39. Zi R, Zhao LS, Fang Q, Qian XH, Fang FY, and Fan CH. Path analysis of the effects of hydraulic conditions, soil properties and plant roots on the soil detachment capacity of karst hillslopes. Catena. (2023) 228:107177. doi: 10.1016/j.catena.2023.107177

Crossref Full Text | Google Scholar

40. Ma F, Han XS, Liu L, Hao YF, and Lu XH. Is “soil anti-scourability” a rational methodology for assessing soil erosion resistance? Insights from a case study and a concise discussion. Soil Till. Res. (2025) 249:106481. doi: 10.1016/j.still.2025.106481

Crossref Full Text | Google Scholar

41. Shahriar MM, Wang JX, Alam S, and Patterson WB. Soil-binding ability of vegetation roots in enhancing erosion resistance of a shallow slope. Int J Geotech. Eng. (2016) 10:409–17. doi: 10.1080/19386362.2016.1168608

Crossref Full Text | Google Scholar

42. Helsen K, Vannoppen W, Honnay O, and Poesen J. Both below-ground and above-ground functional traits can help predict levee grassland root length density as a proxy for flow erosion resistance. J Vegetation. Sci. (2016) 27:1254–63. doi: 10.1111/jvs.12442

Crossref Full Text | Google Scholar