Search

Venezuela has long been identified as an area with large methane emissions and intensive oil exploitation, especially in the Lake Maracaibo region, but production has strongly decreased in recent years. The area is notoriously difficult to observe from space due to its complex topography and persistent cloud cover. We use the unprecedented coverage of the TROPOspheric Monitoring Instrument (TROPOMI) methane observations in analytical inversions with the Integrated Methane Inversion (IMI) framework at the national scale and at the local scale with the Weather Research and Forecasting model with chemistry (WRF-Chem). In the IMI analysis, we find Venezuelan emissions of 7.5 (5.7–9.3) Tg a−1 in 2019, where about half of emissions can be informed by TROPOMI observations, and emissions from oil exploitation are a factor of ∼ 1.6 higher than in bottom-up inventories. Using WRF, we find emissions of 1.2 (1.0–1.5) Tg a−1 from the Lake Maracaibo area in 2019, close to bottom-up estimates. Our WRF estimate is ∼ 40 % lower than the result over the same region from the IMI due to differences in the meteorology used by the two models. We find only a small, non-significant trend in emissions between 2018 and 2020 around the lake, implying the area's methane emission intensity expressed against oil and gas production has doubled over the time period, to ∼ 20 %. This value is much higher than what has previously been found for other oil and gas production regions and indicates that there could be large emissions from abandoned infrastructure.

Abstract Various fire emission estimates for southern Africa during 2019, derived with multiple burned area data sets with resolutions ranging from 500 to 20 m, are evaluated using satellite carbon monoxide (CO) observations. Southern African emissions derived from burned area generated by 20 m Sentinel‐2 satellite imagery are up to 120% higher than other estimates because small fires are better detected with a higher‐resolution satellite instrument. A comprehensive comparison between simulated and observed atmospheric CO indicates that the Sentinel‐2 burned area data significantly improves emission estimates, with up to 15% reduction in CO concentration biases in comparison to emissions based on coarser resolution burned area data. We also found that the temporal lag between emissions and atmospheric CO concentrations during the peak fire month was related to atmospheric transport. These findings emphasize the importance of utilizing higher‐resolution satellite instruments in accurately estimating emissions and understanding the impact of small fires on global climate.

A reduction in anthropogenic methane emissions is vital to limit near-term global warming. A small number of so-called super-emitters is responsible for a disproportionally large fraction of total methane emissions. Since late 2017, the TROPOspheric Monitoring Instrument (TROPOMI) has been in orbit, providing daily global coverage of methane mixing ratios at a resolution of up to 7×5.5 km2, enabling the detection of these super-emitters. However, TROPOMI produces millions of observations each day, which together with the complexity of the methane data, makes manual inspection infeasible. We have therefore designed a two-step machine learning approach using a convolutional neural network to detect plume-like structures in the methane data and subsequently apply a support vector classifier to distinguish the emission plumes from retrieval artifacts. The models are trained on pre-2021 data and subsequently applied to all 2021 observations. We detect 2974 plumes in 2021, with a mean estimated source rate of 44 t h−1 and 5–95th percentile range of 8–122 t h−1. These emissions originate from 94 persistent emission clusters and hundreds of transient sources. Based on bottom-up emission inventories, we find that most detected plumes are related to urban areas and/or landfills (35 %), followed by plumes from gas infrastructure (24 %), oil infrastructure (21 %), and coal mines (20 %). For 12 (clusters of) TROPOMI detections, we tip and cue the targeted observations and analysis of high-resolution satellite instruments to identify the exact sources responsible for these plumes. Using high-resolution observations from GHGSat, PRISMA, and Sentinel-2, we detect and analyze both persistent and transient facility-level emissions underlying the TROPOMI detections. We find emissions from landfills and fossil fuel exploitation facilities, and for the latter, we find up to 10 facilities contributing to one TROPOMI detection. Our automated TROPOMI-based monitoring system in combination with high-resolution satellite data allows for the detection, precise identification, and monitoring of these methane super-emitters, which is essential for mitigating their emissions.

This paper presents the automated plume detection and emission estimation algorithm (APE), developed to detect CO plumes from isolated biomass burning events and to quantify the corresponding CO emission rate. APE uses the CO product of the Tropospheric Monitoring Instrument (TROPOMI) on board the Copernicus Sentinel-5 Precursor (S5P) satellite, launched in 2017, and collocated active fire data from the Visible Infrared Imaging Radiometer Suite (VIIRS), the latter flying 3 min ahead of S5P. After identifying appropriate fire events using VIIRS data, an automated plume detection algorithm based on traditional image processing algorithms selects plumes for further data interpretation. The approach is based on thresholds optimized for data over the United States in September 2020. Subsequently, the CO emission rate is estimated using the cross-sectional flux method, which requires horizontal wind fields at the plume height. Three different plume heights were considered, and the ECMWF Reanalysis v5 (ERA5) data were used to compute emissions. A varying plume height in the downwind direction based on three-dimensional Lagrangian simulation was considered appropriate. APE is verified for observations over Australia and Siberia. For all fire sources identified by VIIRS, only 16 % of the data corresponded to clear-sky TROPOMI CO data with plume signature. Furthermore, the quality filters of APE resulted in emission estimations for 26 % of the TROPOMI CO data with plume signatures. Visual filtering of the APE's output showed a true-positive confidence level of 97.7 %. Finally, we provide an estimate of the emission uncertainties. The greatest contribution of error comes from the uncertainty in Global Fire Assimilation System (GFAS) injection height that leads to emission errors %, followed by systematic errors in the ERA5 wind data. The assumption of constant emission during plume formation and spatial under-sampling of CO column concentration by TROPOMI yields an error of %. The randomized errors from the ensemble ERA5 wind data are found to be less than 20 % for 97 % of the cases.

Carbon monoxide (CO) is an air pollutant that plays an important role in atmospheric chemistry and is mostly emitted by forest fires and incomplete combustion in, for example, road transport, residential heating, and industry. As CO is co-emitted with fossil fuel CO2 combustion emissions, it can be used as a proxy for CO2. Following the Paris Agreement, there is a need for independent verification of reported activity-based bottom-up CO2 emissions through atmospheric measurements. CO can be observed daily at a global scale with the TROPOspheric Monitoring Instrument (TROPOMI) satellite instrument with daily global coverage at a resolution down to 5.5 × 7 km2. To take advantage of this unique TROPOMI dataset, we develop a cross-sectional flux-based emission quantification method that can be applied to quantify emissions from a large number of cities, without relying on computationally expensive inversions. We focus on Africa as a region with quickly growing cities and large uncertainties in current emission estimates. We use a full year of high-resolution Weather Research and Forecasting (WRF) simulations over three cities to evaluate and optimize the performance of our cross-sectional flux emission quantification method and show its reliability down to emission rates of 0.1 Tg CO yr−1. Comparison of the TROPOMI-based emission estimates to the Dynamics–Aerosol–Chemistry–Cloud Interactions in West Africa (DACCIWA) and Emissions Database for Global Atmospheric Research (EDGAR) bottom-up inventories shows that CO emission rates in northern Africa are underestimated in EDGAR, suggesting overestimated combustion efficiencies. We see the opposite when comparing TROPOMI to the DACCIWA inventory in South Africa and Côte d'Ivoire, where CO emission factors appear to be overestimated. Over Lagos and Kano (Nigeria) we find that potential errors in the spatial disaggregation of national emissions cause errors in DACCIWA and EDGAR respectively. Finally, we show that our computationally efficient quantification method combined with the daily TROPOMI observations can identify a weekend effect in the road-transport-dominated CO emissions from Cairo and Algiers.

We quantify weekly methane emissions at 0.25∘ × 0.3125∘ (≈25 × 25 km2) resolution from the Permian Basin, the largest oil production basin in the US, by inverse analysis of satellite observations from the TROPOspheric Monitoring Instrument (TROPOMI) from May 2018 to October 2020. The mean oil and gas emission from the region (± standard deviation of weekly estimates) was 3.7 ± 0.9 Tg a−1, higher than previous TROPOMI inversion estimates that may have used biased prior emissions or background assumptions. We find strong week-to-week variability in emissions superimposed on longer-term trends, and these are consistent with independent inferences of temporal emission variability from tower, aircraft, and multispectral satellite data. New well development and natural gas spot price were significant drivers of variability in emissions over our study period but the concurrent 50 % increase in oil and gas production was not. The methane intensity (methane emitted per unit of methane gas produced) averaged 4.6 % ± 1.3 % and steadily decreased from 5 %–6 % in 2018 to 3 %–4 % in 2020. While the decreasing trend suggests improvement in operator practices during the study period, methane emissions from the Permian Basin remained high, with methane intensity an order of magnitude above the industry target of <0.2 %. Our success in using TROPOMI satellite observations for weekly estimates of emissions from a major oil production basin shows promise for application to near-real-time monitoring in support of climate change mitigation efforts.

Since 2017 the Tropospheric Monitoring Instrument (TROPOMI) on board ESA's Copernicus Sentinel-5 satellite (S5-P) has provided the operational carbon monoxide (CO) data product with daily global coverage on a spatial resolution of 5.5×7 km2. The European Centre for Medium-Range Weather Forecasts (ECMWF) plans to assimilate the retrieved total columns and the corresponding vertical sensitivities in the Copernicus Atmosphere Monitoring Service Integrated Forecasting System (CAMS-IFS) to improve forecasts of the atmospheric chemical composition. The TROPOMI data will primarily constrain the vertical integrated CO field of CAMS-IFS but to a lesser extent also its vertical CO distribution. For clear-sky conditions, the vertical sensitivity of the TROPOMI CO data product is useful throughout the atmosphere, but for cloudy scenes it varies due to cloud shielding and light scattering. To assess the profile information, we deploy an a posteriori profile retrieval that combines an ensemble of TROPOMI CO column retrievals with different vertical sensitivities to obtain a vertical CO profile that is then a representative average for the chosen spatial and temporal domain. We demonstrate the approach on three CO pollution cases. For the so-called “Rabbit Foot Fire” in Idaho on 12 August 2018, we estimate a CO profile showing the pollution at an altitude of about 5 km in good agreement with airborne in situ measurements of the Biomass Burning Flux Measurements of Trace Gases and Aerosol (BB-FLUX) field campaign. The distinct CO enhancement in a plume aloft (length =212 km, width =34 km), decoupled from the ground, is sensed by TROPOMI but is not present in the CAMS-IFS model. For a large-scale event, we analyzed the CO pollution from Siberian wildfires that took place from 14 to 17 August 2018. The TROPOMI data estimate the height of the pollution plume over Canada at 7 km in agreement with CAMS-IFS. However, CAMS-IFS underestimates the enhanced CO vertical column densities sensed by TROPOMI within the plume by more than 100 ppb. Finally, we study the seasonal biomass burning in the Amazon. During the burning season the CO profile retrieved from the TROPOMI measurements (1–15 August 2019) agrees well with the one of CAMS-IFS with a similar vertical shape between ground and 14 km altitude. Hence, our results indicate that assimilating TROPOMI CO retrieval with different vertical sensitivities (e.g., under clear-sky and cloudy conditions) provides information about the vertical distribution of CO.

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.

Satellite remote sensing of methane (CH4) using the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Copernicus Sentinel-5 Precursor (S5-P) satellite is key to monitor and quantify emissions globally. Overall, the S5-P methane data are of satisfying quality, demonstrated by the product validation with ground-based measurements from the Total Carbon Column Observing Network (TCCON). However, analysis of TROPOMI dry-air column mixing ratio (XCH4) data has pointed to a few false methane anomalies up to 20–40 ppb that can potentially be misinterpreted as enhancements due to strong emission sources. These artefacts are caused by spectral features of the underlying surfaces of specific materials (e.g. carbonate rocks), which are not well represented in the forward model of the retrieval algorithm. In this study we show that the observed anomalies are due to the surface model which describes the spectral dependence of the Lambertian albedo by a second-order polynomial in wavelength. By analysing the ECOSTRESS library that contains laboratory spectra for different types of surfaces, we find that a quadratic function might not be the most optimal representation of the surface reflectance spectral dependencies in the short-wave infrared (SWIR) range. Already the use of a third-order polynomial improves the methane data such that the anomalies disappear at several locations (e.g. Siberia, Australia and Algeria) without affecting the data quality elsewhere, and the quality of the fit significantly improves. We also found that the known bias in retrieved methane for low-albedo scenes slightly improves, but still, a posterior correction needs to be applied, leaving open the question about the root cause of the albedo bias. After applying the adjusted surface model globally, we perform the routine validation with TCCON and Greenhouse gases Observing SATellite (GOSAT) data. GOSAT comparison does not significantly improve, while TCCON validation results show a small improvement in some stations of 2–4 ppb, up to a factor of 10 smaller than the artificial XCH4 enhancements. This reflects that TCCON stations are not close to any of the corrected artefacts, hinting at a limitation of the current validation approach of the S5-P XCH4 data product.

The Tropospheric Monitoring Instrument (TROPOMI) on the Copernicus Sentinel 5 Precursor (S5P) satellite, launched in October 2017, provides a wealth of atmospheric composition data, including total columns of carbon monoxide (TCCO) at high horizontal resolution (5.5 km × 7 km). Near-real-time TROPOMI TCCO data have been monitored in the global data assimilation system of the Copernicus Atmosphere Monitoring Service (CAMS) since November 2018 to assess the quality of the data. The CAMS system already routinely assimilates TCCO data from the Measurement of Pollution in the Troposphere (MOPITT) instrument and the Infrared Atmospheric Sounding Interferometer (IASI) outside the polar regions. The assimilation of TROPOMI TCCO data in the CAMS system was tested for the period 6 July to 31 December 2021, i.e. after the TROPOMI algorithm update to version 02.02.00 in July 2021. By assimilating TROPOMI TCCO observations, the CAMS CO columns increase by on average 8 %, resulting in an improved fit to independent observations (IAGOS aircraft profiles and NDACC Fourier transform infrared (FTIR) tropospheric and total-column CO data) compared to a version of the CAMS system where only TCCO from MOPITT and IASI is assimilated. The largest absolute and relative changes from the assimilation of TROPOMI CO are found in the lower and middle troposphere, i.e. that part of the atmosphere that is not already well constrained by the assimilated TIR MOPITT and IASI data. The largest impact near the surface comes from clear-sky TROPOMI data over land, and additional vertical information comes from the retrievals of measurements in cloudy conditions. July and August 2021 saw record numbers of boreal wildfires over North America and Russia, leading to large amounts of CO being released into the atmosphere. The paper assesses the impact of TROPOMI CO assimilation on selected CO plumes more closely. While the CO column can be well constrained by the assimilation of TROPOMI CO data, and the fit to individual IAGOS CO profiles in the lower and middle troposphere is considerably improved, the TROPOMI CO columns do not provide further constraints on individual plumes that are transported across continents and oceans at altitudes above 500 hPa.

The TROPOspheric Monitoring Instrument (TROPOMI), due to its wide swath, performs observations over the ocean in every orbit, enhancing the monitoring capabilities of methane from space. In the short-wave–infrared (SWIR) spectral band ocean surfaces are dark except for the specific sun glint geometry, for which the specular reflectance detected by the satellite provides a signal that is high enough to retrieve methane with high accuracy and precision. In this study, we build upon the RemoTeC full-physics retrieval algorithm for land measurements, and we retrieve 4 years of methane concentrations over the ocean from TROPOMI. We fully assess the quality of the dataset by performing a validation using ground-based measurements of the Total Carbon Column Observing Network (TCCON) from near-ocean sites. The validation results in an agreement of -0.5±0.3 % (-8.4±6.3 ppb) for the mean bias and station-to-station variability, which show that glint measurements comply with the mission requirement of precision and accuracy below 1 %. Comparison to ocean measurements from the Greenhouse gases Observing SATellite (GOSAT) results in a bias of -0.2±0.9 % (-4.4±15.7 ppb), equivalent to the comparison of measurements over land. The full-physics algorithm simultaneously retrieves the amount of atmospheric methane and the physical scattering properties of the atmosphere from measurements in the near-infrared (NIR) and SWIR spectral bands. Based on the scattering properties of the atmosphere and ocean surface reflection we further validate retrievals over the ocean. Using the “upper-edge” method, we identify a set of ocean glint observations where scattering by aerosols and clouds can be ignored in the measurement simulation to investigate other possible error sources such as instrumental errors, radiometric inaccuracies or uncertainties related to spectroscopic absorption cross-sections. With this ensemble we evaluate the RemoTeC forward model via the validation of the total atmospheric oxygen (O2) column retrieved from the O2 A-band, as well as the consistency of XCH4 retrievals using sub-bands from the SWIR band, which show a consistency within 1 %. We discard any instrumental and radiometric errors by a calibration of the O2 absorption line strengths as suggested in the literature.

AirCore in situ vertical profiles sample the atmosphere from near the surface to the lower stratosphere, making them ideal for the validation of satellite tropospheric trace gas data. Here we present intercomparison results of AirCore carbon monoxide (CO) measurements with respect to retrievals from MOPITT (Measurements of Pollution In The Troposphere; version 8) and TROPOMI (TROPOspheric Monitoring Instrument), on board the NASA Terra and ESA Sentinel 5-Precursor satellites, respectively. Mean MOPITT/AirCore total column bias values and their standard deviation (0.4 ± 5.5, 1.7 ± 5.6, and 0.7 ± 6.0 for MOPITT thermal-infrared, near-infrared, and multispectral retrievals, respectively; all in %) are similar to results obtained in MOPITT/NOAA aircraft flask data comparisons from this study and from previous validation efforts. MOPITT CO retrievals are systematically validated using in situ vertical profiles from a variety of aircraft campaigns. Because most aircraft vertical profiles do not sample the troposphere's entire vertical extent, they must be extended upwards in order to be usable in validation. Here we quantify, for the first time, the error introduced in MOPITT CO validation by the use of shorter aircraft vertical profiles extended upwards by analyzing validation results of MOPITT with respect to full and truncated AirCore CO vertical profiles. Our results indicate that the error is small, affects mostly upper tropospheric retrievals (at 300 hPa: ∼ 2.6, 0.8, and 3.2 percent points for MOPITT thermal-infrared, near-infrared, and multispectral, respectively), and may have resulted in the overestimation of MOPITT retrieval biases in that region. TROPOMI can retrieve CO under both clear and cloudy conditions. The latter is achieved by quantifying interfering trace gases and parameters describing the cloud contamination of the measurements together with the CO column; then, the reference CO profiles used in the retrieval are scaled based on estimated above-cloud CO rather than on estimated total CO. We use AirCore measurements as the reference to evaluate the error introduced by this approach in cloudy TROPOMI retrievals over land after accounting for TROPOMI's vertical sensitivity to CO (relative bias and its standard deviation = 2.02 % ± 11.13 %). We also quantify the null-space error, which accounts for differences between the shape of TROPOMI reference profiles and that of AirCore measured profiles (for TROPOMI cloudy enull=0.98 % ± 2.32 %).

The aim of this paper is to highlight how TROPOspheric Monitoring Instrument (TROPOMI) trace gas data can best be used and interpreted to understand event-based impacts on air quality from regional to city scales around the globe. For this study, we present the observed changes in the atmospheric column amounts of five trace gases (NO2, SO2, CO, HCHO, and CHOCHO) detected by the Sentinel-5P TROPOMI instrument and driven by reductions in anthropogenic emissions due to COVID-19 lockdown measures in 2020. We report clear COVID-19-related decreases in TROPOMI NO2 column amounts on all continents. For megacities, reductions in column amounts of tropospheric NO2 range between 14 % and 63 %. For China and India, supported by NO2 observations, where the primary source of anthropogenic SO2 is coal-fired power generation, we were able to detect sector-specific emission changes using the SO2 data. For HCHO and CHOCHO, we consistently observe anthropogenic changes in 2-week-averaged column amounts over China and India during the early phases of the lockdown periods. That these variations over such a short timescale are detectable from space is due to the high resolution and improved sensitivity of the TROPOMI instrument. For CO, we observe a small reduction over China, which is in concert with the other trace gas reductions observed during lockdown; however, large interannual differences prevent firm conclusions from being drawn. The joint analysis of COVID-19-lockdown-driven reductions in satellite-observed trace gas column amounts using the latest operational and scientific retrieval techniques for five species concomitantly is unprecedented. However, the meteorologically and seasonally driven variability of the five trace gases does not allow for drawing fully quantitative conclusions on the reduction in anthropogenic emissions based on TROPOMI observations alone. We anticipate that in future the combined use of inverse modeling techniques with the high spatial resolution data from S5P/TROPOMI for all observed trace gases presented here will yield a significantly improved sector-specific, space-based analysis of the impact of COVID-19 lockdown measures as compared to other existing satellite observations. Such analyses will further enhance the scientific impact and societal relevance of the TROPOMI mission.

The thermal infrared nadir spectra of IASI (Infrared Atmospheric Sounding Interferometer) are successfully used for retrievals of different atmospheric trace gas profiles. However, these retrievals offer generally reduced information about the lowermost tropospheric layer due to the lack of thermal contrast close to the surface. Spectra of scattered solar radiation observed in the near-infrared and/or shortwave infrared, for instance by TROPOMI (TROPOspheric Monitoring Instrument), offer higher sensitivity near the ground and are used for the retrieval of total-column-averaged mixing ratios of a variety of atmospheric trace gases. Here we present a method for the synergetic use of IASI profile and TROPOMI total-column level 2 retrieval products. Our method uses the output of the individual retrievals and consists of linear algebra a posteriori calculations (i.e. calculation after the individual retrievals). We show that this approach has strong theoretical similarities to applying the spectra of the different sensors together in a single retrieval procedure but with the substantial advantage of being applicable to data generated with different individual retrieval processors, of being very time efficient, and of directly benefiting from the high quality and most recent improvements of the individual retrieval processors. We demonstrate the method exemplarily for atmospheric methane (CH4). We perform a theoretical evaluation and show that the a posteriori combination method yields a total-column-averaged CH4 product (XCH4) that conserves the good sensitivity of the corresponding TROPOMI product while merging it with the high-quality upper troposphere–lower stratosphere (UTLS) CH4 partial-column information of the corresponding IASI product. As a consequence, the combined product offers additional sensitivity for the tropospheric CH4 partial column, which is not provided by the individual TROPOMI nor the individual IASI product. The theoretically predicted synergetic effect is verified by comparisons to CH4 reference data obtained from collocated XCH4 measurements at 14 globally distributed TCCON (Total Carbon Column Observing Network) stations, CH4 profile measurements made by 36 individual AirCore soundings, and tropospheric CH4 data derived from continuous ground-based in situ observations made at two nearby Global Atmospheric Watch (GAW) mountain stations. The comparisons clearly demonstrate that the combined product can reliably detect the actual variations of atmospheric XCH4, CH4 in the UTLS, and CH4 in the troposphere. A similar good reliability for the latter is not achievable by the individual TROPOMI and IASI products.

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.

This paper presents an extended scientific HDO/H2O total column data product from short-wave infrared (SWIR) measurements by the Tropospheric Monitoring Instrument (TROPOMI) including clear-sky and cloudy scenes. The retrieval employs a forward model which accounts for scattering, and the algorithm infers the trace gas column information, surface properties, and effective cloud parameters from the observations. Compared to the previous clear-sky-only data product, coverage is greatly enhanced by including scenes over low clouds, particularly enabling data over oceans as the albedo of water in the SWIR spectral range is too low to retrieve under cloud-free conditions. The new dataset is validated against co-located ground-based Fourier transform infrared (FTIR) observations by the Total Carbon Column Observing Network (TCCON). The median bias for clear-sky scenes is 1.4×1021 molec cm−2 (2.9 %) in H2O columns and 1.1×1017 molec cm−2 (−0.3 %) in HDO columns, which corresponds to −17 ‰ (9.9 %) in a posteriori δD. The bias for cloudy scenes is 4.9×1021 molec cm−2 (11 %) in H2O, 1.1×1018 molec cm−2 (7.9 %) in HDO, and −20 ‰ (9.7 %) in a posteriori δD. At low-altitude stations, the bias is small at low and middle latitudes and has a larger value at high latitudes. At high-altitude stations, an altitude correction is required to compensate for different partial columns seen by the station and the satellite. The bias in a posteriori δD after altitude correction depends on sensitivity due to shielding by clouds and on realistic a priori profile shapes for both isotopologues. Cloudy scenes generally involve low sensitivity below the clouds, and since the information is filled up by the prior, a realistic shape of the prior is important for realistic total column estimation in these cases. Over oceans, aircraft measurements with the Water Isotope System for Precipitation and Entrainment Research (WISPER) instrument from a field campaign in 2018 are used for validation, yielding biases of −3.9 % in H2O and −3 ‰ in δD over clouds. To demonstrate the added value of the new dataset, a short case study of a cold air outbreak over the Atlantic Ocean in January 2020 is presented, showing the daily evolution of the event with single-overpass results.

The objective of this study is to derive methane (CH4) emissions from three landfills, which are found to be the most significant CH4 sources in the metropolitan area of Madrid in Spain. We derive CH4 emissions from the CH4 enhancements observed by spaceborne and ground-based instruments. We apply satellite-based measurements from the TROPOspheric Monitoring Instrument (TROPOMI) and the Infrared Atmospheric Sounding Interferometer (IASI) together with measurements from the ground-based COllaborative Carbon Column Observing Network (COCCON) instruments. In 2018, a 2-week field campaign for measuring the atmospheric concentrations of greenhouse gases was performed in Madrid in the framework of Monitoring of the Greenhouse Gases Concentrations in Madrid (MEGEI-MAD) project. Five COCCON instruments were deployed at different locations around the Madrid city center, enabling the observation of total column-averaged CH4 mixing ratios (XCH4). Considering the prevalent wind regimes, we calculate the wind-assigned XCH4 anomalies for two opposite wind directions. Pronounced bipolar plumes are found when applying the method to NO2, which implies that our method of wind-assigned anomaly is suitable to estimate enhancements of trace gases at the urban level from satellite-based measurements. For quantifying the CH4 emissions, the wind-assigned plume method is applied to the TROPOMI XCH4 and to the lower tropospheric CH4 / dry-air column ratio (TXCH4) of the combined TROPOMI+IASI product. As CH4 emission strength we estimate 7.4 × 1025 ± 6.4 × 1024 molec. s−1 from the TROPOMI XCH4 data and 7.1 × 1025 ± 1.0 × 1025 molec. s−1 from the TROPOMI+IASI merged TXCH4 data. We use COCCON observations to estimate the local source strength as an independent method. COCCON observations indicate a weaker CH4 emission strength of 3.7 × 1025 molec. s−1 from a local source (the Valdemingómez waste plant) based on observations from a single day. This strength is lower than the one derived from the satellite observations, and it is a plausible result. This is because the analysis of the satellite data refers to a larger area, covering further emission sources in the study region, whereas the signal observed by COCCON is generated by a nearby local source. All emission rates estimated from the different observations are significantly larger than the emission rates provided via the official Spanish Register of Emissions and Pollutant Sources.

Water vapour is the most abundant natural greenhouse gas in the Earth's atmosphere, and global data sets are required for meteorological applications and climate research. The Tropospheric Monitoring Instrument (TROPOMI) on board Sentinel-5 Precursor (S5P) launched on 13 October 2017 has a high spatial resolution of around 5 km and a daily global coverage. Currently, there is no operational total water vapour product for S5P measurements. Here, we present first results of a new scientific total column water vapour (TCWV) product for S5P using the so-called air-mass-corrected differential optical absorption spectroscopy (AMC-DOAS) scheme. This method analyses spectral data between 688 and 700 nm and has already been successfully applied to measurements from the Global Monitoring Experiment (GOME) on ERS-2, the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on Envisat and GOME-2 on MetOp. The adaptation of the AMC-DOAS method to S5P data requires an additional post-processing procedure to correct the influences of surface albedo, cloud height and cloud fraction. The quality of the new AMC-DOAS S5P water vapour product is assessed by comparisons with data from GOME-2 on MetOp-B retrieved also with the AMC-DOAS algorithm and with four independent data sets, namely reanalysis data from the European Centre for Medium range Weather Forecast (ECMWF ERA5), data obtained by the Special Sensor Microwave Imager and Sounder (SSMIS) flown on the Defense Meteorological Satellite Program (DMSP) platform 16 and two scientific S5P TCWV products derived from TROPOMI measurements. Both are recently published TCWV products for S5P provided by the Max Planck Institute for Chemistry (MPIC) in Mainz and the Netherlands Institute for Space Research (SRON), Utrecht. The SRON TCWV is limited to clear-sky scenes over land. These comparisons reveal a good agreement between the various data sets but also some systematic differences between all of them. On average, the daily derived offset between AMC-DOAS S5P TCWV and AMC-DOAS GOME-2B TCWV is negative (around −1.5 kg m−2) over land and positive over ocean surfaces (more than 1.5 kg m−2). In contrast, SSMIS TCWV is on average lower than AMC-DOAS S5P TCWV by about 3 kg m−2. Monthly averaged ERA5 TCWV and AMC-DOAS S5P TCWV comparison shows spatial features over both land and water surface. Over land, there are systematic spatial structures. There are larger differences between AMC-DOAS S5P TCWV and ERA5 TCWV in tropical regions. Over sea, AMC-DOAS S5P TCWV is slightly lower than ERA5 TCWV by around 2 kg m−2. The AMC-DOAS S5P TCWV and S5P TCWV from MPIC agree on average within 1 kg m−2 over both land and ocean. TCWV from SRON shows daily global averaged differences to AMC-DOAS S5P TCWV of around 1.2 kg m−2. All of these differences are in line with the accuracy of these products and with the typical range of differences of 5 kg m−2 obtained when comparing different TCWV data sets. The AMC-DOAS TCWV product for S5P provides therefore a valuable new and independent data set for atmospheric applications which also has a higher spatial coverage than the other S5P TCWV products.

The TROPOspheric Monitoring Instrument (TROPOMI) provides a daily, spatially resolved (initially 7×7 km2, upgraded to 7×5.6 km2 in August 2019) global dataset of CO columns; however, due to the relative sparseness of reliable ground-based data sources, it can be challenging to characterize the validity and accuracy of satellite data products in remote regions such as the high Arctic. In these regions, satellite intercomparisons can supplement model- and ground-based validation efforts and serve to verify previously observed differences. In this paper, we compare the CO products from TROPOMI, the Atmospheric Chemistry Experiment (ACE) Fourier transform spectrometer (FTS), and a high-Arctic ground-based FTS located at the Polar Environment Atmospheric Research Laboratory (PEARL) in Eureka, Nunavut (80.05∘ N, 86.42∘ W). A global comparison of TROPOMI reference profiles scaled by the retrieved total column with ACE-FTS CO partial columns for the period from 28 November 2017 to 31 May 2020 displays excellent agreement between the two datasets (R=0.93) and a small relative bias of -0.83±0.26% (bias ± standard error of the mean). Additional comparisons were performed within five latitude bands: the north polar region (60 to 90∘ N), northern mid-latitudes (20 to 60∘ N), the equatorial region (20∘ S to 20∘ N), southern mid-latitudes (60 to 20∘ S), and the south polar region (90 to 60∘ S). Latitudinal comparisons of the TROPOMI and ACE-FTS CO datasets show strong correlations ranging from R=0.93 (southern mid-latitudes) to R=0.86 (equatorial region) between the CO products but display a dependence of the mean differences on latitude. Positive mean biases of 7.93±0.61 % and 7.21±0.52 % were found in the northern and southern polar regions, respectively, while a negative bias of -9.41±0.55% was observed in the equatorial region. To investigate whether these differences are introduced by cloud contamination, which is reflected in the TROPOMI averaging kernel shape, the latitudinal comparisons were repeated for cloud-covered pixels and clear-sky pixels only, as well as for the unsmoothed and smoothed cases. Clear-sky pixels were found to be biased higher with poorer correlations on average than clear+cloudy scenes and cloud-covered scenes only. Furthermore, the latitudinal dependence on the biases was observed in both the smoothed and unsmoothed cases. To provide additional context to the global comparisons of TROPOMI with ACE-FTS in the Arctic, both satellite datasets were compared against measurements from the ground-based PEARL-FTS. Comparisons of TROPOMI with smoothed PEARL-FTS total columns in the period of 3 March 2018 to 27 March 2020 display a strong correlation (R=0.88); however, a positive mean bias of 14.7±0.16 % was also found. A partial column comparison of ACE-FTS with the PEARL-FTS in the period from 25 February 2007 to 18 March 2020 shows good agreement (R=0.79) and a mean positive bias of 7.89±0.21 % in the ACE-FTS product relative to the ground-based FTS. The magnitude and sign of the mean relative differences are consistent across all intercomparisons in this work, as well as with recent ground-based validation efforts, suggesting that the current TROPOMI CO product exhibits a positive bias in the high-Arctic region. However, the observed bias is within the TROPOMI mission accuracy requirement of ±15 %, providing further confirmation that the data quality in these remote high-latitude regions meets this specification.

The Sentinel-5 Precursor (S5P) mission with the TROPOspheric Monitoring Instrument (TROPOMI) on board has been measuring solar radiation backscattered by the Earth's atmosphere and surface since its launch on 13 October 2017. In this paper, we present for the first time the S5P operational methane (CH4) and carbon monoxide (CO) products' validation results covering a period of about 3 years using global Total Carbon Column Observing Network (TCCON) and Infrared Working Group of the Network for the Detection of Atmospheric Composition Change (NDACC-IRWG) network data, accounting for a priori alignment and smoothing uncertainties in the validation, and testing the sensitivity of validation results towards the application of advanced co-location criteria. We found that the S5P standard and bias-corrected CH4 data over land surface for the recommended quality filtering fulfil the mission requirements. The systematic difference of the bias-corrected total column-averaged dry air mole fraction of methane (XCH4) data with respect to TCCON data is -0.26±0.56 % in comparison to -0.68±0.74 % for the standard XCH4 data, with a correlation of 0.6 for most stations. The bias shows a seasonal dependence. We found that the S5P CO data over all surfaces for the recommended quality filtering generally fulfil the missions requirements, with a few exceptions, which are mostly due to co-location mismatches and limited availability of data. The systematic difference between the S5P total column-averaged dry air mole fraction of carbon monoxide (XCO) and the TCCON data is on average 9.22±3.45 % (standard TCCON XCO) and 2.45±3.38 % (unscaled TCCON XCO). We found that the systematic difference between the S5P CO column and NDACC CO column (excluding two outlier stations) is on average 6.5±3.54 %. We found a correlation of above 0.9 for most TCCON and NDACC stations. The study shows the high quality of S5P CH4 and CO data by validating the products against reference global TCCON and NDACC stations covering a wide range of latitudinal bands, atmospheric conditions and surface conditions.

In this paper, we compare column-averaged dry-air mole fractions of water vapor (XH2O) retrievals from the COllaborative Carbon Column Observing Network (COCCON) with retrievals from two co-located high-resolution Fourier transform infrared (FTIR) spectrometers as references at two boreal sites, Kiruna, Sweden, and Sodankylä, Finland, from 6 March 2017 to 20 September 2019. In the framework of the Network for the Detection of Atmospheric Composition Change (NDACC), an FTIR spectrometer is operated at Kiruna. The H2O product derived from these observations has been generated with the MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) processor. In Sodankylä, a Total Carbon Column Observing Network (TCCON) spectrometer is operated, and the official XH2O data as provided by TCCON are used for this study. The datasets are in good overall agreement, with COCCON data showing a wet bias of (49.20±58.61) ppm ((3.33±3.37) %, R2=0.9992) compared with MUSICA NDACC and (56.32±45.63) ppm ((3.44±1.77) %, R2=0.9997) compared with TCCON. Furthermore, the a priori H2O volume mixing ratio (VMR) profiles (MAP) used as a priori information in the TCCON retrievals (also adopted for COCCON retrievals) are evaluated with respect to radiosonde (Vaisala RS41) profiles at Sodankylä. The MAP and radiosonde profiles show similar shapes and a good linear correlation of integrated XH2O, indicating that MAP is a reasonable approximation of the true atmospheric state and an appropriate choice for the scaling retrieval methods as applied by COCCON and TCCON. COCCON shows a reduced dry bias (−14.96 %) in comparison with TCCON (−19.08 %) with respect to radiosonde XH2O. Finally, we investigate the quality of satellite data at high latitudes. For this purpose, the COCCON XH2O is compared with retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) generated with the MUSICA processor (MUSICA IASI) and with retrievals from the TROPOspheric Monitoring Instrument (TROPOMI). Both paired datasets generally show good agreement and similar correlations at the two sites. COCCON measures 4.64 % less XH2O at Kiruna and 3.36 % less at Sodankylä with respect to MUSICA IASI, whereas COCCON measures 9.71 % more XH2O at Kiruna and 7.75 % more at Sodankylä compared with TROPOMI. Our study supports the assumption that COCCON also delivers a well-characterized XH2O data product. This emphasizes that this approach might complement the TCCON network with respect to satellite validation efforts. This is the first published study where COCCON XH2O has been compared with MUSICA NDACC and TCCON retrievals and has been used for MUSICA IASI and TROPOMI validation.

The global fire emission inventories depend on ground and airborne measurements of species-specific emission factors (EFs), which translate dry matter losses due to fires to actual trace gas and aerosol emissions. The EFs of nitrogen oxides (NOx) and carbon monoxide (CO) can function as a proxy for combustion efficiency to distinguish flaming from smoldering combustion. The uncertainties in these EFs remain large as they are limited by the spatial and temporal representativeness of the measurements. The global coverage of satellite observations has the advantage of filling this gap, making these measurements highly complementary to ground-based or airborne data. We present a new analysis of biomass burning pollutants using space-borne data to investigate the spatiotemporal efficiency of fire combustion. Column measurements of nitrogen dioxide and carbon monoxide (XNO2 and XCO) from the TROPOspheric Monitoring Instrument (TROPOMI) are used to quantify the relative atmospheric enhancements of these species over different fire-prone regions around the world. We find spatial and temporal patterns in the ΔXNO2 ∕ ΔXCO ratio that point to distinct differences in biomass burning behavior. Such differences are induced by the burning phase of the fire (e.g., high-temperature flaming vs. low-temperature smoldering combustion) and burning practice (e.g., the combustion of logs, coarse woody debris and soil organic matter vs. the combustion of fine fuels such as savanna grasses). The sampling techniques and the signal-to-noise ratio of the retrieved ΔXNO2 ∕ ΔXCO signals were quantified with WRF-Chem experiments and showed similar distinct differences in combustion types. The TROPOMI measurements show that the fraction of surface smoldering combustion is much larger for the boreal forest fires in the upper Northern Hemisphere and peatland fires in Indonesia. These types of fires cause a much larger increase (3 to 6 times) in ΔXCO relative to ΔXNO2 than elsewhere in the world. The high spatial and temporal resolution of TROPOMI also enables the detection of spatial gradients in combustion efficiency at smaller regional scales. For instance, in the Amazon, we found higher combustion efficiency (up to 3-fold) for savanna fires than for the nearby tropical deforestation fires. Out of two investigated fire emission products, the TROPOMI measurements support the broad spatial pattern of combustion efficiency rooted in GFED4s. Meanwhile, TROPOMI data also add new insights into regional variability in combustion characteristics that are not well represented in the different emission inventories, which can help the fire modeling community to improve their representation of the spatiotemporal variability in EFs.

The TROPOspheric Monitoring Instrument (TROPOMI) provides observations of atmospheric methane (CH4) at an unprecedented combination of high spatial resolution and daily global coverage. Hu et al. (2018) reported unexpectedly large methane enhancements over South Sudan in these observations. Here we assess methane emissions from the wetlands of South Sudan using 2 years (December 2017–November 2019) of TROPOMI total column methane observations. We estimate annual wetland emissions of 7.4 ± 3.2 Tg yr−1, which agrees with the multiyear GOSAT inversions of Lunt et al. (2019) but is an order of magnitude larger than estimates from wetland process models. This disagreement may be explained by the underestimation (by up to 4 times) of inundation extent by the hydrological schemes used in those models. We investigate the seasonal cycle of the emissions and find the lowest emissions during the June–August season when the process models show the largest emissions. Using satellite-altimetry-based river water height measurements, we infer that this seasonal mismatch is likely due to a seasonal mismatch in inundation extent. In models, inundation extent is controlled by regional precipitation scaled to static wetland extent maps, whereas the actual inundation extent is driven by water inflow from rivers like the White Nile and the Sobat. We find the lowest emissions in the highest precipitation and lowest temperature season (June–August, JJA) when models estimate large emissions. In general, our emission estimates show better agreement in terms of both seasonal cycle and annual mean with model estimates that use a stronger temperature dependence. This suggests that temperature might be a stronger control for the South Sudan wetlands emissions than currently assumed by models. Our findings demonstrate the use of satellite instruments for quantifying emissions from inaccessible and uncertain tropical wetlands, providing clues for the improvement of process models and thereby improving our understanding of the currently uncertain contribution of wetlands to the global methane budget.

The TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel 5 Precursor (S5-P) satellite provides methane (CH4) measurements with high accuracy and exceptional temporal and spatial resolution and sampling. TROPOMI CH4 measurements are highly valuable to constrain emissions inventories and for trend analysis, with strict requirements on the data quality. This study describes the improvements that we have implemented to retrieve CH4 from TROPOMI using the RemoTeC full-physics algorithm. The updated retrieval algorithm features a constant regularization scheme of the inversion that stabilizes the retrieval and yields less scatter in the data and includes a higher resolution surface altitude database. We have tested the impact of three state-of-the-art molecular spectroscopic databases (HITRAN 2008, HITRAN 2016 and Scientific Exploitation of Operational Missions – Improved Atmospheric Spectroscopy Databases SEOM-IAS) and found that SEOM-IAS provides the best fitting results. The most relevant update in the TROPOMI XCH4 data product is the implementation of an a posteriori correction fully independent of any reference data that is more accurate and corrects for the underestimation at low surface albedo scenes and the overestimation at high surface albedo scenes. After applying the correction, the albedo dependence is removed to a large extent in the TROPOMI versus satellite (Greenhouse gases Observing SATellite – GOSAT) and TROPOMI versus ground-based observations (Total Carbon Column Observing Network – TCCON) comparison, which is an independent verification of the correction scheme. We validate 2 years of TROPOMI CH4 data that show the good agreement of the updated TROPOMI CH4 with TCCON (−3.4 ± 5.6 ppb) and GOSAT (−10.3 ± 16.8 ppb) (mean bias and standard deviation). Low- and high-albedo scenes as well as snow-covered scenes are the most challenging for the CH4 retrieval algorithm, and although the a posteriori correction accounts for most of the bias, there is a need to further investigate the underlying cause.

Measurements of atmospheric column-averaged dry-air mole fractions of carbon dioxide (XCO2), methane (XCH4), and carbon monoxide (XCO) have been collected across the Pacific Ocean during the Measuring Ocean REferences 2 (MORE-2) campaign in June 2019. We deployed a shipborne variant of the EM27/SUN Fourier transform spectrometer (FTS) on board the German R/V Sonne which, during MORE-2, crossed the Pacific Ocean from Vancouver, Canada, to Singapore. Equipped with a specially manufactured fast solar tracker, the FTS operated in direct-sun viewing geometry during the ship cruise reliably delivering solar absorption spectra in the shortwave infrared spectral range (4000 to 11000 cm−1). After filtering and bias correcting the dataset, we report on XCO2, XCH4, and XCO measurements for 22 d along a trajectory that largely aligns with 30∘ N of latitude between 140∘ W and 120∘ E of longitude. The dataset has been scaled to the Total Carbon Column Observing Network (TCCON) station in Karlsruhe, Germany, before and after the MORE-2 campaign through side-by-side measurements. The 1σ repeatability of hourly means of XCO2, XCH4, and XCO is found to be 0.24 ppm, 1.1 ppb, and 0.75 ppb, respectively. The Copernicus Atmosphere Monitoring Service (CAMS) models gridded concentration fields of the atmospheric composition using assimilated satellite observations, which show excellent agreement of 0.52±0.31 ppm for XCO2, 0.9±4.1 ppb for XCH4, and 3.2±3.4 ppb for XCO (mean difference ± SD, standard deviation, of differences for entire record) with our observations. Likewise, we find excellent agreement to within 2.2±6.6 ppb with the XCO observations of the TROPOspheric MOnitoring Instrument (TROPOMI) on the Sentinel-5 Precursor satellite (S5P). The shipborne measurements are accessible at https://doi.org/10.1594/PANGAEA.917240 (Knapp et al., 2020).

The Tropospheric Monitoring Instrument (TROPOMI) on the ESA Copernicus Sentinel-5 satellite (S5-P) measures carbon monoxide (CO) total column concentrations as one of its primary targets. In this study, we analyze TROPOMI observations over Mexico City in the period 14 November 2017 to 25 August 2019 by means of collocated CO simulations using the regional Weather Research and Forecasting coupled with Chemistry (WRF-Chem) model. We draw conclusions on the emissions from different urban districts in the region. Our WRF-Chem simulation distinguishes CO emissions from the districts Tula, Pachuca, Tulancingo, Toluca, Cuernavaca, Cuautla, Tlaxcala, Puebla, Mexico City, and Mexico City Arena by 10 separate tracers. For the data interpretation, we apply a source inversion approach determining per district the mean emissions and the temporal variability, the latter regularized to reduce the propagation of the instrument noise and forward-model errors in the inversion. In this way, the TROPOMI observations are used to evaluate the Inventario Nacional de Emisiones de Contaminantes Criterio (INEM) inventory that was adapted to the period 2017–2019 using in situ ground-based observations. For the Tula and Pachuca urban areas in the north of Mexico City, we obtain 0.10±0.004 and 0.09±0.005 Tg yr−1 CO emissions, which exceeds significantly the INEM emissions of Tg yr−1 for both areas. On the other hand for Mexico City, TROPOMI estimates emissions of 0.14±0.006 Tg yr−1 CO, which is about half of the INEM emissions of 0.25 Tg yr−1, and for the adjacent district Mexico City Arena the emissions are 0.28±0.01 Tg yr−1 according to TROPOMI observations versus 0.14 Tg yr−1 as stated by the INEM inventory. Interestingly, the total emissions of both districts are similar (0.42±0.016 Tg yr−1 TROPOMI versus 0.39 Tg yr−1 adapted INEM emissions). Moreover, for both areas we found that the TROPOMI emission estimates follow a clear weekly cycle with a minimum during the weekend. This agrees well with ground-based in situ measurements from the Secretaría del Medio Ambiente (SEDEMA) and Fourier transform spectrometer column measurements in Mexico City that are operated by the Network for the Detection of Atmospheric Composition Change Infrared Working Group (NDACC-IRWG). Overall, our study demonstrates an approach to deploying the large number of TROPOMI CO data to draw conclusions on urban emissions on sub-city scales for metropolises like Mexico City. Moreover, for the exploitation of TROPOMI CO observations our analysis indicates the clear need for further improvements of regional models like WRF-Chem, in particular with respect to the prediction of the local wind fields.

We have analyzed TROPOspheric Monitoring Instrument (TROPOMI) carbon monoxide (CO) data acquired between November 2017 and March 2019 with respect to other satellite (MOPITT, Measurement Of Pollution In The Troposphere) and airborne (ATom, Atmospheric Tomography mission) datasets to better understand TROPOMI's contribution to the global tropospheric CO record (2000 to present). MOPITT and TROPOMI are two of only a few satellite instruments to ever derive CO from solar-reflected radiances. Therefore, it is particularly important to understand how these two datasets compare. Our results indicate that TROPOMI CO retrievals over land show excellent agreement with respect to MOPITT: relative biases and their SD (i.e., accuracy and precision) are on average -3.73%±11.51%, -2.24%±12.38%, and -3.22%±11.13% compared to the MOPITT TIR (thermal infrared), NIR (near infrared), and TIR + NIR (multispectral) products, respectively. TROPOMI and MOPITT data also show good agreement in terms of temporal and spatial patterns. Despite depending on solar-reflected radiances for its measurements, TROPOMI can also retrieve CO over bodies of water if clouds are present by approximating partial columns under cloud tops using scaled, model-based reference CO profiles. We quantify the bias of TROPOMI total column retrievals over bodies of water with respect to colocated in situ ATom CO profiles after smoothing the latter with the TROPOMI column averaging kernels (AKs), which account for signal attenuation under clouds (relative bias and its SD =3.25%±11.46 %). In addition, we quantify enull (the null-space error), which accounts for differences between the shape of the TROPOMI reference profile and that of the ATom true profile (enull=2.16%±2.23 %). For comparisons of TROPOMI and MOPITT retrievals over open water we compare TROPOMI total CO columns to their colocated MOPITT TIR counterparts. Relative bias and its SD are 2.98 %±15.71 % on average. We investigate the impact of discrepancies between the a priori and reference CO profiles (used by MOPITT and TROPOMI, respectively) on CO retrieval biases by applying a null-space adjustment (based on the MOPITT a priori) to the TROPOMI total column values. The effect of this adjustment on MOPITT and TROPOMI biases is minor, typically 1–2 percentage points.

This study investigates the use of co-located nitrogen dioxide (NO2) and carbon monoxide (CO) retrievals from the TROPOMI satellite to improve the quantification of burning efficiency and emission factors (EFs) over the megacities of Tehran, Mexico City, Cairo, Riyadh, Lahore, and Los Angeles. Efficient combustion is characterized by high NOx (NO+NO2) and low CO emissions, making the NO2∕CO ratio a useful proxy for combustion efficiency (CE). The local enhancement of CO and NO2 above megacities is well captured by TROPOMI at short averaging times compared with previous satellite missions. In this study, the upwind background and plume rotation methods are used to investigate the accuracy of satellite-derived ΔNO2∕ΔCO ratios. The column enhancement ratios derived using these two methods vary by 5 % to 20 % across the selected megacities. TROPOMI-derived column enhancement ratios are compared with emission ratios from the EDGAR v4.3.2 (Emission Database for Global Atmospheric Research v4.3.2) and the MACCity (Monitoring Atmospheric Chemistry and Climate and CityZen) 2018 emission inventories. TROPOMI correlates strongly (r=0.85 and 0.7) with EDGAR and MACCity, showing the highest emission ratio for Riyadh and lowest emission ratio for Lahore. However, inventory-derived emission ratios are 60 % to 85 % higher than TROPOMI column enhancement ratios across the six megacities. The short lifetime of NO2 and the different vertical sensitivity of TROPOMI NO2 and CO explain most of this difference. We present a method to translate TROPOMI-retrieved column enhancement ratios into corresponding emission ratios, thereby accounting for these influences. Except for Los Angeles and Lahore, TROPOMI-derived emission ratios are close (within 10 % to 25 %) to MACCity values. For EDGAR, however, emission ratios are ∼65 % higher for Cairo and 35 % higher for Riyadh. For Los Angeles, EDGAR and MACCity are a factor of 2 and 3 higher than TROPOMI respectively. The air quality monitoring networks in Los Angeles and Mexico City are used to validate the use of TROPOMI. For Mexico City and Los Angeles, these measurements are consistent with TROPOMI-derived emission ratios, demonstrating the potential of TROPOMI with respect to monitoring burning efficiency.

We compare the atmospheric column-averaged dry-air mole fractions of carbon dioxide (XCO2) and methane (XCH4) measured with a pair of COllaborative Carbon Column Observing Network (COCCON) spectrometers at Kiruna and Sodankylä (boreal areas). We compare model data provided by the Copernicus Atmosphere Monitoring Service (CAMS) between 2017 and 2019 with XCH4 data from the recently launched Sentinel-5 Precursor (S5P) satellite between 2018 and 2019. In addition, measured and modeled gradients of XCO2 and XCH4 (ΔXCO2 and ΔXCH4) on regional scales are investigated. Both sites show a similar and very good correlation between COCCON retrievals and the modeled CAMS XCO2 data, while CAMS data are biased high with respect to COCCON by 3.72 ppm (±1.80 ppm) in Kiruna and 3.46 ppm (±1.73 ppm) in Sodankylä on average. For XCH4, CAMS values are higher than the COCCON observations by 0.33 ppb (±11.93 ppb) in Kiruna and 7.39 ppb (±10.92 ppb) in Sodankylä. In contrast, the S5P satellite generally measures lower atmospheric XCH4 than the COCCON spectrometers, with a mean difference of 9.69 ppb (±20.51 ppb) in Kiruna and 3.36 ppb (±17.05 ppb) in Sodankylä. We compare the gradients of XCO2 and XCH4 (ΔXCO2 and ΔXCH4) between Kiruna and Sodankylä derived from CAMS analysis and COCCON and S5P measurements to study the capability of detecting sources and sinks on regional scales. The correlations in ΔXCO2 and ΔXCH4 between the different datasets are generally smaller than the correlations in XCO2 and XCH4 between the datasets at either site. The ΔXCO2 values predicted by CAMS are generally higher than those observed with COCCON with a slope of 0.51. The ΔXCH4 values predicted by CAMS are mostly higher than those observed with COCCON with a slope of 0.65, covering a larger dataset than the comparison between S5P and COCCON. When comparing CAMS ΔXCH4 with COCCON ΔXCH4 only in S5P overpass days (slope = 0.53), the correlation is close to that between S5P and COCCON (slope = 0.51). CAMS, COCCON, and S5P predict gradients in reasonable agreement. However, the small number of observations coinciding with S5P limits our ability to verify the performance of this spaceborne sensor. We detect no significant impact of ground albedo and viewing zenith angle on the S5P results. Both sites show similar situations with the average ratios of XCH4 (S5P/COCCON) of 0.9949±0.0118 in Kiruna and 0.9953±0.0089 in Sodankylä. Overall, the results indicate that the COCCON instruments have the capability of measuring greenhouse gas (GHG) gradients on regional scales, and observations performed with the portable spectrometers can contribute to inferring sources and sinks and to validating spaceborne greenhouse gas sensors. To our knowledge, this is the first published study using COCCON spectrometers for the validation of XCH4 measurements collected by S5P.

Global measurements of atmospheric water vapour isotopologues aid to better understand the hydrological cycle and improve global circulation models. This paper presents a new data set of vertical column densities of H2O and HDO retrieved from short-wave infrared (2.3 µm) reflectance measurements by the Tropospheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite. TROPOMI features daily global coverage with a spatial resolution of up to 7 km×7 km. The retrieval utilises a profile-scaling approach. The forward model neglects scattering, and strict cloud filtering is therefore necessary. For validation, recent ground-based water vapour isotopologue measurements by the Total Carbon Column Observing Network (TCCON) are employed. A comparison of TCCON δD with ground-based measurements by the Multi-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) project for data prior to 2014 (where MUSICA data are available) shows a bias in TCCON δD estimates. As TCCON HDO is currently not validated, an overall correction of recent TCCON HDO data is derived based on this finding. The agreement between the corrected TCCON measurements and co-located TROPOMI observations is good with an average bias of (-0.2±3)×1021 molec cm−2 ((1.1±7.2) %) in H2O and (-2±7)×1017 molec cm−2 ((-1.1±7.3) %) in HDO, which corresponds to a mean bias of (-14±17) ‰ in a posteriori δD. The bias is lower at low- and mid-latitude stations and higher at high-latitude stations. The use of the data set is demonstrated with a case study of a blocking anticyclone in northwestern Europe in July 2018 using single-overpass data.

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.

On 13 October 2017, the Tropospheric Monitoring Instrument (TROPOMI) was launched on the Copernicus Sentinel-5 Precursor satellite in a sun-synchronous orbit. One of the mission's operational data products is the total column concentration of carbon monoxide (CO), which was released to the public in July 2018. The current TROPOMI CO processing uses the HITRAN 2008 spectroscopic data with updated water vapor spectroscopy and produces a CO data product compliant with the mission requirement of 10 % precision and 15 % accuracy for single soundings. Comparison with ground-based CO observations of the Total Carbon Column Observing Network (TCCON) show systematic differences of about 6.2 ppb and single-orbit observations are superimposed by a significant striping pattern along the flight path exceeding 5 ppb. In this study, we discuss possible improvements of the CO data product. We found that the molecular spectroscopic data used in the retrieval plays a key role for the data quality where the use of the Scientific Exploitation of Operational Missions – Improved Atmospheric Spectroscopy Databases (SEOM-IAS) and the HITRAN 2012 and 2016 releases reduce the bias between TROPOMI and TCCON due to improved CH4 spectroscopy. SEOM-IAS achieves the best spectral fit quality (root-mean-square, rms, differences between the simulated and measured spectrum) of 1.5×10-10 mol s−1 m−2 nm−1 sr−1 and reduces the bias between TROPOMI and TCCON to 3.4 ppb, while HITRAN 2012 and HITRAN 2016 decrease the bias even further below 1 ppb. HITRAN 2012 shows the worst fit quality (rms = 2.5×10-10 mol s−1 m−2 nm−1 sr−1) of the tested cross sections and furthermore introduces an artificial bias of about -1.5×1017 molec cm−2 between TROPOMI CO and the CAMS-IFS model in the Tropics caused by the H2O spectroscopic data. Moreover, analyzing 1 year of TROPOMI CO observations, we identified increased striping patterns by about 16 % percent from November 2017 to November 2018. For that, we defined a measure γ, quantifying the relative pixel-to-pixel variation in CO in the cross-track and along-track directions. To mitigate this effect, we discuss two destriping methods applied to the CO data a posteriori. A destriping mask calculated per orbit by median filtering of the data in the cross-track direction significantly reduced the stripe pattern from γ=2.1 to γ=1.6. However, the destriping can be further improved, achieving γ=1.2 by deploying a Fourier analysis and filtering of the data, which not only corrects for stripe patterns in the cross-track direction but also accounts for the variability of stripes along the flight path.

The TROPOspheric Monitoring Instrument (TROPOMI), launched 13 October 2017, has been measuring carbon monoxide (CO) concentrations in the Earth's atmosphere since early November 2017. In the first measurements, TROPOMI was able to measure CO concentrations of the high-pollution event in India of November 2017. In this paper, we studied the extent of the pollution in India, comparing the TROPOMI CO with modeled data from the Weather Research and Forecasting model (WRF) to identify the most important sources contributing to the high pollution, both at ground level and in the total column. We investigated the period 11–19 November 2017. We found that residential and commercial combustion was a much more important source of CO pollution than the post-monsoon crop burning during this period, which is in contrast to what media suggested and some studies on aerosol emissions found. Also, the high pollution was not limited to Delhi and its direct neighborhood but the accumulation of pollution extended over the whole Indo-Gangetic Plain (IGP) due to the unfavorable weather conditions in combination with extensive emissions. From the TROPOMI data and WRF simulations, we observed a buildup of CO during 11–14 November and a decline in CO after 15 November. The meteorological conditions, characterized by low wind speeds and shallow atmospheric boundary layers, were most likely the primary explanation for the temporal accumulation and subsequent dispersion of regionally emitted CO in the atmosphere. This emphasizes the important role of atmospheric dynamics in determining the air quality conditions at ground level and in the total column. Due to its rapidly growing population and economy, India is expected to encounter similar pollution events more often in future post-monsoon and winter seasons unless significant policy measures are taken to reduce residential and commercial emissions.

The Tropospheric Monitoring Instrument (TROPOMI) on the Sentinel-5 Precursor satellite provides measurements of carbon monoxide (CO) total column concentrations based on earthshine radiance measurements in the 2.3 µm spectral range with a spatial resolution of 7 km×7 km and daily global coverage. Due to the high accuracy of the observations, CO pollution can be detected over cities and industrial areas using single orbit overpasses. In this study, we analyzed local CO enhancements in an area around Iran from 1 November to 20 December 2017. We employed the Weather Research and Forecasting (WRF) model v3.8.1 using the EDGAR v4.2 emission inventory and evaluated CO emissions from the cities of Tehran, Yerevan, Urmia, and Tabriz on a spatial resolution comparable to that of TROPOMI. For background conditions, the WRF simulation agrees well with TROPOMI CO, with a mean difference of 5.7 %. However, the emissions for the city area had to be significantly increased in order to match the observations. Moreover, significant differences at the sub-city scale remain. To match the TROPOMI CO observations around the Armenian city of Yerevan, it is necessary to introduce CO emissions along a southeast arterial road of Yerevan. Overall, this hints at deficits in the EDGAR inventory in the region around Iran and indicates TROPOMI's capability to identify localized CO pollution on sub-city scales, which at the same time challenges current atmospheric modeling at high spatial and temporal resolution.

On 13 October 2017, the European Space Agency (ESA) successfully launched the Sentinel-5 Precursor satellite with the Tropospheric Monitoring Instrument (TROPOMI) as its single payload. TROPOMI is the first of ESA's atmospheric composition Sentinel missions, which will provide complete long-term records of atmospheric trace gases for the coming 30 years as a contribution to the European Union's Earth Observing program Copernicus. One of TROPOMI's primary products is atmospheric carbon monoxide (CO). It is observed with daily global coverage and a high spatial resolution of 7×7 km2. The moderate atmospheric resistance time and the low background concentration leads to localized pollution hotspots of CO and allows the tracking of the atmospheric transport of pollution on regional to global scales. In this contribution, we demonstrate the groundbreaking performance of the TROPOMI CO product, sensing CO enhancements above cities and industrial areas and tracking, with daily coverage, the atmospheric transport of pollution from biomass burning regions. The CO data product is validated with two months of Fourier-transform spectroscopy (FTS) measurements at nine ground-based stations operated by the Total Carbon Column Observing Network (TCCON). We found a good agreement between both datasets with a mean bias of 6 ppb (average of individual station biases) for both clear-sky and cloudy TROPOMI CO retrievals. Together with the corresponding standard deviation of the individual station biases of 3.8 ppb for clear-sky and 4.0 ppb for cloudy sky, it indicates that the CO data product is already well within the mission requirement.

A new data set of vertical column densities of the water vapour isotopologues H2O and HDO from the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) instrument for the whole of the mission period from January 2003 to April 2012 is presented. The data are retrieved from reflectance measurements in the spectral range 2339 to 2383 nm with the Shortwave Infrared CO Retrieval (SICOR) algorithm, ignoring atmospheric light scattering in the measurement simulation. The retrievals are validated with ground-based Fourier transform infrared measurements obtained within the Multi-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) project. A good agreement for low-altitude stations is found with an average bias of −3.6×1021 for H2O and −1.0×1018 molec cm−2 for HDO. The a posteriori computed δD shows an average bias of −8 ‰, even though polar stations have a larger negative bias. The latter is due to the large amount of sensor noise in SCIAMACHY in combination with low albedo and high solar zenith angles. To demonstrate the benefit of accounting for light scattering in the retrieval, the quality of the data product fitting effective cloud parameters simultaneously with trace gas columns is evaluated in a dedicated case study for measurements round high-altitude stations. Due to a large altitude difference between the satellite ground pixel and the mountain station, clear-sky scenes yield a large bias, resulting in a δD bias of 125 ‰. When selecting scenes with optically thick clouds within 1000 m above or below the station altitude, the bias in a posteriori δD is reduced from 125 to 44 ‰. The insights from the present study will also benefit the analysis of the data from the new Sentinel-5 Precursor mission.

In the perspective of the upcoming TROPOMI Sentinel-5 Precursor carbon monoxide data product, we discuss the benefit of using CO total column retrievals from cloud-contaminated SCIAMACHY 2.3 µm shortwave infrared spectra to detect atmospheric CO enhancements on regional and urban scales due to emissions from cities and wildfires. The study uses the operational Sentinel-5 Precursor algorithm SICOR, which infers the vertically integrated CO column together with effective cloud parameters. We investigate its capability to detect localized CO enhancements distinguishing between clear-sky observations and observations with low (

We discuss the retrieval of carbon monoxide (CO) vertical column densities from clear-sky and cloud contaminated 2311–2338 nm reflectance spectra measured by the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) from January 2003 until the end of the mission in April 2012. These data were processed with the Shortwave Infrared CO Retrieval algorithm (SICOR) that we developed for the operational data processing of the Tropospheric Monitoring Instrument (TROPOMI) that will be launched on ESA's Sentinel-5 Precursor (S5P) mission. This study complements previous work that was limited to clear-sky observations over land. Over the oceans, CO is estimated from cloudy-sky measurements only, which is an important addition to the SCIAMACHY clear-sky CO data set as shown by NDACC and TCCON measurements at coastal sites. For Ny-Ålesund, Lauder, Mauna Loa and Reunion, a validation of SCIAMACHY clear-sky retrievals is not meaningful because of the high retrieval noise and the few collocations at these sites. The situation improves significantly when considering cloudy-sky observations, where we find a low mean bias b = ±6. 0 ppb and a strong correlation between the validation and the SCIAMACHY results with a mean Pearson correlation coefficient r = 0. 7. Also for land observations, cloudy-sky CO retrievals present an interesting complement to the clear-sky data set. For example, at the cities Tehran and Beijing the agreement of SCIAMACHY clear-sky CO observations with MOZAIC/IAGOS airborne measurements is poor with a mean bias of b = 171. 2 ppb and 57.9 ppb because of local CO pollution, which cannot be captured by SCIAMACHY. For cloudy-sky retrievals, the validation improves significantly. Here the retrieved column is mainly sensitive to CO above the cloud and so not affected by the strong local surface emissions. Adjusting the MOZAIC/IAGOS measurements to the vertical sensitivity of the retrieval, the mean bias adds up to b = 52. 3 ppb and 5.0 ppb for Tehran and Beijing. At the less urbanised region around the airport Windhoek, local CO pollution is less prominent and so MOZAIC/IAGOS measurements agree well with SCIAMACHY clear-sky retrievals with a mean bias of b = 15. 5 ppb, but can be even further improved for cloudy SCIAMACHY observations with a mean bias of b = 0. 2 ppb. Overall the cloudy-sky CO retrievals from SCIAMACHY short-wave infrared measurements present a major extension of the clear-sky-only data set, which more than triples the amount of data and adds unique observations over the oceans. Moreover, the study represents the first application of the S5P algorithm for operational CO data processing on cloudy observations prior to the launch of the S5P mission.

This work presents the operational methane retrieval algorithm for the Sentinel 5 Precursor (S5P) satellite and its performance tested on realistic ensembles of simulated measurements. The target product is the column-averaged dry air volume mixing ratio of methane (XCH4), which will be retrieved simultaneously with scattering properties of the atmosphere. The algorithm attempts to fit spectra observed by the shortwave and near-infrared channels of the TROPOspheric Monitoring Instrument (TROPOMI) spectrometer aboard S5P.The sensitivity of the retrieval performance to atmospheric scattering properties, atmospheric input data and instrument calibration errors is evaluated. In addition, we investigate the effect of inhomogeneous slit illumination on the instrument spectral response function. Finally, we discuss the cloud filters to be used operationally and as backup.We show that the required accuracy and precision of 4 product are met for clear-sky measurements over land surfaces and after appropriate filtering of difficult scenes. The algorithm is very stable, having a convergence rate of 99 %. The forward model error is less than 1 % for about 95 % of the valid retrievals. Model errors in the input profile of water do not influence the retrieval outcome noticeably. The methane product is expected to meet the requirements if errors in input profiles of pressure and temperature remain below 0.3 % and 2 K, respectively. We further find that, of all instrument calibration errors investigated here, our retrievals are the most sensitive to an error in the instrument spectral response function of the shortwave infrared channel.

The Tropospheric Monitoring Instrument (TROPOMI) spectrometer is the single payload of the Copernicus Sentinel 5 Precursor (S5P) mission. It measures Earth radiance spectra in the shortwave infrared spectral range around 2.3 µm with a dedicated instrument module. These measurements provide carbon monoxide (CO) total column densities over land, which for clear sky conditions are highly sensitive to the tropospheric boundary layer. For cloudy atmospheres over land and ocean, the column sensitivity changes according to the light path through the atmosphere. In this study, we present the physics-based operational S5P algorithm to infer atmospheric CO columns satisfying the envisaged accuracy (

The TROPOspheric Monitoring Instrument (TROPOMI) on board the European Space Agency Sentinel-5 Precursor mission is scheduled for launch in the last quarter of 2016. As part of its operational processing the mission will provide CH4 and CO total columns using backscattered sunlight in the shortwave infrared band (2.3 µm). By adapting the CO retrieval algorithm, we have developed a non-scattering algorithm to retrieve total column HDO and H2O from the same measurements under clear-sky conditions. The isotopologue ratio HDO ∕ H2O is a powerful diagnostic in the efforts to improve our understanding of the hydrological cycle and its role in climate change, as it provides an insight into the source and transport history of water vapour, nature's strongest greenhouse gas. Due to the weak reflectivity over water surfaces, we need to restrict the retrieval to cloud-free scenes over land. We exploit a novel 2-band filter technique, using strong vs. weak water or methane absorption bands, to prefilter scenes with medium-to-high-level clouds, cirrus or aerosol and to significantly reduce processing time. Scenes with cloud top heights ≲1 km, very low fractions of high-level clouds or an aerosol layer above a high surface albedo are not filtered out. We use an ensemble of realistic measurement simulations for various conditions to show the efficiency of the cloud filter and to quantify the performance of the retrieval. The single-measurement precision in terms of δD is better than 15–25 ‰ for even the lowest surface albedo (2–4 ‰ for high albedos), while a small bias remains possible of up to ∼ 20 ‰ due to remaining aerosol or up to ∼ 70 ‰ due to remaining cloud contamination. We also present an analysis of the sensitivity towards prior assumptions, which shows that the retrieval has a small but significant sensitivity to the a priori assumption of the atmospheric trace gas profiles. Averaging multiple measurements over time and space, however, will reduce these errors, due to the quasi-random nature of the profile uncertainties. The sensitivity of the retrieval with respect to instrumental parameters within the expected instrument performance is 2, we will greatly increase the amount of useful data on HDO, H2O and their ratio HDO ∕ H2O. We showcase the overall performance of the retrieval algorithm and cloud filter with an accurate simulation of TROPOMI measurements from a single overpass over parts of the USA and Mexico, based on MODIS satellite data and realistic conditions for the surface, atmosphere and chemistry (including isotopologues). This shows that TROPOMI will pave the way for new studies of the hydrological cycle, both globally and locally, on timescales of mere days and weeks instead of seasons and years and will greatly extend the HDO ∕ H2O datasets from the SCIAMACHY and GOSAT missions.

We present a full-mission data product of carbon monoxide (CO) vertical column densities using the 2310–2338 nm SCIAMACHY reflectance measurements over clear-sky land scenes for the period January 2003–April 2012. The retrieval employs the SICOR algorithm, which will be used for operational data processing of the Sentinel-5 Precursor mission. The retrieval approach infers simultaneously carbon monoxide, methane and water vapour column densities together with a Lambertian surface albedo from individual SCIAMACHY measurements employing a non-scattering radiative transfer model. To account for the radiometric instrument degradation including the formation of an ice-layer on the 2.3 µm detector array, we consider clear-sky measurements over the Sahara as a natural calibration target. For these specific measurements, we spectrally calibrate the SCIAMACHY measurements and determine a spectral radiometric offset and the width of the instrument spectral response function as a function of time for the entire operational phase of the mission. We show that the smoothing error of individual clear-sky CO retrievals is less than ±1 ppb and thus this error contribution does not need to be accounted for in the validation considering the much higher retrieval noise. The CO data product is validated against measurements of ground-based Fourier transform infrared spectrometers at 27 stations of the NDACC-IRWG and TCCON network and MOZAIC/IAGOS aircraft measurements at 26 airports worldwide. Overall, we find a good agreement with TCCON measurements with a mean bias b = −1.2 ppb and a station-to-station bias with σ = 7.2 ppb. The negative sign of the bias means a low bias of SCIAMACHY CO with respect to TCCON. For the NDACC-IRWG network, we obtain a larger mean station bias of b = −9.2 ppb with σ = 8.1 ppb and for the MOZAIC/IAGOS measurements we find b = −6.4 ppb with σ = 5.6 ppb. The SCIAMACHY data set is subject to a small but significant bias trend of 1.47 ± 0.25 ppb yr−1. After trend correction, the bias with respect to MOZAIC/IAGOS observation is 2.5 ppb, with respect to TCCON measurements it is −4.6 ppb and with respect to NDACC-IRWG measurements −8.4 ppb. Hence, a discrepancy of 3.8 ppb remains between the global biases with NDACC-IRWG and TCCON, which is confirmed by directly comparing NDACC-IRWG and TCCON measurements. Generally, the scatter of the individual SCIAMACHY CO retrievals is high and dominated by large measurement noise. Hence, for practical usage of the data set, averaging of individual retrievals is required. As an example, we show that monthly mean SCIAMACHY CO retrievals, averaged separately over Northern and Southern Africa, reflect the spatial and temporal variability of biomass burning events in agreement with the global chemical transport model TM5.

We present a sensitivity study of the direct fitting approach to retrieve total ozone columns from the clear sky Global Ozone Monitoring Experiment 2/MetOp-A (GOME-2/MetOp-A) measurements between 325 and 335 nm in the period 2007–2010. The direct fitting of the measurement is based on adjusting the scaling of a reference ozone profile and requires accurate simulation of GOME-2 radiances. In this context, we study the effect of three aspects that introduce forward model errors if not addressed appropriately: (1) the use of a clear sky model atmosphere in the radiative transfer demanding cloud filtering, (2) different approximations of Earth's sphericity to address the influence of the solar zenith angle, and (3) the need of polarization in radiative transfer modeling. We conclude that cloud filtering using the operational GOME-2 FRESCO (Fast Retrieval Scheme for Clouds from the Oxygen A band) cloud product, which is part of level 1B data, and the use of pseudo-spherical scalar radiative transfer is fully sufficient for the purpose of this retrieval. A validation with ground-based measurements at 36 stations confirms this showing a global mean bias of −0.1 % with a standard deviation (SD) of 2.7 %. The regularization effect inherent to the profile scaling approach is thoroughly characterized by the total column averaging kernel for each individual retrieval. It characterizes the effect of the particular choice of the ozone profile to be scaled by the inversion and is part of the retrieval product. Two different interpretations of the data product are possible: first, regarding the retrieval product as an estimate of the true column, a direct comparison of the retrieved column with total ozone columns from ground-based measurements can be done. This requires accurate a priori knowledge of the reference ozone profile and the column averaging kernel is not needed. Alternatively, the retrieval product can be interpreted as an effective column defined by the total column averaging kernel. This interpretation relies much less on the a priori knowledge of the reference ozone profile; however, for its validation, measurements of the vertical ozone distribution are needed. The different manners of data interpretation are demonstrated for simulated and real measurements using on-ground ozone column and ozonesonde measurements for validation.

Insights are given into Tikhonov regularization and its application to the retrieval of vertical column densities of atmospheric trace gases from remote sensing measurements. The study builds upon the equivalence of the least-squares profile-scaling approach and Tikhonov regularization method of the first kind with an infinite regularization strength. Here, the vertical profile is expressed relative to a reference profile. On the basis of this, we propose a new algorithm as an extension of the least-squares profile scaling which permits the calculation of total column averaging kernels on arbitrary vertical grids using an analytic expression. Moreover, we discuss the effective null space of the retrieval, which comprises those parts of a vertical trace gas distribution which cannot be inferred from the measurements. Numerically the algorithm can be implemented in a robust and efficient manner. In particular for operational data processing with challenging demands on processing time, the proposed inversion method in combination with highly efficient forward models is an asset. For demonstration purposes, we apply the algorithm to CO column retrieval from simulated measurements in the 2.3 μm spectral region and to O3 column retrieval from the UV. These represent ideal measurements of a series of spaceborne spectrometers such as SCIAMACHY, TROPOMI, GOME, and GOME-2. For both spectral ranges, we consider clear-sky and cloudy scenes where clouds are modelled as an elevated Lambertian surface. Here, the smoothing error for the clear-sky and cloudy atmosphere is significant and reaches several percent, depending on the reference profile which is used for scaling. This underlines the importance of the column averaging kernel for a proper interpretation of retrieved column densities. Furthermore, we show that the smoothing due to regularization can be underestimated by calculating the column averaging kernel on a too coarse vertical grid. For both retrievals, this effect becomes negligible for a vertical grid with 20–40 equally thick layers between 0 and 50 km.

Trends in the CO andC2H6 partial columns ~0–15 km) have been estimated from four European ground-based solar FTIR (Fourier Transform InfraRed) stations for the 1996–2006 time period. The CO trends from the four stations Jungfraujoch, Zugspitze, Harestua and Kiruna have been estimated to −0.45 ± 0.16% yr−1, −1.00 ± 0.24% yr−1, −0.62 ± 0.19 % yr−1 and −0.61 ± 0.16% yr−1, respectively. The corresponding trends for C2H6 are −1.51 ± 0.23% yr−1, −2.11 ± 0.30% yr−1, −1.09 ± 0.25% yr−1 and −1.14 ± 0.18% yr−1. All trends are presented with their 2-σ confidence intervals. To find possible reasons for the CO trends, the global-scale EMEP MSC-W chemical transport model has been used in a series of sensitivity scenarios. It is shown that the trends are consistent with the combination of a 20% decrease in the anthropogenic CO emissions seen in Europe and North America during the 1996–2006 period and a 20% increase in the anthropogenic CO emissions in East Asia, during the same time period. The possible impacts of CH4 and biogenic volatile organic compounds (BVOCs) are also considered. The European and global-scale EMEP models have been evaluated against the measured CO and C2H6 partial columns from Jungfraujoch, Zugspitze, Bremen, Harestua, Kiruna and Ny-Ålesund. The European model reproduces, on average the measurements at the different sites fairly well and within 10–22% deviation for CO and 14–31% deviation for C2H6. Their seasonal amplitude is captured within 6–35% and 9–124% for CO and C2H6, respectively. However, 61–98% of the CO and C2H6 partial columns in the European model are shown to arise from the boundary conditions, making the global-scale model a more suitable alternative when modeling these two species. In the evaluation of the global model the average partial columns for 2006 are shown to be within 1–9% and 37–50% of the measurements for CO and C2H6, respectively. The global model sensitivity for assumptions made in this paper is also analyzed.

Total columns measured with the ground-based solar FTIR technique are highly variable in time due to atmospheric chemistry and dynamics in the atmosphere above the measurement station. In this paper, a multiple regression model with anomalies of air pressure, total columns of hydrogen fluoride (HF) and carbon monoxide (CO) and tropopause height are used to reduce the variability in the methane (CH4) and nitrous oxide (N2O) total columns to estimate reliable linear trends with as small uncertainties as possible. The method is developed at the Harestua station (60° N, 11° E, 600 m a.s.l.) and used on three other European FTIR stations, i.e. Jungfraujoch (47° N, 8° E, 3600 m a.s.l.), Zugspitze (47° N, 11° E, 3000 m a.s.l.), and Kiruna (68° N, 20° E, 400 m a.s.l.). Linear CH4 trends between 0.13 ± 0.01-0.25 ± 0.02 % yr−1 were estimated for all stations in the 1996-2009 period. A piecewise model with three separate linear trends, connected at change points, was used to estimate the short term fluctuations in the CH4 total columns. This model shows a growth in 1996–1999 followed by a period of steady state until 2007. From 2007 until 2009 the atmospheric CH4 amount increases between 0.57 ± 0.22–1.15 ± 0.17 % yr−1. Linear N2O trends between 0.19 ± 0.01–0.40 ± 0.02 % yr−1 were estimated for all stations in the 1996-2007 period, here with the strongest trend at Harestua and Kiruna and the lowest at the Alp stations. From the N2O total columns crude tropospheric and stratospheric partial columns were derived, indicating that the observed difference in the N2O trends between the FTIR sites is of stratospheric origin. This agrees well with the N2O measurements by the SMR instrument onboard the Odin satellite showing the highest trends at Harestua, 0.98 ± 0.28 % yr−1, and considerably smaller trends at lower latitudes, 0.27 ± 0.25 % yr−1. The multiple regression model was compared with two other trend methods, the ordinary linear regression and a Bootstrap algorithm. The multiple regression model estimated CH4 and N2O trends that differed up to 31 % compared to the other two methods and had uncertainties that were up to 300 % lower. Since the multiple regression method were carefully validated this stresses the importance to account for variability in the total columns when estimating trend from solar FTIR data.

We present an intercomparison of three years of measurements of integrated water vapor (IWV) performed by the mid-infrared solar FTIR (Fourier Transform Infra-Red) instrument on the summit of Mt. Zugspitze (2964 m a.s.l.) and by the nearby near-infrared differential absorption lidar (DIAL) at the Schneefernerhaus research station (2675 m a.s.l.). The solar FTIR was shown to be one of the most accurate and precise IWV sounders in recent work (Sussmann et al., 2009) and is taken as the reference here. By calculating the FTIR-DIAL correlation (22 min coincidence interval, 15 min integration time) we derive an almost ideal slope of 0.996 (10), a correlation coefficient of R = 0.99, an IWV intercept of −0.039 (42) mm (−1.2 % of the mean), and a bias of −0.052 (26) mm (−1.6 % of the mean) from the scatter plot. By selecting a subset of coincidences with an optimum temporal and spatial matching between DIAL and FTIR, we obtain a conservative estimate of the precision of the DIAL in measuring IWV which is better than 0.1 mm (3.2 % of the mean). We found that for a temporal coincidence interval of 22 min the difference in IWV measured by these two systems is dominated by the volume mismatch (horizontal distance: 680 m). The outcome from this paper is twofold: (1) the IWV soundings by FTIR and DIAL agree very well in spite of the differing wavelength regions with different spectroscopic line parameters and retrieval algorithms used. (2) In order to derive an estimate of the precision of state-of-the-art IWV sounders from intercomparison experiments, it is necessary to use a temporal matching on time scales shorter than 10 min and a spatial matching on the 100-m scale.

This paper presents a validation study of SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) carbon monoxide (CO) total column measurements from the Iterative Maximum Likelihood Method (IMLM) algorithm using ground-based spectrometer observations from twenty surface stations for the five year time period of 2003–2007. Overall we find a good agreement between SCIAMACHY and ground-based observations for both mean values as well as seasonal variations. For high-latitude Northern Hemisphere stations absolute differences between SCIAMACHY and ground-based measurements are close to or fall within the SCIAMACHY CO 2σ precision of 0.2 × 1018 molecules/cm2 (∼10%) indicating that SCIAMACHY can observe CO accurately at high Northern Hemisphere latitudes. For Northern Hemisphere mid-latitude stations the validation is complicated due to the vicinity of emission sources for almost all stations, leading to higher ground-based measurements compared to SCIAMACHY CO within its typical sampling area of 8° × 8°. Comparisons with Northern Hemisphere mountain stations are hampered by elevation effects. After accounting for these effects, the validation provides satisfactory results. At Southern Hemisphere mid- to high latitudes SCIAMACHY is systematically lower than the ground-based measurements for 2003 and 2004, but for 2005 and later years the differences between SCIAMACHY and ground-based measurements fall within the SCIAMACHY precision. The 2003–2004 bias is consistent with previously reported results although its origin remains under investigation. No other systematic spatial or temporal biases could be identified based on the validation presented in this paper. Validation results are robust with regard to the choices of the instrument-noise error filter, sampling area, and time averaging required for the validation of SCIAMACHY CO total column measurements. Finally, our results show that the spatial coverage of the ground-based measurements available for the validation of the 2003–2007 SCIAMACHY CO columns is sub-optimal for validation purposes, and that the recent and ongoing expansion of the ground-based network by carefully selecting new locations may be very beneficial for SCIAMACHY CO and other satellite trace gas measurements validation efforts.

We present a method for harmonized retrieval of integrated water vapor (IWV) from existing, long-term, measurement records at the ground-based mid-infrared solar FTIR spectrometry stations of the Network for the Detection of Atmospheric Composition Change (NDACC). Correlation of IWV from FTIR with radiosondes shows an ideal slope of 1.00(3). This optimum matching is achieved via tuning one FTIR retrieval parameter, i.e., the strength of a Tikhonov regularization constraining the derivative (with respect to height) of retrieved water profiles given in per cent difference relative to an a priori profile. All other FTIR-sonde correlation parameters (intercept=0.02(12) mm, bias=0.02(5) mm, standard deviation of coincident IWV differences (stdv)=0.27 mm, R=0.99) are comparable to or better than results for all other ground-based IWV sounding techniques given in the literature. An FTIR-FTIR side-by-side intercomparison reveals a strong exponential increase in stdv as a function of increasing temporal mismatch starting at Δt≈1 min. This is due to atmospheric water vapor variability. Based on this result we derive an upper limit for the precision of the FTIR IWV retrieval for the smallest Δt(=3.75 min) still giving a statistically sufficient sample (32 coincidences), i.e., precision(IWVFTIR) As first application examples long-term IWV data sets are retrieved from the FTIR records of the Zugspitze (47.4° N, 11.0° E, 2964 m a.s.l.) and Jungfraujoch (46.5° N, 8.0° E, 3580 m a.s.l.) NDACC sites. Station-trend analysis comprises a linear fit after subtracting an intra-annual model (3 Fourier components) and constructing an uncertainty interval [95% confidence] via bootstrap resampling. For the Zugspitze a significant trend of 0.79 [0.65, 0.92] mm/decade is found for the time interval [1996–2008], whereas for the Jungfraujoch no significant trend is found. This confirms recent findings that strong variations of IWV trends do occur above land on the local to regional scale (≈250 km) in spite of homogeneous surface temperature trends. This paper provides a basis for future exploitation of more than a dozen existing, multi-decadal FTIR measurement records around the globe for climate studies.

This paper presents extensive {bias determination} analyses of ozone observations from the Atmospheric Chemistry Experiment (ACE) satellite instruments: the ACE Fourier Transform Spectrometer (ACE-FTS) and the Measurement of Aerosol Extinction in the Stratosphere and Troposphere Retrieved by Occultation (ACE-MAESTRO) instrument. Here we compare the latest ozone data products from ACE-FTS and ACE-MAESTRO with coincident observations from nearly 20 satellite-borne, airborne, balloon-borne and ground-based instruments, by analysing volume mixing ratio profiles and partial column densities. The ACE-FTS version 2.2 Ozone Update product reports more ozone than most correlative measurements from the upper troposphere to the lower mesosphere. At altitude levels from 16 to 44 km, the average values of the mean relative differences are nearly all within +1 to +8%. At higher altitudes (45–60 km), the ACE-FTS ozone amounts are significantly larger than those of the comparison instruments, with mean relative differences of up to +40% (about +20% on average). For the ACE-MAESTRO version 1.2 ozone data product, mean relative differences are within ±10% (average values within ±6%) between 18 and 40 km for both the sunrise and sunset measurements. At higher altitudes (~35–55 km), systematic biases of opposite sign are found between the ACE-MAESTRO sunrise and sunset observations. While ozone amounts derived from the ACE-MAESTRO sunrise occultation data are often smaller than the coincident observations (with mean relative differences down to −10%), the sunset occultation profiles for ACE-MAESTRO show results that are qualitatively similar to ACE-FTS, indicating a large positive bias (mean relative differences within +10 to +30%) in the 45–55 km altitude range. In contrast, there is no significant systematic difference in bias found for the ACE-FTS sunrise and sunset measurements.