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

Photochemical model assessment of single source NO2 and O3 plumes using field study data

Single source contribution to ambient O3 and PM2.5 has been estimated with photochemical grid models to support policy demonstrations for National Ambient Air Quality Standards, regional haze, and permit related programs. Limited field data exists to evaluate model representation of the spatial extent and chemical composition of plumes emitted by specific facilities. New tropospheric column measurements of NO2 and in-plume chemical measurements downwind of specific facilities allows for photochemical model evaluation of downwind plume extent, grid resolution impacts on plume concentration gradients, and source attribution methods. Here, photochemical models were applied with source sensitivity and source apportionment approaches to differentiate single source impacts on NO2 and O3 and compare with field study measurements. Source sensitivity approaches (e.g., brute-force difference method and decoupled direct method (DDM)) captured the spatial extent of NO2 plumes downwind of three facilities and the transition of near-source O3 titration to downwind production. Source apportionment approaches showed variability in terms of attributing the spatial extent of NO2 plumes and downwind O3 production. Each of the Community Multiscale Air Quality (CMAQ) source apportionment options predicted large O3 contribution from a large industrial facility in the flight transects nearest the facility when measurements and source sensitivity approaches suggest titration was outpacing production. In general, CMAQ DDM tends to attribute more O3 to boundary inflow and less to within-domain NOX and VOC sources compared to CMAQ source apportionment. The photochemical modeling system was able to capture single source plumes using 1 to 12 km grid resolution with best representation of plume extent and magnitude at the finer resolutions. When modeled at 1 to 12 km grid resolution, primary and secondary PM2.5 impacts were highest at the source location and decrease as distance increases downwind. The use of coarser grid resolution for single source attribution resulted in predicted impacts highest near the source but lower peak source specific concentrations compared to finer grid resolution simulations because impacts were spread out over a larger area.

Source-Receptor Relationships Between Precursor Emissions and O3 and PM2.5 Air Pollution Impacts

Reduced complexity tools that provide a representation of both primarily emitted particulate matter with an aerodynamic diameter less than 2.5 μm (PM2.5), secondarily formed PM2.5, and ozone (O3) allow for a quick assessment of many iterations of pollution control scenarios. Here, a new reduced complexity tool, Pattern Constructed Air Pollution Surfaces (PCAPS), that estimates annual average PM2.5 and seasonal average maximum daily average 8 h (MDA8) O3 for any source location in the United States is described and evaluated. Typically, reduced complexity tools are not evaluated for skill in predicting change in air pollution by comparison with more sophisticated modeling systems. Here, PCAPS was compared against multiple types of emission control scenarios predicted with state-of-the-science photochemical grid models to provide confidence that the model is realistically capturing the change in air pollution due to changing emissions. PCAPS was also applied with all anthropogenic emissions sources for multiple retrospective years to predict PM2.5 chemical components for comparison against routine surface measurements. PCAPS predicted similar magnitudes and regional variations in spatial gradients of measured chemical components of PM2.5. Model performance for capturing ambient measurements was consistent with other reduced complexity tools. PCAPS also did well at capturing the magnitude and spatial features of changes predicted by photochemical transport models for multiple emissions scenarios for both O3 and PM2.5. PCAPS is a flexible tool that provides source-receptor relationships using patterns of air quality gradients from a training data set of generic modeled sources to create interpolated air pollution gradients for new locations not part of the training database. The flexibility provided for both sources and receptors makes this tool ideal for integration into larger frameworks that provide emissions changes and need estimates of air quality to inform downstream analytics, which often includes an estimate of monetized health effects.

Evaluating reduced-form modeling tools for simulating ozone and PM2.5 monetized health impacts

Reduced-form modeling approaches are an increasingly popular way to rapidly estimate air quality and human health impacts related to changes in air pollutant emissions. These approaches reduce computation time by making simplifying assumptions about pollutant source characteristics, transport and chemistry. Two reduced form tools used by the Environmental Protection Agency in recent assessments are source apportionment-based benefit per ton (SA BPT) and source apportionment-based air quality surfaces (SABAQS). In this work, we apply these two reduced form tools to predict changes in ambient summer-season ozone, ambient annual PM2.5 component species and monetized health benefits for multiple sector-specific emission control scenarios: on-road mobile, electricity generating units (EGUs), cement kilns, petroleum refineries, and pulp and paper facilities. We then compare results against photochemical grid and standard health model-based estimates. We additionally compare monetized PM2.5 health benefits to values derived from three reduced form tools available in the literature: the Intervention Model for Air Pollution (InMAP), Air Pollution Emission Experiments and Policy Analysis (APEEP) version 2 (AP2) and Estimating Air pollution Social Impact Using Regression (EASIUR). Ozone and PM2.5 changes derived from SABAQS for EGU scenarios were well-correlated with values obtained from photochemical modeling simulations with spatial correlation coefficients between 0.64 and 0.89 for ozone and between 0.75 and 0.94 for PM2.5. SABAQS ambient ozone and PM2.5 bias when compared to photochemical modeling predictions varied by emissions scenario: SABAQS PM2.5 changes were overpredicted by up to 46% in one scenario and underpredicted by up to 19% in another scenario; SABAQS seasonal ozone changes were overpredicted by 34% to 83%. All tools predicted total PM2.5 benefits within a factor of 2 of the full-form predictions consistent with intercomparisons of reduced form tools available in the literature. As reduced form tools evolve, it is important to continue periodic comparison with comprehensive models to identify systematic biases in estimating air pollution impacts and resulting monetized health benefits.

An Analysis of CMAQ Gas Phase Dry Deposition over North America Through Grid-Scale and Land-Use Specific Diagnostics in the Context of AQMEII4

The fourth phase of the Air Quality Model Evaluation International Initiative (AQMEII4) is conducting a diagnostic intercomparison and evaluation of deposition simulated by regional-scale air quality models over North America and Europe. In this study, we analyze annual AQMEII4 simulations performed with the Community Multiscale Air Quality Model (CMAQ) version 5.3.1 over North America. These simulations were configured with both the M3Dry and Surface Tiled Aerosol and Gas Exchange (STAGE) dry deposition schemes available in CMAQ. A comparison of observed and modeled concentrations and wet deposition fluxes shows that the AQMEII4 CMAQ simulations perform similarly to other contemporary regional-scale modeling studies. During summer, M3Dry has higher ozone (O3) deposition velocities (Vd) and lower mixing ratios than STAGE for much of the eastern U.S. while the reverse is the case over eastern Canada and along the West Coast. In contrast, during winter STAGE has higher O3 Vd and lower mixing ratios than M3Dry over most of the southern half of the modeling domain while the reverse is the case for much of the northern U.S. and southern Canada. Analysis of the diagnostic variables defined for the AQMEII4 project, i.e. grid-scale and land-use (LU) specific effective conductances and deposition fluxes for the major dry deposition pathways, reveals generally higher summertime stomatal and wintertime cuticular grid-scale effective conductances for M3Dry and generally higher soil grid-scale effective conductances (for both vegetated and bare soil) for STAGE in both summer and winter. On a domain-wide basis, the stomatal grid-scale effective conductances account for about half of the total O3 Vd during daytime hours in summer for both schemes. Employing LU-specific diagnostics, results show that daytime Vd varies by a factor of 2 between LU categories. Furthermore, M3Dry vs. STAGE differences are most pronounced for the stomatal and vegetated soil pathway for the forest LU categories, with M3Dry estimating larger effective conductances for the stomatal pathway and STAGE estimating larger effective conductances for the vegetated soil pathway for these LU categories. Annual domain total O3 deposition fluxes differ only slightly between M3Dry (74.4 Tg/year) and STAGE (76.2 Tg/yr), but pathway-specific fluxes to individual LU types can vary more substantially on both annual and seasonal scales which would affect estimates of O3 damages to sensitive vegetation. A comparison of two simulations differing only in their LU classification scheme shows that the differences in LU cause seasonal mean O3 mixing ratio differences on the order of 1 ppb across large portions of the domain, with the differences generally largest during summer and in areas characterized by the largest differences in the fractional coverages of the forest, planted/cultivated, and grassland LU categories. These differences are generally smaller than the M3Dry vs. STAGE differences outside the summer season but have a similar magnitude during summer. Results indicate that the deposition impacts of LU differences are caused both by differences in the fractional coverages and spatial distributions of different LU categories as well as the characterization of these categories through variables like surface roughness and vegetation fraction in look-up tables used in the land-surface model and deposition schemes. Overall, the analyses and results presented in this study illustrate how the diagnostic grid-scale and LU-specific dry deposition variables adopted for AQMEII4 can provide insights into similarities and differences between the CMAQ M3Dry and STAGE dry deposition schemes that affect simulated pollutant budgets and ecosystem impacts from atmospheric pollution.

Characterizing variations in ambient PM2.5 concentrations at the U.S. Embassy in Dhaka, Bangladesh using observations and the CMAQ modeling system

We analyze hourly PM2.5 (particles with an aerodynamic diameter of ≤ 2.5 μm) concentrations measured at the U.S. Embassy in Dhaka over the 2016 - 2021 time period and find that concentrations are seasonally dependent with the highest occurring in winter and the lowest in monsoon seasons. Mean winter PM2.5 concentrations reached ~165-175 μg/m3 while monsoon concentrations remained ~30-35 μg/m3. Annual mean PM2.5 concentration reached ~5-6 times greater than the Bangladesh annual PM2.5 standard of 15 μg/m3. The number of days exceeding the daily PM2.5 standard of 65 μg/m3 in a year approached nearly 50%. Daily-mean PM2.5 concentrations remained elevated (>65 μg/m3) for more than 80 consecutive days. Night-time concentrations were greater than daytime concentrations. The comparison of results obtained from the Community Multiscale Air Quality (CMAQ) model simulations over the Northern Hemisphere using 108-km horizontal grids with observed data suggests that the model can reproduce the seasonal variation of observed data but underpredicts observed PM2.5 in winter months with a normalized mean bias of 13-32%. In the model, organic aerosol is the largest component of PM2.5, of which secondary organic aerosol plays a dominant role. Transboundary pollution has a large impact on the PM2.5 concentration in Dhaka, with an annual mean contribution of ~40 μg/m3.

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 PM_2.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.

Contribution of Organic Nitrates to Organic Aerosol over South Korea during KORUS-AQ

The role of anthropogenic NOx 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 (RONO2) in SOA formation over the Korean Peninsula during the Korea-United States Air Quality field study in Spring 2016 as a model for RONO2 aerosol in cities worldwide. We use aircraft-based measurements of the particle phase and total (gas + particle) RONO2 to explore RONO2 phase partitioning. These measurements show that, on average, one-fourth of RONO2 are in the condensed phase, and we estimate that ≈15% of the organic aerosol (OA) mass can be attributed to RONO2. Furthermore, we observe that the fraction of RONO2 in the condensed phase increases with OA concentration, evidencing that equilibrium absorptive partitioning controls the RONO2 phase distribution. Lastly, we model RONO2 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 RONO2, but there is a large missing source of semivolatile, anthropogenically derived RONO2. 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.

A comparative study of two-way and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US: performance evaluation and impacts of chemistry-meteorology feedbacks on air quality

The two-way coupled Weather Research and Forecasting and Community Multiscale Air Quality (WRF-CMAQ) model has been developed to more realistically represent the atmosphere by accounting for complex chemistry-meteorology feedbacks. In this study, we present a comparative analysis of two-way (with consideration of both aerosol direct and indirect effects) and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US. Long-term (5 years from 2008 to 2012) simulations using WRF-CMAQ with both offline and two-way coupling modes are carried out with anthropogenic emissions based on multiple years of the U.S. National Emission Inventory and chemical initial and boundary conditions derived from an advanced Earth system model (i.e., a modified version of the Community Earth System Model/Community Atmospheric Model). The comprehensive model evaluations show that both two-way WRF-CMAQ and WRF-only simulations perform well for major meteorological variables such as temperature at 2 m, relative humidity at 2 m, wind speed at 10 m, precipitation (except for against the National Climatic Data Center data), and shortwave and longwave radiation. Both two-way and offline CMAQ also show good performance for ozone (O₃) and fine particulate matter (PM_2.5). Due to the consideration of aerosol direct and indirect effects, two-way WRF-CMAQ shows improved performance over offline coupled WRF and CMAQ in terms of spatiotemporal distributions and statistics, especially for radiation, cloud forcing, O₃, sulfate, nitrate, ammonium, elemental carbon, tropospheric O₃ residual, and column nitrogen dioxide (NO₂). For example, the mean biases have been reduced by more than 10 W m^-2 for shortwave radiation and cloud radiative forcing and by more than 2 ppb for max 8 h O₃. However, relatively large biases still exist for cloud predictions, some PM_2.5 species, and PM₁₀ that warrant follow-up studies to better understand those issues. The impacts of chemistry-meteorological feedbacks are found to play important roles in affecting regional air quality in the US by reducing domain-average concentrations of carbon monoxide (CO), O₃, nitrogen oxide (NO _x ), volatile organic compounds (VOCs), and PM_2.5 by 3.1% (up to 27.8%), 4.2% (up to 16.2%), 6.6% (up to 50.9%), 5.8% (up to 46.6%), and 8.6% (up to 49.1%), respectively, mainly due to reduced radiation, temperature, and wind speed. The overall performance of the two-way coupled WRF-CMAQ model achieved in this work is generally good or satisfactory and the improved performance for two-way coupled WRF-CMAQ should be considered along with other factors in developing future model applications to inform policy making.

A comparative study of two-way and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US: performance evaluation and impacts of chemistry-meteorology feedbacks on air quality

The two-way coupled Weather Research and Forecasting and Community Multiscale Air Quality (WRF-CMAQ) model has been developed to more realistically represent the atmosphere by accounting for complex chemistry-meteorology feedbacks. In this study, we present a comparative analysis of two-way (with consideration of both aerosol direct and indirect effects) and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US. Long-term (5 years from 2008 to 2012) simulations using WRF-CMAQ with both offline and two-way coupling modes are carried out with anthropogenic emissions based on multiple years of the U.S. National Emission Inventory and chemical initial and boundary conditions derived from an advanced Earth system model (i.e., a modified version of the Community Earth System Model/Community Atmospheric Model). The comprehensive model evaluations show that both two-way WRF-CMAQ and WRF-only simulations perform well for major meteorological variables such as temperature at 2 m, relative humidity at 2 m, wind speed at 10 m, precipitation (except for against the National Climatic Data Center data), and shortwave and longwave radiation. Both two-way and offline CMAQ also show good performance for ozone (O₃) and fine particulate matter (PM_2.5). Due to the consideration of aerosol direct and indirect effects, two-way WRF-CMAQ shows improved performance over offline coupled WRF and CMAQ in terms of spatiotemporal distributions and statistics, especially for radiation, cloud forcing, O₃, sulfate, nitrate, ammonium, elemental carbon, tropospheric O₃ residual, and column nitrogen dioxide (NO₂). For example, the mean biases have been reduced by more than 10 W m^-2 for shortwave radiation and cloud radiative forcing and by more than 2 ppb for max 8 h O₃. However, relatively large biases still exist for cloud predictions, some PM_2.5 species, and PM₁₀ that warrant follow-up studies to better understand those issues. The impacts of chemistry-meteorological feedbacks are found to play important roles in affecting regional air quality in the US by reducing domain-average concentrations of carbon monoxide (CO), O₃, nitrogen oxide (NO _x ), volatile organic compounds (VOCs), and PM_2.5 by 3.1% (up to 27.8%), 4.2% (up to 16.2%), 6.6% (up to 50.9%), 5.8% (up to 46.6%), and 8.6% (up to 49.1%), respectively, mainly due to reduced radiation, temperature, and wind speed. The overall performance of the two-way coupled WRF-CMAQ model achieved in this work is generally good or satisfactory and the improved performance for two-way coupled WRF-CMAQ should be considered along with other factors in developing future model applications to inform policy making.

Confronting Uncertainties of Simulated Air Pollution Concentrations during Persistent Cold Air Pool Events in the Salt Lake Valley, Utah

Air pollutant accumulations during wintertime persistent cold air pool (PCAP) events in mountain valleys are of great concern for public health worldwide. Uncertainties associated with the simulated meteorology under stable conditions over complex terrain hinder realistic simulations of air quality using chemical transport models. We use the Community Multiscale Air Quality (CMAQ) model to simulate the gaseous and particulate species for 1 month in January 2011 during the Persistent Cold Air Pool Study (PCAPS) in the Salt Lake Valley (SLV), Utah (USA). Results indicate that the temporal variability associated with the elevated NO_x and PM_2.5 concentrations during PCAP events was captured by the model (r = 0.20 for NO_x and r = 0.49 for PM_2.5). However, concentrations were not at the correct magnitude (NMB = -35/12% for PM_2.5 during PCAPs/non-PCAPs), where PM_2.5 was underestimated during PCAP events and overestimated during non-PCAP periods. The underestimated PCAP strength is represented by valley heat deficit, which contributed to the underestimated PM_2.5 concentrations compared with observations due to the model simulating more vertical mixing and less stable stratification than what was observed. Based on the observations, the dominant PM_2.5 species were ammonium and nitrate. We provide a discussion that aims to investigate the emissions and chemistry model uncertainties using the nitrogen ratio method and the thermodynamic ammonium nitrate regime method.

The Chesapeake Bay Program Modeling System: Overview and Recommendations for Future Development

The Chesapeake Bay is the largest, most productive, and most biologically diverse estuary in the continental United States providing crucial habitat and natural resources for culturally and economically important species. Pressures from human population growth and associated development and agricultural intensification have led to excessive nutrient and sediment inputs entering the Bay, negatively affecting the health of the Bay ecosystem and the economic services it provides. The Chesapeake Bay Program (CBP) is a unique program formally created in 1983 as a multi-stakeholder partnership to guide and foster restoration of the Chesapeake Bay and its watershed. Since its inception, the CBP Partnership has been developing, updating, and applying a complex linked modeling system of watershed, airshed, and estuary models as a planning tool to inform strategic management decisions and Bay restoration efforts. This paper provides a description of the 2017 CBP Modeling System and the higher trophic level models developed by the NOAA Chesapeake Bay Office, along with specific recommendations that emerged from a 2018 workshop designed to inform future model development. Recommendations highlight the need for simulation of watershed inputs, conditions, processes, and practices at higher resolution to provide improved information to guide local nutrient and sediment management plans. More explicit and extensive modeling of connectivity between watershed landforms and estuary sub-areas, estuarine hydrodynamics, watershed and estuarine water quality, the estuarine-watershed socioecological system, and living resources will be important to broaden and improve characterization of responses to targeted nutrient and sediment load reductions. Finally, the value and importance of maintaining effective collaborations among jurisdictional managers, scientists, modelers, support staff, and stakeholder communities is emphasized. An open collaborative and transparent process has been a key element of successes to date and is vitally important as the CBP Partnership moves forward with modeling system improvements that help stakeholders evolve new knowledge, improve management strategies, and better communicate outcomes.

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.

Coupled Air Quality and Boundary-Layer Meteorology in Western U.S. Basins during Winter: Design and Rationale for a Comprehensive Study

Wintertime episodes of high aerosol concentrations occur frequently in urban and agricultural basins and valleys worldwide. These episodes often arise following development of persistent cold-air pools (PCAPs) that limit mixing and modify chemistry. While field campaigns targeting either basin meteorology or wintertime pollution chemistry have been conducted, coupling between interconnected chemical and meteorological processes remains an insufficiently studied research area. Gaps in understanding the coupled chemical-meteorological interactions that drive high pollution events make identification of the most effective air-basin specific emission control strategies challenging. To address this, a September 2019 workshop occurred with the goal of planning a future research campaign to investigate air quality in Western U.S. basins. Approximately 120 people participated, representing 50 institutions and 5 countries. Workshop participants outlined the rationale and design for a comprehensive wintertime study that would couple atmospheric chemistry and boundary-layer and complex-terrain meteorology within western U.S. basins. Participants concluded the study should focus on two regions with contrasting aerosol chemistry: three populated valleys within Utah (Salt Lake, Utah, and Cache Valleys) and the San Joaquin Valley in California. This paper describes the scientific rationale for a campaign that will acquire chemical and meteorological datasets using airborne platforms with extensive range, coupled to surface-based measurements focusing on sampling within the near-surface boundary layer, and transport and mixing processes within this layer, with high vertical resolution at a number of representative sites. No prior wintertime basin-focused campaign has provided the breadth of observations necessary to characterize the meteorological-chemical linkages outlined here, nor to validate complex processes within coupled atmosphere-chemistry models.

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).

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).

Furthering a partnership: Air quality modeling and improving public health

Air pollution is one of the top five risk factors for population health globally. In recent years, advances in air pollution data and modeling have occurred simultaneously with advances in data and methods available for health studies. To realize the potential of such advances, the air quality modeling and public health communities should continue to strengthen their engagements and build effective interdisciplinary teams. These partnerships recognize the tight coupling between air quality and health data and methods and the value of expertise from multiple fields to ensure that this information is applied appropriately with a deep understanding of its capabilities and limitations. Building effective multidisciplinary teams takes a sustained commitment to engage with partners with different expertise to establish working partnerships and collaborations to better address public exposures to air pollution. Effective partnerships enable better targeting of research resources to answer important questions and provide essential information to protect public health.Implications: Air quality models are an effective tool that can be used to estimate air pollution exposure in epidemiologic studies and risk assessments. Working together in collaborative multidisciplinary teams will lead to greater advancements in understanding of air pollution impacts and in useful information informing actions to improve public health.

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) O₃ and daily PM_2.5 values from several CMAQ531 simulations to a similarly configured CMAQ521 simulation encompassing 2016. For MDA8 O₃, CMAQ531 has higher O₃ in the winter versus CMAQ521, due primarily to reduced dry deposition to snow, which strongly reduces wintertime O₃ bias (2-4 ppbv monthly average). MDA8 O₃ is lower with CMAQ531 throughout the rest of the year, particularly in spring, due in part to reduced O₃ from the lateral boundary conditions (BCs), which generally increases MDA8 O₃ bias in spring and fall ( 0.5 μg m^-3). For daily 24 h average PM_2.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 O₃ mixing ratios in the regional CMAQ simulation in higher latitudes during winter (due to the decreased O₃ dry deposition to snow in CMAQ53) and lower O₃ mixing ratios in middle and lower latitudes year-round (due to reduced O₃ 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 O₃ mixing ratios and higher PM_2.5 concentrations (0.1-0.25 μg m^-3) throughout the year. Finally, comparisons of the M3Dry and STAGE deposition models showed that MDA8 O₃ is generally higher with M3Dry outside of summer, while PM_2.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 NH₃, 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.

Evaluation of 15 years of modeled atmospheric oxidized nitrogen compounds across the contiguous United States

Atmospheric nitrogen oxide and nitrogen dioxide (NO + NO2, together termed as NO X ) estimates from annual photochemical simulations for years 2002-2016 are compared to surface network measurements of NO X and total gas-phase-oxidized reactive nitrogen (NO Y ) to evaluate the Community Multiscale Air Quality (CMAQ) modeling system performance by U.S. region, season, and time of day. In addition, aircraft measurements from 2011 Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality are used to evaluate how emissions, chemical mechanism, and measurement uncertainty each contribute to the overall model performance. We show distinct seasonal and time-of-day patterns in NO X performance. Summertime NO X is overpredicted with bimodal peaks in bias during early morning and evening hours and persisting overnight. The summertime morning NO X bias dropped from between 28% and 57% for earlier years (2002-2012) to between -2% and 7% for later years (2013-2016). Summer daytime NO X tends to be unbiased or underpredicted. In winter, the evening NO X overpredictions remain, but NO X is unbiased or underpredicted overnight, in the morning, and during the day. NO X overpredictions are most pronounced in the Midwestern and Southern United States with Western regions having more of a tendency toward model underpredictions of NO X . Modeled NO X performance has improved substantially over time, reflecting updates to the emission inputs and the CMAQ air quality model. Model performance improvements are largest for years simulated with CMAQv5.1 or later and for emission inventory years 2014 and later, coinciding with reduced onroad NO X emissions from vehicles with newer emission control technologies and improved treatment of chemistry, deposition, and vertical mixing in CMAQ. Our findings suggest that emissions temporalization of specific mobile source sectors have a small impact on model performance, while chemistry updates improve predictions of NO Y but do not improve summertime NO X bias in the Baltimore/DC area. Sensitivity runs performed for different locations across the country suggest that the improvement in summer NO X performance can be attributed to updates in vertical mixing incorporated in CMAQv5.1.

Examining PM2.5 concentrations and exposure using multiple models

Epidemiologic studies have found associations between fine particulate matter (PM_2.5) exposure and adverse health effects using exposure models that incorporate monitoring data and other relevant information. Here, we use nine PM_2.5 concentration models (i.e., exposure models) that span a wide range of methods to investigate i) PM_2.5 concentrations in 2011, ii) potential changes in PM_2.5 concentrations between 2011 and 2028 due to on-the-books regulations, and iii) PM_2.5 exposure for the U.S. population and four racial/ethnic groups. The exposure models included two geophysical chemical transport models (CTMs), two interpolation methods, a satellite-derived aerosol optical depth-based method, a Bayesian statistical regression model, and three data-rich machine learning methods. We focused on annual predictions that were regridded to 12-km resolution over the conterminous U.S., but also considered 1-km predictions in sensitivity analyses. The exposure models predicted broadly consistent PM_2.5 concentrations, with relatively high concentrations on average over the eastern U.S. and greater variability in the western U.S. However, differences in national concentration distributions (median standard deviation: 1.00 μg m^-3) and spatial distributions over urban areas were evident. Further exploration of these differences and their implications for specific applications would be valuable. PM_2.5 concentrations were estimated to decrease by about 1 μg m^-3 on average due to modeled emission changes between 2011 and 2028, with decreases of more than 3 μg m^-3 in areas with relatively high 2011 concentrations that were projected to experience relatively large emission reductions. Agreement among models was closer for population-weighted than uniformly weighted averages across the domain. About 50% of the population was estimated to experience PM_2.5 concentrations less than 10 μg m^-3 in 2011 and PM_2.5 improvements of about 2 μg m^-3 due to modeled emission changes between 2011 and 2028. Two inequality metrics were used to characterize differences in exposure among the four racial/ethnic groups. The metrics generally yielded consistent information and suggest that the modeled emission reductions between 2011 and 2028 would reduce absolute exposure inequality on average.

Improving Estimates of PM2.5 Concentration and Chemical Composition by Application of High Spectral Resolution Lidar (HSRL) and Creating Aerosol Types from Chemistry (CATCH) Algorithm

Improved characterization of ambient PM2.5 mass concentration and chemical speciation is a topic of interest in air quality and climate sciences. Over the past decades, considerable efforts have been made to improve ground-level PM2.5 using remotely sensed data. Here we present two new approaches for estimating atmospheric PM2.5 and chemical composition based on the High Spectral Resolution Lidar (HSRL)-retrieved aerosol extinction values and types and Creating Aerosol Types from Chemistry (CATCH)-derived aerosol chemical composition. The first methodology (CMAQ-HSRL-CH) improves EPA's Community Multiscale Air Quality (CMAQ) predictions by applying variable scaling factors derived using remotely-sensed information about aerosol vertical distribution and types and the CATCH algorithm. The second methodology (HSRL-CH) does not require regional model runs and can provide atmospheric PM2.5 mass concentration and chemical speciation using only the remotely sensed data and the CATCH algorithm. The resulting PM2.5 concentrations and chemical speciation derived for NASA DISCOVER-AQ (Deriving Information on Surface Conditions from COlumn and VERtically Resolved Observations Relevant to Air Quality) Baltimore-Washington, D.C. Corridor (BWC) Campaign (2011) are compared to surface measurements from EPA's Air Quality Systems (AQS) network. The analysis shows that the CMAQ-HSRL-CH method leads to considerable improvement of CMAQ's predicted PM2.5 concentrations (R2 value increased from 0.37 to 0.63, the root mean square error (RMSE) was reduced from 11.9 to 7.2 μg m-3, and the normalized mean bias (NMB) was lowered from -46.0 to 4.6%). The HSRL-CH method showed statistics (R2=0.75, RMSE=8.6 μgm-3, and NMB=24.0%), which were better than the CMAQ prediction of PM2.5 alone and analogous to CMAQ-HSRL-CH. In addition to mass concentration, HSRL-CH can also provide aerosol chemical composition without specific model simulations. We expect that the HSRL-CH method will be able to make reliable estimates of PM2.5 concentration and chemical composition where HSRL data are available.

Monthly Patterns of Ammonia Over the Contiguous United States at 2-km Resolution

Monthly, high-resolution (∼2 km) ammonia (NH3) column maps from the Infrared Atmospheric Sounding Interferometer (IASI) were developed across the contiguous United States and adjacent areas. Ammonia hotspots (95th percentile of the column distribution) were highly localized with a characteristic length scale of 12 km and median area of 152 km2. Five seasonality clusters were identified with k-means++ clustering. The Midwest and eastern United States had a broad, spring maximum of NH3 (67% of hotspots in this cluster). The western United States, in contrast, showed a narrower midsummer peak (32% of hotspots). IASI spatiotemporal clustering was consistent with those from the Ammonia Monitoring Network. CMAQ and GFDL-AM3 modeled NH3 columns have some success replicating the seasonal patterns but did not capture the regional differences. The high spatial-resolution monthly NH3 maps serve as a constraint for model simulations and as a guide for the placement of future, ground-based network sites.

Evaluating trends and seasonality in modeled PM2.5 concentrations using empirical mode decomposition

Regional-scale air quality models are being used for studying the sources, composition, transport, transformation, and deposition of fine particulate matter (PM_2.5). The availability of decadal air quality simulations provides a unique opportunity to explore sophisticated model evaluation techniques rather than relying solely on traditional operational evaluations. In this study, we propose a new approach for process-based model evaluation of speciated PM_2.5 using improved complete ensemble empirical mode decomposition with adaptive noise (improved CEEMDAN) to assess how well version 5.0.2 of the coupled Weather Research and Forecasting model-Community Multiscale Air Quality model (WRF-CMAQ) simulates the time-dependent long-term trend and cyclical variations in daily average PM_2.5 and its species, including sulfate (SO₄), nitrate (NO₃), ammonium (NH₄), chloride (Cl), organic carbon (OC), and elemental carbon (EC). The utility of the proposed approach for model evaluation is demonstrated using PM_2.5 data at three monitoring locations. At these locations, the model is generally more capable of simulating the rate of change in the long-term trend component than its absolute magnitude. Amplitudes of the sub-seasonal and annual cycles of total PM_2.5, SO₄, and OC are well reproduced. However, the time-dependent phase difference in the annual cycles for total PM_2.5, OC, and EC reveals a phase shift of up to half a year, indicating the need for proper temporal allocation of emissions and for updating the treatment of organic aerosols compared to the model version used for this set of simulations. Evaluation of sub-seasonal and interannual variations indicates that CMAQ is more capable of replicating the sub-seasonal cycles than interannual variations in magnitude and phase.

Meteorological and Air Quality Modeling for Hawaii, Puerto Rico, and Virgin Islands

A photochemical model platform for Hawaii, Puerto Rico, and Virgin Islands predicting O₃, PM_2.5, and regional haze would be useful to support assessments relevant for the National Ambient Air Quality Standards (NAAQS), Regional Haze Rule, and the Prevention of Significant Deterioration (PSD) program. These areas have not traditionally been modeled with photochemical transport models, but a reasonable representation of meteorology, emissions (natural and anthropogenic), chemistry, and deposition could support air quality management decisions in these areas. Here, a prognostic meteorological model (Weather Research and Forecasting) and photochemical transport (Community Multiscale Air Quality) model were applied for the entire year of 2016 at 27, 9, and 3 km grid resolution for areas covering the Hawaiian Islands and Puerto Rico/Virgin Islands. Model predictions were compared against surface and upper air meteorological and chemical measurements available in both areas. The vertical gradient of temperature, humidity, and winds in the troposphere was well represented. Surface layer meteorological model performance was spatially variable, but temperature tended to be underestimated in Hawaii. Chemically speciated daily average PM_2.5 was generally well characterized by the modeling system at urban and rural monitors in Hawaii and Puerto Rico/Virgin Islands. Model performance was notably impacted by the wildfire emission methodology. Model performance was mixed for hourly SO₂, NO₂, PM_2.5, and CO and was often related to how well local emissions sources were characterized. SO₂ predictions were much lower than measurements at monitors near active volcanos on Hawaii, which was expected since volcanic emissions were not included in these model simulations. Further research is needed to assess emission inventory representation of these areas and how microscale meteorology influenced by the complex land-water and terrain interfaces impacts higher time resolution performance.