Atmospheric Environment 43 (2009) 5193–5267
Contents lists available at ScienceDirect
Atmospheric Environment
journal homepage: www.elsevier.com/locate/atmosenv
Review Atmospheric composition change: Ecosystems–Atmosphere interactions
D. Fowler a,*, K. Pilegaard b, M.A. Sutton a, P. Ambus b, M. Raivonen c, J. Duyzer d, D. Simpson e,f, H. Fagerli f, S. Fuzzi g, J.K. Schjoerring h, C. Granier i,j,k, A. Neftel l, I.S.A. Isaksen m,n,P.Lajo,p, M. Maione q, P.S. Monks r, J. Burkhardt s, U. Daemmgen t, J. Neirynck u, E. Personne v, R. Wichink-Kruit w, K. Butterbach-Bahl x, C. Flechard y, J.P. Tuovinen z, M. Coyle a, G. Gerosa aa, B. Loubet v, N. Altimir c, L. Gruenhage ab, C. Ammann l, S. Cieslik ac, E. Paoletti ad, T.N. Mikkelsen b, H. Ro-Poulsen ae, P. Cellier v, J.N. Cape a, L. Horva´th af, F. Loreto ag,U¨ . Niinemets ah, P.I. Palmer ai, J. Rinne aj, P. Misztal a, E. Nemitz a, D. Nilsson ak, S. Pryor al, M.W. Gallagher am, T. Vesala aj, U. Skiba a, N. Bru¨ ggemann x, S. Zechmeister-Boltenstern an, J. Williams ao, C. O’Dowd ap, M.C. Facchini g, G. de Leeuw aq, A. Flossman o, N. Chaumerliac o, J.W. Erisman ar a Centre for Ecology and Hydrology, EH26 0QB Penicuik Midlothian, UK b Risø National Laboratory, Technical University of Denmark, 4000 Roskilde, Denmark c Department of Forest Ecology, University of Helsinki, 00014 Helsinki, Finland d TNO Institute of Environmental Sciences, 3584 CB Utrecht, The Netherlands e Department Radio and Space Science, Chalmers University of Technology, 41296 Gothenburg, Sweden f Norwegian Meteorological Institute, 0313 Oslo, Norway g Istituto di Scienze dell’Atmosfera e del Clima – CNR, 40129 Bologna, Italy h Royal and Veterinary and Agricultural University, 1870 Frederiksberg C, Denmark i UPMC Univ. Paris 06, LATMOS-IPSL; CNRS/INSU, LATMOS-IPSL, 75005 Paris, France j NOAA Earth System Research Laboratory, 80305-3337 Boulder, USA k Cooperative Institute for Research in Environmental Sciences, University of Colorado, 80309-0216 Boulder, USA l Agroscope FAL Reckenholz, Swiss Federal Research Station for Agroecology and Agriculture, 8046 Zurich, Switzerland m Department of Geosciences, University of Oslo, Inst. For Geologibygningen, 0371 OSLO, Norway n Center for International Climate and Environmental Research – Oslo (CICERO), 0349 Oslo, Norway o Laboratoire de Me´te´orologie Physique, Observatoire de Physique du Globe de Clermont-Ferrand, Universite´ Blaise Pascal – CNRS, 63177 Aubie`re, France p Laboratoire de Glaciologie et Ge´ophysique de l’Environnement, Observatoire des Sciences de l’Universite´ de Grenoble, Universite´ J. Fourier – CNRS, 38400 Saint Martin d’Heres, France q Universita’ di Urbino, Istituto di Scienze Chimiche ‘‘F. Bruner’’, 61029 Urbino, Italy r Department of Chemistry, University of Leicester, Leicester LE1 7RH, UK s University of Bonn, Institute of Crop Science and Resource Conservation – Plant Nutrition, 53115 Bonn, Germany t Bundesforschungsanstalt fu¨r Landwirtschaft (FAL) Institut fu¨r Agraro¨kologie, 38116 Braunschweig, Germany u Research Institute for Nature and Forest, 9500 Geraardsbergen, Belgium v INRA, INA PG, UMR Environm & Grandes Cultures, F-78850 Thiverval Grignon, France w Department of Meteorology and Air Quality, Wageningen University and Research Centre, 6700 AA Wageningen, The Netherlands x Institute for Meteorology and Climate Research, Atmospheric Environmental Research (IMK-IFU), Forschungszentrum Karlsruhe GmbH, 82467 Garmisch-Partenkirchen, Germany y Soils, Agronomy and Spatialization (SAS) Unit INRA, 35042 Rennes, France z Finnish Meteorological Institute, 00560 Helsinki, Finland aa Dipartimento di Matematica e Fisica ‘‘Niccolo` Tartaglia’’, Universita` Cattolica del Sacro Cuore, 25121 Brescia, Italy ab Institute for Plant Ecology, Justus-Liebig-University of Giessen, 35392 Giessen, Germany ac Institute for Environment and Sustainability, The European Commission, Joint Research Centre, 21020 Ispra, Italy ad Istituto per la Protezione delle Piante – CNR, 50019 Sesto Fiorentino, Italy ae Botanical Institute, University of Copenhagen, 1353 Copenhagen K, Denmark af Hungarian Meteorological Service, 1675 Budapest, Hungary ag Istituto di Biologia Agroambientale e Forestale – CNR, 00015 Monterotondo Scalo, Italy ah Institute of Agricultural and Environmental Sciences, Estonian University of Life Sciences, 51014 Tartu, Estonia ai School of GeoSciences, University of Edinburgh, EH9 3JN Edinburgh, UK aj Department of Physical Sciences, University of Helsinki, 00014 Helsinki, Finland ak Department of Applied Environmental Science, Atmospher Science Unit, Stockholm University, 10691 Stockholm, Sweden al Atmospheric Science Program, Department of Geography, Indiana University, 47405-7100 Bloomington, USA am School of Earth, Atmospheric and Environmental Sciences, The University of Manchester, M13 9PL Manchester, UK
* Corresponding author. Tel.: þ44 (0) 131 445 4343; fax: þ44 (0) 131 445 3943. E-mail address: [email protected] (D. Fowler)
1352-2310/$ – see front matter Ó 2009 Published by Elsevier Ltd. doi:10.1016/j.atmosenv.2009.07.068 5194 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 an Department of Forest Ecology, Federal Research and Training Centre for Forests, Natural Hazards and Landscape, 1131 Vienna, Austria ao Max-Planck-Institut fu¨r Chemie, 55128 Mainz, Germany ap Department of Experimental Physics and Environmental Change Institute, National University of Ireland, Galway, Ireland aq Climate and Global Change Unit, Research and Development, Finnish Meteorological Institute, 00560 Helsinki, Finland ar Energy Research Centre of The Netherlands, 1755 ZG Petten, The Netherlands article info abstract
Article history: Ecosystems and the atmosphere: This review describes the state of understanding the processes involved in Received 29 January 2009 the exchange of trace gases and aerosols between the earth’s surface and the atmosphere. The gases Received in revised form covered include NO, NO2, HONO, HNO3,NH3,SO2, DMS, Biogenic VOC, O3,CH4,N2O and particles in the size 27 July 2009 range 1 nm–10 mm including organic and inorganic chemical species. The main focus of the review is on the Accepted 29 July 2009 exchange between terrestrial ecosystems, both managed and natural and the atmosphere, although some new developments in ocean–atmosphere exchange are included. The material presented is biased towards Keywords: the last decade, but includes earlier work, where more recent developments are limited or absent. Dry deposition Trace gas fluxes New methodologies and instrumentation have enabled, if not driven technical advances in measure- Resuspension ment. These developments have advanced the process understanding and upscaling of fluxes, especially Biogenic emissions for particles, VOC and NH3. Examples of these applications include mass spectrometric methods, such as Compensation points Aerosol Mass Spectrometry (AMS) adapted for field measurement of atmosphere–surface fluxes using micrometeorological methods for chemically resolved aerosols. Also briefly described are some advances in theory and techniques in micrometeorology. For some of the compounds there have been paradigm shifts in approach and application of both tech- niques and assessment. These include flux measurements over marine surfaces and urban areas using micrometeorological methods and the up-scaling of flux measurements using aircraft and satellite remote sensing. The application of a flux-based approach in assessment of O3 effects on vegetation at regional scales is an important policy linked development secured through improved quantification of fluxes. The coupling of monitoring, modelling and intensive flux measurement at a continental scale within the NitroEurope network represents a quantum development in the application of research teams to address the under- pinning science of reactive nitrogen in the cycling between ecosystems and the atmosphere in Europe. Some important developments of the science have been applied to assist in addressing policy questions, which have been the main driver of the research agenda, while other developments in understanding have not been applied to their wider field especially in chemistry-transport models through deficiencies in obtaining appropriate data to enable application or inertia within the modelling community. The paper identifies applications, gaps and research questions that have remained intractable at least since 2000 within the specialized sections of the paper, and where possible these have been focussed on research questions for the coming decade. Ó 2009 Published by Elsevier Ltd.
Contents
1. Introduction ...... 5196 1.1. Scale...... 5197 1.2. Reactivity of natural surfaces ...... 5197 1.3. Frameworks for analysis and interpretation of trace gas and aerosol exchange ...... 5197 1.4. Bi-directional exchange ...... 5198 1.5. Aerosols ...... 5198 1.6. Ocean–atmosphere exchange ...... 5198 1.7. Wet deposition ...... 5199 2. Reactive gaseous nitrogen compounds – oxidized nitrogen ...... 5199 2.1. Introduction ...... 5199 2.2. Emissions from soils ...... 5199
2.3. Emissions of NOy from plant surfaces ...... 5201 2.4. Canopy atmosphere interactions ...... 5201 2.5. Models and measurements ...... 5201
2.6. Exchange of HNO3,HONO,PAN ...... 5202 2.7. NOx production and emission from snow surfaces ...... 5205 2.8. Up-scaling and regional and global trends ...... 5205 3. Biosphere atmosphere exchange of ammonia ...... 5205 3.1. Introduction ...... 5205 3.2. Advances in measurement methods ...... 5206 3.3. Key controls on biosphere atmosphere exchange of ammonia ...... 5208 3.4. Effects of ecosystem type on ammonia biosphere–atmosphere exchange ...... 5209 3.5. Modelling surface–atmosphere exchange of ammonia ...... 5209 3.6. Dynamic simulation of ecosystem C–N cycling and ammonia fluxes ...... 5210 3.7. Integrating ammonia exchange processes ...... 5211 3.8. Future challenges for ammonia exchange ...... 5212 4. Sulphur dioxide ...... 5212 4.1. Introduction ...... 5212 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5195
4.2. Worldwide advances in SO2 flux monitoring and modelling ...... 5213 4.2.1. Asia ...... 5213 4.2.1.1. Sulphur dioxide deposition to soils ...... 5213 4.2.1.2. Micrometeorological measurements over vegetated areas ...... 5214 4.2.1.3. Long-term deposition studies and inferential modelling ...... 5214 4.2.2. North America ...... 5214 4.2.3. Europe ...... 5215 4.2.3.1. Long-term flux monitoring in the UK ...... 5215 4.2.3.2. Other recent European datasets ...... 5215 4.3. Control of surface uptake rates by leaf cuticular chemistry ...... 5216 4.4. Advances in deposition modelling ...... 5217 4.5. Future challenges ...... 5217 5. Ozone...... 5218 5.1. Introduction ...... 5218 5.2. Deposition rates ...... 5219 5.2.1. European forests ...... 5219 5.2.2. Crops ...... 5220 5.2.3. Grasslands ...... 5221 5.2.4. Other vegetated surfaces ...... 5221 5.2.5. Non-vegetated surfaces ...... 5221 5.2.5.1. Snow ...... 5221 5.2.5.2. Water ...... 5221 5.3. Non-stomatal deposition processes ...... 5222 5.4. Model development and validation ...... 5222 5.5. Risk assessment methods ...... 5223 5.6. Potential effects of climate change ...... 5224 6. Biogenic volatile organic compounds (BVOC) ...... 5225 6.1. Introduction ...... 5225 6.1.1. Volatile isoprenoids ...... 5225 6.1.2. Oxygenated volatile compounds ...... 5225 6.2. Environmental controls on BVOC emissions ...... 5225 6.2.1. Physiological and physico-chemical controls of emissions ...... 5225 6.2.2. Physico-chemical controls of emission in species lacking specific storage structures ...... 5226 6.2.3. Uptake and release of volatile compounds by vegetation ...... 5226
6.2.4. CO2-Dependence of emissions ...... 5227 6.2.5. Induced emissions ...... 5227 6.3. Contemporary difficulties in scaling BVOC emissions from leaf to ecosystem ...... 5227 6.4. BVOC fluxes over Europe, by compound and in relation to the needs of photochemical oxidant models ...... 5227 6.4.1. Flux measurement techniques ...... 5227 6.4.2. Isoprene ...... 5227 6.4.3. Monoterpenes ...... 5228 6.4.4. Sesquiterpenes ...... 5228 6.4.5. Methanol ...... 5228 6.4.6. Acetone and acetaldehyde ...... 5228 6.4.7. Other compounds ...... 5228 6.5. The EU large field campaigns in the Mediterranean area: from BEMA to ACCENT ...... 5228 6.6. Remote sensing of BVOC ...... 5229 7. Deposition and resuspension of aerosols onto and from terrestrial surfaces ...... 5230 7.1. Introduction ...... 5230 7.2. Review of new measurement approaches and instrumentation ...... 5231 7.2.1. Flux measurements of particle numbers (size-resolved or total), without information on chemical composition ...... 5231 7.2.2. Flux measurements of individual aerosol chemical species ...... 5232 7.3. Area sources of particles ...... 5232 7.3.1. Resuspension ...... 5232 7.3.2. Urban emissions of aerosols ...... 5232 7.4. Dry deposition of particles ...... 5233 7.4.1. Dry deposition rates to vegetation ...... 5233
7.4.1.1. Friction velocity (u*) ...... 5234 7.4.1.2. Surface roughness length (z0) and canopy morphology ...... 5234 7.4.1.3. Particle diameter (Dp) ...... 5235 7.4.1.4. Atmospheric stability (z ¼ 1/L)...... 5236 7.4.2. Parameterising and modelling deposition rates ...... 5236 7.4.3. Dry deposition rates to urban areas ...... 5236 7.5. Uncertainties ...... 5236 7.5.1. Uncertainties in the application of micrometeorological flux measurement techniques for deriving the local flux ...... 5236 7.5.2. Relating measured fluxes to surface exchange: flux divergence and the effect of chemical interactions ...... 5237 7.5.3. Interpretation of measurements for model verification ...... 5238 7.6. Future research needs ...... 5238 7.6.1. Deposition measurements and reporting ...... 5238 7.6.1.1. Standardisation of eddy covariance approaches and data analysis procedures ...... 5238 5196 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
7.6.1.2. Improved measurements in the accumulation mode ...... 5238 7.6.1.3. Understanding the effect of stability and leaf properties on deposition velocities ...... 5238 7.6.1.4. Filtering or accounting for chemical interactions and water uptake ...... 5239 7.6.2. Deposition models ...... 5239 7.6.2.1. Migration to a probabilistic approach ...... 5239 7.6.2.2. Improvement of modelling approaches ...... 5239 7.6.2.3. Impact of surface anisotropy on suspension & deposition ...... 5239 7.7. Conclusions – aerosols ...... 5239
8. Ecosystem–atmosphere exchange of the radiatively active gases – N2O and CH4 ...... 5240 8.1. Introduction ...... 5240
8.2. Global budgets of N2O and CH4 ...... 5240 8.3. Biological sources of N2O and CH4 ...... 5240 8.3.1. The biology of production and consumption of N2O and CH4 in soils and sediments ...... 5240 8.3.2. Distribution of active microbial populations in soils ...... 5241
8.3.3. N2O and CH4 fluxes from the main global ecosystems ...... 5241 8.3.4. Plant-mediated transport and production of N2O and CH4 ...... 5241 8.3.4.1. Methane from vegetation ...... 5241 8.3.4.2. Nitrous oxide from vegetation ...... 5241
8.4. New developments in measurements of N2O and CH4 and denitrification ...... 5242 8.4.1. Flux chambers ...... 5242 8.4.2. Micrometeorological methods ...... 5242 8.4.3. Comparison of eddy covariance with chamber methods ...... 5242 8.4.4. Recent methodological advances in measurements of total denitrification rates ...... 5243
8.5. Modelling of N2O and CH4 fluxes at site and regional scales: approaches, applications and uncertainties ...... 5243 8.6. Validation of models by landscape and regional scale measurements ...... 5244 8.7. Conclusions ...... 5244 9. Exchange of trace gases and aerosols over the oceans ...... 5245 9.1. New trace gas interactions at the air–sea interface ...... 5245 9.1.1. Case studies ...... 5245 9.1.1.1. Acetone (ocean uptake) ...... 5245 9.1.1.2. Methanol (ocean uptake) ...... 5246 9.1.1.3. Isoprene (ocean emission) ...... 5247 9.1.1.4. Halogenated organics (ocean emission and uptake) ...... 5247 9.1.1.5. Monoterpenes (ocean emission) ...... 5247 9.1.1.6. Alkyl nitrates (ocean emission) ...... 5247 9.2. Aerosols ...... 5247 9.2.1. Primary marine aerosol (PMA) source functions ...... 5247 9.2.2. Chemical composition of primary sea spray ...... 5248 9.2.3. Secondary aerosol production ...... 5250 10. The processes of wet scavenging of aerosols and trace gases from the atmosphere ...... 5251 10.1. Introduction ...... 5251 10.2. Nucleation scavenging of drops and ice crystals ...... 5251 10.3. Impaction scavenging of aerosol particles ...... 5251 10.4. Scavenging of gases ...... 5252 10.5. Clouds ...... 5252 10.6. Orographic precipitation ...... 5252 10.7. Organic N in air and rain ...... 5253 10.8. Conclusions and some priority areas of future research ...... 5253 11. Ecosystem–atmosphere exchange – concluding remarks ...... 5254 11.1. Policy needs ...... 5254 11.2. Current understanding ...... 5254 11.3. Future developments ...... 5254 Acknowledgements ...... 5255 References ...... 5255
1. Introduction exchange processes is a core activity in understanding the Earth system. The subject of this review is much narrower than the The composition of the earth’s atmosphere is unique in the scope of these opening lines, and is restricted to the trace gases solar system in being largely determined by biological processes and aerosols exchanged between the atmosphere and the earth’s in soils, vegetation and the oceans interacting with physical and surface. However, as is clear from much of the international chemical processes within the atmosphere. The physical surface– assessment of changes in atmospheric composition since the atmosphere exchange of most gases contributing major and trace industrial revolution, these trace atmospheric constituents are constituents of the atmosphere is coupled to biological produc- changing the earth’s climate (IPCC, 2007), global biodiversity tion processes and transferred through the surface–atmosphere (Millenium Ecosystem Assessment, 2005) and the biogeochemical interface. Thus, developing a mechanistic understanding of the cycling of major nutrients including nitrogen, carbon, and production and destruction processes and their interactions with sulphur. D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5197
The earth’s surface is a sink for some atmospheric trace gases the photochemical oxidants in the 1960s and 1970s (Husar et al., and aerosols, and a source for many others, and for most the 1978). These early studies were made to determine the importance of surface–atmosphere interface represents a zone within which most surface removal which is better known as dry deposition (to distin- of the overall control of fluxes occurs. An understanding of the rate guish the process from removal by precipitation) as a sink for reactive controlling processes at this interface is therefore vital in describing trace gases (NO2,SO2,O3). Advances in understanding and computing the exchange process and understanding the global biogeochemical resources have allowed more sophisticated approaches to be adop- cycles. Applications of science in this field are necessary to quantify ted, in which the processes at the surface recognised different sinks and model responses to human perturbation of many of the and interactions with other trace gases, allowing rates of dry depo- biogeochemical cycles (C, N, S, halogens and metals to name but sition to change with time and with surface characteristics. a few). These perturbations include changes in land use or emis- sions of trace gases to the atmosphere, through combustion and 1.1. Scale industrial activities. Taking as an example the global nitrogen cycle, human activity through combustion processes for oxidized nitrogen Emission or deposition schemes to quantify trace gas fluxes and the Haber Bosch process for reduced nitrogen now dominates operate at a range of scales depending on the applications, illustrated the cycling of reactive nitrogen through the atmosphere and back to in Fig. 1.1. For hourly integration the application is primarily terrestrial and marine ecosystems (Galloway et al., 2004). The total for research purposes and mechanistic study at the small scale 2 2 emission of reactive nitrogen (Nr) from human activities at the end (<10 m ). For landscape scale measurements and for assessment of of the 20th century exceeds that from natural processes by a factor the fate of pollutants at the regional scale (106 km2) the application of 4 (20.7 Tg of oxidized and reduced reactive nitrogen Nr from has both research and policy objectives. At this scale, spatial and natural sources within a total of 104 Tg-N in 1993, Galloway et al., temporal integration provides robust parameterisation. The appli- 2004). As nitrogen is a limiting nutrient in many ecosystems, cation in global models to quantify sources and sinks is restricted in these modifications of the natural cycling have profound effects on spatial resolution, typically to 1 1 , and likewise has research and ecosystem function, biodiversity and atmospheric composition policy application Dentener et al., 2006. For the landscape scale, flux (Erisman et al., 2008). The human disturbance of the global carbon measurements may be made directly, using micrometeorological cycle is also extensive, and the quantities involved are very large. methods above canopies of vegetation, soil, or even ocean surfaces Global emissions of CO2 from fossil fuels since 1700 amount to and have become the method of choice for long-term flux approximately 600 Gt-C, which have increased the atmospheric CO2 measurement. These techniques provide, in addition to the target mixing ratio from 280 ppm to 380 ppm in 2006, an increase of about trace gas flux the turbulent exchange characteristics of the under- 30%, (IPCC, 2007). lying surface and the partitioning of available radient energy which These high level indicators of human influence provide essential enables the processes to be investigated at a sufficiently large scale to context for this review paper, but conceal the detailed changes integrate canopy scale fluxes over typically 105 m2 (Baldocchi et al., taking place and the range of chemical species and interactions 2001). involved. The subject area includes many different chemical species, The sections focus mainly on individual trace gases or classes of and it is not possible to be comprehensive in this review for all gases. atmospheric particles, and each considers the surface–atmosphere In particular the subject of the global carbon cycle and CO2 in exchange over specific ecosystems. The exceptions are the ocean particular are much too large to cover in this review. The focus of this surfaces and wet deposition, within which a range of relevant review is on the reactive trace gases and for the greenhouse gases, compounds is considered. CH4 and N2O. The gases are associated with a range of biological sources and have varied chemical reactivity in the atmosphere and at 1.2. Reactivity of natural surfaces surfaces. These differences reveal the range of controls and temporal and spatial variability in rates of exchange, which are the focus of the For many of the short lived gases (<2 days in the boundary review. The review moves through a range of chemical species, layer) there are multiple sinks at the surface and these include identifying the current state of knowledge and, where possible foliar surfaces and soil whose properties as sinks for a range of applications of the new developments in a policy context. gases vary with humidity and the presence of surface water and The gases emitted from terrestrial and ocean ecosystems are influenced, sometimes strongly, by the presence of other gases include all of the major greenhouse gases, H2O, CO2,CH4 and N2O, (Fig. 1.2). The chemical and physical complexity of terrestrial the nitrogen gases (both in reduced and oxidized forms), sulphur surfaces, illustrated in Figs. 1.1 and 1.2 at the microscopic scale is compounds, volatile organic compounds (VOC) and halogens. greatly simplified in the parameterisations used in models. The Quantifying the fluxes of these trace atmospheric components is simplification is necessary in part due to the nature of the flux a prerequisite within an assessmemnt process leading to the devel- measuring systems, which integrate the net fluxes over large areas opment of policy on climate change, eutrophication, acidification of these surfaces, and fail to reveal the microscopic scale of and photochemical oxidant formation. Many research groups have variability of the true exchange. become involved in the measurement and modelling of emission and uptake (deposition) fluxes of trace gases and particles. The mecha- 1.3. Frameworks for analysis and interpretation nistic understanding has developed from two different fields of of trace gas and aerosol exchange study, the first was concerned with the sources of atmospheric trace constitiuents, and the greenhouse gases were among the first The measurements of surface–atmosphere exchange provide at compounds for which surface fluxes were quantified directly by the simplest level the mass exchange per unit area of surface, which field measurements. These included small scale (0.1 m–0.5 m2) may be ground, water or leaf area, per unit time. To extract useful measurements of fluxes from soils and vegetation using chamber information on the underlying processes it is necessary to quantify methods for CO2, CO, CH4,N2O(Junge, 1963). The measurements the contributions each step in the transfer pathway makes to the showed large spatial and temporal variability so that up-scaling to overall exchange between defined points, which in this scheme is regions generated very large uncertainties. The other development simplified to vertical levels between a source and a sink. The most was mainly associated with the atmospheric transport and deposi- widely applied transfer scheme is a resistance analogue (Monteith tion of pollutants, including nitrogen and sulphur compounds and and Unsworth, 2007), in which the flux of trace gas or particle 5198 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
Fig. 1.1. A diagrammatic representation of the scales of measurement of trace gas fluxes for process studies (a transverse section through a Phaseolus vulgaris leaf, showing the palisade and mesophyll cells bounded by epidermal cells and the airspaces for internal exchange between trace gases and intercellular fluid). The field scale at which most of the micrometeorological flux measurements are made and the continental scale where models provide the emission and deposition fluxes. In this case the emission fluxes of oxides of nitrogen over Europe are shown, revealing the importance of international shipping to emissions over continental scales. is treated as an analogue of electrical current flowing through 1.4. Bi-directional exchange a simple network of resistances, which may act in series if there is just one sink at the surface, such as a pure water surface, or may For many of the trace gases, regardless of their reactivity, the have several sinks at the surface, acting in parallel, each repre- exchange fluxes may vary in sign as well as magnitude, with senting a distinct chemical component of the underlying surface. emission and deposition being commonly observed. The most A simple resistance network representing three different sinks widely known example of bi-directional exchange is CO2, which at the surface, and the two atmospheric resistances (Ra and exhibits both deposition and emission fluxes due to photosynthesis Rb, respectively the turbulent transfer resistance and the leaf and respiration respectively. In this case the concept of compen- boundary-layer resistance) are illustrated in Fig. 1.3. sation points as mixing ratios at which no net exchange takes place The atmospheric resistances may be separated from the total is now widely recognised for a range of trace gases (NH3, NO, CO2) resistance using independent measures of the turbulence above the but all controlled by very different processes. vegetation. The overall flux may be measured by a variety of micro- The recognition of bi-directional exchange requires modelling meteorological methods (Monteith and Unsworth, 2007), and thus approaches to simulate the process for application in surface– the total of the surface or canopy resistances to transfer between the atmosphere exchange schemes, as illustrated for NH3 in Fig. 1.4. source and sink may be quantified as the residual, as shown in Fig.1.3. 1.5. Aerosols
The understanding of deposition and emissions of aerosols over terrestrial surfaces has advanced considerably in the last decade, after a long period in which application of a model developed by Slinn (1982) has been a standard for many modelling approaches. Likewise, the emission of aerosols by resuspension from terrestrial surfaces has advanced following innovative new measurement approaches described in this review.
1.6. Ocean–atmosphere exchange
For many years the ocean–atmosphere exchange of trace gases has been treated in a simplistic way (Liss et al., 1981), in part due to the relative simplicity of the ocean surface relative to terrestrial surfaces, but also due to the limited knowledge base to support more complex treatments. However there has been an accelerating interest in ocean– Fig. 1.2. Illustrating the importance of different sinks for reaction of trace gases at the atmosphere exchange as new techniques have become available terrestrial surfaces, notably the external surfaces of vegetation often as in this case to make the flux measurements and as very new issues have been covered by complex layers of epicuticular wax and illustrated in Fig. 1.1, the internal structure of leaves, following uptake through the stomatal apertures and soils greatly identified. Current interest in oceanacidificationandoceaneutro- simplified in this illustration. phication further raise the profile of ocean–atmosphere exchange, and D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5199
deposition footprint of pollutants (e.g. with reference to acidification, χ O3(z-d) eutrophication, metal deposition or photochemical oxidants). These Ra applications require an understanding of the fate and transformations of the emitted trace gases within the atmosphere. The developments
RbO3 in understanding atmospheric chemistry of the trace gases discussed R ++= RRR in this paper are reviewed in the companion paper (Monks et al., this F t cba O3 issue)inthisvolume. R c2 1 R cuticle R R = c3 c1 c 1 1 1 2. Reactive gaseous nitrogen compounds – oxidized nitrogen soil stomata ++ Rc1 Rc2 Rc3 2.1. Introduction
χ ' O3(z 0 )= 0 Developments in understanding surface–atmosphere exchanges of NO and NO2 over the last decade have focussed on three specific issues: the long-term emission of NO from soils; the interaction of Fig. 1.3. A simple resistance analogy for a trace gas with sinks in stomata, on foliar chemical processing of nitrogen oxides in and above plant canopies; surfaces and in soil. and the deposition of NO2 and HNO3 to foliar and soil surfaces. The measurements have been made over different vegetation, but the given that these surfaces cover 71% of the earth’s surface, the relatively recent focus has been on forests, in part because the interactions small section of this review paper on this topic belies its importance in between these processes are greatest for forests, but also because understanding atmospheric composition change. some of the measurements are simpler to make and interpret for mature forests. This section outlines the developments in under- 1.7. Wet deposition standing NOy exchange between terrestrial ecosystems and the atmosphere, concentrating on developments during the last decade. Process understanding in the scavenging of gases and particles by precipitation has continued to advance, with important devel- 2.2. Emissions from soils opments during the last decade. The applications of wet deposition schemes are very important in the Long-Range Transport (LRT) Soil surface emissions of NO are the result of several biological models and increasingly in global chemistry-transport models and abiotic processes in the soil producing and consuming NO. (CTM) (Stevenson et al., 2006). These two applications make very Production and consumption of NO occurs predominantly via the different demands on available knowledge and understanding. In biological nitrification and denitrification processes. Nitrification is þ the case of LRT models in Europe (e.g. EMEP), the applications the oxidation of soil NH4 to NO3 , and denitrification is the anaerobic form part of the integrated assessment process and within the user reduction of soil NO3 to N2O and N2. In nitrification, NO is formed as þ community the pressure to provide ever finer spatial scale esti- a by-product during the oxidation of NH4 to NO2 and possibly also as mates of inputs presents challenges in the capability of LRT models a result of nitrifier reduction of NO2 leading to an NO production of þ and the meteorological models on which they depend. Current 1–4% of the NH4 being oxidized (Skiba et al.,1997). The NO produced demand for assessments of effects at the 1 km 1 km scale allows may be transformed within the soil profile by oxidation to NO3 or it the scale of the input estimate to approach the scale of an individual may be released to the atmosphere following diffusion to the nature reserve, for example. soil surface. In denitrification, NO occurs as an intermediate in the The focus of this paper on surface–atmosphere exchange processes cascade of reductive processes, and in the soil profile NO reduction spans a wide range of trace gases and particles. The motives for studies may contribute to the formation of N2O. Abiotic production of of many of the specific gases were environmental policy related, for NO occurs from oxidation of nitrous acid (HONO) that has been example to determine the atmospheric lifetime, travel distance and produced by protonation of biologically formed NO2 (Venterea et al.,
Fig. 1.4. A diagrammatic representation of bi-directional exchange, for NH3 exchange between the atmosphere and vegetation. 5200 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
2005). Under certain conditions e.g. after application of anhydrous ammonia to agricultural soils or acidic forest soils, the coupled biological-abiotic production of NO may constitute the dominant process for soil NO emissions (Venterea and Rolston, 2000; Go¨dde and Conrad, 1998). Factors that increase nitrification and denitrifi- cation, e.g. substrate and O2 availability, temperature and pH are thus predicted to influence NO formation. Likewise, factors affecting transport processes in the soil are predicted to regulate emissions of NO (and other gases). It has been hypothesized (Davidson,1991)that where WFPS (water-filled pore space) is less than 0.6, nitrification is the dominant process and relatively high emissions of NO may be observed. Under more reducing conditions, 0.6 < WFPS < 0.9, denitrification dominates which has a higher potential for NO production compared to nitrification (Skiba et al., 1997); however under conditions where anoxic conditions are generated by high soil water content or by compaction of fine textured soil the probability of NO being re-consumed by the denitrifying community is greatly enhanced. Soil water may also play a central role in mediating chemical processes leading to NO formation (Venterea et al., 2005). Under most soil conditions, both nitrification and denitrification occur simultaneously and the net flux of NO between soil and atmo- sphere is the result of both processes together. As current views of controls over NO gas emissions are still incomplete and need revision e.g. with emphasis on the role of abiotic formation (Venterea et al., 2005) there is a continuous need to further develop and improve methodologies to identify and characterise the NO formation processes. Go¨dde and Conrad (1998) achieved this by a combined modelling and experimental approach to determine the net NO flux in relation to NO concentration to quantify production and consumption rate constants and compensation concentration. Recent advances in methodological approaches to deepen our understanding of soil based NO emissions have included application of stable isotope techniques. Stark et al. (2002) applied a 15N-isotope pool dilution method- to obtain the simultaneous gross rates of NO forming processes combined with soil emissions, and Russow et al. (2000) adopted a kinetic isotope method (KIM) to study the complex N transformation processes involved in soil NO emissions. Fig. 2.1. Left: NO emission (mgNm 2 h 1) as a function of nitrogen deposition NO and N2O emissions were measured continuously at 15 forest (g N m 2 s 1). Regression lines (solid ¼ significant, dashed ¼ non significant) for 2 1 sites in Europe (Pilegaard et al., 2006) including coniferous and coniferous and deciduous sites, respectively. Right: N2O emission (mgNm h )as deciduous forests in different European climates, ranging from a function of C/N ratio. The full line represents a linear regression and the dotted line boreal to temperate continental forests and from Atlantic to Medi- a logarithmic regression. terranean forests. Furthermore the sites differ in atmospheric N- deposition ranging from low deposition (0.2 g N m2 yr 1) to high 2 1 deposition (4 g N m yr ). rate of O2 supply and thereby controls whether aerobic processes The relationships of the emissions of NO and N2O, with the such as nitrification or anaerobic processes such as denitrification parameters, nitrogen deposition, forest type, age, C/N in the surface dominate within the soil. While N2O emissions are known to horizon, pH, soil temperature and water-filled pore space (WFPS) increase at higher water contents through larger losses from were investigated by means of stepwise multiple regression anal- denitrification the relationship between the NO flux and the soil ysis. NO emission was dependent on forest type and positively water is more complex. Due to limited substrate diffusion at very correlated with nitrogen deposition (Fig. 2.1). WFPS was tested for low water content and limited gas diffusion at high water content, curvature by including a quadratic term, but this was not signifi- nitric oxide emissions are suspected to have a maximum at low to cant. Separately performed regression analyses for deciduous and medium soil water content. coniferous forests showed that the relationship between nitrogen The effects of soil moisture and temperature on NO and N2O deposition and NO emission was only significant for the coniferous emission were studied in laboratory experiment with soil cores forests: (NO (mgNm 2 h 1) ¼ 13.9 þ 25.5 [N deposition from a range of field sites (Schindlbacher et al., 2004). Soil moisture 2 1 2 (mgNm h )], r ¼ 0.82). The N2O emission was significantly and temperature explained most of the variability in NO emission negatively correlated with both the C/N ratio and the age of the (up to 74%) and N2O (up to 86%) emissions for individual soils. NO stands; a logarithmic transformation of N2O emission improved the and N2O were emitted from all soils except from a boreal pine forest significance of the correlation. soil in Finland, where the laboratory experiment showed net NO Soil temperature is a key variable affecting the emission rates of consumption. NO emissions from a German spruce forest ranged both gases (Fig. 2.2). Emissions of both NO (Slemr and Seiler, 1984) from 1.3 to over 600 mg NO–N m 2 h 1 and greatly exceeded and N2O(Skiba, 1998) increase with soil temperature due to the emissions from other soils. Average N2O emissions from this soil 2 1 positive effect of temperature on enzymic processes as long as tended also to be largest (170 40 mgN2O–N m h ), but did not other factors (e.g. substrate or moisture) are not limiting. Soil water differ significantly from other soils. NO and N2O emissions showed þ acts as a transport medium for NO3 and NH4 and influences the a positive exponential relationship with soil temperature. D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5201
400 In general, relationships between nitrogen oxides emission and soil moisture and soil temperature can be found within a single locality when studying short-term variations. However, when comparing ]
-1 annual values from different localities within a large region, parame- h 300 -2 terisations differ appreciably between locations and other factors such as soil properties, stand age, and site hydrological conditions interfere. This scale of variability has slowed the development of improved soil NO emission inventories, and their application in global chemistry- 200 transport models. So while it is disappointing to observe current generation CTMs using 1995 soil NO estimates, it is understandable and identifies a need for better parameterisations and soil emission 100 databases for global application.
2.3. Emissions of NO from plant surfaces NO emission [ μ g NO-N m y
0 Production of NOy on Scots pine branch surfaces by ultraviolet -5 0 5 10 15 20 radiation has been observed in Hyytia¨la¨, southern Finland (Hari Forest floor temperature [°C] et al., 2003)(Fig. 2.4). Other studies have shown that irradiance- dependent NOy emissions from snow and different chamber surfaces Fig. 2.2. Relationship between daily mean forest floor temperature and daily mean NO have been observed to originate from HNO3 or nitrate photolysis. In emissions at the Ho¨glwald Forest (spruce, control) for the observation period January 1, Hyytia¨la¨ forest, Raivonen et al. (2006) investigated whether the NOy 2004–December 31, 2006. For details on measurement and site characteristics see Gasche and Papen (1999). emitted from pine shoots could originate from photolysis of HNO3 attached to the needle surface. Field data of several years from Hyytia¨la¨ were used to test this hypothesis. The HNO3 deposition, estimated for the Hyytia¨la¨ site, has been high enough to account for The results from the annual averages of field data did not show the NOy emission rates observed from the chambers. The particular significant relationships with soil temperature for either NO or for characteristics of the daily pattern of CO2 exchange or stomatal N2O emission. Schindlbacher et al. (2004) showed that N2O emis- control was not reflected in the NOy flux. When a pine branch sions increased with increasing WFPS or decreasing water tension, was rinsed, which reduced the amount of water-soluble nitrogen respectively. Maximum N2O emissions were measured between 80 compounds (e.g. HNO3, nitrates and HONO) from the needle surface, and 95% WFPS or 0 kPa water tension. The optimal moisture for NOy emissions from that branch decreased compared to another NO emission differed significantly between the soils, and ranged non-rinsed branch. Therefore, it was concluded that the results between 15% WFPS in sandy Italian floodplain soil and 65% in loamy support the hypothesis and that HNO3 photolysis on plant surfaces Austrian beech forest soils. For the field data WFPS was not needs to be taken into account both from air chemistry and plant a significant parameter for N2O emission, but had a positive signif- sciences point of view. icant effect on NO emission (Fig. 2.3). The annual average WFPS in the field was higher than the optima found for NO in the laboratory 2.4. Canopy atmosphere interactions experiment, but since not all field sites were studied in the labora- tory it is difficult to provide a general conclusion. The inter-annual The interaction between chemical reactions of nitrogen oxides variation within single sites clearly showed relationships to both taking place in the canopy and trunk space of a forest is a special temperature and soil moisture. An important factor for N2O emis- case because in this area chemical and turbulent timescales change sion is freeze-thaw events which can produce a significant outburst substantially leading to a very complex situation in which even the of N2O(Kitzler et al., 2006). direction of fluxes may change (Duyzer et al., 1995)(Fig. 2.5). This makes it nearly impossible to interpret measurements of the turbulent fluxes of some reactive trace species above the canopy from single point eddy covariance measurements. Several models have been Sandy loam developed to simulate the overall exchange and show the magnitudes Silty loam of the different competing processes. These models describe the Sandy clay loam Loam(1) coupled processes of atmospheric transport and chemical processes Loam(2) above and in canopies in detail. Over forests the situation is even more complex. Flux measurements are usually carried out near the top of rough canopies leading to potential inaccuracies in the K-theory approximation. This theory is relatively easy to combine with vertical atmospheric transport phenomena with fast chemical reactions.
2.5. Models and measurements
Measurementsof small fluxes of NO and NO2 have shown spurious NO emission (% of maximum) results especially at low concentrations due to a lack of specificity of
0 20406080100 monitors and a lack of instrument sensitivity, but other problems may well have contributed, including violation of conditions under which 020406080100 such fluxes may be measured above canopies and the complexity of WFPS(%) interactions; soils and sunlight driven reactions may both be sources Fig. 2.3. The relationship between NO emission and water-filled pore space at different of NOy and these interact with the stomatal sinks and the chemical localities in the NOFRETE project (based on data in Schindlbacher et al., 2004). processing within the canopy trunk space Duyzer et al., 1983. 5202 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
Fig. 2.4. Effect of UV radiation on the NOy concentration in a small Teflon chamber that enclosed a clean pine branch or a branch that had been treated with NH4NO3 solution. The branches were dead and dry, cut from the tree. UV wavelengths were filtered away using a Plexiglas plate (Raivonen et al., 2006).
As a result of these limitations there are only limited data a compensation point above which the flux is directed towards the available for verification of models. Duyzer et al. (2004) described surface and below which the flux is away from the surface. the analyses of a dataset acquired in the framework of a European In summary the flux of NO2 above a forest can be described with project from an experiment carried out in a 20 m high coniferous the following function: forest (Speulderbos, The Netherlands). A 1D multi-layer model of a forest canopy was used to analyse the field data. In each layer 1 FNO2 ¼ f FNO;soil CNO2 vertical transport was described using K-theory; canopy uptake Rc;NO2 was described using a resistance layer model. Simple chemical where all variables have their common meaning and C denotes reactions between ozone and nitric oxide and photolysis of NO2 the concentration of nitrogen dioxide above the canopy. This nitrogen dioxide were described. The coupled differential equa- equation is rather qualitative but indicates the sensitivity of the flux tions were solved numerically. Input to the model calculations were of NO above the canopy. More quantitative model runs are needed, concentrations of nitrogen oxides and ozone at the highest level 2 but these require a large amount of input data and the results are above the forest, levels of radiation, temperature, humidity, wind still uncertain. speed, turbulence parameters and an estimate of the emission of nitric oxide. Output of the model is the concentration and fluxes of A simple resistance model (Duyzer et al., 2005a,b) was tested in the relevant components at the height of each level in and above a deciduous forest (Sorø, Denmark) and is illustrated in Fig. 2.6. the canopy. These may be compared with measured fluxes of these Generally the understanding of the various processes and their components at two levels above and one below the canopy. It is fair interaction is increasing. Nevertheless many uncertainties to say that the comparison between measured and modelled fluxes remain and there is a need for further improvement of models, is not impressive. There are many possible explanations for this especially for Lagrangian models incorporating chemical reac- observation but no clear single cause has been identified. tions. On the other hand, the accuracy of the results of field Depending on the magnitude of this soil flux the NO2 flux is measurements has been rather low. It should be noted that either downward or upward. In the case of the coniferous forest although the interaction between atmospheric chemical reac- described here and the conditions during the experiment the result tions and exchange between the canopy and the atmosphere is was that when the NO flux from the soil exceeded 10 ng m 2 s 1, easy to understand its importance may be limited. In cases where the NO2 emission was upward (i.e. away from the forest). fluxes of nitrogen oxides are small the corrections could be large At high concentrations the NO2 flux is directed towards the in a relative sense but still rather small in an absolute sense. The forest and at small concentrations the flux is more likely to be currently available models could very well be capable of making directed towards the free atmosphere. This may be interpreted as estimates of the magnitude of these effects. In view of all the uncertainties hindering improved estimates in testing of models the limited quality of the description of atmospheric transport processes within the canopy may not be a serious problem here.
2.6. Exchange of HNO3, HONO, PAN
The deposition of HNO3 to terrestrial surfaces has been shown to be primarily controlled by the rates of turbulent exchange in the atmospheric boundary layer and the leaf boundary layer (Huebert and Robert, 1985). The highly reactive and soluble nature of gaseous HNO3 leads to large rates of deposition, approaching the maximum rates of deposition limited by turbulent exchange when each molecule arriving at terrestrial surfaces is immediately absorbed at the surface. In these conditions the surface is considered to be acting as a perfect sink, canopy resistance is zero and the numerical value for the deposition velocity becomes:
Fig. 2.5. A schematic of the various canopy interactions involved in the exchange of V ¼ V ¼ 1=r þ r nitrogen oxides between the free atmosphere and forests. gðNHO3Þ max a b D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5203
Fig. 2.6. Profiles of CO2,NO2,O3 and NO in a beech forest near Sorø, Denmark. The profiles clearly show the effects of stable conditions during night and daytime turbulence mixing the full air column. CO2 is built up during night due to soil respiration; O3 is depleted during night due to deposition and chemical reactions. At the soil surface high concentrations of NO are seen due to emission from the soil.
The values for deposition velocity in these conditions are very US pine forests at Duke Forest, North Carolina, (RH > 75%) and sensitive to wind velocity values and approach several cm s 1 even Blodgett Forest, California, (RH < 30%) (Farmer et al., 2006; Turn- over relatively short vegetation. The consequence of these large rates ipseed et al., 2006; Wolfe et al., 2008). of deposition are that even in areas with small HNO3 concentrations, At Duke Forest fluxes of PAN, PPN and MPAN were measured with dry deposition contributes a substantial quantity if nitrogen. Taking an a CIMS technique (Turnipseed et al., 2006). There were no significant ambient concentration of 0.1 ppb HNO3, the annual deposition of N for differences in the Vd of the three different PAN compounds, but all 1 a forest would be of the order 3 kg N ha annually from HNO3 alone. three species deposited about four times faster than predicted by the The close coupling between rates of turbulent exchange and model of Wesely (1989) during the day, and nearly an order of dry deposition rates for HNO3 also generates substantial spatial magnitude faster during the night, indicating that aqueous solubility variability in N deposition in the landscape, with hotspots for N considerations are insufficient to predict the behaviour of PAN on 1 1 deposition being forests and especially forest edges, hedgerows and surfaces. The average Vd was 2.5 mm s during day and 8 mm s isolated, exposed hills, where wind speeds are larger. during night. In contrast to the considerations of Wesely (1989),wet Several studies have recently attempted to measure total surfaces showed a smaller non-stomatal resistance (Rns ¼ 125 s m 1) 1 oxidized nitrogen (NOy) fluxes or even total reactive nitrogen than dry surfaces (Rns ¼ 250 s m ). (Nr ¼ NOy þ NHx) to ecosystems (Turnipseed et al., 2006). These At the much drier Blodgett forest site, the flux of the sum of all approaches offer the prospect to apply eddy covariance techniques PANs was measured by TD-LIF, based on thermal conversion and NO2 for the robust and relatively cost effective determination of total detection (Farmer et al., 2006). PAN was derived as the difference atmospheric N deposition, but they do not provide the chemical between the ambient temperature and 180 C channel. They found speciation needed to further process understanding. There have, upward fluxes in summer and on average bi-directional exchange however, also been advances in the understanding of individual N with afternoon deposition in winter, when noon-time deposition 1 compounds other than NH3, NO and NO2. velocities averaged 8 mm s . More recently, these measurements A recent lab study (Sparks et al., 2003) has confirmed that PAN have been repeated with the more selective CIMS technique (Wolfe deposition through the stomata can make a significant contribution et al., 2008). Here measurements indicated larger average midday 1 to plant uptake of atmospheric N. In addition, recent instrument values of Vd for PPN (12 mm s ) than for the PAN and MPAN developments in chemical ionisation mass spectroscopy (CIMS) (4 mm s 1), while both compounds deposited slowly at night 1 and thermal-dissociation laser induced fluorescence (TD-LIF) have (Vd < 2mms ). The authors of this study attribute the difference in enabled the application of eddy covariance to the biosphere/ the Vd between compounds to MPAN and PAN production inside the atmosphere exchange of preoxy acyl nitrates (PANs such as PAN, canopy and suggest that the PPN fluxes are a better descriptor of the PPN and MPAN). Measurements were made over two contrasting surface deposition. They suggest that the non-stomatal uptake is 5204 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 dominated, but not fully explained, by thermochemical decompo- governed by the thermodynamic equilibrium with NH4NO3 on leaf sition, and thus strongly linked to canopy temperature. surfaces or fertilizer pellets. Fig. 2.7 shows an example of reduced In summary, this recent measurement evidence suggests that deposition of HNO3 and HCl over a Dutch heathland. deposition rates of PANs in warm climates are at least a factor of 5 The view that HNO3 normally deposits with a near zero canopy larger than predicted by commonly used models and non-stomatal resistance still holds. There is an increasing measurement database deposition is larger to wet and humid surfaces than to dry surfaces. of HNO3 concentrations in national and regional networks suitable During the same TD-LIF study, Farmer et al. (2006) measured for inferential modelling of HNO3 deposition (Tang et al., 2009), fluxes of total alkyl nitrates (gas and aerosol phase), from the differ- which now provides independent confirmation from the model ence between the 180 Cand330 C channels. These compounds results, that HNO3 deposition makes a very significant contribution showed large winter-time midday deposition velocities of 20 mm s 1, to nitrogen deposition across Europe. In addition, Europe-wide 1 1 approaching those of HNO3 (25 mm s ). Even higher Vd of 30 mm s monitoring activities have produced the first hourly monitoring was derived by Horii et al. (2005) for what they interpret as isoprene- datasets of HNO3, which allows for a much more in-depth assess- derived hydroxyalkyl nitrates. ment of the performance of oxidized nitrate chemistry in atmo- Nitric acid (HNO3) has traditionally been believed to deposit at spheric transport models (Tarrason and Nyiri, 2008). the maximum rate possible according to turbulence (Vmax) and its Much development has occurred in measurement techniques flux measurement continues to be used to derive quasi-laminar for nitrous acid (HONO), e.g. based on long path absorption boundary-layer resistances for vegetation (e.g. Pryor and Klemm, photometry (LOPAP) and differential optical absorption spectrom- 2004). This view has been challenged by recent measurements etry (DOAS). This has contributed to the improvement of process that indicated non-negligible canopy resistances in the range of 50 understanding of sources of HONO in the atmosphere, e.g. revealing to >200 s m 1 during midday (Nemitz et al., 2004b, 2008). This larger daytime sources than previously thought and identifying has been attributed to non-zero chemical compensation points NO2 reactions with humic acid as a novel production mechanism
10
] 0 -1 s
-2 -10
-20 HNO3 HCl [ng m [ng χ
F -30
-40 15:00 18:00 21:00 00:00 03:00 06:00 09:00 12:00 15:00 18:00 21:00 00:00 30 V 25 max(HNO3) ]
-1 V 20 max(HCl) V (HNO ) 15 d 3 V d(HCl) [ mm s [ mm 10 d V 5 0 15:00 18:00 21:00 00:00 03:00 06:00 09:00 12:00 15:00 18:00 21:00 00:00 500
400 ] -1 300 HNO
[s m 3
c 200 R HCl 100
0 15:00 18:00 21:00 00:00 03:00 06:00 09:00 12:00 15:00 18:00 21:00 00:00 0.5 120 RH u 100 0.4 *
] 80 -1 0.3 ') [%] 60 0 z ( [m s [m
* 0.2
u 40 RH 0.1 20 0.0 0 15:00 18:00 21:00 00:00 03:00 06:00 09:00 12:00 15:00 18:00 21:00 00:00 GMT
Fig. 2.7. A time series of HNO3 and HCl exchange measured above a Dutch heathland with a denuder gradient system with online analysis by ion chromatography. The panels show: (a) fluxes, (b) deposition velocities of HNO3 and HCl in comparison with their maximum values and (c) Rc for HNO3 and HCl, (d) friction velocity (u*) and relative humidity (h). Data represent 2.5 h running means of 30 min. Vd(HNO3) and Vd(HCl) are reduced compared with their maximum values, presumably due to non-zero chemical compensation points þ originating from deposited NH4 salts. From Nemitz et al. (2004a). D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5205
(Kleffmann et al., 2005; Stemmler et al., 2006). By contrast, appli- based on intensity of rainfall) to estimate soil NO emissions. cations of these approaches to flux gradient measurements are still Yienger and Levy also provide a canopy reduction factor to include rare, which nevertheless confirm surface sources of HONO (Vitel estimates of the chemical conversion and re-deposition of NO as et al., 2002; Kleffmann et al., 2003). NO2 within the canopy. Compared to the methodology by Yienger and Levy, the Skiba-EMEP/CORINAIR approach is more simplistic. 2.7. NOx production and emission from snow surfaces Based on a literature review by Skiba et al. (1997) this approach postulates that 0.3% of any form of nitrogen is volatilized as NO, i.e. Snow lying on the Earth’s surface has traditionally been viewed regardless whether it originates from inorganic or organic fertil- as a chemically inert medium, whose influence on the overlying ization or atmospheric N deposition. Furthermore, a background atmosphere was exerted through its albedo effect, and by emission of 0.1 ng NO–N m 2 s 1 (z0.032 kg NO–N ha 1 a 1)was restricting exchange of gases between the air and land/sea surfaces. assumed (Simpson et al., 1999). In addition, EMEP/CORINAIR also The a priori view was that the boundary layer and troposphere over use a more detailed methodology (BEIS-2), which originates from Antarctica would be somewhat uninteresting, with low concen- the work of Novak and Pierce (1993) and considers soil temperature trations of reactive radicals such as OH, HO2, NO and NO2, and as well as different land use classes. A statistical summary model a composition dominated by longer-lived chemical species. The was developed by Stehfest and Bouwman (2006), which is based equivalent regions of the Arctic atmosphere were assumed domi- on the most extensive literature review currently available. nated by long-range transport of anthropogenic emissions from This methodology for calculating soil NO emissions on global and lower latitudes. However, recent research has shown that this regional scales considers land-use, N fertilization rate [Fertilizer], picture is far from the truth, and that snow is a highly photo- soil N content (three different classes, estimated as 1:10 of soil chemically active medium. Snow-pack impurities, of which there organic carbon content) [SON] and climatic regions. The method- are many, can be photolysed to release reactive trace gases to the ology was recently adapted to calculate a European wide inventory atmosphere. These processes are likely to be active anywhere that of NO emissions from forest soils (Kesik et al., 2006). Kesik et al. sunlight irradiates snow. The importance of these processes to used the process-oriented ecosystem model Forest-DNDC. The boundary layer composition varies with geographical location; in model was extensively tested for its performance to predict NO regions with a high background of radicals, for example arising emissions at the various NOFRETETE field sites, which were located from anthropogenic pollution, emissions from snow are of lesser across Europe and, thus, were covering different climatic condi- importance. But in the remote polar regions, emissions from snow tions (Pilegaard et al., 2006). Regionalisation was finally achieved can be the dominant source of reactive trace gases and have a major by linking the model to a detailed GIS database holding all relevant influence on boundary layer chemical composition. This conclusion information for initializing and driving the model such as data on was first reached for NOx (NO þ NO2), which was measured in vegetation (e.g. forest type) and soil properties (e.g. texture, soil pH, the boundary layer at Summit, Greenland at surprisingly high organic C content) and climate (either present day conditions or concentrations, and with a ratio to NOy that suggested a local projected future climate predictions). This approach demonstrated source. Measurements of NOx within the snow-pack interstitial air for the first time the huge regional differences in NO emissions revealed concentrations that were higher still, suggesting that the from forest soils across Europe as shown in Fig. 2.8, to estimate its snow-pack itself was the source, with a gradient to the atmosphere. significance on a regional scale and to unravel the importance of Subsequent measurements made in Antarctica confirmed that atmospheric N deposition for the magnitude of forest soil NO NOx production within the snow-pack was a feature of both emissions. polar regions (Jones et al., 2000). Additional measurements confirmed that NOx generated within the snow-pack was released 3. Biosphere atmosphere exchange of ammonia to the overlying boundary layer (Jones et al., 2001), contributing to the higher than expected NOx concentrations that were 3.1. Introduction encountered. The air chemistry of polar regions is a rapidly developing field Substantial progress has been made during the last five years in and extends substantially beyond interest in oxidized nitrogen, in understanding ammonia biosphere–atmosphere exchange. Experi- which surface processes clearly provide the major source of reac- mental studies have included controlled laboratory analysis, while tive oxidized nitrogen in Antarctic regions. Discussion of other trace a series of micrometeorological studies have assessed net fluxes gases, including halogen and mercury compounds are beyond the occurring under field conditions. In particular, major advances have scope of this review, however a valuable review of polar halogen been made in modelling the different aspects of ammonia exchange. chemistry and links to oxidized nitrogen chemistry is provided by This has included not just analysis of the drivers of the vertical flux Simpson et al. (2007). densities, but also a consideration of non-stationarities, such as advection effects and chemical interactions. Traditionally, micro- 2.8. Up-scaling and regional and global trends meteorological experiments were designed to avoid these effects, focusing as far as possible on ‘ideal’ micrometeorological conditions, The complexity of processes involved in NO emissions from soils so as to better quantify the vertical exchange processes, and develop has resulted in a significant uncertainty in the regional and global parametrisations for ‘dry deposition schemes’ in regional models source strength of soils for NO. However, different methodologies (Fowler and Duyzer, 1989; Fowler et al., 1998; Sutton et al., 1994; have been developed, e.g. relatively simple statistical models as Simpson et al., 2006). However, for ammonia, it has become well as process-based model approaches, to cope with the problem increasingly clear that these non-stationarities represent important of regionalisation of soil NO fluxes. The most widely used approach effects that are widespread in the real environment and need to be to calculate regional NO emissions from soils is based on the work quantified (Sutton et al., 2007, 2008c). of Yienger and Levy (1995). These authors consider land use and The developments in the last decade have arisen from a wide range respective background emission strengths, nitrogen fertilization of national and international projects. National studies have particu- rate (2.5% loss of applied nitrogen), temperature effects (three larly addressed exchange with key ecosystems of regional importance, classes: cold-linear, exponential and optima) as well as the pulsing such as ammonia losses from agricultural systems (e.g. Milford et al., of NO emissions following prolonged dry periods (four classes, 2001a; Walker et al., 2006; Wichink Kruit et al., 2007)andthe 5206 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
were complemented by the LIFE project, which added long-term ammonia flux data for a number of grassland, moorland and forest ecosystems (e.g. Flechard and Fowler, 1998; Erisman et al., 2001; Spindler et al., 2001). As understanding of ammonia exchange has improved and scientific ambition developed, increased attention has been given to integrating the different drivers of ammonia exchange processes. This has, for example, been reflected in the Braunschweig Integrated Experiment of GRAMINAE (e.g. Sutton et al., 2002, 2008a, 2009a,b), which linked a wide range of biospheric, atmospheric and manage- ment interactions as these control ammonia exchange with managed grassland. This integration has developed substantially under the NitroEurope Integrated Project (Sutton et al., 2007), in which inter- actions between the different components of nitrogen fluxes, including ammonia and oxidized nitrogen and their effect on the net greenhouse gas balance have been investigated. In parallel, major advances have been made in spatial modelling of ammonia fluxes, from individual forest edges to global scales (Theobald et al., 2004; Dentener and Crutzen, 1994; Hertel et al., 2006; Bleeker et al., 2006).
3.2. Advances in measurement methods
Before considering the developments outlined above in more detail, it is important to highlight that the advances have been critically dependent on improvements in measurement technology (See Table 1). At the start of the 1990s, ammonia flux measure- ments were still being made using wet chemistry and manual batch sampling with time integration of typically 2 h (e.g. Sutton et al., 1993a; Duyzer, 1994). The most important advance was the intro- duction of continuous wet chemistry methods for measuring ammonia profiles, including the AMANDA wet rotating denuder (Wyers et al., 1993) and the mini-Wet Effluent Diffusion Denuder (e.g. Blatter et al., 1993; Neftel et al., 1999). Although these tech- Fig. 2.8. Importance of atmospheric N deposition for NO emissions from forests soils. The niques are liable to malfunction, with effort and careful operation map shows the difference in NO emissions between a scenario with zero atmospheric they have produced many key datasets over the last 15 years N deposition and present day atmospheric N deposition. In large parts of central Europe but also Scandinavia forest NO emissions are likely to decrease significantly if atmo- (e.g. Erisman and Wyers, 1993; Sutton et al., 1995, 1997; Fowler spheric N deposition can be reduced to background levels. et al., 1998; Flechard and Fowler, 1998; Neftel et al., 1998; Milford et al., 2001a,b) and still represent the state-of-the-art as regards precise measurement of small ammonia fluxes (Wichink Kruit et al., ammonia inputs into semi-natural ecosystems of conservation value 2007; Neirynck and Ceulemans, 2008; Sutton et al., 2008b). (e.g. Wyers and Erisman, 1998; Neirynck and Ceulemans, 2008). A unique inter-comparison of four continuous wet chemical Collaborative international projects have sought to integrate and systems was made at the GRAMINAE Braunschweig Experiment extend these interests, making the comparison between ecosystem (Sutton et al., 2007, 2009b; Milford et al., 2009), which highlights types and looking at the interactions (e.g. Sutton et al., 2009a). the potential and limitations of the approach. Fig. 3.1 shows the The first European collaborative project dedicated to ammonia ammonia flux measured before and after cutting an agricultural exchange was ‘EXAMINE’. Attention was given to quantifying grassland, and following subsequent fertilization with calcium ammonia exchange with a range of European ecosystems, under both ammonium nitrate. The measurement systems were able to detect experimental and field conditions (e.g. Sutton et al., 1995; Schjoerring the wide range of ammonia fluxes, but the degree of agreement et al., 1998; Neftel et al., 1998; Meixner et al., 1996; Nemitz et al., varied greatly between days. This was a result of varying perfor- 2009b), including analysis of the surface gas-particle interactions mance of the different analysers, highlighting the need for highly between ammonia, nitric acid and hydrochloric acid (Nemitz and intensive instrument maintenance. Sutton, 2004). As part of EXAMINE a major collaborative analysis was Despite the improvements that have been made in the automation made in the North Berwick experiment, which provided a uniquely and reliability of the continuous wet chemical gradient methods (e.g. detailed examination of the processes controlling ammonia exchange Wichink Kruit et al., 2007; Flechard et al., 2007; Sutton et al., 2007), with an oilseed rape canopy (Husted et al., 2000; Nemitz et al., their remain several limitations, which have encouraged researchers 2000a,c; Sutton et al., 2000a,b). to seek alternative ammonia flux measurement approaches. In prin- In the second major European collaboration dedicated to ciple, many refinements have allowed the automated wet chemical ammonia, the GRAMINAE project analyzed the processes controlling methods to become more reliable and comprehensive (such as being ammonia exchange with grassland ecosystems across Europe (Sut- able to measure aerosol and acid gas gradients simultaneously, Trebs ton et al., 2001). This included assessment of both the bi-directional et al., 2006). However, the use of many moving parts can be consid- fluxes of ammonia with agricultural grassland – as these affect ered as inherently liable to faults. Similarly, the response times of atmospheric ammonia balance (e.g. Milford et al., 2001a; Mosquera these instruments are typically >5 min, which means that they et al., 2001), and with semi-natural grasslands as these are impacted are normally limited to the measurement of mean concentration by the atmosphere (e.g. Horva´th et al., 2005). These studies differences and vertical gradients. D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5207
Table 1 Practical suitability of systems to measure ammonia biosphere–atmosphere exchange.
Chemical approach Advantage Disadvantage Application References
Box AGM REA EC Simple, cheap, high Batch filter packs Uncertain gas – aerosol split, batch JJ LLSutton et al. (1993a,b) air volume
Simple, cheap good Batch denuders Low air volume, batch JJ LLDuyzer (1994) gas-aerosol split
Automated batch Automated in field, medium High laboratory processing cost, JJ LLLoubet et al. (2006, 2009) annular denuders cost, precise, high air volume only hourly, need two systems for fluxes
Continuous annular Automatic sensitive, Cost, complexity, fault liable, JJ JJ L L Wyers et al. (1993), Erisman and Wyers (1993), denuders precise, high air volume gradient only Sutton et al. (1995, 2000b, 2001a) and Nemitz et al. (2001b)
Continuous parallel Automatic, sensitive, Cost, complexity, fault liable J J JJ L Nemitz et al. (2001a) and Hensen et al. (2008) plate denuders high air volume REA
Continuous Automatic, Cost, complexity, fault liable JJJJLHensen et al. (2008) mini-WEDD sensitive, precise
Continuous Automatic, sensitive, Cost medium complexity JJJJLFlechard et al. (2007) and Hensen et al. (2008) membrane precise, reliable denuder AIRmonia
Photo-accoustic Automatic, sensitive, Cost, complexity, not reliable JJ J L L Whitehead et al. (2008) in principle reliable
Tunable diode laser Automatic, sensitive, Very high cost, complexity, JJ LJShaw et al. (1998), Famulari et al. (2004), fast response (>10 Hz) maintenance Twigg et al. (2005) and Whitehead et al. (2008)
AGM: Aerodynamic Gradient Method; REA: Relexed Eddy Accumulation, EC: Eddy Covariance.
The benefits of quantifying ammonia fluxes using the gradient technique have been clearly demonstrated by the many papers published using this approach. In terms of informing our under- 50 ) standing of ammonia exchange processes and model development,
-1 Pre-cut s CEH this has almost exclusively been provided by measurements -2 25 FRI using the aerodynamic gradient method (90%), with a few studies 0 (in continental climates) applying the modified Bowen Ratio method (5%). By contrast, the key disadvantage of this method is -25 flux (ng m that it depends on micrometeorological stationarity, with no 3 change in the vertical flux with height. These flux/gradient methods NH -50 are not suitable for the study of exchange fluxes where advection of -75 ammonia from local sources is of interest (e.g. Loubet et al., 2001, )
-1 22/05/00 23/05/00 24/05/00 25/05/00 2006) and where gas-particle ammonia–ammonium interactions s Post-cut -2 CEH are significant (e.g. Brost et al., 1988; Nemitz et al., 1996, 2004b). FRI To address some aspects of advection and air chemistry interac- 500 FAL-D FAL-CH tions, determination of fluxes at a single height offers a way forward. If this can be achieved, in principle, deployment of replicate flux (ng m (ng flux 3 measurement systems at several heights could then be able to
NH 0 determine vertical flux divergences (Sutton et al., 2007). Both the Eddy Covariance (EC) method and Relaxed Eddy Accumulation (REA) allow fluxes to be determined from measurements at one height, and
) 31/05/00 01/06/00 02/06/00 03/06/00 -1 have therefore been the subject of several recent studies. The s 6000 Post-fert -2 CEH advantage of REA is that slow response ammonia measurements can FRI be combined with fast response switching, as has recently been 4000 FAL-D FAL-CH demonstrated in an inter-comparison of 4 REA systems for ammonia (Hensen et al., 2008). A further advantage is that programmed flux (ng m 3 2000 periods of random switching between air up- and down-drafts
NH allows automatic zero checks and the correction of any bias (Nemitz 0 et al., 2001a; Hensen et al., 2008). By contrast, the challenge for REA 06/06/00 07/06/00 08/06/00 09/06/00 and ammonia is that the concentration differential to be measured is typically much smaller than for the gradient method, which to a large Fig. 3.1. Inter-comparison of continuous profile systems for measuring ammonia fluxes extent cancels the benefit of auto-referencing. by the aerodynamic gradient method (AGM), from the GRAMINAE Braunschweig Experiment. Although highly scattered, this flux inter-comparison is unique and repre- Several recent studies have demonstrated the potential of fast sents the current state-of-the-art in chemical detection systems for ammonia fluxes. Increased emissions due to cutting of the underlying grass sward (29 May) and the effect response tunable diode laser absorption spectroscopy (TDLAS) of N fertilization with (100 kg N ha 1, 5 June) are clearly shown (Sutton et al., 2007). for measurement of ammonia fluxes by eddy covariance 5208 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
1200 average concentration in the ecosystem canopy, both of which
TDLAS vary in time and in space (e.g. Sutton et al., 1995; Asman et al., 1000 QC-TDLAS 1998). Within the canopy, several sources and sinks combine together to determine the average ammonia concentration in the )
-1 800 canopy, including exchange with plant tissues through stomata, s
-2 with leaf cuticles and with decomposing leaf litter and the soil 600 surface (e.g. Denmead et al., 1976; Sutton et al., 1993b, 1998). (ng m
3 Ammonia within or immediately above the canopy airspace may 400 undergo chemical reactions, for example forming particulate matter, while depletion of gases within a plant canopy coupled Flux NH Flux 200 with altered microclimate can lead to evaporation of ammonium containing aerosol (Brost et al., 1988; Nemitz et al., 2004b). 0 Finally, the complex nature of ammonia sources and sinks in rural
-200 landscapes means that strong horizontal gradients of ammonia 29/04/2005 29/04/2005 30/04/2005 30/04/2005 30/04/2005 30/04/2005 occur. The result is that ammonia is not simply deposited from 14:24 20:24 02:24 08:24 14:24 20:24 above, but fluxes are often significantly influenced by advection Date, Time (GMT) effects, for example where advection from a ground level source beneath a micrometeorological reference height adds substan- Fig. 3.2. Fluxes of NH3 measured by eddy covariance over intensively managed grassland (Easter Bush, Scotland) several days after the application of liquid manure to tially to a net deposition flux (Loubet et al., 2001, 2008, 2009; the grassland (Sutton et al., 2007). Milford et al., 2001b). It is relevant to summarize the main influences on the primary drivers of exchange, the atmospheric ammonia concentration and the mean concentration of ammonia within the canopy. The first (Shaw et al., 1998; Famulari et al., 2004; Whitehead et al., 2008). of these is influenced partly by dispersion from adjacent In principle, reliable flux measurements can now be made for ammonia sources and partly by exchange with the surface itself. periods of large ammonia fluxes (e.g. after manured application), Over a surface which acts as an ammonia sink, above-canopy as has recently been demonstrated in an inter-comparison of ammonia concentrations are depleted relative to background two laser systems (Fig. 3.2). However, there was little correlation concentrations, while above canopy concentrations may be for fluxes <50 ng m 2 s 1, while the AMANDA systems have been 2 1 significantly enhanced if the surface is a net source (e.g. Sutton shown to be able to measure <10 ng m s (e.g. Sutton et al., et al., 2000a). 1998). Table 1 provides an overview of these and other systems The mean ammonia concentration of the canopy itself results for measuring ammonia fluxes. In principle, TDL and EC has the from the resolution of competing emission and deposition potential to be rated as high as the continuous gradient methods, processes with leaf cuticles, through stomata and with the ground but this still needs to be demonstrated by a more substantial body surface. The concept of ‘compensation point’ concentrations has of published measurements, particularly over longer time often been used to describe these relationships. The earliest view of periods and of a suitable quality for testing of models. a compensation point for ammonia related it to exchange through plant stomata with the leaf apoplast (Lemon and Van Houtte, 1980; 3.3. Key controls on biosphere atmosphere exchange of ammonia Farquhar et al., 1980). Under this interpretation, net ammonia fluxes would depend on the difference between what has since Fig. 3.3 summarizes the main processes affecting the net been termed the ‘stomatal compensation point’ (cs) and the exchange of ammonia with the atmosphere (Sutton et al., 2007). atmospheric concentration (ca). By contrast, subsequent studies The primary driver of ammonia exchange is the difference highlighted the fact that ammonia deposition rates were often between the atmospheric ammonia concentration and the faster than feasible by stomatal uptake, demonstrating the impor- tance of ammonia deposition to leaf cuticles (e.g. Sutton et al., 1993a,b; Duyzer, 1994). The resolution between these positions was provided in the development of the concept of the ‘canopy compensation point’ (cc), which accounts for both bi-directional stomatal exchange and deposition to leaf cuticles (Sutton and Fowler, 1993; Sutton et al., 1995). Such canopy compensation point concepts have since been further developed to include bi-direc- tional exchange with leaf surfaces and exchange with the ground surface under the canopy (e.g. Sutton et al., 1998; Flechard et al., 1999; Nemitz et al., 2001b). These inter-relationships are developed quantitatively in a two- layer canopy compensation point model (Nemitz et al., 2001b). One of the key points to note about ammonia compensation points is that they depend on the net solubility of ammonia in aqueous solution, which is largely dependent on its equililbrium with ammonium ions. By combining the temperature dependence of the Henry equilibrium and the ammonium dissociation equilibrium, Fig. 3.3. Summary of the key issues affecting the net land–atmosphere exchange of the gaseous ammonia concentration can be compared with a given ammonia. Each of these interactions can lead to ammonia fluxes changing with height þ þ ratio of [NH4 ]/[H ], which has been termed G (Nemitz et al., 2000c, above the ground. Ideally, flux measurements, based on e.g. relaxed eddy accumulation 2001a; Sutton et al., 2000b). On this basis, G can be used to provide or eddy covariance, made at several heights above the canopy would be used to quantify these effects, though until now such assessments have had to focus on the use temperature-normalized compensation points, for example þ þ of vertical profiles in mean ammonia concentration. cs ¼ f(T, Gs) where Gs ¼ [NH4 ]apoplast/[H ]apoplast. D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5209
3.4. Effects of ecosystem type on ammonia It is feasible that these emission periods represent events biosphere–atmosphere exchange of desorption of previously deposited ammonia occurring in dry conditions. Conversely, it is also feasible that they represent It has long been established that ecosystem type affects net apparent ‘emissions’, being an artefact whereby horizontal ammonia fluxes (cf., Denmead et al., 1976; Horva´th, 1983; Sutton ammonia concentration gradients away from an adjacent ground- et al., 1993b, 1995; Duyzer, 1994). Overall, unfertilized ecosystems, based ammonia source (e.g. manure spreading, farms etc) lead to such as forest and moorlands are generally sinks for atmospheric an advection error. This would reduce the measured deposition ammonia, while fertilized and grazed agricultural ecosystems tend rate and could explain apparent ammonia upward fluxes in this to show bi-directional fluxes with some periods of deposition and context. This illustration emphasizes the complexity of measuring some periods of emission. Of course, the distinction is not absolute, ammonia exchange processes and highlights the need for further as smaller ammonia emissions may also occur from semi-natural investigation of each option. ecosystems (e.g. Sutton et al., 1995; Flechard and Fowler, 1998). The above example of mainly ammonia deposition to a forest However, such a general difference is clear, and can be explained by ecosystem may be contrasted with recently published measure- the increase in cs and cground that occurs in fertilized and grazed ments of ammonia fluxes over an intensively managed grassland in ecosystems. Two recently published examples of ammonia exchange the Netherlands (Wichink Kruit et al., 2007). The diurnal patterns in provide a useful basis to highlight these differences. ammonia concentration and net exchange flux are illustrated Neirynck et al. (2005) report ammonia flux measurements made in Fig. 3.5. Hourly ammonia concentrations in the air at this site using the AMANDA technique (Wyers et al., 1993) over a coniferous were again very large, 1–50 mgm 3, with an overall mean of around forest in Belgium. Their forest site occurs in an area of intensive 10 mgm 3. In this case, net emission occurred for around 40% of the livestock rearing, so that ammonia concentrations from some wind diurnal period (10:00–20:00), with net deposition at other times. directions are very large (5–25 mgm 3) while for other wind sectors Wichink Kruit et al. (2007) also estimated the canopy 3 0 ammonia concentrations were more moderate (2–4 mgm ). Even compensation point (cc) based on profile estimation of c (zo ). They considering the effects of canopy wetness, in all conditions the then combined this with estimates of surface temperature to esti- 0 mean diurnal profiles show consistent net deposition to the forest mate G(zo )or‘Gc’ from the measurements (Fig. 3.6). Estimated 3 canopy. Curiously, the largest deposition fluxes occurred in dry values of cc were in the range 1–30 mgm , which is comparable conditions, which is unusual, as Rw would be expected to be smaller with other studies for managed grassland (e.g. Milford et al., 2001a; when the canopy is wet (Sutton et al., 1995; 1998; Nemitz et al., Sutton et al., 2001; Loubet et al., 2006), and substantially smaller 2001b). Although this difference is partly explained by different than the upper values implied for the forest in Fig. 3.4. Normalized values of Fmax during conditions of different canopy wetness, this for canopy temperature, the values of Gc were in the range does not to fully explain the difference. Further analysis by Neir- 200–11,000 through a period of May to October 2004, with a mean ynck et al. (2005) showed differences in the overall canopy resis- value of just over 2000. tance (Rc) for ammonia deposition with different canopy wetness and temperature, and with larger values of Rc occurring at higher 3.5. Modelling surface–atmosphere exchange of ammonia ammonia concentrations. Neirynck et al. (2005) did report periods of net ammonia Over recent years, the canopy compensation point approach has emission from their forest canopy (Fig. 3.4). These were recorded become the standard technique to model bi-directional ammonia during periods with winds from the high ammonia wind sector and surface–atmosphere exchange. Starting with the 1-layer models found to only happen at very large ammonia concentrations, which offsetting bi-directional stomatal exchange against deposition to occurred when air temperatures were larger than 15 C and relative leaf surfaces (Sutton and Fowler, 1993; Sutton et al., 1995, 2007), humidity less than 60%. Fig. 3.4 presents an intriguing result, since subsequent models have developed in several directions. The main according to the concepts of ammonia compensation points subsequent developments can be summarized as follows: a different picture should emerge, namely that periods of ammonia Treatment of multiple canopy layers: In addition to ammonia emission occur when atmospheric ammonia concentrations are exchange with the top part of the canopy, leaf litter and the soil small. By contrast, such a relationship is possible when emissions surface have been shown to be important sources of ammonia from a canopy are strong (and not compensation point driven), so emission into the plant canopy (e.g. Nemitz et al., 2000a). For an that it is the emissions from the surface that generate increased oilseed rape canopy Nemitz et al. (2000b) also highlighted the ammonia air concentrations. This phenomenon was also observed importance of an upper and lower part of the main foliage, dis- following harvest of an oilseed rape field (Sutton et al., 2000a,b). tinguishing the main foliage from an over canopy of oilseed ‘siliques’. In practice this three layer model becomes complex to parametrize, and there has now developed consensus that a two-layer model represents an appropriate balance of realistic description while avoiding excess complexity. A recent implementation of the 2-layer model is that of Personne et al. (2009) for the GRAMINAE Integrated Experiment (Sutton et al., 2009b). They used measured bioassay estimates of Gs and Glitter (Mattsson et al., 2008a,b; Herrmann et al., 2008) combined with an energy balance approach to calculate component resistances, showing close agreement with measured ammonia fluxes (Fig. 3.7). Treatment of cuticular fluxes: The initial parametrisations of the cuticular resistance (Rw) allowed only for deposition, depen- dent on relative humidity (Sutton and Fowler, 1993; Sutton et al., 1995) or vapour pressure deficit (Nemitz et al., 2000c, 2001b). As noted above for the forest example, ammonia deposited to a canopy Fig. 3.4. Dependence of ammonia flux on concentration in the high ammonia wind surface may also be re-emitted to the atmosphere, particularly sector during warm daytime conditions with dry canopy (Neirynck et al., 2005). under drying conditions. A first approach to simulate this effect 5210 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
Fig. 3.5. NH3-concentration (upper panel) and NH3-flux (lower panel) measurements above managed grassland in The Netherlands from 18 July until 15 August 2004 (summer 3 2 1 period). The horizontal axis represents time of the day (UTC). Local time is UTCþ2. The vertical axis represents the NH3-concentration (mgm )orNH3-flux (ng m s ). Diamonds A are calculated values for the half-hourly NH3-concentration or NH3-flux; the solid line (dd ) (with vertical 25 and 75 percentile bars) is the median of all half-hourly fluxes for that time. The dashed line (- - -) in the lower panel is the mean leaf wetness signal during this period (Wichink Kruit et al., 2007).
treated the leaf surface as a humidity dependent capacitance 3.6. Dynamic simulation of ecosystem C–N cycling (Qd), which would be in equilibrium with a non-zero leaf surface and ammonia fluxes concentration (cd)(Sutton et al., 1998). In this case an adsorption/ desorption resistance (Rd) is also defined. This first dynamic A disadvantage of the compensation point scheme for simu- approach had the advantage of being able to simulate ammonia lating ammonia fluxes outlined above is that empirical values of charging and discharging of the cuticle, but had the disadvantage G must be provided. The only way forward from this position is to that the leaf surface pH needed to be specified as an input. The develop models of carbon–nitrogen cycling that can simulate approach was further developed by Flechard et al. (1999) who G values for the different pools based on an understanding of the considered the full aqueous chemistry on leaf surfaces, dependent pool dynamics (cf. Massad et al., 2008). To date, the only such on multiple air pollutant inputs and potential leaching of base model to attempt this coupling is the PaSim model of Riedo et al. cations from leaf surfaces. In this model, leaf surface pH is solved by (2002). The model distinguished plant nitrogen pools into struc- ion balance, and the model is able to take account of the effects of tural nitrogen, substrate nitrogen and apoplastic nitrogen (a sub- other trace components such as SO2 on ammonia fluxes. Burkhardt pool of substrate nitrogen), linking these with plant uptake and et al. (2008) have recently extended this model to incorporate the growth processes. The model was parametrised based on measured two-layer approach with bi-directional exchange for each of the fluxes for a Scottish grassland (Milford et al., 2001b) and has leaf surface, stomata and ground surface. The cuticular resistance recently been tested for the Braunschweig Experiment (Sutton clearly responds to the chemistry of the liquid film on vegetation et al., 2009b). Overall, the model was able to simulate the larger and the combination of reactive gases present (Flechard et al., net emissions that occurred after cutting and after fertilization, as 1999). Even in the absence of additional reactive trace gases, the well as the decline in the 10 day period following fertilization. By cuticular resistance declines with increasing NH3 concentration. contrast, the component fluxes were less well described. Bioassays, In a series of chamber experiments Jones et al. (2007) quantified chamber measurements and within-canopy profiles during the the relationships between ambient NH3 concentration and the bulk Braunschweig Experiment (Mattsson et al., 2008a,b; Herrmann canopy resistance for a range of moorland vegetation as shown in et al., 2008; Nemitz et al., 2009b) highlighted leaf litter as being Fig. 3.8 in which the non-stomatal ‘cuticular’ resistance is seen to a key source of emission following cutting. This source is currently increase lineary with NH3 concentration, leading to much smaller not simulated in PaSim, which simulated that increased emissions deposition at high concentrations than if deposition velocity after cutting were due to an increase in apoplastic ammonium. The remained constant with concentration as is usually assumed. bioassays indicated that the foliage was more likely to be a sink of D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5211
0 þ þ 0 Fig. 3.6. Derived canopy compensation points (cc ¼ c(zo )) (upper panel) and ratios between NH4 and H concentration (Gc ¼ G(zo )) (lower panel) from the end of May until the end of October 2004 (diamonds) and the constant value (2200) that is normally assumed for modeling (line) (Wichink Kruit et al., 2007).
soil/litter ammonia emissions, highlighting the need for improved 1. That unreplicated ammonia flux measurements in most studies ecosystem modelling of ammonia exchange that accounts for litter are highly uncertain and need to be considered with caution decomposition processes (Sutton et al., 2000b). when compared with model estimates. 2. That advection effects can significantly influence measured 3.7. Integrating ammonia exchange processes ammonia fluxes, both due to dispersion away from nearby point sources (correction for advection effects increases net The preceding sections have highlighted the many processes deposition) and due to emissions from an emitting field itself and interactions that combine to regulate ammonia fluxes between (correction for advection increases net emission). vegetation and the atmosphere. It thus becomes a major challenge 3. That gas-particle interactions had a minor effect on measured to integrate each of these processes to develop a holistic view. It is ammonia fluxes, though the relative effect on calculated aerosol necessary to quantify the interactions in each case in order that deposition rates was significant (being the cause of apparent valid conclusions can be obtained. This creates a major challenge aerosol emissions). for experimentalists to be able to address all the questions in the 4. That reasonable agreement can be made between relaxed eddy field. For example, in the absence of measurements of horizontal accumulation for ammonia and the aerodynamic gradient concentration profiles, it is difficult to quantify the potential for method, though measurements are not yet sufficiently precise advection effects to have influenced the results presented in to detect flux divergence (except for possible cases of extreme Fig. 3.4. Similarly, it remains an open question in most studies advection errors). whether gas-particle interactions have a significant influence on 5. That net emissions from this grassland canopy are controlled by measured ammonia fluxes. the recapture of leaf litter ammonia emissions by overlying It was with such interactions in mind that the GRAMINAE foliage and the interaction of cuticular exchange pools with Integrated Experiment was designed. A number of findings from mainly stomatal uptake of ammonia from the leaf litter emis- this experiment have already been mentioned, but the experiment sions. Net emissions increase following cutting due to exposure demonstrates both the challenges and the power of developing an of the litter and cutting induced senescence, with a similar integrated approach. Fig. 3.9 summarizes each of the issues and recapture process affecting net emission following fertilization. measurement methods that were investigated during this experi- ment (Sutton et al., 2008a). The key conclusions from this study A range of models is able to simulate the dynamics of net have been summarized by Sutton et al. (2009b) and include: ammonia exchange with the managed grassland, but further 5212 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267
4500 key challenges include the climate dependence of net ammonia Measured flux 4000 emission and deposition, and the characteristic fluxes of other Modelled flux (Gamma Litter) ecosystems in the world.
) 3500 -1
s In principle the models of ammonia exchange incorporate the 3000 -2 main features of environmental conditions and could therefore be 2500 applied in different climates. Here the limitations include the lack of 2000 available data for empirical factors such as G values and the back- 1500 ground data to extrapolate to conditions with different climates. fluxes (ng m fluxes (ng 3 1000 Currently, the estimates of G have mainly been derived for cool NH 500 European conditions and for a very limited number of ecosystems. 0 Although there have been many studies of ammonia emission from -500 fertilized tropical systems, such as rice and maize, there are few 22-May 29-May 05-Jun 12-Jun published studies of ammonia fluxes over semi-natural unfertilized tropical ecosystems. The rates of ammonia deposition or emission Fig. 3.7. Comparison of ammonia fluxes simulated by a two-layer canopy compensa- in these situations are thus highly uncertain. Given the differences tion point model (SURFATM-NH3) with measured fluxes (Fmg) during the GRAMINAE Braunschweig Experiment. For this model scenario, the ground emission is assumed to in biology of these systems, measurements are required to underpin originate from leaf litter based on measured Glitter (Personne et al., 2009). modelling approaches. A modest degree of climate change (e.g. þ2 C) is a much easier matter to simulate, for example based on the analysis of temperature attention is needed to develop dynamic treatments of ammonia effects within existing datasets. The thermodynamics of ammonia emissions from leaf litter decomposition. solubility and dissociation are rather straightforward, indicating for example a doubling in cs every 5 C increase for a given value of 3.8. Future challenges for ammonia exchange G (Sutton et al., 2001). However, caution is needed before making climate change simulations on this basis. Analysis of the PaSim model The results from the GRAMINAE Integrated Experiment provide under different temperature regimes showed that net ammonia a microcosm of some of the key challenges to measure ammonia fluxes for Easter Bush in Scotland (cf. to measured fluxes of Milford fluxes and model the process interactions. In a wider perspective et al., 2001b) were quite insensitive to temperature. For example, increased temperature (in the absence of moisture limitation) led to increased grass growth which diluted available nitrogen pools, thereby reducing G values (Sutton and Milford, unpublished simu- lations). Similarly, increases in wetness, while favouring smaller values of Rw may also lead to increased rates of leaf litter decompo- sition, favouring ammonia emissions. To take another example, in colder conditions, NH3 from manure application to the land surface tends to be emitted at smaller rates, but the emission lasts longer, especially if a waterlogged or frozen soil conditions prevent infiltra- tion. With these illustrations in mind, it becomes a major future challenge to generalize how ammonia fluxes might change in the future under different climatic regimes.
4. Sulphur dioxide
4.1. Introduction
There are three very different spatial scales relevant to the exchange of SO2 at terrestrial surfaces, first the micro-scale, at which the chemical and biological interactions occur (Fig. 1.1). Second is the spatial scale at which most of measurement and interpretation takes place, which is the field scale (averaged over 103–105 m2)for measurements using micrometeorological methods. Lastly, the application of knowledge of the surface exchange process is primarily at regional to continental scales to characterise the fluxes and budgets within chemistry transport models (CTM) and comparisons with the concentration fields observable from satellites. The measurements have mainly been made at the field scale using micrometeorological methods, although there have been some laboratory studies, mainly in the early days of SO2 dry deposition research. The initial measurements were used to esti- mate the regional scales of dry deposition, often using a fixed deposition velocity as a key variable within long-range transport Fig. 3.8. The relationship between ambient ammonia concentrations and the cuticular models (e.g. Fisher et al., 1978). With larger datasets of measure- resistance to deposition for moorland vegetation. A significant difference was found ments covering a wide range of conditions, it is clear that rates between day and night for the bulk canopy resistance (R ), which included both c of dry deposition vary considerably in time and space (Fowler stomatal uptake and deposition to the leaf surfaces. Once the effect of the stomatal resistance (Rs) was accounted for, the cuticular resistance (Rw) was found to approx- and Unsworth, 1979) in particular because the sinks available at imately constant day and night (Jones et al., 2007). terrestrial surfaces, including the stomata in vegetation, leaf D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5213
Estimation of Advection of Quantification of Interactions with acid
farm-scale NH3 NH3 from nearby energy balance & gases and ammonium emissions sources & effects environ controls on particles & effects
from plume on vertical fluxes NH3 exchange on net NH3 fluxes measurements
Continuous
measurement of NH3 fluxes by gradient Effects of dew and REA approaches Effects of cutting & leaf surface & N fertilization chemistry events & choices on NH fluxes 3 NH3 compensation points of foliage Plant bioassay Within-canopy determination of cycling of NH3 fluxes NH3 emission potential
Determination of Effects of leaf senescence
within-canopy NH3 release from and plant species on turbulent exchange litter decomposition NH3 emission potential
Soil chemistry interactions
with plant N uptake & NH3 fluxes
Fig. 3.9. Overview of issues addressed by the GRAMINAE Integrated Experiment (Sutton et al., 2008a).
surfaces and the presence of liquid water on vegetation from dew 4.2. Worldwide advances in SO2 flux monitoring and modelling or rain, present a variable absorbing surface. The data have shown the role of atmospheric composition and surface leaf water chem- 4.2.1. Asia istry in controlling canopy resistance. Sulphur dioxide dry deposition to vegetated surfaces is largely Most dry deposition measurements of sulphur dioxide over the controlled by non-stomatal processes, but in many arid ecosystems last 30 years have been made in N. America and Europe, and have and deserts of the world where vegetation is sparse, the nature and served as a basis for the parameterisation of dry deposition models pH of soils determine the sink strength. In Asia, substantial efforts (Erisman,1994; Smith et al., 2000; Zhang et al., 2002), which in turn have for example gone into the characterization of SO2 uptake by have been applied to ecosystems in different parts of the globe. loess soils, given their large geographical representation in However, most SO2 emission and deposition now occurs outside Northern China, their alkaline nature and their ability to neutralize N. America and Europe. Asia’s contribution in 1985 of 20% to global atmospheric acidity and to serve as an oxidation medium for anthropogenic SO2 emissions has doubled since then, reaching 37% SO2. Both micrometeorological and laboratory- or field-based by the year 2000, of which 23% is emitted by China alone and 5% by flow reactor methods were deployed. New micrometeorological India. measurements over forests and short vegetation have also been Southern China is one of the world’s most sulphur polluted reported over the last 10 years in the region, reflecting the growing areas. Paradoxically, in Northern China, where ambient SO2 concern over increasing sulphur emissions and deposition to concentrations are very large, rainfall is generally alkaline, and the ecosystems. areas polluted by acid rain do not necessarily correspond to the areas of high SO2 emissions. One of the reasons for this discrepancy 4.2.1.1. Sulphur dioxide deposition to soils. Utiyama et al. (2005) is the presence of alkaline soils (yellow sand) distributed over measured dry deposition to loess soil and dead grass in Beijing the arid areas of N.W. China (e.g. the loess plateau and Gobi desert), using the aerodynamic gradient method, though in neutral condi- the windborne erosion of particles with high base cation concen- tions 22% of the time. In stable or unstable thermal stratification, trations can neutralize atmospheric acidity (Utiyama et al., 2005). they used a surface reaction concept for inferring dry deposition. Loess soil, which covers vast areas of the Eurasian continent Two surface kinetics models were considered: either i) the reaction extending from N.E. China to Central Asia, contains Ca in large occurs in soil pores and SO2 molecules diffuse through porosity quantities, and calcium carbonate (CaCO3) reacts with atmospheric while reacting with alkaline sites on the pore surface; or ii) the SO2, to form calcium sulphate (CaSO4). Thus, even bare soil without adsorption mechanism is of Langmuir–Hinshelwood type, where vegetation may be a significant sink (Sorimachi and Sakamoto, the partial pressure of SO2 and its desorption pressure from the site 2007), which may affect the regional SO2 budget if the process is are in equilibrium. The model parameters are then fitted so that the inadequately quantified in dry deposition models. resulting (modelled) vertical SO2 concentration gradient matches In this section we review research and monitoring from the last the observations. Measured deposition velocities (Vd) were in the 1 decade, including SO2 dry deposition measurements from Asia, range 1–12 mm s . North America and Europe, as well as findings from long-term flux Sorimachi and Sakamoto (2007) conducted laboratory-based monitoring experiments. The current state of knowledge concern- flow-reactor measurements of SO2 deposition to soil samples from 12 ing mechanisms of SO2 dry removal from the atmosphere is sites in the arid loess plateau and deserts of Northern China. Canopy reviewed, with consequences for temporal trends in atmospheric resistances in the range 28–650 s m 1 (with a mean of around concentration and deposition, and key future research areas are 200 s m 1) were found to be dependent on RH, as was S(IV) oxidation identified. to S(VI). It was hypothesized that Northern China soils, which are 5214 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 much more alkaline than in Southern China, are a greater sink for SO2 The latter can only provide crude estimates of deposition rates, as and a neutralizing buffer for acidifying atmospheric deposition surrogate surfaces do not adequately account for the complexity of (Sorimachi et al., 2004). By comparison, in modelling SO2 deposition natural surfaces, but they do allow continuous monitoring at to Asia, Xu and Carmichael (1998) used a fixed Rc for deserts of a number of sites and help to detect trends. The use of inferential 500 s m 1, which is clearly too high in the case of Northern China methods requires underpinning measurements of surface resistance deserts. A flow-reactor was also used by Sakamoto et al. (2004) to to representative surfaces to parameterise the models used. It is also determine SO2 dry deposition to yellow sand and soil-mediated SO2 important to note that long-term changes in canopy resistance are oxidation by O3. The deposition velocity for SO2 increased with RH likely, especially in regions in which the relative ambient concen- due to the positive effect of RH on the SO2 oxidation rate. trations of SO2 and NH3 change with time, as for example in Europe during the period 1990–2005 (Fowler et al., 2001c, 2007). 4.2.1.2. Micrometeorological measurements over vegetated Ta et al. (2005) thus provided long-term sulphur dioxide dry areas. Matsuda et al. (2006) reported micrometeorological (aero- deposition estimates across Gansu Province, China, using K2CO3- dynamic gradient) flux measurements of SO2 and O3 over a tropical coated surrogate sulfation plates. Samples were taken monthly for (teak) forest in Northern Thailand in dry and wet seasons. The 11 years at 48 sites distributed across 11 cities in the province. deposition velocity for SO2 in the dry season was rather low, in the The data showed that cumulative SO2 dry deposition fluxes were 1 1 range 1–3.1 mm s in daytime and 0.8–1.1 mm s in night-time. closely related to local SO2 emissions, and had seasonal variations In the wet season, however, Vd was much higher due to enhanced with maxima in winter and minima during summer as a result of non-stomatal uptake in wet conditions, with values in the range higher winter and lower summer SO2 emissions and concentra- 1 1 9.5–13.9 mm s in daytime and 2.6–4.2 mm s in night-time. The tions. Monthly average SO2 deposition velocities, however, peaked data were compared with a recent non-stomatal resistance scheme in April–July at 11–27 mm s 1, and minimum values were observed (Zhang et al., 2003a), and it was concluded that extended experi- in January at 2–10 mm s 1. mental SO2 dry deposition studies are needed in the tropics, while Inferential models (Erisman, 1994; Smith et al., 2000; Zhang Zhang et al. (2003a) recommend more studies to quantify the et al., 2002) may be used to estimate dry deposition at observation different effects of dew and rain on SO2 deposition. sites, where single-height ambient concentration measurements Sulphur dioxide dry deposition was also measured by Matsuda are available together with standard meteorological data. Model et al. (2002) over a red pine forest located in Oshiba Highland, parameters, however, have been largely derived from European and Nagano, Japan, using a Bowen ratio technique. The median daytime N. American studies and may not necessarily be adequate for Asian (12:00 to 14:00) deposition velocity was 9 mm s 1. Measurements vegetation and soils, and numerical evaluations need to be carried compared favourably with estimates by an inferential model for out. Thus Takahashi et al. (2002) simulated the dry deposition of wet conditions, but for dry or mixed wet-dry surfaces there were SO2 to a Japanese cedar (Cryptomeria japonica) forest located in large differences between model and measurements. The authors Gumma Prefecture, based on the results of 1-year’s concentration 1 ascribed the discrepancy to a relative humidity threshold value measurements. The mean modelled Vd at this site was 8.8 mm s used in the inferential scheme to characterise canopy wetness, and (Takahashi et al., 2001). The inferential estimate of the dry sulphur pointed to the need for a refined parameterisation of the cuticle or deposition flux was 11.1 mmol m 2 yr 1 (3.6 kg S ha 1 yr 1), which external leaf surface resistance. compared well with the net throughfall flux (12.4 mmol m 2 yr 1, 1 1 In a study of SO2 and O3 dry deposition to short grassy vegeta- or 4.0 kg S ha yr ). Over a broadleaf forest on typical red soil of tion over an alkaline soil (pH ¼ 9.2) near Beijing, using the Southern China, Xu et al. (2004) simulated Vd for SO2 and partic- 2 aerodynamic gradient method, Sorimachi et al. (2003) measured ulate SO4 , as well as their atmospheric deposition fluxes. The 1 1 mean Vd values of 2 ( 1) mm s and 4 ( 2) mm s in late summer simulations indicated that about 99% of the dry sulphur deposition and early winter, respectively. Although the grass was lush and flux in the forest resulted from SO2, which contributed over 69% of thick in the late summer, and senescent and leafless in the early the total (wet þ dry) annual sulphur deposition. winter observation period, there was no difference in the mean Rc By comparison, Wang et al. (2003) computed dry deposition 1 1 2 (180 270 s m and 180 300 s m , respectively), but the fluxes of SO2 and SO4 for 1 year to agricultural land over red soil uncertainties given reveal a large variability in measured Rc. The (pH ¼ 5.3–5.8) in the Jiangxi province of Central China. The crops difference in Vd stemmed from the higher aerodynamic (Ra) and grown were rice paddies and oilseed rape. Sulphur dioxide 1 1 quasi-laminar sub-layer (Rb) resistances in late summer than in concentrations were measured 8 times day , 7 days month , using early winter. The absence of vegetation and stomatal uptake in a bubbler method. Annual mean modelled estimates of Vd were 1 1 2 early winter, which might otherwise have reduced the SO2 sink 3.8 ( 0.16)mms for SO2 and 0.20 ( 0.12) mm s for SO4 . strength, seems to have been compensated for by the soil alkalinity. Measured monthly mean concentrations ranged from 9 to As the soil was more exposed and the in-canopy aerodynamic 163 mgSm 3 (6.7–121 ppb), with an annual mean of 64 mgSm 3 resistance was reduced, the soil surface offered more adsorption (47 ppb). Estimates of total monthly wet and dry deposition of SO2 2 1 and reaction sites for SO2, with the result that the field was an and SO4 ranged from 2.2 to 20.3 kg S ha with an annual total 1 equally efficient SO2 sink in early winter as in summer. deposition of over 100 kg S ha , of which 83% was via dry deposi- The deposition velocity for SO2 was measured by Jitto et al. tion, accounting for over 90% of total S input to farmland in this area. (2007) during a 1-year experiment over a canopy of irrigated rice paddy in Thailand using the Bowen ratio technique. The deposition 4.2.2. North America velocity was highest around noon and lowest at night. Seasonally- Few long-term datasets of SO2 dry deposition monitoring have 1 averaged values of Vd were 6.7, 12.5, and 15.1 mm s in the winter, emerged over the last 10 years, reflecting the declining importance of summer, and rainy seasons, respectively. SO2 as an acidifying input relative to NOy and NHx. Advances have nonetheless been made in inferential modelling of SO2 uptake, 4.2.1.3. Long-term deposition studies and inferential modelling. As especially regarding the quantification of the non-stomatal (external) alternatives to costly and labour-intensive micrometeorological leaf surface resistance, which serve as a basis for simulating regional measurements of dry deposition, several authors in Asia have esti- patterns of SO2 deposition (Zhang et al., 2002, 2003b). mated long-term SO2 deposition using monitored concentration data Micrometeorological SO2 flux data from 5 sites (2 forests, a corn and inferential models, or long-term artificial collection devices. field, a soybean field and a pasture) in eastern USA (Finkelstein et al., D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5215
2000; Meyers et al., 1998) were compared with modelled data by that the flux remained nearly constant. The measurements of the Zhang et al. (2003a), with the specific objective of evaluating the new atmospheric terms (Ra and Rb) revealed no trend, thus the trend in non-stomatal resistance scheme of the new Canadian model (Zhang Rc is not caused by changes in turbulence, and are clearly a conse- et al., 2002, 2003b). Over the forest sites, Finkelstein et al. (2000) had quence of the chemical affinity of the surface changing with time. noted that wetness tended to increase deposition velocity, but that These dry deposition measurements have proved valuable in the nature of wetness (rain or dew) and its chemistry also controlled explaining the consistently larger decline in ambient SO2 concen- canopy resistance. Non-stomatal surfaces like leaf surface, stem, tration than in emissions in Europe. In the absence of these flux trunk and ground were important sinks for SO2,andtheauthors measurements it would be a matter of speculation as to the under- concluded that a better understanding of surface chemistry and lying cause of the faster decline in ambient concentration than water film chemistry was needed. emission. Even with these measurements there remains the possi- Dew formation has long been recognised as an important sink for bility that SO2 oxidation rates have increased due to the growing SO2 (Fowler and Unsworth, 1974, 1979). In more recent work Meyers oxidizing capacity of the atmosphere and have contributed to the et al. (1998) show that dew is the reason for the relatively high early relative changes in emission and deposition (non-linearity). It will be morning deposition rates at 2 of the 3 low vegetation sites studied in necessary in the further analysis and interpretation of European Eastern USA. Recognizing the weakness of existing North American pollution climate data to carefully examine the relative importance parameterisations (e.g. Meyers et al., 1998) in predicting SO2 depo- of the different contributors to the observed trends in concentration sition rates to non-stomatal surfaces, especially in wet canopies, and deposition and quantify the relative importance of changes in Zhang et al. (2003a) demonstrate that the AURAMS scheme (Zhang dry deposition and oxidation rates in the long-term trends. et al., 2002) performed well at these 5 sites, using different resistance values for dew and rain. The revised non-stomatal resistance scheme 4.2.3.2. Other recent European datasets. The SO2 flux–gradient data (Zhang et al., 2003b) includes a treatment of in-canopy transport, soil obtained over short vegetation by Feliciano et al. (2001), collected and cuticle terms, and is a function of relative humidity, leaf area over a period of 3 years in the mid 1990s at 3 different sites in index and friction velocity. For wet canopies, the cuticular resistance Portugal, were important in providing Rc estimates for the Medi- is treated differently for dew and rain. terranean region of Southern Europe. The 3 sites had contrasting pollution climates, with a coastal, oceanic, humid meadow in 4.2.3. Europe N. Portugal, a hot and semi-arid pseudo-steppe and a site located in 4.2.3.1. Long-term flux monitoring in the UK. Sulphur dioxide fluxes a mostly dry, intensive agricultural area, both in S. Portugal. Median have been monitored continuously since the mid 90s at two rural canopy resistances varied from 140 s m 1 to 200 s m 1 and although sites in the UK, over agricultural land at Sutton Bonnington in the stomatal uptake was important when vegetation was biologically English Midlands, and over moorland at Auchencorth Moss in active, the annual deposition was dominated by non-stomatal S. Scotland (Fowler et al., 2001c, 2005, 2007). The dry deposition mechanisms on wet surfaces. The night-time canopy resistance, measurements have continued to bring surprises over the last 10 a proxy for the non-stomatal resistance, increased with decreasing years. At Sutton Bonnington, the ambient concentrations have relative humidity at all 3 sites. A comparison of nocturnal Rc for the declined from about 2.8 ppb in 1996 to current values close to southern sites showed that, for a given level of relative humidity, the 1.4 ppb and yet the deposition velocity continues to increase due to Rc at the intensive agricultural site was systematically lower than at continued reduction in the canopy resistance (Rc)(Fig. 4.1). Over the pseudo-steppe site, which is used more extensively for grazing the monitoring period the canopy resistance has almost halved and and hay production. Although the authors make no mention of NH3 is now about 70 s m 1. The consequence of the steady decline in being measured at these sites, it might be hypothesized that a higher canopy resistance along with a decline in ambient concentration is NH3 concentration at the intensive agricultural site may have been
10 9 8 7 6
5 NHNH33 SO 4 SO22 3
Mixing ratio (ppb) 2 1 0 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 )
-1 0.6 140
0.5 120 R SO2/NH3SO2/NH3
100 c c2 SO 0.4 Rc (Wheat) 80 2
0.3 m (s 60
0.2 -1 molar ratio (ppb ppb
40 ) 3 0.1 20 /NH 2 0 0 SO 1995 1996 19971998 1999 2000 2001 2002 2003 2004
Fig. 4.1. Changes in the mean concentrations (ppbV) and ratio of ammonia and sulphur dioxide and in the May–July canopy resistance for SO2 deposition on Wheat at Sutton Bonnington between 1996 and 2003. 5216 D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 responsible for the observed lower nocturnal Rc, compared with the extensively-managed, steppe-like grassland. The development of low-cost systems for the long-term moni- toring of SO2,NH3 and other trace gas fluxes holds promise for widening the range of dry deposition datasets for comparison with inferential models. Hole et al. (2008) presentan18-monthdatasetof SO2 fluxes acquired with a conditional time-averaged gradient (COTAG) system (Fowler et al., 2001b; Famulari et al., in press)in a semi-alpine ecosystem in Southern Norway. The mean annual SO2 deposition velocity was 4.0 mm s 1, although the dataset included some negative deposition velocities (upward fluxes), and the annual 1 mean Vd was 13.0 mm s if only the positive values were included. The authors report evidence of enhanced SO2 deposition rates during anepisodeinNovember2005whentheNH3/SO2 ratio was high, and conversely of decreased SO2 uptake and increased NH3 uptake in November 2004 when the NH3/SO2 ratio was low. Comparison with the inferential model by Zhang et al. (2002, 2003b) was satisfactory but the model could not reproduce the large observed variability in exchange rates, which may result from NH3–SO2 co-deposition processes not being included in their resistance scheme. More experimental evidence of the mutual influences of NH3 and SO2 concentrations on their deposition rates was obtained by Derome et al. (2004), though not by micrometeorological measurements but using bulk precipitation collectors and through- fall measurements in Scots pine canopies in SW Finland. The study was conducted over a 6-year period (1993–1998) in the vicinity of a Cu–Ni smelter, which emitted large amounts of gaseous NH3.These emissions were shown to have strongly enhanced the scavenging of atmospheric SO2 by the pine canopy, resulting in increased levels of N and S deposition and increased foliar N and S concentrations. In an NH3 fumigation experiment, Cape et al. (1998) had previously described similar findings over a Scots pine forest in Central Scotland, with the canopy resistance for SO2 decreasing with elevated NH3 concentration. Although NH3 concentrations were not measured in Fig. 4.2. A schematic representation of the dynamic canopy compensation pollution the Finnish study (Derome et al., 2004), they were likely higher than model for SO2 and NH3 exchange over vegetation (from Flechard et al., 1999). normally encountered in the countryside, except near animal housing in areas of intensive agriculture, where such processes could 3 be significant. increased to 2 mgm to provide sufficient NH3 to neutralize the acidity from the ambient SO2 oxidation in solution. The conse- 4.3. Control of surface uptake rates by leaf cuticular chemistry quence of the increase in ambient NH3 is to decrease the canopy resistance for SO2 and increase the deposition rate of SO2 to the A number of authors have addressed the issue of the chemical maximum under the prevailing atmospheric conditions. Demon- control of surface pollutant uptake rates (e.g. Flechard et al., 1999). strating a close link between SO2 deposition and ambient NH3 is not The most important finding for SO2 deposition is that the rates of new, as this was predicted in earlier work in the Netherlands (Van deposition are controlled mainly by the chemistry at the vegeta- Hove et al., 1989). However, this work quantified the process tion–atmosphere interface, and that as the surfaces are wet most of the time, the processes are regulated by chemical processes within the thin film of moisture. In principle, many compounds influence 10 the chemistry of this surface layer, including plant exudates and soil )
2 0 derived compounds, but the key reactant for SO2 is NH3. Thus the SO -10 ambient concentrations of SO and NH essentially regulate the pH -1
2 3 s of the surface moisture and thus control the uptake of SO2. The -2 -20 full surface chemistry of the process has been incorporated into -30 Meas. FSO2 a dynamic mechanistic model shown in Fig. 4.2 (Flechard et al., χ Mod. FSO2 ( NH3 = ambient) -40 χ -3 1999). The chemistry of the surface water film is initialised in the flux (ng m Mod. F ( = 2 μg m 2 SO2 NH3 ) model using measured precipitation chemistry, the model then FSO2, max SO -50 simulates the dynamic responses of the net land–atmosphere -60 exchange of SO2 as the ambient concentrations of the reactive trace gases and meteorological conditions change. The model has been shown to provide good agreement with observed 30 min average fluxes for several days. An example is provided in Fig. 4.3, for a five-
day period at Auchencorth Moss in the Scottish Borders. The 21/03/95 12:00 22/03/95 00:00 22/03/95 12:00 23/03/95 00:00 23/03/95 12:00 24/03/95 00:00 24/03/95 12:00 25/03/95 00:00 25/03/95 12:00 26/03/95 00:00 26/03/95 12:00 general agreement between measured and modelled fluxes is Fig. 4.3. A comparison between measured and modelled SO2 fluxes at Auchencorth excellent during the three-day period, 21st March to 23rd March Moss over the period 21-3-95 to 26-3-95 showing the influence of increasing ambient 1995. From the 24th March, the observed NH3 concentrations are NH3 concentration on SO2 flux (from Flechard et al., 1999). D. Fowler et al. / Atmospheric Environment 43 (2009) 5193–5267 5217 correctly for the first time, demonstrated the effects in field period has meant that the SO2/NH3 ratio has decreased dramatically, conditions at ambient concentrations and provided a mechanistic resulting in a reduced Rc for SO2 (Fig. 4.1). The data for Fig. 4.1 are of model incorporating the full chemical scheme. course site specific and the footprint of the measurements in Not- tinghamshire is only 104–105 m2, but this trend may be regarded as 4.4. Advances in deposition modelling representative for the regions in which ambient SO2 concentrations have declined by up to an order of magnitude since 1970 and The magnitude of dry deposition at the national and regional includes much of central and eastern England and the industrial scales requires that process-based, rather than empirical, parame- regions of Germany, France, the Netherlands and Belgium. While terisations be implemented in atmospheric models, accounting for ambient SO2 was a relatively good proxy for total atmospheric and variations in surface chemical characteristics driven by local pollu- leaf surface acidity 15 or 25 years ago, the relative share of SO2 tion climates. The observation of changes in deposition velocity from compared to other inorganic atmospheric acids (e.g. HNO3 and HCl) the European SO2 deposition studies of the 90s is now widely known is now much smaller. The NitroEurope network of 56 DELTA and is being used by EMEP to explain growing discrepancies in the samplers across the European continent currently provides monthly model-measurement comparisons over Europe. The work has led to mean concentrations of HNO3 and HCl as well as SO2 and NH3 and þ 2- modifications of the EMEP model (Simpson et al., 2003) to simulate aerosol NH4 ,NO3 and SO4 (Tang et al., 2009), with a view to vali- the temporal trends, resulting in an increase in Vd for SO2 over dating European concentration fields of concentration and deposi- many parts of the continent, and driven by the long term, large scale tion for these species. The data show (Fig. 4.5) that the geometric decrease in the SO2/NH3 ratio (Fig. 4.4). The new scheme, for non- mean mixing ratios of SO2, HNO3 and HCl across the network are 0.4, stomatal resistance of both NH3 and SO2, incorporates an acidity to 0.35 and 0.15 ppb, respectively, so that, on average, SO2 makes up alkalinity (SO2/NH3) molar ratio as a scaling factor for resistances. For only about 40% of the sum of acids (SO2 þ HNO3 þ HCl). Further, the SO2, two non-stomatal resistances Rns,wet and Rns,dry are calculated as data indicate that at some sites (e.g. most Danish, French and Italian a function of the SO2/NH3 ratio, and a function of relative humidity is sites), the acidity is largely dominated by HNO3 and HCl, which are used for the transition from dry to wet when the surface cannot be considered in most models (e.g. Simpson et al., 2003) to be deposited w 1 considered fully wet or fully dry. at the maximum rates allowed by turbulence (Rc 0sm ). Under The EMEP non-stomatal deposition scheme has also been used such conditions, the proxy (SO2 þ HNO3 þ HCl)/NH3 would seem in field-scale inferential modelling of N and S dry deposition as part more appropriate to quantify the relative importance of surface of the NitroEurope project, using low-cost, long-term atmospheric acidity and alkalinity in model parameterisations, than the ratio of trace gas and aerosol DELTA samplers (Tang et al., 2009). Another SO2 alone to NH3. Clearly the surface affinity for SO2 uptake will implementation of the parameterisation was made by Zimmer- depend on the presence of fast-depositing, strong acids as the acidity mann et al. (2006) for the simulation of atmospheric deposition to is no longer SO2-dominated, and this needs to be accounted for in Norway spruce, using the SPRUCEDEP SVAT model, and comparison extended surface resistance parameterisations. with throughfall measurements and a canopy base cation budget model. The agreement between (inferential & canopy budget) 4.5. Future challenges modelling and observations was very good for S and oxidized N. Here, the contribution of dry to total (dry þ wet) deposition was The principal controls over SO2 deposition to terrestrial surfaces around 60% for S and for both reduced and oxidized N. have been identified from field, mainly micrometeorological The Dutch IDEM model (Bleeker et al., 2004; Erisman, 1994)also measurements. These studies have enabled the controlling steps in uses an NH3/SO2 molar ratio as a proxy for surface acidity. For NH3, the deposition pathway to be separated and their response to a range of default Rext values are used for 3 classes of the N/S ratio (very environmental variables quantified. In turn the data and responses low, low and high), depending on surface wetness, land-use and time have been used to develop process-based models and applied to 1 of day, while for SO2 the only effect implemented is to add 50 s m to quantify regional deposition budgets at country and continental the non-stomatal resistance when the N/S ratio is very low (<0.02). scales. There have been surprises, notably in the last decade. The The large reduction in European SO2 emissions and ambient largest surprise has been the recognition that long term (w1 year) concentrations over the last 25 years, and the relative stagnation in average deposition velocities change with time due to changes NH3 emissions and concentrations in Western Europe over the same in the chemical climatology at the regional scale. Thus a few