Research and Application of the Calculation Method of River Roughness Coefﬁcient with Vegetation

: The roughness coefﬁcient is a comprehensive parameter reﬂecting river resistance, which is widely used in the planning and design of river regulation and ﬂood control projects. In recent years, as the upstream water conservancy and hydro-power projects have been put into operation, the frequency of low ﬂow in the middle and lower reaches has increased, and the frequency of ﬂood ﬂow has decreased. All kinds of vegetation in the river ﬂoodplain grow luxuriantly, which causes a change in the river resistance and roughness coefﬁcient. The present study was carried out with theoretical analysis and laboratory tests. A formula for the roughness coefﬁcient calculation was derived based on the momentum equilibrium equation and momentum exchange between the vegetation layer and upper layer. The relationship between the depth-averaged velocity within the vegetation layer and depth-averaged velocity of the whole ﬂow was analyzed. The reliability of the formula was veriﬁed by a large amount of previous experimental data. Based on the derived formula, the variation law of the roughness coefﬁcient with vegetation density, vegetation height, and water depth were obtained. For the emerged vegetation ﬂow, the Manning coefﬁcient tended to increase with the increase in the vegetation density and water depth. For the submerged vegetation ﬂow, the Manning coefﬁcient showed a trend of decreasing with the increase in the water depth and increased with the increase in the vegetation height. Finally, the derived formula was applied in the Yueyang reach of the Yangtze River and the Duliujian River. The study can be applied in the ﬁelds of water level-ﬂow discharge relationship analysis and the water surface line calculation of vegetated rivers.


Introduction
Aquatic vegetation is widely found in rivers, wetlands, and coasts.It has the functions of protecting banks, improving the water quality, and maintaining ecosystem diversity.Costanza et al. [1] estimated that the value of ecological services, which include providing nutrients and habitat, controlling erosion, protecting banks, etc., provided by aquatic vegetation was as high as USD 10 trillion per year.Aquatic vegetation is also widely used in river regulation and management.Due to the extensive engineering applications and important ecological functions of aquatic vegetation, studies related to it have received widespread attention in recent years [2][3][4] and have gradually developed into the field of vegetation hydraulics.
Vegetation resistance has become a key field of vegetation hydraulics research because of its importance in modifying flow velocity.Typically, studies of vegetation resistance focus on the Manning coefficient n, friction coefficient f, or shear stress of vegetation flow [5][6][7].In the early stage, researchers established the relationship between the Manning coefficient n and averaged velocity U and hydraulic radius R by field observations or laboratory experiments [8].Kouwen et al. [9] introduced the stiffness parameter MEI to Water 2023, 15, 2638 2 of 12 the relationship, and MEI reflected the effect of vegetation deformation on the roughness coefficient.The U.S. Soil Conservation Service (SCS) has developed tables of roughness coefficients corresponding to various types of vegetated channels to facilitate engineering applications.Some researchers have used the friction coefficient f to describe vegetation resistance.Kouwen and Unny [10], based on flume tests, inferred that the friction coefficient f caused by flexible plastic sheets was related to the relative roughness and oscillation mechanism of plastic sheets.Jarvela [11] obtained the Darcy-Weisbach coefficient f as a function of Reynolds number, flow velocity, water depth, and relative roughness based on flume tests with grass and willow.As measurement techniques advanced, researchers gradually investigated vegetation resistance through systematic model tests.For nearly rigid vegetation such as willow and reed, the vegetation was simplified as a cylindrical or strip shape.Based on the quadratic relationship between the resistance and flow velocity, the variation in the resistance coefficient with flow velocity and vegetation density was analyzed.Stone and Shen [12] used the maximum flow velocity between vegetation as the characteristic velocity.The drag coefficient did not change with the vegetation volume fraction, vegetation diameter, and vegetation Reynolds number when using the averaged velocity within vegetation as the characteristic velocity drag coefficient increased with the increase in the vegetation volume fraction and decreased with the increase in vegetation Reynolds number.Tanino and Nepf [13] conducted a detailed experimental study on the drag coefficient of randomly arranged rigid vegetation and obtained the relationship between the coefficient C D and vegetation Reynolds number and vegetation volume fraction.Yang and Choi [14] calculated the depth-averaged velocity based on solving the velocity distribution in the vegetation layer and the layer above the vegetation, and then obtained the Manning coefficient with vegetation according to the Manning formula.They proposed a coefficient C u in the velocity distribution of the layer above the vegetation, and C u = 1 for a ≤5.0 m −1 , C u = 2 for a >5.0 m −1 .C u was found to vary significantly among different vegetation, which caused a large deviation from the predicted Manning coefficient.Li S et al. [15] divided the vegetated flow into the suspension layer and basal layer, based on which the averaged velocity over the whole flow depth could be derived, and then the Manning coefficient was derived from the Manning formula.They introduced the representative length scale h* into the calculation, which was calculated with energy slope.However, the energy slope is very difficult to obtain in practice.Cheng [16] proposed a representative roughness height to quantify the effect of submerged vegetation on flow resistance in the surface layer, and then developed an approach to estimate the average flow velocity and resistance coefficients for both cases of rigid and flexible vegetation, which was also the calculation of the energy slope.Overall, the methods calculating the Manning coefficient proposed by previous researchers are relatively complicated and need many calculation parameters, especially the energy slope, which is difficult to obtain.Because the research object and derivation process of each method is different, the validation data used are also different, so the resistance coefficients predicted by each method also differ.Wang et al. [17] found that the resistance coefficients of the emerged vegetation calculated by different methods were basically the same, but the calculated resistance coefficients of inundated vegetation varied greatly, and the differences among the results of different methods could be up to 2~4 times.
In flood level prediction and river evolution simulation for rivers containing vegetation, the influence of vegetation is often considered by increasing the resistance.A careful examination and choice of the roughness coefficient is essential.Large deviations in the prediction of vegetation resistance would bring large errors to the aforementioned computational simulations, which would in turn affect the practice of river regulation.For example, the Hec-RAS model result was more accurate when the resistance coefficient was calculated dynamically using the characteristic vegetation [17].Therefore, how to calculate the roughness coefficient of the river containing vegetation easily and scientifically is an importance issue that needs to be solved.In this paper, based on theoretical analysis and Water 2023, 15, 2638 3 of 12 data validation, a relatively convenient and accurate method for calculating the Manning coefficient with vegetation is proposed.

Calculation Method of Roughness Coefficient with Vegetation
The vegetation-induced drag F D can be described by using the drag coefficient C D .
where F D is the vegetation resistance per unit mass of fluid; C D is the vegetation drag coefficient; a is the vegetation density defined by the projected area in the incoming flow plane per unit volume; a = N × d; N is the number of vegetation per unit area; d is the vegetation characteristic length, for cylindrical vegetation, the characteristic length can be used with diameter; U 1 is the mean velocity in the vegetated layer.According to previous research [6,7], the vegetation resistance coefficient varies with the vegetation Reynolds number, vegetation density, etc.For low vegetation density, the drag coefficient C D is generally taken as 1.0.Assuming the vegetation flow is constant and uniform, and the vegetation resistance, bed resistance, and gravity are balanced along the direction of flow, the equation for force balance with submerged vegetation (see Figure 1) is where U is the depth-averaged velocity over the whole flow depth and U 1 is the depthaveraged velocity within the vegetation layer.When vegetation is submerged, U 1 < U, when vegetation is emerged, U 1 = U; a is the vegetation density; h is the vegetation height below the water surface; H is the water depth when the vegetation is emerged, h = H; g is the acceleration of gravity; f is the Darcy-Weisbach coefficient; S 0 is the water surface slope; C D is the vegetation drag coefficient; Φ is the vegetation volume fraction; for cylindrical vegetation, Φ equals adπ/4 and d is the vegetation diameter.
an importance issue that needs to be solved.In this paper, based on theoretical analysis and data validation, a relatively convenient and accurate method for calculating the Man ning coefficient with vegetation is proposed.

Calculation Method of Roughness Coefficient with Vegetation
The vegetation-induced drag FD can be described by using the drag coefficient CD.
where FD is the vegetation resistance per unit mass of fluid; CD is the vegetation drag co efficient; a is the vegetation density defined by the projected area in the incoming flow plane per unit volume; a = N × d; N is the number of vegetation per unit area; d is the vegetation characteristic length, for cylindrical vegetation, the characteristic length can be used with diameter; U1 is the mean velocity in the vegetated layer.According to previous research [6,7], the vegetation resistance coefficient varies with the vegetation Reynolds number, vegetation density, etc.For low vegetation density, the drag coefficient CD is gen erally taken as 1.0.
Assuming the vegetation flow is constant and uniform, and the vegetation resistance bed resistance, and gravity are balanced along the direction of flow, the equation for force balance with submerged vegetation (see Figure 1) is where U is the depth-averaged velocity over the whole flow depth and U1 is the depth averaged velocity within the vegetation layer.When vegetation is submerged, U1 < U when vegetation is emerged, U1 = U; a is the vegetation density; h is the vegetation heigh below the water surface; H is the water depth when the vegetation is emerged, h = H; g is the acceleration of gravity; f is the Darcy-Weisbach coefficient; S0 is the water surface slope; CD is the vegetation drag coefficient; Φ is the vegetation volume fraction; for cylin drical vegetation, Φ equals adπ/4 and d is the vegetation diameter.The water surface slope S0 and Darcy-Weisbach coefficient f can be calculated from the following equation, where n is the Manning coefficient of the river with vegetation and R is the hydraulic radius.Generally, the river width is much larger than the water depth, so the hydraulic radius can be expressed by the mean water depth H; C is the Chezy coefficient; n0 is the roughness coefficient without vegetation.The water surface slope S 0 and Darcy-Weisbach coefficient f can be calculated from the following equation, where n is the Manning coefficient of the river with vegetation and R is the hydraulic radius.Generally, the river width is much larger than the water depth, so the hydraulic radius can be expressed by the mean water depth H; C is the Chezy coefficient; n 0 is the roughness coefficient without vegetation.
Combining Equation (2) with Equation (3) to Equation ( 5), the Manning coefficient n with vegetation was obtained as When the vegetation is emerged, U 1 = U, and the Manning coefficient, Equation ( 6), with vegetation can be simplified as When the vegetation is submerged, the vegetation flow can be divided into twolayer models.The relationship between U 1 and U can be expressed as follows (see Chen et al. [18]).
The coefficient C describes the efficiency of momentum exchange between the vegetation layer and non-vegetation layer.C can be considered as a constant, taking the value of 0.009.Substituting Equation ( 8) into Equation ( 7), the Manning coefficient with submerged vegetation flow can be expressed as

Materials
Previous researchers have conducted many experiments (see Table 1) including Shimizu et al. [19], Dunn et al. [20], Meijer and van Velzen [21], Lopez and Garcia [22], Ghisalberti and Nepf [23], Murphy et al. [24], Nezu and Sanjou [25], Yan [26], Yang [27], and Cheng [16].With these experiments, the calculation method can be validated.The main experimental parameters of the above-mentioned researchers were as follows.The test flume width ranged from 0.3 m to 3.0 m, water depth ranged from 0.05 m to 2.5 m, energy slope ranged from 0.018‰ to 16‰, flow velocity ranged from 0.01 m/s to 1.24 m/s, vegetation height ranged from 0.035 m to 1.5 m, vegetation diameter ranged from 0.001 m to 0.008 m, vegetation density ranged from 0.27 m −1 to 29 m −1 , and the vegetation volume fraction ranged from 0.14% to 11.9%.n 0 was determined by the corresponding test results, and where the model test was not given, it was generally taken as 0.011.

Model Validation
The comparison with the calculated Manning coefficient using Equation ( 9) with the measured Manning coefficient in each test is shown in Figure 2. The validation results showed that the Manning roughness calculated by Equation ( 9) did not differ much from the measured coefficient, with a standard deviation of 17.5%.The comparisons with the calculated Manning coefficient using Yang and Choi's equation [8] and Li et al.'s [9] equation with the measured Manning coefficient in each test are shown in Figures 3 and 4. The R-square of Equation ( 9) was 0.88, the R-square of Yang and Choi's equation was −10.83, and R-square of Li et al.'s equation was 0.85.It can be seen that the model accuracy of Equation ( 9) was much better than Yang and Choi's equation and comparable to that of Li et al.Considering that Equation ( 9) involves the parameters of the natural property of the river and vegetation, it does not utilize the water surface slope, which is commonly used by other researchers [14,15] but is difficult to obtain.We believe that the practical value of Equation ( 9) would be higher than the common formula for calculating the roughness coefficient with vegetation.

Model Validation
The comparison with the calculated Manning coefficient using Equation ( 9) with the measured Manning coefficient in each test is shown in Figure 2. The validation results showed that the Manning roughness calculated by Equation ( 9) did not differ much from the measured coefficient, with a standard deviation of 17.5%.The comparisons with the calculated Manning coefficient using Yang and Choi's equation [8] and Li et al.'s [9] equation with the measured Manning coefficient in each test are shown in Figures 3 and 4 9) involves the parameters of the natural property of the river and vegetation, it does not utilize the water surface slope, which is commonly used by other researchers [14,15] but is difficult to obtain.We believe that the practical value of Equation ( 9) would be higher than the common formula for calculating the roughness coefficient with vegetation.

Variation of Roughness Coefficient
In this section, Equations ( 6) and ( 9), derived in the previous section, were used to analyze the effects of variation in the vegetation density, vegetation height, and water depth on the Manning coefficient of natural river floodplains.Natural river floodplains consist of trees and shrubs, reeds, etc., which present an emerged or submerged state when the river is in a flood stage.
(1) Emerged vegetation flow Trees, which are used to protect embankments, are commonly planted in the floodplain.According to the field survey of the Jingjiang River of the Yangtze River, the averaged tree diameter is about 0.1 m, and the density a is generally between 0.01 and 0.2 m −1 .Based on the experience and the hydraulic calculation manual, the Manning coefficient n of the floodplain without vegetation was taken as 0.03.According to Equation (7), the variation of the Manning coefficient n with the vegetation density and floodplain water depth is shown in Figure 5.It can be seen that for emerged vegetation, the Manning coefficient tends to increase with the increase in the vegetation density and floodplain water depth.The

Variation of Roughness Coefficient
In this section, Equations ( 6) and ( 9), derived in the previous section, were used to analyze the effects of variation in the vegetation density, vegetation height, and water depth on the Manning coefficient of natural river floodplains.Natural river floodplains consist of trees and shrubs, reeds, etc., which present an emerged or submerged state when the river is in a flood stage.
(1) Emerged vegetation flow Trees, which are used to protect embankments, are commonly planted in the floodplain.According to the field survey of the Jingjiang River of the Yangtze River, the averaged tree diameter is about 0.1 m, and the density a is generally between 0.01 and 0.2 m −1 .Based on the experience and the hydraulic calculation manual, the Manning coefficient n of the floodplain without vegetation was taken as 0.03.According to Equation (7), the variation of the Manning coefficient n with the vegetation density and floodplain

Variation of Roughness Coefficient
In this section, Equations ( 6) and ( 9), derived in the previous section, were used to analyze the effects of variation in the vegetation density, vegetation height, and water depth on the Manning coefficient of natural river floodplains.Natural river floodplains consist of trees and shrubs, reeds, etc., which present an emerged or submerged state when the river is in a flood stage.
(1) Emerged vegetation flow Trees, which are used to protect embankments, are commonly planted in the floodplain.According to the field survey of the Jingjiang River of the Yangtze River, the averaged tree diameter is about 0.1 m, and the density a is generally between 0.01 and 0.2 m −1 .Based on the experience and the hydraulic calculation manual, the Manning coefficient n of the floodplain without vegetation was taken as 0.03.According to Equation (7), the variation of the Manning coefficient n with the vegetation density and floodplain (2) Submerged vegetation flow Flexible vegetation such as weeds, crops, and low shrubs are also found in the floodplain.This kind of vegetation bends with the current under the action of water flow and is often in the submerged state.Combined with the relevant literature and field survey, the density of weeds and other flexible vegetation in the floodplain is generally 0.2~5 m −1 .When vegetation is in the submerged state, the height h of vegetation below the water surface is the effective height after bending under the action of water flow.According to Equation ( 9), the relationship of the Manning coefficient with vegetation density, floodplain water depth, and vegetation height is shown in Figures 6 and 7.Under the same vegetation height condition, the greater the water depth, the smaller the Manning coefficient, and the greater the vegetation density, the larger the Manning coefficient; the Manning coefficient increases with the increase in the vegetation density and decreases with the increase in the floodplain water depth.Under the same condition of vegetation height h = 0.5 m, the Manning coefficient was 0.050 when the water depth was 1 m and the vegetation density was 0.2 m −1 , 0.097 when the water depth was 1 m and the vegetation density was 5 m −1 , 0.043 when the water depth was 3 m and the vegetation density was 0.2 m −1 , and 0.063 when the water depth was 1 m and the vegetation density was 5 m −1 .

2) Submerged vegetation flow
Flexible vegetation such as weeds, crops, and low shrubs are also found in the floodplain.This kind of vegetation bends with the current under the action of water flow and is often in the submerged state.Combined with the relevant literature and field survey, the density of weeds and other flexible vegetation in the floodplain is generally 0.2~5 m −1 .When vegetation is in the submerged state, the height h of vegetation below the water surface is the effective height after bending under the action of water flow.According to Equation ( 9), the relationship of the Manning coefficient with vegetation density, floodplain water depth, and vegetation height is shown in Figures 6 and 7.Under the same vegetation height condition, the greater the water depth, the smaller the Manning coefficient, and the greater the vegetation density, the larger the Manning coefficient; the Manning coefficient increases with the increase in the vegetation density and decreases with the increase in the floodplain water depth.Under the same condition of vegetation height h = 0.5 m, the Manning coefficient was 0.050 when the water depth was 1 m and the vegetation density was 0.2 m −1 , 0.097 when the water depth was 1 m and the vegetation density was 5 m −1 , 0.043 when the water depth was 3 m and the vegetation density was 0.2 m −1 , and 0.063 when the water depth was 1 m and the vegetation density was 5 m −1 .
Under the same vegetation density condition, the Manning coefficient also showed a phenomenon of decreasing with increasing water depth and increasing with vegetation height.Under the condition that the vegetation density is constant (a = 1.0 m −1 ), the Manning coefficient increased from 0.03 to 0.257 during the increase in vegetation height from 0 to 1.0 m at a floodplain water depth of 1.0 m, and from 0.03 to 0.064 during the increase in vegetation height from 0 to 1.0 at a vegetation water depth of 3.0 m.Under the same vegetation density condition, the Manning coefficient also showed a phenomenon of decreasing with increasing water depth and increasing with vegetation height.Under the condition that the vegetation density is constant (a = 1.0 m −1 ), the Manning coefficient increased from 0.03 to 0.257 during the increase in vegetation height from 0 to 1.0 m at a floodplain water depth of 1.0 m, and from 0.03 to 0.064 during the increase in vegetation height from 0 to 1.0 at a vegetation water depth of 3.0 m.
The above analysis shows that the Manning coefficient increases significantly when vegetation is present on the floodplain.The Manning coefficient of vegetation generally varies between 0.03 and 0.2.The Manning coefficient of low and sparse vegetation is smaller, while the Manning coefficient of tall and dense vegetation is larger.The above variation pattern can be explained by the vegetation water retention mechanism.The main reason for vegetation water blockage is that the water flow produces vortices of different scales under the action of vegetation including vortices formed by vegetation branches and leaves, and shear vortices formed by submerged vegetation tops.Vortices of different sizes dissipate part of the kinetic energy of the water flow into internal energy, which leads to the slowing down of the water flow and the formation of the water blocking effect.For non-submerged vegetation, the greater the water depth, the more vortices are produced by vegetation, so the greater the energy lost, and the greater the force of vegetation on water flow; the greater the density of the vegetation, the same principle applies, resulting in the phenomenon of increasing the water depth or vegetation density and increasing the Manning coefficient.For submerged vegetation, the water depth of the floodplain increases, therefore the flow velocity of the vegetation layer will decrease, the vortex generated by the vegetation layer decreases, and the energy lost decreases, resulting in the phenomenon of an increase in the water depth and a decrease in the stratification coefficient.The above analysis shows that the Manning coefficient increases significantly when vegetation is present on the floodplain.The Manning coefficient of vegetation generally varies between 0.03 and 0.2.The Manning coefficient of low and sparse vegetation is smaller, while the Manning coefficient of tall and dense vegetation is larger.The above variation pattern can be explained by the vegetation water retention mechanism.The main reason for vegetation water blockage is that the water flow produces vortices of different scales under the action of vegetation including vortices formed by vegetation branches and leaves, and shear vortices formed by submerged vegetation tops.Vortices of different sizes dissipate part of the kinetic energy of the water flow into internal energy, which leads to the slowing down of the water flow and the formation of the water blocking effect.For non-submerged vegetation, the greater the water depth, the more vortices are produced by vegetation, so the greater the energy lost, and the greater the force of vegetation on water flow; the greater the density of the vegetation, the same principle applies, resulting in the phenomenon of increasing the water depth or vegetation density and increasing

Practical Applications
In this section, we applied the roughness calculation method to two typical river sections in the Yangtze and Haihe Basins to illustrate the practicality of the method.

Emerged Vegetation Flow
The Duliujian River is an important flood passage of Tianjin.The research subject is a typical reach of the upper Duliujian River (see Figure 8).Trees are the typical vegetation of the river reach.The average diameter of trees is 6 cm, the number of trees per unit area is 0.18, the height of trees is 10 m, and the density of trees is 0.011 m −1 .The averaged elevation of the river floodplain is about 2.5 m, and the one hundred-year flood level is 5.73 m.During the one hundred-year flood, the averaged depth of floodplain is 3.23 m, and the trees are in an emerged state.The floodplain roughness coefficient without vegetation was taken as 0.025 (the parameters above can be found in Gao et al. [28]).According to Equation (7), the roughness coefficient of the river with trees is 0.057.This result is basically consistent with the roughness coefficient of 0.06 (see Gao et al. [28]) derived from the experiment.

Submerged Vegetation Flow
The middle reaches of the Yangtze River (see Figure 9) are a key area for flood control, and there is much vegetation on the river floodplain.The Yueyang River section is located below the confluence of the Jing River and Dongting Lake, which faces a serious flood control situation in the middle reaches of the Yangtze River.Therefore, the Yueyang River reach was chosen as the study site.This river reach is a compound channel.The width of the river is about 2500 m, the width of the left floodplain is about 400-500 m, and the width of the right floodplain is less than 100 m.Shrubs and crops predominate on the floodplain.vation of the river floodplain is about 2.5 m, and the one hundred-year flood level is 5.73 m.During the one hundred-year flood, the averaged depth of floodplain is 3.23 m, and the trees are in an emerged state.The floodplain roughness coefficient without vegetation was taken as 0.025 (the parameters above can be found in Gao et al. [28]).According to Equation ( 7), the roughness coefficient of the river with trees is 0.057.This result is basically consistent with the roughness coefficient of 0.06 (see Gao et al. [28]) derived from the experiment.

Submerged Vegetation Flow
The middle reaches of the Yangtze River (see Figure 9) are a key area for flood control, and there is much vegetation on the river floodplain.The Yueyang River section is located below the confluence of the Jing River and Dongting Lake, which faces a serious flood control situation in the middle reaches of the Yangtze River.Therefore, the Yueyang River reach was chosen as the study site.This river reach is a compound channel.The width of the river is about 2500 m, the width of the left floodplain is about 400-500 m, and the width of the right floodplain is less than 100 m.Shrubs and crops predominate on the floodplain.
The average diameter of vegetation is about 1 cm, the number of vegetation per unit area is 200 plants/m 2 , the height of vegetation is 50 cm, and the density of vegetation is 2 m −1 .The averaged elevation of the floodplain is about 28 m, and the design flood level of this reach is 32.36 m.Under the design flood level, the average water depth of the floodplain is about 4.36 m.The height of vegetation after bending is 0.4 m under the design flood conditions.According to Equation ( 9), the Manning coefficient of the floodplain is 0.057.Due to the completion and operation of the Three Gorges Project, the chances of where Pi, ni are the wetted perimeter and the Manning coefficient of the ith segmented section, respectively.The typical section of Figure 10 can be divided into two zones, namely, the main channel zone and the floodplain zone.The wetted perimeter of the floodplain and the main channel was 450 m and 2007 m, respectively.Combined with the results of the previous analysis, the Manning coefficients of the floodplain before and after the change in the vegetation density were 0.054 and 0.057, respectively; the Manning coefficient of the main channel remained unchanged at 0.025.The integrated roughness of the section before and after the change in vegetation density was 0.0303 and 0.0309, respectively.The integrated roughness of the section increased after an increase in the vegetation density by about 2.0%.The above study shows that computational parameters of this method are easy to obtain, and the computational process is simple.The method has achieved good results in the water flow of trees, shrubs, and other vegetation.We can see that the calculation method is mainly established based on rigid or near-rigid vegetation, and the momentum exchange coefficient between the vegetation and non-vegetation layers can also be determined based on rigid or near-rigid vegetation.Therefore, we believe that the model proposed in this paper should be mainly used for vegetation such as trees, reeds, and shrubs as the method application may have large errors for flexible vegetation with large deformation under the action of water flow; in cases where there are more leaves or vegetation, density a is difficult to assess accurately.The average diameter of vegetation is about 1 cm, the number of vegetation per unit area is 200 plants/m 2 , the height of vegetation is 50 cm, and the density of vegetation is 2 m −1 .The averaged elevation of the floodplain is about 28 m, and the design flood level of this reach is 32.36 m.Under the design flood level, the average water depth of the floodplain is about 4.36 m.The height of vegetation after bending is 0.4 m under the design flood conditions.According to Equation ( 9), the Manning coefficient of the floodplain is 0.057.Due to the completion and operation of the Three Gorges Project, the chances of flooding in the river reach downstream of the Three Gorges have been significantly reduced.The vegetation coverage on the floodplain has increased compared with its status before the Three Gorges Project, in particular, the density of shrubs and other flexible vegetation has increased significantly.Assuming that the vegetation density before the completion of the Three Gorges was 1 m −1 , other parameters of the vegetation remained unchanged.According to Equation ( 9), the Manning coefficient of the floodplain is 0.054.Therefore, when the vegetation density increases from 1 m −1 to 2 m −1 , the Manning coefficient of the floodplain increases from 0.054 to 0.057, which is an increase of about 5.5%.
Next, the authors analyzed the effect of changes in floodplain resistance on river resistance.A typical cross-section was chosen, as shown in Figures 10 and 11

Conclusions
In this paper, the calculation method of the Manning coefficient was derived based on the momentum equilibrium equation and momentum exchange between the vegetation layer and upper layer.The Manning coefficient is primarily related to the vegetation characteristics and water depth.The variation law of floodplain roughness with vegetation density, vegetation height, and water depth was explored.The study showed that for emerged vegetation, the Manning coefficient tends to increase with the increase in the vegetation density and floodplain water depth.For submerged vegetation, the Manning coefficient showed a trend of decreasing with the increase in the water depth and increasing with the increase in the vegetation height.The Manning coefficient for vegetated flow generally ranges from 0.03 to 0.2.The application of Equation ( 9) in the Yueyang River section of the middle and lower reaches of the Yangtze River showed that the increase in vegetation density after the Three Gorges Project led to an increase in the floodplain Manning coefficient, which increased from 0.054 to 0.057, an increase of about 5.5%, and the comprehensive Manning coefficient of a typical section increased from 0.0303 to 0.0309, an increase of about 2.0%.The results of the study can be applied in the analysis of the water level-flow relationship and water surface line calculation of a vegetated river.It is important to note that further research is needed as to whether the present method is applicable to flexibly submerged vegetation with large deformations.

Conclusions
In this paper, the calculation method of the Manning coefficient was derived based on the momentum equilibrium equation and momentum exchange between the vegetation layer and upper layer.The Manning coefficient is primarily related to the vegetation characteristics and water depth.The variation law of floodplain roughness with vegetation density, vegetation height, and water depth was explored.The study showed that for emerged vegetation, the Manning coefficient tends to increase with the increase in the vegetation density and floodplain water depth.For submerged vegetation, the Manning coefficient showed a trend of decreasing with the increase in the water depth and increasing with the increase in the vegetation height.The Manning coefficient for vegetated flow generally ranges from 0.03 to 0.2.The application of Equation ( 9) in the Yueyang River section of the middle and lower reaches of the Yangtze River showed that the increase in vegetation density after the Three Gorges Project led to an increase in the floodplain Manning coefficient, which increased from 0.054 to 0.057, an increase of about 5.5%, and the comprehensive Manning coefficient of a typical section increased from 0.0303 to 0.0309, an increase of about 2.0%.The results of the study can be applied in the analysis of the water level-flow relationship and water surface line calculation of a vegetated river.It is important to note that further research is needed as to whether the present method is applicable to flexibly submerged vegetation with large deformations.(10) where P i , n i are the wetted perimeter and the Manning coefficient of the ith segmented section, respectively.The typical section of Figure 10 can be divided into two zones, namely, the main channel zone and the floodplain zone.The wetted perimeter of the floodplain and the main channel was 450 m and 2007 m, respectively.Combined with the results of the previous analysis, the Manning coefficients of the floodplain before and after the change in the vegetation density were 0.054 and 0.057, respectively; the Manning coefficient of the main channel remained unchanged at 0.025.The integrated roughness of the section before and after the change in vegetation density was 0.0303 and 0.0309, respectively.The integrated roughness of the section increased after an increase in the vegetation density by about 2.0%.
The above study shows that computational parameters of this method are easy to obtain, and the computational process is simple.The method has achieved good results in the water flow of trees, shrubs, and other vegetation.We can see that the calculation method is mainly established based on rigid or near-rigid vegetation, and the momentum exchange coefficient between the vegetation and non-vegetation layers can also be determined based on rigid or near-rigid vegetation.Therefore, we believe that the model proposed in this paper should be mainly used for vegetation such as trees, reeds, and shrubs as the method application may have large errors for flexible vegetation with large deformation under the action of water flow; in cases where there are more leaves or vegetation, density a is difficult to assess accurately.

Conclusions
In this paper, the calculation method of the Manning coefficient was derived based on the momentum equilibrium equation and momentum exchange between the vegetation layer and upper layer.The Manning coefficient is primarily related to the vegetation characteristics and water depth.The variation law of floodplain roughness with vegetation density, vegetation height, and water depth was explored.The study showed that for emerged vegetation, the Manning coefficient tends to increase with the increase in the vegetation density and floodplain water depth.For submerged vegetation, the Manning coefficient showed a trend of decreasing with the increase in the water depth and increasing with the increase in the vegetation height.The Manning coefficient for vegetated flow generally ranges from 0.03 to 0.2.The application of Equation ( 9) in the Yueyang River section of the middle and lower reaches of the Yangtze River showed that the increase in vegetation density after the Three Gorges Project led to an increase in the floodplain Manning coefficient, which increased from 0.054 to 0.057, an increase of about 5.5%, and the comprehensive Manning coefficient of a typical section increased from 0.0303 to 0.0309, an increase of about 2.0%.The results of the study can be applied in the analysis of the water level-flow relationship and water surface line calculation of a vegetated river.It is important to note that further research is needed as to whether the present method is applicable to flexibly submerged vegetation with large deformations.

Figure 1 .
Figure 1.Sketch of the submerged vegetation flow.

Figure 1 .
Figure 1.Sketch of the submerged vegetation flow.
. The R-square of Equation (9) was 0.88, the R-square of Yang and Choi's equation was −10.83, and R-square of Li et al.'s equation was 0.85.It can be seen that the model accuracy of Equation (9) was much better than Yang and Choi's equation and comparable to that of Li et al.Considering that Equation (
, x FOR PEER REVIEW 7 of 13 water depth is shown in Figure5.It can be seen that for emerged vegetation, the Manning coefficient tends to increase with the increase in the vegetation density and floodplain water depth.The Manning coefficient was about 0.032 for a floodplain water depth of 0.5 m and vegetation density of 0.01 m −1 , 0.071 for a floodplain water depth of 0.5 m and vegetation density of 0.2 m −1 , 0.056 for a floodplain water depth of 3 m and vegetation density of 0.01 m −1 , and 0.214 for a floodplain water depth of 3 m and vegetation density of 0.2 m −1 .

Figure 5 .
Figure 5. Variation in the Manning coefficient with the vegetation density and floodplain water depth for the emerged vegetation flow.

Figure 5 .Figure 6 .
Figure 5. Variation in the Manning coefficient with the vegetation density and floodplain water depth for the emerged vegetation flow.Water 2023, 15, x FOR PEER REVIEW 8 of 13

Figure 7 .
Figure 7. Variation in the Manning coefficient with vegetation height and floodplain water depth for submerged vegetation flow (vegetation density a = 1 m −1 ).

Figure 6 .
Figure 6.Variation in the Manning coefficient with the vegetation density and floodplain water depth for submerged vegetation flow (vegetation height h = 0.5 m).

Figure 6 .
Figure 6.Variation in the Manning coefficient with the vegetation density and floodplain water depth for submerged vegetation flow (vegetation height h = 0.5 m).

Figure 7 .
Figure 7. Variation in the Manning coefficient with vegetation height and floodplain water depth for submerged vegetation flow (vegetation density a = 1 m −1 ).

Figure 7 .
Figure 7. Variation in the Manning coefficient with vegetation height and floodplain water depth for submerged vegetation flow (vegetation density a = 1 m −1 ).

Figure 9 .
Figure 9. Schematic of the middle reaches of Yangtze River from Chai et al. [30].The study area is between Chenglingji and Luoshan.

Figure 9 .
Figure 9. Schematic of the middle reaches of Yangtze River from Chai et al. [30].The study area is between Chenglingji and Luoshan. .

Figure 10 .
Figure 10.Study area image and typical section location.

Figure 11 .
Figure 11.Schematic of the main channel and floodplain zoning of the typical cross section.

Figure 10 . 13 Figure 10 .
Figure 10.Study area image and typical section location.

Figure 11 .
Figure 11.Schematic of the main channel and floodplain zoning of the typical cross section.

Figure 11 .
Figure 11.Schematic of the main channel and floodplain zoning of the typical cross section.on the channel river integrated roughness coefficient (i.e., the Manning coefficient) calculation method:n = P −1 i=N ∑ i=1

Table 1 .
Summary of experiments conducted by previous researchers.