The stirring and mixing of heat and momentum in the ocean surface boundary layer (OSBL) are dominated by 1 to 10 km fluid flows – too small to be resolved in global and regional ocean models. Instead, these processes are parametrized. Two main parametrizations include vertical mixing by surface-forced metre-scale turbulence and overturning by kilometre-scale submesoscale frontal flows and instabilities. In present models, these distinct parametrizations are implemented in tandem, yet ignore meaningful interactions between these two scales that may influence net turbulent fluxes. Using a large-eddy simulation of frontal spin down resolving processes at both scales, this work diagnoses submesoscale and surface-forced turbulence impacts that are the foundation of OSBL parametrizations, following a traditional understanding of these flows. It is shown that frontal circulations act to suppress the vertical buoyancy flux by surface forced turbulence, and that this suppression is not represented by traditional boundary layer turbulence theory. A main result of this work is that current OSBL parametrizations excessively mix buoyancy and overestimate turbulence dissipation rates in the presence of lateral flows. These interactions have a direct influence on the upper ocean potential vorticity and energy budgets with implications for global upper ocean budgets and circulation.
more »
« less
Forcing Space: An Alternative to Regime Diagrams for Predicting Characteristics of Turbulence in the Ocean Surface Mixing Layer
Abstract Various forms of regime diagrams have become an accepted means of identifying the dominant type of forcing of turbulence in the ocean surface layer. However, all of the proposed forms share a number of issues, demonstrated here, that make them an imperfect tool for this purpose. Instead, I suggest a forcing space consisting of surface buoyancy flux (usually dominated by surface heat flux) and a growth rate defined as the inverse of a theoretical time scale for growth of Langmuir circulations in an unstratified water column. Using coastal data, it is demonstrated that, provided forcing conditions are roughly constant for several hours, location in the upper half-plane of this forcing space predicts organizational characteristics of observed turbulence that range in a systematic way between those of “pure” convection and those of full depth Langmuir circulations. In this upper half-plane, where a convective scale velocity exists and the surface Stokes drift velocity can be computed, allowing calculation of a Stokes scale velocity, a linear combination of the two scale velocities provides a consistent estimate of observed rms turbulent vertical velocity. Time dependence is nevertheless a frequent characteristic of ocean surface layer forcing, if only because of the (usually large) diurnal variation in surface heat flux. It is shown that the time scale of response of surface layer turbulence to time variable forcing depends on whether the major change is due to wind/wave or buoyancy forcing. Relevant modeling studies are suggested.
more »
« less
- Award ID(s):
- 1756675
- PAR ID:
- 10364155
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of Physical Oceanography
- Volume:
- 52
- Issue:
- 3
- ISSN:
- 0022-3670
- Format(s):
- Medium: X Size: p. 519-535
- Size(s):
- p. 519-535
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract This work evaluates the fidelity of various upper-ocean turbulence parameterizations subject to realistic monsoon forcing and presents a finite-time ensemble vector (EV) method to better manage the design and numerical principles of these parameterizations. The EV method emphasizes the dynamics of a turbulence closure multimodel ensemble and is applied to evaluate 10 different ocean surface boundary layer (OSBL) parameterizations within a single-column (SC) model against two boundary layer large-eddy simulations (LES). Both LES include realistic surface forcing, but one includes wind-driven shear turbulence only, while the other includes additional Stokes forcing through the wave-average equations that generate Langmuir turbulence. The finite-time EV framework focuses on what constitutes the local behavior of the mixed layer dynamical system and isolates the forcing and ocean state conditions where turbulence parameterizations most disagree. Identifying disagreement provides the potential to evaluate SC models comparatively against the LES. Observations collected during the 2018 monsoon onset in the Bay of Bengal provide a case study to evaluate models under realistic and variable forcing conditions. The case study results highlight two regimes where models disagree 1) during wind-driven deepening of the mixed layer and 2) under strong diurnal forcing.more » « less
-
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
-
Abstract 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
-
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 heatingmore » « less
An official website of the United States government
