Articles | Volume 12, issue 4
Data description paper
14 Dec 2020
Data description paper |  | 14 Dec 2020

A decade of GOSAT Proxy satellite CH4 observations

Robert J. Parker, Alex Webb, Hartmut Boesch, Peter Somkuti, Rocio Barrio Guillo, Antonio Di Noia, Nikoleta Kalaitzi, Jasdeep S. Anand, Peter Bergamaschi, Frederic Chevallier, Paul I. Palmer, Liang Feng, Nicholas M. Deutscher, Dietrich G. Feist, David W. T. Griffith, Frank Hase, Rigel Kivi, Isamu Morino, Justus Notholt, Young-Suk Oh, Hirofumi Ohyama, Christof Petri, David F. Pollard, Coleen Roehl, Mahesh K. Sha, Kei Shiomi, Kimberly Strong, Ralf Sussmann, Yao Té, Voltaire A. Velazco, Thorsten Warneke, Paul O. Wennberg, and Debra Wunch

This work presents the latest release (v9.0) of the University of Leicester GOSAT Proxy XCH4 dataset. Since the launch of the GOSAT satellite in 2009, these data have been produced by the UK National Centre for Earth Observation (NCEO) as part of the ESA Greenhouse Gas Climate Change Initiative (GHG-CCI) and Copernicus Climate Change Services (C3S) projects. With now over a decade of observations, we outline the many scientific studies achieved using past versions of these data in order to highlight how this latest version may be used in the future.

We describe in detail how the data are generated, providing information and statistics for the entire processing chain from the L1B spectral data through to the final quality-filtered column-averaged dry-air mole fraction (XCH4) data. We show that out of the 19.5 million observations made between April 2009 and December 2019, we determine that 7.3 million of these are sufficiently cloud-free (37.6 %) to process further and ultimately obtain 4.6 million (23.5 %) high-quality XCH4 observations. We separate these totals by observation mode (land and ocean sun glint) and by month, to provide data users with the expected data coverage, including highlighting periods with reduced observations due to instrumental issues.

We perform extensive validation of the data against the Total Carbon Column Observing Network (TCCON), comparing to ground-based observations at 22 locations worldwide. We find excellent agreement with TCCON, with an overall correlation coefficient of 0.92 for the 88 345 co-located measurements. The single-measurement precision is found to be 13.72 ppb, and an overall global bias of 9.06 ppb is determined and removed from the Proxy XCH4 data. Additionally, we validate the separate components of the Proxy (namely the modelled XCO2 and the XCH4∕XCO2 ratio) and find these to be in excellent agreement with TCCON.

In order to show the utility of the data for future studies, we compare against simulated XCH4 from the TM5 model. We find a high degree of consistency between the model and observations throughout both space and time. When focusing on specific regions, we find average differences ranging from just 3.9 to 15.4 ppb. We find the phase and magnitude of the seasonal cycle to be in excellent agreement, with an average correlation coefficient of 0.93 and a mean seasonal cycle amplitude difference across all regions of −0.84 ppb.

These data are available at (Parker and Boesch2020).

1 Introduction

Atmospheric methane (CH4) is the second most important greenhouse gas in terms of anthropogenic climate radiative forcing (Myhre et al.2013) with a global warming potential on a 100-year timescale of 28–34 times that of CO2 (Etminan et al.2016) on a mass-per-mass basis. This strong warming potential, when coupled to its short lifetime relative to that of CO2 (Prather et al.2012) makes it of particular interest when considering rapid and achievable mitigation strategies (Nisbet et al.2020).

Scientific debate continues on trying to explain the atmospheric CH4 trend observed over the past couple of decades. Records from surface sites reveal a plateau from 2000 to 2007 and a resumed increase after 2007 (Rigby et al.2008; Dlugokencky et al.2009). Amongst the varied surface sources of CH4, the largest are natural wetlands, agriculture, livestock, biomass burning, waste and fossil fuel production, whereas the primary sink is the OH radical in the atmosphere (Kirschke et al.2013; Saunois et al.2020). Various hypotheses have been offered that attempt to attribute the behaviour in the global growth rate to a particular component or mechanism (Nisbet et al.2016; Schaefer et al.2016; Hausmann et al.2016; McNorton et al.2016b; Buchwitz et al.2017a; Worden et al.2017; Turner et al.2017; Rigby et al.2017) but currently there is no consensus within the community. Many of these studies have utilised satellite observations of atmospheric CH4 and have shown the increasing capability of such measurements to characterise global and regional surface methane fluxes (Jacob et al.2016).

This work presents the most recent update to the University of Leicester GOSAT Proxy XCH4 retrieval. This version (v9.0) now provides over a decade of global total column CH4 observations, from April 2009 to December 2019. A full reprocessing of the entire time series has been performed to ensure consistency throughout the record and to ensure that results utilising the entire record are as robust as possible.

It is the intention of the authors that this study acts as a reference for everyone making use of the data, and as such, we have attempted to provide as much practical detail as possible on the usage of the data.

Section 2 describes the GOSAT observations themselves and highlights any instrument anomalies or data gaps. Section 3 is broken down into several sub-sections detailing the usage of previous versions of these data by the scientific community. Section 4 gives an overview of the retrieval method and details the end-to-end data processing chain including statistics on throughput and data availability. Section 5 shows the validation of the data against the Total Carbon Column Observing Network, characterising not only the final Proxy XCH4 data but also the individual components of the retrieval. Section 6 provides details of the global distribution of the data. Section 7 further characterises the data by performing model comparisons at global and regional scales. Finally we provide a summary and recommendations for future use in Sect. 8.

2 GOSAT TANSO-FTS observations

GOSAT was launched in January 2009 by the Japanese Space Agency (JAXA), as the first satellite mission dedicated to making greenhouse gas observations (Kuze et al.2009). GOSAT is nicknamed “Ibuki”, meaning “breath” in Japanese, highlighting that its mission involves monitoring the breathing of the planet, through measurement of the carbon cycle. In order to achieve this, GOSAT is equipped with a high-resolution Fourier transform spectrometer (TANSO-FTS – Thermal and Near Infrared Sensor for Carbon Observations – Fourier Transform Spectrometer). Shortwave infrared bands at 0.76 µm (O2), 1.6 µm (CO2 and CH4) and 2.0 µm (CO2) all provide near-surface sensitivity while a thermal infrared band between 5.5 and 14.3 µm provides mid-tropospheric sensitivity.

The objective for GOSAT was to provide routine measurements appropriate for regional, and continental-scale flux estimates. Kuze et al. (2009) state a target relative accuracy of 2 % for CH4 over 3-month averages at 100–1000 km spatial scales. The ESA GHG-CCI User Requirements Document (URD) specifies goal (G), breakthrough (B) and threshold (T) requirements, with the goal requirements being the most stringent (Buchwitz et al.2017a). For XCH4 these values are 9, 17 and 34 ppb respectively for precision and 1, 5 and 10 ppb for relative accuracy. We discuss in Sect. 5 how we exceed these breakthrough requirements.

As an FTS acquisition is relatively slow ( 4 s), the GOSAT sampling strategy is tailored to achieve this goal by measuring with a relatively large footprint of 10.5 km, spaced approximately  263 km apart across-track and  283 km along-track. This means that while GOSAT does not “image” the surface, it does return to the same location every 3 d, allowing a long time series of comparable measurements to be obtained. As well as nominally measuring over land in nadir mode, GOSAT is also capable of measuring over the ocean, which is normally too dark in the shortwave infrared (SWIR). This is achieved in the so-called “ocean sun-glint” observation mode, when the sun–satellite angle allows for a sufficiently reflected signal from the glint spot. Discussion of the GOSAT measurement sampling strategy and its evolution over time can be found in Appendix E.

2.1 Instrument anomalies and data gaps

Throughout its 10 years of operation, GOSAT has experienced a number of incidents resulting in instrument anomalies (Kuze et al.2016). These incidents include

  • May 2014 – a solar paddle incident resulting in a temporary instrument shutdown,

  • January 2015 – a switch to the secondary pointing mechanism due to degradation of the primary system,

  • August 2015 – a cryocooler shutdown and restart,

  • May 2018 – a CDMS (command and data management system) incident resulting in GOSAT being inactive for 2 weeks,

  • November 2018 – rotation anomaly of the second solar paddle.

The temporary reduction in observations related to these incidents is discussed in Sect. 4.3 and reflected in Fig. 4.

3 Studies utilising Proxy XCH4 data

The University of Leicester GOSAT Proxy XCH4 data are produced operationally for the ESA Greenhouse Gas Climate Change Initiative (Buchwitz et al.2017b) and the Copernicus Climate Change Service (C3S) (Buchwitz et al.2018) as well as routinely for the UK National Centre for Earth Observation. This work details version 9.0 of the University of Leicester GOSAT Proxy XCH4 data, but previous versions over the past decade have been used for a wide variety of scientific studies. This section details some of these past studies in order to highlight the potential applications for these data.

3.1 Validation of data

Firstly, before any conclusions can be drawn from analysis of the data, the data themselves must be validated to ensure their robustness and reliability. Previous versions of the data have been extensively validated against the TCCON network (Total Carbon Column Observing Network) as part of the ESA Climate Change Initiative (Parker et al.2011; Dils et al.2014; Buchwitz et al.2017a), including extensive validation of the model XCO2 used in the generation of the data (Parker et al.2015). We have also performed validation of the data against aircraft profile observations over the Amazon (Webb et al.2016), one of the most important and challenging regions for the retrieval.

3.2 Comparison to other satellite observations

Although GOSAT was the first satellite mission dedicated to measuring greenhouse gases, successful CH4 retrievals were performed previously from SCIAMACHY and continue to be performed from new missions such as TROPOMI on board Sentinel-5 Precursor and the recently launched GOSAT-2. Furthermore, many thermal infrared missions are capable of measuring CH4 (IASI, AIRS, TES, CrIS), albeit with sensitivity to the mid-troposphere and little sensitivity to the surface. Nevertheless, it is important that these different observations are consistent and their capabilities well-understood if we wish to perform long-term analysis. The ESA Greenhouse Gas Climate Change Initiative (ESA GHG-CCI) (Buchwitz et al.2017a) made substantial efforts to characterise and validate these different observations (Dils et al.2014). The ensemble median algorithm (EMMA) (Reuter et al.2020) homogenises the SCIAMACHY and GOSAT datasets produced via the ESA GHG-CCI project and is intended to be a long time series dataset for climate applications.

Studies such as those of Cressot et al. (2014) and Alexe et al. (2015) have investigated the consistency between flux inversions utilising SCIAMACHY and GOSAT (and IASI) CH4 observations and generally found good consistency in derived emissions. Worden et al. (2015) combined the surface sensitivity of GOSAT with the mid-tropospheric sensitivity of the NASA TES instrument to better estimate the lower tropospheric methane (and hence surface) emissions while Siddans et al. (2017) have compared their IASI CH4 product to GOSAT observations, finding good consistency between the two.

3.3 Investigation of the global growth rate

Perhaps the most important scientific question related to atmospheric CH4 concentrations is understanding the observed long-term behaviour. The cause of the so-called “hiatus” or plateau in atmospheric CH4 between 2000 and 2007 remains unresolved, with various studies speculating on the reason. Although the GOSAT record unfortunately only began in 2009, after the end of the plateau period, it can still help to characterise behaviour and understand the processes that may have contributed to the stalling. GOSAT data have been successfully used to infer long-term global fluxes. As well as contributing to the Global Methane Budget assessments (Saunois et al.2016, 2020), GOSAT data have been used to assess the role of regional wetland emissions (McNorton et al.2016b; Maasakkers et al.2019) and the role of OH variability as a potential cause for the stalling in growth rate (McNorton et al.2016a; Maasakkers et al.2019).

3.4 Regional emissions

GOSAT data have been successfully utilised in regional-scale studies to determine CH4 fluxes over many different regions. These types of studies are of particular interest as they can help inform policy-related discussions on validation and verification of regional- or country-scale emission targets, such as those relevant to the Paris Agreement (Bergamaschi et al.2018a). Fraser et al. (2013) performed regional flux inversions and found large changes over temperate Eurasia and tropical Asia, with the satellite observations providing a significant error reduction over only using surface data. Wecht et al. (2014) performed continental-scale inversions over North America and produced estimates of Californian CH4 emissions and found consistent emission estimates over the Los Angeles Basin between the satellite inversion and that from a dedicated aircraft campaign. Turner et al. (2015) extended this work to the entire US and inferred a US anthropogenic CH4 source over 50 % larger than that from EDGAR and EPA bottom-up inventories. Satellite inversion results from Alexe et al. (2015) showed a redistribution of CH4 emissions in the US from the north-east to south-central. These results are consistent with recent independent studies that suggest that bottom-up estimates of North American fossil fuel emissions (particularly related to natural gas and petroleum production facilities) are systematically underestimated. Ganesan et al. (2017) used GOSAT data to infer India's CH4 emissions between 2010–2015 and found average emissions of 22.0 Tg yr−1 to be consistent with the emissions reported by India to the UNFCCC with no significant trend over time. Sheng et al. (2018) performed a similar study over the US, Canada and Mexico and found that US emissions increased by 2.5 % over the 7-year study period and attributed this to contributions from the oil and gas industry and livestock. In Feng et al. (2017) the individual XCO2 and XCH4 components from the Proxy retrieval are used to infer regional CO2 and CH4 fluxes simultaneously. Finally, Lunt et al. (2019) inferred CH4 emissions over tropical Africa and found a linear increase of between 1.5 and 2.1 Tg yr−1 for 2010–2016, attributing much of this to short-term increase in emissions over the Sudd wetland area in South Sudan.

4 UoL Proxy XCH4 retrieval

The University of Leicester Full-Physics retrieval algorithm (UoL-FP) is based on the original NASA Orbiting Carbon Observatory (OCO) “Full-Physics” retrieval algorithm (Connor et al.2008; Boesch et al.2011; O'Dell et al.2012), which was designed to simultaneously fit the short-wave infrared radiances in the 0.76 µm O2 A-band and the 1.6 µm and 2 µm CO2 bands. This algorithm has been adapted for use on GOSAT observations and modified to perform a variety of different retrievals, including the Proxy method described here. The radiative transfer calculations are accelerated using the low-stream interpolation approach (O'Dell2010).

The concept behind the Proxy XCH4 retrieval approach (Frankenberg et al.2006) is that the majority of atmospheric scattering and instrument effects will be similar for CH4 and CO2 mole fraction retrievals performed in a common absorption band (around 1.6 µm, where both CO2 and CH4 have absorption features). By taking the ratio of the retrieved XCH4∕XCO2, the CO2 acts as a “proxy” for the modifications to the light path induced by scattering (Butz et al.2010) and cancels out those in the CH4 retrieval. This means that moderate aerosol scattering does not adversely impact the retrieval, resulting in a higher number of high-quality XCH4 observations compared to the full-physics approaches where much stricter post-filtering is often required. This is especially useful over the tropics, where moderate aerosol or cirrus effects can limit the coverage of full-physics methods but affect the Proxy approach far less severely.

In order to convert the retrieved XCH4∕XCO2 ratio into a final XCH4 quantity, a model-based estimate of XCO2 is used, according to the equation

(1) XCH 4 (proxy) = XCH 4 XCO 2 × XCO 2 (model) ,

where the relative variability of CO2 in the atmosphere is known to be much lower than that of CH4. This leads to the primary disadvantage of this method, i.e. that the model-based estimates of XCO2 may introduce biases in the retrieved CH4. In an attempt to minimise such biases, CO2 dry-air mole fractions XCO2(model) used in Eq. (1) in the UoL Proxy retrieval scheme are obtained by taking the median of the estimates produced by three atmospheric chemistry transport models which have assimilated surface in situ data: GEOS-Chem (Feng et al.2011), NOAA CarbonTracker (Peters et al.2007) and CAMS (Chevallier et al.2010).

The advantage of the Proxy retrieval approach compared to the “Full-Physics” retrieval as typically used for CO2 (Boesch et al.2011; Cogan et al.2012) is that Proxy retrievals are less sensitive to instrumental effects and require less-strict quality filtering (Schepers et al.2012; Parker et al.2015), thereby ensuring a better coverage of regions (especially in the tropics), where full-physics retrievals are particularly challenging.

Figure 1The average CH4 a priori profile used in the retrievals along with the 25 %–75 % variation (a). The average normalised column averaging kernel produced by the retrieval along with the 25 %–75 % variation, highlighting the strong sensitivity of the retrieval to the surface (b).


4.1 Retrieval inputs and a priori generation

In order to prepare all of the necessary inputs to the retrieval, we use the Leicester Retrieval Preparation Toolset (LRPT) software. The latest version of the GOSAT Level 1B files (version 210.210) is acquired directly from the NIES GDAS Data Server and is processed with the LRPT to extract the measured radiances along with all required sounding-specific ancillary information such as the measurement time, location and geometry. These measured radiances have the recommended radiometric calibration and degradation corrections applied as per Yoshida et al. (2013) with an estimate of the spectral noise derived from the standard deviation of the out-of-band signal. We then format the spectral data for input into the UoL-FP retrieval algorithm and generate a list containing all of the ancillary data necessary to create the retrieval a priori information.

Figure 2Spectra showing the GOSAT radiances for the CO2 and CH4 bands and the resulting residuals (measured  simulated spectral differences). The data shown are the median values for the 25 274 land retrievals that pass the quality filtering for an example month, August 2016.


Sounding-specific a priori information is generated for all individual soundings present in the sounding-selector list above. Atmospheric temperature and water vapour profiles are taken from ECMWF ERA-Interim up to August 2019 and ERA-5 thereafter. CO2 profile information is taken from the 16r1 CAMS atmospheric inversion (Chevallier2019) and incremented by the NOAA estimated global growth rate for recent years. CH4 profiles are taken from a combination of the MACC-II CH4 inversion (v10-S1NOAA – for the troposphere and from a dedicated TOMCAT stratospheric chemistry model simulation (Chipperfield1999) for the stratospheric component. This ensures that the CH4 a priori profiles are sufficiently vertically resolved and capture the sharp decrease in concentration around the tropopause (Fig. 1, top). As the MACC-II data are only available until 2012, the data after this period are repeated each year.

Figure 3Sankey diagram detailing the retrieval throughput at each step of the processing chain for the GOSAT Proxy XCH4 retrieval. As well as the absolute number of soundings, the percentage relative to the initial total is also given.


All of the atmospheric profiles are then interpolated to a sounding-specific retrieval grid. For each sounding a 20-level pressure-based retrieval grid is generated that ranges from the top of the atmosphere (0.1 hPa) to 20 hPa beneath the surface pressure as estimated by ERA-Interim. This 20 hPa buffer allows the surface pressure to be adjusted during the initial cloud-screening process without leading to unphysical extrapolation of the a priori profiles.

In addition to atmospheric a priori information, we also generate sounding-specific a priori information for the spectral dispersion and surface albedo directly from the GOSAT spectra.

The retrieval algorithm also requires the input of spectroscopic parameters for the species being simulated. We use v4.2.0 of the OCO line lists for CO2, H2O and O2, and we take CH4 parameters from the TCCON line lists (Toon2015).

4.2 Cloud screening

Prior to the XCH4 and XCO2 retrievals, cloudy GOSAT soundings are identified and excluded by using the UoL-FP retrieval algorithm to obtain the apparent surface pressure from O2 A-band spectra and comparing it to the surface pressure provided by the ECMWF reanalysis (Dee et al.2011). If the absolute difference between the retrieved and the ECMWF surface pressure is larger than 30 hPa, a sounding is flagged as cloudy and excluded from further processing. The reason why a loose threshold is used for the surface pressure difference is that this procedure only aims to identify and remove soundings which are significantly cloudy. Partially cloudy scenes, or scenes where optically thin clouds are present, are processed by the retrieval algorithm and are dealt with through a post-retrieval quality filtering scheme described later in this section.

4.3XCH4 and XCO2 retrievals

For soundings that pass the cloud screening procedure described above, retrievals with the UoL-FP algorithm for CO2 and CH4 mole fraction profiles are carried out separately. The state vector for these retrievals consists of 20-level profiles for CH4 and CO2 mole fractions along with profile scaling factors for H2O mole fraction and temperature with parameters for surface albedo and spectral dispersion also included, allowing us to explicitly fit the wavelength for each spectra independently.

A post-retrieval quality filtering is then carried out, by selecting the retrievals that meet the following criteria: (1) goodness-of-fit (χ2) parameter between 0.4 and 1.9 for both CH4 and CO2, (2) a posteriori error smaller than 20 ppb for CH4 and 3 ppm for CO2, (3) retrieved XCH4 larger than 1650 ppb and XCO2 larger than 350 ppm, and (4) latitude north of 60 S (to exclude Antarctica).

Figure 2 shows an example of the spectral fits for 1 month of data (August 2016) for all 25 274 successful, quality-filtered data measured over land. The top panels show the averaged measured radiances in the 1.6 µm CO2 (left) and 1.65 µm CH4 (right) retrieval windows, with the bottom panels showing the residual differences to the final simulated spectra. The estimated instrument noise is indicated by the shaded area and the residuals are found to be within the noise.

The retrieved XCH4 and XCO2 satisfying the aforementioned quality criteria are then used in Eq. (1), together with the ensemble median model XCO2 described earlier in this section. Prior to the calculation of XCO2 to be used in Eq. (1), model CO2 profiles are convolved with scene-dependent instrument averaging kernels computed as part of the CO2 retrieval.

Before the final production of the data files, an offset is subtracted from the retrieved XCH4 to remove a residual mean bias to TCCON (see Sect, 5). Currently, a single offset value of 9.06 ppb is used. This offset is applied to all analyses presented here and is built in to the final delivered data.

A summary of the throughput of the whole processing chain described in this section is shown in Fig. 3. This shows that in total between April 2009 and December 2019 we have 19.5 million individual GOSAT soundings which are ingested into the LRPT software. Of these, 95.4 % are successfully preprocessed, with the 4.6 % that fail largely due to incomplete or invalid L1B data. A very small number of successfully preprocessed soundings < 0.1 % are unable to be processed further as we are unable to estimate a noise for those spectra. Of the 18.6 million soundings that continue and are attempted for cloud-clearing, 17.7 million are able to be successfully cloud-cleared. Of these, just over 7 million soundings are found to be cloud-free, with over 10 million determined to be cloudy. A successful CH4 retrieval is performed on the majority of these cloud-free soundings, with just 41 597 failing the retrieval. Of the 7.3 million successful CH4 retrievals, 2.7 million are rejected by our final quality filtering. It should be noted that currently we exclude all retrievals below 60 S (i.e. Antarctica) due to low signal-to-noise ratio and difficulty in distinguishing low cloud from the snow-covered surface. This alone accounts for over 1 million of the 2.7 million rejected retrievals. Finally, we are left with almost 4.6 million successful and quality-controlled XCH4 retrievals (23.5 % of the total measurements performed).

Figure 4The number of GOSAT soundings per month in the final Proxy XCH4 dataset, starting in April 2009 and ending in December 2019. Counts are provided for the global total but also for land and glint observation modes separately.


Figure 4 shows the number of successful retrievals broken down by month and also split into land and glint observation modes. This figure is particularly useful as it highlights any systematic differences in data density over time (e.g. from changes in the GOSAT sampling strategy) and also highlights abrupt data gaps (e.g. from instrument anomalies). In particular, it shows that the increase in monthly data from 2014/2015 onwards is largely a result of an increase in the number of glint observations which is a direct result of the instrument sampling changes that increased the valid glint range. This does highlight that some care should be taken when using the data for some applications as there cannot be assumed consistent temporal–spatial data coverage over the whole data record. Large data gaps, such as in January 2015 and December 2018, are also highlighted and indicate where care may need to be taken when analysing over these periods.

Figure 5Time series plots for each TCCON site comparing GOSAT Proxy XCH4 data to the matching co-located TCCON measurements. Statistics for each site are included, showing the average GOSAT–TCCON difference, the GOSAT–TCCON standard deviation, the correlation coefficients and the number of co-located GOSAT–TCCON measurements.


Figure 6Summary of the statistics from Fig. 5 comparing the GOSAT Proxy XCH4 data to the TCCON measurements for each TCCON site. Panels show the GOSAT–TCCON difference, the GOSAT–TCCON standard deviation, the correlation coefficients between the GOSAT and TCCON data and finally the number of co-located GOSAT–TCCON measurements.


5 Validation against TCCON

Evaluation against the Total Carbon Column Observing Network (TCCON) is the primary mechanism by which satellite-based measurements of XCO2 and XCH4 are validated.

The TCCON network consists of ground-based high-resolution Fourier transform spectrometers, performing direct measurements of solar spectra in the near-infrared. There are currently 27 operational sites located across North America, Europe, Asia and Oceania, including several islands in the Southern Hemisphere and other remote areas. TCCON sites have become operational at different times (see Table F1), and hence the data record length varies between sites, with Burgos (Philippines) (Velazco et al.2017) and Nicosia (Cyprus) being the most recent to come online in 2017 and 2019 respectively. Also note that the Lauder site (Pollard et al.2017) has multiple instruments and we have kept these records separate.

TCCON has been used extensively to validate satellite observations of XCH4 from SCIAMACHY (De Mazière et al.2004; Dils et al.2006), GOSAT (Dils et al.2014; Parker et al.2011; Yoshida et al.2013; Parker et al.2015; Inoue et al.2016) and TROPOMI (Hu et al.2018; Schneising et al.2019; Lambert et al.2019) and allows these measurements to be bias-corrected where necessary. TCCON itself is tied to the World Meteorological Organization (WMO) standard through comparison against integrated aircraft measurements (Wunch et al.2010, 2011).

This work uses the latest available TCCON data, GGG2014. Detailed dataset citations are available for each site in Table F1.

For comparison between TCCON and GOSAT, all GOSAT soundings within ±5 of a TCCON site are taken. For these soundings, the average of the TCCON data within ±2 h of the GOSAT overpass time is calculated, resulting in GOSAT–TCCON pairs when there are TCCON data available. It is these matched GOSAT–TCCON pairs that are then subsequently analysed. In total we use 22 of these sites within our analysis (see Table F1), omitting some sites with insufficient data coverage or high-altitude sites where the total column may not be well-represented or co-located well to the satellite observations. Figure 5 shows the time series of the GOSAT (blue) and TCCON (orange) data between 2009 and 2019 for each individual site. Also provided are the mean GOSAT–TCCON difference (Δ), the standard deviation of the GOSAT–TCCON difference (σ), the correlation coefficient (R) and the total number of GOSAT–TCCON pairs (N). These statistics are also summarised in Fig. 6, which shows that generally the GOSAT–TCCON difference is small (< 5 ppb), the standard deviation (which can be considered to be the single-measurement precision of the GOSAT data) is typically between 10 and 15 ppb, the correlation coefficient is generally high (0.7–0.9) and there are many co-located GOSAT–TCCON measurements with the distribution changing considerably between TCCON sites. Another important validation metric is the relative accuracy, or inter-station bias. This metric is an indication of any spatio-temporal variability of the bias and is defined in Dils et al. (2014) as the standard deviation of the individual site biases. We obtain a value of 3.89 ppb for this metric, again smaller than the estimated TCCON accuracy of ±4 ppb. This meets the “breakthrough” user requirement for the systematic error of 5 ppb as defined by Buchwitz et al. (2017a).

In total across all TCCON sites we find 88 345 matching GOSAT–TCCON data pairs. The correlation between the GOSAT and TCCON data is shown in Fig. 7, presented as a 2-D kernel density estimation (KDE) plot, along with the corresponding marginal 1-D KDE plots on the X and Y axes. An overall difference of 9.06 ppb is removed from the GOSAT data so that, by design, the absolute average difference to TCCON is 0 ppb. The overall standard deviation or single-measurement precision is found to be 13.72 ppb with a correlation coefficient of 0.92. The single-measurement precision of 13.72 ppb comfortably exceeds the precision breakthrough requirement of 17 ppb (Buchwitz et al.2017a), indicating that it “would result in a significant improvement for the targeted application”. Although the data contributing to this plot are from a wide variety of TCCON sites in different locations and at different latitudes, the distribution appears consistent and is tightly aligned to the one-to-one (dashed) line. However, there are signs of a potential hemispheric or latitudinal bias in the data against TCCON, although this is not apparent at all sites; for example Karlsruhe, Lamont, Tsukuba and Lauder all have negligible biases but span a large latitude range. It should also be noted that the uncertainty on the TCCON XCH4 is approximately 4 ppb, and for the majority of sites the GOSAT–TCCON difference is within this uncertainty so care must be taken to not over-interpret any signals at this scale.

Figure 7Correlation between the 88 345 matching TCCON XCH4 data and co-located GOSAT Proxy XCH4 measurements across all TCCON sites. The data are presented as a two-dimensional kernel density estimation (KDE) plot. The distribution sites along the one-to-one line (grey dashed) with a standard deviation (i.e. single-sounding precision) of 13.72 ppb and a correlation coefficient (R) of 0.92. An overall bias to TCCON of 9.06 ppb is removed from the GOSAT data, resulting in an average bias of 0 ppb by design. The individual KDE plots are shown along the upper and right margins.


5.1XCH4∕XCO2 ratio validation

In addition to validation of our final XCH4 Proxy dataset, TCCON data allow us the opportunity to validate the different components in Eq. (1), namely the XCH4∕XCO2 ratio and the model-derived XCO2.

Figure 8 shows the correlation between the retrieved GOSAT XCH4∕XCO2 ratio (ppb ppm−1) (with no bias correction applied) and the corresponding ratio calculated from TCCON. There is an excellent correlation coefficient of 0.89 across the 88 345 matching data points with a standard deviation of just 0.03 ppb ppm−1. An average offset between the two datasets of 0.02 ppb ppm−1 exists and is of a very similar magnitude to the global offset that is removed from the final data of 9.06 ppb. To be clear, the final bias correction which we apply to the Proxy XCH4 is almost entirely attributed to this bias that we identify here in the XCH4∕XCO2 ratio. It should also be noted here that the TCCON data themselves have a bias correction applied to the XCO2 and XCH4 data. This air-mass-independent correction factor derived from airborne calibrations is 1∕0.9898 for XCO2 and 1∕0.9765 for XCH4 (Wunch et al.2010 – Table 5). It is considered that this correction is mainly a result of deficiencies in the spectroscopy, which likely apply to the GOSAT retrievals as well and might go some way to explaining this small difference between TCCON and GOSAT.

Figure 8Correlation between the 88 345 matching TCCON XCH4∕XCO2 ratios and co-located GOSAT XCH4∕XCO2 ratios retrieved as a raw fundamental part of the Proxy XCH4 retrieval (Eq. 1). The data are presented as a two-dimensional kernel density estimation (KDE) plot. The distribution sites along the one-to-one line (grey dashed) with a standard deviation of 0.03 ppb ppm−1 and a correlation coefficient (R) of 0.89. An overall bias to TCCON of 0.02 ppb ppm−1 is present in these raw data. The individual KDE plots are shown along the upper and right margins.


5.2 Validation of XCO2 model

To validate the XCO2 model data used in the generation of the final Proxy data, we evaluate the model median XCO2 mixing ratios against TCCON but also evaluate the three individual models, sampled at the time and location of the GOSAT soundings, with the GOSAT sounding-specific averaging kernel applied. These XCO2 models are all independent of TCCON data but do assimilate NOAA surface site measurements, some of which are close to TCCON sites. Figure 9 shows the correlation between TCCON and the model median XCO2 (top left), GEOS-Chem XCO2 (top right), CAMS XCO2 (bottom left) and CarbonTracker XCO2 (bottom right). In all four cases there is an excellent agreement between TCCON and the model data with correlation coefficients all at 0.99, with average differences ranging from −0.07 to 0.17 ppm and standard deviations between 1.02 and 1.17 ppm. Although the values are all very similar, the model with the smallest difference (GEOS-Chem) has the largest standard deviation, and conversely the model with the largest difference (CAMS) has the smallest standard deviation. Overall the model median performs marginally better than any of the individual models with a standard deviation of 1.02 ppm but retaining a very small difference (0.08 ppm). Parker et al. (2015) provide a detailed assessment of the contribution to the overall uncertainty on the Proxy XCH4 related to the model XCO2.

Figure 9Correlation of the individual XCO2 models (GEOS-Chem, MACC and Carbon-Tracker) used in the generation of the GOSAT Proxy XCH4 data against matching co-located TCCON XCO2 measurements. The first panel shows the median model value which is the quantity directly used in Eq. (1) to generate the final Proxy XCH4 quantities. The data are presented as a two-dimensional kernel density estimation (KDE) plot with the one-to-one line shown as the grey dashed line. Also included are the statistics (difference, standard deviation, correlation coefficient and number of matching measurements) comparing each model to the TCCON XCO2.


Figure 10Global maps showing the GOSAT Proxy XCH4 data from April 2009 to December 2019 separated into seasons – spring (MAM), summer (JJA), autumn (SON) and winter (DJF).

6 Global CH4 distributions

As discussed in Sect. 3, one of the primary applications for the GOSAT Proxy XCH4 data has been as input to global flux inversions. For this reason, it is useful to examine the global spatio-temporal distribution of the data. Figure 10 shows seasonal maps of the GOSAT data from spring (March–April–May) 2009 through to winter (December) 2019. Features of note include a consistent increase in concentration over time; strong regional signals associated with CH4 surface sources, particularly over South America, India, China and Africa; a clear seasonal cycle over many regions; and a significant increase in the number and latitudinal range of GOSAT ocean sun-glint observations from 2014/2015.

Despite changes in the GOSAT sampling pattern and various instrument issues resulting in data gaps (see Sect. 2.1), on a seasonal and global scale there is good data coverage throughout the entire decade of observations.

Figure 11Hovmöller (latitude vs. time) plot of the GOSAT Proxy XCH4 (a), TM5 model XCH4 (b) and GOSAT–TM5 difference (c). The model data have been sampled at the time and location of the GOSAT observation and have had the sounding-specific GOSAT averaging kernel applied.


7 Model comparisons

The purpose of this paper is to present details of the v9.0 Proxy dataset and provide information to facilitate the future use of the data. As such, it is not the intention that this work performs detailed scientific analysis and interpretation. We do not, for example, perform any atmospheric flux inversions using these data as that is a significant study in its own right. However, we do feel that it would be informative to users of the data for us to perform a comparison against existing model XCH4 simulations to give confidence that the data are of sufficient quality to use in such studies.

Figure 12Time series comparison between the GOSAT Proxy XCH4 data and TM5 model XCH4 simulations (Bergamaschi et al.2018b) for individual TransCom regions. The model data have been sampled at the time and location of the GOSAT observation and had the sounding-specific GOSAT averaging kernel applied.


In this section we compare the GOSAT Proxy XCH4 data to a simulation of the TM5 global chemistry transport model (Bergamaschi et al.2013, 2018b) which has assimilated NOAA surface measurements. We have chosen to compare against model simulations that are both widely used within the community and that have already assimilated NOAA surface measurements. The reasoning for this is that any overall differences as might be seen from free-running model simulation are removed and we can clearly state the consistency of our dataset with the NOAA network. By proving good overall agreement to both TCCON measurements (Sect. 5) and the NOAA-constrained model simulations we, believe this indicates the utility of our data for use in further scientific analysis. These model data are the same data as used in the Global Methane Budget 2000–2017 (Saunois et al.2020) and in that study are referred to as “TM5-4DVAR”. For these comparisons, we have sampled the TM5 model at the time and location of the GOSAT measurement, interpolated the model to the GOSAT retrieval grid, applied the sounding-specific GOSAT CH4 averaging kernel (see Fig. 1 and Appendix D), and computed the model total column XCH4 amount.

Figure 11 shows Hovmöller plots (latitudinal mean versus time) for the GOSAT Proxy XCH4 (top), TM5 model simulation (middle) and GOSAT–TM5 difference (bottom). The GOSAT distribution behaves as expected, showing an increase over time between 2009 and 2019 superimposed on top of a north–south gradient and regular seasonal cycle. The TM5 data exhibit very similar characteristics and are in very good agreement to the GOSAT data. The difference between the two datasets (lower panel) shows that although there is a small offset between the two (with GOSAT on average 6.55 ppb larger than the model), there is very good consistency over time. GOSAT and TM5 seem to agree slightly better during the peak of the seasonal cycle, particularly in the tropics, with TM5 exhibiting a shallower trough. At very high northern latitudes, GOSAT is slightly lower than the model, but this relates to observations over Greenland at high altitude and low signal-to-noise ratio for the GOSAT soundings, so care must be taken not to over-interpret this difference.

With Fig. 11 providing confidence that the GOSAT data and TM5 simulation are in broad agreement, it is informative to break these comparisons down to a regional scale. Figure 12 shows time series of the GOSAT and TM5 data over the 16 different TransCom regions (Gurney et al.2002). All regions show good agreement between the modelled and observed data. We compute the de-seasonalised XCH4 over time for each region. The average difference in model and observation ranges from 3.9 ppb (Eurasia Boreal) to 15.4 ppb (Southern Tropical Asia). On average across all regions, the mean difference between the model and observation is 9.8 ppb.

The observed seasonal cycles in each region are very well-represented by the model, with an average correlation coefficient of 0.93 (ranging from 0.84 to 0.98 across all regions). The peak-to-peak seasonal cycle timing and magnitude are very well-reproduced between the two datasets. For example, the average peak-to-peak seasonal cycle amplitude for northern tropical Asia is 61.1 ppb for GOSAT and 62.3 ppb for TM5 whilst for North American Temperate it is 22.7 ppb for GOSAT and 23.6 ppb for TM5. The average GOSAT–model difference between the mean seasonal cycle amplitude across all regions is −0.84 ppb, with the average absolute difference being 5.3 ppb with only a small number of instances where the model and observation strongly disagree (for example in North America Boreal and South America Tropical in 2015, likely indicative of isolated regional emissions).

It is not the purpose of this paper to diagnose or interpret detailed differences between the observations and model, but it is useful to make a few observations relevant for use of the data within future studies. Figure 12 indicates a potential latitude-dependent bias, which is most likely due to model deficiencies in simulating the stratosphere (especially at mid-latitudes and high latitudes; Patra et al.2011; Alexe et al.2015; Saad et al.2016; Wang et al.2017) or inadequate inter-hemispheric mixing but could partly also indicate some latitudinal bias of the satellite retrievals. Accounting for such a latitudinal dependence through the fitting of a second-order polynomial function (as in Bergamaschi et al.2013; Turner et al.2015) may improve the baseline agreement between model and observation and is an approach that users may wish to explore depending upon their application. Furthermore, these model simulations are constrained by the NOAA background observations. Therefore differences between TM5-4DVAR and GOSAT may partly reflect deficiencies of bottom-up inventories used as prior, particularly over strong emission regions (e.g. obvious deficiencies in tropical Africa related to wetlands). When incorporating the GOSAT data into such inversions, this leads to the production of significant increments in the inverted fluxes and better agreement between observation and simulation (as in Alexe et al.2015, and other studies noted in Sect. 3.4).

8 Data availability

The University of Leicester GOSAT Proxy v9.0 XCH4 data are available from the Centre for Environmental Data Analysis data repository at (Parker and Boesch2020). The TCCON data are available from the TCCON Data Archive at (TCCON Data Archive2020; individual data citations are provided in Table F1). CAMS model CO2 (v18r2) data is available from the Copernicus Atmospheric Data Store at!/dataset/cams-global-greenhouse-gas-inversion (ECMWF2020a). MACC model CH4 (v10-S1NOAA) is available from ECMWF at (ECMWF2020b). NOAA CarbonTracker model CO2 (CT2017 and CT2019-NRT) are available from NOAA ESRL, Boulder, Colorado, USA, at (NOAA2020).

9 Summary and outlook

In this work we have presented the latest version of the University of Leicester GOSAT Proxy XCH4 dataset. This dataset now contains over a decade of global CH4 observations, sensitive to surface emissions and hence suited to estimating CH4 fluxes. The capability to estimate global and regional CH4 emissions is vital to improving our understanding of the global methane budget and how this budget may respond and change with respect to a changing future climate.

We begin this work by highlighting the wide variety of studies that previous versions of this dataset have contributed towards, demonstrating the significant utility of this dataset for examining and understanding the global methane budget.

This work provides a thorough description of the data processing chain, explaining in detail how the data are generated and how the high quality of the dataset is ensured. Extensive validation of the data against the TCCON network is performed, validating not only the final Proxy XCH4 data but also the separate components (the XCH4∕XCO2 ratio and the modelled XCO2) that form the final data product.

We also provide global seasonal maps of the data that demonstrate the spatial distribution of the data as well as highlighting particular features and regions that may be of interest for more detailed study.

Finally, as the primary usage of the data is expected to be as input into a flux inversion data assimilation framework in conjunction with atmospheric chemistry transport models and observations from surface networks, it is useful to compare the consistency against existing model simulations. We compare zonally and regionally against TM5 simulations that have assimilated observations from the NOAA surface network. We generally find a high level of consistency whilst identifying the additional utility that the satellite observations should introduce to the system.

Despite GOSAT-1 having a planned mission lifetime of 5 years, it continues to successfully perform measurements 11 years after launch. GOSAT-2 was launched in October 2018 (Suto et al.2020) and will continue the legacy of the GOSAT-1 mission. GOSAT-2 offers several opportunities for development related to the dataset we describe here. Primarily, it ensures that should GOSAT-1 cease operation, the valuable decade-long time series of observations can continue to be extended via GOSAT-2. With a significant overlap in time between the two missions, consistency between the two missions can be assured, albeit with significant future work/development.

In addition, GOSAT-2 has additional capabilities, namely the possibility of measuring carbon monoxide (CO). By measuring CO2, CH4 and CO simultaneously from the same instrument, GOSAT-2 would allow the extension of studies examining biomass burning combustion, leading to constraints on fire emission ratios as have been performed previously for GOSAT-1 (Ross et al.2013; Parker et al.2016).

A strong focus of future CH4-measuring satellites will be to examine anthropogenic emission sources at very high spatial resolution (e.g. PRISMA Pignatti et al.2013; HISUI Matsunaga et al.2017; ENMAP Guanter et al.2015), particularly relating to monitoring of the oil and gas industry. However, many scientific challenges and questions remain regarding the long-term CH4 behaviour and the response to a changing climate. For this reason, a long-term, consistent climate-ready data record as we present here is of continued importance. We expect that these data will be valuable for numerous studies, from regional flux inversions to monitoring long-term trends. With now over a decade of global atmospheric XCH4 observations, this dataset has helped, and will continue to help, us better understand the global methane budget and investigate how it may respond to a future changing climate.

Appendix A: Summary of dataset characteristics

Table A1 summarises the key characteristics of the University of Leicester GOSAT Proxy XCH4 data, including the spatial and temporal extent that the dataset covers, the total number of measurements, and their evaluation against TCCON.

Table A1Table summarising the key characteristics of the University of Leicester GOSAT Proxy XCH4 data.

Download Print Version | Download XLSX

Appendix B: Previous data versions

Table B1 outlines the history and evolution of the University of Leicester GOSAT Proxy XCH4 data product. Entries include the version number, the project that the data were generated for, the version of the GOSAT L1B data used, the time period covered by the data, whether ocean sun-glint data were generated, comments relating to changes/updates from previous versions and peer-reviewed publications that we are aware of that used the data. For the ESA GHG-CCI project, we also indicate which versions were officially delivered as part of the Climate Research Data Packages through the project. All Copernicus C3S versions were delivered to the Copernicus Climate Data Store.

Parker et al. (2011)Ross et al. (2013); Fraser et al. (2013); Dils et al. (2014); Wecht et al. (2014); Cressot et al. (2014); Berchet et al. (2015); Worden et al. (2015)Yoshida et al. (2013)Fraser et al. (2014); Turner et al. (2015); Alexe et al. (2015)Kuze et al. (2014)Parker et al. (2015); Stanevich et al. (2019, 2020)Webb et al. (2016); McNorton et al. (2016b); Siddans et al. (2017); Feng et al. (2017); Buchwitz et al. (2017b, 2018)Ganesan et al. (2017)Parker et al. (2016, 2018); McNorton et al. (2018); Sheng et al. (2018); Maasakkers et al. (2019, 2020); Lunt et al. (2019)Zheng et al. (2019); Reuter et al. (2020); Yin et al. (2020); Parker et al. (2020b); Saunois et al. (2020); Tunnicliffe et al. (2020)Parker et al. (2020a); Lu et al. (2020); Zhang et al. (2020)

Table B1Table showing the evolution of the University of Leicester GOSAT Proxy XCH4 data product. Entries include the version number, the project that the data were generated for, the version of the GOSAT L1B data used, the time period covered by the data, whether ocean sun-glint data were generated, comments relating to changes/updates from previous versions and peer-reviewed publications that we are aware of that used the data. For the ESA GHG-CCI project, we also indicate which versions were officially delivered as part of the Climate Research Data Packages through the project. All Copernicus C3S versions were delivered to the Copernicus Climate Data Store.

Download Print Version | Download XLSX

Appendix C: Data contents and usage notes

This section provides information on the contents and usage of the netCDF data files that we provide containing the Proxy XCH4 data. Whilst we recommend that anyone using the data should discuss their specific usage with the author, the following information is useful to note.

Our data are delivered as daily netCDF files, containing n individual GOSAT soundings. We provide everything in the data files that we believe users would require to make use of our data, including our a priori information (ch4_profile_apriori) and averaging kernels (xch4_averaging_kernel) which are provided on m vertical levels (see Appendix D).

In general, users should only use data that pass our quality checks (i.e. xch4_quality_flag == 0). In some specific use cases, the data that have failed our checks may still be of use, but additional care should be taken in using these data and we strongly recommend discussing such applications with us to determine if that is suitable for your use.

We provide data from both observation modes (nadir land and ocean sun glint) in the same file. While we do not believe that we have any bias between these different modes, for some use cases, users may wish to exclude either of these modes. retr_flag provides information on the mode for each sounding (0 = land, 1 = glint).

The variable named xch4 refers to the final Proxy XCH4 as calculated using Eq. (1). This is the main data product that we provide. In addition to this, we also provide the other components of Eq. (1). raw_xch4 and raw_xco2 refer to the directly retrieved XCH4 and XCO2 quantities. These variables should generally not be used but may be useful for certain applications. For example, some users may wish to use the XCH4∕XCO2 ratio (i.e. raw_xch4/raw_xco2) but replace the model XCO2 that we use (model_xco2) with their own modelled XCO2 which may be more appropriate for their particular application or more consistent with their own model transport. Our retrievals are typically performed on 20 vertical levels, with the first (bottom) level being the surface pressure. However, in a number of instances, especially over high terrain, where the apparent surface pressure from our O2 A-band cloud screening is above the bottom two levels of our pressure profile, this can result in only 19 active retrieval levels. These data are still valid, but variables with a vertical dimension (m) will contain a fill_value of -9999.99 at the first/lowest value. This value should be checked for, and that particular profile should be considered to only have 19 levels, rather than the standard 20.

We identify individual GOSAT soundings by their exposure_id. This may be of use when attempting to match our data to other GOSAT data products. This is a numerical identification that matches the GOSAT L1B file which the sounding was extracted from, appended with an additional three digits (0-indexed) to identify the number of the sounding within that L1B file. Equation (C1) shows the structure of the exposure_id. The exposure_id of 2009080100470440130006 was the seventh (006) sounding originating from the GOSAT L1B file named GOSATTFTS2009080100470440130_1BOB1D21 0210.01. The nomenclature for the GOSAT L1B file includes year (2009), month (08), day (01), hour (00), minute (47), orbit (044) and scene (0130). Note that these times are the times for the start of that GOSAT scene and not the exact measurement time. We provide the time variable as the measurement start time.

Appendix D: Application of GOSAT averaging kernels to model data

In order to correctly compare any model simulation to the satellite observations, the model data must be transformed to be consistent with assumptions made within the retrieval. Ultimately, this requires the satellite averaging kernels to be applied to the model data with any influence from the a priori data taken into account. The theory and methodology to do this are described in detail in Rodgers (2000) and we only briefly outline the method below. Equation (D1) is the equation which should be applied to any CH4 model data and details which variables provided in the data files are required to achieve this. It is assumed that any model data have already been interpolated to the same 20-level pressure grid (pressure_levels) as used in the retrieval. It should be noted here that this interpolation should be done with care to try and ensure that the model XCH4 is conserved via the interpolation process. Once on the same vertical grid as the GOSAT a priori and averaging kernels, Eq. (D1) can be applied to compute the modelled XCH4 by using the pressure_weight, xch4_averaging_kernel and ch4_profile_apriori variables provided in the file. It should be noted that we provide these values for each individual GOSAT sounding and that these are all level (i.e. layer boundary) quantities.

Appendix E: Considerations regarding GOSAT measurement strategy and changes over time

This section provides details on how the GOSAT measurement strategy has evolved over time. GOSAT initially operated primarily on a regular five-point grid, revisiting the same grid point. This changed over time to a three-point grid in order to reduce pointing uncertainty at the extreme angles when using the primary pointing mechanism. However, one consequence of a regular grid is that this resulted in a limited number of observations over some regions, especially islands, where the regular grid point would fall over the ocean. Thanks to the switch to the secondary pointing mechanism, several changes to the sampling strategy were possible.

Firstly, GOSAT is now capable of targeting more specifically and “follows” coastlines in a more efficient manner. The example in Fig. E1 shows the change in sampling location over Indonesia, contrasting 2011 to 2017. Although the exact same grid location is not revisited in the same way, overall there are both more successful measurements and a better geographic coverage.

Figure E1Figure contrasting the routine vs. target observation modes for land-only measurements in 2011 and 2017 over Indonesia. The change in pointing mechanism has allowed a more sophisticated targeting strategy to optimise the number of land measurements over islands and coastlines.

Figure E2Figure contrasting the sampling density pattern over Australia for early in the mission versus recent years. The five-point sampling grid as used in the early years of the missions was updated to a more stable three-point grid which continues to be used. The change in pointing mechanism has allowed the ocean sun-glint measurement range to be extended.

Secondly, GOSAT is now capable of a wider pointing range, and subsequently, the latitudinal range of ocean sun-glint observations has been extended (as observed in Fig. 10). Figure E2 provides example comparisons of the measurement density over Australia between the start of the mission (2009–2010) when the instrument was primarily operating in five-point grid mode to the latter years (2018–2019) where the instrument is operating in three-point grid mode, with extended latitudinal sun-glint coverage.

While we do not believe that these changes are detrimental to the continued consistency of the time series of GOSAT observations, we do feel that it is worth noting as they may have an impact (positively or negatively) on specific applications that a user may wish to use the data for, hence the reason for highlighting them here.

Figure E3Figure showing the number of successful GOSAT XCH4 measurements per 2 latitude–longitude bin from April 2009 to December 2019 separated into seasons – spring (MAM), summer (JJA), autumn (SON) and winter (DJF).

Appendix F: TCCON Data

The section provides details and references for the TCCON data used in this study. Table F1 lists all of the TCCON sites used in the study, along with their latitude, when the data record begins and the citation for each specific dataset.

Strong et al. (2019)Kivi et al. (2014)Wunch et al. (2018)Deutscher et al. (2014)Notholt et al. (2019)Hase et al. (2015)Té et al. (2014)Warneke et al. (2014)Sussmann and Rettinger (2018)Wennberg et al. (2017)Morino et al. (2018c)Wennberg et al. (2016)Goo et al. (2014)Morino et al. (2018a)Shiomi et al. (2014)Morino et al. (2018b)Feist et al. (2014)Griffith et al. (2014a)De Maziere et al. (2017)Griffith et al. (2014b)Sherlock et al. (2014b)Sherlock et al. (2014a)

Table F1The TCCON sites used in this study, along with their latitude, when they were established and the citation for the data used.

Download Print Version | Download XLSX

Author contributions

RJP developed and produced the Proxy XCH4 data, performed the analysis, and wrote the manuscript. AW assisted in the production of the data. AW, PS, ADN, HB, JSA, RBG and NK all contributed to development and analysis at different stages of the processing chain. All authors contributed towards discussion and interpretation of the analysis. PB, FC, PIP and LF provided model data and contributed to the interpretation of the comparisons. All TCCON co-authors provided TCCON data and contributed towards interpretation of the GOSAT–TCCON comparisons.

Competing interests

The authors declare that they have no conflict of interest.


Robert J. Parker, Hartmut Boesch, Alex Webb, Paul I. Palmer and Liang Feng are funded via the UK National Centre for Earth Observation (NE/R016518/1 and NE/N018079/1). Rocio Barrio Guillo and Nikoleta Kalaitzi were funded by a Leicester Institute for Space and EO (LISEO) and ESA-Dragon Programme studentship respectively. JSA was funded by an ESA Living Planet Fellowship. We acknowledge funding from the ESA GHG-CCI and Copernicus C3S projects. We thank the Japanese Aerospace Exploration Agency, National Institute for Environmental Studies and the Ministry of Environment for the GOSAT data and their continuous support as part of the Joint Research Agreement. This research used the ALICE High Performance Computing Facility at the University of Leicester for the GOSAT retrievals and analysis. The TM5-4DVAR CH4 inversions have been supported by ECMWF providing computing resources under the special project “Improve European and global CH4 and N2O flux inversions (2018–2020)”.

TCCON gratefully acknowledges financial support by ESA within the S5P validation programme. Stations at Park Falls, Lamont and Darwin are supported by NASA. Stations at Tsukuba, Rikubetsu and Burgos are supported in part by the GOSAT series project. Burgos is supported in part by the Energy Development Corp. Philippines. Ascension Island and Garmisch stations have been supported by the European Space Agency (ESA) under grant 4000120088/17/I-EF and by the German Bundesministerium für Wirtschaft und Energie (BMWi) under grants 50EE1711C, 50EE1711E and 50EE1711D. We thank the ESA Ariane Tracking Station at North East Bay, Ascension Island, for hosting and local support. The ETL station is funded by CFI/ORF, NSERC, ECCC and the CSA. The Paris station has received funding from Sorbonne Université, the French research center CNRS, the French space agency CNES and Région Île-de-France. The Eureka measurements were made at the Polar Environment Atmospheric Research Laboratory (PEARL) by the Canadian Network for the Detection of Atmospheric Change (CANDAC), primarily supported by the Natural Sciences and Engineering Research Council of Canada, Environment and Climate Change Canada, and the Canadian Space Agency. The Anmyeondo station has received funding from the Korea Meteorological Administration Research and Development Program ”Development and Assessment of IPCC AR6 Climate Change Scenario” under grant 1365003000. The Réunion station 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/IC4) and local activities supported by LACy/UMR8105 – Université de La Réunion. TCCON measurements in Australia are supported by NASA grants NAG5-12247 and NNH05-GD07G, Australian Research Council grants LE0668470, DP089468, DP110103118, DP140101552, DP160101598 and FT180100327, and the GOSAT series project.

Financial support

This research has been supported by the NERC National Centre for Earth Observation (grant nos. NE/N018079/1 and NE/R016518/1).

Review statement

This paper was edited by David Carlson and reviewed by three anonymous referees.


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,, 2015. a, b, c, d, e

Berchet, A., Pison, I., Chevallier, F., Paris, J.-D., Bousquet, P., Bonne, J.-L., Arshinov, M. Y., Belan, B. D., Cressot, C., Davydov, D. K., Dlugokencky, E. J., Fofonov, A. V., Galanin, A., Lavrič, J., Machida, T., Parker, R., Sasakawa, M., Spahni, R., Stocker, B. D., and Winderlich, J.: Natural and anthropogenic methane fluxes in Eurasia: a mesoscale quantification by generalized atmospheric inversion, Biogeosciences, 12, 5393–5414,, 2015. a

Bergamaschi, P., Houweling, S., Segers, A., Krol, M., Frankenberg, C., Scheepmaker, R. A., Dlugokencky, E., Wofsy, S. C., Kort, E. A., Sweeney, C., Schuck, T., Brenninkmeijer, C., Chen, H., Beck, V., and Gerbig, C.: Atmospheric CH4 in the first decade of the 21st century: Inverse modeling analysis using SCIAMACHY satellite retrievals and NOAA surface measurements, J. Geophys. Res.-Atmos., 118, 7350–7369,, 2013. a, b

Bergamaschi, P., Danila, A., Weiss, R. F., Ciais, P., Thompson, R. L., Brunner, D., Levin, I., Meijer, Y., Chevallier, F., Bovensmann, H., Crisp, D., Basu, S., Dlugokencky, E., Engelen, R., Gerbig, C., Günther, D., Hammer, S., Henne, S., Houweling, S., Peylin, P., Pinty, B., Ramonet, M., Reimann, S., Röckmann, T., Schmidt, M., Strogies, M., Sussams, J., Tarasova, O., van Aardenne, J., Vermeulen, A. T., and Vogel, F.: Atmospheric monitoring and inverse modelling for verification of greenhouse gas inventories, EUR 29276 EN, Publications Office of the European Union, Luxembourg,, 2018a. a

Bergamaschi, P., Karstens, U., Manning, A. J., Saunois, M., Tsuruta, A., Berchet, A., Vermeulen, A. T., Arnold, T., Janssens-Maenhout, G., Hammer, S., Levin, I., Schmidt, M., Ramonet, M., Lopez, M., Lavric, J., Aalto, T., Chen, H., Feist, D. G., Gerbig, C., Haszpra, L., Hermansen, O., Manca, G., Moncrieff, J., Meinhardt, F., Necki, J., Galkowski, M., O'Doherty, S., Paramonova, N., Scheeren, H. A., Steinbacher, M., and Dlugokencky, E.: Inverse modelling of European CH4 emissions during 2006–2012 using different inverse models and reassessed atmospheric observations, Atmos. Chem. Phys., 18, 901–920,, 2018b. a, b

Boesch, H., Baker, D., Connor, B., Crisp, D., and Miller, C.: Global Characterization of CO2 Column Retrievals from Shortwave-Infrared Satellite Observations of the Orbiting Carbon Observatory-2 Mission, Remote Sensing, 3, 270–304,, 2011. a, b

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 Mazière, 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,, 2017a. a, b, c, d, e, f

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,, 2017b. a, b

Buchwitz, M., Reuter, M., Schneising, O., Bovensmann, H., Burrows, J. P., Boesch, H., Anand, J., Parker, R., Detmers, R. G., Aben, I., Hasekamp, O. P., Crevoisier, C., Armante, R., Zehner, C., and Schepers, D.: Copernicus Climate Change Service (C3S) global satellite observations of atmospheric carbon dioxide and methane, Proceedings of the International Astronautical Congress, IAC, 2018, 57–60,, 2018. a, b

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.-Atmos., 115, D24302,, 2010. a

Chevallier, F.: Documentation of the CO2 flux service: Description of the CO2 inversion production, Tech. rep., European Centre for Medium-Range Weather Forecasts, 2019. a

Chevallier, F., Feng, L., Bösch, H., Palmer, P. I., and Rayner, P. J.: On the impact of transport model errors for the estimation of CO2 surface fluxes from GOSAT observations, Geophys. Res. Lett., 37, L21803,, 2010. a

Chipperfield, M. P.: Multiannual simulations with a three-dimensional chemical transport model, J. Geophys. Res.-Atmos., 104, 1781–1805,, 1999. a

Cogan, A. J., Boesch, H., Parker, R. J., Feng, L., Palmer, P. I., Blavier, J.-F. L. F., Deutscher, N. M., MacAtangay, R., Notholt, J., Roehl, C., Warneke, T., and Wunch, 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.-Atmos., 117, D21301,, 2012. a

Connor, B. J., Boesch, H., Toon, G., Sen, B., Miller, C., and Crisp, D.: Orbiting Carbon Observatory: Inverse method and prospective error analysis, J. Geophys. Res.-Atmos., 113, D05305,, 2008. a

Cressot, C., Chevallier, F., Bousquet, P., Crevoisier, C., Dlugokencky, E. J., Fortems-Cheiney, A., Frankenberg, C., Parker, R., Pison, I., Scheepmaker, R. A., Montzka, S. A., Krummel, P. B., Steele, L. P., and Langenfelds, R. L.: On the consistency between global and regional methane emissions inferred from SCIAMACHY, TANSO-FTS, IASI and surface measurements, Atmos. Chem. Phys., 14, 577–592,, 2014. a, b

De Mazière, M., Barret, B., Blumenstock, T., Buchwitz, M., De Beek, R., Demoulin, P., Fast, H., Gloudemans, A., Griesfeller, A., Griffith, D., Ionov, D., Janssens, K., Jones, N., Mahieu, E., Mellqvist, J., Mittermeier, R. L., Notholt, J., Rinsland, C., Schrijver, H., Schultz, A., Smale, D., Strandberg, A., Strong, K., Sussmann, R., Warneke, T., and Wood, S.: Comparisons between sciamachy scientific products and ground-based FTIR data for total columns of CO, CH4 and N2O, European Space Agency, (Special Publication) ESA SP, 89–96, 2004. a

De Maziere, 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 GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2017. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, I., Biblot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Greer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Holm, E. V., Isaksen, L., Kallberg, P., Kohler, M., Matricardi, M., McNally, A. P., Mong-Sanz, B. M., Morcette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thepaut, J. N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. a

Deutscher, N., Notholt, J., Messerschmidt, J., Weinzierl, C., Warneke, T., Petri, C., Grupe, P., and Katrynski, K.: TCCON data from Bialystok, Poland, Release GGG2014R2, TCCON data archive, hosted by CaltechDATA,, 2014. a

Dils, B., De Mazière, M., Müller, J. F., Blumenstock, T., Buchwitz, M., de Beek, R., Demoulin, P., Duchatelet, P., Fast, H., Frankenberg, C., Gloudemans, A., Griffith, D., Jones, N., Kerzenmacher, T., Kramer, I., Mahieu, E., Mellqvist, J., Mittermeier, R. L., Notholt, J., Rinsland, C. P., Schrijver, H., Smale, D., Strandberg, A., Straume, A. G., Stremme, W., Strong, K., Sussmann, R., Taylor, J., van den Broek, M., Velazco, V., Wagner, T., Warneke, T., Wiacek, A., and Wood, S.: Comparisons between SCIAMACHY and ground-based FTIR data for total columns of CO, CH4, CO2 and N2O, Atmos. Chem. Phys., 6, 1953–1976,, 2006. a

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,, 2014. a, b, c, d, e

Dlugokencky, E. J., Bruhwiler, L., White, J. W. C., Emmons, L. K., Novelli, P. C., Montzka, S. A., Masarie, K. A., Lang, P. M., Crotwell, A. M., Miller, J. B., and Gatti, L. V.: Observational constraints on recent increases in the atmospheric CH4 burden, Geophys. Res. Lett., 36, L18803,, 2009. a

ECMWF: Copernicus Atmospheric Data Store, available at:!/dataset/cams-global-greenhouse-gas-inversion, last access: 3 December 2020a. a

ECMWF: MACC-III greenhouse gases inversions, available at:, last access: 3 December 2020b. a

Etminan, M., Myhre, G., Highwood, E. J., and Shine, K. P.: Radiative forcing of carbon dioxide, methane, and nitrous oxide: A significant revision of the methane radiative forcing, Geophys. Res. Lett., 43, 12614–12623,, 2016. a

Feist, D. G., Arnold, S. G., John, N., and Geibel, M. C.: TCCON data from Ascension Island, Saint Helena, Ascension and Tristan da Cunha, Release GGG2014R0, TCCON data archive, hosted by CaltechDATA,, 2014. a

Feng, L., Palmer, P. I., Yang, Y., Yantosca, R. M., Kawa, S. R., Paris, J.-D., Matsueda, H., and Machida, T.: Evaluating a 3-D transport model of atmospheric CO2 using ground-based, aircraft, and space-borne data, Atmos. Chem. Phys., 11, 2789–2803,, 2011. a

Feng, L., Palmer, P. I., Bösch, H., Parker, R. J., Webb, A. J., Correia, C. S. C., Deutscher, N. M., Domingues, L. G., Feist, D. G., Gatti, L. V., Gloor, E., Hase, F., Kivi, R., Liu, Y., Miller, J. B., Morino, I., Sussmann, R., Strong, K., Uchino, O., Wang, J., and Zahn, A.: Consistent regional fluxes of CH4 and CO2 inferred from GOSAT proxy XCH4 : XCO2 retrievals, 2010–2014, Atmos. Chem. Phys., 17, 4781–4797,, 2017. a, b

Frankenberg, C., Meirink, J. F., Bergamaschi, P., Goede, A. P. H., Heimann, M., Körner, S., Platt, U., van Weele, M., and Wagner, T.: Satellite chartography of atmospheric methane from SCIAMACHY on board ENVISAT: Analysis of the years 2003 and 2004, J. Geophys. Res.-Atmos., 111, D07303,, 2006. a

Fraser, A., Palmer, P. I., Feng, L., Boesch, H., Cogan, A., Parker, R., Dlugokencky, E. J., Fraser, P. J., Krummel, P. B., Langenfelds, R. L., O'Doherty, S., Prinn, R. G., Steele, L. P., van der Schoot, M., and Weiss, R. F.: Estimating regional methane surface fluxes: the relative importance of surface and GOSAT mole fraction measurements, Atmos. Chem. Phys., 13, 5697–5713,, 2013. a, b

Fraser, A., Palmer, P. I., Feng, L., Bösch, H., Parker, R., Dlugokencky, E. J., Krummel, P. B., and Langenfelds, R. L.: Estimating regional fluxes of CO2 and CH4 using space-borne observations of XCH4 : XCO2, Atmos. Chem. Phys., 14, 12883–12895,, 2014. a

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,, 2017. a, b

Goo, T. Y., Oh, Y. S., and Velazco, V. A.: TCCON data from Anmyeondo (KOR), Release GGG2014R0, TCCON data archive, hosted by CaltechDATA,, 2014. a

Griffith, D. W. T., Deutscher, N., Velazco, V. A., Wennberg, P. O., Yavin, Y., Keppel Aleks, G., Washenfelder, R., Toon, G. C., Blavier, J.-F., Murphy, C., Jones, N., Kettlewell, G., Connor, B., Macatangay, R., Roehl, C., Ryczek, M., Glowacki, J., Culgan, T., and Bryant, G.: TCCON data from Darwin (AU), Release GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2014a. a

Griffith, D. W. T., Velazco, V. A., Deutscher, N., Murphy, C., Jones, N., Wilson, S., Macatangay, R., Kettlewell, G., Buchholz, R. R., and Riggenbach, M.: TCCON data from Wollongong (AU), Release GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2014b. a

Guanter, L., Kaufmann, H., Segl, K., Foerster, S., Rogass, C., Chabrillat, S., Kuester, T., Hollstein, A., Rossner, G., Chlebek, C., Straif, C., Fischer, S., Schrader, S., Storch, T., Heiden, U., Mueller, A., Bachmann, M., Mühle, H., Müller, R., Habermeyer, M., Ohndorf, A., Hill, J., Buddenbaum, H., Hostert, P., Van Der Linden, S., Leitão, P. J., Rabe, A., Doerffer, R., Krasemann, H., Xi, H., Mauser, W., Hank, T., Locherer, M., Rast, M., Staenz, K., and Sang, B.: The EnMAP spaceborne imaging spectroscopy mission for earth observation, Remote Sensing, 7, 8830–8857,, 2015. a

Gurney, K. R., Law, R. M., Denning, A. S., Rayner, P. J., Baker, D., Bousquet, P., Bruhwiler, L., Chen, Y.-H., Clals, P., Fan, S., Fung, I. Y., Gloor, M., Heimann, M., Higuchi, K., John, J., Maki, T., Maksyutov, S., Masarie, K., Peylin, P., Prather, M., Pak, B. C., Randerson, J., Sarmiento, J., Taguchi, S., Takahashi, T., and Yuen, C.-W.: Towards robust regional estimates of CO2 sources and sinks using atmospheric transport models, Nature, 415, 626–630,, 2002. a

Hase, F., Blumenstock, T., Dohe, S., Groß, J., and Kiel, M.: TCCON data from Karlsruhe (DE), Release GGG2014.R1, TCCON Data Archive, hosted by CaltechDATA,, 2015. a

Hausmann, P., Sussmann, R., and Smale, D.: Contribution of oil and natural gas production to renewed increase in atmospheric methane (2007–2014): top–down estimate from ethane and methane column observations, Atmos. Chem. Phys., 16, 3227–3244,, 2016. a

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,, 2018. a

Inoue, M., Morino, I., Uchino, O., Nakatsuru, T., Yoshida, Y., Yokota, T., Wunch, D., Wennberg, P. O., Roehl, C. M., Griffith, D. W. T., Velazco, V. A., Deutscher, N. M., Warneke, T., Notholt, J., Robinson, J., Sherlock, V., Hase, F., Blumenstock, T., Rettinger, M., Sussmann, R., Kyrö, E., Kivi, R., Shiomi, K., Kawakami, S., De Mazière, M., Arnold, S. G., Feist, D. G., Barrow, E. A., Barney, J., Dubey, M., Schneider, M., Iraci, L. T., Podolske, J. R., Hillyard, P. W., Machida, T., Sawa, Y., Tsuboi, K., Matsueda, H., Sweeney, C., Tans, P. P., Andrews, A. E., Biraud, S. C., Fukuyama, Y., Pittman, J. V., Kort, E. A., and Tanaka, T.: Bias corrections of GOSAT SWIR XCO2 and XCH4 with TCCON data and their evaluation using aircraft measurement data, Atmos. Meas. Tech., 9, 3491–3512,, 2016. a

Jacob, D. J., Turner, A. J., Maasakkers, J. D., Sheng, J., Sun, K., Liu, X., Chance, K., Aben, I., McKeever, J., and Frankenberg, C.: Satellite observations of atmospheric methane and their value for quantifying methane emissions, Atmos. Chem. Phys., 16, 14371–14396,, 2016. a

Kirschke, S., Bousquet, P., Ciais, P., et al.: Three decades of global methane sources and sinks, Nat. Geosci., 6, 813–823,, 2013. a

Kivi, R., Heikkinen, P., and Kyro, E.: TCCON data from Sodankyla (FI), Release GGG2014R0, TCCON data archive, hosted by CaltechDATA,, 2014. a

Kuze, A., Suto, H., Nakajima, M., and Hamazaki, T.: Thermal and near infrared sensor for carbon observation Fourier-transform spectrometer on the Greenhouse Gases Observing Satellite for greenhouse gases monitoring, Appl. Opt., 48, 6716–6733,, 2009. a, b

Kuze, A., Taylor, T. E., Kataoka, F., Bruegge, C. J., Crisp, D., Harada, M., Helmlinger, M., Inoue, M., Kawakami, S., Kikuchi, N., Mitomi, Y., Murooka, J., Naitoh, M., O'Brien, D. M., O'Dell, C. W., Ohyama, H., Pollock, H., Schwandner, F. M., Shiomi, K., Suto, H., Takeda, T., Tanaka, T., Urabe, T., Yokota, T., and Yoshida, Y.: Long-Term Vicarious Calibration of GOSAT Short-Wave Sensors: Techniques for Error Reduction and New Estimates of Radiometric Degradation Factors, IEEE T. Geosci. Remote, 52, 3991–4004,, 2014. a

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,, 2016. a

Lambert, J.-C., Keppens, A., Hubert, D., Langerock, B., Eichmann, K.-U., Kleipool, Q., Sneep, M., Verhoelst, T., Wagner, T., Weber, M., Ahn, C., Argyrouli, A., Balis, D., Chan, K., Compernolle, S., Smedt, I. D., Eskes, H., Fjæraa, A., K. Garane, J. G., Gouta, F., and Wang, P.: Sentinel-5 Precursor Mission Performance Centre Quarterly Validation Report of the Copernicus Sentinel-5 Precursor Operational Data Products # 03: July 2018–May 2019, Sentinel-5 Precursor Mission Performance Centre, 2019. a

Lu, X., Jacob, D. J., Zhang, Y., Maasakkers, J. D., Sulprizio, M. P., Shen, L., Qu, Z., Scarpelli, T. R., Nesser, H., Yantosca, R. M., Sheng, J., Andrews, A., Parker, R. J., Boech, H., Bloom, A. A., and Ma, S.: Global methane budget and trend, 2010–2017: complementarity of inverse analyses using in situ (GLOBALVIEWplus CH4 ObsPack) and satellite (GOSAT) observations, Atmos. Chem. Phys. Discuss.,, in review, 2020. a

Lunt, M. F., Palmer, P. I., Feng, L., Taylor, C. M., Boesch, H., and Parker, R. J.: An increase in methane emissions from tropical Africa between 2010 and 2016 inferred from satellite data, Atmos. Chem. Phys., 19, 14721–14740,, 2019. a, b

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,, 2019. a, b, c

Maasakkers, J. D., Jacob, D. J., Sulprizio, M. P., Scarpelli, T. R., Nesser, H., Sheng, J., Zhang, Y., Lu, X., Bloom, A. A., Bowman, K. W., Worden, J. R., and Parker, R. J.: 2010–2015 North American methane emissions, sectoral contributions, and trends: a high-resolution inversion of GOSAT satellite observations of atmospheric methane, Atmos. Chem. Phys. Discuss.,, in review, 2020. a

Matsunaga, T., Iwasaki, A., Tsuchida, S., Iwao, K., Tanii, J., Kashimura, O., Nakamura, R., Yamamoto, H., Kato, S., Obata, K., Mouri, K., and Tachikawa, T.: Current status of Hyperspectral Imager Suite (HISUI) onboard International Space Station (ISS), in: 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), 443–446,, 2017. a

McNorton, J., Chipperfield, M. P., Gloor, M., Wilson, C., Feng, W., Hayman, G. D., Rigby, M., Krummel, P. B., O'Doherty, S., Prinn, R. G., Weiss, R. F., Young, D., Dlugokencky, E., and Montzka, S. A.: Role of OH variability in the stalling of the global atmospheric CH4 growth rate from 1999 to 2006, Atmos. Chem. Phys., 16, 7943–7956,, 2016a. a

McNorton, J., Gloor, E., Wilson, C., Hayman, G. D., Gedney, N., Comyn-Platt, E., Marthews, T., Parker, R. J., Boesch, H., and Chipperfield, M. P.: Role of regional wetland emissions in atmospheric methane variability, Geophys. Res. Lett., 43, 11433–11444,, 2016b. a, b, c

McNorton, J., Wilson, C., Gloor, M., Parker, R. J., Boesch, H., Feng, W., Hossaini, R., and Chipperfield, M. P.: Attribution of recent increases in atmospheric methane through 3-D inverse modelling, Atmos. Chem. Phys., 18, 18149–18168,, 2018. a

Morino, I., Matsuzaki, T., and Horikawa, M.: TCCON data from Tsukuba (JP), 125HR, Release GGG2014.R2, TCCON Data Archive, hosted by CaltechDATA,, 2018a. a

Morino, I., Velazco, V. A., Hori, A., Uchino, O., and Griffith, D. W.: TCCON data from Burgos, Ilocos Norte (PH), Release GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2018b. a

Morino, I., Yokozeki, N., Matzuzaki, T., and Shishime, A.: TCCON data from Rikubetsu (JP), Release GGG2014.R2, TCCON Data Archive, hosted by CaltechDATA,, 2018c. a

Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and Natural Radiative Forcing, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G.-K., Tignor, M., Allen, S., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P., chap. 8, Cambridge University Press, 659–740,, 2013. a

Nisbet, E. G., Dlugokencky, E. J., Manning, M. R., Lowry, D., Fisher, R. E., France, J. L., Michel, S. E., Miller, J. B., White, J. W. C., Vaughn, B., and Bousquet, P., Pyle, J. A., Warwick, N. J., Cain, M., Brownlow, R., Zazzeri, G., Lanoisellé, M., Manning, A. C., Gloor, E., Worthy, D. E. J., Brunke, E.-G., Labuschagne, C., Wolff, E. W., and Ganesan, A. L.: Rising atmospheric methane: 2007–2014 growth and isotopic shift, Glob. Biogeochem. Cy., 30, 1356–1370, 2016. a

Nisbet, E. G., Fisher, R. E., Lowry, D., France, J. L., Allen, G., Bakkaloglu, S., Broderick, T. J., Cain, M., Coleman, M., Fernandez, J., Forster, G., Griffiths, P. T., Iverach, C. P., Kelly, B. F. J., Manning, M. R., Nisbet‐Jones, P. B. R., Pyle, J. A., Townsend‐Small, A., Al‐Shalaan, A., Warwick, N., and Zazzeri, G.: Methane Mitigation: Methods to Reduce Emissions, on the Path to the Paris Agreement, Rev. Geophys., 58, e2019RG000675,, 2020. a

NOAA: NOAA CarbonTracker model CO2 FTP, available at:, last access: 3 December 2020. a

Notholt, J., Petri, C., Warneke, T., Deutscher, N., Buschmann, M., Weinzierl, C., Macatangay, R., and Grupe, P.: TCCON data from Bremen (DE), Release GGG2014R1, TCCON data archive, hosted by CaltechDATA,, 2019. a

O'Dell, C. W.: Acceleration of multiple-scattering, hyperspectral radiative transfer calculations via low-streams interpolation, J. Geophys. Res.-Atmos., 115, D10206,, 2010. a

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,, 2012. a

Parker, R. and Boesch, H.: University of Leicester GOSAT Proxy XCH4 v9.0, Centre for Environmental Data Analysis,, 2020. a, b

Parker, R., Boesch, H., Cogan, A., Fraser, A., Feng, L., Palmer, P. I., Messerschmidt, J., Deutscher, N., Griffith, D. W. T., 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,, 2011. a, b, c

Parker, R. J., Boesch, H., Byckling, K., Webb, A. J., Palmer, P. I., Feng, L., Bergamaschi, P., Chevallier, F., Notholt, J., Deutscher, N., Warneke, T., Hase, F., Sussmann, R., Kawakami, S., Kivi, R., Griffith, D. W. T., and Velazco, V.: Assessing 5 years of GOSAT Proxy XCH4 data and associated uncertainties, Atmos. Meas. Tech., 8, 4785–4801,, 2015. a, b, c, d, e

Parker, R. J., Boesch, H., Wooster, M. J., Moore, D. P., Webb, A. J., Gaveau, D., and Murdiyarso, D.: Atmospheric CH4 and CO2 enhancements and biomass burning emission ratios derived from satellite observations of the 2015 Indonesian fire plumes, Atmos. Chem. Phys., 16, 10111–10131,, 2016. a, b

Parker, R. J., Boesch, H., McNorton, J., Comyn-Platt, E., Gloor, M., Wilson, C., Chipperfield, M. P., Hayman, G. D., and Bloom, A. A.: Evaluating year-to-year anomalies in tropical wetland methane emissions using satellite CH4 observations, Remote Sens. Environ., 211, 261–275,, 2018. a

Parker, R. J., Webb, A., Boesch, H., Somkuti, P., Barrio Guillo, R., Di Noia, A., Kalaitzi, N., Anand, J., Bergamaschi, P., Chevallier, F., Palmer, P. I., Feng, L., Deutscher, N. M., Feist, D. G., Griffith, D. W. T., Hase, F., Kivi, R., Morino, I., Notholt, J., Oh, Y.-S., Ohyama, H., Petri, C., Pollard, D. F., Roehl, C., Sha, M. K., Shiomi, K., Strong, K., Sussmann, R., Te, Y., Velazco, V. A., Warneke, T., Wennberg, P. O., and Wunch, D.: A Decade of GOSAT Proxy Satellite CH4 Observations, Earth Syst. Sci. Data Discuss.,, in review, 2020a. a

arker, R. J., Wilson, C., Bloom, A. A., Comyn-Platt, E., Hayman, G., McNorton, J., Boesch, H., and Chipperfield, M. P.: Exploring Constraints on a Wetland Methane Emission Ensemble (WetCHARTs) using GOSAT Satellite Observations, Biogeosciences Discuss.,, in review, 2020b. a

Patra, P. K., Houweling, S., Krol, M., Bousquet, P., Belikov, D., Bergmann, D., Bian, H., Cameron-Smith, P., Chipperfield, M. P., Corbin, K., Fortems-Cheiney, A., Fraser, A., Gloor, E., Hess, P., Ito, A., Kawa, S. R., Law, R. M., Loh, Z., Maksyutov, S., Meng, L., Palmer, P. I., Prinn, R. G., Rigby, M., Saito, R., and Wilson, C.: TransCom model simulations of CH4 and related species: linking transport, surface flux and chemical loss with CH4 variability in the troposphere and lower stratosphere, Atmos. Chem. Phys., 11, 12813–12837,, 2011. a

Peters, W., Jacobson, A. R., Sweeney, C., Andrews, A. E., Conway, 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., 104, 18925–18930,, 2007. a

Pignatti, S., Palombo, A., Pascucci, S., Romano, F., Santini, F., Simoniello, T., Umberto, A., Vincenzo, C., Acito, N., Diani, M., Matteoli, S., Corsini, G., Casa, R., De Bonis, R., Laneve, G., and Ananasso, C.: The PRISMA hyperspectral mission: Science activities and opportunities for agriculture and land monitoring, in: 2013 IEEE International Geoscience and Remote Sensing Symposium – IGARSS, 4558–4561,, 2013. a

Pollard, D. F., Sherlock, V., Robinson, J., Deutscher, N. M., Connor, B., and Shiona, H.: The Total Carbon Column Observing Network site description for Lauder, New Zealand, Earth Syst. Sci. Data, 9, 977–992,, 2017. a

Prather, M. J., Holmes, C. D., and Hsu, J.: Reactive greenhouse gas scenarios: Systematic exploration of uncertainties and the role of atmospheric chemistry, Geophys. Res. Lett., 39, L09803,, 2012. a

Reuter, M., Buchwitz, M., Schneising, O., Noël, S., Bovensmann, H., Burrows, J. P., Boesch, H., Di Noia, A., Anand, J., Parker, R. J., Somkuti, P., Wu, L., Hasekamp, O. P., Aben, I., Kuze, A., Suto, H., Shiomi, K., Yoshida, Y., Morino, I., Crisp, D., O'Dell, C. W., Notholt, J., Petri, C., Warneke, T., Velazco, V. A., Deutscher, N. M., Griffith, D. W. T., Kivi, R., Pollard, D. F., Hase, F., Sussmann, R., Té, Y. V., Strong, K., Roche, S., Sha, M. K., De Mazière, M., Feist, D. G., Iraci, L. T., Roehl, C. M., Retscher, C., and Schepers, D.: Ensemble-based satellite-derived carbon dioxide and methane column-averaged dry-air mole fraction data sets (2003–2018) for carbon and climate applications, Atmos. Meas. Tech., 13, 789–819,, 2020. a, b

Rigby, M., Prinn, R. G., Fraser, P. J., Simmonds, P. G., Langenfelds, R. L., Huang, J., Cunnold, D. M., Steele, L. P., Krummel, P. B., Weiss, R. F., O'Doherty, S., Salameh, P. K., Wang, H. J., Harth, C. M., Muhle, J., and Porter, L. W.: Renewed growth of atmospheric methane, Geophys. Res. Lett., 35, L22805,, 2008. a

Rigby, M., Montzka, S. A., Prinn, R. G., White, J. W. C., Young, D., O'Doherty, S., Lunt, M. F., Ganesan, A. L., Manning, A. J., Simmonds, P. G., Salameh, P. K., Harth, C. M., Mühle, J., Weiss, R. F., Fraser, P. J., Steele, L. P., Krummel, P. B., McCulloch, A., Park, S., O'Doherty, S., Lunt, M. F., Ganesan, A. L., Manning, A. J., Simmonds, P. G., Salameh, P. K., Harth, C. M., Mühle, J., Weiss, R. F., Fraser, P. J., Steele, L. P., Krummel, P. B., McCulloch, A., Park, S., O'Doherty, S., Lunt, M. F., Ganesan, A. L., Manning, A. J., Simmonds, P. G., Salameh, P. K., Harth, C. M., Mühle, J., Weiss, R. F., Fraser, P. J., Steele, L. P., Krummel, P. B., McCulloch, A., and Park, S.: Role of atmospheric oxidation in recent methane growth, P. Natl. Acad. Sci. USA, 114, 5373–5377,, 2017. a

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding, Series on Atmospheric, Oceanic and Planetary Physics: Volume 2, World Scientific, 256 pp.,, 2000. a

Ross, A. N., Wooster, M. J., Boesch, H., and Parker, R.: First satellite measurements of carbon dioxide and methane emission ratios in wildfire plumes, Geophys. Res. Lett., 40, 4098–4102,, 2013. a, b

Saad, K. M., Wunch, D., Deutscher, N. M., Griffith, D. W. T., Hase, F., De Mazière, M., Notholt, J., Pollard, D. F., Roehl, C. M., Schneider, M., Sussmann, R., Warneke, T., and Wennberg, P. O.: Seasonal variability of stratospheric methane: implications for constraining tropospheric methane budgets using total column observations, Atmos. Chem. Phys., 16, 14003–14024,, 2016. a

Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Brovkin, V., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Curry, C., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., McDonald, K. C., Marshall, J., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prigent, C., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Steele, P., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., van Weele, M., van der Werf, G. R., Weiss, R., Wiedinmyer, C., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: The global methane budget 2000–2012, Earth Syst. Sci. Data, 8, 697–751,, 2016. a

Saunois, M., Stavert, A. R., Poulter, B., Bousquet, P., Canadell, J. G., Jackson, R. B., Raymond, P. A., Dlugokencky, E. J., Houweling, S., Patra, P. K., Ciais, P., Arora, V. K., Bastviken, D., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Carlson, K. M., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., Crill, P. M., Covey, K., Curry, C. L., Etiope, G., Frankenberg, C., Gedney, N., Hegglin, M. I., Höglund-Isaksson, L., Hugelius, G., Ishizawa, M., Ito, A., Janssens-Maenhout, G., Jensen, K. M., Joos, F., Kleinen, T., Krummel, P. B., Langenfelds, R. L., Laruelle, G. G., Liu, L., Machida, T., Maksyutov, S., McDonald, K. C., McNorton, J., Miller, P. A., Melton, J. R., Morino, I., Müller, J., Murguia-Flores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., Parker, R. J., Peng, C., Peng, S., Peters, G. P., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., Riley, W. J., Rosentreter, J. A., Segers, A., Simpson, I. J., Shi, H., Smith, S. J., Steele, L. P., Thornton, B. F., Tian, H., Tohjima, Y., Tubiello, F. N., Tsuruta, A., Viovy, N., Voulgarakis, A., Weber, T. S., van Weele, M., van der Werf, G. R., Weiss, R. F., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The Global Methane Budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623,, 2020. a, b, c, d

Schaefer, H., Fletcher, S. E. M., Veidt, C., Lassey, K. R., Brailsford, G. W., Bromley, T. M., Dlugokencky, E. J., Michel, S. E., Miller, J. B., Levin, I., Lowe, D. C., Martin, R. J., Vaughn, B. H., and White, J. W. C.: A 21st-century shift from fossil-fuel to biogenic methane emissions indicated by 13CH4, Science, 352, 80–84,, 2016. a

Schepers, D., Guerlet, S., Butz, A., Landgraf, J., Frankenberg, C., Hasekamp, O., Blavier, J.-F., Deutscher, N. M., Griffith, D. W. T., Hase, F., Kyro, E., Morino, I., Sherlock, V., Sussmann, R., and Aben, I.: Methane retrievals from Greenhouse Gases Observing Satellite (GOSAT) shortwave infrared measurements: Performance comparison of proxy and physics retrieval algorithms, J. Geophys. Res.-Atmos., 117, D10307,, 2012. a

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,, 2019. a

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,, 2018. a, b

Sherlock, V., Connor, B., Robinson, J., Shiona, H., Smale, D., and Pollard, D.: TCCON data from Lauder (NZ), Release GGG2014R0, TCCON Data Archive, hosted by CaltechDATA,, 2014a. a

Sherlock, V., Connor, B., Robinson, J., Shiona, H., Smale, D., and Pollard, D.: TCCON data from Lauder (NZ), Release GGG2014R0, TCCON data archive, hosted by CaltechDATA,, 2014b. a

Shiomi, K., Kawakami, S., Ohyama, H., Arai, K., Okumura, H., Taura, C., Fukamachi, T., and Sakashita, M.: TCCON data from Saga (JP), Release GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2014. a

Siddans, R., Knappett, D., Kerridge, B., Waterfall, A., Hurley, J., Latter, B., Boesch, H., and Parker, R.: Global height-resolved methane retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) on MetOp, Atmos. Meas. Tech., 10, 4135–4164,, 2017. a, b

Stanevich, I., Jones, D. B. A., Strong, K., Keller, M., Henze, D. K., Parker, R. J., Boesch, H., Wunch, D., Notholt, J., Petri, C., Warneke, T., Sussmann, R., Schneider, M., Hase, F., Kivi, R., Deutscher, N. M., Velazco, V. A., Walker, K. A., and Deng, F.: Characterizing model errors in chemical transport modelling of methane: Using GOSAT XCH4 data with weak constraint four-dimensional variational data assimilation, Atmos. Chem. Phys. Discuss.,, in review, 2019. a

tanevich, I., Jones, D. B. A., Strong, K., Parker, R. J., Boesch, H., Wunch, D., Notholt, J., Petri, C., Warneke, T., Sussmann, R., Schneider, M., Hase, F., Kivi, R., Deutscher, N. M., Velazco, V. A., Walker, K. A., and Deng, F.: Characterizing model errors in chemical transport modeling of methane: impact of model resolution in versions v9-02 of GEOS-Chem and v35j of its adjoint model, Geosci. Model Dev., 13, 3839–3862,, 2020. a

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,, 2019. a

Sussmann, R. and Rettinger, M.: TCCON data from Garmisch (DE), Release GGG2014R2, TCCON Data Archive, hosted by CaltechDATA,, 2018. a

Suto, H., Kataoka, F., Kikuchi, N., Knuteson, R. O., Butz, A., Haun, M., Buijs, H., Shiomi, K., Imai, H., and Kuze, A.: Thermal and near-infrared sensor for carbon observation Fourier-transform spectrometer-2 (TANSO-FTS-2) on the Greenhouse Gases Observing Satellite-2 (GOSAT-2) during its first year on orbit, Atmos. Meas. Tech. Discuss.,, in review, 2020. a

TCCON Data Archive: available at:, last access: 10 December 2020. a

Té, Y., Jeseck, P., and Janssen, C.: TCCON data from Paris (FR), Release GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2014. a

Toon, G. C.: Atmospheric Line List for the 2014 TCCON Data Release,, 2015. a

Tunnicliffe, R. L., Ganesan, A. L., Parker, R. J., Boesch, H., Gedney, N., Poulter, B., Zhang, Z., Lavrič, J. V., Walter, D., Rigby, M., Henne, S., Young, D., and O'Doherty, S.: Quantifying sources of Brazil's CH4 emissions between 2010 and 2018 from satellite data, Atmos. Chem. Phys., 20, 13041–13067,, 2020. a

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,, 2015. a, b, c

Turner, A. J., Frankenberg, C., Wennberg, P. O., and Jacob, D. J.: Ambiguity in the causes for decadal trends in atmospheric methane and hydroxyl, P. Natl. Acad. Sci. USA, 114, 5367–5372,, 2017. a

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.: TCCON Philippines: First measurement results, satellite data and model comparisons in Southeast Asia, Remote Sensing, 9, 1–18,, 2017. a

Wang, Z., Warneke, T., Deutscher, N. M., Notholt, J., Karstens, U., Saunois, M., Schneider, M., Sussmann, R., Sembhi, H., Griffith, D. W. T., Pollard, D. F., Kivi, R., Petri, C., Velazco, V. A., Ramonet, M., and Chen, H.: Contributions of the troposphere and stratosphere to CH4 model biases, Atmos. Chem. Phys., 17, 13283–13295,, 2017. a

Warneke, T., Messerschmidt, J., Notholt, J., Weinzierl, 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 GGG2014.R0, TCCON Data Archive, hosted by CaltechDATA,, 2014. a

Webb, A. J., Bösch, H., Parker, R. J., Gatti, L. V., Gloor, E., Palmer, P. I., Basso, L. S., Chipperfield, M. P., Correia, C. S. C., Domingues, L. G., Feng, L., Gonzi, S., Miller, J. B., Warneke, T., and Wilson, C.: CH4 concentrations over the Amazon from GOSAT consistent with in situ vertical profile data, J. Geophys. Res.-Atmos., 121, 11006–11020,, 2016. a, b

Wecht, K. J., Jacob, D. J., Sulprizio, M. P., Santoni, G. W., Wofsy, S. C., Parker, R., Bösch, H., and Worden, J.: Spatially resolving methane emissions in California: constraints from the CalNex aircraft campaign and from present (GOSAT, TES) and future (TROPOMI, geostationary) satellite observations, Atmos. Chem. Phys., 14, 8173–8184,, 2014. a, b

Wennberg, P. O., Wunch, D., Roehl, C., Blavier, J.-F., Toon, G. C., Allen, N., Dowell, P., Teske, K., Martin, C., and Martin, J.: TCCON data from Lamont (US), Release GGG2014.R1, TCCON Data Archive, hosted by CaltechDATA,, 2016. a

Wennberg, P. O., Roehl, C., Wunch, D., Toon, G. C., Blavier, J. F., Washenfelder, R., Keppel-Aleks, G., Allen, N., and Ayers, J.: TCCON data from Park Falls (USA), Release GGG2014R1, TCCON data archive, hosted by CaltechDATA,, 2017. a

Worden, J. R., Turner, A. J., Bloom, A., Kulawik, S. S., Liu, J., Lee, M., Weidner, R., Bowman, K., Frankenberg, C., Parker, R., and Payne, V. H.: Quantifying lower tropospheric methane concentrations using GOSAT near-IR and TES thermal IR measurements, Atmos. Meas. Tech., 8, 3433–3445,, 2015. a, b

Worden, J. R., Bloom, A. A., Pandey, S., Jiang, Z., Worden, H. M., Walker, 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, 1–11,, 2017. a

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,, 2010.  a, b

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,, 2011. a

Wunch, D., Mendonca, J., Colebatch, O., Allen, N., Blavier, J.-F. L., Roche, S., Hedelius, J., Neufeld, G., Springett, S., Worthy, D., Kessler, R., and Strong, K.: TCCON data from East Trout Lake (CA), Release GGG2014.R1, TCCON Data Archive, hosted by CaltechDATA,, 2018. a

Yin, Y., Chevallier, F., Ciais, P., Bousquet, P., Saunois, M., Zheng, B., Worden, J., Bloom, A. A., Parker, R., Jacob, D., Dlugokencky, E. J., and Frankenberg, C.: Accelerating methane growth rate from 2010 to 2017: leading contributions from the tropics and East Asia, Atmos. Chem. Phys. Discuss.,, in review, 2020. a

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,, 2013. a, b, c

Zhang, Y., Jacob, D. J., Lu, X., Maasakkers, J. D., Scarpelli, T. R., Sheng, J.-X., Shen, L., Qu, Z., Sulprizio, M. P., Chang, J., Bloom, A. A., Ma, S., Worden, J., Parker, R. J., and Boesch, H.: Attribution of the accelerating increase in atmospheric methane during 2010002018 by inverse analysis of GOSAT observations, Atmos. Chem. Phys. Discuss.,, in review, 2020. a

Zheng, B., Chevallier, F., Yin, Y., Ciais, P., Fortems-Cheiney, A., Deeter, M. N., Parker, R. J., Wang, Y., Worden, H. M., and Zhao, Y.: Global atmospheric carbon monoxide budget 2000–2017 inferred from multi-species atmospheric inversions, Earth Syst. Sci. Data, 11, 1411–1436,, 2019. a

Short summary
This work presents the latest release of the University of Leicester GOSAT methane data and acts as the definitive description of this dataset. We detail the processing, validation and evaluation involved in producing these data and highlight its many applications. With now over a decade of global atmospheric methane observations, this dataset has helped, and will continue to help, us better understand the global methane budget and investigate how it may respond to a future changing climate.