Search

Environmentally persistent free radicals (EPFRs) play an important role in aerosol effects on air quality and public health, but their atmospheric abundance and sources are poorly understood. We measured EPFRs contained in PM2.5 collected in Fairbanks, Alaska, in winter 2022. We find that EPFR concentrations were enhanced during surface-based inversion and correlate strongly with incomplete combustion markers, including carbon monoxide and elemental carbon (R2 > 0.75). EPFRs exhibit moderately good correlations with PAHs, biomass burning organic aerosols, and potassium (R2 > 0.4). We also observe strong correlations of EPFRs with hydrocarbon-like organic aerosols, Fe and Ti (R2 > 0.6), and single-particle mass spectrometry measurements reveal internal mixing of PAHs, with potassium and iron. These results suggest that residential wood burning and vehicle tailpipes are major sources of EPFRs and nontailpipe emissions, such as brake wear and road dust, may contribute to the stabilization of EPFRs. Exposure to the observed EPFR concentrations (18 ± 12 pmol m-3) would be equivalent to smoking ∼0.4-1 cigarette daily. Very strong correlations (R2 > 0.8) of EPFR with hydroxyl radical formation in surrogate lung fluid indicate that exposure to EPFRs may induce oxidative stress in the human respiratory tract.

The chemical composition and physical properties of secondary organic aerosol (SOA) generated through OH-initiated oxidation of mixtures containing β-myrcene, an acyclic monoterpene, and d-limonene, a cyclic monoterpene, were investigated to assess the extent of the chemical interactions between their oxidation products. The SOA samples were prepared in an environmental smog chamber, and their composition was analyzed offline using ultraperformance liquid chromatography coupled with electrospray ionization high-resolution mass spectrometry (UPLC-ESI-HRMS). Our results suggested that SOA containing β-myrcene showed a higher proportion of oligomeric compounds with low volatility compared to that of SOA from d-limonene. The formula distribution and signal intensities of the mixed SOA could be accurately predicted by a linear combination of the mass spectra of the SOA from individual precursors. Effects of cross-reactions were observed in the distribution of isomeric oxidation products within the mixed SOA, as made evident by chromatographic analysis. On the whole, β-myrcene and d-limonene appear to undergo oxidation by OH largely independently from each other, with only subtle effects from cross-reactions influencing the yields of specific oxidation products.

The indoor air quality of a residential home during winter in Fairbanks, Alaska, was investigated and contrasted with outdoor levels. Twenty-four-hour average indoor and outdoor filter samples were collected from January 17 to February 25, 2022, in a residential area with high outdoor PM2.5 concentrations. The oxidative potential of PM2.5 was determined using the dithiothreitol-depletion assay (OPDTT). For the unoccupied house, the background indoor-to-outdoor (I/O) ratio of mass-normalized OP (OPmDTT), a measure of the intrinsic health-relevant properties of the aerosol, was less than 1 (0.53 ± 0.37), implying a loss of aerosol toxicity as air was transported indoors. This may result from transport and volatility losses driven by the large gradients in temperature (average outdoor temperature of -19°C/average indoor temperature of 21 °C) or relative humidity (average outdoor RH of 78%/average indoor RH of 11%), or both. Various indoor activities, including pellet stove use, simple cooking experiments, incense burning, and mixtures of these activities, were conducted. The experiments produced PM2.5 with a highly variable OPmDTT. PM2.5 from cooking emissions had the lowest OP values, while pellet stove PM2.5 had the highest. Correlations between volume-normalized OPDTT (OPvDTT), relevant to exposure, and indoor PM2.5 mass concentration during experiments were much lower compared to those in outdoor environments. This suggests that mass concentration alone can be a poor indicator of possible adverse effects of various indoor emissions. These findings highlight the importance of considering both the quantity of particles and sources (chemical composition), as health metrics for indoor air quality.

The oxidative potential (OP) of outdoor PM2.5 in wintertime Fairbanks, Alaska, is investigated and compared to those in wintertime Atlanta and Los Angeles. Approximately 40 filter samples collected in January-February 2022 at a Fairbanks residential site were analyzed for OP utilizing dithiothreitol-depletion (OPDTT) and hydroxyl-generation (OPOH) assays. The study-average PM2.5 mass concentration was 12.8 μg/m3, with a 1 h average maximum of 89.0 μg/m3. Regression analysis, correlations with source tracers, and contrast between cold and warmer events indicated that OPDTT was mainly sensitive to copper, elemental carbon, and organic aerosol from residential wood burning, and OPOH to iron and organic aerosol from vehicles. Despite low photochemically-driven oxidation rates, the water-soluble fraction of OPDTT was unusually high at 77%, mainly from wood burning emissions. In contrast to other locations, the Fairbanks average PM2.5 mass concentration was higher than Atlanta and Los Angeles, whereas OPDTT in Fairbanks and Atlanta were similar, and Los Angeles had the highest OPDTT and OPOH. Site differences were observed in OP when normalized by both the volume of air sampled and the particle mass concentration, corresponding to exposure and the intrinsic health-related properties of PM2.5, respectively. The sensitivity of OP assays to specific aerosol components and sources can provide insights beyond the PM2.5 mass concentration when assessing air quality.

The Alaskan Layered Pollution And Chemical Analysis (ALPACA) field experiment was a collaborative study designed to improve understanding of pollution sources and chemical processes during winter (cold climate and low-photochemical activity), to investigate indoor pollution, and to study dispersion of pollution as affected by frequent temperature inversions. A number of the research goals were motivated by questions raised by residents of Fairbanks, Alaska, where the study was held. This paper describes the measurement strategies and the conditions encountered during the January and February 2022 field experiment, and reports early examples of how the measurements addressed research goals, particularly those of interest to the residents. Outdoor air measurements showed high concentrations of particulate matter and pollutant gases including volatile organic carbon species. During pollution events, low winds and extremely stable atmospheric conditions trapped pollution below 73 m, an extremely shallow vertical scale. Tethered-balloon-based measurements intercepted plumes aloft, which were associated with power plant point sources through transport modeling. Because cold climate residents spend much of their time indoors, the study included an indoor air quality component, where measurements were made inside and outside a house to study infiltration and indoor sources. In the absence of indoor activities such as cooking and/or heating with a pellet stove, indoor particulate matter concentrations were lower than outdoors; however, cooking and pellet stove burns often caused higher indoor particulate matter concentrations than outdoors. The mass-normalized particulate matter oxidative potential, a health-relevant property measured here by the reactivity with dithiothreiol, of indoor particles varied by source, with cooking particles having less oxidative potential per mass than pellet stove particles.

Tea is one of the worlds most popular and widely consumed beverages. It is a common pastime to enjoy a cup of tea in the sunshine. However, little attention has been given to understanding the possible photochemical reactions occurring beneath the calm surface of brewed tea. Epigallocatechin gallate (EGCG), which is widely used in food and beverages, is the most significant active ingredient found in tea. In this study, we investigated the presence of free radicals in both an aqueous EGCG solution and brewed tea under simulated sunlight conditions. To our surprise, we unexpectedly observed the production of hydroxyl radicals (•OH) in brewed tea. It was found that sunlight irradiation played a critical role in the formation of •OH, independent of the presence of metal ions. Furthermore, we demonstrated that the •OH generated from the EGCG aqueous solution induced cell cytotoxicity and DNA damage in vitro. Considering the crucial role of •OH in various fields, including human health and the environment, it is important to further explore the practical implications of •OH production in brewed tea under sunlight. In summary, our study unveils the unexpected formation of •OH in brewed tea and emphasizes the significance of sunlight-induced reactions. The observed cytotoxic and DNA-damaging effects of •OH emphasize the importance of understanding the potential health consequences associated with tea consumption. Further research in this area will contribute to a better understanding of the broader implications of •OH production in brewed tea under sunlight.

BACKGROUND: Although many studies have linked prenatal exposure to PM2.5 to adverse birth outcomes, little is known about the effects of exposure to specific constituents of PM2.5 or mechanisms that contribute to these outcomes. OBJECTIVES: Our objective was to investigate effects of oxidative potential and PM2.5 metal components from non-exhaust traffic emissions, such as brake and tire wear, on the risk of preterm birth (PTB) and term low birth weight (TLBW). METHODS: For a birth cohort of 285,614 singletons born in Los Angeles County, California, in the period 2017-2019, we estimated speciated PM2.5 exposures modeled from land use regression with cokriging, including brake and tire wear related metals (barium and zinc), black carbon, and three markers of oxidative potential (OP), including modeled reactive oxygen species based on measured iron and copper (ROS), OH formation (OPOH), and dithiothreitol (DTT) loss (OPDTT). Using logistic regression, we estimated odds ratios (OR) and 95% confidence intervals (CI) for PTB and TLBW with speciated PM2.5 exposures and PM2.5 mass as continuous variables scaled by their interquartile range (IQR). RESULTS: For both metals and oxidative potential metrics, we estimated increased risks for PTB (ORs ranging from 1.01 to 1.03) and TLBW (ORs ranging from 1.02 to 1.05) per IQR exposure increment that were robust to adjustment for PM2.5 mass. Associations for PM2.5 mass, black carbon, metal components, and oxidative potential (especially ROS and OPOH) with adverse birth outcomes were stronger in Hispanic, Black, and mixed-race or Native American women. DISCUSSION: Our results indicate that exposure to PM2.5 metals from brake and tire wear and particle components that contribute to oxidative potential were associated with an increased risk of PTB and TLBW in Los Angeles County, particularly among Hispanic, Black, and mixed-race or Native American women. Thus, reduction of PM2.5 mass only may not be sufficient to protect the most vulnerable pregnant women and children from adverse effects due to traffic source exposures. https://doi.org/10.1289/EHP12196.

UNLABELLED: Health effects of oxidant gases may be enhanced by components of particulate air pollution that contribute to oxidative stress. Our aim was to examine if within-city spatial variations in the oxidative potential of outdoor fine particulate air pollution (PM2.5) modify relationships between oxidant gases and cardiovascular mortality. METHODS: We conducted a retrospective cohort study of participants in the Canadian Census Health and Environment Cohort who lived in Toronto or Montreal, Canada, from 2002 to 2015. Cox proportional hazards models were used to estimate associations between outdoor concentrations of oxidant gases (Ox, a redox-weighted average of nitrogen dioxide and ozone) and cardiovascular deaths. Analyses were performed across strata of two measures of PM2.5 oxidative potential and reactive oxygen species concentrations (ROS) adjusting for relevant confounding factors. RESULTS: PM2.5 mass concentration showed little within-city variability, but PM2.5 oxidative potential and ROS were more variable. Spatial variations in outdoor Ox were associated with an increased risk of cardiovascular mortality [HR per 5 ppb = 1.028, 95% confidence interval (CI): 1.001, 1.055]. The effect of Ox on cardiovascular mortality was stronger above the median of each measure of PM2.5 oxidative potential and ROS (e.g., above the median of glutathione-based oxidative potential: HR = 1.045, 95% CI: 1.009, 1.081; below median: HR = 1.000, 95% CI: 0.960, 1.043). CONCLUSION: Within-city spatial variations in PM2.5 oxidative potential may modify long-term cardiovascular health impacts of Ox. Regions with elevated Ox and PM2.5 oxidative potential may be priority areas for interventions to decrease the population health impacts of outdoor air pollution.

The photochemical aging of biomass-burning organic aerosols (BBOAs) by exposure to sunlight changes the chemical composition over its atmospheric lifetime, affecting the toxicological and climate-relevant properties of BBOA particles. This study used electron paramagnetic resonance (EPR) spectroscopy with a spin-trapping agent, 5-tert-butoxycarbonyl-5-methyl-1-pyrroline-N-oxide (BMPO), high-resolution mass spectrometry, and kinetic modeling to study the photosensitized formation of reactive oxygen species (ROS) and free radicals in mixtures of benzoquinone and levoglucosan, known BBOA tracer molecules. EPR analysis of irradiated benzoquinone solutions showed dominant formation of hydroxyl radicals (•OH), which are known products of reaction of triplet-state benzoquinone with water, also yielding semiquinone radicals. In addition, hydrogen radicals (H•) were also observed, which were not detected in previous studies. They were most likely generated by photochemical decomposition of semiquinone radicals. The irradiation of mixtures of benzoquinone and levoglucosan led to substantial formation of carbon- and oxygen-centered organic radicals, which became prominent in mixtures with a higher fraction of levoglucosan. High-resolution mass spectrometry permitted direct observation of BMPO-radical adducts and demonstrated the formation of •OH, semiquinone radicals, and organic radicals derived from oxidation of benzoquinone and levoglucosan. Mass spectrometry also detected superoxide radical adducts (BMPO-OOH) that did not appear in the EPR spectra. Kinetic modeling of the processes in the irradiated mixtures successfully reproduced the time evolution of the observed formation of the BMPO adducts of •OH and H• observed with EPR. The model was then applied to describe photochemical processes that would occur in mixtures of benzoquinone and levoglucosan in the absence of BMPO, predicting the generation of HO2• due to the reaction of H• with dissolved oxygen. These results imply that photoirradiation of aerosols containing photosensitizers induces ROS formation and secondary radical chemistry to drive photochemical aging of BBOA in the atmosphere.

Oxidative stress mediated by reactive oxygen species (ROS) is a key process for adverse aerosol health effects. Secondary organic aerosols (SOA) account for a major fraction of fine particulate matter, and their inhalation and deposition into the respiratory tract causes the formation of ROS by chemical and cellular processes, but their relative contributions are hardly quantified and their link to oxidative stress remains uncertain. Here, we quantified cellular and chemical superoxide generation by 9,10-phenanthrenequinone (PQN) and isoprene SOA using a chemiluminescence assay combined with electron paramagnetic resonance spectroscopy as well as kinetic modeling. We also applied cellular imaging techniques to study the cellular mechanism of superoxide release and oxidative damage on cell membranes. We show that PQN and isoprene SOA activate NADPH oxidase in macrophages to release massive amounts of superoxide, overwhelming the superoxide formation by aqueous chemical reactions in the epithelial lining fluid. The activation dose for PQN is 2 orders of magnitude lower than that of isoprene SOA, suggesting that quinones are more toxic. While higher exposures trigger cellular antioxidant response elements, the released ROS induce oxidative damage to the cell membrane through lipid peroxidation. Such mechanistic and quantitative understandings provide a basis for further elucidation of adverse health effects and oxidative stress by fine particulate matter.

Reactive oxygen species (ROS) and environmentally persistent free radicals (EPFR) play an important role in chemical transformation of atmospheric aerosols and adverse aerosol health effects. This study investigated the effects of nitrogen oxides (NOx) during photooxidation of α-pinene and naphthalene on the EPFR content and ROS formation from secondary organic aerosols (SOA). Electron paramagnetic resonance (EPR) spectroscopy was applied to quantify EPFR content and ROS formation. While no EPFR were detected in α-pinene SOA, we found that naphthalene SOA contained about 0.7 pmol μg-1 of EPFR, and NOx has little influence on EPFR concentrations and oxidative potential. α-Pinene and naphthalene SOA generated under low NOx conditions form OH radicals and superoxide in the aqueous phase, which was lowered substantially by 50-80% for SOA generated under high NOx conditions. High-resolution mass spectrometry analysis showed the substantial formation of nitroaromatics and organic nitrates in a high NOx environment. The modeling results using the GECKO-A model that simulates explicit gas-phase chemistry and the radical 2D-VBS model that treats autoxidation predicted reduced formation of hydroperoxides and enhanced formation of organic nitrates under high NOx due to the reactions of peroxy radicals with NOx instead of their reactions with HO2. Consistently, the presence of NOx resulted in the decrease of peroxide contents and oxidative potential of α-pinene SOA.

Secondary organic aerosol (SOA) plays a critical, yet uncertain, role in air quality and climate. Once formed, SOA is transported throughout the atmosphere and is exposed to solar UV light. Information on the viscosity of SOA, and how it may change with solar UV exposure, is needed to accurately predict air quality and climate. However, the effect of solar UV radiation on the viscosity of SOA and the associated implications for air quality and climate predictions is largely unknown. Here, we report the viscosity of SOA after exposure to UV radiation, equivalent to a UV exposure of 6 to 14 d at midlatitudes in summer. Surprisingly, UV-aging led to as much as five orders of magnitude increase in viscosity compared to unirradiated SOA. This increase in viscosity can be rationalized in part by an increase in molecular mass and oxidation of organic molecules constituting the SOA material, as determined by high-resolution mass spectrometry. We demonstrate that UV-aging can lead to an increased abundance of aerosols in the atmosphere in a glassy solid state. Therefore, UV-aging could represent an unrecognized source of nuclei for ice clouds in the atmosphere, with important implications for Earth's energy budget. We also show that UV-aging increases the mixing times within SOA particles by up to five orders of magnitude throughout the troposphere with important implications for predicting the growth, evaporation, and size distribution of SOA, and hence, air pollution and climate.

The reactions of sulfur dioxide (SO2) with surface-bound compounds on atmospheric aerosols lead to the formation of organic sulfur (OS) compounds, thereby affecting the air quality and climate. Here, we show that the heterogeneous reaction of SO2 with authentic urban grime under near-ultraviolet sunlight irradiation leads to a large suite of various organic compounds including OS released in the gas phase. Calculations indicate that at the core area of Guangzhou, building surface uptake of SO2 is 15 times larger than uptake of SO2 on aerosol surfaces, yielding ~20 ng m-3 of OS that represents an important fraction of the observed OS compounds (60 to 200 ng m-3) in ambient aerosols of Chinese megacities. This chemical pathway occurring during daytime can contribute to the observed fraction of OS compounds in aerosols and improve the understanding of haze formation and urban air pollution.

Emerging contaminants are of concern due to their rapidly increasing numbers and potential ecological and human health effects. In this study, the synergistic effects of the presence of multifunctional nitro, amino and carbon-carbon double bond (C═C) groups on the gas phase ozonolysis in O2 or at the air/solid interface were investigated using five simple model compounds. The gas phase ozonolysis rate constants at 296 K were (3.5 ± 0.9) × 10-20 cm3 molecule-1 s-1 for 2-methyl-1-nitroprop-1-ene and (6.8 ± 0.8) × 10-19 cm3 molecule-1 s-1 for 4-methyl-4-nitro-1-pentene, with lifetimes of 134 and 7 days in the presence of 100 ppb ozone in the atmosphere, respectively. The rate constants for gas phase E-N,N-dimethyl-1-propenylamine and N,N-dimethylallylamine reactions with ozone were too fast (>10-18 cm3 molecule-1 s-1) to be measured, implying lifetimes of less than 5 days. A multiphase kinetics model (KM-GAP) was used to probe the gas-solid kinetics of 1-dimethylamino-2-nitroethylene, yielding a rate constant for the surface reaction of 1.8 × 10-9 cm2 molecule-1 s-1 and in the bulk 1× 10-16 cm3 molecule-1 s-1. These results show that a nitro group attached to the C═C lowers the gas phase rate constant by 2-3 orders of magnitude compared to the simple alkenes, while amino groups have the opposite effect. The presence of both groups provides counterbalancing effects. Products with deleterious health effects including dimethylformamide and formaldehyde were identified by FTIR. The identified products differentiate whether the initial site of ozone attack is C═C and/or the amino group. This study provides a basis for predicting the environmental fates of emerging contaminants and shows that both the toxicity of both the parent compounds and the products should be taken into account in assessing their environmental impacts.

Reactive oxygen species (ROS) play a critical role in the chemical transformation of atmospheric secondary organic aerosols (SOA) and aerosol health effects by causing oxidative stress in vivo. Acidity is an important physicochemical property of atmospheric aerosols, but its effects on the ROS formation from SOA have been poorly characterized. By applying the electron paramagnetic resonance spin-trapping technique and the Diogenes chemiluminescence assay, we find highly distinct radical yields and composition at different pH values in the range of 1-7.4 from SOA generated by oxidation of isoprene, α-terpineol, α-pinene, β-pinene, toluene, and naphthalene. We observe that isoprene SOA has substantial hydroxyl radical (•OH) and organic radical yields at neutral pH, which are 1.5-2 times higher compared to acidic conditions in total radical yields. Superoxide (O2 •-) is found to be the dominant species generated by all types of SOAs at lower pH. At neutral pH, α-terpineol SOA exhibits a substantial yield of carbon-centered organic radicals, while no radical formation is observed by aromatic SOA. Further experiments with model compounds show that the decomposition of organic peroxide leading to radical formation may be suppressed at lower pH due to acid-catalyzed rearrangement of peroxides. We also observe 1.5-3 times higher molar yields of hydrogen peroxide (H2O2) in acidic conditions compared to neutral pH by biogenic and aromatic SOA, likely due to enhanced decomposition of α-hydroxyhydroperoxides and quinone redox cycling, respectively. These findings are critical to bridge the gap in understanding ROS formation mechanisms and kinetics in atmospheric and physiological environments.

BackgroundExcess reactive oxygen species (ROS) can cause oxidative stress damaging cells and tissues, leading to adverse health effects in the respiratory tract. Yet, few human epidemiological studies have quantified the adverse effect of early life exposure to ROS on child health. Thus, this study aimed to examine the association of levels of ROS exposure at birth and the subsequent risk of developing common respiratory and allergic diseases in children.Methods1,284 Toronto Child Health Evaluation Questionnaire (T-CHEQ) participants were followed from birth (born between 1996 and 2000) until outcome, March 31, 2016 or loss-to-follow-up. Using ROS data from air monitoring campaigns and land use data in Toronto, ROS concentrations generated in the human respiratory tract in response to inhaled pollutants were estimated using a kinetic multi-layer model. These ROS values were assigned to participants' postal codes at birth. Cox proportional hazards regression models, adjusted for confounders, were then used to estimate hazard ratios (HR) with 95% confidence intervals (CI) per unit increase in interquartile range (IQR).ResultsAfter adjusting for confounders, iron (Fe) and copper (Cu) were not significantly associated with the risk of asthma, allergic rhinitis, nor eczema. However, ROS, a measure of the combined impacts of Fe and Cu in PM2.5, was associated with an increased risk of asthma (HR = 1.11, 95% CI: 1.02-1.21, p

Air pollution is a major risk factor for human health. Chemical reactions in the epithelial lining fluid (ELF) of the human respiratory tract result in the formation of reactive oxygen species (ROS), which can lead to oxidative stress and adverse health effects. We use kinetic modeling to quantify the effects of fine particulate matter (PM2.5), ozone (O3), and nitrogen dioxide (NO2) on ROS formation, interconversion, and reactivity, and discuss different chemical metrics for oxidative stress, such as cumulative production of ROS and hydrogen peroxide (H2O2) to hydroxyl radical (OH) conversion. All three air pollutants produce ROS that accumulate in the ELF as H2O2, which serves as reservoir for radical species. At low PM2.5 concentrations (10 μg m-3), however, Fenton chemistry overwhelms the ROS buffering effect and leads to a tipping point in H2O2 fate, causing a strong nonlinear increase in OH production. This shift in ROS chemistry and the enhanced OH production provide a tentative mechanistic explanation for how the inhalation of PM2.5 induces oxidative stress and adverse health effects.

Rationale: Evidence linking outdoor air pollution with coronavirus disease (COVID-19) incidence and mortality is largely based on ecological comparisons between regions that may differ in factors such as access to testing and control measures that may not be independent of air pollution concentrations. Moreover, studies have yet to focus on key mechanisms of air pollution toxicity such as oxidative stress. Objectives: To conduct a within-city analysis of spatial variations in COVID-19 incidence and the estimated generation of reactive oxygen species (ROS) in lung lining fluid attributable to fine particulate matter (particulate matter with an aerodynamic diameter ⩽2.5 μm [PM2.5]). Methods: Sporadic and outbreak-related COVID-19 case counts, testing data, population data, and sociodemographic data for 140 neighborhoods were obtained from the City of Toronto. ROS estimates were based on a mathematical model of ROS generation in lung lining fluid in response to iron and copper in PM2.5. Spatial variations in long-term average ROS were predicted using a land-use regression model derived from measurements of iron and copper in PM2.5. Data were analyzed using negative binomial regression models adjusting for covariates identified using a directed acyclic graph and accounting for spatial autocorrelation. Measurements and Main Results: A significant positive association was observed between neighborhood-level ROS and COVID-19 incidence (incidence rate ratio = 1.07; 95% confidence interval, 1.01-1.15 per interquartile range ROS). Effect modification by neighborhood-level measures of racialized group membership and socioeconomic status was also identified. Conclusions: Examination of neighborhood characteristics associated with COVID-19 incidence can identify inequalities and generate hypotheses for future studies.

Individual atmospheric particles can contain mixtures of primary organic aerosol (POA), secondary organic aerosol (SOA), and secondary inorganic aerosol (SIA). To predict the role of such complex multicomponent particles in air quality and climate, information on the number and types of phases present in the particles is needed. However, the phase behavior of such particles has not been studied in the laboratory, and as a result, remains poorly constrained. Here, we show that POA+SOA+SIA particles can contain three distinct liquid phases: a low-polarity organic-rich phase, a higher-polarity organic-rich phase, and an aqueous inorganic-rich phase. Based on our results, when the elemental oxygen-to-carbon (O:C) ratio of the SOA is less than 0.8, three liquid phases can coexist within the same particle over a wide relative humidity range. In contrast, when the O:C ratio of the SOA is greater than 0.8, three phases will not form. We also demonstrate, using thermodynamic and kinetic modeling, that the presence of three liquid phases in such particles impacts their equilibration timescale with the surrounding gas phase. Three phases will likely also impact their ability to act as nuclei for liquid cloud droplets, the reactivity of these particles, and the mechanism of SOA formation and growth in the atmosphere. These observations provide fundamental information necessary for improved predictions of air quality and aerosol indirect effects on climate.

This paper presents MATLAB user interfaces for two multiphase kinetic models: the kinetic double-layer model of aerosol surface chemistry and gas--particle interactions (K2-SURF) and the kinetic multilayer model of aerosol surface and bulk chemistry (KM-SUB). Each interface has simple and user-friendly features that allow undergraduate and graduate students in physical, environmental, and atmospheric chemistry classes to learn about multiphase chemistry modeling without prior computer programming or modeling experience. It is easy to input parameters, and the simulation results are promptly displayed in the interface; thus, these model interfaces are particularly suitable for in-classroom and homework teaching applications. The model input parameters include surface and bulk reaction rate coefficients, surface accommodation coefficient, and gas and bulk diffusivities, while model outputs include gas uptake coefficient and surface and bulk concentrations. Students can use the K2-SURF interface to simulate surface processes and the KM-SUB interface to simulate surface and bulk processes. Example simulations were performed for each interface to present atmospherically relevant applications and to demonstrate its versatility for exploring model sensitivity on various kinetic parameters. The K2-SURF interface was used to show how the rate of ozone uptake by an organic surface and temporal evolution of surface concentrations are affected by the surface accommodation coefficient, desorption lifetime, and surface reaction rate coefficient. Additionally, the KM-SUB interface was applied to demonstrate how bulk diffusivities impact the degradation kinetics of oleic acid particles, so that students can learn how the phase state (liquid vs semisolid vs glassy solid) impacts multiphase chemical kinetics. The developed K2-SURF and KM-SUB interfaces are effective tools for modeling surface and bulk reactions in college-level educational settings, helping students to obtain a deeper understanding of the complex behaviors of heterogeneous and multiphase systems.

A critical review of the current state of knowledge of chemical emissions from indoor sources, partitioning among indoor compartments, and the ensuing indoor exposure leads to a proposal for a modular mechanistic framework for predicting human exposure to semivolatile organic compounds (SVOCs). Mechanistically consistent source emission categories include solid, soft, frequent contact, applied, sprayed, and high temperature sources. Environmental compartments are the gas phase, airborne particles, settled dust, indoor surfaces, and clothing. Identified research needs are the development of dynamic emission models for several of the source emission categories and of estimation strategies for critical model parameters. The modular structure of the framework facilitates subsequent inclusion of new knowledge, other chemical classes of indoor pollutants, and additional mechanistic processes relevant to human exposure indoors. The framework may serve as the foundation for developing an open-source community model to better support collaborative research and improve access for application by stakeholders. Combining exposure estimates derived using this framework with toxicity data for different end points and toxicokinetic mechanisms will accelerate chemical risk prioritization, advance effective chemical management decisions, and protect public health.

Historically air constituents have been assumed to be well mixed in indoor environments, with single point measurements and box modeling representing a room or a house. Here we demonstrate that this fundamental assumption needs to be revisited through advanced model simulations and extensive measurements of bleach cleaning. We show that inorganic chlorinated products, such as hypochlorous acid and chloramines generated via multiphase reactions, exhibit spatial and vertical concentration gradients in a room, with short-lived ⋅OH radicals confined to sunlit zones, close to windows. Spatial and temporal scales of indoor constituents are modulated by rates of chemical reactions, surface interactions and building ventilation, providing critical insights for better assessments of human exposure to hazardous pollutants, as well as the transport of indoor chemicals outdoors.

Chemical reactions on indoor surfaces play an important role in air quality in indoor environments, where humans spend 90% of their time. We focus on the challenges of understanding the complex chemistry that takes place on indoor surfaces and identify crucial steps necessary to gain a molecular-level understanding of environmental indoor surface chemistry: (1) elucidate key surface reaction mechanisms and kinetics important to indoor air chemistry, (2) define a range of relevant and representative surfaces to probe, and (3) define the drivers of surface reactivity, particularly with respect to the surface composition, light, and temperature. Within the drivers of surface composition are the roles of adsorbed/absorbed water associated with indoor surfaces and the prevalence, inhomogeneity, and properties of secondary organic films that can impact surface reactivity. By combining laboratory studies, field measurements, and modeling we can gain insights into the molecular processes necessary to further our understanding of the indoor environment.

The neonicotinoid nitenpyram (NPM) is a multifunctional nitroenamine [(R1N)(R2N)C=CHNO2] pesticide. As a nitroalkene, it is structurally similar to other emerging contaminants such as the pharmaceuticals ranitidine and nizatidine. Because ozone is a common atmospheric oxidant, such compounds may be oxidized on contact with air to form new products that have different toxicity compared to the parent compounds. Here we show that oxidation of thin solid films of NPM by gas-phase ozone produces unexpected products, the majority of which do not contain oxygen, despite the highly oxidizing reactant. A further surprising finding is the formation of gas-phase nitrous acid (HONO), a species known to be a major photolytic source of the highly reactive hydroxyl radical in air. The results of application of a kinetic multilayer model show that reaction was not restricted to the surface layers but, at sufficiently high ozone concentrations, occurred throughout the film. The rate constant derived for the O3-NPM reaction is 1 × 10-18 cm3⋅s-1, and the diffusion coefficient of ozone in the thin film is 9 × 10-10 cm2⋅s-1 These findings highlight the unique chemistry of multifunctional nitroenamines and demonstrate that known chemical mechanisms for individual moieties in such compounds cannot be extrapolated from simple alkenes. This is critical for guiding assessments of the environmental fates and impacts of pesticides and pharmaceuticals, and for providing guidance in designing better future alternatives.

Abstract. Information on liquid–liquid phase separation (LLPS) and viscosity (ordiffusion) within secondary organic aerosol (SOA) is needed to improvepredictions of particle size, mass, reactivity, and cloud nucleatingproperties in the atmosphere. Here we report on LLPS and viscosities withinSOA generated by the photooxidation of diesel fuel vapors. Diesel fuelcontains a wide range of volatile organic compounds, and SOA generated bythe photooxidation of diesel fuel vapors may be a good proxy for SOA fromanthropogenic emissions. In our experiments, LLPS occurred over the relativehumidity (RH) range of ∼70 % to ∼100 %,resulting in an organic-rich outer phase and a water-rich inner phase. Theseresults may have implications for predicting the cloud nucleating propertiesof anthropogenic SOA since the presence of an organic-rich outer phase athigh-RH values can lower the supersaturation with respect to water requiredfor cloud droplet formation. At ≤10 % RH, the viscosity was ≥1×108 Pa s, which corresponds to roughly the viscosity of tarpitch. At 38 %–50 % RH, the viscosity was in the range of 1×108 to 3×105 Pa s. These measured viscosities areconsistent with predictions based on oxygen to carbon elemental ratio (O:C)and molar mass as well as predictions based on the number of carbon,hydrogen, and oxygen atoms. Based on the measured viscosities and theStokes–Einstein relation, at ≤10 % RH diffusion coefficients oforganics within diesel fuel SOA is ≤5.4×10-17 cm2 s−1 and the mixing time of organics within 200 nm diesel fuel SOAparticles (τmixing) is 50 h. These small diffusion coefficientsand large mixing times may be important in laboratory experiments, where SOAis often generated and studied using low-RH conditions and on timescales ofminutes to hours. At 38 %–50 % RH, the calculated organic diffusioncoefficients are in the range of 5.4×10-17 to 1.8×10-13 cm2 s−1 and calculated τmixing values arein the range of ∼0.01 h to ∼50 h. These valuesprovide important constraints for the physicochemical properties ofanthropogenic SOA.

Benzo[a]pyrene (BaP), a key polycyclic aromatic hydrocarbon (PAH) often associated with soot particles coated by organic compounds, is a known carcinogen and mutagen. When mixed with organics, the kinetics and mechanisms of chemical transformations of BaP by ozone in indoor and outdoor environments are still not fully elucidated. Using direct analysis in real-time mass spectrometry (DART-MS), kinetics studies of the ozonolysis of BaP in thin films exhibited fast initial loss of BaP followed by a slower decay at long exposure times. Kinetic multilayer modeling demonstrates that the slow decay of BaP over long times can be simulated if there is slow diffusion of BaP from the film interior to the surface, resolving long-standing unresolved observations of incomplete PAH decay upon prolonged ozone exposure. Phase separation drives the slow diffusion time scales in multicomponent systems. Specifically, thermodynamic modeling predicts that BaP phase separates from secondary organic aerosol material so that the BaP-rich layer at the surface shields the inner BaP from ozone. Also, BaP is miscible with organic oils such as squalane, linoleic acid, and cooking oil, but its oxidation products are virtually immiscible, resulting in the formation of a viscous surface crust that hinders diffusion of BaP from the film interior to the surface. These findings imply that phase separation and slow diffusion significantly prolong the chemical lifetime of PAHs, affecting long-range transport of PAHs in the atmosphere and their fates in indoor environments.

Abstract. Secondary organic aerosols (SOA) account for a substantial fraction of airparticulate matter, and SOA formation is often modeled assuming rapidestablishment of gas–particle equilibrium. Here, we estimate thecharacteristic timescale for SOA to achieve gas–particle equilibrium undera wide range of temperatures and relative humidities using astate-of-the-art kinetic flux model. Equilibration timescales werecalculated by varying particle phase state, size, mass loadings, andvolatility of organic compounds in open and closed systems. Modelsimulations suggest that the equilibration timescale for semi-volatilecompounds is on the order of seconds or minutes for most conditions in theplanetary boundary layer, but it can be longer than 1 h if particlesadopt glassy or amorphous solid states with high glass transitiontemperatures at low relative humidity. In the free troposphere with lowertemperatures, it can be longer than hours or days, even at moderate orrelatively high relative humidities due to kinetic limitations of bulkdiffusion in highly viscous particles. The timescale of partitioning oflow-volatile compounds into highly viscous particles is shorter compared tosemi-volatile compounds in the closed system, as it is largely determined bycondensation sink due to very slow re-evaporation with relatively quickestablishment of local equilibrium between the gas phase and thenear-surface bulk. The dependence of equilibration timescales on bothvolatility and bulk diffusivity provides critical insights intothermodynamic or kinetic treatments of SOA partitioning for accuratepredictions of gas- and particle-phase concentrations of semi-volatilecompounds in regional and global chemical transport models.

Indoor surfaces are often coated with organic compounds yet a molecular understanding of what drives these interactions is poorly understood. Herein, the adsorption and desorption of limonene, an organic compound found in indoor environments, on hydroxylated silica (SiO2) surfaces, used to mimic indoor glass surfaces, is investigated by combining vibrational spectroscopy, atomistic computer simulations and kinetic modeling. Infrared spectroscopy shows the interaction involves hydrogen-bonding between limonene and surface O-H groups. Atomistic molecular dynamics (MD) simulations confirm the existence of π-hydrogen bonding interactions, with one or two hydrogen bonds between the silica O-H groups and the carbon-carbon double bonds, roughly one third of the time. The concentration and temperature dependent adsorption/desorption kinetics as measured by infrared spectroscopy were reproduced with a kinetic model, yielding the adsorption enthalpy of ∼55 kJ mol-1, which is consistent with the value derived from the MD simulations. Importantly, this integrated experimental, theoretical and kinetic modeling study constitutes a conceptual framework for understanding the interaction of organic compounds with indoor relevant surfaces and thus provides important insights into our understanding of indoor air chemistry and indoor air quality.

The composition of organic aerosol under different ambient conditions as well as their phase state have been a subject of intense study in recent years. One way to study particle properties is to measure the particle size shrinkage in a diluted environment at isothermal conditions. From these measurements it is possible to separate the fraction of low-volatility compounds from high-volatility compounds. In this work, we analyse and evaluate a method for obtaining particle composition and viscosity from measurements using process models coupled with input optimization algorithms. Two optimization methods, the Monte Carlo genetic algorithm and Bayesian inference, are used together with process models describing the dynamics of particle evaporation. The process model optimization scheme in inferring particle composition in a volatility-basis-set sense and composition-dependent particle viscosity is tested with artificially generated data sets and real experimental data. Optimizing model input so that the output matches these data yields a good match for the estimated quantities. Both optimization methods give equally good results when they are used to estimate particle composition to artificially test data. The timescale of the experiments and the initial particle size are found to be important in defining the range of values that can be identified for the properties from the optimization.

Information on the rate of diffusion of organic molecules within secondary organic aerosol (SOA) is needed to accurately predict the effects of SOA on climate and air quality. Diffusion can be important for predicting the growth, evaporation, and reaction rates of SOA under certain atmospheric conditions. Often, researchers have predicted diffusion rates of organic molecules within SOA using measurements of viscosity and the Stokes-Einstein relation (D ∝ 1/η, where D is the diffusion coefficient and η is viscosity). However, the accuracy of this relation for predicting diffusion in SOA remains uncertain. Using rectangular area fluorescence recovery after photobleaching (rFRAP), we determined diffusion coefficients of fluorescent organic molecules over 8 orders in magnitude in proxies of SOA including citric acid, sorbitol, and a sucrose-citric acid mixture. These results were combined with literature data to evaluate the Stokes-Einstein relation for predicting the diffusion of organic molecules in SOA. Although almost all the data agree with the Stokes-Einstein relation within a factor of 10, a fractional Stokes-Einstein relation (D ∝ 1/ηξ) with ξ = 0:93 is a better model for predicting the diffusion of organic molecules in the SOA proxies studied. In addition, based on the output from a chemical transport model, the Stokes-Einstein relation can overpredict mixing times of organic molecules within SOA by as much as 1 order of magnitude at an altitude of ∼ 3 km compared to the fractional Stokes-Einstein relation with ξ = 0.93. These results also have implications for other areas such as in food sciences and the preservation of biomolecules.

Reactive oxygen species (ROS) play a central role in adverse health effects of air pollutants. Respiratory deposition of fine air particulate matter can lead to the formation of ROS in epithelial lining fluid, potentially causing oxidative stress and inflammation. Secondary organic aerosols (SOA) account for a large fraction of fine particulate matter, but their role in adverse health effects is unclear. Here, we quantify and compare the ROS yields and oxidative potential of isoprene, β-pinene, and naphthalene SOA in water and surrogate lung fluid (SLF). In pure water, isoprene and β-pinene SOA were found to produce mainly OH and organic radicals, whereas naphthalene SOA produced mainly H2O2 and O2•-. The total molar yields of ROS of isoprene and β-pinene SOA were 11.8% and 8.2% in water and decreased to 8.5% and 5.2% in SLF, which can be attributed to ROS removal by lung antioxidants. A positive correlation between the total peroxide concentration and ROS yield suggests that organic (hydro)peroxides may play an important role in ROS formation from biogenic SOA. The total molar ROS yields of naphthalene SOA was 1.7% in water and increased to 11.3% in SLF. This strong increase is likely due to redox reaction cycles involving environmentally persistent free radicals (EPFR) or semiquinones, antioxidants, and oxygen, which may promote the formation of H2O2 and the adverse health effects of anthropogenic SOA from aromatic precursors.

Understanding growth mechanisms for particles in air is fundamental to developing a predictive capability for their impacts on human health, visibility, and climate. In the case of highly viscous semi-solid or solid particles, the likelihood of impinging gases being taken up to grow the particle will be influenced by the initial uptake coefficient and by the residence time of the adsorbed gas on the surface. Here, a new approach that combines Knudsen cell capabilities for gas uptake measurements with temperature programmed desorption (TPD) for binding energy measurements of gases is described. The application of this unique capability to the uptake of organic gases on silica demonstrates its utility and the combination of thermodynamic and kinetic data that can be obtained. Lower limits to the initial net uptake coefficients at 170 K are (3.0 ± 0.6) × 10−3, (4.9 ± 0.6) × 10−3 and (4.3 ± 0.8) × 10−3 for benzene, 1-chloropentane, and methanol, respectively, and are reported here for the first time. The uptake data demonstrated that the ideal gas lattice model was appropriate, which informed the analysis of the TPD data. From the thermal desorption measurements, desorption energies of 34.6 ± 2.5, 45.8 ± 5.5, and 40.0 ± 5.6 kJ mol−1 (errors are 1σ) are obtained for benzene, 1-chloropentane, and methanol, respectively, and show good agreement with previously reported measurements. A multiphase kinetics model was applied to quantify uptake, desorption, and diffusion through the particle multilayers and hence extract desorption kinetics. Implications for uptake of organics on silica surfaces in the atmosphere and the utility of this system for determining relationships between residence times of organic gases and particle surfaces of varying composition are discussed in the context of developing quantitative predictions for growth of aerosol particles in air. [ABSTRACT FROM AUTHOR]

BackgroundCertain metals may play an important role in the adverse health effects of fine particulate air pollution (PM2.5), but few models are available to predict spatial variations in these pollutants.MethodsWe conducted large-scale air monitoring campaigns during summer 2016 and winter 2017 in Toronto, Canada, to characterize spatial variations in iron (Fe) and copper (Cu) concentrations in PM2.5. Information on Fe and Cu concentrations at each site was paired with a kinetic multilayer model of surface and bulk chemistry in the lung epithelial lining fluid to estimate the possible impact of these metals on the production of reactive oxygen species (ROS) in exposed populations. Land use data around each monitoring site were used to develop predictive models for Fe, Cu, and their estimated combined impact on ROS generation.ResultsSpatial variations in Fe, Cu, and ROS greatly exceeded that of PM2.5 mass concentrations. In addition, Fe, Cu, and estimated ROS concentrations were 15, 18, and 9 times higher during summer compared with winter with little difference observed for PM2.5. In leave-one-out cross-validation procedures, final multivariable models explained the majority of spatial variations in annual mean Fe (R 2 = 0.68), Cu (R 2 =0.79), and ROS (R 2 = 0.65).ConclusionsThe combined use of PM2.5 metals data with a kinetic multilayer model of surface and bulk chemistry in the human lung epithelial lining fluid may offer a novel means of estimating PM2.5 health impacts beyond simple mass concentrations.

Secondary organic aerosol (SOA) accounts for a large fraction of submicron particles in the atmosphere. SOA can occur in amorphous solid or semi-solid phase states depending on chemical composition, relative humidity (RH), and temperature. The phase transition between amorphous solid and semi-solid states occurs at the glass transition temperature (Tg). We have recently developed a method to estimate Tg of pure compounds containing carbon, hydrogen, and oxygen atoms (CHO compounds) with molar mass less than 450 g mol-1 based on their molar mass and atomic O: C ratio. In this study, we refine and extend this method for CH and CHO compounds with molar mass up to ∼ 1100 g mol-1 using the number of carbon, hydrogen, and oxygen atoms. We predict viscosity from the Tg-scaled Arrhenius plot of fragility (viscosity vs. Tg •T) as a function of the fragility parameter D. We compiled D values of organic compounds from the literature and found that D approaches a lower limit of ∼ 10 (±1.7) as the molar mass increases. We estimated the viscosity of α-pinene and isoprene SOA as a function of RH by accounting for the hygroscopic growth of SOA and applying the Gordon-Taylor mixing rule, reproducing previously published experimental measurements very well. Sensitivity studies were conducted to evaluate impacts of Tg, D, the hygroscopicity parameter (κ), and the Gordon-Taylor constant on viscosity predictions. The viscosity of toluene SOA was predicted using the elemental composition obtained by high-resolution mass spectrometry (HRMS), resulting in a good agreement with the measured viscosity. We also estimated the viscosity of biomass burning particles using the chemical composition measured by HRMS with two different ionization techniques: electrospray ionization (ESI) and atmospheric pressure photoionization (APPI). Due to differences in detected organic compounds and signal intensity, predicted viscosities at low RH based on ESI and APPI measurements differ by 2-5 orders of magnitude. Complementary measurements of viscosity and chemical composition are desired to further constrain RH-dependent viscosity in future studies.

Polycyclic aromatic hydrocarbons like benzo(a)pyrene (BaP) in atmospheric particulate matter pose a threat to human health because of their high carcinogenicity. In the atmosphere, BaP is mainly degraded through a multiphase reaction with ozone, but the fate and atmospheric transport of BaP are poorly characterized. Earlier modeling studies used reaction rate coefficients determined in laboratory experiments at room temperature, which may overestimate/underestimate degradation rates when applied under atmospheric conditions. Moreover, the effects of diffusion on the particle bulk are not well constrained, leading to large discrepancies between model results and observations. We show how regional and global distributions and transport of BaP can be explained by a new kinetic scheme that provides a realistic description of the temperature and humidity dependence of phase state, diffusivity, and reactivity of BaP-containing particles. Low temperature and humidity can substantially increase the lifetime of BaP and enhance its atmospheric dispersion through both the planetary boundary layer and the free troposphere. The new scheme greatly improves the performance of multiscale models, leading to better agreement with observed BaP concentrations in both source regions and remote regions (Arctic), which cannot be achieved by less-elaborate degradation schemes (deviations by multiple orders of magnitude). Our results highlight the importance of considering temperature and humidity effects on both the phase state of aerosol particles and the chemical reactivity of particulate air pollutants.

The effect of relative humidity (RH) on the chemical composition of secondary organic aerosol (SOA) formed from low-NOx toluene oxidation in the absence of seed particles was investigated. SOA samples were prepared in an aerosol smog chamber at >2%RH and 75%RH, collected on Teflon filters, and analyzed with nanospray desorption electrospray ionization high-resolution mass spectrometry (nano-DESI-HRMS). Measurements revealed a significant reduction in the fraction of oligomers present in the SOA generated at 75%RH compared to SOA generated under dry conditions. In a separate set of experiments, the particle mass concentrations were measured with a scanning mobility particle sizer (SMPS) at RHs ranging from >2 to 90 %. It was found that the particle mass loading decreased by nearly an order of magnitude when RH increased from >2 to 75-90% for low-NOx toluene SOA. The volatility distributions of the SOA compounds, estimated from the distribution of molecular formulas using the "molecular corridor" approach, confirmed that low-NOx toluene SOA became more volatile on average under high-RH conditions. In contrast, the effect of RH on SOA mass loading was found to be much smaller for high-NOx toluene SOA. The observed increase in the oligomer fraction and particle mass loading under dry conditions were attributed to the enhancement of condensation reactions, which produce water and oligomers from smaller compounds in low-NOx toluene SOA. The reduction in the fraction of oligomeric compounds under humid conditions is predicted to partly counteract the previously observed enhancement in the toluene SOA yield driven by the aerosol liquid water chemistry in deliquesced inorganic seed particles.

Terbuthylazine (TBA) is a widely used herbicide, and its heterogeneous reaction with OH radicals is important for assessing its potential to undergo atmospheric long-range transport and to affect the environment and public health. The apparent reaction rate coefficients obtained in different experimental investigations, however, vary by orders of magnitude depending on the applied experimental techniques and conditions. In this study, we used a kinetic multilayer model of aerosol chemistry with reversible surface adsorption and bulk diffusion (KM-SUB) in combination with a Monte Carlo genetic algorithm to simulate the measured decay rates of TBA. Two experimental data sets available from different studies can be described with a consistent set of kinetic parameters resolving the interplay of chemical reaction, mass transport, and shielding effects. Our study suggests that mass transport and shielding effects can substantially extend the atmospheric lifetime of reactive pesticides from a few days to weeks, with strong implications for long-range transport and potential health effects of these substances.

Abstract. Organic compounds present at or near the surface of aqueous droplets can be efficiently oxidized by gas-phase OH radicals, which alter the molecular distribution of the reaction products within the droplet. A change in aerosol composition affects the hygroscopicity and leads to a concomitant response in the equilibrium amount of particle-phase water. The variation in the aerosol water content affects the aerosol size and physicochemical properties, which in turn governs the oxidation kinetics and chemistry. To attain better knowledge of the compositional evolution of aqueous organic droplets during oxidation, this work investigates the heterogeneous OH-radical-initiated oxidation of aqueous methylsuccinic acid (C5H8O4) droplets, a model compound for small branched dicarboxylic acids found in atmospheric aerosols, at a high relative humidity of 85 % through experimental and modeling approaches. Aerosol mass spectra measured by a soft atmospheric pressure ionization source (Direct Analysis in Real Time, DART) coupled with a high-resolution mass spectrometer reveal two major products: a five carbon atom (C5) hydroxyl functionalization product (C5H8O5) and a C4 fragmentation product (C4H6O3). These two products likely originate from the formation and subsequent reactions (intermolecular hydrogen abstraction and carbon–carbon bond scission) of tertiary alkoxy radicals resulting from the OH abstraction occurring at the methyl-substituted carbon site. Based on the identification of the reaction products, a kinetic model of oxidation (a two-product model) coupled with the Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients (AIOMFAC) model is built to simulate the size and compositional changes of aqueous methylsuccinic acid droplets during oxidation. Model results show that at the maximum OH exposure, the droplets become slightly more hygroscopic after oxidation, as the mass fraction of water is predicted to increase from 0.362 to 0.424; however, the diameter of the droplets decreases by 6.1 %. This can be attributed to the formation of volatile fragmentation products that partition to the gas phase, leading to a net loss of organic species and associated particle-phase water, and thus a smaller droplet size. Overall, fragmentation and volatilization processes play a larger role than the functionalization process in determining the evolution of aerosol water content and droplet size at high-oxidation stages.

Abstract. Polycyclic aromatic hydrocarbons (PAHs) are hazardous pollutants, with increasing emissions in pace with economic development in East Asia, but their distribution and fate in the atmosphere are not yet well understood. We extended the regional atmospheric chemistry model WRF-Chem (Weather Research Forecast model with Chemistry module) to comprehensively study the atmospheric distribution and the fate of low-concentration, slowly degrading semivolatile compounds. The WRF-Chem-PAH model reflects the state-of-the-art understanding of current PAHs studies with several new or updated features. It was applied for PAHs covering a wide range of volatility and hydrophobicity, i.e. phenanthrene, chrysene and benzo[a]pyrene, in East Asia. Temporally highly resolved PAH concentrations and particulate mass fractions were evaluated against observations. The WRF-Chem-PAH model is able to reasonably well simulate the concentration levels and particulate mass fractions of PAHs near the sources and at a remote outflow region of East Asia, in high spatial and temporal resolutions. Sensitivity study shows that the heterogeneous reaction with ozone and the homogeneous reaction with the nitrate radical significantly influence the fate and distributions of PAHs. The methods to implement new species and to correct the transport problems can be applied to other newly implemented species in WRF-Chem.