In this paper, some shortcomings involved in the modelling of ozone fluxes in the context of local-scale risk assessment are discussed, especially as related to the data collected within the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests). An enhanced monitoring strategy, that would provide a sounder basis for the development, validation and application of risk assessment modelling tools, is also suggested.
OzoneDoseStomatal fluxForestsMonitoringPollution and Stress FactorsIntroduction
In Europe, the development of air pollution abatement strategies is founded on the effects-based approach, which includes the ozone-induced plant injury as one of the key effects to be minimized. Within this framework, the risk of ozone damage to vegetation is related to numerical exposure and dose indices (UNECE 2004b, Simpson et al. 2007, Tuovinen et al. 2007). Both types of index are presently defined within the risk assessment methodology adopted within the Convention of Long-range Transboundary Air Pollution of the United Nations Economic Commission for Europe (UNECE 2004b). However, there is increasing evidence for the superior biological basis of the dose approach (Matyssek et al. 2007).
The exposure and dose indices differ in that exposure can be evaluated from ozone concentration data alone, while for the calculation of ozone dose the stomatal uptake by vegetation must be modelled. In this paper, we will discuss some shortcomings involved in the modelling of ozone fluxes in the context of local-scale risk assessment, especially as related to the data collected within the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests - UNECE 2007). We will also suggest an enhanced monitoring strategy that would provide a sounder basis for the development, validation and application of risk assessment modelling tools.
Flux modelling principles
The dose-based ozone risk indicator (AFstY, Accumulated stomatal Flux above a threshold Y) introduced within the UNECE risk assessment methodology (UNECE 2004b) can be written as (eqn. 1):
AF_{st}Y = \sum_{i=1}^{N} max(F_{st,i}-Y,0) \cdot \Delta t
where Fst is the stomatal ozone flux per projected leaf area (PLA) to sunlit leaves at the canopy top and Y is the threshold stomatal flux per PLA (in nmol m-2 s-1). AFstY is calculated from hourly values of Fst (denoted by i), so Δt = 1 h; N denotes the number of hours to be included in the calculation period, which corresponds to the growing season. The stomatal flux in eqn. 1 is defined as (eqn. 2):
F_{st,1} = c_i(h_{veg}) \hat{g}_{st,i}(h_{veg})
where ci is the hourly ozone concentration, hveg is vegetation height and gst,i is the effective stomatal conductance (here referred to as “effective” because it also depends on the conductances of the leaf boundary layer and the external plant surfaces, see Tuovinen et al. 2007). The stomatal conductance is modelled using the DO3SE (Deposition of Ozone and Stomatal Exchange) model, which is based on a multiplicative plant species-specific parametrization representing the stomatal responses to environmental and phenological factors (Emberson et al. 2000, UNECE 2004b). Thus the stomatal flux depends on two components, a concentration and a stomatal conductance, both of which are equally significant for the flux.
Measuring concentration
As indicated by eqn. 1 and eqn. 2, the stomatal flux for AFstY, by definition, is to be calculated from hourly-averaged ozone concentration data. These data could be obtained from a standard ozone analyser that is based on UV absorption photometry; this is the reference method defined in the EU Directive on ambient ozone (EU 2002). Within ICP Forests, passive sampling is defined as an option for concentration measurements (UNECE 2000) and, being relatively inexpensive and easy to deploy in the field, is widely used at the ICP Forests Level II monitoring plots across Europe (UNECE 2007). A comparison of passive samplers against the reference method at some ICP Forests sites demonstrated the feasibility of this method but also showed the associated methodological uncertainty (Sanz et al. 2007).
A fundamental property of passive sampling is the time-averaging of the measurement. As a typical sampling time for ozone is two weeks (UNECE 2000), the concentration data obtained have a low temporal resolution as compared to the definition of many air quality indicators, and especially so for AFstY. Thus the hourly data must be derived from the measured mean (14-d or so) concentration by using a statistical technique with additional meteorological (Krupa et al. 2003) or topographical (Loibl et al. 1994) data. This requires a significant amount of prior (hourly) calibration data, and unavoidably further uncertainty is introduced in the modelled hourly values, as exemplified by the results of Gerosa et al. (2007). Moreover, it is questionable to what extent the correlation between high ozone concentrations and environmental factors limiting stomatal uptake, which was a major motivation for the flux-based approach in the first place, can be simulated by this approach.
According to the monitoring recommendations of ICP Forests, the passive samplers are to be located in an open field near, but outside, the forest at a 2-4-m height (UNECE 2000). This contrasts the definition of stomatal flux (eqn. 2), which should be calculated using the tree-top concentration (at hveg ~ 20 m). Consequently, the measured concentration must be transformed to the correct reference height (Fig. 1). This can be accomplished by using a flux-gradient model that relates the vertical concentration profile to wind speed, surface roughness and the bulk surface conductance of ozone deposition (UNECE 2004b). In principle, this kind of model is founded on a micrometeorological theory which assumes horizontal homogeneity. However, this assumption is not fulfilled very well within the present experimental configuration. Furthermore, neutral atmospheric stability, constant surface conductance and independence of the 20-m concentration of the underlying vegetation were employed as necessary assumptions by Schaub et al. (2007), who modelled ozone doses using the ICP Forests monitoring data.
Modelling stomatal conductance
The modelling of stomatal conductance with DO3SE depends on data on wind speed, photosynthetically active radiation, air temperature and air humidity at hveg, and data on soil moisture (Emberson et al. 2000, UNECE 2004b). According to UNECE (2004a), the recommended height for wind measurements is 10 m, while for most of the other meteorological variables it is 2 m. This means that the measured data should be transformed to hveg, using a micrometeorological flux-gradient relationship similar to that employed for the vertical concentration profile. In addition, since soil moisture status is not measured at the ICP Forests monitoring plots, it needs to be modelled based on water budget principles when applying the DO3SE model (Schaub et al. 2007).
A general problem related to flux modelling is the limited validation of models, and here we can note several issues. Firstly, in many cases there are little independent data available for statistically sound cross-validation. This results in conceptual problems concerning the distinction between model calibration and validation, as is evident from the evaluations of the canopy-scale version of DO3SE (Tuovinen et al. 2004).
Secondly, field experiments are seldom designed from the point of view of the characteristics of a certain model, so the data provided by these efforts may be far from optimal for model validation. For example, the observational data may not cover a full range of environmental conditions, as encountered when running the model for a complete growing season, or all the necessary input data, such as soil moisture, are not measured at all (Tuovinen et al. 2004). In particular, the partitioning between stomatal and non-stomatal fluxes would be essential for evaluating flux models. At the canopy scale, this would require measurements of water vapour exchange and/or xylem sap flow. Finally, air quality monitoring programmes, such as that run at the ICP Forests Level II monitoring plots, provide little support for the validation of flux models, as the ozone (or water vapour) flux is not measured.
Discussion
While the passive sampling technique for measuring ozone concentrations offers many advantages over continuous monitoring, especially in a remote and complex forest environment, we have here identified several uncertainty sources specific to the application of this technique to the modelling of ozone fluxes. Some of these issues arise from the monitoring recommendations provided by ICP Forests. All the concentration-related uncertainties discussed above would be avoided, if the ozone concentration was measured at the canopy top using a UV absorption analyser providing hourly-resolved data. In addition, the location of meteorological sensors is not ideal for flux modelling.
In general terms, it can be argued that we are dealing with a trade-off situation as regards the ozone monitoring strategy. On the one hand we have the inexpensive alternative relying on passive sampling; on the other, we could replace this by the costly continuous monitors. The inexpensive alternative makes it possible to run an extensive measurement network. In this case, however, the application of flux-based risk indicators, which should gradually replace the concentration-based indices, entails a large number of additional calculation steps with associated simplifying assumptions. The overall uncertainty involved in these simplifications remains to be quantified. The alternative based on continuous monitors is too costly to be implemented across the existing ICP Forests network. In addition, measurements of ozone deposition fluxes would be needed for model development and validation.
In order to avoid the straightforward trade-off outlined above, an alternative approach could be adopted. We suggest establishing, in addition to the existing Level II sites, a small number of well-equipped measurement sites, “supersites” or “Level III sites”, that would provide data specifically for the flux-based risk assessment purposes. For the calculation of Fst and AFstY, the canopy-top ozone concentration (possibly with vertical within-canopy profiles) and all the input data needed for DO3SE would be measured at these sites on an hourly basis. For further development and validation of DO3SE and other flux models, ozone and water vapour fluxes would be measured above the canopy using the micrometeorological eddy covariance technique, possibly enhanced by sap flow and shoot-scale gas exchange measurements. To overcome the financial and logistic restrictions, all this could be accomplished in practice by collaborating with the existing flux measurement stations (“flux towers”) run across Europe, mainly within large-scale projects such as CarboEurope (2004) and NitroEurope (2006) (Fig. 2). These stations constitute a potential framework for the supersites, providing the necessary infrastructure and expertise, extensive measurement programmes and data bases, as well as a direct connection to flourishing research on atmosphere-biosphere exchanges (e.g., Piao et al. 2008).
Acknowledgements
Discussions with Lisa Emberson (SEI-York), Dave Simpson (EMEP MSC-W) and Mika Aurela (FMI) are highly appreciated.
ReferencesCarboEuropeCarboEurope-IP, Assessment of the European Terrestrial Carbon Balance, Integrated Project, Strategic Plan 2004-2008. Max-Planck-Institute for Biogeochemistry, Jena, Germany, pp. 28.2004Emberson LD, Ashmore MR, Cambridge HM, Simpson D, Tuovinen J-PModelling stomatal ozone flux across Europe. Environmental Pollution 109: 403-413.2000EUDirective 2002/3/EC of the European Parliament and the Council of 12 February 2002 relating to ozone in ambient air. Official Journal of European Communities L 67: 14-30.2002Gerosa G, Ferretti M, Bussotti F, Rocchini DEstimates of ozone AOT40 from passive sampling in forest sites in South-Western Europe. Environmental Pollution 145: 629-635.2007Krupa S, Nosal M, Ferdinand JA, Stevenson RE, Skelly JMA multi-variate statistical model integrating passive sampler and meteorology data to predict the frequency distributions of hourly ambient ozone (O3) concentrations. Environmental Pollution 124: 173-178.2003Loibl W, Winiwarter W, Kopsca A, Zueger J, Baumann REstimating the spatial distribution of ozone concentrations in complex terrain. Atmospheric Environment 28: 2557-2566.1994Matyssek R, Bytnerowicz A, Karlsson P-E, Paoletti E, Sanz M, Schaub M, Wieser GPromoting the O3 flux concept for European forests. Environmental Pollution 146: 587-607.2007NitroEuropeNitroEurope IP brochure. Centre for Ecology and Hydrology, Edinburgh, UK.2006Piao S, Ciais P, Friedlingstein P, Peylin P, Reichstein M, Luyssaert S, Margolis H, Fang J, Barr A, Chen A, Grelle A, Hollinger DY, Laurila T, Lindroth A, Richardson AD, Vesala TNet carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 451: 49-52.2008Sanz MJ, Calatayud V, Sánchez-Peña GMeasures of ozone concentrations using passive sampling in forests of South Western Europe. Environmental Pollution 145: 620-628.2007Schaub M, Emberson L, Büker P, Kräuchi NPreliminary results of modeled ozone uptake for Fagus sylvatica L. trees at selected EU/UN-ECE intensive monitoring plots. Environmental Pollution 145: 636-643.2007Simpson D, Ashmore M, Emberson L, Tuovinen J-PA comparison of two different approaches for mapping potential ozone damage to vegetation. A model study. Environmental Pollution 146: 715-725.2007Tuovinen J-P, Ashmore M, Emberson L, Simpson DTesting and improving the EMEP ozone deposition module. Atmospheric Environment 38: 2373-2385.2004Tuovinen J-P, Simpson D, Emberson L, Ashmore M, Gerosa GRobustness of modelled ozone exposures and doses. Environmental Pollution 146: 578-586.2007UNECEManual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests. Part X: Monitoring of air quality. International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests, United Nations Economic Commission for Europe (UNECE) Convention on Long-range Transboundary Air Pollution, Geneva, pp. 20.2000UNECEManual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests. Part VII: Meteorological Measurements. International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests, United Nations Economic Commission for Europe (UNECE) Convention on Long-range Transboundary Air Pollution, Geneva, Switzerland, pp. 40.2004aUNECEMapping critical levels for vegetation. Manual on methodologies and criteria for modelling and mapping critical loads and levels and air pollution effects, risks and trends. United Nations Economic Commission for Europe (UNECE) Convention on Long-range Transboundary Air Pollution, Geneva, Switzerland, pp. 254.2004bUNECEThe condition of forests in Europe, 2007 executive report. United Nations Economic Commission for Europe (UNECE) Convention on Long-range Transboundary Air Pollution, Geneva, Switzerland, pp. 34.2007
Schematic presentation of the height transformation for estimating the canopy-top concentration.
The flux measurement network of the CarboEurope project (CarboEurope 2004).