胡馨月，孫 倩，古麗古麗·滅迪力哈孜

（新疆農(nóng)業(yè)大學林學與風景園林學院，烏魯木齊 830052）

0 引言

植被是連接自然界水、大氣、土壤內部循環(huán)的重要途徑和橋梁，它也是影響全球環(huán)境及其氣候變化的標志，更是陸地自然生態(tài)系統(tǒng)的重要因素。歸一化植被指數(shù)可以反映區(qū)域植被的生長和其營養(yǎng)資源信息的重要表征性參數(shù)，它還可以用來衡量研究區(qū)域植被的覆蓋度和其生長情況的有效程度度量參數(shù)和重要表征性指數(shù)，廣泛應用于植被的生長和其他植被涵蓋度進行檢測[1]。有學者結合新疆地區(qū)1982—2013年NDVI數(shù)據(jù)，發(fā)現(xiàn)新疆植被覆蓋空間分布南北差異明顯，西北優(yōu)于東南[2]。在區(qū)域和小尺度上，學者研究了不同尺度上對不同地理單元的NDVI植被時空變化和氣候因子[3-9]。歸一化植被指數(shù)是目前最廣泛應用的植被指數(shù)，也是監(jiān)測植被覆蓋狀況和其生長的最佳植被指數(shù)[10]。很多學者對NDVI數(shù)據(jù)做出長期研究，得出了植被覆蓋度變化和影響因子的相關性。近年來，通過NDVI做了大量的科學研究，并且取得了諸多科學研究成果，在全球區(qū)域等宏觀和大尺度的范圍內，對于土地利用覆蓋、植被的動態(tài)變化、植被的分類、作物生長勢、氣候的影響以及監(jiān)測等諸多領域進行了研究，NDVI已經(jīng)具有非常良好的實際適用性[11-13]。

NDVI是最常見且對土壤背景變化很敏感的植被指數(shù)，增強了對植被的響應，它可以消除儀器定標、地形、太陽角，云陰影和大氣條件等對半徑變化的影響，是當前40多種植被指數(shù)中應用最廣的指標之一[14-15]。鑒于此，本研究利用MODIS數(shù)據(jù)產(chǎn)品，通過數(shù)據(jù)預處理獲布爾津縣2010—2014年逐月的MODIS-NDVI數(shù)據(jù)，研究布爾津縣歸一化植被指數(shù)時空演變，分布特征，趨勢分析。本研究結論將定量的呈現(xiàn)布爾津縣植被覆蓋的時空演變規(guī)律，能夠為生態(tài)環(huán)境的區(qū)域保護提供一些科學參考。

1 研究地區(qū)與研究方法

1.1 研究區(qū)域概況

布爾津縣(86°25′—88°06′E，47°22′—49°11′N)屬于新疆阿勒泰地區(qū)，位于阿爾泰山西南麓，與三國接壤是中國重要的邊境縣。布爾津縣呈葫蘆形，南北長200 km，東西寬49～82 km，總面積10345 km2，國界線長218 km。布爾津縣屬大陸性北溫帶寒涼氣候，夏季干燥炎熱，冬季寒冷，降水量小，蒸發(fā)量大，晝夜溫差大，光照充足[16]。

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

遙感數(shù)據(jù)來源于地理空間數(shù)據(jù)云所提供的空間分辨率為1 km的中國2010—2014年5年MODIS-NDVI逐月產(chǎn)品和研究區(qū)矢量邊界。數(shù)據(jù)進行大氣校正、輻射校正，幾何校正，選取布爾津縣范圍，對月度植被指數(shù)空間分布數(shù)據(jù)進行掩膜提取，得到2010—2014年NDVI數(shù)據(jù)，并得到了60幅影像圖。利用ArcGIS軟件，以研究區(qū)影像圖將生成的原始DN數(shù)據(jù)集轉換為標準的逐旬NDVI值，原始NDVI值為擴展NDVI值，采用取值0～10000范圍，DN值和實際值NDVI為DN=NDVI*10000，利用公式得到NDVI實際值。

1.3 研究方法

（1）歸一化植被指數(shù)NDVI

歸一化植被指數(shù)NDVI是根據(jù)通過測量植被強烈反射與植被吸收在近紅外和紅光電磁波段的吸附率之間的值作為數(shù)據(jù)，通過對電磁輻射的特征值進行數(shù)字化分析而得到的歸一化植被覆蓋率的表征指數(shù)，計算如公式(1)。

式中：NIR為700～1000 nm近紅外波段的波譜特征；RED為650 nm紅光波段的波譜特征。

（2）趨勢分析方法

通過一元線性回歸趨勢分析模擬出各個柵格的變動趨勢，可以反映出研究區(qū)植被NDVI空間變化特征[17]，計算如公式(2)。

式中：θslope為NDVI的變化趨勢；n為研究區(qū)的年數(shù)；θslope＞0時，表示NDVI在研究期間內的變化趨勢增加；θslope＜0時，表示NDVI在研究期間內的變化趨勢減小。

2 結果分析

2.1 NDVI的時序變化

布爾津縣屬溫帶干旱氣候區(qū)，3—5月為春季、6—8月為夏季、9—11月為秋季、對NDVI數(shù)據(jù)均值化處理后得到2010—2014年年均各季節(jié)NDVI數(shù)據(jù)集，取均值后制作的NDVI時序變化曲線，如圖1～2。

圖1 2010-2014年布爾津縣三季NDVI趨勢變化圖

圖2 2010-2014年布爾津縣四季NDVI指數(shù)動態(tài)特征圖

由圖1所示，布爾津縣2010—2014年的NDVI平均值在春夏秋三季整體呈現(xiàn)下降的趨勢，2010—2014年的春季的NDVI值波動幅度較小，夏季的NDVI值呈下降趨勢，秋季的NDVI值波動較小。其中2010—2014年的夏季NDVI值都達到最大值，每一年春季NDVI都為最小值。

由圖2所示，NDVI在夏季值最高，2010年最高達到0.53。冬季NDVI值最低，2010年最低值0.01。春秋兩季NDVI值波動幅度緩慢，除2010年春季NDVI值為0.12明顯偏低外，其他年份春季NDVI值均在0.2左右。2010—2014年四季NDVI值波動幅度較平緩。便于理解和分析，人為將生長周期分為不同的階段：春季是植被生長為綠色的時期，且NDVI值呈現(xiàn)平緩增加的現(xiàn)象，NDVI值減少的原因是春季該地區(qū)降水量減少，春季植被覆蓋度不高過渡放牧所致；夏季是植被活力及作物生長最旺盛的時期，NDVI值較為平穩(wěn)；秋季植被新陳代謝減慢，植被的覆蓋率開始下降，NDVI值總體呈下降趨勢；冬季NDVI值沒有明顯的浮動，趨近于0，大部分植被在冬季不能存活。

2.2 NDVI的空間分布特征

本研究使用ArcGIS 10.5對2010—2014年布爾津縣遙感NDVI數(shù)據(jù)進行數(shù)據(jù)提取，為了可以更加直觀的反映布爾津縣5年的NDVI空間分布特征，以5年所監(jiān)測的各月NDVI數(shù)據(jù)為基礎，分析布爾津植被在空間上分布規(guī)律（見圖3）。冬季NDVI數(shù)值相對來說比較穩(wěn)定，變化波動非常微弱，相對而言，春夏秋更加具有研究價值，所以篩選了時空分布上具有明顯特征性的春夏秋三個季節(jié)進行不同年度同季節(jié)的對比。

圖3 2010—2014年布爾津縣NDVI空間分布圖

NDVI值越大，說明植被變化情況劇烈，NDVI值越小，說明植被變化呈下降趨勢。植被指數(shù)較高區(qū)域主要分布地區(qū)中部，植被指數(shù)較低區(qū)域主要分布在西南部，隨著時間的變化推移，NDVI空間布局上發(fā)生了變化，春季分布由北向南逐漸增加，夏季分布沒有明顯變化，秋季分布由北向南先逐年減少。

2010—2014年的春季、夏季、秋季NDVI空間分布數(shù)值，2010年春夏秋對比圖中，夏季的NDVI值達到最高，為0.94；秋季總體上來看，它的數(shù)值與夏季非常接近，最高值為0.92。從圖3可以看出在夏季NDVI高值的區(qū)域面積明顯要大于秋季NDVI數(shù)值高的面積。2011年春夏秋對比圖中，夏季NDVI最高值為0.97；秋季的數(shù)值與春季相接近，最高值為0.77，可以看出夏季NDVI高值的區(qū)域面積分布最廣，春季NDVI高數(shù)值分布區(qū)域面積高于秋季NDVI數(shù)值高值的區(qū)域面積。2012年三季分布圖中，夏季NDVI最高值為0.94；春季NDVI值總體上來說與夏季接近，最高值為0.74，從圖中可以看出夏季NDVI高值分布區(qū)域面積最大，春季高值分布區(qū)域要大于秋季NDVI高值分布區(qū)域面積。2013年的分布圖中，夏季NDVI最高值為0.93，秋季的NDVI總體上來說與夏季接近，最高值為0.89，夏季NDVI高值分布區(qū)域面積大于秋季NDVI高值分布區(qū)域面積。2014年分布圖中，夏季NDVI最高值為0.93，秋季NDVI高值總體上與夏季NDVI高值接近，秋季NDVI低閾值分布范圍最廣。

以不同年份同季節(jié)中，可以看出在春季2011年和2012年的NDVI高值分布區(qū)域較少，其余年份趨于相近，逐年布爾津縣中部植被分布增加，大部分地區(qū)都呈現(xiàn)出植被分布較少，西南地區(qū)邊緣植被分布較少。夏季在2010—2014年中NDVI空間分布變化較小，北部植被分布有所增加。秋季在2011年植被分布最廣，2014年植被分布最少，其余年份植被分布趨于相近，北部植被分布逐漸減少。

2.3 NDVI變化趨勢

變化趨勢SLOPE的計算，能更好的表征出隨著時間變化，長時間序列歸一化植被指數(shù)的動態(tài)變化情況。植被NDVI變化趨勢呈現(xiàn)出夏季高于秋季，秋季高于春季，冬季較穩(wěn)定的特征，數(shù)據(jù)如圖4所示。

圖4 2010—2014年布爾津縣四季NDVI變化趨勢圖

5年間，春季NDVI值最小為-1.06，為負值，表示布爾津縣植被NDVI呈現(xiàn)出了逐漸減小的趨勢，由南向北NDVI呈現(xiàn)出逐步增加的分布特征，到北部邊緣逐漸減少。夏季NDVI值最小為0.01，NDVI值逐漸增加，植被分布廣泛，其夏季NDVI平均變化趨勢較小，但呈現(xiàn)逐漸增加的變化趨勢。夏季冰川融化且雨水充足，地表徑流豐富，人為的種植植物較多，NDVI值越高。

秋季NDVI值最小為0，為正值，呈現(xiàn)出植被指數(shù)呈現(xiàn)出了逐步增加的趨勢。秋季平均降水量降低，農(nóng)作物進入收獲期，耐旱適鹽的植被可以在秋季繼續(xù)生長和存活。冬季NDVI值最小為-0.01，北部NDVI值降低趨勢顯著，其植被變化趨勢減小。冬季11—12月，大部分耐寒耐旱植被仍保持生長，土壤的溫度與含水率對植物的生長有較明顯的影響。這一階段土壤氣候溫度及相對濕度都比較低，在土壤條件的影響下部分植被緩慢地衰敗甚至死亡，植被指數(shù)下降。

3 討論

2010—2014年布爾津縣四季NDVI的變化中，冬季NDVI均值最低，夏季均值最高，造成這一結果的原因是由于冬季溫度較低、降水量較少，部分植被進入落葉期，植被處于生長季末期；而夏季降水量充沛，溫度較高，植被光合作用活躍，處于生長季最旺盛時期。本研究通過分析2010—2014年三季NDVI時序變化，發(fā)現(xiàn)NDVI呈現(xiàn)隨著時間的變化，NDVI的平均值也增加，特別是夏季和秋季NDVI增長較為明顯。

NDVI空間分布表明，研究區(qū)在夏季NDVI平均值達到了最高值，主要分布西南部和東北部逐步減少，可能是由于夏季溫度較高，水分蒸發(fā)速度快，而植物需水量大，不能快速供給植物生長所必需的水分，從而導致部分植物受到干旱脅迫，從而導致植被變化不均。春秋季NDVI均值較低，甚至出現(xiàn)負值，主要分布北部方向植被減少，可能是由于春季氣溫逐漸升高使水汽增多，出現(xiàn)多云天氣。秋季NDVI值與春季NDVI值相接近，秋季植被分布相對沒有明顯變化。

新疆是溫帶大陸性氣候，在全球氣候變化背景下，新疆NDVI時空變化引起了眾多學者的密切關注。相關研究表明[18-20]，近30年來新疆NDVI呈明顯增加趨勢，尤其是在夏秋兩季增加趨勢最為明顯。本研究結果表明，布爾津縣2010—2014年NDVI也表現(xiàn)出與整個新疆相同的變化規(guī)律，布爾津縣2010—2014年NDVI變化趨勢具有更明顯的規(guī)律性，即總體NDVI數(shù)值變化越大，呈現(xiàn)的增加趨勢更為顯著，NDVI數(shù)值變化越小，增加趨勢不明顯。NDVI變化趨勢較小，代表該區(qū)域內的植被分布較穩(wěn)定，植被生長狀況良好，表明該區(qū)域的植被內部環(huán)境較穩(wěn)定不易受到外界氣候或環(huán)境的影響。變化趨勢顯著增加則表示該地區(qū)的植被呈下降趨勢，植被的覆蓋度降低[21-23]。NDVI下降的主要原因有，城市化進程加快，旅游業(yè)發(fā)展，工業(yè)生產(chǎn)加速以及大規(guī)模的城市建設等人為活動引起植被減少[24-26]。

本研究分析了布爾津縣5年的歸一化植被指數(shù)時序變化研究和空間分布特征，但還存在一定局限性，由于研究利用分布率問題可能無法顯示特殊區(qū)域的變化特征，在今后的分析中將通過更高分布率數(shù)據(jù)對特殊區(qū)域進行分析，在變化因素的分析中應考慮人類活動對植被覆蓋度造成的影響。

4 結論

通過對布爾津縣2010—2014年基于時間序列的MODIS-NDVI數(shù)據(jù)處理，從時間和空間角度，采用季節(jié)變化和變化趨勢分析的方法，對NDVI數(shù)據(jù)進行定量分析，得出以下結論：

（1）2010—2014年布爾津縣NDVI從時序變化角度而言，NDVI存在明顯的季節(jié)性差異，夏季NDVI值相對較高，均在0.05～0.6之間波動。春季和秋季年際的NDVI值波動幅度較小，均分別在0.12～0.21和0.29～0.34之間波動。不同季節(jié)NDVI年際變化的總體大小排序為：夏季＞秋季＞春季＞冬季。

（2）夏季NDVI分布的變化趨勢不明顯，夏季植被最為茂盛，空間分布無明顯變化，其中2011年植被生長態(tài)勢最佳；春季NDVI數(shù)值先減少后增加，空間分布上中部和西南部植被分布有所增加，北部植物相對減少，秋季NDVI數(shù)值呈現(xiàn)了先降低后增加的變化趨勢，但變化趨勢相對較小，說明秋季植被生長態(tài)勢變化不劇烈，尤其在研究區(qū)東北方向，植物分布相對減少。

（3）布爾津縣四季5年間的NDVI變化趨勢，春季NDVI呈現(xiàn)了逐步降低的變化趨勢，夏季和秋季NDVI變化趨勢均呈現(xiàn)出了增加趨勢，冬季NDVI變化趨勢微弱。