A new assessment of global and regional budgets, ﬂuxes, and lifetimes of atmospheric reactive N and S gases and aerosols

. We used the EMEP MSC-W (European Monitoring and Evaluation Programme Meteorological Synthesizing Centre – West) model version 4.34 coupled with WRF (Weather Research and Forecasting) model version 4.2.2 meteorology to undertake a present-day (2015) global and regional quantiﬁcation of the concentrations, deposition, budgets, and lifetimes of atmospheric reactive N (N r ) and S (S r ) species. These are quantities that cannot be derived from measurements alone. In areas with high levels of reduced N r (RDN = NH 3 + NH + 4 ), oxidized N r (OXN = NO x + HNO 3 + HONO + N 2 O 5 + NO − 3 + “Other OXN” species), and oxidized S r (OXS = SO 2 + SO 2 −


Introduction
Reactive N (N r ) and S (S r ) gases and particle species are critical determinants of air quality.The principal precursor gas emissions of NO and NO 2 (collectively NO x ), NH 3 , and SO 2 are emitted from a wide range of human activities (Crippa et al., 2018;Hoesly et al., 2018).As the dominant alkaline gas in the atmosphere, NH 3 exerts significant control on the acidity of deposition and on the formation of the secondary inorganic aerosol (SIA) components of fine particles (PM 2.5 ) (Sutton et al., 2020).The latter derive from the reactions between NH 3 and the sulfuric acid (H 2 SO 4 ) and nitric acid (HNO 3 ) that are oxidation products of SO 2 and NO x (Xue et al., 2011;Behera et al., 2013a;Nenes et al., 2020).These particles have significant impacts on human health (Lelieveld et al., 2015;Cohen et al., 2017;Burnett et al., 2018) and on Earth's radiation budget (Bellouin et al., 2011;IPCC, 2021).Sulfur dioxide and sulfate aerosols (collectively referred to as oxidized S and abbreviated as OXS) also have an impact on cloud formation and the acidity of precipitation (Myhre et al., 2017;Aas et al., 2019;Thornhill et al., 2021;McHale et al., 2021).
In addition, oxidized N gases and aerosols such as NO x , HNO 3 , and NO − 3 (collectively abbreviated as OXN), and reduced N gases and aerosol (NH 3 and NH + 4 , collectively abbreviated as RDN), act as powerful plant and microorganism nutrients when deposited to aquatic and terrestrial ecosystems.Excessive anthropogenic input of N can lead to exceedance of N critical loads, which in turn can lead to eutrophication and loss of ecosystem productivity and biodiversity (Erisman et al., 2005;Fowler et al., 2013;Sun et al., 2017).The increase in food production to support a burgeoning human population is increasing agricultural emissions of NH 3 globally, in contrast to the globally declining emissions of SO 2 and NO x (Fowler et al., 2020).A further factor contributing to the different trend in NH 3 emissions is the lack of regulations in most countries specifically targeted at reducing NH 3 emissions.The relative proportions of the N r and S r gas and aerosol species in the atmospheric burdens and the wet and dry deposition of RDN, OXN, and OXS vary substantially around the world because of the different magnitudes of precursor emissions and local meteorology.Knowledge of how these burdens and depositions vary regionally is important for understanding local impacts and mitigation actions.
This information can be derived from atmospheric chemistry transport models.An advantage of these models is that they can simulate concentrations and depositions of air pollutants at greater spatial and temporal resolution and coverage than can be achieved through ambient measurements.Furthermore, as these models simulate the underlying physicalchemical processes, they allow insight into the chemical and meteorological linkages between emissions, concentrations, and deposition of pollutants that measurements cannot reveal: for example, horizontal and vertical pollutant transport as well as the proportion of reduced nitrogen in precipitation that derives from rainout of gaseous NH 3 or particulate NH + 4 .A number of global-scale modelling studies of atmospheric N r and S r have been undertaken previously.Lamarque et al. (2013) presented a multi-model ensemble of N and S deposition from 1850 to 2100, derived from the Atmospheric Chemistry and Climate Model Intercomparison Project (AC-CMIP).The ACCMIP results for "present day" (year 2000 ACCMIP time slice) showed a multi-model mean deposition of 110 TgN yr −1 and 83 TgS yr −1 .Hauglustaine et al. (2014) reported that present day (2000) fine nitrate and ammonium particles make important contributions to aerosol burdens and direct radiative forcing.Bian et al. (2017) conducted an evaluation of global NO − 3 aerosol in 2008 from nine models involved in the Aerosol Comparisons between Observations and Models (AeroCom) phase III and reported a large diversity in global tropospheric burdens of HNO 3 , NO − 3 , NH 3 , and NH + 4 between the models.Tan et al. (2018) calculated the global S and N deposition in 2010 utilizing results of 11 models from the second phase of the Task Force on Hemispheric Transport of Air Pollution (TF HTAP II).The results showed that global S and N depositions in 2010 were 84 TgS yr −1 and 123 TgN yr −1 respectively.Comparisons of their results with those from HTAP I indicated that global distributions of S and N deposition had changed considerably between 2001 and 2010.Globally, there was a reduction in S deposition of 2 TgS (3 %) during this period, with significant reductions in S deposition in Europe, North America, and Russia.In con-
A detailed technical description of EMEP MSC-W rv4.0 is documented in Simpson et al. (2012).An overview of model updates from version rv4.0 to rv4.34 is presented in Simpson et al. (2020).The model has horizontal resolution of 1 • × 1 • and contains 21 terrain-following vertical layers from the surface to 100 hPa.The lowest model layer has a height of ∼ 50 m, and surface concentrations of most species were adjusted to correspond to 3 m above the surface, as described in Simpson et al. (2012Simpson et al. ( , 2020)).
The model aerosol scheme is the Equilibrium Simplified Aerosol Model V4 (EQSAM4clim) (Metzger et al., 2012(Metzger et al., , 2016(Metzger et al., , 2018)), which efficiently parameterizes the aerosol water uptake and full gas-liquid-solid partitioning of mixtures of semi-volatile and non-volatile compounds.In addition to the SO 2− 4 -NO − 3 -NH + 4 scheme, the model includes the reactions between HNO 3 and dust and sea salt to produce coarse nitrate.The aerosol scheme is described in more detail in the Supplement (Sect.S1.1).
Dry deposition of gaseous species is simulated utilizing deposition velocity as described in Simpson et al. (2012Simpson et al. ( , 2020)).Aerosol dry deposition follows a similar formulation to other studies (Lamaud, 1994;Gallagher et al., 1997;Nemitz et al., 2004), with modification for nitrogen compounds (Petroff et al., 2008a, b;Simpson et al., 2012).For wet deposition, all generic fine particles (e.g.SO 2− 4 , NH + 4 , NO − 3 _f) have the same wet scavenging ratio (in cloud) and size-dependent collection efficiency of aerosols by raindrops (below cloud), whilst coarse particles are divided into two groups (coarse sea salt and other coarse particles) with their own sets of parameters (Berge and Jakobsen, 1998;Laakso et al., 2003;Simpson et al., 2012 2021).The terminology "ECLIPSE E " is used here to refer to these monthly emissions.
Dimethyl sulfide (DMS), lightning NO x , soil NO x , and isoprene emissions were configured as reported in Simpson et al. (2017Simpson et al. ( , 2020)), as were the wind-derived emissions of sea-salt particles and dust (Tsyro et al., 2011;Simpson et al., 2012).Emissions from forest and vegetation fires as well as international shipping were also included.et al. (2021) provide a comprehensive evaluation of surface concentrations and wet deposition of N r and S r species from this model configuration against global measurements from 10 monitoring networks.Comparisons of 2010 and 2015 surface concentrations between model and measurement demonstrate that the overall spatial and seasonal variations in the major inorganic pollutants NH 3 , NO 2 , SO 2 , HNO 3 , NH + 4 , NO − 3 , SO 2− 4 , and their wet deposition in East Asia, Southeast Asia, Europe, and North America are captured well by the model.Correlation coefficients between modelled and measured annual mean concentrations of all species are >0.78except for HNO 3 and SO 2 , and they are 0.78 and 0.63 for annual wet deposition of reduced N and oxidized N respectively.Given the intrinsic scale mismatch between a local site measurement and a global-scale chemistry model grid, these comparisons are good and are comparable with model evaluation statistics determined for models of similar resolution (Hauglustaine et al., 2014;Bellouin et al., 2011;Pringle et al., 2010;Xu and Penner, 2012).

Ge
Both model and measurement have uncertainty that constrains the extent to which statistical analyses between mod-elled and measured data can be utilized to assess a model's performance.Model uncertainty comes from uncertainty in emission and meteorology inputs and from parameterizations of chemical/physical processes, whereas measurement uncertainty comes from sampling and analytical procedures.Emission inventory compilation is partially based on reported measurement data and partially on expert estimation, which complicates assigning uncertainties to both emission magnitudes and their spatiotemporal profiles.It is therefore not possible to quantify the contribution of emissions uncertainty to model uncertainty, but it will surely vary regionally and for different species.In Ge et al. (2021), we demonstrated that modelled concentrations of secondary pollutants are less sensitive to the choice of emission inventory than for primary pollutants, as the former are influenced by multiple emissions and the timescales for their formation act to smooth out the differentials in primary emissions.On the other hand, if it is variation in measurement accuracy that contributes to regional variation in model-measurement agreement, then confidence in model output is maintained.Ge et al. (2021) showed better evaluation statistics with measurements in Southeast Asia, Europe, and North America than in East Asia, which suggests shortcomings in regional measurements and/or emissions in the latter network rather than substantial systematic issues with model parameterizations of chemical/physical processes.Ge et al. (2021) also evaluated model response to changes between 2010 and 2015 ECLIPSE E emissions and showed that modelled concentration changes of primary and secondary pollutants and of deposition of reduced N, oxidized N, and oxidized S were consistent with the trends in the corresponding measurements, given realistic uncertainties in emissions and measurements and of meteorological influences between the 2 years.
Nevertheless, considering the fundamental uncertainties in emissions, model, and measurement, all numbers reported in this work should be considered as having underlying uncertainty, albeit that the latest available emission inventory and model version were used to minimize (as far as possible) uncertainties from the first two of these.

Definition of groups of chemical species
The following terminology is used to refer to groups of N r and S r species.RDN represents atmospheric reduced N and is the sum of NH 3 and NH + 4 .OXN represents all oxidized N outputs in the model and consists of HNO 3 , NO x , NO − 3 _f, NO − 3 _c, and "Rest_OXN".The latter comprises N 2 O 5 , HONO, and "Other OXN" species (mainly composed of organic oxidized N species, of which the predominant component is peroxyacetyl nitrate -PAN).Full details are given in Table S1.Total reactive N is the sum of RDN and OXN.OXS represents oxidized S and is the sum of SO 2 and total SO 2− 4 .In this paper, all concentrations, depositions, emissions, and burdens are expressed as mass of the species (e.g. µg m −3 , mg m −2 ) unless otherwise specified (e.g.µgN m −3 , µmol m −3 , TgN yr −1 ).

Definition of world regions
To compare atmospheric budgets and deposition of N r and S r species around the world, we divide the world into the 10 regions illustrated in Fig. 1 (and listed in Table S2).These are based on the 58 reference land and ocean regions defined by the Intergovernmental Panel on Climate Change (IPCC) (Iturbide et al., 2020), each of which represents consistent regional climate features.The rationale behind their definition and their homogeneity is documented in Iturbide et al. (2020).Nine of our regions are continental (albeit with some local surrounding seas); our "Rest of world" region incorporates the majority of the oceans and the remotest continental areas of Antarctica and Greenland.

Results and discussion
3.1 Global N r and S r concentration distributions

Surface distributions
Figure 2 shows the annual mean surface concentrations of RDN, OXN, and OXS and the proportions of each of these that are in the aerosol phase.Maximum concentrations of RDN, OXN, and OXS are observed over the continental regions, particularly eastern China and northern India.These are the regions with the largest population densities and emissions of NH 3 , NO x , and SO x , as shown in Ge et al. (2021).The surface distributions of these N r and S r species are broadly consistent with previous global modelling studies using GEOS-Chem (Zhu et al., 2015), STOCHEM-CRI (Khan et al., 2020), LMDz-INCA (Hauglustaine et al., 2014), and IMPACT (Xu and Penner, 2012).
The maximum RDN concentrations reach 39 µgN m −3 in eastern China; 27 µgN m −3 in northern India; 19 µgN m −3 in southern Europe; and 7-12 µgN m −3 in the eastern United States, south-eastern South America, and eastern Africa.Fertilizer application and livestock farming contribute most to NH 3 emissions in East Asia, South Asia, Europe, and North America, while biomass burning plays an important role in NH 3 emissions in central and eastern Africa (Aksoyoglu et al., 2020;Bray et al., 2018;Dammers et al., 2019;Wang et al., 2018;Whitburn et al., 2015).All of these regions combine both high concentrations of RDN and small proportions of NH + 4 (typically <20 %).Therefore, gaseous NH 3 dominates RDN, indicating that, although reducing NH 3 emissions over these regions will effectively decrease total RDN concentrations, it will have little effect on mitigating SIA formation.A clear increasing gradient of NH + 4 fraction from around 20 % to ∼ 80 % is observed away from these hotspots, indicating conversion of NH 3 to NH + 4 as it disperses from emission sources.The longer lifetime and consequent longer travel distance of NH + 4 compared with NH 3 (see further discussion on this in Sect.3.3) means that NH + 4 accounts for over 90 % of RDN in remote areas, such as over the oceans, deserts, and polar regions.
In general, the concentrations of surface OXN and OXS are also greatest in regions of high RDN concentrations (Fig. 2).This is especially the case over eastern China (maximum of 49 µgN m −3 for OXN and 31 µgS m −3 for OXS) and northern India (maximum of 33 µgN m −3 for OXN and 14 µgS m −3 for OXS).The global maximum OXS concentration (54 µgS m −3 ) surrounds Norilsk in northern Russia where there are very large SO 2 emissions from nickel and palladium mining and smelting operations (Fioletov et al., 2020).As illustration of the magnitude of the SO x emissions from that single model grid (560 Gg yr −1 ), it is double the https://doi.org/10.5194/acp-22-8343-2022Atmos.Chem.Phys., 22, 8343-8368, 2022  S1.
SO x emissions for the whole of the UK (274 Gg yr −1 ) in 2015.
Over continental regions, OXN is predominantly gaseous (i.e.NO x , HNO 3 , HONO, N 2 O 5 , etc.) in areas with both high and low OXN concentrations apart from the remote arid regions of the Sahara, Arabia, and the Tibetan Plateau.The high proportion of gaseous OXN is driven by the widespread emissions of NO x wherever there is population.In contrast, NO − 3 is the major form of OXN over oceans, aside from the major shipping routes (e.g. the sea route from the southern edge of Africa to Southeast Asia) where there are also abundant NO x emissions.The general dominance of aerosol OXN over the oceans is driven both by sea-salt particles acting as an important marine source of coarse NO − 3 (Lee et al., 2014;Schwier et al., 2017) and the reduced depositional sink for gaseous OXN giving time for gas-to-aerosol OXN conversion.
The formation of fine NO − 3 and coarse NO − 3 follows different chemical pathways; therefore, their concentrations and spatial distributions are closely associated with the differing availability of their individual precursors.Figure S1 shows that fine NO − 3 predominates within total NO − 3 over eastern China, Europe, eastern North America, and south-eastern South America, reflecting a relative abundance of gaseous NH 3 in these regions.In contrast, coarse NO − 3 accounts for >80 % of total nitrate over India, Africa, and most of the oceans, which indicates large sources of coarse particles in these regions.Karydis et al. (2016) reported a similar high contribution of coarse NO − 3 to total NO − 3 over India from their modelling of global aerosols.Studies of ambient measurements in India also conclude that the relatively high fraction of coarse nitrate is related to high levels of dust and/or sea salt over this region (Kumar et al., 2010;Boreddy et al., 2021).
OXS is globally dominated by the SO 2− 4 aerosol phase (Fig. 2); SO 2 only accounts for more than a few percent of OXS in the locality of major SO 2 emissions areas such as eastern China, northern Russia, eastern Kazakhstan, and South Africa.This reflects the fast oxidation of SO 2 to SO 2−

Spatial characterization of SO
We apply the concept of chemical domains to characterize the state of ammonia neutralization of both sulfuric and nitric acids (Metzger et al., 2002a, b;Hauglustaine et al., 2014;Xu and Penner, 2012;Ansari and Pandis, 1998).Figure 3 shows the global spatial distribution in relative abundances of SO 2− 4 , fine NO − 3 , and NH + 4 when classified into four such domains.Global maps of the variation in the quantitative species ratios underpinning categorization into the domains shown in Fig. 3 are presented in Fig. S2a for the ratio of total ammonia (T A ) to total sulfate (T S ) and in Fig. S2b for the ratio of free ammonia (T A-free ) to total nitrate (fine; T N-f ).The definitions of these four quantities (in molar concentrations; µmol m −3 ) are given in Eqs. ( 1)-( 4).The SO 2− 4 in T S refers to all of H 2 SO 4 , NH 4 HSO 4 , (NH 4 ) 3 H(SO 4 ) 2 , and (NH 4 ) 2 SO 4 .
The T A-free term quantifies the amount of total ammonia that is in excess of that needed to fully neutralize the sulfate (i.e. the T A remaining when T A T S >2).When T A-free T N-f > 1, sufficient ammonia also exists to fully react with nitrate to form NH 4 NO 3 .An arbitrary threshold for the sum of the concentrations of fine SO 2− 4 , NO − 3 , and NH + 4 to be ≥ 0.5 µg m −3 has been applied to Fig. 3 in order to focus attention away from the most remote areas (Southern Hemisphere oceans and the poles) where the fine SIA concentration is very low and knowledge of the chemical domain is irrelevant for informing emissions mitigation strategies.
In the blue and yellow areas in Fig. 3, T A T S < 2 and sulfuric acid is not fully neutralized; thus, not all SO 2− 4 is in the form of NH 4 (SO 4 ) 2 .Locations where this ratio is between 1 and 2 are coloured blue in Fig. 3 (green colours in Fig. S2a and grey areas in Fig. S2b) and are defined as "sulfate rich".This is the state over Northern Hemisphere oceans and remote continental regions such as the Sahara Desert, Greenland, and Siberia.Over northern Russia and Southern Hemisphere oceans, where T A T S <1 (yellow colour in Fig. 3), conditions can be characterized as "sulfate very rich".As discussed in Sect.3.1.1,northern Russia is an SO 2 emission hotspot, but NH 3 emissions are small.This is also the case around other very large point sources of anthropogenic SO 2 emissions in Kazakhstan, southern Africa, and southern Australia.Over ocean areas, biogenic emission of dimethyl sulfide (DMS) is the major source of sulfur in marine aerosol (Fiddes et al., 2018;Jackson et al., 2020), but NH 3 emission is again very small.A separation between "sulfate rich" and "sulfate very rich" regions arises because there is transport of anthropogenic continental NH 3 emissions over the North-ern Hemisphere oceans in the former areas.In these "sulfate rich" and "sulfate very rich" regions, essentially all NH 3 has been converted to NH + 4 , as illustrated in Fig. 2 where RDN concentrations are almost entirely in the form of NH + 4 .The green and orange areas in Fig. 3 (yellow to red in Fig. S2a) have T A T S >2 and, consequently, T A-free T N-f >0 (non-grey colours in Fig. S2b) and show where sulfuric acid is fully neutralized and excess NH 3 remains to react with HNO 3 to form NH 4 NO 3 .Most continental regions with significant SIA pollution are characterized by these conditions, including East Asia, Southeast Asia, South Asia, Europe, North America, south-eastern South America, western Africa, and eastern Australia.
The green areas in Fig. 3 have a T A-free T N-f ratio between 0 and 1 and indicate where NH 3 is the limiting factor for NH 4 NO 3 formation and, thus, where the chemical domain can be characterized as "NO − 3 rich".The blue colours in Fig. S2b illustrate the quantitative extent of the HNO 3 neutralization in these areas.This situation occurs over coastal areas around East Asia, northern Europe, and the north-eastern United States.These locations experience the transport of oxidized N from the adjacent continent and/or direct shipping emissions of NO x , but NH 3 concentrations are lower because it deposits more quickly to land before transport off the coast.
The orange areas in Fig. 3 are where T A-free T N-f >1 and are characterized as "NH 3 very rich", as the availability of nitrate limits NH 4 NO 3 formation in these areas.This condition prevails over most of the temperate and tropical land surface in the Northern Hemisphere.Concentrations of SIA over a lot of the land in the Southern Hemisphere are (like over the Southern Hemisphere oceans) less than the 0.5 µg m −3 threshold applied to Fig. 3. Figure S2b shows that the most "NH 3 rich" areas are in Central Africa, South Asia, East Asia, Southeast Asia, western Europe, South and Central America, North America, and New Zealand where there are intensive agricultural activities and/or large biomass burning and/or low NO x and SO x emissions.Hauglustaine et al. (2014) and Xu and Penner (2012) reported similar findings, although both of these studies also showed "NH 3 very rich" conditions over tropical oceans, which may be related to oceanic NH 3 emissions if these were included in their simulations.They were not included in our model.Oceanic NH 3 emissions are small (∼ 6 TgN yr −1 , albeit highly uncertain; Fowler et al., 2015), compared with continental NH 3 emissions (53.1 TgN yr −1 ; Table 1), and we deliberately exclude discussion of the chemical domain where the SIA concentration is small.
The continental values of T A-free T N-f calculated here for 2015 are larger than those reported for 2000.For instance, in western Europe, Xu and Penner (2012) and Hauglustaine et al. (2014) reported T A-free T N-f ratios of around 1.0-2.0 and 1.0-4.0respectively, whereas we report values of around 4.0-10.0(Fig. S2b).This temporal difference is driven not only by the increasing NH 3 emissions associated with sustaining a food supply for the growing global population (Sutton et  al., 2013;Behera et al., 2013b;Aneja et al., 2009) but also by the fact that a growing number of countries have been reducing NO x and SO x emissions while most have little or no control policy on NH 3 emissions (Liu et al., 2018;Liu and Matsui, 2021;Aas et al., 2019;de Gouw et al., 2014).These changes in relative emissions of precursors mean an increasingly smaller proportion of gaseous NH 3 is transformed to aerosol NH + 4 , which in turn increases NH 3 concentrations and reduces the effectiveness of PM 2.5 mitigation via NH 3 emission reductions.For example, Liu et al. (2018) reported a significant growth in gaseous NH 3 concentrations over the North China Plain from 2008 to 2016 allied to the rapid SO 2 emission reductions.
Our results also indicate regional differences in terms of the extent of excess ammonia (Fig. S2b).Table S3 lists the maximum, median, and area-weighted mean T A-free T N-f ratios in different world regions.Among the regions with the largest RDN surface concentration, South Asia shows the most "NH 3 rich" condition, with an area-weighted mean ratio of 9.6, followed by South America and Africa, whose areaweighted mean ratios are 8.6 and 7.2 respectively.For regions with an area-weighted mean T A-free T N-f ratio >1, Southeast Asia and Euro_Medi have the lowest area-weighted mean ratios of 1.9 and 2.5 respectively.This lower relative excess of NH 3 of the former is associated with the large proportion of marine area in this region, while the smaller NH 3 abundance in the latter is likely linked to more stringent NH 3 emission targets in Europe than in other world regions via the Gothenburg Protocol and the EU National Emissions Ceilings Directive.
3.2 Global and regional N r and S r deposition An important advantage of an atmospheric model is that it can apportion the different contributions to deposition that measurements cannot readily quantify.For instance, direct measurements of dry deposition are technically difficult and are very limited in terms of both spatial and temporal coverage.Furthermore, whilst measurement of NH + 4 in precipitation permits quantification of wet deposition of RDN, it cannot distinguish how much of it derives from rainout of gaseous NH 3 and how much derives from rainout of particle NH + 4 .Similarly, measurements cannot distinguish between gaseous and particle sources to wet deposition of OXN and OXS.The models therefore provide additional information about the nature of N r and S r deposition.This is relevant because the form of deposited N affects the severity of nitrogen deposition impacts (Sutton et al., 2014;Bobbink et al., 2010).An important effect of deposited N gases and aerosols is the direct toxicity to above-ground parts of various plants.Van Herk et al. (2003) found that lichens are particularly sensitive to direct NH 3 dry deposition, whilst Pescott et al. (2015) reported direct toxic effects of wet-deposited RDN for bryophytes.Additionally, it must be stressed that the lack of dry deposition comparison between model and measurement, due to the lack of measurement data, means that the numbers reported for modelled dry deposition in this section should be regarded as more uncertain.

RDN
Figure 4 shows the global distribution of annual RDN deposition in 2015 separately apportioned into the dry and wet deposition of NH 3 and NH + 4 .The map of total RDN deposition (the sum of all panels of Fig. 4) is provided in Fig. S3.The relative contributions of the different sources of RDN deposition are illustrated in Fig. 5 for the 10 world regions outlined in Fig. 1. Figure 4a-d show similar spatial patterns of the highest depositions over East Asia, South Asia, Europe, eastern North America, south-eastern South America, and central Africa.The greatest RDN deposition (wet + dry) rate occurs in East Asia, exceeding 5300 mgN m −2 yr −1 .The East Asia region receives more than one-fifth of the world's RDN deposition (Fig. 5).The maximum total RDN deposition in South Asia and Southeast Asia is in the range of 3800 to 4600 mgN m −2 yr −1 , whilst in Europe, North America, South America, and Africa maximum total RDN deposition is in the range 1800 to 2200 mgN m −2 yr −1 (Fig. S3).These regions each receive 7 %-11 % of global total RDN deposition (Fig. 5).Around 15 % of RDN deposition occurs over the oceans (Rest of world in Fig. 5).In fact, as most of the other world regions also include areas of water, the proportion of global RDN depositing to seas will be greater than 15 %.
Globally, the total RDN deposition is 55.5 TgN yr −1 , apportioned approximately equally (50 %) between dry and wet deposition.The global deposition can be further apportioned as 40 % NH 3 dry deposition, 9 % NH + 4 dry deposition, 22 % NH 3 wet deposition, and 29 % NH + 4 wet deposition (Table S4).Thus, whilst the large majority (82 %) of RDN dry deposition is in the form of NH 3 , the RDN wet deposition is more evenly split between rainout of gaseous NH 3 (44 %) and rainout of particulate NH + 4 (56 %).Overall, 62 % and 38 % of total global RDN is removed from the atmosphere initially in the form of NH 3 and NH + 4 respectively.However, Fig. 5 shows that there are large regional variations in the contributions of different sources of RDN deposition.In South Asia, South America, and Oceania, more than 50 % of RDN deposition is in the form of NH 3 dry deposition.In East Asia, Euro_Medi, North America, and Africa, dry deposition of NH 3 is in the range of 42 %-49 %.Wet deposition of NH + 4 is dominant only in remote areas such as the Rest of world (54 %) and Rest of Asia (40 %) regions.Overall, the deposition of NH 3 (dry + wet) accounts for well over half of total RDN deposition in all continental regions (except for Rest of Asia), which indicates that reducing NH 3 emissions will be effective for reducing total RDN input from the atmosphere to continental ecosystems.

OXN
Figure 6 shows the global distribution of annual OXN deposition in 2015 apportioned into dry and wet deposition of gaseous OXN (i.e.NO x + HNO 3 + HONO + N 2 O 5 + Other OXN species) and aerosol-phase OXN (NO − 3 _f + NO − 3 _c).The map of total OXN deposition is given in Fig. S4.As described above for RDN deposition, model simulations are able to distinguish in detail the ultimate atmospheric source of the OXN deposition, which measurements cannot; in particular, to distinguish whether dissolved nitrate in precipitation derives from gas-or aerosol-phase species as well as which species.The relative contributions of the different sources of OXN deposition for the 10 world regions outlined in Fig. 1 are shown in Fig. 7.
The regions of high OXN deposition in Fig. 6 are similar to those with high RDN deposition (Fig. 4), although shipping NO x emissions can be seen to have a notable im-pact on deposition of OXN over a large number of marine locations, which is not apparent for RDN deposition.The largest amounts of total OXN deposition are in East Asia (3300 mgN m −2 yr −1 ), followed by Southeast Asia (2800 mgN m −2 yr −1 ) and South Asia (1700 mgN m −2 yr −1 ) (Fig. S4).However, whilst the East Asia region receives 22 % of global RDN deposition (Fig. 5), this region receives only 14 % of global OXN deposition (Fig. 7).In North America, maximum OXN deposition is around 1900 mgN m −2 yr −1 , whereas it is around 900 mgN m −2 yr −1 in Europe.A greater proportion of OXN deposition occurs over the oceans than is the case for RDN deposition (Rest of world proportions are 23 % and 15 % respectively), due both to the shipping emissions of NO x and the longer lifetime (hence, greater transport distance) of NO 2 compared with NH 3 .
The global total OXN deposition in 2015 is 57.9 TgN yr −1 , apportioned as 43 % dry deposition and 57 % wet deposition.The relative contributions of the individual OXN components to global OXN deposition are shown in Fig. S6.The global deposition is further apportioned as 28 % gaseous OXN dry deposition, 15 % aerosol-phase OXN dry deposition, 18 % gaseous OXN wet deposition, and 39 % aerosol-phase OXN wet deposition (Table S4).The global OXN dry deposition is almost equally split among HNO 3 (34 %), coarse NO − 3 (29 %), and NO x (21 %, this is effectively dry deposition of NO 2 ), whereas the rainout of coarse particulate NO − 3 (60 %) is the predominant form in OXN wet deposition, followed by the rainout of HNO 3 (28 %).The dry-deposited fine NO − 3 and the rainout of fine NO − 3 only account for 7 % and 9 % of OXN dry deposition and OXN wet deposition respectively.The total deposition (wet + dry) of Other OXN species (predominantly PAN) contributes 5 % to global total OXN deposition, and the contributions from HONO and N 2 O 5 are negligible.Compared with global RDN deposition, the majority of which derives from wet and dry deposition of gaseous NH 3 (62 %), Fig. S6 shows that it is coarse particulate NO − 3 (46 %) that contributes most to global OXN deposition, followed by HNO 3 (30 %).Overall, slightly greater OXN deposition derives from aerosol species (54 %) than from gaseous species (46 %) globally.
Figure 7 shows that the 10 world regions can be viewed in four groups according to the nature of their OXN deposition.A group comprising Southeast Asia, South America, and Oceania is characterized by large contributions to total OXN deposition derived from wet and dry deposition of gaseous species.In Southeast Asia, wet and dry deposition of HNO 3 constitutes over half (52 %) of total OXN deposition, with coarse NO − 3 (31 %) and NO x (7 %) being the second and third contributors.This indicates that it is the formation of gaseous species (mainly HNO 3 and NO 2 ) that is the most important determinant of OXN deposition in this region.Similar observations are made for South America and Oceania.Fine NO   In the regions of South Asia, Rest of Asia, Africa, and Rest of world, it is coarse NO − 3 (via both wet and dry deposition) that is the most significant contributor to total OXN deposition (wet + dry).This is consistent with the relatively greater abundance of coarse particulate NO − 3 in these regions (Fig. S1).
A third group of regions from Fig. 7 comprises North America and Euro_Medi.In these regions, the proportion contributions to total OXN deposition from wet and dry deposition of HNO 3 and coarse NO − 3 are similar at 31 %-38 %, and there is also significant contribution from dry deposition of NO x (13 %-14 %) and Other OXN (6 % and 3 % for North America and Euro_Medi respectively).The relative contribu-  S1).Note the logarithmic scale.
tions from NO x in these regions are well above their average contribution globally.There is also an important contribution (8 %-10 %) from wet and dry deposition of fine NO − 3 to total OXN deposition in these two regions.
Finally, the East Asia region is characterized by a substantial contribution (24 %) to total OXN deposition from wet and dry deposition of fine NO − 3 , a proportion that is more than double the fine NO − 3 contribution in North America and Euro_Medi.The formation of NH 4 NO 3 aerosol is therefore an important influence on OXN deposition over East Asia, in addition to its important contribution to PM 2.5 and to the latter's health effects.The atmospheric components HNO 3 and coarse NO − 3 contribute an equal amount (30 %) to total wet and dry OXN deposition in East Asia.

OXS
Figure 8 shows the global distribution of annual OXS deposition apportioned into wet and dry deposition of SO 2 and SO 2 4 .
The map of total OXS deposition is provided in Fig. S5.The relative contributions of the different sources of OXS deposition for the 10 world regions outlined in Fig. 1 are shown in Fig. 9. Higher levels of OXS deposition are found near emission sources in East Asia, South Asia, Europe, the eastern United States, and south-eastern Africa for all four deposition components.In addition, both wet and dry deposition of SO 2 is high close to major metal smelting industry such as in northern Russia, where total OXS deposition exceeds 5400 mgS m −2 yr −1 ; eastern Kazakhstan; and South Africa.
In East Asia and South Asia, the highest total (wet + dry) OXS depositions are around 4800 and 2400 mgS m −2 yr −1 respectively.OXS deposition in Europe and North America is much smaller, with the largest values in the range from 900 to 1100 mgS m −2 yr −1 .A substantial proportion of total OXS deposition occurs over the oceans; 34 % occurs over the Rest of world region (Fig. 9), which reflects both deposition of S derived from oxidation of natural DMS emissions and the longer transport distances of anthropogenic aerosol SO 2− 4 .Of the continental regions, Fig. 9 shows that the greatest proportions of OXS deposition occur in the East Asia (16 %) and Rest of Asia (11 %) regions.
The total deposited OXN (TgN) in each region and its fraction of global total OXN deposition are presented below each region name.
tion contributes 11 %, and dry deposition of aerosol SO 2− 4 contributes 10 %.This means that, globally, slightly greater OXS deposition derives from aerosol SO 2− 4 (54 %) than from SO 2 (46 %).This is a similar situation to that for OXN deposition but contrasts with RDN deposition, for which 62 % ultimately derives from the gaseous component (NH 3 ) and only 38 % derives from the aerosol component (NH + 4 ).The total OXS deposition that we calculate for 2015 is around half that reported by Hauglustaine et al. (2014) for 2000 (107 TgS yr −1 ) and 40 % lower than that reported by Tan et al. (2018) for 2010 (84 TgS yr −1 ).It is also lower than the values for 2000 (range 75-100 TgS yr −1 ) from a 23model intercomparison reported by Dentener et al. (2006).The lower OXS deposition value presented here is assumed to be due in part to the considerable reduction in global sulfur emissions from 2000 to 2015 (Fowler et al., 2020;Aas et al., 2019).A direct comparison of modelled S deposition between our work and other studies is confounded by the use of different inventories for both anthropogenic and natural emissions.
As for RDN and OXN, the relative contributions of the different sources of OXS deposition vary across the 10 world regions (Fig. 9).However, the regional composition patterns of OXS deposition are broadly similar to those of RDN deposition, with dry deposition of SO 2 (23 %-46 %) and wet deposition of SO 2− 4 (34 %-56 %) being the two most important contributors in all world regions.The total depositions (wet + dry) of SO 2 and SO 2− 4 contribute almost equal proportions to total OXS deposition in Southeast Asia, South Asia, North America, South America, and Oceania.In the East Asia and Euro_Medi regions, the sum of SO 2 deposition contributes slightly more (55 %-56 %) to total OXS deposition, whereas it is SO 2− 4 that contributes most (57 %-67 %) to total OXS deposition in the Rest of Asia, Africa, and Rest of world regions.
In summary, on a global basis, wet deposition and dry deposition currently contribute equally to RDN deposition.For OXN and OXS, wet deposition contributes slightly more than dry deposition.These proportions vary regionally according to local precursor emission strengths and climatology.For RDN deposition, dry deposition of NH 3 is the predominant contributor in most regions except for Rest of Asia and Rest of world (essentially oceans) where NH 3 emissions are small and RDN deposition is largely determined by the transport of aerosol NH + 4 from other regions (see Sect. 3.2.1).Further reductions in NH 3 would therefore efficiently reduce the deposition of RDN in most world regions.The two largest contributors to OXN deposition in all regions are HNO 3 and coarse NO − 3 (both occurring by both wet and dry deposition).Dry deposition of NO x (i.e.NO 2 ) is the third largest contributor in most regions except for East Asia and Rest of world.The deposition of fine NO − 3 is only important over East Asia and to a lesser extent in the Euro_Medi region.Due to the general dominance of aerosol-phase NO − 3 in OXN deposition, changes in primary NO x emissions may not have a direct and linear influence on the deposition of OXN.Sensitivity experiments focusing on the response of deposition to emission changes are required to gain more insight into this issue.For OXS deposition, SO 2 and SO 2− 4 contribute almost equally to total OXS deposition in most regions.More specifically, the sum of dry deposition of SO 2 and wet deposition of SO 2− 4 accounts for at least 67 % of total OXS deposition in all world regions.
3.3 Global and regional N r and S r budgets

Global budgets
Table 1 summarizes global tropospheric budget data for RDN, OXN, and OXS for 2015.Table 2 does likewise for the individual components NH 3 , NH + 4 , HNO 3 , fine NO − 3 , coarse NO − 3 , and SO 2− 4 .In this study, the top of the troposphere is defined as the top of our model domain at 100 hPa.Table 1 shows that net exchange of chemical species at the top of the model domain ("Flux in" and "Flux out") is negligible.We calculate atmospheric lifetimes based on both sources and sinks.The detailed calculation procedure is described in the Supplement (Sect.S1.2).
The only source of atmospheric RDN is the annual NH 3 surface emission of 53.1 TgN.The annual average burden of RDN is 0.75 TgN, and it has a global lifetime of 4.9-5.2d depending on how this is calculated (Table 1).The total deposition of NH 3 is 34.6 TgN yr −1 (Table 2); this contributes 62 % of the sink of NH 3 and occurs as 65 % dry deposition and 35 % wet deposition.Partitioning of NH 3 into NH + 4 aerosol, which is also the only source of NH + 4 , contributes 21.0 TgN yr −1 (38 %) to the total sink of NH 3 .The NH + 4 loss comes mainly from wet deposition, 15.8 TgN yr −1 (76 %), and to a lesser extent from dry deposition, 5.02 TgN yr −1 (24 %).Net loss to the stratosphere is negligible.The global average tropospheric lifetime of NH + 4 (8.9 d) is much longer than that of NH 3 (1.6 d), which is consistent with the annual average NH + 4 burden (0.51 TgN) being more than double that of NH 3 (0.24 TgN) (Table 2).
Xu and Penner ( 2012 It is important to note that the NH 3 burden across the nine models in Bian et al. (2017) varied by a factor of 17.This indicates the sensitivity of burdens and lifetimes of NH 3 and NH + 4 to the definition of model domain, the deposition pathways, and the interactions between NH 3 and acidic species in different models.
The global atmospheric OXN source and sink fluxes are 56.6 and 57.9 TgN yr −1 respectively (Table 1).The annual average OXN burden is 1.20 TgN, with a global mean lifetime of 7.6 or 7.7 d depending on how this is calcuhttps://doi.org/10.5194/acp-22-8343-2022Atmos.Chem.Phys., 22, 8343-8368, 2022 Table 1.Global tropospheric budgets of total RDN, OXN, and OXS in 2015.Flux (out -in) = Flux out -Flux in.Q is the annual average burden, P is the production pathway (i.e.sources), R is the removal pathway (i.e.sinks), and τ is atmospheric lifetime.The constituents of the three groups are listed in Table S1.
τ = Burden Sources 5.2 7.7 5.9 Lifetime (days): τ = Burden Sinks 4.9 7.6 5.6 lated.Detailed global budgets of HNO 3 and fine and coarse NO − 3 are listed in Table 2.As quantification of the flux of HNO 3 through photolysis and reaction with OH is not available as model output, the HNO 3 budget presented here includes only deposition and aerosol nitrate formation.Xu and Penner (2012) and Hauglustaine et al. (2014) report that 6.9 %-10 % of HNO 3 loss occurs via photolysis and reaction with OH, so these loss routes are relatively small.The total chemical production of HNO 3 in our study is 49.5 TgN yr −1 , which is 17 % higher than the 42.3 TgN yr −1 reported by Xu and Penner (2012) and 2 % higher than the 48.5 TgN yr −1 in Hauglustaine et al. (2014).The total loss of HNO 3 is 49.4 TgN yr −1 , of which 17.6 TgN yr −1 (36 %) is removed from the atmosphere by deposition, almost equally via dry deposition (8.35 TgN yr −1 ) and wet deposition (9.25 TgN yr −1 ).The remaining 64 % of the HNO 3 loss goes into formation of 4.54 TgN yr −1 fine nitrate (14 % of the aerosol nitrate source) and 27.2 TgN yr −1 coarse nitrate (86 % of the aerosol nitrate source).The tropospheric annual average burden of HNO 3 is 0.15 TgN with a lifetime of 1.1 d.
The loss of fine and coarse nitrate is only via deposition, and 72 % of total nitrate loss is due to wet deposition.The calculated nitrate lifetimes range between 4.8 and 7.1 d, which are in line with lifetimes (4.6-6.3 d) reported in other studies (Xu and Penner, 2012;Hauglustaine et al., 2014;Bian et al., 2017).Coarse nitrate (7.0-7.1 d) has a slightly longer lifetime than fine nitrate (4.8 d).Our calculated total nitrate annual average burden is 0.58 TgN, with fine NO − 3 comprising only 10 % of this total.This total nitrate burden is around 3 times greater than the value in Xu and Penner (2012) and Hauglustaine et al. (2014) for 2000 (0.17-0.18TgN), which likely reflects a shift in gas-aerosol partitioning between 2000 and 2015 from gas-phase HNO 3 more towards aerosol-phase NO − 3 , especially given that their HNO 3 burden of 0.30 TgN is double our value.It is also possible that the simulation of nitrate burdens is more sensitive to coarse-particle emissions and the production of coarse nitrate in our model than in the models used by Xu and Penner (2012) and Hauglustaine et al. (2014).
The dominance of coarse nitrate in the global total nitrate burden is also observed in other global modelling studies.Hauglustaine et al. (2014), Xu andPenner (2012), andBauer et al. (2007) reported a relative contribution of coarsemode nitrate of 72 %, 77 %, and 79 % to the total nitrate burden in their respective models for 2000.Additionally, the lifetime of fine nitrate calculated in our work is comparable with other literature values.Park et al. (2004) reported a global fine-nitrate lifetime against deposition of 3.2 d for 2001, Feng et al. (2007) reported a global lifetime of total nitrate (fine + coarse) of 5.0 d for 1997, and Alexander et al. (2020) reported a global lifetime of inorganic nitrate (HNO 3 + total NO − 3 ) for the 2000-2015 period on the order of 3-4 d.
For OXS, emission of SO 2 is the only source, and deposition (wet + dry) is almost the only sink pathway (Table 1).The annual average burden of OXS is 0.78 TgS with a lifetime of 5.6-5.9d depending on how it is calculated.The total source of SO 2− 4 is 27.6 TgS yr −1 (Table 2).The dominant loss of SO 2− 4 is wet deposition (81 %) with the remaining loss (19 %) being dry deposition.The annual average SO 2− 4 burden in the troposphere is 0.63 TgS.The SO 2− 4 lifetime of 8.3-8.5 d is very similar to the NH + 4 lifetime, which can be explained by the close chemical association between these two species in atmospheric SIA.
Differences in species' lifetimes underpin the differences in the vertical and horizontal distributions of the species.For example, as discussed in Sect.3.1.1,gaseous NH 3 , with a short lifetime, only shows a dominant contribution to the total RDN surface concentration in places near emission sources, whereas aerosol-phase NH + 4 , as a RDN reservoir species with a longer lifetime, is the major form of RDN over the remote continental regions and oceans.
The vertical distributions of N r and S r species shown in Fig. 10 also reflect the different species' lifetimes.All species exhibit strong vertical gradients of decreasing concentration from the surface up to 10 km.The gaseous species (NH 3 , NO 2 , and SO 2 ) have peak concentrations at the surface between 20 and 40 • N, as shown in their surface distribution maps (Sect.3.1.1).These gaseous species are emitted at the surface, and their lifetimes are too short for significant mixing to high altitudes.Between the model surface layer (∼ 50 m) and ∼ 2.5 km, the maximum zonal mean of NH 3 and SO 2 decreases from 1.77 to 0.31 µg m −3 and from 1.78 to 0.23 µg m −3 respectively, whereas the concentration gradient of NO 2 is steeper, with the maximum zonal mean decreasing from 4.11 to 0.10 µg m −3 , reflecting that production and removal processes of NO 2 mainly occur in the surface layer.
In contrast, the secondary species, NH + 4 , NO − 3 (fine + coarse), and SO 2− 4 , are mixed to higher altitudes than their precursors.The zonal mean concentrations of the secondary species also show a larger horizontal variation range (10-60 • N), again reflecting their longer transport times.A notable observation from Fig. 10 is that the maximum concentrations of the secondary aerosol species are not at the surface, which is a feature that will not be apparent from surface measurements of these species.The maximum zonal mean of NH + 4 (1.20 µg m −3 ) and SO 2− 4 (2.58 µg m −3 ) occurs in the model level from 88 to 177 m above the surface, whereas the maximum NO − 3 (3.80 µg m −3 ) occurs in the 177-313 m model level.These elevated maxima in secondary components reflect their lack of surface emission and the time taken for the in-atmosphere oxidation and gas-to-particle chemistry.

Regional budgets
Figure 11 shows the differences between annual emissions and depositions for RDN, OXN, and OXS in the 10 world regions illustrated in Fig. 1.The corresponding absolute emission and deposition fluxes for each region are presented in Fig. S7.
For continental regions, the maximum fluxes for RDN are in East Asia (12.5 TgN yr −1 for RDN emission and 12.0 TgN yr −1 for RDN deposition), followed by South Asia (7.67 TgN yr −1 emission and 5.46 TgN yr −1 deposition) (Fig. S7).Both of these regions are characterized by intensive agricultural activities and have RDN fluxes that are larger than the OXN fluxes.Fluxes of RDN in the North America, Euro_Medi, and Southeast Asia regions are lower (3.46-5.14TgN yr −1 ).The spatial distribution of the RDN concentrations presented in Fig. 2 are the consequence of these regional variations in RDN fluxes.Figure 11 shows that all regions are net exporters of RDN except for Southeast Asia and Rest of world, in which marine areas dominate.The net import of RDN in Southeast Asia is, however, very small.In contrast, deposition of RDN in the Rest of world region is nearly 16 times its annual emissions.The difference of 7.89 TgN yr −1 shows that RDN deposition is a notable source of continental anthropogenic pollution to the oceans.The region with greatest export of RDN is South Asia, at 2.21 TgN yr −1 or 29 % of its annual emissions (Fig. 11).Euro_Medi and Africa have the next greatest net exports (1.08 TgN yr −1 (21 % of https://doi.org/10.5194/acp-22-8343-2022Atmos.Chem.Phys., 22, 8343-8368, 2022 τ = Burden Sources 1.6 8.9 1.1 4.8 7.0 8.3 Lifetime (days): emissions) and 0.95 TgN yr −1 (15 % of emissions) respectively).The relative differences between emissions and depositions over North America, South America, and Oceania are in the range of 6.0 %-8.0 %.Although East Asia is the region with the greatest RDN emissions, it only exports 3.8 % (0.48 TgN yr −1 ).This region is geographically generally self-contained from the other regions, but this means that East Asia is largely responsible for its own RDN pollution.
The relatively high availability of acidic species also leads to shorter chemical lifetimes in this region.
For OXN, East Asia and Africa are the two continental regions with the largest annual emission and deposition fluxes (6.72-9.59TgN yr −1 ), followed by South America and North America (5.49-7.02TgN yr −1 ) (Fig. S7).In all of Southeast Asia, Euro_Medi, North America, South America, Africa, Oceania, and Rest of world, the OXN fluxes have larger in-fluence on their total nitrogen budgets than the RDN fluxes, whereas, as noted above, in East Asia and South Asia fluxes of RDN comprise the majority.In contrast to our finding of higher RDN flux (12.0-12.5 TgN yr −1 ) than OXN flux (8.33-9.59TgN yr −1 ) in East Asia, Tan et al. (2018) report higher OXN fluxes (10.5-11.4TgN yr −1 ) than RDN fluxes (7.7-8.5 TgN yr −1 ).Such a difference in RDN and OXN fluxes is partially attributed to differences in emission inventories, model structures, meteorology input, and region definition between the two studies.It may also reflect an increase in NH 3 emissions and a reduction in NO x emissions in this region between 2015 and 2010.
Figure 11 shows that the continental region with the greatest net export of OXN (1.92 TgN yr −1 , 22 % of emissions) is Africa, followed by Euro_Medi (1.61 TgN yr −1 , 26 %), South Asia (1.57TgN yr −1 , 34 %), and South America (1.51 TgN yr −1 , 22 %).In particular, OXN deposition over Rest of world (13.4 TgN yr −1 ) is around 5 times larger than its emission (2.21 TgN yr −1 , from shipping).The difference of 11.2 TgN yr −1 between OXN deposition and emission in this region is greater than the 7.89 TgN yr −1 difference for RDN.It is likely that the oxidation of NO x to reservoir organic nitrate species such as PAN facilitates the longdistance transport of OXN from the continent to the ocean, whereas the short lifetime of NH 3 means less RDN is transported from its continental emission regions.Regardless of the explanation, the very large net deposition of OXN in the Rest of world region again indicates the scale of anthropogenic N deposition to marine areas.
Considering OXS, Fig. S7 shows that the Rest of world region has both the largest emission (10.6 TgS yr −1 ) and deposition (17.4 TgS yr −1 ) fluxes.The former is due to the substantial marine DMS emissions whose atmospheric oxidation ultimately produces SO 2 and H 2 SO 4 (Marandino et al., 2013;Schobesberger et al., 2013;Hoffmann et al., 2016).These OXS emissions largely deposit back to the oceans, but Fig. 11 shows that OXS deposition in the Rest of world is also substantially augmented by long-range transport of OXS pollution from other regions.For continental regions, East Asia has the largest OXS fluxes, with emission and deposition reaching 8.95 and 8.30 TgS yr −1 respectively, followed by Rest of Asia (5.48-6.21TgS yr −1 ) and Euro_Medi (4.05-5.57TgS yr −1 ) (Fig. S7).South Asia is the largest net exporter of OXS with 1.62 TgS yr −1 (37 % of its emission) being deposited outside this region (Fig. 11).Euro_Medi is the second largest net exporter of OXS (1.52 TgS yr −1 , 27 %).All other continental regions are also either positive or negligible exporters of OXS except for Southeast Asia, which receives more OXS deposition than it emits.This is similar to the situation for RDN because, away from major NO x emission areas, NH + 4 is predominantly associated with SO 2− 4 .The tropospheric burden of each species is around 1 or 2 orders of magnitude smaller than its corresponding emission/deposition flux.Figures 12 and S8 respectively show the relative and absolute contributions of individual RDN and OXN species to the total tropospheric N burden and the contributions of the individual OXS species to the total tropospheric S burden, for each of the 10 world regions.In Table 3, we present the regional atmospheric lifetimes of RDN, OXN, and OXS as quantified against their total emission or deposition fluxes.These regional lifetimes are not fully accurate because they do not take the net transport of RDN, OXN, or OXS across the boundary of the region into account.However, given that, aside from the Rest of world region, the net transport in/out of a region is relatively small compared with the emission and deposition fluxes (Fig. 11), they should be reasonable estimates of the regional RDN, OXN, and OXS lifetimes.Table 3 shows that the regional RDN, OXN, and OXS lifetimes vary by up to a factor of 4 across the continental regions, generally being longest over Africa and Rest of Asia, where population and emissions are generally less dense, and shortest over East Asia, Southeast Asia, and South Asia.
The Rest of world region has the highest burdens of all regions with 0.23 TgN RDN, 0.44 TgN OXN, and 0.27 TgS OXS, representing 31 %, 37 %, and 35 % of the global RDN, OXN, and OXS burdens respectively, despite the relatively small emissions for RDN and OXN from this region.This observation is consistent with the conclusions from analyses of regional emission/deposition fluxes in Fig. 11 that the ocean areas are substantially impacted by the transport of N r and S r air pollutants from continental regions.The lifetimes of RDN, OXN, and OXS in Rest of world are also the longest of all regions (Table 3), reflecting that RDN, OXN, and OXS in this region predominantly exist as longer-lived species: SO 2− 4 in the case of S r and NH + 4 , NO − 3 , and Rest OXN (predominantly PAN) in the case of N r .East Asia, Rest of Asia, and Africa have moderate burdens of 0.09-0.11TgN (12 %-14 % of the global burden) for RDN, 0.10-0.16TgN (9 %-13 %) for OXN, and 0.08-0.12TgS (11 %-16 %) for OXS.South Asia, Euro_Medi, and North America have smaller burdens with their contributions to global burdens ranging from 6 % to 9 % for RDN, from 4 % to 7 % for OXN, and from 5 % to 8 % for OXS.
The species contributions to N burdens differ between the regions (Fig. 12).RDN contributes 58 % and 49 % to the total N burden in South Asia and East Asia respectively, which is similarly reflected in the emission and deposition fluxes in the two regions, followed by Rest of Asia (42 %) and Euro_Medi (39 %).The largest relative contributions of OXN to the total N burden are in Oceania (82 %) and South America (70 %).
The largest contribution of NH 3 to the RDN burden is in South America (67 %), and this large proportion of NH 3 gives South America the shortest RDN lifetime (1.8-1.9 d) of any region.In other parts of the world, NH   S1.The corresponding emission and deposition fluxes are shown in Fig. S7.
For OXN, the contribution of HNO 3 to the OXN burden reaches 28 %, 24 %, and 22 % in Southeast Asia, Oceania, and South America respectively, whereas HNO 3 contributes only 5 %-12 % to OXN in other parts of the world (Fig. 12).NO − 3 _f only makes a significant contribution (20 %) to the total OXN burden over East Asia; its contribution is <6 % in other regions.This reflects the strong propensity for NH 4 NO 3 formation in East Asia, which is similarly shown in the composition pattern of OXN deposition in this region as well (Sect.3.2.2).Coarse NO − 3 is a significant contributor to the OXN burden, particularly in South Asia (59 %), Rest of Asia (54 %), and Africa (52 %).This is associated with an abundance of coarse particles like sea salt, dust, and sand in these regions.Therefore, it is critical to correctly account for contributions of coarse nitrate in models, as it overwhelmingly controls the atmospheric nitrate burden and determines the size distribution of nitrate aerosol.In all world regions, SO 2− 4 is the predominant form of the OXS burden (Fig. 12).
To mitigate atmospheric N r and S r burdens and their associated impacts on human health and ecosystems, tackling secondary inorganic species like NH + 4 as well as coarse NO 3and their precursors should be prioritized.The impacts of possible emission controls on atmospheric N r and S r concentrations, depositions, budgets, and lifetimes can be quantified by emission sensitivity tests.
The model evaluation presented in Ge et al. (2021) shows that the model-measurement bias varies from one pollutant to another.For instance, the model underestimates annual global HNO 3 surface concentration measurement by 34 % but overestimates NO 2 by 23 %.Although it is difficult to diagnose the extent to which bias is contributed by bias intrinsic to the model (including emissions input to the model), by bias in measurements, or by model-measurement spatial representation mismatch, it is important to acknowledge that all numbers reported in this work may be biased somewhat low or high.
The intention here is to provide a new overview of how N r and S r budgets vary among different world regions and among different chemical species.As Ge et al. (2021) show that the EMEP model is capable of capturing the overall spatial variation in surface concentrations for major N r and S r pollutants and their wet deposition, and as our results are comparable with similar global modelling studies, we have confidence that we characterize the major regional differences in budgets.However, we also acknowledge that there are varying uncertainties in different regions, as revealed by our model evaluation study (Ge et al., 2021).Therefore, we have greater confidence in regional budgets in North America, Europe, and Southeast Asia than other world regions.Nevertheless, a modelling study cannot resolve whether different regional uncertainties derive from the emissions, model processes, measurements, or a combination of these aspects.Further evaluation studies of emissions and measurements are required to gain more insight.

Conclusions
We used the EMEP MSC-W model 4.34 coupled with WRF 4.2.2 meteorology to undertake a new present-day evaluation of global and regional concentrations, wet and dry deposition, fluxes, and lifetimes of reactive N and S species.These are quantities that cannot be derived from measurements alone.Simulations used the 2015 ECLIPSE annual emissions with EDGAR monthly profiles.Surface concentration and wet deposition outputs from this model configuration have been comprehensively evaluated against several measurement networks.S1).
Table 3. Regional atmospheric lifetimes (in days) for RDN, OXN, and OXS, calculated using either total emission or deposition fluxes.Q is the annual average burden, P is the production pathway (i.e.emissions), R is the removal pathway (i.e.depositions), and τ is atmospheric lifetime.These regional lifetimes are not fully accurate because they do not consider the net transport of RDN, OXN, or OXS across the boundary of the region, although these are relatively small for most continental regions.Please note that the lifetimes derived from the production pathway for Rest of world are not listed in this table because these numbers are highly skewed by the net transport of RDN, OXN, and OXS from other regions (full data are available in Table S5).Our simulations show that the largest surface concentrations and deposition of RDN (NH 3 + NH + 4 ), OXN (NO x + HNO 3 + HONO + N 2 O 5 + NO − 3 _f + NO − 3 _c + Other OXN species), and OXS (SO 2 + SO 2− 4 ) are in eastern China and northern India.In these areas, RDN concentrations are predominantly in the form of NH 3 (NH + 4 is typically <20 %), as also in Europe and North America.NH + 4 comprises >90 % of RDN over oceans and remote continental areas.Gaseous OXN species constitute the majority of OXN over continental regions, whereas NO − 3 dominates OXN over oceans, except for major shipping routes.In eastern China, Europe, and eastern North America, fine NO − 3 dominates nitrate, but coarse NO − 3 dominates in India, as has been noted elsewhere.Aerosol SO 2− 4 is the major form of OXS nearly everywhere except in the vicinity of very large SO 2 emission sources.
Most continental regions are now characterized as "ammonia rich"; this indicates that whilst reducing NH 3 emissions will decrease the RDN concentration and deposition, it will have little effect on mitigating secondary inorganic aerosol.Our results show an increase in ammonia richness compared with model simulations for 2000 in other studies.This is driven not only by increasing NH 3 emissions in parts of the world but also by the fact that a growing number of countries https://doi.org/10.5194/acp-22-8343-2022Atmos.Chem.Phys., 22, 8343-8368, 2022 have been reducing NO x and SO x emissions much faster than NH 3 , for which many countries have yet to set emissions targets.In general, South Asia shows the most ammonia-rich condition.Coastal areas around East Asia, northern Europe, and the north-eastern United States can be characterized as "nitrate rich" where NH 4 NO 3 formation is limited by NH 3 not nitrate.These locations experience the transport of oxidized N from the adjacent continent and/or direct shipping emissions of NO x , but NH 3 concentrations are lower because it deposits more quickly to land before transport off the coast.
Our modelling shows that OXN (57.9 TgN yr −1 , 51 %) and RDN (55.5 TgN yr −1 , 49 %) account for almost equal proportions of the total nitrogen deposition.Total sulfur (OXS) deposition is 50.5 TgS yr −1 .Globally, wet and dry deposition contribute equally to RDN deposition.For OXN and OXS, wet deposition contributes slightly more than dry deposition; however, the contributions of different deposition components vary regionally.An advantage of modelling is that it can distinguish the individual deposition processes, which measurements cannot.For example, the model shows that the majority of the wet deposition of RDN derives from rainout of gaseous NH 3 rather than of particle NH + 4 in many regions.For RDN deposition, dry deposition of NH 3 is the predominant contributor in most regions except for Rest of Asia and Rest of world (essentially oceans), where NH 3 emissions are small and RDN deposition is largely determined by the transport of aerosol NH + 4 from other regions.Thus, further reductions in NH 3 would efficiently reduce the deposition of RDN in most world regions.The two largest contributors to OXN deposition in all regions are HNO 3 and coarse NO − 3 (both occurring by both wet and dry deposition).The deposition of fine NO − 3 is only important over East Asia (24 %) and is <10 % in other regions.For OXS deposition, dry-deposited SO 2 and wet-deposited SO 2− 4 are two largest contributors in all regions.
The annual average tropospheric burden of RDN is 0.75 TgN of which NH 3 and NH + 4 comprise 32 % (0.24 TgN) and 68 % (0.51 TgN) respectively.The average tropospheric lifetime of NH 3 is 1.6 d, whereas that of NH + 4 is more than 5 times longer (8.9 d).The annual average burden of HNO 3 is 0.15 TgN with an average lifetime of only 1.1 d.The total nitrate burden is 0.58 TgN with fine NO − 3 only constituting 10 % of this total.The lifetime of fine nitrate (4.8 d) is shorter than that of coarse nitrate (7.0-7.1 d).The global lifetime of  is longer than that of RDN (4.9-5.2 d), which is consistent with a total OXN burden (1.20 TgN) almost double that of RDN.The tropospheric burden of OXS is 0.78 TgS with a global lifetime of 5.6-5.9d.The SO 2− burden is 0.63 TgS, and its lifetime of 8.3-8.5 d is similar to that of NH + 4 , reflecting the close association between the two species in the atmosphere globally.It is notable that lifetimes of RDN, OXN, and OXS vary by a factor of 4 across different continental regions.East Asia has quite short lifetimes for RDN (2.9-3.0 d), , whereas the lifetimes in Rest of Asia and Africa are around twice as long, which is related to the regional differences in meteorological conditions and the compositions of RDN, OXN, and OXS burdens.
The model confirms that the longer lifetimes of the secondary N r and S r species are consistent with their greater vertical transport, which surface measurements cannot reveal.In addition, whilst the gaseous precursor species have maximum concentrations at the surface, maximum concentrations of the aerosol species are at higher altitudes because of the time taken for their formation.
Comparison of the absolute emission and deposition fluxes of RDN and OXN in different world regions shows that East Asia and South Asia are characterized by larger RDN fluxes than OXN fluxes, whereas most other world regions exhibit larger OXN fluxes than RDN fluxes, which likely reflects greater increases in NH 3 emissions in the two regions.
Comparison of the difference between emission and deposition fluxes shows that the South Asia region is the largest net exporter of RDN (2.21 TgN yr −1 , 29 % of its annual emission), followed by Euro_Medi, Africa, and other continental regions.Despite having the largest RDN emissions and deposition, East Asia has very small net export and is therefore largely responsible for its own RDN pollution.In contrast, Rest of world, which is mainly the ocean regions (as well as Southeast Asia which also contains a lot of marine area), receives more deposition of RDN than it emits.
For OXN, Africa is the largest net exporter of OXN at 22 % (1.92 TgN yr −1 ) of its emission, followed by South Asia.The Rest of world again receives substantial net import of OXN (11.2 TgN yr −1 ) which is even higher than the 7.89 TgN yr −1 net import of RDN.Oxidation of NO x to reservoir organic nitrate species such as PAN likely facilitates the long-distance transport of OXN from the continent to the ocean, whereas the short lifetime of NH 3 means that less RDN is transported from its continental emission regions.The large flux imbalances of OXN and RDN in these areas indicate the scale of anthropogenic N pollution to marine areas.
The two regions with the largest net export of OXS are South Asia (1.62 TgS yr −1 , equivalent to 37 % of its emission) and Euro_Medi (1.52 TgS yr −1 , 27 %).All other continental regions are either positive or zero exporters of OXS except for Southeast Asia, which receives more OXS deposition than it emits.This is similar to the situation for RDN because, away from major NO x emission areas, NH + 4 is predominantly associated with SO 2− 4 .The three largest contributors to the total N burden globally are NH + 4 , coarse NO − 3 , and Rest OXN (mainly PAN).In South Asia, Rest of Asia, and Africa, coarse NO − 3 constitutes over half of the OXN burden, attributable to an abundance of coarse sea-salt and dust particles in these regions.It is critical to account for the contributions of coarse nitrate, as it controls the atmospheric nitrate burden and the size distribution of nitrate aerosol.Aerosol SO 2− 4 is the predominant form of OXS in all regions.
As presented in this work, N r and S r species are abundant air pollutants impacting gas and particle concentrations and deposition globally.The chemical and meteorological linkages between emissions and N r and S r species' concentrations are complex, and the detail can only be revealed by process-based atmospheric modelling, as illustrated here, despite the uncertainties in emissions and model parameterizations.

Figure 1 .
Figure 1.The boundaries of the 10 world regions used in this study (based on the IPCC reference regions described in Iturbide et al., 2020).

Figure 2 .
Figure 2. (a, c, e) Annual mean surface concentrations (log scales) and (b, d, f) inorganic aerosol-phase fractions of (a, b) RDN, (c, d) OXN, and (e, f) OXS species.The full composition of OXN is listed in TableS1.

Figure 3 .
Figure 3.The state of ammonia neutralization of sulfuric and nitric acids for areas where the SIA concentration in the fine-particle fraction (PM 2.5 ) is ≥ 0.5 µg m −3 .Definitions of the four chemical domains are given in the main text.Maps of the quantitative species ratios underpinning categorization into these chemical domains are given in Fig. S2.

Figure 4 .
Figure 4. Annual wet deposition (WDEP) and dry deposition (DDEP) of RDN components for 2015."NH 3 WDEP" and "NH + 4 WDEP" refer to wet deposition of RDN originally present in the atmosphere as gaseous NH 3 or as aerosol NH + 4 respectively.Note the logarithmic scale.

Figure 5 .
Figure 5. Relative contributions to annual RDN deposition for the set of world regions defined in Fig. 1.The two darker colours collectively indicate wet deposition (WDEP), whereas the two lighter colours collectively indicate dry deposition (DDEP).The total deposited RDN (TgN) in each region and its fraction of global total RDN deposition are presented below each region name.

Figure 7 .
Figure 7. Relative contributions to annual OXN deposition for the set of world regions defined in Fig. 1.Labels for individual contributions <2 % are omitted to aid clarity.The darker colours collectively indicate wet deposition (WDEP), whereas the lighter colours collectively indicate dry deposition (DDEP).The "Other" OXN species category predominantly comprises PAN (full compositional details in TableS1).The total deposited OXN (TgN) in each region and its fraction of global total OXN deposition are presented below each region name.

Figure 8 .
Figure 8. Annual wet deposition (WDEP) and dry deposition (DDEP) of OXS components for 2015."SO 2 WDEP" and "SO 2− 4 WDEP" refer to wet deposition of OXS originally present in the atmosphere as gaseous SO 2 or as aerosol SO 2− 4 respectively.Note the logarithmic scale.
), Hauglustaine et al. (2014), and Bian et al. (2017) report similar global total NH 3 emissions (50.5-53.6TgN yr −1 for 2000 and 51.8 TgN yr −1 for 2008) but smaller burdens and shorter lifetimes of NH 3 and NH + 4 than we report here.The respective burdens of NH 3 and NH + 4 are 0.07 TgN (lifetime of 0.46 d) and 0.26 TgN (lifetime of 3.23 d) in Xu and Penner (2012) and 0.09 TgN (lifetime of 0.63 d) and 0.22 TgN (lifetime of 4.52 d) in Hauglustaine et al. (2014) for 2000.Burdens of NH 3 and NH + 4 averaged across nine models in Bian et al. (2017) for 2008 are 0.16 and 0.25 TgN with a corresponding lifetime of 0.72 and 4.3 d respectively.

Figure 9 .
Figure 9. Relative contributions to annual OXS deposition for the set of world regions defined in Fig. 1.The two darker colours collectively indicate wet deposition, whereas the two lighter colours collectively indicate dry deposition.The total deposited OXS (TgS) in each region and its fraction of global total OXS deposition are presented below each region name.

Figure 11 .
Figure 11.Absolute differences (ADs) and relative differences (RDs) between emission and deposition fluxes of RDN (TgN yr −1 ), OXN (TgN yr −1 ), and OXS (TgS yr −1 ) for the 10 world regions defined in Fig. 1.Relative differences are expressed relative to emission fluxes.Difference values for the Rest of world region are much larger than for other regions.The constituents of RDN, OXN, and OXS are listed in TableS1.The corresponding emission and deposition fluxes are shown in Fig.S7.

Figure 12 .
Figure 12.Relative compositions of the tropospheric burdens in 2015 (stacked bar plots, left y axis) of total N (RDN + OXN) and total S (OXS) for the 10 world regions defined in Fig. 1.The total N (RDN + OXN) and S regional burdens are shown as the line plots (right y axis).Rest OXN includes HONO, N 2 O 5 , and Other OXN species (full list given in TableS1).

Table 2 .
Global tropospheric budgets of NH 3 , NH + 4 , HNO 3 , NO − 3 , and SO 2− 4 in 2015.Flux (out -in) = Flux out -Flux in.Q is the annual average burden, P is the production pathway (i.e.sources), R is the removal pathway (i.e.sinks), and τ is atmospheric lifetime.
+ 4 is the majority contributor to the RDN burden.The RDN lifetime in Southeast Asia and East Asia is relatively short (2.6-3.0 d), whereas Rest of Asia has the longest RDN lifetime (7.0-7.3 d) except for Rest of world.