Next Article in Journal
Back from the Past: DNA Barcodes and Morphology Support Ablabesmyia americana Fittkau as a Valid Species (Diptera: Chironomidae)
Previous Article in Journal
Lichen Diversity and Biomonitoring: A Special Issue
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Fate of Endangered Rock Sedge (Carex rupestris) in the Western Carpathians—The Future Perspective of an Arctic-Alpine Species under Climate Change

1
Institute of Botany, Plant Science and Biodiversity Center, Slovak Academy of Sciences, Dúbravská cesta 9, 845 23 Bratislava, Slovakia
2
Department of Biology and General Ecology, Faculty of Ecology and Environmental Sciences, Technical University in Zvolen, Masaryka 24, 960 01 Zvolen, Slovakia
3
Department of Ecosystem Biology, Faculty of Science, University of South Bohemia, Branišovská 1760, 370 05 České Budějovice, Czech Republic
4
Botanical Garden of Comenius University, Botanická 3, 841 04 Bratislava, Slovakia
5
Botanical Garden of Comenius University, Blatnica 315, 038 15 Blatnica, Slovakia
*
Author to whom correspondence should be addressed.
Submission received: 31 July 2019 / Revised: 11 September 2019 / Accepted: 16 September 2019 / Published: 19 September 2019
(This article belongs to the Section Plant Diversity)

Abstract

:
Carex rupestris is an endangered and rare arctic-alpine element of the Western Carpathian flora. Given the geographically isolated and spatially restricted peripheral ranges of arctic-alpine species, there is a good chance that many species of conservation concern irreversibly disappear from the regional flora under the ongoing climate change. In this study, we gathered all existing data on the presence of C. rupestris and focused on its current and future distribution in the Western Carpathians. We found that although the distribution of the species is fragmented and scarce, C. rupestris occurs in several mountain ranges, in four distinct plant community types, which differ considerably in altitude, geological bedrock, and other habitat characteristics. In contrast to the relatively broad range of occupied habitats, C. rupestris shows a narrow temperature niche (mean annual temperature range 0.4–4.0 °C). Ensembles of small models based on climatic characteristics and local topography show that regardless of the climate change scenario (rcp2.6, rcp8.5), many current occurrence sites, mainly in the peripheral zones of the range, will face the excessive loss of suitable environmental conditions. It is expected that the Tatra Mountains will be the only mountain range retaining potentially suitable habitats and providing possible refugia for this cold-adapted species in the future. Such severe shrinkage of distribution ranges and associated geographic isolation raises serious concerns for the fate of the arctic-alpine species in the Western Carpathians.

1. Introduction

Arctic-alpine taxa are mainly distributed in higher latitudes. Due to their adaptation to low temperatures, these species show an extrazonal distribution occupying the highest summits of the European Alpine System (EAS) [1,2], such as the Alps, Carpathians, Sudeten, Pyrenees, Apennines, Dinarides, and Rilo–Rhodopean massifs [3,4,5]. Regarding the extent of their distribution outside the Arctic, there are large differences between taxa. Whereas some arctic-alpine species are common in many mountain ranges of the Northern Hemisphere, others are very rare, and in some cases, only a few fragile populations are known outside the Arctic and the directly connected mountain ranges [6,7]. Small and isolated populations of theses rare cold-adapted species are especially vulnerable to ongoing climate change [8].
Carex rupestris All. (Figure 1A; rock sedge, syn.: Carex petraea Wahlenb.) is an arctic-alpine species, which forms co-dominant or dominant stands in its distributional optimum. This species has a circumpolar distribution and shows an insular distribution pattern in the EAS. It has first reached Western Europe in the Pleistocene from the northern regions of the Holarctic and subsequently migrated eastwards along the chains of the Alps and the Carpathians [9].
C. rupestris is an extremely rare, critically endangered, and protected red-listed species in the whole Carpathian range and adjacent mountain ranges; apart from the Slovak part of the Western Carpathians, two localities are known in Sudeten (Czech Republic) [10,11], a few localities are found in the Polish part of the Eastern Beskids (Poland) [12] and in the Ukraine [13]. Rock sedge is considered to be a relic species [14] as well as an important diagnostic species of the relic arctic-alpine vegetation [3,15]. It has a very specific distribution within the Western Carpathians (Figure 2), what drew the attention of botanists for decades [16,17,18]. It is completely missing in the Polish (northern) part of the Tatra Mountains and its distribution in Slovakia is obviously disjunct, with the habitats differing among the highest mountain range (Vysoké Tatry Mountains) and the other, smaller mountain ranges. However, the exact environmental differentiation of presence sites has not been quantified yet.
The fact that global climate change is already having a significant impact on ecosystems is widely recognized [8,19,20,21]. It might alter natural disturbance regimes, such as drought frequency [22], extreme rainfall events [23], and floods [24], as well as the outbreaks of fire, insect pests, or wind storms in terms of frequency and severity, altering also the corresponding species distribution [25,26]. Arctic and alpine ecosystems are considered to be one of the most vulnerable and susceptible systems to many of the effects of climate change—e.g., increases in temperature, altered precipitation patterns, and changes of seasonal dynamics [27]. There is also an ongoing debate about species range shifts towards higher altitudes and latitudes [28,29,30], with the strongest shifts in species composition in the vegetation of the Tatra Mountains being found in mylonite grasslands of a small plot size [31]. Another striking effect of climate change in mountain ecosystems is the thermophilization of vegetation, where the composition of species accompanying the cold-adapted taxa changed significantly and the cold-adapted plants are recently co-occurring with species preferring warmer conditions than in the past [31,32].
In the Carpathians, long-term climatic data analyses reveal a significant rise of mean, minimum, and maximum temperature over the last decades [33]. Reflecting the geographically isolated and spatially restricted peripheral ranges of arctic-alpine species in the Western Carpathians, there is a good chance that many rare and threatened species irreversibly disappear from the regional flora under ongoing climate change. Declines of alpine species, as well as their local extinctions, have been already noted in the Eastern part of the Carpathians [30], where C. rupestris is listed as extinct at several localities. Almost half of the localities with the occurrence of the cold-adapted rock sedge in the Western Carpathians lie in the peripheral lower mountain ranges, where it is threatened not only by rising temperatures but also by the absence of possibly higher-elevated refugia. The high-risk distribution, as well as the biological traits of the species, such as the non-expansive way of propagation, its relic character, and eye-catching position in the spotlight of botanists over the past decades makes C. rupestris an ideal model system for the study. We hypothesized that under ongoing climate change, many rare plant species, especially cold-adapted arctic-alpine species, on the periphery of their distributional area are in competitive disadvantage against the more competitive, warm-tolerant species, and their life traits predestine their ability to colonize new suitable habitats, what in turn means the restriction of their overall distribution.
The purpose of this paper is: (a) to examine habitat preferences of C. rupestris in the Western Carpathians; (b) to record the current distributional pattern of the species, and (c) to predict potential shifts in the distribution under climate change scenarios.

2. Materials and Methods

2.1. Vegetation Data

The field research was conducted in the Western Carpathians (Slovakia) within the years 2005–2009. We have resurveyed all the known localities with the occurrence of C. rupestris in order to confirm the current presence of the species. All 65 phytosociological relevés from 23 localities were gained following the same standard procedures of the Zürich-Montpellier School [34,35] using the modified 9-degree Braun-Blanquet’s sampling scale [36] and stored in the Slovak Vegetation Database—EU-SK-001 (SVD) [37] or in the private databases of the authors using the TURBOVEG program [38].
The nomenclature of the taxa generally follows the Checklist by Marhold and Hindák [39], taxa not used in this work are quoted with the author citation. The names of syntaxa follow the outline of Kliment et al. [40,41] and Jarolímek and Šibík [42].

2.2. Environmental Data

To examine environmental differences among habitats with C. rupestris occurrence, several environmental variables were recorded or calculated for 62 vegetation plots with the presence of C. rupestris and their known geographic position (Table 1). Data on the elevation, slope, aspect, temperature, precipitations, and solar energy were derived from fine-resolution GIS layers (10 m resolution). A regular grid-based digital terrain model (DTM) was generated from contour lines and height points vectorized from 1:10,000 raster maps in GRASS GIS v6.4 [43]. Aspect values were linearized and rescaled to range from 0 (south), through 2 (east, west) to 4 (north). Solar energy income was calculated according to Hofierka and Šúri [44] using solar geometry [45,46]. Solar energy characteristics included the monthly and yearly sum of insolation time of direct solar radiation, monthly and yearly sum of the direct, diffuse, and reflected solar radiation. Variables were calculated using our r.sunangle Python script created in the GRASS GIS environment with 10-min increments [47]. Air temperature and vertical atmospheric precipitation were calculated from rasters of mean annual precipitations for the years between 1961–1990. The source data were provided by the Slovak Hydrometeorological Institute. Information on geological settings in 1:50,000 scale were acquired from State Geological Institute of Dionýz Štúr and the geological bedrock of sites was classified into three groups: mylonite (includes mylonite and granitic mylonite), marl (includes shale, hornstone, and marl), and lime (includes limestone and dolomite).
A soil sample was collected from 44 plots and they were bulked, homogenized, purified, and air-dried. Subsequently they were analyzed in a certified laboratory (Institute of Botany, Czech Academy of Sciences, Průhonice) following the standard international protocols (N, C, pH); the determination of individual elements was performed using the methods of atomic emission spectroscopy (K), atomic absorption spectroscopy (Mg, Ca); the available phosphorus in soils was determined by photometric method [48].
For the assessment of distribution shifts under climate scenarios, 19 bioclimatic variables with a resolution of ~1km were obtained from the WorldClim database [49]. Current values of the variables (average for 1960–1990) were supplemented with estimates for the years 2050 (average for 2041–2060) and 2070 (average for 2061–2080) based on two possible greenhouse gas scenarios: minimal rcp2.6 [50] and maximal rcp8.5 [51].

2.3. Data Analysis

Phytosociological relevés were processed in JUICE 7.0 [52] and numerical classification was performed using the program SYNTAX 2000 [53] using the ordinal method with Goodman–Kruskal’s coefficient. The resulting vegetation groups were displayed in ordination space using a constrained analysis of principal coordinates (CAP) based on the distance matrix of Goodman–Kruskal’s coefficients [54]. The indicator species analysis [55] was performed to identify taxa with the highest indicator value (IndVal > 0.65), i.e., species with the highest specificity and fidelity to a particular group derived in vegetation classification. Scores of indicator species were displayed in CAP as weighted averages.
We fit a series of simple models to identify environmental differences among vegetation groups. For quantitative environmental variables, a one-way ANOVA followed by a Tukey HSD test were employed. The same approach was used to test for differences in species richness among groups. For qualitative data (geological bedrock), a multinomial model was fitted via neural networks [56] and pairwise post-hoc tests were performed using the same method with Bonferroni adjusted significance levels. Normality, homoscedasticity, and independence of residuals were assessed using diagnostic plots and spatial correlograms. If needed, variables were power- or log-transformed to meet the assumptions of normality and homoscedasticity. In the case of significant autocorrelation patters in the residuals, generalized least square models with a Gaussian spatial correlation structure were employed [57].
Subsequently, we built classification tree models [58] to identify a subset of the environmental variables most strongly associated with a vegetation classification scheme. Trees were fitted using a conditional inference approach with multiple statistical tests to prevent over-fitting [59]. A five-percent significance level was employed for variable selection and as a stopping criterion for growing the trees. The predictive performance of the trees was assessed using cross-validated (leave-one-out) classification accuracy. Separate trees were built for the environmental data from 62 plots and soil data from 44 plots. Since classification trees may be susceptible to collinearity of predictors, we screened the data prior to analysis and removed highly correlated variables—temperature and precipitations (both strongly correlated with altitude), global, direct, diffused, and reflected radiation (all correlated with the duration of solar radiation), pHH2O (strongly correlated with pHKCl), C and Ca content (both strongly correlated with N and P content). Residuals of the trees were checked for autocorrelation using spatial correlograms; no significant spatial autocorrelation pattern was observed.
We built habitat suitability models (HSMs) to estimate the current extent of favorable habitat conditions for the species and to predict the extent of suitable habitats in the year 2050 and 2070 under the two climate change scenarios (see above). Nineteen bioclimatic variables [49] were used as habitat characteristics for HSMs. Since the local topography may create refugia that may not be recognized in purely climate-based models [60], the slope and aspect were therefore included into the modeling. Due to the relatively steep geographic gradient in Slovakia, there was a high degree of collinearity among the habitat characteristics. Therefore, prior to analysis, we removed strongly correlated variables (Pearson |r| > 0.7) in a stepwise manner to diminish redundant information in the dataset. Mean annual temperature, isothermality, precipitation seasonality, slope, and aspect were retained as informative predictors in the HSMs. The data on the presence of C. rupestris were supplemented with true absences obtained from the SVD [37]. Twenty-three sites with known presences of species contrasted sharply with 8697 absences recorded after the year 2000 in the SVD. Due to greatly unbalanced dataset, we used ensembles of small models (ESMs) to calibrate HSMs. This novel approach is especially suitable for the distribution modeling of rare species, as it tends to outperform standard HSMs when a dataset consists of a few occurrences only [61,62]. Briefly, ESMs are built by fitting many small models (typically with 1–3 predictors) and then averaging them into an ensemble model, where the small models are weighted by their cross-validated predictive performance. Here, we used all possible bivariate combinations of the predictor variables and fit ESMs with four frequently used species distribution modeling techniques: artificial neural networks [63], boosted regression trees [64], generalized additive models [65], and random forests [66], yielding 40 bivariate models in total. We followed Breiner et al. [67] and varied the model parameters and complexity separately for each bivariate model and evaluated their predictive performance using a 10-fold cross-validation. The area under the receiver operating characteristic curve (AUC) was used as a quantitative measure of predictive power [68]. For each bivariate model, the parameter setting with the highest cross-validated AUC score was retained as the final optimized bivariate model. These tuned bivariate models were then averaged, weighted by their mean cross-validated AUC scores, to form tuned ESMs. The predictive performance of the ESMs was evaluated using 10-times repeated leave-groups-out cross-validation with 70% of the data used to train the models while 30% was left for validation [68]. Spatial predictions of the ESMs were used to plot habitat suitability maps in QGIS v. 3.6 [69].
The analyses were performed in R [70] using packages ecospat [71], nlme [72] nnet [56], party [59], and vegan [73].

3. Results

3.1. Vegetation Classification of Plots with Carex rupestris

The cluster analysis of all 65 relevés with the occurrence of C. rupestris in the Western Carpathians resulted in classification with four groups of floristically well-differentiated vegetation units, clearly interpretable on the level of higher syntaxa or association level (Figure 3).
Cluster A represents the group of relevés from the alliance Astero-Seslerion (class Elyno-Seslerietea) and encompasses relevés from lower altitudes (Nízke Tatry Mountains, Veľká Fatra Mountains, and Krivánska Malá Fatra Mountains). This group of relevés is typical for presence of (ordered by decreasing indicator value) Minuartia langii, Gentiana clusii, Jovibarba globifera, Sesleria alba, Kernera saxatilis, Thymus pulcherrimus, Primula auricula, Trisetum alpestre, Campanula cochlearifolia, and Festuca tatrae (Figure 4).
Cluster B represents relevés from the alliance Caricion firmae (class Elyno-Seslerietea). It is divided into two sub-clusters, whereas the second one (cluster B2) includes typical relevés of the alliance from the Belianske Tatry Mountains, and the first one (cluster B1) includes relevés from different parts of the Western Carpathians. We did not find any indicator species that would strongly (IndVal > 0.65) distinguish cluster B from the other clusters.
Cluster C comprises clearly differentiated relevés from the alliance Festucion versicoloris (class Carici rupestris-Kobresietea) from high altitudes of the Vysoké Tatry Mountains and one relevé from Západné Tatry Mountains (all of them are from mylonite bedrock). The cluster was typical for presence of Hupersia selago, Festuca supina, Juncus trifidus, Salix reticulata, Saxifraga retusa, Carex fuliginosa, Silene acaulis, Lloydia serotina, Campanula alpina, and Bartsia alpina.
Within cluster D, there are the relevés from the alliance Caricion firmae, which have a very strong relationship to the plant communities of the alliance Oxytropido-Elynion (class Carici rupestris-Kobresietea) that creates a mosaic of both vegetation types, therefore, the species pool for both communities is the same. Arenaria tenella was the only species with sufficient indicator power to distinguish relevés of the cluster D from the other clusters.

3.2. Environmental Discrimination of Vegetation Groups

We tested the environmental differences of vegetation clusters using a series of simple models. Cluster A (alliance Astero-Seslerion) was well separated from all other clusters based on lime bedrock, lower altitudes, and higher temperatures (Table 1). The relevés from the alliance Caricion firmae (cluster B) occur in higher altitudes than those from the alliance Astero-Seslerion. Their soil characteristics are similar except for the content of magnesium, which is lower. Cluster C (alliance Festucion versicoloris) was a well-differentiated group typical for the highest rates of mean annual precipitation and lowest mean annual temperature. The plots occur on the most acidic soils, with the lowest content of nitrogen, carbon, phosphorus, and minerals. Cluster D represents stands from the alliance Caricion firmae from the highest altitudes, which show ecological transitions to the relevés from the alliance Oxytropido-Elynion and often forms contact vegetation. They differ from the typical stands of the alliance Caricion firmae by more acidic soils, milder slopes, higher rates of mean annual precipitation, and lower mean annual temperature [15,74].
The classification tree based on the environmental data showed that the vegetation clusters can be well separated using the information on geological bedrock and altitude (Figure 5). The tree model had a high classification power on out-of-sample data (cross-validated accuracy = 91%). In particular, cluster A was typical for lower altitudes (≤ 1552 m) and cluster C for mylonite bedrock, and were clearly separated. Classification tree based on soil data (not displayed) had considerably lower classification accuracy (68%) and showed a good discrimination power only for cluster C (pHKCl ≤ 6.3).

3.3. Distribution of Carex rupestris

An ensemble of 40 tuned bivariate models showed a high accuracy in predicting the occurrence of C. rupestris (cross-validated AUC ± SD = 0.90 ± 0.02). Predicted suitable habitat (probability of occurrence > 0.5) covers about 1.7% of the area of Slovakia and includes the highest part of the Western Carpathians (Tatra Mountains, Nízke Tatry Mountains, Malá and Veľká Fatra Mountains), which is consistent with the known geographic range of the species (Figure 6). More than 90% of the occurrence records matched with the two highest probability classes for the potential distribution (probability of occurrence > 0.8).
Suitable habitat proportion is expected to decline steeply under all investigated climate change scenarios (Figure 6). Even the mildest rcp2.6 alternative expects a shrinkage of suitable habitat to be 40% in 2050 and the almost complete disappearance of suitable conditions at the known occurrence sites in Nízke Tatry Mountains, Malá, and Veľká Fatra Mountains until 2050. Under rcp8.5, it is expected that there will be an extensive loss of environmental suitability (a decrease to 17% and 6% of the original suitable habitat extent until 2050 and 2070, respectively) and the vanishing of favorable environmental conditions at more than 90% of sites with currently known presence of C. rupestris. The Tatra Mountains is expected to be the only mountain range retaining potentially suitable habitats in 2070.

4. Discussion

4.1. Vegetation Classification and Environmental Differentiation of Plots with Carex rupestris

We found that C. rupestris in the Western Carpathians occurs in four vegetation types with significantly different environmental settings: alliance Astero-Seslerion (cluster A), alliance Caricion firmae (cluster B), transition community from Caricion firmae to the alliance Oxytropido-Elynion (cluster D), and alliance Festucion versicoloris (cluster C). The range of communities span relatively wide gradient of altitudes (1297–2012 m asl), bedrock (from lime to mylonite), and associated soil conditions, such as pH (4.0–7.7) or nitrogen (0.2–2.7%). Apparently, the species is able to tolerate a relatively broad spectrum of habitats. On the other hand, occurrence sites cover a relatively narrow window of temperatures (mean annual temperature 0.4–4.0 °C) suggesting cold-adaptation and a restricted thermal niche of the species.
The geographical distributions of the cold-adapted arctic-alpine species with narrow thermal niches are often very constrained in extrazonal areas. In fact, C. rupestris was considered to be one of the rarest plant species in the former Czechoslovakia [75]. It is believed that the species migrated along the Carpathians from the northern regions of the Holarctic in the Pleistocene [9] and survived the last glacial maximum in a way leading to a disjunctive arctic-alpine occurrence with a distinctive pattern with steep mountain peaks serving as a barrier to migration and often harboring glacial relicts [76].
The distributions of the cold-adapted arctic-alpine species are probably further limited by local evolution, by historical factors, and random events [77], although recent studies based on the measures of palynological richness assume no evident sensitivity of biodiversity at high elevations to past climate changes during the late-glacial/early Holocene period [78]. It seems that the most important limiting factor affecting the distribution of arctic-alpine species is the existence of suitable sites for these species and vegetation during the climatic optimum in the Holocene [79]. On the other hand, the existence of suitable habitats might not be sufficient, as this would assume that species could easily colonize the localities with appropriate conditions. Hence, the species traits controlling its dispersal ability (e.g., the length of pre-flowering and flowering stages, seed productivity, seed size, way of dispersal, clonality, seed germination, mortality of seedlings, need for suitable micro-loci for recruitment) should be considered as well. C. rupestris is a tufted graminoid perennial herb with an extensive, densely branched rhizome, which serves as an efficient local vegetative reproduction. Inflorescence comprises of one terminal spike with wind-pollinated flowers. Seeds show very low germination rates [80]. Moreover, many individuals remain sterile, the size of the populations does not increase and some populations are heavily infected by rust [10,81]. All these characteristics might be considered a competitive disadvantage resulting in a weak ability to colonize new arising habitats.
There was a strong relation between the soil properties and bedrock type of the stands with C. rupestris, which is probably strengthened by the fact that much of the rock in such relatively shallow soils is derived from the underlying bedrock brought near the surface by frost action [82]. The soil reaction is often a key factor influencing the species richness of the high-altitude vegetation and the presence of alpine species, as has been observed in the Alps [83] and other parts of Europe, North America, and Arctic [84,85,86,87,88].

4.2. The Fate of Carex rupestris in the Western Carpathians

The bioclimatic and topologic variables included in ESMs efficiently described the current distribution pattern of C. rupestris and predicted excessive loss of suitable environmental conditions under both of the investigated climate change scenarios. Although some suitable habitats of C. rupestris are predicted to occur even under a pessimistic rcp8.5 scenario, it seems that many current presence sites of the species will lack favorable conditions in the future. Remembering the Ukrainian story [30], even the current elevation of temperatures might inevitably cause the decline or even the extinction of C. rupestris in some current localities, mainly in the peripheral zones of its current distribution. A progressive decline of cold mountain habitats and the thermophilization of mountain plant communities mirroring the degree of recent warming have been already observed also in many European mountain systems [8,31,32].
Although, the highest probability classes of the ESMs further include a number of grid cells for which the occurrence of the species is very likely, the colonization of new suitable habitats is merely hypothetical, reflecting the species traits (see above), isolation of peripheral localities and the size and vitality of the populations. Still, one should bear in mind the limitations of climate-based HSMs. Notably, the models are blind to the ability of a species to track suitable environments, while the understanding of migration processes, the dynamics of shifting populations, and inter-specific interactions may be crucial for accurate spatio-temporal projections of HSMs [89]. Also, direct human activities, such as trampling (Figure 1C), reduce the species richness, diversity, as well as the percentage cover of cold-adapted species in the alpine environment [90]. Heavier disturbances could, therefore, lead to the direct destruction of the alpine environment, significantly eliminate suitable habitats of the species, and increase the risk for upward shifts of warm-tolerant species. The projected decline of C. rupestris in the Western Carpathians may be further exacerbated by the competitive overgrowing of larger plants, which are expected to shift their ranges towards higher altitudes, following elevated temperatures [76]. Pinus mugo (Figure 1C) is considered as one of such plants, significantly changing its frequency and overgrowing sites that were historically covered by high-altitude grasslands [31]. Finally, as stems from the concept of an extinction debt [91,92], the effects of global environmental changes, especially in small populations of long-lived perennials might be underestimated and may cause the delayed extinction of even recently stable-looking populations.
From the spatio-temporal projections of HSM, it is clear that C. rupestris will be restricted to higher elevations and it is expectable that the Tatra Mountains will be the only mountain range retaining potentially suitable habitats and providing possible refugia for this cold-adapted species in the future. Thus, these core refugia are of high conservation concern and deserve strict protection excluding any human actions that might further impair the habitats. Since this is the first attempt to explicitly model the spatial extent of suitable conditions for a rare arctic-alpine species in the Western Carpathians, the temporal trend in other relic arctic-alpine species needs to be estimated to inform the prioritization of conservation targets. Nevertheless, the predicted severe shrinkage of distribution ranges and associated geographic isolation raises serious concerns for the fate of the arctic-alpine species in the Western Carpathians.

Author Contributions

Conceptualization, I.S. and J.Š.; methodology, I.S. and J.Š.; data acquisition, I.S., A.P., D.B., and J.Š.; formal analysis, I.S. and M.S.; data curation, J.Š.; writing—original draft preparation, I.S.; writing—contributed writing the manuscript, M.S., D.S., and J.Š.; writing—review and editing, I.S. and M.S.; visualization, I.S. and M.S.; funding acquisition, M.S. and J.Š.; supervision, M.S. and J.Š.

Funding

This study was supported by the grant agency APVV 16-0431, APVV-15-0210, VEGA 2/0032/17 and VEGA 2/0135/16. The Geographic Resources Analysis Support System (GRASS) and other calculations were performed in the Computing Centre of the Slovak Academy of Sciences using the infrastructure acquired within the projects ITMS 26230120002 and ITMS 26210120002 (Slovak infrastructure for high-performance computing), supported by the Research & Development Operational Program funded by the ERDF.

Acknowledgments

The authors are thankful to Z. Dúbravcová, D. Dítě, M. Perný, and J. Wild for their pleasant fieldwork collaboration and M. Valachovič for his tutorials during the PhD study of the first author. The Institute of Botany of the Czech Academy of Sciences kindly provided the services of their certified soil laboratory for the soil analyses and we would like to thank especially O. Rauch for his support. The authors are grateful also to three anonymous reviewers for their very valuable comments. This work was dedicated to little † Dominik Daniel.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ozenda, P. La Végétation de la Chaîne Alpine dans l’Espace Montagnard Européen, 1st ed.; Masson: Paris, France, 1985; p. 344. [Google Scholar]
  2. Oriolo, G. Naked rush swards (Oxytropido-Elynion Br.-Bl. 1949) on the Alps and the Apennines and their syntaxonomical position. Fitosociologia 2001, 38, 91–101. [Google Scholar]
  3. Petrík, A.; Šibík, J.; Valachovič, M. The class Carici rupestris-Kobresietea bellardii Ohba 1974 also in the Western Carpathians. Hacquetia 2005, 4, 33–51. [Google Scholar]
  4. Hájková, P.; Hájek, M.; Apostolova, I. Diversity of wetland vegetation in the Bulgarian high mountains, main gradients and context-dependence of the pH role. Plant Ecol. 2006, 184, 111–130. [Google Scholar] [CrossRef]
  5. Dítě, D.; Peterka, T.; Dítětová, Z.; Hájková, P.; Hájek, M. Arcto-alpine species at their niche margin: The Western Carpathian refugia of Juncus castaneus and J. triglumis in Slovakia. Ann. Bot. Fennici 2017, 54, 67–82. [Google Scholar] [CrossRef]
  6. Schönswetter, P.; Popp, M.; Brochmann, C. Rare arctic-alpine plants of the European Alps have different immigration histories: The snow bed species Minuartia biflora and Ranunculus pygmaeus. Mol. Ecol. 2006, 15, 709–720. [Google Scholar] [CrossRef] [PubMed]
  7. Hultén, E.; Fries, M. Atlas of North European Vascular Plants, 1st ed.; Koeltz—Scientific Books: Königstein, Germany, 1986; p. 1172. [Google Scholar]
  8. Gottfried, M.; Pauli, H.; Futschik, A.; Akhalkatsi, M.; Barančok, P.; Alonso, J.L.B.; Coldea, G.; Dick, J.; Erschbamer, B.; Calzado, M.R.F.; et al. Continent-wide response of mountain vegetation to climate change. Nat. Clim. Chang. 2012, 2, 111–115. [Google Scholar] [CrossRef]
  9. Pawłowski, B. Dwie ciekawe turzyce z Czarnej Hory. Spraw. Kom. Fizjogr. 1931, 65, 145–152. [Google Scholar]
  10. Holub, J.; Grulich, V. Carex rupestris All. In Červená kniha ohrozených a vzácnych druhov rastlín a živočíchov SR a ČR 4. Vyššie rastliny, 1st ed.; Čeřovský, J., Feráková, V., Holub, J., Maglocký, Š., Procházka, F., Eds.; Príroda: Bratislava, Slovakia, 1999; p. 81. [Google Scholar]
  11. Eliáš jun, P.; Dítě, D.; Kliment, J.; Hrivnák, R.; Feráková, V. Red list of ferns and flowering plants of Slovakia, 5th ed. Biologia 2015, 70, 218–228. [Google Scholar]
  12. Mitka, J.; Michalik, S. Turzyca skalna. In Czerwona Księga Karpat Polskich Rośliny naczyniowe; Mirek, Z., Piękoś-Mirkowa, H., Eds.; Instytut Botaniki im. W. Szafera Polskiej Akademii Nauk: Kraków, Poland, 2008; pp. 526–527. [Google Scholar]
  13. Didukh, Y.P. (Ed.) Chervona Knyha Ukrayiny: Roslynnyi svit; Globalconsulting: Kyiv, Ukraine, 2009; p. 589. [Google Scholar]
  14. Dítě, D.; Hájek, M.; Svitková, I.; Košuthová, A.; Šoltés, R.; Kliment, J. Glacial-relict symptoms in the Western Carpathian flora. Folia Geobot. 2018, 53, 277–300. [Google Scholar] [CrossRef]
  15. Petrík, A.; Dúbravcová, Z.; Jarolímek, I.; Kliment, J.; Šibík, J.; Valachovič, M. Syntaxonomy and ecology of plant communities of the Carici rupestris-Kobresietea bellardii in the Western Carpathians. Biologia 2006, 61, 393–412. [Google Scholar] [CrossRef]
  16. Bernátová, D.; Petrík, A. Ďalšie lokality Carex rupestris All. na Slovensku. Biologia 1983, 38, 491–494. [Google Scholar]
  17. Bernátová, D.; Uhlířová, J. Doplnok k rozšíreniu ostrice skalnej (Carex rupestris All.) vo Veľkej Fatre. Biologia 1986, 41, 75–77. [Google Scholar]
  18. Bernátová, D.; Uhlířová, J. Globulario cordifoliae-Caricetum humilis ass. nova in the rocky part of the Veľká Fatra Mts. Biologia 1994, 49, 1–11. [Google Scholar]
  19. IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, 1st ed.; Cambridge University Press: New York, NY, USA, 2013; p. 1535. [Google Scholar]
  20. Li, D.; Wu, S.; Liu, L.; Zhang, Y.; Li, S. Vulnerability of the global terrestrial ecosystems to climate change. Glob. Chang. Biol. 2018, 24, 4095–4106. [Google Scholar] [CrossRef] [PubMed]
  21. Pederson, N.; Dyer, J.M.; McEwan, R.W.; Hessl, A.E.; Mock, C.J.; Orwig, D.A.; Rieder, H.E.; Cook, B.I. The legacy of episodic climatic events in shaping temperate, broadleaf forests. Ecol. Monogr. 2014, 84, 599–620. [Google Scholar] [CrossRef]
  22. Spinoni, J.; Antofie, T.; Barbosa, P.; Bihari, Z.; Lakatos, M.; Szalai, S.; Szentimrey, T.; Vogt, J. An overview of drought events in the Carpathian Region in 1961–2010. Adv. Sci. Res. 2013, 10, 21–32. [Google Scholar] [CrossRef]
  23. Bartholy, J.; Pongrácz, R.; Torma, C.; Pieczka, I.; Kardos, P.; Hunyady, A. Analysis of regional climate change modelling experiments for the Carpathian basin. Int. J. Glob. Warm. 2009, 1, 238–252. [Google Scholar] [CrossRef]
  24. Palamarchuk, L.; Shpyg, V.; Krakovskaia, S. Floods in the Carpathians: Synoptic analysis and numerical modeling. Geophys. Res. Abstr. 2005, 7, 00967. [Google Scholar]
  25. Navarro, L.; Morin, H.; Bergeron, Y.; Girona, M.M. Changes in spatiotemporal patterns of 20th century spruce budworm outbreaks in eastern Canadian boreal forests. Front. Plant Sci. 2018. [Google Scholar] [CrossRef]
  26. Seidl, R.; Schelhaas, M.J.; Rammer, W.; Verkerk, P.J. Increasing forest disturbances in Europe and their impact on carbon storage. Nat. Clim. Chang. 2014, 4, 806–810. [Google Scholar] [CrossRef] [Green Version]
  27. Ernakovich, J.G.; Hopping, K.A.; Berdanier, A.B.; Simpson, R.T.; Kachergis, E.J.; Steltzer, H.; Wallenstein, M.D. Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. Glob. Chang. Biol. 2014, 20, 3256–3269. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, I.C.; Hill, J.K.; Ohlemuller, R.; Roy, D.B.; Thomas, C.D. Rapid range shifts of species associated with high levels of climate warming. Science 2011, 333, 1024–1026. [Google Scholar] [CrossRef] [PubMed]
  29. Niskanen, A.K.J.; Niittynen, J.; Aalto, J.; Väre, H.; Luoto, M. Lost at high latitudes: Arctic and endemic plants under threat as climate warms. Divers. Distrib. 2019, 25, 809–821. [Google Scholar] [CrossRef]
  30. Kobiv, Y. Response of rare alpine plant species to climate change in the Ukrainian Carpathians. Folia Geobot. 2017, 52, 217–226. [Google Scholar] [CrossRef]
  31. Czortek, P.; Kapfer, J.; Delimat, A.; Eycott, A.; Grytnes, J.A.; Orczewska, A.; Ratyńska, H.; Zięba, A.; Jaroszewicz, B. Plant species composition shifts in the Tatra Mts as a response to environmental change: A resurvey study after 90 years. Folia Geobot. 2018. [Google Scholar] [CrossRef]
  32. Evangelista, A.; Frate, L.; Carranza, M.; Attorre, F.; Pelino, G.; Stanisci, A. Changes in composition, ecology and structure of high-mountain vegetation: A re-visitation study over 42 years. AoB Plants 2016, 8, plw004. [Google Scholar] [CrossRef] [PubMed]
  33. Spinoni, J.; Szalai, S.; Szentimrey, T.; Lakatos, M.; Bihari, Z.; Nagy, A.; Németh, Á.; Kovács, T.; Mihic, D.; Dacic, M.; et al. Climate of the Carpathian Region in the period 1961–2010: Climatologies and trends of 10 variables. Int. J. Clim. 2015, 35, 1322–1341. [Google Scholar] [CrossRef]
  34. Braun-Blanquet, J. Pflanzensoziologie. Grundzüge der Vegetationskunde, 3rd ed.; Springer: Wien, Austria, 1964; p. 631. [Google Scholar]
  35. Westhoff, V.; van den Maarel, E. The Braun-Blanquet approach. In Classification of Plant Communities, 1st ed.; Whittaker, R.H., Ed.; W. Junke: The Hague, The Netherlands, 1978; pp. 289–399. [Google Scholar]
  36. Barkman, J.J.; Doing, H.; Segal, S. Kritische Bemerkungen und Vorschläge zur quantitativen Vegetationsanalyse. Acta Bot. Neerl. 1964, 13, 394–419. [Google Scholar] [CrossRef]
  37. Šibík, J. Slovak Vegetation Database. In Vegetation Databases for the 21st Century; Biodiversity & Ecology; Dengler, J., Oldeland, J., Jansen, F., Chytrý, M., Ewald, J., Finckh, M., Glöckler, F., Lopez-Gonzalez, G., Peet, R.K., Schaminée, J.H.J., Eds.; BEE Press: Hamburg, Germany, 2012; Volume 4, p. 429. [Google Scholar]
  38. Hennekens, S.M.; Schaminée, J.H.J. TURBOVEG, a comprehensive data base management system for vegetation data. J. Veg. Sci. 2001, 12, 589–591. [Google Scholar] [CrossRef]
  39. Marhold, K.; Hindák, F. (Eds.) Zoznam Nižších a Vyšších Rastlín Slovenska, 1st ed.; Veda: Bratislava, Slovakia, 1998; p. 688. [Google Scholar]
  40. Kliment, J.; Bernátová, D.; Jarolímek, I.; Petrík, A.; Šibík, J.; Uhlířová, J. Elyno-Seslerietea. In Rastlinné Spoločenstvá Slovenska. Vysokohorská Vegetácia, 1st ed.; Kliment, J., Valachovič, M., Eds.; Veda: Bratislava, Slovakia, 2007; Volume 4, pp. 147–208. [Google Scholar]
  41. Kliment, J.; Šibík, J.; Šibíková, I.; Dúbravcová, Z.; Jarolímek, I.; Uhlířová, J. High-altitude vegetation of the Western Carpathians—A syntaxonomical review. Biologia 2010, 65, 965–989. [Google Scholar] [CrossRef]
  42. Jarolímek, I.; Šibík, J. (Eds.) Diagnostic, Constant and Dominant Taxa of the Higher Vegetation Units of Slovakia, 1st ed.; Veda: Bratislava, Slovakia, 2008; p. 332. [Google Scholar]
  43. GRASS Development Team. Geographic Resources Analysis Support System (GRASS) Software; Version 6.4.0; Open Source Geospatial Foundation: Chicago, IL, USA, 2010. [Google Scholar]
  44. Hofierka, J.; Šúri, M. The solar radiation model for Open source GIS: Implementation and applications. In Proceedings of the Open Source Free Software GIS-GRASS Users Conference 2002, Trento, Italy, 11–13 September 2002. [Google Scholar]
  45. Krcho, J. Zovšeobecnenie rovnice izalumklín na topografickej ploche a v jej skalárnom poli. Geogr. Časopis 1967, 19, 107–129. [Google Scholar]
  46. Scharmer, K.; Greif, J. (Eds.) Database and exploitation software. In The European Solar Radiation Atlas, 1st ed.; Presses des Mines: Paris, France, 2000; Volume 2, p. 296. [Google Scholar]
  47. Mereďa, P., Jr.; Kučera, J.; Marhold, K.; Senko, D.; Slovák, M.; Svitok, M.; Šingliarová, B.; Hodálová, I. Ecological niche differentiation between tetra- and octoploids of Jacobaea vulgaris. Preslia 2016, 88, 113–136. [Google Scholar]
  48. Olsen, R.S. Estimation of Available Phosphorus in Soils by Extraction with Sodium Bicarbonate; US Dept. of Agriculture: Washington, DC, USA, 1954; p. 19. [Google Scholar]
  49. Hijmans, R.J.; Cameron, S.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Clim. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  50. Weyant, J.; Azar, C.; Kainuma, M.; Kejun, J.; Nakicenovic, N.; Shukla, P.R.; La Rovere, E.; Yohe, G. Report of 2.6 versus 2.9 Watts/m2 RCPP Evaluation Panel; Intergovernmental Panel on Climate Change Secretariat: Geneva, Switzerland, 2009; p. 53. [Google Scholar]
  51. Riahi, K.; Rao, S.; Krey, V.; Cho, C.; Chirkov, V.; Fischer, G.; Kindermann, G.; Nakicenovic, N.; Rafaj, P. RCP8.5—A scenario of comparatively high greenhouse gas emissions. Clim. Chang. 2011, 109, 33–57. [Google Scholar] [CrossRef]
  52. Tichý, L. JUICE, software for vegetation classification. J. Veg. Sci. 2002, 13, 451–453. [Google Scholar] [CrossRef]
  53. Podani, J. Syn-tax 2000. In Computer Program for Data Analysis in Ecology and Systematics for Windows 95, 98 & NT. User’s Manual, 1st ed.; Scientia Publication: Budapest, Hungary, 2001; p. 53. [Google Scholar]
  54. Anderson, M.J.; Willis, T.J. Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology 2003, 84, 511–525. [Google Scholar] [CrossRef]
  55. Dufrêne, M.; Legendre, P. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol. Monogr. 1997, 67, 345–366. [Google Scholar] [CrossRef]
  56. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: New York, NY, USA, 2002; p. 495. [Google Scholar]
  57. Pinheiro, J.C.; Bates, D.M. Mixed-Effects Models in S and S-PLUS. Springer: New York, NY, USA, 2000; p. 528. [Google Scholar]
  58. De’ath, G.; Fabricius, K.E. Classification and regression trees: A powerful yet simple technique for ecological data analysis. Ecology 2000, 81, 3178–3192. [Google Scholar] [CrossRef]
  59. Hothorn, T.; Hornik, K.; Zeileis, A. Unbiased recursive partitioning: A conditional inference framework. J. Comput. Graph. Stat. 2006, 15, 651–674. [Google Scholar] [CrossRef]
  60. Austin, M.P.; Van Niel, K.P. Improving species distribution models for climate change studies: Variable selection and scale. J. Biogeogr. 2010, 38, 1–8. [Google Scholar] [CrossRef]
  61. Lomba, A.; Pellissier, L.; Randin, C.; Vicente, J.; Moreira, F.; Honrado, J.; Guisan, A. Overcoming the rare species modelling paradox: A novel hierarchical framework applied to an Iberian endemic plant. Biol. Conserv. 2010, 143, 2647–2657. [Google Scholar] [CrossRef]
  62. Breiner, F.T.; Guisan, A.; Bergamini, A.; Nobis, M.P. Overcoming limitations of modelling rare species by using ensembles of small models. Methods Ecol. Evol. 2015, 6, 1210–1218. [Google Scholar] [CrossRef]
  63. Ripley, B.D. Pattern Recognition and Neural Networks; Cambridge University Press: New York, NY, USA, 1996; p. 415. [Google Scholar]
  64. Elith, J.; Leathwick, J.R.; Hastie, T. A working guide to boosted regression trees. J. Anim. Ecol. 2008, 77, 802–813. [Google Scholar] [CrossRef] [PubMed]
  65. Wood, S.N. Generalized Additive Models: An Introduction with R, 2nd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2017; p. 496. [Google Scholar]
  66. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  67. Breiner, F.T.; Nobis, M.P.; Bergamini, A.; Guisan, A. Optimizing ensembles of small models for predicting the distribution of species with few occurrences. Methods Ecol. Evol. 2018, 9, 802–808. [Google Scholar] [CrossRef]
  68. Kuhn, M.; Johnson, K. Applied Predictive Modeling; Springer: New York, NY, USA, 2013; p. 600. [Google Scholar]
  69. QGIS Development Team. QGIS Geographic Information System, Version 3.6.0.; Open Source Geospatial Foundation Project: Beaverton, OR, USA, 2019. [Google Scholar]
  70. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  71. Broennimann, O.; Di Cola, V.; Guisan, A. Ecospat: Spatial Ecology Miscellaneous Methods; Version 3.0. R Package; 2018; Available online: https://rdrr.io/cran/ecospat/ (accessed on 18 September 2019).
  72. Pinheiro, J.; Bates, D.; DebRoy, S.; Sarkar, D.; R Core Team. nlme: Linear and Nonlinear Mixed Effects Models; Version 3.1-141; R Package; 2019; Available online: https://cran.r-project.org/web/packages/nlme/citation.html (accessed on 18 September 2019).
  73. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P. vegan: Community Ecology Package; Version 2.3-5; R Package; 2017; Available online: https://rdrr.io/rforge/vegan/ (accessed on 18 September 2019).
  74. Šibík, J.; Petriík, A.; Kliment, J. Syntaxonomical revision of plant communities with Carex firma and Dryas octopetala (alliance Caricion firmae) in the Western Carpathians. Pol. Bot. J. 2004, 49, 181–202. [Google Scholar]
  75. Holub, J. Carex rupestris All. Sborn. Národ. Muz. 1972, 27, 25–26. [Google Scholar]
  76. Fudali, E.; Kučera, J. Bryogeographical elements of moss flora in glacial cirques “Śnieżne Kotły” (Karkonosze Mts) and their threat. Acta Soc. Bot. Pol. 2003, 72, 79–85. [Google Scholar] [CrossRef]
  77. Birks, H.H. The late-quaternary history of arctic and alpine plants. Plant Ecol. Divers. 2008, 1, 135–146. [Google Scholar] [CrossRef]
  78. Feurdean, A.; Tămaş, T.; Tanţău, I.; Fărcaş, S. Elevational variation in regional vegetation responses to late-glacial climate changes in the Carpathians. J. Biogeogr. 2012, 39, 258–271. [Google Scholar] [CrossRef]
  79. Kliment, J.; Šibíková, I.; Šibík, J. On the occurrence of the arctic-alpine and endemic species in the high-altitude vegetation of the Western Carpathians. Thaiszia J. Bot. 2011, 21, 45–60. [Google Scholar]
  80. Alsos, I.G.; Müller, E.; Eidesen, P.B. Germinating seeds or bulbils in 87 of 113 tested Arctic species indicate potential for ex situ seed bank storage. Polar Biol. 2013, 36, 819–830. [Google Scholar] [CrossRef] [Green Version]
  81. Chlebicki, A. Grzyby mikroskopijne na reliktowych stanowiskach turzycy skalnej Carex rupertris v Karpatach. Rocz. Bieszcz. 2000, 9, 61–68. [Google Scholar]
  82. Searcy, K.B.; Wilson, B.F.; Fownes, J.H. Influence of bedrock and aspect on soils and plant distribution in the Holyoke Range, Massachusetts. J. Torrey Bot. Soc. 2003, 130, 158–169. [Google Scholar] [CrossRef]
  83. Vonlanthen, C.M.; Kammer, P.M.; Eugster, W.; Bühler, A.; Veit, H. Alpine vascular plant species richness: The importance of daily maximum temperature and pH. Plant Ecol. 2006, 184, 13–25. [Google Scholar] [CrossRef]
  84. Lunde, T. An investigation into the pH-amplitude of some mountain plants in the county of Troms. Acta Boreal. 1962, 20, 1–103. [Google Scholar]
  85. Gensac, P. Plant and soil groups in the alpine grasslands of the Vanoise Massif, French Alps. Arct. Alp. Res. 1990, 22, 195–201. [Google Scholar] [CrossRef]
  86. Gough, L.; Shaver, G.R.; Carol, J.; Royer, D.L.; Laundre, J.A. Vascular plant species richness in Alaskan arctic tundra: The importance of soil pH. J. Ecol. 2000, 88, 54–66. [Google Scholar] [CrossRef]
  87. Schmidtlein, S.; Ewald, J. Landscape patterns of indicator plants for soil acidity in the Bavarian Alps. J. Biogeogr. 2003, 30, 1493–1503. [Google Scholar] [CrossRef]
  88. Darmody, R.G.; Thorn, C.E.; Schlyter, P.; Dixon, J.C. Relationship of vegetation distribution to soil propertiesin Kärkevagge, Swedish Lapland. Arct. Antarct. Alp. Res. 2004, 36, 21–32. [Google Scholar] [CrossRef]
  89. Thuiller, W.; Albert, C.; Araujo, M.B.; Berry, P.M.; Cabeza, M.; Guisan, A.; Hickler, T.; Midgley, G.F.; Paterson, J.; Schurr, F.M.; et al. Predicting global change impacts on plant species’ distributions: Future challenges. Perspect. Plant Ecol. 2008, 9, 137–152. [Google Scholar] [CrossRef]
  90. Chardon, N.I.; Wipf, S.; Rixen, C.; Beilstein, A.; Doak, D.F. Local trampling disturbance effects on alpine plant populations and communities: Negative implications for climate change vulnerability. Ecol. Evol. 2018, 8, 7921–7935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Tilman, D.; May, R.M.; Lehman, C.L.; Nowak, M.A. Habitat destruction and the extinction debt. Nature 1994, 371, 65–66. [Google Scholar] [CrossRef]
  92. Kuussaari, M.; Bommarco, R.; Heikkinen, R.K.; Helm, A.; Krauss, J.; Lindborg, R.; Öckinger, E.; Pärtel, M.; Pino, J.; Rodà, F.; et al. Extinction debt: A challenge for biodiversity conservation. Trends Ecol. Evol. 2009, 24, 564–571. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Detailed habitus of Carex rupestris in Velická dolina Valley, Vysoké Tatry Mountains (A); typical habitat in Hlinská dolina Valley, Vysoké Tatry Mountains (B); site threatened by trampling on Mt. Siná, Nízke Tatry Mountains (C); and overgrowing by competitively stronger dwarf pine (Pinus mugo) on Mt. Krakova hoľa, Nízke Tatry Mountains (D). Photo: J. Šibík.
Figure 1. Detailed habitus of Carex rupestris in Velická dolina Valley, Vysoké Tatry Mountains (A); typical habitat in Hlinská dolina Valley, Vysoké Tatry Mountains (B); site threatened by trampling on Mt. Siná, Nízke Tatry Mountains (C); and overgrowing by competitively stronger dwarf pine (Pinus mugo) on Mt. Krakova hoľa, Nízke Tatry Mountains (D). Photo: J. Šibík.
Diversity 11 00172 g001
Figure 2. Map of known occurrences of Carex rupestris in the Western Carpathians (Slovakia).
Figure 2. Map of known occurrences of Carex rupestris in the Western Carpathians (Slovakia).
Diversity 11 00172 g002
Figure 3. Dendrogram of the numerical classification of relevés with the occurrence of Carex rupestris in the Western Carpathians based on the ordinal method and Goodman–Kruskal’s coefficient. Abbreviations: AS—Astero-Seslerion (cluster A), CF—Caricion firmae (cluster B; the two sub-clusters reflect geographical differentiation, whereas cluster B2 includes typical relevés of the alliance from Belianske Tatry Mountains, and cluster B1 includes relevés from different parts of the Western Carpathians), FV—Festucion versicoloris (cluster C); CF-OE—transitional relevés between Caricion firmae and Oxytropido-Elynion (cluster D).
Figure 3. Dendrogram of the numerical classification of relevés with the occurrence of Carex rupestris in the Western Carpathians based on the ordinal method and Goodman–Kruskal’s coefficient. Abbreviations: AS—Astero-Seslerion (cluster A), CF—Caricion firmae (cluster B; the two sub-clusters reflect geographical differentiation, whereas cluster B2 includes typical relevés of the alliance from Belianske Tatry Mountains, and cluster B1 includes relevés from different parts of the Western Carpathians), FV—Festucion versicoloris (cluster C); CF-OE—transitional relevés between Caricion firmae and Oxytropido-Elynion (cluster D).
Diversity 11 00172 g003
Figure 4. CAP ordination plot based on distance matrix of Goodman–Kruskal’s coefficients showing differences in species composition among four clusters of relevés (A–B) with the occurrence of C. rupestris in the Western Carpathians. Scores of indicator species (IndVal > 0.65) are overlaid onto the ordination. Confidence ellipses (95%) are displayed to depict expected overlap of the clusters. Abbreviations of species names are as follows: Festtat Festuca tatrae; Joviglo Jovibarba globifera; Kernsax Kernera saxatilis; Seslalb Sesleria albicans; Minulan Minuartia langii; Thympul Thymus pulcherrimus; Campcoc Campanula cochleariifolia; Gentclu Gentiana clusii; Primaur Primula auricula; Trisalp Trisetum alpestre; Arenten Arenaria tenella; Saliret Salix reticulata; Festsup Festuca supina; Junctri Juncus trifidus; Bartalp Bartia alpina; Hupesel Huperzia selago, Saxiret Saxifraga retusa; Sileaca Silene acaulis; Careful Carex fuliginosa; Campalp Campanula alpina; Lloyser Lloydia serotina. For explanation of symbols see Figure 3.
Figure 4. CAP ordination plot based on distance matrix of Goodman–Kruskal’s coefficients showing differences in species composition among four clusters of relevés (A–B) with the occurrence of C. rupestris in the Western Carpathians. Scores of indicator species (IndVal > 0.65) are overlaid onto the ordination. Confidence ellipses (95%) are displayed to depict expected overlap of the clusters. Abbreviations of species names are as follows: Festtat Festuca tatrae; Joviglo Jovibarba globifera; Kernsax Kernera saxatilis; Seslalb Sesleria albicans; Minulan Minuartia langii; Thympul Thymus pulcherrimus; Campcoc Campanula cochleariifolia; Gentclu Gentiana clusii; Primaur Primula auricula; Trisalp Trisetum alpestre; Arenten Arenaria tenella; Saliret Salix reticulata; Festsup Festuca supina; Junctri Juncus trifidus; Bartalp Bartia alpina; Hupesel Huperzia selago, Saxiret Saxifraga retusa; Sileaca Silene acaulis; Careful Carex fuliginosa; Campalp Campanula alpina; Lloyser Lloydia serotina. For explanation of symbols see Figure 3.
Diversity 11 00172 g004
Figure 5. Conditional inference classification tree of vegetation clusters based on environmental data. Bonferroni adjusted probabilities (p), observed classification to clusters (%), and cross-validated sensitivity (SensitivityCV) are displayed for each cluster.
Figure 5. Conditional inference classification tree of vegetation clusters based on environmental data. Bonferroni adjusted probabilities (p), observed classification to clusters (%), and cross-validated sensitivity (SensitivityCV) are displayed for each cluster.
Diversity 11 00172 g005
Figure 6. Habitat suitability maps showing the spatial extent of relatively favorable environmental conditions for C. rupestris (predicted probability of occurrence > 0.5) under current climate and under two climate change scenarios (rcp2.6, rcp8.5) in the year 2050 and 2070. Elevation hillshade is shown in gray. The barplot displays relative proportion of suitable habitat (percentage of the total area of Slovakia) for the current climate and possible future conditions.
Figure 6. Habitat suitability maps showing the spatial extent of relatively favorable environmental conditions for C. rupestris (predicted probability of occurrence > 0.5) under current climate and under two climate change scenarios (rcp2.6, rcp8.5) in the year 2050 and 2070. Elevation hillshade is shown in gray. The barplot displays relative proportion of suitable habitat (percentage of the total area of Slovakia) for the current climate and possible future conditions.
Diversity 11 00172 g006
Table 1. The overview of vegetation and environmental data of 65 C. rupestris plots with known coordinates. The table shows the range of values (min–max), results of simple models testing for differences among clusters of vegetation plots (test statistics—F, or χ2 and associated probabilities—P), average values, and 95% confidence intervals (lower–upper limit) for each cluster. Results of post-hoc tests are displayed by the uppercase letters (clusters with the same letters are not significantly different).
Table 1. The overview of vegetation and environmental data of 65 C. rupestris plots with known coordinates. The table shows the range of values (min–max), results of simple models testing for differences among clusters of vegetation plots (test statistics—F, or χ2 and associated probabilities—P), average values, and 95% confidence intervals (lower–upper limit) for each cluster. Results of post-hoc tests are displayed by the uppercase letters (clusters with the same letters are not significantly different).
VariableUnitRangeF/χ2PCluster ACluster BCluster CCluster D
Vegetation data
 Nr. of relevés 27101711
 Species richness 14–5010.43<0.000123.5 (21.3–25.8) a27.1 (19.1–35.2) a36.1 (32.3–39.8) b29.6 (26.6–32.7) ab
Topology
 Altitudem1297–201289.24<0.00011455 (1425–1485) a1652 (1597–1708) b1862 (1805–1919) c1930 (1862–1997) c
 Slope°10–8014.13<0.000147 (43–51) a43 (34–51) a54 (48–60) a25 (18–33) b
 Aspect 0.04–3.470.150.93191.99 (1.77–2.22)1.39 (0.38–2.40)2.09 (1.66–2.52)1.22 (0.43–2.01)
Geological bedrock% lime/marl/mylonite85.63<0.0001100/0/0 a62/38/0 b0/6/94 c45/55/0 b
Climate
 Temperature°C0.35–4.07111.6<0.00013.0 (2.8–3.2) a2.1 (1.7–2.5) b0.5 (0.4–0.7) c1.3 (1.2–1.4) d
 Precipitationmm1283–2123124.4<0.00011431 (1403–1458) a1463 (1351–1574) a1938 (1899–1977) b1809 (1755–1863) c
Solar radiation
 Duration h·m−2·y−110–1405.590.0019120 (114–125) a96 (84–109) ab100 (79–121) ab83 (74–92) b
 Global radiationMWh·m−2·y−116–764.680.005467 (62–71) a57 (45–69) ab50 (40–60) b53 (46–61) ab
 Direct radiationMWh·m−2·y−10.2–553.750.015648 (44–52) a40 (30–50) ab35 (26–43) b36 (30–43) ab
 Diffused radiationMWh·m−2·y−10.01–19.516.85<0.000118.3 (17.7–18.9) a16.0 (14.6–17.5) b13.8 (11.7–15.8) b15.8 (14.8–16.7) b
 Reflected radiationMWh·m−2·y−10.03–16.933.330.02560.98 (0.80–1.16) ab1.37 (0.85–1.89) ab1.55 (0.36–3.45) a1.34 (1.11–1.58) b
Soil characteristics
 pHH204.30–7.7541.64<0.00017.39 (7.26–7.53) a7.11 (6.77–7.44) ac5.54 (5.21–5.86) b6.53 (5.95–7.11) c
 pHKCl4.03–7.6939.91<0.00017.28 (7.09–7.46) a7.11 (6.76–7.45) ac5.29 (4.93–5.66) b6.24 (5.86–6.62) c
 Cag·kg−10.6–12.49.3890.00016.2 (4.9–7.4) a6.4 (4.1–8.7) a2.3 (1.4–3.1) b6.8 (5.2–8.4) a
 Kmg·kg−147–3406.8870.0008189 (157–221) a145 (106–183) ab103 (76–130) b122 (67–176) ab
 Mgmg·kg−114–214845.44<0.00011210 (956–1463) a615 (58–1172) b48 (23–72) c144 (61–226) b
 Ctot%2.8–30.719.35<0.000122.1 (20.2–24.0) a20.2 (15.7–24.7) a9.3 (6.6–12.0) b18.4 (13.6–23.2) a
 N%0.23–2.658.6650.00021.56 (1.29–1.83) a1.37 (1.00–1.74) a0.70 (0.51–0.89) b1.39 (0.98–1.79) ab

Share and Cite

MDPI and ACS Style

Svitková, I.; Svitok, M.; Petrík, A.; Bernátová, D.; Senko, D.; Šibík, J. The Fate of Endangered Rock Sedge (Carex rupestris) in the Western Carpathians—The Future Perspective of an Arctic-Alpine Species under Climate Change. Diversity 2019, 11, 172. https://0-doi-org.brum.beds.ac.uk/10.3390/d11090172

AMA Style

Svitková I, Svitok M, Petrík A, Bernátová D, Senko D, Šibík J. The Fate of Endangered Rock Sedge (Carex rupestris) in the Western Carpathians—The Future Perspective of an Arctic-Alpine Species under Climate Change. Diversity. 2019; 11(9):172. https://0-doi-org.brum.beds.ac.uk/10.3390/d11090172

Chicago/Turabian Style

Svitková, Ivana, Marek Svitok, Anton Petrík, Dana Bernátová, Dušan Senko, and Jozef Šibík. 2019. "The Fate of Endangered Rock Sedge (Carex rupestris) in the Western Carpathians—The Future Perspective of an Arctic-Alpine Species under Climate Change" Diversity 11, no. 9: 172. https://0-doi-org.brum.beds.ac.uk/10.3390/d11090172

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