Articles | Volume 13, issue 5
https://doi.org/10.5194/essd-13-2307-2021
https://doi.org/10.5194/essd-13-2307-2021
Review article
 | 
28 May 2021
Review article |  | 28 May 2021

The consolidated European synthesis of CH4 and N2O emissions for the European Union and United Kingdom: 1990–2017

Ana Maria Roxana Petrescu, Chunjing Qiu, Philippe Ciais, Rona L. Thompson, Philippe Peylin, Matthew J. McGrath, Efisio Solazzo, Greet Janssens-Maenhout, Francesco N. Tubiello, Peter Bergamaschi, Dominik Brunner, Glen P. Peters, Lena Höglund-Isaksson, Pierre Regnier, Ronny Lauerwald, David Bastviken, Aki Tsuruta, Wilfried Winiwarter, Prabir K. Patra, Matthias Kuhnert, Gabriel D. Oreggioni, Monica Crippa, Marielle Saunois, Lucia Perugini, Tiina Markkanen, Tuula Aalto, Christine D. Groot Zwaaftink, Hanqin Tian, Yuanzhi Yao, Chris Wilson, Giulia Conchedda, Dirk Günther, Adrian Leip, Pete Smith, Jean-Matthieu Haussaire, Antti Leppänen, Alistair J. Manning, Joe McNorton, Patrick Brockmann, and Albertus Johannes Dolman
Abstract

Reliable quantification of the sources and sinks of greenhouse gases, together with trends and uncertainties, is essential to monitoring the progress in mitigating anthropogenic emissions under the Paris Agreement. This study provides a consolidated synthesis of CH4 and N2O emissions with consistently derived state-of-the-art bottom-up (BU) and top-down (TD) data sources for the European Union and UK (EU27 + UK). We integrate recent emission inventory data, ecosystem process-based model results and inverse modeling estimates over the period 1990–2017. BU and TD products are compared with European national greenhouse gas inventories (NGHGIs) reported to the UN climate convention UNFCCC secretariat in 2019. For uncertainties, we used for NGHGIs the standard deviation obtained by varying parameters of inventory calculations, reported by the member states (MSs) following the recommendations of the IPCC Guidelines. For atmospheric inversion models (TD) or other inventory datasets (BU), we defined uncertainties from the spread between different model estimates or model-specific uncertainties when reported. In comparing NGHGIs with other approaches, a key source of bias is the activities included, e.g., anthropogenic versus anthropogenic plus natural fluxes. In inversions, the separation between anthropogenic and natural emissions is sensitive to the geospatial prior distribution of emissions. Over the 2011–2015 period, which is the common denominator of data availability between all sources, the anthropogenic BU approaches are directly comparable, reporting mean emissions of 20.8 Tg CH4 yr−1 (EDGAR v5.0) and 19.0 Tg CH4 yr−1 (GAINS), consistent with the NGHGI estimates of 18.9 ± 1.7 Tg CH4 yr−1. The estimates of TD total inversions give higher emission estimates, as they also include natural emissions. Over the same period regional TD inversions with higher-resolution atmospheric transport models give a mean emission of 28.8 Tg CH4 yr−1. Coarser-resolution global TD inversions are consistent with regional TD inversions, for global inversions with GOSAT satellite data (23.3 Tg CH4 yr−1) and surface network (24.4 Tg CH4 yr−1). The magnitude of natural peatland emissions from the JSBACH–HIMMELI model, natural rivers and lakes emissions, and geological sources together account for the gap between NGHGIs and inversions and account for 5.2 Tg CH4 yr−1. For N2O emissions, over the 2011–2015 period, both BU approaches (EDGAR v5.0 and GAINS) give a mean value of anthropogenic emissions of 0.8 and 0.9 Tg N2O yr−1, respectively, agreeing with the NGHGI data (0.9 ± 0.6 Tg N2O yr−1). Over the same period, the average of the three total TD global and regional inversions was 1.3 ± 0.4 and 1.3 ± 0.1 Tg N2O yr−1, respectively. The TD and BU comparison method defined in this study can be operationalized for future yearly updates for the calculation of CH4 and N2O budgets both at the EU+UK scale and at the national scale. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.4590875 (Petrescu et al., 2020b).

1 Introduction

The global atmospheric concentration of methane (CH4) has increased by 160 % and that of nitrous oxide (N2O) by 22 % since the pre-industrial period (WMO, 2019), and they are well documented as observed by long-term ice-core records (Etheridge et al., 1998; CSIRO). According to the NOAA atmospheric data (https://www.esrl.noaa.gov/gmd/ccgg/trends_ch4/, last access: June 2020) the CH4 concentration in the atmosphere continues to increase and, after a small dip in 2017, has an average growth of 10 ppb yr−1, representing the highest rate observed since the 1980s1 (Nisbet et al., 2016, 2019). This increase was attributed to anthropogenic emissions from agriculture (livestock enteric fermentation and rice cultivation) and fossil-fuel-related activities, combined with a contribution from natural tropical wetlands (Saunois et al., 2020; Thompson et al., 2018; Nisbet et al., 2019). The recent increase in atmospheric N2O is more linked to agriculture in particular due to the application of nitrogen fertilizers and livestock manure on agricultural land (FAO, 2020, 2015; IPCC, 2019b; Tian et al., 2020).

National greenhouse gas inventories (NGHGIs) are prepared and reported on an annual basis by Annex I countries2 based on IPCC Guidelines using national activity data and different levels of sophistication (tiers) for well-defined sectors. These inventories contain annual time series of each country's greenhouse gas (GHG) emissions from the 1990 base year3 until 2 years before the year of reporting and were originally set to track progress towards their reduction targets under the Kyoto Protocol (UNFCCC, 1997). Non-Annex I countries provide some information in biennial update reports (BURs) as well as national communications (NCs), but neither BURs nor NCs report annual time series or use harmonized formats. The IPCC tiers represent the level of sophistication used to estimate emissions, with Tier 1 based on global or regional default values, Tier 2 based on country-specific parameters, and Tier 3 based on more detailed process-level modeling. Uncertainties in NGHGIs are calculated based on ranges in observed (or estimated) emission factors and variability of activity data, using the error propagation method (95 % confidence interval) or Monte Carlo methods, based on clear guidelines (IPCC, 2006).

NGHGIs follow principles of transparency, accuracy, consistency, completeness and comparability (TACCC) under the guidance of the UNFCCC (UNFCCC, 2014). Methodological procedures are taken from the 2006 IPCC Guidelines (IPCC, 2006). The IPCC 2019 Refinement (IPCC, 2019a), which may be used to complement the 2006 IPCC Guidelines, has updated sectors with additional sources and provides guidance on the possible and voluntary use of atmospheric data for independent verification of GHG inventories. Complementary to NGHGIs, research groups and international institutions produce estimates of national GHG emissions, with two families of approaches: atmospheric inversions (top-down, TD) and GHG inventories based on the same principle as NGHGIs but using different methods and input data (bottom-up, BU). These complementary approaches are necessary. First, TD approaches act as an independent check on BU approaches and facilitate a deeper understanding of the scientific processes driving different GHG budgets. Second, NGHGIs only cover a subset of countries, and it is therefore necessary to construct BU estimates independently for all countries. The BU estimates are often used as input data for TD estimates and to track emissions over time, either globally or on the country level, such as in the UNEP Emissions Gap Report (UNEP, 2019). There is no guideline to estimate uncertainties in TD or BU approaches. The uncertainties are usually assessed from the spread of different estimates within the same approach, though some groups or institutions report uncertainties for their estimates using a variety of methods – for instance, by varying parameters or input data. However, this gets complicated when dealing with complex process-based models.

NGHGI official numbers are not always straightforward to compare with other independent estimates. Independent estimates often have different system boundaries and a different focus. BU estimates often have a lot of overlap in terms of methods and other input data, and through harmonization the differences between BU estimates and NGHGIs can be bridged. On the other hand, TD estimates are much more independent and provide the best independent check on NGHGIs. While NGHGIs go through a review process, the UNFCCC procedures do not incorporate mandatory independent, large-scale observation-derived verification but allow the use of atmospheric data for external checks within the data quality control, quality assurance and verification process (2006 IPCC Guidelines, chap. 6: QA/QC procedures). So far, only a few countries (e.g., Switzerland, the UK, New Zealand and Australia) have used atmospheric observations (TD) to complement their national inventory data (Bergamaschi et al., 2018a).

A key priority in the current policy process is to facilitate the global stocktake exercise of the Paris Agreement, the first one coming in 2023, and to assess collective progress towards achieving the near- and long-term objectives, considering mitigation, adaptation and means of implementation. The global stocktake is expected to create political momentum for enhancing commitments in nationally determined contributions (NDCs) under the Paris Agreement. Key components of the global stocktake are the NGHGIs submitted by countries under the enhanced transparency framework of the Paris Agreement. Under the framework, for the first time, developing countries will be required to submit their inventories and also commit to provide regular reports to UNFCCC, alongside developed countries, which will continue to submit also on an annual basis. Some developing countries will face challenges to provide and then update inventories.

The work presented here represents dozens of distinct datasets and models, in addition to the individual country submissions to the UNFCCC for all European countries (NGHGIs), which while following the general guidance laid out in IPCC (2006) still differ in specific approaches, models and parameters, in addition to differences underlying activity datasets. A comprehensive investigation of detailed differences between all datasets is beyond the scope of this paper, though attempts have been previously made for specific subsectors (e.g., agriculture, Petrescu et al., 2020) and in dedicated gas-specific follow-ups to this paper. As this is the most comprehensive comparison of NGHGIs and research datasets (including both TD and BU approaches) for the European continent to date, we focus here on the following rich set of questions that such a comparison raises without necessarily offering detailed solutions. How can we compare the detailed sectoral NGHGI to the observation-based estimates? Which new information are the observation-based estimates likely to bring (mean fluxes, trend, ensemble variability)? What can we expect from such a complex study and how do we proceed going forward?

Table 1Sectors used in this study and data sources providing estimates for these sectors.

a For consistency with the NGHGIs, here we refer to the five reporting sectors as defined by the UNFCCC and the Paris Agreement decision (18/CMP.1), the IPCC Guidelines (IPCC, 2006) and their 2019 Refinement (IPCC, 2019a), with the only exception that the latest IPCC refinement groups together the agriculture and LULUCF sectors in one sector (agriculture, forestry and other land use – AFOLU). b The term “natural” here refers to unmanaged natural CH4 emissions (wetlands, geological, inland waters) not reported under the UNFCCC LULUCF sector. c Anthropogenic (managed) agricultural soils can also have a level of natural emissions. d Natural soils (unmanaged) can have both natural and anthropogenic emissions.

Download Print Version | Download XLSX

Table 2Data sources for CH4 and N2O emissions used in this study.

Download XLSX

We compare official anthropogenic NGHGI emissions with research datasets and wherever needed harmonize research data on total emissions to ensure consistent comparisons of anthropogenic emissions. We analyze differences and inconsistencies between emissions and make recommendations towards future actions to evaluate NGHGI data. While NGHGIs include uncertainty estimates, individual spatially disaggregated research datasets of emissions often lack quantification of uncertainty. Here, we use the median4 and minimum/maximum (min/max) range of different research products of the same type to get a first estimate of uncertainty.

2CH4 and N2O data sources and estimation approaches

We analyze CH4 and N2O emissions in the EU27 + UK from inversions (TD) and anthropogenic emissions from various BU approaches that cover specific sectors. These data (Table 2) span the period from 1990 and 2017, with the same data available for shorter time periods. The data are from peer-reviewed literature and from unpublished research results from the VERIFY project (Table 1 and Appendix A). They are compared with NGHGI official submissions up to 2017 and supplemented by the UNFCCC-NRT inventory to capture 2018 estimates (near real time, EEA 2019). References are given in Table 2 and the detailed description of all products in Appendices A1–A3.

For both CH4 and N2O BU approaches, we used inventories of anthropogenic emissions covering all sectors (EDGAR v5.0 and GAINS) and inventories limited to agriculture (CAPRI and FAOSTAT). For CH4 we used one biogeochemical model of natural peatland emissions (JSBACH–HIMMELI), as well as literature data for geological emissions on land (excluding marine seepage) (Etiope et al., 2019; Hmiel et al., 2020) and for lakes and reservoirs (Del Sontro et al., 2018). Emissions from gas hydrates and termites are not included as they are close to zero in the EU27 + UK (Saunois et al., 2020). Biomass burning emissions of CH4 from the land use, land use change and forestry (LULUCF) sector account for 3 % of the total emissions in the EU27 + UK. These estimates are described in Sect. 2.2. From TD approaches, we used both regional and global inversions, with the latter having a coarser spatial resolution. These estimates are described in Sect. 2.3.

For N2O emissions, we used the same global BU inventories as for CH4 and natural emissions from inland waters (rivers, lakes and reservoirs) (Maavara et al., 2019; Lauerwald et al., 2019). In this study, about 66 % of the N2O emitted by Europe's natural rivers is considered anthropogenic indirect emissions, caused by leaching and runoff of N fertilizers from the agriculture sector. We did not account for natural N2O emissions from unmanaged soils (Tian et al., 2019, estimated pre-industrial soil emissions in Europe at a third of the level of the most recent decade – emissions that in pre-industrial times may have been influenced by human management activities or based on natural processes that have been abolished since). For N2O inversions, we used one regional inversion (FLEXINVERT_NILU) and three global inversions (Friedlingstein et al., 2019; Tian et al., 2020). Agricultural sector emissions of N2O were presented in detail by Petrescu et al. (2020). In this current study these emissions come from the CAPRI model and FAOSTAT, with the latter additionally covering non-CO2 emissions from biomass fires in LULUCF. Fossil-fuel-related and industrial emissions were obtained from GAINS (see Appendix A1). Table A2 in Appendix A presents the methodological differences of the current study with respect to Petrescu et al. (2020).

The units used in this paper are metric tonnes (t) (1 kt= 109g; 1 Mt= 1012g) of CH4 and N2O. The referenced data used for the figures' replicability purposes are available for download at https://doi.org/10.5281/zenodo.4590875 (Petrescu et al., 2020). Upon request, we can provide the codes necessary to plot precisely the style and/or layout of the figures. We focus herein on the EU27 + UK. In the VERIFY project, we have constructed in addition a web tool which allows for the selection and display of all plots shown in this paper (as well as the companion paper on CO2, Petrescu et al., 2021) not only for the regions but for a total of 79 countries and groups of countries in Europe. The data, located on the VERIFY project website (http://webportals.ipsl.jussieu.fr/VERIFY/FactSheets/, last access: February 2021) are free but accessible with a username and password distributed by the project. Figure 1 includes also data from countries outside the EU but located within geographical Europe (Switzerland, Norway, Belarus, Ukraine and the Republic of Moldova).

2.1CH4 and N2O anthropogenic emissions from NGHGIs

UNFCCC NGHGI (2019) emissions are country estimates covering the period 1990–2017. They were kept separate to be compared with other BU and TD data. We supplemented the NGHGI estimates with the member state (MS) near-real-time inventory (MS-NRT; EEA, 2019) to capture one additional year with preliminary estimates5. MS-NRT represents the approximated GHG inventory (also referred to as “proxy estimates”) with an early estimate of the GHG emissions for the preceding year, as required by Regulation (EU) 525/2013 of the European Parliament and of the Council.

Anthropogenic CH4 emissions from the four UNFCCC sectors (Table 1, excluding LULUCF) were grouped together. As anthropogenic NGHGI CH4 emissions from the LULUCF sector are very small for the EU27 + UK (2.6 % in 2017 including biomass burning), we exclude them in Fig. 4 but include them in the total UNFCCC estimates from Figs. 1–3, 5 and 6. Only a few countries6 under the NGHGIs volunteered to report wetland emissions, following the recommendations of the 2014 IPCC wetlands supplement (IPCC, 2014), and these emissions were not included in the NGHGI total, following the IPCC (2006) Guidelines as the reference for NGHGIs and in the absence of a detailed description of what they cover. According to NGHGI data between 2008 and 2017, the wetland emissions in the EU27 + UK reported under LULUCF (CRF Table 4(II) accessible for each EU27 + UK country7) include only managed wetlands which represent one-fourth of the total wetland area in the EU27 + UK (Giacomo Grassi, personal communication, 2019) and sum up to 0.1 Tg of CH4 (Petrescu et al., 2020a).

Anthropogenic N2O emissions are predominantly related to agriculture (for the EU27 + UK, 69 % in 2017) but are also found in the other sectors (Tian et al., 2020). In addition, N2O has natural emissions, which are defined as the pre-industrial background, that is, before the use of synthetic N fertilizers and intensive agriculture, and derive from natural processes in soils but also in lakes, rivers and reservoirs (Maavara et al., 2019; Lauerwald et al., 2019; Tian et al., 2020).

2.2CH4 and N2O anthropogenic and natural emissions from other bottom-up sources

We used four global CH4 and N2O BU anthropogenic emissions inventories: CAPRI, FAOSTAT, GAINS and EDGAR v5.0 (Tables 2 and 3). These estimates are not completely independent from NGHGIs (see Fig. 4 in Petrescu et al., 2020a) as they integrate their own sectorial modeling with the UNFCCC data (e.g., common activity data and IPCC emission factors) when no other source of information is available.

Anthropogenic emissions from these datasets follow and can be matched to Table 1 sectors. The CH4 biomass and biofuel burning emissions are included in NGHGIs under the UNFCCC LULUCF sector, although they are identified as a separate category by the Global Carbon Project CH4 budget synthesis (Saunois et al., 2020). For both CH4 and N2O, CAPRI (Britz and Witzke, 2014; Weiss and Leip, 2012) and FAOSTAT (FAO, 2020) report only agricultural emissions. None of the BU inventories reported uncertainties, except for the 2015 values of EDGAR v5.0 (Solazzo et al., 2021) and for an earlier FAOSTAT dataset only up to 2010 (Tubiello et al., 2013, and Appendix B).

The CH4 natural emissions belong to “peatlands” and “other natural emissions”, with the latter including geological sources and inland waters (lakes and reservoirs), following Saunois et al. (2020). For peatlands, we used the JSBACH–HIMMELI framework and the ensemble of 13 monthly gridded estimates of peatland emissions based on different land surface models as calculated for Saunois et al. (2020), all described in Appendix A2. In EU27 + UK, geological emissions were calculated by scaling up the regional emissions from Etiope et al. (2019) (37.4 Tg CH4 yr−1) to the global ratio of emissions from Hmiel et al. (2020) (see Appendix A2, geological fluxes), obtaining an estimate of 1.3 Tg CH4 yr−1 (marine and land geological). Marine seepage emissions were excluded. This rescaled geological source represents 24 % of the total EU27 + UK natural CH4 emissions. Inland waters (lakes and reservoirs, based on Lauerwald et al., 2019, and Del Sontro et al., 2018) (Appendix A2) are the largest natural component (48 %), with the rest (28 %) being attributed to peatlands. Overall, in the EU27 + UK the natural emissions thus accounted for 5.2 Tg CH4 yr−1.

The N2O anthropogenic emissions from BU datasets belong predominantly to two main categories, as presented in Table 2: (1) direct emissions from the agricultural sector where synthetic fertilizers and manure were applied, as well as from manure management; and (2) indirect emissions on non-agricultural land and water receiving anthropogenic N through atmospheric N deposition, leaching and runoff (also from agricultural land). Furthermore, emissions from industrial processes are declining over time but originate from fossil fuel combustion, air pollution abatement devices, specific chemical reactions, wastewater treatment and land use change. In this study, we do not consider the natural emissions from soils, since these emissions are relatively small for temperate regions, including Europe, and cannot be singled out in landscapes largely dominated by human activities. Therefore, the only “natural” fluxes considered in this study are emissions from inland waters (lakes, rivers and reservoirs, Maavara et al., 2019; Lauerwald et al., 2019, Appendix A3) even if more than half of the emissions (56 % globally, Tian et al., 2020, and 66 % for Europe this study) are due to eutrophication following N fertilizer leaching to inland waters. Emissions from natural soils in this study are considered anthropogenic because, according country-specific national inventory reports (NIRs), all land in the EU27 + UK is considered to be managed, except 5 % of France's EU territory.

2.3CH4 and N2O emission data from top-down inversions

Inversions combine atmospheric observations, transport and chemistry models, and estimates of GHG sources all with their uncertainties to estimate emissions. Emission estimates from inversions depend on the dataset of atmospheric measurements and the choice of the atmospheric model, as well as on other settings (e.g., prior emissions and their uncertainties). Inversion outputs were taken from original publications without evaluation of their performance through specific metrics (e.g., fit to independent cross-validation atmospheric measurements (Bergamaschi et al., 2013, 2018; Patra et al., 2016). Some of the inversions solve explicitly for sectors; others solve for all fluxes in each grid cell and separate sectors using prior grid-cell fractions (see details in Saunois et al., 2020, for global inversions).

For CH4, we use 9 regional TD inversions and 22 global TD inversions listed in Table 2. These inversions are not independent from each other: some are variants from the same modeling group, many use the same transport model and most of them use the same atmospheric data. Different prior data are generally used in models, which produce a greater range of posterior emission estimates (Appendix B3, Table B4). The subset of InGOS inversions (Bergamaschi et al., 2018a) belongs to a project where all models used the same atmospheric data over Europe covering the period 2006–2012. The global inversions from Saunois et al. (2020) were all updated to 2017.

The regional inversions generally use both higher-resolution a priori data and higher-resolution transport models, and, e.g., TM5-JRC runs simultaneously over the global domain at coarse resolution and over the European domain at higher resolution, with atmospheric CH4 concentration boundary conditions taken from global fields. For CH4, 11 global inversions use GOSAT for the period 2010–2017, 8 global inversions use surface stations (SURF) since 2000, 2 global inversions use SURF since 2010 and 1 global inversion uses SURF since 2003 (see Appendix 4 Table in Saunois et al., 2020, and Table 2 below). All regional inversions use observations from SURF stations as a base of their emission calculation.

For N2O, we use one regional inversion (FLEXINVERT_NILU for 2005–2017 period) and three global inversions for the period 1998–2016 from Thompson et al. (2019), listed in Table 2. These inversions are not completely independent from each other since most of them use the same input information (Appendix B, Table B4). The regional inversion uses a higher-resolution transport model for Europe, with atmospheric N2O concentration boundary conditions taken from global fields. As all inversions derived total rather than anthropogenic emissions; emissions from inland waters (lakes, rivers and reservoirs) estimated by Maavara et al. (2019) and Lauerwald et al. (2019) were subtracted from the total emissions. Note that the estimates of Maavara et al. (2019) and Lauerwald et al. (2019) include anthropogenic emissions from N fertilizer leaching accounting for 66 % of the inland water emissions in the EU27 + UK. In 2016, emissions from rivers represent 2.2 % of the total UNFCCC NGHGI (2019) N2O emissions. The natural N2O emissions are small but should be better quantified in the future to allow for a more accurate comparison between BU (anthropogenic sources only) and TD estimates.

The largest share of N2O emissions comes from the agricultural soils (direct and indirect emissions from the applications of fertilizers, whether synthetic or manure) contributing in 2017 69 % of the total N2O emissions (excluding LULUCF) in the EU27 + UK. Table B3 in Appendix B1 presents the allocation of emissions by activity type covering all agricultural activities and natural emissions, following the IPCC classification. We notice that each data product has its own particular way of grouping emissions and does not necessarily cover all emissions activities. Main inconsistencies between models and inventories are observed with activity allocation in the two models (ECOSSE and DayCent). ECOSSE only estimates direct N2O emissions and does not estimate downstream emissions of N2O, for example indirect emissions from nitrate leached into water courses, which also contributes to an underestimation of the total N2O emissions. Field burning emissions are also not included by most of the data sources.

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f01

Figure 1Five-year-average (2011–2015) total CH4 emission estimates (including LULUCF) for the EU27 + UK and five European regions (north, west, central, south and east non-EU). The eastern European region does not include European Russia, and the UNFCCC uncertainty for the Republic of Moldova was not available. Northern Europe includes Norway. Central Europe includes Switzerland. The data belong to the UNFCCC NGHGI (2019) submissions (grey) and base year 1990 (black star), two BU inventories (GAINS and EDGAR v5.0), natural unmanaged emissions (sum of peatland, geological and inland water emissions), and three TD total estimates (regional European inversions – excluding InGOS unavailable for 2013–2015 – and GOSAT and SURF estimates from global inverse models). The relative error on the UNFCCC value represents the UNFCCC NGHGI (2018) reported uncertainty computed with the error propagation method (95 % confidence interval): 9.3 % for the EU27 + UK, 10 % for eastern Europe non-EU, 7.8 % for northern Europe, 10.9 % for southern Europe, 16.1 % for western Europe and 11 % for central Europe. The uncertainty for EDGAR v5.0 was calculated for 2015 and represents the 95 % confidence interval of a lognormal distribution.

3 Results and discussion

3.1 Comparing CH4 emission estimates from different approaches

3.1.1 Estimates of European and regional total CH4 fluxes

We present results of the total CH4 fluxes from the EU27 + UK and five main regions in Europe: north, west, central, east (non-EU) and south. The countries included in these regions are listed in Appendix A, Table A1. Figure 1 shows the total CH4 fluxes from NGHGIs for both the base year 1990 and mean of the 2011–2015 period. This period was the common denominator for which data were available, including 2 years of the Kyoto Protocol first reporting period (2011/12) and reaching the year of the Paris Agreement was adopted. We aim with the selection of this period to bring together all information over a 5-year period for which values are known in 2018. In fact, this can be seen as a reference for what we can achieve in 2023, the year of the first global stocktake, where for most UN Parties the reported inventories will be known until 2021. Given that the global stocktake is only repeated every 5 years, a 5-year average is clearly of interest.

The total NGHGI estimates include emissions from all sectors, and we plot and compare them with fluxes from global datasets, BU models and inversions. We note that for all five regions, the NGHGI-reported CH4 emissions decreased by 21 % in southern Europe, by up to 54 % in eastern Europe and by 35 % for the European Union with respect to the 1990 value. This is encouraging in the context of meeting EU commitments under the PA (at least 50 % and towards 55 % compared with 1990 levels stated by the amended proposal for a regulation of the European Parliament and of the council on establishing the framework for achieving climate neutrality and amending Regulation (EU) 2018/1999 (European Climate Law) (https://ec.europa.eu/clima/sites/clima/files/eu-climate-action/docs/prop_reg_ecl_en.pdf, last access: November 2020) and reaching carbon neutrality by 2050). It also shows that the emissions from BU (anthropogenic and natural) and TD estimates agree well with reported NGHGI data despite the high uncertainty observed in the TD models not only at the EU27 + UK level but also at the regional European level. This uncertainty is represented here by the variability in the model ensembles and denotes the range of the extremes (min and max) of estimates within each model ensemble. From Fig. 1 we clearly note that northern Europe is dominated by natural (wetlands) emissions while western, central and southern European emissions are dominated by anthropogenic sectors (e.g., agriculture).

The EDGAR v5.0 estimate for northern Europe is twice as high when compared to NGHGI and GAINS, and this is because of CH4 emissions from the fuel production and distribution (IPCC sector 1B) and waste sectors. Most Scandinavian countries rely for their power and heat supply on biogenic fuels, which introduces more uncertainty in the use of activity data and emission factors. The allocation of auto-producers as explained in Sect. 3.2 could be another reason for differences. The waste sector emissions for Norway, Sweden, Finland and Estonia are different but still consistent with the landfill emissions from EDGAR v4.3.2, which are known to be up to twice as high as the nationally reported value (Janssens-Maenhout et al., 2019). For eastern Europe we note that BU anthropogenic estimates have the same magnitude as the TD. One possible explanation is linked to the fact that for TD estimates (i.e., using atmospheric inversions) the fluxes are strongly constrained by the density of observations. Where there are few or no observations, the fluxes in the inversion will stay close to the prior estimates, since there is little or no information to adjust them.

In line with Bergamaschi et al. (2018a) we highlight the potential significant contribution from natural unmanaged sources (peatlands, geological and inland water), which for the EU27 + UK accounted for 5.24 Tg CH4 yr−1 (Fig. 1). Taking into account these natural unmanaged CH4 emissions and adding the to the range of the anthropogenic estimates (19–21 Tg CH4 yr−1), the total BU estimates become broadly consistent for all European regions with the range of the TD estimates (23–28 Tg CH4 yr−1).

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f02

Figure 2The contribution of changes (%) in CH4 anthropogenic emissions in the five UNFCCC sectors to the overall change in decadal mean, as reported to UNFCCC NGHGI (2019). The three stacked columns represent the average CH4 emissions from each sector during three periods (1990–1999, 2000–2009 and 2010–2017), and percentages represent the contribution of each sector to the total reduction percentages (black arrows) between periods.

Download

3.1.2 NGHGI sectoral emissions and decadal changes

According to the UNFCCC (2019) NGHGI estimates, in 2017 the EU27 + UK emitted GHGs totaling 3.9 Gt CO2 eq. (including LULUCF); of this total, CH4 emissions accounted for  11 % (0.4 Gt CO2 eq. or 18.1 Mt CH4 yr−1) (Appendix B2, Fig. B1a) with France, the UK and Germany contributing together 36 % of the total CH4 emissions.

The data in Fig. 2 show anthropogenic CH4 emissions and their change from one decade to the next, from UNFCCC NGHGI (2019), with the contribution from different UNFCCC sectors. In 2017, NGHGIs report CH4 from agricultural activities to be 52 % (± 10 %) of the total EU27 + UK CH4 emissions, followed by emissions from waste, 27 % (± 23 %). The large share of agriculture in total anthropogenic CH4 emissions also holds at the global level (IPCC, 2019a). Between the 1990s and the 2000s, the net 17.7 % reduction originates largely from energy and waste, with IPPU (metal and chemical industry) and LULUCF having negligible change. Between the 2000s and 2010–2017, the 15.5 % reduction is distributed more evenly across sectors, with waste having the largest contribution and industry showing no change. The two largest sectors composing total EU27 + UK emission are the agriculture and waste sectors, but the energy and waste sectors have shown higher reductions over the last decade.

The reduction observed in the waste sector is partly due to the adoption of the first EU methane strategy published in 1996 (COM(96), 1996). EU legislation addressing emissions in the waste sector proved to be successful and brought about the largest reductions. Directive 1999/31/EC on the landfill of waste (also referred to as the Landfill Directive) required the members states (MSs) to separate waste, minimizing the amount of biodegradable waste disposed untreated in landfills, and to install landfill gas recovery at all new sites. Based on the 1999 directive, the new 2018/1999 EU regulation on the governance of the Energy Union requires the European Commission to propose a strategic plan for methane, which will become an integral part of the EU's long-term strategy. In the waste sector, the key proposal included the adoption of EU legislation requiring the installation of methane recovery and use systems at new and existing landfills. Other suggested actions included measures aimed at the minimization, separate collection and material recovery of organic waste (Olczak and Piebalgs, 2019).

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f03

Figure 3Total anthropogenic CH4 emissions (excluding LULUCF): (a) of the EU27 + UK and total sectoral emissions as (b) energy, (c) IPPU, (d) agriculture and (e) waste from UNFCCC NGHGI (2019) submissions and MS-NRT 2018 compared to global bottom-up inventory models for agriculture (CAPRI, FAOSTAT) and all sectors excluding LULUCF (EDGAR v5.0, GAINS). CAPRI reports one estimate for Belgium and Luxembourg. The relative error on the UNFCCC value represents the UNFCCC NGHGI (2018) MS-reported uncertainty computed with the error propagation method (95 % confidence interval): 9.4 % for the total EU27 + UK, 23 % for energy and waste, 27 % for IPPU and 10 % for agriculture. The uncertainty for EDGAR v5.0 was calculated for 2015, and the min/max values for all sectors are as follows: EU27 + UK total 15/16, energy 33/37, IPPU 39/34, agriculture 18/18 and waste 32/38; it represents the 95 % confidence interval of a lognormal distribution. The mean values on the right-hand side reflect the values for the common overlapping period 1990–2015. The last reported year in this study refers to 2017 (UNFCCC and FAOSTAT), 2015 (EDGAR v5.0 and GAINS) and 2013 (CAPRI).

Download

3.1.3 NGHGI estimates compared with bottom-up inventories

The data in Fig. 3 present the total anthropogenic CH4 emissions from four BU inventories and UNFCCC NGHGI (2019) excluding those from LULUCF. According to NGHGIs, anthropogenic emissions from the total EU27 + UK of the four UNFCCC sectors (Table 1, excluding LULUCF) amounted to 18.2 Tg of CH4 in the year 2017, which is 10.7 % of the total GHG emissions in CO2 eq. In Fig. 3a, we observe that EDGAR v5.0 and GAINS show consistent trends with NGHGIs (excluding LULUCF), but GAINS reports consistently lower estimates (10 %) and EDGAR v5.0 consistently higher estimates (8 %) compared to NGHGIs. In contrast to the previous version, EDGAR v4.3.2, which was found by Petrescu et al. (2020a) to be consistent with NGHGI (2018) data, EDGAR v5.0 reports higher estimates but within the 9.4 % UNFCCC uncertainty range. The trends in emissions agree better between the two BU inventories and NGHGIs over 1990–2015, with linear trends of 1.5 % yr−1 in NGHGIs compared to 1.5 % yr−1 in GAINS and 1.4 % yr−1 in EDGAR v5.0.

Sectoral time series of anthropogenic CH4 emissions (excluding LULUCF) and their means are shown in Fig. 3b–e. For the energy sector (Fig. 3b), both EDGAR v5.0 and GAINS match the NGHGI trend well thanks to updated methodology that derives bottom-up emission factors and accounts for country-specific information about associated petroleum gas generation and recovery, venting, and flaring (Höglund-Isaksson, 2017). After 2005, GAINS reports consistently lower emissions than UNFCCC due to a phasedown of hard coal production in the Czech Republic, Germany, Poland and the UK; a decline in oil production in particular in the UK; and declining emission factors reflecting reduced leakage from gas distribution networks as old town gas networks are replaced. The consistently higher estimates (+6 % compared to the UNFCCC mean) of EDGAR v5.0 might be due to the use of default emission factors for oil and gas production based on data from the US (Janssens-Maenhout et al., 2019). Next to that, several other reasons could be the cause for the differences (e.g., use of Tier 1 emission factors for coal mines, assumptions for material in the pipelines (in the case of gas transport) and the activity data). EDGAR v5.0, for example, uses the gas pipeline length as a proxy for the activity data; however, this may not be appropriate for the case of the official data, which could consider the total amount of gas being transported or both methods according to the countries. Using pipeline length may overestimate the emissions because the pipeline is not always at 100 % capacity; thus, a larger amount of methane is assumed to be leaked. For coal mining, emissions are a function of the different types of processes being modeled.

The IPPU sector (Fig. 3c), which has only a small share of the total emissions, is not reported in GAINS, while EDGAR v5.0 estimates are less than half of the emissions reported by NGHGI 2019 in this sector. The discrepancy for this sector has a negligible impact on discrepancy for the total CH4 emission. However, we identified that the low bias of EDGAR v5.0 could be explained by fewer activities included in EDGAR v5.0 (e.g., missing solvent, electronics and other manufacturing goods), accounting for 5.5 % of the total IPPU emissions in 2015 reported to UNFCCC. The reason for the remaining difference could be explained by the allocation of emissions from auto-producers8 in EDGAR v5.0 to the energy sector (following the IPCC 1996 guidelines), while in NGHGIs they are reported under the IPPU sector (following the 2006 IPCC Guidelines).

As CAPRI and FAOSTAT just report emissions from agriculture, we only included them in Fig. 3d. The data show that the four data sources (EDGAR v5.0, GAINS, CAPRI and FAOSTAT) show good agreement, with CAPRI at the lower range of emissions (Petrescu et al., 2020a) and on average 3 % lower than UNFCCC and EDGAR v5.0 at the upper range. The reason for EDGAR v5.0 having the highest estimate (contrary to Petrescu et al., 2020a, where NGHGIs were the highest and EDGAR v4.3.2 the second highest) is likely due to the activity data updates in EDGAR v5.0 based on FAOSTAT values, compared to EDGAR v4.3.2. When looking at the time series mean, EDGAR v5.0, GAINS and FAOSTAT show a similar value, +5 % higher than the NGHGI. This shows good consistency between the three BU estimates and UNFCCC likely due to the use of similar activity data and emission factors (EFs); cf. Fig. 4 in Petrescu et al. (2020a).

For the waste sector (Fig. 3e) EDGAR v5.0 shows consistent higher estimates compared to the NGHGI data, while GAINS emissions have an increasing trend after 2000 (mean 1990–2015 value 6 % higher than NGHGIs). The two inventories, EDGAR v5.0 2020 update for landfills and GAINS, used an approach based on the decomposition of waste into different biodegradable streams, with the aim of applying the methodology described in the 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories and the IPCC waste model (IPCC, 2019b) using the first-order-decay (FOD) method. The main differences between the two datasets come from (i) sources for total waste generated per person, (ii) assumption for the fraction composted and (iii) the oxidation. The two inventories may have used different strategies to complete the waste database when inconsistencies were observed in the EUROSTAT database or in the waste trends in UNFCCC.

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f04

Figure 4(a) Comparison of the total CH4 emissions from top-down regional inversions with UNFCCC NGHGI (2019) data and inland water (lakes_reservoirs_ULB, pink), peatland (from JSBACH–HIMMELI, green) and geological emissions (yellow); (b) comparison of anthropogenic CH4 emissions from top-down regional inversions with UNFCCC NGHGI (2019) data. Anthropogenic emissions from these inversions are obtained by removing natural emissions shown in Fig. 4a. The MS-NRT LULUCF estimate does not include the following countries: Austria, Belgium, Estonia, Croatia, Hungary, Luxembourg, Latvia, Malta and Slovenia. UNFCCC NGHGI (2018) reported the uncertainty computed with the error propagation method (95 % confidence interval) is 9.29 % and represents the UNFCCC NGHGI (2018) MS-reported uncertainty for all sectors (including LULUCF). The time series mean was computed for the common period 2006–2012.

Download

3.1.4 NGHGI estimates compared to atmospheric inversions

Regional inversions

Figure 4 compares TD regional estimates with NGHGI anthropogenic data for CH4 and with natural BU emissions. We present TD estimates of the total emissions (Fig. 4a) as well as estimates of anthropogenic emissions only (Fig. 4b), which are calculated by subtracting the natural emissions from the total inversions.

The TD estimates of European CH4 emissions of Fig. 4 use four European regional models (2005–2017) and an ensemble of five different inverse models (InGOS, Bergamaschi et al., 2015) for 2006–2012.

For the common period 2006–2012, the four inverse models give a total CH4 emissions mean of 25.8 (24.0–27.4) Tg CH4 yr−1 compared to anthropogenic total of 20.3 ± 1.9 Tg CH4 yr−1 in NGHGIs (Fig. 4a). The large positive difference between TD and NGHGIs suggests a potentially significant contribution from natural sources (peatlands, geological sources and inland waters), which for the same period report a total mean of 5.2 Tg CH4 yr−1. However, it needs to be emphasized that wetland emission estimates have large uncertainties and show large variability in the spatial (seasonal) distribution of CH4 emissions, but for Europe their inter-annual variability is not very strong (mean of 13 years from JSBACH–HIMMELI peatland emissions 1.4 ± 0.1 Tg CH4 yr−1). Overall, they do represent an important source and could dominate the budget assessments in some regions such as northern Europe (Fig. 1). We also note that the TD trends do not necessarily match those of NGHGIs, and this might be due to strong seasonality of emissions coming from the natural fluxes used as input to the inversions (Saunois et al., 2020).

The natural emissions from inland waters (based on Lauerwald et al., 2019; see Appendix A2) contribute 2.53 Tg CH4 yr−1, or 48 % of the total natural CH4 emissions (sum of lake and reservoir, geological, and peatland emissions). Peatlands (Raivonen et al., 2017, and Susiluoto et al., 2018) account for 1.38 Tg CH4 yr−1, i.e., 27 % of the total natural CH4 emissions, and geological sources sum up to 1.27 Tg CH4 yr−1, i.e., 25 % of the total natural CH4 emissions. It should be noted that geological emissions are an important component of the EU27 + UK emissions budget, although not of concern for climate warming if their source strength has not changed since pre-industrial times (Hmiel et al., 2020). According to the 2006 IPCC Guidelines (IPCC, 2006) CH4 emissions from wetlands are reported by the MS to the NGHGIs under the LULUCF sector and considered anthropogenic. They are included in the total LULUCF values (Figs. 1, 2, 4 and 5), and in 2017 only eight EU countries (Germany, Denmark, Spain, Estonia, Finland, Ireland, Latvia and Sweden) reported CH4 emissions from wetlands, accounting only for 11.2 kt CH4 yr−1.

In an attempt to quantify the anthropogenic CH4 component in the European TD estimates, in Fig. 4b we subtract from the total TD emissions the BU peatland emissions from the regional JSBACH–HIMMELI model and those from geological and inland water sources. It remains however uncertain to perform these corrections due to the prior inventory data allocation of emissions to different sectors (e.g., anthropogenic or natural), which can induce uncertainty of up to 100 % if for example an inventory allocates all emissions to natural emissions and the correction is made by subtracting the natural emissions. The inversion that simulates the closest anthropogenic estimate to the UNFCCC NGHGI (2019) is FLExKF-TM5-4DVAR_EMPA. In 2017, it reports 19.4 Tg CH4 yr−1, while NGHGIs report 18.5 Tg CH4 yr−1. Regarding trends, only FLExKF-TM5-4DVAR_EMPA shows a linear decreasing trend of 2.1 % yr−1, compared to the NGHGI data trend of 1.3 % yr−1 over their overlap period of 2005–2017, while other inversions show no significant trend. From this attempt we clearly note that few of the inversions showed the clear decline of NGHGIs. As NGHGI emissions are dominated by anthropogenic fluxes and decline by almost 30 % compared to 1990, this should also be seen in the corrected anthropogenic inversions. Therefore, we need to further investigate how well the NGHGIs reflect reality or how well the TD estimates capture the trends.

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f05

Figure 5(a) Total CH4 emissions from TD global ensembles based on surface stations (SURF) (yellow) and satellite concentration observations (GOSAT) (green) from 22 global models compared with UNFCCC NGHGI (2019) data (including LULUCF); (b) anthropogenic CH4 emissions from top-down global inversions based on surface stations (SURF) (yellow) and on satellite concentration observations (GOSAT) (green) from different estimates. Anthropogenic emissions from these inversions were obtained by removing the sum of the natural emissions (peatland, inland waters and geological fluxes shown in Fig. 4a) from the total estimates. For consistency with the global data we plot the global wetland emissions from the GCP inversions (blue). The UNFCCC NGHGI (2018) MS-reported uncertainty computed with the error propagation method (95 % confidence interval) is 9.29 % and represents the UNFCCC NGHGI (2018) uncertainty for all sectors (including LULUCF). The time series mean was computed for the common period 2010–2016. Two out of 11 SURF products (GELCA-SURF_NIES and TOMCAT-SURF_UOL) were not available for 2016.

Download

Global inversion estimates

Figure 5 compares TD global estimates with NGHGI data and gives information on the wetland emissions from global wetland models (Saunois et al., 2020). We present TD estimates of the total emissions (Fig. 5a) as well as estimates of anthropogenic emissions (Fig. 5b). The global inversion models were split according to the type of observations used, with 11 of them using satellites (GOSAT) and 11 using surface stations (SURF). Wetlands emissions provided by 22 global TD inversions from the Global Methane Budget (Saunois et al., 2020) are post-processed with prior ratio estimates for wetland CH4 emissions (Appendix B, Table B4).

For the common period 2010–2016 for the EU27 + UK, the two ensembles of regional and global models give a total CH4 emission mean (Fig. 5a) of 22.6 Tg CH4 yr−1 (GOSAT) and 23.7 Tg CH4 yr−1 (SURF) compared to 19.0 ± 1.7 Tg CH4 yr−1 for NGHGIs (Fig. 5a). The mean of the natural wetland emissions from the global inversions is 1.3 Tg CH4 yr−1 and partly explains the positive difference between total emissions from inversions and NGHGI anthropogenic emissions.

In an attempt to quantify the European TD anthropogenic CH4 component, in Fig. 5b we subtract from the total TD CH4 emissions once again the peatland emissions from the regional JSBACH–HIMMELI model and those from geological and inland water sources. The reason for correcting both regional and global inversions with the European peatland emissions from the JSBACH–HIMMELI model lies in the fact that they are in the range of the global wetland emission estimates for Europe (Saunois et al., 2020). Their median for all years (1.43 Tg CH4 yr−1, averaged over 2005–2017) is close to the BU estimates of peatland emissions from the JSBACH–HIMMELI model (1.44 Tg CH4 yr−1, averaged over 2005–2017).

For the 2010–2016 common period, the two ensembles of regional and global models give an anthropogenic CH4 emission mean (Fig. 5b) of 17.4 Tg CH4 yr−1 (GOSAT) and 23.7 Tg CH4 yr−1 (SURF) compared to 19.0 ± 1.7 Tg CH4 yr−1 for NGHGIs (Fig. 5b). For the same period, total CH4 emissions (Fig. 5a) from the SURF and GOSAT ensemble decrease by 0.5 % and 4.6 %, respectively. For anthropogenic CH4 emissions (Fig. 5b), the SURF and GOSAT ensemble show a decrease of 1.1 % and 6.3 %, respectively, compared to 7.3 % for the NGHGIs from 2010 to 2016.

In 2017, the TD ensemble that simulates the closest anthropogenic estimate (Fig. 5b) to the UNFCCC NGHGI (2019) is GOSAT, with the median of GOSAT inversions (16.4 Tg CH4 yr−1) falling within the uncertainty range of UNFCCC (18.4 ± 1.7 Tg CH4 yr−1).

3.2 Comparing N2O emission estimates from different approaches

3.2.1 Estimates of European and regional total N2O fluxes

Similarly, as done for CH4 (Sect. 3.1.1 and Fig. 1), we present results of the total N2O fluxes from the EU27 + UK and five main regions in Europe. Figure 7 summarizes the total N2O fluxes from NGHGIs (including LULUCF) for both the base year 1990 and mean of 2011–2015 period.

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f06

Figure 6Five-year-average (2011–2015) N2O emission estimates for EU27 + UK and five European regions (northern, western central, southern and eastern non-EU). The eastern European region does not include European Russia, and the UNFCCC uncertainty for the Republic of Moldova was not available. Northern Europe includes Norway. Central Europe includes Switzerland. The data belong to UNFCCC NGHGI (2019) submissions (grey) and base year 1990 (black star), two BU inventories (GAINS and EDGAR v5.0), natural emissions (emissions from lakes, rivers and reservoirs), and two TD total estimates (one regional European inversion (FLEXINVERT_NILU) and an average from three global inverse models from GN2OB, Tian et al., 2020). The relative error on the UNFCCC value represents the UNFCCC NGHGI (2018) MS-reported uncertainty computed with the error propagation method (95 % confidence interval): 80.0 % for the EU27 + UK, 50.3 % for eastern Europe non-EU, 26.6 % for northern Europe, 91.6 % for southern Europe, 51.9 % for western Europe and 46.0 % for central Europe.

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f07

Figure 7The contribution of changes in N2O anthropogenic emissions in the five UNFCCC sectors to the overall change in decadal mean, as reported to UNFCCC NGHGI (2019). The emissions follow the atmospheric convention, where positive numbers represent an emission to the atmosphere. The three stacked columns represent the average N2O emissions from each sector during three periods (1990–1999, 2000–2009 and 2010–2017), and percentages represent the contribution of each sector to the total reduction percentages between periods.

Download

The total UNFCCC estimates include emissions from all sectors. We plot these and compare them with fluxes from global datasets, BU models and TD inversions. We note that for all five regions the N2O emissions decreased between 29 % (northern Europe) and 43 % (western Europe) and for the EU27 + UK 37 % with respect to the NGHGI 1990 value. It also shows that at the regional European level the emissions from BU (anthropogenic and natural) and TD estimates agree well with reported NGHGI data within the high uncertainty reported by UNFCCC ( 80 %) or observed in the TD model range (for the EU27 + UK global inversions show a min/max range of 25 %–32 %, while regional inversions show a variability range of 9 %–11 % compared to the mean 2011–2015 value). This TD uncertainty is represented here by the variability in the model ensembles and denotes the range of the extremes (min and max) of estimates within each model ensemble. There is significant uncertainty in northern Europe, where the TD estimates indicate either a source or a sink (Fig. 6). The current observation network is sparse, which currently limits the capability of inverse models to quantify GHG emissions at the country or regional scale.

For all other regions BU anthropogenic emissions agree well with NGHGIs given uncertainties, though we note consistently higher estimates from TD regional and global models estimates. The difference is too high to be attributed to the natural emission, which is related here to inland waters as the only source and which ranges for all five regions between 0.2–1.3 kt N2O yr−1. The blue bar representing the natural emissions has a lower value estimates (Maavara et al., 2019, and Lauerwald et al., 2019), while the maximum value was calculated according to Yao et al. (2020). The higher values in Yao et al. (2020) are primarily due to N2O emissions from small streams, which are not included in Maavara et al. (2019), while both studies agree fairly well regarding N2O emissions from larger rivers (Yao et al., 2020).

3.2.2 NGHGI sectoral emissions and decadal changes

According to the UNFCCC (2019) NGHGI estimates for 2017, the EU27 + UK emitted GHGs totaling 3.9 Gt CO2 eq. (including LULUCF); of this total, N2O emissions accounted for  6 % (0.2 Gt CO2 eq. or 0.8 Mt N2O yr−1) (Fig. 7). France, UK and Germany contributed together 41 % of the total N2O emissions, respectively, which slightly higher than for CH4 (Appendix B2, Fig. B1b).

The data in Fig. 7 show anthropogenic CH4 emissions and their change from one decade to the next, from UNFCCC NGHGI (2019), with the contribution from different UNFCCC sectors. In 2017, NGHGIs reported anthropogenic emissions from the EU27 + UK for the four UNFCCC sectors (excluding LULUCF) (Table 1) to be 0.8 Tg N2O yr−1. The agricultural N2O emissions accounted for 76 % (± 107 %) of the total EU27 + UK emissions followed by emissions from the energy sector with 12 % (± 23 %). We exclude fire emissions as they only account for 1.8 % of the total N2O emissions in the EU27 + UK.

Between the 1990s and the 2000s, the net 17.7 % reduction originates largely from the IPPU and agriculture sectors, which contributed 13.5 % and 4.2 %, respectively. For the period between the 2000s and 2010–2017, the net 15.2 % reduction was again mainly attributed to the IPPU sector (14.1 %), despite very small increases from the LULUCF and waste sectors (+0.6 %).

We note that in 2017 the amount of emissions from the IPPU sector had already decreased by 98 % compared to 1990 and was only 3.5 kt N2O yr−1. Although the IPPU sector contributes in 2017 only 4 % to the total N2O emissions, it was the sector with the largest reduction. IPPU sector emissions are mainly linked to the production of nitric acid (e.g., used in fertilizer production) and adipic acid (e.g., used in nylon production). In the late 1990s and early 2000s the five European adipic acid plants were equipped with efficient abatement technology, cutting emissions by 95 %–99 %, largely through voluntary agreements of the companies. Much of the remaining IPPU emissions, from nitric acid plants, were cut in a similar manner around 2010, a development that has been connected with the introduction of the European Emission Trading System that made it economically interesting for companies to apply emission abatement technologies (catalytic reduction of N2O in the flue gas) to reduce their emissions.

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f08

Figure 8(a) Total EU27 + UK anthropogenic N2O emissions from UNFCCC NGHGI (2019) submissions and MS-NRT 2018 compared to global BU inventories for agriculture (CAPRI, FAOSTAT, DayCent) and all sectors excluding LULUCF (EDGAR v5.0, GAINS). Sectors: (b) energy, (c) IPPU, (d) agriculture and (e) waste. LULUCF not included. CAPRI reports one value for Belgium and Luxembourg. The UNFCCC NGHGI (2018) MS-reported uncertainty was computed with the error propagation method (95 % confidence interval) and is 86 % for the total EU27 + UK (excluding LULUCF), 23 % for energy, 16 % for IPPU, 107 % for agriculture and 626 % (in the figure 100 %) for waste. Uncertainty for EDGAR v5.0 was calculated for 2015, and the min/max values for all sectors are as follows: EU27 + UK total 48/126, energy 12/250, IPPU 13/19, agriculture 74/191 and waste 63/166; it represents the 95 % confidence interval of a lognormal distribution. The mean values on the right-hand side reflect the values for the common overlapping period 1990–2015. The last reported year in this study refers to 2017 (UNFCCC and FAOSTAT), 2015 (EDGAR v5.0 and GAINS every 5 years), 2013 (CAPRI), 2018 ECOSSE and 2011–2015 DayCent.

Download

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f09

Figure 9Total N2O emissions from UNFCCC NGHGI (2019) (including LULUCF) and MS-NRT 2018, compared to the FLEXINVERT_NILU regional inversion over Europe and GCP global inversions (TOMCAT, CAMS-N2O and MIROC4-ACTM). In blue are the natural inland water (lakes_rivers_reservoirs_ULB) N2O emissions. The UNFCCC NGHGI (2018) reported uncertainty computed with the error propagation method (95 % confidence interval) is 80 % and represents the UNFCCC NGHGI (2018) uncertainty for all sectors (including LULUCF). The uncertainty for EDGAR v5.0 was calculated for 2015 and ranged from a minimum of 37 % to a maximum of 73 %; it represents the 95 % confidence interval of a lognormal distribution. The last reported year in this study refers to 2014 (TOMCAT_LEEDS), 2016 (MIROC4-ACTM), 2017 (UNFCCC NGHGI, 2019; FLEXINVERT_NILU; and CAMS-N2O). For MS-NRT, the following countries are missing information from the LULUCF sector: Austria, Estonia, Croatia, Hungary, Latvia, Malta and Slovenia. The time series mean for TD products was computed for the overlapping period 2005–2014.

Download

3.2.3 NGHGI estimates compared with bottom-up inventories

Figure 8 compares the six bottom-up inventories with UNFCCC NGHGI (2019) data and shows that all of them fall on the NGHGI line (Fig. 8a), noting that GAINS only provides emissions every 5 years. Each inventory shows a very good agreement with each other and the NGHGI estimates until 2005. After 2005 the slight increased trend is influenced by the IPPU (Fig. 8c) and waste (Fig. 8e) sectors, with estimates of both EDGAR v5.0 and GAINS for total anthropogenic N2O emissions in the year 2015 being 15.6 % higher than UNFCCC NGHGI (2019). The jump seen in Fig. 8e for EDGAR v5.0 waste emissions in year 2000 is due to updates on activity data for domestic waste water treatment using FAO statistics for 2000–2016, while the previous period 1990–1999 remained unchanged using the data from the previous version of EDGAR v4.3.2. For agriculture (Fig. 8d) five models/inventories show a very good match with the NGHGI. Over 1990–2015, we found linear trends of 0.7 % yr−1 in NGHGIs, GAINS and EDGAR v5.0. This provides further evidence that the sources rely on the same basic activity data from FAOSTAT and follow the IPCC EF Tier 1 or 2 approach (Petrescu et al., 2020a). In contrast, ECOSSE estimates do not use the FAO fertilizer application rate database but instead calculate ideal fertilizer application rates from the nitrogen demand of the crops. This means that it can severely underestimate the applied fertilizer amounts for some areas (e.g., the Netherlands, Denmark or northwestern Germany), and the results are more indicative of emissions under idealized fertilizer application rates. Additionally, as mentioned above, the model simulates only the direct emissions.

In the NGHGI (2018) submissions, the EU27 + UK Tier 1 total uncertainty (based on the IPCC chap. 3 error propagation method described in detail by Petrescu et al., 2020a) for the waste sector was 626 %. The sectoral activity responsible for this high uncertainty was the wastewater treatment and discharge (913 %), and this remains one of the most uncertain sources of N2O having the highest emissions in the waste sector. Emissions are known to vary markedly in space and time even within a single wastewater treatment plant (Gruber et al., 2020), a fact that only recently has been properly accounted for in the inventory guidelines (IPCC, 2019a). However, the total emissions from the waste sector account for only 4.4 % of the total EU27 + UK N2O emissions (excluding LULUCF).

3.2.4 NGHGI estimates compared to atmospheric inversions

Figure 9 compares inversion estimates of the total regional (FLEXINVERT_NILU) and global (three models) N2O inversions with UNFCCC NGHGI (2019). The minimum–maximum range of all inversions is within the 2σ uncertainty of NGHGIs, with the median of global inversions being on average 42 % or 0.4 Tg N2O yr−1 higher than NGHGIs. Over the period 2005–2017, the regional FLEXINVERT_NILU is 65 % higher than UNFCCC NGHGI (2019). From the three global inversions, two show consistently higher estimates (MIROC4-ACTM and TOMCAT) as well as high variability. Regarding trends, only FLEXINVERT_NILU shows a decreasing trend of 2.1 % yr−1 over 2005–2017, compared to UNFCCC NGHGI (2019) with a decreasing trend of 1.2 % yr−1. The global CAMS-N2O inversion agrees best in its mean value (1.0 Tg N2O yr−1) with the UNFCCC estimate (0.9 Tg N2O yr−1) but not in its trend. The higher emissions from TD estimates may be at least in part due to the fact that they include natural emissions of N2O, which are not considered in NGHGI reporting. One estimate (from the O−CN land ecosystem model) is that the natural emissions could amount to 11 % of those reported in NGHGIs for the EU27 + UK region. In addition, the EFs used in NGHGI reporting are very uncertain (up to 300 % for direct agricultural emissions) so there may be a systematic error in these.

For the N2O we do not present the corrected anthropogenic value because the only natural flux, from inland waters, is very low (2.7 kt N2O yr−1), and when subtracted from the four inversions the change is almost negligible. Part of the inland water natural estimate is considered anthropogenic in Europe and is due to the leaching of N fertilizers from agriculture. In this study, it accounts for 66 % of the total inland waters emissions.

4 Data availability

All data files reported in this work which were used for calculations, and figures are available for public download at https://doi.org/10.5281/zenodo.4590875 (Petrescu et al., 2020b). The data we submitted are reachable with one click (without the need for entering login and password) and downloadable with a second click, consistent with the two-click access principle for data published in ESSD (Carlson and Oda, 2018). The data and the DOI number are subject to future updates and only refer to this version of the paper. The raw gridded data, according to the VERIFY consortium governing document, will be made publicly available 12 months after their publication in ESSD. However, they could be released earlier, upon request to the specific author.

Please also see Table 2 for an overview of data sources for CH4 and N2O emissions used in this study.

5 Summary and concluding remarks

This study represents the first comprehensive European verification that compares total and sectoral European CH4 and N2O emission estimates from BU (anthropogenic and natural) with TD estimates in order to assess their use for verification purposes with the UNFCCC NGHGI reporting. Above, in the results sections, we discussed differences between estimates. Identification of source-specific uncertainty is key in understanding these differences and will lead to the reduction of the overall uncertainty in GHG inventories. More specifically, we present the first EU27 + UK and European regional 2011–2015 averaged results for CH4 (Fig. 1) and N2O (Fig. 6) compared to the NGHGI emissions (including LULUCF) for the same period and the one reported in 1990, in the framework of the future global stocktake estimate.

Regarding sources of inconsistencies between CH4 BU estimates and NGHGI data (Fig. 3), at the EU27 + UK level they are mainly caused by the use of different methodologies in calculating emissions as highlighted in Petrescu et al. (2020a). Both BU inventories and the NGHGIs use similar activity data, and the default EFs reported in the 2006 IPCC Guidelines make all data sources agree rather well; thus, inventory spread may not be indicative of the uncertainty. For global consistency purposes, EDGAR v5.0 mostly uses Tier 1 approaches in calculating emissions and uncertainties, a fact which triggers differences with other data sources (GAINS for all sectors and CAPRI for agriculture). Within the UNFCCC reporting process, the two most important emission sectors after agriculture are energy and waste and contribute to the highest reduction percentages (see Fig. 1). Another reason for small inconsistencies between datasets is the allocation of emissions to different sectors, with different data sources using different versions of IPCC Guidelines (e.g., EDGAR v5.0).

For N2O anthropogenic emissions, all BU data sources show good agreement with the UNFCCC NGHGI (2019) data in both trends and means (Fig. 8), and agriculture remains the largest emitter (e.g., urea fertilizers). The uncertainties reported by NGHGIs are very large and will need further improvement.

Regarding the TD estimates, our exercise shows that the comparison between CH4 inversion estimates and NGHGIs is highly uncertain because of the large spread in the inversion results. As TD inversions do not fully distinguish between all emission sectors used by NGHGIs and report either total emissions or a coarse sectorial partitioning, their comparison to NGHGIs is only possible for total emissions. It is also necessary to make an adjustment for natural emissions, which are included in TD inversions but not reported by the NGHGIs. However, the natural N2O emissions do not explain the difference between BU and TD (452 kt N2O), and more research is needed to identify the source of discrepancies.

Some studies (Fronzek et al., 2018) show that model ensembles work well in simulating highly uncertain variables. In general regional inversions show less spread than the global inversions as they used recent updates of transport models and higher-resolution transport. Total CH4 from regional inversions shows a minimum–maximum range of 8.7 Tg CH4 yr−1 compared to 12.4 Tg CH4 yr−1 from global GOSAT inversions and 13.5 Tg CH4 yr−1 from global SURF inversions (Figs. 4 and 5). The global models use fewer observations for Europe compared to the European regional inversions and thus are expected to have larger uncertainties for the European fluxes. In addition, the global models are at coarser resolution and thus likely have larger model representation errors compared to the regional ones, which may contribute to further systematic uncertainty for the European fluxes.

A key challenge for the inversion CH4 community remains the separation of emissions in specific source sectors, as derived total emissions may also include natural emissions (or removals), while in the case of N2O this will not be possible due to the use of definitions (e.g., “natural” N2O emissions are defined as the level of emission in the pre-industrial period). It is therefore not possible to separate between N2O from natural or anthropogenic sources because natural (or unmanaged) soils can have both natural and anthropogenic emissions, while anthropogenic (managed) agricultural soils can also have a level of natural emissions according to the definition of natural. Therefore, the goal of TD inversions for estimating N2O emissions should mainly focus on the trends. Furthermore, the accuracy of derived emissions and the spatial scales at which emissions can be estimated depend on the quality and density of measurements and the quality of the atmospheric models (Bergamaschi et al., 2018b). Significant further developments of the global observation system and the top-down methods would be required to support the implementation of the Paris Agreement.

We provide for the EU27 + UK a consolidated synthesis of the relatively uncertain CH4 and N2O emissions making use of consistently derived BU and TD estimates over the region of Europe, which might illustrate the importance of regional consistent analyses that form the basis of the multilateral facilitative consideration of progress under the enhanced transparency framework of the Paris Agreement. However, the implementation of the Paris Agreement requires accurate quantification of GHG emissions in order to track the progress of all parties with their nationally determined contributions and to assess collective progress towards achieving the purpose of this agreement and its long-term goals (stocktake). As this will be mainly achieved and build upon BU methodologies developed by the IPCC, we need to take into consideration the potential to quantify GHG emissions by using top-down methods (inverse modeling) (Bergamaschi et al., 2018b). One advantage of the inverse estimate is that it provides total emission estimates. Therefore, the capability to quantify anthropogenic emissions depends on the magnitude of natural sources and sinks and the capability to quantify them.

As stated in the introduction, our aim was to identify in this synthesis the issues which cause the differences between NGHGIs, BU and TD to further improve and build a pathway to a verification system (BU use of activity data, emission factor and emission allocation (CH4), very large NGHGI reported uncertainties which need to be reassessed (N2O) and higher TD estimates compared to inventories (CH4 and N2O). Additionally, we advocate the need of analyzing the seasonality of emissions, which is of great importance for CH4 (wetland emission estimates have large uncertainties and show large variability in the spatial (seasonal) distribution) and N2O (agriculture fertilizer application), included mostly in the TD prior information but not included in the reported IPCC Guidelines. In the climate mitigation process these seasonal variations may play an important role for a better quantification of sector-specific uncertainties.

Appendix A: Data sources, methodology and uncertainty descriptions

The country-specific plots are found at http://webportals.ipsl.jussieu.fr/VERIFY/FactSheets/ (upon registration, last access: February 2021) (v1.24).

VERIFY project

VERIFY's primary aim is to develop scientifically robust methods to assess the accuracy and potential biases in national inventories reported by the parties through an independent pre-operational framework. The main concept is to provide observation-based estimates of anthropogenic and natural GHG emissions and sinks as well as associated uncertainties. The proposed approach is based on the integration of atmospheric measurements, improved emission inventories, ecosystem data, and satellite observations, as well as on an understanding of processes controlling GHG fluxes (ecosystem models, GHG emission models).

Two complementary approaches relying on observational data streams will be combined in VERIFY to quantify GHG fluxes:

  1. atmospheric GHG concentrations from satellites and ground-based networks (top-down atmospheric inversion models) and

  2. bottom-up activity data (e.g., fuel use and emission factors) and ecosystem measurements (bottom-up models).

For CO2, a specific effort will be made to separate fossil fuel emissions from ecosystems fluxes. For CH4 and N2O, we will separate agricultural from fossil fuel and industrial emissions. Finally, trends in the budget of the three GHGs will be analyzed in the context of NDC targets.

The objectives of VERIFY are as follows.

  • Objective 1. Integrate the efforts between the research community, national inventory compilers, operational centers in Europe, and international organizations towards the definition of future international standards for the verification of GHG emissions and sinks based on independent observation.

  • Objective 2. Enhance the current observation and modeling ability to accurately and transparently quantify the sinks and sources of GHGs in the land use sector for the tracking of land-based mitigation activities.

  • Objective 3. Develop new research approaches to monitor anthropogenic GHG emissions in support of the EU commitment to reduce its GHG emissions by 40 % by 2030 compared to the year 1990.

  • Objective 4. Produce periodic scientific syntheses of the observation-based GHG balance of EU countries and practical policy-oriented assessments of GHG emission trends and apply these methodologies to other countries.

For more information on project team and products/results check https://verify.lsce.ipsl.fr/ (last access: February 2021).

Table A1Country grouping used for comparison purposes between BU and TD emissions. Countries highlighted in italic are not discussed in the current 2019 synthesis mostly because of unavailability of UNFCCC NGHGI reports (non-Annex I countries*) but are present on the following web portal: http://webportals.ipsl.jussieu.fr/VERIFY/FactSheets/ (last access: February 2021). Results of Annex I countries (NOR, CHE, ISL) and non-EU eastern European countries (EAE) are represented in Fig. 4.

* Non-Annex I countries are mostly developing countries. The reporting to UNFCCC is implemented through national communications (NCs) and biennial update reports (BURs): https://unfccc.int/national-reports-from-non-annex-i-parties (last access: November 2020).

Download XLSX

Table A2Main methodological changes of the current study with respect to Petrescu et al. (2020a). The sectors are highlighted in italics; the changes with respect to AFOLU (Petrescu et al., 2020a) are shown in bold. “NA” means that there are no data available.

Download Print Version | Download XLSX

A1 Anthropogenic CH4 emissions (sectors energy, IPPU, agriculture, LULUCF and waste)

A1.1 Bottom-up CH4 emission estimates

UNFCCC NGHGI (2019)

Under the UNFCCC convention and its Kyoto Protocol national greenhouse gas (GHG) inventories is the most important source of information to track progress and assess climate protection measures by countries. In order to build mutual trust in the reliability of the GHG emission information provided, national GHG inventories are subject to standardized reporting requirements, which have been continuously developed by the Conference of the Parties (COP)9. The calculation methods for the estimation of greenhouse gases in the respective sectors are determined by the methods provided by the 2006 IPCC Guidelines for National Greenhouse Gas Inventories (IPCC, 2006). These guidelines provide detailed methodological descriptions to estimate emissions and removals as well as provide recommendations to collect the activity data needed. As a general overall requirement, the UNFCCC reporting guidelines stipulate that reporting under the convention and the Kyoto Protocol must follow the five key principles of transparency, accuracy, completeness, consistency and comparability (TACCC). The reporting under UNFCCC shall meet the TACCC principles. The three main GHGs are reported in time series from 1990 up to 2 years before the due date of the reporting. The reporting is strictly source category based and is done under the common reporting format tables (CRF), downloadable from the UNFCCC official submission portal: https://unfccc.int/process-and-meetings/transparency-and-reporting/reporting-and-review-under-the-convention/greenhouse-gas-inventories-annex-i-parties/national-inventory-submissions-2019 (last access: January 2021)

The UNFCCC NGHGI anthropogenic CH4 emissions include estimates from four key sectors for the EU27 + UK: (1) energy, (2) industrial processes and product use (IPPU), (3) agriculture, and (5) waste. The tier method a country applies depends on the national circumstances and the individual conditions of the land, which explain the variability of uncertainties among the sector itself as well as among EU countries. The LULUCF CH4 emissions are very small but are included in some figures (see Table 1).

The uncertainty methodology for the NGHGI UNFCCC submissions is based on chap. 3 of the 2006 IPCC Guidelines for National Greenhouse Gas Inventories and is explained in Appendix B of Petrescu et al. (2020a).

EDGAR v5.0

The Emissions Database for Global Atmospheric Research (EDGAR) is an independent global emission inventory of greenhouse gases (GHG) and air pollutants developed by the Joint Research Centre of the European Commission (https://edgar.jrc.ec.europa.eu/index.php, last access: January 2021). The non-CO2 component in EDGAR v5.0 covers a long time series of emissions starting in 1970 until the t−4. CH4 emissions are estimated for all anthropogenic emission sectors with the exception of land use, land use change and forestry (LULUCF) at the country and annual level in a consistent and comparable way for all world countries. Emissions are computed using activity data from international statistics (e.g., IEA (2017), FAO (2017), USGS (2019), etc.); emission factors from IPCC Guidelines (IPCC, 2006); and scientific literature, technology and abatement measures incorporation. Once the emission database is compiled for all countries, sectors and pollutants, annual emission data are disaggregated to monthly emissions applying sector- and country-specific yearly emission profiles (Crippa et al., 2019). In addition, monthly emissions are spatially distributed over global grid maps with a resolution of 0.1× 0.1, making use of sector-specific spatial proxies (Janssens-Maenhout et al., 2019).

The latest version of the EDGAR database EDGAR v5.0 contains estimated CH4 emissions from 1970 until 2015 based fully on statistical data (https://edgar.jrc.ec.europa.eu/overview.php?v=50_GHG, last access: January 2021). The EDGAR v5.0 updated waste emissions were quantified using the first-order-decay method, combining nationally defined inputs (for waste generation rates and compositions) and IPCC's regional default values for parameters associated with waste degradation processes (specific mass of biodegradable organic carbon, the methane volumetric fraction in the obtained landfill gas and the half-life time for each waste component). The total landfilled waste was split into six streams: food and organic waste type, paper and cardboard, textiles, rubber, wood, and sludge and similar effluents. This was done by making use of the EUROSTAT waste database (EUROSTAT, 2020) but also by employing data from waste composition for municipal/household-type waste from Silpa (Silpa et al., 2018).

Uncertainties. The methodological description is explained in detail in Appendix B Petrescu et al. (2020a) and Solazzo et al. (2021).

CAPRI

CAPRI is an economic, partial equilibrium model for the agricultural sector focused on the EU (Britz and Witzke, 201410; Weiss and Leip, 201211. CAPRI stands for Common Agricultural Policy Regionalised Impact analysis, and the name hints at the main objective of the system: assessing the effect of CAP policy instruments not only at the EU or member state level but also at the sub-national level. The model is calibrated for the base year (currently 2012) and then baseline projections are built, allowing the ex ante evaluation of agricultural policies and trade policies on production, income, markets, trade and the environment.

Among other environmental indicators, CAPRI simulates CH4 emissions from agricultural production activities (enteric fermentation, manure management, rice cultivation, agricultural soils). Activity data are mainly based on FAOSTAT and EUROSTAT statistics, and estimation of emissions follows IPCC 2006 methodologies, with a higher or lower level of detail depending on the importance of the emission source. Details on the CAPRI methodology for emissions calculations are referenced in the Annex Table A1.

Uncertainties are not available for the CAPRI estimates.

GAINS

Specific sectors and abatement technologies in GAINS vary by the specific emitted compound, with source sector definition and emission factors largely following the IPCC methodology at the Tier 1 or Tier 2 level. GAINS includes in general all anthropogenic emissions to air but does not cover emissions from forest fires, savannah burning and land use/land use change. Emissions are estimated for 174 countries/regions, with the possibility to aggregate to a global emission estimate and spanning a timeframe from 1990 to 2050 in 5-year intervals. Activity drivers for macroeconomic development, energy supply and demand, and agricultural activities are entered externally; GAINS extends with knowledge required to estimate “default” emissions (emissions occurring due to an economic activity without emission abatement) and emissions and costs of situations under emission control (see Amann et al., 2011).

The GAINS model covers all source sectors of anthropogenic methane (CH4) emissions: agricultural sector emissions from livestock, rice cultivation and agricultural waste burning; energy sector emissions from upstream and downstream sources in fossil fuel extraction and use; and emissions from handling and treatment of solid waste and wastewater source sectors. A description of the modeling of CH4 emissions in GAINS is presented in Höglund-Isaksson et al. (2020). Generation of solid waste and the carbon content of wastewater are derived within the model in consistency with the relevant macroeconomic scenario. The starting point for estimations of anthropogenic CH4 is the methodology recommended in the IPCC (2006 and revision in 2019) Guidelines, for most source sectors using country-specific information to allow for deriving country- and sector/technology-specific emission factors at a Tier 2 level. Consistent methodologies were further developed to estimate emissions from oil and gas systems (Höglund-Isaksson, 2017) and solid waste (Höglund-Isaksson et al., 2018; Gómez-Sanabria et al., 2018). Emission factors are specified in a consistent manner across countries for given sets of technology and with past implementation of emission abatement measures reflected as changes in technology structures. The resulting emission estimates are well comparable across geographic and temporal scales. The GAINS approach to calculate waste emissions is developed in consistency with the first-order-decay method recommended by IPCC (2006 and 2019 revision), applying different decay periods when estimating emissions from flows of different types of organic waste, i.e., food and garden, paper, wood, textile and other. Data on waste generation, composition and treatment are taken from EUROSTAT (https://ec.europa.eu/eurostat/statistics-explained/index.php/Waste_statistics, last access: 2019) and complemented with national information from the UNFCCC (2019) common reporting format tables on the amounts of waste diverted to landfills of various management levels and to treatment, e.g., recycling, composting, biodigestion and incineration.

Uncertainties. Uncertainty is prevalent among many different dimensions in the estimations of emissions, abatement potentials and costs. When constructing global bottom-up emission inventories at a detailed country and source level, it is inevitable that some information gaps will be bridged using default assumptions. As it is difficult to speculate about how such sources of uncertainty affect resulting historical and future emission estimates, we instead address uncertainty in historical emissions by making comparisons to estimates by other publicly available and independently developed bottom-up inventories and various top-down estimates consistent with atmospheric measurements and inverse model results. Although existing publicly available global bottom-up inventories adhere to the recommended guidelines of the IPCC (2006), the flexibility in these is large, and results will depend on the availability and quality of gathered source information. There is accordingly a wide range of possible sources of uncertainty built into estimations in such comprehensive efforts. Having a pool of independently developed inventories, each with its own strengths and weaknesses, can improve the understanding of the scope for uncertainty, in particular when compared against top-down atmospheric measurements.

FAOSTAT

The Food and Agriculture Organization of the United Nations (FAO) provides CH4 emissions from agriculture available at http://www.fao.org/faostat/en/#data/GT/visualize (last access: January 2021) and emissions from land use, available at http://www.fao.org/faostat/en/#data/GL/visualize (last access: January 2021). The FAOSTAT emissions database is computed following Tier 1 2006 IPCC Guidelines for National GHG Inventories (http://www.ipcc-nggip.iges.or.jp/public/2006gl/index.html, last access: December 2019). Country reports to FAO on crops, livestock and agriculture use of fertilizers are the source of activity data. Geospatial data are the source of activity data (AD) for the estimates from cultivation of organic soils, biomass and peat fires. GHG emissions are provided by country, regions and special groups, with global coverage, relative to the period 1961–present (with annual updates, currently to 2017) and with projections for 2030 and 2050, expressed as CO2 eq. for CH4, by underlying agricultural emission sub-domain and by aggregate (agriculture total, agriculture total plus energy, agricultural soils). LULUCF emissions consist of CH4 (methane) associated with biomass and peat fires.

Uncertainties were computed by Tubiello et al. (2013) but are not available in the FAOSTAT database.

A1.2 Top-down CH4 emission estimates

FLEXPART(FLExKF-TM5-4DVAR)

FLExKF applies an extended Kalman filter (Brunner et al., 2012) in combination with backward Lagrangian transport simulations using the model FLEXPART (Stohl et al., 2005; Pisso et al., 2019). It optimizes surface–atmosphere fluxes by assimilating atmospheric observations in a sequential manner, which allows for an analytical solution for relatively large inversion problems (long time periods, number stations O(100)). Since model–observation residuals typically follow a lognormal distribution, the method optimizes log-transformed emissions, which also guarantees a positive solution. Source-receptor matrices (Seibert and Frank, 2004) were computed at 0.25× 0.25 resolution with FLEXPART driven by ECMWF ERA-Interim meteorological fields in the same way as for FLEXINVERT_NILU. Backward simulations were limited to 10 d prior to each observation and to the domain 15 W–35 E, 30–75 N. Fluxes were estimated for this domain on a monthly basis at 0.5× 0.5 resolution. For the version used in this study, FLEXPART(FLExKF-TM5-4DVAR)_EMPA, the background mole fraction was taken from a global TM5-4DVAR assimilation run (Bergamaschi et al., 2018a) where the above domain was cut out following the two-step approach of Rödenbeck et al. (2009).

Uncertainties. The uncertainty in the posterior fluxes is composed of random and systematic errors. The random uncertainties are represented by the posterior error covariance matrix provided by the Kalman filter, which combines errors in the prior fluxes with errors in the observations and model representation. Systematic uncertainties primarily arise from systematic errors in modeled atmospheric transport and in background mole fractions but also include aggregation errors, i.e., errors arising from the way the flux variables are discretized in space and time.

FLEXINVERT_NILU

The FLEXINVERT_NILU framework is based on Bayesian statistics and optimizes surface–atmosphere fluxes using the maximum probability solution (Rodgers, 2000). Atmospheric transport is modeled using the Lagrangian model FLEXPART (Stohl et al., 2005; Pisso et al., 2019) run in the backwards time mode to generate a so-called source-receptor matrix (SRM). The SRM describes the relationship between the change in mole fraction and the fluxes discretized in space and time (Seibert and Frank, 2004) and was calculated for 8 d prior to each observation. For use in the inversions, FLEXPART was driven using ECMWF operational analysis wind fields. The state vector consisted of prior fluxes discretized on an irregular grid based on the SRMs (Thompson and Stohl, 2014). This grid has finer resolution (in this case the finest was 0.25× 0.25) where the fluxes have a strong influence on the observations and coarser resolution where the influence is only weak (the coarsest was 2× 2). The fluxes were solved at 10 d temporal resolution. The state vector also included scalars for the background contribution. The background mixing ratio, i.e., the contribution to the mixing ratio that is not accounted for in the 8 d SRMs, was estimated by coupling the termination points of backwards trajectories (modeled using virtual particles) to initial fields of methane simulated with the Lagrangian FLEXPART-CTM model, which was developed at Empa based on FLEXPART (Stohl et al., 2005; Pisso et al., 2019). In these simulations, we applied the data assimilation method described by Groot Zwaaftink et al. (2018) that constrains modeled fields with surface observations through nudging.

Uncertainties. The posterior fluxes are subject to systematic errors primarily from (1) errors in the modeled atmospheric transport; (2) aggregation errors, i.e., errors arising from the way the flux variables are discretized in space and time; (3) errors in the background methane fields; and (4) the incomplete information from the observations and hence the dependence on the prior fluxes. In addition, there is, to a smaller extent, some error due to calibration offsets between observing instruments. Uncertainties in the observation space were inflated to take into account the model representation errors

InGOS and TM5-4DVAR

The atmospheric models used within the European FP7 project InGOS (Integrated non-CO2 Greenhouse gas Observing System) are described by Bergamaschi et al. (2018a) and their Supplement (https://doi.org/10.5194/acp-18-901-2018-supplement). The models include global Eulerian models with a close-up over Europe (TM5-4DVAR, TM5-CTE, LMDZ), regional Eulerian models (CHIMERE) and Lagrangian dispersion models (STILT, NAME, COMET). The horizontal resolutions over Europe are  1.0–1.2 (longitude) × 0.8–1.0 (latitude) for the global models (close-up) and  0.17–0.56 (longitude) × 0.17–0.5 (latitude) for the regional models. Most models are driven by meteorological fields from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis (Dee et al., 2011). In the case of STILT, the operational ECMWF analyses were used, while for NAME meteorological analyses of the Met Office United Model (UM) were employed. The regional models use boundary conditions (background CH4 mole fractions) from inversions of the global models (STILT from TM3, COMET from TM5-4DVAR, CHIMERE from LMDZ) or estimate the boundary conditions in the inversions (NAME) using baseline observations at Mace Head as prior estimates. In the case of NAME and CHIMERE, the boundary conditions are further optimized in the inversion. The inverse modeling systems applied in this study use different inversion techniques. TM5-4DVAR, LMDZ, and TM3-STILT use 4DVAR variational techniques, which allow optimization of emissions of individual grid cells. These 4DVAR techniques employ an adjoint model in order to iteratively minimize the cost function using a quasi-Newton (Gilbert and Lemaréchal, 1989) or conjugate gradient (Rödenbeck, 2005) algorithm. The NAME model applies a simulated annealing technique, a probabilistic technique for approximating the global minimum of the cost function. In CHIMERE and COMET, the inversions are performed analytically after reducing the number of parameters to be optimized by aggregating individual grid cells before the inversion. TM5-CTE applies an ensemble Kalman filter (EnKF) (Evensen, 2003) with a fixed-lag smoother (Peters et al., 2005).

Uncertainty. In general, the estimated model uncertainties depend on the type of station and for some models (TM5-4DVAR and NAME) also on the specific synoptic situation. In InGOS the uncertainty of the ensemble was calculated as a 1σ estimate. Bergamaschi et al. (2015) showed that the range of the derived total CH4 emissions from north-western and eastern Europe using four different inverse modeling systems was considerably larger than the uncertainty estimates of the individual models because the latter typically use Bayes' theorem to calculate the reduction of assumed a prior emission uncertainties by assimilating measurements (propagating estimated observation and model errors to the estimated emissions). An ensemble of inverse models may provide more realistic overall uncertainty estimates, since estimates of model errors are often based on strongly simplified assumptions and do not represent the total uncertainty.

InTEM-NAME

The Inversion Technique for Emission Modelling (InTEM) (Arnold et al., 2018) uses the NAME (Numerical Atmospheric dispersion Modeling Environment) (Jones et al., 2007) atmospheric Lagrangian transport model. NAME is driven by analysis 3-D meteorology from the UK Met Office Unified Model (Cullen, 1993). The horizontal and vertical resolution of the meteorology has improved over the modeled period from 40 to 12 km (1.5 km over the UK). InTEM is a Bayesian system that minimizes the mismatch between the model and the atmospheric observations given the constraints imposed by the observation and model uncertainties and prior information with its associated uncertainties. The direction- (latitude and longitude) and altitude-varying background concentration and observation station bias are solved for within the inverse system along with the spatial distribution and magnitude of the emissions. The time-varying prior background concentration for the DECC network stations is derived from the MHD observations when they are very largely sensitive only to northern Canada (Arnold et al., 2018). The prior bias (that can be positive or negative) for each station is set to zero with an uncertainty of 1 ppb. The observations from each station are assumed to have an exponentially decreasing 12 h time correlation coefficient and, between stations, a 200 km spatial correlation coefficient. The observations are averaged into 2 h periods. The uncertainty of the observations is derived from the variability of the observations within each 2 h period. The modeling uncertainty for each 2 h period at each station varies and is defined as the larger of the median pollution events in that year at that station or 16.5 % of the magnitude of the pollution event. These values have been derived from analysis of the observations of methane at multiple heights at each station across the DECC network. Each inversion is repeated 24 times; each time 10 % of the observations per year per station are randomly removed in 5 d intervals and the results and uncertainty averaged.

Uncertainty. This random removal of observations allows a greater exploration of the uncertainty, given the potential for some of the emission sources to be intermittent within the time period of the inversion.

CTE-CH4 Europe, CTE-SURF and CTE-GOSAT

CarbonTracker Europe CH4 (CTE-CH4) (Tsuruta et al., 2017) applies an ensemble Kalman filter (Peters et al., 2005) in combination with the Eulerian transport model TM5 (Krol et al., 2005). It optimizes surface fluxes weekly and assimilates atmospheric CH4 observations. TM5 was run at 1× 1 resolution over Europe and 6× 4e resolution globally, constrained by 3-hourly ECMWF ERA-Interim meteorological data. The photochemical sink of CH4 due to tropospheric and stratospheric OH as well as stratospheric Cl and O(1D) were pre-calculated based on Houweling et al. (2014) and Brühl and Crutzen (1993) and not adjusted in the optimization scheme.

Three experiments were conducted, which differ in (1) sets of prior fluxes, (2) sets of assimilated observations and (3) optimization resolution over the Northern Hemisphere. CTE-FMI uses sets of prior fluxes from LPX-Bern DYPTOP (Stocker et al., 2014) for biospheric, EDGAR v4.2 FT2010 (Janssens-Maenhout et al., 2013) for anthropogenic, GFED v4 (Giglio et al., 2013) for biomass burning, Ito and Inatomi (2012) for termites, and Tsuruta et al. (2017) for ocean sources. CTE-SURF and CTE-GOSAT use sets of prior fluxes from the Global Carbon Project (Saunois et al., 2020). CTE-FMI and CTE-SURF assimilated ground-based surface CH4 observations, while CTE-GOSAT assimilated GOSAT XCH4 retrievals from NIES v2.72. CTE-FMI optimized fluxes at 1× 1 resolution over northern Europe, northeast Russia and southeast Canada; 6× 4 resolution over other parts of the Northern Hemisphere land; and region-wise (combined TransCom regions and soil type) over the Southern Hemisphere and ocean. CTE-SURF and CTE-GOSAT fluxes were optimized at 1× 1 resolution over Europe and region-wise elsewhere globally.

Uncertainty. The prior uncertainty is assumed to be a Gaussian probability distribution function, where the error covariance matrix includes errors in prior fluxes, observations and transport model representations. The uncertainty for the prior fluxes was assumed to be 80 % of the fluxes over land and 20 % over ocean, with the correlation between grid cells or regions being 100–500 km over land and 900 km over ocean. The uncertainties for observations and transport model representations vary between observations, with the minimum aggregated uncertainty being 7.5 ppb for surface observations and 15 ppb for GOSAT data. The posterior uncertainty is calculated as the standard deviation of the ensemble members, where the posterior error covariance matrix is driven by the ensemble Kalman filter.

MIROC4-ACTM

The MIROC4-ACTM time-dependent inversions solve for emissions from 53 regions for CH4 and 84 regions for N2O. The inversion framework is based on Bayesian statistics and optimizes surface–atmosphere fluxes using the maximum probability solution. Atmospheric transport is modeled using JAMSTEC's Model for Interdisciplinary Research on Climate, version 4, based on an atmospheric chemistry–transport model (MIROC4-ACTM) (Watanabe et al., 2008; Patra et al., 2018). The source-receptor matrix (SRM) is calculated by simulating unitary emissions from 53 or 84 basis regions, for which the fluxes are optimized. The SRM describes the relationship between the change in mole fraction at the measurement locations for the unitary basis region fluxes. The MIROC4-ACTM meteorology was nudged to the JMA 55-year reanalysis (JRA-55) horizontal wind fields and temperature. The calculation of photochemical losses is performed online. The hydroxyl (OH) radical concentration for reaction with CH4 varies monthly but without any interannual variations. The simulated mole fractions for the total a priori fluxes are subtracted from the observed concentrations before running the inversion calculation (as in Patra et al., 2016, for CH4 inversion). Both inversion results are contributed to the GCP-CH4 and GCP-N2O activities (Saunois et al., 2020; Thompson et al., 2019; Tian et al., 2020).

Uncertainties. The posterior fluxes are subject to systematic errors primarily from (1) errors in the modeled atmospheric transport; (2) aggregation errors, i.e., errors arising from the way the flux variables are discretized in space (84 regions) and time (monthly means); (3) errors in the background mole fractions (assumed to be a minor factor); and (4) the incomplete information from the sparse observational network and hence the dependence on the prior fluxes. In addition, there is, to a much smaller extent, some error due to calibration offsets between observing instruments, which is more pertinent for N2O than for other GHGs. We have validated model transport in the troposphere using SF6 for the inter-hemispheric exchange time using SF6 and CO2 for the age of air in the stratosphere. The simulated N2O concentrations are also compared with aircraft measurements in the upper troposphere and lower stratosphere for evaluating the stratosphere–troposphere exchange rates. Comparisons with ACE-FTS vertical profiles in the stratosphere and mesosphere indicate good parameterization of N2O loss by photolysis and chemical reactions and thus the lifetime, which affect the global total N2O budgets. Random uncertainties are calculated by the inverse model depending on the prior flux uncertainties and the observational data density and data uncertainty. Only 37 sites are used in the inversion, and thus the reduction in priori flux uncertainties has been minimal. The net fluxes from the inversion from individual basis regions are less reliable compared to the anomalies in the estimated fluxes over a period of time.

Global Carbon Project – Global Methane Budget (GMB)

GMB uses an ensemble of 22 top-down global inversions for anthropogenic CH4 emissions presented in Saunois et al. (2020) for the Global Methane Budget. These inversions were simulated by nine atmospheric inversion systems based on various chemistry–transport models, differing in vertical and horizontal resolutions, meteorological forcing, advection and convection schemes, and boundary layer mixing. Surface-based inversions were performed over the period 2000–2017, while satellite-based inversions cover the GOSAT data availability over 2010–2017. The protocol established for these simulations was not stringent as the prior emission flux dataset was not mandatory, and each group selected their constraining observations. More information can be found in Saunois et al. (2020), in particular in their Tables 6 and S6.

Uncertainties. Currently there are no uncertainties reported for the GMB models. This study uses the median and the min/max as uncertainty range estimation from the 22-model ensemble. In general uncertainties might be due to factors like different transport models, physical parameterizations, prior fluxes and observation datasets.

A2 Natural CH4 emissions

A2.1 Bottom-up CH4 emission estimates

CH4 emissions from inland waters

The CH4 estimate from inland waters represents a climatology of average annual diffusive and ebulitive CH4 emissions from lakes and reservoirs at the spatial resolution of 0.1. The climatology is based on five alternative estimates, all relying on the high-resolution HydroLAKES database (Messager et al., 2016), and of which we report the mean and the standard deviation as a measure of uncertainty. Four of these estimates are based on predictions of CH4 emission rates from N and P concentrations. These concentrations were computed for each lake and reservoir of the HydroLAKES dataset (> 1.4 millions), using the mechanistic–stochastic model (MSM) of Maavara et al. (2017, 2019) and Lauerwald et al. (2019); see methodology for inland water N2O emissions for further details. The four estimates result from two empirical equations relating CH4 emissions to chlorophyll a concentrations (Deemer et al., 2016; DelSontro et al., 2018) and two equations relating chlorophyll a concentrations to nutrient concentrations (both from McCauley et al., 1989) in lakes and reservoirs. The fifth estimate is based on direct upscaling from observed CH4 emission rates (155 lakes and reservoirs), which we have classified into rates reported for small lakes (< 0.3 km2), larger (> 0.3 km2) lakes and reservoirs. In addition, we applied a coarse regionalization distinguishing the boreal (> 54 N) from the temperate to sub-tropical (< 54 N) zone.

JSBACH–HIMMELI

The model framework, JSBACH–HIMMELI (Raivonen et al., 2017; Susiluoto et al., 2018), is used to estimate wetland and mineral soil emissions, and an empirical model is used to estimate the emissions from inland water bodies.

JSBACH–HIMMELI is a combination of two models, JSBACH, which is the land-surface model of MPI-ESM (Reick et al., 2013), and HIMMELI, which is a specific model for northern peatland emissions of CH4 (Raivonen et al., 2017). HIMMELI (HelsinkI Model of MEthane buiLd-up and emIssion for peatlands) has been developed especially for estimating CH4 production and transport in northern peatlands. It simulates both CH4 and CO2 fluxes and can be used as a module within different modeling environments (Raivonen et al., 2017; Susiluoto et al., 2018). HIMMELI is driven with soil temperature, water table depth, the leaf area index and anoxic respiration. These parameters are provided to HIMMELI from JSBACH, which models hydrology, vegetation and soil carbon input from litter and root exudates. CH4 emission and uptake of mineral soils are calculated applying the method by Spahni et al. (2011) based on soil moisture estimated by JSBACH.

The distribution of terrestrial vegetation types in JSBACH–HIMMELI is adopted from CORINE land cover data and from native JSBACH land cover for the areas that CORINE does not cover. The HIMMELI methane model is applied for peatlands and the mineral soil approach for the rest. The map of inland water CH4 emissions has been combined with the JSBACH–HIMMELI land use map so that the map of inland waters is preserved and JSBACH grid-based fractions of different land use categories are adjusted accordingly. In order to avoid double-counting the terrestrial CH4 flux estimates have been normalized by the ratio of the two inland water body distributions.

Uncertainties. As in any process modeling, the uncertainties of the bottom-up modeling of CH4 arise from three primary sources: parameters, forcing data (including spatial and temporal resolution) and model structure. An important source of uncertainty in the case of terrestrial CH4 flux modeling is the spatial distribution of peatlands.

The uncertainties of JSBACH–HIMMELI peatland emissions were estimated by comparing the annual totals of measured and simulated methane fluxes at five European observation sites. Two of the sites are located in Finnish Lapland, one in middle Sweden, one in southern Finland and one in Poland.

For the sensitivity of mineral soil fluxes Spahni et al. (2011) tested two soil moisture thresholds, 85 % or 95 % of water holding capacity. Below these two thresholds, mineral soils are assumed to be only CH4 sinks, while above these thresholds they are assumed to be only CH4 sources. We used the higher value, 95 % of water holding capacity. The uncertainty was estimated using CH4 flux simulations of 1 year (2005). We performed two new model runs, using the moisture threshold 95 ± 15 %, and derived the uncertainty from the resulting range in the annual emission sum.

Geological fluxes

To calculate geological CH4 emissions we used literature data for geological emissions on land (excluding marine seepage) (Etiope et al., 2019; Hmiel et al., 2020). From the gridded geological CH4 emissions by Etiope et al. (2019), and using the land–sea mask for the EU27 + UK (to exclude marine seepage), geological CH4 emissions from the land of EU27 + UK are 8.83 Tg CH4 yr−1. Then we scaled this number by the ratio of global geological CH4 emissions estimated by Hmiel et al. (2020) and by Etiope et al. (2019), thus obtaining a value of 8.83×5.4/37.4= 1.3 Tg CH4 yr−1 (marine and land geological). The global total geological CH4 emissions reported by Etiope et al. (2019) and Hmiel et al. (2020) are 37.4 and 5.4 Tg CH4 yr−1, respectively.

A2.2 Top-down CH4 emission estimates

Global Carbon Project – Global Methane Budget (Saunois et al., 2020)

GMB uses an ensemble of 13 monthly gridded estimates of wetland emissions based on different land surface models as calculated for Saunois et al. (2020). Each model conducted a 30-year spin-up and then simulated net methane emissions from wetland ecosystems over 2000–2017. The models were forced by CRU-JRA reconstructed climate fields (Harris, 2019) and by the remote sensing-based wetland dynamical area dataset WAD2M (Wetland Area and Dynamics for Methane Modeling). This dataset provides monthly global areas over 2000–2017 based on a combination of microwave remote sensing data from Schroeder et al. (2015) and various regional inventory datasets. More information is available in Saunois et al. (2020), and more details will be presented in a future publication led by Poulter et al. (2017) and colleagues.

Uncertainty. As described by Saunois et al. (2020) uncertainties are reported as minimum and maximum values of the available studies in brackets. They do not take into account the uncertainty of the individual estimates but rather express the uncertainty as the range of available mean estimates, i.e., the standard error across measurements/methodologies considered.

A3 Anthropogenic and natural N2O emissions

A3.1 Bottom-up N2O emission estimates

UNFCCC NGHGI (2019), EDGAR v5.0 and CAPRI descriptions are found in Appendix A1.

ECOSSE

ECOSSE is a biogeochemical model that is based on the carbon model RothC (Jenkinson and Rayner, 1977; Jenkinson et al., 1987; Coleman and Jenkinson, 1996) and the nitrogen model SUNDIAL (Bradbury et al., 1993; Smith et al., 1996). All processes of the carbon and nitrogen dynamics are considered (Smith et al., 2010a,b). Additionally, in ECOSSE, processes of minor relevance for mineral arable soils (e.g., methane emissions) are implemented in order to represent better processes that are relevant for other soils (e.g., organic soils). ECOSSE can run in different modes and for different time steps. The two main modes are site-specific and limited data. In the later version, basis assumptions/estimates for parameters can be provided by the model. This increases the uncertainty but makes ECOSSE a universal tool that can be applied for large-scale simulations even if the data availability is limited. To increase the accuracy in the site-specific version of the model, detailed information about soil properties, plant input, nutrient application and management can be added as available.

During the decomposition process, material is exchanged between the soil organic matter (SOM) pools according to first-order rate equations, characterized by a specific rate constant for each pool, and modified according to rate modifiers dependent on the temperature, moisture, crop cover and pH of the soil. The N content of the soil follows the decomposition of the SOM, with a stable C:N ratio defined for each pool at a given pH, with N being either mineralized or immobilized to maintain that ratio. Nitrogen released from decomposing SOM as ammonium (NH4+) or added to the soil may be nitrified to nitrate (NO3-).

For spatial simulations the model is implemented in a spatial model platform. This allows us to aggregate the input parameter for the needed resolution. ECOSSE is a one-dimensional model, and the model platform provides the input data in a spatial distribution and aggregates the model outputs for further analysis. While climate data are interpolated, soil data are represented by the dominant soil type or by the proportional representation of the different soil types in the spatial simulation unit (this is a grid cell in VERIFY).

Uncertainties in ECOSSE arise from three primary sources: parameters, forcing data (including spatial and temporal resolution) and model structure.

DayCent

DayCent was designed to simulate soil C dynamics, nutrient flows (N, P, S) and trace gas fluxes (CO2, CH4, N2O, NOx, N2) between soil, plants and the atmosphere at a daily time step. Submodels include soil water content and temperature by layer, plant production and allocation of net primary production (NPP), decomposition of litter and soil organic matter, mineralization of nutrients, N gas emissions from nitrification and denitrification, and CH4 oxidation in non-saturated soils.

The DayCent modeling application at the EU level is a consolidated model framework running on Land Use/Cover Area frame statistical Survey (LUCAS) points (Orgiazzi et al., 2018), which was extensively explained in previous works (Lugato et al., 2017, 2018; Quemada et al., 2020) where a detailed description of numerical and geographical datasets and uncertainty estimations is reported.

Information directly derived from LUCAS (2009–2015) included the soil organic carbon content (SOC), particle size distribution and pH. Hydraulic properties and bulk density was also calculated with an empirically derived pedotransfer. Management information was derived from official statistics (EUROSTAT, https://ec.europa.eu/eurostat/web/agriculture/data, last access: Seprember 2019) and included crop shares at the NUTS2 level. The amount of mineral N was partitioned according to the regional crop rotations and agronomic crop requirements. Organic fertilization and irrigated areas were derived from the Gridded Livestock of the World FAO dataset and the FAO-AQUASTAT product.

Meteorological data were downloaded from the E-OBS gridded dataset (http://www.ecad.eu, last access: June 2019) at 0.1 resolution. For the climatic projection, the gridded data from CORDEX database (https://esgf-node.ipsl.upmc.fr/search/cordex-ipsl/, last access: June 2019) were used. The average annual (2006–2010) atmospheric N deposition from the EMEP model (rv 4.5) was also implemented into the simulations.

Uncertainty. The starting year of the simulation was set in 2009 and projected in the future. The uncertainty analysis, based on the Monte Carlo approach, was done running the model 52 times in each point and, contemporarily, randomly sampling model inputs from probability density functions for SOC pool partition, irrigation, and both mineral and organic fertilization rates. The model outputs (including uncertainties) at point level were upscaled regionally at 1 km resolution by a machine learning approach based on random forest regression.

N2O emissions from inland waters

The N2O estimate represents a climatology of average annual N2O emissions from rivers, lakes, reservoirs and estuaries at the spatial resolution of 0.1. Based on a spatially explicit representation of water bodies and point and non-point sources of N and P, this model quantifies the global-scale spatial patterns in inland water N2O emissions in a consistent manner at 0.5 resolution, which were then downscaled to 0.1 using the spatial distribution of European inland water bodies. The procedure to calculate the cascading loads of N and P delivered to each water body along the river–reservoir–estuary continuum and to topologically connect 1.4 million lakes (extracted from the HydroLAKES database) is described in Maavara et al. (2019) and Lauerwald et al. (2019). The methodology to quantify N2O emissions is based on the application of a mechanistic stochastic model (MSM) to estimate inland water C-N-P cycling as well as N2O production and emission generated by nitrification and denitrification. Using a Monte Carlo analysis, the MSM allows us to generate relationships relating N processes and N2O emissions to N and P loads and water residence time from the mechanistic model outputs, which are subsequently applied for the spatially resolved upscaling. For the estimation of N2O emission, we ran two distinct model configurations relying on EFs scaling to denitrification and nitrification rates: one assuming that N2O production equals N2O emissions and the other taking into account the kinetic limitation on N2O gas transfer and progressive N2O reduction to N2 during denitrification in water bodies with increasing residence time (Maavara et al., 2019). The model outputs from the two scenarios are used to constrain uncertainties in N2O emission estimates.

GAINS

Specific sectors and abatement technologies in GAINS vary by the specific emitted compound, with source sector definition and emission factors largely following the IPCC methodology at the Tier 1 or Tier 2 level. GAINS includes in general all anthropogenic emissions to air but does not cover emissions from forest fires, savannah burning and land use/land use change. Emissions are estimated for 174 countries/regions, with the possibility to aggregate to a global emission estimate and spanning a timeframe from 1990 to 2050 in 5-year intervals. Activity drivers for macroeconomic development, energy supply and demand, and agricultural activities are entered externally; GAINS extends with knowledge required to estimate “default” emissions (emissions occurring due to an economic activity without emission abatement) and emissions and costs of situations under emission control (Amann et al., 2011).

Emissions of nitrous oxide derive from energy, industry, agriculture and waste. Land use change emissions are not included. In the energy sector, certain technologies implemented to improve air quality affect N2O emission factors (like catalytic converters in vehicles), sometimes also negatively. That is also the case for non-selective catalytic reduction devices for NOx abatement in power plants or for fluidized bed combustion. Relevant industrial processes cover nitric acid and adipic acid, with other processes (glyoxal, if relevant, or caprolactam) included. Both processes allow for two different levels of abatement technologies, which both are relatively easily accessible and low cost. The use of N2O in gaseous form, often as an anesthetic for medical purposes, is associated with population numbers and scaled by availability of hospital beds. Marked emission reductions (at low costs) as well as complete phaseout of emissions (high costs) are implemented as technologies. Agricultural emissions in part derive from manure handling, where different management strategies have repercussions on emissions. The larger fraction of emission is from application of nitrogen compounds in different forms to grassland, crops and rice, with rice using a different emission factor. Application of manure and of mineral fertilizer in GAINS can be reduced by advanced computer technology such as automatic steering and variable rate application or by agrochemistry (nitrification inhibitors). The costs of implementation are considered to depend on the size of a farm; hence farm size is an important parameter. In the waste sector, composting and wastewater treatment are considered relevant sources. For wastewater treatment, GAINS also considers a specific emission reduction option when optimizing processes towards N2O reduction (e.g., via favoring the anammox process). All details have been reported by Winiwarter et al. (2018) in their supplementary material.

Uncertainties. The same paper provides full information on the uncertainty of N2O emissions in the GAINS model, which is a consequence of uncertainty provided in the activity data, in the emission factors and in the actual structure of the respective management strategies that also include the share of abatement technology already implemented. Further parameters also described (on uncertainty of future projections and on costs) are not relevant here.

FAOSTAT

The Statistics Division of the Food and Agriculture Organization of the United Nations provides N2O emissions from agriculture (http://www.fao.org/faostat/en/#data/GT/visualize, last access: January 2021) and its sub-domains, as well as N2O emissions from land use linked to biomass burning (metadata: http://fenixservices.fao.org/faostat/static/documents/GT/GT_e_2019.pdf, last access: December 2020, and http://fenixservices.fao.org/faostat/static/documents/GL/GL_e_2019.pdf, last access: December 2020). The FAOSTAT emissions database is computed following Tier 1 of the 2006 IPCC Guidelines for National Greenhouse Gas Inventories (http://www.ipcc-nggip.iges.or.jp/public/2006gl/index.html, last access: December 2020). Country reports to FAO on crops, livestock and agriculture use of fertilizers are the source of activity data. Geospatial data are the source of AD for the estimates from cultivation of organic soils, biomass and peat fires. N2O emissions are provided by country, regions and special groups, with global coverage, relative to the period 1961–present (with annual updates, currently 2017) and with projections for 2030 and 2050 for agriculture only, expressed in both CO2 eq. and N2O by underlying agricultural and land use emission sub-domain and by aggregate (agriculture total, agriculture total plus energy, agricultural soils). The main N2O emissions are reported for the following agricultural activities: manure management, synthetic fertilizers, manure applied to the soils, manure left in pasture, crop residues, cultivation of organic soils and burning crop residues. LULUCF emissions consist of N2O associated with burning biomass and peat fires, as well as from the drainage of organic soils.

Uncertainties were computed by Tubiello et al. (2013) but are not available in the FAOSTAT database.

A3.2 Top-down N2O emission estimates

FLEXINVERT_NILU

The FLEXINVERT_NILU framework is based on Bayesian statistics and optimizes surface–atmosphere fluxes using the maximum probability solution (Rodgers 2000). Atmospheric transport is modeled using the Lagrangian model FLEXPART (Stohl et al., 2005; Pisso et al., 2019) run in the backwards time mode to generate a so-called source-receptor matrix (SRM). The SRM describes the relationship between the change in mole fraction and the fluxes discretized in space and time (Seibert and Frank, 2004) and was calculated for 7 d prior to each observation. For use in the inversions, FLEXPART was driven using ECMWF ERA-Interim wind fields.

The state vector consisted of flux increments (i.e., offsets to the prior fluxes) discretized on an irregular grid based on the SRMs (Thompson et al., 2014). This grid has finer resolution (in this case the finest was 0.5× 0.5) where the fluxes have a strong influence on the observations and coarser resolution where the influence is only weak (the coarsest was 2× 2). The flux increments were solved at 2-weekly temporal resolution. The state vector also included scalars for the background mole fractions. The optimal (posterior) fluxes were found using the conjugate gradient method (e.g., Paige and Saunders, 1975).

The background mole fractions, i.e., the contribution to the modeled mole fractions that is not accounted for in the 7 d SRMs, was estimated by coupling the termination points of backwards trajectories (modeled using virtual particles) to initial fields of mole fractions from the optimized Eulerian model LMDZ (i.e., the CAMS N2O mole fraction product v18r1) following the method of Thompson et al. (2014).

Uncertainties. The posterior fluxes are subject to systematic errors primarily from (1) errors in the modeled atmospheric transport; (2) aggregation errors, i.e., errors arising from the way the flux variables are discretized in space and time; (3) errors in the background mole fractions; and (4) the incomplete information from the observations and hence the dependence on the prior fluxes. In addition, there is, to a smaller extent, some error due to calibration offsets between observing instruments, which is more pertinent for N2O than for other GHGs. Random uncertainties are calculated from a Monte Carlo ensemble of inversions following Chevallier et al. (2007), and uncertainties in the observation space were inflated to take into account the model representation errors.

A3.3 Global N2O Budget – GCP (Tian et al., 2020)

CAMS-N2O

Within the GCP 2019 results, N2O fluxes are estimated using the atmospheric inversion framework, CAMS-N2O. Atmospheric inversions use observations of atmospheric mixing ratios, in this case of N2O, and provide the fluxes that best explain the observations while at the same time being guided by a prior estimate of the fluxes. In other words, the fluxes are optimized to fit the observations within the limits of the prior and observation uncertainties. To produce the optimized (a posteriori) fluxes a number of steps are involved: first, the observations are pre-processed; second, a prior flux estimate is prepared; third, mixing ratios are simulated using the prior fluxes and are used to estimate the model representation error; and fourth, the inversion is performed.

In total 140 ground-based sites, ship and aircraft transects are included in the inversion. The term “site” refers to locations where there is a long-term record of observations and includes ground-based measurements, both from discrete samples (or “flasks”) and quasi-continuous sampling by in situ instruments, as well as aircraft measurements. A prior estimate of the total N2O flux with monthly resolution and inter-annually varying fluxes is prepared from a number of models and inventories. For the soil fluxes (including anthropogenic and natural) an estimate from the land surface model OCN-v1.2 is used, which is driven by observation-based climate data, N fertilizer statistics and modeled N deposition (Zaehle et al., 2011). For the ocean fluxes, an estimate from the ocean biogeochemistry model PlankTOM-v10.2 is used, which is a prognostic model (Buitenhuis et al., 2018). Atmospheric transport is modeled using an offline version of the Laboratoire de Météorologie Dynamique model, LMDZ5, which computes the evolution of atmospheric compounds using archived fields of winds, convection mass fluxes and planetary boundary layer (PBL) exchange coefficients that have been calculated using the online version nudged to ECMWF ERA-Interim winds.

CAMS-N2O uses the Bayesian inversion method to find the optimal fluxes of N2O given prior information about the fluxes and their uncertainty, as well as observations of atmospheric N2O mole fractions. The method is the same as that used in Thompson et al. (2014)

Uncertainty. Uncertainties in CAMS-N2O simulations pertain to observation space and to state space. The uncertainty in the observation space is calculated as the quadratic sum of the measurement and transport uncertainties. The measurement uncertainty is assumed to be 0.3 ppb (approximately 0.1 %) based on the recommendations of data providers. The transport uncertainty includes estimates of uncertainties in advective transport (based on the method of Rödenbeck et al., 2003) and from a lack of subgrid-scale variability (based on the method of Bergamaschi et al., 2010). For the error in each land grid cell, the maximum magnitude of the flux in the cell of interest and its eight neighbors is used, while for ocean grid cells the magnitude of the cell of interest only is used. Posterior flux uncertainties are calculated from a Monte Carlo ensemble of inversions, based on the method of Chevallier et al. (2005).

TOMCAT-INVICAT

TOMCAT-INVICAT (Wilson et al., 2014) is a variational inverse transport model, which is based on the global chemical transport model TOMCAT and its adjoint. It uses a four-dimensional variational (4DVAR) optimization framework based on Bayesian theory which seeks to minimize model–observation differences by altering surface fluxes while allowing for prior knowledge of these fluxes to be retained. TOMCAT (Monks et al., 2017) is an offline chemical transport model in which meteorological data are taken from ECMWF ERA-Interim reanalyses (Dee et al., 2011). The model grid resolution, and therefore the optimized surface flux estimates, has a horizontal resolution of 5.6× 5.6. The model has 60 vertical levels running from the surface to 0.1 hPa. For each individual year's fluxes, which are optimized on a monthly basis, 30 minimization iterations are carried out.

Uncertainty. Uncertainties in TOMCAT-INVICAT N2O inversions are described as follows and further in Thompson et al., (2019). The uncertainty in the observations is calculated as the quadratic sum of the measurement and transport uncertainties. The measurement uncertainty for each observation is assumed to be 0.4 ppb. The transport error for each observation is assumed to be the mean difference between the observation grid cell and its eight neighbors. Prior flux errors are assumed to be 100 % or the prior estimate and are uncorrelated in space and time. Posterior flux uncertainties are not currently able to be calculated.

MIROC4-ACTM

The MIROC4-ACTM time-dependent inversion for 84 regions (TDI84) framework is based on Bayesian statistics and optimizes surface–atmosphere fluxes using the maximum probability solution (Rodgers 2000). Atmospheric transport is modeled using JAMSTEC's Model for Interdisciplinary Research on Climate, version 4, based on an atmospheric chemistry–transport model (MIROC4-ACTM) (Watanabe et al., 2008; Patra et al., 2018). The source-receptor matrix (SRM) is calculated by simulating unitary emissions from 84 basis regions, for which the fluxes are optimized. The SRM describes the relationship between the change in mole fraction at the measurement locations for the unitary basis region fluxes (similar to Rayner et al., 1999). The MIROC4-ACTM meteorology was nudged to the JMA 55-year reanalysis (JRA-55) horizontal wind fields and temperature.

The simulated mole fractions for the total a priori fluxes are subtracted from the observed concentrations before running the inversion calculation (as in Patra et al., 2016, for CH4 inversion).

Uncertainties. The posterior fluxes are subject to systematic errors primarily from (1) errors in the modeled atmospheric transport; (2) aggregation errors, i.e., errors arising from the way the flux variables are discretized in space (84 regions) and time (monthly means); (3) errors in the background mole fractions (assumed to be a minor factor); and (4) the incomplete information from the sparse observational network and hence the dependence on the prior fluxes. In addition, there is, to a much smaller extent, some error due to calibration offsets between observing instruments, which is more pertinent for N2O than for other GHGs. We have validated model transport in the troposphere using SF6 for the inter-hemispheric exchange time and used SF6 and CO2 for the age of air in the stratosphere. The simulated N2O concentrations are also compared with aircraft measurements in the upper troposphere and lower stratosphere for evaluating the stratosphere–troposphere exchange rates. Comparisons with ACE-FTS vertical profiles in the stratosphere and mesosphere indicate good parameterization of N2O loss by photolysis and chemical reactions and thus the lifetime, which affect the global total N2O budgets.

Random uncertainties are calculated by the inverse model depending on the prior flux uncertainties and the observational data density and data uncertainty. Only 37 sites are used in the inversion, and thus the reduction in priori flux uncertainties has been minimal. The net fluxes from the inversion from individual basis regions are less reliable compared to the anomalies in the estimated fluxes over a period of time.

Appendix B

B1 Overview tables

Table B1Comparison of CH4 results from the BU and TD methods for common periods: BU anthropogenic 1990–2015, BU and TD natural 2005–2011, and TD total 2006–2012 representing the common period between all datasets and the last year available. All values are in kt CH4 yr−1, and uncertainties are in kilotons (kt). The UNFCCC NGHGI uncertainties represent the 95 % confidence interval; the uncertainty for EDGAR v5.0 was calculated for 2015, and the min/max values for all sectors are as follows: energy 33/37, IPPU 39/34, agriculture 18/18 and waste 32/38; it represents the 95 % confidence interval of a lognormal distribution. The other uncertainties represent the variability of the model ensembles (TD) and are the min and max of the averaged result over the time period. All values are rounded to the nearest 0.1 kt CH4, and therefore columns do not necessarily add up.

a The three periods were chosen based on the availability of data. The common period between all datasets is 2006–2012. b The last available year is listed in the following: UNFCCC NGHGI 2017, EDGAR v5.0 2015, GAINS 2015, CAPRI 2013, FAOSTAT 2017, JSBACH–HIMMELI 2017, Fluxes from lakes and reservoirs 2011, geological (one value for 2005–2017), GCP-CH4 natural wetlands partition from TD 2015, FLEXPART – FLExKF-TM5-4DVAR 2017, FLEXINVERT_NILU 2017, CTE-CH4 2017, InGOS 2012, GCP ensemble total 2017, total SURF 2017 and total GOSAT 2017. For details on model estimates and yearly values please download the data behind the figures on Zenodo https://doi.org/10.5281/zenodo.4590875. NA: not available.

Download Print Version | Download XLSX

Table B2Comparison of N2O results from the BU and TD methods for different periods: BU anthropogenic 1990–2015, TS total 2005–2014, and the common period between all datasets 2010–2014 and the last year available. All values are in kt N2O yr−1, and uncertainties are in kilotons (kt). The UNFCCC NGHGI uncertainties represent the 95 % confidence interval; the uncertainty for EDGAR v5.0 was calculated for 2015, and the min/max values for all sectors are as follows: energy 12/250, IPPU 13/19, agriculture 74/191 and waste 63/166; it represents the 95 % confidence interval of a lognormal distribution. The other uncertainties represent the variability of the model ensembles (TD) and are the min and max of the averaged result over the time period. All values are rounded to the nearest 0.1 kt N2O, and therefore columns do not necessarily add up.

a The three periods were chosen based on the availability of data. The common period between all datasets is 2010–2014. b The last available year is listed in the following: UNFCCC NGHGI (2019) 2017, EDGAR v5.0 2015, GAINS 2015, CAPRI 2013, FAOSTAT 2017, ECOSSE 2018, DayCent average 2011–2015, Lakes_rivers_reservoirs (one value 2010–2014), FLEXINVERT_NILU 2017, TOMCAT_LEEDS 2014, CAMS-N2O 2017, MIROC4-ACTM 2016. For details on model estimates and yearly values please download the data behind the figures on Zenodo https://doi.org/10.5281/zenodo.4590875. c GAINS reports one value for every 5 years. d UNFCCC uncertainties for agriculture and waste were set to 100 %; the values are much higher: 626 % for waste and 107 % for agriculture. e DayCent 2011–2015. NA: not available.

Download Print Version | Download XLSX

Table B3Adapted from Petrescu et al. (2020a): agriculture and natural N2O emissions – allocation of emissions to different sectors by different data sources.

NA: not available.

Download Print Version | Download XLSX

B2 Overview figures

https://essd.copernicus.org/articles/13/2307/2021/essd-13-2307-2021-f10

Figure B1(a) EU27 + UK total CH4 emissions time series per sector as reported by UNFCCC NGHGI (2019). (b) EU27 + UK total N2O emissions time series per sector as reported by UNFCCC NGHGI (2019).

Download

B3 Source-specific methodology: AD, EFs and uncertainties

Table B4Source-specific activity data (AD), emission factors (EFs) and uncertainty methodology for all current VERIFY and non-VERIFY 2019 data product collection.

NA: not available.

Download XLSX

Table B5Biogeochemical models that computed wetland emissions used in this study. Runs were performed for the whole period 2000–2017. Models run with prognostic (using their own calculation of wetland areas) and/or diagnostic (using WAD2M) wetland surface areas (see Sect. 3.2.1). From Saunois et al. (2020).

Download Print Version | Download XLSX

Table B6Top-down studies used in our new analysis, with their contribution to the decadal and yearly estimates noted. For decadal means, top-down studies have to provide at least 8 years of data over the decade to contribute to the estimate. From Saunois et al. (2020).

Download Print Version | Download XLSX

Table B7List of prior datasets for natural CH4 emissions used by all inverse models.

NA: not available.

Download XLSX

Author contributions

AMRP, AJD and PC designed research and led the discussions; AMRP wrote the initial draft of the paper and edited all the following versions; CQ made the figures; PC and CQ edited and commented to the first versions of the manuscript; GPP provided the Appendix B overview Fig. B1a and b and performed a very detailed review of a previous version; PP and MJM processed the original data submitted to the VERIFY portal; MJM, PP and PB designed and are managing the web portal; ES developed the methodology for the EDGAR uncertainty calculation and provided the CH4 and N2O uncertainties; RLT, GJM, FNT, PB, DBa, DBr, LHI, PR, RL, WW and AJD provided in-depth advice and commented on/edited the initial versions of the manuscript; LP and DG provided advice and edited the text related to the UNFCCC and NGHGI information; RLT, PB, DBr, LHI, PR, RL, AT, WW, PKP, MK, GDO, MC, MS, TM, TA, CDGZ, HT, YY, CW, AL, AL, JMH, AJM and JMN are data providers; all co-authors read and made changes to the subsequent versions of the manuscript and provided specific comments related to their data.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

FAOSTAT statistics are produced and disseminated with the support of its member countries to the FAO regular budget. The views expressed in this publication are those of the author(s) and do not necessarily reflect the views or policies of FAO.

We acknowledge the work of the entire EDGAR group (Marilena Muntean, Diego Guizzardi, Edwin Schaaf and Jos Olivier).

Financial support

Philippe Ciais received support of the European Research Council Synergy project SyG-2013-610028 IMBALANCE-P and from the ANR CLAND Convergence Institute. Prabir Patra received support of the Environment Research and Technology Development Fund (JPMEERF20172001, JPMEERF20182002) of the Environmental Restoration and Conservation Agency of Japan. David Basviken received support of the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement no. 725546). David Bastviken was supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant agreement no. 725546). Tuula Aalto received support from the Caroline Herschel Framework Partnership Agreement under the EU Horizon 2020 program (FPCUP, grant no. 809596) and Academy of Finland (SOMPA, grant no. 312932). Christine D. Groot Zwaaftink received support by the Norwegian Research Council (ICOS-Norway, project 245927). Joe McNorton received financial support from the Horizon2020 CHE project (776186).

Review statement

This paper was edited by Nellie Elguindi and reviewed by two anonymous referees.

References

Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O.: Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039–4072, https://doi.org/10.5194/acp-11-4039-2011, 2011. 

Amann, M., Bertok, I., Borken-Kleefeld, J., Cofala, J., Heyes, C., Höglund-Isaksson, L., Klimont, Z., Nguyen, B., Posch, M., Rafaj, P., Sandler, R., Schöpp, W., Wagner, F., and Winiwarter, W.: Cost-effective control of air quality and greenhouse gases in Europe: modeling and policy applications, Environmental Modelling and Software, 26, 1489–1501, 2011. 

Arnold, T., Manning, A. J., Kim, J., Li, S., Webster, H., Thomson, D., Mühle, J., Weiss, R. F., Park, S., and O'Doherty, S.: Inverse modelling of CF4 and NF3 emissions in East Asia, Atmos. Chem. Phys., 18, 13305–13320, https://doi.org/10.5194/acp-18-13305-2018, 2018. 

Arora, V. K., Melton, J. R., and Plummer, D.: An assessment of natural methane fluxes simulated by the CLASS-CTEM model, Biogeosciences, 15, 4683–4709, https://doi.org/10.5194/bg-15-4683-2018, 2018. 

Beaulieu, J. J., Tank, J. L., Hamilton, S. K., Wollheim, W. M., Hall, R. O., Mulholland, P. J., Peterson, B. J., Ashkenas, L. R., Cooper, L. W., Dahm, C. N., Dodds, W. K., Grimm, N. B., Johnson, S. L., McDowell, W. H., Poole, G. C. , Valett, H. M., Arango, C. P., Bernot, M. J., Burgin, A. J., Crenshaw, C. L., Helton, A. M., Johnson, L. T., O'Brien, J. M., Potter, J. D., Sheibley, R. W., Sobota, D. J., and Thomas, S. M.: Nitrous oxide emission from denitrification in stream and river networks, P. Natl. Acad. Sci. USA, 108, 214–219, https://doi.org/10.1073/pnas.1011464108, 2011. 

Berchet, A., Pison, I., Chevallier, F., Bousquet, P., Bonne, J.-L., and Paris, J.-D.: Objectified quantification of uncertainties in Bayesian atmospheric inversions, Geosci. Model Dev., 8, 1525–1546, https://doi.org/10.5194/gmd-8-1525-2015, 2015a. 

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

Bergamaschi, P., Frankenberg, C., Meirink, J. F., Krol, M., Dentener, F., Wagner, T., Platt, U., Kaplan, J. O., Koerner, S., Heimann, M., Dlugokencky, E. J., andGoede, A.: Satellite chartography of atmospheric methane from SCIAMACHY onboard ENVISAT: (II) Evaluation based on inverse model simulations, J. Geophys. Res., 112, D02304, https://doi.org/10.1029/2006JD007268, 2007. 

Bergamaschi, P., Krol, M., Meirink, J. F., Dentener, F., Segers, A., van Aardenne, J., Monni, S., Vermeulen, A., Schmidt, M., Ramonet, M., Yver, C., Meinhardt, F., Nisbet, E. G., Fisher, R., O'Doherty, S., and Dlugokencky, E. J.: Inverse modeling of European CH4 emissions 2001–2006, J. Geophys. Res., 115, D22309, https://doi.org/10.1029/2010JD014180, 2010. 

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

Bergamaschi, P., Corazza, M., Karstens, U., Athanassiadou, M., Thompson, R. L., Pison, I., Manning, A. J., Bousquet, P., Segers, A., Vermeulen, A. T., Janssens-Maenhout, G., Schmidt, M., Ramonet, M., Meinhardt, F., Aalto, T., Haszpra, L., Moncrieff, J., Popa, M. E., Lowry, D., Steinbacher, M., Jordan, A., O'Doherty, S., Piacentino, S., and Dlugokencky, E.: Top-down estimates of European CH4 and N2O emissions based on four different inverse models, Atmos. Chem. Phys., 15, 715–736, https://doi.org/10.5194/acp-15-715-2015, 2015. 

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

Bergamaschi, P., Danila, A. M., Weiss, R., Ciais, P., Thompson, R. L., Brunner, D., Levin, I., Meijer, Y., Chevallier, F., Janssens-Maenhout, G., Bovensmann, H., Crisp, D., Basu, S., Dlugokencky, E., Engelen, R., Gerbig, C., Günther, D., Hammer, S., Henne, S., Houweling, S., Karstens, U., Kort, E., Maione, M., Manning, A., Miller, J., Montzka, S., Pandey, S., Peters, W., Peylin, P., Pinty, B., Ramonet, M., Reimann, S., Röckmann, T., Schmidt, M., Strogies, M., Sussams, J., Tarasova, O., Van Aardenne, J., Vermeulen, A., and Vogel, F.: Atmospheric monitoring and inverse modelling for verification of greenhouse gas inventories, JRC report, https://doi.org/10.2760/759928, 2018b. 

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, https://doi.org/10.5194/gmd-10-2141-2017, 2017. 

Bousquet, P., Ringeval, B., Pison, I., Dlugokencky, E. J., Brunke, E.-G., Carouge, C., Chevallier, F., Fortems-Cheiney, A., Frankenberg, C., Hauglustaine, D. A., Krummel, P. B., Langenfelds, R. L., Ramonet, M., Schmidt, M., Steele, L. P., Szopa, S., Yver, C., Viovy, N., and Ciais, P.: Source attribution of the changes in atmospheric methane for 2006–2008, Atmos. Chem. Phys., 11, 3689–3700, https://doi.org/10.5194/acp-11-3689-2011, 2011. 

Bouwman, A., Beusen, A. H., and Billen, G.: Human alteration of the global nitrogen and phosphorus soil balances for the period 1970–2050, Global Biogeochem. Cy., 23, GB0A04, https://doi.org/10.1029/2009GB003576, 2009. 

Bradbury, N. J., Whitmore, A. P., Hart, P. B. S., and Jenkinson, D. S.: Modelling the fate of nitrogen in crop and soil in the years following application of 15N-labelled fertilizer to winter wheat, J. Agr. Sci., 121, 363–379, 1993. 

Britz, W. and Witzke, P.: CAPRI model documentation 2014, available at: https://www.capri-model.org/docs/capri_documentation.pdf (last access: September 2020), 2014. 

Brühl, C. and Crutzen, P. J.: MPIC Two-dimensional model, NASA Ref. Publ., 1292, 103–104, 1993. 

Brunner, D., Henne, S., Keller, C. A., Reimann, S., Vollmer, M. K., O'Doherty, S., and Maione, M.: An extended Kalman-filter for regional scale inverse emission estimation, Atmos. Chem. Phys., 12, 3455–3478, https://doi.org/10.5194/acp-12-3455-2012, 2012. 

Brunner, D., Arnold, T., Henne, S., Manning, A., Thompson, R. L., Maione, M., O'Doherty, S., and Reimann, S.: Comparison of four inverse modelling systems applied to the estimation of HFC-125, HFC-134a, and SF6 emissions over Europe, Atmos. Chem. Phys., 17, 10651–10674, https://doi.org/10.5194/acp-17-10651-2017, 2017. 

Buitenhuis, E. T., Suntharalingam, P., and Le Quéré, C.: Constraints on global oceanic emissions of N2O from observations and models, Biogeosciences, 15, 2161–2175, https://doi.org/10.5194/bg-15-2161-2018, 2018. 

Carlson, D. and Oda, T.: Editorial: Data publication – ESSDgoals, practices and recommendations, Earth Syst. Sci. Data, 10, 2275–2278, https://doi.org/10.5194/essd-10-2275-2018, 2018. 

Chevallier, F., Fisher, M., Peylin, P., Serrar, S., Bousquet, P., Bréon, F. M., Chédin, A., and Ciais, P.: Inferring CO2 sources and sinks from satellite observations: Method and application to TOVS data, J. Geophys. Res., 110, D24309, https://doi.org/10.1029/2005jd006390, 2005. 

Chevallier, F., Bréon, F. -M., and Rayner, P. J.: Contribution of the Orbiting Carbon Observatory to the estimation of CO2 sources and sinks: Theoretical study in a variational data assimilation framework, J. Geophys. Res., 112, D09307, https://doi.org/10.1029/2006jd007375, 2007. 

Coleman, K. and Jenkinson, D. S.: RothC-26.3 – A model the turnover of carbon in soil, in: Evaluation of soil organic matter models using existing long-term datasets, edited by: Powlson, D. S., Smith, P., and Smith, J. U., NATO ASI Series I, vol. 38, Springer, Berlin, 237–246, 1996. 

COM(96): 557 Strategy for reducing methane emissions, available at: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:51996DC0557&from=NL (last access: September 2020), 1996. 

CORDEX: Coordinated Regional Climate Downscaling Experiment, available at: https://esgf-node.ipsl.upmc.fr/search/cordex-ipsl/ (last access: June 2020), 2020. 

Crippa, M., Oreggioni, G., Guizzardi, D., Muntean, M., Schaaf, E., Lo Vullo, E., Solazzo, E., Monforti-Ferrario, F., Olivier, J. G. J., and Vignati, E.: Fossil CO2 and GHG emissions of all world countries – 2019 Report, EUR 29849 EN, Publications Office of the European Union, Luxembourg, ISBN 978-92-76-11100-9, https://doi.org/10.2760/687800, JRC117610, 2019. 

Crippa, M., Solazzo, E., Huang, G., Guizzardi, D., Koffi, E., Muntean, M., Schieberle, C., Friedrich, R., and Janssens-Maenhout, G.: High resolution temporal profiles in the Emissions Database for Global Atmospheric Research, Sci. Data, 17, 121, https://doi.org/10.1038/s41597-020-0462-2, 2020. 

CSIRO: Latest Cape Grim greenhouse gas data, available at: https://www.csiro.au/en/Research/OandA/Areas/Assessing-our-climate/Latest-greenhouse-gas-data (last access: October 2020), 2020. 

Cullen, M. J. P.: The unified forecast/climate model, Meteorol. Mag., 122, 81–94, 1993. 

Dee, D. P., Källén, E., Simmons, A. J., and Haimberger, L.: Comments on `Reanalyses suitable for characterizing long-term trends', B. Am. Meteorol. Soc., 92, 65–70, 2011. 

Deemer, B. R., Harrison, J. A., Li, S., Beaulieu, J. J., DelSontro, T., Barros, N., Bezerra-Neto, J. F., Powers, S. M., dos Santos, M. A., and Vonk, J. A.: Greenhouse Gas Emissions from Reservoir Water Surfaces: A New Global Synthesis, Bioscience, 66, 949–964, https://doi.org/10.1093/biosci/biw117, 2016. 

Del Sontro, T., Beaulieu, J. J., and Downing, J. A.: Greenhouse gas emissions from lakes and impoundments: Upscaling in the face of global change, Limnol. Oceanogr. Lett., 3, 64–75, https://doi.org/10.1002/lol2.10073, 2018. 

Directive 1999/31/EC, available at: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:31999L0031&from=EN (last access: June 2020), 1999. 

Dürr, H. H., Laruelle, G. G., van Kempen, C. M., Slomp, C. P., Meybeck, M., and Middlekoop, H.: Worldwide Typology of Nearshore Coastal Systems: Defining the Estuarine Filter of River Inputs to the Oceans, Estuar. Coast., 34, 441–458, https://doi.org/10.1007/s12237-011-9381-y, 2011. 

EEA: Approximated EU GHG inventory: proxy GHG estimates for 2018, EEA Report, available at: https://www.eea.europa.eu/publications/approximated-eu-ghg-inventory-proxy (last access: November 2020), 2019. 

Eisma, R., Vermeulen, A., and Van der Borg, K.: 14CH4 Emissions from Nuclear Power Plants in Northwestern Europe, Radiocarbon, 37, 475–483, 1995. 

Etheridge, D. M., Steele, L. P., Francey, R. J., and Langenfelds, R. L.: Atmospheric methane between 1000 A.D. and present: Evidence of anthropogenic emissions and climatic variability, J. Geophys. Res., 103, 15979–15993, https://doi.org/10.1029/98JD00923, 1998. 

Etiope, G.: Natural Gas Seepage, The Earth’s Hydrocarbon Degassing, Springer International Publishing, 2015. 

Etiope, G., and Milkov, A. V.: A new estimate of global methane flux from onshore and shallow submarine mud volcanoes to the atmosphere, Environ. Geol., 46, 997–1002, 2004. 

Etiope, G., Ciotoli, G., Schwietzke, S., and Schoell, M.: Gridded maps of geological methane emissions and their isotopic signature, Earth Syst. Sci. Data, 11, 1–22, https://doi.org/10.5194/essd-11-1-2019, 2019. 

EUROSTAT: Waste generation and treatment, avialable at: https://ec.europa.eu/eurostat/web/waste/data/database, (last access: June 2019), 2020. 

Evensen, G.: The Ensemble Kalman Filter: theoretical formulation and practical implementation, Ocean Dynam., 53, 343–367, https://doi.org/10.1007/s10236-003-0036-9, 2003. 

FAO: The State of Food and Agriculture: Social protection and agriculture: breaking the cycle of rural poverty, available at: http://www.fao.org/3/a-i4910e.pdf (last access: September 2020), Rome, Italy, 2015. 

FAO: The State of Food Security and Nutrition in the World: Transforming food systems for affordable healthy diets, Rome, Italy, available at: http://www.fao.org/3/ca9692en/CA9692EN.pdf. last access: October 2020. 

FAO-AQUASTAT: FAO's Global Information System on Water and Agriculture, available at: http://www.fao.org/aquastat/en/, last access: September 2020. 

FAOSTAT: Live animal numbers, crop production, total nitrogen fertiliser consumption statistics till 2016, Statistics Division of the Food and Agricultural Organisation of the UN, available at: http://www.fao.org/faostat/en/#home, last access: December, 2018. 

Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Hauck, J., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Bakker, D. C. E., Canadell, J. G., Ciais, P., Jackson, R. B., Anthoni, P., Barbero, L., Bastos, A., Bastrikov, V., Becker, M., Bopp, L., Buitenhuis, E., Chandra, N., Chevallier, F., Chini, L. P., Currie, K. I., Feely, R. A., Gehlen, M., Gilfillan, D., Gkritzalis, T., Goll, D. S., Gruber, N., Gutekunst, S., Harris, I., Haverd, V., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Joetzjer, E., Kaplan, J. O., Kato, E., Klein Goldewijk, K., Korsbakken, J. I., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Marland, G., McGuire, P. C., Melton, J. R., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Neill, C., Omar, A. M., Ono, T., Peregon, A., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Séférian, R., Schwinger, J., Smith, N., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Werf, G. R., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2019, Earth Syst. Sci. Data, 11, 1783–1838, https://doi.org/10.5194/essd-11-1783-2019, 2019. 

Fronzek, S., Pirttioja, N., Carter, T. R., Bindi, M., Hoffmann, H., Palosuo, T., Ruiz-Ramos, M., Tao, F., Trnka, M., Acutis, M., Asseng, S., Baranowski, P., Basso, B., Bodin, P.,Buis, S., Cammarano, D., Deligios, P., Destain, M.-F., Dumont, B., Ewert, F., Ferrise, R., François, L., Gaiser, T., Hlavinka, P., Jacquemin, I., Kersebaum, K. C., Kollas, C.,Krzyszczak, J., Lorite, I. J., Minet, J., Minguez, M. I., Montesino, M., Moriondo, M., Müller, C., Nendel, C., Öztürk, I., Perego, A., Rodríguez, A., Ruane, A. C., Ruget, F., Sanna, M., Semenov, M. A., Slawinski, C., Stratonovitch, P., Supit, I., Waha, K., Wang, E., Wu, L., Zhao, Z., and Rötter, R. P.: Classifying multi-model wheat yield impactresponsesurfacesshowingsensitivitytotemperatureandprecipitationchange, Agr. Syst., 159, 209–224, https://doi.org/10.1016/j.agsy.2017.08.004, 2018. 

GCP: Global Carbon Project: https://www.globalcarbonproject.org/methanebudget/index.htm (last access: December 2019), 2018. 

Gerbig, C., Lin, J. C., Wofsy, S. C., Daube, B. C., Andrews, A. E., Stephens, B. B., Bakwin, P. S., and Grainger, C. A.: Toward constraining regional-scale fluxes of CO2 with atmospheric observations over a continent: 2. Analysis of COBRA data using a receptor-oriented framework, J. Geophys. Res., 108, 4757, https://doi.org/10.1029/2003JD003770, 2003. 

Giglio, L., Randerson, J. T., and van der Werf, G. R.: Analysis of daily, monthly, and annual burned area using the fourth-generation global fire emissions database (GFED4), J. Geophys. Res.-Biogeosci., 118, 317–328, https://doi.org/10.1002/jgrg.20042, 2013. 

Gilbert, J. C. and Lemaréchal, C.: Some numerical experiments with variable storage Quasi-Newton algorithms, Math. Program., 45, 407-435, https://doi.org/10.1007/BF01589113, 1989. 

Gilbert, J. C. and Lemaréchal, C.: Some numerical experiments with variable storage quasi-Newton algorithms, IIASA Working Paper WP-88, A-2361, Laxenburg, 1989. 

Gómez-Sanabria, A., Höglund-Isaksson, L., Rafaj, P., and Schöpp, W.: Carbon in global waste and wastewater flows – its potential as energy source under alternative future waste management regimes, Adv. Geosci., 45, 105–113, https://doi.org/10.5194/adgeo-45-105-2018, 2018. 

Groot Zwaaftink, C. D., Henne, S., Thompson, R. L., Dlugokencky, E. J., Machida, T., Paris, J.-D., Sasakawa, M., Segers, A., Sweeney, C., and Stohl, A.: Three-dimensional methane distribution simulated with FLEXPART 8-CTM-1.1 constrained with observation data, Geosci. Model Dev., 11, 4469–4487, https://doi.org/10.5194/gmd-11-4469-2018, 2018. 

Gruber, W., Villez, K., Kipf, M., Wunderlin, P., Siegrist, H., Vogt, L., Joss, A.: N2O emission in full-scale wastewater treatment: Proposing a refined monitoring strategy, Sci. Total Environ., 699, 134–157, https://doi.org/10.1016/j.scitotenv.2019.134157, 2020. 

Harris, I. C.: CRU JRA v1.1: A forcings dataset of gridded land surface blend of Climatic Research Unit (CRU) and Japanese reanalysis (JRA) data, January 1901–December 2017, University of East Anglia Climatic Research Unit, Centre for Environmental Data Analysis, https://doi.org/10.5285/13f3635174794bb98cf8ac4b0ee8f4ed, 2019. 

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

Heimann, H. and Körner, S.: The global atmospheric tracer model TM3, Technical Report, Max-Planck-Institut für Biogeochemie, 5, 131 pp., 2003. 

Hmiel, B., Petrenko, V. V., Dyonisius, M. N., Buizert, C., Smith, A. M., Place, P. F., Harth, C., Beaudette, R., Hua, Q., Yang, B., Vimont, I., Michel, S. E., Severinghaus, J. P., Etheridge, D., Bromley, T., Schmitt, J., Faïn, X., Weiss, R. F., and Dlugokencky, E.: Preindustrial 14CH4 indicates greater anthropogenic fossil CH4 emissions, Nature, 578, 409–412, 2020. 

Höglund-Isaksson, L.: Global anthropogenic methane emissions 2005–2030: technical mitigation potentials and costs, Atmos. Chem. Phys., 12, 9079–9096, https://doi.org/10.5194/acp-12-9079-2012, 2012. 

Höglund-Isaksson, L.: Bottom-up simulations of methane and ethane from global oil and gas systems, Environ. Res. Lett., 12, 024007, https://doi.org/10.1088/1748-9326/aa583e, 2017. 

Höglund-Isaksson, L., Winiwarter, W., Purohit, P., Gomez-Sanabria, A., Rafaj, P., Schopp, W., and Borken-Kleefeld, J.: Non-CO2 greenhouse gas emissions in the EU-28 from 2005 to 2070 with mitigation potentials and costs – GAINS model methodology, IIASA, 2018. 

Höglund-Isaksson L., Gómez-Sanabria, A., Klimont, Z., Rafaj, P., and Schöpp, W.: Technical potentials and costs for reducing global anthropogenic methane emissions in the 2050 timeframe – results from the GAINS model, Environ. Res. Commun., 2, 025004, https://doi.org/10.1088/2515-7620/ab7457, 2020. 

Hourdin, F. and Armengaud, A.: The use of finite-volume methods for atmospheric advection of trace species. Part I: test of vairious formulations in a general circulation model, Mon. Weather Rev., 127, 822–837, 1999. 

Hourdin, F., Musat, I., Bony, S., Braconnot, P., Codron, F., Dufresne, J. L., Fairhead, L., Filiberti, M. A., Friedlingstein, P., Grandpeix, J. Y., Krinner, G., LeVan, P., and Lott, F..: The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection, Clim. Dynam., 27, 787–813, https://doi.org/10.1007/s00382-006-0158-0, 2006. 

Houweling, S., Krol, M., Bergamaschi, P., Frankenberg, C., Dlugokencky, E. J., Morino, I., Notholt, J., Sherlock, V., Wunch, D., Beck, V., Gerbig, C., Chen, H., Kort, E. A., Röckmann, T., and Aben, I.: A multi-year methane inversion using SCIAMACHY, accounting for systematic errors using TCCON measurements, Atmos. Chem. Phys., 14, 3991–4012, https://doi.org/10.5194/acp-14-3991-2014, 2014. 

IEA: Energy balance statistics for 1970–2015, available at: http://www.iea.org/ (last access: March 2020), 2017. 

IPCC: Guidelines for National Greenhouse Gas Inventories (NGHGI), available at: https://www.ipcc-nggip.iges.or.jp/public/2006gl/ (last access: January 2020), 2006. 

IPCC: Supplement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories: Wetlands, edited by: Hiraishi, T., Krug, T., Tanabe, K., Srivastava, N., Baasansuren, J., Fukuda, M., and Troxler, T. G., IPCC, Switzerland, 2014. 

IPCC: Special Report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems (SRCCL), IPCC report, available at: https://www.ipcc.ch/report/srccl/ (last access: December 2019), 2019a. 

IPCC: 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories, available at: https://www.ipcc.ch/report/2019-refinement-to-the-2006-ipcc-guidelines-for-national-greenhouse-gas-inventories (last access: January 2020), 2019b. 

Ishizawa, M., Mabuchi, K., Shirai, T., Inoue, M., Morino, I., Uchino, O., Yoshida, Y., Maksyutov, S., and Belikov, D.: Inter-annual variability of CO2 exchange in Northern Eurasia inferred from GOSAT XCO2, Environ. Res. Lett., 11, 105001, https://doi.org/10.1088/1748-9326/11/10/105001, 2016. 

Ito, A. and Inatomi, M.: Use of a process-based model for assessing the methane budgets of global terrestrial ecosystems and evaluation of uncertainty, Biogeosciences, 9, 759–773, https://doi.org/10.5194/bg-9-759-2012, 2012. 

Janssens-Maenhout, G., Pagliari, V., Guizzardi, D., and Muntean, M.: Global emission inventories in the Emission Database for Global Atmospheric Research (EDGAR) – Manual (I) I. Gridding: EDGAR emissions distribution on global gridmaps, EUR – Scientific and Technical Research Reports, Publications Office of the European Union, available at: http://publications.jrc.ec.europa.eu/repository/bitstream/JRC78261/edgarv4_manual_i_gridding_pubsy_final.pdf (last access: 7 February 2018), 2013. 

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, https://doi.org/10.5194/essd-11-959-2019, 2019. 

Jenkinson, D. S. and Rayner, J. H.: The turnover of organic matter in some of the Rothamsted classical experiments, Soil Sci., 123, 298–305, 1977. 

Jenkinson, D. S., Hart, P. B. S., Rayner, J. H., and Parry, L. C.: Modelling the turnover of organic matter in long-term experiments at Rothamste, INTECOL Bulletin, 15, 1–8, 1987. 

Jones, A., Thomson, D., Hort, M., and Devenish, B.: The UK Met Office's next-generation atmospheric dispersion model, NAME III, in: Air Pollution Modeling and Its Applications XVII, edited by: Borrego, C. and Norman, A. L., Springer Science+Business Media, 580–589, 2007. 

Kaplan, J. O.: Wetlands at the last glacial maximum: distribution and methane emissions, Geophys. Res. Lett., 29, 3-1–3-4, https://doi.org/10.1029/2001GL013366, 2002. 

Kleinen, T., Brovkin, V., and Schuldt, R. J.: A dynamic model of wetland extent and peat accumulation: results for the Holocene, Biogeosciences, 9, 235–248, https://doi.org/10.5194/bg-9-235-2012, 2012. 

Kleinen, T., Mikolajewicz, U., and Brovkin, V.: Terrestrial methane emissions from the Last Glacial Maximum to the preindustrial period, Clim. Past, 16, 575–595, https://doi.org/10.5194/cp-16-575-2020, 2020. 

Krol, M., Houweling, S., Bregman, B., van den Broek, M., Segers, A., van Velthoven, P., Peters, W., Dentener, F., and Bergamaschi, P.: The two-way nested global chemistry-transport zoom model TM5: algorithm and applications, Atmos. Chem. Phys., 5, 417–432, https://doi.org/10.5194/acp-5-417-2005, 2005. 

Lambert, G. and Schmidt, S.: Reevaluation of the oceanic flux of methane: Uncertainties and long term variations, Chemosphere, 26, 579–589, https://doi.org/10.1016/0045-6535(93)90443-9, 1993. 

Lauerwald, R., Regnier, P., Figueiredo, V., Enrich-Prast, A., Bastviken, D., Lehner, B., Maavara, T., and Raymond, P.: Natural lakes are a minor global source of N2O to the atmosphere, Global Biogeochem. Cy., 33, 1564–1581, https://doi.org/10.1029/2019GB006261, 2019. 

Lehner, B., Verdin, K., and Jarvis, A.: New global hydrography derived from spaceborne elevation data, EOS T. Am. Geophys. Un., 89, 93–94, https://doi.org/10.1029/2008EO100001, 2008. 

Lin, J. C., Gerbig, C., Wofsy, S. C., Andrews, A. E., Daube, B. C., Davis, K. J., and Grainger, C. A.: A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model, J. Geophys. Res., 108, 4493, https://doi.org/10.1029/2002JD003161, 2003. 

Lugato, E., Paniagua, L., Jones, A., De Vries, W., and Leip, A.: Complementing the topsoil information of the Land Use/Land Cover Area Frame Survey (LUCAS) with modelled N2O emissions, PLoS One, 12, https://doi.org/10.1371/journal.pone.0176111, 2017. 

Lugato, E., Leip, A., and Jones, A.: Mitigation potential of soil carbon management overestimated by neglecting N2O emissions, Nat. Clim. Change, 8, 219–223, https://doi.org/10.1038/s41558-018-0087-z, 2018. 

Maavara, T., Lauerwald, R., Regnier, P., and Van Cappellen, P.: Global perturbation of organic carbon cycling by river damming, Nat. Commun., 8, 15347, https://doi.org/10.1038/ncomms15347, 2017. 

Maavara, T., Lauerwald, R., Laruelle, G., Akbarzadeh, Z., Bouskill, N., Van Cappellen, P., and Regnier, P.: Nitrous oxide emissions from inland waters: Are IPCC estimates too high?, Glob. Change Biol., 25, 473–488, https://doi.org/10.1111/gcb.14504, 2019. 

Maksyutov, S., Oda, T., Saito, M., Janardanan, R., Belikov, D., Kaiser, J. W., Zhuravlev, R., Ganshin, A., Valsala, V. K., Andrews, A., Chmura, L., Dlugokencky, E., Haszpra, L., Langenfelds, R. L., Machida, T., Nakazawa, T., Ramonet, M., Sweeney, C., and Worthy, D.: Technical note: A high-resolution inverse modelling technique for estimating surface CO2 fluxes based on the NIES-TM–FLEXPART coupled transport model and its adjoint, Atmos. Chem. Phys., 21, 1245–1266, https://doi.org/10.5194/acp-21-1245-2021, 2021. 

Manning, A. J., O'Doherty, S., Jones, A. R., Simmonds, P. G., and Derwent, R. G.: Estimating UK methane and nitrous oxide emissions from 1990 to 2007 using an inversion modeling approach, J. Geophys. Res.-Atmos., 116, D02305, https://doi.org/10.1029/2010jd014763, 2011. 

Matthews, E. and Fung, I.: Methane emission from natural wetlands: Global distribution, area, and environmental characteristics of sources, Global Biogeochem. Cy., 1, 61–86, https://doi.org/10.1029/GB001i001p00061, 1987. 

McCauley, E., Downing, J. A., and Watson, S.: Sigmoid Relationships between Nutrients and Chlorophyll among Lakes, Can. J. Fish. Aquat. Sci., 46, 1171–1175, https://doi.org/10.1139/f89-152, 1989. 

McGuire, A. D., Christensen, T. R., Hayes, D., Heroult, A., Euskirchen, E., Kimball, J. S., Koven, C., Lafleur, P., Miller, P. A., Oechel, W., Peylin, P., Williams, M., and Yi, Y.: An assessment of the carbon balance of Arctic tundra: comparisons among observations, process models, and atmospheric inversions, Biogeosciences, 9, 3185–3204, https://doi.org/10.5194/bg-9-3185-2012, 2012. 

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, https://doi.org/10.5194/acp-18-18149-2018, 2018. 

Meirink, J. F., Bergamaschi, P., and Krol, M. C.: Four-dimensional variational data assimilation for inverse modelling of atmospheric methane emissions: method and comparison with synthesis inversion, Atmos. Chem. Phys., 8, 6341–6353, https://doi.org/10.5194/acp-8-6341-2008, 2008. 

Melton, J. R. and Arora, V. K.: Competition between plant functional types in the Canadian Terrestrial Ecosystem Model (CTEM) v. 2.0, Geosci. Model Dev., 9, 323–361, https://doi.org/10.5194/gmd-9-323-2016, 2016. 

Menut, L., Bessagnet, B., Khvorostyanov, D., Beekmann, M., Blond, N., Colette, A., Coll, I., Curci, G., Foret, G., Hodzic, A., Mailler, S., Meleux, F., Monge, J.-L., Pison, I., Siour, G., Turquety, S., Valari, M., Vautard, R., and Vivanco, M. G.: CHIMERE 2013: a model for regional atmospheric composition modelling, Geosci. Model Dev., 6, 981–1028, https://doi.org/10.5194/gmd-6-981-2013, 2013. 

Messager, M. L., Lehner, B., Grill, G., Nedeva, I., and Schmitt, O.: Estimating the volume and age of water stored in global lakes using a geo-statistical approach, Nat. Commun., 7, 13603, https://doi.org/10.1038/ncomms13603, 2016. 

Monks, S. A., Arnold, S. R., Hollaway, M. J., Pope, R. J., Wilson, C., Feng, W., Emmerson, K. M., Kerridge, B. J., Latter, B. L., Miles, G. M., Siddans, R., and Chipperfield, M. P.: The TOMCAT global chemical transport model v1.6: description of chemical mechanism and model evaluation, Geosci. Model Dev., 10, 3025–3057, https://doi.org/10.5194/gmd-10-3025-2017, 2017. 

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

Nisbet, E. G., Manning, M. R., Dlugokencky, E. J., Fisher, R. E., Lowry, D., Michel, S. E., und Myhre, C., Platt, S. M., Allen, G., Bousquet, P., Brownlow, R., Cain, M., France, J. L., Hermansen, O., Hossaini, R., Jones, A. E., Levin, I., Manning, A. C., Myhre, G., Pyle, J. A., Vaughn, B. 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, https://doi.org/10.1029/2018GB006009, 2019. 

Niwa, Y., Tomita, H., Satoh, M., Imasu, R., Sawa, Y., Tsuboi, K., Matsueda, H., Machida, T., Sasakawa, M., Belan, B., and Saigusa, N.: A 4D-Var inversion system based on the icosahedral grid model (NICAM-TM 4D-Var v1.0) – Part 1: Offline forward and adjoint transport models, Geosci. Model Dev., 10, 1157–1174, https://doi.org/10.5194/gmd-10-1157-2017, 2017. 

Niwa, Y., Fujii, Y., Sawa, Y., Iida, Y., Ito, A., Satoh, M., Imasu, R., Tsuboi, K., Matsueda, H., and Saigusa, N.: A 4D-Var inversion system based on the icosahedral grid model (NICAM-TM 4D-Var v1.0) – Part 2: Optimization scheme and identical twin experiment of atmospheric CO2 inversion, Geosci. Model Dev., 10, 2201–2219, https://doi.org/10.5194/gmd-10-2201-2017, 2017b. 

NOAA: Atmospheric data, available at: https://www.esrl.noaa.gov/gmd/ccgg/trends_ch4/ (last access: June 2020), 2020. 

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

Olczak, M. and Piebalgs, A: How far should the new EU Methane Strategy go?, European University Institute, https://doi.org/10.2870/621991, 2019. 

Orgiazzi, A., Ballabio, C., Panagos, P., Jones, A., Fernández-Ugalde, O.: LUCAS Soil, the largest expandable soil dataset for Europe: a review, Eur. J. Soil Sci., 69, 140–153, 2018. 

Paige, C. C. and Saunders, M. A.: Solution of sparse indefinite systems of linear equations. Solution of sparse indefinite systems of linear equations, SIAM J. Numer. Anal., 12, 617–629, https://doi.org/10.1137/0712047, 1975. 

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, https://doi.org/10.5194/acp-16-5043-2016, 2016. 

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

Patra, P. K., Saeki, T., Dlugokencky, E. J., Ishijima, K., Umezawa, T., Ito, A., Aoki, S., Morimoto, S., Kort, E. A., Crotwell, A., Ravikumar, K., and Nakazawa, T.: Regional methane emission estimation based on observed atmospheric concentrations (2002–2012), J. Meteorol. Soc. Jpn., 94, 91–113, 2016. 

Patra, P. K., Takigawa, M., Watanabe, S., Chandra, N., Ishijima, K., and Yamashita, Y.: Improved Chemical Tracer Simulation by MIROC4.0-based Atmospheric Chemistry-Transport Model (MIROC4-ACTM), SOLA, 14, 91–96, 2018. 

Peters, W., Miller, J. B., Whitaker, J., Denning, A. S., Hirsch, A., Krol, M. C., Zupanski, D., Bruhwiler, L., and Tans, P. P.: An ensemble data assimilation system to estimate CO2 surface fluxes from atmospheric trace gas observations, J. Geophys. Res., 110, D24304, https://doi.org/10.1029/2005JD006157, 2005. 

Petrescu, A. M. R., Peters, G. P., Janssens-Maenhout, G., Ciais, P., Tubiello, F. N., Grassi, G., Nabuurs, G.-J., Leip, A., Carmona-Garcia, G., Winiwarter, W., Höglund-Isaksson, L., Günther, D., Solazzo, E., Kiesow, A., Bastos, A., Pongratz, J., Nabel, J. E. M. S., Conchedda, G., Pilli, R., Andrew, R. M., Schelhaas, M.-J., and Dolman, A. J.: European anthropogenic AFOLU greenhouse gas emissions: a review and benchmark data, Earth Syst. Sci. Data, 12, 961–1001, https://doi.org/10.5194/essd-12-961-2020, 2020a. 

Petrescu, A. M. R., Qiu, C., Ciais, P., Thompson, R. L., Peylin, P., McGrath, M. J., Solazzo, E., Janssens-Maenhout, G, Tubiello, F. N., Bergamaschi, P., Brunner, D., Peters, G. P., Höglund-Isaksson, L., Regnier, P., Lauerwald, R., Bastviken, D., Tsuruta, A., Winiwarter, W., Patra, P. K., Kuhnert, M., Orregioni, G. D., Crippa, M., Saunois, M., Perugini, L., Markkanen, T., Aalto, T., Groot Zwaaftink, C. D., Yao, Y., Wilson, C., Conchedda, G., Günther, D., Leip, A., Smith, P., Haussaire, J.-M., Leppänen, A., Manning, A. J., McNorton, J., Brockmann, P., and Dolman, A. J.: The consolidated European synthesis of CH4 and N2O emissions for EU27 and UK: 1990–2017, version 2, Zenodo, https://doi.org/10.5281/zenodo.4590875, 2020b. 

Petrescu, A. M. R., McGrath, M. J., Andrew, R. M., Peylin, P., Peters, G. P., Ciais, P., Broquet, G., Tubiello, F. N., Gerbig, C., Pongratz, J., Janssens-Maenhout, G., Grassi, G., Nabuurs, G.-J., Regnier, P., Lauerwald, R., Kuhnert, M., Balcovič, J., Schelhaas, M.-J., Denier van der Gon, H. A. C., Solazzo, E., Qiu, C., Pilli, R., Konovalov, I. B., Houghton, R. A., Günther, D., Perugini, L., Crippa, M., Ganzenmüller, R., Luijkx, I. T., Smith, P., Munassar, S., Thompson, R. L., Conchedda, G., Monteil, G., Scholze, M., Karstens, U., Brokmann, P., and Dolman, A. J.: The consolidated European synthesis of CO2 emissions and removals for the European Union and United Kingdom: 1990–2018, Earth Syst. Sci. Data, 13, 2363–2406, https://doi.org/10.5194/essd-13-2363-2021, 2021. 

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, https://doi.org/10.5194/gmd-12-4955-2019, 2019. 

Poulter, B., Bousquet, P., Canadell, J. G., Ciais, P., Peregon, A., Saunois, M., Arora, V. K., Beerling, D. J., Brovkin, V., Jones, C. D., Joos, F., Gedney, N., Ito, A., Kleinen, T.., Koven, C. D., McDonald, K., Melton, J. R., Peng, C., Prigent, C., Schroeder, R., Riley, W. J., Saito, M., Spahni, R., Tian, H., Taylor, L., Viovy, N., Wilton, D., Wiltshire, A., Xu, X., Zhang, B., Zhang, Z., and Zhu, Q.: Global wetland contribution to 2000–2012 atmospheric methane growth rate dynamics, Environ. Res. Lett., 12, 094013, https://doi.org/10.1088/1748-9326/aa8391, 2017. 

Quemada, M., Lassaletta, L., Leip, A., Jones, A., and Lugato, E.: Integrated management for sustainable cropping systems: looking beyond the greenhouse balance at the field scale, Glob. Change Biol., 26, 2584–2598, 2020. 

Raivonen, M., Smolander, S., Backman, L., Susiluoto, J., Aalto, T., Markkanen, T., Mäkelä, J., Rinne, J., Peltola, O., Aurela, M., Lohila, A., Tomasic, M., Li, X., Larmola, T., Juutinen, S., Tuittila, E.-S., Heimann, M., Sevanto, S., Kleinen, T., Brovkin, V., and Vesala, T.: HIMMELI v1.0: HelsinkI Model of MEthane buiLd-up and emIssion for peatlands, Geosci. Model Dev., 10, 4665–4691, https://doi.org/10.5194/gmd-10-4665-2017, 2017. 

Rayner, P. J., Enting, I. G., Francey, R. J., and Langenfelds, R.: Reconstructing the recent carbon cycle from atmospheric CO2, δ13C and O2/N2 observations, Tellus B, 51, 213–232, 1999. 

Regulation (EU): 2018/1999 (European Climate Law) amended proposal for a regulation of the European parliament and of the council on establishing the framework for achieving climate neutrality, available at: https://ec.europa.eu/clima/sites/clima/files/eu-climate-action/docs/prop_reg_ecl_en.pdf, last access: October 2020. 

Regulation (EU): 525/2013 of the European Parliament and of the Council, available at: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32013R0525&from=EN, last access: November 2020. 

Reick, C., Raddatz, T., Brovkin, V., and Gayler, V.: Representation of natural and anthropogenic land cover change in MPI-ESM, J. Adv. Model. Earth Sy., 5, 459–482, https://doi.org/10.1002/jame.20022, 2013. 

Riley, W. J., Subin, Z. M., Lawrence, D. M., Swenson, S. C., Torn, M. S., Meng, L., Mahowald, N. M., and Hess, P.: Barriers to predicting changes in global terrestrial methane fluxes: analyses using CLM4Me, a methane biogeochemistry model integrated in CESM, Biogeosciences, 8, 1925–1953, https://doi.org/10.5194/bg-8-1925-2011, 2011. 

Ringeval, B., Friedlingstein, P., Koven, C., Ciais, P., de Noblet-Ducoudré, N., Decharme, B., and Cadule, P.: Climate−CH4 feedback from wetlands and its interaction with the climate−CO2 feedback, Biogeosciences, 8, 2137–2157, https://doi.org/10.5194/bg-8-2137-2011, 2011. 

Rödenbeck, C.: Estimating CO2 sources and sinks from atmospheric mixing ratio measurements using a global inversion of atmospheric transport, Tech. Rep. 6, Max-Planck-Institut für Biogeochemie, Jena, available at: http://www.bgc-jena.mpg.de/~christian.roedenbeck/download/2005-Roedenbeck-TechReport6.pdf (last access: June 2020), 2005. 

Rödenbeck, C., Houweling, S., Gloor, M., and Heimann, M.: CO2 flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport, Atmos. Chem. Phys., 3, 1919–1964, https://doi.org/10.5194/acp-3-1919-2003, 2003. 

Rödenbeck, C., Gerbig, C., Trusilova, K., and Heimann, M.: A two-step scheme for high-resolution regional atmospheric trace gas inversions based on independent models, Atmos. Chem. Phys., 9, 5331–5342, https://doi.org/10.5194/acp-9-5331-2009, 2009. 

Rodgers, C: Inverse methods for atmospheric sounding: theory and practice, World Scientific, Singapore, https://doi.org/10.1142/3171, 2000. 

Sanderson, M. G.: Biomass of termites and their emissions of methane and carbon dioxide: A global database, Global Biogeochem. Cy., 10, 543–557, https://doi.org/10.1029/96gb01893, 1996. 

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, https://doi.org/10.5194/essd-8-697-2016, 2016. 

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, https://doi.org/10.5194/essd-12-1561-2020, 2020. 

Schroeder, R., McDonald, K. C., Chapman, B., Jensen, K., Podest, E., Tessler, Z., Bohn, T. J., and Zimmerman, R.: Development and evaluation of a multi-year inundated land surface data set derived from active/passive microwave remote sensing data, Remote Sens., 7, 16668–16732, https://doi.org/10.3390/rs71215843, 2015. 

Schwietzke, S., Sherwood, O. A., Bruhwiler, L. M. P., Miller, J. B., Etiope, G., Dlugokencky, E. J., Michel, S. E., Arling, V. A., Vaughn, B. H., White, J. W. C., and Tans, P. P.: Upward revision of global fossil fuel methane emissions based on isotope database, Nature, 538, 88–91, https://doi.org/10.1038/nature19797, 2016. 

Segers, A. and Houweling, S.: Description of the CH4 Inversion Production Chain, CAMS (Copernicus Atmospheric Monitoring Service) Report, available at: https://atmosphere.copernicus.eu/sites/default/files/2018-11/CAMS73_2015SC3_D73.2.5.5-2018_201811_production_chain_v1_0.pdf (last access: March 2020), 2018. 

Seibert, P. and Frank, A.: Source-receptor matrix calculation with a Lagrangian particle dispersion model in backward mode, Atmos. Chem. Phys., 4, 51–63, https://doi.org/10.5194/acp-4-51-2004, 2004. 

Silpa, K., Yao, L. C., Bhada-Tata, P., and Van Woerden, F.: What a Waste 2.0: A Global Snapshot of Solid Waste Management to 2050. Urban Development, World Bank, Washington, DC, available at: https://openknowledge.worldbank.org/handle/10986/30317 (last access: September 2019), 2018. 

Smith, J. U., Bradbury, N. J., and Addiscott, T. M.: SUNDIAL: A PC-based system for simulating nitrogen dynamics in arable land, Agron J, 88, 38–43, 1996. 

Smith, J. U., Gottschalk, P., Bellarby, J., Chapman, S., Lilly, A., Towers, W., Bell, J., Coleman, K., Nayak, D. R., Richards, M. I., Hillier, J., Flynn, H. C., Wattenbach, M., Aitkenhead, M., Yeluripurti, J. B., Farmer, J., Milne, R., Thomson, A., Evans, C., Whitmore, A. P., Falloon, P., and Smith, P.: Estimating changes in national soil carbon stocks using ECOSSE – a new model that includes upland organic soils. Part I. Model description and uncertainty in national scale simulations of Scotland, Clim. Res., 45, 179–192, https://doi.org/10.3354/cr00899, 2010a. 

Smith, J. U., Gottschalk, P., Bellarby, J., Chapman, S., Lilly, A., Towers, W., Bell, J., Coleman, K., Nayak, D. R., Richards, M. I., Hillier, J., Flynn, H. C., Wattenbach, M., Aitkenhead, M., Yeluripurti, J. B., Farmer, J., Milne, R., Thomson, A., Evans, C., Whitmore, A. P., Falloon, P., and Smith, P.: Estimating changes in national soil carbon stocks using ECOSSE – a new model that includes upland organic soils. Part II. Application in Scotland, Clim. Res., 45, 193–205, https://doi.org/10.3354/cr00902, 2010b. 

Solazzo, E., Crippa, M., Guizzardi, D., Muntean, M., Choulga, M., and Janssens-Maenhout, G.: Uncertainties in the EDGAR emission inventory of greenhouse gases, Atmos. Chem. Phys. Discuss. [preprint], https://doi.org/10.5194/acp-2020-1102, in review, 2020. 

Spahni, R., Wania, R., Neef, L., van Weele, M., Pison, I., Bousquet, P., Frankenberg, C., Foster, P. N., Joos, F., Prentice, I. C., and van Velthoven, P.: Constraining global methane emissions and uptake by ecosystems, Biogeosciences, 8, 1643–1665, https://doi.org/10.5194/bg-8-1643-2011, 2011. 

Spahni, R., Joos, F., Stocker, B. D., Steinacher, M., and Yu, Z. C.: Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century, Clim. Past, 9, 1287–1308, https://doi.org/10.5194/cp-9-1287-2013, 2013. 

Stocker, B. D., Spahni, R., and Joos, F.: DYPTOP: a cost-efficient TOPMODEL implementation to simulate sub-grid spatio-temporal dynamics of global wetlands and peatlands, Geosci. Model Dev., 7, 3089–3110, https://doi.org/10.5194/gmd-7-3089-2014, 2014. 

Stohl, A., Forster, C., Frank, A., Seibert, P., and Wotawa, G.: Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2, Atmos. Chem. Phys., 5, 2461–2474, https://doi.org/10.5194/acp-5-2461-2005, 2005. 

Susiluoto, J., Raivonen, M., Backman, L., Laine, M., Makela, J., Peltola, O., Vesala, T., and Aalto, T.: Calibrating the sqHIMMELI v1.0 wetland methane emission model with hierarchical modeling and adaptive MCMC, Geosci. Model Dev., 11, 1199–1228, https://doi.org/10.5194/gmd-11-1199-2018, 2018. 

Thompson, R. L. and Stohl, A.: FLEXINVERT: an atmospheric Bayesian inversion framework for determining surface fluxes of trace species using an optimized grid, Geosci. Model Dev., 7, 2223–2242, https://doi.org/10.5194/gmd-7-2223-2014, 2014. 

Thompson, R. L., Nisbet, E. G., Pisso, I., Stohl, A., Blake, D., Dlugokencky, E. J., Helmig, D., and White, J. W. C.: Variability in Atmospheric Methane From Fossil Fuel and Microbial Sources Over the Last Three Decades, Geophys. Res. Lett., 112, 11499–11508, https://doi.org/10.1029/2018GL078127, 2018. 

Thompson, R. L., Lassaletta, L., Patra, P. K., Wilson, C., Wells, K. C., Gressent, A., Koffi, E. N., Chipperfield, M. P., Winiwarter, W., Davidson, E. A., Tian, H., and Canadell, J. G.: Acceleration of global N2O emissions seen from two decades of atmospheric inversion, Nat. Clim. Change, 9, 993–998, https://doi.org/10.1038/s41558-019-0613-7, 2019. 

Tian, H., Xu, X., Liu, M., Ren, W., Zhang, C., Chen, G., and Lu, C.: Spatial and temporal patterns of CH4 and N2O fluxes in terrestrial ecosystems of North America during 1979–2008: application of a global biogeochemistry model, Biogeosciences, 7, 2673–2694, https://doi.org/10.5194/bg-7-2673-2010, 2010. 

Tian, H., Chen, G., Lu, C., Xu, X., Ren, W., Zhang, B., Banger, K., Tao, B., Pan, S., Liu, M., Zhang, C., Bruhwiler, L., and Wofsy, S.: Global methane and nitrous oxide emissions from terrestrial ecosystems 3735 due to multiple environmental changes, Ecosystem Health and Sustainability, 1, 1–20, https://doi.org/10.1890/ehs14-0015.1, 2015. 

Tian, H., Yang, J., Xu, R., Lu, C., Canadell, J. G., Davidson, E. A., Jackson, R. B., Arneth, A., Chang, J., Ciais, P., Gerber, S., Ito, A., Joos, F., Lienert, S., Messina, P., Olin, S., Pan, S., Peng, C., Saikawa, E., Thompson, R. L., Vuichard, N., Winiwarter, W., Zaehle, S., and Zhang, B.: Global soil nitrous oxide emissions since the preindustrial era estimated by an ensemble of terrestrial biosphere models: Magnitude, attribution, and uncertainty, Glob. Change Biol., 25, 640–659, https://doi.org/10.1111/gcb.14514, 2019. Nature, 586, 248–256, https://doi.org/10.1038/s41586-020-2780-0, 2020. 

Tian, H., Xu, R., Canadell, J. G., Thompson, R. L., Winiwarter, W., Suntharalingam, P., Davidson, E. A., Ciais, P., Jackson, R. B., Janssens-Maenhout, G., Prather, M. J., Regnier, P., Pan, N., Pan, S., Peters, G. P., Shi, H., Tubiello, F. N., Zaehle, S., Zhou, F., Arneth, A., Battaglia, G., Berthet, S., Bopp, L., Bouwman, A. F., Buitenhuis, E. T., Chang, J., Chipperfield, M. P., Dangal, S. R. S., Dlugokencky, E., Elkins, J. W., Eyre, B. D., Fu, B., Hall, B. D., Ito, A., Joos, F., Krummel, P. B., Landolfi, A., Laruelle, G. G., Lauerwald, R., Li, W., Lienert, S., Maavara, T., Macleod, M., Millet, D. B., Olin, S., Patra, P. K., Prinn, R. G., Raymond, P. A., Ruiz, D. J., van der Werf, G. R., Vuichard, N., Wang, J., Weiss, R. F., Wells, K. C., Wilson, C., Yang, J., and Yao, Y.: A comprehensive quantification of global nitrous oxide sources and sinks, Nature, 586, 248–256, https://doi.org/10.1038/s41586-020-2780-0, 2020. 

Trusilova, K., Rödenbeck, C., Gerbig, C., and Heimann, M.: Technical Note: A new coupled system for global-to-regional downscaling of CO2 concentration estimation, Atmos. Chem. Phys., 10, 3205–3213, https://doi.org/10.5194/acp-10-3205-2010, 2010. 

Tsuruta, A., Aalto, T., Backman, L., Peters, W., Krol, M., van der Laan-Luijkx, I. T., Hatakka, J., Heikkinen, P., Dlugokencky, E. J., Spahni, R., and Paramonova, N.: evaluating atmospheric methane inversion model results for Pallas, northern Finland, Boreal Environ. Res., 20, 506–525, 2015. 

Tsuruta, A., Aalto, T., Backman, L., Hakkarainen, J., van der Laan-Luijkx, I. T., Krol, M. C., Spahni, R., Houweling, S., Laine, M., Dlugokencky, E., Gomez-Pelaez, A. J., van der Schoot, M., Langenfelds, R., Ellul, R., Arduini, J., Apadula, F., Gerbig, C., Feist, D. G., Kivi, R., Yoshida, Y., and Peters, W.: Global methane emission estimates for 2000–2012 from CarbonTracker Europe−CH4 v1.0, Geosci. Model Dev., 10, 1261–1289, https://doi.org/10.5194/gmd-10-1261-2017, 2017. 

Tubiello, F. N.: Greenhouse Gas Emissions Due to Agriculture, Elsevier Encyclopedia of Food Systems, 1, 196–205, https://doi.org/10.1016/B978-0-08-100596-5.21996-3, 2019. 

Tubiello, F. N., Salvatore, M., Rossi, S., Ferrara, A., Fitton, N., and Smith, P.: The FAOSTAT database of greenhouse gas emissions from agriculture, Environ. Res. Lett., 8, 015009, https://doi.org/10.1088/1748-326/8/1/015009, 2013. 

UNEP, United Nations Environment Programme: Emissions Gap Report 2019, available at: https://wedocs.unep.org/bitstream/handle/20.500.11822/30797/EGR2019.pdf?sequence=1&isAllowed=y (last access: September 2020), 2019. 

UNFCCC: Kyoto Climte Change Decision, available at: https://unfccc.int/process-and-meetings/conferences/past-conferences/kyoto-climate-change-conference-december-1997/decisions-kyoto-climate-change-conference-december-1997 (last access: October 2020), 1997. 

UNFCCC: Report of the Conference of the Parties on its nineteenth session, held in Warsaw from 11 to 23 November 2013: Addendum, Part two: Action taken by the Conference of the Parties at its nineteenth session, United Nations Framework Convention on Climate Change, FCCC/CP/2013/10/Add.3, available at: https://unfccc.int/resource/docs/2013/cop19/eng/10a03.pdf (last access: 4 May 2021), 2013 

UNFCCC: Decision 24/CP.19 Revision of the UNFCCC reporting guidelines on annual inventories for Parties included in Annex I to the Convention, FCCC/CP/2013/10/Add.3, 2014. 

UNFCCC: National Inventory Submissions 2019, available at: https://unfccc.int/process-and-meetings/transparency-and-reporting/reporting-and-review-under-the-convention/greenhouse-gas-inventories-annex-i-parties/national-inventory-submissions-2019 (last access: September 2020), 2019. 

UNFCCC NGHGI: NIR reports: UNFCCC: National Inventory Submissions 2018, available at: https://unfccc.int/process-and-meetings/transparencyand-reporting/reporting-and-review-under-theconvention/greenhouse-gas-inventories-annex-i-parties/nationalinventory-submissions-2018 (last access: January 2020), 2018. 

UNFCCC NGHGI: CRFs: available at: https://unfccc.int/process-and-meetings/transparency-and-reporting/reporting-and-review-under-the-convention/greenhouse-gas-inventories-annex-i-parties/national-inventory-submissions-2019 (last access: September 2020), 2019. 

Van Drecht, G., Bouwman, A., Harrison, J., and Knoop, J.: Global nitrogen and phosphate in urban wastewater for the period 1970 to 2050, Global Biogeochem. Cy., 23, GB0A03, https://doi.org/10.1029/2009GB003458, 2009. 

VERIFY: http://webportals.ipsl.jussieu.fr/VERIFY/FactSheets/, last access: January 2021. 

Vermeulen, A. T., Eisma, R., Hensen, A., and Slanina, J.: Transport model calculations of NW-European methane emissions, Environ. Sci. Policy, 2, 315–324, 1999. 

Vermeulen, A. T., Pieterse, G., Hensen, A., van den Bulk, W. C. M., and Erisman, J. W.: COMET: a Lagrangian transport model for greenhouse gas emission estimation – forward model technique and performance for methane, Atmos. Chem. Phys. Discuss., 6, 8727–8779, https://doi.org/10.5194/acpd-6-8727-2006, 2006. 

Wang, X., Jacob, D. J., Eastham, S. D., Sulprizio, M. P., Zhu, L., Chen, Q., Alexander, B., Sherwen, T., Evans, M. J., Lee, B. H., Haskins, J. D., Lopez-Hilfiker, F. D., Thornton, J. A., Huey, G. L., and Liao, H.: The role of chlorine in global tropospheric chemistry, Atmos. Chem. Phys., 19, 3981–4003, https://doi.org/10.5194/acp-19-3981-2019, 2019b. 

Watanabe, S., Miura, H., Sekiguchi, M., Nagashima, T., Sudo, K., Emori, S., and Kawamiya, M.: Development of an atmospheric general circulation model for integrated Earth system modeling on the Earth Simulator, J. Earth Simulator, 9, 27–35, 2008. 

Weiss, F. and Leip, A.: Greenhouse gas emissions from the EU livestock sector: A life cycle assessment carried out with the CAPRI model, Agriculture, Ecosyst. Environ., 149, 124–134, 2012. 

Wilson, C., Chipperfield, M. P., Gloor, M., and Chevallier, F.: Development of a variational flux inversion system (INVICAT v1.0) using the TOMCAT chemical transport model, Geosci. Model Dev., 7, 2485–2500, https://doi.org/10.5194/gmd-7-2485-2014, 2014. 

Winiwarter, W., Höglund-Isaksson, L., Klimont, Z., Schöpp, W., and Amann, M.: Technical opportunities to reduce global anthropogenic emissions of nitrous oxide, Environ. Res. Lett., 13, 014011, https://doi.org/10.1088/1748-9326/aa9ec9, 2018. 

WMO: United in Science Report, available at: https://public.wmo.int/en/resources/united_in_science (last access: January 2020), 2019. 

Yao, Y., Tian, H., Shi, H., Pan, S., Xu, R., Pan, N., and Canadell, J. G.: Increased global nitrous oxide emissions from streams and rivers in the Anthropocene, Nat. Clim. Change, 10, 138–142, https://doi.org/10.1038/s41558-019-0665-8, 2020. 

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

Zaehle, S., Ciais, P., Friend, A. D., and Prieur, V.: Carbon benefits of anthropogenic reactive nitrogen offset by nitrous oxide emissions, Nature Geosci., 4, 601–605, 2011. 

Zhang, Z., Zimmermann, N. E., Kaplan, J. O., and Poulter, B.: Modeling spatiotemporal dynamics of global wetlands: comprehensive evaluation of a new sub-grid TOPMODEL parameterization and uncertainties, Biogeosciences, 13, 1387–1408, https://doi.org/10.5194/bg-13-1387-2016, 2016b. 

Zheng, B., Chevallier, F., Ciais, P., Yin, Y., and Wang, Y.: On the role of the flaming to smoldering transition in the seasonal cycle of African fire emissions, Geophys. Res. Lett., 45, 11998–12007, https://doi.org/10.1029/2018GL079092, 2018a. 

Zheng, B., Chevallier, F., Ciais, P., Yin, Y., Deeter, M., Worden, H., Wang, Y. L., Zhang, Q., and He, K. B.: Rapid decline in carbon monoxide emissions and export from East Asia between years 2005 and 2016, Environ. Res. Lett., 13, 044007, https://doi.org/10.1088/1748-9326/aab2b3, 2018b.  

Zhu, Q., Liu, J., Peng, C., Chen, H., Fang, X., Jiang, H., Yang, G., Zhu, D., Wang, W., and Zhou, X.: Modelling methane emissions from natural wetlands by development and application of the TRIPLEX-GHG model, Geosci. Model Dev., 7, 981–999, https://doi.org/10.5194/gmd-7-981-2014, 2014. 

Zhu, Q., Peng, C., Chen, H., Fang, X., Liu, J., Jiang, H., Yang, Y., and Yang, G.: Estimating global natural wetland methane emissions using process modelling: spatio-temporal patterns and contributions to atmospheric methane fluctuations, Global Ecol. Biogeogr., 24, 959–972, 2015. 

Zhuang, Q., Melillo, J. M., Kicklighter, D. W., Prinn, R. G., McGuire, D. A., Steudler, P. A., Felzer, B. S., and Hu, S.: Methane fluxes between terrestrial ecosystems and the atmosphere at northern high latitudes during the past century: A retrospective analysis with a process-based biogeochemistry model, Global Biogeochem. Cy., 18, GB3010, https://doi.org/10.1029/2004GB002239, 2004. 

1

The rapid development of the gas industry in the 1980s in the former USSR.

2

Annex I Parties include the industrialized countries that were members of the OECD (Organization for Economic Co-operation and Development) in 1992 plus countries with economies in transition (the EIT Parties), including the Russian Federation, the Baltic states, and several central and eastern European states (UNFCCC, https://unfccc.int/parties-observers, last access: February 2020).

3

For most Annex I Parties, the historical base year is 1990. However, parties included in Annex I with an economy in transition during the early 1990s (EIT Parties) were allowed to choose 1 year up to a few years before 1990 as reference because of a non-representative collapse during the breakup of the Soviet Union (e.g., Bulgaria, 1988; Hungary, 1985–1987; Poland, 1988; Romania, 1989; and Slovenia, 1986).

4

The reason for using the median instead of mean for the ensembles is because there is a large spread between global inversions, and we do not want to be biased by outliers/extremes.

5

t−1 refers to an early estimate of the GHG emissions for the preceding year, as required by Regulation (EU) 525/2013 of the European Parliament and of the Council.

6

Denmark, Finland, Germany, Ireland, Latvia, Sweden, France, Estonia and Spain – in total these nine countries report in 2017 11.2 kt of CH4 from managed wetlands (UNFCCC 2019, CRF Table4(II)D: https://unfccc.int/documents/194946, last access: September 2020).

8

Auto-producers of electricity and heat: cogeneration by industries and companies for housing management (central heating and other services) (Olivier et al., 2017, PBL report).

9

The last revision was made by COP 19 in 2013 (UNFCCC, 2013)

Short summary
This study is topical and provides a state-of-the-art scientific overview of data availability from bottom-up and top-down CH4 and N2O emissions in the EU27 and UK. The data integrate recent emission inventories with process-based model data and regional/global inversions for the European domain, aiming at reconciling them with official country-level UNFCCC national GHG inventories in support to policy and to facilitate real-time verification procedures.
Altmetrics
Final-revised paper
Preprint