基于NDVI指數的渾善達克沙地植被時空變化規律分析

朱芳瑩，鹿化煜，弋雙文，王曉勇

(1.南京大學 地理與海洋科學學院 地貌過程與環境實驗室，江蘇 南京 210093；2.中國人民解放軍61243部隊，新疆 烏魯木齊 830000)

基于NDVI指數的渾善達克沙地植被時空變化規律分析

朱芳瑩1,2，鹿化煜1，弋雙文1，王曉勇1

(1.南京大學 地理與海洋科學學院 地貌過程與環境實驗室，江蘇 南京 210093；2.中國人民解放軍61243部隊，新疆 烏魯木齊 830000)

歸一化植被指數(NDVI)；時空變化規律；植被覆蓋度；渾善達克沙地

為從宏觀上認識沙地綠色植被在年、月時間尺度的變化規律，為防沙治沙工作提供理論依據，分析了渾善達克沙地1982—2006年共25年的植被指數數據。研究認為，渾善達克沙地植被生長期為每年的5月下旬至10月，其中7月下旬和整個8月綠色植被生長量達到峰值；1982—2006年的25年間，渾善達克沙地綠色植被覆蓋度以每年0.1百分點的速率緩慢增長，其中沙地西部治沙成果顯著，但沙地中東部，原本的中、高植被覆蓋區，植被的衰退和顯著衰退面積竟然高達總面積的27%。同樣的防沙治沙工作，為何在渾善達克沙地的東、西部出現了截然相反的局面，這需要引起我們的反思和高度重視。

我國北方四大沙地(呼倫貝爾沙地、科爾沁沙地、渾善達克沙地、毛烏素沙地)總面積超12萬km2，是我國農業生產條件最為嚴酷、生態環境最為脆弱、人地矛盾和人水矛盾最為突出、土沙轉換最為頻繁、防治難度最大的地區，而綠色植被，正是這片廣袤土地上最有力的生力軍。研究綠色植被歷年的生長變化特點，有助于我們從宏觀上認識沙地植被在年、月時間尺度上的生長變化規律,有利于我們尋找和發現在防沙治沙過程中存在的問題，為防沙治沙工作提供理論依據。

1 研究區概況

渾善達克沙地位于內蒙古中部錫林郭勒草原南端，總面積約3.1萬km2[1]，是離北京最近的沙源。沙地內部，風沙土上生長的植被主要是草原沙地錦雞兒、柳、蒿灌叢等，河流谷地為禾草、雜草類草甸；沙地西部棕鈣土上是以戈壁針茅為主的草原，在栗鈣土上廣泛分布著克氏針茅草原；沙地的南部零星分布著些許糧田，以耕種耐寒耐旱的大豆、玉米為主；沙地的東部，靠近大興安嶺的區域，植被茂密且種類繁多，主要是以虎榛子、繡線菊、線葉菊灌叢和羊草、貝加爾針茅等植物為主的草原；在各河流流經的區域，或內陸河的發源地、分叉處，有少量的落葉小葉林生長。由于沙地孔隙充分吸收降水和攔截源于南部山地的地表水和地下潛流的補給，因而潛水儲量豐富，沙地中散布著眾多大小不等的沼澤洼地，形成了草原與沙地、沙地與濕地、濕地與湖泊、農田等多種生態系統的綜合交錯區域。

2 數據來源及處理

2.1 數據來源

GIMMS NOAA/AVHRRNDVI全球植被指數變化數據集由全球總量監測和模擬研究組收集處理。該數據集涵蓋1981—2006年共26年的月值植被指數數據，采用ALBERS投影，時間分辨率是15 d，空間分辨率是8 km。由于數據集中1981年的數據不全，故只分析1982—2006年的數據。GIMMS數據集的處理過程包括校正遙感器靈敏度隨時間的變化、太陽高度角、衛星軌道的漂移等影響，消除EL Chichon和Pinatubo火山爆發導致的大氣氣溶膠的影響。GIMMS數據集具有時間序列長、覆蓋范圍廣、時空可比、較強的植被動態變化表征能力等特點[2-3]，已被廣泛應用于全球及區域尺度的植被動態研究之中。

2.2NDVI數據處理

2.2.1 最大值合成法(MVC)

對于NDVI的處理方法，采用最大值合成法MVC(Maximum Value Composites)，用來獲取全年、全月NDVI最大值，該方法可以進一步消除云、大氣、太陽高度角等的部分干擾[4-7]。求取年最大NDVI值的計算公式為

NDVIi=MAX(MNDVI1，MNDVI2，…，MNDVI12)

(1)

式中：i為年序號；MNDVI1，MNDVI2，…，MNDVI12分別為當年的1，2，…，12月的NDVI值。

2.2.2 植被覆蓋度(fg)的計算及分級

根據GUTMAN[8]提出的亞像元模型計算植被覆蓋度，該方法將像元分為均一像元和混合像元兩種。

對于均一像元，也就是植被完全覆蓋的像元，其植被覆蓋度fg=1，像元的NDVI值主要受到植被層狀結構的影響，表達式為

NDVIg=NDVI∝-(NDVI∝-NDVI0)e-kLg

(2)

式中：NDVI∝為Lg無窮大時的NDVI值，在研究中取所有年份中最大的NDVI值為NDVI∝；NDVI0為裸土的NDVI值，在此處設為0.1[9]；k為消光系數，通常認為草地的消光系數取值范圍是0.8

對于混合像元，一個像元是由植被覆蓋部分和非植被覆蓋部分兩方面構成[12]，該像元的NDVI讀數(從數據中直接讀取的值)為植被覆蓋部分NDVIg值和非植被覆蓋部分NDVI0值的加權平均值。此時，可求得該混合像元的植被覆蓋度(fg)為

fg=(NDVI-NDVI0)/(NDVIg-NDVI0)

(3)

式中：NDVIg的計算方法見式(2)，NDVI值可以是月值或年值。

在研究中，將渾善達克沙地分為極低植被覆蓋度(fg≤20%)、低植被覆蓋度(20%80%)共5級[11]。

3 綠色植被覆蓋狀況及生長期分析

使用GIMMSNDVI數據，計算各年份年最大NDVI值，求取多年NDVI的平均值，獲得1982—2006年渾善達克沙地多年平均植被覆蓋度圖。經計算，渾善達克沙地多年平均植被覆蓋度為43%，空間分布上表現為自西部到東部植被覆蓋度逐步增大，呈逐步過渡態勢，如圖1所示。

圖1 渾善達克沙地植被覆蓋度空間分布

研究認為，植物越冬后，由黃色變為綠色并恢復生長，即為返青，一個地區半數植被返青的時間點被認為是返青期；同理，當一個地區60%的植被枯黃后，認為枯黃期隨之到來。返青期與枯黃期之間的時段，稱為綠色植被生長期。所以，當植被覆蓋度大于等于43%×50%=21.5%時，即NDVI月≥0.25時，我們認為植被的返青期到來；當植被覆蓋度小于43%×(1-60%)=17.2%時，即NDVI月<0.22時，植被進入枯黃期。通過累年各月NDVI變化圖可知，渾善達克沙地在5月下旬迎來植被的返青，7月下旬和整個8月綠色植被達到峰值，10月下旬進入枯黃期，沙地植被生長期為每年的5月下旬至10月，如圖2所示。接下來，通過fg的分級，逐月分析渾善達克沙地植被覆蓋度的空間演變，如圖3所示。

圖2 渾善達克沙地累年各月NDVI變化

注：圖中a表示上半月、b表示下半月。

圖3 渾善達克沙地植被生長期覆蓋度空間變化

5月，沙地超過90%的土地綠色植被覆蓋度低于20%，圖上幾乎呈現全黑狀態，僅沙地東部不到10%的土地綠色植被覆蓋度高于20%。

6月，沙地近一半的土地開始轉綠，其中灤河周邊約2.5%的土地綠色植被覆蓋度率先突破40%，綠色植被覆蓋度有較快的提升。

7月，植被快速返青，到7月末，綠色植被覆蓋度超過60%的土地面積已占到總面積的30%。

8月，植被長勢達到巔峰狀態，但沙地西部的沙丘鏈區域，在植被長勢最好的8月仍然地表裸露，是常年不長植被的裸露沙地的典型代表，這片區域約占沙地總面積的5%。

9月開始，植被生長從巔峰態勢進入回落階段，沙地綠色植被覆蓋度逐步下降，至9月末，沙地東部綠色植被幾乎全部開始枯黃，綠色植被覆蓋度逐漸向低于40%的方向轉變。

10月，沙地逐步進入枯黃期，至10月末，沙地的整體綠色植被覆蓋度不足14%，沙地完全進入植被枯黃期。

4 1982—2006年沙地綠色植被的時空變化

從渾善達克沙地植被覆蓋度年際變化曲線看(圖4)，沙地以中、低植被覆蓋度(20%

圖4 渾善達克沙地植被覆蓋度年際變化曲線

分析年植被覆蓋度與不同覆蓋度區域面積比例的相關性(見表1)發現，渾善達克沙地年植被覆蓋度的變化與極低、低覆蓋度區的變化在0.01水平上呈顯著負相關；與中、高、極高覆蓋度區的變化在0.01水平上呈顯著正相關。據此認為，渾善達克沙地25年來植被覆蓋度以波動緩慢提升為主，主要表現為極低、低植被覆蓋度區向中、高植被覆蓋度區的轉變。

表1 渾善達克沙地年植被覆蓋度與不同覆蓋度區域面積比例相關性分析

注：標“**”表示在0.01 水平(雙側)上顯著相關，標“*”表示在0.05水平(雙側)上顯著相關。

選取NDVI值最小的1987—1989年和NDVI值最大的2002—2004年的渾善達克沙地年NDVI值，通過柵格計算器，分析25年來沙地植被覆蓋的空間變化狀況。在分析中，將NDVI變化值大于0.1的區域劃分為植被顯著變化區，包括顯著衰退和顯著增長兩部分，NDVI變化值在-0.1～0.1之間的區域，以0為界限，分為植被衰退區和植被增長區。

由分析可知，渾善達克沙地西部是植被增長和顯著增長的區域，此處位于蘇尼特右旗境內，植被覆蓋度極低，風沙活動劇烈，是近年來防沙治沙的關鍵區域，且治沙效果顯著；相反，沙地中、東部，原本的中、高植被覆蓋區(圖1)，卻出現綠色植被的衰退和顯著衰退，沙地植被衰退和顯著衰退的面積占全區面積的27%(圖5)。同樣的防沙治沙工作，為何在渾善達克沙地的東、西部出現了如此截然相反的局面，這需要引起我們的反思和高度重視。

圖5 渾善達克沙地植被覆蓋度空間變化

5 結 論

(1)渾善達克沙地多年平均植被覆蓋度為43%，空間分布上表現為：沙地自西部到東部，植被覆蓋度逐步增大，呈逐步過渡態勢。

(2)渾善達克沙地植被生長期為每年的5月下旬至10月，其中7月下旬和整個8月綠色植被達到峰值。

(3)25年來，沙地年均植被覆蓋度以每年0.1百分點的速率緩慢增長。

(4)25年來，沙地西部是防沙治沙的關鍵區域，也是綠色植被增長和顯著增長的區域；相反，沙地中、東部(原中、高植被覆蓋區)出現植被的衰退和顯著衰退，沙地植被衰退面積占全區面積的27%。

(致謝：感謝吳霜葉教授、曾琳先生等的幫助！)

[1] 朱芳瑩,鹿化煜,張文超,等.中國北方沙漠(沙地)地理底圖繪制與地表動態過程分析——基于3S技術的新結果[J].第四紀研究,2013,33(2):197-205.

[2] FENSHOLT R,RASMUSSEN K,NIELSEN T T,et al.Evalu-ation of earth observation based long term vegetation trends—Intercomparing NDVI time series trend analysis consistency of Sahel from AVHRR GIMMS,Terra MODIS and SPOT VGT data[J].Remote Sensing of Environment,2009,113(9):1886-1898.

[3] BECK P S A,GOETZ S J.Satellite observations of high nor-thern latitude vegetation productivity changes between 1982 and 2008:ecological variability and regional differences[J].Environmental Research Letters,2011,6(4):45501-45510.

[4] CHEN P Y,SRINIVASAN R,FEDOSEJEVS G,et al.Evalu-ating different NDVI composite techniques using NOAA-14 AVHRR data[J].International Journal of Remote Sensing,2003,24(17):3403-3412.

[5] DU PLESSIS W P.Linear regression relationships betweenNDVI,vegetation and rainfall in Etosha National Park,Namibia.[J].Journal of Arid Environments,1999,42(4):235-260.

[6] 張戈麗,徐興良,周才平,等.近30年來呼倫貝爾地區草地植被變化對氣候變化的響應[J].地理學報,2011,66(1):47-58.

[7] 張月叢,趙志強,李雙成,等.基于SPOT NDVI的華北北部地表植被覆蓋變化趨勢[J].地理研究,2008,27(4):745-754.

[8] GUTMAN G,IGNATOV A.The derivation of the green vegetation fraction from NOAA/AVHRR data for use in numerical weather prediction models[J].International Journal of Remote Sensing,1998,19(8):1533-1543.

[9] 方精云,樸世龍,賀金生,等.近20年來中國植被活動在增強[J].中國科學:生命科學,2003,33(6):554-565.

[10] 陳晉,陳云浩,何春陽,等.基于土地覆蓋分類的植被覆蓋率估算亞像元模型與應用[J].遙感學報,2001,5(6):416-422.

[11] 陳效逑,王恒.1982—2003年內蒙古植被帶和植被覆蓋度的時空變化[J].地理學報,2009,64(1):84-94.

[12] 劉廣峰,吳波,范文義,等.基于像元二分模型的沙漠化地區植被覆蓋度提取——以毛烏素沙地為例[J].水土保持研究,2007,14(2):268-271.

(責任編輯 徐素霞)

庫姆塔格沙漠科學綜合考察項目(2012FY111700)；戈壁生態系統地質本底調查；中國沙漠沙地志編寫項目

S157

A

1000-0941(2017)03-0031-04

朱芳瑩(1989—)，女，江蘇啟東市人，助理工程師，碩士，主要從事沙漠生態調查、生態修復等研究。

2016-04-19