候 靜， 杜靈通， 馬 菁， 張學(xué)儉，

(1.寧夏大學(xué) 西北土地退化與生態(tài)恢復(fù)省部共建國家重點(diǎn)實(shí)驗(yàn)室培育基地/西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點(diǎn)實(shí)驗(yàn)室/西部生態(tài)與生物資源開發(fā)聯(lián)合研究中心， 寧夏 銀川 750021； 2.寧夏農(nóng)林科學(xué)院， 寧夏 銀川 750021)

基于RS與像元二分模型的近20 a寧夏植被覆蓋研究

候 靜1， 杜靈通1， 馬 菁2， 張學(xué)儉1，2

(1.寧夏大學(xué) 西北土地退化與生態(tài)恢復(fù)省部共建國家重點(diǎn)實(shí)驗(yàn)室培育基地/西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點(diǎn)實(shí)驗(yàn)室/西部生態(tài)與生物資源開發(fā)聯(lián)合研究中心， 寧夏 銀川 750021； 2.寧夏農(nóng)林科學(xué)院， 寧夏 銀川 750021)

[目的] 探討近20 a寧夏植被變化特征及其對(duì)氣候的響應(yīng)情況，為了解區(qū)域生態(tài)變化態(tài)勢(shì)，評(píng)估生態(tài)工程成效提供理論基礎(chǔ)和科學(xué)依據(jù)。 [方法] 利用1992，2000，2006和2012年的TM數(shù)據(jù)，基于植被指數(shù)估算植被覆蓋度的原理，運(yùn)用像元二分模型，對(duì)寧夏植被覆蓋度進(jìn)行遙感估算和動(dòng)態(tài)監(jiān)測。 [結(jié)果] 寧夏植被覆蓋呈現(xiàn)南北好，中部差的空間分布特點(diǎn)；近20 a來寧夏的植被覆蓋度表現(xiàn)出明顯的增加趨勢(shì)，即低植被覆蓋區(qū)的面積在逐漸減少，中、高植被覆蓋區(qū)的面積增加了6 434 km2，全區(qū)平均植被覆蓋度提高6%，其中揚(yáng)黃灌區(qū)及南部山區(qū)植被覆蓋增加幅度較大。 [結(jié)論] 寧夏植被覆蓋度的增加主要得益于近年來實(shí)施的一系列植被保護(hù)和生態(tài)恢復(fù)工程，而降水量的變化則是引起植被覆蓋度年際波動(dòng)的主要因素。

遙感； 植被覆蓋度； 像元二分模型； 寧夏

植被覆蓋度是指單位面積內(nèi)植被地上部分的垂直投影面積所占的百分比值。它是衡量地表植被狀況的一個(gè)重要指標(biāo)，直接影響著地表溫度，土壤濕度，地表能量與水循環(huán)等生態(tài)因子[1]。傳統(tǒng)測量植被覆蓋度的方法有目估法[2]、樣方法[3]、儀器法等[4]，這些方法由于存在主觀性強(qiáng)、成本高和操作不便等限制，因此難以大范圍提取植被覆蓋度。遙感技術(shù)以其多平臺(tái)、多波段、短周期、大范圍等特點(diǎn)，已被廣泛的運(yùn)用于環(huán)境評(píng)估、資源開發(fā)、生態(tài)監(jiān)測等諸多領(lǐng)域，為實(shí)時(shí)、連續(xù)、大范圍的監(jiān)測植被覆蓋度提供了技術(shù)支撐[5]。常用于遙感監(jiān)測植被覆蓋度的植被指數(shù)有NDVI，PVI，SAVI，RVI，VCI等，其中植被歸一化指數(shù)NDVI(normalized difference vegetation index)由于檢測靈敏度高，范圍較寬，可消除地形、陰影干擾和削弱太陽高度角及大氣帶來的噪音等，常被用來反映植物生長狀況和空間分布密度等[6-7]。像元二分模型是一種簡單的線性像元分解模型，它通過分解像元信息建立模型，估算植被覆蓋度[8]。Qi J等[9]運(yùn)用此模型研究了美國SanPedro盆地地區(qū)的植被時(shí)空動(dòng)態(tài)變化。李苗苗等[8]對(duì)模型進(jìn)行改進(jìn)，估算了密云水庫上游的植被覆蓋度。

寧夏回族自治區(qū)位于中國從半濕潤向半干旱干旱區(qū)的過渡地帶，從農(nóng)業(yè)生產(chǎn)向牧業(yè)生產(chǎn)的過渡地帶，該區(qū)自然條件和資源具有多樣性、苛刻性和脆弱性[10-11]的特點(diǎn)，因此，該區(qū)的植被變化情況已成為諸多學(xué)者研究的熱點(diǎn)[12-16]。杜靈通等[12]借助SPOT-VGT 1 km NDVI數(shù)據(jù)探討了寧夏植被覆蓋動(dòng)態(tài)變化與氣候因子的關(guān)系；陳曉光等[13]利用GIMMS 8 km NDVI反演了寧夏近20 a的植被覆蓋度；范錦龍等[14]采用PATHFINDER 8 km NDVI數(shù)據(jù)對(duì)該區(qū)植被覆蓋進(jìn)行了時(shí)序分析。而這些研究大多采用低空間分辨率遙感產(chǎn)品。運(yùn)用像元二分模型與植被指數(shù)結(jié)合，對(duì)寧夏地區(qū)植被覆蓋進(jìn)行反演的研究尚不多見。鑒于此，筆者借助空間分辨率為30 m的Landsat 5 TM影像，根據(jù)植被指數(shù)估算植被覆蓋度原理，運(yùn)用改進(jìn)的像元二分模型分解像元信息，計(jì)算植被覆蓋地表與無植被覆蓋地表分別對(duì)像元的貢獻(xiàn)，并利用研究區(qū)土地利用類型確立閾值，估算了寧夏1992，2000，2006和2010年的植被覆蓋度。本研究結(jié)合寧夏三大自然區(qū)地理特征，對(duì)4期植被覆蓋度數(shù)據(jù)進(jìn)行對(duì)比分析，動(dòng)態(tài)監(jiān)測近20 a寧夏植被覆蓋度變化，旨在為把握區(qū)域生態(tài)變化態(tài)勢(shì)、評(píng)估生態(tài)工程成效提供理論基礎(chǔ)和科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

寧夏回族自治區(qū)位于東經(jīng)104°17′—107°39′，北緯35°14′—39°23′，西側(cè)、北側(cè)至東側(cè)由騰格里沙漠、烏蘭布和沙漠和毛烏素沙地相圍，南與黃土高原相連。屬于典型的溫帶大陸性氣候，年均溫7 ℃，有效積溫3 300 ℃，無霜期150 d，年降水量183～677 mm，降水主要集中于7—8月，年蒸發(fā)量1 214～2 803 mm。全區(qū)總面積5.18×104km2，海拔1 090～2 900 m，南北狹長，地勢(shì)南高北低，西部高差較大，東部起伏較緩。土壤類型以灰鈣土、灌淤土、草甸土為主。依據(jù)不同的自然地理?xiàng)l件可將全區(qū)劃分為引黃灌區(qū)、中部干旱帶及南部山區(qū)[12]。引黃灌區(qū)主要集中在北部的銀川平原，是寧夏最為富饒的農(nóng)業(yè)種植區(qū)。寧夏中部多為丘陵山地和山間盆地，地表景觀以荒漠草原為主。南部山區(qū)屬于黃土高原的一部分，特別是六盤山地區(qū)，海拔較高氣候濕潤，植被豐富。

1.2 數(shù)據(jù)準(zhǔn)備

根據(jù)研究的空間尺度和時(shí)間尺度，選取1992，2000，2006和2010年的Landsat 5 TM影像，空間分辨率30 m，由于研究區(qū)范圍較大，需四景影像拼接，軌道號(hào)與獲取時(shí)間如表1所示。TM數(shù)據(jù)具有較高的空間分辨率，但時(shí)間分辨率較低，而研究區(qū)不處于同一軌道，受云雨天氣等因素影響，難以獲取同一天寧夏全境的影像。為減少因時(shí)間差異造成的影響，研究選取植物生長季數(shù)據(jù)，并在影像拼接過程中重疊區(qū)均采用時(shí)序相近的條帶號(hào)(PATH)129的影像數(shù)據(jù)覆蓋。

表1 Landsat 5 TM影像的軌道號(hào)及其獲取時(shí)間

為消除太陽高度角、大氣散射和地形陰影等噪聲的影響，還原影像真實(shí)的光譜信息，首先對(duì)影像進(jìn)行輻射定標(biāo)，即通過影像灰度值(digital number，DN)值來計(jì)算輻射強(qiáng)度，Landsat 5 TM定標(biāo)計(jì)算公式為[17]：

Lλ=gain×Qλ+offset

(1)

式中：λ——波段值；Lλ——遙感器所接收到的輻射強(qiáng)度；Qλ——以DN表示的經(jīng)過量化標(biāo)定的像元值； gain——衛(wèi)星的增益系數(shù)； offset——偏移系數(shù)。在ENVI 5.0環(huán)境下進(jìn)行定標(biāo)時(shí)，軟件可以自動(dòng)從頭文件中讀取增益和偏移系數(shù)。

定標(biāo)之后采用FLAASH大氣校正，設(shè)置研究方法成像中心點(diǎn)經(jīng)緯度、傳感器高度、成像區(qū)域平均高度、成像時(shí)間設(shè)置，大氣模型MLS，氣溶膠模型Rural，氣溶膠反演方法2 Band(K-T)，初始能見度40 km。將計(jì)算結(jié)果鑲嵌，依研究區(qū)矢量數(shù)據(jù)進(jìn)行裁剪。

將寧夏第二次土地利用普查結(jié)果與1∶10萬土地利用類型矢量數(shù)據(jù)[18]生成掩膜文件，數(shù)據(jù)來源于“寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心”(http:∥westdc.westgis.ac.cn)。依照國家土地利用類型現(xiàn)狀分類，保留水體、林地、耕地、城鎮(zhèn)、草地、未利用地將比例較少的利用類型合并。并由中國氣象數(shù)據(jù)共享平臺(tái)獲取當(dāng)年研究區(qū)10個(gè)氣象站的年際降水量數(shù)據(jù)。

1.3 研究方法

歸一化植被指數(shù)(NDVI)是由傳感器所接收的地物光譜信息推算得到的反映地表植被狀況的定量值[8]，如式(2)：

NDVI=(NIR-R)/(NIR+R)

(2)

式中：R——第3波段紅波段(0.63～0.69 μm)； NIR——Landsat第4波段近紅外波段(0.76～0.90 μm)。

像元二分模型假設(shè)像元由植被覆蓋地表和無植被覆蓋地表兩部分組成，這兩個(gè)組分對(duì)像元信息都有貢獻(xiàn)。像元的光譜信息即由這兩部分的面積比例加權(quán)線性組合而成。即：

S=Fc·Sveg+(1-Fc)Ssoil

(3)

將公式(3)變換，則植被覆蓋度Fc可表示為：

Fc=(S-Ssoil)/(Sveg-Ssoil)

(4)

將植被歸一化指數(shù)NDVI，代入公式(4) ，計(jì)算植被覆蓋度Fc：

Fc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

(5)

式中：NDVIsoil——完全是裸土或無植被覆蓋區(qū)域的NDVI值； NDVIveg——完全被植被所覆蓋的像元的NDVI值，即純植被像元的NDVI值。受眾多因素影響，NDVIsoil并不是理論上的零值，而在-0.1～0.2之間變化[19-20]，NDVIveg隨植被類型的差異而變化，因此NDVIveg與NDVIsoil的取值并不固定。實(shí)際運(yùn)用中，相同的土地利用類型植被類型相似，為此，研究借助各土地利用類型掩膜數(shù)據(jù)來確定NDVIsoil，NDVIveg的閾值。在缺少實(shí)測數(shù)據(jù)的情況下，將Fcmax近似取100%，F(xiàn)cmin近似取0%，則：

NDVIsoil=NDVImin； NDVIveg=NDVImax

為了消除噪聲對(duì)影像準(zhǔn)確性的影響，在計(jì)算最大值與最小值時(shí)設(shè)置置信區(qū)間，將像元累計(jì)直方圖拐點(diǎn)作為置信區(qū)間的端點(diǎn)。經(jīng)比較分析，選取NDVI累計(jì)頻率2%，98%的值作為NDVIsoil，NDVIveg的參數(shù)值代入公式(5)計(jì)算植被覆蓋度。

依據(jù)各區(qū)矢量文件分別統(tǒng)計(jì)3大自然區(qū)及全區(qū)的植被覆蓋度，結(jié)合當(dāng)年降水?dāng)?shù)據(jù)建立變量間Pearson相關(guān)系數(shù)。并對(duì)1992和2010年的植被覆蓋情況進(jìn)行專題變化監(jiān)測，更直觀的反映這20 a間研究區(qū)植被覆蓋變化情況。

2 結(jié)果與分析

2.1 植被覆蓋變化特征及其與降水的關(guān)系

根據(jù)植被指數(shù)估算植被覆蓋度原理，選取1992，2000，2006和2010年的TM影像，利用像元二分模型，計(jì)算出寧夏當(dāng)年的植被覆蓋度，剔除無植被覆蓋的水體和異常值，對(duì)影像進(jìn)行密度分割(圖1)。參考前人對(duì)干旱區(qū)植被覆蓋度的劃分標(biāo)準(zhǔn)[14]，將植被覆蓋度劃分為3類，即：0%～30%為低覆蓋，30%～60%為中覆蓋，60%～100%為高覆蓋，據(jù)此統(tǒng)計(jì)4 a不同植被覆蓋度的變化，結(jié)果如圖2所示。

圖1 寧夏1992，2000，2006和2010年植被覆蓋度分級(jí)

圖2 各等級(jí)植被覆蓋區(qū)面積及年均降水量統(tǒng)計(jì)

就空間分布來看，研究區(qū)高植被覆蓋地區(qū)主要分布于南北兩側(cè)，中部地區(qū)植被覆蓋程度較低。北部的銀川平原雖氣候干旱，但受黃河灌溉之利，成為了寧夏最為富饒的農(nóng)業(yè)灌區(qū)，因此植被覆蓋率高。南部的六盤山林區(qū)氣候濕潤，適宜林木生長，以天然次生闊葉林為主，植被覆蓋程度也較高。中部地區(qū)降水稀少且無大的河流流經(jīng)，植被類型以荒漠草原為主，除揚(yáng)黃灌區(qū)外，植被覆蓋都較低。這與寧夏3區(qū)自然地理特征相一致。就時(shí)間序列而言，4 a植被覆蓋度分別38.5%，37.2%，41.2%和44.8%。從1992—2010年全區(qū)植被覆蓋度提高了6%，低覆蓋區(qū)面積下降，中高覆蓋區(qū)面積增加了6 434 km2，植被整體恢復(fù)狀況良好。

研究區(qū)地處內(nèi)陸，氣候干旱，植被的生長狀況深受當(dāng)年降水影響，尤其在荒漠草原與部分旱作農(nóng)業(yè)區(qū)，在降水豐沛的年度，草原植被與農(nóng)作物長勢(shì)良好，地表植被覆蓋也會(huì)增加。為此，筆者運(yùn)用相關(guān)分析法來計(jì)算研究區(qū)植被覆蓋度與降水的關(guān)系。目前中國氣象站站點(diǎn)獲取多為城市氣象數(shù)據(jù)，城市周邊植被覆蓋對(duì)研究區(qū)而言代表性較低，為減少隨機(jī)誤差，本研究以自然區(qū)為單位，平均各區(qū)降水與植被覆蓋數(shù)據(jù)，分別統(tǒng)計(jì)3大自然區(qū)的植被覆蓋度與降水情況(表2)，并進(jìn)行Pearson相關(guān)分析(n=12)。相關(guān)分析結(jié)果通過了α=0.05的顯著性檢驗(yàn)，相關(guān)系數(shù)為0.58，表明研究區(qū)植被覆蓋度與降水呈中度正相關(guān)。

表2 3大自然區(qū)植被覆蓋度及降水量統(tǒng)計(jì)

2.2 植被覆蓋動(dòng)態(tài)監(jiān)測

為了能更好的把握研究區(qū)生態(tài)變化態(tài)勢(shì)、評(píng)估生態(tài)工程成效，本研究對(duì)近20 a間寧夏植被覆蓋進(jìn)行了專題變化監(jiān)測。參考前述的植被覆蓋度劃分標(biāo)準(zhǔn)，對(duì)1992和2010年的植被覆蓋圖做密度分割，將2 a植被覆蓋圖疊加，監(jiān)測植被變化情況(圖3)，并采用轉(zhuǎn)移矩陣數(shù)學(xué)模型，分別統(tǒng)計(jì)3大自然區(qū)的植被覆蓋變化轉(zhuǎn)移矩陣，如表3所示。

由圖3，表3可知，研究區(qū)整體植被覆蓋度有明顯增加。從整體變化面積來看，低覆蓋轉(zhuǎn)向中高覆蓋的區(qū)域面積超過植被退化面積。從地理區(qū)劃來看，中部干旱帶及南部山區(qū)的植被增加更為明顯。這很大程度是是得益于該區(qū)的生態(tài)工程建設(shè)的有效實(shí)施。中部干旱帶的植被變化區(qū)域主要集中于揚(yáng)黃灌溉地區(qū)。1998年，扶貧揚(yáng)黃灌溉工程開始在寧夏中部地區(qū)實(shí)施，開發(fā)灌溉面積5.30×104hm2，使得該區(qū)域農(nóng)業(yè)建設(shè)快速發(fā)展。其中，清水河流域一線工程建設(shè)最多，農(nóng)田得到開墾，植被覆蓋也相應(yīng)提高，由低覆蓋轉(zhuǎn)向高覆蓋。南部山區(qū)的植被恢復(fù)很大程度上得益于國家天然林資源保護(hù)工程的實(shí)施。2000年，六盤山全面禁伐封育，經(jīng)過十幾年間封山育林、人工造林，林木生長量進(jìn)一步加快，植被覆蓋度明顯提升。此外，退牧還草工程、賀蘭山天然林保護(hù)工程、黃河護(hù)岸林等工程的實(shí)施都對(duì)全區(qū)的植被恢復(fù)起到了重要作用，改善了寧夏生態(tài)環(huán)境的狀況。

圖3 寧夏1992與2010年植被覆蓋變化區(qū)域表3 1992與2010年植被覆蓋轉(zhuǎn)移矩陣

km2

3 結(jié)論與討論

本文借助30 m空間分辨率Landsat 5 TM影像對(duì)寧夏近20 a植被覆蓋度進(jìn)行遙感評(píng)估與動(dòng)態(tài)監(jiān)測，較其他低分辨率遙感產(chǎn)品空間精度略高，有利于定量統(tǒng)計(jì)本區(qū)植被覆蓋特征。通過改進(jìn)的像元二分模型分解像元信息，根據(jù)研究區(qū)土地利用類型分別統(tǒng)計(jì)各種地類的NDVIsoil，NDVIveg取值，與一般反演植被覆蓋算法相比較，其結(jié)果更為準(zhǔn)確[10]。研究結(jié)合寧夏3大自然區(qū)具體地理特征，通過對(duì)4期植被覆蓋度數(shù)據(jù)進(jìn)行對(duì)比分析，結(jié)果顯示：

(1) 研究區(qū)植被覆蓋存在南北好，中部差的分布特點(diǎn)。這種格局與該區(qū)的自然地理區(qū)劃相吻合。銀川平原受黃河灌溉之利，農(nóng)業(yè)發(fā)達(dá)，植被覆蓋率高；南部的六盤山林區(qū)氣候濕潤，林木茂盛；中部地區(qū)降水稀少且無大的河流流經(jīng)，除揚(yáng)黃灌區(qū)外，多為典型的荒漠化草原和退化干草原，植被覆蓋較低，生態(tài)環(huán)境脆弱。

(2) 近20 a來寧夏的植被覆蓋度表現(xiàn)出明顯的增加趨勢(shì)，植被整體恢復(fù)狀況良好。1992，2000，2006和2010年植被覆蓋度分別38.5%，37.2%，41.2%和44.8%，植被低覆蓋區(qū)的面積逐漸減少，植被中、高覆蓋區(qū)的面積增加了6 434 km2，全區(qū)平均植被覆蓋度提高6%，其中揚(yáng)黃灌區(qū)及南部山區(qū)植被覆蓋增加幅度較大。

(3) 引起研究區(qū)植被覆蓋度變化的原因包括人為因素和自然因素。其中，人類活動(dòng)是影響植被覆蓋度變化的主要原因。20世紀(jì)寧夏全區(qū)隨著人口的快速增長，農(nóng)牧民過度墾殖、放牧，毀林毀草開荒嚴(yán)重，造成了生態(tài)環(huán)境的極大破壞；20世紀(jì)90年代以后，寧夏實(shí)施了許多生態(tài)建設(shè)工程，例如，揚(yáng)黃灌溉工程，退牧還草工程，六盤山、賀蘭山天然林保護(hù)工程等。這些工程的有效實(shí)施對(duì)于植被恢復(fù)起了巨大的作用，一定程度上遏制了土地荒漠化，水土流失等生態(tài)退化問題，促進(jìn)了生態(tài)環(huán)境的改善。此外，氣候也是影響植被覆蓋度變化的重要因素，研究區(qū)屬于溫帶干旱半干旱氣候，降水不足且波動(dòng)性大，降水量的多少直接制約著植被覆蓋度的變化。

受其他條件限制，本文只選取了1992，2000，2006和2010年4個(gè)時(shí)間點(diǎn)，因此，對(duì)時(shí)段內(nèi)部的精細(xì)變化仍需更深的研究。借助TM 30 m遙感影像反演植被反演結(jié)果在空間分辨率上具有較高精度，但由于影像獲取周期較長，個(gè)別影像時(shí)間差異可能造成研究結(jié)果的不確定性，為此，如何結(jié)合不同時(shí)空精度的影像反演植被覆蓋狀況也是今后研究的一個(gè)重要方向。此外，運(yùn)用像元二分模型決定NDVI閾值時(shí)置信區(qū)間的設(shè)置依圖像大小、圖像清晰度等條件不同取值亦不相同，具有一定的主觀性，因此對(duì)置信區(qū)間的選擇有待進(jìn)一步定量分析。

[1] 方精云，樸世龍，賀金生，等.近20 a來中國植被活動(dòng)在增強(qiáng)[J].中國科學(xué)(C輯)，2003，33(6):554-565.

[2] 章文波，符素華，劉寶元.目估法測量植被覆蓋度的精度分析[J].北京師范大學(xué)學(xué)報(bào):自然科學(xué)版，2001，37(3):402-408.

[3] 張?jiān)葡?，李曉兵，陳云?草地植被蓋度的多尺度遙感與實(shí)地測量方法綜述[J].地球科學(xué)進(jìn)展，2003，18(1):85-93.

[4] 趙春玲，李志剛.中德合作寧夏賀蘭山封山育林育草項(xiàng)目區(qū)植被覆蓋度監(jiān)測[J].寧夏農(nóng)林科技，2000(S):6-14.

[5] 馬娜，胡云鋒，莊大方，等.基于遙感和像元二分模型的內(nèi)蒙古正藍(lán)旗植被覆蓋度格局和動(dòng)態(tài)變化[J].地理科學(xué)，2012，32(2):251-256.

[6] 李曉松，李增元，高志海，等.基于NDVI與偏最小二乘回歸的荒漠化地區(qū)植被覆蓋度高光譜遙感估測[J].中國沙漠，2011，31(1):162-167.

[7] 方紅亮，田慶久.高光譜遙感在植被監(jiān)測中的研究綜述[J].遙感技術(shù)與應(yīng)用，1998，13(1):62-69.

[8] 李苗苗.植被覆蓋度的遙感估算方法研究[D].北京:中國科學(xué)院，2003.

[9] Qi J， Marsett R C， Moran M S， et al. Spatial and temporal dynamics of vegetation in the San Pedro River basin area[J]. Agricultural and forest Meteorology， 2000，105(1):55-68.

[10] 張克斌，王錦林，侯瑞萍，等.我國農(nóng)牧交錯(cuò)區(qū)土地退化研究[J].中國水土保持科學(xué)，2003，1(1):85-90.

[11] 張學(xué)儉，馮銳.近10 a寧夏鹽池縣生態(tài)景觀格局動(dòng)態(tài)變化研究[J].寧夏大學(xué)學(xué)報(bào):自然科學(xué)版，2007，27(4):369-372.

[12] 杜靈通，田慶久.寧夏植被覆蓋動(dòng)態(tài)變化及與氣候因子的關(guān)系[J].中國沙漠，2012，32(5):1479-1485.

[13] 陳曉光，李劍萍，韓穎娟，等.寧夏近20 a來植被覆蓋度及其與氣溫降水的關(guān)系[J].生態(tài)學(xué)雜志，2007，26(9):1375-1383.

[14] 范錦龍，張曉煜.寧夏近18 a植被動(dòng)態(tài)變化遙感時(shí)序分析[J].干旱區(qū)研究，2009，1(1)∶53-58.

[15] 杜靈通，李國旗.基于SPOT-VGT的寧夏鹽池縣近8 a生態(tài)環(huán)境動(dòng)態(tài)監(jiān)測[J].北京林業(yè)大學(xué)學(xué)報(bào)，2008，30(5):46-51.

[16] 杜靈通，田慶久.寧夏回族自治區(qū)NDVI的時(shí)空變化特征研究[J].水土保持通報(bào)，2011，31(5):208-214.

[17] Chander G， Markham B L， Helder D L. Summary of current radiometric calibration coefficients for Landsat MSS， TM， ETM+， and EO-1ALI sensors[J]. Remote sensing of Environment， 2009，113(5):893-903.

[18] Liu Jiyuan， Liu Mingliang， Zhuang Dafang， et al. Study on spatial pattern of land-use change in China during 1995—2000[J]. Science in China(Series D:Earth Sciences)，2003，46(4)：373-384.

[19] Rundquist B C. The influence of canopy green vegetation fraction on spectral measurements over native tall grass prairie[J]. Remote Sensing of Environment， 2002，81(1):129-135.

[20] Carlson T N， Ripley D A. On the relation between NDVI， fractional vegetation cover， and leaf area index[J]. Remote Sensing of Environment， 1997，62(3):241-252.

A Study on Vegetation Coverage in Ningxia Hui Autonomous Region During Past 20 Years Based on Remote Sensing and Dimidiate Pixel Model

HOU Jing， DU Lingtong， MA Jing， ZHANG Xuejian

(1.BreedingBaseforStateKeyLaboratoryofLandDegradationandEcologicalRestorationinNorthwestChina/KeyLabofRestorationandReconstructionofDegradedEcosysteminNorthwestChinaofMinistryofEducation/UnionResearchCenterforEcologyandExploitationofBiologicalResourcesinWesternChina，NingxiaUniversity，Yinchuan，Ningxia750021，China； 2.InstituteofAgriculturalScienceandTechnologyInformation，NingxiaAgricultureandForestryAcademy，Yinchuan，Ningxia750021，China)

[Objective] The aims of this study is to investigate the dynamic change of vegetation coverage and its response to climate change during the past 20 years in Ningxia Hui Autonomous Region in order to provide theoretical foundation and scientific basis for understanding the regional ecological change trend and evaluating of ecological engineering. [Methods] Remote sensing data used in the study was Landsat TM in 1992， 2000， 2006 and 2012. Based on dimidiate pixel model， different vegetation indices for estimating vegetation coverage was calculated， and the dynamic changes of vegetation coverage were monitored. [Results] High vegetation cover was found in the northern and southern parts of Ningxia Hui Autonomous Region， in contrast， the lowest vegetation cover was found in the middle part the region. Vegetation cover increased significantly during past 20 years， with the low vegetation cover area decreased gradually while the middle- and high-vegetation cover area increased about 6 434 km2. The average vegetation cover in the whole region increased by 6%. In particular， vegetation cover in the Yellow River irrigation area of Northern Ningxia and mountainous area of Southern Ningxia increased by a large margin. [Conclusion] The vegetation increase in the study area is benefit from the implementation of a series of vegetation protection and ecological restoration projects in recent years. The annual change in vegetation coverage is mainly caused by the variation of precipitation.

remote sensing； vegetation fraction； dimidiate pixel model； Ningxia Hui Autonomous Region

2014-08-29

2014-12-02

國家自然科學(xué)基金項(xiàng)目“干旱半干旱地區(qū)農(nóng)業(yè)干旱時(shí)空信息演化規(guī)律及其對(duì)氣候變化的響應(yīng)研究：以寧夏為例”(41201438)； 寧夏自然科學(xué)資助項(xiàng)目(NZ14206)； 寧夏高等學(xué)?？茖W(xué)研究重點(diǎn)項(xiàng)目(NGY2014001)； 寧夏大學(xué)研究生創(chuàng)新項(xiàng)目(GIP2015041)

候靜(1989—)，女(漢族)，山西省太原市人，碩士研究生，主要從事3S技術(shù)及其應(yīng)用。 E-mail:weinigod1989@163.com。

張學(xué)儉(1965—)，男(漢族)，寧夏回族自治區(qū)鹽池縣人，碩士，研究員，從事3 S技術(shù)及其應(yīng)用等研究。 E-mail:nxjohn@163.com。

B

1000-288X(2015)05-0127-06

Q149， P94