基于植被覆蓋度的藏羚羊棲息地時(shí)空變化研究

趙海迪，劉世梁,*，董世魁，蘇旭坤，張 翔

(1.北京師范大學(xué)環(huán)境學(xué)院水環(huán)境模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 100875; 2.阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)管理局，庫(kù)爾勒 841000)

基于植被覆蓋度的藏羚羊棲息地時(shí)空變化研究

趙海迪1，劉世梁1,*，董世魁1，蘇旭坤1，張 翔2

(1.北京師范大學(xué)環(huán)境學(xué)院水環(huán)境模擬國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 100875; 2.阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)管理局，庫(kù)爾勒 841000)

植被的分布及動(dòng)態(tài)變化在一定程度上反映了動(dòng)物棲息地的分布與變化?；谥脖贿b感影像數(shù)據(jù)及GIS空間分析，對(duì)阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)植被景觀的變化進(jìn)行分析，同時(shí)通過(guò)植被類(lèi)型要素識(shí)別出適宜藏羚羊生存的主要區(qū)域，對(duì)其植被覆蓋度的時(shí)空變化進(jìn)行分析。結(jié)果表明：保護(hù)區(qū)的植被覆蓋度處于較低水平，低覆蓋度植被區(qū)域占總植被覆蓋區(qū)的50%左右，且植被覆蓋度高的區(qū)域均集中在卡爾墩檢查站的東南部。植被覆蓋度在2000年、2005年、2010年間的增長(zhǎng)趨勢(shì)表明保護(hù)區(qū)植被對(duì)于動(dòng)物的承載能力不斷增大，植被的分布表明藏羚羊的活動(dòng)區(qū)域植被覆蓋度較高，而產(chǎn)羔區(qū)域的植被覆蓋度較低且海拔較高。對(duì)藏羚羊取食植被針茅草原植被覆蓋度的變化分析表明，針茅草原植被覆蓋度隨時(shí)間而有所提高，藏羚羊棲息地的面積可能有所增大，覆蓋度分布的變化表明藏羚羊棲息地有向保護(hù)區(qū)西南方向擴(kuò)散的趨勢(shì)。對(duì)植被覆蓋度與海拔關(guān)系的分析表明，高海拔區(qū)域植被覆蓋度較低，中度海拔區(qū)域植被覆蓋度較高，但是兩者并不存在顯著地相關(guān)性。

藏羚羊；阿爾金山自然保護(hù)區(qū)；歸一化植被指數(shù)(NDVI)；植被覆蓋度；針茅草原

被譽(yù)為‘高原精靈’的藏羚羊是一重要的野生動(dòng)物物種[1]，該物種主要分布于我國(guó)的青海、西藏、新疆三省3700—5500 m的高山荒漠草甸草原和高原草原等環(huán)境中[2- 3]。然而由于環(huán)境的破壞以及非法的狩獵活動(dòng)[4- 5]，使得藏羚羊的數(shù)量和分布在20世紀(jì)均出現(xiàn)明顯的降低，作為國(guó)家I級(jí)保護(hù)動(dòng)物已被列入《瀕危野生動(dòng)植物種國(guó)際貿(mào)易公約》(CITES)公約附錄I[6]。藏羚羊的保護(hù)引起了政府以及科研人員的廣泛關(guān)注[7- 8]，目前關(guān)于藏羚羊的研究主要集中于其行為習(xí)性、生理、季節(jié)性的遷徙及繁殖等領(lǐng)域[9- 12]，然而有關(guān)藏羚羊棲息地的研究卻相對(duì)較少。尤其在阿爾金山保護(hù)區(qū)，藏羚羊棲息地的分布情況至今仍未知。

植被條件決定著動(dòng)物的數(shù)量及品質(zhì)，而且也體現(xiàn)了氣候、土壤等其他環(huán)境要素的狀況，因此植被的分布可以在一定程度上決定棲息地的分布。近些年來(lái)，許多學(xué)者運(yùn)用植被的分布預(yù)測(cè)了動(dòng)物棲息地的分布狀況。Vina等通過(guò)MODIS數(shù)據(jù)指示植被的分布從而成功的預(yù)測(cè)了臥龍保護(hù)區(qū)大熊貓的分布情況[13]。Ito等通過(guò)利用NDVI預(yù)測(cè)蒙古瞪羚的分布，指出NDVI可以作為動(dòng)物棲息地的有效指示指標(biāo)[14]。張軍等通過(guò)分析多個(gè)植被性狀與高原鼢鼠棲息地適宜度的相關(guān)關(guān)系得出植被的覆蓋度可以用于棲息地的預(yù)測(cè)[15]。而且有的學(xué)者指出植被的動(dòng)態(tài)變化可以有效的指示動(dòng)物棲息地的變化情況[16]。遙感歸一化植被指數(shù)(NDVI)與生物量、葉面積指數(shù)有較好的相關(guān)關(guān)系，是反映區(qū)域植被動(dòng)態(tài)變化的最敏感指標(biāo)，能夠相當(dāng)精確地反映植被的季節(jié)和年際變化，常用于大尺度的植被動(dòng)態(tài)監(jiān)測(cè)[17- 19]。因此，本文綜合利用遙感影像NDVI數(shù)據(jù)及GIS空間分析，對(duì)阿爾金山自然保護(hù)區(qū)內(nèi)總體植被及符合藏羚羊食性的植被覆蓋度進(jìn)行年際間的動(dòng)態(tài)變化分析，進(jìn)而初步預(yù)測(cè)棲息地及該區(qū)承載力的變化情況。該研究能夠?yàn)榘柦鹕奖Ｗo(hù)區(qū)藏羚羊的保護(hù)提供理論基礎(chǔ)和科學(xué)依據(jù)，以期更好地保護(hù)瀕危珍稀的藏羚羊。

1 研究區(qū)概況

阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)(36°00′—37°49′N(xiāo)，87°10′—91°18′E)位于我國(guó)新疆維吾爾自治區(qū)、青海省與西藏自治區(qū)三省的交界處，地處青藏高原的西北邊緣地帶，連接可可西里和羌塘自然保護(hù)區(qū)，位置示意圖見(jiàn)圖1。該保護(hù)區(qū)是世界上內(nèi)陸面積最大的保護(hù)區(qū)之一，總面積45000 km2，其南北距離為190 km，東西為360 km，海拔在3748 m至6948 m之間變化。保護(hù)區(qū)于1983年5月正式建立，1985年3月被列為國(guó)家級(jí)自然保護(hù)區(qū)。保護(hù)區(qū)內(nèi)全年無(wú)絕對(duì)無(wú)霜期，無(wú)明顯四季之分，年平均氣溫在0 ℃以下，降水量在300 mm左右，高原多風(fēng)特征明顯，日照輻射強(qiáng)、地溫變化大。由于地處偏遠(yuǎn)、自然條件惡劣的亞歐大陸高寒腹地，加之道路艱險(xiǎn)、高原缺氧，使這里人跡罕至[20]，是我國(guó)四大無(wú)人區(qū)域之一，同時(shí)也使得保護(hù)區(qū)保留了豐富的自然資源及珍稀的野生動(dòng)植物，是我國(guó)少有的高原生態(tài)系統(tǒng)保存完好的地區(qū)，為中外學(xué)者所矚目[5]。

圖1 阿爾金山自然保護(hù)區(qū)位置Fig.1 The location of Altun Nature Reserve

2 研究方法

2.1 數(shù)據(jù)來(lái)源

本研究采用的植被景觀數(shù)據(jù)來(lái)源于研究區(qū)2000年、2005年、2010年3期的LANDSAT TM影像，該影像來(lái)源于美國(guó)聯(lián)邦地質(zhì)調(diào)查局(USGS)和中國(guó)科學(xué)院對(duì)地觀測(cè)與數(shù)字地球科學(xué)中心，影像的條帶號(hào)/行編號(hào)(Path/Row)為139/34、139/35、140/34、140/35、141/34、141/35，成像時(shí)間均在植被覆蓋度較高的7—9月份，并利用ERDAS IMAGE圖像處理軟件結(jié)合野外調(diào)查對(duì)3期影像進(jìn)行人工目視判讀和監(jiān)督分類(lèi)，結(jié)合實(shí)地調(diào)研驗(yàn)證獲得3個(gè)時(shí)期的保護(hù)區(qū)土地利用類(lèi)型圖，其影像的分類(lèi)精度達(dá)到80%。然后提取出各期植被景觀類(lèi)型并進(jìn)行合并，從而對(duì)3期中有植被覆蓋的區(qū)域進(jìn)行植被覆蓋度的變化分析，以預(yù)測(cè)保護(hù)區(qū)整體動(dòng)物棲息地的變化。為了以植被的動(dòng)態(tài)變化更好地反應(yīng)藏羚羊棲息地的變化，該研究提取適合藏羚羊捕食的植被類(lèi)型區(qū)域，進(jìn)一步基于NDVI進(jìn)行藏羚羊適宜植被區(qū)域的植被覆蓋度動(dòng)態(tài)變化分析。植被類(lèi)型數(shù)據(jù)由阿爾金山保護(hù)區(qū)管理中心提供。研究利用2000、2005、2010年的SPOT-VEGETATION逐旬NDVI數(shù)據(jù)，該數(shù)據(jù)由國(guó)家自然科學(xué)基金委員會(huì)“中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心”網(wǎng)站下載，其空間分辨率為1 km×1 km。該數(shù)據(jù)集由比利時(shí)佛萊芒技術(shù)研究所VEGETATION影像處理中心經(jīng)輻射定標(biāo)、大氣校正、云檢測(cè)和幾何校正等預(yù)處理，最后生成10 d最大化合成的NDVI數(shù)據(jù)，再通過(guò)DN=([NDVI]+0.1)/0.004轉(zhuǎn)換為0—250的DN值[21]。

2.2 植被覆蓋度動(dòng)態(tài)變化分析

首先提取3期土地利用圖中的植被景觀部分并對(duì)3期植被進(jìn)行植被覆蓋度的計(jì)算與比較，分析植被覆蓋度的變化情況。然后提取適宜藏羚羊捕食的植被類(lèi)型區(qū)域進(jìn)行主要棲息地的植被動(dòng)態(tài)分析。藏羚羊適宜的海拔為3700—5500 m[22]，其捕食植被類(lèi)型有禾本科、豆科、菊科、薔薇科、玄參科、莎草科[23- 24]，而根據(jù)保護(hù)區(qū)所提供的植被類(lèi)型圖及相關(guān)文獻(xiàn)[20]可知阿爾金山保護(hù)區(qū)的優(yōu)勢(shì)植被類(lèi)型為針茅草原，且根據(jù)實(shí)地野外調(diào)查發(fā)現(xiàn)藏羚羊大多分布于針茅草原區(qū)域內(nèi)。因此，本研究選取了針茅草原進(jìn)行植被覆蓋度的變化分析，進(jìn)而對(duì)藏羚羊棲息地的變化進(jìn)行預(yù)測(cè)，針茅草原在保護(hù)區(qū)的分布見(jiàn)圖2。

圖2 研究區(qū)針茅草原分布Fig.2 Distribution of feather grass steppe in the study area

要使植被指數(shù)能夠定量的反映植被信息，需要賦予NDVI值以相應(yīng)的植被覆蓋含義，將植被指數(shù)轉(zhuǎn)換為植被覆蓋等級(jí)使得對(duì)于植被變化的定量評(píng)價(jià)更加直觀[25]。利用NDVI提取植被覆蓋度的方法主要有經(jīng)驗(yàn)?zāi)Ｐ头?、植被指?shù)法和混合像元分解法[26]，而目前應(yīng)用比較廣泛的是Gutman所提出的混合像元分解法[27]。其基本原理為假定每個(gè)像元所對(duì)應(yīng)的地表單元上只有植被和裸地兩種覆蓋類(lèi)型，每個(gè)像元的NDVI為植被和裸地所對(duì)應(yīng)的NDVI值的加權(quán)和，其權(quán)值即為兩種覆蓋類(lèi)型在像元中的面積百分比[28]。由此可知植被覆蓋度的計(jì)算公式為[29]：

(1)

式中,fc為植被覆蓋度，NDVIi為像元的歸一化植被指數(shù)，NDVIsoil為純裸地覆蓋像元的最小值，NDVIveg為純植被覆蓋像元的最大值。本研究中NDVIsoil值根據(jù)實(shí)際值及相關(guān)文獻(xiàn)[27]，各年均取統(tǒng)一值0.01，NDVIveg為每年所有植被柵格中的最大值。由于年最大NDVI可以較好的反映該年度植被長(zhǎng)勢(shì)最好季節(jié)的植被覆蓋狀況，本研究采用每年各柵格最大化合成的值進(jìn)行植被覆蓋度的分析[30]。

3 結(jié)果與分析

3.1 植被覆蓋度動(dòng)態(tài)變化

利用2000年、2005年、2010年植被覆蓋區(qū)域每個(gè)柵格中年最大合成的NDVI值計(jì)算相應(yīng)年份的植被覆蓋度，然后根據(jù)水利部2008年頒布的《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)》中植被覆蓋度分級(jí)標(biāo)準(zhǔn)，將植被覆蓋度劃分為5個(gè)等級(jí)：<30%(低覆蓋度)、30%—45%(中低覆蓋度)、45%—60%(中等覆蓋度)、60%—75%(中高覆蓋度)和>75%(高覆蓋度)[29]。保護(hù)區(qū)內(nèi)大部分的植被覆蓋區(qū)域具有較低的覆蓋度，低覆蓋度植被占總植被覆蓋區(qū)域的一半左右，而中高覆蓋度區(qū)域所占比例較少，且主要分布于卡爾墩檢查站的西南部及依協(xié)克帕提管護(hù)站附近(圖3)。

圖3 不同植被覆蓋度的空間分布Fig.3 The spatial distribution of different vegetation coverage

植被覆蓋度整體上由低覆蓋度向高覆蓋度轉(zhuǎn)移，呈增加的趨勢(shì)(圖3)，其中，低覆蓋度植被逐年減少，而其余4類(lèi)覆蓋度植被均出現(xiàn)增加的趨勢(shì)(表1)。這表明保護(hù)區(qū)內(nèi)植被面積隨著時(shí)間而增加，即保護(hù)區(qū)對(duì)于動(dòng)物的植被承載能力在增大。這也在另一方面反映出動(dòng)物生存的棲息地有可能在逐年增多。根據(jù)野外調(diào)查及相關(guān)文獻(xiàn)[31]可知，位于保護(hù)區(qū)東北部的阿牙克庫(kù)木湖流域是藏羚羊的主要棲息地，而大部分的中、高植被覆蓋像元均在該流域內(nèi)(圖3)，這表明藏羚羊?qū)τ跅⒌氐倪x擇與植被的狀況具有很大的相關(guān)性。另外藏羚羊在保護(hù)區(qū)內(nèi)的產(chǎn)羔地主要位于保護(hù)區(qū)的西南部[32- 33]，而該區(qū)域內(nèi)的植被覆蓋度較低，這表明藏羚羊選擇的產(chǎn)羔地植被比較稀疏，其對(duì)產(chǎn)羔地的植被要求較低。

表1 各植被覆蓋度所占比例

3.2針茅草原區(qū)植被覆蓋度變化分析

對(duì)比針茅草原分布及不同覆蓋度植被的分布(圖2，圖3)，可以發(fā)現(xiàn)與保護(hù)區(qū)整體植被覆蓋度相似，針茅草原多處于低覆蓋度區(qū)域。為對(duì)比針茅草原在年際間的變化，對(duì)2000年、2005年、2010年針茅草原植被類(lèi)型區(qū)進(jìn)行植被覆蓋度的變化分析。針茅草原中植被覆蓋度相對(duì)較高的區(qū)域集中于卡爾墩檢查站的西南部區(qū)域(圖4)，同時(shí)該區(qū)域也是整個(gè)保護(hù)區(qū)內(nèi)植被覆蓋度較高的區(qū)域，且處于藏羚羊主要的活動(dòng)區(qū)阿牙克庫(kù)木湖流域。這表明藏羚羊?qū)τ跅⒌氐倪x擇與植被的覆蓋情況具有很大相關(guān)性。

圖4 針茅草原植被覆蓋度Fig.4 Vegetation coverage of feather grass steppe

對(duì)比3期針茅草原的植被覆蓋度，可以發(fā)現(xiàn)針茅草原植被覆蓋度整體上隨著時(shí)間呈增加的趨勢(shì)(圖4)，3期平均植被覆蓋度分別為23.32%、24.38%、24.68%，植被覆蓋度非常低的區(qū)域逐漸減少，由針茅草原內(nèi)低覆蓋度向中覆蓋度以及相對(duì)較高的覆蓋度轉(zhuǎn)變。適宜藏羚羊取食的針茅草原覆蓋度的增大表明阿爾金山保護(hù)區(qū)適宜藏羚羊取食的植被面積加大，保護(hù)區(qū)對(duì)于藏羚羊的承載能力增強(qiáng)。同時(shí)針茅草原覆蓋度的增加將使得藏羚羊棲息地的分布范圍擴(kuò)大，藏羚羊的分布將由保護(hù)區(qū)東部向西南部覆蓋度提高區(qū)域擴(kuò)散。該結(jié)果與中國(guó)環(huán)境報(bào)等所報(bào)道的阿爾金山自然保護(hù)區(qū)的科考結(jié)果相吻合。

3.3 植被覆蓋度與海拔關(guān)系分析

為進(jìn)一步通過(guò)植被覆蓋度隨海拔的分布情況預(yù)測(cè)分析動(dòng)物棲息地的海拔分布，對(duì)植被覆蓋度與海拔的關(guān)系進(jìn)行了分析。對(duì)比植被覆蓋分布與海拔分布的情況(圖3，圖5)，可以看出低覆蓋度植被大多處于高海拔或低海拔區(qū)域，而中、高植被覆蓋度區(qū)域的海拔高度大多處于中間水平，這表明動(dòng)物的棲息地在中度海拔區(qū)域分布的可能性較大。由不同植被覆蓋度區(qū)域的平均海拔可以看出(表2)，植被覆蓋度與平均海拔沒(méi)有顯著地相關(guān)性，從總體上看低覆蓋度植被及中、高覆蓋度植被的平均海拔較低，中、高覆蓋度植被的海拔分布比較集中。藏羚羊的主要棲息地阿牙克庫(kù)木湖流域，平均海拔相對(duì)較低，而在藏羚羊的產(chǎn)羔地西南部地區(qū)海拔較高(圖5)。

圖5 植被覆蓋區(qū)域的海拔分布Fig.5 The altitude distribution of areas with vegetation covered

Table2Averageelevationofdifferentvegetationcoveragearea

覆蓋度Coverage年份Year200020052010低Low451844294466中低Mediumlow463746464632中等Medium456846174561中高M(jìn)ediumhigh449544794432高High460145344416

4 結(jié)論與討論

植被的變化在一定程度上可以反映動(dòng)物棲息地的變化，對(duì)棲息地內(nèi)的植被動(dòng)態(tài)進(jìn)行分析，可以預(yù)測(cè)棲息地在空間及面積上的變化。本研究主要從植被覆蓋度的角度分析了阿爾金山保護(hù)區(qū)植被及適宜藏羚羊生存的針茅草原覆蓋度的空間分布及動(dòng)態(tài)變化情況。基于2000、2005、2010年植被覆蓋區(qū)域，對(duì)研究區(qū)植被覆蓋度的分析表明整個(gè)保護(hù)區(qū)的植被覆蓋度處于較低水平，且覆蓋較高的區(qū)域集中分布在卡爾墩檢查站及依協(xié)克帕提管護(hù)站附近。由實(shí)地的野外調(diào)查及相關(guān)資料表明藏羚羊在阿爾金山的阿牙克庫(kù)木湖流域分布較多，而在本文的研究中植被覆蓋度相對(duì)高的區(qū)域均處于該流域內(nèi)，再次驗(yàn)證了通過(guò)植被覆蓋度指示棲息地分布的可行性。研究還發(fā)現(xiàn)藏羚羊主要活動(dòng)區(qū)域的植被覆蓋度較高，而產(chǎn)羔地的植被覆蓋度較低，且海拔較高，表明藏羚羊在產(chǎn)羔時(shí)對(duì)于植被的要求較低而更注重的是安全性。對(duì)植被覆蓋度動(dòng)態(tài)變化的研究表明，保護(hù)區(qū)植被覆蓋度隨時(shí)間呈增大的趨勢(shì)，從一定程度上反映了保護(hù)區(qū)對(duì)動(dòng)物承載能力的增強(qiáng)。針茅草原植被覆蓋度也隨著時(shí)間而增大，表明保護(hù)區(qū)對(duì)于藏羚羊的承載能力有所提高，且由植被覆蓋度的分布變化可以預(yù)測(cè)藏羚羊棲息地可能會(huì)由東部向西南部方向分散。同時(shí)，對(duì)海拔高度與植被覆蓋度的關(guān)系進(jìn)行了分析，結(jié)果表明兩者沒(méi)有顯著地相關(guān)性，中、高植被覆蓋度區(qū)域大多分布于中度海拔區(qū)，高海拔區(qū)域植被覆蓋度較低，從而可以推測(cè)出動(dòng)物偏好的海拔高度應(yīng)為中度海拔。

阿爾金山自然保護(hù)區(qū)藏羚羊的棲息地分布狀況仍然不為人知，通過(guò)藏羚羊適宜生存的環(huán)境要素識(shí)別出其適宜的植被景觀區(qū)域，對(duì)其進(jìn)行動(dòng)態(tài)變化的分析具有十分重要的意義。根據(jù)本研究的結(jié)果，建議保護(hù)區(qū)加強(qiáng)對(duì)卡爾墩檢查站及依協(xié)克帕提管護(hù)站附近的集中保護(hù)工作，同時(shí)由于藏羚羊棲息地可能向西南部擴(kuò)散，因此還需適當(dāng)?shù)財(cái)U(kuò)大保護(hù)的范圍。對(duì)于藏羚羊棲息地的識(shí)別，還需要綜合考慮人類(lèi)活動(dòng)干擾以及坡度、離水源距離等環(huán)境因素，從而與植被覆蓋度的分布相結(jié)合更加準(zhǔn)確的識(shí)別出藏羚羊的棲息地，為基于植被覆蓋度的棲息地動(dòng)態(tài)變化分析提供基礎(chǔ)。另外，實(shí)地的野外監(jiān)測(cè)也需要進(jìn)一步實(shí)施，從而更好地驗(yàn)證植被遙感數(shù)據(jù)以及研究結(jié)果。

[1] Feng Z, Fan B, Li K, Zhang Q D, Yang Q S, Liu B.Molecular characteristics of Tibetan antelope (Pantholopshodgsonii) mitochondrial DNA control region and phylogenetic inferences with related species.Small Ruminant Research, 2008, 75(2/3): 236- 242.

[2] Pan C, Wang P.Dancing elf on the plateau-Tibetan antelope.Forestry of China, 2007, (7A): 10- 13.

[3] Wu X M, Zhang H F.Resources regarding populations of Tibetan Antelope (Pantholopshodgsoni) and the status of its protection.Chinese Journal of Nature, 2011, 33(3): 143- 154.

[4] Li Z C, Yang Q S, Zhang H B, Liu Z H, Yuan G.Population number and distribution of big mammals in eastern Arjin Mountain Nature Reserve, Xinjiang.Sichuan Journal of Zoology, 2006, 25(1): 92- 95.

[5] Ma C.Analysis of human-induced infection in Aerjin Mountain National Nature Reserve.Arid Environmental Monitoring, 2004, 18(2): 101- 102.

[6] Zhou H, Li D Q, Zhang Y G, Yang T, Liu Y.Genetic diversity of microsatellite DNA loci of Tibetan antelope (Chiru,Pantholopshodgsonii) in Hoh Xil National Nature Reserve, Qinghai, China.Journal of Genetics and Genomics, 2007, 34(7): 600- 607.

[7] Liang Q H.The present situation investigation of Tibetan antelope in situ protection.Environmental Protection, 2008, 406(20): 51- 53.

[8] Liu S, Liu F H, Hou G J, A B D.Chasing the footprint of the life-sidelights of two action to protect the Tibetan antelope on Altun mountain.Forestry of Xinjiang, 2005, (6): 18- 19.

[9] Schaller G B, Kang A L, Cai X B, Liu Y L.Migratory and calving behavior of Tibetan antelope population.Acta Theriologica Sinica, 2006, 26(2): 105- 113.

[10] Zhang F, Jiang Z.Mitochondrial phylogeography and genetic diversity of Tibetan gazelle (Procaprapicticaudata): implications for conservation.Molecular Phylogenetics and Evolution, 2006, 41(2): 313- 321.

[11] Lian X M, Zhang T Z, Cao Y F, Su J P, Thirgood S.Group size effects on foraging and vigilance in migratory Tibetan antelope.Behavioural Processes, 2007, 76(3): 192- 197.

[12] Schaller G B, Ren J R, Qiu M J.Observations on the Tibetan antelope (Pantholopshodgsoni).Applied Animal Behaviour Science, 1991, 29(1/4): 361- 378.

[14] Ito T Y, Miura N, Lhagvasuren B, Enkhbileg D, Takatsuki S, Tsunekawa A, Jiang Z.Satellite tracking of Mongolian gazelles (Procapragutturosa) and habitat shifts in their seasonal ranges.Journal of Zoology, 2006, 269(3): 291- 298.

[15] Zhang J, Ge Q Z, Zhang W G, Jiang X L.Correlations between vegetation properties and fitness of plateau zoker habitat.Pratacultural Science, 2011, 28(5): 836- 840.

[16] Stow D, Hamada Y, Coulter L, Anguelova Z.Monitoring shrubland habitat changes through object-based change identification with airborne multispectral imagery.Remote Sensing of Environment, 2008, 112(3): 1051- 1061.

[17] Song F Q, Xing K X, Liu Y, Liu Z C, Kang M Y.Monitoring and assessment of vegetation variation in Northern Shaanxi based on MODIS/NDVI.Acta Ecologica Sinica, 2011, 31(2): 354- 363.

[18] Li H D, Shen W S, Cai B F, Ji D, Zhang X Y.The coupling relationship between variations of NDVI and change of aeolian sandy land in the Yarlung Zangbo River Basin of Tibet, China.Acta Ecologica Sinica, 2013, 33(24): 7729- 7738.

[19] Yuan L H, Jiang W G, Shen W M, Liu Y H, Wang W J, Tao L L, Zheng H, Liu X F.The spatio-temporal variations of vegetation cover in the Yellow River Basin from 2000 to 2010.Acta Ecologica Sinica, 2013, 33(24): 7798- 7806.

[20] Cui D F, Cui N R, Hai Y, Li X S.Plant species diversity in the Arjin Mountain National Nature Reserve areas.Journal of Shihezi University: Natural Science, 1998, 2(4): 281- 288.

[21] Song Y, Ma M G.Study on vegetation cover change in Northwest China based on SPOT VEGETATION data.Journal of Desert Research, 2007, 27(1): 89- 93.

[22] Schallar G B.Wildlife of the Tibetan Steppe.Chicago, USA: University of Chicago Press, 1998: 40- 79.

[23] Cao Y F, Su J P, Lian X M, Zhang T Z, Cui Q H.Food habits of Tibetan antelope (Pantholopshodgsoni) in the Kekexili Nature Reserve.Acta Theriologica Sinica, 2008, 28(1): 14- 19.

[24] Li Z Q, Jiang Z G.Diet analysis for Tibetan gazelle (Procaprapicticaudata) in Tianjun area, Qinghai province, China.Acta Theriologica Sinica, 2007, 27(1): 64- 67.

[25] Zhou Z Y, Chu S L, Wang Z W, Chen Q G.Analysis of vegetation coverage change based on NDVI-A case study in Ganzhou area, Zhangye city, Gansu.Pratacultural Science, 2008, 25(12): 23- 29.

[26] Li D K, Fan J Z, Wang J.Change characteristics and their causes of fractional vegetation coverage (FVC) in Shaanxi Province.Chinese Journal of Applied Ecology, 2010, 21(11): 2896- 2903.

[27] Niu B R, Liu J R, Wang Z W.Remote sensing information extraction based vegetation fraction in semiarid area.Geo-information Science, 2005, 7(1): 84- 97.

[28] Wu Y, Zeng Y, Zhao Y, Wu B F, Wu W B.Monitoring and dynamic analysis of fractional vegetation cover in the Hai River basin based on MODIS data.Resources Science, 2010, 32(7): 1417- 1424.

[29] Ding L Z, Xu G F, Lu J B, Zhang D S, Fang B F.Landscape fragmentation and its effect on biodiversity.Journal of Jiangsu Forestry Science & Technology, 2005, 32(4): 45- 57.

[30] Yang F, He D H.Effects of habitat fragmentation on biodiversity.Ecologic Science, 2006, 25(6): 564- 567.

[31] Gu J H.The ungulates in east Kunlun and Altun mountains region of Xinjiang.Arid Zone Research, 1987, (3): 56- 68.

[32] Buzzard P J, Wong H M, Zhang H B.Population increase at a calving ground of the Endangered Tibetan antelopePantholopshodgsoniiin Xinjiang, China.Oryx, 2012, 46(2): 266- 268.

[33] Bleisch W V, Buzzard P J, Zhang H B, Xü D H, Liu Z H, Li W D, Wong H.Surveys at a Tibetan antelopePantholopshodgsoniicalving ground adjacent to the Arjinshan Nature Reserve, Xinjiang, China: decline and recovery of a population.Oryx, 2009, 43(2): 191- 196.

參考文獻(xiàn):

[2] 潘成，汪平.高原上飛舞的精靈——藏羚羊.中國(guó)林業(yè), 2007, (7A): 10- 13.

[3] 吳曉民，張洪峰.藏羚羊種群資源及其保護(hù).自然雜志, 2011, 33(3): 143- 154.

[4] 李增超，楊奇森，張會(huì)斌，劉志虎，袁剛.新疆阿爾金山自然保護(hù)區(qū)東部大型獸類(lèi)數(shù)量和分布.四川動(dòng)物, 2006, 25(1): 92- 95.

[5] 馬超.阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)人為活動(dòng)的影響簡(jiǎn)析.干旱環(huán)境監(jiān)測(cè), 2004, 18(2): 101- 102.

[7] 梁清華.藏羚羊就地保護(hù)的現(xiàn)狀調(diào)查.環(huán)境保護(hù), 2008, 406(20): 51- 53.

[8] 劉晟, 劉鳳和, 侯國(guó)俊, 阿不都.追逐生命的足跡——阿爾金山保護(hù)藏羚羊二號(hào)行動(dòng)側(cè)記.新疆林業(yè), 2005, (6): 18- 19.

[15] 張軍, 葛慶征, 張衛(wèi)國(guó), 江小雷.植被性狀與高原鼢鼠棲息地適合度的關(guān)系.草業(yè)科學(xué), 2011, 28(5): 836- 840.

[17] 宋富強(qiáng), 邢開(kāi)雄, 劉陽(yáng), 劉志超, 康慕誼.基于MODIS/NDVI的陜北地區(qū)植被動(dòng)態(tài)監(jiān)測(cè)與評(píng)價(jià).生態(tài)學(xué)報(bào), 2011, 31(2): 354- 363.

[18] 李海東, 沈渭?jí)? 蔡博峰, 紀(jì)迪, 張曉勇.雅魯藏布江流域NDVI變化與風(fēng)沙化土地演變的耦合關(guān)系.生態(tài)學(xué)報(bào), 2013, 33(24): 7729- 7738.

[19] 袁麗華, 蔣衛(wèi)國(guó), 申文明, 劉穎慧, 王文杰, 陶亮亮, 鄭華, 劉孝富.2000—2010年黃河流域植被覆蓋的時(shí)空變化.生態(tài)學(xué)報(bào), 2013, 33(24): 7798- 7806.

[20] 崔大方, 崔乃然, 海鷹, 李學(xué)禹.阿爾金山自然保護(hù)區(qū)植物物種多樣性分析.石河子大學(xué)學(xué)報(bào): 自然科學(xué)版, 1998, 2(4): 281- 288.

[21] 宋怡, 馬明國(guó).基于SPOT VEGETATION數(shù)據(jù)的中國(guó)西北植被覆蓋變化分析.中國(guó)沙漠, 2007, 27(1): 89- 93.

[23] 曹伊凡, 蘇建平, 連新明, 張同作, 崔慶虎.可可西里自然保護(hù)區(qū)藏羚羊的食性分析.獸類(lèi)學(xué)報(bào), 2008, 28(1): 14- 19.

[24] 李忠秋, 蔣志剛.青海省天峻地區(qū)藏原羚的食性分析.獸類(lèi)學(xué)報(bào), 2007, 27(1): 64- 67.

[25] 周兆葉, 儲(chǔ)少林, 王志偉, 陳全功.基于NDVI的植被覆蓋度的變化分析——以甘肅省張掖市甘州區(qū)為例.草業(yè)科學(xué), 2008, 25(12): 23- 29.

[26] 李登科, 范建忠, 王娟.陜西省植被覆蓋度變化特征及其成因.應(yīng)用生態(tài)學(xué)報(bào), 2010, 21(11): 2896- 2903.

[27] 牛寶茹, 劉俊蓉, 王政偉.干旱區(qū)植被覆蓋度提取模型的建立.地球信息科學(xué), 2005, 7(1): 84- 97.

[28] 吳云, 曾源, 趙炎, 吳炳方, 武文波.基于MODIS數(shù)據(jù)的海河流域植被覆蓋度估算及動(dòng)態(tài)變化分析.資源科學(xué), 2010, 32(7): 1417- 1424.

[29] 丁立仲, 徐高福, 盧劍波, 章德三, 方炳富.景觀破碎化及其對(duì)生物多樣性的影響.江蘇林業(yè)科技, 2005, 32(4): 45- 57.

[30] 楊芳, 賀達(dá)漢.生境破碎化對(duì)生物多樣性的影響.生態(tài)科學(xué), 2006, 25(6): 564- 567.

[31] 谷景和.新疆東昆侖——阿爾金山地區(qū)的有蹄類(lèi)動(dòng)物.干旱區(qū)研究, 1987, (3): 56- 68.

Studyonspatio-temperalchangeofTibetanAntelope′shabitatbasedonvegetationcoverage

ZHAO Haidi1, LIU Shiliang1,*, DONG Shikui1, SU Xukun1, ZHANG Xiang2

1SchoolofEnvironment,StateKeyLaboratoryofWaterEnvironmentSimulation,BeijingNormalUniversity,Beijing100875,China2NatureReserveAdministrationofAltunNationalNatureReserve,Korla841000,China

The distribution and dynamics of vegetation can reflect the distribution and changes of wildlife habitat to some extent.Tibetan Antelope, well-known as ‘the numen of plateau′, is an endemic animal on the Qinghai-Tibetan Plateau.However, the survival and distribution of Tibetan Antelope have been threatened by environmental deterioration and illegal poaching.This animal has been listed as an Class I species by the Convention on International Trade in Endangered Species (CITES) and as Category I by the Key Protected Wildlife List of China.This situation has attracted a great deal of attentions from government and conservation professionals.The population distribution of Tibetan Antelope remains unkown in Altun National Nature Reserve.The Altun National Nature Reserve, the largest inland nature reserve of China established in May of 1983, is one of the four “no man′s-lands” in China.Due to its unique natural conditions, this reserve is preserving the rich natural resources, rare wild animals and plants.

In this study, we analyzed the dynamics of vegetation in Altun National Nature Reserve based on remote sensing image data and GIS tool.The Normalized Difference Vegetation Index (NDVI) data were collected in 2000, 2005 and 2010.We merged the regions, which were covered by vegetation in the study years and the suitable habitats were identified through interpreting vegetation types, which can be used by Tibetan Antelopes as grazing pastures.Based on NDVI datasets collected in different years, we calculated the vegetation coverage and analyzed the spatial and temporal variations of vegetation coverage.Moreover, the changes of habitat were analyzed on the basis of vegetation coverage variations.The results showed that the vegetation coverage of the reserve was at a low level, and grids with low vegetation coverage took up about 50% of the regions that were covered by vegetation.The grids with high vegetation coverage were mainly distributed in the southwest of Kardun inspection station.The vegetation coverage in the regions where Tibetan Antelopes mainly lived was relatively higher, and that in the regions where the animal lambed was lower.The upward trends of vegetation coverage over time reflected the increase of the reserve′s carrying capacity for animals, and it may also imply that the area of suitable habitat of this animal increased.The vegetation coverage of alpine steppe dominated by feather grass experienced a significant upward trend, also suggesting that carrying capacity of the reserve for Tibetan Antelope increased.The change of spatial distribution of feather grass steppe coverage indicated that habitat of Tibetan Antelopes might spread from east to southwest.The relationship between the altitude and vegetation coverage was not significant correlated.The study of habitat variations based on vegetation coverage is meaningful and can be used for better protecting this wildlife.Furthermore, we suggested that more factors could be integrated with vegetation coverage to precisely identify suitable habitat for animals, and more field surveys and monitoring should be conducted in the future.

tibetan antelope; Altun National Nature Reserve; NDVI; vegetation coverage; feather grass steppe

環(huán)保公益項(xiàng)目(201209033); “十二五”科技支撐項(xiàng)目(2012BAC01B02)

2013- 09- 19;

2014- 04- 01

10.5846/stxb201309192309

*通訊作者Corresponding author.E-mail: shiliangliu@bnu.edu.cn

趙海迪，劉世梁，董世魁，蘇旭坤，張翔.基于植被覆蓋度的藏羚羊棲息地時(shí)空變化研究.生態(tài)學(xué)報(bào),2014,34(12):3285- 3292.

Zhao H D, Liu S L, Dong S K, Su X K, Zhang X.Study on spatio-temperal change of Tibetan Antelope′s habitat based on vegetation coverage.Acta Ecologica Sinica,2014,34(12):3285- 3292.