Rigid Vegetation Affects Slope Flow Velocity

: The mean slope flow velocity is critical in soil erosion models but the mechanism of its variation under rigid vegetation cover remains unclear. On natural slopes, vegetation grows predominantly perpendicular to the horizontal plane (BH), with some growing perpendicularly to the slope surface (BS); however, current research often neglects the effects of these two growth directions on the mean flow velocity. We conducted simulation experiments using different coverage levels, rigid vegetation, slope angles, and flow rates and showed that the flow rate and slope significantly influenced the mean flow velocity. As the coverage of rigid vegetation increased, the mean flow velocity increased more under conditions perpendicular to the horizontal plane (BH) and those perpendicular to the slope (BS). A model for predicting mean flow velocity was developed using vegetation equivalent roughness and the Manning formula, which accurately predicted flow velocity in different conditions. This study contributes to the refinement of slope flow theory and provides data that support soil and water conservation efforts.


Introduction
Soil erosion is a serious ecological and environmental problem that concerns many researchers worldwide [1][2][3].The basic unit of hydraulic erosion is slope runoff, for which the hydraulic properties are intimately linked to soil erosion.Thus, understanding these hydraulic characteristics is essential for studying the soil erosion process.Flow velocity is a fundamental and primary factor in the hydraulic properties of slope runoff and is used in the computation of additional hydraulic factors, including shear force and stream power, which are important parameters for analyzing soil detachment and sediment transport processes [4,5].Therefore, studying the flow velocity from slope runoff not only sheds light on the commencement of erosion, but also facilitates the creation of robust predictive models that can be used toward precisely tracking its progression [6].
The Loess Plateau in Northwestern China is characterized by loose soil [7] and significant precipitation variability [8].In this region, severe degradation of unsustainable land over time and the fragile ecological environment has caused some of the worst soil erosion worldwide.To mitigate the damage caused by soil erosion, the Chinese government implemented the Grain for Green policy in 2000 [9].Vegetation prevents soil erosion and slows desertification through sediment interception and litter accumulation [10,11].The properties of the underlying slope surface are one of the many factors that affect the flow velocity of slope runoff [12,13], making its mechanism complex.Thus, research is required to determine how vegetation affects slope runoff flow velocity.Liu et al. [14] conducted experiments on 30 slopes with varying degrees of vegetation coverage, including a complete mix of grass and shrubs, and showed that slopes with vegetation cover had a noticeably higher resistance coefficient than slopes without vegetation cover, indicating that vegetation reduced the slope flow velocity.Li and Pan [15] also found that Water 2024, 16, 2240 2 of 17 vegetation significantly increased the resistance coefficient of slopes and decreased flow velocity, while Zhang and Hu [16] found that removing the aboveground components of vegetation significantly increased the slope flow velocity.These findings indicate that the aboveground components of vegetation are more effective in reducing flow velocity and increasing hydraulic roughness than the belowground parts.Conversely, vegetation has also been reported to increase slope flow velocity.For example, Cen et al. [17] observed that the average flow velocity initially rose with an increase in the amount of vegetation before an eventual decline.Vegetation also alters surface roughness and affects the water flow [18]: the larger the rough elements, the more concentrated the surface runoff [19], which can also increase the flow velocity.Given the complex relationship between vegetation cover and flow velocity, and the contradicting results in the literature, it is necessary to further study this relationship to refine the theory of slope flow.
Vegetation characteristics such as the shape, stiffness, degree of submergence, and arrangement influence the hydrodynamic properties of slope flow [20].The flow patterns at different positions on the slope are unique, primarily depending on the physical characteristics of the slope [21].Cen et al. [22] studied the hydrodynamic properties of slope flow under combinations of flexible and rigid vegetation and found that different types of vegetation can differentially affect flow resistance.Serio et al. [20] also studied the hydraulic characteristics of flowing rivers under the influence of rigid and flexible vegetation and found that the density and rigidity of vegetation affected the spatiotemporal distribution of the average flow velocity to different degrees.Wu et al. [23] found that water depth increased with the diameter of rigid vegetation, and the flow velocity behind vegetation was lower than that inside its gaps.Teng et al. [24] reported that the distribution of flow velocity under different arrangements of rigid vegetation was lower closer to the channel bottom and in the inner layers of vegetation, with flow velocities along the water depth displaying S-shaped and inverse S-shaped profiles.In natural environments, the growth direction of vegetation is influenced by gravitropism and phototropism.Gravitropism causes roots to grow toward the Earth's center, while stems typically grow vertically upward, perpendicular to the ground plane [25,26].Phototropism directs plants to grow toward light sources, often aligning with vertical growth relative to the horizon.Most plants exhibit this growth pattern in nature [27].However, under certain conditions, such as in high-density vegetation areas or low-light environments, plants may grow perpendicular to the slope to access more light or reduce competition.Currently, most vegetation modeling studies assume that vegetation grows perpendicular to the slope (BS) [22,28,29], overlooking the case where plants grow perpendicular to the horizontal plane (BH).Given the lack of research on the impact of different vegetation growth directions on mean flow velocity, it is essential to differentiate and quantify the effects of these two vegetation types on slope flow velocity.
The prediction of mean velocity has been a focal point of research [30] and it is essential to estimate this accurately in order to forecast runoff and sediment output [12].Nearing et al. [31] hypothesized that mean velocity could be represented by a power function of discharge, while other scholars argue that the prediction of mean velocity should consider coverage in situations with cover [32].Additionally, the Manning formula is widely used to predict slope flow velocity because it includes water depth and roughness in its estimation [33]; however, this formula assumes many constraints, such as a uniform flow, and therefore additional factors must be considered when modifying this formula for practical applications.For example, Cen et al. [17] combined the principle of equal roughness adjustment with the Manning formula to establish a flow velocity prediction model, while Fu et al. [34] reported that the Manning formula did not accurately predict mean velocity under certain vegetation coverage conditions.Overall, further research is needed to develop formulae that predict the mean velocity of slope flow under vegetation.
In summary, the impact of slope vegetation on mean flow velocity is not yet fully understood.Differences in vegetation types used across studies and variations in hydraulic conditions among different experiments contribute to discrepancies in the observed flow velocity patterns.Additionally, most studies have focused solely on flow velocity under rigid vegetation perpendicular to the slope (BS), neglecting the flow velocity patterns under vegetation perpendicular to the horizontal plane (BH).Therefore, it is necessary to explore whether there are similarities or differences in flow velocity patterns under these two vegetation growth orientations.Furthermore, there is a need to refine velocity prediction models and test their applicability under rigid vegetation cover.Based on this, the objectives of this study are to (1) investigate the patterns of average flow velocity under different slopes and flow conditions; (2) analyze the effects of rigid vegetation coverage and the two growth orientations (BS and BH) on mean flow velocity, and understand the mechanisms by which flow velocity differences respond to hydraulic conditions; and (3) develop a new mean flow velocity prediction model based on the Manning equation and compare the accuracy and applicability of different velocity prediction models.Overall, this study contributes to the refinement of slope flow theory and provides data support for soil and water conservation efforts.

Experimental Design
The experimental plot consisted of an acrylic tank with an adjustable slope ranging from 0 • to 30 • , 3.8 m in length, 0.3 m in width, and 0.2 m in depth (Figure 1).A grid-shaped flow stabilizer was installed at the inlet of the tank to maintain the stability of the inflow water.The water supply used for the experiment came from a rectangular tank equipped with a low-flow water pump (QDX100-8-3.5,SRCH, Shanghai, China).The experimental flow was regulated using a glass rotor flow meter (LZB-40, SKJYLEAN, Suzhou, China) designed for a range of 250-2500 L/h.To achieve a uniform roughness scale on the bed, the bottom of the tank was lined with 40-mesh water sandpaper with a 0.38 mm roughness (k sb ).Rigid vegetation, which is a common type of vegetation used in experiments, is typically simulated using cylinders [31].In this experiment, rigid vegetation was simulated using 2 cm diameter PVC pipes, with the vegetation-covered section of the tank extending 2 m with a 20 cm section reserved at the front as a flow stabilization area (Figure 2b).Experiments were conducted using 95, 187, and 286 rigid vegetation stems, as shown in Figure 2a.The rigid vegetation was arranged in a grid pattern, where the front and back rows formed a regular triangular structure.The coverage was defined as the ratio of the total stem cross-sectional area to the area of the covered region, calculated as follows: Equation ( 1) calculates vegetation coverage perpendicular to the slope direction (BS), where Cr is the vegetation coverage, N is the number of rigid vegetation cells in the test area, D is the diameter of the PVC pipe (0.02 m), B is the width of the experimental flume (0.3 m), and L is the length of the area covered by vegetation (2.0 m).For vegetation perpendicular to the horizontal plane (BH), coverage was calculated using the formula for the area of an ellipse (Equation (2)), where θ is the angle of the slope.
Slopes θ were set at 5 • , 10 • , 15 • , and 20 • , which resulted in slope ratios S of 0.0872, 0.1737, 0.2588, and 0.3420, respectively.The designed discharge rates were 5, 10, 20, 30, and 40 L•min −1 , which corresponded to per-width flow rates of 0.000278, 0.000556, 0.001111, 0.001667, and 0.002222 m 3 •m −1 •s −1 , respectively, with five levels of treatment.Three longitudinal observation sections were set along the flume from top to bottom, at 0-0.75 m, 0.75-1.5 m, and 1.25-2.0m, respectively.Water depth h was measured using a water-level gauge (model: SCM60, The Leader, Shanghai, China; accuracy ± 0.01 mm).In each of the three sections, measurements were taken four times, totaling twelve measurements, which were averaged to obtain the final results.If the standard deviation of the results was >5%, re-testing was necessary.

Hydrodynamic Parameters
The formula for calculating mean velocity u is as follows: where Q is the design flow rate (m 3 •s −1 ), h is the average water depth of the cross-section (m), and b 1 is the effective width of the water flow in the experiment (m).Equation ( 4) defines b 1 as follows: Equation ( 5) was used for calculating the Reynolds number Re, as follows: where R is the hydraulic radius (m) and ν 0 is the kinematic viscosity (m 2 •s −1 ), which was calculated with Equation ( 5), as follows: where t is the reading from the thermometer near the flow stabilizer ( • C).Equation ( 7) was used for calculating the Darcy-Weisbach resistance coefficient f, as follows: Equation ( 8) was used to calculate the hydraulic radius R: where A is the cross-sectional area (m 2 ) through which water passes, which is calculated as follows: P is the wetted perimeter (m) and is calculated as follows: Equation ( 10) is the calculation formula for the wetted perimeter of vegetation oriented perpendicularly to the slope (BS), and Equation ( 11) is for vegetation oriented perpendicularly to the horizontal plane (BH).

Data Processing Model Accuracy Assessment Indicators
The accuracy of the mean velocity prediction model was assessed using the adjusted coefficient of determination, Adj.R 2 , as well as the Nash-Sutcliffe efficiency coefficient, NSE; the standard error of the coefficients, SE; and the relative root-mean-square error, RRMSE.The formulae for these three evaluation indicators are as follows: where O i is the observed value, O is the mean of the observed values, P i is the simulated value, P is the mean of the simulated values, σ 2 is the variance of the residuals, n is the number of samples, and k is the number of parameters.

Statistical Analyses
One-way analysis of variance (ANOVA) tests were conducted on the mean velocities under different flow rates, slopes, coverage levels, and types of coverage using SPSS 26 software (IBM SPSS Statistics, Armonk, NY, USA) with a significance level of p < 0.05.The nonlinear regression fitting function in SPSS 26 was used to construct a mean velocity prediction model.R Studio 4.2.0 (R Core Team, Vienna, Austria) and Origin 2023 (OriginLab Corporation, Northampton, MA, USA) were used to create graphics for visualization purposes.

The Impact of Flow Rate and Slope on Mean Velocity
As the flow rate increased, the distribution range of the mean velocity increased, with flow rates of 5, 10, 20, 30, and 40 L•min −1 showing mean velocity distribution ranges of 0.117 to 0.245 m•s −1 , 0.140 to 0.314 m•s −1 , 0.178 to 0.393 m•s −1 , 0.202 to 0.463 m•s −1 , and 0.213 to 0.582 m•s −1 for BS; and 0.158 to 0.307 m•s −1 , 0.183 to 0.419 m•s −1 , 0.204 to 0.482 m•s −1 , 0.222 to 0.463 m•s −1 , and 0.238 to 0.582 m•s −1 for BH, respectively (Figure 3).These results indicated that the flow rate and slope significantly affected the mean velocity distribution range (p < 0.05).A comparison of the two vegetation conditions revealed that at lower flow rates (5, 10, and 20 L•min −1 ), the range of mean velocity under BS was smaller than that under BH.At higher flow rates (30 and 40 L•min −1 ), the mean velocity distribution range under BS conditions was similar to that under BH conditions.As the slope increased, the distribution range of the mean velocity also increased.At slopes of 5 • , 10 • , 15 • , and 20 • , the ranges were 0.117 to 0.475 m•s −1 , 0.179 to 0.498 m•s −1 , 0.167 to 0.557 m•s −1 , and 0.166 to 0.582 m•s −1 for BS conditions; and 0.158 to 0.475 m•s −1 , 0.179 to 0.498 m•s −1 , 0.186 to 0.557 m•s −1 , and 0.166 to 0.582 m•s −1 for BH conditions, respectively, indicating that the range of the mean velocity distributions was similar under both BS and BH conditions.

The Impact of Rigid Vegetation on Mean Velocity
Under rigid vegetation cover conditions, the mean velocity of the slope flow was influenced by a combination of vegetation coverage, flow rate, and slope.At low flow rates, the mean velocity initially increased and then decreased with changes in vegetative coverage (Figure 4).For instance, at a slope of 5° and flow rate of 5 L•min −1 , the mean velocity increased to 0.158, 0.177, 0.145, and 0.117 m•s −1 with increasing coverage (N = 0-286) under BS conditions, while at a slope of 10° and flow rate of 10 L•min −1 , the mean velocity rates were 0.274, 0.312, 0.301, and 0.279 m•s −1 with increasing coverage (N = 0-286) under BH conditions.At a slope of 5° and a flow rate of 30 L•min −1 , the mean velocities with increasing BH coverage conditions were 0.403, 0.274, 0.245, and 0.222 m•s −1 , indicating that at lower slopes, as the flow rate increased, the mean velocity gradually decreased.At medium-to-high slopes (10-20°), the mean velocity gradually decreased and then increased as the flow rate increased, and was most pronounced at a flow rate of 20 L•min −1 .Indeed,

The Impact of Rigid Vegetation on Mean Velocity
Under rigid vegetation cover conditions, the mean velocity of the slope flow was influenced by a combination of vegetation coverage, flow rate, and slope.At low flow rates, the mean velocity initially increased and then decreased with changes in vegetative coverage (Figure 4).For instance, at a slope of 5 • and flow rate of 5 L•min −1 , the mean velocity increased to 0.158, 0.177, 0.145, and 0.117 m•s −1 with increasing coverage (N = 0-286) under BS conditions, while at a slope of 10 • and flow rate of 10 L•min −1 , the mean velocity rates were 0.274, 0.312, 0.301, and 0.279 m•s −1 with increasing coverage (N = 0-286) under BH conditions.At a slope of 5 • and a flow rate of 30 L•min −1 , the mean velocities with increasing BH coverage conditions were 0.403, 0.274, 0.245, and 0.222 m•s −1 , indicating that at lower slopes, as the flow rate increased, the mean velocity gradually decreased.At medium-to-high slopes (10-20 • ), the mean velocity gradually decreased and then increased as the flow rate increased, and was most pronounced at a flow rate of 20 L•min −1 .Indeed, at a slope of 20 • and a flow rate of 20 L•min −1 , the mean velocities with increasing coverage were 0.393, 0.348, 0.288, and 0.358 m•s −1 under BS conditions.At a flow rate of 5 L•min −1 with slopes of 15 • and 20 • , the rates of change in the range of mean velocities were 54.9% and 85.3%, respectively, indicating that as the flow rate increased, the differences in velocity between different coverage levels decreased and the rate of velocity stabilized.However, at a flow rate of 40 L•min −1 with slopes of 15 • and 20 • , the rates of change in the range of mean velocities were only 23.9% and 20.9%, respectively, suggesting that the trend in the mean velocity change at high discharge was gentler than that at low discharge.In addition, the influence of rigid vegetation on velocity at low discharge was more apparent.
at a slope of 20° and a flow rate of 20 L•min −1 , the mean velocities with increasing coverage were 0.393, 0.348, 0.288, and 0.358 m•s −1 under BS conditions.At a flow rate of 5 L•min −1 with slopes of 15° and 20°, the rates of change in the range of mean velocities were 54.9% and 85.3%, respectively, indicating that as the flow rate increased, the differences in velocity between different coverage levels decreased and the rate of velocity stabilized.However, at a flow rate of 40 L•min −1 with slopes of 15° and 20°, the rates of change in the range of mean velocities were only 23.9% and 20.9%, respectively, suggesting that the in the mean velocity change at high discharge was gentler than that at low discharge.In addition, the influence of rigid vegetation on velocity at low discharge was more apparent.a, b, c, d, e, and f) above the bar plots indicate significant differences at the 0.05 significance level: groups with the same letter are not significantly different from each other, while groups with different letters are significantly different.
The ANOVA tests showed that under different slope and flow conditions, the coverage of rigid vegetation and the vegetation type had a significant effect on the mean flow velocity (p < 0.05).Post hoc tests were used to compare the flow velocities between two types of rigid vegetation (BS and BH) under the same hydraulic conditions, and specific results indicated by "a", "b", "c", "d", "e", and "f" labels are shown in Figure 4.For  a, b, c, d, e,  and f) above the bar plots indicate significant differences at the 0.05 significance level: groups with the same letter are not significantly different from each other, while groups with different letters are significantly different.
The ANOVA tests showed that under different slope and flow conditions, the coverage of rigid vegetation and the vegetation type had a significant effect on the mean flow velocity (p < 0.05).Post hoc tests were used to compare the flow velocities between two types of rigid vegetation (BS and BH) under the same hydraulic conditions, and specific results indicated by "a", "b", "c", "d", "e", and "f" labels are shown in Figure 4.For example, under the conditions of a 5 • slope and a flow rate of 10 L•min −1 , the mean flow velocity for BH was significantly higher than that for BS at all three coverage levels (p < 0.05).However, under the conditions of a 20 • slope and a flow rate of 5 L•min −1 , although the mean flow velocity for BH was also higher than that for BS at all three coverage levels, the Water 2024, 16, 2240 9 of 17 difference was not statistically significant where N = 187 (p > 0.05).To compare the effects of the two types of rigid vegetation on the mean velocity, we used a velocity ratio defined as the velocity under BS conditions (u BS ) divided by that under BH conditions (u BH ).A velocity ratio greater than 1 indicated that the velocity under the BS conditions was greater than that under the BH conditions, and vice versa.Correlation tests revealed significant relationships among the velocity ratio and the Reynolds number ratio, resistance coefficient ratio, and flow rate (p < 0.05), with Pearson correlation coefficients of 0.83, -0.95, and 0.61, respectively.As the Reynolds number ratio and flow rate increased, the velocity ratio increased (Figure 5), whereas an increase in the resistance ratio led to a decreased velocity ratio.Figure 5a and 5b, respectively, show the relationships between the velocity ratio and the Reynolds number ratio and between the velocity ratio and the flow rate.The R 2 values of the linear equations are 0.70 and 0.37, respectively.Overall, these results indicate that not only does the difference in rigid vegetation coverage affect the average flow velocity, but it is also affected by the direction of rigid vegetation growth (BS and BH).
example, under the conditions of a 5° slope and a flow rate of 10 L•min −1 , the mean flow velocity for BH was significantly higher than that for BS at all three coverage levels (p < 0.05).However, under the conditions of a 20° slope and a flow rate of 5 L•min −1 , although the mean flow velocity for BH was also higher than that for BS at all three coverage levels, the difference was not statistically significant where N = 187 (p > 0.05).To compare the effects of the two types of rigid vegetation on the mean velocity, we used a velocity ratio defined as the velocity under BS conditions (uBS) divided by that under BH conditions (uBH).A velocity ratio greater than 1 indicated that the velocity under the BS conditions was greater than that under the BH conditions, and vice versa.Correlation tests revealed significant relationships among the velocity ratio and the Reynolds number ratio, resistance coefficient ratio, and flow rate (p < 0.05), with Pearson correlation coefficients of 0.83, -0.95, and 0.61, respectively.As the Reynolds number ratio and flow rate increased, the velocity ratio increased (Figure 5), whereas an increase in the resistance ratio led to a decreased velocity ratio.Figure 5a and 5b, respectively, show the relationships between the velocity ratio and the Reynolds number ratio and between the velocity ratio and the flow rate.The R 2 values of the linear equations are 0.70 and 0.37, respectively.Overall, these results indicate that not only does the difference in rigid vegetation coverage affect the average flow velocity, but it is also affected by the direction of rigid vegetation growth (BS and BH).

Mean Velocity Prediction Model
The prediction of mean velocity has been a major area of research among scholars [35,36], with the mean velocity being expressed as a power function of slope and flow rate, with vegetation coverage being a determinant of exponential size [37].In this study, a nonlinear regression analysis was performed on the experimental data, and the results of this analysis are presented in Table 1.

Mean Velocity Prediction Model
The prediction of mean velocity has been a major area of research among scholars [35,36], with the mean velocity being expressed as a power function of slope and flow rate, with vegetation coverage being a determinant of exponential size [37].In this study, a nonlinear regression analysis was performed on the experimental data, and the results of this analysis are presented in Table 1.The relationship between the predicted and measured mean flow velocities were calculated using Equation (19), which indicated that the Adj.R 2 under BH conditions was 0.77 (Equation ( 18)), while the Adj.R 2 for the prediction model fitted with all data was 0.80 (Equation ( 19); Figure 6c).Because of the more complex dynamics of overland flow compared to open-channel flow, especially when surface runoff is influenced by vegetation, using only flow and slope for prediction is not ideal [38].For slopes covered with vegetation, proposed theories such as vegetation equivalent resistance coefficients [30] and vegetation-based roughness height [39] integrate the effects of vegetation on overland flow with the underlying surface impacts to enhance the practicality of hydraulic models.Rigid vegetation acts as a roughness element on slopes and its impact can be quantified using the resistance coefficient formula for rough areas, as follows [39]: Equations ( 7) and ( 20) calculate the resistance coefficients; thus, by combining and rearranging them, we can derive the vegetation equivalent roughness height k sv (mm): Because the substrate used in this experiment was 40-mesh sandpaper, the roughness height of the slope outside the vegetated area remained constant (k sb = 0.38 mm).Adding the rigid vegetation equivalent roughness height k sv to that of the non-vegetated area k sb yielded the equivalent roughness height k se (mm) of the slope under stiff vegetation cover.The resulting formula is as follows: The equivalent roughness height, k se , of the slope under stiff vegetation was correlated with the slope, flow, and coverage.For easier field observation and application, a nonlinear fit using the equivalent roughness height, k se ; the flow rate, Q; the energy slope, S; and the coverage, Cr, was performed using the following equation: k se = 132.154Q 0.718 S −0.753 e 6.984CrThe Adj.R 2 value of Equation ( 23) was 0.78, indicating that the model structure was reasonable.Manning's resistance coefficient, n, depends on the vegetation coverage and soil type [40].Because this experiment was a fixed-bed test, the Manning's coefficient was primarily influenced by the coverage of rigid vegetation.The formula used for calculating Manning's coefficient, n v , under vegetation cover is as follows: where m 1 and m 2 are the exponents of the hydraulic radius R and energy slope S, respectively.Compared with the basic Manning formula, the use of the hydraulic radius R was more appropriate because of the changes in water depth, h, caused by rigid vegetation cover.Wu et al. [41] found that Manning's coefficient decreases with increasing water depth and increases with increasing slope roughness.This is similar to the findings of Yang et al. [42], who reported that the Manning's coefficient, n b , for slopes without vegetation cover can be expressed as a function of the slope roughness height, k sb , and water depth, h, using the following specific function form: where m 3 and m 4 are the exponents of the slope roughness height k sb and water depth h, respectively.Similarly, by replacing the slope roughness height k sb in Equation ( 25) with the equivalent roughness height k se under a rigid vegetation cover, the following formula for the equivalent vegetation Manning coefficient n ve was derived, as follows: By combining Equations ( 24)-( 26), the model expression for the mean flow velocity u was constructed: where a 1 , a 2 , a 3 , a 4 , and a 5 are the constants.The equivalent roughness height, k se , of the slope in Equation ( 27) can be calculated using Equation ( 23), and the flow velocity prediction model was thus obtained after nonlinear fitting of the experimental data (Table 2).
where  3 and  4 are the exponents of the slope roughness height   and water depth h, respectively.Similarly, by replacing the slope roughness height   in Equation (25) with the equivalent roughness height   under a rigid vegetation cover, the following formula for the equivalent vegetation Manning coefficient   was derived, as follows: By combining Equations ( 24)-( 26), the model expression for the mean flow velocity u was constructed: where  1 ,  2 ,  3 ,  4 , and  5 are the constants.The equivalent roughness height,   , of the slope in Equation ( 27) can be calculated using Equation ( 23), and the flow velocity prediction model was thus obtained after nonlinear fitting of the experimental data (Table 2).

The Impact of Flow Rate and Slope on Mean Velocity
The mean flow velocity, u, a crucial hydraulic characteristic in soil erosion modeling, is influenced by the flow rate [43], slope, and surface conditions.In this study, slope and flow rate increased the mean flow velocity, which suggested that the gravitational potential energy of the water flow rises as the slope becomes steeper, resulting in higher flow velocities on the slope.This outcome aligns with the findings of Zhang et al. [44] and Liu et al. [32].Whether on a smooth bare slope or a slope with surface cover, u remains highly correlated with slope and flow rate [6].For instance, Liu et al. [45] found that on gravel-covered slopes, the correlation coefficients between the mean flow velocity, flow rate, and slope were 0.716 and 0.674, respectively.Similarly, the slope factor and rate of addition could be used to estimate u on slopes with added straw (R 2 = 0.91) [46].However, Ali et al. [47] found that under mobile bed conditions, the slope did not significantly affect u, which may be due to the different substrate conditions.In our study, the substrate and bed conditions were kept constant throughout the experiment to explore how the flow rate and slope affected the morphology and roughness of erosive slopes.Giménez et al. [48] suggested that slope roughness is positively correlated with slope; thus, the greater the slope, the rougher the slope surface.A rougher substrate can offset the kinetic energy converted from the gravitational potential energy [49], thus resulting in no significant change in u.

The Impact of Rigid Vegetation on Mean Velocity
It is commonly believed that vegetation increases surface roughness, and that the energy of the slope flow is mainly consumed by the surface roughness caused by the aboveground parts of vegetation [50].Thus, vegetation reduces the mean flow velocity, u, on slopes.For example, under three different vegetation covers, the flow velocity on slopes decreased by 28-30% [15].Liu et al. [51] added different amounts of straw to slopes to increase surface roughness and found that the flow velocity decreased by 28.44%, 44.09%, and 55.56%.However, in our experiment, there was no monotonic decrease in the association between u and the coverage of rigid vegetation, and u was differentially influenced by the test conditions (Figure 4).At low flow rates and with increases in slope, u fluctuates as follows: it first increases, then decreases, and then finally increases.In comparison, at high flow rates, u fluctuates as follows: it first decreases and then increases in response to an increase in the vegetation cover.The primary reason for the increase in u was that the rigid vegetation changed the cross-sectional area of the flow on the slope, which impacted the spatial distribution of the water [15].Rigid vegetation was the primary roughness element on the slope and reduced the flow velocity due to its obstructive effect, while also accelerating water flow through non-vegetated areas, forming concentrated flows and increased velocity [52].Our findings are supported by those of Wu et al. [23], who observed that the flow velocity behind rigid vegetation was lower than that within vegetation gaps using a three-dimensional laser Doppler velocimeter.Vegetation also affects water depth by increasing the water level upstream and decreasing it downstream of vegetation [53].This backwater-induced local pressure difference increases the local flow velocity and shear stress, increasing the erosion risk near vegetated areas [54].Zhao et al. [28] used PVC pipes to replicate rigid vegetation and showed that the flow velocity increases marginally with greater coverage.At low flow rates, the critical coverage for the increase in the mean flow velocity was low, and this critical coverage increased as the flow rate increased.The main resistance to flow on slopes comes from the particle resistance generated by the slope and the form resistance created by vegetation [55], with particle resistance generally decreasing as the flow rate and slope increase [38].At low flow rates, the slope flow is slower and shallower, and vegetation with low coverage compresses the water flow space, increasing the flow velocity.As the flow rate increases, both the flow velocity and water depth increase, at which time the particle resistance on the slope is gradually replaced by the form resistance of the vegetation.The increased flow velocity owing to the compression of the water cross-sectional area by vegetation at low coverage cannot replace the reduced flow velocity owing to vegetation form resistance [34].Under these conditions, the compression of the water flow space by high-coverage vegetation becomes more pronounced [56], leading to increased flow velocity.The trend in changes in mean flow velocity depends on whether particle resistance or form resistance predominates in the slope flow [57].
For the two different types of rigid vegetation used in our study (BS and BH), the morphological differences resulted in different degrees of obstruction to the water flow, which affected the mean flow velocity.We found that under BS conditions, the average flow velocity (0.331 m•s −1 ) was higher than that under BH conditions (0.311 m•s −1 ).The velocity ratio was positively correlated with both the flow ratio and the Reynolds number ratio, while in contrast, it was found to be negatively correlated with the resistance coefficient ratio.We found that most Reynolds number ratios and velocity ratios were generally less than 1 (Figure 5), indicating that under BS conditions, the rigid vegetation significantly obstructed water flow, inhibited changes in the slope flow state, and lowered the flow velocity.Similarly, an increase in flow suggested greater water depth.Under the BS conditions, vegetation compressed the flow more significantly; therefore, under highflow and BS conditions, the flow velocity gradually exceeded that reported under the BH conditions.These findings are congruent with those of Schoelynck et al. [58], who postulated that vegetation on a slope provides less hindrance to water flow because it is more streamlined.Shan et al. [59] also found that the angle at which vegetation is deflected increases water flow when the water flow contacts the vegetation perpendicularly (e.g., BS)-i.e., it encounters greater obstruction-whereas when the vegetation is angled (e.g., BH), it exerts less drag, making it easier for the flow to bypass [60].In this study, the increase in rigid vegetation coverage altered the connectivity of the water flow.For slope flow under BS conditions, as the number of rigid vegetation N increased from 95 to 286, the coefficient of variation in the mean flow velocity was 30.72%, 35.27%, 29.82%, and 41.61%, respectively.In contrast, for slope flow under BH conditions, the coefficient of variation in the mean flow velocity was 25.68%, 29.82%, and 33.92% as N increased from 95 to 286.In the arrangement of a regular grid, the increase in vegetation coverage made the water flow paths more complex, and the variation in mean flow velocity under different hydraulic conditions became more complicated.Therefore, the coefficient of variation in the flow velocity increased with the increase in coverage.Additionally, the coefficient of variation under BS conditions was consistently higher than that under BH conditions, possibly due to the greater impact of rigid vegetation on water flow in BS conditions.

Evaluation of Mean Flow Velocity Prediction Models
Traditional hydraulic models consider flow velocity as a power function of both the flow rate and slope [32]; however, on vegetated slopes, flow rate and slope are no longer the sole factors affecting flow velocity [61].Equations ( 17)-( 19) reflect the mean flow velocity prediction models fitted using the slope and flow rate as variables.For all data, Figure 6c shows that the data points predicted by Equation ( 19) are close to the 1:1 line, although some data points are more scattered.In comparison, the data points predicted by Equation ( 30) are uniformly distributed near the 1:1 line (Figure 6d), with less dispersion than those predicted by Equation (19).This demonstrates that utilizing the vegetation equivalent roughness height to predict u is rational and yields a highly accurate model.Furthermore, it is unreasonable to predict u solely based on the flow rate and slope on slopes with rigid vegetation; therefore, factors related to vegetation must be considered when predicting mean flow velocities.
To enhance the representativeness of our model results, we examined the flow velocity models from the literature [17,31].Nearing et al. [31] reported experiments that were conducted on soil slopes and their results showed that the model's Adj.R 2 , NSE, RRMSE, and SE were 0.45, −0.01, 0.34, and 0.082, respectively.Our results had a high dispersion of data points near the 1:1 line, which indicated that this model is unsuitable for slopes covered with rigid vegetation (Figure 6a).The model did not consider the slope factor because the roughness of soil slopes increased with slope during erosion, which offset the contribution of the slope to u [48].This is similar to Emmett's [21] study, which concluded through field experiments that slope roughness slows flow velocity and increases water depth, with the increase in water depth being primarily attributed to vegetation and topographic features.Cen et al. [17] reported a model with Adj.R 2 , NSE, RRMSE, and SE values of 0.61, 0.55, 0.22, and 0.082, respectively.This model considered factors such as hydraulic radius, slope, vegetation coverage, and water depth and was more accurate than the previous model [31].However, the equation specifically developed for modeling vegetated slope flow velocities did not assess flow velocities on bare slopes [17].We found that this model was relatively dispersed near the 1:1 line primarily because of differences in the type of simulated vegetation (Figure 6b).Additionally, this model was developed specifically for tufted synthetic grass, which has a stronger disturbance effect on water flow.As a result, the model predicts a higher degree of variability in the average flow velocity compared to Equations ( 19) and (30).The influence of vegetation on the hydrodynamic properties of surface water is affected not only by the above-and underground parts of the vegetation [16] but also by different morphological features of the vegetation [62].As shown in Table 3, Equation (30) demonstrates significant improvements in all metrics, especially in Adj.R 2 , which shows an enhancement ranging from 8.75% to 93.3% compared to other models.This indicates that this model exhibits superior fitting and predictive efficiency.Additionally, the lower RRMSE (0.12) and SE (0.036) values suggest that this model has smaller prediction errors, providing better accuracy and stability.

Conclusions
We conducted indoor simulation experiments to investigate the mechanisms by which rigid vegetation coverage and growth direction affect the mean flow velocity under conditions of four slope gradients, five flow rates, four levels of vegetation coverage, and two vegetation growth directions.First, we analyzed the impact of slope and flow rate on mean flow velocity and found that both significantly affected mean flow velocity (p < 0.05), with the mean flow velocity increasing as both the slope and flow rate increased.Second, we analyzed the impact of rigid vegetation on the mean flow velocity, observing both similarities and differences with previous studies.Our results indicate that rigid vegetation coverage not only reduces the mean flow velocity but can also increase it under certain conditions, primarily owing to the compression of the water flow path by the rigid vegetation.Additionally, the mean flow velocity under BH conditions was generally higher than that under BS conditions, and the velocity ratio was significantly correlated with the Reynolds number ratio, resistance coefficient ratio, and flow rate (p < 0.05).Finally, we established a model to predict the mean flow velocity by combining vegetation equivalent roughness and the Manning formula.The results showed that this model could accurately predict the average flow velocity (Adj.R 2 = 0.87, NSE = 0.86, RRMSE = 0.12, and SE = 0.036).
On natural slopes, the conditions of the underlying surface change due to water flow erosion, and the movement of surface soil particles alters surface roughness.This study was conducted in a fixed-bed flume, which did not account for the effects of bed erosion on flow velocity on slopes.Furthermore, in the field, vegetation growth can exhibit more complex characteristics, such as combinations of different stem diameters, different vegetation cover patterns, and combinations of flexible and rigid vegetation.These vegetation characteristics can all impact the hydrodynamic properties of slope flows.Additionally, the response of flow velocity to slope erosion is also crucial.Therefore, future research should focus on establishing the relationship between slope flow velocity and soil erosion with different vegetation characteristics, which will further contribute to the improvement of soil erosion prediction models.

Figure 2 .Figure 1 . 17 Figure 1 .
Figure 2. Schematic illustration of rigid vegetation arrangements and configurations: (a) The arrangement of rigid vegetation under varying coverage densities, denoted by N = 95, N = 187, and N = 286, with respective measurements of plant spacing.(b) Two configurations of rigid vegetation

Figure 2 .
Figure 2. Schematic illustration of rigid vegetation arrangements and configurations: (a) The ar rangement of rigid vegetation under varying coverage densities, denoted by N = 95, N = 187, and N = 286, with respective measurements of plant spacing.(b) Two configurations of rigid vegetation

Figure 2 .
Figure 2. Schematic illustration of rigid vegetation arrangements and configurations: (a) The arrangement of rigid vegetation under varying coverage densities, denoted by N = 95, N = 187, and N = 286, with respective measurements of plant spacing.(b) Two configurations of rigid vegetation structures: (I) with pipes positioned perpendicular to the slope surface and (II) with pipes positioned perpendicular to the horizontal plane of the slope.
m•s −1 for BS conditions; and 0.158 to 0.475 m•s −1 , 0.179 to 0.498 m•s −1 , 0.186 to 0.557 m•s −1 , and 0.166 to 0.582 m•s −1 for BH conditions, respectively, indicating that the range of the mean velocity distributions was similar under both BS and BH conditions.

Figure 3 .
Figure 3. Distribution of mean flow velocities under different flow rate and slope conditions.(a) Variation in flow velocity as driven by discharge on a sloped surface; (b) variation in flow velocity as driven by discharge on a horizontal plane; (c) variation in flow velocity as driven by slope angle on a sloped surface; (d) variation in flow velocity as driven by slope angle on a horizontal plane.Note: Different letters (a, b, c, and d) above the box plots indicate significant differences at the 0.05 significance level: groups with the same letter are not significantly different from each other, while groups with different letters are significantly different.

Figure 3 .
Figure 3. Distribution of mean flow velocities under different flow rate and slope conditions.(a) Variation in flow velocity as driven by discharge on a sloped surface; (b) variation in flow velocity as driven by discharge on a horizontal plane; (c) variation in flow velocity as driven by slope angle on a sloped surface; (d) variation in flow velocity as driven by slope angle on a horizontal plane.Note: Different letters (a, b, c, and d) above the box plots indicate significant differences at the 0.05 significance level: groups with the same letter are not significantly different from each other, while groups with different letters are significantly different.

Figure 4 .
Figure 4. Variations in mean flow velocity across different degrees of vegetative coverage: Row 1: a slope angle of 5° across flow rates from 5 to 40 L•min −1 ; Row 2: a slope angle of 10° across flow rates from 5 to 40 L•min −1 ; Row 3: a slope angle of 15° across flow rates from 5 to 40 L•min −1 ; Row 4: a slope angle of 20° across flow rates from 5 to 40 L•min −1 .Note: Different letters (a, b, c, d, e, and f) above the bar plots indicate significant differences at the 0.05 significance level: groups with the same letter are not significantly different from each other, while groups with different letters are significantly different.

Figure 4 .
Figure 4. Variations in mean flow velocity across different degrees of vegetative coverage: Row 1: a slope angle of 5 • across flow rates from 5 to 40 L•min −1 ; Row 2: a slope angle of 10 • across flow rates from 5 to 40 L•min −1 ; Row 3: a slope angle of 15 • across flow rates from 5 to 40 L•min −1 ; Row 4: a slope angle of 20 • across flow rates from 5 to 40 L•min −1 .Note: Different letters (a, b, c, d, e, and f) above the bar plots indicate significant differences at the 0.05 significance level: groups with the same letter are not significantly different from each other, while groups with different letters are significantly different.

Figure 5 .
Figure 5. Relationship of the velocity ratio (a) with the Reynolds number ratio and (b) with the flow rate.

Figure 5 .
Figure 5. Relationship of the velocity ratio (a) with the Reynolds number ratio and (b) with the flow rate.

Figure 6 .
Figure 6.Comparison of predicted and measured values across different mean flow velocity prediction models: (a) based on the model developed by Nearing et al. [31]; (b) based on the model developed by Cen et al. [17]; (c,d) based on equations developed in this study, specifically Equations (19) and(30).

Figure 6 .
Figure 6.Comparison of predicted and measured values across different mean flow velocity prediction models: (a) based on the model developed by Nearing et al. [31]; (b) based on the model developed by Cen et al. [17]; (c,d) based on equations developed in this study, specifically Equations (19) and (30).

Table 1 .
Mean flow velocity prediction model for different flow rate and slope conditions.

Table 1 .
Mean flow velocity prediction model for different flow rate and slope conditions.NSE, Nash-Sutcliffe efficiency; RRMSE, relative root-mean-square error; SE, standard error of the coefficients; u, mean velocity; S, slope ratio; Q, flow rate (m 3 •s −1 ).

Table 2 .
The mean velocity prediction model using the equivalent roughness height, k se .
Notes: BS, perpendicular to the slope surface; BH, perpendicular to the horizontal plane; ALL, all data combined; Adj.R 2 , adjusted R square; NSE, Nash-Sutcliffe efficiency; RRMSE, relative root-mean-square error; SE, standard error of the coefficients; u, mean velocity; S, slope ratio; R, hydraulic radius (m); h, average water depth (m); k se , equivalent roughness height (mm).

Table 3 .
Comparison of different mean velocity prediction models.