Next Article in Journal
Post-Monsoon Season Precipitation Reduction over South Asia: Impacts of Anthropogenic Aerosols and Irrigation
Next Article in Special Issue
Spectroscopic BIL-SFG Invariance Hides the Chaotropic Effect of Protons at the Air-Water Interface
Previous Article in Journal
Seasonal Changes of Extremes in Isolated and Mesoscale Precipitation for the Southeastern United States
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The FuGas 2.3 Framework for Atmosphere–Ocean Coupling: Comparing Algorithms for the Estimation of Solubilities and Gas Fluxes

by
Vasco M. N. C. S. Vieira
1,*,
Pavel Jurus
2,3,
Emanuela Clementi
4 and
Marcos Mateus
1
1
MARETEC–Marine Environment and Technology, Instituto Superior Técnico, Universidade de Lisboa, Av Rovisco Pais, 1049-001 Lisboa, Portugal
2
DataCastor, U Svobodarny 1063/6, 190 00 Praha 9 Prague, Czech Republic
3
Institute of Computer Science, Czech Academy of Sciences, 11720 Prague, Czech Republic
4
Istituto Nazionale di Geofisica e Vulcanologia, INGV, 40128 Bologna, Italy
*
Author to whom correspondence should be addressed.
Atmosphere 2018, 9(8), 310; https://doi.org/10.3390/atmos9080310
Submission received: 8 May 2018 / Revised: 19 July 2018 / Accepted: 2 August 2018 / Published: 9 August 2018
(This article belongs to the Special Issue Physical Chemistry of the Air-Water Interface)

Abstract

:
Accurate estimates of the atmosphere–ocean fluxes of greenhouse gases and dimethyl sulphide (DMS) have great importance in climate change models. A significant part of these fluxes occur at the coastal ocean which, although much smaller than the open ocean, have more heterogeneous conditions. Hence, Earth System Modelling (ESM) requires representing the oceans at finer resolutions which, in turn, requires better descriptions of the chemical, physical and biological processes. The standard formulations for the solubilities and gas transfer velocities across air–water surfaces are 36 and 24 years old, and new alternatives have emerged. We have developed a framework combining the related geophysical processes and choosing from alternative formulations with different degrees of complexity. The framework was tested with fine resolution data from the European coastal ocean. Although the benchmark and alternative solubility formulations generally agreed well, their minor divergences yielded differences of up to 5.8% for CH4 dissolved at the ocean surface. The transfer velocities differ strongly (often more than 100%), a consequence of the benchmark empirical wind-based formulation disregarding significant factors that were included in the alternatives. We conclude that ESM requires more comprehensive simulations of atmosphere–ocean interactions, and that further calibration and validation is needed for the formulations to be able to reproduce it. We propose this framework as a basis to update with formulations for processes specific to the air–water boundary, such as the presence of surfactants, rain, the hydration reaction or biological activity.

1. Introduction

The dynamics of the atmosphere-ocean gas exchanges plays a central role in the Earth’s climate, with the ocean acting simultaneously as a sink and source of greenhouse gases and a source of dimethyl sulfide (DMS) to the atmosphere. However, and despite the observed seasonal, inter-annual and regional variability [1,2,3], the open ocean is generally assumed to act as a sink for atmospheric CO2. In the sub-polar regions, for instance, the solubility pump retrieves large amounts of greenhouse gases from the atmosphere, which are then advected to the deep ocean [1,2]. On the other hand, at the surface of coastal waters the balances and fluxes of CO2, CH4, N2O and DMS are very heterogenic, mostly due to factors like upwelling, plankton productivity, and land-originated loads of both natural and anthropogenic cause. Consequently, besides DMS [3,4] the coastal ocean can be, at least seasonally, also a source of CO2 [1,2,5,6,7,8,9,10,11], CH4 [2,10,12], and N2O [2,10,13,14,15] to the atmosphere. Thus, the uncertainty on the atmosphere–coastal ocean exchange of CO2 begs the question about the overall budget of the global ocean [1,3,6,16].
The atmosphere–ocean flux of gases is controlled by the unbalance between the air and water concentrations (the latter mediated by solubility), and by the transfer velocity [17,18,19,20,21,22]. The solubility of gases at the surface can be estimated by algorithms with different physical and chemical backgrounds [20,21,22,23,24,25,26,27]. The transfer velocity depends on numerous processes occurring at the sea surface [17,18,19,20,21,22,28]. However, large-scale simulations of relevant processes for the air–water gas exchange at the coastal ocean and their inherent variability poses serious difficulties for Earth System Models (ESM) and marine biogeochemical models. Constraints related to computational demands explain why such models usually simulate the biosphere at decadal or centennial time scales with daily intervals and spatial resolutions of hundreds to 1000 km. Such are the cases of ESM applications by the Intergovernmental Panel on Climate Change (IPCC), Max Planck Institute (MPI) or Centro Euro-Mediterraneo sui Cambiamenti Climatici (CMCC). Hand in hand with the low resolution for space and time, the atmosphere–ocean gas fluxes are estimated with simple algorithms that disregard the complexity of processes recently unveiled at the coastal ocean. Hence, the standard in current ESM is the parameterization by Wanninkhof [29], relying on wind speed at 10 m elevation under atmospherically neutral conditions (u10N) as the sole driver of transfer velocity. ESM relies on the dominance of strong winds over the global ocean to support the coarse resolution of geophysical process. Contradicting this rationale, Wanninkhof et al. [19] demonstrated that the largest amount of atmosphere–ocean CO2 exchange occurs under low-to-moderate winds. Concomitantly, there emerged empirical formulae best fitting low wind data collected from estuaries. Such are the formulations proposed by Carini et al. [30] or Raymond and Cole [31], accounting only for u10, or by Borges et al. [32], adding current drag with the bottom. Meanwhile, new evidence emerged for the effects of the sea surface cool-skin and warm-layer [33,34,35,36,37,38], atmospheric stability [39,40], wave-breaking [41,42,43,44,45,46,47,48,49,50], water surface renewal [28,51,52,53,54,55,56,57], surfactants [28,58,59,60] and rain [61,62,63,64], just to mention a few examples from an extensive list of published works. In addition to their single influence, these factors also interact. For instance, while surfactants increase surface tension, although decreasing the water surface renewal and micro-scale wave breaking [28,65], there is no report of them affecting whitecap from breaking waves. Surfactants are expected to be more relevant on the coastal ocean due to their supply from land and their effect being more pronounced under low–moderate winds [60,66].
There are many simplistic formulae estimating the gas transfer velocity from the few factors that were determinant for the specific set of data used in their calibration. However, modelling the coastal oceans with fine resolution requires algorithms providing adequate representation of the many processes involved. Developing such an algorithm demands a modular framework able to be updated with the best formulation for each of the mediator processes, and to test and compare different algorithms as well as different degrees of complexity. Hence, we developed the FuGas (Flux of Gases) framework [21,22], currently in its version 2.3. This simulation tool includes 67 alternative formulations to account for factors such as solubility, wind or current-mediated turbulence, atmospheric stability, sea-surface roughness, breaking waves, air and water viscosities, temperature and salinity. This framework is also able to automatically select simpler algorithms in particular locations with insufficient data for more comprehensive solutions. The framework can be used to estimate the air–water fluxes of any atmospheric gas by selecting the respective constants. In this work we rely on the FuGas 2.3 to compare estimates of CO2 transfer velocities obtained by simpler empirical formulations with some more comprehensive alternatives focusing on the sea state. We tested the effect of wave-state with field data from the Baltic Sea. Also, we coupled the Weather Research and Forecasting (WRF) atmospheric model to the WaveWatch III (WW3)–NEMO oceanographic model using simulated data from the European coastal ocean. These new tests are particularly relevant given recent findings that the sea state in fetch-limited locations (as is the coastal ocean) has an influence on the transfer velocity, by contrast with that from the fetch-unlimited open ocean [50,55,57].

2. Methods

Air–water gas fluxes result from the interaction of two factors: (i) the imbalance between the gas concentrations in the air and in the water sets the strength and direction of the flux; and (ii) the resistance the medium does for being crossed by the flow [17,18,19,20,21,22]. Most formulae estimate the flux from F = kw∙kHcp∙Δpgas. The kHcp is the Henry’s constant for gas solubility in its Cw/pa form, where pa is the gas partial pressure in the air and Cw its concentration in the water. The Δpgas is the difference between the partial pressures of the gas in the air (pa) and in the water (pw) i.e., Δpgas = pw − pa. In this form, a positive flux represents uptake by the atmosphere while a negative flux represents uptake by the ocean. The pa must be corrected for the partial pressure of water vapour considering saturation over the sea surface: pmoist = (1 − pH2O/P) pdry. This conversion is detailed in Section 2.1 below. The kw is the transfer velocity of gases across the sub-millimetre thick water surface layer. For moderately soluble gases the airside resistance is not negligible. In these cases, the double layer model [17,18,19,20,21,22] estimates the flux taking into consideration both the water-side and air-side surface layers and thus, F = Kw (Cw − Ca/kH) = Ka (Cw∙kH − Ca). The Ca and Cw are the concentrations of the gas in air and water and the kH is Henry’s constant in its equivalent dimensionless quantity (Ca/Cw). The transfer velocity is estimated from both layers as Kw = (1/kw + 1/(kH∙ka))−1 or its equivalent Ka = (kH/kw + 1/ka)−1. Supplementary Table S1 provides a list of the variables and constants presented in this article. Other variables and constants used in subsidiary calculations are not included. The link to software and data is provided in the Supplementary Materials.

2.1. Solubility

Sarmiento and Gruber [2] compiled the algorithm for the kHcp dependence on temperature and salinity provided by Weiss [23] and Weiss and Price [24]. We converted it to its corresponding dimensionless kH preserving the constants required to estimate Bunsen’s solubility coefficient (β) (Equation (1) or Equation (2)). Solubility coefficients were estimated from the Virial expansion (Equation (3)), where Β was β or β/Vm, depending on which gas it was applied to (Table 3.2.2 in Sarmiento and Gruber [2]). Our software automatically detected the gas from the ai coefficient. When Β = β the kH was estimated from Equation (1), whereas when Β = β/Vm the kH was estimated from Equation (2). This formulation accounted for fugacity (f) of non-ideal gases (Equation (4)) and corrected the gas partial pressure for moisture effects from pmoist = (1 − pH2O/P) pdry considering water vapour saturation over the sea-surface (Equation (5)). Calculations required the air pressure (P), water temperature (Tw), salinity (S), ideal gas law constant (R), molar volume of the specific gas (Vm) and the molar volume of ideal gases (Videal).
k H = ( 1 p H 2 O P ) 101.325 V m RT w β f
k H = 101.325 RT w β f
log ( Β ) = a 1   +   a 2 100 T w   +   a 3 log T w 100   +   a 4 ( T w 100 ) 2   +   S · ( b 1   +   b 2 T w 100   +   b 3 ( T w 100 ) 2 )
f = exp ( 101.325 P ( V m     V ideal ) RT w )
log p H 2 O P = 24.4543     67.4509 ( 100 T w )     4.8489 ln ( T w 100 )     0.000544 S
Johnson [20] developed an algorithm from another approach that corrected for the effects of temperature and salinity on the molecular and thermodynamic properties of the water, its solutes and the specified gas, based on a ideal gas behaviour. His formulae were developed from the compilation by Sander [27] (available on the web since 1999) of the kHcp for nearly all gases in the atmosphere at 25 °C (298.15 K) and 0 ppt. Then, Equation (6) converted the kHcp to kH at a given temperature and 0 ppt salinity. The term −ΔsolnH/R reflected the temperature dependence of solubility, having a value of 2400 for CO2, 1700 for CH4 and 2600 for N2O. The correction to a given salinity (Equation (7)) relied on the empirical Setschenow constants (KS = θ∙logVb) reporting the effect of electrolytes salting-out gases proportionally to their liquid molar volume at boiling point (Vb). The Vb was estimated using the additive Schroeder method, whereas θ was estimated from Equation (8) using a provisional kH# = 0.0409/kHcp.
k H , 0 =   12.1866 P     T w     k H , cp     e Δ Soln H R ( 1 / T w 1 / 298.15 )
k H = k H , 0 · 10 K s S
θ = 7.33532 10   4   +   3.39615 10 5 log ( k H # )     2.40888 10 6 log ( k H # ) 2   +   1.57114 10 7 log ( k H # ) 3

2.2. Transfer Velocity

The simpler transfer velocity formulations rely on the wind velocity 10 m above the sea surface (u10) under atmospherically neutral conditions. Among then, the formulation by Wanninkhof [29] (henceforth also referred to as ‘Wan92’) became the standard used in ESM and satellite data processing (Equation (9a,b)). It includes the Schmidt number of the water (Scw), related to viscosity, raised to an exponent sensitive to the assumed model of surface renewal. For sparingly soluble gases or for CO2 in alkaline waters, and under low winds, the transfer velocity is enhanced by the hydration reaction (bT), which is sensitive to temperature.
k w = ( b T   +   0.31 u 10 2 ) ( Sc w 660 ) 0.5
b T = 2.5 ( 0.5246   +   0.0162 T w   +   0.000499 T w 2 )
Other simple empirical formulations based only on u10 [30,31], or also accounting for bottom–generated turbulence [32], used data collected in estuaries under low wind conditions. However, modelling the coastal ocean at finer resolutions requires an enhanced representation of the multitude of processes involved. Hence, we updated the framework by Vieira et al. [21], with the kw being decomposed into its shear produced turbulence (kwind + kcurrent) and bubbles from whitecapping (kbublle) forcings [32,41,42,44]. Scw was determined from temperature and salinity following Johnson [20]:
k w = ( b T   +   k bubble   +   k wind   +   k current ) ( 600 / Sc w ) 0.5
The formulation by Zhao et al. [43], merged kwind within the kbubble (Equation (11a)), being this term dependent on the wave breaking parameter (RB given by Equation (11b)). The u* was the friction velocity i.e., the velocity of wind dragging on the sea surface, and fp the peak angular frequency of the wind-waves (in rad∙s−1). We estimated the kinematic viscosity of air (υa) from Johnson [20].
k bubble   +   k wind = 0.1315 R B 0.6322
R B = u 2 2 π f p ν a
Later, Zhao and Xie [47] developed a more elaborate algorithm. While the output was still given by Equation (11a), Rb was now estimated from Equation (12a). RHU was a turbulent Reynolds number (Equation (12b)) accounting for the sea state by using the significant wave height (Hs). The wave age (β) was estimated from Equation (12c) using the gravitational acceleration constant (g). While Zhao and Xie [47] tested their algorithm with the drag coefficient (CD) estimated from several empirical formulations, here we estimated CD from its definition (CD = u*2/u102). This was only possible because u* was estimated from u10 using the wind log-linear profile as presented below.
R B =   4.79 C D β 2 3 R HU
R HU = u 10 H s ν a
β = 2.56 ( g H s u 10 2 ) 3 5
The formulation by Woolf [44] split the estimation of kwind and kbubble. The kwind, taken from Jähne et al. [66], related to the transfer mediated by the turbulence generated by wind drag (Equation (13)). The kbubble = 850 W related to the transfer mediated by the bubbles generated by breaking waves. The whitecap cover was estimated from W = 4.02 ∙ 10−7 RH0.96, and the dimensionless form of the Reynolds number for wind waves estimated from RH = u*∙Hsa.
k wind = 1.57 · 10 4 u
The friction velocity (u*) was estimated from the wind log-linear profile (WLLP: Equation (14)) accounting for wind speed at height z (uz), atmospheric stability of the surface boundary layer (through ψm) and sea-surface roughness (through the roughness length z0). Here, κ is von Kármàn’s constant. Historically, the WLLP originated from Monin–Obukhov similarity theory [67].
u = u z κ ln ( z )     ln ( z 0 )   +   ψ m ( z , z 0 , L )
Roughness length (z0) is the theoretical minimal height (most often sub-millimetric) at which wind speed averages zero. It is dependent on surface roughness and often used as its index. It is more difficult to determine over water than over land as there is a strong bidirectional interaction between wind and sea surface roughness. Taylor and Yelland [68] proposed a dimensionless z0 dependency from the wave field, increasing with the wave slope (Equation (15)) as estimated from the significant wave height (Hs) and peak wave period (Lp). Due to the bidirectional nature of the z0 and u* relation, we also tested an iterative solution (iWLP) where Equation (15) was used as a first guess for the z0 and Equation (14) for its subsequent u*. A second iteration re-estimated z0 from the COARE 3.0 [69] adaptation of the Taylor and Yelland [68] formulation, whereby a term for smooth flow was added (Equation (16)), and u* was estimated back again from Equation (14). Applying four iterations were enough for an excellent convergence of the full data array. The coefficients proposed by Taylor and Yelland [68] yielded implausibly high z0 when applied to steep young waves originated by a storm over the central Mediterranean Sea and by strong Mistral blowing off-shore of southern France, corresponding to 2.85% of our observations. These extreme z0 led to absurdly high u* and kw, and thus we imposed a maximum roughness length of 0.01 m. This value roughly matches the maximum expected from the application of alternative formulations for the estimation of z0 from the wave field [70,71,72,73,74].
z 0 H s = 1200 ( H s L p ) 4.5
z 0 = 1200 H s ( H s L p ) 4.5   +   0.11 ν a u
Atmospheric stability characterized the tendency of the surface boundary layer (SBL) to be well mixed (unstable SBL with ψm < 0) or stratified (stable SBL with ψm > 0). The ψm was inferred from the bulk Richardson number (Rib), weighting the air vertical heat gradient and kinetic energy (Equation (17)). Hence, it required the air virtual potential temperature (Tv) and wind speed at two heights (z) within the SBL. Alternatively, the Rib was directly estimated from the air potential temperature (Tp) i.e., replacing Tv and, thus, neglecting humidity [75].
Ri b = g Δ T v Δ z T v u z 2
Tv was estimated from Tp and humidity, following Equation (18) [39] or Equation (19) [67]. Each equation requires different humidity inputs, namely the observed specific humidity (q), relative humidity (hr), liquid water mixing ratio (rL) and water vapour mixing ratio at saturation (rsat):
T v = T p ( 1 + 0.61 q )
T v = T p ( 1 + 0.61 q r sat h r + r L )
Tp and rsat were estimated from the air temperature (Tka) and pressure (P) [67]:
T p = T k , a ( 1000 1013.25 P ) 0.284
r sat = 0.622 e sat 101.32501 P e sat
ln ( e sat ) = ln 0.61078 + 17.2694 T a T k , a 35.86
At surface height z ≈ 0 m, the wind velocity must be set to the theoretical u0 = 0, temperature given by the SST [34,39,69] with rectification for cool-skin and warm-layer effects [34,35,36,37,38], and humidity set to the saturation level at P0 and T0 using q = rsat/(1 + rsat) [69].The Rib was used to estimate the length L from the Monin–Obukhov’s similarity theory, a discontinuous exponential function tending to ±∞ when Rib tends to ±0 and tending to ±0 when Rib tends to ±∞. Rib and L were used to estimate ψm following Stull [67] or Lee [75] algorithms. The transfer velocity of moderately soluble gases, besides taking into consideration the molecular crossing of the water-side surface layer, should also take into consideration the molecular crossing of the air-side surface layer [17,18,19,20,21,22]. We compared values from a single layer model and two-layer “thin film” model [17], the later estimating the air-side transfer velocity (ka) from the COARE formulation as in Equation (23) [46]. CD is the drag coefficient and Sca the Schmidt number of air, which were determined for a given temperature and salinity following Johnson [20].
k a = u 13.3 Sc a 1 / 2   +   CD 1 / 2     5   +   log ( Sc a ) 2 κ

2.3. Validation with Field Data

The algorithms presented in Section 2.2 were applied to observed oceanographic and atmospheric variables (Level 2 data). These are derived geophysical variables at the same resolution and location as the instrument source data (Level 1). For these estimates we used the FuGas 2.2 single processing ensemble. The forecasted CO2 transfer velocities were compared to those inferred from kw = F/Δpgas, where F were the CO2 vertical fluxes measured by eddy–covariance, and the Δpgas were the differences between measured atmospheric and oceanic CO2 concentrations. The field sampling occurred from 22 May 2014 to 26 May 2014 using the atmospheric tower at Östergarnsholm in the Baltic Sea (57°27′ N, 18°59′ E), the Submersible Autonomous Moored Instrument (SAMI–CO2) 1 km away and the Directional Waverider (DWR) 3.5 km away, both south-eastward of the tower [9,76]. The atmospheric measurements were done at 12 m heights. The DWR measured temperatures at 0.5 m depth, taken as representative for the sea-surface. Salinity was obtained from the Asko mooring data provided by the Baltic In–Situ Near–Real–Time Observations available in Copernicus Marine catalogue.
The eddy–covariance (E–C) measurements were subject to quality control. We only used the fluxes for which the wind direction set the SAMI–CO2 and DWR in the footprint of the atmospheric tower (100° < wind direction < 160°). We rejected periods for which there were abrupt changes in wind speed or direction as these represented strong violations of the required homogeneous conditions [77]. An extra filter was applied to remove measurements in which wave age denoted large dominance of swell conditions (β > 2), since they have been identified as biasing measurements of transfer velocity [47,48]. When the ratio is >1 the waves move faster than the wind, indicating the presence of swell [47,48]. Using only the observations passing the quality control check, the air–water CO2 fluxes measured by eddy–covariance were smoothed over 30 min bins and corrected according to the Webb–Pearman–Leuning (WPL) method [78].

2.4. Application within Geoscientific Models

The transfer velocities of CO2, CH4 and N2O were estimated from atmospheric and oceanographic re-analysed (Level 4) data retrieved for the European shores (Mediterranean and North Atlantic to the tip of Scotland) from the 24 May 2014 at 06 h to the 27 May 2014 at 00 h. All variables were interpolated to the same 0.09° grid (roughly 11 km at Europe’s latitudes) and 1 h time steps. This resulted in a data set with 17 variables × 41,776 locations × 66 time instances. For a faster processing of the ≈2 Gb of Level 4 data, we used the multiple processing ensemble of the FuGas 2.3.
The atmospheric variables were retrieved from the standard operational application of the WRF by Meteodata.cz, with model output at 9 km horizontal and 1 h time resolutions. The vertical thickness of the WRF layers varied with space and time. The atmospheric variables where retrieved from the two lowest levels within the SBL. These levels refer to the horizontal layers in the atmospheric model, and are not to be mistaken by the data processing level. Over the ocean, these levels occurred roughly at 0 m and 12 m heights. The second level provided the wind velocity, temperature, pressure and humidity at z meters above sea surface. At surface height z ≈ 0 m, the wind velocity was set to the theoretical u0 = 0, and temperature was given by the SST without rectification for cool-skin and warm-layer effects due to the lack of some required variables. Yet, these effects tend to compensate each other [34,37,38]. Air pressure was given by the WRF at the lower level (≈0 m) and humidity was set to the saturation level using P0 and T0.
The sea surface temperature (SST, here also named Tw) and salinity (S) were estimated by the NEMO modelling system provided in the MyOcean catalogue with 1/12° and 1 day resolutions. The WW3 wave field data for the Mediterranean Sea was supplied by the Istituto Nazionale di Geofisica e Vulcanologia (INGV) using the WW3–NEMO modelling system at 0.0625° and 1 h resolutions [79], and for the North Atlantic by Windguru at roughly 0.5° and 3 h resolutions. The variables included significant wave height (Hs) and peak frequency (fp) for wind sea i.e., disregarding swell. A few aspects did not correspond to the ideal data format for atmosphere–ocean coupling and required further calculations: (i) The peak wave length (Lp) was estimated from the peak frequency assuming the deep-water approximation: Lp = 2πg/fp2; (ii) the Windguru data did not provide a wind sea component where (and when) the wind was too low. For these missing cases, the lowest Hs were attributed and Lp was simulated everywhere else; and (iii) the Windguru and the INGV data overlapped along the Iberian shores, in which case the INGV was given a 2:1 weight over the Windguru data.
SST and S were used to estimate how the solubility algorithms could lead to different estimates of greenhouse gases dissolved below the ocean surface. First, the solubility was converted to the Cw/Ca form corresponding to 1/kH. Then, the differences in solubility between the two algorithms, as estimated from Equation (24), led to the differences in the estimates of greenhouse gases dissolved per cubic meter immediately below the ocean surface (Equation (25)).
Δ s = 1 k H [ Sar 13 ] 1 k H [ Joh 10 ] = C w [ Sar 13 ] C w [ Joh 10 ] C a
Δ g × m 3 = Δ s ρ gas P 101325 M gas 10 6 RT
For these calculations, time variability was disregarded assuming the average Δs over the 66 h time interval. The ∆Cw was converted to mass units using the molar weights of the specific gas (Mgas), respectively 44.01 for CO2, 16.043 for CH4 and 44.013 for N2O. The Ca was estimated from the atmospheric pressure (P in atm) and the partial pressure (ρgas) of CO2, CH4 or N2O, 390 ppm, 1.75 ppm and 0.325 ppm respectively [80], assuming that they were approximately uniform all over the atmospheric surface boundary layer (SBL). The ρgas was divided by 106 to re-scale from ppm to dimensionless. Using the ideal gas law, we divided by R and Ta (assumed Ta = Tw at the sea-surface), and multiplied by 101325 to convert atm back to Pascal.

3. Results and Discussion

3.1. Solubility Estimates

The solubility formulae were compared for the range of environmental conditions commonly found in nature: water temperature (Tw) ranged from −2 °C to 30 °C at 1 °C intervals and a salinity (S) ranged from 0 ppt to 38 ppt at 1 ppt intervals. In general, both formulations estimated similar solubilities despite their distinct chemistry backgrounds, as demonstrated by their ratio kH[20]/kH[2] close to 1 (Figure 1a–c). Nevertheless, differences from 10% to −5% for CO2, 11% to 5% for CH4, and 7% to −5% for N2O were observed in freshwaters. Determining which formulation provides more accurate estimates of greenhouse gas (GHG) solubilities in freshwater is fundamental for Earth System Modelling, because rivers and freshwater reservoirs are important sources of GHG to the atmosphere, even if only seasonably [7,8,9,10,11,12,13,14,16]. In marine waters, our simulations found differences of up to 12% for CO2, 9% for CH4 and 7% for N2O, suggesting that for marine waters is also important to determine which formulation provides more accurate estimates of GHG solubilities. Afterwards, both formulae were applied to the conditions observed at the European coastal ocean during the experiment (the Level 4 SST, S and P). The surface water temperature changed significantly and there were large fresh water inflows from the Black Sea and the Baltic Sea (Figure 2 and Supplementary Video S1). Solubility estimates were averaged over a 66 h time interval and compared at each location by their ratio (Figure 1d–f). Higher divergences reached +4.5% for CO2 solubility in cooler marine waters, +5.65% for CH4 solubility in warmer saltier waters, and −2.1% for N2O solubility in brackish waters. Hence, the algorithms diverge by as much as 0.045 mol∙mol−1 of CO2, 0.0015 mol∙mol−1 of CH4 and 0.012 mol∙mol−1 of N2O (i.e., mol of gas in the ocean surface per mol of gas in the atmosphere) in some of the most sensitive situations for Earth System Modelling and satellite data processing, namely: (i) the cooler marine sub-polar waters where the solubility pump traps greenhouse gases and carries them to the deep ocean [2], and (ii) the warmer waters at the coastal ocean and seas, which have regularly been observed having greenhouse gases and DMS dissolved in concentrations highly unbalanced with those of the atmosphere [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16].
The discrepancies between solubilities estimated by either algorithm led to differences in the estimates of GHG dissolved per cubic meter immediately below the ocean surface (Figure 3). The major differences occurred in the regions with the cooler waters, the warmer waters and the brackish waters. These biases in the estimated total amount of greenhouse gases for the European coastal ocean’s surface during late May 2014, may be an indicator of higher global biases. Hence, similar comparisons with data from other locations and/or periods are required. Performing these comparisons is facilitated the by the FuGas framework.

3.2. Transfer Velocity Estimates from Field Data

The transfer velocity (kw) scaled with u10 (r = 0.91) and with the inverse of the wave age (r = −0.88) i.e., the younger the waves, the larger the kw (Figure 4). This fits the general theory that under low–moderate wind speeds, the kw is set by surface renewal and micro-scale wave breaking (i.e., by kwind), increasing with steeper younger waves [28,51,65], and only under higher wind speeds is kw set by bubbles from breaking waves (i.e., kbubble) [41,42,43,44,45,46,47,48,49,50]. The preponderance of surface renewal and micro-scale wave breaking at the coastal ocean turn it more susceptible to additional factors driving kw. As an example, under the low–moderate wind speeds often verified at the coastal ocean, the hydration reaction enhances kw [19] and surfactants attenuate surface renewal [28,58,59,65,66].
The SBL was stable i.e., 0.016 < Rib < 0.093. The sea surface was slightly to moderately rough (z0 < 0.49 mm). Under these conditions, the kw estimated from the E–C measurements were approximate to the kw estimated by both empirical u10-based parameterizations and extended physically based algorithms (Figure 4 and Figure 5). The small kw fluctuations estimated by empirical u10–based parameterizations were the sole consequence of changes in water viscosity (as estimated by the Scw) driven by changes in water temperature. On the other hand, the physically based algorithms, by accounting for the effects of sea state, breaking waves and atmospheric stability, showed a greater adaptability to local conditions (Figure 5). As an example, depending on the sea state, the formulations by Woolf [44] or Zhao and Xie [47] estimated kw similar to that estimated by the Carini et al. [30] parameterization for the Parker river estuary, or similar to that estimated by the Wanninkhof [29] parameterization for the global ocean. Consequently, the root mean square deviation (RMSD) between the Woolf [44] formulation (henceforth ‘W05va’) and the E–C transfer velocities was 6.34 cm∙h−1, whereas between the Wan92 and the E–C transfer velocities it was 8.4 cm∙h−1. On relative terms, the physically based algorithm showed 75% the deviation of the empirical u10-based algorithm. Nevertheless, the physically based algorithms were still unable to accurately fit E–C observations (Figure 4 and Figure 5).
With the optional coupling of the sea state and wind velocity forcings (through the iWLP joint estimate of z0 and u*), the FuGas enables an enhanced representation of local conditions that was demonstrated to be fundamental for the estimation of transfer velocities in fetch-limited situations [50,55,57]. This coupling agrees with observations that kw scales with the turbulent kinetic energy dissipation at the sea surface (ε), and that ε is better reflected by the sea state [54,56]. Accordingly, the COARE 3.0 included the wave state in the estimation of the roughness parameters essential for the transfer of mass, heat and momentum [69]. We highlight the fact that the roughness Reynolds number is given by R = u*∙z0a i.e., relying precisely on the u* and z0 that are jointly estimated by the FuGas from the SBL and wave field.
The estimation of kw as wind asymptotically tends to zero is generally considered uninteresting for open ocean formulations, although the equatorial global ocean consistently presents such low winds. The bulk of the empirical u10–based parameterizations undervalues the accuracy of the kw estimated under the lowest winds and make it asymptotically tending to zero as well. The parameterizations by Mackay and Yeun [87], McGillis et al. [84], Wanninkhof et al. [19] and Johnson [20] propose eliminating this misrepresentation of reality by making kw asymptotically tend to a constant. Arguably, the accurate estimation of kw under low wind speeds becomes more relevant when studying the coastal ocean and inland waters, and extended physically-based algorithms are better suited for the task. The COARE algorithm added a gustiness term stabilizing the kw in finite velocities under lighter winds [69]. Clayson et al. [51], on the other hand, replaced the gustiness term by a capillary wave parameterization. The FuGas achieved this through the iterative estimation of the wind log–linear profile (the iWLP: Equations (15) and (16)). Under the lowest winds but as long as there were waves, the smooth flow term in the iWLP imposed finite u*, and consequently, finite transfer velocities (Figure 5).

3.3. Transfer Velocity Estimates from Level 4 Data

We performed simulations for the European costal ocean to compare between the ESM standard (the ‘Wan92’) and our extended physically-based alternatives. We show the comparison with the iWLP implemented with the kwind formulation by Jähne et al. [66] and the kbubble formulation by Woolf [44] using the kinematic viscosity of air (the W05va). Their kw estimates diverged mostly under unstable SBL, very rough or very smooth sea-surfaces, or higher friction velocities (Figure 6). The details of the simulations and the differences between kw estimates are presented hereafter.

3.3.1. Atmospheric Stability

Strong winds occurred along the European shores from the 24–26 May 2014. Besides, the air was unusually cold for the season and colder than the sea surface (Supplementary Video S1). Where the wind blew lighter, the rise of the warmer air heated by the sea surface generated turbulent eddies that enhanced mixing within the SBL and, consequently, also enhanced u* and kw. These unstable conditions, identified by Rib < 0, L tending to −0 and ψm < 0 (Supplementary Video S2), occurred more frequently and intensively nearby land masses and often associated to cooler continental breezes blowing off-shore. Its correct simulation required the estimation of the Rib, L and ψm from the algorithms by Grachev and Fairall [39] and Stull [67] that account for humidity considering saturation at 0 m heights. The Rib estimates neglecting humidity, following Lee [75], yielded biased estimates of the SBL conditions as a consequence of biased estimates of the virtual potential temperature. Stull [67] (page 9), Grachev and Fairall [39] and Fairall et al. [69] had already highlighted the importance of accounting for humidity.

3.3.2. Roughness Length

The sea surface agitation was very heterogeneous, particularly at the coastal ocean where it attained both the highest and the lowest estimated roughness lengths (the z0 in Supplementary Video S3). There, the steeper waves, a consequence of shorter fetches, should extract more momentum from the atmosphere under similar u10 conditions [54,55,56,58]. Thus, the rougher coastal ocean surfaces were expected to have relatively more turbulent layers through which gases were transferred at higher rates. The comprehensive formulations simulated this by increasing u* (and consequently kwind) with z0 under similar uz i.e., similar winds generate more drag when blowing over rougher sea surfaces. Closer to the shores, and very often within the Mediterranean, lighter wind blew over smoother sea surfaces. Under these conditions, the iWLP estimated much higher z0 than the WLLP (Supplementary Video S4), demonstrating that the smooth flow was a fundamental driver for the z0 under calmer weather. The increase in z0 led to significantly higher u*, often 1.5 times higher and sometimes more, anticipating a significant impact on the local kwind estimates. These results demonstrated how the iWLP stabilized the kw in finite velocities as wind speed approached zero, and how this becomes a relevant matter when simulating the coastal ocean.

3.3.3. Transfer Velocities

In some rare situations the algorithms estimated unreasonably high kw despite the z0 bounds imposed in the software. To avoid this bias, a 200 cm∙h−1 ceiling was imposed on all kw estimates. The higher transfer velocities were associated to storms, both on the Atlantic side and on the Mediterranean side (Figure 7 and Supplementary Video S5). The root mean square deviation (RMSD) was estimated between the empirical u10-based Wan92 formula and the physically based W05va formula, averaged over time for each location. Results showed that both formulations diverged the most, in absolute terms, under storms and under strong winds near shore (Figure 7). On the other hand, the normalized root mean square deviation (NRMSD)—i.e., the deviation standardized to the Wan92 prediction—showed that, in relative terms, these formulations diverged the most under low and moderate winds and sea states near shore (Figure 7). The physically-based W05va formulation was also used to compare between the kwind and kbubble components of the CO2 transfer velocity. The kbubble term was always lower than the kwind term, and at best 90% of kwind in the fetch-unlimited Atlantic or in the Mediterranean storms (Supplementary Video S5). Over the full domain, the kbubble term was on average 12% and on median 15.9% of the kwind term. Whatever the greenhouse gas, the differences were negligible between estimating kw using the single-layer or the double-layer schemes (Supplementary Video S5). Nevertheless, it is worth noting that the bigger differences were found under Mediterranean storm conditions.

4. Conclusions

The sea state and atmospheric stability were determinant factors for the gas transfer velocities at the coastal ocean, irrespective of current algorithms still needing calibration and validation. Our findings corroborate recent conclusions that the wave state is fundamental for accurate estimates of gas transfer velocities at the fetch-limited coastal ocean [50,57]. Furthermore, our conclusions also agree with Jackson et al. [88] and Shuiqing and Dongliang [49] when inferring about the open ocean. In fact, the COAREG (the gas-dedicated version of the COARE) by Jackson et al. [88] was applied to data from the GasEx cruises on the Northern Atlantic, Equatorial Pacific and Southern Ocean, with their results also suggesting that the sea state and atmospheric stability are important additional factors driving kw in particular situations. Our application of the FuGas to the coastal ocean and the Jackson et al. ([88] Figure 11) application of the COAREG to the global ocean even agree in the magnitude of the enhancement of the transfer velocities, when compared to u10–based empirical parameterizations. Regarding the coarse resolution of Earth System Models, our FuGas 2.2 application demonstrated that physically–based algorithms diverge even more from empirical u10–based parameterizations when applied over the coastal ocean with finer spatial and temporal resolutions. Jackson et al. [88] also implemented a temporal fine resolution (3 h), having observed a wide discrepancy with the kw estimated from monthly means, and concluded that the kw estimated from monthly averaged factors is highly biased due to linearization of highly non-linear processes. The chemical enhancement of the transfer velocity is another factor that may become important under the low wind speeds often verified at the coastal ocean. Yet, it still needs to be implemented in the FuGas and tested regarding its functional form, the choice of gases, and the required environmental conditions [28,89,90].

Supplementary Materials

Software, data and table related to this article are available online at https://www.mdpi.com/2073-4433/9/8/310/s1. The latest FuGas version and videos is available at http://www.maretec.org/en/models/fugas.

Author Contributions

Conceptualization, V.M.N.C.S.V.; Data curation, V.M.N.C.S.V., P.J. and E.C.; Formal analysis, V.M.N.C.S.V.; Funding acquisition, M.M.; Investigation, V.M.N.C.S.V.; Methodology, V.M.N.C.S.V.; Project administration, M.M.; Resources, P.J., E.C. and M.M.; Software, V.M.N.C.S.V.; Supervision, V.M.N.C.S.V.; Validation, V.M.N.C.S.V.; Visualization, V.M.N.C.S.V.; Writing—original draft, V.M.N.C.S.V.; Writing—review and editing, V.M.N.C.S.V., P.J. and M.M.

Funding

This research was funded by ERDF Funds of the Competitiveness Factors Operational Programme—COMPETE and by national funds from the FCT—Foundation for Science and Technology project UID/EEA/50009/2013. The APC was funded by ERDF Funds of the Competitiveness Factors Operational Programme—COMPETE and by national funds from the FCT—Foundation for Science and Technology project UID/EEA/50009/2013.

Acknowledgments

We are thankful to Erik Sahlée for providing the data from Uppsala’s University field station in the Baltic, to Heidi Pettersson for providing the data from the Warider Buoy, and to windguru.cz for the support providing the L4 wave data.

Conflicts of Interest

The authors declare no competing financial interests.

References

  1. Takahashi, T.; Sutherland, S.C.; Sweeney, C.; Poisson, A.; Metzl, N.; Tilbrook, B.; Bates, N.; Wanninkhof, R.; Feely, R.A.; Sabine, C.; et al. Global sea-air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects. Deep Sea Res. Part II Top. Stud. Oceanogr. 2002, 49, 1601–1622. [Google Scholar] [CrossRef]
  2. Sarmiento, J.L.; Gruber, N. Ocean Biogeochemical Dynamics; Princeton University Press: Princeton, NJ, USA, 2013; pp. 73–100. [Google Scholar]
  3. Gypens, N.; Borges, A.V. Increase in dimethylsulfide (DMS) emissions due to eutrophication of coastal waters offsets their reduction due to ocean acidification. Front. Mar. Sci. 2015, 1. [Google Scholar] [CrossRef] [Green Version]
  4. Menzo, Z.M.; Elliott, S.; Hartin, C.A.; Hoffman, F.M.; Wang, S. Climate change impacts on natural sulfur production: Ocean acidification and community shifts. Atmosphere 2018, 9, 167. [Google Scholar] [CrossRef]
  5. Borges, A.V. Do we have enough pieces of the jigsaw to integrate CO2 fluxes in the coastal ocean? Estuaries 2005, 28, 3–27. [Google Scholar] [CrossRef]
  6. Cole, J.J.; Caraco, N.F. Atmospheric exchange of carbon dioxide in a low-wind oligotrophic lake measured by the addition of SF6. Limnol. Oceanogr. 1998, 43, 647–656. [Google Scholar] [CrossRef]
  7. Cole, J.J.; Caraco, N.F. Carbon in catchments: Connecting terrestrial carbon losses with aquatic metabolism. Mar. Freshw. Res. 2001, 52, 101–110. [Google Scholar] [CrossRef]
  8. Borges, A.V.; Delille, B.; Frankignoulle, M. Budgeting sinks and sources of CO2 in the coastal ocean: Diversity of ecosystems counts. Geophys. Res. Lett. 2005, 32, L14601. [Google Scholar] [CrossRef]
  9. Rutgersson, A.; Norman, M.; Schneider, B.; Petterson, H.; Sahlée, E. The annual cycle of carbon dioxide and parameters influencing the air–sea carbon exchange in the Baltic Proper. J. Mar. Syst. 2008, 74, 381–394. [Google Scholar] [CrossRef]
  10. Harley, J.F.; Carvalho, L.; Dudley, B.; Heal, K.V.; Rees, R.M.; Skiba, U. Spatial and seasonal fluxes of the greenhouse gases N2O, CO2 and CH4 in a UK macrotidal estuary. Estuar. Coast. Shelf Sci. 2015, 153, 62–73. [Google Scholar] [CrossRef] [Green Version]
  11. Lansø, A.S.; Bendtsen, J.; Christensen, J.H.; Sørensen, L.L.; Chen, H.; Meijer, H.A.J.; Geels, C. Sensitivity of the air–sea CO2 exchange in the Baltic Sea and Danish inner waters to atmospheric short-term variability. Biogeosciences 2015, 12, 2753–2772. [Google Scholar] [CrossRef]
  12. Dutta, M.K.; Mukherjee, R.; Jana, T.K.; Mukhopadhyay, S.K. Biogeochemical dynamics of exogenous methane in an estuary associated to a mangrove biosphere; The Sundarbans. NE coast of India. Mar. Chem. 2015, 170, 1–10. [Google Scholar] [CrossRef]
  13. Nevison, C.D.; Weiss, R.F.; Erickson, D.J., III. Global oceanic emissions of nitrous oxide. J. Geophys. Res. 1995, 100, 15809–15820. [Google Scholar] [CrossRef]
  14. Nevison, C.D.; Lueker, T.J.; Weiss, R.F. Quantifying the nitrous oxide source from coastal upwelling. Glob. Biogeochem. Cycles 2004, 18, GB1018. [Google Scholar] [CrossRef]
  15. Barnes, J.; Upstill-Goddard, R.C. N2O seasonal distributions and air-sea exchange in UK estuaries: Implications for the tropospheric N2O source from European coastal waters. J. Geophys. Res. 2011, 116, G01006. [Google Scholar] [CrossRef]
  16. Smith, S.V.; Hollibaugh, J.T. Coastal metabolism and the oceanic organic carbon balance. Rev. Geophys. 1993, 31, 75–89. [Google Scholar] [CrossRef]
  17. Liss, P.S.; Slater, P.G. Flux of gases across the air-sea interface. Nature 1974, 247, 181–184. [Google Scholar] [CrossRef]
  18. Upstill-Goddard, R.C. Air-sea gas exchange in the coastal zone. Estuar. Coast. Shelf Sci. 2006, 70, 388–404. [Google Scholar] [CrossRef]
  19. Wanninkhof, R.; Asher, W.E.; Ho, D.T.; Sweeney, C.S.; McGillis, W.R. Advances in quantifying air-sea gas exchange and environmental forcing. Annu. Rev. Mar. Sci. 2009, 1, 213–244. [Google Scholar] [CrossRef] [PubMed]
  20. Johnson, M.T. A numerical scheme to calculate temperature and salinity dependent air-water transfer velocities for any gas. Ocean Sci. 2010, 6, 913–932. [Google Scholar] [CrossRef]
  21. Vieira, V.M.N.C.S.; Martins, F.; Silva, J.; Santos, R. Numerical tools to estimate the flux of a gas across the air–water interface and assess the heterogeneity of its forcing functions. Ocean Sci. 2013, 9, 355–375. [Google Scholar] [CrossRef]
  22. Vieira, V.; Sahlée, E.; Jurus, P.; Clementi, E.; Pettersson, H.; Mateus, M. Improving estimates of atmosphere-ocean greenhouse gas fluxes for Earth-System modelling. In ESA SP, Proceedings of the Living Planet Symposium, Prague, Czeck Republic, 9–13 May 2016; European Space Agency: Paris, France, 2016; Volume 740. [Google Scholar]
  23. Weiss, R.F. Carbon dioxide in water and seawater: The solubility of a non-ideal gas. Mar. Chem. 1974, 2, 203–215. [Google Scholar] [CrossRef]
  24. Weiss, R.F.; Price, B.A. Nitrous oxide solubility in water and seawater. Mar. Chem. 1980, 8, 347–359. [Google Scholar] [CrossRef]
  25. Mackay, D.; Shiu, W.Y. A critical review of Henry’s law constants of environmental interest. J. Phys. Chem. Ref. Data 1981, 10, 1175–1197. [Google Scholar] [CrossRef]
  26. Copin-Montegut, C. A new formula for the effect of temperature on the partial pressure of CO2 in seawater. Mar. Chem. 1988, 25, 29–37. [Google Scholar] [CrossRef]
  27. Sander, R. Compilation of Henry’s law constants (Version 4.0) for water as solvent. Atmos. Chem. Phys. 2015, 15, 4399–4981. [Google Scholar] [CrossRef]
  28. Frew, N.M.; Bock, E.J.; Schimpf, U.; Hara, T.; Haubecker, H.; Edson, J.B.; McGillis, W.R.; Nelson, R.K.; McKenna, S.P.; Uz, B.M.; et al. Air-sea gas transfer: Its dependence on wind stress, small-scale roughness, and surface films. J. Geophys. Res. 2004, 109, C08S17. [Google Scholar] [CrossRef]
  29. Wanninkhof, R. Relationship between wind speed and gas exchange over the ocean. J. Geophys. Res. 1992, 97, 7373–7382. [Google Scholar] [CrossRef]
  30. Carini, S.; Weston, N.; Hopkinson, C.; Tucker, J.; Giblin, A.; Vallino, J. Gas exchange rates in the Parker River estuary, Massachusetts. Biol. Bull. 1996, 191, 333–334. [Google Scholar] [CrossRef] [PubMed]
  31. Raymond, P.A.; Cole, J.J. Gas exchange in rivers and estuaries: Choosing a gas transfer velocity. Estuaries 2001, 24, 312–317. [Google Scholar] [CrossRef]
  32. Borges, A.V.; Vanderborght, J.P.; Schiettecatte, L.S.; Gazeau, F.; Ferron-Smith, S.; Delille, B.; Frankignoulle, M. Variability of the Gas Transfer Velocity of CO2 in a Macrotidal Estuary (the Scheldt). Estuaries 2004, 27, 593–603. [Google Scholar] [CrossRef]
  33. Soloviev, A.V.; Schlüssel, P. Parameterization of the temperature difference across the cool skin of the ocean and of the air–ocean gas transfer on the basis of modelling surface renewal. J. Phys. Oceanogr. 1994, 24, 1339–1346. [Google Scholar] [CrossRef]
  34. Brunke, M.A.; Zeng, X.; Misra, V.; Beljaars, A. Integration of a prognostic sea surface skin temperature scheme into weather and climate models. J. Geophys. Res. 2008, 113, D21117. [Google Scholar] [CrossRef]
  35. Soloviev, A.; Schlüssel, P. Evolution of cool skin and direct air-sea gas transfer coefficient during daytime. Bound. Layer Meteorol. 1996, 77, 45–68. [Google Scholar] [CrossRef]
  36. Soloviev, A.V.; Lukas, R. Observation of spatial variability of diurnal thermocline and rain-formed halocline in the western Pacific warm pool. J. Phys. Oceanogr. 1996, 26, 2529–2538. [Google Scholar] [CrossRef]
  37. Fairall, C.W.; Bradley, E.F.; Godfrey, J.S.; Wick, G.A.; Edson, J.B.; Young, G.S. Cool-skin and warm-layer effects on sea surface temperature. J. Geophys. Res. 1996, 101, 1295–1308. [Google Scholar] [CrossRef]
  38. Zeng, X.; Beljaars, A. A prognostic scheme of sea surface skin temperature for modelling and data assimilation. Geophys. Res. Lett. 2005, 32, L14605. [Google Scholar] [CrossRef]
  39. Grachev, A.A.; Fairall, C.W. Dependence of the Monin–Obukhov Stability Parameter on the Bulk Richardson Number over the Ocean. J. Appl. Meteorol. 1997, 36, 406–414. [Google Scholar] [CrossRef]
  40. Fairall, C.W.; Hare, J.E.; Edson, J.B.; McGillis, W. Parameterization and micrometerological measurements of air-sea gas transfer. Bound. Layer Meteorol. 2000, 96, 63–105. [Google Scholar] [CrossRef]
  41. Asher, W.E.; Farley, P.J. Phase-Doppler anemometer measurement of bubble concentrations in laboratory-simulated breaking waves. J. Geophys. Res. 1995, 100, 7045–7056. [Google Scholar] [CrossRef]
  42. Woolf, D.K. Bubbles and their role in gas exchange. In The Sea Surface and Global Change; Duce, R.A., Liss, P.S., Eds.; Cambridge University Press: New York, NY, USA, 1997; pp. 173–205. [Google Scholar]
  43. Zhao, D.; Toba, Y.; Suzuki, Y.; 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 B Chem. Phys. Meteorol. 2003, 55, 478–487. [Google Scholar] [CrossRef]
  44. Woolf, D.K. Parameterization of gas transfer velocities and sea state-dependent wave breaking. Tellus B Chem. Phys. Meteorol. 2005, 57, 87–94. [Google Scholar] [CrossRef]
  45. Zhang, W.; Perrie, W.; Vagle, S. Impacts of winter storms on air-sea gas exchange. Geophys. Res. Lett. 2006, 33, L14803. [Google Scholar] [CrossRef]
  46. Jeffery, C.; Robinson, I.; Woolf, D. Tuning a physically-based model of the air-sea gas transfer velocity. Ocean Model. 2010, 31, 28–35. [Google Scholar] [CrossRef]
  47. Zhao, D.; Xie, L. A practical bi-parameter formula of gas transfer velocity depending on wave states. J. Oceanogr. 2010, 66, 663–671. [Google Scholar] [CrossRef]
  48. Goddijn-Murphy, L.; Woolf, D.K.; Callaghan, A.H. Parameterizations and algorithms for oceanic whitecap coverage. J. Phys. Oceanogr. 2011, 41, 742–756. [Google Scholar] [CrossRef]
  49. Shuiqing, L.; Dongliang, Z. Gas transfer velocity in the presence of wave breaking. Tellus B Chem. Phys. Meteorol. 2016, 68, 27034. [Google Scholar] [CrossRef] [Green Version]
  50. Gutiérrez-Loza, L.; Ocampo-Torres, F.J.; Nava, H.G. The effect of breaking waves on CO2 air-sea fluxes in the coastal zone. Bound. Layer Meteorol. 2018. [Google Scholar] [CrossRef]
  51. Clayson, C.A.; Fairall, C.W.; Curry, J.A. Evaluation of turbulent fluxes at the ocean surface using surface renewal theory. J. Geophys. Res. 1996, 101, 28503–28513. [Google Scholar] [CrossRef]
  52. Bock, E.J.; Edson, J.B.; Frew, N.M.; Hara, T.; Haussecker, H.; Jähne, B.; McGillis, W.R.; McKenna, S.P.; Nelson, R.K.; Schimpf, U.; et al. Overview of the CoOP experiment: Physical and chemical measurements parameterizing air-sea heat exchange. In Gas Transfer at Water Surfaces; Donelan, M.A., Drennan, W.M., Saltzman, E.S., Wanninkhof, R., Eds.; AGU: Washington, DC, USA, 2002; pp. 39–44. [Google Scholar]
  53. Zappa, C.J.; Asher, W.E.; Jessup, A.T.; Klinke, J.; Long, S.R. Microbreaking and the enhancement of air-water transfer velocity. J. Geophys. Res. Oceans 2004, 109. [Google Scholar] [CrossRef] [Green Version]
  54. Soloviev, A.; Donelan, M.; Graber, H.; Haus, B.; Schlüssel, P. An approach to estimation of near-surface turbulence and CO2 transfer velocity from remote sensing data. J. Mar. Syst. 2007, 66, 182–194. [Google Scholar] [CrossRef]
  55. Turney, D.E.; Banerjee, S. Air-water gas transfer and near-surface motions. J. Fluid Mech. 2013, 733, 588–624. [Google Scholar] [CrossRef]
  56. Wang, B.; Liao, Q.; Fillingham, J.; Bootsma, H.A. On the coefficients of small eddy and surface divergence models for the air-water gas transfer velocity. J. Geophys. Res. Oceans 2015, 120. [Google Scholar] [CrossRef]
  57. Bigdeli, A.; Hara, T.; Loose, B.; Nguyen, A.T. Wave attenuation and gas exchange velocity in marginal sea ice zone. J. Geophys. Res. Oceans 2018, 123. [Google Scholar] [CrossRef]
  58. Frew, N.M.; Nelson, R.K.; McGillis, W.R.; Edson, J.B.; Bock, E.J.; Hara, T. Spatial variations in surface micro-layer surfactants and their role in modulating air-sea exchange. In Gas Transfer at Water Surfaces; Donelan, M.A., Drennan, W.M., Saltzman, E.S., Wanninkhof, R., Eds.; AGU: Washington, DC, USA, 2002; pp. 153–159. [Google Scholar]
  59. McKenna, S.P.; McGillis, W.R. The role of free-surface turbulence and surfactants in air-water gas transfer. Int. J. Heat Mass Transf. 2004, 47, 539–553. [Google Scholar] [CrossRef]
  60. Pereira, R.; Schneider-Zapp, K.; Upstill-Goddard, R.C. Surfactant control of gas transfer velocity along an offshore coastal transect: Results from a laboratory gas exchange tank. Biogeosciences 2016, 13, 3981–3989. [Google Scholar] [CrossRef]
  61. Ho, D.T.; Bliven, L.F.; Wanninkhof, R. The effect of rain on air-water gas exchange. Tellus B Chem. Phys. Meteorol. 1997, 49, 149–158. [Google Scholar] [CrossRef]
  62. Ho, D.T.; Asher, W.E.; Bliven, L.F.; Schlosser, P.; Gordan, E.L. On mechanisms of rain-induced air-water gas exchange. J. Geophys. Res. 2000, 105, 24045–24057. [Google Scholar] [CrossRef] [Green Version]
  63. Frost, T.; Upstill-Goddard, R.C. Meteorological controls of gas exchange at a small English lake. Limnol. Oceanogr. 2002, 47, 1165–1174. [Google Scholar] [CrossRef] [Green Version]
  64. Zappa, C.J.; Ho, D.T.; McGillis, W.R.; Banner, M.L.; Dacey, J.W.H.; Bliven, L.F.; Ma, B.; Nystuen, J. Rain-induced turbulence and air-sea gas transfer. J. Geophys. Res. Oceans 2009, 114. [Google Scholar] [CrossRef] [Green Version]
  65. Zappa, C.J.; Asher, W.E.; Jessup, A.T.; Klinke, J.; Long, S.R. Effect of microscale wave breaking on air-water gas transfer. In Gas Transfer at Water Surfaces; Donelan, M.A., Drennan, W.M., Saltzman, E.S., Wanninkhof, R., Eds.; AGU: Washington, DC, USA, 2002; pp. 23–29. [Google Scholar]
  66. Jähne, B.; Munnich, K.O.; Bosinger, R.; Dutzi, A.; Huber, W.; Libner, P. On the parameters influencing air-water gas exchange. J. Geophys. Res. 1987, 92, 1937–1949. [Google Scholar] [CrossRef]
  67. Stull, R.B. An Introduction to Boundary Layer Meteorology; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1988; pp. 151–195. [Google Scholar]
  68. Taylor, P.K.; Yelland, M.J. The dependence of sea surface roughness on the height and steepness of the waves. J. Phys. Oceanogr. 2001, 31, 572–590. [Google Scholar] [CrossRef]
  69. Fairall, C.W.; Bradley, E.F.; Hare, J.E.; Grachev, A.A.; Edson, J.B. Bulk parameterization of air-sea fluxes: Updates and verification for the COARE algorithm. J. Clim. 2003, 16, 571–591. [Google Scholar] [CrossRef]
  70. Smith, S.D.; Anderson, R.J.; Oost, W.A.; Kraan, C.; Maat, N.; DeCosmo, J.; Katsaros, K.B.; Davidson, K.L.; Bumke, K.; Hasse, L.; et al. Sea surface wind stress and drag coefficients: The HEXOS results. Bound. Layer Meteorol. 1992, 60, 109–142. [Google Scholar] [CrossRef]
  71. Johnson, H.K.; Hoejstrup, J.; Vested, H.J.; Larsen, S.E. Dependence of sea surface roughness on wind waves. J. Phys. Oceanogr. 1998, 28, 1702–1716. [Google Scholar] [CrossRef]
  72. Lange, B.; Johnson, H.K.; Larsen, S.; Hojstrup, J.; Kofoed-Hansen, H.; Yelland, M.J. On detection of a wave age dependency for the sea surface roughness. J. Phys. Oceanogr. 2004, 34, 1441–1458. [Google Scholar] [CrossRef]
  73. Hwang, P. Drag coefficient, dynamic roughness and reference wind speed. J. Oceanogr. 2005, 61, 399–413. [Google Scholar] [CrossRef]
  74. Pan, Y.; Sha, W.; Zhu, S.; Ge, S. A new parameterization for sea surface aerodynamic roughness. Prog. Nat. Sci. 2008, 18, 1365–1373. [Google Scholar] [CrossRef]
  75. Lee, H.N. Improvement of surface flux calculations in the atmospheric surface layer. J. Appl. Meteorol. 1997, 36, 1416–1423. [Google Scholar] [CrossRef]
  76. Högström, U.; Sahlée, E.; Drennan, W.M.; Kahma, K.K.; Smedman, A.-S.; Johansson, C.; Pettersson, H.; Rutgersson, A.; Tuomi, L.; Zhang, F.; et al. Momentum fluxes and wind gradients in the marine boundary layer—A multi-platform study. Boreal Environ. Res. 2008, 13, 475–502. [Google Scholar]
  77. Donelan, M.A.; Drennan, W.M.; Katsaros, K.B. The air-sea momentum flux in conditions of wind sea and swell. J. Phys. Oceanogr. 1997, 27, 2087–2099. [Google Scholar] [CrossRef]
  78. Webb, E.K.; Pearman, G.I.; Leuning, R. Correction of flux measurements for density effects due to heat and water vapour transfer. Quart. J. R. Meteorol. Soc. 1980, 106, 85–100. [Google Scholar] [CrossRef]
  79. Clementi, E.; Oddo, P.; Korres, G.; Drudi, M.; Pinardi, N. Coupled wave-ocean modelling system in the Mediterranean Sea. In Proceedings of the 13th International Workshop on Wave Hindcasting, Banff, AB, Canada, 27 October–1 November 2013. [Google Scholar]
  80. EPA United States Environmental Protection Agency. Climate Change Indicators in the United States—Atmospheric Concentrations of Greenhouse Gases. Available online: https://www.epa.gov/climate-indicators/climate-change-indicators-atmospheric-concentrations-greenhouse-gases (accessed on 29 April 2018).
  81. Jacobs, C.M.J.; Kohsiek, W.; Oost, W.A. Air–sea fluxes and transfer velocity of CO2 over the North Sea: Results from ASGAMAGE. Tellus B Chem. Phys. Meteorol. 1999, 51, 629–641. [Google Scholar] [CrossRef]
  82. Wanninkhof, R.; McGillis, W.R. A cubic relationship between air-sea CO2 exchange and wind speed. Geophys. Res. Lett. 1999, 26, 1889–1892. [Google Scholar] [CrossRef]
  83. Nightingale, P.D.; Malin, G.; Law, C.S.; Watson, A.J.; Liss, P.S.; Liddicoat, M.I.; Boutin, J.; Upstill-Goddard, R.C. In Situ evaluation of air-sea gas exchange parameterizations using novel conservative and volatile tracers. Glob. Biogeochem. Cycles 2000, 14, 373–387. [Google Scholar] [CrossRef]
  84. McGillis, W.R.; Edson, J.B.; Ware, J.D.; Dacey, J.W.H.; Hare, J.E.; Fairall, C.W.; Wanninkhof, R. Carbon dioxide flux techniques performed during GasEx-98. Mar. Chem. 2001, 75, 267–280. [Google Scholar] [CrossRef]
  85. Weiss, A.; Kuss, J.; Peters, G.; Schneider, B. Evaluating transfer velocity-wind speed relationship using a long-term series of direct eddy correlation CO2 flux measurements. J. Mar. Syst. 2007, 66, 130–139. [Google Scholar] [CrossRef]
  86. Ho, D.T.; Law, C.S.; Smith, M.J.; Schlosser, P.; Harvey, M.; Hill, P. Measurements of air-sea gas exchange at high wind speeds in the Southern Ocean: Implications for global parameterizations. Geophys. Res. Lett. 2006, 33, L16611. [Google Scholar] [CrossRef]
  87. Mackay, D.; Yeun, A.T.K. Mass transfer coefficient correlations for volatilization of organic solutes from water. Environ. Sci. Technol. 1983, 17, 211–217. [Google Scholar] [CrossRef] [PubMed]
  88. Jackson, D.L.; Wick, G.A.; Hare, J.E. A comparison of satellite-derived carbon dioxide transfer velocities from a physically based model with GasEx cruise observations. J. Geophys. Res. 2012, 117, G00F13. [Google Scholar] [CrossRef]
  89. Wanninkhof, R.; Knox, M. Chemical enhancement of CO2 exchange in natural waters. Limnol. Oceanogr. 1996, 41, 689–697. [Google Scholar] [CrossRef]
  90. Bade, D.L.; Cole, J.J. Impact of chemically enhanced diffusion on dissolved inorganic carbon stable isotopes in a fertilized lake. J. Geophys. Res. 2006, 111, C01014. [Google Scholar] [CrossRef]
Figure 1. Match–mismatch between estimated solubilities of greenhouse gases. Colour scale is quotient between kH estimated from Johnson [20] and from Sarmiento and Gruber [2]. Solubilities tested for the range of values observable at the Earth’s surface (ac) and at the European coastal ocean from the 24 May 2014 at 06 h to the 27 May 2014 at 00 h (df) with each pixel corresponding to a square sided 11 km and averaged over the 66 h survey period.
Figure 1. Match–mismatch between estimated solubilities of greenhouse gases. Colour scale is quotient between kH estimated from Johnson [20] and from Sarmiento and Gruber [2]. Solubilities tested for the range of values observable at the Earth’s surface (ac) and at the European coastal ocean from the 24 May 2014 at 06 h to the 27 May 2014 at 00 h (df) with each pixel corresponding to a square sided 11 km and averaged over the 66 h survey period.
Atmosphere 09 00310 g001
Figure 2. Sea surface temperature (SST), salinity (S) and air pressure (P) at the European coastal ocean from the 24 May 2014 at 06 h to the 27 May 2014 at 00 h. Each pixel corresponds to a square sided 11 km and averaged over the 66 h survey period.
Figure 2. Sea surface temperature (SST), salinity (S) and air pressure (P) at the European coastal ocean from the 24 May 2014 at 06 h to the 27 May 2014 at 00 h. Each pixel corresponds to a square sided 11 km and averaged over the 66 h survey period.
Atmosphere 09 00310 g002
Figure 3. Bias in the gas mass estimated by each algorithm (Δg∙m−3). Comparing algorithms by Johnson [20] to compilation by Sarmiento and Gruber [2]. Colourscale: Δg m−3.
Figure 3. Bias in the gas mass estimated by each algorithm (Δg∙m−3). Comparing algorithms by Johnson [20] to compilation by Sarmiento and Gruber [2]. Colourscale: Δg m−3.
Atmosphere 09 00310 g003
Figure 4. Transfer velocities (kw) estimated for the Baltic experiment from eddy–covariance observations (E–C). The Jac99 (Jacobs et al. [81]) and WMG99 (Wanninkhof and McGillis [82]) are the empirical u10–based parameterizations estimating higher and lower kw within this wind range.
Figure 4. Transfer velocities (kw) estimated for the Baltic experiment from eddy–covariance observations (E–C). The Jac99 (Jacobs et al. [81]) and WMG99 (Wanninkhof and McGillis [82]) are the empirical u10–based parameterizations estimating higher and lower kw within this wind range.
Atmosphere 09 00310 g004
Figure 5. Comparison among transfer velocities (kw) estimated for the Baltic experiment from empirical u10-based parameterizations and from extended physically based formulations. Were tested the empirical u10-based parameterizations by Wanninkhof et al. [19], Wanninkhof [29], Carini et al. [30], Raymond and Cole [31], Borges et al. [32], Jacobs et al. [81], Wanninkhof and McGillis [82], Nightingale et al. [83], McGillis et al. [84], Weiss et al. [85], Ho et al. [86], Were tested the physically-based parameterizations by Zhao et al. [43], Woolf [44] and Zhao and Xie [47] conjugated with the iWLP algorithm specific of the FuGas.
Figure 5. Comparison among transfer velocities (kw) estimated for the Baltic experiment from empirical u10-based parameterizations and from extended physically based formulations. Were tested the empirical u10-based parameterizations by Wanninkhof et al. [19], Wanninkhof [29], Carini et al. [30], Raymond and Cole [31], Borges et al. [32], Jacobs et al. [81], Wanninkhof and McGillis [82], Nightingale et al. [83], McGillis et al. [84], Weiss et al. [85], Ho et al. [86], Were tested the physically-based parameterizations by Zhao et al. [43], Woolf [44] and Zhao and Xie [47] conjugated with the iWLP algorithm specific of the FuGas.
Atmosphere 09 00310 g005
Figure 6. Comparing the formulation by Wanninkhof [29] (Wan92) to the formulation by Woolf [44] (W05va) using modelled data about the European coastal ocean. Differences in the estimates of transfer velocity (kw) as functions of the bulk Richardson number (Rib), roughness length (z0) and friction velocity (u*). The 1:1 line (in black) represent equal kw estimates.
Figure 6. Comparing the formulation by Wanninkhof [29] (Wan92) to the formulation by Woolf [44] (W05va) using modelled data about the European coastal ocean. Differences in the estimates of transfer velocity (kw) as functions of the bulk Richardson number (Rib), roughness length (z0) and friction velocity (u*). The 1:1 line (in black) represent equal kw estimates.
Atmosphere 09 00310 g006
Figure 7. Comparing transfer velocity (kw) algorithms using modelled data: (Wan92) transfer velocity of CO2 estimated from Wanninkhof [29] and averaged over the 66 h; (RMSD) Root mean square deviation between estimating the transfer velocity following Wanninkhof [29] or Woolf [44] conjugated with the iterative wind log-linear profile. (NRMSD) Root mean square deviation normalized to the Wan92 value at each time and location.
Figure 7. Comparing transfer velocity (kw) algorithms using modelled data: (Wan92) transfer velocity of CO2 estimated from Wanninkhof [29] and averaged over the 66 h; (RMSD) Root mean square deviation between estimating the transfer velocity following Wanninkhof [29] or Woolf [44] conjugated with the iterative wind log-linear profile. (NRMSD) Root mean square deviation normalized to the Wan92 value at each time and location.
Atmosphere 09 00310 g007

Share and Cite

MDPI and ACS Style

Vieira, V.M.N.C.S.; Jurus, P.; Clementi, E.; Mateus, M. The FuGas 2.3 Framework for Atmosphere–Ocean Coupling: Comparing Algorithms for the Estimation of Solubilities and Gas Fluxes. Atmosphere 2018, 9, 310. https://doi.org/10.3390/atmos9080310

AMA Style

Vieira VMNCS, Jurus P, Clementi E, Mateus M. The FuGas 2.3 Framework for Atmosphere–Ocean Coupling: Comparing Algorithms for the Estimation of Solubilities and Gas Fluxes. Atmosphere. 2018; 9(8):310. https://doi.org/10.3390/atmos9080310

Chicago/Turabian Style

Vieira, Vasco M. N. C. S., Pavel Jurus, Emanuela Clementi, and Marcos Mateus. 2018. "The FuGas 2.3 Framework for Atmosphere–Ocean Coupling: Comparing Algorithms for the Estimation of Solubilities and Gas Fluxes" Atmosphere 9, no. 8: 310. https://doi.org/10.3390/atmos9080310

APA Style

Vieira, V. M. N. C. S., Jurus, P., Clementi, E., & Mateus, M. (2018). The FuGas 2.3 Framework for Atmosphere–Ocean Coupling: Comparing Algorithms for the Estimation of Solubilities and Gas Fluxes. Atmosphere, 9(8), 310. https://doi.org/10.3390/atmos9080310

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop