曾秀俐，朱霞，王維，王文杰，張哲

1.湖南科技大學建筑與城鄉(xiāng)規(guī)劃學院，湖南 湘潭 411201 2.青海省環(huán)境監(jiān)測中心站，青海 西寧 810007 3.中國環(huán)境科學研究院環(huán)境信息研究所，北京 100012

植被是氣候和人文因素對環(huán)境影響的敏感指標，也是聯(lián)結(jié)土壤、大氣和水文的自然紐帶[1]。植被指數(shù)是表征地表植被覆蓋、生長狀況的有效度量參數(shù)，其時空變化反映了生物量的動態(tài)演化信息，植被指數(shù)狀況越好，代表生態(tài)環(huán)境質(zhì)量狀況就越好。將植被指數(shù)應(yīng)用于區(qū)域生態(tài)環(huán)境監(jiān)測評估中具有重要的現(xiàn)實意義[2-3]。歸一化植被指數(shù)(normalized difference vegetation index, NDVI)是區(qū)域地表植被覆蓋度與植物長勢的函數(shù)，NDVI是反映植物生長狀況的最有效指標之一，該數(shù)值越高，表明地表葉面積指數(shù)越大，植被覆蓋度越高[4]。

遙感與GIS技術(shù)的發(fā)展為大區(qū)域、長時間序列的生態(tài)環(huán)境監(jiān)測提供了條件。目前利用NDVI進行植被監(jiān)測和評估的傳感器主要有NOAA AVHRR、SPOTVGT、MODIS、Landsat TM等[5-10]，比較而言，MODIS NDVI數(shù)據(jù)具有成像周期短、時間和空間分辨率較高，數(shù)據(jù)記錄完整等特點，對植被覆蓋狀況敏感區(qū)域的生態(tài)環(huán)境監(jiān)測具有明顯優(yōu)勢[11-12]，近年來被廣泛應(yīng)用于中長期植被變化趨勢研究。國內(nèi)外學者利用一元線性回歸[13-14]，Theil-Sen median趨勢分析和Mann-Kendall檢驗[15-17]，小波分析[18-19]和傅里葉變換[20-21]等方法研究植被覆蓋的時空變化特征，或利用相關(guān)分析等方法研究NDVI年內(nèi)年際變化與氣候、降水及人類活動間的關(guān)系[1-2,22-23]。由于一元線性回歸可以在像元尺度上表征一定時間序列植被覆蓋的變化趨勢，在國內(nèi)應(yīng)用相對成熟[24-28]，筆者擬采用該方法來分析研究區(qū)的植被變化趨勢。

邕江流域是我國典型的喀斯特地區(qū)之一，生態(tài)系統(tǒng)十分脆弱，在外部條件作用下很容易發(fā)生逆向演替。邕江流域植被狀況及其變化反映了該流域生態(tài)退化與恢復過程，并在一定程度上影響南寧市的經(jīng)濟社會和生態(tài)環(huán)境可持續(xù)發(fā)展，決定著南寧市的生態(tài)安全。目前研究主要以邕江水質(zhì)及其變化、邕江濱水區(qū)環(huán)境演變等為切入點，鮮有針對該流域長時間序列的植被定量評估和趨勢預(yù)測方面的研究。為此，筆者基于MODIS NDVI數(shù)據(jù)分析邕江流域2000—2010的植被時空特征及變化趨勢，評估了該流域的生態(tài)環(huán)境質(zhì)量狀況，以期為邕江流域植被的定量、動態(tài)監(jiān)測提供技術(shù)與理論支持，為該流域生態(tài)建設(shè)決策和區(qū)域生態(tài)環(huán)境演變研究提供科學依據(jù)，為喀斯特山區(qū)生態(tài)環(huán)境與經(jīng)濟社會可持續(xù)發(fā)展戰(zhàn)略研究提供參考。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

邕江流域位于南寧市南部，上起邕江上游左江的上中、右江的白馬至下游邕寧區(qū)與橫縣交界的六景鄉(xiāng)道莊村之間的水域以及地面水依自然地勢直接流入該水域的集水陸域(圖1)。邕江全長113.8 km，控制流域面積5 688.52 km2。邕江流域地處亞熱帶北側(cè)內(nèi)部邊緣，距北部灣海洋較近，受海洋季風影響，雨量充沛，霜少無雪，年均氣溫21.7 ℃，年均降雨量1 304.2 mm；區(qū)域內(nèi)地形起伏多變，地貌類型復雜多樣，包括平地、低山、石山、丘陵和臺地等；土壤類型以赤紅壤和紫色土為主，還分布有水稻土、石灰土等；植被資源十分豐富，植被類型包括闊葉林、針葉林、栽培植被、灌木和草叢等，流域所處的南寧市素有“綠都”的美譽。

圖1 邕江流域示意Fig.1 Location of Yongjiang Basin

邕江流域共涉及南寧市10個行政區(qū)縣(興寧區(qū)、江南區(qū)、青秀區(qū)、西鄉(xiāng)塘區(qū)、邕寧區(qū)、良慶區(qū)、隆安縣、武鳴縣、賓陽縣和橫縣)，總?cè)丝诩s280萬。近年來經(jīng)濟發(fā)展迅速，截至2010年底，邕江流域生產(chǎn)總值達893億元，為2005的4.5倍；三大產(chǎn)業(yè)所占比例為6.91∶34.80∶58.29，其中，工業(yè)產(chǎn)值在第二產(chǎn)業(yè)中所占比例不斷增加。

1.2 遙感數(shù)據(jù)處理

研究所用的遙感數(shù)據(jù)源自NASA提供的2000—2010年16 d最大合成的MODIS產(chǎn)品數(shù)據(jù)植被指數(shù)MODQ1，該數(shù)據(jù)空間分辨率為250 m。用MRT(MODIS ReProjection Tools)軟件對數(shù)據(jù)進行格式和地圖投影轉(zhuǎn)換，把HDF格式轉(zhuǎn)化為Geotiff格式，把Sinusoidal地圖投影轉(zhuǎn)換為WGS84Albers Equal Area Conic投影，采用最鄰近法進行重采樣；為了有效去除云、大氣和太陽高度角等產(chǎn)生的影響，進一步去除噪聲點并提高數(shù)據(jù)的可靠性，研究區(qū)的逐月NDVI通過對同屬一個月的2期NDVI影像進行最大值合成得到，對同年度的逐月NDVI影像進行加權(quán)平均得出該年度的NDVI[29]。

1.3 研究方法

1.3.1 回歸分析

對一組時間自變量(t)與NDVI因變量(NDVIt)數(shù)據(jù)，可以利用最小二乘法，計算出數(shù)據(jù)集上所有像元的NDVI與時間的回歸系數(shù)(k)：

式中，t為時間，a；NDVIt為t年的NDVI；n為樣本量。若k>0，說明該單元格的植被狀況趨好；若k<0，則說明該單元格的植被狀況趨于惡化。

對t與NDVIt間的線性關(guān)系進行了顯著性檢驗(F檢驗)[30-31]：

1.3.2 相關(guān)系數(shù)分析

采用Pearson相關(guān)系數(shù)建立NDVI與t的相關(guān)關(guān)系，分析NDVI隨t的長時間序列變化特征和趨勢，其公式為：

1.3.3 變異系數(shù)

變異系數(shù)可集中反映NDVI數(shù)據(jù)的離散程度，體現(xiàn)植被的年際波動性，其公式為：

CVNDVI=σNDVI

2 結(jié)果與分析

2.1 植被指數(shù)空間分布特征

2000—2010年邕江流域土地利用與年均NDVI空間分布的相互關(guān)系見圖2和圖3。邕江流域NDVI總體呈現(xiàn)四周高，中部低的空間分布特征，其與邕江流域地形及土地利用植被覆蓋分布格局相吻合；森林分布區(qū)NDVI普遍較高，草本濕地和耕地、草地分布區(qū)次之，建設(shè)用地、水域等其他類型NDVI相對較低。邕江流域中部為南寧中心城市區(qū)，由于市中心人口密集，人為干擾因素不斷加劇，NDVI明顯低于周邊區(qū)域；北部地區(qū)主要植被類型為森林和耕地，高峰嶺—白花山—三狀嶺水源涵養(yǎng)林建設(shè)與保護區(qū)NDVI明顯較高；流域東南部地貌主要為丘陵，林地和農(nóng)作物分布廣泛，尤其是大王灘—鳳亭河水庫飲用水源地保護區(qū)森林覆蓋率較高，NDVI集中在0.4～0.6，其他區(qū)域主要植被類型為耕地和草地，建設(shè)用地零星分布，年均NDVI為0.3～0.5。

圖2 邕江流域土地利用現(xiàn)狀Fig.2 Land use status in Yongjiang Basin

分級統(tǒng)計結(jié)果表明，邕江流域植被狀況總體較好，植被覆蓋區(qū)域占流域總面積的98%以上。其中NDVI小于0.3的僅占流域總面積的1.16%；NDVI為0.3～0.6的占流域總面積的80.32%；大于0.6的占流域總面積的18.52%。

2.2 植被指數(shù)時間變化特征

圖4 2000—2010年邕江流域年際NDVI變化Fig.4 Inter-annual change of NDVI in Yongjiang Basin from 2000-2010

以年均NDVI影像中植被覆蓋區(qū)域的像元NDVI平均值代表當年植被的整體狀態(tài)，研究邕江流域2000—2010年植被年際變化特征(圖4)。由圖4可知，邕江流域年均NDVI為0.42～0.54，2002年達到最大值，為0.54，2004年最小，為0.42；2000—2004年期間，NDVI變化沒有明顯規(guī)律，2004年以后NDVI表現(xiàn)出波動性緩慢增長趨勢。

分析邕江流域氣候因素與人為因素對植被的影響，結(jié)果表明，2000—2005年，邕江流域年均氣溫、土地利用類型及數(shù)量變化不大，年均降雨量是影響該流域植被狀況的主要因素。邕江流域降水量常年均值為1 328.2 mm，2002年邕江流域年平均降雨量達1 412 mm，約為常年均值的1.1倍；2003—2004年降水明顯減少；2004年平均年降水量1 193 mm，較常年均值偏少135 mm。2005年以來，溫度和降水變化幅度不大，人為干擾對植被的影響遠超過氣候因素，植樹造林、創(chuàng)建自然生態(tài)保護區(qū)等生態(tài)治理措施是植被狀況恢復的主要動因。

2.3 植被指數(shù)年際波動特征

年度NDVI的波動是植被受到干擾后在生產(chǎn)力和植被覆蓋度上的年際波動的重要體現(xiàn),是植被群落健康狀況的重要標志[24]。波動值越大說明植被群落越不穩(wěn)定，波動值越小說明植被群落狀態(tài)越穩(wěn)定。用CVNDVI分析邕江流域內(nèi)NDVI的波動特征，根據(jù)計算結(jié)果，借鑒前人的劃分經(jīng)驗[32]，將CVNDVI劃分為5個等級：高波動(CVNDVI≥0.20)，較高波動(0.15≤CVNDVI<0.20)，中波動(0.10≤CVNDVI<0.15)，較低波動(0.05≤CVNDVI<0.10)和低波動(CVNDVI<0.05)。分級統(tǒng)計結(jié)果如表1所示。由表1可知，92.87%的植被覆蓋區(qū)域NDVI呈中低波動，高度波動的區(qū)域僅占7.13%，說明邕江流域植被覆蓋波動性較小，總體穩(wěn)定性較好。

表1 邕江流域NDVI波動特征

邕江流域NDVI的變異系數(shù)分析結(jié)果表明(圖5)，NDVI波動性大小與植被類型密切相關(guān)。高波動區(qū)域主要分布在流域中部和西北部，尤其是位于流域中部的南寧中心城市區(qū)NDVI波動明顯，該區(qū)域人口密集，建設(shè)用地所占比例大，人為干擾劇烈，植被群落穩(wěn)定性較差；低波動地區(qū)主要分布在流域東北部的高峰嶺—白花山—三狀嶺水源涵養(yǎng)林建設(shè)與保護區(qū)和西南部的大王灘—鳳亭河水庫飲用水源地保護區(qū)，是由于該區(qū)域主要植被覆蓋類型為林地，穩(wěn)定性好，受干擾程度低，植被群落穩(wěn)定性較好；中度波動地區(qū)呈零散分布，主要植被類型為耕地和草地。

圖5 邕江流域NDVI變異系數(shù)分布Fig.5 Spatial distribution of variation coefficient of NDVI in Yongjiang Basin

2.4 植被指數(shù)變化趨勢分析

利用一元線性回歸分析邕江流域2000年以來的NDVI變化，并借鑒前人植被指數(shù)退化的劃分標準[27-28]，結(jié)合研究區(qū)實際情況，對邕江流域2000—2010年的植被覆蓋變化程度進行劃分，回歸分析結(jié)果如表2和圖6所示。

圖6 邕江流域趨勢分析結(jié)果Fig.6 Trends of inter-annul NDVI change in Yongjiang Basin

k變化趨勢面積∕km2所占比例∕%<-0.009嚴重退化234.064.13-0.009～-0.003中度退化522.319.23-0.003～0.003基本不變1 734.9430.620.003～0.009中度改善2 450.1943.24≥0.009明顯改善724.3812.78

由表2可知，邕江流域植被狀況有所改善的區(qū)域占流域總面積的56.02%，主要分布在邕江流域東南部的良鳳江及八尺江周邊；30.62%的區(qū)域植被覆蓋基本不變，僅有13.36%的區(qū)域生態(tài)環(huán)境退化，主要分布在南寧中心城區(qū)周邊及流域上游的左江、右江匯入口等區(qū)域，尤其是南寧中心城市區(qū)植被退化嚴重，與近年來該區(qū)域城鎮(zhèn)化加快推進密切相關(guān)。總體而言，植被處于增加趨勢的地區(qū)面積遠遠大于處于減少趨勢的地區(qū)面積。

回歸顯著性檢驗(F檢驗)結(jié)果表明，邕江流域有63.28%的區(qū)域通過了a=0.25的顯著性檢驗(可信度為75%以上)，其中有25.33%的區(qū)域通過了a=0.05的顯著性檢驗(可信度為95%以上)。

邕江流域2000—2010年NDVI與t的Pearson相關(guān)系數(shù)表明，邕江流域75%以上的區(qū)域NDVI與時間呈正相關(guān)關(guān)系，說明邕江流域大部分區(qū)域的植被覆蓋狀況在逐年增強，區(qū)域生態(tài)環(huán)境逐步改善。正相關(guān)區(qū)域主要集中在邕江流域下游的大王灘—鳳亭河水庫飲用水源地保護區(qū)和八尺江周邊，而南寧中心城市區(qū)NDVI與時間呈明顯的負相關(guān)關(guān)系，這與近年來南寧市城鎮(zhèn)化步伐加快導致城區(qū)周邊植被覆蓋度降低密切相關(guān)。Pearson相關(guān)系數(shù)分析結(jié)果與回歸分析結(jié)果保持一致。造成流域內(nèi)整體植被增加、生態(tài)環(huán)境改善的主要動因是近年來，南寧市積極創(chuàng)建國家級生態(tài)示范區(qū)，大規(guī)模植樹造林，強化城市生態(tài)屏障，打造“中國綠城”。

3 結(jié)論

(1)從空間分布上分析，邕江流域植被呈四周高，中部低的分布特征，這種格局與邕江流域的地形特點和土地利用方式相吻合。流域生態(tài)環(huán)境總體較好，植被覆蓋區(qū)域占流域總面積的98%以上。

(2)從時間變化上分析，2000—2010年邕江流域植被覆蓋區(qū)域年均NDVI為0.42～0.54，其中2000—2004年植被覆蓋區(qū)域NDVI波動較大，無明顯的趨勢特征，這與2002—2004年流域降水變化異常有關(guān)；2004年以后，植被覆蓋區(qū)域NDVI呈波動性緩慢增長趨勢，主要是由于南寧市近年來生態(tài)治理取得生態(tài)調(diào)節(jié)效應(yīng)在NDVI變化特征中得以顯現(xiàn)。

(3)從NDVI波動特性上分析，邕江流域植被覆蓋總體穩(wěn)定性較好，其中92.87%的植被覆蓋區(qū)域NDVI呈中低波動；而高波動區(qū)域僅占植被覆蓋區(qū)域總面積的7.13%。植被狀況的穩(wěn)定性與植被類型及數(shù)量密切相關(guān)。

(4)從NDVI變化趨勢上分析，2000—2010年邕江流域地表植被覆蓋改善的區(qū)域遠遠大于植被退化的區(qū)域，其中，改善的區(qū)域占植被覆蓋區(qū)域總面積的56.02%，退化的區(qū)域占13.36%，30.62%的區(qū)域呈穩(wěn)定不變的狀態(tài)。

[1] 劉綠柳,肖風勁.黃河流域植被NDVI與溫度、降水關(guān)系的時空變化[J].生態(tài)學雜志,2006,25(5):477-481.

[2] 孫紅雨,王常耀,牛錚,等.中國植被覆蓋變化及其與氣候因子關(guān)系:基于NOAA數(shù)據(jù)[J].遙感學報,1998,2(3):204-210.

[3] HOPE A S,BOYNTON W L,STOW D C,et al.Interannual growth dynamics of vegetation in the Kuparuk River watershed, Alaska based on the normalized difference vegetation index[J].International Journal of Remote Sensing,2003,24(7):3413-3425.

[4] 趙英時.遙感應(yīng)用分析原理與方法[M].北京:科學出版社,2003.

[5] 王強,張勃,戴聲佩.基于GIMMS AVHRR NDVI數(shù)據(jù)的三北防護林工程區(qū)植被覆蓋動態(tài)變化研究[J].資源環(huán)境與發(fā)展,2011(2):12-19.

[6] JIANG X G,WANG D,TANG L L,et al.Analysing the vegetation cover variation of China from AVHRR-NDVI data[J].International Journal of Remote Sensing,2008,29(17):5301-5311.

[7] ANYAMBA A,TUCKER C J.Analysis of Sahelian vegetation dynamics using NOAA-AVHRR NDVI data from 1981-2003[J].Journal of Arid Environments,2005,63(3):596-614.

[8] 宋富強,邢開雄,劉陽,等.基于MODIS NDVI的陜北地區(qū)植被動態(tài)監(jiān)測與評價[J].生態(tài)學報,2011,31(2):354-363.

[9] HUETE A R,DIDAN K,SHIMABOKURO Y,et al.Regional Amazon basin and global analysis of MODIS vegetation indices:early results and comparisons with AVHRR[C]IGARSS 2000:IEEE 2000 international geoscience and remote sensing symposium - taking the pulse of the planet:the role of remote sensing in managing the environment.Honolulu:IEEE,2000,6(2):536-538.

[10] PETTORELLI N,VIK J O,MYSTERUD A,et al.Using the satellite-derived NDVI to assess ecological responses to environmental change[J].Trends in Ecology & Evolution,2005,20(9):503-510.

[11] 陳燕麗,龍步菊,潘學標,等.MODIS NDVI和AVHRR NDVI對草原植被變化監(jiān)測差異[J]遙感學報,2011,15(4):831-845.

[12] 王正興,劉闖,HUETE A,等.植被指數(shù)研究進展:從AVHRR-NDVI到MODIS-EVI[J].生態(tài)學報,2003,23(5):979-987.

[13] 宋怡,馬明國.基于SPOT VEGETATION數(shù)據(jù)的中國西北植被覆蓋變化分析[J].中國沙漠,2007,27(1):89-93.

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

[15] 王佃來,劉文萍,黃心淵.基于Sen+Mann-Kendall的北京植被變化趨勢分析[J].計算機工程與應(yīng)用,2013,49(5):13-17.

[16] POULIOT D,LATIFOVIC R,OLTHOF I.Trends in vegetation NDVI from 1 km AVHRR data over Canada for the period 1985-2006[J].International Journal of Remote Sensing,2009,30(1):149-168.

[17] FENSHOLT R,LANGANKE T,RUSMUSSEN K,et al.Greenness in semi-arid areas across the globe 1981-2007:an Earth Observing Satellite based analysis of trends and drivers[J].Remote Sensing of Environment,2012,121:144-158.

[18] MARTINEZ B,GILABERT M A.Vegetation dynamics from NDVI time series analysis using the wavelet transform[J].Remote Sensing of Environment,2009,113(9):1823-1842.

[19] 宋濤.基于小波變換的植被覆蓋變化研究[D].福建:福建師范大學,2011.

[20] 黃勇,姚筠,郁凡.基于傅立葉變換的多時次NDVI地表覆被分類方法[J].中國農(nóng)學通報,2009,25(3):260-264.

[21] 王紅說,黃敬峰.基于MODIS NDVI時間序列的植被覆蓋變化特征研究[J].浙江大學學報,2009,35(1):105-110.

[22] 鄭有飛,劉宏舉,吳榮軍,等.貴州省NDVI變化及其與主要氣候因子的相關(guān)性[J].生態(tài)與農(nóng)村環(huán)境學報,2009,25(1):12-17.

[23] 江東,王乃斌,楊小喚,等.植被指數(shù)-地面溫度特征空間的生態(tài)學內(nèi)涵及其應(yīng)用[J].地理科學進展,2001,20(2):146-152.

[24] 王情,劉雪華,呂寶磊.基于SPOT-VGT數(shù)據(jù)的流域植被覆蓋動態(tài)變化及空間格局特征:以淮河流域為例[J].地理科學進展,2013,32(2):270-277.

[25] 毛德華,王宗明,羅玲,等.基于MODIS和AVHRR數(shù)據(jù)源的東北地區(qū)植被NDVI變化及其與氣溫和降水間的相關(guān)分析[J].遙感技術(shù)與應(yīng)用.2012,27(1):77-85.

[26] 安佑志,劉朝順,施潤和,等.基于MODIS時序數(shù)據(jù)的長江三角洲地區(qū)植被覆蓋時空變化分析[J].生態(tài)環(huán)境學報,2012,21(12):1923-1927.

[27] 邱海軍,曹明明.基于SPOT VEGETATION數(shù)據(jù)的中國植被覆蓋時空變化分析[J].資源科學,2011,33(2):335-340.

[28] HERRMANN S M, ANYAMBA A, TUCKER C J.Recent trends in vegetation dynamics in the African sahel and their relationship to climate[J].Global Environmental Change,2005,15(4):394-404.

[29] 何彬方,馮妍,吳文玉,等.安徽省近十年植被指數(shù)時空變化特征[J].生態(tài)學雜志,2010,29(10):1912-1918.

[30] 劉賢趙,張安定,李嘉竹.地理學的數(shù)學方法[M].北京:科學出版社,2009:29-30.

[31] 南英子.簡單線性相關(guān)與回歸分析中各種檢驗及其相互關(guān)系的探討[J].統(tǒng)計與決策,2011(2):26-28.

[32] 張雪艷,胡云峰,莊大方,等.蒙古高原NDVI的空間格局及空間分異[J].地理研究,2009,28(1):10-18. ○