Modelling stomatal ozone flux across Europe

A model has been developed to estimate stomatal ozone flux across Europe for a number of important species. An initial application of this model is illustrated for two species, wheat and beech. The model calculates ozone flux using European Monitoring and Evaluation Programme (EMEP) model ozone concentrations in combination with estimates of the atmospheric, boundary layer and stomatal resistances to ozone transfer. The model simulates the effect of phenology, irradiance, temperature, vapour pressure deficit and soil moisture deficit on stomatal conductance. These species-specific microclimatic parameters are derived from meteorological data provided by the Norwegian Meteorological Institute (DNMI), together with detailed land-use and soil type maps assembled at the Stockholm Environment Institute (SEI). Modelled fluxes are presented as mean monthly flux maps and compared with maps describing equivalent values of AOT40 (accumulated exposure over threshold of 40 ppb or nl l(-1)), highlighting the spatial differences between these two indices. In many cases high ozone fluxes were modelled in association with only moderate AOT40 values. The factors most important in limiting ozone uptake under the model assumptions were vapour pressure deficit (VPD), soil moisture deficit (for Mediterranean regions in particular) and phenology. The limiting effect of VPD on ozone uptake was especially apparent, since high VPDs resulting in stomatal closure tended to co-occur with high ozone concentrations. Although further work is needed to link the ozone uptake and deposition model components, and to validate the model with field measurements, the present results give a clear indication of the possible implications of adopting a flux-based approach for future policy evaluation.     

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.     

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.     

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.     

Variability of ozone concentration in a montane environment,northern Italy

To evaluate the spatial variability of ozone concentrations, two studies were undertaken in the montane environment of Trentino region, northern Italy, in 2007. In the first study, a 225 km² area was considered. Here, a randomized design was used to evaluate the variability of ozone concentration at 1 and 225 km² scale. Measurements were carried out by passive samplers between May and June 2007. In a second study, the whole 6207 km² area of Trentino was considered. The area is covered by five grid cells of the European Monitoring and Evaluation Programme (EMEP). A systematic 15 × 15 km grid was used to allocate 15 passive samplers over the entire province, resulting into 1–4 samplers for each of the 5 EMEP grid cells (2500 km² each) overlapping the study area. Measurements were carried out between June and September 2007. Accuracy of passive samplers was checked by direct comparison with conventional ozone analysers. Significant differences (P = 0.034) were found in ozone concentration among 1 × 1 km quadrates within the 225 km² study area, while variability within the 1 × 1 km grid cells (coefficient of variation, CV′ = 0.12) slightly exceed the measurement error (CV′ = 0.08). At larger scales (225, 2500 and 6207 km²), ozone concentration shows much higher variability (CV′ from 0.18 to 0.28, with peak values at 0.40). Reported differences lead to very different AOT40 estimates even within the same EMEP grid cell. These findings suggest that 1 × 1 km resolution seems appropriate for ozone concentration modelling. On the other hand, significant sub-grid variation may exist at the resolution adopted by the EMEP model. Coupled with the likely variability of other important meteorological, soil and vegetation variables, our findings suggest that ozone risk assessment for vegetation based on large-scale modelled AOT40 and flux needs to be considered with great caution. The evidence reported in this paper asks for more detailed national-scale modelling, and the development of methods to incorporate local scale variations into large-scale models.     

Structural analysis and flux associations of CO2, H2O,heat and ozone over cotton and grape canopies

Micrometeorological tower data, collected over grape and cotton canopies as part of the California ozone deposition experiment (CODE) during the summer of 1991, are used to examine the temporal association between fluxes, and the physical characteristics of the coherent structures which dominate transport for both stable nighttime and unstable daytime conditions. Flux was calculated using the eddy covariance technique and the dominant modes of flux transport determined by quadrant analysis. The mean flux densities for both the cotton and grape site showed the surface acting as a sink for CO₂ and ozone and a source of heat and H₂O during the day, as would be expected, while during the night it became a source for CO₂ and a sink for heat, but remained a sink for ozone and a source of H₂O. The flux association results indicated a single vegetated ozone sink for the grape site, but a vegetated as well as a non-vegetated sink for the cotton site. For both sites, structures simultaneously transporting significant flux contributions of CO₂, H₂O, heat and ozone dominate during unstable conditions, but differed during stable conditions, where unmixed single flux structures dominated over cotton but not over grape. Structure sizes were less than 10 m during nighttime conditions and ranged from 3 to 69 m during the day. The results of this study contribute empirical evidence about the relationship between ozone uptake and the physical and physiological state of vegetation, as well as the limitations placed on eddy scales in simulation models.     

Evaluation of near surface ozone in air quality simulations forced by a regional climate model over Europe for the period 1991–2000

This study aims to evaluate near surface ozone simulated with the modelling system RegCM3/CAMx against ozone measurements from the EMEP database for the recent decade 1991–2000. The RegCM3/CAMx simulations were performed on a 50 km × 50 km grid over Europe driven either by ERA-40 reanalysis (hereafter referred as ERA simulation) or the global circulation model (GCM) ECHAM5 (hereafter referred as ECHAM simulation). A set of statistical metrics is used for the model evaluation, including temporal correlation coefficient, the ratio of the standard deviations and the bias of simulated versus observed values. Overall, a good agreement is found for both ERA and ECHAM simulations at the majority of the selected EMEP stations in all metrics throughout the year based either on monthly or daily ozone values. Based on these results, it is assessed that the modelling system RegCM3/CAMx is suitable to be used for present and future regional climate-air quality simulations with emphasis on near surface ozone. The ERA simulations reproduce more accurately the observed ozone values in comparison to ECHAM simulations because the meteorology of the ERA experiment is closer to real atmospheric conditions than the GCM based experiment. On a seasonal basis, both ERA and ECHAM simulations exhibit a seasonally dependent bias, with winter and spring ozone values being generally under-estimated by the model and summer and autumn values being slightly overestimated. This seasonally dependent bias is also evident from median and peak midday ozone values. However, the highest observed midday ozone peaks in summer, with values higher than 80 ppbv, could not be captured either by ERA or ECHAM simulations. An analysis of day-time and night-time ERA and ECHAM modelled ozone values shows that CAMx performs better during the day-time.     

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.     

Quantifying ozone uptake at the canopy level of spruce,pine and larch trees at the alpine timberline: an approach based on sap flow measurement

Micro-climatic and ambient ozone data were combined with measurements of sap flow through tree trunks in order to estimate whole-tree ozone uptake of adult Norway spruce (Picea abies), cembran pine (Pinus cembra), and European larch (Larix decidua) trees. Sap flow was monitored by means of the heat balance approach in two trees of each species during the growing season of 1998. In trees making up the stand canopy, the ozone uptake by evergreen foliages was significantly higher than by deciduous ones, when scaled to the ground area. However, if expressed per unit of whole-tree foliage area, ozone flux through the stomata into the needle mesophyll was 1.09, 1.18 and 1.40 nmol m(-2) s(-1) in Picea abies, Pinus cembra and Larix decidua, respectively. These fluxes are consistent with findings from measurements of needle gas exchange, published from the same species at the study site. It is concluded that the sap flow-based approach offers an inexpensive, spatially and temporally integrating way for estimating ozone uptake at the whole-tree and stand level, intrinsicly covering the effect of boundary layers on ozone flux.     

Quantifying simultaneous fluxes of ozone, carbon dioxide and water vapor above a subalpine forest ecosystem

Assessing the long-term exchange of trace gases and energy between terrestrial ecosystems and the atmosphere is an important priority of the current climate change research. In this regard, it is particularly significant to provide valid data on simultaneous fluxes of carbon, water vapor and pollutants over representative ecosystems. Eddy covariance measurements and model analyses of such combined fluxes over a subalpine coniferous forest in southern Wyoming (USA) are presented. While the exchange of water vapor and ozone are successfully measured by the eddy covariance system, fluxes of carbon dioxide (CO(2)) are uncertain. This is established by comparing measured fluxes with simulations produced by a detailed biophysical model (FORFLUX). The bias in CO(2) flux measurements is partially attributed to below-canopy advection caused by a complex terrain. We emphasize the difficulty of obtaining continuous long-term flux data in mountainous areas by direct measurements. Instrumental records are combined with simulation models as a feasible approach to assess seasonal and annual ecosystem exchange of carbon, water and ozone in alpine environments. The viability of this approach is demonstrated by: (1) showing the ability of the FORFLUX model to predict observed fluxes over a 9-day period in the summer of 1996; and (2) applying the model to estimate seasonal dynamics and annual totals of ozone deposition and carbon, and water vapor exchange at our study site. Estimated fluxes above this subalpine ecosystem in 1996 are: 195 g C m(-2) year(-1) net ecosystem production, 277 g C m(-2) year(-1) net primary production, 535 mm year(-1) total evapo-transpiration, 174 mm year(-1) canopy transpiration, 2.9 g m(-2) year(-1) total ozone deposition, and 1.72 g O(3) m(-2) year(-1) plant ozone uptake via leaf stomata. Given the large portion of non-stomatal ozone uptake (i.e. 41% of the total annual flux) predicted for this site, we suggest that future research of pollution-vegetation interactions should relate plant response to actively assimilated ozone by foliage rather than to total deposition. In this regard, we propose the Physiological Ozone Uptake Per Unit of Leaf Area (POUPULA) as a practical index for quantifying vegetation vulnerability to ozone damage. We estimate POUPULA to be 0.614 g O(3) m(-2) leaf area year(-1) at our subalpine site in 1996.     

The influence of the photolysis rates on modelled ozone concentrations

The sensitivity of ozone formation in a photochemical model to the photolysis rates has been studied using two photolytic schemes (old EMEP and UKPTM). The results from the comparison of these rates illustrate how the different parameterisation may influence the ozone development. The implementation of the UKPTM photolysis rates into the EMEP chemical scheme resulted in lower ozone levels than those predicted using the old EMEP photolytic scheme.     

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.     

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.     

Summer and autumn ozone fluxes to a forest in the Czech Republic Brdy Mountains

In an effort to examine ozone (O3) deposition over a forest site in the Czech Republic, a low cost eddy flux experiment using slow response ozone and temperature sensors was implemented in July 1993 within the Brdy Mountains. Half-hour 2-Hz ozone and sensible heat measurements made at the Brdy Mountains for 98 days during the period 7 July 1994-20 October 1994 are analyzed and reported. While the Czech Brdy Mountains AOT40 level for the overall 104 day period was 7.6 ppm h (15.1 ppm h for the full 24-h summation), indicating a slight potential for 03 injury, the 1994 summer to autumn'measured forest O3 uptake was 2.4 (+/- 0.9) g m(-2), not unusually high compared to other studies. Average summer midday 03 fluxes and depositidn velocities were -1.0 (+/- 0.6) microg m(-2) s(-1) and 1.1 (+/- 0.7) cm s(-1). and autumn values were -0.36 (+/- 0.4) microg m(-2) s(-1) and 0.7 (+/- 0.5) cm s(-1) respectively. A unique contribution of this study is the first time demonstrated use of slow responding sensors for eddy covariance flux measurements at heights of 20 m above a forest.     

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.     

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.     

A comparison of measured and modeled ozone uptake into plant leaves

Ozone uptake into plant leaves was measured in gas exchange chambers using a mass balance and a variable conductance approach. The variable conductance approach was found to more reliably measure ozone flux through stomata. Measurements using this approach were contrasted with estimates obtained by measuring stomatal conductance g(sw) and modeling ozone uptake using a diffusion equation, assuming a negligible ozone concentration in the substomatal cavity. Actual measurements of uptake were close, but slightly higher than modeled values, providing some support to the idea that substomatal ozone concentrations are close to zero. However, the difference between measured and modeled uptake values suggests either that (i) variable conductance approach measures more ozone uptake than caused by stomatal uptake alone or (ii) ozone conductance is underestimated.     

Dynamic evaluation of a regional air quality model: Assessing the emissions-induced weekly ozone cycle

Air quality models are used to predict changes in pollutant concentrations resulting from envisioned emission control policies. Recognizing the need to assess the credibility of air quality models in a policy-relevant context, we perform a dynamic evaluation of the Community Multiscale Air Quality (CMAQ) modeling system for the “weekend ozone effect” to determine if observed changes in ozone due to weekday-to-weekend (WDWE) reductions in precursor emissions can be accurately simulated. The weekend ozone effect offers a unique opportunity for dynamic evaluation, as it is a widely documented phenomenon that has persisted since the 1970s. In many urban areas of the Unites States, higher ozone has been observed on weekends than weekdays, despite dramatically reduced emissions of ozone precursors (nitrogen oxides [NO_x] and volatile organic compounds [VOCs]) on weekends. More recent measurements, however, suggest shifts in the spatial extent or reductions in WDWE ozone differences. Using 18 years (1988–2005) of observed and modeled ozone and temperature data across the northeastern United States, we re-examine the long-term trends in the weekend effect and confounding factors that may be complicating the interpretation of this trend and explore whether CMAQ can replicate the temporal features of the observed weekend effect. The amplitudes of the weekly ozone cycle have decreased during the 18-year period in our study domain, but the year-to-year variability in weekend minus weekday (WEWD) ozone amplitudes is quite large. Inter-annual variability in meteorology appears to influence WEWD differences in ozone, as well as WEWD differences in VOC and NO_x emissions. Because of the large inter-annual variability, modeling strategies using a single episode lasting a few days or a few episodes in a given year may not capture the WEWD signal that exists over longer time periods. The CMAQ model showed skill in predicting the absolute values of ozone concentrations during the daytime. However, early morning NO_x concentrations were underestimated and ozone levels were overestimated. Also, the modeled response of ozone to WEWD differences in emissions was somewhat less than that observed. This study reveals that model performance may be improved by (1) properly estimating mobile source NO_x emissions and their temporal distributions, especially for diesel vehicles; (2) reducing the grid cell size in the lowest layer of CMAQ; and, (3) using time-dependent and more realistic boundary conditions for the CMAQ simulations.     

A theoretical investigation of the pressure and temperature dependence of atmospheric ozone deposition to trees

Methods for describing the exposure patterns of forests to atmospheric ozone concentrations are compared with special emphasis on the situation at high altitudes, such as the Appalachian Mountains of the eastern USA. Limitations to the use of ozone concentration as mass per unit volume are discussed and a correction for temperature and pressure changes is derived. If identical ozone mass concentrations were measured at two sites separated by 2000 m elevation, the ozone flux at the lower site would exceed the flux at the higher site by 4-8% due to temperature and pressure effects on both air volume and ozone deposition velocity. It is recommended that ozone exposures be described in terms of 'flux-corrected' mass concentrations or volumetric mixing ratios when ambient ozone data from sites at different altitudes are to be compared.     

Operational ozone forecasts for the region of Copenhagen by the Danish Meteorological Institute

The Danish Meteorological Institute (DMI) has developed an operational forecasting system for ozone concentrations in the Atmospheric Boundary Layer; this system is called the Danish Atmospheric Chemistry FOrecasting System (DACFOS). At specific sites where real-time ozone concentration measurements are available, a statistical after-treatment of DACFOS’ results adjusts the next 48 h ozone forecasts. This post-processing of DACFOS’ forecasts is based on an adaptive linear regression model using an optimal state estimator algorithm. The regression analysis uses different linear combinations of meteorological parameters (such as temperature, wind speed, surface heat flux and atmospheric boundary layer height) supplied by the Numerical Weather Prediction model DMI-HIRLAM. Several regressions have been tested for six monitoring stations in Denmark and in England, and four of the linear combinations have been selected to be employed in an automatic forecasting system. A statistical study comparing observations and forecasts shows that this system yields higher correlation coefficients as well as smaller biases and RMSE values than DACFOS; the present post-processing thus improves DACFOS’ forecasts. This system has been operational since June 1998 at the DMI's monitoring station in the north of Copenhagen, for which a new ozone forecast is presented every 6 h on the DMI's internet public homepage.