The turbulent ocean surface boundary layer (OSBL) shoals during daytime solar surface heating, developing a diurnal warm layer (DWL). The DWL significantly influences OSBL dynamics by trapping momentum and heat in a shallow near‐surface layer. Therefore, DWL depth is critical for understanding OSBL transport and ocean‐atmosphere coupling. A great challenge for determining DWL depth is considering wave‐driven Langmuir turbulence (LT), which increases vertical transport. This study investigates observations with moderate wind speeds (4–7 m/s at 10 m height) and swell waves for which breaking wave effects are less pronounced. By employing turbulence‐resolving large eddy simulation experiments that cover observed wind, wave, and heating conditions based on the wave‐averaged Craik‐Lebovich equation, we develop a DWL depth scaling unifying previous approaches. This scaling closely agrees with observed DWL depths from a year‐long mooring deployment in the subtropical North Atlantic, demonstrating the critical role of LT in determining DWL depth and OSBL dynamics.
Measurements collected by a Remote Environmental Monitoring Units (REMUS) 600 autonomous underwater vehicle (AUV) off the coast of southern California demonstrate large-scale coherent wave-driven vortices, consistent with Langmuir turbulence (LT), and played a dominant role in structuring turbulent dissipation within the oceanic surface boundary layer. During a 10-h period with sustained wind speeds of 10 m s−1, Langmuir circulations were limited to the upper third of the surface mixed layer by persistent stratification within the water column. The ensemble-averaged circulation, calculated using conditional averaging of acoustic Doppler dual current profile (AD2CP) velocity profiles using elevated backscattering intensity associated with subsurface bubble clouds, indicates that LT vortex pairs were characterized by an energetic downwelling zone flanked by broader, weaker upwelling regions with vertical velocity magnitudes similar to previous numerical studies of LT. Horizontally distributed microstructure estimates of turbulent kinetic energy dissipation rates were lognormally distributed near the surface in the wave mixing layer with the majority of values falling between wall layer scaling and wave transport layer scaling. Partitioning dissipation rates between downwelling centers and ambient conditions suggests that LT may play a dominant role in elevating dissipation rates in the ocean surface boundary layer (OSBL) by redistributing wave-breaking turbulence.
more » « less- Award ID(s):
- 1829952
- NSF-PAR ID:
- 10535527
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 54
- Issue:
- 9
- ISSN:
- 0022-3670
- Format(s):
- Medium: X Size: p. 1903-1920
- Size(s):
- p. 1903-1920
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Turbulent processes in the ocean surface boundary layer (OSBL) play a key role in weather and climate systems. This study explores a Lagrangian analysis of wave-driven OSBL turbulence, based on a large-eddy simulation (LES) model coupled to a Lagrangian stochastic model (LSM). Langmuir turbulence (LT) is captured by Craik–Leibovich wave forcing that generates LT through the Craik–Leibovich type 2 (CL2) mechanism. Breaking wave (BW) effects are modeled by a surface turbulent kinetic energy flux that is constrained by wind energy input to surface waves. Unresolved LES subgrid-scale (SGS) motions are simulated with the LSM to be energetically consistent with the SGS model of the LES. With LT, Lagrangian autocorrelations of velocities reveal three distinct turbulent time scales: an integral, a dispersive mixing, and a coherent structure time. Coherent structures due to LT result in relatively narrow peaks of Lagrangian frequency velocity spectra. With and without waves, the high-frequency spectral tail is consistent with expectations for the inertial subrange, but BWs substantially increase spectral levels at high frequencies. Consistently, over short times, particle-pair dispersion results agree with the Richardson–Obukhov law, and near-surface dispersion is significantly enhanced because of BWs. Over longer times, our dispersion results are consistent with Taylor dispersion. In this case, turbulent diffusivities are substantially larger with LT in the crosswind direction, but reduced in the along-wind direction because of enhanced turbulent transport by LT that reduces mean Eulerian shear. Our results indicate that the Lagrangian analysis framework is effective and physically intuitive to characterize OSBL turbulence.
-
more » « less
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
Abstract Data from an air–sea interaction tower are used to close the turbulent kinetic energy (TKE) budget in the wave-affected surface layer of the upper ocean. Under energetic wind forcing with active wave breaking, the dominant balance is between the dissipation rate of TKE and the downward convergence in vertical energy flux. The downward energy flux is driven by pressure work, and the TKE transport is upward, opposite to the downgradient assumption in most turbulence closure models. The sign and the relative magnitude of these energy fluxes are hypothesized to be driven by an interaction between the vertical velocity of Langmuir circulation (LC) and the kinetic energy and pressure of wave groups, which is the result of small-scale wave–current interaction. Consistent with previous modeling studies, the data suggest that the horizontal velocity anomaly associated with LC refracts wave energy away from downwelling regions and into upwelling regions, resulting in negative covariance between the vertical velocity of LC and the pressure anomaly associated with the wave groups. The asymmetry between downward pressure work and upward TKE flux is explained by the Bernoulli response of the sea surface, which results in groups of waves having a larger pressure anomaly than the corresponding kinetic energy anomaly, consistent with group-bound long-wave theory.
-
more » « less
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.
