We investigate issues related to volatility and multi‐generational gas‐phase aging parameterizations affecting the formation and evolution of secondary organic aerosol (SOA) in models. We show that ...when assuming realistic values for the mass accommodation coefficient, experimentally observed SOA evaporation rates imply significantly lower “effective volatility” than those derived from SOA growth in smog chambers, pointing to the role of condensed phase processes and suggesting that models need to use different parameters to describe the formation and evolution of SOA. We develop a new, experimentally driven paradigm to represent SOA as a non‐absorbing semi‐solid with very low “effective volatility.” We modify both a box model and a 3D chemical transport model, to include simplified parameterizations capturing the first‐order effects of gas‐phase fragmentation reactions and investigate the implications of treating SOA as a non‐volatile, non‐absorbing semi‐solid (NVSOA). Box model simulations predict SOA loadings decrease with increasing fragmentation, and similar SOA loadings are calculated in the traditional, semi‐volatile (SVSOA) approach and with the new paradigm (NVSOA) before evaporation reduces loadings of SVSOA. Box‐model‐calculated O:C ratios increase with aging in both the SVSOA and the NVSOA paradigms. Consistent with box model results, 3D model simulations demonstrate that predicted SOA loadings decrease with the addition of fragmentation reactions. The NVSOA paradigm predicts higher SOA loadings compared to the SVSOA paradigm over nearly the entire 3D modeling domain, with larger differences close to the surface and in regions where higher dilution favors SVSOA evaporation.
Low effective volatility of SOAGas phase fragmentation reactionBox and regional modeling
It is well known that light absorption from dust and black carbon aerosols has a warming effect on climate while light scattering from sulfate, nitrate, and sea salt aerosols has a cooling effect. ...However, there are large uncertainties associated with light absorption and scattering by different types of organic aerosols, especially in the near‐UV and UV spectral regions. In this paper, we present the results from a systematic laboratory study focused on measuring light absorption by secondary organic aerosols (SOAs) generated from dark α‐pinene + O3 and α‐pinene + NOx + O3 systems in the presence of neutral and acidic sulfate seed aerosols. Light absorption was monitored using photoacoustic spectrometers at four different wavelengths: 355, 405, 532, and 870 nm. Significant light absorption at 355 and 405 nm was observed for the SOA formed from α‐pinene + O3 + NO3 system only in the presence of highly acidic sulfate seed aerosols under dry conditions. In contrast, no absorption was observed when the relative humidity was elevated to greater than 27% or in the presence of neutral sulfate seed aerosols. Organic nitrates in the SOA formed in the presence of neutral sulfate seed aerosols were found to be nonabsorbing, while the light‐absorbing compounds are speculated to be aldol condensation oligomers with nitroxy organosulfate groups that are formed in highly acidic sulfate aerosols. Overall, these results suggest that dark α‐pinene + O3 and α‐pinene + NOx + O3 systems do not form light‐absorbing SOA under typical atmospheric conditions.
Key Points
a‐pinene SOA absorbs UV and near UV radiation
This study examines the uncertainties in simulating mass balance and radiative forcing of mineral dust due to biases in the dust size parameterization. Simulations are conducted quasi-globally (180° ...W–180° E and 60° S–70° N) using the WRF-Chem model with three different approaches to represent dust size distribution (8-bin, 4-bin, and 3-mode). The biases in the 3-mode or 4-bin approaches against a relatively more accurate 8-bin approach in simulating dust mass balance and radiative forcing are identified. Compared to the 8-bin approach, the 4-bin approach simulates similar but coarser size distributions of dust particles in the atmosphere, while the 3-mode approach retains more fine dust particles but fewer coarse dust particles due to its prescribed σg of each mode. Although the 3-mode approach yields up to 10 days of longer dust mass lifetime over the remote oceanic regions than the 8-bin approach, the three size approaches produce a similar dust mass lifetime (3.2 days to 3.5 days) on quasi-global average, reflecting that the global dust mass lifetime is mainly determined by the dust mass lifetime near the dust source regions. With the same global dust emission (~4600 Tg yr−1), the 8-bin approach produces a dust mass loading of 39 Tg, while the 4-bin and 3-mode approaches produce 3% (40.2 Tg) and 25% (49.1 Tg) higher dust mass loading, respectively. The difference in dust mass loading between the 8-bin approach and the 4-bin or 3-mode approaches has large spatial variations, with generally smaller relative difference (<10%) near the surface over the dust source regions. The three size approaches also result in significantly different dry and wet deposition fluxes and number concentrations of dust. The difference in dust aerosol optical depth (AOD) (a factor of 3) among the three size approaches is much larger than their difference (25%) in dust mass loading. Compared to the 8-bin approach, the 4-bin approach yields stronger dust absorptivity, while the 3-mode approach yields weaker dust absorptivity. Overall, on quasi-global average, the three size parameterizations result in a significant difference of a factor of 2~3 in dust surface cooling (−1.02~−2.87 W m−2) and atmospheric warming (0.39~0.96 W m−2) and in a tremendous difference of a factor of ~10 in dust TOA (top of atmosphere) cooling (−0.24~−2.20 W m−2). The impact of different size representations on dust radiative forcing efficiency is smaller. An uncertainty of a factor of 2 is quantified in dust emission estimation due to the different size parameterizations. This study also highlights the uncertainties in modeling dust mass and number loading, deposition fluxes, and radiative forcing resulting from different size parameterizations, and motivates further investigation of the impact of size parameterizations on modeling dust impacts on air quality, climate, and ecosystems.
The authors have decomposed the anthropogenic aerosol radiative forcing into direct contributions from each aerosol species to the planetary energy balance through absorption and scattering of solar ...radiation, indirect effects of anthropogenic aerosol on solar and infrared radiation through droplet and crystal nucleation on aerosol, and semidirect effects through the influence of solar absorption on the distribution of clouds. A three-mode representation of the aerosol in version 5.1 of the Community Atmosphere Model (CAM5.1) yields global annual mean radiative forcing estimates for each of these forcing mechanisms that are within 0.1 W m−2of estimates using a more complex seven-mode representation that distinguishes between fresh and aged black carbon and primary organic matter. Simulating fresh black carbon particles separately from internally mixed accumulation mode particles is found to be important only near fossil fuel sources. In addition to the usual large indirect effect on solar radiation, this study finds an unexpectedly large positive longwave indirect effect (because of enhanced cirrus produced by homogenous nucleation of ice crystals on anthropogenic sulfate), small shortwave and longwave semidirect effects, and a small direct effect (because of cancelation and interactions of direct effects of black carbon and sulfate). Differences between the three-mode and seven-mode versions are significantly larger (up to 0.2 W m−2) when the hygroscopicity of primary organic matter is decreased from 0.1 to 0 and transfer of the primary carbonaceous aerosol to the accumulation mode in the seven-mode version requires more hygroscopic material coating the primary particles. Radiative forcing by cloudborne anthropogenic black carbon is only −0.07 W m−2.
Secondary organic aerosols (SOA) constitute an important fraction of fine-mode atmospheric aerosol mass. Frameworks used to develop SOA parameters from laboratory experiments and subsequently used to ...simulate SOA formation in atmospheric models make many simplifying assumptions about the processes that lead to SOA formation in the interest of computational efficiency. These assumptions can limit the ability of the model to predict the mass, composition, and properties of SOA accurately. In this work, we developed a computationally efficient, process-level model named simpleSOM to represent the chemistry, thermodynamic properties, and microphysics of SOA. simpleSOM simulates multigenerational gas-phase chemistry, phase-state-influenced kinetic gas/particle partitioning, heterogeneous chemistry, oligomerization reactions, and vapor losses to the walls of Teflon chambers. As a case study, we used simpleSOM to simulate SOA formation from the photooxidation of α-pinene. This was done to demonstrate the ability of the model to develop parameters that can reproduce environmental chamber data, to highlight the chemical and microphysical processes within simpleSOM, and discuss implications for SOA formation in chambers and in the real atmosphere. SOA parameters developed from experiments performed in the chamber at the California Institute of Technology (Caltech) reproduced observations of SOA mass yield, O : C, and volatility distribution gathered from other experiments. Sensitivity simulations suggested that multigenerational gas-phase aging contributed to nearly half of all SOA and that in the absence of vapor wall losses, SOA production in the Caltech chamber could be nearly 50% higher. Heterogeneous chemistry did not seem to affect SOA formation over the short timescales for oxidation experienced in the chamber experiments. Simulations performed under atmospherically relevant conditions indicated that the SOA mass yields were sensitive to whether and how oligomerization reactions and the particle phase state were represented in the chamber experiment from which the parameters were developed. simpleSOM provides a comprehensive, process-based framework to consistently model the SOA formation and evolution in box and 3D models.
Secondary organic aerosols (SOA) constitute an important fraction of fine-mode atmospheric aerosol mass.
This paper describes and evaluates a new framework for modeling kinetic gas-particle partitioning of secondary organic aerosol (SOA) that takes into account diffusion and chemical reaction within the ...particle phase. The framework uses a combination of (a) an analytical quasi-steady-state treatment for the diffusion–reaction process within the particle phase for fast-reacting organic solutes, and (b) a two-film theory approach for slow- and nonreacting solutes. The framework is amenable for use in regional and global atmospheric models, although it currently awaits specification of the various gas- and particle-phase chemistries and the related physicochemical properties that are important for SOA formation. Here, the new framework is implemented in the computationally efficient Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) to investigate the competitive growth dynamics of the Aitken and accumulation mode particles. Results show that the timescale of SOA partitioning and the associated size distribution dynamics depend on the complex interplay between organic solute volatility, particle-phase bulk diffusivity, and particle-phase reactivity (as exemplified by a pseudo-first-order reaction rate constant), each of which can vary over several orders of magnitude. In general, the timescale of SOA partitioning increases with increase in volatility and decrease in bulk diffusivity and rate constant. At the same time, the shape of the aerosol size distribution displays appreciable narrowing with decrease in volatility and bulk diffusivity and increase in rate constant. A proper representation of these physicochemical processes and parameters is needed in the next generation models to reliably predict not only the total SOA mass, but also its composition- and number-diameter distributions, all of which together determine the overall optical and cloud-nucleating properties.
Nighttime chemical evolution of aerosol and trace gases in a coal‐fired power plant plume was monitored with the Department of Energy Grumman Gulfstream‐1 aircraft during the 2002 New England Air ...Quality Study field campaign. Quasi‐Lagrangian sampling in the plume at increasing downwind distances and processing times was guided by a constant‐volume balloon that was released near the power plant at sunset. While no evidence of fly ash particles was found, concentrations of particulate organics, sulfate, and nitrate were higher in the plume than in the background air. The enhanced sulfate concentrations were attributed to direct emissions of gaseous H2SO4, some of which had formed new particles as evidenced by enhanced concentrations of nucleation‐mode particles in the plume. The aerosol species were internally mixed and the particles were acidic, suggesting that particulate nitrate was in the form of organic nitrate. The enhanced particulate organic and nitrate masses in the plume were inferred as secondary organic aerosol, which was possibly formed from NO3 radical‐initiated oxidation of isoprene and other trace organic gases in the presence of acidic sulfate particles. Microspectroscopic analysis of particle samples suggested that some sulfate was in the form of organosulfates. Microspectroscopy also revealed the presence of sp2 hybridized C = C bonds, which decreased with increasing processing time in the plume, possibly because of heterogeneous chemistry on particulate organics. Constrained plume modeling analysis of the aircraft and tetroon observations showed that heterogeneous hydrolysis of N2O5 was negligibly slow. These results have significant implications for several issues related to the impacts of power plant emissions on air quality and climate.
Jacobson argues that our statement that "many climate models may overestimate warming by BC" has not been demonstrated. Jacobson challenges our results on the basis that we have misinterpreted some ...model results, omitted optical focusing under high relative humidity conditions and by involatile components, and because our measurements consist of only two locations over short atmospheric time periods. We address each of these arguments, acknowledging important issues and clarifying some misconceptions, and stand by our observations. We acknowledge that Jacobson identified one detail in our experimental technique that places an additional constraint on the interpretation of our observations and reduces somewhat the potential consequences of the stated implications.
A significant fraction of atmospheric particles are composed of inorganic substances that are mixed or coated with organic compounds. The properties and behavior of these particles depend on the ...internal composition and arrangement of the specific constituents in each particle. It is important to know which constituent is on the surface and whether it covers the particle surface partially or entirely. We demonstrate here an instrument consisting of an ultrasensitive single-particle mass spectrometer coupled with a differential mobility analyzer to quantitatively measure in real time individual particle composition, size, density, and shape and to determine which substance is on the surface and whether it entirely covers the particle. For this study, we use NaCl particles completely coated with liquid dioctyl phthalate to generate spherical particles, and NaCl particles partially coated with pyrene, a solid poly aromatic hydrocarbon, to produce aspherical particles with pyrene nodules and an exposed NaCl core. We show that the behavior of the mass spectral intensities as a function of laser fluence yields information that can be used to determine the morphological distribution of individual particle constituents.
The mixing state of soot particles in the atmosphere is of crucial importance for assessing their climatic impact, since it governs their chemical reactivity, cloud condensation nuclei activity, and ...radiative properties. To improve the mixing state representation in models, we present a new approach, the stochastic particle‐resolved model PartMC‐MOSAIC, which explicitly resolves the composition of individual particles in a given population of different types of aerosol particles. This approach tracks the evolution of the mixing state of particles due to emission, dilution, condensation, and coagulation. To make this direct stochastic particle‐based method practical, we implemented a new multiscale stochastic coagulation method. With this method we achieved high computational efficiency for situations when the coagulation kernel is highly nonuniform, as is the case for many realistic applications. PartMC‐MOSAIC was applied to an idealized urban plume case representative of a large urban area to simulate the evolution of carbonaceous aerosols of different types due to coagulation and condensation. For this urban plume scenario we quantified the individual processes that contributed to the aging of the aerosol distribution, illustrating the capabilities of our modeling approach. The results showed for the first time the multidimensional structure of particle composition, which is usually lost in sectional or modal aerosol models.