skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, May 17 until 8:00 AM ET on Saturday, May 18 due to maintenance. We apologize for the inconvenience.

Title: Wave‐Driven Sediment Resuspension Within a Model Eelgrass Meadow

Wave velocity and suspended sediment concentration were measured over a sand bed with and without a model eelgrass meadow. The model meadow was geometrically and dynamically similar to the marine eelgrassZostera marina. Meadows were constructed with three stem densities: 280, 600, and 820 stems/m2. Ripples formed within the meadow only when the spacing between stem rows was larger than the wave excursion. When ripples formed, the ripple geometry was the same as that observed for bare bed. When ripples were present, the near‐bed turbulent kinetic energy (TKE) was dominated by the ripple‐generated turbulence, and both the near‐bedTKEand averaged suspended sediment concentration were similar across all meadow densities and bare bed at the same wave velocity. When ripples were absent, the near‐bedTKEwas dominated by the stem‐generated turbulence, and the averaged suspended sediment concentration was reduced, compared to cases with ripples but at the same wave velocity. For conditions with and without a model meadow, the sediment diffusivity inferred from vertical profiles of suspended sediment concentration increased linearly with distance from the bed.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Page Range / eLocation ID:
p. 1035-1053
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Laboratory experiments examined the impact of model vegetation on wave‐driven resuspension. Model canopies were constructed from cylinders with three diameters (d= 0.32, 0.64, and 1.26 cm) and 12 densities (cylinders/m2) up to a solid volume fraction (ϕ) of 10%. The sediment bed consisted of spherical grains withd50= 85 μm. For each experiment, the wave velocity was gradually adjusted by increasing the amplitude of 2‐s waves in a stepwise fashion. A Nortek Vectrino sampled the velocity atz= 1.3 cm above the bed. The critical wave orbital velocity for resuspension was inferred from records of suspended sediment concentration (measured with optical backscatter) as a function of wave velocity. The critical wave velocity decreased with increasing solid volume fraction. The reduction in critical wave velocity was linked to stem‐generated turbulence, which, for the same wave velocity, increased with increasing solid volume fraction. The measured turbulence was consistent with a wave‐modified version of a stem‐turbulence model. The measurements suggested that a critical value of turbulent kinetic energy was needed to initiate resuspension, and this was used to define the critical wave velocity as a function of solid volume fraction. The model predicted the measured critical wave velocity for stem diametersd= 0.64 to 2 cm. Combining the critical wave velocity with an existing model for wave damping defined the meadow size for which wave damping would be sufficient to suppress wave‐induced sediment suspension within the interior of the meadow.

    more » « less
  2. Abstract

    A new modeling methodology for ripple dynamics driven by oscillatory flows using a Eulerian two‐phase flow approach is presented in order to bridge the research gap between near‐bed sediment transport via ripple migration and suspended load transport dictated by ripple induced vortices. Reynolds‐averaged Eulerian two‐phase equations for fluid phase and sediment phase are solved in a two‐dimensional vertical domain with akεclosure for flow turbulence and particle stresses closures for short‐lived collision and enduring contact. The model can resolve full profiles of sediment transport without making conventional near‐bed load and suspended load assumptions. The model is validated with an oscillating tunnel experiment of orbital ripple driven by a Stokes second‐order (onshore velocity skewed) oscillatory flow with a good agreement in the flow velocity and sediment concentration. Although the suspended sediment concentration far from the ripple in the dilute region was underpredicted by the present model, the model predicts an onshore ripple migration rate that is in very good agreement with the measured value. Another orbital ripple case driven by symmetric sinusoidal oscillatory flow is also conducted to contrast the effect of velocity skewness. The model is able to capture a net offshore‐directed suspended load transport flux due to the asymmetric primary vortex consistent with laboratory observation. More importantly, the model can resolve the asymmetry of onshore‐directed near‐bed sediment flux associated with more intense boundary layer flow speed‐up during onshore flow cycle and sediment avalanching near the lee ripple flank which force the onshore ripple migration.

    more » « less
  3. Abstract

    A large‐scale laboratory experiment was conducted to evaluate cross‐shore sediment transport and bed response on a sandbar under erosive and accretive field‐scale wave conditions (total of 11 cases). Unprecedented vertical resolution of sediment concentration was achieved through the use of conductivity concentration profilers alongside miniature fiber optic backscatter profilers. Observations were made of intrawave (phase‐averaged) and wave‐averaged cross‐shore sediment flux profiles and transport rates in the lower half of the water column on the crest of a sandbar. The net sediment transport rate was partitioned into suspended sediment (SS) and bed load (BL) components to quantify the relative contributions of SS and BL to the total sediment transport rate. Net SS transport rates were greater than net BL transport rates for the positive (wave crest) half‐cycle in 6 of 11 cases, compared to 100% (11 of 11) for the negative (wave trough) half‐cycle. Net (wave‐averaged) BL transport rates were greater, in magnitude, than net SS transport rates for 7 of the 11 cases. The dominant mode of transport was determined from the ratio of net BL to net SS transport rate magnitudes. The net transport rate was negative (offshore‐directed) when SS dominated and positive (onshore‐directed) when BL dominated. Net BL transport rate correlated well with third moments of free‐stream velocity (r2 = 0.72), suggesting that energetics‐type quasi‐steady formulae may be suitable for predicting BL transport under the range of test conditions.

    more » « less
  4. Abstract

    We performed laboratory experiments to investigate the influence of sand content on the dynamics of wave‐supported gravity flows in mud‐dominant environments. The experiments were carried out in an oscillatory water tunnel with a sediment bed of either 1% or 13% sand. Low and high energy regimes are differentiated based on a Stokes Reynolds numberReΔ ≈ 500. In the low energy regime, the sand fraction influences flow dynamics primarily through ripple formation; no ripples form in the 1% sand experiments, whereas ripples form in the 13% experiments that increase turbulence and the wave boundary layer thickness,δm. In the high energy regime, small ripples form in both the 1% and 13% sand experiments and we observe high near‐bed suspended sediment concentrations. The influence of stratification on the boundary layer flow is characterized in terms of the gradient Richardson numberRig. The flow is weakly stratified inside the boundary layer for all runs and critically stratified at or above the top of the boundary layer. In the lower energy regime, the sand content reduces the relative influence of stratification in the boundary layer, shifting the elevation of critical stratification,LB, from approximately 1.3δmto 2.5δmin the 1% and 13% experiments, respectively. In both sets of experimentsLB ≈ δmat the strongest wave energy, indicating a transition to strongly stratified dynamics.

    more » « less
  5. Abstract

    Velocity and forces on individual plants were measured within an emergent canopy with real plant morphology and used to develop predictions for the vertical profiles of velocity and turbulent kinetic energy (TKE). Two common plant species,Typha latifoliaandRotala indica, with distinctive morphology, were considered.Typhahas leaves bundled at the base, andRotalahas leaves distributed over the length of the central stem. Compared to conditions with a bare bed and the same velocity, theTKEwithin both canopies was enhanced. For theTyphacanopy, for which the frontal area increased with distance from the bed, the velocity, integral length‐scale, andTKEall decreased with distance from the bed. For theRotala, which had a vertically uniform distribution of biomass, the velocity, integral length‐scale, andTKEwere also vertically uniform. A turbulence model previously developed for random arrays of rigid cylinders was modified to predict both the vertical distribution and the channel‐average ofTKEby defining the relationship between the integral length‐scale and plant morphology. The velocity profile can also be predicted from the plant morphology. Combining with the new turbulence model, theTKEprofile was predicted from the channel‐average velocity and plant frontal area.

    more » « less