Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ Atmospheric SRef-ID: 1680-7324/acp/2006-6-447 Chemistry European Geosciences Union and Physics
Dry and wet deposition of inorganic nitrogen compounds to a tropical pasture site (Rondonia,ˆ Brazil)
I. Trebs1, L. L. Lara2, L. M. M. Zeri3, L. V. Gatti4, P. Artaxo5, R. Dlugi6, J. Slanina7, M. O. Andreae1, and F. X. Meixner1 1Max Planck Institute for Chemistry, Biogeochemistry Department, P.O. Box 3060, 55020, Mainz, Germany 2Centro de Energia Nuclear na Agricultura (CENA), Laboratorio de Ecologia Isotopica,´ Universidade de Sao˜ Paulo (USP), Av. Centenario, 303 13400-970, Piracicaba, Sao˜ Paulo, SP, Brazil 3Max Planck Institute for Biogeochemistry, Department Biogeochemical Processes, Hans-Knoll-Straße¨ 10, 07745 Jena, Germany 4Instituto de Pesquisas Energeticas´ e Nucleares, CQMA, Atmospheric Chemistry Laboratory, Av. Prof. Lineu Prestes, 2242, Cidade Universitaria, CEP 055508-900, Sao˜ Paulo, SP, Brazil 5Instituto de F´ısica, Universidade de Sao˜ Paulo (USP), Rua do Matao,˜ Travessa R, 187, CEP 05508-900, Sao˜ Paulo, SP, Brazil 6Working Group Atmospheric Processes (WAP), Gernotstrasse 11, 80804 Munich, Germany 7Peking University, College of Environmental Sciences, Beijing 100871, China Received: 28 February 2005 – Published in Atmos. Chem. Phys. Discuss.: 19 May 2005 Revised: 22 November 2005 – Accepted: 13 December 2005 – Published: 8 February 2006
Abstract. The input of nitrogen (N) to ecosystems has grass surface appeared to have a strong potential for daytime increased dramatically over the past decades. While total NH3 emission, owing to high canopy compensation points, (wet + dry) N deposition has been extensively determined which are related to high surface temperatures and to direct in temperate regions, only very few data sets of N wet de- NH3 emissions from cattle excreta. NO2 also significantly position exist for tropical ecosystems, and moreover, reliable accounted for N dry deposition, whereas HNO3, HONO and experimental information about N dry deposition in tropical N-containing aerosol species were only minor contributors. environments is lacking. In this study we estimate dry and Ignoring NH3 emission from the vegetation surface, the an- wet deposition of inorganic N for a remote pasture site in nual net N deposition rate is estimated to be about −11 kgN the Amazon Basin based on in-situ measurements. The mea- ha−1 yr−1. If on the other hand, surface-atmosphere ex- surements covered the late dry (biomass burning) season, a change of NH3 is considered to be bi-directional, the annual transition period and the onset of the wet season (clean condi- net N budget at the pasture site is estimated to range from tions) (12 September to 14 November 2002) and were a part −2.15 to −4.25 kgN ha−1 yr−1. of the LBA-SMOCC (Large-Scale Biosphere-Atmosphere Experiment in Amazonia – Smoke, Aerosols, Clouds, Rain- fall, and Climate) 2002 campaign. Ammonia (NH3), nitric 1 Introduction acid (HNO3), nitrous acid (HONO), nitrogen dioxide (NO2), + nitric oxide (NO), ozone (O3), aerosol ammonium (NH4 ) The supply of reactive nitrogen (N) to global terrestrial − and aerosol nitrate (NO3 ) were measured in real-time, ac- ecosystems has doubled since the 1960s as a consequence companied by simultaneous meteorological measurements. of human activities, such as fertilizer application, cultivation Dry deposition fluxes of NO2 and HNO3 are inferred us- of N fixing legumes and production of nitrogen oxides by ing the “big leaf multiple resistance approach” and particle fossil-fuel burning (Galloway, 1998). The deposition of at- deposition fluxes are derived using an established empiri- mospheric N species constitutes a major nutrient input to the cal parameterization. Bi-directional surface-atmosphere ex- biosphere. On a long-term scale, the increase of N inputs change fluxes of NH3 and HONO are estimated by applying into terrestrial ecosystems may result in (i.) intensified trace a “canopy compensation point model”. N dry and wet depo- gas exchange (ii.) enhanced leaching of nitrate and soil nu- + + 2+ 2+ sition is dominated by NH3 and NH4 , which is largely the trients (e.g., K , Mg , Ca ), (iii.) ecosystem eutrophi- consequence of biomass burning during the dry season. The cation and acidification, (iv.) reduction in biodiversity, and (v.) increased carbon storage (Vitousek et al., 1997). En- Correspondence to: I. Trebs hanced carbon storage due to N deposition has been shown ([email protected]) to increase the terrestrial carbon sink in N-limited temperate
© 2006 Author(s). This work is licensed under a Creative Commons License. 448 I. Trebs et al.: N deposition in the tropics ecosystems, which may have substantial impacts on global rural locations (e.g., North Canada) up to −30 to −70 kgN −1 −1 CO2 concentrations (e.g., Townsend et al., 1996). ha yr for regions that receive N from urban pollution N deposition is considered to be relevant in the tropics and/or agricultural activities (e.g., North Sea, NW Europe due to widespread biomass burning activity and increasing and NE U.S.) (Howarth et al., 1996). In contrast to mod- fertilizer application. It was suggested by e.g., Matson et erately fluctuating air pollution levels that prevail in Eu- al. (1999) and Asner et al. (2001) that in contrast to temper- rope, the United States and Asia throughout the year, tropical ate ecosystems, nitrogen-rich/phosphorus (P)-limited trop- environments such as the Amazon Basin experience every ical rainforest soils may have a reduced productivity fol- year a dramatic change from the “green ocean” clean back- lowing excess N deposition, resulting in a decreased C- ground atmosphere to extremely polluted conditions during storage. Moreover, the humid tropical zone is a major the biomass burning season. Only few studies exist were source area for biogenic nitrous oxide (N2O) and nitric oxide atmospheric N wet removal was determined experimentally (NO) emissions from soils (Reiners et al., 2002). Enhanced in the tropics (Clark et al., 1998; Galloway et al., 1982; N inputs to tropical forests are likely to increase nitrifica- Likens et al., 1987; Srivastava and Ambasht, 1994). The tion/denitrification rates and, hence, the emission of NO and chemical composition of precipitation in the Amazon re- N2O to the atmosphere (Hall and Matson, 1999). The con- gion was determined in previous studies by Andreae et al. version of tropical rainforest into cultivated land and pasture (1990), Stallard and Edmond (1981), Lesack and Melack may lead to a sustained disturbance of the natural N cycle. (1991) and Williams et al. (1997). Galy-Lacaux et al. (2003) During clearing and burning of tropical rainforest, biomass- have estimated N wet + dry deposition for different (trop- associated N is volatilized and a large fraction is emitted in ical) African ecosystems. However, their dry deposition form of gaseous NH3 (Trebs et al., 2004), which may result measurements are based on the application of integrating in considerable N losses of tropical ecosystems (Kauffman et filter methods which are prone to artifacts (cf. Slanina et al., 1998; Kauffman et al., 1995). This is affirmed by the fact al., 2001). Since 1999, N wet + dry deposition are mon- that, in contrast to old growth forests, plant growth in defor- itored on the Caribbean Virgin islands (St. John Island; ested areas is suggested to be limited by N rather than by P see http://www.epa.gov/castnet///sites/vii423.html). Global (Davidson et al., 2004; Oliveira et al., 2001). chemistry and transport models (CTMs) such as MOGUN- The deposition of atmospheric N compounds occurs via TIA have been applied to estimate total N deposition on a dry and wet processes. Nitrogen dioxide (NO2), ammonia global scale (e.g., Holland et al., 1999). Model results sug- (NH3), nitric acid (HNO3) and nitrous acid (HONO) are the gest that net N deposition in the northern temperate latitudes most important contributors to N dry deposition. HNO3 usu- exceeds contemporary tropical N deposition by almost a fac- ally features a rapid downward (net deposition) flux to the tor of two. However, reliable experimental information about surface (Huebert and Robert, 1985). By contrast, the ex- N dry deposition in tropical environments, required to vali- change of NO, NH3, HONO and NO2 between surface and date these model predictions, has been lacking. atmosphere may be bi-directional. The rates of production Kirkman et al. (2002) determined the surface-atmosphere and consumption in vegetation elements and/or soils as well exchange of NOx (NO + NO2) and ozone (O3) at a pasture as the ambient concentration determine whether net emis- site in the Amazon Basin (Rondonia,ˆ Brazil). In this paper, sion or net deposition of these species takes place. Turbu- we complement their study by additionally estimating the lent diffusion controls the transport of gases and particles surface-atmosphere exchange fluxes of NH3, HNO3, HONO, − + from the surface layer to the Earth’s surface. The uptake aerosol NO3 and NH4 at the same pasture site. Our analysis of trace gases by surfaces is considered to be dependent on is based on real-time measurements, supported by simulta- physico-chemical and biological surface properties (Hicks et neous measurements of meteorological quantities covering al., 1987), but also on the solubility and reactivity of the the late dry (biomass burning) season, the transition period, gaseous compound (Wesely, 1989). Hence, soil character- and the onset of the wet season (clean conditions). Fluxes istics, plant stomatal activity and trace gas chemical prop- of NO2, HNO3, NH3 and HONO are estimated by inferen- erties largely determine the deposition velocity. The atmo- tial methods. N wet deposition was determined by collection spheric dry removal of aerosol particles, which may contain of rainwater and subsequent analyses. We estimate the total + − N species such as ammonium (NH4 ) and nitrate (NO3 ), is (wet + dry) annual N deposition at this pasture site and the a function of the particle size (Nicholson, 1988) but also relative contribution of the individual N species. depends on the particle density (e.g., Seinfeld and Pandis, 1998). Dry deposition is enhanced for large particles (espe- cially those larger than a few micrometers) due to the addi- 2 Experimental tional influence of gravitational settling. N wet deposition is a result of in-cloud scavenging (“rain- 2.1 Field site out”) and below-cloud scavenging (“washout”) of atmo- spheric N constituents (Meixner, 1994). The total (wet + Measurements were performed during 12–23 September dry) N deposition ranges from −1 to −2 kgN ha−1 yr−1 for 2002 (dry season, biomass burning), 7–31 October 2002
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 449
(transition period) and 1–14 November 2002 (wet season, clean conditions) as part of the LBA-SMOCC (Large-Scale Biosphere-Atmosphere Experiment in Amazonia – Smoke, Aerosols, Clouds, Rainfall, and Climate) 2002 campaign (Andreae et al., 2004) at a pasture site in the state of Rondonia,ˆ Brazil (Fazenda Nossa Senhora Aparecida, FNS, 10◦04.700 S, 61◦56.020 W, 145 m a.s.l.). The site is located in the south-western part of the Amazon Basin. The location and a simplified sketch of the measurement site are shown in Fig. 1. The old growth rain forest at FNS was cleared by slash and burn activities in 1977. The vegetation at FNS is dominated by C4 grass species (Brachiaria brizantha) with small patches of Brachiaria humidicola and very few iso- lated palms and bushes, and the site is used as a cattle ranch (∼200 “Blanco” cattle, Bos indicus hybrid). The stocking rate at FNS was about one animal per hectare when field measurements took place. The pasture does not receive any fertilizer and is not harvested. FNS is located within a strip of cleared land about 4 km wide and several tens of kilome- ters long (Culf et al., 1996). The towns Ouro Preto do Oeste (∼40 800 inhabitants) and Ji-Parana´ (∼110 000 inhabitants) are situated approximately 8 km and 40 km to the ENE and ESE of the site, respectively. The instrumentation for trace gas/aerosol sampling and online analyses was arranged in an air conditioned wooden house. Rain samples were collected nearby the house (see Fig. 1). An automatic weather station (Met 1) was located in a distance of ∼20 m to the S, and a meteorological tower (Met 2) was situated ∼200 m to the NE of the inlets for trace gas and aerosol measurements. While the sampling site provides a sufficient uniform fetch expanding for 1–2 km from the sampling location in each direction (Andreae et al., 2002), local flow distortions may be caused by the wooden house and some instrument shelters. A more detailed de- Fig. 1. Location of the LBA-SMOCC measurement site Fazenda scription of the measurement site is given in Andreae et Nossa Senhora Aparecida (FNS) in Rondonia,ˆ Brazil. al. (2002) and Kirkman et al. (2002).
2.2 Sampling and analysis 20 min (dry season), 40 min (transition period) and 60 min Table 1 summarizes the specifications of the instruments for (wet season) (Trebs et al., 2004). Aerosol samples of either the measurement of trace gases, aerosol species and meteo- PM2.5 (Dp≤2.5 µm) or total suspended particulate matter rological quantities. Water-soluble N containing trace gases (TSP) were collected. A detailed description and verifica- + tion of the measurement method and of the inlet system can (NH3, HNO3 and HONO) and related aerosol species (NH4 − be found in Slanina et al. (2001) and Trebs et al. (2004). and NO3 ) were measured on-line. Air was taken from a height of 5.3 m above ground through a sophisticated inlet The chemiluminescence NO/NOx analyzer (Thermo En- system, which was designed to reduce wall losses of soluble vironment Instruments) (see Table 1) was equipped with a gases (especially HNO3) and to minimize aerosol losses due molybdenum converter to transform ambient NO2 to NO. to non-isokinetic sampling (see Trebs et al., 2004). Soluble However, the converter basically responds to the sum of NO2 − gases were scavenged with a wet-annular denuder (WAD) + HNO3 + HONO + PAN + aerosol NO3 + organic nitrates. (Wyers et al., 1993), which was combined with a Steam-Jet- Therefore, it is likely that NO2 measurements might be bi- Aerosol Collector (SJAC) (Khlystov et al., 1995) to collect ased by a positive artifact (Fehsenfeld et al., 1990). During particulate N species. For both gaseous and aerosol com- our study, the inlet line for NO/NOx measurements had a pounds, sample collection was followed by subsequent on- length of 25 m (inner diameter = 4.4 mm). Therefore, highly line analysis (ion chromatography (IC)) for anions and flow- soluble and sticky species such as HNO3 and HONO are + injection analyses (FIA) for NH4 . Cycle times were set to assumed to be at least partly removed within the long inlet www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 450 I. Trebs et al.: N deposition in the tropics
Table 1. Specifications of the instrumentation for the measurement of trace gases, aerosol species and meteorological quantities at FNS during LBA-SMOCC 2002 (all heights are above ground).
Parameter/position/ Time Technique or sensor Model, manufacturer Detection limit/ preci- height resolution sion (1) NH3, HNO3, HONO 20–60 min Wet-annular denuder ECN, Petten, Netherlands ≤0.015 ppb (3σ ) for (5.3 m) (WAD), IC, FIA(2) acids, ≤0.118 ppb (3σ ) for NH3 + − ≤ Aerosol (NH4 , NO3 ) 20–60 min Steam-Jet Aerosol Collector ECN, Petten, Netherlands 0.015 ppb (3σ ) for (5.3 m) (SJAC), IC, FIA anions, ≤0.118 ppb + (3σ ) for aerosol NH4 NO concentration 5 min Gas-phase chemilumi- Model 42C TL (trace 0.05±0.025 ppb (10 m) nescense level),Thermo Environment Instruments Inc., USA NO2 concentration 5 min Catalytic conversion of NO2 Model 42C TL (trace level), 0.05±0.025 ppb (10 m) to NO by molybdenum con- Thermo Environment Instru- verter (at 325◦C), gas phase ments Inc., USA chemiluminescence O3 concentration 5 min UV absorption Model 49C 1±0.5 ppb (10 m) Thermo Environment Instruments Inc., USA Air temperature 1 min Pt-100 resistance sensor MP-103A-CG030-W4W ±0.1 K (Met 1, 0.5 m and 5 m) Rotronic, Switzerland Relative humidity 1 min Capacitive sensor MP-103A-CG030-W4W ±1.5% (Met 1, 0.5 m and 5 m) Rotronic, Switzerland Surface wetness 1 min Surface wetness grids at soil 237 WSG, Campbell Scien- – (Met 1) surface tific Ltd., UK Global radiation flux 1 min Pyranometer sensor LI200SZ <±3% (Met 1, 5 m) (LI-COR, Lincoln, Nebraska, USA) Eddy covariance; three di- 10 Hz 3-D ultrasonic anemometer Solent 1012R2, ±5–20% mensional wind and temper- Gill Instruments, UK ature fluctuations (Met 2, 4 m) H2O mixing ratio 10 Hz Infrared closed-path absorp- IRGA LI-COR 6262 ±1% (Met 2, 4 m) tion (LI-COR, Lincoln, Nebraska, USA) Shortwave radiation in and 1 min Pyranometer sensor Kipp & Zonen Pyranometer ±2% out (Met 2, 8.5 m) CM 21 Longwave radiation in and 1 min Pyrgeometer sensor Kipp & Zonen Pyrgeometer ±3% out (Met 2, 8.5 m) CG
(1) 20 min: dry season (12–13 September), 40 min: transition period (7–31 October), 60 min: wet season (1–14 November) (2) IC: ion chromatography, FIA: flow injection analysis
tubing. In addition, mixing ratios of HNO3 and HONO were namely air temperature (T), relative humidity (RH), surface usually below 0.5 ppb (Trebs et al., 2004), indicating that wetness, global radiation flux, momentum, latent and sensi- interferences would be marginal in case any of these gases ble heat flux and ingoing/outgoing short- and longwave ra- would reach the chemiluminescence analyzer. However, el- diation. Eddy covariance measurements were conducted us- evated emissions of isoprene in the area may result in con- ing a Gill 3D-sonic anemometer. The H2O mixing ratio was siderable concentrations of PAN such that we may not rule monitored by a fast LI-COR infrared gas analyzer, and its − out interferences due to this compound. Aerosol NO3 was analog output was directly fed to the ultrasonic anemome- eliminated by the application of an inlet filter. ter A/D converter. Post-processing of the eddy covariance Also listed in Table 1 are those meteorological sensors that data (EDDYWSC, software by Alterra, Wageningen Uni- were used to measure the quantities involved in this study, versity Research, Netherlands) resulted in 30 min averages
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 451 of sensible heat flux, latent heat flux, friction velocity and form (see Thom, 1975). The roughness length z0 of the grass Monin-Obukov length. More details on the eddy covariance surface at the FNS site was taken as 0.11 m (cf. Kirkman et measurements and corresponding data evaluation/calibration al., 2002) and zref was 5.3 m and 10 m for the WAD/SJAC procedures are given in Araujo et al. (2002). and for the NOx measurements, respectively (Table 1). To Precipitation was sampled from 12 September to 14 account for conditions when the reliability of micrometeoro- November 2002 using a wet-only rainwater collector (Ae- logical techniques was low, data were rejected for u∗≤0.01 m −1 rochem Metrics). A total of 23 rainstorm events were col- s and zref/L≥5, i.e. when extremely low turbulence and/or lected representing ∼100% of the precipitation in this period. very high thermal stability was prevailing. Also, data were ◦ Rain samples were stored in the dark at 4 C using polyethy- rejected for zref/L≤−5, which reflects cases of very high lene bottles which were previously cleaned with deionized thermal turbulence production (when Monin-Obukov simi- water and preserved with Thymol. In order to trace possible larity is no longer valid (Ammann, 1999). Thus, about 10% contaminations, the sample pH was measured directly after of the dataset were not used for the flux calculations. + − sampling and before analysis. Analyses of NH4 , NO3 and Rb determines the exchange of gaseous matter by − NO2 were performed for all samples using a Dionex DX600 molecular-turbulent diffusion across the viscous laminar sub- ion chromatograph at the Laboratorio´ de Ecologia Isotopica,´ layer immediately above the vegetation elements and can be CENA/USP (Sao˜ Paulo, Brazil). The detection limit was described by (Hicks et al., 1987):
0.05 µM for all species. More details about sampling and 2 analysis procedures are provided by Lara et al. (2001). 2 Sc 3 Rb = (3) Moreover, a twin Differential Mobility Particle Sizer κ × u∗ P r (DMPS) was employed to measure the dry aerosol parti- cle size distribution in the diameter range from 3 to 850 nm where Sc and Pr are the Schmidt and Prandtl number, respec- (cf. Rissler et al., 2004). The size distribution of particles tively. Pr is 0.72 and Sc is a strong function of the molecular with aerodynamic diameters from 1 to 4 µm was measured diffusivity of the trace gas. Values for Sc were taken from with an Aerodynamic Particle Sizer (TSI APS 3310). Hicks et al. (1987) and Erisman et al. (1994) for the differ- ent trace gas species. The surface resistances Rc could not be directly determined from our field measurements; hence 3 Theory: Estimation of N dry and wet deposition values were adopted from the literature (see Sect. 4.4). The inferential method is valid for trace gases whose mix- 3.1 Trace gas fluxes ing ratio just above the soil and/or vegetation elements is zero. The observation of a net NO2 deposition flux to the Dry deposition fluxes of trace gases have been estimated FNS pasture by Kirkman et al. (2002) justifies the applica- using the inferential method, which is based on the “big tion of the inferential model for NO2 in our study. This is also leaf multiple resistance approach” (Wesely and Hicks, 1977; valid for HNO3, which typically features a rapid downward Hicks et al., 1987). The deposition flux (F) (µg m−2 s−1) of flux with negligible Rc and corresponding high Vd (Hanson a nonreactive trace gas for which the surface is a sink under and Lindberg, 1991). all ambient conditions is defined by: By contrast, NO, HONO and NH3 may be both deposited X(z ) to and emitted from surfaces. Formally, this can be accounted F = −V · X(z ) = − ref (1) d ref R + R + R for by a so-called canopy compensation point concentra- a b c −3 tion Xc(µg m ) that generally refers to the concentration −3 where X(zref) is the trace gas concentration (µg m ) at the of the compound just above the soil and/or vegetation ele- reference height zref (m) and Vd denotes the dry deposition ments (Nemitz et al., 2004a). Xc represents a concentration −1 velocity (m s ), which is the reciprocal of the sum of the analogue of Rc and is the air concentration at which com- −1 turbulent resistance (Ra) (s m ), the quasi-laminar or vis- peting chemical and biological consumption and production −1 cous boundary layer resistance (Rb) (s m ), and the surface processes balance each other (i.e., the net flux is zero) (see −1 resistance (Rc) (s m ). According to Hicks et al. (1987) Ra Sutton et al., 1995): between the reference height (zref) and the roughness length Xc − X(zref) z0 (m) is given by: F = (4) Ra + Rb 1 z = ref − Ra ln 9H zref (2) The net NO emission from the FNS pasture site determined κ × u∗ z 0 L by Kirkman et al. (2002) was very low (0.65 ngN m−2 s−1 where κ denotes the von Karman constant (0.41) and L is the or 0.17 kgN ha−1 yr−1), thus we neglected any contribution Monin-Obukov length (m), a measure of atmospheric stabil- of NO to the surface-atmosphere exchange of N species in ity that is derived from the sensible heat flux and the fric- our study. HONO is generally assumed to be formed by tion velocity u∗ (Garratt, 1992). 9H (zref/L) is the stability heterogeneous reaction of NO2 with surface water (Harri- correction function for heat and inert tracers in its integral son et al., 1996) and it may subsequently be emitted from www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 452 I. Trebs et al.: N deposition in the tropics plant foliar cuticles or soil surfaces. Since there is no indi- Sect. 4.4). Furthermore, Xd (NH3) at time step t is a function cation for any direct HONO emissions by plants (Schimang of the adsorption charge and of the capacitance of the epicu- et al., 2006), the HONO compensation point concentration ticular water film (for details see Sutton et al., 1998). The Xc(HONO) is expected to be a function of the NO2 mixing values for the adsorption charge were adopted from Sutton et ratio (see Sect. 4.4). al. (1998) and the capacitance is a function of the epicuticular To predict the bi-directional surface-atmosphere exchange pH (see Sect. 4.4) of NH3 at the FNS site, we applied a dynamic resistance model proposed by Sutton et al. (1998). Besides uptake and 3.2 Aerosol fluxes emission of NH3 via plant stomata, the dynamic model ac- counts for absorption of NH3 by epicuticular water films un- Up to date, no well established bulk resistance models exist der very humid conditions, and subsequent re-evaporation for the dry deposition of particles. Significant discrepancies (capacitive leaf surface exchange). Since the FNS site is used have been observed between experimental results and model as a cattle ranch and the NH3 flux directly from the soil is predictions (Ruijgrok et al., 1995). The theoretical frame- assumed to be negligible compared to that originating from work proposed by Slinn (1982) is widely used in modeling cattle excreta, we considered a direct NH3 flux from cattle studies to predict particle deposition velocities. However, manure and urine F(NH3)e. The net NH3 flux F(NH3) can Wesely et al. (1985) derived an empirical parameterization −1 be related directly to the NH3 canopy compensation point for the dry deposition velocity Vp (m s ) of submicron sul- concentration Xc(NH3) (Sutton et al., 1998) and is composed fate aerosols (Dp=0.1–1.0 µm) to grass surfaces: of its component fluxes through plant stomata, Fs(NH3), the flux in or out of the epicuticular water film (adsorption ca- Vp = u∗ · 0.002, for L ≥ 0 (8a) pacitor), Fd (NH3), and F(NH3)e: Xs(NH3) − Xc(NH3) " 2 # F(NH3)t = −a 3 Rs(NH3) Vp = u∗ · 0.002 · 1 + , for L < 0 (8b) | {z } L F s(NH3) Xd (NH3)t − Xc(NH3) + +F(NH3)e (5) where a=300 m. This approach generally results in much Rd (NH3) | {z } higher Vp values for submicron particles than predicted by F d(NH3) the Slinn model. Several other studies (e.g., Garland, 2001; X (NH ) − X(NH , z ) = c 3 3 ref Nemitz et al., 2004b; Vong et al., 2004) also showed that Ra + Rb Vp to surfaces of low aerodynamic roughness may be much where Xs(NH3) denotes the NH3 stomatal compensation larger than predicted by the Slinn model, even for particles −3 point concentration (µg m ) and Rs(NH3) denotes the NH3 of different chemical composition than studied by Wesely et −1 stomatal resistance (s m ). Xd (NH3)t is the NH3 adsorp- al. (1985). In this study either PM2.5 or TSP was sampled. −3 + tion concentration (µg m ) associated with the “leaf sur- NH4 is known to be largely attributed to fine mode aerosols face capacitor” at time step t and Rd (NH3) is the charging with Dp≤1.0 µm (Seinfeld and Pandis, 1998), which was resistance of the capacitor (s m−1) (see Sutton et al., 1998). also observed during the SMOCC measurement campaign 1 − Xc(NH3) is then determined by Sutton et al. (1998): (Fuzzi et al., 2005 ). By contrast, aerosol NO3 exhibited a bimodal size distribution (Falkovich et al., 2005). Since we X ( ) = c NH3 t did not find an empirical relationship in the literature to esti- X(NH3, zref)/(Ra + Rb) + Xs(NH3)/Rs(NH3) + Xd (NH3)t /Rd (NH3) + F(NH3)e mate Vp for aerosols with Dp≥1.0 µm and the contribution −1 −1 −1 − (Ra + Rb) + Rs(NH3) + Rd (NH3) of coarse mode aerosol NO3 relative to that of fine mode (6) + aerosol NH4 is presumably small (cf. Sect. 4.2, Table 2), the parameterization by Wesely et al. (1985) is considered as a Xs(NH3) can be parameterized according to (Farquhar et al., 1980; Sutton et al., 1994): reasonable approximation for both PM2.5 and TSP samples.
161 512 1 (−4507.11/TS ) Fuzzi, S., Decesari, S., Facchini, M. C., Cavalli, F., Emblico, Xs(NH3) = · 10 · 0 · 17 000 (7) TS L., Mircea, M., Andreae, M. O. Trebs, I., Hoffer, A., Guyon, P., where T is the surface temperature (K) which was de- Artaxo, P., Rizzo, L. V., Lara, L. L., Pauliquevis, T., Maenhaut, S W., Raes, N., Chi, X., Mayol-Bracero, O. L., Soto, L., Claeys, rived from the outgoing longwave radiation by applying M., Kourtchev, I., Rissler, J., Swietlicki, E., Tagliavini, E., Schkol- the Stefan-Boltzmann law. 0 is the ratio of apoplastic + + nik, G., Falkovich, A. H., Rudich, Y., Fisch, G., and Gatti, L. [NH4 ]/[H ] and was adopted from the literature (Sect. 4.4). V. : Overview of the inorganic and organic composition of size- Rs(NH3) can be calculated from the measured latent heat segregated aerosol in Rondonia,ˆ Brazil, from the biomass burning flux for relatively dry daytime conditions and in the ab- period to the onset of the wet season, J. Geophys. Res., submitted, sence of precipitation according to Nemitz et al. (2004a) (see 2005.
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 453
Table 2. Summary of trace gas and aerosol mixing ratios (∗) during the dry season (12–13 September), the transition period (7–31 October) and the wet season (1–14 November) at FNS during LBA-SMOCC 2002 (conversion factors from ppb to µg m−3 for standard conditions of + − 298.15 K and 1000 hPa: NO: 1.21, NO2: 1.86, O3: 1.94, NH3: 0.69, HNO3: 2.54, HONO: 1.89, aerosol NH4 : 0.73, aerosol NO3 : 2.77).
Dry season Transition period Wet season Species m P 0.25 P 0.75 n m P 0.25 P 0.75 n m P 0.25 P 0.75 n (ppb) (ppb) (ppb) (X) (ppb) (ppb) (ppb) (X) (ppb) (ppb) (ppb) (X) NO 0.09 0.07 0.15 372 0.09 0.07 0.17 295 0.12 0.07 0.13 229 NO2 4.54 2.85 6.46 630 1.78 1.23 2.56 697 1.07 0.82 1.76 268 O3 24.3 11.76 34.87 826 25.59 15.91 32.95 781 14.74 10.77 19.39 316 NH3 1.81 1.10 2.91 298 1.06 0.5 1.74 236 0.55 0.38 0.85 60 HNO3 0.16 0.10 0.25 317 0.06 0.03 0.13 210 0.06 0.04 0.08 52 HONO 0.12 0.08 0.27 323 0.07 0.05 0.09 315 0.06 0.04 0.07 139 + Aerosol NH 1.01 0.73 1.51 291 0.54 0.33 0.85 267 0.47 0.32 0.62 66 −4 Aerosol NO3 0.34 0.17 0.61 297 0.09 0.06 0.15 282 0.06 0.04 0.07 33
(∗) m: median, P0.25: 0.25 percentile, P0.75: 0.75 percentile, n: number of determined data points above the limit of detection (for aerosol species PM2.5 and bulk measurements were included). NOx/O3 data were synchronized to the WAD/SJAC data.
3.3 Determination of characteristic time scales Kirkman et al. (2002). Thereby, j(NO2) was estimated from global radiation data using a relation derived from simultane- The resistance-based approaches presented above to calcu- ous measurements of global radiation and j(NO2) in Amazo- late surface-atmosphere exchange fluxes rely on the “con- nia during LBA-EUSTACH (cf. Kirkman et al., 2002). stant flux layer assumption”, which implies that the trace compounds are considered as chemically-non-reactive trac- HONO is rapidly photolyzed during daylight hours. The ers, such that their flux within the atmospheric surface layer chemical time scale for daytime HONO photolysis is given −1 is constant. However, sufficiently accurate fluxes of com- by τ(HONO)photol.=j(HONO) , whereby the parameteriza- pounds that undergo rapid chemical transformation can be tion provided by Kraus and Hofzumahaus (1998) was used to estimated as long as characteristic chemical time scales are relate j(HONO) to j(NO2). The chemical time scale of het- one order of magnitude larger than turbulent transport times erogeneous HONO formation at nighttime τ(HONO)het. was (Damkohler¨ ratio Dr <0.1) (De Arellano and Duynkerke, derived by considering the HONO production rate PHONO 1992). Following De Arellano and Duynkerke (1992) the (ppb h−1) determined directly from our measurements (see characteristic time of turbulent transport τturb (s) can be cal- Sect. 4.3). Homogeneous daytime HONO formation may culated as: proceed via reaction of NO with OH radicals; however, this u∗ process is very slow and can be neglected compared to day- τ = κ · (z + z ) · (9) turb ref 0 2 time HONO photolysis (Lammel and Cape, 1996). σw where σw denotes the standard deviation of the vertical wind Photochemical reactions involving NH3 and HNO3 are not − component (m s 1). Atmospheric stability strongly deter- considered to contribute to any flux divergence since their mines the time scale of turbulent transport, which typically chemical time scale is much larger than the time scale of ranges from a couple of seconds under unstable conditions turbulent transport (Pandis et al., 1995). By contrast, time up to 2.5 h under stable conditions (for a layer of 10 m height) scales to achieve gas/aerosol equilibrium between gaseous (Dlugi, 1993). NH3, HNO3 and particulate NH4NO3 may occur within the For the characteristic chemical time scale of the NO- time frame of turbulent transport (few seconds for submi- NO2-O3 triad, its photo-stationary equilibrium has to cron particles) and have to be explicitly considered here be considered. NO2 is rapidly photolyzed to NO in (Dlugi, 1993; Meng and Seinfeld, 1996). As reported by the troposphere and it is re-formed by reaction of NO Trebs et al. (2005), the measured concentration product of with O3. The overall chemical time scale of the NO- NH3×HNO3 was persistently below the theoretical equilib- NO2-O3 triad τ(NO-NO2-O3) is given by the combina- rium dissociation constant of the pure NH3/HNO3/NH4NO3 × ] −1 = −1 −1 tion of τ(NO)=(k2 [O3 ) , τ(NO2) k1 =j(NO2) and system during daytime (RH<90%), but approached the the- −1 τ(O3)=(k2×[NO]) (Lenschow, 1982) (where k is the re- oretical equilibrium value during nighttime (RH>90%). The action rate constant) and was calculated in accordance to desired equilibration time scale τ(NH3-HNO3-NH4NO3) www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 454 I. Trebs et al.: N deposition in the tropics can be estimated according to (Wexler and Seinfeld, 1992): derived in a similar way; for instance, the VWM concen- tration obtained during September was considered as repre- −1 τ(NH3 − HNO3 − NH4NO3) sentative of the entire dry season and then multiplied by the ∞ historical mean rainfall for the dry season. The same proce- Z m(Rp)dRp = 3D × (10) dure was performed for October and November representing λ 2 1 + × R × ρp transition period and wet season, respectively. 0 α+Rp p whereby the aerosol mass size distribution m(Rp)dRp was 4 Results and Discussion related to the measured aerosol number size distribution. Rp is the particle radius (m), D is the geometric mean of the In the following, the general patterns of our results will be diffusivity of semi-volatile gaseous species (m2 s−1), m is −3 shown as diel courses of the medians of measured and in- the water-soluble particle mass (kg m ), λ is the mean ferred quantities. Diel courses will be presented either for × −8 free path of air (6.51 10 m at 293.15 K), α denotes the part of the dry season (biomass burning, 12–23 Septem- accommodation coefficient (0.001<α<1) (Wexler and Se- ber) or, where the investigated quantities are relatively in- infeld, 1990), and ρp is the particle density. A value of dependent of season, for the entire measurement campaign ± −3 ρp=1.35 0.15 g cm was used as determined by Reid and (12 September–14 November 2002). To derive surface- Hobbs (1998) for Amazonian biomass burning aerosols. atmosphere exchange fluxes of N-containing compounds, The comparison of τ with τ(NO-NO -O ), turb 2 3 meteorological data and NOx/O3 data were synchronized to τ(HONO)photol., τ(HONO)het. and τ(NH3-HNO3-NH4NO3) the time resolution of the WAD/SJAC system. The conven- provides a test whether reactive species can be treated as tion of negative downward fluxes (net deposition) and posi- a passive tracer during their vertical transport within the tive upward fluxes (net emission) has been adopted. surface layer (Dr <0.1) and, consequently, if the application of the “big leaf multiple resistance approach” and the 4.1 Meteorological conditions “canopy compensation point model” are justifiable. In the Amazon Basin, nighttime radiative cooling usually 3.4 N wet deposition results in the formation of a shallow, decoupled nocturnal boundary layer of high thermodynamic stability character- Storm size influences the chemical composition of the rain- ized by very low wind speeds; while the development of a water, in the sense that larger storms tend to be more dilute. deepening, convectively mixed layer starts with the heating This dependence demands the use of volume weighted means of the surface in the morning (Nobre et al., 1996; Fisch et al., (VWM) for the calculation of monthly and annual concentra- 2004). Median nighttime wind speeds (at height z=5.0 m) tion averages. For each rainwater solute, the concentrations at FNS were low (≤1 m s−1). By contrast, wind speeds measured in the sample were combined to create a VWM were substantially higher during daytime, ranging from 2 −1 concentration for each collection date: to 3 m s . The 24-h average friction velocity (u∗) was ∼0.16 m s−1 and reached maxima of ∼0.4 m s−1 around Pn i=1 cai × vi noon. Differences of u∗ between the dry season, the tran- VWMa = P (11) vi sition period, and the wet season were marginal. As a con- sequence of high global radiation fluxes (500–900 W m−2, where cai is the concentration of species a in sample i, n median) during daytime, ambient temperatures at the FNS is the number of samples and vi is the volume of precipi- site (z=5.0 m) ranged between 30 and 36◦C during sunlight tation solution for sample i. The rain water concentrations hours, while nighttime temperatures were much lower (20– measured during SMOCC were extrapolated to the whole 25◦C). High relative humidities (RHs)(z=5.0 m) were ob- year. Since storm size can vary from year to year we used served during nighttime (90–100%), while daytime RH val- long-term rainfall amounts to obtain the annual N deposition ues usually dropped to 40–50% during afternoon hours. Lo- rate (cf. Lara et al., 2001). Rainfall amounts and air mass cal meteorology changed only marginally from the dry sea- trajectories during the SMOCC period were comparable to son to the wet season. Strong rain events (≥30 mm h−1) oc- conditions during the same period of previous years. Thus, curred in the afternoon hours during the transition period and these months can be considered as representative of the me- the wet season. A more detailed description of diel meteoro- teorological conditions during this period in this area. The N logical conditions is presented by Trebs et al. (2005). wet deposition rates for September, October and November were derived by summation of the VWM concentration de- 4.2 Mixing ratios termined for each precipitation event throughout the month, and subsequent multiplication by the mean of rainfall in the Median diel variations of mixing ratios of NO, NO2,O3, + respective month obtained from historical time series (see NH3, HNO3, HONO and the inorganic aerosol species NH4 − http://www.aneel.gov.br). The annual N wet deposition was and NO3 are shown in Fig. 2a–h for 12 to 23 September (dry
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 455
NH3 (median) P 0.25 - P 0.75 range HNO3 (median) P 0.25 - P 0.75 range 8 0.5 (a) (b) 7 0.4 6 5 0.3 , ppb 3 , ppb 4 3 3 0.2 NH 2 HNO 0.1 1 0 0.0 HONO (median) P 0.25 - P 0.75 range NO (median) P 0.25 - P 0.75 range 0.6 1.25 (c) (d) 0.5 1.00 0.4 0.75 0.3 0.50
0.2 NO, ppb HONO, ppb 0.1 0.25
0.0 0.00
NO2 (median) P 0.25 - P 0.75 range O3 (median) P 0.25 - P 0.75 range 12 60 (e) (f) 10 50
8 40
, ppb 6 30 , ppb 2 3 O
NO 4 20
2 10
0 0 + - aerosol NH4 (median) P 0.25 - P 0.75 range aerosol NO3 (median) P 0.25 - P 0.75 range 2.5 1.4 (g) (h) 1.2 2.0 1.0 , ppb , ppb - + 3 4 1.5 0.8
1.0 0.6 0.4 0.5 Aerosol NO
Aerosol NH 0.2 0.0 0.0 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 time of day [local] time of day [local]
+ − Fig. 2. Diel variations of (a) NH3, (b) HNO3, (c) HONO, (d) NO, (e) NO2, (f) O3, (g) aerosol NH4 (PM2.5) and (h) aerosol NO3 (PM2.5) measured during 12–23 September 2002 (dry season, biomass burning) at FNS during LBA-SMOCC 2002 (conversion factors from ppb −3 to µg m for standard conditions of 298.15 K and 1000 hPa: NH3: 0.69, HNO3: 2.54, HONO: 1.89, NO: 1.21, NO2: 1.86, O3: 1.94, + − aerosol NH4 : 0.73, aerosol NO3 : 2.77). Symbols and grey shading represent medians and interquartile ranges (0.25 to 0.75 percentiles), respectively. season, biomass burning). Table 2 summarizes ambient mix- oxidize significant amounts of NO (see Fig. 2e, f), but may ing ratios measured during the dry season (12–23 Septem- also be partially attributed to rapid HONO photolysis. NO2 ber), the transition period (7–31 October) and the wet season was the most abundant N-containing trace gas during all (1–14 November). A detailed discussion of seasonal and diel three seasons and reached an average mixing ratio of ∼5 ppb + cycles observed for NH3, HNO3, HONO and aerosol NH4 during the dry season. NO2 featured a pronounced diel cy- − cle with nighttime mixing ratios two times higher than dur- and NO3 is given in Trebs et al. (2004, 2005). Despite intensive biomass burning activity during the dry ing daytime (Fig. 2e). Apparently, NO2 was accumulated in season, NO mixing ratios were very low (Fig. 2d). The sharp a shallow nocturnal boundary layer of high thermodynamic peak between 06:00 and 09:00 LT was most likely due to stability due to: (i.) the low water-solubility of NO2 and consequently its low affinity to be taken up by epicuticular rapid photolysis of accumulated nighttime NO2 shortly af- water films and (ii.) chemical production through reaction of ter sunrise, when O3 mixing ratios were still too low to re- www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 456 I. Trebs et al.: N deposition in the tropics
1000000 a)
), s τ (Z = 5.3 m) τ 100000 turb ref
τ 10000 turb (Zref = 10 m)
τ (HONO)het. 1000 τ (HONO)photol. 100 τ (NO-NO2-O3) 10 Characteristic time scale (
1 10000 b) ), s
τ
τ 1000 (NH3-HNO3-NH4NO3), transition period
τ 100 (NH3-HNO3-NH4NO3), dry season
τ (NH3-HNO3-NH4NO3), 10 wet season
Characteristic time scale ( 1 10000 c) ), s
τ
τ 1000 (NH3-HNO3-NH4NO3), transition period
τ 100 (NH3-HNO3-NH4NO3), dry season
τ (NH3-HNO3-NH4NO3), 10 wet season
Characteristic time scale ( 1 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 time of day [local]
Fig. 3. Diel variation of characteristic turbulent time scale τturb for the reference heights zref=10 m and zref=5.3 m in comparison to (a) chemical time scale of the NO-NO2-O3 triad (τ(NO-NO2-O3)), daytime HONO photolysis (τ(HONO)photol.) and heterogeneous HONO formation at night (τ(HONO)het., dry season only), (b) upper estimate of equilibration time scales for the NH3/HNO3/NH4NO3 system τ(NH3-HNO3-NH4NO3) (α=0.1, PMinorg.≤20%) and (c) lower estimate of τ(NH3-HNO3-NH4NO3) (α=1, PM=100%) at FNS during LBA-SMOCC 2002. Except for τ(HONO)het., data from all seasons were used.
NO with O3 in the absence of NO2 photolysis. Through the 4.3 Characteristic timescales transition period until the wet season, NO2 mixing ratios de- clined by a factor of four. O3mixing ratios exhibited a typical Characteristic turbulent time scales (τturb) have been cal- diel variation, which mirrors that of NO2 (high values during culated according to Eq. (9) for zref=10 m (NOx/O3 mea- daytime and lower values during the night; Fig. 2f). This surements) and for zref=5.3 m (WAD/SJAC measurements) was mainly caused by (i.) photochemical daytime produc- (cf. Table 1). To calculate the characteristic time scale tion, (ii.) convective mixing within the boundary layer and for heterogeneous nighttime HONO buildup (τ(HONO)het.) from the free troposphere during daytime and (iii.) dry de- (Fig. 3a), only dry season nighttime HONO production rates position and reaction with NO in a thermally stable stratified (PHONO) were considered. For the transition period and nocturnal boundary layer. the wet season PHONO could not be determined since the HONO diel variation was substantially reduced (cf. Trebs et
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 457 al., 2004). In contrast to other studies (e.g., Alicke et al., Table 3. Ranges of daytime and nighttime surface resistances Rc 2003; Harrison and Kitto, 1994; Lammel and Cape, 1996), for HNO3 and NO2 used to estimate flux scenarios at FNS during our measurements revealed a relatively small average value LBA-SMOCC 2002. −1 of PHONO=0.04 ppb h . Figure 3a shows that τturb was at least two orders of magnitude smaller than τ(HONO)photol. Scenario Rc(HNO3)Rc(NO2) − − and τ(HONO)het., resulting in Dr 0.1. Considering the [s m 1] [s m 1] chemical time scale of the NO-NO2-O3 triad τ(NO-NO2-O3) Day high flux 1 75 (Fig. 3a), largest D values are found between 17:00 LT and r low flux 50 550 20:00 LT, exceeding a value of 0.2 (zref=10 m). However, Night high flux 1 80 during all other periods Dr for the NO-NO2-O3 triad ranged low flux 15 435 between 0.05 and 0.1. Therefore, we conclude that the ap- plication of Eqs. (1) and (4) to calculate surface-atmosphere exchange fluxes of NO2 and HONO are justified, since chem- ical transformations are too slow to affect the vertical con- SO2−, Na+ and K+ were considered as reference for the stancy of turbulent fluxes. 4 relative change of NO− and NH+ while changing thermo- The estimation of the equilibration time scale τ(NH3- 3 4 dynamic conditions. The mass loss of ammonium nitrate HNO3-NH4NO3) was performed by integrating over the measured particle size distribution according to Eq. (10), was recorded as a function of time for particles with a ra- whereby two different cases were considered. Case 1 is an dius of r<0.8 µm and r>1.6 µm. Timescales τ(NH3-HNO3- upper estimate (Fig. 3b), using an accommodation coefficient NH4NO3) were calculated using the algorithm proposed by α=0.1 (Wexler and Seinfeld, 1992) and taking into account (i.) Kramm and Dlugi (1994) and (ii.) Meng and Seinfeld only the inorganic water-soluble aerosol fraction (≤20% of (1996), resulting in τ(NH3-HNO3-NH4NO3)=100–500 s for fine mode particles and ≥880 s for coarse mode particles. PMtot) (cf. Trebs et al., 2005). Case 2, the lower estimate (Fig. 3c), was calculated using α=1 and assuming that the These values are comparable to characteristic times reported entire aerosol mass (water-soluble and non-soluble species) by Harrison et al. (1990) and Meng and Seinfeld (1996) and are in strong favor of our case 1 (upper estimate, Fig. 3b), are available to equilibrate. As shown in Fig. 3b, c τ(NH3- which implies that the equilibration time scale was always HNO3-NH4NO3) increased substantially from the dry sea- son through the transition period to the wet season. This is much larger than that of turbulent transport. Nevertheless, it obviously caused by much higher particle number concentra- should be noted that the influence of the large soluble organic tions measured during the dry season when biomass burning aerosol fraction typical for the Amazon Basin (cf. Trebs et al., 2005) on gas/aerosol partitioning processes is not exactly took place. For case 1 (upper estimate) (Fig. 3b), Dr for known. Taking into account the potential role of WSOC in the NH3-HNO3-NH4NO3 triad during the dry season ranges from 0.1 to 0.17 at nighttime and is significantly smaller than enhancing aerosol water uptake and subsequently the uptake 0.1 during the day. During the transition period and wet sea- of gaseous species, equilibration time scales may be equal or even faster than turbulent transport (cf. Fig. 3c). son, Dr for the NH3-HNO3-NH4NO3 triad is always signifi- cantly smaller than 0.1 (Fig. 3b). However, for case 2 (lower estimate) (Fig. 3c) nighttime Dr substantially exceeds a crit- 4.4 The inferential approach: selection of input parameters ical value of 0.1 during the dry season and the transition pe- riod. Some surface parameters required for the inferential method To verify the theoretically derived values, results from a to estimate of surface-atmosphere exchange fluxes for trace laboratory study will be discussed briefly. Condensation and gases (Eqs. 1–7) could not be directly derived from the re- evaporation of NH3 and HNO3 to/from particles have been sults of our field measurements. Thus, lower and upper investigated under controlled laboratory conditions. Particles scenarios were estimated (except for HONO; see below), were collected during field campaigns in 1991 (Brunnemann comprising a certain range of surface-atmosphere exchange et al., 1996; Seidl et al., 1996) and 1993/1994 in the eastern fluxes. These scenarios were obtained by varying surface pa- part of Germany (Melpitz). The chemical aerosol compo- rameters over a selected range based on results from studies sition was dominated by (NH4)2SO4, NaCl and soot and is in temperate latitudes. Note that for all parameters and quan- comparable to that described in Brunnemann et al. (1996) tities presented, values indicated as “low” and values indi- and Seidl et al. (1996). More than 90% of the NH4NO3 mass cated as “high” correspond to the estimated lower and upper was found in the accumulation mode (Dp≤1 µm). Aerosol fluxes, respectively. samples were exposed to temperature step changes of 1 to Surface resistances Rc (HNO3) and Rc(NO2): Rc(HNO3) 20 K above 278 K, which simulates a moderate diel cycle. was found to be zero in many studies (e.g., Dollard et al., RH was held constant for these cases. A second set of ex- 1987; Huebert and Robert, 1985). Very recently, Nemitz periments was performed where RH was varied between 30 et al. (2004a) showed that non-zero, however, relatively −1 and 96% at constant temperature. In both cases non-volatile small Rc(HNO3) may exist (Rc=15–95 s m ). The ranges www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 458 I. Trebs et al.: N deposition in the tropics
median P 0.25 - P 0.75 range a) 200 120 a) 180 100
160 -1 140 80 s m ),
120 3
= 5.30 m) 60
ref 100
z evaporation (NH ( s 40 of epicuticular -1 80 R water films 60
, s m 20 a
R 40 0 20
0 100000 140 b)
b) Rd high flux
120 10000 Rd low flux
100 -1 -1 1000 s m ), 3 , s m 80 ) 3 100 (NH
60 d R (HNO b
R 40 10
20 1 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 0 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 time of day [local] time of day [local] Fig. 5. Median diel variations of estimated (a) NH3 stomatal resis- R ) tance Rs(NH3) during daytime, using the measured latent heat flux Fig. 4. Diel variations of (a) turbulent resistance ( a and (b) −1 LE (Rs(NH ) was set to 1000 m s for nighttime periods) and (b) quasi-laminar boundary layer resistance (Rb) exemplary for HNO3 3 at FNS during LBA-SMOCC 2002. Symbols and grey shading rep- lower (high flux) and upper (low flux) NH3 charging resistance of resent medians and interquartile ranges (0.25 to 0.75 percentiles), the surface water layer (Rd (NH3)) at FNS during LBA-SMOCC respectively. Data from all seasons were used. 2002. Data from all seasons were used.
of Rc(HNO3) used to estimate lower and upper surface- Since stomata are thought to be closed during nighttime, −1 atmosphere exchange fluxes are given in Table 3. The sur- Rs(NH3) is set to 1000 s m for nighttime periods. The cal- face uptake of HNO3 is considered to be enhanced by the culation of Rd (NH3) is performed in accordance with Sutton presence of epicuticular water films under the humid condi- et al. (1998). tions at the site (lower Rc values were chosen for nighttime). Compensation point concentration Xc(HONO): Up to Rc(NO2) has been determined previously at the FNS site by date, only two studies provide an estimate for Xc(HONO) Kirkman et al. (2002), hence the maximal and minimal diel (Harrison and Kitto, 1994; Stutz et al., 2002). Here, variations of Rc(NO2) were taken directly from Kirkman et the relationship Xc(HONO)=0.03·X(NO2, zref) was used, al. (2002) to estimate lower and upper surface-atmosphere which was found for grassland in the recent study by exchange, respectively. Corresponding averages of maximal Stutz et al. (2002). This corresponds on average to and minimal Rc(NO2) values are given in Table 3. Xc(HONO)=85 ppt at the FNS site. + + Stomatal resistance Rs(NH3) and charging resistance Apoplastic [NH4 ]/[H ] ratio 0 and epicuticular pH: For Rd (NH3): Rs(NH3) is not known for Brachiaria brizantha 0, we have chosen values of 100 (lower estimate) and 200 grass species, therefore the diffusion of H2O through plant (upper estimate), which range at the lower end of data re- stomata is used for the calculation of upper and lower flux ported for grass in the literature (Loubet et al., 2002; Spindler estimates. Thus, Rs(NH3) is estimated according to Nemitz et al., 2001; van Hove et al., 2002). This may be justified by et al. (2004a) using the measured latent heat flux LE. This the poor soil N status at FNS (Kirkman et al., 2002), because approach is only valid for relatively dry daytime conditions N absorbed by the root medium strongly affects the leaf tis- + (10:00–18:00 LT) and in the absence of precipitation. The sue NH4 concentration (Schjoerring et al., 1998a). The ca- transfer of H2O through plant stomata during daytime rep- pacitance of the epicuticular water film is a function of the resents an upper boundary for trace gas stomatal exchange. pH (see Sect. 3.1, Sutton et al., 1998), that is predominantly
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 459
acidic (Flechard et al., 1999) and is taken as 4.5 (lower esti- Xc _high Xc _low Surface wetness (median) mate) and as 4.0 (upper estimate) (cf. Sutton et al., 1998). 8 100 7 90 NH3 flux from cattle excreta F (NH3)e: In order to esti-
80
-3 6 mate the contribution of cattle excreta to the net NH3 flux, 70
g m 5
µ 60 we considered results of Boddey et al. (2004), who investi- ), gated the cycling of N in Brachiaria pastures in the south 3 4 50 3 40 (NH of the Brazilian province Bahia. According to their results, c 30 X 2 − − Surface wetness, % Bos indicus cattle excreted 37 kgN animal 1 yr 1 in manure 20 1 and 49 kgN animal−1 yr−1 in urine when the pastures were 10 0 0 stocked with two animals per hectare. The stocking rate at 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 FNS was one animal per hectare. About 8% of the excreted time of day [local] N may be released as NH3 (A. F. Bouwman, personal com- munication, 2004). Hence, the average F (NH3)e is estimated Fig. 6. Median diel variation of the simulated lower (low flux) and to be ∼10 ngN m−2 s−1, which is applied for the entire mea- upper (high flux) NH3 canopy compensation point Xc(NH3), and surement period (Eqs. 5, 6). the median diel variation of the measured surface wetness, shown exemplarily for a period during the dry season (12–23 September) 4.5 Resistances, NH canopy compensation point, transfer- at FNS during LBA-SMOCC 2002. For details about low and high 3 estimates, see text. and deposition velocities
The turbulent resistance Ra and the quasi-laminar boundary epicuticular NH+ has evaporated, X (NH ) subsequently de- layer resistance Rb (Fig. 4a, b): These resistances feature a 4 c 3 creases (after 09:00 LT). Lower surface temperatures during typical diel variation with lowest values during daytime (Ra (5.3 m)=20 s m−1, R =30 s m−1), because of strong turbu- the transition period and the wet season have caused lower b −3 estimates of Xc(NH3) (0.3–1 µg m during nighttime and lent mixing within the surface layer, and higher values during − −1 −1 1–4 µg m 3 during daytime). nighttime (Ra (5.3 m)=70 s m and Rb=50 s m ). The stomatal resistance Rs(NH3) and the charging re- Transfer- and deposition velocities Vtr.(NH3), sistance Rd (NH3) (Fig. 5a, b): The median diel cycle of Vtr.(HONO), Vd (HNO3), Vd (NO2) (Fig. 7a–d) and Rs(NH3) estimated from LE during daytime (Fig. 5a) reaches Vp(aerosol) (Fig. 8): For compounds featuring a bi- very low values, falling below 50 s m−1between 10:00 and directional surface exchange the concept of a dry deposition = + + −1 14:00 LT. In the afternoon Rs(NH3) increases and attains val- velocity Vd (Ra Rb Rc) (see Eq. 1) is no longer useful −1 ues >100 s m after 17:00 LT. Rs(NH3) is set to a value of (e.g., Kramm and Dlugi, 1994). Thus, for bi-directional NH3 1000 s m−1 for nighttime periods (see above; not shown). and HONO surface-atmosphere exchange the term transfer − − + The simulated Rd (NH3) (Fig. 5b) is very high during day- velocity Vtr.= sgn(Xc X(zref))/(Ra Rb) (see Eqs. 4, 5) time (>1000 s m−1) and remains below 100 s m−1 at night. will be used henceforth. This is obviously due to the fact that grass leaves became dry The estimated median diel variation of Vtr.(NH3) (Fig. 7a) during sunlight hours (see below, Fig. 6), resulting in a very indicates net deposition (= downward flux) at night- low capacity of the foliar cuticle to adsorb NH3. time (Vtr.(NH3)>0) and net emission (= upward flux) The NH3 canopy compensation point concentration (Vtr.(NH3)<0) during the day. Typically, Vtr.(NH3) is pre- −1 Xc(NH3) (Fig. 6): The Xc(NH3) scenario for the dry sea- dicted to range between −2.6 cm s during the day and −1 son calculated according to Eq. (6) lies well within the range 1.5 cm s at nighttime. For cases of NH3 deposition, of values determined for grassland in other studies (Hester- Vtr.(NH3) is in good agreement with values of 0.1–2 cm berg et al., 1996; Meixner et al., 1996; Spindler et al., 2001; s−1 observed by Erisman and Wyers (1993), Hesterberg Sutton et al., 2001) and is a strong function of surface temper- et al. (1996), Nemitz et al. (2004a), Phillips et al. (2004) ature (cf. Eq. 6, 7). Xc(NH3) is predicted to be particularly and Rattray and Sievering (2001). For cases of NH3 emis- high at daytime, although low values of 0 were used to run sion, Vtr.(NH3) agrees well with values reported by e.g., the model. This is caused by prevailing high surface temper- Nemitz et al. (2004a). The application of the dynamic re- atures at the FNS pasture site (35–40◦C at daytime and 20– sistance model from Sutton et al. (1998), which uses non- ◦ 25 C at night during the dry season). NH3 deposits and dis- zero Xc(NH3) values and takes into account that deposited solves in epicuticular water films at high RHs during night- NH3 may (re-)evaporate from surfaces, represents a rather time. After sunrise the increasing surface temperature causes conservative N deposition estimate (lower boundary). Since the NH3 partial pressure above the epicuticular solution to Xc(NH3) was assumed to be zero in several studies (Gould- increase in accordance to Henry’s law. Hence, Xd (NH3) ing et al., 1998; Russell et al., 2003; Singh et al., 2001; increases, which in turn results in higher Xc(NH3) values. Tarnay et al., 2001), which implies that NH3 is persistently This is visible in Fig. 6 as a distinctive peak of Xc(NH3) net deposited to the surface (see Eq. 1), we have comple- after sunrise between 08:00 and 09:00 LT. Once most of the mented the upper and lower bi-directional flux scenarios by www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 460 I. Trebs et al.: N deposition in the tropics
2.0 2.0 a) b) 1.5 1.5 Vd 1.0 1.0 = 5.30 m) ref 0.5 = 5.30 m)
z 0.5 ( ref
z ( -1 0.0
-1 0.0 -0.5 cm s Vtr -0.5 ), cm s 3 -1.0
-1.0
(NH -1.5 d. -2.0 -1.5 (HONO), tr. and V
-2.5 V -2.0 tr. V -3.0 -2.5
2.5 1.4 c) d) 1.2 2,0 1.0 = 10 m) = 5.30 m) ref ref Vd _high (Rc _low) z z ( ( 1.5 0.8 -1 -1 Vd _high (Rc _low) 0.6 cm s cm s 1.0 ), ), 2 3 0.4 (NO d
(HNO 0.5 V d 0.2 V Vd _low (Rc _high) Vd _low (Rc _high) 0.0 0.0 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 time of day [local] time of day [local]
Fig. 7. Median diel variation of estimated transfer- (Vtr.) and deposition velocities (Vd ) for (a) NH3 (zref=5.3 m) (b) HONO (zref=5.3 m), (c) HNO3 (zref=5.3 m) and (d) NO2 (zref=10 m) at FNS during LBA-SMOCC 2002. The bi-directional NH3 flux was complemented by a “deposition only scenario” (cf. Vd in a). Data from all seasons were used, except for Vtr.(HONO), which only represents results from the dry season, since negligible small upward HONO fluxes were predicted during transition period and wet season.
0.9 The estimated median diel variation of Vtr.(HONO)
(Fig. 7b, dry season only) similarly shows net deposition -1 0.8 at night and net emission during the day, with Vtr.(HONO) 0.7
, cm s −1 −1 p Dp = 0.1 - 1.0 µm ranging from −1.7 cm s to 1.5 cm s . For cases of HONO V 0.6 deposition, these values are comparable to Vtr.(HONO) of 0– − 0.5 1.7 cm s 1 reported by Harrison and Kitto (1994) and Stutz 0.4 et al. (2002).
0.3 The estimates of Vd (HNO3) shown in Fig. 7c, reveal max- −1 0.2 imal median values of 2.3 cm s at around 13:00 LT when −1 0.1 Rc(HNO3)=1 s m is applied (see Table 3, upper flux esti- Particle deposition velocity mate), coinciding with the period of highest turbulence. The 0.0 −1 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 lower estimate, when Rc(HNO3) is taken as 50 s m dur- −1 time of day [local] ing daytime and as 15 s m during nighttime (see Table 3), results in median Vd (HNO3) values nearly equal during day − Fig. 8. Median diel course of the particle deposition velocity Vp and night (≤1.1 cm s 1), which is attributed to the compen- (zref=5.3 m) using the empirical parameterization derived by We- sating effect of lower surface uptake and low turbulent ex- sely et al. (1985) for aerosol particles with D =0.1–1 µm at FNS p change during nighttime. Although Vd (HNO3) was found to during LBA-SMOCC 2002. be higher (3–4 cm s−1) in other studies (e.g., Erisman et al., 1988), our values are still in reasonable agreement with depo- sition velocities determined by Dollard et al. (1987), Muller¨ additionally calculating a “deposition only” scenario accord- et al. (1993), Nemitz et al. (2004a) and Rattray and Siever- ing to Eq. (1) (V in Fig. 7a), where R (NH ) was calcu- d c 3 ing (2001) for vegetative canopies similar to that at FNS. It lated from the parallel resistances R (NH ) and R (NH ) s 3 d 3 should be noted here that for a compound that is exchanged (see Fig. 5a, b). Assuming NH3 to be net deposited to the −1 with a very small Rc (as in the case for HNO3), Vd is mainly pasture site results in Vd =0.5–1.5 cm s (Fig. 7a). dependent on Ra +Rb. Thus, this comparison reflects largely the differences in turbulence and surface roughness.
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 461
Median (dry season_high) Median (dry season_low) Median (dry season_Xc = 0) Median (dry season) 80 0.50 a) b) 0.25 60 -1
-1 0.00 s s -2
-2 40 -0.25 -0.50
20 ngN m ngN m ),
3 -0.75 0 -1.00 F(NH -20 F(HONO), -1.25 -40 -1.50 Median (dry season_high) Median (dry season_low) Median (dry season_high) Median (dry season_low) 0.0 0.0 c) d) -0.5 -5.0 -1 -1.0 -1 s s -2
-1.5 -2 -10.0 -2.0 ngN m ngN m ), -2.5 ), 3
2 -15.0 -3.0 F(NO F(HNO -3.5 -20.0 -4.0 -4.5 -25. 0 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 time of day [local] time of day [local]
Fig. 9. Median diel courses of estimated lower (low) and upper (high) surface-atmosphere exchange fluxes for: (a) NH3 (zref=5.3 m) (b) HONO (zref=5.3 m), (c) HNO3 (zref=5.3 m) and (d) NO2 (zref=10 m). Bi-directional NH3 flux scenarios were complemented by a “deposition only” scenario (Xc(NH3)=0). Only one scenario could be calculated for HONO due to the restricted availability of input parameters (cf. Sect. 4.4). Median diel courses are presented exemplarily for a period during the dry season (12–23 September) at FNS during LBA-SMOCC 2002.
The median diel course of Vd (NO2) shown in Fig. 7d ex- RH usually exceeded 90%. However, a quantification of this hibits maxima of ∼1.3 cm s−1 during daytime, while mini- effect would go beyond the scope of this work. mal values are below 0.2 cm s−1. These values compare rel- atively well with ranges of Vd (NO2) reported by Kirkman et 4.6 Surface-atmosphere exchange fluxes al. (2002). However, the dry deposition of NO2 is thought to be mainly featured by uptake through plant stomata, which In this section, the measured mixing ratios (Sect. 4.2), the would imply that much higher dry deposition velocities (and input parameters for the inferential approach (discussed in hence much lower Rc) occur during daytime. As stated by Sect. 4.4) and the results presented in section 4.5 will be used Kirkman et al. (2002), measured day- and nighttime Rc val- to estimate and discuss surface-atmosphere exchange fluxes ues for NO2 at FNS were similar (see also Table 3), which of N-containing trace gases and aerosol particles. was assumed to be mainly a result of stomatal closure due to Bi-directional fluxes of gaseous NH3, F(NH3) (Fig. 9a): high water vapor pressure deficit at lower RH’s during day- Fluxes found at locations in temperate regions with vegeta- time. tive surfaces similar to the grass surface at FNS are compa- rable to our upper bi-directional estimate of F (NH3) for the Deposition velocities Vp (zref=5.3 m) predicted for par- dry season (biomass burning) (Fig. 9a) (Erisman and Wyers, ticles (Fig. 8) are estimated to be highest during daytime 1993; Flechard et al., 1999; Meixner et al., 1996; Spindler −1 (≤0.8 cm s ) in accordance with higher values of u∗. Dur- et al., 2001). The net emission peak between 08:00 and −1 ing nighttime, Vp generally remains below 0.1 cm s . As 09:30 LT in the upper estimate corresponds to the peak of previously indicated, these values are much larger than Vp the predicted NH3 canopy compensation point (cf. Fig. 6) predicted by Slinn (1982) and are therefore considered as an and the highest observed NH3 mixing ratios (see Fig. 2a). upper estimate. It should be pointed out that aerosol wa- As already mentioned above, this feature may be related to ter uptake at high RHs and the resulting particle growth dur- desorption of NH3 from epicuticular water films due to an ing the deposition process may enhance deposition velocities increase of surface temperature after sunrise (see also Trebs (Khlystov, 1998; Gallagher et al., 1997). In our study, this et al., 2005). Relatively high median NH3 net emission would be particularly relevant for nighttime periods when fluxes (5–70 ngN m−2 s−1, upper estimate) during daytime www.atmos-chem-phys.org/acp/6/447/ Atmos. Chem. Phys., 6, 447–469, 2006 462 I. Trebs et al.: N deposition in the tropics
+ - Median (aerosol NH4 ) Median (aerosol NO3 ) -0.35 a) NO 0.0 2 -0.30 Ammonium -1.0 -1 Nitrate -0.25 Nitrite = 5.30 m) -2.0 ref z
( -3.0
-1 -0.20 s
-2 -4.0 -0.15 -5.0 ngN m -6.0 -0.10 N dry depostion, kgN ha -7.0 -0.05 F (aerosol), -8.0 0:00 3:00 6:00 9:00 12:00 15:00 18:00 21:00 0:00 0.00 September October November time of day [local] -0.35 b) Ammonium Fig. 10. Median diel courses of inferred deposition fluxes of aerosol + − -0.30 Nitrate
NH4 and NO3 shown exemplarily for a period during the dry sea- -1 Nitrite son (12–23 September) at FNS during LBA-SMOCC 2002. -0.25
-0.20 (08:00–17:30 LT) are the result of (i.) relatively high NH3 mixing ratios (Fig. 2a), (ii.) direct NH3 emission from cat- -0.15 tle excreta, (iii.) high surface temperatures and (iv.) cor- responding high values of simulated Xc(NH3) (see Fig. 6). -0.10 N wet depostion, kgN ha The simulated nighttime NH3 net deposition (upper esti- mate) (19:00–06:30 LT) is on average −2 to −13 ngN m−2 -0.05 −1 s . The lower NH3 flux estimate shown in Fig. 9a sug- 0.00 September October November gests that the net exchange of NH3 may be significantly smaller when a higher epicuticular pH (4.5) (Sect. 4.4), a + + Fig. 11. Monthly estimates of (a) N dry deposition of NO2, lower apoplastic [NH4 ]/[H ] ratio 0 (Sect. 4.4) and higher + total ammonium (aerosol NH4 +NH3), total nitrate (aerosol Rd (NH3) (Fig. 5b) are applied. The estimated NH3 net de- − NO +HNO ) and nitrite (HONO) (averages of lower and upper position for this case during nighttime is negligible and the 3 3 flux estimates are shown, except for NH3 which represents an av- predicted daytime NH3 net emission varies between 1 and −2 −1 erage of all three estimates (cf. Fig. 9a) and (b) N wet deposition 28 ngN m s . During the transition period and the wet of ammonium, nitrate and nitrite (error bars represent measurement season, however, the lower NH3 flux estimate (not shown) uncertainties) for September (dry season, biomass burning), Octo- features emission during day and night. The “NH3 deposi- ber (transition period) and November (wet season, clean conditions) tion only” scenario (Xc(NH3)=0; Fig. 9a) exhibits highest at FNS during LBA-SMOCC 2002. fluxes during daytime (08:00–15:00 LT), with values rang- ing between −20 and −40 ngN m−2 s−1. At nighttime, net deposition fluxes between −5 and −20 ngN m−2 s−1 were largest uncertainty of all flux estimates since presently nei- estimated. ther Xc(HONO), nor the exact HONO formation mechanism The predicted bi-directional NH3 flux scenarios are most are well known. sensitive to the pH of the epicuticular water film and the Fluxes of gaseous HNO3, F (HNO3) (Fig. 9c): The esti- + + H and NH4 concentrations in the apoplastic fluid (0). If mated median diel HNO3 flux scenario during the dry sea- the pH of the epicuticular water film were >4.5 and 0 con- son is characterized by highest net deposition values from −2 −1 stant (100–200), the upper NH3 flux estimate would show net 09:00–16:30 LT (∼−0.5 to −4 ngN m s , upper esti- emission also during nighttime. On the other hand, increas- mate), coinciding with highest estimated deposition veloci- ing 0 beyond a value of 200 would result in daytime NH3 net ties (cf. Fig. 7c) and the maximal observed HNO3 mixing emissions significantly higher than observed in urban regions ratios (cf. Fig. 2b). Predicted nighttime HNO3 deposition in Europe or North America. fluxes are lower in accordance with lower turbulence (lower Bi-directional fluxes of gaseous HONO, F (HONO) Vd ) and lower HNO3 mixing ratios. The higher values of (Fig. 9b): The estimated median diel HONO flux during the Rc(HNO3) applied for the lower flux estimate (Table 3) re- dry season (12 to 23 September), shows a small net emission sult in daytime net deposition fluxes that are a factor of two during the afternoon but net deposition of up to −1.3 ngN lower than for the upper flux estimate. −2 −1 m s from 17:00–08:40 LT. We like to note, that the Fluxes of gaseous NO2, F (NO2) (Fig. 9d): The diel course HONO flux scenario presented here probably reflects the of the estimated NO2 flux scenario (dry season) is always
Atmos. Chem. Phys., 6, 447–469, 2006 www.atmos-chem-phys.org/acp/6/447/ I. Trebs et al.: N deposition in the tropics 463
a) dry (biomass burning) season b) transition period
3.8% 6.1% 14.7% dry
21.4% 22.4% 10.9% 1.4% 29.4% dry 1.8%