A Delft3D-FLOW model was used to simulate tidal flow in Davis pond marsh in Louisiana, USA. The study area is a freshwater marsh consisting of one main channel and floodplain. Vegetation-induced flow resistance greatly influences tidal flow dynamics in the marsh. This study evaluated eight approaches to estimate vegetation roughness, including two constant Manning's *n* values, four empirical relations for calculating *n*, and two methods for calculating Chezy's *C* values originally embedded in the Delft3D model. Simulated results of water surface elevation (WSE) were compared with the corresponding field observation at eleven stream gauges in the study area. We concluded that the roughness coefficient for vegetated area varies with time as flow depth changes. Among the selected empirical relations for the vegetation roughness, the ones accounting for the effect of the vegetation frontal area and the degree of submergence have closely matched the measurements.

## Notation

The following symbols are used in this paper:

*a*the frontal area of vegetation stem per unit volume (m

^{−1});*B*surface width (m);

*B*,^{V}*B*^{X}volumetric and cross-section blockage factors, respectively (−);

*C*Chezy's coefficient (m

^{1/2}/s);*C*_{b}Chezy's coefficient for non-vegetated channel (m

^{1/2}/s);*C*_{D}drag coefficient (−);

*C**coefficient for the shear stress at the interface between vegetated and non-vegetated zones (−);

*g*gravity acceleration (m/s

^{2});*H*flow depth (m);

*h*_{v}height of vegetation (m);

*m*vegetation density (stem/m

^{2});*n*Manning's coefficient ();

*n*,_{o}*n*_{1}the maximum Manning's coefficient for non-vegetated channels ();

*n*_{v}n_{2}the maximum Manning's coefficient for vegetated channels ();

*N*total number of observed WSE;

*O*_{i}ith value of observed WSE;

average value of the observed WSE;

*S*_{i}ith value of the simulated WSE;

*U*_{vo}average velocity in the vegetation zone (m/s);

*V*averaged velocity across the flow depth (m/s);

*w*vegetation frontal width (m);

*z*_{0}length scale of bed roughness (m);

coefficient (−);

- β
vegetation resistance parameter ()

0.4 (for sparse, low density vegetation H > 0.3 m);

1.6 (for moderately dense vegetation H = 0.3 m);

6.4 (for very dense vegetation H < 0.3 m); and

- η
turbulence length scale (m);

von Karman constant 0.41 (−);

parameter relating to vegetation resistance;

density of water (kg/m

^{3}).

## INTRODUCTION

*h*and its frontal width is

_{v},*w*. Flow depth is H, and surface width is B. The percentage of vegetation occupied cross-section area is the cross sectional blockage factor. The product of vegetation density in stems per unit area and its frontal width is defined as

*α*. It is worth mentioning that Fisher (1992) first examined the differences of using surface area and volume to quantify the blockage factor, and recommended the percentage of surface area for the blockage factor because it is easier to measure. Recently, Green (2005), Nikora

*et al.*(2008), and Luhar & Nepf (2013) used the percentage of cross-sectional area occupied by vegetation as the blockage factor. Since the cross-section based blockage factor better represents the intensity of turbulence wake in vegetation zones, it is more closely correlated with the resistance than the surface area based blockage factor.

Early researchers (e.g. Ree & Palmer 1949; Cowan 1956; Chow 1959; Petryk & Bosmajian 1975) formulated many empirical relations for Manning's roughness coefficient correlating with flow and vegetation properties. However, these relations are only applicable to limited vegetation types and experimental flow conditions (e.g. flow depth) (Klopstra *et al.* 1997; Wu *et al.* 1999; Abood *et al.* 2006). SED2D model, a depth averaged two-dimensional hydrodynamic model, has included an empirical formula to calculate Manning's *n* for channels with vegetation. The *n* value in SED2D model is a function of the maximum *n* value for channel free of vegetation, vegetation height, flow depth, and the maximum *n* value of vegetated water (Letter *et al.* 2011). Through measuring turbulence around vegetation, researchers (e.g. Nepf *et al.* 1997; Nepf 1999) found bed shear stress affected near-bed turbulence production, while far away from the bed, the turbulence wake generated by stems becomes dominated (Nepf *et al.* 1997). The vortex from the stem wake extracts energy from the mean flow, and feeds it into turbulence kinetic energy (Nepf 1999). Consequently, turbulence intensity initially increases with stem density, and then reaches a peak at a given vegetation density, but eventually subsides as the density becomes very high (Nepf *et al.* 1997).

However, to date, there is no consensus on which method is the most feasible for calculating Manning's *n* in freshwater marshes. Vegetation-induced resistance is highly variable in the field due to complex vegetation types, density, and height. To quantify these varieties, we selected the Delft3D-Flow, a widely used open source three-dimensional hydrodynamic model, to study vegetation impact on tidal flow hydrodynamics.

## FLOW MODEL

The Delft3D-FLOW open source program (http://oss.deltares.nl/web/delft3d/source-code) is a three-dimensional (3D) hydrodynamic and sediment transport model capable of simulating unsteady incompressible flow (Deltares 2011). The Reynolds Averaged Navier–Stokes equations, under shallow water (hydrostatic pressure) assumption, were solved using a finite difference scheme on a structured staggered curvilinear grid (Stelling & van Kester 1994). Delft3D-FLOW allows users to choose from two vertical grid systems: *σ*-grid and Z-grid, and four turbulence closure models: constant eddy viscosity coefficient, algebraic eddy viscosity model, *k-L* model, and *k-ɛ* model.

Delft3D-FLOW has a function, Trachytopes, for users to define bed and flow resistance on each sub-grid (Deltares 2011). Three classes are available in the Trachytopes function: area, line, and point class. The area class has three types: the first is a constant coefficient for bed roughness, such as White Colebrook, Chezy, and Manning's coefficients, the second accounts for form resistance resulting from sand dunes, and the third is for calculating roughness coefficient in vegetated channels. In a simulation run, the roughness coefficient often remains a constant with time if choosing the first type of area class, while for the second type, it is determined by dune height, and the third by vegetation properties. The line class of Trachytopes function is used to approximate flow resistance for elements with hedge, bridge piers, and other structures. The point class is used to represent a set of point flow resistance elements, such as groups of individual trees or small plant.

Delft3D-Flow has incorporated vegetation effect through an adjusted bed roughness. For example, an implementation based on Klopstra *et al.* (1996) was added to the Klopstra *et al.* (1997) equation for calculating the *C* value of emergent vegetation (Deltares 2011) (Table 2). Additionally, an artificial term was added to account for extra momentum loss due to vegetation using Baptist's (2005) equation (Table 2). Furthermore, the momentum equations and the *k-ε* turbulent closure model in Delft3D-FLOW have been modified to include vegetation-induced momentum loss as well as influences of vegetation on turbulence generation and dissipation.

Temmerman *et al.* (2005) studied the vegetation impact on flow hydrodynamic and sedimentation processes in a tidal creek within the Paulina salt marsh in the Scheldt estuary, located at the southwest of Netherland, by modifying the momentum and *k-ɛ* equations in Delft3D-FLOW. Lately, Horstman *et al.* (2013) used Delft3D-FLOW for modeling the tidal dynamics in the mangrove forest in Trang Province, Thailand. Both cases are in saltwater marshes where tidal flow is dominant, and vegetation types are distinct from ones in freshwater marshes. None of them have compared the model's performances with results using other methods not embedded in Deflt3D to improve their modelling results.

The objective of this study is to evaluate the accuracy of different methods for calculating vegetation-induced roughness. The methods are not limited to the ones programmed in Delft3D, but all the methods available in literature. Therefore, the first and third types of area class in the Trachytopes function were adopted for incorporating other methods to calculate vegetation-induced roughness. Besides, we selected the orthogonal curvilinear coordinate system, and the vertical plane is *σ*-grid. The *k-ɛ* model was chosen to determine the coefficient of eddy viscosity. The sensitivity of modeling results to the selection of vertical grid and key parameters (e.g. *a* value) was analyzed.

## DESCRIPTION OF STUDY SITE

*et al.*2008). This causes an intrusion of saltwater that threatens the existing freshwater habitat. To reduce the effects from saltwater intrusion, the Davis pond project was constructed to divert freshwater from the Mississippi River into the adjacent estuarine areas. When the project began to operate, unexpected high water levels were observed throughout the study area. To manage the water level, McAlpin

*et al.*(2008) studied the causes by simulating the study site applying RMA2 model. McAlpin

*et al.*(2008) validated the model by comparing the simulated water surface elevations (WSEs) with the observed, and proposed 12 alternatives to reduce water level.

*Panicum hemitomon*(Sasser

*et al.*2008). A previous study (O'Neil 1949) found that

*Panicum hemitomon*is typical in freshwater marshes in the Mississippi River delta. The mature

*Panicum hemitomon*has a stem height about 0.762 m, and a stem diameter of about 1–6 mm (Turner 1994). Its leaf blades are narrow and long, about 1.25 cm wide and ranging from 20 to 30 cm in length. The leaves are rough on the upper side, smooth on the lower side, and grow along the stems (Leithead

*et al.*1971) (Figure 3). In maidencane (

*Panicum hemitomon*) marshes, the averaged stem density is 255 stems per square meter; in deep water marshes, the stem density is only 18 per square meter; and in mixed shallow marshes (

*Eleocharis elongata*,

*Sagittaria lancifolia*,

*Panicum hemitomon*, and

*Pontederia cordata*), the average density is 286 stems per square meter (Turner 1996). The Davis pond marsh is a typical maidencane marsh. Manning's

*n*in the subtropical marsh is 0.26–0.55 for water hyacinth, and 0.16–0.43 for vegetation of mixed species for flow depth about 40–65 cm (Shih & Rahi 1981). The Manning's roughness is inversely calculated by using measured flow depth and velocity in the vegetated area.

## VEGETATION IMPACT

*n*value, or a time-varying

*n*or Chezy's

*C*coefficients for each sub-area. The study area was divided into several sub-areas according to the existence of channels, overbanks, and vegetation height (Figure 4), similar to that in McAlpin

*et al.*(2008). The computational grid is overlaid on the subareas in Figure 4. Table 1 summarizes

*n*values used in the first approach, in which

*n*

_{1}and

*n*

_{2}are the maximum values of Manning's roughness coefficient for unvegetated and vegetated sub-areas, respectively (McAlpin

*et al.*2008). Two options were used to vary

*n*value with time. The equations by Baptist (2005) and Klopstra

*et al.*(1997) in Table 2, already incorporated in Delft3D under the Trachytopes function, were used in the first option to calculate the time-dependent

*C*value. This option treated an individual grass or vegetation as a cylinder with a frontal width equal to its diameter. In this study, the average frontal width was calculated by the number of leaves, the width of individual leaves, and the stem diameter. In the second option, the rest of the equations in Table 2 were programmed into Delft3D-FLOW program using the Trachytopes function. Because of the lack of detailed vegetation distribution map in the study area, we assumed each sub-area is entirely filled with vegetation of uniform height. With this assumption, the volumetric and cross-sectional blockage factor is equal to the ratio of vegetation height to flow depth.

Sub-area . | n_{1}
. | n_{2}
. |
---|---|---|

1 | 0.06 | 0.10 |

2 | 0.06 | 0.40 |

3 | 0.06 | 0.40 |

4 | 0.12 | 0.82 |

5 | 0.12 | 0.82 |

6 | 0.06 | 0.40 |

7 | 0.12 | 0.82 |

Sub-area . | n_{1}
. | n_{2}
. |
---|---|---|

1 | 0.06 | 0.10 |

2 | 0.06 | 0.40 |

3 | 0.06 | 0.40 |

4 | 0.12 | 0.82 |

5 | 0.12 | 0.82 |

6 | 0.06 | 0.40 |

7 | 0.12 | 0.82 |

Authors/reference . | Equations . |
---|---|

Fisher (1992) | |

Reed et al. (1995) | |

SED2D Model | |

Luhar & Nepf (2013) | |

Klopstra et al. (1997) | |

Baptist (2005) | |

For Baptist (2005) equation, a term will be added as an additional term in the momentum equations. For Klopstra *et al.*(1997) equation;
| |

Authors/reference . | Equations . |
---|---|

Fisher (1992) | |

Reed et al. (1995) | |

SED2D Model | |

Luhar & Nepf (2013) | |

Klopstra et al. (1997) | |

Baptist (2005) | |

For Baptist (2005) equation, a term will be added as an additional term in the momentum equations. For Klopstra *et al.*(1997) equation;
| |

## DELFT3D-FLOW MODEL FOR THE STUDY AREA

### Computational mesh

The computational grid of the study area was constructed using available geometric and bed elevation data. The bathymetry of the study area was obtained from USACE, and the elevation is referenced to North American Vertical Datum of 1988 (NAVD88). The computational mesh is a structured grid having 299 points in the main flow direction (M), and 53 points in the direction normal to main flow (N), and 15 layers in the vertical (Figure 4). A finer mesh was used in the inlet canal to capture detailed bathymetry because of rapidly varied bed elevation in this region. The maximum grid sizes in main flow and normal to main flow directions are 177.84 and 462.5 m, respectively, and the minimum are 7.13 and 19.19 m, respectively. The grid aspect ratio (N-grid size/M-grid size) ranges from 1.0 to 20.55. The regions of high aspect ratio are located along the border of the study area, where flow is predominately along the N-direction. Because the velocity gradient close to the bed is very high, the distance between *σ*-layers is smallest at the bottom, and increases toward free surface. The water depth has reached 6.5 m in the inlet canal, and 2.5 m in the other parts during the flood period from 1:00 pm (CDT) on November 30 to 7:00 pm (CDT) on December 3, 2003.

### Boundary conditions

### Results

Eleven stream gauges are available in the study area (Figure 2). The time step is one minute for both approaches to achieve numerically stable solutions. Table 3 is a summary of vegetation roughness parameters for each subarea using the roughness equations in SED2D model (McAlpin *et al.* 2008). In the equations proposed by Luhar & Nepf (2013), Klopstra *et al.* (1997) and Baptist (2005), the average frontal width of an individual plant is calculated as 2.25 times the leaf width (approximately equal to 1.25 cm) plus the mean stem diameter (*d* = 0.4 cm). Since the calculated frontal width, symbolized by *w*, is 3.2 cm, equal to 8*d*, this study assumes 8*d* as the vegetation frontal width. Since the *α* value is the product of vegetation density and its frontal width, the larger the *a* value, the larger the vegetation-induced resistance. In a typical *Panicum hemitomon* grass dominated freshwater marsh in Louisiana, vegetation density is 255 stem/m^{2}, and the average stem diameter is 0.4 cm, thus the *a*-value is 8.160 m^{–1}. This value is within the range, 1–10 m^{–1}, suggested by Luhar *et al.* (2008), Lightbody & Nepf (2006), and Leonard & Luther (1995) for marsh grasses. The drag coefficient for vegetated channel, *C _{D}*, is equal to 1.0. Shear stress at the interface between vegetated and unvegetated regions is quantified by another drag coefficient, denoted as

*C*

_{*}, and is set as 0.10 (Luhar & Nepf 2013).

Sub-area . | n
. _{o} | h
. _{v} | n
. _{v} | . |
---|---|---|---|---|

1 | 0.06 | 0.3048 | 0.10 | 0.05 |

2 | 0.06 | 0.381 | 0.40 | 0.20 |

3 | 0.06 | 0.3048 | 0.40 | 0.20 |

4 | 0.12 | 0.6096 | 0.82 | 0.41 |

5 | 0.12 | 0.762 | 0.82 | 0.41 |

6 | 0.06 | 0.381 | 0.40 | 0.20 |

7 | 0.12 | 0.6856 | 0.82 | 0.41 |

Sub-area . | n
. _{o} | h
. _{v} | n
. _{v} | . |
---|---|---|---|---|

1 | 0.06 | 0.3048 | 0.10 | 0.05 |

2 | 0.06 | 0.381 | 0.40 | 0.20 |

3 | 0.06 | 0.3048 | 0.40 | 0.20 |

4 | 0.12 | 0.6096 | 0.82 | 0.41 |

5 | 0.12 | 0.762 | 0.82 | 0.41 |

6 | 0.06 | 0.381 | 0.40 | 0.20 |

7 | 0.12 | 0.6856 | 0.82 | 0.41 |

*S*and

_{i}*O*are

_{i}*i*th simulated and observed value, respectively, is the average of observed values,

*N*is the total number of observations. The RMSE and NSE values for all the simulations were calculated at all gauges, and their ranges are shown in Figures 8 and 9, respectively. From these figures, one can find the RMSEs using the equation in SED2D model and those using Luhar & Nepf (2013), Klopstra

*et al.*(1997), and Baptist (2005) equations are relatively small, and their NSEs are closer to 1.0, which indicates that the mean square errors by using these roughness equations are much smaller than the variance of the observations, and the simulated results are more accurate than the results from using other equations.

In the first approach, Figure 8 showed that the range of RMSE for the simulated WSE using the constant Manning roughness for unvegetated channels (*n*_{1}) is smaller than the results using Manning's roughness for vegetated channels (*n*_{2}). The NSE value is also closer to 1.0 for the results using *n*_{1} (Figure 9). It can be seen from Figure 7 that WSE were underestimated when using the constant *n*_{1}, while they were overestimated when using the constant *n*_{2}. Although the results using *n*_{1} are slightly better than the ones using *n*_{2}, none of them matched the results well at all gauges. Therefore, several *n* values between *n*_{1} and *n*_{2} values were selected to re-run the model. However, for each value, the WSEs matched the observation data at one gauge in one time interval but not at other times. Also, WSEs matched the observation data better at some gauges but got worse at the rest. Regardless of the *n* value being used, there is no general improvement of modeling results. This excludes the feasibility of using a constant roughness for simulating flow hydrodynamics in a freshwater marsh.

*et al.*(1995) equations are larger (Figure 8) than the results from other equations. The corresponding NSE values using these two equations are also further away from 1.0 than the rest. The equations of Luhar & Nepf (2013), Klopstra

*et al.*(1997) and Baptist (2005) yielded the smallest RMSEs, and NSE values are also the closest to 1.0 (Figures 8 and 9). This attributes to the fact that Fisher (1992) and Reed

*et al.*(1995) equations do not take the degree of submergence into account, which is essential for differentiating the submerged and emerged vegetation. At low flow depth, the vegetation height is comparable to flow depth, or even emergent. In this condition, vegetation-induced roughness is dominant. However, as flow depth increases, vegetation becomes submerged and its influence on flow resistance is diminishing (Nepf 2012). When flow depth is orders of magnitude larger than vegetation height, vegetation-induced resistance will converge to a constant (Augustijn

*et al.*2008). Therefore, this study recommends the use of variable Manning roughness to account for vegetation-induced roughness. Figure 10 shows the results from using variable

*n*values in Luhar & Nepf (2013), Klopstra

*et al.*(1997) and Baptist (2005) equations at gauge 21. These three equations have incorporated the vegetation effect by using a vegetation blockage index (

*a*). In the next section, the sensitivity of results to modeling parameters, such as

*a*value, and vertical grid spacing, will be discussed. Nevertheless, the simulation results showed that these three equations are the best choices for simulating vegetation resistance in this estuarine marsh area.

The simulation was run on a PC with Intel (R) Xeon (R) CPU X5550 (two processors) and 32 GB RAM. The CPU times using the first approach were 74,218.84 and 73,815.20 s for the simulation using constant *n*_{1} and *n*_{2}, respectively, whereas they are 72,790.89 and 77,252.20 s using Reed *et al.* (1995) and Klopstra *et al.* (1997) equations, respectively.

## SENSITIVITY ANALYSIS

### Vegetation blockage index (*a* value)

Using the equations in Luhar & Nepf (2013), Klopstra *et al.* (1997) and Baptist (2005), results showed the best WSE matches with the observation. However, the results can be sensitive to key parameters (e.g. *a* value). To quantify the sensitivity of modeling results to the variation of *a* values, three *a* values were selected based on the average frontal width of an individual grass, equal to *d*, 11*d*, and 22*d*. Since the average vegetation density is 255 stem/m^{2} and the stem diameter is *d* = 0.4 cm, the calculated *a* values are 1.020, 11.220, and 22.440 m^{−1}, respectively.

*a*values were shown with the observed ones at all gauges in Figures 11–13. When

*a*= 1.020 m

^{–1}, the simulated WSEs were underestimated at all gauges (Figure 11), while they were overestimated, especially by the Klopstra

*et al.*(1997) equation, when

*a*= 22.440 m

^{–1}(Figure 13). The simulated WSEs using

*a*= 11.220 m

^{–1}approximately matched the observations (Figure 12). One can find, for

*a*= 1.020 m

^{–1}, the averages of RMSE were the largest, and the NSE values were the most distant from 1.0, for all three equations (Figures 14 and 15). For

*a*= 11.220 m

^{–1}, the averages of RMSE became smaller, and the NSE values were closer to 1.0. For

*a*= 11.220 and 22.440 m

^{–1}, the averages of RMSE using the Klopstra

*et al.*(1997) equation were the largest among three equations, and the NSE values were further distant from 1.0 than the other two equations. This means the Klopstra

*et al.*(1997) equation is very sensitive to the variation of

*a*. Apparently, when

*a*is equal to 1.020 m

^{−1}, none of the three equations gave accurate results of WSE. This implies that the frontal width of an individual grass should be calculated not only by stem diameter but also the number and width of leaves. It can be seen from Figures 8, 9, 14, and 15 that the smallest values of RMSE and the closest NSE values to 1.0 were found when

*a*= 8.160 or 11.220 m

^{–1}. This confirms that the average frontal width of a grass is between 8

*d*and 11

*d*in a typical maidencane marsh with

*Panicum hemitomon*as the dominant species.

### Vertical grid selection

*σ*-grid and Z-grid. The shallow water equations with hydrostatic pressure assumption were solved in both

*σ*and Z-grids. An extension for solving the non-hydrostatic pressure has been added to the Z-grid. All the computational runs conducted until now were using the

*σ*-grid assuming hydrostatic pressure. Since the

*σ*-grid forces the convective transport along the channel bottom, it may allow false fluxes between the main channel and floodplain when their elevation difference prohibits such transport. Therefore, the sensitivity of results to the vertical grid selection was performed in this section using Luhar & Nepf (2013) equation. Additional computational runs were conducted using the same computational grid but 15 Z-layers in the vertical plane assuming both hydrostatic and non-hydrostatic pressure. The time step is 0.1 minute for achieving numerically stable solutions. The CPU times for the runs using the hydrostatic and non-hydrostatic pressure assumptions were 140,968.51 and 211,734.70 s, respectively. The simulated results of WSE at all gauges were compared with the ones using

*σ*-grid and the observations in Figure 16. The results using Z-grid with non-hydrostatic pressure assumption is better than those using the hydrostatic pressure assumption. However, these results still deviate far away from the observations, and those using the

*σ*-grid. Figure 17 showed the ranges of RMSE and NSE at all gauges using

*σ*- and Z-grids, which indicated that the RMSE ranges using the

*σ*-grid are the smallest, and the NSE values are closer to 1.0 than both results using the Z-grid. The CPU time for the simulation using the

*σ*-grid is 76,613.32 s, much less than those using the Z-grid. One major reason is because the

*σ*-grid allows larger time steps (e.g. 1 min) than the Z-grid. Therefore, the Z-grid option did not yield more accurate results but more CPU times than the

*σ*-grid. This study recommends the use of the

*σ*-grid for shallow water marsh application.

## DISCUSSION

Despite significant efforts in quantifying vegetation resistance in estuarine marshes, there is no consensus on which equation is the best. Because of the varieties of vegetation, their complex stem and leaf structures, and their elastic properties in flows, vegetation-induced flow resistance still needs to be determined by using empirical relations in Table 2 or others in the literature. Although this study showed the simulated WSEs best matched the observations using the equations in Luhar & Nepf (2013), Klopstra *et al.* (1997) and Baptist (2005), this conclusion may not be valid for another marsh because vegetation species and ages can alter vegetation's mechanic properties, and thus their influences on flow resistance. Nevertheless, this study found the product of vegetation density and its frontal width is a reliable quantitative measure of vegetation property. Vegetation resistance is dependent on its submergence, and therefore varies with flow depth. Engineers need to validate a selected vegetation resistance equation using observed data of WSE, flow depth, and velocity before adopting it to any field study.

## CONCLUSIONS

This study applied the Delft3D-FLOW model to evaluate the accuracy of various methods for calculating vegetation-induced resistance in a *maidencane* marsh in the Mississippi River delta. Besides two equations embedded in Delft3D model, four other equations were programmed into the Delft3D-FLOW model. The simulated results of WSE were compared with the corresponding observation at eleven USGS gauges. Results showed that vegetation-induced roughness is a spatial and temporal variable that changes with submergence and vegetation blockage index. If treating flow roughness in each sub-area as constant, the results of WSE deviated considerably from the measurements. On the other hand, when using the equations by Luhar & Nepf (2013), Baptist (2005), and Klopstra *et al.* (1997), reasonable matches with the observed WSE were obtained because these equations have taken the degree of submergence and the vegetation blockage index into account. The best results of WSE were obtained by using a constant vegetation density for the *Panicum hemitomon* vegetation type, and the resulting *a* values ranged from 8.160 to 11.220 m^{–1}. However, this does not warrant the universal applicability of these three equations. Further research is needed to better quantify vegetation property and understand their interactions with flow field.

## ACKNOWLEDGEMENT

This project is partially funded by an Iraq Governmental Graduate Student Training Scholarship. In-kind technical and data support from the Coastal Hydraulics Laboratory in Engineering Research and Development Center, the Corps of Engineers, is greatly appreciated.