Explore the vast range of titles available.

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.

A massively parallel large-eddy simulation (LES) code for planetary boundary layers (PBLs) that utilizes pseudospectral differencing in horizontal planes and solves an elliptic pressure equation is described. As an application, this code is used to examine the numerical convergence of the three-dimensional time-dependent simulations of a weakly sheared daytime convective PBL on meshes varying from 323 to 10243 grid points. Based on the variation of the second-order statistics, energy spectra, and entrainment statistics, LES solutions converge provided there is adequate separation between the energy-containing eddies and those near the filter cutoff scale. For the convective PBL studied, the majority of the low-order moment statistics (means, variances, and fluxes) become grid independent when the ratio zi/(CsΔf) > 310, where zi is the boundary layer height, Δf is the filter cutoff scale, and Cs is the Smagorinsky constant. In this regime, the spectra show clear Kolmogorov inertial subrange scaling. The bulk entrainment rate determined from the time variation of the boundary layer height we = dzi/dt is a sensitive measure of the LES solution convergence; we becomes grid independent when the vertical grid resolution is able to capture both the mean structure of the overlying inversion and the turbulence. For all mesh resolutions used, the vertical temperature flux profile varies linearly over the interior of the boundary layer and the minimum temperature flux is approximately −0.2 of the surface heat flux. Thus, these metrics are inadequate measures of solution convergence. The variation of the vertical velocity skewness and third-order moments expose the LES’s sensitivity to grid resolution.

Accurate predictions of the sea state–dependent air–sea momentum flux require a thorough understanding of the wave boundary layer turbulence over surface waves. A set of momentum and energy equations is derived to formulate and analyze wave boundary layer turbulence. The equations are written in wave-following coordinates, and all variables are decomposed into horizontal mean, wave fluctuation, and turbulent fluctuation. The formulation defines the wave-induced stress as a sum of the wave fluctuation stress (because of the fluctuating velocity components) and a pressure stress (pressure acting on a tilted surface). The formulations can be constructed with different choices of mapping. Next, a large-eddy simulation result for wind over a sinusoidal wave train under a strongly forced condition is analyzed using the proposed formulation. The result clarifies how surface waves increase the effective roughness length and the drag coefficient. Specifically, the enhanced wave-induced stress close to the water surface reduces the turbulent stress (satisfying the momentum budget). The reduced turbulent stress is correlated with the reduced viscous dissipation rate of the turbulent kinetic energy. The latter is balanced by the reduced mean wind shear (satisfying the energy budget), which causes the equivalent surface roughness to increase. Interestingly, there is a small region farther above where the turbulent stress, dissipation rate, and mean wind shear are all enhanced. The observed strong correlation between the turbulent stress and the dissipation rate suggests that existing turbulence closure models that parameterize the latter based on the former are reasonably accurate.

Momentum and scalar transport in the marine atmospheric boundary layer (MABL) is driven by a turbulent mix of winds, buoyancy, and surface gravity waves. To investigate the interaction between these processes, a large-eddy simulation (LES) model is developed with the capability to impose a broadband spectrum of time-varying finite-amplitude surface waves at its lower boundary. The LES model adopts a Boussinesq flow model and integrates the governing equations on a time-varying, surface-fitted, nonorthogonal mesh using cell-centered variables with special attention paid to the solution of the pressure Poisson equation near the wavy boundary. Weakly unstable MABLs are simulated with geostrophic winds increasing from 5 to 25 m s−1 and wave age varying from swell-dominated to wind-wave equilibrium. The simulations illustrate cross-scale coupling as wave-impacted near-surface turbulence transitions into shear-convective rolls with increasing distance from the water. In a regime with swell, low winds, and weak heating, wave-induced vertical velocity and pressure signals are readily observed well above the standard reference height ζa = 10 m. At wind-wave equilibrium, the small-scale wave-induced signals are detectable only near the water surface. Below ζa, a nearly-constant-flux layer is observed where the momentum flux carried by turbulence, form stress, and subgrid-scale motions shifts with varying wave age and distance above the water. The spectral content of the surface form stress is wave-age dependent, especially at low wavenumbers. The LES wind profiles deviate from Monin–Obukhov similarity theory in nonequilibrium wind-wave conditions, and entrainment is greatly enhanced by shear-induced engulfment events.

Turbulent flow over strongly forced steep steady and unsteady waves is simulated using large-eddy simulation (LES) with time t and space x varying wave height h ( x , t ) imposed as a lower boundary condition. With steady waves, h ( x , t ) is based on measurements of incipient and active breaking waves collected in a wind-wave flume, while a numerical wave code is used to generate an unsteady evolving wave packet (group). Highly intermittent airflow separation is found in the simulations, and the results suggest separation near a wave crest occurs prior to the onset of wave breaking. The form (pressure) drag is most sensitive to the wave slope, and the form drag can contribute as much as 74% to the total stress. Wind and scalar profiles from the LES display log-linear variations above the wave surface; the LES wind profiles are in good agreement with the measurements. The momentum roughness increases as the water surface changes from wind ripples to incipient breaking to active breaking. However, the scalar roughness decreases as the wave surface becomes rougher. This highlights major differences in momentum and scalar transport over a rough wavy surface. For a rapidly evolving, strongly forced wave group, the form drag is highly correlated with the wave slope, and intermittent separation is found early in the packet evolution when the local wave slope −∂ h /∂ x ( x , t ) ≥ 0.22. The packet root-mean-square wave slope is 0.084, but the form drag fraction is 2.4 times larger than a comparably forced steady wave. Thus, a passing wave group can induce unsteadiness in the wind stress.

This title is part of UC Press's Voices Revived program, which commemorates University of California Press's mission to seek out and cultivate the brightest minds and give them voice, reach, and impact. Drawing on a backlist dating to 1893, Voices Revived makes high-quality, peer-reviewed scholarship accessible once again using print-on-demand technology. This title was originally published in 1987.

The turbulent mixing in thin ocean surface boundary layers (OSBL), which occupy the upper 100 m or so of the ocean, control the exchange of heat and trace gases between the atmosphere and ocean. Here we show that current parameterizations of this turbulent mixing lead to systematic and substantial errors in the depth of the OSBL in global climate models, which then leads to biases in sea surface temperature. One reason, we argue, is that current parameterizations are missing key surface‐wave processes that force Langmuir turbulence that deepens the OSBL more rapidly than steady wind forcing. Scaling arguments are presented to identify two dimensionless parameters that measure the importance of wave forcing against wind forcing, and against buoyancy forcing. A global perspective on the occurrence of wave‐forced turbulence is developed using re‐analysis data to compute these parameters globally. The diagnostic study developed here suggests that turbulent energy available for mixing the OSBL is under‐estimated without forcing by surface waves. Wave‐forcing and hence Langmuir turbulence could be important over wide areas of the ocean and in all seasons in the Southern Ocean. We conclude that surface‐wave‐forced Langmuir turbulence is an important process in the OSBL that requires parameterization. Key Points Climate models have biases in the depth of the ocean surface boundary layer Langmuir turbulence is a key process mixing the ocean surface boundary layer Langmuir turbulence deepens the layer more quickly than wind‐forced turbulence

Large-eddy simulation of atmospheric boundary layers interacting with a coupled and resolved plant canopy reveals the influence of atmospheric stability variations from neutral to free convection on canopy turbulence. The design and implementation of a new multilevel canopy model is presented. Instantaneous fields from the simulations show that organized motions on the scale of the atmospheric boundary layer (ABL) depth bring high momentum down to canopy top, locally modulating the vertical shear of the horizontal wind. The evolution of these ABL-scale structures with increasing instability and their impact on vertical profiles of turbulence moments and integral length scales within and above the canopy are discussed. Linkages between atmospheric turbulence and biological control impact horizontal scalar source distributions. Decreasing spatial correlation between momentum and scalar fluxes with increasing instability results from ABL-scale structures spatially segregating momentum and scalar exchange at canopy top. In combination, these results suggest the need for roughness sublayer parameterizations to incorporate an additional length or time scale reflecting the influence of ABL-scale organized motions.

Generally, ocean waves are thought to act as a drag on the surface wind so that momentum is transferred downward, from the atmosphere into the waves. Recent observations have suggested that when long wavelength waves—which are characteristic of remotely generated swell—propagate faster than the surface wind, momentum can also be transferred upward. This upward momentum transfer acts to accelerate the near-surface wind, resulting in a low-level wave-driven wind jet. Previous studies have suggested that the sign reversal of the momentum flux is well predicted by the inverse wave age, the ratio of the surface wind speed to the speed of the waves at the peak of the spectrum. Data from the 40-yr ECMWF Re-Analysis (ERA-40) have been used here to calculate the global distribution of the inverse wave age to determine whether there are regions of the ocean that are usually in the wind-driven wave regime and others that are generally in the wave-driven wind regime. The wind-driven wave regime is found to occur most often in the midlatitude storm tracks where wind speeds are generally high. The wave-driven wind regime is found to be prevalent in the tropics where wind speeds are generally light and swell can propagate from storms at higher latitudes. The inverse wave age is also a useful indicator of the degree of coupling between the local wind and wave fields. The climatologies presented emphasize the nonequilibrium that exists between the local wind and wave fields and highlight the importance of swell in the global oceans.

The problem is posed and solved for the oceanic surface boundary layer in the presence of wind stress, stable density stratification, equilibrium wind-waves, and remotely generated swell-waves. The addition of swell causes an amplification of the Lagrangian-mean current and rotation toward the swell-wave direction, a fattening of the Ekman velocity spiral and associated vertical Reynolds stress profile, an amplification of the inertial current response, an enhancement of turbulent variance and buoyancy entrainment rate from the pycnocline, and—for very large swell—an upscaling of the coherent Langmuir circulation patterns. Implications are discussed for the parameterization of Langmuir turbulence influences on the mean current profile and the material entrainment rate in oceanic circulation models. In particular, even though the turbulent kinetic energy monotonically increases with wave amplitude inversely expressed by the turbulent Langmuir number La, the Lagrangian shear eddy viscosity profile κL(z) is a nonmonotonic function of La, first increasing with increasing wave amplitude up to approximately the wind-wave equilibrium level, then decreasing with additional swell-wave amplitude. In contrast, the pycnocline entrainment rate is a monotonic function ~La−2.