王醒，方榮杰，張帥普，黃卉，黃朗，吳忠軍

(1.桂林理工大學，廣西環(huán)境污染控制理論與技術重點實驗室，廣西 桂林 541004；2.桂林理工大學，巖溶地區(qū)水污染控制與用水安全保障協(xié)同創(chuàng)新中心，廣西 桂林 541004；3.桂林理工大學，旅游與風景園林學院，廣西 桂林 541004)

植被是陸地生態(tài)系統(tǒng)存在的基礎，也是生物、大氣、水分和土壤等要素相互聯結的樞紐，在維持生態(tài)平衡、調節(jié)氣候、涵養(yǎng)水源、保土固土等方面具有重要作用[1-3]。植被覆蓋度(fractional vegetation cover，Fvc)作為衡量植被覆蓋狀況的重要指標，是了解區(qū)域生態(tài)系統(tǒng)環(huán)境變化的重要途徑[4]，其變化特征可以指示區(qū)域生態(tài)環(huán)境變化狀況，被廣泛應用于氣候、生態(tài)、水文等研究領域[5-7]。

植被覆蓋度時空分布特征受到社會、人文、自然等因素的綜合影響，在不同空間位置對各影響因子的響應程度不同[8]。高程、坡度、坡向等地形因子是影響植被覆蓋變化的重要自然因素，對植被覆蓋度的影響主要體現為靜態(tài)響應[9-10]，即空間分布上植被覆蓋度對地形響應在不同空間地域上往往呈顯著差異[11]。有研究表明，在我國北方延河流域植被覆蓋度隨高程及坡度的升高總體呈先增加后減少的趨勢，陰坡的植被覆蓋度最高[12]；而在西南地區(qū)，高程大于4 000 m時，植被覆蓋度隨高程的升高逐漸下降，而且退化趨勢也最為顯著[13]。秦嶺山地植被覆蓋在低海拔區(qū)域呈減少趨勢，中海拔區(qū)呈明顯的上升趨勢，且北坡的植被覆蓋度在2 000 m以上的高海拔區(qū)域較為穩(wěn)定，而南坡在2 500～3 100 m區(qū)域內有較明顯的減小趨勢[14]。此外，不同地形因子對不同區(qū)域植被覆蓋度的主導作用不同[13]。熊俊楠等[15]和劉堯文等[16]在汶川地震災區(qū)、平潭島等地區(qū)研究發(fā)現高程和坡度比坡向更能影響植被覆蓋度的變化；而蔡宏等[17]在赤水河流域研究表明植被覆蓋度受坡向、坡度的影響明顯高于高程。植被覆蓋度與地形因子間的關系具有復雜性及不確定性，基于不同地區(qū)地形特征的差異性來探究植被覆蓋度的空間分布規(guī)律始終備受關注[9，18]。梯田是我國重點關注的山地地區(qū)農業(yè)耕種模式之一，其獨特的平面和剖面形態(tài)是由人工改變地形而形成。有大規(guī)模梯田的山丘地區(qū)因梯田梯式分布而具有獨特的地形地貌，梯田與周山森林植被的疊加形成穩(wěn)定的景觀格局，該格局可在減少水土流失等生態(tài)保護方面發(fā)揮強大的功能，因此山區(qū)梯田是具有生態(tài)效益和經濟效益并存區(qū)域[19-20]。盡管前人學者已對不同區(qū)域尺度的植被覆蓋度分布特征及其對地形因子的響應機制進行許多研究，但人工梯田地形的特殊性且地表覆蓋受農業(yè)耕種模式的影響嚴重，對具有獨特景觀格局的梯田地區(qū)植被覆蓋度變化及其對地形因子的響應仍需進一步探究。

龍脊梯田旅游開發(fā)始于1992年，是我國AAAA級重點農業(yè)旅游景區(qū)。龍脊梯田獨特的空間結構由森林、村寨及梯田組成，是典型自然-人工復合生態(tài)系統(tǒng)，但近些年來，隨著龍脊梯田旅游產業(yè)的快速發(fā)展，人類活動愈發(fā)增強，多處地表受到流水侵蝕、重力坍塌等影響，境內植被生態(tài)景觀日益破碎，高山區(qū)域的地形結構及植被覆蓋已隨區(qū)域的發(fā)展發(fā)生了明顯變化。因此，本文以龍脊梯田小流域為研究對象，利用1993—2019年3期Landsat遙感影像及DEM地形數據，深入研究龍脊梯田旅游區(qū)自開發(fā)以來植被覆蓋度的時空分布格局及其對地形的響應特征，以期為龍脊梯田生態(tài)環(huán)境建設及保護提供科學依據。

1 研究區(qū)概況

龍脊梯田位于廣西桂林市龍勝各族自治縣龍脊鎮(zhèn)龍脊山(109°32′～110°14′E、25°35′～26°17′N)，分為金坑(大寨)瑤族梯田觀景區(qū)和平安壯族梯田觀景區(qū)。采用高程數據提取龍脊鎮(zhèn)上游小流域作為研究范圍，其面積為126.80 km2。研究區(qū)域主要包含平安、中祿、大寨、江柳、黃江、馬海6個行政村。區(qū)域內崇山峻嶺，溪流縱橫。梯田主要分布在海拔300～1 100 m之間。海拔在700 m以上的山地占3/5以上，最大坡度達50°。該區(qū)域地處中亞熱帶季風氣候區(qū)，四季分明。夏季為東南風，冬季多西北風，風力一般1～3級。年平均氣溫14.4～16.9 ℃，最熱月(7月)平均氣溫25.4 ℃，最高氣溫32 ℃，最冷月(1月)平均氣溫7.1 ℃，最低氣溫-6 ℃。年降雨量1 600～1 700 mm，雨量豐沛，且主要集中在5—9月，年均日照時數1 225.7 h，平均無霜期290 d。植被主要有人工杉木林、竹林、混交林、水稻等，森林覆蓋率達75.6%以上。土壤類型主要為紅壤土、黃壤土及黃棕壤土。

2 數據與方法

2.1 數據來源與處理

2.1.1 數據來源

本次研究選取1993年9月、2006年9月、2019年10月3期影像數據，其中1993年9月和2006年9月為Landsat-5衛(wèi)星TM數據，2019年10月為Landsat-8衛(wèi)星OLI數據，數據均來源于美國地質調查局(United States Geological Survey，USGS)。本研究選取的3期遙感影像時間均處于9—10月，均為龍脊梯田的水稻生長末期，能較好反映出區(qū)域植被狀況。地形數據采用30 m分辨率的DEM數據，來源于中國科學院計算機網絡信息中心地理空間數據云平臺。

2.1.2 數據預處理

為減小光照、大氣等因素的影響，分析數據前先利用ENVI 5.3軟件對所獲取的遙感影像進行輻射定標和大氣校正、影像圖像剪裁等預處理，提高遙感數據的精度和質量，之后計算歸一化植被指數NDVI。基于DEM數據，利用ArcGis 10.2提取龍脊梯田分水嶺作為本次研究流域的矢量邊界。

2.2 研究方法

2.2.1 植被覆蓋度的提取

植被覆蓋度(Fvc)是植被冠層(包括葉、莖、枝)在地面上的垂直投影面積占統(tǒng)計區(qū)面積的百分比，其與歸一化植被指數NDVI之間存在一定的關系[21]。本研究基于標準化植被指數(NDVI)的像元二分模型進行龍脊梯田植被覆蓋度的反演，其原理是假設一個像元只包含純植被和純土壤兩個信息[5]。計算方法如下。

Fvc=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

式中:NDVI為單個像元NDVI值，NDVIsoil為裸土或無植被覆蓋區(qū)域的NDVI值，NDVIveg為完全被植被所覆蓋的像元的NDVI值。理論上NDVIsoil數值應接近于0，NDVIveg數值代表全植被覆蓋單位像元NDVI的最大值，但考慮到植被類型、噪聲、地形及影像質量等多種因素的影響，NDVIsoil與NDVIveg值會與實際值發(fā)生偏差，一般取一定置信度范圍內的最大值與最小值來代表。本文選取研究區(qū)遙感影像中NDVI值累積百分數為2%和98%附近的NDVI值來作為NDVIsoil和NDVIveg值。在ENVI 5.31系統(tǒng)中，通過遙感模型估算研究區(qū)1993年、2006年和2019年的植被覆蓋度。根據研究區(qū)實際植被覆蓋狀況，將每期植被覆蓋度劃分5個級別[22-23]：低植被覆蓋度(0～0.2)、較低植被覆蓋度(0.2～0.4)、中植被覆蓋度(0.4～0.6)、較高植被覆蓋度(0.6～0.8)、高植被覆蓋度(0.8～1)。

2.2.2 地形因子分級

根據研究區(qū)地形實際情況，將研究區(qū)的高程按300～500、500～700、700～900、900～1 100、1 100～1 300、1 300～1 500、1 500～1 700、>1 700 m劃分為8個高程帶?？紤]區(qū)域地形對植被生長的影響，以8°作為緩坡和斜坡界線[24]，按0°～8°、8°～15°、15°～25°、25°～35°、>35°劃分為5個坡度帶。坡向分為北、東北、東、東南、南、西南、西、西北面8個方向[17]。各地形等級面積所占研究區(qū)的比例見表1。

表1 研究區(qū)各地形等級面積占比

3 結果與分析

3.1 土地利用及植被覆蓋度變化

3.1.1 土地利用類型

為更深入了解植被覆蓋度的變化，本文采用監(jiān)督分類和目視解譯相結合的方法對龍脊梯田3期遙感數據進行土地利用類型解譯，根據《土地利用現狀分類》(GB/T21010—2017)，將區(qū)域為林地、耕地、建設用地、未利用地4個土地利用類型，如圖2所示。1993—2019年，林地分布最廣泛，其次是耕地，兩者占地面積高達89%以上；建設用地面積不大，平均面積比例僅為5.78%，但其分布逐漸擴大；未利用地面積最小，由裸土地及局部稀疏草地組成，主要分布在東北方向的高山區(qū)域。

不同土地利用類型變化差異明顯(表1)。26 a間，林地、建設用地面積總體增加，增加的大小表現為：建設用地>林地，而耕地面積總體減少，未利用地面積變化相對平穩(wěn)，1993—2006年，林地、未利用地面積明顯增加，耕地面積明顯減小，建設用地面積變化不明顯；2006—2019年，建設用地明顯增加，未利用地面積明顯減小，林地和耕地面積沒有太大變化。

圖1 1993、2006、2019年龍脊梯田土地利用類型

表2 1993—2019年土地利用類型面積占比

3.1.2 植被覆蓋度時空分布特征

1993—2019年龍脊梯田植被覆蓋整體狀況良好，平均植被覆蓋度達0.76以上(圖2、圖3)。各等級植被覆蓋度在空間上差異較大，較高植被覆蓋度和高植被覆蓋度兩等級分布較廣，幾乎遍及整個區(qū)域，主要集中在無梯田、無村民居住的中部及邊緣的高山區(qū)域，二者占區(qū)域總面積的比例為1993年87.77%、2006年88.96%、2019年87.94%，即龍脊梯田87%以上的區(qū)域植被覆蓋達到了較高水平以上；低、較低、中等植被覆蓋度區(qū)域面積較小，主要分布在人類活動較大、耕地用地較多的區(qū)域。

圖2 1993、2006、2019年龍脊梯田植被覆蓋度

圖3 1993—2019年龍脊梯田植被覆蓋度均值及不同等級面積比例

1993—2019年，植被覆蓋度總體變化呈“低、較低、高植被覆蓋度增加，中等、較高植被覆蓋度減小”的特征；植被覆蓋度在1993—2006年間發(fā)生了較大變化，較高、高級別植被覆蓋度分別呈明顯減小、增加趨勢，變化均達14%以上，低、較低、中植被覆蓋度變化相對較小，其中低、較低、中植被覆蓋度分別呈增加、增加、減小趨勢；各等級植被覆蓋度在2006—2019年間沒有太大變化。

3.2 植被覆蓋度面積轉移

表3展示了龍脊梯田不同等級植被覆蓋度轉移特征。1993—2019年高植被覆蓋度未發(fā)生級別變化的比例最大，高達83.27%；其次是低、較高植被覆蓋度，中等、較高植被覆蓋度未發(fā)生級別變化的比例最小。在植被覆蓋度發(fā)生級別變化中，1993—2006年低植被覆蓋度主要轉移至中等植被覆蓋度，而較低植被覆蓋度主要轉移至中等、較高植被覆蓋度，中等、較高植被覆蓋度主要轉移至比各自高一級別的植被覆蓋度；2006—2019年低、較低、中植被覆蓋度均主要轉移至較高植被覆蓋度，而較高植被覆蓋度主要轉移至高植被覆蓋度；整體而言，1993—2019年低、較低、中等植被覆蓋度均主要轉移至較高植被覆蓋度，較高植被覆蓋度主要轉移至高植被覆蓋度?？梢姡埣固萏锒嗵幍椭脖桓采w度級別有所提高，植被覆蓋度逐漸改善，植被生長整體態(tài)勢良好。

表3 1993—2019年不同植被覆蓋度面積轉移比例

3.3 植被覆蓋度地形響應分析

3.3.1 高程響應

在不同高程范圍內，水、熱、人類活動等是影響植被覆蓋度發(fā)生較大變化的重要因素。圖4展示了研究區(qū)以10 m為間隔高程范圍內的平均植被覆蓋度及各高程帶平均植被覆蓋度變化。1993、2006、2019年植被覆蓋度隨高程的升高總體呈逐漸增加的趨勢，并具有明顯的高程分帶特征。300～600 m范圍內隨高程的升高植被覆蓋度迅速增加，且植被覆蓋度較低；超過600 m范圍植被覆蓋度變化緩慢且復雜，呈波浪式上升的變化趨勢，主要為較高植被覆蓋度等級以上，尤其是高程>1 500 m基本只分布高植被覆蓋度。

圖4 1993—2019年植被覆蓋度隨高程變化特征及各高程帶平均植被覆蓋度變化率

從圖4(b)可知，在各高程帶(表1)上，1993—2019年植被覆蓋度總體呈逐年增加的變化特征。其中植被覆蓋度在300～500 m高程帶上增加的幅度最大；1993—2006年，植被覆蓋度在1 300～1 700 m高程帶減少，在其余高程帶增加； 2006—2019年，植被覆蓋度在500～1 100 m高程帶減少，且減少值隨高程的增高而增大，但在其余高程帶增加。

3.3.2 坡向響應

不同坡向的降雨、太陽輻射等因子強度不同，因而植被生長所需的水、熱、光照等條件不同，導致植被覆蓋度在各坡向上呈現一定的差異[25-27]。圖5展示了研究區(qū)不同坡向上的植被覆蓋度的分布特征。植被覆蓋度分布在各坡向上無明顯差異，3個時相植被覆蓋度最大值和最小值差值最大僅為0.07。1993、2006、2019年植被覆蓋度在東北、北、東面最高(分別為0.80、0.81、0.80)，最小值(分別為0.73、0.78、0.76)分別出現在西北面、東南面、北面。1993—2019年植被覆蓋度在各坡向上整體上呈增加的趨勢，2006年植被覆蓋度總體要高于1993、2019年；2019年植被覆蓋度在東北、東面低于1993年，在其他坡向上均高于1993年。

3.3.3 坡度響應

圖6展示研究區(qū)以1°為間隔坡度范圍內的平均植被覆蓋度及各坡度帶平均植被覆蓋度變化。1993、2006、2019年植被覆蓋度隨坡度的增大呈先增加后減小的趨勢，以35°為分界點。在0°～15°范圍內3期植被覆蓋度間差異不顯著，植被覆蓋度較低，且逐年降低；超過>15°范圍內，植被覆蓋度差異隨坡度的增大逐漸顯著，植被覆蓋度較高，總體呈逐年增加的變化趨勢，特別是1993—2006年植被覆蓋度在>35°范圍上增加的幅度最大，植被覆蓋得到很大的改善。

圖6 1993—2019年植被覆蓋度隨坡度變化特征及各坡帶平均植被覆蓋度變化

4 討論與結論

龍脊梯田以林地、耕地為主要土地利用類型，兩者占地總面積達89%以上，植被覆蓋度變化與人類活動密切相關，1993—2019年各級別植被覆蓋度時空分布特征具有明顯差異，較高、高植被覆蓋度主要分布在林地區(qū)域；中等植被覆蓋度主要分布在以梯田為主的耕地區(qū)域，中等植被覆蓋度是農田農作物生長狀況的鼎盛時期[28]，大面積修造梯田不僅改變地貌特征，對區(qū)域植被覆蓋也產生較大影響；低、較低植被覆蓋度面積最小，主要分布在建設用地區(qū)域，該區(qū)域是旅游開發(fā)的主要場地，植被生長易受人類活動的影響。1993—2006年，植被覆蓋度發(fā)生較大變化，中等、較高、高級別植被覆蓋度面積分別呈減小、減小、增加趨勢，這是由于從1998年政府為了龍脊景區(qū)旅游業(yè)的快速發(fā)展加大資金投入，并于2001年成立景區(qū)的經營管理部門，更加重視自然景觀的建設，加強梯田周邊生態(tài)環(huán)境的修復和保護[29]，實施局部退耕還林措施，使中等、較高植被覆蓋度向更高級別轉換；近幾年隨著區(qū)域旅游業(yè)的快速發(fā)展，游客日益增多，住房及交通道路等建設用地范圍逐漸擴大，造成生態(tài)植被的破環(huán)，導致2019年局部區(qū)域植被覆蓋度有所減小。

不同地形條件下，植被覆蓋度的空間分布格局因水、熱、人類活動等因素不同而存在差異[13]。龍脊梯田植被覆蓋度在各坡向上差異不顯著，但整體隨高程的上升呈增加的趨勢，由于較低高程帶300～1 100 m主要為區(qū)域旅游范圍，梯田分布較多，人類活動頻繁，植被生長狀態(tài)較差；而高程>1 100 m的高山區(qū)域主要分布水源林，溝壑縱橫，土質疏松，植被生長所需的水分、養(yǎng)分條件充足，且人類活動較少，因此，植被覆蓋較高，植被生長狀態(tài)良好。坡度是影響物質和能量在巖石圈的運移與轉換主要因素之一，不同坡度上環(huán)境具有的物質和能量不同，導致植被覆蓋度在各坡度上存在差異[30]。區(qū)域植被覆蓋度隨坡度的增大呈先增加后減小的趨勢，地勢比較平坦區(qū)域是人類活動及生產建設的主要場所，植被受人類活動影響較大，植被覆蓋度較低；隨坡度增大，人類活動逐漸變小，土壤水分及養(yǎng)分能滿足植被生長的需要，植被數量增多，植被覆蓋度增大；而坡度過于陡峭的地勢不利于水土保持，養(yǎng)分容易流失，進而影響植被對水分及養(yǎng)分的吸收，植被覆蓋度有一定程度的降低[16,31-32]。

1993—2019年龍脊梯田植被覆蓋度分布特征主要受人類活動、海拔、坡度等因素的影響，人類通過人工林種植擴大林地面積，提高區(qū)域植被覆蓋度，但人類活動范圍也在變大，低、較低植被覆蓋度不斷擴大。因此，在維持區(qū)域旅游業(yè)快速發(fā)展的同時，應該注重生態(tài)環(huán)境保護，重點治理區(qū)域梯田分布較廣、人類活動較大的低海拔及地勢平坦區(qū)域的生態(tài)環(huán)境，實施局部封山育林的措施，促進植被生長，加大生態(tài)環(huán)境的建設及管理，實現可持續(xù)發(fā)展。

本研究基于像元二分模型反演龍脊梯田不同時期的植被覆蓋度，并分析植被覆蓋度空間分布及其對地形驅動因子的響應特征，結論如下：

(1)龍脊梯田以林地、耕地為主要土地利用類型，1993—2006年平均植被覆蓋度增加，而2006—2019年平均植被覆蓋度呈略減趨勢。植被覆蓋度等級以較高植被覆蓋度和高植被覆蓋度為主，兩者占區(qū)域總面積高達87%以上，主要分布在林地區(qū)域。

(2)1993—2019年低、較低、中等植被覆蓋度發(fā)生等級變化主要向較高植被覆蓋度轉移，較高植被覆蓋度發(fā)生等級變化主要向高植被覆蓋度轉移，高植被覆蓋度有83.27%植被面積未發(fā)生等級變化。

(3)植被覆蓋度對不同地形因子具有不同響應特征。植被覆蓋度在各坡向無明顯差異，但北、東北、東方向的植被覆蓋度較西北、東南方向高。植被覆蓋度在不同高程、坡度等級區(qū)間呈顯著差異，植被覆蓋度整體隨高程增加呈波浪式上升的趨勢，低植被覆蓋度主要出現在人類活動較大、農業(yè)用地較多的300～600 m低海拔帶區(qū)域，高植被覆蓋度主要出現在>1 500 m高程范圍內；植被覆蓋度隨坡度的增加而呈先增加后減小趨勢，在0°～8°坡度范圍內最低，在25°～35°坡度上最高。