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.

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.

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 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.

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.

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.

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 NH3 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.

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.

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.