EGU Journal Logos (RGB) O O O p p p Advances in e Annales e Nonlinear Processes e n n n A A A Geosciencescc Geophysicaecc in Geophysicscc e e e s s s s s s Natural Hazards O Natural Hazards O p p e e and Earth System n A and Earth System n A Sciencesccess Sciencesccess Discussions Atmos.Chem.Phys.Discuss.,13,624A7t–m62o94s,p2h01e3ric O Atmospheric O D www.atmos-chem-phys-discuss.net/13/62C47h/2e0m13i/strypen A Chemistrypen A iscu ACPD d©oAi:1u0th.5o1r(9s4)/2a0c1p3d.-1C3C-6A2t4tr7ib-2u0ti1o3n3.0Laicnends eP.hysicsccess and Physicsccess ssio 13,6247–6294,2013 Discussions n P Atmospheric O Atmospheric O a Thisdiscussionpaperis/hasMbeeeansuunrdeemrerenvtiepen AwforthejournalAtmosphMereicaCsuheremmisetrnytpen A per Multi-model mean c c andPhysics(ACP).PleaserefeTretocthhneiqcoureresscesspondingfinalpaperinACPifTaevcahilnabiqleu.escess | nitrogen and sulfur deposition Discussions Multi-model mean nitrogen and sulfur D deposition fBrioogmeoscitehnceesOpen AAtmospheBrioigceosciencesOpen A iscus J.-F.Lamarqueetal. cc Discussionscc s es es io s s n Chemistry and Climate Model P a TitlePage Climate Ope Climate Ope pe Intercomparison Pron Aject (ACCMIPof t)h:e Pastn A r Abstract Introduction evaluation histoof threi cPaastlccess and projectedDiscussionsccess | Conclusions References D changes Earth System Ope Earth System Ope iscu Tables Figures Dynamicsn Access DynDiascmussiiconssn Access ssion J I P 1 2 3 4 4 4 J.-F. Lamarque , F. Dentener , J. McConnell , C.-U. Ro , M. Shaw , R. Vet , a D. Bergmann5, P. CameronG-Semosictihe5n,tiRfic.OpDoherty6, G. Faluvegi7,GSe.osJc.iGenhtiafinc8Op, per J I Instrumentation e Instrumentation e 9 7 6n 10 7 n B. Josse , Y. H. Lee , I. A. MacKenzie A, D. Plummer , D. T. Shindell , A Back Close D. S. Stevenson6, S. StrodDeMa1et1at,h 1S2o,ydsastn eamdndsGccess. Zeng13 DMaetath Soydsst eamndsccess |D FullScreen/Esc 1NCAREarthSystemLaboratory,NationalCenterforAtmosphericResearchD,isBcuossuioldnser, isc C23EDOue,rsoUepSretAaRnesCeoamrcmhisInssiotMintuo,tJdeoe,iRln DteGRneeoevo,sesNelocaVipr,ecmUnhtSeiCfiAnectnOpen Accesstre,Ispra,Italy Model DGeevoeslocDipiescmnustesiifionncstOpen Access ussionPa PInritnetrearc-tfirvieenDdilsycVuesrssiioonn 4EnvironmentCanada,Toronto,Ontario,Canada p e 5LawrenceLivermoreNationHayldLraobloorgayto aryn,dL iOvermore,CA,USA Hydrology and O r p p e e Earth System6n A247 Earth Systemn A | Sciencesccess Sciencesccess Discussions O O p p Ocean Scienceen A Ocean Scienceen A cce Discussionscce ss ss O O p p Solid Earthen A Solid Earthen A cc Discussionscc e e s s s s O O p p The Cryosphereen A The Cryosphereen A cc Discussionscc e e s s s s 6SchoolofGeoSciences,UniversityofEdinburgh,Edinburgh,UK Dis 7NASAGoddardInstituteforSpaceStudiesandColumbiaEarthInstitute,NewYork,NY,USA cu ACPD s 8PacificNorthwestNationalLaboratory,Richland,WA,USA s io 13,6247–6294,2013 9GAME/CNRM,Me´te´o-France,CNRS–CentreNationaldeRecherchesMe´te´orologiques, n P Toulouse,France a 10CanadianCentreforClimateModelingandAnalysis,EnvironmentCanada,Victoria,British pe Multi-model mean r Columbia,Canada nitrogen and sulfur 11NASAGoddardSpaceFlightCenter,Greenbelt,MD,USA | deposition 12UniversitiesSpaceResearchAssociation,Columbia,MD,USA D 13NationalInstituteofWaterandAtmosphericResearch,Lauder,NewZealand isc J.-F.Lamarqueetal. u s Received:31January2013–Accepted:18February2013–Published:8March2013 sio n Correspondenceto:J.-F.Lamarque([email protected]) P a TitlePage p PublishedbyCopernicusPublicationsonbehalfoftheEuropeanGeosciencesUnion. e r Abstract Introduction | Conclusions References D is c Tables Figures u s s io n J I P a p J I e r Back Close | D FullScreen/Esc is c u s s Printer-friendlyVersion io n P InteractiveDiscussion a p e r 6248 | Abstract D is c ACPD u We present multi-model global datasets of nitrogen and sulfate deposition covering s s time periods from 1850 to 2100, calculated within the Atmospheric Chemistry and ion 13,6247–6294,2013 Climate Model Intercomparison Project (ACCMIP). The computed deposition fluxes P a p 5 are compared to surface wet deposition and ice-core measurements. We use a new er Multi-model mean dataset of wet deposition for 2000–2002 based on critical assessment of the qual- nitrogen and sulfur | ity of existing regional network data. We show that for present-day (year 2000 AC- deposition CMIP time-slice), the ACCMIP results perform similarly to previously published multi- D is model assessments. For this time slice, we find a multi-model mean deposition of c J.-F.Lamarqueetal. u −1 −1 s 50Tg(N)yr from nitrogen oxide emissions, 60Tg(N)yr from ammonia emissions, s 10 and 83Tg(S)yr−1 from sulfur emissions. The analysis of changes between 1980 and ion P 2000 indicates significant differences between model and measurements over the a TitlePage p United States but less so over Europe. This difference points towards misrepresenta- e r Abstract Introduction tion of 1980 NH emissions over North America. Based on ice-core records, the 1850 3 | deposition fluxes agree well with Greenland ice cores but the change between 1850 Conclusions References 15 D and 2000 seems to be overestimated in the Northern Hemisphere for both nitrogen is c Tables Figures and sulfur species. Using the Representative Concentration Pathways to define the u s projected climate and atmospheric chemistry related emissions and concentrations, sio n J I we find large regional nitrogen deposition increases in 2100 in Latin America, Africa P 20 and parts of Asia under some of the scenarios considered. Increases in South Asia ap J I e are especially large, and are seen in all scenarios, with 2100 values more than double r 2000 in some scenarios and reaching >1300mg(N)m−2yr−1 averaged over regional Back Close | to continental scale regions in RCP 2.6 and 8.5, ∼30–50% larger than the values in D FullScreen/Esc anyregioncurrently(2000).ThenewACCMIPdepositiondatasetprovidesnovel,con- is c 25 sistentandevaluatedglobalgriddeddepositionfieldsforuseinawiderangeofclimate us s Printer-friendlyVersion and ecological studies. io n P InteractiveDiscussion a p e r 6249 | 1 Introduction D is c ACPD u Theglobalnitrogencycleisofimportanceforanumberofkey-issues,suchasecology s s and biodiversity (e.g. Phoenix et al., 2006; Bobbink et al., 2010; Butchart et al., 2010), ion 13,6247–6294,2013 eutrophicationandacidification(e.g.Bouwmanetal.,2002;Rodheetal.,2002;Fisher P a p 5 etal.,2011),climatechange-carboncycleinteractions(e.g.Thorntonetal.,2007;Reay er Multi-model mean et al., 2008; Zaehle et al. 2010), food and energy production (Sutton et al., 2011). nitrogen and sulfur | Nitrogen emissions also impact human health through particulate matter and ozone deposition formation. Clearly, nitrogen is central to many aspects of life on Earth (Galloway et al., D is 2008; Fowler et al., 2012). c J.-F.Lamarqueetal. u s Dinitrogen(N )isthemostabundantcomponentintheatmosphere,butischemically s 10 2 io unreactive, only contributing to tropospheric chemistry through lightning and associ- n P ated NO formation (Franzblau and Popp, 1989). More reactive atmospheric nitrogen a TitlePage p components include oxidized NO (=NO+NO +other minor inorganic components e y 2 r Abstract Introduction andorganicnitrogens),thelong-livedgreenhousegasdinitrousoxide(N O;131±10yr, Prather et al., 2012), and reduced nitrogen NH (=NH +NH+) as well2as organic ni- | Conclusions References 15 x 3 4 D trogencomponentssuchasamines(Kanadidouetal.,2012).NOy andNHx arecollec- is c Tables Figures tively identified as reactive nitrogen (N ). Anthropogenic emissions of N components u r r s areestimated(vanAardenneetal.,2001;Lamarqueetal.,2010)tohaveincreasedby sio n J I afactorof5to10since1850.KnowledgeonthelinkbetweenN generationbyhuman r P 20 activities and its subsequent impacts requires an accurate description of largescale ap J I e emissions, atmospheric chemistry transport and removal, all of which occur on rela- r tively fast timescales compared to nitrogen in the other compartments of relevance to Back Close | the Earth system: terrestrial, coastal zones and open ocean. D FullScreen/Esc Globalatmospheric-chemistrytransportmodelsareroutinelyusedasatooltocalcu- is c 25 latetheglobaldispersionofNr.Sinceindividualmodelsarepronetospecificerrorsand uss Printer-friendlyVersion multi-modelmeanstypicallyoutperformindividualmodels(e.g.ReichlerandKim,2008; io n Shindell et al., 2012), it has become common to use model ensemble calculations to P InteractiveDiscussion a improve the quality of the calculations. Previous deposition ensemble studies include p e r 6250 | Lamarqueetal.(2005)using6modelsandfocusingonNO deposition,Denteneretal. D y is (2006) using 23 models and discussing NO , NH and SO deposition under current c ACPD y x x u s and 2 future scenarios, and Sanderson et al. (2008) using 15 models and focusing on s io 13,6247–6294,2013 export of NO and subsequent deposition. n y P The present study uses deposition fields generated by 10 models that participated a 5 p in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP; e Multi-model mean r Lamarque et al., 2012). ACCMIP was designed to inform aspects of the forthcom- nitrogen and sulfur | ing Intergovernmental Panel on Climate Change Assessment Report #5 regarding deposition D the role of long-term changes in atmospheric chemistry on climate and vice-versa. is c J.-F.Lamarqueetal. In this ACCMIP study simulations were provided for time slices representative of the u 10 s periodaround1850,1980,2000,2050and2100(seeSect.2),wherethefutureatmo- sio n spheric chemistry-climate conditions follow the 4 Representative Concentration Path- P ways(Mossetal.,2010).Becauseoftheimportanceofsulfurdeposition(inassociation a TitlePage p e with nitrogen deposition) for the understanding of soil and water acidification (Fisher r Abstract Introduction et al., 2011), we have included an analysis of sulfur deposition to the more traditional 15 | nitrogen deposition analysis. The models used in this analysis represent a combina- Conclusions References D tion of the current generation of Chemistry-Transport Models and Chemistry-Climate is c Tables Figures Models, with somewhat refined horizontal and vertical resolution and more detailed us s descriptions of chemical processes compared to previous studies (namely the results io n J I 20 discussed in Dentener et al., 2006; Sanderson et al., 2008). Rather than represent- P a ing specific meteorological years – as in previous studies – each model meteorology p J I e is rather representative of the decade under consideration (e.g. 1850–1859; or 2000– r 2009) and hence include the effect of climate change. | Back Close Thispaperisstructuredasfollows.Section2givesashortdescriptionofthemodels D FullScreen/Esc 25 and the ACCMIP experiments. We also define there the methodology for computing isc u the multi-model mean (MMM). Section 3 describes the measurement datasets used s s Printer-friendlyVersion in this paper, from surface sites and from ice cores. In Sect. 4, we evaluate the multi- io n model mean for the 2000 period – when extensive measured deposition datasets are P InteractiveDiscussion a available, and which allows us also to make a retrospective analysis of the quality pe r 6251 | of data in this model dataset compared to the earlier and widely used Photocomp D is (Dentener et al., 2006) and HTAP (Sanderson et al., 2008; Vet et al., 2013) deposi- c ACPD u s tion datasets. In Sect. 5, a more limited analysis is given for the changes in deposi- s io 13,6247–6294,2013 tion (and its drivers) since the 1980s and, using ice-cores, for the change between n P decadescenteredaround1850and2000.The1980–2000periodisofparticularinter- a 5 p est since, due to implementation of national and international emission control mea- e Multi-model mean r sures, strong emission reductions have been reported, especially over North America nitrogen and sulfur | andEurope,withcommensurateconsequencesfordeposition.Section6describesthe deposition D overall global structure of total deposition (dry and wet) from 1850 to 2100. A discus- is c J.-F.Lamarqueetal. sion and conclusions follow in Sect. 7. u 10 s s io n P 2 ACCMIP simulations and multi-model mean a TitlePage p e AnoverviewoftheAtmosphericChemistryandClimateModelIntercomparisonProject r Abstract Introduction (ACCMIP) simulations and participating models is given in Lamarque et al. (2012). | Conclusions References Consequently, we focus here on the aspects relevant to the analysis of nitrogen and D is sulfate deposition. c Tables Figures 15 u ACCMIPprovidesanalysisontheroleofatmosphericchemistrychangesinthenear- ss io term-climate forcing included in the CMIP5 (Climate Model Intercomparison Project n J I P Phase 5; Taylor et al., 2012) simulations, including the chemical composition changes a p J I associated with the CMIP5 prescribed forcings. The ACCMIP simulations used in the e r present study (Table 1) consist of time slice experiments (for specific periods span- 20 Back Close | ning 1850 to 2100 with a minimum increment of 10yr) with chemistry diagnostics. Each requested simulation is labeled as Primary (“P”) or Optional (“O”). Although sim- D FullScreen/Esc is ulations using the Representative Concentration Pathway 6.0 (RCP6.0, Masui et al., cu s 2011)emissionswereperformed,notenoughgroupsprovidedthenecessaryfieldsfor s Printer-friendlyVersion io a meaningful analysis of nitrogen and sulfur depositions; therefore, this work focuses n 25 P InteractiveDiscussion on the remaining three RCP projections RCP2.6, RCP4.5 and RCP8.5 (see van Vu- a p uren et al., 2011 and references therein). A primary output of these simulations was er 6252 | nitrogen and sulfate (dry and wet) deposition fields (see Table S1 in Lamarque et al., D is 2012). Although the ACCMIP models were required to specify the anthropogenic and c ACPD u s biomass burning emissions according to Lamarque et al. (2010) for 1850–2000 and s io 13,6247–6294,2013 the RCP emissions (van Vuuren et al., 2011) beyond, there is a range of emissions n P that were used in ACCMIP models, mostly from variations in the treatment of natural a 5 p emissions (Lamarque et al., 2012; Young et al., 2013). e Multi-model mean r As a first step in the quality control of the model calculated depositions, we ana- nitrogen and sulfur | lyze for each model the balance between emission and deposition as reported by the deposition D modeling groups (Tables 2–4) for the year 2000 time slice experiment. This analysis is c J.-F.Lamarqueetal. is done separately for nitrogen originating from nitrogen oxide emissions (deposition u 10 s fields are referred as a group as NOy) or ammonia emissions (NHx) and sulfur emis- sio n sions (SO ). We find that the NO deposition is larger than the emissions (surface and x y P −1 a TitlePage upper-air, including lightning) by approximately 1Tg(N)yr , representing the input of p e nitrogen (mostly in the form of nitric acid) from the stratosphere, in general agreement r Abstract Introduction −1 withobservationalestimatesof0.45Tg(N)yr fromMurphyandFahey(1998),except | 15 −1 Conclusions References for the GISS model where this influx is approximately 5Tg(N)yr . In the case of am- D monium,thebalancebetweenemissionsanddepositionisclearlyattained.Inthecase isc Tables Figures u of sulfur, this balance cannot be fully evaluated due to the lack of diagnostics on the s s formation of sulfate aerosols from the emitted dimethylsulfide (DMS). Boucher et al. ion J I (2003) estimated the yield (DMS to sulfate aerosols) to be 87% when both gas-phase P 20 a and aqueous-phase reactions are taken into account, but this number will be some- pe J I r what model dependent. Within that limitation, it is reasonable to assume that balance Back Close between sulfur deposition and emission is achieved for the listed models. | The focus of this paper is on documenting the multi-model mean (MMM) generated D FullScreen/Esc is for each time slice. The MMM is constructed by linearly interpolating the model gen- c 25 u s erated monthly fields (for example wet deposition combined for all NO species) at s Printer-friendlyVersion ◦ ◦ y io theirnativehorizontalresolution(TableS1)toacommon0.5 ×0.5 grid(finerthanany n P InteractiveDiscussion model grid), identical to the emission grid used in Lamarque et al. (2010). Then, each a p fieldisaveragedacrossmodelsattheoriginalmonthlytemporalresolutiontogenerate e r 6253 | its multi-model mean and standard deviation (Tables S2–S4). As indicated in these ta- D is bles, the largest number of models (up to 9) generated the necessary fields is found c ACPD u s for the historical NO deposition, followed by the RCP8.5 simulations. The number of s y io 13,6247–6294,2013 models performing ammonium chemistry and deposition is however much smaller (2– n P 3 models, depending on the simulated time period) while between 2 and 6 models a 5 p have provided sulfate fields. In all cases, the MMM is constructed using all available e Multi-model mean r model results, therefore leading to variations between the specific models used in the nitrogen and sulfur | average. deposition D is c J.-F.Lamarqueetal. u 3 Description of observational deposition data ss io n Under the auspices of the World Meteorological Organization (WMO) a precipitation P 10 a TitlePage chemistry expert group has performed a critical analysis of available wet deposition p e datafortheyears2000to2002.Whiledrydepositionmaybeavailableatspecificsites, r Abstract Introduction thisisnotadirectlyobservedquantityandwethereforedonotusethisinformation.We | Conclusions References usetheWMO-processedwetdepositiondatasetsinouranalysisoftheperformanceof D is the MMM in the 2000 time slice. c Tables Figures 15 u In the WMO dataset, a careful analysis is made of worldwide reported data of ss io wet precipitation chemistry, and data quality qualifiers are provided to deposition n J I P data obtained mainly from networks in Europe European Monitoring and Evaluation a p J I Programme, EMEP; http://www.nilu.no/projects/ccc/emepdata.html; the United States e r (National Atmospheric Deposition Program, NADP; http://nadp.sws.uiuc.edu/NTN/ 20 Back Close | ntnData.aspx); Canada (Canadian Air and Precipitation Monitoring Network, CAP- MoN; http://www.ec.gc.ca/rs-mn/default.asp?lang=En&n=752CE271-1); Asia (mainly D FullScreen/Esc is fromtheAcidDepositionMonitoringNetworkinEastAsiaEANET;http://www.eanet.cc), cu s and Africa (Deposition of Biogeochemically Important Trace Species (DEBITS), http: s Printer-friendlyVersion io //debits.sedoo.fr). We will focus our analysis on the NADP, EMEP and EANET obser- n 25 P InteractiveDiscussion vations. For more detailed information on data networks, data selection criteria, and a p an application to the HTAP Phase 1 deposition dataset, we refer to Vet et al. (2013), er 6254 | as well as the data networks indicated above. The evaluated wet deposition datasets, D is soon to be available through the World Data Centre for Precipitation Chemistry, only c ACPD u s containasubsetoftheavailablestationsforthatperiod,i.e.stationsthatcorrespondto s io 13,6247–6294,2013 regionallyrepresentativesitesthatfulfilledspecificdatacompletenessandqualitycon- n − P 5 trolmeasures.Inthispaper,weusethewetdepositionmeasurementsofnitrate(NO3), ap ammonium (NH+) and non-sea salt sulfate (SO2−), which were qualified as “good”. er Multi-model mean 4 4 Over Europe (EMEP) and the United States (NADP), deposition measurement sites nitrogen and sulfur | typically started around the 1980s. We use these early measurements to evaluate the deposition D change in deposition over these two decades. To this purpose, we computed from the is c J.-F.Lamarqueetal. raw data (available at the aforementioned web sites), and using the filtering protocols u 10 s s defined by the specific networks, two sets of 6-yr averages (1980–1985 and 1997– io n 2002) for sites that had sufficient observations for both time slices. P a TitlePage Although a quality control similar to the WMO evaluation was not available for the p e 1980s data, comparison of the WMO processed 2000–2002 data with our 1997–2002 r Abstract Introduction 15 indicate very good agreement over the respective regions (not shown), validating the | Conclusions References suitability of our processed dataset to study trends in deposition from the 1980s to the D 2000s. is c Tables Figures u Historical records of nitrate, ammonium, and sulfur were also developed from high- s s depth-resolution measurements in ice cores (see Table S5 for their geographical infor- ion J I mation) using an established continuous ice core analytical system (e.g., McConnell P 20 a and Edwards, 2008). Nitrate and ammonium were measured using spectrophotome- p J I e r try and fluorimetry, respectively, with standard flow through methods (Roethlisberger Back Close et al., 2000). Total sulfur concentrations were measured using magnetic sector Induc- | tively Coupled Plasma Mass Spectrometry (McConnell and Edwards, 2008). At core D FullScreen/Esc siteswithsufficientannualsnowfall,theicecorerecordsweredatedusingannuallayer isc 25 u countingofarangeofseasonallyvaryingchemicalspecies(Sigletal.,2012).Synchro- ss Printer-friendlyVersion io nization to well-known volcanic layers was used for dating at core sites with very low n P InteractiveDiscussion snowfall (e.g., East Antarctica) where annual layers are not preserved (Anschu¨tz et al. a p 2011) or at lower elevation sites (e.g., Akademmi Nauk, McCall Glacier, Flade Isblink) e r 6255 | where surface melting and percolation make annual layer identification difficult (Opel D is et al., 2009). Uncertainty (1σ) in the dating is estimated at ±1yr for the annually dated c ACPD u s sites and ±3yr for all other sites. Decadal averages centered on 1855 and 2000 were s io 13,6247–6294,2013 computed from the high-resolution measurements. n P a p e Multi-model mean r 4 Evaluation of ACCMIP year 2000 time slice deposition 5 nitrogen and sulfur | deposition While observation and model data are available at the monthly time resolution, we fo- D cusouranalysisontheannualmean.Thischoiceismadetolimitthediscussiontothe is c J.-F.Lamarqueetal. u long-termtrendsinnitrogenandsulfatedeposition.Inadditiontothepresentintercom- s s parison project, MMM results from two previous studies are available for comparison io n andanalysisofpotentialimprovements:PhotoComp(Denteneretal.,2006)andHTAP P 10 a TitlePage (Sanderson et al., 2012; Vet et al., 2013). These previous studies are partially inde- p e pendent of the present one, with different emissions and a different sets of models or r Abstract Introduction different versions of the same models. Using the WMO dataset and the model results | Conclusions References interpolated to the location of the observing stations, we can statistically analyze the D is ability of the models to reproduce the observations. c Tables Figures 15 u NHW+e(Fdiigs.p1laby) ianndFiSg.O12−th(Feigg.lo1bca)lwaentdderepgoisointioanl dinisttrhibeuMtioMnMs ocfomNOpa3reNdOto−3 t(hFeigW. 1MaO), ssion J I 4 4 P dataset (over Europe and the United States). Deposition is clearly strongly correlated a p J I withemissions(i.e.largestintheNorthernHemisphere,butalsowithlargeramountsin e r 20 areas of biomass burning and large soil emissions such as Central Africa), albeit with | Back Close significant downwind propagation for each compound. Unlike the nitrogen sources, D FullScreen/Esc sulfur emissions from degassing volcanoes (Andres and Kasgnoc, 1998) lead to the is c formation of deposition hotspots in areas with low anthropogenic emission levels such u s s Printer-friendlyVersion as Central America. io − n In the case of NO deposition, we find that all 3 experiments (ACCMIP, HTAP and 25 3 P InteractiveDiscussion a PhotoComp) perform rather similarly (Fig. 2 and Table 5). However, none of the model p e simulations are able to capture the observed high deposition rates over East Asia or r 6256 |