Explore the vast range of titles available.

We present a review of existing wind–wave coupling models and parameterizations used for large-eddy simulation of the marine atmospheric boundary layer. The models are classified into two main categories: (i) the wave-phase-averaged, sea surface–roughness models and (ii) the wave-phase-resolved models. Both categories are discussed from their implementation, validity, and computational efficiency viewpoints, with emphasis given on their applicability in offshore wind energy problems. In addition to the various models discussed, a review of laboratory-scale and field-measurement databases is presented thereafter. The majority of the presented data have been gathered over many decades of studying air–sea interaction phenomena, with the most recent ones compiled to reflect an offshore wind energy perspective. Both provide valuable data for model validation. We also discuss the modeling knowledge gaps and computational challenges ahead.

Turbulent flow in a weakly convective marine atmospheric boundary layer (MABL) driven by geostrophic winds U g = 10 m s −1 and heterogeneous sea surface temperature (SST) is examined using fine-mesh large-eddy simulation (LES). The imposed SST heterogeneity is a single-sided warm or cold front with temperature jumps Δ θ = (2, −1.5) K varying over a horizontal distance between [0.1, −6] km characteristic of an upper-ocean mesoscale or submesoscale regime. A Fourier-fringe technique is implemented in the LES to overcome the assumptions of horizontally homogeneous periodic flow. Grid meshes of 2.2 × 10 9 points with fine-resolution (horizontal, vertical) spacing ( δx = δy , δz ) = (4.4, 2) m are used. Geostrophic winds blowing across SST isotherms generate secondary circulations that vary with the sign of the front. Warm fronts feature overshoots in the temperature field, nonlinear temperature and momentum fluxes, a local maximum in the vertical velocity variance, and an extended spatial evolution of the boundary layer with increasing distance from the SST front. Cold fronts collapse the incoming turbulence but leave behind residual motions above the boundary layer. In the case of a warm front, the internal boundary layer grows with downstream distance conveying the surface changes aloft and downwind. SST fronts modify entrainment fluxes and generate persistent horizontal advection at large distances from the front.

Sea surface temperature (SST) is characterized by abundant warm and cold structures that influence the overlying atmospheric boundary layer dynamics through two different mechanisms. First, turbulence and large eddies in the lower troposphere are affected by atmospheric stability, which can be modified by local SST, resulting in enhanced vertical mixing and larger surface winds over warmer waters. Second, the thermodynamic adjustment of air density to the underlying SST structures and the subsequent changes in atmospheric pressure drive secondary circulations. This paper aims to disentangle the effects of these processes and explore the environmental conditions that favor them. Two main environmental variables are considered: the large-scale air–sea temperature difference (proxy for stability) and wind speed. Using 5 years of daily reanalyses data, we investigate the 10-m wind response to SST structures. Based on linear regression between wind divergence and SST derivatives, we show that both mechanisms operate over a large spectrum of conditions. Ten-meter wind divergence is strongly impacted by the local SST via its effect on vertical mixing for midwind regimes in slightly unstable to near-neutral conditions, whereas the secondary circulation is important in two distinct regimes: low wind speed with a slightly unstable air column and high background wind speed with a very unstable air column. The first regime is explained by the prolonged Lagrangian time that the air parcel stays over an SST structure while the second one is related to strong heat fluxes at the air–sea interface, which greatly modify the marine atmospheric boundary layer properties. Location and frequency of the environmentally favorable conditions are discussed, as well as the response in low-cloud cover and rainfall.

A mechanism for precipitation enhancement in squall lines moving over mountainous coastal regions is quantified through idealized numerical simulations. Storm intensity and precipitation peak over the sloping terrain as storms descend from an elevated plateau toward the coastline and encounter the marine atmospheric boundary layer (MABL). Storms are most intense as they encounter the deepest MABLs. As the descending storm outflow collides with a moving MABL (sea breeze), surface and low-level air parcels initially accelerate upward, though their ultimate trajectory is governed by the magnitude of the negative nonhydrostatic inertial pressure perturbation behind the cold pool leading edge. For shallow MABLs, the baroclinic gradient across the gust front generates large horizontal vorticity, a low-level negative pressure perturbation, and thus a downward acceleration of air parcels following their initial ascent. A deep MABL reduces the baroclinically generated vorticity, leading to a weaker pressure perturbation and minimal downward acceleration, allowing air to accelerate into a storm’s updraft. Once storms move away from the terrain base and over the full depth of the MABLs, storms over the deepest MABLs decay most rapidly, while those over the shallowest MABLs initially intensify. Though elevated ascent exists above all MABLs, the deepest MABLs substantially reduce the depth of the high- θ e layer above the MABLs and limit instability. This relationship is insensitive to MABL temperature, even though surface-based ascent is present for the less cold MABLs, the MABL thermal deficit is smaller, and convective available potential energy (CAPE) is higher.

Numerical prediction of the interactions between wind and ocean waves is essential for climate modeling and a wide range of offshore operations. Large Eddy Simulation (LES) of the marine atmospheric boundary layer is a practical numerical predictive tool but requires parameterization of surface fluxes at the air–water interface. Current momentum flux parameterizations primarily use wave-phase adapting computational grids, incurring high computational costs, or use an equilibrium model based on Monin–Obukhov similarity theory for rough surfaces that cannot resolve wave phase information. To include wave phase-resolving physics at a cost similar to the equilibrium model, the MOving Surface Drag (MOSD) model is introduced. It assumes ideal airflow over locally piece-wise planar representations of moving water wave surfaces. Horizontally unresolved interactions are still modeled using the equilibrium model. Validation against experimental and numerical datasets with known monochromatic waves demonstrates the robustness and accuracy of the model in representing wave-induced impacts on mean velocity and Reynolds stress profiles. The model is formulated to be applicable to a broad range of wave fields and its ability to represent cross-swell and multiple wavelength cases is illustrated. Additionally, the model is applied to LES of a laboratory-scale fixed-bottom offshore wind turbine model, and the results are compared with wind tunnel experimental data. The LES with the MOSD model shows good agreement in wind–wave–wake interactions and phase-dependent physics at a low computational cost. The model’s simplicity and minimal computational needs make it valuable for studying turbulent atmospheric-scale flows over the sea, particularly in offshore wind energy research.

Lagrangian stochastic models (LSM) are widely used to model the dispersion of sea spray droplets injected from the water surface into the marine atmospheric boundary layer (MABL) and for evaluation of the spray impact on the exchange fluxes between the atmosphere and the ocean. While moving through the MABL the droplets pass through the region of high gradients of air velocity, temperature and humidity occurring in the vicinity of the air–water interface. In this case, the applicability of LSMs constructed under the assumption of weakly inhomogeneous flows is questionable. In this work, we develop a Lagrangian stochastic model taking into account the strongly inhomogeneous structure of the airflow in MABL and, in particular, the anisotropy of turbulence dissipation rate. The model constants and the diffusion matrix coefficients are calibrated by comparison of the LSM prediction for the profiles of droplet concentration and the exchange fluxes of sensible and latent heat against the results of direct numerical simulation of turbulent, droplet-laden airflow over a waved water surface.

The Eulerian multifluid mathematical model is developed to describe the marine atmospheric boundary layer laden with sea spray under the high-wind condition of a hurricane. The model considers spray and air as separate continuous interacting turbulent media and employs the multifluid E –ϵ closure. Each phase is described by its own set of coupled conservation equations and characterized by its own velocity. Such an approach enables us to accurately quantify the interaction between spray and air and pinpoint the effect of spray on the vertical momentum transport much more precisely than could be done with traditional mixture-type approaches. The model consistently quantifies the effect of spray inertia and the suppression of air turbulence due to two different mechanisms: the turbulence attenuation, which results from the inability of spray droplets to fully follow turbulent fluctuations, and the vertical transport of spray against the gravity by turbulent eddies. The results of numerical and asymptotic analyses show that the turbulence suppression by spray overpowers its inertia several meters above wave crests, resulting in a noticeable wind acceleration and the corresponding reduction of the drag coefficient from the reference values for a spray-free atmosphere. This occurs at much lower than predicted previously spray volume fraction values of ∼10 −5 . The falloff of the drag coefficient from its reference values is more strongly pronounced at higher altitudes. The drag coefficient reaches its maximum at spray volume fraction values of ∼10 −4 , which is several times smaller than predicted by mixture-type models.

Air–sea interaction processes over the Gulf Stream have received particular attention over the last decade. It has been shown that sea surface temperature (SST) gradients over the Gulf Stream can alter the near-surface wind divergence through changes in the marine atmospheric boundary layer (MABL). Two mechanisms have been proposed to explain the response: the vertical mixing mechanism (VMM) and the pressure adjustment mechanism (PAM). However, their respective contribution is still under debate. It has been argued that the synoptic perturbations over the Gulf Stream can provide more insight into the MABL response to SST fronts. We analyze the VMM and PAM under different atmospheric conditions obtained from a classification method that is based on the deciles of the statistical distribution of winter turbulent heat fluxes over the Gulf Stream. The lowest deciles are associated with weak air–sea interactions and anticyclonic atmospheric circulation over the Gulf Stream, whereas the highest deciles are related to strong air–sea interactions and a cyclonic circulation. Our analysis includes the low- and high-resolution versions of the ARPEGEv6 atmospheric model forced by observed SST, and the recently released ERA5 global reanalysis. We find that the occurrence of anticyclonic and cyclonic perturbations associated with different anomalous wind regimes can locally modulate the activation of the VMM and the PAM. In particular, the PAM is predominant in anticyclonic conditions, whereas both mechanisms are equally present in most of the cyclonic conditions. Our results highlight the role of the atmospheric circulation and associated anomalous winds in the location, strength, and occurrence of both mechanisms.

The turbulent structure within the marine atmospheric boundary layer is investigated based on four levels of observations at a fixed marine platform. During and before a cold front, the ocean surface is dominated by wind sea and swell waves, respectively, affording the opportunity to test the theory recently proposed in laboratory experiments or for flat land surfaces. The results reveal that the velocity spectra follow a k −1 law within the intermediate wavenumber ( k ) range immediately below inertial subrange during the cold front. A logarithmic height dependence of the horizontal velocity variances is also observed above the height of 20 m, while the vertical velocity variances increase with increasing height following a power law of 2/3. These features confirm the attached eddy model and the “top-down model” of turbulence over the ocean surface. However, the behavior of velocity variances under swell conditions is much different from those during the cold front, although a remarkable k −1 law can be observed in the velocity spectra. The quadrant analysis of the momentum flux also shows a significantly different result for the two conditions.

A new functional form of the neutral drag coefficient for moderate to high wind speeds in the marine atmospheric boundary layer for a range of field measurements as reported in the literature is proposed. This new form is found to describe a wide variety of measurements recorded in the open ocean, coast, fetch-limited seas, and lakes, with almost one and the same set of parameters. This is the result of a reanalysis of the definition of the drag coefficient in the marine boundary layer, which finds that a constant is missing from the traditional definition of the drag coefficient. The constant arises because the neutral friction velocity over water surfaces is not directly proportional to the 10-m wind speed, a consequence of the transition to rough flow at low wind speeds. Within the rough flow regime, the neutral friction velocity is linearly dependent on the 10-m wind speed; consequently, within this rough regime, the new definition of the drag coefficient is not a function of the wind speed. The magnitude of the new definition of the neutral drag coefficient represents an upper limit to the magnitude of the traditional definition.