Abstract
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.
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.
Abstract
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 super(-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 zeta sub(a) = 10 m. At wind-wave equilibrium, the small-scale wave-induced signals are detectable only near the water surface. Below zeta sub(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.
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
Abstract
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.
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
$\unicodeSTIX{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.
Abstract Familial hypercholesterolaemia (FH) is a dominantly inherited disorder present from birth that markedly elevates plasma low-density lipoprotein (LDL) cholesterol and causes premature ...coronary heart disease. There are at least 20 million people with FH worldwide, but the majority remain undetected and current treatment is often suboptimal. To address this major gap in coronary prevention we present, from an international perspective, consensus-based guidance on the care of FH. The guidance was generated from seminars and workshops held at an international symposium. The recommendations focus on the detection, diagnosis, assessment and management of FH in adults and children, and set guidelines for clinical purposes. They also refer to best practice for cascade screening and risk notifying and testing families for FH, including use of genetic testing. Guidance on treatment is based on risk stratification, management of non-cholesterol risk factors, and safe and effective use of LDL lowering therapies. Recommendations are given on lipoprotein apheresis. The use of emerging therapies for FH is also foreshadowed. This international guidance acknowledges evidence gaps, but aims to make the best use of contemporary practice and technology to achieve the best outcomes for the care of FH. It should accordingly be employed to inform clinical judgement and be adjusted for country-specific and local health care needs and resources.
A Global Climatology of Wind-Wave Interaction HANLEY, Kirsty E; BELCHER, Stephen E; SULLIVAN, Peter P
Journal of physical oceanography,
06/2010, Volume:
40, Issue:
6
Journal Article
Peer reviewed
Open access
Abstract
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.
Through interactions between the modes of “waves” and “vortices”, stably stratified turbulence exhibits characteristic features both in 2D and 3D. Using DNS of the Navier–Stokes equations coupled to ...an equation for temperature (the Bousinessq system), we investigate dynamical properties of stably stratified turbulence. After reviewing the importance of three characteristic length scales and their relations in stably stratified turbulence, we present numerical results showing how structures and spectra develop with the growth of the length scales. It is demonstrated that the temperature fluctuations make sharp fronts vertically which result in non-symmetric PDFs of the vertical derivative ∂zθ.