喬雪麗,林 峻,吳建國,趙雪迷,巨喜鋒,鄭江華,3,*

1 新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院, 烏魯木齊 830046 2 新疆治蝗滅鼠指揮部辦公室, 烏魯木齊 830001 3 綠洲生態(tài)教育部重點實驗室, 烏魯木齊 830046

大沙鼠(Rhombomysopimus)是亞洲中部荒漠林以及荒漠草原的典型棲居者,棲息地為亞洲大陸部分國家,是家庭集群動物,其洞穴系統(tǒng)的入住率平均持續(xù)11a[1- 3]。大沙鼠生活在地下洞穴系統(tǒng)中的家庭群體中,該系統(tǒng)包括一個復(fù)雜的隧道和小室網(wǎng)絡(luò),深度達數(shù)米。一個典型的家庭群體包括一個成年雄性,一個或多個成年雌性及其未成熟的后代[4]。

大沙鼠的危害是多方面的。近年來,荒漠林鼠害持續(xù)高發(fā),嚴重危害著荒漠林的健康生長,導(dǎo)致生物結(jié)皮層退化,荒漠化加劇,生態(tài)環(huán)境持續(xù)惡化[5-6]。更為嚴重的是,鼠類大量挖掘洞穴,導(dǎo)致地表植物覆蓋率顯著下降,地區(qū)沙漠化顯著加劇。鼠類挖掘洞穴時所留下的土丘受風(fēng)蝕、徑流的影響,不但使荒漠林和荒漠草原出現(xiàn)大片無植被覆蓋地區(qū),而且也成為了沙塵暴沙塵的重要來源[7]；另一方面,大沙鼠也是亞洲沙漠和草原中鼠疫細菌(耶爾森氏菌)的主要宿主[1]。當(dāng)前,鼠疫病原體主要在非洲、美洲和亞洲的嚙齒類動物中傳播,并主要在非洲造成人類死亡[8]。近年來,我國也陸續(xù)發(fā)現(xiàn)多例青海、內(nèi)蒙等地的人間鼠疫病例。因此,準確地監(jiān)測和預(yù)測大沙鼠的空間分布及危害情況,對及時采取有效措施進行防治具有重要意義。

大沙鼠對于棲息生境的選擇是伴隨著自然選擇的長期進化過程而不斷發(fā)展變化的,這種選擇受諸多條件制約,如不同物種間的競爭或相互利用、自然地理因素、人類社會活動的影響等等[9]。影響大沙鼠棲息地選擇的三個最主要因素是食物條件、地形特點和植被覆蓋度,季節(jié)差異對大沙鼠棲息地的選擇也有明顯的影響作用[9]。戴昆、趙天飆等學(xué)者認為大沙鼠主要選擇丘間地及固定沙丘為筑巢地形[10-11]；張三亮等認為大沙鼠的洞群主要分布在固定沙丘,其洞群主要分布在丘間地,其次是迎風(fēng)坡,大沙鼠不會選擇背風(fēng)坡和丘頂營巢,并發(fā)現(xiàn)大沙鼠分布與優(yōu)勢植物種類、植被蓋度有顯著關(guān)系[12]。

氣候變化影響到世界各地的許多生態(tài)系統(tǒng)和生物群,包括目前許多物種的分布[13- 16]。未來的氣候變化會加速物種的分布,根據(jù)生態(tài)位模型(ENMs)預(yù)測物種的潛在分布[17]可以提醒科學(xué)家和決策者未來氣候變化將對生物多樣性構(gòu)成威脅的潛在風(fēng)險,并幫助他們提出積極的應(yīng)對策略以減輕氣候?qū)ι锒鄻有宰兓挠绊慬18]。在各種 ENM 中,最大熵(MaxEnt)算法具有相對較高的預(yù)測精度,同時使用少量數(shù)據(jù)根據(jù)環(huán)境變量確定物種發(fā)生率[19]。楊楠等使用MaxEnt模型對藍馬雞(Crossoptilonauritum)棲息地適宜性進行了分析和評價[20]；Md.Mamun等利用最大熵模型(MaxEnt)對朝鮮半島入侵外來物種大嘴鱸(Micropterussalmoides)的分布模式進行預(yù)測[21]；袁智文等利用MaxEnt模型構(gòu)建了華南梅花鹿棲息地適宜性模型[22]。

目前,國內(nèi)外大沙鼠的研究大多集中在對其局域分布格局、生理生化等方面[23-24],對于亞洲大陸大沙鼠的潛在分布預(yù)測處于空白狀態(tài),本研究通過收集和篩選大沙鼠地理分布信息,結(jié)合相關(guān)環(huán)境數(shù)據(jù),利用MaxEnt模型對在不同環(huán)境條件下大沙鼠在亞洲的潛在地理分布變化進行預(yù)測,旨在研究以下問題:(1)現(xiàn)代氣候條件下大沙鼠在亞洲的潛在地理分布,并以此為依據(jù)探究大沙鼠的潛在地理分布與環(huán)境因子的關(guān)系,探討不同情景下限制大沙鼠潛在地理分布的重要環(huán)境因子；(2)未來氣候變化情景下,大沙鼠在亞洲各個國家潛在分布區(qū)的地理變化。本文對大沙鼠的空間分布進行預(yù)測、預(yù)報,對有效保護森林資源的可持續(xù)發(fā)展和公共衛(wèi)生監(jiān)測方面起著不可比擬的作用,對亞洲國家制定合理的預(yù)防方案具有重要意義。

1 材料和方法

1.1 物種發(fā)生數(shù)據(jù)

圖1 大沙鼠的分布點數(shù)據(jù)Fig.1 Distribution point data of Rhombomys opimus

大沙鼠的發(fā)生數(shù)據(jù)以下列方式獲得:(1) 2016 至 2020 年的實地調(diào)查和新疆治蝗滅鼠指揮部辦公室提供的新疆大沙鼠發(fā)生數(shù)據(jù)點,共計89個。新疆境內(nèi)的采樣點多為實地考察獲得,調(diào)查點主要位于古爾班通古特沙漠荒漠林、和布克賽爾蒙古自治縣、瑪納斯縣等區(qū)域,使用 GPS 記錄大沙鼠的緯度、經(jīng)度和高度,部分數(shù)據(jù)已發(fā)表[23- 25];(2)國內(nèi)外公開發(fā)表的論文[9,12,26- 33],分布數(shù)據(jù)中僅有地名無經(jīng)緯度信息,通過GPSSPG網(wǎng)站(http://www.gpsspg.com/)進行查詢[34];(3)全球生物多樣性信息交換平臺(GBIF,https://www.gbif.org/)?；谏鲜鍪占?07個分布數(shù)據(jù),本文從空間上篩選了大沙鼠的分布數(shù)據(jù),并利用ArcGIS10.4在研究區(qū)域建立了10′×10′網(wǎng)格,以確保每個網(wǎng)格中只有一條記錄。最終得到123個物種發(fā)生數(shù)據(jù)記錄,如圖1。

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

本文初步選擇了39個可能影響大沙鼠空間分布的環(huán)境變量:氣候數(shù)據(jù)來源于全球氣候數(shù)據(jù)庫網(wǎng)站(http://www.worldclim.org/),獲得19個生物氣候變量,空間分辨率為5arc-minutes,IPCC AR5 分析了4種典型溫室氣體排放濃度(RCP2.6、RCP4.5、RCP6.0、RCP8.5)對未來氣候的影響[35-36],本研究選取 BCC-CSM1.1 全球未來氣候模式 2050 年及 2070 年三種氣候情景(RCP2.6、RCP4.5,RCP8.5)作為未來氣候情景；從世界土壤數(shù)據(jù)庫(http://www.iiasa.ac.at/)獲得17個土壤變量；選擇了3個地形變量(http://www.resdc.cn/)(表1)。

環(huán)境變量被重新采樣為30s的空間分辨率,許多環(huán)境變量在空間上是相關(guān)的,這可能導(dǎo)致預(yù)測過度擬合。因此,用斯皮爾曼的相關(guān)系數(shù)來研究環(huán)境變量中因子之間的相關(guān)性。如果系數(shù)為b<0.75,則保留環(huán)境變量;如果系數(shù)為b>0.75,則保留其中生態(tài)意義較強的環(huán)境變量[37]。

1.3 模型預(yù)測及評估

將大沙鼠分布數(shù)據(jù)和不同場景下的環(huán)境數(shù)據(jù)導(dǎo)入MaxEnt模型中(maxent 3.4.1; http://www.cs.princeton.edu/wschapire/maxent),為了確定影響大沙鼠分布的關(guān)鍵環(huán)境因子,隨機選擇75%的大沙鼠分布數(shù)據(jù)作為訓(xùn)練集,其余25%的數(shù)據(jù)作為測試集,共有10000個背景點,進行了10次重復(fù),并對重復(fù)運行進行了交叉驗證,以保證模型的準確性,開啟刀切法(Jackknife)計算環(huán)境變量對大沙鼠分布的影響,最終得到各氣候環(huán)境因子的響應(yīng)曲線。

利用接受者操作特征曲線(receiver operating characteristic curve,ROC)對MaxEnt模型運算結(jié)果進行檢驗,AUC值從0到1不等,AUC值高表示優(yōu)越,預(yù)測效果的表現(xiàn)分為差(0—0.6)、較差(0.6—0.7)、一般(0.7—0.8)、較好(0.8—0.9)和極好(0.9—1)[38]。

1.4 適生等級劃分和主導(dǎo)環(huán)境因子分析

參照IPCC報告有關(guān)評估可能性劃分方法,并結(jié)合大沙鼠實際分布情況,將大沙鼠潛在地理分布劃分為不同等級及相應(yīng)分布范圍,劃分標(biāo)準為:存在概率<0.05 為不適生區(qū); 0.05≤存在概率<0.33 為低適生區(qū); 0.33≤存在概率<0.66 為中適生區(qū); 存在概率≥0.66 為高適生區(qū)[39]。 綜合Jackknife檢驗、貢獻率和置換重要值分析結(jié)果確定影響大沙鼠適生分布的主導(dǎo)環(huán)境變量,根據(jù)大沙鼠潛在分布概率在不同主導(dǎo)環(huán)境變量中的響應(yīng)曲線,獲得不同分布等級下主導(dǎo)環(huán)境因子的貢獻率。

表1 氣候變量名稱及描述

2 結(jié)果與分析

2.1 模型精度評價

基于MaxEnt模型建模重復(fù)運行10次后得到的訓(xùn)練數(shù)據(jù)AUC值在0.984—0.985之間,測試集AUC值在0.969—0.987,見表2。根據(jù)AUC值的評估標(biāo)準[39],模型預(yù)測的準確性達到“極好”,表明模型能很好的擬合大沙鼠物種分布數(shù)據(jù),其預(yù)測結(jié)果可信。

表2 不同氣候場景下AUC及標(biāo)準差

2.2 環(huán)境變量對大沙鼠分布的影響

利用Jackknife檢驗可以評估環(huán)境變量對預(yù)測結(jié)果的影響程度,從而判斷不同變量對大沙鼠潛在分布的重要性。Jackknife檢驗結(jié)果表明如表3所示,影響大沙鼠在適生區(qū)分布最主要的環(huán)境變量為溫度季節(jié)性變化的標(biāo)準差(BIO4)、土壤基本飽和度(T_BS)、最干季度降水量(BIO17)、最暖季度降水量(BIO18)和土壤可交換鈉鹽(T_ESP),五類環(huán)境變量累計貢獻率達64.5%,其余變量貢獻較小。

表3 影響大沙鼠環(huán)境變量重要性的刀切法檢驗

2.3 當(dāng)前氣候下大沙鼠在亞洲大陸的潛在分布區(qū)

圖2 大沙鼠當(dāng)前在亞洲的適生區(qū)分布范圍 Fig.2 The current suitable distribution range of Rhombomys opimus in Asia

大沙鼠在當(dāng)前氣候情境下的高適生區(qū)面積為940830km2,約占研究區(qū)總面積的1.28%；中適生區(qū)面積為3638189km2,占研究區(qū)總面積的4.95%； 低適生區(qū)面積為10594040km2,約占研究區(qū)總面積的14.43%。

高適生區(qū)(0.66—1)主要分布在:(1)中國境內(nèi):新疆北部和東部地區(qū)的荒漠地區(qū)；甘肅省安西的明水、酒泉、嘉峪關(guān)、張掖和民勤等區(qū)域；內(nèi)蒙古阿拉善等的部分荒漠草原地區(qū)；(2)蒙古國的科布多省、南戈壁省和東戈壁省等地區(qū)；哈薩克斯坦國家的巴爾喀什湖的東部部分地區(qū)、曼吉斯套州和克孜勒奧爾達州等部分荒漠地區(qū)；烏茲別克斯坦的中部地區(qū)；土庫曼斯坦的克茲勒阿爾瓦特和達紹古茲州的部分地區(qū)；土耳其的南部地區(qū)；伊拉克的薩拉赫丁??；其次在吉爾吉斯坦伊塞克湖的西南部呈零星分布、塔吉克斯坦南部與阿富汗交界處、阿富汗的北部薩曼甘省、巴格蘭省、昆都士省和塔哈爾省等部分區(qū)域；伊朗也有少量分布區(qū)域,見圖2。

中適生區(qū)(0.33—0.66)和低適生區(qū)(0.05—0.33)同樣分布在以上國家,其中低適生區(qū)增加了俄羅斯的阿加布里亞特自治區(qū)南部與蒙古國交界處和額爾齊斯河?xùn)|部的少量區(qū)域,阿塞拜疆西部和沙特阿拉伯北部邊疆區(qū)。

2.4 未來氣候下大沙鼠在亞洲大陸的潛在分布區(qū)

從表4中可以得到,Rcp2.6場景下2050年和2070年大沙鼠高適生區(qū)面積和中適生區(qū)面積較當(dāng)前氣候場景下呈增大趨勢,其中2050年和2070年高適生區(qū)面積較當(dāng)前分別增大15.78%和8.32%,中適生區(qū)面積和低適生區(qū)均較當(dāng)前氣候場景下呈明顯增大趨勢。Rcp4.5場景下2050和2070年大沙鼠高適生區(qū)面積較當(dāng)前氣候場景下增大15.10%和13.18%；中適生區(qū)面積較當(dāng)前氣候場景下呈增大趨勢,低適生區(qū)呈明顯增大趨勢。Rcp8.5場景下2050年和2070年高適生區(qū)面積較當(dāng)前增加了13.44%和18.18%,中適生區(qū)面積較當(dāng)前呈下降趨勢；低適生區(qū)面積較當(dāng)前呈明顯增大趨勢。對比4種氣候場景(圖3),大沙鼠的高適生區(qū)和中適生區(qū)空間位置偏移程度較小,但低適生區(qū)空間位置偏移程度較大。

表4 不同氣候變化情境下大沙鼠適生區(qū)面積

3 結(jié)論與討論

1)本研究基于大沙鼠的分布點數(shù)據(jù)和39個環(huán)境變量數(shù)據(jù),利用MaxEnt模型對大沙鼠在亞洲大陸的適生區(qū)進行預(yù)測,ROC 曲線即 AUC 值高于0.9,標(biāo)準差在0.015上下浮動,表明模型達到“極好”的標(biāo)準。

2)Jackknife檢驗結(jié)果表明,影響大沙鼠在適生區(qū)分布最主要的環(huán)境變量為溫度季節(jié)性變化的標(biāo)準差(BIO4)、土壤基本飽和度(T_BS)、最干季度降水量(BIO17)、最暖季度降水量(BIO18)和土壤可交換鈉鹽(T_ESP),5類環(huán)境變量累計貢獻率達64.5%。

3)研究結(jié)果表明,氣候變化影響大沙鼠2050年代和2070年代大沙鼠的分布,并且影響會因地區(qū)而異。對比四種氣候場景下,大沙鼠的高適生區(qū)和中適生區(qū)空間位置偏移程度較小,但低適生區(qū)空間位置偏移程度較大。與其他國家相比,中國適生區(qū)的溫度季節(jié)性變化的標(biāo)準差、土壤基本飽和度、最干季度降水量、最暖季度降水量和土壤可交換鈉鹽更適于大沙鼠生存。中國境內(nèi)的適生范圍為36°—49°N,76°—114°E；高適生區(qū)主要分布于西起新疆裕民縣巴魯克山西麓和中(國)哈(薩克斯坦)邊界的阿拉山口一帶,經(jīng)甘肅省到內(nèi)蒙古蘇尼特右旗的賽汗高畢蘇木。新疆大沙鼠適生區(qū)分布范圍變化較大,3 種情景模式下大沙鼠的適生區(qū)位置向新疆北部擴張；甘肅適生區(qū)位置向西北部擴張；內(nèi)蒙西北部和阿拉善地區(qū)大沙鼠的適生區(qū)位置向四周擴張。

圖3 未來 30 年和 50 年不同氣候變化情景下的大沙鼠的潛在分布區(qū)Fig.3 Potential distribution areas of Rhombomys opimus under different climate change scenarios in the next 30 and 50 years PCR為氣候情景

4)當(dāng)前氣候場景下,大沙鼠高適生區(qū)面積為940 830km2,Rcp2.6、Rcp4.5和Rcp8.5三種氣候場景下2050年高適生區(qū)面積較當(dāng)前分別增長15.78%、15.10%和13.44%；Rcp2.6、Rcp4.5和Rcp8.5三種氣候場景下2070年高適生區(qū)面積較當(dāng)前增長8.32%、13.18%和18.18%。大沙鼠高適生區(qū)主要分布在中國、蒙古國、哈薩克斯坦、烏茲別克斯坦、土庫曼斯坦、土耳其和伊拉克,其次在吉爾吉斯坦國、塔吉克斯坦國、阿富汗國、伊拉克和伊朗國也有少量分布區(qū)域。今后,相關(guān)部門應(yīng)制定合理的預(yù)防方案,防止鼠害對森林資源和公共衛(wèi)生的影響進一步擴大。

本研究運用MaxEnt生態(tài)位模型對大沙鼠的潛在分布進行模擬,模型AUC值達到0.9以上,說明對大沙鼠分布預(yù)測具有較好的效果。大沙鼠的采食和掘洞強烈影響荒漠植物的發(fā)育和外貌,以及荒漠生態(tài)系統(tǒng)的結(jié)構(gòu)和功能,嚴重威脅著天然草地生態(tài)環(huán)境的健康,本研究將模擬結(jié)果與GIS相結(jié)合,對大沙鼠的監(jiān)測、風(fēng)險評估具有重要意義。在建立物種分布模型時,環(huán)境變量和樣本量的選擇會對生態(tài)位模型的預(yù)測結(jié)果產(chǎn)生一定的影響[40]。多數(shù)研究者只考慮氣候因子[41-42]對物種分布的影響,本研究將地形、土壤共39個環(huán)境變量加入其中,剔除貢獻率較小的環(huán)境變量,減少了冗余信息對模擬結(jié)果的影響,模擬效果較好。

結(jié)果顯示,大沙鼠在中國的高適生區(qū)面積在不同場景下達7.9×105—8.7×105km2；同時中國的大沙鼠適生區(qū)范圍隨著碳排放量的增加并沒有表現(xiàn)出規(guī)律性,2050年,在Rcp2.6場景下大沙鼠在中國的高適生區(qū)范圍達到最大為8.7×105km2；2070年,在Rcp8.5場景下大沙鼠在中國的高適生區(qū)范圍達到最大為8.2×105km2；大沙鼠在中國的中適生區(qū)面積在不同場景下達14.0×105—15.4×105km2；大沙鼠在中國的低適生區(qū)面積在不同場景下達35.6×105—38.0×105km2；大沙鼠在中國的高適生區(qū)和低適生區(qū)面積較當(dāng)前呈明顯增長趨勢。

大沙鼠在蒙古國的發(fā)生區(qū)面積較中國次之,2070年大沙鼠在蒙古國的高適生區(qū)范圍隨著碳排放量的增加而增長；中適生區(qū)和低適生區(qū)范圍隨著碳排放量的增加而下降；大沙鼠在蒙古國的高適生區(qū)范圍在Rcp8.5場景下達到最大為2.3×105km2；2050年大沙鼠在蒙古國的高適生區(qū)和中適生區(qū)范圍隨著碳排放量的增加并沒有表現(xiàn)出規(guī)律性,大沙鼠在蒙古國的高適生區(qū)范圍在Rcp4.5場景下達到最大值為2.33×105km2；低適生區(qū)范圍隨著碳排放量的增加而下降。

這表明在未來50年,大沙鼠仍對草地生態(tài)環(huán)境造成嚴重危害,大沙鼠在中國和蒙古國的危險系數(shù)極高,有關(guān)部門應(yīng)密切關(guān)注,嚴防鼠害成災(zāi)。

未來氣候變化下,大沙鼠的適生區(qū)范圍波動較大,高適生區(qū)范圍較當(dāng)前氣候模式呈現(xiàn)明顯增大趨勢,7種氣候模式下的高適生區(qū)面積大小排名為:2070Rcp8.5>2050Rcp2.6>2050Rcp4.5>2050Rcp8.5>2070Rcp4.5>2070Rcp2.6>當(dāng)前氣候,表明未來氣候變化越來越適宜大沙鼠的生存,大沙鼠的危害程度會加重,分布地區(qū)擴大。

劉洋等[34]利用最大熵模型MaxEnt分析了埃及吹綿蚧在中國當(dāng)前和未來(2050年和2070年)的適生區(qū)范圍,結(jié)果表明Rcp8.5氣候模式下的適生區(qū)范圍大于Rcp4.5氣候模式下的適生區(qū)范圍,說明溫室氣體的大量排放導(dǎo)致的全球溫度的持續(xù)升高所帶來的環(huán)境變化可能更加適宜于埃及吹綿蚧的生存,該蟲的危害在未來可能加重；趙佳強等[43]基于 MaxEnt 模型,對刺槐葉癭蚊在中國當(dāng)前和未來(2050年) 的適生區(qū)進行預(yù)測,結(jié)果顯示:3 種外排情景 Rcp2.6、Rcp4.5、Rcp8.5的總適生區(qū)均比當(dāng)前氣候的總適生范圍大,以高度、中度適生區(qū)面積的增大為主,新疆和我國北部區(qū)域面積顯著擴增；唐繼洪等[44]結(jié)合政府間氣候變化專門委員會(IPCC) 第五次評估報告發(fā)布的 BCC-CSM1.1氣候模式數(shù)據(jù),采用最大熵(Maxent) 模型預(yù)測了4種氣候變化(RCP2.6,4.5,6.0和8.5) 情景下2050s及2070s 我國草地螟越冬區(qū)的變化情況,結(jié)果表明:我國草地螟越冬區(qū)將會較當(dāng)前在面積上有所增長并且位置上有向北移動的趨勢。上述研究結(jié)果表明:隨著全球氣候變暖,物種的現(xiàn)代潛在地理分布范圍將會擴大,本研究結(jié)論一致。

MaxEnt模型較其他生態(tài)位模型CLIMEX、GARP具有操作簡便、樣本需求量小和預(yù)測精度高等優(yōu)勢,但不可避免也存在部分局限性,如數(shù)據(jù)量較龐大時運行速度將會滯后等。其他環(huán)境因子變量如:未來的人類活動強度和植被覆蓋等數(shù)據(jù)對大沙鼠的潛在地理分布變化也有影響,由于未來時期的人類強度和植被覆蓋數(shù)據(jù)難以獲取,未加入對大沙鼠潛在地理分布的預(yù)測,在實際應(yīng)用時必須結(jié)合當(dāng)?shù)氐木C合條件。接下來,應(yīng)增加更為全面的環(huán)境變量,探討其對物種分布模型結(jié)果的影響。

致謝：感謝新疆產(chǎn)學(xué)研聯(lián)合培養(yǎng)研究生基地-新疆治蝗滅鼠指揮部辦公室提供平臺和數(shù)據(jù)支持。