Abstract

Measurements of velocity, density, and pressure gradient in the lower Hudson River estuary were used to quantify the dominant terms in the momentum equation and to characterize their variations at tidal and spring–neap timescales. The vertical momentum flux (assumed to be due mainly to turbulent shear stress) was estimated indirectly, based on the residual from the acceleration and pressure gradient terms. The indirect estimates of stress compared favorably to bottom stress estimates using a quadratic drag law, supporting the hypothesis that the tidal momentum equation involves a local balance between tidal acceleration, pressure gradient, and stress divergence.

Estimates of eddy viscosity indicated that there was significant tidal asymmetry, with flood tide values exceeding ebb values by a factor of 2. As a consequence of the asymmetry, the vertical structure of the tidally averaged stress bore no resemblance to the tidally averaged shear. In spite of the asymmetry of vertical mixing, the tidally averaged, estuarine circulation was found to depend simply on the intensity of bottom turbulence, which could be parameterized by a Rayleigh drag formulation based on the tidal velocity magnitude and the tidally averaged near-bottom flow. This seemingly paradoxical result indicates that the estuarine circulation can be modeled without detailed knowledge of the effective eddy viscosity, only requiring an estimate of the bottom drag coefficient, the tidal forcing conditions, and the baroclinic pressure gradient. A notable characteristic of this solution is an inverse dependence of the estuarine circulation on the amplitude of the tides.

Abstract

Measurements just beneath the ocean surface demonstrate that the primary mechanism by which energy from breaking waves is transmitted into the water column is through the work done by the covariance of turbulent pressure and velocity fluctuations. The convergence in the vertical transport of turbulent kinetic energy (TKE) balances the dissipation rate of TKE at first order and is nearly an order of magnitude greater than the sum of the integrated Eulerian and Stokes shear production. The measured TKE transport is consistent with a simple conceptual model that assumes roughly half of the surface flux of TKE by wave breaking is transmitted to depths greater than the significant wave height. During conditions when breaking waves are inferred, the direction of momentum flux is more aligned with the direction of wave propagation than with the wind direction. Both the energy and momentum fluxes occur at frequencies much lower than the wave band, consistent with the time scales associated with wave breaking. The largest instantaneous values of momentum flux are associated with strong downward vertical velocity perturbations, in contrast to the pressure work, which is associated with strong drops in pressure and upward vertical velocity perturbations.

Abstract

Langmuir circulation (LC) is a turbulent upper-ocean process driven by wind and surface waves that contributes significantly to the transport of momentum, heat, and mass in the oceanic surface layer. The authors have previously performed a direct comparison of large-eddy simulations and observations of the upper-ocean response to a wind event with rapid mixed layer deepening. The evolution of simulated crosswind velocity variance and spatial scales, as well as mixed layer deepening, was only consistent with observations if LC effects are included in the model. Based on an analysis of these validated simulations, in this study the fundamental differences in mixing between purely shear-driven turbulence and turbulence with LC are identified. In the former case, turbulent kinetic energy (TKE) production due to shear instabilities is largest near the surface, gradually decreasing to zero near the base of the mixed layer. This stands in contrast to the LC case in which at middepth range TKE production can be dominated by Stokes drift shear. Furthermore, the Eulerian mean vertical shear peaks near the base of the mixed layer so that TKE production by mean shear flow is elevated there. LC transports horizontal momentum efficiently downward leading to an along-wind velocity jet below LC downwelling regions at the base of the mixed layer. Locally enhanced vertical shear instabilities as a result of this jet efficiently erode the thermocline. In turn, enhanced breaking internal waves inject cold deep water into the mixed layer, where LC currents transport temperature perturbation advectively. Thus, LC and locally generated shear instabilities work intimately together to facilitate strongly the mixed layer deepening process.

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.

Abstract

This study makes direct measurements of turbulent fluxes in the mixed layer in order to close heat and momentum budgets across the air–sea interface and to assess the ability of rigid-boundary turbulence models to predict mean vertical gradients beneath the ocean’s wavy surface. Observations were made at 20 Hz at nominal depths of 2.2 and 1.7 m in ∼16 m of water. A new method is developed to estimate the fluxes and the length scales of dominant flux-carrying eddies from cospectra at frequencies below the wave band. The results are compared to independent estimates of those quantities, with good agreement between the two sets of estimates. The observed temperature gradients were smaller than predicted by standard rigid-boundary closure models, consistent with the suggestion that wave breaking and Langmuir circulation increase turbulent diffusivity in the upper ocean. Similarly, the Monin–Obukhov stability function ϕ_h was smaller in the authors’ measurements than the stability functions used in rigid-boundary applications of the Monin–Obukhov similarity theory. The dominant horizontal length scales of flux-carrying turbulent eddies were found to be consistent with observations in the bottom boundary layer of the atmosphere and from laboratory experiments in three ways: 1) in statically unstable conditions, the eddy sizes scaled linearly with distance to the boundary; 2) in statically stable conditions, length scales decreased with increasing downward buoyancy flux; and 3) downwind length scales were larger than crosswind length scales.

The Office of Naval Research's Coupled Boundary Layers and Air–Sea Transfer (CBLAST) program is being conducted to investigate the processes that couple the marine boundary layers and govern the exchange of heat, mass, and momentum across the air–sea interface. CBLAST-LOW was designed to investigate these processes at the low-wind extreme where the processes are often driven or strongly modulated by buoyant forcing. The focus was on conditions ranging from negligible wind stress, where buoyant forcing dominates, up to wind speeds where wave breaking and Langmuir circulations play a significant role in the exchange processes. The field program provided observations from a suite of platforms deployed in the coastal ocean south of Martha's Vineyard. Highlights from the measurement campaigns include direct measurement of the momentum and heat fluxes on both sides of the air–sea interface using a specially constructed Air–Sea Interaction Tower (ASIT), and quantification of regional oceanic variability over scales of O(1–10⁴ mm) using a mesoscale mooring array, aircraft-borne remote sensors, drifters, and ship surveys. To our knowledge, the former represents the first successful attempt to directly and simultaneously measure the heat and momentum exchange on both sides of the air–sea interface. The latter provided a 3D picture of the oceanic boundary layer during the month-long main experiment. These observations have been combined with numerical models and direct numerical and large-eddy simulations to investigate the processes that couple the atmosphere and ocean under these conditions. For example, the oceanic measurements have been used in the Regional Ocean Modeling System (ROMS) to investigate the 3D evolution of regional ocean thermal stratification. The ultimate goal of these investigations is to incorporate improved parameterizations of these processes in coupled models such as the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) to improve marine forecasts of wind, waves, and currents.