西藏藏羚種群數(shù)量動態(tài)研究*

普布頓珠

(西藏自治區(qū)林業(yè)調(diào)查規(guī)劃研究院，西藏 拉薩 850000)

藏羚(Pantholopshodgsonii)屬偶蹄目(Artiodactyla)牛科(Bovidae)藏羚屬(Pantholops)動物，被《國家重點保護野生動物名錄》列為國家一級保護野生動物[1]，也是青藏高原特有的瀕危物種[2]。近年來，由于我國保護力度的加大，藏羚種群數(shù)量逐漸增加，瀕危程度連降兩級，在2016年被世界自然保護聯(lián)盟由“瀕危”降至“近危”[3]。藏羚在我國主要分布在青藏高原的羌塘區(qū)域，僅羌塘自然保護區(qū)種群數(shù)量就超過20×104只[4]。除西藏自治區(qū)外，藏羚在青海省和新疆自治區(qū)也有分布，西藏、青海、新疆3個自治區(qū)(省)藏羚棲息地分別占總棲息地面積的57.7%、2.7%和39.4%[5]。

物種的分布和數(shù)量的變化能直觀反映物種的現(xiàn)狀[6]，對其進行深入分析對該物種的保護與管理提供重要作用。20世紀末因非法偷獵活動猖獗，導(dǎo)致藏羚種群數(shù)量急劇下降，因此藏羚的種群數(shù)量及分布情況備受國內(nèi)外關(guān)注。早在1905年，Rawling在西藏一次就發(fā)現(xiàn)了近20 000只的藏羚群[8]。但到1990年，馮祚建[7]在喀喇昆侖山地區(qū)考察的僅記錄到3 940只藏羚。Fox等[9]于2000年和2002年間秋季在西藏阿魯盆地分別記錄到11 000和10 000只藏羚。目前，保守估計青藏高原地區(qū)的藏羚數(shù)量超過30×104只[10]。

雖然，目前關(guān)于藏羚種群數(shù)量及分布的研究較多[8-10]，但其研究不夠全面。因此，本研究通過文獻檢索、歷史資料收集，結(jié)合從1987—2006年開展的西藏藏羚的專項考察研究，對西藏藏羚種群數(shù)量及分布進行統(tǒng)計分析，從長期的種群數(shù)量及分布數(shù)據(jù)中獲取重要信息，以為未來對藏羚的保護提供科學(xué)依據(jù)。

1 材料和方法

1.1 研究區(qū)概況

西藏自治區(qū)位于中國的西南邊陲，青藏高原的西南部。面積122.84×104km2，約占中國總面積的八分之一，其面積僅次于新疆維吾爾自治區(qū)。南北最寬約1 000 km，東西最長達2 000 km，是世界上面積最大、海拔最高的高原，有“世界屋脊”之稱。羌塘高原是青藏高原的主體，北起昆侖山—可可西里山南麓的新疆、西藏分界線，南至那曲地區(qū)的扎加藏布河谷，東臨唐古拉山一線的西藏、青海分界線，西至念青唐古拉山脈，面積60×104km2，行政范圍隸屬于西藏那曲市和阿里地區(qū)兩地，是藏羚最重要聚集地[11]。

羌塘高原地勢開闊、起伏平緩、生物多樣性較低，植物蓋度也相對較低。但在羌塘高原卻誕生了眾多適應(yīng)高原高寒氣候的特有動物，除藏羚外，這里也是野牦牛(Bosmutus)、西藏野驢(Equuskiang)和雪豹(Pantherauncia)等國家一級保護動物的主要分布區(qū)。藏羚是草食性有蹄類的代表，在羌塘高原地區(qū)具有得天獨厚的生存優(yōu)勢[10]。

1.2 研究方法

根據(jù)西藏藏羚的專項考察研究，對其進行不同年份、不同地區(qū)和不同月份的藏羚種群數(shù)量及分布進行深度探索。通過權(quán)威媒體、會議報道及歷史資料對西藏藏羚種群分布數(shù)據(jù)進行收集。以“藏羚”“西藏”“種群分布”為關(guān)鍵詞搜集藏羚種群的數(shù)據(jù)。

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

基于專項考察研究項目數(shù)據(jù)，運用廣義線性混合模型，評估藏羚種群數(shù)量在空間和時間上的變化趨勢。其中種群數(shù)量為響應(yīng)變量，調(diào)查時間(調(diào)查方位)為解釋變量，路線長度和方位(時間)為隨機變量，利用卡方檢驗分析方位上的差異。所有數(shù)據(jù)整理及分析均在R(4.0.1)軟件中完成，數(shù)據(jù)可視化在“ggplot 2”包中完成。

2 結(jié)果與分析

2.1 藏羚種群的時空分布格局

根據(jù)西藏藏羚的專項考察研究項目調(diào)查數(shù)據(jù)分析，從圖1可以看出，1987年調(diào)查種群密度為0.57只/km2，2006年種群密度為1.73只/km2，年平均密度以0.06只/km2呈上升趨勢，表明總體數(shù)量呈增加趨勢。從空間格局上來看(圖2)，種群密度在0.10～1.72 只/km2范圍內(nèi)，種群密度差異較大，表明藏羚種群空間分布可能受地形地貌、氣候、植被、食物豐富度、人為活動等因素影響較大。

圖1 西藏藏羚種群密度變化——時間尺度

圖2 西藏藏羚種群密度——空間格局

2.2 藏羚種群動態(tài)變化

根據(jù)搜集的數(shù)據(jù)(文獻、新聞和書籍等)可以看出(圖3)，西藏藏羚在20世紀初種群數(shù)量最大，在20世紀80年代末種群數(shù)量急劇銳減到約5×104只，其數(shù)量降低約8倍；到21世紀，種群數(shù)量逐漸恢復(fù)，但也不到最初數(shù)量的一半，恢復(fù)速度以平均每年0.26×104只增長。從1989—2005年平均種群密度上來看，不同方位上密度差異不大，但其東北部種群密度高于其他方向(圖4)。

圖3 藏羚種群數(shù)量長期變化規(guī)律

圖4 1989—2005年間不同方位平均種群密度

2.3 藏羚分布

結(jié)合西藏藏羚相關(guān)數(shù)據(jù)記載，得到1905—2016年的種群數(shù)量及分布信息，不論是整體還是局部其種群數(shù)量差異都較大，其中共有10次完整數(shù)據(jù)，其它都只是局部記載(表1)，但也有一定的規(guī)律，即藏羚種群數(shù)量逐漸增加(密度估計)。根據(jù)記載可知[2]，藏羚主要分布在尼瑪縣、雙湖特區(qū)、申扎縣、班戈縣、改則縣和日土縣，占藏羚分布區(qū)的22.8%，密度大于>0.40只/km2。

表1 西藏藏羚分布情況

3 結(jié)論與討論

本研究通過信息收集方法對西藏藏羚種群數(shù)量動態(tài)進行分析，從局部來看，隨著時間的推移，其種群數(shù)量呈增加趨勢，空間格局也對種群密度有顯著影響；從整體上看，藏羚種群數(shù)量總體呈先降低后上升的趨勢，最大縮減了近90%(瀕臨滅絕)，而不同空間上的其種群密度也有所差異。

本研究雖然未實際監(jiān)測藏羚種群數(shù)量動態(tài)，但根據(jù)可靠數(shù)據(jù)來源，西藏藏羚種群整體呈現(xiàn)出明顯的動態(tài)變化規(guī)律(圖3)，其主要因素是管控程度。從時間尺度來看，1987—2006年局部地區(qū)藏羚種群密度隨時間推移逐漸增加，與整體種群數(shù)量呈現(xiàn)一致趨勢，其中可能由于調(diào)查樣線的長短、調(diào)查季節(jié)、遷徙路線等原因造成的種群密度偏差。20世紀初，藏羚數(shù)量約有100×104只，青藏高原上藏羚群隨處可見。在20世紀下半葉，由于市場對藏羚皮毛的大量需求和國家的不重視，藏羚數(shù)量每年平均以20 000 只驟減[15]，造成其棲息地嚴重片段化；到20世紀80年代末，種群數(shù)量縮減到最低(50 000 只左右)。1988年以后，伴隨著野生動物法律法規(guī)的出臺和國家級自然保護區(qū)的建立，展開有效保護并結(jié)合社區(qū)的宣傳，不斷提高了當(dāng)?shù)卮迕竦纳鷳B(tài)保護意識，藏羚種群數(shù)量得以恢復(fù)，即便如此，當(dāng)前藏羚數(shù)量也才恢復(fù)至20世紀初的一半。根據(jù)1989—2005年的密度值來看，西藏藏羚的種群密度增長值為6.25%，其扣除各種因素致死的藏羚(春夏交接時誤食有毒植物導(dǎo)致死亡；夏季雌性產(chǎn)羔后因中風(fēng)導(dǎo)致死亡；冬季雄性在擇偶交配后由于極度體弱、食物短缺和氣候突變就很容易導(dǎo)致死亡，或者因公藏羚擇偶爭斗死亡；極端氣候下長時間無法獲得食物受饑寒死亡；間隔幾年容易產(chǎn)生較大的疾病，尤其是肺熱病的死亡；天敵傷害；年老體弱而死亡)，估計全區(qū)的藏羚數(shù)量年增長率僅為6.2%，預(yù)計到2023年種群密度將會達到飽和[2]。雖然大的格局下，藏羚種群數(shù)量逐漸增加，但是由于其他因素(氣候、地形、食物、水資源和人為活動等)，也會給藏羚種群數(shù)量和分布造成一定的影響[15]，保護藏羚的棲息地及其遷徙路徑，建立和維護適合的生態(tài)廊道則會降低外在因素造成的影響[17]，對藏羚種群的保護十分有利。