Next Article in Journal
Automated Mapping of Ms 7.0 Jiuzhaigou Earthquake (China) Post-Disaster Landslides Based on High-Resolution UAV Imagery
Next Article in Special Issue
Air-Sea Interactions over Eddies in the Brazil-Malvinas Confluence
Previous Article in Journal
Image Enhancement Driven by Object Characteristics and Dense Feature Reuse Network for Ship Target Detection in Remote Sensing Imagery
Previous Article in Special Issue
Twenty-Seven Years of Scatterometer Surface Wind Analysis over Eastern Boundary Upwelling Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mesoscale Temporal Wind Variability Biases Global Air–Sea Gas Transfer Velocity of CO2 and Other Slightly Soluble Gases

1
Division of Earth and Ocean Sciences, Nicholas School of the Environment, Duke University, Durham, NC 27708, USA
2
College of Oceanography, Hohai University, Nanjing 210098, China
3
Nicholas School of the Environment, Box 90328, Duke University, Durham, NC 27708, USA
4
Department of Civil and Environmental Engineering, Duke University, Durham, NC 27708, USA
5
CNRS, Univ Brest, IRD, Ifremer, LEMAR, F-29280 Plouzané, France
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(7), 1328; https://doi.org/10.3390/rs13071328
Submission received: 7 February 2021 / Revised: 26 March 2021 / Accepted: 29 March 2021 / Published: 31 March 2021
(This article belongs to the Special Issue Remote Sensing of Air-Sea Fluxes)

Abstract

:
The significance of the water-side gas transfer velocity for air–sea CO2 gas exchange (k) and its non-linear dependence on wind speed ( U ) is well accepted. What remains a subject of inquiry are biases associated with the form of the non-linear relation linking k to U (hereafter labeled as f( U ), where f(.) stands for an arbitrary function of U ), the distributional properties of U (treated as a random variable) along with other external factors influencing k, and the time-averaging period used to determine k from U . To address the latter issue, a Taylor series expansion is applied to separate f( U ) into a term derived from time-averaging wind speed (labeled as U , where .   indicates averaging over a monthly time scale) as currently employed in climate models and additive bias corrections that vary with the statistics of U . The method was explored for nine widely used f( U ) parameterizations based on remotely-sensed 6-hourly global wind products at 10 m above the sea-surface. The bias in k of monthly estimates compared to the reference 6-hourly product was shown to be mainly associated with wind variability captured by the standard deviation σ U around U or, more preferably, a dimensionless coefficient of variation I u = σ U / U . The proposed correction outperforms previous methodologies that adjusted k when using U only. An unexpected outcome was that upon setting I u 2 = 0.15 to correct biases when using monthly wind speed averages, the new model produced superior results at the global and regional scale compared to prior correction methodologies. Finally, an equation relating I u 2 to the time-averaging interval (spanning from 6 h to a month) is presented to enable other sub-monthly averaging periods to be used. While the focus here is on CO2, the theoretical tactic employed can be applied to other slightly soluble gases. As monthly and climatological wind data are often used in climate models for gas transfer estimates, the proposed approach provides a robust scheme that can be readily implemented in current climate models.

1. Introduction

Describing air–sea flux of long-lived greenhouse gases such as carbon dioxide (CO2) is of significance for assessing the global carbon cycle and its relation to climate. In climate models, the water-side air–sea flux (F, mol m−2 y−1) is commonly determined using a bulk expression
F = K 0   k   Δ pCO 2
where k is the gas transfer velocity (cm h−1), K0 is the gas solubility (mol L−1 atm−1) in water that is a function of sea surface temperature (SST) and salinity, and ∆pCO2 is the difference in partial pressure of pCO2 between water and air (atm). The dominant factors determining k are governed by a number of physical processes primarily, but not exclusively, associated with wind speed U . For this reason, k is operationally parameterized as a non-linear function of U set at a reference height of 10 m. For comparison purposes, the general formulations (common ones listed in Table 1) take the form of
k = ( Sc / 660 ) 1 / 2   f ( U )
where f( U ) (cm h−1) is a non-linear function of wind speed U , also known as the gas transfer velocity k660 normalized to the dimensionless molecular Schmidt number (Sc) for CO2 in seawater at 20 °C (Sc = 660). The function f( U ) may be quadratic, cubic, or even a higher-order polynomial, and Sc ( 1) is the ratio of the kinematic viscosity (m2 s−1) and the molecular diffusion coefficient (m2 s−1) of CO2 or other gases in seawater. For this reason, Equation (2) is routinely used for slightly soluble gases. The f( U ) can also be derived using turbulent transport theories [1,2,3,4,5], bubbles [6,7], and wave-breaking mechanics [8,9,10]. However, f( U ) cannot be viewed as linking k to an instantaneous U at a point; rather, f( U ) must emerge as an approximation to macroscopic equations derived by averaging gas transfer over space and time scales (analogous to a closure model for turbulent fluxes in Reynolds-averaged Navier–Stokes equations). The spatial scales must be much larger than the largest eddy or wave length impacting gas exchange, whereas the time scales must be sufficiently long to accommodate the effects of turbulent fluctuations (i.e., ensemble of many eddy-turnover times) or wave formation and subsequent breaking, but short enough to resolve mesoscale variations in U . This interval is commensurate with hourly time scales and coincides with time scales associated with the well-known spectral gap in the atmosphere [11]. Fourier power spectra of wind time series sampled from fractions of seconds (turbulent scales) to years support the occurrence of a “gap” in the squared Fourier amplitudes separating mesoscales (longer than few hours) from micro-scales (smaller than minutes). This gap forms the basis of separating U into a micro-scale contribution whose effects on k are to be averaged out and captured by f( U ) and a meso-scale or longer (i.e., larger than hours) contribution [12].
Based on gas transfer velocity parameterizations, modelling and observation-based estimates of the global oceanic CO2 sink vary significantly from −1.18 to −3.1 Pg C yr−1 (negative referring to net flux of CO2 into the ocean) [13,14,15,16,17,18]. The range in these estimates reflects different time periods and uncertainties. Uncertainties result from using various data products, methodological uncertainties in k parameterizations, the relative sparsity of CO2 data coverage in time and space, and thermal and haline effects [19,20,21,22,23,24,25].
Additionally, temporal averaging of wind data substantially contributes to uncertainties in global F estimates due to the non-linearity in f( U ) [26,27,28] and frames the scope of the work here. As an example, if wind blows half the time at a speed of 4 m s−1 and the other half of the time at 16 m s−1 (solid points on the curves), k values estimated from the mean wind speed of 10 m s−1 are biased low by 11.2 and 30.6 cm h−1 relative to the true k (circles on dash lines) for the quadratic and the cubic relations, respectively. Quadratic and cubic equations are taken from [29,30], respectively (Figure 1). Long-term averaged (monthly or longer) wind speeds underestimate gas exchange by 25% and by 50% for quadratic and cubic f( U ), respectively [31]. Such known biases can be handled by: (i) using wind speeds with short temporal intervals (i.e., 6 h) or (ii) applying correction factors when averaging over longer intervals (e.g., monthly) to mitigate these expected biases [26,27,28]. Because using short-term wind speeds (e.g., 6-hourly wind products) globally to evaluate f( U ) is computationally expensive now and in the foreseeable future, bias-corrected methodologies are gaining attention. However, reported biases in k are still pronounced even after applying current correction factors, thereby motivating the development of other approaches. The time is ripe to begin exploring such bias corrections to existing gas exchange formulations given the availability of satellite-based wind products at 6-hourly temporal resolution.
In this study, we compare various published corrections and propose a new method that we test globally and regionally for any function f( U ). The findings here apply for any slightly soluble gas for which its f( U ) is known. For a more accurate correction, wind variance should also be supplied to correct monthly k. In the absence of such information, the work here suggests a constant squared coefficient of variation can be used (Iu2 = 0.15). The manuscript is organized as follows: the datasets and data processing, review of current correction methods, and the proposed new correction method are presented in Section 2. In Section 3, this new method is applied to 29 years of data to obtain the corrected gas transfer velocity and is evaluated by comparing the newly corrected k to results from earlier models and studies. A summary and concluding remarks are presented in Section 4.

2. Data and Methods

2.1. Data and Data Processing

The U and sea surface temperature (SST) data from 1990 to 2018 are used to compute globally averaged k. The 6-hourly and averaged monthly cross-calibrated multiplatform CCMP V2.0 wind data at the 0.25° × 0.25° grid were obtained from the Remote Sensing Systems described elsewhere [42]. The gridded CCMP V2.0 wind products are produced from a combination of satellite (Version-7 remote sensing system radiometer wind speeds, QuikSCAT, and ASCAT scatterometer wind vector wind), moored buoy, and model wind data. The high-resolution CCMP captures wind variability well and is relatively bias free compared to in situ estimates [43,44]. Therefore, the CCMP data is often used to parameterize the gas transfer velocity coefficient and in estimating k [25,33,43,45]. The 6-hourly and monthly SST data at the spatial resolution of 0.25° × 0.25° were derived from the European Centre for Medium-Range Weather Forecasts (ECMWF) fifth generation ERA5 reanalysis product described elsewhere [46]. The wind speeds and SST are linearly interpolated onto a spatial resolution of 0.5° × 0.5°. Within this grid, the statistics of U are assumed to be planar homogeneous. Additionally, both U and SST data were linearly interpolated to 5° × 5° to evaluate the effect of spatial resolution on k estimates. Nine commonly used f( U ) parameterizations (Table 1) were applied to estimate globally averaged gas transfer velocity for CO2. The Sc for CO2 is a function of SST and is determined using a standard formulation [33].

2.2. Review of Prior Correction Methods for the Time-Average Bias

Correction factors are routinely applied to account for biases in k estimates associated with time-averaging of wind speeds. One commonly used method is based on a Reynolds decomposition into a well-defined mean U (e.g., monthly) and fluctuations U   (on the scales of hours) around this average so that U = U + U with U = 0. When f ( U ) = aU 2 , f ( U ) = a U 2 and f ( U ) = aU 2 .   Therefore, the f ( U ) can be evaluated as
f ( U ) = a ( U +   U ) 2 = a ( U 2 + 2 U U +   U 2 ) =   a U 2 [ 1 + ( σ U U ) 2 ] = f ( U ) [ 1 + ( σ U / U ) 2 ]
where σ U is the standard deviation, I u = ( σ U / U ) is as before the coefficient of variation, and the sought correction (as a ratio or bias) can be expressed as
f ( U ) f ( U ) = ( 1 + I u 2 )   and   f ( U ) f ( U ) = a σ U 2
Clearly, this correction depends on the non-linearity of f ( U )   [32,47]. With available 6-hourly wind speed ( U 6   hour ) data, other widely used multiplier corrections [26,27,28,31,40,45,48,49,50,51,52] for the quadratic formulation is expressed in this form as
R 2 = U 6   hour 2 / U 6   hour 2
and for the cubic formulation
R 3 = U 6   hour 3 / U 6   hour 3
These R2 and R3 corrections can be obtained empirically or derived analytically when assuming the distributional properties of U for meso-scale (and longer) variations [26,27,28,31,40,45,48,49,50,51,52]. A Rayleigh distribution, which is commonly used in the evaluation of R2 and R3 [26] arises when the magnitude of the wind velocity is analyzed in two dimensions (usually in the plane parallel to the water surface) (see Text S1 in Supplementary). Assuming that each component is uncorrelated and normally distributed with equal variance in each of the two directions (i.e., planar homogeneous air flow), the overall wind vector magnitude is characterized by a Rayleigh distribution (i.e., a special form of Chi-squared). These R2 and R3 corrections are also simplified using zonally averaged profiles [26,27]. Globally and regionally, the R2 ranges from 1.12 to 1.26 whereas the R3 ranges from 1.35 to 2.17 [26,28,48,53].

2.3. Proposed Correction Based on Taylor Series Expansions

For a wind-only related formulation, f( U ) U n (n   >   1) and upon space-time averaging yields U n   U n . To assess biases arising from setting f ( U ) = f   ( U ) , a Taylor series expansion of any f( U ) form around the space-time averaged value f ( U )   are now introduced and given by
f ( U ) = f ( U ) + df dU | U ( U U ) +   1 2 ! d 2 f dU 2 | U ( U U ) 2   + 1 3 ! d 3 f dU 3 | U ( U U ) 3 +
Applying the space-time averaging operation   . term by term in Equation (7) yields
f ( U ) = f ( U ) + df dU | U U U +   1 2 d 2 f dU 2 | U ( U U ) 2        +   1 6 d 3 f dU 3 | U ( U U ) 3 +
This expression can be arranged as
    f ( U )   = f ( U ) +   k b , with   k b =   f ( U )   f ( U )   Term   1 =   1 2 d 2 f dU 2 | U σ U 2 +   1 6 d 3 f dU 3 | U ( U U ) 3 Term   2  
where k b is the sought bias and σ U 2 is the variance in wind speed around U . In this expression, all derivatives of f( U ) are presumed to be known (e.g., Table 1 or other physics-based formulation) and are being evaluated at U . When n = 1, all the derivative terms are identically zero and k b = 0 . General expression for a quadratic relation (n = 2) f( U ) can be expressed as
f ( U ) =   aU 2
applying Equation (9) to f(U) results in
f ( U ) = a U 2 + k b ,   with   k b = a   σ U 2
suggesting an additive correction (i.e., bias) that only varies with σ U 2 . This expression is consistent with [47], though the approach taken here is more general. For the cubic relations (n = 3),
f ( U ) =   a   U 3 + bU 2 + dU + e
and this results in
f ( U ) = a   U 3 + b U 2 + d U + e + k b ,   with   k b = 3 a U σ U 2 + b σ U 2 + a ( U U ) 3
In this case, the skewness of U as estimated with the Fisher–Pearson coefficient of skewness ( Sk = ( U U ) 3 / σ U 3 ) is also required. In summary, the newly proposed corrections and the three prior correction methods listed in Table 2 are used to adjust for the bias.

3. Results

3.1. Bias in k Induced by Averaging of Wind Data

Globally averaged k for CO2 at various temporal and spatial resolutions were assessed for the parameterizations of f( U ), listed in Table 1. The k computed from maximum spatial (0.5°   ×   0.5°) and temporal resolution (6-hourly) products were used as a reference to illustrate the deviation of k in percentage (Figure 2 and Table S1). As expected, the monthly k underestimates k for all parameterizations. The absolute biases induced by time-averaging wind speed (~10–28% range) are more significant for cubic relations (Serial NO. (2), (4), (5), (7), and (8)) than for quadratic relations (Serial NO. (1), (3), (6), and (9)), with a comparable magnitude in biases at both spatial resolutions of 0.5°   ×   0.5° and 5°   ×   5°. In contrast, uncertainties due to differential spatial resolutions are negligible (less than 1%) at the temporal resolutions of both 6 h and a month (Figure 2). For this reason, we only focus hereafter on the uncertainty induced by differential temporal resolution of wind speed data, though the method can be applied to any type of averaging [28]. As expected, k substantially varies with the choice of f( U ) being used (Figure S2). Undoubtedly, the mechanisms constraining gas transfer velocity must be explored, but this issue is beyond the scope of the present work.

3.2. Assessment of the “Bias Correction Model”

This new “bias correction model” was applied to all the nine parameterizations of k, as shown in Table 1. Measured bias in f(U) (term 1 in Equation (9), Figure S3) and bias kb from the new model (term 2 in Equation (9), Figure S4) were estimated. In term 1, the 6-hourly space-time product was used to evaluate f ( U ) and the monthly space-time product was used to evaluate f ( U ) . Overall, the proposed model reproduces the bias between 6-hourly k and monthly k (Figure 3). Spatially, the differences in the first term and the second term are negligible in quadratic parameterizations ((1), (3), (6), and (9) in Figure S5). In contrast, for the cubic relations such as the k parameterizations of (2), (4), (5), (7), and (8), the differences are small, lower than 0.6 cm h−1 (Figure 3). The higher values in the mid and high latitude of the northern hemisphere (Figure S5) might be associated with large variability in wind speed within a month (Figure S1) due to the occurrence of synoptic high wind events in these regions [54].

3.3. Comparison of Correction Methods

The proposed new bias correction was compared to common correction methodologies. We tested the new corrections presented in Equations (11) and (13) based on two methods. In method 1, the correction was estimated using the annual mean grid-by-grid kb (i.e., with grid-by-grid σ U 2 , U and ( U U ) 3 ). In method 2, the correction was the averaged kb (i.e., with averaged σ u -related terms) (Table 2). As σ u and U both increase in time, the squared coefficient of variation I u 2 =   ( σ U / U ) 2 is a more “conserved” parameter in time with a slowly decreasing (insignificant) trend of 0.002 dec−1 and an average of 0.15 (Figure 4). Thus, for a constant   I u 2 = 0.15 in method 2, the kb of quadratic relations (Equation (11)) can be arranged as
k b = a   U 2 ( σ u U ) 2 = a   U 2 I u 2 = 0.15 a U 2
To simplify further, only the first term in Equation (12) is taken for cubic expressions, so that
f ( U ) = a U 3 + k b
Evidently, kb is also a function of I u 2 for cubic expressions given as
k b = 3 a   U 3 ( σ u U ) 2 = 3 a   U 3 I u 2 = 0.45 a U 3
With constant I u 2 = 0.15 , Equations (14) and (16) can be used to approximate unbiased f( U ), and k considering differences in Sc for CO2 and other gases (see [33] for Sc of other gases). The same f( U ) parameterization for all slightly soluble gases may not be realistic for gases with differing solubilities [55,56,57], but this inquiry is better kept for the future.
The newly proposed corrections were compared to three commonly used methods (listed in Table 2) when adjusting wind variability-induced bias in k for wind speeds sampled at a resolution of 0.5°   ×   0.5° spatially and temporally monthly. Taking CO2 as an example, globally averaged corrected k values were calculated for all the nine parameterizations, and biases in corrected k are estimated in reference to the 6-hourly k.
As expected, deviations of the corrected k from the reference 6-hourly k are significantly reduced using the new approaches (Figure 5). The largest absolute biases yielded by methods 1 and 2 are only 0.6% and 4.46%, respectively. In contrast, the range of absolute biases generated by the other three methods vary from 3.5% to 28% (Table S2). It is worth noting that the bias induced by the time averaging of SST values is negligible (small difference in 6-hourly and corrected k after applying correction only associated with wind variance) due to the small variance of SST (Figure 6) and the smaller effect of SST on k for all parameterizations (Table 3). From the perspective of zonal distribution, the magnitude of the corrected k using method 1 agrees well with the 6-hourly k (Figure 7a). Though the simpler method 2 appears to overestimate k at mid and high latitude regions (within 0–30°N and 30°–60°N), the overall corrected k is generally consistent with the 6-hourly k (especially in the Southern hemisphere), with a smaller root-mean-square errors (RMSE) between corrected k and the 6-hourly k compared to methods 3, 4, and 5 (Figure 7b).
Moreover, the relation between the averaging period Δ t   (spanning from 6 h to a month) versus I u 2 (Figure 8), needed to infer biases in k computed using wind speeds at a coarse temporal resolution, was empirically derived and given as
I u 2 = 0.18 Δ t 0.22 + 0.237 .
Therefore, corrected f( U ) in non-linear wind-only parameterizations for any gases at any temporal resolution from 6-hourly to monthly (or longer, not show here) can now be estimated by substituting   I u 2 in Equations (11) and (14) for quadratic relations, and in Equations (15) and (16) for cubic relations.

3.4. Study Limitation

As shown by the Taylor expansion (Equation (9)), the bias due to time-averaging is a function of the variance in wind speed. Hence, a corollary question to explore is whether the 6-hourly wind speed used as a reference here introduces biases because of missing variances at sub-hourly time scales. There may be a fraction of energy in the scale of 6 h to minutes or seconds commensurate with the time scales, over which the turbulent or wave action impact k but are presumably captured by the functional form of f( U ).
The spectra of 6-hourly global wind is calculated to show the spectral energy distribution from mesoscale to decadal scale (Figure 9). As estimated, the variance in this range is σ d 2   = 0.1 m2 s–2. A f–3 scaling in the spectrum from multi-day to a 12-h range appears to be supported here and implies an enstrophy cascade in quasi-geostrophic flow [58,59]. If the spectrum is extrapolated from 12 h to turbulence scale (seconds) via a Kolmogorov’s –5/3 power law, the “missing variance” in this range is σ m 2 0.001 m2 s–2, which can be ignored. Extrapolations to finer scales via a f–3 scaling would result in an even smaller missing variance. To be clear, this does not imply that the air turbulent time scales (on the order of 10 s and smaller) are minor. The energy contents of these time scales are quite large but are captured by the non-linearity in f( U ) as noted earlier. To illustrate this point, a variance in turbulence scale σ t 2   of 1 m2 s–2 can be estimated from a turbulence similarity relation   σ t 2 = ( 2.3   u * )2 and u *   = U CD1/2, where u * is the air-side friction velocity and CD is a drag coefficient at the reference height of 10 m, derived elsewhere [60]. The energy content in turbulence is clearly an order of magnitude larger than that of the decadal to 0.5 h timescale. However, the effects of these energetic eddies produce water-side eddies (or waves) that are captured by f( U ). Therefore, it is safe to state from this analysis that extrapolating the meso-scale variance to a sub-daily time scale introduces a negligible correction to k (Figure 9) provided the appropriate f( U ) is used.

4. Conclusions

A new bias-correction model was developed based on a Taylor series expansion around the monthly mean wind state that accommodates any non-linearity in f( U ) and explicitly considers wind variability effects on gas transfer velocity. This new model, expressed as a function of the squared coefficient of variation, I u 2 , can be used to adjust k and is derived relative to the 6-hourly time-averaged wind products. For the 6-hourly high resolution CCMP wind data, a simplified bias correction kb as a function of I u 2 can be directly used to correct monthly k for any slightly soluble gas. With increasing wind variability over the last few decades associated with enhanced synoptic-scale high wind events [61,62], it is becoming increasingly necessary to quantify how the trend in wind variability biases k and influences global air–sea fluxes of CO2 and other climate-relevant gases (e.g., N2O, CH4) in the next generation of climate models. In the absence of other wind statistics, a plausible approximation to correct monthly k is to set I u 2 = 0.15 . This correction can be readily accommodated in current climate models. More broadly, the moment expansion approach presented here can be adapted to correct biases associated with averaging non-linear functions so as to accommodate measurements at different temporal or spatial resolutions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs13071328/s1, including Text S1, Tables S1 and S2, Figures S1–S3.

Author Contributions

Conceptualization, Y.G., G.G.K. and N.C.; Data curation, Y.G.; Formal analysis, Y.G., G.G.K. and N.C.; Supervision, G.G.K. and N.C.; Writing—original draft, Y.G.; Writing—review and editing, G.G.K. and N.C. All authors have read and agreed to the published version of the manuscript.

Funding

Y.G. is supported by a scholarship from the China Scholarship Council (CSC) under the grant CSC 201906710071. G.G.K. acknowledges support from the U.S. National Science Foundation (NSF-AGS-1644382, NSF-AGS-2028633, and NSF-IOS-1754893). N.C. is supported by the “Laboratoire d’Excellence” LabexMER (ANR-10-LABX-19) and co-funded by a grant from the French government under the program “Investissements d’Avenir”.

Acknowledgments

The cross-calibrated multiplatform (CCMP-V2) 6-hourly wind speeds are publicly available at www.remss.com/measurements/ccmp, accessed on 7 February 2021. The sea surface temperature data were obtained from https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview, accessed on 7 February 2021.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Katul, G.; Mammarella, I.; Grönholm, T.; Vesala, T. A structure function model recovers the many formulations for air-water gas transfer velocity. Water Resour. Res. 2018, 54, 5905–5920. [Google Scholar] [CrossRef]
  2. Esters, L.; Landwehr, S.; Sutherland, G.; Bell, T.G.; Christensen, K.H.; Saltzman, E.S.; Miller, S.D.; Ward, B. Parameterizing air-sea gas transfer velocity with dissipation. J. Geophys. Res. Oceans 2017, 122, 3041–3056. [Google Scholar] [CrossRef] [Green Version]
  3. Zappa, C.J.; McGillis, W.R.; Raymond, P.A.; Edson, J.B.; Hintsa, E.J.; Zemmelink, H.J.; Dacey, J.W.H.; Ho, D.T. Environmental turbulent mixing controls on air-water gas exchange in marine and aquatic systems. Geophys. Res. Lett. 2007, 34. [Google Scholar] [CrossRef] [Green Version]
  4. Soloviev, A.V. Coupled renewal model of ocean viscous sublayer, thermal skin effect and interfacial gas transfer velocity. J. Mar. Syst. 2007, 66, 19–27. [Google Scholar] [CrossRef]
  5. Vieira, V.M.N.C.; Mateus, M.; Canelas, R.; Leitão, F. The FuGas 2.5 Updated for the Effects of Surface Turbulence on the Transfer Velocity of Gases at the Atmosphere–Ocean Interface. J. Mar. Sci. Eng. 2020, 8, 435. [Google Scholar] [CrossRef]
  6. Liang, J.; Deutsch, C.; McWilliams, J.C.; Baschek, B.; Sullivan, P.P.; Chiba, D. Parameterizing bubble-mediated air-sea gas exchange and its effect on ocean ventilation. Glob. Biogeochem. Cycles 2013, 27, 894–905. [Google Scholar] [CrossRef] [Green Version]
  7. Woolf, D.K.; Leifer, I.S.; Nightingale, P.D.; Rhee, T.S.; Bowyer, P.; Caulliez, G.; de Leeuw, G.; Larsen, S.E.; Liddicoat, M.; Baker, J.; et al. Modelling of bubble-mediated gas transfer: Fundamental principles and a laboratory test. J. Mar. Syst. 2007, 66, 71–91. [Google Scholar] [CrossRef]
  8. Deike, L.; Melville, W.K. Gas transfer by breaking waves. Geophys. Res. Lett. 2018, 45, 410–482. [Google Scholar] [CrossRef]
  9. Zavarsky, A.; Goddijn Murphy, L.; Steinhoff, T.; Marandino, C.A. Bubble-Mediated Gas Transfer and Gas Transfer Suppression of DMS and CO2. J. Geophys. Res. Atmos. 2018, 123, 6624–6647. [Google Scholar] [CrossRef]
  10. Goddijn Murphy, L.; Woolf, D.K.; Callaghan, A.H.; Nightingale, P.D.; Shutler, J.D. A reconciliation of empirical and mechanistic models of the air-sea gas transfer velocity. J. Geophys. Res. Oceans 2016, 121, 818–835. [Google Scholar] [CrossRef] [Green Version]
  11. Van der Hoven, I. Power spectrum of horizontal wind speed in the frequency range from 0.0007 to 900 cycles per hour. J. Meteorol. 1957, 14, 160–164. [Google Scholar] [CrossRef] [Green Version]
  12. Stull, R.B. An Introduction to Boundary Layer Meteorology; Springer: Dordrecht, The Netherlands, 1988. [Google Scholar]
  13. Takahashi, T.; Sutherland, S.C.; Wanninkhof, R.; Sweeney, C.; Feely, R.A.; Chipman, D.W.; Hales, B.; Friederich, G.; Chavez, F.; Sabine, C.; et al. Climatological mean and decadal change in surface ocean pCO2, and net sea–air CO2 flux over the global oceans. Deep Sea Res. Part II Top. Stud. Oceanogr. 2009, 56, 554–577. [Google Scholar] [CrossRef]
  14. Gruber, N.; Gloor, M.; Mikaloff Fletcher, S.E.; Doney, S.C.; Dutkiewicz, S.; Follows, M.J.; Gerber, M.; Jacobson, A.R.; Joos, F.; Lindsay, K.; et al. Oceanic sources, sinks, and transport of atmospheric CO2. Glob. Biogeochem. Cycles 2009, 23, B1005. [Google Scholar] [CrossRef] [Green Version]
  15. Wanninkhof, R.; Park, G.H.; Takahashi, T.; Sweeney, C.; Feely, R.; Nojiri, Y.; Gruber, N.; Doney, S.C.; McKinley, G.A.; Lenton, A.; et al. Global ocean carbon uptake: Magnitude, variability and trends. Biogeosciences 2013, 10, 1983–2000. [Google Scholar] [CrossRef]
  16. Landschützer, P.; Gruber, N.; Bakker, D.C.E.; Schuster, U. Recent variability of the global ocean carbon sink. Glob. Biogeochem. Cycles 2014, 28, 927–949. [Google Scholar] [CrossRef] [Green Version]
  17. Rödenbeck, C.; Bakker, D.C.E.; Gruber, N.; Iida, Y.; Jacobson, A.R.; Jones, S.; Landschützer, P.; Metzl, N.; Nakaoka, S.; Olsen, A.; et al. Data-based estimates of the ocean carbon sink variability—First results of the Surface Ocean pCO2 Mapping intercomparison (SOCOM). Biogeosciences 2015, 12, 7251–7278. [Google Scholar] [CrossRef] [Green Version]
  18. Friedlingstein, P.; Jones, M.W.; O’sullivan, M.; Andrew, R.M.; Hauck, J.; Peters, G.P.; Peters, W.; Pongratz, J.; Sitch, S.; Quéré, C.L.; et al. Global Carbon Budget 2019. Earth Syst. Sci. Data 2019, 11, 1783–1838. [Google Scholar] [CrossRef] [Green Version]
  19. Ward, B.; Wanninkhof, R.; McGillis, W.R.; Jessup, A.T.; DeGrandpre, M.D.; Hare, J.E.; Edson, J.B. Biases in the air-sea flux of CO2 resulting from ocean surface temperature gradients. J. Geophys. Res. Oceans 2004, 109, C8S. [Google Scholar] [CrossRef] [Green Version]
  20. Takahashi, T.; Sutherland, S.C.; Chipman, D.W.; Goddard, J.G.; Ho, C.; Newberger, T.; Sweeney, C.; Munro, D.R. Climatological distributions of pH, pCO2, total CO2, alkalinity, and CaCO3 saturation in the global surface ocean, and temporal changes at selected locations. Mar. Chem. 2014, 164, 95–125. [Google Scholar] [CrossRef] [Green Version]
  21. Bakker, D.C.E.; Pfeil, B.; Smith, K.; Hankin, S.; Olsen, A.; Alin, S.R.; Cosca, C.; Harasawa, S.; Kozyr, A.; Nojiri, Y.; et al. An update to the Surface Ocean CO2; Atlas (SOCAT version 2). Earth Syst. Sci. Data 2014, 6, 69–90. [Google Scholar] [CrossRef] [Green Version]
  22. Woolf, D.K.; Land, P.E.; Shutler, J.D.; Goddijn-Murphy, L.M.; Donlon, C.J. On the calculation of air-sea fluxes of CO2 in the presence of temperature and salinity gradients. J. Geophys. Res. Oceans 2016, 121, 1229–1248. [Google Scholar] [CrossRef] [Green Version]
  23. Roobaert, A.; Laruelle, G.G.; Landschützer, P.; Regnier, P. Uncertainty in the global oceanic CO2 uptake induced by wind forcing: Quantification and spatial analysis. Biogeosciences 2018, 15, 1701–1720. [Google Scholar] [CrossRef] [Green Version]
  24. Woolf, D.K.; Shutler, J.D.; Goddijn Murphy, L.; Watson, A.J.; Chapron, B.; Nightingale, P.D.; Donlon, C.J.; Piskozub, J.; Yelland, M.J.; Ashton, I.; et al. Key uncertainties in the recent air—Sea flux of CO2. Glob. Biogeochem. Cycles 2019, 33, 1548–1563. [Google Scholar] [CrossRef] [Green Version]
  25. Chiodi, A.M.; Dunne, J.P.; Harrison, D.E. Estimating Air-Sea Carbon Flux Uncertainty over the Tropical Pacific: Importance of Winds and Wind Analysis Uncertainty. Glob. Biogeochem. Cycles 2019, 33, 370–390. [Google Scholar] [CrossRef]
  26. Wanninkhof, R. The effect of using time-averaged winds on regional air-sea CO2 fluxes. Geophys. Monogr. Am. Geophys. Union 2002, 127, 351–356. [Google Scholar]
  27. Fangohr, S.; Woolf, D.K.; Jeffery, C.D.; Robinson, I.S. Calculating long-term global air-sea flux of carbon dioxide using scatterometer, passive microwave, and model reanalysis wind data. J. Geophys. Res. 2008, 113, C9032. [Google Scholar] [CrossRef] [Green Version]
  28. Jiang, L.; Cai, W.; Wanninkhof, R.; Wang, Y.; Lüger, H. Air-sea CO2 fluxes on the U.S. South Atlantic Bight: Spatial and seasonal variability. J. Geophys. Res. 2008, 113, C7019. [Google Scholar] [CrossRef]
  29. Wanninkhof, R. Relationship between wind speed and gas exchange. J. Geophys. Res. Atmos. 1992, 97, 7373–7382. [Google Scholar] [CrossRef]
  30. 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]
  31. Feely, R.A.; Wanninkhof, R.; McGillis, W.; Carr, M.E.; Cosca, C.E. Effects of wind speed and gas exchange parameterizations on the air-sea CO2 fluxes in the equatorial Pacific Ocean. J. Geophys. Res. Oceans 2004, 109, C3S–C8S. [Google Scholar] [CrossRef] [Green Version]
  32. Sarmiento, J.L.G.N. Ocean Biogeochemical Dynamics; Princeton University Press: Princeton, NJ, USA, 2006. [Google Scholar]
  33. Wanninkhof, R. Relationship between wind speed and gas exchange over the ocean revisited. Limnol. Oceanogr. Methods 2014, 12, 351–362. [Google Scholar] [CrossRef]
  34. Nightingale, P.M.G.L. 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]
  35. 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]
  36. McGillis, W.R.; Edson, J.B.; Zappa, C.J.; Ware, J.D.; McKenna, S.P.; Terray, E.A.; Hare, J.E.; Fairall, C.W.; Drennan, W.; Donelan, M.; et al. Air-sea CO2 exchange in the equatorial Pacific. J. Geophys. Res. Oceans 2004, 109, C2S–C8S. [Google Scholar] [CrossRef] [Green Version]
  37. 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]
  38. Wanninkhof, R.; Asher, W.E.; Ho, D.T.; Sweeney, C.; 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] [Green Version]
  39. Prytherch, J.; Yelland, M.J.; Pascal, R.W.; Moat, B.I.; Skjelvan, I.; Srokosz, M.A. Open ocean gas transfer velocity derived from long-term direct measurements of the CO2 flux. Geophys. Res. Lett. 2010, 37, L23607. [Google Scholar] [CrossRef]
  40. 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] [Green Version]
  41. Sweeney, C.; Gloor, E.; Jacobson, A.R.; Key, R.M.; McKinley, G.; Sarmiento, J.L.; Wanninkhof, R. Constraining global air-sea gas exchange for CO2 with recent bomb14C measurements. Glob. Biogeochem. Cycles 2007, 21, B2015. [Google Scholar] [CrossRef]
  42. Atlas, R.; Hoffman, R.N.; Ardizzone, J.; Leidner, S.M.; Jusem, J.C.; Smith, D.K.; Gombos, D. A Cross-calibrated, Multiplatform Ocean Surface Wind Velocity Product for Meteorological and Oceanographic Applications. Bull. Am. Meteorol. Soc. 2011, 92, 157–174. [Google Scholar] [CrossRef]
  43. Wanninkhof, R.; Triñanes, J. The impact of changing wind speeds on gas transfer and its effect on global air-sea CO2 fluxes. Glob. Biogeochem. Cycles 2017, 31, 961–974. [Google Scholar] [CrossRef] [Green Version]
  44. Kent, E.C.; Fangohr, S.; Berry, D.I. A comparative assessment of monthly mean wind speed products over the global ocean. Int. J. Climatol. 2013, 33, 2520–2541. [Google Scholar] [CrossRef]
  45. Li, Q.; Guo, X.; Zhai, W.; Xu, Y.; Dai, M. Partial pressure of CO2 and air-sea CO2 fluxes in the South China Sea: Synthesis of an 18-year dataset. Prog. Oceanogr. 2020, 182, 102272. [Google Scholar] [CrossRef]
  46. Service, C.C.C. ERA5: Fifth Generation of ECMWF Atmospheric Reanalyses of the Global Climate. Available online: https://confluence.ecmwf.int/display/CKB/ERA5%3A+data+documentation (accessed on 7 February 2021).
  47. Etcheto, J.B.A.J. Estimating the chemical enhancement effect on the air-sea CO2 exchange using the ERS1 scatterometer wind speeds. In Air-Water Gas Transfer; AEON Verlag & Studio: Hanau, Germany, 1995; pp. 827–841. [Google Scholar]
  48. Wanninkhof, R.; Sullivan, K.F.; Top, Z. Air-sea gas transfer in the Southern Ocean. J. Geophys. Res. Oceans 2004, 109, C8S–C19S. [Google Scholar] [CrossRef]
  49. Lüger, H.; Wanninkhof, R.; Wallace, D.W.R.; Körtzinger, A. CO2 fluxes in the subtropical and subarctic North Atlantic based on measurements from a volunteer observing ship. J. Geophys. Res. 2006, 111, L6024. [Google Scholar] [CrossRef] [Green Version]
  50. Zhang, X.; Cai, W. On some biases of estimating the global distribution of air-sea CO2 flux by bulk parameterizations. Geophys. Res. Lett. 2007, 34, L1608. [Google Scholar] [CrossRef]
  51. Xue, L.; Xue, M.; Zhang, L.; Sun, T.; Guo, Z.; Wang, J. Surface partial pressure of CO2 and air-sea exchange in the northern Yellow Sea. J. Mar. Syst. 2012, 105–108, 194–206. [Google Scholar] [CrossRef]
  52. Astor, Y.M.; Lorenzoni, L.; Guzman, L.; Fuentes, G.; Muller-Karger, F.; Varela, R.; Scranton, M.; Taylor, G.T.; Thunell, R. Distribution and variability of the dissolved inorganic carbon system in the Cariaco Basin, Venezuela. Mar. Chem. 2017, 195, 15–26. [Google Scholar] [CrossRef]
  53. Olsen, A.W.R.T. The effect of wind speed products and wind speed-gas exchange relationships on interannual variability of the air-sea CO2 gas transfer velocity. Tellus B 2005, 57, 95–106. [Google Scholar] [CrossRef] [Green Version]
  54. Sampe, T.A.X.S. Mapping high sea winds from space: A global climatology. Bull. Am. Meteorol. Soc. 2007, 12, 1965–1978. [Google Scholar] [CrossRef] [Green Version]
  55. Liang, J.H.; D’Asaro, E.A.; McNeil, C.L.; Fan, Y.; Harcourt, R.R.; Emerson, S.R.; Yang, B.; Sullivan, P.P. Suppression of CO2 Outgassing by Gas Bubbles Under a Hurricane. Geophys. Res. Lett. 2020, 47, e2020G–e90249G. [Google Scholar] [CrossRef]
  56. Atamanchuk, D.; Koelling, J.; Send, U.; Wallace, D.W.R. Rapid transfer of oxygen to the deep ocean mediated by bubbles. Nat. Geosci. 2020, 13, 232–237. [Google Scholar] [CrossRef]
  57. Leighton, T.G.; Coles, D.G.H.; Srokosz, M.; White, P.R.; Woolf, D.K. Asymmetric transfer of CO2 across a broken sea surface. Sci. Rep. 2018, 8, 1–9. [Google Scholar] [CrossRef]
  58. Charney, J.G. Geostrophic turbulence. J. Atmos. Sci. 1971, 28, 1087–1095. [Google Scholar] [CrossRef] [Green Version]
  59. Nastrom, G.D. A climatology of atmospheric wavenumber spectra of wind and temperature observed by commercial aircraft. J. Atmos. Sci. 1985, 42, 950–960. [Google Scholar] [CrossRef] [Green Version]
  60. Large, W.B. Surface fluxes for practitioners of global ocean data assimilation. In Ocean Weather Forecasting; Springer: Dordrecht, The Netherlands, 2006; pp. 227–270. [Google Scholar]
  61. Young, I.R.; Zieger, S.; Babanin, A.V. Global trends in wind speed and wave height. Science 2011, 332, 451–455. [Google Scholar] [CrossRef]
  62. Young, I.R.; Ribal, A. Multiplatform evaluation of global trends in wind speed and wave height. Science 2019, 364, 548–552. [Google Scholar] [CrossRef]
Figure 1. Conceptual diagram representing the bias in gas transfer velocity (k) estimates associated with averaging wind speed variability (adapted from [32]). The quadratic and cubic relations are in blue and orange, respectively.
Figure 1. Conceptual diagram representing the bias in gas transfer velocity (k) estimates associated with averaging wind speed variability (adapted from [32]). The quadratic and cubic relations are in blue and orange, respectively.
Remotesensing 13 01328 g001
Figure 2. Bias in k of CO2 due to wind speeds at varying spatial resolutions (0.5°   ×   0.5° and 5°   ×   5°) for 6-hourly and monthly gas transfer velocity (k), and temporal bias in k (6 hourly and monthly) at the spatial resolution of 0.5°   ×   0.5° and 5°   ×   5°. The kmon and k6h are gas transfer velocities averaged over all k values estimated from monthly and 6-hourly wind speed records, respectively. k5° and k0.5° are gas transfer velocities averaged over all k values estimated from 5° and 0.5° wind speed, respectively. The bias is estimated as △k*100/k6h,0.5° (k6h,0.5° is k at the resolution of 6-hourly and 0.5°   ×   0.5°).
Figure 2. Bias in k of CO2 due to wind speeds at varying spatial resolutions (0.5°   ×   0.5° and 5°   ×   5°) for 6-hourly and monthly gas transfer velocity (k), and temporal bias in k (6 hourly and monthly) at the spatial resolution of 0.5°   ×   0.5° and 5°   ×   5°. The kmon and k6h are gas transfer velocities averaged over all k values estimated from monthly and 6-hourly wind speed records, respectively. k5° and k0.5° are gas transfer velocities averaged over all k values estimated from 5° and 0.5° wind speed, respectively. The bias is estimated as △k*100/k6h,0.5° (k6h,0.5° is k at the resolution of 6-hourly and 0.5°   ×   0.5°).
Remotesensing 13 01328 g002
Figure 3. Mean bias in gas transfer velocity (k) for CO2 estimated from term 1 (measured bias in f(U)) and term 2 (bias correction kb from new model) of Equation (9) over the period spanning 1990 to 2018 for the parameterizations presented in Table 1.
Figure 3. Mean bias in gas transfer velocity (k) for CO2 estimated from term 1 (measured bias in f(U)) and term 2 (bias correction kb from new model) of Equation (9) over the period spanning 1990 to 2018 for the parameterizations presented in Table 1.
Remotesensing 13 01328 g003
Figure 4. Left panel: time series of global (a) monthly averaged wind speed U (in black) and standard deviation ( σ u ,   in   grey ) around U , (b) monthly squared coefficient of variation I u 2 = ( σ u / U ) 2 from 1990 to 2018 (note the small variations along the ordinate axis). The black and the grey dashed lines in (b) indicate the long-term trend (0.002 dec−1) and average ( I u 2 = 0.15 ), respectively. Right panel: spatial distribution of (c) trends in the wind speed standard deviation ( σ u ) around U , (d) monthly averaged wind speed U , and (e) monthly squared coefficient of variation I u 2 = ( σ u / U ) 2 from 1990 to 2018.
Figure 4. Left panel: time series of global (a) monthly averaged wind speed U (in black) and standard deviation ( σ u ,   in   grey ) around U , (b) monthly squared coefficient of variation I u 2 = ( σ u / U ) 2 from 1990 to 2018 (note the small variations along the ordinate axis). The black and the grey dashed lines in (b) indicate the long-term trend (0.002 dec−1) and average ( I u 2 = 0.15 ), respectively. Right panel: spatial distribution of (c) trends in the wind speed standard deviation ( σ u ) around U , (d) monthly averaged wind speed U , and (e) monthly squared coefficient of variation I u 2 = ( σ u / U ) 2 from 1990 to 2018.
Remotesensing 13 01328 g004
Figure 5. Difference in 6-hourly k and corrected k for CO2 applying five correction methodologies in reference to the 6-hourly k (in %) for k parameterizations listed in Table 1. The bias is estimated as △k*100/k6h,0.5°.
Figure 5. Difference in 6-hourly k and corrected k for CO2 applying five correction methodologies in reference to the 6-hourly k (in %) for k parameterizations listed in Table 1. The bias is estimated as △k*100/k6h,0.5°.
Remotesensing 13 01328 g005
Figure 6. (a) Spatial distribution of averaged variance of sea surface temperature (SST) ( σ SST 2 ) around monthly averaged SST ; (b) Time series of annual averaged variance of SST ( σ SST 2 ) around monthly averaged SST , the dashed line indicates the long-term trend; (c) Spatial pattern of trend in averaged variance of SST ( σ SST 2 ) around monthly averaged SST from 1990 to 2018.
Figure 6. (a) Spatial distribution of averaged variance of sea surface temperature (SST) ( σ SST 2 ) around monthly averaged SST ; (b) Time series of annual averaged variance of SST ( σ SST 2 ) around monthly averaged SST , the dashed line indicates the long-term trend; (c) Spatial pattern of trend in averaged variance of SST ( σ SST 2 ) around monthly averaged SST from 1990 to 2018.
Remotesensing 13 01328 g006
Figure 7. (a) Zonal profiles of corrected k for CO2 using the five correction methodologies in comparison to annual k derived from 6-hourly (red solid curve) and monthly (red dashed curve) wind speed. Zonal variation in k estimated using method 1 (in black) is not visible because it overlaps with the 6-hourly k. Panels (a1a9) show the latitudinal variations in nine k parameterizations listed in Table 1. (b) The RMSE of each method in corrected k from 6-hourly k.
Figure 7. (a) Zonal profiles of corrected k for CO2 using the five correction methodologies in comparison to annual k derived from 6-hourly (red solid curve) and monthly (red dashed curve) wind speed. Zonal variation in k estimated using method 1 (in black) is not visible because it overlaps with the 6-hourly k. Panels (a1a9) show the latitudinal variations in nine k parameterizations listed in Table 1. (b) The RMSE of each method in corrected k from 6-hourly k.
Remotesensing 13 01328 g007
Figure 8. Coefficient of variation I u 2 as a function of the averaging period Δ t   (from 6-hourly to monthly). Circles indicate the results from measurements, and the solid line represents a modelled fit through the measurements. For ∆t > 18 days, Iu2 becomes independent of ∆t. Global climate models operate on a ∆t = 30 days.
Figure 8. Coefficient of variation I u 2 as a function of the averaging period Δ t   (from 6-hourly to monthly). Circles indicate the results from measurements, and the solid line represents a modelled fit through the measurements. For ∆t > 18 days, Iu2 becomes independent of ∆t. Global climate models operate on a ∆t = 30 days.
Remotesensing 13 01328 g008
Figure 9. Energy spectrum of global average 6-hourly wind speed. The spectrum is extrapolated from 12 h to a turbulence scale (seconds) via Kolmogorov’s –5/3 power law (f–5/3, blue dashed line). The resolved spectrum has an exponent of –3 from multi-day to 12 h (f–3, blue solid line) consistent with an enstrophy cascade in a quasi-geostrophic flow. The dashed vertical lines (right to left) indicate frequencies corresponding to the following timescales: sub-hour (=0.5 h), diurnal (=12 h), daily (=24 h), and annual (=8760 h), respectively. The σ d 2 , σ m 2 , and σ t 2   refer to the variance at large (mesoscale to decadal), intermediate (12 h to turbulence), and small (turbulence) scales, respectively.
Figure 9. Energy spectrum of global average 6-hourly wind speed. The spectrum is extrapolated from 12 h to a turbulence scale (seconds) via Kolmogorov’s –5/3 power law (f–5/3, blue dashed line). The resolved spectrum has an exponent of –3 from multi-day to 12 h (f–3, blue solid line) consistent with an enstrophy cascade in a quasi-geostrophic flow. The dashed vertical lines (right to left) indicate frequencies corresponding to the following timescales: sub-hour (=0.5 h), diurnal (=12 h), daily (=24 h), and annual (=8760 h), respectively. The σ d 2 , σ m 2 , and σ t 2   refer to the variance at large (mesoscale to decadal), intermediate (12 h to turbulence), and small (turbulence) scales, respectively.
Remotesensing 13 01328 g009
Table 1. The f( U ) parameterizations used in estimating gas transfer velocity for CO2 (same expressions can be used for other slightly soluble gases [33]). The f( U ) formulations developed from long-term wind speeds (i.e., monthly) were not considered here. In Serial No.9, because the three equations are identical in form with a small difference in their α coefficients, we use the expression f( U ) = 0.251   U 2 [25] in the following analysis as a representative equation of all three models.
Table 1. The f( U ) parameterizations used in estimating gas transfer velocity for CO2 (same expressions can be used for other slightly soluble gases [33]). The f( U ) formulations developed from long-term wind speeds (i.e., monthly) were not considered here. In Serial No.9, because the three equations are identical in form with a small difference in their α coefficients, we use the expression f( U ) = 0.251   U 2 [25] in the following analysis as a representative equation of all three models.
Serial No.Reference f ( U ) Parameterization for CO2
f ( U )   =   k   ( S c / 660 ) 1 / 2
1Wanninkhof (1992) [29] 0.31 U 2
2Wanninkhof and McGillis (1999) [30] 0.0283 U 3
3Nightingale et al. (2000) [34] 0.222 U 2 + 0.333 U
4McGillis et al. (2001) [35]0.026 U 3 + 3.3
5McGillis et al. (2004) [36]0.014 U 3 + 8.2
6Weiss et al. (2007) [37]0.365 U 2 + 0.46 U
7Wanninkhof et al. (2009) [38]0.011 U 3 + 0.064 U 2 + 0.1 U + 3
8Prytherch et al. (2010) [39]0.034 U 3 + 5.3
9Ho et al (2006) [40], Sweeney et al. (2007) [41], Wanninkhof (2014) [33] α U 2
(where α   = 0.266/0.27/0.251)
Table 2. Summary of all the correction methodologies for CO2 and other gases.
Table 2. Summary of all the correction methodologies for CO2 and other gases.
MethodReferenceCorrectionCorrection Details
1This studykb from Equation (11) (for quadratic relations) and Equation (13) (for cubic relations) are added to f( U ) to estimate the corrected k.Grid-by-grid spatially multi-year mean kb
2This studyA simplified method using overall averaged value of kb to fix the bias.
3Wanninkhof (2002) [26](1) The corrected k with multiplier correction R2 (Equation (5)) for the quadratic parameterization is in the form of f ( U ) =   a   R 2 U 2 ,
(2) For the cubic relation with multiplier correction R3 (Equation (6)), the corrected f( U ) is expressed as
f ( U ) =   a   R 3   U 3 + b   U 2 + d   U + e
Assuming a Rayleigh distribution of the 6-hourly wind speeds, R 2 =   Γ ( 2 ) [ Γ ( 3 / 2 ) ] 2 = 1.27 and R 3 =   Γ ( 5 / 2 ) [ Γ ( 3 / 2 ) ] 3 = 1.91 (See Text S1 in Supplementary for details).
4Jiang et al. (2008) [28]Global averaged multiplier correction factors R2 and R3 are estimated using the measured 6-hourly wind speed with R2 = 1.23 and R3 = 1.78.
5Fangohr et al. (2008) [27]Zonal averaged R2 and R3 are used. Large gradients in zonal R2 and R3 are because of the large zonal gradients in wind variance (Figure S1).
Table 3. Parameters used to run scenarios of imposed changes in wind speed ( U ) and sea surface temperature (SST) and their effects on the gas transfer velocity (k) for the nine k parameterizations featured in Table 1. The starting values of U and SST were set to 6.48 m s−1 and 13.73 °C, respectively, according to their climatological global mean. The sensitivities of k to U and SST were assessed from the ratio of the percentage change in k (Y) to percentage change in each factor (X) using the equation: sensitivity = (∆Y/Y)/(∆X/X).
Table 3. Parameters used to run scenarios of imposed changes in wind speed ( U ) and sea surface temperature (SST) and their effects on the gas transfer velocity (k) for the nine k parameterizations featured in Table 1. The starting values of U and SST were set to 6.48 m s−1 and 13.73 °C, respectively, according to their climatological global mean. The sensitivities of k to U and SST were assessed from the ratio of the percentage change in k (Y) to percentage change in each factor (X) using the equation: sensitivity = (∆Y/Y)/(∆X/X).
Serial NOStarting ValueImposed ChangeImposed Change
U SST U SST
U   m · s 1 SST (°C)2%4%8%2%3%4%2%4%8%2%3%4%
Δkk Sensitivity
16.8413.730.491.002.040.090.140.182.022.042.080.380.380.38
20.470.951.980.060.090.113.063.123.250.380.380.38
30.390.791.610.080.120.161.841.851.890.380.380.38
40.430.881.820.070.110.152.192.242.320.380.380.38
50.230.470.980.080.120.161.081.101.150.380.380.38
60.420.861.750.080.120.162.022.042.080.380.380.38
70.430.871.770.080.120.162.022.042.080.380.380.38
80.300.601.240.060.100.131.721.741.800.380.380.38
90.400.811.650.070.110.152.022.042.080.380.380.38
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gu, Y.; Katul, G.G.; Cassar, N. Mesoscale Temporal Wind Variability Biases Global Air–Sea Gas Transfer Velocity of CO2 and Other Slightly Soluble Gases. Remote Sens. 2021, 13, 1328. https://doi.org/10.3390/rs13071328

AMA Style

Gu Y, Katul GG, Cassar N. Mesoscale Temporal Wind Variability Biases Global Air–Sea Gas Transfer Velocity of CO2 and Other Slightly Soluble Gases. Remote Sensing. 2021; 13(7):1328. https://doi.org/10.3390/rs13071328

Chicago/Turabian Style

Gu, Yuanyuan, Gabriel G. Katul, and Nicolas Cassar. 2021. "Mesoscale Temporal Wind Variability Biases Global Air–Sea Gas Transfer Velocity of CO2 and Other Slightly Soluble Gases" Remote Sensing 13, no. 7: 1328. https://doi.org/10.3390/rs13071328

APA Style

Gu, Y., Katul, G. G., & Cassar, N. (2021). Mesoscale Temporal Wind Variability Biases Global Air–Sea Gas Transfer Velocity of CO2 and Other Slightly Soluble Gases. Remote Sensing, 13(7), 1328. https://doi.org/10.3390/rs13071328

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