Articles | Volume 13, issue 2
https://doi.org/10.5194/amt-13-789-2020
https://doi.org/10.5194/amt-13-789-2020
Research article
 | 
19 Feb 2020
Research article |  | 19 Feb 2020

Ensemble-based satellite-derived carbon dioxide and methane column-averaged dry-air mole fraction data sets (2003–2018) for carbon and climate applications

Maximilian Reuter, Michael Buchwitz, Oliver Schneising, Stefan Noël, Heinrich Bovensmann, John P. Burrows, Hartmut Boesch, Antonio Di Noia, Jasdeep Anand, Robert J. Parker, Peter Somkuti, Lianghai Wu, Otto P. Hasekamp, Ilse Aben, Akihiko Kuze, Hiroshi Suto, Kei Shiomi, Yukio Yoshida, Isamu Morino, David Crisp, Christopher W. O'Dell, Justus Notholt, Christof Petri, Thorsten Warneke, Voltaire A. Velazco, Nicholas M. Deutscher, David W. T. Griffith, Rigel Kivi, David F. Pollard, Frank Hase, Ralf Sussmann, Yao V. Té, Kimberly Strong, Sébastien Roche, Mahesh K. Sha, Martine De Mazière, Dietrich G. Feist, Laura T. Iraci, Coleen M. Roehl, Christian Retscher, and Dinand Schepers
Abstract

Satellite retrievals of column-averaged dry-air mole fractions of carbon dioxide (CO2) and methane (CH4), denoted XCO2 and XCH4, respectively, have been used in recent years to obtain information on natural and anthropogenic sources and sinks and for other applications such as comparisons with climate models. Here we present new data sets based on merging several individual satellite data products in order to generate consistent long-term climate data records (CDRs) of these two Essential Climate Variables (ECVs). These ECV CDRs, which cover the time period 2003–2018, have been generated using an ensemble of data products from the satellite sensors SCIAMACHY/ENVISAT and TANSO-FTS/GOSAT and (for XCO2) for the first time also including data from the Orbiting Carbon Observatory 2 (OCO-2) satellite. Two types of products have been generated: (i) Level 2 (L2) products generated with the latest version of the ensemble median algorithm (EMMA) and (ii) Level 3 (L3) products obtained by gridding the corresponding L2 EMMA products to obtain a monthly 5×5 data product in Obs4MIPs (Observations for Model Intercomparisons Project) format. The L2 products consist of daily NetCDF (Network Common Data Form) files, which contain in addition to the main parameters, i.e., XCO2 or XCH4, corresponding uncertainty estimates for random and potential systematic uncertainties and the averaging kernel for each single (quality-filtered) satellite observation. We describe the algorithms used to generate these data products and present quality assessment results based on comparisons with Total Carbon Column Observing Network (TCCON) ground-based retrievals. We found that the XCO2 Level 2 data set at the TCCON validation sites can be characterized by the following figures of merit (the corresponding values for the Level 3 product are listed in brackets) – single-observation random error (1σ): 1.29 ppm (monthly: 1.18 ppm); global bias: 0.20 ppm (0.18 ppm); and spatiotemporal bias or relative accuracy (1σ): 0.66 ppm (0.70 ppm). The corresponding values for the XCH4 products are single-observation random error (1σ): 17.4 ppb (monthly: 8.7 ppb); global bias: −2.0 ppb (−2.9 ppb); and spatiotemporal bias (1σ): 5.0 ppb (4.9 ppb). It has also been found that the data products exhibit very good long-term stability as no significant long-term bias trend has been identified. The new data sets have also been used to derive annual XCO2 and XCH4 growth rates, which are in reasonable to good agreement with growth rates from the National Oceanic and Atmospheric Administration (NOAA) based on marine surface observations. The presented ECV data sets are available (from early 2020 onwards) via the Climate Data Store (CDS, https://cds.climate.copernicus.eu/, last access: 10 January 2020) of the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, last access: 10 January 2020).

1 Introduction

Carbon dioxide (CO2) and methane (CH4) are important greenhouse gases and increasing atmospheric concentrations result in global warming with adverse consequences such as sea level rise (IPCC, 2013). Because of their importance for climate, these gases have been classified as Essential Climate Variables (ECVs) by the Global Climate Observing System (GCOS) (GCOS-154, 2010; GCOS-200, 2016). The generation of XCO2 and XCH4 satellite-derived ECV data products meeting GCOS requirements using European satellite retrieval algorithms started in 2010 in the framework of the GHG-CCI project (http://www.esa-ghg-cci.org/, last access: 10 January 2020) of the European Space Agency's (ESA) Climate Change Initiative (CCI) (Hollmann et al., 2013). Since the end of 2016, this activity continues operationally via the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, last access: 10 January 2020), and the corresponding CO2 and CH4 data products are available via the Copernicus Climate Data Store (CDS, https://cds.climate.copernicus.eu/, last access: 10 January 2020). These ECV data products have been used for a range of applications such as improving our knowledge of CO2 and/or CH4 surface fluxes (e.g., Alexe et al., 2015; Basu et al., 2013; Buchwitz et al., 2017a; Chevallier et al., 2014, 2015; Ganesan et al., 2017; Gaubert et al., 2019; Houweling et al., 2015; Liu et al., 2017; Maasakkers et al., 2019; Miller et al., 2019; Reuter et al., 2014a, b, 2019a; Sheng et al., 2018; Schneising et al., 2014b; Turner et al., 2015, 2019), comparison with climate and other models (e.g., Hayman et al., 2014; Lauer et al., 2017; Schneising et al., 2014a), and for other applications such as computation of CO2 growth rates (e.g., Buchwitz et al., 2018), as well as to better understand changes in the amplitude of the CO2 seasonal cycle (e.g., Yin et al., 2018).

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f01

Figure 1Overview of the presented XCO2 data set. Shown are time series over land for three latitude bands (global, black line; Northern Hemisphere, red; Southern Hemisphere, green) and global maps (half-yearly averages at 1×1 obtained by gridding (averaging) the merged Level 2, i.e., EMMA, product). See Sect. 4 for a detailed discussion.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f02

Figure 2As Fig. 1 but for XCH4.

The C3S satellite greenhouse gas (GHG) data set consists of single-sensor satellite data products and of merged (i.e., combined multi-sensor, multi-algorithm) data products. Here we present the latest version, version 4.1, of the merged Level 2 (L2) and merged Level 3 (L3) XCO2 and XCH4 data products, which cover the time period 2003–2018. The L2 products (XCO2_EMMA and XCH4_EMMA) have been compiled with the ensemble median algorithm (EMMA) originally proposed by Reuter et al. (2013) and recent modifications, which are described in Sect. 3.1. These products contain detailed information for each single observation (i.e., footprint or ground pixel) including time, latitude and longitude, the main parameter (i.e., XCO2 or XCH4), its stochastic uncertainty (e.g., due to instrument noise), an estimate of potential systematic uncertainties (e.g., due to spatial or temporal bias patterns), and its averaging kernel and corresponding a priori profile. The L3 products (XCO2_OBS4MIPS and XCH4_OBS4MIPS) are gridded products at monthly time and 5×5 spatial resolution in Obs4MIPs (Observations for Model Intercomparisons Project, https://www.earthsystemcog.org/projects/obs4mips/, last access: 10 January 2020) format.

Figure 1 provides an overview of the resulting merged XCO2 data product in terms of time series for three latitude bands and global maps and the similarly structured Fig. 2 shows the XCH4 product. As can be seen, XCO2 and XCH4 are both increasing with time and exhibit seasonal fluctuations and spatial variations. The spatiotemporal characteristics of the merged data – e.g., the spatial sampling – reflect the characteristics of the underlying individual sensor satellite data (described in the data section, Sect. 2). Figures 1 and 2 are discussed in detail in the results section, Sect. 4. How these data products have been generated is described in the methods section, Sect. 3. A summary and conclusions are given in Sect. 5.

Table 1Satellite XCO2 Level 2 (L2) data products used as input for the generation of the merged L2 and L3 XCO2 version 4.1 data products. For products which have been generated in the framework of the CCI and C3S projects the corresponding product ID is listed (the other products are external products, which have been obtained from the corresponding websites; see Acknowledgements). Temporal coverage indicates the time coverage of the input data sets.

Download Print Version | Download XLSX

2 Data

In this section, we present an overview about the input data used to generate and validate the new XCO2 and XCH4 data products.

2.1 Satellite data

The input satellite data used to generate the merged satellite data products are individual satellite sensor Level 2 (L2) data products. Table 1 provides an overview about the satellite XCO2 input data sets. As can be seen, in total eight XCO2 L2 data products have been used to generate the merged L2 and Level 3 (L3) XCO2 data products, each corresponding to a different combination of satellite sensor and retrieval algorithm. An overview about the time coverage of these input data products is presented in Fig. 3. As can be seen, the time period 2003 to March 2009 is only covered by one XCO2 product, namely XCO2 retrieved with the Bremen Optimal Estimation DOAS (BESD) algorithm (Reuter et al., 2010, 2011) from the SCIAMACHY/ENVISAT (Bovensmann et al., 1999) instrument. A second SCIAMACHY XCO2 data product is available, which has been retrieved with the Weighting Function Modified Differential Optical Absorption Spectroscopy (WFM-DOAS or WFMD) algorithm (Schneising et al., 2011), but this product is not used because the merging algorithm EMMA (Reuter et al., 2013, described in Sect. 3.1) requires one or more than two input data products (because the median of a set of elements is, according to our definition which avoids averaging, not defined for two elements). Therefore, one of the two products had to be selected, and the choice was the BESD product for XCO2 because of somewhat higher data quality compared to the WFMD product (Buchwitz et al., 2017b) (note however that the WFMD product has the advantage of containing a larger number of observations). As can be seen from Table 1 and Fig. 3, several GOSAT input products have been used from April 2009 onwards and two OCO-2 XCO2 products from September 2014 and May 2015 onwards. Note that additional algorithms/data products are available but have not been used as input, for example the GOSAT BESD XCO2 product (Heymann et al., 2015) and the OCO-2 RemoTeC XCO2 product (Wu et al., 2018). These or other additional products may be added in future versions of the merged XCO2 products. Note also that we always use the bias-corrected version of a data product, if available (some product files contain bias-corrected and uncorrected values).

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f03

Figure 3Individual satellite sensor XCO2 data products contributing to the merged XCO2 data products (see Table 1 for details). The required minimum number of contributing products is shown by the grey area.

Download

All listed satellites perform nadir (down-looking) and glint observations and provide radiance spectra covering the relevant CO2 and CH4 absorption bands located in the shortwave infrared (SWIR) part of the electromagnetic spectrum (around 1.6 and 2 µm) and also cover the O2 A-band spectral region in the near-infrared (NIR, around 0.76 µm). All individual sensor input L2 data products have been generated using retrieval algorithms based on minimizing the difference between a modeled radiance spectrum and the observed spectrum by modifying so-called state vector elements (for details we refer to the references listed in Table 1; for additional information see also the Algorithm Theoretical Basis Documents (ATBDs); Buchwitz et al., 2019b, and Reuter et al., 2019b). The exact definition of the state vector depends on the algorithm, but the general approach is based on the optimal estimation (Rodgers, 2000) formalism or similar approaches (see references in Table 1). Among the state vector elements is a representation of the CO2 vertical profile but also other parameters to consider interfering gases (e.g., water vapor), surface reflection, atmospheric scattering, and other effects and parameters, which have an impact on the (interpretation of the) measured radiance spectrum.

Table 2As Table 1 but for XCH4.

Download Print Version | Download XLSX

Table 3TCCON sites used for the validation of the XCO2 and XCH4 satellite-derived data products.

Download Print Version | Download XLSX

Table 2 and Fig. 4 provide an overview about the satellite XCH4 L2 input data sets. As for XCO2, the time period 2003 to March 2009 is only covered by one SCIAMACHY data product. From April 2009 onwards several GOSAT XCH4 products are available (see Table 2) and have been used to generate the merged XCH4 data L2 and L3 data products. For future updates it is also planned to include XCH4 from the Sentinel-5 Precursor (S5P) satellite (Veefkind et al., 2012), but S5P XCH4 (Hu et al., 2018; Schneising et al., 2019) has not yet been included as the time period covered by these products is currently quite short (less than 2 years). However, we aim to include S5P XCH4 for one of the next updates of the merged methane products.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f04

Figure 4As Fig. 3 but XCH4. For details on each product see Table 2.

Download

2.2 Ground-based data

The satellite data products have been validated by comparison with the XCO2 and XCH4 data products of the TCCON (Wunch et al., 2011). TCCON is a network of ground-based Fourier transform spectrometers (FTSs) recording direct solar spectra in the NIR/SWIR spectral region. From these spectra, accurate and precise column-averaged abundances of CO2, CH4, and a number of other species are retrieved. The TCCON data products (version GGG2014) have been obtained via the TCCON data archive (https://tccondata.org/, last access: 15 July 2019). An overview about the used TCCON sites is presented in Table 3.

In Sect. 4.3, we present annual XCO2 and XCH4 growth rates, which have been derived from the new XCO2 and XCH4 OBS4MIPS data products using the method described in Buchwitz et al. (2018). These growth rates are compared with growth rates derived from marine surface CO2 and CH4 observations, which have been obtained from the National Oceanic and Atmospheric Administration (NOAA) (for details including links and last access see Acknowledgements).

3 Methods

3.1 Merging algorithm EMMA

In order to generate the merged L2 products, the ensemble median algorithm is used, which is described in detail in Reuter et al. (2013). Therefore, we limit the description given here to a short overview of the latest version of the EMMA algorithm. To be specific, we initially describe the EMMA XCO2 algorithm and explain differences relevant for XCH4 at the end of this subsection.

The EMMA XCO2 data product consists of selected individual L2 soundings from the available individual sensor L2 input products (listed in Table 1). The EMMA L2 product is based on selecting the best soundings (i.e., single ground pixel observations) from the ensemble of individual sensor L2 products. Sounding selection is based on monthly time and 10×10 spatial intervals. To decide which individual product is selected for a given month and given grid cell, all input products are first gridded (monthly, 10×10) to consider the fact that the spatiotemporal sampling is different for each individual product (due to different satellite sensors and algorithm-dependent quality-filtering strategies). The selected product is the median in terms of average XCO2 per month and grid cell (note that in case of an even number of products the product which is closest to the mean is selected). The median is used primarily to remove potential outliers. The advantage of the median is also (in contrast to, for example, the arithmetic mean) that no averaging or other modifications to the input data are required. In order for a grid cell to be assigned a valid value, the following criterion has to be fulfilled: a minimum number of data products having a standard error of the mean (SEOM) of less than 1 ppm has to be available (see grey area in Fig. 3). SEOM is defined by 1niσi2, with σi being the (scaled; see below) XCO2 uncertainty of the ith out of n soundings.

This means that EMMA selects for each month and each 10×10 grid cell exactly one product of the available individual L2 input products and then transfers all relevant information (i.e., XCO2 and its uncertainty, related averaging kernels and a priori profile, etc.) from the selected original L2 file into the corresponding daily EMMA L2 product file. This ensures that most of the original information from the selected individual product is also contained in the merged product.

However, some modifications are applied. In order to remove (or at least to minimize) the impact of different a priori assumptions, all products are converted to common a priori CO2 vertical profiles (see Reuter et al., 2013, for details). The new a priori profiles are obtained from the simple empirical CO2 model (SECM, Reuter et al., 2012). SECM is essentially an empirically found function with parameters optimized using a CO2 model (CT2017; see below). The SECM model used here is referred to as SECM2018 and is an update of the SECM model described in Reuter et al. (2012). The main difference is that SECM2018 is using a recent version of NOAA's assimilation system CarbonTracker (Peters et al., 2007, with updates documented at: http://carbontracker.noaa.gov/, last access: 10 January 2020), namely CT2017.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f05

Figure 5Global bias correction as applied by EMMA to the individual satellite XCO2 input data products. Panel (a) shows the difference relative to the SECM2018 model (computed as satellite – model) before the correction and (b) shows the difference after correction.

Download

SECM2018 is also used to correct for potential offsets between the individual data products by adding or subtracting a global offset (i.e., by using one constant offset value for each individual product applied globally and for the full time series). Time series of the individual data products before and after offset correction are shown in Fig. 5. Note that in Fig. 5 all data are relative to SECM2018, which is a very simple CO2 model, and therefore all variations and trends seen in Fig. 5 are at least to some extent model errors. As can be seen from Fig. 5, the correction brings the individual data sets typically closer together without changing any of their other characteristics (e.g., their time dependence). But as can also be seen from Fig. 5, better agreement is only achieved on average, not necessarily for all products during the entire time period. For example, the GOSAT RemoTeC product (blue curve) during 2009–2012 exhibits a somewhat larger difference after the offset correction. The approximately 2 ppm (0.5 %) spike at the beginning of the time series is likely due to a positive bias of the underlying BESD data product, which has not been corrected due to the lack of reference data in this time period (see also the discussion of this aspect in Buchwitz et al., 2018). An obvious issue is also the approximately 1.5 ppm (0.4 %) discontinuity in the first half of 2014 of the PPDF-S (photon path length probability density function/simultaneous) product (light-green curve). Depending on application, this may be an issue when this product is used stand-alone, but this is not a problem for EMMA as EMMA identifies and ignores outliers.

Another modification applied to the individual L2 input products is a potential scaling of their reported uncertainty for the individual L2 soundings. The scaling factor has been chosen such that on average the uncertainty of the reported error is consistent with the standard deviation of satellite minus ground-based validation data differences (see Sect. 4.1 for the validation of the reported uncertainties via the uncertainty ratio).

In order to avoid that an individual input product, which has much more observations than the other products (such as OCO-2 compared to GOSAT), entirely dominates the EMMA product, a method has been implemented to prevent overweighting the contributions from individual L2 input data products. The method is based on limiting the number of L2 data points. For each grid cell and month, we perform the following steps: first, we compute SEOM for each algorithm. From these values, we compute the 25th percentile and divide it by 2. The result is used as the minimum SEOM threshold. If SEOM of an individual algorithm is smaller than this threshold, a subset of soundings is randomly chosen such that SEOM becomes just larger than the threshold. If, for example, all σi are 1 ppm, then SEOM simply becomes 1/n. If in this case, for example, data from four algorithms were available with n1=60, n2=80, n3=100, and n4=1000, the SEOM threshold would become 1/2n3, which would effectively limit the number of soundings of the fourth algorithm to 200 (chosen randomly).

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f06

Figure 6Relative data weight (a) and soundings per month (a) of the individual satellite XCO2 data products contributing to the EMMA XCO2 data product.

Download

In addition to the L2 information of the selected data products, EMMA stores the following diagnostic information for each selected sounding: identifier for the selected L2 algorithm and inter-algorithm spread (IAS) within the grid box of the sounding. Within each grid box, IAS is defined as the algorithm-to-algorithm standard deviation of the grid box averages.

By how much each individual satellite XCO2 data product contributes to the EMMA XCO2 product is shown in Fig. 6. Figure 6a shows the relative data weight (RDW), and Fig. 6b shows the number of soundings per month. How the RDW is defined is explained in detail in Reuter et al. (2013). In short, the RDW is defined as the relative number of soundings weighted with the corresponding (square of the inverse) uncertainty. RDW is high if a (relatively) large number of soundings contribute to the EMMA product and if these soundings have (relatively) low uncertainty compared to the other contributing products. The RDW of a product is a measure of how much information on XCO2 this product contributes to the EMMA product relative to the other contributing products. As can be seen from Fig. 6, the SCIAMACHY BESD product is the only product until early 2009, when the GOSAT time series starts. As can also be seen, OCO-2 dominates in terms of RDW and number of soundings from 2015 onwards. This is because OCO-2 provides much more data with typically better uncertainty compared to the other (GOSAT) product.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f07

Figure 7As Fig. 5 but for XCH4 and using methane model SC4C2018.

Download

The EMMA L2 XCH4 product has been generated similarly to the EMMA L2 XCO2 product, i.e., using essentially the same method as described above. A difference is that the offset correction has been done with a CH4 model instead of SECM2018. This model is the simple CH4 climatological model (SC4C), and we use the year 2018 update referred to as SC4C2018 in the following. The SC4C2018 model is similar to SECM2018 but for XCH4. It is a model-based CH4 climatology adjusted for the annual growth rate (note that this model has also been used as the climatological training and a calibration data set as described in Schneising et al., 2019). The EMMA algorithm SEOM limit controlling the minimum number of data points per grid box, month, and algorithm has been set to 12 ppb for XCH4. The impact of the offset correction for merging the XCH4 products is shown in Fig. 7. Note that in Fig. 7 all data are relative to SC4C2018, which is a very simple CH4 model, and therefore all variations and trends seen in Fig. 7 are at least to some extent model errors. As for CO2 (Fig. 5) the offset correction typically brings the various XCH4 products closer together but does not change any of their other characteristics. The PPDF-S product suffers from a discontinuity (of 8 ppb or 0.4 %) in the first half of 2014 (see above for a similar problem for PPDF-S XCO2).

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f08

Figure 8As Fig. 6 but for XCH4.

Download

Figure 8a shows the RDW, and Fig. 8b shows the number of soundings per month for all individual sensor XCH4 products contributing to the XCH4 EMMA product. Until early 2009, the SCIAMACHY WFMD product is the only product contributing to the EMMA product. Note that the RDW of the SCIAMACHY products drops at the end of 2005 (in contrast to the absolute number of soundings per month). The reason is the increase in the uncertainty of this product due to detector degradation (see, e.g., Schneising et al., 2011, for details). As can also be seen, the two GOSAT proxy (PR) products (i.e., CH4_GOS_OCPR and CH4_GOS_SRPR) dominate the XCH4 EMMA product because they contain more soundings compared to the other (GOSAT) data products.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f09

Figure 9April 2011 XCO2 at 10×10 spatial resolution showing (i) the individual sensor/algorithm input data sets (panels in rows 1–4; see Table 1 for details), (ii) EMMA XCO2 (bottom left), and (iii) the inter-algorithm spread (IAS, 1σ) as computed by EMMA (bottom right; see main text for details). Also shown in the bottom-right panel are the locations of the TCCON sites (pink triangles) and the range of IAS values covered by them (see color bar). Note that the OCO-2 maps (row 4) are empty because this satellite was launched after April 2011 (see Fig. 10 for OCO-2 XCO2).

3.2 Algorithm to generate the Level 3 OBS4MIPS products

The version 4.1 L3 XCO2_OBS4MIPS and XCH4_OBS4MIPS data products have been obtained by gridding (averaging) the version 4.1 L2, i.e., XCO2_EMMA and XCH4_EMMA, products using monthly time and 5×5 spatial resolution. The algorithm for the generation of the OBS4MIPS products is described in Reuter et al. (2019b). Therefore, we here provide only a short overview.

For each individual product, the gridding is based on computing an arithmetic, unweighted average of all soundings falling in a grid box.For each grid box, the standard error of the mean is computed using the uncertainties contained in the corresponding EMMA product files. In order to reduce noise at least two individual observations must be present and the resulting standard error of the mean must be less than 1.6 ppm for XCO2 and less than 12 ppb for XCH4.

Besides XCO2 or XCH4, the final L3 product also includes (per grid box and month) the number of soundings used for averaging; the average column-averaging kernel; the average a priori profile; the standard deviation of the averaged XCO2 or XCH4 values; and an estimate for the total uncertainty computed as the root sum square of two values, where one value is SEOM and the other value is IAS as computed by EMMA. For cases including only one algorithm, the second value is replaced by quadratically adding spatial and seasonal accuracy determined from the TCCON validation.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f10

Figure 10As Fig. 9 but for April 2015. Note that the SCIAMACHY/BESD map (top left) is empty because this product ended in April 2012 (see Fig. 9 for SCIAMACHY/BESD XCO2).

3.3 Validation method

The validation of the merged satellite-derived XCO2 and XCH4 data products is based on comparisons with ground-based XCO2 and XCH4 TCCON observations (using version GGG2014). We present results from two somewhat different validation methods (the EMMA method, Reuter et al., 2013; and the QA/QC method, Buchwitz et al., 2017b; see below), which are similar to other validation methods used in recent years (e.g., Butz et al., 2010; Cogan et al., 2012; Dils et al., 2014; O'Dell et al., 2018; Parker et al., 2011). These methods differ with respect to details such as the chosen colocation criterion, whether the data are brought to a common a priori or not, and if yes which a priori has been used. In the following, we will highlight some of these details as relevant for the two validation methods used for this paper.

Both methods used for the validation of the L2 EMMA products are based on colocating each individual satellite XCO2 (or XCH4) observation with a corresponding value obtained from TCCON using predefined spatial and temporal colocation criteria (see below). The comparisons take into account different a priori assumptions regarding the vertical profiles of CO2 (or CH4) as used for the generation of the L2 input products by converting either the satellite data (QA/QC method) or the TCCON data (EMMA method) to a common a priori. This a priori correction is based on using the satellite averaging kernels and a priori profiles, which are contained (for each single observation) in the EMMA product files. The magnitude of the a priori correction (the explicit formula is shown as Eq. 3 in Dils et al., 2014) depends on the deviation (difference) of the averaging kernel from unity and on the difference of the a priori profiles. Because the averaging kernel profiles are typically close to unity (note that both satellite and the TCCON retrievals correspond to cloud-free conditions) and because the a priori profiles are not totally unrealistic, the a priori correction is typically very small (approximately 0.1 ppm for XCO2 and 1 ppb for XCH4).

Table 4Overview validation results at TCCON sites for data product XCO2_EMMA (version 4.1).

Download Print Version | Download XLSX

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f11

Figure 11Average XCO2 inter-algorithm spread (1σ) during 2003–2018. As can be seen, the scatter is typically around 1 ppm except over parts of the tropics (in particular central Africa), the Himalayas, and at high latitudes, where the scatter can be larger.

The first validation method is the EMMA quality assessment method, which is described in Reuter et al. (2013). Note that EMMA is not only a merging method but also a data quality assessment method, as the assessment of the quality of all satellite input data (listed in Tables 1 and 2) is a key aspect of EMMA. The second method is the quality assessment/quality control (QA/QC) method (Buchwitz et al., 2017b), which is applied to all satellite XCO2 and XCH4 data products generated for the Copernicus Climate Change Service (C3S), i.e., to the merged products but also to all the individual sensor CCI/C3S L2 input products, which are also available via the Copernicus Climate Data Store (CDS) (see products with CCI/C3S product ID listed in Tables 1 and 2).

Key differences between the QA/QC method and the EMMA method are listed as follows.

  • Colocation criteria: QA/QC used ±2 latitude and ±4 longitude as the spatial colocation criterion, but EMMA used 500 km (both methods use the same temporal colocation criterion of 2 h).

  • Filtering criterion surface elevation: EMMA requires a surface elevation difference of less than 250 m between a TCCON site and satellite footprints, whereas the QA/QC does not use this filtering criterion.

  • A priori correction: both methods correct for the use of different a priori CO2 vertical profiles in the various retrieval algorithms, but QA/QC uses the TCCON a priori as common a priori, whereas EMMA uses the SECM2018 model for CO2 and the SC4C2018 model for CH4 (see Sect. 3.1).

  • Approach to quantify seasonal bias and linear bias trend: the EMMA method is based on fitting a trend model, which includes an offset term, a slope term, and a sine term for seasonal fluctuations (see Reuter et al., 2019c) and computes the seasonal bias from the standard deviation of the fitted seasonal fluctuation term and obtains the bias trend and its uncertainty from the fitted slope term. The QA/QC method (Buchwitz et al., 2019a) uses (only) a linear fit to obtain the bias trend and its uncertainty and computes the seasonal bias from the standard deviation of the seasonal biases (as also done by Dils et al., 2014, for their quantity seasonality).

  • Criteria for enough data: both algorithms use several different thresholds for the required minimum number of colocations per TCCON site and minimum length of overlapping TCCON time series.

Despite all these differences, quite similar overall figures of merit have been obtained with both methods (see results section, Sect. 4). This indicates that the overall data quality results do not critically depend on the details of the assessment method (the same conclusion has also been reported for earlier comparisons of results from different assessment methods, e.g., Buchwitz et al., 2015, 2017b).

4 Results and discussion

4.1 Products XCO2_EMMA and XCO2_OBS4MIPS (v4.1)

When generating an EMMA product, a set of standard figures are generated such as Figs. 5 and 6 already discussed but also maps of the EMMA product and of the various input data products for all months of the 2003–2018 time period. Two of these figures are shown here, namely the figures for April 2011 (Fig. 9) and April 2015 (Fig. 10) (note that 2011 is the last full year with data from SCIAMACHY and that 2015 is the first full year with OCO-2 data). The maps in the first four rows of Figs. 9 and 10 show the individual sensor/algorithm L2 input data. As can be seen, the spatial XCO2 patterns are quite similar (e.g., north–south gradient), but there are also significant differences, especially with respect to the spatial coverage. The spatial coverage depends on time and is related to the different satellite instruments but also due to algorithm-dependent quality filtering. The largest differences are between the SCIAMACHY BESD product (top left in Fig. 9) compared to the other products, as the SCIAMACHY product is limited to observations over land, whereas the GOSAT and OCO-2 products also have some ocean coverage due to the ocean-glint mode, which permits the acquisition of an adequate signal (and therefore also signal-to-noise ratio) also over the ocean (note that the reflectivity of water is poor outside of sunglint conditions in the used SWIR spectral regions around 1.6 and 2 µm). The EMMA product is shown in the bottom-left panels of Figs. 9 and 10, and in the bottom-right panel IAS is shown, which quantifies the level of agreement (or disagreement) among the various satellite input data sets. The IAS maps also show the location of the TCCON sites (pink triangles) and the IAS values at the TCCON sites (see pink triangles above the color bar). As can be seen, the TCCON sites are typically located outside of regions where the IAS is highest.

The average IAS for the entire time period 2003–2018 is shown in Fig. 11. As can be seen, the scatter is typically in the range 0.6–1.1 ppm with the exception of parts of the tropics, in particular central Africa, the Himalayas, parts of southeast Asia, and high latitudes. High latitudes typically correspond to large solar zenith angles, which is a challenge for accurate satellite XCO2 retrievals, as this typically corresponds to low signal and therefore low signal-to-noise ratio resulting in enhanced scatter of the retrieved XCO2. In areas with frequent cloud coverage, such as parts of the tropics, sampling is sparse and this may also contribute to a larger scatter.

Detailed validation results for all individual sensors and the EMMA XCO2 Level 2 data products are shown in Appendix A (Fig. A1) for all TCCON sites. The validation results are summarized in Table 4 (per site) and Table 5 (overall) together with the corresponding results of the QA/QC assessment method.

Table 4 lists all TCCON sites, which fulfill either the EMMA method or the QA/QC method criteria with respect to a minimum number of colocations and length of time series. Listed are the numerical values (in ppm), which have been computed for several figures of merit. This includes (i) the overall estimation of the single-observation random error computed as the standard deviation of the satellite minus TCCON differences; (ii) the uncertainty ratio, which is the ratio of the mean value of the reported (1σ) uncertainty to the standard deviation of the satellite–TCCON difference (computed to validate the reported uncertainties); (iii) the overall bias computed as the mean value of the satellite–TCCON differences; and (iv) the seasonal bias, computed as the standard deviation of the biases determined for the four seasons. Also shown in the last two rows are the mean value and the standard deviation of the values listed per TCCON site in the rows above. Several of these values have been used to compute the values listed in Table 5, which shows the overall summary of the quality assessment.

Table 5Validation summary for data product XCO2_EMMA (version 4.1).

Download Print Version | Download XLSX

Table 5 lists (i) the mean value of the single-observation random error, (ii) the global bias computed as the mean value of the biases at the various TCCON sites, (iii) the regional bias computed as the standard deviation of the biases at the various TCCON sites, (iv) the mean seasonal bias, and (v) the spatiotemporal bias computed as the root sum square of the regional and of the seasonal bias. The spatiotemporal bias is used to quantify the achieved performance for relative accuracy, which characterizes the spatially and temporally varying component of the bias (i.e., neglects a possible global bias (global offset), which is reported separately).

The linear bias trend has also been computed by fitting a line to the satellite–TCCON differences (not shown here). The mean value of the linear trend (slope) and its uncertainty (1σ, obtained from the standard deviation of the slope at the various TCCON sites) are -0.05±0.06 ppm yr−1 for the EMMA method and -0.06±0.09 ppm yr−1 for the QA/QC method. This means that no significant long-term bias trend has been detected; i.e., the satellite product is stable.

As can be seen from Table 5, the values computed independently using the EMMA and the QA/QC assessment methods are quite similar, which gives not only confidence in the overall quality assessment summary documented in Table 5 but also in the products and the used validation methods.

Table 6TCCON XCO2 bias in parts per million (ppm; satellite – TCCON). “–” means that the number of available colocations is less than the threshold required by the corresponding assessment method. Note that this table includes only a subset of the 10 sites shown in Fig. 12, namely only those sites with a mean bias being considerably (more than 1.5 times) larger than the standard deviation of the biases.

Assessment method DP is the method used by the data provider. For a see Boesch et al. (2019). For b see Wu et al. (2019).

Download Print Version | Download XLSX

Note however that the quality of the satellite data (at least at TCCON sites) is very likely better than Table 5 suggests (i) because the TCCON retrievals are not free of errors (the 1σ XCO2 uncertainty is about 0.4 ppm; Wunch et al., 2010) and (ii) because of the representation error originating from the (real) spatiotemporal variability of XCO2 around the TCCON sites. The overall error related to this is difficult to quantify, but some indication can potentially be obtained by additional assessment results such as the one shown in Fig. 12. Figure 12 shows the biases as obtained with the EMMA method at the various TCCON sites used for the EMMA method comparisons. Shown are not only the mean satellite–TCCON differences as obtained for the EMMA product but also for all the individual sensor/algorithm input products. The differences are shown as anomalies with respect to the mean; i.e., the sum of the differences in each row is zero. This is equivalent to assuming that for a given satellite product the mean value over all TCCON sites is zero. As can be seen from Fig. 12, the satellite–TCCON differences are dominantly positive (orange and red colors) for higher-latitude TCCON sites and mostly negative (blue colors) for lower-latitude TCCON sites. In order to rule out that this is an artifact of the EMMA assessment method, the overall biases computed with the QA/QC method and biases computed by the individual product data providers (DPs) have also been derived. These biases have been used to compute – for each of the 10 TCCON sites shown in Fig. 12 – the mean bias and the standard deviation of these biases. For 4 of these 10 sites the mean bias is considerably (more than 1.5 times) larger than the standard deviation of the biases, and the corresponding results for these four sites are shown in Table 6. This does not necessarily mean that these sites have the largest biases. This only means that the derived biases at these sites are (independent of their magnitude) the most consistent across all satellite products used for comparison. As can be seen from Table 6, the biases are always positive at Sodankylä, Karlsruhe, and Orléans and always negative at Lamont. Note that this does not imply that all derived biases are significant as some biases are very small, e.g., the FOCAL bias at Sodankylä, which is only 0.02 ppm. Because it is unlikely that all three satellites and several retrieval algorithms produce XCO2 products with similar biases at a given TCCON site, this provides an indication of biases either due to representation errors or due to biases within the TCCON data (Table 6). Note that these biases are within the accuracy stated by TCCON, which is 0.8 ppm (2σ) (Wunch et al., 2010, Hedelius et al., 2017). The accuracy of the TCCON data will be improved for the next data release (planned for 2020). This new TCCON data set will allow for better identification of the causes for the observed biases.

Table 7Overview validation results at TCCON sites for data product XCH4_EMMA (version 4.1).

Download Print Version | Download XLSX

The XCO2_OBS4MIPS product has also been directly compared with TCCON using a comparison method based on the comparison of the monthly satellite product with TCCON monthly mean values. The results are shown in Fig. 13. As can be seen, the mean difference (satellite – TCCON) is 0.18 ppm (which is close to the mean value of the global bias of 0.20 ppm listed in Table 5), the standard deviation is 1.18 ppm (as expected, because of the spatiotemporal averaging, which is somewhat smaller than the value of 1.29 ppm obtained for the XCO2_EMMA product listed in Table 5), and the linear correlation coefficient is 0.99. The spatiotemporal bias, computed as the standard deviation of 3-monthly averages at the TCCON sites listed in Fig. 13, is 0.7 ppm.

Figure 1 presents an overview of the XCO2 data product in terms of time series for three latitude bands and global maps. XCO2 is increasing almost linearly during the 16-year time period (for a discussion of the derived annual growth rates see Sect. 4.3). The main reason for this increase is CO2 emission due to burning of fossil fuels (Le Quéré et al., 2018). The seasonal cycle, which is caused primarily by quasi-regular uptake and release of atmospheric CO2 by the terrestrial vegetation due to photosynthesis and respiration (e.g., Kaminski et al., 2017, Yin et al., 2018), is most pronounced over the Northern Hemisphere. The half-yearly maps for 2003 are based on SCIAMACHY on board ENVISAT (Burrows et al., 1995; Bovensmann et al., 1999) satellite data, and the maps for 2018 contain data from the GOSAT (since 2009) (Kuze et al., 2016) and OCO-2 (since 2014) (Crisp et al., 2004) satellites. GOSAT and OCO-2 also provide good-quality XCO2 retrievals over the oceans due to their sunglint observation mode.

4.2 Products XCH4_EMMA and XCH4_OBS4MIPS (v4.1)

As for XCO2, monthly maps have also been generated for the EMMA XCH4 data product. Two examples are shown in Fig. 14 for September 2010 and in Fig. 15 for September 2018. The individual sensor XCH4 input data are shown in the first four rows, and the EMMA XCH4 product is shown in the bottom-left panel. The bottom-right panel shows the IAS. As can be seen, the spatial patterns of the XCH4 maps are similar but not identical. The IAS shows a quite large variability. The scatter is larger compared to the corresponding XCO2 IAS (Figs. 9 and 10, bottom-right panels), and spatially the grid cells with larger spread are more equally distributed over the globe but with largest differences over the southern part of Asia.

Detailed validation results are shown in Appendix A (Fig. A2), and the validation results are summarized in Tables 7 and 8, which have the same structure as the corresponding XCO2 tables (Tables 4 and 5). These tables also list the results of the QA/QC assessment method, which results in quite similar (within a few ppb) overall quality assessment results (Table 8) as obtained with the EMMA method. The linear bias trend has also been computed by fitting a line to the satellite–TCCON differences (not shown here). The mean value of the linear trend (slope) and its uncertainty (1σ, obtained from the standard deviation of the slope at the various TCCON sites) are -0.1±0.4 ppb yr−1 for the EMMA method and 0.5±0.8 ppb yr−1 for the QA/QC method. As for XCO2, this means that no significant long-term bias trend has been detected; i.e., the satellite product is stable.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f12

Figure 12Average XCO2 differences (satellite – TCCON) for the different satellite XCO2 products at 10 TCCON sites as used by the EMMA assessment method. The differences are shown as anomalies; i.e., the sum of the values corresponding to a given row is zero. Note that here “ACOS” refers to NASA's ACOS algorithm as applied to GOSAT and that “NASA” refers to NASA's ACOS algorithm as applied to OCO-2.

Download

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f13

Figure 13Summary of the comparison of product XCO2_OBS4MIPS with TCCON monthly mean XCO2 (each symbol corresponds to one month and to one TCCON site; each color corresponds to a different TCCON site; TCCON site colors and site IDs (see Table 3) are shown on the right). The comparison is based on 1446 monthly values. The mean difference (satellite – TCCON) is 0.18 ppm and the standard deviation of the difference is 1.18 ppm. The linear correlation coefficient R is 0.99.

Download

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f14

Figure 14September 2010 XCH4 at 10×10 spatial resolution showing (i) the individual sensor/algorithm input data sets (panels in rows 1–4; see Table 2 for details), (ii) EMMA XCH4 (bottom left), and (iii) the inter-algorithm spread (IAS, 1σ) as computed by EMMA (bottom right; see main text for details). Also shown in the bottom-right panel are the locations of the TCCON sites (pink triangles) and the range of IAS values covered by them (see color bar).

Figure 16 shows the TCCON station XCH4 bias anomaly as also shown for XCO2 in Fig. 12; i.e., Fig. 16 shows the biases as obtained with the EMMA method at the various TCCON sites used for the EMMA method comparisons. As for XCO2 not only the mean satellite–TCCON differences as obtained for the EMMA product are shown but also the differences for all the individual sensor/algorithm input products. The differences are shown as anomalies with respect to the mean; i.e., the sum of the differences in each row is zero. As can be seen from Fig. 16, the pattern of satellite–TCCON XCH4 differences has some similarity with the XCO2 difference pattern shown in Fig. 12. For example, the differences are mostly positive at Sodankylä and Garmisch-Partenkirchen and mostly negative at Darwin and Wollongong. But there are also significant differences, for example, with respect to the sign of the bias (e.g., Park Falls, Bremen, Karlsruhe).

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f15

Figure 15As Fig. 14 but for September 2018. Note that the SCIAMACHY/WFMD map (top left) is empty because this product ended in April 2012 (see Fig. 14 for SCIAMACHY/WFMD XCH4). For product GOSAT/PPDF (row 4) no data were available for this month (see Fig. 14 for GOSAT/PPDF XCH4).

Table 8Validation summary for data product XCH4_EMMA (version 4.1).

Download Print Version | Download XLSX

The XCH4_OBS4MIPS product has also been directly compared with TCCON (Fig. 17) using the same method as also used for product XCO2_OBS4MIPS (Fig. 13). As can be seen from Fig. 17, the mean difference (satellite – TCCON) is −2.88 ppb (which is close to the mean value of the global bias of −2.0 ppb of product XCH4_EMMA listed in Table 8), the standard deviation is 8.65 ppb (as expected, because of the averaging, which is somewhat smaller than the value of 17.4 ppb obtained for the XCH4_EMMA product listed in Table 8), and the linear correlation coefficient is 0.97.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f16

Figure 16As Fig. 12 but for XCH4, i.e., average XCH4 differences (satellite – TCCON) for the different satellite XCH4 products at 10 TCCON sites as used by the EMMA assessment method. The differences are shown as anomalies; i.e., the sum of the values corresponding to a given row is zero.

Download

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f17

Figure 17Summary of the comparison of product XCH4_OBS4MIPS with TCCON monthly mean XCH4. The comparison is based on 1518 monthly values. The mean difference (satellite – TCCON) is −2.88 ppb and the standard deviation of the difference is 8.65 ppb. The linear correlation coefficient R is 0.97.

Download

Figure 2 presents an overview of the XCH4 data product in terms of time series for three latitude bands and global maps. As can be seen, XCH4 was nearly constant during 2003–2006 (apart from seasonal fluctuations) but has been increasing since 2007 (for a discussion of the trend and annual growth rates see Sect. 4.3). The reason for this is likely a combination of increasing natural (e.g., wetlands) and anthropogenic (e.g., fossil fuel related) emissions and possibly decreasing sinks (hydroxyl, OH, radical), but it does not seem currently possible to be more definitive (e.g., Worden et al., 2017; Nisbet et al., 2019; Turner et al., 2019; Howarth, 2019; Schaefer, 2019).

4.3 Annual growth rates

Finally, we present an update and extension of the year 2003–2016 annual XCO2 growth rates shown in Buchwitz et al. (2018), using the new OBS4MIPS v4.1 XCO2 data set covering the time period 2003–2018 (Fig. 18). Figure 18a shows the time series of the globally averaged OBS4MIPS version 4.1 XCO2 data product over land. In contrast to Buchwitz et al. (2018), the analysis presented here is based on data over land only as this permits the generation of a time series with better internal consistency (note that the XCO2 OBS4MIPS product is land only for 2003–2008). The average growth rate during 2010–2018, i.e., for the time period where an ensemble of GOSAT and OCO-2 data has been used, is 2.28±0.04 ppm yr−1. As can be seen from Fig. 18b, the year 2017 and 2018 growth rates are less than the growth rates of the years 2015 and 2016, which were years with a strong El Niño. The XCO2 growth rates are in reasonable agreement with the global CO2 growth rates published by National Oceanic and Atmospheric Administration (NOAA) (shown in blue color in Fig. 18b), which are based on marine surface CO2 observations (ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_gr_gl.txt; last access: 30 July 2019). As can be seen from Fig. 18b, the agreement of the satellite-derived XCO2 growth rates with the NOAA surface-CO2-based growth rates is better from year 2010 onwards compared to the time period before when the EMMA data set consists only of one SCIAMACHY data set instead of the full ensemble. For 2018, the XCO2 growth rate is 2.1±0.5 ppm yr−1, which is lower than the NOAA surface CO2 growth rate of 2.43±0.08 ppm yr−1. Note that the 1σ uncertainty ranges of the two growth rate estimates overlap, which indicates that the two growth rate estimates are consistent.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f18

Figure 18(a) Monthly values of the globally averaged XCO2 (over land) as computed from the OBS4MIPS version 4.1 XCO2 data product. The corresponding annual mean XCO2 values are also listed. The increase during 2010–2018 is 2.28±0.04 ppm yr−1 as obtained via a linear fit. (b) Annual XCO2 growth rates (red, with 1σ uncertainties; the corresponding numerical values are also listed with 1σ uncertainty in brackets) and CO2 growth rates from NOAA (shown in blue) obtained from marine surface CO2 observations.

Download

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f19

Figure 19(a) Monthly values of the globally averaged XCH4 (over land) as computed from the OBS4MIPS version 4.1 XCH4 data product. The corresponding annual mean XCH4 values are also listed. The increase during 2010–2018 is 7.9±0.2 ppb yr−1 as obtained via a linear fit. (b) Annual XCH4 growth rates (red, with 1σ uncertainties; the corresponding numerical values are also listed with 1σ uncertainty in brackets) and CH4 growth rates from NOAA (shown in blue) obtained from marine surface CH4 observations.

Download

The growth rate of atmospheric methane is also an important quantity (e.g., Nisbet et al., 2019). The method of Buchwitz et al. (2018) has now also been used to compute annual XCH4 growth rates from satellite XCH4 retrievals. Figure 19a shows the time series of the globally averaged OBS4MIPS version 4.1 XCH4 data product over land. As shown by the linear fit, the average growth rate is 7.9±0.2 ppb yr−1 during 2010–2018, i.e., for the time period where an ensemble of GOSAT data has been used. The annual growth rates are shown in Fig. 19b for the satellite-derived XCH4 (red) and for the NOAA growth rates (ftp://aftp.cmdl.noaa.gov/products/trends/ch4/ch4_gr_gl.txt; last access: 30 July 2019) derived from marine surface CH4 observations. For 2018, the XCH4 growth rate is 10±6 ppb yr−1, which is close to the NOAA surface CH4 growth rate of 9.46±0.56 ppb yr−1.

5 Summary and conclusions

Satellite-derived ensemble XCO2 and XCH4 data products have been generated and validated. These data products are the version 4.1 Level 2 (L2) products XCO2_EMMA and XCH4_EMMA and the Level 3 (L3) products XCO2_OBS4MIPS and XCH4_OBS4MIPS and cover the time period 2003–2018. The data products are freely available for interested users via the Copernicus Climate Data Store (CDS, https://cds.climate.copernicus.eu/, last access: 10 January 2020), where also earlier versions of these data products are accessible. The L2 products have been generated with an adapted version of the EMMA algorithm (Reuter et al., 2013), and the L3 products have been generated by gridding (averaging) the EMMA L2 product to obtain products at monthly time and 5×5 spatial resolution in Obs4MIPS format. The products have been validated by comparisons with TCCON ground-based XCO2 and XCH4 retrievals using TCCON version GGG2014.

From January 2003 to March 2009 the products are based on SCIAMACHY/ENVISAT, and from April 2009 onwards the products use an ensemble of one SCIAMACHY (until early 2012) and several GOSAT products. The XCO2 products contain in addition L2 products from NASA's OCO-2 mission from September 2014 onwards.

The EMMA algorithm selects for each month and each 10×10 grid cell one of the available products, i.e., one from the existing ensemble of L2 input products, and transfers all relevant information (including averaging kernel etc.) from the selected L2 input product into the merged EMMA L2 product. The selected product is the median product. The main purpose of EMMA is to generate a Level 2 product, which covers an as-long-as-possible time series (longer than any of the individual sensor input data sets) with as-high-as-possible accuracy including all information needed, e.g., for surface flux inverse modeling. The median approach helps to reduce the occurrence of potential outliers and thus reduces spatial and temporal biases in the generated data products.

Detailed quality assessment results based on comparisons with TCCON ground-based retrievals have been presented. We found that the XCO2 Level 2 data set at the TCCON validation sites can be characterized by the following figures of merit (the corresponding values for the Level 3 product are listed in brackets) – single-observation random error (1σ): 1.29 ppm (monthly: 1.18 ppm); global bias: 0.20 ppm (0.18 ppm); and spatiotemporal bias or relative accuracy (1σ): 0.66 ppm (0.70 ppm). The corresponding values for the XCH4 products are single-observation random error (1σ): 17.4 ppb (monthly: 8.7 ppb); global bias: −2.0 ppb (−2.9 ppb), spatiotemporal bias (1σ): 5.0 ppb (4.9 ppb). It has also been found that the data products exhibit very good long-term stability as no significant linear bias trends have been identified.

The new data sets have also been used to derive annual XCO2 and XCH4 growth rates, which are in reasonable to good agreement with growth rates from the National Oceanic and Atmospheric Administration (NOAA) based on marine surface observations.

An important application for the EMMA products is to use them together with inverse modeling to obtain improved information on regional-scale CO2 (e.g., Houweling et al., 2015) and CH4 (e.g., Alexe et al., 2015) surface fluxes. Applications for the corresponding OBS4MIPS products are, for example, climate model comparisons (e.g., Lauer et al., 2017) and studies related to annual growth rates (e.g., Buchwitz et al., 2018). It is however important to note that these merged products are not necessarily the most optimal products for all applications as they do not contain all data from a given satellite sensor. For example, users interested primarily in emissions from power plants or other localized CO2 sources will prefer the original OCO-2 Level 2 data product (e.g., Nassar et al., 2017; Reuter et al., 2019a). Especially for users interested in only parts of the time series it is recommended to use the individual sensor products in addition to the merged product as this may significantly increase the robustness, reliability, and uncertainty characterization of key findings.

Appendix A

In this appendix, detailed validation results are shown for the individual sensor and EMMA XCO2 and XCH4 Level 2 data products.

The comparison of the various XCO2 data products with TCCON XCO2 at 10 TCCON sites is shown in Fig. A1. These 10 TCCON sites fulfill the EMMA criteria in terms of a sufficiently large number of colocations as defined to obtain robust conclusions per site. The individual soundings of the EMMA XCO2 product are shown as white circles with a black border. As can be seen, they are located within (mostly close to the center of) the range of values of the individual sensor/algorithm XCO2 values, which is expected.

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f20

Figure A1XCO2 time series at 10 TCCON sites during January 2009–December 2018 as obtained using the EMMA quality assessment method. TCCON GGG2014 XCO2 is shown as thick black dots, the individual satellite L2 input products are shown as colored dots, and the EMMA product is shown as white circles with black borders. The derived numerical values are listed in Table 4.

Download

https://www.atmos-meas-tech.net/13/789/2020/amt-13-789-2020-f21

Figure A2XCH4 time series at 10 TCCON sites during April 2010–December 2018 as obtained using the EMMA quality assessment method. TCCON GGG2014 XCH4 is shown as thick black dots, the individual satellite L2 input products are shown as colored dots, and the EMMA product is shown as a white circles with black borders. The derived numerical values are listed in Table 7.

Download

Figure A2 shows the comparison of the EMMA XCH4 product (white circles with a black border) and of the individual sensor XCH4 input products with TCCON XCH4 originating from the EMMA assessment method. As for the EMMA XCO2 product (Fig. A1), the EMMA XCH4 is located near the center of the clouds of XCH4 values, as expected.

Data availability

The EMMA and OBS4MIPS XCO2 and XCH4 version 4.1 data products (but also several data sets used as input; see data sets with CCI/C3S product ID in Tables 1 and 2) are available (from early 2020 onwards) via the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, ECMWF, 2020a) Climate Data Store (CDS, https://cds.climate.copernicus.eu/, ECMWF, 2020b), including documentation such as the product user guides (http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PUGS/C3S_D312b_Lot2.3.2.3-v1.0_PUGS-GHG_MAIN_v3.1.pdf, Buchwitz et al., 2019c; http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PUGS/C3S_D312b_Lot2.3.2.3-v1.0_PUGS-GHG_ANNEX-D_v3.1.pdf, Reuter et al., 2019d).

Author contributions

MR generated the EMMA and OBS4MIPS XCO2 and XCH4 version 4.1 data sets. MR and MB performed the data analysis. MB wrote the first version of the paper with support of MR. The following authors provided input data or expertise on data sets: MR, MB, OS, SN, HB, JPB, HBoe, ADN, JA, RJP, PS, LW, OPH, IA, AK, HS, KS, YY, IM, DC, CWO'D, JN, CP, TW, VAV, NMD, DWTG, RK, DP, FH, RS, YVT, KS, SR, MKS, MDM, DGF, LTI, CMR, CR, and DS. All authors contributed to significantly improve the paper.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

The generation of the EMMA Level 2 and OBS4MIPS Level 3 data sets and the corresponding data analysis has been funded primarily by the European Union (EU) via the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, last access: 10 January 2020) managed by the European Centre for Medium-Range Weather Forecasts (ECMWF).

The work presented here strongly benefited from additional funding by the European Space Agency (ESA) via ESA's Climate Change Initiative (CCI, http://www.esa-ghg-cci.org/, last access: 10 January 2020) projects GHG-CCI/GHG-CCI+.

The further development of the FOCAL retrieval algorithm used to generate the OCO-2/FOCAL XCO2 input data product would not have been possible without cofunding from the European Union (EU) Horizon 2020 (H2020) research and innovation program projects CHE (grant agreement no. 776186) and VERIFY (grant agreement no. 776810). The generation of the XCO2_OBS4MIPS product also benefited from cofunding from EU H2020 project CCiCC (grant agreement no. 821003).

We thank several space agencies for making available satellite Level 1 (L1) input data: ESA/DLR for SCIAMACHY L1 data, JAXA for GOSAT Level 1B data, and NASA for the OCO-2 L1 data product. We also thank ESA for making the GOSAT L1 product available via the ESA Third Party Mission (TPM) archive.

We thank NIES for the operational GOSAT XCO2 and XCH4 Level 2 products (obtained from: https://data2.gosat.nies.go.jp/, last access: 4- September 2019) and the NASA team for the GOSAT and OCO-2 ACOS Level 2 XCO2 products (the NASA GOSAT L2 data product was obtained from: https://oco2.gesdisc.eosdis.nasa.gov/data/GOSAT_TANSO_Level2/ACOS_L2_Lite_FP.7.3/, last access: 4 September 2019; the NASA OCO-2 data product was obtained from: https://oco2.gesdisc.eosdis.nasa.gov/data/s4pa/OCO2_DATA/OCO2_L2_Lite_FP.9r/, last access: 4 September 2019).

TCCON data were obtained from the TCCON Data Archive, hosted by CaltechDATA, California Institute of Technology (https://tccondata.org/, last access: 15 July 2019).

The TCCON stations Ascension Island, Bremen, Garmisch, Karlsruhe, and Ny-Ålesund have been supported by the German Bundesministerium für Wirtschaft und Energie (BMWi) via the DLR under grants 50EE1711A-E. We thank the ESA Ariane tracking station at North East Bay, Ascension Island, for hosting and local support. Nicholas M. Deutscher is supported by an ARC Future Fellowship, FT180100327. The TCCON site at Réunion island is operated by the Royal Belgian Institute for Space Aeronomy with financial support in 2014, 2015, 2016, 2017, 2018, and 2019 under the EU project ICOS-Inwire and the ministerial decree for ICOS (FR/35/IC2) and local activities supported by LACy/UMR8105 – Université de La Réunion. The TCCON stations at Tsukuba and Burgos are supported in part by the GOSAT series project. Local support for Burgos is provided by the Energy Development Corporation (EDC, Philippines). The Paris TCCON site has received funding from Sorbonne Université, the French research center CNRS, the French space agency CNES, and Région Île-de-France.

We also thank NOAA for access to the surface CO2 (file ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_gr_gl.txt; last access: 30 July 2019) and CH4 (file ftp://aftp.cmdl.noaa.gov/products/trends/ch4/ch4_gr_gl.txt; last access: 30 July 2019) growth rate data sets. Output from NOAA's CarbonTracker has been used as input for the SECM2018 model. CarbonTracker CT2017 results are provided by NOAA ESRL, Boulder, Colorado, USA, from their website at: http://carbontracker.noaa.gov/ (last access: 4 September 2019).

We also thank Peter Bergamaschi for providing MACC-II project inversion system CH4 fields, which have been used as input for the SC4C2018 model.

Financial support

This research has been supported by the European Union (Copernicus Climate Change Service project C3S_312b_Lot2, via contract with DLR, contract no. D/565/67260504; CHE, grant agreement no. 776186; VERIFY, grant agreement no. 776810; CCiCC, grant agreement no. 821003) and the European Space Agency (ESA; grant no. Project GHG-CCI+).

The article processing charges for this open-access publication were covered by the University of Bremen.

Review statement

This paper was edited by John Worden and reviewed by Ray Nassar and two anonymous referees.

References

Alexe, M., Bergamaschi, P., Segers, A., Detmers, R., Butz, A., Hasekamp, O., Guerlet, S., Parker, R., Boesch, H., Frankenberg, C., Scheepmaker, R. A., Dlugokencky, E., Sweeney, C., Wofsy, S. C., and Kort, E. A.: Inverse modelling of CH4 emissions for 2010–2011 using different satellite retrieval products from GOSAT and SCIAMACHY, Atmos. Chem. Phys., 15, 113–133, https://doi.org/10.5194/acp-15-113-2015, 2015. 

Basu, S., Guerlet, S., Butz, A., Houweling, S., Hasekamp, O., Aben, I., Krummel, P., Steele, P., Langenfelds, R., Torn, M., Biraud, S., Stephens, B., Andrews, A., and Worthy, D.: Global CO2 fluxes estimated from GOSAT retrievals of total column CO2, Atmos. Chem. Phys., 13, 8695–8717, https://doi.org/10.5194/acp-13-8695-2013, 2013. 

Boesch, H., Anand, J., Di Noia, A., Buchwitz, M., Somkuti, P., and Parker, R.: Product Quality Assessment Report (PQAR) – ANNEX A for products CO2_GOS_OCFP, CH4_GOS_OCFP & CH4_ GOS_OCPR (v7.2, 2009–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 44 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PQAR/C3S_D312b_Lot2.2.3.2-v1.0_PQAR-GHG_ANNEX-A_v3.1.pdf (last access: 10 January 2020), 2019. 

Bovensmann, H., Burrows, J. P., Buchwitz, M., Frerick, J., Noël, S., Rozanov, V. V., Chance, K. V., and Goede, A. H. P.: SCIAMACHY – Mission objectives and measurement modes, J. Atmos. Sci., 56, 127–150, 1999. 

Bril, A., Oshchepkov, S., and Yokota, T.: Application of a probability density function-based atmospheric light-scattering correction to carbon dioxide retrievals from GOSAT over-sea observations, Remote Sens. Environ., 117, 301–306, 2012. 

Buchwitz, M., Reuter, M., Schneising, O., Boesch, H., Guerlet, S., Dils, B., Aben, I., Armante, R., Bergamaschi, P., Blumenstock, T., Bovensmann, H., Brunner, D., Buchmann, B., Burrows, J. P., Butz, A., Chédin, A., Chevallier, F., Crevoisier, C. D., Deutscher, N. M., Frankenberg, C., Hase, F., Hasekamp, O. P., Heymann, J., Kaminski, T., Laeng, A., Lichtenberg, G., De Mazière, M., Noël, S., Notholt, J., Orphal, J., Popp, C., Parker, R., Scholze, M., Sussmann, R., Stiller, G. P., Warneke, T., Zehner, C., Bril, A., Crisp, D., Griffith, D. W. T., Kuze, A., O'Dell, C., Oshchepkov, S., Sherlock, V., Suto, H., Wennberg, P., Wunch, D., Yokota, T., and Yoshida, Y.: The Greenhouse Gas Climate Change Initiative (GHG-CCI): comparison and quality assessment of near-surface-sensitive satellite-derived CO2 and CH4 global data sets, Remote Sens. Environ., 162, 344–362, https://doi.org/10.1016/j.rse.2013.04.024, 2015. 

Buchwitz, M., Schneising, O., Reuter, M., Heymann, J., Krautwurst, S., Bovensmann, H., Burrows, J. P., Boesch, H., Parker, R. J., Somkuti, P., Detmers, R. G., Hasekamp, O. P., Aben, I., Butz, A., Frankenberg, C., and Turner, A. J.: Satellite-derived methane hotspot emission estimates using a fast data-driven method, Atmos. Chem. Phys., 17, 5751–5774, https://doi.org/10.5194/acp-17-5751-2017, 2017a. 

Buchwitz, M., Reuter, M., Schneising, O., Hewson, W., Detmers, R. G., Boesch, H., Hasekamp, O. P., Aben, I., Bovensmann, H., Burrows, J. P., Butz, A., Chevallier, F., Dils, B., Frankenberg, C., Heymann, J., Lichtenberg, G., De Maziere, M., Notholt, J., Parker, R., Warneke, T., Zehner, C., Griffith, D. W. T., Deutscher, N. M., Kuze, A., Suto, H., and Wunch, D.: Global satellite observations of column-averaged carbon dioxide and methane: The GHG-CCI XCO2 and XCH4 CRDP3 data set, Remote Sens. Environ., 203, 276–295, https://doi.org/10.1016/j.rse.2016.12.027, 2017b. 

Buchwitz, M., Reuter, M., Schneising, O., Noël, S., Gier, B., Bovensmann, H., Burrows, J. P., Boesch, H., Anand, J., Parker, R. J., Somkuti, P., Detmers, R. G., Hasekamp, O. P., Aben, I., Butz, A., Kuze, A., Suto, H., Yoshida, Y., Crisp, D., and O'Dell, C.: Computation and analysis of atmospheric carbon dioxide annual mean growth rates from satellite observations during 2003–2016, Atmos. Chem. Phys., 18, 17355–17370, https://doi.org/10.5194/acp-18-17355-2018, 2018. 

Buchwitz, M., Reuter, M., Schneising-Weigel, O., Aben, I., Wu, L., Hasekamp, O. P., Boesch, H., Di Noia, A., Crevoisier, C., and Armante, R.: Product Quality Assessment Report (PQAR) – Main document for Greenhouse Gas (GHG: CO2 & CH4) data set CDR 3 (2003–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 103 pp., availabel at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PQAR/C3S_D312b_Lot2.2.3.2-v1.0_PQAR-GHG_MAIN_v3.1.pdf (last access: 10 January 2020), 2019a. 

Buchwitz, M., Reuter, M., Schneising-Weigel, O., Aben, I., Wu, L., Hasekamp, O. P., Boesch, H., Di Noia, A., Crevoisier, C., and Armante, R.: Algorithm Theoretical Basis Document (ATBD) – Main document for Greenhouse Gas (GHG: CO2 & CH4) data set CDR 3 (2003–2018), Technical Report Copernicus Climate Change Service (C3S) (main document and 5 Annexes), version 3.1, 03-11-2019, 43 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/ATBD/C3S_D312b_Lot2.1.3.2-v1.0_ATBD-GHG_MAIN_v3.1.pdf (last access: 10 January 2020), 2019b. 

Buchwitz, M., Reuter, M., Schneising-Weigel, O., Aben, I., Wu, L., Hasekamp, O. P., Boesch, H., Di Noia, A., Crevoisier, C., and Armante, R.: Product User Guide and Specification (PUGS) – Main document for Greenhouse Gas (GHG: CO2 & CH4) data set CDR 3 (2003–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 97 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PUGS/C3S_D312b_Lot2.3.2.3-v1.0_PUGS-GHG_MAIN_v3.1.pdf (last access: 10 January 2020), 2019c. 

Burrows, J. P., Hölzle, E., Goede, A. P. H., Visser, H., and Fricke, W.: SCIAMACHY–Scanning Imaging Absorption Spectrometer for Atmospheric Chartography, Acta Astronaut., 35, 445–451, https://doi.org/10.1016/0094-5765(94)00278-t, 1995. 

Butz, A., Hasekamp, O.P., Frankenberg, C., Vidot, J., and Aben, I.: CH4 retrievals from space-based solar backscatter measurements: Performance evaluation against simulated aerosol and cirrus loaded scenes, J. Geophys. Res., 115, D24302, https://doi.org/10.1029/2010JD014514, 2010. 

Butz, A., Guerlet, S., Hasekamp, O., Schepers, D., Galli, A., Aben, I., Frankenberg, C., Hartmann, J.-M., Tran, H., Kuze, A., Keppel-Aleks, G., Toon, G., Wunch, D., Wennberg, P., Deutscher, N., Griffith, D., Macatangay, R., Messerschmidt, J., Notholt, J., and Warneke, T.: Toward accurate CO2 and CH4 observations from GOSAT, Geophys. Res. Lett., 38, L14812, https://doi.org/10.1029/2011GL047888, 2011. 

Chevallier, F.: On the statistical optimality of CO2 atmospheric inversions assimilating CO2 column retrievals, Atmos. Chem. Phys., 15, 11133–11145, https://doi.org/10.5194/acp-15-11133-2015, 2015. 

Chevallier, F., Palmer, P. I., Feng, L., Boesch, H., O'Dell, C. W., and Bousquet, P.: Towards robust and consistent regional CO2 flux estimates from in situ and space-borne measurements of atmospheric CO2, Geophys. Res. Lett., 41, 1065–1070, https://doi.org/10.1002/2013GL058772, 2014. 

Cogan, A. J., Boesch, H., Parker, R. J., Feng, L., Palmer, P. I., Blavier, J.-F. L., Deutscher, N. M., Macatangay, R., Notholt, J., Roehl, C., Warneke, T., and Wunsch, D.: Atmospheric carbon dioxide retrieved from the Greenhouse gases Observing SATellite (GOSAT): Comparison with ground-based TCCON observations and GEOS-Chem model calculations, J. Geophys. Res., 117, D21301, https://doi.org/10.1029/2012JD018087, 2012. 

Crisp, D., Atlas, R. M., Bréon, F.-M., Brown, L. R., Burrows, J. P., Ciais, P., Connor, B. J., Doney, S. C., Fung, I. Y., Jacob, D. J., Miller, C. E., O'Brien, D., Pawson, S., Randerson, J. T., Rayner, P., Salawitch, R. S., Sander, S. P., Sen, B., Stephens, G. L., Tans, P. P., Toon, G. C.,Wennberg, P. O.,Wofsy, S. C., Yung, Y. L., Kuang, Z., Chudasama, B., Sprague, G.,Weiss, P., Pollock, R., Kenyon, D., and Schroll, S.: The Orbiting Carbon Observatory (OCO) mission, Adv. Space Res., 34, 700–709, 2004. 

De Mazière, M., Sha, M. K., Desmet, F., Hermans, C., Scolas, F., Kumps, N., Metzger, J.-M., Duflot, V., and Cammas, J.-P.: TCCON data from Réunion Island (RE), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.reunion01.R1, 2017. 

Deutscher, N. M., Notholt, J., Messerschmidt, J., Weinzierl, C., Warneke, T., Petri, C., Grupe, P., and Katrynski, K.: TCCON data from Bialystok (PL), Release GGG2014R2. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.bialystok01.R2, 2019. 

Dils, B., Buchwitz, M., Reuter, M., Schneising, O., Boesch, H., Parker, R., Guerlet, S., Aben, I., Blumenstock, T., Burrows, J. P., Butz, A., Deutscher, N. M., Frankenberg, C., Hase, F., Hasekamp, O. P., Heymann, J., De Mazière, M., Notholt, J., Sussmann, R., Warneke, T., Griffith, D., Sherlock, V., and Wunch, D.: The Greenhouse Gas Climate Change Initiative (GHG-CCI): comparative validation of GHG-CCI SCIAMACHY/ENVISAT and TANSO-FTS/GOSAT CO2 and CH4 retrieval algorithm products with measurements from the TCCON, Atmos. Meas. Tech., 7, 1723–1744, https://doi.org/10.5194/amt-7-1723-2014, 2014. 

ECMWF: Copernicus Climate Change Service (C3S) website, available at: https://climate.copernicus.eu/, last access: 10 January 2020a. 

ECMWF: Copernicus Climate Data Store (CDS) website, available at: https://cds.climate.copernicus.eu/, last access: 10 January 2020b. 

Feist, D. G., Arnold, S. G., John, N., and Geibel, M. C.: TCCON data from Ascension Island (SH), Release GGG2014R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.ascension01.R0/1149285, 2014. 

Ganesan, A. L., Rigby, M., Lunt, M. F., Parker, R. J., Boesch, H., Goulding, N., Umezawa. T., Zahn, A., Chatterjee, A., Prinn, R. G., Tiwari, Y. K., van der Schoot, M., and Krummel, P. B.: Atmospheric observations show accurate reporting and little growth in India's methane emissions, Nat. Commun., 8, 836, https://doi.org/10.1038/s41467-017-00994, 2017. 

Gaubert, B., Stephens, B. B., Basu, S., Chevallier, F., Deng, F., Kort, E. A., Patra, P. K., Peters, W., Rödenbeck, C., Saeki, T., Schimel, D., Van der Laan-Luijkx, I., Wofsy, S., and Yin, Y.: Global atmospheric CO2 inverse models converging on neutral tropical land exchange, but disagreeing on fossil fuel and atmospheric growth rate, Biogeosciences, 16, 117–134, https://doi.org/10.5194/bg-16-117-2019, 2019. 

GCOS-154: Global Climate Observing System (GCOS), Systematic observation requirements for satellite-based products for climate, Supplemental details to the satellite-based component of the “Implementation Plan for the Global Observing System for Climate in Support of the UNFCCC (2010 update)”, prepared by: World Meteorological Organization (WMO), Intergovernmental Oceanographic Commission, United Nations Environment Programme (UNEP), International Council for Science, Doc.: GCOS 154, availabel at: https://www.wmo.int/pages/prog/gcos/Publications/gcos-154.pdf (last access: 21 February 2019), 2010. 

GCOS-200: The Global Observing System for Climate: Implementation Needs, World Meteorological Organization (WMO), GCOS-200 (GOOS-214), 325 pp., available at: http://unfccc.int/files/science/workstreams/systematic_observation/application/pdf/gcos_ip_10oct2016.pdf (last access: 21 February 2019), 2016. 

Griffith, D. W., Velazco, V. A., Deutscher, N. M., Murphy, C., Jones, N. B., Wilson, S. R., Macatangay, R. C., Kettlewell, G. C., Buchholz, R. R., and Riggenbach, M. O.: TCCON data from Wollongong, (AU), Release GGG2014R0. TCCON data archive, hosted by Caltech-DATA, https://doi.org/10.14291/tccon.ggg2014.wollongong01.R0/1149291, 2014a. 

Griffith, D. W. T., Deutscher, N. M., Velazco, V. A.,Wennberg, P. O., Yavin, Y., Keppel-Aleks, G.,Washenfelder, R. A., Toon, G. C., Blavier, J.-F., Murphy, C., Jones, N. B., Kettlewell, G. C., Connor, B. J., Macatangay, R. C., Roehl, C., Ryczek, M., Glowacki, J., Culgan, T., and Bryant, G. W.: TCCON data from Darwin (AU), Release GGG2014R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.darwin01.R0/1149290, 2014b. 

Hase, F., Blumenstock, T., Dohe, S., Groß, J., and Kiel, M.: TCCON data from Karlsruhe (DE), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.karlsruhe01.R1/1182416, 2015. 

Hayman, G. D., O'Connor, F. M., Dalvi, M., Clark, D. B., Gedney, N., Huntingford, C., Prigent, C., Buchwitz, M., Schneising, O., Burrows, J. P., Wilson, C., Richards, N., and Chipperfield, M.: Comparison of the HadGEM2 climate-chemistry model against in situ and SCIAMACHY atmospheric methane data, Atmos. Chem. Phys., 14, 13257–13280, https://doi.org/10.5194/acp-14-13257-2014, 2014. 

Hedelius, J. K., Parker, H., Wunch, D., Roehl, C. M., Viatte, C., Newman, S., Toon, G. C., Podolske, J. R., Hillyard, P. W., Iraci, L. T., Dubey, M. K., and Wennberg, P. O.: Intercomparability of XCO2 and XCH4 from the United States TCCON sites, Atmos. Meas. Tech., 10, 1481–1493, https://doi.org/10.5194/amt-10-1481-2017, 2017. 

Heymann, J., Reuter, M., Hilker, M., Buchwitz, M., Schneising, O., Bovensmann, H., Burrows, J. P., Kuze, A., Suto, H., Deutscher, N. M., Dubey, M. K., Griffith, D. W. T., Hase, F., Kawakami, S., Kivi, R., Morino, I., Petri, C., Roehl, C., Schneider, M., Sherlock, V., Sussmann, R., Velazco, V. A., Warneke, T., and Wunch, D.: Consistent satellite XCO2 retrievals from SCIAMACHY and GOSAT using the BESD algorithm, Atmos. Meas. Tech., 8, 2961–2980, https://doi.org/10.5194/amt-8-2961-2015, 2015. 

Hollmann, R., Merchant, C. J., Saunders, R., Downy, C., Buchwitz, M., Cazenave, A., Chuvieco, E., Defourny, P., de Leeuw, G., Forsberg, R., Holzer-Popp, T., Paul, F., Sandven, S., Sathyendranath, S., van Roozendael, M., and Wagner, W.: The ESA Climate Change Initiative: satellite data records for essential climate variables, B. Am. Meteorol. Soc., 94, 1541–1552, https://doi.org/10.1175/BAMS-D-11-00254.1, 2013. 

Houweling, S., Baker, D., Basu, S., Boesch, H., Butz, A., Chevallier, F., Deng, F., Dlugokencky, E. J., Feng, L., Ganshin, A., Hasekamp, O., Jones, D., Maksyutov, S., Marshall, J., Oda, T., O'Dell, C. W., Oshchepkov, S., Palmer, P. I., Peylin, P., Poussi, Z., Reum, F., Takagi, H., Yoshida, Y., and Zhuralev, R.: An intercomparison of inverse models for estimating sources and sinks of CO2 using GOSAT measurements, J. Geophys. Res.-Atmos., 120, 5253–5266, https://doi.org/10.1002/2014JD022962, 2015. 

Howarth, R. W.: Ideas and perspectives: is shale gas a major driver of recent increase in global atmospheric methane?, Biogeosciences, 16, 3033–3046, https://doi.org/10.5194/bg-16-3033-2019, 2019. 

Hu, H., Landgraf, J., Detmers, R., Borsdorff, T., Aan de Brugh, J., Aben, I., Butz, A., and Hasekamp, O.: Toward Global Mapping of Methane With TROPOMI: First Results and Intersatellite Comparison to GOSAT, Geophys. Res. Lett., 45, 3682–3689, https://doi.org/10.1002/2018GL077259, 2018. 

IPCC: Climate Change 2013: The Physical Science Basis, Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Report on Climate Change, Cambridge University Press, availabel at: http://www.ipcc.ch/report/ar5/wg1/ (last access: 21 February 2019), 2013. 

Iraci, L., Podolske, J., Hillyard, P., Roehl, C., Wennberg, P. O., Blavier, J.-F., Landeros, J., Allen, N., Wunch, D., Zavaleta, J., Quigley, E., Osterman, G., Albertson, R., Dunwoody, K., and Boyden, H.: TCCON data from Edwards (US), Release GGG2014.R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.edwards01.R1/1255068, 2014. 

Kaminski, T., Scholze, M., Voßbeck, M., Knorr, W., Buchwitz, M., and Reuter, M.: Constraining a terrestrial biosphere model with remotely sensed atmospheric carbon dioxide, Remote Sens. Environ., 203, 109–124, 2017. 

Kiel, M., O'Dell, C. W., Fisher, B., Eldering, A., Nassar, R., MacDonald, C. G., and Wennberg, P. O.: How bias correction goes wrong: measurement of XCO2 affected by erroneous surface pressure estimates, Atmos. Meas. Tech., 12, 2241–2259, https://doi.org/10.5194/amt-12-2241-2019, 2019. 

Kivi, R. and Heikkinen, P.: Fourier transform spectrometer measurements of column CO2 at Sodankylä, Finland, Geosci. Instrum. Method. Data Syst., 5, 271–279, https://doi.org/10.5194/gi-5-271-2016, 2016. 

Kivi, R., Heikkinen, P., and Kyrö, E.: TCCON data from Sodankylä (FI), Release GGG2014R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.sodankyla01.R0/1149280, 2014. 

Kuze, A., Suto, H., Shiomi, K., Kawakami, S., Tanaka, M., Ueda, Y., Deguchi, A., Yoshida, J., Yamamoto, Y., Kataoka, F., Taylor, T. E., and Buijs, H. L.: Update on GOSAT TANSO-FTS performance, operations, and data products after more than 6 years in space, Atmos. Meas. Tech., 9, 2445–2461, https://doi.org/10.5194/amt-9-2445-2016, 2016. 

Lauer, A., Eyring, V., Righi, M., Buchwitz, M., Defourny, P., Evaldsson, M., Friedlingstein, P., de Jeu, R., de Leeuw, G., Loew, A., Merchant, C. J., Müller, B., Popp, T., Reuter, M., Sandven, S., Senftleben, D., Stengel, M., Van Roozendael, M., Wenzel, S., and Willén, U.: Benchmarking CMIP5 models with a subset of ESA CCI Phase 2 data using the ESMValTool, Remote Sens. Environ., 203, 9–39, https://doi.org/10.1016/j.rse.2017.01.007, 2017. 

Le Quéré, C., Andrew, R. M., Friedlingstein, P., Sitch, S., Pongratz, J., Manning, A. C., Korsbakken, J. I., Peters, G. P., Canadell, J. G., Jackson, R. B., Boden, T. A., Tans, P. P., Andrews, O. D., Arora, V. K., Bakker, D. C. E., Barbero, L., Becker, M., Betts, R. A., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Cosca, C. E., Cross, J., Currie, K., Gasser, T., Harris, I., Hauck, J., Haverd, V., Houghton, R. A., Hunt, C. W., Hurtt, G., Ilyina, T., Jain, A. K., Kato, E., Kautz, M., Keeling, R. F., Klein Goldewijk, K., Körtzinger, A., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Lima, I., Lombardozzi, D., Metzl, N., Millero, F., Monteiro, P. M. S., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., Nojiri, Y., Padin, X. A., Peregon, A., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Reimer, J., Rödenbeck, C., Schwinger, J., Séférian, R., Skjelvan, I., Stocker, B. D., Tian, H., Tilbrook, B., Tubiello, F. N., van der Laan-Luijkx, I. T., van der Werf, G. R., van Heuven, S., Viovy, N., Vuichard, N., Walker, A. P., Watson, A. J., Wiltshire, A. J., Zaehle, S., and Zhu, D.: Global Carbon Budget 2017, Earth Syst. Sci. Data, 10, 405–448, https://doi.org/10.5194/essd-10-405-2018, 2018. 

Liu, J., Bowman, K. W., Schimel, D. S., Parazoo, N. C., Jiang, Z., Lee, M., Bloom, A. A., Wunch, D., Frankenberg, C., Sun, Y., O'Dell, C. W., Gurney, K. R., Menemenlis, D., Gierach, M., Crisp, D., and Eldering, A.: Contrasting carbon cycle responses of the tropical continents to the 2015–2016 El Niño, Science, 358, eaam5690, https://doi.org/10.1126/science.aam5690, 2017. 

Maasakkers, J. D., Jacob, D. J., Sulprizio, M. P., Scarpelli, T. R., Nesser, H., Sheng, J.-X., Zhang, Y., Hersher, M., Bloom, A. A., Bowman, K. W., Worden, J. R., Janssens-Maenhout, G., and Parker, R. J.: Global distribution of methane emissions, emission trends, and OH concentrations and trends inferred from an inversion of GOSAT satellite data for 2010–2015, Atmos. Chem. Phys., 19, 7859–7881, https://doi.org/10.5194/acp-19-7859-2019, 2019. 

Miller, S. M., Michalak, A. M., Detmers, R. G., Hasekamp, O. P., Bruhwiler, L. M. P., and Schwietzke, S.: China's coal mine methane regulations have not curbed growing emissions, Nat. Commun., 10, 303, 2019. 

Morino, I., Matsuzaki, T., and Shishime, A.: TCCON data from Tsukuba, Ibaraki (JP), 125HR, Release GGG2014R2. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.tsukuba02.R2, 2018a. 

Morino, I., Velazco, V. A., Hori, A., Uchino, O., and Griffith, D. W. T.: TCCON data from Burgos, Ilocos Norte (PH), Release GGG2014.R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.burgos01.R0, 2018b. 

Nassar, R., Hill, T. G., McLinden, C. A., Wunch, D., Jones, D. B. A., and Crisp, D.: Quantifying CO2 emissions from individual power plants from space, Geophys. Res. Lett., 44, 10045–10053, https://doi.org/10.1002/2017GL074702, 2017. 

Nisbet, E. G., Manning, M. R., Dlugokencky, E. J., Fisher, R. E., Lowry, D., Michel, S. E., Lund Myhre, C., Platt, S. M., Allen, G., Bousquet, P., Brownlow, R., Cain, M., France, J. L., Hermansen, O., Hossaini, R., Jones, A. E., Levin, I., Manning, A. C., Myhre, G., Pyle, J. A., Vaughn, B., Warwick, N. J., and White, J. W. C.: Very strong atmospheric methane growth in the four years 2014–2017: Implications for the Paris Agreement, Global Biogeochem. Cy., 33, 318–342, https://doi.org/10.1029/2018GB006009, 2019. 

Notholt, J., Schrems, O., Warneke, T., Deutscher, N. M., Weinzierl, C., Palm, M., and Buschmann, M.: TCCON data from Ny-Ålesund, Spitzbergen (NO), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.nyalesund01.R1, 2019a. 

Notholt, J., Petri, C., Warneke, T., Deutscher, N., Buschmann, M., Weinzierl, C., Macatangay, R., and Gruppe, P.: TCCON data from Bremen (DE), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.bremen01.R1, 2019b. 

O'Dell, C. W., Connor, B., Bösch, H., O'Brien, D., Frankenberg, C., Castano, R., Christi, M., Eldering, D., Fisher, B., Gunson, M., McDuffie, J., Miller, C. E., Natraj, V., Oyafuso, F., Polonsky, I., Smyth, M., Taylor, T., Toon, G. C., Wennberg, P. O., and Wunch, D.: The ACOS CO2 retrieval algorithm – Part 1: Description and validation against synthetic observations, Atmos. Meas. Tech., 5, 99–121, https://doi.org/10.5194/amt-5-99-2012, 2012. 

O'Dell, C. W., Eldering, A., Wennberg, P. O., Crisp, D., Gunson, M. R., Fisher, B., Frankenberg, C., Kiel, M., Lindqvist, H., Mandrake, L., Merrelli, A., Natraj, V., Nelson, R. R., Osterman, G. B., Payne, V. H., Taylor, T. E., Wunch, D., Drouin, B. J., Oyafuso, F., Chang, A., McDuffie, J., Smyth, M., Baker, D. F., Basu, S., Chevallier, F., Crowell, S. M. R., Feng, L., Palmer, P. I., Dubey, M., García, O. E., Griffith, D. W. T., Hase, F., Iraci, L. T., Kivi, R., Morino, I., Notholt, J., Ohyama, H., Petri, C., Roehl, C. M., Sha, M. K., Strong, K., Sussmann, R., Te, Y., Uchino, O., and Velazco, V. A.: Improved retrievals of carbon dioxide from Orbiting Carbon Observatory-2 with the version 8 ACOS algorithm, Atmos. Meas. Tech., 11, 6539–6576, https://doi.org/10.5194/amt-11-6539-2018, 2018. 

Parker, R., Boesch, H., Cogan, A., Fraser, A., Feng, L, Palmer, P., Messerschmidt, J., Deutscher, N., Griffth, D., Notholt, J., Wennberg, P. O., and Wunch, D.: Methane Observations from the Greenhouse gases Observing SATellite: Comparison to ground-based TCCON data and Model Calculations, Geophys. Res. Lett., 38, L15807, https://doi.org/10.1029/2011GL047871, 2011. 

Peters, W., Jacobson, A. R., Sweeney, C., Andrews, A. E., Con-way, T. J., Masarie, K., Miller, J. B., Bruhwiler, L. M. P. Pétron, G., Hirsch, A. I., Worthy, D. E. J., van der Werf, G. R.,Randerson, J. T., Wennberg, P. O., Krol, M. C., and Tans, P. P.: An atmospheric perspective on North American carbon dioxide exchange: CarbonTracker, P. Natl. Acad. Sci. USA, 104, 18925–18930, https://doi.org/10.1073/pnas.0708986104, 2007. 

Reuter, M., Buchwitz, M., Schneising, O., Heymann, J., Bovensmann, H., and Burrows, J. P.: A method for improved SCIAMACHY CO2 retrieval in the presence of optically thin clouds, Atmos. Meas. Tech., 3, 209–232, https://doi.org/10.5194/amt-3-209-2010, 2010. 

Reuter, M., Bovensmann, H., Buchwitz, M., Burrows, J. P., Connor, B. J., Deutscher, N. M., Griffith, D. W. T., Heymann, J., Keppel-Aleks, G., Messerschmidt, J., Notholt, J., Petri, C., Robinson, J., Schneising, O., Sherlock, V., Velazco, V., Warneke, W., Wennberg, P. O., and Wunch, D.: Retrieval of atmospheric CO2 with enhanced accuracy and precision from SCIAMACHY: Validation with FTS measurements and comparison with model results, J. Geophys. Res., 116, D04301, https://doi.org/10.1029/2010JD015047, 2011. 

Reuter, M., Buchwitz, M., Schneising, O., Hase, F., Heymann, J., Guerlet, S., Cogan, A. J., Bovensmann, H., and Burrows, J. P.: A simple empirical model estimating atmospheric CO2 background concentrations, Atmos. Meas. Tech., 5, 1349–1357, https://doi.org/10.5194/amt-5-1349-2012, 2012. 

Reuter, M., Bösch, H., Bovensmann, H., Bril, A., Buchwitz, M., Butz, A., Burrows, J. P., O'Dell, C. W., Guerlet, S., Hasekamp, O., Heymann, J., Kikuchi, N., Oshchepkov, S., Parker, R., Pfeifer, S., Schneising, O., Yokota, T., and Yoshida, Y.: A joint effort to deliver satellite retrieved atmospheric CO2 concentrations for surface flux inversions: the ensemble median algorithm EMMA, Atmos. Chem. Phys., 13, 1771–1780, https://doi.org/10.5194/acp-13-1771-2013, 2013. 

Reuter, M., Buchwitz, M., Hilker, M., Heymann, J., Schneising, O., Pillai, D., Bovensmann, H., Burrows, J. P., Bösch, H., Parker, R., Butz, A., Hasekamp, O., O'Dell, C. W., Yoshida, Y., Gerbig, C., Nehrkorn, T., Deutscher, N. M., Warneke, T., Notholt, J., Hase, F., Kivi, R., Sussmann, R., Machida, T., Matsueda, H., and Sawa, Y.: Satellite-inferred European carbon sink larger than expected, Atmos. Chem. Phys., 14, 13739–13753, https://doi.org/10.5194/acp-14-13739-2014, 2014a. 

Reuter, M., Buchwitz, M., Hilboll, A., Richter, A., Schneising, O., Hilker, M., Heymann, J., Bovensmann, H., and Burrows, J. P.: Decreasing emissions of NOx relative to CO2 in East Asia inferred from satellite observations, Nat. Geosci., 7, 792–795, https://doi.org/10.1038/ngeo2257, 2014b. 

Reuter, M., Buchwitz, M., Schneising, O., Noël, S., Rozanov, V., Bovensmann, H., and Burrows, J. P.: A Fast Atmospheric Trace Gas Retrieval for Hyperspectral Instruments Approximating Multiple Scattering – Part 1: Radiative Transfer and a Potential OCO-2 XCO2 Retrieval Setup, Remote Sens., 9, 1159, https://doi.org/10.3390/rs9111159, 2017a. 

Reuter, M., Buchwitz, M., Schneising, O., Noël, S., Bovensmann, H., and Burrows, J. P.: A Fast Atmospheric Trace Gas Retrieval for Hyperspectral Instruments Approximating Multiple Scattering – Part 2: Application to XCO2 Retrievals from OCO-2, Remote Sens., 9, 1102, https://doi.org/10.3390/rs9111102, 2017b. 

Reuter, M., Buchwitz, M., Schneising, O., Krautwurst, S., O'Dell, C. W., Richter, A., Bovensmann, H., and Burrows, J. P.: Towards monitoring localized CO2 emissions from space: co-located regional CO2 and NO2 enhancements observed by the OCO-2 and S5P satellites, Atmos. Chem. Phys., 19, 9371–9383, https://doi.org/10.5194/acp-19-9371-2019, 2019a. 

Reuter, M., Buchwitz, M., and Schneising-Weigel, O.: Algorithm Theoretical Basis Document (ATBD) – ANNEX D for products XCO2_EMMA, XCH4_EMMA, XCO2_OBS4MIPS, XCH4_OBS4MIPS (v4.1, 2003–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 37 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/ATBD/C3S_D312b_Lot2.1.3.2-v1.0_ATBD-GHG_ANNEX-D_v3.1.pdf (last access: 10 January 2020), 2019b. 

Reuter, M., Buchwitz, M., and Schneising-Weigel, O.: Product Quality Assessment Report (PQAR) – ANNEX D for products XCO2_EMMA, XCH4_EMMA, XCO2_OBS4MIPS, XCH4_OBS4MIPS (v4.1, 2003–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 37 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PQAR/C3S_D312b_Lot2.2.3.2-v1.0_PQAR-GHG_ANNEX-D_v3.1.pdf (last access: 10 January 2020), 2019c. 

Reuter, M., Buchwitz, M., and Schneising-Weigel, O.: Product User Guide and Specification (PUGS) – ANNEX D for products XCO2_EMMA, XCH4_EMMA, XCO2_OBS4MIPS, XCH4_OBS4MIPS (v4.1, 2003–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 22 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PUGS/C3S_D312b_Lot2.3.2.3-v1.0_PUGS-GHG_ANNEX-D_v3.1.pdf (last access: 10 January 2020), 2019d. 

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practice, World Scientific Publishing, 2000. 

Schaefer, H.: On the Causes and Consequences of Recent Trends in Atmospheric Methane, Curr. Clim. Change Rep., 5, 259–274, https://doi.org/10.1007/s40641-019-00140-z, 2019. 

Schneising, O., Buchwitz, M., Reuter, M., Heymann, J., Bovensmann, H., and Burrows, J. P.: Long-term analysis of carbon dioxide and methane column-averaged mole fractions retrieved from SCIAMACHY, Atmos. Chem. Phys., 11, 2863–2880, https://doi.org/10.5194/acp-11-2863-2011, 2011. 

Schneising, O., Reuter, M., Buchwitz, M., Heymann, J., Bovensmann, H., and Burrows, J. P.: Terrestrial carbon sink observed from space: variation of growth rates and seasonal cycle amplitudes in response to interannual surface temperature variability, Atmos. Chem. Phys., 14, 133–141, https://doi.org/10.5194/acp-14-133-2014, 2014a. 

Schneising, O., Burrows, J. P., Dickerson, R. R., Buchwitz, M., Reuter, M., and Bovensmann, H.: Remote sensing of fugitive methane emissions from oil and gas production in North American tight geologic formations, Earth's Future, 2, https://doi.org/10.1002/2014EF000265, 2014b. 

Schneising, O., Buchwitz, M., Reuter, M., Bovensmann, H., Burrows, J. P., Borsdorff, T., Deutscher, N. M., Feist, D. G., Griffith, D. W. T., Hase, F., Hermans, C., Iraci, L. T., Kivi, R., Landgraf, J., Morino, I., Notholt, J., Petri, C., Pollard, D. F., Roche, S., Shiomi, K., Strong, K., Sussmann, R., Velazco, V. A., Warneke, T., and Wunch, D.: A scientific algorithm to simultaneously retrieve carbon monoxide and methane from TROPOMI onboard Sentinel-5 Precursor, Atmos. Meas. Tech., 12, 6771–6802, https://doi.org/10.5194/amt-12-6771-2019, 2019. 

Sheng, J.-X., Jacob, D. J., Turner, A. J., Maasakkers, J. D., Benmergui, J., Bloom, A. A., Arndt, C., Gautam, R., Zavala-Araiza, D., Boesch, H., and Parker, R. J.: 2010–2016 methane trends over Canada, the United States, and Mexico observed by the GOSAT satellite: contributions from different source sectors, Atmos. Chem. Phys., 18, 12257–12267, https://doi.org/10.5194/acp-18-12257-2018, 2018. 

Sherlock, V., Connor, B., Robinson, J., Shiona, H., Smale, D., and Pollard, D. F.: TCCON data from Lauder (NZ), 125HR, Release GGG2014R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.lauder02.R0/1149298, 2014. 

Shiomi, K., Kawakami, S., Ohyama, H., Arai, K., Okumura, H., Taura, C., Fukamachi, T., and Sakashita, M.: TCCON data from Saga (JP), Release GGG2014R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.saga01.R0/1149283, 2014. 

Strong, K., Roche, S., Franklin, J. E., Mendonca, J., Lutsch, E., Weaver, D., Fogal, P. F., Drummond, J. R., Batchelor, R., and Lindenmaier, R.: TCCON data from Eureka (CA), Release GGG2014R3. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.eureka01.R3, 2019. 

Sussmann, R. and Rettinger, M.: TCCON data from Garmisch (DE), Release GGG2014R2. TCCON data archive, hosted by Caltech-DATA, https://doi.org/10.14291/tccon.ggg2014.garmisch01.R2, 2018. 

Té, Y., Jeseck, P., and Janssen, C.: TCCON data from Paris (FR), Release GGG2014R0. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.paris01.R0/1149279, 2014. 

Turner, A. J., Jacob, D. J., Wecht, K. J., Maasakkers, J. D., Lundgren, E., Andrews, A. E., Biraud, S. C., Boesch, H., Bowman, K. W., Deutscher, N. M., Dubey, M. K., Griffith, D. W. T., Hase, F., Kuze, A., Notholt, J., Ohyama, H., Parker, R., Payne, V. H., Sussmann, R., Sweeney, C., Velazco, V. A., Warneke, T., Wennberg, P. O., and Wunch, D.: Estimating global and North American methane emissions with high spatial resolution using GOSAT satellite data, Atmos. Chem. Phys., 15, 7049–7069, https://doi.org/10.5194/acp-15-7049-2015, 2015. 

Turner, A. J., Frankenberg, C., and Kort, E. A.: Interpreting contemporary trends in atmospheric methane, P. Natl. Acad. Sci. USA, 116, 2805–2813, https://doi.org/10.1073/pnas.1814297116, 2019. 

Veefkind, J. P., Aben, I., McMullan, K., Förster, H., De Vries, J., Otter, G., Claas, J., Eskes, H. J., De Haan, J. F., Kleipool, Q., Van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P.,Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P. F.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83, 2012. 

Velazco V. A., Morino, I., Uchino, O., Hori, A., Kiel, M., Bukosa, B., Deutscher, N. M., Sakai, T., Nagai, T., Bagtasa, G., Izumi, T., Yoshida, Y., and Griffith, D. W. T.: TCCON Philippines: First Measurement Results, Satellite Data and Model Comparisons in Southeast Asia, Remote Sens., 9, 1228, https://doi.org/10.3390/rs9121228, 2017. 

Warneke, T., Messerschmidt, J., Notholt, J.,Weinzierl, 5 C., Deutscher, N. M., Petri, C., Grupe, P., Vuillemin, C., Truong, F., Schmidt, M., Ramonet, M., and Parmentier, E.: TCCON data from Orléans (FR), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.orleans01.R1, 2019. 

Wennberg, P. O., Wunch, D., Roehl, C. M., Blavier, J.-F., Toon, G. C., and Allen, N. T.: TCCON data from Caltech (US), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.pasadena01.R1/1182415, 2015. 

Wennberg, P. O., Wunch, D., Roehl, C. M., Blavier, J.-F., Toon, G. C., Allen, N. T., Dowell, P., Teske, K., Martin, C., and Martin, J.: TCCON data from Lamont (US), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.lamont01.R1/1255070, 2016. 

Wennberg, P. O., Roehl, C. M., Wunch, D., Toon, G. C., Blavier, J.-F., Washenfelder, R., Keppel-Aleks, G., Allen, N. T., and Ayers, J.: TCCON data from Park Falls (US), Release GGG2014R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.parkfalls01.R1, 2017. 

Worden, J. R., Bloom, A. A., Pandey, S., Jiang, Z., Worden, H. M., Walter, T. W., Houweling, S., and Röckmann, T.: Reduced biomass burning emissions reconcile conflicting estimates of the post-2006 atmospheric methane budget, Nat. Commun., 8, 2227, https://doi.org/10.1038/s41467-017-02246-0, 2017. 

Wu, L., Hasekamp, O., Hu, H., Landgraf, J., Butz, A., aan de Brugh, J., Aben, I., Pollard, D. F., Griffith, D. W. T., Feist, D. G., Koshelev, D., Hase, F., Toon, G. C., Ohyama, H., Morino, I., Notholt, J., Shiomi, K., Iraci, L., Schneider, M., de Mazière, M., Sussmann, R., Kivi, R., Warneke, T., Goo, T.-Y., and Té, Y.: Carbon dioxide retrieval from OCO-2 satellite observations using the RemoTeC algorithm and validation with TCCON measurements, Atmos. Meas. Tech., 11, 3111–3130, https://doi.org/10.5194/amt-11-3111-2018, 2018. 

Wu, L., Aben, I., Hasekamp, O. P., and Buchwitz, M.: Product Quality Assessment Report (PQAR) – ANNEX B for products CO2_GOS_ SRFP, CH4_GOS_SRFP (v2.3.8, 2009–2018), Technical Report Copernicus Climate Change Service (C3S), version 3.1, 03-11-2019, 36 pp., available at: http://www.iup.uni-bremen.de/carbon_ghg/docs/C3S/CDR3_2003-2018/PQAR/C3S_D312b_Lot2.2.3.2-v1.0_PQAR-GHG_ANNEX-B_v3.1.pdf (last access: 10 January 2020), 2019.  

Wunch, D., Toon, G. C., Wennberg, P. O., Wofsy, S. C., Stephens, B. B., Fischer, M. L., Uchino, O., Abshire, J. B., Bernath, P., Biraud, S. C., Blavier, J.-F. L., Boone, C., Bowman, K. P., Browell, E. V., Campos, T., Connor, B. J., Daube, B. C., Deutscher, N. M., Diao, M., Elkins, J. W., Gerbig, C., Gottlieb, E., Griffith, D. W. T., Hurst, D. F., Jiménez, R., Keppel-Aleks, G., Kort, E. A., Macatangay, R., Machida, T., Matsueda, H., Moore, F., Morino, I., Park, S., Robinson, J., Roehl, C. M., Sawa, Y., Sherlock, V., Sweeney, C., Tanaka, T., and Zondlo, M. A.: Calibration of the Total Carbon Column Observing Network using aircraft profile data, Atmos. Meas. Tech., 3, 1351–1362, https://doi.org/10.5194/amt-3-1351-2010, 2010. 

Wunch, D., Toon, G. C., Blavier, J.-F. L., Washenfelder, R. A., Notholt, J., Connor, B. J., Griffith, D. W. T., Sherlock, V., and Wennberg, P. O.: The Total Carbon Column Observing Network, Philos. T. R. Soc. A, 369, 2087–2112, https://doi.org/10.1098/rsta.2010.0240, 2011. 

Wunch, D., Mendonca, J., Colebatch, O., Allen, N. T., Blavier, 5 J.-F., Roche, S., Hedelius, J., Neufeld, G., Springett, S., Worthy, D., Kessler, R., and Strong, K.: TCCON data from East Trout Lake, SK (CA), Release GGG2014.R1. TCCON data archive, hosted by CaltechDATA, https://doi.org/10.14291/tccon.ggg2014.easttroutlake01.R1, 2018. 

Yin, Y., Ciais, P., Chevallier, F., Li, W., Bastos, A., Piao, S., Wang, T., and Liu, H.: Changes in the response of the Northern Hemisphere carbon uptake to temperature over the last three decades, Geophys. Res. Lett., 45, 4371–4380, https://doi.org/10.1029/2018GL077316, 2018. 

Yoshida, Y., Kikuchi, N., Morino, I., Uchino, O., Oshchepkov, S., Bril, A., Saeki, T., Schutgens, N., Toon, G. C., Wunch, D., Roehl, C. M., Wennberg, P. O., Griffith, D. W. T., Deutscher, N. M., Warneke, T., Notholt, J., Robinson, J., Sherlock, V., Connor, B., Rettinger, M., Sussmann, R., Ahonen, P., Heikkinen, P., Kyrö, E., Mendonca, J., Strong, K., Hase, F., Dohe, S., and Yokota, T.: Improvement of the retrieval algorithm for GOSAT SWIR XCO2 and XCH4 and their validation using TCCON data, Atmos. Meas. Tech., 6, 1533–1547, https://doi.org/10.5194/amt-6-1533-2013, 2013. 

Download
Short summary
We present new satellite-derived data sets of atmospheric carbon dioxide (CO2) and methane (CH4). The data products are column-averaged dry-air mole fractions of CO2 and CH4, denoted XCO2 and XCH4. The products cover the years 2003–2018 and are merged Level 2 (satellite footprints) and merged Level 3 (gridded at monthly time and 5° x 5° spatial resolution) products obtained from combining several individual sensor products. We present the merging algorithms and product validation results.