# Gridded fossil CO2 emissions and related O2 combustion consistent with national inventories 1959–2018

Matthew William Jones, Robbie M. Andrew, Glen P. Peters, Greet G.A. Janssens-Maenhout, Anthony J. De-Gol, Philippe Ciais, Prabir K. Patra,

Frédéric Chevallier, Corinne Corinne Le quéré

► To cite this version:

Matthew William Jones, Robbie M. Andrew, Glen P. Peters, Greet G.A. Janssens-Maenhout, Anthony J. De-Gol, et al.. Gridded fossil CO2 emissions and related O2 combustion consistent with national inventories 1959–2018. Scientific Data, 2021, 8 (1), 10.1038/s41597-020-00779-6. doi:10.1038/s41597-020-00779-6

HAL Id: hal-03487836

https://hal.science/hal-03487836

Submitted on 19 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# Gridded fossil CO₂ emissions and related O₂ combustion consistent with national inventories 1959–2018

Matthew W. Jones<fnref n="1" />✉, Robbie M. Andrew², Glen P. Peters², Greet Janssens-Maenhout³, Anthony J. De-Gol<fnref n="1" />, Philippe Ciais⁴, Prabir K. Patra⁵, Frederic Chevallier⁴ & Corinne Le Quéré<fnref n="1" />

Quantification of CO₂ fluxes at the Earth's surface is required to evaluate the causes and drivers of observed increases in atmospheric CO₂ concentrations. Atmospheric inversion models disaggregate observed variations in atmospheric CO₂ concentration to variability in CO₂ emissions and sinks. They require prior constraints fossil CO₂ emissions. Here we describe GCP-GridFED (version 2019.1), a gridded fossil emissions dataset that is consistent with the national CO₂ emissions reported by the Global Carbon Project (GCP). GCP-GridFEDv2019.1 provides monthly fossil CO₂ emissions estimates for the period 1959–2018 at a spatial resolution of 0.1°. Estimates are provided separately for oil, coal and natural gas, for mixed international bunker fuels, and for the calcination of limestone during cement production. GCP-GridFED also includes gridded estimates of O₂ uptake based on oxidative ratios for oil, coal and natural gas. It will be updated annually and made available for atmospheric inversions contributing to GCP global carbon budget assessments, thus aligning the prior constraints on top-down fossil CO₂ emissions with the bottom-up estimates compiled by the GCP.

## Background & Summary

Fossil fuel use, cement production and land-use change have perturbed the natural carbon cycle and increased the concentration of carbon dioxide (CO₂) in the Earth's atmosphere by almost 50% since 1750, from 277 ppm in 1750 to 407 ppm in 2018¹⁻³. Routine assessment of the global carbon cycle is required to monitor the ongoing increases in atmospheric CO₂ concentrations, evaluate the causes and drivers of this trend, and quantify the impact of policies that aim to stabilise and reverse it³⁻⁵. The global carbon budget (GCB) was evaluated on multi-year time scales by each of the foregoing Intergovernmental Panel on Climate Change (IPCC) assessment reports⁶⁻¹⁰, while the Global Carbon Project (GCP) has published an annual assessment of the GCB on an annual time-scale for over a decade³,¹¹⁻¹³.

The GCP disaggregates the annual GCB into six components: atmospheric growth (G<sub>ATM</sub>); CO₂ emissions due to fossil fuel combustion, non-combustion uses of fossil fuels, and cement production (E<sub>FF</sub>); CO₂ emissions due to land-use change (E<sub>LUC</sub>); uptake of CO₂ by the global ocean (S<sub>OCEAN</sub>); uptake of CO₂ by the terrestrial biosphere (S<sub>LAND</sub>), the two later fluxes from ocean and land carbon models, respectively; and a budget imbalance term (B<sub>IM</sub>). G<sub>ATM</sub> is the most precisely constrained term of the budget (1σ of 4%)³, while E<sub>FF</sub>, E<sub>LUC</sub>, S<sub>OCEAN</sub> and S<sub>LAND</sub> rely on analysis of national emissions reports¹⁴⁻¹⁶, satellite observations¹⁷,¹⁸, and process-based models³,¹⁹,²⁰ and are more uncertain. If E<sub>FF</sub>, E<sub>LUC</sub>, S<sub>OCEAN</sub> and S<sub>LAND</sub> were perfectly constrained then their sum would be equal to the measured change in the atmospheric stock of CO₂ (G<sub>ATM</sub>). However, the independent analysis of G<sub>ATM</sub>, E<sub>FF</sub>, E<sub>LUC</sub>, S<sub>OCEAN</sub> and S<sub>LAND</sub> using different methodologies results in an unconstrained budget, and a small budget

<footnote n="1">¹Tyndall Centre for Climate Change Research, School of Environmental Sciences, University of East Anglia, Norwich, NR4 7TJ, UK. ²CICERO Center for International Climate Research, Oslo, 0349, Norway. ³European Commission, Joint Research Centre (JRC), Via E. Fermi 2749 (T.P. 123), 21027, Ispra, Varese, Italy. ⁴Laboratoire des Sciences du Climat et de l'Environnement, Institut Pierre-Simon Laplace, CEA-CNRS-UVSQ, CE Orme des Merisiers, 91191, Gif sur Yvette, France. ⁵Research Institute for Global Change, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Yokohama, 236 0001, Japan. ✉e-mail: matthew.w.jones@uea.ac.uk</footnote>
imbalance term (B<sub>IM</sub>) is required to close the budget. The global carbon budget is thus closed as follows (S<sub>OCEAN</sub> and S<sub>LAND</sub> hold negative values)<sup>3</sup>:

$$G_{ATM} = E_{FF} + E_{LUC} + S_{OCEAN} + S_{LAND} + B_{IM} \quad (1)$$

Inversion models use an integrated approach to simultaneously quantify all fluxes of the global carbon budget and they are, by design, constrained by observations of atmospheric CO<sub>2</sub> mole fraction, or satellite derived products of column CO<sub>2</sub>. Inversion models prescribe the fossil carbon emissions (E<sub>FF</sub>) because the current density of the surface network and the sampling of the atmosphere by satellites is too sparse to quantify this flux separately, and then estimate the total land flux (F<sub>LAND</sub> = S<sub>LAND</sub> + E<sub>LUC</sub>) and the ocean sink (S<sub>OCEAN</sub>) using a modelling framework that minimises data-model mismatch across all fluxes according to a cost function (see examples in refs. 21–32 and studies cited therein). By synchronously quantifying E<sub>FF</sub>, F<sub>LAND</sub> and S<sub>OCEAN</sub>, inversion models avoid budget imbalance and hence the global carbon budget equation is closed without a B<sub>IM</sub> term as follows:

$$G_{ATM} = E_{FF}^{inv} + F_{LAND}^{inv} + S_{OCEAN}^{inv} \quad (2)$$

Inversion models require prior constraints on the regional distribution of the CO<sub>2</sub> fluxes that they seek to disaggregate. Here we describe our development of the Global Carbon Budget Gridded Fossil Emissions Dataset (GCP-GridFED; version 2019.1), a new gridded 0.1° × 0.1° global dataset of monthly CO<sub>2</sub> emissions resulting from fossil fuel oxidation and the calcination of limestone during cement production. The gridded nation- and source- specific emissions in GCP-GridFED are consistent with the nation- and source- specific emissions inventories compiled for the GCP’s 2019 GCB assessment<sup>3,33</sup> and for version 2019.1 cover the period 1959-2018. The GCP-GridFED will be updated each year for use by inversion models contributing to the annual updates of the GCB, thus aligning the prior constraints on top-down estimates of fossil CO<sub>2</sub> emissions with the bottom-up estimates used by the GCP.

Gridded estimates of uncertainty in CO<sub>2</sub> emissions are provided as an additional layer of GCP-GridFED and are based on the relative uncertainties (1σ) in fossil CO<sub>2</sub> presented in the uncertainty assessment of the GCB<sup>3</sup> and the relative uncertainties amongst emission sectors<sup>34</sup>. Uncertainties associated with the spatial disaggregation of national emissions are not included (see ‘CO<sub>2</sub> Emissions Uncertainty’). Our approach to uncertainty quantification is broadly representative of the sectoral contributions to total emissions in each grid cell, which changes throughout the time series, and of differences in uncertainty across national emission reports. Inversion models may utilise these uncertainty grids but with the freedom to build more complex covariance structures to suit their requirements.

The global cycles of carbon and oxygen are coupled through their dual involvement in carboxylation reactions (photosynthesis), which consume CO<sub>2</sub> and emit O<sub>2</sub>, and oxidation reactions (respiration and combustion), which consume O<sub>2</sub> and emit CO<sub>2</sub> (refs. 35,36). In addition to CO<sub>2</sub> alone, some inversion models are able to constrain surface fluxes of O<sub>2</sub> or atmospheric potential oxygen (APO ≈ O<sub>2</sub> + 1.1CO<sub>2</sub>)<sup>37–40</sup>. Such models can utilise dual atmospheric measurements of CO<sub>2</sub> and O<sub>2</sub> and dual priors for CO<sub>2</sub> and O<sub>2</sub> surface fluxes and synchronously minimise data-model mismatch with respect to CO<sub>2</sub> and O<sub>2</sub>. Alternatively, O<sub>2</sub> fluxes can be constrained independently using atmospheric O<sub>2</sub> observations and O<sub>2</sub> surface flux priors. GCP-GridFED includes dual estimates of atmospheric O<sub>2</sub> uptake due to the oxidation of fossil fuels, with the aim of supporting the inverse modelling of O<sub>2</sub> or APO and with the view that the data can be used in multi-decadal analyses of the global oxygen budget. Our O<sub>2</sub> uptake estimates are based on the oxidative ratios (OR; uptake of O<sub>2</sub>/emission of CO<sub>2</sub>)<sup>36</sup> applied to the CO<sub>2</sub> emission estimates for coal, oil, natural gas oxidation<sup>36</sup>.

## Methods

**Overview.** GCP-GridFED was produced by scaling monthly gridded emissions for the year 2010, from the Emissions Database for Global Atmospheric Research (EDGAR; version 4.3.2)<sup>41</sup>, to the national annual emissions estimates compiled as part of the 2019 global carbon budget (GCB-NAE) for the years 1959–2018 (ref. 3). EDGAR data for the year 2010 is used because monthly gridded data was only available for this year at the time of product development (new data for 2015 was published recently and will be adopted in future versions of GCP-GridFED)<sup>42</sup>. We describe the key features of the EDGAR and GCB-NAE datasets below (see ‘Input Datasets’).

GCB-NAE and EDGAR provide information regarding the global emission of CO<sub>2</sub> through the combustion of fossil fuels, industrial processes and cement production, and some other minor sources (e.g. consumption of lubricants and paraffin waxes, solvent use, agricultural liming); nonetheless, their merits differ. GCB-NAE provides a consistent long-term dataset of annual national CO<sub>2</sub> emissions (1750–2018), however this dataset is not spatially-explicit below the country level and does not include sub-annual variability in CO<sub>2</sub> emissions. EDGAR provides estimates at high spatial resolution for specific fuels and sectors with a representation of the monthly distribution of emissions. However the EDGARv4.3.2 estimates are only available for 1970–2010 and a constant monthly distribution, matching the year 2010, is used throughout the time series<sup>41</sup>. Our approach merged these two complementary datasets to create a long-term (1959–2018) and gridded (0.1° × 0.1°) dataset of global monthly CO<sub>2</sub> emissions. The start year of 1959 aligns with the period of direct atmospheric measurements of CO<sub>2</sub> concentration<sup>43</sup>. Our approach is to scale EDGAR’s 2010 monthly gridded CO<sub>2</sub> emissions to match the annual gridded CO<sub>2</sub> emissions from GCB-NAE on a nation- and fuel- specific basis (see ‘emissions scaling protocol’, Fig. 1; Table 1).

GCP-GridFED includes additional data layers that are beneficial to inversion models. Gridded uncertainty in CO<sub>2</sub> emissions from each nation and emissions sector is also propagated to our nation-, year- and fuel- specific emissions estimates (Table 2). Gridded estimates of the uptake of O<sub>2</sub> related to oil, coal and natural gas use are also made using the literature-based oxidative ratios presented in the CO<sub>2</sub> release and Oxygen uptake from Fossil Fuel Emission Estimate (COFFEE) dataset<sup>36</sup>.
<figure><img src="image_3.png" /><figcaption>**Fig. 1** A conceptual depiction of the emissions scaling protocol used to produce GCP-GridFED as described in section 2.2. Descriptions of the input datasets are provided in section 2.1. This figure does not depict the procedure used to calculate emissions uncertainty or O₂ combustion; however, the figure indicates the stages at which these additional outputs are produced within the CO₂ emissions scaling protocol (marked as '*' and '#', respectively). Uncertainties in CO₂ concentration are calculated for the EDGAR dataset in the year 2010 and scaled using the same factors as the central estimates for CO₂ emission. O₂ combustion estimates are calculated using oxidative ratios applied to the CO₂ emissions from fuel combustion, prior to the final step of the protocol. Full details of the procedure used to produce gridded CO₂ uncertainties and O₂ combustion estimates can be found in sections 2.3 and 2.4, respectively.</figcaption></figure>

We provide Figs. 2–12, the descriptive details in Tables 3 and 4, the summary statistics in Tables 5–7 and Online-Only Table 1 to outline the key features of GCP-GridFEDv2019.1 and assist with its technical validation.

**Input datasets.** National annual emissions from the global carbon budget 2019 (GCB-NAE). The GCB estimates national annual emissions of CO₂ due to coal, oil and natural gas combustion, the oxidative use of these fuels in non-combustive industrial processes, and the production of cement clinker<sup>3,14–17,44</sup>. National CO₂ emissions are preferentially taken from the country submissions to the United Nations Framework Convention on Climate Change (UNFCCC) for 42 "Annex I" countries over the period 1990–2018<sup>44</sup>. These countries were members of the Organisation for Economic Co-operation and Development (OECD) in 1992, plus 16 non-OECD European countries and Russia, and contributed ~60% of total global emissions in 1990. Emissions in other countries and in Annex I countries prior to 1990 derive from the Carbon Dioxide Information Analysis Center (CDIAC)¹⁵ and are rooted in energy statistics published by the United Nations (UN)¹⁶,⁴⁵. For recent years not covered by either the UNFCCC or CDIAC datasets, the national emissions are predicted using national or regional energy growth rates from the annual BP Statistical Review of World Energy¹⁴. National cement emissions are based on national inventories of cement production and ratios of clinker production from officially reported clinker production data and emission factors, IPCC default emission factors, industry-reported clinker production, and survey-based clinker ratios¹⁶.

*Gridded monthly emissions from EDGAR.* The Emissions Database for Global Atmospheric Research (EDGAR) version 4.3.2<sup>41</sup> is a dataset of global emissions of gases and particulates, including CO₂, based on available national statistics, default emission factors and methods recommended by IPCC<sup>46,47</sup>. EDGAR uses a bottom-up approach that calculates gridded (0.1° × 0.1°) monthly CO₂ emissions for activity sectors based on: statistics that track national levels of each activity; proxy data representing the spatial and temporal distribution of each activity; the mix of technologies used to perform each activity; the fuel mix used by each technology, and; emissions factors for the technology and fuel combinations, which are also corrected for the emission control technologies in place.
<table><thead><tr><th rowspan="2"></th><th colspan="5">GCP-GridFED Source Classes</th></tr><tr><th>Oil</th><th>Natural gas</th><th>Coal</th><th>Bunker Oil</th><th>Cement</th></tr></thead><tbody><tr><td>Power Industry</td><td>1A1a_OIL</td><td>1A1a_GAS</td><td>1A1a_COAL</td><td></td><td></td></tr><tr><td>Fuel Transformation Industries</td><td>1A1b_1A1c_1A5b1_1B1b_1B2a5_1B2a6_1B2b5_2C1b_OIL</td><td>1A1b_1A1c_1A5b1_1B1b_1B2a5_1B2a6_1B2b5_2C1b_GAS</td><td>1A1b_1A1c_1A5b1_1B1b_1B2a5_1B2a6_1B2b5_2C1b_COAL</td><td></td><td></td></tr><tr><td>Manufacturing</td><td>1A2_OIL</td><td>1A2_GAS</td><td>1A2_COAL</td><td></td><td></td></tr><tr><td>Buildings</td><td>1A4_OIL</td><td>1A4_GAS</td><td>1A4_COAL</td><td></td><td></td></tr><tr><td>Transport: road</td><td>1A3b_OIL</td><td></td><td></td><td></td><td></td></tr><tr><td>Transport: Rail, Pipelines, Off-Road</td><td>1A3c_1A3e (a)</td><td></td><td></td><td></td><td></td></tr><tr><td>Fuel Exploitation</td><td>1B1a_1B2a1_1B2a2_1B2a3_1B2a4_1B2c_OIL</td><td>1B1a_1B2a1_1B2a2_1B2a3_1B2a4_1B2c_GAS</td><td></td><td></td><td></td></tr><tr><td>Production of Iron and Steel</td><td></td><td></td><td>2C1a_2C1c_2C1d_2C1e_2C1f_2C2 (b)</td><td></td><td></td></tr><tr><td>Production of Non-ferrous Metals</td><td></td><td></td><td>2C3_2C4_2C5 (b)</td><td></td><td></td></tr><tr><td>Fossil Fuel Fires</td><td></td><td></td><td>7 A (c)</td><td></td><td></td></tr><tr><td>Aviation</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>(Cruising)</td><td>1A3a_CRS_home</td><td></td><td></td><td>1A3a_CRS_international</td><td></td></tr><tr><td>Aviation (Landing and Take-Off)</td><td>1A3a_LTO_home</td><td></td><td></td><td>1A3a_LTO_international</td><td></td></tr><tr><td>Aviation (Climbing & Descent)</td><td>1A3a_CDS_home</td><td></td><td></td><td>1A3a_CDS_international</td><td></td></tr><tr><td>Shipping</td><td></td><td></td><td></td><td>1A3d_1C2 (d)</td><td></td></tr><tr><td>Non-energy Use of Fuels</td><td>2 G (e)</td><td></td><td></td><td></td><td></td></tr><tr><td>Solvents and Product Use</td><td>3 (f)</td><td></td><td></td><td></td><td></td></tr><tr><td>Chemical Processes</td><td></td><td>2B (g)</td><td></td><td></td><td></td></tr><tr><td>Non-metallic Minerals Production</td><td></td><td></td><td></td><td></td><td>2 A (h)</td></tr></tbody></table>

**Table 1.** The relation of GCP-GridFED source classes to EDGAR activity sectors. The filename of the EDGAR grid layers is shown for each EDGAR activity sector. For a full description of each EDGAR sector, see Janssens-Maenhout *et al.* (ref. 41). Assumptions made as to the fuels contributing to unstratified EDGAR sectors are as follows (see section 2.1.2 for further detail): (a) all emissions from off-road, rail and pipeline transport relate to the combustion of oil. (b) all emissions from the production of steel, iron and non-ferrous metals relate to the combustion of coal. (c) all emissions from fossil fuel fires relate to underground coal fires. (d) all emissions from shipping were assumed to relate to bunker oil combustion. (e) all emissions from the non-energy use of fuels sector relate to non-combustion use of oil. (f) all emissions from the solvents and product use sector relate to non-combustion use of oil. (g) all chemical process emissions relate to the non-combustion use of natural gas. (h) all non-metallic minerals production from EDGAR are assumed to be emissions from cement (clinker) production.

A detailed description of EDGAR's gridding procedure is available elsewhere (refs. 41,48) however we summarise below the key features of its design:

* 28 EDGAR activity sectors are based on the 48 sectors defined by IPCC guidelines46,47.

* Activity in each sector is tracked from 1970–2015 using statistics that represent demand and supply of goods and energy, including: fuel-specific energy balances, fuel production, commodity production and cement clinker production and agriculture-related activities.

* Emission factors are taken from the guidelines issued by the IPCC46,47 and are assigned to each country in the following order of preference: national, regional, country group (Annex I/non-Annex I).

* National emissions of CO₂ from each sector are distributed across months using sector-specific or, preferentially, technology-specific monthly shares.

* Emissions are distributed in space using spatial proxy data (that vary stepwise over time 1990–2010), such as population density, point source locations and transport routes.

Some of the uncertainties associated with using proxy data to disaggregate emission in time and space are considered in later sections (see 'CO₂ emissions uncertainty').

**EDGAR sectors included in GCP-GridFED.** Of the 28 EDGAR sectors, the 18 relating to fossil fuel combustion, non-combustion use of fossil fuel and cement production were used in GCP-GridFED. These 18 sectors were selected to correspond as closely as possible with the activities included in the GCB-NAE emission estimates. The 18 activity sectors incorporated from EDGAR into GCP-GridFED are shown in Table 1.

Where possible, emissions from each EDGAR sector were further separated into specific fuels using fuel-specific data from an intermediate processing step of the EDGAR gridding protocol41. Where this was not possible, it was necessary to make the assumptions that follow about the fuels that contribute to emissions in each
<table><caption>**Table 2.** Calculation of uncertainties for each sector in GCP-GridFEDv2019.1. Calculations are based on (i) the ratio of the uncertainty in emissions from each sector (U_TNOs) to the uncertainty in total emissions (U_TNO Tot) and (ii) the relative uncertainty in emissions for Annex I countries (5%) and other countries (10%) from GCB-NAE³. (a) GCP-GridFED source class codes are adopted from EDGAR. See Table 1 for more information. (b) 5% uncertainty in total emissions. (c) 10% uncertainty in total emissions.</caption><thead><tr><th>TNO Sector</th><th>GCP-GridFED source classes included (a)</th><th>U_TNO: Relative Uncertainty Estimate</th><th>U_TNO: Sector-level Uncertainty ÷ Total Uncertainty</th><th>U_GridFED: Annex I Countries (b)</th><th>U_GridFED: Other Countries (c)</th></tr></thead><tbody><tr><td rowspan="3">Public power</td><td>1A1a_OIL</td><td>86.4%</td><td>7.6</td><td>38%</td><td>76%</td></tr><tr><td>1A1a_NATURAL GAS</td><td>86.4%</td><td>7.6</td><td>38%</td><td>76%</td></tr><tr><td>1A1a_COAL</td><td>86.4%</td><td>7.6</td><td>38%</td><td>76%</td></tr><tr><td rowspan="10">Industry</td><td>1A2_OIL</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>1A2_NATURAL GAS</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>1A2_COAL</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>2C1a_2C1c_2C1d_2C1e_2C1f_2C2</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>2C3_2C4_2C5</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>2B</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>1A1b_1A1c_1A5b1_1B1b_1B2a5_1B2a6_1B2b5_2C1b_OIL</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>1A1b_1A1c_1A5b1_1B1b_1B2a5_1B2a6_1B2b5_2C1b_NATURAL GAS</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td>1A1b_1A1c_1A5b1_1B1b_1B2a5_1B2a6_1B2b5_2C1b_COAL</td><td>18.2%</td><td>1.6</td><td>8%</td><td>16%</td></tr><tr><td rowspan="3">Other stationary combustion</td><td>1A4_OIL</td><td>15.4%</td><td>1.4</td><td>7%</td><td>14%</td></tr><tr><td>1A4_NATURAL GAS</td><td>15.4%</td><td>1.4</td><td>7%</td><td>14%</td></tr><tr><td>1A4_COAL</td><td>15.4%</td><td>1.4</td><td>7%</td><td>14%</td></tr><tr><td rowspan="3">Fugitive</td><td>1B1a_1B2a1_1B2a2_1B2a3_1B2a4_1B2c_OIL</td><td>34.1%</td><td>3.0</td><td>15%</td><td>30%</td></tr><tr><td>1B1a_1B2a1_1B2a2_1B2a3_1B2a4_1B2c_NATURAL GAS</td><td>34.1%</td><td>3.0</td><td>15%</td><td>30%</td></tr><tr><td>7A</td><td>34.1%</td><td>3.0</td><td>15%</td><td>30%</td></tr><tr><td rowspan="2">Solvents</td><td>2G</td><td>50%</td><td>4.4</td><td>22%</td><td>44%</td></tr><tr><td>3</td><td>50%</td><td>4.4</td><td>22%</td><td>44%</td></tr><tr><td>Road transport</td><td>1A3b_OIL</td><td>13.9%</td><td>1.2</td><td>6%</td><td>12%</td></tr><tr><td rowspan="2">Shipping</td><td>1A3d_1C2</td><td>7.9%</td><td>0.7</td><td>3%</td><td>7%</td></tr><tr><td>1A3a_CDS</td><td>9.8%</td><td>0.9</td><td>4%</td><td>9%</td></tr><tr><td rowspan="2">Aviation</td><td>1A3a_CRS</td><td>9.8%</td><td>0.9</td><td>4%</td><td>9%</td></tr><tr><td>1A3a_LTO</td><td>9.8%</td><td>0.9</td><td>4%</td><td>9%</td></tr><tr><td>Off-road transport</td><td>1A3c_1A3e</td><td>22.4%</td><td>2.0</td><td>10%</td><td>20%</td></tr><tr><td>Cement (not included in the TNO dataset)</td><td>2A</td><td></td><td></td><td>5%</td><td>10%</td></tr><tr><td>Total Uncertainty</td><td></td><td>11.3%</td><td></td><td>5%</td><td>10%</td></tr></tbody></table>

sector. These assumptions are based on the sector descriptions provided in the IPCC guidelines⁴⁶,⁴⁷ and the major contributing activities and fuel dependencies in each sector. Specifically, we assume that:

* All chemical process emissions relate to the non-combustion use of natural gas.

* All emissions from the non-energy use of fuels sector relate to non-combustion use of oil. This sector chiefly comprises the use of waxes and lubricants.

* All emissions from the solvents and product use sector relate to non-combustion use of oil. This sector chiefly comprises solvents in paint, degreasing and dry cleaning, chemical products and other product use.

* All emissions from the production of steel, iron and non-ferrous metals relate to the oxidation of coal and production of cokes.

* All emissions from fossil fuel fires relate to underground coal fires. This sector also includes oil flaring emissions in Kuwait, however fossil fuel fire emissions were found to be negligible in Kuwait.

* All emissions from off-road, rail and pipeline transport relate to the combustion of oil.

* All emissions from the production of non-metallic minerals relate to cement clinker production.

National CO₂ emissions data were extracted from the EDGAR datasets for the purpose of national annual emissions scaling. National masks were based on the 'countries 2016' dataset of the Geographic Information System of the European Commission (EU-GISCO)⁴⁹.

The appropriate positioning of power plants is key to distributing total emissions accurately because the power sector accounts for ~45% of global emissions⁵⁰. Changes in the available datasets of power plant geolocations are common, and hence we note the importance of recording which datasets are used in each release of gridded
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Time series of (left column) annual CO₂ emissions (Gt CO₂ year⁻¹) and (right column) monthly fossil CO₂ emissions (Mt CO₂ day⁻¹) as estimated by GCP-GridFED. Uncertainties in CO₂ emissions are treated as 5% for Annex I nations and global, following the GCB uncertainty assessment³. (Top row) Total global emissions are disaggregated to (other rows) the top 4 emission regions. Crosses mark input data directly from GCB-NAE³,³³.</figcaption></figure>

emissions products. GridFEDv2019.1 adopts point source geolocations from EDGAR v4.3.2, which are scaled as described below (see ‘GCP-GridFED Protocol’). The EDGAR protocol for geolocating power plant emissions is summarised as follows, with full documentation provided by Janssens-Maenhaut et al.⁴¹. The location, fuel type and seasonality of power plant emissions derives from the CARMAv3.0 dataset⁵¹. The 2010 gridded emissions dataset used here as the scaling basis includes over 60,000 plants mapped globally in CARMAv3.0 in the year 2007. Standard QA/QC screening was applied to the CARMAv3.0 dataset, including gap-filling of missing (0,
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Time series of annual CO₂ emissions as estimated by GCP-GridFEDv2019.1 for the period 1959–2018 (Gt CO₂ year⁻¹). Uncertainties in CO₂ emissions are treated as 5% for Annex I nations and global, following the GCB uncertainty assessment³. (Top row, Left column) Total global emissions are disaggregated to the (other columns) top 4 emission regions and (other rows) source classes. Crosses mark input data directly from GCB-NAE³,³³.</figcaption></figure>

0) plant coordinates, correcting inverted (lon, lat) coordinates and adding some additional points for Russia. National power sector emissions for each fuel type are distributed across plants in proportion to their reported capacities. For larger countries (e.g. USA) with a non-uniform distribution of coal power plants, the fuel-specific distribution of emissions is considered a significant improvement over foregoing approaches. Emissions from each power plant reflect the fuel mix of the plant and the respective carbon intensity of emissions from that fuel mix. However, details of the technologies used by each plant, including carbon capture and storage, are not available. Alternative mappings of point sources can be based on night light detections by satellite⁵² or population data⁵³ but these are least aligned with EDGAR's 'bottom up' approach⁴¹.

*Heating and cooling degree day (HCDD) Correction.* The monthly distribution (seasonality) of global CO₂ emissions is principally determined by seasonality of climate in the Northern Hemisphere, and thus a peak in emissions occurs in the boreal winter months and a trough occurs in the boreal summer months. Although this seasonality is predictable, inter-annual variability in weather influences the distribution of emission across the months. Because the monthly emissions distribution in the EDGAR dataset is derived only from 2010 data, we applied a correction to the EDGAR data to account for the impacts of inter-annual variability on emissions. Specifically, we used a heating and cooling degree day (HCDD) correction to implement inter-annual variability
<table><caption>**Table 3.** Table of groups, variables, dimensions and units of the output files with naming convention GCP\_Global\_\{YYYY\}.nc. Numbers show the length of each dimension for each variable.</caption><thead><tr><th rowspan="2">Groups</th><th rowspan="2">Variables</th><th rowspan="2">Description</th><th rowspan="2">Unit</th><th colspan="3">Dimensions</th></tr><tr><th>Time</th><th>Degrees North of the equator (cell centres)</th><th>Degrees East of the Prime Meridian (cell centres)</th></tr></thead><tbody><tr><td rowspan="5">CO<sub>2</sub></td><td>COAL</td><td>Monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from COAL</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from OIL</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from GAS</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>CEMENT</td><td>Monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from CEMENT</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from BUNKER</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td rowspan="5">CO<sub>2</sub>_uncertainty</td><td>COAL</td><td>Uncertainty in monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from COAL</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Uncertainty in monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from OIL</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Uncertainty in monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from GAS</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>CEMENT</td><td>Uncertainty in monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from CEMENT</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Uncertainty in monthly emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> month-1) from BUNKER</td><td>kg CO<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td rowspan="4">O<sub>2</sub></td><td>COAL</td><td>Monthly uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> month-1) due to COAL use</td><td>kg O<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Monthly uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> month-1) due to OIL use</td><td>kg O<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Monthly uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> month-1) due to GAS use</td><td>kg O<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Monthly uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> month-1) due to BUNKER use</td><td>kg O<sub>2</sub> month-1</td><td>12</td><td>1,800</td><td>3,600</td></tr></tbody></table>

in the monthly distribution of CO₂ emissions from selected EDGAR sectors (power industry, 1A1a; buildings, 1A4; manufacturing, 1A2; and road transport, 1A3b; see Table 1). The HCDD correction approach was implemented as follows.

First, monthly (m) HCDDs were calculated based on gridded (0.5° × 0.5°) daily mean temperature (T) data for the years 1959–2018 from the Climatic Research Unit time-series version 4.03 (CRU-TSv4.03)⁵⁴ and following Spinoni *et al.* (refs. ⁵⁵,⁵⁶). For each 0.5° × 0.5° cell (i_r, j_r) of CRU-TSv4.03, HCDD was calculated as the absolute difference between the daily mean temperature of each month and an upper temperature threshold of 22 °C or a lower temperature threshold of 15.5 °C degrees⁵⁵,⁵⁶, multiplied by the number of days in the month (d).

$$HCDD_{m,i_r,j_r} = \begin{cases} (15.5 - \bar{T}_{m,i_r,j_r}) \cdot d, & \bar{T}_{m,i_r,j_r} < 15.5 \\ \bar{T}_{m,i_r,j_r} - 22.0 \cdot d, & \bar{T}_{m,i_r,j_r} > 22.0 \end{cases} \quad (3)$$

Second, the monthly fraction of annual HCDDs (HCDDfrac) was calculated for each month in the year 2010 as follows.

$$HCDDfrac_{m,i_r,j_r} = \frac{HCDD_{m,i_r,j_r}}{\sum HCDD_{i_r,j_r}} \quad (4)$$
<figure><img src="image_6.png" /><figcaption>Fig. 4 Time series of monthly CO₂ emissions as estimated by GCP-GridFEDv2019.1. for the period 1990–2018 (Mt CO₂ day⁻¹). Uncertainties in CO₂ emissions are treated as 5% for Annex I nations and global, following the GCB uncertainty assessment³. (Top row, Left column) Total global emissions are disaggregated to the (other columns) top 4 emission regions and (other rows) source classes.</figcaption></figure>

Third, the monthly fraction of annual emissions (*Efrac*) for each of the relevant sectors (s; 1A1a, 1A4, 1A2, 1A3b) was calculated on a reduced-resolution grid (0.5° × 0.5°; i_r, j_r) for each month as follows.

$$Efrac_{s,m,i_r,j_r} = \frac{EDGAR_{s,m,i_r,j_r}}{\sum EDGAR_{s,i_r,j_r}} \quad (5)$$

Fourth, a simple linear regression equation of the form below was fitted between monthly HCDDfrac and Efrac in the year 2010.

$$Efrac_{s,m,i_r,j_r} = a_{s,i_r,j_r} + (b_{s,i_r,j_r} \cdot HCDDfrac_{m,i_r,j_r}) + error_{s,i_r,j_r} \quad (6)$$

The HCDD correction was implemented in each year of the time series by predicting *Efrac* based on *HCDDfrac*. The correction was applied only to cells where the R² value of the linear regression equation exceeded 0.66 (where 66% of variation in *Efrac* was explained by variation in *HCDDfrac*). The following conditional approach was applied in all years (1959–2018).
<table><caption>**Table 4.** Table of groups, variables, dimensions and units of the output file GCPalmGlobalAnnual.nc. Numbers show the length of each dimension for each variable.</caption><thead><tr><th rowspan="2">Groups</th><th rowspan="2">Variables</th><th rowspan="2">Description</th><th rowspan="2">Unit</th><th colspan="3">Dimensions</th></tr><tr><th>Time</th><th>Latitude</th><th>Longitude</th></tr><tr><td></td><td></td><td></td><td></td><td>days since 1959-01-01</td><td>Degrees North of the equator (cell centres)</td><td>Degrees East of the Prime Meridian (cell centres)</td></tr></thead><tbody><tr><td rowspan="5">CO<sub>2</sub></td><td>COAL</td><td>Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from COAL</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from OIL</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from GAS</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>CEMENT</td><td>Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from CEMENT</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from BUNKER</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td rowspan="5">CO<sub>2</sub>_uncertainty</td><td>COAL</td><td>Uncertainty in Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from COAL</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Uncertainty in Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from OIL</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Uncertainty in Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from GAS</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>CEMENT</td><td>Uncertainty in Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from CEMENT</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Uncertainty in Annual emissions of CO<sub>2</sub> in each cell (kg CO<sub>2</sub> year-1) from BUNKER</td><td>kg CO<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td rowspan="4">O<sub>2</sub></td><td>COAL</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to COAL use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to OIL use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to GAS use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to BUNKER use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td rowspan="4">O<sub>2</sub>_uncertainty</td><td>COAL</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to COAL use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>OIL</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to OIL use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>GAS</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to GAS use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr><tr><td>BUNKER</td><td>Annual uptake of O<sub>2</sub> in each cell (kg O<sub>2</sub> year-1) due to BUNKER use</td><td>kg O<sub>2</sub> year-1</td><td>60</td><td>1,800</td><td>3,600</td></tr></tbody></table>

$$EDGAR_{s,m,i,j} = \begin{cases} a_{s,i,j} + (b_{s,i,j} \cdot \text{HCDDfrac}_{m,i,j}), & R_{s,i_r,j_r}^2 \ge 0.66 \\ EDGAR_{s,m,i,j}, & \text{else} \end{cases} \quad (7)$$

Here HCDDfrac, a and b were re-gridded by repeating each grid cell in the i_r, j_r dimensions to provide output at the resolution of the EDGAR grid (0.1° × 0.1°; i, j).

**GCP-GridFED Protocol.** CO₂ Emissions. GCP-GridFED was generated using the six-step emissions scaling protocol set out below and applied sequentially for each year in the period 1959–2018 (see Fig. 1):

1. **Group emissions from EDGAR sectors by source class.** The global gridded (i, j) monthly (m) CO₂ emissions were summed across the EDGAR activity sectors (s) in each source class used in this study (S; see Table 1). The monthly distribution of annual emissions was adjusted in advance using Eq. 7.

$$EDGAR_{s,m,i,j} = \sum EDGAR_{s,m,i,j} \quad (8)$$

2. **Extract gridded emissions data from EDGAR for each country.** A subset of gridded monthly CO₂ emissions from each GCP-GridFED source class (see Table 1) was extracted for each country (c) using country masks (True/False) from the EU-GISCO dataset⁴⁹. No subset was extracted for the bunker fuels source
<table><caption>Table 5. Regional summary statistics relating to total annual CO₂ emissions from GridFEDv2019.1.</caption><thead><tr><th>Variable</th><th>Period</th><th>Global</th><th>Northern<br>Extratropics</th><th>Tropics</th><th>Southern<br>Extratropics</th></tr></thead><tbody><tr><td rowspan="7">Mean Annual Emissions<br>(Gt CO<sub>2</sub> year<sup>-1</sup>)</td><td>1960–1969</td><td>11.1</td><td>9.6</td><td>1.3</td><td>0.2</td></tr><tr><td>1970–1979</td><td>17.2</td><td>14.4</td><td>2.5</td><td>0.3</td></tr><tr><td>1980–1989</td><td>20.1</td><td>16.2</td><td>3.5</td><td>0.3</td></tr><tr><td>1990–1999</td><td>23.2</td><td>17.7</td><td>5.1</td><td>0.4</td></tr><tr><td>2000–2009</td><td>28.5</td><td>20.4</td><td>7.5</td><td>0.5</td></tr><tr><td>2010–2018</td><td>34.9</td><td>23.3</td><td>10.9</td><td>0.6</td></tr><tr><td>2018</td><td>36.4</td><td>23.8</td><td>11.9</td><td>0.6</td></tr><tr><td rowspan="7">Mean Fraction of Global<br>Emissions (%)</td><td>1960–1969</td><td></td><td>86.1</td><td>12.1</td><td>1.5</td></tr><tr><td>1970–1979</td><td></td><td>83.9</td><td>14.3</td><td>1.6</td></tr><tr><td>1980–1989</td><td></td><td>80.7</td><td>17.3</td><td>1.7</td></tr><tr><td>1990–1999</td><td></td><td>76.1</td><td>21.8</td><td>1.8</td></tr><tr><td>2000–2009</td><td></td><td>71.5</td><td>26.3</td><td>1.9</td></tr><tr><td>2010–2018</td><td></td><td>66.7</td><td>31.2</td><td>1.7</td></tr><tr><td>2018</td><td></td><td>65.2</td><td>32.8</td><td>1.7</td></tr><tr><td rowspan="7">Mean Annual Emissions<br>Uncertainty (Gt CO<sub>2</sub> year<sup>-1</sup>)</td><td>1960–1969</td><td>2.7</td><td>2.2</td><td>0.4</td><td>0.0</td></tr><tr><td>1970–1979</td><td>4.0</td><td>3.1</td><td>0.8</td><td>0.1</td></tr><tr><td>1980–1989</td><td>5.1</td><td>3.8</td><td>1.2</td><td>0.1</td></tr><tr><td>1990–1999</td><td>6.1</td><td>4.2</td><td>1.8</td><td>0.1</td></tr><tr><td>2000–2009</td><td>7.9</td><td>5.0</td><td>2.7</td><td>0.1</td></tr><tr><td>2010–2018</td><td>10.3</td><td>6.1</td><td>4.1</td><td>0.2</td></tr><tr><td>2018</td><td>10.8</td><td>6.1</td><td>4.5</td><td>0.2</td></tr><tr><td rowspan="7">Mean Annual Emissions<br>Uncertainty (%)</td><td>1960–1969</td><td>24.6</td><td>23.5</td><td>32.4</td><td>26.2</td></tr><tr><td>1970–1979</td><td>23.2</td><td>21.6</td><td>32.7</td><td>24.8</td></tr><tr><td>1980–1989</td><td>25.4</td><td>23.3</td><td>35.3</td><td>25.6</td></tr><tr><td>1990–1999</td><td>26.3</td><td>23.7</td><td>35.7</td><td>25.7</td></tr><tr><td>2000–2009</td><td>27.8</td><td>24.7</td><td>36.5</td><td>25.5</td></tr><tr><td>2010–2018</td><td>29.5</td><td>26.0</td><td>37.2</td><td>25.2</td></tr><tr><td>2018</td><td>29.6</td><td>25.8</td><td>37.5</td><td>24.9</td></tr><tr><td rowspan="7">Mean Annual Emissions<br>Growth (%)</td><td>1960–1969</td><td>4.5</td><td>4.4</td><td>5.0</td><td>5.0</td></tr><tr><td>1970–1979</td><td>3.6</td><td>3.3</td><td>5.7</td><td>3.5</td></tr><tr><td>1980–1989</td><td>1.3</td><td>0.9</td><td>3.2</td><td>2.6</td></tr><tr><td>1990–1999</td><td>0.9</td><td>0.2</td><td>3.5</td><td>2.3</td></tr><tr><td>2000–2009</td><td>2.6</td><td>1.8</td><td>4.7</td><td>2.0</td></tr><tr><td>2010–2018</td><td>1.7</td><td>1.1</td><td>2.9</td><td>0.6</td></tr><tr><td>2018</td><td>2.1</td><td>1.7</td><td>3.1</td><td>-0.2</td></tr></tbody></table>

class; the entire grid layer was scaled globally. Hence, all grid cells were included in an 'international' mask (True in all cells) and treated thereafter in the same way as each country.

$$EDGAR_{S,m,i[c],j[c]} = EDGAR_{S,m,i[True],j[True]} \qquad (9)$$

3. **Sum EDGAR emissions for each GCP-GridFED source class.** Monthly CO₂ emissions were summed both across the months of the year and across the grid extracted for each nation, for each GCP-GridFED source class. The resulting annual emission sub-totals were stored in a tabular format matching the structure of the GCB-NAE data.

$$EDGAR_{S,c} = \sum \sum EDGAR_{S,m,i[c],j[c]} \qquad (10)$$

4. **Calculate scaling factors based on comparison of EDGAR and GCB-NAE emissions estimates.** For each country c and for each source class S, the scaling factor (α) required to convert the annual CO₂ emissions from EDGAR (step 3) to the annual CO₂ emissions estimate from GCB-NAE was derived as follows.

$$GCB_{S,c} = \alpha_{S,c} \cdot EDGAR_{S,c} \qquad (11)$$

5. **Apply annual scaling factors to monthly emission grids.** The scaling factors for each nation and GCP-GridFED source class were applied to the national monthly CO₂ emissions grids generated in step 2.
<table><caption>**Table 6.** Regional summary statistics relating to seasonality of total CO₂ emissions from GridFEDv2019.1.</caption><thead><tr><th>Variable</th><th>Period</th><th>Global</th><th>Northern<br>Extratropics</th><th>Tropics</th><th>Southern<br>Extratropics</th></tr></thead><tbody><tr><td rowspan="6">Mean Daily Emissions<br>(Mt CO<sub>2</sub> day<sup>-1</sup>)</td><td>1960–1969</td><td>30.5</td><td>26.2</td><td>3.7</td><td>0.5</td></tr><tr><td>1970–1979</td><td>47.1</td><td>39.5</td><td>6.7</td><td>0.7</td></tr><tr><td>1980–1989</td><td>55.1</td><td>44.5</td><td>9.5</td><td>0.9</td></tr><tr><td>1990–1999</td><td>63.6</td><td>48.4</td><td>13.9</td><td>1.2</td></tr><tr><td>2000–2009</td><td>78.1</td><td>55.8</td><td>20.5</td><td>1.5</td></tr><tr><td>2010–2018</td><td>95.6</td><td>63.8</td><td>29.8</td><td>1.7</td></tr><tr><td rowspan="6">Mean Seasonal Amplitude<br>(Mt CO<sub>2</sub> day<sup>-1</sup>)</td><td>1960–1969</td><td>6.3</td><td>6.1</td><td>0.1</td><td>0.3</td></tr><tr><td>1970–1979</td><td>10.1</td><td>9.8</td><td>0.2</td><td>0.5</td></tr><tr><td>1980–1989</td><td>11.0</td><td>10.7</td><td>0.2</td><td>0.5</td></tr><tr><td>1990–1999</td><td>12.2</td><td>11.7</td><td>0.2</td><td>0.8</td></tr><tr><td>2000–2009</td><td>13.9</td><td>13.1</td><td>0.3</td><td>1.4</td></tr><tr><td>2010–2018</td><td>15.5</td><td>14.0</td><td>0.4</td><td>1.9</td></tr><tr><td rowspan="6">Mean Seasonal Amplitude (%)</td><td>1960–1969</td><td>20.9</td><td>23.6</td><td>20.7</td><td>7.7</td></tr><tr><td>1970–1979</td><td>21.6</td><td>25.0</td><td>21.2</td><td>7.6</td></tr><tr><td>1980–1989</td><td>20.0</td><td>24.1</td><td>21.0</td><td>5.6</td></tr><tr><td>1990–1999</td><td>19.1</td><td>24.2</td><td>20.9</td><td>6.1</td></tr><tr><td>2000–2009</td><td>18.0</td><td>23.4</td><td>21.0</td><td>6.8</td></tr><tr><td>2010–2018</td><td>16.3</td><td>22.1</td><td>21.6</td><td>6.5</td></tr></tbody></table>

The same scaling factor was used for all months. For the bunker fuels source class, the scaling factor was applied to the equivalent global data.

$$ \text{GridFED}_{S,m,i[c],j[c]} = \alpha_{S,c} \cdot \text{EDGAR}_{S,m,i[c],j[c]} \quad (12) $$

**6. Collate national data to a global output.** Scaled monthly CO₂ emissions grids from all nations were merged into a single grid for each GCP-GridFED source class.

We do not attempt to adjust the EDGARv4.3.2 grids (year 2010) for a range of historical changes to the spatial distribution of emissions, for instance due to the expansion of road networks or flight routes, the commissioning/decommissioning of facilities or large-scale population migration. The resolution of these issues will be prioritised in future developments to the GCP-GridFED protocol. We note that developments introduced in EDGARv5.0 (ref. ⁴²) include refined spatial proxy records and national temporal profiles covering the period 1970–2012, which will support further developments to the GCP-GridFED protocol. Dedicated datasets of fuel-specific monthly CO₂ emissions are also emerging for some countries, including India⁵⁷ and the USA⁵⁸, and could be used preferentially in the GCP-GridFED protocol. Additional sources such as the diffusive coal mine oxidation CO₂, as derived for the dataset CHE-EDGARv4.3.2_FT2015⁴¹,⁵⁹,⁶⁰ will also be considered.

We do not consider emissions of non-CO₂ carbon emissions that later influence atmospheric CO₂ (in particular, CO and CH₄). Here all fossil carbon is assumed to be emitted as fossil CO₂, whereas a fraction is in reality emitted as CO and later represents a diffuse fossil CO₂ source after oxidation to CO₂ (~1.8 Pg CO₂-equivalent year⁻¹)⁶¹. In GridFEDv2019.1, the diffuse nature of this CO₂ source is not considered, and the source is instead placed at the surface at the time and location of oxidation. Meanwhile, fossil CH₄ fugitive emissions represent an additional diffuse source of CO₂ emissions (0.4 Pg CO₂-equivalent year⁻¹) that is not considered here⁶². These diffuse CO₂ sources will also be considered in future developments to the GCP-GridFED protocol.

**CO₂ emissions uncertainty.** We provide gridded uncertainties to complement all gridded layers of the GCP-GridFED dataset, however we note here the incomplete nature of our uncertainty assessment. The gridded uncertainties are based on the total fossil CO₂ emissions uncertainty assessment from the GCB³, combined with variation in relative uncertainties across emission sectors from the recent TNO assessment (Table 2)³⁴ or uncertainties in national total CO₂ emissions, we adopt the values presented in the uncertainty assessment of the GCB; 5% for the 42 Annex I countries that report annually to the UNFCCC⁴⁴ and 10% for other countries³,⁶³ (1σ). Annex I countries are assigned lower uncertainty because for these countries more detailed energy and activity statistics are available, and they are periodically reviewed externally³. We used data presented in the TNO uncertainty assessment to evaluate the ratio of the uncertainties for each sector (U_TNOₐ) to the uncertainty in total emissions (U_TNOₐₜ). We then scaled the ratios to the uncertainties in total emissions that are adopted for Annex I and other countries from GCB-NAE.

$$ U_{\text{GridFEDs}} = \left\{ \frac{U_{\text{TNOs}}}{U_{\text{TNOTot}}} \times 5\%, \text{ UNFCCC Annex I } \frac{U_{\text{TNOs}}}{U_{\text{TNOTot}}} \times 10\%, \text{ other} \right\} \quad (13) $$
<table><caption>**Table 7.** Regional summary statistics relating to annual CO₂ emissions from each source, from GridFEDv2019.1.</caption><thead><tr><th>Variable</th><th>Period</th><th>Global</th><th>Northern Extratropics</th><th>Tropics</th><th>Southern Extratropics</th></tr></thead><tbody><tr><td rowspan="7">Mean Daily Emissions (Mt CO<sub>2</sub> day<sup>-1</sup>)</td><td>1960–1969</td><td>30.5</td><td>26.2</td><td>3.7</td><td>0.5</td></tr><tr><td>1970–1979</td><td>47.1</td><td>39.5</td><td>6.7</td><td>0.7</td></tr><tr><td>1980–1989</td><td>55.1</td><td>44.5</td><td>9.5</td><td>0.9</td></tr><tr><td>1990–1999</td><td>63.6</td><td>48.4</td><td>13.9</td><td>1.2</td></tr><tr><td>2000–2009</td><td>78.1</td><td>55.8</td><td>20.5</td><td>1.5</td></tr><tr><td>2010–2018</td><td>95.6</td><td>63.8</td><td>29.8</td><td>1.7</td></tr><tr><td>2018</td><td>99.8</td><td>65.1</td><td>32.7</td><td>1.7</td></tr><tr><td rowspan="7">Mean Seasonal Amplitude (Mt CO<sub>2</sub> day<sup>-1</sup>)</td><td>1960–1969</td><td>6.3</td><td>6.1</td><td>0.1</td><td>0.3</td></tr><tr><td>1970–1979</td><td>10.1</td><td>9.8</td><td>0.2</td><td>0.5</td></tr><tr><td>1980–1989</td><td>11.0</td><td>10.7</td><td>0.2</td><td>0.5</td></tr><tr><td>1990–1999</td><td>12.2</td><td>11.7</td><td>0.2</td><td>0.8</td></tr><tr><td>2000–2009</td><td>13.9</td><td>13.1</td><td>0.3</td><td>1.4</td></tr><tr><td>2010–2018</td><td>15.5</td><td>14.0</td><td>0.4</td><td>1.9</td></tr><tr><td rowspan="7">Mean Seasonal Amplitude (%)</td><td>1960–1969</td><td>20.9</td><td>23.6</td><td>20.7</td><td>7.7</td></tr><tr><td>1970–1979</td><td>21.6</td><td>25.0</td><td>21.2</td><td>7.6</td></tr><tr><td>1980–1989</td><td>20.0</td><td>24.1</td><td>21.0</td><td>5.6</td></tr><tr><td>1990–1999</td><td>19.1</td><td>24.2</td><td>20.9</td><td>6.1</td></tr><tr><td>2000–2009</td><td>18.0</td><td>23.4</td><td>21.0</td><td>6.8</td></tr><tr><td>2010–2018</td><td>16.3</td><td>22.1</td><td>21.6</td><td>6.5</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Fig. 5** Monthly emissions anomaly relative to the annual mean daily emissions rate from GCP-GridFEDv2019.1. Each grey line represents a year of data in the period 1959–2018. The red line marks the mean value across the time series. (Top row, left column) Global seasonality is disaggregated to (other columns) the top 4 emission regions and (other rows) source classes. No seasonality is present in India, based on the EDGARv4.3.2 gridded input data for the year 2010.</figcaption></figure>

Table 2 shows sectoral TNO uncertainty estimates and the resulting uncertainties adopted in GCP-GridFED for Annex I and other countries, for each sector.

The gridded uncertainty estimates presented here do not include the uncertainties associated with the spatial or temporal (monthly) disaggregation of national emissions, nor do we present a formal assessment of those
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Time series of (left column) annual O₂ uptake (Gt O₂ year⁻¹) and (right column) monthly O₂ uptake (Mt O₂ day⁻¹) through fossil fuel oxidation as estimated by GCP-GridFED. (Top row) Total global uptake is disaggregated to (other rows) the top 4 emission regions. The plotted uptake uncertainties are based on CO₂ emissions uncertainties of 5% for Annex I nations and global³, and they do not include uncertainties in oxidative ratios.</figcaption></figure>

disaggregation uncertainties. We note that spatially-averaged uncertainties resulting from the spatial disaggregation of national emissions estimates to grid cells are on the order of 20–75% (1σ) at spatial resolutions of 1 km to 1° (refs. ⁵³,⁶⁴–⁶⁷). Spatial disaggregation uncertainties occur due to incomplete proxy data coverage (e.g. unmapped or mislocated point sources), poorly constrained nonlinearities (e.g. differences in the emissions intensity between equally dense rural and urban populations), shortcomings in continuous proxy values (e.g. poorly constrained population density) or inappropriate spatial representativeness (e.g. the spatial representativeness of
<figure><img src="image_9.png" /><figcaption>**Fig. 7** Time series of annual O₂ uptake as estimated by GCP-GridFEDv2019.1 for the period 1959–2018 (Gt O₂ year⁻¹). The plotted uptake uncertainties are based on CO₂ emissions uncertainties of 5% for Annex I nations and global³, and they do not include uncertainties in oxidative ratios. (Top row, Left column) Total global uptake is disaggregated to the (other columns) top 4 emission regions and (other rows) source classes.</figcaption></figure>

roadmaps for traffic volume). By construction, these uncertainties are larger for years distant from our reference year 2010 and at monthly resolution. Dedicated analyses of regional emissions at high temporal resolution are yielding new data with which to quantify temporal disaggregation uncertainties⁵⁷,⁶⁸ and to assess the robustness of the temporal profiles employed here and elsewhere⁴².

A full quantitative assessment of these issues, to support the development of comprehensive grid-level uncertainties associated with GCP-GridFED, will be the subject of future work. Overall, our approach to uncertainty quantification is broadly representative of the sectoral contributions to total emissions in each grid cell, which changes throughout the time series. Inversion models may utilise these uncertainty grids but with the freedom to build more complex covariance structures to suit their requirements.

**O₂ Uptake.** The relationship between CO₂ and O₂ fluxes during oxidation reactions can be expressed as an oxidative ratio (OR = flux of O₂ from the atmosphere/flux of CO₂ to the atmosphere, unitless)³⁶,⁶⁹. The OR differs detectably between specific fossil fuel sources, holding a value of −1.17 for coal, −1.44 for oil, and −1.95 for natural gas³⁶,⁶⁹. Uncertainties in OR are thought to be on the order of 2–3%, however variations within fuel classes, such as different grades of coal, have not been studied extensively (ref. ³⁵). Cement clinker production involves a calcination reaction rather than an oxidation reaction, and thus no exchange of oxygen occurs (OR = 0).
<figure><img src="image_10.png" /><figcaption>Fig. 8 Time series of monthly O₂ uptake as estimated by GCP-GridFEDv2019.1 for the period 1990–2018 (Mt O₂ day⁻¹). (Top row, Left column) Total global uptake is disaggregated to the (other columns) top 4 emission regions and (other rows) source classes. The plotted uptake uncertainties are based on CO₂ emissions uncertainties of 5% for Annex I nations and global³, and they do not include uncertainties in oxidative ratios.</figcaption></figure>

GCP-GridFED calculates gridded estimates of the uptake of O₂ during fossil fuel oxidation by applying OR values to the CO₂ emissions estimates for each source.

$$ \text{GridFED\_O2}_{S,i,j} = \text{GridFED}_{S,i,j} \cdot \text{OR}_S \qquad (14) $$

We treat relative uncertainty in O₂ emissions as equal to the relative uncertainty in CO₂ emissions (U_GridFED).

## Data Records

All GCP-GridFEDv2019.1 output grids can be accessed via the Zenodo data repository⁷⁰ (https://doi.org/10.5281/zenodo.3958283).

The data records include 60 files in Network Common Data Form (NetCDF) format with the naming convention GCPGlobal_{YYYY}.nc, where YYYY is the year represented by the contents. Each NetCDF file includes 3 dimensions: time (month of the year expressed as days since the first day of YYYY, n = 12); latitude (Degrees North of the equator [cell centres], n = 1800); longitude (Degrees East of the Prime Meridian [cell centres], n = 3600). Each NetCDF file includes three groups representing CO₂ emissions, CO₂ emissions uncertainty, and O₂ uptake (CO₂, CO₂_uncertainty and O₂, respectively). Each group contains five variables representing
<figure><img src="image_11.png" /><figcaption>**Fig. 9** Spatial distribution of CO₂ emissions as estimated by GCP-GridFED (kg CO₂ year⁻¹). Gridded (0.1° × 0.1°) estimates of total fossil CO₂ emissions are shown for four years of the analytical period.</figcaption></figure>

emissions from each source class (COAL, OIL, GAS, CEMENT, BUNKER) with the units shown in Table 3. Each file contains 1,088,640,000 unique data points. All 60 NetCDF files are contained within a.zip archive named “GCP-GridFEDv2019.1_monthly.zip”.

The data records also include 1 file in NetCDF format, “GCPGlobal_annual.nc”. The NetCDF file includes 3 dimensions: time (year expressed as days since 1959–01–01, n = 60); latitude (Degrees North of the equator [cell centres], n = 1800); longitude (Degrees East of the Prime Meridian [cell centres], n = 3600). Each NetCDF file includes 4 groups representing CO₂ emissions, CO₂ emissions uncertainty, and O₂ uptake, and O₂ uptake uncertainty (CO₂, CO₂_uncertainty, O₂, and O₂_uncertainty respectively). Each group contains 5 variables representing emissions from each source class (COAL, OIL, GAS, CEMENT, BUNKER) with the units shown in Table 4.
<figure><img src="image_12.png" /><figcaption>**Fig. 10** (Left panels) Spatial distribution of O₂ uptake through fossil fuel use as estimated by GCP-GridFEDv2019.1 (kg O₂ year⁻¹). Gridded (0.1° × 0.1°) estimates of total O₂ uptake are shown for four years of the time series (1959–2018). (Right panels) Spatially-explicit oxidative ratios (OR; kg O₂ kg⁻¹ CO₂) for total emissions activities as estimated by GridFEDv2019.1.</figcaption></figure>

The file contains 6,998,400,000 unique data points. The NetCDF file is contained within a.zip archive named "GCP-GridFEDv2019.1_annual.zip".

"GCP-GridFEDv2019.1_monthly.zip" and "GCP-GridFEDv2019.1_annual.zip" can be found within a parent. zip file name "GCP-GridFEDv2019.1.zip". All grids are bottom-left arranged with coordinates referenced to the prime meridian and the equator.

## Technical Validation

We provide Figs. 2–12, the summary statistics in Tables 4–7 and Online-Only Table 1 to outline the key features of GCP-GridFEDv2019.1 and assist with its technical validation.

GCP-GridFED is designed to distribute national annual emissions from GCB-NAE over a spatio-temporal grid based on EDGARv4.3.2. We validated the outputs from GCP-GridFED by comparing the global annual emissions from the output grids (the sum of emissions across the global grid) with the input data supplied to the gridding protocol from GCB-NAE. Throughout the time series of emissions and across all source classes, the global annual emissions totals from GCP-GridFED were always within 0.0077% of the GCB-NAE input data throughout the annual time series (Figs. 2 and 3). The discrepancies were caused by unscalable (zero or NoData) values in sectors of the EDGAR dataset at the national level in 13 countries (EDGAR data summed within the national masks as per Eq. 10). These 13 countries make a small contribution to total global emissions (0.047% in 2018). Online-Only Table 1 provides national-level comparisons of the emissions estimates from GCP-GridFED and GCB-NAE. For the 13 countries where maximum absolute discrepancies exceeded 1% of GCB-NAE emissions, we provide a brief description of the cause of the discrepancy. The GCP-GridFED outputs are robust to within 0.0001% of GCB-NAE values in 195 countries, plus bunker fuels, comprising 99.9% of global emissions in 2018. Hence, we conclude that the GCP-GridFEDv2019.1 is consistent with GCB-NAE emissions estimates for the years 1959–2018.

We also observed a close match between the seasonality seen in the year 2010 in the GCP-GridFED dataset and that seen in the same year of the EDGAR input data, both at the global scale and in large austral and boreal extratropical nations (Fig. 12). This coherence indicates that the seasonality seen in the EDGAR dataset was preserved by the GCP-GridFED protocol. Inter-annual variability in the monthly distribution of emissions can be seen most prominently in the EU27 + UK. Note that EDGARv4.3.2 does not feature monthly variability in emissions for tropical countries⁴¹, and so GCP-GridFED also shows no seasonality in these countries.
<figure><img src="image_13.png" /><figcaption>**Fig. 11** Gridded (0.1° × 0.1°) estimates of relative uncertainty in total CO₂ emissions for four years of the GCP-GridFEDv2019.1 time series. Uncertainty in total emissions is aggregated from the sector-level estimates (see Table 2). The uncertainty estimates account for uncertainty across national emission reports and spatial differences in the sectoral breakdown to total emissions in each grid cell, which changes throughout the time series, however they exclude uncertainties associated with the spatial or temporal (monthly) disaggregation of national emissions (see 'CO₂ emissions uncertainty'). Aggregation of uncertainties to a coarser resolution should account for the non-independence of gridded emissions uncertainties.</figcaption></figure>

**Usage Notes**

The data is intended for use as a prior in inversion model studies, which may wish to incorporate individual priors for each source class or to use total gridded emissions. The data records contain a layer for each source class. Global total emissions can be calculated as the sum of emissions across the 5 source classes. National total
<figure><img src="image_14.png" /><figcaption>**Fig. 12** Monthly emissions anomaly relative to the annual mean daily emissions rate from GCP-GridFEDv2019.1. Each grey line represents a year of data in the period 1959–2018. The red line marks the top row, left column) Global seasonality is disaggregated to (other columns) source classes (other rows) regions.</figcaption></figure>

emissions estimates should be calculated as the sum of coal, oil, gas and cement emissions (bunker fuel emissions should not be included in national emissions totals)⁷¹.

GCP-GridFED will be updated annually and made available for the inversion model runs conducted annually as part of the GCP assessment of the GCB. An updated version of GCP-GridFED (GCP-GridFEDv2020.1) was already made available upon request to support the inversion model runs of the GCP's 2020 GCB assessment⁷²,⁷³ and is now publicly available⁷⁴. GCP-GridFEDv2020.1 is based on the emissions estimates from a preliminary
release of GCB-NAE covering the years 1959–2019 with input data available to June 2020. Further updates will be issued as the GCB-NAE data is updated⁷⁴.

When using GCP-GridFED as a prior in inversion models operating at a coarser resolution, aggregation to the required resolution should account for the non-independence of gridded emissions uncertainties. See ‘CO₂ Emissions Uncertainty’ for further information regarding our treatment of spatial and temporal aggregation in GCP-GridFED.

## Code availability

The code used to perform all steps described here and shown in Fig. 1 can be accessed via the Zenodo dataset repository entry for GCP-GridFEDv2020.1 (https://doi.org/10.5281/zenodo.4277267)⁷⁴. GCP-GridFEDv2020.1 uses the same code and methodology as GCP-GridFEDv2019.1 but includes updated estimates of national annual emissions through to 2019 from the GCP, as discussed in the Usage Notes and also detailed at ref. ⁷⁴.

Received: 29 July 2020; Accepted: 24 November 2020;

Published online: 07 January 2021

## References

1. Joos, F. & Spahni, R. Rates of change in natural and anthropogenic radiative forcing over the past 20,000 years. *Proc. Natl. Acad. Sci. 105*, 1425–1430 (2008).

2. IPCC. *Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change*. (2014).

3. Friedlingstein, P. et al. Global Carbon Budget 2019. *Earth Syst. Sci. Data* **11**, 1783–1838 (2019).

4. Le Quéré, C. et al. Drivers of declining CO₂ emissions in 18 developed economies. *Nat. Clim. Chang.* **9**, 213–217 (2019).

5. Peters, G. P. et al. Towards real-time verification of CO₂ emissions. *Nat. Clim. Chang.* **7**, 848–850 (2017).

6. Ciais, P. et al. Carbon and Other Biogeochemical Cycles. In *Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change* (ed. Intergovernmental Panel on Climate Change) 465–570 (Cambridge University Press, 2013).

7. Denman, K. L. et al. Couplings Between Changes in the Climate System and Biogeochemistry. In *Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change* (eds. Solomon, S. et al.) 499–587 (Cambridge University Press, Cambridge, UK and New York, USA, 2007).

8. Prentice, I. C. et al. The Carbon Cycle and Atmospheric Carbon Dioxide. In *Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change* (eds. Houghton, J. T. et al.) 183–237 (Cambridge University Press, 2001).

9. Schimel, D. et al. Radiative Forcing of Climate Change. In *Climate Change 1995 The Science of Climate Change. Contribution of Working Group I to the Second Assessment Report of the Intergovernmental Panel on Climate Change* (eds. Houghton, J. T. et al.) (Cambridge University Press, 1995).

10. Watson, R. T., Rodhe, H., Oeschger, H., Siegenthaler, U. & Press, C. U. Greenhouse Gases and Aerosols. In *Climate Change: The IPCC Scientific Assessment. Intergovernmental Panel on Climate Change (IPCC)* (eds. Houghton, J. T., Jenkins, G. J. & Ephraums, J. J.) 1–40 (Cambridge University Press, 1990).

11. Le Quéré, C. et al. The global carbon budget 1959–2011. *Earth Syst. Sci. Data* **5**, 165–185 (2013).

12. Canadell, J. G. et al. Contributions to accelerating atmospheric CO₂ growth from economic activity, carbon intensity, and efficiency of natural sinks. *Proc. Natl. Acad. Sci. 104*, 18866–18870 (2007).

13. Friedlingstein, P. et al. Update on CO₂ emissions. *Nat. Geosci.* **3**, 811–812 (2010).

14. BP. *Statistical Review of World Energy*, 2019. https://www.bp.com/content/dam/bp/business-sites/en/global/corporate/pdfs/energy-economics/statistical-review/bp-stats-review-2019-full-report.pdf (2019).

15. Gilfillan, D., Marland, G., Boden, T. & Andres, R. *Global, Regional, and National Fossil-Fuel CO₂ Emissions*. https://energy.appstate.edu/CDIAC (2019).

16. Andrew, R. M. Global CO₂ emissions from cement production, 1928–2018. *Earth Syst. Sci. Data* **11**, 1675–1710 (2019).

17. Hansis, E., Davis, S. J. & Pongratz, J. Relevance of methodological choices for accounting of land use change carbon fluxes. *Global Biogeochem. Cycles* **29**, 1230–1246 (2015).

18. Houghton, R. A. & Nassikas, A. A. Global and regional fluxes of carbon from land use and land cover change 1850–2015. *Global Biogeochem. Cycles* **31**, 456–472 (2017).

19. Sitch, S. et al. Recent trends and drivers of regional sources and sinks of carbon dioxide. *Biogeosciences* **12**, 653–679 (2015).

20. Hauck, J. et al. Consistency and Challenges in the Ocean Carbon Sink Estimate for the Global Carbon Budget. *Front. Mar. Sci.* **7**, 1–33, https://doi.org/10.3389/fmars.2020.571720 (2020).

21. Stephens, B. B. et al. Weak northern and strong tropical land carbon uptake from vertical profiles of atmospheric CO₂. *Science* (80–). **316**, 1732–1735 (2007).

22. Gurney, K. R. et al. Towards robust regional estimates of annual mean [CO]₂ sources and sinks. *Nature* **415**, 626–630 (2002).

23. Rödenbeck, C., Zaehle, S., Keeling, R. & Heimann, M. How does the terrestrial carbon exchange respond to inter-annual climatic variations? A quantification based on atmospheric CO₂ data. *Biogeosciences* **15**, 2481–2498 (2018).

24. Van Der Laan-Luijkx, I. T. et al. The CarbonTracker Data Assimilation Shell (CTDAS) v1.0: Implementation and global carbon balance 2001–2015. *Geosci. Model Dev.* **10**, 2785–2800 (2017).

25. Chevallier, F. et al. Inferring CO₂ sources and sinks from satellite observations: Method and application to TOVS data. *J. Geophys. Res.* **110**, D24309 (2005).

26. Peters, W. et al. Seven years of recent European net terrestrial carbon dioxide exchange constrained by atmospheric observations. *Glob. Chang. Biol.* **16**, 1317–1337 (2010).

27. Basu, S., Miller, J. B. & Lehman, S. Separation of biospheric and fossil fuel fluxes of CO₂ by atmospheric inversion of CO₂ and 14CO₂ measurements: Observation System Simulations. *Atmos. Chem. Phys.* **16**, 5665–5683 (2016).

28. Saeki, T. & Patra, P. K. Implications of overestimated anthropogenic CO₂ emissions on East Asian and global land CO₂ flux inversion. *Geosci. Lett.* **4**, 9 (2017).

29. Gaubert, B. et al. Global atmospheric CO₂ inverse models converging on neutral tropical land exchange, but disagreeing on fossil fuel and atmospheric growth rate. *Biogeosciences* **16**, 117–134 (2019).

30. Rödenbeck, C., Houweling, S., Gloor, M. & Heimann, M. CO₂ flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport. *Atmos. Chem. Phys.* **3**, 1919–1964 (2003).

31. Chevallier, F. et al. CO₂ surface fluxes at grid point scale estimated from a global 21 year reanalysis of atmospheric measurements. *J. Geophys. Res.* **115**, D21307 (2010).

32. Peters, W. et al. An ensemble data assimilation system to estimate CO₂ surface fluxes from atmospheric trace gas observations. *J. Geophys. Res.* **110**, D24304 (2005).
33. Friedlingstein, P. et al. Supplemental data of the Global Carbon Budget 2019. ICOS-ERIC Carbon Portal https://doi.org/10.18160/gcp-2019 (2019).

34. Marshall, J. & Ramirez, T. CHE Project Deliverable 4.3 Attribution Problem Configurations. https://www.che-project.eu/sites/default/files/2020-01/CHE-D4-3-V4-1.pdf (2019).

35. Keeling, R. F. & Manning, A. C. Studies of Recent Changes in Atmospheric O₂ Content. In Treatise on Geochemistry (eds. Holland, H. D. & Turekian, K. K.) 5, 385–404 (2014).

36. Steinbach, J. et al. The CO₂ release and Oxygen uptake from Fossil Fuel Emission Estimate (COFFEE) dataset: Effects from varying oxidative ratios. *Atmos. Chem. Phys.* **11**, 6855–6870 (2011).

37. Gruber, N., Gloor, M., Fan, S. M. & Sarmiento, J. L. Air-sea flux of oxygen estimated from bulk data: Implications for the marine and atmospheric oxygen cycles. *Global Biogeochem. Cycles* **15**, 783–803 (2001).

38. Tohjima, Y., Mukai, H., Nojiri, Y., Yamagishi, H. & Machida, T. Atmospheric O₂/N₂ measurements at two Japanese sites: Estimation of global oceanic and land biotic carbon sinks and analysis of the variations in atmospheric potential oxygen (APO). *Tellus, Ser. B Chem. Phys. Meteorol.* **60 B**, 213–225 (2008).

39. Rödenbeck, C., Le Quéré, C., Heimann, M. & Keeling, R. F. Interannual variability in oceanic biogeochemical processes inferred by inversion of atmospheric O₂/N₂ and CO₂ data. *Tellus, Ser. B Chem. Phys. Meteorol.* **60 B**, 685–705 (2008).

40. Goto, D., Morimoto, S., Aoki, S., Patra, P. K. & Nakazawa, T. Seasonal and short-term variations in atmospheric potential oxygen at Ny-Ålesund, Svalbard. *Tellus, Ser. B Chem. Phys. Meteorol.* **69**, 1–11 (2017).

41. Janssens-Maenhout, G. et al. EDGAR v4.3.2 Global Atlas of the three major greenhouse gas emissions for the period 1970–2012. *Earth Syst. Sci. Data* **11**, 959–1002 (2019).

42. Crippa, M. et al. High resolution temporal profiles in the Emissions Database for Global Atmospheric Research. *Sci. Data* **7**, 1–17 (2020).

43. Keeling, C. D. et al. Atmospheric Carbon-Dioxide Variations at Mauna-Loa Observatory, Hawaii. *Tellus* **28**, 538–551 (1976).

44. UNFCCC. National Inventory Submissions. 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. (2019).

45. UNSD. United Nations Statistics Division: Energy Statistics. http://unstats.un.org/unsd/energy/. (2019).

46. Rypdal, K. & Paciornik, N. 2006 IPCC Guidelines for National Greenhouse Gas Inventories, Volume 1: General Guidance and Reporting, Prepared by the National Greenhouse Gas Inventories Programme. (Intergovernmental Panel on Climate Change, 2006).

47. IPCC. Revised 1996 IPCC Guidelines for National Greenhouse Gas Inventories IPCC/OECD/IEA. (1996).

48. Janssens-Maenhout, G., Pagliari, V., Guizzardi, D. & Muntean, M. Global emission inventories in the Emission Database for Global Atmospheric Research (EDGAR) – Manual (I): Gridding: EDGAR emissions distribution on global grid-maps, JRC Report, EUR 25785 EN. https://publications.jrc.ec.europa.eu/repository (2013).

49. European Commission Eurostat (ESTAT) GISCO. Countries, 2016 - Administrative Units - Dataset ID 5C27B6C0-BC1C-4175-9B0B-783AEEBAAD61. http://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units (2018).

50. Peters, G. P. et al. Carbon dioxide emissions continue to grow amidst slowly emerging climate policies. *Nat. Clim. Chang.* **10**, 3–6 (2020).

51. Ummel, K. Carma Revisited: An Updated Database of Carbon Dioxide Emissions from Power Plants Worldwide. *SSRN Electron. J.* https://doi.org/10.2139/ssrn.2226505 (2012).

52. Oda, T., Maksyutov, S. & Andres, R. J. The Open-source Data Inventory for Anthropogenic CO₂, version 2016 (ODIAC2016): A global monthly fossil fuel CO₂ gridded emissions data product for tracer transport simulations and surface flux inversions. *Earth Syst. Sci. Data* **10**, 87–107 (2018).

53. Andres, R. J., Boden, T. A. & Higdon, D. M. Gridded uncertainty in fossil fuel carbon dioxide emission maps, a CDIAC example. *Atmos. Chem. Phys.* **16**, 14979–14995 (2016).

54. Harris, I., Osborn, T. J., Jones, P. & Lister, D. Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. *Sci. Data* **7**, 109, https://doi.org/10.1038/s41597-020-0453-3 (2020).

55. Spinoni, J. et al. Changes of heating and cooling degree-days in Europe from 1981 to 2100. *Int. J. Climatol.* **38**, e191–e208 (2018).

56. Spinoni, J., Vogt, J. & Barbosa, P. European degree-day climatologies and trends for the period 1951-2011. *Int. J. Climatol.* **35**, 25–36 (2015).

57. Andrew, R. M. Timely estimates of India's annual and monthly fossil CO₂ emissions. *Earth Syst. Sci. Data* **12**, 2411–2421 (2020).

58. EIA. Short-Term Energy Outlook. http://www.eia.gov/forecasts/steo/outlook.cfm (2019).

59. Olivier, J. G. J., Janssens-Maenhout, G., Muntean, M. & Peters, J. A. H.. Trends in global CO₂ emissions: 2016 report, JRC 103425. http://edgar.jrc.ec.europa.eu/overview.php?v=CO₂andGHG1970-2016 (2016).

60. Janssens-Maenhout, G. et al. EDGARv4.3.2_CO₂_FT2015 gridmaps. https://edgar.jrc.ec.europa.eu/overview.php?v=che_h2020. (2019).

61. Zheng, B. et al. Global atmospheric carbon monoxide budget 2000–2017 inferred from multi-species atmospheric inversions. *Earth Syst. Sci. Data* **11**, 1411–1436 (2019).

62. Saunois, M. et al. The Global Methane Budget 2000–2017. *Earth Syst. Sci. Data* **12**, 1561–1623 (2020).

63. Andres, R. J. et al. A synthesis of carbon dioxide emissions from fossil-fuel combustion. *Biogeosciences* **9**, 1845–1871 (2012).

64. Hogue, S., Marland, E., Andres, R. J., Marland, G. & Woodard, D. Uncertainty in gridded CO₂ emissions estimates Earth's Future. *Earth's Fu* **4**, 225–239 (2016).

65. Super, I., Dellaert, S. N. C., Visschedijk, A. J. H. & Van Der Gon, H. A. C. D. Uncertainty analysis of a European high-resolution emission inventory of CO₂ and CO to support inverse modelling and network design. *Atmos. Chem. Phys.* **20**, 1795–1816 (2020).

66. Oda, T. et al. Errors and uncertainties in a gridded carbon dioxide emissions inventory. *Mitig. Adapt. Strateg. Glob. Chang.* **24**, 1007–1050 (2019).

67. Han, P. et al. Evaluating China's fossil-fuel CO₂ emissions from a comprehensive dataset of nine inventories. *Atmos. Chem. Phys.* **4**, 1–21 (2020).

68. EIA. U.S. Energy Information Administration, Short-Term Energy Outlook, http://www.eia.gov/forecasts/steo/outlook.cfm (2020).

69. Keeling, R. F. Development of an interferometric analyzer for precise measurements of the atmospheric oxygen mole fraction. (Harvard University, 1988).

70. Jones, M. W. et al. Gridded fossil CO₂ emissions and related O₂ combustion consistent with national inventories 1959–2018. *Zenodo* https://doi.org/10.5281/zenodo.3958283 (2020).

71. Waldron, C. D., et al. 2006 IPCC Guidelines for National Greenhouse Gas Inventories, Volume 2: Energy in 2006 IPCC Guidelines for National Greenhouse Gas Inventories (ed. Eggleston, S., et al.) (Intergovernmental Panel on Climate Change, 2006).

72. Friedlingstein, P. et al. Global Carbon Budget 2020. *Earth Syst. Sci. Data* **12**, 3269–3340, https://doi.org/10.5194/essd-12-3269-2020 (2020).

73. Friedlingstein, P. et al. Supplemental data of the Global Carbon Budget 2020, available at: https://doi.org/10.18160/gcp-2020, ICOS-ERIC Carbon Portal, last accessed: 16 November 2020. https://doi.org/10.18160/gcp-2020 (2020).

74. Jones, M. W. et al. Gridded fossil CO₂ emissions and related O₂ combustion consistent with national inventories 1959–2019. *Zenodo* https://doi.org/10.5281/zenodo.4277267 (2020).
## Acknowledgements

The authors thank Pierre Friedlingstein, Christian Rödenbeck, Wouter Peters and Ingrid Luijkx for their helpful feedback on the manuscript. We further thank Wouter Peters for his part in securing the funding of the H2020 'CHE' project. We thank Penelope Pickers, Julia Steinbach and Nicolas Mayot for their insights regarding the available data on oxidative ratios. This work was principally funded by European Commission Horizon 2020 (H2020) 'CHE' project (grant no. 776186). M.W.J., F.C., P.C. and G.J.-M. were funded by the H2020 'CHE' project (grant no. 776186). G.P.P., R.M.A, P.C., F.C., G.J.-M. and M.W.J. were funded by the H2020 'VERIFY' project (no. 776810). A.J.D.-G., R.M.A. and G.P.P. were funded by the H2020 '4 C' project (no. 821003). G.P.P. was funded by the H2020 PARIS REINFORCE project (no. 820846). C.L.Q. was funded by the Royal Society (no. RP\R1\191063). We thank the UEA HPC team for support.

## Author contributions

M.W.J., C.L.Q., G.P.P. and R.M.A. designed the methodology with input from G.J.-M. G.J.-M. provided supplemental layers of the EDGAR dataset. R.M.A. and G.P.P. provided the GCB-NAE dataset. C.L.Q. and M.W.J. implemented the HCDD correction. M.W.J. completed all other formal analyses, technical validation and data visualisation and wrote the original draft of the article. M.W.J. wrote the code. A.J.D.-G. optimised the code. This work was made possible by prior funding from the Global Carbon Project. All authors contributed to review and editing of the article.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to M.W.J.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_15.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.

© The Author(s) 2021