Search

Methane (CH4) is an important anthropogenic greenhouse gas, and its rising concentration in the atmosphere contributes significantly to global warming. A comparatively small number of highly emitting persistent methane sources are responsible for a large share of global methane emissions. The identification and quantification of these sources, which often show large uncertainties regarding their emissions or locations, are important to support mitigating climate change. Daily global column-averaged dry air mole fractions of atmospheric methane (XCH4) are retrieved from radiance measurements of the TROPOspheric Monitoring Instrument (TROPOMI) on board on the Sentinel-5 Precursor (S5P) satellite with a moderately high spatial resolution, enabling the detection and quantification of localized methane sources. We developed a fully automated algorithm to detect regions with persistent methane enhancement and to quantify their emissions using a monthly TROPOMI XCH4 dataset from the years 2018–2021. We detect 217 potential persistent source regions (PPSRs), which account for approximately 20 % of the total bottom-up emissions. By comparing the PPSRs in a spatial analysis with anthropogenic and natural emission databases, we conclude that 7.8 % of the detected source regions are dominated by coal, 7.8 % by oil and gas, 30.4 % by other anthropogenic sources like landfills or agriculture, 7.3 % by wetlands, and 46.5 % by unknown sources. Many of the identified PPSRs are in well-known source regions, like the Permian Basin in the USA, which is a large production area for oil and gas; the Bowen Basin coal mining area in Australia; or the Pantanal Wetlands in Brazil. We perform a detailed analysis of the PPSRs with the 10 highest emission estimates, including the Sudd Wetland in South Sudan, an oil- and gas-dominated area on the west coast in Turkmenistan, and one of the largest coal production areas in the world, the Kuznetsk Basin in Russia. The calculated emission estimates of these source regions are in agreement within the uncertainties in results from other studies but are in most of the cases higher than the emissions reported by emission databases. We demonstrate that our algorithm is able to automatically detect and quantify persistent localized methane sources of different source type and shape, including larger-scale enhancements such as wetlands or extensive oil- and gas-production basins.

Global crude steel production is expected to continue to increase in the coming decades to meet the demands of the growing world population. Currently, the dominant steelmaking technology worldwide is the conventional highly CO2-intensive blast furnace–basic oxygen furnace production route (also known as the Linz–Donawitz process), which uses iron ore as raw material and coke as a reducing agent. As a result, large quantities of special gases that are rich in carbon monoxide (CO) are by-products of the various stages of the steelmaking process. Given the challenges associated with satellite-based estimates of carbon dioxide (CO2) emissions at the scale of emitting installations due to significant background levels, co-emitted CO may serve as a valuable indicator of the carbon footprint of steel plants. We show that regional CO release from steel production sites can be monitored from space using 5 years of measurements (2018–2022) from the TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor satellite, benefiting from its relatively high spatial resolution and daily global coverage. We analyse all German steel plants with blast furnaces and basic oxygen furnaces and obtain associated CO emissions in the range of 50–400 kt yr−1 per site. A comparison with the respective CO2 emissions on the level of emitting installations available from emissions trading data of the European Union Emissions Trading System yields a linear relationship with a sector-specific CO/CO2 emission ratio for the analysed steelworks of 3.24 % [2.73–3.89; 1σ], suggesting the feasibility of using CO as a proxy for CO2 emissions from comparable steel production sites. An evaluation at other steel production sites indicates that the derived CO/CO2 emission ratio is also representative of other highly optimised state-of-the-art Linz–Donawitz steelworks outside Germany and that the emission ratio is potentially valuable for estimating sector-specific CO2 emissions from remotely sensed CO emissions, provided that the underlying CO emission estimate is not affected by other sources.

The Tropospheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor (S5P) satellite is part of the latest generation of trace gas monitoring satellites and provides a new level of spatio-temporal information with daily global coverage, which enables the calculation of daily globally averaged CH4 concentrations. To investigate changes in atmospheric methane, the background CH4 level (i.e. the CH4 concentration without seasonal and short-term variations) has to be determined. CH4 growth rates vary in a complex manner and high-latitude zonal averages may have gaps in the time series, and thus simple fitting methods do not produce reliable results. In this paper we present an approach based on fitting an ensemble of dynamic linear models (DLMs) to TROPOMI data, from which the best model is chosen with the help of cross-validation to prevent overfitting. This method is computationally fast and is not dependent on additional inputs, allowing for fast and continuous analysis of the most recent time series data. We present results of global annual methane increases (AMIs) for the first 4.5 years of S5P/TROPOMI data, which show good agreement with AMIs from other sources. Additionally, we investigated what information can be derived from zonal bands. Due to the fast meridional mixing within hemispheres, we use zonal growth rates instead of AMIs, since they provide a higher temporal resolution. Clear differences can be observed between Northern Hemisphere and Southern Hemisphere growth rates, especially during 2019 and 2022. The growth rates show similar patterns within the hemispheres and show no short-term variations during the years, indicating that air masses within a hemisphere are well-mixed during a year. Additionally, the growth rates derived from S5P/TROPOMI data are largely consistent with growth rates derived from Copernicus Atmospheric Monitoring Service (CAMS) global-inversion-optimized (CAMS/INV) data, which use surface observations. In 2019 a reduction in growth rates can be observed for the Southern Hemisphere, while growth rates in the Northern Hemisphere stay stable or increase. During 2020 a strong increase in Southern Hemisphere growth rates can be observed, which is in accordance with recently reported increases in Southern Hemisphere wetland emissions. In 2022 the reduction in the global AMI can be attributed to decreased growth rates in the Northern Hemisphere, while growth rates in the Southern Hemisphere remain high. Investigations of fluxes from CAMS/INV data support these observations and suggest that the Northern Hemisphere decrease is mainly due to the decrease in anthropogenic fluxes, while in the Southern Hemisphere, wetland fluxes continued to rise. While the continued increase in Southern Hemisphere wetland fluxes agrees with existing studies about the causes of observed methane trends, the difference between Northern Hemisphere and Southern Hemisphere methane increases in 2022 has not been discussed before and calls for further research.

The TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor satellite enables the accurate determination of atmospheric methane (CH4) and carbon monoxide (CO) abundances at high spatial resolution and global daily sampling. Due to its wide swath and sampling, the global distribution of both gases can be determined in unprecedented detail. The scientific retrieval algorithm Weighting Function Modified Differential Optical Absorption Spectroscopy (WFMD) has proven valuable in simultaneously retrieving the atmospheric column-averaged dry-air mole fractions XCH4 and XCO from TROPOMI's radiance measurements in the shortwave infrared (SWIR) spectral range. Here we present recent improvements of the algorithm which have been incorporated into the current version (v1.8) of the TROPOMI/WFMD product. This includes processing adjustments such as increasing the polynomial degree to 3 in the fitting procedure to better account for possible spectral albedo variations within the fitting window and updating the digital elevation model to minimise topography-related biases. In the post-processing, the machine-learning-based quality filter has been refined using additional data when training the random forest classifier to further reduce scenes with residual cloudiness that are incorrectly classified as good. In particular, the cloud filtering over the Arctic ocean is considerably improved. Furthermore, the machine learning calibration, addressing systematic errors due to simplifications in the forward model or instrumental issues, has been optimised. By including an additional feature associated with the fitted polynomial when training the corresponding random forest regressor, spectral albedo variations are better accounted for. To remove vertical stripes in the XCH4 and XCO data, an efficient orbit-wise destriping filter based on combined wavelet–Fourier filtering has been implemented, while optimally preserving the original spatial trace gas features. The temporal coverage of the data records has been extended to the end of April 2022, covering a total length of 4.5 years since the start of the mission, and will be further extended in the future. Validation with the ground-based Total Carbon Column Observing Network (TCCON) demonstrates that the implemented improvements reduce the pseudo-noise component of the products, resulting in an improved random error. The XCH4 and XCO products have similar spatial coverage from year to year including high latitudes and the oceans. The analysis of annual growth rates reveals accelerated growth of atmospheric methane during the covered period, in line with observations at marine surface sites of the Global Monitoring Division of NOAA's Earth System Research Laboratory, which reported consecutive annual record increases over the past 2 years of 2020 and 2021.

The Sentinel-5 Precursor (S5P) mission was launched on October 2017 and has since provided data with high spatio-temporal resolution using its remote sensing instrument, the TROPOspheric Monitoring Instrument (TROPOMI). The latter is a nadir viewing passive grating imaging spectrometer. The mathematical inversion of the TROPOMI data yields retrievals of different trace gas and aerosol data products. The column-averaged dry air mole fraction of methane (XCH4) is the product of interest to this study. The daily global coverage of the atmospheric methane mole fraction data enables the analysis of the methane distribution and variation on large scales and also to estimate surface emissions. The spatio-temporal high-resolution satellite data are potentially particularly valuable in remote regions, such as the Arctic, where few ground stations and in-situ measurements are available. In addition to the operational Copernicus S5P total-column averaged dry air mole fraction methane data product developed by SRON, the scientific TROPOMI/WFMD algorithm data product v1.5 (WFMD product) was generated at the Institute of Environmental Physics at the University of Bremen. In this study we focus on the assessment of both S5P XCH4 data products over Greenland and find that spatial maps of both products show distinct features along the coast lines. Anomalies up to and exceeding 100 ppb are observed and stand out in comparison to the otherwise smooth changes in the methane distribution. These features are more pronounced for the operational product compared to the WFMD product. The spatial patterns correlate with the difference between the GMTED2010 digital elevation model (DEM) used in the retrievals and a more recent topography data set indicating that inaccuracies in the assumed surface elevation are the origin of the observed features. These correlations are stronger for the WFMD product. In order to evaluate the impact of the topography dataset on the retrieval we reprocess the WFMD product with updated elevation data. We find a significant reduction of the localized features when GMTED2010 is replaced by recent topography data over Greenland based on ICESat-2 data. This study shows the importance of the chosen topography data on retrieved dry air mole fractions. The use of an accurate and as up-to-date as possible DEM is advised for all S5P data products as well as for future missions which rely on DEM as input data. A modification based on this study is planned to be introduced in the next version of the WFMD data product.

We show new results from an updated version of the Fast atmOspheric traCe gAs retrievaL (FOCAL) retrieval method applied to measurements of the Greenhouse gases Observing SATellite (GOSAT) and its successor GOSAT-2. FOCAL was originally developed for estimating the total column carbon dioxide mixing ratio (XCO2) from spectral measurements made by the Orbiting Carbon Observatory-2 (OCO-2). However, depending on the available spectral windows, FOCAL also successfully retrieves total column amounts for other atmospheric species and their uncertainties within one single retrieval. The main focus of the current paper is on methane (XCH4; full-physics and proxy product), water vapour (XH2O) and the relative ratio of semi-heavy water (HDO) to water vapour (δD). Due to the extended spectral range of GOSAT-2, it is also possible to derive information on carbon monoxide (XCO) and nitrous oxide (XN2O) for which we also show first results. We also present an update on XCO2 from both instruments. For XCO2, the new FOCAL retrieval (v3.0) significantly increases the number of valid data compared with the previous FOCAL retrieval version (v1) by 50 % for GOSAT and about a factor of 2 for GOSAT-2 due to relaxed pre-screening and improved post-processing. All v3.0 FOCAL data products show reasonable spatial distribution and temporal variations. Comparisons with the Total Carbon Column Observing Network (TCCON) result in station-to-station biases which are generally in line with the reported TCCON uncertainties. With this updated version of the GOSAT-2 FOCAL data, we provide a first total column average XN2O product. Global XN2O maps show a gradient from the tropics to higher latitudes on the order of 15 ppb, which can be explained by variations in tropopause height. The new GOSAT-2 XN2O product compares well with TCCON. Its station-to-station variability is lower than 2 ppb, which is about the magnitude of the typical N2O variations close to the surface. However, both GOSAT-2 and TCCON measurements show that the seasonal variations in the total column average XN2O are on the order of 8 ppb peak-to-peak, which can be easily resolved by the GOSAT-2 FOCAL data. Noting that only few XN2O measurements from satellites exist so far, the GOSAT-2 FOCAL product will be a valuable contribution in this context.

Several ambient air quality records corroborate the severe and persistent degradation of air quality over northern India during the winter months, with evidence of a continued, increasing trend of pollution across the Indo-Gangetic Plain (IGP) over the past decade. A combination of atmospheric dynamics and uncertain emissions, including the post-monsoon agricultural stubble burning, make it challenging to resolve the role of each individual factor. Here we demonstrate the potential use of an atmospheric transport model, the Weather Research and Forecasting model coupled with chemistry (WRF–Chem) to identify and quantify the role of transport mechanisms and emissions on the occurrence of the pollution events. The investigation is based on the use of carbon monoxide (CO) observations from the TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor satellite and the surface measurement network, as well as the WRF–Chem simulations, to investigate the factors contributing to CO enhancement over India during November 2018. We show that the simulated column-averaged dry air mole fraction (XCO) is largely consistent with TROPOMI observations, with a spatial correlation coefficient of 0.87. The surface-level CO concentrations show larger sensitivities to boundary layer dynamics, wind speed, and diverging source regions, leading to a complex concentration pattern and reducing the observation-model agreement with a correlation coefficient ranging from 0.41 to 0.60 for measurement locations across the IGP. We find that daily satellite observations can provide a first-order inference of the CO transport pathways during the enhanced burning period, and this transport pattern is reproduced well in the model. By using the observations and employing the model at a comparable resolution, we confirm the significant role of atmospheric dynamics and residential, industrial, and commercial emissions in the production of the exorbitant level of air pollutants in northern India. We find that biomass burning plays only a minimal role in both column and surface enhancements of CO, except for the state of Punjab during the high pollution episodes. While the model reproduces observations reasonably well, a better understanding of the factors controlling the model uncertainties is essential for relating the observed concentrations to the underlying emissions. Overall, our study emphasizes the importance of undertaking rigorous policy measures, mainly focusing on reducing residential, commercial, and industrial emissions in addition to actions already underway in the agricultural sectors.

The COVID-19 pandemic resulted in reduced anthropogenic carbon dioxide (CO2) emissions during 2020 in large parts of the world. To investigate whether a regional-scale reduction of anthropogenic CO2 emissions during the COVID-19 pandemic can be detected using space-based observations of atmospheric CO2, we have analysed a small ensemble of OCO-2 and GOSAT satellite retrievals of column-averaged dry-air mole fractions of CO2, i.e. XCO2. We focus on East China and use a simple data-driven analysis method. We present estimates of the relative change of East China monthly emissions in 2020 relative to previous periods, limiting the analysis to October-to-May periods to minimize the impact of biogenic CO2 fluxes. The ensemble mean indicates an emission reduction by approximately 10 % ± 10 % in March and April 2020. However, our results show considerable month-to-month variability and significant differences across the ensemble of satellite data products analysed. For example, OCO-2 suggests a much smaller reduction (∼ 1 %–2 % ± 2 %). This indicates that it is challenging to reliably detect and to accurately quantify the emission reduction with current satellite data sets. There are several reasons for this, including the sparseness of the satellite data but also the weak signal; the expected regional XCO2 reduction is only on the order of 0.1–0.2 ppm. Inferring COVID-19-related information on regional-scale CO2 emissions using current satellite XCO2 retrievals likely requires, if at all possible, a more sophisticated analysis method including detailed transport modelling and considering a priori information on anthropogenic and natural CO2 surface fluxes.

The switch from the use of coal to natural gas or oil for energy generation potentially reduces greenhouse gas emissions and thus the impact on global warming and climate change because of the higher energy creation per CO2 molecule emitted. However, the climate benefit over coal is offset by methane (CH4) leakage from natural gas and petroleum systems, which reverses the climate impact mitigation if the rate of fugitive emissions exceeds the compensation point at which the global warming resulting from the leakage and the benefit from the reduction of coal combustion coincide. Consequently, an accurate quantification of CH4 emissions from the oil and gas industry is essential to evaluate the suitability of natural gas and petroleum as bridging fuels on the way to a carbon-neutral future. We show that regional CH4 release from large oil and gas fields can be monitored from space by using dense daily recurrent measurements of the TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite to quantify emissions and leakage rates. The average emissions for the time period 2018/2019 from the five most productive basins in the United States, the Permian, Appalachian, Eagle Ford, Bakken, and Anadarko, are estimated to be 3.18±1.13, 2.36±0.88, 1.37±0.63, 0.89±0.56, and 2.74±0.74 Mt yr−1, respectively. This corresponds to CH4 leakage rates relative to the associated production between 1.2 % and 1.4 % for the first four production regions, which are consistent with bottom-up estimates and likely fall below the break-even leakage rate for immediate climate benefit. For the Anadarko Basin, the fugitive emission rate is larger and amounts to 3.9±1.1 %, which likely exceeds the break-even rate for immediate benefit and roughly corresponds to the break-even rate for a 20-year time horizon. The determined values are smaller than previously derived satellite-based leakage rates for the time period 2009–2011, which was an early phase of hydraulic fracturing, indicating that it is possible to improve the climate footprint of the oil and gas industry by adopting new technologies and that efforts to reduce methane emissions have been successful. For two of the world's largest natural gas fields, Galkynysh and Dauletabad in Turkmenistan, we find collective methane emissions of 3.26±1.17 Mt yr−1, which corresponds to a leakage rate of 4.1±1.5 %, suggesting that the Turkmen energy industry is not employing methane emission avoidance strategies and technologies as successfully as those currently widely used in the United States. The leakage rates in Turkmenistan and in the Anadarko Basin indicate that there is potential to reduce fugitive methane emissions from natural gas and petroleum systems worldwide. In particular, relatively newly developed oil and gas plays appear to have larger leakage rates compared to more mature production areas.

Due to proceeding climate change, some regions such as California face rising weather extremes with dry periods becoming warmer and drier, entailing the risk that wildfires and associated air pollution episodes will continue to increase. November 2018 turned into one of the most severe wildfire episodes on record in California, with two particularly destructive wildfires spreading concurrently through the north and the south of the state. Both fires ignited at the wildland–urban interface, causing many civilian fatalities and forcing the total evacuation of several cities and communities. Here we demonstrate that the inherent carbon monoxide (CO) emissions of the wildfires and subsequent transport can be observed from space by analysing radiance measurements of the TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite in the shortwave infrared spectral range. From the determined CO distribution we assess the corresponding air quality burden in major Californian cities caused by the fires and discuss the associated uncertainties. As a result of the prevailing wind conditions, the largest CO load during the first days of the fires is found in Sacramento and San Francisco, with city area averages reaching boundary layer concentration anomalies of about 2.5 mg CO m−3. Even the most polluted city scenes likely comply with the national ambient air quality standards (10 mg CO m−3 with 8 h averaging time). This finding based on dense daily recurrent satellite monitoring is consistent with isolated ground-based air quality measurements.

Glyoxal (CHOCHO) and formaldehyde (HCHO) are intermediate products in the tropospheric oxidation of the majority of volatile organic compounds (VOCs). CHOCHO is also a precursor of secondary organic aerosol (SOA) in the atmosphere. CHOCHO and HCHO are released from biogenic, anthropogenic, and pyrogenic sources. CHOCHO and HCHO tropospheric lifetimes are typically considered to be short during the daytime at mid-latitudes (e.g. several hours), as they are rapidly removed from the atmosphere by their photolysis, oxidation by OH, and uptake on particles or deposition. At night and at high latitudes, tropospheric lifetimes increase to many hours or even days. Previous studies demonstrated that CHOCHO and HCHO vertical column densities (VCDs) are well retrieved from space-borne observations using differential optical absorption spectroscopy (DOAS). In this study, we present CHOCHO and HCHO VCDs retrieved from measurements by TROPOMI (TROPOspheric Monitoring Instrument), launched on the Sentinel-5 Precursor (S5P) platform in October 2017. We observe strongly elevated amounts of CHOCHO and HCHO during the 2018 fire season in British Columbia, Canada, where a large number of fires occurred in August. CHOCHO and HCHO plumes from individual fire hot spots are observed in air masses travelling over distances of up to 1500 km, i.e. much longer than expected for the relatively short tropospheric lifetime expected for CHOCHO and HCHO. Comparison with simulations by the particle dispersion model FLEXPART (FLEXible PARTicle dispersion model) indicates that effective lifetimes of 20 h and more are needed to explain the observations of CHOCHO and HCHO if they decay in an effective first-order process. FLEXPART used in the study calculates accurately the transport. In addition an exponential decay, in our case assumed to be photochemical, of a species along the trajectory is added. We have used this simple approach to test our assumption that CHOCHO and HCHO are created in the fires and then decay at a constant rate in the plume as it is transported. This is clearly not the case and we infer that CHOCHO and HCHO are either efficiently recycled during transport or continuously formed from the oxidation of longer-lived precursors present in the plume, or possibly a mixture of both. We consider the best explanation of the observed CHOCHO and HCHO VCD in the plumes of the fire is that they are produced by oxidation of longer-lived precursors, which were also released by the fire and present in the plume.

Satellite retrievals of column-averaged dry-air mole fractions of carbon dioxide (CO2) and methane (CH4), denoted XCO2 and XCH4, respectively, have been used in recent years to obtain information on natural and anthropogenic sources and sinks and for other applications such as comparisons with climate models. Here we present new data sets based on merging several individual satellite data products in order to generate consistent long-term climate data records (CDRs) of these two Essential Climate Variables (ECVs). These ECV CDRs, which cover the time period 2003–2018, have been generated using an ensemble of data products from the satellite sensors SCIAMACHY/ENVISAT and TANSO-FTS/GOSAT and (for XCO2) for the first time also including data from the Orbiting Carbon Observatory 2 (OCO-2) satellite. Two types of products have been generated: (i) Level 2 (L2) products generated with the latest version of the ensemble median algorithm (EMMA) and (ii) Level 3 (L3) products obtained by gridding the corresponding L2 EMMA products to obtain a monthly 5∘×5∘ data product in Obs4MIPs (Observations for Model Intercomparisons Project) format. The L2 products consist of daily NetCDF (Network Common Data Form) files, which contain in addition to the main parameters, i.e., XCO2 or XCH4, corresponding uncertainty estimates for random and potential systematic uncertainties and the averaging kernel for each single (quality-filtered) satellite observation. We describe the algorithms used to generate these data products and present quality assessment results based on comparisons with Total Carbon Column Observing Network (TCCON) ground-based retrievals. We found that the XCO2 Level 2 data set at the TCCON validation sites can be characterized by the following figures of merit (the corresponding values for the Level 3 product are listed in brackets) – single-observation random error (1σ): 1.29 ppm (monthly: 1.18 ppm); global bias: 0.20 ppm (0.18 ppm); and spatiotemporal bias or relative accuracy (1σ): 0.66 ppm (0.70 ppm). The corresponding values for the XCH4 products are single-observation random error (1σ): 17.4 ppb (monthly: 8.7 ppb); global bias: −2.0 ppb (−2.9 ppb); and spatiotemporal bias (1σ): 5.0 ppb (4.9 ppb). It has also been found that the data products exhibit very good long-term stability as no significant long-term bias trend has been identified. The new data sets have also been used to derive annual XCO2 and XCH4 growth rates, which are in reasonable to good agreement with growth rates from the National Oceanic and Atmospheric Administration (NOAA) based on marine surface observations. The presented ECV data sets are available (from early 2020 onwards) via the Climate Data Store (CDS, https://cds.climate.copernicus.eu/, last access: 10 January 2020) of the Copernicus Climate Change Service (C3S, https://climate.copernicus.eu/, last access: 10 January 2020).

Carbon monoxide (CO) is an important atmospheric constituent affecting air quality, and methane (CH4) is the second most important greenhouse gas contributing to human-induced climate change. Detailed and continuous observations of these gases are necessary to better assess their impact on climate and atmospheric pollution. While surface and airborne measurements are able to accurately determine atmospheric abundances on local scales, global coverage can only be achieved using satellite instruments.The TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite, which was successfully launched in October 2017, is a spaceborne nadir-viewing imaging spectrometer measuring solar radiation reflected by the Earth in a push-broom configuration. It has a wide swath on the terrestrial surface and covers wavelength bands between the ultraviolet (UV) and the shortwave infrared (SWIR), combining a high spatial resolution with daily global coverage. These characteristics enable the determination of both gases with an unprecedented level of detail on a global scale, introducing new areas of application.Abundances of the atmospheric column-averaged dry air mole fractions XCO and XCH4 are simultaneously retrieved from TROPOMI's radiance measurements in the 2.3 µm spectral range of the SWIR part of the solar spectrum using the scientific retrieval algorithm Weighting Function Modified Differential Optical Absorption Spectroscopy (WFM-DOAS). This algorithm is intended to be used with the operational algorithms for mutual verification and to provide new geophysical insights. We introduce the algorithm in detail, including expected error characteristics based on synthetic data, a machine-learning-based quality filter, and a shallow learning calibration procedure applied in the post-processing of the XCH4 data. The quality of the results based on real TROPOMI data is assessed by validation with ground-based Fourier transform spectrometer (FTS) measurements providing realistic error estimates of the satellite data: the XCO data set is characterised by a random error of 5.1 ppb (5.8 %) and a systematic error of 1.9 ppb (2.1 %); the XCH4 data set exhibits a random error of 14.0 ppb (0.8 %) and a systematic error of 4.3 ppb (0.2 %). The natural XCO and XCH4 variations are well-captured by the satellite retrievals, which is demonstrated by a high correlation with the validation data (R = 0.97 for XCO and R = 0.91 for XCH4 based on daily averages).We also present selected results from the mission start until the end of 2018, including a first comparison to the operational products and examples of the detection of emission sources in a single satellite overpass, such as CO emissions from the steel industry and CH4 emissions from the energy sector, which potentially allows for the advance of emission monitoring and air quality assessments to an entirely new level.

Despite its key role in climate change, large uncertainties persist in our knowledge of the anthropogenic emissions of carbon dioxide (CO2) and no global observing system exists that allows us to monitor emissions from localized CO2 sources with sufficient accuracy. The Orbiting Carbon Observatory-2 (OCO-2) satellite allows retrievals of the column-average dry-air mole fractions of CO2 (XCO2). However, regional column-average enhancements of individual point sources are usually small, compared to the background concentration and its natural variability, and often not much larger than the satellite's measurement noise. This makes the unambiguous identification and quantification of anthropogenic emission plume signals challenging. NO2 is co-emitted with CO2 when fossil fuels are combusted at high temperatures. It has a short lifetime on the order of hours so that NO2 columns often greatly exceed background and noise levels of modern satellite sensors near sources, which makes it a suitable tracer of recently emitted CO2. Based on six case studies (Moscow, Russia; Lipetsk, Russia; Baghdad, Iraq; Medupi and Matimba power plants, South Africa; Australian wildfires; and Nanjing, China), we demonstrate the usefulness of simultaneous satellite observations of NO2 and XCO2. For this purpose, we analyze co-located regional enhancements of XCO2 observed by OCO-2 and NO2 from the Sentinel-5 Precursor (S5P) satellite and estimate the CO2 plume's cross-sectional fluxes. We take advantage of the nearly simultaneous NO2 measurements with S5P's wide swath and small measurement noise by identifying the source of the observed XCO2 enhancements, excluding interference with remote upwind sources, allowing us to adjust the wind direction, and by constraining the shape of the CO2 plumes. We compare the inferred cross-sectional fluxes with the Emissions Database for Global Atmospheric Research (EDGAR), the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC), and, in the case of the Australian wildfires, with the Global Fire Emissions Database (GFED). The inferred cross-sectional fluxes range from 31 MtCO2 a−1 to 153 MtCO2 a−1 with uncertainties (1σ) between 23 % and 72 %. For the majority of analyzed emission sources, the estimated cross-sectional fluxes agree, within their uncertainty, with either EDGAR or ODIAC or lie somewhere between them. We assess the contribution of multiple sources of uncertainty and find that the dominating contributions are related to the computation of the effective wind speed normal to the plume's cross section. The flux uncertainties are expected to be reduced by the planned European Copernicus anthropogenic CO2 monitoring mission (CO2M), which will provide not only precise measurements with high spatial resolution but also imaging capabilities with a wider swath of simultaneous XCO2 and NO2 observations. Such a mission, particularly if performed by a constellation of satellites, will deliver CO2 emission estimates from localized sources at an unprecedented frequency and level of accuracy.

The growth rate of atmospheric carbon dioxide (CO2) reflects the net effect of emissions and uptake resulting from anthropogenic and natural carbon sources and sinks. Annual mean CO2 growth rates have been determined from satellite retrievals of column-averaged dry-air mole fractions of CO2, i.e. XCO2, for the years 2003 to 2016. The XCO2 growth rates agree with National Oceanic and Atmospheric Administration (NOAA) growth rates from CO2 surface observations within the uncertainty of the satellite-derived growth rates (mean difference ± standard deviation: 0.0±0.3 ppm year−1; R: 0.82). This new and independent data set confirms record-large growth rates of around 3 ppm year−1 in 2015 and 2016, which are attributed to the 2015–2016 El Niño. Based on a comparison of the satellite-derived growth rates with human CO2 emissions from fossil fuel combustion and with El Niño Southern Oscillation (ENSO) indices, we estimate by how much the impact of ENSO dominates the impact of fossil-fuel-burning-related emissions in explaining the variance of the atmospheric CO2 growth rate. Our analysis shows that the ENSO impact on CO2 growth rate variations dominates that of human emissions throughout the period 2003–2016 but in particular during the period 2010–2016 due to strong La Niña and El Niño events. Using the derived growth rates and their uncertainties, we estimate the probability that the impact of ENSO on the variability is larger than the impact of human emissions to be 63 % for the time period 2003–2016. If the time period is restricted to 2010–2016, this probability increases to 94 %.

Methane is an important atmospheric greenhouse gas and an adequate understanding of its emission sources is needed for climate change assessments, predictions, and the development and verification of emission mitigation strategies. Satellite retrievals of near-surface-sensitive column-averaged dry-air mole fractions of atmospheric methane, i.e. XCH4, can be used to quantify methane emissions. Maps of time-averaged satellite-derived XCH4 show regionally elevated methane over several methane source regions. In order to obtain methane emissions of these source regions we use a simple and fast data-driven method to estimate annual methane emissions and corresponding 1σ uncertainties directly from maps of annually averaged satellite XCH4. From theoretical considerations we expect that our method tends to underestimate emissions. When applying our method to high-resolution atmospheric methane simulations, we typically find agreement within the uncertainty range of our method (often 100 %) but also find that our method tends to underestimate emissions by typically about 40 %. To what extent these findings are model dependent needs to be assessed. We apply our method to an ensemble of satellite XCH4 data products consisting of two products from SCIAMACHY/ENVISAT and two products from TANSO-FTS/GOSAT covering the time period 2003–2014. We obtain annual emissions of four source areas: Four Corners in the south-western USA, the southern part of Central Valley, California, Azerbaijan, and Turkmenistan. We find that our estimated emissions are in good agreement with independently derived estimates for Four Corners and Azerbaijan. For the Central Valley and Turkmenistan our estimated annual emissions are higher compared to the EDGAR v4.2 anthropogenic emission inventory. For Turkmenistan we find on average about 50 % higher emissions with our annual emission uncertainty estimates overlapping with the EDGAR emissions. For the region around Bakersfield in the Central Valley we find a factor of 5–8 higher emissions compared to EDGAR, albeit with large uncertainty. Major methane emission sources in this region are oil/gas and livestock. Our findings corroborate recently published studies based on aircraft and satellite measurements and new bottom-up estimates reporting significantly underestimated methane emissions of oil/gas and/or livestock in this area in EDGAR.

Consistent and accurate long-term data sets of global atmospheric concentrations of carbon dioxide (CO2) are required for carbon cycle and climate-related research. However, global data sets based on satellite observations may suffer from inconsistencies originating from the use of products derived from different satellites as needed to cover a long enough time period. One reason for inconsistencies can be the use of different retrieval algorithms. We address this potential issue by applying the same algorithm, the Bremen Optimal Estimation DOAS (BESD) algorithm, to different satellite instruments, SCIAMACHY on-board ENVISAT (March 2002–April 2012) and TANSO-FTS on-board GOSAT (launched in January 2009), to retrieve XCO2, the column-averaged dry-air mole fraction of CO2. BESD has been initially developed for SCIAMACHY XCO2 retrievals. Here, we present the first detailed assessment of the new GOSAT BESD XCO2 product. GOSAT BESD XCO2 is a product generated and delivered to the MACC project for assimilation into ECMWF's Integrated Forecasting System. We describe the modifications of the BESD algorithm needed in order to retrieve XCO2 from GOSAT and present detailed comparisons with ground-based observations of XCO2 from the Total Carbon Column Observing Network (TCCON). We discuss detailed comparison results between all three XCO2 data sets (SCIAMACHY, GOSAT and TCCON). The comparison results demonstrate the good consistency between SCIAMACHY and GOSAT XCO2. For example, we found a mean difference for daily averages of −0.60 ± 1.56 ppm (mean difference ± standard deviation) for GOSAT–SCIAMACHY (linear correlation coefficient r=0.82), −0.34 ± 1.37 ppm (r = 0.86) for GOSAT–TCCON and 0.10 ± 1.79 ppm (r = 0.75) for SCIAMACHY–TCCON. The remaining differences between GOSAT and SCIAMACHY are likely due to non-perfect collocation (± 2 h, 10° x 10° around TCCON sites), i.e. the observed air masses are not exactly identical but likely also due to a still non-perfect BESD retrieval algorithm, which will be continuously improved in the future. Our overarching goal is to generate a satellite-derived XCO2 data set appropriate for climate and carbon cycle research covering the longest possible time period. We therefore also plan to extend the existing SCIAMACHY and GOSAT data set discussed here by also using data from other missions (e.g. OCO-2, GOSAT-2, CarbonSat) in the future.

The GHG-CCI project (http://www.esa-ghg-cci.org/) is one of several projects of the European Space Agency’s (ESA) Climate Change Initiative (CCI). The goal of the CCI is to generate and deliver data sets of various satellite-derived Essential Climate Variables (ECVs) in line with GCOS (Global Climate Observing System) requirements. The “ECV Greenhouse Gases” (ECV GHG) is the global distribution of important climate relevant gases – namely atmospheric CO2 and CH4 - with a quality sufficient to obtain information on regional CO2 and CH4 sources and sinks. The main goal of GHG-CCI is to generate long-term highly accurate and precise time series of global near-surface-sensitive satellite observations of CO2 and CH4, i.e., XCO2 and XCH4, starting with the launch of ESA’s ENVISAT satellite. These products are currently retrieved from SCIAMACHY/ENVISAT (2002-2012) and TANSO-FTS/GOSAT (2009-today) nadir mode observations in the near-infrared/shortwave-infrared spectral region. In addition, other sensors (e.g., IASI and MIPAS) and viewing modes (e.g., SCIAMACHY solar occultation) are also considered and in the future also data from other satellites. The GHG-CCI data products and related documentation are freely available via the GHG-CCI website and yearly updates are foreseen. Here we present an overview about the latest data set (Climate Research Data Package No. 2 (CRDP#2)) and summarize key findings from using satellite CO2 and CH4 retrievals to improve our understanding of the natural and anthropogenic sources and sinks of these important atmospheric greenhouse gases. We also shortly mention ongoing activities related to validation and initial user assessment of CRDP#2 and future plans.

Wetlands are a major emission source of methane (CH4) globally. In this study, we evaluate wetland emission estimates derived using the UK community land surface model (JULES, the Joint UK Land Earth Simulator) against atmospheric observations of methane, including, for the first time, total methane columns derived from the SCIAMACHY instrument on board the ENVISAT satellite. Two JULES wetland emission estimates are investigated: (a) from an offline run driven with Climatic Research Unit–National Centers for Environmental Prediction (CRU-NCEP) meteorological data and (b) from the same offline run in which the modelled wetland fractions are replaced with those derived from the Global Inundation Extent from Multi-Satellites (GIEMS) remote sensing product. The mean annual emission assumed for each inventory (181 Tg CH4 per annum over the period 1999–2007) is in line with other recently published estimates. There are regional differences as the unconstrained JULES inventory gives significantly higher emissions in the Amazon (by ~36 Tg CH4 yr−1) and lower emissions in other regions (by up to 10 Tg CH4 yr−1) compared to the JULES estimates constrained with the GIEMS product. Using the UK Hadley Centre's Earth System model with atmospheric chemistry (HadGEM2), we evaluate these JULES wetland emissions against atmospheric observations of methane. We obtain improved agreement with the surface concentration measurements, especially at high northern latitudes, compared to previous HadGEM2 runs using the wetland emission data set of Fung et al. (1991). Although the modelled monthly atmospheric methane columns reproduce the large-scale patterns in the SCIAMACHY observations, they are biased low by 50 part per billion by volume (ppb). Replacing the HadGEM2 modelled concentrations above 300 hPa with HALOE–ACE assimilated TOMCAT output results in a significantly better agreement with the SCIAMACHY observations. The use of the GIEMS product to constrain the JULES-derived wetland fraction improves the representation of the wetland emissions in JULES and gives a good description of the seasonality observed at surface sites influenced by wetlands, especially at high latitudes. We find that the annual cycles observed in the SCIAMACHY measurements and at many of the surface sites influenced by non-wetland sources cannot be reproduced in these HadGEM2 runs. This suggests that the emissions over certain regions (e.g. India and China) are possibly too high and/or the monthly emission patterns for specific sectors are incorrect. The comparisons presented in this paper show that the performance of the JULES wetland scheme is comparable to that of other process-based land surface models. We identify areas for improvement in this and the atmospheric chemistry components of the HadGEM Earth System model. The Earth Observation data sets used here will be of continued value in future evaluations of JULES and the HadGEM family of models.

Current knowledge about the European terrestrial biospheric carbon sink, from the Atlantic to the Urals, relies upon bottom-up inventory and surface flux inverse model estimates (e.g. 0.27±0.16 GtC a−1 for 2000–2005 (Schulze et al., 2009), 0.17±0.44 GtC a−1 for 2001–2007 (Peters et al., 2010), 0.45±0.40 GtC a−1 for 2010 (Chevallier et al., 2014), 0.40±0.42 GtC a−1 for 2001–2004 (Peylin et al., 2013)). Inverse models assimilate in situ CO2 atmospheric concentrations measured by surface-based air sampling networks. The intrinsic sparseness of these networks is one reason for the relatively large flux uncertainties (Peters et al., 2010; Bruhwiler et al., 2011). Satellite-based CO2 measurements have the potential to reduce these uncertainties (Miller et al., 2007; Chevallier et al., 2007). Global inversion experiments using independent models and independent GOSAT satellite data products consistently derived a considerably larger European sink (1.0–1.3 GtC a−1 for 09/2009–08/2010 (Basu et al., 2013), 1.2–1.8 GtC a−1 in 2010 (Chevallier et al., 2014)). However, these results have been considered unrealistic due to potential retrieval biases and/or transport errors (Chevallier et al., 2014) or have not been discussed at all (Basu et al., 2013; Takagi et al., 2014). Our analysis comprises a regional inversion approach using STILT (Gerbig et al., 2003; Lin et al., 2003) short-range (days) particle dispersion modelling, rendering it insensitive to large-scale retrieval biases and less sensitive to long-range transport errors. We show that the satellite-derived European terrestrial carbon sink is indeed much larger (1.02±0.30 GtC a−1 in 2010) than previously expected. This is qualitatively consistent among an ensemble of five different inversion set-ups and five independent satellite retrievals (BESD (Reuter et al., 2011) 2003–2010, ACOS (O’Dell et al., 2012) 2010, UoL-FP (Cogan et al., 2012) 2010, RemoTeC (Butz et al., 2011) 2010, and NIES (Yoshida et al., 2013) 2010) using data from two different instruments (SCIAMACHY (Bovensmann et al., 1999) and GOSAT (Kuze et al., 2009)). The difference to in situ based inversions (Peylin et al., 2013), whilst large with respect to the mean reported European carbon sink (0.4 GtC a−1 for 2001–2004), is similar in magnitude to the reported uncertainty (0.42 GtC a−1). The highest gain in information is obtained during the growing season when satellite observation conditions are advantageous, a priori uncertainties are largest, and the surface sink maximises; during the dormant season, the results are dominated by the a priori. Our results provide evidence that the current understanding of the European carbon sink has to be revisited.

Column-averaged dry-air mole fractions of carbon dioxide and methane have been retrieved from spectra acquired by the TANSO-FTS (Thermal And Near-infrared Sensor for carbon Observations-Fourier Transform Spectrometer) and SCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric Cartography) instruments on board GOSAT (Greenhouse gases Observing SATellite) and ENVISAT (ENVIronmental SATellite), respectively, using a range of European retrieval algorithms. These retrievals have been compared with data from ground-based high-resolution Fourier transform spectrometers (FTSs) from the Total Carbon Column Observing Network (TCCON). The participating algorithms are the weighting function modified differential optical absorption spectroscopy (DOAS) algorithm (WFMD, University of Bremen), the Bremen optimal estimation DOAS algorithm (BESD, University of Bremen), the iterative maximum a posteriori DOAS (IMAP, Jet Propulsion Laboratory (JPL) and Netherlands Institute for Space Research algorithm (SRON)), the proxy and full-physics versions of SRON's RemoTeC algorithm (SRPR and SRFP, respectively) and the proxy and full-physics versions of the University of Leicester's adaptation of the OCO (Orbiting Carbon Observatory) algorithm (OCPR and OCFP, respectively). The goal of this algorithm inter-comparison was to identify strengths and weaknesses of the various so-called round- robin data sets generated with the various algorithms so as to determine which of the competing algorithms would proceed to the next round of the European Space Agency's (ESA) Greenhouse Gas Climate Change Initiative (GHG-CCI) project, which is the generation of the so-called Climate Research Data Package (CRDP), which is the first version of the Essential Climate Variable (ECV) "greenhouse gases" (GHGs). For XCO2, all algorithms reach the precision requirements for inverse modelling (< 8 ppm), with only WFMD having a lower precision (4.7 ppm) than the other algorithm products (2.4–2.5 ppm). When looking at the seasonal relative accuracy (SRA, variability of the bias in space and time), none of the algorithms have reached the demanding < 0.5 ppm threshold. For XCH4, the precision for both SCIAMACHY products (50.2 ppb for IMAP and 76.4 ppb for WFMD) fails to meet the < 34 ppb threshold for inverse modelling, but note that this work focusses on the period after the 2005 SCIAMACHY detector degradation. The GOSAT XCH4 precision ranges between 18.1 and 14.0 ppb. Looking at the SRA, all GOSAT algorithm products reach the < 10 ppm threshold (values ranging between 5.4 and 6.2 ppb). For SCIAMACHY, IMAP and WFMD have a SRA of 17.2 and 10.5 ppb, respectively.

The terrestrial biosphere is currently acting as a net carbon sink on the global scale, exhibiting significant interannual variability in strength. To reliably predict the future strength of the land sink and its role in atmospheric CO2 growth, the underlying biogeochemical processes and their response to a changing climate need to be well understood. In particular, better knowledge of the impact of key climate variables such as temperature or precipitation on the biospheric carbon reservoir is essential. It is demonstrated using nearly a decade of SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) nadir measurements that years with higher temperatures during the growing season can be robustly associated with larger growth rates in atmospheric CO2 and smaller seasonal cycle amplitudes for northern mid-latitudes. We find linear relationships between warming and CO2 growth as well as seasonal cycle amplitude at the 98% significance level. This suggests that the terrestrial carbon sink is less efficient at higher temperatures during the analysed time period. Unless the biosphere has the ability to adapt its carbon storage under warming conditions in the longer term, such a temperature response entails the risk of potential future sink saturation via a positive carbon-climate feedback. Quantitatively, the covariation between the annual CO2 growth rates derived from SCIAMACHY data and warm season surface temperature anomaly amounts to 1.25 ± 0.32 ppm yr−1 K−1 for the Northern Hemisphere, where the bulk of the terrestrial carbon sink is located. In comparison, this relationship is less pronounced in the Southern Hemisphere. The covariation of the seasonal cycle amplitudes retrieved from satellite measurements and temperature anomaly is −1.30 ± 0.31 ppm K−1 for the north temperate zone. These estimates are consistent with those from the CarbonTracker data assimilated CO2 data product, indicating that the temperature dependence of the model surface fluxes is realistic.

Urban areas, which are home to the majority of today's world population, are responsible for more than two-thirds of the global energy-related carbon dioxide emissions. Given the ongoing demographic growth and rising energy consumption in metropolitan regions particularly in the developing world, urban-based emissions are expected to further increase in the future. As a consequence, monitoring and independent verification of reported anthropogenic emissions is becoming more and more important. It is demonstrated using SCIAMACHY nadir measurements that anthropogenic CO2 emissions can be detected from space and that emission trends might be tracked using satellite observations. This is promising with regard to future satellite missions with high spatial resolution and wide swath imaging capability aiming at constraining anthropogenic emissions down to the point-source scale. By subtracting retrieved background values from those retrieved over urban areas we find significant CO2 enhancements for several anthropogenic source regions, namely 1.3 &pm; 0.7 ppm for the Rhine-Ruhr metropolitan region and the Benelux, 1.1 &pm; 0.5 ppm for the East Coast of the United States, and 2.4 &pm; 0.9 ppm for the Yangtze River Delta. The order of magnitude of the enhancements is in agreement with what is expected for anthropogenic CO2 signals. The larger standard deviation of the retrieved Yangtze River Delta enhancement is due to a distinct positive trend of 0.3 &pm; 0.2 ppm yr−1, which is quantitatively consistent with anthropogenic emissions from the Emission Database for Global Atmospheric Research (EDGAR) in terms of percentual increase per year. Potential contributions to the retrieved CO2 enhancement by several error sources, e.g. aerosols, albedo, and residual biospheric signals due to heterogeneous seasonal sampling, are discussed and can be ruled out to a large extent.

We analyze an ensemble of seven XCO2 retrieval algorithms for SCIAMACHY (scanning imaging absorption spectrometer of atmospheric chartography) and GOSAT (greenhouse gases observing satellite). The ensemble spread can be interpreted as regional uncertainty and can help to identify locations for new TCCON (total carbon column observing network) validation sites. Additionally, we introduce the ensemble median algorithm EMMA combining individual soundings of the seven algorithms into one new data set. The ensemble takes advantage of the algorithms' independent developments. We find ensemble spreads being often < 1 ppm but rising up to 2 ppm especially in the tropics and East Asia. On the basis of gridded monthly averages, we compare EMMA and all individual algorithms with TCCON and CarbonTracker model results (potential outliers, north/south gradient, seasonal (peak-to-peak) amplitude, standard deviation of the difference). Our findings show that EMMA is a promising candidate for inverse modeling studies. Compared to CarbonTracker, the satellite retrievals find consistently larger north/south gradients (by 0.3–0.9 ppm) and seasonal amplitudes (by 1.5–2.0 ppm).

Global observations of column-averaged dry air mole fractions of carbon dioxide (CO2), denoted by XCO2 , retrieved from SCIAMACHY on-board ENVISAT can provide important and missing global information on the distribution and magnitude of regional CO2 surface fluxes. This application has challenging precision and accuracy requirements. In a previous publication (Heymann et al., 2012), it has been shown by analysing seven years of SCIAMACHY WFM-DOAS XCO2 (WFMDv2.1) that unaccounted thin cirrus clouds can result in significant errors. In order to enhance the quality of the SCIAMACHY XCO2 data product, we have developed a new version of the retrieval algorithm (WFMDv2.2), which is described in this manuscript. It is based on an improved cloud filtering and correction method using the 1.4 μm strong water vapour absorption and 0.76 μm O2-A bands. The new algorithm has been used to generate a SCIAMACHY XCO2 data set covering the years 2003–2009. The new XCO2 data set has been validated using ground-based observations from the Total Carbon Column Observing Network (TCCON). The validation shows a significant improvement of the new product (v2.2) in comparison to the previous product (v2.1). For example, the standard deviation of the difference to TCCON at Darwin, Australia, has been reduced from 4 ppm to 2 ppm. The monthly regional-scale scatter of the data (defined as the mean intra-monthly standard deviation of all quality filtered XCO2 retrievals within a radius of 350 km around various locations) has also been reduced, typically by a factor of about 1.5. Overall, the validation of the new WFMDv2.2 XCO2 data product can be summarised by a single measurement precision of 3.8 ppm, an estimated regional-scale (radius of 500 km) precision of monthly averages of 1.6 ppm and an estimated regional-scale relative accuracy of 0.8 ppm. In addition to the comparison with the limited number of TCCON sites, we also present a comparison with NOAA's global CO2 modelling and assimilation system CarbonTracker. This comparison also shows significant improvements especially over the Southern Hemisphere.

Carbon dioxide (CO2) is the most important greenhouse gas whose atmospheric loading has been significantly increased by anthropogenic activity leading to global warming. Accurate measurements and models are needed in order to reliably predict our future climate. This, however, has challenging requirements. Errors in measurements and models need to be identified and minimised. In this context, we present a comparison between satellite-derived column-averaged dry air mole fractions of CO2, denoted XCO2, retrieved from SCIAMACHY/ENVISAT using the WFM-DOAS (weighting function modified differential optical absorption spectroscopy) algorithm, and output from NOAA's global CO2 modelling and assimilation system CarbonTracker. We investigate to what extent differences between these two data sets are influenced by systematic retrieval errors due to aerosols and unaccounted clouds. We analyse seven years of SCIAMACHY WFM-DOAS version 2.1 retrievals (WFMDv2.1) using CarbonTracker version 2010. We investigate to what extent the difference between SCIAMACHY and CarbonTracker XCO2 are temporally and spatially correlated with global aerosol and cloud data sets. For this purpose, we use a global aerosol data set generated within the European GEMS project, which is based on assimilated MODIS satellite data. For clouds, we use a data set derived from CALIOP/CALIPSO. We find significant correlations of the SCIAMACHY minus CarbonTracker XCO2 difference with thin clouds over the Southern Hemisphere. The maximum temporal correlation we find for Darwin, Australia (r2 = 54%). Large temporal correlations with thin clouds are also observed over other regions of the Southern Hemisphere (e.g. 43% for South America and 31% for South Africa). Over the Northern Hemisphere the temporal correlations are typically much lower. An exception is India, where large temporal correlations with clouds and aerosols have also been found. For all other regions the temporal correlations with aerosol are typically low. For the spatial correlations the picture is less clear. They are typically low for both aerosols and clouds, but depending on region and season, they may exceed 30% (the maximum value of 46% has been found for Darwin during September to November). Overall we find that the presence of thin clouds can potentially explain a significant fraction of the difference between SCIAMACHY WFMDv2.1 XCO2 and CarbonTracker over the Southern Hemisphere. Aerosols appear to be less of a problem. Our study indicates that the quality of the satellite derived XCO2 will significantly benefit from a reduction of scattering related retrieval errors at least for the Southern Hemisphere.

A simple empirical CO2 model (SECM) is presented to estimate column-average dry-air mole fractions of atmospheric CO2 (XCO2) as well as mixing ratio profiles. SECM is based on a simple equation depending on 17 empirical parameters, latitude, and date. The empirical parameters have been determined by least squares fitting to NOAA's (National Oceanic and Atmospheric Administration) assimilation system CarbonTracker version 2010 (CT2010). Comparisons with TCCON (total carbon column observing network) FTS (Fourier transform spectrometer) measurements show that SECM XCO2 agrees quite well with reality. The synthetic XCO2 values have a standard error of 1.39 ppm and systematic station-to-station biases of 0.46 ppm. Typical column averaging kernels of the TCCON FTS, a SCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric CHartographY), and two GOSAT (Greenhouse gases Observing SATellite) XCO2 retrieval algorithms have been used to assess the smoothing error introduced by using SECM profiles instead of CT2010 profiles as a priori. The additional smoothing error amounts to 0.17 ppm for a typical SCIAMACHY averaging kernel and is most times much smaller for the other instruments (e.g. 0.05 ppm for a typical TCCON FTS averaging kernel). Therefore, SECM is well suited to provide a priori information for state-of-the-art ground-based (FTS) and satellite-based (GOSAT, SCIAMACHY) XCO2 retrievals. Other potential applications are: (i) near real-time processing systems (that cannot make use of models like CT2010 operated in delayed mode), (ii) "CO2 proxy" methods for XCH4 retrievals (as correction for the XCO2 background), and (iii) observing system simulation experiments especially for future satellite missions.

SCIAMACHY onboard ENVISAT (launched in 2002) enables the retrieval of global long-term column-averaged dry air mole fractions of the two most important anthropogenic greenhouse gases carbon dioxide and methane (denoted XCO2 and XCH4). In order to assess the quality of the greenhouse gas data obtained with the recently introduced v2 of the scientific retrieval algorithm WFM-DOAS, we present validations with ground-based Fourier Transform Spectrometer (FTS) measurements and comparisons with model results at eight Total Carbon Column Observing Network (TCCON) sites providing realistic error estimates of the satellite data. Such validation is a prerequisite to assess the suitability of data sets for their use in inverse modelling. It is shown that there are generally no significant differences between the carbon dioxide annual increases of SCIAMACHY and the assimilation system CarbonTracker (2.00 &pm; 0.16 ppm yr−1 compared to 1.94 &pm; 0.03 ppm yr−1 on global average). The XCO2 seasonal cycle amplitudes derived from SCIAMACHY are typically larger than those from TCCON which are in turn larger than those from CarbonTracker. The absolute values of the northern hemispheric TCCON seasonal cycle amplitudes are closer to SCIAMACHY than to CarbonTracker and the corresponding differences are not significant when compared with SCIAMACHY, whereas they can be significant for a subset of the analysed TCCON sites when compared with CarbonTracker. At Darwin we find discrepancies of the seasonal cycle derived from SCIAMACHY compared to the other data sets which can probably be ascribed to occurrences of undetected thin clouds. Based on the comparison with the reference data, we conclude that the carbon dioxide data set can be characterised by a regional relative precision (mean standard deviation of the differences) of about 2.2 ppm and a relative accuracy (standard deviation of the mean differences) of 1.1–1.2 ppm for monthly average composites within a radius of 500 km. For methane, prior to November 2005, the regional relative precision amounts to 12 ppb and the relative accuracy is about 3 ppb for monthly composite averages within the same radius. The loss of some spectral detector pixels results in a degradation of performance thereafter in the spectral range currently used for the methane column retrieval. This leads to larger scatter and lower XCH4 values are retrieved in the tropics for the subsequent time period degrading the relative accuracy. As a result, the overall relative precision is estimated to be 17 ppb and the relative accuracy is in the range of about 10–20 ppb for monthly averages within a radius of 500 km. The derived estimates show that the SCIAMACHY XCH4 data set before November 2005 is suitable for regional source/sink determination and regional-scale flux uncertainty reduction via inverse modelling worldwide. In addition, the XCO2 monthly data potentially provide valuable information in continental regions, where there is sparse sampling by surface flask measurements.

Carbon dioxide (CO2) is the most important man-made greenhouse gas (GHG) that cause global warming. With electricity generation through fossil-fuel power plants now being the economic sector with the largest source of CO2, power plant emissions monitoring has become more important than ever in the fight against global warming. In a previous study done by Bovensmann et al. (2010), random and systematic errors of power plant CO2 emissions have been quantified using a single overpass from a proposed CarbonSat instrument. In this study, we quantify errors of power plant annual emission estimates from a hypothetical CarbonSat and constellations of several CarbonSats while taking into account that power plant CO2 emissions are time-dependent. Our focus is on estimating systematic errors arising from the sparse temporal sampling as well as random errors that are primarily dependent on wind speeds. We used hourly emissions data from the US Environmental Protection Agency (EPA) combined with assimilated and re-analyzed meteorological fields from the National Centers of Environmental Prediction (NCEP). CarbonSat orbits were simulated as a sun-synchronous low-earth orbiting satellite (LEO) with an 828-km orbit height, local time ascending node (LTAN) of 13:30 (01:30 p.m. LT) and achieves global coverage after 5 days. We show, that despite the variability of the power plant emissions and the limited satellite overpasses, one CarbonSat has the potential to verify reported US annual CO2 emissions from large power plants (≥5 Mt CO2 yr−1) with a systematic error of less than ~4.9% and a random error of less than ~6.7% for 50% of all the power plants. For 90% of all the power plants, the systematic error was less than ~12.4% and the random error was less than ~13%. We additionally investigated two different satellite configurations using a combination of 5 CarbonSats. One achieves global coverage everyday but only samples the targets at fixed local times. The other configuration samples the targets five times at two-hour intervals approximately every 6th day but only achieves global coverage after 5 days. From the statistical analyses, we found, as expected, that the random errors improve by approximately a factor of two if 5 satellites are used. On the other hand, more satellites do not result in a large reduction of the systematic error. The systematic error is somewhat smaller for the CarbonSat constellation configuration achieving global coverage everyday. Therefore, we recommend the CarbonSat constellation configuration that achieves daily global coverage.

Carbon dioxide (CO2) and methane (CH4) are the two most important anthropogenic greenhouse gases contributing to global climate change. SCIAMACHY onboard ENVISAT (launch 2002) was the first and is now with TANSO onboard GOSAT (launch 2009) one of only two satellite instruments currently in space whose measurements are sensitive to CO2 and CH4 concentration changes in the lowest atmospheric layers where the variability due to sources and sinks is largest. We present long-term SCIAMACHY retrievals (2003–2009) of column-averaged dry air mole fractions of both gases (denoted XCO2 and XCH4) derived from absorption bands in the near-infrared/shortwave-infrared (NIR/SWIR) spectral region focusing on large-scale features. The results are obtained using an upgraded version (v2) of the retrieval algorithm WFM-DOAS including several improvements, while simultaneously maintaining its high processing speed. The retrieved mole fractions are compared to global model simulations (CarbonTracker XCO2 and TM5 XCH4) being optimised by assimilating highly accurate surface measurements from the NOAA/ESRL network and taking the SCIAMACHY averaging kernels into account. The comparisons address seasonal variations and long-term characteristics. The steady increase of atmospheric carbon dioxide primarily caused by the burning of fossil fuels can be clearly observed with SCIAMACHY globally. The retrieved global annual mean XCO2 increase agrees with CarbonTracker within the error bars (1.80±0.13 ppm yr−1 compared to 1.81±0.09 ppm yr−1). The amplitude of the XCO2 seasonal cycle as retrieved by SCIAMACHY, which is 4.3±0.2 ppm for the Northern Hemisphere and 1.4±0.2 ppm for the Southern Hemisphere, is on average about 1 ppm larger than for CarbonTracker. An investigation of the boreal forest carbon uptake during the growing season via the analysis of longitudinal gradients shows good agreement between SCIAMACHY and CarbonTracker concerning the overall magnitude of the gradients and their annual variations. The analysis includes a discussion of the relative uptake strengths of the Russian and North American boreal forest regions. The retrieved XCH4 results show that after years of stability, atmospheric methane has started to rise again in recent years which is consistent with surface measurements. The largest increase is observed for the tropics and northern mid- and high-latitudes amounting to about 7.5±1.5 ppb yr−1 since 2007. Due care has been exercised to minimise the influence of detector degradation on the quantitative estimate of this anomaly.

Total column concentrations of methane have been retrieved from ground-based solar absorption FTIR spectra in the near-infrared recorded in Paramaribo (Suriname). The methane FTIR observations are compared with TM5 model simulations and satellite observations from SCIAMACHY, and represent the first validation of SCIAMACHY retrievals in the inner tropics using ground-based remote sensing techniques. Apart from local biomass burning features, our methane FTIR observations agree well with TM5 model simulations. The comparison of the direct measured CH4/CO2 ratios by FTIR and satellite reveals that the satellite can hardly detect methane emissions of tropical biomass burning due to the used retrieval method.

Carbon dioxide (CO2) is the most important anthropogenic greenhouse gas (GHG) causing global warming. The atmospheric CO2 concentration increased by more than 30% since pre-industrial times – primarily due to burning of fossil fuels – and still continues to increase. Reporting of CO2 emissions is required by the Kyoto protocol. Independent verification of reported emissions, which are typially not directly measured, by methods such as inverse modeling of measured atmospheric CO2 concentrations is currently not possible globally due to lack of appropriate observations. Existing satellite instruments such as SCIAMACHY/ENVISAT and TANSO/GOSAT focus on advancing our understanding of natural CO2 sources and sinks. The obvious next step for future generation satellites is to also constrain anthropogenic CO2 emissions. Here we present a promising satellite remote sensing concept based on spectroscopic measurements of reflected solar radiation and show, using power plants as an example, that strong localized CO2 point sources can be detected and their emissions quantified. This requires mapping the atmospheric CO2 column distribution at a spatial resolution of 2×2 km2 with a precision of 0.5% (2 ppm) or better. We indicate that this can be achieved with existing technology. For a single satellite in sun-synchronous orbit with a swath width of 500 km, each power plant (PP) is overflown every 6 days or more frequent. Based on the MODIS cloud mask data product we conservatively estimate that typically 20 sufficiently cloud free overpasses per PP can be achieved every year. We found that for typical wind speeds in the range of 2–6 m/s the statistical uncertainty of the retrieved PP CO2 emission due to instrument noise is in the range 1.6–4.8 MtCO2/yr for single overpasses. This corresponds to 12–36% of the emission of a mid-size PP (13 MtCO2/yr). We have also determined the sensitivity to parameters which may result in systematic errors such as atmospheric transport and aerosol related parameters. We found that the emission error depends linearly on wind speed, i.e., a 10% wind speed error results in a 10% emission error, and that neglecting enhanced aerosol concentrations in the PP plume may result in errors in the range 0.2–2.5 MtCO2/yr, depending on PP aerosol emission. The discussed concept has the potential to contribute to an independent verification of reported anthropogenic CO2 emissions and therefore could be an important component of a future global anthropogenic GHG emission monitoring system. This is of relevance in the context of Kyoto protocol follow-on agreements but also allows detection and monitoring of a variety of other strong natural and anthropogenic CO2 and CH4 emitters. The investigated instrument is not limited to these applications as it has been specified to also deliver the data needed for global regional-scale CO2 and CH4 surface flux inverse modeling.

An optimal estimation based retrieval scheme for satellite based retrievals of XCO2 (the dry air column averaged mixing ratio of atmospheric CO2) is presented enabling accurate retrievals also in the presence of thin clouds. The proposed method is designed to analyze near-infrared nadir measurements of the SCIAMACHY instrument in the CO2 absorption band at 1580 nm and in the O2-A absorption band at around 760 nm. The algorithm accounts for scattering in an optically thin cirrus cloud layer and at aerosols of a default profile. The scattering information is mainly obtained from the O2-A band and a merged fit windows approach enables the transfer of information between the O2-A and the CO2 band. Via the optimal estimation technique, the algorithm is able to account for a priori information to further constrain the inversion. Test scenarios of simulated SCIAMACHY sun-normalized radiance measurements are analyzed in order to specify the quality of the proposed method. In contrast to existing algorithms for SCIAMACHY retrievals, the systematic errors due to cirrus clouds with optical thicknesses up to 1.0 are reduced to values below 4 ppm for most of the analyzed scenarios. This shows that the proposed method has the potential to reduce uncertainties of SCIAMACHY retrieved XCO2 making this data product potentially useful for surface flux inverse modeling.

Carbon dioxide (CO2) and methane (CH4) are the two most important anthropogenic greenhouse gases. SCIAMACHY on ENVISAT is the first satellite instrument whose measurements are sensitive to concentration changes of the two gases at all altitude levels down to the Earth's surface where the source/sink signals are largest. We have processed three years (2003–2005) of SCIAMACHY near-infrared nadir measurements to simultaneously retrieve vertical columns of CO2 (from the 1.58 μm absorption band), CH4 (1.66 μm) and oxygen (O2 A-band at 0.76 μm) using the scientific retrieval algorithm WFM-DOAS. We show that the latest version of WFM-DOAS, version 1.0, which is used for this study, has been significantly improved with respect to its accuracy compared to the previous versions while essentially maintaining its high processing speed (~1 min per orbit, corresponding to ~6000 single measurements, and per gas on a standard PC). The greenhouse gas columns are converted to dry air column-averaged mole fractions, denoted XCO2 (in ppm) and XCH4 (in ppb), by dividing the greenhouse gas columns by simultaneously retrieved dry air columns. For XCO2 dry air columns are obtained from the retrieved O2 columns. For XCH4 dry air columns are obtained from the retrieved CO2 columns because of better cancellation of light path related errors compared to using O2 columns retrieved from the spectrally distant O2 A-band. Here we focus on a discussion of the XCH4 data set. The XCO2 data set is discussed in a separate paper (Part 1). For 2003 we present detailed comparisons with the TM5 model which has been optimally matched to highly accurate but sparse methane surface observations. After accounting for a systematic low bias of ~2% agreement with TM5 is typically within 1–2%. We investigated to what extent the SCIAMACHY XCH4 is influenced by the variability of atmospheric CO2 using global CO2 fields from NOAA's CO2 assimilation system CarbonTracker. We show that the CO2 corrected and uncorrected XCH4 spatio-temporal pattern are very similar but that agreement with TM5 is better for the CarbonTracker CO2 corrected XCH4. In line with previous studies (e.g., Frankenberg et al., 2005b) we find higher methane over the tropics compared to the model. We show that tropical methane is also higher when normalizing the CH4 columns with retrieved O2 columns instead of CO2. In consistency with recent results of Frankenberg et al. (2008b) it is shown that the magnitude of the retrieved tropical methane is sensitive to the choice of the spectroscopic line parameters of water vapour. Concerning inter-annual variability we find similar methane spatio-temporal pattern for 2003 and 2004. For 2005 the retrieved methane shows significantly higher variability compared to the two previous years, most likely due to somewhat larger noise of the spectral measurements.

Carbon dioxide (CO2) and methane (CH4) are the two most important anthropogenic greenhouse gases. SCIAMACHY on ENVISAT is the first satellite instrument whose measurements are sensitive to concentration changes of the two gases at all altitude levels down to the Earth's surface where the source/sink signals are largest. We have processed three years (2003–2005) of SCIAMACHY near-infrared nadir measurements to simultaneously retrieve vertical columns of CO2 (from the 1.58 μm absorption band), CH4 (1.66 μm) and oxygen (O2 A-band at 0.76 μm) using the scientific retrieval algorithm WFM-DOAS. We show that the latest version of WFM-DOAS, version 1.0, which is used for this study, has been significantly improved with respect to its accuracy compared to the previous versions while essentially maintaining its high processing speed (~1 min per orbit, corresponding to ~6000 single measurements, and per gas on a standard PC). The greenhouse gas columns are converted to dry air column-averaged mole fractions, denoted XCO2 (in ppm) and XCH4 (in ppb), by dividing the greenhouse gas columns by simultaneously retrieved dry air columns. For XCO2 dry air columns are obtained from the retrieved O2 columns. For XCH4 dry air columns are obtained from the retrieved CO2 columns because of better cancellation of light path related errors compared to using O2 columns retrieved from the spectrally distant O2 A-band. Here we focus on a discussion of the XCO2 data set. The XCH4 data set is discussed in a separate paper (Part 2). In order to assess the quality of the retrieved XCO2 we present comparisons with Fourier Transform Spectroscopy (FTS) XCO2 measurements at two northern hemispheric mid-latitude ground stations. To assess the quality globally, we present detailed comparisons with global XCO2 fields obtained from NOAA's CO2 assimilation system CarbonTracker. For the Northern Hemisphere we find good agreement with the reference data for the CO2 seasonal cycle and the CO2 annual increase. For the Southern Hemisphere, where significantly less data are available for averaging compared to the Northern Hemisphere, the CO2 annual increase is also in good agreement with CarbonTracker but the amplitude and phase of the seasonal cycle show systematic differences (up to several ppm) arising partially from the O2 normalization most likely caused by unconsidered scattering effects due to subvisual cirrus clouds. The retrieved XCO2 regional pattern at monthly resolution over various regions show clear correlations with CarbonTracker but also significant differences. Typically the retrieved variability is about 4 ppm (1% of 380 ppm) higher but depending on time and location differences can reach or even exceed 8 ppm. Based on the error analysis and on the comparison with the reference data we conclude that the XCO2 data set can be characterized by a single measurement retrieval precision (random error) of 1–2%, a systematic low bias of about 1.5%, and by a relative accuracy of about 1–2% for monthly averages at a spatial resolution of about 7°×7°. When averaging the SCIAMACHY XCO2 over all three years we find elevated CO2 over the highly populated region of western central Germany and parts of the Netherlands ("Rhine-Main area") reasonably well correlated with EDGAR anthropogenic CO2 emissions. On average the regional enhancement is 2.7 ppm including an estimated contribution of 1–1.5 ppm due to aerosol related errors and sampling.

The reliable prediction of future atmospheric CO2 concentrations and associated global climate change requires an adequate understanding of the CO2 sources and sinks. The sparseness of the existing surface measurement network limits current knowledge about the global distribution of CO2 surface fluxes. The retrieval of CO2 total vertical columns from satellite observations is predicted to improve this situation. Such an application however requires very high accuracy and precision. We report on retrievals of the column-averaged CO2 dry air mole fraction, denoted XCO2, from the near-infrared nadir spectral radiance and solar irradiance measurements of the SCIAMACHY satellite instrument between 2003 and 2005. We focus on northern hemispheric large scale CO2 features such as the CO2 seasonal cycle and show - for the first time - that the atmospheric annual increase of CO2 can be directly observed using satellite measurements of the CO2 total column. The satellite retrievals are compared with global XCO2 obtained from NOAA's CO2 assimilation system CarbonTracker taking into account the spatio-temporal sampling and altitude sensitivity of the satellite data. We show that the measured CO2 year-to-year increase agrees within about 1 ppm/year with CarbonTracker. We also show that the latitude dependent amplitude of the northern hemispheric CO2 seasonal cycle agrees with CarbonTracker within about 2 ppm with the retrieved amplitude being systematically larger. The analysis demonstrates that it is possible using satellite measurements of the CO2 total column to retrieve information on the atmospheric CO2 on the level of a few parts per million.

The three carbon gases carbon monoxide (CO), carbon dioxide (CO2), and methane (CH4) are important atmospheric constituents affecting air quality and climate. The near-infrared nadir spectra measured by SCIAMACHY on ENVISAT contain information on the vertical columns of these gases which we retrieve using a modified DOAS algorithm (WFM-DOAS or WFMD). Our main data products are CO vertical columns and dry-air column averaged mixing ratios of methane (CH4) and CO2 (denoted XCH4 and XCO2). For CO and CH4 we present new results for the year 2003 obtained with an improved version of WFM-DOAS (WFMDv0.5) retrieved from Level 1 version 4 (Lv1v4) spectra. This data set has recently been compared with a network of ground based FTIR stations. Here we describe the WFMDv0.5 algorithm, present global and regional maps, and comparisons with global reference data. We show that major problems of the previous versions (v0.4 and v0.41) related to the varying ice-layer on the SCIAMACHY channel 8 detector have been solved. Compared to MOPITT the SCIAMACHY CO columns are on average higher by about 10–20%. Regionally, however, especially over central South America, differences can be much larger. For methane we present global and regional maps which are compared to TM5 model simulations performed using standard methane emission inventories. We show that methane source regions can be clearly detected with SCIAMACHY. We also show that the methane data product can be significantly further improved using Lv1v5 spectra with improved calibration. For CO2 we present three years of SCIAMACHY CO2 measurements over Park Falls, Wisconsin, USA, retrieved from Lv1v5. We show that the quality of CO2 retrieved from these spectra is significantly higher compared to WFMDv0.4 XCO2 retrieved from Lv1v4.

No abstract available.

Emissions of methane (CH 4 ) in the Permian basin (USA) have been derived for 2019 and 2020 from satellite observations of the Tropospheric Monitoring Instrument (TROPOMI) using the divergence method, in combination with a data driven method to estimate the background column densities. The resulting CH 4 emission data, which have been verified using model data with known emissions, have a spatial resolution of approximately 10 km. The CH 4 emissions show moderate spatial correlation with the locations of oil and gas production and drilling activities in the Permian basin, as well as with emissions of nitrogen oxides (NO x ). Analysis of the emission maps and time series indicates that a significant fraction of methane emissions in the Permian basin is from frequent widespread emissions sources, rather than from a few infrequent very large unplanned releases, which is important considering possible CH 4 emission mitigation strategies. In addition to providing spatially resolved emissions, the divergence method also provides the total emissions of the Permian basin and its main sub-basins. The total CH 4 emission of the Permian is estimated as 3.0 ± 0.7 Tg yr -1 for 2019, which agrees with other independent estimates based on TROPOMI data. For the Delaware sub-basin, it is estimated as 1.4 ± 0.3 Tg yr -1 for 2019, and for the Midland sub-basin 1.2 ± 0.3 Tg yr -1 . In 2020 the emissions are 9% lower compared to 2019 in the entire Permian basin, and respectively 19% and 27% for the Delaware and Midland sub-basins., (© 2023. The Authors.)