## Development phases and structural characteristics of the Penteleu-Viforâta virgin forest in the Curvature Carpathians

iForest - Biogeosciences and Forestry, Volume 13, Issue 5, Pages 389-395 (2020)
doi: https://doi.org/10.3832/ifor3094-013

Research Articles

The stand structure of a virgin forest situated at an average altitude of 1130 m a.s.l. in the Milea Viforâta Nature Reserve (Southern Carpathians, Romania) was investigated to determine the specific development phases of the forest and understand how they influence the stand structure, with the aim of providing optimal solutions and structural models for sustainable forest management. All trees with breast height diameter (dbh) ≥ 8 cm were inventoried in the study plot (1 ha), and the main dendrometrical variables were measured. Radial increment cores were taken from all the trees and were subsequently processed. A total of 317 trees from three species - European beech (Fagus sylvatica), silver fir (Abies alba) and Norway spruce (Picea abies) - were sampled at different development phases (optimum, ageing, breakdown and dieback, rejuvenation). Testing stand structural diversity with the Gini index, a minimal stability was found in the rejuvenation development phase and a maximum stability in the ageing phase. No significant match was found between standard theoretical functions (Normal, Weibull, Gamma and Exponential) and the observed distribution of tree diameter. Also, it was confirmed that dominance of beech in all development phases is a consequence of its high competitive ability and its capacity to endure difficult environmental and biologically stressful conditions. The results revealed a series of structural models specific to these forest ecosystems, which can help managing forests under the selection system.

# Introduction

The Food and Agriculture Organization ([42]) defines virgin forests as “forests developed under natural conditions and undisturbed by human activities”. In Europe, in managed forests, sustainable management and silvicultural practices at the beginning of the 21st century aim to reproduce virgin forest complex structures and their multi-functional outcomes ([44]). The complexity of uneven-aged forest structures ensures economic value through superior timber quality, as well as high biodiversity at all levels ([34], [14]). The dynamics of virgin forests, their response to stress factors, natural competition, species exclusion, species complementarity and niche development have been the subject of research and a reference model, as well as an inspiration for theories and principles of silviculture ([34], [43]). Furthermore, due to the impact of climate change on managed forests and the new upcoming challenges, a different approach in the way forests are managed is needed. New practices and management strategies are required, and virgin forests can provide a model for future silvicultural practices ([22]) and forest management planning to address the challenge of climate change. New approaches should include silvicultural practices aimed at improving the forest stability otherwise reached by unmanaged forests ([31]). However, silvicultural treatments are indispensable, under certain circumstances, to achieve a higher diversity and stability then that attainable following the forest natural dynamics ([6]).

Virgin forests are a natural heritage of extraordinary value and timeless interest. Most of the world’s remaining virgin forests survive thanks to limits on logging and restrictions to human access. However, these limitations are now being removed because of new harvesting techniques and the high value of timber, thus virgin ecosystems are becoming increasingly rare ([34], [14]). Considering the negative impacts that have been recorded on forest succession ([21]), the growing research interest and concern on the survival of virgin forests are justified.

In this study, the term “virgin forest” is used according to the reference literature ([20], [23]), as synonymous of primaveral forest, primary forest and pristine forest ([39]). This study focuses on understanding and modelling virgin forest development phases, defined as theoretical concepts which help characterise the structural dynamics of natural forests ([38]) under the influence of internal and external factors ([47], [11]).

The development phases offer information on virgin forests’ growth dynamics, about how natural processes run without being disturbed by human activities and also how we can reproduce these structures through future management practices ([47], [44]). The aim of this study is to determine the characteristics of each development phase in order to develop solutions for structural models applied in forest management. The main objectives of the study are: (i) to characterise the structure of a virgin stand and its development phases; (ii) to describe the stand structure in relation with diameter (for the entire stand and for each development phase) by theoretical functions; (iii) to determine the relationships among tree species and their respective development phase, age and growth.

# Materials and methods

The study area is situated in the Curvature Carpathians - more specifically, in the Penteleu Massif, which is a part of the Buzau Mountains - in a temperate continental zone at an altitude of 1128 m a.s.l., on a southern slope (Fig. 1). The lithological substrate is composed of sandstone and dune sand with variable granulation. The dominant winds come from NNE and the average annual temperature is 6.1 °C. The inventoried plot is part of the Milea Viforâta Nature Reserve. It is circular in shape, with a radius of 56 m, and has an area of 1 ha. A tape was used to measure trees with a circumference > 25 cm at breast height, corresponding to a diameter of at least 8 cm. The height of each tree in the plot was measured using an ultrasound instrument system (Vertex IV®, Haglof, Sweden). The measured variables were: tree position, circumference, height, species, and storey position. Also, wood core samples were obtained with a Pressler’s increment borer, as radial growth, together with the structure of stand in relation with tree diameter, represent the most important element in describing the stand evolution and its development phases ([15]).

Fig. 1 - The study area (black cross) in the Milea Viforâta Nature Reserve (Romania).

Tree radial growth was measured on cores using the software CooRecorder® v. 7.4 (Cybis Elektronik & Data AB, Saltsjöbaden, Sweden) based on image analysis, and the validation of measurements was done by COFECHA software ([18]). A mean radial growth of the last ten years was computed for each tree (ir10). Nonlinear regression analysis using the “nls” function in R ([4]) was conducted for each species to determine a radial growth (ir) model using the following mathematical expression (eqn. 1):

$${ \bar {ir}}_{10}=a_0 + a_1 dbh +a_2 dbh^2$$

where ir10 is the mean ir of the last ten years for each tree and a0, a1 and a2 are the species-specific parameters to be estimated. The radial increment derived from wood core samples allows us to distinguish different growth rates among groups of trees, and thus to determine the stand’s development phases.

Theoretical distributions (normal - [35]; gamma - [16]; weibull - [28]; exponential - [25]) were fitted using the “fitdristr” function of the R “fitdistrplus” package ([9]). Kolmogorov-Smirnov (KS - [41]), Anderson-Darling (AD -[2]) and chi-square (CHI) tests were used to assess the fit of theoretical distribution to the experimental diameter distribution. Also, to compute the diameter limit of the stand, an exponential distribution of the following form was used ([27] - eqn. 2):

$$y = k \cdot e^{- \alpha x}$$

where y is the number of trees in a narrow diameter interval yx, x is the diameter at breast height, e is the base of the natural logarithm, and k and α are constants characterising a certain structure.

Volume growth was established using a method based on single-inventory information and on wood cores taken from standing trees ([13]). The following equation was used to calculate the volume of each individual tree ([13] - eqn. 3):

\eqalign{\log v = a_0 &+ a_1 \log dbh \\&+ a_2 \log^2 dbh \\ &+a_3 \log h \\&+ a_4 \log^2 h}

where h is the tree height (m), dbh is the diameter at breast height (cm) and a0, …, a4 are the regression coefficients (spruce: a0 = -4.18161, a1 = 2.08131, a2 = -0.11819, a3 = 0.70119, a4 = 0.148181; beech: a0 = -4.11122, a1 = 1.30216, a2 = 0.23636, a3 = 1.26562, a4 = -0.079661; fir: a0= -4.46414, a1 = 2.19479, a2 = -0.12498, a3 = 1.04645, a4 = -0.016848) established for the most important species in Romania’s forests ([13]).

To analyse structural diversity, the Gini index ([12]) was calculated for each development phase ([19]). Each tree was assigned to one of four virgin forest’s development phases ([11]): (i) optimum; (ii) ageing; (iii) breakdown and dieback; (iv) rejuvenantion. The classical method ([26]) was applied for delimitating the development phases on the ground, using the following stand characteristics: number of trees, basal area, volume growth, stand class, vitality, mortality, age structure, stability of the stand and regeneration (see also Tab. S1 in Supplementary material). GPS technology was used for spatial delimitation of the development phases.

All statistical analyses were conducted using the R software ([37]).

# Results

All development phases (optimum, ageing, breakdown and dieback, rejuvenation) were found in the study plot, but with a high variability in the size of the area covered by each phase (Tab. 1). The small area covered by early development phases is the consequence of the lower limit (8 cm) adopted for trees to be inventoried. The number of trees in each development phase was highly positively correlated with the basal area (r= 0.896; p<0.01562). The smallest number of trees was recorded for the optimum phase, followed by the rejuvenation phase, with 90 trees. The breakdown and dieback and ageing phases were well represented within the plot, with 178 and 112 trees, respectively. The optimum phase covered 14 per cent of the total surface area of the plot, with basal areas of 1.05 m2. With the exception of the optimum phase, mean radial growth showed a weak correlation with the average age of trees (r = 0.25, p<0.001), as well as with tree volume (r = 0.35, p<0.001). The volume per tree (overall aboveground tree volume for broadleaves and trunk volume for coniferous) ranged between 0.11 m3 and 3.28 m3.

Tab. 1 - Main structural characteristics of the four development phases recorded at the studied virgin forest.

Characteristics Development phase Total
Optimum Ageing Breakdown
and dieback
Rejuvenation
Number of trees 17 112 178 90 397
Covered area (%) 31 21 41 7 100
Basal area (m2) 1.08 7.01 20.06 15.58 43.73
Gini Index 0.57 0.78 0.69 0.52 -
Min radial growth (ir, mm) 0.42 0.15 0.14 0.15 0.14
Max radial growth (ir, mm) 2.04 3.56 4.36 4.33 4.36
Mean radial growth (ir, mm) 0.97 1.06 1.08 1.42 1.15
Min age (yrs) 32 11 24 27 11
Max age (yrs) 329 365 328 334 365
Average age (yrs) 115.88 89.17 112.46 145.72 113.57
Mean unitary volume (m3) 0.70 0.86 1.74 2.83 1.69

## Characterisation of the stand using dbh theoretical distributions

Theoretical functions (normal, gamma, weibull, exponential) were fitted to the experimental distribution of diameters (Fig. 2). All the statistical tests applied (Kolmogorov-Smirnov, Anderson-Darling and chi-square - Tab. 2) confirmed that none of these distributions corresponds to that observed in the virgin forest stand at Penteleu-Viforâta. Therefore, a mixture of three normal distributions was used to describe the stand structure (Fig. 3). The significant (p<0.001) fitting confirmed the existence of a three-modal distribution of tree diameters, where each peak corresponds to a storey position in the canopy (see also Tab. S2 in Supplementary material).

Fig. 2 - Theoretical (colored lines) vs. experimental (bars) distributions of tree dbh at the Penteleu-Viforâta virgin forest.

Tab. 2 - Goodness-of-fit of theoretical functions to the experimental distribution of diameters at the Penteleu-Viforâta virgin forest.

Theoretical
Function
Statistical tests
CHI Kolmogorov-Smirnov
NO/YES
Kolmogorov-
Smironov
Anderson-Darling
Normal 595.53 rejected 0.20 31.07
Weibull 212.92 rejected 0.15 14.09
Gamma 193.78 rejected 0.12 13.78
Exponential 217.63 rejected 0.24 19.78

Fig. 3 - The mixture of three normal distributions used to describe the distribution of diameters at the Penteleu-Viforâta virgin forest. The significant (p<0.001) fitting confirmed the existence of a three-modal distribution, where each peak corresponds to a storey position in the canopy.

Similar results were obtained when theoretical functions were fitted on observed diameter distribution of each development phase (whenever sufficient samples were available), except for the rejuvenation phase (Tab. 3) where significant (p<0.001) fitting was found for normal, gamma and weibull functions. In the case of all other phases and functions, fitting was rejected.

Tab. 3 - Goodness-of-fit of theoretical functions to the experimental distribution of diameters for each development phase at the Penteleu-Viforâta virgin forest. Results for the optimum phase are not reported due to the limited sample size.

Theoretical function Statistical
test
Phase
Ageing Rejuvenation Breakdown
and dieback
Normal KS 0.24 0.12 0.21
CHI 215.94 27.79 198.81
TEST Reject Fitted Reject
Weibull KS 0.21 0.07 0.15
CHI 94.99 15.57 74.25
TEST Reject Fitted Reject
Gamma KS 0.16 0.06 0.16
CHI 75.68 14.21 65.28
TEST Reject Fitted Reject
Exponential KS 0.31 0.19 0.23
CHI 115.37 24.52 80.02
TEST Reject Reject Reject

## Relation between development phase and species

Beech was the dominant species in all layers and phases, being present in large numbers in every development phase. Ageing, breakdown dieback and rejuvenation were the three phases which mainly cover the study plot, and they were well covered by all three species (beech, fir and spruce). As mentioned above, the limited number of trees falling in the optimum development phase is the consequence of the adoption of a lower limit (8 cm) for trees to be inventoried (Tab. 4).

Tab. 4 - Species distribution in relation with development phase.

Development phase Number of trees
Fir Beech Spruce Total
Ageing 21 82 9 112
Breakdown and dieback 52 108 18 178
Optimum 5 11 1 17
Rejuvenation 18 58 14 90
Total 96 259 42 397

## Relation between age and species

The predominance of trees belonging to lower age class (80-100 years old) was observed in the studied plot. However, a significant number of trees belonging to upper age class were also found, but with a lower frequency. Beech reach an age of 380 years in the Penteleu-Viforâta virgin forest, while fir and spruce are no longer present after 340 years (Tab. S3 in Supplementary material). Most of beech trees, however, are in the lower age classes, with the highest abundance recorded in the age class 80-120 years. Both beech and fir have an uneven-aged distribution characterised by a large number of trees in the younger classes, whereas for spruce, the distribution is platykurtic, with a similar number of trees in every age class.

## Relation between diameter and growth by species

The analysis of mean radial growth based on wood core samples allowed to build species-specific growth models (Fig. 4), which revealed high increment values for all the species. For beech, the highest radial increment (1.82 mm year-1) was recorded for trees with a diameter of 50 cm, while for fir and spruce higher increment values were observed at larger diameters (2.22 mm year-1 at 70 cm diameter and 2.43 mm year-1 at 75 cm diameter, respectively - Fig. 4). However, after reaching these maximum increment levels, radial growth rate decreases and trees enter a senescent overstory stage. This is characterised by high mortality with the stand density curve sloping downward, reaching the diameter limit. In our study, the diameter limit corresponds to the category of 84 cm (Fig. 5).

Fig. 4 - Species-specific growth models based on wood core samples analysis.

Fig. 5 - Stand diameter limit and maximum point of growth for beech, spruce and fir at Penteleu-Viforâta virgin forest.

# Discussion

Compared with even-aged stands, virgin forests are thought to better withstand natural disasters, such as insect outbreaks, windthrows and severe drought, thanks to the co-occurrence of different development phases which ensure long-term sustainability and continuity of the forest ecosystem ([29]). The typical process of identification of development phases in virgin forests can suffer from subjectivity ([11]), whereas the adoption of quantitative methods to investigate stand dynamics based on repetitive measurements combined with GIS data and stand structural characteristics is required to obtain more precise and reliable results.

In the studied forest, the number of trees per ha was 225-400 and the basal area was 32-47 m2. Such values are at the upper limit of the ranges reported in the literature for other old-growth European forests ([17], [47]). Also, the dominance of the ageing and breakdown and dieback phases at Penteleu-Viforâta is likely the result of the protection from wind provided by the surrounding reliefs. Indeed, the analysed plot is situated on a south-facing slope with dominant winds blowing from NNE, and this is similar to other virgin forests reported in the literature ([3]).

The values of the Gini coefficient observed at Penteleu-Viforâta reflect the structural characteristics of each development phase and support the hypothesis of a high stability of the stand ([7]). The lower values of Gini coefficient recorded in the rejuvenation phase indicate a reduced structural diversity. Contrastingly, high values were found in the ageing phase, thus indicating a higher structural diversity and stability. Moreover, the values of Gini coefficient recorded for each development phase were fairly similar to those reported in the literature for virgin forests ([47]), i.e., 0.69-0.76 for the optimum phase, 0.51-0.64 for the ageing phase and 0.51-0.60 for the rejuvenation phase. We believe that such values can be used as reference for further studies focused on virgin beech-fir-spruce forests in Europe. This information, along with the species growth models, could be relevant for the development of a model to manage cultivated forests based on evidences obtained from virgin forests ([32]), thereby increasing long-term functionality and stability of stands ([14]).

The diameter distribution of trees in our study area did not fit any of the most common distributions suited to managed forests, due to the structural complexity of each development phase in virgin forest ecosystems ([14]). It is known that diameters in uneven-aged forests do not follow a general theoretical distribution, although in some cases experimental data showed a fairly good fit to Weibull, Beta and Meyer distributions ([1]). In fact, the structure of a virgin forest stand keeps changing as a result of the balance of inputs and outputs. Inputs are provided by new areas covered by regeneration, and the outputs are the result of the high tree mortality rate due to the strong competition that exists in a virgin forest.

Regarding the relation between development phases and species, the dominance of beech in all development phases, especially in lower-diameter categories, confirms previous findings on this species ([10], [44]). In virgin beech-fir-spruce forests, disturbances are frequent, though balanced by rapid stand growth ([33]). Since virgin ecosystems are not directly influenced by any anthropic disturbances, this growth pattern may occur as a response to climate change ([44]). Furthermore, our results confirm the capacity of beech for regeneration and survival in the understory, showing an increase in growth even under the canopy ([46], [30]).

Leuschner & Meier ([24]) in a study on 42 European forest tree species, reported a maximum longevity of 450, 450 and 300 years for beech, fir and spruce, respectively. In this study, we found several individuals of beech reaching up to 380 years of age, confirming that the Penteleu-Viforâta forest is an optimum site for beech. Although the general climatic conditions suggest a good vegetation site also for the other two species (fir and spruce), their capacity to grow in the understory is lower than that of beech, and this affects their growth dynamics and longevity. Indeed, beech shows a superior capacity of surviving under stress conditions, making it more competitive than fir and spruce. This is reflected by the high-intensity growth of beech in the optimum phase, compared with the two conifers.

Uneven-aged stands are more productive than even-aged stands ([8]); this also holds for mixed forests compared to pure stands ([36]). The most common natural mixed composition in uneven-aged forests of Central Europe consists of beech, spruce and fir, due to their ecological complementarity. Nonetheless, these species show different growth dynamics and competitive ability, particularly in an uneven-aged forest. As mentioned above, beech has a higher capacity to survive in the understory ([5]) due to a well-developed root system reaching deeper soil layers. Moreover, beech growth has a low variability compared to the conifer species ([45]), as a consequence of its capacity to endure a wider range of environmental conditions (including drought). In general, beech shows high growth rates up to approximately 100-120 years of age, while spruce and fir outstrip beech productivity and maintain growth longer than beech and at a higher level, as indicated from the different peaks of mean radial growth for each species observed in this study.

The identification of the maximum growth increment by wood core analysis and the ascertainment of the diameter limit for each species could help define silvicultural practices aimed at optimizing selective cutting in uneven-aged and/or mixed stands. For virgin forests, the diameter limit is taken, from the physiological point of view, as the last category of the diameter distribution a living tree can reach ([40]). In uneven-aged forest management, the maximum point of species growth can be used for predicting the diameter limit of the stand. The diameter limit is dynamic, depending on the site species and the stand history. According to Meyer’s negative exponential distribution, the diameter limit is the diameter category with a frequency of one tree.

# Conclusions

Virgin forests are a natural heritage of great value, offering solutions and structural models for sustainable forest management. The analysis of development phases of virgin forests may provide important information to be included in silvicultural practices aimed at enhance the stability and functionality of these ecosystems.

The structural diversity of the Penteleu-Viforâta virgin forest was analyzed by the Gini index. We observed the dominance of the breakdown and dieback as well as ageing development phases, while maximum and minimum stability were recorded for ageing and rejuvenation phases.

Tree diameter distribution did not match any of the tested theoretical functions, due to the complexity of the stand structure and its dynamics continuously changing in time and space. The study of the relation between development phases and species confirmed the dominance of beech in all phases as a consequence of its competitive ability and capacity to endure a wider range of stressful conditions. Moreover, the wood core analysis revealed a different mean radial increment and maximum radial growth for the three species analyzed, which should be taken into account by forest management planners when establishing the diameter limit.

Virgin forests can be considered as in situ laboratories for forest research, particularly for better understanding the principles underlying the stability and functionality of these ecosystems. These types of forests provide a vast array of resources, offering at the same time a higher level of protection, conservation and biodiversity. In this context, further research on development phases in virgin forests are needed in order to devise sustainable forest management strategies.

# References

(1)
Akhavan R, Sagheb-Talebi K, Zenner EK, Safavimanesh F (2012). Spatial patterns in different forest development stages of an intact old-growth Oriental beech forest in the Caspian region of Iran. European Journal of Forest Research 131: 1355-1366.
CrossRef | Gscholar
(2)
Anderson TW, Darling DA (1954). A test of goodness of fit. Journal of the American Statistical Association (268): 765-769.
CrossRef | Gscholar
(3)
Anić I, Mikac S (2008). Struktura, tekstura i pomladivanje dinarske bukovo jelove paršume Corkova Uvala [Structure, texture and regeneration of Dinaric beech-fir virgin forest of Corkova Uvala]. Šumarski list 11 (12): 505-515. [in Croatian]
Gscholar
(4)
Baty F, Ritz C, Charles S, Brutsche M, Flandrois JP, Delignette-Muller ML (2015). A toolbox for nonlinear regression in R: the package nlstools. Journal of Statistical Software 66 (5): 1-21.
CrossRef | Gscholar
(5)
Bose A, Weiskittel A, Wagner R (2017). Occurrence, pattern of change, and factors associated with American beech-dominance in stands of the northeastern USA forest. Forest Ecology and Management 392: 202-212.
CrossRef | Gscholar
(6)
Brzeziecki B, Bielak K, Bolibok L, Drozdowski S, Zajaczkowski J, Zybura H (2018). Structural and compositional dynamics of strictly protected woodland communities with silvicultural implications, using Bialowieza Forest as an example. Annals of Forest Science. 75: 89.
CrossRef | Gscholar
(7)
Chivulescu S, Badea O, Tomescu R, Silaghi D, Leca S, Turcu DO (2014). Structural features of virgin beech forests in Semenic Mountains. The dynamic structure of virgin beech forest P20 Semenic between 2005-2013. Proceedings of Romanian Academy-Series B: Chemistry, Life Science and Geoscience 16 (2): 115-124.
Gscholar
(8)
Danescu A, Albrecht AT, Bauhus J (2016). Structural diversity promotes productivity of mixed, uneven-aged forests in southwestern Germany. Oecologia 182: 319-333.
CrossRef | Gscholar
(9)
Delignette-muller ML, Dutang C (2015). fitdistrplus: an R package for fitting distributions. Journal of Statistical Software 64 (4): 1-34.
CrossRef | Gscholar
(10)
Diaci J, Rozenbergar D, Mikac S, Anic I, Hartman T, Boncina A (2008). Long-term changes in tree species composition in old-growth Dinaric beech-fir forest. Glasnik za Sumske Pokuse 42: 13-28.
Gscholar
(11)
Feldmann E, Glatthorn J, Hauck M, Leuschner C (2018). A novel empirical approach for determining the extension of forest development stages in temperate old-growth forests. European Journal of Forest Research 137: 321-335.
CrossRef | Gscholar
(12)
Gini C (1912). Variabilità e mutabilità [Variability and mutability]. Cuppini, Bologna, Reprinted in “Memorie di metodologica statistica” (Pizetti E, Salvemini T eds). Libreria Eredi Virgilio Veschi, Rome, Italy, pp. 156. [in Italian]
Gscholar
(13)
Giurgiu V, Decei I, Draghiciu D (2004). Metode si tabele dendrometrice [Forest mensuration methods and tables]. Editura Ceres, Bucharest, Romania, pp. 575. [in Romanian]
Gscholar
(14)
Giurgiu V (2013). Padurile virgine si cvasivirgine ale Romniei [Virgin and quasi-virgin Romanian forests]. Ed. Academiei Romne, Bucharest, Romania, pp. 390. [in Romanian]
Gscholar
(15)
Hassani M, Amani M (2010). Investigation on structure of oriental beech (Fagus orientalis Lipsky) stand at optimal stage in Sangdeh forest. Iranian Journal of Forest and Poplar Research 18 (2): 163-176.
Online | Gscholar
(16)
Hogg RV, Craig AT (1978). Introduction to mathematical statistics (4th edn). Macmillan, New York, USA, pp. 438.
Gscholar
(17)
Holeksa J, Saniga M, Szwagrzyk J, Czerniak M, Staszynska K, Kapusta P (2009). A giant tree stand in the West Carpathians: an exception or a relic of formerly widespread mountain European forests? Forest Ecology and Management 257 (7): 1577-1585.
CrossRef | Gscholar
(18)
Holmes RL (1983). Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin 43: 69-78.
Online | Gscholar
(19)
Klopcic M, Boncina A (2011). Stand dynamics of silver fir (Abies alba Mill.) - European beech (Fagus sylvatica L.) forests during the past century: a decline of silver fir? Forestry 84 (3): 259-271.
CrossRef | Gscholar
(20)
Korpel S (1978). Obnova lesnych porastov v rubáòovom sposobe hospodárenia [Restoration of forest stands in the ruby way of management]. Pesteni lesù, SZN, Prague, Czech Republic, pp. 216-359 [in Slovakian]
Gscholar
(21)
Kral F (1979). Spät- und postglaziale Waldgeschichte der Alpen auf Grund der bisherigen Pollenanalysen [Late and postglacial forest history of the Alps based on the previous pollen analyzes]. Veröffentlichungen des Institutes für Waldbau an der Universität für Bodenkultur, Agrarverlag, Vienna, Austria, pp. 175. [in German]
Gscholar
(22)
Král K, Danek P, Janík D, Kruček M, Vrška T (2018). How cyclical and predictable are Central European temperate forest dynamics in terms of development phases? Journal of Vegetation Science 29 (1): 84-97.
CrossRef | Gscholar
(23)
Leibundgut H (1982). Europäische Urwälder der Bergstufe: dargestellt für Forstleute, Naturwissenschafter und Freunde des Waldes [European primeval forests of the Bergstufe: represented for foresters, naturalists and friends of the forest]. Haupt, Bern-Stuttgart, pp. 308. [in German]
Gscholar
(24)
Leuschner C, Meier IC (2018). The ecology of Central European tree species: trait spectra, functional trade-offs, and ecological classification of adult trees. Perspectives in Plant Ecology, Evolution and Systematics 33: 89-103.
CrossRef | Gscholar
(25)
Marshall AW, Olkin I (1967). A multivariate exponential distribution. Journal of the American Statistical Association 62 (317): 30-44.
CrossRef | Gscholar
(26)
Mayer H (1976). Gebirgswaldbau - Schutzwaldpflege [Mountain forestry, protective forest care]. Gustav Fisher Verlang, Stuttgart, Germany, pp. 435. [in German]
Gscholar
(27)
Meyer HA (1952). Structure, growth, and drain in balanced uneven-aged forests. Journal of Forestry 50 (2): 85-92.
CrossRef | Gscholar
(28)
Nawaz S, Nazrul IM (1980). The Weibull distribution as a general model for forecasting technological change. Technological Forecasting and Social Change 18 (3): 247-56.
CrossRef | Gscholar
(29)
Orman O, Dobrowolska D (2017). Gap dynamics in the Western Carpathian mixed beech old-growth forests affected by spruce bark beetle outbreak. European Journal of Forest Research 136: 571-581.
CrossRef | Gscholar
(30)
Orman O, Dobrowolska D, Szwagrzyk J (2018). Gap regeneration patterns in Carpathian old-growth mixed beech forests - Interactive effects of spruce bark beetle canopy disturbance and deer herbivory. Forest Ecology and Management 430: 451-459.
CrossRef | Gscholar
(31)
O’Hara KL (2016). What is close-to-nature silviculture in a changing world? Forestry 89 (1): 1-6.
CrossRef | Gscholar
(32)
Pach M, Podlaski R (2015). Tree diameter structural diversity in Central European forests with Abies alba and Fagus sylvatica: managed versus unmanaged forest stands. Ecological Research 30 (2): 367-384.
CrossRef | Gscholar
(33)
Paluch JG, Kolodziej Z, Pach M, Jastrzebski R (2015). Spatial variability of close-to-primeval Fagus-Abies-Picea forests in the Western Carpathians (Central Europe): a step towards a generalised pattern. European Journal of Forest Research 134: 235-246.
CrossRef | Gscholar
(34)
Parviainen J (2005). Virgin and natural forests in the temperate zone of Europe. Forest Snow and Landscape Research 18: 9-18.
Gscholar
(35)
Patel JK, Read CB (1996). Handbook of the normal distribution (2nd edn). Statistics Textbooks and Monographs, CRC Press, New York, USA, pp. 1-427.
Online | Gscholar
(36)
Pretzsch H, Bielak K, Block J, Bruchwald A (2013). Productivity of mixed vs. pure stands of oak (Quercus petraea (Matt.) Liebl. and Quercus robur L.) and European beech (Fagus sylvatica L.) along an ecological gradient. European Journal of Forest Research 132 (2): 263-280.
CrossRef | Gscholar
(37)
R Core Team (2014). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
Online | Gscholar
(38)
Remmert H (1991). The mosaic-cycle concept of ecosystems. In: “The Mosaic-Cycle Concept of Ecosystems” (Billings WD, Golley F, Lange OL, Olson JS, Remmert H eds). Ecological Studies, vol. 85, Springer, Berlin, Heidelberg, Germany, pp. 168.
Gscholar
(39)
Schuck A, Parviainen J, Bücking W (1994). A review of approaches to forestry research on structure, succession and biodiversity of undisturbed and semi-natural forests and woodlands in Europe. European Forest Institute, Joensuu, Finland, pp. 42.
Gscholar
(40)
Seceleanu I (2012). Amenajarea Padurilor - Organizare si conducere structurala [Forest management - Organization and management]. Editura Ceres, Bucharest, Romania, pp. 505. [in Romanian]
Gscholar
(41)
Stephens MA (1979). Tests of fit for the logistic distribution based on the empirical distribution function. Biometrika 66 (3): 591-595.
CrossRef | Gscholar
(42)
UN-ECE/FAO (2000). Forest resources of Europe, CIS, North America, Australia, Japan and New Zealand (industrialized temperate/boreal countries): UN-ECE/FAO contribution to the global forest resources Assessment 2000. UN Economic Commission for Europe, FAO, Rome, Italy, pp. 445.
Gscholar
(43)
Vacek Z, Podrázsky V, Bílek L, Bulušek D, Stefančík I, Remeš J, Stícha V, Ambro R (2014). Structural diversity of autochthonous beech forests in Broumovske Steny National Nature Reserve, Czech Republic. Austrian Journal of Forest Science 131: 191-214.
Gscholar
(44)
Visnjić C, Balić B, Vojniković S, Mekić F, Uzunović A (2015). Structural characteristics, dynamics and texture development of virgin forest Ravna Vala on Bjelašnica. Faculty of Forestry, University of Sarajevo, vol. 2015, issue 2, Bosnia-Hercegovina, pp. 100-112.
Gscholar
(45)
Vrška T, Adam D, Hort L, Kolár T, Janík D (2009). European beech (Fagus sylvatica L.) and silver fir (Abies alba Mill.) rotation in the Carpathians - A developmental cycle or a linear trend induced by man? Forest Ecology and Management 258: 347-356.
CrossRef | Gscholar
(46)
Wagner S, Collet C, Madsen P, Nakashizuka T, Nyland RD, Sagheb-Talebi K (2010). Beech regeneration research: from ecological to silvicultural aspects. Forest Ecology and Management 259: 2172-2182.
CrossRef | Gscholar
(47)
Zenner EK, Sagheb-talebi K, Akhavan R, Peck JE (2015). Integration of small-scale canopy dynamics smoothes live-tree structural complexity across development stages in old-growth Oriental beech (Fagus orientalis Lipsky) forests at the multi-gap scale. Forest Ecology and Management 335: 26-36.
CrossRef | Gscholar

# Supplementary Material

#### Authors’ Affiliation

(1)
Serban Chivulescu
Albert Ciceu 0000-0002-4090-8705
Stefan Leca 0000-0002-6110-8379
Bogdan Apostol 0000-0002-5935-4407
National Institute for Research and Development in Forestry ”Marin Dracea” (Romania)
(2)
Albert Ciceu 0000-0002-4090-8705
Transilvania University of Brasov (Romania)
(3)
Octavian Popescu 0000-0002-2597-7155
Institute of Biology Bucharest, Romanian Academy (Romania)
(4)
Octavian Popescu 0000-0002-2597-7155
Molecular Biology Center, Institute of Interdisciplinary Research in Bio-Nano-Sciences, Babes-Bolyai University, Cluj-Napoca (Romania)

#### Corresponding author

Albert Ciceu
albert.ciceu@yahoo.ro

#### Citation

Chivulescu S, Ciceu A, Leca S, Apostol B, Popescu O, Badea O (2020). Development phases and structural characteristics of the Penteleu-Viforâta virgin forest in the Curvature Carpathians. iForest 13: 389-395. - doi: 10.3832/ifor3094-013

#### Paper history

Accepted: Jun 17, 2020

First online: Sep 01, 2020
Publication Date: Oct 31, 2020
Publication Time: 2.53 months

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

#### Breakdown by View Type

(Waiting for server response...)

#### Article Usage

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

Breakdown by View Type
HTML Page Views: 5787
Abstract Page Views: 452

Web Metrics
Days since publication: 391
Overall contacts: 6978
Avg. contacts per week: 124.93

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

(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