Search

We introduce NitroNet, a deep learning model for the prediction of tropospheric NO2 profiles from satellite column measurements. NitroNet is a neural network trained on synthetic NO2 profiles from the regional chemistry and transport model WRF-Chem, which was operated on a European domain for the month of May 2019. This WRF-Chem simulation was constrained by in situ and satellite measurements, which were used to optimize important simulation parameters (e.g. the boundary layer scheme). The NitroNet model receives NO2 vertical column densities (VCDs) from the TROPOspheric Monitoring Instrument (TROPOMI) and ancillary variables (meteorology, emissions, etc.) as input, from which it reproduces NO2 concentration profiles. Training of the neural network is conducted on a filtered dataset, meaning that NO2 profiles showing strong disagreement (>20 %) with colocated TROPOMI column measurements are discarded. We present a first evaluation of NitroNet over a variety of geographical and temporal domains (Europe, the US West Coast, India, and China) and different seasons. For this purpose, we validate the NO2 profiles predicted by NitroNet against satellite, in situ, and MAX-DOAS (Multi-Axis Differential Optical Absorption Spectroscopy) measurements. The training data were previously validated against the same datasets. During summertime, NitroNet shows small biases and strong correlations with all three datasets: a bias of +6.7 % and R=0.95 for TROPOMI NO2 VCDs, a bias of −10.5 % and R=0.75 for AirBase surface concentrations, and a bias of −34.3 % to +99.6 % with R=0.83–0.99 for MAX-DOAS measurements. In comparison to TROPOMI satellite data, NitroNet even shows significantly lower errors and stronger correlation than a direct comparison with WRF-Chem numerical results. During wintertime considerable low biases arise because the summertime/late-spring training data are not fully representative of all atmospheric wintertime characteristics (e.g. longer NO2 lifetimes). Nonetheless, the wintertime performance of NitroNet is surprisingly good and comparable to that of classic regional chemistry and transport models. NitroNet can demonstrably be used outside the geographic and temporal domain of the training data with only slight performance reductions. What makes NitroNet unique when compared to similar existing deep learning models is the inclusion of synthetic model data, which offers important benefits: due to the lack of NO2 profile measurements, models trained on empirical datasets are limited to the prediction of surface concentrations learned from in situ measurements. NitroNet, however, can predict full tropospheric NO2 profiles. Furthermore, in situ measurements of NO2 are known to suffer from biases, often larger than +20 %, due to cross-sensitivities to photooxidants, which other models trained on empirical data inevitably reproduce.

Since 2011, the Global Energy and Water cycle Exchanges (GEWEX) Water Vapor Assessment (G-VAP) has provided performance analyses for state-of-the-art reanalysis and satellite water vapour products to the GEWEX Data and Analysis Panel (GDAP) and the user community in general. A significant component of the work undertaken by G-VAP is to characterise the quality and uncertainty of these water vapour records to (i) ensure full exploitation and (ii) avoid incorrect use or interpretation of results. This study presents results from the second phase of G-VAP, where we have extended and expanded our analysis of total column water vapour (TCWV) from phase 1, in conjunction with updating the G-VAP archive. For version 2 of the archive, we consider 28 freely available and mature satellite and reanalysis data products, remapped to a regular longitude–latitude grid of 2° × 2° and on monthly time steps between January 1979 and December 2019. We first analysed all records for a “common” short period of 5 years (2005–2009), focusing on variability (spatial and seasonal) and deviation from the ensemble mean. We observed that clear-sky daytime-only satellite products were generally drier than the ensemble mean, and seasonal variability/disparity in several regions up to 12 kg m−2 related to original spatial resolution and temporal sampling. For 11 of the 28 data records, further analysis was undertaken between 1988–2014. Within this “long period”, key results show (i) trends between −1.18 ± 0.68 to 3.82 ± 3.94 kg m−2 per decade and −0.39 ± 0.27 to 1.24 ± 0.85 kg m−2 per decade were found over ice-free global oceans and land surfaces, respectively, and (ii) regression coefficients of TCWV against surface temperatures of 6.17 ± 0.24 to 27.02 ± 0.51 % K−1 over oceans (using sea surface temperature) and 3.00 ± 0.17 to 7.77 ± 0.16 % K−1 over land (using surface air temperature). It is important to note that trends estimated within G-VAP are used to identify issues in the data records rather than analyse climate change. Additionally, breakpoints have been identified and characterised for both land and ocean surfaces within this period. Finally, we present a spatial analysis of correlations to six climate indices within the long period, highlighting regional areas of significant positive and negative correlation and the level of agreement among records.

We present a new method for estimating NOx emissions and effective lifetimes from large cities based on NO2 measurements from the TROPOspheric Monitoring Instrument (TROPOMI) (PAL dataset, May 2018–November 2021). As in previous studies, the estimate is based on the downwind plume evolution for different wind directions separately. The novelty of the presented approach lies in the simultaneous fit of downwind patterns for opposing wind directions, which makes the method far more robust (i.e., less prone to local minima with nonphysically high or low lifetimes) than a single exponential decay fit. In addition, the new method does not require the assumption of a city being a “point source” but also derives the spatial distribution of emissions. The method was successfully applied to 100 cities worldwide on a seasonal scale. Fitted emissions generally agree reasonably with the Emissions Database for Global Atmospheric Research (EDGAR) v6 (R=0.72) and are on average 14 % lower, while estimated uncertainties are still rather large (≈ 30 %–50 %). Lifetimes were found to be rather short (2.44 ± 0.68 h) and show no distinct dependency on season or latitude, which might be a consequence of discarding observations at high solar zenith angles (>65°). The main limitations of this and similar methods are the underlying assumptions of steady state (meaning constant emissions, wind fields and chemical conditions) within about 100 km downwind from a city, which is probably a simplification that is too strong in order to reach higher accuracies.

Satellite-derived spatiotemporal patterns of nitrogen oxide (NOx) emissions can improve accuracy of emission inventories to better support air quality and climate research and policy studies. In this study, we develop a new method by coupling the chemical transport Model-Independent SATellite-derived Emission estimation Algorithm for Mixed-sources (MISATEAM) with a divergence method to map high-resolution NOx emissions across US cities using TROPOspheric Monitoring Instrument (TROPOMI) tropospheric nitrogen dioxide (NO2) retrievals. The accuracy of the coupled method is validated through application to synthetic NO2 observations from the NASA-Unified Weather Research and Forecasting (NU-WRF) model, with a horizontal spatial resolution of 4 km × 4 km for 33 large and mid-size US cities. Validation reveals excellent agreement between inferred and NU-WRF-provided emission magnitudes (R= 0.99, normalized mean bias, NMB = −0.01) and a consistent spatial pattern when comparing emissions for individual grid cells (R=0.88±0.06). We then develop a TROPOMI-based database reporting annual emissions for 39 US cities at a horizontal spatial resolution of 0.05° × 0.05° from 2018 to 2021. This database demonstrates a strong correlation (R= 0.90) with the National Emission Inventory (NEI) but reveals some bias (NMB = −0.24). There are noticeable differences in the spatial patterns of emissions in some cities. Our analysis suggests that uncertainties in TROPOMI-based emissions and potential misallocation of emissions and/or missing sources in bottom-up emission inventories both contribute to these differences.

We present WRF-Chem simulations over central Europe with a spatial resolution of 3 km × 3 km and focus on nitrogen dioxide (NO2). A regional emission inventory issued by the German Environmental Agency, with a spatial resolution of 1 km × 1 km, is used as input. We demonstrate by comparison of five different model setups that significant improvements in model accuracy can be achieved by choosing the appropriate boundary layer scheme, increasing vertical mixing strength, and/or tuning the temporal modulation of the emission data (“temporal profiles”) driving the model. The model setup with improved vertical mixing is shown to produce the best results. Simulated NO2 surface concentrations are compared to measurements from a total of 275 in situ measurement stations in Germany, where the model was able to reproduce average noontime NO2 concentrations with a bias of ca. −3 % and R=0.74. The best agreement is achieved when correcting for the presumed NOy cross sensitivity of the molybdenum-based in situ measurements by computing an NOy correction factor from modelled peroxyacetyl nitrate (PAN) and nitric acid (HNO3) mixing ratios. A comparison between modelled NO2 vertical column densities (VCDs) and satellite observations from TROPOMI (TROPOspheric Monitoring Instrument) is conducted with averaging kernels taken into account. Simulations and satellite observations are shown to agree with a bias of +5.5 % and R=0.87 for monthly means. Lastly, simulated NO2 concentration profiles are compared to noontime NO2 profiles obtained from multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements at five locations in Europe. For stations within Germany, average biases of −25.3 % to +12.0 % were obtained. Outside of Germany, where lower-resolution emission data were used, biases of up to +50.7 % were observed. Overall, the study demonstrates the high sensitivity of modelled NO2 to the mixing processes in the boundary layer and the diurnal distribution of emissions.

Bromine monoxide (BrO) is a key radical in the atmosphere, influencing the chemical state of the atmosphere, most notably the abundance of ozone (O3). O3 depletion caused by the release of bromine has been observed and modeled in polar regions, salt pans, and in particular inside volcanic plumes. Furthermore, the molar ratio of BrO and SO2 – which can be detected simultaneously via spectroscopic measurements using the differential optical absorption spectroscopy (DOAS) method – is a proxy for the magmatic composition of a volcano and potentially an eruption forecast parameter. The detection of BrO in volcanic plumes from satellite spectroscopic observations is limited by the precision and sensitivity of the retrieval, which so far only allowed for the detection of BrO during major eruptions. The unprecedented spatial resolution of up to 3.5 km×5.5 km and the high signal-to-noise ratio of the TROPOspheric Monitoring Instrument (TROPOMI) on board Sentinel-5 Precursor (S-5P) enable observing and monitoring volcanic bromine release globally even for minor eruptions or even quiescent degassing. In this study, we investigate how far the BrO retrieval can be improved using TROPOMI data and how well BrO can be detected, even in small eruptions and during quiescent volcanic degassing. There are two steps for which improvements in accuracy are investigated and applied: the improvement and quantitative determination of (1) the detection limit of the DOAS BrO column retrieval and (2) the correction of the non-volcanic background BrO signal. First, the DOAS retrieval settings are varied, and their influence on accuracy and precision is investigated with respect to the detection limit and potential systematic influences. Based on these results, we propose a dedicated DOAS evaluation scheme optimized for the detection of BrO in volcanic plumes. For the DOAS retrieval, we propose the use of a large fit window from 323–360 nm, yielding a statistical uncertainty lower by a factor of 1.8 compared to previous BrO DOAS algorithms while not enhancing systematic influences. Second, the effect of the background BrO is reduced by a latitude-dependent empirical correction scheme correlated to cloud information as well as information on the O3 column. Via these improvements, the combined statistical and systematic uncertainties in the resulting BrO vertical column density is on the order of 7×1012moleculescm-2. We present a new and accurate retrieval algorithm of BrO columns from TROPOMI observations which allows for the detection of even slightly enhanced BrO amounts inside minor eruptive plumes of bromine-rich volcanoes. While designed specifically for TROPOMI observations, the retrieval algorithm is in general also applicable to other hyperspectral satellite observations. However, some parts might require adaptation.

We present an updated (v2) catalog of NOx emissions from point sources as derived from TROPOspheric Monitoring Instrument (TROPOMI) measurements of NO2 (Products Algorithm Laboratory (PAL) product) combined with wind fields from ERA5. Compared to version 1 of the catalog (Beirle et al., 2021), several improvements have been introduced to the algorithm. Most importantly, several corrections are applied, accounting for the effects of plume height on satellite sensitivity, 3D topographic effects, and the chemical loss of NOx, resulting in considerably higher and more accurate NOx emissions. In addition, error estimates are provided for each point source, taking into account the uncertainties of the individual retrieval steps. The v2 catalog is based on a fully automated iterative detection algorithm of point sources worldwide. It lists 1139 locations that have been found to be significant NOx sources. The majority of these locations match power plants listed in the Global Power Plant Database (GPPD). Other NOx point sources correspond to cement plants, metal smelters, industrial areas, or medium-sized cities. The emissions listed in v2 of the catalog show good agreement (within 20 % on average) to emissions reported by the German Environment Agency (Umweltbundesamt, UBA) as well as the United States Environmental Protection Agency (EPA). The data are publicly available at https://doi.org/10.26050/WDCC/No_xPointEmissionsV2 (Beirle et al., 2023).

We present a long-term data set of 1∘×1∘ monthly mean total column water vapour (TCWV) based on global measurements of the Ozone Monitoring Instrument (OMI) covering the time range from January 2005 to December 2020. In comparison to the retrieval algorithm of Borger et al. (2020), several modifications and filters have been applied accounting for instrumental issues (such as OMI's “row anomaly”) or the inferior quality of solar reference spectra. For instance, to overcome issues related to low-quality reference spectra, the daily solar irradiance spectrum is replaced by an annually varying mean earthshine radiance obtained in December over Antarctica. For the TCWV data set, we only consider measurements with an effective cloud fraction less than 20 %, an air mass factor (AMF) greater than 0.1, a snow- and ice-free ground pixel, and an OMI row that is not affected by the row anomaly over the complete time range of the data set. The individual TCWV measurements are then gridded to a regular 1∘×1∘ lattice, from which the monthly means are calculated. The investigation of sampling errors in the OMI TCWV data set shows that these are dominated by the clear-sky bias and cause on average deviations of around −10 %, which is consistent with the findings of previous studies. However, the spatiotemporal sampling errors and those due to the row-anomaly filter are negligible. In a comprehensive intercomparison study, we demonstrate that the OMI TCWV data set is in good agreement with the global reference data sets of ERA5 (fifth-generation ECMWF atmospheric reanalysis), RSS SSM/I (Remote Sensing Systems Special Sensor Microwave Imager), and CM SAF/CCI TCWV-global (COMBI): over ocean the orthogonal distance regressions indicate slopes close to unity with very small offsets and high coefficients of determination of around 0.96. However, over land, distinctive positive deviations of more than +10 kg m−2 are obtained for high TCWV values. These overestimations are mainly due to extreme overestimations of high TCWV values in the tropics, likely caused by uncertainties in the retrieval input data (surface albedo, cloud information) due to frequent cloud contamination in these regions. Similar results are found from intercomparisons with in situ radiosonde measurements from version 2 of the Integrated Global Radiosonde Archive (IGRA2) data set. Nevertheless, for TCWV values smaller than 25 kg m−2, the OMI TCWV data set shows very good agreement with the global reference data sets. Furthermore, a temporal stability analysis proves that the OMI TCWV data set is consistent with the temporal changes in the reference data sets and shows no significant deviation trends. As the TCWV retrieval can be easily applied to further satellite missions, additional TCWV data sets can be created from past missions, such as the Global Ozone Monitoring Experiment-1 (GOME-1) or the SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY); under consideration of systematic differences (e.g. due to different observation times), these data sets can be combined with the OMI TCWV data set in order to create a data record that would cover a time span from 1995 to the present. Moreover, the TCWV retrieval will also work for all missions dedicated to NO2 in the future, such as Sentinel-5 on MetOp-SG. The Max Planck Institute for Chemistry (MPIC) OMI total column water vapour (TCWV) climate data record (CDR) is available at https://doi.org/10.5281/zenodo.7973889 (Borger et al., 2023).

We investigate effects of the three-dimensional (3D) structure of volcanic plumes on the retrieval results of satellite and ground-based UV–Vis observations. For the analysis of such measurements, 1D scenarios are usually assumed (the atmospheric properties only depend on altitude). While 1D assumptions are well suited for the analysis of many atmospheric phenomena, they are usually less appropriate for narrow trace gas plumes. For UV–Vis satellite instruments with large ground pixel sizes like the Global Ozone Monitoring Experiment-2 (GOME-2), the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) or the Ozone Monitoring Instrument (OMI), 3D effects are of minor importance, but usually these observations are not sensitive to small volcanic plumes. In contrast, observations of the TROPOspheric Monitoring Instrument (TROPOMI) on board Sentinel-5P have a much smaller ground pixel size (3.5 × 5.5 km2). Thus, on the one hand, TROPOMI can detect much smaller plumes than previous instruments. On the other hand, 3D effects become more important, because the TROPOMI ground pixel size is smaller than the height of the troposphere and also smaller than horizontal atmospheric photon path lengths in the UV–Vis spectral range. In this study we investigate the following 3D effects using Monte Carlo radiative transfer simulations: (1) the light-mixing effect caused by horizontal photon paths, (2) the saturation effect for strong SO2 absorption, (3) geometric effects related to slant illumination and viewing angles and (4) plume side-effects related to slant illumination angles and photons reaching the sensor from the sides of volcanic plumes. The first two effects especially can lead to a strong and systematic underestimation of the true trace gas content if 1D retrievals are applied (more than 50 % for the light-mixing effect and up to 100 % for the saturation effect). Besides the atmospheric radiative transfer, the saturation effect also affects the spectral retrievals. Geometric effects have a weaker influence on the quantitative analyses but can lead to a spatial smearing of elevated plumes or even to virtual double plumes. Plume side-effects are small for short wavelengths but can become large for longer wavelengths (up to 100 % for slant viewing and illumination angles). For ground-based observations, most of the above-mentioned 3D effects are not important because of the narrow field of view (FOV) and the closer distance between the instrument and the volcanic plume. However, the light-mixing effect shows a similar strong dependence on the horizontal plume extension as for satellite observations and should be taken into account for the analysis of ground-based observations.

In this study, we investigate trends in total column water vapour (TCWV) retrieved from measurements of the Ozone Monitoring Instrument (OMI) for the time range between January 2005 to December 2020. The trend analysis reveals, on global average, an annual increase in the TCWV amount of approximately +0.054 kg m−2 yr−1 or +0.21 % yr−1. After the application of a Z test (to the significance level of 5 %) and a false discovery rate (FDR) test to the results of the trend analysis, mainly positive trends remain, in particular over the northern subtropics in the eastern Pacific. Combining the relative TCWV trends with trends in air temperature, we also analyse trends in relative humidity (RH) on the local scale. This analysis reveals that the assumption of temporally invariant RH is not always fulfilled, as we obtain increasing and decreasing RH trends over large areas of the ocean and land surface and also observe that these trends are not limited to arid and humid regions, respectively. For instance, we find decreasing RH trends over the (humid) tropical Pacific Ocean in the region of the Intertropical Convergence Zone. Interestingly, these decreasing RH trends in the tropical Pacific Ocean coincide well with decreasing trends in precipitation. Moreover, by combining the trends of TCWV, surface temperature, and precipitation, we derive trends for the global water vapour turnover time (TUT) of approximately +0.02 d yr−1. Also, we obtain a TUT rate of change of around 8.4 % K−1, which is 2 to 3 times higher than the values obtained in previous studies.

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.

We present a formalism that relates the vertical column density (VCD) of the oxygen collision complex O2–O2 (denoted as O4 below) to surface (2 m) values of temperature and pressure, based on physical laws. In addition, we propose an empirical modification which also accounts for surface relative humidity (RH). This allows for simple and quick calculation of the O4 VCD without the need for constructing full vertical profiles. The parameterization reproduces the true O4 VCD, as derived from vertically integrated profiles, within -0.7±1.2% (mean ± SD) for Weather Research and Forecasting (WRF) simulations around Germany, 0.2±1.8 % for global reanalysis data (ERA5), and -0.3±1.4% for Global Climate Observing System (GCOS) Reference Upper-Air Network (GRUAN) radiosonde measurements around the world. When applied to measured surface values, uncertainties of 1 K, 1 hPa, and 16 % for temperature, pressure, and RH correspond to relative uncertainties of the O4 VCD of 0.3 %, 0.2 %, and 1 %, respectively. The proposed parameterization thus provides a simple and accurate formula for the calculation of the O4 VCD which is expected to be useful in particular for MAX-DOAS applications.

We present a new method to infer emissions and lifetimes of nitrogen oxides (NOx) based on tropospheric nitrogen dioxide (NO2) observations together with reanalysis wind fields for cities located in polluted backgrounds. Since the accuracy of the method is difficult to assess due to lack of “true values” that can be used as a benchmark, we apply the method to synthetic NO2 observations derived from the NASA-Unified Weather Research and Forecasting (NU-WRF) model at a high horizontal spatial resolution of 4 km × 4 km for cities over the continental United States. We compare the inferred emissions and lifetimes with the values given by the NU-WRF model to evaluate the method. The method is applicable to 26 US cities. The derived results are generally in good agreement with the values given by the model, with the relative differences of 2 % ± 17 % (mean ± standard deviation) and 15 % ± 25 % for lifetimes and emissions, respectively. Our investigation suggests that the use of wind data prior to the satellite overpass time improves the performance of the method. The correlation coefficients between inferred and NU-WRF lifetimes increase from 0.56 to 0.79 and for emissions increase from 0.88 to 0.96 when comparing results based on wind fields sampled simultaneously with satellite observations and averaged over 9 h data prior to satellite observations, respectively. We estimate that uncertainties in NOx lifetime and emissions arising from the method are approximately 15 % and 20 %, respectively, for typical (US) cities. The total uncertainties reach up to 43 % (lifetimes) and 45 % (emissions) by considering the additional uncertainties associated with satellite NO2 observations and wind data. We expect this new method to be applicable to NO2 observations from the TROPOspheric Monitoring Instrument (TROPOMI) and geostationary satellites, such as Geostationary Environment Monitoring Spectrometer (GEMS) or the Tropospheric Emissions: Monitoring Pollution (TEMPO) instrument, to estimate urban NOx emissions and lifetimes globally.

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 %.

Iodine plays a vital role in oxidation chemistry over Antarctica, with past observations showing highly elevated levels of iodine oxide (IO) leading to severe depletion of boundary layer ozone in West Antarctica. Here, we present MAX-DOAS-based (multi-axis differential absorption spectroscopy) observations of IO over three summers (2015–2017) at the Indian Antarctic bases of Bharati and Maitri. IO was observed during all the campaigns with mixing ratios below 2 pptv (parts per trillion by volume) for the three summers, which are lower than the peak levels observed in West Antarctica. This suggests that sources in West Antarctica are different or stronger than sources of iodine compounds in East Antarctica, the nature of which is still uncertain. Vertical profiles estimated using a profile retrieval algorithm showed decreasing gradients with a peak in the lower boundary layer. The ground-based instrument retrieved vertical column densities (VCDs) were approximately a factor of 3 to 5 higher than the VCDs reported using satellite-based instruments, which is most likely related to the sensitivities of the measurement techniques. Air mass back-trajectory analysis failed to highlight a source region, with most of the air masses coming from coastal or continental regions. This study highlights the variation in iodine chemistry in different regions in Antarctica and the importance of a long-term dataset to validate models estimating the impacts of iodine chemistry.

Ground-based Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) is a state-of-the-art remote sensing technique for deriving vertical profiles of trace gases and aerosols. However, MAX-DOAS profile inversions under aerosol pollution scenarios are challenging because of the complex radiative transfer and limited information content of the measurements. In this study, the performances of two inversion algorithms were evaluated for various aerosol pollution scenarios based on synthetic slant column densities (SCDs) derived from radiative transfer simulations. Compared to previous studies, in our study, much larger ranges of aerosol optical depth (AOD) and NO2 vertical column densities (VCDs) are covered. One inversion algorithm is based on optimal estimation; the other uses a parameterized approach. In this analysis, three types of profile shapes for aerosols and NO2 were considered: exponential, Boltzmann, and Gaussian. First, the systematic deviations of the retrieved aerosol profiles from the input profiles were investigated. For most cases, the AODs of the retrieved profiles were found to be systematically lower than the input values, and the deviations increased with increasing AOD. In particular for the optimal estimation algorithm and for high AOD, these findings are consistent with the results in previous studies. The assumed single scattering albedo (SSA) and asymmetry parameter (AP) have a systematic influence on the aerosol retrieval. However, for most cases the influence of the assumed SSA and AP on the retrieval results are rather small (compared to other uncertainties). For the optimal estimation algorithm, the agreement with the input values can be improved by optimizing the covariance matrix of the a priori uncertainties. Second, the aerosol effects on the NO2 profile retrieval were tested. Here, especially for the optimal estimation algorithm, a systematic dependence on the NO2 VCD was found, with a strong relative overestimation of the retrieved results for low NO2 VCDs and an underestimation for high NO2 VCDs. In contrast, the dependence on the aerosol profiles was found to be rather low. Interestingly, the results for both investigated wavelengths (360 and 477 nm) were found to be rather similar, indicating that the differences in the radiative transfer between both wavelengths have no strong effect. In general, both inversion schemes can retrieve the near-surface values of aerosol extinction and trace gas concentrations well.

We present high spatial resolution (up to 2.2×2.2 km2) simulations focussed over south-west Germany using the online coupled regional atmospheric chemistry model system MECO(n) (MESSy-fied ECHAM and COSMO models nested n times). Numerical simulation of nitrogen dioxide (NO2) surface volume mixing ratios (VMRs) are compared to in situ measurements from a network with 193 locations including background, traffic-adjacent and industrial stations to investigate the model's performance in simulating the spatial and temporal variability of short-lived chemical species. We show that the use of a high-resolution and up-to-date emission inventory is crucial for reproducing the spatial variability and resulted in good agreement with the measured VMRs at the background and industrial locations with an overall bias of less than 10 %. We introduce a computationally efficient approach that simulates diurnal and daily variability in monthly-resolved anthropogenic emissions to resolve the temporal variability of NO2. MAX-DOAS (Multiple AXis Differential Optical Absorption Spectroscopy) measurements performed at Mainz (49.99∘ N, 8.23∘ E) were used to evaluate the simulated tropospheric vertical column densities (VCDs) of NO2. We propose a consistent and robust approach to evaluate the vertical distribution of NO2 in the boundary layer by comparing the individual differential slant column densities (dSCDs) at various elevation angles. This approach considers details of the spatial heterogeneity and sensitivity volume of the MAX-DOAS measurements while comparing the measured and simulated dSCDs. The effects of clouds on the agreement between MAX-DOAS measurements and simulations have also been investigated. For low elevation angles (≤8∘), small biases in the range of −14 % to +7 % and Pearson correlation coefficients in the range of 0.5 to 0.8 were achieved for different azimuth directions in the cloud-free cases, indicating good model performance in the layers close to the surface. Accounting for diurnal and daily variability in the monthly-resolved anthropogenic emissions was found to be crucial for the accurate representation of time series of measured NO2 VMR and dSCDs and is particularly critical when vertical mixing is suppressed, and the atmospheric lifetime of NO2 is relatively long.

Clouds impact the radiative transfer of the Earth's atmosphere and strongly influence satellite measurements in the ultraviolet–visible (UV–vis) and infrared (IR) spectral ranges. For satellite measurements of trace gases absorbing in the UV–vis spectral range, particularly clouds ultimately determine the vertical sensitivity profile, mainly by reducing the sensitivity for trace-gas columns below the cloud. The Mainz iterative cloud retrieval utilities (MICRU) algorithm is specifically designed to reduce the error budget of trace-gas retrievals, such as those for nitrogen dioxide (NO2), which strongly depends on the accuracy of small cloud fractions (CFs) in particular. The accuracy of MICRU is governed by an empirical parameterisation of the viewing-geometry-dependent background surface reflectivity taking instrumental and physical effects into account. Instrumental effects are mainly degradation and polarisation effects; physical effects are due to the anisotropy of the surface reflectivity, e.g. shadowing of plants and sun glitter. MICRU is applied to main science channel (MSC) and polarisation measurement device (PMD) data collected between April 2007 and June 2013 by the Global Ozone Monitoring Experiment 2A (GOME-2A) instrument aboard the MetOp-A satellite. CFs are retrieved at different spectral bands between 374 and 758 nm. The MICRU results for MSC and PMD at different wavelengths are intercompared to study CF precision and accuracy, which depend on wavelength and spatial correlation. Furthermore, MICRU results are compared to FRESCO (fast retrieval scheme for clouds from the oxygen A band) and OCRA (optical cloud recognition algorithm) operational cloud products. We show that MICRU retrieves small CFs with an accuracy of 0.04 or better for the entire 1920 km wide swath with a potential bias between −0.01 and −0.03. CFs retrieved at shorter wavelengths are less affected by adverse surface heterogeneities. The comparison to the operational CF algorithms shows that MICRU significantly reduces the dependence on viewing angle, time, and sun glitter. Systematic effects along coasts are particularly small for MICRU due to its dedicated treatment of land and ocean surfaces. The MICRU algorithm is designed for spectroscopic instruments ranging from the GOME to Sentinel-5P/Tropospheric Monitoring Instrument (TROPOMI) but is also applicable to UV–vis imagers like the Advanced Very High Resolution Radiometer (AVHRR), the Moderate Resolution Imaging Spectroradiometer (MODIS), the Visible Infrared Imaging Radiometer Suite (VIIRS), and Sentinel-2.

We present version 1.0 of a global catalog of NOx emissions from point sources, derived from TROPOspheric Monitoring Instrument (TROPOMI) measurements of tropospheric NO2 for 2018–2019. The identification of sources and quantification of emissions are based on the divergence (spatial derivative) of the mean horizontal flux, which is highly sensitive for point sources like power plant exhaust stacks. The catalog lists 451 locations which could be clearly identified as NOx point sources by a fully automated algorithm, while ambiguous cases as well as area sources such as megacities are skipped. A total of 242 of these point sources could be automatically matched to power plants. Other NOx point sources listed in the catalog are metal smelters, cement plants, or industrial areas. The four largest localized NOx emitters are all coal combustion plants in South Africa. About 1/4 of all detected point sources are located in the Indian subcontinent and are mostly associated with power plants. The catalog is incomplete, mainly due to persisting gaps in the TROPOMI NO2 product at some coastlines, inaccurate or complex wind fields in coastal and mountainous regions, and high noise in the divergence maps for high background pollution. The derived emissions are generally too low, lacking a factor of about 2 up to 8 for extreme cases. This strong low bias results from combination of different effects, most of all a strong underestimation of near-surface NO2 in TROPOMI NO2 columns. Still, the catalog has high potential for checking and improving emission inventories, as it provides accurate and independent up-to-date information on the location of sources of NOx and thus also CO2. The catalog of NOx emissions from point sources is freely available at https://doi.org/10.26050/WDCC/Quant_NOx_TROPOMI (Beirle et al., 2020).

In this study, we compare measured and simulated O4 absorptions for conditions of extremely low aerosol optical depth (between 0.034 to 0.056 at 360 nm) on one day during a ship cruise in the tropical Atlantic. For such conditions, the uncertainties related to imperfect knowledge of aerosol properties do not significantly affect the comparison results. We find that the simulations underestimate the measurements by 15 % to 20 %. Even for simulations without any aerosols, the measured O4 absorptions are still systematically higher than the simulation results. The observed discrepancies cannot be explained by uncertainties of the measurements and simulations and thus indicate a fundamental inconsistency between simulations and measurements.

We used a global long-term (1995–2015) data set of total column water vapour (TCWV) derived from satellite observations to quantify to which extent the temporal patterns of various teleconnections can be identified in this data set. To our knowledge, such a comprehensive global TCWV data set was rarely used for teleconnection studies. One important property of the TCWV data set is that it is purely based on observational data. We developed a new empirical method to decide whether a teleconnection index is significantly detected in the global data set. We compared our new method to well-established hypothesis tests and found good agreement with the results of our approach. Based on our empirical method more than 40 teleconnection indices were significantly detected in the global TCWV data set derived from satellite observations. In addition to the satellite data we also applied our method to other global data sets derived from ERA-Interim. One important finding is that the spatial patterns obtained for the ERA TCWV data are very similar to the observational TCWV data set indicating a high consistency between the satellite and ERA data. Moreover, similar results are also found for two selections of ERA data (either all data or mainly clear-sky data). This finding indicates that the clear-sky bias of the satellite data set is negligible for the results of this study. However, for some indices, also systematic differences in the spatial patterns between the satellite and model data set were found probably indicating possible shortcomings in the model data. For most “traditional” teleconnection data sets (surface temperature, surface pressure, geopotential heights and meridional winds at different altitudes) a smaller number of significant teleconnection indices was found than for the TCWV data sets, while for zonal winds at different altitudes, the number of significant teleconnection indices (up to > 50) was higher. The strongest teleconnection signals were found in the data sets of tropospheric geopotential heights and surface pressure. In all global data sets, no “other indices” (solar variability, stratospheric AOD or hurricane frequency) were significantly detected. Since many teleconnection indices are strongly correlated, we also applied our method to a set of orthogonalised indices, which represent the dominant independent temporal teleconnection patterns. The number of significantly detected orthogonalised indices (20) was found to be much smaller than for the original indices (42). Based on the orthogonalised indices we derived the global spatial distribution of the cumulative effect of teleconnections. The strongest effect on the TCWV is found in the tropics and high latitudes.

In urban areas, road traffic is a dominant source of nitrogen oxides (NOx=NO+NO2). Although the emissions from individual vehicles are regulated by the European emission standards, real driving emissions often exceed these limits. In this study, two multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments on opposite sides of the motorway were used to measure the NO2 absorption caused by road traffic at the A60 motorway close to Mainz, Germany. In combination with wind data, the total NOx emissions for the occurring traffic volume can be estimated. Hereto, the ozone-dependent photochemical equilibrium between NO and NO2 is considered. We show that for 10 May 2019 the measured emissions exceed the maximum expected emissions calculated from the European emission standards for standardised test cycles by a factor of 11±7. One major advantage of the method used here is that MAX-DOAS measurements are very sensitive to the integrated NO2 concentration close to the surface. Thus, all emitted NO2 molecules are detected independently from their altitude, and therefore the whole emission plume originating from the nearby motorway is captured, which is a key advantage compared to other approaches such as in situ measurements.

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.

We present comprehensive long-term ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements of aerosols, nitrogen dioxide (NO2), and formaldehyde (HCHO) from Mohali (30.667∘ N, 76.739∘ E; ∼310 m above mean sea level), located in the densely populated Indo-Gangetic Plain (IGP) of India. We investigate the temporal variation in tropospheric columns, surface volume mixing ratio (VMR), and vertical profiles of aerosols, NO2, and HCHO and identify factors driving their ambient levels and distributions for the period from January 2013 to June 2017. We observed mean aerosol optical depth (AOD) at 360 nm, tropospheric NO2 vertical column density (VCD), and tropospheric HCHO VCD for the measurement period to be 0.63 ± 0.51, (6.7 ± 4.1) × 1015, and (12.1 ± 7.5) × 1015 molecules cm−2, respectively. Concerning the tropospheric NO2 VCDs, Mohali was found to be less polluted than urban and suburban locations of China and western countries, but comparable HCHO VCDs were observed. For the more than 4 years of measurements during which the region around the measurement location underwent significant urban development, we did not observe obvious annual trends in AOD, NO2, and HCHO. High tropospheric NO2 VCDs were observed in periods with enhanced biomass and biofuel combustion (e.g. agricultural residue burning and domestic burning for heating). Highest tropospheric HCHO VCDs were observed in agricultural residue burning periods with favourable meteorological conditions for photochemical formation, which in previous studies have shown an implication for high ambient ozone also over the IGP. Highest AOD is observed in the monsoon season, indicating possible hygroscopic growth of the aerosol particles. Most of the NO2 is located close to the surface, whereas significant HCHO is present at higher altitudes up to 600 m during summer indicating active photochemistry at high altitudes. The vertical distribution of aerosol, NO2, and HCHO follows the change in boundary layer height (BLH), from the ERA5 dataset of European Centre for Medium-Range Weather Forecasts, between summer and winter. However, deep convection during the monsoon transports the pollutants at high altitudes similar to summer despite a shallow ERA5 BLH. Strong gradients in the vertical profiles of HCHO are observed during the months when primary anthropogenic sources dominate the formaldehyde production. High-resolution MODIS AOD measurements correlate well but were systematically higher than MAX-DOAS AODs. The ground-based MAX-DOAS measurements were used to evaluate three NO2 data products and two HCHO data products of the Ozone Monitoring Instrument (OMI) for the first time over India and the IGP. NO2 VCDs from OMI correlate reasonably with MAX-DOAS VCDs but are lower by ∼30 %–50 % due to the difference in vertical sensitivities and the rather large OMI footprint. OMI HCHO VCDs exceed the MAX-DOAS VCDs by up to 30 %. We show that there is significant scope for improvement in the a priori vertical profiles of trace gases, which are used in OMI retrievals. The difference in vertical representativeness was found to be crucial for the observed biases in NO2 and HCHO surface VMR intercomparisons. Using the ratio of NO2 and HCHO VCDs measured from MAX-DOAS, we have found that the peak daytime ozone production regime is sensitive to both NOx and VOCs in winter but strongly sensitive to NOx in other seasons.

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.

Total column water vapour has been retrieved from TROPOMI measurements in the visible blue spectral range and compared to a variety of different reference data sets for clear-sky conditions during boreal summer and winter. The retrieval consists of the common two-step DOAS approach: first the spectral analysis is performed within a linearized scheme and then the retrieved slant column densities are converted to vertical columns using an iterative scheme for the water vapour a priori profile shape, which is based on an empirical parameterization of the water vapour scale height. Moreover, a modified albedo map was used combining the OMI LER albedo and scaled MODIS albedo map. The use of the alternative albedo is especially important over regions with very low albedo and high probability of clouds like the Amazon region. The errors of the total column water vapour (TCWV) retrieval have been theoretically estimated considering the contribution of a variety of different uncertainty sources. For observations during clear-sky conditions, over ocean surface, and at low solar zenith angles the error typically is around values of 10 %–20 %, and during cloudy-sky conditions, over land surface, and at high solar zenith angles it reaches values around 20 %–50 %. In the framework of a validation study the retrieval demonstrates that it can well capture the global water vapour distribution: the retrieved H2O vertical column densities (VCDs) show very good agreement with the reference data sets over ocean for boreal summer and winter whereby the modified albedo map substantially improves the retrieval's consistency to the reference data sets, in particular over tropical land masses. However, over land the retrieval underestimates the VCD by about 10 %, particularly during summertime. Our investigations show that this underestimation is likely caused by uncertainties within the surface albedo and the cloud input data: low-level clouds cause an underestimation, but for mid- to high-level clouds good agreement is found. In addition, our investigations indicate that these biases can probably be further reduced by the use of improved cloud input data. For the general purpose we recommend only using VCDs with cloud fraction % and AMF >0.1, which represents a good compromise between spatial coverage and retrieval accuracy. The TCWV retrieval can be easily applied to further satellite sensors (e.g. GOME-2 or OMI) for creating uniform, long-term measurement data sets, which is particularly interesting for climate and trend studies of water vapour.

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.

We present a method to infer CO2 emissions from individual power plants based on satellite observations of co-emitted nitrogen dioxide (NO2), which could serve as complementary verification of bottom-up inventories or be used to supplement these inventories. We demonstrate its utility on eight large and isolated US power plants, where accurate stack emission estimates of both gases are available for comparison. In the first step of our methodology, we infer nitrogen oxides (NOx) emissions from US power plants using Ozone Monitoring Instrument (OMI) NO2 tropospheric vertical column densities (VCDs) averaged over the ozone season (May–September) and a “top-down” approach that we previously developed. Second, we determine the relationship between NOx and CO2 emissions based on the direct stack emissions measurements reported by continuous emissions monitoring system (CEMS) programs, accounting for coal quality, boiler firing technology, NOx emission control device type, and any change in operating conditions. Third, we estimate CO2 emissions for power plants using the OMI-estimated NOx emissions and the CEMS NOx∕CO2 emission ratio. We find that the CO2 emissions estimated by our satellite-based method during 2005–2017 are in reasonable agreement with the US CEMS measurements, with a relative difference of 8 %±41 % (mean ± standard deviation). The broader implication of our methodology is that it has the potential to provide an additional constraint on CO2 emissions from power plants in regions of the world without reliable emissions accounting. We explore the feasibility by comparing the derived NOx∕CO2 emission ratios for the US with those from a bottom-up emission inventory for other countries and applying our methodology to a power plant in South Africa, where the satellite-based emission estimates show reasonable consistency with other independent estimates. Though our analysis is limited to a few power plants, we expect to be able to apply our method to more US (and world) power plants when multi-year data records become available from new OMI-like sensors with improved capabilities, such as the TROPOspheric Monitoring Instrument (TROPOMI), and upcoming geostationary satellites, such as the Tropospheric Emissions: Monitoring Pollution (TEMPO) instrument.

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 Mainz profile algorithm (MAPA) derives vertical profiles of aerosol extinction and trace gas concentrations from MAX-DOAS measurements of slant column densities under multiple elevation angles. This paper presents (a) a detailed description of the MAPA (v0.98), (b) results for the CINDI-2 campaign, and (c) sensitivity studies on the impact of a priori assumptions such as flag thresholds. Like previous profile retrieval schemes developed at MPIC, MAPA is based on a profile parameterization combining box profiles, which also might be lifted, and exponential profiles. But in contrast to previous inversion schemes based on least-square fits, MAPA follows a Monte Carlo approach for deriving those profile parameters yielding best match to the MAX-DOAS observations. This is much faster and directly provides physically meaningful distributions of profile parameters. In addition, MAPA includes an elaborated flagging scheme for the identification of questionable or dubious results. The AODs derived with MAPA for the CINDI-2 campaign show good agreement with AERONET if a scaling factor of 0.8 is applied for O4, and the respective NO2 and HCHO surface mixing ratios match those derived from coincident long-path DOAS measurements. MAPA results are robust with respect to modifications of the a priori MAPA settings within plausible limits.

An improved algorithm for the retrieval of total and tropospheric nitrogen dioxide (NO2) columns from the Global Ozone Monitoring Experiment-2 (GOME-2) is presented. The refined retrieval will be implemented in a future version of the GOME Data Processor (GDP) as used by the EUMETSAT Satellite Application Facility on Atmospheric Composition and UV Radiation (AC-SAF). The first main improvement is the application of an extended 425–497 nm wavelength fitting window in the differential optical absorption spectroscopy (DOAS) retrieval of the NO2 slant column density, based on which initial total NO2 columns are computed using stratospheric air mass factors (AMFs). Updated absorption cross sections and a linear offset correction are used for the large fitting window. An improved slit function treatment is applied to compensate for both long-term and in-orbit drift of the GOME-2 slit function. Compared to the current operational (GDP 4.8) dataset, the use of these new features increases the NO2 columns by ∼1–3×1014 molec cm2 and reduces the slant column error by ∼24 %. In addition, the bias between GOME-2A and GOME-2B measurements is largely reduced by adopting a new level 1b data version in the DOAS retrieval. The retrieved NO2 slant columns show good consistency with the Quality Assurance for Essential Climate Variables (QA4ECV) retrieval with a good overall quality. Second, the STRatospheric Estimation Algorithm from Mainz (STREAM), which was originally developed for the TROPOspheric Monitoring Instrument (TROPOMI) instrument, was optimised for GOME-2 measurements to determine the stratospheric NO2 column density. Applied to synthetic GOME-2 data, the estimated stratospheric NO2 columns from STREAM shows good agreement with the a priori truth. An improved latitudinal correction is introduced in STREAM to reduce the biases over the subtropics. Applied to GOME-2 measurements, STREAM largely reduces the overestimation of stratospheric NO2 columns over polluted regions in the GDP 4.8 dataset. Third, the calculation of AMF applies an updated box-air-mass factor (box-AMF) look-up table (LUT) calculated using the latest version 2.7 of the Vector-LInearized Discrete Ordinate Radiative Transfer (VLIDORT) model with an increased number of reference points and vertical layers, a new GOME-2 surface albedo climatology, and improved a priori NO2 profiles obtained from the TM5-MP chemistry transport model. A large effect (mainly enhancement in summer and reduction in winter) on the retrieved tropospheric NO2 columns by more than 10 % is found over polluted regions. To evaluate the GOME-2 tropospheric NO2 columns, an end-to-end validation is performed using ground-based multiple-axis DOAS (MAXDOAS) measurements. The validation is illustrated for six stations covering urban, suburban, and background situations. Compared to the GDP 4.8 product, the new dataset presents improved agreement with the MAXDOAS measurements for all the stations.

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.

Nitrogen oxides (NOx≡NO+NO2) in the upper troposphere (UT) have a large impact on global tropospheric ozone and OH (the main atmospheric oxidant). New cloud-sliced observations of UT NO2 at 450–280 hPa (∼6–9 km) from the Ozone Monitoring Instrument (OMI) produced by NASA and the Royal Netherlands Meteorological Institute (KNMI) provide global coverage to test our understanding of the factors controlling UT NOx. We find that these products offer useful information when averaged over coarse scales (20∘×32∘, seasonal), and that the NASA product is more consistent with aircraft observations of UT NO2. Correlation with Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD) satellite observations of lightning flash frequencies suggests that lightning is the dominant source of NOx to the upper troposphere except for extratropical latitudes in winter. The NO2 background in the absence of lightning is 10–20 pptv. We infer a global mean NOx yield of 280±80 moles per lightning flash, with no significant difference between the tropics and midlatitudes, and a global lightning NOx source of 5.9±1.7 Tg N a−1. There is indication that the NOx yield per flash increases with lightning flash footprint and with flash energy.

Nitrogen dioxide (NO2) and formaldehyde (HCHO) column data from satellite instruments are used for air quality and climate studies. Both NO2 and HCHO have been identified as precursors to the ozone (O3) and aerosol essential climate variables, and it is essential to quantify and characterise their uncertainties. Here we present an intercomparison of NO2 and HCHO slant column density (SCD) retrievals from four different research groups (BIRA-IASB, IUP Bremen, and KNMI as part of the Quality Assurance for Essential Climate Variables (QA4ECV) project consortium, and NASA) and from the OMI and GOME-2A instruments. Our evaluation is motivated by recent improvements in differential optical absorption spectroscopy (DOAS) fitting techniques and by the desire to provide a fully traceable uncertainty budget for the climate data record generated within QA4ECV. The improved NO2 and HCHO SCD values are in close agreement but with substantial differences in the reported uncertainties between groups and instruments. To check the DOAS uncertainties, we use an independent estimate based on the spatial variability of the SCDs within a remote region. For NO2, we find the smallest uncertainties from the new QA4ECV retrieval (0.8 × 1015 molec. cm−2 for both instruments over their mission lifetimes). Relative to earlier approaches, the QA4ECV NO2 retrieval shows better agreement between DOAS and statistical uncertainty estimates, suggesting that the improved QA4ECV NO2 retrieval has reduced but not altogether eliminated systematic errors in the fitting approach. For HCHO, we reach similar conclusions (QA4ECV uncertainties of 8–12 × 1015 molec. cm−2), but the closeness between the DOAS and statistical uncertainty estimates suggests that HCHO uncertainties are indeed dominated by random noise from the satellite's level 1 data. We find that SCD uncertainties are smallest for high top-of-atmosphere reflectance levels with high measurement signal-to-noise ratios. From 2005 to 2015, OMI NO2 SCD uncertainties increase by 1–2 % year−1, which is related to detector degradation and stripes, but OMI HCHO SCD uncertainties are remarkably stable (increase −1) and this is related to the use of Earth radiance reference spectra which reduces stripes. For GOME-2A, NO2 and HCHO SCD uncertainties increased by 7–9 and 11–15 % year−1 respectively up until September 2009, when heating of the instrument markedly reduced further throughput loss, stabilising the degradation of SCD uncertainty to −1 for 2009–2015. Our work suggests that the NO2 SCD uncertainty largely consists of a random component ( ∼ 65 % of the total uncertainty) as a result of the propagation of measurement noise but also of a substantial systematic component ( ∼ 35 % of the total uncertainty) mainly from stripe effects. Averaging over multiple pixels in space and/or time can significantly reduce the SCD uncertainties. This suggests that trend detection in OMI, GOME-2 NO2, and HCHO time series is not limited by the spectral fitting but rather by the adequacy of assumptions on the atmospheric state in the later air mass factor (AMF) calculation step.

On board the Copernicus Sentinel-5 Precursor (S5P) platform, the TROPOspheric Monitoring Instrument (TROPOMI) is a double-channel, nadir-viewing grating spectrometer measuring solar back-scattered earthshine radiances in the ultraviolet, visible, near-infrared, and shortwave infrared with global daily coverage. In the ultraviolet range, its spectral resolution and radiometric performance are equivalent to those of its predecessor OMI, but its horizontal resolution at true nadir is improved by an order of magnitude. This paper introduces the formaldehyde (HCHO) tropospheric vertical column retrieval algorithm implemented in the S5P operational processor and comprehensively describes its various retrieval steps. Furthermore, algorithmic improvements developed in the framework of the EU FP7-project QA4ECV are described for future updates of the processor. Detailed error estimates are discussed in the light of Copernicus user requirements and needs for validation are highlighted. Finally, verification results based on the application of the algorithm to OMI measurements are presented, demonstrating the performances expected for TROPOMI.

We present time series of the global distribution of water vapor columns over more than 2 decades based on measurements from the satellite instruments GOME, SCIAMACHY, and GOME-2 in the red spectral range. A particular focus is the consistency amongst the different sensors to avoid jumps from one instrument to another. This is reached by applying robust and simple retrieval settings consistently. Potentially systematic effects due to differences in ground pixel size are avoided by merging SCIAMACHY and GOME-2 observations to GOME spatial resolution, which also allows for a consistent treatment of cloud effects. In addition, the GOME-2 swath is reduced to that of GOME and SCIAMACHY to have consistent viewing geometries.Remaining systematic differences between the different sensors are investigated during overlap periods and are corrected for in the homogenized time series. The resulting Climate product v2.2 (https://doi.org/10.1594/WDCC/GOME-EVL_water_vapor_clim_v2.2) allows the study of the temporal evolution of water vapor over the last 20 years on a global scale.

Experience of differential atmospheric absorption spectroscopy (DOAS) shows that a spectral shift between measurement spectra and reference spectra is frequently required in order to achieve optimal fit results, while the straightforward calculation of the optical density proves inferior. The shift is often attributed to temporal instabilities of the instrument but implicitly solved the problem of the tilt effect discussed/explained in this paper. Spectral positions of Fraunhofer and molecular absorption lines are systematically shifted for different measurement geometries due to an overall slope – or tilt – of the intensity spectrum. The phenomenon has become known as the tilt effect for limb satellite observations, where it is corrected for in a first-order approximation, whereas the remaining community is less aware of its cause and consequences. It is caused by the measurement process, because atmospheric absorption and convolution in the spectrometer do not commute. Highly resolved spectral structures in the spectrum will first be modified by absorption and scattering processes in the atmosphere before they are recorded with a spectrometer, which convolves them with a specific instrument function. In the DOAS spectral evaluation process, however, the polynomial (or other function used for this purpose) accounting for broadband absorption is applied after the convolution is performed. In this paper, we derive that changing the order of the two modifications of the spectra leads to different results. Assuming typical geometries for the observations of scattered sunlight and a spectral resolution of 0.6 nm, this effect can be interpreted as a spectral shift of up to 1.5 pm, which is confirmed in the actual analysis of the ground-based measurements of scattered sunlight as well as in numerical radiative transfer simulations. If no spectral shift is allowed by the fitting routine, residual structures of up to 2.5 × 10−3 peak-to-peak are observed. Thus, this effect needs to be considered for DOAS applications aiming at an rms of the residual of 10−3 and below.

In order to promote the development of the passive DOAS technique the Multi Axis DOAS – Comparison campaign for Aerosols and Trace gases (MAD-CAT) was held at the Max Planck Institute for Chemistry in Mainz, Germany, from June to October 2013. Here, we systematically compare the differential slant column densities (dSCDs) of nitrous acid (HONO) derived from measurements of seven different instruments. We also compare the tropospheric difference of SCDs (delta SCD) of HONO, namely the difference of the SCDs for the non-zenith observations and the zenith observation of the same elevation sequence. Different research groups analysed the spectra from their own instruments using their individual fit software. All the fit errors of HONO dSCDs from the instruments with cooled large-size detectors are mostly in the range of 0.1 to 0.3 × 1015 molecules cm−2 for an integration time of 1 min. The fit error for the mini MAX-DOAS is around 0.7 × 1015 molecules cm−2. Although the HONO delta SCDs are normally smaller than 6 × 1015 molecules cm−2, consistent time series of HONO delta SCDs are retrieved from the measurements of different instruments. Both fits with a sequential Fraunhofer reference spectrum (FRS) and a daily noon FRS lead to similar consistency. Apart from the mini-MAX-DOAS, the systematic absolute differences of HONO delta SCDs between the instruments are smaller than 0.63 × 1015 molecules cm−2. The correlation coefficients are higher than 0.7 and the slopes of linear regressions deviate from unity by less than 16 % for the elevation angle of 1°. The correlations decrease with an increase in elevation angle. All the participants also analysed synthetic spectra using the same baseline DOAS settings to evaluate the systematic errors of HONO results from their respective fit programs. In general the errors are smaller than 0.3 × 1015 molecules cm−2, which is about half of the systematic difference between the real measurements.The differences of HONO delta SCDs retrieved in the selected three spectral ranges 335–361, 335–373 and 335–390 nm are considerable (up to 0.57 × 1015 molecules cm−2) for both real measurements and synthetic spectra. We performed sensitivity studies to quantify the dominant systematic error sources and to find a recommended DOAS setting in the three spectral ranges. The results show that water vapour absorption, temperature and wavelength dependence of O4 absorption, temperature dependence of Ring spectrum, and polynomial and intensity offset correction all together dominate the systematic errors. We recommend a fit range of 335–373 nm for HONO retrievals. In such fit range the overall systematic uncertainty is about 0.87 × 1015 molecules cm−2, much smaller than those in the other two ranges. The typical random uncertainty is estimated to be about 0.16 × 1015 molecules cm−2, which is only 25 % of the total systematic uncertainty for most of the instruments in the MAD-CAT campaign. In summary for most of the MAX-DOAS instruments for elevation angle below 5°, half daytime measurements (usually in the morning) of HONO delta SCD can be over the detection limit of 0.2 × 1015 molecules cm−2 with an uncertainty of ∼ 0.9 × 1015 molecules cm−2.

Satellite nitrogen dioxide (NO2) observations have been widely used to evaluate emission changes. To determine trends in nitrogen oxides (NOx) emission over China, we used a method independent of chemical transport models to quantify the NOx emissions from 48 cities and seven power plants over China, on the basis of Ozone Monitoring Instrument (OMI) NO2 observations from 2005 to 2015. We found that NOx emissions over 48 Chinese cities increased by 52 % from 2005 to 2011 and decreased by 21 % from 2011 to 2015. The decrease since 2011 could be mainly attributed to emission control measures in power sector; while cities with different dominant emission sources (i.e., power, industrial, and transportation sectors) showed variable emission decline timelines that corresponded to the schedules for emission control in different sectors. The time series of the derived NOx emissions was consistent with the bottom-up emission inventories for all power plants (r = 0. 8 on average), but not for some cities (r = 0. 4 on average). The lack of consistency observed for cities was most probably due to the high uncertainty of bottom-up urban emissions used in this study, which were derived from downscaling the regional-based emission data to city level by using spatial distribution proxies.

We determined NOx emissions from Paris in summer 2009 and winter 2009/2010 by applying the closed integral method (CIM) to a large set of car multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements performed within the framework of the MEGAPOLI project (http://megapoli.dmi.dk/). MAX-DOAS measurements of the tropospheric NO2 vertical column density (VCD) were performed in large circles around Paris. From the combination of the observed NO2 VCDs with wind fields, the NO2 influx into and the outflux from the encircled area was determined. The difference between the influx and outflux represents the total emission. Compared to previous applications of the CIM, the large number of measurements during the MEGAPOLI campaign allowed the investigation of important aspects of the CIM. In particular, the applicability of the CIM under various atmospheric conditions could be tested. Another important advantage of the measurements during MEGAPOLI is that simultaneous atmospheric model simulations with a high spatial resolution (3 × 3 km2) are available for all days. Based on these model data, it was possible to test the consistency of the CIM and to derive information about favourable or non-favourable conditions for the application of the CIM. We found that in most situations the uncertainties and the variability in the wind data dominate the total error budget, which typically ranges between 30 and 50 %. Also, measurement gaps and uncertainties in the partitioning ratio between NO and NO2 are important error sources. Based on a consistency check, we deduced a set of criteria on whether measurement conditions are suitable or not for the application of the CIM. We also developed a method for the calculation of the total error budget of the derived NOx emissions. Typical errors are between ±30 and ±50 % for individual days (with one full circle around Paris). From the application of the CIM to car MAX-DOAS observations we derive daily average NOx emissions for Paris of 4.0 × 1025 molec s−1 for summer and of 6.9 × 1025 molec s−1 in winter. These values are a factor of about 1.4 and 2.0 larger than the corresponding emissions derived from the application of the CIM to the model data, using the Toegepast Natuurwetenschappelijk Onderzoek (TNO) MEGAPOLI emission inventory, in summer and winter, respectively. Similar ratios (1.5 and 2.3 for summer and winter, respectively) were found for the comparison with the Monitoring Atmospheric composition and climate III (MACC-III) emission inventory. The highest NOx emissions were found during some cold days in February. Enhanced domestic heating and a reduced conversion efficiency of catalytic converters might contribute to these enhanced NOx emissions.

Tropospheric vertical column densities (VCDs) of NO2, SO2 and HCHO derived from the Ozone Monitoring Instrument (OMI) on AURA and the Global Ozone Monitoring Experiment 2 aboard METOP-A (GOME-2A) and METOP-B (GOME-2B) are widely used to characterize the global distributions, trends and dominating sources of these trace gases. They are also useful for the comparison with chemical transport models (CTMs). We use tropospheric VCDs and vertical profiles of NO2, SO2 and HCHO derived from MAX-DOAS measurements from 2011 to 2014 in Wuxi, China, to validate the corresponding products (daily and bi-monthly-averaged data) derived from OMI and GOME-2A/B by different scientific teams. Prior to the comparison, the spatial and temporal coincidence criteria for MAX-DOAS and satellite data are determined by a sensitivity study using different spatial and temporal averaging conditions. Cloud effects on both MAX-DOAS and satellite observations are also investigated. Our results indicate that the discrepancies between satellite and MAX-DOAS results increase with increasing effective cloud fraction and are dominated by the effects of clouds on the satellite products. In comparison with MAX-DOAS, we found a systematic underestimation of all SO2 (40 to 57 %) and HCHO products (about 20 %), and an overestimation of the GOME-2A/B NO2 products (about 30 %), but good consistency with the DOMINO version 2 NO2 product. To better understand the reasons for these differences, we evaluated the a priori profile shapes used in the OMI retrievals (derived from CTM) by comparison with those derived from the MAX-DOAS observations. Significant differences are found for the SO2 and HCHO profile shapes derived from the IMAGES model, whereas on average good agreement is found for the NO2 profile shapes derived from the TM4 model. We also applied the MAX-DOAS profile shapes to the satellite retrievals and found that these modified satellite VCDs agree better with the MAX-DOAS VCDs than the VCDs from the original data sets by up to 10, 47 and 35 % for NO2, SO2 and HCHO, respectively. Furthermore, we investigated the effect of aerosols on the satellite retrievals. For OMI observations of NO2, a systematic underestimation is found for large AOD, which is mainly attributed to effect of the aerosols on the cloud retrieval and the subsequent application of a cloud correction scheme (implicit aerosol correction). In contrast, the effect of aerosols on the clear-sky air mass factor (explicit aerosol correction) has a smaller effect. For SO2 and HCHO observations selected in the same way, no clear aerosol effect is found, probably because for the considered data sets no cloud correction is applied (and also because of the larger scatter). From our findings we conclude that for satellite observations with cloud top pressure (CTP) > 900 hPa and effective cloud fraction (eCF) 2 VCDs can be considerably overestimated if results from different sensors and/or retrievals (e.g. OMI and GOME-2) are used, whereas fewer deviations for HCHO and SO2 VCDs are found.

Knowledge of the field of view (FOV) of a remote sensing instrument is particularly important when interpreting their data and merging them with other spatially referenced data. Especially for instruments in space, information on the actual FOV, which may change during operation, may be difficult to obtain. Also, the FOV of ground-based devices may change during transportation to the field site, where appropriate equipment for the FOV determination may be unavailable. This paper presents an independent, simple and robust method to retrieve the FOV of an instrument during operation, i.e. the two-dimensional sensitivity distribution, sampled on a discrete grid. The method relies on correlated measurements featuring a significantly higher spatial resolution, e.g. by an imaging instrument accompanying a spectrometer. The method was applied to two satellite instruments, GOME-2 and OMI, and a ground-based differential optical absorption spectroscopy (DOAS) instrument integrated in an SO2 camera. For GOME-2, quadrangular FOVs could be retrieved, which almost perfectly match the provided FOV edges after applying a correction for spatial aliasing inherent to GOME-type instruments. More complex sensitivity distributions were found at certain scanner angles, which are probably caused by degradation of the moving parts within the instrument. For OMI, which does not feature any moving parts, retrieved sensitivity distributions were much smoother compared to GOME-2. A 2-D super-Gaussian with six parameters was found to be an appropriate model to describe the retrieved OMI FOV. The comparison with operationally provided FOV dimensions revealed small differences, which could be mostly explained by the limitations of our IFR implementation. For the ground-based DOAS instrument, the FOV retrieved using SO2-camera data was slightly smaller than the flat-disc distribution, which is assumed by the state-of-the-art correlation technique. Differences between both methods may be attributed to spatial inhomogeneities. In general, our results confirm the already deduced FOV distributions of OMI, GOME-2, and the ground-based DOAS. It is certainly applicable for degradation monitoring and verification exercises. For satellite instruments, the gained information is expected to increase the accuracy of combined products, where measurements of different instruments are integrated, e.g. mapping of high-resolution cloud information, incorporation of surface climatologies. For the SO2-camera community, the method presents a new and efficient tool to monitor the DOAS FOV in the field.

Air mass factor (AMF) calculation is the largest source of uncertainty in NO2 and HCHO satellite retrievals in situations with enhanced trace gas concentrations in the lower troposphere. Structural uncertainty arises when different retrieval methodologies are applied within the scientific community to the same satellite observations. Here, we address the issue of AMF structural uncertainty via a detailed comparison of AMF calculation methods that are structurally different between seven retrieval groups for measurements from the Ozone Monitoring Instrument (OMI). We estimate the escalation of structural uncertainty in every sub-step of the AMF calculation process. This goes beyond the algorithm uncertainty estimates provided in state-of-the-art retrievals, which address the theoretical propagation of uncertainties for one particular retrieval algorithm only. We find that top-of-atmosphere reflectances simulated by four radiative transfer models (RTMs) (DAK, McArtim, SCIATRAN and VLIDORT) agree within 1.5 %. We find that different retrieval groups agree well in the calculations of altitude resolved AMFs from different RTMs (to within 3 %), and in the tropospheric AMFs (to within 6 %) as long as identical ancillary data (surface albedo, terrain height, cloud parameters and trace gas profile) and cloud and aerosol correction procedures are being used. Structural uncertainty increases sharply when retrieval groups use their preference for ancillary data, cloud and aerosol correction. On average, we estimate the AMF structural uncertainty to be 42 % over polluted regions and 31 % over unpolluted regions, mostly driven by substantial differences in the a priori trace gas profiles, surface albedo and cloud parameters. Sensitivity studies for one particular algorithm indicate that different cloud correction approaches result in substantial AMF differences in polluted conditions (5 to 40 % depending on cloud fraction and cloud pressure, and 11 % on average) even for low cloud fractions (2 and HCHO retrievals.

We characterize the temporal variation and vertical distribution of nitrogen dioxide (NO2), sulfur dioxide (SO2), formaldehyde (HCHO) and aerosol extinction based on long-term multi-axis differential optical absorption spectroscopy (MAX-DOAS) observations from May 2011 to November 2014 in Wuxi, China. A new inversion algorithm (PriAM) is implemented to retrieve profiles of the trace gases (TGs) and aerosol extinction (AE) from the UV spectra of scattered sunlight recorded by the MAX-DOAS instrument. We investigated two important aspects of the retrieval process. We found that the systematic seasonal variation of temperature and pressure (which is regularly observed in Wuxi) can lead to a systematic bias of the retrieved aerosol profiles (e.g. up to 20 % for the AOD) if it is not explicitly considered. In this study we take this effect into account for the first time. We also investigated in detail the reason for the differences of tropospheric vertical column densities derived from either the geometric approximation or by the integration of the retrieved profiles, which were reported by earlier studies. We found that these differences are almost entirely caused by the limitations of the geometric approximation (especially for high aerosol loads). The results retrieved from the MAX-DOAS observations are compared with independent techniques not only under cloud-free sky conditions, but also under various cloud scenarios. Under most cloudy conditions (except fog and optically thick clouds), the trace gas results still show good agreements. In contrast, for the aerosol results, only near-surface AE could be still well retrieved under cloudy situations. After applying a quality control procedure, the MAX-DOAS data are used to characterize the seasonal, diurnal and weekly variations of NO2, SO2, HCHO and aerosols. A regular seasonality of the three trace gases is found, but not for aerosols. Similar annual variations of the profiles of the trace gases appear in different years. Only NO2 shows a significant seasonality of the diurnal variations. Considerable amplitudes of weekly cycles occur for NO2 and SO2, but not for HCHO and aerosols. The TGs and aerosols show good correlations, especially for HCHO in winter. More pronounced wind direction dependencies, especially for the near-surface concentrations, are found for the trace gases than for the aerosols, which implies that the local emissions from nearby industrial areas (including traffic emissions) dominate the local pollution, while long-distance transport might also considerably contribute to the local aerosol levels.

The instrumental spectral response function (ISRF) is a key quantity in DOAS analysis, as it is needed for wavelength calibration and for the convolution of trace gas cross sections to instrumental resolution. While it can generally be measured using monochromatic stimuli, it is often parameterized in order to merge different calibration measurements and to plainly account for its wavelength dependency. For some instruments, the ISRF can be described appropriately by a Gaussian function, while for others, dedicated, complex parameterizations with several parameters have been developed.Here we propose to parameterize the ISRF as a super-Gaussian, which can reproduce a variety of shapes, from point-hat to boxcar shape, by just adding one parameter to the classical Gaussian. The super-Gaussian turned out to describe the ISRF of various DOAS instruments well, including the satellite instruments GOME-2, OMI, and TROPOMI.In addition, the super-Gaussian allows for a straightforward parameterization of the effect of ISRF changes, which can occur on long-term scales as well as, for example, during one satellite orbit and impair the spectral analysis if ignored. In order to account for such changes, spectral structures are derived from the derivatives of the super-Gaussian, which are afterwards just scaled during spectral calibration or DOAS analysis. This approach significantly improves the fit quality compared to setups with fixed ISRF, without drawbacks on computation time due to the applied linearization. In addition, the wavelength dependency of the ISRF can be accounted for by accordingly derived spectral structures in an easy, fast, and robust way.

Water vapour is known to absorb radiation from the microwave region to the blue part of the visible spectrum with decreasing efficiency. Ab initio approaches to model individual absorption lines of the gaseous water molecule predict absorption lines up to its dissociation limit at 243 nm.We present the first evidence of water vapour absorption near 363 nm from field measurements using data from multi-axis differential optical absorption spectroscopy (MAX-DOAS) and long-path (LP)-DOAS measurements. The identification of the absorptions was based on the recent POKAZATEL line list by Polyansky et al. (2017). For MAX-DOAS measurements, we observed absorption by water vapour in an absorption band around 363 nm with optical depths of up to 2 × 10−3. The retrieved column densities from 2 months of measurement data and more than 2000 individual observations at different latitudes correlate well with simultaneously measured well-established water vapour absorptions in the blue spectral range from 452 to 499 nm (R2 = 0.89), but the line intensities at around 363 nm are underestimated by a factor of 2.6 ± 0.5 by the ab initio model. At a spectral resolution of 0.5 nm, we derive a maximum cross section value of 2.7 × 10−27 cm2 molec−1 at 362.3 nm. The results were independent of the used literature absorption cross section of the O4 absorption, which overlays this water vapour absorption band. Also water vapour absorption around 376 nm was identified. Below 360 nm no water vapour absorption above 1.4 × 10−26 cm2 molec−1 was observed. The newly found absorption can have a significant impact on the spectral retrievals of absorbing trace-gas species in the spectral range around 363 nm. Its effect on the spectral analysis of O4, HONO and OClO is discussed.

The Rann of Kutch (India and Pakistan) is one of the largest salt deserts in the world. Being a so-called "seasonal salt marsh", it is regularly flooded during the Indian summer monsoon. We present 10 years of bromine monoxide (BrO) satellite observations by the Ozone Monitoring Instrument (OMI) over the Great and Little Rann of Kutch. OMI spectra were analysed using Differential Optical Absorption Spectroscopy (DOAS) and revealed recurring high BrO vertical column densities (VCDs) of up to 1.4 × 1014 molec cm−2 during April/May, but no significantly enhanced column densities during the monsoon season (June–September). In the following winter months, the BrO VCDs are again slightly enhanced while the salty surface dries up. We investigate a possible correlation of enhanced reactive bromine concentrations with different meteorological parameters and find a strong relationship between incident UV radiation and the total BrO abundance. In contrast, the second Global Ozone Monitoring Instrument (GOME-2) shows about 4 times lower BrO VCDs over the Rann of Kutch than found by OMI and no clear seasonal cycle is observed. One reason for this finding might be the earlier local overpass time of GOME-2 compared to OMI (around 09:30 vs. 13:30 LT), as the ambient conditions significantly differ for both satellite instruments at the time of the measurements. Further possible reasons are discussed and mainly attributed to instrumental issues. OMI additionally confirms the presence of enhanced BrO concentrations over the Dead Sea valley (Israel/Jordan), as suggested by former ground-based observations. The measurements indicate that the Rann of Kutch salt marsh is probably one of the strongest natural point sources of reactive bromine compounds outside the polar regions and is therefore supposed to have a significant impact on local and regional ozone chemistry.

A method is developed for the calibration of the colour index (CI) and the O4 absorption derived from differential optical absorption spectroscopy (DOAS) measurements of scattered sunlight. The method is based on the comparison of measurements and radiative transfer simulations for well-defined atmospheric conditions and viewing geometries. Calibrated measurements of the CI and the O4 absorption are important for the detection and classification of clouds from MAX-DOAS observations. Such information is needed for the identification and correction of the cloud influence on Multi AXis (MAX-)DOAS profile inversion results, but might be also be of interest on their own, e.g. for meteorological applications. The calibration algorithm was successfully applied to measurements at two locations: Cabauw in the Netherlands and Wuxi in China. We used CI and O4 observations calibrated by the new method as input for our recently developed cloud classification scheme and also adapted the corresponding threshold values accordingly. For the observations at Cabauw, good agreement is found with the results of the original algorithm. Together with the calibration procedure of the CI and O4 absorption, the cloud classification scheme, which has been tuned to specific locations/conditions so far, can now be applied consistently to MAX-DOAS measurements at different locations. In addition to the new threshold values, further improvements were introduced to the cloud classification algorithm, namely a better description of the SZA (solar zenith angle) dependence of the threshold values and a new set of wavelengths for the determination of the CI. We also indicate specific areas for future research to further improve the cloud classification scheme.