On the Variability of the Fluxes of Momentum and Sensible Heat

Direct measurements of the air–sea turbulent fluxes of momentum and heat, along with surface currents, waves and supporting meteorological variables, were acquired during a recent field campaign. Surface currents, measured from a very high frequency radar, were found to steer the stress away from the mean wind direction. Although this effect has been reported in a recent scatterometer study, this is the first time it has been observed in an in situ study with co-located flux, wind and surface current measurements. Data collected during a week of stationary conditions are used to investigate and quantify the sampling variability of the air–sea fluxes of momentum and sensible heat.     

Effects of Roll Vortices on Turbulent Fluxes in the Hurricane Boundary Layer

Boundary-layer secondary circulations or ‘roll vortices’ can have a significant influence on the turbulent exchange of momentum, sensible heat and moisture throughout the hurricane boundary layer. In this study, analyses of data from a WP-3D aircraft of the National Oceanic and Atmospheric Administration (NOAA) are presented. As part of the Coupled Boundary Layer Air-Sea Transfer (CBLAST)-hurricane experiment sponsored through the Office of Naval Research and NOAA’s annual hurricane research program, flights were conducted to investigate energy exchange across the air–sea interface. We present the first in-situ aircraft-based observations of rolls in the hurricane boundary layer and investigate their influence on energy and momentum exchange. The rolls detected in Hurricane Isidore (year 2002) have a characteristic wavelength of about 900 m, in good agreement with analyses of data from a synthetic aperture radar image captured by the Canadian Space Agency’s RADARSAT satellite in the same storm. Our analyses of the airborne data suggest that roll vortices may be a significant factor modulating the air–sea momentum exchange.     

Weakened SST variability in the tropical Atlantic Ocean since 2000

A prominent weakening in equatorial Atlantic sea surface temperature (SST) variability, occurring around the year 2000, is investigated by means of observations, reanalysis products and the linear recharge oscillator (ReOsc) model. Compared to the time period 1982–1999, during 2000–2017 the May–June–July SST variability in the eastern equatorial Atlantic has decreased by more than 30%. Coupled air–sea feedbacks, namely the positive Bjerknes feedback and the negative net heat flux damping are important drivers for the equatorial Atlantic interannual SST variability. We find that the Bjerknes feedback weakened after 2000 while the net heat flux damping increased. The weakening of the Bjerknes feedback does not appear to be fully explainable by changes in the mean state of the tropical Atlantic. The increased net heat flux damping is related to an enhanced response of the latent heat flux to the SST anomalies (SSTa). Strengthened trade winds as well as warmer SSTs are suggested to increase the air–sea specific humidity difference and hence, enhancing the latent heat flux response to SSTa. A combined effect of those two processes is proposed to be responsible for the weakened SST variability in the eastern equatorial Atlantic. The ReOsc model supports the link between reduced SST variability, weaker Bjerknes feedback and stronger net heat flux damping.     

Air–sea coupling in the North Atlantic during summer

In this work, we have investigated the evolution of the summer air–sea interaction in the North Atlantic Ocean and the physical processes involved using reanalysis data and model simulation. It is found that an atmosphere disturbance over the North Atlantic Ocean in the preceding winter favors the build-up of a North Atlantic horseshoe-like sea surface temperature anomaly (SSTA) pattern in the summer through modifying the northeast trade winds and changing ocean upwelling and downwelling. The changed ocean condition (SSTA, upwelling, and downwelling) further intensifies the atmosphere disturbance as a positive feedback. The thermal advection of the atmosphere disturbance weakens the SSTA pattern in the following autumn and winter. The anomalous circulation associated with the air–sea interaction in the observations is characterized by a barotropic structure in the middle and high latitudes of the North Atlantic Ocean. The baroclinic component is enhanced in the model simulation, particularly in the seasons from summer to winter. The life cycle of the air–sea interaction is about 1 year in both the observations and simulations.     

Intercomparison of Air–Sea Interface Fluxes over the Yellow Sea and Korea Strait: Impact of Tsushima Warm Current

Surface-layer meteorological observations obtained from oceanic buoys over the Korean Strait and the Yellow Sea are used to estimate surface-layer turbulent fluxes of heat, moisture and momentum over the East-Asian Marginal Seas. Special emphasis is paid towards explanation of the impact of the Tsushima warm current flowing through the Korean Strait on air–sea interface fluxes. During the active phase of the Tsushima warm current, when the difference in sea surface temperature and air temperature becomes as large as 8°C, the sensible heat flux increases to a value of about 135 W m⁻², while the latent heat flux is around 200 W m⁻². The study attempts to broaden our understanding on the air-sea interaction processes over the Yellow Sea and Korean Strait.     

Aircraft Observations of Sea-Surface Turbulent Fluxes Near the California Coast

Aircraft turbulence data from the Autonomous Ocean Sampling Network project were analyzed and compared to the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk parametrization of turbulent fluxes in an ocean area near the coast of California characterized by complex atmospheric flow. Turbulent fluxes measured at about 35 m above the sea surface using the eddy-correlation method were lower than bulk estimates under unstable and stable atmospheric stratification for all but light winds. Neutral turbulent transfer coefficients were used in this comparison because they remove the effects of mean atmospheric conditions and atmospheric stability. Spectral analysis suggested that kilometre-scale longitudinal rolls affect significantly turbulence measurements even near the sea surface, depending on sampling direction. Cross-wind sampling tended to capture all the available turbulent energy. Vertical soundings showed low boundary-layer depths and high flux divergence near the sea surface in the case of sensible heat flux but minimal flux divergence for the momentum flux. Cross-wind sampling and flux divergence were found to explain most of the observed discrepancies between the measured and bulk flux estimates. At low wind speeds the drag coefficient determined with eddy correlation and an inertial dissipation method after corrections were applied still showed high values compared to bulk estimates. This discrepancy correlated with the dominance of sea swell, which was a usually observed condition under low wind speeds. Under stable atmospheric conditions measured sensible heat fluxes, which usually have low values over the ocean, were possibly affected by measurement errors and deviated significantly from bulk estimates.     

A Lagrangian Model to Predict the Modification of Near-Surface Scalar Mixing Ratios and Air–Water Exchange Fluxes in Offshore Flow

A model was developed to predict the modification with fetch in offshore flow of mixing ratio, air–water exchange flux, and near-surface vertical gradients in mixing ratio of a scalar due to air–water exchange. The model was developed for planning and interpretation of air–water exchange flux measurements in the coastal zone. The Lagrangian model applies a mass balance over the internal boundary layer (IBL) using the integral depth scale approach, previously applied to development of the nocturnal boundary layer overland. Surface fluxes and vertical profiles in the surface layer were calculated using the NOAA COARE bulk algorithm and gas transfer model (e.g., Blomquist et al. 2006, Geophys Res Lett 33:1–4). IBL height was assumed proportional to the square root of fetch, and estimates of the IBL growth rate coefficient, α, were obtained by three methods: (1) calibration of the model to a large dataset of air temperature and humidity modification over Lake Ontario in 1973, (2) atmospheric soundings from the 2004 New England Air Quality Study and (3) solution of a simplified diffusion equation and an estimate of eddy diffusivity from Monin–Obukhov similarity theory (MOST). Reasonable agreement was obtained between the calibrated and MOST values of α for stable, neutral, and unstable conditions, and estimates of α agreed with previously published parametrizations that were valid for the stable IBL only. The parametrization of α provides estimates of IBL height, and the model estimates modification of scalar mixing ratio, fluxes, and near-surface gradients, under conditions of coastal offshore flow (0–50 km) over a wide range in stability.     

Effects of Bubbles and Sea Spray on Air–Sea Exchange in Hurricane Conditions

The lower limit on the drag coefficient under hurricane force winds is determined by the break-up of the air–sea interface due to Kelvin–Helmholtz instability and formation of the two-phase transition layer consisting of sea spray and air bubbles. As a consequence, a regime of marginal stability develops. In this regime, the air–sea drag coefficient is determined by the turbulence characteristics of the two-phase transition layer. The upper limit on the drag coefficient is determined by the Charnock-type wave resistance. Most of the observational estimates of the drag coefficient obtained in hurricane conditions and in laboratory experiments appear to lie between the two extreme regimes: wave resistance and marginal stability.     

Impacts of atmosphere–sea ice–ocean interaction on Southern Ocean deep convection in a climate system model

Deep convection in polar oceans plays a critical role in the variability of global climate. In this study, we investigate potential impacts of atmosphere–sea ice–ocean interaction on deep convection in the Southern Ocean (SO) of a climate system model (CSM) by changing sea ice–ocean stress. Sea ice–ocean stress plays a vital role in the horizontal momentum exchange between sea ice and the ocean, and can be parameterized as a function of the turning angle between sea ice and ocean velocity. Observations have shown that the turning angle is closely linked to the sea-ice intrinsic properties, including speed and roughness, and it varies spatially. However, a fixed turning angle, i.e., zero turning angle, is prescribed in most of the state-of-the-art CSMs. Thus, sensitivities of SO deep convection to zero and non-zero turning angles are discussed in this study. We show that the use of a non-zero turning angle weakens open–ocean deep convection and intensifies continental shelf slope convection. Our analyses reveal that a non-zero turning angle first induces offshore movement of sea ice transporting to the open SO, which leads to sea ice decrease in the SO coastal region and increase in the open SO. In the SO coastal region, the enhanced sea-ice divergence intensifies the formation of denser surface water descending along continental shelf by enhanced salt flux and reduced freshwater flux, combined with enhanced Ekman pumping and weakened stratification, contributing to the occurrence and intensification of continental shelf slope convection. On the other hand, the increased sea ice in the open SO weakens the westerlies, enhances sea-level pressure, and increases freshwater flux, whilst oceanic cyclonic circulation slows down, sea surface temperature and sea surface salinity decrease in the open SO response to the atmospheric changes. Thus, weakened cyclonic circulation, along with enhanced freshwater flux, reduced deep–ocean heat content, and increased stability of sea water, dampens the open–ocean deep convection in the SO, which in turn cools the sea surface temperature, increases sea-level pressure, and finally increases sea-ice concentration, providing a positive feedback. In the CSM, the use of a non-zero turning angle has the capability to reduce the SO warm bias. These results highlight the importance of an accurate representation of sea ice–ocean coupling processes in a CSM.     

The effect of a regional increase in ocean surface roughness on the tropospheric circulation: a GCM experiment

The sensitivity of the atmospheric circulation to an increase in ocean surface roughness in the Southern Hemisphere storm track is investigated in a paired general circulation model experiment. Such a change in sea roughness could be induced by ocean waves generated by storms. Two extended permanent-July runs are made. One with standard sea surface roughness, the other with ten times as a large surface roughness over open sea poleward of 40° S. The regional increase in ocean surface roughness significantly modifies the tropospheric circulation in the Southern Hemisphere. The strongest effect is the reduction of tropospheric winds (by 2 m/s or 10%) above the area with increased roughness. The poleward eddy momentum flux is reduced in the upper troposphere and the meridional eddy sensible heat flux is reduced in the lower troposphere. Zonal mean and eddy kinetic energy are consistently reduced.     

Impact of diurnal atmosphere–ocean coupling on tropical climate simulations using a coupled GCM

The impacts of diurnal atmosphere–ocean (air–sea) coupling on tropical climate simulations are investigated using the SNU coupled GCM. To investigate the effect of the atmospheric and oceanic diurnal cycles on a climate simulation, a 1-day air–sea coupling interval experiment is compared to a 2-h coupling experiment. As previous studies have suggested, cold temperature biases over equatorial western Pacific regions are significantly reduced when diurnal air–sea coupling strategy is implemented. This warming is initiated by diurnal rectification and amplified further by the air–sea coupled feedbacks. In addition to its effect on the mean climatology, the diurnal coupling has also a distinctive impact on the amplitude of the El Nino-Southern Oscillation (ENSO). It is demonstrated that a weakening of the ENSO magnitude is caused by reduced (increased) surface net heat fluxes into the ocean during El Nino (La Nina) events. Primarily, decreased (increased) incoming shortwave radiation during El Nino (La Nina) due to cloud shading is responsible for the net heat fluxes associated with ENSO.     

Regimes of seasonal air–sea interaction and implications for performance of forced simulations

Sea surface temperature (SST) anomalies can induce anomalous convection through surface evaporation and low-level moisture convergence. This SST forcing of the atmosphere is indicated in a positive local rainfall–SST correlation. Anomalous convection can feedback on SST through cloud-radiation and wind-evaporation effects and wind-induced oceanic mixing and upwelling. These atmospheric feedbacks are reflected in a negative local rainfall–SST tendency correlation. As such, the simultaneous rainfall–SST and rainfall–SST tendency correlations can indicate the nature of local air–sea interactions. Based on the magnitude of simultaneous rainfall–SST and rainfall–SST tendency correlations, the present study identifies three distinct regimes of local air–sea interactions. The relative importance of SST forcing and atmospheric forcing differs in these regimes. In the equatorial central-eastern Pacific and, to a smaller degree, in the western equatorial Indian Ocean, SST forcing dominates throughout the year and the surface heat flux acts mainly as a damping term. In the tropical Indo-western Pacific Ocean regions, SST forcing and atmospheric forcing dominate alternatively in different seasons. Atmospheric forcing dominates in the local warm/rainy season. SST forcing dominates with a positive wind-evaporation feedback during the transition to the cold/dry season. SST forcing also dominates during the transition to the warm/rainy season but with a negative cloud-radiation feedback. The performance of atmospheric general circulation model simulations forced by observed SST is closely linked to the regime of air–sea interaction. The forced simulations have good performance when SST forcing dominates. The performance is low or poor when atmospheric forcing dominates.     

Boreal summer intraseasonal oscillations and seasonal Indian monsoon prediction in DEMETER coupled models

Even though multi-model prediction systems may have better skill in predicting the interannual variability (IAV) of Indian summer monsoon (ISM), the overall performance of the system is limited by the skill of individual models (single model ensembles). The DEMETER project aimed at seasonal-to-interannual prediction is not an exception to this case. The reasons for the poor skill of the DEMETER individual models in predicting the IAV of monsoon is examined in the context of the influence of external and internal components and the interaction between intraseasonal variability (ISV) and IAV. Recently it has been shown that the ISV influences the IAV through very long breaks (VLBs; breaks with duration of more than 10 days) by generating droughts. Further, all VLBs are associated with an eastward propagating Madden–Julian Oscillation (MJO) in the equatorial region, facilitated by air–sea interaction on intraseasonal timescales. This VLB-drought–MJO relationship is analyzed here in detail in the DEMETER models. Analyses indicate that the VLB-drought relationship is poorly captured by almost all the models. VLBs in observations are generated through air–sea interaction on intraseasonal time scale and the models’ inability to simulate VLB-drought relationship is shown to be linked to the models’ inability to represent the air–sea interaction on intraseasonal time scale. Identification of this particular deficiency of the models provides a direction for improvement of the model for monsoon prediction.     

寒潮过程中风浪对黄海海气热量通量和动量通量影响研究

采用2009—2013年CFSR(Climate Forecast System Reanalysis)大气和海洋再分析资料对黄海海气间热量通量和动量通量的特征进行统计分析,并通过FVCOMSWAVE浪流耦合模式对典型寒潮过程中风浪的影响效果进行模拟研究与对比分析。统计结果显示,通量受海表大风、海气温差及海洋环流等因子影响,秋冬季节强烈,春夏季节相对较弱,在寒潮活跃的冷季该海域的海流处于弱流期,风浪对海面通量的作用明显增强。海温特征也显示冷季的不稳定性显著强于暖季,因此该海域冷季具有更强的海气热量通量。沿岸站点的比较显示,南部吕泗站面向更开阔的东海海域,其平均波高高出北部20%左右。这与沿海南部通量强于北部特征对应。数值模拟显示,在寒潮过程中,海气界面热量通量和动量通量输送比多年月平均状态显著增强,动量通量增大1~5倍,热量通量增大1~6倍。寒潮过程入海冷锋走向、强度、移动方向显著影响海面热量通量和动量通量大值区的分布。偏北路寒潮纬向型冷锋入海,其强度东部大于西部,造成通量大值区形成在黄海东北部,而偏西路寒潮经向型冷锋入海,其强度南部大于北部,造成通量大值区形成在黄海南部。同时偏北路径寒潮强度大于偏西路径,海气动量通量响应较偏西路径强约25%,热量通量强约50%。耦合风浪作用的模拟显示,海气间热量通量和动量通量明显增大,对不同强度风浪,浪高增加1.5倍,动量通量最大值增大约2倍,热量通量增大10~160 W/m2;浪高减弱至0.5倍,动量通量最大值则减弱约40%,热量通量减小10~55 W/m2。冷锋及其驱动的风浪强烈影响区域海气通量时空特征。     

Local and synoptic mechanisms causing Southern California’s Santa Ana winds

The atmospheric conditions that lead to strong offshore surface winds in Southern California, commonly referred to as Santa Ana winds, are investigated using the North American Regional Reanalysis and a 12-year, 6-km resolution regional climate simulation of Southern California. We first construct an index to characterize Santa Ana events based on offshore wind strength. This index is then used to identify the average synoptic conditions associated with Santa Ana events—a high pressure anomaly over the Great Basin. This pressure anomaly causes offshore geostrophic winds roughly perpendicular to the region’s mountain ranges, which in turn cause surface flow as the offshore momentum is transferred to the surface. We find, however, that there are large variations in the synoptic conditions during Santa Ana conditions, and that there are many days with strong offshore flow and weak synoptic forcing. This is due to local thermodynamic forcing that also causes strong offshore surface flow: a large temperature gradient between the cold desert surface and the warm ocean air at the same altitude creates an offshore pressure gradient at that altitude, in turn causing katabatic-like offshore flow in a thin layer near the surface. We quantify the contribution of “synoptic” and “local thermodynamic” mechanisms using a bivariate linear regression model, and find that, unless synoptic conditions force strongly onshore winds, the local thermodynamic forcing is the primary control on Santa Ana variability.     

Aerodynamic roughness of the sea surface at high winds

The role of the surface roughness in the formation of the aerodynamic friction of the water surface at high wind speeds is investigated. The study is based on a wind-over-waves coupling theory. In this theory waves provide the surface friction velocity through the form drag, while the energy input from the wind to waves depends on the friction velocity and the wind speed. The wind-over-waves coupling model is extended to high wind speeds taking into account the effect of sheltering of the short wind waves by the air-flow separation from breaking crests of longer waves. It is suggested that the momentum and energy flux from the wind to short waves locally vanishes if they are trapped into the separation bubble of breaking longer waves. At short fetches, typical for laboratory conditions, and strong winds the steep dominant wind waves break frequently and provide the major part of the total form drag through the air-flow separation from breaking crests, and the effect of short waves on the sea drag is suppressed. In this case the dependence of the drag coefficient on the wind speed is much weaker than would be expected from the standard parameterization of the roughness parameter through the Charnock relation. At long fetches, typical for the field, waves in the spectral peak break rarely and their contribution to the air-flow separation is weak. In this case the surface form drag is determined predominantly by the air-flow separation from breaking of the equilibrium range waves. As found at high wind speeds up to 60 m s⁻¹ the modelled aerodynamic roughness is consistent with the Charnock relation, i.e. there is no saturation of the sea drag. Unlike the aerodynamic roughness, the geometrical surface roughness (height of short waves) could be saturated or even suppressed when the wind speed exceeds 30 m s⁻¹.     

On the Possible Impact of a Following-Swell on the Atmospheric Boundary Layer

A simple model of the atmospheric boundary layer over the ocean where the swell impact on the atmosphere is explicitly accounted for is suggested. The model is based on Ekman’s equations, where the stress in the wave boundary layer is split into two parts: the turbulent and wave-induced stress. The turbulent stress is parameterized traditionally via the eddy viscosity proportional to the generalized mixing length. The wave-induced stress directed upward (from swell to the atmosphere) is parameterized using the formalism of the wind-over-waves coupling theory. The model can be seen as an extension of the model by Kudryavtsev and Makin (J Phys Oceanogr 34:934–949, 2004) to the scale of the entire atmospheric boundary layer by including the Coriolis force into the momentum conservation equation and generalizing the definition of the mixing length. The regime of low winds for swell propagating along the wind direction is studied. It is shown that the impact of swell on the atmosphere is governed mainly by the swell parameter—the coupling parameter that is the product of the swell steepness and the growth rate coefficient. When the coupling parameter drops below − 1 the impact of swell becomes significant and affects the entire atmospheric boundary layer. The turbulent stress is enhanced near the surface as compared to the no-swell case, and becomes negative above the height of the inner region. The wind profile is characterized by a positive gradient near the surface and a negative gradient above the height of the inner region forming a characteristic bump at the height of the inner region. Results of the model agree at least qualitatively with observations performed in the atmosphere in presence of swell.     

Simulation of the Atlantic meridional overturning circulation in an atmosphere–ocean global coupled model. Part I: a mechanism governing the variability of ocean convection in a preindustrial experiment

A preindustrial climate experiment was conducted with the third version of the CNRM global atmosphere–ocean–sea ice coupled model (CNRM-CM3) for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). This experiment is used to investigate the main physical processes involved in the variability of the North Atlantic ocean convection and the induced variability of the Atlantic meridional overturning circulation (MOC). Three ocean convection sites are simulated, in the Labrador, Irminger and Greenland–Iceland–Norwegian (GIN) Seas in agreement with observations. A mechanism linking the variability of the Arctic sea ice cover and convection in the GIN Seas is highlighted. Contrary to previous suggested mechanisms, in CNRM-CM3 the latter is not modulated by the variability of freshwater export through Fram Strait. Instead, the variability of convection is mainly driven by the variability of the sea ice edge position in the Greenland Sea. In this area, the surface freshwater balance is dominated by the freshwater input due to the melting of sea ice. The ice edge position is modulated either by northwestward geostrophic current anomalies or by an intensification of northerly winds. In the model, stronger than average northerly winds force simultaneous intense convective events in the Irminger and GIN Seas. Convection interacts with the thermohaline circulation on timescales of 5–10 years, which translates into MOC anomalies propagating southward from the convection sites.     

‘Modelling the Arctic Boundary Layer: An Evaluation of Six Arcmip Regional-Scale Models using Data from the Sheba Project’

A primary climate change signal in the central Arctic is the melting of sea ice. This is dependent on the interplay between the atmosphere and the sea ice, which is critically dependent on the exchange of momentum, heat and moisture at the surface. In assessing the realism of climate change scenarios it is vital to know the quality by which these exchanges are modelled in climate simulations. Six state-of-the-art regional-climate models are run for one year in the western Arctic, on a common domain that encompasses the Surface Heat Budget of the Arctic Ocean (SHEBA) experiment ice-drift track. Surface variables, surface fluxes and the vertical structure of the lower troposphere are evaluated using data from the SHEBA experiment. All the models are driven by the same lateral boundary conditions, sea-ice fraction and sea and sea-ice surface temperatures. Surface pressure, near-surface air temperature, specific humidity and wind speed agree well with observations, with a falling degree of accuracy in that order. Wind speeds have systematic biases in some models, by as much as a few metres per second. The surface radiation fluxes are also surprisingly accurate, given the complexity of the problem. The turbulent momentum flux is acceptable, on average, in most models, but the turbulent heat fluxes are, however, mostly unreliable. Their correlation with observed fluxes is, in principle, insignificant, and they accumulate over a year to values an order of magnitude larger than observed. Typical instantaneous errors are easily of the same order of magnitude as the observed net atmospheric heat flux. In the light of the sensitivity of the atmosphere–ice interaction to errors in these fluxes, the ice-melt in climate change scenarios must be viewed with considerable caution.     

Impact of Sea-Spray on the Atmospheric Surface Layer

The feedback effects of sea-spray on the heat and momentum fluxes under equilibrium conditions associated with winds of tropical cyclones are investigated using a one-dimensional coupled sea-spray and atmospheric surface-layer (ASL) model. This model is capable of simulating the microphysical aspects of the evaporation of saline water droplets of various sizes and their dynamic and thermal interaction with the turbulence mixing that is simulated by the Mellor–Yamada 1.5-order closure scheme. Sea-spray droplet generation is described by a state-of-the-art parametrization that predicts the size spectrum of sea-spray droplets for a given surface forcing. The results from a series of simulations indicate the way in which evaporating droplets of various sizes modify the turbulence mixing near the surface, which in turn affects further droplet evaporation. All these results are direct consequences of the effects of sea-spray on the balance of turbulent kinetic energy in the spray-filled surface layer. In particular, the overall impact of sea-spray droplets on the mean wind depends on the wind speed at the level of sea-spray generation. When the wind speed is below 40 m s⁻¹, the droplets are small in size and tend to evaporate substantially and thus cool the spray-filled layer, while for wind speeds above 50 m s⁻¹, the size of the droplets is so large that they do not have enough time to evaporate much before falling back into the sea. The sensible heat carried by the droplets is released to the ambient air, increasing the buoyancy of the surface layer and enhancing the turbulent mixing. The suspension of sea-spray droplets reduces the buoyancy and makes the surface layer more stable, decreasing the friction velocity and the downward turbulent mixing of momentum. The results from the numerical experiments also suggest that, in order not to violate the constant flux assumption critical to the Monin–Obukhov similarity theory, a displacement equal to the mean wave height should be included in the logarithmic profiles of the wind and thermal fields.