朱淑霞 朱弘 程琳 伊賢貴 王賢榮

摘?要：該研究以高盆櫻桃（Cerasus cerasoides）及其近緣種鐘花櫻桃（C. campanulata）為對(duì)象，通過(guò)收集其實(shí)際分布的地理坐標(biāo)，運(yùn)用BIOCLIM模型模擬現(xiàn)代適生區(qū)范圍并預(yù)測(cè)其未來(lái)氣候下（CCM3，2100）潛在分布區(qū)的變化;結(jié)合主成分分析和相關(guān)性分析確定其主導(dǎo)氣候因子，并比較兩者在氣候限制因子方面的差異;利用受試者工作特征曲線(xiàn)（Receiver operation characteristic， ROC）評(píng)估模型預(yù)測(cè)效果。結(jié)果表明：（1）高盆櫻桃和鐘花櫻桃主要分布于我國(guó)長(zhǎng)江以南大部分省區(qū)，兩者現(xiàn)代分布中心分布位于云貴高原，武夷山脈和南嶺山脈。（2）未來(lái)氣候變化情境下（CCM3），高盆櫻桃和鐘花櫻桃的適生范圍將縮小。二者在中國(guó)西南部（高盆櫻桃）和東南部（鐘花櫻桃）的適生范圍可能大幅減小，而鐘花櫻桃在湖南西部的適生區(qū)可能增加。（3）主成分（PCA）及相關(guān)性分析表明，年降水量（bio12）、最濕季降水量（bio16）、最暖季降水量（bio18）、溫度季節(jié)變化方差（bio4）是影響高盆櫻桃及鐘花櫻桃當(dāng)下適生區(qū)的主要?dú)夂蛞蜃?，“熱量變異幅度”是造成二者分布存在差異的最主要環(huán)境因子。（4）鐘花櫻桃（0.816）和高盆櫻桃（0.799）的AUC值均高于隨機(jī)測(cè)試（0.500），說(shuō)明BIOCLIM模型可以準(zhǔn)確預(yù)測(cè)高盆櫻桃及鐘花櫻桃的分布。這為高盆櫻桃和鐘花櫻桃資源保護(hù)、物種鑒定和譜系地理學(xué)的研究提供重要指導(dǎo)。

關(guān)鍵詞： 高盆櫻桃， 鐘花櫻桃， 地理分布， 氣候限制因子， BIOCLIM模型

Modeling geographical distribution pattern and comparison of ecological characteristics between Cerasus cerasoides and C. campanulata

ZHU Shuxia， ZHU Hong， Cheng Lin， YI Xiangui， WANG Xianrong*

（ Collaborative Innovation Center of the Southern Modern Forestry， College of Forest Resources and Environment， Nanjing Forestry University， Nanjing 210037， China ）

Abstract：BIOCLIM model was used to simulate the range of modern suitable areas by collecting the geographical coordinates of its actual distribution and predict the change of their potential distribution under climate change （CCM3） for 2 100 s. The dominant climatic factors were determined by principal component analysis and correlation analysis. The differences between two species in climate limiting factors were analyzed. In addition， the receiver prediction effect was evaluated using the receiver operating characteristic （ROC）. The results were as follows： （1） The actual distributions of C. cerasoides and C. campanulata covered most areas of the South of the Yangtze River in China. The contemporary distribution centers were located in the Yunnan-Guizhou Plateau， the Wuyi Mountain Range and the Nanling Mountains. （2） Under the future climate change situation （CCM3）， the range of potential suitable areas of C. cerasoides and C. campanulata would decrease. The suitable habitat in Southwest （C. cerasoides） and Southeast China （C. campanulata） would be significantly reduced， while that of C. campanulata in West Hunan would increase. （3） Principal component analysis （PCA） and correlation analysis indicated that the annual precipitation （bio12）， precipitation of the driest quarter （bio16）， precipitation of the warmest quarter （bio18） and temperature seasonality （bio4） were the major factors influencing the geographical distribution of C. cerasoides and C. campanulata. “Radiator variability” was the most important environmental factor causing the difference in the distribution of two species. （4） Both of the AUC values of C. campanulata （0.816） and C. cerasoides （0.799） were higher than that of random test （0.500）， which indicated that the BIOCLIM model could accurately predict the distribution of C. cerasoides and C. campanulata. Our results will be an important guide for the study of conservation， identification and phytogeography of C. cerasoides and C. campanulata.

Key words： Cerasus cerasoides， Cerasus campanulata， geographical distribution， climatic limited factor， BIOCLIM model

高盆櫻桃（Cerasus cerasoides）與鐘花櫻桃（C. campanulata）為櫻屬（Cerasus）中特有的具玫紅花色的亞熱帶樹(shù)種，且花期早、抗病蟲(chóng)害能力強(qiáng)。其優(yōu)良的性狀大大拓寬了觀賞櫻花的花期和花色。其中，高盆櫻桃又名云南冬櫻花，多生于溝谷密林中，海拔為1 300～2 200 m，是Don David于1825年發(fā)現(xiàn)的新種;鐘花櫻桃又名福建山櫻花，多生于山谷林中及林緣，海拔為100～600 m，是Maximowiczi于1883年發(fā)表的新種（俞德浚和李朝鑾，1984，1986）。然而，自?xún)蓚€(gè)物種發(fā)現(xiàn)以來(lái)，其分類(lèi)地位一直存在較大爭(zhēng)議，如1910年Koidzumi將鐘花櫻桃作為高盆櫻桃的變種，并將其命名為Prunus cerasoides var. campanulata;王賢榮和向其柏（1998）發(fā)表櫻屬植物分類(lèi)訂正，認(rèn)為鐘花櫻桃應(yīng)為高盆櫻桃分布區(qū)偏東的變種;2003年Flora of China中將高盆櫻桃和鐘花櫻桃定為兩個(gè)種。盡管高盆櫻桃沿海拔和緯度與鐘花櫻桃間存在著明顯的生境分化，但決定其生態(tài)位分化的環(huán)境因子目前并不清楚，鑒別和分析兩者的地理分布格局、生態(tài)因子差異有助于闡明其表型差異以及物種形成的原因（Mallet，2005）。

植被是全球氣候變化的敏感指示器，氣候條件是決定植被物候、生產(chǎn)力、分布格局等特征及其動(dòng)態(tài)變化的主要非生物因素（焦珂?zhèn)サ龋?018）。植物與氣候間的相互關(guān)系以及地理分布格局的模擬和預(yù)測(cè)，是基礎(chǔ)生態(tài)學(xué)和生物地理學(xué)研究的重要工具（Elith et al.，2006）。隨著計(jì)算機(jī)科學(xué)以及地理科學(xué)的不斷發(fā)展，生態(tài)位模型在植物地理分布格局的模擬與預(yù)測(cè)方面得到了廣泛的應(yīng)用，常用的生態(tài)位模型有MaxEnt、GARP、BIOCLIM和DOMAIN等。不同的生態(tài)位模型由于其算法和原理不同，具有各自的優(yōu)缺點(diǎn)（李雪梅等，2007;許仲林等，2015）。BIOCLIM作為物種適生區(qū)預(yù)測(cè)的經(jīng)典模型之一，已經(jīng)在生物多樣性保護(hù)（張藍(lán)月等，2012;王剛等，2013;姜建福等，2014）、氣候變化對(duì)物種分布的影響（李垚等，2014;劉芙蓉等，2017;邱靖等，2018）、外來(lái)入侵植物風(fēng)險(xiǎn)評(píng)估（邁迪娜·吐?tīng)栠d，2017）等方面取得了較好的研究效果。即使在較小的樣本容量時(shí)，BIOCLIM模型也能達(dá)到一個(gè)較為滿(mǎn)意的預(yù)測(cè)結(jié)果（邵慧等，2009）。

BIOCLIM模型預(yù)測(cè)主要基于DIVA-GIS軟件，此軟件操作簡(jiǎn)單，通用性強(qiáng)，結(jié)果直觀（Hijmans et al.，2001;謝春平，2011;朱弘等，2017）。此外，DIVA-GIS網(wǎng)站免費(fèi)提供一整套不同分辨率的全球氣候數(shù)據(jù)資料，以及不同大氣環(huán)流模型產(chǎn)生的氣候數(shù)據(jù)包?；贒IVA-GIS可以快速獲取實(shí)際地理分布點(diǎn)的氣候數(shù)據(jù)并提取未來(lái)氣候信息的預(yù)測(cè)值（田聰?shù)龋?015）。筆者通過(guò)DIVA-GIS軟件獲取高盆櫻桃與鐘花櫻桃氣候數(shù)據(jù)及其實(shí)際地理分布圖，利用BIOCLIM模型預(yù)測(cè)當(dāng)下及其未來(lái)的適生范圍，并進(jìn)行當(dāng)下多度分析。通過(guò)比較兩者在區(qū)域尺度下的氣候響應(yīng)規(guī)律，旨在明確兩者的地理分布格局和生態(tài)特征方面的差異，以期為高盆櫻桃及鐘花櫻桃物種鑒定、分子譜系地理學(xué)分析、野生資源的有效保護(hù)和合理利用提供重要的借鑒。

1?材料與方法

1.1 物種分布數(shù)據(jù)收集及地理分布圖繪制

（1）根據(jù)國(guó)家標(biāo)本資源共享平臺(tái)（http：//www.nsii.org.cn/）、中國(guó)植物圖像庫(kù)（PPBC;http：//www.plantphoto.cn/）、中國(guó)自然標(biāo)本館（CFH;http：//www.cfh.ac.cn/）以及已發(fā)表的相關(guān)文獻(xiàn)資料（陳劍英等，1999;呂月良，2006;蘇倩，2007）獲取分布點(diǎn)信息。刪除錯(cuò)誤及模糊記錄，剔除培育和人工引種的標(biāo)本記錄，最終共獲得186個(gè)精確分布點(diǎn)的信息。其中高盆櫻桃82條記錄，鐘花櫻桃104條記錄。（2）利用百度拾取坐標(biāo)系（http：//api.map.baidu.com/lbsapi/getpoint/index.html）及Google Earth獲取每個(gè)分布點(diǎn)的經(jīng)緯度，錄入到Excel 2013，并轉(zhuǎn)化為文本格式。（3）基于DIVA-GIS（Version7.5）軟件包（http：//www.diva-gis.org），以中國(guó)省級(jí)行政區(qū)劃矢量地圖（下載自國(guó)家基礎(chǔ)地理信息系統(tǒng)網(wǎng)站http：//nfgis.nsdi.gov.cn/）為底圖，導(dǎo)入高盆櫻桃及鐘花櫻桃地理分布數(shù)據(jù)，繪制實(shí)際地理分布圖（圖1）。

1.2 模型建立及氣候數(shù)據(jù)分析

（1）DIVA-GIS網(wǎng)站的全球氣候數(shù)據(jù)資料來(lái)源于WorldClim（http：//www.worldclim.org/），當(dāng)前和未來(lái)CCM3背景下2100年的氣候數(shù)據(jù)分辨率均選用2.5′，將其導(dǎo)入DIVA-GIS軟件，獲取實(shí)際有效分布點(diǎn)所對(duì)應(yīng)的19個(gè)氣候變量數(shù)據(jù)。（2）采用DIVA-GIS耦合的BIOCLIM模型，完成高盆櫻桃與鐘花櫻桃當(dāng)前及其未來(lái)適生區(qū)模擬分布預(yù)測(cè)及適生等級(jí)劃分和當(dāng)下多度分析。（3）首先，利用PAST3（Version 3.14）軟件包（http：//folk.uio.no/ohammer/past/）對(duì)所獲取的氣候變量數(shù)據(jù)進(jìn)行主成分分析（principal component analysis，PCA），通過(guò)研究不同氣候因子在各分布點(diǎn)間的變異大小，選出方差貢獻(xiàn)率較大的候選主要?dú)夂蛞蜃?然后，通過(guò)相關(guān)性分析，僅保留具有低相關(guān)性（r<0.8）的氣候變量作為限制高盆櫻桃和鐘花櫻桃當(dāng)下地理分布的主要?dú)夂蛞蜃?。此外，基于得到的主?dǎo)氣候因子繪制高盆櫻桃和鐘花櫻桃的PCA分布圖，用PAST3繪制兩者的頻率分布直方圖并進(jìn)行正態(tài)曲線(xiàn)擬合，比較兩者在氣候限制因子方面的差異。

1.3 模型精度檢驗(yàn)

受試者工作特征曲線(xiàn)（receiver operation cha-racteristic， ROC）是模型性能評(píng)價(jià)常用的方法。ROC曲線(xiàn)是以真陽(yáng)性率（靈敏度）為縱坐標(biāo)，假陽(yáng)性率（1-特異度）為橫坐標(biāo)所形成的曲線(xiàn)。曲線(xiàn)下的面積即為AUC（areas under curve）值，AUC值通過(guò)曲線(xiàn)上相鄰2點(diǎn)以直線(xiàn)相連與橫軸圍成的各小梯形面積相加得到，即AUC=S1+S2+…+Sn。AUC值越大表示與隨機(jī)分布相距越遠(yuǎn)，環(huán)境變量與預(yù)測(cè)的物種地理分布模型之間相關(guān)性越大，即模型預(yù)測(cè)效果越好（崔麟和魏洪義，2016）。

2?結(jié)果與分析

2.1 地理分布格局

本研究將BIOCLIM模型潛在分布區(qū)預(yù)測(cè)的結(jié)果劃分為非適生區(qū)（Unsuitable distribution range）（0）、低度適生區(qū)（Low distribution range）（0～2.5%）、中度適生區(qū)（Medium distribution range）（2.5%～5%）、高度適生區(qū)（High distribution range）（5%～10%）、極度適生區(qū)（Extremely high distribution range）（10%～20%）和最適生區(qū)（Optimum distribution）（20%～37%）6個(gè)等級(jí)（圖2）。根據(jù)Flora of China的描述，高盆櫻桃主要分布于西藏南部、云南西北部;鐘花櫻桃主要分布于福建、廣東、廣西、海南、湖南、臺(tái)灣、浙江。實(shí)際分布點(diǎn)及預(yù)測(cè)結(jié)果顯示（圖2：A，B），高盆櫻桃和鐘花櫻桃主要分布于亞熱帶地區(qū)，幾乎覆蓋長(zhǎng)江以南各省區(qū)（圖2：A）。

高盆櫻桃主要分布在94.37°—105.83° E、21.58°—29.15° N的地區(qū)，除覆蓋我國(guó)西藏東南部和云南大部分地區(qū)外，在四川南部也有分布。鐘花櫻桃的分布區(qū)為105.78°—121.55° E、21.91°—30.32° N，除Flora of China記錄的分布區(qū)外，江西也是其分布區(qū)，而在海南無(wú)分布點(diǎn)。此外，湖北西南部和東南部、貴州東南部、安徽南部為鐘花櫻桃的潛在分布區(qū)。多度分析被用于分布點(diǎn)的分級(jí)評(píng)估，以分布地點(diǎn)為參數(shù)，每個(gè)地點(diǎn)代表1個(gè)植物群落的存在，在1°×1°的網(wǎng)格內(nèi)統(tǒng)計(jì)的變種豐富度。將高盆櫻桃與鐘花櫻桃集中分布點(diǎn)劃分為6個(gè)等級(jí)，編號(hào)為0、1-1、2-2、3-3、4-4、5-5表示高盆櫻桃與鐘花櫻桃分布程度越來(lái)越集中。綜合當(dāng)下適生區(qū)預(yù)測(cè)表明，高盆櫻桃的現(xiàn)代分布中心位于云貴高原，鐘花櫻桃的現(xiàn)代分布中心位于武夷山脈和南嶺山脈（圖2：C）。

未來(lái)氣候變化情景下，對(duì)比圖2：A和圖2：B可知，高盆櫻桃及鐘花櫻桃的整體分布格局較為相似，而適生區(qū)概率顏色由深到淺，預(yù)測(cè)了未來(lái)高盆櫻桃及鐘花櫻桃的適宜面積縮小、生境破碎。在云南中部，高盆櫻桃的適生范圍將明顯縮小。福建中部、浙江西南部、江西西南部和東南部、廣東北部，鐘花櫻桃原有的適生范圍可能大幅減小。

然而在CCM3情境下，湖南西部可能新增較大面積的鐘花櫻桃適合區(qū)。

2.2 主導(dǎo)氣候因子分析

因19個(gè)氣候因子間可能存在一定相關(guān)性，現(xiàn)通過(guò)PCA和相關(guān)性分析，篩選出相關(guān)性較低且在各分布點(diǎn)間變異豐富的主導(dǎo)氣候因子（劉然等，2018）。最終顯示，年降水量（bio12）、最濕季降水量（bio16）、最暖季降水量（bio18）、溫度季節(jié)變化方差（bio4）是影響高盆櫻桃及鐘花櫻桃當(dāng)下適生區(qū)的主導(dǎo)氣候限制因子，4個(gè)主導(dǎo)氣候因子在高盆櫻桃（鐘花櫻桃）分布地點(diǎn)的波動(dòng)范圍分別是660～2 019（1 223～2 233）mm、380～1 173（578～1 119）mm、380～1 173（446～1 119）mm和326～565（488～924）mm（標(biāo)準(zhǔn)差 × 100）。隨后利用得到的4個(gè)主導(dǎo)氣候因子進(jìn)行的PCA結(jié)果顯示，主成分1和2（PC1和PC2）可解釋96.61%的變異?；赑C1和PC2對(duì)兩物種各分布點(diǎn)進(jìn)行分類(lèi)分析，圖3結(jié)果顯示，高盆櫻桃和鐘花櫻桃各分布點(diǎn)總體分為兩類(lèi)，一類(lèi)主要包括高盆櫻桃的分布點(diǎn)，另一類(lèi)則主要包括鐘花櫻桃的分布點(diǎn)，表明兩個(gè)物種現(xiàn)存分布范圍的生境特點(diǎn)存在明顯差異。

對(duì)兩者分布的限制環(huán)境因子的頻率直方圖進(jìn)行比較分析及其正態(tài)分布進(jìn)行擬合（圖4）。圖4結(jié)果表明，高盆櫻桃的最適年降水量（bio12）范圍為900～1 050 mm，鐘花櫻桃最適年降水量（bio12）范圍為1 500～1 650 mm;最適合的最濕季年降水量（bio16），高盆櫻桃為800～900 mm，鐘花櫻桃為700～800 mm;最適合的最暖季降水量（bio18），高盆櫻桃為600～700 mm，鐘花櫻桃為500～600 mm;最適合的溫度季節(jié)變化方差（bio4），高盆櫻桃為390～420（標(biāo)準(zhǔn)差 × 100），鐘花櫻桃為660～720（標(biāo)準(zhǔn)差 × 100）。這反映了“水分因子”和“熱量變異幅度”對(duì)兩者的影響，“熱量變異幅度”是高盆櫻桃和鐘花櫻桃分布存在差異的最主要的環(huán)境因子（圖4）。

2.3 模型檢驗(yàn)

從圖5可以看出，ROC曲線(xiàn)對(duì)高盆櫻桃及鐘花櫻桃的預(yù)測(cè)結(jié)果均很好。高盆櫻桃AUC的值達(dá)0.799，鐘花櫻桃的AUC值高達(dá)0.816，均顯著高于隨機(jī)測(cè)試AUC（0.500），說(shuō)明BIOCLIM模型的預(yù)測(cè)效果良好，可以準(zhǔn)確預(yù)測(cè)高盆櫻桃及鐘花櫻桃的分布。

3?討論

3.1 生態(tài)因子差異比較及其物種鑒定

王賢榮和向其柏（1998）根據(jù)觀察大量高盆櫻桃和鐘花櫻桃的標(biāo)本發(fā)現(xiàn)，兩者的主要區(qū)別僅在于萼片形狀和開(kāi)花習(xí)性的差別，而形態(tài)的差別可能與局部的小氣候和海拔有關(guān)，如陳璋（2008）通過(guò)研究鐘花櫻桃的表型多樣性發(fā)現(xiàn)，花表型性狀與經(jīng)緯度、海拔、年降水量呈顯著或極顯著相關(guān)，開(kāi)花習(xí)性與年日照時(shí)數(shù)顯著相關(guān);李蒙等（2014）通過(guò)研究山櫻花地理分布與氣候條件因子的關(guān)系時(shí)發(fā)現(xiàn)，山櫻花的分布受溫度影響較大。本研究主導(dǎo)氣候因子分析結(jié)果表明，高盆櫻桃與鐘花櫻桃的平均溫度季節(jié)變化方差（bio4）相差285（標(biāo)準(zhǔn)差 × 100）， 是兩者分布存在差異的最主要的環(huán)境因子。高盆櫻桃位于105° E以西的寒冷干燥的中國(guó)-喜馬拉雅亞區(qū)，主要受印度季風(fēng)影響;鐘花櫻桃位于105° E以東的溫暖濕潤(rùn)的中國(guó)-日本亞區(qū)，主要受太平洋季風(fēng)影響。自西向東具有海拔的梯度性和不同的季風(fēng)環(huán)境， 不同區(qū)域間的異質(zhì)環(huán)境很有可能形成植物的適應(yīng)性分化（葉俊偉等，2017a）。因此，推測(cè)兩者的形態(tài)差別（如萼片性

狀、開(kāi)花習(xí)性等）極有可能是由兩地的溫度差異造成的。

王文采（1992a，b）認(rèn)為中國(guó)亞熱帶地區(qū)植物可能起源于西南地區(qū)，然后向東擴(kuò)散，擴(kuò)散的路線(xiàn)即為北部的秦嶺-大巴山脈、中部的大婁山-武陵山山脈和南部的南嶺山脈。但是這個(gè)假說(shuō)只得到較少的驗(yàn)證。李苗苗（2009）根據(jù)單倍型的地理分布，推測(cè)櫻亞屬的分布中心在第四紀(jì)冰川來(lái)臨之前已經(jīng)形成，且還經(jīng)歷了云貴高原-華中山地-華東江南丘陵的遷移。高盆櫻桃和鐘花櫻桃的現(xiàn)代分布中心云貴高原、武夷山脈和南嶺是冰期時(shí)重要的南方避難所。在冰期和間冰期時(shí)，持續(xù)的長(zhǎng)期隔離使植物保持了較高的種群間遺傳多樣性和遺傳分化（葉俊偉等，2017b）。本研究揭示出兩個(gè)物種分布點(diǎn)環(huán)境因子存在較大差異，且自然分布范圍存在明顯區(qū)別，推測(cè)二者在進(jìn)化過(guò)程中由于氣候等因素產(chǎn)生了明顯分化。但是，高盆櫻桃和鐘花櫻桃是不是一個(gè)物種，還是種下的變種;他們具體的遷移路線(xiàn)是否符合以上假說(shuō);鐘花櫻桃是高盆櫻桃向東遷移之后，為適應(yīng)當(dāng)?shù)氐臍夂驐l件而形成的物種分化嗎？為解釋以上問(wèn)題，后續(xù)可結(jié)合化石、孢粉、親緣地理學(xué)等領(lǐng)域的研究結(jié)果。

3.2 保護(hù)與利用建議

高盆櫻桃與鐘花櫻桃作為典型的亞熱帶樹(shù)種，在秦嶺-淮河以北表現(xiàn)很差，不開(kāi)花或者晚花，只有個(gè)別經(jīng)改良后的鐘花櫻桃品種如陽(yáng)光Cerasus ‘Youkou’和河津櫻C. × kanzakura ‘Kawazu-zakura’等表現(xiàn)良好。國(guó)產(chǎn)櫻花種類(lèi)繁多，但大部分處于未開(kāi)發(fā)的狀態(tài)，品種培育技術(shù)遠(yuǎn)遠(yuǎn)落后于日本。據(jù)記載，目前市場(chǎng)上所應(yīng)用的鐘花櫻桃種系下的品種有13種之多，但基本是日本學(xué)者所培育（王賢榮，2014;王青華等，2014）。高盆櫻桃的品種在云南有所記載，但基本是野生種自然變異所選育的，而且僅限于西南地區(qū)推廣使用。高盆櫻桃及鐘花櫻桃多生長(zhǎng)在林緣或受到較大干擾的次生林中，野生資源破壞嚴(yán)重（段曉梅等，2004;王賢榮等，2006）。因此，對(duì)高盆櫻桃和鐘花櫻桃的資源保護(hù)與品種合理開(kāi)發(fā)利用提出如下建議：（1）通過(guò)天然選優(yōu)、人工雜交、化學(xué)誘變等手段加大新品種選育規(guī)模，培育出更具觀賞價(jià)值和國(guó)人代表性的櫻花新品種。（2）受全球氣候變暖的影響，分布區(qū)處于南緣的幾個(gè)省份的潛在分布區(qū)面積將會(huì)縮小，而在更高緯度的湖南西部會(huì)新增分布區(qū)，總體反映了分布區(qū)北擴(kuò)的趨勢(shì)。這表明高濃度的溫室氣體排放將嚴(yán)重影響高盆櫻桃及鐘花櫻桃的生境，造成生境破碎化、生物多樣性降低等一系列嚴(yán)重后果。因此，為保護(hù)兩者的生境面積，社會(huì)各界應(yīng)采取植樹(shù)造林、節(jié)能減排、開(kāi)發(fā)新能源等有效措施嚴(yán)格控制溫室氣體排放。（3）為加強(qiáng)對(duì)高盆櫻桃及鐘花櫻桃現(xiàn)有資源的保護(hù)力度，應(yīng)將福建中部、浙江西南部、江西西南部和東南部、廣東北部等地的保護(hù)區(qū)作為兩者的重點(diǎn)保護(hù)單元，并將其潛在分布區(qū)納入到新的保護(hù)工作中。通過(guò)移栽、嫁接等手段促進(jìn)物種的再分布。

參考文獻(xiàn)：

CHEN JY， DONG Q， YANG CB， 1999. A study on ecogeo-graphy of Cerasus cerasoides [J]. J SW For Coll， 19（1）： 2-8.[陳劍英， 董瓊， 楊超本， 1999. 冬櫻花生態(tài)地理學(xué)研究 [J]. 西南林學(xué)院學(xué)報(bào)， 19（1）：2-8.]

CHEN Z， 2008. Phenotypic variation in natural populations of Cerasus campanulata Maxim. [J]. J Trop Subtrop Bot， 16（1）： 61-68. [陳璋， 2008. 福建山櫻花天然居群表型變異研究 [J]. 熱帶亞熱帶植物學(xué)報(bào)， 16（1）： 61-68.]

CUI L， WEI HY， 2016.Predicting the geographical distribution of Urochela distincta distant with MaxEnt and DIVA-GIS [J].J Plant Prot， 43（3）： 362-368. [崔麟， 魏洪義， 2016. 基于MaxEnt和DIVA-GIS的亮壯異蝽潛在地理分布預(yù)測(cè) [J]. 植物保護(hù)學(xué)報(bào)， 43（3）：362-368.]

DUAN XM， WANG ZH， FAN GS， 2004. Community characte-ristics and phytogeographical meaning of Cerasus cerasoides in Yunnan Province [J]. J Nanjing For Univ （Nat Sci Ed）， 28（6）： 83-86. [段曉梅， 王自輝， 樊國(guó)盛， 2004. 云南冬櫻花的群落特征及地理分布 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 28（6）： 83-86.]

ELITH J， GRAHAM CH， ANDERSON RP， et al.， 2006.Novel methods improve prediction of species distributions from occurrence data [J].Ecography， 29（2）： 129-151.

Flora of China Editorial Committee， 2003.Flora of China （9）[M].Beijing： Science Press; St.Louis： Missouri Botanical Garden Press.

HACK JJ， KIEHL JT， HURRELL JW， 1998.The hydrologic and thermodynamic characteristics of the NCAR CCM3 [J].J Climate， 11（6）： 1179-1206.

HIJMANS RJ， GUARINO L， CRUZ M， et al.， 2001.Computer tools for spatial analysis of plant genetic resources data： 1.DIVA-GIS [J].Plant Genet Resour Newslett， 127： 15-19.

JIAO KW， GAO JB， WU SH， et al.， 2018.Research progress on the response processes of vegetation activity to climate change [J].Acta Ecol Sin， 38（6）： 2229-2238. [焦珂?zhèn)ィ?高江波， 吳紹洪， 等， 2018. 植被活動(dòng)對(duì)氣候變化的響應(yīng)過(guò)程研究進(jìn)展 [J]. 生態(tài)學(xué)報(bào)， 38（6）：2229-2238.]

JIANG JF， FAN XC， ZHANG Y， et al.， 2014. Modeling the geographic distribution of three endangered Vitis species in China [J].Chin J Ecol， 33（6）： 1615-1622. [姜建福， 樊秀彩， 張穎， 等， 2014. 中國(guó)三種瀕危葡萄屬（Vitis L.）植物的地理分布模擬 [J]. 生態(tài)學(xué)雜志， 33（6）：1615-1622.]

LI M， YI XG， WANG HC， et al.， 2014.Studies on the relationship between Cerasus serrulata distribution region and the environmental factors[J].J Nanjing For Univ （Nat Sci Ed）， 38（S）： 74-80. [李蒙， 伊賢貴， 王華辰， 等， 2014. 山櫻花地理分布與水熱環(huán)境因子的關(guān)系 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 38（S）：74–80.]

LI MM， 2009.Phylogeography of Prunus subgenus Cerasus and genetic diversity of P.pseudocerasus[D].Xian： Northwest University： 1-57. [李苗苗， 2009. 櫻亞屬植物分子親緣地理及中國(guó)櫻桃自然居群遺傳多樣性研究 [D].西安：西北大學(xué)：1-57.]

LI XM， CHENG XQ， 2007.Development of niche theory and its application in each fields of ecology [J].J Beijing For Univ， 29（Suppl. 2）： 294-298. [李雪梅， 程小琴， 2007. 生態(tài)位理論的發(fā)展及其在生態(tài)學(xué)各領(lǐng)域中的應(yīng)用 [J]. 北京林業(yè)大學(xué)報(bào)， 29（增刊2）：294-298.]

LIU FR， LUO JX， YANG MJ， 2017.Potential geographic distribution areas of Toona sinensis and their climatic characte-ristics in Southwest China regions [J].Sichuan For Sci Technol， 38（1）：16-21. [劉芙蓉， 羅建勛， 楊馬進(jìn)， 2017. 西南地區(qū)香椿的潛在地理分布區(qū)及其氣候特點(diǎn) [J]. 四川林業(yè)科技， 38（1）：16-21.]

LIU R， WANG CJ， HE J， et al.， 2018.Analysis of geographical distribution of Abies in China under climate change [J].Bull Bot Res， 38（1）： 37-46. [劉然， 王春晶， 何健， 等， 2018. 氣候變化背景下中國(guó)冷杉屬植物地理分布模擬分析 [J]. 植物研究， 38（1）：37-46.]

L YL， 2006.Genetic diversity and breeding strategy of Cerasus campanulata maxim [D].Nanjing： Nanjing Forest University： 1-119. [呂月良， 2006. 福建山櫻花群體遺傳多樣性、繁育技術(shù)體系和育種策略研究 [D]. 南京：南京林業(yè)大學(xué)：1-119.]

MALLET J， 2005.Hybridization as an invasion of the genome [J].Trends Ecol Evol， 20（5）： 229-237.

QIU J， ZHU H， CHEN X， et al.， 2018.Modeling the suitable areas and ecological characteristics of Sorbus alnifolia using DIVA-GIS software [J].J Beijing For Univ， 40（9）： 25-32. [邱靖， 朱弘， 陳昕， 等， 2018. 基于DIVA-GIS的水榆花楸適生區(qū)模擬及生態(tài)特征 [J]. 北京林業(yè)大學(xué)學(xué)報(bào)， 40（9）：25-32.]

SHAO H， TIAN JQ， GUO K， et al.， 2009.Effects of sample size and species traits on performance of BIOCLIM in predicting geographical distribution of tree species： A case study with 12 deciduous Quercus species indigenous to China[J].Chin J Plant Ecol， 33（5）： 870-877. [邵慧， 田佳倩， 郭柯， 等， 2009. 樣本容量和物種特征對(duì)BIOCLIM模型模擬物種分布準(zhǔn)確度的影響——以12個(gè)中國(guó)特有落葉櫟樹(shù)種為例 [J]. 植物生態(tài)學(xué)報(bào)， 33（5）：870-877.]

SU Q， 2007.Genetic diversity among populations of Cerasus campanulata[D].Nanjing： Nanjing Forest University： 1-65. [蘇倩， 2007. 福建山櫻花種源遺傳多樣性研究 [D]. 南京：南京林業(yè)大學(xué)：1-65.]

TIAN C， MU NR， ZHU ZY， et al.， 2015.Method for the rapid obtaining of climate data based on DIVA-GIS [J].J Agric， 5（5）： 109-113. [田聰， 穆尼熱， 朱忠艷， 等， 2015. 基于DIVA-GIS快速獲取氣候信息的方法 [J]. 農(nóng)學(xué)學(xué)報(bào)， 5（5）：109-113.]

WANG G， MENG YQ， SUN J， et al.， 2013.Modeling the geographic distribution of Rhododendron maculiferum [J].J NW Agric For Univ （Nat Sci Ed）， 41（5）： 173-177. [王剛， 孟艷瓊， 孫婧， 等， 2013. 麻花杜鵑的地理分布模擬 [J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(xué)版）， 41（5）：173-177.]

WANG QT， LIU XH， XU L， 2014.The illustration of main cherry cultivars in China [M].Hangzhou： Zhejiang Science and Technology Press. [王青華， 柳新紅， 徐梁， 2014. 中國(guó)主要栽培櫻花品種圖鑒 [M]. 杭州：浙江科學(xué)技術(shù)出版社]

WANG WT， 1992a.On some distribution patterns and some migration routes found in the eastern Asiatic region [J].J Syst Evol， 30（1）： 1-24. [王文采， 1992a. 東亞植物區(qū)系的一些分布式樣和遷移路線(xiàn) [J]. 植物分類(lèi)學(xué)報(bào)， 30（1）： 1-24.]

WANG WT， 1992b.On some distribution patterns and some migration routes found in the eastern Asiatic region （Cont. ） [J].J Syst Evol， 30（2）： 97-117. [王文采， 1992b. 東亞植物區(qū)系的一些分布式樣和遷移路線(xiàn)（續(xù)） [J]. 植物分類(lèi)學(xué)報(bào)， 30（2）： 97-117.]

WANG XR， 2014.An illustrated monograph of cherry cultivars in China[M].Beijing： Science Press. [王賢榮， 2014. 中國(guó)櫻花品種圖志 [M]. 北京：科學(xué)出版社]

WANG XR， YAN DL， YI XG， 2006. Characteristics of the Cerasus cerasoides var. campanulata community structure and its species diversity [J]. J SW For Coll， 26（4）： 9-13. [王賢榮， 閆道良， 伊賢貴， 2006. 鐘花櫻群落結(jié)構(gòu)特征與物種多樣性研究 [J]. 西南林學(xué)院學(xué)報(bào)， 26（4）： 9-13.]

WANG XR， XIANG QB， 1998.Notes on some species of the genus Cerasus MIQ [J].J Nanjing For Univ， 22（4）：62-64. [王賢榮， 向其柏， 1998. 櫻屬植物分類(lèi)訂正 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)， 22（4）， 62-64.]

XIE CP， 2011.Simple biogeography distribution picture draw based on DIVA-GIS [J].Hubei Agric Sci， 50（11）： 2345-2348.[謝春平， 2011. 基于DIVA-GIS生物地理分布圖的繪制 [J]. 湖北農(nóng)業(yè)科學(xué)， 50（11）：2345-2348.]

XU ZL， PENG HH， PENG SZ，2015.The development and evaluation of species distribution models [J].Acta Ecol Sin， 35（2）： 557-567. [許仲林， 彭煥華， 彭守璋， 2015. 物種分布模型的發(fā)展及評(píng)價(jià)方法 [J]. 生態(tài)學(xué)報(bào)， 35（2）：557-567.]

YE JW， ZHANG Y， WANG XJ， 2017.Phylogeographic breaks and the mechanisms of their formation in the Sino-Japanese floristic region [J].Chin J Plant Ecol， 41（9）： 1003-1019.[葉俊偉， 張陽(yáng)， 王曉娟， 2017. 中國(guó)亞熱帶地闊葉林植物的譜系地理歷史 [J]. 植物生態(tài)學(xué)報(bào)， 41（9）：1003-1019.]

YE JW， ZHANG Y， WANG XJ， 2017.Phylogeographic history of broad-leaved forest plants in subtropical China [J].Acta Ecol Sin， 37（17）： 5894-590. [葉俊偉， 張陽(yáng)， 王曉娟， 2017. 中國(guó)亞熱帶地區(qū)闊葉林植物的譜系地理歷史 [J]. 生態(tài)學(xué)報(bào)， 37（17）：5894-5904.]

YU DJ， LI CL， 1984.Taxa nova Rosacearum sinicarum（V）[J].J Syst Evol， 23（3）：210. [俞德浚， 李朝鑾， 1984. 中國(guó)薔薇科植物分類(lèi)之研究（五） [J]. 植物分類(lèi)學(xué)報(bào)， 23（3）：210.]

YU DJ， LI CL， 1986.Flora Reipublicae Popularis Sinicae [M].Beijing： Science Press， （38）： 46-89. [俞德浚， 李朝鑾， 1986. 中國(guó)植物志 [M]. 北京：科學(xué)出版， （38）：46-89.]

ZHANG LY， YE XB， LIU N， et al.，2012.Research on geographic distribution and potential distribution of Mucuna birdwoodiana[J].Guihaia， 32（1）： 27-32. [張藍(lán)月， 葉向斌， 劉念， 等， 2012. 白花油麻藤的地理分布及適生區(qū)預(yù)測(cè) [J]. 廣西植物， 32（1）：27-32.]

ZHU H， YOU LX， LI YF， et al.， 2017. Modeling the geogra-phical distribution pattern and climatic limited factors of Cerasus schneideriana [J].J Trop Subtrop Bot， 25（4）： 315-322. [朱弘， 尤祿祥， 李涌福，

等， 2017. 浙閩櫻桃地理分布模擬及氣候限制因子分析 [J]. 熱帶亞熱帶植物學(xué)報(bào)， 25（4）：315-322.]