Articles | Volume 21, issue 10
Research article
26 May 2021
Research article |  | 26 May 2021

Uncertainties in eddy covariance air–sea CO2 flux measurements and implications for gas transfer velocity parameterisations

Yuanxu Dong, Mingxi Yang, Dorothee C. E. Bakker, Vassilis Kitidis, and Thomas G. Bell

Air–sea carbon dioxide (CO2) flux is often indirectly estimated by the bulk method using the air–sea difference in CO2 fugacity (ΔfCO2) and a parameterisation of the gas transfer velocity (K). Direct flux measurements by eddy covariance (EC) provide an independent reference for bulk flux estimates and are often used to study processes that drive K. However, inherent uncertainties in EC air–sea CO2 flux measurements from ships have not been well quantified and may confound analyses of K. This paper evaluates the uncertainties in EC CO2 fluxes from four cruises. Fluxes were measured with two state-of-the-art closed-path CO2 analysers on two ships. The mean bias in the EC CO2 flux is low, but the random error is relatively large over short timescales. The uncertainty (1 standard deviation) in hourly averaged EC air–sea CO2 fluxes (cruise mean) ranges from 1.4 to 3.2 mmolm-2d-1. This corresponds to a relative uncertainty of  20 % during two Arctic cruises that observed large CO2 flux magnitude. The relative uncertainty was greater ( 50 %) when the CO2 flux magnitude was small during two Atlantic cruises. Random uncertainty in the EC CO2 flux is mostly caused by sampling error. Instrument noise is relatively unimportant. Random uncertainty in EC CO2 fluxes can be reduced by averaging for longer. However, averaging for too long will result in the inclusion of more natural variability. Auto-covariance analysis of CO2 fluxes suggests that the optimal timescale for averaging EC CO2 flux measurements ranges from 1 to 3 h, which increases the mean signal-to-noise ratio of the four cruises to higher than 3. Applying an appropriate averaging timescale and suitable ΔfCO2 threshold (20 µatm) to EC flux data enables an optimal analysis of K.

1 Introduction

Since the Industrial Revolution, atmospheric CO2 levels have risen steeply due to human activities (Broecker and Peng, 1998). The ocean plays a key role in the global carbon cycle, having taken up roughly one-quarter of anthropogenic CO2 emissions over the last decade (Friedlingstein et al., 2020). Accurate estimates of air–sea CO2 flux are vital to forecast climate change and to quantify the effects of ocean CO2 uptake on the marine biosphere.

Air–sea CO2 flux (F, e.g. in mmolm-2d-1) is typically estimated indirectly by the bulk equation:

(1) F = K 660 ( Sc / 660 ) - 0.5 α ( f CO 2 w - f CO 2 a ) ,

where K660 (in cm h−1) is the gas transfer velocity, usually parameterised as a function of wind speed (e.g. Nightingale et al., 2000); Sc (dimensionless) is the Schmidt number (Wanninkhof, 2014); and α (molL-1atm-1) is the solubility (Weiss, 1974). Sc is equal to 660 for CO2 at 20 C and 35 ‰ saltwater (Wanninkhof et al., 2009). fCO2w and fCO2a are the CO2 fugacity (in µatm) at the sea surface and in the overlying atmosphere, respectively, with fCO2wfCO2a the air–sea CO2 fugacity difference (ΔfCO2). Uncertainties in the K660 parameterisation and limited coverage of fCO2w measurements result in considerable uncertainties in global bulk flux estimates (Takahashi et al., 2009; Woolf et al., 2019).

Eddy covariance (EC) is the most direct method for measuring the air–sea CO2 flux F:

(2) F = ρ w c ,

where ρ is the mean mole density of dry air (e.g. in mol m−3). The dry CO2 mixing ratio c (in ppm or µmol mol−1) is measured by a fast-response gas analyser, and the vertical wind velocity w (in m s−1) is often measured by a sonic anemometer. The prime denotes the fluctuations from the mean, while the overbar indicates time average. Equation (2) does not rely on ΔfCO2 measurements, nor empirical parameters and assumptions of the gas properties (Wanninkhof, 2014). EC flux measurements can therefore be considered useful as an independent reference for bulk air–sea CO2 flux estimates. Furthermore, the typical temporal and spatial scales of EC flux measurements are ca. hourly and 1–10 km2. These scales are much smaller than the temporal and spatial scales of alternative techniques for measuring gas transfer, e.g. by dual tracer methods (daily and 1000 km2) (Nightingale et al., 2000; Ho et al., 2006). EC measurements are thus potentially better suited to capture variations in gas exchange due to small-scale processes at the air–sea interface (Garbe et al., 2014).

The EC CO2 flux method has developed and improved over time. Before 1990, EC was successfully used to measure air–sea momentum and heat fluxes. EC air–sea CO2 flux measurements made during those times were unreasonably high (Jones and Smith, 1977; Wesely et al., 1982; Smith and Jones, 1985; Broecker et al., 1986). After 1990, with the development of the infrared gas analyser, EC became routinely used for terrestrial carbon cycle research (Baldocchi et al., 2001). Development of the EC method was accompanied by improvements in the flux uncertainty analysis, which was generally based on momentum, heat and land–atmosphere gas flux measurements (Lenschow and Kristensen, 1985; Businger, 1986; Lenschow et al., 1994; Wienhold et al., 1995; Mahrt, 1998; Finkelstein and Sims, 2001; Loescher et al., 2006; Rannik et al., 2009, 2016; Billesbach, 2011; Mauder et al., 2013; Langford et al., 2015; Post et al., 2015).

In the late 1990s, the advancement in motion correction of wind measurements (Edson et al., 1998; Yelland et al., 1998) facilitated ship-based EC CO2 flux measurements from a moving platform (McGillis et al., 2001; 2004). After 2000, a commercial open-path infrared gas analyser LI-7500 (Li-COR Inc. USA) became widely used for air–sea CO2 flux measurements (Weiss et al., 2007; Kondo and Tsukamoto, 2007; Prytherch et al., 2010; Edson et al., 2011; Else et al., 2011; Lauvset et al., 2011). The LI-7500 generated extremely large and highly variable CO2 fluxes in comparison to expected fluxes (Kondo and Tsukamoto, 2007; Prytherch et al., 2010; Edson et al., 2011; Else et al., 2011; Lauvset et al., 2011). This problem is generally considered to be an artefact caused by water vapour cross-sensitivity (Kohsiek, 2000; Prytherch et al., 2010; Edson et al., 2011; Landwehr et al., 2014). Mathematical corrections proposed to address this artefact (Edson et al., 2011; Prytherch et al., 2010) were later shown to be unsatisfactory (Else et al., 2011; Ikawa et al., 2013; Blomquist et al., 2014; Tsukamoto et al., 2014) or incorrect (Landwehr et al., 2014).

The most reliable method for measuring EC air–sea CO2 fluxes involves the physical removal of water vapour fluctuations from the sampled air. The simplest approach is to combine a closed-path gas analyser with a physical dryer to eliminate most of the water vapour fluctuation (Miller et al., 2010; Blomquist et al., 2014; Landwehr et al., 2014; Yang et al., 2016; Nilsson et al., 2018). The tuneable-diode-laser-based cavity ring-down spectrometer (CRDS) made by Picarro Inc. (Santa Clara, California, USA) is the most precise closed-path analyser currently available (Blomquist et al., 2014). The closed-path infrared gas analyser LI-7200 (LI-COR Biosciences, Lincoln, Nebraska, USA) is another popular choice.

The advancements in instrumentation and in motion correction methods have significantly improved the quality of air–sea EC CO2 flux observations, but, despite these changes, the flux uncertainties have not been well quantified. The aims of this study are to (1) analyse uncertainties in EC air–sea CO2 flux measurements; (2) propose practical methods to reduce the systematic and random flux uncertainty; and (3) investigate how the EC flux uncertainty influences our ability to estimate and parameterise K660.

Table 1Basic information for all four cruises on the RRS James Clark Ross (JCR) and RRS Discovery that measured air–sea EC CO2 fluxes.

Download Print Version | Download XLSX

2 Experiment and methods

2.1 Instrumental setup

The basic information of four cruises is summarised in Table 1. Appendix A shows the four cruise tracks (Figs. A1 and A2). Data from the Atlantic cruises (AMT28 and AMT29) are limited to 3 N–20 S in order to focus specifically on the performance of two different gas analysers in the same region with low flux signal (tropical zone).

Figure 1EC system (upper panel) and a diagram of the system setup (bottom panel). EC instruments: (1) sonic anemometer, (2) motion sensor, (3) air sample inlet for gas analyser, (4) data logger/gas analyser. Arctic and Atlantic data from 2018 were collected on the RRS James Clark Ross (JCR; upper right) using a Picarro G2311-f, and Atlantic data from 2019 were collected using a LI-7200 on the RRS Discovery (upper left).


The CO2 flux and data logging systems installed on the JCR and Discovery were operated autonomously. The EC systems were approximately 20 ma.m.s.l. on both ships (at the top of the foremasts; Fig. 1) to minimise flow distortion and exposure to sea spray. Computational fluid dynamics (CFD) simulation indicates that the airflow distortion at the top of the JCR foremast is small ( 1 % of the free stream wind speed when the ship is head to wind; Moat and Yelland, 2015). The hull structure of RRS Discovery is nearly identical to that of RRS James Cook. CFD simulation of the James Cook indicates that the airflow at the top foremast is distorted by  2 % for bow-on flows (Moat et al., 2006). The deflection of the streamline from horizontal and effects on the vertical wind component is accounted for by the double rotation (motion correction processes; see Sect. 2.2) prior to the EC flux calculation for both ships.

The EC system on the JCR consists of a three-dimensional sonic anemometer (Metek Inc., Sonic-3 Scientific), a motion sensor (initially Systron Donner Motionpak II, which compared favourably with and was then replaced by a Life Performance-Research LPMS-RS232AL2 in April 2019), and a Picarro G2311-f gas analyser. All instruments sampled at a frequency of 10 Hz or greater, and the data were logged at 10 Hz with a data logger (CR6, Campbell Scientific, Inc.), similar to the setup by Butterworth and Miller (2016). Air is pulled through a long tube (30 m, inner diameter 0.95 cm, Reynolds number 5957) with a dry vane pump at a flow rate of  40 L min−1 (Gast 1023 series). The Picarro gas analyser subsamples from this tube through a particle filter (Swagelok 2 µm) and a dryer (Nafion PD-200T-24M) at a flow of  5 L min−1 (Fig. 1). The dryer is set up in the “re-flux” configuration and uses the lower pressure Picarro exhaust to dry the sample air. This method removes  80 % of the water vapour and essentially all of the humidity fluctuations (Yang et al., 2016). The Picarro internal calculation accounts for the detected residual water vapour and yields a dry CO2 mixing ratio that is used in the flux calculations. A valve controlled by the Picarro instrument injects a “puff” of nitrogen (N2) into the tip of the inlet tube for 30 s every 6 h. This enables estimates of the time delay and high-frequency signal attenuation (Sect. 2.2).

The EC system on RRS Discovery consists of a Gill R3-50 sonic anemometer, a LPMS motion sensor package, and a LI-7200 gas analyser. The LI-7200 gas analyser was mounted within the enclosed staircase, directly underneath the meteorological platform and close to the inlet (inlet length 7.5 m, inner diameter 0.95 cm, Reynolds number 1042). A single pump (Gast 1023) was sufficient to pull air through a particle filter (Swagelok 2 µm), a dryer (Nafion PD-200T-24M), and the LI-7200 at a flow of  7 L min−1. There was no N2 puff system setup on Discovery, but equivalent lab tests confirmed that the delay time was less than on the JCR because of the shorter inlet line. The dryer on the Discovery is set up in the same re-flux configuration as the JCR and uses the lower pressure at the LI-7200 exhaust (limited by an additional 0.08 cm diameter critical orifice) to dry the sample air. This setup removes  60 %–70 % of the water vapour and essentially all of the humidity fluctuations. The dry CO2 mixing ratio, computed by accounting for the LI-7200 temperature, pressure, and residual water vapour measurements, is used in the flux calculations.

Figure 2Flow chart of EC data processing. The raw high-frequency (10 Hz) wind and CO2 data were initially processed separately and then combined to calculate fluxes. CO2 fluxes were filtered by a series of data quality control criteria. The 20 min flux intervals were averaged to longer timescales (hourly or more). The data processing is detailed in the text.


2.2 Flux processing

The EC air–sea CO2 flux calculation steps using the raw data are outlined with a flow chart (Fig. 2) and detailed below. The raw high-frequency wind and CO2 data are processed first, yielding fluxes in 20 min averaging time interval and related statistics. These statistics are then used for quality control of the fluxes. Further averaging of the quality-controlled 20 min fluxes to hourly or longer timescales is then used to reduce random error (Sect. 4.1). Linear detrending was used to identify the turbulent fluctuations (i.e. w and c) throughout the analyses.

To correct the wind data for ship motion, we first generated hourly data files containing the measurements from the sonic anemometer (three-dimensional wind speed components: u, v, and w and sonic temperature Ts), motion sensor (three axis accelerations: accel_x, accel_y, accel_z; and rotation angles: rot_x, rot_y, rot_z), ship heading over ground (HDG; from the gyro compass), and ship speed over ground (SOG; from Global Position System). Spikes larger than 4 standard deviations (SDs) from the median were removed. Secondly, a complementary filtering method using Euler angles (see Edson et al., 1998) was applied to the hourly data files to remove apparent winds generated by the ship movements. The motion-corrected winds were further decorrelated against ship motion to remove any residual motion-sensitivity (Miller et al., 2010; Yang et al., 2013). The motion-corrected winds were double-rotated to account for the wind streamline over the ship, yielding the vertical wind velocity (w) required in Eq. (2). Inspection of frequency spectra showed that the spectral peak at the ship motion frequencies (approximately 0.1–0.3 Hz) had disappeared after the motion correction (Fig. S1 in the Supplement). This indicates that the majority of ship motion had been removed from the measured wind speed. The last step in the wind data processing was the calculation of 20 min average friction velocity, sensible heat flux, and other key variables used for data quality control (Table S1 in the Supplement).

The CO2 data were de-spiked (by removing values > 4 SDs from the median). The Picarro CO2 mixing ratio was further decorrelated against analyser cell pressure and temperature to remove CO2 variations due to the ship's motion. The LI-7200 CO2 mixing ratio was further decorrelated against the LI-7200 H2O mixing ratio and temperature to remove residual air density fluctuations, following Landwehr et al. (2018). CO2 data were also decorrelated against the ship's heave and accelerations because these can produce spurious CO2 variability (Miller et al., 2010; Blomquist et al., 2014).

A lag between CO2 data acquisition and the wind data is created because of the time taken for sample air to travel through the inlet tube. On the JCR, we use the puff system where the lag time is the time difference between the N2 puff start (when the on/off valve is switched) and the time when the diluted signal is sensed by the gas analyser. The lag time can also be estimated by the maximum covariance method, calculated by shifting the time base of the CO2 signal and finding the shift that achieves maximum covariance between the vertical wind velocity (w) signal and the shifted CO2 signal. The lag times estimated by the maximum covariance method agree well with the estimates of the puff procedure (Fig. S2 in the Supplement). These estimates indicate a lag time of 3.3–3.4 s for the Arctic cruises and 3.3 s for cruise AMT28 on the JCR. The lag time on Discovery (AMT29) estimated by the maximum covariance method was 2.6 s, consistent with laboratory test results prior to the cruise.

The inlet tube, particle filter, and dryer cause high-frequency CO2 flux signal attenuation. The N2 puff was also used to assess the response time by considering the e-folding time in the CO2 signal change (similar approaches have been used by Bariteau et al., 2010; Blomquist et al., 2014, Bell et al., 2015). The response time is 0.35 s for the EC system on JCR and 0.25 s for the EC system on Discovery (estimated in the laboratory prior to cruise). These response times were combined with the relative wind speed-dependent, theoretical shapes of the cospectra (Kaimal et al., 1972) to estimate the percentage flux loss due to the inlet attenuation (Yang et al., 2013). The mean attenuation percentage is less than 10 %, with a relative wind speed dependence (Fig. S3 in the Supplement). The attenuation percentage value was applied to the computed flux to compensate for the flux loss due to the high-frequency signal attenuation. Finally, horizontal CO2 fluxes and other statistics such as CO2 range and CO2 trend were computed for quality control purposes (Table S1, Supplement).

The computed 20 min fluxes were filtered for non-ideal ship manoeuvres or violations of the homogeneity/stationary requirement of EC (see Supplement for the quality control criteria).

2.3 Uncertainty analysis methods

2.3.1 Uncertainty components

Uncertainty contains two components: systematic error (δFS) and random error (δFR). According to propagation of uncertainty theory (JCGM, 2008), the total uncertainty in EC CO2 fluxes (from random and systematic errors) can be expressed as

(3) δ F = δ F R 2 + δ F S 2 .

Systematic errors (Sect. 2.3.2) will cause bias in the flux. They thus should be eliminated/minimised with the appropriate system setup and, if needed, effective numerical corrections. Random error results in imprecision (but not bias) and can be reduced by averaging repeated measurements (Sect. 2.3.3). Errors due to insufficient sampling and instrument noise are generally considered most important in EC flux measurements (Lenschow and Kristensen, 1985; Businger, 1986; Mauder et al., 2013; Rannik et al., 2016).

Sampling error is an inherent issue for EC flux measurements and is typically the main source of the CO2 flux uncertainty (Mauder et al., 2013). The sampling error is caused by the difference between the ensemble average and the time average. The calculation of EC flux (Eq. 2) requires the separation between the mean and fluctuating components, which can be represented fully for CO2 mixing ratio c as

(4) c ( x , t ) = c ( x , t ) + c ( x , t ) .

The mean component c represents ensemble average over time (t) and space (x) and does not contribute to the flux. The time average of a stationary turbulent signal and space average of a homogenous turbulent signal theoretically converge on the ensemble average when the averaging time approaches infinity, i.e. T→∞ (Wyngaard, 2010). In practice, Reynolds averaging over a much shorter time interval (10 min to an hour) is typically used for EC flux measurements from a fixed point or from a slow-moving platform such as a ship. This is because the atmospheric boundary layer is only quasi-stationary for a few hours. Non-stationarity (e.g. diurnal variability and synoptic conditions) is an inherent property of the atmospheric boundary layer (Wyngaard, 2010). EC flux observations thus inevitably contain some random error due to insufficient sampling time, and this error is greater at shorter averaging times.

Random error due to instrument noise comes mainly from the white noise of the gas analyser, as the noise from the sonic anemometer is relatively unimportant (Blomquist et al., 2010; Fairall et al., 2000; Mauder et al., 2013). Blomquist et al. (2014) show “pink” noise with a weak spectral slope for their CRDS gas analyser (G1301-f), but the gas analysers on JCR (G2311-f) and Discovery (LI-7200) demonstrate white noise with a constant variance at high frequency (Fig. B2, Appendix B).

Table 2Potential sources of bias in our EC air–sea CO2 flux measurements and the methods used to minimise them.

Download Print Version | Download XLSX

2.3.2 Systematic error

Table 2 details the measures taken during instrument setup and data processing that help eliminate most sources of systematic error in EC CO2 fluxes.

In addition to bias sources related to the instrument setup (Table 2), insufficient sampling time (an inherent issue of EC fluxes) may also generate a systematic error. We use a theoretical method to estimate this systematic error in EC CO2 flux (Lenschow et al., 1994):

(5) | δ F S | 2 σ w σ c a τ w τ c T ,

where σw (m s−1) and σca (ppm) are the standard deviations of the vertical wind velocity and the CO2 mixing ratio due to atmospheric processes, respectively. T is the averaging time interval (s), and τw and τc are integral timescales (s) for vertical wind velocity and CO2 signal, respectively. The definition and estimation of the integral timescale are shown in Appendix B. The sign of δFS could be positive or negative (i.e. under- or overestimation) because of the poor statistics in capturing low-frequency eddies within the flux averaging period (Lenschow et al., 1993). The mean hourly relative systematic error due to insufficient sampling time for four cruises estimated by Eq. (5) is < 5 %. According to propagation of uncertainty theory (JCGM, 2008), the total systematic error is less than 9 % (=7.5%2+5%2).

2.3.3 Random error

Five approaches used to estimate the total random error (A–C) and the random error component due to instrument noise (C–E) in EC CO2 fluxes are discussed below. The random error assessments are empirical (A and D) or theoretical (B, C, and E).

  • A.

    An empirical approach to estimate total random error involves shifting the w data relative to the CO2 data (or vice versa) by a large, unrealistic time shift and then computing the “null fluxes” from the time-desynchronized CO2 and w time series (Rannik et al., 2016). The shift removes any real correlation between CO2 and w due to vertical exchange. The standard deviation of the resultant null fluxes represents the random flux uncertainty (Wienhold et al., 1995). We applied a series of time shifts of  20–60 τw (i.e. using time shifts ranging from 300 to 100 and 100 to 300 s; Rannik et al., 2016). This empirical estimation of total random flux uncertainty will hereafter be referred to as δFR,Wienhold.

  • B.

    Lenschow and Kristensen (1985) derived a rigorous theoretical equation for total random error estimation, which contains both the auto-covariance and cross-covariance functions. The theoretical equation has been numerically approximated by Finkelstein and Sims (2001):

    (6) δ F R,Finkelstein = { 1 n p = - m m r ww ( p ) r cc ( p ) + p = - m m r wc ( p ) r cw ( p ) } 1 / 2 ,

    where n is the number of data points within an averaging time interval, and p is the number of shifting points. The maximum shifting point m can be chosen subjectively ( <n). We found that the random error for m between 1000 and 2000 data points was similar, so for this study we use m=1500 (150 s shift time). The first term in the brackets represents the auto-covariance component, and the second term is the cross-covariance component. rww and rcc are the auto-covariance functions for vertical wind velocity (w) and CO2 mixing ratio (c), respectively. rwc and rcw are the cross-covariance functions for w and c. Here rwc represents shifting w data relative to CO2 data, while rcw represents shifting CO2 data relative to w data.

  • C.

    Blomquist et al. (2010) attributed the sources of CO2 variance σc2 to atmospheric processes (σca2) and white noise (σcn2). The sources of variance are considered to be independent of each other, and the sonic anemometer is assumed to be relatively noise-free. According to propagation of uncertainty theory (JCGM, 2008), the total random flux error can be defined as

    (7) δ F R,Blomquist a σ w T σ c a 2 τ wc + σ c n 2 τ c n 1 / 2 ,

    where the constant a varies from 2 to 2, depending on the relationship between the covariance of the two variables (w and CO2) and the product of their auto-correlations (Lenschow and Kristensen, 1985). Here, τwc is equal to the shorter of τw and τc, which is typically τw (Blomquist et al., 2010), and τcn is the integral timescale of white noise in the CO2 signal. The CO2 variance due to atmospheric processes (σca2) includes two components: variance due to vertical flux (i.e. air–sea CO2 flux), σcav2, and variance due to other atmospheric processes, σcao2 (Fairall et al., 2000). The variance in CO2 due to vertical flux (σcav2) depends on atmospheric stability. σcav2 can be estimated with Monin–Obukhov similarity theory (Blomquist et al., 2010, 2014; Fairall et al., 2000):

    (8) σ c av 2 = 3 w c u f c ( z / L ) 2 ,

    where u is the friction velocity (m s−1), and the similarity function (fc) depends on the stability parameter z/L, where z is the observational height (m), and L is the Obukhov length (m). The expression of fc can be found in Blomquist et al. (2010).

    Equation (7) can be used to assess the random error due to instrument noise by setting σca2=0, referred to hereafter as δFRN,Blomquist. We use the CO2 variance spectra to directly estimate the white noise term σcn2τcn in Eq. (7). The variance is fairly constant at high frequency (1–5 Hz; Fig. B2, Appendix B), which is often referred to as band-limited white noise. The relationship between σcn2τcn and the band-limited noise spectral value φcn is expressed in Blomquist et al. (2010) as

    (9) σ c n 2 τ c n = φ c n 4 .
  • D.

    Billesbach (2011) developed an empirical method to estimate the random error due to instrument noise alone (referred to as ΔFRN,Billesbach). This involves random shuffling of the CO2 time series within an averaging interval and then calculating the covariance of w and CO2. The correlation between w and CO2 is minimised by the shuffling, and any remaining correlation between w and CO2 is due to the unintentional correlations contributed by instrument noise.

  • E.

    Mauder et al. (2013) describe another theoretical approach to estimate the random flux error due to instrument noise:

    (10) δ F RN,Mauder = σ w σ c n n .

    White noise correlates with itself but is uncorrelated with atmospheric turbulence. Thus, the white-noise-induced CO2 variance (σcn) only contributes to the total variance. The value of σcn can be estimated from the difference between the zero-shift auto-covariance value (CO2 variance σc2) and the noise-free variance extrapolated to a time shift of zero (Lenschow et al., 2000):

    (11) σ c n 2 = σ c 2 - σ 2 ( t 0 ) ,

    where σ2(t→0) represents the extrapolation of auto-covariance to a zero shift, which is considered equal to variance due to atmospheric processes (σca2). Figure 3 shows the normalised auto-covariance function curves of w and CO2 as measured by the Picarro G2311-f and the LI-7200. There is a sharp decrease in the CO2 auto-covariance when shifting from 0 s shift to 0.1 s shift for both the Picarro G2311-f and LI-7200 gas analyser. The same sharp decrease is not seen in the vertical wind velocity (w) auto-covariance. The relative difference in the change in normalised auto-covariance shows that white noise makes a much larger relative contribution to the CO2 variance than to the vertical wind velocity variance.

Figure 3Mean normalised auto-covariance functions of CO2 and vertical wind velocity (w) by four different instruments. The magenta line represents a fit to the noise-free auto-covariance function of CO2 (measured by Picarro) extrapolated back to a zero time shift. An example of the white noise and natural variability contributions to the total CO2 (measured by Picarro) variance is indicated by two blue arrows. The sharp decrease of the CO2 auto-covariance between the zero shift and the initial 0.1 s shift corresponds to the large contribution of white noise from the gas analysers. The LI-7200 is the noisier instrument. The noise contributions from the anemometer are relatively small (< 10 %).


Figure 4(a) Air–sea CO2 fluxes (hourly and 6 h averages), (b) random uncertainty in flux (total and due to instrument noise only), and (c) variance in CO2 mixing ratio (total and due to instrument noise only) for two Atlantic cruises.


3 Results

Measurements from AMT28 and AMT29 set the scene for our uncertainty analysis. These two Atlantic cruises transited across the same tropical region (Fig. A2, Appendix A) in October 2018 and September 2019 with different eddy covariance systems (Sect. 2.1). AMT28 and AMT29 show broadly similar latitudinal patterns (Fig. 4a). An obvious question of interest is whether the measured fluxes were the same for the 2 years. To answer this question, the measurement uncertainties must be quantified. The total random uncertainties in CO2 flux (δFR,Finkelstein) are comparable for the two cruises, even though the random error component due to instrument noise (δFRN,Mauder) is about 3 times higher during AMT29 using LI-7200 than during AMT28 using Picarro G2311-f (Fig. 4b; Fig. D1, Appendix D). The similar total random uncertainty in the AMT28 and AMT29 fluxes shows that both gas analysers are equally suitable for air–sea EC CO2 flux measurements. The variance budgets of atmospheric CO2 mixing ratio (used to estimate random flux uncertainty; see Sect. 3.1) are shown in Fig. 4c. Total variance in CO2 mixing ratio is dominated by instrument noise on both cruises. CO2 mixing ratio variance (total and instrument noise) was substantially higher during AMT29.

Table 3Variance in the CO2 mixing ratio estimated using Eqs. (8) and (11) for the Arctic (JR18006/7, Picarro G2311-f) and Atlantic cruises (AMT28, Picarro G2311-f; AMT29, LI-7200). Total CO2 variance (σc2) consists of white noise (σcn2) and atmospheric processes (σca2). The latter can be further broken down to the CO2 variance due to vertical flux (σcav2) and due to other processes (σcao2).

Download Print Version | Download XLSX

3.1 Random uncertainty

Theoretical derivation of flux uncertainty (δFRN,Blomquist, Eq. 7) requires knowledge of the contributions to CO2 mixing ratio variance. Total CO2 variance is made up of instrument noise (σcn2) and atmospheric processes (σca2). Atmospheric processes include vertical flux (σcav2) and other atmospheric processes (σcao2). The variance budgets of CO2 mixing ratio for the four cruises are listed in Table 3. Atmospheric processes contribute a larger CO2 variance in the Arctic (where flux magnitudes are greater) compared to the Atlantic. Vertical flux accounts for  10 % of the variance in CO2 mixing ratio in the Arctic and  1 % of the CO2 variance in the Atlantic. Previous results demonstrate that horizontal transport is a major source of σcao2 for long-lived greenhouse gases (Blomquist et al., 2012). Small changes in CO2 mixing ratio transported horizontally can yield variance that greatly exceeds the variance from vertical flux.

Three quasi-independent methods were used to estimate random uncertainty in EC air–sea CO2 fluxes caused by instrument noise (δFRN; Methods C–E, Sect. 2.3.3). Good agreement was found between all three estimates (Fig. C2, Appendix C) when 2 is used as the constant in Eq. (7) (a). The ΔFRN,Billesbach estimates have more scatter and are slightly higher than the theoretical results, possibly because the random shuffling of data fails to fully exclude the contribution from atmospheric turbulence (Rannik et al., 2016). For the remainder of this study, we use the δFRN,Mauder method to estimate δFRN.

We used three methods to estimate the total random uncertainty (δFR; Methods A–C, Sect. 2.3.3) in the hourly averaged air–sea CO2 fluxes. There is good agreement among the three estimates (r> 0.88; Fig. C1, Appendix C). Again, the constant in Eq. (7) (a) is set to 2, as informed by the instrument noise uncertainty analysis above. We use δFR,Finkelstein (Eq. 6) to estimate the total random flux uncertainty hereafter. Our decision is based on δFR,Finkelstein not requiring the integral timescale (unlike δFR,Blomquist) and showing less scatter than δFR,Wienhold.

Figure 5Relative random uncertainty in hourly CO2 flux and its contribution from noise, vertical flux, and other processes during two Arctic cruises. Relative random uncertainty data are binned into 3 mmolm-2d-1 flux magnitude bins (error bars represent 1 standard deviation).


Figure 5 shows the different relative contributions to the random flux uncertainty for the Arctic cruises (hourly average). Here the uncertainty is normalised by the flux magnitude and then averaged into flux magnitude bins. When the flux magnitude is sufficiently large (> 20 mmolm-2d-1), the total relative random uncertainty in flux asymptotes to about 15 % and is driven by variance associated with both vertical flux and other atmospheric processes. This estimate is similar to uncertainties in air–sea fluxes of other well resolved (i.e. high signal-to-noise ratio) variables (Fairall et al., 2000). At a lower flux magnitude, uncertainty due to atmospheric processes other than vertical flux dominates the total random uncertainty. Uncertainty due to the white noise from the Picarro G2311-f gas analyser is small.

Table 4Summary of hourly average EC CO2 fluxes and associated uncertainties in the mean for the four cruises (mmolm-2d-1). Shown are the mean CO2 flux magnitude (|F|, mmolm-2d-1), upper limitation of the total uncertainty (δF; Eq. 3), upper limitation of the absolute systematic error (|δFS|; propagated from Table 2 and Eq. 5), and random error (δFR; Eq. 6). The random error components are white noise (δFRN; Eq. 10), vertical flux (δFRV; Eqs. 7 and 8), and other atmospheric processes (δFRO=δFR2-δFRN2-δFRV2). The total uncertainty is also expressed as a percentage of the mean flux magnitude (δF/|F|100%).

Download Print Version | Download XLSX

3.2 Summary of systematic and random uncertainties

The total uncertainty δF in the hourly average EC CO2 flux (estimated using Eq. 3) ranges from 1.4 to 3.2 mmolm-2d-1 in the mean for the four cruises (Table 4). Our EC flux system setup was optimal, and subsequent corrections have minimised any bias to < 9 % (Sect. 2.3.2). Systematic error is on average much lower than random error (Table 4). This means the accuracy of the EC CO2 flux measurements is very high, but the precision of hourly averaged EC CO2 air–sea flux measurements is relatively low. In Sect. 4.1, we discuss how the precision can be improved by averaging the observed fluxes for longer.

The theoretical uncertainty estimates above can be compared with a portion of the AMT28 cruise data (15–20 S,  25 W; Fig. 4), when the ship encountered sea surface CO2 fugacity close to equilibrium with the atmosphere (i.e. ΔfCO2 0; Fig. A2, Appendix A). The data from this region are useful for assessing the random and systematic flux uncertainties. The standard deviation of the EC CO2 flux during cruise AMT28 when ΔfCO2 0 is 1.6 mmolm-2d-1, which compares well with the theoretical random flux uncertainty in this region (1.4 mmolm-2d-1). The mean EC CO2 flux from this region was 0.5 mmolm-2d-1, which is indistinguishable from zero considering the random uncertainty. This further confirms the minimal bias in our flux observations.

Figure 6Comparison of relative random uncertainty in hourly CO2 flux and relative standard deviation (RSTD; standard  deviation/|flux  mean|) of the EC CO2 flux from two Arctic cruises. These results are binned in 1 m s−1 wind speed bins.


Figure 6 shows a comparison between the relative uncertainty and the relative standard deviation (RSTD) in the hourly CO2 flux for the two Arctic cruises. Results have been binned into 1 m s−1 wind speed bins. Wind speed was converted to 10 m neutral wind speed (U10N) using the COARE3.5 model (Edson et al., 2013). The relative random error decreases with increasing wind speed. This is partly because the fluxes tend to be larger at higher wind speeds, and so the signal-to-noise ratio in the flux is greater. In addition, at higher wind speeds, a greater number of high-frequency turbulent eddies are sampled by the EC system, providing better statistics of turbulent eddies and lower sampling error.

The RSTD of the flux is greater in magnitude than the estimated flux uncertainty because it also contains environmental variability. The CO2 flux auto-covariance analysis (Sect. 4.1) shows that random error in hourly flux explains  20 % of the flux variance on average for the two Arctic cruises. This implies that the remaining variability in the EC flux ( 80 %) is due to natural phenomena (e.g. changes in ΔfCO2 or wind speed). Similarly, substantial variability is typical in EC-derived CO2 gas transfer velocity at a given wind speed (e.g. Edson et al., 2011; Butterworth and Miller, 2016). K660 is derived from (ECCO2flux)/ΔfCO2, and thus an understanding of EC flux uncertainty can help understand and explain the variability in EC-derived gas transfer velocity estimates (Sect. 4.2).

4 Discussion

4.1 Impact of averaging timescale on flux uncertainty

The random error in flux decreases with increasing averaging time interval T or the number of sampling points n (Eqs. 6, 7, and 10). This is because a longer averaging time interval results in better statistics of the turbulent eddies. However, averaging for too long is also not ideal since the atmosphere is less likely to maintain stationarity. The typical averaging time interval is thus typically between 10 min and 60 min for air–sea flux measurements (20 min intervals were used in this study). The time series of quality controlled 20 min flux intervals can be further averaged over a longer timescale to reduce the random uncertainty. Averaging the 20 min flux intervals assumes that the flux interval data are essentially repeat measurements within a chosen averaging timescale. If the 20 min flux intervals are averaged, one can ask the following question: what is the optimal averaging timescale for interpreting air–sea EC CO2 fluxes?

Figure 7(a) Auto-covariance of the original 20 min fluxes (cruise JR18007) and a fit to the noise-free auto-covariance function extrapolated back to a zero time shift. (b) CO2 flux auto-covariance functions with different averaging timescales. The black line represents the auto-covariance of the original 20 min fluxes. The 20 min fluxes are further averaged at different timescales (1, 2, 3, and 6 h), and the corresponding auto-covariance functions are shown with different colours (dark blue, orange, green, and light blue).


We use an auto-covariance method to determine the optimal averaging timescale. The observed variance in CO2 flux consists of random uncertainty (random noise) as well as natural variability. The random noise component should only contribute to the CO2 flux variance when the data are zero-shifted. After the CO2 flux data are shifted, the noise will not contribute to the auto-covariance function. Figure 7 shows the auto-covariance function of the air–sea CO2 flux with different averaging timescales for Arctic cruise JR18007. For the 20 min fluxes (Fig. 7a), the auto-covariance decreases rapidly between the zero shift and the initial time shift, which indicates that a large fraction of the 20 min flux variance is due to random noise.

Figure 8Effect of the averaging timescale on the noise:signal (random  noise  in  flux/natural  flux  variability) for EC air–sea CO2 flux measurements during four cruises.


The random noise in the CO2 fluxes decreases with a longer averaging timescale, with the greatest effect observed from 20 min to 1 h (Fig. 7b). A fit to the noise-free auto-covariance function extrapolated back to a zero time shift gives us an estimate of the non-noise variability in the natural CO2 flux. Subtracting the extrapolated natural flux variability from the total variance in CO2 flux provides an estimate of the random noise in the flux for each averaging timescale (Fig. 7a). All four cruises consistently demonstrate a non-linear reduction in the noise contribution to the flux measurements when the averaging timescale increases (Fig. 8). The random noise in flux can be expressed relative to the natural variance in flux representing the inverse of the signal-to-noise ratio (i.e. random  noise  in  flux/natural  flux  variability, hereafter referred to as noise:signal).

The noise:signal also facilitates comparison of all four cruises (Fig. 8) and demonstrates the consistent effect that increasing the averaging timescale has on noise:signal. Consistent with Table 4, the Arctic cruises show much lower noise:signal because the flux magnitudes are much larger. Typical detection limits in analytical science are often defined by a 1:3 noise:signal ratio. A 1:3 noise:signal is achieved with a 1 h averaging timescale for the Arctic cruises. The Atlantic cruises encountered much lower air–sea CO2 fluxes, and an averaging timescale of at least 3 h is required to achieve the same 1:3 noise:signal ratio.

The flux measurement uncertainty at a 6 h averaging timescale for the AMT cruises is  0.6 mmolm-2d-1. The analysis presented above permits an answer to the question posed at the beginning of the Results section. The mean difference between the 6 h averaged EC CO2 flux observations on AMT29 and AMT28 (1.3 mmolm-2d-1; Fig. 4a) is much greater than the measurement uncertainty. This significant difference was likely because of the interannual variability in AMT CO2 flux due to changes in the natural environment (e.g. ΔfCO2, sea surface temperature, and physical drivers of interfacial turbulence such as wind speed) during the two cruises.

At a typical research ship speed of  10 knots, the AMT cruises cover  110 km in 6 h, which is equivalent to  1 latitude. Averaging for longer than 6 h is likely to cause substantial loss of real information about the natural variations in air–sea CO2 flux and the drivers of flux variability. For example, the mean flux between 0–20 S during cruise AMT28 is 0.9 mmolm-2d-1. However, the 6 h average EC measurements show that the flux varied between +5 mmolm-2d-1 ( 2–6 S) and 5 mmolm-2d-1 ( 11–13 S; Fig. 4a).

4.2 Effect of CO2 flux uncertainty on the gas transfer velocity K

The uncertainties in the EC CO2 air–sea flux measurement will influence the uncertainty that translates to EC-based estimates of the gas transfer velocity, K. For illustration, K is computed for Arctic cruise JR18007, which had a high flux signal:noise ratio of  5 (Fig. 8). Any data potentially influenced by ice and sea ice melt were excluded using a sea surface salinity filter (data excluded when salinity < 32). Equation (1) is rearranged and used with concurrent measurements of CO2 flux (F), ΔfCO2, and sea surface temperature (SST) to obtain K adjusted for the effect of temperature (K660).

Figure 9Gas transfer velocity (K660) measured on Arctic cruise JR18007 (hourly average, signal:noise  5) vs. 10 m neutral wind speed (U10N). Red squares represent 1 m s−1 bin averages, with error bars representing 1 standard deviation (SD). The red curve represents a quadratic fit using the bin averages: K660= 0.22 U10N2+ 2.46 (R2= 0.76). The grey shaded area represents the standard deviation calculated for each wind speed bin (K660± 1 SD). The cyan region represents the upper and lower bounds in K660 uncertainty computed from the EC flux uncertainty (K660±δK660; see text for detail).


The determination coefficient (R2) of the quadratic fit between wind speed (U10N) and EC-derived K660 (Fig. 9) demonstrates that wind speed explains 76 % of the K660 variance during Arctic cruise JR18007. How much of the remaining 24 % can be attributed to uncertainties in EC CO2 fluxes?

Variability in K660 within each 1 m s−1 wind speed bin can be considered to have minimal wind speed influence. It is thus useful to compare the variability within each wind speed bin (K660± 1 SD) with the upper and lower uncertainty bounds derived from the EC flux measurements. Uncertainty in EC flux-derived K660 (δK660) is calculated from the uncertainty in hourly EC flux (δF) by rearranging Eq. (1) (bulk flux equation) and replacing F with δF. The resultant δK660 is then averaged in wind speed bins. The shaded cyan band in Fig. 9 (K660±δK660) is consistently narrower than the grey shaded band (K660± 1 SD). On average, EC flux-derived uncertainty in K660 can only account for a quarter of the K660 variance within each wind speed bin, and the remaining variance is most likely due to the non-wind speed factors that influence gas exchange (e.g. breaking waves, surfactants).

The analysis above can be extended to assess how EC flux-derived uncertainty affects our ability to parameterise K660 (e.g. as a function of wind speed). To do so, a set of synthetic K660 data is generated (same U10N as the K660 measurements in Fig. 9). The synthetic K660 data are initialised using a quadratic wind speed dependence that matches JR18007 (i.e. K660= 0.22 U10N2+ 2.46). Random Gaussian noise is then added to the synthetic K660 data, with relative noise level corresponding to the relative flux uncertainty values taken from JR18007 (mean of 20 %; Table 4). The relative uncertainty in K660 due to EC flux uncertainty (δK660/K660) shows a wind speed dependence (Fig. S4a in the Supplement), and the artificially generated Gaussian noise incorporates this wind speed dependence (Fig. S4b, Supplement). The R2 of the quadratic fit to the synthetic data as a function of U10N is 0.90 (the rest of the variance is due to uncertainty in K660). Since wind speed explains 76 % of variance in the observed K660, it can be inferred that non-wind speed factors can account for 14 % (i.e. (100–76) %–(100–90) %) of the total variance in K660 from this Arctic cruise. If the synthetic K660 data are assigned a relative flux uncertainty of 50 % (reflective of a region with low fluxes, e.g. AMT28/29), the R2 of the wind speed dependence in the synthetic data decreases to 0.60.

Figure 10Relative uncertainty in EC-estimated hourly K660 (δK660/K660) vs. the magnitude of the air–sea CO2 fugacity difference (fCO2|) during Arctic cruise JR18007 and Atlantic cruises AMT28 and AMT29 (no ΔfCO2 data were collected on JR18006). The data points are colour-coded by wind speed. Blue points are medians of δK660/K660 in 5 µatm bins. Here we use the parameterised K660 (= 0.22 U10N2+ 2.46) to normalise the uncertainty in K660. The dashed line represents the 3:1 signal:noise ratio (δK660/K660=1/3).


The relative uncertainty in EC flux-derived K660 (δK660/K660) is large when fCO2| is small (Fig. 10). Previous EC studies have filtered EC flux data to remove fluxes when the fCO2| falls below a specified threshold (e.g. 20 µatm, Blomquist et al. (2017); 40 µatm, Miller et al. (2010), Landwehr et al. (2014), Butterworth and Miller (2016), Prytherch et al. (2017); 50 µatm, Landwehr et al. (2018)). Analysis of the data presented here suggests that a fCO2| threshold of at least 20 µatm is reasonable for hourly K660 measurements, leading to δK660 of  10 cm h−1 (δK660/K6601/3) or less on average. At very large fCO2| of over 100 µatm, δK660 is reduced to only a few centimetres per hour (cm h−1) (δK660/K660 1/5). At longer flux averaging timescales, it may be possible to relax the minimal fCO2| threshold.

5 Conclusions

This study uses data from four cruises with a range in air–sea CO2 flux magnitude to comprehensively assess the sources of uncertainty in EC air–sea CO2 flux measurements. Data from two ships and two different state-of-the-art CO2 analysers (Picarro G2311-f and LI-7200, both fitted with a dryer) are analysed using multiple methods (Sect. 2.3). Random error accounts for the majority of the flux uncertainty, while the systematic error (bias) is small (Table 4). Random flux uncertainty is primarily caused by variance in CO2 mixing ratio due to atmospheric processes. The random error due to instrument noise for the Picarro G2311-f is 3-fold smaller than for LI-7200 (Table 4 and Fig. D1, Appendix D). However, the contribution of the instrument noise to the total random uncertainty is much smaller than the contribution of atmospheric processes such that both gas analysers are well suited for air–sea CO2 flux measurements.

The mean uncertainty in hourly EC flux is estimated to be 1.4–3.2 mmolm-2d-1, which equates to the relative uncertainty of  20 % in high CO2 flux regions and  50 % in low CO2 flux regions. Lengthening the averaging timescale can improve the signal:noise ratio in EC CO2 flux through the reduction of random uncertainty. Auto-covariance analysis of CO2 flux is used to quantify the optimal averaging timescale (Figs. 7 and 8, Sect. 4.1). The optimal averaging timescale varies between 1 h for regions of large CO2 flux (Arctic in our analysis) and at least 3 h for regions of low CO2 flux (tropical/subtropical Atlantic in our analysis).

The measurement uncertainty in EC CO2 flux contributes directly to scatter in the derived gas transfer velocity, K660. Flux uncertainties determined in this paper are applied to a synthetic K660 dataset. This enables a partitioning of the variance in measured K660 that is due to EC CO2 flux uncertainty, wind speed, and other processes (10 %, 76 %, 14 % for Arctic cruise JR18007). At a given averaging timescale, a fCO2| threshold helps to reduce the scatter in K660. A minimum fCO2| filter of 20 µatm is needed for interpreting hourly K660 data, with the signal:noise ratio in K660 improving further at higher fCO2|.

Appendix A: Cruise tracks

Figure A1Cruise tracks of JR18006 (magenta) and JR18007 (green). The bottom colour bar indicates the CO2 fugacity difference (ΔfCO2) of August 2019 (Bakker et al., 2016; Landschützer et al., 2020), while the right colour bar shows the Arctic sea ice concentrations of 1 August 2019 measured by Advanced Microwave Scanning Radiometer – Earth Observing System Sensor (AMSR-E; Spreen et al., 2008).

Figure A2Cruise tracks of AMT28 (magenta) and AMT29 (green). The ocean is coloured with the ΔfCO2 for October 2018 (Bakker et al., 2016; Landschützer et al., 2020).


Appendix B: Integral timescale and variance spectra of CO2 and vertical wind velocity

Integral timescale is used in the flux uncertainty calculation (Eqs. 5 and 7). The definition of integral timescale τx of variable x is

(B1) τ x = 1 σ x 2 0 r x x ( t ) d t ,

where σx2 is the variance of x, and rxx is the auto-covariance function of x. t is the shifting time of auto-covariance (which is different from the lag time between w and CO2 in the EC flux calculation). We can use Eq. (B1) to estimate the integral timescale of w and CO2 directly. However, integration up to infinity is not practical. Instead we can numerically estimate the timescale by determining the time corresponding to the auto-covariance coefficient function (rxx/σx2) value decaying to 1/e (1/e decaying method) or by integrating the auto-covariance function up to the first zero crossing of the function (zero crossing method) (Rannik et al., 2009).

One can also use similarity theory to estimate the integral timescale theoretically (Blomquist et al., 2010):

(B2) τ w = 2.8 z u r f τ ( z / L ) .

Here, ur is the relative wind speed. The similarity function fτ(z/L) is described by the stability parameter z/L, where z is the observation height (m), and L is the Obukhov length (m) (Blomquist et al., 2010).

Yet another method to estimate the integral timescale is from the peak frequency (fmax) in the w variance spectrum (Kaimal and Finnigan, 1994):

(B3) τ w = 1 2 π f max .

Figure B1Comparison of integral timescales of w estimated by four different methods. Estimated integral timescales from the zero crossing method (integrating the auto-covariance function up to first zero crossing the function) agree well with the estimation of peak frequency method (Eq. B3). However, the similarity method (Eq. B2) overestimates the integral timescale, whereas the 1/e decaying method (determining the time needed for the auto-covariance coefficient function value to decay to 1/e) tends to underestimate the integral timescale.


Figure B2Mean variance spectra for CO2 and w for one Arctic cruise JR18007. The near-constant CO2 variance at high frequency (1–5 Hz) indicates the band-limited noise in the CO2 signal. In contrast, the w spectrum does not show a similar band-limited noise at < 10 Hz.


The integral timescales of w estimated by these four methods for cruise JR18007 are shown in Fig. B1. The integral timescale estimated by the zero crossing method agrees well with the peak frequency estimates using Eq. (B3). The 1/e decaying method tends to underestimate the integral timescale, which is generally observed for turbulent signals (Rannik et al., 2009), whereas the similarity method (Eq. B2) considerably overestimates the integral timescale. Based on the recent analysis (as yet unpublished) of the entire NOAA PSL flux database, the Eq. (B2) formulation is now thought to be an overestimate (review comment for this paper from Blomquist, 2021). In this study we use the integral timescale of w from the zero crossing method to estimate the theoretical flux uncertainty (Eqs. 5 and 7). The theoretical systematic error estimates (Eq. 8) also require the integral timescale of CO2. The integral timescale of CO2 is difficult to evaluate from the above four methods due to instrument noise. Instead, we estimate it by directly integrating the auto-covariance function (Eq. B1) to a shift time of 200 s (we found no significant difference of the integral timescale when integrating the CO2 auto-covariance function for shift times ranging from 150 to 250 s).

Appendix C: Comparison of the uncertainty estimates by different methods

Figure C1Comparison of total random uncertainties in hourly flux estimated by three different methods for the Arctic cruises. The empirical estimates FR,Wienhold agree well with one of the theoretical estimates ΔFR,Finkelstein (r=0.93). The other theoretical estimate ΔFR,Blomquist is slightly higher than the random uncertainties ΔFR,Finkelstein (slope = 1.13) if the constant in Eq. (8) is set equal to 2.


Figure C2Comparison of random error in hourly flux due to instrument white noise, estimated by three different methods for the Arctic cruises. The three uncertainty estimations agree well. The correlation coefficient (r) between δFRN,Mauder and δFRN,Blomquist is 1 if the constant in Eq. (7) (a) is set to 2.


Appendix D: Performance of two gas analysers

Figure D1 shows a comparison between the performance of the Picarro 2311-f and the LI-7200 gas analysers. We estimated that the noise of the LI-7200 is on average 3 times higher than that of the Picarro 2311-f (Table 3). Indeed, random error in the CO2 flux due to the white noise is much higher for the LI-7200 than for the Picarro 2311-f, but the total flux uncertainty of the EC system with the LI-7200 on AMT29 is only slightly higher than that of the EC system with the Picarro 2311-f on AMT28 (Table 4). Again, this is because for both EC systems, sampling error dominates the total random uncertainty, while the contribution of instrument noise (< 30 %) to the total uncertainty is relatively small (Billesbach, 2011; Langford et al., 2015; Mauder et al., 2013; Rannik et al., 2016). Another often used CRDS gas analyser in EC measurements is the Los Gatos Research (LGR) Fast Greenhouse Gas Analyser (FGGA) (Prytherch et al., 2017). Yang et al. (2016) showed that LGR FGGA is ca. 10 times noisier than the Picarro G2311-f, and as a result, the total CO2 flux uncertainty measured by the LGR is 4 times higher than that by the Picarro. From the perspective of measurement noise, Picarro and LI-7200 gas analysers are better suited for air–sea CO2 flux measurements than the LGR FGGA.

Figure D1Comparison of the relative total random uncertainty and the relative random error component due to white noise for different gas analysers. A Picarro G2311-f gas analyser was used on AMT28 and a LI-7200 infrared gas analyser on AMT29.


Data availability

The processed hourly EC CO2 fluxes and uncertainties can be found in the Supplement of this paper. Raw, high-frequency (10 Hz) data are large (tens of gigabytes) and are archived at PML. Please contact the authors directly if you are interested in the raw data.


The supplement related to this article is available online at:

Author contributions

TGB and MY designed and installed the eddy covariance systems on ships and managed the collections of measurements. VK collected and processed the CO2 fugacity data. YD processed and analysed the data with the help of MY and TGB. YD wrote the paper with input from DCEB, TGB, and MY. All authors contributed to and approved the final paper.

Competing interests

The authors declare that they have no conflict of interest.


We thank captains and crew of the RRS James Clark Ross and RRS Discovery and all those who have helped keep the CO2 flux systems running. We are extremely grateful to Brian J. Butterworth (University of Calgary) for his advice on how to set up and run the automated CO2 flux system on JCR and how to code the CR6 data logger, as well as to Tim J. Smyth (PML) for setting up the remote monitoring of flux data. We also greatly appreciate Ian Brown (PML) and Daniel Phillips for fCO2 measurements and Peter S. Liss (UEA) for support and helpful comments.

Financial support

This work is funded by the China Scholarship Council (CSC/201906330072). Air–sea CO2 flux measurements were facilitated by the European Space Agency (ESA AMT4oceanSatFlux project, grant no. 4000125730/18/NL/FF/gp) and support from the Natural Environment Research Council (NERC) via PML's contribution to the ORCHESTRA programme (NE/N018095/1). The Arctic cruises were also supported by NERC, through the DIAPOD (NE/P006280/2) and ChAOS (NE/P006493/1) projects.

Review statement

This paper was edited by Leiming Zhang and reviewed by Byron Blomquist and one anonymous referee.


Bakker, D. C. E., Pfeil, B., Landa, C. S., Metzl, N., O'Brien, K. M., Olsen, A., Smith, K., Cosca, C., Harasawa, S., Jones, S. D., Nakaoka, S., Nojiri, Y., Schuster, U., Steinhoff, T., Sweeney, C., Takahashi, T., Tilbrook, B., Wada, C., Wanninkhof, R., Alin, S. R., Balestrini, C. F., Barbero, L., Bates, N. R., Bianchi, A. A., Bonou, F., Boutin, J., Bozec, Y., Burger, E. F., Cai, W.-J., Castle, R. D., Chen, L., Chierici, M., Currie, K., Evans, W., Featherstone, C., Feely, R. A., Fransson, A., Goyet, C., Greenwood, N., Gregor, L., Hankin, S., Hardman-Mountford, N. J., Harlay, J., Hauck, J., Hoppema, M., Humphreys, M. P., Hunt, C. W., Huss, B., Ibánhez, J. S. P., Johannessen, T., Keeling, R., Kitidis, V., Körtzinger, A., Kozyr, A., Krasakopoulou, E., Kuwata, A., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lo Monaco, C., Manke, A., Mathis, J. T., Merlivat, L., Millero, F. J., Monteiro, P. M. S., Munro, D. R., Murata, A., Newberger, T., Omar, A. M., Ono, T., Paterson, K., Pearce, D., Pierrot, D., Robbins, L. L., Saito, S., Salisbury, J., Schlitzer, R., Schneider, B., Schweitzer, R., Sieger, R., Skjelvan, I., Sullivan, K. F., Sutherland, S. C., Sutton, A. J., Tadokoro, K., Telszewski, M., Tuma, M., van Heuven, S. M. A. C., Vandemark, D., Ward, B., Watson, A. J., and Xu, S.: A multi-decade record of high-quality fCO2 data in version 3 of the Surface Ocean CO2 Atlas (SOCAT), Earth Syst. Sci. Data, 8, 383–413,, 2016. 

Baldocchi, D., Falge, E., Gu, L., Olson, R., Hollinger, D., Running, S., Anthoni, P., Bernhofer, C., Davis, K., and Evans, R.: FLUXNET: A new tool to study the temporal and spatial variability of ecosystem-scale carbon dioxide, water vapor, and energy flux densities, B. Am. Meteorol. Soc., 82, 2415–2434, 2001. 

Bariteau, L., Helmig, D., Fairall, C. W., Hare, J. E., Hueber, J., and Lang, E. K.: Determination of oceanic ozone deposition by ship-borne eddy covariance flux measurements, Atmos. Meas. Tech., 3, 441–455,, 2010. 

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

Billesbach, D. P.: Estimating uncertainties in individual eddy covariance flux measurements: A comparison of methods and a proposed new method, Agr. Forest Meteorol., 151, 394–405,, 2011. 

Blomquist, B. W., Huebert, B. J., Fairall, C. W., and Faloona, I. C.: Determining the sea–air flux of dimethylsulfide by eddy correlation using mass spectrometry, Atmos. Meas. Tech., 3, 1–20,, 2010. 

Blomquist, B.: Interactive comment on “Uncertainties in eddy covariance air-sea CO2 flux measurements and implications for gas transfer velocity parameterisations” by Yuanxu Dong et al., Atmos. Chem. Phys. Discuss.,, 2021. 

Blomquist, B. W., Fairall, C. W., Huebert, B. J., and Wilson, S. T.: Direct measurement of the oceanic carbon monoxide flux by eddy correlation, Atmos. Meas. Tech., 5, 3069–3075,, 2012. 

Blomquist, B. W., Huebert, B. J., Fairall, C. W., Bariteau, L., Edson, J. B., Hare, J. E., and McGillis, W. R.: Advances in Air-Sea CO2 Flux Measurement by Eddy Correlation, Bound.-Lay. Meteorol., 152, 245–276,, 2014. 

Blomquist, B. W., Brumer, S. E., Fairall, C. W., Huebert, B. J., Zappa, C. J., Brooks, I. M., Yang, M., Bariteau, L., Prytherch, J., Hare, J. E., Czerski, H., Matei, A., and Pascal, R. W.: Wind Speed and Sea State Dependencies of Air-Sea Gas Transfer: Results From the High Wind Speed Gas Exchange Study (HiWinGS), J. Geophys. Res.-Oceans, 122, 8034–8062,, 2017. 

Broecker, W. S. and Peng, T. H.: Greenhouse Puzzles Part 1 Keeling's World: is CO2 Greening the Earth, k1–k112, Columbia University, Lamont-Doherty Earth Observatory, 1998. 

Broecker, W. S., Ledwell, J. R., Takahashi, T., Weiss, R., Merlivat, L., Memery, L., Peng, T., Jahne, B., and Munnich, K. O.: Isotopic versus micrometeorologic ocean CO2 fluxes: A serious conflict, J. Geophys. Res.-Oceans, 91, 10517–10527, 1986. 

Businger, J.: Evaluation of the accuracy with which dry deposition can be measured with current micrometeorological techniques, J. Clim. Appl. Meteorol., 25, 1100–1124, 1986. 

Butterworth, B. J. and Miller, S. D.: Air–sea exchange of carbon dioxide in the Southern Ocean and Antarctic marginal ice zone, Geophys. Res. Lett., 43, 7223–7230,, 2016. 

Edson, J. B., Hinton, A. A., Prada, K. E., Hare, J. E., and Fairall, C. W.: Direct covariance flux estimates from mobile platforms at sea, J. Atmos. Ocean. Tech., 15, 547–562,<0547:DCFEFM>2.0.CO;2, 1998. 

Edson, J. B., Fairall, C. W., Bariteau, L., Zappa, C. J., Cifuentes-Lorenzen, A., McGillis, W. R., Pezoa, S., Hare, J. E., and Helmig, D.: Direct covariance measurement of CO2 gas transfer velocity during the 2008 Southern Ocean Gas Exchange Experiment: Wind speed dependency, J. Geophys. Res.-Oceans, 116, C00F10,, 2011. 

Edson, J. B., Jampana, V., Weller, R. A., Bigorre, S. P., Plueddemann, A. J., Fairall, C. W., Miller, S. D., Mahrt, L., Vickers, D., and Hersbach, H.: On the exchange of momentum over the open ocean, J. Phys. Oceanogr., 43, 1589–1610, 2013. 

Else, B. G. T., Papakyriakou, T. N., Galley, R. J., Drennan, W. M., Miller, L. A., and Thomas, H.: Wintertime CO2 fluxes in an Arctic polynya using eddy covariance: Evidence for enhanced air–sea gas transfer during ice formation, J. Geophys. Res.-Oceans, 116, C00G03,, 2011. 

Fairall, C. W., Hare, J. E., Edson, J. B., and McGillis, W.: Parameterization and micrometeorological measurement of air–sea gas transfer, Bound.-Lay. Meteorol., 96, 63–105,, 2000. 

Finkelstein, P. L. and Sims, P. F.: Sampling error in eddy correlation flux measurements, J. Geophys. Res.-Atmos., 106, 3503–3509, 2001. 

Flügge, M., Paskyabi, M. B., Reuder, J., Edson, J. B., and Plueddemann, A. J.: Comparison of direct covariance flux measurements from an offshore tower and a buoy, J. Atmos. Ocean. Tech., 33, 873–890,, 2016. 

Friedlingstein, P., O'Sullivan, M., Jones, M. W., Andrew, R. M., Hauck, J., Olsen, A., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S., Aragão, L. E. O. C., Arneth, A., Arora, V., Bates, N. R., Becker, M., Benoit-Cattin, A., Bittig, H. C., Bopp, L., Bultan, S., Chandra, N., Chevallier, F., Chini, L. P., Evans, W., Florentie, L., Forster, P. M., Gasser, T., Gehlen, M., Gilfillan, D., Gkritzalis, T., Gregor, L., Gruber, N., Harris, I., Hartung, K., Haverd, V., Houghton, R. A., Ilyina, T., Jain, A. K., Joetzjer, E., Kadono, K., Kato, E., Kitidis, V., Korsbakken, J. I., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Liu, Z., Lombardozzi, D., Marland, G., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Niwa, Y., O'Brien, K., Ono, T., Palmer, P. I., Pierrot, D., Poulter, B., Resplandy, L., Robertson, E., Rödenbeck, C., Schwinger, J., Séférian, R., Skjelvan, I., Smith, A. J. P., Sutton, A. J., Tanhua, T., Tans, P. P., Tian, H., Tilbrook, B., van der Werf, G., Vuichard, N., Walker, A. P., Wanninkhof, R., Watson, A. J., Willis, D., Wiltshire, A. J., Yuan, W., Yue, X., and Zaehle, S.: Global Carbon Budget 2020, Earth Syst. Sci. Data, 12, 3269–3340,, 2020. 

Garbe, C. S., Rutgersson, A., Boutin, J., De Leeuw, G., Delille, B., Fairall, C. W., Gruber, N., Hare, J., Ho, D. T., and Johnson, M. T.: Transfer across the air–sea interface, in: Ocean-atmosphere interactions of gases and particles, Springer, Berlin, Heidelberg, pp. 55–112, 2014. 

Ho, D. T., Law, C. S., Smith, M. J., Schlosser, P., Harvey, M., and Hill, P.: Measurements of air–sea gas exchange at high wind speeds in the Southern Ocean: Implications for global parameterizations, Geophys. Res. Lett., 33, L16611,, 2006. 

Ikawa, H., Faloona, I., Kochendorfer, J., Paw U, K. T., and Oechel, W. C.: Air–sea exchange of CO2 at a Northern California coastal site along the California Current upwelling system, Biogeosciences, 10, 4419–4432,, 2013. 

JCGM, J.: Evaluation of measurement data—Guide to the expression of uncertainty in measurement, Int. Organ. Stand., Geneva, 50, 134, 2008. 

Jones, E. P. and Smith, S. D.: A first measurement of sea–air CO2 flux by eddy correlation, J. Geophys. Res., 82, 5990–5992, 1977. 

Kaimal, J. C. and Finnigan, J. J.: Atmospheric boundary layer flows: their structure and measurement, Oxford University Press, New York, 1994. 

Kaimal, J. C., Wyngaard, J. C., Izumi, Y., and Cote, O. R.: Spectral characteristics of surface-layer turbulence, Q. J. Roy. Meteor. Soc., 098, 563–589,, 1972. 

Kohsiek, W.: Water vapor cross-sensitivity of open path H2O/CO2 sensors, J. Atmos. Ocean. Tech., 17, 299–311,<0299:WVCSOO>2.0.CO;2, 2000. 

Kondo, F. and Tsukamoto, O.: Air–sea CO2 flux by eddy covariance technique in the equatorial Indian Ocean, J. Oceanogr., 63, 449–456,, 2007. 

Landschützer, P., Gruber, N., and Bakker, D. C. E.: An observation-based global monthly gridded sea surface pCO2 product from 1982 onward and its monthly climatology (NCEI Accession 0160558), Version 5.5, Dataset, NOAA National Centers for Environmental Information,, 2020. 

Landwehr, S., Miller, S. D., Smith, M. J., Saltzman, E. S., and Ward, B.: Analysis of the PKT correction for direct CO2 flux measurements over the ocean, Atmos. Chem. Phys., 14, 3361–3372,, 2014. 

Landwehr, S., Miller, S. D., Smith, M. J., Bell, T. G., Saltzman, E. S., and Ward, B.: Using eddy covariance to measure the dependence of air–sea CO2 exchange rate on friction velocity, Atmos. Chem. Phys., 18, 4297–4315,, 2018. 

Langford, B., Acton, W., Ammann, C., Valach, A., and Nemitz, E.: Eddy-covariance data with low signal-to-noise ratio: time-lag determination, uncertainties and limit of detection, Atmos. Meas. Tech., 8, 4197–4213,, 2015. 

Lauvset, S. K., McGillis, W. R., Bariteau, L., Fairall, C. W., Johannessen, T., Olsen, A., and Zappa, C. J.: Direct measurements of CO2 flux in the Greenland Sea, Geophys. Res. Lett., 38, L12603,, 2011. 

Lenschow, D. H. and Kristensen, L.: Uncorrelated noise in turbulence measurements, J. Atmos. Ocean. Tech., 2, 68–81, 1985. 

Lenschow, D. H., Mann, J., and Kristensen, L.: How long is long enough when measuring fluxes and other turbulence statistics?, NCAR Tech. Note, NCAR/TN-389, 53, Natl. Cent. for Atmos. Res., Boulder, Colo., 1993. 

Lenschow, D. H., Mann, J., and Kristensen, L.: How long is long enough when measuring fluxes and other turbulence statistics?, J. Atmos. Ocean. Tech., 11, 661–673,<0661:HLILEW>2.0.CO;2, 1994. 

Lenschow, D. H., Wulfmeyer, V., and Senff, C.: Measuring second- through fourth-order moments in noisy data, J. Atmos. Ocean. Tech., 17, 1330–1347,<1330:MSTFOM>2.0.CO;2, 2000. 

Loescher, H. W., Law, B. E., Mahrt, L., Hollinger, D. Y., Campbell, J., and Wofsy, S. C.: Uncertainties in, and interpretation of, carbon flux estimates using the eddy covariance technique, J. Geophys. Res.-Atmos., 111, 1–19,, 2006. 

Mahrt, L.: Flux sampling errors for aircraft and towers, J. Atmos. Ocean. Tech., 15, 416–429,<0416:FSEFAA>2.0.CO;2, 1998. 

Mauder, M., Cuntz, M., Drüe, C., Graf, A., Rebmann, C., Schmid, H. P., Schmidt, M., and Steinbrecher, R.: A strategy for quality and uncertainty assessment of long-term eddy-covariance measurements, Agr. Forest Meteorol., 169, 122–135,, 2013. 

McGillis, W. R., Edson, J. B., Ware, J. D., Dacey, J. W. H., Hare, J. E., Fairall, C. W., and Wanninkhof, R.: Carbon dioxide flux techniques performed during GasEx-98, Mar. Chem., 75, 267–280,, 2001. 

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., and Donelan, M.: Air–sea CO2 exchange in the equatorial Pacific, J. Geophys. Res.-Oceans, 109, C08S02,, 2004. 

Miller, S. D., Marandino, C., and Saltzman, E. S.: Ship-based measurement of air–sea CO2 exchange by eddy covariance, J. Geophys. Res.-Atmos., 115, 1–14,, 2010. 

Moat, B. and Yelland, M.: Airflow distortion at instrument sites on the RRS James Clark Ross during the WAGES project, No. 12, National Oceanography Centre Internal Document, National Oceanography Centre, Southampton, 2015. 

Moat, B. I., Yelland, M. J., and Cooper, E. B.: The airflow distortion at instruments sites on the RRS “James Cook”, National Oceanography Centre Southampton Research and Consultancy Report 11, National Oceanography Centre, Southampton, 44 pp, 2006. 

Nightingale, P. D., Malin, G., Law, C. S., Watson, A. J., Liss, P. S., Liddicoat, M. I., Boutin, J., and Upstill-Goddard, R. C.: In situ evaluation of air–sea gas exchange parameterizations using novel conservative and volatile tracers, Global Biogeochem. Cy., 14, 373–387,, 2000. 

Nilsson, E., Bergström, H., Rutgersson, A., Podgrajsek, E., Wallin, M. B., Bergström, G., Dellwik, E., Landwehr, S., and Ward, B.: Evaluating humidity and sea salt disturbances on CO2 flux measurements, J. Atmos. Ocean. Tech., 35, 859–875,, 2018. 

Post, H., Hendricks Franssen, H. J., Graf, A., Schmidt, M., and Vereecken, H.: Uncertainty analysis of eddy covariance CO2 flux measurements for different EC tower distances using an extended two-tower approach, Biogeosciences, 12, 1205–1221,, 2015. 

Prytherch, J., Yelland, M. J., Pascal, R. W., Moat, B. I., Skjelvan, I., and Neill, C. C.: Direct measurements of the CO2 flux over the ocean: Development of a novel method, Geophys. Res. Lett., 37, L03607,, 2010. 

Prytherch, J., Brooks, I. M., Crill, P. M., Thornton, B. F., Salisbury, D. J., Tjernström, M., Anderson, L. G., Geibel, M. C., and Humborg, C.: Direct determination of the air–sea CO2 gas transfer velocity in Arctic sea ice regions, Geophys. Res. Lett., 44, 3770–3778, 2017. 

Rannik, Ü., Mammarella, I., Aalto, P., Keronen, P., Vesala, T., and Kulmala, M.: Long-term aerosol particle flux observations part I: Uncertainties and time-average statistics, Atmos. Environ., 43, 3431–3439,, 2009. 

Rannik, Ü., Peltola, O., and Mammarella, I.: Random uncertainties of flux measurements by the eddy covariance technique, Atmos. Meas. Tech., 9, 5163–5181,, 2016. 

Smith, S. D. and Jones, E. P.: Evidence for wind-pumping of air–sea gas exchange based on direct measurements of CO2 fluxes, J. Geophys. Res.-Oceans, 90, 869–875, 1985. 

Spreen, G., Kaleschke, L., and Heygster, G.: Sea ice remote sensing using AMSR-E 89-GHz channels, J. Geophys. Res.-Oceans, 113, C02S03,, 2008. 

Takahashi, T., Sutherland, S. C., Wanninkhof, R., Sweeney, C., Feely, R. A., Chipman, D. W., Hales, B., Friederich, G., Chavez, F., Sabine, C., Watson, A., Bakker, D. C. E., Schuster, U., Yoshikawa-Inoue, H., Ishii, M., Midorikawa, T., Nojiri, Y., Körtzinger, A., Steinhoff, T., Hoppema, M., Olafsson, J., Arnarson, T. S., Johannessen, T., Olsen, A., Bellerby, R., Wong, C. S., Delille, B., Bates, N. R., and de Baar, H. J. W.: Climatological mean and decadal change in surface ocean pCO2, and net sea–air CO2 flux over the global oceans, Deep Sea Res. Pt. II, 56, 554–577,, 2009. 

Tsukamoto, O., Kondo, F., and Kamei, Y.: Overestimation of downward air–sea eddy CO2 flux due to optical window contamination of open-path gas analyzer, SOLA, 10, 117–121, 2014. 

Wanninkhof, R.: Relationship between wind speed and gas exchange over the ocean revisited, Limnol. Oceanogr.-Meth., 12, 351–362,, 2014. 

Wanninkhof, R., Asher, W. E., Ho, D. T., Sweeney, C., and McGillis, W. R.: Advances in Quantifying Air-Sea Gas Exchange and Environmental Forcing, Annu. Rev. Mar. Sci., 1, 213–244,, 2009. 

Weiss, R. F.: Carbon dioxide in water and seawater: the solubility of a non-ideal gas, Mar. Chem., 2, 203–215,, 1974. 

Weiss, A., Kuss, J., Peters, G., and Schneider, B.: Evaluating transfer velocity-wind speed relationship using a long-term series of direct eddy correlation CO2 flux measurements, J. Marine Syst., 66, 130–139,, 2007. 

Wesely, M. L., Cook, D. R., Hart, R. L., and Williams, R. M.: Air–sea exchange of CO2 and evidence for enhanced upward fluxes, J. Geophys. Res.-Oceans, 87, 8827–8832, 1982. 

Wienhold, F. G., Welling, M., and Harris, G. W.: Micrometeorological measurement and source region analysis of nitrous oxide fluxes from an agricultural soil, Atmos. Environ., 29, 2219–2227, 1995. 

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., Holding, T., Schuster, U., Girard-Ardhuin, F., Grouazel, A., Piolle, J. F., Warren, M., Wrobel-Niedzwiecka, I., Land, P. E., Torres, R., Prytherch, J., Moat, B., Hanafin, J., Ardhuin, F., and Paul, F.: Key uncertainties in the recent air–sea flux of CO2, Global Biogeochem. Cy., 33, 1548–1563,, 2019. 

Wyngaard, J. C.: Getting to know turbulence, Part I, Chapt. 2, in: Turbulence in the Atmosphere, Cambridge University Press, Cambridge, p. 27–54 2010.  

Yang, M., Nightingale, P. D., Beale, R., Liss, P. S., Blomquist, B., and Fairall, C.: Atmospheric deposition of methanol over the Atlantic Ocean, P. Natl. Acad. Sci. USA, 110, 20034–20039,, 2013. 

Yang, M., Prytherch, J., Kozlova, E., Yelland, M. J., Parenkat Mony, D., and Bell, T. G.: Comparison of two closed-path cavity-based spectrometers for measuring air–water CO2 and CH4 fluxes by eddy covariance, Atmos. Meas. Tech., 9, 5509–5522,, 2016. 

Yelland, M. J., Moat, B. I., Taylor, P. K., Pascal, R. W., Hutchings, J., and Cornell, V. C.: Wind stress measurements from the open ocean corrected for airflow distortion by the ship, J. Phys. Oceanogr., 28, 1511–1526,<1511:WSMFTO>2.0.CO;2, 1998. 

Short summary
Eddy covariance (EC) is the most direct method for measuring air–sea CO2 flux from ships. However, uncertainty in EC air–sea CO2 fluxes has not been well quantified. Here we show that with the state-of-the-art gas analysers, instrumental noise no longer contributes significantly to the CO2 flux uncertainty. Applying an appropriate averaging timescale (1–3 h) and suitable air–sea CO2 fugacity threshold (at least 20 µatm) to EC flux data enables an optimal analysis of the gas transfer velocity.
Final-revised paper