Assessment of Wind and Vegetation Interactions in Constructed Wetlands

Meteorological data from vegetated and un-vegetated wetlands during wet and dry seasons, were collected and analyzed to evaluate the role of wind and vegetation on wetlands’ hydrology. Wind speed diminished by as much as 40%, accompanied by a measurable change in wind directions in the vegetated compared to the open water site. Wind speed and direction means were significantly different (p < 0.001 and <0.01), for vegetated and non-vegetated wetland, respectively. Cattails (Typha sp.) and open water estimates of wind drag coefficients using the log wind profile, were 0.016 and 0.009 for dry season, and 0.012 and 0.005 for wet season, respectively. Wind set up near the wetland outlet was more pronounced at shallow water depth (<20 cm). Measured velocity profile during inflow discharge event with a wind speed of 0.53 ms−1, showed two-layer flows; wind-generated surface water flow opposite to a sub-surface inflow. This opposing surface flow increases hydraulic residence time and improve nutrient uptake. Conversely, wind-generated flows aligned with inflow discharges, accelerates water flow towards the outlet, reduce the duration of water-biotic interactions, and decrease nutrient uptake.


Introduction
Aquatic vegetation is an integral part of all natural and constructed wetlands and plays a vital role in treating both urban and agricultural runoff as well as wastewater. Aquatic plants assimilate nutrients for growth and reproduction, provide surface area for microbial growth, stabilize the soil, act as a filter to trap submerged detritus and sediment, induce resistance to flow, and decrease the ability of the water to carry sediments [1][2][3][4][5].
In south Florida, permanent features such as above ground reservoirs (flow equalization basin (FEB)) and large constructed wetlands (stormwater treatment area (STA)) are used to store and remove phosphorus (P) from stormwater runoff. Presence of emergent aquatic vegetation (EAV) in wetlands in those features reduces prevailing wind speed similar to the presence of trees, forests, and obstacles on land [6]. The key role of aquatic vegetation in both features is to remove phosphorus from agriculture runoff to acceptable levels before water is released southward into the Northern Everglades. However, the presence of EAV in those features induce resistance to flow, which causes delays of water deliveries within and in-between features [7].
Previous research has identified hydraulic factors that have the greatest impact on wetland treatment efficiency, mainly the duration of water-biotic interactions [8], water depth, hydraulic loading rate (HLR), and hydraulic residence time (HRT). However, HRT (how fast/slow water moves through a wetland) is the key factor that determines treatment efficiency and the quality of water at the wetland Water 2020, 12,1937 2 of 20 outlet [9,10]. Water flow through a wetland is affected by the presence of aquatic vegetation (both emergent and submersed), and their associated properties such as vegetation type (rooted, floating, etc.), size, density, and stem flexibility. The mere presence of aquatic vegetation (e.g., emergent and submersed) may cause short circuiting, which is expected to occur particularly in large constructed wetlands (STAs). Short circuiting in a wetland leads to a short contact time with resident biota, and consequently reduces nutrient uptake and performance [1,9].
Despite the vast amount of research in literature regarding the interaction between prevailing winds and aquatic vegetation in general, little information is available to quantify EAV effects on wind speed and directions in large constructed wetlands (i.e., STAs and FEBs), and very little attention has been given to quantifying the interaction between prevailing wind and EAV by calculating roughness length or drag coefficient. Research involving wind effects on wetlands includes both laboratory and field observations. Laboratory experiments showed that the presence of emergent macrophytes significantly enhanced vertical transport, particularly gas transfer across the air-water interface [11]. The gas transfer coefficient also increased as wind speed increased.
In south Florida, previous field studies focused primarily on wetland water flow under hurricane conditions [12]. Sediment resuspension (measured by sediment traps) increased significantly in open water areas compared to areas protected from wind by emergent vegetation [13]. In addition, numerous studies measured drag coefficients of solid objects (e.g., cones, hemisphere, cubes, and pyramids) placed on flat surfaces in wind tunnels and flumes, while others used direct field measurements to calculate wind drag on trees for modeling purposes [14,15]. Furthermore, [16] reported that seasonal variations, up to 5-fold increases in surface roughness lengths, is due to the difference of growing seasons for natural surfaces including wetlands.
Most of the research involving wind effects on wetlands focused largely on high and low frequency of meteorological stressors that impact wetlands in different ways. Severe impacts due to low frequency and high magnitude events such as hurricane and big storms, may last for decades and cover exceptionally large areas [17]. Most high frequency and low magnitude events, such as rainfall, evapotranspiration (ET), inflow discharges and winds, exert a regular, almost predictable, influence on wetlands' ecosystem in terms of primary production, nutrient uptakes, and cycling. However, extraordinarily little is known about the impact of high frequency and low magnitude winds on wetlands in terms of vegetation resistance and wetland hydrology. Therefore, the major goals of this experiment are to answer the following questions, (1) What are the impacts, if any, of winds and EAV on wetland hydrology? (2) Determine if other weather parameters/stressors are being impacted due to the presence of EAV, (3) and if so, which ones are being impacted the most, and what are the implication of such change? Project objectives also included estimating drag coefficients to quantify the direct interaction of wind and EAV in wetlands.
Weather parameters were measured at three locations during wet and dry season in south Florida. Two weather stations were deployed in adjacent vegetated (EAV) and non-vegetated wetlands during wet and dry seasons ( Figure 1). The third location is a permanent station (ENR 308) located at the standard ten-meter height (26 • 37 22.4" N and 80 • 26 21.38" W). Results obtained from the combined dataset of all three stations were used to assess winds and EAV interactions on wetlands hydrology and estimate cattail drag coefficients. Observed wind and water velocity measured in a large constructed wetland during a discharge event, were combined to illustrate the impact of prevailing winds (high frequency and low magnitude) on wetland hydrology.

Data Collection
Metrological data were collected using an Airmar 150WX Ultrasonic weather station (www. fondriest.com, Fairborn, OH, USA). The Airmar 150WX is a compact weather station ideal for moving platforms and wetland environments. The sensor outputs apparent and true wind speed and direction, barometric pressure, air temperature, and GPS location. The Airmar weather station rated for use on buoy platforms (IPX6 waterproof), and measures wind speed and direction with integrated 10 Hz GPS and 2-axis compass, tilt, temperature, pressure, and humidity sensor modules. Sampling intervals for all parameters were set at 15 min.
The initial monitoring period lasted for almost one month (wet season), between 12 June through 10 July 2014, while the second sampling interval covered approximately five months (dry season) between 2 November 2014 through 12 March 2015. The temporary weather stations were located at the center of each Test cell (Figure 1), less than 25 m apart, and were placed at the same exact height above water level (~2 m). The nearby permanent weather station (ENR 308) is set at the standard ten-meter height, approximately eight meters higher than the Airmar weather stations located at the test cells.

Wind Stress
For practical applications and when observed wind speed is less than 10 m s −1 , it is reasonable to use an empirical formula [18] to calculate wind stress as follow: where U 10 = wind speed at 10 m above the water surface; ρ air = 1.22 kg m −3 , air density; and C D = 0.0013 dimensionless drag coefficient; a typical value which gives τ in units of N m −2 , or Pascals (Pa).

Power Law Wind Profile
The wind profile power law describes the relationship between wind speeds at two different heights [19,20]. Most application of the wind profile power law are conducted to assess wind power at specific heights for wind turbines [21], and applications of atmospheric pollution dispersion models [22]. Because the wind profile at atmospheric boundary layer (0-2000 m) is generally logarithmic in nature, it is best approximated by using the log wind profile equation, which accounts for surface roughness and atmospheric stability. However, when surface roughness or atmospheric stability information is lacking, the wind profile power law is commonly used: where u = wind speed (ms −1 ) at height z (meters) above ground level and u r is the wind speed measured at a reference height z r (meters); reference height in this experiment is ten meters (ENR 308 station); and α is an empirically derived coefficient and is approximately = 0.143 for neutral stability in atmospheric conditions. The exponent value is almost constant, particularly if the difference between the two wind speed heights is less than 50 m [21].

Log Wind Profile
The simplest form to represent wind drag is to use a constant value like the one used in Equation (1). However, a more accurate and reliable approach would be to use a logarithmic wind profile above the surface [23]. For this field experiment, the neutral atmospheric stability conditions apply: where U z is the mean wind speed at height z (meters) above the ground surface, u * is the friction or shear velocity (ms −1 ), k is the Von Karman constant (approximately = 0.41), and d is the height in meters above the ground surface that equals to the depth of an air layer trapped in vegetation at which wind speed is zero as a result of flow obstacles. For cattail in this case, it is approximated at 2/3 of the average height (2.0 m) of the obstacles (2/3 × 2.0 = 1.33 m; [23]); z o is the roughness length, which is a corrective measure to account for the effect of the roughness of a surface or an obstacle (EAV in this case) on wind flow. The wind speed at the standard ten-meter height, measured at the ENR 308 station, is determined from the wind speed at 10 m: Water 2020, 12, 1937 5 of 20 While wind speed at the top of vegetation (lower height compared to standard 10-m height) is calculated as follows: Dividing Equation (3) by (4) results in a ratio of wind speeds: Which can be rearranged to derive roughness length Z 0 (m): We can then calculate wind drag coefficient using wind speed at the standard ten-meter heights: Drag coefficients were calculated using Equations (2) and (8) along with wind speed measurements from the ENR 308 (ten-meter height), vegetated wetland (cattail height = 2.0 m), and non-vegetated wetland/control (30 cm water depth).

Wind Setup
In general, surface water moves in the direction of wind causing a rise of water level at the downwind side of an inland water body. However, no previous studies have calculated wind setup in a wetland, presumably due to lack of practical methods, and most likely is the fact that wind friction is commonly considered to be part of variability of overall marsh flow frictions [24]. However, [25] described several methods to calculate wind setup for lakes and coastal areas. A simplified formula that gives a reasonably good estimate is given by [26]: where S = wind setup (m), U = wind speed (m/s), F u = wind fetch (km), and D = average depth over the wind fetch (m).
To assess the impact of wind set up at the wetland outflow, an overland flow equation following a general power law relationship can be used [24,27]: where Q = outflow (m 3 /s), W = mean flow width (m), D = average depth (m), a = empirical parameter, and b = empirical exponent parameter. Walker and Kadlec [27] suggest using 3 to 4 for typical STA or wetland cell. The relative increase of outflow discharge due to wind setup can be estimated based on Equation (11) as: where Q s = outflow with wind setup (m 3 /s), and S = wind setup (m).

Wind Stress
At the experiment site, water flows in a northerly direction (i.e., from south to north) between inflow and outflow structures. Observed wind vectors (speed and directions) were decomposed into two components, one aligned north-south and the other east-west direction to match water-flow alignment between inflow and outflow structures. Wind stress calculated using Equation (1), was also decomposed along these two directions, which resulted in a north-south and east-west stress component, pushing water along those directions (Table 1).   1 -X-axis is oriented north-south direction (north is positive); 2 -Y-axis is oriented east-west direction (east is positive); 3 -units are in Pa; 4 -units m s −1 ; and 5 -statistically significant difference at p < 0.001.
To illustrate wind impacts on wetland hydrology and the generation of wind ripples in open water areas in constructed wetlands, we measured surface water velocities in Cell 3 of STA-2 where submerged aquatic vegetation (SAV) dominated, using a portable current meter (YSI, Handheld ADV FlowTracker©, San Diego, CA, USA). Velocity measurements were taken during a flow discharge event of 14.16 m 3 /s, which is the average discharge observed during previous field experiments.

Data Analysis
Descriptive statistics with graphical techniques were used to identify seasonality or trends, and describe relationships, if any, in time series observations. All analyses were performed with the raw data (i.e., every 15 min), and hourly means calculated from the original observations. Statistics were calculated with JMP (Version12, SAS Institute Inc., Cary, NC, USA) and SigmaPlot (Version 13, SPSS, Inc. Chicago, IL, USA). The level of significance (α) was set at 0.05 or lower as indicated for all analyses. Wind speed probability distributions were calculated to reveal the temporal characteristics of the data sets [28,29]. WRPLOT View software was used to determine wind vector and wind classes from field observations and display wind rose plots (http://www.weblakes.com/company/index.html, Waterloo, Canada) that depict the frequency of occurrence of winds in each of the specified wind-direction sectors and wind-speed classes for a given location and time period.

Wind Speed and Direction
The overall wind speed during the wet season ranged between 0. 10  Observed wind speed and direction in vegetated and non-vegetated sites were significantly different for both dry and wet seasons (Table 2 and Figure 2). Wind speed means of the vegetated wetland during both seasons were half their counterparts in the non-vegetated site ( Table 2). Wind direction in vegetated wetland was significantly different from the open water one ( Figure 2). Open water wetland's 25th and 75th percentiles of wind speed covered a much wider range compared to the vegetated one during both wet and dry season ( Figure 2). However, 25th and 75th percentile of wind direction covered the same range, more or less, during both seasons. The difference in both mean and median values between the two groups is statistically significant ( Table 2).  Wind speed exceedance probability analysis, based on the 15-min observations, indicated that zero occurrence of wind speed above 13 m s −1 during the observation period, at both Cell 12 (EAV) and Cell 13 (open water) locations ( Figure 3). In general, wind speed was higher in the open water compared to the vegetated site and two to three times higher at the standard height (ENR 308) compared to the vegetated site. Median wind speed was 1.  Table 2).
Water 2020, 12, 1937 8 of 19 Wind speed exceedance probability analysis, based on the 15-min observations, indicated that zero occurrence of wind speed above 13 m s −1 during the observation period, at both Cell 12 (EAV) and Cell 13 (open water) locations ( Figure 3). In general, wind speed was higher in the open water compared to the vegetated site and two to three times higher at the standard height (ENR 308) compared to the vegetated site. Median wind speed was 1.  Wind speed also varied during "time-of-day" in both seasons (Figures 4 and 5). Wind speed was higher during daytime (i.e., between 800 through 1900), compared to night-time hours, and decreased during night-and early-day hours (i.e., between 1900 and 700) at both sites and during both seasons. The time-of-day trend was more pronounced during the wet, compared to the dry season ( Figure 4 vs. 5). The maximum speed for an entire day occurred between 1300 and 1800 depending on the season (Figures 4 and 5). Hour-of-day box and whisker plots also demonstrate how the low wind speed values varied in the vegetated, compared to the non-vegetated site. Wind speed values were at their lowest from 0 to 900 h and increased to their highest values starting at hour 900 and peaked around hour 1400-1600 every day (Figures 4 and 5).
The overall average wind direction (resultant vector) for the wet sampling seasons (June-July) were 273° and 205°, for vegetated and open water sites, respectively ( Figure 6). Similarly, during the dry season, wind vectors were 14° and 20°, for vegetated and open water sites, respectively. During the dry season, wind direction was dominated by northerly wind, and westerly and south-westerly direction during the wet season ( Figure 6). Wind speed also varied during "time-of-day" in both seasons (Figures 4 and 5). Wind speed was higher during daytime (i.e., between 800 through 1900), compared to night-time hours, and decreased during night-and early-day hours (i.e., between 1900 and 700) at both sites and during both seasons. The time-of-day trend was more pronounced during the wet, compared to the dry season ( Figure 4 vs. Figure 5). The maximum speed for an entire day occurred between 1300 and 1800 depending on the season (Figures 4 and 5). Hour-of-day box and whisker plots also demonstrate how the low wind speed values varied in the vegetated, compared to the non-vegetated site. Wind speed values were at their lowest from 0 to 900 h and increased to their highest values starting at hour 900 and peaked around hour 1400-1600 every day (Figures 4 and 5).
The overall average wind direction (resultant vector) for the wet sampling seasons (June-July) were 273 • and 205 • , for vegetated and open water sites, respectively ( Figure 6). Similarly, during the dry season, wind vectors were 14 • and 20 • , for vegetated and open water sites, respectively. During the dry season, wind direction was dominated by northerly wind, and westerly and south-westerly direction during the wet season ( Figure 6).
Water 2020, 12,1937; doi:10.3390/w12071937 www.mdpi.com/journal/water Air temperature, relative humidity, and dew point values were un-affected by the presence of emergent vegetation in both wet and dry season, and the observed differences were not statistically significant ( Table 2). Although the differences are not significant, they are consistently higher in the vegetated Cell (12), particularly during the wet season, compared to the dry season (Table 2). At the study site barometric pressure range varied between a minimum of 20 mb to a maximum of 33 mb at the vegetated and un-vegetated sites (Table 2). While, the range in barometric pressure at the ENR 308 station (about 2200 m away and standard ten-meter height) was approximately 35 mb over a fiveyear period.

Wind Stress and Wind Drag Coefficients
The differences in the means and the medians values of the wind-generated surface sheer stress, during wet season for Cell 12 and 13, respectively, are statistically significant (p ≤ 0.001) along both N-S and E-W directions (Table 1). Similarly, wind stress means and medians during the dry season are statistically significant (p ≤ 0.001) and their values along the N-S and E-W directions are an order of magnitude different for dry compared to wet season values (Table 1).
Calculated wind speed profiles (wind speed vs. elevation), using Equations (2) and (8), are different between dry and wet seasons. The range of wind speed is narrower during wet season than during dry season, using the standard height of wind speed measured at ENR 308, Cell 12, and Cell 13 ( Figure 7). For both seasons, the increase in wind speed with height is more gradual at the vegetated site (Cell 12) compared to that at the open water site (Cell 13).

Air Temperature, Relative Humidity, Dew Point, and Barometric Pressure
Air temperature, relative humidity, and dew point values were un-affected by the presence of emergent vegetation in both wet and dry season, and the observed differences were not statistically significant ( Table 2). Although the differences are not significant, they are consistently higher in the vegetated Cell (12), particularly during the wet season, compared to the dry season (Table 2). At the study site barometric pressure range varied between a minimum of 20 mb to a maximum of 33 mb at the vegetated and un-vegetated sites (Table 2). While, the range in barometric pressure at the ENR 308 station (about 2200 m away and standard ten-meter height) was approximately 35 mb over a five-year period.
The standard barometric pressure of 1013.25 mb (millibar) and the less than 2 mb per three hours (mb/3hr) rate of change, are both within expected ranges associated with diurnal oscillation in temperate climates. The barometric pressure in the tropics simply does not change a lot from day to day and barometric pressure (corrected to sea level) varies from 1008.7 to 1018.9 mb. Barometric pressure at Miami, Florida, for example, might change from 1002.0 mbar to 1032.4 mb, a range about 30 mb. In higher latitudes, variation from 981.65 mb to 1049.4 mb, a 68 mb range, is not uncommon. The observed ranges are within the expected values for this region. We also found no major differences in results obtained for the aforementioned parameters (except wind speed and directions), due to presence or absence of emergent vegetation at the two wetland cells. Therefore, no further discussion regarding those parameters will continue due to lack of significant changes observed between data collected in the vegetated and the open water sites.

Wind Speed and Direction
Wind speed and direction means in the vegetated wetland were half their counterparts in the non-vegetated site during wet and dry seasons ( Table 2). The difference in both mean and median values between the vegetated and non-vegetated sites were statistically significantly different for both dry and wet seasons (Table 2 and Figure 2). Similar to vegetation resistance to water movement in the water column induced by the presence of aquatic vegetation [24,30], emergent vegetation in this experiment, induced resistance to air flow (which can be quantified by a drag coefficient) and consequently impacts the frequency of wind classes compared to open water site. The observed reduction in wind speed, is attributed to the mere presence of an obstacle, emergent vegetation in this experiment.
The highest frequencies of wind classes (0 ≤ wind speed ≤ 2.1 m s −1 ) were 92% and 64%, while the second highest (2.1 ≤ wind speed ≤ 4.1 m s −1 ) were 29% and 5%, for vegetated and open water site, respectively, during the wet season sampling period ( Figure 6). The main reason for the observed high frequency (92%) at the lowest wind class in the vegetated site is attributed to the presence of emergent vegetation. Vegetation resistance reduced wind speed and led to the observed highest frequency (92%) to occur at the lowest wind speed class (0 ≤ wind speed ≤ 2.1 m s −1 ). Similar to the wet season trend, the highest frequency of wind speed classes during the dry season was in the range of 0 to 2.1 m s −1 with values of 99% and 75% for vegetated and open water site, respectively. However, the second highest frequency (2.1 ≤ wind speed ≤ 4.1 m s −1 ) was less than 1% for the vegetated site, and 22% for the open water site during the dry season.
Wind speed values were remarkably similar, for both vegetated and open-water sites, at the higher speed (i.e., wind speed > 6 ms −1 ) compared to the low range, which may have greater implications. For example, during hurricanes, one would expect that both locations would be in close agreement due to higher wind speeds that damage existing emergent vegetation (larger impact on existing vegetation). On the other hand, during the current experiment, lower wind speed causes less "bending", prompting higher vegetation resistance, and causing the highest wind speed frequency to occur at the low wind as evident from our results (highest wind speed class occurring at 92 and 99% for wet and dry seasons, respectively). Similar to observed results in this experiment, [31] showed deceleration of flow within the vegetation layer and a maximum Reynolds stress at the tops of the plants were similar to the observed vertical profiles for velocity and Reynolds stress. The higher wind speeds would have greater impacts on emergent vegetation to the point that measured wind speed are equal at top of vegetation and water surface.
Wind directions may have also been impacted by the emergent vegetation bending as a result of the prevailing wind at the study site. Wind direction vectors over the entire wet sampling season (June-July) were 273 • and 205 • , for vegetated and open water sites, respectively, while wind vectors were 14 • and 20 • , for vegetated and open water sites, respectively during the dry season, (Figure 6). The dominant wind direction at the study site is from the north, as indicated by the open water site, during the dry season, and this dominant wind direction change to south-westerly and westerly during the wet season. However, it is clear that the dominant wind direction is from the north in the open water site, while wind direction is a "scatter" around the true north during the dry season in the vegetated site. The observed scatter in wind directions was triggered by the emergent vegetation bending caused by the prevailing winds. We believe that the prevailing wind sways (or bends) the emergent vegetation and creates channel-like patterns in the vegetated site, which not only reduce wind speed but also alter wind direction.

Wind Stress and Drag Coefficients
Wind also induces stress at the water surface, causing the surface layer of the water column to move, in general, in the same direction as wind, due to the small fetch and negligible effects from Coriolis force [32]. Wind stress over surface water produces currents in the top layer aligned with the wind direction, which in turn lead to wind setup and water surface slope. Although wind stress from normal wind conditions may be exceedingly small, the small wind setup may still have an impact on HRT particularly for a shallow wetland.
To demonstrate the effect of wind setup on wetland hydrology, Equations (10) and (11) were used to estimate wind impact set up near the outlet of the open water cell (Figure 7). If the maximum wind speed of 19 m s −1 , observed during this study, were in the same direction as the flow, the wind setup could be as high as 0.6 cm over a depth of 80 cm at the wetland outlet and would increase outflow discharges by 3%. The effect over shallower water is more pronounced. For a water depth of 20 cm, the wind setup is 2.5 cm, and the flow increase by 61%. The same wind condition (i.e., wind speed of 19 m s −1 ) would lead to a wind setup of 1.3 cm over a water depth of 40 cm and increase outflow discharges by 13% (Figure 7).
In addition to the effect of wind set up and increasing water depth at the wetland outlet, wind drag may also increase surface water velocity, and if aligned in the same direction as the surface water flow, decrease HRT, and consequently reduces phosphorus treatment. Several researchers have concluded that decreasing HRT would lead to an overall decrease in TP uptake, while increasing contact time (HRT) with resident biota, leads to improve P performance in a wetland [9,33].
Constructed wetlands by design are deeper near the outlet and shallower near the inlet, which allow gravity to move water between inflow and outflow. Wind set up would increase water depth near the wetland outlet and may cause unnecessary consequences because water depth is a critical factor in determining the health of aquatic vegetation in wetlands. Aquatic vegetation is a major sink in reducing TP in large constructed wetlands, particularly in south Florida due to an extended vegetation growing season. The health of emergent macrophytes, such as cattail, in those large constructed wetlands depends on how deep and how long those water depths persist. Cattail floated island were observed in the Everglades Removal Project (ENRP) when water depth increased to more than four feet and lasted for extended period (>two weeks; [33]). In summary, wind drag increases water velocity and decreases HRT, while wind set up increases water depth near the wetland outlet, and the consequences of both will depend on how long and how fast the wind is blowing.
Wind drag coefficient depends on several factors including vegetation flexibility, leaf density (vegetation porosity or number of stems/plants per m 2 ), prevailing wind speed, and vegetation structure (stem size, height, frontal area, etc.). Drag coefficient is also a function of wind velocity and the vegetation angle of incidence and increases as wind speed increases and also increases as leaf width increases [34]. However, wind drag, in most wetland applications, is considered to be part of the variability of marsh friction and is several times as strong as submerged plant stem drags in open water areas [24].
Vegetation impact on prevailing wind velocity is illustrated and quantified using Equations (2) and (8) (Figure 7). During wet and dry seasons, the increase in wind speed with height is more gradual at the vegetated site (Cell 12) compared to that at the open water site (Cell 13), indicating wind attenuating effect of vegetated site due to higher surface roughness length (Figure 7). Narrow wind speed range vs. height is evident during both wet and dry seasons ( Wind drag coefficients results (  (2) or (8) are almost in perfect agreement with values calculated for cropland, brush, and forest. Open water values from the current experiment are higher than the reported literature values for "smooth open water." However, it was clear during this experiment that the water surface at the open water site was dominated by wind waves and ripples (personal observations), which may not be classified or described as smooth open water. Roughness lengths in wetlands were reported to have a higher value in summer months compared to winter months due to differences in growing seasons [37]. The pattern of growing seasons is different in south Florida. Our data indicates that the roughness length and drag coefficients are lower in the wet season than in the dry season, probably due to lower water levels, which lead to less smooth overall surface.

Winds Impact on Wetland Hydrology
Velocity measurements were taken during a flow discharge event of 14.16 m 3 /s, which is the average discharge observed during previous field experiments. When wind blows across a water body, it pushes the surface water layer via shear stress. This downward transfer of momentum from the air to the water generates a deformed surface current (ripples/wind wave). Those generated waves become bigger as wind speed increases.
The water level responses to wind stress depend on the site location and the wind direction. In this experiment, wind direction is unlikely to change significantly across the full-size wetland (i.e., STA-2 Cell 3; STA2C3), due to the absence of cattail where velocity measurements were taken on 1 March 2014. At station near the inflow, and over the entire STA2C3, wind was blowing from the south, while inflow discharges were moving from north to south (Figure 8). Wind speed was measured at 0.53 m s −1 . The measured velocity profile at that station showed that surface wind-generated water flow, is aligned with wind direction from south to north, and opposite to inflow discharges (north to south; Figure 8). . Surface layer exhibited water flow from north to south (caused by wind stress) coinciding with wind direction from south to north; measured wind speed is approximately (0.53 m s −1 ). A mid-water column layer flowing between inflow and outflow structures from north to south (opposite to wind direction). A near-bottom layer appears to be flowing from north to south, likely due to backscatter from acoustic measurements.
The opposing wind-generated surface flow in the opposite direction of gravity flow may significantly impact HRT in a constructed wetland. For example, during a 14.16 m 3 s −1 (34 cm day −1 or 1.1 ft day −1 ) flow discharge event in STA2C3, wind-generated flow was in the opposite direction to gravity flow. The wind-generated surface water flow (and opposite to inflow discharge direction), during this discharge event, is approximately 0.15-m-deep (Figure 8), and a width of approximately 2000 m. This indicates that up to 12% of the water volume, during this discharge event, could be moving in the opposite direction of gravity discharges (Figure 8). When discharges at the inflow site stop (G333), the wind does not. Instead, wind continues to generate flow in the northerly direction creating a 15-cm surface layer moving back towards the inflow site, while the sub-surface layer continues in a southerly direction due to its momentum, and eventually stops. Those opposing currents generate more velocity shear (−7 to 25 cm s −1 ; Figure 8), which in turn creates more mixing/entrainment and re-cycling of water between the two layers, causing an increase in HRT and improvement of P performance. However, if wind and gravity discharges are aligned in the same direction, wind generated flow would decrease HRT and potentially reduce P performance.
Increasing water velocity between inflow and outflow is likely to decrease HRT, a key factor that impacts and determines P performance in a wetland [1][2][3]. For example, when prevailing wind is opposite to flow direction, wind generated current increases HRT (i.e., contact time) and is likely to increase P performance. In this experiment water moved between inflow and outflow in S-N direction. The calculated wind stress along the x-axis (S-N direction) in the non-vegetated wetland (during the dry seasons) is positive and aligned with the flow discharges in this cell (i.e., S-N). Mean wind stress in vegetated wetland (even smaller compared to non-vegetated wind stress) is also positive and aligned with flow-discharges in the S-N direction for both seasons. However, for the open water during the dry season, wind stress components are negative (i.e., N-S) along the x-axis and opposite to inflow discharges from S-N. Alignment of wind stress with water flow is likely to decrease HRT, which decreases contact time, and consequently P performance in wetlands.
Phosphorus performance in wetlands is the result of too many factors, yet previous researchers pointed out that, in general increasing contact time (HRT) leads to improve P performance. For example, Ruel et al., [33] changed bottom topography by creating baffled wetlands with multiple vertical-scale topographic features to increase HRT, while Dierberg et al., [9] used different HRT in a comparative mesocosm study involving SAV. In both cases increasing HRT leads to improved P performance in those constructed wetlands. Many researchers also identified hydraulic factors that impact wetland treatment efficiency [8]. However, to our knowledge, none reported the effect of a single factor impact such as increasing or decreasing HRT, due to wind action, on nutrient transport and efficiency in large constructed wetland. Most studies focused on the collective effects of several factors combined with HRT (e.g., [3,9]), and literature reported discussions indicated that HRT is expected to impact and change P removal efficiency in a wetland, yet this relationship is complex.
Because no monitoring of P between inflow and outflow took place during this experiment, we were not able to determine how wind-generated flow may have impacted P removal efficiency. In addition, even if P monitoring took place at the inflow and outflow sites, the reduction in P concentrations at the wetland outlet, is the result of the combined effects of too many factors, not individually and solely, by HRT alone. In this case, P performance would be the result of the combined interaction of soil, vegetation, HRT, and wind at a minimum.
In most cases when designing and constructing a wetland for P treatment, little attention has been paid for wind action and its impact on HRT and P performance. Our study showed that wind effect would decrease HRT if wind and flow are in the same directions ( Figure 8). In addition, wind speed and direction are natural phenomena that cannot be controlled. However, one could potentially use a combination of vegetation types in a wetland to change wind speed/direction and optimize HRT and hence P performance. For example, a combination of EAV and SAV strips would be a good choice to achieve this goal. The EAV and SAV in the same wetland would provide high and low vegetation resistance and wind drag coefficients, respectively. However, allowing only strips of SAV within EAV would ensure a small fetch for the wind, resulting in lower speed, and hence, lessen the impact on HRT and P performance.

Conclusions
Meteorological data from vegetated and un-vegetated wetlands during wet and dry seasons, were collected and analyzed to evaluate the role of wind and vegetation on wetland hydrology. Wind speed reductions and changes in wind direction were more pronounced during the wet season compared to the dry season. Vegetation induced resistance reduced wind speed and caused the lowest wind speed class to be the highest frequency class. Water flow, from inflow to outflow, is accelerated when wind direction is aligned with gravity flow and decreased when wind direction (wind-generated flow as a result of wind stress) opposes gravity flow from the wetland inlet.
We demonstrated how wind forces affect surface water flow, by measuring velocity profiles and wind speed/direction in a large constructed wetland. The water velocity profile showed that surface water flow was aligned in the same direction as the wind and was exactly opposite to inflow discharges. Those opposing currents generate more velocity shear, which in turns created more mixing/entrainment and re-cycling of water within the wetland, causing an increase in HRT and improvement of P performance. However, if wind and gravity discharges are aligned in the same direction, wind generated flow would decrease HRT and potentially reduce P performance.
We demonstrated how wind induced stress at the water surface enhanced flow velocity and described how wind stress impact was more pronounced at shallower water depth. We also disclosed how a wind speed of 19 m s −1 resulted in a wind set up, near the wetland outflow, which increased water depth by more than 12% and increased surface flow by as much 61% in a shallow wetland (water depth ≤ 20 cm). However, the impact on a deeper water depth (40 cm) was less pronounced for the same wind speed.
We were also able to calculate roughness lengths and drag coefficients for both vegetated and un-vegetated sites based on field observations. Calculated drag coefficients changed between dry and wet season due to changes in vegetation porosity (space between individual plants and between stems or leaves), caused mainly by prevailing winds. Cattails and open water wind drag coefficients using the log wind profile, were 0.016 and 0.009 for dry season, and 0.012 and 0.005 for wet season, respectively. This study highlights the need to account for changing vegetation resistance (drag coefficients) over time in future model applications and to consider prevailing wind in the design phase of constructed wetlands.
Author Contributions: M.M. and N.W. collaborated fully and equally in all aspect of this research, including conceptualization, methodology, validation, formal analysis, investigation, writing-original draft preparation, writing-review and editing, and visualization. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.