A single-point modeling approach for the intercomparison and evaluation of ozone dry deposition across chemical transport models (Activity 2 of AQMEII4).

A primary sink of air pollutants and their precursors is dry deposition. Dry deposition estimates differ across chemical transport models, yet an understanding of the model spread is incomplete. Here, we introduce Activity 2 of the Air Quality Model Evaluation International Initiative Phase 4 (AQMEII4). We examine 18 dry deposition schemes from regional and global chemical transport models as well as standalone models used for impact assessments or process understanding. We configure the schemes as single-point models at eight Northern Hemisphere locations with observed ozone fluxes. Single-point models are driven by a common set of site-specific meteorological and environmental conditions. Five of eight sites have at least 3 years and up to 12 years of ozone fluxes. The interquartile range across models in multiyear mean ozone deposition velocities ranges from a factor of 1.2 to 1.9 annually across sites and tends to be highest during winter compared with summer. No model is within 50 % of observed multiyear averages across all sites and seasons, but some models perform well for some sites and seasons. For the first time, we demonstrate how contributions from depositional pathways vary across models. Models can disagree with respect to relative contributions from the pathways, even when they predict similar deposition velocities, or agree with respect to the relative contributions but predict different deposition velocities. Both stomatal and nonstomatal uptake contribute to the large model spread across sites. Our findings are the beginning of results from AQMEII4 Activity 2, which brings scientists who model air quality and dry deposition together with scientists who measure ozone fluxes to evaluate and improve dry deposition schemes in the chemical transport models used for research, planning, and regulatory purposes.

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.

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

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

2002-2017 anthropogenic emissions data for air quality modeling over the United States.

The United States Environmental Protection Agency (US EPA) has developed a set of annual North American emissions data for multiple air pollutants across 18 broad source categories for 2002 through 2017. The sixteen new annual emissions inventories were developed using consistent input data and methods across all years. When a consistent method or tool was not available for a source category, emissions were estimated by scaling data from the EPA's 2017 National Emissions Inventory with scaling factors based on activity data and/or emissions control information. The emissions datasets are designed to support regional air quality modeling for a wide variety of human health and ecological applications. The data were developed to support simulations of the EPA's Community Multiscale Air Quality model but can also be used by other regional scale air quality models. The emissions data are one component of EPA's Air Quality Time Series Project which also includes air quality modeling inputs (meteorology, initial conditions, boundary conditions) and outputs (e.g., ozone, PM2.5 and constituent species, wet and dry deposition) for the Conterminous US at a 12 km horizontal grid spacing.

4D-Var Inversion of European NH3 Emissions Using CrIS NH3 Measurements and GEOS-Chem Adjoint With Bi-Directional and Uni-Directional Flux Schemes.

We conduct the first 4D-Var inversion of NH3 accounting for NH3 bi-directional flux, using CrIS satellite NH3 observations over Europe in 2016. We find posterior NH3 emissions peak more in springtime than prior emissions at continental to national scales, and annually they are generally smaller than the prior emissions over central Europe, but larger over most of the rest of Europe. Annual posterior anthropogenic NH3 emissions for 25 European Union members (EU25) are 25% higher than the prior emissions and very close (<2% difference) to other inventories. Our posterior annual anthropogenic emissions for EU25, the UK, the Netherlands, and Switzerland are generally 10%-20% smaller than when treating NH3 fluxes as uni-directional emissions, while the monthly regional difference can be up to 34% (Switzerland in July). Compared to monthly mean in-situ observations, our posterior NH3 emissions from both schemes generally improve the magnitude and seasonality of simulated surface NH3 and bulk NH x wet deposition throughout most of Europe, whereas evaluation against hourly measurements at a background site shows the bi-directional scheme better captures observed diurnal variability of surface NH3. This contrast highlights the need for accurately simulating diurnal variability of NH3 in assimilation of sun-synchronous observations and also the potential value of future geostationary satellite observations. Overall, our top-down ammonia emissions can help to examine the effectiveness of air pollution control policies to facilitate future air pollution management, as well as helping us understand the uncertainty in top-down NH3 emissions estimates associated with treatment of NH3 surface exchange.

Technical note: AQMEII4 Activity 1: evaluation of wet and dry deposition schemes as an integral part of regional-scale air quality models.

We present in this technical note the research protocol for phase 4 of the Air Quality Model Evaluation International Initiative (AQMEII4). This research initiative is divided into two activities, collectively having three goals: (i) to define the current state of the science with respect to representations of wet and especially dry deposition in regional models, (ii) to quantify the extent to which different dry deposition parameterizations influence retrospective air pollutant concentration and flux predictions, and (iii) to identify, through the use of a common set of detailed diagnostics, sensitivity simulations, model evaluation, and reduction of input uncertainty, the specific causes for the current range of these predictions. Activity 1 is dedicated to the diagnostic evaluation of wet and dry deposition processes in regional air quality models (described in this paper), and Activity 2 to the evaluation of dry deposition point models against ozone flux measurements at multiple towers with multiyear observations (to be described in future submissions as part of the special issue on AQMEII4). The scope of this paper is to present the scientific protocols for Activity 1, as well as to summarize the technical information associated with the different dry deposition approaches used by the participating research groups of AQMEII4. In addition to describing all common aspects and data used for this multi-model evaluation activity, most importantly, we present the strategy devised to allow a common process-level comparison of dry deposition obtained from models using sometimes very different dry deposition schemes. The strategy is based on adding detailed diagnostics to the algorithms used in the dry deposition modules of existing regional air quality models, in particular archiving diagnostics specific to land use-land cover (LULC) and creating standardized LULC categories to facilitate cross-comparison of LULC-specific dry deposition parameters and processes, as well as archiving effective conductance and effective flux as means for comparing the relative influence of different pathways towards the net or total dry deposition. This new approach, along with an analysis of precipitation and wet deposition fields, will provide an unprecedented process-oriented comparison of deposition in regional air quality models. Examples of how specific dry deposition schemes used in participating models have been reduced to the common set of comparable diagnostics defined for AQMEII4 are also presented.

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

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

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

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

Atmospheric nitrogen deposition in the Chesapeake Bay watershed: A history of change.

The Chesapeake Bay watershed has been the focus of pioneering studies of the role of atmospheric nitrogen (N) deposition as a nutrient source and driver of estuarine trophic status. Here, we review the history and evolution of scientific investigations of the role of atmospheric N deposition, examine trends from wet and dry deposition networks, and present century-long (1950-2050) atmospheric N deposition estimates. Early investigations demonstrated the importance of atmospheric deposition as an N source to the Bay, providing 25%-40% among all major N sources. These early studies led to the unprecedented inclusion of targeted decreases in atmospheric N deposition as part of the multi-stakeholder effort to reduce N loads to the Bay. Emissions of nitrogen oxides (NOx) and deposition of wet nitrate, oxidized dry N, and dry ammonium ( NH 4 + ) sharply and synchronously declined by 60%-73% during 1995-2019. These decreases largely resulted from implementation of Title IV of the 1990 Clean Air Act Amendments, which began in 1995. Wet NH 4 + deposition shows no significant trend during this period. The century-long atmospheric N deposition estimates indicate an increase in total atmospheric N deposition in the Chesapeake watershed from 1950 to a peak of ~15 kg N/ha/yr in 1979, trailed by a slight decline of <10% through the mid-1990s, and followed by a sharp decline of about 40% thereafter through 2019. An additional 21% decline in atmospheric N deposition is projected from 2015 to 2050. A comparison of the Potomac River and James River watersheds indicates higher atmospheric N deposition in the Potomac, likely resulting from greater emissions from higher proportions of agricultural and urban land in this basin. Atmospheric N deposition rose from 30% among all N sources to the Chesapeake Bay watershed in 1950 to a peak of 40% in 1973, and a decline to 28% by 2015. These data highlight the important role of atmospheric N deposition in the Chesapeake Bay watershed and present a potential opportunity for decreases in deposition to contribute to further reducing N loads and improving the trophic status of tidal waters.

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

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

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

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

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

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

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

Per- and polyfluoroalkyl substances (PFASs) have been released into the environment for decades, yet contributions of air emissions to total human exposure, from inhalation and drinking water contamination via deposition, are poorly constrained. The atmospheric transport and fate of a PFAS mixture from a fluoropolymer manufacturing facility in North Carolina were investigated with the Community Multiscale Air Quality (CMAQ) model applied at high resolution (1 km) and extending ∼150 km from the facility. Twenty-six explicit PFAS compounds, including GenX, were added to CMAQ using current best estimates of air emissions and relevant physicochemical properties. The new model, CMAQ-PFAS, predicts that 5% by mass of total emitted PFAS and 2.5% of total GenX are deposited within ∼150 km of the facility, with the remainder transported out. Modeled air concentrations of total GenX and total PFAS around the facility can reach 24.6 and 8500 ng m-3 but decrease to ∼0.1 and ∼10 ng m-3 at 35 km downwind, respectively. We find that compounds with acid functionality have higher deposition due to enhanced water solubility and pH-driven partitioning to aqueous media. To our knowledge, this is the first modeling study of the fate of a comprehensive, chemically resolved suite of PFAS air emissions from a major manufacturing source.

Linking multi-media modeling with machine learning to assess and predict lake Chlorophyll a concentrations.

Eutrophication and excessive algal growth pose a threat on aquatic organisms and the health of the public, environment, and the economy. Understanding what drives excessive algal growth can inform mitigation measures and aid in advance planning to minimize impacts. We demonstrate how simulated data from weather, hydrological, and agroecosystem numerical prediction models can be combined with machine learning (ML) to assess and predict Chlorophyll a (Chl a) concentrations, a proxy for lake eutrophication and algal biomass. The study area is Lake Erie for a 16-year period, 2002-2017. A total of 20 environmental variables from linked and coupled physical models are used as input features to train the ML model with Chl a observations from 16 measuring stations. Included are meteorological variables from the Weather Research and Forecasting (WRF) model, hydrological variables from the Variable Infiltration Capacity (VIC) model, and agricultural management practice variables from the Environmental Policy Integrated Climate (EPIC) agroecosystem model. The consolidation of these variables is conducive to a successful prediction of Chl a. Aside from the synergistic effects that weather, hydrology, and fertilizers have on eutrophication and excessive algal growth, we found that the application of different forms of both P and N fertilizers are highly ranked for the prediction of Chl a concentration. The developed ML model successfully predicts Chl a with a coefficient of determination of 0.81, bias of -0.12 µg/l and RMSE of 4.97 µg/l. The developed ML-based modeling approach can be used for impact assessment of agriculture practices in a changing climate that affect Chl a concentrations in Lake Erie.

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.

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 PM2.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 O3 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 NOx 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.

Impacts of tiled land cover characterization in the Model for Prediction Across Scales-Atmosphere (MPAS-A).

Parameterization of subgrid-scale variability of land cover characterization (LCC) is an active area of research, and can improve model performance compared to the dominant (i.e., most abundant tile) approach. The "Noah" land surface model implementation in the global Model for Predictions Across Scales-Atmosphere (MPAS-A), however, only uses the dominant LCC approach that leads to oversimplification in regions of highly heterogeneous LCC (e.g., urban/suburban settings). Thus, in this work we implement a subgrid tiled approach as an option in MPAS-A, version 6.0, and assess the impacts of tiled LCC on meteorological predictions for two gradually refining meshes (92-25 and 46-12 km) focused on the conterminous U.S for January and July 2016. Compared to the dominant approach, results show that using the tiled LCC leads to pronounced global changes in 2-m temperature (July global average change ~ -0.4 K), 2-m moisture, and 10-m wind speed for the 92-25 km mesh. The tiled LCC reduces mean biases in 2-m temperature (July U.S. average bias reduction ~ factor of 4) and specific humidity in the central and western U.S. for the 92-25 km mesh, improves the agreement of vertical profiles (e.g., temperature, humidity, and wind speed) with observed radiosondes; however, there is increased bias and error for incoming solar radiation at the surface. The inclusion of subgrid LCC has implications for reducing systematic temperature biases found in numerical weather prediction models, particularly those that employ a dominant LCC approach.

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

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

A Measurement-Model Fusion Approach for Improved Wet Deposition Maps and Trends.

Air quality models provide spatial fields of wet deposition (WD) and dry deposition that explicitly account for the transport and transformation of emissions from thousands of sources. However, many sources of uncertainty in the air quality model including errors in emissions and meteorological inputs (particularly precipitation) and incomplete descriptions of the chemical and physical processes governing deposition can lead to bias and error in the simulation of WD. We present an approach to bias correct Community Multiscale Air Quality model output over the contiguous United States using observation-based gridded precipitation data generated by the Parameter-elevation Regressions on Independent Slopes Model and WD observations at the National Atmospheric Deposition Program National Trends Network sites. A cross-validation analysis shows that the adjusted annual accumulated WD for NO3 -, NH4 +, and SO4 2- from 2002 to 2012 has less bias and higher correlation with observed values than the base model output without adjustment. Temporal trends in observed WD are captured well by the adjusted model simulations across the entire contiguous United States. Consistent with previous trend analyses, WD NO3 - and SO4 2- are shown to decrease during this period in the eastern half of the United States, particularly in the Northeast, while remaining nearly constant in the West. Trends in WD of NH4 + are more spatially and temporally heterogeneous, with some positive trends in the Great Plains and Central Valley of CA and slightly negative trends in the south.

Updates to the Noah Land Surface Model in WRF-CMAQ to Improve Simulated Meteorology, Air Quality, and Deposition.

Regional, state, and local environmental regulatory agencies often use Eulerian models to investigate the potential impacts on pollutant deposition and air quality from changes in land use, anthropogenic and natural emissions, and climate. The Noah land surface model (LSM) in the Weather Research and Forecasting (WRF) model is widely used with the Community Multiscale Air Quality (CMAQ) model for such investigations, but there are many inconsistencies that need to be changed so that they are consistent with dry deposition and emission processes. In this work, the Noah LSM in WRFv3.8.1 is improved in its linkage to CMAQv5.2 by adding important parameters to the WRF/Noah output, updating the WRF soil and vegetation reference tables that influence CMAQ wet and dry photochemical deposition processes, and decreasing WRF/Noah's top soil layer depth to be consistent with CMAQ processes (e.g., windblown dust and bidirectional ammonia exchange). The modified WRF/Noah-CMAQ system (both off-line and coupled) impacts meteorological predictions of 2-m temperature (T2; increases and decreases), 2-m mixing ratio (Q2; decreases), and 10-m wind speed (WSPD10; decreases) in the United States. These changes are mostly driven by leaf area index values and aerodynamic roughness lengths updated in the vegetation tables based on satellite data, with additional impacts from soil tables updated based on recent soil data. Improvements in the consistency in the treatment of land surface processes between CMAQ and WRF resulted in improvements in both estimated meteorological (e.g., T2, WSPD10, and latent heat fluxes) and chemical (e.g., ozone, sulfur dioxide, and windblown dust) model estimates.