Abstract. The results of a comparison exercise of radiative transfer models (RTM) of various international research groups for Multiple AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) viewing geometry are presented. Besides the assessment of the agreement between the different models, a second focus of the comparison was the systematic investigation of the sensitivity of the MAX-DOAS technique under various viewing geometries and aerosol conditions. In contrast to previous comparison exercises, box-air-mass-factors (box-AMFs) for different atmospheric height layers were modelled, which describe the sensitivity of the measurements as a function of altitude. In addition, radiances were calculated allowing the identification of potential errors, which might be overlooked if only AMFs are compared. Accurate modelling of radiances is also a prerequisite for the correct interpretation of satellite observations, for which the received radiance can strongly vary across the large ground pixels, and might be also important for the retrieval of aerosol properties as a future applicationCorrespondence to: T. Wagner (thomas.wagner@iup.uni-heidelberg.de) of MAX-DOAS. The comparison exercises included different wavelengths and atmospheric scenarios (with and without aerosols). The strong and systematic influence of aerosol scattering indicates that from MAX-DOAS observations also information on atmospheric aerosols can be retrieved. During the various iterations of the exercises, the results from all models showed a substantial convergence, and the final data sets agreed for most cases within about 5%. Larger deviations were found for cases with low atmospheric optical depth, for which the photon path lengths along the line of sight of the instrument can become very large. The differences occurred between models including full spherical geometry and those using only plane parallel approximation indicating that the correct treatment of the Earth's sphericity becomes indispensable. The modelled box-AMFs constitute an universal data base for the calculation of arbitrary (total) AMFs by simple convolution with a given trace gas concentration profile. Together with the modelled radiances and the specified settings for the various exercises, they can serve as test cases for future RTM developments.Published by Copernicus GmbH on behalf of the European Geosciences Union.

Abstract. We conducted long-term network observations using standardized Multi-Axis Differential optical absorption spectroscopy (MAX-DOAS) instruments in Russia and ASia (MADRAS) from 2007 onwards and made the first synthetic data analysis. At seven locations (Cape Hedo, Fukue and Yokosuka in Japan, Hefei in China, Gwangju in Korea, and Tomsk and Zvenigorod in Russia) with different levels of pollution, we obtained 80 927 retrievals of tropospheric NO2 vertical column density (TropoNO2VCD) and aerosol optical depth (AOD). In the technique, the optimal estimation of the TropoNO2VCD and its profile was performed using aerosol information derived from O4 absorbances simultaneously observed at 460–490 nm. This large data set was used to analyze NO2 climatology systematically, including temporal variations from the seasonal to the diurnal scale. The results were compared with Ozone Monitoring Instrument (OMI) satellite observations and global model simulations. Two NO2 retrievals of OMI satellite data (NASA ver. 2.1 and Dutch OMI NO2 (DOMINO) ver. 2.0) generally showed close correlations with those derived from MAX-DOAS observations, but had low biases of up to ~50%. The bias was distinct when NO2 was abundantly present near the surface and when the AOD was high, suggesting a possibility of incomplete accounting of NO2 near the surface under relatively high aerosol conditions for the satellite observations. Except for constant biases, the satellite observations showed nearly perfect seasonal agreement with MAX-DOAS observations, suggesting that the analysis of seasonal features of the satellite data were robust. Weekend reduction in the TropoNO2VCD found at Yokosuka and Gwangju was absent at Hefei, implying that the major sources had different weekly variation patterns. While the TropoNO2VCD generally decreased during the midday hours, it increased exceptionally at urban/suburban locations (Yokosuka, Gwangju, and Hefei) during winter. A global chemical transport model, MIROC-ESM-CHEM (Model for Interdisciplinary Research on Climate–Earth System Model–Chemistry), was validated for the first time with respect to background NO2 column densities during summer at Cape Hedo and Fukue in the clean marine atmosphere.

[1] This study compares the limb scattered radiances calculated by six radiative transfer models for a variety of viewing conditions. Atmospheres that include molecular scattering, aerosol scattering, and ozone absorption are considered. All models treat single scattering accurately in full spherical geometry. Two ''approximate spherical'' models (CDI and LIMBTRAN) rely on the plane-parallel atmosphere approximation to calculate the diffuse radiance field; the remaining four ''spherical'' models (Siro, MCC++, GSLS, and CDIPI) treat multiple scattering in a spherical atmosphere. Only three of the models (Siro, MCC++, and GSLS) have vector treatment with polarization. A brief comparison of vector radiances with the limb scattered radiances measured by the SOLSE and LORE instruments demonstrates agreement usually within 15% and always within 30%. The inclusion of polarization appears to have little effect on the level of agreement among the models (which agree to within 2% for this sample case). A more general comparison among calculated scalar radiances follows, including four solar zenith angles (20°, 60°, 80°, and 90°), three relative azimuth angles (20°, 90°, and 160°), and two surface albedos (0 and 0.95). The single scattered radiances agree to within 1% for almost every case. Comparisons of the total radiance show larger differences, with 2-4% spread among the results of the spherical models. The approximate spherical models show a positive radiance difference relative to the other models that increases with tangent height, reaching as much as 8% at 60 km. The rule used to divide the model atmosphere into discrete layers is shown to affect the calculated radiance, causing a height-dependent difference of up to 1% for 1 km layer thickness.

Context. Twilight studies have proved to be important tools for analyzing the atmospheric structure with interesting consequences on the characterization of astronomical sites. Active discussions of this topic have started again recently in connection with the evaluation of Dome C, Antarctica, as a potential astronomical site and several site-testing experiments, including twilight brightness measurements, are being prepared. Aims. The present work provides the first absolute photometric measurements of twilight sky brightness for ESO-Paranal (Chile). They are meant both as a contribution to the site monitoring and as reference values in the analysis of other sites, including Dome C. Methods. The UBVRI twilight sky brightness was estimated on more than 2000 FORS1 archival images, which include both flats and standard star observations taken in twilight and covering a Sun zenith distance range 94• -112• . Results. The comparison with a low-altitude site shows that Paranal V twilight sky brightness is about 30% lower, implying that some fraction of multiple scattering has to take place at an altitude of a few km above sea level.