Acidity and the multiphase chemistry of atmospheric aqueous particles and clouds

The acidity of aqueous atmospheric solutions is a key parameter driving both the partitioning of semi-volatile acidic and basic trace gases and their aqueous-phase chemistry. In addition, the acidity of atmospheric aqueous phases, e.g., deliquesced aerosol particles, cloud, and fog droplets, is also dictated by aqueous-phase chemistry. These feedbacks between acidity and chemistry have crucial implications for the tropospheric lifetime of air pollutants, atmospheric composition, deposition to terrestrial and oceanic ecosystems, visibility, climate, and human health. Atmospheric research has made substantial progress in understanding feedbacks between acidity and multiphase chemistry during recent decades. This paper reviews the current state of knowledge on these feedbacks with a focus on aerosol and cloud systems, which involve both inorganic and organic aqueous-phase chemistry. Here, we describe the impacts of acidity on the phase partitioning of acidic and basic gases and buffering phenomena. Next, we review feedbacks of different acidity regimes on key chemical reaction mechanisms and kinetics, as well as uncertainties and chemical subsystems with incomplete information. Finally, we discuss atmospheric implications and highlight the need for future investigations, particularly with respect to reducing emissions of key acid precursors in a changing world, and the need for advancements in field and laboratory measurements and model tools.

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.

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.

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.

Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1

The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency's (EPA) Office of Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O₃) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM_2.5) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O₃ mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM_2.5 bias (PM_2.5 is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O₃ and PM_2.5 on average in January and July. Overall, the seasonal variation in simulated PM_2.5 generally improves in CMAQv5.1 (when considering all model updates), as simulated PM_2.5 concentrations decrease in the winter (when PM_2.5 is generally overestimated by CMAQ) and increase in the summer (when PM_2.5 is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O₃ mean bias, as O₃ tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O₃ is low); however, O₃ correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NO _x (NO + NO₂), VOC and SO _x (SO₂ + SO₄) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O₃ due to large, widespread reductions in observed emissions.

Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1

The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency's (EPA) Office of Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O₃) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM_2.5) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O₃ mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM_2.5 bias (PM_2.5 is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O₃ and PM_2.5 on average in January and July. Overall, the seasonal variation in simulated PM_2.5 generally improves in CMAQv5.1 (when considering all model updates), as simulated PM_2.5 concentrations decrease in the winter (when PM_2.5 is generally overestimated by CMAQ) and increase in the summer (when PM_2.5 is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O₃ mean bias, as O₃ tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O₃ is low); however, O₃ correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NO _x (NO + NO₂), VOC and SO _x (SO₂ + SO₄) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O₃ due to large, widespread reductions in observed emissions.

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

This paper describes the development and implementation of an extendable aqueous-phase chemistry option (AQCHEM -KMT(I)) for the Community Multiscale Air Quality (CMAQ) modeling system, version 5.1. Here, the Kinetic PreProcessor (KPP), version 2.2.3, is used to generate a Rosenbrock solver (Rodas3) to integrate the stiff system of ordinary differential equations (ODEs) that describe the mass transfer, chemical kinetics, and scavenging processes of CMAQ clouds. CMAQ's standard cloud chemistry module (AQCHEM) is structurally limited to the treatment of a simple chemical mechanism. This work advances our ability to test and implement more sophisticated aqueous chemical mechanisms in CMAQ and further investigate the impacts of microphysical parameters on cloud chemistry. Box model cloud chemistry simulations were performed to choose efficient solver and tolerance settings, evaluate the implementation of the KPP solver, and assess the direct impacts of alternative solver and kinetic mass transfer on predicted concentrations for a range of scenarios. Month-long CMAQ simulations for winter and summer periods over the US reveal the changes in model predictions due to these cloud module updates within the full chemical transport model. While monthly average CMAQ predictions are not drastically altered between AQCHEM and AQCHEM-KMT, hourly concentration differences can be significant. With added in-cloud secondary organic aerosol (SOA) formation from biogenic epoxides (AQCHEM-KMTI), normalized mean error and bias statistics are slightly improved for 2-methyltetrols and 2-methylglyceric acid at the Research Triangle Park measurement site in North Carolina during the Southern Oxidant and Aerosol Study (SOAS) period. The added in-cloud chemistry leads to a monthly average increase of 11-18 % in "cloud" SOA at the surface in the eastern United States for June 2013.

Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1

The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency's (EPA) Office of Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O₃) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM_2.5) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O₃ mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM_2.5 bias (PM_2.5 is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O₃ and PM_2.5 on average in January and July. Overall, the seasonal variation in simulated PM_2.5 generally improves in CMAQv5.1 (when considering all model updates), as simulated PM_2.5 concentrations decrease in the winter (when PM_2.5 is generally overestimated by CMAQ) and increase in the summer (when PM_2.5 is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O₃ mean bias, as O₃ tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O₃ is low); however, O₃ correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NO _x (NO + NO₂), VOC and SO _x (SO₂ + SO₄) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O₃ due to large, widespread reductions in observed emissions.

Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1

The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency's (EPA) Office of Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O3) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM2.5) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O3 mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM2.5 bias (PM2.5 is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O3 and PM2.5 on average in January and July. Overall, the seasonal variation in simulated PM2.5 generally improves in CMAQv5.1 (when considering all model updates), as simulated PM2.5 concentrations decrease in the winter (when PM2.5 is generally overestimated by CMAQ) and increase in the summer (when PM2.5 is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O3 mean bias, as O3 tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O3 is low); however, O3 correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NO x (NO + NO2), VOC and SO x (SO2 + SO4) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O3 due to large, widespread reductions in observed emissions.

Near-road modeling and measurement of cerium-containing particles generated by nanoparticle diesel fuel additive use

Cerium oxide nanoparticles (nCe) are used as a fuel-borne catalyst in diesel engines to reduce particulate emissions, yet the environmental and human health impacts of the exhaust particles are not well understood. To bridge the gap between emission measurements and ambient impacts, size-resolved measurements of particle composition and mass concentration have been performed in Newcastle-upon-Tyne, United Kingdom, where buses have used an nCe additive since 2005. These observations show that the noncrustal cerium fraction thought to be associated with the use of nCe has a mass concentration ∼ 0.3 ng m(-3) with a size distribution peaking at 100-320 nm in aerodynamic diameter. Simulations with a near-roadway multicomponent sectional aerosol dynamic model predict that the use of nCe additives increases the number concentration of nuclei mode particles (<50 nm in diameter) while decreasing the total mass concentration. The near-road model predicts a downwind mass size distribution of cerium-containing particles peaking at 150 nm in aerodynamic diameter, a value similar to that measured for noncrustal cerium in Newcastle. This work shows that both the emission and atmospheric transformation of cerium-containing particles needs to be taken into account by regional modelers, exposure scientists, and policymakers when determining potential environmental and human health impacts.