Variables related to nitrogen deposition improve defoliation models for European forests

Abstract

Key message

Atmospheric deposition of nitrogen compounds and soil and foliar variables related to N deposition resulted important factors accounting for the variability of defoliation in European forest plots.

Context

Nitrogen (N) deposition has increased in the northern hemisphere because of anthropogenic-related emission of N compounds. Increased N availability may have an adverse impact on forest sustainability.

Aims

This study aims to test the importance of throughfall N (N_thr) deposition in explaining the variability of the frequency of trees with defoliation >25 % (F₂₅), an indicator of forest condition.

Methods

A pan-European data set (71 plots) with enhanced quality control was considered. The additive effect of N_thr-related predictors (identified conceptually and by rank correlation) in explaining F₂₅ was investigated by partial least square regression in comparison with a reference model based on site, stand, management and climate data. Reported damage to foliage, N_thr deposition, foliar N ratios and mineral top-soil variables were added stepwise to the reference model.

Results

N-related variables improved defoliation models. Higher N_thr deposition led to higher F₂₅ for beech and Norway spruce, while the effect was opposite for Scots pine. Higher foliar N ratios led to higher F₂₅ for all species.

Conclusion

N_thr deposition, damage to foliage, foliar N/P, N/Ca, N/Mg, N/K, top-soil pH, C/N and exchangeable base cation resulted important factors (although with possible diverse effect) in explaining the variability of F₂₅ among plots.

1 Introduction

The human-induced alteration of the nitrogen (N) cycle has caused a three to five times increase in atmospheric inorganic N deposition (largely nitrate and ammonium) during the last century, with a further increase by a factor of 2.5 projected at a global level by the end of this century (Lamarque et al. 2005). N availability controls the productivity of many forest ecosystems, and, under N limitation, increased N deposition stimulates growth and enhances carbon (C) sequestration. On the other hand, when N is not a limiting factor, further N addition by deposition may not stimulate growth because other nutrients become limiting (e.g. Erisman et al. 2011). In both cases, however, detrimental effects of increased N deposition on tree health have been reported due to nutrient imbalance, increased sensitivity to pests/pathogens (e.g. Roelofs et al. 1985), storms (e.g. Braun et al. 2003) and extreme weather conditions (e.g. Dobbertin 2005).

The response of forest ecosystems to N deposition can be evaluated by controlled experiments with N addition or by observational studies, e.g. long-term, large-scale monitoring. The experimental approach permits the identification of the effect of N inputs independently of other factors. Such an approach, however, is quite expensive and impractical with mature forests for which it has been generally restricted to a few sites or a limited period of time. Further, N treatments may not reflect the effect of chronic N deposition in terms of frequency and amount of N added. The analysis of observational data from long-term, large-scale concurrent monitoring of N deposition, tree nutrition, health and growth as well as important covariates such as soil and meteorological variables can be an alternative and complementary way of investigating the effects of N deposition on forests (Ferretti et al. 2014a). The potential for causal inference is lower in comparison with manipulative experiments, but observational approaches have the advantage of maintaining a high degree of realism with regard to environmental conditions and N deposition regimes. In this respect, the data collected at the intensive monitoring plots of the International Cooperative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests) are an important asset for two main reasons (Ferretti and Fischer 2013): Plots cover a wide range of ecological conditions (including a marked N deposition gradient), and—despite some inconsistencies that may always happen in large-scale programs—data are collected according to comparable methods (Ferretti and Fischer 2013). In these plots, monitoring of tree condition is based on the assessment of “defoliation”, a popular (and incorrect) term adopted to mean the reduction of foliage density on trees as compared to a reference standard. Despite several criticisms (e.g. Innes et al. 1993), the frequency of trees with defoliation >25 % (hereafter referred to as F₂₅, Ferretti and Fischer 2013) is largely used to report forest condition data and adopted as a Sustainable Forest Management (SFM) indicator under Criterion 2, Forest Health and Vitality (Forest Europe 2011). A number of factors, including stand age, meteorology, pests, and diseases were reported to explain most of the variation in defoliation (e.g. de Vries et al. 2014). The objective of this paper is to test the null hypothesis (H₀) of no effect of throughfall N deposition (N_thr) and related variables (soil and foliar nutrients) in defoliation models. To do this, we evaluated the possible, additive influence of N-related predictors and damage to foliage in explaining the spatial variability of F₂₅ once other important drivers (i.e. geography, stand, local climate and management) are taken into account.

2 Material and methods

2.1 Monitoring sites and quality assurance

Data were collected at the ICP Forests intensive monitoring (level II) permanent plots. Each plot consists of about 2500 m² in which several monitoring activities are carried out (Ferretti and Fischer 2013). To minimise the effect of random annual variations and anomalies, we concentrated on mean values over a defined time window and considered plots for which deposition, meteorological and defoliation data were available for all the 3 years 2007–2009 (see below).

A number of quality assurance and quality control (QA/QC) procedures are routinely carried out in order to improve accuracy and comparability over time and space (Ferretti and Fischer 2013).

Additional QA/QC checks were performed for this study, and only plots for which data passed the QA/QC ring tests for total N in deposition carried out in 2005 and 2009 were considered.

From this reduced data set, we considered plots with three main tree species (MTS) for which at least 15 plots were retained, namely beech (Fagus sylvatica L., n = 20), Norway spruce [Picea abies (L.) Karst., n = 33] and Scots pine plots (Pinus sylvestris L., n = 18). The number of trees evaluated ranged from 18 to 304 per plot and year. All in all, 71 plots were selected in nine countries covering a large latitudinal gradient across Europe (Fig. 1).

Fig. 1

Geographical distribution of the sampling plots

2.2 Study concept

The study concentrated on factors explaining the spatial variability of mean F₂₅ and was based on mean values per plot over a 3-year time window (2007–2009). Therefore, we did not consider time trends and were not concerned with time dependency. Within such a context, a stepwise procedure was adopted (Ferretti et al. 2014b): First, a reference model was considered, with geographical (latitude, longitude, elevation), stand (tree density) and meteorological variables (annual precipitation and per cent of precipitation occurring in summer) as predictors of F₂₅ (dependent variable). Secondly, individual sets of predictors were added stepwise to the reference model to evaluate their additional role in explaining F₂₅ (Table 1).

Table 1 Plotwise aggregated response and predictor variables selected for the PLS regression models

Since the number of plots with available data was not the same for all variables, we used three complete, internally consistent, aligned designs, each of them with a different number of plots, i.e. n = 71 for testing the role of N_thr deposition and foliar damage; n = 65 for testing the role of N-foliar ratios; and n = 56 for testing the role of soil variables. For each design, we compare the reference model (see above) against other models including additional variables. In this way, we were able to compare models based always on the same plots.

2.3 Response and predictor variables

The frequency of trees with defoliation higher than 25 % (F₂₅) for the MTS in each plot was adopted as response variable. Defoliation is the reduction of foliage density on a tree as compared to a reference standard, and it is visually assessed by trained observers according to semi-quantitative classes: 0 %, 5 % (>0–5 %), 10 % (>5–10 %), …, 99 % (>95–100 %, when a tree is still alive) and 100 % (dead tree) at 5 % steps (see details in Ferretti and Fischer 2013). As any other ecological data based on visual observation, defoliation assessment is prone to observer error (Innes et al. 1993), and QA/QC activities have been routinely implemented to control such an error (e.g. Ferretti and Fischer 2013). It is worth noting, however, that, since our study is aimed at identifying the additional role of N-related variables in a step-wise procedure that compares different models over the same plots, differences in defoliation assessment that could exist among surveyors are likely to have a limited impact on the overall results.

Predictors were chosen on the basis of their availability and include geographical and meteorological variables (in the reference model), reported foliar damage due to reasons other than air pollution (e.g. biotic factors such as pests, or abiotic factors such as storms) (hereafter referred as BADoAP), annual total N_thr deposition and a set of variables known to be directly or indirectly related to N deposition, namely foliar N ratios (N/Ca, N/P, N/K, N/Mg) and soil pH, exchangeable base cations (BCE) and C/N ratio in the mineral top soil (Ferretti et al. 2014b; Cools and De Vos 2011; Cools et al. 2014) (Table 1). For both predictors and response variables, assessment and measurement methods were those reported by Ferretti and Fischer (2013). BADoAP on foliage were evaluated on an ordinal scale for each tree and averaged for each MTS and plot. Throughfall deposition was sampled by 3–10 gutters (France and Slovenia) or 10–50 funnels located on a systematic basis (other countries). For most compounds, throughfall deposition is usually larger than open field deposition because the former also includes the dry deposition of compounds collected by the tree canopy. However, in some of the plots used in this study, throughfall deposition of N compounds (mainly ammonium) was lower than open field deposition because of the active uptake by the foliage (e.g. Ferretti et al. 2014b). Thus, we decided to use N_thr deposition of total N as a predictor because it represents the total amount of N carried to the forest soil through atmospheric deposition and canopy processes. Total N_thr includes both inorganic and organic N, which has different ecological effects; however, in most of Europe except for northern Fennoscandia, the contribution of organic N to total N_thr is relatively low.

Precipitation amounts were measured in a clearing close to the permanent plot. Chemical analyses of foliar tissues were available for 65 plots, while soil analyses were available for 56 plots, and plot average of the topmost mineral layers was considered. Further details on methods for sampling and analysis are reported by Ferretti and Fisher (2013; see also http://icp-forests.net/page/icp-forests-manual).

2.4 Statistical methods

Correlations among variables were explored using Spearman regression. Partial least square (PLS) regression (Wold et al. 2001; Mehmood et al. 2012) was used for statistical modelling of defoliation (Ferretti et al. 2014a). PLS regression generalises and combines features from principal component analysis (PCA) and multiple linear regression (MLR). It is particularly useful when the need is to predict a set of dependent variables (Y) from a large set of independent variables (X). Among the data mining techniques, PLS regression is considered effective for providing information on the relative importance of predictors (Mehmood et al. 2012), but it does not provide a robust way of testing significance of the coefficients for the predictors. The variable importance in the projections (VIP) scores are used to rank the importance of the predictors, and—since the average of the squared VIP scores equals 1—variables with scores >1 were considered to be the most important (Ferretti et al. 2014a). PLS regression was carried out by the non-linear estimation by iterative partial least squares (NIPALS) algorithm (Geladi and Kowalski 1986) for different set of predictors and separately for each tree species, using the pls package (Mevik and Wehrens 2007). All variables were standardised to zero mean and unit variance, and all analyses were performed in the R statistical environment (R Core Team 2013).

3 Results

3.1 Descriptive statistics and univariate relationships

Individual 2007–2009 plot average F₂₅ ranged from 0 to 43.3 % in beech, from 0 to 72.0 % for spruce and from 0 to 32.7 % for pine (Table 1). Annual throughfall N deposition ranged from 1.5 to 32.3 kg ha⁻¹ and showed significant latitudinal and longitudinal (p < 0.01) gradients (Table 2), with the lowest values in northern Finland and maximum values in Germany. Due to the different geographical distribution of tree species, however, the latitudinal gradient was negative for spruce and pine plots (N deposition decreases from central Europe to the Nordic countries) and positive for beech (N deposition decreases from central Europe to southern Italy and France). This has an effect on the N_thr deposition gradient experienced by the different species (see Table 1). The more southerly spatial distribution of beech plots, including calcareous areas in Italy and Switzerland (Fig. 1), is reflected in a larger range of pH and BCE values for this species (Table 1).

Table 2 Spearman rank order correlation coefficient between N deposition and the other variables considered in the study

Significant (p < 0.05) positive linear relationships were found between N_thr deposition and the foliar N/P, N/Ca, N/Mg and N/K ratios, across all data as well as for one or more of the studied tree species (Table 2). A negative relationship was detected between N_thr deposition and soil pH, measured in the topmost mineral soil, across all data and for pine and spruce plots. A positive relationship between soil BCE and N_thr deposition was evident only for pine plots.

3.2 Multivariate analysis

Factors contributing to explaining the variability of F₂₅ were identified using one reference plus five N-related PLS-regression models for each species (see “Material and methods”). As foliar and soil variables were correlated among themselves and with N deposition, they were added separately to the reference model. When interpreting the results, emphasis should be placed on the effect of the addition of a given group of variables (e.g. foliar N ratios), while the importance of each variable within the model can be biased by their mutual correlation.

The performance of each model is reported in Table 3. The results showed marked differences among species. The reference model was able to explain 4.6 % (beech), 8.1 % (Norway spruce) and 29 % (Scots pine) of the total variation in F₂₅ (Table 3). The most important predictors in the reference models were latitude (for Norway spruce and Scots pine), longitude (for Norway spruce), tree density (for all species), precipitation (for beech) and summer precipitation (for Scots pine) (Tables 4, 5 and 6). Model performance always improved when a set of N-related predictors was added. This was particularly obvious for beech (up to 41.7 % additional variance explained than in the reference model) and Scots pine (up to 40.6 %), and less marked for Norway spruce (up to 30.9 %) (Table 3).

Table 3 Model performance

Table 4 Variable importance in the projections (VIP) score for beech

Table 5 Variable importance in the projections (VIP) score for Norway spruce

Table 6 Variable importance in the projections (VIP) score for Scots pine

For beech, the contribution of the various sets of N-related predictors in improving the reference model was as follows: foliar > N_thrdep + BADoAP > soil > BADoAP > N_thrdep (Table 3). In particular, F₂₅ increases with increasing values of BADoAP, N_thrdep, the combination of N deposition and damage, foliar N ratios (N/Ca, N/K) and soil C/N, and decreases with increasing values soil pH and BCE (Table 4).

For Norway spruce, the contribution of the various sets of predictors in improving the reference model was as follows: foliar > N_thrdep + BADoAP > BADoAP ~ soil > N_thrdep (Table 3). In particular, F₂₅ increases with N deposition and damage (single and in combination), foliar N/P, N/Ca, N/Mg and N/K, and soil BCE, and decreases with soil C/N (Table 5).

For Scots pine, the improvement of the reference model was due to soil > N_thrdep + BADoAP ~ N_thrdep > foliar > BADoAP (Table 3). In particular, F₂₅ increases with foliar N/P and N/Mg, and decrease with N_thrdep and soil pH and C/N (Table 6).

Table 7 shows the number of models for which each variable obtained a VIP > 1. N_thrdep and N-related variable obtained VIP > 1 in almost all models where they were included. In the case of Scots pine, summer precipitation was as important as latitude (also important for Norway spruce) and tree density (all species).

Table 7 Summary of the PLS regression models

4 Discussion

High levels of inorganic N deposition on forest ecosystems due to anthropogenic activities have been reported for decades in Europe and North America. Most studies have shown how N deposition affects soil and foliar chemistry, tree growth, carbon sequestration and plant species richness and diversity (e.g. de Vries et al. 2014), but only a few studies have discussed the effects of elevated N deposition on forest health in terms of defoliation (e.g. Lomsky et al. 2012; Klap et al. 2000 and references therein). Our results suggest that—to different extents—variables related to N deposition can all contribute to explain the frequency of trees with defoliation >25 %. Positive correlation between N deposition and F₂₅ was estimated for beech (most of the beech plots have a suboptimal foliar N status), while a negative one was observed for pine, whose foliar N status was low for half of the plots located in Nordic countries where also N deposition is low. These apparently contradictory effects of N deposition on defoliation may be explained when considering the different ranges across the investigated N deposition and saturation gradient covered by the investigated tree species.

In fact, the degree of N saturation (Aber et al. 1998) is an underlying factor/process influencing forest health. For all tree species considered, F₂₅ increases with increasing foliar N ratios, which are good indicators of N-induced nutritional imbalances. In particular, the foliar N/Ca and N/K ratios were the most influential in the case of beech; foliar N/Mg, N/Ca and N/P for spruce, and foliar N/P and N/Mg for pine, indicating that crown condition in all these tree species is related to imbalance in foliar nutrients (e.g. Veresoglou et al. 2014). N-induced nutrient imbalances can be either caused by direct foliar N uptake (e.g. Fleischer et al. 2013) or mediated through soil processes (e.g. Moore and Houle 2013). In particular, N deposition effects on soil may include soil acidification, due to nitrification and leaching of nitrate (with possible release of toxic aluminum compounds) and leaching of base cations (for a synthesis, see de Vries et al. 2014). In our data set, soil pH was significantly and negatively correlated with N_thr deposition (p < 0.01), which can probably be explained by similar, though unrelated, geographical gradients. In the Southern Mediterranean part of Europe, soil types tend to be more calcareous and have higher soil pH values compared to Central or Northern Europe. Simultaneously, N deposition values tend to be lower in Southern and Northern Europe.

The effect of soil pH on F₂₅ was indeed important for beech plots (where the soil pH ranges between 3.10 and 6.95), but not for Norway spruce and Scots pine (soil pH range, 2.9–4.7). Similarly, the fact that latitude had a positive effect on beech F₂₅ and a negative one on spruce and pine F₂₅ is probably also due to a “range effect”. The effect of latitude in the reference model and of soil pH in the model including soil variables probably accounts for part of the effect of the deposition and N saturation gradient in Europe.

The impact of leaching of base cations on tree nutrition, coupled to an increased availability of N, leads to positive relationships between N deposition and foliar N ratios, as shown in many observational and fertilisation studies (e.g. de Vries et al. 2014). In our data, almost all foliar ratios increase with increasing N deposition. In particular, the significant relationship between N_thr deposition and foliar N/K for all tree species supports the hypothesis of K leached from leaf tissue in exchange with ammonium ions (Moreno et al. 2001) as K does not easily leach from soil. More generally, N deposition may alter both terms of the foliar N ratios. On the one hand, N deposition improves N nutrition, and on the other hand, it can contribute to deteriorated nutrition of other nutrients by acting as a fertiliser promoting tree growth and therefore increasing nutrient demand by trees (Lukac et al. 2010), by reducing nutrient uptake capacity of trees due to reduced carbon allocation to roots and detrimental effects on mycorrhiza (Kjøller et al. 2012) and by reducing nutrient availability in the soil due to acidification.

In summary, N deposition and N-deposition-related variables at soil and foliar level improved substantially (3.9–41.7 %, according to the species and the predictor) explanatory models of the frequency of defoliated trees in beech, Norway spruce and Scots pine across Europe.