Hydraulic Features of Flow through Emergent Bending Aquatic Vegetation in the Riparian Zone

Vegetation in riparian zones has a significant influence on resistance, velocity distribution and turbulence intensity. This study experimentally investigated the effect of emergent bending riparian zone vegetation on the flow. The results showed that the frond and stem parts of Acorus calami had different influences on hydraulic features and that the relative depth ratio of water depth h to stem height h s was a key determinant of those influences. Manning coefficient n varied greatly with the variation of vegetation densities, relative depth ratio of water depth h to stem height h s , Re and Fr. Manning coefficient n increased with increasing vegetation density, particularly in cases when h/h s > 1. The velocity distributions did not follow logarithmic profiles, but they instead exhibited double logarithmic profiles. In addition, vegetation characteristics were shown to influence the height of maximum velocity. The position of maximum velocity is further away from the bed in cases with denser vegetation distribution. Finally, turbulence intensity showed more significant variation in the stem part and peaked near the middle of the stem, at z/h s = 0.5, where z was the distance from the bottom.


Introduction
A riparian zone has a complex ecotone; it is a transitional semiterrestrial area regularly influenced by fresh water that usually extends from the edges of water bodies to the edges of upland communities [1].Riparian vegetation is a main focus of river restoration schemes, and the importance of its preservation to river ecology has also been noted [2].Riparian vegetation has many ecological, aesthetic and economic benefits, such as providing terrestrial wildlife habitat; improving water quality; stabilizing streambanks and floodplains; supplying energy subsidies for aquatic and terrestrial ecosystems [3]; and serving as an important buffer in a healthy river ecosystem.Riparian zone revegetation has also been identified in recent years as one of the most common river restoration practices to obtain a more natural state [4,5].For instance, Acorus calamus, which flourishes in riparian zones, has been extensively used in river restoration projects in China.Vegetation preservation is significant in the ecology of natural and artificial systems [6].However, from an engineering point of view, riparian vegetation practices have not been often encouraged because of increased flow resistance, sediment transport effects and decreased flood discharge efficiency compared to unvegetated regions [7].Thus, a better understanding of the physical processes governing flow resistances in vegetated areas could help resolve these conflicting engineering and ecological considerations [8].
With growing interest in river restoration and flood management requirements, several studies on the hydraulic properties of riverine vegetation have been conducted.Stephan and Gutknecht (2002) reported that Gessner who emphasized the importance of the alteration of flow velocity due to macrophyte communities was one of the first biologists [9].Since then, many laboratory, numerical, analytical and field studies have been conducted to account for the complex dynamic interactions between vegetation and moving fluids [5,10].For example, Kouwen et al. [11,12] studied the interaction between flexible vegetation and flow and reported that vegetation increases flow resistance, changes backwater profiles, and modifies sediment transport and deposition [11,12].López and García (2001) proposed a κ-ε model to compute the mean velocity profile and turbulence characteristics in open channel vegetated flows [6].Defina and Bixio [13] explored the capability of two different mathematical models to predict fully developed one-dimensional open channel flow in the presence of rigid, complex-shaped vegetation with submerged or emergent leaves.
Many previous studies concentrated on determining roughness coefficients to develop design methods [2].In terms of submerged degree, aquatic plants in riparian zones fall into two categories: submerged plants and non-submerged (emergent) plants.Accordingly, flow resistance problems are usually classified into two groups: flow over submerged, short vegetation and flow through non-submerged (emergent), tall vegetation [14].Most efforts to study vegetal resistance have concentrated on submerged plants and rigid roughness [14].For simplification, plants have been simulated by a group of cylinders of the same height and diameter with regular spacing.Stone [15] conducted laboratory experiments to explore the hydraulic flow in an open channel with cylindrical roughness.Wilson, et al. [2] performed laboratory experiments on the effects of two forms (rods or stipes) of submerged flexible vegetation on the turbulence structure within the canopy and above the canopy region based on stiffness calculations.These studies would seem to characterize the flow regardless of the vegetation's rigidity [16,17].However, in terms of stiffness, aquatic plants are commonly classified into three groups: flexible, rigid and semi-flexible.Flexible and rigid plants have been more extensively studied, while few investigations have explored the mechanical response of semi-flexible plants.To address the new objectives of river restoration and environmental flood management, a better understanding of hydraulic flow features through emergent semi-flexible aquatic plants in riparian zones is required.
In this study, we focus on hydraulic flow features through Acorus calami, a typical emergent semi-flexible macrophyte, which can be divided into frond and stem.For simplification, an artificial A. calamus is used to conduct the experiments in a flume.The frond and stem have different influences on drag force, mean velocity profile and turbulence structure.The primary objective is to investigate the features of Manning coefficient n, mean velocity profiles and turbulence structures in the frond and stem of the emergent bending macrophytes.

Materials
Acorus calamus is a tall perennial monocot of the Acoraceae family and is prevalent in riparian zones and in sunny locations.Mature Acorus calami are, on average, 50-120 cm high.Its leaves are sword shaped and are from 55 to 80 cm long and between 0.7 and 1.7 cm wide.In natural riparian zones, emergent vegetations exhibit three different primary motions under the impact of flow: erecting, compressing and bending [18].A. calami stand rigid in water with a low flow velocity.As the discharge or water depth increases, the drag force on the plants increases as well, so the stems cannot withstand the flow and start to bend.This serves as a self-protection reaction, in that the vegetal element decreases its frontal area to decrease the drag force [18,19].Individual natural A. calami have particular properties (e.g., height, stiffness), which hinder the generalization of the hydraulic features of natural A. calami.To investigate the hydraulic flow features of A. calami, it is necessary to use plant forms that are simplified and generalized compared to the complexity of real vegetation assemblages [2].This is a similar approach to that adopted by other authors (e.g., Nepf and Vivoni 1999;Wilson, et al. 2003) [2,17].In the present study, artificial plants manufactured from PE plastic were used to simulate natural A. calami.All parameters, including geometric, physical and kinematic parameters, were scaled using Froude's law.In general, the dimensional parameters are 1:1 scale replicas of A. calami.The simulated A. calami are from 0.55 m to 0.7 m high and have 6 leaves from approximately 0.005 m to 0.018 m wide and between 0.65 m and 0.75 m long.They have basal stems with a mean diameter of 0.0115 m and mean height of 0.0612 m.The geometry and notation of simulated A. calami are shown in Figure 1.In addition, the bending stiffness (defined as the gradient of the force-deflection curves) of real and simulated A. calami are evaluated in the laboratory using the approach of Wilson (2003) [2].The bending stiffnesses of natural and manufactured A. calamus are 14.35 N/m and 13.42 N/m, respectively.All geometric and physical parameters of natural and manufactured A. calami satisfy the Froude similarity rule.

Experimental Setup and Measurement Technique
The experiments were conducted in a 26-m-long, 0.5-m-wide and 0.7-m-deep glass-walled flume (Figure 2).The slope of the flume is fixed at 0.07692%.Water was pumped into a stilling basin and flowed into the flume.The test region (12 m in length) is located 9 m from the entrance section.Upstream and downstream of the test segment, there were two transition segments (12 m and 5 m in length, respectively).The water surface profiles were controlled by the downstream tailgate, which was installed at the end of flume and could be adjusted using a gear system that allowed the opening size to be prescribed with a high degree of accuracy.The two transition segments and the tailgate prevented large-scale disturbances from the inlet and outlet and allowed for the development of a quasi-constant water flow by depth.The outlet and inlet of the flume were both connected to a tank, which enables the continuous recirculation of steady discharges.To ensure uniform conditions, the flow depths were controlled and measured using three water level gauges at three cross sections along the flume (10 m, 15 m, and 20 m from the entrance section), and the tailgate was adjusted as needed.Discharge was measured using the overflow weir and a water level gauge in the water-receiving channel.3D mean velocities and velocity fluctuations were measured using a 3D sideways-looking Acoustic Doppler Velocimeter (ADV) manufactured by SonTek Inc., (San Diego, CA, USA).The ADV was mounted in a wood frame across the center section of the test segment, and it could be easily moved upstream or downstream, so that all sampling points were vertically aligned.The ADV technology is based on the pulse-to-pulse coherent measuring method.The instrument consists of three modules: a measuring probe, a conditioning module, and a processing module [20].Mean velocities and velocity fluctuations can be measured from a small sampling volume (approximately 0.125 cm 3 ), which is 5 cm away from the sensing elements.The probe was connected to the computer via a processing card.Real-time data could be recorded using a data acquisition program installed in the sampling computer.For each measurement, the sample-reporting rate was 25 Hz, and the acoustic frequency was 10 MHz.The velocity components, u, v, and w, correspond to the stream-wise (x), lateral (y), and vertical (z) directions, respectively.A 180-second sampling period was set for each test run.The WinADV-program, a post-processing program, was used to filter and post-process the sampled data.Data with average correlations less than or equal to 70% were filtered out.

Test Series Description
Plants were placed into a PVC board platform with dimensions of 12 m × 0.5 m × 0.02 m (length, width and thickness in the principal flow direction, respectively) with prepared holes.The platform was fitted onto the bed of the test segment, 9 m away from the flume entrance segment (shown in Figure 2).The plants were arranged in a staggered pattern, with spacing of 15 cm, 30 cm and 45 cm, by sticking the end of the plant into the drilled holes.The arrays of plant arrangements are shown as V1, V2, V3 and V4 in Figure 3  Depending on the vegetation set-up, 8-9 h-Q combinations have been tested for each vegetation distribution pattern to investigate the influence of the ratio of water depth h on the stem height h s and h/h s , as well as the influence of the vegetation density, δ, on Manning coefficient n and the velocity distribution.Data for the 5 test series, with a total of 44 test runs, have been obtained in this study.A summary of the test series is presented in Table 1.In all the test series, the Reynolds numbers, Re, ranged from approximately 5600-48,000, indicating that all of the test runs are within the range of turbulent flow.The Froude numbers, Fr, ranged from 0.070 to 0.604, which means that all test runs can be considered subcritical flow.The plants exhibited no bending but instead oscillated slightly with low discharge.Conversely, because the leaves are flexible, the stems bend somewhat, and the leaves are streamlined with higher discharge.Furthermore, for each vegetation distribution pattern, the greatest discharge run (circled data in Table 1) has been selected to investigate flow velocity distribution and turbulence by water depth.

Vegetation Density
Most previous experiments have been conducted with arrays of rigid emergent cylindrical elements of constant diameter, where the momentum absorbing areas of the plant were constant with plant height.In these cases, the density was defined as the projected area of rod per volume [2,19].This definition can be extended to quantify the vegetation density for emergent semi-flexible plants.Vegetation density can be defined as the momentum absorbing area per unit water volume.Determining the momentum absorbing area is a topic of significant concern.In some previous studies, the total momentum absorbing area is represented by the total frontal area of plants [21,22].However, because wetted leaves spread in the flow direction in water, considerable foliage area is hidden behind the frontal area.Thus, it is important to use the total wetted surface area instead of the total frontal area in these calculations.In practice, the wetted surface area consists of the wetted stem surface area and the wetted frond surface area.When the water depth is deeper than the stem height, both the wetted stem surface area and the wetted frond surface area should be considered.Otherwise, only the stem area should be considered.Vegetation density can be described as in Equation ( 1): where δ is the vegetation density, m −1 ; N is the total number of plants; D is the diameter of stem, m; A is the area of testing reach, m 2 ; h is the water depth, m; h s is the height of stem, m; and L w is the width of wetted foliage, m.
In Equation ( 1), all parameters (aside from L w ) can be easily obtained.However, because L w mainly depends on the bending degree, it cannot be directly determined.In this study, L w is determined statistically using the geometric mean after measuring the wetted foliage widths of 20 plants for each test run.

Resistance Coefficient
Estimating the flow resistance of vegetation is of great importance in river management because of its potentially significant effects on channel conveyance [14].Efforts to quantify hydraulic roughness in vegetated channels date to the 1950s and 1960s [5].The methods for quantifying resistance vary from one-dimensional approaches (e.g., Manning coefficient or Darcy-Weisbach friction factor) to three-dimensional approaches based on the drag force equation.Järvelä (2002) calculated the Darcy-Weisbach friction coefficient based on measured head losses [14].In practice, Manning coefficient n is the most widely used resistance quantifying metric in river hydraulics.Many previous studies have presented Manning coefficient n as a function of hydraulic radius, vegetation density, vegetation stiffness, and plant frontal area.Flow resistance is dominated by vegetation resistance, and the vegetation drag coefficient is used as a fitting parameter [21].In this study, Manning coefficient n is calculated using Equation (2) proposed by Fischenich (2000), which is based on the drag force concept [23].
where R represents hydraulic radius; C d is the effective vegetal drag coefficient; A veg is the frontal area of the vegetation; and k n is the units term.

Turbulence Intensity
The quantitative analysis of turbulent flow is based on measurements of velocity fluctuations at a single point in the flow (or at multitude points, from which the flow fields can be inferred).In this study, the mean flow velocity components (u, v, and w) and velocity fluctuation components in turbulent flow (u', v', w') correspond to the streamwise, lateral and vertical directions, respectively.Velocity fluctuations can be defined as the deviation from the mean velocity.In general, root-mean-square values (RMS) of velocity fluctuations are considered measures of turbulence intensities.Turbulence intensities corresponding to streamwise, lateral and vertical directions are defined as follows [24]: 2, and where N is the total number of observations in a given series; u', v' and w' are velocity fluctuations of streamwise flow, lateral flow and vertical flow, respectively; and RMS u , RMS v and RMS w are turbulence intensities corresponding to streamwise, lateral and vertical flow, respectively.

Vegetation Density
The profiles of vegetation densities δ against water depth h for the four vegetation distribution patterns are plotted in Figure 4. Figure 4 illustrates that vegetation density increases with rising water level in all cases.However, when the water level is lower than stem height, only minimal increases in vegetation densities can be observed.Conversely, much more significant increases in vegetation densities can be observed at greater h.Furthermore, the V4 cases show the largest increasing rate of vegetation density for h/h s > 1, while the V1 cases have the least increase rate of δ.

Manning Coefficient n
In this study, Manning coefficient n is used as a measure of flow resistance.The total resistance of the flume includes sidewall and vegetation bottom resistance, designated as n w and n b , respectively.Because the flume wall is glass, its resistance can be considered negligible compared to that of vegetation.In other words, the flow resistance in a vegetated flume is dominated by the vegetation bottom resistance.The flow resistance n = n b can be evaluated by applying Equation (2).After investigating various test runs, it may be concluded how Manning coefficient n varies with the relative depth ratio of water depth h to stem height h s , the vegetation density and the Reynolds number.These are shown in Figure 5a-d.

Manning coefficient n h/h s
Figure 5 shows that the Manning coefficient n of vegetation flow (e.g., V1, V2, V3 and V4) is greater than that of non-vegetation flow (e.g., V0).Additionally, Figure 5 illustrates that the Manning coefficient has great variation with the distance between two plants.When the distance between two plants is less under the same discharge situation, the Manning n becomes greater.The sequence is V4 > V3 > V2 > V1.It is similar to Hall's results [25].Furthermore, Figure 5a shows that the Manning coefficient n increases with increasing the ratio of water depth to plant stem height h/h s .However, the increment of change in n below the stem height (h/h s = 1.0) is less than that above the stem height.When the ratio of flow depth to plant stem height, h/h s < 1, the n increment is approximately 3%.Conversely, when h/h s > 1, the Manning n increments are approximately 5.7%, 7.7%, 16.4%, 23.2% at V1, V2, V3 and V4, respectively.It illustrates that the contribution of foliage with flexible leaves is significantly greater than that of plant stems without leaves.Thus, variations in Manning n depend on the position of the separating point between the stem and foliage.The leaf is found to be an important factor in the analyses, and properties such as rigidity are found to influence the drag exerted by leaves.This can also be illustrated by Figure 5b.When the vegetation density is less than approximately 6.5 m −1 , an increasing trend of Manning n is not apparent.However, when the vegetation density is greater than 6.5 m −1 , the Manning n increases more significantly with increasing vegetation density.Manning coefficient is almost uniform up to the critical plant density and then increases [6].
Figure 5c also shows that the general tendency is an increase in the n value with increasing Reynolds number Re.Additionally, in high dense distribution (V3 and V4), the increase tendency is very dramatical when Re > 10,000.It implies that Re = 10,000 would be a critical Renolds number.When the Re < 10,000, the flow would be laminar flow.Otherwise, it is transit flow or turbulent flow.Figure 5c shows that the Renolds number Re the distance between two plants Figure 5d indicates that the Manning n decreases with increasing Fr.Therefore, the separating point and water depth are the definitive factors which have most important influence on the Manning coefficient n.

Velocity Distribution
Velocities at any point can be measured in three mutually perpendicular directions: streamwise (parallel to the boundary, u), lateral (normal to the boundary, v) and vertical (w).The test runs V1-9, V2-9, V3-9 and V4-9, which are the largest discharge runs for different vegetation density patterns (0.02391 m 3 /s for V1, 0.02056 m 3 /s for V2, 0.01290 m 3 /s for V3 and 0.01090 m 3 /s for V4), are used to investigate velocities and turbulence intensities in x, y and z directions.The sampling vertical line stands at the center of four neighboring plants.The surveyed velocity profiles are shown on the plane z/h s against u/u* (Figure 6), where z is the distance from the bed, h s is the plant stem height and u* is the shear velocity estimated using ghJ u = * , where J is the flume slope [26].
For all experimental test runs with Re > 12,500 and Fr < 1, vegetation flows are considered to be turbulent subcritical flows.Figure 6 shows that the mean flow with plants is greatly slowed and no longer follows a logarithmic vertical profile.From the velocity profiles, two zones can be distinguished: the stem part (z/h s < 1) and the foliage part (z/h s > 1).However, because both profiles in the stem and foliage parts almost follow logarithmic behavior, the whole velocity profile can be characterized as a double logarithmic profile.Furthermore, both the stem and foliage parts have significant variations in mean velocities and exhibit the generation of a horizontal shear layer.The magnitude of the mean velocity in the foliage part is less than that in the stem part because the foliage consumes much more energy and momentum from the flow through the generation of turbulence.A local velocity maximum can be observed near the bed in the stem part.In addition, Figure 6 shows that vegetation densities have some influence on the vertical velocity distribution.The vertical velocities of V3 and V4 are much lower than those of V1 and V2, as expected.Figure 6 also shows that the position of maximum velocity is farther from the bed in the denser vegetation distributions.This result is similar to the report by Afzalimehr, et al. [27].

Turbulence
In each vegetation distribution pattern, the largest discharge run was investigated to measure the turbulence intensities and Reynolds stress ( ' 'w u , instantaneous streamwise velocity fluctuations multiplied by the instantaneous vertical velocity fluctuations, then time averaged).For example, the turbulence intensities (RMS u , RMS v , RMS w ) and Reynolds stress ( ' 'w u ) for V3 (Q = 0.0129 m 3 /s) and V4 (Q = 0.0109 m 3 /s) are shown in Figure 7a,b, respectively.).This finding indicates that the turbulent transport of momentum in the vertical direction is so small that the Reynolds stress can be neglected.
Figure 7 also shows that the turbulence intensity distribution is anisotropic.In general, turbulence intensity in the z direction is much smaller than that in the other two directions.As z/h s > 1, turbulence intensity has little variation and can be considered constant.However, the turbulence intensity varies dramatically for 0 < z/h s < 1.The turbulence intensity peaks near the middle of the stem.In another words, turbulence intensity increases for z/h s < 0.5 and decreases for 0.5 < z/h s < 1. Afzalimhr, et al. [27] reported similar results of approximately constant turbulence intensity above the stem height.When the water depth is higher than the stem height, the flexible leaves bend over in the flow when the discharge is high.As the leaves become parallel to the flow, the momentum flux is absorbed, decreasing the turbulence intensity of the flow within the canopy.While the foliage induces a greater drag force, the turbulence generated by shear is reduced due to the inhibition of momentum exchange by the plant frontal area.In conclusion, the stem and foliage of emergent bending aquatic vegetation have quite different effects on the turbulence intensity of the flow, with more intense turbulence in the stem regime.

Conclusions
This paper presented an experimental study that explored the effect of emergent bending vegetations on the Manning coefficient n, the velocity distribution and turbulence structure within the stem and foliage for various distribution combinations.In this study, simulated A. calami, 1:1 scaled with a real plant, were arranged into four density patterns in a laboratory flume.An ADV was used to measure the mean velocities and velocity fluctuations.
Compared with non-vegetation flow, vegetation flow in riparian zones was characterized as having larger Manning coefficients n, a non-log law of mean velocity profiles and obvious turbulence.Furthermore, these features depended significantly on the vegetation properties, densities and distribution patterns.Particularly, because the frond and stem parts of A. calami had different influences on hydraulic features, the ratio of water depth to stem height h/h s was a key factor for investigating the hydraulic features of this species.For instance, Manning coefficient n varied greatly with the variation of vegetation densities δ, relative depth ratio h/h s , Re and Fr.Manning coefficient n increased with the increasing vegetation density and increased more rapidly in cases with h/h s > 1.
Vertical velocity profiles did not follow logarithmic profiles.According to the vertical velocity distribution, the flow could be separated into two zones: a stem zone and a foliage zone.In general, the velocity distribution showed double logarithmic profiles.The magnitude of the mean velocity in the foliage zone was less than that in the stem zone.Vegetation densities were shown to be influential on the vertical velocity distribution.The peak point of velocity was further away from the bed in cases with denser vegetation distributions.Therefore, vegetation played an important role in the boundary layer thickness and the height of maximum velocity.
For emergent vegetation, the Reynolds stresses within vegetations were quite small, and the turbulence intensities were found to be anisotropic.Turbulence intensities peaked at around z/h s = 0.5.This result indicates the presence of strong turbulence intensity around the stem and that bending flexible foliage weakened the turbulence intensities.
This work focuses on a single type of aquatic plant in riparian zone.However, riparian zones actually include erecting, bending and submerged plants.Further studies should investigate the interaction between flow and different vegetation types and combinations.Field surveys and additional laboratory experiments are planned.

Figure 1 .
Figure 1.Geometry and notation of simulated A. calami in flow.
. V1 represents one or two plants in a row, with a spacing of 45 cm and a total of 40 plants in the test region.V2 shows two plants in a row, with a spacing of 30 cm and a total of 60 plants in the test region.V3 indicates three plants in a row, with a spacing of 15 cm and a total of 240 plants in the test region.V4 shows four plants in a row, with a spacing of 15 cm and a total of 280 plants in the test region.In addition, we explored the flow features of a non-vegetation flume, which is symbolized by V0.

Figure 4 .
Figure 4. Profiles of vegetation densities against water depth in different vegetation distribution patterns, where h is the water depth and δ is the vegetation density.

Figure 5 .
Figure 5. Manning coefficient n for each series: (a) n against relative depth ratio, h/h s , where h is the water depth and h s is the plant stem height; (b) n against vegetation density, δ(m −1 ); (c) n against Reynolds number, Re; (d) n against Froude number, Fr.

Figure 7 .b
Figure 7. Turbulence intensities and Reynolds stress against relative depth ratio of vertical distance from bed z to stem height h s : (a) represents V3 pattern and Q = 0.0128 m 3 /s; (b) represents V4 pattern and Q = 0.0109 m 3 /s.

Table 1 .
Summary of experiments for vegetation density series.