Introduction

Trends in environmental conditions in Europe

Environmental conditions are rapidly changing on a global scale since the beginning of the industrial times (Schöpp et al. 2003; IPCC 2007). The human development has not only altered the chemical composition of the atmosphere, but also of water bodies, soils, flora and fauna. Since the late 19th century until its peak in the mid-1980s, sulphur (S) emissions in Europe increased by a factor of seven, while nitrogen (N) emissions increased by a factor near five, correlated with a similar increase in acid deposition (Schöpp et al. 2003). The initial concern for negative effects of acid deposition started in Scandinavia in the seventies, in particular in south Norway, where slow weathering bedrocks and shallow soils made the ecosystems sensitive, and widespread death of trout and salmon was attributed to the enhanced leaking of aluminium (Al) which damaged gills (Baker and Schofield 1982). Low pH and increased Al levels cause chronic stress can also lead to lower body weight and smaller size and makes fish less able to compete for food and habitat (Lien et al. 1996). Shortly afterwards, in the beginning of the eighties, the health and vitality of forest ecosystems became a subject of wide public and political concern due to the extensive forest damage observed in rural areas in Central Europe, which was connected to air pollution and acid deposition (e.g Schütt et al. 1983; Lammel 1984). The extensive forest damage observed in rural areas in Central Europe in the beginning of the 1980’s (e.g. Schütt et al. 1983; Lammel 1984) was in particular the case for the Ore Mountains on the border between Germany, Poland and the Czech Republic, which were heavily exposed to sulphur dioxide (SO2).

Of more recent concern are the increasing concentrations of tropospheric ozone (O3), which is not only a greenhouse gas with the third strongest radiative forcing on climate (Forster et al. 2007), but is also the air pollutant considered to be causing the most damage to plants in Europe and the US today (Ashmore 2005; Karnosky et al. 2007a; Matyssek et al. 2007b). Tropospheric O3 pollution has shifted from a regional to a global issue because of its intercontinental transport (Derwent et al. 2004; Vingarzan 2004). In fact, tropospheric O3 has been recognized as an important factor within “global change” (IPCC 2001; Ashmore 2005), that has the capacity of reducing carbon sink strength of forest ecosystems (e.g. Karnosky et al. 2003; Sitch et al. 2007; Pretzsch et al. 2010) and modify metabolic responses under elevated atmospheric CO2 (Karnosky et al. 2007b; Matyssek et al. 2010a). Emissions associated with fossil fuel and biomass burning have acted to approximately double the global mean tropospheric O3 concentration (Gauss et al. 2006), and further increases are expected over the twenty-first century (Gauss et al. 2003; IPCC 2007). However, major uncertainties in modelled O3 response exist, mainly due to changes in methane (Wild et al. 2012). Trend analyses of O3 measurements are mainly restricted to northern and western parts of Europe. Overviews of reported trends are given in several publications, such as Monks et al. (2003) and Mol et al. (2010). The temporal trends in tropospheric O3 do not show a uniform picture, apart from sites at the western coast of Europe which are measuring background values. Most sites show a substantial downward trend in peak O3 concentrations (98 and 95 percentiles) over the past 10 to 15 years, which is due to European abatement on nitrogen oxides and volatile organic compound emissions. However, despite this downward trend in peak O3 concentrations, there has been a slight increase in the accumulated dose over a threshold of 40 ppb (AOT40; Fuhrer et al. 1997) in rural stations (Fiala et al. 2003).

Another more recent concern is the rise in CO2 concentrations by 15 % since the late 19th century. Over the 20th century, this has been accompanied by an increase in global average surface temperature by 0.6 ± 0.2 °C (IPCC 2007) and it is expected to increase by 1.4 to 5.8 °C until 2100 (IPCC 2007). While past temperature change was most pronounced in winter months, current temperature increases in Europe also include spring and summer temperatures (Rebetez and Reinhard 2008). In the European Alps, temperature has increased by more than twice the global average in the 20th century (Christensen and Christensen 2003). Most of this increase has occurred during the last 20 years (Rebetez and Reinhard 2008). For future precipitation, an increase in extreme events is predicted (Christensen and Christensen 2003; Frei et al. 2006), while more hot spells are expected in summer due to increasing variability in temperature (Schär et al. 2004).

A comparison of estimated trends in S and N emission, CO2 concentration and temperature for Europe between 1880 and 2005 is shown in Fig. 1. It illustrates why S emission was of highest concern in the early 1980s. By that time, S emission had increased by 700 % of the pre-industrial time and total N emission by roughly 540 %. CO2 concentrations had risen by 18 % and temperature by 0.4 °C. Since then, S and N emissions have declined whereas CO2 and temperature is constantly rising, causing an increasing interest in climate change impacts on forest ecosystems.

Fig. 1
figure 1

Relative changes in: (i) SO2 and N emissions in Europe between 1880 and 2007 (data until 1990 based on (Schöpp et al. 2003) and EEA data from 1990–2007), (ii) CO2 concentrations since 1960 (station Mauna Loa) in comparison to estimated values from 1880 1880 base 292 ppm and (iii) mean annual temperature deviation with respect to the 1961–1990 reference period for the Northern hemisphere (University of East Anglia data set) averaged using a 10-year moving window. The figure also shows the ICP Level and Level II forest monitoring networks in Europe (big grey bars), that started in 1986 and 1994

Impacts of air quality and acid deposition on forest ecosystems

Forest condition is influenced by a multitude of stress factors, including air pollution and climatic factors. The hypotheses on air pollution, which have been subject to considerable controversy over the years (Skelly and Innes 1994), include direct, above-ground impacts of sulphur dioxide (SO2) and ozone (O3) on foliage and indirect, soil-mediated impacts of N and S deposition on roots that may cause nutrient deficiencies and aggravate natural stress, such as physiological drought and the occurrence of pests and diseases (e.g. Landmann and Bonneau 1995). The relative contribution of a given factor may vary on its spatial and temporal distribution, i.e. between regions in Europe and from year to year.

Especially in the eighties and nineties of the previous century, there has been a large research effort to gain insight in the mechanisms behind air-pollution impacts on forest condition by a combination of experiments, monitoring and modelling. Based on this research, the following air-pollution related impacts on forest condition and forest growth have been hypothesized (see e.g. De Vries et al. 2000b):

  • Elevated concentrations of the gases SO2 and O3, causing (i) leaf stomata disturbances and premature senescence, leading to water regulation stress and physiological drought, (ii) reduction or shift of carbon allocation leading to a weakened root system and (iii) accelerated foliar leaching affecting the nutrient status.

  • Eutrophication by elevated N (NO3, NH4) deposition inputs, causing (i) a shift in deficiency from N to base cation nutrients (Bc), where Bc stands for the sum of calcium (Ca), magnesium (Mg) and potassium (K), due to elevated demands induced by an initial growth increase, (ii) enhanced drought stress since an elevated N input favours growth of canopy biomass, whereas root growth is relatively unaffected and (iii) an increased sensitivity to natural stresses, such as frost and fungal diseases.

  • Soil acidification by S (SO4) and N deposition, including (i) the loss of Bc from the soil causing deficiency of these nutrients (notably Mg), (ii) the release of toxic Al affecting fine root growth and inhibiting Bc uptake and (iii) a decrease in pH that may increase the mobility of heavy metals.

Elevated atmospheric N deposition can also cause changes in plant species composition and the following main mechanisms for effects are considered (Bobbink et al. 2010):

  • Direct toxicity of nitrogen gases and aerosols may result in changes in physiology and reduce growth at high concentrations,

  • Increased availability of nitrogen may result in an increase in plant productivity, litter production and N mineralisation, and thus cause competitive exclusion of characteristic species by fast-growing nitrophilic species,

  • High concentrations of ammonium and ammonia in the soil solution may be toxic to sensitive plant species, especially in habitats with moderately acidic conditions

  • N deposition may result in soil acidification causing acid-resistant plant species to become dominant and species typical of intermediate pH to disappear.

  • High availability of nitrogen may increase susceptibility to secondary stress and disturbance factors, such as plant pathogens and insect pests.

At present, the concern is related to the (long-term) resilience of forest ecosystems to cope with the interacting impacts of air pollution and climate change, and the possibility to counteract these impacts by forest management. Below, we discuss effects of N and acid deposition, O3 exposure, and the impacts of climatic factors.

Aims of the large-scale forest monitoring networks in Europe

As a consequence of the concern of acid deposition effects, a large-scale monitoring of forest condition began under the umbrella of the United Nations Economic Commission for Europe (UN-ECE) and EU resulting in the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests) network in the mid-1980s (Lorenz 1995) and a related EU Scheme on the Protection of Forests against Atmospheric Pollution. ICP-Forests is part of the UN-ECE Long-Range Transboundary Air Pollution (LRTAP) convention. From 1986 onwards, data on tree crown condition, i.e. defoliation and discoloration, were gathered annually at forested plots that were placed in a systematic fashion (in most countries on a 16 × 16 km grid), following harmonized methods and centrally stored (e.g. Fischer et al. (2010; e.g. EC-UN/ECE 2007). The major aim of this so-called Level I Monitoring Programme was to gain insight in the geographic and temporal variations in forest condition and its possible relationship with stress factors, including air pollution. As it became more and more clear that large-scale monitoring was not able to answer questions regarding the cause-effect relationships, a “Pan- European Programme for Intensive and Continuous Monitoring of Forest Ecosystems” (the so-called Level II Monitoring Programme) started in 1994 under the umbrella of the UN-ECE/ICP Forests and EU, with currently approximately 760 permanent observation plots, located in 30 participating countries, including 500 plots with data on both atmospheric deposition and forests ecosystem impacts (De Vries et al. 2003c; Lorenz and Becher 2012). In the sequel both networks are referred to as the ICP Forests Level I and Level II networks. An overview of the measurements carried out and the methods used in both monitoring networks is given in Ferretti and Fischer (2013).

The ICP Forest monitoring programmes initially focused on air pollution impacts on forest nutrition and forest health, since the results should be useful for the validation and further development of protocols on air pollution control strategies. However, over time, natural stress factors and climate change gained increasing attention, including impacts on forest growth and plant species diversity of the forest ground vegetation. This change in focus follows from the changes in environmental conditions shown in Fig. 1. In 1985, when the time that the first large-scale forest health monitoring began at so-called Level I plots, S and N emissions and related acidification effects were top priority in forest impact research. However, since then, S emissions have declined by 75 % to only 160 % of pre-industrial times, while total N emissions have only declined by 25 % and are still roughly four times as high as in the 19th century. CO2 rise had slightly increased from between 0.5 ppm per year until 1985 to 1.7 ppm per year between 1980 and 2007. Temperature has risen more dramatically with an increase of almost 0.8 °C since 1980 (double the value of the previous 100 years). Thus, during the short period of large-scale forest monitoring in Europe, environmental conditions have drastically changed, illustrating the increasing interest in the combined interacting impacts of N deposition, NOx induced O3 exposure and climate change.

Whereas the ICP Forests Level I network predominantly aims to gain insight in the geographic and temporal variations in forest condition, the aim of the ICP Forests Level II network is to clarify cause-effect relationships. The ICP Forests Level II network includes 15 different surveys, with an assessment of crown condition and damages, tree stem increment and the chemical composition of foliage and soil on all plots, whereas atmospheric deposition, meteorological parameters, soil solution chemistry, ground vegetation composition litterfall, O3 injury, ambient air quality, ground vegetation biomass, leaf area index and soil water are conducted on a number of selected permanent observation plots spread over Europe (Ferretti et al. 2013). The Level II network is now following trends in stresses and responses for approximately 760 sites spread over 40 European countries, with 15-year time series. These data allow an evaluation of the impacts due to changes in weather conditions, including extremes, O3 exposures and deposition of atmospheric S and N compounds.

Aim and contents of this paper

In this paper, field evidence is presented of various impacts of air pollution and climate change on forest ecosystems based on monitoring data from the ICP Forests Level I and Level II networks, with a focus on the ICP Forests Intensive Monitoring Level II network. First, a summary is given of impacts of air pollution and climate on forests ecosystems as reported in the literature on the basis of laboratory and field research, including hypotheses on the overall interacting impacts of air pollution and climate change. Next, results from monitoring studies, both at a European wide scale and related national studies, are presented in terms of measured trends and geographic variation in N and S deposition and O3 concentrations/exposure, followed by results on effects of elevated N and S inputs and/or O3 exposure and/or climate on:

  • Nutrient status: element budgets, soil solution chemistry and foliar chemistry

  • Forest health in terms of crown condition and crown transparency

  • Forest growth and related tree carbon sequestration

  • Species diversity of the ground vegetation

We conclude with an evaluation of the monitoring results in view of the knowledge based on results of laboratory and field studies. Rather than giving an in-depth review of the interactions between air quality and climate on forest ecosystems, we present a broad overview of empirical, field data based, forest responses to many drivers.

Impacts as derived from laboratory and field studies

An overview of the potential impacts of air quality and climate on forest ecosystems in terms of the soil and solution quality, nutrition, condition and growth of forests (trees) and the species composition of forest ground vegetation is given in Table 1. The various impacts are discussed further in view of results of laboratory and field studies reported in literature. Experimental studies under controlled conditions are valuable contributions for the mechanistic understanding of effects of drivers on forests, but it is often difficult to upscale their findings to adult trees, grown under real field conditions. In addition, field studies do give insight at field scale, but results are often influenced with other confounding factors and results are only valid for the site where the investigation took place. However, both sources of information give insight in, and allow hypotheses of, the interacting impacts of air quality and climate on forests ecosystems.

Table 1 Tentative overview of key stress factors for different effects in various ecosystem compartments. A ‘+’ signifies that an impact is expected, whereas a ‘−‘ implies no impact. A ‘+/−’ signifies that the impact is likely to be small (after De Vries et al. 2003d)

Nitrogen and sulphur deposition impacts

Overall impacts at increasing nitrogen saturation

Atmospheric deposition has an impact on forest ecosystems through “eutrophication by nitrogen” and soil acidification. Nitrogen is an essential plant nutrient and many terrestrial ecosystems are adapted to conditions of low N availability, a situation that often leads to plant communities with high species diversity (Bobbink et al. 1998). Aber et al. (1989) launched the theory on ecosystem nitrogen saturation and the different stages that can be identified in view of: (i) impacts on soil chemical processes such as mineralization, immobilization, nitrification, affecting N leaching, (ii) plant nutrition and forest growth and (iii) plant species diversity. Until a certain physiological optimum is reached, forests will react to additional N inputs by an increased biomass production. Below the threshold level for growth, however, changes in the ecosystem are already observed. Especially the forest vegetation may gradually change towards more nitrophilic species. Changes in the plant species diversity of the forest undergrowth towards more nitrophilic species have been observed at levels where the growth of forest may still respond to additional N (Ellenberg 1985; Bobbink et al. 1998; Bobbink and Hettelingh 2011). In forested plots with a continuous elevated N input, the ecosystem may approach “N saturation” (Aber et al. 1989). In this stage, the leaching of N, mostly as nitrate (NO3) will increase. At very high NH4 deposition levels, elevated NH4 accumulation may occur in the topsoil (e.g. Roelofs et al. 1985). Elevated leaching of NO3, in combination with sulphate (SO4), is associated with elevated leaching of Al and/or base cations (BC). This cation release, either by weathering or by exchange of protons against Al and BC, implies a decrease in the acid neutralizing capacity of the soil, which is defined as soil acidification (Van Breemen et al. 1984). Soil acidification is associated with a decrease in pH and base cation saturation and an increase of the concentration of Al3+ in the soil solution, (Cronan and Grigal 1995; Marschner 1990; Mengel 1991), which in turn causes a decrease in the ratio of nutrient base cations (Bc), i.e. calcium (Ca), magnesium (Mg) and potassium (K) to Al (e.g. Sverdrup and Warfvinge 1993; Cronan and Grigal 1995). In this stage, a decrease in forest condition and even forest growth is assumed to occur (Aber et al. 1989, 1998).

Forest condition

Enhanced dissolution of Al by deposition of acidifying N and S compounds has been considered a probable threat to forest vitality (Ulrich 1984; Ulrich et al. 1980). Hypothesized mechanisms of Al toxicity are hampered root growth and inhibition of uptake of nutrient base cations (Bc) (Matzner and Murach 1995; Schulze 1989), which may be aggravated by a loss of mycorrhiza or root damage (e.g. Roelofs et al. 1985) or increased levels of dissolved NH4 (Boxman et al. 1988).

An excess input of N may furthermore increase the N content in foliage, which in turn may cause a nutritional imbalance, i.e. causing a shift in deficiency from N to the macro nutrients K, P, Mg and Ca. The chemical composition of the foliage of forest trees is an important indicator for tree nutrition. It provides information on deficiency or excess of nutrients, either in absolute values or relative to the content of other elements. Increased growth rate and elevated N concentrations in foliage may dilute the pool of other nutrients in absolute or relative terms. Furthermore, strong accumulation of N in foliage (e.g. as amino acids) may cause an increased sensitivity to climatic factors, such as frost and drought (Aronsson 1980; De Visser 1994) and diseases and plagues, such as attacks by fungi on (Roelofs et al. 1985; Van Dijk et al. 1992; Flückiger and Braun 1998). It may also cause water stress as a result of increased canopy size, increased shoot/root ratio, and loss of mycorrhizal infection. In this context, a critical N content of 1.8 % in needles has often been mentioned in the literature (Aronsson 1980). Apart from laboratory experiments, the hypothesis of acid deposition induced forest condition decline was based on field observations and foliage analyses, showing e.g. relations between acid deposition, deficiencies of Mg and K based on foliar analysis and yellowing of needles (Zöttl and Mies 1983). In the eighties, several authors considered soil acidification, as the main cause for forest decline (for example Ulrich and Pankrath 1983; Hutchinson et al. 1986)

Forest growth

As summarized by Kozlowski and Pallardy (1996), the requirements for tree growth are carbon dioxide, water and minerals for raw materials, light as energy resource, oxygen, and favourable temperature for growth processes. The capacity for photosynthetic processes (i.e. foliar biomass) and the competition for resources are constraining tree growth. The impact of the deposition of acid N and S compounds on forest growth can thus be positive, due to a fertilizing effect of N deposition, and negative due to a weakening of the forest condition by mechanisms discussed above. Since the productivity of many temperate ecosystems is N limited, adding N via deposition has the potential to increase growth, and therefore to sequester CO2 from the atmosphere. A range of studies has shown positive forest growth and C accumulation responses under low to moderate N additions (Vitousek and Howarth 1991; Aber et al. 1995; Bergh et al. 1999; Franklin et al. 2003). Observations of increased tree growth of European forests (Spiecker et al. 1996) are thus associated with the effect of increased N inputs. Results of 15N experimental studies on the fate of N, combined with C/N ratios in forest ecosystem compartments (e.g. Nadelhoffer et al. 1999) and results of long-term (8–30 year) low dose N fertilizer experiments on the C pool in biomass and soil (e.g. Högberg et al. 2006; Pregitzer et al. 2008; Hyvönen et al. 2008) showed a clear N fertilization effect varying between approximately 20–40 kg C per kg N. Several meta-analysis of impacts of N addition experiments on carbon responses, such net primary production (NPP) and net ecosystem CO2 exchange, indicated similar positive responses (e.g. LeBauer and Treseder 2008; Wamelink et al. 2009; Liu and Greaver 2009).

The nitrogen saturation hypothesis (Aber et al. 1989, 1998) predicts that the final stages of N saturation lead to tree decline and even death. Some studies indeed report decreased growth in response to high N-loads (e.g. Magill et al. 1997, 2004). Furthermore, Boxman et al. (1998) observed a growth improvement in a highly N-saturated Scots pine (Pinus sylvestris L.) stand in which the N input to the forest floor was reduced from ≈ 60 kg N ha−1 yr−1 to <5 kg N ha−1 yr−1 by means of a transparent roof and application of clean, artificial throughfall water. The mechanism for a decline in growth includes the loss of the nutrient base cations Ca, Mg and K from the soil in response to acid deposition, since their (lowered) availability may limit tree growth. Field studies have shown for example that low foliar Ca concentrations, associated with low soil Ca concentrations, were related to elevated rates of dark respiration, and lower rates of net photosynthesis (McLaughlin and Wimmer 1999). Roberts et al. (1989) showed that chlorosis of older needles, associated with reduced photosynthesis and carbohydrate translocation, of Norway spruce (Picea abies) and silver fir (Abies alba) in mountainous areas of central Europe resulted directly from foliar Mg deficiency. These changes can result in lower rates of root growth and wood development. They indicated that foliar Mg deficiency resulted primarily from low rates of Mg uptake due to the low availability of Mg in some acidified soils in central Europe.

Species diversity of ground vegetation

An important aspect of forest biodiversity is the species composition of the ground vegetation (Fischer et al. 2009). In many countries a shift in species composition has been reported, usually from a cryptogam-dominated to a grass-dominated undergrowth, and an increase in species indicative for nitrogen-rich circumstances (Van Dobben et al. 1999). These changes are now mostly ascribed to deposition of nitrogen compounds originating from combustion processes and agriculture (e.g. Tyler 1987; Bobbink et al. 1998; Thimonier et al. 1992; Gilliam 2006). In response, many experimental nitrogen addition studies have been carried out to establish empirical ‘critical loads’ for nitrogen, i.e. the maximum deposition tolerated by a given ecosystem without a change in species composition (e.g. Achermann and Bobbink 2003) and models have been developed to assess critical loads by simulation (Van Dobben et al. 2006; De Vries et al. 2010). Generalisation of the impacts of N deposition on the plant species diversity of forest ecosystems across Europe is, however, difficult because the effects of air-borne N deposition depends on (i) the duration and total amount of the inputs, (ii) the form of the N input, (iii) the intrinsic sensitivity of the present (plant) species and (iv) the abiotic conditions in the ecosystem (Achermann and Bobbink 2003). A way to gain insight into N deposition impacts is to relate monitoring data of the ground vegetation with environmental factors, including N deposition, as described further.

Ozone effects

As a strong oxidant, O3 causes several types of visible symptoms including chlorosis and necrosis. Compounds resulting from oxidation by O3 interfere with the cell’s energy production, decrease the photosynthesis, resulting in slower plant growth and they impair water use efficiency and other functions. Plants weakened by O3 may also be more susceptible to pests, disease, and drought. Ambient O3 concentrations in rural and forested areas have been shown to reach concentrations high enough to produce phytotoxic effects in native vegetation (e.g. Bussotti and Ferretti 1998; Sanz and Millán 1998; Skelly et al. 1999; Sanz et al. 2000; Innes et al. 2001; Ashmore 2005; Ferretti et al. 2007a). For example, from 2000 to 2002, a total of 65 species representing 52 genuses were observed to manifest O3 symptoms in Spain, Italy and Switzerland. Symptomatic species were found at 47 % of the common monitoring plots in 2001 and 38 % in 2002 (Ferretti et al. 2007a).

Forest condition

Several authors have hypothesized that ambient O3 is one of the key factors resulting in spatial and temporal changes of tree crown defoliation, interacting with acid deposition (Reich 1987; Guderian 1985; Sandermann 1996; Takemoto et al. 2001). It is assumed that their combined effects differ significantly from the sum of individual effects due to their complex synergistic or antagonistic interactions (Bytnerowicz et al. 2007). Relatively few studies have, however, investigated O3 effects on crown condition. Datasets describing O3 effects on adult trees from long-term studies are scarce, given the fact that most of the available knowledge about O3 effects on trees has been gained from juvenile trees (e.g. Baumgarten et al. 2000; Bortier et al. 2000; Kolb and Matyssek 2003; Matyssek and Sandermann 2003; Nunn et al. 2005; Schaub et al. 2005; Gerosa et al. 2009; Vapaavuori et al. 2009) or from intensive studies which focused on a better mechanistic understanding of the physiological and bio-chemical processes within the canopy of several trees (e.g. Matyssek et al. 2007a; Wieser et al. 2013).

The complex interactions between environmental conditions as they occur inside the canopy of adult trees and O3 are thus not fully clarified and the relationship between forest tree crown condition (to describe the condition of entire forest ecosystems) and O3 effects are still not well established (Manning 2005; Ollinger et al. 1997; Percy and Ferretti 2004; Paoletti 2006). Response validation under actual forest conditions with adult trees has become feasible upon recent developments of free-air O3 fumigation systems (Musselman and Hale 1997; Karnosky et al. 2007b; Werner and Fabian 2002). Also, due to the expected increase in future ambient O3 concentrations (Fowler et al. 1999; Percy et al. 2003) it is necessary to determine if exposure to ambient O3 levels actually affects crown condition in natural forests (Manning 2005).

Until approximately 2000, most O3 research focused impacts of the external O3 exposure on forest ecosystems, using either an average or peak O3 concentrations or an accumulated dose over a threshold O3 concentration of 40 ppb (AOT40). The AOT40 index was developed within the context of the UN-ECE as a metric to derive a ‘critical level’, defined as ‘the concentration in the atmosphere above which direct adverse effects on receptors such as plants and ecosystems may occur, according to present knowledge’ (Ashmore and Fuhrer 2000). The currently set AOT 40 threshold of 5 ppmh to protect forests from ozone (De Santis 1999) is frequently exceeded in many parts of Europe (Karlsson et al. 2004b). From 2000 to 2002, this critical level was exceeded at 77–100 % of the 81 monitored forest sites across France, Italy, Spain and Switzerland (Gerosa et al. 2007). However, already for more than a decade, the uncertainties associated with the metric and the currently applied critical level of 5 ppm h are clearly articulated under the LRTAP Convention and individual countries are encouraged to modify the approach for application to their local conditions (Ashmore and Fuhrer 2000). Today, there is consensus that the so-called Phytotoxic Ozone Dose (i.e. the accumulated stomatal flux of ozone above a flux threshold of Y, PODY) exposure index provides much stronger relationships with O3 effects compared to external concentration based exposure indices (Emberson et al. 2000; Simpson et al. 2003; Karlsson et al. 2003; Ashmore et al. 2004; Matyssek et al. 2007a; Mills et al. 2010; Grünhage et al. 2012). Karlsson et al. (2007), for example, has demonstrated that the O3 impact on biomass, as well as on visible injury could be better explained by the accumulated stomatal O3 uptake flux as compared to the concentration based AOT40 index, for young trees. For both, young and adult beech (Fagus sylvatica L.) trees, Baumgarten et al. (2000) showed an improved association between injury and cumulative O3 uptake, as compared to the exposure-based approach. Within the UN-ECE, the PODY is therefore now suggested for regional ozone risk assessment instead of AOT40 (Mills et al. 2010).

The Phytotoxic Ozone Dose is calculated as the sum over time of the differences between hourly mean values of the stomatal flux of O3 based on the projected leaf area (PLA), Fst (in nmol m−2 PLA s−1) and a stomatal flux threshold of Y nmol m−2 s−1, for an appropriate time-window for the periods when Fst exceeds The value of PODY is formally named AFstY; Accumulated stomatal Flux over thresholds of Y nmol m−2 s−1). The POD approach takes into account the varying influences of air temperature, water vapour pressure deficit (VPD) of the surrounding leaves, light (irradiance), soil water potential (SWP) or plant available water (PAW), O3 concentration and plant development (phenology) on the stomatal flux of O3. It therefore provides an estimate of the critical amount of O3 entering through the stomata and reaching the sites of action inside the plant and is species-specific. The stomatal flux-based critical level of O3, CLef mmol m−2 PLA, is then the cumulative stomatal flux of ozone, PODY, above which direct adverse effects may occur according to present knowledge (Mills et al. 2010). Simpson et al. (2007) applied the EMEP chemical transport model to map both indicators of O3 damage across Europe metrics, the AOT40 and the AFstY for two vegetation types, wheat (Triticum spp.) and beech forests. The results show that exceedances of critical levels for either type of indicator are widespread, but that the indicators give very different spatial patterns of hot spots across Europe (Fig. 2).

Fig. 2
figure 2

Calculated geographic variation over Europe of the ozone exposure index AOT40f [ppb h], calculated as the accumulated exposure over a threshold of 40 ppb, (left), and the ozone exposure index AFstY [mmol m−2], calculated as the stomatal flux Fst over a threshold Y of 1.6 nmol m−2 s−1 (right). Source: Simpson et al. 2007

The decoupling between external O3 exposure and internal O3 flux becomes particular apparent under limited light conditions or drought (Matyssek et al. 2004; Nunn et al. 2006). Schaub et al. (2007) found a divergence between AOT40 and O3 flux indices under conditions where water vapour deficit (VPD) and soil moisture deficit may limit physiological activity, stomatal conductance and hence O3 uptake into the leaves. The advantage of O3 risk assessment by means of PODy can be demonstrated for Fennoscandia, where the potentially most important effect of a future climate change on stomatal O3 uptake would be mediated through the earlier onset of the start of the growing season, therefore leading to a substantial and increasing risk for negative impacts on vegetation in Northern Sweden (Karlsson et al. 2007). Currently, the flux-based O3 uptake approach has been applied for adult trees of beech (Matyssek et al. 2004; Baumgarten et al. 2009; Schaub et al. 2007; Braun et al. 2010; Gerosa et al. 2008), conifer species such as spruce, larch (Larix) and stone pine (Pinus pinea) (Wieser et al. 2013; Wieser et al. 2002; Nunn et al. 2006; Zapletal et al. 2012) and Mediterranean species such as quercus species (Calatayud et al. 2011; Gerosa et al. 2008; Manes et al. 2007).

Forest growth

Increasing O3 levels have the potential of limiting carbon sink strength of forest ecosystems (IPCC 2007; Pretzsch et al. 2010; Sitch et al. 2007; Wittig et al. 2009), due to reduced net photosynthesis, accelerated leaf senescence and increased dark respiration (Matyssek et al. 2010b).

The negative impacts of O3 on physiological functions and biomass production have been investigated in a number of experiments. Most has focused on physiological or biochemical effects of O3 exposure, in attempts to determine mechanisms of toxicity, dose response work for risk assessment, or development of diagnostic tests for O3 injury. Many of these demonstrated relationships between O3 exposure and reductions in both growth and leaf gas exchange (Matyssek et al. 1995; Sandermann et al. 1997; Novak et al. 2005). Experimental studies, however, also show high variability of different growth parameters in response to elevated O3 exposures. An example is Kranzberg Forest, Germany, where 50- to 60-year-old mixed stand of Norway spruce and European beech tree cohorts have been exposed to double ambient O3 from 2000 through 2007 (Werner and Fabian 2002). Within this stand, Pretzsch et al. (2010) found that O3 fumigation induced a shift in resource allocation into height growth at the expense of diameter growth. Norway spruce showed a non-significant loss of stem diameter and basal area increment at breast height of about 12 %. Similar reductions were reported for beech with −12 % for stem increment and −17 % for basal area. The height growth of spruce, however, significantly increased by 70 %, whereas this was not the case for beech. In terms of total volume growth, twice ambient O3 concentrations thus induced a 2.7 % growth increase for spruce and a 43.5 % growth reduction for beech (Pretzsch et al. 2010).

Many experiments have been carried out to assess how much exposure to O3 actually affects tree growth (Sandermann et al. 1997; Chappelka and Samuelson 1998; Dizengremel 2001; Samuelson and Kelly 2001; Ferretti et al. 2002; Karlsson et al. 2003; Karlsson et al. 2004a). Despite the high variability of different tree species, the overall impact of experimental results shows a clear decline in growth in response to elevated O3 exposures. With a meta-analysis, including 313 articles covering the period 1970 through 2006, Wittig et al. (2009) addressed the question as how current ambient O3 concentrations decrease productivity relative to preindustrial O3 concentrations. They found that current ambient O3 concentrations of 40 ppb on average significantly reduced total biomass (and therefore carbon storage) of young trees by 7 % compared to trees grown under preindustrial O3 levels. Elevated O3 concentrations of 64 ppb reduced total biomass by 11 % relative to trees grown under ambient O3 concentrations.

While results of experimental studies under controlled conditions, as described above, are valuable contributions for the mechanistic understanding of O3 effects on forest trees, it is difficult to upscale their findings to adult trees, grown under real field conditions (Kolb et al. 1997; Baumgarten et al. 2000; Schaub et al. 2005). For example Nunn et al. (2005) applied a correlative O3 uptake model (Emberson et al. 2000) to relate O3 flux to radial growth and showed that despite similar O3 regimes, young beech seedling grown in phytotrons were more sensitive to O3 relative to adult beech trees due to higher O3 uptake. On the other hand, Kitao et al. (2012) found that O3-induced changes in leaf-level photosynthetic parameters had the capacity of limiting the seasonal photosynthetic carbon gain in the same adult beech trees, but the reduction in leaf photosynthesis could not immediately be related to the extent of growth reduction in stems.

Much has been written about the presumed or probable effects of O3 on growth of trees in the forest (Manning 2005). In Europe, estimates of average forest yield losses around 10 % were suggested (Broadmeadow 1998), while values as high as 42 % were suggested for beech (Pretzsch et al. 2010). At the same time, in the eastern USA, growth losses of 0–10 % were predicted (Chappelka and Samuelson 1998). Growth reductions in eastern hardwood trees, ranging from 4 to 12 %, were predicted for north-eastern forests (Ollinger et al. 1997). Verifying these loss predictions is difficult, as the bases for most of them are laboratory experiments that may not be relevant for actual field conditions. In fact, much of the available information about O3 effects on trees does not apply well to actual forest conditions and large trees (Matyssek and Sandermann 2003). Due to the complex interrelationships among the parameters affecting tree growth over the life span of a tree, establishing a cause and effect relationship for ambient O3 exposure and tree growth in forests has proven to be a difficult task. Manning (2005) suggests, that smaller, better-defined studies in the forest should lead to a better understanding of factor interaction with O3 in affecting tree growth. While it is difficult to work in natural ambient conditions, studies under these conditions are essential to develop a better understanding of how O3 and environmental factors interact to affect radial growth in trees. In this context, Matyssek et al. (2010b) promotes free-air fumigation studies, which allow whole-plant assessment of adult trees to gain a better understanding of the modifying impact of biotic agents in the presence of natural mycorrhization.

Species diversity of ground vegetation

Novak et al. (2008) demonstrated that plant competition is an important factor affecting plant responses to O3 stress, but the direction and severity of these effects depend on the interacting species and micro site conditions. Therefore, competition in combination with elevated O3 exposures may alter species composition and diversity over the long run, but quantitative information on this potential impact is yet lacking.

Climate impacts

Crown condition

Crown transparency, a measure of reduced foliage mass has been shown to be highly influenced by summer drought (Dobbertin and Rigling 2006; Seidling 2007; Carnicer et al. 2011). Drought may affect crown condition and mortality in many ways, i.e. by increasing needle or leaf loss during or after the drought (Christiansen 1992; Gruber 1990; Bussotti et al. 2005; Heiniger and Schmid 1986, 1989; Livsey and Barklund 1993), dying of buds and branches (Aronsson et al. 1978), root mortality and predisposition to attack by Armillaria spp. or other root pathogens (Whitney and Timmer 1983) and reduced foliar nutrient contents (Wehrmann 1961). Experimentally introduced drought stress to stands of spruce and pine have shown a considerable defoliation, measured as a reduction in foliage mass (Christiansen and Nilsen 1990), and yellowing caused by reduced foliar nutrient concentrations (Dambrine et al. 1993). In addition, a dry or hot summer reduces growth (Spiecker 1986), increases attacks by bark beetles and other insects (Worrell 1983), and induces heavy flowering (Owens and Blake 1985), partly at the expense of vegetative shoot formation (Lüscher 1989).

Spring or early summer drought have the potential to reduce current year leaf or needle size (Dobbertin et al. 2010). Extreme drought may either lead to premature leaf or needle shedding causing high crown transparency during the summer, or for conifers it may cause needle loss during autumn which will up show as reduced foliage mass in the following year (Solberg 2004). Most commonly for temperate trees with only one annual flush, however, is that late summer drought causes a decrease in primordial, which subsequently results in reduced foliage mass and shorter shoots in the following season. Water stress is therefore very important with respect to the monitored forest condition. Innes (1993) mentioned that the most alarming and frequent observations of a decrease in forest condition in Central Europe coincided with the dry years 1982 and 1983. Landmann (1995) mentioned that defoliation appears to be highest in soils poorly supplied with water and/or in stands in which trees, at some stage of development, have suffered from competition for water.

Mortality of forest and woodland plants has also been increasingly observed as a response to extreme drought and heat events (Allen and Breshears 1998; Allen et al. 2010; Hanson and Weltzin 2000; Breshears et al. 2005; Shaw et al. 2005; Berg et al. 2006; Gitlin et al. 2006; Carnicer et al. 2011). Drought and heat events, such as in 2003 in Europe, reduced net primary productivity (Ciais et al. 2005; Bréda et al. 2006) and increased mortality (Martínez-Vilalta and Piñol 2002; Lloret et al. 2005; Bigler et al. 2006; Dobbertin and Rigling 2006). A widespread increase in tree mortality in the Western United States has also recently been observed and attributed mainly to warming and consequent increases in water stress (Van Mantgem et al. 2009).

Forest growth

Change in water supply can severely alter above- and below-ground tree growth, as well as affect the vitality of trees (Dobbertin 2005; Waring 1987). Increasing temperature under equal or reduced precipitation will increase the evaporative pressure on plants (Rebetez and Dobbertin 2004). Plants react to drought stress with stomata closure to reduce water loss or leaf wilting (Zweifel et al. 2007). While the reduction of stomatal conductance is an adaptation to drought stress, closed stomata reduces the uptake of atmospheric carbon dioxide, and in turn limits the carbon dioxide fertilization of photosynthesis and therefore reduces carbon uptake (Körner 2003). Reduced photosynthesis may not only result in reduced tree growth and related carbon allocation (Zweifel et al. 2006; Adams et al. 2009; McDowell et al. 2008; Zweifel et al. 2009), but it may also reduce the resistance of trees against insects and pathogens (Kolb et al. 1998). As different species are differently adapted to drought, the inter-tree competition will also be altered.

It is well known that reduced water availability results in decreased tree growth and reduced foliage mass, either by premature leaf shedding or reduced foliage production in the following year (Pallardy 2008). The extreme drought of 2003, for example, led to a reduction of forest growth in low-land forests in central Europe of up to 60 % (Dobbertin 2005; Jolly et al. 2005; Leuzinger et al. 2005). Although respiration loss was also reduced, many forests released more carbon in 2003 than they sequestrated via tree stem growth (Ciais et al. 2005) resulting in a net loss of carbon. Many tree-ring studies have successfully linked ring width reduction to reduced precipitation or increased temperature of various months during or before the current the tree growing season (Schweingruber 1996).

Impacts as derived from monitoring studies

Geographic variation and trends in acidic atmospheric deposition and ozone exposure at ICP Forests Level II plots

Atmospheric deposition of nitrogen and sulphur

In the framework of the activities carried out on the ICP Forests Level II plots, atmospheric deposition is measured on the basis of a combination of (i) bulk deposition in the open field in the vicinity (within 1 kilometre) of the Intensive Monitoring plot and (ii) throughfall below the forest canopy within the plot. For monitoring of throughfall data, use was mostly made of funnels whereas gutters were used at a minority of the plots. The sampling frequencies for bulk deposition, throughfall and stemflow measurements were mostly weekly or fortnightly and the number of bulk deposition and throughfall samplers used mostly ranged from 2–6 and 10–15, respectively (De Vries et al. 2003c).

The geographic variation in acidic components over Europe is illustrated by maps of the average annual S and N deposition in the period 2005–2007, as presented in Fig. 3. The results show that high S deposition, and to a lesser extent NO3-N deposition, mainly occur in central Europe (Poland, Germany, Czech Republic, Slovak Republic, Slovenia and Romania). High nitrate (N-NO3) and ammonium (N-NH4) deposition also occurs in the Netherlands, Belgium, parts of France, and north of Italy and Spain. The lowest nitrogen deposition was observed on plots in Scandinavia with values mostly below 2 kg N ha−1 yr−1 for both N-NO3 and N-NH4.

Fig. 3
figure 3

Mean annual average sulphur (SO4-S; left), nitrate-nitrogen (NO3-N; middle) and ammonium-nitrogen in (NH4-N; right) throughfall deposition in 2005 to 2007. Source: Fischer et al. (2010)

Trends in S and N deposition at Intensive Monitoring plots, assessed over the last 10–15 years, show a clear decline for S and a less significant decline for N. This follows from recent studies by, Fischer et al. (2010) and Lorenz and Becher (2012), who presented results of trends in S and N deposition between 1998–2007 and 2001–2010, respectively, based on 87 plots with throughfall measurements and 55 plots with bulk deposition measurements. Their results showed an overall strong decrease of annual S deposition both in bulk deposition and throughfall. For example, mean S throughfall deposition decreased from 10 kg S ha−1 yr−1 in 1998 to 5 kg S ha−1 yr−1 in 2010, implying an average annual decline of 0.42 kg S ha−1 yr−1 over 12 years. The mean S bulk deposition decreased from about 7 kg S ha−1 yr−1 in 1998 to 4 kg S ha−1 yr−1 in 2010, implying an average annual decline of 0.25 kg S ha−1 yr−1. The mean throughfall deposition of inorganic N on plots with continuous measurements from 1998 to 2010 decreased from approximately 11 to 9 kg N ha−1 yr−1, i.e. circa 0.16 kg N ha−1 yr−1. Trends were clearer for NO3 than for NH4. At 83 % of the intensive monitoring plots, no significant decrease in atmospheric throughfall deposition was detected (Fischer et al. 2010), illustrating the limited reduction in ammonia (NH3) emissions since 1998.

Maps of trends in S and N deposition for the period 1998–2007 (Fig. 4) show that decreases in S deposition, and to a lesser extent NO3-N deposition, mainly occur in central Europe, although there is no consistent pattern in the results due to many local factors affecting the deposition. The deposition of NH4-N, being mainly due to agricultural activities, mostly decrease in high deposition areas with intensive agriculture in the Netherlands, Belgium and Germany in response to emission reduction measures taken in those countries.

Fig. 4
figure 4

Trends in sulphate-sulphur (SO4-S; left), nitrate-nitrogen (NO3-N; middle) and ammonium-nitrogen in (NH4-N; right) throughfall deposition from 1998 to 2007. Source: Fischer et al. (2010)

Tropospheric ozone concentrations

In National Air Quality Networks, O3 concentrations are measured with continuous monitors that are costly and mostly concentrated in urban or suburban areas, making it very difficult to infer the O3 spatial patterns, especially in complex mountainous areas or along the sea shore (Bacci et al. 1990; Millán et al. 2000). Passive samplers are considered an alternative, or a complement to continuous monitors, e.g. when a higher density in the measuring network is required. Their use in assessing the ecological effects of tropospheric O3 has significantly increased in recent years (e.g. Manning et al. 1996; Blum et al. 1997; Krupa and Legge 2000; Sanz et al. 2001). In the framework of the activities carried out on the ICP Forests Level II plots, mainly O3 passive sampler measurements were thus carried out since 2000 to (i) gain knowledge into the spatial and temporal exposure trends of O3 concentrations and (ii) to evaluate the potential risk to forest ecosystems (Sanz and Krause 2001).

Biweekly O3 measurements at the ICP Forests Level II network in 2002 showed concentrations ranging between 30 and 45 ppb at most of the stations. Patterns in Germany and France show increasing concentrations going from North to South and towards high-elevation areas, with the seasonal average concentrations rarely exceeding 45 ppb, but there was no distinct geographical pattern in the UK. In Switzerland and Italy, about 30 % of the plots ranged between 45 and 60 ppb in 2002 (Lorenz et al. 2005). Passive sampler measurements for the period 2000–2004 show that the highest O3 concentrations in forest areas across Europe occur in the south (Fig. 5), in particular Italy, Switzerland and Spain (Sanz et al. 2007). Highest measured values occurred in central and northern Italy and in southern Switzerland. These results are consistent with the output from continuous monitors (De Leeuw and Bogman 2001; Fiala et al. 2002, 2003). The hot-spots in Italy and Switzerland are under the influence of the industry and traffic emissions originating from the Po Valley, including the city of Milan (e.g. Novak et al. 2003; Staffelbach et al. 1997b; Staffelbach et al. 1997a). Spring and summer O3 concentrations are high enough to be of potential risk to sensitive plants in nearly all regions of Europe. The currently recognized AOT40 thresholds for damage to forest trees are frequently exceeded in many areas of Europe.

Fig. 5
figure 5

Mean ozone concentrations for April–September in 2000–2004 at 91 Intensive Monitoring plots. Due to the higher solar radiation, concentrations typically increase from North to South Europe. Source: EC-UN/ECE 2007

The annual changes in O3 concentrations in the 2000–2004 period observed at the monitoring plots (Fig. 6) were mainly related to different climatic conditions (Sanz et al. 2007), since O3 concentrations increase with higher air temperature and solar radiation and decrease in cloudy and rainy periods (e.g. Amoriello et al. 2003). Warm and sunny weather enhances the O3 concentration because of the emission of volatile hydrocarbons (including vegetation emissions) increases with temperature, higher solar radiation increases photochemical processes, and high temperature results in more rapid chemical O3 formation (Fiala et al. 2003). Over the 2000–2004 period, highest concentrations were measured in 2003, a year with one of the hottest summers on record for Europe, producing the highest O3 values for Central Europe since the 80s. The high air temperatures measured during the summer of 2001 are also reflected in the higher O3 levels during that year. In the five-year study period, concentrations were lowest in 2004, due to the relatively low solar radiation recorded during April–September 2004 (Fig. 6).

Fig. 6
figure 6

Shares of plots belonging to 8 classes of ozone concentrations based on the mean ozone values (April–September) for the period 2000–2004. Only plots with at least 70 % data cover during the observation period were considered. Source: (EC-UN/ECE 2007)

Water and element budgets and nutrient-acidity status

Water budgets

Water influences the availability of nutrients by affecting bio/geochemical processes and the loss of nutrients from the rooting zone by leaching. To understand the element cycles in a forest ecosystem, a quantification of the hydrological situation is indispensable. The hydrological budget is mainly determined by the input of water by rainfall and the loss of water from the forest by interception evaporation, soil evaporation, transpiration and leaching. The first and only water flux study at a European wide scale has been carried out by Van der Salm et al. (2007b), using data from hundreds of intensive forest monitoring plots for the period 1996-1998. They calculated interception evaporation on the basis of measured bulk deposition and throughfall, while soil evaporation, transpiration and leaching were calculated with a hydrological simulation model called SWATRE. Van der Salm et al. (2007b) concluded that yearly but also monthly and two weekly measured throughfall data were quite well simulated for 245 plots passing a calibration procedure and various quality checks.

Evapotranspiration and leaching fluxes varied considerably for the investigated plots (Fig. 7). Transpiration fluxes are generally highest (>450 mm) in central Europe where both rainfall and radiation are relatively high (Fig. 7 left). At most plots in southern Europe, transpiration is somewhat lower due to the limited rainfall, except for plots in mountainous areas. In northern Europe transpiration fluxes are lower than in central Europe due to the decrease in radiation with increasing latitude and the relatively low precipitation in part of Sweden and Finland. The leaching fluxes (Fig. 7 right) closely resemble the precipitation and throughfall maps. High leaching fluxes are found in North-western Scandinavia, Britain and Ireland and on some plots in mountainous areas in central and southern Europe.

Fig. 7
figure 7

Average annual transpiration fluxes (left) and average leaching fluxes (right) for 245 European forest monitoring plots for the period 1996–1998. Source: Van der Salm et al. (2007b)

Element budgets

At 121 intensive forest monitoring plots, the inputs and outputs of S and N compounds (total N, NH4 and NO3), Al and BC (Ca, Mg, K and Na) have been assessed for the same period 1996–1998, as presented in De Vries et al. (2007). Total N and S deposition was derived from both bulk deposition and throughfall while accounting for forest canopy exchange, while N and S leaching was based on a combination of measured dissolved element concentrations and the water balance model SWAP, using methods described in detail in De Vries et al. (2007). Results showed that on average SO4 behaves as a tracer, since average SO4 leaching fluxes are comparable to SO4 deposition, whereas N is strongly retained in the soil and/or denitrified since N leaching is generally much lower than N deposition (Fig. 8a, b). The higher leaching of SO4 implies that S deposition is still the dominant source of actual soil acidification, despite the fact that N deposition is higher on average.

Fig. 8
figure 8

Relations between total deposition and leaching fluxes of S (a) and N (b) at 121 forest monitoring sites for the period 1996–1998. Source: De Vries et al. (2007)

Although average S leaching is close to S deposition, there is a large variation in both fluxes that can partly be attributed to errors in both the input and output assessment (Fig. 8a). At several sites the leaching of SO4 is considerably higher than the deposition indicating a strong release of SO4, which has been adsorbed during previous decades. Results of the leaching of total N and NO3 against the total N deposition showed that the leaching of N is generally negligible below a total N input of 10 kg ha−1 yr−1 (Fig. 8b). These results are in accordance with those found by e.g. Gundersen et al. (1998), Dise et al. (1998a, b, 2009), MacDonald et al. (2002) and Van der Salm et al. (2007a), who found those results, largely based on a review of literature data while using NO3 leaching and throughfall N inputs. At N inputs between 10 and 20 kg ha−1 yr−1, leaching of N is generally elevated, although lower than the input indicating N retention at the plots. At N inputs above 20 kg ha−1 yr−1, N leaching is also mostly elevated and in several cases (seven plots), it is near or even above (for two plots) the N deposition (Fig. 8b). The latter situation indicates a clear disturbance in the N cycle in response to the elevated N input. NO3 dominates N leaching and the increase in N leaching in response to elevated N loads (above 20 kg ha−1 yr−1) appears to be larger for deciduous trees than for conifers (De Vries et al. 2007).

In acidic soils, atmospheric deposition of S and N compounds do lead to elevated Al concentrations, in response to elevated concentrations of SO4 and NO3. The range in Al leaching fluxes at the same 121 Intensive Monitoring plots is quite comparable to the range in S and N deposition (Fig. 9a). At several sites, Al leaching fluxes are even higher than the acid deposition. This occurred at sites where the leaching of SO4 is higher than the input (S release from the soil causing acidification) and acid input is almost completely buffered by the release of Al. Sites where the Al leaching flux is lower imply that part of the potential acid input is buffered by S and/or N retention and/or BC release. In soils with a pH above 5.0, the release of Al is generally negligible, independent of the S and N input, since BC release by weathering and cation exchange buffers the incoming net acidity in those soils, but at sites below pH 4.5 there is high correlation between Al concentration and the sum of SO4 and NO3 concentration (Fig. 9b). Comparable results were found by Dise et al. (2001).

Fig. 9
figure 9

Relations between Al leaching fluxes and the total deposition flux of S+N (a) and the concentration of total Al against total SO4+NO3 in the subsoil of Intensive Monitoring plots with a pH < 4.5 at 121forest monitoring sites (b). The solid line in A represents a regression line being equal to: Al = −95 + 0.74 (SO4+NO3) (R 2 = 0.86). Source: De Vries et al. (2007)

Soil solution and foliar chemistry

Soil chemistry data at a single point in time do not give adequate information on the impacts of acid atmospheric deposition data due to natural soil acidification that may also lead to low pH values and low base saturation under natural circumstances. De Vos and Cools (2011) summarize the results of a second European Forest Soil Condition survey in which they evaluate carbon and nutrient stocks, acidification and contamination levels, as compared to the first European forest soil condition survey in the 1990s on the ICP Forests Level I network (Vanmechelen et al. 1997). About 44 % of the number of Level I plots assessed during the first forest soil condition survey were revisited during the BioSoil survey (2326 plots), but only a limited set of 127 permanent Level II plots were investigated, with about 75 % of them having historical soil data thus enabling assessment of temporal changes. Unfortunately, since methods were not always consistent, De Vos and Cools (2011) did not present any changes in soil chemistry at European scale. Unlike soil chemistry, one time monitoring data of soil solution chemistry and foliar chemistry already shows effects of atmospheric N and S deposition by elevated N and S concentrations in soil solution and forest foliage. Below the results are discussed of both one-time surveys and trends at European scale.

Soil solution chemistry

At the above mentioned 121 Intensive Monitoring plots, De Vries et al. (2007) showed that in acidic sub-soils with a pH below 4.5, more than 80 % of the variation in Al concentration could be explained by a variation in SO4 and NO3 concentrations (Fig. 9b), which in turn were strongly related to the deposition of S and N, respectively (see Fig. 8). In addition, in coastal areas peak Al values in the range 200–800 umol l−1 are mostly caused by sea salts blown inland by storms (Lange et al. 2006). Although SO4 is important in releasing Al, results showed that NO3 concentrations were mostly higher, reflecting the increasing role of N in soil acidification. Results of the Al concentrations and of Al/Bc ratios (Bc is Ca + Mg + K) in the 121 Intensive Monitoring plots in the year 2000 in different layers are given in Fig. 10. Results show a clear increase in Al concentration and Al/Bc ratio going downward from the organic layer to the mineral soil (Fig. 5a, b) due to Bc uptake and elevated Al release with depth. Differences between the mineral topsoil (0–20 cm) and the mineral subsoil (20–80 cm) are smaller than the difference between the organic layer and mineral topsoil. Although different numbers of samplers at different depths influence the results, it does indicate the large impact of Al release and nitrification in the mineral topsoil. More details on the results are given in De Vries et al. (2003b).

Fig. 10
figure 10

Cumulative frequency distributions of the concentrations of Al (a) and of the ratios of Al against the base cations Ca, Mg and K (b) in the organic layer, the mineral topsoil (0–20 cm) and the mineral subsoil (20–80 cm) of 121 Intensive Monitoring plots in the period 1996–1998. Source: De Vries et al. (2003b)

Recently, Iost et al. (2012) examined trends in two soil solution parameters, i.e. the base cation/aluminium (Bc/Al) ratio and the dissolved inorganic N concentration, in terms of exceedances of possible critical limits for these parameters mentioned in literature (Sverdrup and Warfvinge 1993). Data were evaluated for the organic layer, mineral topsoil (0–40 cm soil depth) and mineral subsoil (below 40 cm soil depth) from the early 1990s to 2006. For the time trend analysis, they selected all plots and soil depths that provided data for at least five consecutive years. Their results showed no trends in critical limit exceedance for the Bc/Al ratio in nearly all plots. For the N concentration no trends were detected in most of the plots, but in cases where a trend was detected, it was usually decreasing.

The constant Bc/Al ratio may seem strange in view of the reduced S deposition. The most likely explanation is that this reduction has led to both a reduced Al and Bc concentration because of less buffering due to the decrease in potential acid deposition. Reported monitoring data for the Lehstenbach site in Germany experienced a significant decrease in acid deposition load between 1988 and 1999 and as a consequence, a significant decrease in soil solution Al concentration. At the same time, however, soil solution Ca and Mg decreased even stronger resulting in a decline in soil solution Ca to Al ratio (Alewell et al. 2000). Leeters et al. (2007) also found a decrease of both Al, K and Mg concentration in response to a decrease in S and N deposition for 14 intensive monitoring plots In the Netherlands, in the period 1990–2001.

Foliar nutrient status

Insight in the possible impact of N deposition on a nutritional imbalance in foliage has been derived for 109 Intensive Monitoring plots with information on both the chemical composition of the foliage and total N deposition. Ranges in N deposition at plots with a balanced and unbalanced ratio of the base cations K, Ca or Mg to N are given in Table 2 (De Vries et al. 2003d). Criteria for deficiency (low) or excess (high) in nutrient contents and nutrient ratios are based on Stefan et al. (1997). The number of plots with a clearly unbalanced nutrition was approximately 50 %, with an unbalanced P/N ratio being most common, followed by Mg/N. When considering all elements, including P K, Ca and Mg, there was an unbalanced ratio of one or more of those elements at 57 of the 109 plots. The results indicated a larger N deposition at the plots with an unbalanced ratio. Considering all base cations, the median N deposition was 10 kg ha−1 yr−1 (range of 2–32 kg ha−1 yr−1) at the plots with a balanced nutrition and 21 kg ha−1 yr−1 (range of 8–34 kg ha−1 yr−1) at the plots with an unbalanced nutrition. Results obtained on the relation between N contents in first year needles of Scots pine and total N deposition at 68 Intensive Monitoring plots in Europe (Fig. 11) indicate that N contents in foliage may exceed the previously mentioned critical N content of 1.8 % related to drought and frost stress at N deposition levels above 20 kg ha−1 yr−1.

Table 2 Ranges in N deposition at 109 Intensive Monitoring plots in Europe with a balanced and unbalanced ratio of P, K, Ca or Mg to N
Fig. 11
figure 11

Relationship between N contents in first year needles of Scots pine and total N deposition at 68 plots in Europe. Source: De Vries et al. (2003d)

Regarding foliar chemistry, no trend assessments have yet been made at European scale, but results are available at country level. Leeters et al. (2007), for example, evaluated trends in foliar N and S contents in the foliage of 9 intensive monitoring plots In the Netherlands in the period 1990–2005. Their analysis included three plots for Douglas-fir (Pseudotsuga menziesii), Pine and Oak (Quercus), respectively. Results showed that N and S contents in the foliage decreases significantly over this period 1990–2005 for the coniferous tree species (Douglas-fir and Pine) but no trends were found for Oak.

Forest crown condition

A possible method to explore relationships between air pollution and crown condition is by means of correlative studies, i.e. statistical analysis of data that have been collected within the framework of monitoring programmes, even without an ad-hoc experimental design (Ferretti et al. 2007a, b). This is the case of many studies carried out to evaluate relationships between air pollution and tree crown condition (Innes and Whittaker 1993; Mather 1994; Dobbertin et al. 1997; Klap et al. 2000; Zierl 2002; Ferretti et al. 2003; Ozolincius et al. 2005). These studies have a number of inherent limitations (Ferretti et al. 2003; Karlsson et al. 2003; Manning 2005): (i) high variability of the forests being examined (in experimental terms: non-homogeneous experimental material), (ii) unbalanced and virtually biased allocation of plots to air pollution levels (non-random allocation of treatments, absence of replication), (iii) effects of confounding factors and covariates, (iv) absence of a true control, (v) potential inadequacy of sampling numbers for proper statistical analysis and interpretation, (vi) co-linearity between multi data, (vii) constraints in data availability and (viii) nature of the available response indicators. However, while these limitations must be taken into account at the interpretation stage, in many cases correlative studies are one of the few available options to gain an insight into the possible relationships between O3 exposure and forest conditions across wide geographical areas. In this section, we report results based on monitoring studies using a correlative approach.

European scale studies

Acid deposition

Since 1986, the crown condition variables defoliation and discolouration have been monitored by assessments on permanent plots throughout Europe, as indicators of forest health. This effort has been carried out on plots at systematically sampled locations in a grid over Europe (Level I), subjectively located intensive monitoring plots (Level II) as well as a number of additional, national sets of plots. These large data sets together with data on stress factors have been used in various research projects in order to test hypotheses on cause-effect mechanisms of forest health.

The gathering of crown condition data into international data bases, have enabled studies at the European scale. Unfortunately, although the participating countries have co-operated closely on harmonization of methods, methodological differences apparently accounted for more than 30 % of the variation in defoliation on the Level I-plots in the period 1987–1997 (Klap et al. 2000). Main results of the analyses are shown in Table 3. Apart from the variable ‘country’, which explained variance 32 to 39 % of the variation (except for Holm oak (Quercus ilex)), the most important predictor variable was age (except for Q. ilex), which adds 2 to 14 % to the explained variance. Further step-wise completion of the statistical model revealed drought (relative transpiration) and various air pollution variables (acid deposition, air concentrations of SO2 and NO2 and AOT60 of O3) as significant predictors for the variation in crown condition, but the additional variance explained by the various selected predictor variables was generally very small: around 1 % (Table 3).

Table 3 Explained variance by country, age, altitude and environmental variables (acid deposition, relative transpiration, air concentrations of SO2 and NO2 and AOT60) by different ‘best’ selected statistical models (After Klap et al. 2000)

However, despite this, statistical relationships have been found in a number of studies. In the above mentioned 10-year study of the Level-I plots regression analyses, significant correlations were revealed between crown condition and factors considered as natural, i.e. tree age, altitude and drought stress, while indications of effects by pathogenic fungi and insects were also found. In addition, significant relationships were also found for air pollution, most notable for O3, but also for NOx, SO2 and acid deposition, suggesting negative effects in parts of central Europe, and in particular for deciduous tree species. The relationships were weaker for conifers, in particular for spruce (Klap et al. 1997, 2000; Müller-Edzards et al. 1997).

In a multiple regression analysis based on measured data at 222 Level II plots, De Vries et al. (2000a) found that the measured SO4 throughfall deposition was significantly related to increased defoliation of Norway spruce. In a similar data set study, De Vries et al. (2003d) used 262 Level II plots and found that 30–50 % of the variation in defoliation in pine, spruce, oak and beech could be explained by the variation in age, soil type, precipitation, N and S deposition and foliar chemistry (Table 4). The relationship with age was generally strong, except for pine. Positive correlations between defoliation and deposition of S and N were found for spruce and oak, while this was weak or absent for pine and beech. The defoliation of spruce and oak appeared to be larger in poorly buffered sandy soils compared to well-buffered clayey soils, while relationships with other factors such as precipitation and N deposition was variable.

Table 4 Overview of predictor variables explaining defoliation of the four most represented tree species of the Intensive Monitoring Plots with the number of plots (N) and the percentage accounted for (R2 adj) (after De Vries et al. 2003d)

In studies by Lorenz et al. (2003) and Fischer et al. (2003), based on about 3000 Level I plots, sulphur deposition was related to either high or increasing defoliation. However, significant relationships were found for Scots pine only. In addition, the studies indicated negative effects of insect attacks, drought and age. In a 20-year study Fischer et al. (2007) found some common, European-wide temporal variation of defoliation for most tree species, with a peak in the mid-1990s and an increase in 2003 and 2004. This might indicate that climatic conditions have a considerable influence, with particular attention to the summer droughts in the early 1990s and in 2003. Until now, O3 exposure has not been included in European scale crown condition studies due to a lack of data on the European wide scale.

Ozone exposure

From 2001 to 2005, Schaub and Sanz (2005) confirmed 108 symptomatic species at 136 EU and UN-ECE/ICP Forests Level II monitoring plots across 10 countries, including 35 tree-, 27 shrub- and 46 perennial herb species (Schaub et al. 2002). Temperate plots located in humid areas such as the Italian Po Plane and Canton Ticino in southern Switzerland had the highest number of symptomatic species but no relationship was found between O3 exposure, frequency of symptomatic sites and frequency of species with symptoms (Lorenz et al. 2005).

Ferretti et al. (2007b) evaluated whether O3 is a significant predictor of defoliation of beech by relating crown defoliation data with measured O3 concentration data, collected at a number of long-term monitoring plots in France, Italy, Spain and Switzerland within the framework of the Level II monitoring network (see also Ferretti et al. 2007a; Sanz et al. 2007; Gerosa et al. 2007). Ozone seems to have had a significant role (p = 0.008) in determining mean defoliation of beech plots. When year-by-year deviations from predicted defoliation are considered, a possible effect of O3 was obvious only at high exposure levels. Therefore, this raises the question whether one O3 critical level to protect vegetation is appropriate for all European countries (e.g. De Santis 1999). Ozone was selected as a significant predictor of defoliation (multivariate approach) and high O3 exposure was associated with deviation from predicted crown defoliation (temporal autocorrelation approach). This suggests that O3 may have some role in determining the overall defoliation status of the sites (which however is mostly controlled by site characteristics) and in determining annual deviations from predictions (which is mostly controlled by other biotic and abiotic factors) (Ferretti et al. 2007b).

National scale studies

National scale studies relating crown condition variables to environmental conditions were mostly limited to a selected set of variables, focusing on acid deposition, O3 exposure or climatic effects, as summarized below.

Acid deposition

West-Germany was the first country to carry out a nationwide assessment of crown condition, in 1983. The results of this were interpreted as evidence for air pollution effects (Schöpfer and Hradetzky 1984). This was, however, based on weak and indirect relationships such as a slightly higher defoliation in multi-layered stands than in one-layered stands, because the former stand structure was believed to intercept more with the moving air masses and hence capture more air pollution. Later, based on 61 Level-II plots with 4 tree species in Germany, Seidling (2004) found no clear indication of acid deposition effects on crown condition, by showing that high defoliation during 1995–2000 was statistically related both to low concentrations of soil Ca as well as low concentrations of SO4 in soil solution. In a study on 1800 Level-I plots in Germany, defoliation (taken as the fraction of trees with >25 % defoliation) was positively correlated to a number of S variables (S deposition; exceedance of critical loads of S and foliar S concentration) and foliar Ca concentration, while it was negatively correlated to N deposition (Augustin et al. 2005).

Early studies from Great Britain did not indicate negative effects of acid deposition on crown condition (Innes and Whittaker 1993; Mather et al. 1995), and the highest defoliation was found on sites having the lowest Al/BC ratios in soil water (Freer-Smith and Read 1995). However, a cause-effect relationship was considered unlikely, since plots with the highest loads of acid deposition were mainly located at high elevations exposed to mechanical damage from snow and wind (not quantified in the study) (Mather et al. 1995).

In Norway, two apparently contradictory results were found concerning the effect of acid deposition on Norway spruce defoliation. Both studies worked in southern Norway, which is a very appropriate region for correlative studies, because it contains a gradient from the highest deposition loads in Norway down to almost background values. Nellemann and Frogner (1994) found a relationship between defoliation and modelled (PROFILE) exceedance of critical loads based on about 100 Level-I plots, while Solberg and Tørseth (1997), who used a stratum of about 100 old forest officers plots, found no relationships with either deposition or soil chemical variables. While Solberg and Tørseth (1997) used standard statistical analyses methods, Nellemann and Frogner (1994) used an unusual statistical approach: Their 100 observations were first aggregated to 10 decile values, and for these 10 data points they fitted an S-shaped line and presented an R2 value of 0.94. The S-shaped fit was the only way to resolve the fact that the plots with the highest exceedance had defoliation similar to the plots with the lowest exceedance; however, they did not provide a plausible explanation for this curve. Later, Solberg et al. (2002) carried out modelling of exceedance of critical loads of acidity for the same 100 forest officers’ plots, and found the exceedance variables to not correlate with crown condition variables. Their soil solution Al/Ca were mostly below the criterion of 1.0, even though critical acid loads were exceeded at many plots, as calculated by the MAGIC and PROFILE models, demonstrating that modelled exceedance does not mean that unfavourable soil solution chemistry is present today but only that it might be reached in the future. Hence, exceedance may not be an appropriate variable in correlative studies searching for acid deposition effects.

Augustaitis et al. (2007) analysed Scots pine defoliation and air pollution data from 45 monitoring plots during 1994–2004 in Lithuania. During these years, air pollution decreased considerably. For example the SO2 concentrations decreased by 80 %. At the same time defoliation decreased considerably. This co-incidence between two decreasing trends produced positive correlations between defoliation and air pollution variables, and they concluded from this that air pollution were among the key factors affecting defoliation. However, this seems questionable as the SO2 concentrations were only in the range 1.6–6 μg m−3, which is clearly below critical levels, while it is unlikely that the ecosystem could recover from an eventual soil acidification stress during only 10 years. In addition, their hypothesis testing was questionable, since they treated all data from 45 stands during 11 consecutive years as 495 independent observations, i.e. not taking temporal autocorrelation into account.

In the Netherlands defoliation has been found to correlate mainly to tree species and age, and to a much lesser extent to foliar N content and pH of the soil solution (Hendriks et al. 1994). In a later study Hendriks et al. (2000) found that a range of factors such as deposition, climatic, stand, site and biotic factors was correlated to defoliation, however, the relationships were weak.

Ozone exposure

Augustaitis and Bytnerowicz (2008) investigated the possible effect of natural and anthropogenic environmental factors on Scots pine defoliation and stem growth in Lithuania, including acid deposition and O3 exposure by applying a linear multiple regression technique. Results confirmed earlier findings that acidifying air compounds and their deposition are key factors resulting in pine defoliation changes (Augustaitis et al. 2003, 2005, 2007). They explained from 23 % to 28 % of the variance of residual defoliation of pine trees. The effect of peak O3 concentrations was less significant (19.3 %), however, they increased the explanation rate of defoliation residual variability by air concentration of acidifying compounds, their concentration in precipitation and deposition from 3 to 8 % (Augustaitis and Bytnerowicz 2008).The findings of that study revealed that the contribution of peak O3 concentrations to the integrated impact of the other environmental factors on Scots pine defoliation and basal growth increment remain statistically significant, confirming results obtained in Finland (Utriainen and Holopainen 2000) and the Carpathian forest (Muzika et al. 2004). Also, the effects of peak O3 concentrations have higher significance for spruce trees in Norway than the mean diurnal O3 concentration or the AOT40 index in southern Sweden (Karlsson et al. 2006). Therefore, it seems that peak O3 concentration is one of the key factors affecting forest ecosystems in North-eastern and Northern Europe (Augustaitis and Bytnerowicz 2008).

With emphasis on the simulation of the local water balance, and based on forest condition data for 1983–1985 from the Swiss Forests Inventory (Stierlin et al. 1994), Zierl (2002) investigated possible relationships between leaf/needle loss and O3 concentrations in Switzerland by means of statistical analysis, based on the hypothesis about the possible impact of O3 on crown defoliation. All results from the temporal and spatial statistical analysis support the hypothesis that tropospheric O3 contributes to the forest decline observed during the last decades in Switzerland. However, under different climatic conditions plants might respond differently to O3, which results in regional and yearly differences in sensitivities of forest stands. Especially drought stress during the summer months can provide an effective protection from O3 damage by enforcing stomatal narrowing or even closure. At dry forest sites the risk of ozone damage even seems to be highest in years with moderate O3 concentrations, but small drought stress.

Baumgarten et al. (2009) compared exposure and flux-based indices of O3 risk at 13 Level II monitoring sites across Bavaria, Southern Germany from 2002 to 2005. Both, the exposure-based indices AOT40 and MPOC (Maximum Permissible Ozone Concentration, Grünhage et al. 2001; Krause et al. 2005) as well as the flux-based POD indicator suggested that Bavarian forests are at risk from ozone. However, O3 visible symptoms were found only to a small extent (<3 % of total leaf area, >1–5 % of leaves) on beech trees throughout the study years, and not on all plots.

In a study dedicated on the impact of O3 on forest ecosystems in Italy, O3 was found to be a predictor of crown transparency residuals in beech sites over two consecutive years, but the variance explained amounts to less than 10 % (Bussotti and Ferretti 2009). In a more recent study, Ferretti et al. (2012a) investigated the relation between defoliation data (2007–2009) collected at ICP Forests Level I (n = 15) and Level II (n = 1) plots in Trentino, Italy and site and environmental factors including O3, using multiple regression models and Linear Mixed Models (LMM). Both approaches showed that defoliation values increased with increasing frequency of reported damage (biotic and abiotic) and with decreasing level of foliar N:K. Despite the recorded high values, O3 was never reported as a significant predictor of defoliation at the plots in Trentino. These results were further confirmed for the 2000–2009 period, taking into account AOT40 and stomatal flux, at the Level II plot.

Climatic factors

Weather conditions, especially drought, are major drivers for year to year variations in crown condition (e.g. Lorenz et al. 1998). In Britain, Mather et al. (1995) analysed data from 294 monitoring plots together with modelled climatic and air pollution variables between 1989 and 1992, and concluded that drought was the major factor affecting the crown condition of trees. Using multivariate analyses techniques defoliation could be attributed to soil moisture deficit and potential evapotranspiration, and the cause-effect inference was supported by observations of drought stress effects on trees in the very dry 1989 summer, particularly in southern Britain. However, drought and other factors together explained only 5–10 % of the variation in defoliation. Solberg (2004) found that summer drought was followed by increases in defoliation, discolouration, more cones and mortality, working on 450 forest officers spruce plots in southeast Norway. Increased defoliation occurred the year after the drought, which was interpreted as a result of increased needle fall in the autumn after the dry summer. Seidling (2007) studied German defoliation data from 1990–2004, and found increased defoliation in pine, spruce, beech and oak after dry summers, mainly with a delay of one year. The response was particularly evident after the extremely dry summer of 2003, which is also supported by the annual German forest condition report from 2004 which showed that their 20 year time series of spruce defoliation peaked in the following year 2004 (Anon 2005). In the Swiss party of the Rhone Valley, a considerable decline of Scots pine has been recorded, where 59 % of the trees died from 1994 to 2004, most of following the drought periods 1996–1998 and 2003–2004 (Dobbertin and Rigling 2006). Their analyses indicated that the decline was caused by an interaction between drought and mistletoe (Viscum album) infection, because the mistletoes caused increased water use. In Italy, Bussotti et al. (1995) demonstrated closely related temporal variations in litterfall and defoliation in a Holm oak monitoring plot, and concluded that summer drought was a major driver for defoliation, possibly in combination with other stress factors.

Biotic agents

Many biotic agents are likely to be responsible for discolouration and defoliation, and some well-known examples of insects include Elatobius sp. lice attacks on Sitka spruce (Picea sitchensis) in Britain, oak caterpillars (Tortrix spp.) in central Europe, various moths on mountain birch (Betula pubescens) in Fennoscandia, and pine sawflies (Diprion pini and Neodiprion sertifer) in Fennoscandia. An example of fungal diseases is Gremmeniella abietina, which had an epidemic in Sweden and Norway in the years 2001–2003, causing considerable browning and defoliation. The extent of the disease was well described using data from grids of plot in the national forest inventory and the national forest damage inventory, even though these grids were sparse (Wulff et al. 2006). And as mentioned above, a widespread mistletoe infection apparently contributed to a Scots pine decline in the Swiss alps by causing excessive water use and more drought stress (Dobbertin and Rigling 2006).

Other stress factors

In general, correlative studies often indicated a multitude of stress factors. For example, based on the longest available time series of beech and oak defoliation in several of the German federal states, Eichhorn et al. (2005) used logistic regression analyses which demonstrated that defoliation was related to a number of factors, i.e. age, crown spacing, stand composition, fruit bearing and insect defoliation. Tree age is often the dominating explanatory variable, as for example found in a German study: Zirlewagen et al. (2007) found that age explained 35–64 % of the variance in defoliation, which is in line with the findings of Seidling (2007). However, for cause-effect inferences, age cannot unambiguously be held responsible for defoliation, and the often identified relationship between defoliation and age may in fact represent an interaction between stress factors and age: In a study on the relationships between defoliation and age in Norway, it was demonstrated that the relationship varied from region to region, and in one region of Norway there was almost no increase in defoliation with increasing age (Solberg et al. 2009). This was the western part of Norway, where Norway spruce has been introduced, and where the soils are free from the severe root rot Heterobasidiuon annosum and where drought stress is rare due to abundant rainfall. Hence, when a large fraction of the variation in defoliation often is attributed to age, this may in fact camouflage stress affects that hit older trees harder than younger trees.

Forest growth and tree carbon sequestration

European scale studies

The influence of site and environmental factors on measured forest growth data was investigated at tree level by Laubhann et al. (2009) and at stand level by Solberg et al. (2009), using five-year growth data for the period 1994–1999 for nearly 400 intensively monitored forest plots. Evaluations focused on the influence of nitrogen and acid deposition in interaction with forest stand characteristics and weather variables. The impact of O3 exposure was not included. The study included the main tree species common beech, sessile or pedunculate oak (Quercus petraea and Q. robur), Scots pine and Norway spruce. Requirements for plot selection were different for both methods, resulting in 382 plots for the individual tree-based approach and 363 plots for the stand-growth model approach. The individual tree-based model had measured basal area increment of each individual tree as a growth response variable and tree size (diameter at breast height), tree competition (basal area of larger trees and stand density index), site factors (e.g. soil C/N ratio, temperature), and environmental factors (e.g. temperature change compared to long-term average N and S deposition) as influencing parameters. In the stand-growth model, stem volume increment was used as the growth response variable, after filtering out the expected growth. Expected growth was modelled as a function of site productivity, stand age and a stand density index. Relative volume growth was then calculated as actual growth in % of expected growth. Concerning the meteorological variables temperature, precipitation, and drought the deviations from the long-term (30 years) means were used.

Results of the multivariate analyses at individual tree level, shown in Table 5, indicated a 1.2–1.5 % increase in basal area increment (coefficients varying between 0.012 and 0.015 relative increase), depending on tree species in response to a fertilizing N deposition effect of 1 kg N ha−1 yr−1. In this case, the response was significant for all included tree species. Laubhann et al. (2009) also showed that volume increment was proportionally related to basal area increment. Referring to a total carbon uptake for European forests of 1729 kg C ha−1 yr−1 by De Vries et al. (2006), they then estimated the average response in terms of C sequestration between 20.7 and 25.8 kg C per kg N, depending on tree species composition. Increasing temperature was also found to have a positive influence on forest growth, but this effect was less significant. Results of the multivariate analyses at stand level are shown in Table 6. The results indicated roughly a 1 % increase in site productivity in response to a fertilizing effect of N deposition of 1 kg of N ha−1 yr−1 for Scots pine and 2 % for Norway spruce. Stronger responses were seen for N sensitive sites (high soil C/N ratio), having roughly a 1.3 and 2.2 % increase in growth for pine and spruce, respectively, in response to a fertilizing effect of 1 kg N ha−1 yr−1. These responses for pine and spruce were recalculated in terms of C sequestration, by taking the product of the measured mean annual volume increment times the mean wood density times the estimated growth increase, assuming a C content of 50 %. In converting volume increment to tree C changes, the increase in branches and woody roots is accounting for. The means of the modelled C sequestration were 19 and 38 kg C/kg N. More information on the approach used and results obtained is given in Solberg et al. (2009).

Table 5 Multivariate regression results at tree level indicating the relative change in basal area increment per unit change in influencing factor of the four most represented tree species at the Intensive Monitoring Plots. Values for N deposition are given in bold. The regression results indicate the relative change in stem volume growth per unit change in influencing factor. For example, a value of 0.013 for N deposition implies an increase of stem growth of 1.3 % per kg N deposition Note: - implies that the effect was statistically insignificant (p > 0.05) (after Laubhann et al. 2009)
Table 6 Multivariate regression results at stand level indicating the relative change in stem volume growth per unit change in influencing factor of the four most represented tree species at the Intensive Monitoring Plots. Values for N deposition are given in bold. Note: - implies that the effect was statistically insignificant (p > 0.05). (after Solberg et al. 2009)

S deposition and acidic deposition both had mostly positive slopes, but these two explanatory variables were highly correlated to N deposition. Signs for the parameter estimates for temperature deviation gave significant positive parameter estimates for pine and spruce and partially for beech, while estimates were not significant for oak. Estimates for drought were generally negative for spruce pine and beech and again mostly insignificant for oak.

National scale studies

As with crown condition, national scale studies relating forest growth to environmental conditions were mostly limited to either air quality (acid deposition and/or O3 exposure) or climatic effects, as summarized below.

Air quality

The possible influence of nitrogen and acid deposition on forest growth was evaluated in a study with nationwide data-sets for Norway (Solberg et al. 2004). Results indicated a 1–2 % increase in growth per kg N annual deposition for pine and spruce, but for beech and oak, the response was not significant. Augustaitis and Bytnerowicz (2008) investigated the impact of acid deposition and O3 exposure on both defoliation (see before) and stem growth of Scots pine stands in Lithuania. Results showed that the integrated impact of acidifying compounds and their deposition accounted for 18.5 % of variability in stem increment residuals, while O3 increased the degree of the explanation by 8.7 %–27.2 %.

The impact of O3 on forest ecosystems is specifically monitored and evaluated in Italy, where O3 levels are measured in 30 plots using passive samplers together with forest response parameters, including crown condition (transparency), BAI (basal area increment) and visible symptoms on spontaneous vegetation (Bussotti and Ferretti 2009). Results showed that levels of AOT40 are above the concentration-based critical level of 5 ppmh in all sites, but the evidence of impact on forest vegetation remains limited. As mentioned before, O3 was found to be a predictor of crown transparency, but the variance explained was less than 10 % and the relation between BAI reduction and O3 was even less certain (Bussotti and Ferretti 2009).

As with defoliation, Ferretti et al. (2012a) also investigated the relation between growth data, reported as relative Basal Area Increment (BAI), collected in the period 2005–2009 at ICP Forests Level I (n = 15) and Level II (n = 1) plots in Trentino, Italy and site and environmental factors, using multiple regression. Results showed that relative BAI increased with foliar N:Mg and decreased with tree diameter. As with defoliation, O3 was not a significant predictor of BAI at the sites in Trentino. These results were further confirmed for the 2000–2009 period, taking into account AOT40 and stomatal flux at the Level II plot. Large variations in exposure and flux did not result into consistent variation in growth. In a recent investigation carried out in Italy (Ferretti et al. 2012b), variables related to management, climate, tree health, and N in soil, foliage and throughfall were the most important predictors of relative growth of Level II plots.

Braun et al. (2007) applied an epidemiological analysis across 83 beech and 61 Norway spruce plots across Switzerland and associated a 7.4 % reduction in shoot length with an AOT40 level of 0–10 ppm h for beech but shoot growth of Norway spruce was not affected by O3. Karlsson et al. (2006) used relative annual basal area increment of mature Norway spruce trees in south-central Sweden during 9 years and found a decline in basal area increment of 4.6 % in 19- to 35-year old Norway spruce trees that could be correlated with differences in ambient O3 exposures between years in the range of 2–9 ppm h. At the same site, O3-specific symptoms were demonstrated by means of light electron microscopy from mature Norway spruce trees (Kivimäenpää et al. 2004).

Climatic factors

Various national scale studies showed that soil water availability and also temperature change was correlated with forest growth. By combining tree-ring analysis with models for soil water availability, that had been locally calibrated, Lebourgeois et al. (2005) showed that soil water deficit, in particular in early summer, could well explain annual and early wood growth fluctuations for common beech on French Level II plots. For Norway spruce, early and latewood growth depended mainly on current summer soil water availability, while for silver fir early wood growth was also influenced by the water supply of the previous year’s summer and autumn (Lebourgeois 2007). Using annual stem circumference measurements at Swiss Level II plots, a tree growth reduction was found at low altitude in 2003 as compared to the wet years 2002, while high altitude forests showed no growth change or increased growth (Jolly et al. 2005). The reason, was that at low altitude soil water availability was highly reduced in 2003 and correlated well with the observed growth reductions, while water was not limiting at high altitude (Graf Pannatier et al. 2007). At high altitude, the temperature during the summer of 2003 was increased favouring increased tree growth (Graf Pannatier et al. 2007).

Species diversity of the ground vegetation

European scale studies

Relationships between species composition of the ground vegetation and environmental factors were evaluated at 366 Intensive monitoring plots, as presented in van Dobben and de Vries (2010). Included environmental factors were: (i) site characteristics, such as climate zone, altitude, tree species and soil type, (ii) climatic data, such as temperature and precipitation, (iii) chemical soil data, including element contents in the humus layer and the mineral topsoil both in soil solid and liquid phase, and (iv) bulk deposition data of NH4, NO3, SO4, Ca, Mg, K, Na and Cl. The data used in the analysis are annual mean values over the period 1995–1999, being the period in which the ground vegetation data were assessed. Bulk deposition collected in nearby open field precipitation samplers was used as an estimator for atmospheric deposition. Although throughfall may give a better estimation of the deposition to the forest floor, these data were not used because they were available for a limited number of plots only. A separate analysis was carried out for these plots to check the effect of throughfall, which showed that impacts inferred on the basis of bulk deposition data were very similar to those inferred on the basis of throughfall data.

Table 7 is a summary of an ordination analysis of the first observations on all Intensive Monitoring plots available up to 1999. Out of the 14.5 % of the variance in the species abundances that can be explained by the available predictors, by far the largest part is due to the ‘traditional’ factors soil, climate and tree species, and only 0.6 % is due to bulk deposition. Na and NO3 are the only ions in bulk deposition that have a significant effect, and of these two only NO3 is of anthropogenic origin. In addition, 7 % variance is due to the countries themselves, probably as a result of methodological differences. It should be noted that in ordination (i.e., multivariate statistics), percentages explained variance are intrinsically low. The ordination diagram (not shown here) indicates the following main causes of variation in the vegetation (in order of decreasing importance): (1) soil fertility, (2) altitude, (3) tree species (coniferous vs. deciduous), (4) climate (oceanic vs. continental). The effect of deposition is small but statistically significant and ecologically understandable, i.e. there is a slight tendency of nitrophilic species to increase with increasing NO3 deposition.

Table 7 Percentage explained variance of plant species abundances at Intensive Monitoring Plots that could be ascribed to four main groups of variables in a European and Dutch study (after Van Dobben and de Vries 2010)

Conclusions similar to the above, based on the same one time survey of ground vegetation data, are drawn by de Vries et al. (2003a). These authors stress the bias in monitoring results caused by methodological differences between the countries, especially if biodiversity indicators like the number of species per plot are used. Such bias can be due to differences in e.g. plot size, number and distance of subplots, number of yearly visits, size of the field team and experience of fieldworkers. In general, analyses based on composite measures like mean indicator values per plot or measures derive from multivariate statistics are considered more robust than analyses based on individual species. Seidling and Fischer (2008) and Seidling et al. (2008) present an analysis of a subset of the European Intensive Monitoring plots based on Ellenberg indicator values (Ellenberg et al. 1991). Also their analysis leads to the conclusion that there is a slight but statistically significant increase of nitrophilic species (i.e., an increase in mean Ellenberg nutrient indicator) with increasing nitrogen deposition.

It should be noted that small temporal shifts in the vegetation due to slow changes in soil chemistry caused by deposition will not be detected in transversal studies over a wide geographical range like the ones presented above. Van Dobben et al. (2010) presented a new analysis of the European Intensive Monitoring plots, based on repeated data on ground vegetation collected up to 2006. In this extended data set the effect of N deposition (in this case derived from EMEP modelling results) on ground vegetation is more pronounced than in the earlier analyses. Table 8 is a summary of an ordination analysis of the repeated observations on Intensive Monitoring plots available up to 2006, indicating the explanation of temporal changes in the vegetation and environmental variables. The total explained variance was 8.2 % only, but there appeared to be a significant effect of N deposition, entailing a clear shift in species composition toward nitrophilic species, explaining 1.9 % of the total explained variation of 8.2 %. This shift occurs both in space (i.e., when moving from ‘clean’ areas to areas with a high deposition), and over time (i.e. when comparing the first and the last observation of plots where time series are available).

Table 8 Explained variation of the change in vegetation at Intensive Monitoring Plots, using EMEP estimates to quantify N deposition. Change is determined as the difference last - first relevee of each species in each plot where the time lag between the first and last relevee is >6 years (after Van Dobben et al. 2010)

National scale studies

Although many countries have forest monitoring schemes at the national level, only few of them have their results published in scientific literature, the most recent examples being from the Netherlands (Van Dobben and de Vries 2010) and Norway (Økland et al. 2004). The results of these two national evaluations appear to be very comparable to the above one on the European level. The result of the study in the Netherlands is also given in Table 8, which shows that again there is a significant but small contribution of deposition to the explained variance. In this case the components of deposition with a significant effect are SO2 and Mg. Of these, only SO2 is of anthropogenic origin, and Mg is probably just an indicator for the distance to the coast. The Norwegian study (Økland et al. 2004) also found the traditional factors (soil, climate, humidity) to be the most important determinants of the vegetation. In this study soil acidification due to atmospheric deposition is considered as one of the main drivers of vegetation changes (mainly a loss of species in areas with high deposition). However, this conclusion is only based on inference as the vegetation was not directly related to measured or simulated deposition.

All results referred to so far are based on the analysis of generalised indicators, i.e. derived from multivariate approaches (ordination axes) or indicator values. An alternative approach may be the analysis of the response of single species to environmental factors using logistic regression (e.g. De Vries et al. 2003d). However, with studies over a large geographical extent the distribution area of each species should also be taken into account. Failure to do so (like in the above studies) may lead to spurious results because outside its distribution area a species will be absent anyway, even if abiotic conditions are favourable. The results of such studies are therefore hard to interpret.

Discussion and conclusions

Explicit assessments of the impacts of different air quality variables (deposition of N and acidity and O3 exposure) and weather variables (temperature and precipitation) on forest ecosystems are difficult to derive on the basis of monitoring data because of the various interacting effects. This holds less for the impacts on forest nutrient status and soil (solution) chemistry but is specifically true for the impacts on forest health, forest growth and species diversity of the ground vegetation. Below, we discuss the strengths, weaknesses and perspectives of the European forest monitoring programme, followed by a discussion of the main results of this programme in view of those strengths and weaknesses and knowledge obtained from field and laboratory research and we end with main conclusions that can be drawn from the monitoring results.

Strengths, weaknesses and perspectives of the European forest monitoring programme

The main strength of the ICP Forests monitoring programme is the large number of variables that has been measured at each plot, the large spatial extent and the temporal coverage over the last approximately 20 years. Furthermore, compared to other monitoring programmes, ICP Forests has invested great efforts into harmonization. This has led to various manuals for the data assessment, including quality assurance and quality control procedures.

An important weakness of the forest monitoring programme is the use of non-specific response indicators, such as tree defoliation, which are affected by various interacting air quality and climate parameters. In combination with correlations between those drivers, such as those between NOx deposition and O3 exposure, this causes the problem of confounding factors in regression relations. Furthermore, despite the focus and efforts to come to harmonized approaches, there is still a lack of standardization for various assessments, specifically for crown condition and ground vegetation data, eventually causing substantial observer bias. Finally, since the programme is extremely labour intensive, there has until recently been a very strong focus on the actual field measurements all the way through to data management and data storage at the expense of actual data analyses and dissemination of results.

Despite the limitations of confounding factors due to non-specific response indicators, correlations between drivers and methodological inconsistencies, there are large potentials to explore the wealth of real-field data obtained from of large-scale long-term monitoring networks such as the Intensive Monitoring Level II network of ICP Forests. Potentials and perspectives, in combination with more advanced approaches and in interaction with other programmes, include:

  • Integration and further data analysis: The ICP Forests Monitoring Programme in Europe has strong links to other ICPs such as ICP Modelling & Mapping, ICP Integrated Monitoring and ICP Vegetation allowing for more integrated and interdisciplinary analyses across Europe. Another source for future potential lies in a stronger collaboration with programmes such as ICOS, EMEP, LTER, etc. Until now, integrated analyses of the data are still limited while the wealth of data allows for such an analysis, especially considering the current time coverage.

  • Model validation: The monitoring data set also becomes more and more attractive for validation purposes of regional and global scale models. This is particularly true for physiological and environmental parameters describing the water, carbon and nutrient cycles.

  • Introduction of new techniques and approaches: The technical advances allow more and more parameters to be monitored at higher, i.e. hourly resolution which is another great potential for model validation purposes. Furthermore, there are ideas to include not only Level II & I but also so-called Level III (super) sites which are highly equipped long-term monitoring plots that allow integrated studies and analyses from cell to ecosystem level and validation of detailed process based models. Finally the monitoring systems can be used as a ground validation systems for remote sensing techniques, that are used to assess various forestry aspects, such as growth, crown condition and tree nutrition.

Impacts on forest ecosystems as derived from monitoring studies

Nitrogen deposition effects on forest growth and plant species diversity

The two presented European scale studies, although being different in their approaches, using a different set of plots and different methods to estimate the N induced carbon uptake in stem wood, resulted in similar results. In summary, the results of this wide European growth dataset in the period 1993–2000 indicated a clear fertilization effect of N deposition on European forests varying between approximately 20–40 kg C per kg N (see also De Vries et al. 2008) although higher growth responses to N deposition (up to 4 %) have also been observed in selected Swiss observation plots (Braun et al. 1999). The results are in line with approaches focusing on the fate of N in forest ecosystems (e.g. Nadelhoffer et al. 1999), results of N fertilizer experiments and model simulations (De Vries et al. 2009), as well meta- analysis of various N addition experiments. In the correlative study, S deposition and acidic deposition were also positively correlated with forest growth, but both variables were highly correlated to nitrogen deposition. It has to be realized that the overall positive effect of N deposition on forest growth does not preclude negative effects in extreme cases of acidification. Especially the removal of base cations from the adsorption complex of moderately acid soil may have negative effects on the long term. Lawrence et al. (2005) compared changes in soil chemistry, based on archived soil samples in 1926, 1964, and 2001, with patterns of tree growth in the same period near St. Petersburg, Russia. Results showed a depletion of exchangeable Ca and Mg in the period between 1926 and 1964 and an increase in exchangeable Al between 1964 and 2001. The onset of Al mobilization coincided with decreased diameter growth and a suppression of climate-tree growth relationships in Norway spruce, suggesting that the trees were responding to a decline in Ca and Mg availability. However, overall, the impact of acid deposition on growth does not seem to be negative at European scale.

Even though large impacts of N deposition on forest ground vegetation have been found in local N addition studies, such an effect is hard to assess in large-scale (i.e. Europe-wide or at least country-wide) monitoring studies relating variation in environmental factors, including N deposition, to variation in plant species diversity. These studies lead to the conclusion that factors that are considered to be key indicators for the vegetation in phyto-sociological literature (e.g. Braun-Blanquet 1964; Grime et al. 1988; Oberdorfer 1979; Ellenberg et al. 1991) are most important. These factors are soil chemistry, tree layer characteristics (most notably the tree species), climate and hydrology. Effects of N deposition may be significant but are always small.

A disclaimer with large-scale monitoring studies is that they may be strongly hampered by methodological inconsistencies, especially if they span more than one country. The European ICP Forest network has large differences in plot size and number of subplots between the countries and sometimes even within a country. Observers may have different levels of species knowledge (Ringvall et al. 2005), and in transversal studies their species knowledge tends to increase over time (De Vries et al. 2003a). Floristic views may differ per country or even per observer. Plots may be fenced to exclude the effect of grazing or unfenced, observation dates relative to the growing season may be different and there may be different number of visits per year. On the abiotic side analytical techniques may be different. Anyway evaluations over the complete geographical extent are only possible for variables that have been measured in all plots, and this requires a level of standardisation that is hardly attainable. Mosses and lichens, which are probably the best indicators for environmental conditions (Van Dobben et al. 1999), are not recorded by all countries within ICP Forests. These methodological problems make the large-scale monitoring rather insensitive, i.e. only environmental variables with a rather large effect are likely to be significant in such a study, despite large numbers of plots and a large geographical extent. Large-scale monitoring studies with repeated data, however, are likely to pick up relationships between N deposition and trends in plant species diversity indicators.

Acid deposition effects on soil and forest condition

Results of monitoring studies at the European scale, presented above, show that elevated S and N deposition is reflected by an increase in S and N concentration in forest foliage and elevated SO4 and NO3 concentrations in soil solution. This in turn leads to soil acidification in terms of elevated Al and/or Bc leaching from the forest ecosystem and increased dissolved Al and/or Bc concentrations. Although on average N deposition is higher than S deposition, it still leads to less cation mobilization (acidification) because SO4 generally behaves conservative, whereas N is strongly retained in the soil and/or denitrified. In combination with results from country scale monitoring, the studies also show that reductions in S and N deposition observed over the last 10–15 years are reflected in a decrease of foliar S and N concentrations and dissolved SO4, NO3, Al and Bc concentrations. Since both Al and Bc decrease with a decrease in acid deposition, the Bc/Al ratio stays mostly rather constant.

Despite the impacts on forest nutrient status and soil (solution) chemistry, acidic deposition was found to explain only about 1 % of the variation in crown condition on a European scale, whereas most monitoring results at national scale do not even show any convincing statistical relationship between crown condition and acid deposition or sulphur dioxide exposure. These results imply that there is no clear evidence that acid deposition negatively affects forest health, in terms of crown condition, on a large regional scale. This result is in line with the observation that defoliation tends to oscillate, which indicates that acute stress and recovery are drivers in forest health, rather than chronic air pollution stress (Landmann 1993; Mather et al. 1995; Kandler 1988). In the past several regional studies have concluded that air pollution effects on crown condition are clearly evident, but the provided evidence for this is generally weak and questionable. For example, the first nationwide assessments of crown condition in West Germany in 1983 was interpreted as evidence for air pollution effects based on weak and speculative relationships (Schöpfer and Hradetzky 1984), and the strong increase in defoliation observed in the next year was seen as evidence for the increasing damage by air pollutants, but this increase could well be explained by a change in the assessment methods (Kandler 1992).

The limited evidence of acidification impacts on forest condition from monitoring study results is in line with the earlier presented evidence that the linkage between Bc/Al ratios and root or shoot growth found in laboratory experiments (e.g. Sverdrup and Warfvinge 1993), is generally not found in the field situation (De Wit et al. 2001a, b, 2010). Effects are even lacking in extreme situations. For example, the Uusikaarlepyy monitoring site in Finland has a very acidic soil with extremely high concentrations of inorganic Al, while the spruce trees growing there show no defoliation (Rautjärvi et al. 2002). Also in studies from moderately polluted areas North America, indications of air pollution effects on crown condition have been ruled out in favour of other factors, such as climate, stand and site condition, and pests, with special regard to drought and insect attacks (Brooks 1994; Burton et al. 1991). Despite the large body of data gathered and a widespread research effort, the effects of long-range transported acid deposition on forest health remain indecisive (e.g. Skelly and Innes 1994; Solberg and Tørseth 1997; Thomas et al. 2002; Solberg et al. 2004; Binkley and Högberg 1997; Kandler and Innes 1995; Klap et al. 1997).

Altogether, an evaluation of correlative studies based on European crown condition assessments together with related studies in Europe and North America, indicate that negative effects from air pollution on forests in Europe are weak or cover minor parts of Europe. Widespread, negative effects from acid deposition and soil acidification are not detectable. If occurring, negative effects of air pollution on forests in Europe seems to be local (“classical smoke damage”) due to industry plumes containing high concentrations of SO2 (Lomský et al. 2002) but the extent of such damage has decreased considerably during the last 30–40 years and this SO2 damage type is today rare.

Ozone exposure effects on forest condition and forest growth

The assessment of visible O3-induced injury at monitoring sites across Europe has turned out to be a useful and inexpensive method to show spatial and temporal trends for potentially harmful O3 exposures (Sanz et al. 2007; Gerosa et al. 2007). However, species-specific differences in O3 sensitivity, differences in microclimatic conditions and species composition make it difficult to relate visible injury to O3 exposure and determine critical levels of O3 risk (Schaub and Sanz 2005). Most monitoring results both at national and European scale mostly showed a convincing statistical relationship between crown condition and O3 exposure, but generally, the effects seem to be limited. Until now, the possible impact of O3 on forest growth has hardly been evaluated in either national or European field scale monitoring studies with the exception of Augustaitis and Bytnerowicz (2008), who investigated the impact of acid deposition and O3 exposure on both defoliation and stem growth of Scots pine stands in Lithuania. Their results suggested that peak O3 concentration is one of the key factors affecting forest ecosystems in North-eastern and Northern Europe.

The limited results of correlative studies based on monitoring results do not yet confirm the hypothesis that ambient O3 is a key factor explaining spatial and temporal changes of tree crown defoliation and/or forest growth. The problem with routine forest monitoring is, however, the use on non-specific response indicators, such as tree defoliation and tree growth, which are subject to many other stressors than O3 (e.g. Percy and Ferretti 2004). Furthermore, there is generally a correlation between NOx deposition and O3 exposure, while severe drought conditions may limit O3 uptake by plants, therefore preventing the development of visible injury (Fischer et al. 2005). Consequently, field evidence of O3 effects on forests is limited and extrapolation of experimental results to the real forest condition is problematic (Ferretti et al. 2007c).

In this context, it has also to be realized that all correlative studies were made using either peak or average O3 concentrations or air quality standards such as AOT40 (Kärenlampi and Skärby 1996; Fuhrer et al. 1997). These O3 exposure metrics, however, do not account for environmental and biological variables that affect ozone uptake and plant injury (Matyssek et al. 2007a; Karlsson et al. 2007) and are thus quite inappropriate to predict ozone effects on regional scales. Until now, there are no publications relating monitored forest growth changes to changes in the phytotoxic ozone dose. This is a potential promising approach in correlative field studies. Critical levels in terms of phytotoxic ozone dose (POD) have until now been based on the results of experiments with seedlings under controlled condition (e.g. Karlsson et al. 2003) with relative reduction of biomass used as a response indicator without any results from field gradient studies or long-term monitoring studies. Such monitoring results are highly relevant since the extrapolation of results obtained from seedlings to mature trees under real forest condition has been severely challenged (Kolb and Matyssek 2003). Currently, a study is carried out in which the combined effects of climatic factors, nitrogen deposition and phytotoxic ozone dose (POD) are related to changes in basal area increment, using fifteen-year growth data for the period 1994–2009 for nearly 400 intensively monitored forest plots, as a follow up of the study by Solberg et al. (2009). Preliminary results (not yet published), however, show a large correlation between N and POD, thus obscuring the negative effects of cumulative O3 uptake.

Climate effects on forest condition and forest growth

Correlative studies of measured or modelled climatic and air pollution variables with monitored changes in crown condition at both European and national scale showed that climatic variables (especially drought) were generally the most significant predictors for the variation in crown condition. This result is in line with field studies and experimental studies, showing that drought can reduce foliage production and foliage size, increase foliage shedding, cause yellowing due to reduced foliar nutrient concentrations or cause increased attacks by pests and diseases. For example, it has been demonstrated that oak decline, which has received great attention during the last decades, is not new, but rather that it has occurred repeatedly during the past three centuries, and can be attributed to climatic extremes (winter frost, summer drought), as well as defoliating insects and pathogenic fungi, while air pollution and nutritional imbalance do not seem to be important (Thomas et al. 2002).

Both the European scale studies and various national scale studies also showed that soil water availability and drought stress are correlated with forest growth. At European scale, an increase in temperature deviation gave significant positive relationships with the growth of pine, spruce and beech, while increased drought had a significant negative impact on the growth of these tree species. The various national scale studies showed specifically a large impact of summer soil water availability on growth. The results are fully in line with the experimental observations that reduced water availability results in decreased tree growth and reduced foliage mass.

Overall conclusions

As mentioned above, there are various limitations to the analyses of large scale and long-term data sets obtained from the Intensive Monitoring Level II network of ICP Forests, such as the problem of confounding factors due to non-specific response indicators, correlations between drivers and methodological inconsistencies. Nevertheless, results of monitoring studies enabled some conclusions on large scale relations between forest ecosystem impacts and drivers of change that are in line with other field studies i.e.:

  • Reductions in N and S deposition observed over the last 10–15 years are reflected in a decrease of foliar N and S contents and dissolved SO4, NO3, Al and Bc concentrations. Since both Al and Bc decrease with a decrease in acid deposition, the Bc/Al ratio stays mostly rather constant.

  • Climatic factors and in particularly drought stress, appear to be primary drivers for forest crown condition variations, while air pollution effects appear limited in time and space.

  • Drought usually increases crown transparency and reduces tree growth, while temperature rise increases forest growth.

  • Ozone, especially peak O3 concentrations, are related to a decrease in crown condition and growth of forest ecosystems in parts of Europe.

  • Negative impacts of acid (N and S) deposition on crown condition and forest growth are hardly or not derived. Inversely, N deposition has significantly increased above ground tree growth in European forests.

  • N deposition has negatively affected the species diversity of the ground vegetation although the impact is limited compared to traditional factors such as tree species and climate.