賀曉慧 高 健 朱 麗 郝瑞敏 黃 磊 朱 晉程 莉 周 潔

（1.內(nèi)蒙古科技大學(xué)包頭師范學(xué)院生態(tài)環(huán)境學(xué)院，包頭 014030； 2.中國科學(xué)院西北生態(tài)環(huán)境資源研究院 內(nèi)陸河流域生態(tài)水文重點(diǎn)實(shí)驗(yàn)室， 蘭州 730000）

氣候變化是目前全球所面臨的主要環(huán)境問題，已成為21 世紀(jì)威脅生物多樣性的主要因素之一［1-2］。在全球氣候變暖與人類活動(dòng)共同作用下，全球植物蓋度、豐度、優(yōu)勢種群、物種分布范圍、物候，甚至植被演替和生態(tài)系統(tǒng)穩(wěn)定性都發(fā)生了明顯變化［2-5］。眾多研究表明，氣候變化可能會(huì)導(dǎo)致物種未來適宜生境的縮減或擴(kuò)張，適宜生境向高海拔或高緯度地區(qū)遷移，對(duì)未來物種多樣性保護(hù)和利用、生態(tài)安全的維護(hù)提出了嚴(yán)峻挑戰(zhàn)［6-8］。因此，研究物種地理分布與氣候的關(guān)系是傳統(tǒng)植物生態(tài)學(xué)和生物地理學(xué)的研究熱點(diǎn)，也是全球生態(tài)學(xué)和全球變化生物學(xué)的重點(diǎn)研究領(lǐng)域之一。

物種分布模型（SDMs）是以生態(tài)位理論為基礎(chǔ)，利用物種分布數(shù)據(jù)和環(huán)境變量數(shù)據(jù)，模擬和預(yù)測物種潛在分布的有力工具，通常分為回歸算法（廣義線性模型（GLM）、廣義相加模型（GAM）、多適應(yīng)回歸樣條模型（MARS）），包絡(luò)算法（表面分布區(qū)分室模型（SRE）、棲息地模型（HABITAT）），分類算法（柔性判別分析（FDA）、分類回歸樹分析（CTA）），常用機(jī)器學(xué)習(xí)算法（人工神經(jīng)網(wǎng)絡(luò)（ANN）、支持向量機(jī)（SVM）、最大似然法（MAXLIKE）、隨機(jī)森林（RF）），模糊數(shù)學(xué)（模糊物元模型（FME）、模糊神經(jīng)網(wǎng)絡(luò)（FNN））和最大熵模型（MaxEnt）［5，7，9-10］。其中，回歸算法的精度解釋能力較好，但對(duì)樣本數(shù)量依賴較大；包絡(luò)算法簡單易用，但氣候變量作用等同，精度較低［5，10］；常用機(jī)器學(xué)習(xí)算法盡管精確度較高，使用較容易，但缺乏生態(tài)學(xué)意義［10］；MaxEnt 模型由于界面友好，方法成熟，在樣本量較小的情況下仍能取得高的模擬精度，被學(xué)者們廣泛使用［9-11］，但也造成了模型的濫用。因此，僅靠單一物種分布模型通常難以準(zhǔn)確地模擬物種適宜生境，也容易出現(xiàn)過度擬合［10，12］。隨著研究的不斷深入及對(duì)模擬精度要求的提高，綜合適宜生境模型（CHS）應(yīng)運(yùn)而生［10，12-13］。該模型不僅利用氣候因素創(chuàng)建了氣候適宜性模型，還可以利用土壤等因素創(chuàng)建限制性模型［5，9，13］。研究表明，綜合適宜生境模型較單一模型具有更高的建模精度，降低了模型不確定性［13］，近年來在生物地理學(xué)、保護(hù)生物學(xué)等領(lǐng)域的研究中得到越來越多的應(yīng)用［7，9-10］。

珍稀瀕危物種蒙古扁桃（Amygdalus mongolica）是蒙古高原第三紀(jì)孑遺的代表物種，因其具有藥用、觀賞、水土保持等多種經(jīng)濟(jì)、社會(huì)和生態(tài)價(jià)值而被廣泛關(guān)注［14］。蒙古扁桃為強(qiáng)旱生小灌木，主要分布于我國西北地區(qū)荒漠及荒漠草原區(qū)的低山丘陵地、石質(zhì)坡地及干河床等生境條件較為嚴(yán)酷的地區(qū)，是維持當(dāng)?shù)厣鷳B(tài)環(huán)境穩(wěn)定的重要物種［15］。隨著氣候環(huán)境變化、人類活動(dòng)及野生動(dòng)物干擾，蒙古扁桃種群遭到嚴(yán)重破壞，使其處于瀕危狀態(tài)。近年來，對(duì)蒙古扁桃的研究主要集中在遺傳多樣性［15］、群落特征［16］、藥用價(jià)值［17］、地理分布［18-20］，等。而氣候變化下蒙古扁桃地理分布變化研究主要采用單一模型（MaxEnt）進(jìn)行模擬，同時(shí)受物種種群樣點(diǎn)數(shù)較少、環(huán)境指標(biāo)選取單一等制約［16-18］，模擬精度有待明確。因此，本研究基于蒙古扁桃的種群分布點(diǎn)位數(shù)據(jù)，結(jié)合環(huán)境氣候因子和環(huán)境土壤因子，通過構(gòu)建物種綜合適宜生境模型（CHS），模擬當(dāng)前和預(yù)測未來（2050S 和2090S）共3 種 不 同 氣 候 情 景 下（SSP1-2.6、SSP3-7.0 和SSP5-8.5）蒙古扁桃的潛在地理分布，探究其在氣候變化下適宜生境的變化趨勢及遷移方向，以期為蒙古扁桃的種質(zhì)資源管理和保護(hù)提供理論依據(jù)。

1 數(shù)據(jù)來源與研究方法

1.1 蒙古扁桃物種分布數(shù)據(jù)來源

蒙古扁桃當(dāng)前在中國的地理分布數(shù)據(jù)主要通過以下幾種途徑獲?。海?）中國數(shù)字植物標(biāo)本館（CVH，https：//www.cvh.ac.cn）；（2）全球生物多樣性信息網(wǎng)絡(luò)（GBIF，https：//www.gbif.org）；（3）公開發(fā)表的文獻(xiàn)［15-16，18-20］；（4）2020 年和2021 年的野外調(diào)查數(shù)據(jù)等。對(duì)只有位置信息沒有坐標(biāo)信息的種群點(diǎn)，利用Google Earth 定位獲取。并利用ENMTools 軟件，結(jié)合2.5′（～5 km×5 km）的環(huán)境數(shù)據(jù)篩選分布數(shù)據(jù)，最終確定121個(gè)蒙古扁桃有效分布點(diǎn)進(jìn)行建模（見圖1）。

1.2 環(huán)境數(shù)據(jù)獲取及預(yù)處理

選用氣候、土壤兩類環(huán)境變量。氣候數(shù)據(jù)源于世界氣候數(shù)據(jù)庫（WORLD-CLIM version 2.1，http：//www.worldclim.org），包括當(dāng)前、2050S（2041—2060）和2090S（2081—2100）3 個(gè)時(shí)期的氣候和降水等19 個(gè)環(huán)境氣候數(shù)據(jù)。當(dāng)前氣候數(shù)據(jù)為1970～2000年的平均值［21］，未來數(shù)據(jù)采用CMIP 6的BCCCSM2-MR 全球氣候模式下的環(huán)境氣候數(shù)據(jù)，空間分辨率為2.5′（～5 km×5 km）。未來氣候變化情景包括：SSP1-2.6（低輻射強(qiáng)迫情景，假設(shè)2100 年輻射強(qiáng)迫穩(wěn)定在2.6 W·m-2）、SSP3-7.0（中等至高等輻射強(qiáng)迫情景，假設(shè)2100 年輻射強(qiáng)迫穩(wěn)定在7.0 W·m-2）、SSP5-8.5（高輻射強(qiáng)迫情景，假設(shè)2100年輻射強(qiáng)迫穩(wěn)定在8.5 W·m-2）。土壤數(shù)據(jù)采用世界和諧土壤數(shù)據(jù)庫（HWSD）數(shù)據(jù)，共17 個(gè)土壤變量（見表1）。

表1 蒙古扁桃地理分布模擬所需環(huán)境變量Table 1 Environment variables used for predicting spatial distribution of A.mongolica

為避免環(huán)境變量多重共線性導(dǎo)致的模型過度擬合［7］，本研究對(duì)19 個(gè)環(huán)境氣候變量進(jìn)行相關(guān)分析，選擇相關(guān)系數(shù)小于0.7 的變量。最終確定參與建模的氣候變量為：Bio1（年平均氣溫）、Bio2（平均日較差）、Bio4（溫度季節(jié)性變化）、Bio12（年降水量）、Bio15（降水量季節(jié)性變化）、Bio19（最冷季降水量）（見表1）。

1.3 模型構(gòu)建

采用綜合適宜生境模型（CHS）模擬不同氣候條件下蒙古扁桃在中國的潛在分布，該模型主要包括2 個(gè)獨(dú)立模型，即基于EM 策略模擬氣候變化情景下遷移趨勢的氣候適宜性模型和基于MaxEnt模擬擴(kuò)散限制的分布限制模型（Li）［5，9］。

氣候適宜性模型基于R 語言的Biomod2 構(gòu)建［5］。首先，生成3 套偽采樣點(diǎn)數(shù)據(jù)集，每個(gè)數(shù)據(jù)集包括500個(gè)偽采樣點(diǎn)數(shù)據(jù)，以減少由隨機(jī)生成點(diǎn)位引起的模擬不確定性。其次，選取10 種模型分別建模，包括：廣義線性模型、廣義相加模型、多元自適應(yīng)回歸樣條模型、柔性判別分析、廣義增強(qiáng)模型（GBM）、分類樹分析、人工神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)、隨機(jī)森林、最大熵模型。隨機(jī)選擇80%的樣本數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)，剩余20%的樣本數(shù)據(jù)用于模型驗(yàn)證，算法重復(fù)5 次，對(duì)樣本數(shù)據(jù)集采用真實(shí)技能統(tǒng)計(jì)（TSS）和受試者工作特征曲線（ROC）進(jìn)行精度評(píng)估。TSS 和AUC 值越大表示模型精度越高［7，9］，本研究選擇TSS＞0.70，ROC＞0.80 的模型構(gòu)建EM，模型權(quán)重由TSS值決定：

式中：EM，i為EM 評(píng)價(jià)單元i的值，Wj為模型j的 權(quán) 重，Xi，j為 模 型j結(jié) 果 中 評(píng) 價(jià) 單 元i的 值。EM為0～1 的連續(xù)值，使用TSSmax為閾值對(duì)其進(jìn)行二值化，劃分為適宜生境與不適宜生境［7，13］。

分布限制模型利用MaxEnt 模擬蒙古扁桃土壤適宜性。輸入17個(gè)環(huán)境土壤變量及蒙古扁桃分布數(shù)據(jù)，在MaxEnt 3.3.3 軟件中進(jìn)行5 倍交叉驗(yàn)證來訓(xùn)練模型，并用ROC 評(píng)估模型精度。把能夠確保蒙古扁桃分布點(diǎn)的90%正確納入適宜生境，作為閾值，將蒙古扁桃土壤適宜性劃分為土壤適宜生境和土壤不適宜生境［5，7］。

基于氣候適宜性模型與分布限制模型，建立CHS 模型模擬蒙古扁桃潛在分布，并將每個(gè)評(píng)價(jià)單元的CHS值定義為：

式 中：CHS，i表 示 評(píng) 價(jià) 單 元i的 綜 合 生 境 適 宜性，EM，i是評(píng)價(jià)單元i在不同氣候變化情景下的氣候適宜性模型結(jié)果，Li是評(píng)價(jià)單元i的分布限制模型結(jié)果。

2 結(jié)果與分析

2.1 模型精度

結(jié)果表明，不同模型算法的精度存在差異，SRE 模型的TSS 值略小于0.70，其余模型的TSS 值均大于0.70，ROC 值均大于0.80（見圖2）；且氣候適宜性模型（EM）的TSS 值大于0.75，分布限制模型的ROC值為0.939，證明模型結(jié)果較好。

圖2 EM 中10 種模型算法的真實(shí)技能統(tǒng)計(jì)量（TSS）和受試者工作特征曲線（ROC）Fig.2 The true skill statistic（TSS） and receiver operating characteristic curve（ROC） for the ensemble model（EM） for A. mongolica

2.2 蒙古扁桃當(dāng)前和未來氣候情景下的適宜生境分布

綜合適宜生境模型精細(xì)地模擬了蒙古扁桃在全國的適宜生境（見表2），當(dāng)前氣候情景下，蒙古扁桃適宜生境面積為64.00 萬km2，占全國總面積的6.67%（見表2）。由圖3 可見，蒙古扁桃在內(nèi)蒙古自治區(qū)主要分布于阿拉善盟（阿拉善左旗和阿拉善右旗）、烏海市、鄂爾多斯市、巴彥淖爾市、包頭市、呼和浩特市、集寧市、錫林郭勒盟西南部及赤峰市西部，適宜生境面積大且成片分布；在甘肅省主要分布于張掖市、金昌市、武威市、蘭州市、白銀市及臨夏市，酒泉市有零星分布；在寧夏自治區(qū)主要分布于銀川市、吳忠市、固原縣及西峰市等偏北部地區(qū)；在青海省主要分布于西寧市、平安縣、共和縣、烏蘭縣東部等東北部地區(qū)；在新疆自治區(qū)主要分布于阿瓦提縣東部、庫爾勒市北部、哈密市東部等，呈條帶狀分布；在陜西省主要分布于榆林市、延安市；在山西省主要分布于大同市，朔州市也有少量破碎狀分布；在河北省張家口市有少量片狀分布。

圖3 當(dāng)前氣候情景下蒙古扁桃適宜生境分布Fig.3 Distribution of suitable habitat of A. mongolica under current climate scenarios

表2 不同氣候情景下蒙古扁桃在中國的適宜生境面積Table 2 Suitable habitat area of A. mongolica under different climate scenarios

未來蒙古扁桃的適宜生境面積整體呈減少趨勢。 在2050S 和2090S，SSP1-2.6、SSP3-7.0 和SSP5-8.5的氣候情景下，蒙古扁桃適宜生境面積分別 為62.62 萬km2、59.63 萬km2、58.44 萬km2和60.85 萬km2、57.92 萬km2、50.80 萬km2。蒙古扁桃適宜生境面積隨著輻射強(qiáng)迫的增大呈現(xiàn)出縮減的趨勢（見表2）。此外，從時(shí)間尺度看，2090S 較2050S蒙古扁桃適宜生境面積也呈現(xiàn)出縮減趨勢。從圖4 可見，隨著時(shí)間推移和輻射強(qiáng)迫增大，內(nèi)蒙古中西部，甘肅中部、河西走廊，寧夏中北部一直為蒙古扁桃的適生分布區(qū)，而內(nèi)蒙古東部赤峰市的蒙古扁桃適宜生境面積逐漸消失，內(nèi)蒙古與陜西和山西接壤處，陜西、山西、河北等地的蒙古扁桃適宜生境面積也在逐步喪失。

圖4 未來不同氣候情景下蒙古扁桃適宜生境分布Fig.4 Distribution of suitable habitat of A. mongolica under different climate scenarios in the future

2.3 蒙古扁桃未來氣候情景下的適宜生境分布變化

在綜合適宜生境模型的優(yōu)勢下，本研究發(fā)現(xiàn)，在全球氣候變化下，未來蒙古扁桃適宜生境退化區(qū)域和適宜性增加區(qū)域都分布在穩(wěn)定適生區(qū)的邊緣，穩(wěn)定適生區(qū)主要分布在內(nèi)蒙古阿拉善盟（阿拉善左旗和阿拉善右旗）、巴彥淖爾市、烏海市、包頭市中北部、鄂爾多斯市，寧夏石嘴山市、銀川市和吳忠市，甘肅蘭州、白銀、武威、金昌和張掖市（見圖5）。適宜性退化區(qū)域主要位于內(nèi)蒙古赤峰市、錫林浩特市、集寧市南部、呼和浩特市、鄂爾多斯東南部，山西省大同市、朔州市，陜西省的榆林市、延安市，河北張家口市，甘肅西峰市等（見圖5）。適宜性增加區(qū)域主要位于內(nèi)蒙古阿拉善盟，青海烏蘭縣，新疆哈密市、庫爾勒市等，且有向西方向的擴(kuò)張。蒙古扁桃未來不同氣候情景下適宜穩(wěn)定性比例大于75%，2050S的平均適宜穩(wěn)定性比例高于2090S（見表3）。

圖5 未來不同氣候情景下蒙古扁桃適宜生境變化Fig.5 The change of suitable habitat of A. mongolica under different climate scenarios in the future

表3 未來不同氣候情景下蒙古扁桃適宜生境穩(wěn)定比例Table 3 Stable habitat percentage of A. mongolica under different climate scenarios in the future

蒙古扁桃從當(dāng)前到未來2050S、2090S 氣候情景下，適宜生境質(zhì)心整體上有逐步向西和向高緯度地區(qū)遷移的趨勢。但不同時(shí)期不同氣候情景下質(zhì)心的遷移方向和距離具有差異性。當(dāng)前氣候情景下，蒙古扁桃適宜生境質(zhì)心位于內(nèi)蒙古自治區(qū)阿拉善左旗（39°56′37″N，105°2′2″E）。在2050S，SSP1-2.6 氣候情景下，適宜生境質(zhì)心向東南方向遷移，遷移距離較短。但隨著輻射強(qiáng)迫的增強(qiáng)，適宜生境質(zhì)心逐步向西遷移，且遷移距離逐步加大（見圖6）。在2090S，SSP1-2.6 氣候情景下，適宜生境質(zhì)心較當(dāng)前向西遷移，遷移距離較短且向高緯度地區(qū)遷移，隨后隨著輻射強(qiáng)迫增強(qiáng)，適宜生境質(zhì)心向西南方向遷移。總體上看，隨著時(shí)間和氣候情景的變化，蒙古扁桃適宜生境質(zhì)心從當(dāng)前的阿拉善左旗遷移到了阿拉善右旗（見圖6）。

圖6 不同氣候情景下蒙古扁桃適宜生境質(zhì)心的變化Fig.6 The centroid changes of suitable habitat of A. mongolica under different climate scenarios

3 討論

3.1 模型、指標(biāo)及種群點(diǎn)位對(duì)模擬精度的影響

蒙古扁桃生態(tài)位模型模擬結(jié)果的精確度受多方面因素影響。首先是模型的選擇，本研究采用物種綜合適宜生境模型［5，7，9］模擬和預(yù)測蒙古扁桃在不同氣候情景的適宜生境，發(fā)現(xiàn)不同模型算法的精度存在差異，SRE 模型的TSS 值較小，RF 模型的TSS 值較大（見圖2）；本研究氣候適宜性模型的TSS 值大于0.75，分布限制模型（MaxEnt）的ROC值為0.939，表明本研究選用綜合適宜生境模型取得了較高的模擬精度，模擬結(jié)果可信。同時(shí)該模型也更加適合于特殊生境物種地理分布及其對(duì)氣候變化的響應(yīng)研究［10，13］。而眾多研究也表明，構(gòu)建組合模型可以避免單一模型產(chǎn)生擬合不足或擬合過度的現(xiàn)象［12，22］，使模擬精度更高、通用性更強(qiáng)、穩(wěn)健性更好［7，9，12，23］。其次，關(guān)于環(huán)境變量的選擇，并不是參與建模的環(huán)境變量種類和數(shù)量越多越能提高預(yù)測精度［8，24］，例如：Wordclim 數(shù)據(jù)具有較強(qiáng)的共線性，使用19 個(gè)氣候變量會(huì)造成建模時(shí)信息冗余，影響擬合結(jié)果的準(zhǔn)確性［21］。一些看似對(duì)物種分布具有很強(qiáng)影響，且在生態(tài)位模擬中也具有很大貢獻(xiàn)率的環(huán)境變量（如海拔）并不適用于生態(tài)位建模［25-26］。而選擇合適的環(huán)境因子至關(guān)重要，例如：荒漠植物的生長發(fā)育對(duì)土壤環(huán)境依賴性較強(qiáng)，把對(duì)荒漠植物分布格局和擴(kuò)散具有重要影響的土壤因子納入生態(tài)指標(biāo)就顯得尤為重要［27］。這是本研究未選擇地形（海拔、坡度和坡向）因子，而增加土壤因子作為限制性變量的原因，也可能是與前人研究結(jié)果中適宜生境面積產(chǎn)生差異的主要原因［18-20］。馬松梅等［18-19］發(fā)現(xiàn)，人類對(duì)土地的開發(fā)、盲目的采礦、城市化加快等嚴(yán)重干擾和破壞了蒙古扁桃的生境，使得蒙古扁桃的生境分布區(qū)面積不斷縮小和破碎化，甚至瀕臨滅絕。因此，在未來需要不斷挖掘人類活動(dòng)與干擾等影響指標(biāo)，定量化研究其對(duì)物種生態(tài)位的影響，使物種生態(tài)位研究更加準(zhǔn)確［10，28］。再次，物種采樣點(diǎn)的位置誤差、采樣點(diǎn)的數(shù)量、均勻度及樣本的選擇性偏差均會(huì)影響模型的模擬精度［5，29］。采樣點(diǎn)越接近物種實(shí)際分布，模型模擬精度可能越高。本研究利用ENMTools 軟件，對(duì)采樣點(diǎn)數(shù)據(jù)進(jìn)行了篩選，確保每5 km×5 km 的范圍內(nèi)有至少1 個(gè)采樣點(diǎn)，使采樣點(diǎn)相對(duì)均勻；為了保證建模的精度，采用存在點(diǎn)位數(shù)據(jù)和偽不存在點(diǎn)位數(shù)據(jù)同時(shí)進(jìn)行建模，這樣建模優(yōu)于只采用存在點(diǎn)位數(shù)據(jù)［10，13］。此外，本研究模型精度效果采用了TSS 和AUC 同時(shí)評(píng)價(jià)，也有效規(guī)避了單純使用AUC 引起的不準(zhǔn)確性［30］。綜上所述，本研究經(jīng)過環(huán)境變量降維、物種分布點(diǎn)篩選等，利用綜合適宜生境模型取得了較高的模擬精度，模擬結(jié)果相對(duì)比先前的研究應(yīng)該準(zhǔn)確性更高，更能說明蒙古扁桃的適宜生境及未來氣候變化趨勢。

3.2 蒙古扁桃潛在適生區(qū)的變化格局

蒙古扁桃是亞洲戈壁荒漠區(qū)孑遺的旱生落葉灌木，具有耐高溫、耐干旱等特性，主要生長于降水少、氣溫高的西北礫質(zhì)荒漠、低山丘陵、石質(zhì)坡地及干河床等地，分布范圍相對(duì)較為狹窄［14-15］。研究表明，當(dāng)前氣候情景下蒙古扁桃的適宜生境面積僅占全國總面積的6.67%（見表2），且主要分布于蒙甘寧地區(qū)（見圖3），這基本與馬松梅等［18］、段義忠等［20］的結(jié)果相似。生物地理分布受到諸多自然地理要素的影響與制約，其中氣候要素是一個(gè)重要的影響因素［31］。隨著全球氣候變暖，未來蒙古扁桃的潛在地理分布范圍隨著時(shí)間推移及溫室氣體排放濃度的增加呈現(xiàn)出逐漸減小的趨勢（見圖4）。其中，2090S的SSP5-8.5情景下適宜生境面積較當(dāng)前減少約20.63%（見表2）。全球氣候變化對(duì)高海拔地區(qū)物種的影響相對(duì)顯著［32］，蒙古扁桃主要分布于內(nèi)蒙古高原、黃土高原等地區(qū)，該區(qū)域范圍的地形異質(zhì)性低，對(duì)未來氣候變暖的緩沖作用較小，所以增溫是必然趨勢［33］。由于蒙古扁桃所能適應(yīng)溫度的生態(tài)幅較小［19］，在未來全球氣溫上升的大背景下，溫度變化可能成為制約蒙古扁桃生存的重要因素，從而導(dǎo)致物種分布區(qū)面積的下降。具體來看，未來蒙古扁桃分布區(qū)范圍收縮主要集中于蒙陜晉接壤處，擴(kuò)張出的新分布區(qū)主要分布于相對(duì)較高緯度（見圖5）。氣溫變暖迫使現(xiàn)今分布于蒙陜晉接壤處的蒙古扁桃種群往北拓展生存空間，陰山山脈由于海拔、緯度均比蒙陜晉接壤處高，在氣候變暖的背景下可能成為蒙古扁桃種群尋求合適生境道路上較為理想的“避難”場所。此外，陰山山脈長期以來便是蒙古扁桃分布的核心區(qū)域［14］。較多北方干旱地區(qū)植物在未來氣候變暖的背景下也出現(xiàn)了分布區(qū)縮小的現(xiàn)象，例如：我國荒漠草原防風(fēng)固沙的先驅(qū)植物——黑沙蒿（Artemisia ordosica），在未來情景下的潛在高適生區(qū)分布面積較當(dāng)前減少了5.2%～26.8%［34］；我國一級(jí)瀕危保護(hù)孑遺植物四合木（Tetraena mongolica）的未來分布區(qū)范圍將出現(xiàn)縮減趨勢［35］；第三紀(jì)古地中海孑遺植物綿刺（Potaninia mongolica）的未來高適生區(qū)出現(xiàn)收縮，尤其是在2070S 時(shí)期RCP8.5 情景下收縮更為嚴(yán)重［36］；亞洲中部荒漠特有種半日花（Helianthemum songaricum）未來的潛在適生區(qū)面積較當(dāng)代有所減少［29］。由此可見，在未來氣候變化下，干旱區(qū)植物可能存在適生區(qū)面積喪失的風(fēng)險(xiǎn)，從而可能影響物種（尤其影響瀕危保護(hù)物種）［37］的生存。

隨著未來全球溫度的不斷升高，物種為尋求與升溫前類似的生境，將出現(xiàn)向高緯度或高海拔地區(qū)遷移的趨勢［34，38-40］。例如：白刺屬物種未來均有向西北方向青海等高海拔地區(qū)遷移的趨勢［39］；四合木受溫度影響顯著，未來（2050 年）適宜生境面積將會(huì)減少，且適生區(qū)呈現(xiàn)破碎化分布，向東北、西北高緯度地區(qū)遷移［35］；未來氣候變化背景下孑遺植物桫欏（Alsophila spinulosa）的高適宜生境中心向西北地區(qū)和高緯度地區(qū)轉(zhuǎn)移［40］。本研究表明，從當(dāng)前到未來（2050S、2090S）不同氣候情景下，蒙古扁桃分布區(qū)面積邊緣出現(xiàn)收縮與擴(kuò)展之外，適宜生境的質(zhì)心也逐步由阿拉善左旗遷移至阿拉善右旗（往西約400 km）（見圖6）。有研究表明未來氣候條件下，中亞干旱區(qū)降水量會(huì)出現(xiàn)下降趨勢，但中亞偏東地區(qū)可能出現(xiàn)少量的增加［41］，這可能是導(dǎo)致蒙古扁桃出現(xiàn)適宜生境質(zhì)心遷移的因素之一。根據(jù)蒙古扁桃種群遺傳結(jié)構(gòu)分析顯示，現(xiàn)今蒙古扁桃實(shí)際的種群主要分為陰山山脈及周邊地區(qū)、賀蘭山及河西走廊兩個(gè)種群集合區(qū)［42］。通識(shí)意義上來講，物種一般會(huì)通過種群擴(kuò)散來維持或發(fā)展種群數(shù)量并保證遺傳多樣性，進(jìn)而實(shí)現(xiàn)種群的延續(xù)。對(duì)于蒙古扁桃賀蘭山及河西走廊種群集合區(qū)而言，阿拉善左旗以西，河西走廊以北的阿拉善右旗地區(qū)是較為理想的種群擴(kuò)展區(qū)。此外，阿拉善左旗緊鄰內(nèi)蒙古烏海市、寧夏自治區(qū)，頻繁的人類活動(dòng)（例如放牧、砍柴、土地開發(fā)）可能會(huì)對(duì)蒙古扁桃的生境產(chǎn)生破壞［18］，從而影響其交配系統(tǒng)及種群的基因流，導(dǎo)致物種分布區(qū)發(fā)生變化［43］。

4 結(jié)論

本研究表明物種綜合適宜生境模型能夠準(zhǔn)確地模擬和預(yù)測當(dāng)前和未來蒙古扁桃的適宜生境。蒙古扁桃作為我國西北干旱區(qū)重要的瀕危物種，當(dāng)前氣候情景下主要分布于蒙甘寧等地，少量分布于青海、新疆、陜西、山西、河北等地；未來隨著全球氣候變暖，蒙古扁桃的適宜生境面積隨著時(shí)間推移及碳排放濃度增加均呈現(xiàn)出減少趨勢，但適宜生境穩(wěn)定比例在68.49%～87.79%，暫無滅絕風(fēng)險(xiǎn)。蒙古扁桃適宜生境質(zhì)心整體上向西和向高緯度地區(qū)遷移，分布格局呈現(xiàn)“西北擴(kuò)”與“東南縮”的變化趨勢。相關(guān)研究結(jié)果能夠?yàn)槊晒疟馓乙吧N質(zhì)資源保護(hù)和撫育、人工合理種植及種群恢復(fù)與重建等提供科學(xué)參考與理論依據(jù)，對(duì)于維持西北干旱地區(qū)生物多樣性和維持生態(tài)環(huán)境穩(wěn)定具有重要意義。