## Growth dynamics of the Norway spruce and silver fir understory in continuous cover forestry

iForest - Biogeosciences and Forestry, Volume 13, Issue 1, Pages 56-64 (2020)
doi: https://doi.org/10.3832/ifor3183-012

Research Articles

The conversion to uneven-aged, mixed-species stands represents one possible way to mitigate the consequences of disturbances in Norway spruce forests in Central Europe. A better understanding of the establishment and growth dynamics of the understory can contribute to a more effective conversion process. Here we investigate the structure of understory, light climate and growth of natural regeneration of Norway spruce and silver fir in two forest stands undergoing conversion to continuous cover forestry. Stand-wise forest inventory was conducted in 1993 and 2013. The natural regeneration was surveyed, and the light conditions and inter-tree competition were quantified in 51 sample plots established across the stands in 2013. Our results suggest that the diffuse radiation strongly affects the height growth of fir and spruce natural regeneration. We do not confirm the effect of local sapling density on the regeneration dynamics. The results further show that fir trees grow faster than spruce under less intensity of diffuse light. Most of the spruce and fir trees reached the upper limit of the lower overstory (DBH 12 cm) at approximately 50 years of age. Thus, more substantial reductions in stand density can lead to a well-differentiated structure in less than five decades. To control the prospective representation of spruce and fir in mixed-species forests undergoing the conversion, managing of light conditions is crucial.

# Introduction

Most of the stands dominated by Norway spruce (Picea abies [L.] Karst.) in Central Europe originated from traditional even-aged management. They are productive yet rather unstable forest ecosystems that are often threatened by windthrows, bark beetle outbreaks ([44], [32]), nutrient leaching and accelerated soil acidification ([1]) and recently also by climate change ([2], [4], [45]). The negative impact of these factors together with current trends emphasizing the advantages of incorporating close-to-nature principles into silvicultural management encourages the conversion of pure Norway spruce forests into uneven-aged, mixed-species stands ([22]). It is assumed that structurally diverse and mixed-species forests have higher adaptation potential when exposed to unpredictable stress or disturbances ([5]). As a traditional multi-aged silvicultural system, single-tree selection (also known as Plenterwald) represents one of the appropriate alternatives for this transformation ([55]). Stand structures of this kind are considered stable and able to continuously produce high-quality timber. In comparison with even-aged plantations, intensively managed selection forests are usually formed by silver fir (Abies alba Mill.) with an admixture of Norway spruce and European beech (Fagus sylvatica L.) mixed at small spatial scales ([41]). In conifer-dominated mixed stands, silver fir is a species known to be more tolerant to higher temperatures and drought ([52], [54]) and to have greater resistance to wind and insect damage than Norway spruce ([13], [19]); thus, its increased proportion in converted stands can promote forest stability. However, the key assumption is the elimination of browsing by ungulates ([27]).

The process of conversion from homogenous to multi-aged stand structures is a complicated and long-lasting task. Its realization depends on the initial condition of the respective forest stand and, according to Schütz ([41]), the conversion has to consider stand stability, the life expectancy of the canopy trees, structural diversification and the prospects for regeneration establishment. Natural regeneration plays a crucial role in the whole process and is the necessary precondition for forming the overstory. One of the main factors regulating the species composition, spatial distribution and development of natural regeneration is the light environment. Variability in the light microclimate in canopy gaps allows the establishment and growth of tree species with diverse ecological requirements ([20], [50], [47], [15]). In stands undergoing conversion to single-tree selection, low light intensity favours the natural regeneration of semi-shade- (Norway spruce) and shade-tolerant species (silver fir and European beech - [42]). Seedlings and saplings of these tree species are able to survive and grow under low light intensities ([10], [15]). In shaded conditions, silver fir has a competitive advantage because it is able to achieve a positive carbon balance at a lower radiation intensity than Norway spruce ([24], [37]). This advantage is attributed, in particular, to the superior morphological plasticity of fir saplings’ crowns and to the more efficient photosynthetic capacity of fir needles growing under low light levels. Crown morphology therefore significantly influences the utilization of diffuse light, resulting in faster height growth and the differences in competitiveness of spruce and fir in the understory ([39], [26]).

Unlike the regeneration cuts applied in even-aged stands (e.g., regular shelterwood cut), where the light conditions are much more homogenous and bound to certain stand area, the distribution of light in stands undergoing conversion to single-tree selection is spatially diverse, dominated by small sunflecks of short duration that supply almost all light required for natural regeneration ([42]). Heterogeneous vertical canopy structure and irregular distribution of natural regeneration cause high variability of light conditions and modify the strength of inter- and intraspecific competition within the cohorts of natural regeneration ([18], [25], [31], [12]). On the other hand, deliberate regulation of light and competition in understory has a significant impact on subsequent changes in stand structure and tree species composition ([12]). Therefore, detailed insights into the growth processes of understory trees in the specific conditions of stands undergoing conversion may help to manage the process of conversion following principles of close-to-nature forest management more effectively. Under the ongoing climate change, maintaining the appropriate spruce-fir (beech) ratio in regeneration as well as in the mature stands represents the crucial task of the silvicultural treatments. Moreover, for the planning of conversions in vulnerable spruce stands it is important to quantify also the time-frames needed for the achievement of irregular structure ([55]).

The goal of the present study was to examine and better understand the differences between the growth dynamics of Norway spruce and silver fir in the understory of forest stands undergoing conversion to selection forests. The particular objectives of the study were to (i) assess the tree species composition of regeneration and analyse the sapling density in relation to light conditions; (ii) quantify the effect of site, light conditions and neighbourhood characteristics on height growth and selected morphological attributes of Norway spruce and silver fir saplings; and (iii) quantify the time needed for regeneration to grow into the overstory.

# Materials and methods

## Study area

The study was conducted in two conifer forest stands (D1A and D2C) in an advanced phase of conversion to selection forests, growing on different sites. The study area is situated in the southwestern part of the Nízke Tatry Mountains (48° 52′ N, 19° 14′ E) in Central Slovakia. Since the 1980s, the stands have been part of the demonstration area Donovaly-Mistríky (50.3 ha), established by the Technical University in Zvolen, to follow the effects of silvicultural interventions on the structure and development of forest stands in various phases of conversion to selection forests ([40]). The investigated stands, covering 3.58 ha (D1A) and 2.50 ha (D2C), are located at an altitude of 900-950 m a.s.l. on a NW slope. The parent material is represented by granite, and the dominant soil types are shallow, acidic, nutrient-poor umbric leptosols ([53]). The mean annual temperature reaches 4.2-4.5 °C, and the mean annual precipitation is 1000-1100 mm. The D1A and D2C differ in slope (25% and 85%, respectively) and soil skeleton content (30% and 65% on average, respectively), resulting in more severe growth conditions and lower productive potential in the D2C. The site index (defined as the height of dominant trees at the age of 100 years) reaches 28 for the spruce and 26 for the fir trees in the D1A, and 24 for the spruce and 20 for the fir trees in the D2C. The ground vegetation is dominated by Vaccinium myrtillus L., with a cover of 36% in the D1A and 48% in the D2C ([53]). In the whole area, the forests are classified as Piceion excelsae ([36]). The dominant tree species in the investigated stands is Norway spruce. Before the start of the conversion process, the proportion of Norway spruce reached in the D1A 75% and in the D2C 99% of the growing stock (data from the local Forest Management Plan from 1978). The rest of the growing stock was largely dominated by silver fir, with a minor proportion (1-3%) of European larch (Larix decidua Mill.), Scots pine (Pinus sylvestris L.) and rowan (Sorbus aucuparia L.). European beech, once the natural admixture in fir-spruce forests, was eliminated from the stands due to charcoal production at the end of the 18th century. Intending to restore its presence, small groups of beech were underplanted in the stands in 1981-85. The groups were located mostly within existing gaps, and 3-year-old seedlings were used for the underplanting. The number of underplanted beech seedlings was 350 ha-1 in the D1A and 400 ha-1 in the D2C ([40]). The initial structure of the investigated stands before the start of the conversion process in the early 1980s was already partly differentiated. At the beginning of the study in 1993, the stands were already in an advanced stage of conversion. The last deliberate silvicultural intervention was applied in 1995. The selection cutting with intensities of 78 m3 ha-1 (D1A) and 32 m3 ha-1 (D2C) was mainly focused on the reduction of stems in the overstory. In the following period, the silvicultural interventions in both stands were limited to sanitation cuttings of windthrown Norway spruce, reaching an intensity of 53 m3 ha-1 in the D1A and 4 m3 ha-1 in the D2C (data according to local Reports on the Annual Harvest from 1993-2013). Regarding the damage by ungulates, the experimental setting was representative of standard forest management conditions with no special protection measures (e.g., fencing) and browsing controlled by hunting.

## Data collection

The diameter distribution data (overstory trees with DBH > 8 cm) originated from full callipering conducted in the studied stands in 1993 and 2013. Regeneration (trees with DBH ≤ 8 cm) was surveyed in 31 (D1A) and 20 (D2C) sample plots established in a regular grid of 35 × 35 m in 2013. Each sample plot consisted of two concentric circles (250 m2 and 100 m2) and a square (25 m2) located in the centre. In the central 5 × 5 m square, for the trees in the “seedlings” (height ≤ 20 cm, excluding one-year-old seedlings) and “younger saplings” (height 20-130 cm) categories, we recorded the tree species and height. In the circle of 100 m2, the trees classified as “older saplings” (height > 130 cm, DBH ≤ 4 cm) were sampled, and we recorded the tree species, DBH and height. In the circle of 250 m2, the trees in the category “undergrowth” (DBH 4-8 cm) were measured using the same parameters as those used for the older saplings. In all the regeneration categories, we also registered the presence of mechanical damage (broken terminal shoot) and damage by ungulates, i.e., browsing and bark stripping. For a detailed analysis of tree growth, we selected 65 Norway spruce and 65 silver fir trees from the older saplings category in each stand (260 trees in total). The saplings were chosen to cover the entire range of light and competition conditions and were randomly distributed throughout the investigated stands. Saplings with visible damage and located in the vicinity of recent cuttings were excluded from the analysis. If the saplings were growing in a cluster, we selected the highest tree. For this set of older saplings, we recorded the tree species, DBH, height, the length of the terminal and lateral shoots located on the upper three whorls and the length of the whole live crown from the top to the last live branch. The shoot lengths were measured using a telescopic measuring rod mEssfix-S® (Nedo, Dornstetten, Germany) with a range of 8 m. Subsequently, we calculated the following characteristics (sensu [50]): the mean annual height increment of the last 3 years, the apical dominance ratio (ratio between the lengths of the terminal and lateral shoot located on the third whorl), the relative growth ratio (mean annual height increment of the last 3 years/ tree height × 100) and the relative crown length (length of live crown/tree height × 100). Moreover, for each selected older sapling, the intensity of competition was quantified using a competition index ([18], [26]) that considered the competing trees inside a 4 m2 large circular area (1.13 m radius) surrounding the sapling. The area was divided into four quadrants, and the competition index (I) was calculated according to the following formula (eqn. 1):

$$I={\frac{1}{H_{OS}}} \cdot {\frac{1}{4}}\sum_{j=1}^{4} {\frac{C_{j}H_{j}} {D_{j}}}$$

where HOS is the height of a selected older sapling (in cm), Cj is the percent surface covered by the competing saplings in quadrant j, Hj is the mean height of the competing saplings in quadrant j (in cm), and Dj is the mean distance between the competing saplings in quadrant j and the selected older sapling (in cm).

In the centre of each sample plot, hemispherical photographs were taken at 1.3 m and 3.5 m above the ground (upper limits for categories of younger and older saplings) to assess the light environment. The amount of natural regeneration quantified by sapling density was subsequently analysed in relation to the light conditions. Additionally, hemispherical photographs were made above the terminal shoots of all 260 selected older saplings. We used a Nikon CoolPix® P5000 10 MP camera (Nikon Inc., Tokyo, Japan) with an 8 mm fish-eye lens converter FC-E8 and a self-levelling mount (Regent Instruments Inc., Sainte-Foy, Quebec, Canada) installed on a telescopic tripod Manfrotto with a maximum height of 7.3 m. Hemispherical photographs were taken from July-September 2013 under overcast skies or after sunset to avoid the influence of direct solar radiation. The acquired images were processed using the software WinSCANOPY 2006b (Regent Instruments Inc.) using an automatic pixel classification based on grey scale. The local light conditions were quantified by the ratio between the diffuse photosynthetic photon flux density under and over the canopy that resulted from the analysis of hemispherical photographs (“indirect site factor” in WinSCANOPY 2006b).

In addition, we also analysed the age structure and past growth of a sample of lower overstory stems (DBH = 8-12 cm). For 50 randomly selected trees (28 Norway spruces and 22 silver firs) with a DBH near the upper limit of the first diameter class (10-12 cm), we carried out a complete stem analysis. The trees were cut directly above the ground, their total height was recorded and 2 cm-thick cross-sections were taken at 0.5 m intervals along the stem, with the first cross-section taken at the stem base of the felled tree. In the laboratory, the samples were prepared according to standard methodology ([11]) and scanned at 1200 dpi using an Epson Expression 10000XL® scanner. The digitized cross-sections were subsequently analysed using WinDENDRO 2009b software (Regent Instruments Inc.), and the number of rings was counted for each disc. Based on the method of Carmean ([8]), we reconstructed the annual height growth of the analysed trees, which was subsequently used for the construction of the growth curves. The average annual height growth rate of the spruce and fir from the lower overstory was calculated according to the fitted species-specific age-height curves as the ratio of height and respective age.

## Statistical analyses

All of the statistical analyses were performed using the software package R ([38]) and STATISTICA® v.10.0 ([48]). A finite mixture of two Weibull functions (seven-parameter form - [58]) was fitted to the empirical diameter distributions of each measurement year and forest stand by the “mixdist” package of R software. The differences between the empirical diameter distributions and their equilibrium curves ([53]) were examined using a two-sample Kolmogorov-Smirnov test. For testing differences between the relative radiation in the D1A and D2C, we used a t-test for independent samples. The relationships between diffuse light and sapling density were analysed using simple linear regression ([57]).

Multiple linear regression was used to test the effect of multiple explanatory variables, including stand type, species, diffuse light and competition on the individual growth traits (y): tree height increment (both absolute and relative), apical dominance ratio and relative crown length. To quantify the competition status, two indices were tested: (i) the distance-dependent competition index I (compInd, eqn. 2); and (ii) the simple sum of the heights of the competing trees (sumHeight, eqn. 3). To quantify the effect of a variable within the levels of the other variables, we added linear interactions in the model. We developed the following model variants (eqn. 2, eqn. 3):

\eqalign {y_{ij} &= b_0 + b_1 \cdot standType_{ij} + b_2 \cdot species_{ij} \\ &+ b_3 \cdot compInd_{ij} + b_4 \cdot difLight_{ij} \\ &+ b_5 \cdot compInd_{ij} \cdot difLight_{ij} + \varepsilon}
\eqalign{y_{ij} &= b_0 + b_1 \cdot standType_{ij} + b_2 \cdot species_{ij} \\ &+ b_3 \cdot sumHeight_{ij} + b_4 \cdot difLight_{ij} \\ &+ b_5 \cdot sumHeight_{ij} \cdot difLight_{ij} + \varepsilon}

where y is the growth trait of interest, b0-b5 are regression coefficients to be estimated, the indices i and j denote the location in which a plot is placed and the observational plot itself, respectively, and ε is the error term.

However, we found a significant correlation between compInd and difLight (R2 = 0.17, p < 0.05) as well as between sumHeight and difLight (R2 = 0.11, p < 0.05). The variance inflation factor (calculated using the “vif” R function of the “car” package - [21]) indicated that multicollinearity can strongly affect the estimation of the regression coefficients (VIF > 4 for both predictors). Therefore, with the aim of keeping all the predictor variables in the models, we first removed the relationship among these variables. The residuals from the relationship between diffuse light (predictor) and the competition index/sum of the heights (dependent variable) were then used in the model to represent the pure competition effect (the residuals replaced compInd and sumHeight in eqn. 2 and eqn. 3). We used “ggplot2” in R ([56]) to estimate the effects of the predictors with the respective 95% confidence intervals.

The age-height curves of the spruce and fir from the lower overstory were fitted using the Chapman-Richards function (eqn. 4):

$$h=a(1-{e}^{-bt})^{c}$$

where h is the tree height, t is the age, and a, b, c are coefficients to be estimated.

For the parameterization of the growth function we used the nonlinear method of least squares, and for the optimization, we used the derivate Rosenbrock and quasi-Newton method in STATISTICA® v. 10.0.

# Results

## Stand history and natural regeneration

Between 1993 and 2013, the density of stems with a DBH > 8 cm decreased in both the studied stands, from 707 ha-1 to 633 ha-1 in the D1A and from 631 ha-1 to 622 ha-1 in the D2C (Tab. 1). At the same time, an increase in the basal area and growing stock was recorded. Regarding the tree species composition, we observed a substantial reduction in the spruce proportion, especially when considering the stem numbers (from 68% to 52% in the D1A and from 96% to 81% in the D2C). The decrease in spruce was compensated by the increased proportion of other tree species, especially fir in the lower diameter classes (from 25% to 39% in the D1A and from 3% to 17% in the D2C). Although not of the same magnitude, similar changes were observed when considering the basal area and volume.

Tab. 1 - Basic stand characteristics (trees of DBH > 8 cm) in investigated forests in 1993 and 2013. (*): European larch, Scots pine, rowan.

Characteristics Species D1A D2C
1993 2013 1993 2013
Stem density Total (N ha-1) 707 633 631 622
Spruce (%) 68 52 96 81
Fir (%) 25 39 3 17
Beech (%) 5 8 0 1
Others* (%) 2 1 1 1
Basal area Total (m2 ha-1) 45 47 33 39
Spruce (%) 72 61 97 87
Fir (%) 24 33 3 13
Beech (%) 1 2 0 0
Others* (%) 3 4 0 0
Growing stock Total (m3 ha-1) 495 519 363 438
Spruce (%) 71 62 97 87
Fir (%) 26 33 3 13
Beech (%) 1 1 0 0
Others* (%) 2 4 0 0

The diameter distributions of the investigated stands in 1993 as well as those in 2013 approached S-shaped curves (Fig. 1). The differences between the empirical and optimal (reverse J-shaped) distributions in both the measurement years and study sites remained significant (p < 0.05). During the study period, we recorded moderately increased stem numbers in diameter classes ≥ 38 cm and underrepresented numbers of trees in the lower diameter classes.

Fig. 1 - Equilibrium curves and diameter distributions recorded in the D1A (left) and the D2C (right).

The overall density of the regeneration (i.e., all the individuals, from seedlings to undergrowth, with a DBH ≤ 8 cm) in the studied stands was 21.461 ± 2.380 ha-1 in the D1A and 15.016 ± 2.697 ha-1 in the D2C. In the regeneration, four tree species were present: naturally regenerated spruce, fir, rowan as well as underplanted beech (Tab. 2). In both stands and in all but one regeneration category (younger saplings in the D2C), the proportion of fir was significantly higher than in the overstory. Regarding the regeneration categories, the share of spruce in the seedlings and younger saplings categories was higher than in the older saplings and undergrowth categories. An opposite trend of increasing proportion with regeneration age was observed for the fir, which showed higher proportions in the categories of older saplings and undergrowth. The most frequently damaged tree species in both the investigated stands was fir, with the damage concentrated above all to the categories of older saplings and undergrowth. However, the damage in any category of regeneration did not exceed 15%. The dominant type of damage was caused by ungulates (browsing and bark stripping). The regeneration of spruce, beech and rowan was damaged only sporadically (< 1%).

Tab. 2 - Density and species composition of regeneration (mean ± standard error) in investigated forests in 2013.

Stand Species Seedlings Younger saplings Older saplings Undergrowth
N ha-1 % N ha-1 % N ha-1 % N ha-1 %
D1A Spruce 6338 ± 1623 51 2723 ± 844 41 457 ± 102 25 138 ± 29 22
Fir 5476 ± 801 44 3246 ± 387 49 1254 ± 231 69 451 ± 82 70
Beech - - 62 ± 36 1 112 ± 29 6 51 ± 17 8
Rowan 538 ± 203 5 615 ± 126 9 - - - -
Total 12352 ± 1812 100 6646 ± 943 100 1823 ± 235 100 640 ± 79 100
D2C Spruce 4000 ± 965 67 6080 ± 1655 86 810 ± 165 53 250 ± 47 57
Fir 1760 ± 905 29 700 ± 139 10 515 ± 117 34 174 ± 43 40
Beech - - 40 ± 40 1 195 ± 103 13 12 ± 8 3
Rowan 260 ± 113 4 220 ± 93 3 - - - -
Total 6020 ± 1360 100 7040 ± 1729 100 1520 ± 205 100 436 ± 59 100

The light microclimate in the sample plots did not differ significantly between the investigated sites (t = 0.581, p = 0.568 in height 1.3 m and t = 0.447, p = 0.660 in height 3.5 m). The median values of the relative diffuse light reached 13.2% (D1A) and 12.0% (D2C) at 1.3 m height and 14.7% (D1A) and 13.0% (D2C) at 3.5 m height. For Norway spruce, significant correlations between the relative diffuse light and sapling density were found in both categories (younger saplings: R = 0.47, p = 0.008; older saplings: R = 0.57, p = 0.001), but only in the D1A (Fig. S1 in Supplementary material). In silver fir, no significant relationship between relative light and the sapling density was revealed in any category or site.

## Growth of older saplings

The intensity of relative diffuse light measured above the terminal shoot of the analysed spruces and firs ranged from 7.2 to 26.7% in the studied stands. Under these light conditions, the values of the relative growth ratio varied between 0.3% and 10.7% for the spruce and between 0.2% and 13.2% for the fir (Tab. 3, Fig. S2 in Supplementary material).

Tab. 3 - Basic characteristics of older saplings (mean ± standard error; the range is given in parentheses).

Parameter D1A D2C
Spruce Fir Spruce Fir
Height
(cm)
189 ± 6 210 ± 6 186 ± 6 206 ± 6
(131-308) (130-324) (131-311) (130-310)
Mean height increment
(cm year-1)
8.9 ± 0.7 12.9 ± 0.9 5.7 ± 0.5 8.6 ± 0.6
(0.7-30.0) (1.0-30.5) (0.8-16.5) (0.5-21.2)
Relative growth ratio
(%)
4.6 ± 0.4 6.1 ± 0.4 3.1 ± 0.3 3.8 ± 0.3
(0.3-10.7) (0.8-13.2) (0.5-10.3) (0.2-10.6)
Apical dominance ratio 0.9 ± 0.1 0.9 ± 0.1 0.7 ± 0.1 0.7 ± 0.1
(0.1-1.8) (0.3-1.6) (0.3-1.5) (0.1-1.6)
Relative crown length
(%)
68.7 ± 1.5 74.5 ± 1.8 58.6 ± 2.5 70.7 ± 1.8
(27.3-90.1) (20.7-95.2) (17.2-88.2) (40.1-93.3)
Competition index 0.5 ± 0.1 0.4 ± 0.1 0.5 ± 0.1 0.4 ± 0.1
(0.0-0.9) (0.0-0.7) (0.1-0.8) (0.0-0.8)

The competition index (compInd) performed better than the sum of the heights of the competitors (sumHeight) regarding the explained variability in the growth traits (Tab. S1). Therefore, in the following analysis, we interpreted the results of the first variant (eqn. 2). The proportion of the variability explained by the predictors ranged from 36% (apical dominance ratio) to 52% (height increment). The diffuse light and stand type variables most significantly affected the growth traits tested in the study (p < 0.001 - Fig. 2). However, diffuse light explained most of the variability in all the growth traits, ranging from 22% (apical dominance ratio) to 34% (height increment). The competition status showed a negative correlation, but only for the relative crown length (proportion of the variability explained was 7.4%, p < 0.001). The interaction between competition and diffuse light was only significant for the relative crown length and showed a positive effect of increased diffuse light at the same competition level (2.3% of the total variability, p < 0.001). Except for the apical dominance ratio, we found that the tree species significantly contributed to the explanation of the total variability in the growth traits (5-8% of the total variability, p < 0.001).

Fig. 2 - Estimation of coefficients of the multiple linear regression model (fitted by eqn. 2). The coefficients express the effects of predictors on the individual growth traits (height increment, relative growth ratio, apical dominance ratio and relative crown length). The inner graph zooms in the predictors that have the coefficients values close to zero. The whiskers denote 95% confidence intervals. (standType): stand type; (diffLight): diffuse light; (compInd): competition index; (*): p < 0.05; (**): p < 0.01; (***): p < 0.001; (ns): not significant.

## Age structure and height growth of the lower overstory trees

The age of the analysed spruce and fir trees (DBH 10-12 cm) varied between 40 and 99 years (median value 52 years) in the D1A and between 32 and 107 years (median value 47 years) in the D2C. The majority of the trees (83% in the D1A and 88% in the D2C) reached the upper limit of the first diameter class (i.e., 12 cm) at the age of 30-60 years (hereafter “normal-growing trees”). However, among the analysed samples, we found several trees that reached the given threshold at ages distinctively exceeding 60 years (hereafter “slow-growing trees” - Tab. 4, Fig. 3; Fig. S3 in Supplementary material). The slow-growing trees were almost exclusively spruce individuals, and only one slow-growing fir was found in the D1A. The average heights of the slow-growing trees reaching the upper limit of the first diameter class were 25-35% lower than those of the normal-growing trees.

Tab. 4 - Parameters of age-height curves fitted with the Chapman-Richards function.

Density Species No a b c R
D1A Normal-growing spruce 9 13.329 -0.070 5.555 0.954
Slow-growing spruce 3 11.563 -0.032 4.928 0.955
Normal-growing fir 11 14.354 -0.045 3.477 0.969
Slow-growing fir 1 127.990 -0.001 1.078 0.992
D2C Normal-growing spruce 13 12.272 -0.059 3.636 0.906
Slow-growing spruce 3 19.223 -0.011 2.169 0.862
Normal-growing fir 10 19.763 -0.040 3.717 0.942

Fig. 3 - Growth curves of individual lower overstory trees (DBH 10-12 cm) with fitted species-specific age-height curves

For the calculation of the average annual height increment, we used the values of the fitted heights at the age of 60 years for both the normal- and slow-growing trees. The height growth rate of the normal-growing trees reached 18.4 cm in the D2C and 20.5 cm in the D1A per year for the spruce and 18.7 to 23.1 cm per year for the fir. Regarding the slow-growing trees, the rates were considerably lower, with a range of 6.3-9.0 cm per year.

# Discussion

## Stand structure, natural regeneration and effect of light on sapling density

The results of this study indicate that the investigated stands still show certain deviations from the optimal diameter distributions and tree species composition (e.g., overrepresented spruce, absence of natural regeneration of beech, etc.). Low intensity of silvicultural interventions in the last two decades resulted in an increase in stem density and volume in the upper diameter classes and a decrease in the lower diameter classes of the stands (Tab. 1, Fig. 1). The progress of the ongoing conversion process was slowed, and therefore, the achievement of the optimal selection structure was delayed as well. In both of the studied stands, the share of the trees with a DBH > 36 cm in the growing stock (80% in the D1A and 67% in the D2C) was considerably higher compared with the theoretical structure (61% and 54%). As a consequence, the current growing stock significantly exceeds the optimal levels (410 m3 ha-1 for the D1A and 320 m3 ha-1 for the D2C), by 26% (D1A) or 37% (D2C). Nevertheless, due to underrepresented numbers of trees in the lower diameter classes, the amount of transmitted solar radiation (median values 13.2% in the D1A and 12% in the D2C) was higher in comparison with similar studies performed in spruce-fir (beech) selection forests with diffuse light levels between 3.9 and 6.4% ([3], [15]). This structure allowed the advanced natural regeneration (i.e., older saplings and undergrowth) to survive in large numbers, which exceeds the threshold value of 600 ha-1 reported by Duc ([17]), the value necessary for the maintenance of selection structure in spruce-fir forests. Regarding the future development of the two analysed conversion stands, the regular application of silvicultural interventions that gradually reduces the density of large-diameter stems can increase the transition rate between the lower and upper overstory and significantly improve the differentiation of stand structure ([42], [7]).

In the investigated stands, we observed a progressive change in the representation of the main tree species, spruce and fir. During the previous two decades, a visible increase in the fir proportion at the expense of spruce was observed in the parent stand, especially in the lower overstory. Moreover, fir became established also in the understory (Tab. 1, Tab. 2). The main reason for the increasing proportion of fir was probably the low intensity of diffuse light in the stand. Shaded conditions do not allow spruce to realize its high growth potential, and under such conditions, fir is a stronger competitor ([24], [47], [12]). Another important factor to be considered is the impact of ungulates, as the overbrowsing of regeneration was identified by several studies ([30], [16]) as one of the main reasons for fir decline in Europe. In the stands surveyed in our study, the low density of ungulates (personal observation), causing only insignificant damage to fir regeneration, was very likely another reason for the current successful development of fir regeneration.

The natural regeneration of fir prefers mostly microsites under stand canopies or at the edges of canopy gaps ([25], [14], [34]). Here, fir survival is limited by light, with the minimal required values from 2% for seedlings to 10% for 5- to 15-year-old saplings ([23], [28]). The lower ability of fir saplings to utilize high-light conditions compared to spruce was confirmed also by Grassi & Bagnaresi ([24]). The results of our study, similar to those of Cescatti ([10]), indicate that the density of fir saplings was not significantly affected by sunlight. The same situation was observed for spruce in the D2C, while the only positive correlation between diffuse light and regeneration density was confirmed for spruce in the D1A (Fig. S1 in Supplementary material). We assume that in the D2C, the influence of diffuse light was confounded by less favourable soil conditions ([53]) that can play a decisive role in the regeneration processes of spruce as well as fir ([14], [34]). Another major factor could be the competition of ground vegetation ([20], [51]), represented by large patches of Vaccinium myrtillus in the D2C ([53]). However, the dynamics of regeneration processes can be affected also by other factors, such as seed dispersal ([51]), micro-relief ([14]), seedbed characteristics ([46]) or the tree species composition of the canopy layer ([34], [9]).

## Impact of site conditions and neighbourhood characteristics on saplings’ growth

Recently, the analysis of the height growth of spruce and fir natural regeneration in uneven-aged, mixed forests was the object of several studies ([50], [47], [15]). A direct comparison of the results was nevertheless problematic because there were considerable differences in the height range of the surveyed samples. In general, the relative height increments reported for the natural regeneration of spruce and fir in the mentioned studies varied between 5% and 11%. In our study, comparable results were observed for the older saplings in the D1A; however, in the less favourable site conditions of the D2C, the relative height increment along the examined light gradient was distinctively lower (Tab. 3). Similar to our results (Fig. 2), the positive impact of diffuse light on the height growth of spruce as well as fir with heights < 2.2 m was confirmed also by the studies of Filipiak et al. ([20]) and Grassi & Giannini ([26]). On the contrary, in the study of Diaci & Firm ([15]) no correlation between height increment and diffuse light was revealed for spruce and fir seedlings. The authors observed a significant relationship only in the sapling category (height 20-130 cm), with significantly larger height increments of fir. Under average diffuse light levels of 6-8%, the annual height increment of fir exceeded that of spruce by approximately 1 cm, corresponding to the results of our study (Fig. 2 and Fig. S2 in Supplementary material). Similar results were confirmed by Stancioiu & O’Hara ([47]) from mixed mountain forests in Romania. In low-light environments (relative irradiance < 35%), the height increments of the fir and beech saplings reached significantly higher values than those of spruce. In contrast, in intermediate light conditions (relative irradiance 35-70%), the authors registered comparable height growth of all the tree species and in open conditions (relative irradiance > 70%), even the tendency of spruce to outgrow beech and fir.

The ability of saplings to effectively use solar radiation is crucially affected by the morphologies of their assimilatory apparatuses ([33]). Therefore, the apical dominance ratio and the relative crown length of conifer saplings could be regarded as simple indicators of light conditions ([35], [39], [50]), especially for spruce, as was documented by our results (Fig. S2). In general, the lower the apical dominance ratio, the more seedling growth is affected by light deficiency ([26]). In our study, the threshold value of 1 was reached at relative radiation levels of 18% (fir) and 19% (spruce) in the D1A; however, in the D2C, the threshold value was reached only after the relative light exceeded 22% (fir) and 26% (spruce - Fig. S2 in Supplementary material). This indicates that the more severe a site is, the higher are the saplings’ demands for light. For silviculture in sites with less favourable growth conditions, it is therefore recommended to adequately reduce the period of shading ([6]). Under such conditions, in the stands undergoing conversion, it is useful to reduce the canopy density more strongly and keep the growing stock at a lower level. Nevertheless, a too strong intervention can result in the establishment of spatially continuous, even-aged natural regeneration that, regarding stand differentiation, should be avoided ([42], [40]). As registered in our study, the spruce and fir reacted to changes in the light microclimate, showing differences in height growth and relative crown length. The observed differences between these two species in height growth can most likely be attributed to the physiological characteristics of the needles, with better adaptation of fir to low light environments. The differences could also be attributed to the better ability of fir saplings to adjust their crowns to changing light conditions ([24], [39]).

Along the examined light gradient (relative radiation 7-27%), we observed a negative influence of local sapling density (competition index) only on the relative crown length of the older saplings of spruce and fir (Fig. 2). In contrast, no significant relationship was confirmed between local sapling density and height growth or the apical dominance ratio. Several studies ([29], [26]) have suggested that in shaded conditions, light has by far the highest influence on tree growth, and light usually overlay the impact of the other factors including local sapling density. Intraspecific competition was confirmed to affect the height increment only in non-limiting light conditions ([35], [18]). We therefore assume, that lower values of diffuse light registered in the conversion stands in our study were the main reason why the local sapling density did not play any important role in the height growth of older saplings. Nevertheless, the generalization of these results is limited by the fact that we analysed only the height growth of the highest saplings in the clusters.

## Long-term growth dynamics of the lower overstory trees

Continuous ingrowth of natural regeneration into the smallest diameter classes is the basic feature of a system’s self-regulation that, in conversion stands, begins to function only after a relatively long time period. Based on the analysis of 100 years of conversion from regular afforestation to currently differentiated, spruce-dominated forest stands in Swiss Jura, Schütz ([41]) concluded that the formation of irregular vertical structures requires at least 60-70 years. This estimation more or less corresponds with the findings from our stands, where the majority of the young spruces and firs (“normal-growing trees”) reached the upper limit of the first diameter class at approximately 50 years, with an average height growth rate of 18-23 cm per year (Fig. 3 and Fig. S3 in Supplementary material). Despite the different site conditions, the lower growing stock and thus the more favourable light microclimate for the growth of the regeneration and lower overstory in the D2C resulted in the similar ages of the normal-growing trees, reaching the DBH threshold in both the investigated stands. In particular, the light-activated trees that are continuously growing up to the higher stand layers are critical for the maintenance of the equilibrium state of a selection forest ([41], [49]). In addition to the trees of average growth, in the studied stands were rarely present also slow-growing spruce and fir trees with average height increment less than 10 cm (mainly in the parts with complicated access and therefore a lower cutting intensity). Such trees can survive under the stand canopy for a very long period (in special cases more than 200 years) and form the potential reserve for the transitions into higher tree layers. After the release, they are able to reach dimensions comparable to the normal-growing trees regardless of their age ([31]).

Considering the decreasing height increment rate during the previous 20 years (Fig. 3), we assume that a more regular application of silvicultural interventions in both of the investigated stands during this period, especially a more substantial reduction in large overstory stems, would considerably promote the growth of fir and especially spruce regeneration and ingrowth to the overstory. As confirmed by several studies ([40], [9], [7]), stand structures with excessive growing stock negatively affect the density of natural regeneration, the transition of saplings into the lower overstory and the further transition into the higher stand layers. If the stand density exceeds a certain threshold, regeneration is no longer able to provide enough saplings for continuous height transitions, and the selection structure will be progressively lost ([43]). Moreover, low light levels can also influence the tree species composition of the understory, as this represents a competitive disadvantage for spruce against fir ([12]). In conversion stands, where forest managers usually have to address this kind of problem, the recommended solution could be the regular application of classical differentiation thinning (sensu [42]) focused on the refinement of the structure and improvements towards the equilibrium state of a selection forest.

# Conclusions

Our results showed that the formation of lower overstory in stands undergoing the conversion to single-tree selection requires a time period of about 50 years. In these conditions, height growth of fir and spruce natural regeneration was significantly influenced by the intensity of diffuse radiation, but not the local sapling density. Natural regeneration of fir showed better dynamics of height growth than spruce. Considering the current negative impacts of climate change, especially on spruce, lower growing stocks in the stands could provide more suitable conditions for spruce regeneration and therefore contribute to its maintenance in spruce-fir mixtures. Results on conversion of monospecific forest stands to mixed diverse forest structures are an important step toward an understanding and cultivating stable forest ecosystems resilient to disturbances.

# Authors’ contributions statement

The study was conducted with contributions from all authors. JV, MS and SK designed and planned the study. JV, DS, JP and ZP conducted the field measurements and laboratory work. MB and JV performed the statistical analyses. JV, SK and PJ analyzed and interpreted data. JV prepared the draft of the manuscript.

# Acknowledgements

This research was funded by the Ministry of Education, Science, Research and Sports of the Slovak Republic - Scientific Grant Agency (VEGA), projects VEGA 1/0492/17 and VEGA 1/0021/18.

# Conflict of interest

The authors declare that they have no conflict of interest.

# References

(1)
Augusto L, Ranger J, Binkley D, Rothe A (2002). Impact of several common tree species of European temperate forests on soil fertility. Annals of Forest Science 59: 233-253.
CrossRef | Gscholar
(2)
Boden S, Kahle HP, Wilpert K, Spiecker H (2014). Resilience of Norway spruce (Picea abies (L.) Karst) growth to changing climatic conditions in Southwest Germany. Forest Ecology and Management 315: 12-21.
CrossRef | Gscholar
(3)
Boncina A, Diaci J, Cencic L (2002). Comparison of the two main types of selection forests in Slovenia: distribution, site conditions, stand structure, regeneration and management. Forestry 75: 365-373.
CrossRef | Gscholar
(4)
Bošela M, Sedmák R, Sedmáková D, Marušák R, Kulla L (2014). Temporal shifts of climate-growth relationships of Norway spruce as an indicator of health decline in the Beskids, Slovakia. Forest Ecology and Management 325: 108-117.
CrossRef | Gscholar
(5)
Brang P, Spathelf P, Larsen JB, Bauhus J, Bončina A, Chauvin C, Drössler L, Güemes CG, Heiri C, Kerr G, Lexer MJ, Mason B, Mohren F, Mühlethaler U, Nocentini S, Svoboda M (2014). Suitability of close-to-nature silviculture for adapting temperate European forests to climate change. Forestry 87: 492-503.
CrossRef | Gscholar
(6)
Burschel P, Huss J (1997). Grundriß des Waldbaus [Silviculture]. Parey Buchverlag, Berlin, Germany, pp. 487. [in German]
Gscholar
(7)
Cameron AD, Hands MOR (2010). Developing a sustainable irregular structure: an evaluation of three inventories at 6-year intervals in an irregular mixed-species stand in Scotland. Forestry 83: 469-475.
CrossRef | Gscholar
(8)
Carmean WH (1972). Site index curves for upland oaks in the Central States. Forest Science 18: 102-120.
Online | Gscholar
(9)
Cavlović J, Bozić M, Boncina A (2006). Stand structure of an uneven-aged fir-beech forest with an irregular diameter structure: modeling the development of the Belevine forest, Croatia. European Journal of Forest Research 125: 325-333.
CrossRef | Gscholar
(10)
Cescatti A (1996). Selective cutting radiative regime and natural regeneration in a mixed coniferous forest. In: “Modelling Regeneration Success and Early Growth of Forest Stands” (Skovsgaard JP, Johannsen P eds). Danish Forest and Landscape Research Institute, Horsholm, Denmark, pp. 474-483.
Gscholar
(11)
Cook ER, Kairiukstis LA (1990). Methods of dendrochronology: applications in the environmental sciences. Kluwer Academic Publishers, Dordrecht, Netherlands, pp. 394.
Gscholar
(12)
Danescu A, Kohnle U, Bauhus J, Weiskittel A, Albrecht AT (2018). Long-term development of natural regeneration in irregular, mixed stands of silver fir and Norway spruce. Forest Ecology and Management 430: 105-116.
CrossRef | Gscholar
(13)
Dengler A (1980). Waldbau [Silviculture]. Verlag Paul Parey, Hamburg, Berlin, Germany, pp. 350. [in German]
Gscholar
(14)
Diaci J, Pisek R, Boncina A (2005). Regeneration in experimental gaps of subalpine Picea abies forest in the Slovenian Alps. European Journal of Forest Research 124: 29-36.
CrossRef | Gscholar
(15)
Diaci J, Firm D (2011). Long-term dynamics of a mixed conifer stand in Slovenia managed with a farmer selection system. Forest Ecology and Management 262: 931-939.
CrossRef | Gscholar
(16)
Diaci J, Rozenbergar D, Anic I, Mikac S, Saniga M, Kucbel S, Visnjic C, Ballian D (2011). Structural dynamics and synchronous silver fir decline in mixed old-growth mountain forests in Eastern and Southeastern Europe. Forestry 84: 479-491.
CrossRef | Gscholar
(17)
Duc P (1991). Investigations on recruitment dynamics in Emmental selection forest plots. Schweizerische Zeitschrift für Forstwesen 142: 299-319. [in German]
Gscholar
(18)
Duchesneau R, Lesage I, Messier C, Morin H (2001). Effects of light and intraspecific competition on growth and crown morphology of two size classes of understory balsam fir saplings. Forest Ecology and Management 140: 215-225.
CrossRef | Gscholar
(19)
Dvorak L, Bachmann P, Mandallaz D (2001). Sturmschäden in ungleichförmigen Beständen [Storm damage in irregular stands]. Schweizerische Zeitschrift für Forstwesen 152: 445-452. [in German]
CrossRef | Gscholar
(20)
Filipiak M, Iszkulo G, Korybo J (2005). Relation between photosynthetic photon flux density (PPFD) and growth of silver fir (Abies alba Mill.) seedlings in a forest stand dominated by spruce [Picea abies (L.) Karst.] in the Sudety Mts (SW Poland). Polish Journal of Ecology 53: 177-184.
Online | Gscholar
(21)
Fox J, Weisberg S (2011). An R companion to applied regresion (2nd edn). Sage, Thousand Oaks, CA, USA, pp. 608.
Online | Gscholar
(22)
Gamborg C, Larsen JB (2003). “Back to nature” - a sustainable future for forestry? Forest Ecology and Management 179: 559-571.
CrossRef | Gscholar
(23)
Givnish TJ (1988). Adaptation to sun and shade: a whole-plant perspective. Australian Journal of Plant Physiology 15: 63-92.
CrossRef | Gscholar
(24)
Grassi G, Bagnaresi U (2001). Foliar morphological and physiological plasticity in Picea abies and Abies alba saplings along a natural light gradient. Tree Physiology 21: 959-967.
CrossRef | Gscholar
(25)
Grassi G, Minotta G, Tonon G, Bagnaresi U (2004). Dynamics of Norway spruce and silver fir natural regeneration in a mixed stand under uneven-aged management. Canadian Journal of Forest Research 34: 141-149.
CrossRef | Gscholar
(26)
Grassi G, Giannini R (2005). Influence of light and competition on crown and shoot morphological parameters of Norway spruce and silver fir saplings. Annals of Forest Science 62: 269-274.
CrossRef | Gscholar
(27)
Heuze P, Schnitzler A, Klein F (2005). Is browsing the major factor of silver fir decline in the Vosges Mountains of France? Forest Ecology and Management 430: 219-228.
CrossRef | Gscholar
(28)
Jaworski A (1995). Silvicultural characteristic of forest trees. Gutenberg, Kraków, Poland, pp. 237. [in Polish]
Gscholar
(29)
Klinka K, Wang Q, Kayahara GJ, Carter RE, Blackwell BA (1992). Light-growth response relationships in Pacific silver fir (Abies amabilis) and subalpine fir (Abies lasiocarpa). Canadian Journal of Forest Research 70: 1919-1930.
CrossRef | Gscholar
(30)
Klopcic M, Jerina K, Boncina A (2010). Long-term changes of structure and tree species composition in Dinaric uneven-aged forests: are red deer an important factor? European Journal of Forest Research 129: 277-288.
CrossRef | Gscholar
(31)
Lin CJ, Laiho O, Lähde E (2012). Norway spruce (Picea abies L.) regeneration and growth of understory trees under single-tree selection silviculture in Finland. European Journal of Forest Research 131: 683-691.
CrossRef | Gscholar
(32)
Mezei P, Grodzki W, Blaenec M, Jakuš R (2014). Factors influencing the wind-bark beetles’ disturbance system in the course of an Ips typographus outbreak in the Tatra Mountains. Forest Ecology and Management 312: 67-77.
CrossRef | Gscholar
(33)
Oliver CD, Larson BC (1996). Forest stand dynamics. John Wiley, New York, USA, pp. 520.
Gscholar
(34)
Paluch J (2005). Spatial distribution of regeneration in West-Carpathian uneven-aged silver fir forests. European Journal of Forest Research 124: 47-54.
CrossRef | Gscholar
(35)
Parent S, Messier C (1995). Effets d’un gradient de lumière sur la croissance en hauteur et la morphologie de la cime du sapin baumier régénéré naturellement [Effects of a light gradient on height growth and the morphology of the crown of naturally regenerated balsam fir]. Canadian Journal of Forest Research 25: 878-885. [in French]
CrossRef | Gscholar
(36)
Pawlowski B, Sokolowski M, Wallisch K (1928). Die Pflanzenassoziationen des Tatra-Gebirges. Teil 7. Die Pflanzenassoziationen und die Flora des Morskie Oko-Tales [Plant associations of the Tatra Mountains. Part 7. Plant associations and flora of the Morskie Oko Valley]. Extrait de Bullettin de l’Academie Polonaise des Sciences et des Lettres 2: 205-272. [in German]
Gscholar
(37)
Prpić B, Seletković Z (2001). Ecological constitution of the silver fir. In: “Silver Fir in Croatia” (Matić S, Vukelić J, Prpić B, Vidaković M, Klepac D, Glavaš M, Ljuljka B, Figurić M, Jakovac H eds). Tiskara Meić, Zagreb, Croatia, pp. 269-276.
Gscholar
(38)
R Development Core Team (2018). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
Online | Gscholar
(39)
Robakowski P, Wyka T, Samardakiewicz S, Kierzkowski D (2004). Growth, photosynthesis, and needle structure of silver fir (Abies alba Mill.) seedlings under different canopies. Forest Ecology and Management 201: 211-227.
CrossRef | Gscholar
(40)
Saniga M, Vencurik J (2007). Dynamics of structure and regeneration processes of the forests in various phases of the conversion to a selection forest in forest management unit Korytnica. Vydavatelstvo Technickej univerzity, Zvolen, Slovakia, pp. 82. [in Slovak]
Gscholar
(41)
Schütz JP (2001a). Opportunities and strategies of transforming regular forests to irregular forests. Forest Ecology and Management 151: 87-94.
CrossRef | Gscholar
(42)
Schütz JP (2001b). Der Plenterwald und weitere Formen strukturierter und gemischter Wälder [Selection forest and other forms of structured and mixed forests]. Parey Buchverlag, Berlin, Germany, pp. 220. [in German]
Gscholar
(43)
Schütz JP (2006). Modelling the demographic sustainability of pure beech plenter forests in Eastern Germany. Annals of Forest Science 63: 93-100.
CrossRef | Gscholar
(44)
Schütz JP, Götz M, Schmid W, Mandallaz D (2006). Vulnerability of spruce (Picea abies) and beech (Fagus sylvatica) forest stands to storms and consequences for silviculture. European Journal of Forest Research 125: 291-302.
CrossRef | Gscholar
(45)
Seidl R, Thom D, Kautz M, Martin-Benito D, Peltoniemi M, Vacchiano G, Wild J, Ascoli D, Petr M, Honkaniemi J, Lexer MJ, Trotsiuk V, Mairota P, Svoboda M, Fabrika M, Nagel TA, Reyer CPO (2017). Forest disturbances under climate change. Nature Climate Change 7: 395-402.
CrossRef | Gscholar
(46)
Simard MJ, Bergeron Y, Sirois L (1998). Conifer seedling recruitment in a southwestern Canadian boreal forests: the importance of substrate. Journal of Vegetation Science 9: 575-582.
CrossRef | Gscholar
(47)
Stancioiu PT, O’Hara KL (2006). Regeneration growth in different light environments of mixed species, multiaged, mountainous forests of Romania. European Journal of Forest Research 125: 151-162.
CrossRef | Gscholar
(48)
StatSoft Inc. (2010). Electronic statistics textbook. Tulsa, OK, USA.
Online | Gscholar
(49)
Sterba H, Zingg A (2001). Target diameter harvesting - a strategy to convert even-aged forests. Forest Ecology and Management 151: 95-105.
CrossRef | Gscholar
(50)
Szymura TH (2005). Silver fir saplings bank in seminatural stand: individuals architecture and vitality. Forest Ecology and Management 212: 101-108.
CrossRef | Gscholar
(51)
Szymura TH, Dunajski A, Aman I, Makowski M, Szymura M (2007). The spatial pattern and microsites requirements of Abies alba natural regeneration in the Karkonosze Mountains. Dendrobiology 58: 51-57.
Online | Gscholar
(52)
Van Der Maaten-Theunissen M, Kahle HP, Van Der Maaten E (2013). Drought sensitivity of Norway spruce is higher than that of silver fir along an altitudinal gradient in southwestern Germany. Annals of Forest Science 70: 185-193.
CrossRef | Gscholar
(53)
Vencurik J, Kucbel S, Saniga M, Jaloviar P (2016). Structure and height growth of natural regeneration in selection forests with different Norway spruce proportion in Slovakia. Technical university, Zvolen, Slovakia, pp. 64. [in Slovak]
Gscholar
(54)
Vitali V, Büntgen U, Bauhus J (2018). Seasonality matters - The effects of past and projected seasonal climate change on the growth of native and exotic conifer species in Central Europe. Dendrochronologia 48: 1-9.
CrossRef | Gscholar
(55)
Von Lüpke B, Ammer C, Bruciamacchie M, Brunner A, Ceitel J, Collet C, Deuleuze C, Placido JD, Huss J, Jankovič J, Kantor P, Larsen JB, Lexer M, Löf M, Longauer R, Madsen P, Modrzynski J, Mosandl R, Pampe A, Pommering A, Tesar V, Thompson R, Zientarski J (2004). Silvicultural strategies for conversion. In: “Norway spruce conversion: options and consequences” (Spiecker H, Hansen J, Klimo E, Skovsgaard JP, Sterba H, von Teuffel K eds). Brill, Leiden, Netherlands, pp. 121-164.
Gscholar
(56)
Wickham H (2016). ggplot2: Elegant graphics for data analysis. Springer Verlag, New York, USA, pp. 276.
Online | Gscholar
(57)
Zar JH (1999). Biostatistical analysis. Prentice Hall, New Jersey, USA, pp. 663.
Gscholar
(58)
Zhang L, Gove JH, Liu C, Leak WB (2001). A finite mixture of two Weibull distributions for modeling the diameter distributions of rotated-sigmoid, uneven-aged stands. Canadian Journal of Forest Research 31 (9): 1654-1659.
CrossRef | Gscholar

#### Authors’ Affiliation

(1)
Jaroslav Vencurik
Stanislav Kucbel
Milan Saniga
Peter Jaloviar
Denisa Sedmáková 0000-0001-5556-0839
Ján Pittner
Zuzana Parobeková
Department of Silviculture, Faculty of Forestry, Technical University in Zvolen, Masarykova 24, 960 53 Zvolen (Slovak Republic)
(2)
Michal Bosela 0000-0001-6706-8614
Department of Forest Management and Geodesy, Faculty of Forestry, Technical University in Zvolen, Masarykova 24, 960 53 Zvolen (Slovak Republic)

#### Corresponding author

Jaroslav Vencurik
vencurik@tuzvo.sk

#### Citation

Vencurik J, Kucbel S, Saniga M, Jaloviar P, Sedmáková D, Pittner J, Parobeková Z, Bosela M (2020). Growth dynamics of the Norway spruce and silver fir understory in continuous cover forestry. iForest 13: 56-64. - doi: 10.3832/ifor3183-012

#### Paper history

Accepted: Nov 28, 2019

First online: Feb 05, 2020
Publication Date: Feb 29, 2020
Publication Time: 2.30 months

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

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

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

Breakdown by View Type
HTML Page Views: 6494
Abstract Page Views: 479

Web Metrics
Days since publication: 443
Overall contacts: 7955
Avg. contacts per week: 125.70

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

(No citations were found up to date. Please come back later)

#### iForest Database Search

Search By Author

Search By Keyword

Citing Articles

Search By Author

Search By Keywords

#### PubMed Search

Search By Author

Search By Keyword