Next Article in Journal
Dynamics of a Predator–Prey Model with the Effect of Oscillation of Immigration of the Prey
Next Article in Special Issue
Diversity and Origin of the Central Mexican Alpine Flora
Previous Article in Journal
Characteristics of Eukaryotic Plankton Communities in the Cold Water Masses and Nearshore Waters of the South Yellow Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Altitudinal Vascular Plant Richness and Climate Change in the Alpine Zone of the Lefka Ori, Crete

1
Department of Geoinformation in Environmental Management, CIHEAM–Mediterranean Agronomic Institute of Chania, 73100 Chania, Greece
2
Department of Natural Resources, Faculty of Geo-Information Science and Earth Observation (ITC), University of Twente, 7500 AE Enschede, The Netherlands
3
School of Economics and Management, Open University of Cyprus, Latsia 2252, Nicosia, Cyprus
4
School of Pure and Applied Sciences, Open University of Cyprus, Latsia 2252, Nicosia, Cyprus
*
Author to whom correspondence should be addressed.
Submission received: 15 November 2020 / Revised: 5 January 2021 / Accepted: 7 January 2021 / Published: 9 January 2021
(This article belongs to the Special Issue Diversity, Ecology and Conservation of Alpine Plants)

Abstract

:
High mountain zones in the Mediterranean area are considered more vulnerable in comparison to lower altitudes zones. Lefka Ori massif, a global biodiversity hotspot on the island of Crete is part of the Global Observation Research Initiative in Alpine Environments (GLORIA) monitoring network. The paper examines species and vegetation changes with respect to climate and altitude over a seven-year period (2001–2008) at a range of spatial scales (10 m Summit Area Section-SAS, 5 m SAS, 1 m2) using the GLORIA protocol in a re-survey of four mountain summits (1664 m–2339 m). The absolute species loss between 2001–2008 was 4, among which were 2 endemics. At the scale of individual summits, the highest changes were recorded at the lower summits with absolute species loss 4 in both cases. Paired t-tests for the total species richness at 1 m2 between 2001–2008, showed no significant differences. No significant differences were found at the individual summit level neither at the 5 m SAS or the 10 m SAS. Time series analysis reveals that soil mean annual temperature is increasing at all summits. Linear regressions with the climatic variables show a positive effect on species richness at the 5 m and 10 m SAS as well as species changes at the 5 m SAS. In particular, June mean temperature has the highest predictive power for species changes at the 5 m SAS. Recorded changes in species richness point more towards fluctuations within a plant community’s normal range, although there seem to be more significant diversity changes in higher summits related to aspects. Our work provides additional evidence to assess the effects of climate change on plant diversity in Mediterranean mountains and particularly those of islands which remain understudied.

1. Introduction

Global warming represents one of the most significant threats to biodiversity, worldwide [1] with some studies [2] tentatively forecasting the loss of 1 million living species by 2040. During the last two decades, evidence of global warming has accumulated and has been documented in the periodical IPCC reports [3,4] with the most recent report [5] suggesting a warming of 0.85 (0.65 to 1.06) °C, over the past two centuries. In alpine areas, climate drives plant species distribution [6,7] and this is the main reason why they are considered sensitive to any future climate changes [6,7] and why alpine ecosystems may serve as case studies to detect early indications of such changes [8]. Empirical studies in the Alps [7,9,10,11] provided already evidence that many plants have extended their altitudinal range as a result of recent climate warming. Species richness in high alpine mountains has also revealed considerable alterations [7,11,12,13,14,15]. According to future projections for 2100, at low altitude mountain ranges, half of the total species might be at risk [16], with rare and endemic species being particularly vulnerable to climatic changes [17], a general concern also for subnival–nival species [7].
For the Mediterranean Basin, future warming is expected to exceed global rates by 25% [18] with regional climate models predicting a sharp decrease in precipitation coupled with increasing summer temperatures [19,20,21]. In the case of mountain areas, worldwide changes in temperature and precipitation are projected higher than other regions [4], with Southern European mountains being more vulnerable to climate change [22]. For Mediterranean mountains, projections of future climate, under different emission scenarios and Atmosphere-Ocean-Coupled General Circulation Models, predict warming over 1.4 °C and a reduction of precipitation over −4.8% until 2085 (which vary according to the employed scenarios) [22].
The overall picture of several global change scenarios [23] is that significant biodiversity loss is likely in the Mediterranean mountains, due to habitat tracking problems and interspecific competition. Higher temperatures in mountain regions will lead to an upward shift of biotic zones and a possible decrease in endemic species [3,24]. In Europe, there is already evidence for shifting of vegetation belts and species distribution ranges [9,25,26,27] which are likely to continue in the future [28,29]. European-wide mountain studies report an increase in plant species richness [30] over a 145-year period, with significant abundance changes for species in highest elevations [31]. In particular, a synthesis on plant diversity changes across European summits showed that on Mediterranean summits species, gains were outstripped by losses of cryophilic species, resulting in a net species loss of on average 1.4 species over a seven-year period [32].
Evidence on how climate change affects the flora of Mediterranean is relatively limited [33,34,35] with few field based monitoring studies [36,37]. The largest Mediterranean islands are in their majority mountainous and considered global biodiversity hotspots [38,39], for which there is a substantial body of climatic and ecological biological evidence on the impacts of climate change (see reviews in [40,41]). On those islands, mountain massifs are very important in the context of national and international biodiversity and the Lefka Ori in Crete is such an example [42,43].
The vegetation formations of the Lefka Ori above the tree line are typical of the oro-Mediterranean and alti-Mediterranean zone [44,45] of the Mediterranean mountains and comprise either low prickly scrub formations or communities of spiny cushion-shaped dwarf shrubs [45,46]. Recognizing its ecological importance, Lefka Ori has been included in the Global Observation Research Initiative in Alpine Environments (GLORIA) network (see www.gloria.ac.at). GLORIA is an international long-term observation network in alpine environments to detect impacts of climate change on mountain ecosystems. In 2001, the first field baseline data in GLORIA-Europe were collected from 18 target regions, amongst which is the Lefka Ori massif in Crete [35]. Variation among summits in species richness and measures of similarity have been recorded for GLORIA summits in Europe [11,34,35,47], in New Zealand [48] and in Australia [49]. Changes in species richness and composition, paralleling changes in temperature, have also been identified in many European mountain regions [7,11,15].
This paper is part of the GLORIA project methodological framework [8] which consists of re-visitation of summit-areas worldwide, including the Lefka Ori, every seven years. Therefore, the focus herein is on the vegetation and species patterns of four high summits on the massif. The importance of this work is twofold. This is the first consistent monitoring scheme of this kind in Greek mountains and the second in a Mediterranean island mountain. Moreover, the work provides additional evidence to the effects of climate change on plant biodiversity in Mediterranean mountains and particularly those of islands which remain understudied. The aim of this study is to detect species’ and vegetation change across four summits, which were first surveyed in 2001 [35] and re-surveyed in 2008 using the same GLORIA protocol. In order to achieve this aim, the paper addresses the following questions:
  • What are the changes in sub-surface temperature?
  • What are the changes in species richness, species, and vegetation turnover?
  • Are changes scale specific?
  • Do different life forms and endemics show strong changes over time?
  • What are the effects of temperature on species richness?
  • Are observed changes fluctuations or signals of climate change?

2. Materials and Methods

2.1. Study Area

The study area is the Lefka Ori massif situated on the western part of Crete, Greece (Figure 1), which culminates at the Pachnes summit (2453 m).
It is a rugged marble and dolomite massif rich in rock debris and karstic formations and the wettest place on the island with over 2000 mm mean annual precipitation [50]. Above the tree line, the vegetation is dominated by low prickly scrub formations or spiny cushion-shaped dwarf shrubs (such as Berberis cretica L., Euphorbia acanthothamnos Heldr. & Sart. ex Boiss., Juniperus oxycedrus L. subsp oxycedrus Acantholimon androcaceum (Jaub. & Spach) Boiss., and Astragalus angustifolius Lam) [51]. Grazing, by far the major human impact on Lefka Ori, decreases in intensity along altitude. The four selected summits were Low summit (1664 m), Chorafas (1965 m), South-East Kakovoli (2160 m), and Sternes (2339 m), reported hereafter as LOW, CHO, SEK, and STR, respectively.

2.2. Vegetation Sampling

In May–June 2008, we revisited four mountain summits of the Lefka Ori massif (GR-LEO study site) which were sampled in 2001. Summit selection and sampling were based on the GLORIA project multi-summit approach [8] (Figure S1). Two Summit Area Sections referred to herein as SAS were defined per summit within which floristic survey took place. The first SAS in each summit was defined as a polygon with four corners at each cardinal direction at 10 m lower from the summit top (referred to as 10 m SAS) where a complete list of plants was compiled. Within each 10 m SAS, a smaller SAS was defined (referred to as 5 m SAS) in a similar manner at 5 m lower from the summit top in order to evaluate the quantitative floristic composition for each principal exposure (N-E-S-W). On each summit, the four corner 1 m × 1 m quadrats of a 3 m × 3 m grid were established as permanent plots in each cardinal direction (i.e., east, south, west, north) 5 m below the highest summit point. In each of the 1 m × 1 m quadrats, we recorded all vascular plant species and the percent cover of each species was visually estimated following [8]. The position of every plot was recorded using a GPS receiver, accompanied by photo documentation.
For species identification, the Flora Europaea [52] and Mountain Flora of Greece [53] were used. Nomenclature of the plant taxa given in this paper is according to Turland, Chilton and Press [54]; Chilton and Turland [55]. For the nomenclature on chorotypes and life forms, we refer to Jahn and Schönfelder [56]; while bioclimatic zonation is according to Quézel [45].

2.3. Climatic Variation

In the center of each 3 m × 3 m grid, a data logger (StowAway Tidbit v2; Onset Corporation, Bourne, MA, USA) recorded the soil temperature 10 cm below the surface at hourly intervals continuously from June 2001 to July 2008. During this time period, temperature data were used to calculate annual values of absolute minimum soil temperature, annual daily mean soil temperature, absolute maximum soil temperature, temperature sums, growing degree days, and length of the growing season across the sampled years. The assumption made herein is that the defrosting period starts when the mean daily temperature is constantly above 0 °C while the freezing period starts in the autumn when the mean daily temperatures are constantly below 0 °C [35].

2.4. Data Analysis

2.4.1. Climate

We used the lm function in R [57] to build up a linear regression time series model for each data set. The explanatory variables were “Time t”, which refers to the number of days elapsed since the origin (t = 1, day one), the average temperature of each month, as a dummy variable and the cyclical effect. We modeled the daily changes on the average daily temperature (based on four measurements). In order to measure the monthly seasonal effect, we calculated the regression coefficient for each month, which could be interpreted as the average increase in the temperature due to specified month, which is constant over the years. All the monthly effects were statistically significant. The shape of the seasonal trend is taken into account with the cosine element, which is dependent on time t, i.e., how far away from the origin we are. In this way, we managed to explain about 95% of the variability of the daily average temperature.

2.4.2. Species Richness

We have compiled species richness data at the 1 m2 quadrats, on each aspect of the SAS at the 5 and 10 m contour line and took the means at the 5 m and 10 m SAS and at the ‘whole of summit’ species (5 + 10 m SAS) (Table 1 to provide a summit overview. We used paired t-tests, following a Shapiro-Wilk test for normality, to assess changes in total species richness at the 1 m2 quadrats and at the 5 m SAS and 10 m SAS for all summits, and linear regression models to assess temperature effects on species richness at various scales as well as changes in species richness. Prior to employing regression, we checked that all assumptions were met (including continuous variables, linear relationship of variables independence of observations, homoscedacity, outliers, and residuals). These statistical analyses were performed with PAST software version 4.03 [58].

2.4.3. Species Diversity and Turnover

We use Shannon diversity index (Shannon H’ Log Base 10) to measure diversity at every cardinal direction per summit as well as the overall diversity per summit. We used two indices to assess turnover at the summit level for all scales as follows [59]:
  • Species turnover based on individual species frequency:
    Tsp = (A + D)/(A + D + U),
    where A is the frequency of quadrats where the species appeared in 2008; D is the frequency of quadrats where the species disappeared; U is the frequency of quadrats where the species’ frequency was unchanged.
  • Vegetation turnover index at different spatial scales and life forms (at summit level):
    Tveg = (A + D)/(A + D + B)
    where A is the number of new species in 2008; D is the number of disappearing species; B is the number of species present in both years of comparison.
For both indices, low turnover values are close to 0.01 and complete turnover is 1.0.

3. Results

3.1. Climate

In all four models (Table 2), there is a constant increase on the average temperature per year. More specifically, it increases by 0.118 degrees per year (LOW), 0.035 degrees per year (CHO), 0.047 degrees per year (SEK), and 0.041 degrees per year (STR). The highest average temperature usually occurs in July and August. However, the average temperature between years remains constant without any significant change on its variability.

3.2. General Patterns

Overall, at the study area level (GR-LEO–4 summits), 66 species were recorded (Table S1), 18 of which were endemics. The absolute species loss was 4 (8 species were lost and 4 new gained) among which were 2 endemics (3 lost, 1 new gained) (Table 2). At the scale of individual summits, the highest changes were recorded at the lower summits (LOW, CHO) with absolute species loss 4 in both cases (Table 3).
Mean species richness at the 5 + 10 m SAS was highest at the lower summit (LOW) and lowest at the highest summit (demonstrate a trend in decreasing richness with altitude). This pattern was also observed at the 10 m SAS, 5 m SAS, and the 1 m2 level (Table 1).
A two tailed t-test comparing the total species richness at 1 m2 (64 quadrats) between 2001–2008, revealed no significant differences at the 0.05 level (t-value = 0.03541, p-value = 0.971807). No significant differences were found at the individual summit level neither at the 5 m SAS or the 10 m SAS between 2001 and 2008 (Table 4).

3.3. Species and Vegetation Turnover

Mean Tsp was higher at CHO summit (0.46 ± 0.12 SE) and lower at the SEK summit (0.18 ± 0.05 SE). Endemic species turnover was between 0.17 (at LOW 5 m SAS) to 0.33 (at STR 10 m SAS) (Figure 2 but also Table S2). At the site level, chamaephytes turnover was between 0 to 0.29 and hemicryptophytes turnover was between 0.16 to 0.44, with the exception of Euphorbia acanthothamnos Boiss., which at the site level had a turnover of 0.5 (Table S2). The other major cushion-shaped plants recorded such as Astragalus angustifolius (Willd.) Hayek subsp. angustifolius and Prunus prostrara Labill. had very low or no turnover at all. Phanerophytes recorded behave similarly, with Berberis cretica L. having the highest turnover 0.5; whereas Rhamnus lycioides subsp. oleoides (L.) Jahand. & Maire and Acer sempervirens L. show no changes at the site level. Vegetation turnover (Tveg) across the four summits, within the 5 m SAS, were between 0.11 and 0.32; whereas within the 10 m SAS (not including the 5 m SAS), Tveg values were between 0.12 and 0.31 (Table S3). Overall vegetation turnover ranged from 0.06 to 0.28. At the 1 m2 scale, turnover values range from no change to 0.25 at CHO summit.

3.4. Species Richness and Climate

Linear regressions with the climatic variables used can explain both species richness at the 5 and 10 m SAS as well as species changes at the 5 m SAS. They all have moderate to high predictive power (27 to 90%). Temperature has a positive effect on species richness at 5 SAS and 10 SAS with high predictive power (over 70%) and a positive effect on changes at 5 SAS but with moderate predictive power. June mean temperature, in particular, has the highest predictive power for species changes at the 5 m SAS. Linear regressions with the climatic variables used had no explanatory power for changes observed at the 10 m SAS (Table 4).

3.5. Species Diversity

Species diversity decreases with altitude for all summits and years with the exception of the highest summit (STR: see Figure 3). STR summit shows the highest variation in species diversity when different aspects and years are examined, while the lowest summit LOW shows the lowest. The highest species diversity is consistently observed in the East slopes for the 2001 data and in the South slopes for 2008 data.

4. Discussion

4.1. Species Richness and Climatic Parameters

The GLORIA initiative responds to the need for a consistent long-term monitoring, an activity timely and essential for detecting and predicting mountain species’ responses to local climate change [8]. The value of experimental field sites for climate change studies is promoted in conservation biology, particularly for tracing/documenting aberrant responses [60], and to document the effects of elevational gradients on species richness often contested in biogeography [61,62,63]. Compared to the changes reported by similar studies emerging from Europe [15,36,37], changes from Lefka Ori seem to be variable. The increase in species number on individual summit level was very small compared to what has been reported elsewhere (max three species). A comparison between and within summits for the two monitoring years (2001 and 2008) did not reveal any significant changes in species richness (Table 1 and Table 3) throughout sampling scales but evident changes in different cardinal directions were recorded (Figure 3). Temperature explains species richness at the spatial scales examined (5 SAS, 10 SAS) (Table 4), showing a positive effect. This is consistent with the other similar studies from GLORIA network [59] alluding to an increase of warm-adapted species as temperature increases continues [64]. Temperature also explains changes at the 5 m SAS and in particular June mean temperature has the highest predictive power (Table 4). At the coarser spatial scale (10 m SAS), temperature had no effect on recorded changes. This is possibly due to the highest contribution/presence of larger perennial shrubs at that scale suggesting that changes are perhaps related more to biotic interactions than temperature changes [65].
In this study, the sub-surface soil temperature data over a seven-year period for the GR-LEO study site showed that there is a constant increase in the average temperature per year ranging from 0.047 (SEK) to 0.118 degree (LOW) (Table 2). Based on observed data for the periods 1970–2000 in Crete, Tsanis et al. [66] give a monthly temperature increase of 0.5 degree per decade and a monthly average precipitation decrease of 26 mm per decade. Recent research for Crete suggests that the known trends for the eastern Mediterranean will be mirrored on the island’s climate with the optimistic estimates of 12% lower temperature and 2 degrees’ temperature increase until 2040 further aggravating, following the same trend in the decades to come [66].

4.2. Changes in Species Composition and Abundance

Overall, there are less pronounced species changes at the smallest sampling scales while it seems that changes are more evident in middle altitudes, particularly the increase of therophytes and the highest turnover at the CHO summit (Figure 2). Similar trends were also reported during the first survey of the area [35], while high species turnover in middle altitudes is also suggested by modeling studies for Mediterranean mountains [25].
There is a higher turnover in chamaephytes–hemicryptophytes compared to other life forms across summits and sampling scales. Increase in chamaephytes has also been reported for other Mediterranean mountains [36,37]. Since chamaephytes and hemicryptophytes constitute the majority of the flora in the study area [51], this result is expected. These changes of course may indicate changes in biotic relations between life forms i.e., a competitive advantage of some life forms at the expense of others. Although not a pronounced change, a moderate increase was recorded in therophytes in the middle summit CHO, which has also the highest species turnover (Figure 2).
As evidence on the responses of mountain ecosystems on climate change are mounting, it becomes clearer that there are many aspects which need to be taken into consideration in order to draw reliable conclusions (see review in [67]). For example, plants’ response will be too slow to keep pace with the speed of climate change [68] or that invasions of non-native plants on mountainous areas cannot be excluded [69]. In fact, a summit’s response will be ‘individualist’ since vegetation composition particularly between the lowest and the highest summits in studies similar [36,37] to this one usually differs. Therefore, the traits of the constituent species in every summit will determine their ability to compete, disperse, and adapt. Alexander et al. [70] also highlight the importance of community dynamics in determining mountain species responses. This can be either related to new composition but also to facilitation of existing plants. For example, the role of cushion plants in facilitating range expansion has been demonstrated for mountain environments [65,71]. Mediterranean mountains in general and Lefka Ori in particular have a high number of cushion plants such as Astragalus angustifolius Lam., Acantholimon androcaceum (Jaub. & Spach) Boiss., but also shrubs such as Berberis cretica L. which act as ‘eco-engineers’ in a harsh environment. Their persistence and role in the future may be critical for the survival of other mountain plants. This study did not record any significant changes of the major cushion-shaped plants across summits.
Although species diversity decreases with altitude and across years for three out of four summits, there is an increase in overall species diversity from 2001–2008 in the highest summit (STR) (Table 3). In addition, STR summit shows the highest variation in species diversity when different aspects between the two years were surveyed. Differences in species numbers among cardinal directions in European mountain summits have been recently documented [72]. This study suggests that species richness on Mediterranean summits tend to be highest on the West facing and lowest on the East and North facing summit slopes, although no significant differences were recorded. What is emerging from the GR-LEO summit data is that species diversity is higher in the East slopes for the 2001 data and the South slopes for the 2008 data, compared to the other two directions. This might point towards a more individualistic response per site and the effect of summer drought, which may limit species numbers on the warmer slopes of Mediterranean mountains.

4.3. Changes in Endemic Species

Mountain areas with high endemism level have received much attention in the literature owing to the vulnerability of endemic species [8,73,74]. The flora of the Cretan mountains comprises relict species, many of which are endemic derivatives of lowland species and species that also occur in the continental Greek mountains such as Scutellaria hirta Boiss. from Scutellaria sieberi Benth and Bellis longifolia Boiss. & Heldr. from Bellis sylvestris Cirillo [75]. According to [42], regional and local endemism in Greece increase in a southerly direction, culminating in the Lefka Ori, which displays one of the highest rates of narrow endemism in the Mediterranean area. Narrow endemics seem to be highly threatened on the low mountain summits of the north-eastern Alps [17,73] and this was also reported for the Lefka Ori in Crete [35], compared to endemics with wider distributions [76]. Our results show that endemics’ turnover was very small across summits and overall in the area where a small decrease was recorded in absolute numbers. Still the highest turnover was recorded at the middle altitude (CHO) but also at the highest summit (STR). The first finding is consistent with the overall turnover vegetation pattern recorded herein. In the case of STR, a high value is due to the fact that contribution of the endemics to the overall flora of the Lefka Ori increases with altitude (i.e., endemics richness along an elevation gradient) [51]. In similar fashion, a small decrease in endemics is reported also for the Apennines [36] while in Sierra Nevada it seems that endemics were more affected [77]. These findings might support the view of the mountains’ resilience at least as endemic flora is concerned although findings should be interpreted with caution [78].
Mediterranean mountains acted as refugia in the Mediterranean Basin during the Pleistocene cycles preventing species extinctions and likely supporting speciation processes [22,79]. Medail and Diadema [80] identified 18 out of 52 of these refugia in mountains, among which are the Pyrenees, South Apennines, Middle and High Atlas, Taurus and the mountains of Crete. In addition to differences in biogeography, glaciation pattern seems to be very different in the Apennines and the Sierra Nevada compared to Crete with evidence still contested for the Lefka Ori (see [81]).
Lefka Ori still host important populations of tertiary elements [51]. As in the case of several other Mediterranean mountains, some of these elements are in isolated occurrences i.e., linked to specific azonal habitats resulting from geomorphological and hydrological processes [82]. Micro-habitat diversity influences species distribution at small scales [83] with some micro-habitats reported to have lower than soil surface temperatures [84] acting as a potential buffer to climate change [83].
Even if endemics persist, there might be other consequences for them including changes in their phenology, reported already in other mountain regions of the world [85], and disproportionate range losses [23,73]. Range boundaries of many narrow endemics in mountain environments are associated with localized environments (geology, geomorphology). Since these species already exhibit realized climatic niches much more restricted compared to their potential climatic niche, their distributions might not be immediately affected by climate change [86] nor influenced by other species’ altitudinal shift due to biotic interactions [35]. Out of the 18 endemics recorded in this study, the ones with the highest turnover rate included Crepis sibthorpiana Boiss. & Heldr., Centaurea idea Boiss. & Heldr. and Alyssum fragillimum (Bald.) Rech. f. These three species are locally quite common i.e., widely distributed along the massif over 1500 m and are not restricted to a single specialized habitat. Although turnover might be due to interannual fluctuations, it might imply that widely distributed species will be more affected by changes in climate, compared to narrow endemics associated with localized habitats. In addition, a recent study in Sierra Nevada showed that interspecific hybridization between narrowly endemic alpine plants and widely distributed lowland plants is an overlooked consequence of altitudinal shift of the latter [87]. Therefore, studying endemic plants population viabilities and extinction probabilities under interspecific competition have been suggested [11].

5. Conclusions

Recorded changes in species richness in Lefka Ori point more towards fluctuations within a plant community’s normal range, although there seem to be more significant diversity changes in higher summits related to aspects. Empirical studies like this one may serve as a validation and refinement of climate species modeling related techniques in an attempt to draw meaningful conclusion on the effects of climate change on biodiversity [88]. For example, a recent synthesis [89] predicted small extinction risks in Europe and the Mediterranean area compared to other geographical areas of the world. Subsequent re-surveys will determine whether the trends and patterns recorded herein will hold. In addition, this study did not find any evidence showing that endemics have been influenced by changes in climate parameters, although the pace of current global warming and its interactions with other drivers of change has no precedent [90]. Mediterranean mountains might continue to hold as biodiversity refugia for the years to come under the current climate change projections.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-2818/13/1/22/s1, Figure S1: Scheme of the Multi-Summit sampling design (from Pauli, H.; Gottfried, M.; Lamprecht, A.; Niessner, S.; Rumpf, S.; Winkler, M.; Steinbauer, K. and Grabherr, G., coordinating authors and editors (2015). The GLORIA field manual—standard Multi-Summit approach, supplementary methods and extra approaches. 5th edition. GLORIA-Coordination, Austrian Academy of Sciences & University of Natural Resources and Life Sciences, Vienna). Table S1: List of species present across the GR-LEO summits (LOW, CHO, SEK and STR) during the 2001 and 2008 surveys, their families, life forms (P = Phanerophyte, Ch = Chamaephyte, H = Hemicryptophyte, G = Geophyte, T = Therophyte) and their endemism to Cretan flora. Table S2: Species turnover (Tsp) for species recorded in the 1 m2 quadrats for each summit (LOW, CHO, SEK AND STR) and for the four summits combined (GR-LEO); A = Number of quadrats the species appeared in 2008, D = Number of quadrats the species disappeared in 2008, U = Number of quadrats the species unchanged. Table S3: The species richness in 2001 and 2008, the number of new and lost species and also the number of species found in both surveys; used for calculating the vegetation turnover (Tveg) between 2001 and 2008 in summit area sections and for different groups of species and life forms.

Author Contributions

Conceptualization, G.K. and D.G.; methodology, G.K., D.G., P.N. and I.N.V.; fieldwork, G.K., D.G., I.R., and P.N., formal analysis, I.N.V. and M.A.T.; writing—original draft preparation, I.N.V.; writing—review and editing, All authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the project EU-GLORIA (2001–2003) No. EVK2-CT-2000-00056. Funding for the second monitoring cycle came from the Austrian Academy of Sciences.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available in the GLORIA network repository (gloria.ac.at).

Acknowledgments

We would like to acknowledge the support of the GLORIA office at BOKU.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Malcolm, J.R.; Liu, C.; Neilson, R.P.; Hansen, L.; Hannah, L. Global Warming and Extinctions of Endemic Species from Biodiversity Hotspots. Conserv. Biol. 2006, 20, 538–548. [Google Scholar] [CrossRef]
  2. Thomas, C.D.; Cameron, A.; Green, R.E.; Bakkenes, M.; Beaumont, L.J.; Collingham, Y.C.; Erasmus, B.F.N.; de Siqueira, M.F.; Grainger, A.; Hannah, L.; et al. Extinction Risk from Climate Change. Nature 2004, 427, 145–148. [Google Scholar] [CrossRef] [PubMed]
  3. McCarthy, J.J.; Intergovernmental Panel on Climate Change (Eds.) Climate Change 2001: Impacts, Adaptation, and Vulnerability: Contribution of Working Group II to the Third Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2001; ISBN 978-0-521-80768-5. [Google Scholar]
  4. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Tignor, M.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, V.; Midgley, P.M. (Eds.) IPCC, 2013: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 1535. [Google Scholar]
  5. Hock, R.; Rasul, G.; Adler, C.; Cáceres, B.; Gruber, S.; Hirabayashi, Y.; Jackson, M.; Kääb, A.; Kang, S.; Kutuzov, S.; et al. High Mountain Areas; IPCC—Intergovernmental Panel on Climate Change: Geneva, Switzerland, 2019; pp. 131–202. [Google Scholar]
  6. Körner, C. Alpine Plant. Life: Functional Plant. Ecology of High. Mountain Ecosystems, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2003; ISBN 978-3-540-00347-2. [Google Scholar]
  7. Pauli, H.; Gottfried, M.; Reiter, K.; Klettner, C.; Grabherr, G. Signals of Range Expansions and Contractions of Vascular Plants in the High Alps: Observations (1994–2004) at the GLORIA Master Site Schrankogel, Tyrol, Austria. Glob. Chang. Biol. 2007, 13, 147–156. [Google Scholar] [CrossRef]
  8. Pauli, H.; European Commission; Directorate General for Research. The Gloria Field Manual: Multi-Summit Approach; Directorate-General for Research: Luxembourg, 2004; ISBN 978-92-894-4737-9. [Google Scholar]
  9. Holzinger, B.; Hülber, K.; Camenisch, M.; Grabherr, G. Changes in Plant Species Richness over the Last Century in the Eastern Swiss Alps: Elevational Gradient, Bedrock Effects and Migration Rates. Plant Ecol. 2008, 195, 179–196. [Google Scholar] [CrossRef]
  10. Vittoz, P.; Bodin, J.; Ungricht, S.; Burga, C.A.; Walther, G.-R. One Century of Vegetation Change on Isla Persa, a Nunatak in the Bernina Massif in the Swiss Alps. J. Veg. Sci. 2008, 19, 671–680. [Google Scholar] [CrossRef]
  11. Erschbamer, B.; Kiebacher, T.; Mallaun, M.; Unterluggauer, P. Short-Term Signals of Climate Change along an Altitudinal Gradient in the South Alps. Plant Ecol. 2009, 202, 79–89. [Google Scholar] [CrossRef]
  12. Pauli, H.; Gottfried, M.; Reiter, K.; Grabherr, G. High Mountain Summits as Sensitive Indicators of Climate Change Effects on Vegetation Patterns: The “Multi Summit-Approach” of GLORIA (Global Observation Research Initiative in Alpine Environments). In Global Change and Protected Areas; Visconti, G., Beniston, M., Iannorelli, E.D., Barba, D., Eds.; Advances in Global Change Research; Springer: Dordrecht, The Netherlands, 2001; Volume 9, pp. 45–51. ISBN 978-90-481-5686-3. [Google Scholar]
  13. Virtanen, R.; Dirnböck, T.; Dullinger, S.; Grabherr, G.; Pauli, H.; Staudinger, M.; Villar, L. Patterns in the Plant Species Richness of European High Mountain Vegetation. In Alpine Biodiversity in Europe; Nagy, L., Grabherr, G., Körner, C., Thompson, D.B.A., Eds.; Ecological Studies; Springer: Berlin/Heidelberg, Germany, 2003; Volume 167, pp. 149–172. ISBN 978-3-642-62387-5. [Google Scholar]
  14. Parolo, G.; Rossi, G. Upward Migration of Vascular Plants Following a Climate Warming Trend in the Alps. Basic Appl. Ecol. 2008, 9, 100–107. [Google Scholar] [CrossRef]
  15. Michelsen, O.; Syverhuset, A.O.; Pedersen, B.; Holten, J.I. The Impact of Climate Change on Recent Vegetation Changes on Dovrefjell, Norway. Diversity 2011, 3, 91–111. [Google Scholar] [CrossRef] [Green Version]
  16. Dullinger, S.; Gattringer, A.; Thuiller, W.; Moser, D.E.; Zimmermann, N.; Guisan, A.; Willner, W.; Plutzar, C.; Leitner, M.; Mang, T.; et al. Extinction debt of high-mountain plants under twenty-first-century climate change. Nat. Clim. Chang. 2012, 2, 619–622. [Google Scholar] [CrossRef]
  17. Grabherr, G.; Gottfried, M.; Gruber, A.; Pauli, H. Patterns and Current Changes in Alpine Plant Diversity. Ecol. Stud. Anal. Synth. (USA) 1995, 113, 167–181. [Google Scholar]
  18. Cramer, W.; Guiot, J.; Fader, M.; Garrabou, J.; Gattuso, J.-P.; Iglesias, A.; Lange, M.A.; Lionello, P.; Llasat, M.C.; Paz, S.; et al. Climate Change and Interconnected Risks to Sustainable Development in the Mediterranean. Nat. Clim. Chang. 2018, 8, 972–980. [Google Scholar] [CrossRef] [Green Version]
  19. Fronzek, S.; Carter, T.R.; Jylhä, K. Representing Two Centuries of Past and Future Climate for Assessing Risks to Biodiversity in Europe. Glob. Ecol. Biogeogr. 2012, 21, 19–35. [Google Scholar] [CrossRef]
  20. Kirtman, B.; Power, S.B.; Adedoyin, A.J.; Boer, G.J.; Bojariu, R.; Doblas-Reyes, F.; Fiore, A.M.; Kimoto, M.; Meehl, G.; Prather, M.; et al. Near-Term Climate Change: Projections and Predictability. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 76. [Google Scholar]
  21. Polade, S.D.; Pierce, D.W.; Cayan, D.R.; Gershunov, A.; Dettinger, M.D. The Key Role of Dry Days in Changing Regional Climate and Precipitation Regimes. Sci. Rep. 2015, 4, 4364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Nogués-Bravo, D.; Araújo, M.B.; Lasanta, T.; Moreno, J.I.L. Climate Change in Mediterranean Mountains during the 21st Century. Ambi 2008, 37, 280–285. [Google Scholar] [CrossRef]
  23. Thuiller, W.; Lavorel, S.; Araujo, M.B.; Sykes, M.T.; Prentice, I.C. Climate Change Threats to Plant Diversity in Europe. Proc. Natl. Acad. Sci. USA 2005, 102, 8245–8250. [Google Scholar] [CrossRef] [Green Version]
  24. Mooney, H.A.; Kalin Arroyo, M.T.; Bond, W.J.; Canadell, J.; Hobbs, R.J.; Lavorel, S.; Neilson, R.P. Mediterranean-Climate Ecosystems. In Global Biodiversity in a Changing Environment: Scenarios for the 21st Century; Chapin, F.S., Sala, O.E., Huber-Sannwald, E., Eds.; Ecological Studies; Springer: New York, NY, USA, 2001; pp. 157–199. ISBN 978-1-4613-0157-8. [Google Scholar]
  25. Ruiz-Labourdette, D.; Nogués-Bravo, D.; Ollero, H.S.; Schmitz, M.F.; Pineda, F.D. Forest Composition in Mediterranean Mountains Is Projected to Shift along the Entire Elevational Gradient under Climate Change: Forest Dynamics under Climate Change. J. Biogeogr. 2012, 39, 162–176. [Google Scholar] [CrossRef]
  26. Kapralov, D.S.; Shiyatov, S.G.; Moiseev, P.A.; Fomin, V.V. Changes in the Composition, Structure, and Altitudinal Distribution of Low Forests at the Upper Limit of Their Growth in the Northern Ural Mountains. Russ. J. Ecol. 2006, 37, 367–372. [Google Scholar] [CrossRef]
  27. Kullman, L. Tree Line Population Monitoring of Pinus Sylvestris in the Swedish Scandes, 1973?2005: Implications for Tree Line Theory and Climate Change Ecology. J. Ecol. 2007, 95, 41–52. [Google Scholar] [CrossRef]
  28. Lenoir, J.; Gegout, J.C.; Marquet, P.A.; de Ruffray, P.; Brisse, H. A Significant Upward Shift in Plant Species Optimum Elevation During the 20th Century. Science 2008, 320, 1768–1771. [Google Scholar] [CrossRef]
  29. Engler, R.; Randin, C.F.; Thuiller, W.; Dullinger, S.; Zimmermann, N.E.; Araújo, M.B.; Pearman, P.B.; Lay, G.L.; Piedallu, C.; Albert, C.H.; et al. 21st Century Climate Change Threatens Mountain Flora Unequally across Europe. Glob. Chang. Biol. 2011, 17, 2330–2341. [Google Scholar] [CrossRef]
  30. Steinbauer, M.J.; Grytnes, J.-A.; Jurasinski, G.; Kulonen, A.; Lenoir, J.; Pauli, H.; Rixen, C.; Winkler, M.; Bardy-Durchhalter, M.; Barni, E.; et al. Accelerated Increase in Plant Species Richness on Mountain Summits Is Linked to Warming. Nature 2018, 556, 231–234. [Google Scholar] [CrossRef] [PubMed]
  31. Rumpf, S.B.; Hülber, K.; Klonner, G.; Moser, D.; Schütz, M.; Wessely, J.; Willner, W.; Zimmermann, N.E.; Dullinger, S. Range Dynamics of Mountain Plants Decrease with Elevation. Proc. Natl. Acad. Sci. USA 2018, 115, 1848–1853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Pauli, H.; Gottfried, M.; Dullinger, S.; Abdaladze, O.; Akhalkatsi, M.; Alonso, J.L.B.; Coldea, G.; Dick, J.; Erschbamer, B.; Calzado, R.F.; et al. Recent Plant Diversity Changes on Europe’s Mountain Summits. Science 2012, 336, 353–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Peñuelas, J.; Boada, M. A Global Change-Induced Biome Shift in the Montseny Mountains (NE Spain). Glob. Chang. Biol. 2003, 9, 131–140. [Google Scholar] [CrossRef] [Green Version]
  34. Stanisci, A.; Pelino, G.; Blasi, C. Vascular Plant Diversity and Climate Change in the Alpine Belt of the Central Apennines (Italy). Biodivers. Conserv. 2005, 14, 1301–1318. [Google Scholar] [CrossRef]
  35. Kazakis, G.; Ghosn, D.; Vogiatzakis, I.N.; Papanastasis, V.P. Vascular Plant Diversity and Climate Change in the Alpine Zone of the Lefka Ori, Crete. Biodivers. Conserv. 2007, 16, 1603–1615. [Google Scholar] [CrossRef]
  36. Stanisci, A.; Frate, L.; Morra Di Cella, U.; Pelino, G.; Petey, M.; Siniscalco, C.; Carranza, M.L. Short-Term Signals of Climate Change in Italian Summit Vegetation: Observations at Two GLORIA Sites. Plant Biosyst. An. Int. J. Deal. All Asp. Plant Biol. 2016, 150, 227–235. [Google Scholar] [CrossRef] [Green Version]
  37. Calzado, M.R.F.; Mesa, J.M.; Merzouki, A.; Porcel, M.C. Vascular Plant Diversity and Climate Change in the Upper Zone of Sierra Nevada, Spain. Plant Biosyst. An. Int. J. Deal. All Asp. Plant Biol. 2012, 146, 1044–1053. [Google Scholar] [CrossRef]
  38. Medail, F.; Quezel, P. Biodiversity Hotspots in the Mediterranean Basin: Setting Global Conservation Priorities. Conserv. Biol. 1999, 13, 1510–1513. [Google Scholar] [CrossRef]
  39. Davis, S.D.; Heywood, V.H.; Hamilton, A.C.; World Wide Fund for Nature; International Union for Conservation of Nature and Natural Resources (Eds.) Centres of Plant Diversity: A Guide and Strategy for Their Conservation; World Wide Fund for Nature (WWF) and ICUN-World Conservation Union: Cambridge, UK, 1994; ISBN 978-2-8317-0197-4. [Google Scholar]
  40. Vogiatzakis, I.N.; Mannion, A.M.; Sarris, D. Mediterranean Island Biodiversity and Climate Change: The Last 10,000 Years and the Future. Biodivers. Conserv. 2016, 25, 2597–2627. [Google Scholar] [CrossRef]
  41. Médail, F. The Specific Vulnerability of Plant Biodiversity and Vegetation on Mediterranean Islands in the Face of Global Change. Reg. Environ. Chang. 2017, 17, 1775–1790. [Google Scholar] [CrossRef] [Green Version]
  42. Strid, A. Phytogeographical Aspects of the Greek Mountain Flora. Fragm. Florist. Et Geobot. Suppl. 1993, 2, 411–433. [Google Scholar]
  43. Phitos, D. (Ed.) The Red Data Book of Rare and Threatened Plants of Greece; 1. Publ.; World Wide Fund for Nature: Athen, Greece, 1995; ISBN 978-960-7506-04-7. [Google Scholar]
  44. Quézel, P. The study of groupings in the countries surrounding the Mediterranean: Some methodological aspects. In Mediterranean-Type Shrublands; Di Castri, F., Goodall, D.W., Sprecht, R.L., Eds.; Elsevier: Amsterdam, The Netherlands, 1981; pp. 87–93. [Google Scholar]
  45. Quézel, P. Floristic composition and phytosociological structure of sclerophyllous mattoral around the Mediterranean. In Mediterranean-Type Shrublands; Di Castri, F., Goodall, D.W., Sprecht, R.L., Eds.; Elsevier: Amsterdam, The Netherlands, 1981; pp. 107–121. [Google Scholar]
  46. Zafran, J. Contributions à La Flore et à Vegétation de La Crète; Université de Provence: Aix en Provence, France, 1990. [Google Scholar]
  47. Coldea, G.; Pop, A. Floristic Diversity in Relation to Geomorphological and Climatic Factors in the Subalpinealpine Belt of the Rodna Mountains (the Romanian Carpathians). Pirineos 2004, 158–159, 61–72. [Google Scholar] [CrossRef] [Green Version]
  48. Mark, A.F.; Dickinson, K.J.M.; Maegli, T.; Halloy, S.R.P. Two GLORIA Long-term Alpine Monitoring Sites Established in New Zealand as Part of a Global Network. J. R. Soc. N. Z. 2006, 36, 111–128. [Google Scholar] [CrossRef]
  49. Pickering, C.; Hill, W.; Green, K. Vascular Plant Diversity and Climate Change in the Alpine Zone of the Snowy Mountains, Australia. Biodivers. Conserv. 2008, 17, 1627–1644. [Google Scholar] [CrossRef] [Green Version]
  50. Rackham, O.; Moody, J.A. The Making of the Cretan Landscape; Manchester University Press, Distributed Exclusively in the USA and Canada by St. Martin’s Press: Manchester, UK; New York, NY, USA, 1996; ISBN 978-0-7190-3646-0. [Google Scholar]
  51. Vogiatzakis, I.N.; Griffiths, G.H.; Mannion, A.M. Environmental Factors and Vegetation Composition, Lefka Ori Massif, Crete, S. Aegean. Environ. Factors Veg. Compos. Crete. Glob Ecol. Biogeogr. 2003, 12, 131–146. [Google Scholar] [CrossRef]
  52. Tutin, T.G.; Heywood, V.H.; Burges, N.A.; Moore, D.M.; Valentine, D.H.; Walters, S.M.; Webb, D.A. Flora Europaea Volumes 1–5; Cambridge University Press: Cambridge, UK; New York, NY, USA, 1964. [Google Scholar]
  53. Strid, A.; Tan, K. (Eds.) Mountain Flora of Greece; Cambridge University Press: Cambridge, UK; New York, NY, USA, 1986; ISBN 978-0-521-25737-4. [Google Scholar]
  54. Turland, N.J.; Chilton, J.R.; Press, J.R. Flora of the Cretan Area: Annotated Checklist and Atlas; HMSO Publications: London, UK, 1993. [Google Scholar]
  55. Turland, N.; Chilton, L. Flora of Crete: Supplement II, Additions 1997–2008; Marengo Publications: Retford, UK, 2008. [Google Scholar]
  56. Jahn, R.; Schönfelder, P. Exkursionsflora Für Kreta; Ulmer: Stuttgart, Germany, 1995. [Google Scholar]
  57. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria; Available online: https://www.r-project.org/ (accessed on 21 December 2020).
  58. Hammer, O.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontol. Electron. 2001, 4, 9. [Google Scholar]
  59. Venn, S.; Pickering, C.; Green, K. Short-Term Variation in Species Richness across an Altitudinal Gradient of Alpine Summits. Biodivers. Conserv. 2012, 21, 3157–3186. [Google Scholar] [CrossRef]
  60. Parmesan, C.; Hanley, M.E. Plants and Climate Change: Complexities and Surprises. Ann. Bot. 2015, 116, 849–864. [Google Scholar] [CrossRef]
  61. Rahbek, C. The Elevational Gradient of Species Richness: A Uniform Pattern? Ecography 1995, 18, 200–205. [Google Scholar] [CrossRef]
  62. Lomolino, M.V. Elevation Gradients of Species-Density: Historical and Prospective Views. Glob. Ecol. Biogeogr. 2001, 10, 3–13. [Google Scholar] [CrossRef]
  63. Korner, C.; Spehn, E.M.; Spehn, E.M. Mountain Biodiversity: A Global Assessment; Routledge: Abingdon, UK, 2019; ISBN 978-0-429-34258-5. [Google Scholar]
  64. Gottfried, M.; Pauli, H.; Futschik, A.; Akhalkatsi, M.; Barančok, P.; Benito Alonso, J.L.; Coldea, G.; Dick, J.; Erschbamer, B.; Fernández Calzado, M.R.; et al. Continent-Wide Response of Mountain Vegetation to Climate Change. Nat. Clim. Chang. 2012, 2, 111–115. [Google Scholar] [CrossRef]
  65. Cavieres, L.A.; Brooker, R.W.; Butterfield, B.J.; Cook, B.J.; Kikvidze, Z.; Lortie, C.J.; Michalet, R.; Pugnaire, F.I.; Schöb, C.; Xiao, S.; et al. Facilitative Plant Interactions and Climate Simultaneously Drive Alpine Plant Diversity. Ecol. Lett. 2014, 17, 193–202. [Google Scholar] [CrossRef] [PubMed]
  66. Tsanis, I.K.; Koutroulis, A.G.; Daliakopoulos, I.N.; Jacob, D. Severe Climate-Induced Water Shortage and Extremes in Crete. Clim. Chang. 2011, 106, 667–677. [Google Scholar] [CrossRef]
  67. Malanson, G.P.; Resler, L.M.; Butler, D.R.; Fagre, D.B. Mountain Plant Communities: Uncertain Sentinels? Prog. Phys. Geogr. Earth Environ. 2019, 43, 521–543. [Google Scholar] [CrossRef]
  68. Essl, F.; Dullinger, S.; Rabitsch, W.; Hulme, P.E.; Pyšek, P.; Wilson, J.R.U.; Richardson, D.M. Delayed Biodiversity Change: No Time to Waste. Trends Ecol. Evol. 2015, 30, 375–378. [Google Scholar] [CrossRef]
  69. Petitpierre, B.; McDougall, K.; Seipel, T.; Broennimann, O.; Guisan, A.; Kueffer, C. Will Climate Change Increase the Risk of Plant Invasions into Mountains? Ecol. Appl. 2016, 26, 530–544. [Google Scholar] [CrossRef]
  70. Alexander, J.M.; Diez, J.M.; Levine, J.M. Novel Competitors Shape Species’ Responses to Climate Change. Nature 2015, 525, 515–518. [Google Scholar] [CrossRef]
  71. le Roux, P.C.; McGeoch, M.A. Spatial Variation in Plant Interactions across a Severity Gradient in the Sub-Antarctic. Oecologia 2008, 155, 831–844. [Google Scholar] [CrossRef]
  72. Winkler, M.; Lamprecht, A.; Steinbauer, K.; Hülber, K.; Theurillat, J.-P.; Breiner, F.; Choler, P.; Ertl, S.; Gutiérrez Girón, A.; Rossi, G.; et al. The Rich Sides of Mountain Summits—A Pan-European View on Aspect Preferences of Alpine Plants. J. Biogeogr. 2016, 43, 2261–2273. [Google Scholar] [CrossRef]
  73. Dirnböck, T.; Essl, F.; Rabitsch, W. Disproportional Risk for Habitat Loss of High-Altitude Endemic Species under Climate Change. Glob. Chang. Biol. 2011, 17, 990–996. [Google Scholar] [CrossRef]
  74. Grabherr, G.; Gottfried, M.; Pauli, H. GLORIA: A Global Observation Research Initiative in Alpine Environments. Mt. Res. Dev. 2000, 20, 190–191. [Google Scholar] [CrossRef]
  75. Greuter, W. The relict element of the flora of Crete and its evolutionary significance. In Taxonomy, Phytogeography and Evolution; Academic Press: London, UK, 1972. [Google Scholar]
  76. Theurillat, J. Climate change and the alpine flora: Some perspectives. In Potential Ecological Impacts of Climate Change in the Alps and Fennoscandian Mountains; Jardin Botanique: Geneve, Switzerland, 1995; pp. 121–127. [Google Scholar]
  77. Fernández Calzado, M.R.; Molero, J. Changes in the Summit Flora of a Mediterranean Mountain (Sierra Nevada, Spain) as a Possible Effect of Climate Change. LAZA 2013, 34, 65–75. [Google Scholar] [CrossRef] [Green Version]
  78. Harrison, S.; Noss, R. Endemism Hotspots Are Linked to Stable Climatic Refugia. Ann. Bot. 2017, 119, 207–214. [Google Scholar] [CrossRef] [PubMed]
  79. Vogiatzakis, I.N. (Ed.) Mediterranean Mountain Environments: Vogiatzakis/Mediterranean Mountain Environments; John Wiley & Sons, Ltd.: Chichester, UK, 2012; ISBN 978-1-119-94115-6. [Google Scholar]
  80. Médail, F.; Diadema, K. Glacial Refugia Influence Plant Diversity Patterns in the Mediterranean Basin. J. Biogeogr. 2009, 36, 1333–1345. [Google Scholar] [CrossRef]
  81. Hughes, P.D.; Woodward, J.C. Quaternary Glaciation in the Mediterranean Mountains: A New Synthesis. Geol. Soc. Lond. Spec. Publ. 2017, 433, 1–23. [Google Scholar] [CrossRef]
  82. Nagy, L.; Grabherr, G. The Biology of Alpine Habitats; OUP Oxford: Oxford, UK, 2009; ISBN 978-0-19-856703-5. [Google Scholar]
  83. Scherrer, D.; Körner, C. Topographically Controlled Thermal-Habitat Differentiation Buffers Alpine Plant Diversity against Climate Warming. J. Biogeogr. 2011, 38, 406–416. [Google Scholar] [CrossRef]
  84. Kulonen, A.; Imboden, R.A.; Rixen, C.; Maier, S.B.; Wipf, S. Enough Space in a Warmer World? Microhabitat Diversity and Small-Scale Distribution of Alpine Plants on Mountain Summits. Divers. Distrib. 2018, 24, 252–261. [Google Scholar] [CrossRef] [Green Version]
  85. Munson, S.M.; Sher, A.A. Long-Term Shifts in the Phenology of Rare and Endemic Rocky Mountain Plants. Am. J. Bot. 2015, 102, 1268–1276. [Google Scholar] [CrossRef] [Green Version]
  86. Thomas, C.D. Climate, Climate Change and Range Boundaries. Divers. Distrib. 2010, 16, 488–495. [Google Scholar] [CrossRef]
  87. Gómez, J.M.; González-Megías, A.; Lorite, J.; Abdelaziz, M.; Perfectti, F. The Silent Extinction: Climate Change and the Potential Hybridization-Mediated Extinction of Endemic High-Mountain Plants. Biodivers. Conserv. 2015, 24, 1843–1857. [Google Scholar] [CrossRef]
  88. Lenoir, J.; Svenning, J.-C. Climate-Related Range Shifts—A Global Multidimensional Synthesis and New Research Directions. Ecography 2015, 38, 15–28. [Google Scholar] [CrossRef]
  89. Urban, M.C. Accelerating Extinction Risk from Climate Change. Science 2015, 348, 571–573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Sandel, B.; Monnet, A.-C.; Govaerts, R.; Vorontsova, M. Late Quaternary Climate Stability and the Origins and Future of Global Grass Endemism. Ann. Bot. 2017, 119, 279–288. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Topographic profile of South Lefka Ori massif and the four Gloria summits (GR-LEO). LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Figure 1. Topographic profile of South Lefka Ori massif and the four Gloria summits (GR-LEO). LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Diversity 13 00022 g001
Figure 2. Vegetation turnover over the two monitoring periods. LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Figure 2. Vegetation turnover over the two monitoring periods. LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Diversity 13 00022 g002
Figure 3. Species diversity per summit per aspect for the two monitoring periods. (N: North, S: South, E: East, W: West, T: Total Summit). LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Figure 3. Species diversity per summit per aspect for the two monitoring periods. (N: North, S: South, E: East, W: West, T: Total Summit). LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Diversity 13 00022 g003
Table 1. Summary information for the four summits: Location, altitude, and size of the SAS *, mean species richness in 2001 and the mean difference between 2008 and 2001.
Table 1. Summary information for the four summits: Location, altitude, and size of the SAS *, mean species richness in 2001 and the mean difference between 2008 and 2001.
LOW **CHOSEKSTR
Location E24.0873 N35.2686E24.0813 N35.2734E24.0892 N35.2819E24.0549 N35.2940
Altitude (m a.s.l.)1664196521602339
Mean species richness per area20012008200120082001200820012008
1 m2 quadrats mean ± 1 SE11.87 ± 0.87+0.063.75 ± 0.59−0.061.50 ± 0.22−0.121.44 ± 0.44 +0.44
5 m SAS (mean of aspects)42.00 ± 1.47−2.2516.50 ± 0.86–2.008.75 ± 0.75+0.506.50 ± 1.32+0.75
10 m SAS (mean of aspects)39.25 ± 0.62−4.0017.00 ± 0.70+1.0010.00 ± 1.29+0.258.50 ± 1.19+2.75
5 + 10 m SAS (mean of aspects)46.75 ± 0.47–2.2520.50 ± 0.64–1.0012.75 ± 0.47–0.759.50 ± 1.50+1.75
* SAS: Summit Area Section, ** LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Table 2. Time series regression for the temperatures in the four summits.
Table 2. Time series regression for the temperatures in the four summits.
Explanatory VariableLOW *CHOSEKSTR
Time t0.000320.0000960.000130.000112
January10.518.2227.7056.266
February10.517.9957.6846.553
March11.398.5856.9725.964
April11.258.3005.7205.422
May12.119.5987.8027.351
June12.059.5378.5147.504
July12.5610.038.7168.117
August12.6910.058.8337.863
September11.418.7697.4816.205
October10.827.9406.7585.488
November9.3096.5315.3294.193
December9.5887.0066.3505.427
sin((2/365)*π*t)−4.976−4.030−4.099−3.848
cos((2/365)*π*t)8.3717.4757.3707.339
* LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Table 3. Number of species in GR-LEO (total) and per summit.
Table 3. Number of species in GR-LEO (total) and per summit.
GR-LEO20012008Trend *LostNew
Total species7066 Diversity 13 00022 i00184
Endemics2018 Diversity 13 00022 i00131
Summits
LOW **5955 Diversity 13 00022 i00162
CHO3228 Diversity 13 00022 i00173
SEK1818--22
STR1415 Diversity 13 00022 i00201
* red arrow indicates decreasing trend while green arrow increasing trend, ** LOW: low summit, CHO: Chorafas, SEK: South East Kakovoli, STR: Sternes.
Table 4. Statistical tests on species richness data comparing 2001 and 2008 temperature.
Table 4. Statistical tests on species richness data comparing 2001 and 2008 temperature.
Site/Sampling ScaleModelt-Valuep-Value
LOW 5 m SAS *Pair t-test1.09440.133
LOW 10 m SASPair t-test3.70330.003 *
CHO 5 m SASPair t-test2.82840.0633
CHO 10 m SASPair t-test0.56240.6134
SEK 10 m SASPair t-test1.73210.1817
SEK 10 m SASPair t-test0.26410.8089
STR 5 m SASPair t-test1.56670.2151
STR 10 m SASPair t-test1.96120.1447
Site/Sampling ScaleModelR2sig.level
SAS 5–All sitesLinear regression against temp. means0.780.001
SAS 5–All sitesLinear regression against min temp. means0.850.001
SAS 5–All sitesLinear regression against max temp. means0.630.001
SAS 10–All sitesLinear regression against temp. means0.860.001
SAS 10–All sitesLinear regression against min temp. means0.90.001
SAS 10–All sitesLinear regression against max temp. means0.760.001
SAS–5 SR changes–all sitesLinear regression against temp. means0.340.05
SAS–5 SR changes–all sitesLinear regression against min temp. means0.360.05
SAS–5 SR changes–all sitesLinear regression against max temp. means0.270.05
SAS–5 SR changes–all sitesLinear regression against June temp. means0.460.05
SAS–10 SR changes–all sitesLinear regression against temp. means0.11n.s
SAS–10 SR changes–all sitesLinear regression against min temp. means0.16n.s
SAS–10 SR changes–all sitesLinear regression against max temp. means0.13n.s
SAS–10 SR changes–all sitesLinear regression against June temp. means0.14n.s
* SAS: Summit Area Section.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kazakis, G.; Ghosn, D.; Remoundou, I.; Nyktas, P.; Talias, M.A.; Vogiatzakis, I.N. Altitudinal Vascular Plant Richness and Climate Change in the Alpine Zone of the Lefka Ori, Crete. Diversity 2021, 13, 22. https://0-doi-org.brum.beds.ac.uk/10.3390/d13010022

AMA Style

Kazakis G, Ghosn D, Remoundou I, Nyktas P, Talias MA, Vogiatzakis IN. Altitudinal Vascular Plant Richness and Climate Change in the Alpine Zone of the Lefka Ori, Crete. Diversity. 2021; 13(1):22. https://0-doi-org.brum.beds.ac.uk/10.3390/d13010022

Chicago/Turabian Style

Kazakis, George, Dany Ghosn, Ilektra Remoundou, Panagiotis Nyktas, Michael A. Talias, and Ioannis N. Vogiatzakis. 2021. "Altitudinal Vascular Plant Richness and Climate Change in the Alpine Zone of the Lefka Ori, Crete" Diversity 13, no. 1: 22. https://0-doi-org.brum.beds.ac.uk/10.3390/d13010022

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