Success of the Paris Agreement relies on rapid reductions in fossil fuel CO2 (ffCO2) emissions. Atmospheric data can verify the ffCO2 reductions pledged by nations in their nationally determined contributions. However, estimating ffCO2 from atmospheric CO2 is challenging due to natural fluxes and varying backgrounds. One approach is to combine with nitrogen oxides (NOx = NO + NO2), which are co-emitted with CO2 during combustion. A key challenge in using NOx to estimate ffCO2 is the computational cost of modelling atmospheric photochemistry. Additionally, the NO2 : NO column ratio must be well understood to convert model NOx columns to NO2 columns for comparison with satellite data. We use random forest regression to parameterise NOx chemistry, relying only on meteorological parameters and NOx concentration. The regression is trained on outputs from a nested GEOS (Goddard Earth Observing System)-Chem model simulation for mainland Europe in 2019. We develop a monthly NOx chemistry parameterisation that performs well when tested on perturbed emission runs (R2 > 0.95) and on unseen meteorology for 2021 (R2 > 0.79). We also parameterise the NO2 : NO ratio (R2 > 0.99 on perturbed outputs, R2 > 0.92 on unseen meteorology). Additionally, we present an alternative method to predict NOx rates by scaling baseline NOx rates with changes in NOx concentration (R2 = 1.0 on perturbed outputs). Our models reproduce NO2 columns with minimal deviation from full-chemistry models, with reconstruction error smaller than the TROPOspheric Monitoring Instrument (TROPOMI) precision in over 99.9 % of cases, supporting robust ffCO2 inversion efforts. These results provide a robust framework for accurately estimating fossil fuel CO2 emissions from atmospheric data, enabling more reliable monitoring and verification of global emissions reductions.

The Copernicus Atmosphere Monitoring Service (CAMS) has recently produced a greenhouse gas reanalysis (version egg4) that covers almost 2 decades from 2003 to 2020 and which will be extended in the future. This reanalysis dataset includes carbon dioxide (CO2) and methane (CH4). The reanalysis procedure combines model data with satellite data into a globally complete and consistent dataset using the European Centre for Medium-Range Weather Forecasts' Integrated Forecasting System (IFS). This dataset has been carefully evaluated against independent observations to ensure validity and to point out deficiencies to the user. The greenhouse gas reanalysis can be used to examine the impact of atmospheric greenhouse gas concentrations on climate change (such as global and regional climate radiative forcing), assess intercontinental transport, and serve as boundary conditions for regional simulations, among other applications and scientific uses. The caveats associated with changes in assimilated observations and fixed underlying emissions are highlighted, as is their impact on the estimation of trends and annual growth rates of these long-lived greenhouse gases.

Abstract The impact of urbanization on local weather patterns affects over half the global population. Global numerical weather prediction systems have reached a resolution at which urban conurbations can be spatially resolved, justifying their representation within land surface parameterizations with the aim of improving local predictions. Additionally, real‐time atmospheric monitoring of trace gas emissions can utilize weather variables relevant for urban areas. We investigated whether a simple single‐layer urban canopy scheme can be used within a global forecast model to jointly improve predictions of near‐surface weather variables and residential CO2 emissions. The scheme has been implemented in the Integrated Forecast System used operationally at the European Centre for Medium‐Range Weather Forecasts running at ∼9 km horizontal resolution. First, we selected a suitable urban land cover map (ECOCLIMAP‐SG) based on comparisons with regional data and land surface temperature MODIS retrievals. The urban scheme is verified by providing improved 2 m temperature (∼10%) and 10 m wind (∼17%) RMSE values for both summer and winter months around urban environments. The influence of the scheme was most noticeable at night. Additionally, we have implemented a simple temperature‐dependent residential emissions model to calculate real‐time CO2 heating emissions. These were validated against existing offline products, national reporting and by comparing atmospheric simulations with total column CO2 observations. The results show an improved temporal variability of emissions, which arise from synoptic scale temperature changes. Given the improved predictability from the urban scheme for both weather and emissions, it will be operationally implemented in an upcoming model cycle.

The Copernicus Atmosphere Monitoring Service has recently produced a greenhouse gases reanalysis (version egg4) that covers almost two decades from 2003 to 2020 and will be extended in the future. This reanalysis dataset includes carbon dioxide (CO2) and methane (CH4). The reanalysis procedure combines model data with satellite data into a globally complete and consistent dataset using the European Centre for Medium-range Weather Forecasts’ Integrated Forecasting System (IFS). This dataset has been carefully evaluated against independent observations to ensure validity and point out deficiencies to the user. The greenhouse gas reanalysis can be used to examine the impact of atmospheric greenhouse gases concentrations on climate change, such as global and regional climate radiative forcing, assess intercontinental transport, and also serve as boundary conditions for regional simulations, among other applications and scientific studies. The caveats associated with changes in assimilated observations and fixed underlying emissions are highlighted, as well as their impact on the estimation of trends and annual growth rates of these long-lived greenhouse gases.

The CO2 Human Emissions project has generated realistic high-resolution 9 km global simulations for atmospheric carbon tracers referred to as nature runs to foster carbon-cycle research applications with current and planned satellite missions, as well as the surge of in situ observations. Realistic atmospheric CO2, CH4 and CO fields can provide a reference for assessing the impact of proposed designs of new satellites and in situ networks and to study atmospheric variability of the tracers modulated by the weather. The simulations spanning 2015 are based on the Copernicus Atmosphere Monitoring Service forecasts at the European Centre for Medium Range Weather Forecasts, with improvements in various model components and input data such as anthropogenic emissions, in preparation of a CO2 Monitoring and Verification Support system. The relative contribution of different emissions and natural fluxes towards observed atmospheric variability is diagnosed by additional tagged tracers in the simulations. The evaluation of such high-resolution model simulations can be used to identify model deficiencies and guide further model improvements.

Abstract. Concentrations of atmospheric methane (CH4), the second most important greenhouse gas, continue to grow. In recent years this growth rate has increased further (2020: +14.7 ppb), the cause of which remains largely unknown. Here, we demonstrate a high-resolution (~80 km), short-window (24-hour) 4D-Var global inversion system based on the ECMWF Integrated Forecasting System (IFS) and newly available satellite observations. The largest national disagreement found between prior (63.1 Tg yr−1) and posterior (59.8 Tg yr−1) CH4 emissions is from China, mainly attributed to the energy sector. Emissions estimated form our global system agree well with previous basin-wide regional studies and point source specific studies. Emission events (leaks/blowouts) >10 t hr−1 were detected, but without accurate prior uncertainty information, were not well quantified. Our results suggest that global anthropogenic CH4 emissions for 2020 were 5.7 Tg yr−1 (+1.6 %) higher than for 2019, mainly attributed to the energy and agricultural sectors. Regionally, the largest 2020 increases were seen from China (+2.6 Tg yr−1, 4.3 %), with smaller increases from India (+0.8 Tg yr−1, 2.2 %) and Indonesia (+0.3 Tg yr−1, 2.6 %). Results show the rise in emissions, and subsequent atmospheric growth, would have occurred with or without the COVID-19 slowdown. During the onset of the global slowdown (March–April, 2020) energy sector CH4 emissions from China increased; however, during later months (May–June, 2020) emissions decreased below expected pre-slowdown levels. The accumulated impact of the slowdown on CH4 emissions from March–June 2020 is found to be small. Changes in atmospheric chemistry, not investigated here, may have contributed to the observed growth in 2020. Future work aims to develop the global IFS inversion system and to extend the 4D-Var window-length using a hybrid ensemble-variational method.

The Copernicus Atmosphere Monitoring Service (CAMS), part of the European Union's Earth observation program Copernicus, entered operations in July 2015. Implemented by the European Centre for Medium-Range Weather Forecasts (ECMWF) as a truly European effort with over 23,500 direct data users and well over 200 million end users worldwide as of March 2022, CAMS delivers numerous global and regional information products about air quality, inventory-based emissions and observation-based surface fluxes of greenhouse gases and from biomass burning, solar energy, ozone and UV radiation, and climate forcings. Access to CAMS products is open and free of charge via the Atmosphere Data Store. The CAMS global atmospheric composition analyses, forecasts, and reanalyses build on ECMWF's Integrated Forecasting System (IFS) and exploit over 90 different satellite data streams. The global products are complemented by coherent higher-resolution regional air quality products over Europe derived from multisystem analyses and forecasts. CAMS information products also include policy support such as quantitative impact assessment of short- and long-term pollutant-emission mitigation scenarios, source apportionment information, and annual European air quality assessment reports. Relevant CAMS products are cited and used for instance in IPCC Assessment Reports. Providing dedicated support for users operating smartphone applications, websites, or TV bulletins in Europe and worldwide is also integral to the service. This paper presents key achievements of the CAMS initial phase (2014-21) and outlines some of its new components for the second phase (2021-28), e.g., the new Copernicus anthropogenic CO2 emissions Monitoring and Verification Support capacity that will monitor global anthropogenic emissions of key greenhouse gases.

The Copernicus Atmosphere Monitoring Service (CAMS) provides near-real-time forecast and reanalysis of aerosols using the ECMWF Integrated Forecasting System with atmospheric composition extension, constrained by the assimilation of MODIS and the Polar Multi-Sensor Aerosol Optical Properties (PMAp) aerosol optical depth (AOD). The objective of this work is to evaluate two new near-real-time AOD products to prepare for their assimilation into CAMS, namely the Copernicus AOD (collection 1) from the Sea and Land Surface Temperature Radiometer (SLSTR) on board Sentinel 3-A/B over ocean and the NOAA EPS AOD (v2.r1) from VIIRS on board S-NPP and NOAA20 over both land and ocean. The differences between MODIS (C6.1), PMAp (v2.1), VIIRS (v2.r1), and SLSTR (C1) AOD as well as their departure from the modeled AOD were assessed at the model grid resolution (i.e., level-3) using the 3-month AOD average (December 2019–February 2020 and March–May 2020). VIIRS and MODIS show the best consistency across the products, which is explained by instrument and retrieval algorithm similarities. VIIRS AOD is frequently lower over the ocean background and higher over biomass burning and dust source land regions compared to MODIS. VIIRS shows larger spatial coverage over land and resolves finer spatial structures such as the transport of Australian biomass burning smoke over the Pacific, which can be explained by the use of a heavy aerosol detection test in the retrieval algorithm. Our results confirm the positive offset over ocean (i) between Terra/MODIS and Aqua/MODIS due to the non-corrected radiometric calibration degradation of Terra/MODIS in the Dark Target algorithm and (ii) between SNPP/VIIRS and NOAA20/VIIRS due to the positive bias in the solar reflective bands of SNPP/VIIRS. SLSTR AOD shows much smaller level-3 values than the rest of the products, which is mainly related to differences in spatial representativity at the IFS grid spatial resolution due to the stringent cloud filtering applied to the SLSTR radiances. Finally, the geometry characteristics of the instrument, which drive the range of scattering angles sampled by the instrument, can explain a large part of the differences between retrievals such as the positive offset between PMAp datasets from MetOp-B and MetOp-A.

The societal benefits of numerical weather prediction (NWP) forecasts are most evident in populated areas. An urban representation within NWP models should provide improved forecast accuracy. Here, we present the preliminary implementation of an urban scheme within the Integrated Forecasting System (IFS) using a simplified single‐layer urban canopy model. The scheme makes assumptions of canyon geometry and considers fluxes from roads, walls, and roofs. Temperature observations were used to optimize single‐column model (SCM) parameters using the Gauss‐Newton method. Observation comparisons over six European cities, show a 2‐m temperature root‐mean‐squared error reduction from 1.85 to 1.75 K with the urban scheme. Optimized parameters were used globally at kilometric scale in a land surface model. A sensitivity experiment assuming a 100% urban world showed spatially averaged northern hemisphere 2‐m temperatures increased by 0.54 K (January) and 0.42 K (July) at night caused by changes in the albedo, emissivity, roughness, and thermal and hydrological properties. Global ∼1‐km resolution simulations using ancillary urban mapping information produce an urban heat island effect over major and minor conurbations. Only major conurbations were well represented at ∼9‐km resolution. Results from SCM simulations show a heightening of the planetary boundary layer over city sites, with the largest enhancements occurring at night in July (84 ± 48 m) caused by an increased sensible heat flux. These initial developments show the importance of a high‐resolution urban representation within NWP models. Improved parameterization and mapping will enable an online representation of energy, water, and trace gas fluxes over residential areas.

In this study, we present a novel monitoring methodology that combines satellite retrievals and forecasts to detect local CH4 concentration anomalies worldwide. These anomalies are caused by rapidly changing anthropogenic emissions that significantly contribute to the CH4 atmospheric budget and by biases in the satellite retrieval data. The method uses high-resolution (7 km × 7 km) retrievals of total column CH4 from the TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel 5 Precursor satellite. Observations are combined with high-resolution CH4 forecasts (∼ 9 km) produced by the Copernicus Atmosphere Monitoring Service (CAMS) to provide departures (observations minus forecasts) at close to the satellite's native resolution at appropriate time. Investigating these departures is an effective way to link satellite measurements and emission inventory data in a quantitative manner. We perform filtering on the departures to remove the synoptic-scale and meso-alpha-scale biases in both forecasts and satellite observations. We then apply a simple classification scheme to the filtered departures to detect anomalies and plumes that are missing (e.g. pipeline or facility leaks), underreported or overreported (e.g. depleted drilling fields) in the CAMS emissions. The classification method also shows some limitations to detect emission anomalies only due to local satellite retrieval biases linked to albedo and scattering issues.

The growth in anthropogenic carbon dioxide (CO2) emissions acts as a major climate change driver, which has widespread implications across society, influencing the scientific, political, and public sectors. For an increased understanding of the CO2 emission sources, patterns, and trends, a link between the emission inventories and observed CO2 concentrations is best established via Earth system modelling and data assimilation. Bringing together the different pieces of the puzzle of a very different nature (measurements, reported statistics, and models), it is of utmost importance to know their level of confidence and boundaries well. Inversions disaggregate the variation in observed atmospheric CO2 concentration to variability in CO2 emissions by constraining the regional distribution of CO2 fluxes, derived either bottom-up from statistics or top-down from observations. The level of confidence and boundaries for each of these CO2 fluxes is as important as their intensity, though often not available for bottom-up anthropogenic CO2 emissions. This study provides a postprocessing tool CHE_UNC_APP for anthropogenic CO2 emissions to help assess and manage the uncertainty in the different emitting sectors. The postprocessor is available under https://doi.org/10.5281/zenodo.5196190 (Choulga et al., 2021). Recommendations are given for regrouping the sectoral emissions, taking into account their uncertainty instead of their statistical origin; for addressing local hot spots; for the treatment of sectors with small budget but uncertainties larger than 100 %; and for the assumptions around the classification of countries based on the quality of their statistical infrastructure. This tool has been applied to the EDGARv4.3.2_FT2015 dataset, resulting in seven input grid maps with upper- and lower-half ranges of uncertainty for the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System. The dataset is documented and available under https://doi.org/10.5281/zenodo.3967439 (Choulga et al., 2020). While the uncertainty in most emission groups remains relatively small (5 %–20 %), the largest contribution (usually over 40 %) to the total uncertainty is determined by the OTHER group (of fuel exploitation and transformation but also agricultural soils and solvents) at the global scale. The uncertainties have been compared for selected countries to those reported in the inventories submitted to the United Nations Framework Convention on Climate Change and to those assessed for the European emission grid maps of the Netherlands Organisation for Applied Scientific Research. Several sensitivity experiments are performed to check (1) the country dependence (by analysing the impact of assuming either a well- or less well-developed statistical infrastructure), (2) the fuel type dependence (by adding explicit information for each fuel type used per activity from the Intergovernmental Panel on Climate Change), and (3) the spatial source distribution dependence (by aggregating all emission sources and comparing the effect against an even redistribution over the country). The first experiment shows that the SETTLEMENTS group (of energy for buildings) uncertainty changes the most when development level is changed. The second experiment shows that fuel-specific information reduces uncertainty in emissions only when a country uses several different fuels in the same amount; when a country mainly uses the most globally typical fuel for an activity, uncertainty values computed with and without detailed fuel information are the same. The third experiment highlights the importance of spatial mapping.

Urban areas make up only a small fraction of the Earth’s surface; however, they are home to over 50% of the global population. Accurate numerical weather prediction (NWP) forecasts in these areas offer clear societal benefits; however, land-atmosphere interactions are significantly different between urban and non-urban environments. Forecasting urban weather requires higher model resolution than the size of the urban domain, which is often achievable by regional but not global NWP models. Here we present the preliminary implementation of an urban scheme within the land surface component of the global Integrated Forecasting System (IFS), at recently developed ~1km horizontal resolution. We evaluate the representation error of fluxes and NWP variables at coarser resolutions (~9 km and ~31 km), using the high resolution as truth. We evaluate the feasibility of the scheme and its urban representation at ~1km scales. Availability of urban mapping data limit the affordable complexity of the global scheme; however, using generalisations model performance is improved over urban sites, even adopting simple schemes, and the modelled Urban Heat Island effects show broad agreement with observations. Several directions for future work are explored including a more complex urban representation, restructuring of the urban tiling and the introduction of an urban emissions model for trace gas emissions.

The CO2 Human Emissions (CHE) project has been tasked by the European Commission to prepare the development of a European capacity to monitor anthropogenic CO2 emissions. The monitoring of fossil fuel CO2 emissions has to come with a sufficiently low uncertainty in order to be useful for policymakers. In this context, the main approaches to estimate fossil fuel emissions, apart from bottom-up inventories, are based on inverse transportmodeling either on its own or within a coupled carbon cycle fossil fuel data assimilation system. Both approaches make use of atmospheric CO2 and other tracers (e.g., CO and NOx) and rely on the availability of prior fossil fuel CO2 emission estimates and uncertainties (as well as biogenic fluxes for the transport inverse modeling). For a robust estimate of the uncertainty, information from different sources needs to be brought together.A methodology to calculate yearly and monthly anthropogenic CO2 emission uncertainties based on IPCC guidelines (2006 IPCC Guidelines for National Greenhouse Gas Inventories + its 2019 Refinements) has been developed. Emission uncertainties are calculated for all world countries, under the assumption of two categories of world countries, depending on whether the country’s statistical infrastructure is well or less developed. For well-developed statistical infrastructure, emission uncertainties are lower, while less developed statistical infrastructure countries have higher emission uncertainties. A sensitivity analysis is investigating the impact of the well or less developed infrastructure assumption for several countries on the global emission uncertainty. Sensitivity experiments with different anthropogenic CO2 sources distributions, as well as the first results on using these prior anthropogenic CO2 uncertainties in ensemble perturbation runs will be presented.

Effective monitoring of global wildfire activity requires comprehensive knowledge of changing environmental (including atmospheric and hydrological) conditions, fuel availability and routine observations of fire locations and intensity. The European Centre for Medium-Range Weather Forecasts (ECMWF) through its operation of, and contribution to, different Copernicus Services is in a unique position to provide detailed information on the conditions leading to wildland fire activity, the evolution of wildfires, and their potential impacts, when they occur. Fire weather forecasts from the Copernicus Emergency Management Service, and surface climate anomalies from the Copernicus Climate Change Service both provide context to the environmental conditions required for wildfires to persist. Analyses based on observations of fire radiative power, along with analyses and forecasts of associated atmospheric pollutants, from the Copernicus Atmosphere Monitoring Service aid in quantifying the scale and intensity in near-real-time and the subsequent atmospheric impacts. During 2019, regions of anomalously hot and dry surface conditions in Arctic Siberia and southeast Australia experienced large-scale, long-duration wildfires which burned thousands of square kilometres with a total intensity that was significantly above the average of the previous 16 years of data in those regions. We present an overview of the evolution of fire activity in Siberia between June-August 2019, and Australia between September 2019-January 2020, based on ECMWF/Copernicus data for fire weather, climate anomalies and active fires. We will show that the different datasets, while being relatively independent, show a strong correspondence and provide a wealth of information vital to understanding global wildfires, their underlying causes and environmental impacts.

Abstract. In this study we present a novel monitoring methodology to detect local CH4 concentration anomalies worldwide that are related to rapidly changing anthropogenic emissions that significantly contribute to the CH4 atmospheric budget. The method uses high resolution (7 km × 7 km) retrievals of total column CH4 from the Tropospheric Monitoring Instrument (TROPOMI) onboard the Sentinel 5 Precursor satellite. Observations are combined with high resolution CH4 forecasts (~ 9 km) produced by the Copernicus Atmosphere Monitoring Service (CAMS) to provide departures (observations minus forecasts) close to the native satellite resolution at appropriate time. Investigating the departures is an effective way to link satellite measurements and emission inventory data in a quantitative manner. We perform filtering on the departures to remove the large-scale biases on both forecasts and satellite observations. We then use a simple classification on the filtered departures to detect anomalies and plumes coming from CAMS emissions that are missing (e.g. pipeline or facility leaks), under-reported or over-reported (e.g. depleted drilling fields). Additionally, the classification helps to detect local satellite retrieval errors due to land surface albedo issues.

Incremental 4D‐Var is a data assimilation algorithm used routinely at operational numerical weather prediction (NWP) centres worldwide. The algorithm solves a series of quadratic minimization problems (inner‐loops) obtained from linear approximations of the forward model around nonlinear trajectories (outer‐loops). Since most of the computational burden is associated with the inner‐loops, many studies have focused on developing computationally efficient algorithms to solve the least‐square quadratic minimization problem, in particular through time parallelization. This paper presents the first implementation and testing of a recently proposed method for parallelizing incremental 4D‐Var, the Randomized Incremental Optimal Technique (RIOT), which replaces the traditional sequential conjugate gradient (CG) iterations in the inner‐loop of the minimization with fully parallel randomized singular value decomposition (RSVD) of the preconditioned Hessian of the cost function. RIOT is tested using the standard Lorenz‐96 model (L‐96) as well as two realistic high‐dimensional atmospheric source inversion problems based on aircraft observations of black carbon concentrations. A new outer‐loop preconditioning technique tailored to RSVD was introduced to improve convergence stability and performance. Results obtained with the L‐96 system show that the performance improvement from RIOT compared to standard CG algorithms increases significantly with nonlinearities. Overall, in the realistic black carbon source inversion experiments, RIOT reduces the wall‐clock time of the 4D‐Var minimization by a factor of 2 to 3, at the cost of a factor of 4 to 10 increase in energy cost due to the large number of parallel cores used. Furthermore, RIOT enables reduction of the wall‐clock time computation of the analysis‐error covariance matrix by a factor of 40 compared to a standard iterative Lanczos approach. Finally, as evidenced in this study, implementation of RIOT in an operational NWP system will require a better understanding of its convergence properties as a function of the Hessian characteristics and, in particular, the degree of freedom for signal of the inverse problem.

The European Union’s Copernicus Atmosphere Monitoring Service (CAMS) operationally provides daily forecasts of global atmospheric composition. It uses the ECMWF Integrated Forecasting System (IFS), which includes meteorological and atmospheric composition variables, such as reactive gases, greenhouse gases and aerosols, for its global forecasts and reanalyses. The current green-house gases operational suite monitors CH4 and CO2 and assimilates TANSO and IASI retrievals for both species. The TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor (S5P) satellite launched in October 2017 yields a wealth of atmospheric composition data, including CH4 retrievals at unprecedented high horizontal resolution (7km) and up to daily revisit time. We used the IFS to perform monitoring experiments at different horizontal resolutions (25 km and 9 km). We also performed first data assimilation experiments at 25 km horizontal resolution. This first set of monitoring experiments shows the potential of the TROPOMI CH4 retrievals to correct known biases that exist in the current CAMS analyses and forecasts. Assimilation experiments of TROPOMI CH4 shows that adding the instrument in the operational chain would significantly improve the analysis and forecasts. Detection of CH4 sources seen by TROPOMI compared to CAMS also shows the potential of the instrument to inform on and infer anthropogenic and natural sources. For example, discrepancies between TROPOMI retrievals and CAMS fields in the CH4 levels associated with oil and gas extraction activities show very promising perspectives for monitoring and analysis of CH4 concentration and emissions. We will finally discuss the challenges and progress made towards performing inversions using the IFS operational system.

Atmospheric flux inversions use observations of atmospheric CO2 to provide anthropogenic and biogenic CO2 flux estimates at a range of spatio-temporal scales. Inversions require prior flux, a forward model and observation errors to estimate posterior fluxes and uncertainties. Here, we investigate the forward transport error and the associated biogenic feedback in an Earth system model (ESM) context. These errors can occur from uncertainty in the initial meteorology, the analysis fields used, or the advection schemes and physical parameterisation of the model. We also explore the spatio-temporal variability and flow-dependent error covariances. We then compare the error with the atmospheric response to uncertainty in the prior anthropogenic emissions. Although transport errors are variable, average total-column CO2 (XCO2) transport errors over anthropogenic emission hotspots (0.1–0.8 ppm) are comparable to, and often exceed, prior monthly anthropogenic flux uncertainties projected onto the same space (0.1–1.4 ppm). Average near-surface transport errors at three sites (Paris, Caltech and Tsukuba) range from 1.7 to 7.2 ppm. The global average XCO2 transport error standard deviation plateaus at ∼0.1 ppm after 2–3 d, after which atmospheric mixing significantly dampens the concentration gradients. Error correlations are found to be highly flow dependent, with XCO2 spatio-temporal correlation length scales ranging from 0 to 700 km and 0 to 260 min. Globally, the average model error caused by the biogenic response to atmospheric meteorological uncertainties is small (<0.01 ppm); however, this increases over high flux regions and is seasonally dependent (e.g. the Amazon; January and July: 0.24±0.18 ppm and 0.13±0.07 ppm). In general, flux hotspots are well-correlated with model transport errors. Our model error estimates, combined with the atmospheric response to anthropogenic flux uncertainty, are validated against three Total Carbon Observing Network (TCCON) XCO2 sites. Results indicate that our model and flux uncertainty account for 21 %–65 % of the total uncertainty. The remaining uncertainty originates from additional sources, such as observation, numerical and representation errors, as well as structural errors in the biogenic model. An underrepresentation of transport and flux uncertainties could also contribute to the remaining uncertainty. Our quantification of CO2 transport error can be used to help derive accurate posterior fluxes and error reductions in future inversion systems. The model uncertainty diagnosed here can be used with varying degrees of complexity and with different modelling techniques by the inversion community.

The CO2 Human Emission (CHE) project is an European initiative bringing together a consortium of 22 European partners to build a prototype global CO2 source inversion system that can provide policy-relevant information on the spatiotemporal characteristics of anthropogenic CO2 emissions. This prototype shall evolve toward a new Copernicus CO2 service, which will provide a Monitoring and Verification Support (MVS) capacity that can address the challenge of the global stocktake (GST) devised under the Paris Agreement. The global inversion system will build on existing operational infrastructures (CAMS, C3S) at the European Centre for Medium-range Weather Forecast (ECMWF) to exploit ground-based measurements as well as space-based observations from current and future satellite missions (e.g., Sentinel 5p and CO2M). We will present ongoing efforts at ECMWF to develop a source inversion capability in the current operational Integrated Forecasting System (IFS), which will serve as the basis for the future global CO2 inversion prototype. Preliminary results will be discussed, that include model transport error estimations based on Monte-Carlo ensemble simulations as well as the first chemical source optimization experiments performed with the IFS 4D-Var system.

As part of the Copernicus Atmosphere Monitoring Service (CAMS), operated by ECMWF on behalf of the European Commission, global analyses and forecasts of atmospheric composition have been produced operationally since 2015. These were built on many years of previous work under the GEMS and MACC projects, which began producing regular forecasts in 2007. Since the transition to an operational service, there have continued to be many new developments and improvements to the system in five major upgrades, including increased horizontal and vertical resolution, updated emissions and paramterisations, additional species such as nitrate aerosol, as well as updates to the underlying meteorological model and data assimilation. The components of this system (aerosols, gas-phase chemistry, meteorology and the ocean) are also now coupled more tightly via active feedbacks then ever before. In this interactive presentation, we will demonstrate the impact of a number of these developments on the performance of the resulting global air quality forecasts, alongside the continuing evolution of our approaches to assessing model improvement against independent in-situ and remote-sensing observations from a variety of platforms. Because the continuing evolution of an operational system can make the analysis of long-term trends problematic, we will also contrast this with the CAMS global reanalysis product, which (while not using the very latest version of the model) do provide a consistent long-term dataset from 2003 onwards.