Articles | Volume 23, issue 14
Research article
19 Jul 2023
Research article |  | 19 Jul 2023

East Asian methane emissions inferred from high-resolution inversions of GOSAT and TROPOMI observations: a comparative and evaluative analysis

Ruosi Liang, Yuzhong Zhang, Wei Chen, Peixuan Zhang, Jingran Liu, Cuihong Chen, Huiqin Mao, Guofeng Shen, Zhen Qu, Zichong Chen, Minqiang Zhou, Pucai Wang, Robert J. Parker, Hartmut Boesch, Alba Lorente, Joannes D. Maasakkers, and Ilse Aben

We apply atmospheric methane column retrievals from two different satellite instruments (Greenhouse gases Observing SATellite – GOSAT; TROPOspheric Monitoring Instrument – TROPOMI) to a regional inversion framework to quantify East Asian methane emissions for 2019 at 0.5× 0.625 horizontal resolution. The goal is to assess if GOSAT (relatively mature but sparse) and TROPOMI (new and dense) observations inform consistent methane emissions from East Asia with identically configured inversions. Comparison of the results from the two inversions shows similar correction patterns to the prior inventory in central northern China, central southern China, northeastern China, and Bangladesh, with less than 2.6 Tg a−1 differences in regional posterior emissions. The two inversions, however, disagree over some important regions, particularly in northern India and eastern China. The methane emissions inferred from GOSAT observations are 7.7 Tg a−1 higher than those from TROPOMI observations over northern India but 6.4 Tg a−1 lower over eastern China. The discrepancies between the two inversions are robust against varied inversion configurations (i.e., assimilation window and error specifications). We find that the lower methane emissions from eastern China inferred by the GOSAT inversion are more consistent with independent ground-based in situ and total column (TCCON) observations, indicating that the TROPOMI retrievals may have high XCH4 biases in this region. We also evaluate inversion results against tropospheric aircraft observations over India during 2012–2014 by using a consistent GOSAT inversion of earlier years as an intercomparison platform. This indirect evaluation favors lower methane emissions from northern India inferred by the TROPOMI inversion. We find that in this case the discrepancy in emission inference is contributed by differences in data coverage (almost no observations by GOSAT vs. good spatial coverage by TROPOMI) over the Indo-Gangetic Plain. The two inversions also differ substantially in their posterior estimates for northwestern China and neighboring Kazakhstan, which is mainly due to seasonally varying biases between GOSAT and TROPOMI XCH4 data that correlate with changes in surface albedo.

1 Introduction

Methane (CH4) is a powerful greenhouse gas, with a global warming potential ∼80 times that of carbon dioxide (CO2) on a 20-year timescale and ∼30 times on a 100-year timescale (Forster et al., 2021). In 2020, the atmospheric methane concentration increased to 1889±2 ppbv, 262 % of pre-industrial levels in 1750, driven primarily by increasing anthropogenic emissions (WMO, 2021). The last decade has seen rapid growth of atmospheric methane (∼8.6 ppbv a−1), after a brief period of stabilization in the early 2000s (Dlugokencky et al., 2011; Fletcher and Schaefer, 2019; Rigby et al., 2008; Yin et al., 2021; Zhang et al., 2021). Rising methane concentrations, if continued at current rates in coming decades, may negate benefits of CO2 emission reduction, and therefore curbing methane emissions in the 2020s is vital for the success of the Paris Agreement (Ganesan et al., 2019; Nisbet et al., 2019).

Information on methane emissions is required at global, national, and regional levels to guide climate actions on methane. Current bottom-up inventories are often inadequate for this purpose because of their large uncertainties in emission factors and lack of information on emission activities (Saunois et al., 2020). Independent measurements of atmospheric methane, including those from satellite remote sensing, are thus used to evaluate and improve these bottom-up inventories (Jacob et al., 2016). This is generally done through an inversion of atmospheric observations with a chemical transport model (CTM) to characterize the relationship between emissions and concentrations. Atmospheric methane is measured by two classes of satellite instruments, point source imagers and area flux mappers. While point sources imagers (e.g., Sentinel-2, Landsat, GHGSat) specialize in detecting large emissions from point sources, area flux mappers provide high-precision measurements that can be used to constrain methane fluxes on regional and global scales (Jacob et al., 2022). Area flux mappers that are currently in operation include the TANSO-FTS instrument on board the Greenhouse gases Observing SATellite (GOSAT) launched in 2009 (Kuze et al., 2016) and the more recent TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor (S5P) satellite launched in 2017 (Hu et al., 2016; Lorente et al., 2021a; Veefkind et al., 2012). Satellite observations made by these area flux mappers are especially valuable in constraining methane emissions over regions with no or only sparse ground networks, including Africa, South America, and East and South Asia (Lu et al., 2021).

Both GOSAT and TROPOMI operate in sun-synchronous orbits and retrieve column-averaged dry-air methane mole fractions (XCH4) from backscattered solar shortwave infrared radiation. TROPOMI continuously images the land surface at a pixel resolution of 7 km × 7 km (5.5 km × 7 km after August 2019) with daily global coverage (Hu et al., 2018; Lorente et al., 2021a; Sha et al., 2021), while GOSAT in its standard viewing mode measures with a 3 d return time in 10 km diameter circular footprints that are typically spaced ∼250 km apart (Butz et al., 2011; Kuze et al., 2009, 2016; Yokota et al., 2009). As a result of differing sampling strategies, TROPOMI generates much higher observation density than GOSAT, which in principle should benefit fine-resolution inversions. The two instruments also measure at different wavelengths: GOSAT at the 1.65 µm band and TROPOMI at the 2.3 µm band. This affects the algorithm that can be applied to retrieve XCH4. Operational TROPOMI retrievals use the RemoTeC full-physics method (Hu et al., 2018). The method is prone to spatially and temporally variable biases owing to scattering artifacts (Hu et al., 2018; Lorente et al., 2021a; Sha et al., 2021). These biases in general are not reducible with more observations and, if not corrected, can translate into biases in emission estimates in an inversion. Because of spectrally adjacent CO2 and CH4 absorption in the 1.65 µm band, GOSAT retrievals can alternatively use the CO2 proxy method, in which XCH4 is derived from directly retrieved CH4 to CO2 column ratios and independently specified (simulated or assimilated) CO2 columns (Alexe et al., 2015; Frankenberg et al., 2005, 2006; Parker et al., 2015, 2020). The proxy method usually results in reduced variable biases, as scattering artifacts largely cancel out in retrieving CH4 to CO2 column ratios. It is, however, subject to any errors in specified CO2 columns. The proxy method also leads to a better retrieval success rate over regions with high aerosol loadings or thin clouds, as the method is less sensitive to this interference compared to the full-physics approach.

A number of studies have applied GOSAT data in inversions on a range of scales (Alexe et al., 2015; Cressot et al., 2014; Feng et al., 2022; Lu et al., 2021; Maasakkers et al., 2019; Monteil et al., 2013; Pandey et al., 2016; Turner et al., 2015; Zhang et al., 2021). TROPOMI data have also been applied in several regional inversion studies (Chen et al., 2022; McNorton et al., 2022; Shen et al., 2021, 2022; Zhang et al., 2020), often with the focus on resolving fine-scale emission hotspots. Qu et al. (2021) performed global inversions of GOSAT and TROPOMI observations at 2× 2.5 resolution in a comparative analysis, and they showed that methane emissions inferred from the two inversions are generally consistent on the global scale but with significant regional discrepancies including over China.

Here we present high-resolution (0.5× 0.625) inversions of GOSAT and TROPOMI observations over East Asia. The main objective is to assess the consistency of methane fluxes inferred from the two sets of satellite data that differ in their data coverage and regional bias (Qu et al., 2021), adding information to the uncertainty characterization of satellite-based methane emission accounting. We perform the analyses with identically configured inversions to isolate the effects of observation data, and we further use independent ground-based observations to evaluate the discrepancies between the two inversions and discuss the cause of differences. This study focuses on East Asia (including China and northern India), which is one of the world's major methane-emitting regions and accounts for more than 20 % of global emissions (UNFCCC, 2020). The region has been an important contributor to global increases in methane emissions, but the magnitude of the trend and its sectoral attributions are debated (Ganesan et al., 2017; Gao et al., 2021; Liu et al., 2021; Miller et al., 2019; Sheng et al., 2021; Zhang et al., 2021).

2 Observation data

2.1 Satellite observations

We used XCH4 observations from GOSAT and TROPOMI for 2019 in regional inversions over East Asia. For GOSAT, we use the University of Leicester Proxy XCH4 v9.0 retrievals (Parker and Boesch, 2020). Our inversion assimilates only high-quality GOSAT retrievals flagged as “xch4_quality_flag = 0” over both land and ocean (glint mode). This GOSAT product is based on the CO2 proxy method, which use the ratio between simulated (XCO2model) and retrieved (XCO2raw) CO2 columns to correct for retrieved methane columns (XCH4raw) that are sensitive to aerosol and surface interference:

(1) XCH 4 = XCH 4 raw XCO 2 raw × XCO 2 model .

This limits variable biases because both XCH4raw and XCO2raw are similarly affected by scattering artifacts, but the method is subject to any biases in specified CO2 columns (Parker et al., 2015).

The University of Leicester Proxy XCH4 v9.0 retrieval takes the median CO2 columns from three atmospheric chemistry transport models as XCO2model, and the range of the three models characterizes the XCO2model uncertainty (Fig. S1). The disagreement among these three models is ∼1 ppm over remote regions, ∼2 ppm over eastern China, and 2–4 ppm in India, Bangladesh, and southwestern China, which roughly translates to uncertainties of 0.3, 0.4, and 0.5 %–1.0 %, respectively, in retrieved XCH4.

For TROPOMI, we use the SRON RemoTeC-S5P XCH4 scientific product from Lorente et al. (2021a). The improved algorithm by Lorente et al. (2021a) was later implemented in the official operational product (v2.02.00) in July 2021 (Lorente et al., 2022). They derived an empirical correction formula to improve surface-reflectance-dependent biases identified in TROPOMI full-physics retrievals. The correction in general improves data quality over scenes with low (e.g., snow cover) and high surface albedo (e.g., deserts), which are challenging for a full-physics algorithm. Large corrections are made in eastern China, Xinjiang (China), Southeast Asia, and Siberia (Fig. S2). Bias-corrected TROPOMI retrievals flagged with “qa_value = 1” are used for inversion. This version of the TROPOMI product does not provide glint-mode ocean retrievals.

Figure 1 shows the spatial distributions of XCH4 measured by GOSAT and TROPOMI, annually averaged on the 0.625× 0.5 grid. Both datasets show high XCH4 in eastern China and northern India and low XCH4 over the Mongolian and Tibetan plateaus, although TROPOMI provides much better spatial coverage than GOSAT over most regions. There are in total 45 018 observations for GOSAT and 8 860 722 for TROPOMI. We take averages when multiple measurements fall within a 0.625×0.5 grid cell on any individual day (this procedure affects primarily dense TROPOMI data), and the resulting gridded daily observations are used in the inversion. The spatial distribution of gridded daily observation numbers is shown in Fig. S3.

We refer to the XCH4 retrieval products used in this study as GOSAT or TROPOMI observations and corresponding inversions as GOSAT or TROPOMI inversions for simplicity. There are other operational and scientific retrieval products available from both GOSAT and TROPOMI measurements (e.g., the operational GOSAT XCH4 retrieval, Yoshida et al., 2013; the scientific TROPOMI/WFMD XCH4 retrieval, Schneising et al., 2019). Our analyses and conclusions are specific to the two retrieval products used here, though we expect that some of them can also apply to other retrievals.

Figure 1The 2019 annual average methane column mole fractions over the East Asia domain for GOSAT (University of Leicester Proxy v9.0 retrieval) and TROPOMI (Lorente et al., 2021a, full-physics retrieval), presented on the 0.5× 0.625 GEOS-Chem grid.

2.2 Independent evaluation data

We use a suite of independent high-quality methane observations to evaluate the posterior emissions inferred from satellite observations, including surface in situ observations, ground-based remote sensing observations, and tropospheric in situ measurements from commercial airlines. Table S1 provides a descriptive list of these surface sites, and Fig. 2 shows the locations of surface sites and a representative flight path. These suborbital observations are of good accuracy and precision compared to satellite observations.

Surface in situ observations are available through World Data Centre for Greenhouse Gases (WDCGG) or the CH4 GLOBALVIEWplus v4.0 ObsPack (Schuldt et al., 2021). The five sites are Anmyeon-do (AMY), Pha Din (PDI), Lulin (LLN), Ulaan Uul (UUM), and Waliguan (WLG) (Dlugokencky et al., 1994, 2021; Lee et al., 2019; Nguyen and Steinbacher, 2021). Observations are done with either continuous (hourly) online instruments or weekly collected flasks (Table S1). Daytime measurements are used for evaluating simulations. Most of these sites are continental or subcontinental background sites (PDI, LLN, UUM, and WLG), and their observations are insensitive to local methane emissions. An exception is AMY, which is affected by local Korean emissions as well as upwind eastern China emissions.

Total methane column observations by ground-based Fourier transform infrared (FTIR) spectrometers are available at two TCCON sites located in eastern China, Hefei (HF) and Xianghe (XH) (Liu et al., 2023; Yang et al., 2020), and their observations are sensitive to methane emissions from eastern China. We note that a previous evaluation of GOSAT and TROPOMI against TCCON did not include data from these two sites, as their data were not available then (Qu et al., 2021). We use only measurements with solar zenith angles <60 to ensure high data quality.

All the above surface sites are far from northern India, which is a major methane-emitting region in the study domain. The only relevant dataset available to us in this area comes from the Civil Aircraft for the Regular Investigation of the atmosphere Based on an Instrument Container (CARIBIC) project (available via the CH4 GLOBALVIEWplus v4.0 ObsPack, Schuldt et al., 2021), which includes regular flights in the troposphere over northern India. However, these data were collected in earlier years between 2012 and 2014 before the time of TROPOMI. In the absence of better observation data, we compare these 2012–2014 aircraft observations to a simulation driven by a similarly configured GOSAT inversion for an earlier period (2010–2017) (Zhang et al., 2022). By doing so, we assume that any systematic bias derived from this comparison should still be representative of the 2019 GOSAT inversion.

Figure 2Spatial distribution of prior emissions. Locations of independent data for evaluation (seven surface sites and aircraft route) are shown. Circles represent background sites and triangles source-region sites. Total column measurements are coded in red and in situ measurements in blue. The solid purple line shows a CARIBIC aircraft route that measured tropospheric methane over India on 22 November 2012.

3 Inverse analysis

3.1 Forward model and prior emissions

We use GEOS-Chem v12.9.3 as the forward model for the inversion. The simulation is conducted for 2019 over East Asia (15–55 N, 60–140 E) on a 0.5×0.625 horizontal grid with 47 vertical layers and is driven by MERRA-2 meteorological fields from the NASA Global Modeling and Assimilation Office (GMAO) (Gelaro et al., 2017). The initial concentration fields on 1 January 2019 and 3-hourly boundary conditions for the nested domain are taken from a global inversion of TROPOMI data for 2019 (Qu et al., 2021). We find that the boundary conditions from this global inversion still have biases over East Asia (discussed further in Sect. 4.3.3), which may partly be due to the fact that Qu et al. (2021) used an earlier version of TROPOMI retrievals. In our inversion, we optimize for systematic biases at four lateral boundaries together with methane emissions.

Prior emissions (Fig. 2) used in GEOS-Chem simulations are compiled from bottom-up sectoral inventories (Table S2). In brief, we use EDGAR v4.3.2 (Janssen-Maenhout et al., 2019) for anthropogenic methane emissions, with those from fossil fuel exploitation replaced by Scarpelli et al. (2020) (oil and gas; coal outside of China) and Sheng et al. (2019) (coal in China). A comparison with a more recent inventory EDGAR v6 shows no large revisions of anthropogenic methane emissions over the study region that we expect to have a great impact on the inversion results (Fig. S4). For natural emissions, we use the ensemble average of the WetCHARTs version 1.0 inventory for wetlands (Bloom et al., 2017), the Quick Fire Emissions Dataset (QFED) v2.4r8 for biomass burning, Fung et al. (1991) for termite emissions, and Maasakkers et al. (2019) for geological sources.

While methane sinks are not optimized in our regional inversion, they are explicitly simulated in GEOS-Chem simulations. We use monthly OH fields from a full-chemistry GEOS-Chem simulation (Wecht et al., 2014) and soil absorption from Murguia-Flores et al. (2018).

3.2 Inversion procedure

We perform analytical Bayesian inversions to optimize a state vector x containing annual methane emissions from 600 clusters and average methane column biases at four model boundaries. We optimize emissions on 600 spatial clusters instead of the native 0.5×0.625 grid (Fig. S5), which are generated based on a Gaussian mixed model (GMM) algorithm proposed by Turner and Jacob (2015). This strategy significantly reduces the computation of an analytical inversion while accounting for major patterns in the distribution of methane emissions. We also optimize for biases in boundary conditions on four sides of our domain (east, south, west, north). Examination of our prior simulation finds domain-wide biases against either GOSAT or TROPOMI observations that can only be attributed to biased boundary conditions. The optimization is done annually for our main result. In addition, we also perform a seasonal optimization in a sensitivity inversion.

Assuming a Gaussian distribution of error, the optimal estimate of x is obtained by minimizing the cost function (Brasseur and Jacob, 2017; Rodgers, 2000):

(2) J x = x - x A T S A - 1 x - x A + y - F x T S O - 1 ( y - F ( x ) ) ,

where xA represents prior estimates for x and y is the observation vector containing either TROPOMI or GOSAT observations, and F is a function of x representing the forward model. SA and SO are respectively prior and observation error covariance matrices, and their specification is described and discussed in Sect. 3.3.

The forward model (GEOS-Chem) can be described by a linear equation:

(3) F x = K x ,

where K=∇xF is the Jacobian matrix, which describes the sensitivity of observations to the state vector. The cost function is minimized at xJ(x)=0, which yields the optimal estimate (x^),

(4) x ^ = x A + K T S O - 1 K + S A - 1 - 1 K T S O - 1 y - K x A ,

with the posterior error covariance matrix S^,

(5) S ^ = K T S O - 1 K + S A - 1 - 1 ,

and the averaging kernel matrix A that describes the sensitivity of the optimal solution to the true value:

(6) A = x ^ x = I n - S ^ S A - 1 .

The trace of A is referred to as the degrees of freedom for signal (DOFS), which represents the number of independent pieces of information constrained by an observing system.

We apply a transformation vector w to aggregate the posterior estimate regionally (x^r=wTx^). The corresponding posterior error covariance for the region (σ^r2) is then computed as

(7) σ ^ r 2 = w T S ^ w .

3.3 Error specification

The observation error covariance matrix SO represents total random errors from both the methane retrieval (y) and the forward model (F(x)). It can be decomposed as SO=ΣCΣ, where Σ is the diagonal standard deviation matrix and C is the error correlation matrix. In general, inverting SO (which has a dimension of 10 000–10 000 000) in Eqs. (4) and (5) is computationally difficult if C is non-diagonal. The computational challenge can be eased by omitting error correlations (SO=Σ2), but this assumption of error independence unrealistically increases the power of individual observations, leading to overfitting (highly unlikely departure of the posterior solution from the prior estimate) (Zhang et al., 2018). To remedy this issue, previous studies introduce a scalar factor γ (SO=Σ2γ), which serves to enlarge the observation error (γ is usually <1) and thus de-weight individual observations (e.g., Maasakkers et al., 2019; Qu et al., 2021; Zhang et al., 2018). The γ value, which plays the same role as the regularization parameter in Tikhonov methods, can be determined through the graph-based L-curve method (Hansen, 1998; Lu et al., 2021); however, results are sometimes ambiguous and often difficult to interpret physically.

Here, we propose an alternative method. We first determine the diagonal matrix Σ following the residual error method (Heald et al., 2004), which yields observation error standard deviations that average 16 ppbv for TROPOMI and 18 ppbv for GOSAT. Then, we specify a full error correlation matrix C. We parametrize the entry Cij as a function of the distance (Δdij) and the time (Δtij) between ith and jth observations:

(8) C i j = exp - Δ d i j ρ d exp - Δ t i j ρ t ,

where ρd and ρt are correlation scales in space and time, respectively. The values of ρd and ρt can be determined empirically by analyzing spatial and temporal correlations in prior residual errors (Fig. S6). In our case, we find ρt=7 d and ρd=400 km. Finally, we find C̃-1, a computationally tractable (diagonal) approximation to C−1, and replace SO-1 in Eqs. (4) and (5) with Σ-1C̃-1Σ-1. See Appendix A for the derivation of C̃-1. Compared to the traditional γ factor, this method provides better interpretability by explicitly specifying error correlations. Moreover, C̃-1 can be unequivocally determined once C is specified. For comparison, we also include a sensitivity inversion in which SO is specified as Σ2γ with γ=0.6 for GOSAT observations and γ=0.09 for TROPOMI observations following the procedure by Lu et al. (2021) (Fig. S7).

For the prior error covariance matrix SA, we take it as a diagonal matrix and assume a 50 % standard deviation for prior emissions and a 1 % standard deviation for boundary conditions. We also test two alternative configurations for SA in sensitivity inversions: (1) the relative error standard deviation for prior emissions is enlarged to 100 %; and (2) the error standard deviation for prior emissions is specified as 50 % or 1×10-10 kg m−2 s−1, whichever is larger. The latter SA specification gives the inversion more freedom to adjust at locations where prior methane emissions are small or nonexistent.

4 Results and discussion

4.1 Comparison of methane emissions from TROPOMI and GOSAT inversions

Figure 3 shows the correction patterns of methane emissions (posterior–prior emissions) respectively inferred from TROPOMI and GOSAT inversions. Both inversions find that the prior inventory underestimates methane emissions from northeastern China (NEC) and Bangladesh (BAN) and overestimates emissions from central southern China (CSC). The two inversions also find similar correction patterns in central northern China (CNC) with upward adjustments over central Shanxi and downward adjustments over neighboring Henan province. These agreements reflect some consistencies between TROPOMI and GOSAT inversions at the regional level.

TROPOMI and GOSAT inversions show large differences over important source regions, including eastern China (EC) and northern India (IND) (Fig. 3). While the GOSAT inversion suggests that methane emissions over IND should be increased and those from EC decreased relative to prior estimates, the TROPOMI inversion finds the opposite. As a result, regional total methane emissions inferred by the two inversions differ by 7.7 Tg a−1 or 27 % over IND (TROPOMI: 24.6±0.6 Tg a−1, GOSAT: 32.3±0.8 Tg a−1) (errors reported for regional estimates are 1σ standard deviations derived from posterior error covariance matrices) and by 6.4 Tg a−1 or 29 % over EC (TROPOMI: 28.0±0.8 Tg a−1, GOSAT: 21.6±1.0 Tg a−1) (Fig. 3c). In addition, the two inversions also disagree over the northwestern part of the domain (NWD, including parts of Kazakhstan and northern Xinjiang in China, and SXJC, including mainly southern Xinjiang), where TROPOMI indicates large upward adjustments while GOSAT finds agreement with the prior inventory.

Figure 3Spatial distributions of methane emission corrections (posterior–prior) inferred by (a) GOSAT and (b) TROPOMI inversions. Panel (c) shows emissions aggregated by region as defined in blue rectangles in (a) and (b). Error bars represent the standard deviation of regional estimates derived from posterior error covariance matrices (Eq. 7). These errors do not include systematic uncertainties due to inversion setups and are thus optimistic, but they are relevant for comparing results from two identically configured inversions.

Table S2 summarizes methane emission estimates from TROPOMI and GOSAT inversions over the entire East Asia domain and over China. The two inversions find consistent posterior methane emissions from East Asia (TROPOMI: 142.7±1.3 Tg a−1; GOSAT: 142.6±1.5 Tg a−1), with differences in China (TROPOMI: 73.7±0.9 Tg a−1; GOSAT: 66.4±1.1 Tg a−1) largely canceled out by differences in northern India. For China, we attribute 67.9 Tg a−1 for the TROPOMI inversion and 61.6 Tg a−1 for the GOSAT inversion to anthropogenic emissions based on prior sectoral fractions in each spatial cluster. These values are at the high end of previous inversion-based estimates of 43–62 Tg a−1 (Deng et al., 2022; Lu et al., 2021; Miller et al., 2019; Qu et al., 2021; Saunois et al., 2020; Sheng et al., 2021; Stavert et al., 2022; Wang et al., 2021; Zhang et al., 2021, 2022) and are higher than China's latest submission to the United Nations Framework Convention on Climate Change (UNFCCC, 55 Tg a−1) for 2014 (UNFCCC, 2020). These previous inversions mainly used GOSAT observations but differ greatly in their inversion setups (e.g., time, domain coverage, spatial resolution, transport model), thus resulting in a considerable range of estimates. In contrast, the differences in inversions presented in this work are fully due to satellite observations. Our TROPOMI inversion results are consistent with a recent TROPOMI inversion study by Chen et al. (2022), who reported estimates of China's total, anthropogenic, and natural methane emissions of 70.0 (61.6–79.9), 65.0 (57.7–68.4), and 5.0 (3.9–11.6) Tg a−1.

In addition to the main inversion, we also perform a series of sensitivity inversions. The objective is to test whether the comparison between the GOSAT and TROPOMI inversions (e.g., Fig. 3) is affected by the configurations such as the assimilation window and error specifications. There are four sensitivity tests including (1) optimizing emissions seasonally instead of annually, (2) increasing prior error standard deviations from 50 % to 100 %, (3) assigning a minimum prior error standard deviation equivalent to 1×10-10 kg m−2 s−1, and (4) applying a traditional regularization factor following Lu et al. (2021) (γ=0.6 for GOSAT and γ=0.09 for TROPOMI) to account for error correlations in SO instead of the method proposed in Sect. 3.3.

Figures S8 and S9 show that the major findings from our comparison of the GOSAT and TROPOMI inversions in Fig. 3 (agreement in NEC, BAN, and CSC; disagreement in EC, IND, NWD, SXJC) are robust against these perturbed inversion configurations. Consistent with the main inversion, the sensitivity tests find good agreement between the GOSAT and TROPOMI inversions in posterior methane emissions from NEC (upward adjustment), BAN (upward adjustment), and CSC (downward adjustment) but find that discrepancies range 4.1–8.2 Tg a−1 for EC and 5.1–8.8 Tg a−1 for IND (Fig. S9). These results indicate that the effects of inversion configurations are only moderate on systematic differences between the GOSAT and TROPOMI inversions.

4.2 Evaluation of inversion results with independent observations

Both TROPOMI and GOSAT posterior simulations can reduce errors against their respective “training” data relative to the prior simulation (Fig. 4), which is expected for successful inversions. However, concentration fields from the two simulations show varied degrees of agreement across the domain (Fig. 5a). In this section, we use independent high-quality observations to evaluate whether GOSAT and TROPOMI inversion results are consistent, and in the case that they are not, which one is in better agreement with independent data.

Table 1 summarizes performance metrics against these independent observations. Figure S10 plots the time series of these observations in comparison with prior and posterior simulations. GOSAT and TROPOMI inversions perform similarly at background sites such as PDI, UUM, WLG, and LLN. Both posterior simulations considerably reduce biases against in situ observations at WLG and PDI and achieve reasonable agreement at PDI, UUM, and WLG (absolute biases <8 ppbv and R2 between 0.40–0.73). Among these sites, WLG show a relatively low posterior R2 (GOSAT: 0.40; TROPOMI: 0.41) due to inability to capture sub-seasonal variability (Fig. S10). Seasonal optimization done in one of the sensitivity inversions only improves R2 at WLG marginally (Fig. S10). An exception is LLN (a high-mountain background site in the southeast of the domain) where biases grow larger in both posterior simulations (10.8 ppbv for GOSAT and 16.7 ppbv for TROPOMI). This is mainly caused by large seasonal biases in the eastern boundary (Fig. 5c) (see Sect. 4.3.3 for more discussion). The bias is the largest during the monsoon season (May to August) (Fig. S10).

Figure 4Differences in XCH4 between simulations and satellite observations from GOSAT (a, c) and TROPOMI (b, d). Panels (a) and (b) show results for the prior simulation; panel (c) is for the posterior simulation driven by the GOSAT inversion, and panel (d) is for the posterior simulation driven by the TROPOMI inversion. Root mean square errors (RMSEs, in ppbv) and mean biases (MBs, in ppbv, simulation – observation) are inset.

Figure 5Differences in tropospheric methane concentrations (TROPOMI–GOSAT) between GOSAT and TROPOMI posterior simulations. Panel (a) shows the total differences, while panels (b) and (c) decompose the differences to methane emissions and boundary condition bias corrections. The corrections of boundary conditions (in ppbv) by the two inversions are shown.

Table 1Evaluation of simulated methane concentrations against independent observationsa.

a Five sites report surface in situ measurements, with PDI, LLN, UUM, and WLG being continental-scale background sites and AMY a regional site. Two sites (XH and HF) located in eastern China report ground-based total column measurements. The aircraft measurements (CARIBIC) are taken over northern India. b Main inversions are unable to improve the performance for temporal variability, as the optimization of methane emissions is done only annually. Seasonal inversions improve the performance at sites AMY, WLG, XH, LLN (only GOSAT), and HF but in most cases only slightly. Other factors that affect the R2 metric include model transport errors and observation representativeness. c Small prior biases at XH and HF should not be interpreted as evidence for unbiased prior emissions from EC because the prior simulation has substantial low biases in background concentrations as shown by data at WLG (upwind of EC). d Large biases between simulations and observations occur on five days (22 July, 30 September, 3 November, 23 November, and 3 December) at site HF (Fig. S10). Relatively low R2 values are largely affected by these data. Excluding this subset of observations results in correlation coefficients of ∼0.8 for all simulations and mean biases of -2.9±1.4, 5.4±1.3, and 15.9±1.7 ppbv for prior, GOSAT, and TROPOMI simulations, respectively. e Indirect evaluation is performed for CARIBIC data. The value in the “GOSAT” column represents the mean bias between the posterior simulation of the 2010–2017 GOSAT inversion and 2012–2014 CARIBIC aircraft observations. We assume that GOSAT inversions are consistent between years so that the 2012–2014 bias is representative for the 2019 condition. The value in the “TROPOMI” column is computed by subtracting the mean difference along aircraft paths between 2019 GOSAT and TROPOMI posterior simulations (∼4.9 ppbv) (Fig. 5a) from the 2012–2014 GOSAT bias. R2 is not reported for this indirect comparison.

Download Print Version | Download XLSX

On the other hand, methane concentrations from the TROPOMI and GOSAT posterior simulations differ by  10–20 ppbv at sites in methane source regions (i.e., XH and HF within EC and AMY in Korea downwind of EC) (Fig. 5a). Their differences in concentrations are due mainly to higher methane emissions inferred by the TROPOMI inversion than GOSAT over EC (by 6.4 Tg a−1) and Korea (Fig. 3). Our evaluation against in situ measurements at AMY and total column measurements at XH and HF shows consistently high biases of  15–27 ppbv by the TROPOMI posterior simulation and comparatively better agreement (bias ∼8 ppbv) with the GOSAT posterior simulation (Table 1). Smaller mean biases are achieved by the prior simulation at XH and HF (Table 1), but this is largely because of the low background concentration caused by biases in prior boundary conditions (as indicated by the large negative prior bias at the upwind background site WLG; Table 1). The ability to capture temporal variations can be further improved by seasonal optimization of emissions, especially for HF where the influence of the seasonal cycles in rice emissions is strong (Fig. S10). Overall, our results at AMY, XH, and HF support the lower methane emissions from EC inferred by the GOSAT inversion over the TROPOMI inference and indicate that TROPOMI XCH4 retrievals may have regional high biases over EC (more discussion in Sect. 4.3.1).

Methane concentrations from the TROPOMI and GOSAT posterior simulations differ by 4.9 ppbv on average along the CARIBIC flight tracks over the Indo-Gangetic Plain (Fig. 5a). This difference is mainly due to different IND methane emissions between the two inversions (Fig. 5b) with minor contributions from boundary condition bias inference (Fig. 5c). In the absence of concurrent independent observations over IND, we use CARIBIC aircraft observations that are only available from 2012 to 2014 to evaluate the inversions. Since these observations predate TROPOMI, we can only indirectly evaluate by using a simulation driven by methane emissions from a GOSAT inversion for earlier years as an intercomparison platform. We take inversion results from a previous study (Zhang et al., 2022), which performed an East Asia inversion also using GOSAT Proxy XCH4 retrievals. Their inversion is almost identically configured as this study except that it was for 2010–2017. Consistent with our GOSAT results, the GOSAT inversion from Zhang et al. (2022) also found that IND methane emissions should be adjusted upward.

Comparison with these aircraft observations indicates that the 2012–2014 simulation driven by GOSAT-optimized emissions from Zhang et al. (2022) overestimates the aircraft observations by ∼14.9 ppbv (Table 1). On the other hand, the 2019 posterior simulation from the GOSAT inversion is about 4.9 ppbv higher than that from the TROPOMI inversion along flight tracks (Fig. 5a). Assuming that our 2019 GOSAT inversion is consistent with the 2010–2017 GOSAT inversion by Zhang et al. (2022) (mean bias 14.9 ppbv), it thus suggests that the TROPOMI inversion likely agrees better with the CARIBIC observations (mean bias 10.0 ppbv) than the GOSAT inversion. Unlike the EC case, we find relatively small systematic differences in TROPOMI and GOSAT XCH4 retrievals over IND (Fig. 6). Our analysis suggests that good data coverage of TROPOMI over IND is likely responsible for its better performance in constraining methane emissions (see Sect. 4.3.2 for more discussion).

4.3 Attribution of TROPOMI and GOSAT inversion differences

4.3.1 Regional differences in XCH4 retrievals

To understand the cause of differences in the inferred methane emissions, we first compare coincident TROPOMI and GOSAT XCH4 retrievals. The comparison is done following Zhang et al. (2010) where a CTM simulation is used as an intercomparison platform to account for differences in prior profiles and vertical sensitivity between TROPOMI and GOSAT retrievals. TROPOMI XCH4 is on average higher than GOSAT XCH4 over EC by ∼6 ppbv, SXJC by ∼10 ppbv, and NWD by ∼10 ppbv (Fig. 6b), which leads to higher methane emissions inferred by the TROPOMI inversion over these regions (Fig. 3). These differences persist throughout the year in EC and SXJC but appear to be highly seasonal in NWD. The largest TROPOMI–GOSAT differences in NWD ( 30–40 ppbv) occur between December and March. In other regions of interest, the annual averaged TROPOMI–GOSAT XCH4 differences are in general less than 5 ppbv including IND where the two inversions find large discrepancies in posterior methane emissions.

Independent ground-based observations are more consistent with the GOSAT inversion and thus do not support high emissions from EC inferred by the TROPOMI inversion, which indicates that TROPOMI retrievals have systematic regional high biases over EC. In addition, even with enhanced methane emissions in EC, SXJC, and NWD from the TROPOMI inversion, the posterior simulation cannot fully capture these high XCH4 concentrations (Fig. 4d). This is also a hint toward retrieval biases, as it indicates that the inversion finds it difficult to reconcile these high XCH4 patterns with known methane sources and wind information, given our specification of error parameters (SA and SO).

Figure 6Differences in XCH4 between GOSAT and TROPOMI (defined as TROPOMI–GOSAT) shown on the 0.5× 0.625 grid (a) and by region (b). Panel (a) shows annual averages for each grid cell, and panel (b) shows time series of regional averages. Regions are defined in blue rectangles in Fig. 3a.

In addition to EC, large XCH4 differences between GOSAT and TROPOMI are also found in the northwestern part of the domain (SXJC and NWD). Although we do not have independent observations over these regions, we speculate that TROPOMI retrievals have positive biases. SXJC is featured with high surface albedo (desert), while in NWD large TROPOMI and GOSAT differences occur during December and March when surface albedo is low (snow and/or ice cover) (Fig. S11). High and low surface albedo scenes are known to be challenging for the full-physics retrieval. We suggest applying the “blended albedo” filter to TROPOMI observations over these regions before inversion (Chen et al., 2022; Wunch et al., 2011).

In our study, we use the TROPOMI science product from Lorente et al. (2021a), who applied a posterior correction for surface-albedo-dependent biases identified in originally retrieved TROPOMI data. We find that this bias correction scheme does improve the agreement between TROPOMI and GOSAT overall in both their methane column concentrations (Fig. S12) and posterior methane emissions (Fig. S13). However, the agreement is not improved in EC, SXJC, and NWD.

Previous studies have reported decreased accuracy of GOSAT CO2 Proxy retrievals in India owing to errors in the specified CO2 field (Parker et al., 2015; Schepers et al., 2012), which is consistent with a large uncertainty in modeled XCO2 applied to GOSAT CH4  retrievals in India (Fig. S1). The range of modeled XCO2 used in the GOSAT product is equivalent to an XCH4 uncertainty of 0.7 % (∼13 ppbv) in India and Bangladesh. Our result shows that TROPOMI XCH4 is lower than GOSAT XCH4 in the western Indo-Gangetic Plain (around Delhi) and higher in a few locations outside the Indo-Gangetic Plain (Fig. 6a), but the regional difference between the two retrievals is overall small (<5 ppbv) in IND compared to those in EC, SXJC, and NWD (Fig. 6b). Exceptions are November and December when the differences are up to 20 ppbv in IND.

4.3.2 Spatial coverage of observations

Although methane emissions from IND inferred by the GOSAT inversion are considerably larger than those inferred by the TROPOMI inversion, we find relatively small differences in coincident XCH4 retrievals there (Fig. 6), indicating that retrieval biases are unlikely the only cause of discrepancies. Moreover, the two satellite products differ greatly in their data density over the Indo-Gangetic Plain (blue ellipse in Fig. 7) where the discrepancy in inferred methane emissions is the largest. GOSAT has almost no observations over the region, while TROPOMI samples the region fairly well (Fig. S3). We have shown above that indirect comparison with CARIBIC tropospheric aircraft measurements favors lower emissions from IND estimated by the TROPOMI inversion (Table 1). In this section, we explore whether differences in data coverage between TROPOMI and GOSAT may contribute to the discrepancies in inferred emissions.

Figure 7Averaging kernel sensitivities for GOSAT (a) and TROPOMI (b) inversions. Values represent the ability of observations to constrain methane emissions (0: not at all, 1: perfectly). The eastern Indo-Gangetic Plain is marked by blue rectangles. Panel (c) compares the DOFS of regional emissions constrained by TROPOMI and GOSAT inversions.

Figure 7 compares the ability of TROPOMI and GOSAT inversions to constrain the distribution of methane emissions measured by averaging kernel sensitivities (diagonal elements of the averaging kernel matrix). The sum of averaging kernel sensitivities over a region represents the number of pieces of independent information (also known as degrees of freedom for signal, DOFS) constrained by an observation system. Figure 7 shows that the TROPOMI inversion has a larger DOFS value (74) than the GOSAT inversion (46) over the East Asia domain. A large difference in DOFS between the two inversions is found in IND (TROPOMI: 23 vs. GOSAT: 13), indicating a weak observational constraint on emissions from IND by the GOSAT inversion, even though it infers a large emission correction.

We further investigate why this correction is inferred by the GOSAT inversion by examining the contribution of individual observations to the correction. This analysis indicates that the correction is primarily driven by observations in Bangladesh (Fig. 8a). Low XCH4 biases are found over Bangladesh when we compare the prior simulation to either GOSAT or TROPOMI observations (Fig. 4). In the absence of GOSAT observations over the Indo-Gangetic Plain, the inversion partly attributes these XCH4 biases to emissions from the Indo-Gangetic Plain, which is upwind of Bangladesh most of the time, leading to a substantial upward correction of emissions from IND. In contrast, the XCH4 bias over Bangladesh is corrected locally by the TROPOMI inversion. In this case, only small corrections are inferred for emissions from the Indo-Gangetic Plain and the corrections are informed mainly by observations over the Indo-Gangetic Plain (Fig. 8b).

Figure 8Contribution of individual observations to the correction of emissions from the Indo-Gangetic Plain by the (a) GOSAT and (b) TROPOMI inversions. This is done by decomposing the computation of Eq. (4). Results are aggregated on the inversion grid. The scales are different between the two panels.

4.3.3 Regional boundary conditions

Our evaluation against surface observations shows improved agreement at background sites (i.e., PDI, UUM, and WLG) by both inversions (Table 1). This is achieved through simultaneous optimization for biases in boundary conditions together with emissions. As WLG, UUM, and PDI are respectively sensitive to the western, northern, and southern boundaries, this result suggests that satellite observations can correct biases along these boundaries, supporting our inversion configuration. Furthermore, we find that a sensitivity inversion not optimizing for boundary condition biases (S0) cannot reduce large prior biases at WLG and PDI and leads to unrealistically high methane emissions over East Asia (222 Tg a−1) including China (102 Tg a−1).

An exception in Table 1 is LLN (a high-mountain background site in the southeast of the domain) where biases are increased by both inversions. Although the site AMY is also close to the eastern boundary, it has little influence from the southeast monsoon (Fig. 5c). The biases show strong seasonality, with the largest occurring in summer, consistent with ocean-to-land (southeast to northwest) transport by the summer monsoon. Our analysis suggests that this increase in biases is caused by large adjustments at the eastern boundary (GOSAT: 4.5 ppbv; TROPOMI: 25.4 ppbv) rather than changes in methane emissions (Fig. 5). This result indicates that satellite observations that are mainly over land are insufficient to constrain the eastern boundary, which consists mainly of ocean.

We then assess the impact of biases along the eastern boundary on inferred methane emissions. We perform sensitivity inversions using varied levels of fixed (not optimized by the inversion) eastern boundary conditions and find relatively small effects on quantifying annual emissions as expected from prevailing westerlies in midlatitudes. A positive bias of 10 ppbv would result in a reduction of annual methane emissions by 3.3 Tg a−1 (∼2 %) over the East Asia domain, 1.8 Tg a−1 (∼2 %) over China, and 0.75 Tg a−1 (∼3 %) over EC (the most affected region) (Fig. 9). Although the inversion has a weak constraint on the eastern boundary conditions, it does not have a great influence on the posterior emissions.

Figure 9Impact of biases in the eastern boundary on quantification of annual methane emissions. Inversions are performed by using fixed eastern boundary conditions. Sensitivity results are computed from perturbing these fixed eastern boundary conditions by 10 (S10), 20 (S20), and 30 (S30) ppbv.


5 Conclusions

We estimate methane emissions from East Asia for 2019 by applying atmospheric methane column retrievals from two different satellite instruments (GOSAT and TROPOMI) to a high-resolution regional inversion framework, in which methane emissions are optimized on 600 spatial clusters with up to about 0.5 horizontal resolution. Our objective is to assess if consistent methane emissions from East Asia are inferred from inversion of GOSAT and TROPOMI observations. This information adds to the uncertainty characterization of satellite-data-based methane emission quantification.

The two inversions estimate a consistent magnitude of methane emissions from East Asia (TROPOMI: 142.7 Tg a−1; GOSAT: 142.6 Tg a−1) compared to the prior estimate (130 Tg a−1) but differ by ∼10 % in China (TROPOMI: 73.7 Tg a−1; GOSAT: 66.4 Tg a−1). Comparisons at the regional scale show that the GOSAT and TROPOMI inversions find consistent results over central northern China, central southern China, northeastern China, and Bangladesh, where the inferred emissions differ by less than 2.6 Tg a−1. However, the two inversions show large differences over some of the important regions including northern India and eastern China. The inferred methane emissions by GOSAT observations are 7.7 Tg a−1 higher than those by TROPOMI over northern India but 6.4 Tg a−1 lower over eastern China. Large differences in inferred emissions are also found in northwestern China and Kazakhstan (SXJC and NWD). These findings from the comparison of the GOSAT and TROPOMI inversions are robust against varied inversion configurations.

We evaluate the inversion results by comparing GOSAT and TROPOMI posterior simulations with independent observations. We find that independent ground-based in situ observations at AMY and total column observations at XH and HF are more compatible with lower methane emissions from eastern China inferred by the GOSAT inversion than those by the TROPOMI inversion. We also indirectly evaluate against tropospheric aircraft observations over India during 2012–2014 by using a consistent GOSAT inversion of earlier years as an intercomparison platform, which favors lower methane emissions from northern India inferred by the TROPOMI inversion over those by the GOSAT inversion.

The fact that high eastern China emissions inferred from TROPOMI are inconsistent with independent observations suggests high regional biases in TROPOMI retrievals over eastern China. Large retrieval differences between GOSAT and TROPOMI are also found in northwestern China and Kazakhstan, which also leads to substantially higher methane emissions inferred by the TROPOMI inversion. Unfortunately, we do not have independent observations to evaluate the results in these two regions. However, we note that large TROPOMI XCH4 variations in Kazakhstan and northern Xinjiang are coincident with seasonal changes in surface albedo, suggesting possible overcorrection of surface-albedo-dependent biases in TROPOMI retrievals at the regional level.

The two inversions show large discrepancies in emissions over northern India along the Indo-Gangetic Plain, although GOSAT and TROPOMI XCH4 values agree reasonably well. We find that the discrepancy in emissions from the Indo-Gangetic Plain is related to differences in data coverage. In the absence of GOSAT observations over the Indo-Gangetic Plain, the inversion attributes the model–observation differences in XCH4 over Bangladesh partly to its upwind region. In contrast, the TROPOMI inversion finds little emission correction based on the observations over the Indo-Gangetic Plain and attributes the XCH4 differences over Bangladesh primarily to local emissions.

Both inversions show improved agreement at background sites, supporting our optimization of boundary condition biases. An exception is LLN where both inversions show large positive concentration biases against in situ measurements, which results from overcorrections at the eastern boundary by inversions. However, our simulations demonstrate that methane concentration biases at the eastern boundary have relatively small impacts on annual emission inference. The newer version of the TROPOMI methane product includes glint-mode ocean observations, which may improve the optimization of eastern boundary conditions.

Appendix A: Approximation to the inverse of the error correlation matrix

The observational error covariance matrix is decomposed as SO=ΣCΣ, where Σ is a diagonal matrix, while C is generally non-diagonal. The inverse of SO can then be written as SO-1=Σ-1C-1Σ-1. However, the computations of C−1 and SO-1 quickly become intractable as the dimension of the C matrix (m) grows. We therefore seek C̃-1 that approximates C−1 but is easy to compute. To do so, we assume that C̃-1 is a diagonal matrix.

For clarity, we denote C−1 as X and C̃-1 as X̃. We have the following linear system by definition:

(A1) CX = I ,

where I is an identity matrix. To find X is to find its column vectors xi such that

(A2) C x i = e i , i = 1 , 2 , , m .

Here ei=0,,1,,0T is a unit vector with its ith element being 1 and the rest 0.

By assuming that X̃ is diagonal, we impose the condition that its column vectors x̃ispanei. We apply the oblique projection technique to find the solution for x̃i such that the residual vector ei-Cx̃i is orthogonal to the one-dimensional subspace spanned by Cei (Saad, 2003). Hence, we have

(A3) C e i T e i - C x ̃ i = 0 .

Solving the equation yields

(A4) x ̃ i i = C i i | | C e i | | 2 2 = 1 j = 1 m C i j 2 ,

where x̃ii is the ith element of x̃i and ||||2 represents the L-2 norm. Because C is a correlation matrix, its diagonal element Cii is equal to 1.

Consequently, we obtain

(A5) C ̃ - 1 = X ̃ = diag 1 j = 1 m C 1 j 2 , 1 j = 1 m C 2 j 2 , , 1 j = 1 m C m j 2 .

Note that computation of C̃-1 can be readily parallelized for speed-up.

The diagonal elements of C̃-1 can be interpreted as the weight for individual observations. The weight is 1 for an independent observation i uncorrelated with any other observations (Cii=1 and Cij=0 for ij), while the weight can be substantially smaller than 1 for an observation with strong correlation with others (many nonzero Cij terms or large Cij terms for ij).

Data availability

The TROPOMI methane observations are from (Lorente et al., 2021b). The GOSAT methane observations are the University of Leicester GOSAT Proxy XCH4 v9.0 accessible through (Parker and Boesch, 2020). Surface observations at PDI are downloaded from (Nguyen and Steinbacher, 2021). Surface observations at AMY, LLN, UUM, and WLG and aircraft observations from the CARIBIC project are available via the NOAA ObsPack CH4 product (Schuldt et al., 2021). The Xianghe FTIR CH4 data are accessible through (Yang et al., 2019). The Hefei FTIR CH4 from the TCCON network can be accessed by contacting Cheng Liu at the University of Science and Technology of China.


The supplement related to this article is available online at:

Author contributions

RL and YZ designed the study. RL performed the inverse modeling with contributions from YZ, WC, PZ, JL, ZQ, and ZC. RL analyzed and interpreted results with contributions from YZ, CC, HM, GS, ZQ, MZ, RJP, HB, AL, JDM, and IA. RJP and HB provided the GOSAT methane retrievals. AL, JDM, and IA provided the TROPOMI methane retrievals. MZ and PW provided ground-based FTIR methane retrievals at the Xianghe site. RL and YZ wrote the paper with inputs from all authors.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We thank the team that realized TROPOMI and its data products, consisting of the partnership between Airbus Defense and Space Netherlands, KNMI, SRON, and TNO, commissioned by NSO and ESA. Sentinel-5 Precursor is part of the EU Copernicus program; Copernicus (modified) Sentinel-5P data (2018–2020) have been used. 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. The research on GOSAT retrievals used the ALICE High-Performance Computing Facility at the University of Leicester. We thank the TCCON community and Weidong Nan (IAP) for supporting the Xianghe FTIR measurements and the Department of Precision Machinery and Precision Instrumentation, University of Science and Technology of China (Cheng Liu's group), for providing ground-based remote sensing data (Hefei) with their own independent intellectual property rights. We thank the Korea Meteorological Administration, Vietnam Meteorological and Hydrological Administration, China Meteorological Administration, and NOAA for providing surface measurements through GLOBALVIEWplus CH4 ObsPack and WDCGG. We thank the High-Performance Computing Center of Westlake University and the National Supercomputing Center at Wuxi for facility support and technical assistance.

Financial support

This research has been supported by the National Key Research and Development Program of China (grant no. 2021YFB3901000), the National Natural Science Foundation of China (grant no. 42007198), the Environmental Defense Fund, the TROPOMI national program from the NSO, the UK National Centre for Earth Observation funded by the National Environment Research Council (NE/R016518/1 and NE/N018079/1), and the Copernicus Climate Change Service (C3S2_312a_Lot2).

Review statement

This paper was edited by Gabriele Stiller and reviewed by Julia Marshall and two 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. 

Bloom, A. A., Bowman, K. W., Lee, M., Turner, A. J., Schroeder, R., Worden, J. R., Weidner, R., McDonald, K. C., and Jacob, D. J.: A global wetland methane emissions and uncertainty dataset for atmospheric chemical transport models (WetCHARTs version 1.0), Geosci. Model Dev., 10, 2141–2156,, 2017. 

Brasseur, G. P. and Jacob, D. J.: Modeling of Atmospheric Chemistry, Cambridge University Press, Cambridge, USA,, 2017. 

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

Chen, Z., Jacob, D. J., Nesser, H., Sulprizio, M. P., Lorente, A., Varon, D. J., Lu, X., Shen, L., Qu, Z., Penn, E., and Yu, X.: Methane emissions from China: a high-resolution inversion of TROPOMI satellite observations, Atmos. Chem. Phys., 22, 10809–10826,, 2022. 

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. 

Deng, Z., Ciais, P., Tzompa-Sosa, Z. A., Saunois, M., Qiu, C., Tan, C., Sun, T., Ke, P., Cui, Y., Tanaka, K., Lin, X., Thompson, R. L., Tian, H., Yao, Y., Huang, Y., Lauerwald, R., Jain, A. K., Xu, X., Bastos, A., Sitch, S., Palmer, P. I., Lauvaux, T., d'Aspremont, A., Giron, C., Benoit, A., Poulter, B., Chang, J., Petrescu, A. M. R., Davis, S. J., Liu, Z., Grassi, G., Albergel, C., Tubiello, F. N., Perugini, L., Peters, W., and Chevallier, F.: Comparing national greenhouse gas budgets reported in UNFCCC inventories against atmospheric inversions, Earth Syst. Sci. Data, 14, 1639–1675,, 2022. 

Dlugokencky, E. J., Steele, L. P., Lang, P. M., and Masarie, K. A.: The growth rate and distribution of atmospheric methane, J. Geophys. Res., 99, 17021–17043,, 1994. 

Dlugokencky, E. J., Nisbet, E. G., Fisher, R. E., and Lowry, D.: Global atmospheric methane: budget, changes and dangers, Philos. T. Roy. Soc. A, 369, 2058–2072,, 2011. 

Dlugokencky, E. J., Crotwell, A. M., Mund, J. W., Crotwell, M. J., and Thoning, K. W.: Atmospheric Methane Dry Air Mole Fractions from the NOAA GML Carbon Cycle Cooperative Global Air Sampling Network, 1983–2020, Version: 2021-07-30 [data set],, 2021. 

Feng, L., Palmer, P. I., Zhu, S., Parker, R. J., and Liu, Y.: Tropical methane emissions explain large fraction of recent changes in global atmospheric methane growth rate, Nat. Commun., 13, 1378,, 2022. 

Fletcher, S. E. M. and Schaefer, H.: Rising methane: A new climate challenge, Science, 364, 932–933,, 2019. 

Forster, P., Storelvmo, T., Armour, K., Collins, W., Dufresne, J.-L., Frame, D., Lunt, D. J., Palmer, T. M. M. D., Watanabe, M., Wild, M., and Zhang, H.: The Earth's Energy Budget, Climate Feedbacks, and Climate Sensitivity, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK, New York, NY, USA,, 2021. 

Frankenberg, C., Meirink, J. F., van Weele, M., Platt, U., and Wagner, T.: Assessing Methane Emissions from Global Space-Borne Observations, Science, 308, 1010–1014,, 2005. 

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., 111, D07303,, 2006. 

Fung, I. Y., John, J. G., Lerner, J., Matthews, E., Prather, M. J., Steele, L. P., and Fraser, P. J.: Three dimensional model synthesis of the global methane cycle, J. Geophys. Res., 96, 13033–13065,, 1991. 

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. 

Ganesan, A. L., Schwietzke, S., Poulter, B. I., Arnold, T., Lan, X., Rigby, M. L., Vogel, F. R., Werf, G. R., Janssens-Maenhout, G., Boesch, H., Pandey, S., Manning, A. J., Jackson, R. B., Nisbet, E. G., and Manning, M. R.: Advancing Scientific Understanding of the Global Methane Budget in Support of the Paris Agreement, Global Biogeochem. Cy., 33, 1475–1512,, 2019. 

Gao, J., Guan, C., Zhang, B., and Li, K.: Decreasing methane emissions from China's coal mining with rebounded coal production, Environ. Res. Lett., 16, 124037,, 2021. 

Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., da Silva, A. M., Gu, W., Kim, G.-K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. 

Hansen, P. C.: Rank-deficient and discrete ill-posed problems: numerical aspects of linear inversion, Society for Industrial and Applied Mathematics (SIAM), Philadelphia,, 1998. 

Heald, C. L., Jacob, D. J., Jones, D. B. A., Palmer, P. I., Logan, J. A., Streets, D. G., Sachse, G. W., Gille, J., Hoffman, R. N., and Nehrkorn, T.: Comparative inverse analysis of satellite (MOPITT) and aircraft (TRACE-P) observations to estimate Asian sources of carbon monoxide, J. Geophys. Res., 109, D23306,, 2004. 

Hu, H., Hasekamp, O., Butz, A., Galli, A., Landgraf, J., Aan de Brugh, J., Borsdorff, T., Scheepmaker, R., and Aben, I.: The operational methane retrieval algorithm for TROPOMI, Atmos. Meas. Tech., 9, 5423–5440,, 2016. 

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

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. 

Jacob, D. J., Varon, D. J., Cusworth, D. H., Dennison, P. E., Frankenberg, C., Gautam, R., Guanter, L., Kelley, J., McKeever, J., Ott, L. E., Poulter, B., Qu, Z., Thorpe, A. K., Worden, J. R., and Duren, R. M.: Quantifying methane emissions from the global scale down to point sources using satellite observations of atmospheric methane, Atmos. Chem. Phys., 22, 9617–9646,, 2022. 

Janssens-Maenhout, G., Crippa, M., Guizzardi, D., Muntean, M., Schaaf, E., Dentener, F., Bergamaschi, P., Pagliari, V., Olivier, J. G. J., Peters, J. A. H. W., van Aardenne, J. A., Monni, S., Doering, U., Petrescu, A. M. R., Solazzo, E., and Oreggioni, G. D.: EDGAR v4.3.2 Global Atlas of the three major greenhouse gas emissions for the period 1970–2012, Earth Syst. Sci. Data, 11, 959–1002,, 2019. 

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. Optics, 48, 6716–6733,, 2009. 

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. 

Lee, H., Han, S.-O., Ryoo, S.-B., Lee, J.-S., and Lee, G.-W.: The measurement of atmospheric CO2 at KMA GAW regional stations, its characteristics, and comparisons with other East Asian sites, Atmos. Chem. Phys., 19, 2149–2163,, 2019. 

Liu, C., Sun, Y., Shan, C., Wang, W., Notholt, J., Palm, M., Yin, H., Tian, Y., Gao, J., and Mao, H.: Long-Term Observations of Atmospheric Constituents at the First Ground-Based High-Resolution Fourier-Transform Spectrometry Observation Station in China, Engineering, 22, 201–214,, 2023. 

Liu, G., Peng, S., Lin, X., Ciais, P., Li, X., Xi, Y., Lu, Z., Chang, J., Saunois, M., Wu, Y., Patra, P. K., Chandra, N., Zeng, H., and Piao, S.: Recent Slowdown of Anthropogenic Methane Emissions in China Driven by Stabilized Coal Production, Environ. Sci. Technol. Lett., 8, 739–746,, 2021. 

Lorente, A., Borsdorff, T., Butz, A., Hasekamp, O., aan de Brugh, J., Schneider, A., Wu, L., Hase, F., Kivi, R., Wunch, D., Pollard, D. F., Shiomi, K., Deutscher, N. M., Velazco, V. A., Roehl, C. M., Wennberg, P. O., Warneke, T., and Landgraf, J.: Methane retrieved from TROPOMI: improvement of the data product and validation of the first 2 years of measurements, Atmos. Meas. Tech., 14, 665–684,, 2021a. 

Lorente, A., Borsdorff, T., aan de Brugh, J., Landgraf, J., and Hasekamp, O.: SRON S5P – RemoTeC scientific TROPOMI XCH4 [data set], (last access: 29 December 2021), 2021b. 

Lorente, A., Borsdorff, T., Landgraf, J., and the SRON L2 team: SRON RemoTeC-S5P scientific XCH4 data product Product User Guide, (last access: 1 December 2022), 2022. 

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., Boesch, 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., 21, 4637–4657,, 2021. 

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. 

McNorton, J., Bousserez, N., Agustí-Panareda, A., Balsamo, G., Cantarello, L., Engelen, R., Huijnen, V., Inness, A., Kipling, Z., Parrington, M., and Ribas, R.: Quantification of methane emissions from hotspots and during COVID-19 using a global atmospheric inversion, Atmos. Chem. Phys., 22, 5961–5981,, 2022. 

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

Monteil, G., Houweling, S., Butz, A., Guerlet, S., Schepers, D., Hasekamp, O. P., Frankenberg, C., Scheepmaker, R. A., Aben, I., and Röckmann, T.: Comparison of CH4 inversions based on 15 months of GOSAT and SCIAMACHY observations, J. Geophys. Res.-Atmos., 118, 11807–11823,, 2013. 

Murguia-Flores, F., Arndt, S., Ganesan, A. L., Murray-Tortarolo, G., and Hornibrook, E. R. C.: Soil Methanotrophy Model (MeMo v1.0): a process-based model to quantify global uptake of atmospheric methane by soil, Geosci. Model Dev., 11, 2009–2032,, 2018. 

Nguyen, N. A. and Steinbacher, M.: Atmospheric CH4 at Pha Din by Viet Nam Meteorological and Hydrological Administration, dataset published as CH4_PDI_surface-insitu_VNMHA_data1 at WDCGG, ver. 2020-07-13-2137 [data set],, 2021. 

Nisbet, E. G., Manning, M. R., Dlugokencky, E. J., Fisher, R. E., Lowry, D., Michel, S. E., Myhre, C. L., Platt, S. M., Allen, G., Bousquet, P., Brownlow, R., Cain, M., France, J. L., Hermansen, O., Hossaini, R., Jones, A. E., Levin, I., Manning, A. C., Myhre, G., Pyle, J. A., Vaughn, B. H., Warwick, N. J., and White, J. W. C.: Very Strong Atmospheric Methane Growth in the 4 Years 2014–2017: Implications for the Paris Agreement, Global Biogeochem. Cy., 33, 318–342,, 2019. 

Pandey, S., Houweling, S., Krol, M., Aben, I., Chevallier, F., Dlugokencky, E. J., Gatti, L. V., Gloor, E., Miller, J. B., Detmers, R., Machida, T., and Röckmann, T.: Inverse modeling of GOSAT-retrieved ratios of total column CH4 and CO2 for 2009 and 2010, Atmos. Chem. Phys., 16, 5043–5062,, 2016. 

Parker, R. J. and Boesch, H.: University of Leicester GOSAT Proxy XCH4 v9.0, Centre for Environmental Data Analysis [data set],, 2020. 

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. 

Parker, R. J., Webb, A., Boesch, H., Somkuti, P., Barrio Guillo, R., Di Noia, A., Kalaitzi, N., Anand, J. S., 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., Té, Y., Velazco, V. A., Warneke, T., Wennberg, P. O., and Wunch, D.: A decade of GOSAT Proxy satellite CH4 observations, Earth Syst. Sci. Data, 12, 3383–3412,, 2020. 

Qu, Z., Jacob, D. J., Shen, L., Lu, X., Zhang, Y., Scarpelli, T. R., Nesser, H., Sulprizio, M. P., Maasakkers, J. D., Bloom, A. A., Worden, J. R., Parker, R. J., and Delgado, A. L.: Global distribution of methane emissions: a comparative inverse analysis of observations from the TROPOMI and GOSAT satellite instruments, Atmos. Chem. Phys., 21, 14159–14175,, 2021. 

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., Mühle, J., and Porter, L. W.: Renewed growth of atmospheric methane, Geophys. Res. Lett., 35, L22805,, 2008. 

Rodgers, C. D.: Inverse methods for atmospheric sounding: theory and practice, World Scientific, Singapore, ISBN 978-981-02-2740-1, 2000. 

Saad, Y.: Iterative methods for sparse linear systems, 2nd ed., Society for Industrial and Applied Mathematics (SIAM), Philadelphia,, 2003. 

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. 

Scarpelli, T. R., Jacob, D. J., Maasakkers, J. D., Sulprizio, M. P., Sheng, J.-X., Rose, K., Romeo, L., Worden, J. R., and Janssens-Maenhout, G.: A global gridded (0.1×0.1) inventory of methane emissions from oil, gas, and coal exploitation based on national reports to the United Nations Framework Convention on Climate Change, Earth Syst. Sci. Data, 12, 563–575,, 2020. 

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. 

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. 

Schuldt, K. N., Aalto, T., Andrews, A., Aoki, S., Arduini, J., Baier, B., Bergamaschi, P., Biermann, T., Biraud, S. C., Boenisch, H., Brailsford, G., Chen, H., Colomb, A., Conil, S., Cristofanelli, P., Cuevas, E., Daube, B., Davis, K., Mazière, M. D., Delmotte, M., Desai, A., DiGangi, J. P., Dlugokencky, E., Elkins, J. W., Emmenegger, L., Fischer, M. L., Gatti, L. V., Gehrlein, T., Gerbig, C., Gloor, E., Goto, D., Haszpra, L., Hatakka, J., Heimann, M., Heliasz, M., Hermanssen, O., Hintsa, E., Holst, J., Ivakhov, V., Jaffe, D., Joubert, W., Kang, H.-Y., Karion, A., Kazan, V., Keronen, P., Ko, M.-Y., Kominkova, K., Kort, E., Kozlova, E., Krummel, P., Kubistin, D., Labuschagne, C., Langenfelds, R., Laurent, O., Laurila, T., Lauvaux, T., Lee, J., Lee, H., Lee, C.-H., Lehner, I., Leppert, R., Leuenberger, M., Lindauer, M., Loh, Z., Lopez, M., Machida, T., Mammarella, I., Manca, G., Marek, M. V., Martin, M. Y., Matsueda, H., McKain, K., Miles, N., Miller, C. E., Miller, J. B., Moore, F., Morimoto, S., Munro, D., Myhre, C. L., Mölder, M., Müller-Williams, J., Nichol, S., Niwa, Y., O’Doherty, S., Obersteiner, F., Piacentino, S., Pichon, J. M., Pittman, J., Plass-Duelmer, C., Ramonet, M., Richardson, S., Rivas, P. P., Saito, K., Santoni, G., Sasakawa, M., Scheeren, B., Schuck, T., Schumacher, M., and Seifert, T.: Multi-laboratory compilation of atmospheric methane data for the period 1983–2020; obspack_ch4_ 1_ GLOBALVIEWplus_v4.0_2021-10-14 [data set],, 2021. 

Sha, M. K., Langerock, B., Blavier, J.-F. L., Blumenstock, T., Borsdorff, T., Buschmann, M., Dehn, A., De Mazière, M., Deutscher, N. M., Feist, D. G., García, O. E., Griffith, D. W. T., Grutter, M., Hannigan, J. W., Hase, F., Heikkinen, P., Hermans, C., Iraci, L. T., Jeseck, P., Jones, N., Kivi, R., Kumps, N., Landgraf, J., Lorente, A., Mahieu, E., Makarova, M. V., Mellqvist, J., Metzger, J.-M., Morino, I., Nagahama, T., Notholt, J., Ohyama, H., Ortega, I., Palm, M., Petri, C., Pollard, D. F., Rettinger, M., Robinson, J., Roche, S., Roehl, C. M., Röhling, A. N., Rousogenous, C., Schneider, M., Shiomi, K., Smale, D., Stremme, W., Strong, K., Sussmann, R., Té, Y., Uchino, O., Velazco, V. A., Vigouroux, C., Vrekoussis, M., Wang, P., Warneke, T., Wizenberg, T., Wunch, D., Yamanouchi, S., Yang, Y., and Zhou, M.: Validation of methane and carbon monoxide from Sentinel-5 Precursor using TCCON and NDACC-IRWG stations, Atmos. Meas. Tech., 14, 6249–6304,, 2021. 

Shen, L., Zavala-Araiza, D., Gautam, R., Omara, M., Scarpelli, T. R., Sheng, J., Sulprizio, M. P., Zhuang, J., Zhang, Y., Qu, Z., Lu, X., Hamburg, S. P., and Jacob, D. J.: Unravelling a large methane emission discrepancy in Mexico using satellite observations, Remote Sens. Environ., 260, 112461,, 2021. 

Shen, L., Gautam, R., Omara, M., Zavala-Araiza, D., Maasakkers, J. D., Scarpelli, T. R., Lorente, A., Lyon, D., Sheng, J., Varon, D. J., Nesser, H., Qu, Z., Lu, X., Sulprizio, M. P., Hamburg, S. P., and Jacob, D. J.: Satellite quantification of oil and natural gas methane emissions in the US and Canada including contributions from individual basins, Atmos. Chem. Phys., 22, 11203–11215,, 2022. 

Sheng, J., Song, S., Zhang, Y., Prinn, R. G., and Janssens-Maenhout, G.: Bottom-Up Estimates of Coal Mine Methane Emissions in China: A Gridded Inventory, Emission Factors, and Trends, Environ. Sci. Tech. Lett., 6, 473–478,, 2019. 

Sheng, J., Tunnicliffe, R. L., Ganesan, A. L., Maasakkers, J. D., Shen, L., Prinn, R. G., Song, S., Zhang, Y., Scarpelli, T. R., Anthony Bloom, A., Rigby, M. L., Manning, A. J., Parker, R. J., Boesch, H., Lan, X., Zhang, B., Zhuang, M., and Lu, X.: Sustained methane emissions from China after 2012 despite declining coal production and rice-cultivated area, Environ. Res. Lett., 16, 104018,, 2021. 

Stavert, A. R., Saunois, M., Canadell, J. G., Poulter, B., Jackson, R. B., Regnier, P., Lauerwald, R., Raymond, P. A., Allen, G. H., Patra, P. K., Bergamaschi, P., Bousquet, P., Chandra, N., Ciais, P., Gustafson, A., Ishizawa, M., Ito, A., Kleinen, T., Maksyutov, S., McNorton, J., Melton, J. R., Müller, J., Niwa, Y., Peng, S., Riley, W. J., Segers, A., Tian, H., Tsuruta, A., Yin, Y., Zhang, Z., Zheng, B., and Zhuang, Q.: Regional trends and drivers of the global methane budget, Glob. Change Biol., 28, 182–200,, 2022. 

Turner, A. J. and Jacob, D. J.: Balancing aggregation and smoothing errors in inverse models, Atmos. Chem. Phys., 15, 7039–7048,, 2015. 

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. 

UNFCCC: Greenhouse Gas Inventory Data [data set], by_ party (last access: 30 May 2022), 2020. 

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

Wang, F., Maksyutov, S., Janardanan, R., Tsuruta, A., Ito, A., Morino, I., Yoshida, Y., Tohjima, Y., Kaiser, J. W., Janssens-Maenhout, G., Lan, X., Mammarella, I., Lavric, J. V., and Matsunaga, T.: Interannual variability on methane emissions in monsoon Asia derived from GOSAT and surface observations, Environ. Res. Lett., 16, 024040,, 2021. 

Wecht, K. J., Jacob, D. J., Frankenberg, C., Jiang, Z., and Blake, D. R.: Mapping of North American methane emissions with high spatial resolution by inversion of SCIAMACHY satellite data, J. Geophys. Res.-Atmos., 119, 7741–7756,, 2014. 

WMO: WMO Greenhouse Gas Bulletin – The State of Greenhouse Gases in the Atmosphere Based on Global Observations through 2020, (last access: 6 July 2022), 2021. 

Wunch, D., Wennberg, P. O., Toon, G. C., Connor, B. J., Fisher, B., Osterman, G. B., Frankenberg, C., Mandrake, L., O'Dell, C., Ahonen, P., Biraud, S. C., Castano, R., Cressie, N., Crisp, D., Deutscher, N. M., Eldering, A., Fisher, M. L., Griffith, D. W. T., Gunson, M., Heikkinen, P., Keppel-Aleks, G., Kyrö, E., Lindenmaier, R., Macatangay, R., Mendonca, J., Messerschmidt, J., Miller, C. E., Morino, I., Notholt, J., Oyafuso, F. A., Rettinger, M., Robinson, J., Roehl, C. M., Salawitch, R. J., Sherlock, V., Strong, K., Sussmann, R., Tanaka, T., Thompson, D. R., Uchino, O., Warneke, T., and Wofsy, S. C.: A method for evaluating bias in global measurements of CO2 total columns from space, Atmos. Chem. Phys., 11, 12317–12337,, 2011. 

Yang, Y., Zhou, M., Langerock, B., Sha, M. K., Hermans, C., Wang, T., Ji, D., Vigouroux, C., Kumps, N., Wang, G., De Mazière, M., and Wang, P.: Ground-based FTIR CO2, CH4 and CO measurements at Xianghe, China, Royal Belgian Institute for Space Aeronomy [data set],, 2019. 

Yang, Y., Zhou, M., Langerock, B., Sha, M. K., Hermans, C., Wang, T., Ji, D., Vigouroux, C., Kumps, N., Wang, G., De Mazière, M., and Wang, P.: New ground-based Fourier-transform near-infrared solar absorption measurements of XCO2, XCH4 and XCO at Xianghe, China, Earth Syst. Sci. Data, 12, 1679–1696,, 2020. 

Yin, Y., Chevallier, F., Ciais, P., Bousquet, P., Saunois, M., Zheng, B., Worden, J., Bloom, A. A., Parker, R. J., Jacob, D. J., 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., 21, 12631–12647,, 2021. 

Yokota, T., Yoshida, Y., Eguchi, N., Ota, Y., Tanaka, T., Watanabe, H., and Maksyutov, S.: Global Concentrations of CO2 and CH4 Retrieved from GOSAT: First Preliminary Results, Sola, 5, 160–163,, 2009. 

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.  

Zhang, L., Jacob, D. J., Liu, X., Logan, J. A., Chance, K., Eldering, A., and Bojkov, B. R.: Intercomparison methods for satellite measurements of atmospheric composition: application to tropospheric ozone from TES and OMI, Atmos. Chem. Phys., 10, 4725–4739,, 2010. 

Zhang, Y., Jacob, D. J., Maasakkers, J. D., Sulprizio, M. P., Sheng, J.-X., Gautam, R., and Worden, J.: Monitoring global tropospheric OH concentrations using satellite observations of atmospheric methane, Atmos. Chem. Phys., 18, 15959–15973,, 2018. 

Zhang, Y., Gautam, R., Pandey, S., Omara, M., Maasakkers, J. D., Sadavarte, P., Lyon, D. R., Nesser, H., Sulprizio, M. P., Varon, D. J., Zhang, R., Houweling, S., Zavala-Araiza, D., Alvarez, R. A., Lorente, A., Hamburg, S. P., Aben, I., and Jacob, D. J.: Quantifying methane emissions from the largest oil-producing basin in the United States from space, Sci. Adv., 6, eaaz5120,, 2020. 

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 2010–2018 by inverse analysis of GOSAT observations, Atmos. Chem. Phys., 21, 3643–3666,, 2021. 

Zhang, Y., Fang, S., Chen, J., Lin, Y., Chen, Y., Liang, R., Jiang, K., Parker, R. J., Boesch, H., Steinbacher, M., Sheng, J.-X., Lu, X., Song, S., and Peng, S.: Observed changes in China's methane emissions linked to policy drivers, P. Natl. Acad. Sci. USA, 119, e2202742119,, 2022. 

Short summary
We compare and evaluate East Asian methane emissions inferred from different satellite observations (GOSAT and TROPOMI). The results show discrepancies over northern India and eastern China. Independent ground-based observations are more consistent with TROPOMI-derived emissions in northern India and GOSAT-derived emissions in eastern China.
Final-revised paper