Title: Langmuir Turbulence in the Arctic Ocean: Insights From a Coupled Sea Ice –Wave Model
Upper ocean mixing governs the vertical transport of heat, momentum, and tracers in the ocean surface boundary layer (OSBL), yet large-scale climate models often misrepresent its underlying processes, leading to significant uncertainty in sea ice and ocean predictions. Langmuir turbulence (LT) is one of the primary mechanisms of mixing in the open ocean and is generated by the interaction of wind stress and wave-induced Stokes drift. Observations have confirmed LT activity in leads, polynyas, and the marginal ice zone (MIZ), however its spatial and seasonal variability remains poorly constrained. In this study, we conduct the first Arctic-wide assessment of LT potential using a coupled sea ice–wave model that integrates neXtSIM and WAVEWATCH III. We analyze the spatiotemporal variability of LT by examining model-resolved turbulent dissipation and vertical kinetic energy within the OSBL. Our analysis reveals that LT potential is higher in the MIZ during melt and freeze-up, when partial sea ice cover allows intermittent wave propagation. Under these conditions, LT commonly coexists with wind-driven shear, forming a mixed-forcing regime that shapes upper-ocean energetics in response to evolving sea ice and wave states. Sea ice concentration and wind–wave alignment strongly influence the intensity and distribution of LT-driven mixing. On average, LT contributes roughly 15 % of the total upper-ocean dissipation in the Arctic MIZ, with episodic wave-driven events during transitional ice periods doubling local mixing rates compared to wind-only conditions. This analysis highlights the energetic role of wave-induced mixing in the upper ocean, with potential implications for vertical momentum transport, mixed layer structure, and sea ice–ocean interactions in the Arctic.  more » « less
Award ID(s):
2146889
PAR ID:
10682004
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
The Cryosphere
Date Published:
Format(s):
Medium: X
Institution:
Brown University
Sponsoring Org:
National Science Foundation
More Like this
  1. Measurements collected by a REMUS 600 AUV off the coast of southern California demonstrate large-scale coherent wave-driven vortices, consistent with Langmuir turbulence (LT), played a dominant role in structuring turbulent dissipation within the oceanic surface boundary layer. During a 10-hour period with sustained wind speeds of 10 m/s, 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 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 OSBL by redistributing wave-breaking turbulence. 
    more » « less
  2. Abstract Measurements collected by a REMUS 600 AUV off the coast of southern California demonstrate large-scale coherent wave-driven vortices, consistent with Langmuir turbulence (LT), played a dominant role in structuring turbulent dissipation within the oceanic surface boundary layer. During a 10-hour 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 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 OSBL by redistributing wave breaking turbulence. 
    more » « less
  3. The Arctic Ocean has experienced significant sea ice loss over recent decades, shifting towards a thinner and more mobile seasonal ice regime. However, the impacts of these transformations on the upper ocean dynamics of the biologically productive Pacific Arctic continental shelves remain underexplored. Here, we quantified the summer upper mixed layer depth and analyzed its interannual to decadal evolution with sea ice and atmospheric forcing, using hydrographic observations and model reanalysis from 1996 to 2021. Before 2006, a shoaling summer mixed layer was associated with sea ice loss and surface warming. After 2007, however, the upper mixed layer reversed to a generally deepening trend due to markedly lengthened open water duration, enhanced wind-induced mixing, and reduced ice meltwater input. Our findings reveal a shift in the primary drivers of upper ocean dynamics, with surface buoyancy flux dominant initially, followed by a shift to wind forcing despite continued sea ice decline. These changes in upper ocean structure and forcing mechanisms may have substantial implications for the marine ecosystem, potentially contributing to unusual fall phytoplankton blooms and intensified ocean acidification observed in the past decade 
    more » « less
  4. null (Ed.)
    Abstract Turbulence driven by wind and waves controls the transport of heat, momentum, and matter in the ocean surface boundary layer (OSBL). For realistic ocean conditions, winds and waves are often neither aligned nor constant, for example, when winds turn rapidly. Based on a Large Eddy Simulation (LES) method, which captures shear-driven turbulence (ST) and Langmuir turbulence (LT) driven by the Craik-Leibovich vortex force, we investigate the OSBL response to abruptly turning winds. We design idealized LES experiments, whose winds are initially constant to equilibrate OSBL turbulence before abruptly turning 90° either cyclonically or anticyclonically. The transient Stokes drift for LT is estimated from a spectral wave model. The OSBL response includes three successive stages that follow the change in direction. During stage 1, turbulent kinetic energy (TKE) decreases due to reduced TKE production. Stage 2 is characterized by TKE increasing with TKE shear production recovering and exceeding TKE dissipation. Transient TKE levels may exceed their stationary values due to inertial resonance and non-equilibrium turbulence. Turbulence relaxes to its equilibrium state at stage 3, but LT still adjusts due to slowly developing waves. During stages 1 and 2, greatly misaligned wind and waves lead to Eulerian TKE production exceeding Stokes TKE production. A Reynolds stress budget analysis and Reynolds-averaged Navier-Stokes equation models indicate that Stokes production furthermore drives the OSBL response. The Coriolis effects result in asymmetrical OSBL responses to wind turning directions. Our results suggest that transient wind conditions play a key role in understanding realistic OSBL dynamics. 
    more » « less
  5. Abstract This study investigates how Langmuir turbulence (LT) driven by Stokes drift shear affects the heated ocean surface boundary layer (OSBL) based on turbulence-resolving large-eddy simulations (LES) and assesses an analytic vertical mixing parameterization based on a simplified second-moment closure (SMC) approach. Diurnal solar heating forces OSBL shoaling to generate a diurnal warm layer (DWL) in which heat and momentum are trapped. Without LT, relatively weak turbulent mixing results in a near-surface jet that is associated with enhanced turbulent kinetic energy (TKE) production of shear-driven turbulence (ST), which approximately balances TKE dissipation rates. Conversely, LT maintains strong mixing, delaying the DWL formation and preventing the TKE dissipation enhancement by generating a less sheared jet. However, sufficiently strong heating destroys TKE to ultimately reduce mixing and to create more sheared jets, which effectively shifts the LT to an ST-dominated regime. A second-moment turbulence budget analysis suggests that 1) the near-surface OSBL responds rapidly to the surface forcing, 2) Stokes drift impacts heat and momentum budgets in profoundly different ways, and 3) buoyancy terms are to leading order negligible. Building on these findings and introducing a physics-based mixing length, we develop a simplified SMC model that can be solved for near-surface expressions for key turbulent variables and mixing coefficients in terms of known variables. For ST, these expressions are consistent with the Monin–Obukhov similarity theory. For LT, these expressions reveal a fundamental dependence of turbulent variables on Stokes drift shear. 
    more » « less