Particle-Phase Diffusion Modulates Partitioning of Semivolatile Organic Compounds to Aged Secondary Organic Aerosol

Model Framework for Predicting Semivolatile Organic Material Emissions Indoors from Organic Aerosol Measurements: Applications to HOMEChem Stir-Frying

Cooking activities emit myriad low-volatility, semivolatile, and highly volatile organic compounds that together form particles that can accumulate to large indoor concentrations. Absorptive partitioning thermodynamics governs the particle-phase organic aerosol concentration mainly via temperature and sorbing mass impacts. Cooking activities can increase the organic sorbing mass by 1-2 orders of magnitude, increasing particle-phase concentrations and affecting emission rate calculations. Although recent studies have begun to probe the volatility characteristics of indoor cooking particles, parametrizations of cooking particle mass emissions have largely neglected these thermodynamic considerations. Here, we present an improved thermodynamics-based model framework for estimating condensable organic material emission rates from a time series of observed concentrations, given that adequate measurements or assumptions can be made about the volatility of the emitted species. We demonstrate the performance of this methodology by applying data from stir-frying experiments performed during the House Observations of Microbial and Environmental Chemistry (HOMEChem) campaign to a two-zone box model representing the UTest House. Preliminary estimates of organic mass emitted on a per-stir-fry basis for three types of organic aerosol factors are presented. Our analysis highlights that using traditional nonvolatile particle models and emission characterizations for some organic aerosol emitting activities can incorrectly attribute concentration changes to emissions rather than thermodynamic effects.

Nonequilibrium Behavior in Isoprene Secondary Organic Aerosol

Recent studies have shown that instantaneous gas-particle equilibrium partitioning assumptions fail to predict SOA formation, even at high relative humidity (∼85%), and photochemical aging seems to be one driving factor. In this study, we probe the minimum aging time scale required to observe nonequilibrium partitioning of semivolatile organic compounds (SVOCs) between the gas and aerosol phase at ∼50% RH. Seed isoprene SOA is generated by photo-oxidation in the presence of effloresced ammonium sulfate seeds at <1 ppbv NOx, aged photochemically or in the dark for 0.3-6 h, and subsequently exposed to fresh isoprene SVOCs. Our results show that the equilibrium partitioning assumption is accurate for fresh isoprene SOA but breaks down after isoprene SOA has been aged for as short as 20 min even in the dark. Modeling results show that a semisolid SOA phase state is necessary to reproduce the observed particle size distribution evolution. The observed nonequilibrium partitioning behavior and inferred semisolid phase state are corroborated by offline mass spectrometric analysis on the bulk aerosol particles showing the formation of organosulfates and oligomers. The unexpected short time scale for the phase transition within isoprene SOA has important implications for the growth of atmospheric ultrafine particles to climate-relevant sizes.

Effective mass accommodation for partitioning of organic compounds into surface films with different viscosities

Indoor surfaces can act as reservoirs and reaction media influencing the concentrations and type of species that people are exposed to indoors. Mass accommodation and partitioning are impacted by the phase state and viscosity of indoor surface films. We developed the kinetic multi-layer model KM-FILM to simulate organic film formation and growth, but it is computationally expensive to couple such comprehensive models with indoor air box models. Recently, a novel effective mass accommodation coefficient (αeff) was introduced for efficient and effective treatments of gas-particle partitioning. In this study, we extended this approach to a film geometry with αeff as a function of penetration depth into the film, partitioning coefficient, bulk diffusivity, and condensed-phase reaction rate constant. Comparisons between KM-FILM and the αeff method show excellent agreement under most conditions, but with deviations before the establishment of quasi-equilibrium within the penetration depth. We found that the deposition velocity of species and overall film growth are impacted by bulk diffusivity in highly viscous films (Db ∼<10-15 cm2 s-1). Reactions that lead to non-volatile products can increase film thicknesses significantly, with the extent of film growth being dependent on the gas-phase concentration, rate coefficient, partitioning coefficient and diffusivity. Amorphous semisolid films with Db > ∼10-17-10-19 cm2 s-1 can be efficient SVOC reservoirs for compounds with higher partitioning coefficients as they can be released back to the gas phase over extended periods of time, while glassy solid films would not be able to act as reservoirs as gas-film partitioning is impeded.

Vapors Are Lost to Walls, Not to Particles on the Wall: Artifact-Corrected Parameters from Chamber Experiments and Implications for Global Secondary Organic Aerosol

Atmospheric models of secondary organic aerosol (OA) (SOA) typically rely on parameters derived from environmental chambers. Chambers are subject to experimental artifacts, including losses of (1) particles to the walls (PWL), (2) vapors to the particles on the wall (V2PWL), and (3) vapors to the wall directly (VWL). We present a method for deriving artifact-corrected SOA parameters and translating these to volatility basis set (VBS) parameters for use in chemical transport models (CTMs). Our process involves combining a box model that accounts for chamber artifacts (Statistical Oxidation Model with a TwO-Moment Aerosol Sectional model (SOM-TOMAS)) with a pseudo-atmospheric simulation to develop VBS parameters that are fit across a range of OA mass concentrations. We found that VWL led to the highest percentage change in chamber SOA mass yields (high NO_x: 36-680%; low NO_x: 55-250%), followed by PWL (high NO_x: 8-39%; low NO_x: 10-37%), while the effects of V2PWL are negligible. In contrast to earlier work that assumed that V2PWL was a meaningful loss pathway, we show that V2PWL is an unimportant SOA loss pathway and can be ignored when analyzing chamber data. Using our updated VBS parameters, we found that not accounting for VWL may lead surface-level OA to be underestimated by 24% (0.25 μg m^-3) as a global average or up to 130% (9.0 μg m^-3) in regions of high biogenic or anthropogenic activity. Finally, we found that accurately accounting for PWL and VWL improves model-measurement agreement for fine mode aerosol mass concentrations (PM_2.5) in the GEOS-Chem model.

Comparison of Phase States of PM2.5 over Megacities, Seoul and Beijing, and Their Implications on Particle Size Distribution

Although the particle phase state is an important property, there is scant information on it, especially, for real-world aerosols. To explore the phase state of fine mode aerosols (PM_2.5) in two megacities, Seoul and Beijing, we collected PM_2.5 filter samples daily from Dec 2020 to Jan 2021. Using optical microscopy combined with the poke-and-flow technique, the phase states of the bulk of PM_2.5 as a function of relative humidity (RH) were determined and compared to the ambient RH ranges in the two cities. PM_2.5 was found to be liquid to semisolid in Seoul but mostly semisolid to solid in Beijing. The liquid state was dominant on polluted days, while a semisolid state was dominant on clean days in Seoul. These findings can be explained by the aerosol liquid water content related to the chemical compositions of the aerosols at ambient RH; the water content of PM_2.5 was much higher in Seoul than in Beijing. Furthermore, the overall phase states of PM_2.5 observed in Seoul and Beijing were interrelated with the particle size distribution. The results of this study aid in a better understanding of the fundamental physical properties of aerosols and in examining how these are linked to PM_2.5 in polluted urban atmospheres.

Phase state of organic aerosols may limit temperature-driven thermodynamic repartitioning following outdoor-to-indoor transport

Ambient aerosols often experience temperature and humidity gradients following outdoor-to-indoor transport, causing organic aerosols (OA) to either gain or lose mass via gas-particle repartitioning. Recent models have sought to quantify these effects using equilibrium partitioning thermodynamics. However, evidence suggests some indoor OA may possess glassy or semisolid phase states with higher viscosities than liquid OA. Characteristic partitioning timescales of higher-viscosity particles are significantly longer than for liquid particles, which may either fully or partially inhibit repartitioning. For outdoor OA experiencing a temperature change during transport indoors, the ultimate repartitioning state depends on the relationship between the gas-particle partitioning rate coefficient (k_gp) of semivolatile organics and the indoor particle loss rate coefficient (l_p). That is, thermodynamic equilibrium partitioning may occur when semivolatile k_gp ≫ l_p, no repartitioning when semivolatile k_gp ≪ l_p, and partial repartitioning when their magnitudes are similar. Longer indoor particle lifetimes, higher particle number, and larger particle sizes all raise k_gp (driving repartitioning towards equilibrium). For simulated U.S. residences, equilibrium condensation was likely reached in humid climate zones during warm meteorological conditions. In colder regions, the degree of evaporative repartitioning depended on whether organics could repartition before the particle phase state adjusts to indoor conditions, which is uncertain. When an appreciable temperature gradient exists, this study not only confirmed that all outdoor-originating OA that is liquid indoors will reach thermodynamic equilibrium, but also concluded that a plurality (46% for this domain) of such OA that is semisolid may also achieve thermodynamic equilibrium during its indoor lifetime.

Monoterpene Photooxidation in a Continuous-Flow Chamber: SOA Yields and Impacts of Oxidants, NOx, and VOC Precursors

Monoterpene photooxidation plays an important role in secondary organic aerosol (SOA) formation in the atmosphere. The low-volatility products can enhance new particle formation and particle growth and thus influence climate feedback. Here, we present the results of α-pinene and Δ-3-carene photooxidation experiments conducted in continuous-flow mode in an environmental chamber under several reaction conditions. The roles of oxidants, addition of NO, and VOC molecular structure in influencing SOA yield are illustrated. SOA yield from α-pinene photooxidation shows a weak dependence on H₂O₂ concentration, which is a proxy for HO₂ concentration. The high O/C ratios observed in the α-pinene photooxidation products suggest the production of highly oxygenated organic molecules (HOM). Addition of ozone to the chamber during low-NO_x photooxidation experiments leads to higher SOA yield. With the addition of NO, the production of N-containing HOMs is enhanced and the SOA yield shows a modest, nonlinear dependence on the input NO concentration. Carene photooxidation leads to higher SOA yield than α-pinene under similar reaction conditions, which agrees with the lower volatility retrieved from evaporation kinetics experiments. These results improve the understanding of SOA formation from monoterpene photooxidation and could be applied to refine the representation of biogenic SOA formation in models.

Process-Level Modeling Can Simultaneously Explain Secondary Organic Aerosol Evolution in Chambers and Flow Reactors

Secondary organic aerosol (SOA) data gathered in environmental chambers (ECs) have been used extensively to develop parameters to represent SOA formation and evolution. The EC-based parameters are usually constrained to less than one day of photochemical aging but extrapolated to predict SOA aging over much longer timescales in atmospheric models. Recently, SOA has been increasingly studied in oxidation flow reactors (OFRs) over aging timescales of one to multiple days. However, these OFR data have been rarely used to validate or update the EC-based parameters. The simultaneous use of EC and OFR data is challenging because the processes relevant to SOA formation and evolution proceed over very different timescales, and both reactor types exhibit distinct experimental artifacts. In this work, we show that a kinetic SOA chemistry and microphysics model that accounts for various processes, including wall losses, aerosol phase state, heterogeneous oxidation, oligomerization, and new particle formation, can simultaneously explain SOA evolution in EC and OFR experiments, using a single consistent set of SOA parameters. With α-pinene as an example, we first developed parameters by fitting the model output to the measured SOA mass concentration and oxygen-to-carbon (O:C) ratio from an EC experiment (<1 day of aging). We then used these parameters to simulate SOA formation in OFR experiments and found that the model overestimated SOA formation (by a factor of 3-16) over photochemical ages ranging from 0.4 to 13 days, when excluding the abovementioned processes. By comprehensively accounting for these processes, the model was able to explain the observed evolution in SOA mass, composition (i.e., O:C), and size distribution in the OFR experiments. This work suggests that EC and OFR SOA data can be modeled consistently, and a synergistic use of EC and OFR data can aid in developing more refined SOA parameters for use in atmospheric models.

Solid organic-coated ammonium sulfate particles at high relative humidity in the summertime Arctic atmosphere

Physical and chemical properties of individual atmospheric particles determine their climate impacts. Hygroscopic inorganic salt particles mixed with trace amounts of organic material are predicted to be liquid under typical tropospheric conditions in the summertime Arctic. Yet, we unexpectedly observed a significant concentration of solid particles composed of ammonium sulfate with an organic coating under conditions of high relative humidity and low temperature. These particle properties are consistent with marine biogenic-derived new particle formation and growth, with particle collision hypothesized to result in the solid phase. This particle source is predicted to have increasing relevance in the context of declining Arctic sea ice and increasing open water, with impacts on clouds, and therefore climate.

The ability of atmospheric aerosols to impact climate through water uptake and cloud formation is fundamentally determined by the size, composition, and phase (liquid, semisolid, or solid) of individual particles. Particle phase is dependent on atmospheric conditions (relative humidity and temperature) and chemical composition and, importantly, solid particles can inhibit the uptake of water and other trace gases, even under humid conditions. Particles composed primarily of ammonium sulfate are presumed to be liquid at the relative humidities (67 to 98%) and temperatures (−2 to 4 °C) of the summertime Arctic. Under these atmospheric conditions, we report the observation of solid organic-coated ammonium sulfate particles representing 30% of particles, by number, in a key size range (<0.2 µm) for cloud activation within marine air masses from the Arctic Ocean at Utqiaġvik, AK. The composition and size of the observed particles are consistent with recent Arctic modeling and observational results showing new particle formation and growth from dimethylsulfide oxidation to form sulfuric acid, reaction with ammonia, and condensation of marine biogenic sulfate and highly oxygenated organic molecules. Aqueous sulfate particles typically undergo efflorescence and solidify at relative humidities of less than 34%. Therefore, the observed solid phase is hypothesized to occur from contact efflorescence during collision of a newly formed Aitken mode sulfate particle with an organic-coated ammonium sulfate particle. With declining sea ice in the warming Arctic, this particle source is expected to increase with increasing open water and marine biogenic emissions.

Rapid growth of anthropogenic organic nanoparticles greatly alters cloud life cycle in the Amazon rainforest

Aerosol-cloud interactions remain uncertain in assessing climate change. While anthropogenic activities produce copious aerosol nanoparticles smaller than 10 nanometers, they are too small to act as efficient cloud condensation nuclei (CCN). The mechanisms responsible for particle growth to CCN-relevant sizes are poorly understood. Here, we present aircraft observations of rapid growth of anthropogenic nanoparticles downwind of an isolated metropolis in the Amazon rainforest. Model analysis reveals that the sustained particle growth to CCN sizes is predominantly caused by particle-phase diffusion-limited partitioning of semivolatile oxidation products of biogenic hydrocarbons. Cloud-resolving numerical simulations show that the enhanced CCN concentrations in the urban plume substantially alter the formation of shallow convective clouds, suppress precipitation, and enhance the transition to deep convective clouds. The proposed nanoparticle growth mechanism, expressly enabled by the abundantly formed semivolatile organics, suggests an appreciable impact of anthropogenic aerosols on cloud life cycle in previously unpolluted forests of the world.

Heterogeneous Nucleation Drives Particle Size Segregation in Sequential Ozone and Nitrate Radical Oxidation of Catechol

Secondary organic aerosol formation via condensation of organic vapors onto existing aerosol transforms the chemical composition and size distribution of ambient aerosol, with implications for air quality and Earth's radiative balance. Gas-to-particle conversion is generally thought to occur on a continuum between equilibrium-driven partitioning of semivolatile molecules to the pre-existing mass size distribution and kinetic-driven condensation of low volatility molecules to the pre-existing surface area size distribution. However, we offer experimental evidence in contrast to this framework. When catechol is sequentially oxidized by O₃ and NO₃ in the presence of (NH₄)₂SO₄ seed particles with a single size mode, we observe a bimodal organic aerosol mass size distribution with two size modes of distinct chemical composition with nitrocatechol from NO₃ oxidation preferentially condensing onto the large end of the pre-existing size distribution (∼750 nm). A size-resolved chemistry and microphysics model reproduces the evolution of the two distinct organic aerosol size modes─heterogeneous nucleation to an independent, nitrocatechol-rich aerosol phase.

Humidity Dependence of the Condensational Growth of α-Pinene Secondary Organic Aerosol Particles

The influence of relative humidity (RH) on the condensational growth of organic aerosol particles remains incompletely understood. Herein, the RH dependence was investigated via a series of experiments for α-pinene ozonolysis in a continuously mixed flow chamber in which recurring cycles of particle growth occurred every 7 to 8 h at a given RH. In 5 h, the mean increase in the particle mode diameter was 15 nm at 0% RH and 110 nm at 75% RH. The corresponding particle growth coefficients, representing a combination of the thermodynamic driving force and the kinetic resistance to mass transfer, increased from 0.35 to 2.3 nm² s^-1. The chemical composition, characterized by O:C and H:C atomic ratios of 0.52 and 1.48, respectively, and determined by mass spectrometry, did not depend on RH. The Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) was applied to reproduce the observed size- and RH-dependent particle growth by optimizing the diffusivities D_b within the particles of the condensing molecules. The D_b values increased from 5 α^-1 × 10^-16 at 0% RH to 2 α^-1 × 10^-12 cm^-2 s^-1 at 75% RH for mass accommodation coefficients α of 0.1 to 1.0, highlighting the importance of particle-phase properties in modeling the growth of atmospheric aerosol particles.

A core-shell box model for simulating viscosity dependent secondary organic aerosol (CSVA) and its application

Secondary organic aerosol (SOA) plays a key role in air pollution and global climate change. However, the understanding and modelling of SOA properties and evolution are still limited. In this paper, we developed a novel kinetic Core-Shell box model for Viscosity dependent SOA simulation (CSVA), which includes explicit gas-phase reactions (MCM), homogeneous nucleation by H₂SO₄-NH₃-H₂O, viscosity dependent mass transfer between gas and particle phases (organic and aqueous phases) and particle-phase reactions. The gas-particle mass transfer is represented by chainlike reactions analogizing to electrical resistance. The CSVA model is verified and applied to chamber experiments of toluene oxidation systems. The monomers and dimers of SOA are determined by coupling the high-resolution Orbitrap mass spectra and MCM mechanism. The majority of dimers are confirmed to be peroxyhemiacetals formed by reactions of hydroperoxides with aldehydes in the particle phase. The results show that CSVA can well capture the following processes: (1) relative humidity (RH) dependent nucleation of the H₂SO₄-NH₃-H₂O system, (2) particle size-dependent hygroscopic growth of inorganics (e.g., NaCl and (NH₄)₂SO₄) and organics (levoglucosan and SOA), (3) NO_x dependent SOA formation, (4) viscosity-induced evolution of particle size distribution, and (5) effect of RH on SOA formation. In particular, our model reproduces the phenomenon that the evolution of SOA particle size distribution from a one-peak mode into a two-peak mode is due to viscosity.

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (Statistical Oxidation Model coupled to the TwO Moment Aerosol Sectional (SOM-TOMAS)), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (D_b) of SOA produced from α-pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the D_b for SOA can be constrained ((1-7) × 10^-15 cm² s^-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass-transition temperature, viscosity, and, ultimately, D_b (∼10^-15 cm² s^-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size-resolved SOA model to constrain the particle phase state of SOA using historical measurements of the evolution of the particle size distribution.

Indoor aerosol water content and phase state in U.S. residences: impacts of relative humidity, aerosol mass and composition, and mechanical system operation

Hygroscopic particulate matter (PM) constituents promote uptake of aerosol water (AW), depending on relative humidity (RH), which can constrain qualities such as organic aerosol (OA) phase state and inorganic aerosol (IA) deliquescence and efflorescence. This work provides a first incorporation of AW predictions into residential indoor PM simulations. The indoor model, IMAGES, which simulates factored OA concentrations and thermodynamics using the 2D-volatility basis set, was expanded to predict speciated IA concentrations, AW with κ-Köhler theory of hygroscopic growth, and OA phase state with glass transition temperatures. Since RH is the largest driver of AW and varies with meteorology, simulations were conducted using a database of historical ambient weather and pollution records spanning the sixteen U.S. climate zones, facilitating assessment of seasonal and regional trends. Over this diverse simulation set, the residential indoor AW mass was ∼10 to 100 times smaller than dry PM mass. This relative AW amount indoors was about ∼10 times smaller than outdoors, since indoor-emitted aerosol is likely less hygroscopic. The indoor OA phase state was typically semisolid, suggesting kinetic limitations might inhibit thermodynamic OA partitioning equilibrium from being established indoors. Residences in hot and humid climates during the summertime may have liquid indoor OA, while amorphous solid indoor OA can exist in cold climates. Deliquescence and efflorescence of recirculated IA within HVAC systems during cooling or heating, respectively, was also modeled. Oftentimes, two IA populations with different histories existing as wet or dry aerosol were generated by HVAC operation depending on indoor and outdoor environmental conditions and the HVAC operating mode.

Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model

Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T _g), oxygen to carbon (O : C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core-shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~ 70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10²-10¹² Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity < 10² Pa s. Phase separation while in a liquid phase state, i.e., liquid-liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM_2.5) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.