植被恢復對洪雅縣近15年景觀格局的影響

王 鵬，李賢偉，趙安玖，黃從德，范 川，賴家明

(四川農業大學林學院，長江上游生態林業工程四川省重點實驗室，雅安 625014)

植被恢復對洪雅縣近15年景觀格局的影響

王 鵬，李賢偉*，趙安玖，黃從德，范 川，賴家明

(四川農業大學林學院，長江上游生態林業工程四川省重點實驗室，雅安 625014)

退耕還林和天然林保護等重大林業生態工程的實施對我國西部地區生態環境建設起到了積極的推動作用。植被的重建、恢復和保護，不僅改變了景觀的類型組成，也強烈地影響著景觀的空間構造。以四川盆地山地丘陵區邊緣的四川省洪雅縣為研究對象，采用景觀生態學原理和方法，對工程實施前后15年來(1994—2009)洪雅縣景觀格局的變化進行分析，探討了植被恢復對景觀格局的影響。結果顯示，林地增加、耕地和草地減少，是15年間洪雅縣景觀變化的主要特點。耕地的減少和林地的增加主要集中于人為活動影響更為強烈的海拔1000 m以下區域，反映了退耕還林對景觀類型組成改變的巨大推動作用；海拔1500 m以上區域，更多地受到天然林保護和自然因素影響，面積變化相對穩定，景觀格局空間構造方面的變化更為顯著；海拔1000—1500 m區域，兩種工程實施的影響交錯，過程相對復雜。坡度≤5°的平緩地中，人為活動影響強烈，林地增加明顯；而在坡度25°附近的坡地中，耕地較多，退耕還林潛力較大，耕地—林地轉換相對其他坡度區段更為明顯。總體上，林地向低海拔、緩坡度區域擴張，斑塊間空隙還林面積較多，有連片發展趨勢；而耕地在面積減少的同時，向低海拔、緩坡度區域收縮，斑塊規模更小，格局更加破碎。

植被恢復；退耕還林；天然林保護；洪雅縣；景觀變化

景觀生態學強調景觀格局與生態過程在不同尺度上的相互作用[1]，景觀格局變化受到來自生物和非生物因素的共同影響，反映景觀要素在空間上的分布變化和時間上的變化趨勢，并影響與之相對應的生態過程[2]。土地利用/覆蓋變化是景觀格局變化的重要決定因素，其對景觀的結構、功能以及變化的影響已成為景觀生態學研究中的熱點問題[3]。植被恢復作為土地利用/覆蓋變化的一種表現，對于改善生態系統結構和功能，維持生態系統平衡起到至關重要的作用，也是恢復生態學研究和生態環境建設的重點領域[4- 5]。

自1998年以來，國家先后實施了包括天然林保護、退耕還林在內的六大重點林業生態工程，已經取得了階段性的成果。四川盆地山地丘陵區地處長江上游，是我國實施生態環境保護和治理的重點區域，區內坡地眾多，長期以來土地利用方式粗放，濫砍濫伐，毀林開荒，陡坡開墾，導致植被退化、水土流失、土壤退化嚴重[6]。工程實施十多年來，生態退耕以及植被恢復和保護成效顯著，森林覆蓋率穩步提高。林業生態工程的實施，植被的重建、恢復和保護，包括自然和人類活動在內的驅動因子對景觀格局和生態過程的影響積極且深遠[7]。如何從景觀尺度評價工程實施背景下植被恢復對景觀格局的影響，對于工程后續政策的制定與實施，具有重要的現實意義。洪雅縣地處四川盆地向青藏高原過渡地帶，位于“華西雨屏區”核心區域，區內雨量充沛，氣候適宜，植被恢復自然條件良好，是全國退耕還林示范縣和先進縣。同時，洪雅縣海拔跨度大，退耕類型復雜，植被恢復過程多樣，且受氣象及地形條件影響，遙感影像獲取及解譯難度較大，研究相對滯后。本文選擇四川省洪雅縣為研究區域，采用景觀生態學原理和方法，借助景觀格局指數，對國家重大林業生態工程實施前后15年來(1994—2009年)洪雅縣景觀格局進行分析，獲得景觀組成變化的主要特點，以及植被恢復對景觀格局的影響，期望對該地區后續工程實施和生態環境建設提供科學依據。

1 研究區概況

四川省洪雅縣位于四川盆地西南緣(102°49′—103°32′ E，29°24′—30°00′ N)，地處四川盆地與涼山山原過渡地帶，地勢由東北丘陵平壩向西南高山峽谷逐漸抬升，以山地、丘陵地貌為主。研究區總面積1955.07 km2，人口約34 萬，以農業人口為主。境內氣候溫和，雨量充沛，年平均氣溫16.8 ℃，年平均降水量1493.8 mm，屬中亞熱帶濕潤氣候。土壤類型以山地黃壤和紫色土為主，地帶性植被類型為亞熱帶常綠闊葉林。洪雅縣生態環境良好，森林動植物資源豐富，有四川省最大的國有林場——洪雅林場。自1998年先后開始實施天然林保護工程和退耕還林工程(試點)以來，洪雅縣植被恢復成效顯著，森林覆蓋率逐步上升，主要造林樹種有柳杉(Cryptomeriafortunei)、杉木(Cunninghamialanceolata)、光皮樺(Betulaluminifera)、雜交竹(Hybridbamboo)、苦竹(Sinobambusatootsik)、慈竹(Neosinocalamusaffinis)等。

2 數據處理與分析方法

2.1 數據來源與處理

研究采用1994年6月26日和2009年6月3日成像的兩期Landsat TM遙感影像作為主要數據源，其他數據包括研究區森林資源分布圖、林相圖、1∶5萬地形圖以及由地形圖數字化形成的10 m分辨率數字高程模型(DEM)。TM影像采用最大似然法進行監督分類，經類型歸并后，獲得研究區1994年和2009年兩期景觀類型圖(圖1)，包括耕地、林地、草地、水體、建筑用地、裸地六類。其中，耕地包括旱地和水田，林地包括針葉林、闊葉林、灌木林、和竹林，裸地既包括高海拔極端氣候條件下自然形成的裸露地塊，也包括林場采伐經營后的采伐跡地。

圖1 研究區1994年和2009年景觀類型圖Fig.1 Landscape maps of the study area in 1994 and 2009

為盡量減少由影像解譯引起的不確定性，采用ERDAS中的ATCOR3模塊分別對兩期TM影像做輻射和地形校正處理，減輕由大氣和地形引起的輻射誤差，獲取近似地表反射率的待解譯影像；而在解譯過程中，對兩期影像采用相同的處理方式，如相同的最小制圖單元，減少景觀格局分析中的不確定性[8- 9]，使景觀格局的對比具有更高的可比性。經野外實地驗證，兩期影像總體解譯精度均在85%以上。

2.2 海拔和坡度梯度的劃分

地形地貌對山區景觀組分的分布有決定性的影響。為研究洪雅縣各景觀類型在空間上的分布變化特征，從研究區DEM數據中獲取海拔和坡度圖，并采用等步長的方式，在像元尺度上將海拔和坡度梯度分別劃分為7個和10個區段(表1)，并與景觀類型圖疊加統計，獲得各區段景觀類型的面積和變化信息。

表1 海拔和坡度分區表

2.3 景觀格局分析方法

景觀格局指數作為景觀格局分析常用的方法，能夠高度概況景觀鑲嵌體的組成和空間構造方面的特征，為景觀生態學研究提供簡單而有效的定量信息[10]。但多數景觀格局指數之間存在高度的相關性[11- 12]，并且對于指數的解釋需要根據不同的格局特征以及指數隨格局的變化規律來進行，特別是在格局指數與類型組分之間經常存在的非線性關系[13- 15]，給指數的解釋帶來困難[16]。基于此，本文主要選擇幾個描述格局基本特征的指數，從斑塊面積、斑塊數目、斑塊形狀、斑塊周長面積關系方面，對景觀格局的變化及其生態學意義進行闡述，包括平均斑塊面積、斑塊密度、最大斑塊指數、邊緣密度、平均形狀指數、平均周長-面積比等6個指數，并采用Fragstats 3.4[17]軟件從類型和景觀水平分別進行計算，指數計算方式和意義請參閱軟件說明。

3 結果與分析

3.1 景觀格局總體特征

圖2 1994年和2009年景觀類型面積比例及轉移 Fig.2 Area percentage of each landscape type in 1994 and 2009, and gains and losses as a percentage of study area from 1994 to 2009 柱端數字為各類型轉移面積占總轉入和總轉出的比例; 耕地 farmland；林地 forestland；草地 grassland；水體 water；建筑用地 building land；裸地 bare land

在1994年和2009年這兩個時期，耕地和林地均為洪雅縣主要的景觀類型，二者的面積之和占景觀總面積的85%以上(圖2)。林地作為最大的景觀類型，在整個景觀中占據支配地位，面積比例由60.58%上升到66.30%，而耕地類型則由27.39%減少到22.27%，草地作為第三大類型同樣呈減少趨勢，面積比例由7.32%減少到5.81%，其他類型變化相對較小。從1994年到2009年，洪雅縣有近1/3的景觀發生了類型轉變，其中耕地和林地轉出最多，分別占景觀總面積的12.85%和10.49%，其次為草地(6.14%)，這3種類型發生的變化占洪雅縣景觀變化面積的90%以上，是景觀變化的主體類型。其中，林地是轉入面積最多的類型，占轉入總面積的50%以上，其次為耕地(23.97%)和草地(14.35%)。因而，耕地、草地減少，林地增加，是1994—2009年間洪雅縣主要的景觀變化特征。

從1994年到2009年，洪雅縣景觀整體平均斑塊面積略有增大，由5.57 hm2增加到6.08 hm2；斑塊密度則由17.96 個/km2降低到16.45 個/km2；最大斑塊指數由47.31%增大到56.27%；邊緣密度由116.17 m/hm2減少到105.53 m/hm2。這些指數的變化表明景觀總體似乎具有這樣一種變化趨勢，即景觀中小的斑塊聚合成為更大的斑塊、斑塊間的邊界減少、景觀由多個小斑塊分散分布向集中聚合的方式轉變、大斑塊在景觀中的比重更大，景觀總體有逆破碎化的變化趨勢。

3.2 耕地、林地空間分布變化

對不同海拔和坡度梯度上景觀總面積、耕地和林地面積，以及耕地和林地轉入和轉出面積進行統計，得到耕地和林地在不同海拔和坡度區段上的分布及其變化，可以反映景觀變化主體類型在空間上的分布變化特征。從統計結果得出(圖3)，洪雅縣85%以上的林地分布在海拔500—2500 m的區域，并且隨著海拔升高，在2500 m以下，林地在各區段內的比例不斷增加，而在2500 m以上區域，接近森林分布的垂直邊界，比例出現下降；林地在坡度梯度上主要分布在gt;10°的區域)，且隨著坡度的增大，林地在各坡度區段內的比例不斷升高；耕地則主要分布在海拔1500 m以下區域)，且隨著海拔升高，比例不斷下降；而在坡度梯度上，85%以上的耕地分布在≤25°的區域，并且坡度越大，耕地在坡度區段內的比例越低。總體上，洪雅縣林地主要分布在中高海拔的陡坡區域，而耕地則主要存在于低海拔的平地、緩坡區域。

從1994年到2009年，林地在海拔2500 m以下區域均表現為增加的趨勢，而在海拔2500 m以上呈減少趨勢。這其中，林地的凈增加以海拔500—1000 m的區段最多，占所有區段林地凈增加總面積的58.50%，是林地凈增加的主要區域；其次為500 m以下區段，比例達到凈增加總面積的18.72%；而在海拔1000 m以上區域，林地的面積變化則相對較小。在林地凈變化率方面，海拔500 m以下區域林地增加最為明顯(圖4)，達到了87.66%；其次為500—1000 m區域，增加了29.40%；而在1500—2500 m區域，變化率相對較小；海拔2500 m以上區域，隨著海拔升高，林地的減少趨勢更為明顯。林地在各坡度梯度上均表現為增加趨勢(圖3)，除≤5°的平緩地增加較多外，在20—25°區段存在一個面積增加的峰值。而在凈變化率上(圖4)，林地以≤30°區域增加相對明顯，除≤5°的平緩地之外，在20—25°區段也存在一個峰值，但相對于海拔梯度來說，差異相對較小。

圖3 1994年和2009年林地、耕地面積比例隨海拔和坡度梯度分布對比Fig.3 Comparison of forest and farmland distributions along elevation and slope gradients in 1994 and 2009

從1994年到2009年，耕地在所有的海拔區段均呈減少趨勢(圖3)。其中以海拔500—1000 m區段減少最多，占耕地凈減少總面積的近60%，是耕地減少的主要分布區段；其次為500 m以下區域，耕地減少面積占凈減少總面積的20.47%；其他區段耕地面積變化相對較小。在耕地凈變化率方面(圖4)，低海拔的500 m以下區域減少的比例最少，且總體上海拔越高，減少愈加明顯。但海拔1000—1500 m的區域比較特殊，這一區域耕地的轉入和轉出大體相抵，與林地在這一區段的特征相對應，均沒有表現出明顯的變化，這可能與退耕還林等林業工程實施之前，該區域存在過毀林開荒有關。耕地在各坡度區段上也呈減少趨勢(圖3)，以≤30°的區域減少較多，占耕地凈減少總面積的88.28%。在凈變化率方面，坡度20—25°區段存在一個峰值(圖4)，耕地的減少最為明顯。

圖4 1994—2009年耕地、林地隨海拔和坡度梯度凈變化率Fig.4 Net change rates of forest and farmland along elevation and slope gradients from 1994 to 2009

由此，從1994年到2009年，林地在不同海拔梯度上的變化趨勢可以分為3種類型：海拔1000 m以下，是人口和經濟活動分布的主要區域，受人為活動影響，如退耕還林、宜林荒山荒地造林等，林地增加明顯，且低海拔的500 m以下區域，受到的影響更為強烈；海拔1000—2500 m區域，人口分布相對較少，受人為活動影響較弱，盡管也受到諸如林場經營采伐活動的影響，但總體上這一區域表現為近自然的森林恢復，林地變化較為穩定；海拔2500 m以上區域，人口稀少，主要受到自然環境影響，林地中灌木林有向草地轉化的趨勢，但同時受到影像成像時植被的物候影響，這一結論有待進一步驗證。耕地在低海拔的1000 m以下區域與林地的變化相對應，耕地的減少對應著林地的增加，這與退耕還林的實施密切相關。而在高海拔的1000 m以上區域，耕地面積較少，林地的變化更多地受到天然林保護、荒山造林，以及自然環境影響，耕地與林地之間的轉換關系不明顯。

在坡度梯度上，林地的增加除≤5°的平緩地面積較多，受人為影響強烈外，在5—25°區域隨著坡度升高，增加越為明顯。一方面，中等坡度的25°附近面積較多，造林潛力大；另一方面，與gt;25°的陡坡區域相比，受到人為活動影響更大，退耕還林在其中起到了關鍵作用，與耕地在25°附近的減少率峰值相對應。gt;25°的陡坡區域林地的增加隨著坡度升高逐漸減少，盡管退耕還林工程以gt;25°的陡坡耕地為實施的重點區域，但這一區域原有的林地比例較高，因而林地并沒有明顯的增加趨勢。

3.3 類型水平格局特征

在類型水平，景觀中各類型的格局變化特征不相一致。從1994年到2009年，耕地類型的平均斑塊面積有所減小，由4.94 km2減少到3.71 km2(圖5)，斑塊規模介于水體和建筑用地之間，遠小于林地；而林地則一直具有最大的斑塊規模，并且平均斑塊面積從18.80 km2增加到26.13 km2；草地的平均斑塊面積相對較小，且變化不大。從斑塊密度上來看(圖5)，相對于其他類型，耕地的斑塊密度一直處于較高的水平，并且由5.50 個/km2增加到6.00 個/km2；而林地的斑塊密度在三大類型當中最低，且由3.22 個/km2降低到2.54 個/km2，與耕地的變化特征相反；草地的斑塊密度降低明顯，由6.09 個/km2降低到4.61 個/km2。最大斑塊指數中(圖5)，耕地和林地也表現出相反的變化特征，耕地中最大的斑塊受到分割，面積減少，而林地中的最大斑塊得到補充，面積增大，且增加顯著。在描述邊界特征的指數邊緣密度中(圖5)，耕地和林地一直具有較高的邊緣密度，其次為草地，表現為具有較多的邊界。從平均形狀指數來看(圖5)，林地較耕地具有更加復雜的形狀，而在平均周長-面積比中(圖5)，耕地則具有更大的周長-面積比。

圖5 1994年和2009年類型水平格局指數Fig.5 Class level landscape pattern indices in both 1994 and 2009

總體上，對于耕地類型來說，受退耕還林工程實施的影響，類型總面積減少，類型中大的斑塊受到其他類型的分割，平均斑塊面積減小，斑塊數目增多，大斑塊在景觀中的比例被削弱，表現出破碎化的傾向；而林地類型則與耕地相反，受退耕還林和天然林保護工程實施的影響，斑塊間空隙得到補充，總面積增大，斑塊規模增大，斑塊數目減少，大斑塊在景觀中的影響力得到增強，森林總體呈現連片式發展趨勢，連通性增大，有明顯的逆破碎化傾向。

4 討論

本文采用兩期TM遙感影像，對洪雅縣近15年來植被恢復對景觀格局的影響進行了研究，反映了重大林業工程的實施前后景觀格局變化的主要特點。退耕還林工程實施引起的耕地減少、林地增加主要集中于海拔1000 m以下區域，且對景觀類型組成的影響強于格局的空間構造。而天然林保護工程實施對景觀格局的影響主要針對海拔1500 m以上的林地，林地作為這一區域景觀的主體，面積變化較少，但空間構造方面卻發生了顯著的變化，這一區域林地的斑塊規模更大，邊界更少，景觀的完整性得到恢復和保護，天然林保護對景觀格局空間構造方面的影響更為顯著。海拔1000—1500 m區域，雖然耕地和林地的面積均沒有表現出明顯的變化趨勢，但仍然存在較多的轉化。一方面，在工程實施之前，這一區域可能存在毀林開荒等森林砍伐活動，并且瓦屋山水電站的建設也會占用林地，導致林地的轉出；另一方面，退耕還林、宜林荒山荒地造林、天然林保護以及林場經營活動，也會導致林地的轉入。因而，海拔1000—1500 m區域實際上是多種因素共同作用的結果。海拔1000 m以下，景觀受到退耕還林、宜林荒山荒地造林等因素的驅動，植被的恢復主要是以重建的方式對格局的類型組成帶來影響，政策性的人為活動的影響更為強烈；而海拔1500 m以上區域，景觀主要受到天然林保護等因素的驅動，植被的恢復更多的是以保護、維持的方式對格局的空間構造帶來影響，自然因素的影響相對更為明顯；海拔1000—1500 m區域作為兩種影響的交錯區域，退耕還林和天然林保護的影響并存，變化過程相對復雜。

盡管研究的區域、尺度，以及數據來源不同，但耕地、草地與林地之間緊密的轉換關系，在近幾年的相關研究中都有所報道[18- 21]，退耕還林和天然林保護工程的實施在其中起到了巨大的推動作用。工程的實施不僅提高了研究區森林覆蓋率，減少了水土流失，提高了森林的水源涵養能力，也使森林景觀的格局發生了顯著的變化，景觀的破碎化得到緩解和修復，給生物多樣性的維持和保護也帶來了積極的影響。如果將植被的重建作為植被恢復的第一階段，那么植被的保護和維持、格局的優化將是低海拔區域退耕還林之后植被恢復的重點。低海拔區域人口集中，社會經濟需要持續發展，人為干擾影響巨大，如何在經濟發展與生態保護之間取得平衡，實現退得住、保得住，提高人工林穩定性，使格局進一步優化，生態系統服務功能逐步提升，將是今后退耕還林等生態環境建設實現可持續發展的重要方向。而在高海拔區域，在排除各種人為干擾之后，格局趨于優化，以原有植被為本底的景觀自我維持能力得到提升，后續的植被恢復應當仍然以保護和維持為主，促進植被向頂級類型演替，逐步提高抵御外界干擾和脅迫的能力，實現生態的完全修復。

本文采用幾個描述景觀格局基本特征的指數對格局進行分析，反映了不同類型的格局變化趨勢，但對指數結果的解釋仍然存在困難，特別是幾個描述斑塊邊界特征的指數。耕地由于面積減少，邊界總數量降低，并且格局中面積減少對邊界數量的影響強于格局的空間構造，因為斑塊數量的增加和更大的周長-面積比也可能會導致更多的邊界；而林地由于面積得到補充，小的斑塊被更大的斑塊替代，斑塊聚合程度更高，邊界數量減少，格局的空間構造對邊界數量的影響強于面積的增加。但景觀格局分析當中，格局的組成成分和空間構造之間關系比較多樣，且多為非線性關系，在運用格局指數進行景觀格局研究時難以將二者厘清。有必要在以后的景觀格局分析中，將不同的變化過程進行區分，將指數的解釋與特定的變化過程結合，研究指數對不同過程的響應，可能會對景觀格局的分析帶來積極的影響。

5 結論

(1)從1994年到2009年，洪雅縣植被恢復效果顯著，林地面積比例由60.58%上升到66.30%，森林覆蓋率上升了5.72個百分點，耕地、草地則分別減少了5.12和1.51個百分點。耕地、草地減少，林地增加，這3種類型的變化是洪雅縣景觀變化的主要特點。

(2)重大林業工程的實施不僅帶來了植被的重建與恢復，也給景觀的格局帶來了積極的影響。從1994年到2009年，洪雅縣景觀總體有逆破碎化的趨勢。其中林地作為景觀的主體，空間分布向低海拔、緩坡度方向擴展，格局總體呈連片式發展趨勢；與林地的格局變化相反，耕地的空間分布向低海拔、緩坡度方向收縮，格局總體更加趨于破碎。低海拔的1000 m以下區域，景觀的變化以類型轉換為主，耕地—林地轉化強度更大；而高海拔的1500 m以上區域，面積變化較小，景觀格局的空間構造方面的變化更為顯著；海拔1000—1500 m區域，退耕還林和天然林保護的影響交錯，多種因素共同作用，變化過程相對復雜。坡度≤5°的平緩地面積較多，人為活動影響強烈，林地增加明顯；而中等坡度的25°附近區域坡耕地較多，退耕還林的潛力較大，工程實施對耕地—林地轉化的推動更為顯著。

(3)植被的恢復不僅改變了景觀中類型的組成，也給景觀的空間構造帶來影響。本文選擇了幾個描述格局基本特征的指數，可以說明格局中存在的基本變化趨勢，但景觀中存在著多樣的土地利用變化過程，不同的區域其變化過程或階段可能不同。發展新的格局分析方法，將景觀格局分析與變化過程結合，研究指數對不同過程的響應，可能會使指數的生態學意義更加明確。

致謝：感謝洪雅縣林業局在數據收集中給予的大力支持；感謝同門陳浩、王鈞、馮帥、賴元長在數據處理中提供的幫助。

[1] Turner M G. Landscape ecology: What is the state of the science?. Annual Review of Ecology, Evolution, and Systematics, 2005, 36(1): 319- 344.

[2] Turner M G. Landscape ecology: The effect of pattern on process. Annual Review of Ecology, and Systematics, 1989, 20(1): 171- 197.

[3] Wu J G, Hobbs R. Key issues and research priorities in landscape ecology: An idiosyncratic synthesis. Landscape Ecology, 2002, 17(4): 355- 365.

[4] Peng S L, Lu H F. Some key points of restoration ecology. Acta Ecologica Sinica, 2003, 23(7): 1249- 1257.

[5] Ren H, Peng S L, Lu H F. The restoration of degraded ecosystems and restoration ecology. Acta Ecologica Sinica, 2004, 24(8): 1756- 1764.

[6] Li X W, Zhang J, Hu T X, Luo C D. Theory of the Conversion of Farmland to Forest and Application of Forest-Grass Mode. Beijing: Science Press, 2009: 45- 51.

[7] Qiu Y, Zhang Y, Han J, Wang J, Meng Q H. Spatiotemporal pattern of ecological cropland conversion and vegetation succession: A review. Chinese Journal of Ecology, 2008, 27(11): 2002- 2009.

[8] Shao G F, Wu J G. On the accuracy of landscape pattern analysis using remote sensing data. Landscape Ecology, 2008, 23(5): 505- 511.

[9] Lechner A M, Langford W T, Bekessy S A, Jones S D. Are landscape ecologists addressing uncertainty in their remote sensing data? Landscape Ecology, 2012, 27(9): 1249- 1261.

[10] Li H B, Wu J G. Application of spatial pattern analysis in landscape Ecology // Wu J, ed. Lectures in Modern Ecology (III): Advances and Key Topics. Beijing: Higher Education Press, 2007: 228- 250.

[11] Li H B, Reynolds J F. A simulation experiment to quantify spatial heterogeneity in categorical maps. Ecology, 1994, 75(8): 2446- 2455.

[12] Riitters K H, O′Neill R V, Hunsaker C T, Wickham J D, Yankee D H, Timmins S P, Jones K B, Jackson B L. A factor analysis of landscape pattern and structure metrics. Landscape Ecology, 1995, 10(1): 23- 39.

[13] Gustafson E J, Parker G R. Relationships between landcover proportion and indices of landscape spatial pattern. Landscape Ecology, 1992, 7(2): 101- 110.

[14] Hargis C D, Bissonette J A, David J L. The behavior of landscape metrics commonly used in the study of habitat fragmentation. Landscape Ecology, 1998, 13(3): 167- 186.

[15] Long J A, Nelson T A, Wulder M A. Characterizing forest fragmentation: Distinguishing change in composition from configuration. Applied Geography, 2010, 30(3): 426- 435.

[16] Li H B, Wu J G. Use and misuse of landscape indices. Landscape Ecology, 2004, 19(4): 389- 399.

[17] McGarigal K, Cushman S A, Neel M C, Ene E. FRAGSTATS v3: Spatial Pattern Analysis Program for Categorical Maps. [2013- 3- 18]. http://www.umass.edu/landeco/research/fragstats/fragstats.html

[18] Zhao F, Ju H B, Huang J W. Landscape changes after conversion of cropland to forest project implementation in Xichang Area. Scientia Silvae Sinicae, 2006, 42(S1): 56- 61, 171- 171.

[19] Zhu Z Q, Liu L M, Zhang J L. Impact of grain for green project on landscape pattern in hilly loess region in Southern Ningxia: landscape evolution process assessment of Zhong-zhuang-cun small watershed in 1993—2005. Acta Ecologica Sinica, 2010, 30(1): 146- 154.

[20] Zhang G K, Lu J H, Song K S, Zhang S Q, Liang Y H. The change of landscape pattern in Zhenlai Xian, Jilin Province in recent ten years. Acta Ecologica Sinica, 2012, 32(12): 3958- 3965.

[21] Zhou D C, Zhao S Q, Zhu C. The impact of the grain for green project on the land use/cover change in the Northern Farming-pastoral Ecotone, China — A case study of Kezuohouqi County. Scientia Geographica Sinica, 2012, 22(4): 442- 449.

參考文獻：

[4] 彭少麟, 陸宏芳. 恢復生態學焦點問題. 生態學報, 2003, 23(7): 1249- 1257.

[5] 任海, 彭少麟, 陸宏芳. 退化生態系統恢復與恢復生態學. 生態學報, 2004, 24(8): 1756- 1764.

[6] 李賢偉, 張健, 胡庭興, 羅承德. 退耕還林理論基礎及林草模式的實踐應用. 北京: 科學出版社, 2009: 45- 51.

[7] 邱揚, 張英, 韓靜, 王軍, 孟慶華. 生態退耕與植被演替的時空格局. 生態學雜志, 2008, 27(11): 2002- 2009.

[10] 李哈濱, 鄔建國. 空間分析方法在景觀生態學中的應用 // 鄔建國. 現代生態學講座 (III): 學科進展與熱點論題. 北京: 高等教育出版社, 2007: 228- 250.

[18] 趙峰, 鞠洪波, 黃建文. 西昌地區實施退耕還林工程后的景觀格局變化. 林業科學, 2006, 42(S1): 56- 61,171- 171.

[19] 朱戰強, 劉黎明, 張軍連. 退耕還林對寧南黃土丘陵區景觀格局的影響——以中莊村典型小流域為例. 生態學報, 2010, 30(1): 146- 154.

[20] 張國坤, 盧京花, 宋開山, 張樹清, 梁雨華. 吉林省鎮賚縣近10年景觀格局變化. 生態學報, 2012, 32(12): 3958- 3965.

[21] 周德成, 趙淑清, 朱超. 退耕還林還草工程對中國北方農牧交錯區土地利用/覆被變化的影響——以科爾沁左翼后旗為例. 地理科學, 2012, 22(4): 442- 449.

EffectsofvegetationrestorationonlandscapepatternofHongyaCountryinrecent15years

WANG Peng, LI Xianwei*, ZHAO Anjiu, HUANG Congde, FAN Chuan, LAI Jiaming

KeyLaboratoryofForestryEcologicalEngineeringinUpperReachesofYangtzeRiverofSichuanProvince,CollegeofForestry,SichuanAgriculturalUniversity,Ya′an625014,China

Great forestry ecological engineerings, i.e. Conversion of Farmland to Forest (CFF) and Natural Forest Conservation (NFC), have played positive roles in promoting the ecological environment constructions in western China. The restoration of vegetation in these projects changed not only the composition but also the configuration of landscape pattern. In this paper, the landscape pattern of Hongya Country in 1994 and 2009 were compared to examine the effects of vegetation restoration on landscape pattern in the mountain-hilly region of the edge of Sichuan Basin, China. Landscape patterns at the two periods were both derived from Landsat TM images using a supervised classification procedure. The area of vegetation types were compared to examine the spatial distribution changes at different elevation and slope zones, and six landscape pattern indices were calculated to quantify the changes of landscape pattern at both landscape and class level. The elevation and slope zones were classified with the step sizes of 500m and 5° respectively, and the indices were calculated in Fragstats 3.4.

From 1994 to 2009, the total area of forest in Hongya increased by 5.72% from 60.58% to 66.30%, and the farmland and grassland decreased by 5.12% and 1.51% respectively. The increase of forestland and decrease of farmland and grassland are the major effects of vegetation restoration. Most of the conversions from farmland to forestland occurred under the elevation of 1000m, indicating a remarkable promotion of CFF to the change of landscape composition. Furthermore, as the dominant type, forest occupies more than 80% of total area in the region higher than 1500m. The vegetation restoration here was mainly influenced by the NFC and natural conditions, exhibiting little change of amount, but more change of landscape configuration. In region between the elevations of 1000m and 1500m, the amounts of both forestland and farmland were relatively stable, but remained some farm-forest conversions, and the landscape changes here were more complex, partly for the overlapping influences of CFF and NFC. At the slope aspect, forest in the region ≤ 5° increased obviously for the intense impact of human activities; moreover, in the region of near 25° forest also exhibited a little increasing trend, indicating more potential conversions from farmland to forest in the implement of CFF.

The overall landscape pattern change of Hongya exhibited a little defragmentation trend from 1994 to 2009. For farmland, influenced by CFF, patches were cut apart and converted to forest, resulted in more patches and smaller mean patch size, and the patches were more fragmented. While, forest was influenced by both the CFF and NFC, the gaps between forest patches were more filled, resulted in fewer patches and larger mean patch size, and the patches were more expanded and aggregated. In summary, as the consequences of vegetation restoration, forest was more continuously expanded to areas with lower elevation; while farmland was more shrunk, and its pattern was more fragmented.

At last, different change processes can result in different composition and configuration changes of landscape pattern. The method of traditional landscape pattern analysis cannot distinguish various processes of landscape change, in which the landscape and class level indices are indistinct and uninterpretable. There is a need to develop novel approaches to analysis the changes of landscape pattern for various change processes, especially for the interpretation of landscape indices.

vegetation restoration; conversion of farmland to forest; natural forest conservation; Hongya country; landscape change

國家“十一五”、“十二五”科技支撐計劃項目(2006BAD23B05, 2011BAC09B05)

2013- 03- 18;

2013- 08- 02

*通訊作者Corresponding author.E-mail: lxw@sicau.edu.cn

10.5846/stxb201303180447

王鵬，李賢偉，趙安玖，黃從德，范川，賴家明.植被恢復對洪雅縣近15年景觀格局的影響.生態學報,2013,33(20):6721- 6729.

Wang P, Li X W, Zhao A J, Huang C D, Fan C, Lai J M.Effects of vegetation restoration on landscape pattern of Hongya Country in recent 15 years.Acta Ecologica Sinica,2013,33(20):6721- 6729.