ALBERT

Description: 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.

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

Description: 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.

Description: This paper presents the retrieval algorithm for the operational Ozone Monitoring Instrument (OMI) total bromine monoxide (BrO) data product (OMBRO) developed at the Smithsonian Astrophysical Observatory (SAO) and shows comparisons with correlative measurements and retrieval results. The algorithm is based on direct nonlinear least squares fitting of radiances from the spectral range 319.0–347.5 nm. Radiances are modeled from the solar irradiance, attenuated by contributions from BrO and interfering gases, and including rotational Raman scattering, additive and multiplicative closure polynomials, correction for Nyquist undersampling and the average fitting residual spectrum. The retrieval uses albedo- and wavelength-dependent air mass factors (AMFs), which have been pre-computed using a single mostly stratospheric BrO profile. The BrO cross sections are multiplied by the wavelength-dependent AMFs before fitting so that the vertical column densities (VCDs) are retrieved directly. The fitting uncertainties of BrO VCDs typically vary between 4 and 7×1012 molecules cm−2 (∼10 %–20 % of the measured BrO VCDs). Additional fitting uncertainties can be caused by the interferences from O2-O2 and H2CO and their correlation with BrO. AMF uncertainties are estimated to be around 10 % when the single stratospheric-only BrO profile is used. However, under conditions of high tropospheric concentrations, AMF errors due to this assumption of profile can be as high as 50 %. The retrievals agree well with GOME-2 observations at simultaneous nadir overpasses and with ground-based zenith-sky measurements at Harestua, Norway, with mean biases less than -0.22±1.13×1013 and 0.12±0.76×1013 molecules cm−2, respectively. Global distribution and seasonal variation of OMI BrO are generally consistent with previous satellite observations. Finally, we confirm the capacity of OMBRO retrievals to observe enhancements of BrO over the US Great Salt Lake despite the current retrieval setup considering a stratospheric profile in the AMF calculations. OMBRO retrievals also show significant BrO enhancements from the eruption of the Eyjafjallajökull volcano, although the BrO retrievals are affected under high SO2 loading conditions by the sub-optimum choice of SO2 cross sections.

Description: 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.

Description: The spatial and temporal changes in SO2 regimes over China during 2005 to 2016 and their associated driving mechanism are investigated based on a state-of-the-art retrieval dataset. Climatological SO2 exhibits pronounced seasonal and regional variations, with higher loadings in wintertime and two prominent maxima centered in the North China Plain and the Cheng-Yu District. In the last decade, overall SO2 decreasing trends have been reported nationwide, with spatially varying downward rates according to a general rule – the higher the SO2 loading, the more significant the decrease. However, such decline is in fact not monotonic, but instead four distinct temporal regimes can be identified by empirical orthogonal function analysis. After an initial rise at the beginning, SO2 in China undergoes two sharp drops in the periods 2007–2008 and 2014–2016, amid which 5-year moderate rebounding is sustained. Despite spatially coherent behaviors, different mechanisms are tied to North China and South China. In North China, the same four regimes are detected in the time series of emission that is expected to drive the regime of atmospheric SO2, with a percentage of explained variance amounting to 81 %. Out of total emission, those from the industrial sector dominate SO2 variation throughout the whole period, while the role of household emission remains uncertain. In contrast to North China, SO2 emissions in South China exhibit a continuous descending tendency, due to the coordinated cuts of industrial and household emissions. As a result, the role of emissions only makes up about 45 % of the SO2 variation, primarily owing to the decoupled pathways of emission and atmospheric content during 2009 to 2013 when the emissions continue to decline but atmospheric content witnesses a rebound. Unfavorable meteorological conditions, including deficient precipitation, weaker wind speed and increased static stability, outweigh the effect of decreasing emissions and thus give rise to the rebound of SO2 during 2009 to 2013.

Description: 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.

Description: 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 (

Description: The differential optical absorption spectroscopy (DOAS) method is a well-known remote sensing technique that is nowadays widely used for measurements of atmospheric trace gases, creating the need for harmonization and characterization efforts. In this study, an intercomparison exercise of DOAS retrieval codes from 17 international groups is presented, focusing on NO2 slant columns. The study is based on data collected by one instrument during the Multi-Axis DOAS Comparison campaign for Aerosols and Trace gases (MAD-CAT) in Mainz, Germany, in summer 2013. As data from the same instrument are used by all groups, the results are free of biases due to instrumental differences, which is in contrast to previous intercomparison exercises.While in general an excellent correlation of NO2 slant columns between groups of  〉  99.98 % (noon reference fits) and  〉  99.2 % (sequential reference fits) for all elevation angles is found, differences between individual retrievals are as large as 8 % for NO2 slant columns and 100 % for rms residuals in small elevation angles above the horizon.Comprehensive sensitivity studies revealed that absolute slant column differences result predominantly from the choice of the reference spectrum while relative differences originate from the numerical approach for solving the DOAS equation as well as the treatment of the slit function. Furthermore, differences in the implementation of the intensity offset correction were found to produce disagreements for measurements close to sunrise (8–10 % for NO2, 80 % for rms residual). The largest effect of  ≈  8 % difference in NO2 was found to arise from the reference treatment; in particular for fits using a sequential reference. In terms of rms fit residual, the reference treatment has only a minor impact. In contrast, the wavelength calibration as well as the intensity offset correction were found to have the largest impact (up to 80 %) on rms residual while having only a minor impact on retrieved NO2 slant columns.