Search

Instruments for air quality observations on geostationary satellites provide multiple observations per day and allow for the analysis of the diurnal variation in important air pollutants such as nitrogen dioxide (NO2). The South Korean instrument GEMS (Geostationary Environmental Monitoring Spectrometer), launched in February 2020, is the first geostationary instrument that is able to observe the diurnal variation in NO2. The measurements have a spatial resolution of 3.5 km × 8 km and cover a large part of Asia. This study compares 1 year of tropospheric NO2 vertical column density (VCD) observations from the operational GEMS L2 product, the scientific GEMS IUP-UB (Institute of Environmental Physics at the University of Bremen) product, the operational TROPOspheric Monitoring Instrument (TROPOMI) product, and ground-based differential optical absorption spectroscopy (DOAS) measurements in South Korea. The GEMS L2 tropospheric NO2 VCDs overestimate the ground-based tropospheric NO2 VCDs with a median relative difference of +61 % and a correlation coefficient of 0.76. The median relative difference is −2 % for the GEMS IUP-UB product and −16 % for the TROPOMI product, with correlation coefficients of 0.83 and 0.89, respectively. The scatter in the GEMS products can be reduced when observations are limited to the TROPOMI overpass time. Diurnal variations in tropospheric NO2 VCDs differ by the pollution level of the analyzed site but with good agreement between the GEMS IUP-UB and ground-based observations. Low-pollution sites show weak or almost no diurnal variation. In summer, the polluted sites show a minimum around noon, indicating the large influence of photochemical loss. Most variation is seen in spring and autumn, with increasing NO2 in the morning, a maximum close to noon, and a decrease towards the afternoon. Winter observations show rather flat or slightly decreasing NO2 throughout the day. Winter observations under low-wind-speed conditions at high-pollution sites show enhancements of NO2 throughout the day. This indicates that under calm conditions, dilution and the less effective chemical loss in winter do not balance the accumulating emissions. Diurnal variation observed at a low-pollution site follows seasonal wind patterns. A weekday–weekend effect analysis shows good agreement between the different products. However, the GEMS L2 product, while agreeing with the other data sets on weekdays, shows significantly less reduction on weekends. The influence of the stratospheric contribution and the surface reflectivity product on the satellite tropospheric NO2 VCD products is investigated. While the TM5 model's stratospheric VCDs, used in the TROPOMI product, are too high, resulting in tropospheric NO2 VCDs that are too low and even negative, when used in the GEMS IUP-UB retrieval, the GEMS L2 stratospheric VCD is too low. Surface reflectivity comparisons indicate that the GEMS L2 reflectivity makes a large contribution to the observed overestimation and scatter.

Bromine monoxide (BrO) is relevant to atmospheric oxidative capacity, affecting the lifetime of greenhouse gases (i.e., methane, dimethylsulfide) and mercury oxidation. However, measurements of BrO radical vertical profiles are rare, and BrO is highly variable. As a result, the few available aircraft observations in different regions of the atmosphere are not easily reconciled. Autonomous multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments placed at remote mountaintop observatories (MT-DOAS) present a cost-effective alternative to aircraft, with the potential to probe the climate-relevant yet understudied free troposphere more routinely. Here, we describe an innovative full-atmosphere BrO and formaldehyde (HCHO) profile retrieval algorithm using MT-DOAS measurements at Mauna Loa Observatory (MLO – 19.536° N, 155.577° W; 3401 m a.s.l.). The retrieval is based on time-dependent optimal estimation and simultaneously inverts 190+ individual BrO (and formaldehyde, HCHO) SCDs (slant column densities; SCD = dSCD + SCDRef) from solar stray light spectra measured in the zenith and off-axis geometries at high and low solar zenith angles (92° > SZA > 30°) to derive BrO concentration profiles from 1.9 to 35 km with 7.5 degrees of freedom (DoFs). Two case study days are characterized by the absence (26 April 2017, base case) and presence of a Rossby-wave-breaking double tropopause (29 April 2017, RW-DT case). Stratospheric-BrO vertical columns are nearly identical on both days (VCD = (1.5 ± 0.2) × 1013 molec. cm−2), and the stratospheric-BrO profile peaks at a lower altitude during the RW-DT (1.6–2.0 DoFs). Tropospheric-BrO VCDs increase from (0.70 ± 0.14) × 1013 molec. cm−2 (base case) to (1.00 ± 0.14) × 1013 molec. cm−2 (RW-DT) owing to a 3-fold increase in BrO in the upper troposphere (1.7–1.9 DoFs). BrO at MLO increases from (0.23 ± 0.03) pptv (base case) to (0.46 ± 0.03) pptv (RW-DT) and is characterized by an added time resolution (∼ 3.8 DoFs). Up to (0.9 ± 0.1) pptv BrO is observed above MLO in the lower free troposphere in the absence of the double tropopause. We validate the retrieval using aircraft BrO profiles and in situ HCHO measurements aboard the NSF/NCAR GV aircraft above MLO (11 January 2014) that establish BrO peaks around 2.4 pptv above 13 km in the upper troposphere–lower stratosphere (UTLS) during a similar RW-DT event (0.83 × 1013 molec. cm2 tropospheric-BrO VCD above 2 km). The tropospheric-BrO profile measured using MT-DOAS (RW-DT case) and using the aircraft agree well (after averaging-kernel smoothing). Furthermore, these tropospheric-BrO profiles over the central Pacific Ocean are found to closely resemble those over the eastern Pacific Ocean (2–14 km) and are in contrast to those over the western Pacific Ocean, where a C-shaped tropospheric-BrO profile shape has been observed.

The Geostationary Environment Monitoring Spectrometer (GEMS) on board GEO-KOMPSAT-2B was launched in February 2020 and has been monitoring atmospheric chemical compositions over Asia. We present the first evaluation of the operational GEMS formaldehyde (HCHO) vertical column densities (VCDs) during and after the in-orbit test (IOT) period (August–October 2020) by comparing them with the products from the TROPOspheric Monitoring Instrument (TROPOMI) and Fourier-transform infrared (FTIR) and multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments. During the IOT, the GEMS HCHO VCDs reproduced the observed spatial pattern of TROPOMI VCDs over the entire domain (r= 0.62) with high biases (10 %–16 %). We found that the agreement between GEMS and TROPOMI was substantially higher in Northeast Asia (r= 0.90), encompassing the Korean Peninsula and east China. GEMS HCHO VCDs captured the seasonal variation in HCHO, primarily driven by biogenic emissions and photochemical activities, but showed larger variations than those of TROPOMI over coastal regions (Kuala Lumpur, Singapore, Shanghai, and Busan). In addition, GEMS HCHO VCDs showed consistent hourly variations with MAX-DOAS (r= 0.77) and FTIR (r= 0.86) but were 30–40 % lower than ground-based observations. Different vertical sensitivities of GEMS and ground-based instruments caused these biases. Utilizing the averaging kernel smoothing method reduces the low biases by approximately 10 % to 15 % (normalized mean bias (NMB): −47.4 % to −31.5 % and −38.6 % to −26.7 % for MAX-DOAS and FTIR, respectively). The remaining discrepancies are due to multiple factors, including spatial collocation and different instrumental sensitivities, requiring further investigation using inter-comparable datasets.

Volatile organic compounds (VOCs) are key precursors of particulate matter and tropospheric ozone. Although the terrestrial biosphere is by far the largest source of VOCs into the atmosphere, the emissions of biogenic VOCs remain poorly constrained at the regional scale. In this work, we derive top-down biogenic emissions over Europe using weekly averaged TROPOMI formaldehyde (HCHO) data from 2018 to 2021. The systematic bias of the TROPOMI HCHO columns is characterized and corrected for based on comparisons with FTIR data at seven European stations. The top-down fluxes of biogenic, pyrogenic, and anthropogenic VOC sources are optimized using an inversion framework based on the MAGRITTEv1.1 chemistry transport model and its adjoint. The inversion leads to strongly increased isoprene emissions with respect to the MEGAN–MOHYCAN inventory over the model domain (from 8.1 to 18.5 Tg yr−1), which is driven by the high observed TROPOMI HCHO columns in southern Europe. The impact of the inversion on biomass burning VOCs (+13 %) and anthropogenic VOCs (−17 %) is moderate. An evaluation of the optimized HCHO distribution against ground-based remote sensing (FTIR and MAX-DOAS) and in situ data provides generally improved agreement at stations below about 50∘ N but indicates overestimated emissions in northern Scandinavia. Sensitivity inversions show that the top-down emissions are robust with respect to changes in the inversion settings and in the model chemical mechanism, leading to differences of up to 10 % in the total emissions. However, the top-down emissions are very sensitive to the bias correction of the observed columns, as the biogenic emissions are 3 times lower when the correction is not applied. Furthermore, the use of different a priori biogenic emissions has a significant impact on the inversion results due to large differences among bottom-up inventories. The sensitivity run using CAMS-GLOB-BIOv3.1 as a priori emissions in the inversion results in 30 % lower emissions with respect to the optimization using MEGAN–MOHYCAN. In regions with large temperature and cloud cover variations, there is strong week-to-week variability in the observed HCHO columns. The top-down emissions, which are optimized at weekly increments, have a much improved capability of representing these large fluctuations than an inversion using monthly increments.

Ozone depletion over the polar regions is monitored each year by satellite- and ground-based instruments. In this study, the vortex-averaged ozone loss over the last 3 decades is evaluated for both polar regions using the passive ozone tracer of the chemical transport model TOMCAT/SLIMCAT and total ozone observations from Système d'Analyse par Observation Zénithale (SAOZ) ground-based instruments and Multi-Sensor Reanalysis (MSR2). The passive-tracer method allows us to determine the evolution of the daily rate of column ozone destruction and the magnitude of the cumulative column loss at the end of the winter. Three metrics are used in trend analyses that aim to assess the ozone recovery rate over both polar regions: (1) the maximum ozone loss at the end of the winter, (2) the onset day of ozone loss at a specific threshold, and (3) the ozone loss residuals computed from the differences between annual ozone loss and ozone loss values regressed with respect to sunlit volume of polar stratospheric clouds (VPSCs). This latter metric is based on linear and parabolic regressions for ozone loss in the Northern Hemisphere and Southern Hemisphere, respectively. In the Antarctic, metrics 1 and 3 yield trends of −2.3 % and −2.2 % per decade for the 2000–2021 period, significant at 1 and 2 standard deviations (σ), respectively. For metric 2, various thresholds were considered at the total ozone loss values of 20 %, 25 %, 30 %, 35 %, and 40 %, all of them showing a time delay as a function of year in terms of when the threshold is reached. The trends are significant at the 2σ level and vary from 3.5 to 4.2 d per decade between the various thresholds. In the Arctic, metric 1 exhibits large interannual variability, and no significant trend is detected; this result is highly influenced by the record ozone losses in 2011 and 2020. Metric 2 is not applied in the Northern Hemisphere due to the difficulty in finding a threshold value in enough of the winters. Metric 3 provides a negative trend in Arctic ozone loss residuals with respect to the sunlit VPSC fit of −2.00 ± 0.97 (1σ) % per decade, with limited significance at the 2σ level. With such a metric, a potential quantitative detection of ozone recovery in the Arctic springtime lower stratosphere can be made.

We present a database of MAX-DOAS (Multi-AXis Differential Optical Absorption Spectroscopy) ground-based observations of NO2 and H2CO tropospheric vertical column densities (VCDtropo) performed for the first time in the city of Kinshasa. These measurements were conducted between November 2019 and July 2021 and processed using the standardized inversion tools developed in the ESA FRM4DOAS (Fiducial Reference Measurements for Ground-Based DOAS Air-Quality Observations) project. The retrieved geophysical quantities are used to validate column observations from the TROPOspheric Monitoring Instrument (TROPOMI) over Kinshasa. In the validation, we experiment with three different comparison cases of increasing complexity. In the first case, a direct comparison between MAX-DOAS observations (hourly average of MAX-DOAS VCDtropo at overpass) and TROPOMI shows an underestimation of TROPOMI with a median bias of −38 % for NO2 and −39 % for H2CO based on monthly comparison. The second case takes into account the different vertical sensitivities of the two instruments and the a priori profile. We note significant changes in the median bias for both compounds: −12 % for NO2 and +11 % for H2CO. The third case builds on the second case by considering also the direction of sight of the MAX-DOAS. For this third case, we find a median bias of +44 % for NO2 and a median bias of +4 % for H2CO. However this case is impacted by low sampling and is considered less reliable. The findings from this study underscore the significance of employing a realistic a priori profile in TROPOMI column extraction, particularly within heavily polluted urban zones like Kinshasa. The investigation also highlights the necessity for prudence when integrating the MAX-DOAS line of sight due to the noise generated during subsampling and the limited horizontal sensitivity of MAX-DOAS observations. Importantly, the study further reveals the pronounced pollution levels of NO2, H2CO and aerosols in both the city of Kinshasa and its adjacent regions, underscoring the imperative for consistent monitoring and effective regulatory measures by local authorities.

Tropospheric vertical column densities (VCDs) of nitrogen dioxide (NO2) retrieved from sun-synchronous satellite instruments have provided abundant NO2 data for environmental studies, but such data are limited by retrieval uncertainties and insufficient temporal sampling (e.g., once a day). The Geostationary Environment Monitoring Spectrometer (GEMS) launched in February 2020 monitors NO2 at an unprecedented hourly resolution during the daytime. Here we present a research product for tropospheric NO2 VCDs, referred to as POMINO–GEMS (where POMINO is the Peking University OMI NO2 algorithm). We develop a hybrid retrieval method combining GEMS, TROPOMI (TROPOspheric Monitoring Instrument) and GEOS-CF (Global Earth Observing System Composition Forecast) data to generate hourly tropospheric NO2 slant column densities (SCDs). We then derive tropospheric NO2 air mass factors (AMFs) with explicit corrections for surface reflectance anisotropy and aerosol optical effects through parallelized pixel-by-pixel radiative transfer calculations. Prerequisite cloud parameters are retrieved with the O2–O2 algorithm by using ancillary parameters consistent with those used in NO2 AMF calculations. The initial retrieval of POMINO–GEMS tropospheric NO2 VCDs for June–August 2021 exhibits strong hotspot signals over megacities and distinctive diurnal variations over polluted and clean areas. POMINO–GEMS NO2 VCDs agree with the POMINO–TROPOMI v1.2.2 product (R=0.98; NMB = 4.9 %) over East Asia, with slight differences associated with satellite viewing geometries and cloud and aerosol properties affecting the NO2 retrieval. POMINO–GEMS also shows good agreement with the following: OMNO2 (Ozone Monitoring Instrument (OMI) NO2 Standard Product) v4 (R=0.87; NMB = −16.8 %); and GOME-2 (Global Ozone Monitoring Experiment-2) GDP (GOME Data Processor) 4.8 (R=0.83; NMB = −1.5 %) NO2 products. POMINO–GEMS shows small biases against ground-based MAX-DOAS (multi-axis differential optical absorption spectroscopy) NO2 VCD data at nine sites (NMB = −11.1 %), with modest or high correlation in diurnal variation at six urban and suburban sites (R from 0.60 to 0.96). The spatiotemporal variation in POMINO–GEMS correlates well with mobile car MAX-DOAS measurements in the Three Rivers source region on the Tibetan Plateau (R=0.81). Surface NO2 concentrations estimated from POMINO–GEMS VCDs are consistent with measurements from the Ministry of Ecology and Environment of China for spatiotemporal variation (R=0.78; NMB = −26.3 %) and diurnal variation at all, urban, suburban and rural sites (R≥0.96). POMINO–GEMS data will be made freely available for users to study the spatiotemporal variations, sources and impacts of NO2.

We introduce the new Global Ozone Monitoring Experiment-2 (GOME-2) daily and monthly level-3 product of total column ozone (O3), total and tropospheric column nitrogen dioxide (NO2), total column water vapour, total column bromine oxide (BrO), total column formaldehyde (HCHO), and total column sulfur dioxide (SO2) (daily products https://doi.org/10.15770/EUM_SAF_AC_0048, AC SAF, 2023a; monthly products https://doi.org/10.15770/EUM_SAF_AC_0049, AC SAF, 2023b). The GOME-2 level-3 products aim to provide easily translatable and user-friendly data sets to the scientific community for scientific progress as well as to satisfy public interest. The purpose of this paper is to present the theoretical basis as well as the verification and validation of the GOME-2 daily and monthly level-3 products. The GOME-2 level-3 products are produced using the overlapping area-weighting method. Details of the gridding algorithm are presented. The spatial resolution of the GOME-2 level-3 products is selected based on the sensitivity study. The consistency of the resulting level-3 products among three GOME-2 sensors is investigated through time series of global averages, zonal averages, and bias. The accuracy of the products is validated by comparison to ground-based observations. The verification and validation results show that the GOME-2 level-3 products are consistent with the level-2 data. Small discrepancies are found among three GOME-2 sensors, which are mainly caused by the differences in the instrument characteristic and level-2 processor. The comparison of GOME-2 level-3 products to ground-based observations in general shows very good agreement, indicating that the products are consistent and fulfil the requirements to serve the scientific community and general public.

In this paper, we present the updated SAGE-CCI-OMPS+ climate data record of monthly zonal mean ozone profiles. This dataset covers the stratosphere and combines measurements by nine limb and occultation satellite instruments – SAGE II (Stratospheric Aerosol and Gases Experiment II), OSIRIS (Optical Spectrograph and InfraRed Imaging System), MIPAS (Michelson Interferometer for Passive Atmospheric Sounding), SCIAMACHY (SCanning Imaging Spectrometer for Atmospheric CHartographY), GOMOS (Global Ozone Monitoring by Occultation of Stars), ACE-FTS (Atmospheric Chemistry Experiment Fourier Transform Spectrometer), OMPS-LP (Ozone Monitor Profiling Suite Limb Profiler), POAM (Polar Ozone and Aerosol Measurement) III, and SAGE III/ISS (Stratospheric Aerosol and Gases Experiment III on the International Space Station). Compared to the original version of the SAGE-CCI-OMPS dataset (Sofieva et al., 2017b), the update includes new versions of MIPAS, ACE-FTS, and OSIRIS datasets and introduces data from additional sensors (POAM III and SAGE III/ISS) and retrieval processors (OMPS-LP). In this paper, we show detailed intercomparisons of ozone profiles from different instruments and data versions, with a focus on the detection of possible drifts in the datasets. The SAGE-CCI-OMPS+ dataset has a better coverage of polar regions and of the upper troposphere and the lower stratosphere (UTLS) than the previous dataset. We also studied the influence of including new datasets on ozone trends, which are estimated using multiple linear regression. The changes in the merged dataset do not change the overall morphology of post-1997 ozone trends; statistically significant trends are observed in the upper stratosphere. The largest changes in ozone trends are observed in polar regions, especially in the Southern Hemisphere. The updated SAGE-CCI-OMPS+ dataset contains profiles of deseasonalized anomalies and ozone concentrations from 1984 to 2021, in 10∘ latitude bins from 90∘ S to 90∘ N and in the altitude range from 10 to 50 km. The dataset is open access and available at https://climate.esa.int/en/projects/ozone/data/ (last access: 9 March 2023) and at ftp://cci_web@ftp-ae.oma.be/esacci (ESA Climate Office; last access: 9 March 2023).

Airborne imaging differential optical absorption spectroscopy (DOAS), ground-based stationary DOAS, and car DOAS measurements were conducted during the S5P-VAL-DE-Ruhr campaign in September 2020. The campaign area is located in the Rhine-Ruhr region of North Rhine-Westphalia, western Germany, which is a pollution hotspot in Europe comprising urban and large industrial sources. The DOAS measurements are used to validate spaceborne NO2 tropospheric vertical column density (VCD) data products from the Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI). Seven flights were performed with the airborne imaging DOAS instrument for measurements of atmospheric pollution (AirMAP), providing measurements that were used to create continuous maps of NO2 in the layer below the aircraft. These flights cover many S5P ground pixels within an area of 30 km × 35 km and were accompanied by ground-based stationary measurements and three mobile car DOAS instruments. Stationary measurements were conducted by two Pandora, two Zenith-DOAS, and two MAX-DOAS instruments. Ground-based stationary and car DOAS measurements are used to evaluate the AirMAP tropospheric NO2 VCDs and show high Pearson correlation coefficients of 0.88 and 0.89 and slopes of 0.90 ± 0.09 and 0.89 ± 0.02 for the stationary and car DOAS, respectively. Having a spatial resolution of about 100 m × 30 m, the AirMAP tropospheric NO2 VCD data create a link between the ground-based and the TROPOMI measurements with a nadir resolution of 3.5 km × 5.5 km and are therefore well suited to validate the TROPOMI tropospheric NO2 VCD. The observations on the 7 flight days show strong NO2 variability, which is dependent on the three target areas, the day of the week, and the meteorological conditions. The AirMAP campaign data set is compared to the TROPOMI NO2 operational offline (OFFL) V01.03.02 data product, the reprocessed NO2 data using the V02.03.01 of the official level-2 processor provided by the Product Algorithm Laboratory (PAL), and several scientific TROPOMI NO2 data products. The AirMAP and TROPOMI OFFL V01.03.02 data are highly correlated (r=0.87) but show an underestimation of the TROPOMI data with a slope of 0.38 ± 0.02 and a median relative difference of −9 %. With the modifications in the NO2 retrieval implemented in the PAL V02.03.01 product, the slope and median relative difference increased to 0.83 ± 0.06 and +20 %. However, the modifications resulted in larger scatter and the correlation decreased significantly to r=0.72. The results can be improved by not applying a cloud correction for the TROPOMI data in conditions with high aerosol load and when cloud pressures are retrieved close to the surface. The influence of spatially more highly resolved a priori NO2 vertical profiles and surface reflectivity are investigated using scientific TROPOMI tropospheric NO2 VCD data products. The comparison of the AirMAP campaign data set to the scientific data products shows that the choice of surface reflectivity database has a minor impact on the tropospheric NO2 VCD retrieval in the campaign region and season. In comparison, the replacement of the a priori NO2 profile in combination with the improvements in the retrieval of the PAL V02.03.01 product regarding cloud heights can further increase the tropospheric NO2 VCDs. This study demonstrates that the underestimation of the TROPOMI tropospheric NO2 VCD product with respect to the validation data set has been and can be further significantly improved.

Clouds have a strong impact on satellite measurements of tropospheric trace gases in the ultraviolet, visible, and near-infrared spectral ranges from space. Therefore, trace gas retrievals rely on information on cloud fraction, cloud albedo, and cloud height from cloud products. In this study, the cloud parameters from different cloud retrieval algorithms for the Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) are compared: the Optical Cloud Recognition Algorithm (OCRA) a priori cloud fraction, the Retrieval Of Cloud Information using Neural Networks (ROCINN) CAL (Clouds-As-Layers) cloud fraction and cloud top and base height, the ROCINN CRB (Clouds-as-Reflecting-Boundaries) cloud fraction and cloud height, the Fast Retrieval Scheme for Clouds from the Oxygen A-band (FRESCO) cloud fraction, the interpolated FRESCO cloud height from the TROPOMI NO2 product, the cloud fraction from the NO2 fitting window, the O2–O2 cloud fraction and cloud height, the Mainz Iterative Cloud Retrieval Utilities (MICRU) cloud fraction, and the Visible Infrared Imaging Radiometer Suite (VIIRS) cloud fraction. Two different versions of the TROPOMI cloud products OCRA/ROCINN, FRESCO, and the TROPOMI NO2 product are included in the comparisons (processor version 1.x and 2.x). Overall, the cloud parameters retrieved by the different algorithms show qualitative consistency in version 1.x and good agreement in version 2.x with the exception of the VIIRS cloud fraction, which cannot be directly compared to the other data. Differences between the cloud retrievals are found especially for small cloud heights with a cloud fraction threshold of 0.2, i.e. clouds that are particularly relevant for tropospheric trace gas retrievals. The cloud fractions of the different version 2 cloud products primarily differ over snow- and ice-covered pixels and scenes with sun glint, for which only MICRU includes an explicit treatment. All cloud parameters show some systematic problems related to the across-track dependence, where larger values are found at the edges of the satellite view. The consistency between the cloud parameters from different algorithms depends strongly on how the data are filtered for the comparison, for example, what quality value is used or whether snow- and ice-covered pixels are excluded from the analysis. In summary, clear differences were found between the results of various algorithms, but these differences are reduced in the most recent versions of the cloud data.

Operational retrievals of tropospheric trace gases from space-borne spectrometers are based on one-dimensional radiative transfer models. To minimize cloud effects, trace gas retrievals generally implement a simple cloud model based on radiometric cloud fraction estimates and photon path length corrections. The latter relies on measurements of the oxygen collision pair (O2–O2) absorption at 477 nm or on the oxygen A-band around 760 nm to determine an effective cloud height. In reality however, the impact of clouds is much more complex, involving unresolved sub-pixel clouds, scattering of clouds in neighbouring pixels, and cloud shadow effects, such that unresolved three-dimensional effects due to clouds may introduce significant biases in trace gas retrievals. Although clouds have significant effects on trace gas retrievals, the current cloud correction schemes are based on a simple cloud model, and the retrieved cloud parameters must be interpreted as effective values. Consequently, it is difficult to assess the accuracy of the cloud correction only based on analysis of the accuracy of the cloud retrievals, and this study focuses solely on the impact of the 3D cloud structures on the trace gas retrievals. In order to quantify this impact, we study NO2 as a trace gas example and apply standard retrieval methods including approximate cloud corrections to synthetic data generated by the state-of-the-art three-dimensional Monte Carlo radiative transfer model MYSTIC. A sensitivity study is performed for simulations including a box cloud, and the dependency on various parameters is investigated. The most significant bias is found for cloud shadow effects under polluted conditions. Biases depend strongly on cloud shadow fraction, NO2 profile, cloud optical thickness, solar zenith angle, and surface albedo. Several approaches to correct NO2 retrievals under cloud shadow conditions are explored. We find that air mass factors calculated using fitted surface albedo or corrected using the O2–O2 slant column density can partly mitigate cloud shadow effects. However, these approaches are limited to cloud-free pixels affected by surrounding clouds. A parameterization approach is presented based on relationships derived from the sensitivity study. This allows measurements to be identified for which the standard NO2 retrieval produces a significant bias and therefore provides a way to improve the current data flagging approach.

A TROPOspheric Monitoring Instrument (TROPOMI) validation campaign was held in the Netherlands based at the CESAR (Cabauw Experimental Site for Atmospheric Research) observatory during September 2019. The TROpomi vaLIdation eXperiment (TROLIX-19) consisted of active and passive remote sensing platforms in conjunction with several balloon-borne and surface chemical (e.g., ozone and nitrogen dioxide) measurements. The goal of this joint NASA-KNMI geophysical validation campaign was to make intensive observations in the TROPOMI domain in order to be able to establish the quality of the L2 satellite data products under realistic conditions, such as non-idealized conditions with varying cloud cover and a range of atmospheric conditions at a rural site. The research presented here focuses on using ozone lidars from NASA's Goddard Space Flight Center to better evaluate the characterization of ozone throughout TROLIX-19. Results of comparisons to the lidar systems with balloon, space-borne and ground-based passive measurements are shown. In addition, results are compared to a global coupled chemistry meteorology model to illustrate the vertical variability and columnar amounts of both tropospheric and stratospheric ozone during the campaign period.

Knowledge of sulfur dioxide layer height (SO2 LH) is important to understand volcanic eruption processes, the climate impact of SO2 emissions and to mitigate volcanic risk for civil aviation. However, the estimation of SO2 LH from ground-based instruments is challenging in particular for rapidly evolving and sustained eruptions. Satellite wide-swath nadir observations have the advantage to cover large-scale plumes and the potential to provide key information on SO2 LH. In the ultraviolet, SO2 LH retrievals leverage the fact that, for large SO2 columns, the light path and its associated air mass factor (AMF) depends on the SO2 absorption (and therefore on the vertical distribution of SO2), and SO2 LH information can be obtained from the analysis of measured back-scattered radiances coupled with radiative transfer simulations. However, existing algorithms are mainly sensitive to SO2 LH for SO2 vertical columns of at least 20 DU. Here we develop a new SO2 LH algorithm and apply it to observations from the high-spatial-resolution TROPOspheric Monitoring Instrument (TROPOMI). It is based on an SO2 optical depth look-up table and an iterative approach. The strength of this scheme lies in the fact that it is a Covariance-Based Retrieval Algorithm (COBRA; Theys et al., 2021). This means that the SO2-free contribution of the measured optical depth is treated in an optimal way, resulting in an improvement of the SO2 LH sensitivity to SO2 columns as low as 5 DU, with a precision better than 2 km. We demonstrate the value of this new data through a number of examples and comparison with satellite plume height estimates (from IASI and CALIOP), and back-trajectory analyses. The comparisons indicate an SO2 LH accuracy of 1–2 km, except for some difficult observation conditions, in particular for optically thick ash plumes or partially SO2-filled scenes.

Dual-scan ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements of tropospheric nitrogen dioxide (NO2) and aerosols were carried out in Uccle (50.8∘ N, 4.35∘ E; Brussels region, Belgium) for 2 years from March 2018 to February 2020. The MAX-DOAS instrument operated in both UV and visible wavelength ranges in a dual-scan configuration consisting of two submodes: (1) an elevation scan in a fixed viewing azimuthal direction and (2) an azimuthal scan in a fixed low elevation angle (2∘). By analyzing the O4 and NO2 differential slant column density (dSCD) at six different wavelength intervals along every azimuthal direction and by applying a new optimal-estimation-based inversion approach (the so-called mapping MAX-DOAS technique), the horizontal distribution of the NO2 near-surface concentrations and vertical column densities (VCDs) as well as the aerosol near-surface extinction coefficients are retrieved along 10 azimuthal directions. The retrieved horizontal NO2 concentration profiles allow the identification of the main NO2 hotspots in the Brussels area. Correlative comparisons of the retrieved horizontal NO2 distribution were conducted with airborne, mobile, air quality model, and satellite datasets, and overall good agreement is found. The comparison with TROPOMI observations from operational and scientific data products reveals that the characterization of the horizontal distribution of tropospheric NO2 VCDs by ground-based measurements and an adequate a priori NO2 profile shape in TROPOMI retrievals lead to better consistency between satellite and ground-based datasets.

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.

This paper reports on ground-based validation of the atmospheric OClO data record produced within the framework of EUMETSAT's Satellite Application Facility on Atmospheric Chemistry Monitoring (AC SAF) using the Global Ozone Monitoring Experiment (GOME)-2A and GOME-2B instrument measurements, covering the 2007–2016 and 2013–2016 periods, respectively. OClO slant column densities are compared to correlative measurements collected from nine Zenith-Scattered-Light Differential Optical Absorption Spectroscopy (ZSL-DOAS) instruments from the Network for the Detection of Atmospheric Composition Change (NDACC) distributed in both the Arctic and Antarctic. Sensitivity tests are performed on the ground-based data to estimate the impact of the different OClO DOAS analysis settings. On this basis, we infer systematic uncertainties of about 25 % (i.e., about 3.75×1013 molec. cm−2) between the different ground-based data analyses, reaching total uncertainties ranging from about 26 % to 33 % for the different stations (i.e., around 4 to 5×1013 molec. cm−2). Time series at the different sites show good agreement between satellite and ground-based data for both the inter-annual variability and the overall OClO seasonal behavior. GOME-2A results are found to be noisier than those of GOME-2B, especially after 2011, probably due to instrumental degradation effects. Daily linear regression analysis for OClO-activated periods yield correlation coefficients of 0.8 for GOME-2A and 0.87 for GOME-2B, with slopes with respect to the ground-based data ensemble of 0.64 and 0.72, respectively. Satellite minus ground-based offsets are within 8×1013 molec. cm−2, with some differences between GOME-2A and GOME-2B depending on the station. Overall, considering all the stations, a median offset of about -2.2×1013 molec. cm−2 is found for both GOME-2 instruments.

Three-dimensional (3D) cloud structures may impact atmospheric trace gas products from ultraviolet–visible (UV–Vis) sounders. We used synthetic and observational data to identify and quantify possible cloud-related bias in NO2 tropospheric vertical column density (TVCD). The synthetic data were based on high-resolution large eddy simulations which were input to a 3D radiative transfer model. The simulated visible spectra for low-earth-orbiting and geostationary geometries were analysed with standard retrieval methods and cloud correction schemes that are employed in operational NO2 satellite products. For the observational data, the NO2 products from the TROPOspheric Monitoring Instrument (TROPOMI) were used, while the Visible Infrared Imaging Radiometer Suite (VIIRS) provided high-spatial-resolution cloud and radiance data. NO2 profile shape, cloud shadow fraction, cloud top height, cloud optical depth, and solar zenith and viewing angles were identified as the metrics being the most important in identifying 3D cloud impacts on NO2 TVCD retrievals. For a solar zenith angle less than about 40∘ the synthetic data show that the NO2 TVCD bias is typically below 10 %, while for larger solar zenith angles the NO2 TVCD is low-biased by tens of percent. The horizontal variability of NO2 and differences in TROPOMI and VIIRS overpass times make it challenging to identify a similar bias in the observational data. However, for optically thick clouds above 3000 m, a low bias appears to be present in the observational data.

We present an updated perspective on near-global total ozone trends for the period 1995–2020. We use the GOME-type (Global Ozone Monitoring Experiment) Total Ozone Essential Climate Variable (GTO-ECV) satellite data record which has been extended and generated as part of the European Space Agency's Climate Change Initiative (ESA-CCI) and European Union Copernicus Climate Change Service (EU-C3S) ozone projects. The focus of our work is to examine the regional patterns and seasonal dependency of the ozone trend. In the Southern Hemisphere we found regions that indicate statistically significant positive trends increasing from 0.6 ± 0.5(2σ) % per decade in the subtropics to 1.0 ± 0.9 % per decade in the middle latitudes and 2.8 ± 2.6 % per decade in the latitude band 60–70∘ S. In the middle latitudes of the Northern Hemisphere the trend exhibits distinct regional patterns, i.e., latitudinal and longitudinal structures. Significant positive trends (∼ 1.5 ± 1.0 % per decade) over the North Atlantic region, as well as barely significant negative trends (−1.0 ± 1.0 % per decade) over eastern Europe, were found. Moreover, these trends correlate with long-term changes in tropopause pressure. Total ozone trends in the tropics are not statistically significant. Regarding the seasonal dependence of the trends we found only very small variations over the course of the year. However, we identified different behavior depending on latitude. In the latitude band 40–70∘ N the positive trend maximizes in boreal winter from December to February. In the middle latitudes of the Southern Hemisphere (35–50∘ S) the trend is maximum from March to May. Further south toward the high latitudes (55–70∘ S) the trend exhibits a relatively strong seasonal cycle which varies from 2 % per decade in December and January to 3.8 % per decade in June and July.

In this study we focus on the retrieval of aerosol and trace gas vertical profiles from multi-axis differential optical absorption spectroscopy (MAX-DOAS) observations for the first time over Thessaloniki, Greece. We use two independent inversion algorithms for the profile retrievals: the Mexican MAX-DOAS Fit (MMF) and the Mainz Profile Algorithm (MAPA). The former is based on the optimal estimation method (OEM), while the latter follows a parameterization approach. We evaluate the performance of MMF and MAPA, and we validate their retrieved products with ancillary data measured by other co-located reference instruments. The trace gas differential slant column densities (dSCDs), simulated by the forward models, are in good agreement, except for HCHO, where larger scatter is observed due to the increased spectral noise of the measurements in the UV. We find an excellent agreement between the tropospheric column densities of NO2 retrieved by MMF and MAPA (slope=1.009, Pearson's correlation coefficient R=0.982) and a good correlation for the case of HCHO (R=0.927). For aerosols, we find better agreement for the aerosol optical depths (AODs) in the visible (i.e., at 477 nm) compared to the UV (at 360 nm), and we show that the agreement strongly depends on the O4 scaling factor that is used in the analysis. The agreement for NO2 and HCHO near-surface concentrations is similar to the comparison of the integrated columns with slightly decreased correlation coefficients. The seasonal mean vertical profiles that are retrieved by MMF and MAPA are intercompared, and the seasonal variation in all species along with possible sources is discussed. The AODs retrieved by the MAX-DOAS are validated by comparing them with AOD values measured by a CIMEL sun photometer and a Brewer spectrophotometer. Four different flagging schemes were applied to the data in order to evaluate their performance. Qualitatively, a generally good agreement is observed for both wavelengths, but we find a systematic bias from the CIMEL sun photometer and Brewer spectrophotometer measurements, due to the limited sensitivity of the MAX-DOAS in retrieving information at higher altitudes, especially in the UV. An in-depth validation of the aerosol vertical profiles retrieved by the MAX-DOAS is not possible since only in very few cases is the true aerosol profile known during the period of study. However, we examine four cases, where the MAX-DOAS provided a generally good estimation of the shape of the profiles retrieved by a co-located multi-wavelength lidar system. The NO2 near-surface concentrations are validated against in situ observations, and the comparison of both MMF and MAPA revealed good agreement with correlation coefficients of R=0.78 and R=0.73, respectively. Finally, the effect of the O4 scaling factor is investigated by intercomparing the integrated columns retrieved by the two algorithms and also by comparing the AODs derived by MAPA for different values of the scaling factor with AODs measured by the CIMEL sun photometer and the Brewer spectrophotometer.

Retrievals of trace gas concentrations from satellite observations are mostly performed for clear regions or regions with low cloud coverage. However, even fully clear pixels can be affected by clouds in the vicinity, either by shadowing or by scattering of radiation from clouds in the clear region. Quantifying the error of retrieved trace gas concentrations due to cloud scattering is a difficult task. One possibility is to generate synthetic data by three-dimensional (3D) radiative transfer simulations using realistic 3D atmospheric input data, including 3D cloud structures. Retrieval algorithms may be applied on the synthetic data, and comparison to the known input trace gas concentrations yields the retrieval error due to cloud scattering. In this paper we present a comprehensive synthetic dataset which has been generated using the Monte Carlo radiative transfer model MYSTIC (Monte Carlo code for the phYSically correct Tracing of photons In Cloudy atmospheres). The dataset includes simulated spectra in two spectral ranges (400–500 nm and the O2A-band from 755–775 nm). Moreover it includes layer air mass factors (layer-AMFs) calculated at 460 nm. All simulations are performed for a fixed background atmosphere for various sun positions, viewing directions and surface albedos. Two cloud setups are considered: the first includes simple box clouds with various geometrical and optical thicknesses. This can be used to systematically investigate the sensitivity of the retrieval error on solar zenith angle, surface albedo and cloud parameters. Corresponding 1D simulations are also provided. The second includes realistic three-dimensional clouds from an ICON large eddy simulation (LES) for a region covering Germany and parts of surrounding countries. The scene includes cloud types typical of central Europe such as shallow cumulus, convective cloud cells, cirrus and stratocumulus. This large dataset can be used to quantify the trace gas concentration retrieval error statistically. Along with the dataset, the impact of horizontal photon transport on reflectance spectra and layer-AMFs is analysed for the box-cloud scenarios. Moreover, the impact of 3D cloud scattering on the NO2 vertical column density (VCD) retrieval is presented for a specific LES case. We find that the retrieval error is largest in cloud shadow regions, where the NO2 VCD is underestimated by more than 20 %. The dataset is available for the scientific community to assess the behaviour of trace gas retrieval algorithms and cloud correction schemes in cloud conditions with 3D structure.

Iodine compounds destroy ozone (O3) in the global troposphere and form new aerosols, thereby affecting the global radiative balance. However, few reports have described the latitudinal distribution of atmospheric iodine compounds. This work reports iodine monoxide (IO) measurements taken over unprecedented sampling areas from the Arctic to the Southern Hemisphere and spanning sea surface temperatures (SSTs) of approximately 0 to 31.5 ∘C. The highest IO concentrations were observed over the Western Pacific warm pool (WPWP), where O3 minima were also measured. There, a negative correlation was found between O3 and IO mixing ratios at extremely low O3 concentrations. This correlation is not explained readily by the O3-dependent oceanic fluxes of photolabile inorganic iodine compounds, which is the dominant source in recent global-scale chemistry transport models representing iodine chemistry. Actually, the correlation rather implies that O3-independent pathways can be similarly important in the WPWP. The O3-independent fluxes result in a 15 % greater O3 loss than that estimated for O3-dependent processes alone. The daily O3 loss rate related to iodine over the WPWP is as high as approximately 2 ppbv (parts per billion by volume) despite low O3 concentrations of approximately 10 ppbv, with the loss being up to 100 % greater than that without iodine. This finding suggests that warming SST driven by climate change might affect the marine atmospheric chemical balance through iodine–ozone chemistry.

We present the first global glyoxal (CHOCHO) tropospheric column product derived from the TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor satellite. Atmospheric glyoxal results from the oxidation of other non-methane volatile organic compounds (NMVOCs) and from direct emissions caused by combustion processes. Therefore, this product is a useful indicator of VOC emissions. It is generated with an improved version of the BIRA-IASB scientific retrieval algorithm relying on the differential optical absorption spectroscopy (DOAS) approach. Among the algorithmic updates, the DOAS fit now includes corrections to mitigate the impact of spectral misfits caused by scene brightness inhomogeneity and strong NO2 absorption. The product comes along with a full error characterization, which allows for providing random and systematic error estimates for every observation. Systematic errors are typically in the range of 1 ×1014–3 ×1014 molec. cm−2 (∼30 %–70 % in emission regimes) and originate mostly from a priori data uncertainties and spectral interferences with other absorbing species. The latter may be at the origin, at least partly, of an enhanced glyoxal signal over equatorial oceans, and further investigation is needed to mitigate them. Random errors are large (>6×1014 molec. cm−2) but can be reduced by averaging observations in space and/or time. Benefiting from a high signal-to-noise ratio and a large number of small-size observations, TROPOMI provides glyoxal tropospheric column fields with an unprecedented level of detail. Using the same retrieval algorithmic baseline, glyoxal column data sets are also generated from the Ozone Monitoring Instrument (OMI) on Aura and from the Global Ozone Monitoring Experiment-2 (GOME-2) on board Metop-A and Metop-B. Those four data sets are intercompared over large-scale regions worldwide and show a high level of consistency. The satellite glyoxal columns are also compared to glyoxal columns retrieved from ground-based Multi-AXis DOAS (MAX-DOAS) instruments at nine stations in Asia and Europe. In general, the satellite and MAX-DOAS instruments provide consistent glyoxal columns both in terms of absolute values and variability. Correlation coefficients between TROPOMI and MAX-DOAS glyoxal columns range between 0.61 and 0.87. The correlation is only poorer at one mid-latitude station, where satellite data appear to be biased low during wintertime. The mean absolute glyoxal columns from satellite and MAX-DOAS generally agree well for low/moderate columns with differences of less than 1×1014 molec. cm−2. A larger bias is identified at two sites where the MAX-DOAS columns are very large. Despite this systematic bias, the consistency of the satellite and MAX-DOAS glyoxal seasonal variability is high.

The purpose of the EUNADICS-AV (European Natural Airborne Disaster Information and Coordination System for Aviation) prototype early warning system (EWS) is to develop the combined use of harmonised data products from satellite, ground-based and in situ instruments to produce alerts of airborne hazards (volcanic, dust, smoke and radionuclide clouds), satisfying the requirement of aviation air traffic management (ATM) stakeholders (https://cordis.europa.eu/project/id/723986, last access: 5 November 2021). The alert products developed by the EUNADICS-AV EWS, i.e. near-real-time (NRT) observations, email notifications and netCDF (Network Common Data Form) alert data products (called NCAP files), have shown significant interest in using selective detection of natural airborne hazards from polar-orbiting satellites. The combination of several sensors inside a single global system demonstrates the advantage of using a triggered approach to obtain selective detection from observations, which cannot initially discriminate the different aerosol types. Satellite products from hyperspectral ultraviolet–visible (UV–vis) and infrared (IR) sensors (e.g. TROPOMI – TROPOspheric Monitoring Instrument – and IASI – Infrared Atmospheric Sounding Interferometer) and a broadband geostationary imager (Spinning Enhanced Visible and InfraRed Imager; SEVIRI) and retrievals from ground-based networks (e.g. EARLINET – European Aerosol Research Lidar Network, E-PROFILE and the regional network from volcano observatories) are combined by our system to create tailored alert products (e.g. selective ash detection, SO2 column and plume height, dust cloud, and smoke from wildfires). A total of 23 different alert products are implemented, using 1 geostationary and 13 polar-orbiting satellite platforms, 3 external existing service, and 2 EU and 2 regional ground-based networks. This allows for the identification and the tracking of extreme events. The EUNADICS-AV EWS has also shown the need to implement a future relay of radiological data (gamma dose rate and radionuclides concentrations in ground-level air) in the case of a nuclear accident. This highlights the interest of operating early warnings with the use of a homogenised dataset. For the four types of airborne hazard, the EUNADICS-AV EWS has demonstrated its capability to provide NRT alert data products to trigger data assimilation and dispersion modelling providing forecasts and inverse modelling for source term estimate. Not all of our alert data products (NCAP files) are publicly disseminated. Access to our alert products is currently restricted to key users (i.e. Volcanic Ash Advisory Centres, national meteorological services, the World Meteorological Organization, governments, volcano observatories and research collaborators), as these are considered pre-decisional products. On the other hand, thanks to the EUNADICS-AV–SACS (Support to Aviation Control Service) web interface (https://sacs.aeronomie.be, last access: 5 November 2021), the main part of the satellite observations used by the EUNADICS-AV EWS is shown in NRT, with public email notification of volcanic emission and delivery of tailored images and NCAP files. All of the ATM stakeholders (e.g. pilots, airlines and passengers) can access these alert products through this free channel.

Sensitive and accurate detection of sulfur dioxide (SO2) from space is important for monitoring and estimating global sulfur emissions. Inspired by detection methods applied in the thermal infrared, we present here a new scheme to retrieve SO2 columns from satellite observations of ultraviolet back-scattered radiances. The retrieval is based on a measurement error covariance matrix to fully represent the SO2-free radiance variability, so that the SO2 slant column density is the only retrieved parameter of the algorithm. We demonstrate this approach, named COBRA, on measurements from the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Sentinel-5 Precursor (S-5P) satellite. We show that the method reduces significantly both the noise and biases present in the current TROPOMI operational DOAS SO2 retrievals. The performance of this technique is also benchmarked against that of the principal component algorithm (PCA) approach. We find that the quality of the data is similar and even slightly better with the proposed COBRA approach. The ability of the algorithm to retrieve SO2 accurately is further supported by comparison with ground-based observations. We illustrate the great sensitivity of the method with a high-resolution global SO2 map, considering 2.5 years of TROPOMI data. In addition to the known sources, we detect many new SO2 emission hotspots worldwide. For the largest sources, we use the COBRA data to estimate SO2 emission rates. Results are comparable to other recently published TROPOMI-based SO2 emissions estimates, but the associated uncertainties are significantly lower than with the operational data. Next, for a limited number of weak sources, we demonstrate the potential of our data for quantifying SO2 emissions with a detection limit of about 8 kt yr−1, a factor of 4 better than the emissions derived from the Ozone Monitoring Instrument (OMI). We anticipate that the systematic use of our TROPOMI COBRA SO2 column data set at a global scale will allow missing sources to be identified and quantified and help improve SO2 emission inventories.

Launched in October 2017, the TROPOspheric Monitoring Instrument (TROPOMI) aboard Sentinel-5 Precursor provides the potential to monitor air quality over point sources across the globe with a spatial resolution as high as 5.5 km × 3.5 km (7 km × 3.5 km before 6 August 2019). The DLR nitrogen dioxide (NO2) retrieval algorithm for the TROPOMI instrument consists of three steps: the spectral fitting of the slant column, the separation of stratospheric and tropospheric contributions, and the conversion of the slant column to a vertical column using an air mass factor (AMF) calculation. In this work, an improved DLR tropospheric NO2 retrieval algorithm from TROPOMI measurements over Europe is presented. The stratospheric estimation is implemented using the STRatospheric Estimation Algorithm from Mainz (STREAM), which was developed as a verification algorithm for TROPOMI and does not require chemistry transport model data as input. A directionally dependent STREAM (DSTREAM) is developed to correct for the dependency of the stratospheric NO2 on the viewing geometry by up to 2×1014 molec./cm2. Applied to synthetic TROPOMI data, the uncertainty in the stratospheric column is 3.5×1014 molec./cm2 in the case of significant tropospheric sources. Applied to actual measurements, the smooth variation of stratospheric NO2 at low latitudes is conserved, and stronger stratospheric variation at higher latitudes is captured. For AMF calculation, the climatological surface albedo data are replaced by geometry-dependent effective Lambertian equivalent reflectivity (GE_LER) obtained directly from TROPOMI measurements with a high spatial resolution. Mesoscale-resolution a priori NO2 profiles are obtained from the regional POLYPHEMUS/DLR chemistry transport model with the TNO-MACC emission inventory. Based on the latest TROPOMI operational cloud parameters, a more realistic cloud treatment is provided by a Clouds-As-Layers (CAL) model, which treats the clouds as uniform layers of water droplets, instead of the Clouds-As-Reflecting-Boundaries (CRB) model, in which clouds are simplified as Lambertian reflectors. For the error analysis, the tropospheric AMF uncertainty, which is the largest source of NO2 uncertainty for polluted scenarios, ranges between 20 % and 50 %, leading to a total uncertainty in the tropospheric NO2 column in the 30 %–60 % range. From a validation performed with ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements, the new DLR tropospheric NO2 data show good correlations for nine European urban/suburban stations, with an average correlation coefficient of 0.78. The implementation of the algorithm improvements leads to a decrease of the relative difference from −55.3 % to −34.7 % on average in comparison with the DLR reference retrieval. When the satellite averaging kernels are used to remove the contribution of a priori profile shape, the relative difference decreases further to ∼ −20 %.

The TROPOspheric Monitoring Instrument (TROPOMI), launched in October 2017 on board the Sentinel-5 Precursor (S5P) satellite, monitors the composition of the Earth's atmosphere at an unprecedented horizontal resolution as fine as 3.5 × 5.5 km2. This paper assesses the performances of the TROPOMI formaldehyde (HCHO) operational product compared to its predecessor, the OMI (Ozone Monitoring Instrument) HCHO QA4ECV product, at different spatial and temporal scales. The parallel development of the two algorithms favoured the consistency of the products, which facilitates the production of long-term combined time series. The main difference between the two satellite products is related to the use of different cloud algorithms, leading to a positive bias of OMI compared to TROPOMI of up to 30 % in tropical regions. We show that after switching off the explicit correction for cloud effects, the two datasets come into an excellent agreement. For medium to large HCHO vertical columns (larger than 5 × 1015 molec. cm−2) the median bias between OMI and TROPOMI HCHO columns is not larger than 10 % (< 0.4 × 1015 molec. cm−2). For lower columns, OMI observations present a remaining positive bias of about 20 % (< 0.8 × 1015 molec. cm−2) compared to TROPOMI in midlatitude regions. Here, we also use a global network of 18 MAX-DOAS (multi-axis differential optical absorption spectroscopy) instruments to validate both satellite sensors for a large range of HCHO columns. This work complements the study by Vigouroux et al. (2020), where a global FTIR (Fourier transform infrared) network is used to validate the TROPOMI HCHO operational product. Consistent with the FTIR validation study, we find that for elevated HCHO columns, TROPOMI data are systematically low (−25 % for HCHO columns larger than 8 × 1015 molec. cm−2), while no significant bias is found for medium-range column values. We further show that OMI and TROPOMI data present equivalent biases for large HCHO levels. However, TROPOMI significantly improves the precision of the HCHO observations at short temporal scales and for low HCHO columns. We show that compared to OMI, the precision of the TROPOMI HCHO columns is improved by 25 % for individual pixels and by up to a factor of 3 when considering daily averages in 20 km radius circles. The validation precision obtained with daily TROPOMI observations is comparable to the one obtained with monthly OMI observations. To illustrate the improved performances of TROPOMI in capturing weak HCHO signals, we present clear detection of HCHO column enhancements related to shipping emissions in the Indian Ocean. This is achieved by averaging data over a much shorter period (3 months) than required with previous sensors (5 years) and opens new perspectives to study shipping emissions of VOCs (volatile organic compounds) and related atmospheric chemical interactions.

In this paper, we present the MErged GRIdded Dataset of Ozone Profiles (MEGRIDOP) in the stratosphere with a resolved longitudinal structure, which is derived from data from six limb and occultation satellite instruments: GOMOS, SCIAMACHY and MIPAS on Envisat, OSIRIS on Odin, OMPS on Suomi-NPP, and MLS on Aura. The merged dataset was generated as a contribution to the European Space Agency Climate Change Initiative Ozone project (Ozone_cci). The period of this merged time series of ozone profiles is from late 2001 until the end of 2018. The monthly mean gridded ozone profile dataset is provided in the altitude range from 10 to 50 km in bins of 10∘ latitude × 20∘ longitude. The merging is performed using deseasonalized anomalies. The created MEGRIDOP dataset can be used for analyses that probe our understanding of stratospheric chemistry and dynamics. To illustrate some possible applications, we created a climatology of ozone profiles with resolved longitudinal structure. We found zonal asymmetry in the climatological ozone profiles at middle and high latitudes associated with the polar vortex. At northern high latitudes, the amplitude of the seasonal cycle also has a longitudinal dependence. The MEGRIDOP dataset has also been used to evaluate regional vertically resolved ozone trends in the stratosphere, including the polar regions. It is found that stratospheric ozone trends exhibit longitudinal structures at Northern Hemisphere middle and high latitudes, with enhanced trends over Scandinavia and the Atlantic region. This agrees well with previous analyses and might be due to changes in dynamical processes related to the Brewer–Dobson circulation.

The second Cabauw Intercomparison of Nitrogen Dioxide measuring Instruments (CINDI-2) took place in Cabauw (the Netherlands) in September 2016 with the aim of assessing the consistency of multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements of tropospheric species (NO2, HCHO, O3, HONO, CHOCHO and O4). This was achieved through the coordinated operation of 36 spectrometers operated by 24 groups from all over the world, together with a wide range of supporting reference observations (in situ analysers, balloon sondes, lidars, long-path DOAS, direct-sun DOAS, Sun photometer and meteorological instruments). In the presented study, the retrieved CINDI-2 MAX-DOAS trace gas (NO2, HCHO) and aerosol vertical profiles of 15 participating groups using different inversion algorithms are compared and validated against the colocated supporting observations, with the focus on aerosol optical thicknesses (AOTs), trace gas vertical column densities (VCDs) and trace gas surface concentrations. The algorithms are based on three different techniques: six use the optimal estimation method, two use a parameterized approach and one algorithm relies on simplified radiative transport assumptions and analytical calculations. To assess the agreement among the inversion algorithms independent of inconsistencies in the trace gas slant column density acquisition, participants applied their inversion to a common set of slant columns. Further, important settings like the retrieval grid, profiles of O3, temperature and pressure as well as aerosol optical properties and a priori assumptions (for optimal estimation algorithms) have been prescribed to reduce possible sources of discrepancies. The profiling results were found to be in good qualitative agreement: most participants obtained the same features in the retrieved vertical trace gas and aerosol distributions; however, these are sometimes at different altitudes and of different magnitudes. Under clear-sky conditions, the root-mean-square differences (RMSDs) among the results of individual participants are in the range of 0.01–0.1 for AOTs, (1.5–15) ×1014molec.cm-2 for trace gas (NO2, HCHO) VCDs and (0.3–8)×1010molec.cm-3 for trace gas surface concentrations. These values compare to approximate average optical thicknesses of 0.3, trace gas vertical columns of 90×1014molec.cm-2 and trace gas surface concentrations of 11×1010molec.cm-3 observed over the campaign period. The discrepancies originate from differences in the applied techniques, the exact implementation of the algorithms and the user-defined settings that were not prescribed. For the comparison against supporting observations, the RMSDs increase to a range of 0.02–0.2 against AOTs from the Sun photometer, (11–55)×1014molec.cm-2 against trace gas VCDs from direct-sun DOAS observations and (0.8–9)×1010molec.cm-3 against surface concentrations from the long-path DOAS instrument. This increase in RMSDs is most likely caused by uncertainties in the supporting data, spatiotemporal mismatch among the observations and simplified assumptions particularly on aerosol optical properties made for the MAX-DOAS retrieval. As a side investigation, the comparison was repeated with the participants retrieving profiles from their own differential slant column densities (dSCDs) acquired during the campaign. In this case, the consistency among the participants degrades by about 30 % for AOTs, by 180 % (40 %) for HCHO (NO2) VCDs and by 90 % (20 %) for HCHO (NO2) surface concentrations. In former publications and also during this comparison study, it was found that MAX-DOAS vertically integrated aerosol extinction coefficient profiles systematically underestimate the AOT observed by the Sun photometer. For the first time, it is quantitatively shown that for optimal estimation algorithms this can be largely explained and compensated by considering biases arising from the reduced sensitivity of MAX-DOAS observations to higher altitudes and associated a priori assumptions.

Ozone data derived from the Tropospheric Monitoring Instrument (TROPOMI) sensor on board the Sentinel-5 Precursor satellite show exceptionally low total ozone columns in the polar region of the Northern Hemisphere (Arctic) in spring 2020. Minimum total ozone column values around or below 220 Dobson units (DU) were seen over the Arctic for 5 weeks in March and early April 2020. Usually the persistence of such low total ozone column values in spring is only observed in the polar Southern Hemisphere (Antarctic) and not over the Arctic. These record low total ozone columns were caused by a particularly strong polar vortex in the stratosphere with a persistent cold stratosphere at higher latitudes, a prerequisite for ozone depletion through heterogeneous chemistry. Based on the ERA5, which is the fifth generation of the European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis, the Northern Hemisphere winter 2019/2020 (from December to March) showed minimum polar cap temperatures consistently below 195 K around 20 km altitude, which enabled enhanced formation of polar stratospheric clouds. The special situation in spring 2020 is compared and discussed in context with two other Northern Hemisphere spring seasons, namely those in 1997 and 2011, which also displayed relatively low total ozone column values. However, during these years, total ozone columns below 220 DU over several consecutive days were not observed in spring. The similarities and differences of the atmospheric conditions of these three events and possible explanations for the observed features are presented and discussed. It becomes apparent that the monthly mean of the minimum total ozone column value for March 2020 (221 DU) was clearly below the respective values found in March 1997 (267 DU) and 2011 (252 DU), which highlights the special evolution of the polar stratospheric ozone layer in the Northern Hemisphere in spring 2020. A comparison with a typical ozone hole over the Antarctic (e.g., in 2016) indicates that although the Arctic spring 2020 situation is remarkable, with total ozone column values around or below 220 DU observed over a considerable area (up to 0.9 million km2), the Antarctic ozone hole shows total ozone columns typically below 150 DU over a much larger area (of the order of 20 million km2). Furthermore, total ozone columns below 220 DU are typically observed over the Antarctic for about 4 months.

This paper reports on consolidated ground-based validation results of the atmospheric NO2 data produced operationally since April 2018 by the TROPOspheric Monitoring Instrument (TROPOMI) on board of the ESA/EU Copernicus Sentinel-5 Precursor (S5P) satellite. Tropospheric, stratospheric, and total NO2 column data from S5P are compared to correlative measurements collected from, respectively, 19 Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS), 26 Network for the Detection of Atmospheric Composition Change (NDACC) Zenith-Scattered-Light DOAS (ZSL-DOAS), and 25 Pandonia Global Network (PGN)/Pandora instruments distributed globally. The validation methodology gives special care to minimizing mismatch errors due to imperfect spatio-temporal co-location of the satellite and correlative data, e.g. by using tailored observation operators to account for differences in smoothing and in sampling of atmospheric structures and variability and photochemical modelling to reduce diurnal cycle effects. Compared to the ground-based measurements, S5P data show, on average, (i) a negative bias for the tropospheric column data, of typically −23 % to −37 % in clean to slightly polluted conditions but reaching values as high as −51 % over highly polluted areas; (ii) a slight negative median difference for the stratospheric column data, of about −0.2 Pmolec cm−2, i.e. approx. −2 % in summer to −15 % in winter; and (iii) a bias ranging from zero to −50 % for the total column data, found to depend on the amplitude of the total NO2 column, with small to slightly positive bias values for columns below 6 Pmolec cm−2 and negative values above. The dispersion between S5P and correlative measurements contains mostly random components, which remain within mission requirements for the stratospheric column data (0.5 Pmolec cm−2) but exceed those for the tropospheric column data (0.7 Pmolec cm−2). While a part of the biases and dispersion may be due to representativeness differences such as different area averaging and measurement times, it is known that errors in the S5P tropospheric columns exist due to shortcomings in the (horizontally coarse) a priori profile representation in the TM5-MP chemical transport model used in the S5P retrieval and, to a lesser extent, to the treatment of cloud effects and aerosols. Although considerable differences (up to 2 Pmolec cm−2 and more) are observed at single ground-pixel level, the near-real-time (NRTI) and offline (OFFL) versions of the S5P NO2 operational data processor provide similar NO2 column values and validation results when globally averaged, with the NRTI values being on average 0.79 % larger than the OFFL values.

Sentinel-5 Precursor (S-5P), launched in October 2017, carrying the TROPOspheric Monitoring Instrument (TROPOMI) nadir-viewing spectrometer, is the first mission of the Copernicus Programme dedicated to the monitoring of air quality, climate, and ozone. In the presented study, the TROPOMI tropospheric nitrogen dioxide (NO2) level-2 (L2) product (OFFL v1.03.01; 3.5 km × 7 km at nadir observations) has been validated over strongly polluted urban regions by comparison with coincident high-resolution Airborne Prism EXperiment (APEX) remote sensing observations (∼ 75 m × 120 m). Satellite products can be optimally assessed based on (APEX) airborne remote sensing observations, as a large amount of satellite pixels can be fully mapped at high accuracy and in a relatively short time interval, reducing the impact of spatiotemporal mismatches. In the framework of the S-5P validation campaign over Belgium (S5PVAL-BE), the APEX imaging spectrometer has been deployed during four mapping flights (26–29 June 2019) over the two largest urban regions in Belgium, i.e. Brussels and Antwerp, in order to map the horizontal distribution of tropospheric NO2. For each flight, 10 to 20 TROPOMI pixels were fully covered by approximately 2700 to 4000 APEX measurements within each TROPOMI pixel. The TROPOMI and APEX NO2 vertical column density (VCD) retrieval schemes are similar in concept. Overall, for the ensemble of the four flights, the standard TROPOMI NO2 VCD product is well correlated (R = 0.92) but biased negatively by −1.2 ± 1.2 × 1015 molec cm−2 or −14 ± 12 %, on average, with respect to coincident APEX NO2 retrievals. When replacing the coarse 1∘ × 1∘ the massively parallel (MP) version of the Tracer Model version 5 (TM5) a priori NO2 profiles by NO2 profile shapes from the Copernicus Atmospheric Monitoring Service (CAMS) regional chemistry transport model (CTM) ensemble at 0.1∘ × 0.1∘, R is 0.94 and the slope increases from 0.82 to 0.93. The bias is reduced to −0.1 ± 1.0 × 1015 molec cm−2 or −1.0 ± 12 %. The absolute difference is on average 1.3 × 1015 molec cm−2 (16 %) and 0.7 × 1015 molec cm−2 (9 %), when comparing APEX NO2 VCDs with TM5-MP-based and CAMS-based NO2 VCDs, respectively. Both sets of retrievals are well within the mission accuracy requirement of a maximum bias of 25 %–50 % for the TROPOMI tropospheric NO2 product for all individual compared pixels. Additionally, the APEX data set allows the study of TROPOMI subpixel variability and impact of signal smoothing due to its finite satellite pixel size, typically coarser than fine-scale gradients in the urban NO2 field. For a case study in the Antwerp region, the current TROPOMI data underestimate localized enhancements and overestimate background values by approximately 1–2 × 1015 molec cm−2 (10 %–20 %).

Multi-axis differential optical absorption spectroscopy (MAX-DOAS) and direct sun NO2 vertical column network data are used to investigate the accuracy of tropospheric NO2 column measurements of the GOME-2 instrument on the MetOp-A satellite platform and the OMI instrument on Aura. The study is based on 23 MAX-DOAS and 16 direct sun instruments at stations distributed worldwide. A method to quantify and correct for horizontal dilution effects in heterogeneous NO2 field conditions is proposed. After systematic application of this correction to urban sites, satellite measurements are found to present smaller biases compared to ground-based reference data in almost all cases. We investigate the seasonal dependence of the validation results as well as the impact of using different approaches to select satellite ground pixels in coincidence with ground-based data. In optimal comparison conditions (satellite pixels containing the station) the median bias between satellite tropospheric NO2 column measurements and the ensemble of MAX-DOAS and direct sun measurements is found to be significant and equal to −34 % for GOME-2A and −24 % for OMI. These biases are further reduced to −24 % and −18 % respectively, after application of the dilution correction. Comparisons with the QA4ECV satellite product for both GOME-2A and OMI are also performed, showing less scatter but also a slightly larger median tropospheric NO2 column bias with respect to the ensemble of MAX-DOAS and direct sun measurements.

Ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements of aerosols and tropospheric nitrogen dioxide (NO2) were carried out in Uccle (50.8∘ N, 4.35∘ E), Brussels, during 1 year from March 2018 until March 2019. The instrument was operated in both the UV and visible wavelength ranges in a dual-scan configuration consisting of two sub-modes: (1) an elevation scan in a fixed viewing azimuthal direction (the so-called main azimuthal direction) pointing to the northeast and (2) an azimuthal scan in a fixed low elevation angle (2∘). By applying a vertical profile inversion algorithm in the main azimuthal direction and a parameterization technique in the other azimuthal directions, near-surface NO2 volume mixing ratios (VMRs) and vertical column densities (VCDs) were retrieved in 10 different azimuthal directions. The dual-scan MAX-DOAS dataset allows for partly resolving the horizontal distribution of NO2 around the measurement site and studying its seasonal variations. Furthermore, we show that measuring the tropospheric NO2 VCDs in different azimuthal directions improves the spatial colocation with measurements from the Sentinel-5 Precursor (S5P), leading to a reduction of the spread in validation results. By using NO2 vertical profile information derived from the MAX-DOAS measurements, we also resolve a systematic underestimation in S5P NO2 data due to the use of inadequate a priori NO2 profile shape data in the satellite retrieval.

The Airborne ROmanian Measurements of Aerosols and Trace gases (AROMAT) campaigns took place in Romania in September 2014 and August 2015. They focused on two sites: the Bucharest urban area and large power plants in the Jiu Valley. The main objectives of the campaigns were to test recently developed airborne observation systems dedicated to air quality studies and to verify their applicability for the validation of space-borne atmospheric missions such as the TROPOspheric Monitoring Instrument (TROPOMI)/Sentinel-5 Precursor (S5P). We present the AROMAT campaigns from the perspective of findings related to the validation of tropospheric NO2, SO2, and H2CO. We also quantify the emissions of NOx and SO2 at both measurement sites. We show that tropospheric NO2 vertical column density (VCD) measurements using airborne mapping instruments are well suited for satellite validation in principle. The signal-to-noise ratio of the airborne NO2 measurements is an order of magnitude higher than its space-borne counterpart when the airborne measurements are averaged at the TROPOMI pixel scale. However, we show that the temporal variation of the NO2 VCDs during a flight might be a significant source of comparison error. Considering the random error of the TROPOMI tropospheric NO2 VCD (σ), the dynamic range of the NO2 VCDs field extends from detection limit up to 37 σ (2.6×1016 molec. cm−2) and 29 σ (2×1016 molec. cm−2) for Bucharest and the Jiu Valley, respectively. For both areas, we simulate validation exercises applied to the TROPOMI tropospheric NO2 product. These simulations indicate that a comparison error budget closely matching the TROPOMI optimal target accuracy of 25 % can be obtained by adding NO2 and aerosol profile information to the airborne mapping observations, which constrains the investigated accuracy to within 28 %. In addition to NO2, our study also addresses the measurements of SO2 emissions from power plants in the Jiu Valley and an urban hotspot of H2CO in the centre of Bucharest. For these two species, we conclude that the best validation strategy would consist of deploying ground-based measurement systems at well-identified locations.

We present the inter-comparison of delta slant column densities (SCDs) and vertical profiles of nitrous acid (HONO) derived from measurements of different multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments and using different inversion algorithms during the Second Cabauw Inter-comparison campaign for Nitrogen Dioxide measuring Instruments (CINDI-2) in September 2016 at Cabauw, the Netherlands (51.97∘ N, 4.93∘ E). The HONO vertical profiles, vertical column densities (VCDs), and near-surface volume mixing ratios are compared between different MAX-DOAS instruments and profile inversion algorithms for the first time. Systematic and random discrepancies of the HONO results are derived from the comparisons of all data sets against their median values. Systematic discrepancies of HONO delta SCDs are observed in the range of ±0.3×1015 molec. cm−2, which is half of the typical random discrepancy of 0.6×1015 molec. cm−2. For a typical high HONO delta SCD of 2×1015 molec. cm−2, the relative systematic and random discrepancies are about 15 % and 30 %, respectively. The inter-comparison of HONO profiles shows that both systematic and random discrepancies of HONO VCDs and near-surface volume mixing ratios (VMRs) are mostly in the range of ∼±0.5×1014 molec. cm−2 and ∼±0.1 ppb (typically ∼20 %). Further we find that the discrepancies of the retrieved HONO profiles are dominated by discrepancies of the HONO delta SCDs. The profile retrievals only contribute to the discrepancies of the HONO profiles by ∼5 %. However, some data sets with substantially larger discrepancies than the typical values indicate that inappropriate implementations of profile inversion algorithms and configurations of radiative transfer models in the profile retrievals can also be an important uncertainty source. In addition, estimations of measurement uncertainties of HONO dSCDs, which can significantly impact profile retrievals using the optimal estimation method, need to consider not only DOAS fit errors, but also atmospheric variability, especially for an instrument with a DOAS fit error lower than ∼3×1014 molec. cm−2. The MAX-DOAS results during the CINDI-2 campaign indicate that the peak HONO levels (e.g. near-surface VMRs of ∼0.4 ppb) often appeared in the early morning and below 0.2 km. The near-surface VMRs retrieved from the MAX-DOAS observations are compared with those measured using a co-located long-path DOAS instrument. The systematic differences are smaller than 0.15 and 0.07 ppb during early morning and around noon, respectively. Since true HONO values at high altitudes are not known in the absence of real measurements, in order to evaluate the abilities of profile inversion algorithms to respond to different HONO profile shapes, we performed sensitivity studies using synthetic HONO delta SCDs simulated by a radiative transfer model with assumed HONO profiles. The tests indicate that the profile inversion algorithms based on the optimal estimation method with proper configurations can reproduce the different HONO profile shapes well. Therefore we conclude that the features of HONO accumulated near the surface derived from MAX-DOAS measurements are expected to represent the ambient HONO profiles well.

The QA4ECV (Quality Assurance for Essential Climate Variables) version 1.1 stratospheric and tropospheric NO2 vertical column density (VCD) climate data records (CDRs) from the OMI (Ozone Monitoring Instrument) satellite sensor are validated using NDACC (Network for the Detection of Atmospheric Composition Change) zenith-scattered light differential optical absorption spectroscopy (ZSL-DOAS) and multi-axis DOAS (MAX-DOAS) data as a reference. The QA4ECV OMI stratospheric VCDs have a small bias of ∼0.2 Pmolec.cm-2 (5 %–10 %) and a dispersion of 0.2 to 1 Pmolec.cm-2 with respect to the ZSL-DOAS measurements. QA4ECV tropospheric VCD observations from OMI are restricted to near-cloud-free scenes, leading to a negative sampling bias (with respect to the unrestricted scene ensemble) of a few peta molecules per square centimetre (Pmolec.cm-2) up to −10 Pmolec.cm-2 (−40 %) in one extreme high-pollution case. The QA4ECV OMI tropospheric VCD has a negative bias with respect to the MAX-DOAS data (−1 to −4 Pmolec.cm-2), which is a feature also found for the OMI OMNO2 standard data product. The tropospheric VCD discrepancies between satellite measurements and ground-based data greatly exceed the combined measurement uncertainties. Depending on the site, part of the discrepancy can be attributed to a combination of comparison errors (notably horizontal smoothing difference error), measurement/retrieval errors related to clouds and aerosols, and the difference in vertical smoothing and a priori profile assumptions.

TROPOMI (the TROPOspheric Monitoring Instrument), on board the Sentinel-5 Precursor (S5P) satellite, has been monitoring the Earth's atmosphere since October 2017 with an unprecedented horizontal resolution (initially 7 km2×3.5 km2, upgraded to 5.5 km2×3.5 km2 in August 2019). Monitoring air quality is one of the main objectives of TROPOMI; it obtains measurements of important pollutants such as nitrogen dioxide, carbon monoxide, and formaldehyde (HCHO). In this paper we assess the quality of the latest HCHO TROPOMI products versions 1.1.(5-7), using ground-based solar-absorption FTIR (Fourier-transform infrared) measurements of HCHO from 25 stations around the world, including high-, mid-, and low-latitude sites. Most of these stations are part of the Network for the Detection of Atmospheric Composition Change (NDACC), and they provide a wide range of observation conditions, from very clean remote sites to those with high HCHO levels from anthropogenic or biogenic emissions. The ground-based HCHO retrieval settings have been optimized and harmonized at all the stations, ensuring a consistent validation among the sites. In this validation work, we first assess the accuracy of TROPOMI HCHO tropospheric columns using the median of the relative differences between TROPOMI and FTIR ground-based data (BIAS). The pre-launch accuracy requirements of TROPOMI HCHO are 40 %–80 %. We observe that these requirements are well reached, with the BIAS found below 80 % at all the sites and below 40 % at 20 of the 25 sites. The provided TROPOMI systematic uncertainties are well in agreement with the observed biases at most of the stations except for the highest-HCHO-level site, where it is found to be underestimated. We find that while the BIAS has no latitudinal dependence, it is dependent on the HCHO concentration levels: an overestimation (+26±5 %) of TROPOMI is observed for very low HCHO levels (<2.5×1015 molec. cm−2), while an underestimation (-30.8%±1.4 %) is found for high HCHO levels (>8.0×1015 molec. cm−2). This demonstrates the great value of such a harmonized network covering a wide range of concentration levels, the sites with high HCHO concentrations being crucial for the determination of the satellite bias in the regions of emissions and the clean sites allowing a small TROPOMI offset to be determined. The wide range of sampled HCHO levels within the network allows the robust determination of the significant constant and proportional TROPOMI HCHO biases (TROPOMI =+1.10±0.05 ×1015+0.64±0.03 × FTIR; in molecules per square centimetre). Second, the precision of TROPOMI HCHO data is estimated by the median absolute deviation (MAD) of the relative differences between TROPOMI and FTIR ground-based data. The clean sites are especially useful for minimizing a possible additional collocation error. The precision requirement of 1.2×1016 molec. cm−2 for a single pixel is reached at most of the clean sites, where it is found that the TROPOMI precision can even be 2 times better (0.5–0.8×1015 molec. cm−2 for a single pixel). However, we find that the provided TROPOMI random uncertainties may be underestimated by a factor of 1.6 (for clean sites) to 2.3 (for high HCHO levels). The correlation is very good between TROPOMI and FTIR data (R=0.88 for 3 h mean coincidences; R=0.91 for monthly means coincidences). Using about 17 months of data (from May 2018 to September 2019), we show that the TROPOMI seasonal variability is in very good agreement at all of the FTIR sites. The FTIR network demonstrates the very good quality of the TROPOMI HCHO products, which is well within the pre-launch requirements for both accuracy and precision. This paper makes suggestions for the refinement of the TROPOMI random uncertainty budget and TROPOMI quality assurance values for a better filtering of the remaining outliers.

In September 2016, 36 spectrometers from 24 institutes measured a number of key atmospheric pollutants for a period of 17 d during the Second Cabauw Intercomparison campaign for Nitrogen Dioxide measuring Instruments (CINDI-2) that took place at Cabauw, the Netherlands (51.97∘ N, 4.93∘ E). We report on the outcome of the formal semi-blind intercomparison exercise, which was held under the umbrella of the Network for the Detection of Atmospheric Composition Change (NDACC) and the European Space Agency (ESA). The three major goals of CINDI-2 were (1) to characterise and better understand the differences between a large number of multi-axis differential optical absorption spectroscopy (MAX-DOAS) and zenith-sky DOAS instruments and analysis methods, (2) to define a robust methodology for performance assessment of all participating instruments, and (3) to contribute to a harmonisation of the measurement settings and retrieval methods. This, in turn, creates the capability to produce consistent high-quality ground-based data sets, which are an essential requirement to generate reliable long-term measurement time series suitable for trend analysis and satellite data validation. The data products investigated during the semi-blind intercomparison are slant columns of nitrogen dioxide (NO2), the oxygen collision complex (O4) and ozone (O3) measured in the UV and visible wavelength region, formaldehyde (HCHO) in the UV spectral region, and NO2 in an additional (smaller) wavelength range in the visible region. The campaign design and implementation processes are discussed in detail including the measurement protocol, calibration procedures and slant column retrieval settings. Strong emphasis was put on the careful alignment and synchronisation of the measurement systems, resulting in a unique set of measurements made under highly comparable air mass conditions. The CINDI-2 data sets were investigated using a regression analysis of the slant columns measured by each instrument and for each of the target data products. The slope and intercept of the regression analysis respectively quantify the mean systematic bias and offset of the individual data sets against the selected reference (which is obtained from the median of either all data sets or a subset), and the rms error provides an estimate of the measurement noise or dispersion. These three criteria are examined and for each of the parameters and each of the data products, performance thresholds are set and applied to all the measurements. The approach presented here has been developed based on heritage from previous intercomparison exercises. It introduces a quantitative assessment of the consistency between all the participating instruments for the MAX-DOAS and zenith-sky DOAS techniques.

Multi-axis differential optical absorption spectroscopy (MAX-DOAS) tropospheric NO2 column retrievals from four European measurement stations are compared to simulations from five regional air quality models which contribute to the European regional ensemble forecasts and reanalyses of the operational Copernicus Atmosphere Monitoring Service (CAMS). Compared to other observational data usually applied for regional model evaluation, MAX-DOAS data are closer to the regional model data in terms of horizontal and vertical resolution, and multiple measurements are available during daylight, so that, for example, diurnal cycles of trace gases can be investigated. In general, there is good agreement between simulated and retrieved NO2 column values for individual MAX-DOAS measurements with correlations between 35 % and 70 % for individual models and 45 % to 75 % for the ensemble median for tropospheric NO2 vertical column densities (VCDs), indicating that emissions, transport and tropospheric chemistry of NOx are on average well simulated. However, large differences are found for individual pollution plumes observed by MAX-DOAS. Most of the models overestimate seasonal cycles for the majority of MAX-DOAS sites investigated. At the urban stations, weekly cycles are reproduced well, but the decrease towards the weekend is underestimated and diurnal cycles are overall not well represented. In particular, simulated morning rush hour peaks are not confirmed by MAX-DOAS retrievals, and models fail to reproduce observed changes in diurnal cycles for weekdays versus weekends. The results of this study show that future model development needs to concentrate on improving representation of diurnal cycles and associated temporal scalings.

We present different methods for in-field elevation calibration of MAX-DOAS (Multi AXis Differential Optical Absorption Spectroscopy) instruments that were applied and inter-compared during the second Cabauw Intercomparison campaign for Nitrogen Dioxide measuring Instruments (CINDI-2). One necessary prerequisite of consistent MAX-DOAS retrievals is a precise and accurate calibration of the elevation angles of the different measuring systems. Therefore, different methods for this calibration were applied to several instruments during the campaign, and the results were inter-compared. This work first introduces and explains the different methods, namely far- and near-lamp measurements, white-stripe scans, horizon scans and sun scans, using data and results for only one (mainly the Max Planck Institute for Chemistry) instrument. In the second part, the far-lamp measurements and the horizon scans are examined for all participating groups. Here, the results for both methods are first inter-compared for the different instruments; secondly, the two methods are compared amongst each other. All methods turned out to be well-suited for the calibration of the elevation angles of MAX-DOAS systems, with each of them having individual advantages and drawbacks. Considering the results of this study, the systematic uncertainties of the methods can be estimated as ±0.05∘ for the far-lamp measurements and the sun scans, ±0.25∘ for the horizon scans, and around ±0.1∘ for the white-stripe and near-lamp measurements. When comparing the results of far-lamp and horizon-scan measurements, a spread of around 0.9∘ in the elevation calibrations is found between the participating instruments for both methods. This spread is of the order of a typical field of view (FOV) of a MAX-DOAS instrument and therefore affecting the retrieval results. Further, consistent (wavelength dependent) offsets of 0.32∘ and 0.40∘ between far-lamp measurements and horizon scans are found, which can be explained by the fact that, despite the flat topography around the measurement site, obstacles such as trees might mark the visible horizon during daytime. The observed wavelength dependence can be explained by surface albedo effects. Lastly, the results are discussed and recommendations for future campaigns are given.

An improved tropospheric nitrogen dioxide (NO2) retrieval algorithm from the Global Ozone Monitoring Experiment-2 (GOME-2) instrument based on air mass factor (AMF) calculations performed with more realistic model parameters is presented. The viewing angle dependency of surface albedo is taken into account by improving the GOME-2 Lambertian-equivalent reflectivity (LER) climatology with a directionally dependent LER (DLER) dataset over land and an ocean surface albedo parameterisation over water. A priori NO2 profiles with higher spatial and temporal resolutions are obtained from the IFS (CB05BASCOE) chemistry transport model based on recent emission inventories. A more realistic cloud treatment is provided by a clouds-as-layers (CAL) approach, which treats the clouds as uniform layers of water droplets, instead of the current clouds-as-reflecting-boundaries (CRB) model, which assumes that the clouds are Lambertian reflectors. On average, improvements in the AMF calculation affect the tropospheric NO2 columns by ±15 % in winter and ±5 % in summer over largely polluted regions. In addition, the impact of aerosols on our tropospheric NO2 retrieval is investigated by comparing the concurrent retrievals based on ground-based aerosol measurements (explicit aerosol correction) and the aerosol-induced cloud parameters (implicit aerosol correction). Compared with the implicit aerosol correction utilising the CRB cloud parameters, the use of the CAL approach reduces the AMF errors by more than 10 %. Finally, to evaluate the improved GOME-2 tropospheric NO2 columns, a validation is performed using ground-based multi-axis differential optical absorption spectroscopy (MAXDOAS) measurements at different BIRA-IASB stations. At the suburban Xianghe station, the improved tropospheric NO2 dataset shows better agreement with coincident ground-based measurements with a correlation coefficient of 0.94.

In October 2017, the Sentinel-5 Precursor (S5P) mission was launched, carrying the TROPOspheric Monitoring Instrument (TROPOMI), which provides a daily global coverage at a spatial resolution as high as 7 km × 3.5 km and is expected to extend the European atmospheric composition record initiated with GOME/ERS-2 in 1995, enhancing our scientific knowledge of atmospheric processes with its unprecedented spatial resolution. Due to the ongoing need to understand and monitor the recovery of the ozone layer, as well as the evolution of tropospheric pollution, total ozone remains one of the leading species of interest during this mission. In this work, the TROPOMI near real time (NRTI) and offline (OFFL) total ozone column (TOC) products are presented and compared to daily ground-based quality-assured Brewer and Dobson TOC measurements deposited in the World Ozone and Ultraviolet Radiation Data Centre (WOUDC). Additional comparisons to individual Brewer measurements from the Canadian Brewer Network and the European Brewer Network (Eubrewnet) are performed. Furthermore, twilight zenith-sky measurements obtained with ZSL-DOAS (Zenith Scattered Light Differential Optical Absorption Spectroscopy) instruments, which form part of the SAOZ network (Système d'Analyse par Observation Zénitale), are used for the validation. The quality of the TROPOMI TOC data is evaluated in terms of the influence of location, solar zenith angle, viewing angle, season, effective temperature, surface albedo and clouds. For this purpose, globally distributed ground-based measurements have been utilized as the background truth. The overall statistical analysis of the global comparison shows that the mean bias and the mean standard deviation of the percentage difference between TROPOMI and ground-based TOC is within 0 –1.5 % and 2.5 %–4.5 %, respectively. The mean bias that results from the comparisons is well within the S5P product requirements, while the mean standard deviation is very close to those limits, especially considering that the statistics shown here originate both from the satellite and the ground-based measurements. Additionally, the TROPOMI OFFL and NRTI products are evaluated against already known spaceborne sensors, namely, the Ozone Mapping Profiler Suite, on board the Suomi National Polar-orbiting Partnership (OMPS/Suomi-NPP), NASA v2 TOCs, and the Global Ozone Monitoring Experiment 2 (GOME-2), on board the Metop-A (GOME-2/Metop-A) and Metop-B (GOME-2/Metop-B) satellites. This analysis shows a very good agreement for both TROPOMI products with well-established instruments, with the absolute differences in mean bias and mean standard deviation being below +0.7 % and 1 %, respectively. These results assure the scientific community of the good quality of the TROPOMI TOC products during its first year of operation and enhance the already prevalent expectation that TROPOMI/S5P will play a very significant role in the continuity of ozone monitoring from space.

We describe a formaldehyde (HCHO) retrieval algorithm for the Geostationary Environment Monitoring Spectrometer (GEMS) that will be launched by the Korean Ministry of Environment in 2019. The algorithm comprises three steps: preprocesses, radiance fitting, and postprocesses. The preprocesses include a wavelength calibration, as well as interpolation and convolution of absorption cross sections; radiance fitting is conducted using a nonlinear fitting method referred to as basic optical absorption spectroscopy (BOAS); and postprocesses include air mass factor calculations and bias corrections. In this study, several sensitivity tests are conducted to examine the retrieval uncertainties using the GEMS HCHO algorithm. We evaluate the algorithm with the Ozone Monitoring Instrument (OMI) Level 1B irradiance/radiance data by comparing our retrieved HCHO column densities with OMI HCHO products of the Smithsonian Astrophysical Observatory (OMHCHO) and of the Quality Assurance for Essential Climate Variables project (OMI QA4ECV). Results show that OMI HCHO slant columns retrieved using the GEMS algorithm are in good agreement with OMHCHO, with correlation coefficients of 0.77–0.91 and regression slopes of 0.94–1.04 for March, June, September, and December 2005. Spatial distributions of HCHO slant columns from the GEMS algorithm are consistent with the OMI QA4ECV products, but relatively poorer correlation coefficients of 0.52–0.76 are found compared to those against the OMHCHO products. Also, we compare the satellite results with ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) observations. OMI GEMS HCHO vertical columns are 9 %–25 % lower than those of MAX-DOAS at Haute-Provence Observatory (OHP) in France, Bremen in Germany, and Xianghe in China. We find that the OMI GEMS retrievals have less bias than the OMHCHO and OMI QA4ECV products at OHP and Bremen in comparison with MAX-DOAS.

In this study the consistency between MAX-DOAS measurements and radiative transfer simulations of the atmospheric O4 absorption is investigated on 2 mainly cloud-free days during the MAD-CAT campaign in Mainz, Germany, in summer 2013. In recent years several studies indicated that measurements and radiative transfer simulations of the atmospheric O4 absorption can only be brought into agreement if a so-called scaling factor () is applied to the measured O4 absorption. However, many studies, including those based on direct sunlight measurements, came to the opposite conclusion, that there is no need for a scaling factor. Up to now, there is no broad consensus for an explanation of the observed discrepancies between measurements and simulations. Previous studies inferred the need for a scaling factor from the comparison of the aerosol optical depths derived from MAX-DOAS O4 measurements with that derived from coincident sun photometer measurements. In this study a different approach is chosen: the measured O4 absorption at 360 nm is directly compared to the O4 absorption obtained from radiative transfer simulations. The atmospheric conditions used as input for the radiative transfer simulations were taken from independent data sets, in particular from sun photometer and ceilometer measurements at the measurement site. This study has three main goals: first all relevant error sources of the spectral analysis, the radiative transfer simulations and the extraction of the input parameters used for the radiative transfer simulations are quantified. One important result obtained from the analysis of synthetic spectra is that the O4 absorptions derived from the spectral analysis agree within 1 % with the corresponding radiative transfer simulations at 360 nm. Based on the results from sensitivity studies, recommendations for optimised settings for the spectral analysis and radiative transfer simulations are given. Second, the measured and simulated results are compared for 2 selected cloud-free days with similar aerosol optical depths but very different aerosol properties. On 18 June, measurements and simulations agree within their (rather large) uncertainties (the ratio of simulated and measured O4 absorptions is found to be 1.01±0.16). In contrast, on 8 July measurements and simulations significantly disagree: for the middle period of that day the ratio of simulated and measured O4 absorptions is found to be 0.82±0.10, which differs significantly from unity. Thus, for that day a scaling factor is needed to bring measurements and simulations into agreement. Third, recommendations for further intercomparison exercises are derived. One important recommendation for future studies is that aerosol profile data should be measured at the same wavelengths as the MAX-DOAS measurements. Also, the altitude range without profile information close to the ground should be minimised and detailed information on the aerosol optical and/or microphysical properties should be collected and used. The results for both days are inconsistent, and no explanation for a O4 scaling factor could be derived in this study. Thus, similar but more extended future studies should be performed, including more measurement days and more instruments. Also, additional wavelengths should be included.

Multi-axis differential optical absorption spectroscopy (MAX-DOAS) is a widely used measurement technique for the detection of a variety of atmospheric trace gases. Using inverse modelling, the observation of trace gas column densities along different lines of sight enables the retrieval of aerosol and trace gas vertical profiles in the atmospheric boundary layer using appropriate retrieval algorithms. In this study, the ability of eight profile retrieval algorithms to reconstruct vertical profiles is assessed on the basis of synthetic measurements. Five of the algorithms are based on the optimal estimation method, two on parametrised approaches, and one using an analytical approach without involving any radiative transfer modelling. The synthetic measurements consist of the median of simulated slant column densities of O4 at 360 and 477 nm, as well as of HCHO at 343 nm and NO2 at 477 nm, from seven datasets simulated by five different radiative transfer models. Simulations are performed for a combination of 10 trace gas and 11 aerosol profiles, as well as 11 elevation angles, three solar zenith, and three relative azimuth angles. Overall, the results from the different algorithms show moderate to good performance for the retrieval of vertical profiles, surface concentrations, and total columns. Except for some outliers, the root-mean-square difference between the true and retrieved state ranges between (0.05–0.1) km−1 for aerosol extinction and (2.5–5.0) ×1010 molec cm−3 for HCHO and NO2 concentrations.

The TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5 Precursor (S5P) satellite launched in October 2017 yields a wealth of atmospheric composition data, including retrievals of total column ozone (TCO3) that are provided in near-real-time (NRT) and off-line. The NRT TCO3 retrievals (v1.0.0–v1.1.2) have been included in the data assimilation system of the Copernicus Atmosphere Monitoring Service (CAMS), and tests to monitor the data and to carry out first assimilation experiments with them have been performed for the period 26 November 2017 to 30 November 2018. The TROPOMI TCO3 data agree to within 2 % with the CAMS analysis over large parts of the globe between 60∘ N and 60∘ S and also with TCO3 retrievals from the Ozone Monitoring Instrument (OMI) and the Global Ozone Monitoring Experiment-2 (GOME-2) that are routinely assimilated by CAMS. However, the TCO3 NRT data from TROPOMI show some retrieval anomalies at high latitudes, at low solar elevations and over snow/ice (e.g. Antarctica and snow-covered land areas in the Northern Hemisphere), where the differences with the CAMS analysis and the other data sets are larger. These differences are particularly pronounced over land in the NH during winter and spring (when they can reach up to 40 DU) and come mainly from the surface albedo climatology that is used in the NRT TROPOMI TCO3 retrieval. This climatology has a coarser horizontal resolution than the TROPOMI TCO3 data, which leads to problems in areas where there are large changes in reflectivity from pixel to pixel, e.g. pixels covered by snow/ice or not. The differences between TROPOMI and the CAMS analysis also show some dependency on scan position. The assimilation of TROPOMI TCO3 has been tested in the CAMS system for data between 60∘ N and 60∘ S and for solar elevations greater than 10∘ and is found to have a small positive impact on the ozone analysis compared to Brewer TCO3 data and an improved fit to ozone sondes in the tropical troposphere and to IAGOS aircraft profiles at West African airports. The impact of the TROPOMI data is relatively small because the CAMS analysis is already well constrained by several other ozone retrievals that are routinely assimilated. When averaged over the periods February–April and September–October 2018, differences between experiments with and without assimilation of TROPOMI data are less than 2 % for TCO3 and less than 3 % in the vertical for seasonal mean zonal mean O3 mixing ratios, with the largest relative differences found in the troposphere.

Satellite retrieval of vertical column densities (VCDs) of tropospheric nitrogen dioxide (NO2) is critical for NOx pollution and impact evaluation. For regions with high aerosol loadings, the retrieval accuracy is greatly affected by whether aerosol optical effects are treated implicitly (as additional “effective” clouds) or explicitly, among other factors. Our previous POMINO algorithm explicitly accounts for aerosol effects to improve the retrieval, especially in polluted situations over China, by using aerosol information from GEOS-Chem simulations with further monthly constraints by MODIS/Aqua aerosol optical depth (AOD) data. Here we present a major algorithm update, POMINO v1.1, by constructing a monthly climatological dataset of aerosol extinction profiles, based on level 2 CALIOP/CALIPSO data over 2007–2015, to better constrain the modeled aerosol vertical profiles. We find that GEOS-Chem captures the month-to-month variation in CALIOP aerosol layer height (ALH) but with a systematic underestimate by about 300–600 m (season and location dependent), due to a too strong negative vertical gradient of extinction above 1 km. Correcting the model aerosol extinction profiles results in small changes in retrieved cloud fraction, increases in cloud-top pressure (within 2 %–6 % in most cases), and increases in tropospheric NO2 VCD by 4 %–16 % over China on a monthly basis in 2012. The improved NO2 VCDs (in POMINO v1.1) are more consistent with independent ground-based MAX-DOAS observations (R2=0.80, NMB = −3.4 %, for 162 pixels in 49 days) than POMINO (R2=0.80, NMB = −9.6 %), DOMINO v2 (R2=0.68, NMB = −2.1 %), and QA4ECV (R2=0.75, NMB = −22.0 %) are. Especially on haze days, R2 reaches 0.76 for POMINO v1.1, much higher than that for POMINO (0.68), DOMINO v2 (0.38), and QA4ECV (0.34). Furthermore, the increase in cloud pressure likely reveals a more realistic vertical relationship between cloud and aerosol layers, with aerosols situated above the clouds in certain months instead of always below the clouds. The POMINO v1.1 algorithm is a core step towards our next public release of the data product (POMINO v2), and it will also be applied to the recently launched S5P-TROPOMI sensor.

We present an intercomparison study of four airborne imaging DOAS instruments, dedicated to the retrieval and high-resolution mapping of tropospheric nitrogen dioxide (NO2) vertical column densities (VCDs). The AROMAPEX campaign took place in Berlin, Germany, in April 2016 with the primary objective to test and intercompare the performance of experimental airborne imagers. The imaging DOAS instruments were operated simultaneously from two manned aircraft, performing synchronised flights: APEX (VITO–BIRA-IASB) was operated from DLR's DO-228 D-CFFU aircraft at 6.2 km in altitude, while AirMAP (IUP-Bremen), SWING (BIRA-IASB), and SBI (TNO–TU Delft–KNMI) were operated from the FUB Cessna 207T D-EAFU at 3.1 km. Two synchronised flights took place on 21 April 2016. NO2 slant columns were retrieved by applying differential optical absorption spectroscopy (DOAS) in the visible wavelength region and converted to VCDs by the computation of appropriate air mass factors (AMFs). Finally, the NO2 VCDs were georeferenced and mapped at high spatial resolution. For the sake of harmonising the different data sets, efforts were made to agree on a common set of parameter settings, AMF look-up table, and gridding algorithm. The NO2 horizontal distribution, observed by the different DOAS imagers, shows very similar spatial patterns. The NO2 field is dominated by two large plumes related to industrial compounds, crossing the city from west to east. The major highways A100 and A113 are also identified as line sources of NO2. Retrieved NO2 VCDs range between 1×1015 molec cm−2 upwind of the city and 20×1015 molec cm−2 in the dominant plume, with a mean of 7.3±1.8×1015 molec cm−2 for the morning flight and between 1 and 23×1015 molec cm−2 with a mean of 6.0±1.4×1015 molec cm−2 for the afternoon flight. The mean NO2 VCD retrieval errors are in the range of 22 % to 36 % for all sensors. The four data sets are in good agreement with Pearson correlation coefficients better than 0.9, while the linear regression analyses show slopes close to unity and generally small intercepts.

Global observations of tropospheric nitrogen dioxide (NO2) columns have been shown to be feasible from space, but consistent multi-sensor records do not yet exist, nor are they covered by planned activities at the international level. Harmonised, multi-decadal records of NO2 columns and their associated uncertainties can provide crucial information on how the emissions and concentrations of nitrogen oxides evolve over time. Here we describe the development of a new, community best-practice NO2 retrieval algorithm based on a synthesis of existing approaches. Detailed comparisons of these approaches led us to implement an enhanced spectral fitting method for NO2, a 1°  ×  1° TM5-MP data assimilation scheme to estimate the stratospheric background and improve air mass factor calculations. Guided by the needs expressed by data users, producers, and WMO GCOS guidelines, we incorporated detailed per-pixel uncertainty information in the data product, along with easily traceable information on the relevant quality aspects of the retrieval. We applied the improved QA4ECV NO2 algorithm to the most current level-1 data sets to produce a complete 22-year data record that includes GOME (1995–2003), SCIAMACHY (2002–2012), GOME-2(A) (2007 onwards) and OMI (2004 onwards). The QA4ECV NO2 spectral fitting recommendations and TM5-MP stratospheric column and air mass factor approach are currently also applied to S5P-TROPOMI. The uncertainties in the QA4ECV tropospheric NO2 columns amount to typically 40 % over polluted scenes. The first validation results of the QA4ECV OMI NO2 columns and their uncertainties over Tai'an, China, in June 2006 suggest a small bias (−2 %) and better precision than suggested by uncertainty propagation. We conclude that our improved QA4ECV NO2 long-term data record is providing valuable information to quantitatively constrain emissions, deposition, and trends in nitrogen oxides on a global scale.