Modelling stomatal ozone flux and deposition to grassland communities across Europe

Regional scale modelling of both ozone deposition and the risk of ozone impacts is poorly developed for grassland communities. This paper presents new predictions of stomatal ozone flux to grasslands at five different locations in Europe, using a mechanistic model of canopy development for productive grasslands to generate time series of leaf area index and soil water potential as inputs to the stomatal component of the DO(3)SE ozone deposition model. The parameterisation of both models was based on Lolium perenne, a dominant species of productive pasture in Europe. The modelled seasonal time course of stomatal ozone flux to both the whole canopy and to upper leaves showed large differences between climatic zones, which depended on the timing of the start of the growing season, the effect of soil water potential, and the frequency of hay cuts. Values of modelled accumulated flux indices and the AOT40 index showed a five-fold difference between locations, but the locations with the highest flux differed depending on the index used; the period contributing to the accumulation of AOT40 did not always coincide with the modelled period of active ozone canopy uptake. Use of a fixed seasonal profile of leaf area index in the flux model produced very different estimates of annual accumulated total canopy and leaf ozone flux when compared with the flux model linked to a simulation of canopy growth. Regional scale model estimates of both the risks of ozone impacts and of total ozone deposition will be inaccurate unless the effects of climate and management in modifying grass canopy growth are incorporated.     

Modelling of stomatal conductance and ozone flux of Norway spruce: comparison with field data

It has been proposed that stomatal flux of ozone would provide a more reliable basis than ozone exposure indices for the assessment of the risk of ozone damage to vegetation across Europe. However, implementation of this approach requires the development of appropriate models which need to be rigorously tested against actual data collected under field conditions. This paper describes such an assessment of the stomatal component of the model described by Emberson et al. (2000. Modelling stomatal ozone flux across Europe. Environmental Pollution 110). Model predictions are compared with field measurements of both stomatal conductance (g(s)) and calculated ozone flux for shoots of mature Norway spruce (Picea abies) growing in the Tyrol Mountains in Austria. The model has been developed to calculate g(s) as a function of leaf phenology and four environmental variables: photosynthetic flux density (PFD), temperature, vapour pressure deficit (VPD) and soil moisture deficit (SMD). The model was run using climate data measured on site, although the SMD component was omitted since the necessary data were not available. The model parameterisation for Norway spruce had previously been collected from the scientific literature and therefore established independently from the measurement study. Overall, strong associations were found between model predictions and measured values of stomatal conductance to ozone (GO(3)) and calculated stomatal ozone flux (FO(3)). Average diurnal profiles of GO(3) and FO(3) showed good agreement between the field data and modelled values except during the morning period of 1990. The diurnal pattern of ozone flux was determined primarily by PFD and VPD, as there was little diurnal variation in ozone concentration. In general, the model predicted instances of high ozone flux satisfactorily, indicating its potential applicability in identifying areas of high ozone risk for this species.     

Comparison of different stomatal conductance algorithms for ozone flux modelling

A multiplicative and a semi-mechanistic, BWB-type [Ball, J.T., Woodrow, I.E., Berry, J.A., 1987. A model predicting stomatal conductance and its contribution to the control of photosynthesis under different environmental conditions. In: Biggens, J. (Ed.), Progress in Photosynthesis Research, vol. IV. Martinus Nijhoff, Dordrecht, pp. 221-224.] algorithm for calculating stomatal conductance (g(s)) at the leaf level have been parameterised for two crop and two tree species to test their use in regional scale ozone deposition modelling. The algorithms were tested against measured, site-specific data for durum wheat, grapevine, beech and birch of different European provenances. A direct comparison of both algorithms showed a similar performance in predicting hourly means and daily time-courses of g(s), whereas the multiplicative algorithm outperformed the BWB-type algorithm in modelling seasonal time-courses due to the inclusion of a phenology function. The re-parameterisation of the algorithms for local conditions in order to validate ozone deposition modelling on a European scale reveals the higher input requirements of the BWB-type algorithm as compared to the multiplicative algorithm because of the need of the former to model net photosynthesis (A(n)).     

Modelling annual pasture dynamics: Application to stomatal ozone deposition

Modelling ozone (O₃) deposition for impact risk assessment is still poorly developed for herbaceous vegetation, particularly for Mediterranean annual pastures. High inter-annual climatic variability in the Mediterranean area makes it difficult to develop models characterizing gas exchange behaviour and air pollutant absorption suitable for risk assessment. This paper presents a new model to estimate stomatal conductance (g_s) of Trifolium subterraneum, a characteristic species of dehesa pastures. The MEDPAS (MEDiterranean PAStures) model couples 3 modules estimating soil water content (SWC), vegetation growth and g_s. The g_s module is a reparameterized version of the stomatal component of the EMEP DO₃SE O₃ deposition model. The MEDPAS model was applied to two contrasting years representing typical dry and humid springs respectively and with different O₃ exposures. The MEDPAS model reproduced realistically the g_s seasonal and inter-annual variations observed in the field. SWC was identified as the major driver of differences across years. Despite the higher O₃ exposure in the dry year, meteorological conditions favoured 2.1 times higher g_s and 56 day longer growing season in the humid year compared to the dry year. This resulted in higher ozone fluxes absorbed by T. subterraneum in the humid year. High inter-family variability was found in gas exchange rates, therefore limiting the relevance of single species O₃ deposition flux modelling for dehesa pastures. Stomatal conductance dynamics at the canopy level need to be considered for more accurate O₃ flux modelling for present and future climate scenarios in the Mediterranean area.     

Modelling atmospheric mercury transport and deposition across Europe and the UK

There are inadequate measurements of surface ambient concentrations of mercury species and their deposition rates for the UK deposition budget to be characterized. In order to estimate the overall mercury flux budget for the UK, a simple long-term 1D Lagrangian trajectory model was constructed that treats emissions (1998), atmospheric transformation and deposition across Europe. The model was used to simulate surface concentrations of mercury and deposition across Europe at a resolution of 50 km×50 km and across the UK at 20 km×20 km. The model appeared to perform adequately when compared with the few available measurements, reproducing mean concentrations of elemental gaseous mercury at particular locations and the magnitude of regional gradients. The model showed that 68% of the UK's mercury emissions are exported and 32% deposited within the UK. Of deposition to the UK, 25% originates from the Northern Hemisphere/global background, 41% from UK sources and 33% from other European countries. The total mercury deposition to the UK is in good agreement with other modelling, 9.9 tonne yr⁻¹ cf. 9.0 tonne yr⁻¹, for 1998. However, the attribution differs greatly from the results of other coarser-scale modelling, which allocates 55% of the deposition to the UK from UK sources, 4% from other European countries and 60% from the global background atmosphere. The model was found to be sensitive to the speciation of emissions and the dry deposition velocity of elemental gaseous mercury. The uncertainties and deficiencies are discussed in terms of model parameterization and input data, and measurement data with which models can be validated. There is an urgent requirement for measurements of removal terms, concentrations, and deposition with which models can be parameterized and validated.     

Implications of climate change for the stomatal flux of ozone: a case study for winter wheat

Climate change factors such as elevated CO2 concentrations, warming and changes in precipitation affect the stomatal flux of ozone (O3) into leaves directly or indirectly by altering the stomatal conductance, atmospheric O3 concentrations, frequency and extent of pollution episodes and length of the growing season. Results of a case study for winter wheat indicate that in a future climate the exceedance of the flux-based critical level of O3 might be reduced across Europe, even when taking into account an increase in tropospheric background O3 concentration. In contrast, the exceedance of the concentration-based critical level of O3 will increase with the projected increase in tropospheric background O3 concentration. The influence of climate change should be considered when predicting the future effects of O3 on vegetation. There is a clear need for multi-factorial, open-air experiments to provide more realistic information for O3 flux-effect modelling in a future climate.     

A comparison of two stomatal conductance models for ozone flux modelling using data from two Brassica species

In this study we tested and compared a multiplicative stomatal model and a coupled semi-empirical stomatal-photosynthesis model in their ability to predict stomatal conductance to ozone (g_st) using leaf-level data from oilseed rape (Brassica napus L.) and broccoli (Brassica oleracea L. var. italica Plenck). For oilseed rape, the multiplicative model and the coupled model were able to explain 72% and 73% of the observed g_st variance, respectively. For broccoli, the models were able to explain 53% and 51% of the observed g_st variance, respectively. These results support the coupled semi-empirical stomatal-photosynthesis model as a valid alternative to the multiplicative stomatal model for O₃ flux modelling, in terms of predictive performance.     

Ozone risk assessment for agricultural crops in Europe: Further development of stomatal flux and flux–response relationships for European wheat and potato

Applications of a parameterised Jarvis-type multiplicative stomatal conductance model with data collated from open-top chamber experiments on field grown wheat and potato were used to derive relationships between relative yield and stomatal ozone uptake. The relationships were based on thirteen experiments from four European countries for wheat and seven experiments from four European countries for potato. The parameterisation of the conductance model was based both on an extensive literature review and primary data. Application of the stomatal conductance models to the open-top chamber experiments resulted in improved linear regressions between relative yield and ozone uptake compared to earlier stomatal conductance models, both for wheat (r²=0.83) and potato (r²=0.76). The improvement was largest for potato. The relationships with the highest correlation were obtained using a stomatal ozone flux threshold. For both wheat and potato the best performing exposure index was AF_st6 (accumulated stomatal flux of ozone above a flux rate threshold of 6 nmol ozone m⁻² projected sunlit leaf area, based on hourly values of ozone flux). The results demonstrate that flux-based models are now sufficiently well calibrated to be used with confidence to predict the effects of ozone on yield loss of major arable crops across Europe. Further studies, using innovations in stomatal conductance modelling and plant exposure experimentation, are needed if these models are to be further improved.     

Moving toward effective ozone flux assessment

We present a comment about "Ozone risk assessment for plants: central role of metabolism-dependent changes in reducing power" by Dizengremel, Le Thiec, Bagard, and Jolivet. As tools for summarizing plant O(3) sensitivity in simple indices, Dizengremel et al. suggest: reducing power, as antioxidant regeneration through the Halliwell/Asada cycle requires NADPH from the photosynthetic light reaction; Rubisco/PEPc ratio, as an index of the energy balance between anabolic and catabolic reactions; and water-use efficiency as a time-integrated approximation of the carbon gain to stomatal O(3) uptake ratio. The scientific background is solid, and simple enough (although expensive) to be translated into modelling and routine use. In the last decade, several approaches have been developed, mostly by using photosynthesis as a metric of defence. All these approaches should be experimentally tested in different and realistic conditions, before the results are transferred to the field and used in effective O(3) flux modelling and assessment.     

Modelling stomatal responses of spring wheat (Triticum aestivum L. cv. Turbo) to ozone and different levels of water supply

Spring wheat (Triticum aestivum L. cv. Turbo) was exposed to different levels of ozone and water supply in open-top chambers in 1991. Air was charcoal filtered (CF), non-filtered (NF) and CF plus proportional addition of ambient or twice ambient ozone (CF1, CF2). Seasonal means of O(3), taken over 24 h, were 2.3, 20.6, 17.3, and 34.5 nl litre(-1) for CF, NF, CF1 and CF2 treatments, respectively. A split-plot design was used to obtain two levels of water supply: one-half of the pots was irrigated sufficiently not to show any symptoms of drought stress; the others were exposed to low water supply and received 50% of these amounts. Using a steady-state porometer approximately 800 measurements of stomatal conductance (g(s)) were made on flag leaves from 68 to 106 days after sowing. The measurements yielded only small differences of maximum conductance between the two levels of water supply. Therefore, low water supply did not protect wheat plants against ozone injury via reduced stomatal uptake in this experiment. To describe the effects of environmental variables on the stomatal behaviour, boundary-line analysis and non-linear regression analysis were used. Besides microclimatic parameters, the ozone dose of flag leaves was introduced as an independent variable affecting stomatal aperture. A well-defined boundary line for ozone dose was found, suggesting that increasing ozone dose caused stomatal closure in wheat flag leaves. But at high ozone doses, co-acting senescence seems also responsible for the decrease in stomatal conductance. A multiplicative boundary-line model was used to predict stomatal conductance from combinations of environmental variables. In the test carried out with the measurements of stomatal conductance, the model accounted only for 40% of the variation of g(s). Generalized stomatal response patterns of the herbaceous growth form, the dependence of the variables' age and ozone dose and the lack of an important factor influencing stomatal response (water status of the plant) in the model, are suggested as explanations of the poor results of the test.     

Validation of the stomatal flux approach for the assessment of ozone visible injury in young forest trees. Results from the TOP (transboundary ozone pollution) experiment at Curno, Italy

This paper summarises some of the main results of a two-year experiment carried out in an Open-Top Chambers facility in Northern Italy. Seedlings of Populus nigra, Fagus sylvatica, Quercus robur and Fraxinus excelsior have been subjected to different ozone treatments (charcoal-filtered and non-filtered air) and soil moisture regimes (irrigated and non-irrigated plots). Stomatal conductance models were applied and parameterised under South Alpine environmental conditions and stomatal ozone fluxes have been calculated.The flux-based approach provided a better performance than AOT40 in predicting the onset of foliar visible injuries. Critical flux levels, related to visible leaf injury, are proposed for P. nigra and F. sylvatica (ranging between 30 and 33 mmol O₃ m⁻²). Soil water stress delayed visible injury appearance and development by limiting ozone uptake. Data from charcoal-filtered treatments suggest the existence of an hourly flux threshold, below which may occur a complete ozone detoxification.     

Simulation of stomatal conductance for Aleppo pine to estimate its ozone uptake

The data from a previous experiment carried out in open-top chambers to assess the effects of ozone (O3) exposure on growth and physiology of Aleppo pine (Pinus halepensis Mill.) were re-assessed to test the performance of the EMEP O3 stomatal conductance model used to estimate tree O3 uptake at a European scale. Aleppo pine seedlings were exposed during three consecutive years to three different O3 treatments: charcoal filtered air, non-filtered air and non-filtered air supplemented with 40 nl l(-1). The results of the model using the default parameterisation already published for Mediterranean conifers showed a poor performance when compared to measured data. Therefore, modifications of g(max), f(min), and new f(VPD), f(temp) and f(phen) functions were developed according to the observed data. This re-parameterisation resulted in a significant improvement of the performance of the model when compared to its original version.     

Long-term changes in the tropospheric column ozone from the ozone soundings over Europe

The tropospheric column of ozone is analyzed from the measurements of the vertical profile of ozone by balloon-born ozonesondes. The data base includes ∼16,000 ozone profiles collected above six European stations—three of them have begun the ozonesoundings since 1970. We present a trend analysis (with data up to 2005) focusing on detection of the long-term tropospheric ozone variability over the European network. The ozone time series have been examined separately for each station and season (DJF, MAM, JJA, SON) using a flexible trend model. A trend component of the model is taken as a smooth curve without a priori defined shape. A large increase in the European tropospheric ozone since the beginning of the 1970s (net change of ∼10% in summers and ∼30% in winters) and a kind of stabilization in the early 1990s have been corroborated by the study. This pattern comes from the most extensive data set of ozonesoundings over Hohenpeissenberg and Payern. The decadal differences in the trend pattern between these and other European stations are disclosed. The results of a stepwise regression model using various characteristics of the ozone and temperature profiles as explanatory variables for the tropospheric column ozone (TCO₃) variations show that the ozone changes may be reconstructed using the ozone mixing ratio at 500 hPa, the thermal tropopause (TT) height, and the difference between ozonepause and TT heights. It appears that the last two factors induce 20–30% of the net TCO₃ change over Hohenpeissenberg in the 1970–2004 period.     

Injury response of Phaseolus vulgaris to ozone flux density

This study describes a quantitative relationship between mean O₃ flux density and the length of exposure needed for the occurrence of visual injury to Phaseolus vulgaris L. Similar relationships were found for 14 day old and 6 week old plants using a whole leaf gas exchange cuvette system. Cultivars Seafarer (O₃ sensitive) and Gold Crop (O₃ resistant) exhibited similar responses at flux densities > 3 mg m⁻² h⁻¹ but only Seafarer was injured below this flux density. O₃ concentration and length of exposure period alone did not contain sufficient information to describe the onset of visual foliar injury. The use of O₃ concentrations in excess of normal ambient conditions compensated for low leaf conductances so that flux densities in the cuvette were similar to those found in the field.     

Critical levels for ozone effects on vegetation in Europe

The evidence of detrimental effects of ozone on vegetation in Europe, and the need to develop international control policies to reduce ozone exposures which are based on the effects of the pollutant, has led to attempts to define so-called critical levels of ozone above which adverse effects on trees, crops and natural vegetation may occur. This review is a critical assessment of the scientific basis of the concepts used to define critical levels for ozone and identifies the key limitations and uncertainties involved. The review focuses on the Level I critical level approach, which provides an environmental standard or threshold to minimise the effects of ozone on sensitive receptors, but does not seek to quantify the impacts of exceeding the critical level under field conditions. The concept of using the AOT (accumulated exposure over a threshold) to define long-term ozone exposure is demonstrated to be appropriate for several economically important species. The use of 40 ppb (giving the AOT40 index) as a threshold concentration gives a good linear fit to experimental data from open-top chambers for arable crops, but it is less certain that it provides the best fit to data for trees or semi-natural communities. Major uncertainties in defining critical level values relate to the choice of response parameter and species; the absence of data for many receptors, especially those of Mediterranean areas; and extrapolation to field conditions from relatively short-term open-top chamber experiments. The derivation of critical levels for long-lived organisms, such as forest trees, may require the use of modelling techniques based on physiological data from experimental studies. The exposure-response data which have been applied to derive critical levels should not be used to estimate the impacts of ozone over large areas, because of the uncertainties associated with extrapolation from the open-top chamber method, especially for forest trees, and because of spatial variation in atmospheric and environmental conditions, which may alter ozone uptake.     

Reactive uptake of ozone at simulated leaf surfaces: Implications for ‘non-stomatal’ ozone flux

The reaction of ozone (O₃) with α-pinene has been studied as a function of temperature and relative humidity and in the presence of wax surfaces that simulate a leaf surface. The objective was to determine whether the presence of a wax surface, in which α-pinene could dissolve and form a high surface concentration, would lead to enhanced reaction with O₃. The reaction of O₃ itself with the empty stainless steel reactor and with aluminium and wax surfaces demonstrated an apparent activation energy of around 30 kJ mol^?1 for all the surfaces, similar to that observed in long-term field measurements of O₃ fluxes to vegetation. However, the absolute reaction rate was 14 times greater for aluminium foil and saturated hydrocarbon wax surfaces than for stainless steel, and a further 5 times greater for beeswax than hydrocarbon wax. There was no systematic dependence on either relative or absolute humidity for these surface reactions over the range studied (20–100% RH). Reaction of O₃ with α-pinene occurred at rates close to those predicted for the homogeneous gas-phase reaction, and was similar for both the empty reactor and in the presence of wax surfaces. The hypothesis of enhanced reaction at leaf surfaces caused by enhanced surface concentrations of α-pinene was therefore rejected. Comparison of surface decomposition reactions on different surfaces as reported in the literature with the results obtained here demonstrates that the loss of ozone at the earth's surface by decomposition to molecular oxygen (i.e. without oxidative reaction with a substrate) can account for measured ‘non-stomatal’ deposition velocities of a few mm s^?1. In order to quantify such removal, the effective molecular surface area of the vegetation/soil canopy must be known. Such knowledge, combined with the observed temperature-dependence, provides necessary input to global-scale models of O₃ removal from the troposphere at the earth's surface.     

Leaf conductance response of phaseolus vulgaris to ozone flux density

The effect of ozone flux density on leaf conductance to ozone in Phaseolus vulgaris was examined. The change in conductance was measured within the first two hours of fumigation for mature, fruiting 6-week-old plants of an ozone sensitive cultivar (Seafarer); for young, 14-day-old plants of the same cultivar; and for an ozone resistant cultivar (Gold Crop). Young Seafarer plants showed no change in conductance to ozone over a wide range of ozone flux densities. Gold Crop showed a decrease in conductance of −3.1 % /(mgO₃ m⁻² h⁻¹) whereas mature Seafarer plants exhibited a stronger decrease of −7.7% /(mgO₃ m⁻² h⁻¹). Diffusion porometer measurements taken on fruiting Seafarer plants in the field illustrated that a decrease in leaf diffusive conductance to water is related to visual ozone injury.     

Simulations of stomatal conductance and ozone uptake to Norway spruce saplings in open-top chambers

A simulation model was developed to estimate the stomatal conductance and ozone flux to Norway spruce saplings in open-top chambers. The model was parameterized against needle conductance measurements that were made on 4-6-year-old spruce saplings, grown in open-top chambers, in July-September during three different seasons. The spruce saplings were either maintained well watered or subject to a 7-8 week drought period in July-September each year. The simulated conductance showed a good agreement with the measured conductance for the well-watered as well as the drought stress-treated saplings. The simulations were significantly improved when different vapour pressure deficit (VPD) functions were applied for well-watered and drought-stressed spruce saplings. The cumulated ozone uptake which was calculated from the conductance simulations showed less variation between years, compared to the cumulative ozone exposure index AOT40 (accumulated exposure over a threshold of 40 ppb or nl l(-1)) for the corresponding time periods. Measurements in May 1995 demonstrated the occurrence of long-term 'memory-effects' from the drought stress treatments on the conductance. Memory-effects need to be considered when simulation models for stomatal conductance are to be applied to long-lived forest trees under a multiple stress situation.     

Exposure to moderate concentrations of tropospheric ozone impairs tree stomatal response to carbon dioxide

With rising concentrations of both atmospheric carbon dioxide (CO₂) and tropospheric ozone (O₃), it is important to better understand the interacting effects of these two trace gases on plant physiology affecting land-atmosphere gas exchange. We investigated the effect of growth under elevated CO₂ and O₃, singly and in combination, on the primary short-term stomatal response to CO₂ concentration in paper birch at the Aspen FACE experiment. Leaves from trees grown in elevated CO₂ and/or O₃ exhibited weaker short-term responses of stomatal conductance to both an increase and a decrease in CO₂ concentration from current ambient level. The impairement of the stomatal CO₂ response by O₃ most likely developed progressively over the growing season as assessed by sap flux measurements. Our results suggest that expectations of plant water-savings and reduced stomatal air pollution uptake under rising atmospheric CO₂ may not hold for northern hardwood forests under concurrently rising tropospheric O₃.