Next Article in Journal
Modelling the Effect of Weed Competition on Long-Term Volume Yield of Eucalyptus globulus Labill. Plantations across an Environmental Gradient
Next Article in Special Issue
Landscape Planning—Paving the Way for Effective Conservation of Forest Biodiversity and a Diverse Forestry?
Previous Article in Journal
Effects of Leaf Age and Exogenous Hormones on Callus Initiation, Rooting Formation, Bud Germination, and Plantlet Formation in Chinese Fir Leaf Cuttings
Previous Article in Special Issue
Spatial and Temporal Variation of NDVI in Response to Climate Change and the Implication for Carbon Dynamics in Nepal
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 9
Issue 8
10.3390/f9080479
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Five Decades of Structural and Compositional Changes in Managed and Unmanaged Montane Stands: A Case Study from South-East Europe

by
Srđan Keren
1,*,
Milan Medarević
2,
Snežana Obradović
2 and
Brane Zlokapa
3
1
Department of Biometry and Forest Productivity, Faculty of Forestry, University of Agriculture in Krakow, al. 29-Listopada 46, 31-425 Krakow, Poland
2
Faculty of Forestry, University of Belgrade, 11000 Belgrade, Serbia
3
Faculty of Forestry, University of Banja Luka, 78000 Banja Luka, Bosnia and Herzegovina
*
Author to whom correspondence should be addressed.
Forests 2018, 9(8), 479; https://doi.org/10.3390/f9080479
Submission received: 3 July 2018 / Revised: 29 July 2018 / Accepted: 4 August 2018 / Published: 7 August 2018
(This article belongs to the Special Issue Forest Landscape Ecology: Linking Past, Present, and Future Data)
Browse Figures
Versions Notes

Abstract

:
The recent research has indicated that restoration of old-growth attributes such as large-sized living trees and snags contributes to sustaining biodiversity on the landscape level. The extent to which these attributes are restored, maintained, or diminished by total salvage logging, selection silviculture, and strict protection has been partly investigated in the past. However, studies examining the influence of partial salvage logging are largely absent. Thus, we compared long-term structural and compositional changes in three montane beech-fir-spruce stands in Serbia that were exposed to different management regimes for five decades (partial salvage logging, selection silviculture, and strict protection). Tree species composition of partly salvaged stand and selection stand significantly differed from that in the adjacent unmanaged stand. However, the diameter distributions of compared stands often exhibited the same structural forms in certain periods, despite the greater share of large-size trees in the unmanaged stand. The study indicated that managing for old-growth attributes such as large trees may be possible by applying not only rotated sigmoid and negative exponential structures, but also the increasing-q diameter structure as high basal areas in studied beech-fir-spruce stands did not impair the ingrowth of young trees when conifers dominated the upperstory. The study further revealed that partial salvaging may serve as a sound alternative to promoting old-growth attributes such as large veteran trees and snags.

1. Introduction

Strictly protected or otherwise created areas of old-growth forests make up a very small percentage of the overall forest area worldwide [1,2,3]. Therefore, strict protection alone is considered to be insufficient for sustaining old-growth attributes and species diversity over larger landscapes [2,4,5]. Parviainen et al. [2] indicated that the more management activities follow natural processes in production forests, the less there is a need for the total conservation of forests. In this context, single-tree selection silviculture has been classified as a close-to-nature system, but this approach aims at choosing mostly those natural processes that best achieve management goals such as target species composition, optimal growth, or high-quality timber [6]. Therefore, some authors took a step further by proposing a more comprehensive natural forest management that, in addition to sustainable diameter structure, natural regeneration and tree species composition, also includes the presence of large living veteran trees and large snags [4,7,8]. However, the problem of the lack of old-growth attributes was not acknowledged until relatively recently [9]. As a consequence, most studies focusing on the restoration of old-growth characteristics very rarely cover longer time spans [10]. For this reason, knowledge on restoration of old-growth attributes is still lacking.
The importance of rotated sigmoid, negative exponential, and increasing-q structural forms [4,11,12,13] lies in the fact that they can support high basal areas and allow adequate recruitment of large trees without endangering the overall tree-size demographic equilibrium. Moreover, the creation of stand structures with the full-range of tree sizes, including live veteran trees and standing dead trees, is crucial for sustaining biodiversity [14]. However, such structures are deficient in many areas, as forests have been mostly used for commercial reasons [15]. Nevertheless, despite high demands for wood products, many countries are now making an effort to restore some of the extinct old-growth attributes such as snags and large-sized living trees as they represent indispensible habitats for a large number of species [4,7].
Information about the structural and compositional changes of forests that are subject to different harvesting regimes may serve as a guideline for restoring old-growth conditions in managed forests, that is, at the larger landscape scale [10]. For management purposes, it is thus vital to investigate and understand how forest ecosystems respond to different silvicultural activities. Several studies in North America have investigated the restoration of old-growth attributes in managed forests through the application of a structural complexity enhancement approach and selection systems [4,8,16,17] as well as extended rotation forestry [18,19]. In contrast, there have been only few similar studies in Europe [10,13,14]. Moreover, long-term studies of structural dynamics and compositional changes have been rare for many ecologically and economically important forest types around the world [19].
The diversity of forest structures and natural tree species composition in most parts of Europe was greatly reduced by clear-cutting after the late 1700s [14,20,21]. European beech (Fagus sylvatica L.; hereafter beech) was especially affected, whereas conifers such as Norway spruce (Picea abies L. Karst; hereafter spruce) and Scots pine (Pinus sylvestris L.) were favored. However, forestry activities in some countries have recently led to increased structural and compositional diversity [22,23]. This indicates that the paradigm of close-to-nature silviculture has become partly accepted in practice, at least with respect to certain structural attributes related to live trees.
Since forests in South-East Europe are often located on limestone and steep slopes, and thus susceptible to erosion, the history of forest management in this part of Europe was mostly characterized by gradual conversion of old-growth forests to managed uneven-aged forests, usually by the application of selection silviculture [24]. This approach may have contributed to better conservation of trees having diameters at breast height (dbh) between 50 cm and 80 cm [25]. However, the policies in South-East Europe related to conserving standing dead trees (snags) and very large living trees (>80 cm dbh) have focused mostly on strictly protected areas.
A mixture of beech, spruce, and silver fir (Abies alba Mill.; hereafter fir) is often the dominant forest type in the mountainous regions of Europe [26]. In Serbia, mixed stands composed of these species cover relatively small areas [27,28] because both fir and spruce in this country are close to the southern edge of their natural range, which makes them more susceptible to climate changes and disturbance events. For this reason, these stands were given special attention and were often thoroughly monitored. The basis of our research consists of long-term data from partly salvaged stand, selection stand, and an unmanaged stand on Tara mountain in western Serbia. Salvage logging may have both positive and negative effects on forest ecosystems [29]. Recently, however, especially negative effects are studied after all damaged trees are removed following large or medium natural disturbances [30,31,32,33]. In the Dinaric forests, such events occur also occasionally, but small-scale gaps are predominant, whereas the wind tends to be the major disturbance agent [34]. In this study, partial salvage logging was conducted in response to low severity disturbances rather than stand-replacing disturbance (details are provided in the Section 2.1).
While the studies on the effects of selection silviculture and strict protection have been moderately frequent, the studies on the influence of partial salvage logging on stand structural characteristics have been almost nonexistent. Consequently, the major goals of the present study were to (i) compare whether three differently managed beech-fir-spruce stands (partial salvaging, selection management, and strict protection) differed across five decades with respect to the form of the diameter distribution and tree species composition, (ii) investigate how different management regimes influenced the recruitment of large-sized living trees, and (iii) test for differences between partly salvaged stand and unmanaged stand in terms of the density and basal area of snags.

2. Materials and Methods

2.1. Study Area and Management Regimes

The study area is located in the Tara National Park (NP) in western Serbia. The Tara mountain is situated in the range of Dinaric Mountains between 43°51′ and 43°57′ N and 17°03′ and 17°11′ E. Different types of management are applied in this area, including the selection system, partial salvaging, and strict protection. At the altitude between 900 m to 1350 m above sea level, the dominant forest type is a mixture of beech, fir, and spruce on limestone covered mostly with brown soil types. The climate for the period 1975–2005 was assessed to be moderately continental with mean yearly precipitation of 991 mm and mean yearly temperature of 7.9 °C [28]. The study area has not experienced ungulate over-browsing, as has been the case in the mixed forests of some countries in South-East Europe [35,36].
The first registered cutting activities in the Tara mountain date back to 1852. Mostly conifers with dbh above 40 cm were cut [37], while beech was spared until the mid 20th century [38]. Both conifers and beech were used to satisfy the needs of local wood processing mills for timber, with some being exported. Nevertheless, the landscape did not experience dramatic changes in tree species composition, since mostly single-tree harvesting was performed [24,39]. The first organized management plan was made in 1960. Since then, inventories in selection stand and an adjacent unmanaged stand (locally called Racanska Sljivovica) have been conducted every 10 years. The unmanaged stand was declared as a strict reserve in 1957, but regular measurements began simultaneously with the inventories of selection stand in 1960.
A few years prior to national inventory, a smaller area (8.0 ha) that was only partly salvaged was first measured in 1955, although with somewhat less frequent measurements (see Data collection). This stand has been only partly salvaged in order to leave standing dead trees that were deemed not to jeopardize the stability and health of surrounding live trees. This means that single-tree mortality (e.g., due to competition or aging) has been partly tolerated; however, if a tree was broken by agents such as wind or snow, or if it was determined to suffer from beetle or fungi attack, it was felled and removed from the stand. This approach of low-intensity cutting was part of an experiment with the intention to increase the share of large trees and snags, but it differs from structural complexity enhancement approach described in the USA [4] as the latter had certain constraints. For instance, the rotated sigmoid diameter distribution was prescribed as one of the objectives, and target diameter was set to 90 cm. In Tara, the initial objective in 1955 was to establish a negative exponential form with high basal area, no target diameter, and a certain amount of snags were allowed if they did not disturb a demographic equilibrium of live trees. Such an approach is not traditional in the region but was used in the study area solely as an experimental treatment.
On the other hand, no harvesting was performed in the unmanaged stand Tara, at least from 1957 to the present [40]. All of the mentioned stands are located in close proximity to each other and have been exposed to virtually the same climate and natural disturbance events, which enabled us to differentiate and examine the influence of the different types of forest management.

2.2. Data Collection

The stands in NP Tara were given special attention in the frame of the national forest inventory in Serbia where relatively large areas were fully callipered over time [41]. Since unmanaged stand and partly salvaged stand had similar site conditions (soil type, inclination, altitude, etc.), for proper comparison we also chose one stand from selection forest that had similar site conditions to those found in adjacent unmanaged stand and partly salvaged stand. The harvests in selection stand were conducted once every ten years, whereas target basal area was not set as a single value but as the range from 31 m2·ha−1 to 38 m2·ha−1. The area of selection stand that was studied amounted to 11.2 ha, whereas the area of the unmanaged stand was 15.4 ha. Both selection stand and unmanaged stand were fully callipered every ten years from 1960 to 2010. On the other hand, the measurements in partly salvaged stand started five years earlier and they were somewhat less frequent (1955, 1975, and 2005). Nevertheless, we find them important for characterizing the long-term structural and compositional dynamics under partial salvage logging. The data for partly salvaged stand were collected on the total area of 8.0 ha where merely partial salvaging has been conducted.
At every inventory, the dbh of all stems >10.0 cm were measured to the nearest 1.0 mm. The stands were not sampled or replicated, but complete census of all trees was conducted in each of three studied stands. The total numbers of trees that were measured during every census are provided in the Appendix A, Table A1. For meaningful comparison of partly salvaged stand with selection stand and unmanaged stand in Tara, we used only the data from selection stand and unmanaged stand that were collected in the years 1960, 1980, and 2010, as those inventory years were closest to the measurement years in partly salvaged stand.
Standing dead trees with dbh >10.0 cm and height >1.30 m were callipered from 1975 to 2005 in partly salvaged stand, whereas snags in unmanaged stand were inventoried only in 2010. Since repetitive measurements on snags were not available, the mortality rate could not be determined. Nevertheless, the one-time measurement in unmanaged stand provided useful information with respect to the long-term outcomes of strict protection, whereas the two measurements (in 1975 and 2005) gave us an insight into the density of snags in partly salvaged stand. Concerning selection stand in Tara, snags were not inventoried, since managers regularly remove damaged and dead trees from selection stands. Thus, their occurrence in selection stands in this part of Europe may be only temporary and accidental rather than the result of planned activities.

2.3. Data Analysis

A long-term experiment spanning 50 years was used to examine the effects of partial salvage logging, selection management, and strict protection on live-tree size distributions, species composition, and the abundance of standing dead trees in the study area. Diameter distributions were constructed using 10 cm diameter classes as applied by Silver et al. [19]. For detailed quantification that provides clear insight into the changes in dbh distribution shapes at the species- and stand level, we applied the method proposed by Janowiak et al. [42]. Mixed (linear, polynomial, cubic) regression models were calculated and best-fitting models were selected based on lowest Akaike information criterion (AIC) [43], whereas log10 of tree density +1 was regressed on all combinations of size class midpoints to the third exponent. Depending on the signs of coefficients, diameter distribution shapes could be classified into: Negative exponential (−dbh); increasing-q (−dbh2 or −dbh3); rotated sigmoid (−dbh + dbh2 − dbh3); concave (−dbh + dbh2, −dbh + dbh3); or unimodal (+dbh − dbh2, +dbh − dbh3, +dbh2 − dbh3) (see [42,44]). If best model indicated variable shape (−dbh2 + dbh3 or +dbh − dbh2 + dbh3), then the second best model was selected. The changes in live-tree size distribution shapes in selection stand, unmanaged stand, and partly salvaged stand were observed over the period of five decades, whereas partly salvaged stand was compared to selection stand and unmanaged stand in different time periods.
In order to test for differences in tree species composition among these stands, we applied a two-variable χ2-test, whereas species frequencies rather than their proportions or percentages were used. To avoid biased χ2 estimation for low expected frequencies (see [45]), we grouped the initial 10 cm classes into four size classes matching the following size categories: poles (dbh of 10.1 cm–20.0 cm), small trees (dbh of 20.1 cm–30.0 cm), medium trees (dbh of 30.1 cm–50.0 cm), and large trees (dbh > 50.0 cm). A one-variable χ2-test was applied for comparison of snag density distributions in partly salvaged stand and unmanaged stand.
The intensity of basal area dynamics for large trees >50 cm dbh was analyzed by using the large tree basal area index (LTIba), which was adopted after Klopcic and Boncina [46], who employed a similar index for the total stand volume. In our study, the LTIba represents a measure of the relative change in basal area of large trees (>50 cm dbh), calculated as the quotient between this basal area at a particular point in time and the respective basal area of large trees at the first forest inventory. The index was computed for each tree species and also for total trees. The possible index values may range theoretically from −1 (all large trees removed from the stand before the final measurement) to +∞.

3. Results

3.1. Structural Characteristics of Live Trees

In the first observation period (1955/1960), partly salvaged stand had the highest density, with 513 stems·ha−1, selection stand had 497 stems·ha−1, and unmanaged stand had the lowest density with 355 stems·ha−1. A similar relationship was also observed two decades later (1975/1980). Nevertheless, by the end of the observation period (2005/2010), selection stand had the highest stem density, while the number of pole trees in partly salvaged stand decreased significantly, to the extent that it ended up with the lowest overall stem density (Figure 1).
Diameter distribution shapes were compared among partly salvaged stand, selection stand, and unmanaged stand in different time periods. At the beginning stage of observation (1955/1960), fir, spruce, and total trees in partly salvaged stand as well as total trees in selection stand had negative exponential distribution, whereas total trees in unmanaged stand exhibited an increasing-q distribution. At this time, selection stand exhibited higher demographic stability compared to unmanaged stand as the dominant species of fir and spruce had rotated sigmoid shapes in the first, and increasing-q and unimodal distributions in the latter, respectively. In the first period, only beech had the same distribution (increasing-q) in all three sites. In 1975/1980, fir in selection stand and unmanaged stand had the same shape (rotated sigmoid) as well as total trees in selection stand and partly salvaged stand (negative exponential), whereas the structural forms of all other comparative pairs differed by species and also total trees. Eventually, in the last observation period in 2005/2010, in selection stand all distributions by species and for total trees were rotated sigmoid, and in the unmanaged stand three out of four possible cases were also rotated sigmoid (only beech had the unimodal shape). On the other hand, partly salvaged stand exhibited more varied shapes, with fir having the rotated sigmoid shape, beech the unimodal shape, and spruce and total trees having the increasing-q distribution form (Table 1; Figure 1, Figure 2 and Figure 3).
From the perspective of temporal changes, the rotated sigmoid distribution showed stability at the species level for 50 years in selection stand (fir) and for 30 years in selection stand (spruce) and unmanaged stand (fir). The negative exponential and increasing-q shapes were less stable, as they were maintained for only two decades. The negative exponential shape, after initially being determined for fir and spruce in partly salvaged stand, eventually transformed into the rotated sigmoid (fir) and increasing-q (spruce) shapes. A similar demographic trajectory was determined for total trees in selection stand. Only beech in partly salvaged stand and unmanaged stand showed a different structural trajectory by switching from the initial increasing-q form into an unimodal distribution (Table 1; Figure 1, Figure 2 and Figure 3).

3.2. Tree Species Composition

At the beginning of observation (1955/1960), significant differences in tree species composition were determined among partly salvaged stand, selection stand, and unmanaged stand. Despite different cutting strategies, however, no significant difference was found five decades later between selection stand and partly salvaged stand, whereas these two stands differed significantly from adjacent unmanaged stand throughout whole study period (χ2-tests, p < 0.05). Partly salvaged stand exhibited the largest decrease in overall stem density over time (Appendix A, Table A2). Although we were not able to investigate spatial interactions among tree species based on the available dataset, it is important to point out that the share of beech was higher in partly salvaged stand compared to selection stand and unmanaged stand throughout the whole study period (Table 2; Figure 4). In the final observation year, the share of large beech trees (>50 cm dbh) amounted to 40% of the total basal area of large trees in partly salvaged stand, whereas in selection stand and unmanaged stand the share of such beech trees amounted to only 25% and 16%, respectively (Appendix A, Table A3 and Table A4). With respect to total basal area, it remained highest in unmanaged stand throughout the whole period of five decades (42.4 m2·ha−1–56.7 m2·ha−1), moderate in partly salvaged stand (38.1 m2·ha−1–49.9 m2·ha−1), and smallest in selection stand (34.0 m2·ha−1–37.0 m2·ha−1) (Figure 4).
Whereas in partly salvaged stand all species experienced a decrease in the density of pole trees by 43%, their recruitment in selection stand and unmanaged stand simultaneously increased by 2% and 102%, respectively. In selection stand, only fir experienced a relatively small decrease in pole trees (17%), while in unmanaged stand, beech poles almost disappeared from 1960 to 2010. On the other hand, the density of fir poles in this unmanaged stand increased in the same period by 74%. Regarding the density of large trees >50 cm dbh, their relative recruitment over five decades was highest in partly salvaged stand (91%), while unmanaged stand and selection stand followed with 72% and 68%, respectively. In all studied stands, conifers had higher recruitment rates than beech (for details see Appendix A). The increase in the share of large beech trees in unmanaged stand was relatively high; however, the absolute share of beech in the species composition of this stand remained low (Table 2; Figure 4).
The share of sycamore maple (Acer pseudoplatanus L.), a more light-demanding species than fir, beech, and spruce, increased in selection stand from 1.5% to 3.6% and in partly salvaged stand from 0.3% to 0.7%. In unmanaged stand, the share of maple was lowest, amounting to only 0.2% of the total stem density throughout the whole study period (Table 2).
In contrast to the species shares based on stem density, the composition based on basal area showed a much higher share of conifers in unmanaged stand compared to partly salvaged stand and selection stand (Figure 4). The dynamics of basal area of large trees expressed with the LTIba index showed similar recruitment of this parameter between partly salvaged stand and unmanaged stand (0.94 and 0.98, respectively), whereas the index value for basal area of large trees in selection stand amounted to 0.60. However, on the species level, large fir trees in partly salvaged stand exhibited the highest value (LTIba = 6.12) due to an absolute basal area increase from 2.4 m2·ha−1 to 16.8 m2·ha−1 over the 50-year period (Figure 4). All LTIba values are provided in the Appendix A, Table A5.

3.3. Standing Dead Trees

With respect to standing dead trees in partly salvaged stand, there were a total of 12 snags·ha−1 in 1975 and 18 snags·ha−1 in 2005. Between the two periods, the number of snags with diameters between 20–40 cm increased most noticeably (Figure 5), which probably caused the two distributions to differ significantly (χ2 = 22.1, p = 0.0047). Despite this increase in snag density in partly salvaged stand over time, there was a significant difference in snag diameter distributions between this stand and adjacent unmanaged stand in the final period 2005/2010 (χ2 = 38.07, p = 0.0003). Unmanaged stand was characterized by a total of 68 snags·ha−1, of which large snags with dbh >50 cm amounted to 16 snags·ha−1. On the other hand, partly salvaged stand had on average 2 large snags·ha−1 (Figure 5). The basal area of snags in partly salvaged stand increased between the two periods from 0.9 m2·ha−1 to 1.4 m2·ha−1, whereas the basal area of snags in unmanaged stand was noticeably higher, amounting to 7.9 m2·ha−1.

4. Discussion

4.1. Structure of Living Trees

The enhancement of old-growth attributes over large forested landscapes is currently the subject of debate. While there are many arguments for strict protection (e.g., [17,47]), a number of active silvicultural techniques have been recently suggested to accelerate large tree recruitment [4,48] or to naturally regenerate and preserve certain tree species [49,50]. Our study contributes to previous research by showing that partial salvage logging could also serve as a useful technique for restoring old-growth attributes such as large-sized trees and snags while maintaining tree-size demographic equilibrium. Under partial salvage logging, the negative exponential diameter distribution was maintained for two decades (1955–1975) before eventually being transformed into the increasing-q form in 2005. While this form occasionally led to unimodal structure, e.g., beech, in unmanaged stand and partly salvaged stand, it more often transformed into more sustainable rotated sigmoid or negative exponential shapes, as shown in the examples from selection stand and unmanaged stand.
Although our dataset did not allow us to investigate small-scale spatial interactions among tree species, other studies from beech-conifer forests have reported that the recruitment of young trees decreases when the share of beech increases in total basal area or total growing stock [25,51,52]. Given that the share of beech was highest in partly salvaged stand, this may explain why partial salvaging failed to maintain the negative exponential distribution or establish the rotated sigmoid distribution, whereas the cumulative dbh structure of adjacent unmanaged stand, even under strict protection, changed from the increasing-q into the rotated sigmoid form. It is also important to note that very high total basal area in unmanaged stand did not impede the maintenance of the rotated sigmoid structure. Similar outcomes were reported for the old-growth forest Janj in Bosnia-Herzegovina [25]. Thus, partial salvaging combined with removal of some large beech trees could result in more large fir and spruce trees making up a greater share of total basal area, which would create conditions for establishing the rotated sigmoid structure. Such an approach might be a sound solution when managing for old-growth conditions [4,7], at least in mixed beech-fir-spruce stands.
Negative exponential dbh distribution usually transforms into rotated sigmoid shape in selection forests as a result of the increase in basal area of large trees (see [53]). This was exactly what happened in selection stand in Tara mountain, primarily due to the increase in basal area of conifers by 2010. Although initially the negative exponential form was one of the management goals in selection stand [28], the outcome of applied selection system was not a targeted negative exponential form but actually the rotated sigmoid form. This result corroborates the statements by Schütz [54] that negative exponential form is difficult to establish and maintain. In addition, the rotated sigmoid shape proved to be sustainable for a longer period than the negative exponential form. Therefore, the rotated sigmoid distribution may be considered as a sound alternative to the negative exponential structure in the context of selection management. Schütz [54] indicated that the rotated sigmoid shape occurred in Swiss selection forests because medium-sized trees experienced accelerated diameter growth. In selection stand in Tara, however, this form occurred as a consequence of a slight decrease in medium-sized trees (probably due to cutting) and simultaneous increase in the recruitment of large-sized trees >50 cm in dbh. Thus, tree-size demographic equilibrium was not put at risk, as the rotated sigmoid form proved to be sustainable over five decades (in the case of fir in selection stand). Our results hence confirm the findings from northern hardwoods that demonstrated similar elasticity in the rotated sigmoid structure [42,55].
Although the negative exponential form may have certain limitations, it may be applicable when restoring some old-growth attributes if the target diameter is increased in mixed beech-fir-spruce forests. This is at least one of the conclusions from recent research in the old-growth forest Lom that is also located in the Dinaric Mountains (see [56]). On the other hand, the present study indicated that besides the rotated sigmoid and negative exponential structural forms, alternative options for management towards higher degree of old-growthness could also be pursued in the increasing-q form, which eventually occurred in partly salvaged stand. Similar outcomes were obtained by Gronewold et al. [16] in the northern hardwoods of Michigan (USA). Moreover, the same authors reported that the increasing-q structure may be sustained over time under single-tree selection. Several other studies have also determined the increasing-q shape in managed and unmanaged northern hardwoods, although it was less frequently reported than rotated sigmoid and negative exponential distributions [4,42,55,57].
In our study, the increasing-q shape was sustained for a maximum of two decades, which in certain circumstances may be considered as a drawback. Nevertheless, switching between rotated sigmoid, negative exponential, and increasing-q shapes depends on tree species composition, and the structural form may be changed in a short period of time (two to three decades) even by the application of only partial salvaging. Therefore, any of these distribution forms may be acceptable when managing towards old-growth conditions, at least in mixed beech-fir-spruce forests.
Some authors have already proposed the modification of conventional selection systems to promote old-growth attributes for northern hardwoods (e.g., [4,16]). We believe such modifications might be also possible in beech-fir-spruce stands in South-East Europe. Namely, our results showed a decrease in stem density and an increase in basal area in selection stand without disruption of demographic equilibrium over a period of five decades. Similar results were reported by O’Hara et al. [58] for Swiss plenter forests composed of the same tree species. This confirms the statement that high wood productivity is rather stable in selection stands, even when tree density is reduced [59]. Consequently, modifying selection systems towards a higher degree of old-growthness might be feasible in practice. Please note, however, that our case study did not include the replication of stands, and additional studies are necessary to extrapolate the outcomes to the broader scale. Thus, the question of the acceptance of alternative silvicultural techniques on a wider scale goes beyond the scope of our study and remains a challenge for future research and forest management.

4.2. The Influence of Management Regimes on Tree Species Composition

The study showed that different management regimes may create similar structural forms, but that the effects on tree species composition were varied. From a statistical viewpoint, both selection stand and partly salvaged stand had significantly different tree species compositions than adjacent unmanaged stand. Similar outcomes were reported when comparing old-growth forests and managed forests in Bosnia-Herzegovina [25]. This indicates that, despite climatic changes, management practices still have a large influence on the structural and compositional dynamics of mixed beech-fir-spruce stands.
With respect to tree species replacement patterns, a number of studies [12,60,61,62,63] have reported long-term fir/spruce decline and beech progression in European old-growth forests, which have sometimes been ascribed to the effects of climate change. Likewise, in managed forests across Europe, where beech was practically eliminated in 19th and 20th century, we now face a boomerang effect as beech naturally fights to reoccupy its areas. The present study revealed, however, that forest management still has a crucial impact on tree species replacement patterns, at least with respect to mixed beech-fir-spruce stands. The results also showed that in partly salvaged stand and unmanaged stand in Tara, beech is in decline and silver fir dominates all stand layers. However, this result should be treated with caution as we do not have the detailed information on possible human activities in the unmanaged stand prior to official strict protection in 1957. It is yet important to note that Petritan et al. [52] also reported successful recruitment of silver fir in old-growth forest Sinca in Romania, so further studies are needed to check if these are the isolated cases.
Given that every silvicultural approach shortens the lifespan of trees, our results support the conclusions of Vrška et al. [64] on the importance of cyclic replacement among tree species. Although fir and spruce are sensitive species in Serbia because they are close to the southern edge of their natural range, it seems that beech-fir-spruce stands are still rather stable compositionally (see also [56]). This finding may be important for the future management of these mixed forests; however, long-term research is needed to clarify the influence of climate change on each constituent species (see [65]).

4.3. Standing Dead Trees in Partly Salvaged Stand and Unmanaged Stand in Tara

Total snag density in unmanaged stand Tara was higher compared to other Dinaric old-growth forests as well as old-growth forests in northern hardwoods [55], but a common characteristic of all Dinaric old-growth forests is the similar density of large snags >50 cm dbh (16 ha−1–19 ha−1) (e.g., [66,67]). In certain situations, the discrepancy in snag densities between forests managed with conventional silvicultural systems and unmanaged references might not be evident [16]. Nonetheless, significant differences between managed and unmanaged forests are typically reported [68]. Although we observed the effects of an unconventional partial salvaging approach, the stem density and basal area of snags in partly salvaged stand were also significantly lower compared to the unmanaged stand.
Despite this, partial salvage logging led to a noticeable increase in snag density between 1975 and 2005. The largest increase in snags occurred in the 20 cm–40 cm dbh range, but large snags are also of special interest as they contain a number of microhabitats [7]. So far, however, few studies have proposed threshold values for large snag density in managed European stands. While some authors have suggested five snags with dbh >30 cm per hectare (e.g., [69]), others have proposed up to nine snags with dbh >43 cm per hectare [68]. However, the conventional selection systems that are usually applied in montane beech-conifer stands usually do not achieve these values [66,68]. In our case, even partly salvaged stand has heretofore failed to meet the levels proposed above. Nevertheless, considering the significant recruitment of live veteran trees in the last few decades under partial salvaging, there is obviously great potential to increase the density of large snags in partly salvaged stand in the near future.

5. Conclusions

Different management regimes resulted in significant differences in tree species composition between studied managed and unmanaged stands over a period of five decades. Contrary to a number of recent studies that reported the issues with the recruitment of young fir trees in European old-growth forests, in the unmanaged stand in Tara there was actually a problem with beech recruitment. However, this result should be treated with caution, as there is no information on possible human activities before strict protection in 1957. If such activities took place, they certainly influenced cyclic replacement trajectory among tree species (see [64]). Similar outcomes regarding fir and beech recruitment were determined for selection stand in Tara, however, the influence of forestry in this managed stand clearly showed that silvicultural practices still have major influence on tree species composition despite climate change.
On the other hand, despite the greater share of large-size trees in unmanaged stand, the diameter distributions of studied managed and unmanaged stands had the same structural forms in certain time periods. Since increasing-q and negative exponential forms most often transformed into rotated sigmoid form over time and with increase of basal area, the application of rotated sigmoid structure seems to be most practical if the management goal is to increase the share of large living trees.
Considering very low share of old-growth forests in Europe and worldwide, the promotion of old-growth attributes should be pursued. To this end, significant progress could be made if partial salvaging (instead of total salvaging) were to be included in conventional selection systems, at least on designated land units. Contrary to total salvaging, the study showed that partial salvaging may serve as a good alternative to promoting old-growth attributes such as large veteran trees and a small but valuable share of large snags. Forest managers could partly tolerate very large snags as their limited presence did not endanger demographic equilibrium of living trees in studied beech-fir spruce stands. Nonetheless, additional studies are necessary to investigate the full range of effects of partial salvage logging.

Author Contributions

Data curation, all authors; Formal analysis, S.K.; Investigation, S.K., M.M. and S.O.; Methodology, S.K.; Project administration, M.M.; Writing-original draft, S.K.; Writing-review & editing, all authors.

Funding

The research was funded by the Ministry of Education, Science and Technological Development (Serbia) in the frame of the project “Sustainable management of the total forest potentials in the Republic of Serbia”-EVBR 37008. We are also grateful to the Ministry of Science and Higher Education (MNiSW, Poland) for supporting open access publishing and the English proofreading of the article in the frame of Statutory funds: DS-3418/ZBiPL/2018.

Acknowledgments

For continuing cooperation, the authors would like to express special thanks to the administration of the National Park Tara.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table
Table A1. Total numbers of trees for which diameters at breast height (dbh, 1.30 m above ground) were measured during every inventory (complete census) in partly salvaged stand (area 8.0), selection-managed stand (area 11.2 ha), and unmanaged stand (area 15.4 ha).
Table A1. Total numbers of trees for which diameters at breast height (dbh, 1.30 m above ground) were measured during every inventory (complete census) in partly salvaged stand (area 8.0), selection-managed stand (area 11.2 ha), and unmanaged stand (area 15.4 ha).
Partly Salvaged Stand
Inventory year195519752005
Total number of measured trees410444403296
Selection Stand
Inventory year196019802010
Total number of measured trees556756225298
Unmanaged Stand
Inventory year196019802010
Total number of measured trees546763146591
Table
Table A2. Stem density and basal area (BA) for different size categories in partly salvaged-logged stand in Tara Mountain across five decades.
Table A2. Stem density and basal area (BA) for different size categories in partly salvaged-logged stand in Tara Mountain across five decades.
Partly Salvaged Stand
Beech 1955Beech 1975Beech 2005
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees200.35130.2350.09
small trees211.04150.7160.31
medium trees354.46323.99202.72
large trees3713.003713.893212.07
Ʃ11318.859618.826315.19
Fir 1955Fir 1975Fir 2005
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees2103.702053.631282.26
small trees844.13994.84623.05
medium trees657.419110.878110.12
large trees92.36277.595216.81
Ʃ36717.6042226.9332332.24
Spruce 1955Spruce 1975Spruce 2005
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees190.34160.2890.16
small trees70.3590.4540.18
medium trees60.6881.0160.69
large trees00.0020.4531.05
Ʃ321.36352.20222.09
Total trees 1955Total trees 1975Total trees 2005
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees2494.392344.141432.53
small trees1135.521236.02733.57
medium trees10612.6213115.9110713.58
large trees4615.566622.178830.18
Ʃ51338.0955548.2341249.86
Table
Table A3. Stem density and basal area (BA) for different size categories in selection managed stand in Tara Mountain across five decades.
Table A3. Stem density and basal area (BA) for different size categories in selection managed stand in Tara Mountain across five decades.
Selection Stand
Beech 1960Beech 1980Beech 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees250.43320.57470.83
small trees190.94150.75140.67
medium trees263.15222.77172.13
large trees196.43134.17134.23
Ʃ8810.95838.26917.86
Fir 1960Fir 1980Fir 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees1853.261672.961542.72
small trees783.81743.62492.41
medium trees728.55819.76597.36
large trees113.06246.57329.55
Ʃ34618.6834622.9029422.03
Spruce 1960Spruce 1980Spruce 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees230.41230.41270.48
small trees130.63130.62110.54
medium trees131.50161.97141.81
large trees20.4551.2592.80
Ʃ502.98574.25625.63
Total trees 1960Total trees 1980Total trees 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees2384.192334.122434.30
small trees1125.481045.12793.87
medium trees11413.6112214.849411.70
large trees3410.694312.435717.09
Ʃ49733.9750236.5047336.96
Table
Table A4. Stem density and basal area (BA) for different size categories in unmanaged stand in Tara Mountain across five decades.
Table A4. Stem density and basal area (BA) for different size categories in unmanaged stand in Tara Mountain across five decades.
Unmanaged Stand
Beech 1960Beech 1980Beech 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees90.1640.0730.05
small trees150.75120.5770.32
medium trees151.70182.18141.84
large trees82.56124.21186.30
Ʃ475.17457.03428.52
Fir 1960Fir 1980Fir 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees691.21831.471202.12
small trees432.09442.14311.54
medium trees658.24587.28475.93
large trees288.324514.155016.64
Ʃ20419.8622925.0424826.23
Spruce 1960Spruce 1980Spruce 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees70.12460.81480.86
small trees90.46110.52150.76
medium trees456.18263.63212.67
large trees308.524213.124716.28
Ʃ9215.2812518.0813220.57
Total trees 1960Total trees 1980Total trees 2010
Size categorystems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)stems ha−1BA (m2 ha−1)
pole trees851.501342.371723.03
small trees683.35663.24532.62
medium trees13317.2210713.868410.72
large trees6920.3710332.7611940.34
Ʃ35542.4341052.2242856.71
Table
Table A5. The values of the LTIba index for the assessment of basal area dynamics of large trees >50 cm dbh, calculated for each tree species and for all species combined at stand level.
Table A5. The values of the LTIba index for the assessment of basal area dynamics of large trees >50 cm dbh, calculated for each tree species and for all species combined at stand level.
LTIba Index Values
Partly Salvaged StandSelection StandUnmanaged Stand
Beech−0.07−0.341.46
Fir6.122.121.00
Sprucena5.190.91
All species0.940.600.98
na—not available in case when basal area of large trees >50 cm dbh was equal to 0 in the first inventory.

References

  1. Davis, M.B. Extent and location. In Eastern Old-Growth Forests: Prospects for Rediscovery and Recovery; Davis, M.B., Ed.; Island Press: Washington, DC, USA, 1996; pp. 18–36. [Google Scholar]
  2. Parviainen, J.; Bücking, W.; Vandekerkhove, K.; Schuck, A.; Päivinen, R. Strict forest reserves in Europe: Efforts to enhance biodiversity and research on forests left for free development in Europe (EU-COST-Action E4). Forestry 2000, 73, 107–118. [Google Scholar] [CrossRef]
  3. Global Forest Resource Assessment 2010. Available online: http://www.fao.org/docrep/013/i1757e/i1757e.pdf (accessed on 17 May 2018).
  4. Keeton, W.S. Managing for late-successional/old-growth characteristics in northern hardwood-conifer forests. For. Ecol. Manag. 2006, 235, 129–142. [Google Scholar] [CrossRef]
  5. Côté, P.; Tittler, R.; Messier, C.; Kneeshaw, D.D.; Fall, A.; Fortin, M.J. Comparing different forest zoning options for landscape-scale management of the boreal forest: Possible benefits of the TRIAD. For. Ecol. Manag. 2010, 259, 418–427. [Google Scholar] [CrossRef] [Green Version]
  6. Schütz, J.P.; Saniga, M.; Diaci, J.; Vrška, T. Comparing close-to-nature silviculture with processes in pristine forests: Lessons from Central Europe. Ann. For. Sci. 2016, 73, 911–921. [Google Scholar] [CrossRef]
  7. Bauhus, J.; Puettmann, K.; Messier, C. Silviculture for old-growth attributes. For. Ecol. Manag. 2009, 258, 525–537. [Google Scholar] [CrossRef] [Green Version]
  8. Ford, S.E.; Keeton, W.S. Enhanced carbon storage through management for old-growth characteristics in northern hardwood-conifer forests. Ecosphere 2017, 8, e01721. [Google Scholar] [CrossRef] [Green Version]
  9. Kimmins, J.P.; Rempel, R.; Welham, C.; Seely, B.; Rees, V.K. Biophysical sustainability, process-based monitoring and forest ecosystem management decision support systems. For. Chron. 2007, 83, 502–514. [Google Scholar] [CrossRef] [Green Version]
  10. Paillet, Y.; Pernot, C.; Boulanger, V.; Debaive, N.; Fuhr, M.; Gilg, O.; Gosselin, F. Quantifying the recovery of old-growth attributes in forest reserves: A first reference for France. For. Ecol. Manag. 2015, 346, 51–64. [Google Scholar] [CrossRef]
  11. Alessandrini, A.; Biondi, F.; Di, A.; Ziaco, E.; Piovesan, G. Tree size distribution at increasing spatial scales converges to the rotated sigmoid curve in two old-growth beech stands of the Italian Apennines. For. Ecol. Manag. 2011, 262, 1950–1962. [Google Scholar] [CrossRef]
  12. Diaci, J.; Rozenbergar, D.; Anic, I.; Mikac, S.; Saniga, M.; Kucbel, S.; Visnjic, C.; Ballian, D. Structural dynamics and synchronous silver fir decline in mixed old-growth mountain forests in Eastern and Southeastern Europe. Forestry 2011, 84, 479–491. [Google Scholar] [CrossRef] [Green Version]
  13. Motta, R.; Garbarino, M.; Berretti, R.; Meloni, F.; Nosenzo, A.; Vacchiano, G. Development of old-growth characteristics in uneven-aged forests of the Italian Alps. Eur. J. For. Res. 2015, 134, 19–31. [Google Scholar] [CrossRef]
  14. Dieler, J.; Uhl, E.; Biber, P.; Müller, J.; Rötzer, T.; Pretzsch, H. Effect of forest stand management on species composition, structural diversity, and productivity in the temperate zone of Europe. Eur. J. For. Res. 2017, 136, 739–766. [Google Scholar] [CrossRef]
  15. Spies, T.A.; Franklin, J.F. The Diversity and Maintenance of Old-Growth Forests. In Biodiversity in Managed Landscapes: Theory and Practice; Szaro, R.C., Johnston, D., Eds.; Oxford University Press: New York, NY, USA, 1996; pp. 296–314. ISBN 0-19-507958-2. [Google Scholar]
  16. Gronewold, C.A.; D’Amato, A.W.; Palik, B.J. The influence of cutting cycle and stocking level on the structure and composition of managed old-growth northern hardwoods. For. Ecol. Manag. 2010, 259, 1151–1160. [Google Scholar] [CrossRef]
  17. Hanson, J.J.; Lorimer, C.G.; Halpin, C.R.; Palik, B.J. Ecological forestry in an uneven-aged, late-successional forest: Simulated effects of contrasting treatments on structure and yield. For. Ecol. Manag. 2012, 270, 94–107. [Google Scholar] [CrossRef]
  18. Curtis, R.; Marshall, D.D. Douglas-Fir Rotations-Time for Reappraisal? West. J. Appl. For. 1993, 8, 81–85. [Google Scholar]
  19. Silver, E.J.; D’Amato, A.W.; Fraver, S.; Palik, B.J.; Bradford, J.B. Structure and development of old-growth, unmanaged second-growth, and extended rotation Pinus resinosa forests in Minnesota, USA. For. Ecol. Manag. 2013, 291, 110–118. [Google Scholar] [CrossRef]
  20. Bengtsson, J.; Nilsson, S.G.; Franc, A.; Menozzi, P. Biodiversity, disturbances, ecosystem function and management of european forests. For. Ecol. Manag. 2000, 132, 39–50. [Google Scholar] [CrossRef]
  21. Pretzsch, H. Diversity and productivity in forests: Evidence from long-term experimental plots. In Forest Diversity and Function: Temperate and Boreal Systems; Scherer-Lorenzen, M., Korner, C., Schulze, E.-D., Eds.; Springer: Berlin, Germany, 2005; Volume 176, pp. 41–64. ISBN 978-3-540-22191-3. [Google Scholar]
  22. Fischer, C.; Mölder, A. Trend to increasing structural diversity in German forests: Results from National Forest Inventories 2002 and 2012. Ann. For. Sci. 2017, 74, 80. [Google Scholar] [CrossRef]
  23. Spathelf, P.; Bolte, A.; van der Maaten, E. Is Close-to-Nature Silviculture (CNS) an adequate concept to adapt forests to climate change? Landbauforsch. App. Agric. For. Res. 2015, 65, 161–170. [Google Scholar] [CrossRef]
  24. Bončina, A. History, current status and future prospects of uneven-aged forest management in the Dinaric region: An overview. Forestry 2011, 84, 467–478. [Google Scholar] [CrossRef]
  25. Keren, S.; Diaci, J.; Motta, R.; Govedar, Z. Stand structural complexity of mixed old-growth and adjacent selection forests in the Dinaric Mountains of Bosnia and Herzegovina. For. Ecol. Manag. 2017, 400, 531–541. [Google Scholar] [CrossRef]
  26. Bosela, M.; Tobin, B.; Šeben, V.; Petráš, R.; Larocque, G.R. Different mixtures of Norway spruce, silver fir, and European beech modify competitive interactions in central European mature mixed forests. Can. J. For. Res. 2015, 45, 1577–1586. [Google Scholar] [CrossRef]
  27. Bankovic, S.; Medarevic, M.; Pantic, D.; Petrovic, N.; Sljukic, B.; Obradovic, S. The growing stock of the Republic of Serbia: State and problems. Bull. Fac. For. 2009, 100, 7–30. [Google Scholar] [CrossRef]
  28. Medarević, M.; Banković, S.; Pantić, D.; Obradović, S. Effects of the control method (Goč variety) in selection forest management in western Serbia. Arch. Biol. Sci. 2010, 62, 407–418. [Google Scholar] [CrossRef]
  29. Diaci, J.; Rozenbergar, D.; Fidej, G.; Nagel, T.A. Challenges for uneven-aged silviculture in restoration of post-disturbance forests in central Europe: A synthesis. Forests 2017, 8, 378. [Google Scholar] [CrossRef]
  30. Lang, K.D.; Schulte, L.A.; Guntenspergen, G.R. Windthrow and salvage logging in an old-growth hemlock-northern hardwoods forest. For. Ecol. Manag. 2009, 259, 56–64. [Google Scholar] [CrossRef] [Green Version]
  31. Fidej, G.; Rozman, A.; Diaci, J. Drivers of regeneration dynamics following salvage logging and different silvicultural treatments in windthrow areas in Slovenia. For. Ecol. Manag. 2018, 409, 378–389. [Google Scholar] [CrossRef]
  32. Thorn, S.; Bässler, C.; Svoboda, M.; Müller, J. Effects of natural disturbances and salvage logging on biodiversity—Lessons from the Bohemian Forest. For. Ecol. Manag. 2017, 388, 113–119. [Google Scholar] [CrossRef]
  33. Thorn, S.; Bässler, C.; Brandl, R.; Burton, P.J.; Cahall, R.; Campbell, J.L.; Castro, J.; Choi, C.Y.; Cobb, T.; Donato, D.C.; et al. Impacts of salvage logging on biodiversity: A meta-analysis. J. Appl. Ecol. 2018, 55, 279–289. [Google Scholar] [CrossRef] [PubMed]
  34. Nagel, T.A.; Mikac, S.; Dolinar, M.; Klopcic, M.; Keren, S.; Svoboda, M.; Diaci, J.; Boncina, A.; Paulic, V. The natural disturbance regime in forests of the Dinaric Mountains: A synthesis of evidence. For. Ecol. Manag. 2017, 388, 29–42. [Google Scholar] [CrossRef]
  35. Klopčič, M.; Simončič, T.; Bončina, A. Comparison of regeneration and recruitment of shade-tolerant and light-demanding tree species in mixed uneven-aged forests: Experiences from the Dinaric region. Forestry 2015, 88, 552–563. [Google Scholar] [CrossRef]
  36. Klopčič, M.; Mina, M.; Bugmann, H.; Bončina, A. The prospects of silver fir (Abies alba Mill.) and Norway spruce (Picea abies (L.) Karst) in mixed mountain forests under various management strategies, climate change and high browsing pressure. Eur. J. For. Res. 2017, 136, 1071–1090. [Google Scholar] [CrossRef]
  37. Vasić, M. Tara planina, Šumarska Studija; N. štamparija „Davidović“: Belgrade, Serbia, 1908. [Google Scholar]
  38. Milojković, D. Razvoj gazdovanja i uređivanja šuma GJ Tara. Šumarstvo 1986, 1–2, 11–27. [Google Scholar]
  39. Obradović, S. Stanje i Razvoj Sastojina Bukve, jele i Smrče Prašumskog Porekla u Srbiji kao Osnov za Planiranje i obezbeđivanje Prirodi Bliskog Gazdovanja. Ph.D. Thesis, University of Belgrade, Belgrade, Serbia, 2015. [Google Scholar]
  40. Pantic, D.; Medarevic, M.; Bankovic, S.; Obradovic, S.; Sljukic, B.; Pesic, B. Structural, production and dynamic characteristics of the strict forest reserve “Račanska šljivovica” on mt. Tara. Bull. Fac. For. 2011, 93–114. [Google Scholar] [CrossRef]
  41. Bončina, A.; Čavlović, J.; Čurović, M.; Govedar, Z.; Klopčič, M.; Medarević, M. A comparative analysis of recent changes in Dinaric uneven-aged forests of the NW Balkans. Forestry 2014, 87, 71–84. [Google Scholar] [CrossRef]
  42. Janowiak, M.K.; Nagel, L.M.; Webster, C.R. Spatial Scale and Stand Structure in Northern Hardwood Forests: Implications for Quantifying Diameter Distributions. For. Sci. 2008, 54, 497–506. [Google Scholar]
  43. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach; Springer: New York, NY, USA, 2002; ISBN 0-387-95364-7. [Google Scholar]
  44. Zenner, E.K.; Peck, J.L.E. Floating neighborhoods reveal contribution of individual trees to high sub-stand scale heterogeneity. For. Ecol. Manag. 2018, 412, 29–40. [Google Scholar] [CrossRef]
  45. Lafond, V.; Cordonnier, T.; De Coligny, F.; Courbaud, B. Reconstructing harvesting diameter distribution from aggregate data. Ann. For. Sci. 2012, 69, 235–243. [Google Scholar] [CrossRef]
  46. Klopcic, M.; Boncina, A. Stand dynamics of silver fir (Abies alba Mill.)-European beech (Fagus sylvatica L.) forests during the past century: A decline of silver fir? Forestry 2011, 84, 259–271. [Google Scholar] [CrossRef]
  47. Jaroszewicz, B.; Bobiec, A.; Eycott, A.E. Lack of demographic equilibrium indicates natural, large-scale forest dynamics, not a problematic forest conservation policy—A reply to Brzeziecki et al. J. Veg. Sci. 2017, 28, 218–222. [Google Scholar] [CrossRef]
  48. Fassnacht, K.S.; Bronson, D.R.; Palik, B.J.; Amato, A.W.D.; Lorimer, C.G.; Martin, K.J. Accelerating the Development of Old-Growth Characteristics in Second-Growth Northern Hardwoods; General Technical Report; U.S. Depoartment of Agriculture, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2015.
  49. Brzeziecki, B.; Pommerening, A.; Miścicki, S.; Drozdowski, S.; Zybura, H. A common lack of demographic equilibrium among tree species in Bialowieża National Park (NE Poland): Evidence from long-term plots. J. Veg. Sci. 2016, 27, 460–469. [Google Scholar] [CrossRef]
  50. Dobrowolska, D.; Niemczyk, M.; Olszowska, G. The influence of stand structure on European yew Taxus baccata populations in its natural habitats in central Poland. Pol. J. Ecol. 2017, 65, 369–384. [Google Scholar] [CrossRef]
  51. Čavlović, J.; Božić, M.; Bončina, A. Stand structure of an uneven-aged fir-beech forest with an irregular diameter structure: Modeling the development of the Belevine forest, Croatia. Eur. J. For. Res. 2006, 125, 325–333. [Google Scholar] [CrossRef]
  52. Petritan, I.C.; Commarmot, B.; Hobi, M.L.; Petritan, A.M.; Bigler, C.; Abrudan, I.V.; Rigling, A. Structural patterns of beech and silver fir suggest stability and resilience of the virgin forest Sinca in the Southern Carpathians, Romania. For. Ecol. Manag. 2015, 356, 184–195. [Google Scholar] [CrossRef]
  53. O’Hara, K.L.; Gersonde, R.F. Stocking control concepts in uneven-aged silviculture. Forestry 2004, 77, 131–143. [Google Scholar] [CrossRef] [Green Version]
  54. Schütz, J.P. Die Plenterung und ihre Unterschiedlichen Formen Skript zu Vorlesung Waldbau II und Waldbau IV; ETH Zürich: Zürich, Switzerland, 2002. [Google Scholar]
  55. Goodburn, J.M.; Lorimer, C.G. Population structure in old-growth and managed northern hardwoods: An examination of the balanced diameter distribution concept. For. Ecol. Manag. 1999, 118, 11–29. [Google Scholar] [CrossRef]
  56. Motta, R.; Berretti, R.; Castagneri, D.; Dukić, V.; Garbarino, M.; Govedar, Z.; Lingua, E.; Maunaga, Z.; Meloni, F. Toward a definition of the range of variability of central European mixed Fagus–Abies–Picea forests: The nearly steady-state forest of Lom (Bosnia and Herzegovina). Can. J. For. Res. 2011, 41, 1871–1884. [Google Scholar] [CrossRef]
  57. Schwartz, J.W.; Nagel, L.M.; Webster, C.R. Effects of uneven-aged management on diameter distribution and species composition of northern hardwoods in Upper Michigan. For. Ecol. Manag. 2005, 211, 356–370. [Google Scholar] [CrossRef]
  58. O’Hara, K.L.; Hasenauer, H.; Kindermann, G. Sustainability in multi-aged stands: An analysis of long-term plenter systems. Forestry 2007, 80, 163–181. [Google Scholar] [CrossRef]
  59. Pretzsch, H.; Rötzer, T. Indicating Forest Ecosystem and Stand Productivity: From Deductive to Inductive Concepts. In Ecological Forest Management Handbook, Series of Applied Ecology and Environmental Management; Laroque, G.R., Ed.; CRC Press Taylor & Francis Group: Boca Raton, FL, USA, 2016; pp. 415–444. [Google Scholar]
  60. Jaloviar, P.; Saniga, M.; Kucbel, S.; Pittner, J.; Vencurik, J.; Dovciak, M. Seven decades of change in a European old-growth forest following a stand-replacing wind disturbance: A long-term case study. For. Ecol. Manag. 2017, 399, 197–205. [Google Scholar] [CrossRef]
  61. Keren, S.; Motta, R.; Govedar, Z.; Lucic, R.; Medarevic, M.; Diaci, J. Comparative structural dynamics of the Janj mixed old-growth mountain forest in Bosnia and Herzegovina: Are conifers in a long-term decline? Forests 2014, 5, 1243–1266. [Google Scholar] [CrossRef]
  62. Marinšek, A.; Diaci, J. A comparison of structural characteristics and ecological factors between forest reserves and managed silver fir: Norway spruce forests in Slovenia. Ekológia (Bratislava) 2011, 30, 51–66. [Google Scholar] [CrossRef]
  63. Szwagrzyk, J.; Szewczyk, J.; Maciejewski, Z. Shade-tolerant tree species from temperate forests differ in their competitive abilities: A case study from Roztocze, south-eastern Poland. For. Ecol. Manag. 2012, 282, 28–35. [Google Scholar] [CrossRef]
  64. Vrška, T.; Adam, D.; Hort, L.; Kolář, T.; Janík, D. 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? For. Ecol. Manag. 2009, 258, 347–356. [Google Scholar] [CrossRef]
  65. Dobrowolska, D.; Bončina, A.; Klumpp, R. Ecology and silviculture of silver fir (Abies alba Mill.): A review. J. For. Res. 2017, 22, 326–335. [Google Scholar] [CrossRef]
  66. Keren, S.; Diaci, J. Comparing the quantity and structure of deadwood in selection managed and old-growth forests in South-East Europe. Forests 2018, 9, 76. [Google Scholar] [CrossRef]
  67. Motta, R.; Garbarino, M.; Berretti, R.; Bjelanovic, I.; Borgogno Mondino, E.; Čurović, M.; Keren, S.; Meloni, F.; Nosenzo, A. Structure, spatio-temporal dynamics and disturbance regime of the mixed beech–silver fir–Norway spruce old-growth forest of Biogradska Gora (Montenegro). Plant Biosyst. 2015, 149, 966–975. [Google Scholar] [CrossRef]
  68. Marage, D.; Lemperiere, G. The management of snags: A comparison in managed and unmanaged ancient forests of the Southern French Alps. Ann. For. Sci. 2005, 62, 135–142. [Google Scholar] [CrossRef] [Green Version]
  69. Mason, F.; Nardi, G.; Whitmore, D. Recherches sur la restauration des habitats du bois mort: L’exemple du LIFE “Bosco della Fontana” (Italie). In Bois Mort et à Cavités, Une clé Pour Des Forêts Vivantes; Vallauri, D., André, J., Dodelin, B., Eynard-Machet, R., Rambaud, D., Eds.; Tec & Doc: Paris, France, 2005; pp. 285–291. [Google Scholar]
Forests 09 00479 g001 550
Figure 1. Diameter distributions of live trees with diameter at breast height (dbh) >10 cm in partly salvaged stand in Tara Mountain in 1955, 1975, and 2005 per species and for total trees.
Figure 1. Diameter distributions of live trees with diameter at breast height (dbh) >10 cm in partly salvaged stand in Tara Mountain in 1955, 1975, and 2005 per species and for total trees.
Forests 09 00479 g001
Forests 09 00479 g002 550
Figure 2. Diameter distributions of live trees with diameter at breast height (dbh) >10 cm in selection-managed stand in Tara Mountain in 1960, 1980, and 2010 per species and for total trees.
Figure 2. Diameter distributions of live trees with diameter at breast height (dbh) >10 cm in selection-managed stand in Tara Mountain in 1960, 1980, and 2010 per species and for total trees.
Forests 09 00479 g002
Forests 09 00479 g003 550
Figure 3. Diameter distributions of live trees with diameter at breast height (dbh) >10 cm in unmanaged stand in Tara Mountain in 1960, 1980, and 2010 per species and for total trees.
Figure 3. Diameter distributions of live trees with diameter at breast height (dbh) >10 cm in unmanaged stand in Tara Mountain in 1960, 1980, and 2010 per species and for total trees.
Forests 09 00479 g003
Forests 09 00479 g004 550
Figure 4. Species-specific and total basal area changes across time in partially salvaged stand, selection-managed stand, and unmanaged stand in the study area.
Figure 4. Species-specific and total basal area changes across time in partially salvaged stand, selection-managed stand, and unmanaged stand in the study area.
Forests 09 00479 g004
Forests 09 00479 g005 550
Figure 5. Snags in partly salvaged stand and unmanaged stand in the study area.
Figure 5. Snags in partly salvaged stand and unmanaged stand in the study area.
Forests 09 00479 g005
Table
Table 1. The summary of best-fitting models and corresponding shapes of diameter distributions per species and total trees in differently managed stands across time in the study area.
Table 1. The summary of best-fitting models and corresponding shapes of diameter distributions per species and total trees in differently managed stands across time in the study area.
Partly Salvaged Stand (1955)Selection Stand (1960)Unmanaged Stand (1960)
FirSpruceBeechTotalFirSpruceBeechTotalFirSpruceBeechTotal
RMSE0.140.060.090.070.030.030.060.080.070.130.090.06
Adjusted R20.970.980.960.990.980.990.990.940.990.910.950.98
p0.00000.00540.00000.00000.00120.00320.00000.00150.00000.00150.00000.0000
AIC−10.3−8.5−17.4−22.9−12.0−15.9−39.3−25.2−19.5−3.1−13.7−16.0
ShapeNENEIQNERSNEIQNEIQUNIIQIQ
Partly Salvaged Stand (1975)Selection Stand (1980)Unmanaged Stand (1980)
FirSpruceBeechTotalFirSpruceBeechTotalFirSpruceBeechTotal
RMSE0.080.070.050.050.060.030.060.070.060.240.130.05
Adjusted R20.990.980.980.990.980.990.980.930.980.800.860.98
p0.00000.00000.00040.00000.00010.00010.00000.00180.00010.0230.00040.0001
AIC−19.0−13.7−27.7−29.2−11.3−16.1−23.9−18.4−14.53.8−6.2−9.5
ShapeNENEUNINERSRSNENERSIQIQRS
Partly Salvaged Stand (2005)Selection Stand (2010)Unmanaged Stand (2010)
FirSpruceBeechTotalFirSpruceBeechTotalFirSpruceBeechTotal
RMSE0.110.080.070.090.030.060.090.090.100.070.060.06
Adjusted R20.980.940.950.980.990.990.970.980.950.970.950.99
p0.00000.00040.00000.00000.00010.00010.00050.00010.00170.00010.00000.0000
AIC−13.8−11.5−21.6−18.3−8.4−19.5−11.3−12.2−13.3−16.7−21.8−21.1
ShapeRSIQUNIIQRSRSRSRSRSRSUNIRS
RMSE—root mean square error; AIC—Akaike information criterion; Structural shapes: UNI—convex or unimodal shape, RS—rotated sigmoid, IQ—increasing-q, NE—negative exponential, (sensu Janowiak et al. [42]).
Table
Table 2. Long-term changes in tree species composition (%) of stands in Tara Mountain under different management regimes such as (i) partial salvaging, (ii) selection system, and (iii) strict protection in unmanaged stand.
Table 2. Long-term changes in tree species composition (%) of stands in Tara Mountain under different management regimes such as (i) partial salvaging, (ii) selection system, and (iii) strict protection in unmanaged stand.
Stem Density Per Hectare (%)
Stand CategoryYearBeech (%)Fir (%)Spruce (%)Maple (%)Other sp. (%)
Partly Salvaged Stand195521.971.56.30.30.0
197517.476.06.40.20.0
200515.378.65.40.70.0
Selection Stand196017.869.510.11.51.1
198016.568.811.31.71.7
201019.262.213.13.61.9
Unmanaged Stand196013.257.325.90.23.4
198011.056.030.50.22.3
20109.758.130.90.21.1

Share and Cite

MDPI and ACS Style

Keren, S.; Medarević, M.; Obradović, S.; Zlokapa, B. Five Decades of Structural and Compositional Changes in Managed and Unmanaged Montane Stands: A Case Study from South-East Europe. Forests 2018, 9, 479. https://doi.org/10.3390/f9080479

AMA Style

Keren S, Medarević M, Obradović S, Zlokapa B. Five Decades of Structural and Compositional Changes in Managed and Unmanaged Montane Stands: A Case Study from South-East Europe. Forests. 2018; 9(8):479. https://doi.org/10.3390/f9080479

Chicago/Turabian Style

Keren, Srđan, Milan Medarević, Snežana Obradović, and Brane Zlokapa. 2018. "Five Decades of Structural and Compositional Changes in Managed and Unmanaged Montane Stands: A Case Study from South-East Europe" Forests 9, no. 8: 479. https://doi.org/10.3390/f9080479

APA Style

Keren, S., Medarević, M., Obradović, S., & Zlokapa, B. (2018). Five Decades of Structural and Compositional Changes in Managed and Unmanaged Montane Stands: A Case Study from South-East Europe. Forests, 9(8), 479. https://doi.org/10.3390/f9080479

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop