Organic nitrates contribute significantly to the total organic aerosol burden. However, the formation and loss mechanisms of particulate organic nitrates (PONs) remain poorly understood. In this study, with the CMAQ modeling system, we implement a detailed biogenic volatile organic carbon gas phase oxidation mechanism and an explicit representation of multiphase organic nitrate formation and loss, including both aqueous-phase uptake and vapor-pressure driven partitioning into organic aerosol as well as condensed-phase reactions. We find vapor-pressure dependent partitioning is the leading mechanism for formation of PONs and hydrolysis is a major loss mechanism for PON resulting in substantial amounts of organic aerosol that originate as an organic nitrate. Partitioning and hydrolysis together can produce high concentrations (up to 5 μg/m³) of PON-derived aerosols over the southeast United States. The main source of PON-derived aerosols is monoterpene nitrates that have been chemically processed to lose their nitrate functionality through aqueous chemistry. In contrast, the major portion of aqueous aerosol and in-cloud PON, which retains its nitrate moiety, are soluble isoprene nitrates. We evaluate the model using the observations from the Southern Oxidant and Aerosol Study (SOAS) campaign in the Southeast US in summer 2013 and show implementing aerosol-phase pathways for organic nitrates dramatically improves the magnitude of total alkyl nitrates (ANs) in CMAQ. The contribution of PONs to the total ANs at the SOAS site is estimated to be ~20%, a value in the range of the measurements. The predicted AN composition is shifted from monoterpene to isoprene and anthropogenic organic nitrates.

Atmospheric models that accurately describe the fate and transport of trace species for the right reasons aid in the development of effective air-quality management strategies that safeguard human health. Controllable emissions facilitate the formation of biogenic secondary organic aerosol (BSOA) to enhance the atmospheric fine particulate matter (PM_2.5) burden. Previous modeling with the EPA's Community Multiscale Air Quality (CMAQ) model predicted that anthropogenic primary organic aerosol (POA) emissions had the greatest impact on BSOA. That experiment included formation processes involving semivolatile partitioning but not aerosol liquid water (ALW), a ubiquitous PM constituent. We conduct 17 summertime CMAQ simulations with updated chemistry and evaluate changes in BSOA due to the removal of individual pollutants and source sectors for the contiguous U.S. CMAQ predicts SO₂ from electricity generating units, and mobile source NO_X emissions have the largest impacts on BSOA. The removal of anthropogenic NO_X, SO₂, and POA emissions during the simulation reduces the nationally averaged BSOA by 23, 14, and 8% and PM_2.5 by 9.2, 14, and 5.3%, respectively. ALW mass concentrations decrease by 10 and 35% in response to the removal of NO_X and SO₂ emissions. This work contributes chemical insight into ancillary benefits of Federal NO_X and SO₂ rules that concurrently reduce organic PM_2.5 mass.

The lack of statistically robust relationships between IEPOX (isoprene epoxydiol)-derived SOA (IEPOX SOA) and aerosol liquid water and pH observed during the 2013 Southern Oxidant and Aerosol Study (SOAS) emphasizes the importance of modeling the whole system to understand the controlling factors governing IEPOX SOA formation. We present a mechanistic modeling investigation predicting IEPOX SOA based on Community Multiscale Air Quality (CMAQ) model algorithms and a recently introduced photochemical box model, simpleGAMMA. We aim to (1) simulate IEPOX SOA tracers from the SOAS Look Rock ground site, (2) compare the two model formulations, (3) determine the limiting factors in IEPOX SOA formation, and (4) test the impact of a hypothetical sulfate reduction scenario on IEPOX SOA. The estimated IEPOX SOA mass variability is in similar agreement (r² ∼ 0.6) with measurements. Correlations of the estimated and measured IEPOX SOA tracers with observed aerosol surface area (r² ∼ 0.5-0.7), rate of particle-phase reaction (r² ∼ 0.4-0.7), and sulfate (r² ∼ 0.4-0.5) suggest an important role of sulfate in tracer formation via both physical and chemical mechanisms. A hypothetical 25% reduction of sulfate results in ∼70% reduction of IEPOX SOA formation, reaffirming the importance of aqueous phase chemistry in IEPOX SOA production.

Several models were used to describe the partitioning of ammonia, water, and organic compounds between the gas and particle phases for conditions in the southeastern US during summer 2013. Existing equilibrium models and frameworks were found to be sufficient, although additional improvements in terms of estimating pure-species vapor pressures are needed. Thermodynamic model predictions were consistent, to first order, with a molar ratio of ammonium to sulfate of approximately 1.6 to 1.8 (ratio of ammonium to 2× sulfate, R_N/2S ≈ 0.8 to 0.9) with approximately 70% of total ammonia and ammonium (NH _x ) in the particle. Southeastern Aerosol Research and Characterization Network (SEARCH) gas and aerosol and Southern Oxidant and Aerosol Study (SOAS) Monitor for AeRosols and Gases in Ambient air (MARGA) aerosol measurements were consistent with these conditions. CMAQv5.2 regional chemical transport model predictions did not reflect these conditions due to a factor of 3 overestimate of the nonvolatile cations. In addition, gas-phase ammonia was overestimated in the CMAQ model leading to an even lower fraction of total ammonia in the particle. Chemical Speciation Network (CSN) and aerosol mass spectrometer (AMS) measurements indicated less ammonium per sulfate than SEARCH and MARGA measurements and were inconsistent with thermodynamic model predictions. Organic compounds were predicted to be present to some extent in the same phase as inorganic constituents, modifying their activity and resulting in a decrease in [H⁺]_air (H⁺ in μgm^-3 air), increase in ammonia partitioning to the gas phase, and increase in pH compared to complete organic vs. inorganic liquid-liquid phase separation. In addition, accounting for nonideal mixing modified the pH such that a fully interactive inorganic-organic system had a pH roughly 0.7 units higher than predicted using traditional methods (pH = 1.5 vs. 0.7). Particle-phase interactions of organic and inorganic compounds were found to increase partitioning towards the particle phase (vs. gas phase) for highly oxygenated (O : C≥0.6) compounds including several isoprene-derived tracers as well as levoglu-cosan but decrease particle-phase partitioning for low O: C, monoterpene-derived species.

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.

Mounting evidence from field and laboratory observations coupled with atmospheric model analyses shows that primary combustion emissions of organic compounds dynamically partition between the vapor and particulate phases, especially as near-source emissions dilute and cool to ambient conditions. The most recent version of the Community Multiscale Air Quality model version 5.2 (CMAQv5.2) accounts for the semivolatile partitioning and gas-phase aging of these primary organic aerosol (POA) compounds consistent with experimentally derived parameterizations. We also include a new surrogate species, potential secondary organic aerosol from combustion emissions (pcSOA), which provides a representation of the secondary organic aerosol (SOA) from anthropogenic combustion sources that could be missing from current chemical transport model predictions. The reasons for this missing mass likely include the following: (1) unspeciated semivolatile and intermediate volatility organic compound (SVOC and IVOC, respectively) emissions missing from current inventories, (2) multigenerational aging of organic vapor products from known SOA precursors (e.g., toluene, alkanes), (3) underestimation of SOA yields due to vapor wall losses in smog chamber experiments, and (4) reversible organic compounds-water interactions and/or aqueous-phase processing of known organic vapor emissions. CMAQ predicts the spatially averaged contribution of pcSOA to OA surface concentrations in the continental United States to be 38.6 and 23.6 % in the 2011 winter and summer, respectively. Whereas many past modeling studies focused on a particular measurement campaign, season, location, or model configuration, we endeavor to evaluate the model and important uncertain parameters with a comprehensive set of United States-based model runs using multiple horizontal scales (4 and 12 km), gas-phase chemical mechanisms, and seasons and years. The model with representation of semivolatile POA improves predictions of hourly OA observations over the traditional nonvolatile model at sites during field campaigns in southern California (CalNex, May-June 2010), northern California (CARES, June 2010), the southeast US (SOAS, June 2013; SEARCH, January and July, 2011). Model improvements manifest better correlations (e.g., the correlation coefficient at Pasadena at night increases from 0.38 to 0.62) and reductions in underprediction during the photochemically active afternoon period (e.g., bias at Pasadena from -5.62 to -2.42 μg m^-3). Daily averaged predictions of observations at routine-monitoring networks from simulations over the continental US (CONUS) in 2011 show modest improvement during winter, with mean biases reducing from 1.14 to 0.73μg m^-3, but less change in the summer when the decreases from POA evaporation were similar to the magnitude of added SOA mass. Because the model-performance improvement realized by including the relatively simple pcSOA approach is similar to that of more-complicated parameterizations of OA formation and aging, we recommend caution when applying these more-complicated approaches as they currently rely on numerous uncertain parameters. The pcSOA parameters optimized for performance at the southern and northern California sites lead to higher OA formation than is observed in the CONUS evaluation. This may be due to any of the following: variations in real pcSOA in different regions or time periods, too-high concentrations of other OA sources in the model that are important over the larger domain, or other model issues such as loss processes. This discrepancy is likely regionally and temporally dependent and driven by interferences from factors like varying emissions and chemical regimes.

During the 2013 Southern Oxidant and Aerosol Study, Fourier Transform Infrared Spectroscopy (FTIR) and Aerosol Mass Spectrometer (AMS) measurements of submicron mass were collected at Look Rock (LRK), Tennessee, and Centreville (CTR), Alabama. Carbon monoxide and submicron sulfate and organic mass concentrations were 15-60% higher at CTR than at LRK but their time series had moderate correlations (r~0.5). However, NO_x had no correlation (r=0.08) between the two sites with nighttime-to-early-morning peaks 3~10 times higher at CTR than at LRK. Organic mass (OM) sources identified by FTIR Positive Matrix Factorization (PMF) had three very similar factors at both sites: Fossil Fuel Combustion (FFC) related organic aerosols, Mixed Organic Aerosols (MOA), and Biogenic Organic Aerosols (BOA). The BOA spectrum from FTIR is similar (cosine similarity > 0.6) to that of lab-generated particle mass from the photochemical oxidation of both isoprene and monoterpenes under high NO_x conditions from chamber experiments. The BOA mass fraction was highest during the night at CTR but in the afternoon at LRK. AMS PMF resulted in two similar pairs of factors at both sites and a third nighttime NO_x-related factor (33% of OM) at CTR but a daytime nitrate-related factor (28% of OM) at LRK. NO_x was correlated with BOA and LO-OOA for NO_x concentrations higher than 1 ppb at both sites, producing 0.5 ± 0.1 μg m^-3 for CTR-LO-OOA and 1.0 ± 0.3 μg m^-3 for CTR-BOA above 1 ppb additional biogenic OM for each 1 ppb increase of NO_x.

The role of anthropogenic NO_x emissions in secondary organic aerosol (SOA) production is not fully understood but is important for understanding the contribution of emissions to air quality. Here, we examine the role of organic nitrates (RONO₂) in SOA formation over the Korean Peninsula during the Korea-United States Air Quality field study in Spring 2016 as a model for RONO₂ aerosol in cities worldwide. We use aircraft-based measurements of the particle phase and total (gas + particle) RONO₂ to explore RONO₂ phase partitioning. These measurements show that, on average, one-fourth of RONO₂ are in the condensed phase, and we estimate that ≈15% of the organic aerosol (OA) mass can be attributed to RONO₂. Furthermore, we observe that the fraction of RONO₂ in the condensed phase increases with OA concentration, evidencing that equilibrium absorptive partitioning controls the RONO₂ phase distribution. Lastly, we model RONO₂ chemistry and phase partitioning in the Community Multiscale Air Quality modeling system. We find that known chemistry can account for one-third of the observed RONO₂, but there is a large missing source of semivolatile, anthropogenically derived RONO₂. We propose that this missing source may result from the oxidation of semi- and intermediate-volatility organic compounds and/or from anthropogenic molecules that undergo autoxidation or multiple generations of OH-initiated oxidation.

Organic aerosol (OA) constitutes a significant fraction of atmospheric fine particle mass. However, the precursors and chemical processes responsible for a majority of OA are rarely conclusively identified. We use online observations of hundreds of simultaneously measured molecular components obtained from 15 laboratory OA formation experiments with constraints on their effective saturation vapor concentrations to attribute the VOC precursors and subsequent chemical pathways giving rise to the vast majority of OA mass measured in two forested regions. We find that precursors and chemical pathways regulating OA composition and volatility are dynamic over hours to days, with their variations driven by coupled interactions between multiple oxidants. The extent of physical and photochemical aging, and its modulation by NO_x, were key to a uniquely comprehensive combined composition-volatility description of OA. Our findings thus provide some of the most complete mechanistic-level guidance to the development of OA descriptions in air quality and Earth system models.

Organic compounds and liquid water are major aerosol constituents in the southeast United States (SE US). Water associated with inorganic constituents (inorganic water) can contribute to the partitioning medium for organic aerosol when relative humidities or organic matter to organic carbon (OM/OC) ratios are high such that separation relative humidities (SRH) are below the ambient relative humidity (RH). As OM/OC ratios in the SE US are often between 1.8 and 2.2, organic aerosol experiences both mixing with inorganic water and separation from it. Regional chemical transport model simulations including inorganic water (but excluding water uptake by organic compounds) in the partitioning medium for secondary organic aerosol (SOA) when RH > SRH led to increased SOA concentrations,· particularly at night. Water uptake to the organic phase resulted in even greater SOA concentrations as a result of a positive feedback in which water uptake increased SOA, which further increased aerosol water and organic aerosol. Aerosol properties· such as the OM/OC and hygroscopicity parameter (κ_org), were captured well by the model compared with measurements during the Southern Oxidant and Aerosol Study (SOAS) 2013. Organic nitrates from monoterpene oxidation were predicted to be the least water-soluble semivolatile species in the model, but most biogenically derived semivolatile species in the Community Multiscale Air Quality (CMAQ) model were highly water soluble and expected to contribute to water-soluble organic carbon (WSOC). Organic aerosol and SOA precursors were abundant at night, but additional improvements in daytime organic aerosol are needed to close the model-measurement gap. When taking into account deviations from ideality, including both inorganic (when RH > SRH) and organic water in the organic partitioning medium reduced the mean bias in SOA for routine monitoring networks and improved model performance compared to observations from SOAS. Property updates from this work will be released in CMAQ v5.2.