喻素芳+佘光輝+羅葉紅+潘婷+李麗娟

摘要：利用月合成MODIS-NDVI時(shí)間序列數(shù)據(jù)、氣象數(shù)據(jù)、DEM等數(shù)據(jù)，分析研究了2001-2010年廣西省10年間植被覆蓋變化的特征。結(jié)果表明，廣西植被覆蓋率較高且總體呈上升趨勢(shì)，從區(qū)域角度分析，桂南地區(qū)增加趨勢(shì)最顯著，其次是桂中地區(qū);從海拔梯度角度分析，海拔200 m以下的區(qū)域植被覆蓋增加趨勢(shì)最顯著，但此區(qū)域的植被覆蓋率較低，NDVI最高分布在海拔400～800 m的區(qū)域。植被覆蓋變化受氣溫、降水及人類活動(dòng)影響較大，植被覆蓋梯度變化除受這些因素影響外，也與植被垂直地理分布特征相關(guān)。

關(guān)鍵詞：MODIS-NDVI;植被覆蓋;氣溫;顯著

中圖分類號(hào)：S812 ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A ? ? ? ?文章編號(hào)：0439-8114（2015）02-0321-05

DOI：10.14088/j.cnki.issn0439-8114.2015.02.017

Analysis of Characteristics of Vegetation Cover Change in Guangxi Province Based on MODIS-NDVI

YU Su-fang1，2，SHE Guang-hui1，LUO Ye-hong2，PAN Ting2， LI Li-juan2

（1. School of Forest Resources and Environment ，Nanjing Forestry University ，Nanjing 210037， China;

2. School of Forestry，Guangxi University ，Nanning 530004， China）

Abstract： Base on the monthly time series MODIS-NDVI， meteorologic and digital elevation model（DEM） data from 2001 to 2010， analysed the change characters of surface vegetation coverage of Guangxi Province during this 10 years. The results showed that the vegetation coverage of Guangxi was at a high level with a growing tendency. The growing tendency was the most remarkable in southern Guangxi and followed by the middle Guangxi area. Vegetation coverage change with the vary classified elevation levels. The 200 m below level hold the most obvious growing trend while its NDVI was lower. NDVI was highest at altitude of 400～800 m range. The temperature， rainfall and human activities could influence the vegetation coverage. Beside all these three factors，the vegetation coverage was also associated with the vertical characteristics of geographical distribution of vegetation.

Key words： MODIS-NDVI; vegetation cover; temperature; obvious

地表植被的空間分布是植被對(duì)氣候長(zhǎng)期適應(yīng)的結(jié)果，任何地區(qū)的植被均能反映出該地區(qū)的氣候類型[1]，并且植被覆蓋狀況能直接反映該地區(qū)生態(tài)環(huán)境狀況，對(duì)全球環(huán)境變化具有指示器作用[2]。研究某一區(qū)域植被變化對(duì)研究該區(qū)域生態(tài)系統(tǒng)具有重要意義[3]。

植被指數(shù)（Vegetation Index）是指從多光譜遙感數(shù)據(jù)中提取的有關(guān)地球表面植被狀況的定量數(shù)值，用以表征地表植被覆蓋、生長(zhǎng)狀況和生物量等，常用的植被指數(shù)有：RVI、GVI、DDVI、PVI、EVI及NDVI等[4，5]。歸一化植被指數(shù)（NDVI），又稱標(biāo)準(zhǔn)化植被指數(shù)，是植物生長(zhǎng)狀態(tài)以及植被空間分布密度的最佳指示因子，與植被分布密度呈線性相關(guān)。同時(shí)，NDVI對(duì)土壤背景的變化較為敏感;與綠葉的葉面積系數(shù)（LAI）有很好的正相關(guān)性，可直接表征某地區(qū)的植被覆蓋狀況[3，5]。

一些學(xué)者基于不同的遙感數(shù)據(jù)源研究了不同地區(qū)地表植被的時(shí)空分布、覆被變化及其驅(qū)動(dòng)力等內(nèi)容[6-8];宋怡等[7]、崔曉臨等[8]利用SPOT、MODIS遙感數(shù)據(jù)對(duì)中國(guó)西北地區(qū)的植被覆蓋變化進(jìn)行了分析研究;李秀花等[9]、劉綠柳等[10]、李雙成等[11]、張笑鶴[12]研究了NDVI的動(dòng)態(tài)演變規(guī)律與氣候之間的相關(guān)性，結(jié)果表明，全球氣候變暖的顯著變化影響了地表植被覆蓋;崔曉臨等[13]研究得出不同海拔NDVI與氣溫具有不同的相關(guān)性;何勇等[14]發(fā)現(xiàn)在眾多氣候因子中，氣溫及降水對(duì)植被NDVI具有最重要和最直接的影響;潘薔等[15]通過(guò)對(duì)北京植被指數(shù)變化及影響因素分析得出，除了氣候變化影響NDVI外，人類活動(dòng)包括土地利用類型的改變、農(nóng)業(yè)生產(chǎn)水平的提高以及植被建設(shè)的管理對(duì)植被覆蓋均具有不同程度的影響。目前，國(guó)內(nèi)對(duì)植被覆蓋變化的研究主要集中在長(zhǎng)江以北地區(qū)[7-15]，華南地區(qū)極少，本研究基于2000-2010年MODIS NDVI時(shí)間序列數(shù)據(jù)，以地處亞熱帶地區(qū)的華南西部的廣西省為研究區(qū)域，分析該區(qū)域10年間植被覆蓋變化特征，并試圖從氣候變化、植被分布等方面分析和探討植被覆蓋變化的原因，以期為區(qū)域植被恢復(fù)和管理提供理論支持。endprint

1 ?研究區(qū)概況

廣西壯族自治區(qū)地處中國(guó)南疆，位于東經(jīng)104°26′-112°04′，北緯20°54′-26°24′之間，北回歸線橫貫全區(qū)中部，南瀕熱帶海洋，北接南嶺山地，西延云貴高原，屬云貴高原向東南沿海丘陵過(guò)渡地帶，地形特點(diǎn)為四周高中部低，形似盆地，且山地多、平原少（圖1）。廣西地處中、南亞熱帶季風(fēng)氣候區(qū)，在太陽(yáng)輻射、大氣環(huán)流和地理環(huán)境的共同作用下，形成了氣候溫暖、熱量豐富（各地年平均氣溫在16.5～23.1 ℃之間），降水豐沛、干濕分明（各地年均降水量為1 080～2 760 mm，大部分地區(qū)為1 300～2 000 mm），日照適中、冬少夏多（各地年均日照時(shí)數(shù)為1 169～2 219 h），災(zāi)害頻繁、旱澇突出，沿海、山地風(fēng)能資源豐富的氣候特點(diǎn)。

2 ?材料與方法

2.1 ?試驗(yàn)材料

中分辨率成像光譜儀（MODerate-resolution Imaging Spectroradiometer-MODIS）是Terra和Aqua衛(wèi)星上搭載的主要傳感器之一，兩顆衛(wèi)星相互配合每1～2 d可重復(fù)觀測(cè)整個(gè)地球表面，本研究使用了從中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站（http：//www.gscloud.cn）下載的分辨率為1 km的MODIS月合成NDVI數(shù)據(jù)（TERRA星），時(shí)間序列數(shù)據(jù)由2001年1月至2010年12月共120期，該產(chǎn)品是由分辨率為250 m的MODIS NDVI每天的數(shù)據(jù)產(chǎn)品計(jì)算得到，計(jì)算方法為取月內(nèi)最大值。其他材料有廣西壯族自治區(qū)行政區(qū)劃圖（圖1）、從廣西氣象局獲取的2001至2010年平均氣溫和年平均降水量、廣西2001至2011年統(tǒng)計(jì)年鑒及從中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)（http：/datamirror.csdb.cn）下載的分辨率為90 m的DEM數(shù)據(jù)。

2.2 ?試驗(yàn)方法

利用軟件ERDAS 9.1和ARCGIS 10.0對(duì)獲取的120期月平均NDVI數(shù)據(jù)進(jìn)行解譯和分析，發(fā)現(xiàn)在這10年中，每年上半年的數(shù)據(jù)（1～6月份）均有1～2個(gè)月的數(shù)據(jù)不同程度的受到云層的影響，為了使各年份間的數(shù)據(jù)具有更可靠的比較性，在計(jì)算每個(gè)月的10年平均值時(shí)，剔除了受到云層影響的數(shù)據(jù)，用無(wú)云層覆蓋的數(shù)據(jù)計(jì)算取平均值;年平均值采用各年下半年（7～12月份）的數(shù)據(jù)進(jìn)行計(jì)算取平均值，以其作為年平均值。

為分析過(guò)去幾年研究區(qū)域每個(gè)像元的NDVI變化趨勢(shì)，本研究利用Stow等[16]研究的以時(shí)間t為自變量用每年的平均NDVI和時(shí)間序列建立的一元線性回歸方程的分析方法來(lái)進(jìn)行分析，公式如下：

Slope=■ （1）

式（1）中：Slope表示線性回歸方程的斜率，n表示累計(jì)總監(jiān)測(cè)年數(shù)，i表示第1～n年，■■表示第i年的平均值，這條直線不是簡(jiǎn)單的描述第一年與最后一年的植被指數(shù)的關(guān)系，而是這10年總的變化趨勢(shì)，如果Slope=0，說(shuō)明NDVI基本保持不變;如果Slope<0，說(shuō)明NDVI呈遞減趨勢(shì);反之，NDVI則呈遞增趨勢(shì)。

3 ?結(jié)果與分析

3.1 ?平均NDVI值及年際變化

2001-2010年10年間月平均NDVI值變化情況如圖2所示。由圖2可知，全年中NDVI值最小值出現(xiàn)在2月份，其值為0.553，隨月份增加，氣溫回升，NDVI不斷增大，到9月份，出現(xiàn)峰值0.861，9月過(guò)后又逐漸下降，至次年的2月份。

由圖3可知，NDVI值總體呈增加的趨勢(shì)，但在這10年中，2002年的NDVI出現(xiàn)了一個(gè)小峰值，由廣西省2001-2010年年平均氣溫和年平均降水量圖（圖4）可知，2001年雨量充沛、氣溫適宜，這些綜合因素使得2002年植被長(zhǎng)勢(shì)良好;2009年NDVI出現(xiàn)了一個(gè)較低值拐點(diǎn)，為導(dǎo)致這一年NDVI下降的主要影響因素;2008年，廣西遭遇了特大雨水冰凍災(zāi)害，低溫導(dǎo)致植被、作物、莊稼均遭到極其嚴(yán)重的損失（圖4），從而NDVI值下降，說(shuō)明氣溫與NDVI值存在明顯的相關(guān)性，這與李蕙敏等[4]的研究結(jié)果相似。

3.2 ?10年間變化趨勢(shì)分析

根據(jù)地理位置將廣西全區(qū)分為五個(gè)區(qū)：桂北地區(qū)（包括桂林、柳州、河池）、桂東地區(qū)（包括賀州、梧州、貴港）、桂西北地區(qū)（包括百色）、桂中地區(qū)（包括南寧、來(lái)賓、崇左）和桂南地區(qū)（包括北海、防城港、欽州、玉林）。利用Stow等[16]研究的回歸方程對(duì)研究區(qū)域每個(gè)像元的變化趨勢(shì)進(jìn)行分析， 根據(jù)計(jì)算得到的Slope， 將其分為五個(gè)等級(jí)： 退化（<-0.000 4）、 基本不變（-0.000 4～0.000 4）、 輕微改善（0.000 4～0.004 0）、中度改善（0.004 0～0.009 0）及明顯改善（>0.009 0）。

2001-2010年的10年植被覆蓋變化趨勢(shì)圖（圖5）和面積分布圖（圖6）表明，廣西全區(qū)的植被覆蓋整體呈增長(zhǎng)趨勢(shì)，全區(qū)83%以上的區(qū)域植被覆蓋情況得到了不同程度的改善，其中中度改善和明顯改善的面積分別為43.4%和11.9%，合計(jì)占總面積的一半以上，而退化的區(qū)域僅占全區(qū)總面積的12.0%。

廣西為全國(guó)巖溶分布面積大、石漠化嚴(yán)重的?。▍^(qū)）之一，石漠化已成為廣西最大的生態(tài)問(wèn)題。長(zhǎng)期以來(lái)，林業(yè)部門(mén)積極采取各種措施，改善石漠化地區(qū)生態(tài)環(huán)境，將石漠化治理與當(dāng)?shù)靥厣a(chǎn)業(yè)發(fā)展、農(nóng)民增收和脫貧致富相結(jié)合，多年來(lái)，通過(guò)大力實(shí)施造林綠化工程、退耕還林工程、珠江流域防護(hù)林工程、森林生態(tài)效益補(bǔ)償、中央財(cái)政造林補(bǔ)貼試點(diǎn)項(xiàng)目、新造油茶示范林項(xiàng)目等，加快植樹(shù)造林和封山育林，大幅度增加森林植被，這一系列的措施使得植被的恢復(fù)取得了良好的效果，各區(qū)域植被覆蓋情況有所改善，如圖7所示。由圖7可知，從區(qū)域分布來(lái)看，桂南、桂中及桂西北地區(qū)植被恢復(fù)明顯，其中桂南地區(qū)的植被恢復(fù)程度最高，恢復(fù)面積占桂南地區(qū)總面積的93.4%，并且北海、欽州、崇左和來(lái)賓等地市的Slope值較高的像元比較集中。植被退化主要集中在桂北和桂東地區(qū)，其中桂林、賀州、貴港北部和柳州北部區(qū)域較為突出，桂北地區(qū)2008年遭遇了特大雨雪冰凍災(zāi)害，植被遭到破壞且恢復(fù)較慢，桂東地區(qū)近年經(jīng)濟(jì)發(fā)展較快，城市化進(jìn)程加速，說(shuō)明氣溫的變化、經(jīng)濟(jì)的發(fā)展對(duì)植被的恢復(fù)有較大的影響。endprint

3.3 ?不同海拔NDVI年際變化

廣西的地勢(shì)北高南低（圖1），海拔從沿海-25 m至2 140 m， 平均海拔394 m， 利用空間分辨率為90 m的DEM數(shù)據(jù)將廣西全區(qū)高程分為五個(gè)等級(jí)區(qū)域：200 m以下、200～400 m、400～800 m、800～1 200 m和1 200 m以上。利用arcgis統(tǒng)計(jì)各等級(jí)高程范圍的10年平均NDVI值，統(tǒng)計(jì)結(jié)果（表1）表明，平均NDVI值最大的區(qū)域?yàn)橹械群０胃叨龋?00～800 m）區(qū)域，最大均值為0.823 7;平均NDVI值最小的區(qū)域?yàn)楹０?200 m的區(qū)域，最小均值為0.761 3;而高山地區(qū)（海拔在1 200 m以上）的平均NDVI值僅高于海拔<200 m的區(qū)域，海拔200～1 200 m范圍平均NDVI值差異不大。利用年份作為橫坐標(biāo)，NDVI值作為縱坐標(biāo)，建立各等級(jí)高程的NDVI變化趨勢(shì)圖（圖8），由圖8可見(jiàn)，不同等級(jí)的高程區(qū)域，10年間NDVI的變化均呈上升趨勢(shì)，利用Excel軟件對(duì)各高程范圍的NDVI進(jìn)行線性回歸分析（表1），結(jié)果表明，海拔<200 m的區(qū)域范圍NDVI增加趨勢(shì)最顯著，其他依次為：200～400 m、800～1 200 m、400～800 m、1 200 m以上。

平均NDVI值海拔分布特征與廣西植被的垂直地帶性分布特征有比較密切的關(guān)系，海拔<200 m的區(qū)域，主要以農(nóng)田耕地植被為主，建設(shè)區(qū)和水域主要分布在這一海拔范圍內(nèi)，人類活動(dòng)頻繁，NDVI值低。隨著國(guó)家一系列恢復(fù)植被的措施的實(shí)施，該區(qū)域人類活動(dòng)受到限制，減弱了對(duì)植被覆蓋的影響，建城區(qū)的綠化覆蓋面積也在逐年遞增（2010年是2001年2.66倍），使得該區(qū)域的植被覆蓋增加迅速;海拔200～1 200 m的低山丘陵和中山區(qū)域，植被主要為季節(jié)性雨林、季風(fēng)常綠闊葉林和典型常綠闊葉林[17]，闊葉林郁閉度較高。另外，此區(qū)域植被的熱量、光照、水分等條件均比較充足，植物的生長(zhǎng)環(huán)境良好，植被覆蓋率較大，NDVI較高;海拔1 200 m以上的區(qū)域?yàn)槌＞G闊葉和針葉混交林，在風(fēng)力較大的山頂和山脊則形成山頂苔蘚矮林，針葉林和苔癬矮林郁閉度均相對(duì)較低，且隨著海拔的升高，氣溫降低，熱量減少，植被長(zhǎng)勢(shì)降低，NDVI值相對(duì)較低，但隨著全球氣候變暖，高海拔區(qū)域的植被覆蓋出現(xiàn)了逐年遞增的趨勢(shì)。

4 ?小結(jié)與討論

4.1 ?結(jié)論

本文采用2001-2010年MODIS NDVI時(shí)間序列數(shù)據(jù)和空間分辨率90 m的DEM數(shù)據(jù)研究廣西全區(qū)10年植被覆蓋變化特征，研究結(jié)果表明：（1）廣西土地面積2 374.08萬(wàn)hm2，研究大區(qū)域范圍植被覆蓋率的連續(xù)動(dòng)態(tài)變化，需要有足夠的空間和時(shí)間分辨率的觀測(cè)數(shù)據(jù)，MODIS數(shù)據(jù)是一個(gè)合理的選擇，其重復(fù)觀察周期短（1～2 d），空間分辨率適中（250 ～1 000 m）;（2）廣西總體植被覆蓋率較高，每年9月NDVI值最高，秋季到冬季逐漸下降，2月最低，春季到夏季逐漸升高，全年變化呈單峰曲線。2001-2010年由于受到氣溫、降水以及人類活動(dòng)等影響，年平均NDVI的變化出現(xiàn)小幅波動(dòng)，但總體呈上升趨勢(shì);（3）研究區(qū)域植被恢復(fù)效果明顯，明顯改善占11.9%、中度改善占43.4%、輕度改善占28.6%，基本不變占4.1%、退化區(qū)域所占12.0%，其中，桂南和桂中地區(qū)改善明顯的區(qū)域占比例較大;（4）植被覆蓋具有明顯梯度分布特征，植被覆蓋率從高到低依次為400～800 m、200～400 m、800～1 200 m、1 200 m以上、200 m以下， 增加趨勢(shì)從高到低依次為200 m以下、200～400 m、800～1 200 m、400～800 m、1 200 m以上。

綜上研究結(jié)果，NDVI值與氣候因子、植被分布特點(diǎn)以及地區(qū)建設(shè)發(fā)展有關(guān)，可以比較準(zhǔn)確地反映地區(qū)植被的覆蓋狀況。長(zhǎng)期進(jìn)行NDVI值觀測(cè)，可以將這些方面作為植被覆蓋變化驅(qū)動(dòng)力進(jìn)行研究，發(fā)現(xiàn)其內(nèi)在的作用機(jī)制，為區(qū)域植被恢復(fù)和管理提供理論支持。

4.2 ?建議

1）宋怡等[7]、崔曉臨等[8]采用一年中的最大NDVI值作為年平均值對(duì)研究區(qū)域的植被覆蓋進(jìn)行分析研究，而本研究利用7～12月份NDVI值求平均作為年平均值，兩種計(jì)算方法的結(jié)果對(duì)分析植被覆蓋變化特征是否存在差異，有待進(jìn)一步論證。

2）植被覆蓋的變化受諸多因素的影響，本研究?jī)H從氣溫、降水量、植被分布特征等進(jìn)行了分析，在后續(xù)的研究中可以考慮廣西土地利用類型的變化與植被覆蓋變化的相關(guān)性。

3）中國(guó)現(xiàn)今主要進(jìn)行的是城市帶動(dòng)周邊鄉(xiāng)鎮(zhèn)和農(nóng)村發(fā)展方針，由于區(qū)域發(fā)展也是植被變化的一個(gè)主要影響因子，因而未來(lái)可進(jìn)行城區(qū)內(nèi)外的對(duì)比研究，以期能探究城市發(fā)展進(jìn)程對(duì)周邊區(qū)域NDVI值的動(dòng)態(tài)影響。

參考文獻(xiàn)：

[1] 侯光雷，劉德贏，張正祥，等.松嫩平原不同氣候區(qū)植被NDVI對(duì)氣候變化的響應(yīng)[J].中國(guó)農(nóng)業(yè)氣象，2012，33（2）：271-277.

[2] 孫紅雨，王長(zhǎng)耀，牛 ?錚，等.中國(guó)地表植被覆蓋變化及其與氣候因子關(guān)系-基于NOAA時(shí)間序列數(shù)據(jù)分析[J].遙感學(xué)報(bào)，1998，2（3）：204-210.

[3] 徐 ?茜，任志遠(yuǎn)，楊 ?忍.黃土高原地區(qū)歸一化植被指數(shù)時(shí)空動(dòng)態(tài)變化及其與氣候因子的關(guān)系[J].陜西師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，40（1）：82-87.

[4] 李蕙敏，劉洪斌，武 ?偉.近10年重慶市歸一化植被指數(shù)變化分析[J].地理科學(xué)，2010，30（1）：119-123.

[5] 張 ?波，郭晉平，劉艷紅.太原市城市綠地斑塊植被特征和形態(tài)特征的熱環(huán)境效應(yīng)研究[J].中國(guó)園林，2010（1）：92-96.

[6] LUCIANO T， ROSA L. Quantifying intra-annual persistent behaviour in SPOT-VEGETATION NDVI data for Mediterranean ecosystems of southern Italy[J]. Remote Sensing of Environment，2006，101：95-103.endprint

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

[8] 崔曉臨，白紅英，尚小清.基于MODIS NDVI的秦嶺地區(qū)植被覆蓋變化研究[J].西北大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，42（6）：1021-1026.

[9] 李秀花，師慶東，常順利，等.1981-2001年中國(guó)西北干旱區(qū)NDVI變化分析[J]. 干旱區(qū)地理，2008，31（6）：940-944.

[10] 劉綠柳，許紅梅.黃河流域主要植被類型NDVI變化規(guī)律及其與氣象因子的關(guān)系[J].中國(guó)農(nóng)業(yè)氣象，2007，28（3）：334-337.

[11] 李雙成，劉逢媛，戴爾阜，等.NDVI與氣候因子耦合關(guān)系及其地域差異的定量遞歸分析-以云南省縱向嶺谷區(qū)為例[J].北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，44（3）：483-489.

[12] 張笑鶴.西南地區(qū)NDVI和NPP時(shí)空動(dòng)態(tài)及其與氣候因子相關(guān)性分析[D].北京：中國(guó)林業(yè)科學(xué)研究院，2011.

[13] 崔曉臨，白紅英，王 ?濤.秦嶺地區(qū)植被NDVI海拔梯度差異及其氣溫響應(yīng)[J].資源科學(xué)，2013，35（3）：618-626.

[14] 何 ?勇，董文杰，郭曉寅，等.基于MODIS的中國(guó)陸地植被生長(zhǎng)及其與氣候的關(guān)系[J].生態(tài)學(xué)報(bào)，2007，27（12）：5086-5092.

[15] 潘 ?薔，范文義，于海群，等.北京市植被指數(shù)變化與影響因素分析[J].北京林業(yè)大學(xué)學(xué)報(bào)，2012，34（2）：26-33.

[16] STOW D， DAESCHNER S， HOPE A，et al.Variability of the seas onally integrated normalized difference vegetation index a cross the north slope of Alaska in the 1990s[J].International Journal o f Remote Sensing，2003，24（5）：1111-1117.

[17] 蘇宗明.廣西天然植被類型分類系統(tǒng)[J]. 廣西植物，1998，18（3）：237-246.endprint