4.2. Rainfall Effects on Soil Aggregates Stability
After the rainfall for more than one month, the stability of water-stable aggregates was lower than that before the rainfall, and Shi et al. [
34] thought that the stability of soil aggregates will not be improved again after the drying process, which was also confirmed in this study. Among the seven natural rainfalls, every two rainfalls were spaced by several days, and the aggregates would be fully exposed to the sun and dried. The final result showed that the stability of aggregates was lower than that before the rainfall (
Table 6). Dimoyiannis et al. [
35] proved through continuous monitoring for two years that among many factors, rainfall and air temperature can strongly affect the dynamic change of soil aggregate stability, which is ultimately ascribed to the dynamic change of soil moisture in the final analysis. He et al. [
36] also showed that the dynamic change of MWD with time is negatively correlated with the dynamic change of soil moisture, and the change in the MWD value of water-stable aggregates in this study also follows this law. With the repeated rainfall, soil aggregates are frequently in wetted-dry-wetted cycles, and their stability finally declines. The spatial change of soil aggregate stability, an index used to evaluate the damage of aggregate resistance against external forces to the aggregate [
37], on the slope may result from the spatial heterogeneity of soil particles induced by water erosion. Soil particles are redistributed along the slope after being splashed by raindrops and washed by runoffs, thus changing the soil properties at different positions on the slope. In the research of Zhang et al. [
38], it was found that the relative content of clay in the surface layer of the slope increased along the runoff direction, but the deformation and destruction of the slope led to the decrease in the relative content of clay in sediments located at the lower part of the slope. Clay particles have large specific surface area and strong cation exchange ability, which can effectively gather stable aggregates [
39], and their reduction will inevitably weaken the stability of aggregates at the lower part of the slope.
As shown in
Figure 7, the difference values of aggregate stability parameters before and after rainfall under various erosion intensities could be observed. The GMD of water-stable aggregates increased by 0.005 and 0.002 in the non-eroded area and slightly eroded area, respectively, which could be regarded as unchanged compared with that before rainfall, and decreased significantly under other erosion intensities. The maximum value of MWD of water-stable aggregates was 0.433, and the minimum value was 0.16, which decreased significantly under all erosion intensities. With the increase of erosion intensity, the MWD decreased by 0.017, 0.071, 0.118, 0.125, 0.123, and 0.189, respectively, which was the same as that studied by Xia et al. [
40]. This may be because the aggregate is decomposed by the rainfall into differently sized particles, which block soil pores and form crusts, thus increasing surface runoffs [
41], and the more severely eroded areas are more seriously scoured. As shown in
Table 6, the PAD after rainfall was higher than that before rainfall, meaning that soil aggregates were easily dispersed when encountering rainfall or runoffs, thus reducing the stability of soil aggregates. The parameters of mechanically stable aggregates increased under different erosion intensities. Teh, C. B. S. [
10] pointed out that the stability of soil aggregates depends on the size of individual aggregates. In
Table 6, the large aggregates on the whole slope after rainfall were higher than those before rainfall, so the parameters of mechanically stable aggregates would increase under different erosion intensities. Based on the changes in parameters of mechanically stable aggregates and water-stable aggregates under different erosion intensities, it could be found that with the increase of erosion intensity, the stability parameters of aggregates changed extremely obviously, especially under the severe erosion intensity.
4.3. Relationship between Micro-Topographic Factors and Soil Aggregate Stability
The correlation between micro-topographic factors and stability parameters of soil aggregates was analyzed, and the results are shown in
Figure 8. The micro-topographic factors were all correlated with the stability parameters of soil aggregates. The slope was significantly negatively correlated with water-stable GMD (
p < 0.05) with a correlation coefficient of −0.21, and presented significant positive correlations with the surface roughness and the surface cutting degree (
p < 0.05), with correlation coefficients of 0.65 and 0.58, respectively. The surface roughness showed a significant positive correlation with the surface cutting degree (
p < 0.05) and a significant negative correlation with the plane curvature (
p < 0.05), with correlation coefficients of 0.81 and −0.33, respectively, and a significant negative correlation with water-stable GMD with a correlation coefficient of −0.17 (
p < 0.05). The surface cutting degree was significantly negatively correlated with the plane curvature and mechanically stable GMD (
p < 0.05), and the correlation coefficients were −0.41 and −0.16, respectively. The water-stable MWD showed a significant negative correlation with water-stable GMD (
p < 0.05), and it was significantly negatively correlated with both mechanically stable MWD and GMD (
p < 0.05). There was a significant positive correlation between mechanically stable MWD and mechanically stable GMD (
p < 0.05), and the correlation coefficient was 0.95.
According to correlation results between soil aggregate stability parameters and topographic factors in
Figure 8, it was found that every two parameters between micro-topographic factors and soil aggregate stability parameters were correlated but insignificantly. Hence, there might be a relationship between micro-topographic factors and soil stability parameters that cannot be explained by a single factor, and linear relationship and nonlinear relationship might exist at the same time. GAM was used to fit the relationship between topographic factors and soil stability parameters. Four soil aggregate stability parameters (water-stable aggregate MWD, water-stable aggregate GMD, mechanically stable aggregate MWD, and mechanically stable aggregate GMD) were used as response variables and four micro-topographic factors (slope, surface roughness, surface cutting degree, and plane curvature) were used to construct models, respectively, for multi-factor fitting (
Table 9). The results showed that the variance explained rate R
2 of the model fitting value to the response variable was 0.378–0.519, and the fitting result was of certain reference significance.
For water-stable aggregate MWD, the response variable was nonlinearly related to the slope and surface cutting degree, and linearly related to the surface roughness and plane curvature, among which the slope was the most important influencing factor. The water-stable aggregate GMD was also nonlinearly correlated with the slope and surface cutting degree, and linearly correlated with the surface roughness and plane curvature, but it was influenced by the surface roughness most intensely, and then by the slope. The mechanically stable aggregate MWD and the mechanically stable aggregate GMD were both nonlinearly related to the plane curvature and linearly related to other topographic factors, and the plane curvature was the biggest influencing factor.
Because the relationship between variables is nonparametric in GAM and cannot be described by simple mathematical formulas, the graphical smooth curve is a common way to observe the dependence between variables in GAM (
Figure 9). It could be seen from
Figure 9a,b that both MWD and GMD of water-stable aggregates would decrease with the increase of slope, and the relationship between GMD and slope was extremely significant (
p < 0.01), which might be attributed to the fact that the erosion and loss of surface soil increased with the increase of slope [
42], which led to the increase of runoff velocity, the enhanced scouring of aggregates, the increase of fragmentation rate, and the decrease of aggregation degree. Both MWD and GMD of water-stable aggregates would decrease with the increase of surface roughness, and the relationship between GMD and surface roughness was significant (
p < 0.05). Decomposition of soil surface aggregates not only results from the external force of external rainfall and runoffs but also from the wetting process of expansion and explosion caused by the difference between air entrained in aggregates and atmospheric pressure [
43]. Generally, the increase of surface roughness will reduce the flow velocity, increase the surface water storage capacity, and enhance infiltration [
44], which accelerates the wetting process of aggregates, leading to accelerated decomposition of aggregates and their decreased stability.
Figure 9b shows the relationship between the surface roughness and the stability parameters of soil aggregates. The change of the surface roughness value was significantly correlated with the MWD of water-stable aggregates and the MWD of mechanically stable aggregates (
p < 0.01). The changes in MWD and GMD of water-stable aggregates first increased and then decreased with the increase in the surface cutting degree, meaning that the stability of water-stable aggregates was high at the subsidence and weak at the uplift. The changes in MWD and GMD of water-stable aggregates also showed the same characteristics with the plane curvature, and the stability of aggregates was negatively correlated with the plane curvature, which was also reflected in the study of Nimmod et al. [
45]. This is because organic carbon can be more easily accumulated in the concave surface, and can facilitate soil particles to form aggregates and enhance the stability. In
Figure 9c,d, the parameters of mechanically stable aggregates were almost not different from the changes of topographic factors but nonlinearly correlated with the plane curvature, namely, it would decline with the increase of the plane curvature. As a whole, however, mechanical stability was of not too much reference significance to the evaluation of soil aggregate stability.
It is worth noting, however, that the small-scale micro-topography on the slope has no significant influence on the stability of soil aggregates, but it can reflect the erosion-deposition process on the slope, and it can also influence the factors that can affect the stability of soil aggregates, such as soil moisture, nutrients (SOM, SOC), soil mechanical composition, and microorganisms. Therefore, the correlation between the stability index and topographic factors discussed in this study can be understood as the correlation between the stability index and the spatial distribution complex of soil properties that affect the stability of soil aggregates, so multi-factor analysis should be strengthened in the follow-up research.