Strong light scattering of highly oxygenated organic aerosols impacts significantly on visibility degradation

Secondary organic aerosols (SOAs) account for a large fraction of atmospheric aerosol mass and play significant roles in visibility impairment by scattering solar radiation. However, comprehensive evaluations of SOA scattering abilities under ambient relative humidity (RH) conditions on the basis of field measurements are still lacking due to the difficulty of simultaneously direct quantifications of SOA scattering efficiency in dry state and SOA water uptake abilities. In this study, field measurements of aerosol chemical and physical properties were conducted in winter in Guangzhou (lasting about 3 months) using a humidified nephelometer system and aerosol chemical speciation monitor. A modified multilinear regression model was proposed to retrieve dry-state mass scattering efficiencies (MSEs, defined as scattering coefficient per unit aerosol mass) of aerosol components. The more oxidized oxygenated organic aerosol (MOOA) with an O/C ratio of 1.17 was identified as the most efficient light scattering aerosol component. On average, 34 % mass contribution of MOOA to total submicron organic aerosol mass contributed 51 % of dry-state organic aerosol scattering. The overall organic aerosol hygroscopicity parameter κOA was quantified directly through hygroscopicity closure, and hygroscopicity parameters of SOA components were further retrieved using a multilinear regression model by assuming hydrophobic properties of primary organic aerosols. The highest water uptake ability of MOOA among organic aerosol factors was revealed with κMOOA reaching 0.23, thus further enhancing the fractional contribution of MOOA in ambient organic aerosol scattering. In particular, the scattering abilities of MOOA were found to be even higher than those of ammonium nitrate under RH of < 70%, which was identified as the most efficient inorganic scattering aerosol component, demonstrating that MOOA had the strongest scattering abilities in ambient air (average RH of 57 %) during winter in Guangzhou. During the observation period, secondary aerosols contributed dominantly to visibility degradation (∼ 70%), with substantial contributions from MOOA (16 % on average), demonstrating significant impacts of MOOA on visibility degradation. The findings of this study demonstrate that more attention needs to be paid to SOA property changes in future visibility improvement investigations. Also, more comprehensive studies on MOOA physical properties and chemical formation are needed to better parameterPublished by Copernicus Publications on behalf of the European Geosciences Union. 7714 L. Liu et al.: Strong light scattering of highly oxygenated organic aerosols ize its radiative effects in models and implement targeted control strategies on MOOA precursors for visibility improvement.


Introduction
Atmospheric aerosols directly scatter and absorb solar radiation and thus have significant radiative effects on both the Earth-atmosphere radiative budget and the atmospheric environment.Aerosols represent the dominant contributor to atmospheric visibility impairment in polluted regions (Liu et al., 2017a).With rapid industrialization and urbanization, China has experienced severe haze pollution in the last decade, and frequent low visibility events have aroused public attention and concern, especially since 2013.In recent years, the Chinese government has implemented stringent control policies called "Blue Sky Actions" to lower aerosol mass concentration and improve atmospheric visibility.However, Xu et al. (2020) revealed that the visibility improvement in southern China was less than expected; it was especially poor in the Pearl River Delta region, due to the nonlinear responses of the visibility improvement to PM 2.5 (particulate matter with aerodynamic diameter less than 2.5 µm) mass concentration reduction.Several recent literature reports also proved that visibility was less improved than PM 2.5 mass concentrations.Results of Liu et al. (2020) demonstrated that increased aerosol extinction efficiency associated with nitrate was responsible for the less improved visibility in eastern China.Hu et al. (2021) raised the challenge of visibility improvement due to increased nitrate contribution in Beijing area.However, results of Xu et al. (2020) demonstrate that this situation was likely associated with both increased aerosol scattering efficiency and aerosol hygroscopicity and particularly pointed out that other than changes of inorganic aerosol components, special attention should be paid to scattering efficiency and hygroscopicity changes of secondary organic aerosol (SOA).
Organic aerosols including primary and secondary organic aerosols (POAs and SOAs) represent a large and sometimes even dominant fraction of submicron aerosol mass (Jimenez et al., 2009).In particular, SOA was found to contribute dominantly to total organic aerosol mass under polluted hazy conditions (Huang et al., 2014;Kuang et al., 2020b).Wang et al. (2019) reported increased contributions of both secondary organic and inorganic aerosol mass across China due to clean air actions, and the nonlinear responses of secondary aerosol mass concentration to emission reductions were further confirmed during the COVID lockdown period as reported by Huang et al. (2020).Xu et al. (2019) also reported substantial changes of SOA properties such as enhanced oxidation state.However, most previous studies only paid attention to influences of nitrate increase on visibility degradation, whereas synergistic effects of SOA on visibility have never been the focus due to the complexity of SOA hygroscopicity and scattering efficiency.Organic aerosol evolves in the atmosphere including their sizes and chemical structures and thus also their optical properties and hygroscopicity (Jimenez et al., 2009), leading to the difficulty of quantifying contributions of organic aerosol in visibility degradation.Li et al. (2022) reported that nitrate and SOA dominated particle extinction in dry state in Beijing due to clean air actions; however, evaluations in ambient air were lacking, further stressing the importance of comprehensive evaluations of SOA scattering abilities under ambient relative humidity (RH) conditions to better elucidate the roles of SOA in visibility degradation and long-term visibility changes.
In this study, we comprehensively quantified the dry-state mass scattering efficiencies (MSEs) of both primary and secondary organic aerosol components and organic aerosol hygroscopicity, thus systematically evaluating contributions of SOA factors to aerosol scattering and visibility degradation in ambient air.

Campaign information
Aerosol physical and chemical properties were simultaneously measured during winter from 13 December 2020 to 25 February 2021 in the Guangzhou urban area.Instruments were housed in an air-conditioned container in Haizhu wetland park (Sect.S1 in the Supplement for site information).A PM 2.5 inlet (BGI, SCC 2.354) was used for aerosol sampling, and a sample flow of 8-9 L min −1 was maintained during the observation period, thus generally satisfying the flow requirement (8 L min −1 ) of 2.5 µm cutting diameter.A Nafion dryer was used to lower sample RH to less than 40 %.A humidified nephelometer system with a total flow of about 5 L min −1 and a quadrupole aerosol chemical speciation monitor (Q-ACSM) with a flow of 3 L min −1 were placed downstream of the dryer to measure aerosol scattering abilities under controlled RH conditions and aerosol chemical compositions.An AE33 aethalometer (Drinovec et al., 2015) set up with a flow rate of 5 L min −1 was separately operated under another inlet (PM 2.5 , BGI SCC 1.829) to measure aerosol absorptions and thus indirectly black carbon (BC) mass concentration.Measurements of meteorological parameters such as temperature, wind speed and direction, and RH were made by an automatic weather station.

Humidified nephelometer measurements
The humidified nephelometer system is a laboratoryassembled one, including two Aurora 3000 nephelometers, one measuring aerosol scattering abilities (aerosol scattering and back scattering coefficients at 450, 525 and 635 nm) under low RH conditions (mostly less than 30 %, dry nephelometer) and another one measuring aerosol scattering abilities under controlled RH conditions (wet nephelometer).The humidified nephelometer system can operate either in fixed-RH mode or in scanning-RH mode; details about techniques of fixed-RH mode and scanning-RH mode were introduced in detail in several previous studies (Kuang et al., 2017(Kuang et al., , 2020c)).In this study, the humidified nephelometer system was operated in scanning-RH mode before 26 January 2021 and in fixed-RH mode (80 % RH) from 26 January to 9 February.The RH ranges of scanning-RH mode were 75 %-90 % from 11 December 2020 to 5 January 2021 and 60 %-90 % from 13 to 26 January 2021.The humidified nephelometer system provides direct measurements of the aerosol light scattering enhancement factor f (RH, λ) = σ sp (RH,λ) σ sp (dry λ) , where σ sp (RH, λ) is the aerosol scattering coefficient of light wavelength λ at condition of RH (Titos et al., 2016;C. Zhao et al., 2019), and f (RH, 525) is referred directly to as f (RH) was usually used to derive the optical hygroscopicity parameter κ sca through f (RH) = 1 + κ sca × RH 100−RH (Brock et al., 2016;Kuang et al., 2017Kuang et al., , 2018Kuang et al., , 2020cKuang et al., , 2021b;;Xu et al., 2020).The nephelometer measurements are associated with truncation error and nonideality of the light source.The dry-state aerosol scattering coefficients were corrected using the empirical formula provided by Qiu et al. (2021).RH 0 in the dry nephelometer was in the range of 6 %-49 % with an average of 22 %; thus dry-state aerosol scattering coefficient at 525 nm (σ sp,525 ) was further corrected using measured aerosol optical hygroscopicity through σ sp,525 = σ sp,525,measured / 1 + κ sca × RH 0 100−RH 0 .

Q-ACSM measurements and positive matrix factorization (PMF) analysis
The Q-ACSM was deployed to routinely characterize and measure the mass concentrations of nonrefractory submicron aerosol components at a time resolution of 15 min, including organics, sulfate, nitrate, ammonium and chloride; for details about the Q-ACSM setup, please refer to Ng et al. (2011).
The mass concentrations and mass spectra were processed using ACSM data analysis software (ACSM Local 1.5.10.0Released 6 July 2015) based on Igor Pro (version 6.37).The detailed procedures were described in Ng et al. (2011) and Sun et al. (2013).A composition-dependent CE value consistent with Middlebrook et al. (2012) and Sun et al. (2013) was chosen considering that aerosol samples were dried before entering the ACSM instrument.).The mass spectra of these factors (Fig. S9 in the Supplement) and more details about the factor analysis can be found in Sect.S2.

Organic aerosol hygroscopicity derivation
Organic aerosol was usually treated as nearly hydrophobic in many previous studies when considering environmental effects of organic aerosol (Cheng et al., 2016); however it is quantified in this study based on the most recently developed organic aerosol hygroscopicity quantification method by Kuang et al. (2020a).On the basis of field measurements, the organic aerosol hygroscopicity parameter κ (κ OA ) can only be estimated through closure between directly measured overall aerosol hygroscopicity and aerosol chemical compositions using the volume mixing rule (Kuang et al., 2020d).Kuang et al. (2020c) developed an optical method to calculate κ OA based on the combination of the humidified nephelometer system measurements and bulk PM 1 aerosol chemical-composition measurements, and the application of this method was further manifested and discussed in Kuang et al. (2021b) and thus used in this study to estimate κ OA .The humidified nephelometer system provides direct measurements of the optical hygroscopicity parameter κ sca and aerosol scattering Ångström exponent, which can be used together to derive a κ value referred to as κ f (RH) (Kuang et al., 2017), which can be treated as the overall aerosol hygroscopicity parameter in the hygroscopicity closure (Kuang et al., 2021a).In the closure, ions were paired using the scheme proposed by Gysel et al. (2007) as listed in Table S1 in the Supplement.The same as in Kuang et al. (2021b)  can be estimated using the following formula by assuming volume additivity and zero κ of BC: where ε represents the volume fraction whose calculation needs total aerosol volume concentrations, and the subscript represents name of an aerosol component.The total PM 1 aerosol volume concentration was calculated using measurements of the dry nephelometer following the machine learning method proposed by Kuang et al. (2018).Organic aerosol density varies over a wide range, and previous studies demonstrate that it is tightly associated with aerosol oxidation state (Kuwata et al., 2012) and that a higher O/C ratio usually corresponds to higher aerosol density.Following Wu et al. (2016), densities of primary organic aerosol (POA) components including HOA and COA were assumed to be 1 g cm −3 , and density of MOOA was assumed to be 1.4 g cm −3 due to its highly oxygenated feature with O/C of 1.17; however, the density of LOOA was assumed to be 1.2 g cm −3 due to its moderate O/C value of 0.7.In addition, the difference between predicted volume concentration from nephelometer measurements and the total aerosol volume concentrations summed from known aerosol components was attributed to aerosol components unidentified by the mass spectrometer; thus its volume fraction and hygroscopicity were named ε X and κ X in Eq. ( 1), where κ X was assumed to be 0.05 since the unidentified parts are usually metals and dust with low aerosol hygroscopicity.For more comprehensive discussions on κ OA estimations as well as κ X treatment, please refer to Kuang et al. (2021b).
Using the ZSR (Zdanovskii, Stokes and Robinson) mixing rule, the κ OA derived at RH of 80 % can be further expressed as where ε represents volume fractions of primary and secondary organic aerosol components in total organic aerosols.Assuming κ values of HOA and COA as zero, Eq. ( 2) can be simplified as κ OA = κ MOOA ×ε MOOA +κ LOOA ×ε LOOA .Considering the noisy characteristics of derived κ OA as shown in Fig. 4a, this simplified formula was not directly used to fit all derived κ OA values.Instead, average diurnal variations of derived κ OA were firstly acquired and then fitted using κ OA = κ MOOA × ε MOOA + κ LOOA × ε LOOA with average diurnal volume fractions of MOOA and LOOA (ε MOOA and ε LOOA ) in total organic aerosol as inputs, which yields average κ MOOA and κ LOOA .
3 Results and discussions

Overview of the pollution characteristics during winter in Guangzhou
The time series of corrected aerosol scattering coefficients in dry state (σ sp,525 ), nonrefractory PM 1 (NR PM 1 ) mass concentrations and resolved organic aerosol factors, as well as meteorological parameters, are shown in Fig. 1.The scattering coefficient at 525 nm of aerosols in dry state (σ sp,525 ) varied over a wide range of 8 to 688 Mm −1 , with an average of 118 Mm −1 .The average NR PM 1 is 20 µg m −3 , with the highest NR PM 1 mass concentrations reaching beyond 160 µg m −3 .This suggests relatively clean conditions compared to aerosol pollution in other polluted regions in China (Zhou et al., 2020); however, severe pollution episodes occurred occasionally.Three haze pollution episodes characterized by relatively high aerosol mass loading and scattering coefficients were observed before February (gray shaded areas in Fig. 1).The evolution and formation of these three episodes were tightly associated with the stagnant meteorological conditions, as indicated by the low wind speed (< 2 m s −1 ) and increasing RH during the last two severe pollution episodes.As shown in Fig. 1c, organic aerosol and nitrate contributed dominantly to the increase of aerosol mass, while sulfate remained almost flat or increased slightly during these pollution episodes.For example, for the episode in which the worst pollution conditions occurred, the daily average NR PM 1 increased from 19 to 143 µg m −3 from 12 to 16 February, with the organic aerosol increasing from 9.3 to 69.8 µg m −3 and nitrate increasing from 5.5 to 44.2 µg m −3 ; however, sulfate only increased from 1.4 to 8.5 µg m −3 .This phenomenon is quite different from the results reported by Guo et al. (2020) and W. Chen et al. ( 2021) that organic aerosol dominated the aerosol mass increase with a significant increase of both sulfate and nitrate in pollution episodes in autumn of Guangzhou urban area; however, it is generally consistent with the increasing characteristics reported by C. Chen et al. (2021).These observations suggest that the aerosol pollution differs much among seasons and years due to the highly variable characteristics of meteorological conditions.As for the organic aerosol mass increase, the time series of resolved organic aerosol factors are also shown in Fig. 1d and e.For the three observed pollution episodes, both increases of primary and secondary organic aerosol (represented by summation of MOOA and LOOA) were observed, with LOOA contributing dominantly to SOA increase and HOA contributing dominantly to POA increase.However, the accumulation of POA (summation of HOA and COA) contributed almost twice as much as the increase of SOA, suggesting that primary emissions, especially vehicle emissions, played significant roles in aerosol mass increase during pollution episodes of winter in Guangzhou.

Strong scattering ability of MOOA in dry state
Traditional multiple linear regression models were usually applied to determine MSEs of different aerosol components using simultaneously measured aerosol scattering coefficients and mass concentrations of aerosol components (Hand and Malm, 2007;Han et al., 2015;Chan et al., 1999).However, the traditional model failed in this study due to co-variations of input variables and impactor inconsistencies (PM 2.5 versus PM 1 ) between aerosol chemical compositions and aerosol scattering measurements.Details about this failure are discussed in Sect.S4.A new methodology was proposed to lower correlations between variables and reduce the impacts of measurement inconsistency of aerosol populations between the nephelometer and the mass spec-  creases of LOOA cases were almost always accompanied by the spontaneously quick HOA increase because most LOOA rapid increase cases happened during pollution episodes and start near dusk when accumulation of vehicle emissions and nitrate formation occurred.Thus, the average increase of LOOA was even smaller than that of AN and HOA.Cases of slight but significant MOOA increase were also identified, and average results are also shown in Fig. 2d, showing that MOOA increases were usually accompanied by significant nitrate formation.These results demonstrate that MSEs of aerosol components cannot be quantified directly from responses of aerosol scattering to aerosol emission or formation processes.However, for these cases, mass increases of aerosol components and corresponding changes in aerosol scattering matter most, and impacts of unidentified aerosol components are reduced substantially through differentially considering that the average time change for these identified cases is only 4 h.In addition, as listed in Table S3, the correlations between changes of most variables for all identified cases are much smaller than their time series correlations shown in Table S1.Thus, the modified multiple linear regression model σ sp,525 = HOA × MSE HOA + COA × MSE COA + LOOA×MSE LOOA + MOOA×MSE MOOA + AS × MSE AS + AN × MSE AN + BC × MSE BC was applied to retrieve MSEs of aerosol components.The derived MSEs for HOA, COA, LOOA, MOOA, AN, AS and BC were 2.1, 3.9, 3.4, 9.9, 7.1, 5.5 and 3.3 m 2 g −1 , respectively.The fitted σ sp,525 correlated highly (R 2 = 0.98, average ratio 1.0) with observed σ sp,525 as shown in Fig. 2e.Derived MSEs were used to calculate σ sp,525 during the whole observation period using the formula σ sp,525 = 2.1 × HOA + 3.9 × COA+3.4×LOOA+9.9×MOOA+5.5×AS+7.1×AN+3.3 × BC and compared with observed σ sp,525 .Good consistency (R 2 = 0.93 and average ratio of 1.05; Fig. 2f) was achieved between calculated and observed σ sp,525 values.In addition, the retrieved MSE AN (7.1) using the modified multilinear regression model was quite consistent with the estimated one (7 m 2 g −1 ) based on average changes shown in Fig. 2a, which indirectly confirms the reliability of the modified method.Tao et al. (2019) quantified MSEs of fine-mode AS and AN as well as elemental carbon (EC) using size-resolved filter measurements in four seasons of the Guangzhou urban area.Their results demonstrate that MSEs of AN and AS bear relatively large standard deviations and variations among seasons; however, the MSE of EC varied little among seasons with small standard deviations (2.6±0.1 m 2 g −1 ).The derived MSE BC of 3.3 m 2 g −1 was close to the MSE EC reported in Tao et al. (2019).The derived MSE AN and MSE AS were significantly higher than those reported in previous studies in which MSE AN and MSE AS were estimated through Mie theory of size-resolved filter measurements (Tao et al., 2019;Chen et al., 2020).For example, Tao et al. (2019) reported MSEs of 4.4 ± 1.3 m 2 g −1 for AN and 4.3 ± 0.9 m 2 g −1 for AS in winter of urban Guangzhou for fine-mode aerosols (< 2.1 µm).The reason explaining this inconsistency is that the derived MSE AN using the multiple regression method here is MSE * AN = σ sp,525 (AN,PM 2.5 ) [AN](PM 1 ) ; however, MSE AN derived for example in Tao et al. (2019) in the fine mode is MSE AN,PM 2.1 = σ sp,525 (AN,PM 2.1 ) [AN](PM 2.1 ) .The MSE AN,PM 2.1 , MSE AN,PM 1 and MSE * AN values of 4.4, 5.3 and 7.5 m 2 g −1 are simulated using the reported average AN mass size distributions reconstructed from size-resolved filter measurements in winter of urban Guangzhou by Tao et al. (2019) (as shown in Fig. 3a) as inputs of the Mie model.The simulated MSE * AN of 7.5 is very close to the retrieved MSE * AN and is much higher than simulated MSE AN,PM 1 due to substantial mass contributions of 1 to 2.1 µm as shown in Fig. 3a, demonstrating that good consistency between results of the multiple regression model and results of Tao et al. (2019) was achieved.
If using the simulated ratio MSE AN,PM 1 /MSE * AN to approximate MSE AN,PM 1 , MSE AS,PM 1 and MSE MOOA,PM 1 , values of 5.0, 3.9 and 6.9 m 2 g −1 for MSE AN,PM 1 , MSE AS,PM 1 and MSE MOOA,PM 1 would be obtained, which fall in the reported ranges of MSE AS and MSE AN (Tao et al., 2019); however, the high MSE MOOA,PM 1 needs further investigation.MSE is determined by the aerosol volume scattering efficiency (VSE) defined as aerosol scattering per unit aerosol volume and aerosol density (ρ a ) with MSE = VSE ρ a .High MSE * MOOA of 9.9 m 2 g −1 was retrieved; however, most of the difference between MSE * MOOA and MSE * AN might be explained by their density differences.Based on the Mie theory, aerosol scattering is in general proportional to aerosol volume (Kuang et al., 2018); thus the volume size distribution is a determining factor in VSE variations.The VSE of PM 1 under different unimodal volume lognormal distribution conditions with a refractive index of 1.53-10 −7 i was simulated and is shown in Fig. 3b.The approximated MSE AN,PM 1 and MSE MOOA,PM 1 of 5.0 and 6.9 m 2 g −1 corresponds to VSE AN,PM 1 and VSE MOOA,PM 1 of 8.6 and 9.7 m 2 µm −3 according to the aerosol densities discussed in Sect.2.4, falling within the VSE ranges of geometric mean diameter (D gv ) near 500 nm and geometric standard deviation (σ gv ) of 1.3 to 1.5.This result is consistent with conclusions of several previous studies that the MOOA with the highest oxygen state that has experienced complex chemical aging such as aqueous-phase reactions likely shares a similar shape of mass and volume size distribution with inorganic secondary aerosols (Kuang et al., 2020b;Wang et al., 2021), and this result rationalizes using the MSE AN,PM 1 /MSE * AN ratio to derive MSE AS,PM 1 and MSE MOOA,PM 1 .In addition, the aerosol refractive index also played a significant role in aerosol VSE variations (G.Zhao et al., 2019;Liu et al., 2013), and the high MSE of MOOA might also be related with the high real part of its refractive index.Laboratory results of Li et al. (2017) revealed enhanced light scattering of SOA formed through multiphase reactions due to an increase of the real part of the refractive index.Zhao et al. (2021) reported that the real part of the aerosol refractive index varied over a wide range (1.36 to 1.78) and in general increased with the mass fraction increase of organic aerosol, suggesting a generally higher real part of the refractive index of organic aerosol.In general, these results revealed strong scattering abilities of MOOA under dry state; however the size distribution and refractive index of MOOA need further comprehensive investigations.
Moreover, effective densities of HOA, COA and LOOA are near 1 g cm −3 , suggesting that VSEs of HOA, COA and LOOA are around or slightly higher than their corresponding MSEs.As shown in Fig. 3b, VSE PM 1 < 4 m 2 µm −3 means that D gv was generally lower than 250 nm; thus more than 99 % of aerosol mass resided in PM 1 under σ gv < 1.8.Therefore, derived MSEs of HOA, COA and LOOA can be treated as their corresponding MSE PM 1 values.Cai et al. (2020) reported average HOA and COA volume size distributions in urban Beijing using PMF techniques.They found (Fig. 7 in Cai et al., 2020) that HOA volume peaked near 200 nm, and COA volume size distribution showed bimodal characteristics, with the first mode peaking near 90 nm and the second mode peaking near 350 nm, yielding MSEs that share similar magnitude with the retrieved MSEs of HOA and COA.These results further confirmed the reliability of the newly proposed multiple regression method.

Water uptake abilities of organic aerosols
The time series of derived κ OA is shown in Fig. 4a, where the estimated κ OA ranged from −0.08 to near 0.35 with an average of 0.09, which is generally consistent with the values reported in other regions (Kuang et al., 2020d).It was found that variations of derived κ OA correlated tightly with mass fractions of MOOA (R = 0.6) and POA (R = −0.52),https://doi.org/10.5194/acp-22-7713-2022Atmos.Chem.Phys., 22, 7713-7726, 2022  as shown in Fig. 4b and c.The MOOA enhanced the overall κ OA , and POA lowered κ OA , which is consistent with conclusions of previous studies; however, no correlations were found between κ OA and mass fractions of LOOA (R = 0.06).
As shown in Fig. 4d and e, a drastic increase of POA before dusk would bring a drastic decrease of κ OA to around 0.05, which are in accordance with reported results in previous literature that most POA components are hydrophobic with a κ value of almost zero.As introduced in Sect.2.4, the average diurnal variations of volume fractions of MOOA and LOOA were used as variables to fit average diurnal variations of derived κ OA shown in Fig. 4e to retrieve κ MOOA and κ LOOA .The fitted results are shown in Fig. 4e (R 2 = 0.89, average ratio of 1), yielding κ MOOA and κ LOOA values of 0.23 and 0.13.The relatively lower value of κ LOOA and the phenomenon of sometimes co-increasing LOOA and POA during the dusk period (as shown in Figs. 1 and 3d) explain the weak correlations between κ OA and LOOA mass fraction.Most previous field-based organic aerosol hygroscopicity studies focused only on the overall κ OA characterization of entire organic aerosol population and were rarely specific to secondary organic aerosol factors (Kuang et al., 2020d).Considering the estimated O/C ratios of 1.17 and 0.7 for MOOA and LOOA, the retrieved κ LOOA and κ MOOA values are consistent with previous findings that organic aerosol oxidation state impacts significantly on organic aerosol hygroscopicity and usually higher hygroscopicity of more oxygenated organic aerosol (Kuang et al., 2020d).In particular, Lambe et al. (2011) investigated the relationships between organic aerosol hygroscopicity of laboratory-generated SOA of varying kinds of volatile organic compounds using cloud condensation nuclei activity measurements and reported a linear relationship of κ OA = (0.18 ± 0.04) × O/C + 0.03, yielding κ OA values of 0.24 and 0.16 for O/C of 1.17 and 0.7, which is slight higher but generally consistent with retrieved ones of this study.Overall, the retrieved κ LOOA and κ MOOA are further verified indirectly through comparison with previous literature results, confirming the strongest water uptake abilities of MOOA among OA factors.

Dominant contribution of MOOA to organic aerosol scattering ability
High scattering efficiency and water uptake abilities of MOOA resulted in the strongest light scattering abilities of MOOA among all organic aerosol factors.As shown in Fig. 5a and b, on average, a 34 % mass contribution of MOOA to entire PM 1 organic aerosol populations, however, contributed 51 % of organic aerosol scattering in dry state.
The dominant role of MOOA in organic aerosol scattering would be further enhanced under ambient RH conditions due to the highest water uptake abilities of MOOA among organic aerosol factors.Results of Kuang et al. (2017) demonstrate that the hygroscopicity parameter κ can be directly linked to the optical hygroscopicity parameter κ sca through κ sca = κ × R sca and thus aerosol light scattering enhancement factor f (RH).Particle number size distributions plays the most important role in R sca variations, with κ playing the smaller role.Yu et al. (2018) investigated R sca variations from measurements of several field campaigns and found that R sca varied within the range of 0.55 to 0.8 with an average of 0.66 and parameterized R sca with scattering Ångström exponent.Here, the relationship between VSE PM 1 and R sca was further simulated using Mie theory through varying the volume geometric mean D gv of lognormal size distributions from 100 to 700 nm considering that aerosol size distributions also play the dominant role in VSE PM 1 variations.Simulated results are shown in Fig.
The problem remains as to whether a continuous increase of MSE as a function of ambient RH increase can be applied because aerosol phase states are also crucial in determining the responses of aerosol scattering to RH increases (Kuang et al., 2016).Many factors such as ambient RH (Liu et al., 2017b(Liu et al., , 2016) ) and the RH history that aerosol particles have experienced (Kuang et al., 2016) and the deliquescent and crystalline properties determined by aerosol chemical compositions and mixing states are important in determining aerosol phase state (Li et al., 2021).The ambient RH ranged from 20 to 94, with an average of 57 %, with the histogram of ambient RH also shown in Fig. S12.Results of Liu et al. (2017b) found a transition from a semisolid to a liquid state at RH of about 60 %, suggesting ambient aerosol particle might exist as semisolid or solid phase at low RH conditions.The phase state of ambient aerosols depends not only on ambient RH, but also on the RH history that aerosols have experienced (Kuang et al., 2016).For instance, aerosols under ambient RH of 40 % in the afternoon would remain liquid if their crystalline RH is lower than 40 % and they had experienced high RH conditions (such as > 80 %) during the morning.Therefore, the lowest RH that aerosols have experienced in the afternoon and highest RH they have experienced in the morning are crucial for their phase state for hydrophilic aerosols.In addition, the deliquescent RH and crystalline RH are another two crucial parameters; however, they are quite complex for ambient multicomponent aerosols (Kuang et al., 2016;Li et al., 2021)  for the humidified nephelometer system from 13 to 26 January, and a continuous increase of the aerosol light scattering enhancement factor was always observed (Fig. S13), with the RH in the dry nephelometer always lower than 35 % (Fig. 1a), suggesting that aerosols were either not crystalized under RH of < 35 % or were deliquesced under RH of < 60 %.The lines of 35 % and 60 % are plotted in Fig. 1a, and most days (> 85 %) either had its lowest RH > 35 % or its highest RH > 60 %, suggesting a liquid state most of the time for internally mixed ambient secondary aerosols and a continuous increase of aerosol light scattering as ambient RH changes.Therefore, continuous increases of MSEs of LOOA, MOOA, AN and AS as a function of RH were considered, and results are shown in Fig. 5d.MSE PM 1 of MOOA changed from 6.9 m 2 g −1 under dry conditions to 11 and 16 m 2 g −1 (corresponding to f (RH) values of 1.6 and 2.3) under RH of 80 % and 90 %, revealing that scattering abilities of MOOA would be largely enhanced by aerosol hygroscopic growth.The MSE PM 1 of LOOA was enhanced from 3.4 under dry conditions to 7 m 2 g −1 of 90 % RH; however, MSE PM 1 values of HOA and COA remained constant due to their hydrophobic properties.Both MSE PM 1 values of AS and AN increased quickly as a function of RH, and their f (RH) values reached as high as 2.5 and 2.2 at RH of 80 %.The MSE PM 1 of AN exceeded that of MOOA at RH near 70 %; however, MSE PM 1 of MOOA was always higher than that of AS for RH < 90 %.Average MSE PM 1 values of secondary aerosol components considering water uptake under ambient RH conditions during the observation period are 6.8, 8.5, 8.9 and 4.2 m 2 g −1 for AS, AN, MOOA and LOOA, respectively, demonstrating the strongest scattering abilities of MOOA under ambient air conditions during the observation period.Thus, the average contribution of MOOA to PM 1 organic aerosol scattering are further enhanced to 54 % under ambient conditions, as shown in Fig. 5c (campaign average RH ∼ 57 %).

Implications for visibility improvement and aerosol radiative effect simulations
Atmospheric visibility measures the maximum distance that people can see, which is determined by air extinction.The strong light scattering abilities of MOOA might have significant effects on atmospheric visibility and direct aerosol radiative effects and thus have broad implications for both aerosol environmental and climate effects.The contributions of different aerosol components to visibility degradation under different visibility conditions were estimated as fractional contributions to ambient atmospheric extinction caused by both aerosols and air molecules.The results are shown in Fig. 5e, and a detailed estimation method is introduced in Sect.S6.It shows that AN contributed most to visibility degradation especially under polluted conditions, which is consistent with findings of several recent studies that nitrate plays an increasing and even dominant role in visibility degradation (Liu et al., 2020;Li et al., 2022) in several regions of China.MOOAs contributed slightly less than ammonium sulfate, and their contributions increased to ∼ 20 % for a visibility range of 10-20 km.Contribution of MOOA to total organic aerosol scattering under ambient conditions was slightly higher than that under dry conditions due to the relatively low RH conditions during winter.However, the contributions of MOOA might be much higher during other seasons due to much higher RH conditions, for example, yearly average RH of > 70 % in Guangzhou (Xu et al., 2020).Overall, secondary aerosols contributed dominantly to visibility degradation (∼ 70 % on average), and MOOA represented the third contributor among secondary aerosol components (16 % on average), demonstrating significant impacts of MOOA on visibility degradation.Thus, more attention should be paid to property changes of SOA regarding visibility improvement investigations and policymaking.Moreover, MOOA with high scattering abilities would likely contribute substantially to aerosol optical depth, and the accurate estimations of organic aerosol radiative effects in models need accurate representations of MOOA, even its mass contribution to organic aerosol is small.However, constant scattering efficiencies for organic aerosols derived from fixed size distributions and the refractive index were usually assumed in current chemical models (Latimer and Martin, 2019), and unsatisfactory performances of current models in SOA simulations were generally reported, which hinders the accurate representations of direct aerosol radiative effects.The role of MOOA will likely become more important in future due to stringent control policies of precursors of inorganic aerosol components such as sulfur dioxide and nitrogen oxides.Both formation pathways and precursor sources of MOOA are complex; however, current understanding of physical and chemical properties as well as formation pathways of MOOA remains limited.Therefore, we call for more comprehensive studies on formation, evolution and physical properties of MOOA to better parameterize optical properties of MOOA in models and implement targeted control strategies on MOOA precursors of volatile organic compounds in the future.
Data availability.All data needed are presented in figures, tables and supplementary figures.Raw datasets of this study are available from the first author and corresponding author upon request.
Supplement.Supporting information on site information, Q-ACSM PMF analysis, the traditional multiple linear regression model analysis and the visibility contribution estimation method, including two supporting tables and 13 supporting figures, is given in the Supplement.The supplement related to this article is available online at: https://doi.org/10.5194/acp-22-7713-2022-supplement.
Author contributions.YK and LL designed the aerosol experiments.YK conceived and led this research.LL and YK wrote the manuscript.MMZ, LL and CQY conducted the Q-ACSM measurements.MMZ and YH performed the PMF analysis.BX, BL, HBX, FL and YK performed aerosol optical property measurements.XJD and TD obtained funding for this campaign and provided instruments.JT, JCT, WYX, HBT and MS provided insights into data analysis, and all authors contributed to discussions and revisions of this paper.
Competing interests.The contact author has declared that neither they nor their coauthors have any competing interests.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Review statement.This paper was edited by Sergey A. Nizkorodov and reviewed by two anonymous referees.

Figure 1 .
Figure 1.Time series of (a) RH; (b) aerosol scattering coefficient in dry state, and NR PM 1 concentration; (c) secondary inorganic aerosols (SIAs), including sulfate and nitrate, and organic aerosol in right axis; (d) MOOA and LOOA; (e) HOA and COA; and (f) wind speed and wind direction.Gray shaded areas are three identified pollution episodes, and the pink shaded area is the period of the spring festival.
trometer.This method mainly considers the responses of the aerosol scattering coefficient to quick mass concentration increases of aerosol components.Using AN as an example, cases of significant increases of AN were identified, and average changes of aerosol components as well as σ sp,525 for these cases are shown in Fig.2a.On average, AN dominated the aerosol mass increase (> 90 %) in these cases; however, changes of other aerosol components differed much among cases, as indicated by large standard deviations.The MSE AN can be roughly estimated as around 7 m 2 g −1 if assuming σ sp,525 was solely contributed by AN increase.As shown in Fig.1d and e, prominent increases of HOA, COA and LOOA were frequently observed.Average changes of aerosol components and σ sp,525 for identified cases of HOA increase or COA increase are shown in Fig. 2b.This shows that increases of HOA or COA were usually accompanied by significant increases of BC and LOOA; thus the impacts of HOA or COA increases on observed aerosol scattering increases cannot be isolated.Similar results were obtained with LOOA and MOOA.As shown in Fig. 1, remarkable inhttps://doi.org/10.5194/acp-22-7713-2022Atmos.Chem.Phys., 22, 7713-7726, 2022

Figure 2 .
Figure 2. Average changes of aerosol components and σ sp,525 (right axis) for identified increase cases of (a) nitrate,; (b) HOA, (c) LOOA and (d) MOOA.Black error bars represent standard deviations.(e) Comparisons between observed σ sp,525 changes for identified cases and multiple linear fitted values.(f) The comparison between observed σ sp,525 and calculated σ sp,525 using retrieved MSEs of aerosol components.Dashed red lines represent 20 % deviation lines.

Figure 3 .
Figure 3. (a) AN mass size distributions derived by Tao et al. (2019) based on size-resolved filter measurements, Da is the aerodynamic diameter.(b) Simulated aerosol volume scattering efficiency (VSE) under different volume size distributions through varying volume geometric mean D gv and standard deviation σ gv of lognormal size distributions.The blue shaded area corresponds to PM 1 , and the pink shaded area corresponds to PM 1−2.1 .

Figure 4 .
Figure 4. Time series of derived κ OA (a).Correlations between derived κ OA and mass fractions of MOOA (b) and HOA+COA (c) in total organic aerosols.(d) Average diurnal variations of derived κ OA (blue) and corresponding fitted ones.(e) Average diurnal variations of resolved mass concentrations of organic aerosol factors.

Figure 5 .
Figure 5. (a) Average mass fractions of organic aerosol components in total organic aerosol.Contributions of different organic aerosol components to total organic aerosol scattering coefficients under dry (b) and ambient (c) conditions for PM 1 .(d) MSEs of aerosol components under different RH conditions.(e) Contributions of different aerosol components to visibility degradation under different visibility conditions.SA corresponds to summation of secondary aerosol components.

Acknowledgements.
This work is supported by the National Key Research and Development Program of China (2019YFCO214605 and 2016YFC0202000), the National Natural Science Foundation of China (41805109 and 42105092), the Guangdong Basic and Applied Basic Research Foundation (2019A1515110791), the Science and Technology Innovation Team Plan of Guangdong Meteorological Bureau (grant no.GRMCTD202003), and the Natural Science Foundation of Fujian Province, China (2021J01463).

Financial support .
This research has been supported by the National Key Research and Development Program of China (grant nos.2019YFC0214605 and 2016YFC0202000), the National Natural Science Foundation of China (grant nos.41805109 and 42105092), the Guangdong Basic and Applied Basic Research Foundation (grant no.2019A1515110791), the Science and Technology Innovation Team Plan of Guangdong Meteorological Bureau (grant no.GRMCTD202003) and the Natural Science Foundation of Fujian Province, China (grant no.2021J01463).