Next Article in Journal
The Terebelliformia-Recent Developments and Future Directions
Next Article in Special Issue
Flora and Vegetation Characteristics of the Natural Habitat of the Endangered Plant Pterygopleurum neurophyllum
Previous Article in Journal
Western Pacific Zooplankton Community along Latitudinal and Equatorial Transects in Autumn 2017 (Northern Hemisphere)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity of Ostrya carpinifolia Forests in Ravine Habitats of Serbia (S-E Europe)

by
Dimitrije Sekulić
1,*,
Branko Karadžić
1,
Nevena Kuzmanović
2,
Snežana Jarić
1,
Miroslava Mitrović
1 and
Pavle Pavlović
1
1
Department of Ecology, Institute for Biological Research “Siniša Stanković”—National Institute of Republic of Serbia, University of Belgrade, Bulevar despota Stefana 142, 11000 Belgrade, Serbia
2
Institute of Botany and Botanical Garden, Faculty of Biology, University of Belgrade, Takovska 43, 11000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Diversity 2021, 13(2), 59; https://doi.org/10.3390/d13020059
Submission received: 31 December 2020 / Revised: 27 January 2021 / Accepted: 29 January 2021 / Published: 3 February 2021

Abstract

:
We investigated vegetation in ravine habitats of Serbia, in order to classify hop hornbeam (Ostrya carpinifolia Scop.) forests in syntaxonomic terms, assess the effects of environmental factors on their floristic differentiation, and detect the biodiversity components of the analyzed communities. Both K-means clustering and Bayesian classification revealed five ecologically interpretable groups of forests that belong to the alliances Ostryo carpinifoliae-Fagion sylvaticae, Ostryo carpinifoliae-Tilion platyphylli, Fraxino orni-Ostryion carpinifoliae, Pseudofumario albae-Ostryion carpinifoliae, and Achilleo ageratifoliae-Ostryion carpinifoliae. Canonical correspondence analysis indicated that these alliances are clearly differentiated along a combined light–moisture gradient (from shade and mesic to sunny and xeric variants). The alpha diversity increases from xeric to mesic alliances. A lower alpha diversity in xeric forests may be explained by the stress conditions that prevent mesic species from colonizing the saxatile habitats. Extremely high—almost the greatest possible—values of both the species turnover and beta diversity were detected in all variants of the analyzed forests. Such high diversity may be the result of the strong environmental gradients in ravine habitats. The investigated forests represent an important pool of rare, paleo-endemic species that survived Quaternary glaciations in ravine refugia.

1. Introduction

Ostrya carpinifolia Scop., commonly known as hop hornbeam, dominates in xeric amphi-Adriatic alliances Carpinion orientalis Horvat 1954 and Fraxino orni-Ostryion Tomažić 1940 of the class Quercetea pubescentis Doing-Kraft ex Scamoni et Passarge 1959 [1,2,3,4]. Čarni et al. [4] and Stupar et al. [5] analyzed the ecological, floristic, and chorological differences between these alliances. In addition to the zonal belt of sub-Mediterranean xeric forests, the hop hornbeam forms mesic communities that are extrazonally distributed in inland parts of the Apennines and the Balkan Peninsula [2,3,4,5,6,7,8].
Ostrya carpinifolia Scop. is a species that is adapted to a warm to moderate climate, but it avoids a strictly Mediterranean climate with a marked summer drought. This species has wide tolerance limits with respect to both light and moisture gradients [9]. Due to its high tolerance to adverse soil conditions, it is most commonly located in canyons [2,3,8]. The stress conditions found in skeletal soils, and on screes and bare rocks, prevent stronger competitors from replacing the hop hornbeam in ravine habitats.
Ravine refugia of the Balkan Peninsula were important for species survival during periodic glaciations and for the postglacial recolonization of northern Europe during interglacial stages of the Quaternary [10,11,12,13,14].
Despite numerous studies having been conducted on ravine vegetation in central parts of the Balkan Peninsula [3,4,5,6,7,8,15,16,17,18,19,20,21,22,23], available information on the biodiversity components of Ostrya carpinifolia forests is incomplete. Therefore, we performed investigations of ravine vegetation in Serbia, in order to (1) classify hop hornbeam forests, (2) assess the effects of environmental factors on their floristic differentiation, and (3) detect patterns of biodiversity components (alpha and beta diversity) in investigated communities.

2. Materials and Methods

2.1. Study Area and Sampling Procedures

Using the Braun-Blanquet [24] sampling procedure, we collected 144 relevés from canyons of the Gradac, Trešnjica, Derventa, Djetinja, Ljutina, Mileševka, and Lim rivers (Figure 1). The set of originally collected relevés is stored in the FLORA database, and it is available upon request. The geographic coordinates of the study area are presented in Table 1. The size of the sampling sites varied from 15 × 15 m in mesic forests to 20 × 20 m in xeric communities. The importance of a species at a site was assessed using the numeric scale of combined cover-abundance values [25]. The taxonomic nomenclature was harmonized with the nomenclature for vascular plants and bryophytes in the Plant List Database (http://www.theplantlist.org).

2.2. Statistical Analyses

The classification of relevés may be performed using a wide spectrum of clustering methods. Hierarchical classification methods (different variants of agglomerative and divisive clustering) may produce a high number of misclassifications. To avoid the risk of misclassifications, we used both K-means clustering [26] and Bayesian classification [27]. These non-hierarchical methods enable the allocation of misclassified relevés and the formation of clusters that are as homogeneous as possible. Subjectivity in the initial estimation of the number of clusters is the main drawback of both methods. Marinković et al. [28] proposed a simple procedure to avoid this problem. It selects the number of clusters by maximizing the variance ratio:
V R = σ B 2 σ W 2 ,
where σ2B denotes the between-group variance (i.e., variance of cluster centroids), and σ2W is the within-group variance (the sum of variances within each of the k clusters). The fidelity of each species to different clusters was assessed using the phi coefficient [29].
The effects of environmental variables on the floristic differentiation of communities were assessed using canonical correspondence analysis [30]. Environmental variables included topographic parameters (altitude, slope, and aspect) and a set of variables (light, moisture, soil acidity, soil nutrients, temperature, and continentality) that were estimated indirectly, using the weighted average of Ellenberg’s indicator values [31]. The stepwise forward selection procedure [27] was used to detect environmental variables with statistically significant effects on vegetation. During each step of the procedure, we expanded the multiple regression model by adding an environmental variable explaining most of the residual variance (i.e., the variance of vegetation data, not explained by previously selected environmental variables). The statistical significance of the hypothesis that vegetation is independent of selected environmental variables was assessed using the non-parametric Monte Carlo permutation test (3000 permutations, p < 0.05).
The alpha diversity was assessed using Shannon’s entropy:
H = i = 1 s p i log p i ,
where pi is the proportion of species i within a site, and s denotes the number of species within the site. Pielou [32] proved that the equitability component of alpha diversity can be calculated using the following equation:
E = H H m a x = i = 1 s p i log p i / log s ,
where Hmax is the greatest possible entropy.
The beta diversity represents the variability in species composition among communities [33,34]. Floristic difference between a pair of communities is attributable to both species turnover and species richness [35,36,37,38,39,40,41]. Two communities maximally differ if they have no species in common (Figure 2a). Such a situation corresponds to complete species replacement from one to another community. However, a pair of communities may differ not because of species replacement, but due to a difference in species richness. In that case, the species-richer community may be considered a nest for the species-poorer community (Figure 2c). Therefore, the beta diversity component related to differences in species richness is frequently denoted as “nestedness”. Floristic differences between two communities are usually caused by both species turnover and nestedness (Figure 2b).
Components of the beta diversity were determined using the procedures described by Baselga [35] and Podani et al. [37]. Statistical analyses were performed using the updated version of the FLORA software [42].

3. Results

We recorded 353 species of vascular plants in the investigated forests. A high percentage of phanerophytes (trees and shrubs) was characteristic for the analyzed communities. The screes and cliffs are covered by sparsely distributed individuals of Ostrya carpinifolia Scop. and other xeric trees. On more favorable soils, the hop hornbeam occurs in combination with numerous xeric and mesic species. Compared to the phanerophytes, the herbaceous plants were more diverse. A high percentage of herbaceous species belonged to the group of endemic Tertiary relicts. A selection of the huge amount of paleo-endemic species occurring in the investigated forests is shown in Figure 3.

3.1. Classification of Communities

The results of K-means clustering and Bayesian classification showed that the analyzed forests can be grouped into five homogeneous clusters (Table 2).
For each cluster, a set of diagnostic species was detected using the phi coefficient. A synoptic table of the analyzed forests is presented in Appendix A.
Cluster 1 comprises mesophilous, shady forests on deep soil. Diagnostic species of the cluster are mesic species (Asperula taurina L., Epimedium alpinum L., Fraxinus excelsior L., Evonymus europaeus L., Alliaria petiolata (M.Bieb.) Cavara & Grande, Pulmonaria officinalis L., Melampyrum hoermannianum K. Malý, Angelica sylvestris L., Parietaria officinalis L., Athyrium filix-femina (L.) Roth, Cardamine bulbifera (L.) Crantz, Cyclamen purpurascens Mill., Acer pseudoplatanus L., and Brachypodium sylvaticum (Huds.) P.Beauv.). The forests included in this cluster may be assigned to the mesic alliance Ostryo carpinifoliae-Fagion sylvaticae Borhidi 1963 of the order Fagetalia sylvaticae Pawłowski 1928.
Cluster 2 consists of mesic forests on both deep and skeletal soils. Its diagnostic species included mesic species (Salvia glutinosa L., Hedera helix L., Asarum europaeum L., Tilia platyphyllos Scop., Scrophularia nodosa L., Aremonia agrimonoides (L.) DC., Juglans regia L., Campanula rapunculoides L., Lilium martagon L., Campanula trachelium L., Sanicula europaea L., Arum maculatum L., and Sorbus torminalis (L.) Crantz.) and more xeric species (Helleborus odorus Waldst. & Kit. ex Willd., Scutellaria altissima L., Digitalis grandiflora Mill., Calamintha grandiflora (L.) Moench., Arabis turrita L., and Galium schultesii Vest). These forests belong to the alliance Ostryo carpinifoliae-Tilion platyphylli Košir et al. 2008.
Cluster 3 includes the xeric hop hornbeam forests that belong to the alliance Fraxino orni-Ostryion carpinifoliae Tomažič 1940. Its diagnostic species are Potentilla micrantha Ram. ex DC., Carex cariophyllea Latourr., Glechoma hirsuta Waldst. & Kit., Clinopodium nepeta (L.) Kuntze, and Fraxinus ornus L.
Cluster 4 comprises xerothermophilous forests on screes and shallow skeletal soils. Diagnostic species of the cluster are Pseudofumaria alba subsp. acaulis (Wulfen) Lidén, Stachys recta L., Teucrium montanum L., Prunus spinosa L., Draba lasiocarpa Rochel, Allium flavum L., Evonymus verrucosus Scop., Melica ciliata L., Campanula lingulata Waldst. & Kit., Fritillaria montana Hoppe ex W.D.J.Koch, Minuartia bosniaca (Beck) K.Malý, Clinopodium thymifolium (Scop.) Kuntze, Stipa calamagrostis (L.) Wahlenb., Rhamnus saxatilis Jacq., Coronilla emerus subsp. emeroides (Boiss. & Spuner) Holmboe, and Galium purpureum L. These communities, which are floristically and ecologically clearly separated from other variants of investigated hop hornbeam forests, may be included in a new alliance of Pseudofumario albae-Ostryion carpinifoliae all. nova.
Cluster 5 is represented by the chasmophytic xeric communities in stony habitats. Its diagnostic species are Achillea ageratifolia (Sibth. & Sm.) Boiss, Seseli rigidum Waldst. & Kit., Euphorbia glabriflora Vis., Frangula rupestris Schur, Amphoricarpos neumayeri Vis., Globularia cordifolia L., Hieracium waldsteinii Tausch, Geranium macrorrhizum L., Campanula secundiflora Vis. & Pancic, Genista radiata (L.) Scop., Saxifraga crustata Vest., Valeriana montana L., Cerastium decalvans Schloss. & Vuk., Jurinea mollis (Torn.) Rchb., Saxifraga tridactylites L., and Silene pusilla Waldst. & Kit. We have included communities of the cluster in a new alliance of Achilleo ageratifoliae-Ostryion carpinifoliae all. nova.
These variants of the analyzed forests are distributed in clearly distinguishable habitat types (Figure 4). The mesic alliance Ostryo carpinifoliae-Fagion sylvaticae occurs in shady habitats, on north-facing slopes or at the bottom of canyons, on deep soil (Figure 4a). Forests of the alliance Ostryo carpinifoliae-Tilion platyphylli are also distributed in shady habitats but, unlike the previous alliance, on more skeletal soils (Figure 4b).
Xeric variants of the investigated forests usually occur in dry, sunny habitats. They differ with respect to soil conditions. Forests that belong to the alliance Fraxino orni-Ostryion carpinifoliae occupy sunny habitats, either on gentle slopes, or at the foot of steep slopes, where the soil is relatively deep (Figure 4c). The new alliances proposed in this article (Pseudofumario albae-Ostryion carpinifoliae and Achilleo ageratifoliae-Ostryion carpinifoliae) occur in adverse scree and saxatile habitats (Figure 4d,e).

3.2. Patterns in the Vegetation–Environment Relationship

Forward selection analysis (Table 3) indicates that all of the environmental variables, except for the altitude, are significant for floristic differentiation of the analyzed forests.
A combined light–moisture gradient separates the analyzed forests from the shady and mesic variants (the alliances Ostryo carpinifoliae-Fagion sylvaticae and Ostryo carpinifoliae-Tilion platyphylli) to the sunny and xeric variants (the alliances Fraxino orni-Ostryion carpinifoliae, Pseudofumario albae-Ostryion carpinifoliae, and Achilleo ageratifoliae-Ostryion carpinifoliae). Light and soil pH gradients are correlated. Sunny and xeric variants of forests occur in screes and rocky cliffs. Due to the shallow soil on carbonate and serpentine rock, the soil pH in xeric variants of the forests (the alliances Pseudofumario albae-Ostryion carpinifoliae and Achilleo ageratifoliae-Ostryion carpinifoliae) is alkaline. In shady communities (the alliances Ostryo carpinifoliae-Fagion sylvaticae and Ostryo carpinifoliae-Tilion platyphylli), which are developed on relatively deep soil, the soil reaction is slightly acidic to neutral. The gradients of soil nutrients and soil pH are negatively correlated (Figure 5).
Both forward selection and canonical correspondence analysis indicate that the variants of the analyzed forests are clearly differentiated with respect to shade, moisture, and soil nutrient gradients. The environmental conditions of the investigated communities (Figure 6) additionally support these findings.

3.3. Alpha Diversity

The highest alpha diversity (expressed as Shannon’s entropy) was recorded in mesic alliance Ostryo carpinifoliae-Tilion platyphylli, while the lowest entropy was detected in xeric alliance Fraxino orni-Ostryion carpinifoliae (Figure 7). A similar trend was observed for species richness. Contrary to the alpha diversity and species richness, the values of species equitability were almost identical (values near 1) in all variants of the analyzed forests.

3.4. Beta Diversity

Extremely high values for both species turnover and beta diversity were detected in all variants of the analyzed forests (Figure 8). The ternary graphs, presented in Figure 8, were obtained using Baselga’s method [35]. We performed analyses of beta diversity using an alternative method, described by Podani et al. [37], but obtained essentially the same results, which are not presented in this article.

4. Discussion

The classification of relevés may be performed using a wide spectrum of clustering methods [43,44,45]. Hill et al. [46] emphasized that the number of misclassifications is a key parameter in assessing the analytical power of clustering methods. Agglomerative classification procedures usually produce a huge number of misclassifications [43,46]. The most popular polythetic divisive methods perform well when the first principal axis accounts for most of the variation in a data set. However, if variances of principal axes are similar, these methods may produce serious misclassifications [47,48]. The inability to correct misclassifications is the main drawback of hierarchical classification methods. Non-hierarchical clustering methods, however, enable the allocation of misclassified relevés to their most similar cluster. Therefore, we performed classification of the analyzed communities using K-means clustering and Bayesian classification. These methods are the most powerful variants of non-hierarchical clustering methods [27]. The variance ratio indicates that K-means clustering produces more acceptable results than Bayesian classification. The main drawback of Bayesian classification is a rigid assumption that all variables must be normally distributed. Austin et al. [49] and Karadžić et al. [50] emphasized that the distribution of species along gradients usually deviates from the Gaussian distribution. Contrary to Bayesian classification, K-means clustering is not restricted by the normality assumption.
Using both methods, we revealed that the analyzed forests may be divided into five clusters. In addition to a clear floristic distinction (Table A1), these variants of hop hornbeam forests differ ecologically. The mesic alliances Ostryo carpinifoliae-Fagion sylvaticae and Ostryo carpinifoliae-Tilion platyphylli occur in shady habitats, on north-facing slopes or at the bottom of canyons. They differ with respect to soil conditions. While forests of Ostryo carpinifoliae-Fagion sylvaticae inhabit sites with deep soils, forests of Ostryo carpinifoliae-Tilion platyphyll occur on relatively shallow soils.
The syntaxonomy of the alliances can be represented by the following classification units:
Class: Carpino-Fagetea sylvaticae Jakucs ex Passarge 1968;
Order: Fagetalia sylvaticae Pawłowski 1928;
Alliance: Ostryo carpinifoliae-Fagion sylvaticae Borhidi 1963;
Order: Aceretalia pseudoplatani Moor 1976;
Alliance: Ostryo carpinifoliae-Tilion platyphylli (Košir et al. 2008) Čarni in Willner et al., 2016.
Mucina et al. [1] included the basiphilous beech and mixed fir-beech forests of the Balkan Peninsula in the alliances Aremonio-Fagion (Horvat 1950) Borhidi in Török et al. 1989 (with the suballiances Ostryo-Fagenion and Lonicero alpigenae-Fagenion), Geranio striati-Fagion Gentile 1970, and Fagion sylvaticae Luquet 1926 (all basiphilous beech forests lacking numerous diagnostic species of Aremonio-Fagion and Geranio-Fagion). Such taxonomy may be questioned, because it excludes the Peri-Pannonian submontane lime-beech forests Tilio tomentose-Fagion sylvaticae (Marinšek, Čarni et Šilc 2013) Karadžić 2018. Moreover, it excludes the ravine forests of Fagus sylvatica and Corylus colurna that occur in eastern parts of Serbia and in Bulgaria. The geographically vicariant alliances Ostryo carpinifoliae-Fagion sylvaticae Borhidi 1963 and Fago sylvaticae-Colurnion colurnae Borhidi 1964 represent the most diverse variants of basiphilous beech forests in Serbia [51]. Therefore, we included the mesic forests of hop hornbeam and Fagus sylvatica in the alliance Ostryo carpinifoliae-Fagion sylvaticae, rather than in the alliance Aremonio-Fagion, which includes mesic forests from both canyon and non-ravine habitats.
Xeric variants of the investigated forests usually occur in dry, sunny habitats. They differ with respect to soil conditions. Forests that belong to the alliance Fraxino orni-Ostryion carpinifoliae are distributed on skeletal, but relatively deep, soil. The new alliances proposed in this article (Pseudofumario albae-Ostryion carpinifoliae and Achilleo ageratifoliae-Ostryion carpinifoliae) occur in adverse scree and saxatile habitats. It is the stress gradient (from relatively deep but skeletal soils to screes and bare rocks) that clearly separates these three variants of xeric hop hornbeam forests. The classification of xeric hop hornbeam forests involves the following syntaxa:
Class: Quercetea pubescentis Doing-Kraft ex Scamoni et Passarge 1959;
Order: Quercetalia pubescenti-petraeae Klika 1933;
Alliances: Fraxino orni-Ostryion Tomažič 1940,
Pseudofumario albae-Ostryion carpinifoliae all. nova, and
Achilleo ageratifoliae-Ostryion carpinifoliae all. nova.
Species-rich hop hornbeam forests are widely distributed in central parts of the Balkan Peninsula [3,4,5,6,7,8,16,17,18,19,20,21,22,23]. Due to their extremely heterogenous structure, the syntaxonomy of these forests is still unresolved. Our article clearly indicates that hop hornbeam forests in Serbia may be divided into five ecologically and floristically distinguishable alliances. Similar communities are extended outside the study area. Therefore, a detailed chorological and ecological analysis of associations within each of these alliances requires further vegetation studies. The forward selection analysis indicates that all environmental variables, except for the altitude, are significant for floristic differentiation of the analyzed forests. Most of the analyzed communities are distributed within a narrow altitudinal range (from 300 to 550 m). Therefore, the effects of altitude on the differentiation of communities is insignificant. The specific topography of canyons forms strong thermal and moisture gradients [8]. As canonical correspondence analysis indicates, the investigated ravine forests are clearly differentiated along the moisture and temperature gradients. A strong light gradient is a consequence of the different density of dominating trees in the analyzed forests. The screes and cliffs are covered by sparsely distributed individuals of hop hornbeam and other xeric trees. The shade increases from the extremely xeric alliances Pseudofumario albae-Ostryion carpinifoliae and Achilleo ageratifoliae-Ostryion carpinifoliae to forests of the alliance Ostryo carpinifoliae-Fagion sylvaticae, which are dominated by dense populations of Fagus sylvatica and other broad-leaved trees.
Analyses of the species diversity involved investigations of both alpha (within-community) and beta (between-community) diversity. The alpha diversity depends on both species richness and species equitability (evenness of species abundance) [52,53]. The high alpha diversity of the investigated forests is attributable to the diverse environmental conditions that enable the coexistence of numerous species with different ecological requirements [8,50]. A heterogeneous, patchy environment prevents the dominance of one species and promotes a polydominant community structure [6].
The beta diversity can be partitioned into two additive components: Species turnover and nestedness (the difference in species richness) [35,37,42]. Extremely high—almost the greatest possible—values of both species turnover and beta diversity were found in all variants of the analyzed forests. Such high diversity can be explained by the heterogeneous microclimatic conditions along strong thermal and moisture gradients [8].
The high diversity of phanerophytes, the presence of numerous Tertiary relicts, and a high proportion of paleo-endemic species are the main characteristics of ravine vegetation in the Balkan Peninsula [6,7,8]. The high diversity of phanerophytes suggests that broad-leaved ravine forests are remnants of Tertiary subtropical vegetation [8].
The Balkan Peninsula is one of the most important biodiversity hotspots in Europe [54,55,56]. According to Velčev et al. [57], Balkan endemic taxa can be divided into two main categories: Paleo-endemics and neo-endemics. Paleo-endemics include relict species from the Tertiary period, which dominated Europe at the beginning of the so-called Quaternary period. The high percentage of paleo-endemic species in canyons can be explained by the ‘chasmophytic divergence’ hypothesis [7,8,58]. The hypothesis is based on three main assumptions: (1) A strong selective pressure in adverse rocky habitats; (2) reduced gene flow between populations; and (3) genetic drift.
The expansion of chasmophytes into adjacent non-ravine habitats is prevented by more efficient competitors. Due to dispersion barriers, the chasmophytic populations are isolated and their gene flow is reduced. Even within the same saxatile habitat, the gene flow may be considerably inhibited due to obstructions in the formation of progeny. The failure of seed to establish new individuals occurs frequently, because most of the holes and crevices that are suitable for ecesis are already occupied. Such a situation reduces gene flow significantly. A combination of reduced gene flow, a strong selective pressure, and genetic drift may contribute to the fast speciation of endemic taxa. Examples of the amazing diversification of chasmophytes are genera Edraianthus A. DC. and Amphoricarpos Vis. [59,60].
The effective size of chasmophyte populations is usually very low. In such a situation, the genetic drift stochastically eliminates alleles. The fixation of (often non-adaptive) alleles and increased proportion of homozygous individuals reduce the ecological plasticity of chasmophytes. Therefore, most of the species that are adapted to saxatile habitats belong to the group of critically endangered taxa.
Due to the high diversity, the protection of broad-leaved ravine forests should be the top priority in biodiversity conservation projects [61].

5. Conclusions

We used both K-means clustering and Bayesian classification to obtain maximally homogeneous clusters of relevés. The optimal number of clusters was determined by maximizing the ratio of between-group to within-group variances [28]. The results of both the K-means clustering and Bayesian classification indicate that the analyzed forests may be divided into five variants (alliances Ostryo carpinifoliae-Fagion sylvaticae, Ostryo carpinifoliae-Tilion platyphylli, Fraxino orni-Ostryion carpinifoliae, Pseudofumario albae-Ostryion carpinifoliae, and Achilleo ageratifoliae-Ostryion carpinifoliae). Despite a great overlap of species distributions, these alliances are clearly distinguishable floristically and ecologically.
A combined light–moisture gradient clearly separates the analyzed forests. Due to shallow soil on carbonate and serpentine rock, the soil pH in xeric variants of the forests (the alliances Pseudofumario albae-Ostryion carpinifoliae and Achilleo ageratifoliae-Ostryion carpinifoliae) is alkaline. In shady communities that are developed on relatively deep soil, the soil reaction is slightly acidic to neutral.
The greatest alpha diversity (expressed as Shannon’s entropy) was recorded in mesic alliance Ostryo carpinifoliae-Tilion platyphylli, while the lowest entropy was detected in the xeric alliance Fraxino orni-Ostryion carpinifoliae. The extremely high species turnover and beta diversity in the analyzed forests can be explained by the heterogeneous microclimatic conditions along strong thermal and moisture gradients.
The investigated forests represent an important pool of rare, paleo-endemic species that survived the Quaternary glaciations in ravine refugia.

Author Contributions

Conceptualization, D.S., N.K., and B.K.; methodology, D.S. and B.K.; software, B.K.; validation, D.S., M.M., S.J., P.P., and N.K.; formal analysis, D.S.; investigation, D.S., S.J., N.K., M.M., P.P., and B.K.; resources, D.S.; data curation, D.S., N.K., and S.J.; writing—original draft preparation, D.S.; writing—review and editing, D.S., N.K., and B.K.; visualization, D.S.; supervision, B.K. and N.K.; project administration, M.M.; funding acquisition, M.M. and P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia (No. 451–03–68/2020–14/200007).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The set of originally collected relevés is stored in the FLORA database, and it is available upon request.

Acknowledgments

Editor and three of the anonymous reviewers gave valuable comments for the improvement of this paper and are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Synoptic table of the analyzed forests. Each row of the table specifies the fidelity (expressed as the phi coefficient, multiplied by 100) of a species to the alliances Ostryo carpinifoliae-Fagion sylvaticae (A), Ostryo carpinifoliae-Tilion platyphylli (B), Fraxino orni-Ostryion carpinifoliae (C), Pseudofumario albae-Ostryion carpinifoliae (D), and Achilleo ageratifoliae-Ostryion carpinifoliae (E). Diagnostic species are shaded.
Table A1. Synoptic table of the analyzed forests. Each row of the table specifies the fidelity (expressed as the phi coefficient, multiplied by 100) of a species to the alliances Ostryo carpinifoliae-Fagion sylvaticae (A), Ostryo carpinifoliae-Tilion platyphylli (B), Fraxino orni-Ostryion carpinifoliae (C), Pseudofumario albae-Ostryion carpinifoliae (D), and Achilleo ageratifoliae-Ostryion carpinifoliae (E). Diagnostic species are shaded.
SpeciesAlliances
ABCDE
Asperula taurina L.38....
Epimedium alpinum L.47....
Fraxinus excelsior L.39....
Evonymus europaeus L.39....
Alliaria petiolata (M.Bieb.) Cavara & Grande39....
Heracleum sphondylium L.39....
Primula veris L.30....
Pulmonaria officinalis L.39....
Melampyrum hoermannianum K. Malý44....
Angelica sylvestris L.35....
Galeopsis speciosa Mill.35....
Parietaria officinalis L.373...
Convallaria majalis L.35....
Athyrium filix-femina (L.) Roth43....
Cardamine bulbifera (L.) Crantz291...
Equisetum telmateia Ehrh.30....
Cyclamen purpurascens Mill.53....
Acer pseudoplatanus L.30....
Brachypodium sylvaticum (Huds.) P.Beauv.5011...
Salvia glutinosa L.3822...
Hedera helix L.4537...
Asarum europaeum L.2818...
Tilia platyphyllos Scop. 10227..
Helleborus odorus Waldst. & Kit. ex Willd.1155...
Scutellaria altissima L.1539...
Scrophularia nodosa L. 447...
Aremonia agrimonoides (L.) DC.10264..
Juglans regia L..41...
Cornus mas L. .2992.
Campanula rapunculoides L..2764.
Digitalis grandiflora Mill..362..
Calamintha grandiflora (L.) Moench..2616..
Arabis turrita L..271.1
Lilium martagon L..37...
Rubus hirtus Wald. et Kit..42...
Campanula trachelium L..3314..
Sanicula europaea L..37...
Arum maculatum L..32...
Sorbus torminalis (L.) Crantz..37...
Galium schultesii Vest .40...
Cornus sanguinea L..3510..
Potentilla micrantha Ram. ex DC..729..
Carex cariophyllea Latourr..526..
Glechoma hirsuta Waldst. et Kit..3255.
Lathyrus sylvestris L...292.
Clinopodium nepeta (L.) Kuntze..34..
Geranium sanguineum L..721..
Fraxinus ornus L..153017.
Pseudofumaria alba (Mill.) Lidén 53010
Stachys recta L...1926.
Teucrium montanum L...1129.
Prunus spinosa L...1224.
Draba lasiocarpa Rochel..228.
Allium flavum L...622.
Evonymus verrucosus Scop...32518
Melica ciliata L....3212
Campanula lingulata Waldst. & Kit..8.33.
Fritillaria montana Hoppe ex W.D.J.Koch...27.
Minuartia bosniaca (Beck) K.Malý...27.
Micromeria thymifolia (Scop.) Fritsch ...329
Stipa calamagrostis (L.) P.Beauv. ...29.
Rhamnus saxatilis Jacq....34.
Galium purpureum L....2424
Globularia cordifolia L....1321
Hieracium waldsteinii Tausch....50
Geranium macrorrhizum L.....48
Campanula secundiflora Vis. & Pancic ....39
Genista radiata (L.) Scop.....37
Saxifraga crustata Vest.....34
Valeriana montana L. ....34
Cerastium decalvans Schloss. & Vuk.....30
Jurinea mollis (Torn.) Rchb.....31
Saxifraga tridactylites L....232
Silene pusilla Waldst. et Kit.....31
Achillea ageratifolia (Sibth. & Sm.) Boiss....38
Seseli rigidum Waldst. & Kit.....35
Edraianthus graminifolius (L.) A.DC....437
Euphorbia glabriflora Vis....734
Frangula rupestris Schur....34
Amphoricarpos neumayeri Vis....221

References

  1. Mucina, L.; Bültmann, H.; Dierssen, K.; Theurillat, J.P.; Raus, T.; Čarni, A.; Šumberová, K.; Willner, W.; Dengler, J.; García, R.G.; et al. Vegetation of Europe: Hierarchical floristic classification system of vascular plant, bryophyte, lichen, and algal communities. Appl. Veg. Sci. 2016, 19, 3–264. [Google Scholar] [CrossRef]
  2. Blasi, C.; Filibeck, G.; Rosati, L. Classification of southern Italy Ostrya carpinifolia woods. Fitosociologia 2006, 43, 3–23. [Google Scholar]
  3. Košir, P.; Čarni, A.; di Pietro, R. Classification and phytogeographical differentiation of broad-leaved ravine forests in southeastern Europe. J. Veg. Sci. 2008, 19, 331–342. [Google Scholar] [CrossRef]
  4. Čarni, A.; Košir, P.; Karadžić, B.; Matevski, V.; Redžić, S.; Škvorc, Ž. Thermophilous deciduous forests in Southeastern Europe. Plant Biosyst. 2009, 143, 1–13. [Google Scholar] [CrossRef]
  5. Stupar, V.; Brujić, J.; Škvorc, Ž.; Čarni, A. Vegetation types of thermophilous deciduous forests (Quercetea pubescentis) in the Western Balkans. Phytocoenologia 2016, 46, 49–68. [Google Scholar] [CrossRef]
  6. Karadžić, B.; Mijović, A.; Popović, R.; Perišić, S.; Marinković, S. Forest vegetation in West-Serbian canyons: Biodiversity hotspots. In Forest Research: A Challenge for an Integrated European Approach; Ragodlou, K., Ed.; NAGREF-Forest Research Institute: Thessaloniki, Greece, 2001; pp. 513–518. [Google Scholar]
  7. Karadžić, B. Chasmophytic forests of Ostrya carpinifolia in west-Serbian canyons. Biol. Nyssana 2017, 8, 73–81. [Google Scholar]
  8. Karadžić, B.; Bulić, Z.; Jarić, S.; Mitrović, M.; Pavlović, P. Vegetation in Ravine Habitats of Montenegro. In The Rivers of Montenegro; Pešić, V., Paunović, M., Kostianoy, A.G., Eds.; Springer Nature Switzerland AG: Cham, Switzerland, 2020; pp. 201–230. [Google Scholar]
  9. Popović, R.; Kojić, M.; Karadžić, B. Ecological characteristics of six important submediterranean tree species in Serbia. Bocconea 1996, 5, 431–438. [Google Scholar]
  10. Birks, H.J.B.; Willis, K.J. Alpines, trees, and refugia in Europe. Plant Ecol. Divers. 2008, 1, 147–160. [Google Scholar] [CrossRef]
  11. Hewitt, G.M. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol. J. Linn. Soc. 1996, 58, 247–276. [Google Scholar] [CrossRef]
  12. Hewitt, G.M. Postglacial recolonization of European Biota. Biol. J. Linn. Soc. 1999, 68, 87–112. [Google Scholar] [CrossRef]
  13. Hewitt, G.M. The genetic legacy of the quaternary ice ages. Nature 2000, 405, 907–913. [Google Scholar] [CrossRef] [PubMed]
  14. Hewitt, G.M. Mediterranean peninsulas: The evolution of hotspots. In Biodiversity Hotspots. Distribution and Protection of Conservation Priority Areas; Zachos, F.E., Habel, J.C., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 123–147. [Google Scholar]
  15. Puncer, I.; Zupančič, M. Die ökologische und wirtschaftliche Bedeutung der Ostrya carpinifolia Scop. in Slowenien. Studia Geobot. 1982, 2, 25–32. [Google Scholar]
  16. Tomić, Z. Cenoareal crnog graba (Ostrya carpinifolia Scop.) u Srbiji. Bot. Serbica 1994, 28, 173–182. [Google Scholar]
  17. Jovanović, B. Šuma crnog graba u okolini Titovog Užica (Sesleria variae-Ostryetum ass. n.). Šumarstvo 1972, 25, 3–9. [Google Scholar]
  18. Trinajstić, I.; Cerovečki, Z. O cenoarealu crnoga graba, Ostrya carpinifolia Scop. (Corylaceae) u Hrvatskoj [On coenoareal of hop-hornbeam, Ostrya carpinifolia Scop. (Corylaceae) in Croatia]. Biosistematika 1978, 4, 57–65. [Google Scholar]
  19. Tomić, Z. Fitocenoze crnoga graba (Ostrya carpinifolia Scop.) u Srbiji [Phytocenoses of Hop hornbeam (Ostrya carpinifolia Scop.) in Serbia]. Ph.D. Thesis, University in Belgrade, Belgrade, Serbia, 1980. [Google Scholar]
  20. Lakušić, R.; Pavlović, D.; Redžić, S. Horološko-ekološka i floristička diferencijacija šuma i šikara sa bjelograbićem (Carpinus orientalis Mill.) i crnim grabom (Ostrya carpinifolia Scop.) na prostoru Jugoslavije. Glas. Republičkog Zavoda za Zaštitu Prir. Prir. Muz. Titogr. 1982, 15, 103–116. [Google Scholar]
  21. Tomić, Z. Specijski diverzitet u crnograbovim šumama sveze Orno-Ostryon Tomž. 1940. i njegove karakteristike. Zaštita Prir. 1998, 50, 57–61. [Google Scholar]
  22. Tomić, Z. Zajednica Orno-Ostryetum Aich. 1933 u refugijumima jugozapadne Srbije. Glas. Šumarskog Fak. 2000, 82, 177–189. [Google Scholar]
  23. Tomić, Z. Šume crnog graba (Hop hornbeam forests). In Vegetacija Srbije 2—Šumske zajednice 2; Škorić, D.M., Ed.; Srpska akademija nauka i umetnosti, Odeljenje hemijskih i bioloških nauka: Beograd, Srbija, 2006; pp. 29–68. [Google Scholar]
  24. Braun-Blanquet, J. Plant Sociology: The Study of Plant Communities; Hafner: London, UK, 1965. [Google Scholar]
  25. Van der Maarel, E. Transformation of coverabundance values in phytosociology and its effects on community similarity. Vegetatio 1979, 39, 97–114. [Google Scholar]
  26. MacQueen, J. Some Methods for Classification and Analysis of Multivariate Observations, Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability, Oakland, CA, USA, 27 December 1965–7 January 1966; University of California Press: Los Angeles, CA, USA, 1967; pp. 281–297. [Google Scholar]
  27. James, G.; Witten, D.; Hastie, T.; Tibshirani, R. An Introduction to Statistical Learning: With Applications in R; Springer: New York, NY, USA, 2013. [Google Scholar]
  28. Marinković, N.; Karadžić, B.; Pešić, V.; Gligorović, B.; Grosser, C.; Paunović, M.; Nikolić, V.; Raković, M. Faunistic patterns and diversity components of leech assemblages in karst springs of Montenegro. Knowl. Manag. Aquat. Ecosyst. 2019, 420, 26. [Google Scholar]
  29. Chytrý, M.; Tichý, L.; Holt, J.; Botta-Dukát, Z. Determination of diagnostic species with statistical fidelity measures. J. Veg. Sci. 2002, 13, 79–90. [Google Scholar] [CrossRef]
  30. Ter Braak, C.J. Canonical correspondence analysis: A new eigenvector technique for multivariate direct gradient analysis. Ecology 1986, 67, 1167–1179. [Google Scholar] [CrossRef] [Green Version]
  31. Ellenberg, H.; Leuschner, C. Vegetation Mitteleuropas mit den Alpen, 6th ed.; Ulmer: Stuttgart, Germany, 2010. [Google Scholar]
  32. Pielou, E.C. Population and Community Ecology: Principles and Methods; Gordon and Breach: New York, NY, USA, 1974. [Google Scholar]
  33. Whittaker, R.H. Evolution and measurement of species diversity. Taxon 1972, 21, 213–251. [Google Scholar] [CrossRef] [Green Version]
  34. Anderson, M.J.; Crist, T.O.; Chase, J.M.; Vellend, M.; Inouye, B.D.; Freestone, A.L.; Sanders, N.J.; Cornell, H.V.; Comita, L.S.; Davies, K.F.; et al. Navigating the multiple meanings of beta diversity: A roadmap for the practicing ecologist. Ecol. Lett. 2011, 14, 19–28. [Google Scholar] [CrossRef] [PubMed]
  35. Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 2010, 19, 134–143. [Google Scholar] [CrossRef]
  36. Baselga, A. The relationship between species replacement, dissimilarity derived from nestedness, and nestedness. Glob. Ecol. Biogeogr. 2012, 21, 1223–1232. [Google Scholar] [CrossRef]
  37. Podani, J.; Ricotta, C.; Schmera, D. A general framework for analyzing beta diversity, nestedness and related community-level phenomena based on abundance data. Ecol. Complex. 2013, 15, 52–61. [Google Scholar] [CrossRef] [Green Version]
  38. Schmera, D.; Podani, J. Comments on separating components of beta diversity. Community Ecol. 2011, 12, 153–160. [Google Scholar] [CrossRef]
  39. Cardoso, P.; Rigal, F.; Carvalho, J.C.; Fortelius, M.; Borges, P.A.V.; Podani, J.; Schmera, D. Partitioning taxon, phylogenetic and functional beta diversity into replacement and richness difference components. J. Biogeogr. 2014, 41, 749–761. [Google Scholar] [CrossRef] [Green Version]
  40. Carvalho, J.C.; Cardoso, P.; Gomes, P. Determining the relative roles of species replacement and species richness differences in generating beta-diversity patterns. Glob. Ecol. Biogeogr. 2012, 21, 760–771. [Google Scholar] [CrossRef]
  41. Legendre, P. Interpreting the replacement and richness difference components of beta diversity. Glob. Ecol. Biogeogr. 2014, 23, 1324–1334. [Google Scholar] [CrossRef]
  42. Karadžić, B. FLORA: A software package for statistical analysis of ecological data. Water Res. Manag. 2013, 3, 45–54. [Google Scholar]
  43. Pielou, E.C. The Interpretation of Ecological Data: A Primer on Classification and Ordination; Wiley: New York, NY, USA, 1984. [Google Scholar]
  44. Jongman, R.H.G.; ter Braak, C.J.F.; van Tongeren, O.F.R. Data Analysis in Community and Landscape Ecology; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  45. Legendre, P.; Legendre, L. Numerical Ecology; Elsevier: New York, NY, USA, 2012. [Google Scholar]
  46. Hill, M.O.; Bunce, R.G.H.; Shaw, M.W. Indicator species analysis, a divisive polythetic method of classification and its application to a survey of native pinewoods in Scotland. J. Ecol. 1975, 63, 597–613. [Google Scholar] [CrossRef]
  47. Van Groenewoud, H. The robustness of Correspondence, Detrended Correspondence, and TWINSPAN Analysis. J. Veg. Sci. 1992, 3, 239–246. [Google Scholar] [CrossRef]
  48. Belbin, L.; McDonald, C. Comparing Three Classification Strategies for Use in Ecology. J. Veg. Sci. 1993, 4, 341–348. [Google Scholar] [CrossRef] [Green Version]
  49. Austin, M.P.; Nicholls, A.O.; Doherty, M.D.; Meyers, J.A. Determining species response functions to an environmental gradient by means of a β-function. J. Veg. Sci. 1994, 5, 215–228. [Google Scholar] [CrossRef]
  50. Karadžić, B.; Marinković, S.; Kataranovski, D. Use of the β-function to estimate the skewness of species responses. J. Veg. Sci. 2003, 14, 799–805. [Google Scholar] [CrossRef]
  51. Karadžić, B. Beech forests (order Fagetalia sylvaticae Pawlowski 1928) in Serbia. Bot. Serbica 2018, 42, 91–107. [Google Scholar]
  52. Whittaker, R.H. Communities and Ecosystems; MacMillan: New York, NY, USA, 1970. [Google Scholar]
  53. Lloyd, M.; Ghelardi, R.J. A Table for Calculating the Equitability Component of Species Diversity. J. Anim. Ecol. 1964, 33, 217–225. [Google Scholar] [CrossRef]
  54. Stevanović, V. Analysis of the Central European and Mediterranean orophytic element on the mountains of the West and Central Balkan Peninsula, with special reference to endemics. Bocconea 1996, 5, 77–97. [Google Scholar]
  55. Stevanović, V.; Tan, K.; Petrova, A. Mapping the endemic flora of the Balkans—A progress report. Bocconea 2007, 21, 131–137. [Google Scholar]
  56. Tomović, G.; Niketić, M.; Lakušić, D.; Ranđelović, V.; Stevanović, V. Balkan endemic plants in Central Serbia and Kosovo regions: Distribution patterns, ecological characteristics, and centres of diversity. Bot. J. Linn. Soc. 2014, 176, 173–202. [Google Scholar] [CrossRef] [Green Version]
  57. Velčev, V.; Kožuharov, S.; Ančev, M. Atlas of the Endemic Plants in Bulgaria; Publishing House of the Bulgarian Academy of Sciences: Sofia, Bulgaria, 1992. [Google Scholar]
  58. Davis, H.P. Cliff vegetation in the Eastern Mediterranean. J. Ecol. 1951, 39, 63–93. [Google Scholar] [CrossRef]
  59. Surina, B.; Schönswetter, P.; Schneeweiss, G.M. Quaternary range dynamics of ecologically divergent species (Edraianthus serpyllifolius and E. tenuifolius, Campanulaceae) within the Balkan refugium. J. Biogeogr. 2011, 38, 1381–1393. [Google Scholar] [CrossRef]
  60. Caković, D.; Stešević, D.; Schönswetter, P.; Frajman, B. How many taxa? Spatiotemporal evolution and taxonomy of Amphoricarpus (Asteraceae, Carduoideae) on the Balkan Peninsula. Org. Divers. Evol. 2015, 15, 429–455. [Google Scholar] [CrossRef]
  61. Council Directive 92/43/EEC of 21 May 1992 on the Conservation of Natural Habitats and of Wild Fauna and Flora. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:31992L0043&from=EN (accessed on 1 February 2021).
Figure 1. A map of the study area. The canyons of the Gradac (a), Trešnjica (b), Derventa (c), Djetinja (d), Ljutina (e), Mileševka (f) and Lim (g) rivers are circled.
Figure 1. A map of the study area. The canyons of the Gradac (a), Trešnjica (b), Derventa (c), Djetinja (d), Ljutina (e), Mileševka (f) and Lim (g) rivers are circled.
Diversity 13 00059 g001
Figure 2. Venn diagrams of overlapping species in two communities. These diagrams may explain the concept of beta diversity components. The greatest beta diversity implies the complete replacement of a species from one to another community (a). Two communities may differ not because of species replacement, but due to differences in species richness (c). The beta diversity usually depends on both species turnover and nestedness (b).
Figure 2. Venn diagrams of overlapping species in two communities. These diagrams may explain the concept of beta diversity components. The greatest beta diversity implies the complete replacement of a species from one to another community (a). Two communities may differ not because of species replacement, but due to differences in species richness (c). The beta diversity usually depends on both species turnover and nestedness (b).
Diversity 13 00059 g002
Figure 3. Important endemic species that occur in hop hornbeam forests: Onosma stellulatum Wald. et Kit (a), Lathyrus binatus Pančić (b), Achillea ageratifolia (Sibth. & Sm.) Boiss. (c), Daphne malyana Blečić (d), Athamanta turbith (L.) Brot. subsp. haynaldii (Borbas & Uechtr.) Tutin. (e), Aquilegia grata Zimmeter (f), Campanula secundiflora Vis. & Pančić (g), Euphorbia glabriflora Vis. (h), and Hieracium waldsteinii Tausch (i).
Figure 3. Important endemic species that occur in hop hornbeam forests: Onosma stellulatum Wald. et Kit (a), Lathyrus binatus Pančić (b), Achillea ageratifolia (Sibth. & Sm.) Boiss. (c), Daphne malyana Blečić (d), Athamanta turbith (L.) Brot. subsp. haynaldii (Borbas & Uechtr.) Tutin. (e), Aquilegia grata Zimmeter (f), Campanula secundiflora Vis. & Pančić (g), Euphorbia glabriflora Vis. (h), and Hieracium waldsteinii Tausch (i).
Diversity 13 00059 g003
Figure 4. Habitats of the alliances Ostryo carpinifoliae-Fagion sylvaticae (a), Ostryo carpinifoliae-Tilion platyphylli (b), Fraxino orni-Ostryion carpinifoliae (c), Pseudofumario albae-Ostryion carpinifoliae (d), and Achilleo ageratifoliae-Ostryion carpinifoliae (e).
Figure 4. Habitats of the alliances Ostryo carpinifoliae-Fagion sylvaticae (a), Ostryo carpinifoliae-Tilion platyphylli (b), Fraxino orni-Ostryion carpinifoliae (c), Pseudofumario albae-Ostryion carpinifoliae (d), and Achilleo ageratifoliae-Ostryion carpinifoliae (e).
Diversity 13 00059 g004
Figure 5. Canonical correspondence analysis of the investigated forests.
Figure 5. Canonical correspondence analysis of the investigated forests.
Diversity 13 00059 g005
Figure 6. Environmental conditions in hop hornbeam forests. The blue bars are mean values and the lines are standard deviations of the light intensity, moisture, and soil nutrients in the investigated communities.
Figure 6. Environmental conditions in hop hornbeam forests. The blue bars are mean values and the lines are standard deviations of the light intensity, moisture, and soil nutrients in the investigated communities.
Diversity 13 00059 g006
Figure 7. Alpha diversity components of the analyzed forests. The blue bars are mean values and the lines are standard deviations of the species richness, Shannon’s entropy, and equitability in different variants of the analyzed forests.
Figure 7. Alpha diversity components of the analyzed forests. The blue bars are mean values and the lines are standard deviations of the species richness, Shannon’s entropy, and equitability in different variants of the analyzed forests.
Diversity 13 00059 g007
Figure 8. Beta diversity of the investigated forests. Each point (red circle) represents a pair of communities. Its position is determined by values of the beta diversity and its additive components (species turnover and nestedness). The blue rectangle corresponds to the centroid (average value) of the points.
Figure 8. Beta diversity of the investigated forests. Each point (red circle) represents a pair of communities. Its position is determined by values of the beta diversity and its additive components (species turnover and nestedness). The blue rectangle corresponds to the centroid (average value) of the points.
Diversity 13 00059 g008
Table 1. The geographic coordinates of the study area.
Table 1. The geographic coordinates of the study area.
RiverGeographic CoordinatesNumber of Relevés
Gradac44.188436 N, 19.868344 E–44.233333 N, 19.867778 E23
Trešnjica44.141667 N, 19.531667 E–44.150278 N, 19.540278 E28
Derventa43.954167 N, 19.356944 E–43.963056 N, 19.358333 E21
Djetinja43.848611 N, 19.808611 E–43.844444 N, 19.772778 E18
Ljutina43.475833 N, 19.405556 E–43.501389 N, 19.422500 E20
Mileševka43.364444 N, 19.728889 E–43.358889 N, 19.739722 E15
Lim43.190278 N, 19.766111 E–43.213056 N, 19.764167 E19
Table 2. Dependence of classification results on a pre-selected number of clusters. The greatest ratio of between-group variance (B) to within-group variance (W) assures that the overlap of homogeneous clusters is minimized. In our data set, the greatest variance ratio (bold font) was obtained for five clusters.
Table 2. Dependence of classification results on a pre-selected number of clusters. The greatest ratio of between-group variance (B) to within-group variance (W) assures that the overlap of homogeneous clusters is minimized. In our data set, the greatest variance ratio (bold font) was obtained for five clusters.
Classification Method
K-Means ClusteringBayesian Classification
Number of GroupsBWB/WBWB/W
20.03591.47600.02430.03441.47210.0234
30.05342.14730.02490.05132.14350.0239
40.07142.78890.02560.07212.78650.0259
50.10263.28760.03120.08903.42570.0260
60.10973.85990.02840.08664.08110.0212
70.10944.57250.02390.09344.72420.0198
80.11155.25540.02120.10895.32050.0205
90.11235.95420.01890.11125.96120.0187
100.12216.51670.01870.12016.52230.0184
Table 3. Results of the forward selection analysis.
Table 3. Results of the forward selection analysis.
VariableEigenvalueF StatisticProbability
Soil nutrients0.42145.1010.000 *
Light0.23652.8180.000 *
Temperature0.22522.6820.000 *
Continentality0.17302.0510.000 *
Soil acidity0.13241.5640.000 *
Moisture0.12341.4570.000 *
Aspect0.10871.2820.002 *
Slope0.10251.2080.006 *
Altitude0.07990.9390.326
* Signs indicate the environmental variables significant for florstic differentiation of the ana-lyzed forests.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sekulić, D.; Karadžić, B.; Kuzmanović, N.; Jarić, S.; Mitrović, M.; Pavlović, P. Diversity of Ostrya carpinifolia Forests in Ravine Habitats of Serbia (S-E Europe). Diversity 2021, 13, 59. https://doi.org/10.3390/d13020059

AMA Style

Sekulić D, Karadžić B, Kuzmanović N, Jarić S, Mitrović M, Pavlović P. Diversity of Ostrya carpinifolia Forests in Ravine Habitats of Serbia (S-E Europe). Diversity. 2021; 13(2):59. https://doi.org/10.3390/d13020059

Chicago/Turabian Style

Sekulić, Dimitrije, Branko Karadžić, Nevena Kuzmanović, Snežana Jarić, Miroslava Mitrović, and Pavle Pavlović. 2021. "Diversity of Ostrya carpinifolia Forests in Ravine Habitats of Serbia (S-E Europe)" Diversity 13, no. 2: 59. https://doi.org/10.3390/d13020059

APA Style

Sekulić, D., Karadžić, B., Kuzmanović, N., Jarić, S., Mitrović, M., & Pavlović, P. (2021). Diversity of Ostrya carpinifolia Forests in Ravine Habitats of Serbia (S-E Europe). Diversity, 13(2), 59. https://doi.org/10.3390/d13020059

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