1.Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871, China 2.Institute for Environmental and Climate Research, Jinan University, Guangzhou 510632, China Manuscript received: 2018-11-19 Manuscript revised: 2019-03-27 Manuscript accepted: 2019-05-16 Abstract:Assessment of the radiative forcing of aerosols in models still lacks sufficient input data for aerosol hygroscopicity. The light scattering enhancement factor [f( RH,Λ)] is a crucial parameter for describing aerosol hygroscopic growth properties. In this paper, we provide a survey of f( RH,Λ) studies in China for the past seven years, including instrument developments of humidified nephelometers, ambient f( RH,Λ) measurements in China, f( RH,Λ) parameterization schemes, and f( RH,Λ) applications in aerosol measurements. Comparisons of different f( RH,Λ) parameterizations are carried out to check their performance in China using field measurement datasets. We also summary the parameterization schemes for predicting f( RH,Λ) with aerosol chemical compositions. The recently developed methods to observe other aerosol properties using f( RH,Λ) measurements, such as calculating the aerosol hygroscopicity parameter, cloud condensation nuclei number concentration, aerosol liquid water content, and aerosol asymmetry factor, are introduced. Suggestions for further research on f( RH,Λ) in China are given. Keywords: aerosol, hygroscopicity, light-scattering enhancement factor, humidified nephelometer 摘要:气溶胶吸湿性是模式中计算气溶胶辐射强迫的关键参数. 光散射增强因子 [f(RH,Λ)] 是描述气溶胶吸湿增长特性的主要参数之一. 本文对中国过去7年的f(RH,Λ)研究做了全面的综述,包括加湿浊度计的仪器研发, 中国地区外场f(RH,Λ)的观测, f(RH,Λ)的参数化方案, 和f(RH,Λ)在气溶胶特性观测方面的应用. 利用外场观测数据对于各种f(RH,Λ)的参数化方案进行了检验, 同时总结了各种利用气溶胶化学组成进行f(RH,Λ)参数化的方案. 本文介绍了最近利用f(RH,Λ)发展的观测其他气溶胶变量的新方法, 如气溶胶吸湿性因子, 云凝结核浓度, 气溶胶液态含水量和气溶胶不对称因子. 最后本文给出了未来中国f(RH,Λ) 研究的建议. 关键词:气溶胶, 吸湿性, 光散射增强因子, 加湿浊度计
HTML
--> --> --> -->
4.1. Applications of f( RH,λ) parameterizations in China in recent years
Fitting results in China are described here. One-parameter equations are simple and commonly used. (Yu et al., 2018) compared the performances of Eq. (3) and Eq. (6) to fit the measured f( RH,525) curves. The results showed that the deviations of the fitted f(55%,525), f(70%,525), f(85%,525) by Eq. (6) were all within 10%, while the fitted deviations of Eq. (3) can reach 20%, suggesting Eq. (6) should be used to fit f( RH,λ) curves. With additional parameters in the equations, two-parameter equations are introduced to describe f( RH,λ) curves more precisely. (Yang et al., 2015) used Eq. (9) to fit f( RH,550) curves, and the square of the correlation coefficient (R2) between the fitted and measured f( RH,550) values was 0.87. (Liu et al., 2013) introduced a two-parameter equation, Eq. (9), to fit f( RH,525) curves, and the R2 between the fitted and measured f( RH,525) values was 0.93. (Kuang et al., 2016) proved that f( RH,550) curves without deliquescent phenomena and f( RH,550) curves with RHs higher than deliquescence points, can be fitted well by Eq. (9). (Zhang et al., 2015) applied Eq. (7) and Eq. (9) to f( RH,550) curves. The fitted parameters of Eq. (7) and Eq. (9) changed along with local pollution and dust episodes. For Eq. (7), the results showed that larger f( RH,λ) values resulted in bigger values of parameters c and γ. Moreover, the fitted parameters of the measured f( RH,550) were found to be related to the mass fractions of nitrate. Similarly, (Qi et al., 2018) also utilized both Eq. (7) and Eq. (9) to fit measured f( RH,550). The results showed that the R2 of Eq. (7) was 0.42 and the R2 of Eq. (9) was 0.35. They found that parameter γ in Eq. (7) was more sensitive to the variations of f( RH,550), and the values of parameters a and b in Eq. (9) grew with the increase in aerosol hygroscopicity. (Zhao et al., 2018b) assessed the ability of Eqs. (5), (7), (8) and (9) to describe f( RH,525) curves and, ultimately, according to the calculated R2 between the fitted and measured f( RH,525) values, Eq. (8) was selected to fit f( RH,525) curves. In winter, the R2 values between the fitted and measured f( RH,525) values were 0.88, 0.88 and 0.89 under very clean, moderately clean and polluted conditions, respectively. Meanwhile, the R2 values were 0.89, 0.94 and 0.96 for summer, but 0.89, 0.94 and 0.98 for autumn, respectively. (Deng et al., 2016) and (Wu et al., 2017) both utilized the three-parameter Eq. (10) to fit measured f( RH,525) and f( RH,520) curves, respectively. All the fitting results are illustrated in detail in Table 1.
2 4.2. Comparisons of different f( RH,λ) parameterizations -->
4.2. Comparisons of different f( RH,λ) parameterizations
In order to access the ability of these proposed parameterizations to fit f( RH,λ) curves, the curves of f( RH,450), f( RH,525) and f( RH,635) measured at Gucheng in 2016 and illustrated in (Yu et al., 2018) and Zhangqiu in 2017, expressed in section 3, are also introduced here. The fitting results of the data from Gucheng and Zhangqiu are illustrated in Tables 2 and 3, respectively. The two-parameter equations can better express f( RH,λ) curves than the single-parameter equations, while the three-parameter equation does not show any obvious advantage. As shown in Table 2, for the Gucheng data introduced in (Yu et al., 2018), Eq. (8) describes the f( RH,λ) curves with the highest R2 and lowest root-mean-square error (RMSE) for all three wavelengths. The corresponding R2 values at 450 nm, 525 nm and 635 nm are 0.993, 0.994 and 0.995, respectively. The R2 values of the three-parameter Eq. (10) at 450 nm, 525 nm and 635 nm are 0.990, 0.991 and 0.992, respectively, which are lower than those of Eq. (7) or (8). The R2 values of single-parameter schemes are the smallest, but still higher than 0.92. Equation (5) is the best single-parameter equation, with corresponding R2 values at 450 nm, 525 nm and 635 nm of 0.975, 0.974 and 0.972, respectively. Table 3 shows that, for the data from Zhangqiu, Eq. (10) is the best, and the R2 of the two-parameter schemes is larger than that of the single-parameter schemes. For the single-parameter schemes, the R2 values are all higher than 0.92, and Eq. (3) is the best single-parameter scheme, with R2 values of 0.970, 0.967 and 0.970 at 450 nm, 525 nm and 635 nm, respectively. For the two-parameter schemes, the R2 values are all higher than 0.98. Equation (7) is proven to be the best, and the R2 values at 450 nm, 525 nm and 635 nm are 0.988, 0.989 and 0.989, respectively. The observation sites at Gucheng and Zhangqiu are both in the NCP, surrounded by farmland and residential areas. The degree of pollution at these two sites is severe. As the dataset from Gucheng was collected during autumn and the dataset from Zhangqiu was collected during summer, these selected two datasets represent the background conditions of the NCP in different seasons. Therefore, conclusions drawn from these two datasets have a certain applicability in the NCP, but not all locations in China. As for other locations in China, further studies are needed to evaluate the performance of different fitting equations.
-->
5.1. Parameterizing f( RH,λ) at a certain RH with aerosol chemical compositions
(Zhang et al., 2015) sought the relationships between f( RH,λ) and the mass fraction of organic matter (OM) (Forg) or inorganic matter (Fio). Results showed that the statistical relationship between f(85%,550) and Forg was f(85%,550)=(2.05 0.02)-(1.20 0.04)Forg, while the statistical relationship between f(85%,550) and Fio was f(85%)=(1.10 0.01)+(0.96 0.02)Fio. The study of (Wu et al., 2017) indicated that f(80%,520) increased with an increase in secondary inorganic aerosols mass fraction (Fsia), and the R between them was 0.758. The corrsponding statistical relaitionship was calculated to be f(80%,520)=0.981+3.502Fsia. Meanwhile, f(80%,520) decreased with an increase in Forg, with an R of -0.679, and the relationship was f(80%,520)=2.813-1.955Forg. (Qi et al., 2018) quantified the relationship between f(80%,550) and different aerosol chemical compositions. The results showed that f(80%,550) was negatively related to OM, and the R2 between them was 0.55. The calculated relationship was found to be f(80%,550)=-(0.99 0.12)Forg+(1.70 0.05). The f(80%,550) was also positively correlated with Fio or the nitrate mass fraction (Fni), and the R2 between both was 0.70. The corresponding relationships were f(80%,550)=(1.36 0.12)Fio+(0.82 0.05) and f(80%,550)=(2.465 0.217)Fni+(1.027 0.024).
2 5.2. Parameterizing the fitted parameters of f( RH,λ) with aerosol chemical compositions -->
5.2. Parameterizing the fitted parameters of f( RH,λ) with aerosol chemical compositions
Section 5.1 represents the calculated relationships between f( RH,λ) at a fixed RH condition with aerosol chemical compositions. However, it is critical to obtain the aerosol scattering properties at any RH condition to calculate the ambient aerosol optical properties. As the f( RH,λ) curves are usually described by some appropriate equations, several studies have been carried out to explore the relationships between aerosol chemical compositions and the fitted parameters of f( RH,λ). (Zhang et al., 2015) employed the method of (Quinn et al., 2005) to study the relationship between Forg and the parameter γ in Eq. (3), where the Forg was calculated by Forg= OM/( OM+ SO42-). The resulting statistical relationship was γ=-0.25Forg+0.48. However, this scheme is only applicable to the aerosol types dominated by sulfates and organics, and the R2 between γ and Forg was only 0.14. A further study by (Zhang et al., 2015) found that nitrate also made an important contribution to the hygroscopic growth of aerosol. The statistical relationship was γ=-0.26Forg+0.50, and the R2 between Forg and γ was 0.56 when the Forg was described as Forg= OM/( OM+ SO42-). Furthermore, if the Forg was set as Forg= OM/( OM+ SO3-+ SO42-), the relationship was γ=-0.42Forg+0.54 and the R2 between Forg and γ was 0.68. (Yu et al., 2018) newly proposed a scheme to bridge the gap between aerosol chemical compositions and f( RH,λ) through Eq. (5). They first linked the aerosol hygroscopicity parameter $\kappa$ and chemical compositions with three dominant inorganic ions ( NH4+, NO3-, SO42-) and OM. Then, the relationship among $\kappa$, the fitted parameter $\kappa_\rm sca$ of f( RH,λ), and the scattering ?ngstr?m exponent (SAE), which represents the size of aerosol, was established. Thus, the hygroscopic parameter $\kappa$ can be determined through aerosol chemical compositions and then the fitted parameter $\kappa_\rm sca$ can be obtained by $\kappa$ and SAE. The resulting empirical equation to parameterize f( RH,λ) with aerosol chemical compositions is shown as \begin{eqnarray} f({\rm RH},525)&=&1+(0.01+0.70 f_{{\rm NH}_{4}}+0.44f_{{\rm NO}_{3}^-}+0.62f_{{\rm SO}_{4}^{2-}}+ \\ &&0.06f_{{\rm OM}})(0.45+0.15{\rm SAE})\frac{{\rm RH}}{1-{\rm RH}} . \ \ (11)\end{eqnarray} The verification results showed that the correlation coefficient of f(80%,525) obtained by observation and parameterization was 0.81. In addition, the method can further consider the influence of the aerosol size on f( RH,λ). This proposed scheme links chemical compositions to f( RH,λ) and can be applied to chemical transport models to reduce assessment errors of aerosol direct radiative forcing and atmospheric visibility.
-->
6.1. Calculating the aerosol hygroscopicity parameter $\kappa$
The RH dependence of the growth of an aerosol particle owing to water uptake can be parameterized in a good approximation by a one-parameter equation——proposed, for example, by (Petters and Kreidenweis, 2007). The parameter $\kappa$ is a simple measure of the particle's hygroscopicity and captures all solute properties. Traditionally, an overall hygroscopicity parameter $\kappa$ can be retrieved from measured f( RH,λ), hereinafter referred to as $\kappa_f(\rm RH)$, by combining concurrently measured PNSDs and mass concentrations of black carbon. (Kuang et al., 2017) proposed a new method to directly derive $\kappa_f(\rm RH)$ based only on measurements from a three-wavelength humidified nephelometer system. The advantage of this newly proposed approach is that $\kappa_f(\rm RH)$ can be estimated without any additional information on PNSDs and black carbon. (Kuang et al., 2017) verified this method with measurements from different field campaigns conducted in the NCP, and their results demonstrated that this method of deriving $\kappa_f(\rm RH)$ is applicable at different sites and in different seasons of the NCP and might also be applicable in other regions around the world. This work directly links f( RH,λ) to $\kappa$ for the first time, which is a breakthrough for studying the impacts of aerosol hygroscopic growth on aerosol optical properties and connecting aerosol optical properties to aerosol water content. Thus, their results should make the humidified nephelometer system more convenient when it comes to aerosol hygroscopicity research, as well as facilitate other research on the roles of aerosol water in aerosol radiative effects.
2 6.2. Calculating ambient ALWC -->
6.2. Calculating ambient ALWC
Aerosol liquid water plays significant roles in the atmospheric environment, atmospheric chemistry, and climate. Unfortunately, until now, no instruments have been available for real-time monitoring of ambient ALWC. The challenging issue hindering the development of ambient ALWC measurement techniques is that the amount of ambient ALWC is quite small and very sensitive to RH changes, which makes it unable to be measured directly. (Kuang et al., 2018) proposed a novel method to calculate ambient ALWC based on measurements of a three-wavelength humidified nephelometer system. The proposed ALWC calculation method includes two steps. The first step is calculating the dry state total volume concentration of ambient aerosol particles, Va( RHdry), with a machine learning model based on measurements of the "dry" nephelometer. The second step is calculating the volume growth factor Vg( RH) of ambient aerosol particles due to water uptake, using measured f( RH,λ) and SAE. Then, ambient ALWC can be derived based on the calculated Va( dry) and Vg( RH). (Kuang et al., 2018) validated the proposed ALWC calculation method using ambient ALWC calculated from the aerosol thermodynamic model of ISORROPIA with measured aerosol chemistry data. Their results demonstrated that a good agreement was achieved between the ALWC calculated from measurements of the humidified nephelometer system and from those estimated using the ISORROPIA model, with a slope and intercept of 1.14 and -8.6 μm3 cm-1 (R2=0.92), respectively. The thermodynamic model needs chemistry data as inputs, which requires expensive chemical instruments, and has a relatively low temporal resolution (1 h). The thermodynamic model only accounts for the ALWC contribution from inorganic aerosol components and is unable to take into account that from organic matter. Even when measurements of aerosol organic matter are available, accurate estimation of their hygroscopicity still remains unresolved. The advantage of this newly proposed method is that the required measurement data can be obtained solely from the humidified nephelometer system, which is quite stable and has a high temporal resolution. Also, this method considers the contributions of both inorganic and organic aerosol components to ambient ALWC, facilitating the real-time monitoring of ambient ALWC and promoting the study of aerosol liquid water and its role in atmospheric chemistry, secondary aerosol formation, and climate change.
2 6.3. Calculating number concentrations of CCN -->
6.3. Calculating number concentrations of CCN
Based on measurements of a humidified nephelometer system, (Tao et al., 2018) proposed a new method to calculate the number concentration of CCN (NCCN), which is a key parameter of cloud microphysics and the indirect radiative effect of aerosol. In general, NCCN is directly measured under supersaturated conditions in CCN chambers, which are complex and costly. As accumulation-mode aerosols contribute most to both aerosol optical properties and the aerosol CCN activity, NCCN can be predicted by its relationship with the aerosol scattering coefficient. In this new method, a look-up table that involves the scattering coefficient σsp, SAE, and hygroscopicity parameter $\kappa$, is established to derive NCCN based on measurements of a three-wavelength humidified nephelometer system that can measure the required three parameters (i.e., σsp, SAE and $\kappa$). This method has been validated by comparing the measurements of a humidified nephelometer system and a CCN counter in Gucheng in 2016 (Tao et al., 2018). For the calculated and measured NCCN, good agreements can be achieved. Relative deviations are within 30%, and the slope and R of the regression are 1.03 and 0.966, respectively. Because the humidified nephelometer system is simply operated and stable, the real-time monitoring of NCCN, especially on aircraft, can be facilitated by this new method. Furthermore, this method is more applicable for studies of aerosol-cloud interaction, due to its applicability at lower supersaturations than 0.1%.
2 6.4. Calculating the asymmetry factor of aerosol -->
6.4. Calculating the asymmetry factor of aerosol
In addition to aerosol optical depth and aerosol single-scattering albedo, the aerosol phase function is the most important factor for assessing direct aerosol radiative forcing. However, little attention has been paid to the measurements and parameterization of the asymmetry factor g. (Zhao et al., 2018a) proposed a novel method to calculate g based on measurements from the humidified nephelometer system. This method constrains the uncertainty of g within 2.56% for dry aerosol populations and 4.02% for ambient aerosols, where the aerosol hygroscopic growth has been taken into account. The total uncertainty of the calculation of g using the Random Forest machine learning model is 4.47%. Sensitivity studies show that aerosol hygroscopicity plays a vital role in the accuracy of predicting g. This new method for calculating g has been validated by comparing the values of g from the Random Forest machine learning model and those from field-measured phase function (Zhao et al., 2018a). The g values with these two methods show good consistency, with 95% of the data within the relative difference of 6.5%.