Abstract. In this paper, we present long-term observations of atmospheric nitrogen dioxide (NO2) and formaldehyde (HCHO) in Nanjing using a Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) instrument. Ground-based MAX-DOAS measurements were performed from April 2013 to February 2017. The MAX-DOAS measurements of NO2 and HCHO vertical column densities (VCDs) are used to validate ozone monitoring instrument (OMI) satellite observations over Nanjing. The comparison shows that the OMI observations of NO2 correlate well with the MAX-DOAS data with Pearson correlation coefficient (R) of 0.91. However, OMI observations are on average a factor of 3 lower than the MAX-DOAS measurements. Replacing the a priori NO2 profiles by the MAX-DOAS profiles in the OMI NO2 VCD retrieval would increase the OMI NO2 VCDs by ∼30 % with correlation nearly unchanged. The comparison result of MAX-DOAS and OMI observations of HCHO VCD shows a good agreement with R of 0.75 and the slope of the regression line is 0.99. An age-weighted backward-propagation approach is applied to the MAX-DOAS measurements of NO2 and HCHO to reconstruct the spatial distribution of NO2 and HCHO over the Yangtze River Delta during summer and winter time. The reconstructed NO2 fields show a distinct agreement with OMI satellite observations. However, due to the short atmospheric lifetime of HCHO, the backward-propagated HCHO data do not show a strong spatial correlation with the OMI HCHO observations. The result shows that the MAX-DOAS measurements are sensitive to the air pollution transportation in the Yangtze River Delta, indicating the air quality in Nanjing is significantly influenced by regional transportation of air pollutants. The MAX-DOAS data are also used to evaluate the effectiveness of air pollution control measures implemented during the Youth Olympic Games 2014. The MAX-DOAS data show a significant reduction of ambient aerosol, NO2 and HCHO (30 %–50 %) during the Youth Olympic Games. Our results provide a better understanding of the transportation and sources of pollutants over the Yangtze River Delta as well as the effect of emission control measures during large international events, which are important for the future design of air pollution control policies.

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

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

Abstract. Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) is a well-established ground-based measurement technique for the detection of aerosols and trace gases particularly in the boundary layer and the lower troposphere: ultraviolet- and visible radiation spectra of skylight are analysed to obtain information on different atmospheric parameters, integrated over the light path from space to the instrument. An appropriate set of spectra recorded under different viewing geometries ("Multi-Axis") allows retrieval of tropospheric aerosol and trace gas vertical distributions by applying numerical inversion methods. 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 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, sun photometer and others). 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. The profiles were found to be in good qualitative agreement: most participants obtained the same features in the retrieved vertical trace gas and aerosol distributions, however sometimes at different altitudes and of different intensity. Under clear sky conditions, the root-mean-square differences of aerosol optical thicknesses, trace gas (NO2, HCHO) vertical columns and surface concentrations among the results of individual participants vary between 0.01–0.1, (1.5–15) x 1014 molec cm-2 and (0.3–8) x 1010 molec cm-3, respectively. For the comparison against supporting observations, these values increase to 0.02–0.2, (11–55) x 1014 molec cm-2 and (0.8–9) x 1010 molec cm-3. It is likely that a large part of this increase is caused by imperfect spatio-temporal overlap of the different observations. In contrast to what is often assumed, the MAX-DOAS vertically integrated extinction profiles and the sun photometer total aerosol optical thickness were found to not necessarily being comparable quantities, unless information on the real aerosol vertical distribution is available to account for the low sensitivity of MAX-DOAS observations at higher altitudes.

Abstract. In September 2016, 36 spectrometers from 24 institutes measured a number of key atmospheric pollutants for a period of 17 days 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 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, to discuss the performance of the various types of instruments and 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 dimer (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. 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 an unprecedented 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 reference, 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 measurement performance of all the participating instruments for the MAX-DOAS and zenith-sky DOAS techniques.

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

We present a new aerosol extinction profile retrieval algorithm for multi-axis differential optical absorption spectrometer (MAX-DOAS) measurements at high-altitude sites. The algorithm is based on the lookup table method. It is applied to retrieve aerosol extinction profiles from the longterm MAX-DOAS measurements (February 2012 to February 2016) at the Environmental Research Station Schneefernerhaus (UFS), Germany (47.417 • N, 10.980 • E), which is located near the summit of Zugspitze at an altitude of 2650 m. The lookup table consists of simulated O 4 differential slant column densities (DSCDs) corresponding to numerous possible aerosol extinction profiles. The sensitivities of O 4 absorption to several parameters were investigated for the design and parameterization of the lookup table. In the retrieval, simulated O 4 DSCDs for each possible profile are derived by interpolating the lookup table to the observation geometries. The cost functions are calculated for each aerosol profile in the lookup table based on the simulated O 4 DSCDs, the O 4 DSCD observations, and the measurement uncertainties. Valid profiles are selected from all the possible profiles according to the cost function, and the optimal solution is defined as the weighted mean of all the valid profiles. A comprehensive error analysis is performed to better estimate the total uncertainty. Based on the assumption that the lookup table covers all possible profiles under clear-sky conditions, we determined a set of O 4 DSCD scaling factors for differ-ent elevation angles and wavelengths. The profiles retrieved from synthetic measurement data can reproduce the synthetic profile. The results also show that the retrieval is insensitive to measurement noise, indicating the retrieval is robust and stable. The aerosol optical depths (AODs) retrieved from the long-term measurements were compared to coinciding and co-located sun photometer observations. High correlation coefficients (R) of 0.733 and 0.798 are found for measurements at 360 and 477 nm, respectively. However, especially in summer, the sun photometer AODs are systematically higher than the MAX-DOAS retrievals by a factor of ∼ 2. The discrepancy might be related to the limited measurement range of the MAX-DOAS and is probably also related to the decreased sensitivity of the MAX-DOAS measurements at higher altitudes. The MAX-DOAS measurements indicate the aerosol extinction decreases with increasing altitude during all seasons, which agrees with the co-located ceilometer measurements. Our results also show maximum AOD and maximum Ångström exponent in summer, which is consistent with observations at an AERONET station located ∼ 43 km from the UFS.

<p><strong>Abstract.</strong> In this paper, we present long term observations of atmospheric nitrogen dioxide (NO<sub>2</sub>) and formaldehyde (HCHO) in Nanjing using a Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) instrument. Ground based MAX-DOAS measurements were performed from April 2013 to February 2017. The MAX-DOAS measurements of NO<sub>2</sub> and HCHO vertical column densities (VCDs) are used to validate OMI satellite observations over Nanjing. The comparison shows that the OMI observations of NO<sub>2</sub> correlate well with the MAX-DOAS data with Pearson correlation coefficient (<i>R</i>) of 0.91. However, OMI observations are on average a factor of 3 lower than the MAX-DOAS measurements. Replacing the a priori NO<sub>2</sub> profiles by the MAX-DOAS profiles in the OMI NO2 VCD retrieval would increase the OMI NO<sub>2</sub> VCDs by ~ 30 % with correlation nearly unchanged. The comparison result of MAX-DOAS and OMI observations of HCHO VCD shows a good agreement with <i>R</i> of 0.75 and the slope of the regression line is 0.99. We developed a new technique to assemble the source contribution map using backward trajectory analysis. The age weighted backward propagation approach is applied to the MAX-DOAS measurements of NO<sub>2</sub> and HCHO to reconstruct the spatial distribution of NO<sub>2</sub> and HCHO over the Yangtze River Delta during summer and winter time. The reconstructed NO<sub>2</sub> fields show a distinct agreement with OMI satellite observations. However, due to the short atmospheric lifetime of HCHO, the backward propagated HCHO data does not show a strong spatial correlation with the OMI HCHO observations. The result shows the MAX-DOAS measurements are sensitive to the air pollution transportation in the Yangtze River Delta, indicating the air quality in Nanjing is significantly influenced by regional transportation of air pollutants. The MAX-DOAS data are also used to evaluate the effectiveness of air pollution control measures implemented during the Youth Olympic Games 2014. The MAX-DOAS data show a significant reduction of ambient aerosol, NO<sub>2</sub> and HCHO (30 %–50 %) during the Youth Olympic Games. Our results provide a better understanding of the transportation and sources of pollutants in over the Yangtze River Delta as well as the effect of emission control measures during large international event, which are important for the future design of air pollution control policies.</p>