A framework for expanding aqueous chemistry in the Community Multiscale Air Quality (CMAQ) model version 5.1

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.

Nonagricultural emissions enhance dimethylamine and modulate urban atmospheric nucleation

Gas-phase dimethylamine (DMA) has recently been identified as one of the most important vapors to initiate new particle formation (NPF), even in China's polluted atmosphere. Nevertheless, there remains a fundamental need for understanding the atmospheric life cycle of DMA, particularly in urban areas. Here we pioneered large-scale mobile observations of the DMA concentrations within cities and across two pan-region transects of north-to-south (∼700 km) and west-to-east (∼2000 km) in China. Unexpectedly, DMA concentrations (mean ± 1σ) in South China with scattered croplands (0.018 ± 0.010 ppbv) were over three times higher than those in the north with contiguous croplands (0.005 ± 0.001 ppbv), suggesting that nonagricultural activities may be an important source of DMA. Particularly in non-rural regions, incidental pulsed industrial emissions led to some of the highest DMA concentration levels in the world (>2.3 ppbv). Besides, in highly urbanized areas of Shanghai, supported by direct source-emission measurements, the spatial pattern of DMA was generally correlated with population (R² = 0.31) due to associated residential emissions rather than vehicular emissions. Chemical transport simulations further show that in the most populated regions of Shanghai, residential DMA emissions can contribute for up to 78% of particle number concentrations. Shanghai is a case study for populous megacities, and the impacts of nonagricultural emissions on local DMA concentration and nucleation are likely similar for other major urban regions globally.

Extreme Emission Reduction Requirements for China to Achieve World Health Organization Global Air Quality Guidelines

A big gap exists between current air quality in China and the World Health Organization (WHO) global air quality guidelines (AQG) released in 2021. Previous studies on air pollution control have focused on emission reduction demand in China but ignored the influence of transboundary pollution, which has been proven to have a significant impact on air quality in China. Here, we develop an emission-concentration response surface model coupled with transboundary pollution to quantify the emission reduction demand for China to achieve WHO AQG. China cannot achieve WHO AQG by its own emission reduction for high transboundary pollution of both PM_2.5 and O₃. Reducing transboundary pollution will loosen the reduction demand for NH₃ and VOCs emissions in China. However, to meet 10 μg·m^-3 for PM_2.5 and 60 μg·m^-3 for peak season O₃, China still needs to reduce its emissions of SO₂, NO_x, NH₃, VOCs, and primary PM_2.5 by more than 95, 95, 76, 62, and 96% respectively, on the basis of 2015. We highlight that both extreme emission reduction in China and great efforts in addressing transboundary air pollution are crucial to reach WHO AQG.

Pronounced increases in nitrogen emissions and deposition due to the historic 2020 wildfires in the western U.S

Wildfire outbreaks can lead to extreme biomass burning (BB) emissions of both oxidized (e.g., nitrogen oxides; NO_x = NO+NO₂) and reduced form (e.g., ammonia; NH₃) nitrogen (N) compounds. High N emissions are major concerns for air quality, atmospheric deposition, and consequential human and ecosystem health impacts. In this study, we use both satellite-based observations and modeling results to quantify the contribution of BB to the total emissions, and approximate the impact on total N deposition in the western U.S. Our results show that during the 2020 wildfire season of August-October, BB contributes significantly to the total emissions, with a satellite-derived fraction of NH₃ to the total reactive N emissions (median ~ 40%) in the range of aircraft observations. During the peak of the western August Complex Fires in September, BB contributed to ~55% (for the contiguous U.S.) and ~ 83% (for the western U.S.) of the monthly total NO_x and NH₃ emissions. Overall, there is good model performance of the George Mason University-Wildfire Forecasting System (GMU-WFS) used in this work. The extreme BB emissions lead to significant contributions to the total N deposition for different ecosystems in California, with an average August - October 2020 relative increase of ~78% (from 7.1 to 12.6 kg ha^-1 year^-1) in deposition rate to major vegetation types (mixed forests + grasslands/shrublands/savanna) compared to the GMU-WFS simulations without BB emissions. For mixed forest types only, the average N deposition rate increases (from 6.2 to 16.9 kg ha^-1 year^-1) are even larger at ~173%. Such large N deposition due to extreme BB emissions are much (~6-12 times) larger than low-end critical load thresholds for major vegetation types (e.g., forests at 1.5-3 kg ha^-1 year^-1), and thus may result in adverse N deposition effects across larger areas of lichen communities found in California's mixed conifer forests.

Using wildland fire smoke modeling data in gerontological health research (California, 2007-2018)

Widespread population exposure to wildland fire smoke underscores the urgent need for new techniques to characterize fire-derived pollution for epidemiologic studies and to build climate-resilient communities especially for aging populations. Using atmospheric chemical transport modeling, we examined air quality with and without wildland fire smoke PM_2.5. In 12-km gridded output, the 24-hour average concentration of all-source PM_2.5 in California (2007-2018) was 5.16 μg/m³ (S.D. 4.66 μg/m³). The average concentration of fire-PM_2.5 in California by year was 1.61 μg/m³ (~30% of total PM_2.5). The contribution of fire-source PM_2.5 ranged from 6.8% to 49%. We define a "smokewave" as two or more consecutive days with modeled levels above 35 μg/m³. Based on model-derived fire-PM_2.5, 99.5% of California's population lived in a county that experienced at least one smokewave from 2007 to 2018, yet understanding of the impact of smoke on the health of aging populations is limited. Approximately 2.7 million (56%) of California residents aged 65+ years lived in counties representing the top 3 quartiles of fire-PM_2.5 concentrations (2007-2018). For each year (2007-2018), grid cells containing skilled nursing facilities had significantly higher mean concentrations of all-source PM_2.5 than cells without those facilities, but they also had generally lower mean concentrations of wildland fire-specific PM_2.5. Compared to rural monitors in California, model predictions of wildland fire impacts on daily average PM_2.5 carbon (organic and elemental) performed well most years but tended to overestimate wildland fire impacts for high-fire years. The modeling system isolated wildland fire PM_2.5 from other sources at monitored and unmonitored locations, which is important for understanding exposures for aging population in health studies.

Projections of future wildfires impacts on air pollutants and air toxics in a changing climate over the western United States

Wildfires emit smoke particles and gaseous pollutants that greatly aggravate air quality and cause adverse health impacts in the western US (WUS). This study evaluates how wildfire impacts on air pollutants and air toxics evolve from the present climate to the future climate under a high anthropogenic emission scenario at regional and city scales. Through employing multiple climate and chemical transport models, small changes in domain-averaged air pollutant concentrations by wildfires are simulated over WUS. However, such changes significantly increase future city-scale pollutant concentrations by up to 53 ppb for benzene, 158 ppb for formaldehyde, 655 μg/m³ for fine particulate matter (PM_2.5), and 102 ppb for ozone, whereas that for the present climate are 104 ppb for benzene, 332 ppb for formaldehyde, 1,378 μg/m³ for PM_2.5, and 140 ppb for ozone. Despite wildfires induce smaller changes in the future, the wildfire contribution ratios can increase by more than tenfold compared to the present climate, indicating wildfires become a more critical contributor to future air pollution in WUS. In addition, additional 6 exceedance days/year for formaldehyde and additional 3 exceedance days/year for ozone suggest increasing health impacts by wildfires in the future.

Nitrogen burden from atmospheric deposition in East Asian oceans in 2010 based on high-resolution regional numerical modeling

East Asian oceans are possibly affected by a high nitrogen (N) burden because of the intense anthropogenic emissions in this region. Based on high-resolution regional chemical transport modeling with horizontal grid scales of 36 and 12 km, we investigated the N burden into East Asian oceans via atmospheric deposition in 2010. We found a high N burden of 2-9 kg N ha^-1 yr^-1 over the Yellow Sea, East China Sea (ECS), and Sea of Japan. Emissions over East Asia were dominated by ammonia (NH₃) over land and nitrogen oxides (NO_x) over oceans, and N deposition was dominated by reduced N over most land and open ocean, whereas it was dominated by oxidized N over marginal seas and desert areas. The verified numerical modeling identified that the following processes were quantitatively important over East Asian oceans: the dry deposition of nitric acid (HNO₃), NH₃, and coarse-mode (aerodynamic diameter greater than 2.5 μm) NO₃^-, and wet deposition of fine-mode (aerodynamic diameter less than 2.5 μm) NO₃^- and NH₄⁺. The relative importance of the dry deposition of coarse-mode NO₃^- was higher over open ocean. The estimated N deposition to the whole ECS was 390 Gg N yr^-1; this is comparable to the discharge from the Yangtze River to the ECS, indicating the significant contribution of atmospheric deposition. Based on the high-resolution modeling over the ECS, a tendency of high deposition in the western ECS and low deposition in the eastern ECS was found, and a variety of deposition processes were estimated. The dry deposition of coarse-mode NO₃^- and wet deposition of fine-mode NH₄⁺ were the main factors, and the wet deposition of fine-mode NO₃^- over the northeastern ECS and wet deposition of coarse-mode NO₃^- over the southeastern ECS were also found to be significant processes determining N deposition over the ECS.

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.

Study on the variation of air pollutant concentration and its formation mechanism during the COVID-19 period in Wuhan

To prevent the spread of COVID-19 (2019 novel coronavirus), from January 23 to April 8 in 2020, the highest Class 1 Response was ordered in Wuhan, requiring all residents to stay at home unless absolutely necessary. This action was implemented to cut down all unnecessary human activities, including industry, agriculture and transportation. Reducing these activities to a very low level during these hard times meant that some unprecedented naturally occurring measures of controlling emissions were executed. Ironically, however, after these measures were implemented, ozone levels increased by 43.9%. Also worthy of note, PM_2.5 decreased 31.7%, which was found by comparing the observation data in Wuhan during the epidemic from 8th Feb. to 8th Apr. in 2020 with the same periods in 2019. Utilizing CMAQ (The Community Multiscale Air Quality modeling system), this article investigated the reason for these phenomena based on four sets of numerical simulations with different schemes of emission reduction. Comparing the four sets of simulations with observation, it was deduced that the emissions should decrease to approximately 20% from the typical industrial output, and 10% from agriculture and transportation sources, attributed to the COVID-19 lockdown in Wuhan. More importantly, through the CMAQ process analysis, this study quantitatively analyzed differences of the physical and chemical processes that were affected by the COVID-19 lockdown. It then examined the differences of the COVID-19 lockdown impact and determined the physical and chemical processes between when the pollution increased and decreased, determining the most affected period of the day. As a result, this paper found that (1) PM_2.5 decreased mainly due to the reduction of emission and the contrary contribution of aerosol processes. The North-East wind was also in favor of the decreasing of PM_2.5. (2) O₃ increased mainly due to the slowing down of chemical consumption processes, which made the concentration change of O₃ pollution higher at about 4 p.m.-7 p.m. of the day, while increasing the concentration of O₃ at night during the COVID-19 lockdown in Wuhan. The higher O₃ concentration in the North-East of the main urban area also contributed to the increasing of O₃ with unfavorable wind direction.

Study on the variation of air pollutant concentration and its formation mechanism during the COVID-19 period in Wuhan

To prevent the spread of COVID-19 (2019 novel coronavirus), from January 23 to April 8 in 2020, the highest Class 1 Response was ordered in Wuhan, requiring all residents to stay at home unless absolutely necessary. This action was implemented to cut down all unnecessary human activities, including industry, agriculture and transportation. Reducing these activities to a very low level during these hard times meant that some unprecedented naturally occurring measures of controlling emissions were executed. Ironically, however, after these measures were implemented, ozone levels increased by 43.9%. Also worthy of note, PM_2.5 decreased 31.7%, which was found by comparing the observation data in Wuhan during the epidemic from 8th Feb. to 8th Apr. in 2020 with the same periods in 2019. Utilizing CMAQ (The Community Multiscale Air Quality modeling system), this article investigated the reason for these phenomena based on four sets of numerical simulations with different schemes of emission reduction. Comparing the four sets of simulations with observation, it was deduced that the emissions should decrease to approximately 20% from the typical industrial output, and 10% from agriculture and transportation sources, attributed to the COVID-19 lockdown in Wuhan. More importantly, through the CMAQ process analysis, this study quantitatively analyzed differences of the physical and chemical processes that were affected by the COVID-19 lockdown. It then examined the differences of the COVID-19 lockdown impact and determined the physical and chemical processes between when the pollution increased and decreased, determining the most affected period of the day. As a result, this paper found that (1) PM_2.5 decreased mainly due to the reduction of emission and the contrary contribution of aerosol processes. The North-East wind was also in favor of the decreasing of PM_2.5. (2) O₃ increased mainly due to the slowing down of chemical consumption processes, which made the concentration change of O₃ pollution higher at about 4 p.m.-7 p.m. of the day, while increasing the concentration of O₃ at night during the COVID-19 lockdown in Wuhan. The higher O₃ concentration in the North-East of the main urban area also contributed to the increasing of O₃ with unfavorable wind direction.

Role of sea fog over the Yellow Sea on air quality with the direct effect of aerosols

In this study, we investigate the impact of sea fog over the Yellow Sea on air quality with the direct effect of aerosols for the entire year of 2016. Using the WRF-CMAQ two-way coupled model, we perform four model simulations with the up-to-date emission inventory over East Asia and dynamic chemical boundary conditions provided by hemispheric model simulations. During the spring of 2016, prevailing westerly winds and anticyclones caused the formation of a temperature inversion over the Yellow Sea, providing favorable conditions for the formation of fog. The inclusion of the direct effect of aerosols enhanced its strength. On foggy days, we find dominant changes of aerosols at an altitude of 150-200 m over the Yellow Sea resulted by the production through aqueous chemistry (~12.36% and ~3.08% increases in sulfate and ammonium) and loss via the wet deposition process (~-2.94% decrease in nitrate); we also find stronger wet deposition of all species occurring in PBL. Stagnant conditions associated with reduced air temperature caused by the direct effect of aerosols enhanced aerosol chemistry, especially in coastal regions, and it exceeded the loss of nitrate. The transport of air pollutants affected by sea fog extended to a much broader region. Our findings show that the Yellow Sea acts as not only a path of long-range transport but also as a sink and source of air pollutants. Further study should investigate changes in the impact of sea fog on air quality in conjunction with changes in the concentrations of aerosols and the climate.

Unexpected air quality impacts from implementation of green infrastructure in urban environments: A Kansas City case study

Green infrastructure (GI) implementation can benefit an urban environment by reducing the impacts of urban stormwater on aquatic ecosystems and human health. However, few studies have systematically analyzed the biophysical effects on regional meteorology and air quality that are triggered by changes in the urban vegetative coverage. In this study we use a state-of-the-art high-resolution air quality model to simulate the effects of a hypothetically feasible vegetation-focused GI implementation scenario in Kansas City, MO/KS on regional meteorology and air quality. Full year simulations are conducted for both the base case and GI land use scenarios using two different land surface models (LSMs) schemes inside the meteorological model. While the magnitudes of the changes in air quality due to the GI implementation differ using the two LSMs, the model outputs consistently showed increases in summertime PM_2.5 (1.1 μg m^-3, approximately 10% increase using NOAH LSM), which occurred mostly during the night and arose from the primary components, due to the cooler surface temperatures and the decreased planetary boundary layer height (PBLH). Both the maximum daily 8-hour average ozone and 1 h daily maximum O₃ during summertime, decreased over the downtown areas (maximum decreases of 0.9 and 1.4 ppbv respectively). The largest ozone decreases were simulated to happen during the night, mainly caused by the titration effect of increased NO_x concentration from the lower PBLH. These results highlight the region-specific non-linear process feedback from GI on regional air quality, and further demonstrate the need for comprehensive coupled meteorological-air quality modeling systems and necessity of accurate land surface model for studying these impacts.

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.

A multiphase CMAQ version 5.0 adjoint

We present the development of a multiphase adjoint for the Community Multiscale Air Quality (CMAQ) model, a widely used chemical transport model. The adjoint model provides location- and time-specific gradients that can be used in various applications such as backward sensitivity analysis, source attribution, optimal pollution control, data assimilation, and inverse modeling. The science processes of the CMAQ model include gas-phase chemistry, aerosol dynamics and thermodynamics, cloud chemistry and dynamics, diffusion, and advection. Discrete adjoints are implemented for all the science processes, with an additional continuous adjoint for advection. The development of discrete adjoints is assisted with algorithmic differentiation (AD) tools. Particularly, the Kinetic PreProcessor (KPP) is implemented for gas-phase and aqueous chemistry, and two different automatic differentiation tools are used for other processes such as clouds, aerosols, diffusion, and advection. The continuous adjoint of advection is developed manually. For adjoint validation, the brute-force or finite-difference method (FDM) is implemented process by process with box- or column-model simulations. Due to the inherent limitations of the FDM caused by numerical round-off errors, the complex variable method (CVM) is adopted where necessary. The adjoint model often shows better agreement with the CVM than with the FDM. The adjoints of all science processes compare favorably with the FDM and CVM. In an example application of the full multiphase adjoint model, we provide the first estimates of how emissions of particulate matter (PM2.5) affect public health across the US.

The Acidity of Atmospheric Particles and Clouds

Acidity, defined as pH, is a central component of aqueous chemistry. In the atmosphere, the acidity of condensed phases (aerosol particles, cloud water, and fog droplets) governs the phase partitioning of semi-volatile gases such as HNO3, NH3, HCl, and organic acids and bases as well as chemical reaction rates. It has implications for the atmospheric lifetime of pollutants, deposition, and human health. Despite its fundamental role in atmospheric processes, only recently has this field seen a growth in the number of studies on particle acidity. Even with this growth, many fine particle pH estimates must be based on thermodynamic model calculations since no operational techniques exist for direct measurements. Current information indicates acidic fine particles are ubiquitous, but observationally-constrained pH estimates are limited in spatial and temporal coverage. Clouds and fogs are also generally acidic, but to a lesser degree than particles, and have a range of pH that is quite sensitive to anthropogenic emissions of sulfur and nitrogen oxides, as well as ambient ammonia. Historical measurements indicate that cloud and fog droplet pH has changed in recent decades in response to controls on anthropogenic emissions, while the limited trend data for aerosol particles indicates acidity may be relatively constant due to the semi-volatile nature of the key acids and bases and buffering in particles. This paper reviews and synthesizes the current state of knowledge on the acidity of atmospheric condensed phases, specifically particles and cloud droplets. It includes recommendations for estimating acidity and pH, standard nomenclature, a synthesis of current pH estimates based on observations, and new model calculations on the local and global scale.

Evaluation of Regional Air Quality Models over Sydney, Australia: Part 2, Comparison of PM2.5 and Ozone

Accurate air quality modelling is an essential tool, both for strategic assessment (regulation development for emission controls) and for short-term forecasting (enabling warnings to be issued to protect vulnerable members of society when the pollution levels are predicted to be high). Model intercomparison studies are a valuable support to this work, being useful for identifying any issues with air quality models, and benchmarking their performance against international standards, thereby increasing confidence in their predictions. This paper presents the results of a comparison study of six chemical transport models which have been used to simulate short-term hourly to 24 hourly concentrations of fine particulate matter less than and equal to 2.5 µm in diameter (PM2.5) and ozone (O3) for Sydney, Australia. Model performance was evaluated by comparison to air quality measurements made at 16 locations for O3 and 5 locations for PM2.5, during three time periods that coincided with major atmospheric composition measurement campaigns in the region. These major campaigns included daytime measurements of PM2.5 composition, and so model performance for particulate sulfate (SO42−), nitrate (NO3−), ammonium (NH4+) and elemental carbon (EC) was evaluated at one site per modelling period. Domain-wide performance of the models for hourly O3 was good, with models meeting benchmark criteria and reproducing the observed O3 production regime (based on the O3/NOx indicator) at 80% or more of the sites. Nevertheless, model performance was worse at high (and low) O3 percentiles. Domain-wide model performance for 24 h average PM2.5 was more variable, with a general tendency for the models to under-predict PM2.5 concentrations during the summer and over-predict PM2.5 concentrations in the autumn. The modelling intercomparison exercise has led to improvements in the implementation of these models for Sydney and has increased confidence in their skill at reproducing observed atmospheric composition.

Differences in Model Performance and Source Sensitivities for Sulfate Aerosol Resulting from Updates of the Aqueous- and Gas-Phase Oxidation Pathways for a Winter Pollution Episode in Tokyo, Japan

During the Japanese intercomparison study, Japan’s Study for Reference Air Quality Modeling (J-STREAM), it was found that wintertime SO42– concentrations were underestimated over Japan with the Community Multiscale Air Quality (CMAQ) modeling system. Previously, following two development phases, model performance was improved by refining the Fe- and Mn-catalyzed oxidation pathways and by including an additional aqueous-phase pathway via NO2 oxidation. In a third phase, we examined a winter haze period in December 2016, involving a gas-phase oxidation pathway whereby three stabilized Criegee intermediates (SCI) were incorporated into the model. We also included options for a kinetic mass transfer aqueous-phase calculation. According to statistical analysis, simulations compared well with hourly SO42– observations in Tokyo. Source sensitivities for four domestic emission sources (transportation, stationary combustion, fugitive VOC, and agricultural NH3) were investigated. During the haze period, contributions from other sources (overseas and volcanic emissions) dominated, while domestic sources, including transportation and fuel combustion, played a role in enhancing SO42– concentrations around Tokyo Bay. Updating the aqueous phase metal catalyzed and NO2 oxidation pathways lead to increase contribution from other sources, and the additional gas phase SCI chemistry provided a link between fugitive VOC emission and SO42– concentration via changes in O3 concentration.

Vapor-pressure pathways initiate but hydrolysis products dominate the aerosol estimated from organic nitrates

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

A Feasible Methodological Framework for Uncertainty Analysis and Diagnosis of Atmospheric Chemical Transport Models

The current state of quantifying uncertainty in chemical transport models (CTM) is often limited and insufficient due to numerous uncertainty sources and inefficient or inaccurate uncertainty propagation methods. In this study, we proposed a feasible methodological framework for CTM uncertainty analysis, featuring sensitivity analysis to filter for important model inputs and a new reduced-form model (RFM) that couples the high-order decoupled direct method (HDDM) and the stochastic response surface model (SRSM) to boost uncertainty propagation. Compared with the SRSM, the new RFM approach is 64% more computationally efficient while maintaining high accuracy. The framework was applied to PM_2.5 simulations in the Pearl River Delta (PRD) region and found five precursor emissions, two  pollutants in lateral boundary conditions (LBCs), and three meteorological inputs out of 203 model inputs to be important model inputs based on sensitivity analysis. Among these selected inputs, primary PM_2.5 emissions, PM_2.5 concentrations of LBCs, and wind speed were identified as key uncertainty sources, which collectively contributed 81.4% to the total uncertainty in PM_2.5 simulations. Also, when evaluated against observations, we found that there were systematic underestimates in PM_2.5 simulations, which can be attributed to the two-product method that describes the formation of secondary organic aerosol.

Influence of bromine and iodine chemistry on annual, seasonal, diurnal, and background ozone: CMAQ simulations over the Northern Hemisphere

Bromine and iodine chemistry has been updated in the Community Multiscale Air Quality (CMAQ) model to better capture the influence of natural emissions from the oceans on ozone concentrations. Annual simulations were performed using the hemispheric CMAQ model without and with bromine and iodine chemistry. Model results over the Northern Hemisphere show that including bromine and iodine chemistry in CMAQ not only reduces ozone concentrations within the marine boundary layer but also aloft and inland. Bromine and iodine chemistry reduces annual mean surface ozone over seawater by 25%, with lesser ozone reductions over land. The bromine and iodine chemistry decreases ozone concentration without changing the diurnal profile and is active throughout the year. However, it does not have a strong seasonal influence on ozone over the Northern Hemisphere. Model performance of CMAQ is improved by the bromine and iodine chemistry when compared to observations, especially at coastal sites and over seawater. Relative to bromine, iodine chemistry is approximately four times more effective in reducing ozone over seawater over the Northern Hemisphere (on an annual basis). Model results suggest that the chemistry modulates intercontinental transport and lowers the background ozone imported to the United States.

Fire behavior and smoke modeling: Model improvement and measurement needs for next-generation smoke research and forecasting systems

There is an urgent need for next-generation smoke research and forecasting (SRF) systems to meet the challenges of the growing air quality, health, and safety concerns associated with wildland fire emissions. This review paper presents simulations and experiments of hypothetical prescribed burns with a suite of selected fire behavior and smoke models and identifies major issues for model improvement and the most critical observational needs. The results are used to understand the new and improved capability required for the next-generation SRF systems and to support the design of the Fire and Smoke Model Evaluation Experiment (FASMEE) and other field campaigns. The next-generation SRF systems should have more coupling of fire, smoke, and atmospheric processes to better simulate and forecast vertical smoke distributions and multiple sub-plumes, dynamical and high-resolution fire processes, and local and regional smoke chemistry during day and night. The development of the coupling capability requires comprehensive and spatially and temporally integrated measurements across the various disciplines to characterize flame and energy structure (e.g., individual cells, vertical heat profile and the height of well mixing flaming gases), smoke structure (vertical distributions and multiple sub-plumes), ambient air processes (smoke eddy, entrainment and radiative effects of smoke aerosols), fire emissions (for different fuel types and combustion conditions from flaming to residual smoldering), as well as night-time processes (smoke drainage and super-fog formation).

Model Performance Differences in Sulfate Aerosol in Winter over Japan Based on Regional Chemical Transport Models of CMAQ and CAMx

Sulfate aerosol (SO42−) is a major component of particulate matter in Japan. The Japanese model intercomparison study, J-STREAM, found that although SO42− is well captured by models, it is underestimated during winter. In the first phase of J-STREAM, we refined the Fe- and Mn-catalyzed oxidation and partly improved the underestimation. The winter haze in December 2016 was a target period in the second phase. The results from the Community Multiscale Air Quality (CMAQ) and Comprehensive Air quality Model with eXtentions (CAMx) regional chemical transport models were compared with observations from the network over Japan and intensive observations at Nagoya and Tokyo. Statistical analysis showed both models satisfied the suggested model performance criteria. CMAQ sensitivity simulations explained the improvements in model performance. CMAQ modeled lower SO42− concentrations than CAMx, despite increased aqueous oxidation via the metal catalysis pathway and NO2 reaction in CMAQ. Deposition explained this difference. A scatter plot demonstrated that the lower SO42− concentration in CMAQ than in CAMx arose from the lower SO2 concentration and higher SO42− wet deposition in CMAQ. The dry deposition velocity caused the difference in SO2 concentration. These results suggest the importance of deposition in improving our understanding of ambient concentration behavior.

Spatial variation of modelled total, dry and wet nitrogen deposition to forests at global scale

Forests are an important biome that covers about one third of the global land surface and provides important ecosystem services. Since atmospheric deposition of nitrogen (N) can have both beneficial and deleterious effects, it is important to quantify the amount of N deposition to forest ecosystems. Measurements of N deposition to the numerous forest biomes across the globe are scarce, so chemical transport models are often used to provide estimates of atmospheric N inputs to these ecosystems. We provide an overview of approaches used to calculate N deposition in commonly used chemical transport models. The Task Force on Hemispheric Transport of Air Pollution (HTAP2) study intercompared N deposition values from a number of global chemical transport models. Using a multi-model mean calculated from the HTAP2 deposition values, we map N deposition to global forests to examine spatial variations in total, dry and wet deposition. Highest total N deposition occurs in eastern and southern China, Japan, Eastern U.S. and Europe while the highest dry deposition occurs in tropical forests. The European Monitoring and Evaluation Program (EMEP) model predicts grid-average deposition, but also produces deposition by land use type allowing us to compare deposition specifically to forests with the grid-average value. We found that, for this study, differences between the grid-average and forest specific could be as much as a factor of two and up to more than a factor of five in extreme cases. This suggests that consideration should be given to using forest-specific deposition for input to ecosystem assessments such as critical loads determinations.

Photochemical model evaluation of 2013 California wild fire air quality impacts using surface, aircraft, and satellite data

The Rim Fire was one of the largest wildfires in California history, burning over 250,000 acres during August and September 2013 affecting air quality locally and regionally in the western U.S. Routine surface monitors, remotely sensed data, and aircraft based measurements were used to assess how well the Community Multiscale Air Quality (CMAQ) photochemical grid model applied at 4 and 12 km resolution represented regional plume transport and chemical evolution during this extreme wildland fire episode. Impacts were generally similar at both grid resolutions although notable differences were seen in some secondary pollutants (e.g., formaldehyde and peroxyacyl nitrate) near the Rim fire. The modeling system does well at capturing near-fire to regional scale smoke plume transport compared to remotely sensed aerosol optical depth (AOD) and aircraft transect measurements. Plume rise for the Rim fire was well characterized as the modeled plume top was consistent with remotely sensed data and the altitude of aircraft measurements, which were typically made at the top edge of the plume. Aircraft-based lidar suggests O₃ downwind in the Rim fire plume was vertically stratified and tended to be higher at the plume top, while CMAQ estimated a more uniformly mixed column of O₃. Predicted wildfire ozone (O₃) was overestimated both at the plume top and at nearby rural and urban surface monitors. Photolysis rates were well characterized by the model compared with aircraft measurements meaning aerosol attenuation was reasonably estimated and unlikely contributing to O₃ overestimates at the top of the plume. Organic carbon was underestimated close to the Rim fire compared to aircraft data, but was consistent with nearby surface measurements. Periods of elevated surface PM_2.5 at rural monitors near the Rim fire were not usually coincident with elevated O₃.

Refinement of Modeled Aqueous-Phase Sulfate Production via the Fe- and Mn-Catalyzed Oxidation Pathway

We refined the aqueous-phase sulfate (SO42−) production in the state-of-the-art Community Multiscale Air Quality (CMAQ) model during the Japanese model inter-comparison project, known as Japan’s Study for Reference Air Quality Modeling (J-STREAM). In Japan, SO42− is the major component of PM2.5, and CMAQ reproduces the observed seasonal variation of SO42− with the summer maxima and winter minima. However, CMAQ underestimates the concentration during winter over Japan. Based on a review of the current modeling system, we identified a possible reason as being the inadequate aqueous-phase SO42− production by Fe- and Mn-catalyzed O2 oxidation. This is because these trace metals are not properly included in the Asian emission inventories. Fe and Mn observations over Japan showed that the model concentrations based on the latest Japanese emission inventory were substantially underestimated. Thus, we conducted sensitivity simulations where the modeled Fe and Mn concentrations were adjusted to the observed levels, the Fe and Mn solubilities were increased, and the oxidation rate constant was revised. Adjusting the concentration increased the SO42− concentration during winter, as did increasing the solubilities and revising the rate constant to consider pH dependencies. Statistical analysis showed that these sensitivity simulations improved model performance. The approach adopted in this study can partly improve model performance in terms of the underestimation of SO42− concentration during winter. From our findings, we demonstrated the importance of developing and evaluating trace metal emission inventories in Asia.

Southeast Atmosphere Studies: learning from model-observation syntheses

Concentrations of atmospheric trace species in the United States have changed dramatically over the past several decades in response to pollution control strategies, shifts in domestic energy policy and economics, and economic development (and resulting emission changes) elsewhere in the world. Reliable projections of the future atmosphere require models to not only accurately describe current atmospheric concentrations, but to do so by representing chemical, physical and biological processes with conceptual and quantitative fidelity. Only through incorporation of the processes controlling emissions and chemical mechanisms that represent the key transformations among reactive molecules can models reliably project the impacts of future policy, energy and climate scenarios. Efforts to properly identify and implement the fundamental and controlling mechanisms in atmospheric models benefit from intensive observation periods, during which collocated measurements of diverse, speciated chemicals in both the gas and condensed phases are obtained. The Southeast Atmosphere Studies (SAS, including SENEX, SOAS, NOMADSS and SEAC4RS) conducted during the summer of 2013 provided an unprecedented opportunity for the atmospheric modeling community to come together to evaluate, diagnose and improve the representation of fundamental climate and air quality processes in models of varying temporal and spatial scales. This paper is aimed at discussing progress in evaluating, diagnosing and improving air quality and climate modeling using comparisons to SAS observations as a guide to thinking about improvements to mechanisms and parameterizations in models. The effort focused primarily on model representation of fundamental atmospheric processes that are essential to the formation of ozone, secondary organic aerosol (SOA) and other trace species in the troposphere, with the ultimate goal of understanding the radiative impacts of these species in the southeast and elsewhere. Here we address questions surrounding four key themes: gas-phase chemistry, aerosol chemistry, regional climate and chemistry interactions, and natural and anthropogenic emissions. We expect this review to serve as a guidance for future modeling efforts.