基于Maxent模型的檵木分布格局模擬

孟 影,馬姜明,*,王永琪,莫燕華

1 廣西師范大學(xué)可持續(xù)發(fā)展創(chuàng)新研究院, 桂林 5410062 廣西師范大學(xué)珍稀瀕危動(dòng)植物生態(tài)與環(huán)境保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室, 桂林 5410063 廣西師范大學(xué)生命科學(xué)學(xué)院, 桂林 541006

理解和預(yù)測(cè)物種如何應(yīng)對(duì)全球氣候變化是生態(tài)學(xué)研究的核心問(wèn)題之一[1],基于Maxent模型可以檢測(cè)物種潛在分布區(qū)對(duì)于歷史上氣候變化的應(yīng)對(duì)措施,從而預(yù)測(cè)未來(lái)全球氣候變暖對(duì)于不同物種可能造成的潛在影響。Maxent模型應(yīng)用于預(yù)測(cè)保護(hù)物種、入侵種和藥用植物的潛在分布格局,不但具有良好的模型預(yù)測(cè)效果,而且具有較好的穩(wěn)定性[2]。模型推出后廣泛應(yīng)用于竹葉花椒(Zanthoxylumarmatum)[3]、 紅松(Pinuskoraiensis)[4]、美味獼猴桃(Actinidiadeliciosa)[5]、假臭草(Praxelisclematidea)[6]、黃杉(Pseudotsugasinensis)林[7]、互花米草(Spartinaalterniflora)[8]、雙花木屬(Disanthus)[9]等物種的潛在分布研究,但是對(duì)于分布較廣的優(yōu)勢(shì)物種與環(huán)境之間的關(guān)系及其分布格局的研究較少。

氣候變化下物種的潛在分布格局發(fā)生變化,適生區(qū)也發(fā)生變遷。古氣候記錄顯示,金縷梅科(Hamamelidaceae)的化石證據(jù)發(fā)掘很豐富,最早的化石記錄普遍發(fā)現(xiàn)于晚白壟紀(jì)[10]。在末次盛冰期和末次間冰期氣候不斷旋回中,極大影響了現(xiàn)代植物的分布[11]。末次間冰期以來(lái),全球氣候波動(dòng)事件頻繁,強(qiáng)烈地影響了植物的歷史動(dòng)態(tài)變化,塑造了現(xiàn)代分布格局[12]。末次盛冰期以來(lái),全球氣候環(huán)境發(fā)生了從低溫—升溫—高溫的轉(zhuǎn)變,涵蓋了不同時(shí)間段的植被格局及演替特征[13]。例如：張愛(ài)平等準(zhǔn)確地模擬了3種云杉屬(Picea)植物末次間冰期、末次盛冰期、全新世中期和現(xiàn)代的適生分布區(qū)分布,推測(cè)了青扦(Piceawilsonii)和紫果云杉(Piceapurpurea)的生物避難所[12]。李璇等模擬白櫟(Quercusfabri)在末次盛冰期的結(jié)果顯示：高度適生區(qū)集中在華中地區(qū),隨氣候的變暖逐漸向北移動(dòng)[14]。Bertrand 等已經(jīng)證實(shí)了氣候變暖會(huì)影響植物遷徙向高海拔或高緯度地區(qū)[15]。

檵木(Loropetalumchinense)為金縷梅科檵木屬(Loropetalum)一種木本植物[16],廣泛分布于我國(guó)中部、南部及西南各省,喜生于向陽(yáng)的丘陵及山地,在廣西北部及中部,無(wú)論是“土山”生境還是喀斯特“石山”生境,在景觀上通常能形成以檵木為優(yōu)勢(shì)種的天然林斑塊,這尤其對(duì)于喀斯特石山生境的植被恢復(fù)具有非常重要的意義[17]。目前,有關(guān)檵木的研究主要為其功能成分研究[18]、藥理價(jià)值[19]、種群空間格局及關(guān)聯(lián)性[20],以及檵木群落的土壤生態(tài)化學(xué)計(jì)量學(xué)特征[21]、凋落葉的分解研究[22]、群落老齡林植物葉性狀[23]等,而其潛在分布及其與氣候的關(guān)系研究尚不清楚。因此,本研究將以檵木為研究對(duì)象,基于最大熵模型(Maxent)模擬當(dāng)代、末次間冰期、末次冰期等不同時(shí)期的潛在分布,研究歷史氣候變化對(duì)檵木分布格局的影響,揭示檵木對(duì)冰期和暖期的響應(yīng)變化,對(duì)于探討檵木發(fā)揮其生態(tài)優(yōu)勢(shì)提供現(xiàn)實(shí)意義。本研究試圖回答以下問(wèn)題：(1)氣候因子是如何影響檵木的地理分布,其分布的主導(dǎo)氣候因子是什么？(2)明確檵木的氣候適宜性,劃分各等級(jí)的潛在適宜級(jí)別。(3)不同時(shí)期檵木潛在分布區(qū)如何應(yīng)對(duì)氣候變化？這些問(wèn)題的解決不僅能夠?yàn)闄嬆咀鳛閮?yōu)勢(shì)樹(shù)種的選擇提供有力支撐,也為石漠化地區(qū)生態(tài)治理提供理論基礎(chǔ)。

1 數(shù)據(jù)來(lái)源及研究方法

1.1 樣品數(shù)據(jù)

檵木地理分布信息：檢索全球生物多樣性信息平臺(tái)(GBIF)和中國(guó)數(shù)字植物標(biāo)本館(http://www.cvh.org.cn/)等數(shù)據(jù)庫(kù)獲得檵木標(biāo)本的采集地信息,檵木最新分布點(diǎn)參考中國(guó)知網(wǎng)、讀秀、萬(wàn)方等國(guó)內(nèi)已公開(kāi)發(fā)表的論文數(shù)據(jù)。然后,對(duì)這些具體到縣的地理信息借助高德拾取器(https://lbs.amap.com/console/show/picker)來(lái)確定它的地理坐標(biāo),并轉(zhuǎn)化為WGS84地理坐標(biāo)[24]。本研究除去經(jīng)緯度重復(fù)和采集地不詳?shù)臉?biāo)本信息,共選取了227個(gè)檵木的地理位點(diǎn)。并用Excel編輯檵木現(xiàn)有分布的經(jīng)緯度信息,格式為.csv；用ArcGIS轉(zhuǎn)換圖層的格式并剪切需要的區(qū)域。

1.2 環(huán)境數(shù)據(jù)

從WorldClim數(shù)據(jù)庫(kù)(http://www.worldclim.org)下載,本文通過(guò) ArcGIS 軟件提取了19個(gè)歷史氣候數(shù)據(jù),其中包括11個(gè)溫度數(shù)據(jù)和8個(gè)降雨數(shù)據(jù)[25]。它們總體反映了溫度及降水的每年每季每月的平均值、最高溫、最低溫、最大降雨量、最小降雨量等數(shù)值(表1)。

表1 Maxent模型中的19個(gè)環(huán)境變量

1.3 地圖數(shù)據(jù)地圖數(shù)據(jù)

中國(guó)國(guó)界和省級(jí)行政區(qū)劃圖矢量底圖來(lái)源于國(guó)家基礎(chǔ)地理信息中心(http://ngcc.sbsm.gov.cn/)。

1.4 研究方法

A. 數(shù)據(jù)準(zhǔn)備: 用Excel編輯的物種現(xiàn)有分布的經(jīng)緯度信息,格式為. csv。

B. 打開(kāi)Maxent. bat, 將檵木的分布點(diǎn)數(shù)據(jù)和19個(gè)環(huán)境數(shù)據(jù)導(dǎo)入 MaxEnt 軟件,在右下角“Output directory”中設(shè)定輸出文件夾的位置,文件格式選擇. asc。在相關(guān)參數(shù)設(shè)置中,25%是用于驗(yàn)證模型,75%是軟件將隨機(jī)選取,最大迭代次數(shù)為500。為了使運(yùn)行結(jié)果精確性高,次數(shù)設(shè)置為10次。

C.點(diǎn)擊“Run”開(kāi)始運(yùn)行,運(yùn)算完成后會(huì)在已設(shè)定的輸出文件夾中生成多個(gè)文件。然后根據(jù)文件夾的數(shù)據(jù)結(jié)果進(jìn)行分析。

通過(guò)Maxent自身提供的模塊,即刀切法測(cè)定不同環(huán)境變量對(duì)檵木潛在適生區(qū)預(yù)測(cè)的貢獻(xiàn)。該方法分別計(jì)算模擬3個(gè)變量的訓(xùn)練得分情況。判定標(biāo)準(zhǔn)為：“僅此變量”時(shí)分值比例較高,說(shuō)明該因子的預(yù)測(cè)能力精確度較高,對(duì)物種分布貢獻(xiàn)較大,“除此變量”時(shí)的分值比例降低時(shí),說(shuō)明該變量具有較多的獨(dú)特信息,對(duì)物種分布較為重要。

Maxent模型采用ROC(receiver operating characteristic)曲線分析法對(duì)風(fēng)險(xiǎn)預(yù)測(cè)結(jié)果進(jìn)行精度檢驗(yàn)[26-27]。ROC曲線的下面積(AUC值)作為評(píng)價(jià)模型準(zhǔn)確性的指標(biāo),該值范圍為[0,1],越接近數(shù)值1表示與隨機(jī)分布距離越遠(yuǎn),獲得最大 AUC 值的方案作為推薦的預(yù)測(cè)研究方案。

2 結(jié)果與分析

2.1 全國(guó)檵木的分布圖

將227個(gè)檵木分布點(diǎn)導(dǎo)入ArcGIS生成的地理分布圖(圖1),可以看出檵木樣點(diǎn)分布東至浙江、西至四川、南至廣西、北至湖北。從標(biāo)本的分布密度來(lái)看,主要集中分布于中國(guó)的中部、南方以及西南省市等地區(qū),但這種分布與采集標(biāo)本的強(qiáng)度、頻率和質(zhì)量是密切相關(guān)的,不能反映其客觀地理分布格局和可能的潛在分布區(qū)域。

圖1 全國(guó)檵木分布點(diǎn)Fig.1 The distribution point of Lorpetalum chinense in China

2.2 檵木空間分布格局

檵木空間分布的經(jīng)度范圍在102—122° E,分布的平均經(jīng)度是113.07° E,標(biāo)準(zhǔn)差是4.279°；檵木分布的緯度范圍23—34° N,分布的平均緯度是27.03° N,標(biāo)準(zhǔn)差是2.235°(圖2)；這表明檵木在經(jīng)度分布范圍較寬,而緯度分布范圍較窄。數(shù)據(jù)離散程度的變異系數(shù)表明經(jīng)度、緯度的集中程度,變異系數(shù)值越小代表穩(wěn)定性越高(或者不改變性)。經(jīng)度因子的變異系數(shù)值為0.038,緯度因子的變異系數(shù)值為0.082,即經(jīng)度因子的穩(wěn)定性較大,較緯度因子而言,我國(guó)檵木分布格局主要受經(jīng)度因子的影響。

圖2 檵木經(jīng)度、緯度頻率分布Fig.2 Longitude and latitude frequency distribution of Lorpetalum chinense

2.3 降水量分析

根據(jù)檵木物種分布數(shù)據(jù)和氣候因子的描述統(tǒng)計(jì)表(表2),8個(gè)降水量因子變異系數(shù)絕對(duì)值大小順序：濕度變化方差(bio15)<年降水量(bio12)<暖季降水量(bio18)<雨季降水量(bio16)<最濕月降水量(bio13)<干季降水量(bio17)<最干月降水量(bio14)<冷季降水量(bio19)。因此,濕度變化方差、年降水量、暖季降水量為限制檵木分布的主要降水因子。

2.4 溫度因子分析

表1可知,11個(gè)溫度因子中變異系數(shù)的絕對(duì)值大小順序：最熱月份極高溫(bio5)<最熱季均溫(bio10)<溫差日較差值(bio2)<雨季度均溫(bio8)<年溫度較差(bio7)<年均氣溫(bio1)<晝夜溫差與年溫差百分比(bio3)<溫度變化方差(bio4)<最冷季均溫(bio11)<最干季均溫(bio9)<最冷月份極低溫(bio6)。因此,最熱月份極高溫、最熱季均溫、溫差日較差值為限制檵木分布的主要溫度因子。

表2 氣候因子的描述統(tǒng)計(jì)結(jié)果

2.5 基于刀切法的環(huán)境變量重要性分析

圖3 刀切法檢測(cè)主要環(huán)境變量對(duì)檵木在中國(guó)分布影響的重要程度 Fig.3 Importance of major environmental variables on distribution of Loropetalum chinense in China by Jackknife test bio1: 年平均氣溫；bio2: 晝夜溫差月均值；bio3: 晝夜溫差與年溫差百分比；bio4: 溫度變化方差；bio5: 極熱月份最高溫；bio6: 極冷月份最低溫；bio7: 年溫度變化范圍；bio8: 雨季均溫；bio9: 極干季度均溫；bio10: 極熱季度均溫；bio11: 極冷季度均溫；bio12: 年均降雨量；bio13: 極濕月降雨量；bio14: 極干月降雨量；bio15: 濕度變化方差；bio16: 極濕季降雨量；bio17: 極干季降雨量；bio18: 極熱季平均降雨量；bio19: 極冷季平均降雨量；Alt:海拔

受到極端氣候條件影響下,預(yù)測(cè)物種潛在分布格局的限制更為顯著[11,17],因此,按照影響數(shù)值大的氣候因子中,擇優(yōu)選擇反映極值氣候因子的原則。由圖3可知,用刀切法(Jackknife)檢測(cè)在基于氣候溫度、水分等要素所選擇的19個(gè)環(huán)境變量對(duì)檵木適生性分布的貢獻(xiàn)大小。bio14(最干月份降雨量)對(duì)應(yīng)的深藍(lán)條最長(zhǎng),3個(gè)時(shí)期預(yù)測(cè)百分比均達(dá)到了75%以上,表明它對(duì)預(yù)測(cè)檵木的分布有重要影響。按照氣候因子得分值的大小,影響檵木潛在分布的環(huán)境因子還包括bio17(最干季度降雨量),其次為bio19(最冷季度降雨量),bio15(濕度變化方差),bio2(晝夜溫差月均值)。溫度的高低和水分的有效性都會(huì)影響檵木的生理活動(dòng)和生化過(guò)程。

模型生成的各氣候因子響應(yīng)曲線顯示預(yù)測(cè)分布概率隨環(huán)境變量的變化趨勢(shì),響應(yīng)曲線反映各環(huán)境變量對(duì)模型預(yù)測(cè)的影響。一般認(rèn)為,當(dāng)分布概率值>0.5時(shí),其對(duì)應(yīng)的生態(tài)因子的值適合植物生長(zhǎng)[4]。圖4可見(jiàn),bio14(最干月降雨量)、bio17(極干季降雨量)、bio19(極冷季平均降雨量)和bio15(濕度變化方差)的響應(yīng)曲線均呈單峰曲線,說(shuō)明檵木對(duì)這些環(huán)境變量有明顯的適應(yīng)性。bio14(最干月降雨量)小于15 mm時(shí),檵木分布概率幾乎為零；降雨量大于15 mm分布概率急劇上升,降雨量約45 mm時(shí),分布概率達(dá)到峰值,此時(shí)降雨量范圍是檵木最適的生存條件。隨后分布概率又急劇降低,當(dāng)降雨量達(dá)到170 mm左右時(shí),分布概率幾乎減小為0。當(dāng)分布概率>0.5時(shí)稱(chēng)為生存概率,即生存概率在0.5—1之間為適生范圍,檵木最干月降雨量的適應(yīng)范圍約為23—93 mm。除bio14外,對(duì)當(dāng)代檵木生態(tài)適宜區(qū)預(yù)測(cè)貢獻(xiàn)最大bio17(極干季降雨量),降雨量達(dá)到160 mm時(shí),分布概率達(dá)到峰值,此時(shí)降雨量達(dá)到檵木最適生存條件；若以生存概率>0.5適應(yīng)范圍,降雨量的適應(yīng)范圍約98—300 mm。bio19(極冷季平均降雨量)在200 mm時(shí)達(dá)到峰值,為檵木最適的生存條件；若以生存概率>0.5為適應(yīng)范圍,降雨量的適應(yīng)范圍約110—350 mm。bio15(濕度變化方差)濕度范圍在52%時(shí)達(dá)到峰值；生存概率>0.5為適應(yīng)范圍,濕度變化方差范圍在42%—65%。

圖4 主要?dú)夂蛞蜃拥捻憫?yīng)曲線Fig.4 Response curves of major climate factors

2.6 模型評(píng)估

從Maxent軟件輸出的遺漏率、發(fā)生概率與累積閾值關(guān)系函數(shù)圖(圖5)中可以看出測(cè)試集遺漏率與預(yù)測(cè)遺漏率基本吻合,說(shuō)明模型與訓(xùn)練數(shù)據(jù)的擬合越好,測(cè)試數(shù)據(jù)和培訓(xùn)數(shù)據(jù)是獨(dú)立的。

基于最大熵模型,對(duì)中國(guó)地區(qū)金縷梅科檵木在末次冰期(LGM)、末次間冰期(LIG)和當(dāng)代這3個(gè)不同時(shí)期的潛在分布區(qū)進(jìn)行了模擬,結(jié)果如圖6所示。在所有模擬過(guò)程中,訓(xùn)練集和驗(yàn)證集AUC數(shù)值分別為0.947和0.954,模擬結(jié)果數(shù)值大于隨機(jī)模型數(shù)值(0.5),說(shuō)明投影分布與所收集數(shù)據(jù)吻合。

圖5 測(cè)試集遺漏率Fig.5 Omission on test sample

圖6 檵木預(yù)測(cè)結(jié)果AUC圖Fig.6 Lorpetalum chinense prediction the result of AUC chart AUC：曲線下的面積 Area Under Curve

2.7 檵木分布格局模擬分析

圖7 LIG時(shí)期檵木在中國(guó)的潛在地理分布圖Fig.7 LIG period potential geographical distribution of Lorpetalum chinense in China LIG: 末次間冰期last interglacial

檵木適生等級(jí)劃分應(yīng)用ArcGIS軟件將MaxEnt軟件輸出的ASCⅡ格式文件轉(zhuǎn)為Raster格式文件,利用ArcGIS的統(tǒng)計(jì)分析功能計(jì)算在各等級(jí)的分布面積[8, 28]。MaxEnt模型默認(rèn)適生等級(jí)為10級(jí),根據(jù)檵木在全國(guó)的實(shí)際分布情況,利用ArcGIS的“重分類(lèi)工具”將檵木的生境適宜性按照“自然間斷點(diǎn)分級(jí)法”劃分為4個(gè)等級(jí)[29]：不適生區(qū)：P<0.23；低適生區(qū)：0.23≤P<0.54；中適生區(qū)：0.54≤P<0.85；高適生區(qū)：0.85≤P<1。

在LIG時(shí)期,檵木總面積為 43.48 萬(wàn)km2,占國(guó)土總面積的4.53%(表3),其中高、中、低度適生區(qū)總面積分別是0.56 萬(wàn)km2,3.51 萬(wàn)km2,39.41 萬(wàn)km2,占全國(guó)總面積的0.05%、0.37%和4.1%。分布區(qū)主要在廣西、廣東、福建浙江南部以及臺(tái)灣北部區(qū)域,檵木的高度適生區(qū)集中在臺(tái)灣北部、中適生區(qū)分布于福建東南部(圖7)。與當(dāng)代和LGM時(shí)期相比,LIG時(shí)期適生區(qū)向南移,有明顯的收縮以及片段化現(xiàn)象,適生區(qū)總面積占比最小。

末次冰期LGM,檵木的潛在地理分布區(qū)總面積為153.76 萬(wàn)km2,占國(guó)土總面積的16.01%。其中中、低度適生區(qū)總面積分別是103.83 萬(wàn)km2,49.93 萬(wàn)km2,占全國(guó)總面積的10.82%和5.2%(表3)。見(jiàn)圖8可知檵木中適生區(qū)的區(qū)域主要集中分布在福建、浙江南部、重慶西部、貴州東部、湖南南部以及江西有零星分布。末次冰期(LGM)沒(méi)有明顯的高適宜區(qū)分布,低適生區(qū)占比最大。相比末次間冰期,檵木在LGM時(shí)期的潛在分布區(qū)向北遷移。

圖8 LGM(末次冰期)檵木在中國(guó)的潛在地理分布圖Fig.8 LGM period potential geographical distribution of Lorpetalum chinense in ChinaLGM: 末次冰期last glacial maximum

圖9 當(dāng)前氣候條件下檵木在中國(guó)的潛在地理分布圖Fig.9 Current potential geographical distribution of Lorpetalum chinense in China

表3 當(dāng)前、LIG、LGM檵木在我國(guó)的適生面積

當(dāng)代我國(guó)檵木的潛在地理分布主要在北回歸線以北的區(qū)域,總面積為162.55 萬(wàn)km2,占國(guó)土總面積的16.93%,其中高、中、低度適生區(qū)總面積分別是8.36 萬(wàn)km2,102.4 萬(wàn)km2,51.79 萬(wàn)km2,占全國(guó)總面積的0.87%、10.67%和5.39%(表3)。高度適生區(qū)主要分布于廣西北部、江西東部、以及湖南省南部、福建東部等地區(qū)。低度適生區(qū)主要分布于湖北、貴州西部、四川南部、廣西南部、安徽南部。中度適生區(qū)主要分布于重慶、貴州、湖南、廣東、江西、福建以及湖南東部(圖7)。當(dāng)代適生區(qū)總面積相比于LIG,LGM兩個(gè)時(shí)期的潛在分布區(qū)總面積有明顯擴(kuò)增的現(xiàn)象。在氣候變暖的條件下,檵木適生區(qū)域增大,適生程度也有所提高。

3 結(jié)論與討論

已有的關(guān)于檵木生境及適宜性的研究,地域記載于1979年《中國(guó)植物志》北回歸以南未見(jiàn)蹤跡,根據(jù)全球生物多樣性信息平臺(tái)和中國(guó)數(shù)字植物標(biāo)本館(http://www.cvh.org.cn/)等數(shù)據(jù)庫(kù)獲得當(dāng)代檵木標(biāo)本的采集地信息,預(yù)測(cè)結(jié)果表明：檵木當(dāng)代主要分布區(qū)為江西、貴州、四川、重慶、安徽、湖南、廣西、廣東、福建等地區(qū),集中在亞熱帶季風(fēng)性濕潤(rùn)性氣候,因而具有連續(xù)性的分布特性,這與檵木陽(yáng)性喜光、喜溫暖、濕潤(rùn)的微酸性土壤等環(huán)境需求一致[9]。檵木經(jīng)度分布點(diǎn)范圍在102°—123°E、緯度分布點(diǎn)范圍在23°—34°N,與中國(guó)植物志中對(duì)該物種的分布范圍描述基本一致[30]。

通過(guò)AUC檢驗(yàn)的結(jié)果,可以看出Maxent模型對(duì)檵木的預(yù)測(cè)結(jié)果是比較可靠的,Maxent模型是基于當(dāng)前存在數(shù)據(jù)進(jìn)行分析的模型,模型結(jié)果可用于當(dāng)前存在的檵木種群。全國(guó)尺度下對(duì)檵木適宜性分布的主導(dǎo)環(huán)境因子有最干月份平均降雨量、極干季降雨量、極冷季平均降雨量、濕度變化方差,其中降水因子起主導(dǎo)作用。本研究表明,分布區(qū)降水變化特征可作為研究檵木生生物學(xué)特征的重要?dú)夂蛞蛩?。?項(xiàng)環(huán)境因子的適宜范圍可知,檵木適宜生長(zhǎng)在濕潤(rùn)的環(huán)境中[9]。洪文君等的研究發(fā)現(xiàn),檵木具有較強(qiáng)的適應(yīng)性以及耐陰性特征、自身對(duì)水分條件有較高的依賴(lài)性是一致的[31]?；贛axent模型得出的檵木潛在分布區(qū)域均具有相似的氣候特點(diǎn),作為生態(tài)恢復(fù)力較強(qiáng)的物種[12],可為分布區(qū)進(jìn)行植被恢復(fù)時(shí)提供樹(shù)種選擇。

末次間冰期與末次冰期檵木適宜性劃分發(fā)現(xiàn)：LGM時(shí)期適宜區(qū)的分布明顯發(fā)生了北移,有分布擴(kuò)增的現(xiàn)象,與方炎明等研究的變遷分析表明：末次間冰期開(kāi)始轉(zhuǎn)暖,氣候的波動(dòng)影響植被帶的遷移,降雨量增多,植被帶呈現(xiàn)北移狀態(tài)的結(jié)果是一致的[32]。Maxent模型對(duì)檵木在間冰期至當(dāng)代的適生區(qū)進(jìn)行預(yù)測(cè)時(shí)發(fā)現(xiàn),其高適生區(qū)范圍向南進(jìn)行擴(kuò)張。與 Gong等研究的中國(guó)亞熱帶常綠闊葉林中的檵木的種群模擬,表明氣候變化導(dǎo)致亞熱帶植被物種的地理分布發(fā)生南移的結(jié)果是一致的[33],在喀斯特地貌石漠化恢復(fù)過(guò)程中,檵木具備適應(yīng)資源貧瘠的能力和較高的生產(chǎn)力,屬于優(yōu)勢(shì)種和先鋒種[34]。有助于改善初期的生態(tài)環(huán)境,為石漠化地區(qū)進(jìn)行植被恢復(fù)提供理論依據(jù)。因此需要相關(guān)部門(mén)緊密配合,加強(qiáng)檵木的保護(hù)和推廣工作。