Using observations from the Martha's Vineyard Coastal Observatory, this paper investigates how momentum flux in the marine atmospheric surface layer over the coastal ocean varies with fetch, wave age, and wave slope and assesses the performance of the COARE 3.5 bulk parameterization. Long-fetch (at least 300 km) and short-fetch (3-6 km away from land) conditions have very similar momentum flux, with the latter being just 15% higher. The COARE 3.5 wind speed-dependent formulation closely matches the observations. The sea state dependence of wave age and wave slope is analyzed by considering both peak frequency and mean frequency in a wave spectrum. The slope of regression lines between normalized roughness and wave age is sensitive to wind speed ranges and the scatter of momentum flux, potentially explaining why earlier studies did not find a universal formula to characterize the wave dependence. Although the observed momentum flux exhibits an obvious dependence on wave age, a robust quadratic fit between momentum flux and the 10-m neutral wind speed exists only for young waves. For the short-fetch conditions, the momentum flux does not increase as wave slope increases. This may be because waves generated by local wind are still weak and swell that propagates from other areas dominates the wave spectra. In other words, the wave slope computed by considering the spectra does not properly reflect the local wind-wave interaction.

Large-scale offshore wind farms are expected to influence surface waves by modifying local wind forcing through wake effects. We use regional coupled ocean-atmosphere-wave model simulations to investigate a realistic large-scale offshore wind development scenario in the northeastern U.S. during boreal summer. Near-surface wind speeds are reduced by 10% over lease areas and within downstream wake regions, leading to decreases in significant wave height (3%) and wave-supported momentum flux (30%). This further leads to reductions in surface roughness length (16%) and near-surface ocean turbulent kinetic energy (20%). Spectral analysis shows a clear reduction in wind-sea energy, indicating suppressed local wind-wave growth near the wind farms. Weaker winds favor the development of longer-period waves, increasing dominant wave phase speed by 3% and suggesting a transition to an older sea state. Modern bulk flux algorithms often parameterize surface roughness using inverse wave age and/or wave slope. This raises the question of whether wake-driven reductions in inverse wave age and wave height impact air-sea momentum exchange. To assess this, we compare fully coupled simulations with an atmosphere-only run excluding wave coupling. Results show that about one-third of the reduction in roughness length can be attributed to sea state changes, while two-thirds result from lower friction velocity due to lower wind speeds. However, the impact of sea state on the drag coefficient and momentum flux is negligible (1%), suggesting that wake-induced wind speed reductions are the primary driver, with sea state changes playing a secondary role.

Energy is transferred from the atmosphere to the ocean primarily through ocean surface waves, and the majority is dissipated locally in the near‐surface ocean. Observations of turbulent kinetic energy (TKE) in the upper ocean have shown dissipation rates exceeding law‐of‐the‐wall theory by an order of magnitude. The excess near‐surface ocean TKE dissipation rate is thought to be driven primarily by wave breaking, which limits wave growth and transfers energy from the surface wave field to the wave‐affected layer of the ocean. Here, the statistical properties of breaking wave dynamics in a coastal area are extracted from visible imagery and used to estimate TKE dissipation rates due to breaking waves. The statistical properties of whitecap dynamics are quantified with Λ(c), a distribution of total whitecap crest length per unit area as a function of crest speed, and used to compute energy dissipation by breaking waves, S ds . S ds approximately balances elevated subsurface dissipation in young seas but accounts for only a fraction of subsurface dissipation in older seas. The wind energy input is estimated from wave spectra from polarimetric imagery and laser altimetry. S ds balances the wind energy input except under high winds. Λ(c)‐derived estimates of TKE dissipation rates by breaking waves compare well with the atmospheric deficit in TKE dissipation, a measure of energy input to the wave field (Cifuentes‐Lorenzen et al., 2024). These results tie the observed atmospheric dissipation deficit and enhancement in subsurface TKE dissipation to wave driven energy transport, constraining the TKE dissipation budget near the air‐sea interface.

In this work, we seek to address the validity of Monin-Obukhov similarity theory (MOST) in the wave-affected surface boundary layer of the atmosphere. While bulk flux formulas which rely on MOST have been tested with and applied to measurements and models of air/sea interaction for several decades, the influences of surface wave-mediated fluxes on MOST have not been thoroughly quantified. We assess several months of direct covariance data from a stationary tower deployed with instruments inside the wave-affected surface layer. These measurements are analyzed in the context of the turbulent kinetic energy (TKE) equation and MOST, extending previous work due to the inclusion of directly estimated wave-coherent energy fluxes. Scaled TKE dissipation rates are reduced from what is predicted by MOST during events with large wave-coherent surface fluxes, resulting in a dissipation deficit in the energy budget (roughly 30%). However, we find that shear is much less impacted by these wave events showing much smaller deviations from baselines (less than 10%). During much of the experiment, the dissipation deficit is balanced by the wave-coherent pressure work, suggesting a general understanding of the combined turbulent and wave-driven energetics. However, several large storms in the fall of 2022 yielded much larger dissipation deficits than can be explained by the wave-coherent pressure work, highlighting that more work is needed to understand energetics in the wave-affected surface layer more generally.

Hurricane Epsilon (2020) was a late-season, category-3 tropical cyclone that underwent extratropical transition and became Extratropical Cyclone Epsilon on October 26. The upper ocean response to the passage of the storm was observed by three types of autonomous platforms: an eXpendable Spar buoy, an Air-Launched Autonomous Micro-Observer profiling float, and two Seagliders. Taken together, this array enabled the rare collection of contemporaneous observations of the upper ocean, air-sea interface, and atmospheric boundary layer before, during, and after the passage of the storm. The evidence presented suggests that Extratropical Cyclone Epsilon contributed to breaking down the residual North Atlantic summer stratification regime and accelerated the shift to the prolonged ocean cooling associated with winter. The synergistic capabilities of the observational array are significant for two reasons: (1) by enabling the comparison of complementary atmosphere and ocean observations, taken from different platforms, they permit a comprehensive approach to better understand how storm-induced momentum, heat, and moisture fluxes alter upper ocean structure, and (2) they demonstrate the ability of future, targeted deployments of similar observational arrays to assess the fidelity of coupled ocean-atmosphere-wave numerical prediction models.

Ocean surface radiation measurement best practices have been developed as a first step to support the interoperability of radiation measurements across multiple ocean platforms and between land and ocean networks. This document describes the consensus by a working group of radiation measurement experts from land, ocean, and aircraft communities. The scope was limited to broadband shortwave (solar) and longwave (terrestrial infrared) surface irradiance measurements for quantification of the surface radiation budget. Best practices for spectral measurements for biological purposes like photosynthetically active radiation and ocean color are only mentioned briefly to motivate future interactions between the physical surface flux and biological radiation measurement communities. Topics discussed in these best practices include instrument selection, handling of sensors and installation, data quality monitoring, data processing, and calibration. It is recognized that platform and resource limitations may prohibit incorporating all best practices into all measurements and that spatial coverage is also an important motivator for expanding current networks. Thus, one of the key recommendations is to perform interoperability experiments that can help quantify the uncertainty of different practices and lay the groundwork for a multi-tiered global network with a mix of high-accuracy reference stations and lower-cost platforms and practices that can fill in spatial gaps.

Surface waves grow through a mechanism in which atmospheric pressure is offset in phase from the wavy surface. A pattern of low atmospheric pressure over upward wave orbital motions and high pressure over downward wave orbital motions travels with the water wave, leading to a pumping of kinetic energy from the atmospheric boundary layer into the waves. This pressure pattern persists above the air/water interface, modifying the turbulent kinetic energy in the atmospheric wave-affected boundary layer. Here, we present field measurements of the transfer of energy from wind to waves through wave-coherent atmospheric pressure work. Measured pressure work cospectra are consistent with an existing model for atmospheric pressure work. Measured pressure work energy fluxes reach 0.1-0.2 W m$^{-2}$ during the largest measured wind event (winds reaching 16.5 m s$^{-1}$). The implications for these measurements and their importance to the turbulent kinetic energy budget are discussed.

This work serves as an observation-based exploration into the role of wave-driven turbulence at the air-sea interface by measuring Turbulent Kinetic Energy (TKE) dissipation rates above and below the sea surface. Subsurface ocean measurements confirm a TKE dissipation rate enhancement relative to the predicted law-of-the-wall (epsilon(obs) > epsilon(p)), which appears to be fully supported by wave breaking highlighting the role of the transport terms in balancing the subsurface TKE budget. Simultaneous measurements of TKE dissipation rates on the atmospheric side capture a deficit relative to the law-of-the-wall (epsilon(obs) < epsilon(p)). This deficit is explained in terms of wave-induced perturbations, with observed convergence to the law-of-the-wall at 14 m above mean sea level. The deficit on the atmospheric side provides an estimate of the energy flux divergence in the wave boundary layer. An exponential function is used to integrate in the vertical and provide novel estimates of the amount of energy going into the wave field. These estimates correlate well with classic spectral input parameterizations and can be used to derive an effective wave-scale, capturing wind-wave coupling purely from atmospheric observations intimately tied to wave-induced perturbations of the air-flow. These atmospheric and oceanic observations corroborate the commonly assumed input-dissipation balance for waves at wind speeds in the 8-14 ms(-1) range in the presence of developed to young seas. At wind speeds above 14 ms(-1) under young seas (U10cp>1.2)observations suggest a deviation from the TKE input-dissipation balance in the wave field.

Atmospheric fronts embedded in extratropical cyclones are high-impact weather phenomena, contributing significantly to midlatitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves on air-sea fluxes remain undocumented. Using the multi-year in situ near-surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hours. The fully-coupled model simulations indicate that the misaligned waves significantly increase wave roughness length (300%), drag coefficient (30%), and momentum flux (20%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the compensating effect between the decrease in temperature and humidity scaling parameters and the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave-based bulk flux algorithm. Yet, the suggested modification to the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind-wave coupling in the bulk formula of the numerical models may help improve the air-sea interaction simulations under the passing atmospheric fronts in the midlatitudes.

In winter, the Northwest Tropical Atlantic Ocean can be characterized by various wave age-based interactions among ocean current, surface wind and surface waves, which are critical for accurately describing surface wind stress. In this work, coupled wave-ocean-atmosphere model simulations are conducted using two different wave roughness parameterizations within COARE3.5, including one that relies solely on wind speed and another that uses wave age and wave slope as inputs. Comparisons with the directly measured momentum fluxes during the ATOMIC/EUREC(4)A experiments in winter 2020 show that, for sea states dominated by short wind waves under moderate to strong winds, the wave-based formulation (WBF) increases the surface roughness length in average by 25% compared to the wind-speed-based approach. For sea states dominated by remotely generated swells under moderate to strong wind intensity, the WBF predicts significantly lower roughness length and surface stress (asymptotic to 15%), resulting in increased near-surface wind speed above the constant flux layer (asymptotic to 5%). Further investigation of the mixed sea states in the model and data indicates that the impact of swell on wind stress is over-emphasized in the COARE3.5 WBF, especially under moderate wind regimes. Various approaches are explored to alleviate this deficiency by either introducing directional alignment between wind and waves or using the mean wave period instead of the wave period corresponding to the spectral peak to compute the wave age. The findings of this study are likely to be site-dependent, and mostly concern specific regimes of wind and waves where the original parameterization was deficient.

High-frequency wind measurements from Saildrone autonomous surface vehicles are used to calculate wind stress in the tropical east Pacific. Comparison between direct covariance (DC) and bulk wind stress estimates demonstrates very good agreement. Building on previous work that showed the bulk input data were reliable, our results lend credibility to the DC estimates. Wind flow distortion by Saildrones is comparable to or smaller than other platforms. Motion correction re-sults in realistic wind spectra, albeit with signatures of swell-coherent wind fluctuations that may be unrealistically strong. Fractional differences between DC and bulk wind stress magnitude are largest at wind speeds below 4 m s21. The size of this effect, however, depends on choice of stress direction assumptions. Past work has shown the importance of using current -relative (instead of Earth-relative) winds to achieve accurate wind stress magnitude. We show that it is also important for wind stress direction.

Bubble plume penetration depths have been identified as a key parameter linking subsurface turbulent kinetic energy (TKE) dissipation rates and whitecaps. From data collected in the Atlantic sector of the Southern Ocean, nominally 50 degrees S 40 degrees W, bubble plume penetration depths were estimated from Acoustic Doppler Current Profiler measurements of the acoustic backscatter anomaly. Bubble presence at depth was corroborated using independent measurements of optical scattering. Here, an effective wavelength, observations of significant wave height and atmospheric forcing were used to scale penetration depths of breaking waves under open ocean conditions. The parameterization was developed assuming a correlation between the observed penetration depth and an estimate of the TKE dissipation rate enhancement under breaking waves. The effective wavelength was defined from the effective phase speed based on a momentum and energy balance across the atmospheric wave boundary layer and was considered to be the largest actively wind-coupled wave and representative of large-scale breaking for wave ages ranging from 15 to 35 (i.e., 15 = = 35). This yields a dimensional penetration depth parameterization in terms of inverse wave age and the length scales under consideration. The parameterization captures the bubble plume penetration depth with stronger forcing leading to deeper injections, reaching up to 9 m. Both length scales are effective at defining the depth of a wave-affected layer in terms of bubble presence with the effective wavelength better collapsing the data under mixed conditions with deeper plumes associated to larger fractional whitecap coverage.

The Observing Air-Sea Interactions Strategy (OASIS) is a new United Nations Decade of Ocean Science for Sustainable Development programme working to develop a practical, integrated approach for observing air-sea interactions globally for improved Earth system (including ecosystem) forecasts, CO2 uptake assessments called for by the Paris Agreement, and invaluable surface ocean information for decision makers. Our "Theory of Change" relies upon leveraged multi-disciplinary activities, partnerships, and capacity strengthening. Recommendations from >40 OceanObs'19 community papers and a series of workshops have been consolidated into three interlinked Grand Ideas for creating #1: a globally distributed network of mobile air-sea observing platforms built around an expanded array of long-term time-series stations; #2: a satellite network, with high spatial and temporal resolution, optimized for measuring air-sea fluxes; and #3: improved representation of air-sea coupling in a hierarchy of Earth system models. OASIS activities are organized across five Theme Teams: (1) Observing Network Design & Model Improvement; (2) Partnership & Capacity Strengthening; (3) UN Decade OASIS Actions; (4) Best Practices & Interoperability Experiments; and (5) Findable-Accessible-Interoperable-Reusable (FAIR) models, data, and OASIS products. Stakeholders, including researchers, are actively recruited to participate in Theme Teams to help promote a predicted, safe, clean, healthy, resilient, and productive ocean.

Two decades of high-resolution satellite observations and climate modeling studies have indicated strong ocean–atmosphere coupled feedback mediated by ocean mesoscale processes, including semipermanent and meandrous SST fronts, mesoscale eddies, and filaments. The air–sea exchanges in latent heat, sensible heat, momentum, and carbon dioxide associated with this so-called mesoscale air–sea interaction are robust near the major western boundary currents, Southern Ocean fronts, and equatorial and coastal upwelling zones, but they are also ubiquitous over the global oceans wherever ocean mesoscale processes are active. Current theories, informed by rapidly advancing observational and modeling capabilities, have established the importance of mesoscale and frontal-scale air–sea interaction processes for understanding large-scale ocean circulation, biogeochemistry, and weather and climate variability. However, numerous challenges remain to accurately diagnose, observe, and simulate mesoscale air–sea interaction to quantify its impacts on large-scale processes. This article provides a comprehensive review of key aspects pertinent to mesoscale air–sea interaction, synthesizes current understanding with remaining gaps and uncertainties, and provides recommendations on theoretical, observational, and modeling strategies for future air–sea interaction research. Recent high-resolution satellite observations and climate models have shown a significant impact of coupled ocean–atmosphere interactions mediated by small-scale (mesoscale) ocean processes, including ocean eddies and fronts, on Earth’s climate. Ocean mesoscale-induced spatial temperature and current variability modulate the air–sea exchanges in heat, momentum, and mass (e.g., gases such as water vapor and carbon dioxide), altering coupled boundary layer processes. Studies suggest that skillful simulations and predictions of ocean circulation, biogeochemistry, and weather events and climate variability depend on accurate representation of the eddy-mediated air–sea interaction. However, numerous challenges remain in accurately diagnosing, observing, and simulating mesoscale air–sea interaction to quantify its large-scale impacts. This article synthesizes the latest understanding of mesoscale air–sea interaction, identifies remaining gaps and uncertainties, and provides recommendations on strategies for future ocean–weather–climate research.

Rainfall alters the physical and chemical properties of the surface ocean, and its effect on ocean skin temperature and surface heat fluxes is poorly represented in many air-sea interaction models. We present radiometric observations of ocean skin temperature, near-surface (5 cm) temperature from a towed thermistor, and bulk atmospheric and oceanic variables, for 69 rain events observed over the course of 4 months in the Indian Ocean as part of the DYNAMO project. We test a state-of-the-art prognostic model developed by Bellenger et al. (2017, https://doi.org/10.1002/2016JC012429) to predict ocean skin temperature in the presence of rain, and demonstrate a physically motivated modification to the model that improves its performance with increasing rain rate. We characterize the vertical skin-bulk temperature gradient induced by rain and find that it levels off at high rain rates, suggestive of a transition in skin-layer physics that has been previously hypothesized in the literature. We also quantify the small bias that will be present in turbulent sensible heat fluxes parameterized from ocean temperature measurements made at typical "bulk" depths during a rain event. Finally, a wind threshold is observed above which the surface ocean remains well-mixed during a rain event; however, the skin temperature is observed to decrease at all wind speeds in the presence of rain.

In winter, the Northwest Tropical Atlantic Ocean can be characterized by various regimes of interactions among ocean current, surface wind, and wind waves, which are critical for accurately describing surface wind stress. In this work, coupled wave-ocean-atmosphere model simulations are conducted using two different wave roughness parameterizations within COARE3.5, including one that relies solely on wind speed and another that uses wave age and wave slope as inputs. Comparisons with the directly measured momentum fluxes during the ATOMIC/EUREC4A experiments in winter 2020 show that, for sea states dominated by short wind waves under moderate to strong winds, the wave-based formulation increases the surface roughness length by 40\% compared to the wind-speed-based approach. For sea states dominated by remotely generated swells under moderate to strong wind intensity, the wave-based formulation predicts significantly lower roughness length and surface stress (~20%), resulting in increased near-surface wind speed above the constant flux layer (~5%). Further investigation of the mixed sea states in the model and data indicates that the impact of swell on wind stress is over-emphasized in the COARE3.5 wave-based formulation, especially under moderate wind regimes. Various approaches are explored to alleviate this deficiency by either introducing directional alignment between wind and waves or using the mean wave period instead of the wave period corresponding to the spectral peak to compute the wave age.

Physical understanding and modeling of Ka-band ocean surface backscatter is challenging due to a lack of measurements. In the framework of the NASA Earth Ventures Suborbital-3 Submesoscale Ocean Dynamics Experiment (S-MODE) mission, a Ka-Band Ocean continuous wave Doppler Scatterometer (KaBODS) built by the University of Massachusetts, Amherst (UMass) was installed on the Woods Hole Oceanographic Institution (WHOI) Air-Sea Interaction Tower. Together with ASIT anemometers, a new data set of Ka-band ocean surface backscatter measurements along with surface wind/wave and weather parameters was collected. In this work, we present the KaBODS instrument and an empirical Ka-band wind Geophysical Model Function (GMF), the so-called ASIT GMF, based on the KaBODS data collected over a period of three months, from October 2019 to January 2020, for incidence angles ranging between 40° and 68°. The ASIT GMF results are compared with an existing Ka-band wind GMF developed from data collected during a tower experiment conducted over the Black Sea. The two GMFs show differences in terms of wind speed and wind direction sensitivity. However, they are consistent in the values of the standard deviation of the model residuals. This suggests an intrinsic geophysical variability characterizing the Ka-band surface backscatter. The observed variability does not significantly change when filtering out swell-dominated data, indicating that the long-wave induced backscatter modulation is not the primary source of the KaBODS backscatter variability. We observe evidence of wave breaking events, which increase the skewness of the backscatter distribution in linear space, consistent with previous studies. Interestingly, a better agreement is seen between the GMFs and the actual data at an incidence angle of 60° for both GMFs, and the statistical analysis of the model residuals shows a reduced backscatter variability at this incidence angle. This study shows that the ASIT data set is a valuable reference for studies of Ka-band backscatter. Further investigations are on-going to fully characterize the observed variability and its implication in the wind GMF development.

The warm Gulf Stream sea surface temperatures strongly impact the evolution of winter clouds behind atmospheric cold fronts. Such cloud evolution remains challenging to model. The Gulf Stream is too wide within the ERA5 and MERRA2 reanalyses, affecting the turbulent surface fluxes. Known problems within the ERA5 boundary layer (too-dry and too-cool with too strong westerlies), ascertained primarily from ACTIVATE 2020 campaign aircraft dropsondes and secondarily from older buoy measurements, reinforce surface flux biases. In contrast, MERRA2 winter surface winds and air-sea temperature/humidity differences are slightly too weak, producing surface fluxes that are too low. Reanalyses boundary layer heights in the strongly forced winter cold-air-outbreak regime are realistic, whereas late-summer quiescent stable boundary layers are too shallow. Nevertheless, the reanalysis biases are small, and reanalyses adequately support their use for initializing higher-resolution cloud process modeling studies of cold-air outbreaks.