A simple lightning assimilation technique for improving retrospective WRF simulations

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.

Establishing the Suitability of the Model for Prediction Across Scales for Global Retrospective Air Quality Modeling

The U.S. EPA is leveraging recent advances in meteorological modeling to construct an air quality modeling system to allow consistency from global to local scales. The Model for Prediction Across Scales-Atmosphere (MPAS-A or MPAS) has been developed by the National Center for Atmospheric Research (NCAR) as a global complement to the Weather Research and Forecasting model (WRF). Patterned after a regional coupled system with WRF, the Community Multiscale Air Quality (CMAQ) modeling system has been coupled within MPAS to explore global-to-local chemical transport modeling. Several options were implemented into MPAS for retrospective applications. Nudging-based data assimilation was added to support continuous simulations of past weather to minimize error growth that exists with a weather forecast configuration. The Pleim-Xiu land-surface model, the Asymmetric Convective Model 2 boundary layer scheme, and the Pleim surface layer scheme were added as the preferred options for retrospective air quality applications with WRF. Annual simulations were conducted using this EPA-enhanced MPAS configuration on two different mesh structures and compared against WRF. MPAS generally compares well with WRF over the conterminous United States. Errors in MPAS surface meteorology are comparable to WRF throughout the year. Precipitation statistics indicate MPAS performs slightly better than WRF. Solar radiation in MPAS is higher than WRF and measurements, suggesting fewer clouds in MPAS than WRF. Upper-air meteorology is well-simulated by MPAS, but errors are slightly higher than WRF. These comparisons lend confidence to use MPAS for retrospective air quality modeling and suggest ways it can be further improved in the future.

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.

Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model

Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T _g), oxygen to carbon (O : C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core-shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~ 70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10²-10¹² Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity < 10² Pa s. Phase separation while in a liquid phase state, i.e., liquid-liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM_2.5) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.

Significant ground-level ozone attributed to lightning-induced nitrogen oxides during summertime over the Mountain West States

Using lightning flash data from the National Lightning Detection Network with an updated lightning nitrogen oxides (NOx) emission estimation algorithm in the Community Multiscale Air Quality (CMAQ) model, we estimate the hourly variations in lightning NOx emissions for the summer of 2011 and simulate its impact on distributions of tropospheric ozone (O3) across the continental United States. We find that typical summer-time lightning activity across the U.S. Mountain West States (MWS) injects NOx emissions comparable to those from anthropogenic sources into the troposphere over the region. Comparison of two model simulation cases with and without lightning NOx emissions show that significant amount of ground-level O3 in the MWS during the summer can be attributed to the lightning NOX emissions. The simulated surface-level O3 from a model configuration incorporating lightning NOx emissions showed better agreement with the observed values than the model configuration without lightning NOx emissions. The time periods of significant reduction in bias in simulated O3 between these two cases strongly correlate with the time periods when lightning activity occurred in the region. The inclusion of lightning NOx increased daily maximum 8 h O3 by up to 17 ppb and improved model performance relative to measured surface O3 mixing ratios in the MWS region. Analysis of model results in conjunction with lidar measurements at Boulder, Colorado during July 2014 corroborated similar impacts of lightning NOx emissions on O3 emissions estimated for other summers is comparable to the 2011 air quality. The magnitude of lightning NOx estimates suggesting that summertime surface-level O3 levels in the MWS region could be significantly influenced by lightning NOx.

An Integrated Agriculture, Atmosphere, and Hydrology Modeling System for Ecosystem Assessments

We present a regional-scale integrated modeling system (IMS) that includes Environmental Policy Integrated Climate (EPIC), Weather Research and Forecast (WRF), Community Multiscale Air Quality (CMAQ), and Soil and Water Assessment Tool (SWAT) models. The centerpiece of the IMS is the Fertilizer Emission Scenario Tool for CMAQ (FEST-C), which includes a Java-based interface and EPIC adapted to regional applications along with built-in database and tools. The SWAT integration capability is a key enhanced feature in the current release of FEST-C v1.4. For integrated modeling demonstration and evaluation, FEST-C EPIC is simulated over three individual years with WRF/CMAQ weather and N deposition. Simulated yearly changes in water and N budgets along with yields for two major crops (corn grain and soybean) match those inferred from intuitive physical reasoning and survey data given different-year weather conditions. Yearlong air quality simulations with an improved bidirectional ammonia flux modeling approach directly using EPIC-simulated soil properties including NH₃ content helps reduce biases of simulated gas-phase NH₃ and NH₄ ⁺ wet deposition over the growing season. Integrated hydrology and water quality simulations applied to the Mississippi River Basin show that estimated monthly streamflow and dissolved N near the outlet to the Gulf of Mexico display similar seasonal patterns as observed. Limitations and issues in different parts of the integrated multimedia simulations are identified and discussed to target areas for future improvements.

PLAIN LANGUAGE SUMMARY

Computer modeling tools with land-water-air processes are important for understanding nutrient cycling and its negative impacts on air and water quality. We have developed an integrated modeling system that includes agriculture, atmosphere, and hydrology components. The centerpiece of the system is a computer system that includes an agricultural ecosystem model and tools used to connect different modeling components. The agricultural system can conduct simulations for 42 types of grassland and cropland with the influence of site, soil, and management information along with weather and nitrogen deposition from the atmosphere component. An air quality computer model then uses information from the agricultural model, such as how much ammonia is in the soil, to predict how much ammonia gets in the air. Then, the watershed hydrology and water quality model uses the information from the agricultural and atmospheric models to understand the influence of agriculture and atmosphere on water quality. The paper demonstrates and evaluates the integrated modeling system on issues mainly related to N cycling. The system performs reasonably well in comparison with survey and observation data given the configured modeling constraints. The paper also identifies and discusses the advantages and limitations in each part of the system for future applications and improvements.

New Bidirectional Ammonia Flux Model in an Air Quality Model Coupled With an Agricultural Model

Ammonia surface flux is bidirectional; that is, net flux can be either upward or downward. In fertilized agricultural croplands and grasslands there is usually more emission than deposition especially in midday during warmer seasons. In North America, most of the ammonia emissions are from agriculture with a significant fraction of that coming from fertilizer. A new bidirectional ammonia flux modeling system has been developed in the Community Multiscale Air Quality (CMAQ) model, which has close linkages with the Environmental Policy Integrated Climate (EPIC) agricultural ecosystem model. Daily inputs from EPIC are used to calculate soil ammonia concentrations that are combined with air concentrations in CMAQ to calculate bidirectional surface flux. The model is evaluated against surface measurements of NH₃ concentrations, NH₄ ⁺ and SO₄ ^2- aerosol concentrations, NH₄ ⁺ wet deposition measurements, and satellite retrievals of NH₃ concentrations. The evaluation shows significant improvement over the base model without bidirectional ammonia flux. Comparisons to monthly average satellite retrievals show similar spatial distribution with the highest ammonia concentrations in the Central Valley of California (CA), the Snake River valley in Idaho, and the western High Plains. In most areas the model underestimates, but in the Central Valley of CA, it generally overestimates ammonia concentration. Case study analyses indicate that modeled high fluxes of ammonia in CA are often caused by anomalous high soil ammonia loading from EPIC for particular crop types. While further improvements to parameterizations in EPIC and CMAQ are recommended, this system is a significant advance over previous ammonia bidirectional surface flux models.

Simulating lightning NO production in CMAQv5.2: evolution of scientific updates

This work describes the lightning nitric oxide (LNO) production schemes in the Community Multiscale Air Quality (CMAQ) model. We first document the existing LNO production scheme and vertical distribution algorithm. We then describe updates that were made to the scheme originally based on monthly National Lightning Detection Network (mNLDN) observations. The updated scheme uses hourly NLDN (hNLDN) observations. These NLDN-based schemes are good for retrospective model applications when historical lightning data are available. For applications when observed data are not available (i.e., air quality forecasts and climate studies that assume similar climate conditions), we have developed a scheme that is based on linear and log-linear parameters derived from regression of multiyear historical NLDN (pNLDN) observations and meteorological model simulations. Preliminary assessment for total column LNO production reveals that the mNLDN scheme overestimates LNO by over 40% during summer months compared with the updated hNLDN scheme that reflects the observed lightning activity more faithfully in time and space. The pNLDN performance varies with year, but it generally produced LNO columns that are comparable to hNLDN and mNLDN, and in most cases it outperformed mNLDN. Thus, when no observed lightning data are available, pNLDN can provide reasonable estimates of LNO emissions over time and space for this important natural NO source that influences air quality regulations.

Experimental and model estimates of the contributions from biogenic monoterpenes and sesquiterpenes to secondary organic aerosol in the southeastern United States

Atmospheric organic aerosol (OA) has important impacts on climate and human health but its sources remain poorly understood. Biogenic monoterpenes and sesquiterpenes are important precursors of secondary organic aerosol (SOA), but the amounts and pathways of SOA generation from these precursors are not well constrained by observations. We propose that the less-oxidized oxygenated organic aerosol (LO-OOA) factor resolved from positive matrix factorization (PMF) analysis on aerosol mass spectrometry (AMS) data can be used as a surrogate for fresh SOA from monoterpenes and sesquiterpenes in the southeastern US. This hypothesis is supported by multiple lines of evidence, including lab-in-the-field perturbation experiments, extensive ambient ground-level measurements, and state-of-the-art modeling. We performed lab-in-the-field experiments in which the ambient air is perturbed by the injection of selected monoterpenes and sesquiterpenes, and the subsequent SOA formation is investigated. PMF analysis on the perturbation experiments provides an objective link between LO-OOA and fresh SOA from monoterpenes and sesquiterpenes as well as insights into the sources of other OA factors. Further, we use an upgraded atmospheric model and show that modeled SOA concentrations from monoterpenes and sesquiterpenes could reproduce both the magnitude and diurnal variation of LO-OOA at multiple sites in the southeastern US, building confidence in our hypothesis. We estimate the annual average concentration of SOA from monoterpenes and sesquiterpenes in the southeastern US to be roughly 2 μg m-3.

Coupling of organic and inorganic aerosol systems and the effect on gas-particle partitioning in the southeastern US

Several models were used to describe the partitioning of ammonia, water, and organic compounds between the gas and particle phases for conditions in the southeastern US during summer 2013. Existing equilibrium models and frameworks were found to be sufficient, although additional improvements in terms of estimating pure-species vapor pressures are needed. Thermodynamic model predictions were consistent, to first order, with a molar ratio of ammonium to sulfate of approximately 1.6 to 1.8 (ratio of ammonium to 2× sulfate, RN/2S ≈ 0.8 to 0.9) with approximately 70% of total ammonia and ammonium (NH x ) in the particle. Southeastern Aerosol Research and Characterization Network (SEARCH) gas and aerosol and Southern Oxidant and Aerosol Study (SOAS) Monitor for AeRosols and Gases in Ambient air (MARGA) aerosol measurements were consistent with these conditions. CMAQv5.2 regional chemical transport model predictions did not reflect these conditions due to a factor of 3 overestimate of the nonvolatile cations. In addition, gas-phase ammonia was overestimated in the CMAQ model leading to an even lower fraction of total ammonia in the particle. Chemical Speciation Network (CSN) and aerosol mass spectrometer (AMS) measurements indicated less ammonium per sulfate than SEARCH and MARGA measurements and were inconsistent with thermodynamic model predictions. Organic compounds were predicted to be present to some extent in the same phase as inorganic constituents, modifying their activity and resulting in a decrease in [H+]air (H+ in μgm-3 air), increase in ammonia partitioning to the gas phase, and increase in pH compared to complete organic vs. inorganic liquid-liquid phase separation. In addition, accounting for nonideal mixing modified the pH such that a fully interactive inorganic-organic system had a pH roughly 0.7 units higher than predicted using traditional methods (pH = 1.5 vs. 0.7). Particle-phase interactions of organic and inorganic compounds were found to increase partitioning towards the particle phase (vs. gas phase) for highly oxygenated (O : C≥0.6) compounds including several isoprene-derived tracers as well as levoglu-cosan but decrease particle-phase partitioning for low O: C, monoterpene-derived species.

Improving the simulation of convective dust storms in regional-to-global models

Convective dust storms have significant impacts on atmospheric conditions and air quality and are a major source of dust uplift in summertime. However, regional-to-global models generally do not accurately simulate these storms, a limitation that can be attributed to (1) using a single mean value for wind speed per grid box, i.e., not accounting for subgrid wind variability and (2) using convective parametrizations that poorly simulate cold pool outflows. This study aims to improve the simulation of convective dust storms by tackling these two issues. Specifically, we incorporate a probability distribution function for surface wind in each grid box to account for subgrid wind variability due to dry and moist convection. Furthermore, we use lightning assimilation to increase the accuracy of the convective parameterization and simulated cold pool outflows. This updated model framework is used to simulate a massive convective dust storm that hit Phoenix, AZ, on 6 July 2011. The results show that lightning assimilation provides a more realistic simulation of precipitation features, including timing and location, and the resulting cold pool outflows that generated the dust storm. When those results are combined with a dust model that accounts for subgrid wind variability, the prediction of dust uplift and concentrations are considerably improved compared to the default model results. This modeling framework could potentially improve the simulation of convective dust storms in global models, regional climate simulations, and retrospective air quality studies.