氣候變化對(duì)荒漠區(qū)子午沙鼠(Meriones meridianus)分布格局的影響

王志成， 王 艷， 康宇坤， 張彩軍， 安 康， 蘇軍虎,3*

(1.甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院/草業(yè)生態(tài)系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室/甘肅省草業(yè)工程實(shí)驗(yàn)室/中美草地畜牧業(yè)可持續(xù)發(fā)展研究中心，甘肅 蘭州 730070； 2.甘肅農(nóng)業(yè)大學(xué)-新西蘭梅西大學(xué)草地生物多樣性研究中心， 甘肅 蘭州 730070； 3.甘肅省祁連山草原生態(tài)系統(tǒng)野外科學(xué)觀測(cè)研究站， 甘肅 天祝 733200)

氣候是影響生物多樣性的主要自然因素[1]。近年來，全球氣候變暖已經(jīng)成為21世紀(jì)人類生存和生物種群延續(xù)所面臨的重大挑戰(zhàn)，這種變化對(duì)物種的生存環(huán)境、分布范圍、種群大小、種間關(guān)系、生態(tài)系統(tǒng)結(jié)構(gòu)和功能等均已產(chǎn)生不同程度的影響[2]。根據(jù)目前溫室氣體排放趨勢(shì)，預(yù)計(jì)到2100年，全球平均溫度可能上升5.6℃[3]，這將會(huì)導(dǎo)致一些物種棲息地?cái)U(kuò)展且種群數(shù)量增多，另外一些種群的分布范圍和豐富度則會(huì)減少，如果物種對(duì)氣候變化做出的響應(yīng)遠(yuǎn)遠(yuǎn)慢于氣候變化速率，也可能會(huì)導(dǎo)致物種滅絕[4-7]。已有研究表明，物種為了響應(yīng)氣候變暖，會(huì)向更高海拔或更高緯度遷移[8]，但由于地理障礙造成的遷移限制，可能會(huì)加速一些物種的分化，甚至滅絕[9-10]。因此，準(zhǔn)確預(yù)測(cè)氣候變化對(duì)生物種群的影響，是當(dāng)前和未來物種管理和生物多樣性保護(hù)的重要內(nèi)容，同時(shí)也對(duì)維持生態(tài)系統(tǒng)平衡和生物多樣性穩(wěn)定具有重要意義。

物種分布模型(Species distribution model)主要是利用物種分布信息和環(huán)境數(shù)據(jù)，根據(jù)一定的算法運(yùn)算來構(gòu)建模型，以預(yù)測(cè)物種的實(shí)際分布(Actual distribution)和潛在分布(Potential distribution)[11]。通過物種分布模型，可以獲得關(guān)于物種與區(qū)域之間的決策信息數(shù)據(jù)，如：(1)物種的潛在發(fā)生區(qū)域大??；(2)物種當(dāng)前與未來區(qū)域之間的差異；(3)在環(huán)境變化的相互作用下適應(yīng)能力較強(qiáng)的物種；(4)在未來氣候情景中生物多樣性較高的區(qū)域等[12]。目前，有較多的物種分布模型被廣泛運(yùn)用于生境潛在分布研究，如生態(tài)位因子分析(Ecological niche factor analysis)模型、基于統(tǒng)計(jì)概率和規(guī)則集的遺傳算法(Genetic algorithm for ruleset production)模型、生物氣候(Biology climate)模型和最大熵(Maximum entropy,MaxEnt)模型等[13]。其中MaxEnt是一種基于最大熵理論的密度估計(jì)和物種地理分布預(yù)測(cè)模型，該模型通過使用已知的物種分布數(shù)據(jù)，確定物種當(dāng)前適宜生境、預(yù)測(cè)物種在未來氣候的地理分布變化[14]。該模型可以在分布點(diǎn)較少的情況下完成預(yù)測(cè)，且具有相對(duì)較高的預(yù)測(cè)精度和較強(qiáng)的預(yù)測(cè)能力[15]，目前廣泛運(yùn)用于全球變暖背景下物種的分布和氣候之間的關(guān)系、區(qū)域氣候變化對(duì)植物群落和功能的影響、生態(tài)系統(tǒng)功能群和關(guān)鍵種的監(jiān)測(cè)和預(yù)測(cè)、生態(tài)系統(tǒng)不同尺度多樣性的管理和保護(hù)、外來物種入侵區(qū)域的預(yù)測(cè)、生態(tài)系統(tǒng)的關(guān)鍵物種的潛在分布預(yù)測(cè)和保護(hù)區(qū)規(guī)劃等方面[4,16-17]。

子午沙鼠(MerionesmeridianusPallas,1773)隸屬嚙齒目(Rodentia)鼠科(Muridae)沙鼠亞科(Gerbillinae)沙鼠屬(Meriones)，是分布于歐亞大陸中部荒漠、半荒漠地區(qū)的小型嚙齒類動(dòng)物，常棲居于沙質(zhì)松土地帶和固定、半固定沙丘中，主要以環(huán)境中耐旱植物的莖、葉和種子等為食，如泡泡刺(Nitrariasphaerocarpa)、沙拐棗(Calligonummongolicm)、豬毛菜(Salsolapasserina)等，是我國荒漠地區(qū)的優(yōu)勢(shì)物種之一[18]。當(dāng)子午沙鼠以自然密度存在時(shí)，其作為“生態(tài)系統(tǒng)工程師”對(duì)維持群落的穩(wěn)定、更新和發(fā)展有顯著的作用，如影響氮含量、改善草原土壤，甚至形成“沃島效應(yīng)”[19]；然而，當(dāng)種群密度過大時(shí)，它們會(huì)破壞植被、危害農(nóng)區(qū)作物，加速荒漠草原的沙化[20]。此外，子午沙鼠還是典型的鼠疫傳播者和病菌攜帶者，對(duì)當(dāng)?shù)厝丝诮】荡嬖谳^大威脅[21]。因此，準(zhǔn)確預(yù)測(cè)未來子午沙鼠的分布范圍變化，對(duì)維持荒漠區(qū)生態(tài)系統(tǒng)平衡以及預(yù)防公共衛(wèi)生安全具有重要意義。本研究使用MaxEnt模型模擬子午沙鼠的空間分布范圍，并評(píng)估氣候變化對(duì)子午沙鼠分布格局的影響，重點(diǎn)比較了1970—2000年與2050RCP2.6(低排放情景)和2050RCP8.5(高排放情景)[22]兩種氣候情景下子午沙鼠棲息地的變化趨勢(shì)；通過確定影響當(dāng)前子午沙鼠分布的環(huán)境參數(shù)，預(yù)測(cè)子午沙鼠2050s(2041—2060年)兩種排放情景下的分布和范圍變化，為鼠疫疫情防控和風(fēng)險(xiǎn)評(píng)估工作以及子午沙鼠的防治和物種管理工作提供科學(xué)依據(jù)。

1 材料與方法

1.1 子午沙鼠地理分布數(shù)據(jù)的獲取與處理

本研究所用子午沙鼠分布位點(diǎn)主要有以下三種來源：(1)中國動(dòng)物主題數(shù)據(jù)庫(http://www.zoology.csdb.cn);(2)查閱相關(guān)書籍，包括《中國動(dòng)物志獸綱》和《中國獸類野外手冊(cè)》[23-24]；(3)文獻(xiàn)查詢獲取[25-28]。通過以上來源查詢子午沙鼠分布位置的GPS位點(diǎn)，并導(dǎo)入ArcGIS 10.9中，剔除重復(fù)出現(xiàn)和距離較近的GPS位點(diǎn)，避免空間自相關(guān)性過高引起的負(fù)面影響[28]，最終保留51個(gè)GPS位點(diǎn)(圖1)，以“.csv”格式保存。

圖1 子午沙鼠分布點(diǎn)數(shù)據(jù)Fig.1 Distribution points of midday gerbil

1.2 生物氣候變量的獲取與處理

從Worldclim(http://www.worldclim.org)獲取了代表溫度和降水參數(shù)的氣候數(shù)據(jù)，該數(shù)據(jù)提供了基準(zhǔn)(1970—2000年)全球和區(qū)域氣象站的平均氣候數(shù)據(jù)，并選取了BCC-CSM-1大氣環(huán)流模型作為未來2050s(2041—2060年)的氣候模式。此大氣環(huán)流模型是由中國國家氣候中心開發(fā)，其對(duì)溫度和降水各變量的模擬和預(yù)測(cè)都有較高的準(zhǔn)確性，已被廣泛應(yīng)用[30]。選取在未來氣候變化中溫室排放氣體最低和最高兩種排放情景，即RCP2.6和RCP8.5進(jìn)行分析，所有生物氣候變量的空間分辨率均為2.5 弧分。本文所用中國行政區(qū)劃矢量數(shù)據(jù)來自自然資源部標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站(http://bzdt.ch.mnr.gov.cn)。所有數(shù)據(jù)經(jīng)重采樣統(tǒng)一坐標(biāo)系和分辨率后轉(zhuǎn)換為“asc”格式進(jìn)行后續(xù)分析[31]。同時(shí)，通過ENMTools軟件對(duì)環(huán)境變量之間具有高度相關(guān)的生物氣候變量進(jìn)行剔除[32]，在去除11個(gè)高度相關(guān)的變量(|r|>0.80)后，對(duì)年平均氣溫(Bio1)、晝夜溫差月均值(Bio2)、最冷月最低溫(Bio6)、溫度年較差(Bio7)、最干季度均溫(Bio9)、年降水量(Bio12)、降水季節(jié)性(Bio15)、最干季度降水量(Bio17)8個(gè)生物氣候變量進(jìn)行分析(表1)。

表1 用于子午沙鼠潛在地理分布建模的生物氣候變量Table 1 Bioclimatic variables for modeling the potential geographic distribution of midday gerbil

1.3 最大熵建模和模型評(píng)估

將子午沙鼠的分布點(diǎn)數(shù)據(jù)導(dǎo)入MaxEnt v3.4.4模型中，將75%的分布點(diǎn)作為訓(xùn)練集，25%的分布點(diǎn)作為測(cè)試集，用于對(duì)模型的預(yù)測(cè)效果進(jìn)行評(píng)價(jià)[33]。獲得的測(cè)試集數(shù)據(jù)輸入受試者工作(Receiver operating characteristic，ROC)曲線并計(jì)算曲線下面積(Area under curve，AUC)。AUC評(píng)價(jià)是一種閾值獨(dú)立的評(píng)價(jià)模型表現(xiàn)力的方法，其取值范圍一般為0.5到1.0，其評(píng)判標(biāo)準(zhǔn)為：AUC在0.5～0.6為不及格；0.6～0.7為較差；0.7～0.8為一般；0.8～0.9為良好；0.9～1為優(yōu)秀[34]。AUC值越接近1，環(huán)境因子與預(yù)測(cè)的分布范圍相關(guān)性越大。在環(huán)境參數(shù)設(shè)置中選擇刀切法(Jackknife)、繪制響應(yīng)曲線(Response curves)、制作預(yù)測(cè)圖(Make picture of predictions)，輸出結(jié)果為Logistic格式的連續(xù)柵格數(shù)據(jù)，模型迭代次數(shù)為5 000次，置10次重復(fù)，其余參數(shù)為軟件默認(rèn)值。對(duì)于10次重復(fù)結(jié)果，取其平均值作為最終模擬結(jié)果[34]。

1.4 子午沙鼠適生區(qū)劃分及其面積變化

本研究將子午沙鼠基準(zhǔn)氣候和未來潛在分布區(qū)域劃分為4個(gè)等級(jí)：<0.1為非適生區(qū)；0.1～0.3為低適生區(qū)；0.3～0.5為中適生區(qū)；0.5～0.8為高適生區(qū)。通過分別計(jì)算各等級(jí)對(duì)應(yīng)分布區(qū)面積比例，以明確氣候變化對(duì)子午沙鼠實(shí)際分布的影響。另外為了對(duì)比子午沙鼠在基準(zhǔn)氣候與未來氣候情景下分布格局差異，本研究通過基于MaxEnt軟件自動(dòng)生成平等的訓(xùn)練敏感性和特異性(Equal training sensitivity and specificity)閾值(平均值為0.15)將邏輯概率圖轉(zhuǎn)化為二進(jìn)制圖(僅指示適生區(qū)與非適生區(qū)兩種狀態(tài))[35]。在ArcGIS軟件中SDM工具箱中的“Distribution changes between binary”工具中操作。具體步驟如下：輸入基準(zhǔn)氣候情景下子午沙鼠的SDM層命名為“gerbil_thresholded”，輸入2050 s物種的SDM層并命名為“gerbil_2050_thresholded”，建立兩個(gè)柵格文件之間的連接；選擇輸出層的位置，選擇合適的等面積投影，文檔的輸出格式為TIFF格式(.tif)。

2 結(jié)果與分析

2.1 MaxEnt模型準(zhǔn)確性評(píng)估

采用受試者工作曲線ROC分析法對(duì)MaxEnt模型的預(yù)測(cè)結(jié)果進(jìn)行精度驗(yàn)證，在本次模型對(duì)不同氣候情景下子午沙鼠預(yù)測(cè)結(jié)果的平均AUC值均大于0.83(表2)，模型表現(xiàn)力良好，故滿足精度要求，并且模型在10次重復(fù)的穩(wěn)定性相對(duì)較好。因此，在本研究中MaxEnt模型預(yù)測(cè)的不同氣候情景下子午沙鼠的分布范圍具有較高的準(zhǔn)確性和可靠性。

表2 MaxEnt模型在不同氣候下的AUC值Table 2 AUC value of MaxEnt model under different climates

2.2 環(huán)境變量貢獻(xiàn)率分析

本研究結(jié)果表明，8種生物氣候變量對(duì)子午沙鼠的分布有不同的貢獻(xiàn)率(表3)。其中年降水量(Bio12)和年平均氣溫(Bio1)的相對(duì)貢獻(xiàn)率較高，分別為65.50%和15.00%；溫度年較差(Bio7)、最干季度降水量(Bio17)和降水季節(jié)性(Bio15)的相對(duì)貢獻(xiàn)率次之，分別為5.10%，5.10%和5.00%；最冷月最低溫(Bio6)，晝夜溫差月均值(Bio2)和最干季度均溫(Bio9)的相對(duì)貢獻(xiàn)率較小，分別為1.70%，1.40%和1.20%?；诘肚蟹z驗(yàn)結(jié)果可知(圖2)：8個(gè)生物氣候變量對(duì)子午沙鼠物種構(gòu)建模型的貢獻(xiàn)率存在差異。

圖2 MaxEnt模型刀切法檢驗(yàn)結(jié)果Fig.2 Results of Jackknife test in MaxEnt model

表3 生物氣候變量因子貢獻(xiàn)率Table 3 Relative contribution of bioclimatic variables

氣候因子與物種存在概率之間的響應(yīng)曲線可以表明環(huán)境與物種之間的關(guān)系。由圖3可知，隨著年平均氣溫的升高，子午沙鼠的存在概率先升高后降低，年平均氣溫在2.75～10.55℃時(shí)，子午沙鼠存在概率達(dá)50%以上，其中，年平均氣溫為5.10～8.10℃時(shí)，子午沙鼠存在概率達(dá)60%以上；年平均氣溫低于—1.10℃或高于14.30℃時(shí)，子午沙鼠存在概率低于20%。隨著年降水量的升高，子午沙鼠的存在概率整體呈下降趨勢(shì)，年降水量在229.60 mm以下時(shí)，子午沙鼠的存在概率達(dá)50%以上，其中，年降水量在24.40 mm以下時(shí)，子午沙鼠的存在概率達(dá)60%以上；年降水量在高于369.60 mm時(shí)，子午沙鼠的存在概率低于20%。

圖3 主要生物氣候變量響應(yīng)曲線Fig.3 Response curves of major bioclimatic variables

2.3 子午沙鼠在中國的潛在分布預(yù)測(cè)

MaxEnt運(yùn)行結(jié)果表明，在基準(zhǔn)氣候情景下子午沙鼠的中高適生區(qū)主要分布在新疆塔里木盆地和準(zhǔn)噶爾盆地、甘肅省西北部及河西走廊地區(qū)、青海省柴達(dá)木盆地、寧夏與內(nèi)蒙古巴丹吉林沙漠地區(qū)(圖4)；其中高適生區(qū)面積占比為14.20%，中適生區(qū)面積占10.01%；與基準(zhǔn)氣候相比，2050時(shí)段下子午沙鼠高適生區(qū)面積均呈現(xiàn)縮減的趨勢(shì)。其中在RCP2.6氣候情景下，子午沙鼠高適生區(qū)占比縮減到1.55%，中適生區(qū)面積增加到18.52%；在RCP8.5氣候情景下其高適生區(qū)面積占2.30%，中適生區(qū)面積占18.17%(表4)。

表4 不同氣候情景子午沙鼠各適生區(qū)面積占中國面積百分比Table 4 Percentage of suitable areas for midday gerbil in different climate scenarios in China

圖4 不同氣候情景下子午沙鼠在中國的潛在分布Fig.4 The potential distribution of midday gerbil in China under different climate scenarios

2.4 未來氣候情景下子午沙鼠分布范圍面積變化

通過對(duì)基準(zhǔn)氣候和未來氣候下子午沙鼠的分布區(qū)域進(jìn)行對(duì)比發(fā)現(xiàn)：在2050RCP2.6氣候情景下，適生區(qū)面積預(yù)計(jì)增加了3.50×105km2，而適生區(qū)預(yù)計(jì)縮減1.13×106km2；未變化的適生區(qū)面積預(yù)計(jì)在4.65×105km2，在2050RCP8.5氣候情景下，適生區(qū)面積預(yù)計(jì)將增加1.74×105km2，適生區(qū)面積預(yù)計(jì)將縮減1.59×106km2，未發(fā)生變化的適生區(qū)面積預(yù)計(jì)為1.82×102km2(表5，圖5)。

表5 未來氣候情景下子午沙鼠分布范圍面積變化Table 5 Changes in the distribution range of midday gerbil under future climate scenarios

3 討論

3.1 模型的合理性

MaxEnt模型是評(píng)估物種潛在分布的一種較好的選擇[36]，該模型具備對(duì)環(huán)境變量共線性的敏感性較低、位點(diǎn)較少時(shí)仍然可以保持穩(wěn)健、能夠擬合復(fù)雜的變量關(guān)系等優(yōu)勢(shì)[37]。因此，近年來也常常廣泛應(yīng)用于保護(hù)生物學(xué)、入侵生物學(xué)以及其他涉及物種分布和遺傳多樣性的全球氣候變化的研究[38]。特別是在保護(hù)生物學(xué)相關(guān)的研究中，該模型被廣泛的應(yīng)用于評(píng)估物種棲息地適宜性，研究物種的生態(tài)位漂移，并預(yù)測(cè)全球氣候變化對(duì)特定物種地理分布的影響[39]。本研究通過參考以往文獻(xiàn)中記載的參數(shù)模型并結(jié)合子午沙鼠的實(shí)際分布區(qū)，對(duì)模型參數(shù)進(jìn)行了設(shè)置，結(jié)果發(fā)現(xiàn)子午沙鼠MaxEnt模型的測(cè)試集和訓(xùn)練集的AUC值均達(dá)到0.83以上，表現(xiàn)力良好，表明該模型具有較高的準(zhǔn)確值和可靠性。然而，由于物種分布的GPS位點(diǎn)、數(shù)據(jù)集的分辨率以及使用的閾值、運(yùn)行次數(shù)、軟件參數(shù)的設(shè)置、環(huán)境預(yù)測(cè)因子等因素容易導(dǎo)致過度擬合[40]。因此在建模之前，必須剔除相關(guān)性較高的生物氣候變量，同時(shí)對(duì)參數(shù)設(shè)置進(jìn)行優(yōu)化。本研究為解決這一問題，首先通過ENMTools工具初步篩選了生物氣候變量，其次剔除了相同的物種分布點(diǎn)和貢獻(xiàn)率為0的生物氣候變量，并且嘗試了不同的參數(shù)設(shè)置，最終確定了最優(yōu)的模型參數(shù)和相關(guān)性較低的8個(gè)生物氣候變量。

此外，預(yù)測(cè)潛在分布區(qū)范圍的目的是為了給管理者和研究學(xué)者提供更為精確、更有價(jià)值的理論數(shù)據(jù)。但模型預(yù)測(cè)所考慮的環(huán)境因子帶有一定的局限性，不能包含實(shí)際物種生境的所有影響因子及其他控制因子[41]。本研究中未考慮人類干擾、植被、地形因子、天敵影響以及食物資源的競(jìng)爭(zhēng)等因素，這些因素對(duì)子午沙鼠分布范圍的變化也可能會(huì)有不同程度的影響，盡管本研究基于當(dāng)前較為科學(xué)的模型，利用溫度、降雨指標(biāo)得到了子午沙鼠未來分布的預(yù)測(cè)結(jié)果，但我們的預(yù)測(cè)區(qū)域覆蓋地域同《世界自然保護(hù)聯(lián)盟瀕危物種紅色名錄》中記錄的子午沙鼠分布區(qū)相比[45]，一致性僅為68.43%，結(jié)果準(zhǔn)確性還有待進(jìn)一步提升。因?yàn)橛绊懳锓N分布的因素各不相同，影響子午沙鼠的因素也有很多，植被類型、蓋度，地表基質(zhì)、地形和人為干擾等因素影響其棲息地選擇[40-42]，在后續(xù)的數(shù)據(jù)分析中，應(yīng)克服部分指標(biāo)量化困難、數(shù)據(jù)難以收集等問題，納入更多影響物種分布的因素，來提高物種分布模擬的準(zhǔn)確性。

3.2 影響子午沙鼠潛在分布的環(huán)境因子

本研究結(jié)果表明，影響子午沙鼠潛在分布的主要生物氣候變量為年平均氣溫和年降水量，這與武曉東等[46]的研究結(jié)果一致：當(dāng)物種棲息地周圍環(huán)境溫度過高且不能滿足其生物學(xué)和生態(tài)學(xué)特性要求時(shí)，會(huì)導(dǎo)致物種的擴(kuò)散、物種數(shù)量的減少甚至滅絕。雖然子午沙鼠主要分布于荒漠、半荒漠地區(qū)，但溫度大幅升高對(duì)其分布的影響較大。在基準(zhǔn)氣候情景下，子午沙鼠的適生區(qū)主要是在準(zhǔn)噶爾盆地、塔里木盆地、柴達(dá)木盆地、河西走廊以及巴丹吉林沙漠及其周圍的荒漠地區(qū)，但在未來氣候情景下，其分布范圍大量縮減，主要分布至周圍的阿爾泰山脈、天山山脈、昆侖山脈、青藏高原、祁連山脈、內(nèi)蒙古高原、黃土高原、太行山脈等地區(qū)。未來子午沙鼠的物種分布范圍發(fā)生如此變化一方面是因?yàn)樯郎睾徒邓黾雍蟾淖兞嘶哪桶牖哪貐^(qū)的植被結(jié)構(gòu)[47-48]，從而影響子午沙鼠的食物資源的分布；另一方面是由于溫度和降水增加造成了與其同域分布的具有競(jìng)爭(zhēng)關(guān)系的物種數(shù)量增加[49]，使其向周圍競(jìng)爭(zhēng)較小的地區(qū)遷移擴(kuò)散。

3.3 子午沙鼠物種管理的啟示

氣候變化影響子午沙鼠未來分布格局的變化，預(yù)計(jì)在未來氣候變化情景下其地理分布范圍進(jìn)一步縮減，作為荒漠區(qū)的優(yōu)勢(shì)物種，這對(duì)未來荒漠地區(qū)的生態(tài)安全和公共衛(wèi)生安全等均產(chǎn)生重要影響。已有研究證實(shí)，荒漠嚙齒類動(dòng)物對(duì)植被具有顯著的正相關(guān)效應(yīng)，即“沃島效應(yīng)”[50]，子午沙鼠的活動(dòng)促進(jìn)了荒漠植被對(duì)營養(yǎng)物質(zhì)的吸收，進(jìn)而加速了植被群落的更新和演替，當(dāng)子午沙鼠的分布范圍縮減后，這種效應(yīng)也隨之減小，加之人類活動(dòng)等因素的影響，會(huì)加劇荒漠區(qū)的沙化；在未來氣候情景下子午沙鼠高適生區(qū)的荒漠草原、農(nóng)田等生態(tài)系統(tǒng)的危害增加，同時(shí)還會(huì)增加鼠疫等疾病的威脅[21]。因此在未來對(duì)子午沙鼠的管理中，應(yīng)重點(diǎn)關(guān)注子午沙鼠未來適生區(qū)的鼠害防控工作，同時(shí)，還應(yīng)該提高對(duì)子午沙鼠物種如何響應(yīng)氣候變化的認(rèn)識(shí)、改進(jìn)評(píng)估有害生物發(fā)生風(fēng)險(xiǎn)的方法。

4 結(jié)論

本研究基于國內(nèi)的子午沙鼠51個(gè)分布記錄區(qū)域，從19個(gè)生物氣候變量中選取8個(gè)變量，通過MaxEnt模型等預(yù)測(cè)其潛在分布范圍，發(fā)現(xiàn)年平均氣溫和年降水量是影響子午沙鼠分布的主要生物環(huán)境變量(總貢獻(xiàn)率達(dá)到80%)，未來氣候情景下子午沙鼠的適宜生境范圍總體呈現(xiàn)縮減趨勢(shì)，未來氣候變暖、降水增加會(huì)顯著影響子午沙鼠的分布范圍，建議今后子午沙鼠的防控策略和生態(tài)系統(tǒng)管理應(yīng)考慮氣候變化因素。