摘要：研究高原鼠兔的地理分布及其對氣候變化的響應，對進一步了解高原鼠兔在生態(tài)系統(tǒng)中作用意義重大。本研究利用青海省果洛州高原鼠兔存在點和環(huán)境因子數(shù)據(jù)，采用物種分布模型預測在3種排放背景下，當前與未來高原鼠兔宜生分布區(qū)。結果表明：極度梯度提升樹（eXtreme Gradient Boosting，XGBoost）的ROC和TSS值最高，模型表現(xiàn)最佳；最濕季降水量對高原鼠兔分布的影響最大，其次為氣溫年較差和最濕月降水量，且高原鼠兔當前生境適宜度隨海拔和降水量的增加而降低，而隨氣溫的增加呈先增加后降低的趨勢；當前氣候條件下高原鼠兔主要分布在果洛州西北部，未來氣候影響下高原鼠兔適宜生境面積逐漸減少，且適宜生境集中在瑪多縣西北部；在3種排放情景下，高原鼠兔適宜度隨海拔和降雨量的增加而降低，而隨著氣溫的升高，低排放背景下高原鼠兔適宜度增加，高排放背景下其適宜度降低。

關鍵詞：高寒草甸；高原鼠兔；物種分布模型；空間分布；氣候因子

中圖分類號：S812.6""" 文獻標識碼：A"""" 文章編號：1007-0435（2024）06-1902-11

The Distribution of Plateau Pika （Ochotona Curzoniae） in Guoluo Prefecture，

Qinghai Province and Its Response to Climate Change

YANG Li-yan1， SONG Mei-ling1，2， WANG Yu-qin1，2， WANG Hong-sheng1，2， ZHOU Rui1，2*

（1. College of Animal Husbandry and Veterinary Science， Qinghai University， Xining，Qinghai Province 810016， Qinghai;

2. State Key Laboratory of Sanjiangyuan Ecology and Plateau Agriculture and Animal Husbandry， Provincial-Ministerial

Co-construction， Qinghai University， Xining， Qinghai Province 810016， Qinghai）

Abstract：Studying the geographical distribution of plateau pikas and their responses to climate change is of great guiding significance for further understanding its role in the ecosystem. In this study，based on the information of existence points and environmental factor data of plateau pikas in Guoluo Tibetan Autonomous Prefecture of Qinghai Province，distribution models were used to predict the distribution range of plateau pikas in Guoluo Prefecture at present and future under three emission scenarios. The results showed that eXtreme Gradient Boosting （XGBoost） had the highest ROC and TSS values. The climate factors that significantly affect the distribution of plateau pika are the wettest season precipitation，followed by annual temperature range and wettest month precipitation. Moreover，the current habitat suitability of plateau pika decreases with the increase of altitude and precipitation，and shows a trend of first increasing and then decreasing with the increase of temperature;Under the current climate conditions，the plateau pika is mainly distributed in the northwest of Guoluo Prefecture. In the future，under the influence of climate，the suitable habitat area for plateau pika will gradually decrease，while the suitable habitats will gradually concentrate in the northwest of Maduo County;In three greenhouse gas emission scenarios，the suitability of plateau pika decreases with increasing altitude and rainfall，while in the low emission scenario，the suitability of plateau pika increases with increasing temperature. In the medium to high emission scenario，the suitability of plateau pika decreases with increasing temperature.

Key words：Alpine meadow;Plateau pika;Species distribution model;Spatial distribution;Climatic factors

動物空間分布格局反映了個體聚集、隨機或離散的位置，是生態(tài)系統(tǒng)結構的重要組成部分［1］。了解動物種群的空間分布有助于理解動物習性與環(huán)境因子的疊加效應，對物種保護具有重要價值［2］。環(huán)境因子作為生態(tài)系統(tǒng)中最重要的因素之一，不僅影響著生物的生存、生長和繁殖，還影響著整個生態(tài)系統(tǒng)的結構和功能［3］。近年來，物種分布模型（Species distribution models，SDMs）廣泛應用于氣候變化下預測物種的潛在分布［4］，通過對物種和環(huán)境的相關數(shù)據(jù)分析，建立一個以統(tǒng)計或理論為基礎的實證模型，是探討全球變化下物種-環(huán)境相互關系的重要手段之一［5］。近十幾年來，由于處理數(shù)據(jù)算法與生態(tài)位物種分布模型的發(fā)展，學者們已經建立了數(shù)十個不同的分布預測模型，它們的應用領域和理論方法都不盡相同，但都對最大熵模型（Maxent）給予了極高的評價，認為Maxent操作簡單，采用存在點進行預測，可以在同一時間處理大量的數(shù)據(jù)［6］，但不存在點未納入該模型。Thuiller等［7］指出，隨著輸入數(shù)據(jù)的變化，物種分布模式的穩(wěn)定性會下降，只用單一模型進行分布預測的方法并不穩(wěn)定可靠。目前物種分布研究使用的物種分布模型也逐漸由單一模型轉向組合模型［8］。2003年BIOMOD的出現(xiàn)［9］，得到了許多生物學家的認可并被多次采用，BIOMOD是R語言中的一個軟件包，可以利用存在點和不存在點來預測物種的分布，解決許多模型的不確定性以及檢驗物種-環(huán)境的相互關系。Aguirre-Gutierrez等［10］比較了BIOMOD中Maxent的預測準確性以及5種模型的合成效果，結果表明，后者比Maxent模型在預測準確性與合成效果都稍微高一些，Biomod2幾乎涵蓋了Maxent模型，更精確地劃分出該屬的分布區(qū)與未分布的分布區(qū)，適用于對不同物種的潛在分布進行比較，推導出制約該地理分布的主要因素。目前該模型的最新版本為BIOMOD2軟件包（Version：4.2-4），其中有11種物種分布模型算法。

青藏高原被稱為“地球第三極”，具有“亞洲水塔”之稱，是世界上最高的獨立地貌單元，也是生物多樣性的重要保護地，在維持氣候系統(tǒng)穩(wěn)定、生物多樣性保護以及水資源供給等方面，起重要的生態(tài)安全和穩(wěn)定屏障的作用［11-13］。在氣候變化、人類活動等因素的影響下，青藏高原退化的高寒草甸面積已達到1.62×105 km2，主要表現(xiàn)為草地生產力和覆蓋度下降、物種豐富度減少、土壤斑塊化和毒雜草比例增加等方面［14］。草地退化增加草原鼠類棲息地適合度、存活率以及繁殖率，使其種群密度逐步增大，從而使草地退化加?。?5］。高原鼠兔（Ochotona curzoniae）是青藏高原生態(tài)系統(tǒng)中的重要物種之一，主要通過采食和掘洞等行為對草地植被和土壤造成破壞［16］，種群密度過大可對草地生產力和畜牧業(yè)的發(fā)展造成不可估量的損失［17］。由于草地鼠害具有危害面積廣、暴發(fā)速度快和治理難度大等特點，持續(xù)時間長且不易消除，對我國草原生態(tài)安全及牧區(qū)糧食安全構成重大威脅［18］。因此，掌握高原鼠兔等草原鼠類的空間分布及影響其棲息地選擇的因子從而控制草地害鼠種群數(shù)量，對維護草地生態(tài)系統(tǒng)安全有重要意義，也對草原畜牧業(yè)的可持續(xù)發(fā)展意義重大。

本研究利用R語言的BIOMOD2軟件包，通過實地調查果洛藏族自治州高原鼠兔存在點和不存在點，利用11種物種分布模型分別預測果洛州高原鼠兔當前分布圍，采用不同溫室氣體排放濃度和社會發(fā)展水平下（SSP126、SSP245和SSP585）的氣候預測結果作為從低到高3個等級的氣候情景代表未來全球氣候變暖的3 種趨勢，進一步預測2050 s、2070 s和2090 s時間段高原鼠兔分布范圍變化。以期為青藏高原鼠害的科學防治和草原生態(tài)保護提供理論基礎。

1 材料與方法

1.1 研究區(qū)選擇

研究區(qū)域選擇在果洛藏族自治州，位于青藏高原的東部，青海省東南部，東經96°54′～101°51′，北緯32°31′～35°37′。年平均氣溫為-4℃，冷季最低月份歷年平均氣溫為-12.1℃，低限氣溫達到-48.1℃，暖季最熱月份歷年平均氣溫為9℃度，極限高溫為28.1℃，年降水量在400～760 mm之間。果洛藏族自治州平均海拔4 200米以上，植物優(yōu)勢種為矮生嵩草（Kobresia humilis）、線葉嵩草（Kobresia capillifolia）、異針茅（Stipa aliena）、虎耳草（Saxifraga stolonifera）等;主要伴生種有唐松草（Thalictrum aquilegifolium）、毛茛（Ranunculus japonicus）、乳白香青（Anaphalis lactea）等，優(yōu)質的牧草有華扁穗草（Blysmus sinocompressus）、線葉嵩草（K. capillifolia）、垂穗披堿草（Elymus nutans）等20余種。優(yōu)勢鼠種為高原鼠兔。

1.2 研究內容及方法

2022年3月對果洛藏族自治州的高原鼠兔分布區(qū)使用樣帶法進行調查，有高原鼠兔鼠洞即為一個存在點，反之為不存在點。為了消除分布點之間的空間自相關，調查點之間相距至少5 km。采用Excel 2021進行數(shù)據(jù)錄入和初步分析，利用R語言BIOMOD2軟件包的11種模型分別預測果洛州高原鼠兔當前分布范圍，依據(jù)ROC和TSS值選取最佳預測結果，將結果導入Arc GIS10.8 軟件，通過Conversion工具將導入的原始數(shù)據(jù)轉換為柵格數(shù)據(jù)，再利用Reclassify工具，以模型生成的靈敏度和邏輯閾值作為高原鼠兔不同適宜生境的劃分閾值對物種的宜生區(qū)等級進行劃分［19］。采用不同溫室氣體排放濃度和社會發(fā)展水平下的SSP126、SSP245和SSP585氣候預測結果作為從低到高3個等級的氣候情景代表未來全球氣候變暖的3 種趨勢，進一步預測2050 s、2070 s和2090 s時間段高原鼠兔分布范圍變化，利用ArcGIS 10.8對模擬的高原鼠兔當前和未來分布范圍進行地圖繪制。

1.3 高原鼠兔分布點調查

2022年3月對果洛藏族自治州地區(qū)的鼠兔分布區(qū)進行了實地考察。由于研究區(qū)域內不同地區(qū)的地理位置差別很大，因此高原鼠兔調查點的合理布設會受限于地形、距離等多種因素。為提高調查的效率和準確性，本實驗使用地面樣帶法調查，有高原鼠兔鼠洞即為一個存在點，反之為不存在點，調查點之間相距至少5 km，共取得590個存在點和446個不存在點的數(shù)據(jù)（圖1），對應的空間位置由全球定位系統(tǒng)記錄。

1.4 環(huán)境變量選擇及其數(shù)據(jù)獲取

本研究選取的環(huán)境變量包括地形、土壤、氣候3類數(shù)據(jù)，涉及3個地形變量、3個土壤理化性質變量、19個氣候變量，共計25個環(huán)境變量。

數(shù)字高程模型（Digital elevation model，DEM）數(shù)據(jù)來自地理空間數(shù)據(jù)云網站（http：//www. gscloud. cn）、土壤數(shù)據(jù)選自中國科學院資源環(huán)境數(shù)據(jù)云平臺（http：//www.resdc. cn）、未來氣候數(shù)據(jù)來源于聯(lián)合國政府間氣候變化專門委員會（IPCC）發(fā)布的第5次排放情景數(shù)據(jù)（http：∥www. worldclim.org）以及北京氣候中心獨立開發(fā)的大氣環(huán)流模GCM（Beijing climate center climates system model version1.1，BCC CSM1.1）。

氣候數(shù)據(jù)采用世界氣候數(shù)據(jù)庫（WORD-CLIM，version 1.4，http：//www. Worldclim.org），數(shù)據(jù)中包含19個環(huán)境變量（表3），空間分辨率為1 km［20-21］。本研究選取2050s（2041—2060年）、2070s（2061—2080年）和2090s（2081—2100年）3個時間段以及3種溫室氣體排放濃度和社會發(fā)展水平下（SSP126、SSP245和SSP585）的氣候預測的數(shù)據(jù)參與建模，空間分辨率為1 km。

1.5 BIOMOD2物種分布模型及特點

BIOMOD2是一種基于 R語言的新的計算架構，其中包含11種物種分布模型［22］，利用不同的模型方法來分析物種-環(huán)境之間的關系，本研究分別用11種模型來預測高原鼠兔的分布區(qū)，其中11種模型簡介如下（表1）。

1.6 最優(yōu)模型選取

為了對模型精度進行評價，以70%實地調查樣點為訓練樣本、30%實地調查樣點為檢驗樣本。重復運行10次并隨機選取運行過程中的訓練值和測試值，采用受試者工作特征曲線（Receiver operator characteristic curve，ROC）與真實技巧統(tǒng)計值（True skill statistic，TSS）來評價模型精度。TSS的取值范圍一般在0～1，ROC的取值范圍一般在0.5～1（表2），在運行過程中，當TSS與ROC的值越接近于1時，表明模型的模擬效果好［34］。

1.7 數(shù)據(jù)處理

采用Excel 2021進行數(shù)據(jù)錄入和初步分析，使用R語言中的BIOMOD2包對高原鼠兔分布范圍進行模擬，利用GraphPad Prism8繪制箱線圖，以模型生成的靈敏度和邏輯閾值作為高原鼠兔不同適宜生境的劃分閾值，最終將高原鼠兔生境劃分為適宜生境（適宜度指數(shù)≥0.22）、中適宜生境（0.185≤適宜度指數(shù)lt;0.22）、低適宜生境（0.15≤分類閾值lt;0.185）和不適宜生境（適宜度指數(shù)lt;0.15）。依據(jù)分類閾值利用ArcGIS 10.8對模擬的高原鼠兔當前和未來分布范圍進行地圖繪制，利用Origin2022繪制當前與未來氣候條件下果洛州高原鼠兔適宜區(qū)分布變化的趨勢圖。

2 結果與分析

2.1 不同模型精度對比

在11種模型運行中，ANN、GAM運行失敗，只有9種模型運行成功。其中，XGBoost的擬合效果最佳，ROC和TSS分別為0.988和0.874，這說明，XGBoost在10次運行的過程中能較準確地對當前高原鼠兔分布進行預測，并且具有一定的穩(wěn)定性。模擬效果最差為SRE，ROC和TSS分別為0.509和0.017，其基本上不能準確預測高原鼠兔的分布區(qū)范圍，其他模型模擬效果均一般（圖2）。

2.2 影響高原鼠兔分布的環(huán)境因子

依據(jù)運行10次XGBoost模型的結果可以得到影響高原鼠兔分布的不同環(huán)境因子重要值（表3）。其中排名前五的因子依次是最濕季降水量（bio16）、氣溫年較差（bio7）、最濕月降水量（bio13）、最干月降水量（bio14）和降水的季節(jié)性（bio15），因此將這5個因子納入最終模型中。

2.3 當前氣候條件下高原鼠兔潛在適宜分布區(qū)

基于XGBoost最優(yōu)模型重復運行了10次，從中選取效果最好的結果，利用最優(yōu)模型生成的適宜度指數(shù)矢量數(shù)據(jù)進行繪圖，得到果洛州在當前氣候條件下高原鼠兔適宜區(qū)分布情況（圖3）?？傮w來看果洛州高原鼠兔在西北部分布多，在東南部分布少。果洛州高原鼠兔適宜區(qū)面積占總面積的48.92%，低適宜生境的面積占總面積的6.96%，占比面積最小；適宜生境中瑪多縣占比最高為84.57%，其次為瑪沁縣、久治縣、甘德縣和達日縣，分別為38.86%，31.81%，24.53%和11.67%，占比最低的班瑪縣為3.30%；不適宜生境占比最高為班瑪縣（35.57%），其次為久治縣、瑪多縣、達日縣和瑪沁縣，分別為7.67%，7.63%，2.70%和0.94%，占比最低的甘德縣（0.21%）；果洛州高原鼠兔中適宜區(qū)面積中甘德縣占比最大（67.06%），瑪多縣占比面積最小（4.67%）；果洛州高原鼠兔低適宜區(qū)面積中久治縣占比最大（18.95%），瑪多縣面積占比最?。?.13%）（表4）。

2.4 未來氣候條件下高原鼠兔適宜區(qū)分布變化

本研究通過預測3種排放情景下2050s、2070s和2090s三個時間段果洛州高原鼠兔的適宜生境的變化趨勢，可以看出隨著溫室氣體排放濃度增加和未來時間的變化，高原鼠兔適宜生境逐漸縮小，中低適宜生境和不適宜生境逐漸擴大（圖4）。其中，到21世紀50年代，與當前分布預測相比，在ssp126排放情境下，高原鼠兔適宜區(qū)面積增加7.04%，在ssp245、ssp585排放情景下，高原鼠兔適宜區(qū)面積分別減少24.35%，22.28%，中低適宜生境和不適宜生境均有所增大。到21世紀70年代，隨3種排放情景的變化，高原鼠兔適宜區(qū)面積與當前分布預測相比分別減少18.33%，30.44%和34.78%，中低適宜生境顯著增大（10.04%，23.29%和25.87%），而低適宜和不適宜生境比例無明顯變化。到21世紀90年代，隨3種排放情景的變化，高原鼠兔適宜區(qū)面積與當前氣候潛在分布區(qū)相比分別減少20.05%、26.39%和45.3%，低適宜和不適宜生境比例增大，而中適宜生境比例變化幅度不大（表5）。

2.5 當前和未來氣候條件下高原鼠兔適宜度與環(huán)境因子的關系

本研究利用影響高原鼠兔分布的未來環(huán)境因素與高原鼠兔適宜度建立擬合曲線，結果表明，在當前氣候條件下，高原鼠兔適宜度隨海拔的升高呈現(xiàn)先升高后下降的趨勢，且高原鼠兔最優(yōu)適宜度在海拔為4 360 m的區(qū)域（圖5）。在3種排放情景下，未來2050 s、2070 s和2090 s三個時間段高原鼠兔的適宜度變化趨勢也呈現(xiàn)先增高后降低的趨勢，且2050 s時間段內高原鼠兔的適宜度指數(shù)在3種排放情景下均高于2070 s和2090 s這2個時間段。綜合比較當前和未來氣候條件下高原鼠兔適宜度隨海拔的變化趨勢，本研究發(fā)現(xiàn)高原鼠兔最優(yōu)適宜度隨排放濃度和時間的增大而降低。

高原鼠兔適宜度隨降水量的增加呈現(xiàn)下降的趨勢，且3種排放濃度下，當前氣候條件高原鼠兔適宜度均高于未來三個時間段。高原鼠兔對溫度的響應依然顯著，在低排放濃度（SSP126）下，當前氣候條件下高原鼠兔適宜度呈現(xiàn)先增加后降低的趨勢，并在溫度為19.1℃時適宜度最大，而未來三個時間段模擬的氣候條件下高原鼠兔適宜度隨氣溫的升高增大，且2050s、2070s和2090s三個時間段達到高原鼠兔最優(yōu)適宜度的溫度分別為19.7℃，20.3℃和19.1℃；在中等排放濃度（SSP245）下，當前和未來氣候條件中，高原鼠兔適宜度均呈現(xiàn)先增加后降低的趨勢，但當前氣候條件下高原鼠兔最優(yōu)適宜度所在的溫度（19.1℃）遠高于2050s（16.5℃）、2070s（16.8℃）和2090s（17.2℃）三個時間段；在高排放濃度（SSP585）下，也呈現(xiàn)與中等排放濃度下類似的規(guī)律，而未來2050s、2070s和2090s三個時間段高原鼠兔最優(yōu)適宜度達到的溫度分別為16.4℃，16.2℃和15.4℃。

3 討論

本研究中的11個模型中，除了GAM和 ANN的運行失敗之外，其他9個模型都能較好地預測高原鼠兔當前的分布范圍。9種模型中XGBoost的模擬效果最佳。花蕊等［34］通過對青藏高原東部瑪曲縣高原鼠兔現(xiàn)今及未來的時空分布的研究發(fā)現(xiàn)，在10個模型中RF模型表現(xiàn)最佳。吳藝楠等［35］使用BIOMOD在青海湖流域開展了高原鼠兔分布的模擬，發(fā)現(xiàn) GBM 和Maxent模型模擬的效果最好。杜嘉星等［36］在黃河源高原鼠兔分布模擬中使用的 10 種獨立模型中，RF模型表現(xiàn)最好。值得注意的是，XGBoost是BIOMOD2包在2023年2月更新后新加入的一種物種分布模型算法，是一款基于梯度提升算法的機器學習功能庫，它具有良好的學習性能和快速的學習能力，受到了廣泛的重視。XGBoost在連續(xù)型特征的處理、數(shù)據(jù)的預排序和分塊存儲等方面有較好的表現(xiàn)，在另一個常用機器學習庫scikit-learn中實現(xiàn)，XGBoost的性能通常會提高10倍［37］，其致力于讓提升樹突破自身計算極限，以達到快速運算、提升性能的目的［38］。本研究中，由于氣候環(huán)境變量和物種分布數(shù)據(jù)龐大，而且環(huán)境因子變量復雜，需要XGBoost這種簡單、高效和預測速度快的算法，有助于處理分類回歸等問題的模型。

氣候的變化表現(xiàn)為降雨和氣溫，并對物種的存活和繁殖有重要的影響。在本研究中，當前氣候條件下，最濕季降水量、氣溫年較差、最濕月降水量、最干月降水量與降水的季節(jié)性對當前高原鼠兔分布有重要影響，排在前5的指標中4個為降雨量指標，說明降雨量是影響高原鼠兔分布的重要因素。郭新磊等［39］從氣候條件方面進行了分析，認為其洞穴密度主要受降雨和輻射的影響，張衛(wèi)國等［40］的研究也表明，高原鼠兔種群數(shù)量與第3季降水量的年際變率有較好的一致性，特別是生長季的降水量對鼠兔種群密度有重要影響。高原鼠兔繁殖需要充足的食物用于哺育幼崽，植被的營養(yǎng)狀況就直接影響繁殖期的長度與幼崽的存活率，而高寒草地的植被生長主要受氣溫和降水的影響［41］。Zhou等［42］對高原鼠兔活動與環(huán)境因素關系的研究也發(fā)現(xiàn)，適度的降水量會增加高原鼠兔的活動量，但較高的降水量會顯著降低高原鼠兔的活動，這主要是由于高原鼠兔是一種高體溫和高新陳代謝的小型哺乳動物，適度的降水會有效減弱由于環(huán)境高溫對高原鼠兔帶來的不利影響，而當降水量較大時，可能會影響高原鼠兔的采食，導致其營養(yǎng)不良，甚至死亡。另外極端降水通過提高植物群落的高度和蓋度，通過改變高原鼠兔棲息地群落結構，不利于高原鼠兔的生存［43］。隨著降雨量增加，高原鼠兔適宜生境呈現(xiàn)降低的趨勢。高原鼠兔生活在海拔較高的地方，并且對溫度比較敏感［44-46］。我們還發(fā)現(xiàn)，當前氣候條件下，高原鼠兔的適宜度隨溫度的升高呈現(xiàn)先增加后降低的趨勢。我們推測這主要是由于在適當?shù)臏囟确秶鷥龋咴笸每梢杂行ЬS持自身的體溫，使其能在惡劣的高原環(huán)境中生存，當環(huán)境溫度過高時，高原鼠兔會減少在洞外活動的時間，以防止體溫過高導致死亡［47］。本研究也發(fā)現(xiàn)，各環(huán)境因子中，海拔和坡度的重要值較低，這說明海拔和坡度這兩個地形因素對高原鼠兔分布的影響并不明顯。而趙健赟研究發(fā)現(xiàn)海拔和坡度均顯著影響高原鼠兔的分布，這主要是由于本研究高原鼠兔存在點的取樣尺度較小，并且各存在點之間海拔差異較小導致［48］。

本研究采用SSP126、SSP245和SSP585等3個不同排放情景下的氣候變化，對果洛州鼠兔生境未來的適宜度進行模擬，隨著排放濃度增加，高原鼠兔適宜生境逐漸縮小，中低適宜區(qū)生境和不適宜生境增大，高原鼠兔是一種對降水、氣溫等環(huán)境因子非常敏感的小型哺乳動物，其地理分布格局可能受到降水和氣溫變化的影響。本研究通過預測2050 s、2070 s和2090 s三個時間段的3種排放情景下，果洛藏族自治州高原鼠兔的適宜生境隨海拔、溫度與降水量的變化趨勢，可看出當溫室氣體排放逐漸增高時，高原鼠兔的適宜度隨海拔的增加先升高后降低，而當海拔lt;5 000 m時，當前氣候背景下高原鼠兔適宜度高于未來3個時間段，當海拔gt;5 000 m時，當前氣候背景下高原鼠兔適宜度低于未來3個時間段。由此我們推測，隨著未來氣候的變化，高原鼠兔逐漸向海拔較高的棲息地遷移，導致其生境適宜度降低，是高原鼠兔響應氣候變暖的一種策略。大氣溫度和降水量在時間尺度上的變化也進一步證實了這一推測，當前氣候背景下高原鼠兔適宜度顯著高于未來三個時間段，且隨著排放濃度的增大，高原鼠兔適宜度下降速度隨降水量的增加逐漸變快。高原鼠兔最佳適宜度的溫度范圍隨著排放濃度的增加呈現(xiàn)逐漸降低的趨勢。為應對全球氣候變化，高原鼠兔會放棄原有的平坦開闊的棲息地，逐漸向高海拔平地或坡度較高的山坡擴散。我們在對其中一處研究區(qū)進行連續(xù)2年的高原鼠兔種群變動調查過程中也初步發(fā)現(xiàn)這一現(xiàn)象。但這一現(xiàn)象的出現(xiàn)也不完全由氣候變化主導，人為的滅鼠活動和放牧也會強烈干擾高原鼠兔的棲息地選擇，導致高原鼠兔遷移至干擾較少的區(qū)域活動。然而，現(xiàn)有研究并沒有充足的證據(jù)來證實我們的這一推測，需要開展進一步的實驗進行驗證。

4 結論

本研究采用地面調查對果洛藏族自治州高原鼠兔分布點進行匯總，并采用BIOMOD2的11個模型對果洛藏族自治州高原鼠兔分布區(qū)進行模擬，其中極度梯度提升樹（eXtreme Gradient Boosting，XGBoost）模型表現(xiàn)最佳，氣候因子中最濕季降水量（bio16）對高原鼠兔生境選擇的影響顯著；當前氣候條件下高原鼠兔主要分布在果洛州西北部，在未來氣候的影響下高原鼠兔適宜生境面積逐漸減少，中低適宜生境和不適宜生境隨之增大，且適宜生境逐漸集中在瑪多縣西北部地區(qū)；在未來時間尺度下，高原鼠兔的適宜度隨海拔和降雨量的增加而降低，在低排放濃度下，高原鼠兔適宜度隨氣溫的增加而升高，在中高排放濃度下，高原鼠兔適宜度隨氣溫的增加而降低。

參考文獻

［1］ 張莉. 中國北方大尺度植被格局研究［D］. 北京：華北電力大學，2011，19-21

［2］ 何倩蕓，王小祎，羅概，等. 不同放牧模式下高原鼠兔洞口的空間分布格局變［J］. 獸類學報，2021，41（4）：431-440

［3］ 祁應蓮，陳志，張慧武，等. 高原鼠兔種群密度與生境因子作用關系研究［J/OL］. 草業(yè)科學，http：//kns.cnki.net/kcms/detail/62.1069.S.20240125.1100.003.html，2024-01-26/2024-02-27

［4］ 許仲林，彭煥華，彭守璋. 物種分布模型的發(fā)展及評價方法［J］. 生態(tài)學報，2015，35（2）：557-567

［5］ 畢雅瓊，張明旭，陳元，等. 基于Biomod2組合模型的中國野生芍藥Paeonia lactiflora適宜生境分布［J］. 中國中藥雜志，2022，47（2）：376-384

［6］ 劉濤，劉玉萍，呂婷，等. 基于Biomod2組合模型預測青藏高原特有屬扇穗茅屬物種的潛在分布［J］. 草地學報，2020，28（6）：1650-1656

［7］ THUILLER W. Biomod：optimizing predictions of species distributions and projecting potential future shifts under global change［J］. Global Change Biology，2003（32）：369-373

［8］ 楊穎倫. 基于MaxEnt模型的滇藏玉蘭適生區(qū)分布預測研究［J］. 現(xiàn)代農業(yè)科技，2024（3）：98-102

［9］ FIELDING A H， BELL J F. A review of methods for the assessment of prediction errorsin conservation presence/ absence models［J］. Environmental Conservation，1997，24（1）：38-49

［10］AGUIRRE-GUTIERREZ J，CARVALHEIRO L G，POLCE C，et al. Fit-for purpose：Species distribution model performance depends on evaluation criteria：Dutch Hover flies as a case study［J］. Plos One，2013，8（5）：e63708

［11］高明龍. 基于Biomod2組合模型的我國山楊潛在分布區(qū)研究［J/OL］. 南京林業(yè)大學學報（自然科學版），http：//kns.cnki.net/kcms/detail/32.1161.S.20230313.0929.002.html，2023-03-13/2023-12-20

［12］曹曉云，周秉榮，周華坤，等. 氣候變化對青藏高原植被生態(tài)系統(tǒng)的影響研究進展［J］. 干旱氣象，2022，40（6）：1068-1080

［13］孫鴻烈，鄭度，姚檀棟，等. 青藏高原國家生態(tài)安全屏障保護與建設［J］. 地理學報，2012，67（1）：3-12

［14］王英成. 三江源區(qū)退化高寒草甸土壤微生物多樣性及空間變化規(guī)律研究［D］. 西寧：青海大學，2021：23-24

［15］YU C Q，ZHANG Y J，CLAUS H，et al. Ecological and environmental issues faced by a developing Tibet［J］. Environmental Science amp; Technology，2012，46：1979-1980

［16］高艷美，吳鵬飛. 高寒草甸退化對土壤昆蟲多樣性的影響［J］. 生態(tài)學報，2016，36（8）：2327-2336

［17］花蕊，周睿，包達爾罕，等. 基于旋翼無人機低空遙感的高原鼠兔危害等級劃分技術研究［J］. 草業(yè)學報，2022，31（4）：165-176

［18］龐曉攀. 高寒草甸植物生產力與土壤主要養(yǎng)分對高原鼠兔干擾響應的研究［D］.蘭州：蘭州大學，2020：67-68

［19］衛(wèi)萬榮. 高原鼠兔和高原鼢鼠種群消長規(guī)律及其與植被關系的研究［D］. 蘭州：蘭州大學，2018：40-41

［20］楚彬，包達爾罕，葉國輝，等. 青藏高原東緣高原鼢鼠（Eospalax baileyi）生境適宜性研究［J］. 中國草地學報，2023，45（8）：100-108

［21］FICK S E，HIJMANS R J. World Clim 2：New 1 km spatial resolution climate surfaces for global land areas. International Journal of Climatology，2017，37：4302-4315

［22］劉炳蓮，苗潤澤，陳旭升，等. 基于組合模型預測吉林東部山區(qū)原麝的適宜分布區(qū)［J］. 四川動物，2023，42（2）：121-130

［23］NELDER J A，WEDDERBURN R W M. Generalized Linea Models［J］. Journal of the Royal Statistical Society，1972，135（3）：370-384

［24］HASTIE T，TIBSHIRANI R. Generalized Additive Models［J］. Statistical Science，1986，1（3）：297-310

［25］BUSBY J R. BIOCLIM：a bioclimate analysis and prediction system［J］. Plant Protection Quarterly，1991，6（4）：115-131

［26］FRIEDMAN J H. Multivariate adaptive regression splines［J］. The Annals of Statistics，1991，19（1）：1-67

［27］HASTIE T，BUJA T A. Flexible Discriminant Analysis by Optimal Scoring［J］. Journal of the American Statistical Association，1994，89（428）：1255-1270

［28］BREIMAN L，F(xiàn)RIEDMAN J H，OLSHEN R A，et al. Classification and Regression Trees［M］. New York：Chapman and Hall，1984，103-105

［29］RIDGEWAY G. The state of boosting［J］. Computing Science and Statistics，1999，31：172-181

［30］FRIEDNAN J，HASTIE T，TIBSHIRANI R. Additive logistic regression：a statistical view of boosting （With discussion and are joinder by the authors） ［J］. The Annalsof Statistics，2000，28（2）：337-407

［31］RIPLEY B D. Pattern Recognition and Neural Networks［M］. Cambridge：Cambridge University Press，2007：1-15

［32］BREIMAN L. Random Forests［J］. Machine Learning，2001，45（1）：5-32

［33］徐凱，張會妨. 基于RFE-BXGBoost的軸承套圈溝道表面缺陷識別方法［J］. 機電工程，2023，40（11）：1691-1699

［34］花蕊，周睿，包達爾罕，等. 瑪曲縣高原鼠兔地理分布預測及其對氣候變化的潛在響應［J］. 草原與草坪，2020，40（3）：1-8

［35］吳藝楠，馬育軍，劉文玲，等. 基于BIOMOD的青海湖流域高原鼠兔分布模擬［J］. 動物學雜志，2017，52（3）： 390-402

［36］杜嘉星，孫義，向波，等. 基于BIOMOD的黃河源區(qū)高原鼠兔潛在分布及其影響因子［J］. 草業(yè)科學，2019，36（4）：1074-1083，917

［37］TIANQI C， CARLOS G. XGBoost： A Scalable Tree Boosting System［J］. CoRR，2016，16（3）：4681-4690

［38］鄧威，梅玉杰，李勇，等. 基于BOXGBoost的配變日峰值負荷預測及重過載預警方法［J/OL］. 電力系統(tǒng)及其自動化學報， https：//doi.org/10.19635/j.cnki.csu-epsa.001295，2023-07-03/2023-12-16

［39］郭新磊，宜樹華，秦彧，等. 基于無人機的青藏高原鼠兔潛在棲息地環(huán)境初步研究［J］. 草業(yè)科學，2017，34（6）： 1306-1313

［40］張衛(wèi)國，丁連生，韓天虎. 降水對高原鼠兔種群消長的影響［J］. 草業(yè)科學，1999，（6）：20-22，25

［41］宋文杰，張海蘭，洛藏昂毛，等. 多年鼠類防控對高原鼠兔種群數(shù)量及植物群落結構的影響［J］. 草地學報，2023，31（9）：2853-2859

［42］ZHOU R， HUA R， TANG Z S， et al. Daily and Seasonal Activity Patterns of Plateau Pikas （Ochotona curzoniae） on the Qinghai-Tibet Plateau， China， and Their Relationship with Weather Condition［J］. Animals， 2023，13（10）： 1689-1701

［43］宜樹華，曹文達，張建國，等. 氣候變化和人類活動背景下高原鼠兔動態(tài)分布預估研究進展［J］. 南通大學學報（自然科學版），2020，19（4）：16-30

［44］馬蘭，格日力. 高原鼠兔低氧適應分子機制的研究進展［J］. 生理科學進展，2007（2）：143-146

［45］葛慶征. 極端降水對高原鼠兔（Ochotona curzoniae）種群動態(tài)的影響［D］. 蘭州：蘭州大學，2013：147-148

［46］馬波，王小明，劉曉慶，等. 高原鼠兔洞穴數(shù)量與其棲息地植被分布格局的GIS分析［J］. 生物多樣性，2011， 19（1）：71-78

［47］WAN D， SUN R， WANG Z， et al. Effects of temperature and photoperiod on thermogenesis in plateau pikas （Ochotona curzoniae） and root voles （Microtus oeconomus）［J］. Journal of comparative physiology，1999，169 （1）：77-83

［48］趙健赟. 基于BIOMOD2集成模型的黃河源高原鼠兔潛在分布與干擾強度分析［J/OL］. 生態(tài)學雜志，http：//kns.cnki.net/kcms/detail/21.1148.Q.20230830.1848.002.html，2023-08-21/2023-11-08

（責任編輯 彭露茜）