Articles | Volume 22, issue 19
Atmos. Chem. Phys., 22, 12749–12767, 2022
Atmos. Chem. Phys., 22, 12749–12767, 2022
Research article
30 Sep 2022
Research article | 30 Sep 2022

Long-term regional trends of nitrogen and sulfur deposition in the United States from 2002 to 2017

Long-term regional trends of nitrogen and sulfur deposition in the United States from 2002 to 2017
Sarah E. Benish1,a, Jesse O. Bash1, Kristen M. Foley1, K. Wyat Appel1, Christian Hogrefe1, Robert Gilliam1, and George Pouliot1 Sarah E. Benish et al.
  • 1US Environmental Protection Agency, Research Triangle Park, NC 27711, USA
  • aformerly at: Oak Ridge Institute for Science and Education (ORISE), US Environmental Protection Agency, Research Triangle Park, NC 27711, USA

Correspondence: Sarah E. Benish ( and Jesse O. Bash (


Atmospheric deposition of nitrogen (N) and sulfur (S) compounds from human activity has greatly declined in the United States (US) over the past several decades in response to emission controls set by the Clean Air Act. While many observational studies have investigated spatial and temporal trends of atmospheric deposition, modeling assessments can provide useful information over areas with sparse measurements, although they usually have larger horizontal resolutions and are limited by input data availability. In this analysis, we evaluate wet, dry, and total N and S deposition from multiyear simulations within the contiguous US (CONUS). Community Multiscale Air Quality (CMAQ) model estimates from the EPA's (Environmental Protection Agency) Air QUAlity TimE Series (EQUATES) project contain important model updates to atmospheric deposition algorithms compared to previous model data, including the new Surface Tiled Aerosol and Gaseous Exchange (STAGE) bidirectional deposition model which contains land-use-specific resistance parameterization and land-use-specific deposition estimates needed to estimate the differential impacts of N deposition to different land use types. First, we evaluate model estimates of wet deposition and ambient concentrations, finding underestimates of SO4, NO3, and NH4 wet deposition compared to National Atmospheric Deposition Program observations and underestimates of NH4 and SO4 and overestimates of SO2 and TNO3 (HNO3+NO3) compared to the Clean Air Status and Trends Network (CASTNET) ambient concentrations. Second, a measurement–model fusion approach employing a precipitation and bias correction to wet-deposition estimates is found to reduce model bias and improve correlations compared to the unadjusted model values. Model agreement of wet deposition is poor over parts of the West and Northern Rockies, due to errors in precipitation estimates caused by complex terrain and uncertainty in emissions at the relatively coarse 12 km grid resolution used in this study. Next, we assess modeled N and S deposition trends across climatologically consistent regions in the CONUS. Total deposition of N and S in the eastern US is larger than the western US with a steeper decreasing trend from 2002–2017; i.e., total N declined at a rate of approximately −0.30 kg N ha−1 yr−1 in the Northeast and Southeast and by −0.02 kg N ha−1 yr−1 in the Northwest and Southwest. Widespread increases in reduced N deposition across the Upper Midwest, Northern Rockies, and West indicate evolving atmospheric composition due to increased precipitation amounts over some areas, growing agricultural emissions, and regional NOx/SOx emission reductions shifting gas–aerosol partitioning; these increases in reduced N deposition are generally masked by the larger decreasing oxidized N trend. We find larger average declining trends of total N and S deposition between 2002–2009 than 2010–2017, suggesting a slowdown of the rate of decline likely in response to smaller emission reductions. Finally, we document changes in the modeled total N and S deposition budgets. The average annual total N deposition budget over the CONUS decreases from 7.8 in 2002 to 6.3 kg N ha−1 yr−1 in 2017 due to declines in oxidized N deposition from NOx emission controls. Across the CONUS during the 2002–2017 time period, the average contribution of dry deposition to the total N deposition budget drops from 60 % to 52 %, whereas wet deposition dominates the S budget rising from 45 % to 68 %. Our analysis extends upon the literature documenting the growing contribution of reduced N to the total deposition budget, particularly in the Upper Midwest and Northern Rockies, and documents a slowdown of the declining oxidized N deposition trend, which may have consequences on vegetation diversity and productivity.

1 Introduction

Human activity doubled the amount of reactive nitrogen (Nr) in the environmental globally over the past century (Fowler et al., 2013). Major sources responsible for the increase include fossil fuel combustion from vehicles and electric utilities emitting nitrogen oxides (NOx=NO + NO2) and sulfur dioxide (SO2), as well as agricultural activities releasing ammonia (NH3) (Galloway et al., 2008). After entering the atmosphere, the major nitrogen (N) and sulfur (S) removal pathways occur by precipitation (wet deposition) or uptake by surfaces, such as terrestrial and aquatic vegetation (dry deposition). Consequences of deposition include risks to public health (Townsend et al., 2003) and damages to ecosystems, such as acidification of soil and waterways, changes to species composition (Bobbink et al., 2010), and massive cyanobacteria and algal blooms (Jaworski, 1990), which may cause physiological stress and compromise immune function among wildlife (Refsnider et al., 2021). Despite adverse effects of deposition to humans and the environment, the total (wet + dry) N and S budgets remain uncertain, particularly for sensitive ecosystems in critical regions due to observational gaps.

Deposition monitoring and assessment played a decisive role informing the United States (US) Clean Air Act Amendments (CAAA) of 1990 (United States Congress, 1990). Wet deposition, sampled weekly in rain or snow, has been measured by the National Atmospheric Deposition Program (NADP) National Trends Network (NTN) since 1978. Dry-deposition measurements are inferred by combining weekly measured concentrations from the Clean Air Status and Trends Network (CASTNET) and modeled deposition velocities from the multilayer model (MLM) (Meyers et al., 1998). Dry-deposition modeling is still uncertain, particularly for land use and dry-deposition schemes in models and emission data. Despite providing critical deposition information, the limited number of NADP and CASTNET sites in essential locations, such as areas with complex terrain, near urban centers, at high elevation, or in forest ecosystems, restrict a thorough understanding on the amount and consequences of deposition. For instance, strong concentration gradients in N deposition have been documented over urban areas such as Boston, Massachusetts (Rao et al., 2013), and near Baltimore, Maryland (Bettez et al., 2013), and along coastlines like the Chesapeake Bay (Loughner et al., 2016). Additionally, networks are unable to capture the full budget of Nr. For instance, organic N species are estimated to contribute ∼25 % to total Nr (Jickells et al., 2013) but are not easily quantified by the NTN due to limitations in field and laboratory techniques (Walker et al., 2012).

Chemical transport models (CTMs) can be used to study deposition relationships and trends for locations without measurements and compounds that are challenging to quantify (e.g., organic N). However, to provide reliable estimates, CTMs estimates must first be evaluated with measurement data. Various studies have compared N and S wet deposition and concentrations estimated by the Community Multiscale Air Quality (CMAQ) model (, last access: 11 March 2022) with observed values in the US (Xing et al., 2015; Zhang et al., 2018; Appel et al., 2011), providing a wide range of agreement with measurements. Zhang et al. (2018) found adjusting CMAQv5.0.2 wet-deposition estimates with observed precipitation results in a larger negative normalized mean bias (NMB) than the original CMAQ estimates for TNO3 (NO3+HNO3, −31.6 % unadjusted vs. −35.6 % adjusted), NHx (NH3+NH4, −30.9 % unadjusted vs. −35.1 % adjusted), and total S (−5.1 % unadjusted vs. −10.5 % adjusted). Similar model agreement was published using CMAQv4.7 (Appel et al., 2011) and attributed to the coarse grid resolution (36 km). In Europe, Theobald et al. (2019) characterized the model bias of wet deposition of NOx, NHx, and SOx from six CTMs between 1990 and 2010 to range from 30 %–40 %. Biased model simulations may lead to incorrect scientific conclusions, such as critical-load designations (Williams et al., 2017). Methods to fuse model estimates with measured values can therefore be used to provide more reliable insights into “biogeochemical cycles and assessments of ecosystems and human health effects” (World Meteorological Organization, 2017). As such, the NADP Total Deposition Science Committee (TDep; see, last access: 11 March 2022) advances methods to improve estimates of atmospheric deposition from CMAQ (see Schwede and Lear, 2014). TDep products only employ a fusion approach to dry deposition and currently use an older version of CMAQ, although efforts to update the model version and incorporate wet-deposition fusion described herein are ongoing.

In this study, we assess deposition from the EPA's (Environmental Protection Agency) Air QUAlity TimE Series (EQUATES) project (, last access: 11 March 2022), which utilized the CMAQ version 5.3.2 modeling platform. The EQUATES project used consistent methodology (e.g., input data, processing methods, model versions) with deposition science updates compared to previous model versions, including the new Surface Tiled Aerosol and Gaseous Exchange (STAGE) bidirectional deposition model and improvements to organic N chemistry. Long-term air quality simulations like EQUATES are used by the scientific community in ecological and epidemiological assessments, including critical-load analyses. By employing these science updates in the model, consistent methodology, finer grid spacing (12 km), and a bias adjustment to wet deposition, we evaluate our best understanding of simulated deposition trends in the US. We evaluate the model agreement with observations of wet deposition and concentrations used to infer dry deposition and describe a measurement–model fusion technique in Sect. 2. In Sect. 3, we focus on the modeled trend and budget of N and S deposition over nine climatologically consistent regions within the contiguous US (CONUS) over a period of large emission reductions designated by the CAAA. Section 4 ends with Conclusions.

2 Methods and materials

2.1 EQUATES model configuration

The long-term EQUATES simulations from 2002–2017 were run using the CMAQ model version 5.3.2. Emissions were from different national emission inventories (NEIs) or used a base year NEI and estimated emissions for other years by creating scaling factors based on activity data and emission control information (Xing et al., 2013; Foley et al., 2020). A summary of recent CMAQ version 5.3 model updates and their impact on modeled concentrations is provided in Appel et al. (2021). The Weather Research and Forecasting version 4.1.1 (WRFv4.1.1) model was used for meteorology. Lateral boundary conditions for the 12 km grid spacing CONUS domain used in this study were provided by a 108 km grid spacing Northern Hemispheric simulation. Hemispheric and North American emission inventories were specifically prepared for EQUATES to ensure consistent input data and methods across all years. The Surface Tiled Aerosol and Gaseous Deposition (STAGE) option in CMAQv5.3.2 (Galmarini et al., 2021) was used to estimate atmospheric dry-deposition rates. The STAGE dry-deposition option generally performs similar to M3Dry (the other dry-deposition option in CMAQv5.3.2) when compared to network observations of ambient gaseous and aerosol pollutant concentrations (Appel et al., 2021) while providing additional land-use-specific deposition data useful for assessments of ecosystem exposure (Hood et al., 2021).

Figure 1Locations of the 200 National Atmospheric Deposition Program (NADP) National Trends Network (NTN, white circles) and 75 Clean Air Status and Trends Network (CASTNET, triangles) sites examined in this study. NADP NTN sites shown in black circles did not meet completeness criteria thresholds and therefore were not included in this analysis. Color-coded US climate regions shown in this map are referred to throughout this analysis. The black-bordered white circles indicate NADP NTN sites that meet annual completeness criteria for 13 years of the time series and are examined in the model evaluation presented in Sect. 3.1.

2.2 US deposition observations

We assess CMAQ's ability to reproduce annual accumulated wet deposition and ambient concentrations of N and S species measured by the NADP NTN and EPA's CASTNET networks (Fig. 1). Since anthropogenic sources dominate deposition and their area of influence is largely regional (Paulot et al., 2013), we evaluate the agreement over CONUS climate regions from 2002–2017 following the NOAA definition defined by Karl and Koss (1984), although these groupings do not imply the regions are isolated from each other. The nine climate regions include the Northwest, West, Northern Rockies, Southwest, South, Ohio Valley, Southeast, and Northeast as shown in Fig. 1.

The NADP's NTN (, last access: 11 March 2022) currently collects weekly precipitation samples at 263 locations across the US. Samples are analyzed by ion chromatography or flow injection analysis to quantify wet deposition of nitrate (NO3), ammonium (NH4), and sulfate (SO4), as well as several other compounds. First, modeled wet deposition of NO3, NH4, and SO4 is calculated using the approach by Appel et al. (2011) that accounts for chemical transformations of several species in the aqueous phase. Then, the model and observations are paired in time and space and annually accumulated. Following similar completeness criteria thresholds as in Zhang et al. (2019), NTN sites must be available for at least 13 of the 16 years simulated and have at least 60 % annual coverage each year to be considered for analysis. Applying these criteria, 200 NTN sites meet the minimum requirements for our analysis (Table S1 in the Supplement). The CASTNET measures concentrations and estimates dry deposition at 99 sites in the US (, last access: 11 March 2022). Weekly ambient concentrations of gases and particles are collected with an open-face three-stage filter pack. Since the reported dry-deposition values estimated by CASTNET are model derived, we instead evaluate the CMAQ performance simulating the ambient air concentrations of sulfur dioxide (SO2), sulfate (SO4), total nitrate (TNO3=NO3+HNO3), and ammonium (NH4). CASTNET sites available for at least 13 of the 16 years simulated with 75 % annual coverage each year are included in this analysis, resulting in 75 valid sites meeting the minimum requirements (Table S2). The EQUATES simulations are evaluated by comparison with available monitoring data using several model performance statistics (Simon et al., 2012).

2.3 Measurement–model fusion technique

The modeled wet-deposition fields are adjusted to account for input biases and uncertainty in the chemical and physical processes governing deposition. No corrections are applied to dry deposition due to limited dry-deposition measurements. Since model performance is improved for annual instead of seasonal values (Table S3; refer to Zhang et al., 2019, for detailed seasonal model evaluation), we apply a measurement–model fusion technique previously described by Zhang et al. (2019) to adjust the modeled annual wet-deposition fields of inorganic N (NO3+NH4) and S, briefly described here. First, the modeled wet-deposition fields from EQUATES are adjusted using observation-based 4 km precipitation fields generated by the Parameter-elevation Regressions on Independent Slopes Model (PRISM,, last access: 11 March 2022), annually accumulated and regridded to the 12 km CMAQ domain (EQUATESprecip-adj). While wet deposition relies on the season and precipitation rate (Schichtel et al., 2019), the relationship between precipitation and wet deposition is more linear when annually accumulated. Since PRISM incorporates observed elevation and climatic factors to estimate precipitation, NTN and TDep maps also used these data to improve estimates of wet deposition over complex terrain (Schwede and Lear, 2014). In each grid cell, the EQUATES modeled annual total wet deposition (WDmod) is adjusted by the ratio of annually accumulated precipitation from PRISM (Precipobs) to the WRF estimated annual accumulated rainfall (Precipmod) (Eq. 1):

(1) EQUATES precip-adj = Precip obs Precip mod × WD mod .

After adjusting simulated wet deposition by precipitation, an additional bias adjustment (EQUATESbias-adj) is applied using all NTN observations that meet annual data completeness, which varies year to year and includes sites of all classifications. For each NTN site, the median ratio between the EQUATESprecip-adj wet deposition and the NTN observations is calculated for all sites within a 300 km radius. Then, the median biases are interpolated using universal kriging with a linear trend in the spatial coordinates and an exponential covariance structure using Python's PyKrige module (, last access: 11 March 2022). A cross-validation analysis employing this measurement–model fusion technique with CMAQv5.0.2 is presented in Zhang et al. (2019). The final wet-deposition bias adjustment (Eq. 2) is computed by multiplying the precipitation-adjusted wet deposition by the inverse of the bias field, b:

(2) EQUATES bias-adj = 1 b × EQUATES precip-adj .

Figure 2Taylor plot comparing annual accumulated wet deposition (kg ha−1) of NH4 (a), NO3 (b), and SO4 (c) collected at NTN sites (black circles, along the x axis) with model output. The symbols differentiate between the various models (EQUATES and ECODEP) and wet-deposition corrections as described in the text. The azimuthal angle denotes the Pearson correlation coefficient (r2); the gold dashed radial distance shows the standard deviation (kg ha−1); and the dotted semicircles centered at the observation marker (black circle) denote the root-mean-square error.


Table 1EQUATES model performance metrics of annual (2002–2017) accumulated wet deposition of NH4, NO3, and SO4, including the Pearson correlation coefficient (r2), mean bias (MB, kg ha−1), and normalized mean bias (NMB, %).

Download Print Version | Download XLSX

Figure 3Scatterplots of annual accumulated bias-adjusted modeled and NTN observed wet deposition (kg ha−1) of ammonium (a, NH4), nitrate (b, NO3), and sulfate (c, SO4) from 2002 to 2017 colored by the climate region. Panel (d) shows NADP NTN observed and modeled precipitation (cm).


3 Results and discussion

3.1 Model performance

In this section, we investigate the EQUATES model performance, estimating wet deposition and ambient concentrations of N and S compounds throughout the US. The impact of the precipitation and bias correction technique on reproducing 2002–2017 accumulated measured wet deposition of NH4, NO3, and SO4 is shown in the Taylor diagram in Fig. 2 and summarized in Table 1. In addition, since previous modeling studies demonstrated moderate skill simulating wet deposition (Appel et al., 2011; Zhang et al., 2018), we also compare to an earlier time series utilizing CMAQv5.0.2 (ECODEP; 2002–2012 only). Little difference was found between EQUATES- and ECODEP-simulated wet-deposition statistics over the CONUS, although the improvement in the EQUATES estimates is more noticeable over regions with improved precipitation estimates. For example, a high precipitation bias in the WRFv3.4 meteorology used for ECODEP is reduced or eliminated in the WRFv4.1.1 meteorology used for EQUATES. The decrease in bias is likely attributed to the implementation of a hybrid coordinate system in WRFv4.1.1 which has been found to improve precipitation estimates at high-elevation sites in the western US (Beck et al., 2020), although the EQUATES precipitation is still biased low on average relative to PRISM, particularly in the Northwest (Fig. S1 in the Supplement). The close grouping of markers in Fig. 2a compared to the other species presented in Fig. 2 reflect the smaller standard deviation in the NH4 wet-deposition measurements (1.5 kg ha−1), which are comparable and slightly underestimated among the various models and adjustments. Applying the precipitation adjustment (EQUATESprecip-adj) reduces the underestimation of NH4 and SO4 wet deposition and increases the overestimation of NO3 wet deposition compared to the unadjusted model output (Table 1). The bias-adjusted NH4, NO3, and SO4 wet-deposition estimates have a similar standard deviation as the observations (falling near the dashed arc) and also feature a lower root-mean-square error (RMSE) (dotted semicircles) and higher correlation compared to the other estimates. Interannually, the bias correction to EQUATES output reduces the NMB of NH4 and SO4 wet deposition by 10 %–30 % and 5 %–20 %, respectively, compared to the unadjusted model output, with the largest improvement occurring from 2010–2017 (Fig. S2). Corrections applied to modeled NO3 wet deposition reduce the NMB by ∼20 % in the beginning of the time series (2002–2009) and <10 % at the end (2010–2017), possibly related to more uncertainty in emissions earlier in the earlier time period. The bias corrections to NH4 wet deposition are particularly noteworthy in the Upper Midwest, Northern Rockies, Ohio Valley, and South (Fig. S3), improving the agreement between the modeled and observed standard deviation, increasing the correlation, and lowering the RMSE compared to the unadjusted model values. The bias adjustment greatly improves the agreement of NO3 and SO4 wet deposition to observations in the Northeast, Northwest, and Ohio Valley by similarly improving the model performance metrics (Figs. S4 and S5). Due to the large improvements in estimating wet deposition using the measurement–model fusion approach, we use these corrections when assessing wet and total deposition throughout the rest of the analysis.

Table 2Evaluation of bias-adjusted SO4, NO3, and NH4 wet deposition and precipitation compared to NTN measurements in the contiguous US and nine climate regions that meet the completeness criteria (obs/mod: observed/modeled).

Download Print Version | Download XLSX

Figure 3 provides an overall comparison between the bias-corrected EQUATES wet deposition and NTN observations, along with an evaluation of simulated and observed precipitation. Table 2 provides summary statistics of wet deposition and precipitation for each climate region. A slight negative bias of modeled precipitation and bias-corrected NH4, NO3, and SO4 wet deposition is found for the entire CONUS, although agreement across the climate regions varies widely. The model generally exhibits better performance in the eastern climate regions, especially the Ohio Valley and Upper Midwest, compared to the rest of the CONUS. Agreement in the West, Northwest, and Northern Rockies is particularly poor, likely due to the errors simulating precipitation in the complex terrain (Table 2). In these western climate regions, precipitation is generally overestimated, consistent with previous model versions (Zhang et al., 2018).

The 16-year annual accumulated NH4 wet deposition is highest in the Upper Midwest, seen in both the model and measurements, and lowest in the Northwest. The model generally underestimates NH4 wet deposition for all regions except the Ohio Valley, with MB values ranging from 0.03 kg ha−1 in the Ohio Valley to −0.15 kg ha−1 in the Southeast. The correlation coefficients for wet NH4 deposition are on average stronger in the northern (r2>0.66) than southern US (r2<0.60). The weaker correlation between the measured and modeled trends of NH4 wet deposition (r2=0.43) is impacted by the smaller magnitude of the observed and modeled NH4 trends and reflects moderate model performance reproducing decreasing trends in many parts of the southern US and increasing trends in the northern US (Fig. S6). The 16-year annual accumulated NO3 wet deposition is also highest in the Upper Midwest but lowest in the West. The simulation overestimates NO3 wet deposition in the West, South, Ohio Valley, and Southwest with MB values all approximately 0.10 kg ha−1, while all other regions are underestimated with a MB ranging from −0.01 in the Northwest to −0.39 kg ha−1 in the Northern Rockies. Correlation coefficients for NO3 wet deposition are greater than 0.60 in all regions except the Southwest and Northern Rockies (r2∼0.50). Similar to NO3 wet deposition, the 16-year annual accumulated wet SO4 deposition is also highest in the Upper Midwest and lowest in the West. The simulation overestimates SO4 wet deposition across most of the southern US, where the MB ranges from 0.01 (Southwest) to 0.16 kg ha−1 (West). Most other regions show an underestimation of SO4 wet deposition with MB ranging from −0.03 (Ohio Valley) to −0.32 kg ha−1 (Northeast). Correlation coefficients are strongest for SO4 wet deposition than the other species, particularly in the northern US where r2 values are generally greater than 0.80. Like NH4 and NO3 wet deposition, correlation coefficients for SO4 wet deposition are weakest in the Southwest and West. The model also captures the clear downward trends of NO3 and SO4 wet deposition (r2=0.79 and 0.94, respectively), although the magnitudes of the decreasing modeled trends are slightly underestimated (Fig. S6).

Figure 4Scatterplots comparing 2002–2017 annual average concentrations (µg m−3) of sulfur dioxide (a, SO2), sulfate (b, SO4), total nitrate (c, TNO3), and ammonium (d, NH4) between EQUATES output and CASTNET observations (N=1175). The marker color corresponds to the climate region. Data from CASTNET sites with at least 75 % valid data available for each of the 13 years analyzed are included in the analysis shown.


To indirectly evaluate the ability of the EQUATES simulations to reproduce dry deposition, we compare the model-simulated annual average concentrations of SO2, SO4, TNO3, and NH4 with observations collected at CASTNET sites in Fig. 4 and in Table 3. Correlations between measured and modeled mean concentrations show good agreement for NH4 (r2=0.94), TNO3 (r2=0.81), SO2 (r2=0.90), and SO4 (r2=0.96). Overall, the model shows negatives biases for NH4 and SO4 concentrations (−17.4 % and −15.4 %, respectively) and positive biases for TNO3 and SO2 concentrations (2.8 % and 26.0 %, respectively) but generally improved from CMAQv5.0 (Xing et al., 2015). The highest 16-year average concentrations for all species are found in the Upper Midwest for both the simulations and observations. The model underestimates concentrations of NH4 across all regions, with MB ranging from −0.02 (Northern Rockies) to −0.22µg m−3 (Ohio Valley). Correlation coefficients for NH4 concentrations are generally strong (r2>0.80) everywhere, except the West and Southwest (r2=0.09 and 0.39, respectively). Average modeled TNO3 concentrations are mostly overestimated in the northern US, with MB ranging from 0.02 (Northern Rockies) to 0.26 µg m−3 (Ohio Valley), and underestimated in the southern US, with MB ranging from −0.01 (Southwest) to −0.64µg m−3 (West). Concentrations of SO2 are generally overestimated except in the western regions (West, Southwest, and Northern Rockies), and strong correlations are found throughout all regions. The model underestimates SO4 concentrations in most regions, with MB ranging from −0.08 (West) to −0.58µg m−3 (Southeast), while it overestimates it in the Northern Rockies (MB =0.07µg m−3) and Southwest (MB =0.01µg m−3). Correlation coefficients for SO4 are also generally strong (r2>0.80), with exceptions in the Southwest and West (r2=0.59 and 0.31, respectively). The EQUATES simulations reproduce the decreasing trend in observed concentrations with strong correlations (r2>0.79) (Fig. S7) but overpredict trends of SO2 and TNO3 concentrations and underpredicts trends of NH4 and SO4 concentrations.

Table 3Evaluation of EQUATES modeled average concentrations (µg m−3) of SO2, SO4, TNO3, and NH4 compared to CASTNET measurements in the CONUS and eight climate regions (no CASTNET data are available in the Northwest climate region; obs/mod: observed/modeled).

Download Print Version | Download XLSX

3.2 Trend analysis

3.2.1 Total deposition trends

Figure 5 shows the spatial distribution and overall trend of modeled total deposition of N (i.e., the sum of oxidized and reduced N) and S between 2002 and 2017. For the trend analysis shown here and throughout, we calculate a linear least-squares regression with significance examined at the 95 % confidence level using a Wald test, with any insignificant trends set to zero. The eastern US has higher total N deposition amounts than the western US in both 2002 and 2017, particularly in parts of the Ohio Valley and Upper Midwest, regions associated with large N emissions from agriculture (Dammers et al., 2019) and energy consumption (Kim et al., 2006). The major N deposition region shifts from the eastern to the central US between 2002 and 2017, with highest deposition amounts found in parts of Iowa, North Carolina, and Indiana, a pattern consistent with N inventories compiled by Sabo et al. (2019) indicating an increase in agricultural fertilizer and livestock waste in the Midwest between 2002 and 2012. Increases in agricultural fertilizer have previously been linked to cropland expansion in the Great Plains (Lark et al., 2015; Wright et al., 2017) and increasing demand for domestically sourced biofuels (Donner and Kucharik, 2008). Urban regions in the central and eastern US indicate a substantial amount of N deposition compared to nearby rural areas (Fig. S8), consistent with previous findings that bulk N deposition in urban areas is twice as much as rural and remote sites (Decina et al., 2019).

Similar to total N deposition, total S deposition has a noticeable spatial gradient in the east compared to the west, particularly in the beginning of the time series, following larger changes in emissions in these regions (Aas et al., 2019; Holland et al., 1999). In 2002, total S deposition across the Southeast, Ohio Valley, and Northeast is broadly greater than 10 kg S ha−1, whereas S deposition for most of the western and central US climate regions is less than 5 kg S ha−1. By 2017, total S deposition continues to be highest in the eastern US, although total S deposition has decreased to less than 7.5 kg ha−1 on average. As the primary ion in acid deposition, S has been shown to reduce the biodiversity of natural ecosystems (Clark et al., 2019) but is also an essential macro-nutrient required for crop yield and health (Hinckley et al., 2020). Model projections spanning a range of future climate emission scenarios at the end of the century project global S deposition to agricultural soils to decrease by 70 %–90 % compared to 2005–2009 (Feinberg et al., 2021). Consequences of declining S deposition under growing food demand include risk of S deficiencies unless mitigating fertilizer strategies are developed and implemented. The spatial changes in decreasing total N and S deposition track changes in SO2 and NOx emissions since 1970 (Nopmongcol et al., 2019), demonstrating the success of air quality mitigation policies implemented to achieve the National Ambient Air Quality Standards (NAAQS) under the CAAA (United States Congress, 1990) in also reducing nutrient and acid deposition. Concentrations of SO2, a primarily directly emitted compound, indicate greater decreasing trends than the secondary formed SO4. The non-linear relationship between SO2 and SO4 trends is shown in both the measurements and model output, as well as noted in other studies using several global CTMs (Aas et al., 2019) and monitoring data (Sickles and Shadwick, 2015), which attributed the increasing oxidative capacity of the atmosphere more efficiently converting SO2 to SO4 (Sickles and Shadwick, 2015). As SO2 emissions decline from coal-burning power plants due to controls implemented to meet the NAAQS (United States Congress, 1990), additional oxidants are available to oxidize SO2. The decline in SO2 emissions accompanied by near-constant NH3 emissions results in less acidic cloud droplets (Redington et al., 2009; Pye et al., 2020), which increases the oxidation rate of SO2 via the ozone pathway (Paulot et al., 2017). Lastly, the trends in wet deposition of SO4 (Fig. S6) are larger than the trends in concentration (Fig. S7), as SO4 is more efficiently scavenged by precipitation.

Figure 5Spatial distribution of total N (a, b, c) and S (d, e, f) deposition in 2002 (a and d, kg ha−1) and 2017 (b and e, kg ha−1) and the 2002–2017 annual trend (c and f, kg ha−1 yr−1) with significance at the 95 % confidence level. Grey areas in panels (c) and (f) indicate where the trend is unavailable due to a lack of PRISM data or not significant (i.e., p value of the Wald test is greater than 0.05).

Figure 6 summarizes average total N and S trends across the climate regions and CONUS based on three different time bins: 2002–2017, 2002–2009, and 2010–2017. While there are significant regional deposition composition changes from 2002–2017, smaller changes in the overall deposition trends are observed for the CONUS. From 2002 to 2017, the largest average trend in decreasing total N deposition (−0.27 to −0.31 kg N ha−1 yr−1) occurs in the Ohio Valley, Northeast, and Southeast. Similarly, the Ohio Valley, Northeast, and Southeast experience the fastest rate of total S deposition decline from 2002 to 2017 (−0.68 to −0.86 kg S ha−1 yr−1), while the West, Northwest, Northern Rockies, and Southwest indicate near-zero (−0.03 kg S ha−1 yr−1) rates of decline. In many regions, especially the Northeast, Ohio Valley, and Southeast, a larger decreasing trend in total N and S deposition is found between 2002–2009 than 2010–2017. The high amount of total N and S deposition found in these regions combined with the slower rate of decline in recent years suggests additional emission controls may be necessary to continue deposition declines. On the other hand, the South is the only region with a very small increasing trend of ∼0.01 kg N ha−1 yr−1 of total N deposition with large variability from 2010–2017 but also indicates a larger declining trend of −0.08 kg N ha−1 yr−1 from 2002–2009. Similarly, a large decreasing trend in total N deposition is found in the Midwest from 2002–2009, while the trend is near zero from 2010–2017. Fertilizer use in the Midwest has only grown modestly (∼1.3 % yr−1), so increases in total N deposition have been largely attributed to increasing reduced N deposition. Growing-reduced N deposition is a result of NH3 emissions increasing exponentially with temperature (except below freezing and where emissions are near 0) (Riddick et al., 2016) and by increasing the partitioning of NH3 remaining in the gas phase due to NOx and SO2 emission reductions (Warner et al., 2017) and increasing importance of on-road mobile emissions of NH3 (Fenn et al., 2018). As the world's demand for food grows and NH3 emissions from agriculture increase, models project increasing Nr deposition in many areas across the globe by 2050 (Paulot et al., 2013) and 2100 (Lamarque et al., 2011).

Figure 6Comparison of mean total N (a) and total S (b) deposition trends in the nine CONUS climate regions from 2002–2017, 2002–2009, and 2010–2017. Lines on the bars extend to the 5th and 95th percentiles of the trends. Trends include all grid cell locations (insignificant trends are set to 0).


To elucidate the chemical drivers responsible for the varying changes in total N deposition across the US, we examine trends in total oxidized and reduced N in Fig. 7. Overall, most regions show larger decreasing trends in oxidized N than increasing trends in reduced N. Regions with insignificant total trends (see grey areas in Fig. 5) indicate similar magnitudes of oxidized and reduced N trends. The largest decreasing trends in total oxidized N are found in the eastern US and along the Pacific coast. The large declines in total oxidized N in the Upper Midwest, Ohio Valley, Northeast, and Southeast showcase both the coordinated application of emission controls in major contributing regions and increasingly lower amounts of NOx transported from source to receptor areas (Lloret and Valiela, 2016). Unlike trends of total N deposition (Fig. 5), statistically significant declining trends for total oxidized N are found throughout most of the CONUS, including the western US. No east–west spatial gradient is observed for trends in total reduced N deposition. Instead, reduced N deposition increases throughout most of the US, particularly in the Upper Midwest and Ohio Valley at an average rate of 0.09 kg N ha−1 yr−1, but is masked by the larger decreasing oxidized N trends. Declines in total reduced N deposition are located in Los Angeles, California; Tulsa, Oklahoma; and Baton Rouge, Louisiana, suggesting decreases in NH3 emissions near these cities. High NH3 deposition in Baton Rouge (>10 kg ha−1) has been attributed to ammonia plant operations in the area (Guo et al., 2018), while dairy farms (Nowak et al., 2012) and vehicles (Cao et al., 2021) are major sources of NH3 in Los Angeles.

Figure 7Maps of 2002–2017 trends of total oxidized (a) and reduced N deposition (b, kg N ha−1 yr−1) and comparison of average trends across nine CONUS climate regions (c). Lines in panel (c) show the 5th and 95th percentiles of the trends for each region. Grey colors on the map indicate where the trend is unavailable or insignificant.

3.2.2 Wet- and dry-deposition trends

Figure 8 shows trends of the wet and dry components of oxidized and reduced N and S deposition. Average dry deposition of oxidized N is found to decline faster than average wet deposition of oxidized N for all climate regions, as oxidized N is more efficiently removed by dry deposition due to high deposition velocities (Zhang et al., 2012). Larger amounts of dry vs. wet deposition across the eastern US in 2002 (Fig. S9a and d) and more comparable amounts in 2017 (Fig. S9b and e) suggest a shift over time in the contributions of oxidized wet and dry deposition, although the portioning is highly dependent upon annual rainfall and changes to climate. The Southeast, Northeast, and Ohio Valley have the steepest simulated average dry and wet oxidized N deposition trends (approximately −0.3 and −0.1 kg N ha−1 yr−1, respectively), consistent with the observed trends (Fig. S10). The NTN sites with the highest mean NO3 deposition, particularly those in the Northeast, experience the largest decreasing trends between 2002–2009 (−2.2 kg N ha−1 yr−1) compared to 2010–2017 (−1.0 kg N ha−1 yr−1), consistent with the trends shown in Fig. 6.

Trends of reduced N are increasing and overall smaller in magnitude compared to the oxidized N trends (Fig. 7b). Regions of elevated wet and dry reduced N deposition have expanded and increased in magnitude across the CONUS (Fig. S11) compared to oxidized N, also observed in the NTN NH4 measurements. The increasing trends of NH4 wet deposition accompany statistically significant increasing trends in annual precipitation predominantly over the Upper Midwest, Southwest, and Northern Rockies (Fig. S12). Similar increasing trends of NH4 wet deposition between 1989 and 2016 observed in the Midwest and Mid-Atlantic have also been linked to rising annual precipitation amounts (Feng et al., 2021). Similarly, NTN sites in the Upper Midwest with the highest mean NH4 wet deposition experience larger increasing trends from 2002 to 2017 (Fig. S10). Similar statistically significant increases in NH4 wet-precipitation concentrations have been reported, although increases are largest for the period prior to our model simulations (1985–1999) (Mchale et al., 2021). The decreasing reduced N deposition in the three urban centers discussed in Sect. 3.2.1 is present in the dry-deposition trends but not in wet deposition. The Southeast and Northeast have a near-zero trend in wet reduced N deposition but a larger increasing trend in dry reduced N deposition. Parts of the Southeast, including southwestern Florida, experience a decrease in precipitation across the entire period, which may partially explain the decreasing trend of wet reduced N deposition.

Decreasing trends in wet and dry S deposition are found across the CONUS, although the amount varies by climate region (Figs. 8c, S13). Similar to the oxidized N trends, the Northeast, Southeast, and Ohio Valley experience the largest decreasing trend of approximately −0.5 kg S ha−1 yr−1 for dry deposition and −0.3 kg S ha−1 yr−1 for wet deposition, consistent with the NTN SO4 observations (Fig. S10). Sites with the largest mean SO4 wet deposition in the Ohio Valley and Northeast also experience the largest decreasing trends in wet SO4 deposition (−1.9 kg S ha−1 yr−1, Fig. S10), consistent with monthly mean precipitation-weighted concentrations from NTN sites from 1985 to 2017 (Mchale et al., 2021). Large decreasing trends in dry S deposition are clearly seen over large urban centers in the western US, while the eastern US experiences a significant declining trend everywhere in response to decreases in regional SO2 emissions (Likens et al., 2001).

Figure 8Annual average trend (2002–2017, kg ha−1 yr−1) of wet (blue) and dry (green) oxidized nitrogen (a), reduced nitrogen (b), and sulfur (c) deposition throughout the CONUS and nine climate regions. Whiskers extend to the 5th and 95th percentiles of the trends.


Figure 9Average total (wet + dry) N (a) and S (b) deposition in 2002 (left bars) and 2017 (right bars) throughout the nine climate regions and CONUS. Colors denote the wet and dry component of each species. The numbers in each bar denote the percentage contribution to the annual total.


3.3 Total nitrogen and sulfur deposition budgets

The average N and S deposition budgets across all CONUS climate regions generally decrease from 2002 to 2017 (Fig. 9). The average total N budget for the CONUS decreases from 7.8 to 6.3 kg N ha−1 yr−1 between 2002 and 2017. The Upper Midwest, Ohio Valley, Northeast, and Southeast have the largest average total N, decreasing from 12–14 in 2002 to 8–10 kg N ha−1 yr−1 in 2017. In contrast, the Northwest, West, and Southwest have the smallest average total N, which decreases only slightly from 3.1–3.9 to 2.8–3.5 kg N ha−1 yr−1 over the same time period. The decrease in total N deposition is largely due to reductions in oxidized N, particularly from dry deposition, and declines in NOx emissions (Nopmongcol et al., 2019). The average wet oxidized N deposition for the CONUS decreases from 1.6 in 2002 to 1.1 kg N ha−1 yr−1 in 2017, while the dry component decreases from 3.7 to 1.7 kg N ha−1 yr−1. Dry deposition constitutes most of the total N budget (slightly decreasing from 50 %–82 % in 2002 to 42 %–76 % in 2017 across the climate regions), with the largest average contributions found in the western US. Spatial differences in the dry N deposition flux, inferred using Ozone Monitoring Instrument (OMI) satellite observations of NO2 and ground measurements from 2005 to 2014, indicate an increasing flux across the western US and decreasing flux over the eastern US due to reductions in NOx and NH3 emissions (Jia et al., 2016). Unlike oxidized N, the wet reduced N deposition is larger on average than the dry deposition for all regions except the West and Northwest, regions that typically experience less rainfall than the other regions. Similar to Du et al. (2014), we find an increasing contribution of wet reduced N deposition to the total across the CONUS from 19 % to 30 % across the time period (Fig. S14).

Figure 10Changes in the percent of total N deposited as reduced nitrogen from 2002 to 2017 throughout the NOAA climate regions and CONUS.


The relative proportion of reduced N to the total budget has shifted over the time period of 2002 to 2017 (Fig. 10). On average across the CONUS, oxidized N deposition dominates the total budget until 2016, when reduced N comprises >50 % of the total, consistent with previous studies (Du et al., 2014; Kanakidou et al., 2016; Ackerman et al., 2019; Li et al., 2016; Nopmongcol et al., 2019) and attributable to the absence of NH3 emission controls from agricultural sources combined with the success of NOx regulations. The reductions in NOx and SOx emissions due to policy controls result in less NHx partitioning to the aerosol phase and more to the gas phase, favoring dry deposition closer to source regions. Despite spatial variations in agricultural NH3 sources and a lack of an emission trend during this study period over some locations, all climate regions indicate an increasing amount of reduced N to the total budget, agreeing with results from Nopmongcol et al. (2019). In order to protect critical ecosystems within US national parks from future N deposition, Ellis et al. (2013) project at least a 50 % reduction in anthropogenic NH3 emissions is required relative to RCP-projected (Representative Concentration Pathway) emission scenario levels in 2050. The fraction of reduced N to the total is largest in regions with larger NH3 emissions (Dammers et al., 2019; Behera et al., 2013), such as the Upper Midwest and Northern Rockies, throughout the time period studied here. These regions exceed 50 % of reduced N as a fraction of the total budget as early as 2006 for the Upper Midwest and 2011 for the Northern Rockies. Approximately 80 % of the total global NH3 emissions in 2005 were attributable to agriculture, including animal feedlot operations, and are anticipated to increase in the US and globally (Behera et al., 2013). Present-day (2015) global simulations indicate substantial regional variation in the Nr budget, reflecting the impact of regional activities and meteorology (Ge et al., 2022). The increasing importance of reduced N is likely to impact future N deposition budgets and the competitive nature among plants with varying affinities for the different forms of nitrogen (Choudhary et al., 2016; Kahmen et al., 2006).

The total S budget also decreases across the CONUS on average, falling from 5.3 in 2002 to 1.8 kg S ha−1 yr−1 in 2017. Decreasing SO2 emissions have been previously shown to be responsible for decreasing S deposition (Nopmongcol et al., 2019). Similar to the total N budget, the largest average amounts of total S are located in the Ohio Valley, Northeast, and Southeast, decreasing from 11–14 in 2002 to 8–10 kg S ha−1 yr−1 in 2017. In the West, Ohio Valley, Northeast, and Southeast, dry deposition dominates the budget for 2002, while transitioning wet deposition dominated by 2017. Wet deposition comprises a greater portion of the total S budget in 2017 across all climate regions. The primary driver behind the shift in the total S budget may be emission reductions in SO2 implemented to meet the NAAQS (leading to decreasing ambient concentrations of SO2 and SO4) that more efficiently decrease the dry component than the wet (Sickles and Shadwick, 2015).

4 Conclusions

In this study, we examine CMAQ model simulations from the EPA's Air QUAlity TimE Series (EQUATES) project to investigate spatial and temporal trends of nitrogen (N) and sulfur (S) deposition over a period with substantial emission reductions implemented to meet requirements set by the Clean Air Act Amendments (CAAA). We assess changes in modeled dry deposition and precipitation- and bias-adjusted wet deposition across nine climatologically consistent regions within the US. EQUATES CMAQ simulations included important science updates to previous modeling studies, such as a revised model for bidirectional air–surface exchange of NH3 and a mechanistic representation of organic N (Pye et al., 2015), necessary to improve modeling of Nr budgets. The measurement–model fusion technique for modeled wet-deposition estimates reduces the model bias and improves correlations compared to the unadjusted model values and can thus improve the accuracy of critical-load assessments based on model data. Comparisons to wet-deposition observations from the NADP network indicate the model generally underestimates wet deposition of NH4, NO3, and SO4 across the CONUS. Agreement is poor in regions with complex terrain like the Northern Rockies and West, potentially due to errors constraining emissions and/or reproducing precipitation at 12 km grid spacing, as sub-grid variability in precipitation estimates can differ from sub-grid variability in concentrations across the country. Similarly, the model underestimates ambient concentrations of NH4 and SO4 compared to CASTNET measurements across the CONUS, while SO2 and TNO3 concentrations are overestimated. Additional measurements in regions with high bias, such as the West, Northwest, and Southwest, as well as model improvements to precipitation estimates, particularly in areas of complex terrain, could greatly improve agreement with the NADP network observations.

The modeled total N and S deposition amounts and trends are larger across the eastern than the western US, particularly in the Northeast, Southeast, Ohio Valley, and Upper Midwest climate regions. While the average CONUS trend indicates a decrease of −0.11 kg N ha−1 yr−1 in total N deposition from 2002–2017, significant increasing trends are found in localized parts of the Upper Midwest, Northern Rockies, and West. Such increases in total N deposition are dominated by increases in reduced N that are detected across all climate regions but are masked by the larger declining trends of oxidized N deposition as a result of policies targeting NOx emissions enacted under the CAAA. The increasing trend in reduced N from 2002 to 2017 is a result of growing agricultural emissions combined with NOx/SOx emission reductions and higher precipitation amounts over time concentrated in the Upper Midwest, Southwest, and Northern Rockies. Despite the overall average declining N trend from 2002 to 2017, the simulations indicate a larger average decreasing trend of total N and S at the beginning of the time series (2002 to 2009) than the end of the period (2010 to 2017). Total N trends in the South reveal a slight increasing trend (0.01 kg N ha−1 yr−1) from 2010 to 2017, contrary to the larger declining trends found between 2002–2009. The modeled slowing and slightly increasing trend for total N across the climate regions in the EQUATES simulations hints at future changes to the N budget that may impact plant biodiversity and productivity (Kahmen et al., 2006).

The average total N deposition budget over the CONUS has decreased from 7.8 in 2002 to 6.3 kg N ha−1 yr−1 in 2017, largely due to decreases in oxidized N deposition as a result of NOx emission controls implemented under the CAAA. The major contributor to the N budget is dry deposition, although the contribution across the climate regions has declined from 50 %–82 % in 2002 to 42 %–76 % in 2017, in part due to changes in precipitation amounts. The regions with the largest average total N deposition, including the Ohio Valley, Northeast, Upper Midwest, and South, also experience the largest reduction in total N deposition. Similar to the total N budget, the average total S budget over the CONUS has declined from 5.3 in 2002 to 1.8 kg S ha−1 yr−1 in 2017. The largest amount of total S occurs in the eastern US, where total deposition is more impacted by dry deposition. By the end of the time series, wet deposition dominates the S budget for all climate regions, comprising 68 % of the total budget in the CONUS on average.

Our analysis, in addition to the analyses of many others (Zhang et al., 2018; Li et al., 2016; Du et al., 2014), highlights the increasing contribution of reduced compounds to the N budget in all climate regions, particularly the Northern Rockies and Upper Midwest. Both the model and observations indicate statistically significant increasing trends in NH4 deposition across the CONUS due to changes in precursor emissions. Reductions in NOx and SO2 emissions, and therefore their oxidation products, due to policy controls result in less NHx partitioning to the aerosol phase and more to the gas phase, favoring dry deposition of NH3 closer to source regions. Additionally, growth of NH3 emissions from agricultural sources (Behera et al., 2013) allows for additional NH4 aerosol formation in the atmosphere given sufficient NOx and SOx. This study further confirms the growth of reduced N in recent years across all climate regions and additionally suggests an attenuation of the declining oxidized N deposition trend. However, since N deposition over urban areas across the CONUS is likely already underestimated (Rao et al., 2013), with increasing urbanization only expected to further increase N deposition amounts (Joyce et al., 2020), addressing modeling uncertainty in emissions and chemistry at relevant spatial and temporal resolutions is imperative. For instance, a more robust characterization of dry deposition in the US, including low-cost dry deposition and ambient measurement methods, is needed, particularly in agricultural areas and regions with transitions from urban to rural environments, as well as for species not routinely measured, such as organic compounds. In this regard, additional ambient NH3 measurements are of critical importance for improving the modeling and analytical methods used in air quality (and related ecosystem) management for particulate matter and N deposition.


The supplement related to this article is available online at:

Author contributions

KMF and GP were the EQUATES project leads and responsible for project administration, methodology, validation, formal analysis, data curation, and visualization. RG conducted the meteorological runs, and KWA and CH performed the model simulations. RG, CH, KWA, and JOB provided the methodology, formal analysis, and validation. SEB performed the data analysis and created all figures and tables. SEB wrote the paper with comments from all authors.

Competing interests

The contact author has declared that none of the authors has any competing interests.


The views expressed in this paper are those of the authors and do not necessarily represent the views or policies of the US Environmental Protection Agency.

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The authors would like to thank Chris Allen, Michael Aldridge, Megan Beardsley, James Beidler, David Choi, Alison Eyth, Caroline Farkas, Janice Godfrey, Barron Henderson, Shannon Koplitz, Rich Mason, Rohit Mathur, Chris Misenis, Norm Possiel, Havala Pye, Lara Reynolds, Matthew Roark, Sarah Roberts, Donna Schwede, Karl Seltzer, Darrell Sonntag, Kevin Talgo, Claudia Toro, and Jeff Vukovich for their contribution to the development of emission and meteorology inputs for the EQUATES simulations used in this study. This research was supported in part by an appointment to the Research Participation Program at the US EPA Office of Research and Development (ORD), administered by the Oak Ridge Institute for Science and Education (ORISE) through an interagency agreement between the US Department of Energy and the US EPA.

Review statement

This paper was edited by Qiang Zhang and reviewed by Ruth Heindel and Yuqiang Zhang.


Aas, W., Mortier, A., Bowersox, V., Cherian, R., Faluvegi, G., Fagerli, H., Hand, J. L., Klimont, Z., Galy-Lacaux, C., Lehmann, C. M. B., Myhre, C. L., Myhre, G., Olivie, D., Sato, K., Quaas, J., Rao, P. S. P., Schulz, M., Shindell, D. T., Skeie, R. B., Stein, A., Takemura, T., Tsyro, S., Vet, R., and Xu, X.: Global and regional trends of atmospheric sulfur, Sci. Rep.-UK, 9, 953,, 2019. 

Ackerman, D., Millet, D. B., and Chen, X.: Global Estimates of Inorganic Nitrogen Deposition Across Four Decades, Global Biogeochem. Cy., 33, 100–107,, 2019. 

Appel, K. W., Foley, K. M., Bash, J. O., Pinder, R. W., Dennis, R. L., Allen, D. J., and Pickering, K.: A multi-resolution assessment of the Community Multiscale Air Quality (CMAQ) model v4.7 wet deposition estimates for 2002–2006, Geosci. Model Dev., 4, 357–371,, 2011. 

Appel, K. W., Bash, J. O., Fahey, K. M., Foley, K. M., Gilliam, R. C., Hogrefe, C., Hutzell, W. T., Kang, D., Mathur, R., Murphy, B. N., Napelenok, S. L., Nolte, C. G., Pleim, J. E., Pouliot, G. A., Pye, H. O. T., Ran, L., Roselle, S. J., Sarwar, G., Schwede, D. B., Sidi, F. I., Spero, T. L., and Wong, D. C.: The Community Multiscale Air Quality (CMAQ) model versions 5.3 and 5.3.1: system updates and evaluation, Geosci. Model Dev., 14, 2867–2897,, 2021. 

Beck, J., Brown, J., Dudhia, J., Gill, D., Hertneky, T., Klemp, J., Wang, W., Williams, C., Hu, M., James, E., Kenyon, J., Smirnova, T., and Kim, J.-H.: An Evaluation of a Hybrid, Terrain-Following Vertical Coordinate in the WRF-Based RAP and HRRR Models, Weather Forecast., 35, 1081–1096,, 2020. 

Behera, S. N., Sharma, M., Aneja, V. P., and Balasubramanian, R.: Ammonia in the atmosphere: a review on emission sources, atmospheric chemistry and deposition on terrestrial bodies, Environ. Sci. Pollut. R., 20, 8092–8131,, 2013. 

Bettez, N. D., Marino, R., Howarth, R. W., and Davidson, E. A.: Roads as nitrogen deposition hot spots, Biogeochemistry, 114, 149–163,, 2013. 

Bobbink, R., Hicks, K., Galloway, J. N., Spranger, T., Alkemade, R., Ashmore, M., Bustamante, M., Cinderby, S., Davidson, E., Dentener, F., Emmett, B., Erisman, J. W., Fenn, M. E., Gilliam, F. S., Nordin, A., Pardo, L. H., and De Vries, W.: Global assessment of nitrogen deposition effects on terrestrial plant diversity: a synthesis, Ecol. Appl., 20, 30–59,, 2010. 

Cao, H., Henze, D. K., Cady-Pereira, K., McDonald, B. C., Harkins, C., Sun, K., Bowman, K. W., Fu, T.-M., and Nawaz, M. O.: COVID-19 Lockdowns Afford the First Satellite-Based Confirmation That Vehicles Are an Under-recognized Source of Urban NH3 Pollution in Los Angeles, Environ. Sci. Tech. Let., 9, 3–9, 10.1021/acs.estlett.1c00730, 2021. 

Choudhary, S., Blaud, A., Osborn, A. M., Press, M. C., and Phoenix, G. K.: Nitrogen accumulation and partitioning in a High Arctic tundra ecosystem from extreme atmospheric N deposition events, Sci. Total Environ., 554–555, 303–310,, 2016. 

Clark, C. M., Simkin, S. M., Allen, E. B., Bowman, W. D., Belnap, J., Brooks, M. L., Collins, S. L., Geiser, L. H., Gilliam, F. S., Jovan, S. E., Pardo, L. H., Schulz, B. K., Stevens, C. J., Suding, K. N., Throop, H. L., and Waller, D. M.: Potential vulnerability of 348 herbaceous species to atmospheric deposition of nitrogen and sulfur in the United States, Nat. Plants, 5, 697–705,, 2019. 

Dammers, E., McLinden, C. A., Griffin, D., Shephard, M. W., Van Der Graaf, S., Lutsch, E., Schaap, M., Gainairu-Matz, Y., Fioletov, V., Van Damme, M., Whitburn, S., Clarisse, L., Cady-Pereira, K., Clerbaux, C., Coheur, P. F., and Erisman, J. W.: NH3 emissions from large point sources derived from CrIS and IASI satellite observations, Atmos. Chem. Phys., 19, 12261–12293,, 2019. 

Decina, S. M., Hutyra, L. R., and Templer, P. H.: Hotspots of nitrogen deposition in the world's urban areas: a global data synthesis, Front. Ecol. Environ., 18, 92–100,, 2019. 

Donner, S. D. and Kucharik, C. J.: Corn-based ethanol production compromises goal of reducing nitrogen export by the Mississippi River, P. Natl. Acad. Sci. USA, 105, 4513–4518,, 2008. 

Du, E., De Vries, W., Galloway, J. N., Hu, X., and Fang, J.: Changes in wet nitrogen deposition in the United States between 1985 and 2012, Environ. Res. Lett., 9, 095004,, 2014. 

Ellis, R. A., Jacob, D. J., Sulprizio, M. P., Zhang, L., Holmes, C. D., Schichtel, B. A., Blett, T., Porter, E., Pardo, L. H., and Lynch, J. A.: Present and future nitrogen deposition to national parks in the United States: critical load exceedances, Atmos. Chem. Phys., 13, 9083–9095,, 2013. 

EPA (United States Environmental Protection Agency): EQUATES: EPA's Air QUAlity TimE Series Project,, last access: 11 September 2022. 

Feinberg, A., Stenke, A., Peter, T., Hinckley, E.-L. S., Driscoll, C. T., and Winkel, L. H. E.: Reductions in the deposition of sulfur and selenium to agricultural soils pose risk of future nutrient deficiencies, Commun. Earth Environ., 2, 101,, 2021. 

Feng, J., Vet, R., Cole, A., Zhang, L., Cheng, I., O'Brien, J., and Macdonald, A.-M.: Inorganic chemical components in precipitation in the eastern U.S. and Eastern Canada during 1989–2016: Temporal and regional trends of wet concentration and wet deposition from the NADP and CAPMoN measurements, Atmos. Environ., 254, 118367,, 2021. 

Fenn, M. E., Bytnerowicz, A., Schilling, S. L., Vallano, D. M., Zavaleta, E. S., Weiss, S. B., Morozumi, C., Geiser, L. H., and Hanks, K.: On-road emissions of ammonia: An underappreciated source of atmospheric nitrogen deposition, Sci. Total Environ., 625, 909–919,, 2018. 

Foley, K. M., Possiel, N., Aldridge, M., Allen, M., Bash, J. O., Beardsley, M., Beidler, J., Choi, D., Eder, B., Farkas, C., Gilliam, R., Godfrey, J., Henderson, B. H., Hogrefe, C., Koplitz, S. N., Mason, R., and Vukovich, J.: EQUATES: EPA's Air QUAlity TimE Series Project, 19th Annual CMAS Conference, Virtual, 28 October 2020, (last access: 11 March 2022), 2020. 

Fowler, D., Coyle, M., Skiba, U., Sutton, M. A., Cape, J. N., Reis, S., Sheppard, L. J., Jenkins, A., Grizzetti, B., Galloway, J. N., Vitousek, P., Leach, A., Bouwman, A. F., Butterbach-Bahl, K., Dentener, F., Stevenson, D., Amann, M., and Voss, M.: The global nitrogen cycle in the twenty-first century, Philos. T. Roy. Soc. B, 368, 20130164,, 2013. 

Galloway, J. N., Townsend, A. R., Esirman, J. W., Bekunda, M., Cai, Z., Freney, J. R., Martinelli, L. A., Seitzinger, S. P., and Sutton, M. A.: Transformation of the Nitrogen Cycle: Recent Trends, Questions, and Potential Solutions, Science, 320, 889–892,, 2008. 

Galmarini, S., Makar, P., Clifton, O. E., Hogrefe, C., Bash, J. O., Bellasio, R., Bianconi, R., Bieser, J., Butler, T., Ducker, J., Flemming, J., Hodzic, A., Holmes, C. D., Kioutsioukis, I., Kranenburg, R., Lupascu, A., Perez-Camanyo, J. L., Pleim, J., Ryu, Y.-H., San Jose, R., Schwede, D., Silva, S., and Wolke, R.: Technical note: AQMEII4 Activity 1: evaluation of wet and dry deposition schemes as an integral part of regional-scale air quality models, Atmos. Chem. Phys., 21, 15663–15697,, 2021. 

Ge, Y., Vieno, M., Stevenson, D. S., Wind, P., and Heal, M. R.: A new assessment of global and regional budgets, fluxes, and lifetimes of atmospheric reactive N and S gases and aerosols, Atmos. Chem. Phys., 22, 8343–8368,, 2022. 

Guo, H., Han, F., Wang, Z., Pardue, J., and Zhang, H. L.: Deposition of sulfur and nitrogen components in Louisiana in August, 2011, Sci. Total Environ., 636, 124–133,, 2018. 

Hinckley, E. L. S., Crawford, J. T., Fakhraei, H., and Driscoll, C. T.: A shift in sulfur-cycle manipulation from atmospheric emissions to agricultural additions, Nat. Geosci., 13, 597–604,, 2020. 

Holland, D. M., Principe, P. P., and Sickles, J. E.: Trends in atmospheric sulfur and nitrogen species in the eastern United States for 1989–1995, Atmos. Environ., 33, 37–49,, 1999. 

Hood, R. R., Shenk, G. W., Dixon, R. L., Smith, S. M. C., Ball, W. P., Bash, J. O., Batiuk, R. A., Boomer, K., Brady, D. C., Cerco, C., Claggett, P., de Mutsert, K., Easton, Z. M., Elmore, A. J., Friedrichs, M. A. M., Harris, L. A., Ihde, T. F., Lacher, L., Li, L., Linker, L. C., Miller, A., Moriarty, J., Noe, G. B., Onyullo, G., Rose, K., Skalak, K., Tian, R., Veith, T. L., Wainger, L., Weller, D., and Zhang, Y. J.: The Chesapeake Bay Program Modeling System: Overview and Recommendations for Future Development, Ecol. Modell., 465, 109635,, 2021. 

Jaworski, N. A.: Retrospective study of the water quality issues of the upper Potomac estuary, Rev. Aquat. Sci., 3, 11–40, 1990. 

Jia, Y., Yu, G., Gao, Y., He, N., Wang, Q., Jiao, C., and Zuo, Y.: Global inorganic nitrogen dry deposition inferred from ground- and space-based measurements, Sci. Rep.-UK, 6, 19810,, 2016. 

Jickells, T., Baker, A. R., Cape, J. N., Cornell, S. E., and Nemitz, E.: The cycling of organic nitrogen through the atmosphere, Philos. T. Roy. Soc. B, 368, 20130115,, 2013. 

Joyce, E. E., Walters, W. W., Le Roy, E., Clark, S. C., Schiebel, H., and Hastings, M. G.: Highly concentrated atmospheric inorganic nitrogen deposition in an urban, coastal region in the US, Environ. Res. Commun., 2, 081001,, 2020. 

Kahmen, A., Renker, C., Unsiker, S. B., and Buchmann, N.: Niche Complentarity for Nitrogen: An Explanation for the Biodiversity and Ecosystem Functioning Relationship?, Ecology, 87, 1244–1255,[1244:NCFNAE]2.0.CO;2, 2006. 

Kanakidou, M., Myriokefalitakis, S., Daskalakis, N., Fanourgakis, G., Nenes, A., Baker, A. R., Tsigaridis, K., and Mihalopoulos, N.: Past, Present and Future Atmospheric Nitrogen Deposition, J. Atmos. Sci., 73, 2039–2047,, 2016. 

Karl, T. R. and Koss, W. J.: Regional and National Monthly, Seasonal, and Annual Temperature Weighted by Area, 1895–1983, National Climatic Data Center, Asheville, NC, 1984. 

Kim, S. W., Heckel, A., McKeen, S. A., Frost, G. J., Hsie, E. Y., Trainer, M. K., Richter, A., Burrows, J. P., Peckham, S. E., and Grell, G. A.: Satellite-observed U.S. power plant NOx emission reductions and their impact on air quality, Geophys. Res. Lett., 33, L22812,, 2006. 

Lamarque, J.-F., Kyle, G. P., Meinshausen, M., Riahi, K., Smith, S. J., van Vuuren, D. P., Conley, A. J., and Vitt, F.: Global and regional evolution of short-lived radiatively-active gases and aerosols in the Representative Concentration Pathways, Climatic Change, 109, 191–212,, 2011. 

Lark, T. J., Meghan Salmon, J., and Gibbs, H. K.: Cropland expansion outpaces agricultural and biofuel policies in the United States, Environ. Res. Lett., 10, 044003,, 2015. 

Li, Y., Schichtel, B. A., Walker, J. T., Schwede, D. B., Chen, X., Lehmann, C. M. B., Puchalski, M. A., Gay, D. A., and Collett, J. L., Jr.: Increasing importance of deposition of reduced nitrogen in the United States, P. Natl. Acad. Sci. USA, 113, 5874–5879,, 2016. 

Likens, G. E., Butler, T. J., and Buso, D. C.: Long- and short-term changes in sulfate deposition: Effects of the 1990 Clean Air Act Amendments, Biogeochemistry, 52, 1–11, 2001. 

Lloret, J. and Valiela, I.: Unprecedented decrease in deposition of nitrogen oxides over North America: the relative effects of emission controls and prevailing air-mass trajectories, Biogeochemistry, 129, 165–180,, 2016. 

Loughner, C. P., Tzortziou, M., Shroder, S., and Pickering, K. E.: Enhanced dry deposition of nitrogen pollution near coastlines: A case study covering the Chesapeake Bay estuary and Atlantic Ocean coastline, J. Geophys. Res.-Atmos., 121, 14221–14238,, 2016. 

McHale, M. R., Ludtke, A. S., Wetherbee, G. A., Burns, D. A., Nilles, M. A., and Finkelstein, J. S.: Trends in precipitation chemistry across the U.S. 1985–2017: Quantifying the benefits from 30 years of Clean Air Act amendment regulation, Atmos. Environ., 247, 118219,, 2021. 

Meyers, T. P., Finkelstein, P., Clarke, J., Ellestad, T. G., and Sims, P. F.: A multilayer model for inferring dry deposition using standard meteorological measurements, J. Geophys. Res.-Atmos., 103, 22645–22661,, 1998. 

Nopmongcol, U., Beardsley, R., Kumar, N., Knipping, E., and Yarwood, G.: Changes in United States deposition of nitrogen and sulfur compounds over five decades from 1970 to 2020, Atmos. Environ., 209, 144–151,, 2019. 

Nowak, J. B., Neuman, J. A., Bahreini, R., Middlebrook, A. M., Holloway, J. S., McKeen, S. A., Parrish, D. D., Ryerson, T. B., and Trainer, M.: Ammonia sources in the California South Coast Air Basin and their impact on ammonium nitrate formation, Geophys. Res. Lett., 39, L07804,, 2012. 

Paulot, F., Jacob, D. J., and Henze, D. K.: Sources and processes contributing to nitrogen deposition: an adjoint model analysis applied to biodiversity hotspots worldwide, Environ. Sci. Technol., 47, 3226–3233,, 2013. 

Paulot, F., Fan, S., and Horowitz, L. W.: Contrasting seasonal responses of sulfate aerosols to declining SO2 emissions in the Eastern U.S.: Implications for the efficacy of SO2 emission controls, Geophys. Res. Lett., 44, 455–464,, 2017. 

Pye, H. O., Luecken, D. J., Xu, L., Boyd, C. M., Ng, N. L., Baker, K. R., Ayres, B. R., Bash, J. O., Baumann, K., Carter, W. P. L., Edgerton, E. S., Fry, J. L., Hutzell, W. T., Schwede, D. B., and Shepson, P. B.: Modeling the Current and Future Roles of Particulate Organic Nitrates in the Southeastern United States, Environ. Sci. Technol., 49, 14195–14203,, 2015. 

Pye, H. O. T., Nenes, A., Alexander, B., Ault, A. P., Barth, M. C., Clegg, S. L., Collett Jr., J. L., Fahey, K. M., Hennigan, C. J., Herrmann, H., Kanakidou, M., Kelly, J. T., Ku, I.-T., McNeill, V. F., Riemer, N., Schaefer, T., Shi, G., Tilgner, A., Walker, J. T., Wang, T., Weber, R., Xing, J., Zaveri, R. A., and Zuend, A.: The acidity of atmospheric particles and clouds, Atmos. Chem. Phys., 20, 4809–4888,, 2020. 

Rao, P., Hutyra, L. R., Raciti, S. M., and Templer, P. H.: Atmospheric nitrogen inputs and losses along an urbanization gradient from Boston to Harvard Forest, MA, Biogeochemistry, 121, 229–245,, 2013. 

Redington, A. L., Derwent, R. G., Witham, C. S., and Manning, A. J.: Sensitivity of modelled sulphate and nitrate aerosol to cloud, pH and ammonia emissions, Atmos. Environ., 43, 3227–3234,, 2009. 

Refsnider, J. M., Garcia, J. A., Holliker, B., Hulbert, A. C., Nunez, A., and Streby, H. M.: Effects of harmful algal blooms on stress levels and immune functioning in wetland-associated songbirds and reptiles, Sci. Total Environ., 788, 147790,, 2021. 

Riddick, S., Ward, D., Hess, P., Mahowald, N., Massad, R., and Holland, E.: Estimate of changes in agricultural terrestrial nitrogen pathways and ammonia emissions from 1850 to present in the Community Earth System Model, Biogeosciences, 13, 3397–3426,, 2016. 

Sabo, R. D., Clark, C. M., Bash, J., Sobota, D., Cooter, E., Dobrowolski, J. P., Houlton, B. Z., Rea, A., Schwede, D., Morford, S. L., and Compton, J. E.: Decadal Shift in Nitrogen Inputs and Fluxes Across the Contiguous United States: 2002–2012, J. Geophys. Res.-Biogeo., 124, 3104–3124,, 2019. 

Schichtel, B. A., Gebhart, K. A., Morris, K. H., Cheatham, J. R., Vimont, J., Larson, R. S., and Beachley, G.: Long-term trends of wet inorganic nitrogen deposition in Rocky Mountain National Park: Influence of missing data imputation methods and associated uncertainty, Sci. Total Environ., 687, 817–826,, 2019. 

Schwede, D. B. and Lear, G. G.: A novel hybrid approach for estimating total deposition in the United States, Atmos. Environ., 92, 207–220,, 2014. 

Sickles II, J. E. and Shadwick, D. S.: Air quality and atmospheric deposition in the eastern US: 20 years of change, Atmos. Chem. Phys., 15, 173–197,, 2015. 

Simon, H., Baker, K. R., and Phillips, S.: Compilation and interpretation of photochemical model performance statistics published between 2006 and 2012, Atmos. Environ., 61, 124–139,, 2012. 

Theobald, M. R., Vivanco, M. G., Aas, W., Andersson, C., Ciarelli, G., Couvidat, F., Cuvelier, K., Manders, A., Mircea, M., Pay, M.-T., Tsyro, S., Adani, M., Bergström, R., Bessagnet, B., Briganti, G., Cappelletti, A., D'Isidoro, M., Fagerli, H., Mar, K., Otero, N., Raffort, V., Roustan, Y., Schaap, M., Wind, P., and Colette, A.: An evaluation of European nitrogen and sulfur wet deposition and their trends estimated by six chemistry transport models for the period 1990–2010, Atmos. Chem. Phys., 19, 379–405,, 2019. 

Townsend, A. R., Howarth, R. W., Bazzaz, F. A., Booth, M. S., Cleveland, C. C., Collinge, S. K., Dobson, A. P., Mallin, M. A., Rogers, C. A., Wayne, P., and Wolfe, A. H.: Human health effects of a changing global nitrogen cycle, Front. Ecol. Environ., 1, 240–246,[0240:HHEOAC]2.0.CO;2, 2003. 

United States Congress: Clean Air Act Amendments of 1990, 1990. 

Walker, J. T., Dombek, T. L., Green, L. A., Gartman, N., and Lehmann, C. M. B.: Stability of organic nitrogen in NADP wet deposition samples, Atmos. Environ., 60, 573–582,, 2012. 

Warner, J., Dickerson, R. R., Wei, Z., Strow, L. L., Wang, Y., and Liang, Q.: Increased atmospheric ammonia over the world's major agricultural areas detected from space, Geophys. Res. Lett., 44, 2875–2884,, 2017. 

Williams, J. J., Chung, S. H., Johansen, A. M., Lamb, B. K., Vaughan, J. K., and Beutel, M.: Evaluation of atmospheric nitrogen deposition model performance in the context of U.S. critical load assessments, Atmos. Environ., 150, 244–255,, 2017. 

World Meteorological Organization: Global Atmosphere Watch Workshop on Measurement-Model Fusion for Global Total Atmospheric Deposition (MMF-GTAD), (last access: 11 March 2022), 2017.  

Wright, C. K., Larson, B., Lark, T. J., and Gibbs, H. K.: Recent grassland losses are concentrated around U.S. ethanol refinerie, Environ. Res. Lett., 12, 044001,, 2017. 

Xing, J., Pleim, J., Mathur, R., Pouliot, G., Hogrefe, C., Gan, C.-M., and Wei, C.: Historical gaseous and primary aerosol emissions in the United States from 1990 to 2010, Atmos. Chem. Phys., 13, 7531–7549,, 2013. 

Xing, J., Mathur, R., Pleim, J., Hogrefe, C., Gan, C.-M., Wong, D. C., Wei, C., Gilliam, R., and Pouliot, G.: Observations and modeling of air quality trends over 1990–2010 across the Northern Hemisphere: China, the United States and Europe, Atmos. Chem. Phys., 15, 2723–2747,, 2015. 

Zhang, L., Jacob, D. J., Knipping, E. M., Kumar, N., Munger, J. W., Carouge, C. C., van Donkelaar, A., Wang, Y. X., and Chen, D.: Nitrogen deposition to the United States: distribution, sources, and processes, Atmos. Chem. Phys., 12, 4539–4554,, 2012. 

Zhang, Y., Mathur, R., Bash, J. O., Hogrefe, C., Xing, J., and Roselle, S. J.: Long-term trends in total inorganic nitrogen and sulfur deposition in the US from 1990 to 2010, Atmos. Chem. Phys., 18, 9091–9106,, 2018. 

Zhang, Y., Foley, K. M., Schwede, D. B., Bash, J. O., Pinto, J. P., and Dennis, R. L.: A Measurement-Model Fusion Approach for Improved Wet Deposition Maps and Trends, J. Geophys. Res.-Atmos., 124, 4237–4251,, 2019. 

Short summary
We assess Community Multiscale Air Quality (CMAQ) model simulations of nitrogen and sulfur deposition over US climate regions to evaluate the model ability to reproduce long-term deposition trends and total deposition budgets. A measurement–model fusion technique is found to improve estimates of wet deposition. Emission controls set by the Clean Air Act successfully decreased oxidized nitrogen deposition across the US; we find increasing amounts of reduced nitrogen to the total nitrogen budget.
Final-revised paper