CONTEXT: Crop models are essential tools for assessing the impact of climate change on national or regional agricultural production. Starting from meteorology, soil and crop management, fertilization and irrigation practices, they predict the yield of specific crop varieties. For long term assessments, climate models are the source of primary information. To make climate model results usable in a specific time frame context, bias adjustment (BA) is required. In fact, climate models tend to deviate from day-to-day values of the physical parameters while conserving the climate variability signal. BA brings the climatic signal to the actual values observed in a specific location and period, and to be representative of a specific period in absolute terms. BA techniques come in different flavours. The broadest categorization is univariate and multivariate methods. Multivariate methods adjust the variables considering possible cross -correlations while univariate methods treat the variables one by one without accounting for possible dependence on one another. OBJECTIVE: The hypothesis tested in this paper is that since crop models require as input climate variables that are in most of the cases cross -correlated, the multi-variate bias adjustment of the latter is likely to improve performance compared to univariate bias adjusted climate model results. METHODS: To verify this hypothesis, 14 BA methods were applied to 9 variables from 8 climate models at 21 locations across Europe and Northern Africa for a period of 5 years. Twelve crop models, from the AgMIP Wheat community, were run using the climate model results. All crop models, except one, were restarted at every growing season. The crop models were also run using the AgMERRA re -analysis. The latter were used as reference to compare the results when using the other climate models treated with the various sets of biasadjustment methods. RESULTS AND CONCLUSIONS: The results show that multivariate BA treatment should be preferred to univariate ones. The error obtained by comparing crop simulation obtained with AgMERRA with those obtained with multivariate bias -adjusted climate prediction is systematically lower. The error reduction varies as a function of the variable, the location, the crop model, and the climate model though the tendency is for smaller errors when multivariate methods are used to treat the latter. The results are attributed to the nature of crop models and the fact that multivariate methods consider more adequately the correlation existing between the meteorological variables. SIGNIFICANCE: The study shows the importance of considering the nature of a model and the selection of input data that best suited to the former. In this case the improvements produced when using multivariate data appears to be significant especially in the light of the variety of crop models used and the similar response obtained and it is therefore recommended.

Abstract. Quantification of land surface-atmosphere fluxes of carbon dioxide (CO2) fluxes and their trends and uncertainties is essential for monitoring progress of the EU27+UK bloc as it strives to meet ambitious targets determined by both international agreements and internal regulation. This study provides a consolidated synthesis of fossil sources (CO2 fossil) and natural sources and sinks over land (CO2 land) using bottom-up (BU) and top-down (TD) approaches for the European Union and United Kingdom (EU27+UK), updating earlier syntheses (Petrescu et al., 2020, 2021b). Given the wide scope of the work and the variety of approaches involved, this study aims to answer essential questions identified in the previous syntheses and understand the differences between datasets, particularly for poorly characterized fluxes from managed ecosystems. The work integrates updated emission inventory data, process-based model results, data-driven sectoral model results, and inverse modeling estimates, extending the previous period 1990–2018 to the year 2020 to the extent possible. BU and TD products are compared with European National Greenhouse Gas Inventories (NGHGIs) reported by Parties including the year 2019 under the United Nations Framework Convention on Climate Change (UNFCCC). The uncertainties of the EU27+UK NGHGI were evaluated using the standard deviation reported by the EU Member States following the guidelines of the Intergovernmental Panel on Climate Change (IPCC) and harmonized by gap-filling procedures. Variation in estimates produced with other methods, such as atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), originate from within-model uncertainty related to parameterization as well as structural differences between models. By comparing NGHGIs with other approaches, key sources of differences between estimates arise primarily in activities. System boundaries and emission categories create differences in CO2 fossil datasets, while different land use definitions for reporting emissions from Land Use, Land Use Change and Forestry (LULUCF) activities result in differences for CO2 land. The latter has important consequences for atmospheric inversions, leading to inversions reporting stronger sinks in vegetation and soils than are reported by the NGHGI. For CO2 fossil emissions, after harmonizing estimates based on common activities and selecting the most recent year available for all datasets, the UNFCCC NGHGI for the EU27+UK accounts for 3392 ± 49 Tg CO2 yr-1 (926 ± 13 Tg C yr-1), while eight other BU sources report a mean value of 3340 [3238,3401] [25th,75th percentile] Tg CO2 yr-1 (948 [937,961] Tg C yr-1). The sole top-down inversion of fossil emissions currently available accounts for 3800 Tg CO2 yr-1 (1038 Tg C yr-1), a value close to that of the NGHGI, but for which uncertainty estimates are not yet available. For the net CO2 land fluxes, during the most recent five-year period including the NGHGI estimates, the NGHGI accounted for -91 ± 32 Tg C yr-1 while six other BU approaches reported a mean sink of -62 [-117,-49] Tg C yr-1 and a 15-member ensemble of dynamic global vegetation models (DGVMs) reported -69 [-152,-5] Tg C yr-1. The five-year mean of three TD regional ensembles combined with one non-ensemble inversion of -73 Tg C yr-1 has a slightly smaller spread (0th–100th percentile of [-135,45] Tg C yr-1), and was calculated after removing land-atmosphere CO2 fluxes caused by lateral transport of carbon (crops, wood trade and inland waters) resulting in increased agreement with the the NGHGI and bottom-up approaches. Results at the sub-sector level (Forestland, Cropland, Grassland) show generally good agreement between the NGHGI and sub-sector-specific models, but results for a DGVM are mixed. Overall, for both CO2 fossil and net CO2 land fluxes, we find current independent approaches are consistent with the NGHGI at the scale of the EU27+UK. We conclude that CO2 emissions from fossil sources have decreased over the past 30 years in the EU27+UK, while large uncertainties on net uptake of CO2 by the land surface prevent trend identification. In addition, a gap on the order of 1000 Tg C yr-1 between CO2 fossil emissions and net CO2 uptake by the land exists regardless of the type of approach (NGHGI, TD, BU), falling well outside all available estimates of uncertainties. However, uncertainties in top-down approaches to estimate CO2 fossil emissions remain uncharacterized and are likely substantial. The data used to plot the figures are available at https://doi.org/10.5281/zenodo.7365863.

Abstract. This study, performed under the umbrella of the Task Force on Hemispheric Transport of Air Pollution (TF-HTAP), responds to the need of the global and regional atmospheric modelling community of having a mosaic emission inventory of air pollutants that conforms to specific requirements: global coverage, long time series, spatially distributed emissions with high time resolution, and a high sectoral resolution. The mosaic approach of integrating official regional emission inventories based on locally reported data, with a global inventory based on a globally consistent methodology, allows modellers to perform simulations of a high scientific quality while also ensuring that the results remain relevant to policymakers. HTAP_v3, an ad-hoc global mosaic of anthropogenic inventories, has been developed by integrating official inventories over specific areas (North America, Europe, Asia including Japan and Korea) with the independent Emissions Database for Global Atmospheric Research (EDGAR) inventory for the remaining world regions. The results are spatially and temporally distributed emissions of SO2, NOx, CO, NMVOC, NH3, PM10, PM2.5, Black Carbon (BC), and Organic Carbon (OC), with a spatial resolution of 0.1 x 0.1 degree and time intervals of months and years covering the period 2000–2018 (DOI 10.5281/zenodo.7516361, https://edgar.jrc.ec.europa.eu/dataset_htap_v3). The emissions are further disaggregated to 16 anthropogenic emitting sectors. This paper describes the methodology applied to develop such an emission mosaic, reports on source allocation, differences among existing inventories, and best practices for the mosaic compilation. One of the key strengths of the HTAP_v3 emission mosaic is its temporal coverage, enabling the analysis of emission trends over the past two decades. The development of a global emission mosaic over such long time series represents a unique product for global air quality modelling and for better-informed policy making, reflecting the community effort expended by the TF-HTAP to disentangle the complexity of transboundary transport of air pollution.

The food system, spanning from pre-production processes to post-production stages, is responsible for about one third of global greenhouse gas emissions and requires significant mitigation efforts to prevent dangerous levels of global warming. This article summarises trends and drivers of global food system emissions from 1990 to 2018. We highlight regional diversity in patterns of food system emissions and identify the highest global emitters. While food system emissions have stabilised in some regions and countries, global emissions are increasing, with growth in certain sectors and countries outweighing the handful of cases where sustained emissions reductions have been realised. Emissions from livestock rearing account for a large portion of global emissions, and the contribution of post-production emissions is steadily increasing in all regions. We also provide an overview of food system policies at the national level, mapping them to each emissions segment. This highlights the significant shortfall in policy activity required to address the challenge of climate change mitigation in general, and the impacts of livestock and post-production emissions in particular. Our work lays the groundwork for addressing specific country-level questions on optimal policy pathways to achieve emission reductions.

Abstract. This study, performed under the umbrella of the Task Force on Hemispheric Transport of Air Pollution (TF-HTAP), responds to the global and regional atmospheric modelling community's need of a mosaic emission inventory of air pollutants that conforms to specific requirements: global coverage, long time series, spatially distributed emissions with high time resolution, and a high sectoral resolution. The mosaic approach of integrating official regional emission inventories based on locally reported data, with a global inventory based on a globally consistent methodology, allows modellers to perform simulations of high scientific quality while also ensuring that the results remain relevant to policymakers. HTAP_v3, an ad hoc global mosaic of anthropogenic inventories, has been developed by integrating official inventories over specific areas (North America, Europe, Asia including Japan and South Korea) with the independent Emissions Database for Global Atmospheric Research (EDGAR) inventory for the remaining world regions. The results are spatially and temporally distributed emissions of SO2, NOx, CO, non-methane volatile organic compounds (NMVOCs), NH3, PM10, PM2.5, black carbon (BC), and organic carbon (OC), with a spatial resolution of 0.1∘ × 0.1∘ and time intervals of months and years, covering the period 2000–2018 (https://doi.org/10.5281/zenodo.7516361, Crippa, 2023, https://edgar.jrc.ec.europa.eu/dataset_htap_v3, last access: June 2023). The emissions are further disaggregated into 16 anthropogenic emitting sectors. This paper describes the methodology applied to develop such an emission mosaic, reports on source allocation, differences among existing inventories, and best practices for the mosaic compilation. One of the key strengths of the HTAP_v3 emission mosaic is its temporal coverage, enabling the analysis of emission trends over the past 2 decades. The development of a global emission mosaic over such long time series represents a unique product for global air quality modelling and for better-informed policymaking, reflecting the community effort expended by the TF-HTAP to disentangle the complexity of transboundary transport of air pollution.

Knowledge of the spatial distribution of the fluxes of greenhouse gases (GHGs) and their temporal variability as well as flux attribution to natural and anthropogenic processes is essential to monitoring the progress in mitigating anthropogenic emissions under the Paris Agreement and to inform its global stocktake. This study provides a consolidated synthesis of CH4 and N2O emissions using bottom-up (BU) and top-down (TD) approaches for the European Union and UK (EU27 + UK) and updates earlier syntheses (Petrescu et al., 2020, 2021). The work integrates updated emission inventory data, process-based model results, data-driven sector model results and inverse modeling estimates, and it extends the previous period of 1990–2017 to 2019. BU and TD products are compared with European national greenhouse gas inventories (NGHGIs) reported by parties under the United Nations Framework Convention on Climate Change (UNFCCC) in 2021. Uncertainties in NGHGIs, as reported to the UNFCCC by the EU and its member states, are also included in the synthesis. Variations in estimates produced with other methods, such as atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), arise from diverse sources including within-model uncertainty related to parameterization as well as structural differences between models. By comparing NGHGIs with other approaches, the activities included are a key source of bias between estimates, e.g., anthropogenic and natural fluxes, which in atmospheric inversions are sensitive to the prior geospatial distribution of emissions. For CH4 emissions, over the updated 2015–2019 period, which covers a sufficiently robust number of overlapping estimates, and most importantly the NGHGIs, the anthropogenic BU approaches are directly comparable, accounting for mean emissions of 20.5 Tg CH4 yr−1 (EDGARv6.0, last year 2018) and 18.4 Tg CH4 yr−1 (GAINS, last year 2015), close to the NGHGI estimates of 17.5±2.1 Tg CH4 yr−1. TD inversion estimates give higher emission estimates, as they also detect natural emissions. Over the same period, high-resolution regional TD inversions report a mean emission of 34 Tg CH4 yr−1. Coarser-resolution global-scale TD inversions result in emission estimates of 23 and 24 Tg CH4 yr−1 inferred from GOSAT and surface (SURF) network atmospheric measurements, respectively. The magnitude of natural peatland and mineral soil emissions from the JSBACH–HIMMELI model, natural rivers, lake and reservoir emissions, geological sources, and biomass burning together could account for the gap between NGHGI and inversions and account for 8 Tg CH4 yr−1. For N2O emissions, over the 2015–2019 period, both BU products (EDGARv6.0 and GAINS) report a mean value of anthropogenic emissions of 0.9 Tg N2O yr−1, close to the NGHGI data (0.8±55 % Tg N2O yr−1). Over the same period, the mean of TD global and regional inversions was 1.4 Tg N2O yr−1 (excluding TOMCAT, which reported no data). The TD and BU comparison method defined in this study can be operationalized for future annual updates for the calculation of CH4 and N2O budgets at the national and EU27 + UK scales. Future comparability will be enhanced with further steps involving analysis at finer temporal resolutions and estimation of emissions over intra-annual timescales, which is of great importance for CH4 and N2O, and may help identify sector contributions to divergence between prior and posterior estimates at the annual and/or inter-annual scale. Even if currently comparison between CH4 and N2O inversion estimates and NGHGIs is highly uncertain because of the large spread in the inversion results, TD inversions inferred from atmospheric observations represent the most independent data against which inventory totals can be compared. With anticipated improvements in atmospheric modeling and observations, as well as modeling of natural fluxes, TD inversions may arguably emerge as the most powerful tool for verifying emission inventories for CH4, N2O and other GHGs. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.7553800 (Petrescu et al., 2023).

We present a near-real-time global gridded daily CO 2 emissions dataset (GRACED) throughout 2021. GRACED provides gridded CO 2 emissions at a 0.1° × 0.1° spatial resolution and 1-day temporal resolution from cement production and fossil fuel combustion over seven sectors, including industry, power, residential consumption, ground transportation, international aviation, domestic aviation, and international shipping. GRACED is prepared from the near-real-time daily national CO 2 emissions estimates (Carbon Monitor), multi-source spatial activity data emissions and satellite NO 2 data for time variations of those spatial activity data. GRACED provides the most timely overview of emissions distribution changes, which enables more accurate and timely identification of when and where fossil CO 2 emissions have rebounded and decreased. Uncertainty analysis of GRACED gives a grid-level two-sigma uncertainty of value of ±19.9% in 2021, indicating the reliability of GRACED was not sacrificed for the sake of higher spatiotemporal resolution that GRACED provides. Continuing to update GRACED in a timely manner could help policymakers monitor energy and climate policies’ effectiveness and make adjustments quickly.

Timely, fine-grained gridded carbon emission datasets are particularly important for global climate change research. Often, fine-grained datasets are challenging to visualize over the globe, and clear visualization tools are also needed. Therefore, we present a near-real-time global gridded daily CO2 emissions dataset (GRACED). GRACED provides gridded CO2 emissions at a 0.1° × 0.1° spatial resolution and 1-day temporal resolution from cement production and fossil fuel combustion over seven sectors, including power, industry, residential consumption, ground transportation, domestic aviation, international aviation, and international shipping. GRACED is prepared from the near-real-time daily national CO2 emissions estimates (Carbon Monitor), multi-source spatial activity data and satellite NO2 data for time variations of those spatial activity data. Here, we examined the spatial patterns of sectoral CO2 emission changes from January 1, 2019, to December 31, 2021. In 2021, most regions showed rapid rebounds in carbon emissions compared with 2020, reflecting the continuing challenges to accelerate climate mitigation in the post-COVID era. GRACED provides the most timely and more refined overview than any other previously published datasets, which enables more accurate and timely identification of when and where fossil CO2 emissions have rebounded and decreased as the world recovers from COVID-19 and witnesses contrasted efforts to decarbonize energy systems. Uncertainty analysis of GRACED gives a grid-level two-sigma uncertainty of value of ±19.9%, indicating the reliability of GRACED was not sacrificed for the sake of higher spatiotemporal resolution that GRACED provides. In addition, we also examined the distribution of emission in a grid-wise perspective for major emission datasets, and compared it with GRACED. The similarity in emission distribution was observed in GRACED and other datasets. One of the advantages of our dataset is that it provides worldwide near-real-time monitoring of CO2 emissions with different fine spatial scales at the sub-national level, such as cities, thus enhancing our comprehension of spatial and temporal changes in CO2 emissions and anthropogenic activities. With the continued extension of GRACED time series, we present crucial daily-level input to analyze CO2 emission changes in the post-COVID era, which will ultimately facilitate and aid in designing more localized and adaptive management policies for the purpose of climate change mitigation in the post-COVID era.

Quantification of land surface–atmosphere fluxes of carbon dioxide (CO2) and their trends and uncertainties is essential for monitoring progress of the EU27+UK bloc as it strives to meet ambitious targets determined by both international agreements and internal regulation. This study provides a consolidated synthesis of fossil sources (CO2 fossil) and natural (including formally managed ecosystems) sources and sinks over land (CO2 land) using bottom-up (BU) and top-down (TD) approaches for the European Union and United Kingdom (EU27+UK), updating earlier syntheses (Petrescu et al., 2020, 2021). Given the wide scope of the work and the variety of approaches involved, this study aims to answer essential questions identified in the previous syntheses and understand the differences between datasets, particularly for poorly characterized fluxes from managed and unmanaged ecosystems. The work integrates updated emission inventory data, process-based model results, data-driven categorical model results, and inverse modeling estimates, extending the previous period 1990–2018 to the year 2020 to the extent possible. BU and TD products are compared with the European national greenhouse gas inventory (NGHGI) reported by parties including the year 2019 under the United Nations Framework Convention on Climate Change (UNFCCC). The uncertainties of the EU27+UK NGHGI were evaluated using the standard deviation reported by the EU member states following the guidelines of the Intergovernmental Panel on Climate Change (IPCC) and harmonized by gap-filling procedures. Variation in estimates produced with other methods, such as atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), originate from within-model uncertainty related to parameterization as well as structural differences between models. By comparing the NGHGI with other approaches, key sources of differences between estimates arise primarily in activities. System boundaries and emission categories create differences in CO2 fossil datasets, while different land use definitions for reporting emissions from land use, land use change, and forestry (LULUCF) activities result in differences for CO2 land. The latter has important consequences for atmospheric inversions, leading to inversions reporting stronger sinks in vegetation and soils than are reported by the NGHGI. For CO2 fossil emissions, after harmonizing estimates based on common activities and selecting the most recent year available for all datasets, the UNFCCC NGHGI for the EU27+UK accounts for 926 ± 13 Tg C yr−1, while eight other BU sources report a mean value of 948 [937,961] Tg C yr−1 (25th, 75th percentiles). The sole top-down inversion of fossil emissions currently available accounts for 875 Tg C in this same year, a value outside the uncertainty of both the NGHGI and bottom-up ensemble estimates and for which uncertainty estimates are not currently available. For the net CO2 land fluxes, during the most recent 5-year period including the NGHGI estimates, the NGHGI accounted for −91 ± 32 Tg C yr−1, while six other BU approaches reported a mean sink of −62 [-117,-49] Tg C yr−1, and a 15-member ensemble of dynamic global vegetation models (DGVMs) reported −69 [-152,-5] Tg C yr−1. The 5-year mean of three TD regional ensembles combined with one non-ensemble inversion of −73 Tg C yr−1 has a slightly smaller spread (0th–100th percentiles of [-135,+45] Tg C yr−1), and it was calculated after removing net land–atmosphere CO2 fluxes caused by lateral transport of carbon (crop trade, wood trade, river transport, and net uptake from inland water bodies), resulting in increased agreement with the NGHGI and bottom-up approaches. Results at the category level (Forest Land, Cropland, Grassland) generally show good agreement between the NGHGI and category-specific models, but results for DGVMs are mixed. Overall, for both CO2 fossil and net CO2 land fluxes, we find that current independent approaches are consistent with the NGHGI at the scale of the EU27+UK. We conclude that CO2 emissions from fossil sources have decreased over the past 30 years in the EU27+UK, while land fluxes are relatively stable: positive or negative trends larger (smaller) than 0.07 (−0.61) Tg C yr−2 can be ruled out for the NGHGI. In addition, a gap on the order of 1000 Tg C yr−1 between CO2 fossil emissions and net CO2 uptake by the land exists regardless of the type of approach (NGHGI, TD, BU), falling well outside all available estimates of uncertainties. However, uncertainties in top-down approaches to estimate CO2 fossil emissions remain uncharacterized and are likely substantial, in addition to known uncertainties in top-down estimates of the land fluxes. The data used to plot the figures are available at https://doi.org/10.5281/zenodo.8148461 (McGrath et al., 2023).

In 2015 and 2016, atmospheric transport modeling challenges were conducted in the context of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) verification, however, with a more limited scope with respect to emission inventories, simulation period and number of relevant samples (i.e., those above the Minimum Detectable Concentration (MDC)) involved. Therefore, a more comprehensive atmospheric transport modeling challenge was organized in 2019. Stack release data of Xe-133 were provided by the Institut National des Radioéléments/IRE (Belgium) and the Canadian Nuclear Laboratories/CNL (Canada) and accounted for in the simulations over a three (mandatory) or six (optional) months period. Best estimate emissions of additional facilities (radiopharmaceutical production and nuclear research facilities, commercial reactors or relevant research reactors) of the Northern Hemisphere were included as well. Model results were compared with observed atmospheric activity concentrations at four International Monitoring System (IMS) stations located in Europe and North America with overall considerable influence of IRE and/or CNL emissions for evaluation of the participants' runs. Participants were prompted to work with controlled and harmonized model set-ups to make runs more comparable, but also to increase diversity. It was found that using the stack emissions of IRE and CNL with daily resolution does not lead to better results than disaggregating annual emissions of these two facilities taken from the literature if an overall score for all stations covering all valid observed samples is considered. A moderate benefit of roughly 10% is visible in statistical scores for samples influenced by IRE and/or CNL to at least 50% and there can be considerable benefit for individual samples. Effects of transport errors, not properly characterized remaining emitters and long IMS sampling times (12-24 h) undoubtedly are in contrast to and reduce the benefit of high-quality IRE and CNL stack data. Complementary best estimates for remaining emitters push the scores up by 18% compared to just considering IRE and CNL emissions alone. Despite the efforts undertaken the full multi-model ensemble built is highly redundant. An ensemble based on a few arbitrary runs is sufficient to model the Xe-133 background at the stations investigated. The effective ensemble size is below five. An optimized ensemble at each station has on average slightly higher skill compared to the full ensemble. However, the improvement (maximum of 20% and minimum of 3% in RMSE) in skill is likely being too small for being exploited for an independent period.

Atmospheric nitrous oxide (N2O) contributes to global warming and stratospheric ozone depletion, so reducing uncertainty in estimates of emissions from different sources is important for climate policy. In this study, we simulate atmospheric N2O using an atmospheric chemistry-transport model (ACTM), and the results are first compared with the in situ measurements. Five combinations of known (a priori) N2O emissions due to natural soil, agricultural land, other human activities, and sea–air exchange are used. The N2O lifetime is 127.6 ± 4.0 yr in the control ACTM simulation (range indicates interannual variability). Regional N2O emissions are optimized using Bayesian inverse modeling for 84 partitions of the globe at monthly intervals, using measurements at 42 sites around the world covering 1997–2019. The best estimated global land and ocean emissions are 12.99 ± 0.22 TgN yr−1 and 2.74 ± 0.27 TgN yr−1, respectively, for 2000–2009, and 14.30 ± 0.20 TgN yr−1 and 2.91 ± 0.27 TgN yr−1, respectively, for 2010–2019. On regional scales, we find that the most recent ocean emission estimation, with lower emissions in the Southern Ocean regions, fits better with that predicted by the inversions. Marginally higher (lower) emissions than the inventory/model for the tropical (extratropical) land regions are estimated and validated using independent aircraft observations. Global land and ocean emission variabilities show a statistically significant correlation with El Niño Southern Oscillation (ENSO). Analysis of regional land emissions shows increases over America (Temperate North, Central, and Tropical), Central Africa, and Asia (South, East, and Southeast) between the 2000s and 2010s. Only Europe as a whole recorded a slight decrease in N2O emissions due to the chemical industry. Our inversions suggest revisions to seasonal emission variations for three of the 15 land regions (East Asia, Temperate North America, and Central Africa), and the Southern Ocean region. The terrestrial ecosystem model (Vegetation Integrative SImulator for Trace Gases) can simulate annual total emissions in agreement with the observed N2O growth rate since 1978, but the lag-time scales of N2O emissions from nitrogen fertilizer application may need to be revised.

Abstract. Knowledge of the spatial distribution of the fluxes of greenhouse gases and their temporal variability as well as flux attribution to natural and anthropogenic processes is essential to monitoring the progress in mitigating anthropogenic emissions under the Paris Agreement and to inform its Global Stocktake. This study provides a consolidated synthesis of CH4 and N2O emissions using bottom-up (BU) and top-down (TD) approaches for the European Union and UK (EU27+UK) and updates earlier syntheses (Petrescu et al., 2020, 2021). The work integrates updated emission inventory data, process-based model results, data-driven sector model results, inverse modelling estimates, and extends the previous period 1990–2017 to 2020. BU and TD products are compared with European National GHG Inventories (NGHGI) reported by Parties under the United Nations Framework Convention on Climate Change (UNFCCC) in 2021. The uncertainties of NGHGIs were evaluated using the standard deviation obtained by varying parameters of inventory calculations, reported by the EU Member States following the guidelines of the Intergovernmental Panel on Climate Change (IPCC) and harmonized by gap-filling procedures. Variation in estimates produced with other methods, such as atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), arise from diverse sources including within-model uncertainty related to parameterization as well as structural differences between models. By comparing NGHGIs with other approaches, the activities included are a key source of bias between estimates e.g. anthropogenic and natural fluxes, which, in atmospheric inversions are sensitive to the prior geospatial distribution of emissions. For CH4 emissions, over the updated 2015–2019 period, which covers a sufficiently robust number of overlapping estimates, and most importantly the NGHGIs, the anthropogenic BU approaches are directly comparable, accounting for mean emissions of 20.5 Tg CH4 yr−1 (EDGAR v5v6.0, last year 2018) and 18.4 Tg CH4 yr−1 (GAINS, 2015), close to the NGHGI estimates of 17.5 ± 2.1 Tg CH4 yr−1. TD inversions estimates give higher emission estimates, as they also detect natural emissions. Over the same period, high resolution regional TD inversions report a mean emission of 34 Tg CH4 yr−1. Coarser-resolution global-scale TD inversions result in emission estimates of 23 Tg CH4 yr−1 and 24 Tg CH4 yr−1 inferred from GOSAT and surface (SURF) network atmospheric measurements, respectively. The magnitude of natural peatland and mineral soils emissions from the JSBACH-HIMMELI model, natural rivers, lakes and reservoirs emissions, geological sources and biomass burning together could account for the gap between NGHGI and inversions and account for 8 Tg CH4 yr−1. For N2O emissions, over the 2015–2019 period, both BU products (EDGAR v5v6.0 and GAINS) report a mean value of anthropogenic emissions of 0.9 Tg N2O yr−1, close to the NGHGI data (0.8 ± 55 % Tg N2O yr−1). Over the same period, the mean of TD global and regional inversions was 1.4 Tg N2O yr−1 (excluding TOMCAT which reported no data). The TD and BU comparison method defined in this study can be "operationalized" for future annual updates for the calculation of CH4 and N2O budgets at the national and EU27+UK scales. Future comparability will be enhanced with further steps involving analysis at finer temporal resolutions and estimation of emissions over intra-annual timescales, of great importance for CH4 and N2O, which may help identify sector contributions to divergence between prior and posterior estimates at the annual/inter-annual scale. Even if currently comparison between CH4 and N2O inversions estimates and NGHGIs is highly uncertain because of the large spread in the inversion results, TD inversions inferred from atmospheric observations represent the most independent data against which inventory totals can be compared. With anticipated improvements in atmospheric modelling and observations, as well as modelling of natural fluxes, TD inversions may arguably emerge as the most powerful tool for verifying emissions inventories for CH4, N2O and other GHGs. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.6992472 (Petrescu et al., 2022).

Food systems are important contributors to global emissions of air pollutants. Here, building on the EDGAR-FOOD database of greenhouse gas emissions, we estimate major air pollutant compounds emitted by different stages of the food system, at country level, during the past 50 years, resulting from food production, processing, packaging, transport, retail, consumption and disposal. Air pollutant estimates from food systems include total nitrogen and its components (N2O, NH3 and NOx), SO2, CO, non-methane volatile organic compounds (NMVOC) and particulate matter (PM10, PM2.5, black carbon and organic carbon). We show that 10% to 90% of air pollutant emissions come from food systems, resulting from steady increases over the past five decades. In 2018, more than half of total N (and 87% of ammonia) emissions come from food systems and up to 35% of particulate matter. Food system emissions are responsible for about 22.4% of global mortality due to poor air quality and 1.4% of global crop production losses. The exact contribution of food systems to air pollution is unknown. On the basis of the European Commission’s EDGAR-FOOD database, a global emission inventory of air pollutants from the food systems, this study quantifies historic emissions of major pollutant compounds at each stage of the food supply chain at country level.

The effect of the 2018 extreme meteorological conditions in Europe on methane (CH 4 ) emissions is examined using estimates from four atmospheric inversions calculated for the period 2005–2018. For most of Europe, we find no anomaly in 2018 compared to the 2005–2018 mean. However, we find a positive anomaly for the Netherlands in April, which coincided with positive temperature and soil moisture anomalies suggesting an increase in biogenic sources. We also find a negative anomaly for the Netherlands for September–October, which coincided with a negative anomaly in soil moisture, suggesting a decrease in soil sources. In addition, we find a positive anomaly for Serbia in spring, summer and autumn, which coincided with increases in temperature and soil moisture, again suggestive of changes in biogenic sources, and the annual emission for 2018 was 33 ± 38% higher than the 2005–2017 mean. These results indicate that CH 4 emissions from areas where the natural source is thought to be relatively small can still vary due to meteorological conditions. At the European scale though, the degree of variability over 2005–2018 was small, and there was negligible impact on the annual CH 4 emissions in 2018 despite the extreme meteorological conditions. This article is part of a discussion meeting issue ‘Rising methane: is warming feeding warming? (part 2)’.

Comprehensive and reliable information on anthropogenic sources of greenhouse gas emissions is required to track progress towards keeping warming well below 2°C as agreed upon in the Paris Agreement. Here we provide a dataset on anthropogenic GHG emissions 1970-2019 with a broad country and sector coverage. We build the dataset from recent releases from the “Emissions Database for Global Atmospheric Research” (EDGAR) for CO2 emissions from fossil fuel combustion and industry (FFI), CH4 emissions, N2O emissions, and fluorinated gases and use a well-established fast-track method to extend this dataset from 2018 to 2019. We complement this with information on net CO2 emissions from land use, land-use change and forestry (LULUCF) from three available bookkeeping models.

Precise and high-resolution carbon dioxide (CO2) emission data is of great importance in achieving carbon neutrality around the world. Here we present for the first time the near-real-time Global Gridded Daily CO2 Emissions Dataset (GRACED) from fossil fuel and cement production with a global spatial resolution of 0.1° by 0.1° and a temporal resolution of 1 day. Gridded fossil emissions are computed for different sectors based on the daily national CO2 emissions from near-real-time dataset (Carbon Monitor), the spatial patterns of point source emission dataset Global Energy Infrastructure Emissions Database (GID), Emission Database for Global Atmospheric Research (EDGAR), and spatiotemporal patters of satellite nitrogen dioxide (NO2) retrievals. Our study on the global CO2 emissions responds to the growing and urgent need for high-quality, fine-grained, near-real-time CO2 emissions estimates to support global emissions monitoring across various spatial scales. We show the spatial patterns of emission changes for power, industry, residential consumption, ground transportation, domestic and international aviation, and international shipping sectors from January 1, 2019, to December 31, 2020. This gives thorough insights into the relative contributions from each sector. Furthermore, it provides the most up-to-date and fine-grained overview of where and when fossil CO2 emissions have decreased and rebounded in response to emergencies (e.g., coronavirus disease 2019 [COVID-19]) and other disturbances of human activities of any previously published dataset. As the world recovers from the pandemic and decarbonizes its energy systems, regular updates of this dataset will enable policymakers to more closely monitor the effectiveness of climate and energy policies and quickly adapt.