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.

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.

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, PM_2.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 NH₃ accounting for NH₃ bi-directional flux, using CrIS satellite NH₃ observations over Europe in 2016. We find posterior NH₃ 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 NH₃ 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 NH₃ 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 NH₃ emissions from both schemes generally improve the magnitude and seasonality of simulated surface NH₃ 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 NH₃. This contrast highlights the need for accurately simulating diurnal variability of NH₃ 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 NH₃ emissions estimates associated with treatment of NH₃ surface exchange.

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.

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 (¹³⁷Cs) 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 ¹³⁷Cs 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 ¹³⁷Cs 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 ¹³⁷Cs 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 ¹³⁷Cs re-emission.

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 Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module in the Community Multiscale Air Quality (CMAQ) modeling system version 5.3.2

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

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

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

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.

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

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.

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

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

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.

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.

Projections of Atmospheric Nitrogen Deposition to the Chesapeake Bay Watershed

Atmospheric deposition is among the largest pathways of nitrogen loading to the Chesapeake Bay Watershed (CBW). The interplay between future climate and emission changes in and around the CBW will likely shift the future nutrient deposition abundance and chemical regime (e.g., oxidized vs. reduced nitrogen). In this work, a Representative Concentration Pathway (RCP) from the Community Earth System Model is dynamically downscaled using the Weather Research and Forecasting (WRF) and Community Multiscale Air Quality (CMAQ) model coupled to the agro-economic Environmental Policy Integrated Climate (EPIC) model. The relative impacts of emission and climate changes on atmospheric nutrient deposition are explored for a recent historical period and a period centered on 2050. The projected regional emissions in CMAQ reflect current federal and state regulations, which use baseline and projected emission years 2011 and 2040, respectively. The historical simulations of 2-m temperature and precipitation have cool and dry biases, and temperature and precipitation are projected to both increase. Ammonium wet deposition agrees well with observations, but nitrate wet deposition is underpredicted. Climate and deposition changes increase simulated future ammonium fertilizer application. In the CBW at 2050, these changes (along with widespread decreases in anthropogenic nitrogen oxide and sulfur oxide emissions, and relatively constant NH3 emissions) decrease total nitrogen deposition by 21%, decrease annual average oxidized nitrogen deposition by 44%, and increase reduced nitrogen deposition by 10%. These results emphasize the importance of decreased anthropogenic emissions on the control of future nitrogen loading to the Chesapeake Bay in a changing climate.

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.

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.

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.

Long-term trends in total inorganic nitrogen and sulfur deposition in the U.S. from 1990 to 2010

Excess deposition (including both wet and dry deposition) of nitrogen and sulfur are detrimental to ecosystems. Recent studies have investigated the spatial patterns and temporal trends of nitrogen and sulfur wet deposition, but few studies have focused on dry deposition due to the scarcity of dry deposition measurements. Here, we use long-term model simulations from the coupled WRF-CMAQ model covering the period from 1990 to 2010 to study changes in spatial distribution as well as temporal trends in total (TDEP), wet (WDEP) and dry deposition (DDEP) of total inorganic nitrogen (TIN) and sulfur (TSO4). We first evaluate the model's performance in simulating WDEP over the U.S. by comparing the model results with observational data from the U.S. National Atmospheric Deposition Program. The coupled model generally underestimates the WDEP of both TIN (including both the oxidized nitrogen deposition-TNO3, and the reduced nitrogen deposition-NHX) and TSO4, with better performance in the eastern U.S. than the western U.S. TDEP of both TIN and TSO4 show significant decreases over the U.S., especially in the east due to the large emission reductions that occurred in that region. The decreasing trends of TIN TDEP are caused by decrease of TNO3, and the increasing trends of TIN deposition over the Great Plains and Tropical Wet Forests regions are caused by increases in NH3 emissions although it should be noted that these increasing trends are not significant. TIN WDEP shows decreasing trends throughout the U.S., except for the Marine West Coast Forest region. TIN DDEP shows significant decreasing trends in the region of Eastern Temperate Forests, Northern Forests, Mediterranean California and Marine West Coast Forest, and significant increasing trends in the region of Tropical Wet Forests, Great Plains and Southern Semi-arid Highlands. For the other three regions (North American Deserts, Temperate Sierras and Northwestern Forested Mountains), the decreasing or increasing trends were not significant. Both the WDEP and DDEP of TSOx have decreases across the U.S., with a larger decreasing trend in the DDEP than that in the WDEP. Across the U.S. during the 1990-2010 period, DDEP of TIN accounted for 58-65% of TDEP of TIN. TDEP of TIN over the U.S. was dominated by deposition of TNO3 during the first decade, which then shifts to reduced nitrogen (NHX) dominance after 2003 resulting from combination of NOx emission reductions and NH3 emission increases. The sulfur DDEP is usually higher than the sulfur WDEP until recent years, as the sulfur DDEP has a larger decreasing trend than WDEP.

Modeling NH4NO3 Over the San Joaquin Valley During the 2013 DISCOVER-AQ Campaign

The San Joaquin Valley (SJV) of California experiences high concentrations of particulate matter NH4NO3 during episodes of meteorological stagnation in winter. A rich data set of observations related to NH4NO3 formation was acquired during multiple periods of elevated NH4NO3 during the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) field campaign in SJV in January and February 2013. Here NH4NO3 is simulated during the SJV DISCOVER-AQ study period with the Community Multiscale Air Quality (CMAQ) model, diagnostic model evaluation is performed using the DISCOVER-AQ data set, and integrated reaction rate analysis is used to quantify HNO3 production rates. Simulated NO3- generally agrees well with routine monitoring of 24-hr average NO3-, but comparisons with hourly average NO3- measurements in Fresno revealed differences at higher time resolution. Predictions of gas-particle partitioning of total nitrate (HNO3 + NO3-) and NHx (NH3 + NH4+) generally agree well with measurements in Fresno, although partitioning of total nitrate to HNO3 is sometimes overestimated at low relative humidity in afternoon. Gas-particle partitioning results indicate that NH4NO3 formation is limited by HNO3 availability in both the model and ambient. NH3 mixing ratios are underestimated, particularly in areas with large agricultural activity, and additional work on the spatial allocation of NH3 emissions is warranted. During a period of elevated NH4NO3, the model predicted that the OH + NO2 pathway contributed 46% to total HNO3production in SJV and the N2O5 heterogeneous hydrolysis pathway contributed 54%. The relative importance of the OH + NO2 pathway for HNO3 production is predicted to increase as NOx emissions decrease.

Long-term trends in total inorganic nitrogen and sulfur deposition in the U.S. from 1990 to 2010

Excess deposition (including both wet and dry deposition) of nitrogen and sulfur are detrimental to ecosystems. Recent studies have investigated the spatial patterns and temporal trends of nitrogen and sulfur wet deposition, but few studies have focused on dry deposition due to the scarcity of dry deposition measurements. Here, we use long-term model simulations from the coupled WRF-CMAQ model covering the period from 1990 to 2010 to study changes in spatial distribution as well as temporal trends in total (TDEP), wet (WDEP) and dry deposition (DDEP) of total inorganic nitrogen (TIN) and sulfur (TSO₄). We first evaluate the model's performance in simulating WDEP over the U.S. by comparing the model results with observational data from the U.S. National Atmospheric Deposition Program. The coupled model generally underestimates the WDEP of both TIN (including both the oxidized nitrogen deposition-TNO₃, and the reduced nitrogen deposition-NH_X) and TSO₄, with better performance in the eastern U.S. than the western U.S. TDEP of both TIN and TSO₄ show significant decreases over the U.S., especially in the east due to the large emission reductions that occurred in that region. The decreasing trends of TIN TDEP are caused by decrease of TNO₃, and the increasing trends of TIN deposition over the Great Plains and Tropical Wet Forests regions are caused by increases in NH₃ emissions although it should be noted that these increasing trends are not significant. TIN WDEP shows decreasing trends throughout the U.S., except for the Marine West Coast Forest region. TIN DDEP shows significant decreasing trends in the region of Eastern Temperate Forests, Northern Forests, Mediterranean California and Marine West Coast Forest, and significant increasing trends in the region of Tropical Wet Forests, Great Plains and Southern Semi-arid Highlands. For the other three regions (North American Deserts, Temperate Sierras and Northwestern Forested Mountains), the decreasing or increasing trends were not significant. Both the WDEP and DDEP of TSO_x have decreases across the U.S., with a larger decreasing trend in the DDEP than that in the WDEP. Across the U.S. during the 1990-2010 period, DDEP of TIN accounted for 58-65% of TDEP of TIN. TDEP of TIN over the U.S. was dominated by deposition of TNO₃ during the first decade, which then shifts to reduced nitrogen (NH_X) dominance after 2003 resulting from combination of NO_x emission reductions and NH₃ emission increases. The sulfur DDEP is usually higher than the sulfur WDEP until recent years, as the sulfur DDEP has a larger decreasing trend than WDEP.

Long-term trends in the ambient PM2.5- and O3-related mortality burdens in the United States under emission reductions from 1990 to 2010

Concentrations of both fine particulate matter (PM2.5) and ozone (O3) in the United States (US) have decreased significantly since 1990, mainly because of air quality regulations. Exposure to these air pollutants is associated with premature death. Here we quantify the annual mortality burdens from PM2.5 and O3 in the US from 1990 to 2010, estimate trends and inter-annual variability, and evaluate the contributions to those trends from changes in pollutant concentrations, population, and baseline mortality rates. We use a fine-resolution (36 km) self-consistent 21-year simulation of air pollutant concentrations in the US from 1990 to 2010, a health impact function, and annual county-level population and baseline mortality rate estimates. From 1990 to 2010, the modeled population-weighted annual PM2.5 decreased by 39 %, and summertime (April to September) 1 h average daily maximum O3 decreased by 9 % from 1990 to 2010. The PM2.5-related mortality burden from ischemic heart disease, chronic obstructive pulmonary disease, lung cancer, and stroke steadily decreased by 54% from 123 700 deaths year-1 (95% confidence interval, 70 800-178 100) in 1990 to 58 600 deaths year-1 (24 900-98 500) in 2010. The PM2.5-related mortality burden would have decreased by only 24% from 1990 to 2010 if the PM2.5 concentrations had stayed at the 1990 level, due to decreases in baseline mortality rates for major diseases affected by PM2.5. The mortality burden associated with O3 from chronic respiratory disease increased by 13% from 10 900 deaths year-1 (3700-17 500) in 1990 to 12 300 deaths year-1 (4100-19 800) in 2010, mainly caused by increases in the baseline mortality rates and population, despite decreases in O3 concentration. The O3-related mortality burden would have increased by 55% from 1990 to 2010 if the O3 concentrations had stayed at the 1990 level. The detrended annual O3 mortality burden has larger inter-annual variability (coefficient of variation of 12%) than the PM2.5-related burden (4%), mainly from the inter-annual variation of O3 concentration. We conclude that air quality improvements have significantly decreased the mortality burden, avoiding roughly 35 800 (38%) PM2.5-related deaths and 4600 (27%) O3-related deaths in 2010, compared to the case if air quality had stayed at 1990 levels (at 2010 baseline mortality rates and population).

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.

Persistence of initial conditions in continental scale air quality simulations

This study investigates the effect of initial conditions (IC) for pollutant concentrations in the atmosphere and soil on simulated air quality for two continental-scale Community Multiscale Air Quality (CMAQ) model applications. One of these applications was performed for springtime and the second for summertime. Results show that a spin-up period of ten days commonly used in regional-scale applications may not be sufficient to reduce the effects of initial conditions to less than 1% of seasonally-averaged surface ozone concentrations everywhere while 20 days were found to be sufficient for the entire domain for the spring case and almost the entire domain for the summer case. For the summer case, differences were found to persist longer aloft due to circulation of air masses and even a spin-up period of 30 days was not sufficient to reduce the effects of ICs to less than 1% of seasonally-averaged layer 34 ozone concentrations over the southwestern portion of the modeling domain. Analysis of the effect of soil initial conditions for the CMAQ bidirectional NH₃ exchange model shows that during springtime they can have an important effect on simulated inorganic aerosols concentrations for time periods of one month or longer. The effects are less pronounced during other seasons. The results, while specific to the modeling domain and time periods simulated here, suggest that modeling protocols need to be scrutinized for a given application and that it cannot be assumed that commonly-used spin-up periods are necessarily sufficient to reduce the effects of initial conditions on model results to an acceptable level. What constitutes an acceptable level of difference cannot be generalized and will depend on the particular application, time period and species of interest. Moreover, as the application of air quality models is being expanded to cover larger geographical domains and as these models are increasingly being coupled with other modeling systems to better represent air-surface-water exchanges, the effects of model initialization in such applications needs to be studied in future work.

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.

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.

Modeling the Current and Future Roles of Particulate Organic Nitrates in the Southeastern United States

Organic nitrates are an important aerosol constituent in locations where biogenic hydrocarbon emissions mix with anthropogenic NOx sources. While regional and global chemical transport models may include a representation of organic aerosol from monoterpene reactions with nitrate radicals (the primary source of particle-phase organic nitrates in the Southeast United States), secondary organic aerosol (SOA) models can underestimate yields. Furthermore, SOA parametrizations do not explicitly take into account organic nitrate compounds produced in the gas phase. In this work, we developed a coupled gas and aerosol system to describe the formation and subsequent aerosol-phase partitioning of organic nitrates from isoprene and monoterpenes with a focus on the Southeast United States. The concentrations of organic aerosol and gas-phase organic nitrates were improved when particulate organic nitrates were assumed to undergo rapid (τ = 3 h) pseudohydrolysis resulting in nitric acid and nonvolatile secondary organic aerosol. In addition, up to 60% of less oxidized-oxygenated organic aerosol (LO-OOA) could be accounted for via organic nitrate mediated chemistry during the Southern Oxidants and Aerosol Study (SOAS). A 25% reduction in nitrogen oxide (NO + NO2) emissions was predicted to cause a 9% reduction in organic aerosol for June 2013 SOAS conditions at Centreville, Alabama.

Differences between magnitudes and health impacts of BC emissions across the United States using 12 km scale seasonal source apportionment

Recent assessments have analyzed the health impacts of PM2.5 from emissions from different locations and sectors using simplified or reduced-form air quality models. Here we present an alternative approach using the adjoint of the Community Multiscale Air Quality (CMAQ) model, which provides source-receptor relationships at highly resolved sectoral, spatial, and temporal scales. While damage resulting from anthropogenic emissions of BC is strongly correlated with population and premature death, we found little correlation between damage and emission magnitude, suggesting that controls on the largest emissions may not be the most efficient means of reducing damage resulting from anthropogenic BC emissions. Rather, the best proxy for locations with damaging BC emissions is locations where premature deaths occur. Onroad diesel and nonroad vehicle emissions are the largest contributors to premature deaths attributed to exposure to BC, while onroad gasoline emissions cause the highest deaths per amount emitted. Emissions in fall and winter contribute to more premature deaths (and more per amount emitted) than emissions in spring and summer. Overall, these results show the value of the high-resolution source attribution for determining the locations, seasons, and sectors for which BC emission controls have the most effective health benefits.

Towards a climate-dependent paradigm of ammonia emission and deposition

Existing descriptions of bi-directional ammonia (NH3) land-atmosphere exchange incorporate temperature and moisture controls, and are beginning to be used in regional chemical transport models. However, such models have typically applied simpler emission factors to upscale the main NH3 emission terms. While this approach has successfully simulated the main spatial patterns on local to global scales, it fails to address the environment- and climate-dependence of emissions. To handle these issues, we outline the basis for a new modelling paradigm where both NH3 emissions and deposition are calculated online according to diurnal, seasonal and spatial differences in meteorology. We show how measurements reveal a strong, but complex pattern of climatic dependence, which is increasingly being characterized using ground-based NH3 monitoring and satellite observations, while advances in process-based modelling are illustrated for agricultural and natural sources, including a global application for seabird colonies. A future architecture for NH3 emission-deposition modelling is proposed that integrates the spatio-temporal interactions, and provides the necessary foundation to assess the consequences of climate change. Based on available measurements, a first empirical estimate suggests that 5°C warming would increase emissions by 42 per cent (28-67%). Together with increased anthropogenic activity, global NH3 emissions may increase from 65 (45-85) Tg N in 2008 to reach 132 (89-179) Tg by 2100.

Sensitivity of continental United States atmospheric budgets of oxidized and reduced nitrogen to dry deposition parametrizations

Reactive nitrogen (Nr) is removed by surface fluxes (air-surface exchange) and wet deposition. The chemistry and physics of the atmosphere result in a complicated system in which competing chemical sources and sinks exist and impact that removal. Therefore, uncertainties are best examined with complete regional chemical transport models that simulate these feedbacks. We analysed several uncertainties in regional air quality model resistance analogue representations of air-surface exchange for unidirectional and bi-directional fluxes and their effect on the continental Nr budget. Model sensitivity tests of key parameters in dry deposition formulations showed that uncertainty estimates of continental total nitrogen deposition are surprisingly small, 5 per cent or less, owing to feedbacks in the chemistry and rebalancing among removal pathways. The largest uncertainties (5%) occur with the change from a unidirectional to a bi-directional NH3 formulation followed by uncertainties in bi-directional compensation points (1-4%) and unidirectional aerodynamic resistance (2%). Uncertainties have a greater effect at the local scale. Between unidirectional and bi-directional formulations, single grid cell changes can be up to 50 per cent, whereas 84 per cent of the cells have changes less than 30 per cent. For uncertainties within either formulation, single grid cell change can be up to 20 per cent, but for 90 per cent of the cells changes are less than 10 per cent.

Estimation of in-canopy ammonia sources and sinks in a fertilized Zea mays field

An analytical model was developed to describe in-canopy vertical distribution of ammonia (NH(3)) sources and sinks and vertical fluxes in a fertilized agricultural setting using measured in-canopy mean NH(3) concentration and wind speed profiles. This model was applied to quantify in-canopy air-surface exchange rates and above-canopy NH(3) fluxes in a fertilized corn (Zea mays) field. Modeled air-canopy NH(3) fluxes agreed well with independent above-canopy flux estimates. Based on the model results, the urea fertilized soil surface was a consistent source of NH(3) one month following the fertilizer application, whereas the vegetation canopy was typically a net NH(3) sink with the lower portion of the canopy being a constant sink. The model results suggested that the canopy was a sink for some 70% of the estimated soil NH(3) emissions. A logical conclusion is that parametrization of within-canopy processes in air quality models are necessary to explore the impact of agricultural field level management practices on regional air quality. Moreover, there are agronomic and environmental benefits to timing liquid fertilizer applications as close to canopy closure as possible. Finally, given the large within-canopy mean NH(3) concentration gradients in such agricultural settings, a discussion about the suitability of the proposed model is also presented.

A note on elevated total gaseous mercury concentrations downwind from an agriculture field during tilling

Elevated mercury concentrations were measured at the University of Connecticut's mercury forest flux tower during spring agricultural field operations on an adjacent corn field. Concentrations at the tower were elevated, a peak of 7.03 ng m(-3) over the background concentration of 1.74+/-0.26 ng m(-3), during times when the prevailing wind was from the direction of the corn field and during periods when the soil was disturbed by tilling. Strong deposition to the forest was recorded at the point of measurement when atmospheric mercury concentrations were elevated. The strongest deposition rate was a 1 hour maximum of -4011 ng m(-2) h(-1) following the initial peak in atmospheric concentrations, Analyses of the meteorological conditions and mercury content in agricultural soil, manure and the diesel consumed in the tilling operation indicate that the source of the mercury was from the agricultural tilling operations and it was advected over the tower enriching the atmospheric concentrations above the forest canopy leading to deposition. These results indicate that agriculture operations resulting in a disturbed soil surface may be a source of atmospheric mercury originating from the pool of mercury bound in the soil. This represents a previously undocumented source of mercury emissions resulting from anthropogenic activities.