Evaluating critical hydrodynamic characteristics of rill formation on saturated soil slopes

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 2 = 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. • A new analytical model V c(NAM) was built to estimate critical hydrodynamic parameters on saturated soil slopes • The critical hydrodynamic of rill formation on saturated soil slopes was explored. • Shear stress was optimal hydrodynamic parameter to predicting erosion rates of rill initiation.