世界自然遺產地荔波縣植被覆蓋時空變化研究

胡向紅　許玉鳳　陳如霞

摘要 基于2000—2015年MODIS NDVI數據，在GIS技術支持下，采用最大值合成法、趨勢分析和變異系數等方法，研究了世界自然遺產保護區荔波縣植被覆蓋的時空變化特征。結果表明：荔波縣植被覆蓋較好的區域主要分布在東部和西部區域，城鎮及交通線沿線的植被覆蓋較差。研究期間植被覆蓋總體呈上升趨勢，其中2000—2006年間植被覆蓋變化比較緩慢，2006—2012年間植被上升趨勢最顯著，2013—2015年間植被上升趨勢變緩。研究區植被覆蓋變化空間差異性顯著，自然保護區的植被覆蓋較好，植被覆蓋比較穩定；南部區域的植被退化應該引起重視。

關鍵詞 遙感數據；植被覆蓋；時空變化；世界自然遺產保護區；貴州荔波

中圖分類號 P966；TP79 文獻標識碼 A 文章編號 1007-5739（2017）04-0223-03

Abstract Based on MODIS NDVI data from 2000 to 2015，using the maximum synthesis，trend analysis and variation coefficient methods，the spatial and temporal change of vegetation cover in the world natural heritage site Libo County were studied under the GIS support. The results showed that the vegetation coverage in Libo County was mainly distributed in the eastern and western regions，and the vegetation coverage along the urban and transportation lines was poor. The vegetation coverage generally showed a rising trend during the study period. The rising trend was significant from 2006 to 2012，and relatively slow from 2000 to 2006，and slower from 2013 to 2015. The spatial variation of vegetation coverage was significant in the study area，and the vegetation coverage was good and relatively stable in the nature reserves. It should be paid more attention to the vegetation degradation in the southern region.

Key words remote sensing data；vegetation coverage；temporal and spatial variation；the world natural heritage site；Libo Guizhou

地表植被在調節區域氣候、涵養水源、保持水土和維持生物多樣性等方面具有重要的支撐與促進作用，對地方、區域乃至全球生態系統的狀況、環境壓力和景觀變化有指示作用[1-2]。植被變化是自然因素和社會因素共同作用的結果，是全球變化的主要表現形式[3-4]。植被覆蓋變化是區域水土流失、生態環境惡化的根源，因此有必要對植被動態變化進行監測研究[5]。

遙感技術已經成為植被動態變化監測研究的主要手段，主要用于土地覆蓋變化[6]、植被指數監測[7]等。大區域范圍內的植被變化監測主要采用植被指數，其中歸一化植被指數（NDVI）能夠準確反映地表植被的生長狀況以及覆蓋程度等特征，應用最為廣泛[5-8]。國內外利用植被指數監測植被變化的研究成果較多[9-14]，在定量監測的研究方法和研究體系等方面有了長足的發展，但是大多數是針對北方平原地區、河流流域地區的研究，而對于西南山地的相關研究較少[15]。

荔波縣地處貴州南部，旅游資源豐富，隨著“山地公園省、多彩貴州風”宏觀構想目標的建立，對荔波縣世界自然遺產保護區的保護和開發尤為重要，遙感生態監測作用顯得尤為突出，但目前利用遙感數據研究荔波縣植被覆蓋變化的文獻很少。因此，本文利用2000—2015年的遙感數據對荔波縣植被覆蓋狀況進行了動態監測分析，目的在于了解荔波縣植被資源的動態變化及其變化趨勢，為世界自然遺產地荔波縣的生態保護提供科學依據。

1 研究區概況

荔波縣地處貴州省南部，隸屬于黔南布依族苗族自治州（圖1），范圍在東經107°37′～108°18′、北緯25°07′～25°09′之間，面積2 431.8 km2，包括1個街道、5個鎮、2個鄉；少數民族人口占總人口的90%以上，人口較多的少數民族有布依族、水族、瑤族、苗族等，布依族是荔波縣第一大民族，占全縣總人口的65%（2010年數據），是典型的少數民族聚集的縣區。

荔波縣東北與黔東南苗族侗族自治州的從江縣、榕江縣接壤，東南與廣西壯族自治區的環江縣、南丹縣毗鄰，西與獨山縣相連，北與三都水族自治縣交界。植物資源豐富，境內有樟江國家級風景名勝區和國家級茂蘭自然保護區。

荔波縣屬于中亞熱帶季風濕潤氣候區，氣溫分布的總趨勢是南高北低，全縣年平均氣溫18.3 ℃，常年最熱月為7月，最冷月為1月。7月平均氣溫為27 ℃以下，極端最高氣溫不超過40.0 ℃；1月平均氣溫各地均在5.5 ℃以上，極端最低氣溫在-10 ℃以上?？h內霜期少，無霜期在270 d以上。降水主要集中在夏季。6—8月各月雨量在2 000 mm以上，占全年總雨量的50%左右；冬季（12月至翌年2月）僅占全年總雨量的5%左右；秋季（9—11月）占15%左右；春季（3—5月）海洋季風逐漸增強，降水占全年雨量的30%左右；于4月下旬進入雨季，4—10月降水量占全年總降水量的81%。

荔波縣地處貴州廣西省區的交界處，地形聚集了貴州高原和廣西丘陵的綜合優勢，既有高山也有丘陵，屬于喀斯特地形地貌。荔波境內河流均屬龍江水系，歸珠江流域。荔波縣河流落差大，水能理論蘊藏量豐富。

2 數據與方法

2.1 數據來源

采用美國NASA提供的MODIS NDVI數據，分辨率為16 d，空間分辨率為250 m，時間跨度為2000—2015年，該數據已經通過了輻射校正、云體掩膜、大氣校正等處理。首先利用MRT（MODIS Reprojection Tools）軟件對遙感數據進行格式轉換、投影轉換和鑲嵌工作，統一采用WGS_84地理坐標系、UTM投影；利用行政區劃圖邊界為掩膜裁剪出研究區邊界圖，通過批處理得到荔波縣2000—2015年的柵格圖像數據集。

2.2 最大值合成法

最大值合成法MVC（maximum value composite）是將全年23期遙感影像取最大值，可以在一定程度上去除太陽高度角、云、大氣等方面的影響[16-18]。

2.3 趨勢分析

一元線性回歸趨勢分析可以模擬每個柵格單元的變化趨勢，被廣泛應用于基于柵格的植被綠度變化速率[12]。對各像元上的NDVI值進行一元線性回歸運算，可以了解多年植被空間變化規律，計算公式為：

k=■（1）

式（1）中，k為變化趨勢斜率；n為監測年數；fi為第i年的植被覆蓋狀況，本文中是指NDVI值。k>0表示研究期間NDVI呈上升趨勢，植被覆蓋有所改善；k<0表明研究期間NDVI值呈下降趨勢，植被覆蓋有所退化；k=0為穩定狀態，植被覆蓋狀況基本沒有改變。

2.4 變異系數（穩定性分析）

植被NDVI的年際波動是植被受到干擾后在生產力和植被覆蓋度上波動的重要體現，是植被群落健康狀況的重要標志[19]。波動值大說明植被群落受到干擾強度大、不穩定（植被狀況趨于恢復或退化）；波動值小說明植被群落狀況穩定[20]。

利用研究期間植被NDVI的變異系數來反映研究區域內植被變化的波動程度，其計算公式為：

S=■（2）

CV=■（3）

式（2）、（3）中，CV為變異系數，S為標準差，NDVI為多年植被NDVI的平均值。

3 結果與分析

3.1 植被覆蓋的年際變化

2000—2015年荔波縣年均NDVI值呈波動上升趨勢，增長率為0.018/10年（R2=0.574 7），見圖2（a）。其中，2000—2006年期間植被覆蓋增加緩慢，幾乎沒有增長（R2=0.046 1），見圖2（b）；2006—2012年期間植被覆蓋增加迅速，增長率為0.049/10年（R2=0.691 4），見圖2（c）；2013—2015年期間植被覆蓋增加速度稍有減緩，增長率為0.015/10年（R2=0.994 1），見圖2（d）。近16年間NDVI值隨時間的變化呈顯著上升趨勢，通過了置信度99%的顯著性檢驗，說明研究區經過近10多年的退耕還林、還草等生態保護建設工程，植被狀況顯著改善。2000—2006年研究區的植被覆蓋基本上無明顯變化，這一研究結論與吳昌廣等人的研究結果一致，可能與西南地區的嚴重干旱有關[21]。

3.2 植被覆蓋空間分布及其變化

3.2.1 荔波縣植被覆蓋空間分布。荔波縣的植被覆蓋總體狀況較好，最低值在0.5以上，但是有明顯的空間差異，植被覆蓋較好的區域主要分布在南部、東北部和西部的部分山地，在城市、鄉鎮、水域等分布的區域及交通沿線地區的植被覆蓋狀況較低（圖3）。

3.2.2 2000—2015年植被覆蓋的空間變化。荔波縣植被NDVI變化趨勢存在明顯的空間差異。2000—2015年研究區植被退化面積占研究區面積的23.09%，主要分布在荔波中部、南部、東北和西北等地區，在城鎮地區、交通沿線表現更加明顯。中部地區受建設用地增加的影響，植被退化比較明顯；南部、東北和西北部主要是農業用地分布范圍較大。

植被覆蓋增加的地區占研究區面積的75.06%，主要分布在荔波縣東部和西部地區；穩定不變的區域占研究區面積的1.31%（圖4）。

樟江和茂蘭等自然保護區的植被變化呈上升趨勢，而南部植被覆蓋較好的區域，植被覆蓋呈退化趨勢，需要引起重視。

3.3 植被穩定性的空間格局

多年植被NDVI的變異系數反映了植被的穩定性，變異系數越大則植被穩定性越差?？傮w來看荔波縣植被穩定性較好。植被較好的地區，植被變異系數較小，較穩定；植被覆蓋較差的地區，變異系數較大，不穩定（圖5）。

在水域附近區域的植被變異系數較大，表明這些地區受濕地環境的影響，對環境脅迫比較敏感。

人口增加及城鎮化發展的加快、城市人口的增加對植被覆蓋產生了影響，針對人為干涉產生的結果，建議主要采用自然力進行生態修復，特別是出現輕微退化的生態脆弱區，要繼續加大封山育林的管理力度。對于嚴重退化的區域要結合人工修復，開展生態修復工程來改善其生態環境。

4 結論

本文利用2000—2015年MODIS數據，采用最大值合成法、趨勢分析和變異系數等方法，研究了世界自然遺產保護區荔波縣植被覆蓋的時空變化及變化趨勢，得出以下研究結論：

（1）荔波縣植被覆蓋較好的區域主要分布在東部和西部區域，城鎮、交通線附近分布區的植被覆蓋較差。

（2）研究期間荔波縣植被覆蓋呈上升趨勢，不同階段的變化趨勢有所區別。2000—2006年期間，植被覆蓋變化比較緩慢，可能與這一階段的干旱狀況有關；2006—2012年間植被上升趨勢最顯著；2013—2015年間植被變化趨勢開始變緩。

（3）研究期間，研究區植被覆蓋變化有顯著的空間差異性，自然保護區分布區的植被覆蓋較好，植被覆蓋比較穩定；南部區域的植被退化應該引起重視。

（4）造成植被退化的原因主要是人類活動，地形地貌也是造成植被變化的原因之一，氣溫、降水等變化也對植被變化的影響作用十分顯著。具體影響原因及其影響程度是今后進一步研究的重點。

5 參考文獻

[1] 孫紅雨，王長耀，牛錚，等.中國地表植被覆蓋變化及其與氣候因子關系[J].遙感學報，2009，2（3）：204-210.

[2] 張善紅.延安市植被變化及其對氣溫變化的響應[J].現代農業科技，2014（17）：269-272.

[3] O′CALLAGHAN J F，MARK D M.The extraction of drainage networks from digital[J].ComputerVision，Graphics，and Image Processing，1984， 28：323-344.

[4] JENSON S K，DOMINGUE J O.Extraction topographic structure from elevation data for geographic information system analysis[J].Photogramm-etric Enginerimg and Remote Sensing，1988，54（11）：1593-1600.

[5] 楊嘯.基于時序NDVI的湖北省植被覆蓋動態變化監測分析[J].長江流域資源與環境 2013，22（2）：226-231.

[6] 李琳，譚炳香，馮秀蘭.北京郊區植被覆蓋變化動態遙感監測：以懷柔區為例[J].農業網絡信息，2008（6）：38-42.

[7] 郭鈮.植被指數及其研究進展[J].干旱氣象，2003，21（4）：71-75.

[8] 張思詩，杜曉初.基于NDVI的三峽庫區植被時空變化規律研究[J].湖南大學學報（自然科學版），2012，34（3）：317-319.

[9] 杜子濤，占玉林，王長耀.基于NDVI序列影像的植被覆蓋變化研究[J].遙感技術與應用，2008，23（1）：47-50.

[10] 王蕊，李虎.2001—2010年蒙古國MODIS-NDVI時空變化監測分析[J].地球信息科學學報，2011，13（5）：665-671.

[11] 邱海軍，曹明明.基于SPOT VEGETATION數據的中國植被覆蓋時空變化分析[J].資源科學，2011，33（2）：335-340.

[12] HOPE A，BOYNTON W，STOW D，et al.Inter-annual growth dynamics of vegetation in the Kuparuk River watershed based on the normalized difference vegetation vegetation index[J].International Journal of Remote Sensing，2003，24（17）：3413-3425.

[13] DURO D，COOPS N C，WULDER M A，et al.Development of a large area biodiversity monitoring system driven by remote sensing[J].Progress in Physical Geography，2007，31（3）：235-260.

[14] 馬保東，陳紹杰，吳立新，等.基于SPOT VEGETATION的礦區植被遙感監測方法[J].地理與地理信息科學，2009，25（1）：84-86.

[15] 李威，李月臣，唐誼娟，等.基于MODIS-NDVI的貴州省2001—2010年植被覆蓋變化分析及分異特征[J].重慶師范大學學報（自然科學版），2015，32（4）：146-151.

[16] HOLBEN B N.Characteristics of maximum-value composite images for temporal AVHRR data[J].Remote Sensing，1986，7：1435-1445.

[17] TUCKER C J，NEWCOMB W W，DREGNE H E.AVHRR datasets for determination of desert spatital extend international[J].Journal of Remote Sensing，1994，17：354-356.

[18] 趙英時.遙感應用分析原理與方法[M].北京：科學出版社，2003.

[19] 侯美婷，趙海燕，王箏，等.基于衛星遙感的植被NDVI對氣候變化響應的研究進展[J].氣候與環境研究，2013，18（3）：353-364.

[20] 于泉洲，梁春玲，劉煜杰，等.基于MODIS的山東省植被覆蓋時空變化及其原因分析[J].生態環境學報，2015，24（11）：1799-1807.

[21] 吳昌廣，周志翔，肖文發，等.基于MODIS NDVI的三峽庫區植被覆蓋動態監測[J].林業科學，2012，48（1）：22-28.