Journal cover Journal topic
Ocean Science An interactive open-access journal of the European Geosciences Union
Journal topic
Ocean Sci., 15, 235-247, 2019
https://doi.org/10.5194/os-15-235-2019
Ocean Sci., 15, 235-247, 2019
https://doi.org/10.5194/os-15-235-2019

Research article 08 Mar 2019

Research article | 08 Mar 2019

# Measurements of air–sea gas transfer velocities in the Baltic Sea

Gas transfer velocities in the Baltic Sea
Leila Nagela, Kerstin E. Krall1, and Bernd Jähne1,2 Leila Nagel et al.
• 1Institute of Environmental Physics, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
• 2Heidelberg Collaboratory for Image Processing, Heidelberg University, Berliner Straße 43, 69120 Heidelberg, Germany
• apreviously at: Institute of Environmental Physics, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
Abstract

Heat transfer velocities measured during three different campaigns in the Baltic Sea using the active controlled flux technique (ACFT) with wind speeds ranging from 5.3 to 14.8 m s−1 are presented. Careful scaling of the heat transfer velocities to gas transfer velocities using Schmidt number exponents measured in a laboratory study allows us to compare the measured transfer velocities to existing gas transfer velocity parameterizations, which use wind speed as the controlling parameter. The measured data and other field data clearly show that some gas transfer velocities are much lower than those based on the empirical wind speed parameterizations. This indicates that the dependencies of the transfer velocity on the fetch, i. e., the history of the wind and the age of the wind-wave field, and the effects of surface-active material need to be taken into account.

1 Introduction

The transfer of a trace gas across the air–sea interface is commonly characterized by the gas transfer velocity k, which links the gas flux j with the concentration difference across the interface, Δc:

$\begin{array}{}\text{(1)}& j=k\mathrm{\Delta }c.\end{array}$

Traditionally, k is parameterized with the wind speed measured at a height of 10 m, u10, since wind speed is the most readily available parameter. Different authors proposed different functional dependencies between k and u10, for example a gradual transition from a smooth to a wavy regime (Jähne1982) or piecewise linear , linear and quadratic , quadratic , or cubic terms .

gives an overview of the most commonly used techniques to measure the gas transfer velocity. In the last decades, the dual-tracer technique, especially with the tracer pair 3He−SF6, as well as eddy covariance measurements of the gases CO2 and dimethylsulfide (DMS), has become state of the art for measuring the gas transfer velocity in situ. A recent review article by proposed

$\begin{array}{}\text{(2)}& {k}_{\mathrm{600}}\phantom{\rule{0.33em}{0ex}}\left[\mathrm{cm}\phantom{\rule{0.125em}{0ex}}{\mathrm{h}}^{-\mathrm{1}}\right]=\mathrm{0.262}±\mathrm{0.022}{u}_{\mathrm{10}}^{\mathrm{2}}\phantom{\rule{0.33em}{0ex}}\left[{u}_{\mathrm{10}}\phantom{\rule{0.25em}{0ex}}\mathrm{in}\phantom{\rule{0.25em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{-\mathrm{1}}\right]\end{array}$

as the best fit to all available 3He−SF6 dual-tracer data points, where k600 denotes the transfer velocity scaled to a CO2-equivalent transfer velocity at 20 C. However, mass balance techniques such as the dual-tracer method have a large time constant of up to weeks and large spatial scales of a few tens of kilometers, smoothing away varying micrometeorological and surface conditions (e.g., the degree of surface contamination by surface-active material).

In contrast, the eddy covariance method provides measurements of the gas transfer velocity on timescales below 1 h and spatial scales of a few kilometers. However, bin averaging over wind speed intervals is frequently necessary, since even under idealized conditions, not all realizations of the turbulent field can be measured, so that each single flux measurement obtained during a 30 min time period is still uncertain .

In this study, the active controlled flux technique (ACFT), a thermographic technique, is used, which is capable of measuring the heat transfer velocity with a temporal resolution of about 20 min, which can then be scaled to gas transfer velocities. This technique is described in Sect. 3.1. The ACFT was deployed during three cruises in the Baltic Sea to investigate the variability of the transfer velocities under field conditions.

Earlier measurements of the gas transfer velocity in the Baltic Sea are sparse. used the eddy covariance technique to measure the transfer of CO2 in the Arkona Basin and used the same technique in the Gotland Sea. Both studies found a very high variability of the gas transfer velocity.

2 Factors influencing air–sea gas exchange

The common approach is to parameterize the gas transfer velocity with wind speed alone. However, a wealth of studies have shown that a multitude of factors influence gas transfer, for example the contamination of the water surface with surface-active material (e.g., ), bubble entrainment (e.g., ), fetch (e.g., ), rain (e.g., ) and convective mixing (e.g., ).

Since the method discussed in this paper is insensitive to bubble contributions and can only be used to measure the interfacial part of the air–sea gas transfer, and no measurements were performed in rain conditions, only the influence of surface-active material and fetch will be discussed here.

## 2.1 Surfactants

One factor contributing to the disagreement between gas transfer velocities measured at the same wind speed even with the same measuring technique are surface-active materials (surfactants), which reduce the gas transfer velocity. This reduction in the gas transfer velocity in the presence of surfactants is not caused by the additional diffusion of the gas through the monomolecular surfactant layer at the water surface but by hydrodynamic effects in the mass boundary layer. Surfactant presence at the water surface inhibits eddy motion close to the surface and reduces fluid velocities. Upwelling at the surface is hindered by a reduction in the surface divergence due to the viscoelastic properties of the surfactant . Vertical velocity fluctuations near the interface are considered vital to gas transfer enhancement. Decreased vertical transport of fresh fluid towards the water surface results in a thicker boundary layer and thus a reduced transfer velocity .

Surfactants are enriched in the sea surface microlayer in the world's oceans over a wide range of wind speeds as high as u10=13 m s−1 . ln the Baltic Sea, high surface activities were measured , with a seasonal dependency at a near-shore location. The reduction in the gas transfer velocity due to surfactants has been observed in studies, where the gas transfer velocity was measured in laboratory setups in fresh water with added artificial surfactants , in water sampled from the ocean , during field studies and during field studies where artificial surfactants were released on the ocean surface . Gas transfer is found to be highly variable, with a reduction of up to 60 % under surfactant influence.

The gas transfer velocity k of sparingly soluble gases is commonly parameterized with the friction velocity u*, a measure for momentum input,

$\begin{array}{}\text{(3)}& k=\frac{\mathrm{1}}{\mathit{\beta }}{u}_{*}{\mathit{Sc}}^{-n},\end{array}$

with the momentum transfer resistance parameter β and the Schmidt number exponent n . Both the momentum transfer resistance β and the Schmidt number exponent n depend on the hydrodynamic properties of the water surface. For a hydrodynamically smooth water surface, e.g., at very low wind speeds or under surfactant influence, the Schmidt number exponent is found to be $n=\mathrm{2}/\mathrm{3}$, while for a wavy water surface, $n=\mathrm{1}/\mathrm{2}$. For increasing friction velocity, this change from $n=\mathrm{2}/\mathrm{3}$ to 1∕2 is found to be smooth rather than sudden . In addition, this change in the Schmidt number exponent also depends on the contamination of the water surface with surface-active material, with the change starting at higher friction velocities and being steeper for a surfactant-covered water surface (Krall2013).

## 2.2 Fetch and wave age

Another factor influencing the gas transfer velocity, which is disregarded in the widely used wind speed only parameterizations, is the dependency on fetch or the age of the wave field. The earliest indications that the fetch is an important parameter were seen by , who used an 18 m long wind-wave tank and found almost a doubling of the gas transfer velocity compared to the earlier work by Liss (1973), who used a tank of only 4.5 m length. pointed out that the differences observed between gas transfer measurements in lakes and the ocean might be caused by growing wave fields and thus increasing near-surface turbulence over distances as great as a few hundreds of kilometers offshore. and developed a parameterization for the transfer velocity based on the breaking-wave parameter and the whitecap coverage, both of which depend on the fetch. The considerations above indicate that there should be a dependency of the gas transfer velocity on the fetch. But unfortunately there is no solid knowledge because more detailed measurements and theories are lacking.

3 Measuring technique

## 3.1 Active thermography

The active controlled flux technique (ACFT) can be used to measure gas transfer velocities under laboratory as well as under field conditions with a high temporal (minutes) and spatial (meters) resolution, using heat as a proxy tracer. A carbon dioxide laser with a scanning optic is used to deposit energy directly to the water surface. An infrared camera measures the resulting heating. For this study the system theory approach proposed in was used. In this approach, the laser is switched on and off with changing frequencies. At low laser forcing frequencies the water surface will reach the thermal equilibrium, resulting in constant heating. At higher forcing frequencies this equilibrium is not reached and the measured amplitude is damped. Using Fourier analysis to determine this amplitude damping depending on the laser forcing frequency, the time to reach the thermal equilibrium, which corresponds to the response time of the system, is calculated. It is linked to the transfer velocity by

$\begin{array}{}\text{(4)}& {k}_{\mathrm{heat}}=\sqrt{\frac{{D}_{\mathrm{heat}}}{\mathit{\tau }}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\text{or}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathit{\tau }=\frac{{D}_{\mathrm{heat}}}{{k}_{\mathrm{heat}}^{\mathrm{2}}}\end{array}$

. This analysis technique is particularly suitable for field measurements as it requires no absolute calibration. A more detailed description of the analysis method, the necessary correction for the penetration depth of the infrared camera and the error estimation can be found in .

## 3.2 Scaling heat transfer velocities to gas transfer velocities

To compare the measured transfer velocities of heat to the transfer velocities of a gas like CO2, Schmidt number scaling is applied:

$\begin{array}{}\text{(5)}& {k}_{\mathrm{gas}}={k}_{\mathrm{heat}}{\left(\frac{\mathit{Sc}}{\mathit{Pr}}\right)}^{-n},\end{array}$

where kgas and kheat are the transfer velocities for the gas and heat, respectively. The Schmidt number $\mathit{Sc}=\mathit{\nu }/{D}_{\mathrm{gas}}$ and the Prandtl number $\mathit{Pr}=\mathit{\nu }/{D}_{\mathrm{heat}}$ are given by the kinematic viscosity of the water divided by the diffusion coefficient of the gas and of heat in water, respectively. The Schmidt number exponent n varies between $n=\mathrm{2}/\mathrm{3}$ for a flat and $n=\mathrm{1}/\mathrm{2}$ for a wavy water surface . Schmidt number scaling is used to provide a value for the gas transfer velocity, which is independent of the specific measurement technique or tracer.

However, using heat as a proxy for a gas tracer has one significant drawback. Diffusion of heat is about 100 times faster than diffusion of a dissolved gas in water. Because of this, any uncertainty in the Schmidt number exponent n leads to a relatively large uncertainty for the heat transfer velocity scaled to a gas transfer velocity. It is generally given by

$\begin{array}{}\text{(6)}& \frac{\mathrm{\Delta }k}{{k}_{\mathrm{gas}}}=\mathrm{ln}\left(\frac{\mathit{Sc}}{\mathit{Pr}}\right)\mathrm{\Delta }n,\end{array}$

where Δk and Δn are the absolute uncertainties for the transfer velocity and the Schmidt number exponent, respectively. For the whole expected range of $n=\mathrm{2}/\mathrm{3}$ to 1∕2, ${\mathrm{\Delta }}_{n}=±\mathrm{0.083}$ (Fig. 1) and $\mathit{Sc}/\mathit{Pr}\approx \mathrm{600}/\mathrm{9}$, the relative scaling error is ±35 %. This is quite a large uncertainty.

In the past decade, several studies found deviations between heat scaled by the Schmidt number and the simultaneously measured gas transfer velocities . However, a more recent study by showed that using a model independent analysis method, as proposed by and the correct Schmidt number exponent results in good agreement.

For field measurements, the importance of using a Schmidt number exponent, depending on the water surface condition, is also highlighted in , who relate the gas transfer velocity to the turbulent energy dissipation rate.

Figure 1Possible ranges of Schmidt number exponents for a clean and surfactant-covered water surface as a function of the wind speed as inferred from experiments in the Heidelberg Aeolotron wind-wave tank (Krall2013) for the wind speeds encountered during this study. Friction velocities measured in the Aeolotron were taken from Bopp (2011) and converted to the wind speed at 10 m height using the drag coefficient parameterization by . To scale the heat transfer velocities measured in the present work, the mean values of the Schmidt number exponent were used.

Currently, there are no measurement techniques available to measure the Schmidt number exponent in the field with the same temporal resolution as the heat transfer measurements. Therefore, the scaling in the present work was done using Schmidt number exponents measured in the Heidelberg Aeolotron wind-wave tank (Krall2013), as opposed to , who used a fixed Schmidt number exponent of 1∕2. In , Schmidt number exponents were measured with different concentrations of the surface-active material (surfactant) Triton X-100. The mean of the Schmidt number exponent of the two extreme cases presented in , corresponding to clean water and water with 167 µg l−1 Triton X-100, respectively, was used to scale the heat transfer velocities to gas transfer velocities (see Fig. 1) to account for possible contamination of the water surface with surface-active material. The difference between the mean and both extreme values of the Schmidt number exponent was used as the uncertainty of the Schmidt number exponent. Since the Aeolotron wind-wave tank is an annular facility, it has virtually unlimited fetch, comparable with open-ocean conditions. Due to the lack of simultaneously measured Schmidt number exponents in the field, this approach is more realistic than using $n=\mathrm{1}/\mathrm{2}$ for all encountered wind conditions disregarding a potentially smooth condition ($n=\mathrm{2}/\mathrm{3}$) of the water surface. The approach used here reduces the uncertainty of Δn from ±0.083 to $<±\mathrm{0.030}$ (Fig. 1). The resulting relative uncertainty of k is then $\mathrm{\Delta }k/k<±\mathrm{13}$ %. Another source of uncertainty lies in transferring the lab measurements of the Schmidt number exponent to the field conditions, since the friction velocity u* is measured in the lab (Bopp2011) as opposed to the wind speed at 10 m height, which is commonly measured in the field. To convert lab measurements to field conditions, the drag coefficient, ${C}_{\mathrm{D}}={u}_{*}^{\mathrm{2}}{u}_{\mathrm{10}}^{-\mathrm{2}}$, taken from was used.

4 Measurements

## 4.1 Baltic Sea campaigns 2009 and 2010

Three ship campaigns were conducted in 2009 and 2010. Figure 2 show the tracks of these three cruises. The first one (Alkor Cruise 336, ) took place from 25 April until 7 May 2009 on the German RV Alkor. It included measurements northwest of Rügen and the Gotland Sea. The second cruise on the same vessel (Alkor Cruise 356, ) between 30 June and 13 July 2010 included measurement stations spread across the whole Baltic Sea. The third cruise took place on the Finnish research vessel RV Aranda from 14 until 19 September 2010. Due to the stormy weather conditions, most measurements were conducted in the Finnish archipelago and only two measurements were conducted under open-ocean conditions in the Gulf of Finland.

Figure 2Map of the Baltic Sea. The tracks of the three cruises are shown.

## 4.2 Experimental setup on ship

To use the ACFT method described in Sect. 3.1, a CO2 laser (Firestar f200, Synrad, Inc.) was used to heat the water surface. A scanning system (Micro Max 671, Cambridge Technology, Inc.) was used to widen the laser to create a heated patch on the water surface. The temperature response of the water surface was recorded with an infrared camera (CMT 256, Thermosensorik). Laser, scanner and camera are synchronized by custom electronics. A watertight box, including the infrared (IR) laser, the IR camera and the electronics, was installed on rails on top of an aluminum cradle at the bow of the research vessels. During transit times the box was retracted and fixed over the vessel, while it was moved over the ocean during measurement times. A more detailed description of all instruments used is given in .

Measurements were only conducted at stations where the vessel was standing in one position. Nevertheless due to currents the water surface moved relative to the ship. As direct sun irradiation disturbs the infrared signals, most measurement were conducted during nighttime or on cloudy days. Nevertheless, reflections of the thermal signature of the sky and the ship itself cannot be avoided. However, the periodic forcing of the heat flux as described in Sect. 3.1 suppresses these effects (lock-in technique).

Wind speed measured at 10 m height was provided by the weather station of each vessel. On RV Alkor, 1 min mean wind speeds were stored only for the times during which measurements with the ACFT were performed. On RV Aranda, 10 s mean values were stored for the whole duration of the cruise. During data processing, averages of the stored values were calculated for the times during which the respective ACFT measurements were performed.

5 Results

## 5.1 Measured transfer velocities

The first results of the cruise in 2009 are already published in . For this study a reevaluation with slight differences in the correction of the penetration depth of the infrared camera was done. Also, the improved Schmidt number scaling described in Sect. 3.2 was used, while used $n=\mathrm{1}/\mathrm{2}$ for all conditions. The obtained heat transfer velocities are given in Table A1. Figure 3 shows the measured transfer velocities, scaled to a Schmidt number of 600. To compare the results with other field measurements, the parameterization by , which parameterizes the transfer velocity with the wind speed is also shown. This parameterization was chosen for comparison, since it is one of the few in which a margin of uncertainty is included (gray band in Fig. 3).

Figure 3Measured k600 transfer velocities plotted against the wind speed of the RV Alkor spring 2009 cruise. For comparison the best fit of , Eq. (2), is added.

Figure 4 shows the measured heat transfer velocities against the wind speed for the Alkor campaign in 2010 in comparison with the parameterization by . Schmidt number scaling was done with the same method as for the Alkor 2009 data set. During most of the RV Alkor campaign in 2010 the wind speeds were rather low. At low wind speeds, the response time of the water surface is very long, as it increases with the square of the inverse transfer velocity (Eq. 4). The time a water parcel stays in the heated patch (residence time) is limited due to surface currents and the movement of the ship relative to the water surface. In the thermal equilibrium, the heat energy deposited on the water surface by the laser equals the energy removed from the surface by processes driving heat transfer, which results in a constant water surface temperature. Only if the residence time is longer than the response time does the water surface reach the thermal equilibrium. Otherwise a lower temperature and therefore a higher amplitude damping will be observed, which leads to an overestimation of the measured transfer velocities. The residence times were estimated from the infrared images themselves by measuring the time a single structure stayed in the heated patch. All measurements with wind speeds of 4 m s−1 and below are not reliable because the estimated residence times were found to be too long. Therefore, they will be excluded from further analysis.

Figure 4Measured k600 transfer velocities plotted against the wind speed of the RV Alkor summer 2010 cruise. Conditions for which the measured transfer velocity is likely overestimated are marked with open circles and will not be used for further analysis. For comparison the wind speed parameterization taken from is added.

This highlights the difficulties of measuring gas transfer velocities at very low wind speeds. However, difficulties also exist with other approaches to measure the gas transfer velocity in the field, such as dual-tracer studies, where the timescales required for measurements are very long at low wind speeds and sufficiently long periods of low winds are rarely encountered.

Figure 5Measured k600 transfer velocities plotted against the wind speed of the RV Aranda fall 2010 cruise. The filled circles show the open-ocean measurements, while the open circles are data from the archipelago. For comparison, the wind speed parameterization by is also shown.

The heat transfer velocities scaled to Sc=600 measured on RV Aranda in 2010 are shown in Fig. 5. The transfer velocities measured in the shielded archipelago are significantly lower than the ones measured under open-ocean conditions.

## 5.2 Comparison with other field and laboratory data

Figure 6 shows a comparison between the measured transfer velocities and the empirical parameterization of . The measurements from the Alkor 2009 and Alkor 2010 cruises coincide within the error margins with the empirical parameterization by Ho, except for the value at the highest wind speed, which is approx. 40 % lower. The two open-ocean measurements during the RV Aranda cruise 2010 are slightly lower than the empirical parameterization but still close to it.

This is, however, not the case for the RV Aranda cruise measurements in the shielded archipelago. The measured values are significantly lower. On average, the values are only about one half of the transfer velocities predicted by the empirical parameterization. There are three possible explanations for this finding: bubble-mediated transfer, fetch (or wave age) and surfactants. In the following sections, these possible explanations will be discussed in detail.

Figure 6Comparison of scaled heat transfer velocities measured in the Baltic Sea and gas transfer velocities measured in the Heidelberg Aeolotron wind-wave facility with a clean water surface (green shaded area). The measurements on RV Aranda in 2010 which were made under open-ocean conditions (i.e., with a virtually unlimited fetch) are marked with filled circles, while the fetch limited measurements in the archipelago are marked with open circles. Also shown is the lower limit for a smooth water surface; Eq. (3). The region between the transfer velocities measured with a clean water surface as the upper boundary and the values for a smooth water surface as the lower boundary for possible transfer velocities is shaded in magenta. Also shown are the data set from the North Atlantic of (K&R1983) using the radon deficit method, DMS eddy covariance measurements and the parameterization of previous Baltic Sea gas transfer measurements by . The individual data points in and scatter too strongly to be shown here. Also shown is the parameterization by .

### 5.2.1 Bubble-mediated transfer

It is known that active thermography misses the contribution by bubbles to the transfer; see Sect. 2. Because of the high solubility of dimethylsulfide (DMS), this tracer's gas transfer has almost no bubble-induced component, and the transfer velocities of DMS measured by do indeed have values very similar to ours (Fig. 6). Another observation, which supports this argument, are the higher CO2 gas exchange values (CO2 has a significantly lower solubility than DMS with a higher expected bubble-induced contribution) measured in the Baltic Sea by and . We only show the combined linear or quadratic parameterization by , since does not give a parameterization.

A very helpful indication comes, however, from laboratory experiments, which suggest that this explanation is not correct. No evidence for a significant bubble contribution to gas transfer was found in a laboratory study (Krall2013) up to the highest wind speed used in that study (≈12 m s−1), although tracers with solubilities much lower than CO2 (dimensionless solubility α≈0.7) and DMS (α≈11.2) were used, including N2O (α≈0.5), trifluoromethane (α≈0.26) and pentafluoroethane (α≈0.07). In another study, found no differences between gas transfer velocities of N2O and heat transfer velocities for wind speeds as high as 12 m s−1, which indicates that bubble contribution for both the transfer heat and that of N2O is not significant.

### 5.2.2 Fetch and wave-age effects

A second explanation would be the effect of fetch or, equivalently, quickly varying wind conditions with young wave ages. This effect has almost not be studied so far. Only recently, using active thermography measurements in the Heidelberg Aeolotron, showed that, at very short fetches and low wind speed, the gas transfer velocity is significantly higher than at infinite fetch. This finding is supported by an old data set, which constitutes the most diligently measured gas transfer velocities using the radon deficit method . One part of this data set was measured during the Joint Air Sea Interaction (JASIN) cruise in the North Atlantic with highly varying wind speeds. The measured gas transfer velocities are higher or as high as predicted by the empirical parameterization. However, the transfer velocities measured during the First GARP (Global Atmospheric Research Program) Global Experiment (FGGE) cruise with constantly blowing trade winds are significantly lower. One value is 3 times lower than predicted by the empirical parameterization of Ho. These measurements clearly indicate that even on the open ocean (i. e. without fetch limitations), there will be significant differences in the gas transfer velocity. The data suggest that this effect may be as large as a factor of 5.

Surprisingly, the thermographic measurements in the Baltic Sea show just the opposite dependency. In the shielded archipelago with fetches that are probably short, the transfer velocities are lower and not higher. Thus fetch dependency does not seem to be the correct explanation in this case at rather high wind speeds, where the Aeolotron data by also show no significant fetch dependency.

### 5.2.3 Surfactants

The third and most likely reason for the lower gas exchange rates during part of the Aranda 2010 cruise is a reduction in the transfer velocity by surface films. The reduction of about a factor of 2 is consistent with earlier measurements discussed in Sect. 2.1.

At this point it is helpful to compare the field data with laboratory data augmented by physical arguments about the mechanisms of air–sea gas transfer. However, a direct comparison is physically not valid because the conditions concerning the wave field and surface contamination will be different. Despite that, laboratory data are suited to explore the possible upper and lower limits of the gas transfer velocity at a given wind speed. The Heidelberg Aeolotron laboratory has a virtually unlimited fetch due to its annular shape, so it may resemble the ocean conditions in the best possible way. The gas transfer velocities measured when the water surface in the Aeolotron was carefully cleaned by skimming the top layer of the water before the start of each measurement to remove surface-active material can be considered to be the upper limit (green shaded area in Fig. 6). Those gas transfer velocities were measured with the method described in and are published in . In the green shaded area, the increase in the gas transfer velocities at low wind speeds and short fetches observed by is taken into account, too.

The lower limit for possible gas transfer velocities is given by the prediction of (Eq. 3 with $n=\mathrm{2}/\mathrm{3}$ and β=12.1) for a smooth water surface. These values have been confirmed by measurements in a small annular wind or wave facility, when the water surface was covered by surfactants . The highest friction velocity in water at which the water surface remained smooth and without wind waves in this facility was 1.4 cm s−1 corresponding to a smooth water surface up to a wind speed of u10≈13 m s−1. This is supported by the findings of , who measured surfactant enrichment in the sea surface microlayer up to u10≈13 m s−1 as well.

The region between these upper and lower bounds for gas transfer is shaded in a magenta color in Fig. 6. This difference between highest and lowest possible gas transfer velocities alone indicates that the gas transfer is highly variable and not only dependent on wind speed alone.

All field data shown and based on mass balance methods, eddy covariance and active thermography are compatible with this shaded region of possible gas transfer velocities. The parameterization of CO2 measured with the eddy covariance technique in the Baltic Sea according to is slightly higher than the upper limit resulting from laboratory measurements. Because of the high scatter of these data, some individual measurements are even much higher. This means that we still see discrepancies between measurements based on mass balance (now including also active thermography) and eddy covariance measurements, although they are not as bad as in the early days of eddy covariance measurements .

6 Conclusions and outlook

Heat exchange measurements were conducted in the Baltic Sea during three different campaigns using the active controlled flux technique. The measured heat transfer velocities, scaled to gas transfer velocities using realistic Schmidt number exponents, show high variability even at the same wind speed. It is new that, even at high wind speeds in the range of 8 to 15 m s−1, significantly lower gas transfer velocities were measured, which were about a factor of 2 lower than the average transfer velocities measured by the dual-tracer technique and parameterized by the relation of . Based on arguments from several lab studies, the influence of surfactants is the most likely reason for variability of the gas transfer velocity under the environmental conditions for the thermographic measurements in the Baltic Sea. But a possible influence of fetch and bubbles on these measurements cannot completely be ruled out.

Therefore this study clearly indicates that a better understanding of air–sea gas transfer urgently requires more systematic measurements of the effects of bubbles, fetch (or the age of the wave field) and surfactants. In the field, the most promising approach is eddy covariance measurements together with active thermography. For laboratory measurements some serious limitations must be overcome. One is the fetch gap. In linear facilities only very short fetches can be studied, which are no longer than the maximum length of the water tunnel in the facility. Even at these short fetches, significant variations in the gas transfer rate can be measured. This has recently been demonstrated by using active thermography.

In order to increase the fetch range available in the lab, gas exchange measurements could be performed in annular facilities under unsteady wind speed conditions. In the Heidelberg Aeolotron it is possible to switch on the wind in a few seconds, while it takes several minutes for the wave field to develop to a stationary state. Unfortunately, it is very hard to make gas exchange measurements with a temporal resolution of below 1 min using conventional mass balance techniques.

A very promising technique for fast measurements of gas transfer is the recently developed mass boundary layer imaging technique . Using this technique will enable the measurement of the gas transfer velocity simultaneously and in the same footprint as the heat transfer velocity. This will allow a direct comparison as well as in-depth studies of the physical mechanisms governing air–sea gas and heat transfer.

Data availability
Data availability.

Data have been uploaded to https://pangaea.de/. All third-party data sets used are cited in the text.

Appendix A: Numerical values of the measured transfer velocities

Tables A1, A2 and A3 give the numerical values of the measurements conducted during the cruises in the Baltic Sea.

Table A1Measured heat transfer velocities kheat depending on time, position, wind speed, and water and air temperature for the measurements on RV Alkor in 2009. Furthermore the Prandtl number Pr, the Schmidt number exponent n and the scaled transfer velocity k600 are given. The times given are approximate starting times in UTC. Each measurement lasted about 20 min.

Table A2Measured heat transfer velocities kheat depending on time, position, wind speed, and water and air temperature for the measurements on RV Alkor in 2010. Furthermore the Prandtl number Pr, the Schmidt number exponent n and the scaled transfer velocity k600 are given. The times given are approximate starting times in UTC. Each measurement lasted about 20 min.

Measurements marked with an asterisk (*) are considered unreliable; see Sect. 5.1.

Table A3Measured heat transfer velocities kheat depending on time, position, wind speed, and water and air temperature for the measurements on RV Aranda in 2010. Furthermore the Prandtl number Pr, the Schmidt number exponent n and the scaled transfer velocity k600 are given. The times given are approximate starting times in UTC. Each measurement lasted about 20 min. All measurements were conducted in a fetch-limited position with the exception of the two conditions marked with an asterisk (*).

Author contributions
Author contributions.

LN was involved in designing the experiments, performed the experiments and calculated the heat transfer velocities from the IR images. KEK calculated the k600 values, prepared all figures and tables, and served as communicating author. BJ was involved in designing the experiments and the analytical methods and outlined the main conclusions. All authors contributed to writing the manuscript.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We would like to thank Kimmo Kahma and Heidi Pettersson, Finish Meteorological Institute, for the possibility of participating in the Aranda CO2_WAVE10_CTD10/2010 cruise. We would also like to thank Robert Schmidt and Bernd Schneider, Institut für Ostseeforschung, Warnemünde, for their chief scientist work during the cruises Alkor 336 and Alkor 356. We are grateful for the assistance onboard by the captains and crews of RV Aranda and RV Alkor. We would like to thank Uwe Schimpf and Günther Balschbach for help with preparing and conducting the measurements as well as logistical support. Financial support for this work by the German Federal Ministry of Education and Research (BMBF) joint project “Surface Ocean Processes in the Anthropocene” (SOPRAN, FKZ 03F0462F, 03F0611F and 03F0622F) within the international SOLAS project is gratefully acknowledged. We acknowledge financial support by Deutsche Forschungsgemeinschaft within the funding programme Open Access Publishing, by the Baden-Württemberg Ministry of Science, Research and the Arts and by Ruprecht-Karls-Universität Heidelberg.

Edited by: Mario Hoppema
Reviewed by: three anonymous referees

References

Asher, W. E., Jessup, A. T., and Atmane, M. A.: Oceanic application of the active controlled flux technique for measuring air-sea transfer velocities of heat and gases, J. Geophys. Res., 109, C08S12, https://doi.org/10.1029/2003JC001862, 2004. a

Atmane, M. A., Asher, W., and Jessup, A. T.: On the use of the active infrared technique to infer heat and gas transfer velocities at the air-water free surface, J. Geophys. Res., 109, C08S14, https://doi.org/10.1029/2003JC001805, 2004. a

Bell, T. G., De Bruyn, W., Miller, S. D., Ward, B., Christensen, K. H., and Saltzman, E. S.: Air-sea dimethylsulfide (DMS) gas transfer in the North Atlantic: evidence for limited interfacial gas exchange at high wind speed, Atmos. Chem. Phys., 13, 11073–11087, https://doi.org/10.5194/acp-13-11073-2013, 2013. a, b

Bell, T. G., De Bruyn, W., Marandino, C. A., Miller, S. D., Law, C. S., Smith, M. J., and Saltzman, E. S.: Dimethylsulfide gas transfer coefficients from algal blooms in the Southern Ocean, Atmos. Chem. Phys., 15, 1783–1794, https://doi.org/10.5194/acp-15-1783-2015, 2015. a, b

Bopp, M.: Messung der Schubspannungsgeschwindigkeit am Heidelberger Aeolotron mittels der Impulsbilanzmethode, Bachelor thesis, Institut für Umweltphysik, Fakultät für Physik und Astronomie, Univ. Heidelberg, Heidelberg, Germany, 2011. a, b

Brockmann, U. H., Hühnerfuss, H., Kattner, G., Broecker, H.-C., and Hentzschel, G.: Artificial surface films in the sea area near Sylt, Limnol. Oceanogr., 27, 1050–1058, https://doi.org/10.4319/lo.1982.27.6.1050, 1982. a

Broecker, H. C., Siems, W., and Petermann, J.: The influence of wind on CO2 exchange in a wind wave tunnel, including the effects of monolayers, J. Mar. Res., 36, 595–610, 1978. a

Broecker, W. S., Ledwell, J. R., Takahashi, T., Weiss, R., Merlivat, L., Memery, L., Jähne, B., and Münnich, K. O.: Isotopic versus micrometeorologic ocean CO2 fluxes: A serious conflict, J. Geophys. Res., 91, 10517–10528, https://doi.org/10.1029/JC091iC09p10517, 1986. a

Coantic, M.: A model of gas transfer across air–water interfaces with capillary waves, J. Geophys. Res., 91, 3925–3943, https://doi.org/10.1029/JC091iC03p03925, 1986. a

Crosswell, J. R.: Bubble Clouds in Coastal Waters and Their Role in Air-Water Gas Exchange of CO2, Journal of Marine Science and Engineering, 3, 866–890, https://doi.org/10.3390/jmse3030866, 2015. a

Csanady, G. T.: The role of breaking wavelets in air-sea gas transfer, J. Geophys. Res., 95, 749–759, https://doi.org/10.1029/JC095iC01p00749, 1990. a

Deacon, E. L.: Gas transfer to and across an air-water interface, Tellus, 29, 363–374, https://doi.org/10.3402/tellusa.v29i4.11368, 1977. a, b

Edson, J. B., Jampana, V., Weller, R. A., Bigorre, S. P., Plueddemann, A. J., Fairall, C. W., Miller, S. D., Mahrt, L., Vickers, D., and Hersbach, H.: On the exchange of momentum over the open ocean, J. Phys. Oceanogr., 43, 1589–1610, https://doi.org/10.1175/JPO-D-12-0173.1, 2013. a, b

Esters, L., Landwehr, S., Sutherland, G., Bell, T. G., Christensen, K. H., Saltzman, E. S., Miller, S. D., and Ward, B.: Parameterizing air-sea gas transfer velocity with dissipation, J. Geophys. Res., 122, 3041–3056, https://doi.org/10.1002/2016JC012088, 2017. a

Frew, N., Bock, E., Schimpf, U., Hara, T., Haussecker, H., Edson, J., McGillis, W., Nelson, R., McKenna, S., Uz, B., and Jähne, B.: Air-sea gas transfer: Its dependence on wind stress, small-scale roughness, and surface films, J. Geophys. Res., 109, C08S17, https://doi.org/10.1029/2003JC002131, 2004. a, b

Frew, N. M., Goldman, J. C., Denett, M. R., and Johnson, A. S.: Impact of phytoplankton-generated surfactants on air-sea gas-exchange, J. Geophys. Res., 95, 3337–3352, https://doi.org/10.1029/JC095iC03p03337, 1990. a, b

Frew, N. M., Bock, E. J., McGillis, W. R., Karachintsev, A. V., Hara, T., Münsterer, T., and Jähne, B.: Variation of air–water gas transfer with wind stress and surface viscoelasticity, in: Air-water Gas Transfer, Selected Papers from the Third International Symposium on Air-Water Gas Transfer, edited by: Jähne, B. and Monahan, E. C., 529–541, AEON, Hanau, https://doi.org/10.5281/zenodo.10405, 1995. a

Garbe, C. S., Rutgersson, A., Boutin, J., Delille, B., Fairall, C. W., Gruber, N., Hare, J., Ho, D., Johnson, M., de Leeuw, G., Nightingale, P., Pettersson, H., Piskozub, J., Sahlee, E., Tsai, W., Ward, B., Woolf, D. K., and Zappa, C.: Transfer across the air-sea interface, in: Ocean-Atmosphere Interactions of Gases and Particles, edited by: Liss, P. S. and Johnson, M. T., 55–112, Springer, Heidelberg, https://doi.org/10.1007/978-3-642-25643-1_2, 2014. a

Goldman, J. C., Dennet, M. R., and Frew, N. M.: Surfactant effects on air-sea gas exchange under turbulent conditions, Deep-Sea Res, 35, 1953–1970, https://doi.org/10.1016/0198-0149(88)90119-7, 1988. a

Harrison, E. L., Veron, F., Ho, D. T., Reid, M. C., Orton, P., and McGillis, W. R.: Nonlinear interaction between rain- and wind-induced air-water gas exchange, J. Geophys. Res.-Oceans, 117, C03034, https://doi.org/10.1029/2011JC007693, 2012. a

Ho, D. T., Wanninkhof, R., Schlosser, P., Ullman, D. S., Hebert, D., and Sullivan, K. F.: Toward a universal relationship between wind speed and gas exchange: Gas transfer velocities measured with 3He∕SF6 during the Southern Ocean Gas Exchange Experiment, J. Geophys. Res., 116, C00F04, https://doi.org/10.1029/2010JC006854, 2011. a, b, c, d, e, f, g, h, i

Jessup, A. T., Asher, W. E., Atmane, M., Phadnis, K., Zappa, C. J., and Loewen, M. R.: Evidence for complete and partial surface renewal at an air-water interface, Geophys. Res. Lett., 36, 1–5, https://doi.org/10.1029/2009GL038986, 2009. a

Jähne, B.: Trockene Deposition von Gasen über Wasser (Gasaustausch), in: Austausch von Luftverunreinigungen an der Grenzfläche Atmosphäre/Erdoberfläche, Zwischenbericht für das Umweltbundesamt zum Teilprojekt 1: Deposition von Gasen, BleV-R-64.284-2, edited by: Flothmann, D., Battelle Institut, Frankfurt, https://doi.org/10.5281/zenodo.10278, 1982. a

Jähne, B., Münnich, K. O., and Siegenthaler, U.: Measurements of gas exchange and momentum transfer in a circular wind-water tunnel, Tellus, 31, 321–329, https://doi.org/10.1111/j.2153-3490.1979.tb00911.x, 1979. a, b

Jähne, B., Münnich, K. O., Bösinger, R., Dutzi, A., Huber, W., and Libner, P.: On the parameters influencing air-water gas exchange, J. Geophys. Res., 92, 1937–1950, https://doi.org/10.1029/JC092iC02p01937, 1987. a, b, c

Jähne, B., Libner, P., Fischer, R., Billen, T., and Plate, E. J.: Investigating the transfer process across the free aqueous boundary layer by the controlled flux method, Tellus, 41, 177–195, https://doi.org/10.3402/tellusb.v41i2.15068, 1989. a, b, c

Krall, K. E.: Laboratory Investigations of Air-Sea Gas Transfer under a Wide Range of Water Surface Conditions, Dissertation, Institut für Umweltphysik, Fakultät für Physik und Astronomie, Univ. Heidelberg, Heidelberg, Germany, https://doi.org/10.11588/heidok.00014392, 2013. a, b, c, d, e, f, g, h, i

Kräuter, C.: Visualization of air-water gas exchange, Dissertation, Institut für Umweltphysik, Fakultät für Physik und Astronomie, Univ. Heidelberg, Heidelberg, Germany, https://doi.org/10.11588/heidok.00018209, 2015. a

Kräuter, C., Trofimova, D., Kiefhaber, D., Krah, N., and Jähne, B.: High resolution 2-D fluorescence imaging of the mass boundary layer thickness at free water surfaces, J. Europ. Opt. Soc. Rap. Public., 9, 14016, https://doi.org/10.2971/jeos.2014.14016, 2014. a

Kromer, B. and Roether, W.: Field measurements of air-sea gas exchange by the radon deficit method during JASIN 1978 and FGGE 1979, Meteor Forschungsergebnisse, Deutsche Forschungsgemeinschaft, Reihe A/B, Gebrüder Bornträger, Berlin, Stuttgart, Germany, 24, 55–76, https://doi.org/10.1594/PANGAEA.604844, 1983. a, b

Kunz, J. and Jähne, B.: Investigating small scale air-sea exchange processes via thermography, Front. Mech. Eng., 4, 4, https://doi.org/10.3389/fmech.2018.00004, 2018. a, b, c, d

Lee, R. J. and Saylor, J. R.: The effect of a surfactant monolayer on oxygen transfer across an air/water interface during mixed convection, Int. J. Heat Mass Tran., 53, 3405–3413, https://doi.org/10.1016/j.ijheatmasstransfer.2010.03.037, 2010. a

Liss, P. S.: Processes of gas exchange across an air–water interface, Deep-Sea Res., 20, 221–238, https://doi.org/10.1016/0011-7471(73)90013-2, 1973. a

Liss, P. S. and Merlivat, L.: Air-sea gas exchange rates: Introduction and synthesis, in: The Role of Air-Sea Exchange in Geochemical Cycling, edited by: Buat-Menard, P., 113–129, Reidel, Boston, MA, https://doi.org/10.1007/978-94-009-4738-2_5, 1986. a

McKenna, S. P. and Bock, E. J.: Physicochemical effects of the marine microlayer on air-sea gas transport, in: Marine Surface Films: Chemical Characteristics, Influence on Air-Sea Interactions, and Remote Sensing, edited by: Gade, M., Hühnerfuss, H., and Korenowski, G. M., 77–91, Springer Berlin Heidelberg, Germany, https://doi.org/10.1007/3-540-33271-5_9, 2006. a

McKenna, S. P. and McGillis, W. R.: The role of free-surface turbulence and surfactants in air-water gas transfer, Int. J. Heat Mass Tran., 47, 539–553, https://doi.org/10.1016/j.ijheatmasstransfer.2003.06.001, 2004. a

Mesarchaki, E., Kräuter, C., Krall, K. E., Bopp, M., Helleis, F., Williams, J., and Jähne, B.: Measuring air–sea gas-exchange velocities in a large-scale annular wind–wave tank, Ocean Sci., 11, 121–138, https://doi.org/10.5194/os-11-121-2015, 2015. a, b

Nagel, L.: Active Thermography to Investigate Small-Scale Air-Water Transport Processes in the Laboratory and the Field, Dissertation, Institut für Umweltphysik, Fakultät für Chemie und Geowissenschaften, Univ. Heidelberg, Heidelberg, Germany, https://doi.org/10.11588/heidok.00016831, 2014. a, b

Nagel, L., Krall, K. E., and Jähne, B.: Comparative heat and gas exchange measurements in the Heidelberg Aeolotron, a large annular wind-wave tank, Ocean Sci., 11, 111–120, https://doi.org/10.5194/os-11-111-2015, 2015. a, b

Nightingale, P. D., Malin, G., Law, C. S., Watson, A. J., Liss, P. S., Liddicoat, M. I., Boutin, J., and Upstill-Goddard, R. C.: In situ evaluation of air-sea gas exchange parameterization using novel conservation and volatile tracers, Global Biogeochem. Cy, 14, 373–387, https://doi.org/10.1029/1999GB900091, 2000. a

Pereira, R., Ashton, I., Sabbaghzadeh, B., Shutler, J. D., and Upstill-Goddard, R. C.: Reduced air–sea CO2 exchange in the Atlantic Ocean due to biological surfactants, Nat. Geosci., 11, 492–496, https://doi.org/10.1038/s41561-018-0136-2, 2018. a

Richter, K. and Jähne, B.: A laboratory study of the Schmidt number dependency of air-water gas transfer, in: Gas Transfer at Water Surfaces 2010, edited by: Komori, S., McGillis, W., and Kurose, R., Kyoto Univ. Press, Kyoto, Japan, 322–332, https://doi.org/10.5281/zenodo.14955, 2011. a, b

Roether, W. and Kromer, B.: Optimum application of the radon deficit method to obtain air–sea gas exchange rates, in: Gas transfer at water surfaces, edited by: Brutsaert, W. and Jirka, G. H., 447–457, Reidel, Hingham, MA, https://doi.org/10.1007/978-94-017-1660-4_41, 1984. a

Rutgersson, A., Norman, M., Schneider, B., Pettersson, H., and Sahlée, E.: The annual cycle of carbon-dioxide and parameters influencing the air-sea carbon exchange in the Baltic Proper, J. Marine Syst., 74, 381–394, https://doi.org/10.1016/j.jmarsys.2008.02.005, 2008. a, b, c, d

Rutgersson, A., Smedman, A., and Sahlée, E.: Oceanic convective mixing and the impact on air-sea gas transfer velocity, Geophys. Res. Lett., 38, L02602, https://doi.org/10.1029/2010GL045581, 2011. a

Sabbaghzadeh, B., Upstill-Goddard, R. C., Beale, R., Pereira, R., and Nightingale, P. D.: The Atlantic Ocean surface microlayer from 50 N to 50 S is ubiquitously enriched in surfactants at wind speeds up to 13 m s−1, Geophys. Res. Lett., 44, 2852–2858, https://doi.org/10.1002/2017GL072988, 2017GL072988, 2017. a, b

Salter, M., Upstill-Goddard, R., Nightingale, P., Archer, S., Blomquist, B., Ho, D., Huebert, B., Schlosser, P., and Yang, M.: Impact of an artificial surfactant release on air-sea gas fluxes during Deep Ocean Gas Exchange Experiment II, J. Geophys. Res., 116, C11016, https://doi.org/10.1029/2011JC007023, 2011. a, b

Schimpf, U., Nagel, L., and Jähne, B.: The 2009 SOPRAN active thermography pilot experiment in the Baltic Sea, in: Gas Transfer at Water Surfaces 2010, edited by: Komori, S., McGillis, W., and Kurose, R., Kyoto Univ. Press, Kyoto, Japan, 358–367, https://doi.org/10.5281/zenodo.14956, 2011. a, b, c

Schmidt, R.: Cruise Report of Alkor Cruise 336, available at: https://www2.bsh.de/aktdat/dod/fahrtergebnis/2009/20100107.htm (last access: 4 March 2019), 2009. a

Schmidt, R. and Schneider, B.: The effect of surface films on the air–sea gas exchange in the Baltic Sea, Mar. Chem., 126, 56–62, https://doi.org/10.1016/j.marchem.2011.03.007, 2011. a, b

Schneider, B.: Cruise Report of Alkor Cruise 356, available at: https://www2.bsh.de/aktdat/dod/fahrtergebnis/2010/20110109.htm (last access: 4 March 2019), 2010. a

Toba, Y. and Koga, M.: A parameter describing overall conditions of wave breaking, whitecapping, sea-spray production, and wind stress, in: Oceanic Whitecaps and Their Role in Air-Sea Exchange Processes, 37–47, D. Reidel Publishing Company, Springer, Dordrecht, https://doi.org/10.1007/978-94-009-4668-2_4, 1986. a

Wanninkhof, R.: Relationship between wind speed and gas exchange over the ocean, J. Geophys. Res., 97, 7373–7382, https://doi.org/10.1029/92JC00188, 1992.  a, b

Wanninkhof, R. and McGillis, W. R.: A cubic relationship between gas transfer and wind speed., Geophys. Res. Lett., 26, 1889–1892, https://doi.org/10.1029/1999GL900363, 1999. a

Wanninkhof, R., Asher, W. E., Ho, D. T., Sweeney, C., and McGillis, W. R.: Advances in quantifying air-sea gas exchange and environmental forcing, Annu. Rev. Mar. Sci., 1, 213–244, https://doi.org/10.1146/annurev.marine.010908.163742, 2009. a

Weiss, A., Kuss, J., Peters, G., and Schneider, B.: Evaluating transfer velocity-wind speed relationship using a long-term series of direct eddy correlation CO2 flux measurements, J. Marine Syst., 66, 130–139, https://doi.org/10.1016/j.jmarsys.2006.04.011, 2007. a, b, c, d, e, f

Woolf, D., Leifer, I., Nightingale, P., Rhee, T., Bowyer, P., Caulliez, G., de Leeuw, G., Larsen, S., Liddicoat, M., Baker, J., and Andreae, M.: Modelling of bubble-mediated gas transfer: Fundamental principles and a laboratory test, J. Marine Syst., 66, 71–91, https://doi.org/10.1016/j.jmarsys.2006.02.011, 2007. a

Woolf, D. K.: Parametrization of gas transfer velocities and sea-state-dependent wave breaking, Tellus B, 57, 87–94, https://doi.org/10.1111/j.1600-0889.2005.00139.x, 2005. a, b

Wurl, O., Miller, L., and Vagle, S.: Production and fate of transparent exopolymer particles in the ocean, J. Geophys. Res., 116, C00H13, https://doi.org/10.1029/2011JC007342, 2011. a

Zappa, C., Ho, D., McGillis, W., Banner, M., Dacey, J. W. H., Bliven, L., Ma, B., and Nystuen, J.: Rain-induced turbulence and air-sea gas transfer, J. Geophys. Res., 114, C07009, https://doi.org/10.1029/2008JC005008, 2009. a

Zappa, C. J., Asher, W. E., Jessup, A. T., Klinke, J., and Long, S. R.: Microbreaking and the enhancement of air-water transfer velocity, J. Geophys. Res., 109, C08S16, https://doi.org/10.1029/2003JC001897, 2004. a

Zhao, D., Toba, Y., Suzuki, Y., and Komori, S.: Effect of wind waves on air-sea gas exchange: proposal of an overall CO2 transfer velocity formula as a function of breaking-wave parameter, Tellus, 55B, 478–487, https://doi.org/10.1034/j.1600-0889.2003.00055.x, 2003. a, b