Seasonal variations of throughfall chemistry in pure and mixed stands of Oriental beech (Fagus orientalis Lipsky) in Hyrcanian forests (Iran)

Abstract

Key message

Throughfall nutrient fluxes were generally higher in a mixed stand than in a pure stand of Oriental beech. Throughfall fluxes were higher than bulk precipitation fluxes except for nitrogen and indicate higher canopy uptake of this element in the pure stand compared to the mixed stand.

Context

Oriental beech is an economically relevant tree species for Iran and adjacent countries. Yet little is known about nutrient cycling in Oriental beech stands and the influence of the degree of mixture with other species.

Aim s

We assessed the effect of seasons on nutrient fluxes in precipitation and whether throughfall chemistry differed between pure and mixed stands.

Methods

Bulk precipitation in the open field and throughfall were sampled during one whole month within each season from August 2013 to May 2014 in a pure (81 % of beech trees on average) and a mixed stand (57 % of beech trees) of Oriental beech. Samples were analysed for pH, nitrate (NO₃ ⁻), ammonium (NH₄ ⁺), phosphorus (P), calcium (Ca²⁺), magnesium (Mg²⁺) and potassium (K⁺).

Results

Nutrient concentrations were generally higher in the growing season than in the dormant season, both in bulk precipitation and in throughfall. Nutrient fluxes were higher in fall and these peaks coincided with higher amounts of precipitation. The concentrations and fluxes of NH₄ ⁺, NO₃ ⁻, P, Ca²⁺, Mg²⁺ and K⁺ were generally higher in the mixed stand than in the pure stand in all seasons. Compared to the open field, throughfall fluxes were usually higher, except for NO₃ ⁻ and NH₄ ⁺, indicating direct canopy uptake of nitrogen.

Conclusion

Canopy composition in Oriental beech stands (owing to differences in foliage chemistry) and seasons (owing to differences in precipitation regime and phenological stages) have a significant effect on throughfall nutrient fluxes.

1 Introduction

Forest nutrient cycles are tightly connected to the hydrological cycle, because rainfall and soil water act as the main transportation and solution media (Bhat et al. 2011; Bruijnzeel 1989). Forest canopy has a major impact on soil nutrient availability through the amount and chemical composition of litter and throughfall (Carnol and Bazgir 2013; De Schrijver et al. 2007; Hansen et al. 2009; Prescott 2002). Throughfall is that part of incident precipitation which passes through the forest canopy, either directly in gaps or after interacting with foliage and branches (Parker 1983).

In forest ecosystems, before rain water reaches the forest floor as throughfall or stemflow (i.e. the throughfall fraction that flows down tree stems), the chemical content of rain water is strongly affected by exchanges within tree canopies (Herrmann et al. 2006). Rainwater is thereby influenced by three main processes (Chiwa et al. 2004): (1) the washing of deposits accumulated on the canopy between precipitation events, (2) the leaching of material from internal plant tissues and (3) the uptake by foliage of solutes. Tree species with different foliage properties exhibit different patterns of throughfall water and ion fluxes (Beier et al. 1993; Hansen 1995; Raat et al. 2002; Whelan et al. 1998). Furthermore, as dry deposition is also a function of foliage properties (Lovett 1994), it may cause species-related differences (Hojjati et al. 2009; Rothe et al. 2002).

As opposed to litterfall, throughfall and stemflow add nutrients directly to the available nutrient pool without the interference of any process of decomposition (Abbasian et al. 2014). Atmospheric deposition is difficult to measure directly. However, the nutrient content of throughfall is often used as a measure of deposition for ions that show minimal foliage uptake or leaching (Rothe and Binkley 2001).

In recent years, mixed stands have received growing attention due to their greater biodiversity, ecological stability and their greater productivity in comparison to pure stands (Pretzsch et al. 2010). Few authors, however, have assessed the impact of tree species composition on atmospheric deposition by studying throughfall fluxes both in pure and mixed stands. Berger et al. (2008) compared stands of pure beech, pure spruce and a mixture of beech and spruce and found that nitrogen and sulphur throughfall fluxes were lowest in the beech stand and highest in the spruce stand, while the mixed stand showed intermediate values. In this same study, the leaching of base cations tended to be lowest in the beech stand, intermediate in the mixed stand and highest in the spruce stand. Talkner et al. (2010) found that the throughfall and stemflow fluxes of K⁺, Ca²⁺, Mg²⁺ and PO₄ ³⁻ were higher in mixed deciduous species plots than in pure beech plots due to higher canopy leaching rates in the mixed species plots.

The temperate deciduous forests of northern Iran, known as the Hyrcanian forests (also called Caspian forests), cover about two million hectares in a narrow strip over 800 km long and 20–70 km wide (Hosseini Ghaleh Bahmani et al. 2012) and encompass different forest types thanks to their 80 tree and shrub species (Sagheb Talebi et al. 2004). Oriental beech (Fagus orientalis Lipsky) is an economically relevant tree species for Iran and adjacent countries. Beech stands in northern Iran are of special concern since Oriental beech is one of the most valuable and important timber species, covering 17.5 % of the surface area and accounting for 30 % of the volume of the industrial forests in Iran (Sagheb Talebi and Eslami 2008).

Forest productivity is limited by the supply of one or more nutrients in almost all forests, and forest nutrition management is a key factor in the management of commercial forests (Rouhi-Moghaddam et al. 2008; Stape et al. 2006). While there are many studies comparing the mineral nutrition and nutrient cycling in beech stands (Berger et al. 2009a; Staelens et al. 2007), most of them have focused on European beech (Fagus sylvatica L.) and information about throughfall chemistry below Oriental beech canopy is still rare (Abbasian et al. 2014; Moslehi et al. 2011). To date, most studies on Oriental beech stands have been restricted to the quantification of interception loss (Ahmadi et al. 2009; Mattaji et al. 2012; Rahmani et al. 2011). Because of the significant lack of data on throughfall chemistry under Oriental beech canopy, our focus in the present study is throughfall chemistry in beech stands in the Caspian forest. We hypothesized that the chemical composition of throughfall differs in the growing season (spring and summer) compared to the dormant season (late fall and winter) and in pure compared to mixed stands of Oriental beech. Our study addresses the impact of tree composition on throughfall chemistry by asking the following questions:

• 1. Does tree composition in mixed and pure stands affect atmospheric input?

• 2. Do seasons affect nutrient fluxes in beech stands?

In addition to these specific research questions, owing to the lack of deposition data for this region and these types of stands, we attempted to estimate annual rainfall and throughfall fluxes of the main nutrients for the period July 2013 to June 2014, based on the four monthly fluxes that we measured.

2 Material and methods

2.1 Site description

This study was conducted from August 2013 to May 2014 in intact natural stands of F. orientalis, within the Gorazbon district at the Kheyrud Forest Research Station of the University of Tehran, located approximately 7 km east of Nowshahr, Mazandaran province, northern Iran (36°33′N, 50°33′E) (Fig. 1). The Gorazbon district covers about 1000 ha and its elevation ranges from 550 to 1380 m a.s.l. Ground slopes range from 0 to 115 %, and the average slope in the whole district is about 15 %. This area is dominated by natural temperate forests containing native mixed deciduous tree species such as F. orientalis Lipsky, Carpinus betulus L. and Acer velutinum Boiss. (Jourgholami et al. 2013). The mean annual precipitation is 1146 mm, with the highest amount in October (200 mm) and the lowest amount in July (49 mm). Two sample parcels (Fig. 1) were selected within the Gorazbon district, one covered by an almost pure Oriental beech stand (with beech representing 70 % of the standing volume), the other covered by a mixed stand (52 % of standing volume; Table 1).

Fig. 1

The experimental area located at the Forest Research Station of the University of Tehran, Iran

Table 1 Stand characteristics in the two forest parcels in which the study plots were selected

Both parcels are in the same elevation class (1150–1350 m a.s.l.) and lie on similar parent material (limestone). Mean tree density and volume were 225 trees ha⁻¹ and 463 m³ ha⁻¹ standing volume in the mixed stand and 266 trees ha⁻¹ and 500 m³ ha⁻¹ standing volume in the pure stand (Table 1). In each of these two parcels, 30 rectangular plots, 500 m² each, homogeneous with respect to slope and aspect, were delimited. In addition, four open areas (two in each parcel) were selected for the sampling of gross rainfall (Fig. 2). We studied the seasonal variation of precipitation chemistry in the pure and mixed stands by intensive sampling in mid-spring, mid-summer, mid-fall and mid-winter at the plots.

Fig. 2

Location of the open areas and plots in pure and mixed stands of Oriental beech trees

2.2 Bulk precipitation and throughfall measurement

Precipitation in the open field (bulk precipitation, BP) and throughfall (TF) were sampled during one whole month within each season: in summer (August 2013), fall (November 2013), winter (February 2014) and spring (May 2014). One manual plastic collector, 9 cm in diameter and 25 cm in height, was installed in each of the four open areas for bulk precipitation sampling. For the sampling of throughfall, six collectors were installed in each plot. A nylon mesh was placed at the opening of the collectors to prevent contamination by coarse debris. The collectors were emptied twice, in the middle and at the end of the sampling months. In the middle of the sampling month, after sample collection, the collectors were rinsed with distilled water before being replaced. The volumes collected were measured using a graduated cylinder and the contents of the six collectors were pooled in the field to one composite sample per plot, resulting in 2 × 30 samples for each collection period. The samples were then taken to the laboratory and stored at 4 °C. One exception to this general procedure concerned the first 2 weeks of February, during which precipitation fell in the form of snow. After the 2-week exposure, the collectors were replaced by clean ones and pooling took place after melting in the laboratory.

The samples of the two collection periods within a month were pooled before analysis to obtain monthly samples. The pH of the rainfall was measured with a pH meter. The nitrate (NO₃ ⁻) and ammonium (NH₄ ⁺) concentrations were measured using the Kjeldahl method. The NH₄ ⁺ fraction in the sample was first determined by distilling the sample with sodium hydroxide. The amount of ammonia into which NH₄ ⁺ was converted was then determined by back titration. The NO₃ ⁻ fraction was determined using the remaining solution, into which Devarda’s alloy was added to reduce NO₃ ⁻ into NH₄ ⁺, which in turn was analysed as described above. Phosphorus (P) was measured with a UV/V spectrophotometer (SQ-2800, US). Concentrations of potassium (K⁺), calcium (Ca²⁺) and magnesium (Mg²⁺) were determined using a flame photometer (Abbasian et al. 2014).

We calculated monthly fluxes of NO₃ ⁻, NH₄ ⁺, P, Ca²⁺, Mg²⁺ and K⁺ in bulk precipitation and throughfall for August 2013, November 2013, February 2014 and May 2014 by multiplying concentrations by the precipitation volume. In order to estimate fluxes on an annual scale, we further used the monthly amounts of precipitation measured at the Nowshahr meteorological station from July 2013 to June 2014 as an estimate of bulk precipitation for the months for which we had no data. To estimate throughfall volumes for the missing months, we calculated an interception factor for each season as the ratio of measured throughfall to bulk precipitation volumes and applied this factor to the bulk precipitation estimates. Within each season, the measured monthly nutrient concentrations were applied to the other 2 months for which no values were available. Measured and estimated monthly fluxes were then totalled to obtain estimates of the annual nutrient fluxes in bulk precipitation and throughfall for the period July 2013 to June 2014.

As a second estimate, given that our measurements were well distributed over the year, we extrapolated the sum of measured fluxes to annual values by simply dividing this sum by the number of measurement days and by multiplying by 365 days.

2.3 Statistical analysis

The normality of the variables was checked using the Kolmogorov-Smirnov test, and Levene’s test was used to examine the equality of variances. The effect of stand type, season and the interactions between stand type and season were evaluated by repeated measures ANOVA. Differences between stand types and between seasons were tested using the Duncan test. For all statistical analyses, the SPSS software (version 22) was used.

3 Results

3.1 Spatial variability of the open areas

In a preliminary analysis, in order to check the comparability of the forest parcels, we assessed the spatial variability of our open areas. All four open areas showed very similar precipitation volume and chemistry. There were no significant differences between the two sets of open areas located in each of the two forest parcels (not shown).

3.2 Quantity of bulk precipitation and throughfall

During the period from August 2013 until May 2014, the amount of bulk precipitation was significantly higher than the amount of throughfall at the forest stands (P < 0.05; Fig. 3). Bulk precipitation amounts ranged from 56.7 in summer to 119.2 mm month⁻¹ in fall. Monthly throughfall amounts in the pure and mixed stands were higher in the dormant season compared to the growing season due to higher precipitation in the fall and winter months. They ranged from 44.4 to 84.8 in the pure stand and from 45.4 to 85.3 mm month⁻¹ in the mixed stand. Interception loss was slightly higher in the pure stand than in the mixed stand. The difference was significant (P < 0.05) in spring, summer and winter.

Fig. 3

Mean (±SD) monthly precipitation volumes and nutrient concentrations in rainfall (measured in open areas) and throughfall in pure and mixed beech stands. Different lowercase letters indicate significant (P < 0.05) differences among stands and different capital letters indicate significant (P < 0.05) differences among seasons

3.3 Nutrient concentrations and fluxes

The concentrations of nutrients were generally higher in the growing season than in the dormant season, both in bulk precipitation and in throughfall (Fig. 3). In the growing season, under mixed canopy, the concentration of N in the form of NO₃ ⁻ (NO₃ ⁻-N) was significantly higher in throughfall than in bulk precipitation, while the opposite was true for the pure stand in spring. In the dormant season, bulk precipitation showed similar or higher NO₃ ⁻-N concentrations than throughfall. The concentration of N in the form of NH₄ ⁺ (NH₄ ⁺-N) tended to be higher than that of NO₃ ⁻-N. In the mixed stand, the NH₄ ⁺-N concentration was significantly higher in throughfall than in bulk precipitation in spring, summer and winter. In the pure stand, the only significant difference in comparison to bulk precipitation was observed in spring, with lower NH₄ ⁺-N in throughfall (as also observed for NO₃ ⁻-N). With few exceptions, NO₃ ⁻-N and NH₄ ⁺-N concentrations were significantly higher in the mixed stand than in the pure stand (TF mixed ≥ TF pure). P concentration increased after passing through the canopy. The P concentrations under the pure stand, under the mixed stand and in the bulk precipitation were significantly different (TF mixed > TF pure > BP). The concentrations of Ca²⁺, Mg²⁺ and K⁺ followed a similar pattern of TF mixed ≥ TF pure ≥ BP in all seasons. The pH value of throughfall under both stands was significantly higher than in bulk precipitation in the growing season and tended to decrease from the growing season to the dormant season. The pH of throughfall under mixed canopy was significantly higher than under pure canopy in the growing season.

Fluxes of nutrients were generally higher in fall than in other seasons, and these peaks coincided with higher precipitation amounts (Fig. 4). This was true for all nutrients in bulk precipitation (except K⁺ and Ca²⁺), as well as for P, Ca²⁺, Mg²⁺ and K⁺ in throughfall for both stands. The same pattern was apparent for NO₃ ⁻-N and NH₄ ⁺-N; however, the differences were not always significant. The flux of NO₃ ⁻-N tended to be higher in bulk precipitation than in throughfall in each season except in the mixed stand in summer. NH₄ ⁺-N fluxes tended to be higher in bulk precipitation compared to throughfall in both types of stands in the dormant season and compared to throughfall in the pure stand in spring. NO₃ ⁻-N and NH₄ ⁺-N fluxes were generally higher under mixed canopy than under pure canopy (TF mixed ≥ TF pure). Monthly fluxes of P, Mg²⁺, Ca²⁺ and K⁺ mostly followed the order: TF mixed ≥ TF pure ≥ BP. This pattern was observed in all seasons for P and K⁺ fluxes.

Fig. 4

Mean (±SD) monthly nutrient fluxes in rainfall (measured in open areas) and throughfall in pure and mixed beech stands. Different lowercase letters indicate significant (P < 0.05) differences among stands and different capital letters indicate significant (P < 0.05) differences among seasons

Our estimation of annual fluxes (Table 2) showed that although the annual throughfall amount was similar under both canopy types, throughfall fluxes of all nutrients were higher in the mixed stand than in the pure stand by approximately 5 to 40 % depending on the method and the nutrient. For example, the estimated N fluxes varied between <10 kg^-1 ha^-1 year^-1 in the pure stand and 11 to 12 kg^-1 ha^-1 year^-1 in the mixed stand, while P fluxes lay between <2 kg^-1 ha^-1 year^-1 in the pure stand and >2.5 kg^-1 ha^-1 year^-1 in the mixed stand. Compared to bulk precipitation fluxes, throughfall fluxes were higher except for NO₃ ⁻-N and NH₄ ⁺-N.

Table 2 Approximate annual precipitation and fluxes of nutrients in beech stands and in the open area for the period July 2013 to June 2014. Method 1 uses precipitation measured by the nearby Nowshahr meteorological station; method 2 extrapolates the measured values by weighting by the number of days (see text for a more detailed description)

4 Discussion

Throughfall chemistry and nutrient fluxes were investigated from August 2013 to May 2014 in pure and mixed stands of Oriental beech. To date, while considerable research has been done comparing throughfall chemistry in European beech stands, to our knowledge, there is no published study on the effect of season and degree of mixture with other species on throughfall chemistry under Oriental beech canopies.

4.1 General effect of the tree canopy on precipitation chemistry and seasonal variation

Our results regarding the general effect of Oriental beech canopy on precipitation chemistry and its seasonal variation are consistent with those widely found in the literature for other deciduous tree species. In our study, throughfall water was less acidic than bulk precipitation, especially in the growing season, which is due to the capacity of deciduous canopies to neutralize precipitation (Devlaeminck et al. 2005). The throughfall fluxes of almost all the analysed nutrients differed significantly from the bulk precipitation fluxes, both in the dormant and the growing season. These differences reflect the wash-off of dry deposited particles and the exchange with leaf surfaces (absorption and leaching; Parker 1983). The lower NO₃ ⁻ and NH₄ ⁺ concentrations and fluxes observed under canopy compared to bulk precipitation, including in the dormant season, show that these two forms of inorganic N were absorbed by foliage, also in the senescent stage, and by branches. Adriaenssens et al. (2011) likewise found a significant uptake of nitrogen, mostly of NH₄ ⁺ but also of NO₃ ⁻, by the senescent leaves of European beech. Nitrogen uptake at leaf senescence was assumed by these authors to be linked to N retranslocation by trees (see also Houle et al. 1999). The lower N fluxes in throughfall compared to bulk precipitation also indicate that dry deposition of inorganic N compounds on foliage, in the form of gases or aerosols, is probably quite low in our study site. In many European temperate broadleaved forests, throughfall fluxes exceed bulk precipitation fluxes due to the wash-off of dry-deposited compounds (e.g., Kristensen et al. 2004; Thimonier et al. 2005), even though N is most likely taken up by the canopy as well (Harrison et al. 2000).

The strongly enhanced throughfall K⁺ concentrations and fluxes that we observed at our beech stands is also well documented in the literature for tree canopies (Stachurski and Zimka 2002). Potassium is known to be readily leached from foliage during rainwater passage (Samecka-Cymerman et al. 2010). Together with Mg²⁺ and Ca²⁺, K⁺ is known to be leached especially in spring during leaf sprouting (Herrmann et al. 2006) and during leaf senescence. Phosphorus, mostly in the form of phosphate, is likewise known to be leached from the canopy (Potter et al. 1991). Accordingly, we observed higher throughfall concentrations and fluxes of this element at our plots. Our study confirms the significant effect of seasons on the nutrient input to forest soils via throughfall, which is due both to variations in bulk precipitation and to stand phenology.

4.2 Effect of stand type (pure or mixed)

Although our pure and mixed stands are located in two distinct forest parcels, the four open areas (two open areas in each parcel) showed highly similar rainfall volumes and chemistry, which indicates that both types of stands are subjected to the same deposition climate. Even though this was expected from the careful selection of sites with similar elevation ranges and topographies, these results confirm the comparability of the two types of stand.

In the present study, we did not measure the leaf area index of our two types of stands, which is an important stand parameter that influences the correlation between canopy cover and deposition processes (Staelens et al. 2006). However, tree density and tree volume data are good indicators of stand structure. Tree density and standing volume were slightly higher in the pure stand compared to the mixed stand at the parcel level, which is reflected by slightly higher interception losses in the pure stand. However, the differences in interception are very small (see Table 2 on the annual scale). This suggests that the differences in throughfall chemistry that we observed between pure and mixed stands are to be ascribed to differences in foliage properties (translating into, for example, different leaching or absorption rates of nutrients) rather than to differences in throughfall volumes.

The concentrations and fluxes of all nutrients were significantly higher in the mixed stands compared to the pure stands across all seasons, with few exceptions. Even though throughfall collectors were installed below the canopies of beech trees, the crown of beech trees was intertwined with that of other broadleaved species in the mixed stand. The differences between pure and mixed stands can thus be explained by differences in the foliage properties of beech trees themselves or by the influence of the other tree species.

To our knowledge, there are no studies that compare the foliage chemistry of Oriental beech to that of any of the other broadleaved species in mixture with beech as found at our sites. However, European beech trees are known to have lower foliage nutrient concentrations than European hornbeam (C. betulus L.) (Krauss and Heinsdorf, 2005). Furthermore, foliar analyses performed on the very beech trees under which throughfall was sampled showed that foliage concentrations of N, Ca and Mg were higher in the mixed stand than in the pure stand (unpublished results). The higher foliar concentrations in the mixed stand were related to higher measured nutrient concentrations in the soil of our plots in the mixed stand compared with the pure stand, P being the exception (higher P concentrations in both beech foliage and soil in the pure stand; unpublished results). The higher throughfall nutrient fluxes observed in our mixed stand are thus likely to be related to the higher foliar concentrations in the beech trees of the mixed stand and in the other species in mixture, as well as the higher availability of nutrients in the soils under mixed canopy. Our findings are thus similar to those of Talkner et al. (2010). These authors also observed higher leaching fluxes of K⁺, Mg²⁺ and Ca²⁺ under mixed beech canopy compared to pure stands of European beech, which they related to differences in foliage and soil properties. The lower N throughfall fluxes in our pure stand, on the other hand, may be ascribed to a higher uptake of this element directly from precipitation, in relation to a lower N availability as indicated by lower foliar and soil N concentrations, and a higher carbon to nitrogen ratio (C:N) in the upper soil layer in the pure stand. For P, however, the higher fluxes under mixed canopy are inconsistent with the lower P concentrations in the beech leaves and the soil of the mixed stand. Higher P concentration in throughfall could be related to the overlapping of canopies of different trees species in the mixed stand.

4.3 Estimation of annual fluxes

Using our precipitation measurements in the open areas from four separate months in the year, complemented by the precipitation measured by the Nowshahr meteorological station, we estimated an annual precipitation of 1114 mm year^-1 for the period July 2013 to June 2014. The simple extrapolation method (weighting by the number of days) resulted in a value of 1009 mm year^-1. These two estimates are close to the long-term mean of 1146 mm year^-1. Our measurement period seems therefore to be representative of a normal year in terms of precipitation, the quantity of which has been shown to be a factor in the variation in annual bulk and throughfall fluxes (e.g., Thimonier et al. 2005).

The simple extrapolation method by weighting by the number of days resulted in somewhat lower but on the whole similar fluxes (Table 2). The annual fluxes of Mg²⁺, NH₄ ⁺-N and NO₃ ⁻-N estimated in this study were within the range of typical values for temperate European beech forests (e.g., Berger et al. 2009b; Thimonier et al. 2005). Even though our estimated throughfall flux of Ca²⁺ was rather low in comparison to other findings in the literature, Berger et al. (2009b) found similar low values. The estimated throughfall flux of K⁺, on the other hand, was rather high in comparison to values in Parker (1983), Thimonier et al. (2005) and Hojjati et al. (2009). We cannot exclude the possibility that we over- or underestimated throughfall fluxes for some of the nutrients. However, the fact that our two methods resulted in quite consistent estimates indicates that these values are robust.

With about 12 to 13 kg ha^-1 year^-1 for bulk precipitation fluxes, 9 to 12 kg ha^-1 year^-1 in the forest stands, the values are within the range of empirical critical loads for nitrogen as a nutrient (10-20 kg ha^-1 year^-1 for broadleaved stands; Bobbink and Hettelingh 2011). Below these values, no negative effects on sensitive parts of the ecosystem are expected to occur. The fact that the canopy visibly takes up N might indicate that the ecosystem is still a sink for this element.

5 Conclusion

The results of this study show that canopy composition in Oriental beech stands has a great effect on nutrient fluxes via throughfall. The chemistry of throughfall in pure and mixed stands appears to be influenced by differences in the foliage chemistry of beech trees themselves and possibly that of accompanying tree species, as well as by differences in soil properties among the two types of stands. Characterizing the nutrient status and fluxes for the tree species that are in mixture with Oriental beech in the Gorazbon forest would be helpful to interpret the differences between mixed and pure stands. There is also a need to investigate deposition in these stands for a longer time period in order to better estimate annual fluxes of nutrients and their year to year variability.