skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Dispersion of Buoyant and Sinking Particles in a Simulated Wind‐ and Wave‐Driven Turbulent Coastal Ocean
Abstract In shallow coastal oceans, turbulent flows driven by surface winds and waves and constrained by a solid bottom disperse particles. This work examines the mechanisms driving horizontal and vertical dispersion of buoyant and sinking particles for times much greater than turbulent integral time scales. Turbulent fields are modeled using a wind‐stress driven large eddy simulation (LES), incorporating wave‐driven Langmuir turbulence, surface breaking wave turbulent kinetic energy inputs, and a solid bottom boundary. A Lagrangian stochastic model is paired to the LES to incorporate Lagrangian particle tracking. Within a subset of intermediate buoyant rise velocities, particles experience synergistic vertical mixing in which breaking waves (BW) inject particles into Langmuir downwelling velocities sufficient to drive deep mixing. Along‐wind dispersion is controlled by vertical shear in mean along‐wind velocities. Wind and bottom friction‐driven vertical shear enhances dispersion of buoyant and sinking particles, while energetic turbulent mixing, such as from BW, dampens shear dispersion. Strongly rising and sinking particles trapped at the ocean surface and bottom, respectively, experience no vertical shear, resulting in low rates of along‐wind dispersion. Crosswind dispersion is shaped by particle advection in wind‐aligned fields of counter‐rotating Langmuir and Couette roll cells. Langmuir cells enhance crosswind dispersion in neutrally to intermediately buoyant particles through enhanced cell hopping. Surface trapping restricts particles to Langmuir convergence regions, strongly inhibiting crosswind dispersion. In shallow coastal systems, particle dispersion depends heavily on particle buoyancy and wave‐dependent turbulent effects.  more » « less
Award ID(s):
1634578
PAR ID:
10449897
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Oceans
Volume:
126
Issue:
3
ISSN:
2169-9275
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract This study investigates the dynamics of velocity shear and Reynolds stress in the ocean surface boundary layer for idealized misaligned wind and wave fields using a large-eddy simulation (LES) model based on the Craik–Leibovich equations, which captures Langmuir turbulence (LT). To focus on the role of LT, the LES experiments omit the Coriolis force, which obscures a stress–current-relation analysis. Furthermore, a vertically uniform body force is imposed so that the volume-averaged Eulerian flow does not accelerate but is steady. All simulations are first spun-up without wind-wave misalignment to reach a fully developed stationary turbulent state. Then, a crosswind Stokes drift profile is abruptly imposed, which drives crosswind stresses and associated crosswind currents without generating volume-averaged crosswind currents. The flow evolves to a new stationary state, in which the crosswind Reynolds stress vanishes while the crosswind Eulerian shear and Stokes drift shear are still present, yielding a misalignment between Reynolds stress and Lagrangian shear (sum of Eulerian current and Stokes drift). A Reynolds stress budgets analysis reveals a balance between stress production and velocity–pressure gradient terms (VPG) that encloses crosswind Eulerian shear, demonstrating a complex relation between shear and stress. In addition, the misalignment between Reynolds stress and Eulerian shear generates a horizontal turbulent momentum flux (due to correlations of along-wind and crosswind turbulent velocities) that can be important in producing Reynolds stress (due to correlations of horizontal and vertical turbulent velocities). Thus, details of the Reynolds stress production by Eulerian and Stokes drift shear may be critical for driving upper-ocean currents and for accurate turbulence parameterizations in misaligned wind-wave conditions. 
    more » « less
  2. Abstract Dispersion processes in the ocean surface boundary layer (OSBL) determine marine material distributions such as those of plankton and pollutants. Sheared velocities drive shear dispersion, which is traditionally assumed to be due to mean horizontal currents that decrease from the surface. However, OSBL turbulence supports along-wind jets; located in near-surface convergence and downwelling regions, such turbulent jets contain strong local shear. Through wind-driven idealized and large-eddy simulation (LES) models of the OSBL, this study examines the role of turbulent along-wind jets in dispersing material. In the idealized model, turbulent jets are generated by prescribed cellular flow with surface convergence and associated downwelling regions. Numeric and analytic model solutions reveal that horizontal jets substantially contribute to along-wind dispersion for sufficiently strong cellular flows and exceed contributions due to vertical mean shear for buoyant surface-trapped material. However, surface convergence regions also accumulate surface-trapped material, reducing shear dispersion by jets. Turbulence resolving LES results of a coastal depth-limited ocean agree qualitatively with the idealized model and reveal long-lived coherent jet structures that are necessary for effective jet dispersion. These coastal results indicate substantial jet contributions to along-wind dispersion. However, jet dispersion is likely less effective in the open ocean because jets are shorter lived, less organized, and distorted due to spiraling Ekman currents. 
    more » « less
  3. This study utilizes a large-eddy simulation (LES) approach to systematically assess the directional variability of wave-driven Langmuir turbulence (LT) in the ocean surface boundary layer (OSBL) under tropical cyclones (TCs). The Stokes drift vector, which drives LT through the Craik–Leibovich vortex force, is obtained through spectral wave simulations. LT’s direction is identified by horizontally elongated turbulent structures and objectively determined from horizontal autocorrelations of vertical velocities. In spite of a TC’s complex forcing with great wind and wave misalignments, this study finds that LT is approximately aligned with the wind. This is because the Reynolds stress and the depth-averaged Lagrangian shear (Eulerian plus Stokes drift shear) that are key in determining the LT intensity (determined by normalized depth-averaged vertical velocity variances) and direction are also approximately aligned with the wind relatively close to the surface. A scaling analysis of the momentum budget suggests that the Reynolds stress is approximately constant over a near-surface layer with predominant production of turbulent kinetic energy by Stokes drift shear, which is confirmed from the LES results. In this layer, Stokes drift shear, which dominates the Lagrangian shear, is aligned with the wind because of relatively short, wind-driven waves. On the contrary, Stokes drift exhibits considerable amount of misalignments with the wind. This wind–wave misalignment reduces LT intensity, consistent with a simple turbulent kinetic energy model. Our analysis shows that both the Reynolds stress and LT are aligned with the wind for different reasons: the former is dictated by the momentum budget, while the latter is controlled by wind-forced waves. 
    more » « less
  4. Abstract Mixing processes in the upper ocean play a key role in transferring heat, momentum, and matter in the ocean. These mixing processes are significantly enhanced by wave‐driven Langmuir turbulence (LT). Based on a paired analysis of observations and simulations, this study investigates wind fetch and direction effects on LT at a coastal site south of the island Martha’s Vineyard (MA, USA). Our results demonstrate that LT is strongly influenced by wind fetch and direction in coastal oceans, both of which contribute to controlling turbulent coastal transport processes. For northerly offshore winds, land limits the wind fetch and wave development, whereas southerly winds are associated with practically infinite fetch. Observed and simulated two‐dimensional wave height spectra reveal persistent southerly swell and substantially more developed wind‐driven waves from the south. For oblique offshore winds, waves develop more strongly in the alongshore direction with less limited fetch, resulting in significant wind and wave misalignments. Observations of coherent near‐surface crosswind velocities indicate that LT is only present for sufficiently developed waves. The fetch‐limited northerly winds inhibit wave developments and the formation of LT. In addition to limited fetch, strong wind–wave misalignments prevent LT development. Although energetic and persistent, swell waves do not substantially influence LT activity during the observation period because these relatively long swell waves are associated with small Stokes drift shear. These observational results agree well with turbulence‐resolving large eddy simulations (LESs) based on the wave‐averaged Navier–Stokes equation, validating the LES approach to coastal LT in the complex wind and wave conditions. 
    more » « less
  5. Langmuir turbulence consists of Langmuir circulation (LC) generated at the surface of rivers, lakes, bays, and oceans by the interaction between the wind-driven shear and surface gravity waves. In homogeneous shallow water, LC can extend to the bottom of the water column and interact with the bottom boundary layer. Large-eddy simulation (LES) of LC in shallow water was performed with the finite volume method and various forms of subgrid-scale (SGS) model characterized by different near-wall treatments of the SGS eddy viscosity. The wave forcing relative to wind forcing in the LES was set following the field measurements of full-depth LC during the presence of LC engulfing a water column 15 m in depth in the coastal ocean, reported in the literature. It is found that the SGS model can greatly impact the structure of LC in the lower half of the water column. Results are evaluated in terms of (1) the Langmuir turbulence velocity statistics and (2) the lateral (crosswind) length scale and overall cell structure of LC. LES with an eddy viscosity with velocity scale in terms of S and Ω (where S is the norm of the strain rate tensor and Ω is the norm of the vorticity tensor) and a Van Driest wall damping function (referred to as the S-Omega model) is found to provide best agreement with pseudo-spectral LES in terms of the lateral length scale and overall cell structure of LC. Two other SGS models, namely the dynamic Smagorinsky model and the wall-adapting local-eddy viscosity model are found to provide less agreement with pseudo- spectral LES, for example, as they lead to less coherent bottom convergence of the cells and weaker associ ated upward transport of slow downwind moving fluid. Finally, LES with the S-Omega SGS model is also found to lead to good agreement with physical measurements of LC in the coastal ocean in terms of Langmuir turbulence decay during periods of surface heating 
    more » « less