Overland Flow Resistance Law under Sparse Stem Vegetation Coverage

To explore the characteristics of overland flow resistance under the condition of sparse vegetative stem coverage and improve the basic theoretical research of overland flow, the resistance characteristics of overland flow were systematically investigated under four slope gradients (S), seven flow discharges (Q), and six degrees of vegetation coverage (Cr). The results show that the Manning roughness coefficient (n) changes with the ratio of water depth to vegetation height (h/hv) while the Reynolds number (Re), Froude number (Fr), and slope (S) are closely related to vegetation coverage. Meanwhile, h/hv, Re, and Cr have strong positive correlations with n, while Fr and S have strong negative correlations with n. Through data regression analysis, a power function relationship between n and hydraulic parameters was observed and sensitivity analysis was performed. It was concluded that the relationship between n and h/hv, Re, Cr, Q, and S shows the same law; in particular, for sparse stem vegetation coverage, Cr is the dominant factor affecting overland flow resistance under zero slope condition, while Cr is no longer the first dominant factor affecting overland flow resistance under non-zero slope condition. In the relationship between n and Fr, Cr has the least effect on overland flow resistance. This indicates that when Manning roughness coefficient is correlated with different hydraulic parameters, the same vegetation coverage has different effects on overland flow resistance. Therefore, it is necessary to study overland flow resistance under the condition of sparse stalk vegetation coverage.


Introduction
Overland flow is one of the causes of soil erosion and water loss. It is also the main driving factor for the migration of pollutants, sediment production, and topographic evolution [1][2][3][4]. Overland flow characteristics are mainly affected by surface coverage, slope, and flow discharge [5,6]. Vegetation plays an important role in the benign cycle of the ecological environment and has become an indispensable part of the ecological environment, so it is for the ecological slope [7][8][9][10][11][12]. Vegetation modifies the mixing effect of turbulence and dissipates the energy of overland flow, affects sediment transport and settlement, and then reduces soil erosion on slopes. Meanwhile, only when vegetation cover reaches a certain degree can the vegetation effect be prominent; if the vegetation cover is too large, the effect of vegetation will not only not increase, but the surface soil moisture and nutrients will decrease greatly due to the sharp increase of vegetation consumption, leading to soil drying and increasing soil burden. Therefore, it is necessary to clarify the change characteristics of the vegetation slope runoff resistance coefficient, reveal the hydrodynamic characteristics and dynamic erosion processes of slope runoff, and provide a scientific basis for the treatment of vegetation as it relates to soil erosion.
The presence of vegetation has significant effects on the resistance of land to overland flow [8,[13][14][15][16][17][18][19]. At present, a great amount of research has been conducted on the influence of vegetation on the Manning roughness coefficient. Previous studies have mainly focused on the effects of vegetation type [20][21][22], vegetation flexibility [23], vegetation relative submergence height [24,25], and vegetation stem diameter [19,26] on the parameters of flow resistance. However, the effect of vegetation coverage on flow resistance and characteristics has also been studied, but the research premise is mostly based on larger vegetation coverage [27][28][29][30][31].
In fact, plant density plays an important role in influencing flow resistance [32]. With the continuous strengthening of research on the characteristics of resistance to overland flow, the effects of vegetation cover on flow resistance have been given an increasing amount of attention recently. Zhao et al. [19] considered that the Reynolds number (Re) for vegetated slopes was much higher than it was for bare slopes due to the influence of vegetation stems on the effective flow width and unit flow discharge. Pan and Shangguan [33] and Wu et al. [34] found that on natural slope, the vegetation stems could increase the infiltration and reduce the runoff and then weaken the soil erosion. Therefore, the Re decreased as vegetative cover increased under simulated rainfall on mobile beds. In addition, the Froude number (Fr) decreased and Manning roughness coefficient (n) increased with increased vegetative cover [19,33,35]. Nepf [36] believed that even with similar arrays and Reynolds numbers, the coefficient of flow resistance may decrease with increasing vegetation density. The experiment of Liu and Zeng [37] revealed that the Drag coefficient logarithmically increases with the vegetation density under supercritical flow, while for subcritical flow, the Drag coefficient logarithmically decreases with the vegetation density. This has also been presented by Wu et al. [38] based on limited experimental results, whom concluded that the roughness coefficients for sub-and supercritical flows have opposite trends of variation. However, Stoesser et al. [39] believed that the Drag coefficient increased with increasing vegetation density.
Therefore, the effects of the density of vegetation cover on the flow resistance coefficient and hydraulic parameters remain unclear. In addition, as far as the study of vegetation cover is concerned, many researchers of studied the types of vegetation cover [8,19,40,41], with most research concentrating on grass cover [42][43][44][45][46][47][48]. However, some studies have shown that many areas are still sparsely vegetated [49,50], especially in the wind and water erosion transitional belt of the Loess Plateau, where the vegetation is relatively vulnerable to erosion and most areas have <10% vegetation cover [40]. In regions with sparse vegetation cover, the vegetation type mainly consists of stalk and woody vegetation such as shrubs and trees [51]. However, relatively few studies have addressed the characteristics of overland flow resistance in areas covered by sparsely vegetated areas with stem-type and woody vegetation, which is also the focus of this article.
Traditional studies of vegetative cover have considered the combined roles of leaves and stems. However, the presence of vegetation stems will affect the longitudinal turbulent mixing effect and energy dissipation of flow [52][53][54]. In addition, Zhou and Shangguan [55,56] conducted research on loess erosion through simulated rainfall experiments, and the results showed that vegetation stalks have a greater contribution to reducing runoff. Pan et al. [33] showed that vegetation stems can significantly reduce the sediment yield of overland flow. Because water and sediments mainly pass between and directly act on vegetation stems, vegetative stems behave as the dominant roughness element in overland flow [19]. Thus, this paper simulated vegetation water discharge and used fluid mechanics theory to clarify the characteristics of overland flow resistance under the condition of sparse stem cover within a certain range of flow and slope; quantify the intrinsic relationship between the Manning roughness coefficient and the characteristic hydraulic parameters of the overland flow (h/hv, Cr, S, Re, Fr, Q). Especially in the fragile ecological environment areas where stem vegetation is mainly distributed, the vegetation coverage on the slope is highly sparse. However, few studies have been conducted on the effect of sparse vegetation coverage on the characteristics of overland flow. Therefore, it is necessary to clarify the effect of sparse vegetation coverage on the characteristics of overland flow resistance. It provides a more comprehensive theoretical framework for future studies on the effect of vegetation coverage on overland flow. The results help to simplify the calculation of hydraulic parameters used to predict sediment transport capacity, so as to better understand the role of plants in soil and water conservation.

Experimental Structure
The overland flow study reported here mainly used the basic theory of flow in an open channel [57,58]. Many studies on the effects of vegetation on a slope during water erosion have used variable slopes involving rectangular flume experiments. To enhance the experimental effect and improve the experimental efficiency, single-plant vegetation has sometimes been generalized into an effective water-blocking cylinder in these experiments [59,60]. The structure of the test tank used in the present study is shown in Figure 1. The following design points were determined based on the purpose of the test: (1) The test area was a polymethyl methacrylate flume. Four slopes (0, 1.0%, 2.0%, and 3.0%) were selected. The test flume was 5.0 m long with a width of 0.4 m and a height of 0.3 m; this flume was divided into three parts: an upstream flat water section (1 m), a vegetation coverage section (3 m), and a downstream water measuring section (1 m).
(2) Plastic cylindrical rods were used to simulate the stems of rigid vegetation. The height of the simulated vegetation (hv) was set at 0.15 m with a stem diameter (d) of 0.004 m to match the scale of the model. Simulated vegetation was arranged in a regular pattern; the distances between two rows of cylinders were 0.04 m, 0.06 m, 0.08 m, 0.10 m, 0.12 m, and 0.14 m, respectively. The ratios of stem diameter to adjacent plant spacing (d/l) were 10.00%, 6.67%, 5.00%, 4.00%, 3.33%, and 2.86%, respectively. The total number of cylindrical rods in each entire vegetation segment was 688, 420, 344, 239, 172, and 134, respectively. The stem cover (Cr) values assessed in this study were approximately 0.72%, 0.37%, 0.20%, 0.13%, 0.11%, and 0.07%, respectively (Figure 2), the corresponding stem cover (Cr) values can be calculated by using Equation (1) in the "Data Analysis" section.
(3) A quantitative analysis of the Manning roughness coefficient and hydraulic parameters was conducted to systematically study the effects of the change of stem vegetation on the internal relationship between the coefficient and parameters. The final design of seven discharge flows (Q) were 0.0313 m 3 /min, 0.1205 m 3 /min, 0.2190 m 3 /min, 0.3019 m 3 /min, 0.4241 m 3 /min, 0.5203 m 3 /min, and 0.6319 m 3 /min (0.0005 m 3 /s, 0.0020 m 3 /s, 0.0037 m 3 /s, 0.0050 m 3 /s, 0.0071 m 3 /s, 0.0087 m 3 /s, 0.0105 m 3 /s). In this study, a total of 168 water discharge tests were conducted, and 1176 sets of test data were obtained.
(4) In order to avoid the phenomenon of a high level of water jump in the upstream of the vegetation section and the downstream wake depth of the vegetation section, the influence of the flow instability on the test results was eliminated as much as possible. A water leveling device was set up upstream of the vegetation section, and observations of the length of the water depth and high area were performed downstream. Because the length of the upstream flat water section and the downstream water measuring section are both 1 m, it was found through observation that the upstream water jump and the downstream water flooding phenomenon did not affect the vegetation section. The length of the vegetation coverage section was 3 m, and 5 fixed longitudinal observation sections were set at equal intervals (As shown in Figure 1). The distance between adjacent sections is 0.75 m. In order to further ensure the full development of upstream water flows, this experiment used Sections 2 and 4 as the observation sections for water depth measurement. The vegetation cover observation length was 1.5 m. As a result, the influence of the instability of the upstream water flow on the accuracy of the observation data of the vegetation section was further reduced.

Data Analysis
The ratio of the total area of stem vegetation to the total area of the flume bed of the vegetation test section is used to quantify the vegetative stem cover. Its expression is as follows:

Data Analysis
The ratio of the total area of stem vegetation to the total area of the flume bed of th vegetation test section is used to quantify the vegetative stem cover. Its expression is a follows:

Data Analysis
The ratio of the total area of stem vegetation to the total area of the flume bed of the vegetation test section is used to quantify the vegetative stem cover. Its expression is as follows: where Cr is effective vegetative stem coverage (%), d is the diameter of the stem vegetation (m). N is the number of stem vegetation with diameter d. W is the width of the test flume (0.4 m), and L is the vegetation paving test section (3 m).
Because of the existence of stem vegetation, the cross-sectional area of overland flow is reduced; thus, the flume width cannot represent the effective flow width. The effective flow width is calculated as follows [61]: where w is the effective flow width of the section (m). The hydraulic radius R is the ratio of the interface water area to the wet circumference (m). It is expressed as follows: P is the wet perimeter (m). A is the effective cross-sectional area of the flow (m 2 ); it is expressed as follows: where h is the mean flow depth (m). The mean inter-stem velocity was estimated as follows: where V is the mean velocity (m s −1 ). Q is the discharge (m 3 s −1 ), which is determined by the electromagnetic flowmeter (Model: LDG-MIK-DN100).
During the experiment, we measured the pressure with piezometer tubes in Sections 2 and 4 and recorded the flow depths as h 2 and h 4 , respectively. R 2 , R 4 and V 2 , V 4 are respectively calculated by Equations (3) and (6). The flow depth h, current velocity V, and the hydraulic radius R were calculated using the mean values of cross Sections 2 and 4 (i.e., h = (h 2 + h 4 )/2; R = (R 2 + R 4 )/2; V= (V 2 + V 4 )/2).
The flow pattern is the basic parameter of the hydrodynamic characteristics. Generally, the Re and Fr are used to distinguish the flow pattern parameters. The Reynolds number is a parameter that reflects the motion state of a liquid and measures the ratio of inertial force to viscous force in the process of liquid flow. The Reynolds number can be expressed as where R is the hydraulic radius, and ν is the motion viscosity coefficient (m 2 /s). Using the Poiseuille formula, ν was calculated as: where T is the temperature of the water ( • C).
In the mechanical sense, the Froude number represents the ratio between the action of the inertial force and the action of gravity in the flow process of water. The expression of the Froude number is where g is the acceleration due to gravity (m/s 2 ). Manning roughness coefficient, n, reflects the retardation effect of bed roughness on the flow process. The flow resistance of a vegetation slope is generally believed to mainly consist of boundary resistance or particle resistance and additional resistance caused by water flowing through vegetation, which can be obtained by Equation (10).
where J is the hydraulic gradient. The hydraulic gradient is the head loss per unit distance along the water flow path, which can be obtained by Equation (11).
where l is the length of water along the course (m). In this experiment, l is the distance between Sections 2 and 4, with a length of 1.5 m. h f is the frictional head loss (m) where z 2 and z 4 are the head height of the initial (Section 2) and final positions (Section 4) is the flow head difference between the initial and final positions. The following statistical parameters and equations were used to evaluate the performance of the Manning roughness coefficient. The Nash-Sutcliffe model efficiency (NSE) was calculated as where O i , P i , and O av are the measured, analogue, and average values of the measured values, respectively, and N is the number of samples. Values of NSE closer to 1 indicate a better simulation effect of the model. Through the calculation of the above formulas, the values of relevant hydraulic parameters under different stem cover conditions in four gentle slope states were obtained, which were used to systematically study the relationship between the slope flow resistance and water depth, flow pattern, stem density, and slope. The specific experimental data are shown in Tables 1-4. Table 1. Experimental data of hydraulic parameters of six different vegetation stem cover under the condition of S = 0.0% (Cr is 0.72%, 0.37%, 0.20%, 0.13%, 0.11% and 0.07%, respectively).

Statistical Analysis
All statistical analyses were performed in First Optimization software (version 5.0). To analyze the correlations between the resistance coefficient and hydraulic parameters, gravel coverage degrees, slope gradients, and flow discharges, a regression analysis was implemented to quantify the relationship between independent and dependent variables.

The Relationship between n and h/hv under Different Density and Slope Conditions
It is generally believed that the resistance caused by vegetation to water flow acts on the water body along the water depth [62,63]; then, the flow resistance coefficient is positively related to the water depth. It should be pointed out that this experiment is related to the non-submerged state of vegetation. In order to make the experiment and the results show better physical meaning, this study did not directly choose the water depth but the ratio h/hv of the water depth to the height of the vegetation, and the study was carried out in a dimensionless form. Variation characteristics of water flow resistance coefficient under different coverage density and slope conditions. In order to further analyze the variation law of the flow resistance coefficient for overland flow with non-submerged vegetation, the test data were organized as shown in Figure 3.
is the dominant factor affecting the resistance coefficient, followed by Cr (0.1986) and S (0.0873), the latter having the smallest effect.   Figure 3 shows the trend of variation in the vegetation Manning roughness coefficient n with the values of h/hv under six coverage conditions. Under six different stem cover densities, n increased with h/hv as a power function (Figure 3). This is consistent with the conclusion that the relationship between flow resistance and water depth in Zhang et al. [64] and Juha [65] was based on the influence of the distribution of vegetation on slope flow resistance. In terms of the coefficient of resistance and values of h/hv, as the vegetation coverage increased, the correlation between the coefficient of resistance and values of h/hv strengthened, which is consistent with the results of previous literature [66,67]. When Cr ≥ 0.13%, its regression correlation coefficient R 2 was always greater than 0.8835 (mean was approximately 0.9266); when Cr ≤ 0.11%, the correlation between the resistance coefficient and h/hv was weak, and its regression correlation coefficient R 2 was less than 0.5769 (mean was approximately 0.5109). In addition, the n value of all test points did not increase from zero, which was caused by the wall friction from the plexiglass plate laid on the test floor; Jin et al. [25] and Jeon et al. [68] also had similar ideas related to the flow resistance of partially submerged stems.
The dn/d (h/hv) ratio increased with increasing vegetation coverage (Figure 3). This indicates that under high coverage, when Re is held constant and the increase in velocity is small, the water depth is bound to increase which causes the drag caused by vegetation in increase along the water depth. The drag caused by the wall friction was far less than the additional resistance of the vegetation. This is consistent with the conclusion of James et al. [69] on the additional resistance of vegetation. In addition, with a constant values of h/hv, the coverage and the resistance coefficient were positively correlated. Under the same coverage, the slope has almost no effect on the resistance coefficient.
The factors affecting resistance to overland flow in areas covered by vegetation mainly include ratio of water depth to vegetation height, coverage, slope, flow, and so on. The power function formula between h/hv, Cr, S and n can be obtained using data regression analysis. The NSE values in Equations (14) and (15) were 0.8571 and 0.8528, respectively, which are close to 1, indicating that the equation is very good. The absolute value of the index of each variable can reflect the sensitivity of the change to the resistance coefficient. It can be seen that under the condition of zero slope, Cr (0.2968) has a slightly greater influence on the drag coefficient than h/hv (0.2941), so it is the dominant factor affecting the overland flow resistance. Under the condition of non-zero slope, h/hv (0.3472) is the dominant factor affecting the resistance coefficient, followed by Cr (0.1986) and S (0.0873), the latter having the smallest effect.

The Relationship between n and Re under Different Density and Slope Conditions
For the open channel flow of submerged vegetation, n is a constant, which does not change with a change in water depth, Reynolds number, and relative roughness, but is related to the boundary material. Due to the influence of surface conditions such as vegetation shape, stiffness, density, and distribution form, the relationship between n and Re of overland flow with vegetation is more complicated than that of open channel flow, and its change trend is closely related to the composition of rough vegetation units. Figure 4 plots the variation law of n and Re under different coverage, which can be used to explore the characteristics of vegetation resistance.

The Relationship between n and Fr under Different Density and Slope Conditions
The resistance coefficient is closely related to the Reynolds number, but also to its Froude number. The Fr is used to indicate whether the overland flow is in a rapid-flow or slow-flow state and represents the ratio of inertia force to gravity. When the flow rate is  Figure 4 shows the range of the Re value was 1218-24,155; the research flow patterns were transition flow and turbulent flow. Under six coverage conditions, n increased with the increase of Re, which is consistent [11,70] or inconsistent [8,[71][72][73] with the research conclusions of various scholars. This is true because the former group of scholars studied turbulent flow, while the latter group studied laminar, from which it can be inferred that the flow pattern of the overland flow is the main factor affecting the relationship between the resistance coefficient and Re. This phenomenon may occur because the composition of the flow resistance is different in different conditions. In the case where Re is very small, because the flow is not easily spread across the entire bottom slope, the boundary resistance generated by the inner wall of the flume has an absolute advantage over the morphological resistance due to the presence of vegetation; as Re increases, the boundary resistance is continuously decreasing while the morphological resistance is prominent.
When Cr is relatively large, the rate of change of n with Re is relatively large; that is, the larger Re is, the greater the increment of n is. This is true because an increase in vegetation density reduces the cross section of water area, increases the wet cycle of vegetation flow, and increases the solid-liquid contact area, so the vegetation change rate is greater. It can also be seen from Figure 4 that when Cr ≥ 0.13%, under the condition of the same coverage density and the same Reynolds number, a larger S corresponds to a smaller n value. When Cr decreases further ≤ 0.11%, this phenomenon gradually disappears.
The occurs because as the slope continues to increase, the potential energy of the flow and the inertial force of the flow both increased continuously; the turbulence of the flow also increased continuously, causing under the same Re conditions. The thickness of the water layer on the slope is lessened, the solid-liquid contact area is continuously reduced, the Manning roughness coefficient decreases, and the flow resistance capacity of slope decreases correspondingly. However, when the amount of vegetation coverage is relatively small, the influence of slope on the resistance coefficient is difficult to reflect because the morphological resistance of vegetation plays only a small role.
The power function formula between Re, Cr, S and n can be obtained based on the regression analysis data. Using Equations (16) and (17)

The Relationship between n and Fr under Different Density and Slope Conditions
The resistance coefficient is closely related to the Reynolds number, but also to its Froude number. The Fr is used to indicate whether the overland flow is in a rapid-flow or slow-flow state and represents the ratio of inertia force to gravity. When the flow rate is basically held unchanged, with a larger Froude number the average flow velocity of overland runoff will be larger, runoff will have a stronger ability to cause erosion, the water depth will be shallower, and the runoff on the slope will be smaller. Previous studies have considered that there is a power function relationship between the resistance coefficient and the Froude number, but the influence of coverage and slope has seldom been considered. In order to fully consider the relationship between the resistance coefficient and the Froude number under different coverage and slope conditions, the test data were organized as shown in Figure 5.
within the slope range of this test, Fr was the leading factor affecting the resistance coefficient, followed by S and finally Cr.

The Relationship between n and Cr and S under Different Flow Discharge Conditions
The density of vegetation coverage and the slope are important factors affecting the change of resistance of sloping surface. In order to further study the influence of coverage density and slope on the Manning roughness coefficient, the experimental data are organized as shown in Figures 6 and 7. Under the four slope conditions, n gradually increased with an increase of Cr, and dn/dCr gradually decreased with an increase of slope, especially when the S is 2.0% or 3.0%; with Cr continuously increasing from 0.07% to 0.37%, n was basically in a steady change ( Figure 6). In addition, under the same coverage conditions, n increased with an increase of Q and the two were positively correlated. This finding was consistent with some previous results showing that the Manning roughness coefficient increased with an increasing degree of coverage [19,[33][34][35]. This occurred because under the same flow conditions with greater vegetation coverage the double-layer advantages of the number of plants and the thickness of the backwater increased the solidliquid contact area, and the water blocking ability was bound to increase.  Figure 5). With an increase of Fr, n showed a decreasing trend, which is similar to the conclusions of Li et al. [70] and Zhang et al. [11]; Fr basically presents the change rule of a power function. When Fr is constant, a larger S resulted in a greater resistance coefficient. As coverage decreases, the correlation between n and Fr weaker, and the range of n values is shranking ( Figure 5). Under the same coverage conditions, as the slope increased, the flow from slow to rapid development improved. In addition, when the slope S ≤ 1.0% the flow state under this study conditions was slow flow; meanwhile, when the slope S > 1.0% the flow changed to rapids flow.
Using the data regression analysis, the power function formula for Fr, Cr, S and n could be obtained. Compared with other hydraulic parameters, the correlation was slightly worse, and the NSE values in Equations (18) and (19) were 0.5587 and 0.7613, respectively. According to the fitting formula, under a constant Fr condition, Cr and S were positively correlated with n, while Fr and n were negatively correlated. In addition, within the slope range of this test, Fr was the leading factor affecting the resistance coefficient, followed by S and finally Cr.

The Relationship between n and Cr and S under Different Flow Discharge Conditions
The density of vegetation coverage and the slope are important factors affecting the change of resistance of sloping surface. In order to further study the influence of coverage density and slope on the Manning roughness coefficient, the experimental data are organized as shown in Figures 6 and 7. Under the four slope conditions, n gradually increased with an increase of Cr, and dn/dCr gradually decreased with an increase of slope, especially when the S is 2.0% or 3.0%; with Cr continuously increasing from 0.07% to 0.37%, n was basically in a steady change ( Figure 6). In addition, under the same coverage conditions, n increased with an increase of Q and the two were positively correlated. This finding was consistent with some previous results showing that the Manning roughness coefficient increased with an increasing degree of coverage [19,[33][34][35]. This occurred because under the same flow conditions with greater vegetation coverage the double-layer advantages of the number of plants and the thickness of the backwater increased the solid-liquid contact area, and the water blocking ability was bound to increase.
However, some inconsistent opinions have been expressed related to the relationship between the resistance coefficient and slope gradient. Wang et al. [74] study showed that under the steep slope test conditions where the flow regime is laminar flow, when the vegetation coverage is greater than 4.65%, the resistance coefficient decreases with the increase of the slope, that is, the two are in a negative correlation. When the vegetation coverage is less than 4.65%, the slope has little effect on resistance coefficient, and there is basically no change with the increase in slope. Pan and Shangguan [73,75] showed that under the condition of laminar flow steep slope, for bare slope, the surface water resistance coefficient increased with the increase of slope, while the slope vegetation coverage was 70%, showing an opposite trend. It can be seen that the vegetation coverage is an important factor affecting the trend change between the resistance coefficient and the slope. However, under the condition of gentle slope with transitional flow and turbulent flow, the influence of sparse stem vegetation coverage on the relationship between overland flow resistance coefficient and slope remains to be studied. In order to further explore the inherent relationship between the two, the test data is sorted out as shown in Figure 7.
When the stem coverage density was 0.72% to 0.13%, n decreased with an increase of S; that is, a negative correlation existed between n and S. As the coverage decreased further (Cr ≤ 0.11%), the slope had little effect on the roughness coefficient. This indicates that with the decrease of vegetation cover density to a certain extent, particle resistance plays a dominant role in water flow resistance, while vegetation form resistance plays a small role, which is similar to the conclusion of Zhang et al. [76] that "particle resistance is not or rarely affected by slope S", so Manning roughness coefficient tends to be stable with the increase of slope. However, some inconsistent opinions have been expressed related to the relationship between the resistance coefficient and slope gradient. Wang et al. [74] study showed that under the steep slope test conditions where the flow regime is laminar flow, when the  When the stem coverage density was 0.72% to 0.13%, n decreased with an increase of S; that is, a negative correlation existed between n and S. As the coverage decreased further (Cr ≤ 0.11%), the slope had little effect on the roughness coefficient. This indicates that with the decrease of vegetation cover density to a certain extent, particle resistance plays a dominant role in water flow resistance, while vegetation form resistance plays a small role, which is similar to the conclusion of Zhang et al. [76] that "particle resistance is not or rarely affected by slope S", so Manning roughness coefficient tends to be stable with the increase of slope.
Because Figures 6 and 7 both reflect the relationship between Cr, S, and n under different flow discharges, the power function formula between Cr, S, Q, and n could be obtained based on the data regression analysis. In Equations (20) and (21), the NSE values were 0.8550 and 0.8095, respectively, and was close to 1 showing equation was optimal. According to the fitting formula, Cr, Q and n are positively correlated under the same flow condition, negative correlation existed between S and n. It can be seen that under the condition of zero slope, the influence of Cr (0.2693) on the flow resistance of slope is obviously greater than that of Q (0.1994), so it is the dominant factor affecting the overland flow resistance. Under the condition of non-zero slope, Q (0.2997) is the dominant factor affecting the resistance coefficient, followed by Cr (0.2540) and S (0.1529), the latter having smallest effect. Because Figures 6 and 7 both reflect the relationship between Cr, S, and n under different flow discharges, the power function formula between Cr, S, Q, and n could be obtained based on the data regression analysis. In Equations (20) and (21), the NSE values were 0.8550 and 0.8095, respectively, and was close to 1 showing equation was optimal. According to the fitting formula, Cr, Q and n are positively correlated under the same flow condition, negative correlation existed between S and n. It can be seen that under the condition of zero slope, the influence of Cr (0.2693) on the flow resistance of slope is obviously greater than that of Q (0.1994), so it is the dominant factor affecting the overland flow resistance. Under the condition of non-zero slope, Q (0.2997) is the dominant factor affecting the resistance coefficient, followed by Cr (0.2540) and S (0.1529), the latter having smallest effect.

Conclusions
For overland flows with discharges ranging from 0.0313 to 0.6319 m 3 /min and a slope from 0% to 3.0%, a number of hydraulic parameters were measured with a vegetative stem cover from 0.07% to 0.72% to investigate the effects of the resistance to overland flow under the condition of sparse stem cover. The following conclusions can be drawn based on this experimental study: The value of n showed a power function increasing with an increase of h/hv; the correlation between n and h/hv was stronger with an increase in vegetation cover. Meanwhile, dn/d (h/hv) increased with an increase of Cr. Under the condition of zero slope, Cr was the dominant factor affecting the overland flow resistance. Under the condition of non-zero slope, h/hv was the dominant factor affecting the resistance coefficient, followed by Cr and S, the latter having smallest effect.
In addition, n increased with an increase of Re; when Cr was greater, the rate of change of n with Re was larger. Also, Cr and Re were positively correlated with n. Under the same Fr both Cr and S were positively correlated with n, and Fr and n were negatively correlated with each other.
The value of n gradually increased with an increase of Cr; the dn/dCr ratio gradually decreased with an increase of S. Under the same Cr, n had a positive correlation with Q. Meanwhile, n decreased with increasing S, so the two had a negative correlation. As the coverage decreased further (Cr ≤ 0.11%), S had little effect on the roughness coefficient.
In short, the relationship between n and h/hv, Re, Cr, and S all shows the same law; that is, for sparse stem vegetation coverage, Cr is the dominant factor affecting overland flow resistance under zero slope condition, while Cr is no longer the first dominant factor affecting overland flow resistance under non-zero slope condition. Therefore, Cr is an important factor affecting slope flow resistance and hydraulic parameters.

Data Availability Statement:
The authors confirm that the data supporting the findings of this study are available within the supporting information.