*
 

iForest - Biogeosciences and Forestry

iForest - Biogeosciences and Forestry
*

Species interactions in pure and mixed-species stands of silver fir and European beech in Mediterranean mountains

iForest - Biogeosciences and Forestry, Volume 14, Issue 1, Pages 1-11 (2021)
doi: https://doi.org/10.3832/ifor3476-013
Published: Jan 02, 2021 - Copyright © 2021 SISEF

Research Articles

Interactions between tree species determine the dynamics of forest communities. Spatial and temporal changes in resource availability, variation in species composition and spatial distribution of trees may alter competitive interactions between species and, therefore, affect tree growth and forest productivity. In this study, we analyzed the intra and inter-specific interactions between European beech (Fagus sylvatica L.) and silver fir (Abies alba Mill.) in southern Italy (Molise and Calabria regions), and how these interactions affect basal area increments in mixed-species and pure stands. Results showed that intra-specific interactions have a negative effect on the basal area increment, both in pure and mixed-species stands of Molise and Calabria. Basal area increment was higher influenced by intra-specific interactions in pure stands than in mixed-species stands. Silver fir in Molise showed higher basal area increment in mixed-species stand, probably in relation with stand structure and space occupation that resulted in less competition between individual trees. European beech showed high values of intra-specific interactions in pure stands, likely related to the low self-tolerance of this species and to the spatial arrangement of trees, due to canopy closure. The absence of inter-specific interactions in mixed-species stands could be explained by the sub-dominant position of European beech, which may have limited the benefit derived from niche separation and complementarity for silver fir.

Abies alba Mill., Fagus sylvatica L., Intra-specific Interactions, Inter-specific Interactions, Stand Productivity, Tree Growth

  Introduction 

Interactions between trees in forest stands are often dynamics due to spatial and temporal changes in resource availability ([64], [12], [20], [21]). Competitive interactions may also increase due to changes in species composition and to different spatial patterning of trees (clustered or uniform - [37]), strongly influencing tree growth. Therefore, interactions between tree species are considered the main mechanisms responsible for the effects of niche complementarity and the most important drivers for the dynamics of forest communities.

Several studies have shown that higher tree species richness may lead to greater forest stand productivity ([67], [35]). Mixed-species stands showed higher degree of resistance and resilience to biotic and abiotic disturbances ([32], [60]), and wider range of ecosystem services in comparison to corresponding monocultures ([26]). Nevertheless, further experimental insights on the ecosystem functions of mixed-species stands and the ecological mechanisms that cause their over-yielding is required for implementing resource-efficient management strategies ([21]), particularly, in drought-prone environments ([29], [15]). Indeed, under limited resource availability, as the number of individuals in a tree population tends to fill the available growing space and the density continues to increase, competition for water (besides nutrients and light) increases as well ([63]).

Interactions between individual trees in a forest stand can be studied and identified through the use of competition indices ([16], [6], [7]). In mixed-species populations, in order to better understand the interactions that occur between tree species, a distinction between intra- and inter-specific competition indices is recommended ([18]). Several studies have found that in presence of a reduction of competition or of a facilitating effect ([18]) inter-specific interactions can advantage tree growth in mixed-species stands ([3]). Although in mixed-species stands the interactions are often less severe than in the corresponding monocultures ([13], [45]), the positive effects of neighborhood complementarity on tree growth may depend on several factors, such as stand density ([14]), site quality ([49]) and changes in growth conditions ([18]). It must be pointed out, however, that distinguishing between competition reduction and facilitation effect is almost impossible and, therefore, these two processes are often defined as complementarity ([23], [1]).

Competitive interactions between tree species, but also complementarity effects, may occur above- and/or below-ground. In forests with tree species that differ intrinsically in crown architecture and leaf phenology (e.g., deciduous vs. evergreen species, conifers vs. broadleaved trees), species-specific occupation of the growing space and partition of the intercepted radiation may play a major role in the use of complementary resources ([22], [53]). Moreover, in mixed-species stands, with shade-tolerant species and a heterogeneous spatial structure, the use of limiting resources (e.g., light) can be more efficient, since the trees in the dominated plane may take advantage of irradiance transmitted through the crown of the dominant trees, so as to survive and grow ([52], [38]). Further, also co-occurrence of species that differ in their root architecture may improve the uptake of nutrients and water ([11]). Below-ground competition does not only consider the interactions between dominant or/and sub-dominant trees with the subject tree, but also those trees whose roots occupy the root distribution area of the subject tree ([10]). For example, silver fir (Abies alba Mill.) develops a deeper root system than European beech (Fagus sylvatica L.) and adopts a strategy that allows the absorption of water during the growing season, namely in periods of high drought stress ([65]). Species-specific root uptake, but also capacity stomatal adjustment, may affect the hydraulic redistribution of soil water, particularly in drought-prone Mediterranean environments. Though unclear distinction between silver fir and European beech in the sensitivity of stomatal conductance to water potential has been reported ([33]), silver fir promptly closes the stomata to maintain leaf water potential during drought events (i.e., isohydric behavior - [39]). European beech, instead, continues to transpire even when soil water content diminishes (i.e., anisohydric behavior - [39]).

European beech has been often studied mixed with the comparatively drought-intolerant Norway spruce (Picea abies Karst. - [48]) or drought-tolerant Scots pine (Pinus sylvestris L. - [15]). However, mixed-species forests that include silver fir and European beech have recently earned attention and praise among forest researchers ([18], [10], [54]). Renewed interest for silver fir can be related to recent findings of enhanced productivity of silver fir in central Europe ([9]) and to the relatively high drought tolerance of this species, particularly in southern Italy and mountain environments ([2]). Silver fir has a geographical distribution similar to that of European beech, though restricted to mountain ranges, and grow well in the southernmost locations in the mountain ranges of Italy ([41]). Toward its southern distribution limit, silver fir can be competitive with European beech (moderately drought tolerant) and, in central Europe, a good replacement for the more drought-sensitive Norway spruce at higher elevations ([59]).

This study aimed to analyze intra- and inter-specific interactions between European beech and silver fir, and the influence these interactions may have on the growth of trees in mixed-species stands and corresponding pure stands. In order to better understand the interactions between these two major forest species, we used competition indices at the individual tree level. Individual competition indices reflect the local density of the competitors, which interact with the reference tree ([66]), giving the possibility to quantify the influence of neighboring individuals on the growth of the reference tree. We hypothesized that: (i) intra-specific competition is higher in pure stands than in mixed-species stands; and (ii) tree species mixing and stand structural complexity mitigate the negative effect of inter-specific interaction.

  Materials and methods 

Study sites

The study sites included two triplets located in southern Italy (Molise and Calabria regions). In particular, the following forest sites were sampled: (i) “Abeti Soprani” (Molise) and (ii) Sila National Park (“Sila Piccola”, Calabria - Fig. 1).

Fig. 1 - Location of the study sites in the southern Italy. Red triangle indicates silver fir in the pure stands (S); green circle indicates European beech in the pure stands (B); yellow square indicates the mixed-species stands (M).

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Each triplet is composed of three rectangular stands, including two pure stands of silver fir and European beech, and one mixed-species stand of these species. Stand size varied between 0.12 and 0.28 ha. Stands were selected according to the minimum number of individuals (≥ 50 trees of silver fir and/or European beech per stand) and on the elevation range (<250 m between the triplets). Furthermore, mixed and pure stands were identified on the basis of the basal area occupied by the tree species investigated inside the stand, according to the classification of Mina et al. ([43]). Stands were considered pure when the basal area of the studied tree species was at least 85%. Mixed-species stands, on the other hand, were classified as a mixture of the two species when, both species represented 85% or more of the total basal area and each of the two species contributed at least 20% of the basal area (Tab. 1). The horizontal structure was different in the two regions. Stands in Molise showed a random spatial arrangement of the trees, while in Calabria the structure was mainly clustered (Fig. 2).

Tab. 1 - Location of the study sites and structural parameters in the pure and mixed-species stands of European beech and silver fir in Molise and Calabria. (BA): stand basal area; (BAs): basal area percentage of the main species; (dg): quadratic mean diameter; (hg): quadratic mean height.

Stand
Characteristics
Molise Calabria
beechpure firpure beechmixed firmixed beechpure firpure beechmixed firmixed
Latitude (N) 41°52′ 13.65″ 41°52′ 11.02″ 41°52′ 14.35″ 39°09′ 13.49″ 39°09′ 19.31″ 39°09′ 07.96″
Longitude (E) 14°16′ 37.85″ 14°17′ 26.03″ 14°16′ 51.26″ 16°39′ 53.05″ 16°41′ 2.30″ 16°40′ 12.33″
Aspect (°) 354 317 327 - 343 319 307 -
Altitude (m a.s.l.) 1182 1230 1182 - 1289 1220 1332 -
Slope (°) 21 20 27 - 24 37 28 -
BA (m2 ha-1) 33.6 54.4 45.2 - 59.7 61.5 74.6 -
BAs (%) 92 92 20 78 86 85 33 61
No. trees ha-1 989 486 134 202 641 750 293 293
dg (cm) 21 37 28 58 35 32 35 47
hg (m) 15 21 19 27 22 21 22 27
Soil depth (cm) 20 50 40 - 50 70 90 -

  Enlarge/Reduce  Open in Viewer

Fig. 2 - Spatial arrangement of trees in pure and mixed-species stands in Molise (A: silver fir; B: European beech; C: mixed-species) and in Calabria (D: silver fir; E: European beech; F: mixed-species). Green indicate silver fir and violet European beech. A different colour is assigned to each species (minor species included).

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Stands have similar exposure (North, North-West) and slope, between 20° and 28°, except for pure stand of silver fir in Calabria (37°), but differ for soil type and depth. In Molise, the main soil type is Rendzina, whereas in Calabria is Distrudept. Soils in pure stands of European beech are relatively shallow and rocky in comparison with those of the mixed-species stands and the pure stands of silver fir (Tab. 1). All triplets have not been thinned since 1979-1980.

Field measurements

At the beginning of summer 2017, the structural attributes and dendrometric parameters were determined in each stand, for all living trees with a diameter ≥ 7 cm, through the Field-Map® technology and software (Institute of Forest Ecosystem Research Ltd. Jílové u Prahy, Czech Republic - ⇒ http:/­/­www.­fieldmap.­cz/­). For each tree, the stem diameter at 1.30 m, the tree height, the crown radius in the four cardinal directions, the tree species, and the spatial position were recorded. Cores were sampled only from dominant and sub-dominant trees of European beech and silver fir.

Two cores per tree were collected at breast height (1.3 m) perpendicularly to the slope, in order to avoid compression wood (for silver fir) and tension wood (for European beech). Ring-width series were measured with a LINTAB measuring table with an accuracy of 0.01 mm, equipped with a Leica MS5® stereoscope (Leica Microsystems, Wetzlar, Germany). The analysis was performed with the TSAP software package ([57]) and subsequently were cross-dated using COFECHA ([31], [28]) to assess the measurement and dating accuracy. We successfully cross-dated 526 cores in total.

The raw chronologies were converted into tree basal area increment (BAI), according to the following standard equation (eqn. 1):

\begin{equation} BAI= \pi \left (r_{t}^{2} - r_{t-1}^{2} \right ) \end{equation}

where r is the radius of the tree and t is the year of tree-ring formation.

Individual-based competition indices

To analyze the competitive dynamics within the stands, we used distance-dependent competition indices. Generally, although not for stands with stationary structures ([70]), these indices are considered useful for explaining competitive interactions between individual trees ([62]). Distance-dependent indices, in addition to the individual tree sizes, may also provide spatial information on the competitive status of an individual tree ([68]).

To calculate the distance-dependent indices, several methods can be applied ([30], [17], [66]). Three distance-dependent competition indices (CI) based on stem diameter at breast height (DBH) were used in this study: Hegyi ([30] - eqn. 2), Martin & Ek ([40] - eqn. 3) and Rouvinen & Kuuluvainen ([58] - eqn. 4).

The competition indices were calculated using the following formulas:

\begin{equation} CI_{Heg} =\sum_{j=1}^{n} {\frac{DBH_{j} / DBH_{i}} {dist_{ij} +1}} \end{equation}
\begin{equation} CI_{ME} =\sum_{j=1}^{n} { \left [\frac{DBH_{j}} {DBH_{i}} \cdot { e} ^ {- \left (\frac{16 dist_{ij}} {DBH_{i} + DBH_{j}} \right )} \right ]} \end{equation}
\begin{equation} CI_{RK}=\sum_{j=1}^{n} {\frac{DBH_{j} / DBH_{i}} {{ \left (dist_{ij} \right )} ^ {2}}} \end{equation}

In the above equations, DBHi is the diameter at breast height of the subject tree i, DBHj the diameter at breast height of the competitor tree j, distij the distance between the subject tree i and the competitor tree j, and n is the number of competitors in the neighborhood zone.

In order to calculate the competition indices only European beech and silver fir were considered, being the overwhelming species in the sampled stands. The competition indices were calculated at the beginning of the growing period (2011), which allowed us to build the mixed model (see below). In order to obtain the DBHs for 2011, the basal area increments of all European beech and silver fir trees in mixed-species and pure stands were reconstructed. In detail, basal area increments in the last 5 years (2011-2016) and the thickness of the bark (in some cores also the beginning of growth in 2017) of all cored trees were considered. The 5-year basal area increment (BAI5years) and bark thickness were adapted to DBH measured in 2017 of cored trees, using linear models for each site and each species. This allowed to estimate the 5-year basal area increment (BAI5years) and the bark thickness of non-cored trees, and values were sequentially subtracted from the DBH measured in 2017 to reconstruct DBH in 2011 of non-cored European beech and silver fir trees.

The selection of competitors was determined by the fixed radius method. The competition indices (CI) values were obtained using neighborhood radii ranging 1-10 m (in 1 m increments) from the subject tree. Linear regression models were used to identify the optimal search radius and the corresponding CI was used in the subsequent analyses. Linear regression models were tested and evaluated on the basis of the constancy of variances (p > 0.05), adjusted r2, and significance of the correlation (p < 0.1). The BAI5years-CI relationship for silver fir trees was calculated only for trees with DBH ≥ 20 cm, given that the cored trees had a DBH ≥ 20 cm. European beech trees, however, had both a DBH < 20 cm and DBH ≥ 20 cm. Therefore, to assess whether the optimal search radius in European beech varied between the considered DBH values, the BAI5years-CI relationship was calculated both for all cored trees (20 cm > DBH ≥ 20 cm) and for cored trees with a DBH ≥ 20 cm only.

For each tree, the CI was divided into inter-specific and intra-specific competition CI, when trees of different species from the focal tree species (CIinter) and trees of the focal species (CIintra) were considered, respectively.

Modelling approach

To investigate the effects of competitive interactions between mixed-species and pure stands, a mixed-effects model was applied. Specifically, the basal area growth of the individual trees (cm2 year-1) was chosen as the dependent variable, being less influenced by the age of the tree than the ring-width time series ([8]). A natural logarithmic transformation of the original dependent variable was performed, usually being one of the best variables reflecting the non-linear tree growth curve ([71]).

In the model, we considered the hypothesis that the basal area increment of the subject tree is a function of the tree size (SIZE), of the competition (COMP), of site characteristics (SITE), and of species in admixture (MIXTURE - eqn. 5 - eqn. 9:

\begin{equation} \ln (BAI_{5years} )= a + b\, SIZE + c\,COMP + d\,SITE + e\,MIXTURE \end{equation}

where the dependent variable is the 5-year basal area increment (BAI5years, cm2). The four independent variables were defined as:

\begin{equation} SIZE = b_1 \ln (DBH) + b_2 {DBH} ^ {2} + b_3 \left ( \frac{1}{Age} \right ) \end{equation}

where the variable DBH (cm) is that measured at the beginning of the growing period (2011), and Age is the tree age in the year 2011;

\begin{equation} COMP=c_1 CI_{intra} + c_2 CI_{inter} + c_3 mix-pure \end{equation}

where the competition variables are the intra-specific (CIintra) and inter-specific (CIinter) competition indices calculated at the beginning of the growing period (2011) of European beech and silver fir trees in mixed-species and pure stands; mix-pure is the dummy variable, with mix-pure = 0 for pure species stands and mix-pure = 1 for mixed-species stands;

\begin{equation} SITE= d_1 ALT+ d_2 SLOPE+ d_3 DEPTH \end{equation}

where the site variables are topographic features of the stand, including elevation (ALT), and slope and depth of the soil;

\begin{equation} MIXTURE= e_1 BA_{beech} + e_2 BA_{fir} \end{equation}

where BAbeech and BAfir are the basal area (measured at the beginning of the growing period - 2011) of European beech and silver fir, respectively.

The models (silver fir mixed-pure and European beech mixed-pure) were developed with competition variables conveying intra- and inter-specific species interactions. The models for silver fir in mixed-species and pure stands were built with 55 trees, while the models for European beech in mixed-species and pure stands were built with 50 trees. To test whether there was variability in the basal area increment between stands, we fitted a linear mixed model, including the stand (intercept) as a random effect.

Prior to fit the models through the “lmer” function of the R package “lme4” ([5]), the incidence of multicollinearity between the independent variables was evaluated, with the “VIF” function of the R package “car” ([24]). The threshold used to assess multicollinearity was VIF = 2.5 ([19]). After sub-setting the independent variables according to multicollinearity, models were fitted through the likelihood criterion (ML - [46]) and compared using the Akaike’s information criterion (AIC - see Appendix 1 in Supplementary material). We discarded redundant explanatory variables through the “chisq.test” function of the standard “statistics” package in R. The best models were fitted using the limited maximum likelihood method (REML), with a significance threshold of p < 0.05. Normality and homogeneity of the residuals of the final models were visually inspected by plotting the residuals against fitted values ([72] - see Fig. S1 in Supplementary material). Conditional R2 was obtained through the “r.squared GLMM” function of the “MuMIn” package ([4]).

Statistical analyses

The differences between mixed-species and pure stands, in terms of competition and basal area increment, were assessed through the analysis of variance (ANOVA). Analyses were performed for each stand to identify significant correlations between the competition indices calculated at the beginning of the growing period (inter- and intra-specific) and basal area increment. For multiple comparisons of subgroups of variables, Tukey’s post-hoc test was performed. All statistical analyses were performed using the R software package (⇒ https:/­/­www.­r-project.­org/­).

  Results 

Stand structure

Stands showed differences both in tree age (Tab. 2) and stem diameter distribution (Fig. 3). Age differences were observed between Molise and Calabria, as generally younger trees were found in Calabria (Tab. 2).

Tab. 2 - Tree age (range), average diameter at breast height (DBH) of sampled trees, 5-year basal area increment (BAI5years) with the standard deviation (SD), and number of trees and cores used in the study.

Site Stand Tree age
(yrs)
DBH ± SD
(cm)
BAI5years ± SD
(cm2 yr-1)
Cores/trees
Molise beechpure 50-140 22 ± 6 6.71 ± 5.32 110/60
firpure 57-171 42 ± 8 13.58 ± 18.55 86/44
beechmixed 50-119 27 ± 11 17.26 ± 17.32 46/31
firmixed 76-161 59 ± 14 31.24 ± 18.26 38/19
Calabria beechpure 50-109 41 ± 15 13.21 ± 13.35 78/46
firpure 39-82 37 ± 12 15.57 ± 14.52 82/43
beechmixed 53-128 40 ± 14 14.53 ± 14.45 38/25
firmixed 51-108 52 ± 16 27.53 ± 17.12 46/28

  Enlarge/Reduce  Open in Viewer

Fig. 3 - DBH distribution classes of all species in pure and mixed-species stands in Molise (A: silver fir; B: European beech; C: mixed-species) and in Calabria (D: silver fir; E: European beech; F: mixed-species). The size of DBH classes was 10 cm both in pure and mixed-species stands. Tree species are referred to by colours in the legend.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Considering all the species in each stand, diameter distribution showed a negative exponential pattern in almost all sites, except for pure silver fir in Molise (Fig. 3). DBH had a continuous distribution, with a maximum value of 80 cm, except for silver fir (maximum DBH of 60 cm - Fig. 3A) and European beech (40 cm - Fig. 3B) in pure stands, in Molise. Variation in tree height also showed differences between the stands in Molise and those in Calabria (Fig. S2 in Supplementary material). In Molise, trees of silver fir were taller than those of European beech. In Calabria, instead, only in mixed-species stand, trees of silver fir grew taller than those of European beech, the first species occupying the dominant layer while the second the sub-dominated layer.

Selection of competition index and identification of optimal search radius

Among the competition indices that were tested, the Hegyi index (CIHeg), based on the relationship between the stem radial diameters and the distances between the subject tree and the competitors, explained the highest proportion of variability in basal area increment (Tab. 3). The Martin and Ek index (CIME) was the second best, while the Rouvinen and Kuuluvainen (CIRK) was the least appropriate for these stands. CIHeg explained a high proportion of variability in basal area increment for European beech (up to 72% in mixed-species and pure stands) than for silver fir, except for European beech in mixed-species stand in Calabria (from 26% to 39%).

Tab. 3 - Range of R2 calculated for relationships between the 5-years basal area increment (BAI5years) of species and individual competition indices. (CIHeg): Hegyi index; (CIME): Martin and Ek index; (CIRK): Rouvinen and Kuuluvainen index.

Site Competition
Index
firpure beechpure firmixed beechmixed
Molise CI Heg 0.15-0.40 0.08-0.50 0.16-0.47 0.44-0.72
CI ME 0.05-0.39 0.00-0.41 0.14-0.40 0.28-0.63
CI RK 0.03-0.21 0.00-0.20 0.01-0.19 0.02-0.14
Calabria CI Heg 0.12-0.45 0.29-0.58 0.00-0.53 0.26-0.39
CI ME 0.03-0.44 0.13-0.21 0.00-0.39 0.08-0-33
CI RK 0.00-0.43 0.04-0.08 0.01-0.20 0.00-0.03

  Enlarge/Reduce  Open in Viewer

The optimal search radii were selected according to the highest prediction accuracy of the linear regression models (i.e., constancy of variances, adjusted r2, and significance of the correlation - Tab. 4) obtained from the selected competition index (CIHeg). The optimal search radii varied according to the species and stands examined (Tab. 4). In Molise, the highest accuracies of the linear regression model for silver fir were reached for optimal search radii of 10 m. In Calabria, silver fir in pure and mixed-species stands, however, reached the highest prediction accuracy of the linear regression model for an optimal search radius of 9 m and 8 m, respectively.

Tab. 4 - Optimal search radii selected on the basis of the best values of linear regression models, considering DBH ≥ 20 cm. (CV): Constancy of variances (p > 0.05); (r2): adjusted r2; (p): significance of the correlation (p < 0.1); (n): number of tree samples; R (m): optimal radius.

Sites Stands CV r 2 p n R (m)
Molise beechpure 0.80 0.47 < 0.01 22 7
firpure 1.76 0.36 <0.01 18 10
beechmixed 0.13 0.69 <0.01 13 9
firmixed 0.45 0.43 <0.01 17 10
Calabria beechpure 0.78 0.38 0.03 10 7
firpure 3.14 0.39 <0.01 14 9
beechmixed 3.38 0.30 0.07 9 5
firmixed 1.12 0.46 0.02 9 8

  Enlarge/Reduce  Open in Viewer

Considering only trees with DBH ≥ 20 cm, in Molise the optimal search radius for European beech in pure and mixed-species stands was 7 m and 9 m, respectively (Tab. 4). In Calabria, the optimal search radii for European beech in pure and mixed-species stands was 7 m and 5 m, respectively (Tab. 4). The same optimal search radii were identified also for all sampled trees inside the stand (Tab. S1).

Individual-level competition and BAI of mixed-species vs. pure stands

Tab. 5shows the competition indices (CIHeg) and the 5-years basal area increment calculated for each stand. Both intra-specific (CI2011_intra) and inter-specific (CI2011_inter) indices showed higher levels of competition in pure stands of Molise than in pure stands of Calabria. European beech in mixed-species stands had higher values of intra- and inter-specific competition in Molise (0.44 ± 0.25 and 0.99 ± 0.32, respectively) than in Calabria (0.30 ± 0.33 and 0.48 ± 0.35, respectively). In contrast, silver fir in mixed-species stands showed higher values of intra- and inter-specific competition in Calabria (0.83 ± 0.43 and 0.49 ± 0.26, respectively) than in Molise (0.49 ± 0.39 and 0.34 ± 0.11, respectively).

Tab. 5 - Intra- and inter-specific CIHeg in the mixed-species stands for European beech and silver fir at the beginning of the growing period (CI2011_intra and CI2011_inter) and basal area increment of the last 5 years (BAI5years) for the analyzed species.

Site Stand CI2011_intra CI2011_inter BAI5years
(cm2 yr-1)
Molise beechpure 1.83 ± 0.57 0.09 ± 0.14 8.57 ± 4.59
firpure 1.89 ± 0.57 0.10 ± 0.12 10.12 ± 8.13
beechmixed 0.44 ± 0.25 0.99 ± 0.32 19.52 ± 10.93
firmixed 0.49 ± 0.39 0.34 ± 0.11 33.10 ± 19.55
Calabria beechpure 0.79 ± 0.54 0.02 ± 0.02 22.98 ± 13.01
firpure 1.94 ± 0.47 0 14.23 ± 9.42
beechmixed 0.30 ± 0.33 0.48 ± 0.35 17.13 ± 17.36
firmixed 0.83 ± 0.43 0.49 ± 0.26 28.13 ± 17.17

  Enlarge/Reduce  Open in Viewer

The 5-years basal area increment (BAI5years) in pure stands was higher in Calabria than in Molise. Silver fir (33.10 ± 19.55) and European beech (19.52 ± 10.93) in mixed-species stands of Molise showed higher BAI5years values than in pure stands. In Calabria, instead, only silver fir in mixed-species stand (28.13 ± 17.17) had higher BAI5years values than in pure stand (14.23 ± 9.42).

Growth models

The selection of the best intra-specific and inter-specific interaction models was based on the lowest values of AIC (Tab. S2 and Tab. S3 in Supplementary material).

Intra-specific interactions models

The best intra-specific interaction models for silver fir mixed-pure, four independent variables were used: CI2011_intra, slope, DBH2 and mix-pure (Tab. S4). In the model for European beech mixed-pure, six independent variables were used: CI2011_intra, BAbeech, age, soil depth, DBH2 and mix-pure (Tab. S4). In the final model for silver fir mixed-pure, only CI2011_intra affected the basal area increment significantly (Tab. 6); the basal area increment decreased with increasing CI2011_intra. Furthermore, the model showed differences between stands, with R2 of 0.51, when all effects were included, and R2 of 0.50, when only fixed effects were considered.

Tab. 6 - Parameter estimates and fitting statistics for the best intra-specific interaction models through the limited maximum likelihood method (REML). (CI2011_intra): intra-specific interactions; (DBH2): diameter at breast height squared (cm2); (Soil_depth): depth of the soil; (VarPlot): Plot (intercept) variance; (VarRes): Residual variance; (SD): standard deviation; (AIC); Akaike’s information criterion (mixed effect).

Stand
Type
Fixed effects Random effects
Factor Est. SE Prob. VarPlot
(SD)
VarRes
(SD)
R2m R2c AIC
Silver fir
mixed-pure
CI2011_intra -0.693 0.102 <0.001 0.011
(0.106)
0.465
(0.682)
0.500 0.510 126.537
European beech
mixed-pure
CI2011_intra -0.667 0.122 <0.001 0.000
(0.000)
0.513
(0716)
0.600 0.600 147.014
DBH 2 0.595 0.12 <0.001
Soil_depth -0.45 0.126 <0.001

  Enlarge/Reduce  Open in Viewer

In the final model for European beech mixed-pure, CI2011_intra, DBH2 and Soil_ depth affected the basal area increment significantly (Tab. 6); the basal area increment decreased with increasing CI2011_intra and Soil_depth. DBH2 showed direct responses to basal area increment. Furthermore, the model did not show differences between stands, with an R2 of 0.60, both when all effects were included and if only fixed effects were considered.

Basal area increment (BAI5years) revealed that only silver fir in mixed-species stand of Molise differed significantly from the corresponding pure stand (p < 0.05 - Fig. 4A). In both Molise and Calabria, the basal area increments of European beech did not show significant differences between mixed-species and pure stands (Fig. 5A). The competition indices (CI2011_intra), on the other hand, differed significantly between pure and mixed-species stands of Molise and Calabria (p < 0.05 - Fig. 4B, Fig. 5B). In Molise and Calabria, silver fir in mixed-species stands differed significantly from silver fir in pure stands (p < 0.05 - Fig. 4B). Values of CI2011_intra were higher in mixed-species and pure stands of European beech in Molise than in Calabria (Fig. 5B).

Fig. 4 - Basal area increment (logBAI) and competition indices (CI2011_intra) of silver fir in mixed-species and pure stands, considering the last 5 years and the beginning of the growing period. Smixed refers to mixed-species stands (Calabria and Molise) and Spure to pure stands (Calabria and Molise). Capital letters indicate differences between the sites for basal area increment and competition indices, and small letters indicate differences between species (ANOVA and post-hoc tests, p < 0.05). The vertical lines represent the “whiskers” for the 5th and 95th percentiles of data distribution.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Fig. 5 - Basal area increment (logBAI) and competition indices (CI2011_intra) of European beech in mixed-species and pure stands, considering the last 5 years and the beginning of the growing period. Bmixed refers to mixed-species stands (Calabria and Molise) and Bpure to pure stands (Calabria and Molise). Capital letters indicate differences between sites for each competition index, intra- and inter-specific (ANOVA and post-hoc tests, p < 0.05). The vertical lines represent the “whiskers” for the 5th and 95th percentiles of data distribution.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Inter-specific interactions models

Even for inter-specific interactions, the best models for silver fir mixed-pure, four independent variables were used: CI2011_inter, slope, DBH2 and mix-pure (Tab. S5). In the model for European beech mixed-pure, six independent variables were used: CI2011_inter, BAbeech, age, Soil_depth, DBH2 and mix-pure (Tab. S5). In the final model for silver fir mixed-pure, only DBH2 significantly affected the basal area increment (Tab. 7). DBH2 showed direct responses to tree basal area increment. The model showed also differences between stands, with an R2 of 0.46, when all effects were included, and R2 of 0.36, when only fixed effects were considered.

Tab. 7 - Parameter estimates and fitting statistics for the best inter-specific interaction models through the limited maximum likelihood method (REML). (DBH2): diameter at breast height squared (cm2); (Age): tree age (years); (Soil_depth): depth of the soil; (BAbeech): basal area occupied by competitors of European beech; (VarPlot): Plot (intercept) variance; (VarRes): Residual variance; (SD): standard deviation; (AIC); Akaike’s information criterion (mixed effect).

Stand
Type
Fixed effects Random effects
Factor Est. SE Prob. VarPlot
(SD)
VarRes
(SD)
R2m R2c AIC
Silver fir
mixed-pure
DBH 2 0.562 0.114 <0.001 0.093
(0.305)
0.479
(0.692)
0.360 0.460 130.837
European beech
mixed-pure
DBH 2 0.678 0.125 <0.001 0.000
(0.000)
0.564
(0.751)
0.570 0.570 141.972
Age 0.192 0.113 0.097
Soil_depth -0.456 0.135 0.001
BA beech -0.541 0.122 <0.001

  Enlarge/Reduce  Open in Viewer

In the final model for European beech mixed-species and pure stands, BAbeech, DBH2 and Soil_depth affected the basal area increment significantly (Tab. 7); the basal area increment decreased with increasing BAbeech and Soil_depth. DBH2 showed direct responses to basal area increment. The model did not show differences between stands, with an R2 of 0.57, both when all effects were included and if only fixed effects were considered.

  Discussion 

Influence of intra-specific interactions on trees basal area increment in mixed-species and pure stands

We analyzed the effect of interactions between tree species on the basal area increment of uneven-aged mixed-species and pure stands of European beech and silver fir, implementing growth models and competition indices at the individual tree level. The final models showed strong influence of intra-specific interactions (CI2011_intra) on the basal area increment (BAI5years) of both species (Tab. 6). The dummy variable mix-pure (Tab. S4 in Supplementary material) indicates the influence of variables related to the pure and mixed-species stands on the basal area increment (BAI5years). The absence of significance of this variable, suggests that intra-specific interactions (CI2011_intra) had high influence on the basal area increment (BAI5years) of both species in mixed-species and pure stands. Therefore, European beech and silver fir in mixed-species and pure stands were more affected by competition with conspecific than allospecific neighbors. Several studies indicated that European beech is a strong self-competitor ([49], [51], [10]), with high intra-specific competition for above-ground ([55]) and below-ground ([56], [18], [44]) resources. In literature studies on silver fir focused mainly on drought response ([34]), while are rather scarce in terms of competitive interactions between mixed-species and pure stands.

Analysis of variance (ANOVA) showed that the influence of intra-specific interactions (CI2011_intra) was higher in pure stands of silver fir and European beech than in mixed-species stands (Fig. 4B, Fig. 5B), supporting our first hypothesis. In Molise, ANOVA highlighted significantly higher increase in the basal area increment (BAI5years) of silver fir in mixture than in monoculture (Fig. 4A). In accordance with our results, Toïgo et al. ([65]) observed an increase in the basal area increment of silver fir more marked in mixed-species stands than in pure stands. Crown architecture, vertical positioning within the canopy and root morphology may lead to differences in resources accessibility between pure and mixed-species stands, and, therefore, a more efficient use of water and light ([36], [23]). However, in our case, higher basal area increments of silver fir trees in the mixed-species stand in Molise, compared to the corresponding pure stand (Fig. 4A) might depend on differences in stand structure (Fig. 3). Silver fir trees in the mixed-species stand were bigger (DBH > 50 cm) than in the pure stand (DBH < 50 cm), which might explain the reduction in competitive interactions and, consequently the higher basal area increment of silver fir trees in the mixed-species stand.

European beech, instead, showed high tree-to-tree competition in pure stands, most probably related to the low self-tolerance of this species ([47]), which is attenuate in mixed-species stands ([50]). For example, Mölder & Leuschner ([45]) found a lower impact of intra-specific interactions on the basal area increment of European beech in mixed-species forests than in pure stands. This less negative impact of intra-specific interactions was related to the lower number of European beech trees among competitors, confirming that European beech is a strong self-competitor. High levels of intra-specific interactions may also depend on high degree of canopy closure in pure stands of European beech ([49]). In our case, significantly higher intra-specific competition levels (CI2011_intra) were observed in pure stand of European beech, in Molise, which had closed canopy (Fig. 2B), in comparison with the corresponding one in Calabria (Fig. 2E). A closed-canopy may derive from different spatial arrangement of trees, random (in Molise) vs. clustered (in Calabria). In fact, mixed-species and pure stands of European beech showed higher influence of intra-specific interactions (CI2011_intra) in Molise than in Calabria (Fig. 5B). Gap recruitment of small competing trees leading to clustered regeneration may occur more frequently in drier sites ([61]). Edaphic heterogeneity and local aridity may also cause clumped spatial distribution of trees.

Several studies found that the stand-level spatial arrangement of trees may influence the competitive intensity between species ([25]), especially in uneven-aged and uneven-structured stands ([69]), as well as in those stands with trees differing in crown morphology, crown efficiency and space occupation ([42]). This would suggest that the forest structure (e.g., spatial arrangement, crown morphology, etc.) may regulate tree-to-tree competition and tree growth in structurally diverse forests ([25]).

Intra- and inter-specific interactions in mixed-species stand

Inter-specific interaction (Tab. 7), was not the most important variable to explain the competitive interactions in mixed stands. This suggests that inter-specific interactions (CI2011_inter) in mixed-species stands did not significantly influence the basal area increment, supporting our second hypothesis. Gazol et al. ([27]) observed higher growth rates in silver fir admixed with European beech. Other studies of intra- and inter-specific interactions on fir-beech mixtures indicated either no effects of mixing on basal area increment of either species ([18], [59]) or positive mixing effects for beech, whereas positive effects for silver fir could be restricted to sites of lower productivity ([65]). A negative influence of intra-specific interactions (in comparison with inter-specific interactions) on the growth of European beech was confirmed in the present study. However, the effect of these interactions on the growth of silver fir needs to be further explored, in view of contradicting results of previous studies ([18], [10], [44]). Indeed, Mina et al. ([44]) found that the growth of silver fir was negatively affected by the presence of European beech. In contrast, Del Río et al. ([18]) revealed that silver fir was more affected by intra-specific than inter-specific competition, suggesting the presence of some degree of complementarity between European beech and silver fir. Bošela et al. ([10]) noted that the growth of silver fir was positively influenced by admixture with European beech, showing facilitating effect or a reduction in competition of European beech on silver fir. Evidence of inter-specific facilitation or reduction in tree-to-tree competition could be linked to differences in functional traits and phenological patterns of silver fir vs. European beech ([18]). For example, compartmentalization of resource utilization may allow a different strategy for root water uptake ([34]). Complementary crown architectures and contrasting physiological adaptations may facilitate canopy expansion (European beech) and improve crown light interception in spring and autumn (silver fir - [38]). This suggests that, in mixed-species stands, tree species with complementary traits may increase the exploitation of resources in comparison with those species with similar traits ([55]). According to Pretzsch & Schütze ([55]), improved resource use efficiency (e.g., light, water and nutrients) in mixed-species stands is attained when a complex vertical structure (small and large individuals of both species) occurs. In our case, the absence of inter-specific interactions (CI2011_inter) can be related to the smaller diameter in European beech than in silver fir, in the mixed-species stands of Molise and Calabria (Fig. 3). These differences between European beech and silver fir were also observed for the variation in tree height, silver fir being taller than European beech (Fig. S2 in Supplementary material). Therefore, the sub-dominant structural position of European beech within mixed-species stands may determine a reduction in competition for silver fir.

Overall, forest management that promotes multi-layered uneven-sized species mixtures might allow maintaining these Mediterranean mountain forests continuously, increasing their resilience and productivity. Although the importance of species mixture is recognized, management of mixed-species uneven-aged forests at the range edge of species distribution needs to be optimized for effective application. Indeed, the sustainability of species mixture might become an issue in drought-prone environments ([15]), in which a stable species distribution could be negatively affected by increasing disturbance frequency and intensity. For a continuous delivery of ecosystem services of mixed-species mountain forests, the selective felling might regulate the occupation of the growing space of the species and enhance the development of the stands towards a multi-layered structure, making them more stable against natural hazards and ensuring the protective functions of the forests, for soils and biodiversity, as well as the maintenance of carbon stocks.

  Conclusion 

The growth models showed that intra-specific interactions influenced the basal area increment more than inter-specific interactions. In particular, intra-specific interactions negatively influenced the basal area increment of European beech and silver fir, both in pure and mixed-species stands of Molise and Calabria. In accordance with our first hypothesis, higher intra-specific interactions were observed in pure stands of both species than in corresponding mixed-species stands. Silver fir showed relatively lower levels of intra-specific competition, with higher basal area increment in mixed-species stands. This higher basal area increment of individual trees in mixed-species stands can be explained by the distinct stand structure, which may have increased space occupation and improved resource utilization. For European beech, high values of intra-specific interactions could be linked to its low self-tolerance and to the spatial arrangement of trees within the stands, following canopy closure. The absence of inter-specific interactions in mixed-species stands was probably related to the sub-dominant position of European beech, which may have reduced the benefit derived from niche separation and complementarity for silver fir. These results support our second hypothesis, providing useful insights on the structural and functional dynamics that occur in uneven-aged mixed-species forests vs. pure stands of European beech and silver fir. This suggest that mixtures of European beech and silver fir, a natural species combination in mountain regions, should be promoted as an option for mountain ranges in southern Europe. However, further studies with a larger number of samples are needed to draw more definitive conclusions, so as to help evaluate forest adaptation measures under global change scenarios.

  Acknowledgements 

The authors are grateful to Carabinieri Forestali of Petilia Policastro (KR, Italy) for their help in selecting the study sites. This study generated from the COST (European Cooperation in Science and Technology) Action CLIMO (“Climate-Smart Forestry in Mountain Regions” - CA15226) financially supported by the EU Framework Programme for Research and Innovation HORIZON 2020.

  References

(1)
Ammer C (2019). Diversity and forest productivity in a changing climate. New Phytologist 221 (1): 50-66.
CrossRef | Gscholar
(2)
Antonucci S, Rossi S, Lombardi F, Marchetti M, Tognetti R (2019). Influence of climatic factors on silver fir xylogenesis along the Italian Peninsula. IAWA journal 40 (2): 259-275.
CrossRef | Gscholar
(3)
Aussenac R, Bergeron Y, Gravel D, Drobyshev I (2019). Interactions among trees: a key element in the stabilising effect of species diversity on forest growth. Functional Ecology 33 (2): 360-367.
CrossRef | Gscholar
(4)
Barton K (2013). MuMln: multi-model inference. R Package Version 1.9.5.
Online | Gscholar
(5)
Bates D, Maechler M, Bolker B, Walker S (2014). lme4: Linear mixed-effects models using Eigen and S4. R package version 1: 1-23.
Online | Gscholar
(6)
Biging GS, Dobbertin M (1992). A comparison of distance-dependent competition measures for height and basal area growth of individual conifer trees. Forest Science 38 (3): 695-720.
Online | Gscholar
(7)
Biging GS, Dobbertin M (1995). Evaluation of competition indices in individual tree growth models. Forest Science 41 (2): 360-377.
Online | Gscholar
(8)
Biondi F (1999). Comparing tree-ring chronologies and repeated timber inventories as forest monitoring tools. Ecological Application 9 (1): 216-227.
CrossRef | Gscholar
(9)
Bošela M, Petráš R, Sitková Z, Priwitzer T, Pajtík J, Hlavatá H, Sedmák R, Tobin B (2014). Possible causes of the recent rapid increase in the radial increment of silver fir in the Western Carpathians. Environmental Pollution 184: 211-221.
CrossRef | Gscholar
(10)
Bošela M, Tobin B, Seben V, Petráš R, Larocque GR (2015). Different mixtures of Norway spruce, silver fir, and European beech modify competitive interactions in central European mature mixed forests. Canadian Journal of Forest Research 45 (11): 1577-1586.
CrossRef | Gscholar
(11)
Brassard BW, Chen HY, Cavard X, Laganière J, Reich PB, Bergeron Y, Paré D, Yuan Z (2013). Tree species diversity increases fine root productivity through increased soil volume filling. Journal of Ecology 101 (1): 210-219.
CrossRef | Gscholar
(12)
Brooker RW (2006). Plant-plant interactions and environmental change. New Phytologist 171 (2): 271-284.
CrossRef | Gscholar
(13)
Coates KD, Lilles EB, Astrup R (2013). Competitive interactions across a soil fertility gradient in a multispecies forest. Journal of Ecology 101 (3): 806-818.
CrossRef | Gscholar
(14)
Condés S, Del Río M, Sterba H (2013). Mixing effect on volume growth of Fagus sylvatica and Pinus sylvestris is modulated by stand density. Forest of Ecology and Management 292: 86-95.
CrossRef | Gscholar
(15)
Conte E, Lombardi F, Battipaglia G, Palombo C, Altieri S, La Porta N, Marchetti M, Tognetti R (2018). Growth dynamics, climate sensitivity and water use efficiency in pure vs. mixed pine and beech stands in Trentino (Italy). Forest of Ecology and Management 409: 707-718.
CrossRef | Gscholar
(16)
Daniels RF (1976). Simple competition indices and their correlation with annual loblolly pine tree growth. Forest Science 22 (4): 454-456.
Online | Gscholar
(17)
Daniels RF, Burkhart HE, Clason TR (1986). A comparison of competition measures for predicting growth of loblolly pine trees. Canadian Journal of Forest Research 16 (6): 1230-1237.
CrossRef | Gscholar
(18)
Del Río M, Condés S, Pretzsch H (2014). Analyzing size-symmetric vs. size-asymmetric and intra-vs. inter-specific competition in beech (Fagus sylvatica L.) mixed stands. Forest of Ecology and Management 325: 90-98.
CrossRef | Gscholar
(19)
Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, García Marquéz JR, Gruber B, Lafourcade B, Leitão PJ, Münkemüller T, McClean C, Osborne PE, Reineking B, SchröDer B, Skidmore AK, Zurell D, Lautenbach S (2013). Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography 36 (1): 27-46.
CrossRef | Gscholar
(20)
Forrester DI, Vanclay JK, Forrester RI (2011). The balance between facilitation and competition in mixtures of Eucalyptus and Acacia changes as stands develop. Oecologia 166 (1): 265-272.
CrossRef | Gscholar
(21)
Forrester DI (2014). The spatial and temporal dynamics of species interactions in mixed-species forests: from pattern to process. Forest of Ecology and Management 312: 282-292.
CrossRef | Gscholar
(22)
Forrester DI, Albrecht AT (2014). Light absorption and light-use efficiency in mixtures of Abies alba and Picea abies along a productivity gradient. Forest of Ecology and Management 328: 94-102.
CrossRef | Gscholar
(23)
Forrester DI, Bauhus J (2016). A review of processes behind diversity productivity relationships in forests. Current Forestry Reports 2 (1): 4-61.
CrossRef | Gscholar
(24)
Fox J, Weisberg S, Adler D, Bates D, Baud-Bovy G, Ellison S, Firth D, Friendly M. Gorjanc G, Graves S (2012). Package “car”. R Foundation for Statistical Computing, Vienna, Austria.
Online | Gscholar
(25)
Fraver S, D’Amato AW, Bradford JB, Jonsson BG, Jönsson M, Esseen PA (2014). Tree growth and competition in an old-growth Picea abies forest of boreal Sweden: influence of tree spatial patterning. Journal of Vegetation Science 25 (2): 374-385.
CrossRef | Gscholar
(26)
Gamfeldt L, Snäll T, Bagchi R, Jonsson M, Gustafsson L, Kjellander P, Riuz-Jaen MC, Fröberg M, Stendahl J, Philipson CD, Mikusinski G, Andersson E, Westerlund B, Andrén H, Moen J, Bengtsson J (2013). Higher levels of multiple ecosystem services are found in forests with more tree species. Nature Communications 4 (1): 59.
CrossRef | Gscholar
(27)
Gazol A, Camarero J, Gomez-Aparicio JL (2016). Functional diversity enhances silver fir growth resilience to an extreme drought. Journal of Ecology 104: 1063-1075.
CrossRef | Gscholar
(28)
Grissino-Mayer HD (2001). Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA. Tree Ring Research 57: 205-221.
Online | Gscholar
(29)
Grossiord C, Granier A, Ratcliffe S, Bouriaud O, Bruelheide H, Chećko E, Forrester DI, Dawud SM, Finér L, Pollastrini M, Scherer-Lorenzen M, Valladares F, Bonal D, Gessler A (2014). Tree diversity does not always improve resistance of forest ecosystems to drought. Proceedings of the National Academy of Sciences USA 111 (41): 14812-14815.
CrossRef | Gscholar
(30)
Hegyi F (1974). A simulation model for managing jack-pine stands simulation (Fries J ed). Royal College of Forestry, Stockholm, Sweden, pp. 74-90.
Gscholar
(31)
Holmes RL, Adams RK, Fritts HC (1986). Tree-ring chronologies of western North America: California, eastern Oregon and northern Great Basin with procedures used in the chronology development work including users manuals for computer programs COFECHA and ARSTAN. Chronology Series VI, University of Arizona, Tucson, AZ, pp. 185.
Online | Gscholar
(32)
Jactel H, Brockerhoff EG (2007). Tree diversity reduces herbivory by forest insects. Ecology Letters 10 (9): 835-848.
CrossRef | Gscholar
(33)
Klein T (2014). The variability of stomatal sensitivity to leaf water potential across tree species indicates a continuum between isohydric and anisohydric behaviours. Functional Ecology 28: 1313-1320.
CrossRef | Gscholar
(34)
Lebourgeois F, Gomez N, Pinto P, Mérian P (2013). Mixed stands reduce Abies alba tree-ring sensitivity to summer drought in the Vosges mountains, western Europe. Forest of Ecology and Management 303: 61-71.
CrossRef | Gscholar
(35)
Liang J, Crowther TW, Picard N, Zhou M, Alberti G, Schulze E-D, McGuire AD, Bozzato F, Pretzsch H, De-Miguel S, Paquette A, Hérault B, Scherer-Lorenzen M, Barrett CB, Glick HB, Hengeveld GM, Nabburs G-J, Pfautsch S, Viana H, Vibrans AC, Ammer C, Schall P, Verbyla D, Tchebakova N, Fischer M, Watson JV, Chen HYH, Lei X, Schelhaas M-J, Lu H, Gianelle D, Parfenova EI, Salas C, Lee E, Lee B, Kim HS, Bruelheide H, Coomes DA, Piotto D, Sunderland T, Schmid B, Gourlet-Fleury S, Sonké B, Tavani R, Zhu J, Brandl S, Vayreda J, Kitahara F, Searle EB, Neldner VJ, Ngugi MR, Baraloto C, Frizzera L, Balazy R, Oleksyn J, Zawila-Niedzwiecki T, Bouriaud O, Bussotti F, Finér L, Jaroszewicz B, Jucker T, Valladares F, Jagodzinski AM, Peri PL, Gonmadje C, Marthy W, O’Brien T, Martin EH, Marshall AR, Rovero F, Bitariho R, Niklaus PA, Alvarez-Loayza P, Chamuya N, Valencia R, Mortier F, Wortel V, Engone-Obiang NL, Ferreira LV, Odeke DE, Vasquez RM, Lewis SL, Reich PB (2016). Positive biodiversity-productivity relationship predominant in global forests. Science 354 (6309): aaf8957.
CrossRef | Gscholar
(36)
Linares JC, Camarero JJ, Carreira JA (2010). Competition modulates the adaptation capacity of forests to climatic stress: insights from recent growth decline and death in relict stands of the Mediterranean fir Abies pinsapo. Journal of Ecology 98 (3): 592-603.
CrossRef | Gscholar
(37)
Lorimer CG, Dahir SE, Nordheim EV (2001). Tree mortality rates and longevity in mature and old-growth hemlock-hardwood forests. Journal of Ecology 89 (6): 960-971.
CrossRef | Gscholar
(38)
Lu H, Mohren GM, Den Ouden J, Goudiaby V, Sterck FJ (2016). Overyielding of temperate mixed forests occurs in evergreen-deciduous but not in deciduous-deciduous species mixtures over time in the Netherlands. Forest of Ecology and Management 376: 321-332.
CrossRef | Gscholar
(39)
Magh RK, Bonn B, Grote R, Burzlaff T, Pfautsch S, Rennenberg H (2019). Drought superimposes the positive effect of silver fir on water relations of European beech in mature forest stands. Forests 10 (10): 897.
CrossRef | Gscholar
(40)
Martin GL, Ek AR (1984). A comparison of competition measures and growth models for predicting plantation red pine diameter and height growth. Forest Science 30 (3): 731-743.
Online | Gscholar
(41)
Marziliano PA, Tognetti R, Lombardi F (2019). Is tree age or tree size reducing height increment in Abies alba Mill. at its southernmost distribution limit? Annals of Forest Science 76: 17.
CrossRef | Gscholar
(42)
Metz J, Seidel D, Schall P, Scheffer D, Schulze ED, Ammer C (2013). Crown modeling by terrestrial laser scanning as an approach to assess the effect of aboveground intra-and interspecific competition on tree growth. Forest of Ecology and Management 310: 275-288.
CrossRef | Gscholar
(43)
Mina M, Huber MO, Forrester DI, Thürig E, Rohner B (2018a). Multiple factors modulate tree growth complementarity in central European mixed forests. Journal of Ecology 106 (3): 1106-1119.
CrossRef | Gscholar
(44)
Mina M, Del Río M, Huber MO, Thürig E, Rohner B (2018b). The symmetry of competitive interactions in mixed Norway spruce, silver fir and European beech forests. Journal of Vegetation Science 29 (4): 775-787.
CrossRef | Gscholar
(45)
Mölder I, Leuschner C (2014). European beech grows better and is less drought sensitive in mixed than in pure stands: tree neighbourhood effects on radial increment. Trees 28 (3): 777-792.
CrossRef | Gscholar
(46)
Pinheiro JC, Bates DM (2000). Statistics and computing. Mixed-effects models in S and S-PLUS Statistics and Computing book series, Springer-Verlag, New York, USA, pp. 3-56.
CrossRef | Gscholar
(47)
Pretzsch H (2006). Species-specific allometric scaling under self-thinning: evidence from long-term plots in forest stands. Oecologia 146 (4): 572-583.
CrossRef | Gscholar
(48)
Pretzsch H, Schütze G (2009). Transgressive overyielding in mixed compared with pure stands of Norway spruce and European beech in Central Europe: evidence on stand level and explanation on individual tree level. European Journal of Forest Research 128 (2): 183-204.
CrossRef | Gscholar
(49)
Pretzsch H, Block J, Dieler J, Dong PH, Kohnle U, Nagel J, Spellmann H, Zingg A (2010). Comparison between the productivity of pure and mixed stands of Norway spruce and European beech along an ecological gradient. Annals of Forest Science 67 (7): 712.
CrossRef | Gscholar
(50)
Pretzsch H, Matthew C, Dieler J (2012). Allometry of tree crown structure. Relevance for space occupation at the individual plant level and for self-thinning at the stand level. In: “Growth and Defence in Plants”. Springer, Berlin, Heidelberg, Germany, pp. 287-310.
CrossRef | Gscholar
(51)
Pretzsch H, Schütze G, Uhl E (2013). Resistance of European tree species to drought stress in mixed versus pure forests: evidence of stress release by inter-specific facilitation. Plant Biology 15: 483-495.
CrossRef | Gscholar
(52)
Pretzsch H (2014). Canopy space filling and tree crown morphology in mixed-species stands compared with monocultures. Forest Ecology and Management 327: 251-264.
CrossRef | Gscholar
(53)
Pretzsch H, Biber P, Schütze G, Bielak K (2014). Changes of forest stand dynamics in Europe. Facts from long-term observational plots and their relevance for forest ecology and management. Forest Ecology and Management 316: 65-77.
CrossRef | Gscholar
(54)
Pretzsch H, Biber P (2016). Tree species mixing can increase maximum stand density. Canadian Journal of Forest Research 46 (10): 1179-1193.
CrossRef | Gscholar
(55)
Pretzsch H, Schütze G (2016). Effect of tree species mixing on the size structure, density, and yield of forest stands. European Journal of Forest Research 135 (1): 1-22.
CrossRef | Gscholar
(56)
Rewald B, Leuschner C (2009). Belowground competition in a broad-leaved temperate mixed forest: pattern analysis and experiments in a four-species stand. European Journal of Forest Research 128 (4): 387-398.
CrossRef | Gscholar
(57)
Rinn F (2003). TSAP-win. Time series analysis and presentation for dendrochronology and related applications. User Reference, Heidelberg, Germany.
Online | Gscholar
(58)
Rouvinen S, Kuuluvainen T (1997). Structure and asymmetry of tree crowns in relation to local competition in a natural mature Scots pine forest. Canadian Journal of Forest Research 27 (6): 890-902.
CrossRef | Gscholar
(59)
Schwarz JA, Bauhus J (2019). Benefits of mixtures on growth performance of silver fir (Abies alba) and European beech (Fagus sylvatica) increase with tree size without reducing drought tolerance. Frontiers in Forests and Global Change 2: 1483.
CrossRef | Gscholar
(60)
Seidl R, Spies TA, Peterson DL, Stephens SL, Hicke JA (2016). Searching for resilience: addressing the impacts of changing disturbance regimes on forest ecosystem services. Journal of Applied Ecology 53 (1): 120-129.
CrossRef | Gscholar
(61)
Spies TA, Franklin JF, Klopsch M (1990). Canopy gaps in Douglas-fir forests of the Cascade Mountains. Canadian Journal of Forest Research 20 (5): 649-658.
CrossRef | Gscholar
(62)
Stadt KJ, Huston C, Coates KD, Feng Z, Dale MR, Lieffers VJ (2007). Evaluation of competition and light estimation indices for predicting diameter growth in mature boreal mixed forests. Annals of Forest Science 64 (5): 477-490.
CrossRef | Gscholar
(63)
Tognetti R, Cherubini P, Marchi S, Raschi A (2007). Leaf traits and tree rings suggest different water-use and carbon assimilation strategies by two co-occurring Quercus species in a Mediterranean mixed-forest stand in Tuscany, Italy. Tree Physiology 27: 1741-1751.
CrossRef | Gscholar
(64)
Tielbörger K, Kadmon R (2000). Temporal environmental variation tips the balance between facilitation and interference in desert plants. Ecology 81 (6): 1544-1553.
CrossRef | Gscholar
(65)
Toïgo M, Vallet P, Perot T, Bontemps JD, Piedallu C, Courbaud B (2015). Overyielding in mixed forests decreases with site productivity. Journal of Ecology 103 (2): 502-512.
CrossRef | Gscholar
(66)
Tomé M, Burkhart HE (1989). Distance-dependent competition measures for predicting growth of individual trees. Forest Science 35 (3): 816-831.
Online | Gscholar
(67)
Vilà M, Carrillo-Gavilán A, Vayreda J, Bugmann H, Fridman J, Grodzki W, Haase J, Kunstler G, Schelhaas M, Trasobares A (2013). Disentangling biodiversity and climatic determinants of wood production. PLoS One 8(2): e53530.
CrossRef | Gscholar
(68)
Wang D, Tang SC, Hsieh HC, Chung CH, Lin CY (2012). Distance-dependent competition measures for individual tree growth on a Taiwania plantation in the Liuguei area. Taiwan Journal of Forest Science 27: 215-227.
Gscholar
(69)
Weber P, Bugmann H, Fonti P, Rigling A (2008). Using a retrospective dynamic competition index to reconstruct forest succession. Forest of Ecology and Management 254 (1): 96-106.
CrossRef | Gscholar
(70)
Weiskittel AR, Hann DW, Kershaw JAJ, Vanclay JK (2011). Forest growth and yield modeling. John Wiley and Sons, Chichester, UK, pp. 344.
Online | Gscholar
(71)
Wykoff WR (1990). A basal area increment model for individual conifers in the northern Rocky Mountains. Forest Science 36 (4): 1077-1104.
Online | Gscholar
(72)
Zuur A, Ieno EN, Walker N, Saveliev AA, Smith GM (2009). Mixed effects models and extensions in ecology with R. Statistics for Biology and Health, Springer Science and Business Media, New York, USA, pp. 574.
Online | Gscholar

Authors’ Affiliation

(1)
Soraya Versace
Vittorio Garfì 0000-0002-6980-2920
Di Febbraro Mirko 0000-0001-8898-7046
Department of Biosciences and Territory, University of Molise, c.da Fonte Lappone, 86090 Pesche, IS (Italy)
(2)
Soraya Versace
Michele Dalponte 0000-0001-9850-8985
Lorenzo Frizzera 0000-0003-3150-8395
Damiano Gianelle 0000-0001-7697-5793
Department of Sustainable Agro-ecosystems and Bioresources, Research and Innovation Centre, Fondazione E. Mach, v. E. Mach 1, 38010 San Michele all’Adige, TN (Italy)
(3)
Damiano Gianelle 0000-0001-7697-5793
Roberto Tognetti 0000-0002-7771-6176
The EFI Project Centre on Mountain Forests (MOUNTFOR), v. Edmund Mach 1, 38010, San Michele all’Adige, TN (Italy)
(4)
Roberto Tognetti 0000-0002-7771-6176
Department of Agricultural, Environmental and Food Sciences, University of Molise, v. Francesco De Sanctis, 86100 Campobasso (Italy)

Corresponding author

 

Citation

Versace S, Garfì V, Dalponte M, Febbraro Mirko D, Frizzera L, Gianelle D, Tognetti R (2021). Species interactions in pure and mixed-species stands of silver fir and European beech in Mediterranean mountains. iForest 14: 1-11. - doi: 10.3832/ifor3476-013

Academic Editor

Giorgio Vacchiano

Paper history

Received: Apr 20, 2020
Accepted: Oct 26, 2020

First online: Jan 02, 2021
Publication Date: Feb 28, 2021
Publication Time: 2.27 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2021

  Open Access

This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Creative Commons Licence

Breakdown by View Type

(Waiting for server response...)

Article Usage

Total Article Views: 25527
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 21169
Abstract Page Views: 1707
PDF Downloads: 2317
Citation/Reference Downloads: 4
XML Downloads: 330

Web Metrics
Days since publication: 1172
Overall contacts: 25527
Avg. contacts per week: 152.47

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Feb 2023)

Total number of cites (since 2021): 3
Average cites per year: 1.00

 
 

Publication Metrics

by Dimensions ©

List of the papers citing this article based on CrossRef Cited-by.

 

iForest Similar Articles

 

This website uses cookies to ensure you get the best experience on our website. More info