王秋生,溫 璐,蘇旭坤

1 內(nèi)蒙古大學(xué)生態(tài)與環(huán)境學(xué)院,呼和浩特 010021

2中國(guó)科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100085

3 中國(guó)科學(xué)院大學(xué),北京 100049

根據(jù)《IPCC第六次評(píng)估報(bào)告》的評(píng)估結(jié)果顯示,在未來(lái)幾十年里,所有地區(qū)的氣候變化都將加劇且仍得不到有效遏制,而不可遏制的氣候變化仍被認(rèn)為是21世紀(jì)全球生物多樣性面臨的最主要威脅之一[1]。人為原因引起的氣候變化正以前所未有的速度改變著地球上物種的分布格局,最明顯的表現(xiàn)就是物種分布范圍向高緯度或高海拔移動(dòng)[2—3]。被譽(yù)為“世界屋脊”“地球第三極”的青藏高原[4],是珍稀野生動(dòng)物的天然棲息地和高原物種基因庫(kù),是中國(guó)乃至亞洲重要的生態(tài)安全屏障,是中國(guó)生態(tài)文明建設(shè)的重點(diǎn)地區(qū)之一。在青藏高原的氣候環(huán)境不斷變化的背景下,地表環(huán)境和生態(tài)系統(tǒng)變化對(duì)氣候變化的響應(yīng)與反饋過(guò)程將產(chǎn)生一系列影響[5—6]。地表環(huán)境和生態(tài)系統(tǒng)變化也直接影響著物種的分布格局。青藏高原的物種豐富度在氣候變化的影響下,從東到西逐漸下降；雖然東南山區(qū)擁有較高的物種豐富度,卻是生物多樣性損失最嚴(yán)重的地區(qū)。21世紀(jì)很有可能看到在青藏高原北部有蹄類物種增加的現(xiàn)象,青藏高原野生動(dòng)物的分布區(qū)域減少,物種豐富度格局也將發(fā)生顯著變化[7]。只根據(jù)當(dāng)前的分布情況建立自然保護(hù)區(qū)和國(guó)家公園來(lái)保護(hù)目標(biāo)物種是遠(yuǎn)遠(yuǎn)不夠的[8—9]。將自然氣候因素和人為氣候因素結(jié)合以獲得更準(zhǔn)確的分布預(yù)測(cè)結(jié)果,此研究使用的氣候數(shù)據(jù)采用CMIP6模式采用的共享社會(huì)經(jīng)濟(jì)路徑(Shared Socioeconomic Pathways,SSPs)包括SSP126、SSP245、SSP370和SSP585四種共享社會(huì)經(jīng)濟(jì)路徑情景,能夠反映輻射強(qiáng)迫和社會(huì)經(jīng)濟(jì)發(fā)展間的關(guān)聯(lián),體現(xiàn)氣候因素與社會(huì)經(jīng)濟(jì)因素間的綜合作用,更有利氣候變化預(yù)估與影響的研究[10]。因此,采用最新的氣候數(shù)據(jù),針對(duì)不同的氣候情景來(lái)評(píng)估物種分布,將保護(hù)行動(dòng)的規(guī)模和自然保護(hù)地的優(yōu)化和管理與氣候變化預(yù)測(cè)的結(jié)果相匹配便顯得極為重要[11—12]。

藏羚羊(Pantholopshodgsonii)被譽(yù)為“高原精靈”,其已有百萬(wàn)年的演化歷史,是青藏高原的特有種與旗艦種[13]。藏羚羊是《國(guó)家重點(diǎn)保護(hù)野生動(dòng)物名錄》中的國(guó)家一級(jí)保護(hù)野生動(dòng)物,《瀕危野生動(dòng)植物種國(guó)際貿(mào)易公約》中被列為瀕危野生動(dòng)物[14],2016年被世界自然保護(hù)聯(lián)盟(IUCN)瀕危野生動(dòng)物保護(hù)紅色目錄《Red List of Threatened Species》列為近危(NT,Near Threatened)物種(IUCN,2016)[15]。隨著保護(hù)政策實(shí)施和保護(hù)力度加強(qiáng),藏羚羊的種群數(shù)量逐漸恢復(fù)。粗略估計(jì),目前青藏高原藏羚羊數(shù)量大約為30萬(wàn)只[16],但近年來(lái),青藏高原的氣候變化、基礎(chǔ)設(shè)施建設(shè)以及放牧等人為干擾日益加劇,藏羚羊賴以生存的草場(chǎng)和棲息地被破壞,并導(dǎo)致生境破碎化和喪失,直接威脅著藏羚羊等有蹄類野生動(dòng)物的生存和種群延續(xù)[17]。對(duì)生活在生態(tài)脆弱地區(qū)的藏羚羊來(lái)說(shuō),活動(dòng)范圍或生活習(xí)慣都受到限制,極小的氣候事件都有可能對(duì)藏羚羊的棲息地產(chǎn)生不成比例的影響[18]。另外,我國(guó)青藏高原建立的自然保護(hù)區(qū)大多存在保護(hù)與經(jīng)濟(jì)發(fā)展的沖突,缺乏統(tǒng)籌規(guī)劃,藏羚羊等珍稀瀕危野生動(dòng)物存在保護(hù)空缺等問(wèn)題[4]。最大熵模型(MaxEnt)相比其他物種分布模型,結(jié)果更精準(zhǔn),在時(shí)空外推性方面也表現(xiàn)更好,被廣泛應(yīng)用于珍稀瀕危物種的地理分布及氣候變化背景下的生境時(shí)空格局演變的研究[19—23],如Luo等利用物種分布模型評(píng)估未來(lái)氣候變化對(duì)青藏高原32種有蹄類動(dòng)物的潛在影響[7]；晏婷婷利用物種分布模型評(píng)估未來(lái)氣候變化對(duì)我國(guó)特有物種大熊貓的潛在影響[24]。

因此,本研究中應(yīng)用物種分布模型,在當(dāng)前和未來(lái)氣候變化情景下模擬藏羚羊生境時(shí)空格局演變并確定影響其分布的最主要因素,明確百年尺度藏羚羊保護(hù)空缺變化。研究藏羚羊生境時(shí)空變化,對(duì)青藏高原生物多樣性保護(hù)與維持具有重要意義,對(duì)青藏高原自然保護(hù)地體系優(yōu)化提供科學(xué)支撐并助力青藏高原生態(tài)安全屏障建設(shè)。

1 材料與方法

1.1 研究區(qū)概況

青藏高原平均海拔4000米以上,地勢(shì)復(fù)雜,氣候獨(dú)特，常與南極和北極作比,形象地被稱為“世界第三極”[25]。青藏高原地理位置為73°19′—104°47′E,26°00′—39°47′N,總面積約為2.5×107km2。青藏高原特殊的地理環(huán)境和地表特征形成了極其復(fù)雜的氣候條件,擁有大量的藏羚羊適宜生境[13,16]。本文中的研究區(qū)邊界(圖1)依照中國(guó)科學(xué)院生態(tài)環(huán)境研究中心的中國(guó)生態(tài)功能區(qū)劃并結(jié)合保護(hù)區(qū)及青藏高原邊界。研究區(qū)內(nèi)已建立1個(gè)國(guó)家公園、3個(gè)國(guó)家級(jí)自然保護(hù)區(qū)和1個(gè)省級(jí)自然保護(hù)區(qū)用于保護(hù)藏羚羊等有蹄類野生動(dòng)物及其生境,分別是:三江源國(guó)家公園(面積1.5×105km2)、新疆阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)(面積4.68×104km2)、可可西里國(guó)家級(jí)自然保護(hù)區(qū)(面積5.03×104km2)、羌塘國(guó)家級(jí)自然保護(hù)區(qū)(面積3.02×105km2)和新疆中昆侖自治區(qū)級(jí)保護(hù)區(qū)(面積3.2×104km2)。

圖1 研究區(qū)位置示意圖

1.2 研究方法

1.2.1藏羚羊出現(xiàn)點(diǎn)野外調(diào)查

藏羚羊在不同季節(jié)的大致生境范圍主要是通過(guò)查閱文獻(xiàn)和報(bào)告獲取,并于2012年12月、2013年6月和2021年8—9月,開(kāi)展野外樣線調(diào)查,野外調(diào)查獲取阿爾金山國(guó)家級(jí)保護(hù)區(qū)和羌塘國(guó)家級(jí)保護(hù)區(qū)的藏羚羊出現(xiàn)點(diǎn)位信息。根據(jù)提前布設(shè)的樣線,沿樣線以10—20 km/h車速,通過(guò)肉眼觀察并輔以蔡斯8倍雙筒望遠(yuǎn)鏡,用手持GPS記錄觀測(cè)出現(xiàn)點(diǎn)坐標(biāo),并記錄距藏羚羊距離、雌雄亞幼、相對(duì)角度、地形和距水源地距離等指標(biāo),通過(guò)分析計(jì)算得到藏羚羊?qū)嶋H出現(xiàn)點(diǎn)位信息。通過(guò)樣線調(diào)查法共得到420個(gè)藏羚羊點(diǎn)位數(shù)據(jù),另外從相關(guān)文獻(xiàn)中獲取部分藏羚羊出現(xiàn)點(diǎn)位坐標(biāo)并進(jìn)行數(shù)字化,通過(guò)生境信息和ArcGIS的地理配準(zhǔn)對(duì)比校正,獲取了70個(gè)藏羚羊點(diǎn)位數(shù)據(jù)作為補(bǔ)充[26]。將所有樣點(diǎn)通過(guò)ArcGIS緩沖區(qū)分析法篩選,以0.5 km作為緩沖半徑,若1 km 范圍出現(xiàn)多點(diǎn)則僅保留1個(gè)。最終獲得青藏高原藏羚羊有效點(diǎn)位386個(gè)。

1.2.2氣象數(shù)據(jù)和地形數(shù)據(jù)獲取

本研究涉及9個(gè)環(huán)境因子,包括6個(gè)氣候因子、3個(gè)地形因子(海拔、坡度和坡向)。利用R語(yǔ)言的dismo安裝包,對(duì)Wordcilm2.1(http://www.worldclim.org/)網(wǎng)站下載的當(dāng)前氣候數(shù)據(jù)(2000—2018年)進(jìn)行處理。根據(jù)每月最高氣溫、最低氣溫、降水量三個(gè)數(shù)據(jù)項(xiàng),計(jì)算出2000—2018年的19個(gè)生物氣候變量(Bio1—Bio19)。利用SPSS對(duì)19個(gè)生物氣候因子進(jìn)行相關(guān)性分析,剔除相關(guān)性大于0.8的變量,最終確定溫度季節(jié)性變化(Bio4),溫度年變幅(Bio7),最干季均溫(Bio9),最干月降雨(Bio14),最濕季降雨(Bio16)和最干季降雨(Bio17)作為氣候變量,分辨率為30″(面積約0.83km2)[27—29]。海拔數(shù)據(jù)來(lái)源于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心,基于GIS技術(shù)生成坡度(slope)及坡向(aspect),分辨率統(tǒng)一為30″(面積約0.83km2)。

未來(lái)生物氣候變量(Bio1—Bio19)數(shù)據(jù)采用國(guó)際耦合模式計(jì)劃第六階段(CMIP6)的MIROC6模式模擬結(jié)果(https://esgf.node.llnl.gov/projects/cmip6/),模擬時(shí)段為2021—2100年。本文選取2021—2040年,2041—2061年,2061—2080年,2081—2100年四個(gè)時(shí)段。氣候壓迫情景選擇SSP126、SSP585兩種極端的路徑情景和SSP245中間路徑情景?？臻g分辨率與當(dāng)前氣候條件下的氣候因子數(shù)據(jù)保持一致。

1.2.3潛在生境識(shí)別

藏羚羊的潛在生境識(shí)別應(yīng)用最大熵模型(MaxEnt)完成,將已知的藏羚羊點(diǎn)位數(shù)據(jù)和利用ArcGIS 10.6軟件轉(zhuǎn)化工具生成的ASCII格式的氣候圖層導(dǎo)入MaxEnt軟件進(jìn)行運(yùn)算[29]。運(yùn)算過(guò)程中,設(shè)置訓(xùn)練數(shù)據(jù)集為75%,測(cè)試驗(yàn)證數(shù)據(jù)集為25%,迭代次數(shù)為10000,設(shè)置調(diào)控倍頻為1,重復(fù)訓(xùn)練數(shù)為10,進(jìn)行數(shù)據(jù)模擬。在MaxEnt軟件中,檢驗(yàn)各生物氣候變量對(duì)青藏高原藏羚羊地理分布的限制強(qiáng)度的生物氣候變量貢獻(xiàn)率是由Jackknife模塊計(jì)算得到[30]。

進(jìn)行氣候變化的影響分析時(shí),通過(guò)最大化卡帕值(maximum Kappa value)的閾值選取方法,將MaxEnt模型預(yù)測(cè)結(jié)果在ArcGIS 10.6軟件中通過(guò)重分類轉(zhuǎn)化成適宜生境(0)和適宜生境(1)的柵格數(shù)據(jù)[31—33]。通過(guò)ArcGIS柵格計(jì)算器統(tǒng)計(jì)基準(zhǔn)氣候和未來(lái)各種氣候條件下藏羚羊潛在適宜生境柵格數(shù)量得出對(duì)應(yīng)氣候條件下的藏羚羊潛在適宜生境面積。對(duì)基準(zhǔn)氣候條件下藏羚羊潛在適宜生境與未來(lái)不同氣候條件下的藏羚羊潛在適宜生境進(jìn)行疊加分析,公式為:

C=HN×10+HW,C∈(0,1,10,11)

式中,HN為基準(zhǔn)氣候條件下的預(yù)測(cè)結(jié)果；HW為未來(lái)某一氣候路徑情景。X=10表示中高適宜生境→低適宜生境；X=1表示低適宜生境→中高適宜生境；X=0表示不變低適宜生境；X=11不變中高適宜生境。

MaxEnt模型預(yù)測(cè)的結(jié)果是0到1之間的ASCII數(shù)據(jù),而且適生值越高,代表物種存在的概率就越高[34]。本研究中,我們選擇模型敏感度和特異性數(shù)值之和最大MSS (maximum sensitivity plus specificity)的閾值作為重分類的閾值,閾值為0.195和0.380。對(duì)模型結(jié)果中的13個(gè)MSS和ESS值進(jìn)行均值化處理,結(jié)合其他學(xué)者的閾值判定方法將預(yù)測(cè)結(jié)果劃分為三個(gè)等級(jí)[21,22,27]:適生值< 0.2為低適生區(qū)；0.2 ≤適生值< 0.4為中適生區(qū)；0.4 ≤適生值< 1為高適生區(qū)。中適生區(qū)和高適生區(qū)所在的區(qū)域是藏羚羊的適宜生境的概率更高。因此,本研究主要關(guān)注“中高適生區(qū)”的中適生區(qū)和高適生區(qū)。

2 結(jié)果與分析

2.1 模型精度驗(yàn)證

由圖2可知,預(yù)測(cè)分布結(jié)果的ROC曲線下面積(AUC)的訓(xùn)練集和測(cè)試集的均值為0.898,這表明MaxEnt模型預(yù)測(cè)得到的青藏高原藏羚羊的潛在分布區(qū)具有很好的精度,其預(yù)測(cè)結(jié)果可靠度高,不具有隨機(jī)性。

圖2 MaxEnt模型的受試者操作特性曲線

通過(guò)Jackknife模塊對(duì)各氣候因子的貢獻(xiàn)率進(jìn)行評(píng)估可知,溫度季節(jié)性變化(Bio4)訓(xùn)練得分值為最高,貢獻(xiàn)率為37.2%,其次是海拔,貢獻(xiàn)率為17.9%,隨后是最干月降雨(Bio14),貢獻(xiàn)率是15%。坡度,最干季降雨(Bio17),最濕季降雨(Bio16)貢獻(xiàn)率分別是10.7%,8.9%,8.3%。而溫度年變幅(Bio7),最干季均溫(Bio9),坡向貢獻(xiàn)率較小。這表明最干月降雨(Bio4)和海拔是影響藏羚羊分布的主導(dǎo)環(huán)境因子(圖3)。

圖3 刀切法檢測(cè)環(huán)境因子對(duì)藏羚羊在青藏高原分布的重要程度

2.2 T1時(shí)段藏羚羊生境與保護(hù)空缺識(shí)別

藏羚羊在青藏高原的中高適生面積為6.11×105km2,占青藏高原總面積的24.44%,其中適生區(qū)面積為2.70×105km2,占10.80%,高適生區(qū)面積為3.41×105km2,占13.64%。中高適生區(qū)主要是在30.38°N—37.30°N、74.38°E—82.31°E之間,主要分布在羌塘高原的東部,可可西里,阿爾金山南側(cè)且一直延伸到三江源國(guó)家公園附近,在昆侖山中部也有部分分布。中高適生區(qū)幾乎涵蓋整個(gè)可可西里國(guó)家級(jí)自然保護(hù)區(qū),也分布在羌塘國(guó)家級(jí)自然保護(hù)區(qū)的東南部分,新疆阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)和三江源國(guó)家公園的北部也有大面積的分布；中高適生區(qū)主要分布在西藏的北部,青海的西部和新疆的東南部。中適生區(qū)在新疆中昆侖自然保護(hù)區(qū)有部分點(diǎn)狀分布,祁連山西側(cè)也有少部分的中適生區(qū)分布(圖4)。青藏高原的國(guó)家級(jí)和省級(jí)保護(hù)區(qū)保護(hù)了青藏高原藏羚羊中高適生區(qū)的56.92%,而剩下43.08%的中高適生區(qū)約2.70×105km2在保護(hù)區(qū)外,未受到有效保護(hù)。未受到保護(hù)的中高適生區(qū)主要圍繞在羌塘國(guó)家級(jí)自然保護(hù)區(qū)、可可西里國(guó)家級(jí)自然保護(hù)區(qū)與三江源國(guó)家公園的周圍區(qū)域。另外,在青海湖周圍區(qū)域也存在部分保護(hù)空缺區(qū)。

圖4 當(dāng)前氣候情景下藏羚羊適宜生境分布

2.3 T2—T5時(shí)段藏羚羊生境變化

通過(guò)疊加分析可知,在SSP126,SSP245,SSP585三種壓迫情景下,中高適宜區(qū)向高緯度和高海拔轉(zhuǎn)移趨勢(shì)的明顯。棲息地的受影響程度較小。藏羚羊的核心分布區(qū)域(中高適生區(qū))仍然以西藏的中北部,青海的西部和新疆的東南部為主。從時(shí)間尺度上看,在三種不同壓迫情景下,中高適生區(qū)的變化趨勢(shì)基本一致,在整體上呈現(xiàn)出西南部分減少,東北部分增加的趨勢(shì),特別是在未來(lái)的T4和T5中高適生區(qū)的面積變化明顯。

與當(dāng)前氣候條件相比,T2階段的(SSP126,SSP245,SSP585)的情境下,三種情景下藏羚羊的中高適生區(qū)新增分別是6.72×104km2,6.75×104km2,6.53×104km2；喪失的中高適生區(qū)有6.20×104km2,6.19×104km2,5.69×104km2。喪失的中高適生區(qū)主要主要分布在羌塘保護(hù)區(qū)西南邊界處,新疆中昆侖自然保護(hù)區(qū)的西部和中部,和三江源國(guó)家公園西南緣與羌塘保護(hù)區(qū)交界處。新增的中高適生區(qū)主要在阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)的北部,三江源國(guó)家公園的中北部在保護(hù)區(qū)外圍也有大量零星分布。青海湖附近區(qū)域也有少量新增的藏羚羊中高適宜區(qū)。

與當(dāng)前氣候條件相比,T4階段的(SSP126,SSP245,SSP585)的情境中,三種情景下藏羚羊的中高適生區(qū)喪失的面積分別是7.68×104km2,6.16×104km2,7.99×104km2；新增的中高適生區(qū)有6.59×104km2,7.82×104km2,7.29×104km2。喪失的中高適生區(qū)主要是在喪失的中高適生區(qū)也呈現(xiàn)出點(diǎn)狀零星分布。隨著強(qiáng)迫程度的增加,中高適生區(qū)喪失的面積也會(huì)增加。T5階段的SSP585喪失中高適宜區(qū)面積最多,達(dá)到9.22×104km2,喪失的中高適生區(qū)呈現(xiàn)出部分面狀分布,如新疆中昆侖自然保護(hù)區(qū)中部及南部,羌塘高原的南部邊緣,新疆的南部的邊緣,西藏的北部,祁連山的西南部,青海湖西北部和三江源國(guó)家公園的中北部。羌塘自然保護(hù)區(qū)的西北部出現(xiàn)較多低適宜區(qū)轉(zhuǎn)變?yōu)橹懈哌m宜區(qū)的區(qū)域而可可西里國(guó)家級(jí)自然保護(hù)區(qū)中高適生區(qū)對(duì)于氣候變化的響應(yīng)程度略低。

圖5 不同時(shí)期不同氣候變化情景下藏羚羊適宜生境變化

2.4 T2—T5時(shí)段藏羚羊保護(hù)空缺識(shí)別

在三種氣候情景下,保護(hù)空缺面積整體減小幅度不明顯。青藏高原西南緣的保護(hù)空缺面積會(huì)隨氣候變化減少,青藏高原東北部分出現(xiàn)新的適宜棲息地。新疆中昆侖自治區(qū)級(jí)保護(hù)區(qū)在SSP585的壓迫情景下棲息地面積明顯減少(圖6)。羌塘國(guó)家級(jí)自然保護(hù)區(qū)的南緣部分保護(hù)空缺將向更高緯度擴(kuò)展。SSP3245情景下,保護(hù)空缺的面積受氣候變化的影響最小。

圖6 藏羚羊生境潛在分布與保護(hù)空缺分布

3 討論

本研究發(fā)現(xiàn),氣溫和海拔是對(duì)藏羚羊適生區(qū)分布的貢獻(xiàn)率最高的兩個(gè)生物氣候因子。其中,溫度季節(jié)性變化(Bio4)和海拔總貢獻(xiàn)率達(dá)到55.1%。過(guò)去的研究發(fā)現(xiàn)藏羚羊喜好生活在全年3/4以上的月份平均氣溫不足0℃,海拔介于3250—5500 m的高寒荒漠地區(qū),本次研究很好地驗(yàn)證了這一點(diǎn)。氣溫上升會(huì)影響野生動(dòng)物的生理行為,例如,顯著的溫度升高和氣溫的變化可能會(huì)引起熱應(yīng)激并影響適應(yīng)高山地區(qū)物種的生命率[35]。溫度升高可能會(huì)限制物種的分布,因?yàn)樗鼤?huì)限制覓食等重要活動(dòng)的可用時(shí)間,或者限制基因交流所需的時(shí)間[36—37]。藏羚羊的適宜生境的變化也同多數(shù)前人的研究一致,隨著氣候變暖,其分布范圍將逐漸縮小,并向高緯度和高海拔轉(zhuǎn)移[38]。本研究結(jié)果同樣表明青藏高原西南部部分地區(qū)的海拔和氣候條件已經(jīng)不滿足藏羚羊生存,藏羚羊可能會(huì)為了躲避低海拔及低緯度生境的濕熱化,往高海拔和高緯度地方遷移。藏羚羊損失的適宜生境的主要出現(xiàn)在青藏高原西南部,和中北部的低海拔地區(qū)；新出現(xiàn)的適宜區(qū)域整體上集中于青藏高原面的東部和東北部地區(qū),整體表現(xiàn)出向青藏高原東部和東北部遷移的趨勢(shì)。另外,藏羚羊適宜生境還表現(xiàn)出向核心地帶縮減的趨勢(shì)[23](圖7)。而在同一時(shí)期,藏羚羊適宜生境的面積并沒(méi)有隨強(qiáng)迫程度的增加,而呈現(xiàn)出一定的規(guī)律性。本研究針對(duì)此問(wèn)題做出了兩種假設(shè):(1)文中所用的基于CMIP6模式的未來(lái)氣候數(shù)據(jù)的準(zhǔn)確性和完整性有待于提高,CMIP6模式雖然擁有更精確的空間分辨率、改進(jìn)的云微物理過(guò)程參數(shù)化方案和更復(fù)雜的地球系統(tǒng)過(guò)程并在未來(lái)氣候模擬采用新的情景框架,但是CMIP6模式對(duì)青藏高原降水和氣溫的模擬能力仍然有限,對(duì)于青藏高原這樣的復(fù)雜地形區(qū)域,氣候模式的水平分辨率仍較為粗糙,冷偏差與濕偏差依然存在,并且偏差會(huì)隨著時(shí)間的推移而變大[39]；(2)目前關(guān)于氣候變化下高海拔的有蹄類動(dòng)物的潛在范圍變化研究還很少,這可能限制了我們對(duì)氣候變化有蹄類野生動(dòng)物潛在影響的認(rèn)識(shí)和對(duì)制定保護(hù)策略的指導(dǎo)。

圖7 未來(lái)80年藏羚羊適宜生境轉(zhuǎn)換

三江源國(guó)家公園、可可西里國(guó)家級(jí)自然保護(hù)區(qū)、羌塘國(guó)家級(jí)自然保護(hù)區(qū)、新疆阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)和新疆中昆侖自治區(qū)級(jí)自然保護(hù)區(qū)雖保護(hù)了56.92%藏羚羊的中高適生區(qū),但仍有43.18%的保護(hù)空缺區(qū)在保護(hù)區(qū)外,未受到有效保護(hù)。且在低適生區(qū)內(nèi)仍有大量零星的藏羚羊中高適宜生境。盡管藏羚羊在世界自然保護(hù)聯(lián)盟(IUCN)紅色名錄中的保護(hù)等級(jí)已經(jīng)下降,但在氣候變化的影響下,藏羚羊的滅絕率仍將有上升的風(fēng)險(xiǎn)[24]。對(duì)于藏羚羊適宜生境喪失嚴(yán)重性、藏羚羊的分布的復(fù)雜性以及當(dāng)前保護(hù)空缺的急迫性。建議相關(guān)部門在以國(guó)家公園為主體的自然保護(hù)地體系優(yōu)化的大背景下,綜合考慮藏羚羊及其他有蹄類野生動(dòng)物的生物學(xué)特性,調(diào)整保護(hù)邊界,并在保護(hù)地之間布設(shè)生態(tài)廊道,保障藏羚羊等瀕危有蹄類野生動(dòng)物棲息地完整性和連續(xù)性。

4 結(jié)論

本文利用386個(gè)藏羚羊?qū)嶋H分布點(diǎn)位數(shù)據(jù)和9個(gè)環(huán)境變量數(shù)據(jù),采用最大熵模型(MaxEnt)模擬青藏高原地區(qū)藏羚羊當(dāng)前分布格局以及氣候變化背景下青藏高原藏羚羊生境時(shí)空格局演變,并得到影響其分布的主導(dǎo)環(huán)境變量主要是溫度季節(jié)性變化(Bio4),海拔,最干月降雨(Bio14)。通過(guò)空間疊加分析,得到不同情景下的中高適生區(qū)和保護(hù)空缺分布。在三種氣候情境下藏羚羊的適宜生境都具有向高緯度、高海拔遷移趨勢(shì)。在五個(gè)自然保護(hù)地外仍然有近44%的中高適生區(qū)未受到保護(hù)。將不同氣候情景下的適宜生境有效聯(lián)系,綜合考慮藏羚羊遷徙廊道,采取相應(yīng)的保護(hù)措施,如擴(kuò)建和增建自然保護(hù)地、實(shí)施長(zhǎng)期監(jiān)測(cè)以獲得藏羚羊轉(zhuǎn)移和分布趨勢(shì)的有效信息、加強(qiáng)畜牧儲(chǔ)備平衡等。建議在今后的自然保護(hù)地管理中應(yīng)對(duì)藏羚羊保護(hù)空缺區(qū)域加強(qiáng)關(guān)注。