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.     

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.     

Visible leaf injury in young trees of Fagus sylvatica L. and Quercus robur L. in relation to ozone uptake and ozone exposure. An Open-Top Chambers experiment in South Alpine environmental conditions

An Open-Top Chambers experiment on Fagus sylvatica and Quercus robur seedlings was conducted in order to compare the performance of an exposure-based (AOT40) and a flux-based approaches in predicting the appearance of ozone visible injuries on leaves. Three different ozone treatments (charcoal-filtered; non-filtered; and open plots) and two soil moisture treatments (watered and non-watered plots) were performed. A Jarvisian stomatal conductance model was drawn up and parameterised for both species and typical South Alpine environmental conditions, thus allowing the calculation of ozone stomatal fluxes for every treatment. A critical ozone flux level for the onset of leaf visible injury in beech was clearly identified between 32.6 and 33.6 mmolO3 m(-2). In contrast, it was not possible to identify an exposure critical level using the AOT40 index. Water stress delayed the onset of the leaf visible injuries, but the flux-based approach was able to take it into account accurately.     

Assessing the risk caused by ground level ozone to European forest trees: a case study in pine, beech and oak across different climate regions

Two different indices have been proposed for estimation of the risk caused to forest trees across Europe by ground-level ozone, (i) the concentration based AOT40 index (Accumulated Over a Threshold of 40 ppb) and (ii) the recently developed flux based AFstY index (Accumulated stomatal Flux above a flux threshold Y). This paper compares the AOT40 and AFstY indices for three forest trees species at different locations in Europe. The AFstY index is estimated using the DO(3)SE (Deposition of Ozone and Stomatal Exchange) model parameterized for Scots pine (Pinus sylvestris), beech (Fagus sylvatica) and holm oak (Quercus ilex). The results show a large difference in the perceived O(3) risk when using AOT40 and AFstY indices both between species and regions. The AOT40 index shows a strong north-south gradient across Europe, whereas there is little difference between regions in the modelled values of AFstY. There are significant differences in modelled AFstY between species, which are predominantly determined by differences in the timing and length of the growing season, the periods during which soil moisture deficit limits stomatal conductance, and adaptation to soil moisture stress. This emphasizes the importance of defining species-specific flux response variables to obtain a more accurate quantification of O(3) risk.     

Use of sap flow measurements to validate stomatal functions for mature beech (Fagus sylvatica) in view of ozone uptake calculations

For a quantitative estimate of the ozone effect on vegetation reliable models for ozone uptake through the stomata are needed. Because of the analogy of ozone uptake and transpiration it is possible to utilize measurements of water loss such as sap flow for quantification of ozone uptake. This technique was applied in three beech (Fagus sylvatica) stands in Switzerland. A canopy conductance was calculated from sap flow velocity and normalized to values between 0 and 1. It represents mainly stomatal conductance as the boundary layer resistance in forests is usually small. Based on this relative conductance, stomatal functions to describe the dependence on light, temperature, vapour pressure deficit and soil moisture were derived using multivariate nonlinear regression. These functions were validated by comparison with conductance values directly estimated from sap flow. The results corroborate the current flux parameterization for beech used in the DO₃SE model.     

Comparison of seasonal variations of ozone exposure and fluxes in a Mediterranean Holm oak forest between the exceptionally dry 2003 and the following year

Ozone and energy fluxes have been measured using the eddy covariance technique, from June to December 2004 in Castelporziano near Rome (Italy), and compared to similar measurements made in the previous year. The studied ecosystem consisted in a typical Mediterranean Holm oak forest. Stomatal fluxes have been calculated using the resistance analogy and by inverting the Penmann-Monteith equation. Results showed that the average stomatal contribution accounts for 42.6% of the total fluxes. Non-stomatal deposition proved to be enhanced by increasing leaf wetness and air humidity during the autumnal months. From a comparison of the two years, it can be inferred that water supply is the most important limiting factor for ozone uptake and that prolonged droughts alter significantly the stomatal conductance, even 2 months after the soil water content is replenished. Ozone exposure, expressed as AOT40, behaves similarly to the cumulated stomatal flux in dry conditions whereas a different behaviour for the two indices appears in wet autumnal conditions. A difference also occurs between the two years.     

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.     

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 critical levels of ozone for the forested area of austria modifications of the aot40 concept

GOAL, SCOPE AND BACKGROUND: Ozone is the most important air pollutant in Europe for forest ecosystems and the increase in the last decades is significant. The ozone impact on forests can be calculated and mapped based on the provisional European Critical Level (AOT40 = accumulated exposure over a threshold of 40 ppb, 10,000 ppb x h for 6 months of one growing season calculated for 24 h day(-1)). For Norway spruce, the Austrian main tree species, the ozone risk was assessed in a basis approach and because the calculations do not reflect the health status of forests in Austria, the AOT40 concept was developed. METHODS: Three approaches were outlined and maps were generated for Norway spruce forests covering the entire area of Austria. The 1st approach modifies the AOT40 due to the assumption that forests have adapted to the pre-industrial levels of ozone, which increase with altitude (AOTalt). The 2nd approach modifies the AOT40 according to the ozone concentration in the sub-stomata cavity. This approach is based on such factors as light intensity and water vapour saturation deficit, which affect stomatal uptake (AOTsto). The 3rd approach combines both approaches and includes the hemeroby. The pre-industrial ozone level approach was applied for autochthonous ('natural') forest areas, the ozone-uptake approach for non-autochthonous ('altered') forest areas. RESULTS AND DISCUSSION: The provisional Critical Level (AOT40) was established to allow a uniform assessment of the ozone risk for forested areas in Europe. In Austria, where ozone risk is assessed with utmost accuracy due to the dense grid of monitoring plots of the Forest Inventory and because the continuously collected data from more than 100 air quality measuring stations, an exceedance up to the five fold of the Critical Level was found. The result could lead to a yield loss of up to 30-40% and to a severe deterioration in the forest health status. However, the data of the Austrian Forest Inventory and the Austrian Forest Damage Monitoring System do not reflect such an ozone impact. Therefore, various approaches were outlined including the tolerance and avoidance mechanisms of Norway spruce against ozone impact. Taking into consideration the adaptation of forests to the pre-industrial background level of ozone, the AOT40 exceedances are markedly reduced (1st approach). Taking into account the stomatal uptake of ozone, unrealistic high amounts of exceedances up to 10,000 ppb x h were found. The modelled risk does not correspond with the health status and the wood increment of the Austrian forests (2nd approach). Consolidating the forgoing two approaches, a final map including the hemeroby was generated. It became clear that the less natural ('altered') forested regions are highly polluted. This means, that more than half of the spruce forests are endangered by ozone impact and AOT40 values of up to 30,000 ppb x h occur (3rd approach). CONCLUSIONS: The approaches revealed that a plausible result concerning the ozone impact on spruce forests in Austria could only be reached by combining pre-industrial ozone levels, ozone flux into the spruce needles and the hemeroby of forests.     

Comparison of modelled and measured ozone concentrations and meteorology for a site in south-west Sweden: implications for ozone uptake calculations

Measurements of ground-level ozone concentrations and meteorology (temperature, vapour pressure deficit (VPD), solar radiation) at the monitoring site Ostad (south-west Sweden) were compared to data from the corresponding grid in the EMEP photo-oxidant model for 1997, 1999 and 2000. The influence of synoptic weather on the agreement between model and measurements was studied. Implications of differences between modelled and observed inputs for ozone flux calculations for wheat and potato were investigated. The EMEP model output of ozone, temperature and VPD correlated well with measurements during daytime. Deviations were larger during the night, especially in calm conditions, attributed to local climatological conditions at the monitoring site deviating from average conditions of the grid. These differences did not lead to significant differences in calculated ozone uptake, which was reproduced remarkably well. The uptake calculations were sensitive to errors in the ozone and temperature input data, especially when including a flux threshold.     

Development of a multi-factor model for predicting the effects of ambient ozone on the biomass of white clover

Results are presented from the UN/ECE ICP Vegetation (International Cooperative Programme on effects of air pollution on natural vegetation and crops) experiments in which ozone(O(3))-resistant (NC-R) and -sensitive (NC-S) clones of white clover (Trifolium repens cv. Regal) were exposed to ambient O(3) episodes at 14 sites in eight European countries in 1996, 1997 and 1998. The plants were grown according to a standard protocol, and the forage was harvested every 28 days for 4-5 months per year by excision 7 cm above the soil surface. Biomass ratio (NC-S/NC-R) was related to the climatic and pollutant conditions at each site using multiple linear regression (MLR) and artificial neural networks (ANNs). Twenty-one input parameters [e.g. AOT40, 7-h mean O(3) concentration, daylight vapour pressure deficit (VPD), daily maximum temperature] were considered individually and in combination with the aim of developing a model with high r(2) and simple structure that could be used to predict biomass change in white clover. MLR models were generally more complex, and performed less well for unseen data than non-linear ANN models. The ANN model with the best performance had five inputs with an r(2) value of 0.84 for the training data, and 0.71 for previously unseen data. Two inputs to the model described the O(3) conditions (AOT40 and 24-h mean for O(3)), two described temperature (daylight mean and 24-h mean temperature), and the fifth input appeared to be differentiating between semi-urban and rural sites (NO concentration at 17:00). Neither VPD nor harvest interval was an important component of the model. The model predicted that a 5% reduction in biomass ratio was associated with AOT40s in the range 0.9-1.7 ppm x h (microl l(-1) h) accumulated over 28 days, with plants being most sensitive in conditions of low NO(x), medium-range temperature, and high 24-h mean O(3) concentration.     

Establishing ozone flux–response relationships for winter wheat: Analysis of uncertainties based on data for UK and Polish genotypes

The work outlined in this paper had three objectives. The first was to explore the effects of ozone pollution on grain yield and quality of commercially-grown winter wheat cultivars. The second was to derive a stomatal ozone flux model for winter wheat and compare with those already developed for spring wheat. The third was to evaluate exposure- versus flux–response approaches from a risk assessment perspective, and explore the implications of genetic variation in modelled ozone flux.Fifteen winter wheat cultivars were grown in open-top chambers where they were exposed to four levels of ozone. During fumigation, stomatal conductance measurements were made over the lifespan of the flag leaf across a range of environmental conditions. Although significant intra-specific variation in ‘ozone sensitivity’ (in terms of impacts on yield) was identified, yield was inversely related (R² = 0.63, P < 0.001) to the accumulated hourly averaged ozone exposure above 40 ppb during daylight hours (AOT40) across the dataset. The adverse effect of ozone on yield was principally due to a decline in seed weight. Algorithms defining the influence of environmental variables on stomatal uptake were subtly different from those currently in use, based on data for spring wheat, to map ozone impacts on pan-European cereal yield. Considerable intra-specific variation in phenological effects was identified. This meant that an ‘average behaviour’ had to be derived which reduced the predictive capability of the derived stomatal flux model (R² = 0.49, P < 0.001, 15 cultivars included). Indeed, given the intra-specific variability encountered, the flux model that was derived from the full dataset was no better in predicting O₃ impacts on wheat yield than was the AOT40 index. The study highlights the need to use ozone risk assessment tools appropriate to specific vegetation types when modelling and mapping ozone impacts at the regional level.     

High-resolution modelling of AOT40 and stomatal ozone uptake in wheat and grassland: a comparison between 2000 and the hot summer of 2003 in Switzerland

The aim was to compare the ozone risk for agricultural crops in Switzerland during the hot and dry year 2003 with the more 'normal' situation in 2000. An improved version of the Ozone DEposition Model ODEM was used at a 2 x 2 km resolution. The distribution of the index AOT40 was compared with the accumulated stomatal ozone flux, AF(st). Averaged AOT40 at 2 m and at canopy height was much higher in 2003 than in 2000, but inter-annual differences in AF(st) for wheat and grasslands were small due to the limiting effect of low soil water contents in 2003. AOT40 suggested larger potential yield losses in wheat in 2003, while using AF(st) with a threshold of 6 nmol m(-2) s(-1) (AF(st)6) yielded similar estimates for both years. The data show that modelling of AF(st) can be used to differentiate ozone risks between regions and years at a national scale.     

Relations between crown condition and ozone and its dependence on environmental factors

It is now widely accepted that stomatal ozone uptake provides a more reliable indicator for the assessment of risk of ozone damage to forest trees than ozone exposure. Thus, environmental factors influencing stomatal opening strongly affect the sensitivity of forest trees towards high ozone concentrations. In this study the hydrological model WAWAHAMO is used to simulate stomatal opening dependent on air temperature, solar radiation, vapor pressure deficit and soil water potential. Then, a statistical analysis was carried out to investigate relations between crown defoliation and ozone concentrations dependent on stomatal opening. The results support the hypothesis that tropospheric ozone contributes to the observed crown defoliation. However, environmental factors strongly affect tree response to high ozone concentrations.     

An evaluation of ozone exposure metrics for a seasonally drought-stressed ponderosa pine ecosystem

Ozone stress has become an increasingly significant factor in cases of forest decline reported throughout the world. Current metrics to estimate ozone exposure for forest trees are derived from atmospheric concentrations and assume that the forest is physiologically active at all times of the growing season. This may be inaccurate in regions with a Mediterranean climate, such as California and the Pacific Northwest, where peak physiological activity occurs early in the season to take advantage of high soil moisture and does not correspond to peak ozone concentrations. It may also misrepresent ecosystems experiencing non-average climate conditions such as drought years. We compared direct measurements of ozone flux into a ponderosa pine canopy with a suite of the most common ozone exposure metrics to determine which best correlated with actual ozone uptake by the forest. Of the metrics we assessed, SUM0 (the sum of all daytime ozone concentrations > 0) best corresponded to ozone uptake by ponderosa pine, however the correlation was only strong at times when the stomata were unconstrained by site moisture conditions. In the early growing season (May and June). SUM0 was an adequate metric for forest ozone exposure. Later in the season, when stomatal conductance was limited by drought. SUM0 overestimated ozone uptake. A better metric for seasonally drought-stressed forests would be one that incorporates forest physiological activity, either through mechanistic modeling, by weighting ozone concentrations by stomatal conductance, or by weighting concentrations by site moisture conditions.     

Quantification of ozone uptake at the stand level in a Pinus canariensis forest in Tenerife, Canary Islands: an approach based on sap flow measurements

Ozone uptake was studied in a pine forest in Tenerife, Canary Islands, an ecotone with strong seasonal changes in climate. Ambient ozone concentration showed a pronounced seasonal course with high concentrations during the dry and warm period and low concentrations during the wet and cold season. Ozone uptake by contrast showed no clear seasonal trend. This is because canopy conductance significantly decreased with soil water availability and vapour pressure deficit. Mean daily ozone uptake averaged 1.9 nmol m(-2) s(-1) during the wet and cold season, and 1.5 nmol m(-2) s(-1) during the warm and dry period. The corresponding daily mean ambient ozone concentrations were 42 and 51 nl l(-1), respectively. Thus we conclude that in Mediterranean type forest ecosystems the flux based approach is more capable for risk assessment than an external, concentration based approach.     

Ozone fluxes above and within a pine forest canopy in dry and wet conditions

The physiological and physical processes controlling ozone dry deposition to vegetated surfaces are still not fully understood. In particular, the role of the understorey and the possible action of dew on ozone deposition have not received much attention so far. This paper presents the results of an experiment aimed at quantifying ozone dry deposition to a maritime pine forest in the “Les Landes” area in France. Ozone deposition fluxes were measured using the eddy-covariance technique above and within the canopy. We investigate the factors acting on ozone deposition in both dew-wetted and dry conditions. The values obtained for the ozone deposition velocity are well in the range of previously published measurements over coniferous forests. For the present forest, ozone uptake by the understorey is a significant portion of ozone deposition to the whole pine stand. The understorey contributes more to the overall ozone flux than to the other measured scalar fluxes (sensible heat and water vapour). During dry nights the surface conductance for ozone and the friction velocity are strongly correlated, showing that ozone deposition is largely controlled by dynamical processes. During the day, in dry conditions, the canopy stomatal conductance is the major parameter controlling ozone deposition. However, in winter, when the stomatal conductance is low, the influence of dynamical processes persists during day-time. It is also found that surface wetness associated with dew significantly enhances ozone deposition, during the night as well as in the morning.     

Modeling coupled interactions of carbon,water, and ozone exchange between terrestrial ecosystems and the atmosphere. I: model description

A new biophysical model (FORFLUX) is presented to study the simultaneous exchange of ozone, carbon dioxide, and water vapor between terrestrial ecosystems and the atmosphere. The model mechanistically couples all major processes controlling ecosystem flows trace gases and water implementing recent concepts in plant eco-physiology, micrometeorology, and soil hydrology. FORFLUX consists of four interconnected modules-a leaf photosynthesis model, a canopy flux model, a soil heat-, water- and CO2- transport model, and a snow pack model. Photosynthesis, water-vapor flux and ozone uptake at the leaf level are computed by the LEAFC3 sub-model. The canopy module scales leaf responses to a stand level by numerical integration of the LEAFC3model over canopy leaf area index (LAI). The integration takes into account (1) radiative transfer inside the canopy, (2) variation of foliage photosynthetic capacity with canopy depth, (3) wind speed attenuation throughout the canopy, and (4) rainfall interception by foliage elements. The soil module uses principles of the diffusion theory to predict temperature and moisture dynamics within the soil column, evaporation, and CO2 efflux from soil. The effect of soil heterogeneity on field-scale fluxes is simulated employing the Bresler-Dagan stochastic concept. The accumulation and melt of snow on the ground is predicted using an explicit energy balance approach. Ozone deposition is modeled as a sum of three fluxes- ozone uptake via plant stomata, deposition to non-transpiring plant surfaces, and ozone flux into the ground. All biophysical interactions are computed hourly while model projections are made at either hourly or daily time step. FORFLUX represents a comprehensive approach to studying ozone deposition and its link to carbon and water cycles in terrestrial ecosystems.     

Environmental control of stomatal conductance in forest trees of the Great Smoky Mountains National Park

To determine if stomatal conductance (g(s)) of forest trees could be predicted from measures of leaf microclimate, diurnal variability in in situ g(s) was measured in black cherry (Prunus serotina), red maple (Acer rubrum), and northern red oak (Quercus rubra). Relative to overstory trees, understory saplings exhibited little diurnal variability in g(s) and ozone uptake. Depending on species and site, up to 30% of diurnal and seasonal variation in g(s )of overstory trees was explained by photosynthetically active radiation and vapor pressure deficit. Daily maximum g(s) was significantly related to soil moisture in overstory northern red oak and black cherry (R(2) ranged from 33 to 65%). Although g(s) was not fully predicted using instantaneous measures of leaf microclimate, ozone uptake of large forest trees was reduced by low soil moisture.     

Increasing risk for negative ozone impacts on vegetation in northern Sweden

Trends were found for increasing surface ozone concentrations during April-September in northern Sweden over the period 1990-2006 as well as for an earlier onset of vegetation growing season. The highest ozone concentrations in northern Sweden occurred in April and the ozone concentrations in April showed a strong increasing trend. A model simulation of ozone flux for Norway spruce indicated that the provisional ozone flux based critical level for forests in Europe is exceeded in northern Sweden. Future climate change would have counteracting effects on the stomatal conductance and needle ozone uptake, mediated on the one hand by direct effect of increasing air temperatures and on the other through increasing water vapour pressure difference between the needles and air. Thus, there is a substantial and increasing risk for negative impacts of ozone on vegetation in northern Sweden, related mainly to increasing ozone concentrations and an earlier onset of the growing season.