A multi-axis differential optical absorption spectroscopy aerosol profile retrieval algorithm for high-altitude measurements: application to measurements at Schneefernerhaus (UFS), Germany
- 1Remote Sensing Technology Institute, German Aerospace Center (DLR), Oberpfaffenhofen, Germany
- 2Faculty of Civil, Geo and Environmental Engineering, Technical University of Munich (TUM), Munich, Germany
- 3Max Planck Institute for Chemistry, Mainz, Germany
- 4German Meteorological Office (DWD), Hohenpeißenberg, Germany
- 5Meteorological Institute, Ludwig Maximilian University of Munich (LMU), Munich, Germany
Correspondence: Ka Lok Chan (firstname.lastname@example.org)
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 long-term 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 O4 differential slant column densities (DSCDs) corresponding to numerous possible aerosol extinction profiles. The sensitivities of O4 absorption to several parameters were investigated for the design and parameterization of the lookup table. In the retrieval, simulated O4 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 O4 DSCDs, the O4 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 O4 DSCD scaling factors for different 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.
Atmospheric aerosols play an important role in atmospheric physics and chemistry. They affect the atmospheric radiation budget by absorbing and scattering radiation, as well as providing nuclei for the formation of clouds (Haywood and Boucher, 2000; Bellouin et al., 2005; Li and Kou, 2011; Heald et al., 2014). Aerosols also have significant impacts on global climate change, local air quality, and visibility (Bäumer et al., 2008; Levy et al., 2013; Viana et al., 2014). Moreover, exposure to atmospheric aerosols can be harmful to human health (Valavanidis et al., 2008; Brook et al., 2010; Karanasiou et al., 2012). Besides primary aerosols which are directly introduced into the atmosphere, aerosols can also be secondarily formed through chemical reactions (Hinds, 2012). A significant increasing amount of anthropogenic aerosols and precursors have been released into the atmosphere since the industrial revolution (Liu et al., 1991; Junker and Liousse, 2008), which has become a far-reaching environmental problem in recent years. Aerosols can be long-range transported and hence influence regions far from the sources (Wiegner et al., 2011; Almeida-Silva et al., 2013; Lee et al., 2013; Zhang et al., 2014; Chan and Chan, 2017; Chan, 2017; Chan et al., 2018). The properties and vertical distribution of aerosols vary strongly with time and location. Therefore, it is important to measure the spatial and temporal variations in aerosols for better understanding of the role of aerosols in atmospheric processes. In addition, anthropogenic contribution to atmospheric aerosol load is one of the largest uncertainties in climate forcing assessments. Accurate measurements of aerosol optical properties are necessary for the further assessment of environmental and radiative effects of aerosols (Stocker et al., 2013).
Methodologies for aerosol monitoring are mature and well-established: the backbone is certainly the AERONET network of sun photometers (Holben et al., 1998) providing the spectral aerosol optical depth (AOD) from direct sun observations. They might be complemented by active lidar remote sensing to provide range-resolved information. The latter includes research lidars (e.g., Pappalardo et al., 2014) and networks of ceilometers (e.g., Wiegner et al., 2014; Cazorla et al., 2017). These measurements provide – depending on the complexity of the system – the vertical distribution of the particle backscatter and extinction coefficient at typically one to three wavelengths with a very high vertical resolution on the order of 10 m. However, the uncertainty of the retrieved AOD is larger than that of sun photometers due to the restrictions of the measurement range. In the case of ceilometers, inherent assumptions of the data evaluation further add to the uncertainty. Recently, the potential of a multi-axis differential optical absorption spectrometer (MAX-DOAS) for range-resolved aerosol retrievals was investigated as well (Platt and Stutz, 2008; Wagner et al., 2004; Frieß et al., 2006).
Ground-based multi-axis differential optical absorption spectroscopy is a remote sensing technique for measuring atmospheric aerosols and trace gases. MAX-DOAS instruments measure the spectra of scattered sunlight at several different viewing directions, and information of trace gas absorption along the light paths can be obtained by applying the differential optical absorption spectroscopy (DOAS) method to the ultraviolet–visible (UV–VIS) band. The retrieval of aerosol extinction profiles from MAX-DOAS measurements typically relies on the absorption signal of oxygen collision complex (O4). As the vertical distribution profile of O4 is well-known and stable, it is an ideal indicator of the atmospheric distribution of photon paths. Photon paths of scattered sunlight can be influenced by aerosols and hence change the measured O4 slant columns. Therefore, aerosol vertical extinction profiles can be retrieved by fitting the O4 observations to radiative transfer simulations. Since the experimental setup is relatively simple and inexpensive, MAX-DOAS instruments have been widely used to measure the vertical distribution of atmospheric aerosols and trace gases in the past 2 decades (e.g., Hönninger et al., 2004; Irie et al., 2008, 2011; Li et al., 2010, 2013; Clémer et al., 2010; Frieß et al., 2011; Halla et al., 2011; Vlemmix et al., 2011; Wagner et al., 2011; Ma et al., 2013; Wang et al., 2014a, 2016; Chan et al., 2015, 2017; Jin et al., 2016).
In the retrieval of vertical profile information from MAX-DOAS measurements, the aerosol profile is usually regarded as the state vector (x), and the measured O4 differential slant column densities (DSCDs) of each scanning cycle are regarded as the measurement vector (y). The radiative transfer model used to simulate the O4 DSCDs is regarded as the forward model (F). As the radiative transfer in the atmosphere is nonlinear, the retrieval is a nonlinear problem. Moreover, the retrieval is ill-posed, which means the information contained in the observation is insufficient to determine a unique solution. In many of the other MAX-DOAS studies (e.g., Frieß et al., 2006, 2011; Clémer et al., 2010; Irie et al., 2011; Wang et al., 2014a, 2016; Chan et al., 2017), aerosol profiles are retrieved using the optimal estimation method (OEM) (Rodgers, 2000). The inversion of the aerosol profile is solved iteratively by minimizing the cost function. Vertical profile information can also be retrieved from MAX-DOAS observations using parameterized approaches (e.g., Lee et al., 2009; Li et al., 2010; Vlemmix et al., 2011; Wagner et al., 2011; Sinreich et al., 2013). These methods simplify aerosol profiles as limited parameters, e.g., aerosol optical depth (AOD), layer height, shape parameter (Wagner et al., 2011; Hartl and Wenig, 2013). The optimal solution is usually determined by minimizing the difference between simulations and measurements.
However, as the retrieval is ill-posed and errors exist in both measurements and simulations, the profile with the lowest cost function may not be the one closest to the true profile. Moreover, in the typical OEM-based algorithms, the iteration stops as soon as the cost function is smaller than a certain threshold. Therefore, the retrieved profile is not necessarily the one with the smallest cost function. At high-altitude sites, the aerosol profile retrieval is more challenging, as the O4 concentration and the aerosol load are both much lower than that at low-altitude sites. The vertical gradient of the aerosol extinction is also much smaller and the relative contribution from aerosols above the retrieval height to the total AOD is more significant. As a result, the signal-to-noise ratio (SNR) of high-altitude MAX-DOAS measurements is much lower and hence affects the retrieval quality.
In this paper, we present a new MAX-DOAS aerosol profile retrieval algorithm suitable for high-altitude measurements. It is based on an O4 DSCD lookup table. The lookup table includes simulated O4 DSCDs corresponding to a very large number of aerosol extinction profiles. Our retrieval algorithm is applied to MAX-DOAS observations at the Environmental Research Station Schneefernerhaus (Umweltforschungsstation Schneefernerhaus, UFS). The UFS is located close to the summit of Zugspitze (2962 m above sea level), the highest mountain of Germany, at an altitude of 2650 m. The O4 concentration at Zugspitze is ∼40 % lower compared to sea level. As the measurement site is surrounded by the mountainous area of the Alps and far from polluted areas, the aerosol load is much lower than at low-altitude sites. The annual averaged AOD measured by the sun photometer at the UFS is around 0.1 at 350–500 nm. Moreover, the surface around the UFS is very complex, which complicates the radiative transfer simulation. As a result, the model errors are larger compared to the flat and simple surfaces. In the study, we first analyzed the simulation uncertainty caused by the simplification of topography definition (see Sect. 3.3). Then we studied the sensitivity of O4 absorption to several parameters (see Sect. 3.4 and Appendix B). Based on the results, we designed the O4 DSCD lookup table and the inversion method (see Sect. 3.5 to 3.8). In Sect. 3.9, we present our method for determining the O4 DSCD scaling factors based on the lookup table. Discussions of the retrieved aerosol profiles from the long-term measurements at the UFS are presented in Sect. 4.
2.1 MAX-DOAS measurements
The MAX-DOAS instrument is set up on the platform on the fifth floor of the UFS (47.417∘ N, 10.980∘ E), about 20 m above ground level, which is about 2650 m above sea level. The instrument consists of a scanning telescope, a stepping motor which controls the viewing zenith angle of the telescope, and two spectrometers covering the ultraviolet (UV) and visible (VIS) wavelength bands. Incoming sunlight is redirected by a prism reflector and a quartz fiber bundle to the spectrometers for spectral analysis. The field of view (FOV) of the instrument is ∼0.98∘. Two spectrometers (OMT Instruments, OMT ctf-60) each equipped with a charge-coupled device (CCD) detector are used to measure the spectra of both UV (320–478 nm) and VIS (427–649 nm) wavelength ranges. The full width at half maximum (FWHM) spectral resolutions of the UV and VIS spectrometers are about 1.0 and 0.6 nm, respectively. The scanning direction of the telescope is controlled by the stepping motor.
As the measurement geometry is limited by the topography, the viewing azimuth angle of the telescope was adjusted to due south (180∘) with the lowest elevation angle of 1∘. Each scanning cycle consists of measurements at elevation angles (α) of 90∘ (zenith), 30, 20, 10, 5, 2, and 1∘. A single measurement at each elevation angle lasts for ∼1 min, and a full scanning cycle takes about 10 min. The recorded spectrum of each measurement is the sum of the CCD readouts within ∼1 min. In order to optimize the measurement SNR, avoid saturation, and achieve a constant signal level, the data acquisition software automatically adjusts the exposure time of each readout to make the maximum count close to 70 % of saturation level (65 535 counts). Depending on the intensity of received light, the exposure time of each readout varies from tens of milliseconds to a few seconds. The measurements of the UV and VIS bands are taken by the two spectrometers simultaneously, but their exposure times are adjusted individually. The instrument takes measurements continuously during daytime (solar zenith angle (SZA) <85∘), but during the noon ( solar azimuth angle (SAA) <185∘) and twilight periods ( SZA <92∘), the instrument takes only zenith measurements.
The MAX-DOAS instrument has been running since February 2012. However, the measurement was interrupted between February 2013 and July 2013 due to instrument maintenance. In February 2016, the measurement was interrupted again, and the VIS spectrometer was found to be degraded. In this paper, we present 4 years of MAX-DOAS measurements from February 2012 to February 2016.
2.2 Sun photometer measurements
Next to the MAX-DOAS instrument, a sun photometer is installed at the UFS, which provides measurements of radiances at 12 wavelengths between 340 and 1640 nm with a temporal resolution of 1 s. The instrument was developed at the Meteorological Institute of Ludwig Maximilian University of Munich (LMU) based on a system operated in the framework of the SAMUM campaigns (Toledano et al., 2009, 2011) but with improved electronics and data acquisition developed by Physikalische Messsysteme Ltd. In this study, the AODs derived from sun photometer measurements applying the well-established Rayleigh calibration method were used for the intercomparison with the MAX-DOAS retrieval. For this purpose, AOD measurements at 340 and 380 nm were interpolated to 360 nm while AODs at 477 nm were interpolated from the measurements at 440 and 500 nm. The interpolation followed the Ångström exponent method. Measurements were given as hourly averages. Due to the reduced accuracy under large SZA, only the measurements between 10:00 and 14:00 UTC each day were used. In order to ensure the data quality, only cloud-free conditions and periods of stable aerosol abundance (variability of radiances below 5 % within 1 h) were considered. These requirements reduce the number of available sun photometer measurements considerably. Note that the AOD is often below 0.02 at the relevant wavelengths with an uncertainty on the order of ±0.015 due to calibration errors, Rayleigh correction and radiometric accuracy. As the uncertainty of the AOD measured by the sun photometer is relatively large, the uncertainty of the Ångström exponent would be further amplified. Consequently they are not used in this study.
The aerosol optical properties which are not available from the UFS measurements but required for our MAX-DOAS inversion scheme (single-scattering albedo and phase function) were estimated from the AERONET measurements at Hohenpeißenberg, which is located at an altitude of 980 m and approximately 43 km north of the UFS. The AERONET data were available at 440, 675, 870, and 1020 nm. Therefore, the data at 360 nm were extrapolated, and the data at 477 nm were interpolated. As Hohenpeißenberg and UFS are located at different altitudes, the aerosol optical properties might be slightly different. Therefore, we have analyzed the uncertainties caused by the differences in single-scattering albedo and phase functions through a sensitivity analysis. The results show that the influences of aerosol optical properties are in general less than 3 %; see Appendix B3. Some other MAX-DOAS studies also found that aerosol optical properties show only small impacts on aerosol profile retrieval (e.g., Chan et al., 2019).
2.3 Ceilometer measurements
The UFS is also equipped with a Lufft (previously Jenoptik) ceilometer (model CHM15kx; see Wiegner and Geiß, 2012) operated by the German Weather Service (DWD). Ceilometers are single-wavelength backscatter lidars, and the received signals follow the well-known lidar equation (Wiegner et al., 2014). The CHM15kx is eye-safe and fully automated, which allows unattended 24/7 operation. It can be used to monitor aerosol layers (e.g., volcanic ash; see Schäfer et al., 2011) and validate meteorological and chemistry transport models (see, e.g., Emeis et al., 2011) and is foreseen for model assimilation (e.g., Wang et al., 2014b; Warren et al., 2018; Chan et al., 2018).
The CHM15kx ceilometer is equipped with a diode-pumped Nd:YAG laser emitting pulses at 1064 nm. The received backscatter signals are stored in 1024 range bins with a resolution of 15 m. The temporal resolution is set to 15 s. The signals are corrected for incomplete overlap by a correction function provided by the manufacturer.
A strict retrieval of the particle extinction coefficient from ceilometer measurements is not possible due to the unknown lidar ratio; furthermore, exploitation of the signal in the range of incomplete overlap is subject to errors. Thus, in order to convert the ceilometer measurements to aerosol extinction profiles, we followed an approach mentioned in Wagner et al. (2019). The range-corrected attenuated backscatter data from July 2016 to December 2017 were seasonally averaged. Data of the altitude between 500 m and 5 km above the instrument were averaged with a vertical grid resolution of 500 m. Data below 500 m were assumed to be constant, following the values at 500 m. The extinction coefficients were first calculated by scaling the attenuated backscatter profiles (β∗) to the seasonal average AODs at 360 and 477 nm obtained from the sun photometer. The extinction profiles were then used to correct for the attenuation of the backscatter profiles following the lidar equation (Klett, 1981; Fernald, 1984). The corrected backscatter profiles (β) were then scaled to the AODs at 360 and 477 nm measured by the sun photometer to obtain the extinction profiles; see Fig. 1. Note that the ceilometer measures at 1064 nm and the optical properties of aerosols depend on the wavelength. Therefore, the uncertainties of these profiles are very large and they should be considered qualitative only.
The results shown in Fig. 1 indicate that the aerosol load at the UFS is highest in summer (June, July, and August) and lowest in winter (December, January, and February). The seasonal results also indicate large variations in the aerosol load from the surface up to 2 km. Above 2 km the variability is smaller; however, their contribution to the total column is still substantial (∼30 %–50 %).
In this study, we developed an aerosol profile retrieval algorithm for MAX-DOAS measurements based on the lookup table method. According to the measurement sensitivity, we first parameterized the aerosol profile as the aerosol extinction coefficients of three altitude layers and defined a profile set which is assumed to include all possible profiles. O4 DSCDs corresponding to each profile in the set were simulated and stored in the lookup table. In the retrieval, O4 DSCDs are calculated from the measured spectra and then compared to the simulated ones corresponding to each profile of the set using a cost function. According to the cost function, valid profiles are selected from the set, and the optimal solution is defined as the weighted mean of all the valid profiles.
3.1 O4 DSCD calculation
The DSCDs of O4 were derived from both UV and VIS spectra using the DOAS technique (Platt et al., 1979; Platt and Stutz, 2008). In the retrieval, DSCD is defined as the difference between the slant column density (SCD) of each off-zenith spectrum () and the corresponding zenith reference spectrum (). The QDOAS spectrum analysis software (version 3.2) developed by BIRA-IASB (http://uv-vis.aeronomie.be/software/QDOAS/, last access: 1 April 2020) was used for the spectral fitting analysis. The calibration of the spectrometers was performed by fitting the measured solar spectra to the literature solar reference (Chance and Kurucz, 2010). All the measured spectra were first corrected for offset and dark current.
Details of the DOAS fit settings for both bands are listed in Table 1. The fitting windows were determined according to both the absorption signal of O4 and the SNR of the spectrometers. For UV spectra, the fitting window is 338–370 nm, which is the same for most of the other MAX-DOAS studies (e.g., Clémer et al., 2010; Wang et al., 2014a; Kreher et al., 2019), and it covers the strong absorption peak at 360.8 nm and a weak absorption peak at 344 nm. For VIS spectra, because the spectral range of the spectrometer begins at 427 nm and the SNR close to the spectral edges is low, we therefore adapted a smaller fitting window of 440–490 nm, which is a bit narrower than the fitting window of 425–490 nm commonly used in other MAX-DOAS studies (e.g., Clémer et al., 2010; Chan et al., 2017; Kreher et al., 2019). The VIS fitting window covers the strong absorption peak at 477 nm and a weak absorption peak at 446.5 nm. As the temperature at the UFS typically varies between 263 and 279 K (Risius et al., 2015), trace gas absorption cross sections measured at 273 K were used in the DOAS fit. Absorption cross sections of several trace gases as well as a synthetic ring spectrum were included in the DOAS fit. For each scanning cycle, the zenith spectra before and after the cycle were temporally interpolated to the measurement time of each off-zenith spectrum. The broadband spectral structures caused by Rayleigh and Mie scattering were removed by including a low-order polynomial in the DOAS fit. Small shift and squeeze of the wavelengths were allowed in the wavelength mapping process in order to compensate for small uncertainties caused by the instability of the spectrograph.
The root mean square (rms) of fit residual was used to evaluate the performance of the DOAS fit. DSCDs with a residual rms larger than were not considered in the following analysis. Under cloud-free conditions, the residual rms of most of the UV spectra varies between and , while the residual rms of most of the VIS spectra varies between and . This is because both the light intensity and the O4 absorption are stronger at the VIS band; hence the measurement SNR is higher.Volkamer et al. (2005)Chance and Orphal (2011)Rothman et al. (2010)Bogumil et al. (2003)Bogumil et al. (2003)Serdyuchenko et al. (2014)Serdyuchenko et al. (2014)Thalman and Volkamer (2013)Chance and Spurr (1997)
3.2 Cloud screening
The aerosol profile retrieval requires the forward simulation of the radiative transfer in the atmosphere. As the radiative transfer is rather complicated for cloudy-sky conditions, the forward simulation usually assumes a cloud-free atmosphere. The aerosol retrieval might result in large uncertainty under cloudy or foggy conditions. Therefore, it is important to filter out the measurements taken under cloudy or foggy conditions. In this study, a cloud screening approach based on a color index (CI) (Wagner et al., 2014, 2016) was applied to filter out cloudy measurements. The CI is defined as the ratio of radiative intensities at 330 and 390 nm in this study. Larger CI indicates the UV–VIS intensity ratio is higher; hence, the sky is more blue. Our cloud screening method is presented in Appendix A. The cloud screening results during the entire measurement period are summarized in Table 2. Among the four seasons, the percentage of cloudy measurements is highest in summer and lowest in winter. In total, about 60 % of the zenith measurements were determined as cloudy scenes, and the corresponding scanning cycles were not used in the aerosol profile retrieval.
3.3 Topography effect and the simplification in radiative transfer model
The topography around the UFS is quite complex, which complicates the radiative transfer simulations. As shown in Fig. 2, the surface altitude varies between 600 and 2800 m a.s.l. along the viewing direction of the MAX-DOAS instrument. Figure 2 also shows the type of surface in different colors which includes forests, meadows, rocks, etc. Some parts of the surface are seasonally or permanently covered by snow, while some steep slopes cannot be covered by snow even in winter.
Three-dimensional radiative transfer models (RTMs) can consider such a complex terrain, but they are computationally expensive and unaffordable for retrieval. Due to the limitation of the two-dimensional RTM LIDORT (Spurr et al., 2001; Spurr, 2008) used in the study, we simplified the ground topography to a flat surface at an altitude of 2650 m a.s.l. in the radiative transfer simulations. In order to estimate the error caused by this simplification, we investigated using the three-dimensional RTM TRACY-2.
TRACY-2 is a full spherical Monte Carlo atmospheric RTM (Deutschmann, 2008; Wagner et al., 2007), which allows the simulation of three-dimensional radiative transport as well as two-dimensional variation in the surface height. The model was compared to other RTMs and very good agreement was found (Wagner et al., 2007). We also performed an intercomparison with LIDORT. The result shows that with the same definition of topography and atmosphere, the difference between the O4 DSCDs simulated by the two RTMs is less than 3 %.
For the three-dimensional simulations carried out in this study, a pseudo-reality topography was defined with the exact ground altitude (obtained from Google Earth) in the azimuth direction of the MAX-DOAS measurements taken into account, whereas in the dimension orthogonal to this direction the surface altitude was set constant. This simplification was chosen to reduce the computational effort. We feel that this approach is justified since the atmospheric light paths in the viewing direction of the instruments can be very large (up to several tens of kilometers). It is most important to take this variation in the surface altitude along this direction into account, whereas the influence of the orography perpendicular to this direction is expected to be small.
Simulations were performed with all the combinations of three different SZAs (30, 50, and 70∘), three different relative solar azimuth angles (RAAs) (30, 60, and 90∘), and two different aerosol extinction profiles (an aerosol-free profile and a box-shape profile with AOD =0.12 and box height =3 km), i.e., altogether 18 cases. For each case, O4 DSCDs at 360 and 477 nm were simulated with both the flat surface at 2650 m and the pseudo-reality topography using TRACY-2. The relative errors of O4 DSCDs simulated with the flat surface compared to those simulated with the pseudo-reality topography are calculated. A fixed surface albedo of 0.07 was used in the simulations. For both wavelengths, the single-scattering albedo was set to 0.93 and the phase function was defined as a Henyey–Greenstein phase function with the asymmetry parameter set to 0.68. The atmospheric profile was defined as the US standard midlatitude atmosphere (Anderson et al., 1986). Figure 3 shows the results of some of the cases: panels (a) and (b) show the results of six cases with SZA =50∘ and different RAAs and both aerosol extinction profiles; panels (c) and (d) show the results of six cases with RAA =60∘ and different SZAs and also both aerosol extinction profiles.
As shown in all the panels of Fig. 3 as well as in all the other cases not shown, O4 DSCDs simulated with the flat surface are in general slightly underestimated compared to the pseudo-reality topography. The difference could be explained by the scattering in the valleys where the concentration of O4 is higher. For the flat surface at 2650 m, the light paths below 2650 m would not be taken into account, and hence the O4 DSCDs would be underestimated. Moreover, the relative error has no obvious correlation with elevation angle, SZA, RAA, and aerosol load. This is because the light path below 2650 m is influenced by the topography, and the influence differs with the observation geometry. In addition, the light path is also influenced by the aerosols both below and above 2650 m. Since only a pseudo-reality surface and a constant surface albedo are used in the study, the actual error caused by the topography simplification is expected to be much more complicated.
In order to make the compensation feasible, we consider the error to be the combination of a systematic error and a random error. Based on the results of all 18 cases of this study, the mean bias for each elevation angle and each wavelength is considered to be the systematic error, while the standard deviation of the relative difference is considered to be the random error; see Table 3. In the aerosol profile retrieval, systematic errors are first corrected from the measured O4 DSCDs, while random errors are included in the error budget in the calculation of cost functions (see Sect. 3.7.2). In the following text, measured O4 DSCDs refer to the values corrected by the systematic error unless otherwise mentioned.
3.4 Sensitivity analysis
In order to make full use of the measurement sensitivity and reduce unnecessary computational efforts, our retrieval algorithm was designed according to the sensitivity of O4 absorption. We performed several sensitivity analyses to determine the optimal vertical grid, step size of the aerosol extinction for each layer, and maximum aerosol extinction. In addition, these sensitivity analyses also help to estimate the measurement and model errors, which are very important for the retrieval. The sensitivity analyses are based on forward simulations of O4 DSCDs using LIDORT. We tested the sensitivities of O4 absorption to surface albedo, aerosol optical properties, and the aerosol vertical profile. The results of the sensitivity analyses are shown in Appendix B. The extreme and median values of the parameters are also discussed in that section.
3.5 Parameterization of the aerosol extinction profile
As discussed in Appendix B4, O4 absorption is insensitive to aerosols above 2 km. Therefore, our retrieval only focuses on aerosols between 0 and 2 km above the MAX-DOAS instrument (i.e., 2650–4650 m a.s.l.). In order to limit the complexity of the retrieval, avoid unreasonable results, and make full use of the measurement sensitivity, we parameterize the aerosol extinction profile as aerosol extinctions in three layers. The thicknesses of the two lower layers are defined as 0.5 km. Due to the lower sensitivity at high altitude, the thickness of the third layer is set to 1 km. The aerosol profile is denoted as a three-dimensional state vector x,
where σ1 is the aerosol extinction coefficient between 0 and 0.5 km (2650–3150 m a.s.l.), σ2 is the aerosol extinction coefficient between 0.5 and 1 km (3150–3650 m a.s.l.), and σ3 is the aerosol extinction coefficient between 1 and 2 km (3650–4650 m a.s.l.). The definition of x is illustrated in Fig. 4a. The vertical resolution of our retrieval grid is lower compared to that of many other studies (e.g., Clémer et al., 2010; Chan et al., 2017; Tirpitz et al., 2020); however, the vertical gradient of aerosol extinction at such a high-altitude site is expected to be small and this is also proved by the ceilometer measurements. Therefore, the coarse resolution is considered to be sufficient for the retrieval of the UFS MAX-DOAS measurements.
In order to formulate the lookup table, we defined a profile set (denote as XLUT) which is assumed to include all possible aerosol extinction profiles under cloud-free conditions. XLUT is a finite set of x, and the variation steps of σ1, σ2, and σ3 were determined according to the sensitivity and accuracy of measurement. XLUT includes only the profiles with reasonable shapes, and the variation range of σ1, σ2, and σ3 covers the actual aerosol load at the UFS. In this way, unreasonable and unrealistic retrieval results can be avoided.
As discussed in Appendix B6, the measurement sensitivity decreases with increasing surface aerosol extinction, and the sensitivity is very low when the surface aerosol extinction coefficient exceeds 0.3 km−1. Therefore, σ1 is defined to vary between 0 and 0.3 km−1. The variation step increases from 0.001 km−1 per step to 0.02 km−1 per step with increasing aerosol extinction, so that the difference of O4 DSCD per step is similar to the average spectral fitting error (∼2 %). In total, we define 65 values for σ1; see Table 4.
As illustrated in Fig. 4b, the values of σ2 and σ3 are defined as a tree, which means we define different values of σ2 for each σ1, and the values of σ3 are also defined depending on σ2. According to the ceilometer observations at the UFS, strong elevated aerosol layers are unlikely to exist under cloud-free conditions. Therefore we allow only weak elevated layers in designing the profile set. We assume that for reasonable profiles, σ2 should not exceed σ1 by more than 30 %, and σ3 should not exceed σ2 by more than 30 %, either. According to the sensitivity, for each value of σ1 (σ1>0), we define 14 possible values for σ2 which varies from 0 to 1.3σ1 with a step size of 0.1σ1. In the case σ1=0, elevated layers are not considered, and then σ2 and σ3 can only be 0. Similarly, σ3 varies between 0 and 1.3σ2. Due to the lower measurement sensitivity at high altitude, we define nine possible ratios between σ3 and σ2 (see Table 4). In the case σ2=0, σ3 can only be 0.
XLUT includes the profiles with all the combinations of σ1, , and . For each of the 64 nonzero values of σ1, there are corresponding profiles. For σ1=0, there is only one profile with . Therefore, the profile set consists of aerosol extinction profiles in total.
3.6 Definitions of other dimensions of the lookup table
The basic idea of the lookup table method is to replace the repetitive time-consuming computation with a pre-calculated array. In this study, we replace the forward simulation of O4 DSCDs with a lookup table, so that all the possible aerosol extinction profiles can be considered in the retrieval of each measurement cycle with an affordable computational effort.
Besides the parameterized aerosol extinction profile x, we consider another four input parameters for the forward simulation which can be described as a function,
where ΔSs refers to the simulated O4 DSCD, λ represents the wavelength, α indicates the elevation angle, θ is the SZA, and ϕ is the RAA. All the input parameters are well-known in the retrieval.
In order to formulate the lookup table, the input parameters need to be parameterized as a grid with finite nodes. As already presented in Sect. 3.5, the aerosol extinction profile (x) is parameterized as a profile set which consists of 7553 possible profiles. As the simulated O4 DSCDs are used to fit to the measured ones, only the data at 360 and 477 nm and at the six non-zenith elevation angles of the measurement cycles are included in the lookup table. SZA (θ) and RAA (ϕ) are parameterized as a grid with resolution. The grid includes 5005 combinations of SZA and RAA, which can cover all possible solar positions for the daytime measurements at the UFS. When we obtain data from the lookup table, as the input SZA and RAA are not integers, the output ΔSs is interpolated from the data of the four adjacent nodes of the SZA–RAA grid. In total, the five input parameters are parameterized as a grid with nodes. Details of the parameterization of the input parameters are summarized in Table 4.
As discussed in Appendix B, besides the input parameters we defined, O4 DSCDs can also be affected by other parameters such as the ground albedo, aerosol optical properties, and others. Since accurate measurements of these parameters are not available and their influence is relatively small, they are considered uncertainties. In creating the lookup table, these parameters were fixed to the median values. Details of the simulation settings are listed in Table 5. O4 DSCDs corresponding to all the nodes of the lookup table were simulated using LIDORT.
As discussed in Appendix B5, the influence from the aerosols above 2 km is also considered a kind of uncertainty and treated in a similar way as the other unknown parameters. In the simulations for creating the lookup table, the aerosol extinction coefficient between 2 and 4 km was defined as 0.5σ3, so that this so-called parameter is fixed to the “median” value. Note that the aerosol extinction coefficient above 2 km is neither considered a part of the retrieved profile nor counted in the retrieved AOD.
3.7 Error estimation
Most of the other MAX-DOAS studies only consider the spectral fitting error in their retrieval. However, this fitting error only contributes to a small part of the total error. In addition, the total error is not directly proportional to the spectral fitting error. As the measurement and simulation uncertainties play an important part in our inversion method, we perform a comprehensive error analysis for the MAX-DOAS measurement and radiative transfer simulation of O4 DSCDs. In this study, errors from seven major sources are taken into account in estimating the total uncertainty.
3.7.1 Error in measured O4 DSCDs
Two error sources related to measured O4 DSCDs are taken into account in the total uncertainty estimation, and they are the DOAS fitting error (ϵfit) and the error caused by temperature variation (ϵtemp).
ϵfit is the byproduct of the DSCD calculation, derived from the fit residual and the absorption cross section of O4. It is proportional to the rms of the fit residual. For low elevation angles (1, 2, 5∘), the percentage of ϵfit compared to the DSCD typically varies between 1 % and 3 % at the UV band and between 0.3 % and 0.7 % at the VIS band, which is rather small compared to other sources of error. However, for the elevation angle of 30∘, as the absolute DSCD value is much smaller, the percentage of ϵfit can be up to ∼25 % and ∼10 % at the UV and VIS bands, respectively.
As discussed in Sect. 3.1, the O4 absorption cross section measured at 273 K was used in the DOAS fitting. However, the effective temperature of the MAX-DOAS measurements could be significantly different from 273 K. Previous studies show that O4 absorption has a strong and systematic dependence on temperature (Thalman and Volkamer, 2013; Wagner et al., 2019). In order to estimate ϵtemp, we compared the O4 DSCDs retrieved using the cross sections measured at 253 and 293 K to those retrieved with the cross sections measured at 273 K. The comparison shows that the O4 DSCDs are underestimated by 5.1 % at the UV band and 2.5 % at the VIS band when the effective temperature is 293 K. On the other hand, the O4 DSCDs are overestimated by 6.9 % at the UV band and 3.9 % at the VIS band when the effective temperature is 253 K. These systematic errors are almost constant, regardless of the observation geometry. Between 253 and 293 K, the average variation rate of O4 DSCD at the UV band is 0.3 % K−1. This result is in general agreement with Wagner et al. (2019). They found that with the fitting window of 352–387 nm, O4 DSCDs retrieved using the cross section at 203 K are reported to be 30 % smaller than those retrieved using the cross section at 293 K, i.e., 0.33 % K−1 on average. Based on the fact that the temperature at the measurement site varies between ∼258 and 288 K during daytime in most cases, we estimate the ϵtemp of all measurements as 4.5 % and 2.4 % of the O4 DSCD at the UV and VIS bands, respectively.
3.7.2 Error in simulated O4 DSCDs
Five error sources related to simulated O4 DSCDs are taken into account in estimating the total uncertainty. They are the random error caused by the simplification of the topography definition (ϵtopo), the error caused by surface albedo (ϵSA), the error caused by single-scattering albedo (ϵSSA), the error caused by phase function (ϵPF), and the error caused by aerosols above retrieval height (ϵ2–4 km).
As discussed in Sect. 3.3, the random error caused by the simplification of the topography definition (ϵtopo) of each elevation angle and each wavelength is derived from the standard deviation of the relative errors of the 18 cases simulated using the three-dimensional RTM TRACY-2. Values of ϵtopo are listed in Table 3.
For the uncertainties from the other four sources (ϵSA, ϵSSA, ϵPF, and ϵ2–4 km), as discussed in Appendix B, they can be estimated by radiative transfer simulations. Since they differ under different observation geometries and different aerosol loads, we determine them using simple lookup tables in the retrieval. In order to simplify the error estimation process, we assume that the uncertainties from the four sources are only influenced by the AOD, while the influence from different vertical distribution of aerosols is neglected. In addition, from the O4 DSCD lookup table, we found that O4 DSCD at 5∘ is almost negatively correlated with AOD, while it is insensitive to the shape of the profile. Therefore, we use the O4 DSCD measured at 5∘ as the indicator for estimating the AOD in deriving uncertainty values from the error lookup tables.
The error lookup tables consist of the values of ϵSA, ϵSSA, ϵPF, and ϵ2–4 km for all the combinations of SZA and RAA (with resolution) and 65 profiles of the XLUT with . The calculation of the error lookup tables was similar to the sensitivity study. In order to estimate the uncertainty caused by each parameter, O4 DSCDs were simulated under both median and extreme values, while all the other parameters were fixed as the median settings listed in Table 5. The relative difference between the two simulations is treated as the uncertainty and stored in the lookup table.
As discussed in Appendix B1, the uncertainty caused by surface albedo (ϵSA) was derived from the relative difference of the O4 DSCDs simulated with the surface albedo set to 0.2 (extreme value) and 0.1 (median value).
As discussed in Appendix B2, in the estimation of the uncertainty caused by single-scattering albedo (ϵSSA), the extreme value was chosen as 0.997 for both the UV and VIS bands, while the median value was chosen as 0.92 and 0.93 for the UV and VIS bands, respectively.
As discussed in Appendix B3, from all the phase functions measured by the AERONET station in Hohenpeißenberg during the period of 2013–2014, the phase function with which the simulated O4 DSCDs at all elevation angles are closest to the median values was chosen as the so-called median phase function. The phase function with which the simulated O4 DSCDs are closest to the rank of 95 % (i.e., 2σ) was chosen as the “extreme” phase function. ϵPF was derived from the relative difference between O4 DSCDs simulated with median and extreme phase functions.
As discussed in Appendix B5, the error caused by aerosols above 2 km (ϵ2–4 km) is treated similarly to ϵSA, ϵSSA, and ϵPF in the study. The so-called median O4 DSCDs were simulated with profiles with the aerosol extinction coefficient between 2 and 4 km equal to 0.5σ3 (50 % of the aerosol extinction coefficient between 1 and 2 km), while the extreme values were simulated with the aerosol extinction coefficient between 2 and 4 km set equal to σ3. ϵ2–4 km was derived from the relative difference between the extreme and median results.
3.7.3 Total uncertainty
We assume that the seven kinds of errors mentioned in Sect. 3.7.1 and 3.7.2 follow the normal distribution, and the total uncertainty of each band and each elevation angle can be determined by the root mean square of the seven errors as
Examples of the error budgets of two measurement cycles for both wavelength bands are shown in Figs. 5 and 6. The cycle shown in Fig. 5 was measured in summer under relatively high aerosol load (AOD at 440 nm measured by the sun photometer around noon of that day was ∼0.2), while the cycle shown in Fig. 6 was measured in winter under relatively low aerosol load (AOD at 440 nm measured by the sun photometer around noon of that day was ∼0.015). In addition, the former cycle was measured under a smaller SZA compared to the latter one (64 and 79∘, respectively), while the RAA was much larger than that of the latter (97 and 39∘, respectively). The results show that contributions from different error sources are quite different in different measurement cycles, at different wavelengths, and at different elevation angles.
3.7.4 Other possible error sources
Besides the seven abovementioned error sources, there are still some other sources of error which are difficult to estimate and hence not included in the error estimation.
Error in O4 DSCD scaling factors. In this study, we found that an elevation-dependent O4 DSCD scaling factor is needed to bring measurements and modeled results into agreement. We determined the factors based on the statistical analysis of the long-term measurements; see Sect. 3.9. However, as it is still difficult to estimate the uncertainties of the scaling factors, they are currently not taken into account in calculating the total uncertainty.
Error caused by horizontal gradients of the aerosol extinction. Besides its direct effect on the measurements, the complex topography might also cause systematic horizontal gradients of the aerosol extinction. For example polluted air masses from the valleys might be transported to higher altitudes according to the vertical mixing and the prevailing wind direction. Such effects can be especially important for the measurements discussed here because of the rather low AOD. Further quantification of the effects of possible horizontal gradients is beyond the scope of this study but might be one reason for the observed elevation dependence of the O4 DSCD scaling factor.
Error caused by the variation in atmospheric profile. The O4 DSCD lookup table was calculated using the US standard climatology data, but the change of atmospheric temperature and pressure can slightly affect the O4 absorption. However, since it is difficult to estimate the accurate uncertainty, and real-time measurements of temperature and pressure profiles are not available, the error caused by the variation in the atmospheric profile is not taken into account in calculating the total uncertainty.
Systematic effect of the surface albedo on the measurements at the high-altitude station. Due to the dependence of the snow coverage on altitude, the surface albedo close to the instrument is typically higher than at locations far away. Since the measurements at high elevation angles are usually more sensitive to air masses closer to the instrument, they are probably more strongly affected by snow and ice than measurements at low elevation angles. In this study, this effect cannot be further quantified, but it might be one reason for the need for different O4 DSCD scaling factors for different elevation angles; see Sect. 3.9.
In order to avoid the underestimation of the measurement uncertainty, we set a relatively relaxed threshold of cost functions for choosing valid profiles; see Sect. 3.8.
3.8 Inversion method
Aerosol extinction profiles are retrieved from the measured O4 DSCDs of each scanning cycle. The measurements of the UV and VIS bands are retrieved separately. The measured O4 DSCDs at the UV and VIS bands are fitted to the simulated O4 DSCDs at 360 and 477 nm, respectively. In the retrieval, we assume the state of atmosphere to be stable during a scanning cycle and the distribution of aerosols to be homogeneous in horizontal direction. For a single scanning cycle, the measured O4 DSCDs at the wavelength λ are denoted as a measurement vector
where M is the number of off-zenith measurements in each scanning cycle and is six in this study. ΔSλ,1, ΔSλ,2, …, ΔSλ,6 are the O4 DSCDs measured at O4 wavelength band λ with the viewing elevation angles of 1, 2, 5, 10, 20, and 30∘, respectively.
The simulated O4 DSCDs corresponding to each possible aerosol extinction profile in XLUT can be obtained from the lookup table. Similar to ym, the simulation vector ys for each possible profile x is denoted as
Aerosol extinction profiles can be derived by fitting the forward simulation to the measured O4 DSCDs. Typically, the optimal solution can be determined by minimizing the cost function, which is defined as
where Sϵ is the uncertainty covariance matrix. Assuming the measurements of each viewing elevation angle are independent, Sϵ is a diagonal matrix and its diagonal elements are equal to the square of the total uncertainties of each elevation angle defined in Eq. (3),
Our cost function definition is similar to the cost functions used in many of the MAX-DOAS studies based on the OEM (e.g., Clémer et al., 2010; Frieß et al., 2016; Wang et al., 2016; Chan et al., 2017), but it only includes the item related to measurement error, while the item related to the a priori profile is omitted. This is because the a priori profile is not needed in our retrieval algorithm.
χ2 indicates the difference between ys and ym. However, as the retrieval is ill-posed and the SNR of the measurements at the UFS is low, the single profile with the lowest χ2 is not necessarily the one closest to the true profile. In order to overcome this limitation, we consider all the profiles in XLUT with χ2(x) ≤ 1.5M (nine in this study) to be valid profiles and calculate the weighted mean profile as the optimal result. A profile with χ2 ≤ M indicates that the measured and simulated O4 DSCDs agree within the measurement errors, but in order to avoid underestimation of the measurement errors, we define the threshold as 1.5M. The weight of each valid profile for the calculation of the optimal solution is defined as
and the optimal solution can be calculated as
3.9 O4 DSCD correction
Discrepancies between measured and simulated O4 DSCDs are found in many other MAX-DOAS studies (Wagner et al., 2009, 2019; Clémer et al., 2010; Chan et al., 2015, 2017; Wang et al., 2016). The discrepancies are often explained by the systematic errors of the absorption cross section of O4 as well as the radiative transfer simulation, and a correction is therefore necessary. Previous studies suggested multiplying a constant scaling factor (usually between 0.75 and 0.9) with the measured O4 DSCD for all elevations to correct for the systematic error (e.g., Wagner et al., 2009; Clémer et al., 2010; Chan et al., 2015; Wang et al., 2016). Some recent studies suggested elevation-dependent scaling factors. Irie et al. (2015) suggested a set of scaling factors for 477 nm which gradually decreases with increasing elevation angle, varying from 0.984 for 1∘ to 0.667 for 30∘. Zhang et al. (2018) suggested a set of scaling factors for 360 nm, which also decreases with increasing elevation angle, varying from 1.02 for 1∘ to 0.909 for 30∘. Chan et al. (2017) derived a set of elevation-dependent scaling factors for 477 nm by comparing modeled and measured (relative) intensity, varying from 0.792 for 1∘ to 0.957 for 30∘. On the other hand, some other MAX-DOAS studies did not find it necessary to apply any correction to O4 DSCDs. For example, Frieß et al. (2011) reported that for the MAX-DOAS measurements in an Arctic area, the measured and simulated O4 DSCDs are in good agreement without any correction. Note that the scaling factor mentioned here refers to the ratio between simulated and measured O4 DSCDs, which is opposite to some other studies.
In order to assess whether the O4 DSCD correction is necessary for the MAX-DOAS measurements at the UFS, we compared the measured O4 DSCDs to the simulated ones in the lookup table. Assuming our profile set (XLUT) covers all possible aerosol profiles under cloud-free conditions, we derived the O4 scaling factor for each elevation angle and each wavelength based on the statistical analysis. The AODs measured by the sun photometer were used to restrict the range of possible profiles.
Figure 7 shows the scattered plots of measured and simulated O4 DSCDs of the scanning cycle on 7 December 2015 at ∼13:55 UTC. The measurements of both (a) UV and (b) VIS bands are shown. According to the cloud screening as well as the skycam images, this day was absolutely cloud free. Total AOD measured by the sun photometer at that time is 0.02 and 0.017 for the 360 and 477 nm bands, respectively. In each plot, the x axis indicates the O4 DSCDs measured (or simulated) at the elevation angle of 1∘, while the y axis represents the O4 DSCDs at the other elevation angles. Different colors indicate measurements at different elevation angles. The simulated O4 DSCDs (ys(x)) of all the possible profiles in XLUT are shown as colored dots. We assume the MAX-DOAS measurement of AOD between 0 and 2 km (denoted as τ2k, ) varies between 50 % and 100 % of the total AOD measured by the sun photometer (denoted as τsp,λ) in most cases, and the data points of the profiles fulfilling this assumption are highlighted. The measured O4 DSCDs (already corrected for the systematic errors caused by the topography simplification) are plotted as square markers with error bars showing the total uncertainties. It is obvious that at most of the elevation angles, the measured O4 DSCD does not agree with the simulations within the total error. As a result, at both the UV and VIS bands, no profiles in XLUT satisfy the selection requirement (χ2 ≤ 9; see dashed curves in Fig. 8). No profiles matching the measurement is unlikely to happen under such clear-sky conditions, hence implying a systematic error, and correction of the error is necessary.
In order to determine whether the O4 scaling factor is constant for all elevations or dependent on the viewing elevation angles, we first assume it is constant and plot the corrected O4 DSCD measurements in Fig. 7. The plus signs indicate the measured O4 DSCDs corrected with constant scaling factors of 0.8, 0.9, 1.1, and 1.2. Furthermore, the corrected O4 DSCDs should vary along the colored dashed lines if any other constant scaling factor is applied to the measurements. However, the forward simulation of O4 DSCDs does not overlap with the dashed lines in most of the cases (especially for 5 and 10∘ of the UV band), indicating that a constant O4 scaling factor for all viewing elevation angles could not resolve the systematic error. Therefore, different scaling factors should be applied to different elevation angles.
In this study, the O4 DSCD scaling factors for each viewing elevation angle and wavelength were determined through the statistical analysis of the long-term observations. We assume the scaling factor mainly depends on the viewing elevation angle, while being less sensitive to other factors such as solar geometry, aerosol load, temperature, etc.
Figure 7 shows that the simulated O4 DSCDs at high elevation angles (e.g., 20 and 30∘) vary in a very narrow range. Based on the assumption that XLUT covers all possible aerosol profiles, the measured O4 DSCDs should lie within the range. The scaling factor can be derived by taking the ratio of the simulated and measured values. As the simulated value varies in a narrow range, the uncertainty of the derived scaling factor should also be low. In order to have better statistics of the scaling factors, this method was applied to the long-term measurements. In addition, only the measurements taken under cloud-free and low-aerosol-load () conditions were used, so as to avoid accounting for data contaminated by clouds in the analysis. Here it should be noted that measurements with AOD <0.03 are almost entirely found during winter due to the strong seasonal variation in aerosol load at the UFS. Subsequently, for the wavelength λ and the ith elevation angle of each scanning cycle, we calculate the variation range of the simulated O4 DSCDs for all the profiles fulfilling , which can be described as a set,
Only if , was the scanning cycle then taken into account. In most cases, measured O4 DSCDs at high elevation angles are lower than simulated ones. Therefore we calculate the scaling factor from the minimum value in to avoid overestimation of the scaling factor. The scaling factor derived from this scanning cycle is denoted as
where ΔSλ,i is the measured O4 DSCD (already corrected for the systematic errors caused by the topography). For the elevation angles of 5, 10, 20, and 30∘ at the UV band and 10, 20, and 30∘ at the VIS band, numerous scanning cycles from the long-term measurements fulfill the selection criterion, and hence there are sufficient samples of for statistical analysis. We analyzed the frequency distribution of of each elevation and each wavelength band. The result shows that the distributions of follow the normal distribution function with small standard deviation. For instance, for the elevation angle of 20∘, the standard deviations of UV and VIS bands are both ∼0.16. Subsequently, with the maximum frequency was derived by Gaussian fit. The peak value was used as the scaling factor, which is denoted as .
For the low elevation angles (1 and 2∘ at the UV band and 1, 2, and 5∘ at the VIS band), as O4 DSCD varies in a wide range, it is impossible to determine the scaling factor with the method mentioned above. However, it is found that in many scanning cycles, within the possible profiles in XLUT, the simulated O4 DSCDs at low elevation angles are well-correlated to those at the neighboring elevation angle. Therefore, once the scaling factor of the higher elevation angle is determined, we can derive an expected value of the O4 DSCD at the lower elevation angle from the corrected O4 DSCD at the higher one, and the scaling factor can be derived by taking the ratio of the expected value and the measured value.
For the wavelength λ and for each scanning cycle, a subset of XLUT is defined as
and the elements of X† are denoted as . The corresponding simulated O4 DSCD at the ith elevation angle is denoted as
A third-order polynomial regression is applied between and . The regression function is denoted as g. Only if the correlation coefficient R2 ≥ 0.98 is this scanning cycle taken into account. As the scaling factor of the (i+1)th elevation () is already determined, the expected value of the O4 DSCD at the ith elevation angle can be derived with the regression function
and the scaling factor derived from this scanning cycle is
Similar to the high elevation angles, the frequency distribution of from all the available samples was analyzed by fitting to a Gaussian function. The peak value of is used as . The scaling factor of the (i−1)th elevation is then derived in the same way. The scaling factors of 1 and 2∘ at the UV band and 1, 2, and 5∘ at the VIS band were determined using this method.
* The O4 DSCDs which are already corrected for the systematic errors caused by the topography simplification.
The determined scaling factors are listed in Table 6. The corrected O4 DSCDs are indicated as triangles in Fig. 7. The result shows that except for the elevation angle of 1∘, the simulated O4 DSCDs are overestimated compared to the measured ones. It should be noted that the determination of the scaling factors is based on the measured O4 DSCDs which are already corrected for the systematic errors caused by the topography simplification (discussed in Sect. 3.3). Comparing to the original measurements, the result still indicates that the simulated O4 DSCDs at high elevation angles are overestimated. This result is opposite to the results of the other studies. At the moment we have no clear explanation for this finding, it might be related to the specific properties of the high-altitude station, e.g., the highly structured topography, horizontal gradients of the aerosol extinction, and systematic dependence of the surface albedo on altitude.
Figure 8 shows the cumulative distribution of χ2 of all the profiles in XLUT for the scanning cycle shown in Fig. 7. The distribution of χ2 before and after the DSCD correction is shown as dashed and solid curves, respectively. The result indicates that for both UV (blue curves) and VIS (red curves) bands, the χ2 values of most profiles in XLUT are significantly lower after the correction. As a result, a number of profiles fulfill the selection criterion (χ2 ≤ 9). Note that the AODs measured by the MAX-DOAS are still expected to be lower than the sun photometer observations due to the fact that the MAX-DOAS only reports the AOD below 2 km while the sun photometer covers the entire atmosphere.
Our retrieval algorithm was applied to the long-term measurement data of the UFS MAX-DOAS from February 2012 to February 2013 and from July 2013 to February 2016. The results are also compared to sun photometer measurements. This section presents the results as well as their discussion.
4.1 Dependency of retrieval result on the threshold of cost function
As presented in Sect. 3.8, we consider all the profiles with χ2≤9 to be valid profiles, and the retrieved profile is defined as the weighted mean of all the possible profiles. In this section, we investigate the dependency of the retrieval result on the threshold of χ2 by comparing the results calculated with different χ2 thresholds. Taking the two measurement cycles mentioned in Figs. 5 and 6 for example, Fig. 9 (5 July 2015 at ∼16:26 UTC) and Fig. 10 (7 December 2015 at ∼13:55 UTC) show the weighted mean profiles, the variation range of valid profiles and the number of valid profiles corresponding to different χ2 thresholds. The profiles are shown as colored curves which indicate the aerosol extinction coefficients in the three layers (i.e., σ1, σ2, and σ3).
The results of both scanning cycles show that the retrieved profile is not sensitive to the threshold of χ2 when there is a sufficient number of valid profiles (number of profiles exceeds ∼800 and ∼400 for UV and VIS, respectively; see the grey curves in Figs. 9 and 10). This is because profiles with larger χ2 have lower weight (w). In addition, when the threshold value is increased, more profiles with both higher and lower aerosol extinction coefficients are taken into account. As a result, the variation range of valid profiles becomes larger but the weighted mean remains similar. The result shows that the retrieval with a χ2 threshold of 9 is stable; therefore, it is used in the study.
4.2 Estimation of the uncertainties of retrieved profiles
Still taking the two measurement cycles mentioned in Sect. 4.1 as examples, we analyzed the weight distribution of valid profiles; see Figs. 11 and 12. The distributions of aerosol extinction coefficients in the three layers (σ1, σ2, and σ3) are shown as solid curves. For each layer, aerosol extinction coefficients of all the valid profiles are grouped, and the y axis refers to the total weight of each group. The three vertical dashed lines indicate the weighted mean aerosol extinction coefficient of each layer (i.e., σ1, σ2, and σ3 of ). The result shows that the distributions of σ1, σ2, and σ3 are all asymmetric for both the UV and VIS bands. In particular for the layer of 1–2 km (σ3) at the UV band, the weight decreases monotonically with increasing aerosol extinction in both of the two cycles. Taking the cycle shown in Fig. 12 (7 December 2015 at ∼13:55 UTC) as an example, there are altogether 205 (12.8 %) and 120 (12.6 %) valid profiles with σ3=0 in the UV and VIS bands, respectively. These profiles contribute total weights of 0.122 and 0.101 for the UV and VIS retrievals, respectively.
In order to estimate the uncertainty of , we calculate the weighted standard deviations of σ1, σ2, and σ3 of all the valid profiles. Due to the asymmetric distribution, the weighted standard deviations are calculated separately for the left (negative) and right (positive) sides. For the lth () layer, denote the aerosol extinction coefficient of each profile as σl(x), then the weighted standard deviation of the left side is calculated from all the valid profiles with ,
and the weighted standard deviation of the right side is calculated from all the valid profiles with ,
The uncertainties of are indicated as error bars in Figs. 11 and 12. For each layer, the total weight of the profiles covered by the error bar is labeled in the charts. At the UV band, the total weight of the valid profiles covered by the uncertainties is 59 %–66 %, which is close to the standard normal distribution. However, the percentage can be up to 90 % at the VIS band. This is because the SNR of the measurement at the VIS band is higher. Therefore the retrieval of the VIS band has higher selectivity, and the weight is more concentrated to the mean value.
4.3 Retrieval of synthetic measurement data
In order to test the effectiveness of our retrieval algorithm, we generated some synthetic measurement data for the application to our algorithm. Figure 13 shows the results of three representative synthetic profiles at 360 and 477 nm. In each chart, the true profile is shown as the grey curve. Profile 1 is a tangent curve with aerosols distributed between 0 and 6 km above the instrument. The aerosol extinction decreases with increasing altitude, which is 0.04 km−1 at surface level, ∼89 % at 2 km, and 50 % at 3 km. The total AOD is 0.12, of which ∼92 % is contributed from the altitude below 3 km. Profile 2 has a similar shape as Profile 1, but the aerosol extinction between 0.5 and 1 km above the instrument was enhanced. The aerosol extinction peaks at 0.75 km, and the average aerosol extinction coefficient between 0.5 and 1 km is larger than the bottom layer by ∼10 %. In addition, the aerosol extinction coefficients at other altitudes are increased by a factor of 2 compared to Profile 1. Profile 3 is an exponential profile. The total AOD is 0.12, the scaling height is 1.5 km, and the surface aerosol extinction coefficient is 0.03 km−1.
We first simulated O4 DSCDs at 360 and 477 nm with each profile. The solar position was set as SZA =60∘ and RAA =60∘, and the other parameters followed the settings used in calculating the lookup table listed in Table 5 (excluding the aerosol extinction coefficients above 2 km). In order to test the stability of the retrieval, we also generated a set of noisy data for each profile and each wavelength by adding random noise to the simulated O4 DSCDs. We assume the measurement noise at all elevation angles is the same and follows a normal distribution with a standard deviation of 2 % of the DSCD of the lowest elevation angle. This noise level is realistic for the measurements at the UFS.
Aerosol profiles were then retrieved from both the original and noisy synthetic data using our algorithm. In the error estimation, the DOAS fitting error (ϵfit) was defined as the average values of the UFS measurements, while the other six kinds of errors followed the common settings presented in Sect. 3.7. O4 DSCD correction was not applied. The solid and dashed blue curves in Fig. 13 show the profiles retrieved from the original and noisy data, respectively, and the error bars indicate the uncertainties calculated by Eqs. (16) and (17). The results show that for Profile 1 and Profile 3 our retrieval algorithm can reproduce the true profiles well from not only the original data but also the noisy data. For Profile 2, the retrieved profile cannot reproduce the elevated layer, but the error bar covers the aerosol extinction of the true profile. This is because the retrieval is ill-posed, which means the limited input information does not correspond to a unique profile with an elevated layer. Instead, many other profiles without the elevated layer can also fit the input information. Adding noise to the synthetic data can affect the retrieved aerosol extinction coefficients. However the influence is small in most cases. In addition, the noise can amplify the uncertainty of the retrieved profile. The results indicate that our LUT-based retrieval is stable.
We also retrieved the synthetic data using the bePRO profiling tool developed by BIRA-IASB (Clémer et al., 2010; Hendrick et al., 2014). It is an OEM-based algorithm and uses LIDORT as the forward model. In the retrieval of all six cases, the a priori profile was defined as an exponential profile with AOD =0.12 and scaling height =1.5 km, shown as the dotted orange curve in each panel of Fig. 13. The vertical grid was defined as 20 layers of 200 m thickness each. For Profile 1 and Profile 2, the uncertainty covariance matrix of the a priori profile (Sa) was defined as in Clémer et al. (2010) and Wang et al. (2014a): the diagonal elements corresponding to the bottom layer, Sa (1, 1), were set as the square of a scaling factor β (β=0.2) times the maximum partial AOD of the profiles; the other diagonal elements decrease linearly with increasing altitude to 0.2×Sa (1, 1); the off-diagonal elements of Sa were defined using Gaussian functions with correlation length γ=0.05 km. For Profile 3, as the difference between the true and a priori profiles is quite large, we set β=0.4 and γ=0.1 km, so that the constraint from the a priori profile is weaker. The measurement uncertainty covariance matrix (Sϵ) was also defined as in most of the other MAX-DOAS studies so that Sϵ is a diagonal matrix with variances equal to the square of the DOAS fitting error (). We defined ϵfit the same as in the LUT retrieval, but the other six error sources were not included. The retrieval parameters related to the radiative transfer simulation followed the settings of our LUT-based retrieval.
The results retrieved from the data with and without noise are shown in Fig. 13 as solid and dashed red curves, respectively. In all 12 retrieval cases, the O4 DSCDs simulated with retrieved profiles are well-correlated to the input values (the relative root-mean-square error varies between 0.7 % and 4.7 %). However, as the retrieval is ill-posed, the retrieved profiles cannot reproduce the true profile well. Especially at high altitudes (above 1 km), the retrieved profiles are mostly dominated by the a priori profile. The OEM retrieval is also sensitive to measurement noise, which can be seen from the large variations in profile shape and aerosol extinction. The results indicate that the LUT-based algorithm is much more suitable for measurements with low SNR.
4.4 Comparison to sun photometer measurements
Figure 14 shows the comparison of AODs measured by the MAX-DOAS and sun photometer during the entire study period. The seasonally averaged AODs measured by both instruments are listed in Table 7. As the AOD measured by the MAX-DOAS refers to the AOD between 0 and 2 km while the AOD measured by sun photometer refers to the total AOD, the sun photometer results should be larger. Despite the difference, the time series (panels (a) and (c) of Fig. 14) show that the AODs measured by both instruments have a similar seasonal variation with higher AOD in summer and lower AOD in winter. The monthly average data show that the difference between the AODs measured by the MAX-DOAS and sun photometer is much larger in summer, which coincides with the ceilometer profiles shown in Fig. 1 which indicate much higher aerosol extinction coefficients above 2 km in summer. The underestimation of the MAX-DOAS may also be related to the decreased sensitivity of measurement at higher altitudes.
The correlation between hourly averaged AODs measured by the MAX-DOAS and sun photometer is shown in Fig. 14b, d. AODs show a general agreement at the UV and the VIS bands with correlation coefficients of R=0.733 and 0.798, respectively. However, AODs from the MAX-DOAS are lower; consequently the slope of the regression lines is 0.5308 and 0.3556 for the UV and VIS bands, respectively. As the MAX-DOAS only reports AODs below 2 km while the sun photometer measures the total AODs, the MAX-DOAS AODs are indeed expected to be lower. This is in particular true in cases of large AODs due to very strong convection of polluted air masses from the valley and/or the presence of Saharan dust layers. Then, particles are often transported beyond the range of the MAX-DOAS measurements and the disagreement is largest. This feature might be strengthened by the decreased sensitivity of the MAX-DOAS measurement at higher altitudes, so that the upper part of an aerosol layer is missed. In addition, a few data points lie above the 1:1 reference lines. This might be explained by the inhomogeneous distribution of aerosols in the horizontal direction. The light paths of the MAX-DOAS and the sun photometer are different. The MAX-DOAS measures scattered sunlight while the sun photometer derives the AOD from direct sun measurements. Therefore, when the aerosol load along the light path of the MAX-DOAS is higher than that of the direct sun measurement, the AOD measured by the MAX-DOAS may exceed the AOD measured by the sun photometer. For most of these points, the difference between the results of the two instruments is within their uncertainty ranges; i.e., the disagreement is probably due to the measurement and retrieval errors.
4.5 Temporal variation in aerosol characteristics
The seasonally averaged aerosol extinction profiles derived from the long-term measurements are shown in Fig. 15. The result indicates that the aerosol load is high in summer and low in winter, which coincides with the ceilometer results shown in Fig. 1. The seasonal pattern can be explained by the higher biogenic emissions from vegetation in summer. Moreover, the mixing layer is higher in summer; thus anthropogenic aerosols are more likely dispersed to upper altitudes. The shape of the profiles also agrees with the ceilometer results that the averaged aerosol extinction decreases with increasing altitude in all seasons – taking into account the coarse vertical resolution of the MAX-DOAS. In addition, Fig. 15 shows a much larger vertical gradient at 360 nm in summer. This might be explained by the lower sensitivity of the UV measurement for higher altitudes due to the more decreased visibility at shorter wavelengths.
We compared the seasonally averaged aerosol extinction coefficients at 360 and 477 nm in the bottom layer (0–0.5 km above the instrument, σ1); see Fig. 16. The averaged aerosol extinction coefficients are shown as bar charts. The ratio between the aerosol extinction coefficients at 360 and 477 nm is indicated by the grey curve. The result shows that the aerosol extinction coefficient ratio between 360 and 477 nm is significantly higher in summer than in the other seasons.
From these ratios the Ångström exponent (AE) can be calculated using the seasonally averaged surface aerosol extinction coefficients at 360 and 477 nm. The results are listed in Table 8. The seasonal averaged AEs of 380–500 nm from the AERONET measurements at Hohenpeißenberg from April 2013 to February 2016 are also listed for comparison. The result shows that both the UFS and Hohenpeißenberg measured the highest AE in summer and the lowest in winter. The AE at the UFS is in general lower than that measured at Hohenpeißenberg with a smaller difference in summer. This can be explained by the different altitude of the two sites. As the AERONET station at Hohenpeißenberg is located at ∼950 m a.s.l., a larger contribution of anthropogenic aerosols is expected. The extremely low AE at the UFS in spring, autumn, and winter agrees with the result measured at a plateau site (Lhasa, China, 3688 m a.s.l.) reported by Xin et al. (2007). The annual mean AE at that site is reported to be 0.06±0.31, which is significantly lower than that measured at low-altitude sites, especially urban and forest sites. In general, a smaller AE implies larger aerosol particle sizes (Dubovik et al., 2002). The increased AE at UFS in summer indicates a larger contribution of fine particles. The result is consistent with the fact that the particle size of biogenic secondary aerosols is in general smaller than ice particles transported from the lower altitudes to upper altitudes in summer.
We have developed a new MAX-DOAS aerosol profile retrieval algorithm based on a parameterized O4 DSCD lookup table. The algorithm is applied to the long-term MAX-DOAS measurements at the UFS, Germany, a high-altitude site located at 2650 m a.s.l.
Observations of O4 absorptions at both 360 and 477 nm were analyzed. We first investigated the sensitivities of O4 absorption to several parameters. According to the sensitivity analysis result, we defined an aerosol profile set which consists of 7553 possible profiles and then simulated O4 DSCDs with all the profiles and all possible observation geometries. In the retrieval of each measurement cycle, the simulated O4 DSCDs corresponding to all the possible profiles are obtained from the lookup table. The cost functions (χ2) are calculated for each possible profile according to the simulated and measured O4 DSCDs as well as the measurement uncertainties. A comprehensive error analysis is performed to estimate the total uncertainty. Valid profiles are selected from the profile set according to the cost function. The optimal solution is defined as the weighted mean of the valid profiles.
One key result of our study is that an elevation-dependent O4 DSCD scaling factor is needed to bring measured and simulated O4 DSCDs into agreement. Assuming the lookup table covers all possible aerosol profiles under clear-sky conditions, we determined the scaling factors based on the statistical analysis of the long-term measurements. The agreement between measured and simulated O4 DSCDs is greatly improved by this correction.
In addition, we developed a simple cloud screening method which is based on the statistical analysis of the color index. The developed cloud screening method is applied to the long-term measurements to filter out data taken under cloudy conditions.
In order to test the effectiveness of the algorithm, we retrieved profiles from synthetic data. The results indicate that our algorithm can reproduce the true profile well, and the retrieval is stable to measurement noise.
The AOD retrieved from the long-term MAX-DOAS measurements was compared to the sun photometer observations at the UFS. The results show reasonable agreement with each other. However, especially in summer, the sun photometer results are systematically larger (by about a factor of 2) than the MAX-DOAS results. This discrepancy is due to the different definitions of AOD measured by the MAX-DOAS and sun photometer. The larger difference in summer also coincides with the ceilometer measurements at the UFS which indicate larger aerosol extinctions at high altitude in summer. The long-term observation results show that the aerosol load at the UFS is higher in summer and lower in winter. Higher AOD in summer is mainly related to a higher frequency of extended mixing layers that allows particles to disperse from lower to upper altitudes. According to the MAX-DOAS measurements, the mean aerosol extinction decreases with increasing altitude for all seasons, which agrees with the ceilometer measurements. The Ångström exponent derived from MAX-DOAS surface measurement is higher in summer and extremely low in winter, which implies a smaller particle size in summer. This might be due to a significant contribution from biogenic sources in summer.
The study demonstrated that the developed method is effective for MAX-DOAS measurements at the UFS. Since the profile set only consists of reasonable profiles and the final solution is calculated from the weighted mean of all valid profiles, and because the retrieval does not rely on a priori profiles, many of the limitations of retrieval algorithms based on the optimal estimation method are overcome. In addition, as the O4 DSCDs of all possible profiles are pre-calculated, our method significantly reduces the computational time, so that real-time retrievals should be possible.
In this study, the color index (CI) is defined as the ratio of radiative intensities at 330 and 390 nm. Measured CIs (denoted as CImeas) were calculated from the zenith UV spectra (offset and dark current corrected) by taking the ratio of the counts at 330 and 390 nm. Figure A1 shows the time series of CImeas calculated from all the zenith spectra with SZA (solar zenith angle) during the entire study. The result shows that the variation range of CImeas is stable within the two periods. However, the optical throughput of the instrument in the UV spectral range has been enhanced after regular maintenance of the optical system in 2013. Hence, the CI increased systematically in the second period. Therefore, calibration of CImeas is necessary in order to make the CImeas values measured during the two periods comparable to each other. The calibration was performed following the method suggested in Wagner et al. (2016). CImeas values measured under overcast skies were fitted to the simulated minimum CI. The correction factor was determined to be 2.70 and 2.06 for the periods of February 2012–January 2013 and August 2013–February 2016, respectively. CImeas was subsequently converted to CIcal (calibrated CI) by multiplying the corresponding correction factor.
Figure A2 shows the frequency distribution of CIcal measured with different SZAs. The CIcal values from the long-term measurements were grouped by their SZA with a step size of 2∘. The relative frequency distributions are color coded. The result shows a bimodal frequency distribution of CIcal for all SZAs. The peaks with lower and higher CIs correspond to measurements under overcast and clear skies, respectively. This pattern is similar to the CI measured on Jungfraujoch, Switzerland (3570 m a.s.l.), reported by Gielen et al. (2014) and different from the results measured at the low-altitude sites reported by Gielen et al. (2014) and Wagner et al. (2016). This is because the high-altitude sites are seldom influenced by anthropogenic aerosols; hence the sky is either clear or covered by cloud or fog most of the time. Based on this feature, we defined the threshold for cloud screening as the CIcal with the minimum probability between the two peaks (denoted as CIcal,valley). The CIcal,valley was determined by fitting the probability density function to a Gaussian function. The circle markers shown in Fig. A2 indicate the determined CIcal,valley. In order to minimize the noise, the CIcal,valley was fitted to a fourth-order polynomial. The resulting smoothed CIcal,valley was used as the threshold (indicated as dashed curve in Fig. A2). Based on this approach, ∼60 % of the zenith measurements were determined to be cloudy scenes, and the corresponding scanning cycles were not used in the following analysis.
We investigated the sensitivity of O4 absorption to surface albedo, single-scattering albedo (SSA), scattering phase function, aerosol extinction at different altitudes, aerosol extinction above retrieval height, and surface aerosol extinction. In the test for each parameter, O4 DSCDs at 360 and 477 nm and at the six off-zenith elevations were simulated with the parameter being tested set as different values, while all the other parameters were fixed. In this section, we only present the results of the sensitivity tests under the common settings listed in Table B1. In the following subsections, all the unmentioned simulation parameters followed the common settings. The extreme and median values of each parameter are also discussed in the following subsections.
B1 Sensitivity to surface albedo
It is difficult to estimate the surface albedo around the measurement site. In other studies, the surface albedo at low-altitude sites was usually estimated to be 0.05–0.1 (e.g., Irie et al., 2008; Ma et al., 2013; Wagner et al., 2011; Chan et al., 2017; Li et al., 2010, 2013; Clémer et al., 2010; Wang et al., 2016), while at a high-altitude site it was estimated to be 0.2 (Franco et al., 2015). As for the UFS, on the one hand, the snow cover and naked rocks are more reflective than typical urban and rural surfaces; on the other hand, the deep valleys close to the site can significantly decrease the surface albedo. In addition, the measurements at different elevation angles are sensitive to different parts of the surface. The effective surface albedo also depends on the observation geometry. The forming and melting of the snow cover can affect the surface albedo as well. However, the RTM can only assume a constant surface albedo. Therefore, we have to estimate a variation range of the surface albedo and consider the possible uncertainty in the retrieval. In this study, we empirically estimate that the surface albedo varies between 0.025 and 0.2 with a median value of 0.1 for both 360 and 477 nm.
In order to estimate the uncertainty of simulated O4 DSCD due to the surface albedo, we simulated O4 DSCDs with extreme surface albedo values (0.025 and 0.2) and the median value (0.1), while the other parameters were fixed as the settings listed in Table B1. Besides the box-shape profile with AOD =0.12, we also did a test with an aerosol-free profile. The relative differences of the O4 DSCDs simulated with extreme surface albedo values compared to those simulated with the median value are shown in Fig. B1.
The result shows that at both 360 and 477 nm, O4 DSCDs of all elevation angles slightly decrease with increasing surface albedo, and the variation rate differs with different elevation angles and different aerosol loads. Based on our estimation of the variation range of surface albedo, if the estimated median value (0.1) is used in the forward simulation, the uncertainty caused by the surface albedo assumption would be less than 3 %, and the positive and negative errors are nearly equal. Our further simulations also show that the uncertainty caused by surface albedo depends on the observation geometry. In the aerosol profile retrieval, we use a simple lookup table to determine the simulation error caused by surface albedo (see Sect. 3.7.2).
B2 Sensitivity to single-scattering albedo
As aerosol optical property data at the UFS are not available, we use the AERONET data at Hohenpeißenberg instead. According to the long-term data, for the single-scattering albedo (SSA) at 360 nm, 90 % of the data vary between 0.87 and 0.997, and the median value is 0.93; for the SSA at 477 nm, 90 % of the data vary between 0.85 and 0.997, and the median value is 0.92.
In order to estimate the uncertainty of simulated O4 DSCD due to the SSA, O4 DSCDs were simulated with the median and extreme SSA values (0.87, 0.93, and 0.997 for 360 nm; 0.85, 09.2, and 0.997 for 477 nm), while the other parameters were fixed as the settings listed in Table B1. The relative differences between the O4 DSCDs simulated with extreme and median SSA values are shown in Fig. B2.
The result indicates that using the median SSA in the forward simulation would result in less than 1 % error in O4 DSCDs in 90 % of the cases. In addition, the positive and negative errors are mostly equal. Although the measurements of SSA were taken at a much lower-altitude site, the sensitivity result shows the error attributed to SSA is rather small (<1 %). Therefore, using the SSA values from Hohenpeißenberg should not have a big influence on the retrieval. Since the simulation error caused by SSA can be influenced by the aerosol load as well as the observation geometry, it is determined using a simple lookup table in our aerosol profile retrieval (see Sect. 3.7.2).
B3 Sensitivity to scattering phase function
The estimation of the uncertainty of simulated O4 DSCD due to scattering phase function is also based on the AERONET data at Hohenpeißenberg. Unlike most of the other simulation parameters which can be defined by a single number, the parameter of scattering phase function is defined by function values at different scattering angles. In order to estimate the uncertainty, we simulated O4 DSCDs with all the phase functions from 2013 to 2014 (altogether 179 available data), while the other parameters were fixed as the settings listed in Table B1. The frequency distributions of simulated O4 DSCDs are shown in Fig. B3. For each elevation angle, the percentage standard deviation is indicated beside the curve, and the grey dashed line shows the median value. The results indicate that the distribution of the simulated O4 DSCDs follows the normal distribution, and the standard deviation at 477 nm is larger compared to that at 360 nm.
Based on the simulation results, the phase function with which the simulated O4 DSCDs at all elevation angles are closest to the median values is chosen as the so-called median phase function for each wavelength. In our aerosol profile retrieval, the error caused by scattering phase function is also determined using a simple lookup table (see Sect. 3.7.2).
B4 Sensitivity to aerosols at different altitudes
The sensitivity of O4 DSCD to aerosol extinction at different altitudes was estimated by simulating O4 DSCDs with box-shape aerosol profiles with the same aerosol extinction coefficient of 0.04 km−1 and different box heights varying from 0 to 8 km. The other parameters were fixed as the settings listed in Table B1. Figure B4 shows the simulated O4 DSCDs at 360 and 477 nm for each elevation angle. The result indicates that the sensitivities of O4 DSCDs at all elevation angles decrease rapidly with increasing box height (and also increasing AOD). Furthermore, O4 DSCDs at all elevation angles are almost constant when the box height varies between 2 and 8 km, which indicates that O4 absorption is almost insensitive to the aerosols above 2 km. Taking the O4 DSCD measured at 360 nm with an elevation angle of 2∘ as an example, the sensitivity to aerosols at 2 km is lower than that at surface level by a factor of ∼40. In addition, measurements at lower elevation angles are more sensitive to aerosols close to the surface compared to those at higher elevations. According to the result, our retrieval of aerosol profiles would only focus on aerosols below 2 km.
This result coincides with the results reported in the MAX-DOAS studies based on the OEM (e.g., Frieß et al., 2006, 2016; Clémer et al., 2010; Bösch et al., 2018). In these studies, the averaging kernels – which indicate the measurement sensitivity to aerosols at different altitudes – are all close to zero at the altitudes above 2 km.
B5 Sensitivity to aerosols above retrieval height
As discussed in Sect. B4, our aerosol profile retrieval would only focus on aerosols below 2 km. However, the aerosol load on Zugspitze is usually very low and the aerosol extinction coefficient above 2 km is usually of the same order of magnitude as the one below 2 km. We estimate that the aerosol extinction coefficient between 2 and 4 km (denoted as AE2–4) varies from 0 % to 100 % of the aerosol below 2 km (denoted as AE0–2), and the median value is 50 % of AE0–2. In order to estimate the sensitivity of O4 absorption to AE2–4, O4 DSCDs were simulated with profiles with the same aerosol extinction coefficient below 2 km (AE0–2=0.04 km−1), and AE2–4 equals 0 %, 50 %, and 100 % of AE0–2, and the other parameters were fixed as the settings listed in Table B1. The differences between the O4 DSCDs simulated with extreme and median AE2–4 are shown in Fig. B5. The result indicates that the aerosols above 2 km can affect the O4 DSCDs by up to ∼3 %, which is similar to the surface albedo. Therefore, we consider the influence from the aerosols above 2 km to be a kind of measurement uncertainty and treat it in the same way as the errors caused by surface albedo, single-scattering albedo, and phase function uncertainties. Similarly, the uncertainty caused by aerosols above retrieval height is determined using a simple lookup table in our profile retrieval (see Sect. 3.7.2).
B6 Sensitivity to surface aerosol extinction
In order to estimate the sensitivity of O4 DSCD to surface aerosol extinction, O4 DSCDs were simulated with box-shape profiles with the constant box height of 2 km and aerosol extinction coefficient vary from 0 to 1 km−1. The other parameters were fixed as the settings listed in Table B1. The simulated O4 DSCDs are shown in Fig. B6. The result indicates that the sensitivities of O4 absorption at all elevation angles and both wavelength bands decrease with increasing aerosol extinction. Furthermore, the sensitivity is very low when the surface aerosol extinction coefficient exceeds 0.3 km−1. The O4 DSCDs at all elevation angles and both wavelengths decrease monotonically with increasing aerosol extinction. In addition, measurements at lower elevation angles are much more sensitive.
The MAX-DOAS data are available upon request from the corresponding author (email@example.com).
ZW and KLC designed the retrieval algorithm. ZW, KLC, and RH contributed to the measurement and instrument maintenance. KLC, KPH, AD, TW, and MW provided useful comments for the discussion and support for the article writing. TW provided support for the radiative transfer calculation. MW and KLC provided auxiliary ceilometer data. ZW analyzed the measurement data and prepared the article, with contributions from all co-authors.
The authors declare that they have no conflict of interest.
This work is funded by the DLR-DAAD Research Fellowships 2014 (50019750) program with reference number 91549461. We are thankful for the help from the colleagues of the University of Heidelberg in installing the MAX-DOAS instrument at the UFS. We thank Ina Mattis and her colleagues at DWD for their effort in establishing and maintaining the AERONET site in Hohenpeißenberg. The authors would also like to thank the Royal Belgian Institute for Space Aeronomy (BIRA-IASB) for the provision of the QDOAS software used in this study. We are grateful to Till Rehm and his colleagues at the UFS for the maintenance work. We thank the Bayerisches Umweltministerium for supporting the UFS as part of their mission.
This research has been supported by the DLR-DAAD Research Fellowships (grant no. 50019750).
The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.
This paper was edited by Alexander Kokhanovsky and reviewed by three anonymous referees.
Aliwell, S. R., Roozendael, M. V., Johnston, P. V., Richter, A., Wagner, T., Arlander, D. W., Burrows, J. P., Fish, D. J., Jones, R. L., Tørnkvist, K. K., Lambert, J.-C., Pfeilsticker, K., and Pundt, I.: Analysis for BrO in zenith-sky spectra: An intercomparison exercise for analysis improvement, J. Geophys. Res.-Atmos., 107, ACH 10-1–ACH 10-20, https://doi.org/10.1029/2001JD000329, 2002. a, b
Almeida-Silva, M., Almeida, S. M., Freitas, M. C., Pio, C. A., Nunes, T., and Cardoso, J.: Impact of Sahara Dust Transport on Cape Verde Atmospheric Element Particles, J. Toxicol. Env. Heal. A, 76, 240–251, https://doi.org/10.1080/15287394.2013.757200, 2013. a
Anderson, G. P., Clough, S. A., Kneizys, F., Chetwynd, J. H., and Shettle, E. P.: AFGL atmospheric constituent profiles (0.120 km), Tech. rep., AIR FORCE GEOPHYSICS LAB HANSCOM AFB MA, 1986. a
Bäumer, D., Vogel, B., Versick, S., Rinke, R., Möhler, O., and Schnaiter, M.: Relationship of visibility, aerosol optical thickness and aerosol size distribution in an ageing air mass over South-West Germany, Atmos. Environ., 42, 989–998, https://doi.org/10.1016/j.atmosenv.2007.10.017, 2008. a
Bellouin, N., Boucher, O., Haywood, J., and Reddy, M. S.: Global estimate of aerosol direct radiative forcing from satellite measurements, Nature, 438, 1138–1141, https://doi.org/10.1038/nature04348, 2005. a
Bogumil, K., Orphal, J., Homann, T., Voigt, S., Spietz, P., Fleischmann, O., Vogel, A., Hartmann, M., Kromminga, H., Bovensmann, H., Frerick, J., and Burrows, J.: Measurements of molecular absorption spectra with the SCIAMACHY pre-flight model: instrument characterization and reference data for atmospheric remote-sensing in the 230–2380 nm region, J. Photoch. Photobio. A, 157, 167–184, https://doi.org/10.1016/S1010-6030(03)00062-5, 2003. a, b
Bösch, T., Rozanov, V., Richter, A., Peters, E., Rozanov, A., Wittrock, F., Merlaud, A., Lampel, J., Schmitt, S., de Haij, M., Berkhout, S., Henzing, B., Apituley, A., den Hoed, M., Vonk, J., Tiefengraber, M., Müller, M., and Burrows, J. P.: BOREAS – a new MAX-DOAS profile retrieval algorithm for aerosols and trace gases, Atmos. Meas. Tech., 11, 6833–6859, https://doi.org/10.5194/amt-11-6833-2018, 2018. a
Brook, R. D., Rajagopalan, S., Pope, C. A., Brook, J. R., Bhatnagar, A., Diez-Roux, A. V., Holguin, F., Hong, Y., Luepker, R. V., Mittleman, M. A., Peters, A., Siscovick, D., Smith, S. C., Whitsel, L., and Kaufman, J. D.: Particulate Matter Air Pollution and Cardiovascular Disease, Circulation, 121, 2331–2378, https://doi.org/10.1161/CIR.0b013e3181dbece1, 2010. a
Cazorla, A., Casquero-Vera, J. A., Román, R., Guerrero-Rascado, J. L., Toledano, C., Cachorro, V. E., Orza, J. A. G., Cancillo, M. L., Serrano, A., Titos, G., Pandolfi, M., Alastuey, A., Hanrieder, N., and Alados-Arboledas, L.: Near-real-time processing of a ceilometer network assisted with sun-photometer data: monitoring a dust outbreak over the Iberian Peninsula, Atmos. Chem. Phys., 17, 11861–11876, https://doi.org/10.5194/acp-17-11861-2017, 2017. a
Chan, K., Wiegner, M., Wenig, M., and Pöhler, D.: Observations of tropospheric aerosols and NO2 in Hong Kong over 5 years using ground based MAX-DOAS, Sci. Total Environ., 619–620, 1545–1556, https://doi.org/10.1016/j.scitotenv.2017.10.153, 2017. a, b, c, d, e, f, g, h
Chan, K. L., Hartl, A., Lam, Y. F., Xie, P. H., Liu, W. Q., Cheung, H. M., Lampel, J., Pöhler, D., Li, A., Xu, J., Zhou, H. J., Ning, Z., and Wenig, M. O.: Observations of tropospheric NO2 using ground based MAX-DOAS and OMI measurements during the Shanghai World Expo 2010, Atmos. Environ., 119, 45–58, https://doi.org/10.1016/j.atmosenv.2015.08.041, 2015. a, b, c
Chan, K. L., Wiegner, M., Flentje, H., Mattis, I., Wagner, F., Gasteiger, J., and Geiß, A.: Evaluation of ECMWF-IFS (version 41R1) operational model forecasts of aerosol transport by using ceilometer network measurements, Geosci. Model Dev., 11, 3807–3831, https://doi.org/10.5194/gmd-11-3807-2018, 2018. a, b
Chan, K. L., Wang, Z., Ding, A., Heue, K.-P., Shen, Y., Wang, J., Zhang, F., Shi, Y., Hao, N., and Wenig, M.: MAX-DOAS measurements of tropospheric NO2 and HCHO in Nanjing and a comparison to ozone monitoring instrument observations, Atmos. Chem. Phys., 19, 10051–10071, https://doi.org/10.5194/acp-19-10051-2019, 2019. a
Chance, K. and Kurucz, R.: An improved high-resolution solar reference spectrum for earth's atmosphere measurements in the ultraviolet, visible, and near infrared, J. Quant. Spectrosc. Ra., 111, 1289–1295, https://doi.org/10.1016/j.jqsrt.2010.01.036, 2010. a
Chance, K. and Orphal, J.: Revised ultraviolet absorption cross sections of H2CO for the HITRAN database, J. Quant. Spectrosc. Ra., 112, 1509–1510, https://doi.org/10.1016/j.jqsrt.2011.02.002, 2011. a
Chance, K. V. and Spurr, R. J. D.: Ring effect studies: Rayleigh scattering, including molecular parameters for rotational Raman scattering, and the Fraunhofer spectrum, Appl. Optics, 36, 5224–5230, https://doi.org/10.1364/AO.36.005224, 1997. a
Clémer, K., Van Roozendael, M., Fayt, C., Hendrick, F., Hermans, C., Pinardi, G., Spurr, R., Wang, P., and De Mazière, M.: Multiple wavelength retrieval of tropospheric aerosol optical properties from MAXDOAS measurements in Beijing, Atmos. Meas. Tech., 3, 863–878, https://doi.org/10.5194/amt-3-863-2010, 2010. a, b, c, d, e, f, g, h, i, j, k, l
Deutschmann, T.: Atmospheric radiative transfer modelling with Monte Carlo methods, Institute of environmental physics University of Heidelberg, p. 3935, 2008. a
Dubovik, O., Holben, B., Eck, T. F., Smirnov, A., Kaufman, Y. J., King, M. D., Tanré, D., and Slutsker, I.: Variability of Absorption and Optical Properties of Key Aerosol Types Observed in Worldwide Locations, J. Atmos. Sci., 59, 590–608, https://doi.org/10.1175/1520-0469(2002)059<0590:VOAAOP>2.0.CO;2, 2002. a
Emeis, S., Forkel, R., Junkermann, W., Schäfer, K., Flentje, H., Gilge, S., Fricke, W., Wiegner, M., Freudenthaler, V., Groβ, S., Ries, L., Meinhardt, F., Birmili, W., Münkel, C., Obleitner, F., and Suppan, P.: Measurement and simulation of the 16/17 April 2010 Eyjafjallajökull volcanic ash layer dispersion in the northern Alpine region, Atmos. Chem. Phys., 11, 2689–2701, https://doi.org/10.5194/acp-11-2689-2011, 2011. a
Franco, B., Hendrick, F., Van Roozendael, M., Müller, J.-F., Stavrakou, T., Marais, E. A., Bovy, B., Bader, W., Fayt, C., Hermans, C., Lejeune, B., Pinardi, G., Servais, C., and Mahieu, E.: Retrievals of formaldehyde from ground-based FTIR and MAX-DOAS observations at the Jungfraujoch station and comparisons with GEOS-Chem and IMAGES model simulations, Atmos. Meas. Tech., 8, 1733–1756, https://doi.org/10.5194/amt-8-1733-2015, 2015. a
Frieß, U., Monks, P. S., Remedios, J. J., Rozanov, A., Sinreich, R., Wagner, T., and Platt, U.: MAX-DOAS O4 measurements: A new technique to derive information on atmospheric aerosols: 2. Modeling studies, J. Geophys. Res.-Atmos., 111, D14, https://doi.org/10.1029/2005JD006618, 2006. a, b, c
Frieß, U., Sihler, H., Sander, R., Pöhler, D., Yilmaz, S., and Platt, U.: The vertical distribution of BrO and aerosols in the Arctic: Measurements by active and passive differential optical absorption spectroscopy, J. Geophys. Res.-Atmos., 116, D14, https://doi.org/10.1029/2011JD015938, 2011. a, b, c
Frieß, U., Klein Baltink, H., Beirle, S., Clémer, K., Hendrick, F., Henzing, B., Irie, H., de Leeuw, G., Li, A., Moerman, M. M., van Roozendael, M., Shaiganfar, R., Wagner, T., Wang, Y., Xie, P., Yilmaz, S., and Zieger, P.: Intercomparison of aerosol extinction profiles retrieved from MAX-DOAS measurements, Atmos. Meas. Tech., 9, 3205–3222, https://doi.org/10.5194/amt-9-3205-2016, 2016. a, b
Gielen, C., Van Roozendael, M., Hendrick, F., Pinardi, G., Vlemmix, T., De Bock, V., De Backer, H., Fayt, C., Hermans, C., Gillotay, D., and Wang, P.: A simple and versatile cloud-screening method for MAX-DOAS retrievals, Atmos. Meas. Tech., 7, 3509–3527, https://doi.org/10.5194/amt-7-3509-2014, 2014. a, b
Halla, J. D., Wagner, T., Beirle, S., Brook, J. R., Hayden, K. L., O'Brien, J. M., Ng, A., Majonis, D., Wenig, M. O., and McLaren, R.: Determination of tropospheric vertical columns of NO2 and aerosol optical properties in a rural setting using MAX-DOAS, Atmos. Chem. Phys., 11, 12475–12498, https://doi.org/10.5194/acp-11-12475-2011, 2011. a
Hartl, A. and Wenig, M. O.: Regularisation model study for the least-squares retrieval of aerosol extinction time series from UV/VIS MAX-DOAS observations for a ground layer profile parameterisation, Atmos. Meas. Tech., 6, 1959–1980, https://doi.org/10.5194/amt-6-1959-2013, 2013. a
Heald, C. L., Ridley, D. A., Kroll, J. H., Barrett, S. R. H., Cady-Pereira, K. E., Alvarado, M. J., and Holmes, C. D.: Contrasting the direct radiative effect and direct radiative forcing of aerosols, Atmos. Chem. Phys., 14, 5513–5527, https://doi.org/10.5194/acp-14-5513-2014, 2014. a
Hendrick, F., Müller, J.-F., Clémer, K., Wang, P., De Mazière, M., Fayt, C., Gielen, C., Hermans, C., Ma, J. Z., Pinardi, G., Stavrakou, T., Vlemmix, T., and Van Roozendael, M.: Four years of ground-based MAX-DOAS observations of HONO and NO2 in the Beijing area, Atmos. Chem. Phys., 14, 765–781, https://doi.org/10.5194/acp-14-765-2014, 2014. a
Hinds, W. C.: Aerosol technology: properties, behavior, and measurement of airborne particles, John Wiley & Sons, 2012. a
Holben, B., Eck, T., Slutsker, I., Tanre, D., Buis, J., Setzer, A., Vermote, E., Reagan, J., Kaufman, Y., Nakajima, T., Lavenu, F., Jankowiak, I., and Smirnov, A.: AERONET – A Federated Instrument Network and Data Archive for Aerosol Characterization, Remote Sens. Environ., 66, 1–16, https://doi.org/10.1016/S0034-4257(98)00031-5, 1998. a
Irie, H., Kanaya, Y., Akimoto, H., Iwabuchi, H., Shimizu, A., and Aoki, K.: First retrieval of tropospheric aerosol profiles using MAX-DOAS and comparison with lidar and sky radiometer measurements, Atmos. Chem. Phys., 8, 341–350, https://doi.org/10.5194/acp-8-341-2008, 2008. a, b
Irie, H., Takashima, H., Kanaya, Y., Boersma, K. F., Gast, L., Wittrock, F., Brunner, D., Zhou, Y., and Van Roozendael, M.: Eight-component retrievals from ground-based MAX-DOAS observations, Atmos. Meas. Tech., 4, 1027–1044, https://doi.org/10.5194/amt-4-1027-2011, 2011. a, b
Irie, H., Nakayama, T., Shimizu, A., Yamazaki, A., Nagai, T., Uchiyama, A., Zaizen, Y., Kagamitani, S., and Matsumi, Y.: Evaluation of MAX-DOAS aerosol retrievals by coincident observations using CRDS, lidar, and sky radiometer inTsukuba, Japan, Atmos. Meas. Tech., 8, 2775–2788, https://doi.org/10.5194/amt-8-2775-2015, 2015. a
Jin, J., Ma, J., Lin, W., Zhao, H., Shaiganfar, R., Beirle, S., and Wagner, T.: MAX-DOAS measurements and satellite validation of tropospheric NO2 and SO2 vertical column densities at a rural site of North China, Atmos. Environ., 133, 12–25, https://doi.org/10.1016/j.atmosenv.2016.03.031, 2016. a
Junker, C. and Liousse, C.: A global emission inventory of carbonaceous aerosol from historic records of fossil fuel and biofuel consumption for the period 1860–1997, Atmos. Chem. Phys., 8, 1195–1207, https://doi.org/10.5194/acp-8-1195-2008, 2008. a
Karanasiou, A., Moreno, N., Moreno, T., Viana, M., de Leeuw, F., and Querol, X.: Health effects from Sahara dust episodes in Europe: Literature review and research gaps, Environ. Int., 47, 107–114, https://doi.org/10.1016/j.envint.2012.06.012, 2012. a
Kreher, K., Van Roozendael, M., Hendrick, F., Apituley, A., Dimitropoulou, E., Frieß, U., Richter, A., Wagner, T., Abuhassan, N., Ang, L., Anguas, M., Bais, A., Benavent, N., Bösch, T., Bognar, K., Borovski, A., Bruchkouski, I., Cede, A., Chan, K. L., Donner, S., Drosoglou, T., Fayt, C., Finkenzeller, H., Garcia-Nieto, D., Gielen, C., Gómez-Martín, L., Hao, N., Herman, J. R., Hermans, C., Hoque, S., Irie, H., Jin, J., Johnston, P., Khayyam Butt, J., Khokhar, F., Koenig, T. K., Kuhn, J., Kumar, V., Lampel, J., Liu, C., Ma, J., Merlaud, A., Mishra, A. K., Müller, M., Navarro-Comas, M., Ostendorf, M., Pazmino, A., Peters, E., Pinardi, G., Pinharanda, M., Piters, A., Platt, U., Postylyakov, O., Prados-Roman, C., Puentedura, O., Querel, R., Saiz-Lopez, A., Schönhardt, A., Schreier, S. F., Seyler, A., Sinha, V., Spinei, E., Strong, K., Tack, F., Tian, X., Tiefengraber, M., Tirpitz, J.-L., van Gent, J., Volkamer, R., Vrekoussis, M., Wang, S., Wang, Z., Wenig, M., Wittrock, F., Xie, P. H., Xu, J., Yela, M., Zhang, C., and Zhao, X.: Intercomparison of NO2, O4, O3 and HCHO slant column measurements by MAX-DOAS and zenith-sky UV-Visible spectrometers during the CINDI-2 campaign, Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2019-157, in review, 2019. a, b
Lee, H., Irie, H., Kim, Y. J., Noh, Y., Lee, C., Kim, Y., and Chun, K. J.: Retrieval of Aerosol Extinction in the Lower Troposphere Based on UV MAX-DOAS Measurements, Aerosol Sci. Tech., 43, 502–509, https://doi.org/10.1080/02786820902769691, 2009. a
Lee, Y. C., Lam, Y. F., Kuhlmann, G., Wenig, M. O., Chan, K. L., Hartl, A., and Ning, Z.: An integrated approach to identify the biomass burning sources contributing to black carbon episodes in Hong Kong, Atmos. Environ., 80, 478–487, https://doi.org/10.1016/j.atmosenv.2013.08.030, 2013. a
Levy II, H., Horowitz, L. W., Schwarzkopf, M. D., Ming, Y., Golaz, J.-C., Naik, V., and Ramaswamy, V.: The roles of aerosol direct and indirect effects in past and future climate change, J. Geophys. Res.-Atmos., 118, 4521–4532, https://doi.org/10.1002/jgrd.50192, 2013. a
Li, X., Brauers, T., Shao, M., Garland, R. M., Wagner, T., Deutschmann, T., and Wahner, A.: MAX-DOAS measurements in southern China: retrieval of aerosol extinctions and validation using ground-based in-situ data, Atmos. Chem. Phys., 10, 2079–2089, https://doi.org/10.5194/acp-10-2079-2010, 2010. a, b, c
Li, X., Brauers, T., Hofzumahaus, A., Lu, K., Li, Y. P., Shao, M., Wagner, T., and Wahner, A.: MAX-DOAS measurements of NO2, HCHO and CHOCHO at a rural site in Southern China, Atmos. Chem. Phys., 13, 2133–2151, https://doi.org/10.5194/acp-13-2133-2013, 2013. a, b
Liu, S. C., McKeen, S. A., and Madronich, S.: Effect of anthropogenic aerosols on biologically active ultraviolet radiation, Geophys. Res. Lett., 18, 2265–2268, https://doi.org/10.1029/91GL02773, 1991. a
Ma, J. Z., Beirle, S., Jin, J. L., Shaiganfar, R., Yan, P., and Wagner, T.: Tropospheric NO2 vertical column densities over Beijing: results of the first three years of ground-based MAX-DOAS measurements (2008–2011) and satellite validation, Atmos. Chem. Phys., 13, 1547–1567, https://doi.org/10.5194/acp-13-1547-2013, 2013. a, b
Pappalardo, G., Amodeo, A., Apituley, A., Comeron, A., Freudenthaler, V., Linné, H., Ansmann, A., Bösenberg, J., D'Amico, G., Mattis, I., Mona, L., Wandinger, U., Amiridis, V., Alados-Arboledas, L., Nicolae, D., and Wiegner, M.: EARLINET: towards an advanced sustainable European aerosol lidar network, Atmos. Meas. Tech., 7, 2389–2409, https://doi.org/10.5194/amt-7-2389-2014, 2014. a
Platt, U., Perner, D., and Patz, H. W.: Simultaneous measurement of atmospheric CH2O, O3, and NO2 by differential optical absorption, J. Geophys. Res.-Oceans, 84, 6329–6335, https://doi.org/10.1029/JC084iC10p06329, 1979. a
Risius, S., Xu, H., Di Lorenzo, F., Xi, H., Siebert, H., Shaw, R. A., and Bodenschatz, E.: Schneefernerhaus as a mountain research station for clouds and turbulence, Atmos. Meas. Tech., 8, 3209–3218, https://doi.org/10.5194/amt-8-3209-2015, 2015. a
Rodgers, C. D.: Inverse methods for atmospheric sounding: Theory and practice, vol. 2, World scientific, 2000. a
Rothman, L., Gordon, I., Barber, R., Dothe, H., Gamache, R., Goldman, A., Perevalov, V., Tashkun, S., and Tennyson, J.: HITEMP, the high-temperature molecular spectroscopic database, J. Quant. Spectrosc. Ra., 111, 2139–2150, https://doi.org/10.1016/j.jqsrt.2010.05.001, 2010. a
Schäfer, K., Thomas, W., Peters, A., Ries, L., Obleitner, F., Schnelle-Kreis, J., Birmili, W., Diemer, J., Fricke, W., Junkermann, W., Pitz, M., Emeis, S., Forkel, R., Suppan, P., Flentje, H., Gilge, S., Wichmann, H. E., Meinhardt, F., Zimmermann, R., Weinhold, K., Soentgen, J., Münkel, C., Freuer, C., and Cyrys, J.: Influences of the 2010 Eyjafjallajökull volcanic plume on air quality in the northern Alpine region, Atmos. Chem. Phys., 11, 8555–8575, https://doi.org/10.5194/acp-11-8555-2011, 2011. a
Serdyuchenko, A., Gorshelev, V., Weber, M., Chehade, W., and Burrows, J. P.: High spectral resolution ozone absorption cross-sections – Part 2: Temperature dependence, Atmos. Meas. Tech., 7, 625–636, https://doi.org/10.5194/amt-7-625-2014, 2014. a, b
Sinreich, R., Merten, A., Molina, L., and Volkamer, R.: Parameterizing radiative transfer to convert MAX-DOAS dSCDs into near-surface box-averaged mixing ratios, Atmos. Meas. Tech., 6, 1521–1532, https://doi.org/10.5194/amt-6-1521-2013, 2013. a
Spurr, R.: LIDORT and VLIDORT: linearized pseudo-spherical scalar and vector discrete ordinate radiative transfer models for use in remote sensing retrieval problems, Light scattering reviews, 3, 229–275, 2008. a
Spurr, R., Kurosu, T., and Chance, K.: A linearized discrete ordinate radiative transfer model for atmospheric remote-sensing retrieval, J. Quant. Spectrosc. Ra., 68, 689–735, https://doi.org/10.1016/S0022-4073(00)00055-8, 2001. a
Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, B., and Midgley, B.: IPCC, 2013: climate change 2013: the physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change, 2013. a
Thalman, R. and Volkamer, R.: Temperature dependent absorption cross-sections of O2−O2 collision pairs between 340 and 630 nm and at atmospherically relevant pressure, Phys. Chem. Chem. Phys., 15, 15371–15381, https://doi.org/10.1039/C3CP50968K, 2013. a, b
Tirpitz, J.-L., Frieß, U., Hendrick, F., Alberti, C., Allaart, M., Apituley, A., Bais, A., Beirle, S., Berkhout, S., Bognar, K., Bösch, T., Bruchkouski, I., Cede, A., Chan, K. L., den Hoed, M., Donner, S., Drosoglou, T., Fayt, C., Friedrich, M. M., Frumau, A., Gast, L., Gielen, C., Gomez-Martín, L., Hao, N., Hensen, A., Henzing, B., Hermans, C., Jin, J., Kreher, K., Kuhn, J., Lampel, J., Li, A., Liu, C., Liu, H., Ma, J., Merlaud, A., Peters, E., Pinardi, G., Piters, A., Platt, U., Puentedura, O., Richter, A., Schmitt, S., Spinei, E., Stein Zweers, D., Strong, K., Swart, D., Tack, F., Tiefengraber, M., van der Hoff, R., van Roozendael, M., Vlemmix, T., Vonk, J., Wagner, T., Wang, Y., Wang, Z., Wenig, M., Wiegner, M., Wittrock, F., Xie, P., Xing, C., Xu, J., Yela, M., Zhang, C., and Zhao, X.: Intercomparison of MAX-DOAS vertical profile retrieval algorithms: studies on field data from the CINDI-2 campaign, Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2019-456, in review, 2020. a
Toledano, C., Wiegner, M., Garhammer, M., Seefeldner, M., Gasteiger, J., Müller, D., and Koepke, P.: Spectral aerosol optical depth characterization of desert dust during SAMUM 2006, Tellus B, 61, 216–228, https://doi.org/10.1111/j.1600-0889.2008.00382.x, 2009. a
Toledano, C., Wiegner, M., Groß, S., Freudenthaler, V., Gasteiger, J., Müller, D., Müller, T., Schladitz, A., Weinzierl, B., Torres, B., and O'neill, N. T.: Optical properties of aerosol mixtures derived from sun-sky radiometry during SAMUM-2, Tellus B, 63, 635–648, https://doi.org/10.1111/j.1600-0889.2011.00573.x, 2011. a
Valavanidis, A., Fiotakis, K., and Vlachogianni, T.: Airborne Particulate Matter and Human Health: Toxicological Assessment and Importance of Size and Composition of Particles for Oxidative Damage and Carcinogenic Mechanisms, J. Environ. Sci. Heal. C, 26, 339–362, https://doi.org/10.1080/10590500802494538, 2008. a
Viana, M., Pey, J., Querol, X., Alastuey, A., de Leeuw, F., and Lükewille, A.: Natural sources of atmospheric aerosols influencing air quality across Europe, Sci. Total Environ., 472, 825–833, https://doi.org/10.1016/j.scitotenv.2013.11.140, 2014. a
Vlemmix, T., Piters, A. J. M., Berkhout, A. J. C., Gast, L. F. L., Wang, P., and Levelt, P. F.: Ability of the MAX-DOAS method to derive profile information for NO2: can the boundary layer and free troposphere be separated?, Atmos. Meas. Tech., 4, 2659–2684, https://doi.org/10.5194/amt-4-2659-2011, 2011. a, b
Volkamer, R., Spietz, P., Burrows, J., and Platt, U.: High-resolution absorption cross-section of glyoxal in the UV-Vis and IR spectral ranges, J. Photoch. Photobio. A, 172, 35–46, https://doi.org/10.1016/j.jphotochem.2004.11.011, 2005. a
Wagner, T., Dix, B., Friedeburg, C. v., Fries, U., Sanghavi, S., Sinreich, R., and Platt, U.: MAX-DOAS O4 measurements: A new technique to derive information on atmospheric aerosols – Principles and information content, J. Geophys. Res.-Atmos., 109, D22, https://doi.org/10.1029/2004JD004904, 2004. a
Wagner, T., Burrows, J. P., Deutschmann, T., Dix, B., von Friedeburg, C., Frieß, U., Hendrick, F., Heue, K.-P., Irie, H., Iwabuchi, H., Kanaya, Y., Keller, J., McLinden, C. A., Oetjen, H., Palazzi, E., Petritoli, A., Platt, U., Postylyakov, O., Pukite, J., Richter, A., van Roozendael, M., Rozanov, A., Rozanov, V., Sinreich, R., Sanghavi, S., and Wittrock, F.: Comparison of box-air-mass-factors and radiances for Multiple-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) geometries calculated from different UV/visible radiative transfer models, Atmos. Chem. Phys., 7, 1809–1833, https://doi.org/10.5194/acp-7-1809-2007, 2007. a, b
Wagner, T., Deutschmann, T., and Platt, U.: Determination of aerosol properties from MAX-DOAS observations of the Ring effect, Atmos. Meas. Tech., 2, 495–512, https://doi.org/10.5194/amt-2-495-2009, 2009. a, b
Wagner, T., Beirle, S., Brauers, T., Deutschmann, T., Frieß, U., Hak, C., Halla, J. D., Heue, K. P., Junkermann, W., Li, X., Platt, U., and Pundt-Gruber, I.: Inversion of tropospheric profiles of aerosol extinction and HCHO and NO2 mixing ratios from MAX-DOAS observations in Milano during the summer of 2003 and comparison with independent data sets, Atmos. Meas. Tech., 4, 2685–2715, https://doi.org/10.5194/amt-4-2685-2011, 2011. a, b, c, d
Wagner, T., Apituley, A., Beirle, S., Dörner, S., Friess, U., Remmers, J., and Shaiganfar, R.: Cloud detection and classification based on MAX-DOAS observations, Atmos. Meas. Tech., 7, 1289–1320, https://doi.org/10.5194/amt-7-1289-2014, 2014. a
Wagner, T., Beirle, S., Remmers, J., Shaiganfar, R., and Wang, Y.: Absolute calibration of the colour index and O4 absorption derived from Multi AXis (MAX-)DOAS measurements and their application to a standardised cloud classification algorithm, Atmos. Meas. Tech., 9, 4803–4823, https://doi.org/10.5194/amt-9-4803-2016, 2016. a, b, c
Wagner, T., Beirle, S., Benavent, N., Bösch, T., Chan, K. L., Donner, S., Dörner, S., Fayt, C., Frieß, U., García-Nieto, D., Gielen, C., González-Bartolome, D., Gomez, L., Hendrick, F., Henzing, B., Jin, J. L., Lampel, J., Ma, J., Mies, K., Navarro, M., Peters, E., Pinardi, G., Puentedura, O., Puķīte, J., Remmers, J., Richter, A., Saiz-Lopez, A., Shaiganfar, R., Sihler, H., Van Roozendael, M., Wang, Y., and Yela, M.: Is a scaling factor required to obtain closure between measured and modelled atmospheric O4 absorptions? An assessment of uncertainties of measurements and radiative transfer simulations for 2 selected days during the MAD-CAT campaign, Atmos. Meas. Tech., 12, 2745–2817, https://doi.org/10.5194/amt-12-2745-2019, 2019. a, b, c, d
Wang, S., Cuevas, C. A., Frieß, U., and Saiz-Lopez, A.: MAX-DOAS retrieval of aerosol extinction properties in Madrid, Spain, Atmos. Meas. Tech., 9, 5089–5101, https://doi.org/10.5194/amt-9-5089-2016, 2016. a, b, c, d, e, f
Wang, T., Hendrick, F., Wang, P., Tang, G., Clémer, K., Yu, H., Fayt, C., Hermans, C., Gielen, C., Müller, J.-F., Pinardi, G., Theys, N., Brenot, H., and Van Roozendael, M.: Evaluation of tropospheric SO2 retrieved from MAX-DOAS measurements in Xianghe, China, Atmos. Chem. Phys., 14, 11149–11164, https://doi.org/10.5194/acp-14-11149-2014, 2014a. a, b, c, d
Wang, Y., Sartelet, K. N., Bocquet, M., Chazette, P., Sicard, M., D'Amico, G., Léon, J. F., Alados-Arboledas, L., Amodeo, A., Augustin, P., Bach, J., Belegante, L., Binietoglou, I., Bush, X., Comerón, A., Delbarre, H., García-Vízcaino, D., Guerrero-Rascado, J. L., Hervo, M., Iarlori, M., Kokkalis, P., Lange, D., Molero, F., Montoux, N., Muñoz, A., Muñoz, C., Nicolae, D., Papayannis, A., Pappalardo, G., Preissler, J., Rizi, V., Rocadenbosch, F., Sellegri, K., Wagner, F., and Dulac, F.: Assimilation of lidar signals: application to aerosol forecasting in the western Mediterranean basin, Atmos. Chem. Phys., 14, 12031–12053, https://doi.org/10.5194/acp-14-12031-2014, 2014b. a
Warren, E., Charlton-Perez, C., Kotthaus, S., Lean, H., Ballard, S., Hopkin, E., and Grimmond, S.: Evaluation of forward-modelled attenuated backscatter using an urban ceilometer network in London under clear-sky conditions, Atmos. Environ., 191, 532–547, https://doi.org/10.1016/j.atmosenv.2018.04.045, 2018. a
Wiegner, M., Groß, S., Freudenthaler, V., Schnell, F., and Gasteiger, J.: The May/June 2008 Saharan dust event over Munich: Intensive aerosol parameters from lidar measurements, J. Geophys. Res.-Atmos., 116, d23213, https://doi.org/10.1029/2011JD016619, 2011. a
Wiegner, M., Madonna, F., Binietoglou, I., Forkel, R., Gasteiger, J., Geiß, A., Pappalardo, G., Schäfer, K., and Thomas, W.: What is the benefit of ceilometers for aerosol remote sensing? An answer from EARLINET, Atmos. Meas. Tech., 7, 1979–1997, https://doi.org/10.5194/amt-7-1979-2014, 2014. a, b
Xin, J., Wang, Y., Li, Z., Wang, P., Hao, W. M., Nordgren, B. L., Wang, S., Liu, G., Wang, L., Wen, T., Sun, Y., and Hu, B.: Aerosol optical depth (AOD) and Ångström exponent of aerosols observed by the Chinese Sun Hazemeter Network from August 2004 to September 2005, J. Geophys. Res.-Atmos., 112, D5, https://doi.org/10.1029/2006JD007075, 2007. a
Zhang, J., Wang, S., Guo, Y., Zhang, R., Qin, X., Huang, K., Wang, D., Fu, Q., Wang, J., and Zhou, B.: Aerosol vertical profile retrieved from ground-based MAX-DOAS observation and characteristic distribution during wintertime in Shanghai, China, Atmos. Environ., 192, 193–205, https://doi.org/10.1016/j.atmosenv.2018.08.051, 2018. a
Zhang, Z., Wenig, M., Zhou, W., Diehl, T., Chan, K. L., and Wang, L.: The contribution of different aerosol sources to the Aerosol Optical Depth in Hong Kong, Atmos. Environ., 83, 145–154, https://doi.org/10.1016/j.atmosenv.2013.10.047, 2014. a