Search for: All records

Recent high-resolution large-eddy simulations (LES) of a stable atmospheric boundary layer (SBL) with mesh sizes N=(5123,10243,20483) or mesh spacings ▵=(0.78,0.39,0.2) m are analyzed. The LES solutions are judged to be converged based on the good collapse of vertical profiles of mean winds, temperature, and low-order turbulence moments, i.e., fluxes and variances, with increasing N. The largest discrepancy is in the stably stratified region above the low-level jet. Subfilter-scale (SFS) motions are extracted from the LES with N=20483 and are compared to sonic anemometer fields from the horizontal array turbulence study (HATS) and its sequel over the ocean (OHATS). The results from the simulation and observations are compared using the dimensionless resolution ratio Λw/▵f where ▵f is the filter width and Λw is a characteristic scale of the energy-containing eddies in vertical velocity. The SFS motions from the observations and LES span the ranges 0.1<Λw/▵f<20 and are in good agreement. The small, medium, and large range of Λw/▵f correspond to Reynolds-averaged Navier–Stokes (RANS), the gray zone (a.k.a. “Terra Incognita”), and fine-resolution LES. The gray zone cuts across the peak in the energy spectrum and then flux parameterizations need to be adaptive and account for partially resolved flux but also “stochastic” flux fluctuations that represent the turbulent correlation between the fluctuating rate of strain and SFS flux tensors. LES data with mesh 20483 will be made available to the research community through the web and tools provided by the Johns Hopkins University Turbulence Database.

 

A persistent spatial organization of eddies is identified in the lowest portion of the stably stratified planetary boundary layer. The analysis uses flow realizations from published large-eddy simulations (Sullivan et al. in J Atmos Sci 73(4):1815–1840, 2016) ranging in stability from near-neutral to almost z-less stratification. The coherent turbulent structure is well approximated as a series of uniform momentum zones (UMZs) and uniform temperature zones (UTZs) separated by thin layers of intense gradients that are significantly greater than the mean. This pattern yields stairstep-like instantaneous flow profiles whose shape is distinct from the mean profiles that emerge from long-term averaging. However, the scaling of the stairstep organization is closely related to the resulting mean profiles. The differences in velocity and temperature across the thin gradient layers remain proportional to the surface momentum and heat flux conditions regardless of stratification. The vertical thickness of UMZs and UTZs is proportional to height above the surface for near-neutral and weak stratification, but becomes thinner and less dependent on height as the stability increases. Deviations from the logarithmic mean profiles for velocity and temperature observed under neutral conditions are therefore predominately due to the reduction in eddy size with increasing stratification, which is empirically captured by existing Monin–Obukhov similarity relations for momentum and heat. The zone properties are additionally used to explain trends in the turbulent Prandtl number, thus providing a connection between the eddy organization, mean profiles, and turbulent diffusivity in stably stratified conditions.

 

Air–sea momentum and scalar fluxes are strongly influenced by the coupling dynamics between turbulent winds and a spectrum of waves. Because direct field observations are difficult, particularly in high winds, many modeling and laboratory studies have aimed to elucidate the impacts of the sea state and other surface wave features on momentum and energy fluxes between wind and waves as well as on the mean wind profile and drag coefficient. Opposing wind is common under transient winds, for example, under tropical cyclones, but few studies have examined its impacts on air–sea fluxes. In this study, we employ a large-eddy simulation for wind blowing over steep sinusoidal waves of varying phase speeds, both following and opposing wind, to investigate impacts on the mean wind profile, drag coefficient, and wave growth/decay rates. The airflow dynamics and impacts rapidly change as the wave age increases for waves following wind. However, there is a rather smooth transition from the slowest waves following wind to the fastest waves opposing wind, with gradual enhancement of a flow perturbation identified by a strong vorticity layer detached from the crest despite the absence of apparent airflow separation. The vorticity layer appears to increase the effective surface roughness and wave form drag (wave attenuation rate) substantially for faster waves opposing wind.

Surface waves increase friction at the sea surface and modify how wind forces upper-ocean currents and turbulence. Therefore, it is important to include effects of different wave conditions in weather and climate forecasts. We aim to inform more accurate forecasts by investigating wind blowing over waves propagating in the opposite direction using large-eddy simulation. We find that when waves oppose wind, they decay as expected, but also increase the surface friction much more drastically than when waves follow wind. This finding has important implications for how waves opposing wind are represented as a source of surface friction in forecast models.

 

The coupled dynamics of turbulent airflow and a spectrum of waves are known to modify air–sea momentum and scalar fluxes. Waves traveling at oblique angles to the wind are common in the open ocean, and their effects may be especially relevant when constraining fluxes in storm and tropical cyclone conditions. In this study, we employ large-eddy simulation for airflow over steep, strongly forced waves following and opposing oblique wind to elucidate its impacts on the wind speed magnitude and direction, drag coefficient, and wave growth/decay rate. We find that oblique wind maintains a signature of airflow separation while introducing a cross-wave component strongly modified by the waves. The directions of mean wind speed and mean wind shear vary significantly with height and are misaligned from the wind stress direction, particularly toward the surface. As the oblique angle increases, the wave form drag remains positive, but the wave impact on the equivalent surface roughness (drag coefficient) rapidly decreases and becomes negative at large angles. Our findings have significant implications for how the sea-state-dependent drag coefficient is parameterized in forecast models. Our results also suggest that wind speed and wind stress measurements performed on a wave-following platform can be strongly contaminated by the platform motion if the instrument is inside the wave boundary layer of dominant waves.

Surface waves increase friction at the sea surface and modify how wind forces upper-ocean currents and turbulence. Therefore, it is important to include effects of different wave conditions in weather and climate forecasts. We aim to inform more accurate forecasts by investigating wind blowing over waves propagating in oblique directions using large-eddy simulation. We find that waves traveling at a 45° angle or larger to the wind grow as expected, but do not increase or even decrease the surface friction felt by the wind—a surprising result that has significant implications for how oblique wind-waves are represented as a source of surface friction in forecast models.

 

Ocean fronts are an important submesoscale feature, yet frontogenesis theory often neglects turbulence – even parameterized turbulence – leaving theory lacking in comparison with observations and models. A perturbation analysis is used to include the effects of eddy viscosity and diffusivity as a first-order correction to existing strain-induced inviscid, adiabatic frontogenesis theory. A modified solution is obtained by using potential vorticity and surface conditions to quantify turbulent fluxes. It is found that horizontal viscosity and diffusivity tend to be readily frontolytic – reducing frontal tendency to negative values under weakly non-conservative perturbations and opposing or reversing front sharpening, whereas vertical viscosity and diffusivity tend to only weaken frontogenesis by slowing the rate of sharpening of the front even under strong perturbations. During late frontogenesis, vertical diffusivity enhances the rate of frontogenesis, although perturbation theory may be inaccurate at this stage. Surface quasi-geostrophic theory – neglecting all injected interior potential vorticity – is able to describe the first-order correction to the along-front velocity and ageostrophic overturning circulation in most cases. Furthermore, local conditions near the front maximum are sufficient to reconstruct the modified solution of both these fields. 

Large-eddy simulation (LES) is used to model turbulent winds in a nominally neutral atmospheric boundary layer at varying mesh resolutions. The boundary layer is driven by wind shear with zero surface heat flux and is capped by a stable inversion. Because of entrainment the boundary layer is in a weakly stably stratified regime. The simulations use meshes varying from 128²× 64 to 1024²× 512 grid points in a fixed computational domain of size (2560, 2560, 896) m. The subgrid-scale (SGS) parameterizations used in the LES vary with the mesh spacing. Low-order statistics, spectra, and structure functions are compared on the different meshes and are used to assess grid convergence in the simulations. As expected, grid convergence is primarily achieved in the middle of the boundary layer where there is scale separation between the energy-containing and dissipative eddies. Near the surface second-order statistics do not converge on the meshes studied. The analysis also highlights differences between one-dimensional and two-dimensional velocity spectra; differences are attributed to sampling errors associated with aligning the horizontal coordinates with the vertically veering mean wind direction. Higher-order structure functions reveal non-Gaussian statistics on all scales, but are highly dependent on the mesh resolution. A generalized logarithmic law and a k⁻¹spectral scaling regime are identified with mesh-dependent parameters in agreement with previously published results.

 

A common technique for estimating the sea surface generation functions of spray and aerosols is the so-called flux–profile method, where fixed-height concentration measurements are used to infer fluxes at the surface by assuming a form of the concentration profile. At its simplest, this method assumes a balance between spray emission and deposition, and under these conditions the concentration profile follows a power-law shape. It is the purpose of this work to evaluate the influence of waves on this power-law theory, as well as investigate its applicability over a range of droplet sizes. Large-eddy simulations combined with Lagrangian droplet tracking are used to resolve the turbulent transport of spray droplets over moving, monochromatic waves at the lower surface. The wave age and the droplet diameter are varied, and it is found that droplets are highly influenced both by their inertia (i.e., their inability to travel exactly with fluid streamlines) and the wave-induced turbulence. Deviations of the vertical concentration profiles from the power-law theory are found at all wave ages and for large droplets. The dynamics of droplets within the wave boundary layer alter their net vertical fluxes, and as a result, estimates of surface emission based on the flux–profile method can yield significant errors. In practice, the resulting implication is that the flux–profile method may unsuitable for large droplets, and the combined effect of inertia and wave-induced turbulence is responsible for the continued spread in their surface source estimates.