A new method to determine the aerosol optical properties from multiple wavelength O 4 absorptions by MAX-DOAS observation

Ground based Multi-AXis Differential Optical Absorption Spectroscopy (MAX-DOAS) observation was carried out from November 2016 to February 2017 in Beijing, China to measure the O4 absorptions in UV and visible bands and further to illustrate its relationship with aerosol optical properties (AOPs) under different the weather types. According to 20 relative humidity, visibility and PM2.5, we classified the observation periods into clear, non-haze, haze, heavy-haze, fog and rainy five different weather conditions. There are obvious differences for measured AOPs under different weather conditions, especially scattering coefficient (σsca) and absorption coefficient (σabs). It was also found that both the O4 Differential Slant Column Densities (DSCDs) at UV and visible bands varied in the order of clear days > non-haze days > haze days > heavyhaze days > fog days. The correlation coefficients (R) between O4 DSCDs at 360.8 and 477.1 nm mainly varied in the order 25 of clear days > non-haze days > haze days > heavy-haze days. Based on the statistics of O4 DSCDs at elevation angle 1 with the corresponding linear regression between UV and visible bands of segmental periods, the relationships between O4 DSCDs and AOPs were established. It mainly should be clear or non-haze days when the correlation slope is greater than 1.0, correlation coefficient (R) greater than 0.9 and O4 DSCDs mainly greater than 2.5 × 10 molec cm. Meanwhile, σsca and σabs are less than 45 and 12 Mm , respectively. For haze or heavy-haze days, the correlation slope is less than 0.6, R less 30 than 0.8 and O4 DSCDs mainly less than 1.3× 10 molec cm, under which σsca and σabs are mainly located at 200-900 and 20-60 Mm. Additionally, the determination method was well validated based on another MAX-DOAS measurement at Gucheng from 19 to 27 November 2016. For more precise and accurate inversion of AOPs, more detailed look-up tables for Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2019-15 Manuscript under review for journal Atmos. Meas. Tech. Discussion started: 6 February 2019 c © Author(s) 2019. CC BY 4.0 License.


Introduction
Atmospheric aerosols influence the radiative budget by scattering and absorbing solar radiation directly.They also affect global climate change, cloud formation, regional air quality, and human health (Seinfeld and Pandis, 2006;Kim and Ramanathan, 2008;Karanasiou et al., 2012;Levy et al., 2013;Viana et al., 2014).It is important to obtain comprehensive knowledge on the spatial distributions, temporal variations in aerosol loading, and aerosol optical properties (AOPs).Different aerosol types have different optical properties.For example, black carbon (BC) aerosols are characterized by strong light absorption.Recent studies indicated that it can heat the air and contributes to global warming (Ramanathan et al., 2007;Galdos et al., 2013;Ramana et al., 2010;Fyfe et al., 2013;Allen et al., 2012).It can also change the atmospheric vertical temperature profile, causing variations in the planetary boundary layer (PBL) structure (Ding et al., 2016;Wilcox et al., 2016;Wang et al., 2018).However, dust aerosol and some secondary heterogeneousreaction aerosols, playing an important role during pollution episodes in China, are mainly based on scattering optical characteristics (Huang et al., 2014;Wang et al., 2018).
Measurements of AOPs, e.g., aerosol extinction coefficient (AEC), aerosol optical depth (AOD), single-scattering albedo (SSA), asymmetry factor, and Ångström exponent (AE), could provide more comprehensive information for a better understanding of the role of aerosols in atmospheric processes.AOD is an important parameter to evaluate the ability of aerosol particles to attenuate the solar radiation, which is defined as the integration of AEC from the surface to the top of atmosphere.The AE is the sum of aerosol scattering and absorption coefficients.Moreover, SSA is defined as the ratio of scattering efficiency to the total extinction, which is the dominant intensive parameter determining aerosol direct radiative forcing.The asymmetric factor is used to evaluate the aerosol forward scattering ability, while the Ångström exponent is a parameter to evaluate the aerosol particle size.Previous measurements of AOPs indicated that different aerosol types (such as biomass burning, urban-industrial, dust, and sea-salt aerosols) exhibit significant differences in optical properties.The differences of the optical properties of these kinds of aerosols are used to clarify the mechanisms of aerosol radiative forcing (Dubovik et al., 2001).For biomass burning aerosol, the Ångström exponent is mainly distributed between 1.1 and 2.1 at wavelength bands of 440-870 nm and SSA mainly ranging from ∼ 0.88 to 0.99 at 440 nm (Eck et al., 2003;Bergstrom et al., 2007;Weinzierl et al., 2017).The SSA of urban-industrial aerosol tends to be ∼ 0.95 in cleaner conditions and ∼ 0.85 in industrial conditions (Liousse et al., 1996;Remer and Kaufman, 1998;Garland et al., 2009;He et al., 2009;Shen et al., 2018).Dust exhibits a SSA of ∼ 0.92 to 0.93 in the blue spectral range at 440 nm, but ∼ 0.96-0.99 in wavelengths greater than 550 nm (Kaufman et al., 2001;Dubovik et al., 2001;Bergstrom et al., 2007;Weinzierl et al., 2017).The SSA in oceanic aerosol is mainly greater than 0.97 due to sea salt and water-soluble particles with high relative humidity (Tanré et al., 1999;Dubovik et al., 2001;Hess et al., 1998;Eck et al., 2005).
Multi-axis differential optical absorption spectroscopy (MAX-DOAS) remote sensing is an effective tool for atmospheric aerosol measurements based on O 4 molecular ultraviolet-visible light absorption (Platt and Stutz, 2008).O 4 is the collision complex of O 2 and its concentration is proportional to the square of the O 2 concentration.Due to the vertical O 4 profile being well known and nearly constant, it can be used as an indicator for the atmospheric distribution photon paths due to its nearly constant characteristics (Wagner et al., 2004;Frieß et al., 2006Frieß et al., , 2016)).The O 4 cross sections exhibit four main absorption bands in the UV-visible region around 360.8, 477.1, 577.1, and 630.8 nm (Thalman and Volkamer, 2013).By collecting the scattered sunlight spectra at zenith and different elevation angles closed to the horizon by MAX-DOAS, the O 4 absorptions in four bands (338-370, 425-490, 540-588, and 602-645 nm) can be estimated, and aerosol vertical profiles can be further derived (Honninger and Platt, 2002;Hytch et al., 2003;Hönninger et al., 2004;Wagner et al., 2004;Wittrock et al., 2004;Clémer et al., 2010).The sunlight at different wavelength bands has different abilities to traverse the atmosphere; thus the light path length at different wavelength bands is diverse, which can change the corresponding O 4 absorptions.Conversely, the correlation analysis between O 4 absorptions retrieved at UV range and visible range could also provide information about the impacts of aerosol scattering on photon paths (Lee et al., 2011).Besides the extinction coefficient profile and AOD, there are no previous detailed studies on the other AOPs' retrieval based on MAX-DOAS measurements.
In this paper, we try to establish a new method to determine several different aerosol optical properties from multiple-wavelength O 4 absorptions observed by MAX-DOAS measurement.The measurement of UV and visible O 4 absorptions was performed by a MAX-DOAS instrument in Beijing from November 2016 to February 2017.Combining the O 4 absorptions and measured AOPs, some empirical relationships between them can be found under different weather conditions, which are fundamental to determine the AOPs from MAX-DOAS-observed O 4 absorptions at different wavelength bands.Furthermore, another short measurement campaign was used to validate the feasibility and reliability of the new method to infer the AOPs under different weather conditions based on the O 4 absorptions.

The MAX-DOAS measurements
The MAX-DOAS instrument was installed on the roof of the Chinese Academy of Meteorological Sciences building (CAMS, 39.9475 • N, 116.3273 • E) for continuous measurements of O 4 absorptions from November 2016 to February 2017.The MAX-DOAS instrument consists of three major parts: a telescope unit, two spectrometers with temperature stabilized at 20 • , and a computer acting as the controlling and data acquisition unit.The viewing elevation angle of the telescope is controlled by a stepping motor.Scattered sunlight collected by the telescope is redirected by a prism reflector and a quartz fiber bundle to the spectrometers.Two spectrometers (Acton Spectrapro 300i Czerny-Turner optical spectrometer) are used to cover both the UV (300-460 nm) and visible (400-560 nm) wavelength ranges.The full width at half maximum (FWHM) spectral resolutions of these two spectrometers are all 0.6 nm, or 7.2 detector pixels.Moreover, the optical spectrometer is equipped with a chargecoupled device (CCD) detector camera (model DU 440-BU) with 2048 pixels.The field of view (FOV) of the instrument is estimated to be less than 0.3 • .
A full measurement scanning sequence consists of 11 elevation angles, i.e., 1, 2, 3, 4, 5, 6, 8, 10, 15, 30, and 90 • .The instrument azimuth angle is 138 • and the exposure time is fixed to 60 000 ms for each elevation angle measurement.A full measurement sequence takes about 11 min.Dark current and offset spectra were measured by blocking incoming light using a mechanical shutter and were subtracted from the measurement spectra before spectral analysis.The routine measurements were continuously repeated as long as the solar zenith angle (SZA) was lower than 80 • .

O 4 absorptions in UV and visible ranges
The O 4 absorptions were derived in the fitting windows of 339 to 387 nm in the UV range and 425 to 490 nm in the visible spectral interval.The measured spectra were analyzed using the QDOAS software developed by BIRA-IASB (http://uv-vis.aeronomie.be/software/QDOAS/,last access: 16 June 2019).The corresponding zenith spectrum was taken as a reference spectrum for off-zenith elevation angles during each measurement scan.The DOAS fitting generates the differential slant column density (DSCD) of O 4 between the measured spectra and reference spectrum.Details of DOAS fit settings are listed in Table 1.We used corrected I 0 (Aliwell et al., 2002).Figure 1 shows a typical DOAS retrieval for the O 4 absorptions at 360.8 and 477.1 nm.Afterwards, DOAS fit results with a root mean square (rms) larger than 5 × 10 −4 were filtered out, and about 99.07 % of all O 4 measurements remained for further discussion.

Ancillary data
Quality-assured level 2.0 sun photometer AODs, asymmetric factor, and Ångström exponent at the Beijing_CAMS AERONET site (http://aeronet.gsfc.nasa.gov/,last access: 16 June 2019) were employed, which is collocated with the MAX-DOAS instrument just 2 m away.A sun photometer (CE-318) collects direct sunlight about every 15 min during non-rainy daytime hours.These aerosol optical parameters at multiple wavelengths were normalized to 450 nm according to Wang et al. (2016).In addition, the scattering coefficients (σ sca ) were measured at three wavelengths (λ = 450, 520, and 700 nm) using an integrating nephelometer (Aurora 4000, Ecotech) at Peking University Urban Atmosphere Environment Monitoring Station (PKUERS, 39.9892 • N, 116.3069 • E).The absorption coefficients (σ abs ) were measured using a seven-wavelength Aethalometer (AE-31, Magee Scientific) at λ = 370, 470, 520, 660, 880, and 950 nm also located at PKUERS.Both σ sca and σ abs have the temporal sampling resolution of 5 min.In order to ensure the accuracy of the data, the corrections for σ sca and σ abs were used as referred to in Shen et al. (2018).The SSA was calculated by the measured σ sca at 450 nm and estimated σ abs at 450 nm using the following equation: The visibility and the relative humidity (RH) information were collected from the weather history data at Beijing international airport (http://www.wunderground.com/, last access: 16 June 2019) about 26 km from CAMS.In situ data of PM 2.5 concentrations were obtained from Guanyuan station (39.9425 • N, 116.3610 • E), belonging to the national environmental monitoring network (http://beijingair.sinaapp.com/, last access: 16 June 2019), which is about ∼ 2 km from the CAMS site.All these data are normalized to hourly averages for further discussion.

Wintertime aerosol optical properties
The time series of PM 2.5 concentrations, σ sca , σ abs , SSA, AOD, Ångström exponent, asymmetry factor, and the corresponding meteorological data, i.e., RH and visibility, from November 2016 to February 2017 are presented in Fig. 2. The typical meteorological conditions of high RH and low visibility coincided with significantly high PM 2.5 concentration and high AOD.As indicated in the gray areas of Fig.  value within a shorter period.PM 2.5 concentrations, σ sca , and σ abs increased to exceed the maximum values 465 µg cm −3 , 1331.151, and 123.402Mm −1 within 1-3 h (with increments up to 200 µg cm −3 , 600, and 100 Mm −1 ) during episode I, respectively.In episode II, the maximum values of PM 2.5 concentrations, σ sca , and σ abs are 585 µg cm −3 , 1473.523 and 153.431Mm −1 , respectively.In addition, SSA mainly remained over 0.85 during the entire wintertime; except it was observed to be less than 0.8 in late November 2016 and several days during January and February 2017.Generally, the high values of SSA were always accompanied by the peak of PM 2.5 concentrations, which suggests that the scattering properties of atmospheric aerosols were enhanced during the explosive increase stage of particle concentrations.Meanwhile, it is also associated with the decrease in Ångström exponent and the increase in asymmetry factors simultaneously.This is typically related to the particle size growth process (Guo et al., 2014;Yu et al., 2011Yu et al., , 2016)).
With the increase in pollution level indicated by AQI (except fog and rainy days), AOD increased dramatically from 0.311 under clear days to 1.338 on heavy-haze days.There are no obvious changes for σ sca and σ abs between clear days and light-haze days.Nevertheless, the σ sca increased sharply from light-haze days to heavy-haze days with the averaged value from 44.524 to 449.741 Mm −1 .The averaged value of σ abs is 8.257 Mm −1 in light-haze days and it increased as much as 5 times on heavy-haze days.Moreover, the averaged SSA was about 0.847 on light-haze days and similar to that on haze days, but it increased about 3.53 % from haze days to heavy-haze days with averaged values of 0.846 to 0.878.This suggests that the aerosol scattering and absorption abilities have changed but the ratio of scattering to extinction changed slightly as the processes of particle pollution became severe.No obvious variations in Ångström exponent were observed among clear days to heavy-haze days, but it decreased by about 2.83 % on fog days.In a previous study, the Ångström exponents were usually higher than 0.80 when AOD was greater than 0.60 in Beijing, which reveals the major contribution of small particles to higher aerosol loading (Che et al., 2015).However, our study demonstrates that small particles made a major contribution to the aerosols throughout the whole winter in Beijing.The obvious decrease in Ångström exponent on fog days is attributed to the increase in water vapor in particles.In addition, the averaged asymmetry factor was about 0.697 on fog days and 8.52 % higher than for other weather conditions.This indicates the increased forward scattering on fog days (Yoon and Kim, 2006).DSCDs at 477.1 nm are obviously higher than at 360.8 nm on clear and light-haze days, and slightly larger than that at 360.8 nm on haze and heavy-haze days, which can be explained by the fact that the observable light path length in the visible range is longer than in the UV range.Even in UV bands, the observed O 4 DSCDs at 353 nm were reported to be lower than those at 380 nm for most of the elevations under haze conditions during winter in Beijing (Lee et al., 2011).This phenomenon revealed that O 4 absorptions in the short wavelength range were more significantly affected by light diffusion under hazy conditions.However, we found there are no obvious differences between O 4 DSCDs at 360.8 and 477.1 nm on fog days, during which the high contents of water vapor decreased the visibility and the atmospheric absorption paths from the UV to visible ranges.We further analyzed the relationship of O 4 absorptions between the UV and visible bands.As shown in the right column of Fig. 3, the correlation coefficient (R 2 ) of O 4 DSCDs between 360.8 and 477.1 nm varied in the order of clear days > light-haze days > haze days > heavy-haze days.Strong correlation between UV and visible O 4 absorptions (R 2 > 0.9) was achieved for clear and light-haze days.Under haze and heavy-haze conditions, R 2 was 0.81 and 0.74, respectively, which is much lower than that on clear and lighthaze days.That is because the increase in light-absorbing and light-scattering aerosols can more obviously result in reduced light path lengths in shorter wavelength bands than longer wavelength bands during haze and heavy-haze days.

UV and visible O 4 absorptions under different weather conditions
The changes in AOPs, especially aerosol scattering and absorption properties, are mainly manifested in the variations in O 4 absorptions at different wavelength bands.The correlation information of O 4 DSCDs at different bands will also    be affected by the variation in AOPs.For more details, i.e., 21 February 2017, was chosen to exhibit the influence of AOP changes on O 4 DSCDs in Fig. 4. Comparing Fig. 4a to b, it can be found that σ sca and σ abs have similar variation trends; a slight turning and an abrupt decrease occurred at ∼ 09:05 and ∼ 12:00 (UTC+8 for all times throughout) (especially for σ sca ), respectively, while the time-indicated O 4 DSCDs seem to be three segments with a higher correlation coefficient divided by the breaks at 10:00 and 12:00.
In order to explore the relationship between the O 4 DSCDs at different wavelength bands and the variations in σ sca and σ abs , we defined the change speed (v sca ), acceleration (a sca ), and change rate of acceleration (a sca ) of σ sca (Fig. 4b, e, and  f) as the following three equations.
Accordingly, the relevant time series of v sca , a sca , and a sca are displayed in Fig. 4d to f.In this case, we can find two time break points, defined as t 1 and t 2 (t 1 = 10:00 and t 2 = 12:00), at which σ sca and σ abs have significant variations based on the calculated v sca , a sca , and a sca .In addition, we found the indicator of a sca can describe the specific moment of change (increasing or decreasing) in σ sca more clearly than v sca and a sca in this case.|a scat 1 | and |a scat 2 | are all higher than 20 000.Consequently, the O 4 DSCDs at 360.8 and 477.1 nm can be divided into three segments for the periods of 08:00-10:00, 10:00-12:00, and 12:00-16:00 and the correlation between UV and visible O 4 DSCDs was further analyzed individually.As shown in Fig. 4c, the R 2 during 08:00-10:00 and 10:00-12:00 is obviously larger than that of all day in Fig. 4a; however, it is smaller for the segment of 12:00-16:00.Moreover, there were huge divergences among the correlation slopes among these three segments due to the change in aerosol scattering and absorption properties.Therefore, it can be concluded that the diurnal variations in O 4 DSCDs provide the information of the light path length impacted by aerosol loading.Further, the varied relationship between O 4 DSCDs at UV and visible ranges implies change in the aerosol scattering and absorption properties.
Using the method discussed above, we have defined the time break points with aerosol property changes and further classified the observations into several segmental periods with the criteria of |v sca | > 1000, |a sca | > 10 000, or |a sca | > 20 000.The summary of time break points and corresponding change speed (v sca ), acceleration (a sca ), and change rate of acceleration (a sca ) of σ sca are listed in Table S1 in the Supplement.

Implications of O 4 absorptions for aerosol optical properties
In order to derive the aerosol optical properties from multiple-wavelength O 4 absorptions, the complete 4 months of observational O 4 and AOPs data were used for discussion under different weather conditions.Hourly data of O 4 DSCDs were divided into five weather conditions and made the linear regression between UV and visible O 4 DSCDs.In total, there were about 218 segments (776 h in 97 d), including 67, 31, 61, 44, and 15 segments for clear, light-haze, haze, heavy-haze, and fog days, respectively.Figure 5 illustrates the statistics of O 4 DSCDs in UV and visible bands, and the slope and R 2 of correlation analysis between them, as well as the O 4 DSCD ratio of UV to visible bands for different weather conditions.In general, the O 4 DSCDs in UV wavelengths mainly ranged in 3.00-4.00× 10 43 , 2.50-3.50× 10 43 , 0.50-1.10× 10 43 , 0.25-0.80× 10 43 , and 0.20-0.40× 10 43 molec.cm −2 in clear, light-haze, haze, heavyhaze, and fog days, respectively.And the O 4 DSCDs in visible wavelengths are mainly distributed between 4.00 and 6.50 × 10 43 , 3.00 and 5.50 × 10 43 , 0.50 and 1.30 × 10 43 , 0.25 and 0.60 × 10 43 , and 0.25 and 0.60 × 10 43 molec.cm −2 under the above five different weather conditions, which are higher than those in UV wavelengths, especially for clear and light-haze days.Moreover, the corresponding ratio of visible to UV O 4 DSCDs is 1.45-1.70,1.45-1.65,1.00-1.65,0.85-1.35,and 0.80-1.35under these five weather conditions, respectively.The linear regression results show that the correlation slopes between UV and visible O 4 DSCDs are greater than 1.00 (mainly greater than 1.40) and the correlation R 2 values are greater than 0.93 mostly on clear days.Under light-haze conditions, the correlation slopes are greater than 1.00 (mainly greater than 1.20) and the correlation R 2 values are mainly greater than 0.90.The correlation slopes are mainly less than 0.60 and the correlation R 2 values have a wider range (maximum value < 0.80 and occasional fitting failure) on haze days.On heavy-haze days, the correlation slopes are less than 0.60-0.80 and the correlation R 2 values are 0.50-0.80mostly (some fitting failure cases appeared).On fog days, the correlation slopes are floated around 1.00 and the correlation R 2 values are mainly 0.75-0.85.
Meanwhile, the statistical characteristics of AOPs under different weather conditions are shown in Fig. 6.Similar to the results in Table 2, σ sca shows the increasing trend and was mainly distributed between 21. 83 and 47.01, 28.33 and 57.25, 134.72 and 349.47, 228.42 and 649.11, and 450.00 and 1004.88Mm −1 on clear, light-haze, haze, heavy-haze and fog days, respectively.The σ abs values were mainly distributed between 2. 61 and 8.26, 3.99 and 11.89, 22.25 and 57.31, 34.84 and 72.22, and 70.01 and 115.26Mm −1 under the five different weather conditions above.The AODs were mainly distributed between 0.11 and 0.35, 0.12 and 0.37, 0.76 and 1.70, 1.37 and 2.38, and 0.69 and 1.38 under these five weather conditions.The Ångström exponents were more disperse for clear and light-haze days than the other weather conditions.Except the fog days, the asymmetry factor in other weather conditions is not much different.
Combining the statistical information on O 4 absorptions and AOPs, we could conclude some empirical relationships as follows.
Under the condition that the correlation slopes between UV and visible O 4 DSCDs are greater than 1.0 and the correlation R 2 values are greater than 0.9, and simultaneously, the UV and visible O 4 DSCDs are mainly greater than 2.5 × 10 43 molec.cm −2 and 3.0 × 10 43 molec.cm −2 , we know the weather should mainly be clear or light-haze days.It can be suspected that σ sca and σ abs are less than 45 and 12 Mm −1 , and AODs are less than 0.4.
If the correlation slope is around 1.0 and with a correlation R 2 of 0.75-0.85,we know the weather should mainly be fog days.σ sca and σ abs are located at 450-1200 and 60-90 Mm −1 , while AODs are greater than 0.7.Therefore, the obtained empirical relationships represent the potential ability to determine the basic aerosol loading conditions from the MAX-DOAS observed O 4 absorptions.

Discussion
To investigate the feasibility and reliability, another short MAX-DOAS measurement campaign operating in Gucheng, Hebei Province (39.1382 • N, 115.7163 • E) from 19 to 27 November 2016 was used for the application of the new method to determine AOPs from O 4 absorptions.The MAX-DOAS instrument is the same as that installed in CAMS.Due to the absence of a sun photometer instrument, AODs at 450 nm were obtained by profiling the aerosol extinction coefficient based on MAX-DOAS measurements by utilizing the optimal estimation method (Frieß et al., 2006(Frieß et al., , 2016;;Xing et al., 2017).In addition, σ sca and σ abs were acquired using the same co-located integrating nephelometer (Aurora 1000, Ecotech) and seven-wavelength Aethalometer (AE-31, Magee Scientific), respectively.Moreover, the temporal res-olution of σ sca and σ abs measurements is 1 min and 1 respectively.
Figure 7a and b show the diurnal variations and segmental correlation of O 4 DSCDs in the UV and visible bands during this campaign.According to the empirical relationships discussed in Sect.3.3, it can be inferred that the period segment during 09:00-11:00 on 25 November should be haze or heavy-haze conditions because the UV and visible O 4 DSCDs are all less than 0.5 × 10 43 molec.cm −2 , and simultaneously the correlation slope and R 2 between them are 0.42 and 0.59, which is in line with the determination conditions that UV and visible O 4 DSCDs are mainly less than 1.1 × 10 43 and 1.3 × 10 43 molec.cm −2 and simultaneously the correlation slope and R 2 between them are mainly less than 0.6 and 0.8.Similarly, the periods of 09:00-12:00 on 21 November, 10:50-16:00 on 22 November, 10:00-15:00 on 23 November, 08:00-15:00 on 26 November and 11:00-15:00 on 27 November are mainly clear or light-haze weather.And 09:00-10:00 on 19 November, 09:00-12:00 on 20 November, and 09:00-10:50 on 22 November can be mainly regarded as haze or heavy-haze weather types.
Furthermore, the time series of in situ σ sca , σ abs , and MAX-DOAS-retrieved AOD are shown in Fig. 7c and d.According to the empirical relationships summarized above, the σ sca , σ abs , and AOD are mainly located at 200-900, 20-60 Mm −1 , and 0.9-2.5 under the haze segment of 09:00-11:00 of 25 November.Simultaneously, the in situ measured σ sca , σ abs , and MAX-DOAS-retrieved AOD during the above periods range in 588.30-730.77, 58.19-67.63Mm −1 , and 1.39-2.22.The inferred results are in good agreement with the measured results.This indicates that the concluded empirical relationships can be used as the criterion to accurately determine the ranges of aerosol optical parameters of σ sca , σ abs , and AOD.Nevertheless, we found two segments with correlation slopes > 1.0 and R 2 < 0.9 during 12:00-15:00 on 19 November and 11:00-15:00 on 25 November, which is not included in cases of the empirical relationships.This suggests that more refined and quantitative relationships between aerosol optical parameters and O 4 absorptions need to be further achieved with the increases in the measured data, which can be established as a look-up table to retrieve the aerosol optical properties in the future.
Moreover, in order to illustrate the variations in the O 4 absorptions due to the change in aerosol loadings, we used a radiative transfer model of SCIATRAN to simulate O 4 DSCDs in the UV and visible bands under conditions with different aerosol optical properties, but under a fixed given series of observation geometry with elevation angle, solar zenith angle, and relative azimuth angle.In total, 11 different aerosol scenarios were simulated, and the linear-regression analysis for the simulated UV and visible O 4 DSCDs under different aerosol conditions were performed.The forward radiative transfer model (RTM) simulation results presented in the Supplement (Table S2 and Fig. S1) demonstrate that the O 4 absorptions information, including the value of UV and Vis- ible O 4 DSCDs, the corresponding linear-regression slope and R 2 between them, could greatly reflect the variation in aerosol optical properties.More details can be found in the Supplement.

Summary and conclusions
Ground-based MAX-DOAS measurements for O 4 DSCDs in the UV and visible wavelength bands were carried out in Beijing from November 2016 to February 2017.Combined with the measured σ sca and σ abs and AOD, we have summarized the characteristics of O 4 absorptions and parameters of AOPs under different weather conditions during autumn-winter seasons.It was found that the averaged AOD increased from 0.311 on clear days to 1.338 on heavy-haze days.The averaged σ sca changed dramatically from 44.524 Mm −1 on lighthaze days to 449.741 Mm −1 on heavy-haze days.Moreover, the averaged σ abs also obviously increased from 8.257 Mm −1 on light-haze days to 53.257 Mm −1 for heavy-haze conditions.Both the measured UV and visible O 4 DSCDs varied in the order of clear days > light-haze days > haze days > heavyhaze days > fog days.The corresponding correlation information (slope and R 2 ) between O 4 DSCDs at UV and visible wavelength bands also changed synchronously when σ sca and σ abs have varied.
Considering the simultaneous variation in O 4 absorptions and AOPs, the segmental period correlation analysis between UV and visible O 4 DSCDs was performed.Afterwards, the empirical relationships between O 4 absorptions and AOPs could be concluded for different aerosol load-ings.Clear and light-haze days occur under linear-regression slopes greater than 1.0 and R 2 mainly greater than 0.9, and simultaneously UV and visible O 4 DSCDs mainly greater than 2.5 × 10 43 molec.cm −2 .σ sca , σ abs , and AOD are mainly less than 45, 12 Mm −1 , and 0.4 under these conditions, respectively.When the correlation slopes, R 2 , and O 4 DSCDs are less than 0.6, 0.8, and 1.3 × 10 43 molec.cm −2 , respectively, it should mainly be haze or heavy-haze days.Under this condition, the σ sca , σ abs , and AOD can be inferred to be mainly located at 200-900, 20-60 Mm −1 , and 0.9-2.5, respectively.In addition, it should be fog days if the correlation slopes are around 1.0 and R 2 is 0.75-0.85.Another MAX-DOAS measurement campaign carried out at Gucheng from 19 to 27 November 2016 was used to validate the proposed new method, which could determine the AOPs from the observed O 4 absorptions well.
In this paper, we present a new method to directly deduce the parameters of aerosol optical properties from the observed UV and visible O 4 absorptions, which expands the usage of the MAX-DOAS technique to quickly semi-quantify the aerosol scattering and absorption properties.With the improvement of the look-up table, more precise and accurate inversion of aerosol optical properties can be achieved.Since only the O 4 DSCDs at an elevation angle of 1 • were employed to obtain the aerosol scattering and absorption at the surface, it can be expected that vertical spatially resolved of aerosol scattering and absorption can be retrieved by using O 4 DSCDs at different elevation angles in a future study.

Figure 1 .
Figure 1.Typical DOAS spectral fittings for O 4 absorptions in the (a) UV and (b) visible bands.Black lines represent the absorption signal and the red lines represent the sum of the absorption signal and the fit residual.

Figure 2 .
Figure 2. Time series of PM 2.5 and AOPs (σ sca , σ abs , SSA, AOD, Ångström exponent, and asymmetry factor) and meteorological parameters (relative humidity and visibility) during the observation in Beijing from November 2016 to February 2017.

Figure 3
Figure3shows the examples of diurnal pattern and the corresponding correlation of UV and visible O 4 DSCDs (elevation angle = 1 • ) at 360.8 and 477.1 nm under five different weather conditions, except for the rainy days.In view of the absolute strength of O 4 absorption, both the O 4 DSCDs in the UV and visible bands varied in the order of clear days > lighthaze days > haze days > heavy-haze days > fog days.This shows the dependency of O 4 absorption on the scattering sunlight path impacted by the aerosol loading.Moreover, O 4 DSCDs at 477.1 nm are obviously higher than at 360.8 nm on clear and light-haze days, and slightly larger than that at 360.8 nm on haze and heavy-haze days, which can be explained by the fact that the observable light path length in the visible range is longer than in the UV range.Even in UV bands, the observed O 4 DSCDs at 353 nm were reported to be lower than those at 380 nm for most of the elevations under haze conditions during winter in Beijing(Lee et al., 2011).This phenomenon revealed that O 4 absorptions in the short wavelength range were more significantly affected by light diffusion under hazy conditions.However, we found there are no obvious differences between O 4 DSCDs at 360.8 and 477.1 nm on fog days, during which the high contents of water vapor decreased the visibility and the atmospheric absorption paths from the UV to visible ranges.We further analyzed the relationship of O 4 absorptions between the UV and visible bands.As shown in the right column of Fig.3, the correlation coefficient (R 2 ) of O 4 DSCDs between 360.8 and 477.1 nm varied in the order of clear days > light-haze days > haze days > heavy-haze days.Strong correlation between UV and visible O 4 absorptions (R 2 > 0.9) was achieved for clear and light-haze days.Under haze and heavy-haze conditions, R 2 was 0.81 and 0.74, respectively, which is much lower than that on clear and lighthaze days.That is because the increase in light-absorbing and light-scattering aerosols can more obviously result in reduced light path lengths in shorter wavelength bands than longer wavelength bands during haze and heavy-haze days.The changes in AOPs, especially aerosol scattering and absorption properties, are mainly manifested in the variations in O 4 absorptions at different wavelength bands.The correlation information of O 4 DSCDs at different bands will also

Figure 3 .
Figure 3.Diurnal variation and correlation analysis of O 4 DSCDs at 360.8 and 477.1 nm under different weather conditions: (a) and (b) clear day, (c) and (d) light-haze day, (e) and (f) haze day, (g) and (h) heavy-haze day, (i) and (j) fog day.The color bar represents the time sequence.

Figure 4 .
Figure 4.An example day on 21 February 2017: (a) the correlations between O 4 DSCDs at 360.8 and 477.1 nm.The color bar represents time sequence.Panel (b) shows the time series of aerosol scattering and absorption coefficients.The correlation information between O 4 DSCDs at 360.8 and 477.1 nm at 08:00-10:00, 10:00-12:00, and 12:00-16:00 on 21 February 2017 is shown in (c).Panels (d) to (f) show the time series of v sca , a sca , and a sca of scattering coefficients, respectively.

Figure 5 .
Figure 5. Box plots of statistics for the O 4 DSCD under different weather conditions: (a) in the UV band, (b) in the visible band, (c) the ratio of visible / UV O 4 DSCD, (d) correlation slope, and (e) R 2 between the UV and visible ranges (visible O 4 DSCD = slope • UV O 4 DSCD + intercept).

Table 2 .
Statistics of AQI and several aerosol optical properties under different weather conditions.