Evaluating critical hydrodynamic characteristics of rill formation on saturated soil slopes

Abstract

As a specific near surface hydrological condition, soil saturation can significantly affect the critical hydrodynamic characteristic and soil erosion rate of rill formation, leading to severe rill erosion. Nevertheless, few studies have investigated the characteristics of critical hydrodynamic parameters and their relationships with rill erosion rate under critical hydrodynamic conditions of rill formation on saturated soil slopes. Consequently, the quantification of critical hydrodynamic parameters and their effects on rill erosion rate under critical hydrodynamic conditions of rill formation on saturated soil slopes is of great significance for understanding the dynamic mechanism of rill formation and evolution and for predicting and controlling soil loss. In this study, indoor simulated rainfall experiments were performed and a new analytical model (V_c(NAM)) was applied to calculate the aforementioned critical parameters under a wide range of hydraulic conditions comprising five slope gradients (S_G) (2°, 5°, 10°, 15°, and 20°) and three rainfall intensities (R_I) (30, 60, and 90 mm/h). The results indicated that the new analytical model (V_c(NAM)) was suitable for estimating critical hydrodynamic parameters on saturated soil slopes. The critical flow velocity (V_c), the critical shear stress (τ_c), and the critical stream power (ω_c) apparently increased, whereas the critical rill length (L_c) decreased with the increase of slope gradients and rainfall intensities. Moreover, the erosion rate at the critical condition increased with decreasing L_c and increasing V_c, τ_c, and ω_c. Pearson correlation analysis indicated that τ_c and ω_c were significantly positively correlated, whereas L_c was negatively correlated with erosion rate under the critical conditions. Stepwise regression analysis revealed that the erosion rate under critical hydrodynamic conditions of rill formation could be well predicted by τ_c (R² = 0.83) with the linear model. The results provide an accurate model for evaluating critical conditions of rill formation and a basis for further understanding the intrinsic dynamic mechanism of rill formation on saturated soil slopes.