Consumer, industrial, and commercial product usage is a source of exposure to potentially hazardous chemicals. In addition, cleaning agents, personal care products, coatings, and other volatile chemical products (VCPs), evaporate and react in the atmosphere producing secondary pollutants. Here, we show high air emissions from VCP usage (≥ 14 kg person −1 yr −1 , at least 1.7× higher than current operational estimates) are supported by multiple estimation methods and constraints imposed by ambient levels of ozone, hydroxyl radical (OH) reactivity, and the organic component of fine particulate matter (PM 2.5 ) in Pasadena, California. A near-field model, which estimates human chemical exposure during or in the vicinity of product use, indicates these high air emissions are consistent with organic product usage up to ~75 kg person −1 yr −1 , and inhalation of consumer products could be a non-negligible exposure pathway. After constraining the PM 2.5 yield to 5% by mass, VCPs produce ~41% of the photochemical organic PM 2.5 (1.1 ± 0.3 μ g m −3 ) and ~17% of maximum daily 8-hr average ozone (9 ± 2 ppb) in summer Los Angeles. Therefore, both toxicity and ambient criteria pollutant formation should be considered when organic substituents are developed for VCPs in pursuit of safer and sustainable products and cleaner air.

Abstract. We describe simulations using an updated version of the Community Multiscale Air Quality model version 5.3 (CMAQ v5.3) to investigate the contribution of intermediate-volatility organic compounds (IVOCs) to secondary organic aerosol (SOA) formation in southern California during the CalNex study. We first derive a model-ready parameterization for SOA formation from IVOC emissions from mobile sources. To account for SOA formation from both diesel and gasoline sources, the parameterization has six lumped precursor species that resolve both volatility and molecular structure (aromatic versus aliphatic). We also implement new mobile-source emission profiles that quantify all IVOCs based on direct measurements. The profiles have been released in SPECIATE 5.0. By incorporating both comprehensive mobile-source emission profiles for semivolatile organic compounds (SVOCs) and IVOCs and experimentally constrained SOA yields, this CMAQ configuration best represents the contribution of mobile sources to urban and regional ambient organic aerosol (OA). In the Los Angeles region, gasoline sources emit 4 times more non-methane organic gases (NMOGs) than diesel sources, but diesel emits roughly 3 times more IVOCs on an absolute basis. The revised model predicts all mobile sources (including on- and off-road gasoline, aircraft, and on- and off-road diesel) contribute ∼1 µg m−3 to the daily peak SOA concentration in Pasadena. This represents a ∼70 % increase in predicted daily peak SOA formation compared to the base version of CMAQ. Therefore, IVOCs in mobile-source emissions contribute almost as much SOA as traditional precursors such as single-ring aromatics. However, accounting for these emissions in CMAQ does not reproduce measurements of either ambient SOA or IVOCs. To investigate the potential contribution of other IVOC sources, we performed two exploratory simulations with varying amounts of IVOC emissions from nonmobile sources. To close the mass balance of primary hydrocarbon IVOCs, IVOCs would need to account for 12 % of NMOG emissions from nonmobile sources (or equivalently 30.7 t d−1 in the Los Angeles–Pasadena region), a value that is well within the reported range of IVOC content from volatile chemical products. To close the SOA mass balance and also explain the mildly oxygenated IVOCs in Pasadena, an additional 14.8 % of nonmobile-source NMOG emissions would need to be IVOCs (assuming SOA yields from the mobile IVOCs apply to nonmobile IVOCs). However, an IVOC-to-NMOG ratio of 26.8 % (or equivalently 68.5 t d−1 in the Los Angeles–Pasadena region) for nonmobile sources is likely unrealistically high. Our results highlight the important contribution of IVOCs to SOA production in the Los Angeles region but underscore that other uncertainties must be addressed (multigenerational aging, aqueous chemistry and vapor wall losses) to close the SOA mass balance. This research also highlights the effectiveness of regulations to reduce mobile-source emissions, which have in turn increased the relative importance of other sources, such as volatile chemical products.

Allergic airway diseases represent a complex health problem which can be exacerbated by the synergistic action of pollen particles and air pollutants such as ozone. Understanding human exposures to aeroallergens requires accurate estimates of the spatial distribution of airborne pollen levels as well as of various air pollutants at different times. However, currently there are no established methods for estimating allergenic pollen emissions and concentrations over large geographic areas such as the United States. A mechanistic modeling system for describing pollen emissions and transport over extensive domains has been developed by adapting components of existing regional scale air quality models and vegetation databases. First, components of the Biogenic Emissions Inventory System (BEIS) were adapted to predict pollen emission patterns. Subsequently, the transport module of the Community Multiscale Air Quality (CMAQ) modeling system was modified to incorporate description of pollen transport. The combined model, CMAQ-pollen, allows for simultaneous prediction of multiple air pollutants and pollen levels in a single model simulation, and uses consistent assumptions related to the transport of multiple chemicals and pollen species. Application case studies for evaluating the combined modeling system included the simulation of birch and ragweed pollen levels for the year 2002, during their corresponding peak pollination periods (April for birch and September for ragweed). The model simulations were driven by previously evaluated meteorological model outputs and emissions inventories for the eastern United States for the simulation period. A semi-quantitative evaluation of CMAQ-pollen was performed using tree and ragweed pollen counts in Newark, NJ for the same time periods. The peak birch pollen concentrations were predicted to occur within two days of the peak measurements, while the temporal patterns closely followed the measured profiles of overall tree pollen. For the case of ragweed pollen, the model was able to capture the patterns observed during September 2002, but did not predict an early peak; this can be associated with a wider species pollination window and inadequate spatial information in current land cover databases. An additional sensitivity simulation was performed to comparatively evaluate the dispersion patterns predicted by CMAQ-pollen with those predicted by the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model, which is used extensively in aerobiological studies. The CMAQ estimated concentration plumes matched the equivalent pollen scenario modeled with HYSPLIT. The novel pollen modeling approach presented here allows simultaneous estimation of multiple airborne allergens and other air pollutants, and is being developed as a central component of an integrated population exposure modeling system, the Modeling Environment for Total Risk studies (MENTOR) for multiple, co-occurring contaminants that include aeroallergens and irritants.

Per- and polyfluoroalkyl substances (PFASs) have been released into the environment for decades, yet contributions of air emissions to total human exposure, from inhalation and drinking water contamination via deposition, are poorly constrained. The atmospheric transport and fate of a PFAS mixture from a fluoropolymer manufacturing facility in North Carolina were investigated with the Community Multiscale Air Quality (CMAQ) model applied at high resolution (1 km) and extending ∼150 km from the facility. Twenty-six explicit PFAS compounds, including GenX, were added to CMAQ using current best estimates of air emissions and relevant physicochemical properties. The new model, CMAQ-PFAS, predicts that 5% by mass of total emitted PFAS and 2.5% of total GenX are deposited within ∼150 km of the facility, with the remainder transported out. Modeled air concentrations of total GenX and total PFAS around the facility can reach 24.6 and 8500 ng m–3 but decrease to ∼0.1 and ∼10 ng m–3 at 35 km downwind, respectively. We find that compounds with acid functionality have higher deposition due to enhanced water solubility and pH-driven partitioning to aqueous media. To our knowledge, this is the first modeling study of the fate of a comprehensive, chemically resolved suite of PFAS air emissions from a major manufacturing source.