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.     

Effect of climate change on fluxes of nitrogen from the Tovdal River basin, Norway, to adjoining marine areas

The mass transport model TEOTIL was used to project nitrate (NO3) fluxes from the Tovdal River basin, southernmost Norway, given four scenarios of climate change. Forests, uplands, and open water currently account for 90% of the NO3 flux. Climate scenarios for 2071-2100 suggest increased temperature by 2-4 degrees C and precipitation by 3-11%. Climate experiments and long-term monitoring were used to estimate future rates of nitrogen (N) leaching. More water will run through the terrestrial catchments during the winter but less will run in the spring. The annual NO3 flux from the Tovdal River to the adjoining Topdalsfjord is projected to remain unchanged, but with more NO3 delivered in the winter and less in the spring. Algal blooms in coastal waters can be expected to occur earlier in the year. Major sources of uncertainty are in the long-term fate of N stored in soil organic matter and the impacts of forest management.     

Ammonia emission factors for N fertilizers applied to two contrasting grassland soils

Ammonia volatilization from nitrogen (N) fertilizer applied throughout the year to two soil types was measured using a system of small wind tunnels. Losses from urea ranged from 12 to 46% of the applied N. Small losses, averaging <1%, were measured from ammonium nitrate (AN) and calcium nitrate applications. Factors influencing these losses are discussed. Using these results and those from other workers, emission factors for urea and AN applications to grassland in the UK were determined as 23.0 and 1.6% of the applied N, respectively. Emission factors for these fertilizers when applied to arable land were estimated as 11.8 and 0.8%, respectively. The emission factor for all other applied N (as straight and compound fertilizers) was assumed to be similar to that for AN. Calculations showed that fertilizer applications to agricultural land in the UK contributes 34 kt NH3-N per year, equivalent to 17% of the total annual NH3 emission.     

Predicting emissions from oil and gas operations in the Uinta Basin,Utah

In this study, emissions of ozone precursors from oil and gas operations in Utah’s Uinta Basin are predicted (with uncertainty estimates) from 2015–2019 using a Monte-Carlo model of (a) drilling and production activity, and (b) emission factors. Cross-validation tests against actual drilling and production data from 2010–2014 show that the model can accurately predict both types of activities, returning median results that are within 5% of actual values for drilling, 0.1% for oil production, and 4% for gas production. A variety of one-time (drilling) and ongoing (oil and gas production) emission factors for greenhouse gases, methane, and volatile organic compounds (VOCs) are applied to the predicted oil and gas operations. Based on the range of emission factor values reported in the literature, emissions from well completions are the most significant source of emissions, followed by gas transmission and production. We estimate that the annual average VOC emissions rate for the oil and gas industry over the 2010–2015 time period was 44.2E+06 (mean) ± 12.8E+06 (standard deviation) kg VOCs per year (with all applicable emissions reductions). On the same basis, over the 2015–2019 period annual average VOC emissions from oil and gas operations are expected to drop 45% to 24.2E+06 ± 3.43E+06 kg VOCs per year, due to decreases in drilling activity and tighter emission standards.

Implications: This study improves upon previous methods for estimating emissions of ozone precursors from oil and gas operations in Utah’s Uinta Basin by tracking one-time and ongoing emission events on a well-by-well basis. The proposed method has proven highly accurate at predicting drilling and production activity and includes uncertainty estimates to describe the range of potential emissions inventory outcomes. If similar input data are available in other oil and gas producing regions, then the method developed here could be applied to those regions as well.     

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.     

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.     

Impact of biogenic emission uncertainties on the simulated response of ozone and fine particulate matter to anthropogenic emission reductions

The role of emissions of volatile organic compounds and nitric oxide from biogenic sources is becoming increasingly important in regulatory air quality modeling as levels of anthropogenic emissions continue to decrease and stricter health-based air quality standards are being adopted. However, considerable uncertainties still exist in the current estimation methodologies for biogenic emissions. The impact of these uncertainties on ozone and fine particulate matter (PM2.5) levels for the eastern United States was studied, focusing on biogenic emissions estimates from two commonly used biogenic emission models, the Model of Emissions of Gases and Aerosols from Nature (MEGAN) and the Biogenic Emissions Inventory System (BEIS). Photochemical grid modeling simulations were performed for two scenarios: one reflecting present day conditions and the other reflecting a hypothetical future year with reductions in emissions of anthropogenic oxides of nitrogen (NOx). For ozone, the use of MEGAN emissions resulted in a higher ozone response to hypothetical anthropogenic NOx emission reductions compared with BEIS. Applying the current U.S. Environmental Protection Agency guidance on regulatory air quality modeling in conjunction with typical maximum ozone concentrations, the differences in estimated future year ozone design values (DVF) stemming from differences in biogenic emissions estimates were on the order of 4 parts per billion (ppb), corresponding to approximately 5% of the daily maximum 8-hr ozone National Ambient Air Quality Standard (NAAQS) of 75 ppb. For PM2.5, the differences were 0.1-0.25 microg/m3 in the summer total organic mass component of DVFs, corresponding to approximately 1-2% of the value of the annual PM2.5 NAAQS of 15 microg/m3. Spatial variations in the ozone and PM2.5 differences also reveal that the impacts of different biogenic emission estimates on ozone and PM2.5 levels are dependent on ambient levels of anthropogenic emissions.     

Carbon dioxide fluxes over a grazed prairie and seeded pasture in the Northern Great Plains

Temperate grasslands are vast terrestrial ecosystems that may be an important component of the global carbon (C) cycle; however, annual C flux data for these grasslands are limited. The Bowen ratio/energy balance (BREB) technique was used to measure CO2 fluxes over a grazed mixed-grass prairie and a seeded western wheatgrass [Pascopyrum smithii (Rybd) L?ve] site at Mandan, ND from 24 April to 26 October in 1996, 1997, and 1998. Above-ground biomass and leaf area index (LAI) were measured about every 21 days throughout the season. Root biomass and soil organic C and N content were determined to 110 cm depth in selected increments about mid-July each year. Peak above-ground biomass and LAI coincided with peak fluxes and occurred between mid-July to early August. Biomass averaged 1227 and 1726 kg ha(-1) and LAI 0.44 and 0.59, for prairie and western wheatgrass, respectively. Average CO2 flux for the growing season was 279 g CO2 m(-2) for prairie and 218 g CO2 m(-2) for western wheatgrass (positive flux is CO2 uptake and negative flux is CO2 loss to the atmosphere). Using prior measured dormant season CO2 fluxes from the prairie sites gave annual flux estimates that ranged from -131 to 128 g CO2 m(-2) for western wheatgrass and from -70 to 189 g CO2 m(-2) for the prairie. This wide range in calculated annual fluxes suggests that additional research is required concerning dormant season flux measurements to obtain accurate estimates of annual CO2 fluxes. These results suggest Northern Great Plains mixed-grass prairie grasslands can either be a sink or a source for atmospheric CO2 or near equilibrium, depending on the magnitude of the dormant season flux.     

A comparison of two different approaches for mapping potential ozone damage to vegetation. A model study

Two very different types of approaches are currently in use today for indicating risk of ozone damage to vegetation in Europe. One approach is the so-called AOTX (accumulated exposure over threshold of Xppb) index, which is based upon ozone concentrations only. The second type of approach entails an estimate of the amount of ozone entering via the stomates of vegetation, the AFstY approach (accumulated stomatal flux over threshold of Y nmol m(-2) s(-1)). The EMEP chemical transport model is used to map these different indicators of ozone damage across Europe, for two illustrative vegetation types, wheat and beech forests. The results show that exceedences of critical levels for either type of indicator are widespread, but that the indicators give very different spatial patterns across Europe. Model simulations for year 2020 scenarios suggest reductions in risks of vegetation damage whichever indicator is used, but suggest that AOT40 is much more sensitive to emission control than AFstY values.     

Meteorologically adjusted trends of daily maximum ozone concentrations in Taipei,Taiwan

The meteorological conditions exert large impacts on ozone concentrations, and may mask the long-term trends in ozone concentrations resulting from precursor emissions. Estimation of long-term trends of ozone concentrations due to the changes in precursor emissions is important for corresponding control strategy. Multiple linear regression (method I), multilayer perceptron (MLP) neural network (method II) and Komogorov-Zurbenko (KZ) filter method plus MLP methodology (method III), are used to estimate the meteorologically adjusted long-term trends of daily maximum ozone concentrations by removing the masking effects of meteorological conditions in this study. The daily maximum ozone concentrations and relative meteorological variables were extracted from six air-monitoring stations in Taipei area from 1994 to 2001. The data collected during 1994–2000 period were used as modeling set and utilized to estimate the meteorologically adjusted trends, and the data of 2001 were used as the validation data. The meteorologically adjusted trends of ozone for these three methods were calculated and compared. The results show that both MLP and KZ filter +MLP models are more suitable than multiple linear regression for estimating the long-term trends of ozone in Taipei, Taiwan. The long-term linear trends of meteorologically adjusted ozone concentrations due to the precursor emissions show an increase trend at all stations, and the percent changes per year range from 1.0% to 2.25% during the modeling period in Taipei area.     

Stress responses of Calluna vulgaris to reduced and oxidised N applied under 'real world conditions'

Effects and implications of reduced and oxidised N, applied under 'real world' conditions, since May 2002, are reported for Calluna growing on an ombrotrophic bog. Ammonia has been released from a 10 m line source generating monthly concentrations of 180-6 microg m(-3), while ammonium chloride and sodium nitrate are applied in rainwater at nitrate and ammonium concentrations below 4mM and providing up to 56 kg N ha(-1) year(-1) above a background deposition of 10 kg N ha(-1) year(-1). Ammonia concentrations, >8 microg m(-3) have significantly enhanced foliar N concentrations, increased sensitivity to drought, frost and winter desiccation, spring frost damage and increased the incidence of pathogen outbreaks. The mature Calluna bushes nearest the NH3 source have turned bleached and moribund. By comparison the Calluna receiving reduced and oxidised N in rain has shown no significant visible or stress related effects with no significant increase in N status.     

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.     

Halogenated compounds and climate change: future emission levels and reduction costs

OBJECTIVES: This work assesses the contribution to climate change resulting from emissions of the group of halogenated greenhouse gases. METHODS: A bottom-up emission model covering 22 technological sectors in four major regions is described. Emission estimates for 1996 and projection for 2010 and 2020 are presented. The costs for deep cuts into projected emission levels are calculated. RESULTS: The substances covered by this study have contributed emissions of 1100 +/- 800 MT CO2 equivalents per year in 1996. In terms of their relative contribution to emissions of CO2 equivalents, this corresponds to 3 +/- 2% of global emissions of all anthropogenic greenhouse gases. The wide range of uncertainty is due to the poorly quantified net global warming potential of the ozone depleting substances, which have an indirect cooling effect on climate through the destruction of stratospheric ozone. For annual emissions of HFCs, PFCs and SF6 (which are regulated under the Kyoto Protocol and for which global warming potentials are well defined), the relative contribution is projected to increase to 2% (600 MT CO2 eq.) of global greenhouse gas emissions by 2010. This trend is expected to continue, emissions are projected to grow to a contribution of roughly 3% (870 MT CO2 eq.) in 2020 compared to 0.9% (300 MT CO2 eq.) in 1996. For HFCs, PFCs and SF6, this study identifies global emission reduction potentials of 260 MT CO2 eq. per year in 2010 and 640 MT CO2 eq. per year in 2020 at below US$ 50 per ton. These values correspond to roughly 40% and 75% of projected emissions in 2010 and 2020, respectively.     

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.     

Use of (15)N-labelled nitrogen deposition to quantify the source of nitrogen in runoff at a coniferous-forested catchment at Gårdsjön, Sweden

To determine the source of dissolved inorganic nitrogen (N) in runoff, approx. 35kg N enriched with the stable isotope (15)N (2110 per thousand delta(15)N) was added to a mature coniferous forested catchment for one whole year. The total N input was approx. 50kg ha(-1) year(-1). The enrichment study was part of a long-term whole-catchment ammonium nitrate addition experiment at G?rdsj?n, Sweden. The (15)N concentrations in precipitation, throughfall, runoff and upper forest floor were measured prior to, during, and 3-9years following the (15)N addition. During the year of the (15)N addition the delta(15)N level in runoff largely reflected the level in incoming N, indicating that the leached NO(3)(-) came predominantly from precipitation. Only 1.1% of the incoming N was lost during the year of the tracer addition. The cumulative loss of tracer N over a 10-year period was only 3.9% as DIN and 1.1% as DON.     

Boreal bog plants: nitrogen sources and uptake of recently deposited nitrogen

The main goals of this study were to determine the delta15N signature of quantitatively important boreal bog plants as basis for discussing their N sources, and to assess the effects of five different 3 year N treatments (i.e. 0, 5, 10, 20 and 40 kg N ha(-1) year(-1)) on the bog plants and surface peat at different depths (i.e. 0, 5, 10, 20 and 40 cm) by using 15N as tracer. Plants and peat were analyzed for N concentration, 15N natural abundance and 15N at.%. From the results we draw three main conclusions: First, the relative importance of different N sources is species-specific among bog plants. Second, an annual addition of 5 kg N ha(-1) year(-1) was sufficient to significantly increase the N concentration in Sphagnum mosses, liverworts and shallow rooted vascular plants, and an annual addition of 40 kg N ha(-1) year(-1) during 3 years was not sufficient to increase the N concentration in deep rooted plants, although the 15N content increased continuously, indicating a possible longer term effect. Third, an annual addition of 40 kg N ha(-1) year(-1) during 3 years increased the N content in surface peat at depths of 5 and 10 cm, but not at depths of 20 and 40 cm, indicating the capacity of the living Sphagnum mosses and the surface peat to take up deposited N, and thereby function as a filter.     

The effects of excess nitrogen deposition on young Norway spruce trees. Part I the soil

The effects of wet-deposited nitrogen on soil acidification and the health of Norway spruce were investigated in a pot experiment using an open-air spray/drip system. Nitrogen was applied as ammonium ((NH(4))(2)SO(4)) or nitrate (HNO(3)/NaNO(3)) in simulated rain to either the soil or the foliage for a period of two years five months. Symptoms of forest decline were not reproduced. Adverse effects relating to soil acidification and N saturation were observed and depended on the chemical form of N. The plant-soil system absorbed most of the soil-applied NH(+)(4) at doses of up to 65 kgN ha(-1) year(-1) but only 54% at a dose of 125 kgN ha(-1) year(-1). About 60% of soil-applied NO(-)(3) was absorbed in all treatments. Ammonium treatments acidified the soil, increased base cation leaching, and mobilised acidic cations. Nitrification was not the major source of acidity, however. Nitrate inputs increased soil pH. Critical loads calculated using current criteria were 60-120 and 30-60 kgN ha(-1) year(-1) for ammonium and nitrate, respectively. Ammonium is more likely to damage forest ecosystems, however, illustrating the need for care in the definition of critical loads.     

Ozone and UV-B variations at Ispra from 1993 to 1997

An analysis of the variability of the total ozone column at Ispra (Italy) has been performed to ascertain if, even in a short-time interval of 5 years (1993–1997), a decline of the monthly mean ozone values could be demonstrated. A linear fit of the data displays a decrease of 0.21% per year with a mean value equal to 319±2 D.U. and an amplitude of the annual cycle of about 10% of the mean. A linear regression of the surface monthly mean ozone values has also been performed showing a decreasing trend (−1% per year) that could contribute, even if for a very small amount, to the decline of the total ozone values. Ispra monthly mean total ozone data have been compared with those of three stations located within 2° latitude and 3° longitude from Ispra (Haute Provence, Hohenpeissenberg and Arosa). A linear fit of the data shows some discrepancies in the ozone changes, which can be attributed to the limited length of the observational period.An analysis has been performed to verify if the variation of ozone at Ispra is in agreement with that of the solar UV measured at a wavelength (305 nm) where the ozone absorption is still remarkable. The results, taken at a fixed solar zenith angle of 68°, show a clear anticorrelation between the monthly mean values of UV and the corresponding values of the total ozone column; the linear fit of the UV data displays an increase of 2.0% per year, much higher than expected from the ozone decrease, and a mean value of 1.4±0.1 mW m^-2 nm^-1.     

Long-term trends and spatial variability in nitrate leaching from alpine catchment-lake ecosystems in the Tatra Mountains (Slovakia-Poland)

Relationships between catchment characteristics of 31 alpine lakes and observed trends in lake water concentrations of nitrate were evaluated in the Tatra Mountains. Nitrate concentrations increased from background levels <4 microeql(-1) in the 1930s to maxima (up to 55 microeql(-1)) in the 1980s, after which they declined to 4-44 microeql(-1) by the late 1990s. In-lake nitrate concentrations correlated negatively with parameters characterising catchment-weighted mean pools (CWM; kgm(-2)) of soil, i.e. with percent land cover with meadow and soil depth, and positively with grade of terrain, annual precipitation, and the highest elevation in the catchment. The CWM pool of soil and annual precipitation explained together 65% of the current spatial variability in nitrate concentrations. Denitrification and direct N deposition on surface area explained 14% of the variability. Increased atmospheric N deposition and declining net N retention in soils were responsible for long-term changes in nitrate concentrations. Long-term decline in %N retention in soils decreased along with the estimated decline in C:N ratios (from 21 to 18 on average during the last 70 years). An empirical model linking nitrate concentrations in different types of alpine Tatra Mountain lakes to four independent variables (CWM soil pool, annual precipitation, increased N deposition, and average trend in soil C:N ratios) explained 80% of the observed spatial and temporal nitrate variability over the period 1937-2000.     

Increasing interannual and altitudinal ozone mixing ratios in the Catalan Pyrenees

Interannual, seasonal, daily and altitudinal patterns of tropospheric ozone mixing ratios, as well as ozone phytotoxicity and the relationship with NO_x precursors and meteorological variables were monitored in the Central Catalan Pyrenees (Meranges valley and Forest of Guils) over a period of 5 years (2004–2008). Biweekly measurements using Radiello passive samplers were taken along two altitudinal transects comprised of thirteen stations ranging from 1040 to 2300 m a.s.l. Visual symptoms of ozone damage in Bel-W3 tobacco cultivars were evaluated biweekly for the first three years (2004–2006). High ozone mixing ratios, always above forest and vegetation protection AOT40 thresholds, were monitored every year. In the last 14 years, the AOT40 (Apr–Sept.) has increased significantly by 1047 μg m^?3 h per year. Annual means of ozone mixing ratios ranged between 38 and 67 ppb_v (38 and 74 ppb_v during the warm period) at the highest site (2300 m) and increased at a rate of 5.1 ppb_v year^?1. The ozone mixing ratios were also on average 35–38% greater during the warm period and had a characteristic daily pattern with minimum values in the early morning, a rise during the morning and a decline overnight, that was less marked the higher the altitude. Whereas ozone mixing ratios increased significantly with altitude from 35 ppb_v at 1040 m–56 ppb_v at 2300 m (on average for 2004–2007 period), NO₂ mixing ratios decreased with altitude from 5.5 ppb_v at 1040 m–1 ppb_v at 2300 m. The analysis of meteorological variables and NO_x values suggests that the ozone mainly originated from urban areas and was transported to high-mountain sites, remaining aloft in absence of NO. Ozone damage rates increased with altitude in response to increasing O₃ mixing ratios and a possible increase in O₃ uptake due to more favorable microclimatic conditions found at higher altitude, which confirms Bel-W3 as a suitable biomonitor for ozone concentrations during summer time. Compared to the valley-bottom site the annual means of ozone mixing ratios are 37% larger in the higher sites. Thus the AOT40 for the forest and vegetation protection threshold is greatly exceeded at higher sites. This could have substantial effects on plant life at high altitudes in the Pyrenees.