Langmuir turbulence and filament frontogenesis in the oceanic surface boundary layer

2019 ◽  
Vol 879 ◽  
pp. 512-553 ◽  
Author(s):  
Peter P. Sullivan ◽  
James C. McWilliams

Submesoscale currents, small-scale turbulence and surface gravity waves co-exist in the upper ocean and interact in complex ways. To expose the couplings, the frontogenetic life cycle of an idealized cold dense submesoscale filament interacting with upper ocean Langmuir turbulence is investigated in large-eddy simulations (LESs) based on the incompressible wave-averaged equations. The simulations utilize large domains and fine meshes with $6.4\times 10^{9}$ grid points. Case studies are made with surface winds or surface cooling with waves oriented in across-filament (perpendicular) or down-filament (parallel) directions relative to the two-dimensional filament axis. The currents $u$, $v$ and $w$ are aligned with the across-filament, down-filament and vertical directions, respectively. Frontogenesis is induced by across-filament Lagrangian secondary circulations in the boundary layer, and it is shown to be strongly impacted by surface waves, in particular the propagation direction relative to the filament axis. In a horizontally heterogeneous boundary layer, surface waves induce both mean and fluctuating Stokes-drift vortex forces that modify a linear, hydrostatic turbulent thermal wind (TTW) approximation for momentum. Down-filament winds and waves are found to be especially impactful, they significantly reduce the peak level of frontogenesis by fragmenting the filament into primary and secondary down-welling sites in a broad frontal zone over a width ${\sim}500~\text{m}$. At peak frontogenesis, opposing down-filament jets $\langle v\rangle$ overlie each other resulting in a vigorous vertical shear layer $\unicode[STIX]{x2202}_{z}\langle v\rangle$ with large vertical momentum flux $\langle v^{\prime }w^{\prime }\rangle$. Filament arrest is induced by a lateral shear instability that generates horizontal momentum flux $\langle u^{\prime }v^{\prime }\rangle$ at low wavenumbers. The turbulent vertical velocity patterns, indicative of coherent Langmuir cells, change markedly across the horizontal domain with both across-filament and down-filament winds under the action of submesoscale currents.

2017 ◽  
Vol 837 ◽  
pp. 341-380 ◽  
Author(s):  
Peter P. Sullivan ◽  
James C. McWilliams

The evolution of upper ocean currents involves a set of complex, poorly understood interactions between submesoscale turbulence (e.g. density fronts and filaments and coherent vortices) and smaller-scale boundary-layer turbulence. Here we simulate the lifecycle of a cold (dense) filament undergoing frontogenesis in the presence of turbulence generated by surface stress and/or buoyancy loss. This phenomenon is examined in large-eddy simulations with resolved turbulent motions in large horizontal domains using${\sim}10^{10}$grid points. Steady winds are oriented in directions perpendicular or parallel to the filament axis. Due to turbulent vertical momentum mixing, cold filaments generate a potent two-celled secondary circulation in the cross-filament plane that is frontogenetic, sharpens the cross-filament buoyancy and horizontal velocity gradients and blocks Ekman buoyancy flux across the cold filament core towards the warm filament edge. Within less than a day, the frontogenesis is arrested at a small width,${\approx}100~\text{m}$, primarily by an enhancement of the turbulence through a small submesoscale, horizontal shear instability of the sharpened filament, followed by a subsequent slow decay of the filament by further turbulent mixing. The boundary-layer turbulence is inhomogeneous and non-stationary in relation to the evolving submesoscale currents and density stratification. The occurrence of frontogenesis and arrest are qualitatively similar with varying stress direction or with convective cooling, but the detailed evolution and flow structure differ among the cases. Thus submesoscale filament frontogenesis caused by boundary-layer turbulence, frontal arrest by frontal instability and frontal decay by forward energy cascade, and turbulent mixing are generic processes in the upper ocean.


2012 ◽  
Vol 39 (18) ◽  
Author(s):  
Stephen E. Belcher ◽  
Alan L. M. Grant ◽  
Kirsty E. Hanley ◽  
Baylor Fox-Kemper ◽  
Luke Van Roekel ◽  
...  

Nature ◽  
1985 ◽  
Vol 318 (6046) ◽  
pp. 519-522 ◽  
Author(s):  
S. A. Thorpe

2017 ◽  
Vol 47 (12) ◽  
pp. 2863-2886 ◽  
Author(s):  
Qing Li ◽  
Baylor Fox-Kemper

AbstractLarge-eddy simulations (LESs) with various constant wind, wave, and surface destabilizing surface buoyancy flux forcing are conducted, with a focus on assessing the impact of Langmuir turbulence on the entrainment buoyancy flux at the base of the ocean surface boundary layer. An estimate of the entrainment buoyancy flux scaling is made to best fit the LES results. The presence of Stokes drift forcing and the resulting Langmuir turbulence enhances the entrainment rate significantly under weak surface destabilizing buoyancy flux conditions, that is, weakly convective turbulence. In contrast, Langmuir turbulence effects are moderate when convective turbulence is dominant and appear to be additive rather than multiplicative to the convection-induced mixing. The parameterized unresolved velocity scale in the K-profile parameterization (KPP) is modified to adhere to the new scaling law of the entrainment buoyancy flux and account for the effects of Langmuir turbulence. This modification is targeted on common situations in a climate model where either Langmuir turbulence or convection is important and may overestimate the entrainment when both are weak. Nevertheless, the modified KPP is tested in a global climate model and generally outperforms those tested in previous studies. Improvements in the simulated mixed layer depth are found, especially in the Southern Ocean in austral summer.


2017 ◽  
Vol 47 (10) ◽  
pp. 2429-2454 ◽  
Author(s):  
Eric D. Skyllingstad ◽  
Jenessa Duncombe ◽  
Roger M. Samelson

AbstractGeneration of ocean surface boundary layer turbulence and coherent roll structures is examined in the context of wind-driven and geostrophic shear associated with horizontal density gradients using a large-eddy simulation model. Numerical experiments over a range of surface wind forcing and horizontal density gradient strengths, combined with linear stability analysis, indicate that the dominant instability mechanism supporting coherent roll development in these simulations is a mixed instability combining shear instability of the ageostrophic, wind-driven flow with symmetric instability of the frontal geostrophic shear. Disruption of geostrophic balance by vertical mixing induces an inertially rotating ageostrophic current, not forced directly by the wind, that initially strengthens the stratification, damps the instabilities, and reduces vertical mixing, but instability and mixing return when the inertial buoyancy advection reverses. The resulting rolls and instabilities are not aligned with the frontal zone, with an oblique orientation controlled by the Ekman-like instability. Mean turbulence is enhanced when the winds are destabilizing relative to the frontal orientation, but mean Ekman buoyancy advection is found to be relatively unimportant in these simulations. Instead, the mean turbulent kinetic energy balance is dominated by mechanical shear production that is enhanced when the wind-driven shear augments the geostrophic shear, while the resulting vertical mixing nearly eliminates any effective surface buoyancy flux from near-surface, cold-to-warm, Ekman buoyancy advection.


2016 ◽  
Vol 46 (2) ◽  
pp. 371-394 ◽  
Author(s):  
Jacob O. Wenegrat ◽  
Michael J. McPhaden

AbstractOcean currents in the surface boundary layer are sensitive to a variety of parameters not included in classic Ekman theory, including the vertical structure of eddy viscosity, finite boundary layer depth, baroclinic pressure gradients, and surface waves. These parameters can modify the horizontal and vertical flow in the near-surface ocean, making them of first-order significance to a wide range of phenomena of broad practical and scientific import. In this work, an approximate Green’s function solution is found for a model of the frictional ocean surface boundary layer, termed the generalized Ekman (or turbulent thermal wind) balance. The solution admits consideration of general, more physically realistic forms of parameters than previously possible, offering improved physical insight into the underlying dynamics. Closed form solutions are given for the wind-driven flow in the presence of Coriolis–Stokes shear, a result of the surface wave field, and thermal wind shear, arising from a baroclinic pressure gradient, revealing the common underlying physical mechanisms through which they modify currents in the ocean boundary layer. These dynamics are further illustrated by a case study of an idealized two-dimensional front. The solutions, and estimates of the global distribution of the relative influence of surface waves and baroclinic pressure gradients on near-surface ocean currents, emphasize the broad importance of considering ocean sources of shear and physically realistic parameters in the Ekman problem.


2009 ◽  
Vol 39 (5) ◽  
pp. 1077-1096 ◽  
Author(s):  
Gregory P. Gerbi ◽  
John H. Trowbridge ◽  
Eugene A. Terray ◽  
Albert J. Plueddemann ◽  
Tobias Kukulka

Abstract Observations of turbulent kinetic energy (TKE) dynamics in the ocean surface boundary layer are presented here and compared with results from previous observational, numerical, and analytic studies. As in previous studies, the dissipation rate of TKE is found to be higher in the wavy ocean surface boundary layer than it would be in a flow past a rigid boundary with similar stress and buoyancy forcing. Estimates of the terms in the turbulent kinetic energy equation indicate that, unlike in a flow past a rigid boundary, the dissipation rates cannot be balanced by local production terms, suggesting that the transport of TKE is important in the ocean surface boundary layer. A simple analytic model containing parameterizations of production, dissipation, and transport reproduces key features of the vertical profile of TKE, including enhancement near the surface. The effective turbulent diffusion coefficient for heat is larger than would be expected in a rigid-boundary boundary layer. This diffusion coefficient is predicted reasonably well by a model that contains the effects of shear production, buoyancy forcing, and transport of TKE (thought to be related to wave breaking). Neglect of buoyancy forcing or wave breaking in the parameterization results in poor predictions of turbulent diffusivity. Langmuir turbulence was detected concurrently with a fraction of the turbulence quantities reported here, but these times did not stand out as having significant differences from observations when Langmuir turbulence was not detected.


2018 ◽  
Vol 48 (9) ◽  
pp. 2057-2080 ◽  
Author(s):  
Jinliang Liu ◽  
Jun-Hong Liang ◽  
James C. McWilliams ◽  
Peter P. Sullivan ◽  
Yalin Fan ◽  
...  

AbstractA large-eddy simulation (LES) model is configured to investigate the effect of the horizontal (northward) component of Earth’s rotation on upper-ocean turbulence. The focus is on the variability of the effect with latitude/hemisphere in the presence of surface gravity waves and when capped by a stable stratification beneath the surface layer. When is included, the mean flow, turbulence, and vertical mixing depend on the wind direction. The value and effect of are the largest in the tropics and decrease with increasing latitudes. The variability in turbulent flows to wind direction is different at different latitudes and in opposite hemispheres. When limited by stable stratification, the variability in turbulence intensity to wind direction reduces, but the entrainment rate changes with wind direction. In wave-driven Langmuir turbulence, the variability in mean current to wind direction is reduced, but the variability of turbulence to wind direction is evident. When there is wind-following swell, the variability in the mean current to wind direction is further reduced. When there is strong wind-opposing swell so that the total wave forcing is opposite to the wind, the variability in the mean current to wind direction is reduced, but the variability of turbulence to wind direction is enhanced, compared to in Ekman turbulence. The profiles of eddy viscosity, including its shape and its value, show a strong wind direction dependence for both stratified wind-driven and wave-driven Langmuir turbulence. Our study demonstrates that wind direction is an important parameter to upper-ocean mixing, though it is overlooked in existing ocean models.


Sign in / Sign up

Export Citation Format

Share Document