Sensitivity of air quality to vehicle ammonia emissions in the United States.

The US Environmental Protection Agency (EPA) estimates on-road vehicles emissions using the Motor Vehicle Emission Simulator (MOVES). We developed updated ammonia emission rates for MOVES based on road-side exhaust emission measurements of light-duty gasoline and heavy-duty diesel vehicles. The resulting nationwide on-road vehicle ammonia emissions are 1.8, 2.1, 1.8, and 1.6 times higher than the MOVES3 estimates for calendar years 2010, 2017, 2024, and 2035, respectively, primarily due to an increase in light-duty gasoline vehicle NH3 emission rates. We conducted an air quality simulation using the Community Multi-Scale Air Quality (CMAQv5.3.2) model to evaluate the sensitivity of modeled ammonia and fine particulate matter (PM2.5) concentrations in calendar year 2017 using the updated on-road vehicle ammonia emissions. The average monthly urban ammonia ambient concentrations increased by up to 2.3 ppbv in January and 3.0 ppbv in July. The updated on-road NH3 emission rates resulted in better agreement of modeled ammonia concentrations with 2017 annual average ambient ammonia measurements, reducing model bias by 5.8 % in the Northeast region. Modeled average winter PM2.5 concentrations increased in urban areas, including enhancements of up to 0.5 µg/m3 in the northeast United States. The updated ammonia emission rates have been incorporated in MOVES4 and will be used in future versions of the NEI and EPA's modeling platforms.

Volatile Chemical Product Enhancements to Criteria Pollutants in the United States.

Volatile chemical products (VCPs) are a significant source of reactive organic carbon emissions in the United States with a substantial fraction (>20% by mass) serving as secondary organic aerosol (SOA) precursors. Here, we incorporate a new nationwide VCP inventory into the Community Multiscale Air Quality (CMAQ) model with VCP-specific updates to better model air quality impacts. Model results indicate that VCPs mostly enhance anthropogenic SOA in densely populated areas with population-weighted annual average SOA increasing 15-30% in Southern California and New York City due to VCP emissions (contribution of 0.2-0.5 µg m-3). Annually, VCP emissions enhance total population-weighted PM2.5 by ∼5% in California, ∼3% in New York, New Jersey, and Connecticut, and 1-2% in most other states. While the maximum daily 8 h ozone enhancements from VCP emissions are more modest, their influence can cause a several ppb increase on select days in major cities. Printing Inks, Cleaning Products, and Paints and Coatings product use categories contribute ∼75% to the modeled VCP-derived SOA and Cleaning Products, Paints and Coatings, and Personal Care Products contribute ∼81% to the modeled VCP-derived ozone. Overall, VCPs enhance multiple criteria pollutants throughout the United States with the largest impacts in urban cores.

National Exposure Models for Source-Specific Primary Particulate Matter Concentrations Using Aerosol Mass Spectrometry Data.

This paper investigates the feasibility of developing national empirical models to predict ambient concentrations of sparsely monitored air pollutants at high spatial resolution. We used a data set of cooking organic aerosol (COA) and hydrocarbon-like organic aerosol (HOA; traffic primary organic PM) measured using aerosol mass spectrometry across the continental United States. The monitoring locations were selected to span the national distribution of land-use and source-activity variables commonly used for land-use regression modeling (e.g., road length, restaurant count, etc.). The models explain about 60% of the spatial variability of the measured data (R2 0.63 for the COA model and 0.62 for the HOA model). Extensive cross-validation suggests that the models are robust with reasonable transferability. The models predict large urban-rural and intra-urban variability with hotspots in urban areas and along the road corridors. The predicted national concentration surfaces show reasonable spatial correlation with source-specific national chemical transport model (CTM) simulations (R2: 0.45 for COA, 0.4 for HOA). Our measured data, empirical models, and CTM predictions all show that COA concentrations are about two times higher than HOA. Since COA and HOA are important contributors to the intra-urban spatial variability of the total PM2.5, our results highlight the potential importance of controlling commercial cooking emissions for air quality management in the United States.

Characterizing the Air Emissions, Transport, and Deposition of Per- and Polyfluoroalkyl Substances from a Fluoropolymer Manufacturing Facility.

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.

Assessing PM2.5 Model Performance for the Conterminous U.S. with Comparison to Model Performance Statistics from 2007-2015.

Previous studies have proposed that model performance statistics from earlier photochemical grid model (PGM) applications can be used to benchmark performance in new PGM applications. A challenge in implementing this approach is that limited information is available on consistently calculated model performance statistics that vary spatially and temporally over the U.S. Here, a consistent set of model performance statistics are calculated by year, season, region, and monitoring network for PM2.5 and its major components using simulations from versions 4.7.1-5.2.1 of the Community Multiscale Air Quality (CMAQ) model for years 2007-2015. The multi-year set of statistics is then used to provide quantitative context for model performance results from the 2015 simulation. Model performance for PM2.5 organic carbon in the 2015 simulation ranked high (i.e., favorable performance) in the multi-year dataset, due to factors including recent improvements in biogenic secondary organic aerosol and atmospheric mixing parameterizations in CMAQ. Model performance statistics for the Northwest region in 2015 ranked low (i.e., unfavorable performance) for many species in comparison to the 2007-2015 dataset. This finding motivated additional investigation that suggests a need for improved speciation of wildfire PM2.5emissions and modeling of boundary layer dynamics near water bodies. Several limitations were identified in the approach of benchmarking new model performance results with previous results. Since performance statistics vary widely by region and season, a simple set of national performance benchmarks (e.g., one or two targets per species and statistic) as proposed previously are inadequate to assess model performance throughout the U.S. Also, trends in model performance statistics for sulfate over the 2007 to 2015 period suggest that model performance for earlier years may not be a useful reference for assessing model performance for recent years in some cases. Comparisons of results from the 2015 base case with results from five sensitivity simulations demonstrated the importance of parameterizations of NH3 surface exchange, organic aerosol volatility and production, and emissions of crustal cations for predicting PM2.5 species concentrations.

Impacts of Future European Emission Reductions on Aerosol Particle Number Concentrations Accounting for Effects of Ammonia, Amines, and Organic Species.

Although they are currently unregulated, atmospheric ultrafine particles (<100 nm) pose health risks because of, e.g., their capability to penetrate deep into the respiratory system. Ultrafine particles, often minor contributors to atmospheric particulate mass, typically dominate aerosol particle number concentrations. We simulated the response of particle number concentrations over Europe to recent estimates of future emission reductions of aerosol particles and their precursors. We used the chemical transport model PMCAMx-UF, with novel updates including state-of-the-art descriptions of ammonia and dimethylamine new particle formation (NPF) pathways and the condensation of organic compounds onto particles. These processes had notable impacts on atmospheric particle number concentrations. All three emission scenarios (current legislation, optimized emissions, and maximum technically feasible reductions) resulted in substantial (10-50%) decreases in median particle number concentrations over Europe. Consistent reductions were predicted in Central Europe, while Northern Europe exhibited smaller reductions or even increased concentrations. Motivated by the improved NPF descriptions for ammonia and methylamines, we placed special focus on the potential to improve air quality by reducing agricultural emissions, which are a major source of these species. Agricultural emission controls showed promise in reducing ultrafine particle number concentrations, although the change is nonlinear with particle size.

Predictions of PFAS regional-scale atmospheric deposition and ambient air exposure.

Per- and polyfluoroalkyl substances (PFAS) are a large class of human-made compounds that have contaminated the global environment. One environmental entry point for PFAS is via atmospheric emission. Air releases can impact human health through multiple routes, including direct inhalation and contamination of drinking water following air deposition. In this work, we convert the reference dose (RfD) underlying the United States Environmental Protection Agency's GenX drinking water Health Advisory to an inhalation screening level and compare to predicted PFAS and GenX air concentrations from a fluorochemical manufacturing facility in Eastern North Carolina. We find that the area around the facility experiences ~15 days per year of GenX concentrations above the inhalation screening level we derive. We investigate the sensitivity of model predictions to assumptions regarding model spatial resolution, emissions temporal profiles, and knowledge of air emission chemical composition. Decreasing the chemical specificity of PFAS emissions has the largest impact on deposition predictions with domain-wide total deposition varying by as much as 250 % for total PFAS. However, predicted domain-wide mean and median air concentrations varied by <18 % over all scenarios tested for total PFAS. Other model features like emission temporal variability and model spatial resolution had weaker impacts on predicted PFAS deposition.

Reactive organic carbon air emissions from mobile sources in the United States.

Mobile sources are responsible for a substantial controllable portion of the reactive organic carbon (ROC) emitted to the atmosphere, especially in urban environments of the United States. We update existing methods for calculating mobile source organic particle and vapor emissions in the United States with over a decade of laboratory data that parameterize the volatility and organic aerosol (OA) potential of emissions from on-road vehicles, nonroad engines, aircraft, marine vessels, and locomotives. We find that existing emission factor information from Teflon filters combined with quartz filters collapses into simple relationships and can be used to reconstruct the complete volatility distribution of ROC emissions. This new approach consists of source-specific filter artifact corrections and state-of-the-science speciation including explicit intermediate-volatility organic compounds (IVOCs), yielding the first bottom-up volatility-resolved inventory of US mobile source emissions. Using the Community Multiscale Air Quality model, we estimate mobile sources account for 20 %-25 % of the IVOC concentrations and 4.4 %-21.4 % of ambient OA. The updated emissions and air quality model reduce biases in predicting fine-particle organic carbon in winter, spring, and autumn throughout the United States (4.3 %-11.3 % reduction in normalized bias). We identify key uncertain parameters that align with current state-of-the-art research measurement challenges.

Simulations of smog-chamber experiments using the two-dimensional volatility basis set: linear oxygenated precursors.

We use a two-dimensional volatility basis set (2D-VBS) box model to simulate secondary organic aerosol (SOA) mass yields of linear oxygenated molecules: n-tridecanal, 2- and 7-tridecanone, 2- and 7-tridecanol, and n-pentadecane. A hybrid model with explicit, a priori treatment of the first-generation products for each precursor molecule, followed by a generic 2D-VBS mechanism for later-generation chemistry, results in excellent model-measurement agreement. This strongly confirms that the 2D-VBS mechanism is a predictive tool for SOA modeling but also suggests that certain important first-generation products for major primary SOA precursors should be treated explicitly for optimal SOA predictions.

Emissions of Carbonaceous Particulate Matter and Ultrafine Particles from Vehicles-A Scientific Review in a Cross-Cutting Context of Air Pollution and Climate Change.

Airborne particulate matter (PM) is a pollutant of concern not only because of its adverse effects on human health but also on visibility and the radiative budget of the atmosphere. PM can be considered as a sum of solid/liquid species covering a wide range of particle sizes with diverse chemical composition. Organic aerosols may be emitted (primary organic aerosols, POA), or formed in the atmosphere following reaction of volatile organic compounds (secondary organic aerosols, SOA), but some of these compounds may partition between the gas and aerosol phases depending upon ambient conditions. This review focuses on carbonaceous PM and gaseous precursors emitted by road traffic, including ultrafine particles (UFP) and polycyclic aromatic hydrocarbons (PAHs) that are clearly linked to the evolution and formation of carbonaceous species. Clearly, the solid fraction of PM has been reduced during the last two decades, with the implementation of after-treatment systems abating approximately 99% of primary solid particle mass concentrations. However, the role of brown carbon and its radiative effect on climate and the generation of ultrafine particles by nucleation of organic vapour during the dilution of the exhaust remain unclear phenomena and will need further investigation. The increasing role of gasoline vehicles on carbonaceous particle emissions and formation is also highlighted, particularly through the chemical and thermodynamic evolution of organic gases and their propensity to produce particles. The remaining carbon-containing particles from brakes, tyres and road wear will still be a problem even in a future of full electrification of the vehicle fleet. Some key conclusions and recommendations are also proposed to support the decision makers in view of the next regulations on vehicle emissions worldwide.

Human-health impacts of controlling secondary air pollution precursors.

Exposure to ozone and fine particle (PM2.5) air pollution results in premature death. These pollutants are predominantly secondary in nature and can form from nitrogen oxides (NOX), sulfur oxides (SOX), and volatile organic compounds (VOCs). Predicted health benefits for emission reduction scenarios often incompletely account for VOCs as precursors as well as the secondary organic aerosol (SOA) component of PM2.5. Here, we show that anthropogenic VOC emission reductions are more than twice as effective as equivalent fractional reductions of SOX or NOX at reducing air pollution-associated cardiorespiratory mortality in the United States. A 25% reduction in anthropogenic VOC emissions from 2016 levels is predicted to avoid 13,000 premature deaths per year, and most (85%) of the VOC-reduction benefits result from reduced SOA with the remainder from ozone. While NOX (-5.7 ± 0.2 % yr-1) and SOX (-12 ± 1 % yr-1) emissions have declined precipitously across the U.S. since 2002, anthropogenic VOC emissions (-1.8 ± 0.3 % yr-1) and concentrations of non-methane organic carbon (-2.4 ± 1.0 % yr-1) have changed less. This work indicates preferentially controlling VOCs could yield significant benefits to human health.

Predicting the Nonlinear Response of PM2.5 and Ozone to Precursor Emission Changes with a Response Surface Model.

Reducing PM2.5 and ozone concentrations is important to protect human health and the environment. Chemical transport models, such as the Community Multiscale Air Quality (CMAQ) model, are valuable tools for exploring policy options for improving air quality but are computationally expensive. Here, we statistically fit an efficient polynomial function in a response surface model (pf-RSM) to CMAQ simulations over the eastern U.S. for January and July 2016. The pf-RSM predictions were evaluated using out-of-sample CMAQ simulations and used to examine the nonlinear response of air quality to emission changes. Predictions of the pf-RSM are in good agreement with the out-of-sample CMAQ simulations, with some exceptions for cases with anthropogenic emission reductions approaching 100%. NOX emission reductions were more effective for reducing PM2.5 and ozone concentrations than SO2, NH3, or traditional VOC emission reductions. NH3 emission reductions effectively reduced nitrate concentrations in January but increased secondary organic aerosol (SOA) concentrations in July. More work is needed on SOA formation under conditions of low NH3 emissions to verify the responses of SOA to NH3 emission changes predicted here. Overall, the pf-RSM performs well in the eastern U.S., but next-generation RSMs based on deep learning may be needed to meet the computational requirements of typical regulatory applications.

Reactive organic carbon emissions from volatile chemical products.

Volatile chemical products (VCPs) are an increasingly important source of anthropogenic reactive organic carbon (ROC) emissions. Among these sources are everyday items, such as personal care products, general cleaners, architectural coatings, pesticides, adhesives, and printing inks. Here, we develop VCPy, a new framework to model organic emissions from VCPs throughout the United States, including spatial allocation to regional and local scales. Evaporation of a species from a VCP mixture in the VCPy framework is a function of the compound-specific physiochemical properties that govern volatilization and the timescale relevant for product evaporation. We introduce two terms to describe these processes: evaporation timescale and use timescale. Using this framework, predicted national per capita organic emissions from VCPs are 9.5 kg per person per year (6.4 kg C per person per year) for 2016, which translates to 3.05 Tg (2.06 Tg C), making VCPs a dominant source of anthropogenic organic emissions in the United States. Uncertainty associated with this framework and sensitivity to select parameters were characterized through Monte Carlo analysis, resulting in a 95 % confidence interval of national VCP emissions for 2016 of 2.61-3.53 Tg (1.76-2.38 Tg C). This nationwide total is broadly consistent with the U.S. EPA's 2017 National Emission Inventory (NEI); however, county-level and categorical estimates can differ substantially from NEI values. VCPy predicts higher VCP emissions than the NEI for approximately half of all counties, with 5 % of all counties having greater than 55 % higher emissions. Categorically, application of the VCPy framework yields higher emissions for personal care products (150 %) and paints and coatings (25 %) when compared to the NEI, whereas pesticides (-54 %) and printing inks (-13 %) feature lower emissions. An observational evaluation indicates emissions of key species from VCPs are reproduced with high fidelity using the VCPy framework (normalized mean bias of -13 % with r =0.95). Sector-wide, the effective secondary organic aerosol yield and maximum incremental reactivity of VCPs are 5.3 % by mass and 1.58 gO3 g-1, respectively, indicating VCPs are an important, and likely to date underrepresented, source of secondary pollution in urban environments.

Modeling secondary organic aerosol formation from volatile chemical products.

Volatile chemical products (VCPs) are commonly-used consumer and industrial items that are an important source of anthropogenic emissions. Organic compounds from VCPs evaporate on atmospherically relevant time scales and include many species that are secondary organic aerosol (SOA) precursors. However, the chemistry leading to SOA, particularly that of intermediate volatility organic compounds (IVOCs), has not been fully represented in regional-scale models such as the Community Multiscale Air Quality (CMAQ) model, which tend to underpredict SOA concentrations in urban areas. Here we develop a model to represent SOA formation from VCP emissions. The model incorporates a new VCP emissions inventory and employs three new classes of emissions: siloxanes, oxygenated IVOCs, and nonoxygenated IVOCs. VCPs are estimated to produce 1.67 µg m-3 of noontime SOA, doubling the current model predictions and reducing the SOA mass concentration bias from -75% to -58% when compared to observations in Los Angeles in 2010. While oxygenated and nonoxygenated intermediate volatility VCP species are emitted in similar quantities, SOA formation is dominated by the nonoxygenated IVOCs. Formaldehyde and SOA show similar relationships to temperature and bias signatures indicating common sources and/or chemistry. This work suggests that VCPs contribute up to half of anthropogenic SOA in Los Angeles and models must better represent SOA precursors from VCPs to predict the urban enhancement of SOA.

Evaluation of the offline-coupled GFSv15-FV3-CMAQv5.0.2 in support of the next-generation National Air Quality Forecast Capability over the contiguous United States.

As a candidate for the next-generation National Air Quality Forecast Capability (NAQFC), the meteorological forecast from the Global Forecast System with the new Finite Volume Cube-Sphere dynamical core (GFS-FV3) will be applied to drive the chemical evolution of gases and particles described by the Community Multiscale Air Quality modeling system. CMAQv5.0.2, a historical version of CMAQ, has been coupled with the North American Mesoscale Forecast System (NAM) model in the current operational NAQFC. An experimental version of the NAQFC based on the offline-coupled GFS-FV3 version 15 with CMAQv5.0.2 modeling system (GFSv15-CMAQv5.0.2) has been developed by the National Oceanic and Atmospheric Administration (NOAA) to provide real-time air quality forecasts over the contiguous United States (CONUS) since 2018. In this work, comprehensive region-specific, time-specific, and categorical evaluations are conducted for meteorological and chemical forecasts from the offline-coupled GFSv15-CMAQv5.0.2 for the year 2019. The forecast system shows good overall performance in forecasting meteorological variables with the annual mean biases of -0.2 °C for temperature at 2 m, 0.4% for relative humidity at 2 m, and 0.4 m s-1 for wind speed at 10 m compared to the METeorological Aerodrome Reports (METAR) dataset. Larger biases occur in seasonal and monthly mean forecasts, particularly in spring. Although the monthly accumulated precipitation forecasts show generally consistent spatial distributions with those from the remote-sensing and ensemble datasets, moderate-to-large biases exist in hourly precipitation forecasts compared to the Clean Air Status and Trends Network (CASTNET) and METAR. While the forecast system performs well in forecasting ozone (O3) throughout the year and fine particles with a diameter of 2.5 µm or less (PM2.5) for warm months (May-September), it significantly overpredicts annual mean concentrations of PM2.5. This is due mainly to the high predicted concentrations of fine fugitive and coarse-mode particle components. Underpredictions in the southeastern US and California during summer are attributed to missing sources and mechanisms of secondary organic aerosol formation from biogenic volatile organic compounds (VOCs) and semivolatile or intermediate-volatility organic compounds. This work demonstrates the ability of FV3-based GFS in driving the air quality forecasting. It identifies possible underlying causes for systematic region- and time-specific model biases, which will provide a scientific basis for further development of the next-generation NAQFC.

The Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module in the Community Multiscale Air Quality (CMAQ) modeling system version 5.3.2.

Air quality modeling for research and regulatory applications often involves executing many emissions sensitivity cases to quantify impacts of hypothetical scenarios, estimate source contributions, or quantify uncertainties. Despite the prevalence of this task, conventional approaches for perturbing emissions in chemical transport models like the Community Multiscale Air Quality (CMAQ) model require extensive offline creation and finalization of alternative emissions input files. This workflow is often time-consuming, error-prone, inconsistent among model users, difficult to document, and dependent on increased hard disk resources. The Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module, a component of CMAQv5.3 and beyond, addresses these limitations by performing these modifications online during the air quality simulation. Further, the model contains an Emission Control Interface which allows users to prescribe both simple and highly complex emissions scaling operations with control over individual or multiple chemical species, emissions sources, and spatial areas of interest. DESID further enhances the transparency of its operations with extensive error-checking and optional gridded output of processed emission fields. These new features are of high value to many air quality applications including routine perturbation studies, atmospheric chemistry research, and coupling with external models (e.g., energy system models, reduced-form models).

Predicting wildfire particulate matter and hypothetical re-emission of radiological Cs-137 contamination incidents.

Radiological release incidents can potentially contaminate widespread areas with radioactive materials and decontamination efforts are typically focused on populated areas, which means radionuclides may be left in forested areas for long periods of time. Large wildfires in contaminated forested areas have the potential to reintroduce these radionuclides into the atmosphere and cause exposure to first responders and downwind communities. One important radionuclide contaminant released from radiological incidents is radiocesium (137Cs) due to high yields and its long half-life of 30.2 years. An Eulerian 3D photochemical transport model was used to estimate potential ambient impacts of 137Cs re-emission due to wildfire following hypothetical radiological release scenarios. The Community Multiscale Air Quality (CMAQ) model did well at predicting levels and periods of increased PM2.5 carbon due to wildfire smoke at routine surface monitors in California during the summer of 2016. The model also did well at capturing the extent of the surface mixing layer compared to aerosol lidar measurements. Emissions from a large hypothetical wildfire were introduced into the wildland-urban interface (WUI) impacted by a hypothetical radiological release event. While ambient concentrations tended to be highest near the fire, the highest population committed effective dose equivalent by inhalation to an adult from 137Cs over an hour was downwind where wind flows moved smoke to high population areas. Seasonal variations in meteorology (wind flows) can result in differential population impacts even in the same metropolitan area. Modeled post-incident ambient levels of 137Cs both near these wildfires and further downwind in nearby urban areas were well below levels that would necessitate population evacuation or warrant other protective action recommendations such as shelter-in-place. These results suggest that 1) the modeling system captures local to regional scale transport and levels of PM2.5 from wildfire and 2) first responders and downwind population would not be expected to be at elevated risk from the initial inhalathion exposure of 137Cs re-emission.

The Community Multiscale Air Quality (CMAQ) model versions 5.3 and 5.3.1: system updates and evaluation.

The Community Multiscale Air Quality (CMAQ) model version 5.3 (CMAQ53), released to the public in August 2019 and followed by version 5.3.1 (CMAQ531) in December 2019, contains numerous science updates, enhanced functionality, and improved computation efficiency relative to the previous version of the model, 5.2.1 (CMAQ521). Major science advances in the new model include a new aerosol module (AERO7) with significant updates to secondary organic aerosol (SOA) chemistry, updated chlorine chemistry, updated detailed bromine and iodine chemistry, updated simple halogen chemistry, the addition of dimethyl sulfide (DMS) chemistry in the CB6r3 chemical mechanism, updated M3Dry bidirectional deposition model, and the new Surface Tiled Aerosol and Gaseous Exchange (STAGE) bidirectional deposition model. In addition, support for the Weather Research and Forecasting (WRF) model's hybrid vertical coordinate (HVC) was added to CMAQ53 and the Meteorology-Chemistry Interface Processor (MCIP) version 5.0 (MCIP50). Enhanced functionality in CMAQ53 includes the new Detailed Emissions Scaling, Isolation and Diagnostic (DESID) system for scaling incoming emissions to CMAQ and reading multiple gridded input emission files. Evaluation of CMAQ531 was performed by comparing monthly and seasonal mean daily 8 h average (MDA8) O3 and daily PM2.5 values from several CMAQ531 simulations to a similarly configured CMAQ521 simulation encompassing 2016. For MDA8 O3, CMAQ531 has higher O3 in the winter versus CMAQ521, due primarily to reduced dry deposition to snow, which strongly reduces wintertime O3 bias (2-4 ppbv monthly average). MDA8 O3 is lower with CMAQ531 throughout the rest of the year, particularly in spring, due in part to reduced O3 from the lateral boundary conditions (BCs), which generally increases MDA8 O3 bias in spring and fall ( 0.5 µg m-3). For daily 24 h average PM2.5, CMAQ531 has lower concentrations on average in spring and fall, higher concentrations in summer, and similar concentrations in winter to CMAQ521, which slightly increases bias in spring and fall and reduces bias in summer. Comparisons were also performed to isolate updates to several specific aspects of the modeling system, namely the lateral BCs, meteorology model version, and the deposition model used. Transitioning from a hemispheric CMAQ (HCMAQ) version 5.2.1 simulation to a HCMAQ version 5.3 simulation to provide lateral BCs contributes to higher O3 mixing ratios in the regional CMAQ simulation in higher latitudes during winter (due to the decreased O3 dry deposition to snow in CMAQ53) and lower O3 mixing ratios in middle and lower latitudes year-round (due to reduced O3 over the ocean with CMAQ53). Transitioning from WRF version 3.8 to WRF version 4.1.1 with the HVC resulted in consistently higher (1.0-1.5 ppbv) MDA8 O3 mixing ratios and higher PM2.5 concentrations (0.1-0.25 µg m-3) throughout the year. Finally, comparisons of the M3Dry and STAGE deposition models showed that MDA8 O3 is generally higher with M3Dry outside of summer, while PM2.5 is consistently higher with STAGE due to differences in the assumptions of particle deposition velocities to non-vegetated surfaces and land use with short vegetation (e.g., grasslands) between the two models. For ambient NH3, STAGE has slightly higher concentrations and smaller bias in the winter, spring, and fall, while M3Dry has higher concentrations and smaller bias but larger error and lower correlation in the summer.

Simulation of organic aerosol formation during the CalNex study: updated mobile emissions and secondary organic aerosol parameterization for intermediate-volatility organic compounds.

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 µgm-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.

Criteria pollutant impacts of volatile chemical products informed by near-field modeling.

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 (PM2.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 PM2.5 yield to 5% by mass, VCPs produce ~41% of the photochemical organic PM2.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.