甘肅白龍江流域植被覆蓋度及景觀格局變化

高彥凈，謝余初，錢大文，張玲玲，鞏 杰

（蘭州大學(xué) 西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，西部環(huán)境與氣候變化研究院，蘭州730000）

植被是生態(tài)環(huán)境的重要構(gòu)成部分，又是維持生態(tài)環(huán)境，發(fā)揮有效生態(tài)效能的功能體，是衡量自然生態(tài)環(huán)境狀況和性質(zhì)的主要指示物［1］。植被覆蓋度是指植被冠層在地面上的垂直投影面積與土地總面積的百分比，是反映地表植被分布特征的定量指標(biāo)［2］，對(duì)衡量地表植被及區(qū)域生態(tài)環(huán)境變化有著重要指示作用［3－4］。隨著遙感技術(shù)的發(fā)展，植被覆蓋度被廣泛應(yīng)用于植被動(dòng)態(tài)監(jiān)測(cè)、自然災(zāi)害監(jiān)測(cè)以及氣候和水文等研究工作［5］，國(guó)內(nèi)外許多學(xué)者基于植被覆蓋度進(jìn)行了大量不同尺 度 的 研 究［6－12］，其 中，李 登 科 等［7］利 用 MODIS NDVI數(shù)據(jù)分析了陜西省封山育林，退耕還林等一系列生態(tài)建設(shè)工程實(shí)施后植被覆蓋度時(shí)空變化特征。陳效逑等［8］基于AVHRR NDVI數(shù)據(jù)對(duì)內(nèi)蒙古植被帶植被空間分布格局進(jìn)行研究，并分析了典型草原植被帶的年季時(shí)空變化及其與水熱因子的關(guān)系。甘春英等［9］以TM影像為數(shù)據(jù)源，分析了連江流域1988—2006年植被演變特點(diǎn)及空間分布特征，發(fā)現(xiàn)連江流域植被覆蓋度變化顯著。以上研究包括了干旱和濕潤(rùn)地區(qū)的植被覆蓋度變化研究，且多從時(shí)間尺度上考慮年際和（或）月間的植被覆蓋度動(dòng)態(tài)變化，而從景觀生態(tài)學(xué)角度對(duì)流域植被景觀格局及其時(shí)空變化過程研究相對(duì)較少，尤其是對(duì)于處在黃土高原、青藏高原和秦巴山區(qū)三大地形交匯過渡帶的植被覆蓋及其景觀格局變化研究就更少。因此，加強(qiáng)對(duì)生態(tài)過渡帶在流域尺度上的景觀格局演變過程研究，對(duì)西部生態(tài)環(huán)境的恢復(fù)和生態(tài)可持續(xù)發(fā)展的進(jìn)程具有重要意義。

白龍江流域地處甘肅南部，是長(zhǎng)江上游重要的水源地和生態(tài)屏障，同時(shí)也是我國(guó)滑坡、泥石流災(zāi)害四大高發(fā)區(qū)之一［13］。流域內(nèi)地形復(fù)雜，地質(zhì)災(zāi)害頻繁、人類活動(dòng)干擾劇烈，生態(tài)環(huán)境十分脆弱。尤其是在八十年代以來，由于不合理的土地利用方式衍生出了一系列生態(tài)環(huán)境問題，諸如植被破壞、水土流失、滑坡、泥石流等災(zāi)害；而植被是制約滑坡、泥石流等地質(zhì)災(zāi)害的重要因子和生態(tài)環(huán)境的關(guān)鍵要素。1998年長(zhǎng)江洪澇災(zāi)害后，我國(guó)相繼推行天然林資源保護(hù)、退耕還林還草等一系列重大生態(tài)建設(shè)工程，白龍江流域也率先開展了生態(tài)恢復(fù)的植被建設(shè)工作。然而，由于這些工程涉及面廣，工程量大，如何客觀、準(zhǔn)確地了解植被生態(tài)恢復(fù)效果成為難題。遙感與地理信息系統(tǒng)技術(shù)對(duì)植被變化監(jiān)測(cè)具有觀測(cè)尺度大、觀測(cè)時(shí)間序列長(zhǎng)、間隔周期短的特點(diǎn)，因此基于歸一化植被指數(shù)NDVI（Normal Differential Vegetation Index）可以在大區(qū)域尺度上進(jìn)行地表植被及其景觀格局變化研究。基于此，本文擬以甘肅白龍江流域?yàn)槔肧POT NDVI數(shù)據(jù)以像元二分模型、景觀格局分析為方法基礎(chǔ)，分析1998—2012年流域植被覆蓋動(dòng)態(tài)變化及其景觀格局變化特征，旨在為流域植被恢復(fù)、生態(tài)環(huán)境改善和流域水土保持供科學(xué)依據(jù)。

1 研究區(qū)概況

甘肅白龍江流經(jīng)迭部縣、舟曲縣、宕昌縣、武都區(qū)和文縣（圖1），全長(zhǎng)475km，流域面積約1.84萬km2。地理位置32°36′—34°24′N，102°48′—105°30′E。流域內(nèi)山巒疊嶂、溝壑縱橫、坡陡流急，素以“山大溝深”著稱，是我國(guó)水土保持與水源涵養(yǎng)的重要功能區(qū)。流域內(nèi)氣候類型復(fù)雜多樣，屬于暖溫帶濕潤(rùn)氣候和北亞熱帶濕潤(rùn)氣候交錯(cuò)帶，夏季高溫多雨，冬季溫暖少雨，年均降水量400～850mm，且集中在5—10月，年均氣溫6～14.9℃。研究區(qū)植被覆蓋較好，森林廣布、生物資源薈萃，是甘肅唯一的油橄欖、茶葉、銀杏等亞熱帶作物產(chǎn)地。由于海拔高，植被的垂直分異明顯：海拔在800～3 700m之間的區(qū)域?yàn)閬喐呱结樔~林；海拔在3 700～4 100m間的區(qū)域多為由杜鵑和高山柳組成的高山灌木林；海拔4 100m以上的高山通常分布有高山草甸或者基巖。研究區(qū)總?cè)丝?30多萬人，其中武都和宕昌的人口數(shù)最多，有漢族、藏族、回族等多個(gè)民族。

圖1 研究區(qū)地理位置

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源與處理

研究采用的數(shù)據(jù)由兩部分組成：（1）比利時(shí)弗萊芒技術(shù)研究所網(wǎng)站（http：∥free.vgt.vito.be）提供的1998—2012年的SPOT NDVI數(shù)據(jù)，空間分辨率為1km，時(shí)間分辨率為逐旬。獲得的SPOT NDVI遙感數(shù)據(jù)已經(jīng)過大氣校正、輻射校正等預(yù)處理，通過MVC（maximum value composite）獲得年最大 NDVI，以減少云、大氣等噪聲對(duì)數(shù)據(jù)的影響。（2）氣象數(shù)據(jù)為白龍江流域及周邊9個(gè)地面氣象站點(diǎn)（武都區(qū)、宕昌縣、舟曲縣、迭部縣、岷縣、禮縣、卓尼縣、康縣、文縣）的1998—2012年逐月降水和氣溫?cái)?shù)據(jù)。運(yùn)用克里格插值法對(duì)氣象數(shù)據(jù)進(jìn)行插值，并運(yùn)用偏相關(guān)分析法，在ArcGIS柵格計(jì)算器中分別計(jì)算植被覆蓋度與降水和氣溫的相關(guān)系數(shù)。

2.2 植被覆蓋度計(jì)算與分級(jí)

研究區(qū)植被覆蓋度提取是基于像元二分模型原理，假設(shè)像元由裸土和植被兩部分構(gòu)成，通過NDVI計(jì)算，進(jìn)而推導(dǎo)出植被覆蓋度［14］，計(jì)算公式如下：

式中：FC——一個(gè)像元內(nèi)的植被覆蓋度；NDVI——該像元上的歸一化植被指數(shù)；NDVIveg，NDVIsoil——影像上植被和裸土部分對(duì)應(yīng)的NDVI值。理論上，大多數(shù)裸土所對(duì)應(yīng)的NDVIsoil應(yīng)該接近0，但實(shí)際上由于受眾多因素（大氣、地表濕度、土壤類型、地表粗糙度等）的影響，NDVIsoil也會(huì)隨著空間而變化［14］。NDVIveg也會(huì)因植被類型的不同等因素而改變。因此，取給定置信度區(qū)間的最大值與最小值，分別選擇累積的5%和95%為置信區(qū)間，反復(fù)對(duì)比分析后選定NDVIveg和NDVIsoil的數(shù)值，進(jìn)而計(jì)算植被覆蓋度分布狀況。同時(shí)，參考中華人民共和國(guó)水利部2008年頒布的《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》（SL190—2007），并結(jié)合研究區(qū)實(shí)際情況，以TM影像為基礎(chǔ)，通過目視解譯得到的2002年、2010年土地利用數(shù)據(jù)，將植被覆蓋度劃分為：FC＜30%（低植被覆蓋度）、30%≤FC＜45%（中低植被覆蓋度）、45%≤FC＜60%（中植被覆蓋度）、60%≤FC＜75%（中高覆蓋度）和FC≥75%（高覆蓋度）5個(gè)等級(jí)［11－12，15］。

2.3 植被覆蓋景觀格局分析

景觀指數(shù)不僅高度濃縮景觀格局信息，而且可以反映景觀結(jié)構(gòu)和空間配置特征變化，進(jìn)而建立景觀結(jié)構(gòu)與過程或現(xiàn)象的聯(lián)系，更好地解釋與理解景觀功能［16］。根據(jù)景觀指數(shù)的生態(tài)學(xué)意義和指數(shù)間的相關(guān)性，按照既可以描述景觀要素又可以描述景觀總體特征的原則［17］，并結(jié)合白龍江流域植被覆蓋特點(diǎn)，擬從植被覆蓋的破碎度、多樣性和景觀空間構(gòu)型角度分別選取斑塊密度（PN）、斑塊占景觀面積的比例（PLAND）、景觀形狀指數(shù)（LSI）、散布與并列指數(shù)（IJI）、香農(nóng)多樣性指數(shù)（SHDI）、香農(nóng)均勻度指數(shù)（SHEI）、面積加權(quán)分維數(shù)（AWMPFD）等指數(shù)進(jìn)行分析，各指標(biāo)表達(dá)公式及其生態(tài)意義見文獻(xiàn)［16，18］。

3 結(jié)果與分析

3.1 1998－2012年流域植被覆蓋度時(shí)空變化

由圖2a和圖3可知，高植被覆蓋和中高植被覆蓋是白龍江流域的主要植被覆蓋類型，二者面積之和均占流域總面積的66.67%以上。1998—2005年，面積增加最大的是高植被覆蓋度（650.56km2），中植被覆蓋度次之（342.14km2）；而面積減少最明顯的是低植被覆蓋度（－621.64km2），中高植被覆蓋度次之（－475.7km2）。從不同等級(jí)植被覆蓋變化程度上看，低植被覆蓋動(dòng)態(tài)度最大（－2.87%），中植被次之（2.73%），高植被最小（1.09%）。2005—2012年，中高植被覆蓋度的面積及動(dòng)態(tài)度最大（分別為482.58km2和2.09%），中低植被覆蓋次之（分別為－192.76km2和－2.02%），高植被覆蓋變化最小（分別為－86.74km2和－0.16%）。整個(gè)研究期間，高植被覆蓋度面積增加最大（737.30km2），而低被覆蓋度面積減少最多（－814.40km2），表明流域內(nèi)植被生態(tài)恢復(fù)明顯，植被覆蓋程度不斷增強(qiáng)。

圖2 1998－2012年各植被覆蓋度等級(jí)的面積變化及其動(dòng)態(tài)度

從空間分布上來看，高、中高植被覆蓋度主要分布在文縣、拱壩河上中游、迭部縣南部、岷江西岸至迭山之間的區(qū)域，同時(shí)在武都區(qū)東南面有少量分布。低和中低植被覆蓋度緊密相連，且呈現(xiàn)多核地域分布格局，具體集中在舟曲—武都區(qū)段白龍江兩岸及其以北區(qū)域、宕昌縣北部及岷江東部、迭部縣北部等區(qū)域，這些地區(qū)除迭部縣北部因海拔高出露大片裸巖之外，主要是坡耕地分布區(qū)域。中植被覆蓋區(qū)分布相對(duì)比較零散，大多分布在低和中低植被覆蓋度的外圍，以舟曲—武都區(qū)段白龍江兩岸及其以北區(qū)域分布面積最大（圖3）。在整個(gè)研究期間，白龍江流域植被覆蓋度分布格局變化不大，總體呈現(xiàn)植被轉(zhuǎn)好的趨勢(shì)。

圖3 甘肅白龍江流域植被覆蓋度空間分布

3.2 1998－2012年流域不同等級(jí)植被覆蓋度變化情況

從不同植被覆蓋度等級(jí)之間面積動(dòng)態(tài)轉(zhuǎn)化來看（表1），1998—2005年，白龍江流域植被覆蓋度提升，首先因?yàn)榇罅恐懈咧脖桓采w轉(zhuǎn)化為高植被覆蓋，轉(zhuǎn)化面積為1 460.31km2；其次為低植被覆蓋向中低植被覆蓋轉(zhuǎn)化、中植被覆蓋向中高植被覆蓋轉(zhuǎn)化，轉(zhuǎn)化面積分別為618.51，535.31km2。這說明研究區(qū)的天然林保護(hù)、退耕還林還草政策實(shí)施初見成效。2005—2012年，植被覆蓋度繼續(xù)提高，而且與1998—2005年植被覆蓋變化規(guī)律一致，即低一等級(jí)的植被覆蓋向高一等級(jí)的植被覆蓋轉(zhuǎn)化的面積最多，其中，中高植被向高植被覆蓋轉(zhuǎn)化的面積最大（1 282.28km2），其次為中植被覆蓋向中高植被覆蓋轉(zhuǎn)化（804.21 km2）。1998—2012年，中高植被覆蓋向較高植被覆蓋轉(zhuǎn)化的面積最大（達(dá)1 676.94km2），其次是中、低植被覆蓋分別向中高、中低植被覆蓋轉(zhuǎn)化；年平均轉(zhuǎn)化幅度最大的是中高植被覆蓋轉(zhuǎn)化為高植被覆蓋，轉(zhuǎn)化幅度約0.65%。這種大面積由低級(jí)植被覆蓋向高級(jí)植被覆蓋轉(zhuǎn)移的現(xiàn)象，對(duì)應(yīng)土地利用變化分析可以發(fā)現(xiàn)，主要是未利用地轉(zhuǎn)變?yōu)楣唷⒉莸兀赂睾臀蠢玫叵蛄值剞D(zhuǎn)移，說明研究區(qū)植被生態(tài)恢復(fù)工程效果顯著，生態(tài)環(huán)境有所改善。

然而，2012年低和中低植被覆蓋度面積比重達(dá)18.63%，說明近15a來白龍江流域森林保護(hù)與植被恢復(fù)工作依然很嚴(yán)峻。

表1 各類型植被覆蓋度轉(zhuǎn)移矩陣

為了更好地反映流域內(nèi)不同等級(jí)植被覆蓋度的變化程度，在進(jìn)行植被覆蓋度等級(jí)變化轉(zhuǎn)移矩陣分析的基礎(chǔ)上，將植被覆蓋度向高等級(jí)的覆蓋度增加2級(jí)以上（含2級(jí)）定義為顯著增加，增加1級(jí)的為輕度增加；反之，若向低等級(jí)的覆蓋度降低2級(jí)以上則為顯著減少，降低1級(jí)的為輕度減少［10，19］。

由圖4和表1可知，15a來白龍江流域植被生態(tài)恢復(fù)明顯，1998—2005年植被覆蓋度顯著增加區(qū)域主要分布在宕昌縣的北部，2005—2012年主要分布在宕昌縣北部、岷江東岸、舟曲—武都段白龍江兩岸及其以北區(qū)域。這可能與這些區(qū)域原來植被覆蓋度比較低，且是水土保持與退耕還林還草工程的重點(diǎn)實(shí)施治理區(qū)有關(guān)。從面積變化上看，1998—2012年低植被覆蓋度向中植被覆蓋度以上轉(zhuǎn)化是植被覆蓋度顯著增加的主要原因，其轉(zhuǎn)化面積約占總轉(zhuǎn)化面積的36.58%。植被覆蓋度顯著減少的區(qū)域集中于西秦嶺和迭山北部，其高、中高植被覆蓋度分別向中、低植被覆蓋度轉(zhuǎn)化是主導(dǎo)因素。

圖4 甘肅白龍江流域不同等級(jí)植被覆蓋度變化空間分布

3.3 植被覆蓋度景觀格局變化

3.3.1 景觀破碎度 斑塊數(shù)能夠較好地反映斑塊之間的破碎化程度，是反映景觀異質(zhì)性的重要指標(biāo)。1998—2012年，白龍江流域整體景觀破碎度逐漸增強(qiáng)，且不同等級(jí)植被覆蓋度的斑塊數(shù)呈現(xiàn)單一增長(zhǎng)的趨勢(shì)，斑塊數(shù)整體增加了58.23%（圖5），表明流域內(nèi)景觀異質(zhì)性增強(qiáng)，也間接地反映出人類活動(dòng)對(duì)區(qū)域景觀格局有著重要的影響。其中，中高植被覆蓋度破碎化最嚴(yán)重，斑塊數(shù)變化最大，年均增加幅度高達(dá)6.02%，從圖3可以看出主要是高植被覆蓋度覆蓋區(qū)退化成很多中高植被覆蓋度造成的；其次為中低植被覆蓋度（增加61.91%），主要是因?yàn)樯鷳B(tài)工程實(shí)施后，植樹造林，退耕還林還草等措施打破了中低植被覆蓋度連成一片的局面；高植被覆蓋度增加的最少（增加12.5%）。

圖5 1998－2012年甘肅白龍江流域各斑塊類型水平上的景觀格局指數(shù)

3.3.2 景觀多樣性 在景觀水平尺度上，景觀多樣性指標(biāo)直接反映了植被景觀要素的均衡化和空間分布均勻化。Shannon′s多樣性指數(shù)（SHDI）從1998年的1.391略減至1.342，其均勻度指數(shù)（SHEI）也由0.864減至0.834，說明研究區(qū)不同等級(jí)植被覆蓋度趨向于不均衡化，植被景觀空間分布格局向不均勻發(fā)展，高和中高植被覆蓋度的面積比重逐漸增大，且分布范圍廣（表3）。

3.4 植被覆蓋變化原因分析

氣溫與降水是植被生長(zhǎng)的重要影響因素，而近年來，隨著城鎮(zhèn)化、修路和采礦等活動(dòng)的加劇，人類干擾對(duì)植被覆蓋度的影響在逐漸增強(qiáng)。降水量的增加有利于植被覆蓋度的提高，適度的氣溫升高也有利于植被生長(zhǎng)，1998—2012年研究區(qū)降水量和氣溫都呈增加趨勢(shì)（圖6）。

圖7顯示，植被覆蓋度與降水顯著相關(guān)地區(qū)主要分布在文縣的中東部、武都區(qū)的中南部和宕昌縣的北部，而植被覆蓋度和氣溫顯著相關(guān)地區(qū)零星分布在宕昌縣、舟曲縣和迭部縣。此外，近年來，白龍江流域長(zhǎng)江中上游防護(hù)林體系建設(shè)工程、天然林資源保護(hù)和退耕還林還草工程的實(shí)施對(duì)植被覆蓋度的增加起到了促進(jìn)作用。如1989—2008年宕昌縣在實(shí)施“長(zhǎng)江中上游防護(hù)林體系建設(shè)工程”中人工造林、封山育林、撫育幼林等552.6km2，1999—2009武都區(qū)人工造林1122.06km2，退耕還林349.13km2。這說明研究區(qū)植被覆蓋度的增加是自然因素和人為因素共同作用的結(jié)果。

表3 1998－2012年甘肅白龍江流域景觀水平上的景觀格局指數(shù)

圖6 1998－2012年甘肅白龍江流域年平均氣溫和年平均降水量變化趨勢(shì)

從圖6可知，植被覆蓋度顯著增加地區(qū)的植被覆蓋度與降水、氣溫的相關(guān)性均不明顯，說明生態(tài)工程植被是研究區(qū)植被面積增加的主要原因。此外，隨著人口增加，經(jīng)濟(jì)發(fā)展，城鎮(zhèn)化速度加快、工礦建設(shè)及交通基礎(chǔ)設(shè)施的建設(shè)對(duì)研究區(qū)植被覆蓋斑塊數(shù)增加和景觀破碎化程度加重等產(chǎn)生了重要影響［20－21］。

圖7 植被覆蓋度與氣溫和降水的空間相關(guān)性

4 結(jié)論與討論

本文基于遙感和地理信息系統(tǒng)技術(shù)平臺(tái)，利用SPOT NDVI數(shù)據(jù)對(duì)植被覆蓋進(jìn)行了分析，客觀地反映了甘肅白龍江流域植被覆蓋度及其景觀格局變化，揭示了該流域近15a來植被覆蓋度整體提高，植被恢復(fù)明顯，同時(shí)景觀破碎化加重。

從植被覆蓋度面積變化來看，1998—2012年白龍江流域植被覆蓋度明顯增加，生態(tài)環(huán)境有了明顯的改善。2012年植被覆蓋和長(zhǎng)勢(shì)明顯好于1998年，低植被覆蓋度面積減少最多（814.40km2），中高和高植被覆蓋度面積約占研究區(qū)總面積的比例達(dá)70.7%，植被覆蓋狀況明顯好轉(zhuǎn)。

從植被覆蓋度空間變化來看，1998—2012年白龍江流域植被覆蓋的空間分布呈現(xiàn)出一定的規(guī)律性，低和中低植被覆蓋度多集中在人類活動(dòng)相對(duì)頻繁的河谷兩岸及石質(zhì)性山區(qū)與黃土地貌交錯(cuò)過渡的土石山區(qū)，舟曲—武都段白龍江兩岸及其以北區(qū)域尤為突出，其植被景觀格局穩(wěn)定性較差，存在較大的潛在生態(tài)風(fēng)險(xiǎn)，這與趙彩霞［20］、謝余初［21］等的研究結(jié)果類似。可見，流域內(nèi)植被生態(tài)恢復(fù)工作仍然很緊迫，建議研究區(qū)結(jié)合自身實(shí)際情況，堅(jiān)持可持續(xù)發(fā)展戰(zhàn)略，不可以環(huán)境破壞為代價(jià)來發(fā)展經(jīng)濟(jì)。

從植被覆蓋度的景觀格局來看，研究期間，流域植被景觀逐漸破碎化，結(jié)構(gòu)復(fù)雜化，異質(zhì)性和幾何形狀特征逐漸增強(qiáng)，不同覆蓋度等級(jí)的植被景觀混雜程度和受干擾程度持續(xù)增強(qiáng)，且趨向于不均衡化發(fā)展。說明15a來白龍江流域森林保護(hù)與植被恢復(fù)工作依然很嚴(yán)峻，因此，在實(shí)施生態(tài)建設(shè)工程的同時(shí)，應(yīng)注意加強(qiáng)流域生態(tài)系統(tǒng)管理和生態(tài)建設(shè)成果保護(hù)，減少破壞，避免形成“邊恢復(fù)，邊破壞”的現(xiàn)象。

自然因素和人為因素是研究區(qū)植被覆蓋度及景觀格局發(fā)生變化的主要驅(qū)動(dòng)力。研究區(qū)植被覆蓋度顯著增加地區(qū)與氣溫和降水的相關(guān)性不大，說明近15a來，自然因素對(duì)研究區(qū)植被覆蓋度的增加影響較小，而人為因素起到了重要作用。盡管研究區(qū)較早開展了植被建設(shè)工作，“長(zhǎng)江中上游防護(hù)林體系建設(shè)工程”、“天然林資源保護(hù)工程”和“退耕還林還草工程”相繼實(shí)施，但其工程建設(shè)對(duì)植被恢復(fù)的影響及其貢獻(xiàn)率還有待進(jìn)一步研究。

［1］ 孫紅雨，王長(zhǎng)耀，牛錚，等.中國(guó)地表植被覆蓋變化及其與氣候因子關(guān)系［J］.遙感學(xué)報(bào)，1998，2（3）：204－210.

［2］ 胡良軍，邵明安.論水土流失研究中的植被覆蓋度量指標(biāo)［J］.西北林學(xué)院學(xué)報(bào)，2001，16（1）：40－43.

［3］ Sellers P J，Los S O，Tucker C J.A revised land surface parameterization（SiB2）for atmospheric GCMs.PartⅡ：the generation of global fields of terrestrial biophysical parameters from satellite data［J］.Journal of Climate，1996，9（4）：706－737.

［4］ 陳云浩，李曉兵，史培軍，等.北京海淀區(qū)植被覆蓋的遙感動(dòng)態(tài)研究［J］.植物生態(tài)學(xué)報(bào)，2001，25（5）：588－593.

［5］ 吳云，曾源，吳炳方，等.基于MODIS數(shù)據(jù)的三北防護(hù)林工程區(qū)植被覆蓋度提取與分析［J］.生態(tài)學(xué)雜志，2009，28（9）：1712－1718.

［6］ Samir M A，Mohamed A H.Vegetation dynamic in semi arid Butana Plain，Sudan［J］.Journal of Geography and Regional Planning，2011（4）：183－198.

［7］ 李登科，范建忠，王娟.陜西省植被覆蓋度變化特征及其成因［J］.應(yīng)用生態(tài)學(xué)報(bào)，2010，21（11）：2896－2903.

［8］ 陳效逑，王恒.1982—2003年內(nèi)蒙古植被帶和植被覆蓋度的時(shí)空變化［J］.地理學(xué)報(bào)，2009，64（1）：84－94.

［9］ 甘春英，王兮之，李保生，等.連江流域近18a來植被覆蓋度變化分析［J］.地理科學(xué)，2011，31（8）：1019－1024.

［10］ 鄧飛，全占軍，于云江，等.20a來烏蘭木倫河流域植被蓋度變化及影響因素［J］.水土保持研究，2011，18（3）：137－140.

［11］ 王兮之，陸龍，梁釗雄，等.青海湟水流域植被覆蓋度時(shí)空變化分析［J］.水土保持研究，2010，17（6）：172－176.

［12］ 龐吉林，張克斌，喬娜，等.基于RS的鹽池縣近10a植被覆蓋度動(dòng)態(tài)變化研究［J］.水土保持研究，2012，19（4）：112－115，121.

［13］ 陳冠，孟興民，喬良，等.“7·22”岷縣漳縣地震地質(zhì)災(zāi)害分布、特征及與影響因子間關(guān)系分析［J］.工程地質(zhì)學(xué)報(bào)，2013，21（5）：750－760.

［14］ 李苗苗，吳炳方，顏長(zhǎng)珍，等.密云水庫上游植被覆蓋度的遙感估算［J］.資源科學(xué)，2004，26（4）：153－159.

［15］ 王冬梅，孟興民，邢釗，等.基于RS的武都區(qū)植被覆蓋度動(dòng)態(tài)變化及其驅(qū)動(dòng)力分析［J］.干旱區(qū)資源與環(huán)境，2012，26（11）：92－97.

［16］ 鄔建國(guó).景觀生態(tài)學(xué)：格局、過程、尺度與等級(jí)［M］.北京：高等教育出版社，2007.

［17］ 陳文波，肖篤寧，李秀珍，等.景觀指數(shù)分類、應(yīng)用及構(gòu)建研究［J］.應(yīng)用生態(tài)學(xué)報(bào)，2002，13（1）：121－125.

［18］ McGarigal K，Cushman S A，Neel M C，et al.FRAGSTATS：Spatial Pattern Analysis Program for Categorical Maps［M］.Massachusetts：University of Massachusetts，Amherst，2002.

［19］ 劉憲鋒，楊勇，任志遠(yuǎn)，等.2000—2009年黃土高原地區(qū)植被覆蓋度時(shí)空變化［J］.中國(guó)沙漠，2013，33（4）：1244－1249.

［20］ 趙彩霞.甘肅白龍江流域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［D］.蘭州：蘭州大學(xué)，2013.

［21］ 謝余初，鞏杰，趙彩霞.甘肅白龍江流域水土流失的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［J］.生態(tài)學(xué)雜志，2014，33（3）：702－708.