Distributional responses to climate change for alpine species of Cyananthus and Primula endemic to the Himalaya-Hengduan Mountains
Xie Hea,b,d, Kevin S. Burgessc, Lian-Ming Gaob*, De-Zhu Lia,b,d**     
a. Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China;
b. CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China;
c. Department of Biology, College of Letters & Sciences, Columbus State University, University System of Georgia, Columbus, GA, 31907-5645, USA;
d. Kunming College of Life Sciences, University of Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
Abstract: Global warming increases the vulnerability of plants, especially alpine herbaceous species, to local extinction. In this study, we collected species distribution information from herbarium specimens for ten selected Cyananthus and Primula alpine species endemic to the Himalaya-Hengduan Mountains (HHM). Combined with climate data from WorldClim, we used Maximum Entropy Modeling (MaxEnt) to project distributional changes from the current time period to 2070. Our predictions indicate that, under a wide range of climate change scenarios, the distributions of all species will shift upward in elevation and northward in latitude; furthermore, under these scenarios, species will expand the size of their range. For the majority of the species in this study, habitats are available to mitigate upward and northward shifts that are projected to be induced by changing climate. If current climate projections, however, increase in magnitude or continue to increase past our projection dates, suitable habitat for future occupation by alpine species will be limited as we predict range contraction or less range expansion for some of the species under more intensified climate scenarios. Our study not only underscores the value of herbarium source information for future climate model projections but also suggests that future studies on the effects of climate change on alpine species should include additional biotic and abiotic factors to provide greater resolution of the local dynamics associated with species persistence under a warming climate.
Keywords: Global warming    MaxEnt    Range expansion    Alpine species    Himalaya-Hengduan Mountains    
1. Introduction

The global mean temperature has increased by approximately 0.74 ℃ in the last one hundred years and this warming trend is projected to accelerate in the coming decades, increasing 2.8-5.3 ℃ by 2085 (Nogués-Bravo et al., 2007; Pacifici et al., 2015). Numerous studies have shown that in response to such warming trends the distribution of plant species tends to shift towards suitable habitat, specifically, upward in elevation and northward in latitude (Hickling et al., 2006; Chen et al., 2011). Given that the evolution of climatic niche partitioning tends to be much slower than rates of climate change (Quintero and Wiens, 2013), climate change-induced contractions in range sizes threaten many species with local extinction (Wiens, 2016). However, in response to similar climate changes, range expansion has been documented (Berry et al., 2003; Hamann and Wang, 2006). Thus, understanding how species ranges respond to climate change remains a major challenge.

Alpine species are considered especially sensitive to climate change (Lenoir et al., 2008). Due to their relatively low biotic complexity, climate change is suspected to be a dominant factor affecting range shifts in alpine species, which are often accompanied with marginal extinction at the warm edge and expansion at the cold edge (Parmesan, 2006; Pauli et al., 2007). In addition, because species distributed at high elevations are usually endemic to specific regions and have relatively narrow climate tolerances, theyoften experience a "nowhere-to-go" scenario, i.e., cold habitats are limited to the summits of mountainous regions (Nogués-Bravo et al., 2007). Models for Primula and Gentiana species predict decreased niches in response to climate change and grazing enforcement (Wehn and Johansen, 2015; Fois et al., 2016). Five comprehensive model projections have also shown that approximately 1/3 to 1/2 of the species distributed at higher elevations will lose more than 80% of their suitable habitats by 2070-2100, and cold-adapted mountain flora are at the highest risk for local extinction (Engler et al., 2011).

The Himalaya-Hengduan Mountains (HHM) region stretches over two global biodiversity hotspots, the Mountains of Southwest China and the Himalayas, both of which harbor a diversity of habitats and a high level of plant biodiversity (Mittermeier, 2004). A large proportion of alpine and sub-nival plant species are endemic to the region (Xu et al., 2014) and the area has experienced extensive glacial retreat and tree line advancement in the last century (Baker and Moseley, 2007). Two alpine genera in particular, Cyananthus and Primula, are suspected to be at risk of extinction due to climate change in the HHM (Zhou et al., 2013; Yan et al., 2015). The genus Cyananthus (Campanulaceae) comprises approximately 26 species with the highest level of diversity in the HHM (Zhou et al., 2013) and Primula (Primulaceae) contains about 500 species, among which 200 species are mainly distributed in this region (Yan et al., 2015). Most of the species in these two genera are perennial herbs, and although climate change is suspected to be a significant threat to remaining populations, there are few empirical studies that have investigated how these species may respond to climate change.

In this study, we explore how the distribution of selected Cyanathus and Primula alpine species endemic to the HHM may respond to climate change. Specifically, we extract species distribution information from herbarium specimens collected in the last century and, combined with climate data from WorldClim, we use species distribution models (SDM) to project species distributional dynamics under future climate warming conditions.

2. Materials and methods 2.1. Herbarium data

Species occurrence data were obtained from the Chinese Virtual Herbarium (CVH: http://www.cvh.ac.cn/). To assess possible changes in distribution for alpine species, we targeted Primula and Cyananthus species that had at least 100 specimen records distributed in the HHM: five species from each genus satisfied this criterion; their elevations ranged from ~2800 m to 4000 m. The original dataset contained 2, 866 specimen records, ranging from 123 specimens (Cyananthus lobatus) to 614 specimens (Primula sonchifolia). The following specimens were removed from the dataset: duplicate records with the same collector and collection number; specimens with elevation or collection records located far beyond the species distribution ranges based on the Flora of China; specimens with obvious identification errors; specimens that lacked spatial coordinates for georeferencing in Google Earth; and, to reduce the effect of spatial autocorrelation and the consequent overfitting, records within five kilometers of each other (Veloz, 2009). Our final dataset contained 708 records, which were used to conduct species distribution modeling (Table 1, Fig. 1).

Table 1 Number of original Cyananthus and Primula specimens and the number of their occurrences with spatial coordinates after data filtration. Included are the Threshold values used to distinguish the area for which each species was present and the AUC (Area Under Receiver Operator Curve) value of each species' model.
Species Number of specimens Number of occurrences Threshold values AUC
C. hookeri 306 63 0.2057 0.988 ± 0.006
C. incanus 254 66 0.3074 0.980 ± 0.023
C. inflatus 301 57 0.2206 0.991 ± 0.006
C. lobatus 123 26 0.3904 0.978 ± 0.037
C. macrocalyx 289 39 0.0980 0.991 ± 0.009
P. poissonii 259 16 0.2723 0.994 ± 0.006
P. secundifolia 231 16 0.3488 0.984 ± 0.009
P. sikkimensis 285 62 0.0792 0.988 ± 0.008
P. sonchifolia 614 307 0.1482 0.993 ± 0.005
P. tangutica 204 56 0.2150 0.991 ± 0.004

Fig. 1 Model projections were based on the occurrences of five species of Cynananthus (a) and five species of Primula (b) in the Himalaya-Hengduan Mountains.
2.2. Environmental data

Climate data were downloaded from WorldClim ver. 1.4 (http://www.worldclim.org/). We downloaded 19 bioclimatic variables (Table S1) for both current (Current: the average for 1960-1990) and future (2070: the average for 2060-2080) conditions at the highest available spatial resolution (30 arc-seconds; ~1 km). We chose the Global Circulation Model (GCM) ACCESS1-0, which has been shown to perform the best in the HHM region (Zhang et al., 2015; Wu et al., 2017). We used two greenhouse gas concentration trajectories for future climate conditions: an optimistic scenario whereby emissions peak around 2040 and then decline (Representative Concentration Pathway (RCP 4.5), and a pessimistic scenario whereby emissions continue to rise throughout the century (RCP 8.5).

We performed a Pearson's correlation test for the 19 climatic variables (Table S2) and the variables with correlation coefficients above 0.90 were removed to reduce the effect of multicollinearity. After this procedure, eight bioclimatic variables remained in our analyses, among which four variables were associated with temperature (bio1, bio2, bio3, bio4) and another four variables were associated with precipitation (bio12, bio14, bio15, bio18). Additional static environmental variables previously shown to contribute to model performance were used in the analysis (Stanton et al., 2012), including a land-cover layer (resolution of 10 arc-seconds) downloaded from the European Space Agency (ESA: http://due.esrin.esa.int/page_globcover.php), a soil layer (resolution of 30 arc-seconds) downloaded from the International Institute for Applied Systems Analysis (IIASA: http://webarchive.iiasa.ac.at/Research/LUC/External-World-soil-database/HTML/index.html?sb=1), and an elevation layer (resolution of 30 arc-seconds) downloaded from the Geospatial Information Authority of Japan (GIAJ: https://globalmaps.github.io/el.html). These variables were assumed to be constant under current and future scenarios. All layers were cut to restrict the projection area to the HHM, standardized to the same resolution as the climatic layers (30 arcseconds) with the same coordinate system (WGS1984), and transferred to ASCII format to be prepared for modeling operations.

2.3. MaxEnt modeling and output data analysis

We used MaxEnt to project species distribution under current and 2070 climate conditions (Phillips et al., 2006; Elith et al., 2011). Occurrence data and environmental layers for the current time period and the 2070 projection were used as input data in the model. We used the default settings for MaxEnt parameters and set the iteration time as 10 for each species, then we used the average logistic output layer (based on the ten replicates of cross-validation runs for each species) for subsequent calculations. Values for area under the receiver operating characteristic curves (AUC) were used to evaluate the performance of the models, where the value 0.5 indicates a random prediction and 1.0 a perfect model prediction (Bai et al., 2018). We reclassified the MaxEnt output file using the 10-percentile training presence logistic (10TPL) threshold value (Table 1) to define a species potential distribution region, above which species were considered "present" in the region, a method widely recognized for distinguishing suitable from unsuitable regions (Escalante et al., 2013; Kramer-Schadt et al., 2013; Radosavljevic et al., 2014; Hughes, 2017). We then calculated the longitude, latitude, elevation and range size for each pixel of potential presence and estimated the mean values for the current time period and 2070 for each species. The average value of longitude, latitude and elevation between the current time period and 2070 were compared within and among species. The shifts in total range size from the current time period to 2070 were also estimated. We then overlaid the distribution maps for the two time periods to detect potential changes in the occupied area.

All analyses were performed using MaxEnt 3.3.3k (Phillips et al., 2006), R ×64 3.3.3 (R Core Team, 2016), and ArcGIS10.2 (Environmental Systems Resource Institute (ESRI), 2014).

3. Results

All models in this study had AUC values > 0.978 (Table 1), indicating high performance. The mean 10TPL value was 0.229 (Table 1), which means that on average, we considered the cell with occurrence possibility above 0.229 to be suitable for species survival.

For all species, our models predict that, in response to climate change, distributions will shift upward in elevation, northward in latitude, and westward in longitude. Furthermore, these changes in distribution are predicted to be stronger in more pessimistic scenarios such as RCP 8.5. The mean shifts in elevation across the ten species were 532.49 m in the RCP 4.5 scenario (range: 348.11 m for Primula sikkimensis to 871.38 m for Primula poissonii) and 661.43 m in the RCP 8.5 scenario (range: 419.17 m for P. sikkimensis to 1042.38 m for P. poissonii). Mean shifts in latitude were 1.36° under the RCP 4.5 scenario (range: 0.82° for C. lobatus to 1.83° for P. sikkimensis) and 1.61° under the RCP 8.5 scenario (range: 1.03° for C. lobatus to 2.28° for P. sikkimensis) (Tables 2 and 3, Fig. 2).

Table 2 Projected changes in range size and elevation for Cyananthus and Primula species in the Himalaya-Hengduan Mountains. Listed are the estimated Range sizes (% change in range size) and mean Elevations (m ± standard error) in current and future (2070) projections under the optimistic (RCP 4.5) and the pessimistic (RCP 8.5) models for climate change.
Species Range size Elevation
Current 2070 RCP 4.5 2070 RCP 8.5 Current 2070 RCP 4.5 2070 RCP 8.5
C. hookeri 520818.64 865587.08 (66.2) 958649.68 (84.1) 3760.17 ± 612.13 4119.43 ± 626.8 4251.78 ± 596.61
C. incanus 402835.87 928989.23 (130.6) 1007652.52 (150.1) 4028.67 ± 662.1 4417.64 ± 576.73 4576.02 ± 514.01
C. inflatus 315840.63 361409.43 (14.4) 280491.98 (-11.2) 2896.07 ± 808.34 3561.06 ± 793.91 3690.23 ± 754.99
C. lobatus 191350.21 373843.47 (95.4) 399883.09 (109.0) 3736.28 ± 573.75 4151.39 ± 662.57 4184.38 ± 657.76
C. macrocalyx 340121.21 675029.92 (98.5) 584848.57 (72.0) 3595.87 ± 855.39 4388.12 ± 695.65 4596.41 ± 624.4
P. poissonii 214198.85 379857.37 (77.3) 403875.45 (88.6) 2848.15 ± 909.2 3719.53 ± 977.33 3890.53 ± 943.81
P. secundifolia 401837.26 675402.84 (68.1) 705069.54 (75.5) 3739.06 ± 779.75 4201.97 ± 731.96 4289.41 ± 713.41
P. sikkimensis 351578.17 1029850.46 (192.9) 1141228.36 (224.6) 3920.62 ± 580 4268.73 ± 648.09 4339.79 ± 688.74
P. sonchifolia 277936.42 402046.28 (44.7) 309688.54 (11.4) 3276.31 ± 818.39 3937.41 ± 720.87 4150.97 ± 657.89
P. tangutica 426533.12 476937.34 (11.8) 376277.5 (-11.8) 3549.12 ± 641.05 3909.9 ± 665.11 3995.07 ± 631.35

Table 3 Projected changes in longitude and latitude for Cyananthus and Primula species in the Himalaya-Hengduan Mountains. Listed are the estimated mean Longitude (±standard error) and mean Latitude (±standard error) in current and future (2070) projections under the optimistic (RCP 4.5) and the pessimistic (RCP 8.5) models for climate change.
Species Longitude Latitude
Current 2070 RCP 4.5 2070 RCP 8.5 Current 2070 RCP 4.5 2070 RCP 8.5
C. hookeri 97.44 ± 5.51 96.19 ± 5.75 95.47 ± 5.72 30.18 ± 2.25 31.13 ± 2.56 31.26 ± 2.5
C. incanus 96.12 ± 4.91 94.86 ± 5.58 93.78 ± 5.55 29.38 ± 1.39 30.93 ± 2.03 31.23 ± 2.11
C. inflatus 100.12 ± 3.98 98.47 ± 5.06 98.02 ± 4.89 27.29 ± 1.95 28.87 ± 1.93 28.97 ± 1.71
C. lobatus 94.22 ± 7.02 92.62 ± 7.01 92.9 ± 6.85 28.55 ± 1.46 29.37 ± 1.47 29.58 ± 1.57
C. macrocalyx 99.17 ± 4.44 95.44 ± 5.94 94.56 ± 5.89 28.76 ± 1.8 30.22 ± 1.86 30.28 ± 1.68
P. poissonii 100.65 ± 2.64 97.84 ± 4.09 97.35 ± 4.28 26.96 ± 1.76 28.63 ± 1.7 29.05 ± 1.53
P. secundifolia 98.21 ± 4.81 95.89 ± 6.06 95.57 ± 5.98 29.33 ± 1.88 30.25 ± 1.97 30.38 ± 1.9
P. sikkimensis 97.17 ± 5.87 95.32 ± 6.8 94.58 ± 7.53 30.41 ± 2.26 32.24 ± 2.83 32.69 ± 3.06
P. sonchifolia 100.27 ± 4.24 99.02 ± 5.64 98.31 ± 6.12 28.57 ± 2.11 30.4 ± 2.35 30.5 ± 2.21
P. tangutica 100.67 ± 3.31 99.23 ± 2.89 98.9 ± 2.57 34.32 ± 2.53 35.34 ± 2.29 35.86 ± 2.45

Fig. 2 The projected distributional maps for the ten studied species. Shown is the area to be lost (red region), area to remain unchanged (blue region) and area to be gained (purple region) from the current time period to 2070, under two climate scenarios RCP 4.5 and RCP 8.5.

In the RCP 4.5 scenario, our model predicts that from now until 2070 the range for all species will expand (range: 11.8% for Primula tangutica to 192.9% for P. sikkimensis); moreover, the newly occupied area will be larger than the area lost. Under the RCP 8.5 scenario, six of the 10 species showed greater range expansion than under the RCP 4.5 scenario (range: 75.5% for Primula secundifolia to 224.6% for P. sikkimensis), two species showed less range expansion (72.0% for Cyananthus macrocalyx and 11.4% for P. sonchifolia), while Cyananthus inflatus and P. tangutica were projected to lose area by 11.2% and 11.8%, respectively (Table 2, Fig. 2). Loss in distribution area was mainly confined to the southern portion of species distribution ranges, whereas gains in distribution area were mainly confined to the northwest on the Qinghai-Tibet Plateau.

4. Discussion 4.1. Changes in species distribution

Our models predict that the distribution of all Cyanathus and Primula species in our study will show a consistent response to a warming climate; specifically, the distributions of these species will shift upward in elevation and northward in latitude so that species will occupy living conditions similar to that of their previous habitats. Our predictions are consistent with previous studies that modeled distributional changes in response to climate change for other species in the region. For example, northward and upward shifts in distribution have been detected in wild soybean (He et al., 2016) and yews (Poudel et al., 2014). Furthermore, the westward shifts in distribution projected by our models are not surprising. Because high mountain habitats converge in the western portion of the HHM, when species distributions shift upwards in elevation, they will likely shift westward in longitude. This same trend has been projected for other alpine species in the HHM (Liang et al., 2018; He et al., unpublished).

4.2. Range shifts under future scenarios of climate change

All the species in our study are projected to expand their ranges under the more optimistic RCP 4.5 scenario but some of the species are projected to show greater range expansion under the more pessimistic RCP 8.5 scenario. These projections are in stark contrast to those of previous studies that have projected decreases in species ranges due to climate change. For example, the distribution of beech trees is expected to decrease by 2100 due to increased climate variability (Saltre et al., 2015), and by 2070 a decline in suitable climate space is projected for the two threatened dipterocarp trees (Deb et al., 2017) and the Tertiary relict tree species Davidia involucrate (Tang et al., 2017). Under the global warming scenarios projected in our study, however, suitable habitats may not be limited for most of the alpine species, which appear to be able to expand their distribution under our model projections. Similarly, a high proportion of plant species (63.6%) in the HHM have shown range expansion under future climate warming scenarios, but only a small proportion of plant species (11.9%) have been projected to experience range contraction (Liang et al., 2018). According to our model, the range of two species (C. inflatus and P. tangutica), which expand in the RCP 4.5 scenario, are projected to contract under the RCP 8.5 scenario. Under this same scenario, the ranges of another two species (C. macrocalyx and P. sonchifolia) are expected to expand less. Unsurprisingly, the species (P. tangutica) that occupies the highest latitude in our study is projected to encounter the greatest limitations to range expansion, likely due to the lack of suitable habitat in the north. While other three species (C. inflatus, C. macrocalyx and P. sonchifolia) that occupy relatively low latitudes and low elevations are projected to lose the largest areas in the southern margin of their distribution. Collectively, our projections indicate that although there is suitable habitat for future range expansion, if more extreme climatic changes persist, some of the alpine species in our study may decrease in range size, possibly resulting in local extinction.

We also found that the expansion area was in most cases larger than the contraction area, indicating that most of the species in our study will actually increase their area of occupied sites in response to climate change. Other models have projected that in response to a warming climate the land surface area suitable for plant occupation in some mountain habitats would increase rather than decrease (Elsen and Tingley, 2015). For example, Liang et al. (2018) concluded that, in the Hengduan Mountains, species would expand their ranges because the complex topography in high elevation areas provide sufficient land surfaces, which would promote not only upward shifts in elevation, but also expansion into the adjacent Qinghai-Tibet Plateau, which can accommodate northward shifting species and mitigate range contractions. Similarly, most Rhodiola species in the HHM are expected to expand their range sizes through upward shifts in elevation, as climate warming is expected to reduce permafrost and increase potential habitats for upshifting species (You et al., 2018). In this study, northwestern Yunnan, western Sichuan, and southeastern Tibet appear to be stable areas with suitable habitats that can accommodate alpine herbaceous species in response to climate change. These regions may be particularly useful for future in situ conservation efforts (Fois et al., 2017; Heywood, 2018), whereas immediate actions to minimize the effects of grazing and tourism may help to provide future habitats for species affected by climate change.

4.3. Future considerations

In this study, we assumed that potential habitats could be fully occupied by species with full dispersal capability when only considered the impacts of climate change. Our projections, however, did not consider several biotic factors that may impact our results, such as limited dispersal capability, pollination syndrome, and competition among species. For example, a loss of distribution area of up to 61.10% (full dispersal assumption) and 68.91% (no dispersal) was observed in SDM projections of crop wild relatives (Aguirre-Gutiérrez et al., 2017); species range sizes may also be controlled by mating systems, where self-pollinated species typically have larger geographic ranges than those of close relatives that are out-crossing (Grossenbacher et al., 2015). Although these sources of variation may affect the outcomes of our projections on a local scale, we anticipate that the general trends shown here may not change as climate is likely a dominant factor leading species distributional changes at a broader regional scale (Pearson and Dawson, 2003).

The model projections in this study span ten species across two genera and reveal almost the same distributional changes across all species. The consistency of these projections provides confidence in the validity of our models. However, although our modeling process was previously shown to perform well in the HHM region, our projections rely on this single GCM. To expand the confidence in model projections, future studies should include more GCMs to evaluate the magnitude of variation that may be associated with different climate models.

The species selected in our study are mainly those with an abundant number of herbarium specimens in order to have enough valid data for modeling projections. As such, the species included in our study are those with a relatively wide distribution and abundance. Species with a more localized distribution and less abundance are largely ignored in this analysis due to the limited number of herbarium specimens available for data collection. Such species, however, may be the most susceptible to climate change due to their narrow niche occupation and habitat requirements. This limitation highlights the need for future studies to include fieldbased investigations on the distribution of locally endemic species and their fate under future climate change scenarios.

5. Conclusion

Our study highlights the usefulness of herbarium specimen data for evaluating species distribution and predicting species distributional changes under the effects of climate change. The results are consistent with previous studies that show species will shift upward in elevation, northward in latitude and westward in longitude under a warming climate scenario, however, our species are also projected to expand their ranges as the loss of previously occupied habitat is less than that of newly acquired areas. Great opportunities for the exploration of species distributional dynamics under climate change remain to be explored, including using digitalized herbarium specimens, citizen scientists, and targeted field surveys. These new approaches to examining responses to climate change may provide insights for targeted species protection. The utilization of more comprehensive models that consider both the biotic and abiotic factors contributing to species persistence will provide valuable data for making conservation decisions under future climate scenarios.

This study was funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB31010000), the National Key Basic Research Program of China (2014CB954100), and the Program of Science and Technology Talents Training of Yunnan Province (2017HA014). We appreciate the Chinese Virtual Herbarium for providing specimen data of Cyananthus and Primula species.

Supplementary data to this article can be found online at https://doi.org/10.1016/j.pld.2019.01.004.

References
Aguirre-Gutiérrez J. van Treuren R., Hoekstra R., van Hintum T.J.L., Vaclavik T., 2017. Crop wild relatives range shifts and conservation in Europe under climate change. Divers. Distrib, 23(7): 739-750. DOI:10.1111/ddi.12573
Bai Y., Wei X., Li X., 2018. Distributional dynamics of a vulnerable species in response to past and future climate change:a window for conservation prospects. PeerJ, 6: e4287. DOI:10.7717/peerj.4287
Baker B.B., Moseley R.K., 2007. Advancing treeline and retreating glaciers:Implications for conservation in Yunnan, P. R. China. Arctic Antarct. Alpine Res, 39(2): 200-209. DOI:10.1657/1523-0430(2007)39[200:ATARGI]2.0.CO;2
Berry P.M., Dawson T.P., Harrison P.A., Pearson R., Butt N., 2003. The sensitivity and vulnerability of terrestrial habitats and species in Britain and Ireland to climate change. J. Nat. Conserv, 11(1): 15-23. DOI:10.1078/1617-1381-00030
Chen I.C., Hill J.K., Ohlemuller R., Roy D.B., Thomas C.D., 2011. Rapid range shifts of species associated with high levels of climate warming. Science, 333(6045): 1024-1026. DOI:10.1126/science.1206432
Deb J.C., Phinn S., Butt N., McAlpine C.A., 2017. The impact of climate change on the distribution of two threatened Dipterocarp trees. Ecol. Evol, 7(7): 2238-2248. DOI:10.1002/ece3.2846
Elith J., Phillips S.J., Hastie T., Dudík M., Chee Y.E., Yates C.J., 2011. A statistical explanation of MaxEnt for ecologists. Divers. Distrib, 17(1): 43-57. DOI:10.1111/j.1472-4642.2010.00725.x
Elsen P.R., Tingley M.W., 2015. Global mountain topography and the fate of montane species under climate change. Nat. Clim. Change, 5(6): 772-776. DOI:10.1038/nclimate2656
Engler R., Randin C.F., Thuiller W., Dullinger S., Zimmermann N.E., AraÚJo M.B., et al, 2011. 21st century climate change threatens mountain flora unequally across Europe. Glob. Chang. Biol, 17(7): 2330-2341. DOI:10.1111/j.1365-2486.2010.02393.x
Environmental Systems Resource Institute (ESRI), 2014. ArcMap 10.2. ESRI, Redlands, California.
Escalante T., Rodríguez-Tapia G., Linaje M., Illoldi-Rangel P., González-López R., 2013. Identification of areas of endemism from species distribution models:threshold selection and Nearctic mammals. Tip, 16(1): 5-17. DOI:10.1016/s1405-888x(13)72073-4
Fois M., Bacchetta G., Cogoni D., Fenu G., 2017. Current and future effectiveness of the Natura 2000 network for protecting plant species in Sardinia:a nice and complex strategy in its raw state? J. Environ. Plann. Manag, 61(2): 332-347. DOI:10.1080/09640568.2017.1306496
Fois M., Cuena-Lombra#241;a A. Fenu G., Cogoni D., Bacchetta G., 2016. The reliability~of conservation status assessments at regional level:Past, present and future perspectives on Gentiana lutea L. ssp. lutea in Sardinia. J. Nat. Conserv, 33: 1-9. DOI:10.1016/j.jnc.2016.06.001
Grossenbacher D., Briscoe Runquist R., Goldberg E.E., Brandvain Y., 2015. Geographic range size is predicted by plant mating system. Ecol. Lett, 18(7): 706-713. DOI:10.1111/ele.12449
Hamann A., Wang T., 2006. Potential effects of climate change on ecosystem and tree species distribution in British Columbia. Ecology, 87(11): 2773-2786. DOI:10.1890/0012-9658(2006)87[2773:peocco]2.0.co;2
He S.L., Wang Y.S., Li D.Z., Yi T.S., 2016. Environmental and historical determinants of patterns of genetic differentiation in wild soybean (Glycine soja Sieb. et Zucc).. Sci. Rep, 6: 22795. DOI:10.1038/srep22795
Heywood V.H., 2018. Conserving plants within and beyond protected areas e still problematic and future uncertain. Plant. Divers.. DOI:10.1016/j.pld.2018.10.001
Hickling R., Roy D.B., Hill J.K., Fox R., Thomas C.D., 2006. The distributions of a wide range of taxonomic groups are expanding polewards. Glob. Chang. Biol, 12(3): 450-455. DOI:10.1111/j.1365-2486.2006.01116.x
Hughes A.C., 2017. Mapping priorities for conservation in Southeast Asia. Biol.Conserv, 209: 395-405. DOI:10.1016/j.biocon.2017.03.007
Kramer-Schadt S., Niedballa J., Pilgrim J.D., Schröder B., Lindenborn J., Reinfelder V., et al, 2013. The importance of correcting for sampling bias in MaxEnt species distribution models. Divers. Distrib, 19(11): 1366-1379. DOI:10.1111/ddi.12096
Lenoir J., Gegout J.C., Marquet P.A., de Ruffray P., Brisse H., 2008. A significant upward shift in plant species optimum elevation during the 20th century. Science, 320(5884): 1768-1771. DOI:10.1126/science.1156831
Liang Q., Xu X., Mao K., Wang M., Wang K., Xi Z., et al, 2018. Shifts in plant distributions in response to climate warming in a biodiversity hotspot, the Hengduan Mountains. J. Biogeogr, 45(6): 1334-1344. DOI:10.1111/jbi.13229
Mittermeier, R.A., 2004. Hotspots revisited. Cemex.
Nogués-Bravo D. Araújo M.B., Errea M.P., Martínez-Rica J.P., 2007. Exposure of global mountain systems to climate warming during the 21st century. Glob. Environ.Change, 17(3-4): 420-428. DOI:10.1016/j.gloenvcha.2006.11.007
Pacifici M., Foden W.B., Visconti P., Watson J.E.M., Butchart S.H.M., Kovacs K.M., et al, 2015. Assessing species vulnerability to climate change. Nat. Clim. Change, 5(3): 215-224. DOI:10.1038/nclimate2448
Parmesan C., 2006. Ecological and evolutionary responses to recent climate change. Annu. Rev. Ecol. Evol. Syst, 37(1): 637-669. DOI:10.1146/annurev.ecolsys.37.091305.110100
Pauli H., Gottfried M., Reiter K., Klettner C., Grabherr G., 2007. 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, 13(1): 147-156. DOI:10.1111/j.1365-2486.2006.01282.x
Pearson R.G., Dawson T.P., 2003. Predicting the impacts of climate change on the distribution of species:are bioclimate envelope models useful? Glob. Ecol.Biogeogr, 12(5): 361-371. DOI:10.1046/j.1466-822X.2003.00042.x
Phillips S.J., Anderson R.P., Schapire R.E., 2006. Maximum entropy modeling of species geographic distributions. Ecol. Model, 190(3-4): 231-259. DOI:10.1016/j.ecolmodel.2005.03.026
Poudel R.C., Möller M., Liu J., Gao L.-M., Baral S.R., Li D.-Z., et al, 2014. Low genetic diversity and high inbreeding of the endangered yews in Central Himalaya:implications for conservation of their highly fragmented populations. Divers. Distrib, 20(11): 1270-1284. DOI:10.1111/ddi.12237
Quintero I., Wiens J.J., 2013. Rates of projected climate change dramatically exceed past rates of climatic niche evolution among vertebrate species. Ecol. Lett, 16(8): 1095-1103. DOI:10.1111/ele.12144
R Core Team, 2016. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/.
Radosavljevic A., Anderson R.P., Araújo M., 2014. Making better Maxent models of species distributions:complexity overfitting and evaluation. J. Biogeogr, 41(4): 629-643. DOI:10.1111/jbi.12227
Saltre F., Duputie A., Gaucherel C., Chuine I., 2015. How climate, migration ability and habitat fragmentation affect the projected future distribution of European beech. Glob. Chang. Biol, 21(2): 897-910. DOI:10.1111/gcb.12771
Stanton J.C., Pearson R.G., Horning N., Ersts P., Reşit Akçakaya H., 2012. Combining static and dynamic variables in species distribution models under climate change. Methods Ecol. Evol, 3(2): 349-357. DOI:10.1111/j.2041-210X.2011.00157.x
Tang C.Q., Dong Y.F., Herrando-Moraira S., Matsui T., Ohashi H., He L.Y., et al, 2017. Potential effects of climate change on geographic distribution of the Tertiary relict tree species Davidia involucrata in China. Sci. Rep, 7: 43822. DOI:10.1038/srep43822
Veloz S.D., 2009. Spatially autocorrelated sampling falsely inflates measures of accuracy for presence-only niche models. J. Biogeogr, 36(12): 2290-2299. DOI:10.1111/j.1365-2699.2009.02174.x
Wehn S., Johansen L., 2015. The distribution of the endemic plant Primula scandinavica, at local and national scales, in changing mountainous environments. Biodiversity, 16(4): 278-288. DOI:10.1080/14888386.2015.1116408
Wiens J.J., 2016. Climate-related local extinctions are already widespread among plant and animal species. PLoS Biol, 14(12): e2001104. DOI:10.1371/journal.pbio.2001104
Wu Q., Jiang X.W., Xie J., 2017. Evaluation of surface air temperature in Southwestern China simulated by the CMIP5 models. Plateau Meteorol, 36(2): 358-370. DOI:10.7522/j.issn.1000-0534.2016.00046
Xu B., Li Z.M., Sun H., 2014. Plant diversity and floristic characters of the alpine subnival belt flora in the Hengduan Mountains, SW China. J. Systemat. Evol, 52(3): 271-279. DOI:10.1111/jse.12037
Yan H.F., Liu Y.J., Xie X.F., Zhang C.Y., Hu C.M., Hao G., et al, 2015. DNA barcoding evaluation and its taxonomic implications in the species-rich genus Primula L.in China. PLoS One, 10(4): e0122903. DOI:10.1371/journal.pone.0122903
You J., Qin X., Ranjitkar S., Lougheed S.C., Wang M., Zhou W., et al, 2018. Response to climate change of montane herbaceous plants in the genus Rhodiola predicted by ecological niche modelling. Sci. Rep, 8(1): 5879. DOI:10.1038/s41598-018-24360-9
Zhang W.L., Zhang J.Y., Fan G.Z., 2015. Evaluation and projection of dry- and wetseason precipitation in southwestern China using CMIP5 models. Chin. J. Atmos.Sci, 39(3): 559-570.
Zhou Z., Hong D., Niu Y., Li G., Nie Z., Wen J., et al, 2013. Phylogenetic and biogeographic analyses of the Sino-Himalayan endemic genus Cyananthus(Campanulaceae) and implications for the evolution of its sexual system. Mol.Phylogenet. Evol, 68(3): 482-497. DOI:10.1016/j.ympev.2013.04.027