1 Introduction
The direct and indirect climate effects of atmospheric aerosols greatly depend on the particles' spatial distribution in the atmosphere and their climate-relevant properties, including their hygroscopicity, optical properties, and their ability to act as cloud condensation nuclei (CCN) and ice nuclei . These properties, in turn, are closely related to the aerosol mixing state . Aerosol mixing state refers to the way in which different aerosol chemical species are distributed among and within the aerosol particles . As shown in many observational field studies, atmospheric aerosols have complex mixing states , ranging between the two extremes of an “internal mixture”, where the composition of all particles within the population is identical (and equal to the bulk composition of the aerosol), and an “external mixture”, where each particle in a population consists of only a single species (which may be different for each particle).
This poses a unique challenge for the modeling of aerosols in Earth system models, which, for the sake of computational efficiency, represent aerosols by simplifying the true aerosol mixing state using various mixing-state-related assumptions. For example, bulk aerosol models predict the abundance of individual aerosol chemical species by tracking the species' mass concentrations, inherently treating the aerosol as external mixtures of, e.g., sulfate, black carbon, organic carbon, sea salt, and dust . Univariate sectional models are able to represent size-resolved composition but cannot resolve the diversity of the aerosol within a certain size range. For modal models, the ability to resolve mixing state depends on the definition and the placement of the modes. Different approaches for modal models have been developed, ranging from a small number of internally mixed, non-overlapping modes (e.g., three modes in MAM3, ; or CMAQv5.2, ) to a larger number of modes that may overlap in a given size range and separate out different aerosol mixtures (e.g., nine modes in MADE3, ; or 16 modes in MATRIX, ). For these multi-modal models, the aerosol processes of gas–aerosol partitioning and coagulation make it necessary to define rules for how the modes interact . Condensation of secondary aerosol on a mode reserved for a pure species (e.g., black carbon or dust) requires moving mass over to a mixed mode when a critical mass fraction of secondary aerosol is exceeded. Transfer terms due to coagulation of particles in different modes can be calculated analytically , and rules need to be defined regarding the destination mode after coagulation. Generally, the transfer of aerosol mass from smaller modes to larger modes during growth can lead to inaccuracies. The removal of particles due to scavenging by cloud activation is another issue that is difficult to reconcile. Hence, the choice of the number of modes, their compositions, and the criteria for transfer between modes are user-defined, which introduces structural uncertainty in aerosol simulations that still needs to be quantified.
Given that modal models are to some extent mixing-state-aware, the following question arises: how well do modal models represent mixing state? Due to the scarcity of relevant observational data, we are not yet at the point where we can comprehensively validate model output of aerosol mixing state as is done for other aerosol-related quantities, such as bulk mass concentrations or aerosol optical depth. However, higher-detail models can serve as benchmarks to perform a verification of simulated aerosol mixing state. This paper aims to verify the global distribution of aerosol mixing state represented by a modal model by using benchmark simulations from the particle-resolved stochastic aerosol model PartMC-MOSAIC . Our usage of the term “aerosol representation” in this paper encompasses the representation of processes that go along with the aerosol representation itself, since the two are in practice tightly coupled.
We used the aerosol mixing state index as a metric to quantify aerosol mixing state. The mixing state index can be interpreted as a label for particle populations to rigorously characterize where the population lies on the spectrum from external ( %) to internal ( %) mixture. This concept has been successfully applied to observational data and for error quantification studies . Particularly relevant for this work is the study by , which showed that assuming an internal mixture when the aerosol is actually not completely internally mixed can result in errors of up to 150 % in CCN predictions.
PartMC-MOSAIC tracks the composition of individual particles and therefore resolves aerosol mixing state explicitly . However, this modeling approach is computationally very expensive and therefore not practical for large-scale simulations of several months or years of simulation time. To estimate the global spatial distribution of mixing state, we recently developed a machine-learned (ML) model based on high-detail particle-resolved simulations that uses inputs that are known from global model simulations to predict . In this paper, we use this ML model to predict the spatial distribution of the mixing index and then compare the results with values that are derived from the Community Earth System Model version 2
This paper is organized as follows. In Sect. we introduce the setup of the Earth system model simulations. The definition of mixing state indices and the derivation of aerosol mixing state indices for modal models are given in Sect. . Section briefly describes the ML model generated with machine learning and particle-resolved modeling for estimating the benchmark aerosol mixing state indices. Section focuses on the comparison of mixing state indices from the particle-resolved and modal models, and Sect. summarizes our findings.
2 Global model simulations
Here we employed CESM2 to provide the global model simulation data. Specifically, we used the component set FHIST to set up the global simulations with aerosols. This component set represents a typical historical simulation in the Community Atmospheric Model
MAM4 is the default aerosol module of this component set, which represents the aerosol size distribution with four lognormal modes (Aitken, accumulation, coarse, and primary carbon modes; ). MAM4 tracks six aerosol species, and these are distributed over the four modes as follows. The Aitken mode consists of dust, sulfate, secondary organic aerosol (SOA), and sea salt. The accumulation mode includes sulfate, SOA, sea salt, primary organic matter (POM), black carbon (BC), and dust. The coarse mode contains sulfate, dust, and sea salt. The primary carbon mode contains only BC and POM, which are supplied by primary aerosol emissions.
The choice of modes in MAM4 is motivated by the desire to treat the microphysical aging of the primary carbonaceous aerosols in the atmosphere similar to other modal models used in regional or global models . In MAM4, mass and number concentrations of BC and POM in the primary carbon mode are transferred to the accumulation mode by the processes of intermodal coagulation and condensation of SOA and sulfuric acid onto the primary carbon mode. The accumulation mode then represents aged BC and POM, as these species are internally mixed with other aerosol species. The MAM4 treatment of aging is critical for improving the long-range transport of carbonaceous aerosols to remote regions such as the polar region, which suffered from a low bias in a prior version of the model when only three internally mixed modes were used .
We ran the model for the year 2011 with 6 years (2005–2010) of spinup. The simulation was conducted at a resolution of 0.9 latitude by 1.25 longitude along with emission inventories from CMIP6 emissions . We stored the instantaneous outputs every 3 h during the simulation, which yields 2920 timestamps for each surface-layer grid cell for the entire year of simulation time. The surface layer was chosen to be in line with the PartMC-MOSAIC model scenarios that were used as training data for the ML models of mixing state indices (see Sect. ) and which were designed to represent conditions in the planetary boundary layer.
3 Aerosol mixing state indices: definition and calculation
3.1 Particle-based aerosol mixing state index
The mixing state index quantifies where an aerosol population lies on the continuum from external to internal mixing – that is, how spread out the chemical species are over an aerosol population. We will focus here on the mixing state of submicron aerosols (PM) due to their relevance for light scattering and absorption and their contribution to CCN formation .
To summarize, the mixing state index is given by the affine ratio of the average particle species diversity, , and bulk population species diversity, , as
1
The diversities and are calculated as follows. First, the per-particle mixing entropies are determined for each particle by 2 Here, is the number of distinct aerosol species and is the mass fraction of species in particle . These values are then averaged (mass-weighted) over the entire population to obtain the average particle species diversity by where is the total number of particles in the population and is the mass fraction of particle in the population. Finally, the bulk diversity is calculated as where is the bulk mass fraction of species in the population.
Note that the definition of “species” for calculating is based on application needs. It can be based on operationally defined chemical species , elemental composition , or species groups such as volatile and nonvolatile species or hygroscopic and non-hygroscopic species . Other possibilities include the propensity for aerosols to undergo heterogeneous reactions, quantified by the heterogeneous reaction rate coefficient for a specific reaction. In this paper we consider three different definitions of , which we explain in more detail in Sect. .
3.2 Mode-based aerosol mixing state indexThe framework laid out in Sect. can be easily generalized to a modal modeling framework (see Fig. ). The bulk mixing entropy, , and the bulk diversity, , can be calculated using the bulk mass fractions, , of species from the MAM4 simulation and Eqs. () and (). To calculate the average particle mixing entropy, , and the average particle species diversity, , we use where is the mass fraction of species in mode , is the mass fraction of mode in the population, and represents the per-mode mixing entropies. Finally, the mixing state index, , can be calculated using Eq. (). Note that Eqs. () and () are analogous to Eqs. () and (). A detailed derivation of these equations is provided in the Appendix .
Figure 1
Illustration of the mode-based calculation of the aerosol mixing state index. The coarse mode is removed because only modes dominated by submicron particles are used for calculations. Note that the Aitken mode mass fraction is very low compared to the other modes and the caption does not obscure any data.
[Figure omitted. See PDF]
In this study, we consider the mixing states of submicron aerosols including the Aitken, accumulation, and primary carbon modes, and we do not include the coarse mode because the coarse particles are above 1 m. Since the mixing entropies are mass-weighted (rather than number-weighted), the mixing state index is more representative of the modes with the larger particles, i.e., the accumulation and primary carbon modes.
3.3 Grouped surrogate speciesHere we compare and contrast the aerosol mixing state indices defined in three different ways, namely based on the mixing of optically absorbing and non-absorbing species (), based on the mixing of primary carbonaceous and non-primary carbonaceous species (), and based on the mixing of hygroscopic and non-hygroscopic species (). Table shows the definitions of these aerosol mixing state indices.
Table 1
Aerosol mixing state index definitions. Six aerosol species (bc: black carbon, dst: dust, ncl: sea salt, pom: primary organic matter, soa: secondary organic aerosol, so4: sulfate) are used in calculating the aerosol mixing state indices based on different species groupings. The mixing state indices , , and are based on two grouped surrogate species.
Aerosol mixing state index (symbol) | Grouped species | |
---|---|---|
Optical property () | (bc) | (pom, dst, ncl, soa, so4) |
Primary carbon () | (bc, pom) | (dst, ncl, soa, so4) |
Hygroscopicity () | (bc, pom, dst) | (ncl, soa, so4) |
For , we considered two surrogate species: black carbon (strongly absorbing, assigned a mass absorption coefficient in CESM2 at 533 nm and 0 % RH of 8.144 m g) and the five other aerosol species grouped together (less absorbing or non-absorbing, with mass absorption coefficients in CESM2 at 533 nm and 0 % RH of 0.1442, , , , and m g for POM, SOA, dust, sea salt, and sulfate, respectively). Thus, a lower value in refers to the case where the strongly absorbing species black carbon and the sum of the other species (termed “non-absorbing” here for convenience) are more externally mixed.
The index is motivated by the primary carbon treatment of MAM4, where the primary particulate organic matter and black carbon are assigned to a separate primary carbon mode . A lower value in refers to the situation where the primary carbonaceous species and all other species exist separately in different particles.
Similarly, was also calculated from two surrogate species. We combined black carbon, primary organic matter, and dust as one surrogate species, given their comparatively lower hygroscopicities (kappa values of , , and 0.068, respectively). Accordingly, NaCl (1.16), SOA (0.14), and sulfate (0.507) were grouped as the other surrogate species. Here, a lower value in represents the case where hygroscopic and non-hygroscopic species tend to be present in separate particles.
4 Machine-learned models of mixing state indicesAerosol mixing state indices can be calculated directly using particle-resolved modeling, but this comes with large computational costs. Alternatively, developed ML models, which integrate machine learning and particle-resolved aerosol simulations to estimate aerosol mixing state indices. To generate the training and testing data sets for developing such ML models, an ensemble of particle-resolved model scenarios was created using the particle-resolved model PartMC-MOSAIC . In brief, PartMC-MOSAIC simulates individual aerosol particles within a representative volume of air, including stochastic coagulation, particle-phase thermodynamics, gas- and particle-phase chemistries, and dynamic gas–particle mass transfer. Thus, the composition of the individual particles within a population evolves dynamically, and assumptions about mixing state are not necessary.
The strategy to generate the data was to vary the input parameters (45 in total) for the PartMC-MOSAIC model, including primary emissions of different aerosol types (e.g., carbonaceous aerosol and dust emissions, including contribution from Aitken mode, accumulation mode, and coarse mode size ranges), primary emissions of gas phase species (e.g., SO, NO, and various volatile organic compounds), and meteorological parameters (see Table 1 in , for more information). For instance, to vary the gas emissions, scaling factors were sampled from 0 % to 200 % for different gas species, based on the emission rates in . A Latin hypercube sampling approach was employed to sample the parameter space efficiently for the training and testing data sets. We note that new particle formation and growth was not simulated explicitly, but Aitken mode sulfate particles were introduced into the simulation by emission for a subset of scenarios as a proxy for having particles present that originate from new particle formation. While PartMC-MOSAIC includes the process of new particle formation , the reason for this simplification was that considerable uncertainty exists regarding the subsequent growth of the freshly nucleated particles , which poses a challenge for a highly detailed aerosol model such as PartMC-MOSAIC. Errors in representing this particle type adequately may result in underestimating the abundance of BC-free particles in some regions and thereby overestimating the degree of internal mixture. This would imply that the error in the MAM4 simulations is even larger than currently indicated. Other processes that are not explicitly included in generating training data are aerosol removal by nucleation-scavenging and other cloud processes. However, for the purpose of this study, the emphasis is on the aerosol state, i.e., having a sufficiently comprehensive set of aerosol populations that can serve as training data, not necessarily that all the processes are included.
The ML models were derived by the machine learning algorithm eXtreme Gradient Boosting
10 where is the mixing state index (, , or ) at location in the model layer nearest the surface at time , and denotes the function for calculating the corresponding mixing state index . The set names (aerosol), (gas), and (environmental) represent the predictors (features) used for predicting the mixing state index. The choice of features is determined by the overlap of variables that are present in both PartMC-MOSAIC and CESM2. Aerosol species include black carbon, mineral dust, sea salt, primary organic aerosol, secondary organic aerosol, and sulfate. Of note is that we used the bulk (not the per-mode) concentrations of submicron aerosol species as the features. The gas species include dimethyl sulfide, hydrogen peroxide, sulfuric acid, ozone, semi-volatile organic gas, and sulfur dioxide. The environmental variables are air temperature, relative humidity, and solar zenith angle. Table shows the performance of the ML models when predicting the mixing state indices. The mixing state calculation in this study was purely based on the above six aerosol species (excluding other aerosol species) for a fair comparison with the mode-based aerosol mixing state index, which resulted in slightly different performance of the ML model compared to . The average error of the ML model (using the hold-out testing samples) is about 5 % for and 8 % for and (measured by mean absolute error).
Table 2Predictive performance of the ML models using the testing data set. Metrics include the mean absolute error (MAE), root-mean-square error (RMSE), median absolute deviation (MAD), index of agreement (; ), Pearson correlation coefficient (PCC), and coefficient of determination ().
MAE | RMSE | MAD | PCC | |||
---|---|---|---|---|---|---|
XGBoost ML models | ||||||
0.048 | 0.072 | 0.030 | 0.974 | 0.953 | 0.906 | |
0.079 | 0.107 | 0.056 | 0.955 | 0.916 | 0.836 | |
0.082 | 0.112 | 0.057 | 0.955 | 0.916 | 0.835 |
We would like to emphasize that this ML modeling framework cannot compensate for any biases that the global model (here CESM2) might have in simulating the quantities that serve as the features. Instead, what we can expect from this approach is that it provides the most likely mixing state associated with the species concentrations that CESM2 simulates.
5 Results5.1 Quantitative comparison of mode-based and particle-based mixing state indices
Let and denote the mixing state indices computed by the ML model and by the MAM4 model for each grid cell at timestamp , respectively. The corresponding time-averaged values for a certain time interval and for each grid cell are and . Here we consider the full year as the time-averaging interval. An analysis of the seasonal variation of mixing state indices can be found in .
To compare the annual mean values, we calculated the mean difference () and the mean absolute difference () for each grid cell of the layer closest to the surface: where the subscript refers to the mixing state index (o, c, or h), and the total number of timestamps is . Since it only makes sense to quantify mixing state when at least two species are present in a given location, areas where the mass fraction of any one surrogate species was higher than 99 % for (due to the low mass fraction of black carbon) and 97.5 % for and were ignored for the calculation and appear as hatched areas in Fig. . We will first discuss the overall probability density functions of these quantities (Fig. ) and then their spatial distributions (Fig. ).
Figure 2
Probability density functions of annual averaged mixing state indices using the MAM4 model and ML model. The thin black lines refer to their mean values.
[Figure omitted. See PDF]
Figure 3
Global distribution of annually averaged mixing state indices (, , and ) using the ML model, MAM4 model, their mean difference (), and mean absolute difference (). Areas are hatched where the mass fraction of any one surrogate species was higher than 99 % for (due to the low mass fraction of black carbon) and 97.5 % for and .
[Figure omitted. See PDF]
Figure shows the probability density functions of the annual averaged mixing state indices computed by the ML model (), by MAM4 (), their average difference (), and their average absolute difference () for each surface-layer grid cell. The results show large discrepancies in mixing state indices between the ML model and the MAM4 model, without a clear relationship between them (see Fig. d–f).
The annual average of the mixing state index estimated by the ML model, , ranged between 55 % and 96 %, with a mean of 73 %. Calculated by the MAM4 model, varied spatially from 46 % to 99.76 %, with a higher mean of 86 %. The similar mean values of (14 %) and (18 %) were caused by higher values in compared to , which is confirmed below with Fig. . The averaged mixing state index ranged between 31 % and 84 % with a mean of 54 %, while had a wider range (from 9 % to 99.81 %) with a mean (of 58 %). Similarly, ranged from 21 % to 81 % with a mean of 58 %, while varied between 10 % and 99.85 % with a mean of 63 %. The large discrepancy between the mean difference (4.8 % for and 4.7 % for ) and mean absolute difference (30 % for and 38 % for ) indicates that the errors in and were symmetric (positive and negative) but large. The maximal errors in and between the two methods were up to 59 and 76 percentage points, respectively.
The implications of these discrepancies are more easily discussed with Fig. , which illustrates the global spatial distribution of annually averaged mixing state indices predicted by the ML model (first column), MAM4 (second column), their mean difference (third column), and their mean absolute difference (fourth column). The differences in mixing state indices between the ML model and MAM4 varied strongly across the globe.
High values of occurred in the continental regions (77 %) compared to oceans (69 %). Specifically, the ML model predicted high values for in central Africa (20 S–15 N, 12–30 E), the Arctic (66.5–90 N), and southern Asia (5–38 N, 60–90 E). These are also the regions with relatively larger mass fractions of black carbon ( %; see Fig. ). The mixing state index showed a higher degree of internal mixing over the globe (with a median of 90 %) compared to the ML model. The only exceptions were oceans in the Northern Hemisphere at the mid-latitudes (45–60 N, dominated by sea salt, sulfate, and secondary organic aerosol in the accumulation mode) and Antarctica (66.5–90 S, dominated by sea salt and sulfate in the accumulation mode as well as sulfate in Aitken mode), where was 75 %. Qualitatively, the MAM4 model captured the trend that areas with high black carbon concentration (defined here as concentrations above the 95 % percentile) tended to have higher values.
The ML model estimate suggested a rather homogeneous spatial distribution of the annually averaged mixing state, with values of approximately 50 %. Compared to , values were lower (primary carbonaceous aerosol more externally mixed) at high latitudes and higher at low and mid-latitudes (primary carbonaceous aerosol more internally mixed). Note that, while values were similar in the Arctic and Antarctic, the abundance of primary carbonaceous species was predicted to be higher in the Arctic compared to the Antarctic (see Fig. ).
The spatial distributions of were similar to . That is, the MAM4 model predicted that the hygroscopic species and non-hygroscopic species were more externally mixed at high latitudes and more internally mixed at low latitudes. In contrast, the spatial distribution of shows qualitative differences compared to in two aspects. First, was higher than at high latitudes, meaning that hygroscopic species and non-hygroscopic appeared more internally mixed than primary carbonaceous and non-carbonaceous species in this region. Second, areas over the North Atlantic Ocean (0–20 N, 20–45 W), southern Africa (5–32 S, 5–20 E), and Australia (10–30 S, 100–140 E) appeared rather externally mixed. These are areas where mineral dust is the dominant aerosol species (see Fig. ).
These two facts lead to the overall finding that exhibits the largest differences between the two methods. This applies especially to regions where mineral dust was the dominant aerosol species, which points to an important structural issue of the four-mode setup used in MAM4. While the ML model predicted a more external mixture in these regions (dust externally mixed from sea salt and other species), the MAM4 model could not represent this because the accumulation mode included all six aerosol species in an internal mixture. Figure illustrates the relationship of the mean absolute difference of and the mass fraction of dust for all model grid points. It confirms that grid points with large dust mass fractions were associated with larger mean absolute differences in . These results confirm the tradeoff discussed in : MAM3 (and MAM4 in ) intentionally combines dust and sea salt in the same mode to reduce the computational burden; however, this simplification does not always realistically reflect the aerosol mixing state in the ambient atmosphere.
It is interesting to note that the areas where sea salt is present, but not dust, are not associated with large errors, even though sea salt – just like mineral dust – is a primary aerosol type. The reason for this lies in our surrogate species definitions (Table ) for computing the mixing state index. Based on our mixing state definitions, sea salt, secondary organic aerosol, and sulfate are always grouped together. Therefore, none of the mixing state indices as defined here tell us how externally mixed sea salt is when it is considered as a single aerosol type.
Figure 4
Dependence of mean absolute difference of on dust mass fractions for all model grid points.
[Figure omitted. See PDF]
Figure further demonstrates the zonal mean annual aerosol mixing state indices, highlighting that differences between and tended to be zonally structured, where the MAM4 model overestimated at low latitudes, while it underestimated at high latitudes relative to the ML model. In contrast, the MAM4 model overestimated at all latitudes north of 60 S.
Figure 5
Zonal mean annual aerosol mixing state indices (a) , (b) , and (c) using the MAM4 model and ML model. The bands refer to the standard deviation.
[Figure omitted. See PDF]
5.2 Interpretation of findingsFrom Sect. , the following picture emerges: MAM4 overestimates the mixing state index except in regions at high latitudes in the Southern Hemisphere. At the same time, and are overestimated at low latitudes to mid-latitudes and underestimated at high latitudes. These findings point towards too rapid a transfer from the carbonaceous mode to the accumulation mode at low latitudes to mid-latitudes and too slow a transfer at high latitudes.
To conceptually illustrate these relationships, here we use and as examples and contrast the conditions for high and low latitudes. Figure a–f show conditions representative of high latitudes. A grid cell sampled from the CESM2/MAM4 simulation (73 N, 151 W) contains 15 % BC and 37 % POM, distributed over the accumulation and primary carbon mode as shown in Fig. a and d. The corresponding value for is 80 %. Figure b depicts particle population that was sampled from the MAM4 population in Fig. a. All particles, except for the smallest ones (corresponding to Aitken mode particles), contain BC, which results in the relatively high mixing state index value for . Note that in MAM4 BC is not included in the Aitken mode by definition. Considering the same particle population, but now evaluating the mixing state metric , which quantifies the degree of mixing of primary carbon and other species, yields the following observation. The entire primary carbon mode, by definition, consists of POM and BC, which results in an appreciable number of particles that contain only primary carbon (BC POM), giving a mixing state index of only 27 %.
We now compare the MAM4-sampled particle populations above to particle populations that were sampled from our PartMC scenario library. We searched for populations with similar mass fractions of BC and POM as in the MAM4 populations and that were simulated at a similar latitude as the grid point location of the CESM2/MAM4 model output. Figure c shows that the PartMC results have comparatively more BC-free particles, and Fig. f shows that comparatively more particles are mixtures of primary carbon and other species. Overall, this means that in MAM4 BC appears too internally mixed (because irrespective of whether BC is placed in the primary carbon or accumulation mode, it is by design mixed with other species) and that at high latitudes the primary carbon mode is not transferring mass to the accumulation mode as quickly as is the case in PartMC simulations.
The reason why MAM4 behaves in this way can be explained by the aging process treatment in MAM4. Aging in MAM4 is formulated using a threshold criteria. That is, BC and POA mass is transferred from the primary carbon mode to the accumulation mode when a certain threshold of sulfate and SOA has condensed. In MAM4 this threshold is set to a relatively large value. This is done to prevent BC from being removed too quickly by wet deposition – because the primary carbon mode has a lower hygroscopicity than the accumulation mode and thus a lower wet scavenging efficiency – thereby counteracting a low bias in BC concentrations in the Arctic regions. From we already know that using such a high threshold may not be appropriate. However, the global model also has biases in other processes that contribute to the low BC bias in the Arctic, and setting the threshold to a high value compensates for these errors. Our results are a reflection of this fact. While adjusting the threshold criteria in MAM4 to a lower value may improve the agreement with the ML simulations in some regions, it may deteriorate the overall results in other areas. This is a good example how structural uncertainty manifests itself, namely by the fact that adjusting a parameter does not fundamentally fix the issue.
Figure g–l show conditions representative of low latitudes. A grid cell sampled from the CESM2/MAM4 simulation (20 N, 120 E) contains 11 % BC and 24 % POM, distributed over the accumulation and primary carbon mode as shown in Fig. g and j, with most of the mass in the accumulation mode. The corresponding value for is therefore 99 %, an almost complete internal mixture. For the same reason, is also very high. Similarly to the high-latitude case, Fig. i and l show that the comparable PartMC population has comparatively more BC-free particles and more particles that contain very low amounts of primary carbonaceous material, leading to lower values of both and compared to the MAM4 results.
Figure 6
Illustration to explain the differences in mixing state representation between MAM4 and the ML model at high and low latitudes.
[Figure omitted. See PDF]
5.3 Comparison to observational dataThe question that arises from Sect. and is of course the following: which spatial distribution of aerosol mixing state reflects reality more closely? The validation of simulated mixing state indices with observational data is still challenging since per-particle mass fractions of species are required for calculating the mixing state indices (see Sect. ). These are in principle obtainable from in situ deployments of single-particle mass spectrometers or by using electron microscopy techniques, but their quantitative derivation comes with challenges and is not routinely done, so that only very few data sets exist that allow for a meaningful comparison . Keeping these limitations in mind, reported a qualitative comparison of available measurements of mixing state metrics in locations in developed countries (Paris, France; Pittsburgh, USA; various locations in Japan) with seasonally averaged results from the ML model based on particle-resolved simulations. This showed that the ML model was able to capture the range of values that is consistent with the observations.
We further compared the ML model estimates using recent observations from China. Specifically, we compared and with values from Taizhou and Beijing derived from Single Particle Soot Photometer (SP2) measurements. For both locations, overestimated the observed values, while was in the range of the observations. Specifically, the measured at a suburban site Taizhou from 26 May to 18 June 2017 ranged from 62 % to 82 %. During the same time period (but in the year 2011), the values of were between 63 % and 84 %, while was between 84 % and 96 %. The values at the urban site of Beijing ranged between 55 % and 70 % in winter (from 10 November to 10 December 2016) and varied between 60 % and 75 % in summer (from 18 May to 25 June 2017). Using our simulations of the year 2011, varied from 60 % to 88 % in winter and from 59 % to 83 % in summer. As a comparison, ranged from 92 % to 97 % in winter and from 87 % to 95 % in summer. A caveat when comparing and , respectively, with the observations reported in and is that the definition of and included BC-free particles, while the values in the measurements by and were calculated only considering the subpopulation of BC-containing particles. This might introduce a bias in the mixing state index between the index used in this paper and the observations (depending on the fraction of the BC-free particles present at any given location).
We can also relate our index qualitatively to the SP2 measurements in the Finnish Arctic during winter 2011–2012 . Although this study did not provide quantitative mixing state index calculations, it is an important finding that BC-containing particles (with various amounts of coatings) co-existed with BC-free particles. As we saw in Sect. , this condition can easily be represented with a particle-resolved approach. However, the modal model with modes configured as in MAM4 puts black carbon in all accumulation-sized particles (Fig. ), which is not consistent with the observations.
6 Conclusions
In this paper we present a framework for evaluating the error in submicron aerosol mixing state induced by aerosol representation assumptions, which is one of the important contributors to structural uncertainty in aerosol models. We quantitatively compared mixing state indices for submicron aerosol predicted by the modal model MAM4 within the global model CESM to a machine-learned model based on high-detail particle-resolved simulations. We focused on the mixing of optically absorbing and non-absorbing species (), the mixing of primary carbonaceous with other aerosol species (), and the mixing of hygroscopic and non-hygroscopic species ().
For , the MAM4 modal representation generally overestimated the degree of mixing of BC with other aerosol species. This overestimation is due to the fact that MAM4's choice of modes does not allow for representing BC-free particles in the accumulation and primary carbon modes. This is in contrast to field observations by , which showed that BC and POM may be externally mixed near sources. The implication of this is that, if optical properties are calculated based on the aerosol composition, absorption will be overestimated.
For and , the error tended to be zonally structured, where the MAM4 model overestimated the mixing state indices at low latitudes and underestimated them at high latitudes compared to the ML model. This behavior could be explained by modeling choices in MAM4, in particular that (1) BC is always emitted with POM, (2) no BC-free particles exist in the submicron modes, and (3) dust is always internally mixed with other aerosol species.
Mixing state is an important emergent property that affects the aerosol radiative forcing and aerosol–cloud interactions, but it is not easy to constrain this property globally. To the best of our knowledge, this is the first study that evaluated the spatial distribution of aerosol mixing state as predicted by a global model. Since errors in mixing state predictions propagate into errors in aerosol climate impacts, our findings provide a framework and reference for Earth system model developers and users regarding simulation reliability. For example, this framework can be used to (1) quantify model bias in simulating mixing state in different regions, identifying model structural deficiencies, and (2) provide insights into potential improvements of model process representations for a more realistic simulation of aerosols.
Appendix A Derivation of mode-based aerosol mixing state index
Table details the notation for aerosol mass and mass fractions to calculate using modal information.
To explain how to obtain Eqs. () and () from Eqs. () and (), let us assume that each mode contains particles and the number of species in the population is . The mixing entropy of particle in mode , , is given by A1 The average particle mixing entropy of the entire population (summed over all modes), , is A2
Given that each mode is assumed to be internally mixed, particles within the same mode have the same composition, and we have A3
This results in A4
Therefore, based on Eq. () and the fact that , Eq. () can be rewritten as A5 With the mode-based , the other mixing state quantities can be computed as described in Sect. .
Table A1Aerosol mass and mass fraction definition and notation. The number of modes is ( for MAM4 without the coarse mode), the number of particles in mode is , and the number of species is .
Quantity | Meaning |
---|---|
mass of species in particle frommode | |
total mass of particle from mode | |
total mass of species from mode | |
total mass of mode | |
total mass of species in population | |
total mass of the population | |
mass fraction of species in particle (within mode ) | |
mass fraction of species in mode | |
mass fraction of particle from mode in population | |
mass fraction of mode in population | |
mass fraction of species in population |
Figure A1
Aerosol species mixing ratio (g kg). Accumulation mode: a1; Aitken mode: a2; primary carbon mode: a4. The coarse mode (a3) is not used in this study and therefore omitted in this figure. Black carbon: bc; dust: dst; sea salt: ncl; primary organic matter: pom; secondary organic aerosol: soa; sulfate: so4.
[Figure omitted. See PDF]
Figure A2
Fraction of aerosol species mixing ratio (%). Accumulation mode: a1; Aitken mode: a2; primary carbon mode: a4. The coarse mode (a3) is not used in this study and therefore omitted in this figure. Black carbon: bc; dust: dst; sea salt: ncl; primary organic matter: pom; secondary organic aerosol: soa; sulfate: so4.
[Figure omitted. See PDF]
Code and data availability
Notebooks and data to reproduce the global mixing state index analysis are available at
Author contributions
ZZ, MW, and NR conceptualized the analysis and wrote the manuscript with input from the co-authors. ZZ developed the code, carried out the simulations, and performed the analysis. LZ, PLM, and XL provided scientific suggestions for the manuscript. All authors were involved in helpful discussions and contributed to the manuscript.
Competing interests
Some authors are members of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors have also no other competing interests to declare.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
We would like to acknowledge high-performance computing support from Cheyenne (10.5065/D6RX99HX) provided by NCAR's Computational and Information Systems Laboratory, sponsored by the National Science Foundation. The CESM project is supported primarily by the National Science Foundation. This research is part of the Blue Waters sustained-petascale computing project, which is supported by the National Science Foundation (awards OCI-0725070 and ACI-1238993), the State of Illinois, and as of December 2019 the National Geospatial-Intelligence Agency. Blue Waters is a joint effort of the University of Illinois at Urbana-Champaign and its National Center for Supercomputing Applications. Po-Lun Ma and Xiaohong Liu were supported by the Enabling Aerosol-cloud interactions at GLobal convection-permitting scalES (EAGLES) project (74358), funded by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research, Earth System Model Development program. The Pacific Northwest National Laboratory is operated for the U.S. Department of Energy by Battelle Memorial Institute under contract DE-AC05-76RL01830.
Financial support
This research has been supported by the Office of Biological and Environmental Research (grant no. DE-SC0019192), the NSF Division of Atmospheric and Geospace Sciences (grant no. AGS-1254428), the Office of Biological and Environmental Research (Enabling Aerosol-cloud interactions at GLobal convection-permitting scalES (EAGLES) project (grant no. 74358)), the Office of Advanced Cyberinfrastructure (grant no. OCI-0725070), and the Division of Advanced Cyberinfrastructure (grant no. ACI-1238993).
Review statement
This paper was edited by Qiang Zhang and reviewed by three anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2021. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Aerosol mixing state is an important emergent property that affects aerosol radiative forcing and aerosol–cloud interactions, but it has not been easy to constrain this property globally. This study aims to verify the global distribution of aerosol mixing state represented by modal models. To quantify the aerosol mixing state, we used the aerosol mixing state indices for submicron aerosol based on the mixing of optically absorbing and non-absorbing species (
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details





1 Department of Civil and Environmental Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
2 Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
3 Department of Civil and Environmental Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA; National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL, USA
4 Atmospheric Sciences and Global Change Division, Pacific Northwest National Laboratory, Richland, WA, USA
5 Department of Atmospheric Sciences, Texas A&M University, College Station, TX, USA
6 Department of Atmospheric Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, USA