Vegetation exists in many natural rivers and woodlands. No matter what type the vegetation is, it will alter the velocity field of flow, consequently increasing figflow resistance. The influence of vegetation on the velocity of flow depends on its type (rigid or flexible) and condition (submerged or emergent). As a prerequisite for analysis of flow resistance, pollutant mixing process, etc., the velocity profile of vegetated flows has drawn the interest of many researchers. For example, Temple (1986) studied velocity distribution in channels with grass; turbulent structure of vegetated layer was studied (Shimzu & Tsujimoto, 1994; Nepf & Koch, 1999; Nepf & Vivoni, 2000; Lopez & Garcia, 2001; Dimitris & Panayotis, 2011 ). There are extensive studies on the velocity distribution via experimental and analytical methods (Tsujimoto & Kitamur 1990, Klopstra et al., 1997; Meijer & Van Velzen, 1999; Nepf & Koch, 1999; Ghisalberti & Nepf, 2004; Defina & Bixio, 2005; Baptist et al., 2007; Kubrak, et al. 2008; Huai et al. 2009; Yang & Choi, 2010; Nepf, 2012; Tang, 2018) .
However, due to different types of vegetation and flow conditions, various methods are proposed for predicting velocity profiles. For the flows with submerged vegetation, the most common method is a two-layer approach, in which different analytical models are applied in the lower vegetation layer and the upper surface layer, using a mixing-layer analogy (Klopstra et al., 1997; Meijer & Van Velzen, 1999; Defina & Bixio, 2005; Baptist, et al., 2007; Huai, et al. 2009; Yang & Choi, 2010; Nepf, 2012; Singh et al. 2019; Tang, 2018) . Klopstra et al. (1997) proposed a two-layer model of velocity: one in the vegetation layer and one above it called the surface layer. In the vegetation layer, the turbulent eddy viscosity is expressed by the product of a characteristic length (λ) and velocity, following Boussinesq hypothesis. λ was found to empirically relate to H/h (H is the flow depth and h denotes the height of vegetation). After further investigation on this model through more data, Meijer & Van Velzen (1999) found that λ could be approximated by . Similarly, Defina & Bixio (2005) established a different form of analytical solution for velocity in the vegetation layer. Yang & Choi (2010) and Nepf (2012) proposed different empirical models for velocity in the vegetated layer.
Considering only certain data used in the validation of the above-mentioned models, the present study aims at comparing the four models of Klopstra et al. (1997) , Defina & Bixio (2005) , Yang & Choi (2010) and Nepf (2012) for the prediction of velocity profiles in submerged rigid vegetation against a wide range of independent datasets. Thus it is to demonstrate the capability and robust of each model, consequently making some recommendation for the future application. In this study, the 19 datasets used cover different submergence ratios (H/h) ranging from 1.25 to 3.4, different vegetation densities (defined as a, the frontal area of the vegetation per unit volume) (a = 1.1 - 10 m−1) and bed slopes So (1.8 × 10−6 - 4.0 × 10−3).
2. Description of Models
In open-channel flow with vegetation, the bed and wall boundary stress are both negligible compared with the drag force on the vegetation (Nepf & Vivoni, 2000; Stone & Shen, 2002) . Thus, the momentum equation of steady fully-developed 1-D vegetated flow may be described as (Klopstra et al., 1997; Ghisalberti & Nepf, 2004; Defina & Bixio, 2005; Kubrak et al., 2008; Tang, 2018) :
where τ is the shear stress, g the gravity, z the vertical coordinate above the bed, So the bed slope and Fv is the drag force per unit mass generated by the vegetation, see Figure 1. The drag force Fv is given by:
where h is the height of vegetation, u the local stream wise velocity, CD the drag coefficient, a the frontal area of vegetation (Av) per unit volume, representing the density of vegetation, and m is the number of vegetation per unit area.
where νT is the total eddy viscosity of vegetated flow, and λ is a mixing length.
Under steady flow condition, inserting Equations (2) & (3) into (1) gives:
For given a and CD, one can obtain an analytical solution for u2 in Equation (4) under appropriate boundary conditions (e.g. Klopstra et al., 1997; Defina & Bixio, 2005 ). The boundary conditions are as follows:
At the bed (z = 0), where the bed shear stress is neglected compared with drag force of vegetation, the local equilibrium between vegetation drag and gravity force leads to:
At the top of the vegetation (z = h), the Reynolds stress is described as:
Figure 1. Velocity profile in a channel with submerged vegetation.
2.1. Model by Klopstra et al. (1997)
For the vegetation layer, Equation (4) has the following analytical solution for the velocity under boundary conditions (5) & (6):
For the surface layer, the velocity is described by the well-known logarithmic profile:
where κ is von Karman’s constant (0.40), zm (=h − hs) is the zero-plane displacement of the logarithmic profile, hs is the distance between the top of vegetation and the virtual bed of the surface layer (see Figure 1), zo is the equivalent height of bed roughness, and u* is the shear velocity, given by:
2.2. Model by Defina and Bixio (2005)
Based on Equation (4) with similar assumptions, and using the mixing length concept, Defina & Bixio (2005) proposed another form of velocity profile by introducing a parameter β, which depends on λ/h. The velocity profile is described as follows:
For the vegetation layer:
where λ is recommended to be .
For the surface layer, the velocity is described by the same Equation (11), where hs and zo are respectively given by:
2.3. Model by Yang and Choi (2010)
Yang and Choi (2010) proposed a two-layer model that predicts the velocity profile for flows with both flexible and rigid submerged vegetation. The flow through the vegetation is assumed to be driven by gravity and the shear force generated at the top of the vegetation. This is a result of the co-flowing between the low speed vegetation layer and relatively higher speed surface layer, together with the drag force of vegetation. Thus, the flow velocity in the vegetation layer is given by:
The flow in the upper non-vegetated layer is assumed to be driven by gravity force and Reynolds stresses, leading to a logarithmic profile as:
where cu is recommended as 2 when a > 5; otherwise it is 1.
2.4. Model by Nepf (2012)
Nepf (2012) proposed a model using an empirical equation to describe the velocity profile in the vegetation layer. Although the model is supposed to be applicable to both flexible and rigid vegetation, most of its parameters are only proposed for rigid vegetation. Therefore, only the parameters for the rigid vegetation are introduced here. For the surface layer, the velocity profile is assumed as Equation (11), where zm and z0 are given by:
where z0 is recommended for ah > 0.1; otherwise, z0 is estimated as 0.1h (Nepf & Ghisalberti, 2008) .
The vegetation layer is divided into two zones, which depends on the penetration of the turbulence stresses from the surface layer into the vegetation layer. An empirical Equation (24) is proposed for estimating the depth of turbulent shear penetration (δ) for (ah) values of (0.2 - 3):
Therefore in the upper zone of the vegetation layer (h − δ < z < h), the flow is driven by gravity and turbulent stresses and balanced by vegetation drag. The velocity profile is assumed to be exponential:
where uh = u at the top of the vegetation, and u2 = velocity in the lower region of vegetation (z < h − δ), which is described by Equation (5).
3. Data for Study
To compare the four models in Section 2, we used a wide range of different experimental data from the literature for submerged rigid vegetation. Table 1
Table 1. The dataset used for evaluating the models of submerged rigid vegetation.
summarizes a total of 19 datasets used in this study. These datasets include different submergence ratio, i.e. H/h = 1.25 - 3.4, various vegetation densities a = 1.1 - 10 m−1 and bed slopes So between 1.8 × 10−6 and 4.0 × 10−3. Note that the data of Dunn et al. (1996) and runs 34 and 36 of Meijer & Van Velzen (1999) are extracted from Dimitris & Panayotis (2011) . The values of CD in the dataset of emergent vegetation were estimated for all runs, taken from the original papers.
For simplicity of analysis, the models by Klopstra et al. (1997) , Defina & Bixio (2005) , Yang & Choi (2010) and Nepf (2012) are denoted as Klopstra, Defina, Yang and Nepf models respectively in the subsequent sections.
4.1. Comparison of Results between the Models
The comparisons between the four models are shown in Figures 2-4, which demonstrate that all the models, except the Yang model, can predict profiles reasonably well for the lower zone of vegetation. However, none of them can predict the velocity profile well in the surface layer for all the data tested. For the
surface layer, both Klopstra and Yang models under-estimate, particularly the Yang model significantly under-estimates the prediction for all the data, whereas the Defina and Nepf models both over-estimate the velocity in most cases although the two models can predict reasonably well in certain cases. For example, the Defina model predicts reasonably well for the data by Huai et al. (2009) , Meijer & Van Velzen (1999) , and Shimizu & Tsujimoto (1994) , whilst the Nepf model does well for the data by Dunn et al. (1996), Huai et al. (2009) and Shimizu & Tsujimoto (1994) in Figure 3 and Figure 4. It is not surprising that the Defina model predicts well for the data of Meijer & Van Velzen (1999) as the same value of λ was used in the Defina model. It appears that both models performs reasonably well when H/h is not very large (e.g. H/h < 2.0).
Despite the Klopstra and Defina models having similar assumptions, the Defina model performs better than the Klopstra model in most cases. The reason appears from the different λ models used. Close examination shows that if the same λ model is used, the predicted velocity profiles are almost the same although z0 and zm proposed by the two models are different. Note that the Klopstra model did not consider the effect of the 1st term on the right side of Equation (18).
The Nepf model does not perform well in some cases. This is mainly due to the insufficiency of the formulas of velocity profile in the vegetation layer. The proposed formula for calculating the velocity in the upper zone of vegetation depends on two velocity predictions: one (u2) at the top of vegetation and the other (uo) in the lower part of the vegetation. Therefore, if these two velocities deviate from the measured ones significantly then the formula for the upper zone of vegetation will underperform. The proposed formula for the lower zone of vegetation neglects the dispersive, turbulent and viscous stresses and, as mentioned before, in the above experiments the vegetation is deeply penetrated by the turbulent stress and in such cases the dispersive stress might become significant. Bearing in mind that this model is based on the empirical equations for estimation of zm and zo, which is dependent on CD and vegetation density (a), further study is needed about the empirical Equations (22) & (23).
The predicted velocity by the Yang model significantly deviates from the measured ones. It can be argued that assuming a uniform distribution of velocity profile over the whole vegetation depth is not a realistic assumption. The measured velocity profile in the vegetation in most of the cases tested herein is exponential for a large part of the vegetation depth. Moreover, in some cases illustrated here, this model has predicted the velocity very close to the measured one but the trend is not well captured because the assumption of a uniform velocity does not reflect reality. The inaccuracy of the logarithmic part might be explained from the following two aspects. First, the equation for the surface layer depends on the predicted velocity in the vegetation layer. Therefore, the accuracy of the model in the surface layer is partly dependent on the precision of prediction in the vegetation layer. Second, examining the measured data indicates that the penetration of turbulent stresses into the vegetation layer is significant, so the effective roughness should be much smaller than the vegetation height, which is taken as the effective roughness in this model. Therefore, in all cases the velocity in the upper region is under-estimated.
4.2. Further Discussion on Klopstra & Defina Models
Further examination of the Defina and Klopstra models shows that their differences are very small if the same λ value is used. Because both z0 and hs are important to describe the velocity in the surface layer, they relate to the λ value in the models. Since both flow depth and vegetation height will have impact on λ (characteristic length of eddy), λ may be represented by , where k is a constant. The predicted velocity by both models decreases as increasing k value, because the larger the k value, the bigger λ becomes, i.e. the stronger the eddy, indicating a smaller velocity. Close studies on all the data tested suggest that when k = 1/40 both models have good agreement with the data, as shown in Figures 5-8.
In a river with submerged vegetation, the velocity profile in the lower vegetation layer is significantly different from the upper non-vegetation surface layer. Comparison has been made on four commonly used two-layer models of velocity profile: Klopstra et al. (1997) , Defina & Bixio (2005) , Yang & Choi (2010) and Nepf (2012) . This study has demonstrated that all the models except the Yang model can predict the velocity reasonably well in the vegetation layer near the bed. However, in the surface layer none of the models using their default parameter
values can predict velocities well for a wide range of vegetation (H/h = 1.25 - 3.4). Both the Klopstra and Yang models significantly under-estimate the velocity. The Nepf model can predict reasonably well for certain cases with lower density of vegetation, i.e. ah is up to 0.35. However, the Defina model is more capable of prediction for a relatively large range of cases.
Close examination of λ parameters in the Defina and Klopstra models shows that the two models are very close when the same parameter λ value is used. This study shows that when λ is described by with k having an optimum value of 1/40, the Defina and Klopstra models can both predict the velocity profiles well for a wide range of vegetated flows. More data are needed for examining the new recommended value of k in future research.
The author would like to acknowledge the financial support by the National Natural Science Foundation of China (11772270) and the Research Development Funding of XJTLU (RDF-15-01-10).
 Baptist, M. J., Babovic, V., Rodriguez Uthubrubu, J., Keijzer, M., Uittenbogaadr, R. E., Mynett, A., & Verwey, A. (2007). On Inducing Equations for Vegetation Resistance. Journal of Hydraulic Research, 45, 435-450.
 Dimitris, S., & Panayotis, P. (2011). Macroscopic Turbulence Models and Their Application in Turbulent Vegetated Flows. Journal of Hydraulic Engineering, 137, 315-332.
 Hao, W. L., Zhu, C. J., & Chang, X. P. (2014). Research on Vertical Distribution of Longitudinal Velocity in the Flow with Submerged Vegetation. Journal of Heibei University of Engineering, 31, 64-67. (In Chinese)
 Huai, W. X., Chen, Z., & Han, J. (2009). Mathematical Model for the Flow with Submerged and Emerged Rigid Vegetation. Journal of Hydrodynamics, 21, 722-729.
 Klopstra, D., Barneveld, H. J., Noortwijk, J. M., & Velzen, E. H. (1997). Analytical Model for Hydraulic Roughness of Submerged Vegetation. Theme A. Proceedings of 27th IAHR Congress, San Francisco: ASCE, 775-780.
 Kubrak, E., Kubrak, J., & Rowinski, P. M. (2008). Vertical Velocity Distribution through and above Submerged, Flexible Vegetation. Hydrological Sciences Journal, 53, 905-920.
 Lopez, F., & Garcia, M. H. (2001). Mean Flow and Turbulence Structure of Open-Channel Flow through Non-Emergent Vegetation. Journal of Hydraulic Engineering, 127, 392-402.
 Shimizu, Y., & Tsujimoto, T. (1994). Numerical Analysis of Turbulent Open-Channel Flow over a Vegetation Layer Using k-Turbulence Model. Journal of Hydroscience and Hydraulic Engineering, 11, 57-67.
 Singh, P., Rahimi, H., & Tang, X. (2019). Parameterization of the Modeling Variables in Velocity Analytical Solutions of Open-Channel Flows with Double-Layered Vegetation. Environmental Fluid Mechanics.
 Stone, B. M., & Shen, H. T. (2002). Hydraulic Resistance of Flow in Channels with Cylindrical Roughness. Journal of Hydraulic Engineering, 128, 500-506.
 Tang, X. (2018). A Mixing-Length-Scale-Based Analytical Model for Predicting Velocity Profiles of Open Channel Flows with Submerged Rigid Vegetation. Water and Environment Journal, 32, 1-10.
 Yang, W., & Choi, S. (2010). A Two-Layer Approach for Depth-Limited Open Channel Flows with Submerged Vegetation. Journal of Hydraulic Research, 48, 466-475.