Articles | Volume 20, issue 21
Research article
07 Nov 2020
Research article |  | 07 Nov 2020

Quantifying sources of Brazil's CH4 emissions between 2010 and 2018 from satellite data

Rachel L. Tunnicliffe, Anita L. Ganesan, Robert J. Parker, Hartmut Boesch, Nicola Gedney, Benjamin Poulter, Zhen Zhang, Jošt V. Lavrič, David Walter, Matthew Rigby, Stephan Henne, Dickon Young, and Simon O'Doherty

Brazil's CH4 emissions over the period 2010–2018 were derived for the three main sectors of activity: anthropogenic, wetland and biomass burning. Our inverse modelling estimates were derived from GOSAT (Greenhouse gases Observing SATellite) satellite measurements of XCH4 combined with surface data from Ragged Point, Barbados, and the high-resolution regional atmospheric transport model NAME (Numerical Atmospheric-dispersion Modelling Environment). We find that Brazil's mean emissions over 2010–2018 are 33.6±3.6Tgyr-1, which are comprised of 19.0±2.6Tgyr-1 from anthropogenic (primarily related to agriculture and waste), 13.0±1.9Tgyr-1 from wetlands and 1.7±0.3Tgyr-1 from biomass burning sources. In addition, between the 20112013 and 2014–2018 periods, Brazil's mean emissions rose by 6.9±5.3Tgyr-1 and this increase may have contributed to the accelerated global methane growth rate observed during the latter period. We find that wetland emissions from the western Amazon increased during the start of the 2015–2016 El Niño by 3.7±2.7Tgyr-1 and this is likely driven by increased surface temperatures. We also find that our estimates of anthropogenic emissions are consistent with those reported by Brazil to the United Framework Convention on Climate Change. We show that satellite data are beneficial for constraining national-scale CH4 emissions, and, through a series of sensitivity studies and validation experiments using data not assimilated in the inversion, we demonstrate that (a) calibrated ground-based data are important to include alongside satellite data in a regional inversion and that (b) inversions must account for any offsets between the two data streams and their representations by models.

1 Introduction

Methane (CH4) is the second most important anthropogenic greenhouse gas behind carbon dioxide due to its radiative properties and atmospheric abundance (IPCC2013). After a brief plateau period around the turn of the century (Cunnold2002; Dlugokencky et al.2003), CH4 mole fractions began rising again globally after 2007 (Rigby et al.2008; Dlugokencky et al.2009; Frankenberg et al.2011; Nisbet et al.2016) with some of the strongest growth rates occurring from 2014 onward (Nisbet et al.2019). This increase in CH4 growth rate was accompanied by a shift in the δ13C-CH4 isotopic ratios to more negative values, suggesting a change in the global makeup of sources and/or sinks. The drivers responsible for this shift are presently not well understood, and proposals include increases from tropical wetlands or agriculture, decreases in biomass burning, and changes in fossil fuel emissions or in the hydroxyl radical sink (e.g. Monteil et al.2011; Schaefer et al.2016; Schwietzke et al.2016; Nisbet et al.2016, 2019; Rigby et al.2017; Worden et al.2017; McNorton et al.2018; Turner et al.2019). Quantifying the CH4 budget and understanding how major sources and sinks have evolved is key to designing emission pathways that limit global warming due to the importance of CH4 in meeting global climate targets (Ganesan et al.2019; Nisbet et al.2019, 2020).

The Paris Agreement pledges to limit warming to less than 2 C with an aspiration for less than 1.5 C warming from pre-industrial levels (UNFCCC2015). The mitigation action taken by each country is dependent on their own Nationally Determined Contributions and accounting for national emissions will occur through inventory or “bottom-up” methods. To assess whether these self-determined targets are being met, independent estimates can be derived using “top-down” strategies that use atmospheric measurements to quantify sector-level emissions estimates at near real time and at high resolution (e.g. Ganesan et al.2019). Using both top-down and bottom-up methods together for national-scale greenhouse gas estimation is considered to be best practice (Calvo Buendia et al.2019) and allows for the greatest process-level understanding of changes in the atmosphere.

Brazil is thought to be a major contributor to global CH4 emissions due to its variety of natural and human-made sources. Anthropogenic emissions arise from agriculture, waste and biomass burning (Ministry of Foreign Affairs et al.2019). Brazil's 2018 Biennial Update Report to the United Framework Convention on Climate Change (UNFCCC) states that 17.6±4.0 Tg of CH4 was emitted from anthropogenic sources in 2015. The majority of these emissions were from agricultural processes (70 % from enteric fermentation, manure management and crop residue burning) with the remainder coming from waste (16 %), energy (4 %) and land-use change (6 %) (Ministry of Foreign Affairs et al.2019).

Around 60 % of the Amazon basin and 80 % of the Pantanal wetland region (Ministry of Science and Innovation2016; Schulz et al.2019) exist within Brazil in the northern and central-western regions of the country, respectively. The primary areas of agricultural activity are in central and southern provinces and include cattle ranching and sugar cane production, while waste and fossil fuel emissions are focused in population centres along the eastern coast (Ministry of Foreign Affairs et al.2019). Biomass burning occurs along the “arc of deforestation” along the southern edge of the Amazon rainforest during and after the dry season (July–October). This is in contrast to Amazon wetland emissions which peak during and after the wet season (December–March).

Current top-down estimates of CH4 emissions from Brazil, the Amazon and tropical South America vary depending on the method, source of data and area considered. In the synthesis of Saunois et al. (2016), across the Tropical South America region, emissions estimates derived using different datasets and top-down methods span the large range of 63119 Tg yr−1 (2369 Tg yr−1 from wetlands) for 2012. Across the Amazon basin, estimates of total emissions derived from aircraft measurements are between ∼16 and 72 Tg yr−1 derived for May 2009 (Beck et al.2013) and 3143 Tg yr−1 for 20102013 (Wilson et al.2016; Pangala et al.2017). A recent study that performed a regional analysis using satellite data by Janardanan et al. (2020) found Brazil's emissions alone, on average, to be 56.2 Tg yr−1 (39.8±12.4Tgyr-1 from wetlands) across 20112017. In addition, many previous studies have estimated emissions globally using satellite data (e.g. Bergamaschi et al.2009; Feng et al.2017). The wide range of estimates indicate that large uncertainties exist and these uncertainties are exacerbated when estimating emissions over smaller scales such as the Amazon basin or when quantifying individual sources.

Through use of a high-resolution regional inversion framework coupled with satellite measurements of CH4, we inferred spatial and temporal distributions of Brazil's CH4 emissions from 2010 to 2018. The regional inversion approach provides the benefit that uncertainties in the hydroxyl radical CH4 sink (Rigby et al.2017; Turner et al.2017; Nguyen et al.2020), a limitation in global approaches, can be neglected. Owing to a spatial and temporal difference in Brazil's major CH4 sources, these emissions are further partitioned into source sectors (Sect. 3.1) and are also presented for different wetland regions (Sect. 3.2). We demonstrate the importance of the inversion setup when using satellite data to estimate country- and basin-scale emissions. Independent validation using in situ data is shown in Sect. 3.3, and sensitivity studies, testing a range of different input factors, are discussed in Sect. 3.4.

2 Methods

2.1CH4 measurements

We used data from three sources: (1) the University of Leicester v7.2 total column CH4 product from the Thermal And Near-infrared Sensor for carbon Observation Fourier Transform Spectrometer (TANSO-FTS) instrument on board the Greenhouse gases Observing SATellite (GOSAT) from April 2010 to November 2018, (2) surface data from Ragged Point, Barbados (RPB; coordinates: 13.17 N, 59.43 W), and (3) surface data from the Amazon Tall Tower Observatory (ATTO; coordinates: 2.15 S, 59.01 W) for external validation of the inversion. Figure 1 shows the positions of GOSAT points for a 1-year period and the locations of the two surface stations.

Figure 1GOSAT measurements grouped into February–April (FMA), May–July (MJJ), August–October (ASO) and November–January (NDJ). Measurements are shown here for the February 2014–January 2015 period and averaged in 1 bins for clarity. FMA and ASO are representative of the wet and dry seasons over the Amazon, respectively. Positions for the ATTO and RPB sites are also shown. The Brazilian border is highlighted in blue.

Dry-air column-averaged CH4 mole fractions (XCH4) were derived using the CO2 proxy method, which multiplies the XCH4∕XCO2 ratio by a model XCO2 field (Parker et al.2011, 2015). The model XCO2 is based on the median of three global models which all assimilated surface site measurement data: GEOS-Chem (Feng et al.2011), Carbon Tracker (Peters et al.2007) and LMDZ (MACC/CAMS) (Chevallier et al.2010). This GOSAT product was previously compared to aircraft measurements over the Amazon basin by extrapolating the aircraft profiles through the troposphere and using a stratospheric model. Differences ranged from −1.9 to 9.7 nmol mol−1 (Webb et al.2016).

We used level 2 GOSAT measurements that were taken in nadir mode within an area that extended from 35.8 S to 7.3 N and from 76.0 to 32.8 W and that passed the quality threshold. We only used nadir measurements to minimise the effect of any unquantified biases between nadir and glint mode. In addition, we filtered data where the surface pressure deviated from the retrieval grid by enough to reduce the number of retrieval levels to less than 20. The remaining data points were averaged across a 0.23×0.35 grid to match the lowest resolution of the atmospheric transport model grid cell (see Sect. 2.3) across the 2010–2018 time period leaving ∼1300 data points on average per month.

Data from RPB were used alongside the satellite measurements in the inversion to provide additional constraints on the boundary conditions. RPB is part of the Advanced Global Atmospheric Gases Experiment (AGAGE) network (Prinn et al.2018b) and predominantly measures well-mixed background air. Measurements up to 2017 were made using GC-FID (gas chromatograph with a flame ionisation detector) and beyond this with a CRDS (cavity ring-down spectrometer). All data were averaged into hourly samples.

Measurements from ATTO (Andreae et al.2015; Botía B. et al.2020) were used for external validation of the inversion results. ATTO is located near Manaus within the Amazon rainforest. The position and predominant north-easterly wind direction means that this site is particularly sensitive to CH4 emitted from wetlands but may also receive air masses from regions of biomass burning and other human activity (Andreae et al.2015; Pöhlker et al.2019). CH4 mole fractions from 2014 to 2018 derived from CRDS instrumentation have been used in this study. Hourly-mean measurements from the highest inlet on the tower, at 79 m, were used as they are assumed to be the most representative of regional air masses.

2.2 Atmospheric transport model

To provide the relationship between atmospheric mole fractions at a receptor and a surface emissions field, we used the high-resolution Lagrangian atmospheric transport model NAME (Numerical Atmospheric-dispersion Modelling Environment) (Jones et al.2007). Model particles were released for each GOSAT and surface measurement time and location and tracked backward in time for 30 d. The model tracked the interaction of these particles with the surface (defined as 0–40 m above modelled ground level) to quantify the sensitivity to regional emissions. The times and locations that particles left the model domain were recorded to quantify the sensitivity to boundary conditions. NAME was driven by meteorological inputs from the Unified Model (UM) spanning resolutions between 0.23 to 0.09 latitude and 0.35 to 0.14 longitude over the 2010–2018 period. The annual mean sensitivity to GOSAT measurements used in this study are shown in Fig. A1 in the appendix.

Satellite measurements require footprints of the total atmospheric column and model particles were released at multiple heights based on the pressure levels defined within the GOSAT product (see Ganesan et al.2017, for a description of how NAME was used to simulate XCH4 by applying averaging kernels, pressure weights and a priori information for satellite data). The main modification in the NAME setup from Ganesan et al. (2017) made here is that surface pressure in GOSAT was corrected to match the surface pressure from the UM. Occasionally, the corrected surface pressure level was lower than the first model level, and in these cases, the retrievals were discarded. This ensured consistency between the model defining the GOSAT pressure levels and NAME.

2.3 Inversion method

Top-down emissions estimates were inferred using a hierarchical Bayesian trans-dimensional inversion using reversible jump Markov chain Monte Carlo (MCMC). A full description of the method can be found in Ganesan et al. (2014) and Lunt et al. (2016). The hierarchical component employs a set of hyperparameters that define the model–measurement and prior emissions uncertainties, which were explored as part of the inversion. Inclusion of these additional model parameters allows for uncertainties in the system to be more accurately captured. The trans-dimensional component of the inversion allowed for the spatial inversion grid to be estimated as part of the inversion rather than being defined a priori.

The a priori inputs to the inversion are described in Sect. 2.4. The emissions PDF (probability density function) was defined as log-normal to prevent non-physical negative solutions from being reached. The standard deviation of this PDF was allowed to vary between 0.05 and 20.0 (with a value of one being equivalent to the prior emissions magnitude). The model–measurement uncertainty was governed by a Gaussian distribution centred on zero with a standard deviation that was a hyper-parameter in the inversion. The standard deviation hyper-parameter was described by a uniform distribution with a range of 0.2 to 200 nmol mol−1.

Each month, we estimated emissions from within the NAME domain (at the resolution explored by the trans-dimensional method), as well as offsets to a priori boundary condition “curtains” on each edge of the domain (Sect. 2.4). In addition, an offset parameter was included to account for any differences between the satellite and the calibrated ground-based measurements and their representation by models. A normal PDF was defined for both of these types of offsets, centred around zero, and where the standard deviations of the PDFs were governed by hyper-parameters. The standard deviations of the boundary condition offsets were allowed to vary up to 100.0 nmol mol−1 and up to 50.0 nmol mol−1 for the offset between surface and satellite data, respectively. The necessity of this parameter to produce the most robust results is discussed in Sect. 3.3.

The Metropolis–Hastings MCMC sampler was run with 500 000 iterations with the initial 100 000 samples discarded as burn-in. Every 500th iteration was saved and used to build posterior PDFs for each parameter. The mean and 2.5–97.5 percentiles were used to produce posterior estimates and 95 % confidence intervals.

2.4 A priori fields

A priori emissions and boundary condition fields are summarised in Table 1. Emissions were inferred for the three major source sectors in Brazil: anthropogenic, biomass burning and wetlands. Maps for two representative months in the wet (January) and dry (September) seasons for 2014 are shown for each sector in Fig. 2.

Anthropogenic emissions, excluding biomass burning, were from the EDGAR (Emission Database for Global Atmospheric Research) v4.3.2 database (Janssens-Maenhout et al.2019). Annual emissions were available up to 2012 and then assumed to be equal to the 2012 emissions thereafter. The biomass burning contribution was from GFED (Global Fire Emissions Database) v4.1 (van Der Werf et al.2017) at monthly resolution to the year 2015 and assumed to be held at 2015 values thereafter.

Wetland emissions were based on the output from the JULES land surface model (Clark et al.2011), which was modified to use the wetland fractional map from Surface WAter Microwave Product Series (SWAMPS). We used a version of SWAMPS that was updated from Schroeder et al. (2015) to include wetlands occurring under dense canopies, to remove rice agriculture and to include any inland water. Wetland emissions across South America were scaled to 44 Tg yr−1 based on the mean bottom-up estimate for Tropical South America from Saunois et al. (2016).

Figure 2A priori emissions for wetlands, anthropogenic and biomass burning sectors for January and September 2014. These months are representative of peak wetland extent in the wet season (January) and biomass burning activity (September). Note that the EDGAR anthropogenic inventory is annual resolution and is shown here for 2012.

A priori mole fractions at the boundaries of the domain were derived from the CAMS CH4 flux inversion product v17r1 (accessible at, last access: 22 October 2020). This version assimilated the global surface measurement network and did not use satellite data. This product was only available up to 2017, so to extend the analysis to 2018, the climatological mean of the 2010–2017 period was used.

Table 1Summary of a priori emissions and boundary conditions used in this study. All maps have been regridded to 0.23 latitude by 0.35 longitude resolution.

a Repeats this year thereafter.
b For 2018, a climatological mean of the 2011–2017 period was used.

Download Print Version | Download XLSX

2.5 Sector attribution

The total emissions estimated from the inversion were partitioned into each of the three major source sectors using the fraction of each source in the a priori emission fields in each grid cell. Due to the largely distinct spatial or temporal distributions of the sectors as shown in Fig. 3, the fractional map of each source is not overly dependent on the inventories used. The influence of the a priori distributions on the robustness of the sector partitioning is discussed in Sect. 3.4.

Figure 3Fraction of wetland, anthropogenic and biomass burning emissions for January (a–c) and September (d–f) 2014 derived from the a priori emissions. Note that the EDGAR anthropogenic inventory is annual resolution and is shown here for 2012.

2.6 Validation with ATTO

To provide a validation of the inversion results, we compared a model prediction of mole fractions at ATTO derived from hourly NAME sensitivities convolved with our posterior emissions maps and boundary conditions against measured values. Four tests were run using different configurations of the inversion. The first three estimates were from inversions that used variants of the GOSAT and RPB dataset. The first inversion utilised GOSAT data alone. The second inversion used both GOSAT and RPB measurements but did not include an offset parameter between satellite and surface data in the inversion. The third inversion used GOSAT and RPB measurements and included an offset parameter that was estimated in the inversion (our main results). These tests and the resulting comparisons with ATTO data allowed us to determine the factors that are most important when using satellite data to constrain country-scale emissions. We performed a final test which scaled our posterior emissions map so that emissions from the Brazilian Amazon matched those derived by Wilson et al. (2016) using four aircraft sites. For the whole Amazon basin, the lowest value in the range presented in Wilson et al. (2016) of 31.6 Tg yr−1 was used, with 19Tgyr-1 coming from the Brazilian Amazon, based on wetland extent. This test allowed us to investigate the fit of previous results against ATTO data.

In addition to these experiments, we simulated the model prediction at ATTO using a second regional Lagrangian model, the FLEXible PARTicle dispersion model (FLEXPART), for 2014–2017 (Pisso et al.2019). The setup for FLEXPART was the same as NAME, except the surface was defined as 0–50 m above ground level and the meteorological drivers were 1 resolution from the European Centre for Medium-Range Weather Forecasts (ECMWF). Particles were tracked backwards for 30 d. This test allowed us to assess whether results are significantly impacted by systematic uncertainties in NAME.

2.7 Sensitivity Studies

Sensitivity tests against a range of inputs to the inversion were performed to assess the robustness of our results. Three categories of inputs were tested: a priori emissions, a priori boundary conditions and the model XCO2 fields used to derive XCH4. In most cases, comparisons were performed for 2014 only, but if differences were seen, the analysis was expanded across the entire time range of 2010–2018. The sensitivity study details are summarised in Table 2.

To test the sensitivity to a priori emissions, we ran a set of inversions where emissions were perturbed one at a time from each source sector. We changed the magnitudes of emissions from each sector and tested variations of wetland extent maps. For the latter, three additional wetland distributions were used: two using JULES emissions either with Bergamaschi et al. (2007) (hereafter referred to as Kaplan, which is based on land cover maps from optical imagery) or the high-resolution Tropical and Sub-Tropical Wetland Distribution v2.0 (Gumbricht et al.2017, hereafter referred to as Gumbricht). The final variation used the Wetland Methane Emissions and Uncertainty dataset for atmospheric chemistry an transport modelling (WetCHARTs) dataset based on an ensemble of wetland models (Bloom et al.2017). We did not modify these wetland distributions to include any emissions that might occur when the water table is below the surface. These four wetland distribution maps are shown in Fig. A2 in the appendix.

To test the sensitivity to a priori boundary conditions, we used a variation of the global mole fractions used to generate the boundary condition curtains. We used the climatological mean of the MOZART global model (Emmons et al.2010) over the 2010–2014 time period. The setup for MOZART is described in Palmer et al. (2018).

To test the sensitivity to the model XCO2 used to derive XCH4, we generated 10 variations of XCH4 for each measurement. These were created by randomly selecting between the median (the main results) and the extremes in the ensemble members that are included with the data product. We reran the inversion for each of the 10 datasets for the full 2010–2018 time period, which allowed us to investigate random errors in XCO2. However, additional uncertainties could nevertheless remain due to sparse CO2 observations in the region.

Bergamaschi et al. (2007)Gumbricht et al. (2017)(Bloom et al.2017)(Emmons et al.2010; Palmer et al.2018)

Table 2Sensitivity studies performed in this study for a priori fields or for model XCO2. For the a priori emissions, each sector was varied one at a time over the full South America domain with other sectors kept in their original configurations.

 Based on bottom-up emissions estimates from Saunois et al. (2016) for Tropical South America.

Download Print Version | Download XLSX

3 Results

3.1 Annual and seasonal emissions by sector

Mean emissions from 2010 to 2018 for Brazil are 33.6±3.6Tgyr-1 (Fig. 4). These emissions correspond to mean anthropogenic emissions of 19.0±2.6Tgyr-1, mean wetland emissions of 13.0±1.9Tgyr-1, and mean biomass burning emissions of 1.7±0.3Tgyr-1. Maps of these posterior emissions and the difference from the a priori inputs are shown for each season in Fig. A3 in the appendix.

Both our 2012 and 2015 estimates of anthropogenic emissions of 16.2±3.0 and 18.3±2.5Tgyr-1 are consistent within uncertainties with Brazil's Third Biennial Update Report to the UNFCCC, which estimates 15.6 Tg yr−1 in 2012 and 16.3±3.8Tgyr-1 in 2015 (Ministry of Foreign Affairs et al.2019), when LULUCF (land use, land-use change and forestry) contributions are removed.

The overall rise in emissions over the 2010–2018 period generally occurred in late 2013 and early 2014 and was sustained thereafter. Average emissions during 2014–2018 rose over 20112013 levels by 6.9±5.3Tgyr-1, and this is driven by changes in anthropogenic, wetland and biomass burning emissions of 3.3±3.7, 2.6±2.8 and 1.0±0.4Tgyr-1, respectively.

Across 2010–2018, we find that total emissions maximise in April and minimise in October, and the overall seasonality reflects the net effect of different seasonal patterns in the three sectors. Anthropogenic emissions (the largest sector) peak in April, are lowest in August–October (dry season) and could be a result of seasonality in cattle, manure management (e.g. Cardoso et al.2019) or landfill emissions (e.g. Machado et al.2009; Imbiriba et al.2020). Anthropogenic emissions are only estimated annually in EDGAR and in reports to the UNFCCC and thus do not capture this important feature. Wetland emissions peak during the wet season between February and April and are lowest in October, and this seasonality is more pronounced in our estimates than in the a priori emissions. Anthropogenic and wetland emissions are discussed further for different regions of Brazil in Sect. 3.2. Biomass burning emissions maximise in September and the seasonality is consistent with GFED.

Our analysis shows that individual years show some differences from the bottom-up estimates. We find the largest biomass burning emissions in 2010, a year with strong drought and intensive burning due to high Atlantic sea surface temperatures (Lewis et al.2011; van der Laan-Luijkx et al.2015); annual mean emissions in 2010 were 5.5±0.5Tgyr-1 (based on April–December due to a change in model resolution of the UM and the inability to model GOSAT measurements from January–March) but with a monthly value in September at 22.3-2.7+2.6Tgyr-1, a value that is 6.3 Tg yr−1 larger than reflected in GFED. Our estimates are consistent with GFED at most other times. Wetland emissions are highest in 2015, which corresponds to a strong El Niño year. The a priori model emissions do not capture the increase in 2015 but do simulate a decrease from 2016. This feature is discussed further in Sect. 3.2.

Figure 4Brazil's CH4 emissions derived from GOSAT and RPB measurements between 2010 and 2018. Total emissions (grey) are split into the three major sectors: wetlands (cyan), anthropogenic (yellow) and biomass burning (red). Prior and posterior emissions are dashed and solid lines, respectively. Shading indicates the 95 % confidence interval. (a) Monthly emissions, (b) monthly emissions smoothed with a 12-month rolling mean and (c) monthly means across the 2010 to 2018 period. Errors in mean values assume a 50 % correlation between individual months.


The performance of the inversion is demonstrated through a comparison of modelled mole fractions derived from the posterior emissions and boundary conditions with the measurements used in inversion. We show this fit for both GOSAT and RPB in Fig. A4 in the appendix, and we find both datasets to be represented well by the inversion.

3.2 Sub-national emissions

In addition to the Brazilian totals presented above, we aggregated our posterior emissions for the major regions of Brazil: the Amazon basin, the Pantanal and the remainder of the country (Figs. 5 and 6 for wetland and anthropogenic sectors, respectively). The Amazon basin was defined using the TransCom definition for Tropical South America (Saunois et al.2016), and the Pantanal region was defined using the TRIP River Routing Model output (Oki et al.1999). These regions were further masked to only include the area within Brazil using the public domain Natural Earth database (, last access: 22 October 2020).

Figure 5Brazil's wetland emissions aggregated over Amazon and Pantanal regions. (a) Monthly emissions with an inset map showing the masks used to delineate between the Amazon and Pantanal regions, (b) monthly emissions smoothed with a 12-month rolling mean and (c) seasonal means across the 2010 to 2018 period. Errors for mean values assume a 50 % correlation between individual months.


We aggregated wetland emissions (Fig. 5) into mean values for the 2010–2018 period, changes between 20112013 and 2014–2018, and means for each month. Mean wetland emissions from the Brazilian Amazon and Pantanal regions across the 2010–2018 period are 9.2±1.8 and 1.9±0.5Tgyr-1, respectively. Wetland emissions in the Amazon and Pantanal comprise 65 % and 26 % of total emissions, respectively, with emissions from the Pantanal being dominated by the anthropogenic sector. While emissions from the Pantanal are not significantly different from the a priori emissions, these results are found to be robust in our sensitivity studies as discussed in Sect. 3.4. There is only a small change in wetland emissions over the two regions between 20112013 and 2014–2018. Differences are 1.5±2.6 and 1.0±0.7Tgyr-1, for the Amazon and Pantanal, respectively. Only the small Pantanal change is significant within the 95 % confidence interval. We also find that there is an offset in peak emissions between the Amazon and the Pantanal regions. Amazon wetland emissions peak around February–March, whereas the Pantanal peaks in April. The seasonality for the Amazon is earlier than reflected in the a priori emissions.

Because wetland emissions from the Pantanal exhibit a similar seasonal pattern to the seasonality in anthropogenic emissions across Brazil, we analysed the regions that are responsible for driving the anthropogenic seasonal cycle. Figure 6 shows the anthropogenic emissions aggregated over the Amazon and Pantanal regions and over the remaining Brazilian territory. We find that the seasonal cycle is dominated by the emissions outside of the Amazon and the Pantanal. Therefore, while Pantanal wetland and anthropogenic emissions have the same seasonal pattern, this is not due to a misattribution between sectors based on the current configuration of the wetland and anthropogenic prior emissions. However, it is important to note the difficulty of wetland models in capturing the full seasonal cycle in the Pantanal due to overbank inundation, so there could be some uncertainty in the fractional partitioning due to uncertainty in the wetland models used in the main results and in the sensitivity studies (Parker et al.2018).

Figure 6Brazil's anthropogenic emissions aggregated over the Amazon basin, Pantanal region and the rest of Brazilian territory. (a) Monthly emissions with an inset map showing the masks used to delineate different regions; (b) monthly emissions smoothed with a 12-month rolling mean and (c) seasonal means across the 2010 to 2018 period. Errors for mean values assume a 50 % correlation between individual months.


Wetland emissions are 3.9±3.0Tgyr-1 larger in the 2015 wet season relative to 20112014, a feature that is not present in the a priori emissions. We show in Fig. 7 that this increase is driven from the Amazon and not by Pantanal wetlands. We investigated changes in some of the major environmental influences to understand what could drive this pattern. Figure 7 shows our derived emission maps, changes in surface temperature from the WFDEI (WATCH Forcing Data methodology applied to ERA-Interim) meteorological dataset (, last access: 22 October 2020) and changes in SWAMPS inundation for the wet season, defined as February–April (FMA). We show differences between 2015 and 20112014 and between 2016 and 2015 for the Amazon and the Pantanal regions.

We find that the increase in 2015 originates mainly from the Western Amazon (defined as the Brazilian Amazon area west of 55 W) with a rise of 3.7±2.7Tgyr-1. This coincides with increased surface temperatures from this region. Wetland extent did not significantly change in the Western Amazon between 2015 and preceding years. Emissions then decrease after mid-2015 to levels that are sustained from 2016 to 2018. We find that this decrease is correlated with both lower soil temperatures and decreased inundation. The a priori emissions may be simulating the decrease after 2016 because the a priori emissions are constrained to the observational inundation fields. However, these results suggest that there may be uncertainties in the wetland model temperature sensitivity.

Figure 7Differences during the wet season (FMA). (a–c) 2015 minus 20112014 average and (d–f) 2016 minus 2015. (a, d) CH4 emissions, (b, e) surface temperature from the WFDEI meteorological dataset and (c, f) SWAMPS wetland fraction. All regions outside of the Brazilian Amazon and Pantanal have been masked for clarity.

3.3 Validation against ground-based data

We used independent data from ATTO to assess the robustness of our inversion results and to understand what factors are important for the inversion setup. The results of these tests can be seen in Fig. 8. Other datasets besides ATTO exist, such as aircraft data from the Amazon (Wilson et al.2016; Pangala et al.2017) but were not available for use.

An inversion using only GOSAT data produced a mean difference between modelled ATTO data and measurements of 42.5 nmol mol−1 (Fig. 8a and e). This difference can largely be attributed to modelled boundary conditions that are consistently elevated throughout the year above the lowest ATTO data. Introduction of the surface baseline station of RPB (Fig. 8b and f) improved the boundary condition estimation, with the modelled boundary conditions now consistent with ATTO data in most months and lower than ATTO data in months with significant regional emissions (i.e. times when ATTO may not be representative of boundary conditions). Despite consistency with boundary conditions, this setup produced the highest mean difference with ATTO, 67.4 nmol mol−1, due to large regional emissions being estimated. The third case, the setup of our main results, which allowed for an offset between the GOSAT and RPB measurements to be estimated in the inversion, resulted in the best fit to ATTO (Fig. 8c and g). The model achieved consistent boundary conditions and the smallest mean difference with ATTO (18.9 nmol mol−1). In our inversions from 2010 to 2018, we estimate a mean offset parameter between GOSAT and RPB data of 22±8nmolmol-1. The numbers presented for the offsets are a combination of any bias between the data themselves but also in the model's interpretation of these datasets. The model simulates the three-dimensional atmospheric fields necessary to combine these two datasets together. However, the interpretation of these tests show that near-surface data that help to constrain boundary conditions are required because GOSAT data alone do not have enough resolving power to partition boundary conditions and emissions. An offset parameter should then be included to account for a combination of any differences between in situ data and satellite data and any offsets due to the atmospheric model.

When our posterior emissions estimates were scaled to match previous results derived by Wilson et al. (2016) (but keeping the posterior boundary conditions fixed from our main results), a larger offset from ATTO of 45.4 nmol mol−1 (Fig. 8d and h) again resulted. This test indicates that larger emissions from the Amazon are inconsistent with ATTO and its representation by the NAME model.

To assess the possibility of large systematic uncertainties in NAME, we show a comparison of the validation at ATTO generated using NAME (as in Fig. 8c) with those generated using FLEXPART. This comparison is shown in Fig. A5 in the appendix and shows that the posterior emissions and boundary conditions derived here are consistent with ATTO across both models. These results provide additional confidence in the magnitude of emissions that we derive.

Figure 8(a–d) Comparison of modelled and measured mole fractions at ATTO for different posterior emissions with (e–h) histograms of the differences. Posterior emissions derived in inversions using (a, e) GOSAT measurements only, (b, f) both GOSAT and RPB data but no offset parameter in the inversion, and (c, g) both GOSAT and RPB data with an offset parameter between the two datasets in the inversion. (d, h) Posterior emission distribution of our main results but scaled so that total CH4 emissions in the Brazilian Amazon were equivalent to those derived in Wilson et al. (2016).


3.4 Sensitivity studies

3.4.1 Sensitivity to a priori emissions

Sensitivity tests of the effect of different wetland distributions are shown in Fig. 9. Total emissions do not change significantly between these sensitivity tests, despite the large seasonal cycle in the a priori WetCHARTs emissions that is not reflected in the other wetland distributions. This suggests that the inversion is well constrained by the atmospheric data and is not significantly influenced by the prior emissions. There are small differences in wetland and anthropogenic partitioning, but emissions are consistent within uncertainties.

In addition to different a priori wetland distributions, sensitivity tests to perturb a priori emissions from each source sector are shown in Figs. A6, A7 and A8 in the appendix. The main impact of perturbing a priori emissions from any source sector comes in the partitioning of total emissions into the sources, particularly between anthropogenic and wetland emissions. This is due to a small overlap between anthropogenic and wetland sources (Figs. 2 and 3). However, the trade-off between these two sectors is smaller than the initial perturbation to the prior and emissions are still consistent within confidence intervals, suggesting that the sectoral partitioning is robust. The largest sensitivity to the a priori emissions is shown when doubling a priori biomass burning emissions (Fig. A8), and the resulting posterior biomass burning estimate is not consistent within uncertainties to the unperturbed case. Overall, these tests show that our results and the associated sectoral partitioning, with the exception of some influence from the biomass burning prior emissions, are robust to the a priori emissions used.

Figure 9Brazil's CH4 emissions in a sensitivity inversion using perturbed a priori wetland distributions. All other components were held at the original configuration. Total emissions for all of South America were scaled to 44 Tg yr−1 for all distributions apart from WetCHARTs, which used its derived emissions. (a) Total, (b) wetland, (c) anthropogenic and (d) biomass burning emissions. Prior emissions are only shown in (a, b) as they only vary for these components in this test.


3.4.2 Influence of a priori boundary conditions

Results of using a different global model of a priori boundary conditions are shown in Fig. A9 in the appendix. Due to differences in the seasonal cycle when comparing the CAMS and MOZART boundary conditions in 2014, this analysis was run for the full 2010–2018 period to provide a longer comparison. While there is some month-to-month variability, the overall patterns are consistent between the two inversion setups, suggesting that the inversion is robust to the a priori boundary conditions. As demonstrated in Sect. 3.3, it is important to include data that can help the inversion constrain the boundary conditions, through, for example, surface measurements from remote background stations.

3.4.3 Influence of model XCO2 on XCH4

We generated 10 variations of the XCH4 dataset used in the inversion based on different model XCO2. Figure 10 shows the emissions estimates that result when perturbing XCH4 by random values of the model XCO2 used to generate XCH4 using the CO2 proxy method. Because of some differences in 2014, this analysis was run for the full 2010–2018 period.

Across the 10 variations, mean emissions over 2010–2018 range from 33.834.8 Tg yr−1 in total, corresponding to 19.019.4 Tg yr−1 for anthropogenic, 13.013.4 Tg yr−1 for wetlands and 1.71.8 Tg yr−1 for biomass burning. Individual months can exhibit larger ranges in the 10 variants, in some cases spanning >10Tgyr-1. The differences based on model XCO2 does not exhibit any particular seasonality. The change between the 20112013 and 2014–2018 periods across these 10 inversions produces a range of 5.97.0 Tg yr−1. Thus, the increase in emissions is robust to uncertainties in XCO2.

Figure 10Brazil's CH4 emissions in a sensitivity inversion using 10 variants of GOSAT data generated with different model XCO2. All other components were held at the original configuration. Each line represents a different inversion run. (a) Total, (b) wetland, (c) anthropogenic and (d) biomass burning emissions. Prior emissions do not vary in this test.


4 Discussion

We find that Brazil's emissions increased during 2014–2018 over 20112013 levels by 6.9±5.3Tgyr-1, and this coincides with a large increase in global CH4 mole fraction growth rate in 2014 (Nisbet et al.2019). The increase in Brazil's emissions is primarily driven by anthropogenic and wetland sources. Brazil's anthropogenic emissions are dominated by agriculture and mainly cattle, which is likely to be the main source for the inferred anthropogenic change. However, we did not have sufficient information with which to robustly separate total anthropogenic emissions into individual sub-sectors. Future work should couple measurements of δ13C-CH4 from Brazil along with campaigns to sample representative isotopic source signatures (Ganesan et al.2018) to better understand whether changes in these sources are consistent with isotopic constraints.

The increased wetland emissions that we derive in the wet season of 2015 primarily originates from the Western Amazon. Previous studies have found that changes in wetland CH4 emissions exhibit complex dynamics during El Niño years. Zhang et al. (2018) found through model simulations that the 20152016 El Niño led to larger instantaneous growth in CH4 emissions than previous El Niño periods. This study also showed that there was a large increase in the Western Amazon due to increased soil respiration from high soil temperatures, despite a decline in wetland extent due to drought. This pattern is consistent with the results that we have derived from atmospheric data rather than from model simulations. In contrast to Zhang et al. (2018), who also found a 20152016 El Niño effect on Western Amazon emissions, we find that emissions increased during the 2015 wet season rather than the 2016 wet season. We instead show a decline in 2016 emissions, surface temperature and wetland extent compared to 2015 levels. This discrepancy in temporal response from Zhang et al. (2018) suggests that the dynamics of the wetland response to climatic perturbations may require further investigation.

Our results show that emissions can be derived for a country of the size of Brazil from satellite data coupled with high-resolution atmospheric transport modelling, but careful consideration needs to be paid to the setup of the inversion. Regional inversions use atmospheric data to estimate boundary conditions and regional emissions. Due to the lower signal-to-noise ratio of GOSAT data (which are sensitive to surface emissions that are mixed through the entire atmospheric column) compared to ground-based data (although the reduced surface sensitivity and precision of satellite data needs to be weighed against the greater geographical coverage), we find that additional surface data are required to better constrain the boundary conditions. However, we find that when combining satellite data with calibrated surface data in an inversion, it is critical to incorporate an offset parameter between the two datasets in the inversion. The GOSAT product used here has been previously corrected by 7.7 nmol mol−1 as a global average offset to independent ground-based measurements from the Total Carbon Column Observing Network (TCCON) (Wunch et al.2011). However, large regional variations can still exist (Dils et al.2014). The estimated offset in an inversion is due to biases between the different datasets as well as their representation by the atmospheric transport model. Other inversion studies have imposed latitude-dependent bias corrections on other GOSAT data (Bergamaschi et al.2009; Turner et al.2015).

Janardanan et al. (2020) estimated Brazil's CH4 emissions using a coupled global Eulerian–Lagrangian model from 20112017 using GOSAT and surface data and find total emissions to be 56.2 Tg yr−1 compared with 33.3±3.7Tgyr-1 derived in this study. The difference between our results can be attributed to the natural wetland emissions estimates for which Janardanan et al. (2020) derive 39.8±12.4Tgyr-1 compared to 13.1±1.9Tgyr-1 presented here. Anthropogenic estimates (excluding biomass burning) are similar at 16.5 Tg yr−1 compared with our estimate of 18.8±2.6Tgyr-1. One factor in this difference could be the differing GOSAT retrieval products used which were derived using different algorithms (CO2 proxy vs. full physics retrievals). Another reason for the discrepancy could stem from Janardanan et al. (2020) not allowing for an offset parameter between the surface/aircraft and satellite data within their inversion. In the case where we similarly set up our inversion to not include an offset parameter (as shown by the ATTO comparison in Fig. 8b), we also derive larger total emissions of 51.7±3.5Tgyr-1 for 2014. Not allowing for this offset produces the poorest comparison to the independent ATTO measurements. However, it is important to note that our validation is based on only one site because of the availability of data.

Studies deriving Amazon basin CH4 emissions using aircraft data from within the Amazon are also higher than our estimates at 49 Tg yr−1 (Miller et al.2007; Wilson et al.2016). However, these higher estimates, as shown in Fig. 8d, when simulated with NAME, are less consistent when compared with CH4 mole fractions measured at the ATTO tower. The wetland results presented here are most consistent with the lower bound estimates from Saunois et al. (2016) which range from 23.463.7 within Tropical South America. As discussed in Sect. 3.4.1, neither varying the magnitude of the prior input for wetlands nor the wetland extent map used significantly altered our posterior estimates.

We propose one reason for the difference from aircraft-based estimates could be that the studies using aircraft data may not be able to constrain emissions over the whole of the Amazon basin and furthermore at the country scale, though our comparison at present has only been validated by one in situ measurement station. Future work should perform a detailed comparison between aircraft-derived estimates and those derived from satellites, investigating the inversion setup and the degree of constraint by the datasets to understand the reasons for this discrepancy. The main benefit of using satellite data is in their widespread coverage, which allows for country-scale emissions to be derived (albeit with inclusion of calibrated near-surface data in the inversion).

Overall, we derive lower emissions than previous studies. We show the validation of our results at ATTO using two models, NAME and FLEXPART. The consistency between the two models in simulating the magnitude of mole fractions at ATTO provides some confidence in the lower emissions we derive over previous studies. In future, performing a full set of inversion results using a large range of models with different physical parameterisations could help to quantify the magnitude of any systematic uncertainties.

5 Conclusions

We estimated Brazil's CH4 emissions from 2010 to 2018 using a combination of GOSAT satellite data and surface data from Ragged Point, Barbados. Due to the spatial and temporal separation in the three main sources of Brazil's emissions (anthropogenic, wetland and biomass burning), we were able to derive emissions estimates by sector.

We find mean emissions from 2010 to 2018 to be 33.6±3.6Tgyr-1, corresponding to 19.0±2.6Tgyr-1 from anthropogenic, 13.0±1.9Tgyr-1 from wetland and 1.7±0.3Tgyr-1 from biomass burning. We find a rise of 6.9±5.3Tgyr-1 occurring between the 20112013 and 2014–2018 periods. Both anthropogenic and wetland sources drive the increase in emissions over the period. This rise in emissions occurred during a period of accelerated global CH4 growth, suggesting that Brazil's CH4 sources have a significant influence on changes in the atmosphere.

We find that wetland emissions from the Western Amazon increased by 3.7±2.7Tgyr-1 in the 2015 wet season, at the start of the 20152016 El Niño, and decreased subsequently from 2016. We show that the increase is likely to be driven by increased surface temperatures (and thus respiration rates) rather than through changes in inundation.

Our study demonstrates that satellite data, with its enhanced coverage compared to surface data, can be used to infer country-scale emissions. This is beneficial for independently comparing top-down estimates with national reports to the UNFCCC. However, we show that satellite data must be used in conjunction with calibrated surface data, which provide critical constraints on boundary conditions in regional inversions. It is also necessary to account for any offsets between datasets which can result from either biases between satellite data and surface data or from the atmospheric transport model used to simulate these data. Otherwise the resulting emissions estimates may be biased. Our sensitivity studies show that our emissions estimates are insensitive to most inputs, but the largest differences are driven by uncertainties in the model XCO2 used to derive XCH4.

Appendix A

Figure A1Annual mean NAME sensitivity map for GOSAT measurements in nadir mode within the area from 35.8 S to 7.3 N and from 76.0 to 32.8 W over Brazil for 2014.

Figure A2CH4 emissions for each of the four a priori wetland emissions used in the wetland extent sensitivity study. Panels (a)(c) SWAMPS, Kaplan and Gumbricht fractional maps are combined with the JULES emissions output. Details of these inversion setups are described in Table 2. This is shown for April 2014 which is a wet season month with high emissions in the Amazon basin and the Pantanal region.

Figure A3Emissions maps across the 20112018 time period grouped into February–April (FMA), May–July (MJJ), August–October (ASO) and November–January (NDJ). Panels (a)(d) CH4 posterior emissions maps and (e)(h) difference between the posterior and the a priori emissions input.

Figure A4Modelled mole fractions derived from the posterior emissions estimate compared to measurements from (a) GOSAT and (b) RPB. The orange line shows posterior boundary conditions and the blue line shows the total modelled mole fraction. Measurements are displayed as red dots.


Figure A5Comparison of modelled and measured mole fractions at ATTO with the same posterior emissions and boundary conditions convolved with sensitivity maps derived from two different models. (a) NAME model (main results) and (b) FLEXPART model where (c)(d) shows histograms of the difference. Posterior emissions were derived from our inversion setup using both GOSAT and RPB data with an offset parameter between the two datasets allowed within the inversion.


Figure A6Brazil's CH4 emissions in a sensitivity inversion using a perturbed a priori wetlands emissions magnitude. All other components were held at the original configuration. A priori wetland emissions for all of South America were scaled to 32 (low), 44 (mean) or 50 (high) Tg yr−1 (as defined for Tropical South America in Saunois et al.2016). (a) Total, (b) wetland, (c) anthropogenic and (d) biomass burning emissions. Prior emissions are only shown in (a, b) as they only vary for these components in this test.


Figure A7Brazil's CH4 emissions in a sensitivity inversion using perturbed a priori anthropogenic emissions. All other components were held at the original configuration. Anthropogenic emissions for all of South America were doubled from EDGAR. (a) Total, (b) wetland, (c) anthropogenic and (d) biomass burning emissions. Prior emissions are only shown in (a, c) as they only vary for these components in this test.


Figure A8Brazil's CH4 emissions in a sensitivity inversion using perturbed a priori biomass burning emissions. All other components were held at the original configuration. Biomass burning emissions for all of South America were doubled from GFED. (a) Total, (b) wetland, (c) anthropogenic and (d) biomass burning emissions. Prior emissions are only shown in (a, d) as they only vary for these components in this test.


Figure A9Brazil's CH4 emissions in a sensitivity inversion using perturbed a priori boundary conditions from the MOZART model. All other components were held at the original configuration. (a) Total, (b) wetland, (c) anthropogenic and (d) biomass burning emissions. A priori emissions do not vary in this test.


Code and data availability

University of Leicester GOSAT Proxy XCH4 data can be accessed via the Copernicus Climate Data Store or by contacting Rob Parker. RPB data can be accessed from (Prinn et al.2018a) and by contacting Dickon Young ( ATTO data can be accessed from (last access: September 2018) and by contacting Jošt Lavrič ( The inversion code and NAME footprints used in this study can be accessed by contacting Rachel Tunnicliffe ( and Anita Ganesan (

Author contributions

RLT and ALG designed the methodology and wrote the article. RLT performed the analysis. RJP and HB provided the GOSAT data. NG provided JULES CH4 emission fields. BP and ZZ provided SWAMPS wetland extent maps. JVL and DW provided data from the Amazon Tall Tower Observatory. MR advised on methodology. DY and SO'D provided data from Ragged Point, Barbados.

Competing interests

The authors declare that they have no conflict of interest.


This work and its contributors (Rachel L. Tunnicliffe) were supported by the Newton Fund through the Met Office Climate Science for Service Partnership Brazil (CSSP Brazil) and the Natural Environment Research Council (NERC) Methane Observations and Yearly Assessments programme (MOYA, NE/N016548/1). Anita L. Ganesan was funded by the NERC Independent Research Fellowship NE/L010992/1. Robert J. Parker and Hartmut Boesch were funded via the UK National Centre for Earth Observation (NCEO grant no. nceo020005).

We thank the Japanese Aerospace Exploration Agency, National Institute for Environmental Studies and the Ministry of Environment for the GOSAT L1B data and their continuous support as part of the Joint Research Agreement. We thank Alistair Manning (Met Office) for useful and insightful discussions about this work.

The operation of the Ragged Point site was funded by the National Aeronautical and Space Administration (NASA, USA) (grants NAG5-12669, NNX07AE89G and NNX11AF17G to MIT; grants NAG5-4023, NNX07AE87G, NNX07AF09G, NNX11AF15G and NNX11AF16G to SIO) under the AGAGE programme and the National Oceanic and Atmospheric Administration (NOAA, USA) (contract RA-133R-15-CN-0008 to the University of Bristol). ATTO data were supported by the Max Planck Society (MPG), the German Federal Ministry of Education and Research (contracts 01LB1001A and 01LK1602A) and the Brazilian Ministério da Ciência, Tecnologia e Inovação (MCTI/FINEP contract 01.11.01248.00) as well as the Amazon State University (UEA) (FAPEAM, LBA/INPA and SDS/CEUC/RDS-Uatumã). We acknowledge the Swiss Federal Office for Meteorology and Climatology (MeteoSwiss) for providing access to ECMWF ERA-Interim reanalysis products for use with the FLEXPART model. This work was carried out using the computational facilities of the Advanced Computing Research Centre, University of Bristol – (last access: 22 October 2010). GOSAT retrievals used the ALICE high-performance computing facility at the University of Leicester.

Financial support

This research has been supported by the Newton Fund (grant no. CSSP Brazil) and the Natural Environment Research Council (NERC) Methane Observations and Yearly Assessments programme (MOYA, NE/N016548/1).

Review statement

This paper was edited by Ilse Aben and reviewed by two anonymous referees.


Andreae, M. O., Acevedo, O. C., Araùjo, A., Artaxo, P., Barbosa, C. G. G., Barbosa, H. M. J., Brito, J., Carbone, S., Chi, X., Cintra, B. B. L., da Silva, N. F., Dias, N. L., Dias-Júnior, C. Q., Ditas, F., Ditz, R., Godoi, A. F. L., Godoi, R. H. M., Heimann, M., Hoffmann, T., Kesselmeier, J., Könemann, T., Krüger, M. L., Lavric, J. V., Manzi, A. O., Lopes, A. P., Martins, D. L., Mikhailov, E. F., Moran-Zuloaga, D., Nelson, B. W., Nölscher, A. C., Santos Nogueira, D., Piedade, M. T. F., Pöhlker, C., Pöschl, U., Quesada, C. A., Rizzo, L. V., Ro, C.-U., Ruckteschler, N., Sá, L. D. A., de Oliveira Sá, M., Sales, C. B., dos Santos, R. M. N., Saturno, J., Schöngart, J., Sörgel, M., de Souza, C. M., de Souza, R. A. F., Su, H., Targhetta, N., Tóta, J., Trebs, I., Trumbore, S., van Eijck, A., Walter, D., Wang, Z., Weber, B., Williams, J., Winderlich, J., Wittmann, F., Wolff, S., and Yáñez-Serrano, A. M.: The Amazon Tall Tower Observatory (ATTO): overview of pilot measurements on ecosystem ecology, meteorology, trace gases, and aerosols, Atmos. Chem. Phys., 15, 10723–10776,, 2015. a, b

Beck, V., Gerbig, C., Koch, T., Bela, M. M., Longo, K. M., Freitas, S. R., Kaplan, J. O., Prigent, C., Bergamaschi, P., and Heimann, M.: WRF-Chem simulations in the Amazon region during wet and dry season transitions: evaluation of methane models and wetland inundation maps, Atmos. Chem. Phys., 13, 7961–7982,, 2013. a

Bergamaschi, P., Frankenberg, C., Meirink, J. F., Krol, M., Dentener, F., Wagner, T., Platt, U., Kaplan, J. O., Körner, S., Heimann, M., Dlugokencky, E. J., and Goede, A.: Satellite chartography of atmospheric methane from SCIAMACHY on board ENVISAT: 2. Evaluation based on inverse model simulations, J. Geophys. Re., 112, D02304,, 2007. a, b

Bergamaschi, P., Frankenberg, C., Meirink, J. F., Krol, M., Villani, M. G., Houweling, S., Dentener, F., Dlugokencky, E. J., Miller, J. B., Gatti, L. V., Engel, A., and Levin, I.: Inverse modeling of global and regional CH4 emissions using SCIAMACHY satellite retrievals, J. Geophys. Res., 114, D22301,, 2009. a, b

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

Botía, S., Gerbig, C., Marshall, J., Lavric, J. V., Walter, D., Pöhlker, C., Holanda, B., Fisch, G., de Araújo, A. C., Sá, M. O., Teixeira, P. R., Resende, A. F., Dias-Junior, C. Q., van Asperen, H., Oliveira, P. S., Stefanello, M., and Acevedo, O. C.: Understanding nighttime methane signals at the Amazon Tall Tower Observatory (ATTO), Atmos. Chem. Phys., 20, 6583–6606,, 2020. a

Calvo Buendia, E., Tanabe, K., Kranjc, A., Baasansuren, J., Fukuda, M., S., N., Osako, A., Pyrozhenko, Y., Shermanau, P., and Federici, S.: 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories, available at: (last access: 22 October 2020), 2019. a

Cardoso, A. d. S., Oliveira, S. C., Janusckiewicz, E. R., Brito, L. F., Morgado, E. D. S., Reis, R. A., and Ruggieri, A. C.: Seasonal effects on ammonia, nitrous oxide, and methane emissions for beef cattle excreta and urea fertilizer applied to a tropical pasture, Soil Till. Res., 194, 104341,, 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

Clark, D. B., Mercado, L. M., Sitch, S., Jones, C. D., Gedney, N., Best, M. J., Pryor, M., Rooney, G. G., Essery, R. L. H., Blyth, E., Boucher, O., Harding, R. J., Huntingford, C., and Cox, P. M.: The Joint UK Land Environment Simulator (JULES), model description – Part 2: Carbon fluxes and vegetation dynamics, Geosci. Model Dev., 4, 701–722,, 2011. a

Cunnold, D. M.: In situ measurements of atmospheric methane at GAGE/AGAGE sites during 1985–2000 and resulting source inferences, J. Geophys. Res., 107, 4225,, 2002. 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

Dlugokencky, E. J., Houweling, S., Bruhwiler, L., Masarie, K. A., Lang, P. M., Miller, J. B., and Tans, P. P.: Atmospheric methane levels off: Temporary pause or a new steady-state?, Geophys. Res. Lett., 30, 1992–1995,, 2003. a

Dlugokencky, E. J., Bruhwiler, L., White, J. W., 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 CH 4 burden, Geophys. Res. Lett., 36, 3–7,, 2009. a

Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G., Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando, J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.: Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67,, 2010. a, b

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

Frankenberg, C., Aben, I., Bergamaschi, P., Dlugokencky, E. J., Van Hees, R., Houweling, S., Van Der Meer, P., Snel, R., and Tol, P.: Global column-averaged methane mixing ratios from 2003 to 2009 as derived from SCIAMACHY: Trends and variability, J. Geophys. Res.-Atmos., 116, 1–12,, 2011. a

Ganesan, A. L., Rigby, M., Zammit-Mangion, A., Manning, A. J., Prinn, R. G., Fraser, P. J., Harth, C. M., Kim, K.-R., Krummel, P. B., Li, S., Mühle, J., O'Doherty, S. J., Park, S., Salameh, P. K., Steele, L. P., and Weiss, R. F.: Characterization of uncertainties in atmospheric trace gas inversions using hierarchical Bayesian methods, Atmos. Chem. Phys., 14, 3855–3864,, 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

Ganesan, A. L., Stell, A. C., Gedney, N., Comyn-Platt, E., Hayman, G., Rigby, M., Poulter, B., and Hornibrook, E. R.: Spatially Resolved Isotopic Source Signatures of Wetland Methane Emissions, Geophys. Res. Lett., 45, 3737–3745,, 2018. a

Ganesan, A. L., Schwietzke, S., Poulter, B., Arnold, T., Lan, X., Rigby, M., Vogel, F. R., van der 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. Cycles, 33, 1475–1512,, 2019. a, b

Gumbricht, T., Roman-Cuesta, R. M., Verchot, L., Herold, M., Wittmann, F., Householder, E., Herold, N., and Murdiyarso, D.: An expert system model for mapping tropical wetlands and peatlands reveals South America as the largest contributor, Global Change Biol., 23, 3581–3599,, 2017. a, b

Imbiriba, B. C. O., Ramos, J. R. d. S., de Sousa Silva, R., Cattanio, J. H., do Couto, L. L., and Mitschein, T. A.: Estimates of methane emissions and comparison with gas mass burned in CDM action in a large landfill in Eastern Amazon, Waste Manag., 101, 28–34,, 2020. a

IPCC: 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. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1535 pp., 2013. a

Janardanan, R., Maksyutov, S., Tsuruta, A., Wang, F., Tiwari, Y. K., Valsala, V., Ito, A., Yoshida, Y., Kaiser, J. W., and Janssens-Maenhout, G.: Country-Scale Analysis of Methane Emissions With a High-Resolution Inverse Model Using GOSAT and Surface Observations, Remote Sens.-Basel, 12, 1–27,, 2020. a, b, c, d

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

Jones, A., Thomson, D., Hort, M., and Devenish, B.: The U.K. Met Office's Next-Generation Atmospheric Dispersion Model, NAME III, pp. 580–589,, 2007. a

Lewis, S. L., Brando, P. M., Phillips, O. L., Van Der Heijden, G. M., and Nepstad, D.: The 2010 Amazon drought, Science, 331, 554,, 2011. a

Lunt, M. F., Rigby, M., Ganesan, A. L., and Manning, A. J.: Estimation of trace gas fluxes with objectively determined basis functions using reversible-jump Markov chain Monte Carlo, Geosci. Model Dev., 9, 3213–3229,, 2016. a

Machado, S. L., Carvalho, M. F., Gourc, J. P., Vilar, O. M., and do Nascimento, J. C.: Methane generation in tropical landfills: Simplified methods and field results, Waste Manag., 29, 153–161,, 2009. a

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

Miller, J. B., Gatti, L. V., D'Amelio, M. T., Crotwell, A. M., Dlugokencky, E. J., Bakwin, P., Artaxo, P., and Tans, P. P.: Airborne measurements indicate large methane emissions from the eastern Amazon basin, Geophys. Res. Lett., 34, L10809,, 2007. a

Ministry of Foreign Affairs, Ministry of Science, and Communications: Brazil's Third Biennial Update Report to the United Nations Framework Convention on Climate Change, available at: (last access: 22 October 2020), 2019. a, b, c, d

Ministry of Science and Innovation: Third National Communication of Brazil to the United Nations Framework Convention on Climate Change – Volume I, available at: (last access: 22 October 2020), 2016. a

Monteil, G., Houweling, S., Dlugockenky, E. J., Maenhout, G., Vaughn, B. H., White, J. W. C., and Rockmann, T.: Interpreting methane variations in the past two decades using measurements of CH4 mixing ratio and isotopic composition, Atmos. Chem. Phys., 11, 9141–9153,, 2011. a

Nguyen, N. H., Turner, A. J., Yin, Y., Prather, M. J., and Frankenberg, C.: Effects of Chemical Feedbacks on Decadal Methane Emissions Estimates, Geophys. Res. Lett., 47, 1–13,, 2020. a

Nisbet, E., Fisher, R., Lowry, D., France, J., Allen, G., Bakkaloglu, S., Broderick, T., Cain, M., Coleman, M., Fernandez, J., Forster, G., Griffiths, P., Iverach, C., Kelly, B., Manning, M., Nisbet‐Jones, P., Pyle, J., 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

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., Vaughn, B., Bousquet, P., Pyle, J. A., Warwick, N. J., Cain, M., Brownlow, R., Zazzeri, G., Lanoisellé, M., Manning, A. C., Gloor, E., Worthy, D. E., Brunke, E. G., Labuschagne, C., Wolff, E. W., and Ganesan, A. L.: Rising atmospheric methane: 2007–2014 growth and isotopic shift, Global Biogeochem. Cycles, 30, 1356–1370,, 2016. a, b

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. Cycles, 33, 318–342,, 2019. a, b, c, d

Oki, T., Nishimura, T., and Dirmeyer, P.: Assessment of Annual Runoff from Land Surface Models Using Total Runoff Integrating Pathways (TRIP), J. the Meteorol. Soc. Japan Ser. II, 77, 235–255,, 1999. a

Palmer, P. I., O'Doherty, S., Allen, G., Bower, K., Bösch, H., Chipperfield, M. P., Connors, S., Dhomse, S., Feng, L., Finch, D. P., Gallagher, M. W., Gloor, E., Gonzi, S., Harris, N. R. P., Helfter, C., Humpage, N., Kerridge, B., Knappett, D., Jones, R. L., Le Breton, M., Lunt, M. F., Manning, A. J., Matthiesen, S., Muller, J. B. A., Mullinger, N., Nemitz, E., O'Shea, S., Parker, R. J., Percival, C. J., Pitt, J., Riddick, S. N., Rigby, M., Sembhi, H., Siddans, R., Skelton, R. L., Smith, P., Sonderfeld, H., Stanley, K., Stavert, A. R., Wenger, A., White, E., Wilson, C., and Young, D.: A measurement-based verification framework for UK greenhouse gas emissions: an overview of the Greenhouse gAs Uk and Global Emissions (GAUGE) project, Atmos. Chem. Phys., 18, 11753–11777,, 2018. a, b

Pangala, S. R., Enrich-Prast, A., Basso, L. S., Peixoto, R. B., Bastviken, D., Hornibrook, E. R. C., Gatti, L. V., Ribeiro, H., Calazans, L. S. B., Sakuragui, C. M., Bastos, W. R., Malm, O., Gloor, E., Miller, J. B., and Gauci, V.: Large emissions from floodplain trees close the Amazon methane budget, Nature, 552, 230,, 2017. 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, 2–7,, 2011. a

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

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 CH4observations, Remote Sens. Environ., 211, 261–275,, 2018. a

Peters, W., Jacobson, A. R., Sweeney, C., Andrews, A. E., Conway, T. J., Masarie, K., Miller, J. B., Bruhwiler, L. M., Pétron, G., Hirsch, A. I., Worthy, D. E., Van Der Werf, G. R., Randerson, J. T., Wennberg, P. O., Krol, M. C., and Tans, P. P.: An atmospheric perspective on North American carbon dioxide exchange: CarbonTracker, P. Natl. Acad. Sci. USA, 104, 18925–18930,, 2007. a

Pisso, I., Sollum, E., Grythe, H., Kristiansen, N. I., Cassiani, M., Eckhardt, S., Arnold, D., Morton, D., Thompson, R. L., Groot Zwaaftink, C. D., Evangeliou, N., Sodemann, H., Haimberger, L., Henne, S., Brunner, D., Burkhart, J. F., Fouilloux, A., Brioude, J., Philipp, A., Seibert, P., and Stohl, A.: The Lagrangian particle dispersion model FLEXPART version 10.4, Geosci. Model Dev., 12, 4955–4997,, 2019. a

Pöhlker, C., Walter, D., Paulsen, H., Könemann, T., Rodríguez-Caballero, E., Moran-Zuloaga, D., Brito, J., Carbone, S., Degrendele, C., Després, V. R., Ditas, F., Holanda, B. A., Kaiser, J. W., Lammel, G., Lavrič, J. V., Ming, J., Pickersgill, D., Pöhlker, M. L., Praß, M., Löbs, N., Saturno, J., Sörgel, M., Wang, Q., Weber, B., Wolff, S., Artaxo, P., Pöschl, U., and Andreae, M. O.: Land cover and its transformation in the backward trajectory footprint region of the Amazon Tall Tower Observatory, Atmos. Chem. Phys., 19, 8425–8470,, 2019. a

Prinn, R., Weiss, R., Arduini, J., Arnold, T., DeWitt, H., Fraser, P., Ganesan, A., Gasore, J., Harth, C., Hermansen, O., Kim, J., Krummel, P., Li, S., Loh, Z., Lunder, C., Maione, M., Manning, A., Miller, B., Mitrevski, B., Muhle, J., O'Doherty, S., Park, S., Reimann, S., Rigby, M., Saito, T., Salameh, P., Schmidt, R., Simmonds, P., Steele, L., Vollmer, M., Wang, H. J., Yao, B., Yokouchi, Y., Young, D., and Zhou, L.: History of chemically and radiatively important atmospheric gases from the Advanced Global Atmospheric Gases Experiment (AGAGE), Carbon Dioxide Information Analysis Center (CDIAC), Oak Ridge National Laboratory (ORNL), Oak Ridge, TN, USA,, 2018a. a

Prinn, R. G., Weiss, R. F., Arduini, J., Arnold, T., DeWitt, H. L., Fraser, P. J., Ganesan, A. L., Gasore, J., Harth, C. M., Hermansen, O., Kim, J., Krummel, P. B., Li, S., Loh, Z. M., Lunder, C. R., Maione, M., Manning, A. J., Miller, B. R., Mitrevski, B., Mühle, J., O'Doherty, S., Park, S., Reimann, S., Rigby, M., Saito, T., Salameh, P. K., Schmidt, R., Simmonds, P. G., Steele, L. P., Vollmer, M. K., Wang, R. H., Yao, B., Yokouchi, Y., Young, D., and Zhou, L.: History of chemically and radiatively important atmospheric gases from the Advanced Global Atmospheric Gases Experiment (AGAGE), Earth Syst. Sci. Data, 10, 985–1018,, 2018b. a

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, 2–7,, 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., and Park, S.: Role of atmospheric oxidation in recent methane growth, P. Natl. Acad. Sci. USA, 114, 5373–5377,, 2017. a, b

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, b, c, d, e, f

Schaefer, H., Fletcher, S. E., 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.: A 21st-century shift from fossil-fuel to biogenic methane emissions indicated by 13CH4, ScienceMag,, 2016. a

Schroeder, R., McDonald, K. C., Chapman, B. D., Jensen, K., Podest, E., Tessler, Z. D., Bohn, T. J., and Zimmermann, R.: Development and evaluation of a multi-year fractional surface water data set derived from active/passive microwave remote sensing data, Remote Sens.-Basel, 7, 16688–16732,, 2015. a

Schulz, C., Whitney, B. S., Rossetto, O. C., Neves, D. M., Crabb, L., de Oliveira, E. C., Terra Lima, P. L., Afzal, M., Laing, A. F., de Souza Fernandes, L. C., da Silva, C. A., Steinke, V. A., Torres Steinke, E., and Saito, C. H.: Physical, ecological and human dimensions of environmental change in Brazil's Pantanal wetland: Synthesis and research agenda, Sci. Total Environ., 687, 1011–1027,, 2019. a

Schwietzke, S., Sherwood, O. A., Bruhwiler, L. M., Miller, J. B., Etiope, G., Dlugokencky, E. J., Michel, S. E., Arling, V. A., Vaughn, B. H., White, J. W., and Tans, P. P.: Upward revision of global fossil fuel methane emissions based on isotope database, Nature, 538, 88–91,, 2016. 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

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

Turner, A. J., Frankenberg, C., and Kort, E. A.: Interpreting contemporary trends in atmospheric methane, P. Natl. Acad. Sci. USA, 116, 2805–2813,, 2019. a

UNFCCC: United Nations Framework Convention on Climate Change Paris Agreement, available at: (last access: 22 October 2020), 2015. a

van der Laan-Luijkx, I. R., Krol, M. C., Gatti, L. V., Domingues, L. G., Correia, C. S. C., Miller, J. B., Gloor, M., Leeuwen, T. T., Kaiser, J. W., Wiedinmyer, C., Basu, S., Clerbaux, C., and Peters, W.: Global Biogeochemical Cycles drought derived with CarbonTracker South America, Global Biogeochem. Cycles, 29, 1092–1108,, 2015.  a

van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720,, 2017. 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., 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., 121, 11006–11020,, 2016. a

Wilson, C., Gloor, M., Gatti, L. V., Miller, J. B., Monks, S. A., McNorton, J., Bloom, A. A., Basso, L. S., and Chipperfield, M. P.: Contribution of regional sources to atmospheric methane over the Amazon Basin in 2010 and 2011, Global Biogeochem. Cycles, 30, 400–420,, 2016. a, b, c, d, e, f, g

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

Wunch, D., Toon, G. C., Blavier, J. F. L., Washenfelder, R. A., Notholt, J., Connor, B. J., Griffith, D. W., Sherlock, V., and Wennberg, P. O.: The total carbon column observing network, Philosophical Transactions of the Royal Society A: Mathematical, Phys. Eng. Sci., 369, 2087–2112,, 2011. a

Zhang, Z., Zimmermann, N. E., Calle, L., Hurtt, G., Chatterjee, A., and Poulter, B.: Enhanced response of global wetland methane emissions to the 2015–2016 El Niño-Southern Oscillation event, Environ. Res. Lett., 13, 7,, 2018. a, b, c

Short summary
This study quantifies Brazil’s emissions of a potent atmospheric greenhouse gas, methane. This is in the field of atmospheric modelling and uses remotely sensed data and surface measurements of methane concentrations as well as an atmospheric transport model to interpret the data. Because of Brazil’s large emissions from wetlands, agriculture and biomass burning, these emissions affect global methane concentrations and thus are of global significance.
Final-revised paper