Quantifying the effects of mixing state on aerosol optical properties

Abstract

Calculations of the aerosol direct effect on climate rely on simulated aerosol fields. The model representation of aerosol mixing state potentially introduces large uncertainties into these calculations, since the simulated aerosol optical properties are sensitive to mixing state. In this study, we systematically quantified the impact of aerosol mixing state on aerosol optical properties using an ensemble of 1800 aerosol populations from particle-resolved simulations as a basis for Mie calculations for optical properties. Assuming the aerosol to be internally mixed within prescribed size bins caused overestimations of aerosol absorptivity and underestimations of aerosol scattering. Together, these led to errors in the populations' single scattering albedo of up to $-\mathrm{22.3}$ % with a median of $-\mathrm{0.9}$ %. The mixing state metric $\mathit{\chi }$ proved useful in relating errors in the volume absorption coefficient, the volume scattering coefficient and the single scattering albedo to the degree of internally mixing of the aerosol, with larger errors being associated with more external mixtures. At the same time, a range of errors existed for any given value of $\mathit{\chi }$. We attributed this range to the extent to which the internal mixture assumption distorted the particles' black carbon content and the refractive index of the particle coatings. Both can vary for populations with the same value of $\mathit{\chi }$. These results are further evidence of the important yet complicated role of mixing state in calculating aerosol optical properties.

1 Introduction

Particles scatter and absorb incoming solar radiation, thereby impacting the global radiative balance and temperatures on Earth . Black carbon (BC), commonly emitted from combustion, has a direct radiative forcing of $+\mathrm{0.11}$ $\mathrm{W}{\mathrm{m}}^{-\mathrm{2}}$ . At the same time, the overall global average aerosol direct radiative forcing in the clear-sky environment is $-\mathrm{0.22}$ $\mathrm{W}{\mathrm{m}}^{-\mathrm{2}}$ because of the presence of other non-absorbing aerosol species, which exert a cooling impact .

Radiative effects of aerosols depend on their optical properties, which, as a whole, are determined by the individual particles that the aerosol consists of. As observed in field campaigns, particles are mixtures of inorganic and organic species and exhibit significant spacial and temporal variation in their abundance and composition , with considerable diversity in composition existing within individual aerosol populations. The topic of this paper is to quantify the importance of diversity in composition for aerosol optical properties.

Aerosol composition impacts aerosol optical properties for several reasons. First, aerosol species differ in their complex refractive index. While inorganic species and many organic species have a purely real refractive index for wavelength of visible sunlight (i.e., only scatter radiation), black carbon and some organic carbon species have a nonzero imaginary part of the refractive index and hence also absorb radiation . Second, aerosol species differ in their hygroscopicity. This governs aerosol water uptake in a humidified environment, which is important for scattering .

Lastly, the arrangement of the different aerosol species within a particle is important for determining their scattering and absorption. For mixed particles without strongly absorbing species, i.e., BC, a volume-mixing rule can be used to calculate the refractive index of the entire particle. When the particle contains BC, assuming a core–shell configuration was proven to be more accurate (still assuming sphericity as particle shape). The absorption enhancement of BC-containing particles due to their non-absorbing coatings has been widely investigated . Taking the nonspherical shapes of BC-containing particles into account complicates matter considerably since Mie calculations cannot be applied and more sophisticated optical models need to be used, which are computationally much more expensive. Using the discrete dipole approximation model, found that the absorption coefficients' enhancement of BC / NaCl mixtures is higher for compact BC particles completely embedded in NaCl than for lacy BC particles.

To understand the importance of aerosol composition in calculating aerosol optical properties, it is useful to define the term aerosol mixing state, that is, the distribution of aerosol species among the particles in a population . Aerosol mixing state in the ambient atmosphere ranges between the two idealized extremes of an external mixture on the one hand, where each particle is composed of a single species, and an internal mixture on the other hand, where all particles consist of the same mixture of species. Aerosols close to emission sources tend to be more (although not completely) externally mixed . After emission, aging processes, such as coagulation between particles and condensation of gas species on the particles, transform aerosol populations towards more internal mixtures . Past studies quantified the importance of mixing state for aerosol optical properties through optical closure studies. For example, using measured aerosol size distributions and aerosol composition observed over the East China Sea, found that the internal mixture assumption for fine particles increased the absorption aerosol optical depth by a factor of 2 or more.

Aerosol mixing state is challenging to represent in 3D chemical transport models, which usually rely on simplifying assumptions for computational efficiency. These assumptions then influence the magnitude of calculated aerosol optical properties. Optical properties are here understood by three widely used parameters: the absorption cross section, the scattering cross section and the asymmetry parameter . Many 3D models use a modal approach to represent aerosols, such as the Community Multiscale Air Quality Modeling System (CMAQ) and the modal aerosol module (MAM) . The modes are externally mixed from each other, whereas within each mode, the aerosol is assumed to be internally mixed. For BC-containing modes, sphericity and a core–shell configuration are assumed, so that Mie calculations can be applied to calculate optical properties. found that neglecting the diversity in coating thickness for BC-containing particles (a result of the internal mixture assumption) leads to overestimated absorption enhancement by up to 200 %. Another approach is the sectional model representation, which tracks size-resolved composition but not particle composition diversity within a certain size, such as TwO-Moment Aerosol Sectional (TOMAS) and the GLObal Model of Aerosol Processes (GLOMAP) . Still, mixing state assumptions need to be invoked for each size bin. Recently, aerosol modules with more detailed BC mixing state representation were implemented in global climate models . These approaches better represent the evolution of BC aging processes within a size bin by adding a second dimension for BC mass fraction. However, this two-dimensional bin approach still does not capture the mixing state information of other, non-BC aerosol species.

The uncertainties in optical properties introduced by mixing state assumptions were also evaluated through model sensitivity studies. Using the AQMEII-2 model inter-comparison framework, quantified the sensitivity of aerosol optical properties to several parameters, including aerosol mixing state and size distribution. They found that aerosol mixing state is the dominant factor introducing uncertainties, explaining 30 %–35 % of the uncertainty in aerosol optical depth and single scattering albedo (SSA). found that the direct radiative forcing (DRF) can vary from $-\mathrm{1.65}$ to $-\mathrm{1.34}$ W m${}^{-\mathrm{2}}$ over the pan-Arctic region depending on the assumption of an internal or external mixture. The variation is similar when the assumptions are used to calculate DRF at the top of the atmosphere . These sensitivity studies have in common that no benchmark simulations exist that represent the real mixing state, and therefore the importance of mixing state can only be assessed based on differences between varied idealized assumptions. By applying a detailed particle-resolved benchmark model, found that simple mixing state assumptions can result in an erroneous distribution of BC cores and coating material and lead to errors in absorption. This effect was further confirmed to be the main source of the discrepancies between simulated and experimentally determined particle optical properties .

The goal of this study is to systematically quantify the errors in optical properties due to simplified assumptions for mixing state, here quantified with the mixing state metric $\mathit{\chi }$ . A similar framework was used to quantify the error in cloud condensation nuclei (CCN) concentration , showing that CCN error ranges from $-\mathrm{40}$ % to 150 % when assuming the aerosol was internally mixed. The error depended on the supersaturation level that CCN concentrations were evaluated at and also aerosol mixing state. In this work, we want to answer the following questions: given the aerosol mixing state, what is the error in aerosol optical values when assuming internal mixture, and what are the leading causes for this error?

The paper is structured as follows: model description, scenario design and the definition of metrics are given in Sect. . Section  shows the relation between the errors in aerosol scattering and absorption and mixing state for dry aerosol populations, and Sect.  further analyzes the errors for the aerosol populations at different levels of ambient relative humidity. The errors in single scattering albedo and its implications for aerosol direct radiative forcing are analyzed in Sect. . Section  summarizes the main findings.

2 Model description, scenario libraries and metrics

2.1 The stochastic particle-resolved model PartMC-MOSAIC

The model used for this study is the particle-resolved model PartMC-MOSAIC (Particle Monte Carlo model–Model for Simulating Aerosol Interactions and Chemistry). A comprehensive description of the model can be found in and for PartMC and in for MOSAIC. PartMC is a Lagrangian box model that tracks the evolution of particles in a fully mixed computational volume. The processes of emission, coagulation and dilution are simulated stochastically. Gas-phase chemistry and gas–aerosol partitioning are incorporated by coupling with the deterministic model MOSAIC. Specifically, MOSAIC uses the carbon-bond-based mechanism CBM-Z for gas-phase photochemical reactions , a multicomponent Taylor expansion method (MTEM) for calculating electrolyte activity coefficients in aqueous inorganic mixtures and a multicomponent equilibrium solver for aerosols (MESA) for calculating the phase states of the particles . The secondary organic aerosol (SOA) treatment follows the Secondary Organic Aerosol Model (SORGAM) . Aerosol water uptake is calculated using the Zdanovskii–Stokes–Robinson (ZSR) method based on the composition of the inorganic portion of the particles. By this method, organic species are treated as hydrophobic and do not contribute to water uptake. The impact of this assumption on optical properties was quantified by , who found that errors in single scattering albedo can be up to 6 % if neglecting the water uptake of organic compounds.

2.2 Scenario library design

Following the strategy in and , we created a scenario library of PartMC-MOSAIC simulations for this study, with a focus on the aging of carbonaceous aerosol. To produce particle populations with a wide range of compositions and mixing states, we varied the model input parameters within the ranges shown in Table . We used Latin hypercube sampling to create input parameter combinations for a total of 100 model simulations. The simulation time for each simulation was 24 h, beginning at 06:00 local time with hourly output. This yielded a total of 2500 particle populations. All scenarios were run with 10 000 computational particles. To create aerosol initial conditions with realistic mixing states, we adopted the approach described in : we carried out a first set of simulations, starting with the aerosol initial concentrations set to zero for all simulations (the “initial runs”). We then repeated the same set of simulations but replaced the aerosol initial condition with a randomly sampled population from the initial runs (the “restart runs”). For the analysis in this paper, we only used the results from the restart runs. Within our ensemble or aerosol populations, some were found with higher species concentrations than what would be expected in the ambient atmosphere. We applied upper thresholds to eliminate those which were calculated as the sum of the 75th percentile and 1.5 times the IQR (interquartile range) for each of the aerosol species. After this procedure, 1809 out of 2500 populations were used for the error analysis presented in the remainder of the paper.

Table 1

Baseline and range for the input variables.

Input parameters	Baseline	Range
Environment variables
Relative humidity (RH)		[0.1, 1) or [0.4, 1)
Latitude		($\mathrm{70}{}^{\circ }$ S, $\mathrm{70}{}^{\circ }$ N) or ($\mathrm{90}{}^{\circ }$ S, $\mathrm{90}{}^{\circ }$ N)
Day of year		[1, 365]
Temperature		Based on latitude and day of year
Gas emission rates (mol m${}^{-\mathrm{2}}$ s${}^{-\mathrm{1}}$)
Sulfur dioxide (SO${}_{\mathrm{2}}$)	$\mathrm{8.5}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Nitrogen dioxide (NO${}_{\mathrm{2}}$)	$\mathrm{3.0}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Nitrogen oxide (NO)	$\mathrm{5.7}\times {\mathrm{10}}^{-\mathrm{8}}$	[0 %–200 %]
Ammonia (NH${}_{\mathrm{3}}$)	$\mathrm{8.9}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Carbon oxide (CO)	$\mathrm{7.8}\times {\mathrm{10}}^{-\mathrm{7}}$	[0 %–200 %]
Methanol (CH3OH)	$\mathrm{2.3}\times {\mathrm{10}}^{-\mathrm{10}}$	[0 %–200 %]
Acetaldehyde (ALD2)	$\mathrm{1.7}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Ethanol (ANOL)	$\mathrm{5.3}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Acetone (AONE)	$\mathrm{7.8}\times {\mathrm{10}}^{-\mathrm{10}}$	[0 %–200 %]
Dimethyl sulfide (DMS)	$\mathrm{3.8}\times {\mathrm{10}}^{-\mathrm{11}}$	[0 %–200 %]
Ethene (ETH)	$\mathrm{1.8}\times {\mathrm{10}}^{-\mathrm{8}}$	[0 %–200 %]
Formaldehyde (HCHO)	$\mathrm{4.1}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Isoprene (ISOP)	$\mathrm{2.4}\times {\mathrm{10}}^{-\mathrm{10}}$	[0 %–200 %]
Internal olefin carbons (OLEI)	$\mathrm{5.9}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Terminal olefin carbons (OLET)	$\mathrm{5.9}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Paraffin carbon (PAR)	1$.7\times {\mathrm{10}}^{-\mathrm{7}}$	[0 %–200 %]
Toluene (TOL)	$\mathrm{6.1}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Xylene (XYL)	$\mathrm{5.6}\times {\mathrm{10}}^{-\mathrm{9}}$	[0 %–200 %]
Carbonaceous aerosol emission (single mode)
Geometric mean diameter ($D_{\mathrm{g}}$)		[25, 250] $\mathrm{nm}$
Geometric standard deviation of diameter (${\mathit{\sigma }}_{\mathrm{g}}$)		[1.4, 2.5]
$\mathrm{BC}/\mathrm{OC}$ mass ratio		[0, 100 %]
Particle emission flux		[0, $\mathrm{1.6}\times {\mathrm{10}}^{\mathrm{7}}$] m${}^{-\mathrm{2}}$ s${}^{-\mathrm{1}}$

2.3 Optical property calculations

We calculated the optical properties of the particle populations using Mie calculations . These properties included the asymmetry parameter $g$, scattering cross section ${\mathit{\sigma }}_{\mathrm{scat}}$ and absorption cross section ${\mathit{\sigma }}_{\mathrm{abs}}$ for each particle. Particles were assumed to be spherical, and when BC was present, a core–shell configuration was assumed, with BC as the core and non-BC species as the shell. In PartMC-MOSAIC, each chemical species was assigned a refractive index, and the values were the same as , as listed in Table . The shell refractive index of the particle was the volume average of all the shell species, including aerosol water. The absorptivity of brown carbon has been of great interest in recent years ; however, this was not considered in the current work. We used the values for wavelength $\mathit{\lambda }$ of 550 $\mathrm{nm}$ for our analysis.

Figure 1

Conceptual framework of composition averaging. The colors indicate different aerosol species. Light blue stands for water, and black stands for black carbon. The other colors conceptually represent other chemical species. In total, we track 18 aerosol species in addition to BC and water. Composition averaging is applied to the dry populations, and then water uptake is recalculated for RH $=\mathrm{50}\mathit{\%}$ and RH $=\mathrm{90}\mathit{\%}$. The composition averaging procedure conserves bulk mass concentration of each species, the total number concentration and the particle diameters within each bin.

[Figure omitted. See PDF]

In PartMC-MOSAIC, all particles are tracked individually in a well-mixed computational volume, and we obtained the ensemble optical property values by summing over all particles in the volume. The ensemble scattering coefficients ${\mathit{\beta }}_{\mathrm{scat}}\left(\mathit{\lambda }\right)$, ensemble extinction coefficients ${\mathit{\beta }}_{\mathrm{ext}}\left(\mathit{\lambda }\right)$ and ensemble absorption coefficients ${\mathit{\beta }}_{\mathrm{abs}}\left(\mathit{\lambda }\right)$ at wavelength $\mathit{\lambda }$ are given as $$\begin{array}{}\text{1}& {\displaystyle }{\mathit{\beta }}_{\mathrm{scat}}\left(\mathit{\lambda }\right)=\sum _i^N{\mathit{\sigma }}_{\mathrm{scat},i}\left(\mathit{\lambda }\right)n_i,\text{2}& {\displaystyle }{\mathit{\beta }}_{\mathrm{ext}}\left(\mathit{\lambda }\right)=\sum _i^N{\mathit{\sigma }}_{\mathrm{ext},i}\left(\mathit{\lambda }\right)n_i,\text{3}& {\displaystyle }{\mathit{\beta }}_{\mathrm{abs}}\left(\mathit{\lambda }\right)={\mathit{\beta }}_{\mathrm{ext}}\left(\mathit{\lambda }\right)-{\mathit{\beta }}_{\mathrm{scat}}\left(\mathit{\lambda }\right),\end{array}$$ where $i$ is the particle index, $n_i$ is the number concentration associated with particle $i$ and $N$ is the number of computational particles in the population. We determined the optical properties of all particle populations of our scenario libraries using these equations.

Two additional derived quantities of interest are the absorption enhancement and the mass absorption coefficient. The absorption enhancement of BC-containing particles due to coatings is defined as 4$$E_{\mathrm{abs}}\left(\mathit{\lambda }\right)={\mathit{\beta }}_{\mathrm{abs}}\left(\mathit{\lambda }\right)/{\mathit{\beta }}_{\mathrm{abs},\mathrm{BC}}\left(\mathit{\lambda }\right),$$ where ${\mathit{\beta }}_{\mathrm{abs},\mathrm{BC}}\left(\mathit{\lambda }\right)$ is the volume absorption coefficient when the particle coatings are removed from the BC cores.

We can also calculated the BC-specific mass absorption coefficients ${\mathrm{MAC}}_{\mathrm{BC}}$ (m${}^{\mathrm{2}}$ g${}^{-\mathrm{1}}$) using 5$${\mathrm{MAC}}_{\mathrm{BC}}={\displaystyle \frac{\sum _i^N{\mathit{\sigma }}_{\mathrm{abs},i}\left(\mathit{\lambda }\right)n_i}{\sum _i^N{m}_{\mathrm{BC},i}{n}_i}},$$ where $m_{\mathrm{BC},i}$ is the BC mass in particle $i$.

2.4 Quantifying the impact of mixing state through composition averaging

To quantify the impacts of mixing state on aerosol optical properties, we employed the strategy of “composition averaging” similar to to create sensitivity scenarios. The technique is shown conceptually in Fig. . For each population in our reference scenario library, we averaged the dry particle compositions within prescribed size bins. We chose eight size bins between 0.039 and 10 $\mathrm{\mu }$m, consistent with the bin structure of the sectional aerosol module MOSAIC used in WRF-Chem .

Table 2

Refractive indices of aerosol species at $\mathit{\lambda }=\mathrm{550}\mathrm{nm}$.

Compounds	Refractive index
H${}_{\mathrm{2}}$SO${}_{\mathrm{4}}$	1.43
(NH${}_{\mathrm{4}}$)${}_{\mathrm{2}}$SO${}_{\mathrm{4}}$	1.52
(NH${}_{\mathrm{4}}$)HSO${}_{\mathrm{4}}$	1.47
NH${}_{\mathrm{4}}$NO${}_{\mathrm{3}}$	1.5
H${}_{\mathrm{2}}$O	1.33
BC	1.82 + 0.74$i$
SOA	1.45
OC	1.45

The composition averaging procedure preserves the bulk mass concentration of each species, the total number concentration and the particle diameters within each bin Appendix B1; i.e., after composition averaging, each bin still contains particles of different sizes. It changes the per-particle compositions so that each bin becomes internally mixed; however the composition can vary between bins. This mimics the assumption frequently made in sectional models, namely that each size bin contains an internally mixed aerosol. PartMC-MOSAIC represents particles outside the MOSAIC bin range, especially for the lower boundary, and we used an extra bin (bin 0) to preserve the total number and mass concentrations. Since the aerosol water content plays an important role in aerosol optical properties, we further calculated water uptake for the reference populations and for the composition-averaged populations for 50 % ($P_{\mathrm{2}}$) and for 90 % ($P_{\mathrm{3}}$) relative humidity, respectively. At RH $=\mathrm{50}$ %, depending on the exact composition, some particles take up water, and at RH $=\mathrm{90}$ %, most particles take up water, except particles that only contain hydrophobic species, such as pure black carbon or primary organic carbon. Note that while the dry aerosol mass was conserved by the composition averaging procedure, the water content was recalculated after composition averaging and could change compared to the reference population.

Figure  illustrates the changes of two important parameters for aerosol optical properties due to composition averaging, BC mass fraction and the real part of the refractive index. In the reference case, a wide range of BC mass fractions exists within the same size bin (Fig. a). After composition averaging, all particles within a size bin have the same BC mass fraction (Fig. b). Since composition averaging preserves the particle diameters, BC and other species are redistributed so that all particles within a size bin are assigned the same mass fractions. Specifically, if a particle has lower BC mass fraction than the average level in the same size bin, BC is added to this particle from those with higher BC content. The coating species are also redistributed after composition averaging, which causes the refractive index of the coating to change (Fig. c and d). Hence, comparing optical properties before and after composition averaging in the dry population $P_{\mathrm{1}}$ isolates the impact of mixing state on aerosol optical properties. We will discuss the impact of composition averaging for dry conditions in Sect.  and the impact of water uptake in Sect. .

Figure 2

Two-dimensional number distributions of BC mass fraction and dry diameter: (a) reference and (b) composition-averaged. Real part of the refractive index and dry diameter: (c) reference and (d) composition-averaged. The population is taken from scenario 76 at $t=\mathrm{1}$ h, with $\mathit{\chi }=\mathrm{36}\mathit{\%}$. Red numbers and vertical gray lines represent the size bin ranges. The two red rectangles are for the analysis in Sect. .

[Figure omitted. See PDF]

2.5 Mixing state metrics

We quantified the optical properties error introduced by a simplified mixing state representation using the metrics developed by . These metrics include the single-particle diversity $D_i$, the average particle species diversity $D_{\mathit{\alpha }}$ and bulk population species diversity $D_{\mathit{\gamma }}$. For a population with $N$ particles, total mass $\mathit{\mu }$ and $A$ species, we can calculate these metrics from the total mass of particle $i$, ${\mathit{\mu }}_i$, total mass of species $a$ in the population, ${\mathit{\mu }}^a$, and mass of species $a$ in particle $i$, ${\mathit{\mu }}_i^a$, for $i=\mathrm{1}$, $\mathrm{\dots }$, $N$ and $a=\mathrm{1}$, $\mathrm{\dots }$, $A$. The mass fraction of species $a$ in particle $i$, $p_i^a$, mass fraction of particle $i$ in the population, $p_i$, and mass fraction of species $a$ in the population, $p^a$, are given by

6$$p_i^a={\displaystyle \frac{{\mathit{\mu }}_i^a}{{\mathit{\mu }}_i}},p_i={\displaystyle \frac{{\mathit{\mu }}_i}{\mathit{\mu }}},p^a={\displaystyle \frac{{\mathit{\mu }}^a}{\mathit{\mu }}}.$$

The single-particle diversity $D_i$ describes the effective species number in each particle and is defined as 7$$D_i=\prod _{a=\mathrm{1}}^A(p_i^a{)}^{-p_i^a}.$$ As shown in Fig. S1 in the Supplement, if a particle only contains one species, $D_i$ is 1. If the chemical species are present in equal amounts in the particle, $D_i$ equals the number of species. If the species are unevenly distributed, $D_i$ is a real number ranging between 1 and the number of species in the particle.

Based on $D_i$, we can construct $D_{\mathit{\alpha }}$ and $D_{\mathit{\gamma }}$, which describes the average effective species number in each particle and bulk population, respectively: $$\begin{array}{}\text{8}& {\displaystyle }{D}_{\mathit{\alpha }}=\prod _{i=\mathrm{1}}^N(D_i{)}^{p_i},\text{9}& {\displaystyle }{D}_{\mathit{\gamma }}=\prod _{i=\mathrm{1}}^A(p^a{)}^{-p^a}.\end{array}$$ Finally, the mixing state metric $\mathit{\chi }$ is defined as the affine ratio of $D_{\mathit{\alpha }}$ and $D_{\mathit{\gamma }}$: 10$$\mathit{\chi }={\displaystyle \frac{D_{\mathit{\alpha }}-\mathrm{1}}{D_{\mathit{\gamma }}-\mathrm{1}}}.$$ The values of $\mathit{\chi }$ vary between 0 % and 100 %. Take the three particle populations in Fig. S1 as an example. All three populations have the same bulk species mass concentration. Thus, they have the same bulk effective species diversity $D_{\mathit{\gamma }}$. However, the species are distributed differently within the populations. When the particles are externally mixed, each particle only contains one species and $D_{\mathit{\alpha }}$ is 1, which results in $\mathit{\chi }=\mathrm{0}\mathit{\%}$. When all particles have the same species mass fractions, $D_{\mathit{\alpha }}$ equals to $D_{\mathit{\gamma }}$, and we obtain $\mathit{\chi }$ of 100 %. The population is fully internally mixed. For the intermediately mixed population, $\mathit{\chi }$ ranges between 0 % and 100 %. In many applications, $\mathit{\chi }$ is used to describe the mixing state of chemical species, and we can therefore also refer to it as the chemical abundance metric ${\mathit{\chi }}_{\mathrm{chem}}$ .

For this work, our focus is the optical properties of the particles. Differing from ${\mathit{\chi }}_{\mathrm{chem}}$, we grouped the aerosol species by absorbing and non-absorbing species, i.e., BC and non-BC species, and defined a new index accordingly, ${\mathit{\chi }}_{\mathrm{BC}}$. It still ranges between 0 % and 100 % and signifies the degree to which BC and non-BC species are mixed. Since we only consider two (surrogate) aerosol species, the maximum value of $D_i$, $D_{\mathit{\alpha }}$ and $D_{\mathit{\gamma }}$ is 2. The same metric was chosen by to characterize the mixing state of BC-containing aerosol in Beijing and by to understand the role of mixing state in aerosol light absorption enhancement. For the remainder of the paper, unless otherwise noted, we will refer to the ${\mathit{\chi }}_{\mathrm{BC}}$ simply as $\mathit{\chi }$.

Figure  shows the range of bulk chemical species concentrations, the mixing state metric and optical properties within the selected scenario library. The simulated aerosol bulk species mass concentration in the library covered a wide range of urban conditions (Fig. a), and the values were comparable to the measurements in different locations . Both mixing state metrics ${\mathit{\chi }}_{\mathrm{chem}}$ and ${\mathit{\chi }}_{\mathrm{BC}}$ were larger than 30 %, with a median value of 85 %. The fact that mixing state metric values smaller than 30 % did not occur in our scenario library is consistent with the notion that aerosol species rarely exist in a completely external mixture but rather form some degree of internal mixtures already at the time of emission. Additionally, in urban environments, particles age quickly, forming internal mixtures with secondary species .

Our range of ${\mathit{\chi }}_{\mathrm{BC}}$ values encompasses the range observed in field measurements at Taizhou, China, where ${\mathit{\chi }}_{\mathrm{BC}}$ ranged between 68 % and 79 % for a period in May–June , and at Beijing, China, where ${\mathit{\chi }}_{\mathrm{BC}}$ ranged between 55 % and 70 % in winter and between 60 % and 75 % in summer . The range of ${\mathit{\chi }}_{\mathrm{chem}}$ is also consistent with the daily range of 37 %–73 % of particle samples collected in Paris, France .

Figure c shows that the single scattering albedo (SSA) was larger than 0.4 for all populations, with a median value of 0.88. While SSA values lower than 0.5 are considered extremely low (4 %), most populations (72 %) had a SSA larger than 0.85, which is consistent with fine-mode SSA observations from AERONET . The distribution of simulated total number concentration (Fig. d) is consistent with the observed number concentration of particles in the accumulation-mode size range . Note that the simulations presented here do not include the process of new particle formation. As a result, the simulated particle populations are more representative of accumulation-mode particles in a range of different environments.

Figure 3

Distribution of (a) bulk species concentration, (b) mixing state metric ${\mathit{\chi }}_{\mathrm{BC}}$ and ${\mathit{\chi }}_{\mathrm{chem}}$, (c) SSA and (d) total number concentration in the scenario library. Error bars in (a) are for $\pm$1 interquartile range (IQR), and numbers are the species median concentration in micrograms per cubic meter ($\mathrm{\mu }\mathrm{g}{\mathrm{m}}^{-\mathrm{3}}$).

[Figure omitted. See PDF]

3 Errors in aerosol absorptivity and scattering for dry particles

This section describes how we quantified the error introduced by composition averaging assumptions and how this error depends on mixing state. Similar to the approach used by , we stratified the populations by the black carbon mixing state metric $\mathit{\chi }$. To isolate the impacts of mixing state (in the sense of how the chemical species except for aerosol water are distributed across the population) from the impacts of water uptake, we first analyzed the results for the dry population scenarios $P_{\mathrm{1}}$. Particles were partially or fully deliquescent in scenarios $P_{\mathrm{2}}$ (RH $=\mathrm{50}$ %) and $P_{\mathrm{3}}$ (RH $=\mathrm{90}$ %). These populations will be further analyzed in Sect.  to quantify the water uptake effects on aerosol optical properties resulting from internally mixing hygroscopic and hydrophobic species.

The errors in aerosol optical properties due to the internal mixture assumption were defined by comparing the values of reference and composition-averaged populations. The relative error $\mathit{ϵ}$ for the aerosol populations was calculated as

11$$\mathit{ϵ}(v,\mathit{\chi })={\displaystyle \frac{v^{\prime }\left(\mathit{\chi }\right)-v\left(\mathit{\chi }\right)}{v\left(\mathit{\chi }\right)}},$$ where $v$ stands for ${\mathit{\beta }}_{\mathrm{abs}}$, ${\mathit{\beta }}_{\mathrm{scat}}$, $E_{\mathrm{abs}}$ or single scattering albedo in the reference library, and $v^{\prime }$ is for the same parameters in the sensitivity library. These parameters are stratified by the mixing state metric $\mathit{\chi }$.

3.1 Errors in aerosol absorptivity due to composition averaging

Absorption was overestimated universally after composition averaging, and, as expected, the error was higher for more externally mixed populations (low $\mathit{\chi }$ values), with $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ reaching up to $+\mathrm{70}\mathit{\%}$ for $\mathit{\chi }$ of 30 % (Fig. a). Each dot in Fig. a represents a particle population from the scenario library. As shown in the box plot inset, the mean overestimation was 18 %, and the maximum reached over 80 %. The figure further contains information of BC bulk mass concentration and relative average BC core size changes, which are the two main factors in determining absorptivity , as represented by marker size and color, respectively. The relative average BC core size change for a population is defined as

12$$\mathrm{\Delta }{D}^{\mathrm{core}}={\displaystyle \frac{{\sum }_{i=\mathrm{1}}^N{n}_i{D_i^{\mathrm{core}}}^{\prime }-{\sum }_{i=\mathrm{1}}^N{n}_i{D}_i^{\mathrm{core}}}{{\sum }_{i=\mathrm{1}}^N{n}_i{D}_i^{\mathrm{core}}}},$$ where $i$ is the particle index, $n_i$ and $D_i^{\mathrm{core}}$ are the associated number concentration and BC core diameter in the reference scenario, and ${D_i^{\mathrm{core}}}^{\prime }$ is the BC core diameter in the sensitivity scenario. The number concentration $n_i$ is always greater than 0, and if there is no core for particle $i$, $D_i^{\mathrm{core}}$ is 0. It is interesting to note that $\mathrm{\Delta }{D}^{\mathrm{core}}$ is always positive; that is, the average core diameter after composition averaging is larger than the average core diameter before composition averaging. This is a result of particle mass being a convex function of particle diameter (assuming spherical particles). Calculating the new core diameters after composition averaging will therefore always lead to larger core diameters on average than averaging the core diameters before composition averaging, as shown in Fig. S2.

The decreasing error with increasing $\mathit{\chi }$ can be explained by the magnitude of $\mathrm{\Delta }{D}^{\mathrm{core}}$. Evidently, composition averaging caused larger changes of BC core sizes when the populations were more externally mixed. For example, for $\mathit{\chi }=\mathrm{30}$ %, the change in core sizes was as large than $+\mathrm{25}$ %, while for $\mathit{\chi }=\mathrm{95}$ %, the change in core sizes was less than 5 %. We also noticed a range of errors for populations with $\mathit{\chi }$ between 60 % and 70 %, i.e., partially internally mixed populations. In fact, the highest overestimation of 82 % was reached at $\mathit{\chi }=\mathrm{63}$ %. As indicated by the circle size, these populations contained very little BC (0.01 $\mathrm{\mu }\mathrm{g}{\mathrm{m}}^{-\mathrm{3}}$), and even small changes in core sizes can lead to large relative errors in the volume absorption coefficient.

Figure 4

Relative error in absorption coefficients (a) $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ and (b) $\mathit{ϵ}\left(E_{\mathrm{abs}}\right)$ after composition averaging for dry particles. Each marker represents an aerosol population. The color denotes the change of BC diameter due to composition averaging, and the marker size represents BC bulk mass in the population. The box plot inset shows the distribution of the error. The red line shows the median, and the edges of the dashed lines are the minimum and maximum values. Red numbers are for the minimum, first quartile, median, third quartile and maximum values.

[Figure omitted. See PDF]

Given the constraint that composition averaging preserves the particle number concentration and sizes, it follows that, for some particles, this operation increases the sizes of BC cores (while at the same time decreasing the coating thickness), whereas for other particles it decreases the BC cores sizes (while increasing the coating thickness). It is therefore not immediately clear that composition averaging consistently causes overestimation of aerosol absorption coefficients.

At a per-particle scale, for particles of the same diameter, ${\mathit{\sigma }}_{\mathrm{abs}}$ increases with increasing BC core, even though the coating thickness (and hence the absorption enhancement) decreases (Fig. S3). However, $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ is determined by the entire population. The internal mixture in each size bin is reached by moving species from a group of particles to another group of particles. As the BC mass fraction distribution in Fig.  shows, there are two major groups of particles in the population: group 1 are particles with higher BC mass fraction, and group 2 are particles with lower BC mass. Particles in group 1 experience decreased absorbing ability because they are losing BC and vice versa for particles in group 2.

To further illustrate the effects at the population level, we show the effects of composition averaging on the volume absorption coefficient for a simplified case of five monodisperse populations of different sizes, starting out with completely externally mixed populations consisting of BC and ammonium bisulfate (Fig. ). Absorption coefficients are normalized by the absorption coefficient for $f_{\mathrm{BC}}=\mathrm{1}$ (pure BC). The black line shows the normalized volume absorption coefficient for populations when all particles are externally mixed for bulk BC mass fractions $f_{\mathrm{BC}}$ varying between 0 % and 100 %. For external mixtures, absorption increases linearly with increasing BC mass fraction (black line). The linear relationship applies for all five externally mixed populations with different diameters, so we can only see one black line in the figure.

Figure 5

Normalized absorption coefficient as a function of BC mass fraction for five monodisperse populations with different sizes. The coating species is ammonium bisulfate with refractive index 1.47. Absorption coefficients are normalized by ${\mathit{\beta }}_{\mathrm{abs}}$ of the population with $f_{\mathrm{BC}}=\mathrm{1}$ (pure BC). The black line is for BC in an external mixture. Colored lines are for BC in an internal mixture of different sizes. The table on the right sketches three 300 $\mathrm{nm}$ internal and external populations with BC mass fraction of 0 %, 50 % and 100 %. Black is for BC and yellow for coating species.

[Figure omitted. See PDF]

The colored lines represent the internally mixed monodisperse populations (i.e., after composition averaging) for different diameters. These populations all have higher absorption coefficients compared to the corresponding externally mixed populations. The effect is more pronounced for larger particles and intermediate BC mass fractions because the maximum ${\mathrm{\Delta }}_{\mathrm{core}}$ is reached. As the table (Fig. ) shows, for a 300 $\mathrm{nm}$ population, the normalized absorption is 0.76 when the particles are internally mixed, higher than an external mixed population (0.5). Although this example is an idealized case since our populations lie between external and internal mixtures before composition averaging and are polydisperse, this illustrates that assuming internal mixture will lead to absorption overestimation.

The coating redistribution after composition averaging also changes the absorption enhancement. As shown in Fig. S4, the median $E_{\mathrm{abs}}$ is 1.88 for the reference populations with BC mass fraction less than 10 %, while it is 1.98 for the corresponding populations of the sensitivity library. The absorption enhancement decreases as the bulk BC mass fraction decreases. These values are within the range of previous studies . Similar to the error in volume absorption coefficient $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$, the errors are larger for the populations for lower mixing state metric (Fig. b). Lastly, since composition averaging conserves the bulk species mass concentrations, the denominator in Eq. () (total BC mass concentration) remains unchanged, and the errors for ${\mathrm{MAC}}_{\mathrm{BC}}$ are the same as for ${\mathit{\beta }}_{\mathrm{abs}}$.

3.2 Error in aerosol scattering due to composition averaging

Considering the volume scattering coefficient, composition averaging resulted in a negative relative error (Fig. a). Similar to what we found for $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$, the magnitudes of $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ decreased with increasing $\mathit{\chi }$ but were overall smaller, with the largest underestimation of $-\mathrm{32}\mathit{\%}$ for a population with $\mathit{\chi }=\mathrm{40}\mathit{\%}$ and a median of $-\mathrm{1.2}\mathit{\%}$.

Figure 6

(a) Same as Fig.  but for $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$. The color is for refractive index relative change, and the marker size represents BC bulk mass in the population. The red box is the population analyzed in (b). (b) Size-resolved scattering coefficients for the reference and sensitivity (composition-averaged scenario library). The population is from scenario 77 at $t=\mathrm{2}$ h.

[Figure omitted. See PDF]

Two factors affect the particle scattering ability by composition averaging, the change of the BC core size (and the corresponding change in coating thickness) and the change in the refractive index of the coating. As Fig.  shows, adding a BC core decreases the scattering ability for particles with diameters less than 1200 $\mathrm{nm}$, which is the typical size range considered in our study. This explains the larger scattering underestimation with higher BC mass concentration in Fig. a.

Figure 7

Relation between scattering cross section, refractive index and diameter for a wavelength of 550 nm. Blue lines are for non-absorbing particles, and symbols indicate the refractive index. The red line is for absorbing particles containing a BC core of $\mathrm{0.2}{D}_{\mathrm{dry}}$ (constant BC volume fraction across the size range).

[Figure omitted. See PDF]

To further explore the effects of coating volume changes, Fig. b shows the size-resolved scattering coefficients before and after composition averaging for the aerosol populations from scenario 77 at $t=\mathrm{2}$ h, which produced the largest scattering coefficients underestimation ($-\mathrm{32}\mathit{\%}$). There is a significant decrease of ${\mathit{\sigma }}_{\mathrm{scat}}$ in the size range of 400–800 $\mathrm{nm}$ in the sensitivity populations, and the core mass ratio increment in bin 4 is responsible for this decrease (Fig. S5).

The blue lines in Fig.  show the scattering cross sections for two different real refractive indices. For particles with diameters between 800 and 1200 $\mathrm{nm}$, a lower refractive index leads to a larger scattering cross section, although the difference is smaller than the change caused by adding a BC core. Similar to the BC core size change $\mathrm{\Delta }{D}^{\mathrm{core}}$ in Fig. , we defined a volume-weighted refractive index change, $\mathrm{\Delta }{m}^{\mathrm{real}}$, to help understand the changes in scattering. The index change is defined as 13$$\mathrm{\Delta }{m}^{\mathrm{real}}={\displaystyle \frac{{\sum }_{i=\mathrm{1}}^N{V}_i{m^{\mathrm{real}}}^{\prime }-{\sum }_{i=\mathrm{1}}^N{V}_i{m}^{\mathrm{real}}}{{\sum }_{i=\mathrm{1}}^N{V}_i{m}^{\mathrm{real}}}},$$ where $i$ is the particle index, $V_i$ is the particle volume, $m^{\mathrm{real}}$ is the real part of the coating refractive index of the particles in the reference library and ${m^{\mathrm{real}}}^{\prime }$ is for particles in the sensitivity library. Here we applied the total particle volume $V_i$ in the equation to focus on the relation between the changes in scattering and changes in the refractive index. As shown in Fig. , aerosol populations with small errors in scattering tend to be associated with small $\mathrm{\Delta }{m}^{\mathrm{real}}$. For more externally mixed populations (with lower $\mathit{\chi }$), $\mathrm{\Delta }{m}^{\mathrm{real}}$ tended to be larger.

For the effects of composition averaging for particle scattering, we conclude that at a given value of $\mathit{\chi }$, the magnitude of $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ was determined by the change in core/coating volumes and by changes in the coating refractive index. The increase of BC core sizes after composition averaging is the major factor for the decrease of the scattering coefficients. Populations with large underestimation are those with higher BC mass concentrations and large refractive index changes. It is worth emphasizing that we did not consider the absorption of organic carbon that might be present in the coating .

4 The effects of water uptake on aerosol optical properties

The analysis so far has been based on dry aerosol populations. In this section we investigate the impact of water uptake on the errors in absorption and scattering by considering RH values of 50 % and 90 %. As a reminder, we performed composition averaging on the dry population first and then calculated water uptake based on the averaged composition for RH $=\mathrm{50}\mathit{\%}$ and RH $=\mathrm{90}\mathit{\%}$, respectively.

Figure 8

Box plot of (a) scattering relative error $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ and (b) absorption relative error $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ at three RH levels (0 %, 50 % and 90 %). Dots are the populations with values outside Q3 $+$ 1.5 times the IQR.

[Figure omitted. See PDF]

Considering all populations, the range of relative errors in ${\mathit{\beta }}_{\mathrm{scat}}$ decreased with increasing RH, with the median error over all populations decreasing from $-\mathrm{1.2}\mathit{\%}$ (RH $=\mathrm{0}\mathit{\%}$) to $-\mathrm{1.0}\mathit{\%}$ (RH $=\mathrm{50}\mathit{\%}$) and $-\mathrm{0.2}\mathit{\%}$ (RH $=\mathrm{90}\mathit{\%}$) (Fig. a). In contrast, the range of relative errors in ${\mathit{\beta }}_{\mathrm{abs}}$ remained approximately the same (Fig. b), with a median of approximately 13 %.

The different response of $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ and $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ after the populations became humidified was due to the scattering coefficients increasing strongly at higher relative humidities (Fig. S6a). The enhancement ratio, defined by the ${\mathit{\beta }}_{\mathrm{scat}}$ values for the higher RH cases and the dry case, had a median of 1.33 at RH $=\mathrm{50}$ % and 3.35 at RH $=\mathrm{90}$ % in our scenario populations. These values are in accordance with previous studies .

As for the absorption coefficients in the humidified environments, the differences between reference and sensitivity cases remained almost the same for both ${\mathit{\beta }}_{\mathrm{abs}}$ and ${\mathrm{MAC}}_{\mathrm{BC}}$ (Fig. S6b and c), indicating that the errors in absorptivity introduced by composition averaging were not sensitive to RH.

5 Errors in single scattering albedo and implications for directive radiative forcing

The changes of scattering and absorption coefficients lead to changes in SSA, which is an important quantity that determines radiative forcing. With the definition of SSA, we can calculate the absolute error $\mathrm{\Delta }$SSA as

14$$\begin{array}{rl}\mathrm{\Delta }\mathrm{SSA}& ={\displaystyle \frac{{\mathit{\beta }}_{\mathrm{scat}}^{\prime }}{{\mathit{\beta }}_{\mathrm{scat}}^{\prime }+{\mathit{\beta }}_{\mathrm{abs}}^{\prime }}}-{\displaystyle \frac{{\mathit{\beta }}_{\mathrm{scat}}}{{\mathit{\beta }}_{\mathrm{scat}}+{\mathit{\beta }}_{\mathrm{abs}}}}\\ & ={\displaystyle \frac{{\mathit{\beta }}_{\mathrm{scat}}^{\prime }{\mathit{\beta }}_{\mathrm{abs}}-{\mathit{\beta }}_{\mathrm{scat}}{\mathit{\beta }}_{\mathrm{abs}}^{\prime }}{({\mathit{\beta }}_{\mathrm{scat}}^{\prime }+{\mathit{\beta }}_{\mathrm{abs}}^{\prime })({\mathit{\beta }}_{\mathrm{scat}}+{\mathit{\beta }}_{\mathrm{abs}})}},\end{array}$$ where ${\mathit{\beta }}_{\mathrm{scat}}^{\prime }$ and ${\mathit{\beta }}_{\mathrm{abs}}^{\prime }$ refer to the scattering and absorption coefficients after composition averaging. Based on the previous analysis, we know that ${\mathit{\beta }}_{\mathrm{scat}}^{\prime }$ tends to be lower than ${\mathit{\beta }}_{\mathrm{scat}}$ and ${\mathit{\beta }}_{\mathrm{abs}}^{\prime }$ greater than ${\mathit{\beta }}_{\mathrm{abs}}$. Combining these changes with Eq. (), these variations will result in negative values for $\mathrm{\Delta }$SSA and the relative error $\mathit{ϵ}\left(\mathrm{SSA}\right)$, which is confirmed by Fig. .

Figure 9

Same as Fig. a but for $\mathit{ϵ}\left(\mathrm{SSA}\right)$. Red numbers in the inset box plot are for the minimum, median and maximum values.

[Figure omitted. See PDF]

Figure  shows that $\mathit{ϵ}\left(\mathrm{SSA}\right)$ was negative for all the dry aerosol populations, with a median value of $-\mathrm{0.9}\mathit{\%}$ and a largest value of $-\mathrm{22.3}\mathit{\%}$. The dependence of $\mathit{ϵ}\left(\mathrm{SSA}\right)$ on the mixing state metric $\mathit{\chi }$ shows a similar pattern as for the volume scattering coefficient $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$. The errors decreased with increasing $\mathit{\chi }$, indicating the SSA underestimation was smaller for more internally mixed populations. For the populations with the same mixing state metric $\mathit{\chi }$, errors were higher for the populations with more BC mass concentrations. Aerosol populations with higher SSA errors were also associated with higher refractive index changes.

In order to further connect $\mathit{ϵ}\left(\mathrm{SSA}\right)$ with $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ and $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ and investigate the effects of RH, we sorted the populations by $\mathit{ϵ}({\mathit{\beta }}_{\mathrm{scat})}$ and $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ ranges and calculated the averaged $\mathit{ϵ}\left(\mathrm{SSA}\right)$ for each $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)-\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ bin for the three RH levels, as shown in Fig. . For all three RH levels, $\mathit{ϵ}\left(\mathrm{SSA}\right)$ was negative, meaning that composition averaging causes an underestimation of SSA. The largest $\mathit{ϵ}\left(\mathrm{SSA}\right)$ ($-\mathrm{22.3}$ %) occurred for the largest underestimation in $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ in the RH $=\mathrm{0}$ % environment. Populations with $\mathit{ϵ}\left(\mathrm{SSA}\right)$ lower than $-\mathrm{10}$ % were related to populations with large negative magnitudes of $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$. Relative errors in SSA decreased in a more humidified environment, accompanied by decreasing errors in scattering coefficients. The median underestimation of SSA decreased from $\mathrm{0.9}$ % (RH $=\mathrm{0}$ %) to $\mathrm{0.7}$ % (RH $=\mathrm{50}$ %) and $\mathrm{0.4}$ % (RH $=\mathrm{90}$ %).

Figure 10

Relation between errors in SSA, scattering and absorption coefficients. Color represents the averaged $\mathit{ϵ}$(SSA) in the corresponding $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ and $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ histograms.

[Figure omitted. See PDF]

The underestimation of SSA can have significant impacts in calculating direct radiative forcing. evaluated the response of directive radiative forcing to changes of several quantities, including aerosol optical depth and single scattering albedo. They found that the total uncertainties in directive radiative forcing ranged from 0.2 to 3.1 $\mathrm{W}{\mathrm{m}}^{-\mathrm{2}}$, and SSA introduced the largest uncertainties. Through perturbation analysis, also found the SSA to be the dominant factor for direct radiative forcing uncertainties. They perturbed SSA by $\pm \mathrm{3}\mathit{\%}$ over land, which resulted in uncertainties in direct aerosol radiative forcing between $-\mathrm{0.59}$ and $+\mathrm{0.69}$ $\mathrm{W}{\mathrm{m}}^{-\mathrm{2}}$. The median SSA errors for our simulations were on the order of the perturbations imposed in the study by , and we therefore conclude that mixing state effects can have impacts on radiative forcing similar to the ones reported in . Furthermore, the spatial and temporal variations of relative humidity imply that the errors in optical properties for a population with a given mixing state may vary depending on location, season and time of day. Radiative transfer calculations would be required for a more in-depth analysis of radiative forcing impacts.

6 Conclusions

Simplified representation of aerosol mixing state used in current regional or global models may introduce errors in simulating aerosol optical properties, thus leading to uncertainties in calculating directive radiative forcing. In this study, the errors introduced by internal mixture assumptions used in sectional aerosol models were systematically quantified. We created a reference scenario library with 1800 aerosol populations by performing particle-resolved aerosol model simulations with PartMC-MOSAIC. We constructed a sensitivity library where particles were internally mixed in a prescribed set of size bins by applying composition averaging. This operation has the properties of conserving number concentration and particle sizes, and hence differences in any quantity can be solely attributed to mixing state impacts. Aerosol populations from the reference and sensitivity library were then exposed to three different RH levels to understand the relative role of chemical species and water redistribution introduced by the internal mixture assumption.

The internal mixture assumption generally led to an overestimation of the volume absorption coefficients and an underestimation of the volume scattering coefficients. The relative errors for $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ and $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ reached up to 70 % and $-\mathrm{32}\mathit{\%}$, respectively. The relative errors generally increased for more externally mixed populations, although at a given value of $\mathit{\chi }$ a range of errors could be found, especially for the error in the scattering coefficient. The range of error in the absorption coefficient can be explained by the magnitude of BC core size changes that are induced by composition averaging. The finding of overestimation of BC absorption due to simplified mixing state representation was consistent with many other studies, including the works by , , and . The error in the scattering coefficient can be explained by the magnitude of the changes in the refractive index of the coating that are induced by the composition averaging.

For the cases with RH of 50 % and 90 %, the bulk aerosol water content was almost identical for the aerosol populations in reference and sensitivity libraries. The relative error in the volume absorption coefficient $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{abs}}\right)$ displayed a similar pattern for RH of 50 % and 90 % compared to the dry environment. The relative error in the volume scattering coefficient $\mathit{ϵ}\left({\mathit{\beta }}_{\mathrm{scat}}\right)$ decreased for higher relative humidities because of the enhanced scattering cross section through hygroscopic growth.

The absorption overestimation and scattering underestimation resulted in an consistent underestimation of SSA, with median errors of $-\mathrm{0.9}$ % (RH $=\mathrm{0}$ %), $-\mathrm{0.7}\mathit{\%}$ (RH $=\mathrm{50}$ %) and $-\mathrm{0.4}\mathit{\%}$ (RH $=\mathrm{90}$ %). Populations with the largest underestimation of SSA ($-\mathrm{22.3}$ %) were associated with populations with the largest underestimation in scattering.

It is worth emphasizing that we used Mie theory with a core–shell configuration to calculate optical properties assuming spherical particle shapes. Our results are therefore most representative of BC-containing populations where the BC core is collapsed rather than a fractal aggregate . More accurate methods, such as discrete dipole approximation (DDA) and dynamic shape factor should be used to represent these more irregular particle shapes .

Code and data availability

The simulation data and codes are available at 10.13012/B2IDB-8157303_V1 (). PartMC v2.6.0 is archived at 10.5281/zenodo.5644422 ().

The supplement related to this article is available online at: https://doi.org/10.5194/acp-22-9265-2022-supplement.

Author contributions

YY, ZZ and NR designed the particle-resolved scenario libraries. JHC developed codes for calculating per-particle optical properties, and JC contributed to interpret results. YY and NR performed the analysis and prepared the manuscript, with edits from co-authors.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

Joseph Ching is an International Research Fellow of Japan Society for the Promotion of Science (JSPS) and acknowledges the financial support from the JSPS Postdoctoral Fellowships for Research in Japan (Standard). Zhonghua Zheng is funded by the NCAR Advanced Study Program Postdoctoral Fellowship. This material is based upon work supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation under cooperative agreement no. 1755088.

Financial support

This research has been supported by the National Science Foundation (grant no. 1254428); the U.S. Department of Energy (grant no. DESC0022130); the Japan Society for the Promotion of Science (grant no. JP19F19402); the Research Institute for Humanity and Nature (grant no. 14200133); the Ministry of Education, Culture, Sports, Science and Technology (grant no. JPMXD1420318865); and the Meteorological Research Institute (MRI), Japan (grant no. Fundamental Technology Research, M5 and P5).

Review statement

This paper was edited by Dantong Liu and reviewed by two anonymous referees.