Next Article in Journal
Mapping Forest Biomass Using Remote Sensing and National Forest Inventory in China
Previous Article in Journal
Vulnerability of Plantation Carbon Stocks to Defoliation under Current and Future Climates
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 5
Issue 6
10.3390/f5061243
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

Comparative Structural Dynamics of the Janj Mixed Old-Growth Mountain Forest in Bosnia and Herzegovina: Are Conifers in a Long-Term Decline?

by
Srdjan Keren
1,*,
Renzo Motta
2,
Zoran Govedar
1,
Radovan Lucic
3,
Milan Medarevic
4 and
Jurij Diaci
5
1
Faculty of Forestry, University in Banja Luka, Bulevar Stepa Stepanovic 75, Banja Luka 78000, Bosnia and Herzegovina
2
Department of Agriculture, Forestry and Food sciences, University of Turin, via Leonardo da Vinci, Grugliasco 44, Italy
3
Research and Project Center, Jovana Ducica 28, Banja Luka 78000, Bosnia and Herzegovina
4
Faculty of Forestry, University of Belgrade, Kneza Viseslava 1, Belgrade 11000, Serbia
5
Department of Forestry and Renewable Forest Resources, Biotechnical Faculty, University in Ljubljana, Vecna pot 83, Ljubljana 1000, Slovenia
*
Author to whom correspondence should be addressed.
Forests 2014, 5(6), 1243-1266; https://0-doi-org.brum.beds.ac.uk/10.3390/f5061243
Submission received: 16 April 2014 / Revised: 12 May 2014 / Accepted: 20 May 2014 / Published: 10 June 2014
Browse Figures
Versions Notes

Abstract

:
Regression of conifers in European mixed old-growth mountain forests has been observed for a long period and studied from different aspects. Old-growth (OG) forests in Bosnia and Herzegovina (BiH) have not experienced heavy air pollution and chronic overbrowsing that have affected many other European OG forests, while climatic and anthropogenic disturbances have been well documented. We analysed stand structure in the Janj OG forest, compared it with inventories of Lom and Perucica OG forests (BiH) and with earlier inventories of the same reserves. At present, OG forest Janj is characterized by a high growing stock (1215 m3∙ha−1). This is due to good site quality, prevalence of conifers (84%) and dominant endogenous processes in recent decades. In all three OG forests, indicators of structural change exhibited progression of European beech over time. Historical evidence revealed the occurrence of warm summers and droughts followed by bark beetle outbreaks in the 1920s, 1940s and early 1950s, which in turn influenced a marked conifer decline. It seems likely that repeated canopy opening released waves of European beech regeneration. These stand structural changes have delayed the rejuvenation of conifers and can help explain the early observations of conifer decline.

1. Introduction

Mixed forests composed of Fagus sylvatica L. (hereafter beech), Abies alba Mill. (fir) and Picea abies (L.) Karst. (Norway spruce) represent one of the most ecologically and economically important forest categories of European mountainous regions [1]. They grow in the relatively humid and cold climates of the mountains [2]. These forest sites have never been intensively used for agriculture or settlements and are still well preserved. Many European old-growth forests are situated in this vegetation belt [3]. However, there are many reports of synchronous regression of conifers, especially of silver fir in old-growth forests (e.g., [4,5,6,7]) as well as in managed forests [8,9,10,11]. Many factors, natural and anthropogenic, influence the coexistence of beech, fir and spruce, making this topic challenging to study [12,13,14,15].
It is well established that, among the main species of this altitudinal belt, spruce is the most tolerant to the cold continental climate, followed by fir, while beech may better tolerate a warmer and dryer climate [16]. Both conifers are more successful on sites with less competition from ground vegetation (rocky outcrops, CWD) and on silicate bedrock [2]. Fir is the strongest competitor in low light levels under dense stands, while the competitive strength of beech increases within canopy gaps [17,18,19] and that of spruce after intermediate disturbances [6]. Silver fir is very vulnerable to browsing [20] and polluted air, especially SO2 pollution [21]. Regarding species coexistence, not only are the direct effects of environmental changes important, such as higher temperatures or changes in the disturbance regime, but also indirect effects such as higher frequency of seed years in beech [22] or changes in ground vegetation structure due to nitrogen deposition and a prolonged growing season [23].
Neighbourhood effects [24] form an important group of coexistence factors. They include processes such as the changes in microsite conditions induced by the presence of certain canopy species. This variability may include light transmission [25], soil characteristics [4,26], nutrient availability [27,28] and biological factors and agents (e.g., allelopathy, species specific pathogens and predators) [29]. Such processes can lead toward reciprocal, random or self-replacement of different species in the canopy [12,30,31]. Thus, coexistence is driven by many factors, many of which are difficult to separate from the management regime in managed forests [3]. For this reason, the study of coexistence should involve research of old-growth forests. Since these are rare in Europe, comparative and integrative studies are important.
Most European old-growth forests are relatively small remnants surrounded by a buffer zone and located inside a matrix of managed forests. A common perception is that they have been well preserved from human impact. However, in past centuries rural areas were more densely inhabited than today and unregulated management that included felling of individual trees for fuel wood or construction on site, forest grazing and charcoal burning was common [7,9,15,32,33]. On the other hand, in the forest matrix surrounding old-growth forests, conifers were favoured for economic reasons. This enhanced their overall seed potential and likely positively influenced conifer regeneration in old-growth forests because of their limited size.
Moreover, in some regions deer became extinct, which additionally favoured regeneration of fir [14]. After WWII, deer densities in Europe sharply increased, and overbrowsing negatively influenced fir abundance [20]. Between the 1960s and 1980s, fir was heavily affected by air pollution [21]. Therefore, local direct historical influences and regional conditions (e.g., deer density, air pollution) are very important for understanding the recent structure of old-growth forests and their long-term dynamics [34].
Old-growth forests in Bosnia (Perucica, Lom and Janj) have not experienced chronic overbrowsing by deer as have other European old-growth forests [15,35], and environmental pollution has also been less acute [36]. These forests therefore represent an excellent resource for the study of interactions of climatic fluctuations and indirect anthropogenic disturbances, both of which are well documented in the literature [33,37,38,39,40,41]. The Janj old-growth forest has received less scientific attention [42,43] than the forests of Perucica (e.g., [5,19,44,45,46]) and Lom [32,47,48]. The aims of this study were (1) to describe the structural changes over a 60-year period in the Janj old-growth forest and (2) to compare and contrast the changes with those in the Lom and Perucica (Bosnia and Herzegovina) old-growth forests and other similar old-growth forests in Europe.

2. Materials and Methods

2.1. Study Site

The research was carried out in the Janj old-growth forest located on the mountainous massif of Stolovas near Sipovo, Bosnia and Herzegovina (BiH). The geographic position of the core area of Janj is 44°08'N, 17°17'E, and the altitude ranges from 1260–1400 m a.s.l. (Figure 1). The parent rock is made of dolomite bedrock on which deep rendzina and brown soil types of various depths are formed. Mean annual rainfall amounts to 1200 mm, and mean annual temperature is around 5 °C. The core area of old-growth forest at Janj has been classified as plant association Piceo-Abieti-Fagetum [49].
This virgin forest was protected by the decision of the State Institute for Protection of Cultural Monuments and Natural Rarities in Bosnia and Herzegovina (No. 245/54 from July 17, 1954). It is believed that even before the official protection of the Janj old-growth forest, this forest had not been exposed to organized logging activities. Unlike most European old-growth forests, Janj is a unique site for studying natural regeneration processes since previous studies have not reported the impact of ungulates as significant. The original purpose of the reserve was to serve primarily as an object for scientific research. Due to its inaccessibility, Janj has also remained undisturbed by human activities. However, for this reason combined with adverse historical circumstances (economic crisis, war, etc.), Janj has been little studied. In this old-growth forest there is a core area (57.2 ha; 1280–1400 m a.s.l.) surrounded by a buffer zone of 237.8 ha. Only low-intensity cuttings have been performed in the buffer zone, mostly in the form of salvage cuttings.
Forests 05 01243 g001 1024
Figure 1. Geographic location, buffer zone and core area with research plot network in the Janj old-growth forest. The insertion at the upper right-hand side of the picture represents the location of the Lom, Janj and Perucica old-growth forests in Bosnia and Herzegovina.
Figure 1. Geographic location, buffer zone and core area with research plot network in the Janj old-growth forest. The insertion at the upper right-hand side of the picture represents the location of the Lom, Janj and Perucica old-growth forests in Bosnia and Herzegovina.
Forests 05 01243 g001
Inventory data from Janj was compared with recent inventories of the Lom [47] and Perucica [50] old-growth forests. Perucica is located approx. 150 km southeast of Janj, while Lom is located about 80 km northwest (Figure 1). All three old-growth forests are situated approx. 90 km from the Adriatic coast. Lom (1250–1522 m a.s.l.), with a total protected area of 297.8 ha and a core area of 55.8 ha (1280–1350 m a.s.l.), is on prevalent limestone bedrock with numerous karst phenomena such as sinkholes and rock outcrops. The annual average precipitation is about 1600 mm, and the mean annual temperature is 7.8 °C (Drinic station, 730 m a.s.l.). Perucica (1000–1600 m a.s.l.), with total protected area of 1434.6 ha, is characterised by silicate parent material in the lower parts, while in the higher parts limestone is frequent [44,45]. The mean annual air temperature is between 5.9 and 8.6 °C and mean annual precipitation ranges from 1400–1900 mm (Suha and Cemerno station, 690–1329 m a.s.l., respectively).

2.2. Field Measurements

Measurements and recordings in Janj were carried out during summer 2011. To allow direct comparison with Lom, we followed the methodology described in Motta et al. [47]. A regular 100 m grid with 40 sampling points in the core area was superimposed (Figure 1). Each grid intersection defined the centre of a sampling plot where the following inventories were carried out: in a 452 m2 circular plot (radius = 12 m) species dbh (to the nearest 0.01 m) for all living trees above 7.5 cm was measured; heights of living trees were measured (to the nearest 0.5 m) for a sample of 100 trees for each species (spruce, fir and beech); in a 78.5 m2 circular plot (radius = 5 m) species and height of regeneration individuals from 10 cm height to 7.5 cm diameter at breast height (dbh) were tallied.
With respect to course woody debris (CWD), the following measurements were carried out: on two 50 m line intersects oriented northward (the first) and eastward (the second) from the centre of the sampling point, logs crossing the line were measured (see [51]); in two 50 × 4 m rectangular plots centred on the previous line, stumps (diameter at the ground and at the top) and snags (dbh) were measured. For each element of CWD, decay classes were recorded (class 1 fresh, class 5 very old) (see [52]); CWD was grouped into snags (standing dead trees, dbh ≥ 7.5 cm and height ≥ 1.3 m), downed logs (fallen stems or branches ≥ 7.5 cm diameter and length > 1 m) and stumps (short, vertical remains from cutting or windthrow, top diameter ≥ 7.5 cm and height < 1.30 m). The separation of snags from logs was established at a 45° leaning angle.
Data for the recent inventory of old-growth Lom originates from Motta et al. [47], while data for Perucica was measured in 2011 on 12 square sampling plots (0.50 ha) installed on a regular 1000 × 1000 m grid covering the whole core area of the reserve. Within plots, all trees with dbh ≥ 10 cm were tallied according to species and dbh was measured by Lucic [47]. The data for comparison of structural changes in the three old-growth forests originates from the study of Drinic [44]. In 1952/53, he installed 14 approximately one-hectare plots within the old-growth forests of Perucica (six plots), Janj (four plots) and Lom (four plots). On these plots, dbh and height of all trees with dbh ≥ 10 cm and some additional structural parameters were measured. Inventories in 1952/1953 were not carried out on the same plots as in the recent inventories; therefore, some bias may be expected. However, the ecological and terrain characteristics of the Lom and Janj old-growth forests are similar across the core area and Drinic’s sample size in each old-growth forest was large (ca. 4 ha). Therefore, we believe this first inventory provides good insight into the overall structural characteristics at that time. Perucica is larger and more variable. For all these reasons, special care was taken in the analyses to account for possible bias. Data was summarised and we focused only on the most persistent and clear trends.
For the Janj old-growth forest, data on dbh structure from the state forest inventory in 1965 published by Matic et al. [53] and the inventory in 2003 published by Maunaga et al. [43] were also available. Both studies followed official procedures for forest inventory in BiH, meaning that concentric circles were superimposed in a regular 100 × 100 m grid and the inventory threshold was set with dbh ≥ 10 cm.

2.3. Analyses

To determine site quality by tree species, we compared height curves from this study with standardized site quality curve lines for spruce, fir and beech valid in BiH [54]. For determining the average height of the upper stand story, we used the heights of the 20 largest diameter trees for each species. Volume of living trees and snags was calculated according to local volume tables. For estimating volume of logs, stumps and broken snags, the same method as described in Motta et al. [55] was applied.
The shape of diameter distributions was analysed using the methodology presented in Leak [56] and Janowiak et al. [57]. For cumulative dbh frequency distribution and for each tree species, a series of multiple regressions between the base 10 logarithm of trees per hectare (dependent variable) and all possible combinations of the dbh class midpoint, the midpoint squared and the midpoint cubed (independent variables) were calculated. The highest adjusted R2 and lowest root mean square error values were used as a basis for selection of the best-fitting model from all significant models (p < 0.05). Based on the sign of the coefficient, dbh distribution shapes could be classified into rotated sigmoid (RS), unimodal (UNI), negative exponential (NE), increasing-q (IQ) and concave (CO). If shapes were inconsistent or variable, the second best-fitting model was selected.
We studied neighbourhood effects based on transition probabilities between distinct stories of old-growth forest. We compared relative tree species abundance in each sampling plot in three layers: (1) regeneration from 10 cm in height to 7.5 cm dbh (understory); (2) poles with dbh from 7.5–27.5 cm (middle-story) and (3) canopy trees with dbh over 27.5 cm (upper-story). We created three transition matrices to avoid possible oversimplification of the data as observed by White et al. [58]. Transition matrices were based on the proportions of understory to upper-story trees, middle-story to upper-story trees and understory to middle-story trees. Further, we calculated relative abundances of tree species after 20 generations, after species composition was stabilized, using present abundances and transition probabilities as described by Stevens [59]. In these calculations, Sycamore maple (Acer pseudoplatanus L.) was not considered since it was not present in the middle or upper-story. Because of the many assumptions of these models, which are difficult to verify, and because we did not sample regeneration below each tree on sampling plots, we focused only on the most obvious relationships.
Changes in relative tree species composition of the three old-growth forests according to tree density and basal area in plots for the last six decades was tested with a t-test. Species proportions per plot were obtained from species counts, thus the arcsine transformation was applied to meet the assumptions of normality.

3. Results

3.1. Structure and Site Quality of the Janj Old-Growth Forest

In the Janj old-growth forest, the mean tree density was 516 ha−1 (Table 1). Density varied from 243–995 trees ha−1. Species composition in regard to the number of trees was strongly in favour of beech with 63%, followed by fir and spruce with 23% and 14%, respectively. The total number of trees with dbh over 50 cm amounted to 115 (fir 64 trees ha−1, spruce 37 trees ha−1, beech 14 trees ha−1). Only fir (27 trees ha−1) and spruce (15 trees ha−1) were present above 80 cm dbh.
Table
Table 1. Structural characteristics of the Janj old-growth forest (core area).
Table 1. Structural characteristics of the Janj old-growth forest (core area).
Mean diameter (cm)Density (trees ha−1)%Basal area (m2∙ha−1)%Volume (m3∙ha−1)
European beech22.53276311.718192
Silver fir60.61172334.752648
Norway spruce58.2721420.230376
Sycamore maple–<1<1–<1–
Total46.7516 66.7 1215
Range9.0–100.1243–995 26.1–102.5 419–1946
St. deviation23.0155 17.7 347
Basal area on sample plots varied from 26.1–102.5 m2. Average basal area was 34.7 m2, 20.2 m2 and 11.7 m2 for fir, spruce and beech, respectively. The share of sycamore maple in basal area was negligible although the site conditions suit this species well. With respect to basal area, fir dominated with 52%, followed by spruce and beech with 30% and 18%, respectively. Thus, species composition was significantly different when estimated according to the tree density and basal area.
The frequency distribution of diameters was dominated by beech within the small diameter classes, but its share sharply declined towards the medium and large diameter classes (Figure 2). Large diameter classes were dominated by fir and spruce. For 20 cm dbh class widths, the best-fitting multiple regression model of dbh frequency distribution indicated an UNI shape for fir, IQ shape for spruce, CO shape for beech, while for the cumulative distribution an RS shape was indicated (Table 5).
Forests 05 01243 g002 1024
Figure 2. Diameter distribution of living trees in the core area of Janj.
Figure 2. Diameter distribution of living trees in the core area of Janj.
Forests 05 01243 g002
The relation between dbh and height of trees suggested that site quality for all three species (beech, fir, spruce) could be assigned to the first class, indicating that these species grow in the most favourable conditions for their development. In younger stages, beech was taller than its competitors of the same diameter, while within large dbh classes (≥70 cm) spruce and fir were taller. The largest diameter and tallest individuals measured were firs with a dbh of 124 cm and a height of 50.5 m, respectively. Spruce reached 110 cm in dbh and 47.5 m in height, while beech grew to 92 cm in dbh and 44.7 m in height, respectively. There were no occurrences in any of the species of the largest diameter tree concurrently being the tallest. The mean values of height and corresponding dbh for the 20 largest diameter trees of each species are presented in Table A1.

3.2. Coarse Woody Debris

The total volume of CWD in Janj was 386.5 m3∙ha−1, which was 32% of the growing stock (Table 2). Most of the CWD volume was in logs, followed by snags and stumps with 79%, 19% and 2%, respectively. Logs were also more frequent than snags and stumps combined.
All five CWD decay classes were present only in logs, where the fifth decay class was prevalent. Most snags were in the first decay class, while the fifth was not present. Stumps were mostly in the fourth decay class. The first and fifth decay classes of stumps were completely absent in the core area (Table 3). On the other hand, large standing snags (≥50 cm dbh) were present with 19 trees ha−1, contributing to biodiversity and the survival of other species, especially nesting birds that depend on them. Interestingly, the number of smaller diameter snags (7.5–50 cm dbh) was lower, amounting to 17 trees ha−1. We also tested for relations between CWD volume and growing stock, but the results were not significant.
Table
Table 2. Volume of coarse woody debris types in the core area.
Table 2. Volume of coarse woody debris types in the core area.
Snags
(m3∙ha−1)		Logs
(m3∙ha−1)		Stumps
(m3∙ha−1)		Total
(m3∙ha−1)
Volume70.9306.78.8386.5
Range 0–510.2103.8–582.50–26.1162.1–892.6
SD107.8125.25.5169.7
Table
Table 3. Different decay classes of coarse woody debris in the core area.
Table 3. Different decay classes of coarse woody debris in the core area.
Decay classesSnags
(m3∙ha−1)		Logs
(m3∙ha−1)		Stumps
(m3∙ha−1)
120.12.20
216.230.30.5
318.864.23.0
415.888.85.3
50121.20

3.3. Replacement Patterns and Neighbourhood Effects

There were several indices that suggested a prevalent endogenous disturbance regime in the Janj old-growth forest in recent decades. This offered favourable conditions for investigation of neighbourhood effects of tree species. In the understory and middle-story, beech was by far the most frequent species, at 69% and 83%, respectively (Table 4). Fir was the second most abundant species in the understory and middle-story at 15% and 10% and prevailed in the upper-story at 48%. Between the understory and upper-story, beech was the most abundant species under its conspecific, while between the understory and middle-story, as well as between the middle-story and upper-story, it was more frequent under conifers. Fir indicated some signs of reciprocal replacement since in all three matrices it was more frequent under beech or spruce. Also, spruce was not the most frequent species under its conspecific in any of the comparisons. Models of future relative abundance indicated a strong increase of beech in the upper story and a moderate decrease in the middle story and vice versa for both conifers.

3.4. Structural Changes of the Janj, Lom and Perucica Old-Growth Forests over Six Decades

In Janj, the share of beech in tree density was similar in the first two inventories of 1952 and 1965, at 38% and 37%, respectively. It substantially increased to 54% and 63% in the last two inventories of 2003 and 2011, respectively. From 1952–2011, spruce decreased from 30%–14%, while the share of fir decreased less intensively from 32%–23%, respectively. A similar pattern of beech progression and conifer decline was indicated in Lom and Perucica (Figure 3). This trend was also comparable regarding the shares of tree species in relation to basal area. Here, the only exception was an increase in the share of fir in Janj, from 43%–52% from 1952–2011, respectively. The increase in the share of beech in density was statistically significant for all three old-growth forests, while its increase in basal area was only significant for Perucica.
Table
Table 4. Transition probabilities of species between three canopy stories, based on 40 sample plots. Present species mixture in individual stories is given in parentheses next to the species. Relative abundances of the species after 20 generations are listed in the far right column.
Table 4. Transition probabilities of species between three canopy stories, based on 40 sample plots. Present species mixture in individual stories is given in parentheses next to the species. Relative abundances of the species after 20 generations are listed in the far right column.
Upper story Relative
Understorybeech (0.24)fir (0.48)spruce (0.28)abundances
beech (0.69)0.730.690.680.74
fir (0.15)0.140.110.200.16
spruce (0.05)0.070.090.050.10
maple (0.11)0.060.110.06/*
Middle story
Understorybeech (0.83)fir (0.10)spruce (0.07)
beech0.630.750.680.69
fir0.160.130.160.19
spruce0.100.040.090.12
maple0.110.080.07/*
Upper story
Middle storybeechfirspruce
beech0.750.830.830.77
fir0.150.100.090.14
spruce0.100.070.080.09
* Maple was not present in the middle or upper layer.
Forests 05 01243 g003 1024
Figure 3. Change in tree species composition according to tree density (above) and basal area (below) in the last six decades. Three, two and one stars denote a significance code of p < 0.05, p < 0.01 and p < 0.001, respectively.
Figure 3. Change in tree species composition according to tree density (above) and basal area (below) in the last six decades. Three, two and one stars denote a significance code of p < 0.05, p < 0.01 and p < 0.001, respectively.
Forests 05 01243 g003
In all inventories, tree species composition according to number of individuals in extended dbh size classes in Janj and Lom showed a similar overall pattern (Figure 4). Beech significantly prevailed in the first extended dbh class, while spruce and fir dominated in the third class. The exception was Perucica in 1952 where silver fir was dominant across all dbh classes, while spruce and beech showed a similar pattern according to increasing dbh classes as in the other sites, but with a smaller share.
Forests 05 01243 g004 1024
Figure 4. Changes in tree species composition with respect to density within extended dbh size classes in the Janj old-growth forest at different inventory years (a–h) and comparison with the Lom and Perucica old-growth forests. (a) Share of silver fir, Norway spruce and European beech in Janj old-growth forest in 1952 according to extended dbh classes; (b–d) Same data for Janj old-growth forest in 1965, 2003 and 2011, respectively; (e,f) Same data for Lom old-growth forest in 1952 and 2011, respectively; (g,h) same data for Perucica old-growth forests in 1952 and 2011, respectively.
Figure 4. Changes in tree species composition with respect to density within extended dbh size classes in the Janj old-growth forest at different inventory years (a–h) and comparison with the Lom and Perucica old-growth forests. (a) Share of silver fir, Norway spruce and European beech in Janj old-growth forest in 1952 according to extended dbh classes; (b–d) Same data for Janj old-growth forest in 1965, 2003 and 2011, respectively; (e,f) Same data for Lom old-growth forest in 1952 and 2011, respectively; (g,h) same data for Perucica old-growth forests in 1952 and 2011, respectively.
Forests 05 01243 g004
During the last half of the century, the proportion of beech in the first extended dbh class in Janj continually increased, from 53% in 1952 to 88% in 2011. On the other hand, the share of both spruce and fir in the same size class was in constant decline. In recent decades, beech has become prevalent in the second extended dbh class. Its share also increased in the third extended dbh class, from a modest 5% in 1952 to 12% in 2011. The relative abundance of fir in the second and third extended diameter class was relatively unchanged, while that of spruce moderately decreased (Figure 5a–d). Also, in Lom the share of beech from 1952–2011 increased from 71%–74%, 43%–62% and 9%–12% within the first, second and third extended diameter classes, respectively. During the same period, the share of spruce moderately decreased over all classes, while the share of silver fir decreased only in the mean dbh class (Figure 5e,f). In Perucica, the increase in beech and decrease in conifers were similar to Janj and Lom, but were more pronounced (Figure 5g,h).
Forests 05 01243 g005 1024
Figure 5. Comparison of diameter distributions for the Janj old-growth forest in 1952 (a) and 2011 (b); Lom in 1952 (c) and 2011 (d); and Perucica in 1952 (e) and 2011 (f).
Figure 5. Comparison of diameter distributions for the Janj old-growth forest in 1952 (a) and 2011 (b); Lom in 1952 (c) and 2011 (d); and Perucica in 1952 (e) and 2011 (f).
Forests 05 01243 g005
For cumulative tree dbh distributions (20 cm dbh classes), the best-fitting multiple regression models changed from IQ to RS in Janj, from NE to RS in Lom and from VAR to NE in Perucica (Table 5, Figure 5). In all three old-growth forests, the total number of trees in the first diameter class (7.5–17.5 cm) increased, while it decreased through medium diameters. In 1952, the best-fitting model for both conifers in all three old-growth forests was IQ. This indicates increasing reduction rate in the number of trees from one diameter class to the next with increasing tree size and could be indicative of past disturbance within the group of large diameter trees. In the following 60 years, the fir model changed to UNI in Janj, to RS in Lom and to CO in Perucica. The UNI model may indicate problems with recruitment. The spruce model changed to RS in Lom and remained of the same IQ shape in Janj and Perucica. Changes from IQ to RS and CO may indicate improved conditions for large diameter trees. During the observed period, the shape of beech frequency distributions changed from NE, VAR and CO to CO. This indicates a shift to larger diameter trees.
Table
Table 5. Parameters of the best-fitting model for multiple regressions between the base 10 logarithm of trees per hectare and all possible combinations of dbh.
Table 5. Parameters of the best-fitting model for multiple regressions between the base 10 logarithm of trees per hectare and all possible combinations of dbh.
Janj 1952Janj 2011
FirSpruceBeechCumulativeFirSpruceBeechCumulative
RMSE0.120.110.090.120.120.230.130.06
Adj. R20.920.930.990.950.950.780.990.99
N (dbh class)55456666
ShapeIQIQNEIQUNIIQCORS
Lom 1952Lom 2011
RMSE0.090.090.150.100.080.020.090.04
Adj. R20.970.960.910.980.970.990.990.99
N (dbh class)66566666
ShapeIQIQVARNERSRSCORS
Perucica 1952Perucica 2011
RMSE0.180.090.080.140.030.100.080.09
Adj. R20.940.950.990.970.990.860.990.99
N (dbh class)77777777
ShapeIQIQCOIQCOIQCONE
RMSE—root mean square error; UNI—convex or unimodal shape; RS—rotated sigmoid shape; IQ—increasing-q; CO—concave; VAR—variable (sensu Janowiak, et al. [57]).
Graphical comparison of empirical dbh distributions revealed several similarities among old-growth forests and inventories (Figure 5). In addition to empirical dbh distributions, curves of the best regression models were also inspected (Figure A1). This enabled a direct comparison despite different dbh class thresholds between years of measurements. In all cases, beech dbh distribution follows an approximate straight line if densities are plotted in a logarithmic scale, indicating the closest match with NE. In the recent inventory, beech was more frequent across all dbh classes and also shifted to larger dbh classes. However, the shape remained relatively similar between years and old-growth forests. Beech distribution had the steepest slope of all tree species, which may be attributed to a higher mortality rate.
Distributions of conifers were more irregular in shape and indicated for both species a non-constant q ratio across dbh classes. In all inventories and old-growth forests, the spruce curve was consistently the most “flattened” within smaller dbh classes among all species. Between inventories, spruce density decreased over all classes in Janj and Perucica, but this was more pronounced within small dbh classes. In Lom, spruce density decreased across medium dbh classes and slightly increased in smaller and larger classes (Figure 5c,d). The same was true for fir, while in Perucica fir decreased over all classes.

4. Discussion

4.1. Recent Stand Structure in Janj and Comparison with Similar Old-Growth Forests

Tree species composition according to growing stock in Janj comprised one of the highest shares of conifers (84%) when compared to other old-growth forests in the region [15,44,47,60]. A similar low percentage of beech (18%) was found in the Igman old-growth forest in BiH, which is located about 92 km in the southeast direction. In Janj, spruce was present with an especially high percentage (30%), which can be compared with more continental Carpathian or Alpine old-growth forests, for example Lucka Bela, which also has a 30% share [61]. It is very likely that continentality and high altitude positively influenced the share of spruce in Janj [62]. However, past natural or anthropogenic disturbance events cannot be ruled out. Near the Janj core area is the geographical name “Pozar” in the Serbian language, which translates to fire in English (Figure 1).
In Lom and Perucica conifers were also dominant in BA, with 68% and 65%, respectively [47,50]. With respect to the tree density in all three old-growth forests, beech was the dominant species with 63%, 60% and 50% in Janj, Lom and Perucica, respectively. The average growing stock of 1215 m3∙ha−1 and basal area of 67 m2∙ha−1 in Janj were also among the highest when measured over larger areas and not on a particular selected permanent plot in this type of forest [5,6]. Average growing stock in Lom was lower with 763 m3∙ha−1, and Perucica was between these two values with 937 m3∙ha−1.
The high growing stock in Janj is likely a result of the high percentage of conifers, good site quality and a relatively long period (up to 60 years) without significant natural disturbances. The latter was also indicated by the forest structure, which is absent of large gaps or cohorts of small trees. In Janj, the average CWD volume of 387 m3∙ha−1 was above the average value in Lom (327 m3∙ha−1). Motta [63] observed comparable values of 406 and 420 m3∙ha−1 in mixed Dinaric old-growth forests in Perucica and Biogradska gora, respectively. CWD values from Bosnian mixed old-growth forests were well above the European average of 220 m3∙ha−1 [64] for beech old-growth forests. This may be attributed to the slower decomposition rate of conifers compared to beech [65], their higher overall growing stock and productive sites. The CWD volume to GS ratio in Janj (32%) was lower than in Lom (43%) and also below the European average ratio of 37%, which again may indicate lesser intensity of recent disturbances.
In Janj, the understory and the middle-story were dominated by beech in regard to tree density, while the upper-story was dominated by fir and spruce. Despite beech dominance in lower stories and frequent replacement patterns by beech, there was some indication of reciprocal replacement between species. Moderate reciprocal replacement was often reported for mixed forests with silver fir (e.g., [12,32,66]). Transition matrices indicate the future decrease of beech in the middle-story and increase in the upper-story if the assumptions of the model are to be met [58,59]. The probability that canopy fir and spruce will be replaced by middle story beech does not depend solely on its frequency since both conifers attain a higher position within the canopy. Thus, the assessment of the future species mixture with transition probabilities is above all an indication of potential beech progression.

4.2. Six Decades of Structural Changes in the Janj, Lom and Perucica Old-Growth Forests

Comparison of recent inventories (Motta et al. [47], Lucic [50] and this study) with studies from Drinic [44], Matic et al. [53] and Maunaga et al. [43] may be biased since data were not collected on the same plots. Therefore, emphasis was on changes that were strongly expressed, preferably synchronous across old-growth forests and consistent with the data from recent, more intensive inventories. During the past six decades, the average tree density in Janj and Lom increased by 25% and 10%, respectively, while in Perucica a decrease was indicated. The density increase in Janj and Lom may be attributed to the recent history of endogenous disturbance, but partially also to the 2.5 cm lower dbh limit at the recent inventory. Total basal area increased in Janj by 13%, while it decreased in Lom and Perucica by 6% and 26%, respectively.
The share of spruce and fir in tree density decreased in all reserves, while the share of beech increased. Also, the share of beech in BA increased, while the share of spruce decreased. The share of fir in BA increased in Janj, remained stable in Lom and decreased in Perucica. It seems that fir better compensated for the decrease in density by diameter growth when compared to spruce. Another reason for the more stable fir structure might lie in the lower intensity of indirect anthropogenic disturbance in Bosnia compared to Slovenia, Croatia and Slovakia. Namely, lower SO2 emissions [36] and lower deer densities [35] may be responsible for the less pronounced silver fir decline in Bosnia. Diaci et al. [15] indicated a synchronous silver fir decline in most studied southeast European mixed old-growth mountain forests, while the share of fir in BA increased in the Igman and Trstionica mixed old-growth mountain forests in Bosnia. The data from this study indicated a more pronounced spruce decline in comparison to fir, which may be attributed to its greater susceptibility to drought and high temperatures at the limit of its geographic range towards the Adriatic Sea [62], or its successional stage due to past (anthropogenic) disturbance. This is also in line with the current spruce decline in the northern Velebit National Park in Croatia, but see Battipaglia et al. [67] for different results.
For all inventories and old-growth forests, comparison of tree species composition showed beech dominance in the first extended dbh class, its decrease along extended dbh classes and an increase in conifers. This may be due to the inferior competitive ability of beech in larger dbh classes (see Motta et al. [47]), but since the species mixture changed with time in favour of beech, it may also indicate the onset of beech progression before the 1950s. The lower density of conifers compared with beech may, to some extent, not be in conflict with the demographic equilibrium due to lower mortality of conifers. This was also indicated by the less steeply sloped curves of their dbh distributions (comp. Shimano [68]). However, stable dbh frequency shapes would be expected in time if the different natural mortality rates were the main reason for the dbh distribution shape patterns.
The best fitting curves for dbh frequency of conifers indicated an IQ shape in 1952, while in the recent inventory their shapes were more variable. For conifers, IQ distributions may indicate, due to increasing Q, a higher mortality rate in large-diameter trees and potentially past disturbance, while a unimodal shape (Janj) may indicate heavier past disturbance, a succession or insufficient recruitment [56,57]. NE and RS shapes are often reported as common in old-growth forests [15,56]. In old-growth northern hardwoods, Ducey et al. [69] found deviations from the classic “balanced” exponential distribution on plots heavily stocked with beech; however, in our study beech was sub-dominant to conifers and exhibited a tendency towards NE and CO shapes in all inventories. CO patterns in beech were often observed due to a few widely scattered large trees. It seems that they have filled some of the space made available by the conifer decline.

4.3. Historical Evidence and Possible Causes of Beech Progression

Coexistence of spruce, fir and beech is a complex subject [12,13,17,19,46]. It is driven by the interaction of many factors, anthropogenic and natural. While small scale endogenous disturbance favours fir, intermediate disturbance promotes beech, spruce and admixed broadleaved species [5,19]. Regarding site conditions, conifers are for example more susceptible to competition from ground vegetation than broadleaves and are therefore more successful in the regeneration on rockier sites or acidic soils and on CWD with less abundant ground flora [16]. Besides the natural drivers of stand composition in old-growth forests, the influence of indirect anthropogenic factors may also be strong [3,7,33,70]. Air pollution, over-browsing and position within the agricultural landscape may change natural processes [71].
Factors influencing coexistence act at different spatial scales from tree neighbourhood effects, through the regional impact of diverse deer densities and air pollution to the larger scale impact of climate change. Moreover, species coexistence in mixed mountain forests should also be observed from a long-term perspective since all three species are long-lived and reach ages up to five centuries and more. Many dominant trees in old-growth forests originated two to four centuries ago [42,47,72]. At that time, the densities of ungulates in the Dinaric Mountains were significantly lower, with some species even going extinct due to poaching [14], and the climate was cooler (little ice age, 16th–mid 19th century period; see [73,74,75]). In managed forests, conifers were favoured and beech treated as a non-commercial species. It was often felled as firewood or for charcoal and potash production [9]. This type of forest exploitation did not require sophisticated logging infrastructure and has affected some remote old-growth forests [15]. Before the last century, a cooler climate, unregulated logging of beech, and low densities of ungulates favoured regeneration of conifers.
During the last century, a new constellation of factors has evolved which have given an indirect competitive advantage to beech. Compared to the mid 19th century, warmer climate conditions were already being reported at the beginning of the 20th century [75,76]. Warm spells induced a diffuse decline of conifers in mixed mountain forests on extreme sites and at the limits of their geographic range [4]. After WWII, modern silviculture was introduced in former Yugoslavia, which included beech as a valuable species, while forest grazing was prohibited [77]. In the 1950s, with the onset of accelerated industrialization, air pollution began to be a factor that negatively affected conifers, especially silver fir [21,78]. This species was also weakened by increasing ungulate densities [14]. Higher temperatures and longer vegetation periods coupled with nitrogen deposition may have also contributed to more exuberant ground vegetation [23,79], which represents a greater obstacle for the development of conifer seedlings than those of beech [16].
In BiH, air pollution was less pronounced compared to other parts of Yugoslavia [36], and game densities were never considered to be excessive. Red deer in particular was extremely rare, while roe deer and chamois densities even then rarely exceeded one animal per square kilometre [15,35]. On the other hand, the consequences of climate fluctuations at the beginning of the 20th century in BiH received a great deal of attention in the scientific literature (comp. Safar [4] and references therein). Bark beetle outbreaks on fir and spruce in the Dinaric Mountains were already reported in 1923 and 1924 [38]. A new wave of outbreaks started after the extremely hot and dry summer of 1928, which was followed by a harsh winter (compared to meteorological data in Brunetti et al. [76], Böhm et al. [80] and Kress et al. [81]). This resulted in a patchy and diffuse decline of conifers on south-exposed and shallow soil sites, which was later followed by a large-scale bark beetle outbreak on conifers: spruce, fir and Scots pine (Pinus sylvestris L.), with the exception of black pine (Pinus nigra J.F. Arnold; see [32,37]). Some authors shared the opinion that bark beetle outbreaks were partially caused by unregulated felling and improper handling of felling residues [38,82]. According to Popovic [39], in 1928 only individual conifers at the most exposed sites were weakened by the drought and consequently infested by bark beetles, but at the beginning of 1929, groups of 10–15 trees were already being attacked, and by the end of 1929, this increased to circular patches of 100 trees or more. In 1930, the size of patches reached 1000–5000 trees. Healthy trees were also attacked, but not small trees and saplings. At the start of the outbreak, its upper borderline was at about 1200 m, but by 1930 it reached an altitude of 1400 m a.s.l. Tregubov [39], an early researcher of old-growth forests in the Lom region, also mentioned these events.
The next important period of high summer temperatures and drought was in the late 1940s and early 1950s (comp. to Brunetti et al. [76], Böhm et al. [80] and Kress et al. [81]). For these events, reports of bark beetle outbreaks in the Perucica old-growth forest are available. According to Eic [41] and Golubovic [33], in the fir dominated part of Perucica, 20,000 m3 of infested conifer logs were cut down, which represented about 10% of the exploitable timber. Logs were not properly handled, increasing the potential for beetle infestations. Eic [41] further stated that the increased infection was also due to fires in neighbouring forests after heavy forest exploitation.
The combination of the bark beetle outbreaks in the early 1920s, the severe climatic event in 1928 followed by the droughts of the late 1940s and early 1950s, and anthropogenic disturbances [4] had weakened the competitive strength of conifers and probably induced a large-scale synchronous regeneration wave of beech (Figure 6), which had already been observed by the early researchers of old-growth forests [4,46].
Forests 05 01243 g006 1024
Figure 6. A dense, ubiquitous small tree layer of beech with scattered conifer seedling clusters below is typical for a large part of the Janj reserve (November 1, 2012; photo: J. Diaci).
Figure 6. A dense, ubiquitous small tree layer of beech with scattered conifer seedling clusters below is typical for a large part of the Janj reserve (November 1, 2012; photo: J. Diaci).
Forests 05 01243 g006
Beech seed crops are not as frequent as those of conifers [83]; however, beech mast years have become more frequent in the last century, possibly due to temperature increases, climatic extremes and increased atmospheric nitrogen [22]. Regarding disturbance events in BiH from the 1920s to the early 1950s, it is also important that early summer drought is a strong masting predictor for beech [84]. Large-scale synchronous regeneration of beech has probably delayed the regeneration of conifers, lasting for several decades [19]. Documented climatic events, likely coupled with anthropogenic disturbances, are indicated in this study by the following: (1) beech dominance among small trees in all inventories and its increase overall dbh classes in recent inventories; (2) the future decrease of the beech middle-story modelled with transition probabilities; (3) the IQ shape of dbh frequency distribution in the 1950s for conifers in all old-growth forests and (4) changes in dbh frequency distributions. Other recent studies have also suggested these events, e.g., the disturbance study in Motta et al. [47] and gap age analysis in Bottero et al. [48].

5. Conclusions

The recent forest structure of the Janj old-growth forest in comparison to other Dinaric and European mixed mountain forests is characterized by a high growing stock and a large proportion of conifers, particularly spruce. The former may be due to site quality, dominant conifers, and the absence of large disturbances in recent decades. The high share of conifers in comparison to other Dinaric mixed old-growth forests in Slovenia and Croatia could be attributed to the higher altitude and continental climate. However, other Bosnian mixed old-growth forests (Igman, Lom and partly Perucica) are in a similar altitudinal belt and exhibit fewer conifers; therefore, past natural or anthropogenic disturbances may have also influenced the high share of conifers. This seems likely since spruce dbh distribution changed in time and the recent distribution does not show sustainable recruitment.
With respect to tree height and diameter size, beech is slightly inferior to both coniferous trees, but many changes in structural parameters indicate its progression (e.g., increasing abundance within small trees; changes in dbh distributions, GS and BA; neighbourhood effects). The causes of conifer regression are likely associated with climatic extremes, especially high summer temperatures and droughts, which have weakened their vitality and influenced bark beetle infestations. However, many anthropogenic disturbances may have accelerated this process (unregulated and sanitary fellings in old-growth forests and in the near vicinity, burning of slash and clearing of land for fields and pasturing, building of logging infrastructure). Some evidence suggests that the sequence of droughts, bark beetle infestations, and anthropogenic disturbances during the 1920s, 1940s and early 1950s had diffusely opened the canopy layer of the three analysed old-growth forests in BiH, reduced conifer and increased beech seed potential, and thus triggered a wave of beech regeneration. Conifers did regenerate below the dense beech seedling and sapling layer, but with a certain delay, thus their density and competitive ability probably decreased with each wave of drought and proliferous beech regeneration.
The results of this study indicate conifer decline in some Bosnian old-growth forests. It seems that the decline was not continuous but rather followed extreme climatic events (summer droughts and harsh winters), most likely in conjunction with bark beetle outbreaks and anthropogenic disturbances. This hypothesis deserves in-depth research in the future. If confirmed, it could contribute to a better understanding of future changes in forest composition due to climate change.

Acknowledgments

The authors gratefully acknowledge the constructive comments and statistical assistance provided by Gary Kerr. We would also like to thank Thomas A. Nagel for his advice on calculating transition probabilities. This research was funded by the Slovenian Research Agency, Research Program P4-0059 and Ministry of Science and Technology, the Republic of Srpska, Bosnia-Herzegovina.

Author Contributions

Initial ideas and financing of the research: Jurij Diaci, Renzo Motta and Srdjan Keren; data compilation and evaluation: all authors; manuscript drafting: Srdjan Keren and Jurij Diaci; all authors contributed to revising successive drafts.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix

Table
Table A1. Heights of thickest trees of major tree species in old-growth Janj.
Table A1. Heights of thickest trees of major tree species in old-growth Janj.
SpeciesEuropean beechSilver firNorway spruce
Nodbhheightdbhheightdbhheight
14831.27238.66837.9
24830.47340.37036.3
35131.57434.77035.5
45330.57438.47137.5
55440.87541.07437.4
65633.48042.37540.1
75732.18145.07541.8
85830.08542.07543.0
95833.08635.37736.3
105839.09439.47734.6
115936.09649.57839.3
125937.29839.08346.2
136032.59944.58543.5
146039.210850.59141.9
156935.611048.59244.9
167035.611046.29544.5
177235.511747.59643.4
187935.411744.89944.3
198344.712046.510045.5
209241.612444.811045.0
Mean62.235.394.742.983.141.0
Forests 05 01243 g007 1024
Figure A1. Temporal changes in species-specific diameter distributions of old-growth forests (a) Janj, (b) Lom, (c) Perucica. Curves and lines represent the best regression models. The first letter indicates old-growth forest (e.g., l denotes Lom), the next two letters represent Latin species names (e.g., aa denotes Abies alba) and the following two numbers indicate the year of measurement (e.g., 11 denotes 2011).
Figure A1. Temporal changes in species-specific diameter distributions of old-growth forests (a) Janj, (b) Lom, (c) Perucica. Curves and lines represent the best regression models. The first letter indicates old-growth forest (e.g., l denotes Lom), the next two letters represent Latin species names (e.g., aa denotes Abies alba) and the following two numbers indicate the year of measurement (e.g., 11 denotes 2011).
Forests 05 01243 g007

References

  1. Bohn, U. General Map of Natural Vegetation of Europe; German Federal Agency for Natural Conservation: Bonn, Germany, 1996. [Google Scholar]
  2. Ellenberg, H. Vegetation Ecology of Central Europe; Cambridge University Press: Cambridge, UK, 1988. [Google Scholar]
  3. Peterken, G.F. Natural Woodland: Ecology and Conservation in Northern Temperate Regions; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  4. Safar, J. Ugibanje i Obnavljanje Jele u Prebornim Šumama Gorskog Kotara. Sumarski List 1951, 75, 299–303. [Google Scholar]
  5. Leibundgut, H. Europäische Urwälder der Bergstufe; Haupt Verlag AG: Bern, Switzerland, 1982. [Google Scholar]
  6. Korpel, S. Die Urwälder der Westkarpaten; Gustav Fischer Verlag: Stuttgart, Germany, 1995. [Google Scholar]
  7. Vrska, T.; Adam, D.; Hort, L.; Kolár, 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]
  8. Malek, J. Jedle v Pralesovitých Rezervacích a Vybernych Lesich. Lesnicka Prace 1981, 60, 255–259. [Google Scholar]
  9. Matic, S. Utjecaj Ekoloskih i Strukturnih Cinilaca na Prirodno Pomladivanje Prebornih Suma Jele i Bukve u Gorskom Kotaru. Glasnik za Sumske Pokuse 1983, 21, 223–400. [Google Scholar]
  10. Ficko, A.; Poljanec, A.; Boncina, A. Do changes in Spatial Distribution, Structure and Abundance of Silver Fir (Abies alba Mill.) Indicate Its Decline? For. Ecol. Manag. 2011, 261, 844–854. [Google Scholar] [CrossRef]
  11. Boncina, A.; Cavlovic, J.; Curovic, M.; Govedar, Z.; Klopcic, M.; Medarevic, M. A Comparative Analysis of Recent Changes in Dinaric Uneven–Aged Forests of the NW Balkans. Forestry 2014, 87, 71–84. [Google Scholar] [CrossRef]
  12. Safar, J. Problemi Izmjene Vrsta u Šumama. Sumarski List 1952, 76, 89–100. [Google Scholar]
  13. Paluch, J. The Spatial Pattern of a Natural European Beech (Fagus sylvatica L.)–Silver Fir (Abies alba Mill.) Forest: A Patch-Mosaic Perspective. For. Ecol. Manag. 2007, 253, 161–170. [Google Scholar] [CrossRef]
  14. Klopcic, M.; Jerina, K.; Boncina, A. Long-Term Changes of Structure and Tree Species Composition in Dinaric Uneven-Aged Forests: Are Red Deer an Important Factor? Eur. J. For. Res. 2010, 129, 277–288. [Google Scholar] [CrossRef]
  15. 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, 5, 479–491. [Google Scholar]
  16. Ott, E.; Frehner, M.; Frey, H.-U.; Lüscher, P. Gebirgsnadelwälder: Praxisorientierter Leitfaden Für Eine Standortgerechte Waldbehandlung; Verlag Paul Haupt: Bern, Switzerland, 1997. [Google Scholar]
  17. Mlinsek, D. Pomlajevanje in Nekatere Razvojne Znacilnosti Bukovega in Jelovega Mladovja v Pragozdu na Rogu. Zbornik Biotehniske Fakultete 1967, 15, 7–32. [Google Scholar]
  18. Wraber, M. Topografski, Ekoloski in Socioloski Podatki o Slovenskih Pragozdovih. In Juznoevropske Prasume i Visokoplaninska Flora i Vegetacija Istocnoalpsko-Dinarskog Prostora; Fukarek, P., Ed.; Academy of Sciences Bosnia and Herzegovina: Sarajevo, Bosnia and Herzegovina, 1970; pp. 91–109. [Google Scholar]
  19. Nagel, T.A.; Svoboda, M.; Rugani, T.; Diaci, J. Gap Regeneration and Replacement Patterns in an Old-Growth Fagus-Abies Forest of Bosnia and Herzegovina. Plant Ecol. 2010, 208, 307–318. [Google Scholar] [CrossRef]
  20. Motta, R. Impact of Wild Ungulates on Forest Regeneration and Tree Composition of Mountain Forests in the Western Italian Alps. For. Ecol. Manag. 1996, 88, 93–98. [Google Scholar]
  21. Elling, W.; Dittmar, C.; Pfaffelmoser, K.; Rötzer, T. Dendroecological Assessment of the Complex Causes of Decline and Recovery of the Growth of Silver Fir (Abies alba Mill.) in Southern Germany. For. Ecol. Manag. 2009, 257, 1175–1187. [Google Scholar] [CrossRef]
  22. Drobyshev, I.; Overgaard, R.; Saygin, I.; Niklasson, M.; Hickler, T.; Karlsson, M.; Sykes, M.T. Masting Behaviour and Dendrochronology of European Beech (Fagus sylvatica L.) in Southern Sweden. For. Ecol. Manag. 2010, 259, 2160–2171. [Google Scholar] [CrossRef]
  23. Kellner, O.; Redbo-Torstensson, P. Effects of Elevated Nitrogen Deposition on the Field-Layer Vegetation in Coniferous Forests. Ecol. Bull. 1995, 227–237. [Google Scholar]
  24. Frelich, L.E.; Reich, P.B. Neighborhood Effects, Disturbance Severity and Community Stability in Forests. Ecosystems 1999, 2, 151–166. [Google Scholar] [CrossRef]
  25. Canham, C.D.; Denslow, J.S.; Platt, W.J.; Runkle, J.R.; Spies, T.A.; White, P.S. Light Regimes Beneath Closed Canopies and Tree-Fall Gaps in Temperate and Tropical Forests. Can. J. For. Res. 1990, 20, 620–631. [Google Scholar] [CrossRef]
  26. Zinke, P.J. The pattern of Influence of Individual Forest Trees on Soil Properties. Ecology 1962, 43, 130–133. [Google Scholar] [CrossRef]
  27. Mladenoff, D.J. Dynamics of Nitrogen Mineralization and Nitrification in Hemlock and Hardwood Treefall Gaps. Ecology 1987, 68, 1171–1180. [Google Scholar] [CrossRef]
  28. Finzi, A.C.; Canham, C.D.; van Breemen, N. Canopy Tree Soil Interactions within Temperate Forests: Species Effects on pH and Cations. Ecol. Appl. 1998, 8, 447–454. [Google Scholar]
  29. Packer, A.; Clay, K. Soil Pathogens and Spatial Patterns of Seedling Mortality in a Temperate Tree. Nature 2000, 404, 278–281. [Google Scholar] [CrossRef]
  30. Fox, J.F. Alternation and Coexistence of Tree Species. Am. Nat. 1977, 111, 69–89. [Google Scholar]
  31. Woods, K.D. Patterns of Tree Replacement—Canopy Effects on Understory Pattern in Hemlock—Northern Hardwood Forests. Vegetatio 1984, 56, 87–107. [Google Scholar]
  32. Tregubov, V. Les Forets Vierges Montagnardes des Alps Dinariques. Massif de Klekovatcha—Guermetch. Ph.D. Thesis, Montpellier University, Montpellier, France, 1941. [Google Scholar]
  33. Golubovic, U. Da li je Peručica Prašuma? Šumarski List 1963, 11/12, 484–486. [Google Scholar]
  34. Uotila, A.; Kouki, J.; Kontkanen, H.; Pulkkinen, P. Assessing the Naturalness of Boreal Forests in Eastern Fennoscandia. For. Ecol. Manag. 2002, 161, 257–277. [Google Scholar] [CrossRef]
  35. Anonymous. Census of Game Density in Hunting Area Sipovo, Sume Republike Srpske. Unpublished Report. Ministry of Agriculture, Forestry and Water Resources of the Republika Srpska: Banja Luka, Bosnia and Herzegovina, 2013. [Google Scholar]
  36. Berge, E.; Bartnicki, J.; Olendrzynski, K.; Tsyro, S.G. Long-Term Trends in Emissions and Transboundary Transport of Acidifying Air Pollution in Europe. J. Environ. Manag. 1999, 57, 31–50. [Google Scholar] [CrossRef]
  37. Bambulovic, P. Susenje Cetinjavih Suma u Bosni. Sumarski List 1930, 54, 446–450. [Google Scholar]
  38. Batic, J. Susenje Jelovih i Smrekovih Suma u Lici. Sumarski List 1930, 54, 345–347. [Google Scholar]
  39. Popovic, J. Susenje Cetinara Zbog Napada Potkornjaka u Šumama Drinske i Vrbaske Banovine; Radovi fitopatoloskog zavoda u Sarajevu: Sarajevo, Bosnia and Herzegovina, 1931. [Google Scholar]
  40. Frölich, J. Urwaldpraxis; Neumann Verlag: Radebeul/Berlin, Germany, 1954. [Google Scholar]
  41. Eic, N. Nase presume. Narodni Sumar 1951, 2/3, 56–67. [Google Scholar]
  42. Schütz, J.-P. Etude des Phénomenes de la Croissance en Hauteur et en Diametre du Sapin (Abies alba Mill.) et de l’épicéa (Picea abies Karst.) Dans Deux Peuplements Jardinés et une Foret Vierge. Ph.D. Thesis, ETH, Zurich, Switzerland, 1969. [Google Scholar]
  43. Maunaga, Z.; Burlica, C.; Pavlovic, B.; Bozalo, G.; Koprivica, M.; Lazarev, V.; Mataruga, M.; Rapajic, Z.; Govedar, Z.; Stanivukovic, Z.; et al. Studija Prasume Janj i Lom; Sumarski Fakultet Univerziteta u Banjaluci: Banjaluka, Bosnia and Herzegovina, 2005. [Google Scholar]
  44. Drinic, P. Taksacioni Elementi Sastojina Jele, Smrce i Bukve Prasumskog Tipa u Bosni. Radovi Poljoprivredno-Sumarskog Fakulteta u Sarajevu 1956, 1, 107–160. [Google Scholar]
  45. Fukarek, P. Sumske Zajednice Prasumskog Rezervata Perucice u Bosni. In Juznoevropske Prasume i Visokoplaninska Flora i Vegetacija Istocnoalpsko-Dinarskog Prostora; Fukarek, P., Ed.; Academy of Sciences Bosnia and Herzegovina: Sarajevo, Bosnia and Herzegovina, 1970; pp. 157–262. [Google Scholar]
  46. Pintaric, K. Urwald Perucica als Naturliches Forschungslaboratorium. Algebr. For. Z. 1978, 33, 691–693. [Google Scholar]
  47. Motta, R.; Berretti, R.; Castagneri, D.; Dukic, 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 Forests: The Nearly Steady-State Forest of Lom (Bosnia-Herzegovina). Can. J. For. Res. 2011, 41, 1871–1884. [Google Scholar] [CrossRef]
  48. Bottero, A.; Garbarino, M.; Dukic, V.; Govedar, Z.; Lingua, E.; Nagel, T.A.; Motta, R. Gap-Phase Dynamics in the Old-Growth Forest of Lom, Bosnia and Herzegovina. Silva Fenn. 2011, 45, 875–887. [Google Scholar]
  49. Govedar, Z. Nacini Prirodnog Obnavljanja Mjesovitih Suma Jele i Smrce (Abieti-Piceetum illyricum) na Podrucju Zapadnog Djela Republike Srpske. Ph.D. thesis, Belgrade University, Belgrade, Serbia, 2005. [Google Scholar]
  50. Lucic, R. Structural and Productivity Characteristics of Forests with Old-Growth Character in the National Park Sutjeska. Master Thesis, Belgrade University, Belgrade, Serbia, 2012. [Google Scholar]
  51. Van Wagner, C.E. The Line Intersect Method in Forest Fuel Sampling. For. Sci. 1968, 14, 20–26. [Google Scholar]
  52. Goedhart, P.W.; van Hees, A.F.M.; van der Werf, B.; Kramer, K. Biodiversity on Coarse Woody Debris, Model Description and Parameterisation. In EU 5.OP Project Nat Man Working Report; Wageningen UR: Wageningen, The Netherlands, 2004. [Google Scholar]
  53. Matic, V.; Drinic, P.; Stefanovic, V.; Ciric, M. Stanje Šuma u SR Bosni i Hercegovini Prema Inventuri Šuma na Velikim Površinama u 1964–1968. Godini, Šumarski Fakultet i Institut za Sumarstvo u Sarajevu 1971, 7, 1–639. [Google Scholar]
  54. Matic, V. Prirast i Prinos Suma; Faculty of Forestry: Sarajevo, Bosnia and Herzegovina, 1980. [Google Scholar]
  55. Motta, R.; Berretti, R.; Lingua, E.; Piussi, P. Coarse Woody Debris, Forest Structure and Regeneration in the Valbona Forest Reserve, Paneveggio, Italian Alps. For. Ecol. Manag. 2006, 235, 155–163. [Google Scholar] [CrossRef]
  56. Leak, W.A. Long-Term Structural Change in Uneven-Aged Northern Hardwoods. For. Sci. 1996, 42, 160–165. [Google Scholar]
  57. 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]
  58. White, P.S.; Mackenzie, M.D.; Busing, R.T. A Critique on Overstory Understory Comparisons Based on Transition-Probability Analysis of an Old Growth Spruce-Fir Stand in the Appalachians. Vegetatio 1985, 64, 37–45. [Google Scholar] [CrossRef]
  59. Stevens, M.H.H. A Primer of Ecology with R; Springer: Dordrecht, The Netherlands; New York, NY, USA, 2009. [Google Scholar]
  60. Anic, I.; Mikac, S. Struktura, Tekstura i Pomlađivanje Dinarske Bukovo-Jelove Prasume Čorkova uvala. Sumarski List 2008, 11–12, 505–515. [Google Scholar]
  61. Firm, D.; Nagel, T.A.; Diaci, J. Disturbance History and Dynamics of an Old-Growth Mixed Species Mountain Forest in the Slovenian Alps. For. Ecol. Manag. 2008, 257, 1893–1901. [Google Scholar] [CrossRef]
  62. Horvat, I.; Glavac, V.; Ellenberg, H. Vegetation Sudosteuropas; G. Fischer: Stuttgart, Germany, 1974. [Google Scholar]
  63. 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. 2014. submitted. [Google Scholar]
  64. Christensen, M.; Hahn, K.; Mountford, P.E.; Odor, P.; Standovar, T.; Rozenbergar, D.; Diaci, J.; Wijdeven, S.; Meyer, P.; Winter, S.; et al. Dead Wood in European Beech (Fagus sylvatica) Forest Reserves. For. Ecol. Manag. 2005, 210, 267–282. [Google Scholar] [CrossRef]
  65. Harmon, M.E.; Franklin, J.F.; Swanson, F.J.; Sollins, P.; Gregory, S.V.; Lattin, J.D.; Anderson, N.H.; Cline, S.P.; Aumen, N.G.; Sedell, J.R.; et al. Ecology of Coarse Woody Debris in Temperate Ecosystems. Adv. Ecol. Res. 1986, 15, 133–302. [Google Scholar] [CrossRef]
  66. Simak, M. Untersuchungen Über den Natürlichen Baumartenwechsel in Schweizerischen Plenterwäldern. Mitteilungen Schweiz. Anst. Forstl. Vers. 1951, 27, 406–468. [Google Scholar]
  67. Battipaglia, G.; Saurer, M.; Cherubini, P.; Siegwolf, R.T.W.; Cotrufo, M.F. Tree Rings Indicate Different Drought Resistance of Native (Abies alba Mill.) and a Nonnative (Picea abies (L.) Karst.) Species Co-occurring at a Dry Site in Southern Italy. For. Ecol. Manag. 2009, 257, 820–828. [Google Scholar] [CrossRef]
  68. Shimano, K. A Power Function for Forest Structure and Regeneration Pattern of Pioneer and Climax Species in Patch Mosaic Forests. Plant Ecol. 2000, 146, 207–220. [Google Scholar] [CrossRef]
  69. Ducey, J.M.; Gunn, S.J.; Whitman, A.A. Late-Successional and Old-Growth Forests in the Northeastern United States: Structure, Dynamics, and Prospects for Restoration. Forests 2013, 4, 1055–1086. [Google Scholar] [CrossRef]
  70. Diaci, J.; Rozenbergar, D.; Boncina, A. Stand Dynamics of Dinaric Old-Growth Forest in Slovenia: Are Indirect Human Influences Relevant? Plant Biosyst. 2010, 144, 194–201. [Google Scholar] [CrossRef]
  71. Nowacki, G.J.; Abrams, M.D. Forest Composition, Structure, and Disturbance History of the Alan-Seeger-Natural-Area, Huntington County, Pennsylvania. Bull. Torrey Bot. Club 1994, 121, 277–291. [Google Scholar]
  72. Mlinsek, D. Rast in Sposobnost Reagiranja Pragozdne Bukve. Zbornik Biotehniške Fakultete 1967, 15, 63–79. [Google Scholar]
  73. Sercelj, A. The Origins and Development of Forests in Slovenia; Academia Scientiarum et Artium Slovenica, Classis IV: Historia Naturalis, Ljubljana, Slovenia, 1996. [Google Scholar]
  74. Mann, M.E.; Bradley, R.S.; Hughes, M.K. Northern Hemisphere Temperatures During the Past Millennium: Inferences, Uncertainties, and Limitations. Geophys. Res. Lett. 1999, 26, 759–762. [Google Scholar]
  75. Buntgen, U.; Kyncl, T.; Ginzler, C.; Jacks, D.S.; Esper, J.; Tegel, W.; Heussner, K.U.; Kyncl, J. Filling the Eastern European Gap in Millennium-Long Temperature Reconstructions. Proc. Natl. Acad. Sci. USA 2013, 110, 1773–1778. [Google Scholar]
  76. Brunetti, M.; Maugeri, M.; Monti, F.; Nanni, T. Temperature and Precipitation Variability in Italy in the Last Two Centuries from Homogenised Instrumental Time Series. Int. J. Climatol. 2006, 26, 345–381. [Google Scholar] [CrossRef]
  77. Mlinsek, D. From Clear-Cutting to a Close-to-Nature Silvicultural System. IUFRO News 1996, 25, 6–8. [Google Scholar]
  78. Prpic, B. Propadanje Suma u SR Hrvatskoj i Jugoslaviji. Sumarksi List 1989, 113, 235–243. [Google Scholar]
  79. Sebesta, J.; Samonil, P.; Lacina, J.; Oulehle, F.; Houska, J.; Bucek, A. Acidification of Primeval Forests in the Ukraine Carpathians: Vegetation and Soil Changes over Six Decades. For. Ecol. Manag. 2011, 262, 1265–1279. [Google Scholar] [CrossRef]
  80. Böhm, R.; Auer, I.; Schöner, W.; Ganekind, M.; Gruber, C.; Jurkovic, A.; Orlik, A.; Ungersböck, M. Eine neue Webseite mit Instrumentellen Qualitäts-Klimadaten für den Grossraum Alpen zurück bis 1760. Wiener Mitteilungen 2009, 216, 7–20. [Google Scholar]
  81. Kress, A.; Saurer, M.; Siegwolf, R.T.; Frank, D.C.; Esper, J.; Bugmann, H. A 350 Year Drought Reconstruction from Alpine Tree Ring Stable Isotopes. Glob. Biogeochem. Cycles 2010, 24, 1–16. [Google Scholar]
  82. Omanovic, S. Izvjestaji—Pregled Radova na Suzbijanju Potkornjaka na Podrucju Direkcije sSuma u Sarajevu. Sumarski List 1932, 6, 374–380. [Google Scholar]
  83. Kantorowicz, W. Half a Century of Seed Years in Major Tree Species of Poland. Silvae Genetica 2000, 49, 245–249. [Google Scholar]
  84. Piovesan, G.; Adams, J.M. Masting Behaviour in Beech: Linking Reproduction and Climatic Variation. Can. J. Bot.-Rev. Can. Bot. 2001, 79, 1039–1047. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

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. https://0-doi-org.brum.beds.ac.uk/10.3390/f5061243

AMA Style

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(6):1243-1266. https://0-doi-org.brum.beds.ac.uk/10.3390/f5061243

Chicago/Turabian Style

Keren, Srdjan, Renzo Motta, Zoran Govedar, Radovan Lucic, Milan Medarevic, and Jurij Diaci. 2014. "Comparative Structural Dynamics of the Janj Mixed Old-Growth Mountain Forest in Bosnia and Herzegovina: Are Conifers in a Long-Term Decline?" Forests 5, no. 6: 1243-1266. https://0-doi-org.brum.beds.ac.uk/10.3390/f5061243

Article Metrics

Back to TopTop