[1] Direct measurements of air-sea heat, momentum, and mass (including CO 2 , DMS, and water vapor) fluxes using the direct covariance method were made over the open ocean from the NOAA R/V Ronald H. Brown during the Southern Ocean Gas Exchange (SO GasEx) program. Observations of fluxes and the physical processes associated with driving air-sea exchange are key components of SO GasEx. This paper focuses on the exchange of CO 2 and the wind speed dependency of the transfer velocity, k, used to model the CO 2 flux between the atmosphere and ocean. A quadratic dependence of k on wind speed based on dual tracer experiments is most frequently encountered in the literature. However, in recent years, bubble-mediated enhancement of k, which exhibits a cubic relationship with wind speed, has emerged as a key issue for flux parameterization in high-wind regions. Therefore, a major question addressed in SO GasEx is whether the transfer velocities obey a quadratic or cubic relationship with wind speed. After significant correction to the flux estimates (primarily due to moisture contamination), the direct covariance CO 2 fluxes confirm a significant enhancement of the transfer velocity at high winds compared with previous quadratic formulations. Regression analysis suggests that a cubic relationship provides a more accurate parameterization over a wind speed range of 0 to 18 m s −1 . The Southern Ocean results are in good agreement with the 1998 GasEx experiment in the North Atlantic and a recent separate field program in the North Sea.

Concurrent wavefield and turbulent flux measurements acquired during the Southern Ocean (SO) Gas Exchange (GasEx) and the High Wind Speed Gas Exchange Study (HiWinGS) projects permit evaluation of the dependence of the whitecap coverage W on wind speed, wave age, wave steepness, mean square slope, and wind-wave and breaking Reynolds numbers. The W was determined from over 600 high-frequency visible imagery recordings of 20 min each. Wave statistics were computed from in situ and remotely sensed data as well as from a WAVEWATCH III hindcast. The first shipborne estimates of W under sustained 10-m neutral wind speeds U10N of 25 m s−1 were obtained during HiWinGS. These measurements suggest that W levels off at high wind speed, not exceeding 10% when averaged over 20 min. Combining wind speed and wave height in the form of the wind-wave Reynolds number resulted in closely agreeing models for both datasets, individually and combined. These are also in good agreement with two previous studies. When expressing W in terms of wavefield statistics only or wave age, larger scatter is observed and/or there is little agreement between SO GasEx, HiWinGS, and previously published data. The wind speed–only parameterizations deduced from the SO GasEx and HiWinGS datasets agree closely and capture more of the observed W variability than Reynolds number parameterizations. However, these wind speed–only models do not agree as well with previous studies than the wind-wave Reynolds numbers.

Predicting future climate hinges on our understanding of and ability to quantify air‐sea gas transfer. The latter relies on parameterizations of the gas transfer velocity k, which represents physical mass transfer mechanisms and is usually parameterized as a nonlinear function of wind forcing. In an attempt to reduce uncertainties in k, this study explores empirical parameterizations that incorporate both wind speed and sea state dependence via wave‐wind and breaking Reynolds numbers, RH and RB. Analysis of concurrent eddy covariance gas transfer and measured wavefield statistics supplemented by wave model hindcasts shows for the first time that wave‐related Reynolds numbers collapse four open ocean data sets that have a wind speed dependence of CO2 transfer velocity ranging from lower than quadratic to cubic. Wave‐related Reynolds number and wind speed show comparable performance for parametrizing dimethyl sulfide (DMS) which, because of its higher solubility, is less affected by bubble‐mediated exchange associated with wave breaking.

We investigate the momentum and energy exchange across the wave boundary layer (WBL). Directly at the air-sea interface, we test three wave-growth parametrizations by comparing estimates of the wave-induced momentum flux derived from wave spectra with direct covariance estimates of the momentum flux. An exponential decay is used to describe the vertical structure of the wave-induced momentum in the atmospheric WBL through use of a decay rate, a function of the dimensionless decay rate and wavenumber (A α k). The decay rate is varied to minimize the difference between the energy extracted from the WBL and the energy flux computed from wave spectra using our preferred wave-growth parametrization. For wave ages (i.e. the peak phase speed to atmospheric friction velocity ratio) in the range 15 < c p /u * < 35 we are able to balance these two estimates to within 10%. The decay rate is used to approximate the WBL height as the height to which the wave-induced flux is 0.1 of its surface value and the WBL height determined this way is found to be between 1-3 m. Finally, we define an effective phase speed with which to parametrize the energy flux for comparison with earlier work, which we ultimately attempt to parametrize as a function of wind forcing.

Bubble size distributions ranging from 0.5 to 125 lm radius were measured optically during high winds of 13 m s 21 and large-scale wave breaking as part of the Southern Ocean Gas Exchange Experiment.Very small bubbles with radii less than 60 mm were measured at 6-9 m depth using optical measurements of the near-forward volume scattering function and critical scattering angle for bubbles (80 ). The bubble size distributions generally followed a power law distribution with mean slope values ranging from 3.6 to 4.6. The steeper slopes measured here were consistent with what would be expected near the base of the bubble plume. Bubbles, likely stabilized with organic coatings, were present for time periods on the order of 10-100 s at depths of 6-9 m. Here, relatively young seas, with an inverse wave age of approximately 0.88 and shorter characteristic wave scales, produced lower bubble concentrations, shallower bubble penetration depths, and steep bubble size distribution slopes. Conversely, older seas, with an inverse wave age of 0.70 and longer characteristic wave scales, produced relatively higher bubble concentrations penetrating to 15 m depth, larger bubble sizes, and shallower bubble size distribution slopes. When extrapolated to 4 m depth using a previously published bubble size distribution, our estimates suggest that the deeply penetrating small bubbles measured in the Southern Ocean supplied 36% of the total void fraction and likely contributed to the transfer and supersaturation of low-solubility gases.