Surface Flux Modeling for Air Quality Applications

Meteorological effects on sources and future projection of nitrogen deposition to lakes in China

Lake ecosystems are extremely sensitive to nitrogen growth, which leads to water quality degradation and ecosystem health decline. Nitrogen depositions, as one of the main sources of nitrogen in water, are expected to change under future climate change scenarios. However, it remains not clear how nitrogen deposition to lakes respond to future meteorological conditions. In this study, a source-oriented version of Community Multiscale Air Quality (CMAQ) Model was used to estimate nitrogen deposition to 263 lakes in 2013 and under three RCP scenarios (4.5, 6.0 and 8.5) in 2046. Annual total deposition of 58.2 Gg nitrogen was predicted for all lakes, with 23.3 Gg N by wet deposition and 34.9 Gg N by dry deposition. Nitrate and ammonium in aerosol phase are the major forms of wet deposition, while NH3 and HNO3 in gas phase are the major forms of dry deposition. Agriculture emissions contribute to 57% of wet deposition and 44% of dry deposition. Under future meteorological conditions, wet deposition is predicted to increase by 5.5% to 16.4%, while dry deposition would decrease by 0.3% to 13.0%. Changes in wind speed, temperature, relative humidity (RH), and precipitation rates are correlated with dry and wet deposition changes. The predicted changes in deposition to lakes driven by meteorological changes can lead to significant changes in aquatic chemistry and ecosystem functions. Apart from future emission scenarios, different climate scenarios should be considered in future ecosystem health evaluation in response to nitrogen deposition.

Satellite Observations of Atmospheric Ammonia Inequalities Associated with Industrialized Swine Facilities in Eastern North Carolina

Industrialized swine facilities adversely affect the health and well-being of Eastern North Carolina residents in the U.S. and are an issue of environmental racism. Concentrated animal feeding operations (CAFOs) emit various harmful and noxious air pollutants, including ammonia (NH3). There are limited measurements of CAFO-related air quality, contributing to disputes around its severity. We use NH3 vertical column densities from the space-based Infrared Atmospheric Sounding Interferometer (IASI) to report systematic, distributive inequalities in NH3 column enhancements (ΔNH3 columns), equal to NH3 columns less an observationally determined tropospheric background. Population-weighted block group-scale ΔNH3 columns are higher by 27 ± 3% for Black and African Americans, 35 ± 3% for Hispanics and Latinos, and 49 ± 3% for American Indians compared to non-Hispanic/Latino whites in Eastern North Carolina (April-August 2016-2021). Surface winds and air temperature influence block group-scale NH3 distributions, with higher absolute NH3 inequalities for all groups on calm days and for Black and African Americans and Hispanics and Latinos on hot days, consistent with effects from NH3 volatization downfield of facilities from, e.g., manure-covered fields, particles, and other surfaces. ΔNH3 columns correspond spatially with permitted swine facilities, with residents living multiple kilometers from swine CAFOs chronically exposed to elevated NH3. Trends in NH3 columns over 2008-2023 are driven by regional-scale atmospheric processes rather than localized NH3 changes in CAFO emissions. Results are discussed in local decision-making contexts that have broad relevance for air quality issues without protective federal regulatory standards.

An updated modeling framework to simulate Los Angeles air quality - Part 1: Model development, evaluation, and source apportionment

This study describes a modeling framework, model evaluation, and source apportionment to understand the causes of Los Angeles (LA) air pollution. A few major updates are applied to the Community Multiscale Air Quality (CMAQ) model with a high spatial resolution (1 km × 1 km). The updates include dynamic traffic emissions based on real-time, on-road information and recent emission factors and secondary organic aerosol (SOA) schemes to represent volatile chemical products (VCPs). Meteorology is well predicted compared to ground-based observations, and the emission rates from multiple sources (i.e., on-road, volatile chemical products, area, point, biogenic, and sea spray) are quantified. Evaluation of the CMAQ model shows that ozone is well predicted despite inaccuracies in nitrogen oxide (NO x ) predictions. Particle matter (PM) is underpredicted compared to concurrent measurements made with an aerosol mass spectrometer (AMS) in Pasadena. Inorganic aerosol is well predicted, while SOA is underpredicted. Modeled SOA consists of mostly organic nitrates and products from oxidation of alkane-like intermediate volatility organic compounds (IVOCs) and has missing components that behave like less-oxidized oxygenated organic aerosol (LO-OOA). Source apportionment demonstrates that the urban areas of the LA Basin and vicinity are NO x -saturated (VOC-sensitive), with the largest sensitivity of O3 to changes in VOCs in the urban core. Differing oxidative capacities in different regions impact the nonlinear chemistry leading to PM and SOA formation, which is quantified in this study.

Quantifying anthropogenic emission of iron in marine aerosol in the Northwest Pacific with shipborne online measurements

Anthropogenic emissions are recognized as significant contributors to atmospheric soluble iron (Fe) in recent years, which may affect marine primary productivity, especially in Fe-limited areas. However, the contribution of different emission sources to Fe in marine aerosol has been primarily estimated by modeling approaches. Quantifying anthropogenic Fe based on field measurements remains a great challenge. In this study, online multi-element measurements and Positive Matrix Factorization (PMF) were combined for the first time to quantify sources of atmospheric Fe and soluble Fe in the Northwest Pacific during a cruise in spring 2015. Fe concentration in 624 atmospheric PM2.5 samples measured online was 74.58 ± 90.87 ng/m3. The PMF results showed anthropogenic activities, including industrial coal combustion, biomass burning, and maritime transport, were important in this region, contributing 31.4 % of atmospheric Fe on average. In addition, anthropogenic Fe concentration resolved by PMF was comparable to the simulation results of the CMAQ (Community Multiscale Air Quality) and GEOS-Chem (Goddard Earth Observing System-Chemical transport) models, with better correlation to CMAQ (r = 0.76) than GEOS-Chem (r = 0.26). This study developed a new method to estimate atmospheric soluble Fe, which integrates Fe source apportionment results and Fe solubility from different sources. Soluble Fe concentration was estimated as 3.93 ± 5.14 ng/m3, of which 87.0 % was attributed to anthropogenic emissions. Notably, ship emission alone contributed 27.5 % of soluble Fe, though its contribution to total Fe was only 2.2 %. Finally, the total deposition fluxes of atmospheric Fe (37.11 ± 38.43 μg/m2/day) and soluble Fe (1.85 ± 2.13 μg/m2/day) were estimated. This study developed a new methodology for quantifying contribution of anthropogenic emissions to Fe in marine aerosol, which could greatly help the assessment of impacts of human activities on marine environment.

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.

Ambient ozone at a rural Central European site and its vertical concentration gradient close to the ground

The representativeness of ambient air quality of an in situ measurement is key in the use and correct interpretation of the measured concentration values. Though the horizontal representativeness aspect is generally not neglected in air pollution studies, a detailed, high-resolution vertical distribution of ambient air pollutant concentrations is rarely addressed. The aim of this study is twofold: (i) to explore the vertical distribution of ground-level ozone (O3) concentrations measured at four heights above the ground-namely at 2, 8, 50, and 230 m-and (ii) to examine in detail the vertical O3 concentration gradient in air columns between 2 and 8, 8 and 50, and 50 and 230 m above the ground. We use the daily mean O3 concentrations measured continuously at the Košetice station, representing the rural Central European background ambient air quality observed during 2015-2021. We use the semiparametric GAM (generalised additive model) approach (with complexity or roughness-penalised splines implementation) to analyse the data with sufficient flexibility. Our models for both O3 concentrations and O3 gradients use (additive) decomposition into annual trend and seasonality (plus an overall intercept). The seasonal and year-to-year patterns of the modelled O3 concentrations look very similar at first glance. Nevertheless, a more detailed look through O3 gradients shows that they differ substantially with respect to their seasonal and long-term dynamics. The vertical O3 concentration gradient in 2-230 m is not uniform but changes substantially with increasing height and shows by far the highest dynamics near the ground between 2 and 8 m, differing in both the seasonal and annual aspects for all the air columns inspected. We speculate that non-linear changes of both seasonal and annual components of vertical O3 gradients are due to atmospheric-terrestrial interactions and to meteorological factors, which we will explore in a future study.

Comparison of ozone formation attribution techniques in the northeastern United States

The Integrated Source Apportionment Method (ISAM) has been revised in the Community Multiscale Air Quality (CMAQ) model. This work updates ISAM to maximize its flexibility, particularly for ozone (O3) modeling, by providing multiple attribution options, including products inheriting attribution fully from nitrogen oxide reactants, fully from volatile organic compound (VOC) reactants, equally from all reactants, or dynamically from NO x or VOC reactants based on the indicator gross production ratio of hydrogen peroxide (H2O2) to nitric acid (HNO3). The updated ISAM has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3.2 and beyond). This study's primary objective is to document these ISAM updates and demonstrate their impacts on source apportionment results for O3 and its precursors. Additionally, the ISAM results are compared with the Ozone Source Apportionment Technology (OSAT) in the Comprehensive Air-quality Model with Extensions (CAMx) and the brute-force method (BF). All comparisons are performed for a 4 km horizontal grid resolution application over the northeastern US for a selected 2 d summer case study (9 and 10 August 2018). General similarities among ISAM, OSAT, and BF results add credibility to the new ISAM algorithms. However, some discrepancies in magnitude or relative proportions among tracked sources illustrate the distinct features of each approach, while others may be related to differences in model formulation of chemical and physical processes. Despite these differences, OSAT and ISAM still provide useful apportionment data by identifying the geographical and temporal contributions of O3 and its precursors. Both OSAT and ISAM attribute the majority of O3 and NO x contributions to boundary, mobile, and biogenic sources, whereas the top three contributors to VOCs are found to be biogenic, boundary, and area sources.

Atmospheric deposition of reactive nitrogen to a deciduous forest in the southern Appalachian Mountains

Assessing nutrient critical load exceedances requires complete and accurate atmospheric deposition budgets for reactive nitrogen (Nr). The exceedance is the total amount of Nr deposited to the ecosystem in excess of the critical load, which is the amount of Nr input below which harmful effects do not occur. Total deposition includes all forms of Nr (i.e., organic and inorganic) deposited to the ecosystem by wet and dry pathways. Here we present results from the Southern Appalachian Nitrogen Deposition Study (SANDS), in which a combination of measurements and field-scale modeling was used to develop a complete annual Nr deposition budget for a deciduous forest at the Coweeta Hydrologic Laboratory. Wet deposition of ammonium, nitrate, nitrite, and bulk organic N were measured directly. The dry deposited Nr fraction was estimated using a bidirectional resistance-based model driven with speciated measurements of Nr air concentrations (e.g., ammonia, ammonium aerosol, nitric acid, nitrate aerosol, bulk organic N in aerosol, total alkyl nitrates, and total peroxy nitrates), micrometeorology, canopy structure, and biogeochemistry. Total annual deposition was ~6.7 kg N ha-1 yr-1, which is on the upper end of Nr critical load estimates recently developed for similar ecosystems in the nearby Great Smoky Mountains National Park. Of the total (wet + dry) budget, 51.1% was contributed by reduced forms of NrNH x = ammonia + ammonium ) , with oxidized and organic forms contributing ~41.3% and 7.6%, respectively. Our results indicate that reductions inNH x deposition would be needed to achieve the lowest estimates (~3.0 kg N ha-1 yr-1) of Nr critical loads in southern Appalachian forests.

Dry Deposition Methods Based on Turbulence Kinetic Energy: Part 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single-Point Model

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. To address these issues with dry deposition, firstly we take advantage of three-dimensional (3-D) turbulence information in estimating resistances by proposing and validating a 3-D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3-D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (1) present a set of resistance formulations for ozone (O3) based on the 3-D turbulence velocity scale; (2) intercompare estimations of such resistances with those obtained using existing formulations; and, (3) evaluate simulated O3 fluxes using a single-point dry deposition model against long-term observations of O3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

Spatial heterogeneity of ammonia fluxes in a deciduous forest and adjacent grassland

Gas-phase ammonia (NH3), emitted primarily from agriculture, contributes significantly to reactive nitrogen (Nr) deposition. Excess deposition of Nr to the environment causes acidification, eutrophication, and loss of biodiversity. The exchange of NH3 between land and atmosphere is bidirectional and can be highly heterogenous when underlying vegetation and soil characteristics differ. Direct measurements that assess the spatial heterogeneity of NH3 fluxes are lacking. To this end, we developed and deployed two fast-response, quantum cascade laser-based open-path NH3 sensors to quantify NH3 fluxes at a deciduous forest and an adjacent grassland separated by 700 m in North Carolina, United States from August to November, 2017. The sensors achieved 10 Hz precisions of 0.17 ppbv and 0.23 ppbv in the field, respectively. Eddy covariance calculations showed net deposition of NH3 (-7.3 ng NH3-N m-2 s-1) to the forest canopy and emission (3.2 ng NH3-N m-2 s-1) from the grassland. NH3 fluxes at both locations displayed diurnal patterns with midday peaks and smaller peaks in the afternoons. Concurrent biogeochemistry data showed over an order of magnitude higher NH3 emission potentials from green vegetation at the grassland compared to the forest, suggesting a possible explanation for the observed flux differences. Back trajectories originating from the site identified the upwind urban area as the main source region of NH3. Our work highlights that adjacent natural ecosystems sharing the same airshed but different vegetation and biogeochemical conditions may differ remarkably in NH3 exchange. Such heterogeneities should be considered when upscaling point measurements, downscaling modeled fluxes, and evaluating Nr deposition for different natural land use types in the same landscape. Additional in-situ flux measurements accompanied by comprehensive biogeochemical and micrometeorological records over longer periods are needed to fully characterize the temporal variabilities and trends of NH3 fluxes and identify the underlying driving factors.

Improving the particle dry deposition scheme in the CMAQ photochemical modeling system

Dry deposition of atmospheric aerosols in large-scale models is a critical, but highly uncertain, sink process with a strong dependence on particle size, meteorological conditions, and land surface properties. This study investigates the particle dry deposition scheme implemented in the standard Community Multiscale Air Quality (CMAQ) model v5.2.1, characterizes its underlying parameterized components with comparison to a similar scheme in a contemporary regional-scale model, and proposes two updated schemes that are then evaluated with available ambient particle deposition velocity (Vd) measurements. Both updated schemes reduce the surprisingly strong dependence of deposition velocity on the aerosol mode width, with one scheme further introducing a dependence on vegetation coverage that is broadly consistent with variability in observations between vegetated and non-vegetated surfaces. Compared to the base scheme, the updated scheme with vegetation dependence increases Vd for submicron particles and decreases it for larger particles by an average of 37% and -66%, respectively. This scheme performs statistically better than the base scheme, reducing fractional biases by 56%-97% for vegetated land-use types and has roughly equivalent performance over water. The base and updated schemes are tested with three annual CMAQ (v5.2.1) simulations for the year 2011; predicted ambient aerosol concentrations are evaluated with routine monitoring network observations and predicted dry deposition fluxes are evaluated with data from the Clean Air Status and Trends Network (CASTNET). The updated scheme with vegetation dependence reduces negative fractional biases for PM10 by 41% and positive fractional biases for PM2.5 organic carbon by 15%. This scheme has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3 and beyond) to replace the scheme used in previous versions of CMAQ (v4.5 through v5.2.1).

Dry Deposition Methods Based on Turbulence Kinetic Energy: Part 2. Extension to Particle Deposition Using a Single-Point Model

Magnitude of atmospheric turbulence, a key driver of several processes that contribute to aerosol (i.e., particle) deposition, is underrepresented in current models. Various formulations have been developed to model particle dry deposition; all these formulations typically rely on friction velocity and some use additional ad hoc factors to represent enhanced impacts of turbulence. However, none were formally linked with the three-dimensional (3-D) turbulence. Here, we propose a set of 3-D turbulence-dependent resistance formulations for particle dry deposition simulation and intercompare the performance of new resistance formulations with that obtained from using the existing formulations and measured dry deposition velocity. Turbulence parameters such as turbulence velocity scale, turbulence factor, intensity of turbulence, effective sedimentation velocity, and effective Stokes number are newly introduced into two different particle deposition schemes to improve turbulence strength representation. For an assumed particle size distribution, the newly proposed schemes predict stronger diurnal variation of particle dry deposition velocity and are comparable to corresponding measurements while existing formulations indicate large underpredictions. We also find that the incorporation of new turbulence parameters either introduced or added stronger diurnal variability to sedimentation velocity and collection efficiencies values, making the new schemes predict higher deposition values during daytime and nighttime when compared to existing schemes. The findings from this research may help improve the capability of dry deposition schemes in regional and global models.

Modeling past and future spatiotemporal distributions of airborne allergenic pollen across the contiguous United States

Exposures to airborne allergenic pollen have been increasing under the influence of changing climate. A modeling system incorporating pollen emissions and atmospheric transport and fate processes has been developed and applied to simulate spatiotemporal distributions of two major aeroallergens, oak and ragweed pollens, across the contiguous United States (CONUS) for both historical (year 2004) and future (year 2047) conditions. The transport and fate of pollen presented here is simulated using our adapted version of the Community Multiscale Air Quality (CMAQ) model. Model performance was evaluated using observed pollen counts at monitor stations across the CONUS for 2004. Our analysis shows that there is encouraging consistency between observed seasonal mean concentrations and corresponding simulated seasonal mean concentrations (oak: Pearson = 0.35, ragweed: Pearson = 0.40), and that the model was able to capture the statistical patterns of observed pollen concentration distributions in 2004 for most of the pollen monitoring stations. Simulation of pollen levels for a future year (2047) considered conditions corresponding to the RCP8.5 scenario. Modeling results show substantial regional variability both in the magnitude and directionality of changes in pollen metrics. Ragweed pollen season is estimated to start earlier and last longer for all nine climate regions of the CONUS, with increasing average pollen concentrations in most regions. The timing and magnitude of oak pollen season vary across the nine climate regions, with the largest increases in pollen concentrations expected in the Northeast region.

The sensitivities of ozone and PM2.5 concentrations to the satellite-derived leaf area index over East Asia and its neighboring seas in the WRF-CMAQ modeling system

Vegetation plays an important role as both a sink of air pollutants via dry deposition and a source of biogenic VOC (BVOC) emissions which often provide the precursors of air pollutants. To identify the vegetation-driven offset between the deposition and formation of air pollutants, this study examines the responses of ozone and PM_2.5 concentrations to changes in the leaf area index (LAI) over East Asia and its neighboring seas, using up-to-date satellite-derived LAI and green vegetation fraction (GVF) products. Two LAI scenarios that examine (1) table-prescribed LAI and GVF from 1992 to 1993 AVHRR and 2001 MODIS products and (2) reprocessed 2019 MODIS LAI and 2019 VIIRS GVF products were used in WRF-CMAQ modeling to simulate ozone and PM_2.5 concentrations for June 2019. The use of up-to-date LAI and GVF products resulted in monthly mean LAI differences ranging from -56.20% to 96.81% over the study domain. The increase in LAI resulted in the differences in hourly mean ozone and PM_2.5 concentrations over inland areas ranging from 0.27 ppbV to -7.17 ppbV and 0.89 μg/m³ to -2.65 μg/m³, and the differences of those over the adjacent sea surface ranging from 0.69 ppbV to -2.86 ppbV and 3.41 μg/m³ to -7.47 μg/m³. The decreases in inland ozone and PM_2.5 concentrations were mainly the results of dry deposition accelerated by increases in LAI, which outweighed the ozone and PM_2.5 formations via BVOC-driven chemistry. Some inland regions showed further decreases in PM_2.5 concentrations due to reduced reactions of PM_2.5 precursors with hydroxyl radicals depleted by BVOCs. The reductions in sea surface ozone and PM_2.5 concentrations were accompanied by the reductions in those in upwind inland regions, which led to less ozone and PM_2.5 inflows. The results suggest the importance of the selective use of vegetation parameters for air quality modeling.

Simulating PM2.5 Concentrations during New Year in Cuenca, Ecuador: Effects of Advancing the Time of Burning Activities

Fine particulate matter (PM2.5) is dangerous to human health. At midnight on 31 December, in Ecuadorian cities, people burn puppets and fireworks, emitting high amounts of PM2.5. On 1 January 2022, concentrations between 27.3 and 40.6 µg m-3 (maximum mean over 24 h) were measured in Cuenca, an Andean city located in southern Ecuador; these are higher than 15 µg m-3, the current World Health Organization guideline. We estimated the corresponding PM2.5 emissions and used them as an input to the Weather Research and Forecasting with Chemistry (WRF-Chem 3.2) model to simulate the change in PM2.5 concentrations, assuming these emissions started at 18:00 LT or 21:00 LT on 31 December 2021. On average, PM2.5 concentrations decreased by 51.4% and 33.2%. Similar modeling exercises were completed for 2016 to 2021, providing mean decreases between 21.4% and 61.0% if emissions started at 18:00 LT. Lower mean reductions, between 2.3% and 40.7%, or even local increases, were computed for emissions beginning at 21:00 LT. Reductions occurred through better atmospheric conditions to disperse PM2.5 compared to midnight. Advancing the burning time can help reduce the health effects of PM2.5 emissions on 31 December.

Diurnal landfill methane flux patterns across different seasons at a landfill in Southeastern US

Diurnal patterns of methane flux are examined at a landfill in the Southeastern US. Methane fluxes are measured by an eddy covariance (EC) tower during representative one-week periods in three seasons: summer, fall, and winter. Measured methane fluxes are compared with atmospheric pressure, temporal variation of atmospheric pressure, wind shear velocity, and air temperature. Landfill methane flux varies significantly with shear velocity and temporal changes in atmospheric pressure when the atmosphere is neutral. Under unstable atmospheric conditions, air temperature correlates best with methane flux, which is corroborated with an independent dataset of tracer correlation method (TCM) measurements for similar measurement periods. These field data support a mathematical model previously proposed to describe atmospheric effects on methane flux from landfills. The field data also indicate significant diurnal methane flux variations, with daytime fluxes up to 23 times greater than nighttime fluxes. Because the majority of historical TCM measurements of whole landfill methane flux are between 12 pm and 6 pm at this landfill, when daily emissions are highest because of atmospheric effects, average diurnal fluxes might have been overestimated by as much as 73%. Methane emissions are most representative of diurnal average emissions when atmospheric stability is near-neutral, which occurs in the late morning (∼11 am) and in the early evening (∼5 pm) at this site.

Development and evaluation of an advanced National Air Quality Forecasting Capability using the NOAA Global Forecast System version 16

A new dynamical core, known as the Finite-Volume Cubed-Sphere (FV3) and developed at both NASA and NOAA, is used in NOAA's Global Forecast System (GFS) and in limited-area models for regional weather and air quality applications. NOAA has also upgraded the operational FV3GFS to version 16 (GFSv16), which includes a number of significant developmental advances to the model configuration, data assimilation, and underlying model physics, particularly for atmospheric composition to weather feedback. Concurrent with the GFSv16 upgrade, we couple the GFSv16 with the Community Multiscale Air Quality (CMAQ) model to form an advanced version of the National Air Quality Forecasting Capability (NAQFC) that will continue to protect human and ecosystem health in the US. Here we describe the development of the FV3GFSv16 coupling with a "state-of-the-science" CMAQ model version 5.3.1. The GFS-CMAQ coupling is made possible by the seminal version of the NOAA-EPA Atmosphere-Chemistry Coupler (NACC), which became a major piece of the next operational NAQFC system (i.e., NACC-CMAQ) on 20 July 2021. NACC-CMAQ has a number of scientific advancements that include satellite-based data acquisition technology to improve land cover and soil characteristics and inline wildfire smoke and dust predictions that are vital to predictions of fine particulate matter (PM_2.5) concentrations during hazardous events affecting society, ecosystems, and human health. The GFS-driven NACC-CMAQ model has significantly different meteorological and chemical predictions compared to the previous operational NAQFC, where evaluation of NACC-CMAQ shows generally improved near-surface ozone and PM_2.5 predictions and diurnal patterns, both of which are extended to a 72 h (3 d) forecast with this system.

Extension of a gaseous dry deposition algorithm to oxidized volatile organic compounds and hydrogen cyanide for application in chemistry transport models

The dry deposition process refers to flux loss of an atmospheric pollutant due to uptake of the pollutant by the Earth's surfaces, including vegetation, underlying soil, and any other surface types. In chemistry transport models (CTMs), the dry deposition flux of a chemical species is typically calculated as the product of its surface layer concentration and its dry deposition velocity (V _d); the latter is a variable that needs to be highly empirically parameterized due to too many meteorological, biological, and chemical factors affecting this process. The gaseous dry deposition scheme of Zhang et al. (2003) parameterizes V _d for 31 inorganic and organic gaseous species. The present study extends the scheme of Zhang et al. (2003) to include an additional 12 oxidized volatile organic compounds (oVOCs) and hydrogen cyanide (HCN), while keeping the original model structure and formulas, to meet the demand of CTMs with increasing complexity. Model parameters for these additional chemical species are empirically chosen based on their physicochemical properties, namely the effective Henry's law constants and oxidizing capacities. Modeled V _d values are compared against field flux measurements over a mixed forest in the southeastern US during June 2013. The model captures the basic features of the diel cycles of the observed V _d. Modeled V _d values are comparable to the measurements for most of the oVOCs at night. However, modeled V _d values are mostly around 1 cm s^-1 during daytime, which is much smaller than the observed daytime maxima of 2-5 cm s^-1. Analysis of the individual resistance terms and uptake pathways suggests that flux divergence due to fast atmospheric chemical reactions near the canopy was likely the main cause of the large model-measurement discrepancies during daytime. The extended dry deposition scheme likely provides conservative V _d values for many oVOCs. While higher V _d values and bidirectional fluxes can be simulated by coupling key atmospheric chemical processes into the dry deposition scheme, we suggest that more experimental evidence of high oVOC V _d values at additional sites is required to confirm the broader applicability of the high values studied here. The underlying processes leading to high measured oVOC V _d values require further investigation.

Dry Deposition of Atmospheric Aerosols: Approaches, Observations, and Mechanisms

Aerosols are liquid or solid particles suspended in the atmosphere, typically with diameters on the order of nanometers to microns. These particles impact air quality and the radiative balance of the planet. Dry deposition is a key process for the removal of aerosols from the atmosphere and plays an important role in controlling the lifetime of atmospheric aerosols. Dry deposition is driven by turbulence and shows a strong dependence on particle size. This review summarizes the mechanisms behind aerosol dry deposition, including measurement approaches, field observations, and modeling studies. We identify several gaps in the literature, including deposition over the cryosphere (i.e., snow and ice surfaces) and the ocean; in addition, we highlight new techniques to measure black carbon fluxes. While recent advances in aerosol instrumentation have enhanced our understanding of aerosol sources and chemistry, dry deposition and other loss processes remain poorly investigated.

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.

Dry deposition of reactive nitrogen to different ecosystems across eastern China: A comparison of three community models

The assessment of nitrogen ecosystem loads mostly use the method of sampling observation combined with numerical model to estimate the spatial distribution pattern of nitrogen dry deposition flux. The selection of models is important which directly affects the reliability of the deposition flux results. In this study, the performance of three widely used models (WRF-Chem, EMEP, CMAQ) are compared. The dry deposition fluxes of typical active nitrogen components over eastern China showed uncertainties by a factor of 0.5 ~ 2 between the oxidized nitrogen (OXN) results of the three models and the observation network while the reduced nitrogen (RDN) simulation results are underestimated by a quarter of the observation reports. These three models show different results on four typical ecosystems: simulation of EMEP got the highest for OXN dry deposition flux on each ecosystem (urban 14.94 ± 4.92kgN ⋅ ha^-1 ⋅ yr^-1, cropland/grassland 5.53 ± 5.11kgN ⋅ ha^-1 ⋅ yr^-1, forest 4.75 ± 4.32kgN ⋅ ha^-1 ⋅ yr^-1, water bodies 1.48 ± 1.53kgN ⋅ ha^-1 ⋅ yr^-1); WRF-Chem has the highest value of RDN on the urban (8.91 ± 6.44kgN ⋅ ha^-1 ⋅ yr^-1) and water bodies (1.01 ± 1.44kgN ⋅ ha^-1 ⋅ yr^-1) while EMEP is highest in cropland/grassland (3.42 ± 3.43kgN ⋅ ha^-1 ⋅ yr^-1) and forest (2.34 ± 1.94kgN ⋅ ha^-1 ⋅ yr^-1). CMAQ is in medium range for both OXN and RDN simulations on each ecosystem. Compare with the critical loads, CMAQ generates more exceeded critical load areas than WRF-Chem and EMEP on cropland/grassland and forests ecosystem. For water bodies, WRF-chem and CMAQ showed higher exceeding critical load areas than EMEP. In summary, EMEP generally underestimates while the CMAQ and WRF-Chem model would overestimate the impacts on the ecosystems. So, policy implementation needs special attention accounting the difference of simulation effect with different models.

Dry Deposition of Ozone over Land: Processes, Measurement, and Modeling

Dry deposition of ozone is an important sink of ozone in near surface air. When dry deposition occurs through plant stomata, ozone can injure the plant, altering water and carbon cycling and reducing crop yields. Quantifying both stomatal and nonstomatal uptake accurately is relevant for understanding ozone's impact on human health as an air pollutant and on climate as a potent short-lived greenhouse gas and primary control on the removal of several reactive greenhouse gases and air pollutants. Robust ozone dry deposition estimates require knowledge of the relative importance of individual deposition pathways, but spatiotemporal variability in nonstomatal deposition is poorly understood. Here we integrate understanding of ozone deposition processes by synthesizing research from fields such as atmospheric chemistry, ecology, and meteorology. We critically review methods for measurements and modeling, highlighting the empiricism that underpins modeling and thus the interpretation of observations. Our unprecedented synthesis of knowledge on deposition pathways, particularly soil and leaf cuticles, reveals process understanding not yet included in widely-used models. If coordinated with short-term field intensives, laboratory studies, and mechanistic modeling, measurements from a few long-term sites would bridge the molecular to ecosystem scales necessary to establish the relative importance of individual deposition pathways and the extent to which they vary in space and time. Our recommended approaches seek to close knowledge gaps that currently limit quantifying the impact of ozone dry deposition on air quality, ecosystems, and climate.

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

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

PLAIN LANGUAGE SUMMARY

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

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

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

Characteristics of Atmospheric Reactive Nitrogen Deposition in Nyingchi City

Atmospheric reactive nitrogen (N) deposition has been proven to be an important nutrient input from external environments to forest ecosystems. However, the magnitude of atmospheric N deposition in the Tibetan region of China is not well known. In this study, multi-year (between 2005 and 2016) measurements of dry and wet N deposition were carried out in Nyingchi (NC) city, southeastern Tibet. Bulk deposition was collected by the rain gauge method; dry deposition was calculated by the inferential method, namely, multiplying ambient N concentrations by dry deposition velocity (V_d) of the N species. During the entire period, annual bulk and dry N deposition fluxes averaged 2.19 and 1.85 kg N ha^-1 yr^-1, respectively. Total N deposition fluxes (the sum of reduced and oxidized N species in dry and bulk deposition) showed an obvious increasing trend, especially for oxidized N species. Both bulk and dry N deposition showed a consistent seasonal pattern, with the highest fluxes in summer and the lowest in winter. Our findings suggest that N deposition to the urban environment in southeast Tibet has recently shifted from ammonium-dominated to nitrate-dominated conditions.

The impact of air pollutant deposition on solar energy system efficiency: An approach to estimate PV soiling effects with the Community Multiscale Air Quality (CMAQ) model

Deposition and accumulation of aerosol particles on photovoltaics (PV) panels, which is commonly referred to as "soiling of PV panels," impacts the performance of the PV energy system. It is desirable to estimate the soiling effect at different locations and times for modeling the PV system performance and devising cost-effective mitigation. This study presents an approach to estimate the soiling effect by utilizing particulate matter (PM) dry deposition estimates from air quality model simulations. The Community Multiscale Air Quality (CMAQ) modeling system used in this study was developed by the U.S. Environmental Protection Agency (U.S. EPA) for air quality assessments, rule-making, and research. Three deposition estimates based on different surface roughness length parameters assumed in CMAQ were used to illustrate the soling effect in different land-use types. The results were analyzed for three locations in the U.S. for year 2011. One urban and one suburban location in Colorado were selected because there have been field measurements of particle deposition on solar panels and analysis on the consequent soiling effect performed at these locations. The third location is a coastal city in Texas, the City of Brownsville. These three locations have distinct ambient environments. CMAQ underestimates particle deposition by 40% to 80% when compared to the field measurements at the two sites in Colorado due to the underestimations in both the ambient PM₁₀ concentration and deposition velocity. The estimated panel transmittance sensitivity due to the deposited particles is higher than the sensitivity obtained from the measurements in Colorado. The final soiling effect, which is transmittance loss, is estimated as 3.17 ± 4.20% for the Texas site, 0.45 ± 0.33%, and 0.31 ± 0.25% for the Colorado sites. Although the numbers are lower compared to the measurements in Colorado, the results are comparable with the soiling effects observed in U.S.

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.

Attributing differences in the fate of lateral boundary ozone in AQMEII3 models to physical process representations

Increasing emphasis has been placed on characterizing the contributions and the uncertainties of ozone imported from outside the US. In chemical transport models (CTMs), the ozone transported through lateral boundaries (referred to as LB ozone hereafter) undergoes a series of physical and chemical processes in CTMs, which are important sources of the uncertainty in estimating the impact of LB ozone on ozone levels at the surface. By implementing inert tracers for LB ozone, the study seeks to better understand how differing representations of physical processes in regional CTMs may lead to differences in the simulated LB ozone that eventually reaches the surface across the US. For all the simulations in this study (including WRF/CMAQ, WRF/CAMx, COSMO-CLM/CMAQ, and WRF/DEHM), three chemically inert tracers that generally represent the altitude ranges of the planetary boundary layer (BC1), free troposphere (BC2), and upper troposphere-lower stratosphere (BC3) are tracked to assess the simulated impact of LB specification. Comparing WRF/CAMx with WRF/CMAQ, their differences in vertical grid structure explain 10 %-60 % of their seasonally averaged differences in inert tracers at the surface. Vertical turbulent mixing is the primary contributor to the remaining differences in inert tracers across the US in all seasons. Stronger vertical mixing in WRF/CAMx brings more BC2 downward, leading to higher BCT (BCT = BC1+BC2+BC3) and BC2/BCT at the surface in WRF/CAMx. Meanwhile, the differences in inert tracers due to vertical mixing are partially counteracted by their difference in sub-grid cloud mixing over the southeastern US and the Gulf Coast region during summer. The process of dry deposition adds extra gradients to the spatial distribution of the differences in DM8A BCT by 5-10 ppb during winter and summer. COSMO-CLM/CMAQ and WRF/CMAQ show similar performance in inert tracers both at the surface and aloft through most seasons, which suggests similarity between the two models at process level. The largest difference is found in summer. Sub-grid cloud mixing plays a primary role in their differences in inert tracers over the southeastern US and the oceans in summer. Our analysis of the vertical profiles of inert tracers also suggests that the model differences in dry deposition over certain regions are offset by the model differences in vertical turbulent mixing, leading to small differences in inert tracers at the surface in these regions.

Spatial variation of modelled total, dry and wet nitrogen deposition to forests at global scale

Forests are an important biome that covers about one third of the global land surface and provides important ecosystem services. Since atmospheric deposition of nitrogen (N) can have both beneficial and deleterious effects, it is important to quantify the amount of N deposition to forest ecosystems. Measurements of N deposition to the numerous forest biomes across the globe are scarce, so chemical transport models are often used to provide estimates of atmospheric N inputs to these ecosystems. We provide an overview of approaches used to calculate N deposition in commonly used chemical transport models. The Task Force on Hemispheric Transport of Air Pollution (HTAP2) study intercompared N deposition values from a number of global chemical transport models. Using a multi-model mean calculated from the HTAP2 deposition values, we map N deposition to global forests to examine spatial variations in total, dry and wet deposition. Highest total N deposition occurs in eastern and southern China, Japan, Eastern U.S. and Europe while the highest dry deposition occurs in tropical forests. The European Monitoring and Evaluation Program (EMEP) model predicts grid-average deposition, but also produces deposition by land use type allowing us to compare deposition specifically to forests with the grid-average value. We found that, for this study, differences between the grid-average and forest specific could be as much as a factor of two and up to more than a factor of five in extreme cases. This suggests that consideration should be given to using forest-specific deposition for input to ecosystem assessments such as critical loads determinations.

Impacts of transportation sector emissions on future U.S. air quality in a changing climate. Part II: Air quality projections and the interplay between emissions and climate change

In Part II of this work we present the results of the downscaled offline Weather Research and Forecasting/Community Multiscale Air Quality (WRF/CMAQ) model, included in the "Technology Driver Model" (TDM) approach to future U.S. air quality projections (2046-2050) compared to a current-year period (2001-2005), and the interplay between future emission and climate changes. By 2046-2050, there are widespread decreases in future concentrations of carbon monoxide (CO), nitrogen oxides (NO_x = NO + NO₂), volatile organic compounds (VOCs), ammonia (NH₃), sulfur dioxide (SO₂), and particulate matter with an aerodynamic diameter ≤ 2.5 μm (PM_2.5) due mainly to decreasing on-road vehicle (ORV) emissions near urban centers as well as decreases in other transportation modes that include non-road engines (NRE). However, there are widespread increases in daily maximum 8-hr ozone (O₃) across the U.S., which are due to enhanced greenhouse gases (GHG) including methane (CH₄) and carbon dioxide (CO₂) under the Intergovernmental Panel on Climate Change (IPCC) A1B scenario, and isolated areas of larger reduction in transportation emissions of NO_x compared to that of VOCs over regions with VOC-limited O₃ chemistry. Other notable future changes are reduced haze and improved visibility, increased primary organic to elemental carbon ratio, decreases in PM_2.5 and its species, decreases and increases in dry deposition of SO₂ and O₃, respectively, and decreases in total nitrogen (TN) deposition. There is a tendency for transportation emission and CH₄ changes to dominate the increases in O₃, while climate change may either enhance or mitigate these increases in the west or east U.S., respectively. Climate change also decreases PM_2.5 in the future. Other variable changes exhibit stronger susceptibility to either emission (e.g., CO, NO_x, and TN deposition) or climate changes (e.g., VOC, NH₃, SO₂, and total sulfate deposition), which also have a strong dependence on season and specific U.S. regions.

A Comparison of Simulated and Field-Derived Leaf Area Index (LAI) and Canopy Height Values from Four Forest Complexes in the Southeastern USA

Vegetative leaf area is a critical input to models that simulate human and ecosystem exposure to atmospheric pollutants. Leaf area index (LAI) can be measured in the field or numerically simulated, but all contain some inherent uncertainty that is passed to the exposure assessments that use them. LAI estimates for minimally managed or natural forest stands can be particularly difficult to develop as a result of interspecies competition, age and spatial distribution. Satellite-based LAI estimates hold promise for retrospective analyses, but we must continue to rely on numerical models for alternative management analysis. Our objective for this study is to calculate and validate LAI estimates generated from the USDA Environmental Policy Impact Climate (EPIC) model (a widely used, field-scale, biogeochemical model) on four forest complexes spanning three physiographic provinces in Virginia and North Carolina. Measurements of forest composition (species and number), LAI, tree diameter, basal area, and canopy height were recorded at each site during the 2002 field season. Calibrated EPIC results show stand-level temporally resolved LAI estimates with R² values ranging from 0.69 to 0.96, and stand maximum height estimates within 20% of observation. This relatively high level of performance is attributable to EPIC's approach to the characterization of forest stand biogeochemical budgets, stand history, interspecies competition and species-specific response to local weather conditions. We close by illustrating the extension of this site-level approach to scales that could support regional air quality model simulations.

Atmospheric modeling to assess wind dependence in tracer dilution method measurements of landfill methane emissions

The short-term temporal variability of landfill methane emissions is not well understood due to uncertainty in measurement methods. Significant variability is seen over short-term measurement campaigns with the tracer dilution method (TDM), but this variability may be due in part to measurement error rather than fluctuations in the actual landfill emissions. In this study, landfill methane emissions and TDM-measured emissions are simulated over a real landfill in Delaware, USA using the Weather Research and Forecasting model (WRF) for two emissions scenarios. In the steady emissions scenario, a constant landfill emissions rate is prescribed at each model grid point on the surface of the landfill. In the unsteady emissions scenario, emissions are calculated at each time step as a function of the local surface wind speed, resulting in variable emissions over each 1.5-h measurement period. The simulation output is used to assess the standard deviation and percent error of the TDM-measured emissions. Eight measurement periods are simulated over two different days to look at different conditions. Results show that standard deviation of the TDM- measured emissions does not increase significantly from the steady emissions simulations to the unsteady emissions scenarios, indicating that the TDM may have inherent errors in its prediction of emissions fluctuations. Results also show that TDM error does not increase significantly from the steady to the unsteady emissions simulations. This indicates that introducing variability to the landfill emissions does not increase errors in the TDM at this site. Across all simulations, TDM errors range from -15% to 43%, consistent with the range of errors seen in previous TDM studies. Simulations indicate diurnal variations of methane emissions when wind effects are significant, which may be important when developing daily and annual emissions estimates from limited field data.

Evaluation and Intercomparison of Five North American Dry Deposition Algorithms at a Mixed Forest Site

To quantify differences between dry deposition algorithms commonly used in North America, five models were selected to calculate dry deposition velocity (V _d) for O₃ and SO₂ over a temperate mixed forest in southern Ontario, Canada, where a 5-year flux database had previously been developed. The models performed better in summer than in winter with correlation coefficients for hourly V _d between models and measurements being approximately 0.6 and 0.3, respectively. Differences in mean V _d values between models were on the order of a factor of 2 in both summer and winter. All models produced lower V _d values than the measurements of O₃ in summer and SO₂ in summer and winter, although the measured V _d may be biased. There was not a consistent tendency in the models to overpredict or underpredict for O₃ in winter. Several models produced magnitudes of the diel variation of V _d (O₃) comparable to the measurements, while all models produced slightly smaller diel variations than the measurements of V _d (SO₂) in summer. A few models produced larger diel variations than the measurements of V _d for O₃ and SO₂ in winter. Model differences were mainly due to different surface resistance parameterizations for stomatal and nonstomatal uptake pathways, while differences in aerodynamic and quasi-laminar resistances played only a minor role. It is recommended to use ensemble modeling results for ecosystem impact assessment studies, which provides mean values of all the used models and thus can avoid too much overestimations or underestimations.

Evaluation of the Community Multiscale Air Quality Model for Simulating Winter Ozone Formation in the Uinta Basin

The Weather Research and Forecasting (WRF) and Community Multiscale Air Quality (CMAQ) models were used to simulate a 10 day high-ozone episode observed during the 2013 Uinta Basin Winter Ozone Study (UBWOS). The baseline model had a large negative bias when compared to ozone (O₃) and volatile organic compound (VOC) measurements across the basin. Contrary to other wintertime Uinta Basin studies, predicted nitrogen oxides (NO _x ) were typically low compared to measurements. Increases to oil and gas VOC emissions resulted in O₃ predictions closer to observations, and nighttime O₃ improved when reducing the deposition velocity for all chemical species. Vertical structures of these pollutants were similar to observations on multiple days. However, the predicted surface layer VOC mixing ratios were generally found to be underestimated during the day and overestimated at night. While temperature profiles compared well to observations, WRF was found to have a warm temperature bias and too low nighttime mixing heights. Analyses of more realistic snow heat capacity in WRF to account for the warm bias and vertical mixing resulted in improved temperature profiles, although the improved temperature profiles seldom resulted in improved O₃ profiles. While additional work is needed to investigate meteorological impacts, results suggest that the uncertainty in the oil and gas emissions contributes more to the underestimation of O₃. Further, model adjustments based on a single site may not be suitable across all sites within the basin.

Emulation and Sensitivity Analysis of the Community Multiscale Air Quality Model for a UK Ozone Pollution Episode

Gaussian process emulation techniques have been used with the Community Multiscale Air Quality model, simulating the effects of input uncertainties on ozone and NO₂ output, to allow robust global sensitivity analysis (SA). A screening process ranked the effect of perturbations in 223 inputs, isolating the 30 most influential from emissions, boundary conditions (BCs), and reaction rates. Community Multiscale Air Quality (CMAQ) simulations of a July 2006 ozone pollution episode in the UK were made with input values for these variables plus ozone dry deposition velocity chosen according to a 576 point Latin hypercube design. Emulators trained on the output of these runs were used in variance-based SA of the model output to input uncertainties. Performing these analyses for every hour of a 21 day period spanning the episode and several days on either side allowed the results to be presented as a time series of sensitivity coefficients, showing how the influence of different input uncertainties changed during the episode. This is one of the most complex models to which these methods have been applied, and here, they reveal detailed spatiotemporal patterns of model sensitivities, with NO and isoprene emissions, NO₂ photolysis, ozone BCs, and deposition velocity being among the most influential input uncertainties.

A Photosynthesis-based Two-leaf Canopy Stomatal Conductance Model for Meteorology and Air Quality Modeling with WRF/CMAQ PX LSM

A coupled photosynthesis-stomatal conductance model with single layer sunlit and shaded leaf canopy scaling is developed for the Pleim-Xiu land surface model (LSM) option in the meteorology and air quality modeling system - WRF/CMAQ (Weather Research and Forecast model and Community Multiscale Air Quality model). The photosynthesis-based model for the PX LSM (PX PSN) is implemented and evaluated in a diagnostic box model that has evapotranspiration and ozone deposition components taken directly from WRF/CMAQ. We evaluate PX PSN for latent heat (LH) estimation at four FLUXNET sites with different vegetation types and landscape characteristics and at one FLUXNET site with ozone flux measurements against the simple Jarvis approach used in the current PX LSM. Overall, the PX PSN simulates LH as well as the PX Jarvis approach. The PX PSN, however, shows distinct advantages over the PX Jarvis approach on grassland that likely results from its treatment of C₃ and C₄ plants for CO₂ assimilation estimation. Simulations using Moderate Resolution Imaging Spectroradiometer (MODIS) LAI rather than LAI observations assess how the model would perform with the grid averaged data available in the Eulerian grid model (WRF/CMAQ). While MODIS LAI generally follows the seasonality of the observed LAI, it cannot capture the extreme highs and lows of the site measurements. MODIS LAI estimates degrade model performance at all sites but one site having old and tall trees. Ozone deposition velocity and ozone flux along with LH are simulated especially well by PX PSN as compared to significant PX Jarvis overestimation.

Evaluation and error apportionment of an ensemble of atmospheric chemistry transport modeling systems: multivariable temporal and spatial breakdown

Through the comparison of several regional-scale chemistry transport modeling systems that simulate meteorology and air quality over the European and North American continents, this study aims at (i) apportioning error to the responsible processes using timescale analysis, (ii) helping to detect causes of model error, and (iii) identifying the processes and temporal scales most urgently requiring dedicated investigations. The analysis is conducted within the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII) and tackles model performance gauging through measurement-to-model comparison, error decomposition, and time series analysis of the models biases for several fields (ozone, CO, SO2, NO, NO2, PM10, PM2.5, wind speed, and temperature). The operational metrics (magnitude of the error, sign of the bias, associativity) provide an overallsense of model strengths and deficiencies, while apportioning the error to its constituent parts (bias, variance, and covariance) can help assess the nature and quality of the error. Each of the error components is analyzed independently and apportioned to specific processes based on the corresponding timescale (long scale, synoptic, diurnal, and intraday) using the error apportionment technique devised in the former phases of AQMEII. The application of the error apportionment method to the AQMEII Phase 3 simulations provides several key insights. In addition to reaffirming the strong impact of model inputs (emission and boundary conditions) and poor representation of the stable boundary layer on model bias, results also highlighted the high interdependencies among meteorological and chemical variables, as well as among their errors. This indicates that the evaluation of air quality model performance for individual pollutants needs to be supported by complementary analysis of meteorological fields and chemical precursors to provide results that are more insightful from a model development perspective. This will require evaluaion methods that are able to frame the impact on error of processes, conditions, and fluxes at the surface. For example, error due to emission and boundary conditions is dominant for primary species (CO, particulate matter (PM)), while errors due to meteorology and chemistry are most relevant to secondary species, such as ozone. Some further aspects emerged whose interpretation requires additional consideration, such as the uniformity of the synoptic error being region- and model-independent, observed for several pollutants; the source of unexplained variance for the diurnal component; and the type of error caused by deposition and at which scale.

Advanced error diagnostics of the CMAQ and Chimere modelling systems within the AQMEII3 model evaluation framework

The work here complements the overview analysis of the modelling systems participating in the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) by focusing on the performance for hourly surface ozone by two modelling systems, Chimere for Europe and CMAQ for North America. The evaluation strategy outlined in the course of the three phases of the AQMEII activity, aimed to build up a diagnostic methodology for model evaluation, is pursued here and novel diagnostic methods are proposed. In addition to evaluating the "base case" simulation in which all model components are configured in their standard mode, the analysis also makes use of sensitivity simulations in which the models have been applied by altering and/or zeroing lateral boundary conditions, emissions of anthropogenic precursors, and ozone dry deposition. To help understand of the causes of model deficiencies, the error components (bias, variance, and covariance) of the base case and of the sensitivity runs are analysed in conjunction with timescale considerations and error modelling using the available error fields of temperature, wind speed, and NO x concentration. The results reveal the effectiveness and diagnostic power of the methods devised (which remains the main scope of this study), allowing the detection of the timescale and the fields that the two models are most sensitive to. The representation of planetary boundary layer (PBL) dynamics is pivotal to both models. In particular, (i) the fluctuations slower than ~ 1.5 days account for 70-85 % of the mean square error of the full (undecomposed) ozone time series; (ii) a recursive, systematic error with daily periodicity is detected, responsible for 10-20 % of the quadratic total error; (iii) errors in representing the timing of the daily transition between stability regimes in the PBL are responsible for a covariance error as large as 9 ppb (as much as the standard deviation of the network-average ozone observations in summer in both Europe and North America); (iv) the CMAQ ozone error has a weak/negligible dependence on the errors in NO2, while the error in NO2 significantly impacts the ozone error produced by Chimere; (v) the response of the models to variations of anthropogenic emissions and boundary conditions show a pronounced spatial heterogeneity, while the seasonal variability of the response is found to be less marked. Only during the winter season does the zeroing of boundary values for North America produce a spatially uniform deterioration of the model accuracy across the majority of the continent.

On the implications of aerosol liquid water and phase separation for organic aerosol mass

Organic compounds and liquid water are major aerosol constituents in the southeast United States (SE US). Water associated with inorganic constituents (inorganic water) can contribute to the partitioning medium for organic aerosol when relative humidities or organic matter to organic carbon (OM/OC) ratios are high such that separation relative humidities (SRH) are below the ambient relative humidity (RH). As OM/OC ratios in the SE US are often between 1.8 and 2.2, organic aerosol experiences both mixing with inorganic water and separation from it. Regional chemical transport model simulations including inorganic water (but excluding water uptake by organic compounds) in the partitioning medium for secondary organic aerosol (SOA) when RH > SRH led to increased SOA concentrations,· particularly at night. Water uptake to the organic phase resulted in even greater SOA concentrations as a result of a positive feedback in which water uptake increased SOA, which further increased aerosol water and organic aerosol. Aerosol properties· such as the OM/OC and hygroscopicity parameter (κorg), were captured well by the model compared with measurements during the Southern Oxidant and Aerosol Study (SOAS) 2013. Organic nitrates from monoterpene oxidation were predicted to be the least water-soluble semivolatile species in the model, but most biogenically derived semivolatile species in the Community Multiscale Air Quality (CMAQ) model were highly water soluble and expected to contribute to water-soluble organic carbon (WSOC). Organic aerosol and SOA precursors were abundant at night, but additional improvements in daytime organic aerosol are needed to close the model-measurement gap. When taking into account deviations from ideality, including both inorganic (when RH > SRH) and organic water in the organic partitioning medium reduced the mean bias in SOA for routine monitoring networks and improved model performance compared to observations from SOAS. Property updates from this work will be released in CMAQ v5.2.

Turbulent Fluxes and Pollutant Mixing during Wintertime Air Pollution Episodes in Complex Terrain

Cold air pools (CAPs) are stagnant stable air masses that form in valleys and basins in the winter. Low wintertime insolation limits convective mixing, such that pollutant concentrations can build up within the CAP when pollutant sources are present. In the western United States, wintertime CAPs often persist for days or weeks. Atmospheric models do not adequately capture the strength and evolution of CAPs. This is in part due to the limited availability of data quantifying the local turbulence during the formation, maintenance, and destruction of persistent CAPs. This paper presents observational data to quantify the turbulent mixing during two CAP episodes in Utah's Salt Lake Valley during February of 2004. Particulate matter (PM) concentration data and turbulence measurements for CAP and non-CAP time periods indicate that two distinct types of mixing scenarios occur depending on whether the CAP is dry or cloudy. Where cloudy, CAPs have enhanced vertical mixing due to top-down convection from the cloud layer. A comparison between the heat and momentum fluxes during 5 days of a dry CAP episode in February to those of an equivalent 5 day time period in March with no CAP indicates that the average turbulent kinetic energy during the CAP was suppressed by approximately 80%.

Preliminary evaluation of a regional atmospheric chemical data assimilation system for environmental surveillance

We report the progress of an ongoing effort by the Air Resources Laboratory, NOAA to build a prototype regional Chemical Analysis System (ARLCAS). The ARLCAS focuses on providing long-term analysis of the three dimensional (3D) air-pollutant concentration fields over the continental U.S. It leverages expertise from the NASA Earth Science Division-sponsored Air Quality Applied Science Team (AQAST) for the state-of-science knowledge in atmospheric and data assimilation sciences. The ARLCAS complies with national operational center requirement protocols and aims to have the modeling system to be maintained by a national center. Meteorology and chemistry observations consist of land-, air- and space-based observed and quality-assured data. We develop modularized testing to investigate the efficacies of the various components of the ARLCAS. The sensitivity testing of data assimilation schemes showed that with the increment of additional observational data sets, the accuracy of the analysis chemical fields also increased incrementally in varying margins. The benefit is especially noted for additional data sets based on a different platform and/or a different retrieval algorithm. We also described a plan to apply the analysis chemical fields in environmental surveillance at the Centers for Disease Control and Prevention.

Bidirectional exchange of biogenic volatiles with vegetation: emission sources, reactions, breakdown and deposition

Biogenic volatile organic compound (BVOC) emissions are widely modelled as inputs to atmospheric chemistry simulations. However, BVOC may interact with cellular structures and neighbouring leaves in a complex manner during volatile diffusion from the sites of release to leaf boundary layer and during turbulent transport to the atmospheric boundary layer. Furthermore, recent observations demonstrate that the BVOC emissions are bidirectional, and uptake and deposition of BVOC and their oxidation products are the rule rather than the exception. This review summarizes current knowledge of within-leaf reactions of synthesized volatiles with reactive oxygen species (ROS), uptake, deposition and storage of volatiles, and their oxidation products as driven by adsorption on leaf surface and solubilization and enzymatic detoxification inside leaves. The available evidence indicates that because of the reactions with ROS and enzymatic metabolism, the BVOC gross production rates are much larger than previously thought. The degree to which volatiles react within leaves and can be potentially taken up by vegetation depends upon compound reactivity, physicochemical characteristics, as well as upon their participation in leaf metabolism. We argue that future models should be based upon the concept of bidirectional BVOC exchange and consider modification of BVOC sink/source strengths by within-leaf metabolism and storage.

Toward verifying fossil fuel CO2 emissions with the CMAQ model: motivation, model description and initial simulation

Motivated by the question of whether and how a state-of-the-art regional chemical transport model (CTM) can facilitate characterization of CO2 spatiotemporal variability and verify CO2 fossil-fuel emissions, we for the first time applied the Community Multiscale Air Quality (CMAQ) model to simulate CO2. This paper presents methods, input data, and initial results for CO2 simulation using CMAQ over the contiguous United States in October 2007. Modeling experiments have been performed to understand the roles of fossil-fuel emissions, biosphere-atmosphere exchange, and meteorology in regulating the spatial distribution of CO2 near the surface over the contiguous United States. Three sets of net ecosystem exchange (NEE) fluxes were used as input to assess the impact of uncertainty of NEE on CO2 concentrations simulated by CMAQ. Observational data from six tall tower sites across the country were used to evaluate model performance. In particular, at the Boulder Atmospheric Observatory (BAO), a tall tower site that receives urban emissions from Denver CO, the CMAQ model using hourly varying, high-resolution CO2 fossil-fuel emissions from the Vulcan inventory and Carbon Tracker optimized NEE reproduced the observed diurnal profile of CO2 reasonably well but with a low bias in the early morning. The spatial distribution of CO2 was found to correlate with NO(x), SO2, and CO, because of their similar fossil-fuel emission sources and common transport processes. These initial results from CMAQ demonstrate the potential of using a regional CTM to help interpret CO2 observations and understand CO2 variability in space and time. The ability to simulate a full suite of air pollutants in CMAQ will also facilitate investigations of their use as tracers for CO2 source attribution. This work serves as a proof of concept and the foundation for more comprehensive examinations of CO2 spatiotemporal variability and various uncertainties in the future.

IMPLICATIONS

Atmospheric CO2 has long been modeled and studied on continental to global scales to understand the global carbon cycle. This work demonstrates the potential of modeling and studying CO2 variability at fine spatiotemporal scales with CMAQ, which has been applied extensively, to study traditionally regulated air pollutants. The abundant observational records of these air pollutants and successful experience in studying and reducing their emissions may be useful for verifying CO2 emissions. Although there remains much more to further investigate, this work opens up a discussion on whether and how to study CO2 as an air pollutant.

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.