曾春陽，曾廣宇 ，程麗華 ，劉志斌 ，廖 源 ，黎 波

林地變化及景觀斑塊特征的演變

曾春陽1，曾廣宇2，程麗華3，劉志斌1，廖 源4，黎 波4

（1. 廣西林業(yè)勘測設計院，廣西 南寧 530011； 2. 廣西林業(yè)科學研究院，廣西 南寧 530002；3. 廣西國土測繪院， 廣西 南寧 530023; 4. 岑溪市林業(yè)局，廣西 岑溪 543200）

利用廣西岑溪市1999年、2009年兩期森林資源規(guī)劃設計調查數據，借助GIS分析工具及景觀生態(tài)學原理，分析了間隔期內林地變化及林地景觀斑塊特征演變情況。研究結果表明：在10年中，岑溪市林地面積有所增加，林地的增加尤其是有林地較大幅度的增加導致了林地整體質量的提升；研究區(qū)內林地景觀斑塊形狀指數有所下降，形狀趨于規(guī)則；斑塊分維數有所下降，形狀趨于規(guī)律；斑塊密度指數有所上升，受干擾程度有所增加。全市林地景觀格局呈現(xiàn)多樣化，破碎化程度增加，林地景觀格局的演變主要受人類活動的干擾影響。

林地；斑塊特征；景觀生態(tài)學；破碎化；岑溪市

林地是森林資源的載體和重要組成部分，是林業(yè)發(fā)展的基礎，是控制森林資源消耗、增強森林生態(tài)防護效益、實現(xiàn)國民經濟可持續(xù)發(fā)展的根本與保障[1-3]。研究林地的變化，對揭示林地利用規(guī)律、形成合理保護和利用林地的決策具有極其重要的意義[4]，而林地景觀斑塊特征的演變影響著景觀中許多空間格局和生態(tài)過程[5]，反映了林地景觀受人類活動影響的情況。本文利用廣西岑溪市1999年、2009年兩期森林資源規(guī)劃設計調查數據，借助GIS分析工具及景觀生態(tài)學原理，分析了間隔期內林地變化及林地景觀斑塊特征演變情況，對岑溪市以后的林地保護利用、景觀規(guī)劃和管理及林地資源的可持續(xù)利用具有一定的理論和實踐作用。

1 研究區(qū)概況

廣西壯族自治區(qū)岑溪市位于東經110°3′～111°22′， 北 緯 22°37′～ 23°13′， 屬 南 亞 熱 帶 濕潤季風氣候區(qū)，地處云開大山北麓東段的丘陵山區(qū)，地勢東南高，西北低，境內山脈為云開大山及勾漏山余脈延伸。全市植被類型豐富，主要樹種有馬尾松Pinus massoniana、杉木Cunninghamia Lanceolata、尾葉桉Eucalyptus urophylla、巨尾桉 Eucalyptus grandis、 紅 錐 Castanopsis Hystrix、大葉櫟Castanopsis fissa、青岡Cyclobalanopsis glauca、栲樹Castanopsis fargesii、八角Lllicium verum、肉桂Cinnamomum cassia等。

2 研究方法

2.1 數據來源與分類系統(tǒng)

本文以廣西壯族自治區(qū)岑溪市為例，研究中所采用的岑溪市1999年、2009年林地空間分布數據均來自于兩期森林資源規(guī)劃設計調查成果。根據《中華人民共和國森林法實施條例》及國家林業(yè)局2003年4月頒布的《森林資源規(guī)劃設計調查主要技術規(guī)定》，結合《廣西壯族自治區(qū)森林資源規(guī)劃設計調查主要技術規(guī)程》（2009年）及實際情況，本次研究把林地類型分為：①有林地(woodland)；②疏林地(open woodland)；③灌木林地(shrub land)；④未成林造林地(yong afforestation land)；⑤其它林地(other forestland)，其它林地中含無立木林地、宜林地、輔助生產用地及苗圃地。非林地統(tǒng)一歸為一類：⑨非林地(non forestland)。

2.2 林地變化及斑塊信息提取

將1999年及2009年兩期森林資源規(guī)劃設計調查數據統(tǒng)一到同一坐標系和投影下，把矢量數據轉化為20 m×20 m的ArcGIS柵格數據，利用GIS的圖層疊置功能求算出間隔期內林地變化信息，并以此為基礎，計算林地轉移矩陣及其流向百分比矩陣[6-9]。同時將兩期柵格數據導入Fragstats3.3中，得到林地斑塊特征有關數據。

各林地類型變化率指數可以定量描述區(qū)域林地類型變化的速度及區(qū)域林地變化中的類型差異[10]。公式表達為：

式中： K為研究時段內變化率，Ua，Ub分別為該林地類型在研究初期和末期的面積，T為研究時段。

2.3 林地景觀斑塊特征研究指標的選取

關于斑塊特征的評價指標有很多，根據研究區(qū)域內獲取數據的特點。本文采用目前相對成熟的方法計算和分析，從斑塊的破碎化特征和形狀特征對林地景觀斑塊進行破碎化評價。所選取主要指標的名稱，計算公式及生態(tài)意義詳見表1[11-14]。

3 研究結果及分析

3.1 林地變化分析

1999～2009年林地轉移數量及流向百分比矩陣詳見表2。10年間，不同林地的變化情況有所不同，岑溪市1999～2009年林地變化的主要特征為：

有林地、疏林地、其它林地面積有所增加，其中以有林地和其它林地增加面積最多，分別增加了6 547.5 hm2、2 688.4 hm2；灌木林地、未成林造林地、非林地面積有所減少，分別減少6 399.6 hm2、2 176.7 hm2、900.0 hm2。就有林地而言，有林地與各林地均存在轉入轉出的關系。有林地轉出到其它林地面積最大，轉出到非林地的面積最?。还嗄玖值丶拔闯闪衷炝值剞D入到有林地的面積較大，疏林地轉入有林地面積最小?？傮w上，凈減少的有林地主要是轉換成了其它林地，其次是灌木林地，再次是疏林地。減少的灌木林地主要轉換為有林地和其它林地，凈增的灌木林地主要來源是有林地和非林地。疏林地、未成林造林地、其它林地主要是與有林地之間轉入轉出。從各林地之間的交換量（轉入加轉出）來看，有林地、灌木林地、其它林地居多；疏林地和非林地較少。

表1 景觀生態(tài)特征指標及生態(tài)意義Table 1 Landscape ecological index and ecological significance

表2 1999～2009年岑溪市各地類轉換信息Table 2 Land-use changing information in Cenxi City from 1999 to 2009

1999～2009年，岑溪市林地總面積增加了900.0 hm2，年平均增加0.12%，而有林地面積凈增加6 547.5 hm2，年平均增加0.35%，這與退耕還林工程的實施及征占用林地補償平衡政策相關，使林地得到了有效的保護。按照Costanza估算的森林服務價值為302美元/hm2計算[15]，僅由于林地面積的增加，所增加的生態(tài)服務功能價值就達到了27.18萬美元。岑溪市林地的增加尤其是有林地的較大幅度增加導致林地整體質量有所提升。

3.2 林地景觀斑塊特征演變分析

3.2.1 林地景觀斑塊的總體特征演變分析

岑溪市土地總面積為277 122.0 hm2，其中林地面積由1999年的204 210.2 hm增加到2009年的205 110.2 hm2，林地景觀處于岑溪市景觀本底或基質的位置。研究區(qū)內1999年、2009年林地景觀分別有1 514、3 617個斑塊，各個景觀類型所擁有的斑塊總面積和總周長也極不均衡（見表3）。兩個年度中，有林地在整個林地景觀中都占有絕對優(yōu)勢，斑塊總面積分別占91.81%、94.60%，斑塊數分別占20.48%、26.51%，是研究區(qū)內的主導景觀類型?？梢?，有林地一直是整個林地景觀的基質，控制著整個景觀的物質和能量流動，其它林地類型景觀組分鑲嵌分布其中，呈小面積零星分布。灌木林地斑塊總面積分別占5.50%、2.36%，但是斑塊數所占比重兩個年度都超過了40%，說明灌木林地破碎化程度越來越大，受人工干擾的程度很深，這主要是因為岑溪市的灌木林地大部分都為經濟林樹種，屬于國家特別規(guī)定的灌木林。未成林造林地斑塊數、斑塊總面積、斑塊總周長都大幅下降，說明岑溪市在間隔期內生態(tài)環(huán)境趨于良好，未成林造林地自然成長后進入有林地。疏林地、其它林地斑塊數、斑塊總面積、斑塊總周長比重有所增加，這主要是部分林農采用“拔大毛”的方式對有林地進行采伐導致向疏林地轉化，2008年的雨雪冰凍災害及采伐后沒有及時更新造林使有林地向其它林地轉化。

表3 岑溪市林地景觀斑塊數量變化特征Table 3 Patch quantity change characteristics of forest landscape in Cenxi city

3.2.2 林地景觀斑塊的形狀及破碎化特征演變分析

斑塊的平均面積越大，景觀連通性就越好，反之就越差。由表4可見，景觀總體的平均值由1999年的134.88 hm2下降為2009年的56.71 hm2，各林地景觀類型的斑塊平均值均有所下降，這表明各林地景觀類型連通性均有所下降，但有林地還是一直保持著相對較高的連通性。各景觀類型斑塊面積的標準差反應了同一景觀類型中斑塊面積的變動程度大小。兩期數據各景觀斑塊的面積標準差大小順序分別為：有林地＞灌木林地＞其它林地＞未成林造林地＞疏林地、有林地＞未成林造林地＞其它林地＞疏林地＞灌木林地，除灌木林地外指標的大小順序變化不大，顯示了有林地景觀的斑塊面積大小差異一直非常顯著，異質性特征最為強烈，而其它林地類型景觀斑塊間面積大小差異不十分明顯。有林地景觀斑塊間面積大小的顯著差異顯示了森林植被的自然屬性，而其它幾種林地景觀類型面積差異相對較小，說明受到人類活動的深刻影響和干擾。

斑塊的形狀指數越接近于1，說明該類型景觀的斑塊形狀越接近于圓形，數值越大說明形狀越復雜。整個林地景觀類型的形狀指數都有所下降，形狀趨于規(guī)則。而因為有林地的斑塊總面積和平均面積在整個景觀中占有絕對優(yōu)勢，故其斑塊形狀較為復雜，其它景觀類型的平均面積都很小，形狀反而比有林地更加接近于圓形，形狀破碎化程度也較低。

表4 岑溪市林地景觀斑塊形狀及破碎化變化特征Table 4 Patch shape and heterogeneity change characteristics of forest landscape in Cenxi city

分維數用來測定斑塊形狀對斑塊內部生態(tài)過程的影響，如動物遷徙，物質交流等。其理論取值范圍為1.0～2.0，分維數愈接近于1，表明斑塊的自我相似度越強，形狀越有規(guī)律，斑塊受人為干擾的程度也愈大；分維數越接近于2，則表明斑塊具有愈為復雜的形狀[16]。間隔期內，各林地景觀類型的分維數都有所下降，形狀越趨于規(guī)律。受地形等環(huán)境因子控制的景觀類型分維數較大，如其它林地分別為1.48、1.28，而受人類經營活動影響較大的其它景觀類型分維數相對較小。

斑塊密度指數（亦稱孔隙度），其可以揭示出景觀基質被該類型斑塊分割的程度，孔隙度越高，表明某一類型在景觀中的分布就越廣、影響就越大，從而可以識別不同景觀類型受干擾的程度[17]。間隔期內，岑溪市各林地景觀類型的密度指數都有所上升，受干擾程度有所增加，其中有林地的密度分別為0.11塊/km2、0.15塊/km2，其一直是岑溪市的基質景觀，面積大，連通性好，因此孔隙度一直要比其它景觀類型的低，而其它景觀類型的林地斑塊小而分散，表明受人為作用影響較大。

4 結 論

研究結果表明：岑溪市1999～2009年間，林地面積有所增加，林地的增加尤其是有林地較大幅度的增加導致了林地整體質量的提升；林地景觀斑塊平均面積減小，破碎化程度增加，全市林地景觀格局呈現(xiàn)多樣化，林地景觀格局的演變主要受人類活動的干擾影響。

（1）通過林地矩陣轉移分析，岑溪市林地總面積增加了900.0 hm2，年平均增加0.12%，而有林地面積凈增加6 547.5 hm2，年平均增加0.35%，這與退耕還林工程的實施及征占用林地補償平衡政策相關，使林地得到了有效的保護。

（2）研究區(qū)內各個景觀類型的斑塊數、所擁有的斑塊總面積和總周長極不均衡，有林地一直作為主導景觀類型，其斑塊面積、周長大小差異大，一直是整個林地景觀的基質，控制著整個景觀的物質和能量流動。其它林地類型景觀組分鑲嵌分布其中，呈小面積零星分布。

（3）研究區(qū)內整個林地景觀類型的形狀指數都有所下降，形狀趨于規(guī)則，有林地斑塊形狀較為復雜，其它景觀類型的比有林地更加接近于圓形，形狀破碎化程度也較低。各林地景觀類型的分維數都有所下降，形狀越趨于規(guī)律，受地形等環(huán)境因子控制的景觀類型分維數較大，而受人類經營活動影響較大的其它景觀類型分維數相對較小。各景觀類型的密度指數都有所上升，受干擾程度有所增加，有林地一直是岑溪市的基質景觀，面積大，連通性好，而其它景觀類型的林地斑塊小而分散，受人為作用影響較大。

林地的變化及林地景觀斑塊特征的演變不僅體現(xiàn)了林業(yè)發(fā)展的趨勢，而且表征了國土狀態(tài)的演變趨勢。這個發(fā)展趨勢合理與否，直接關系到林業(yè)生產的規(guī)模和布局，也對整個國計民生起著舉足輕重的作用。通過對岑溪市林地變化和林地景觀斑塊特征演變研究，其今后的林地發(fā)展和保護策略應注意以下幾個方面的內容：林地必須以保護為主，尤其是生態(tài)公益林區(qū)內的林地，盡量減少各種不利的干擾，控制林地景觀的破碎程度；合理對林地進行開發(fā)和利用，增加林地景觀的功效性，促進林地系統(tǒng)生產力的提高；協(xié)調好利用和開發(fā)之間的關系，同時加強區(qū)域內物質循環(huán)、能量流動、生物的遷徙及保護生物多樣性的監(jiān)測。

[1] 徐新良，劉紀遠，莊大方，等. 中國林地資源時空動態(tài)特征及驅動力分析[J]. 北京林業(yè)大學學報，2004，26(1)：41 - 46.

[2] 杜國明，張樹文. 20世紀末期大興安嶺北段林地變化研究[J]. 資源科學，2006,28(3):59-65.

[3] 曾春陽，唐代生，唐嘉鍇. 森林立地指數的地統(tǒng)計學空間分析[J]. 生態(tài)學報，2010,30(13):3465-3471.

[4] 楊存建，張 洋，程曦冉，等. 基于遙感和GIS的20世紀90年代后半期四川林地動態(tài)變化[J]. 生態(tài)學報，2006,26(12):4113-4119.

[5] 封 斌，高保山，麻保林，等. 陜北榆林風沙區(qū)農田防護林結構配置與效益研究[J]. 西北林學院學報，2005,20(2):118-124.

[6] Emma D，Sara H. Ontogenetic changes in social behaviour in the forest tent caterpillar [J] . Malacosoma disstria. Behavioral Ecology and Sociobiology，2004，56(2)：177-184.

[7] Stephen L R，Bryan B. Modeling and spatial prediction of pre-settlement patterns of forest distribution using witness tree data[J]. Environmental and Ecological Statistics，2006，13(4):427-448.

[8] Hao Q Y，Zhou Y P，Wang L H，et al. Optimization models of stand structure and selective cutting cycle for large diameter trees of broadleaved forest in Changbai Mountain [J] . Journal of Forestry Research，2006，17(2):135-140.

[9] 唐寬金，鄭新奇，閆弘文，等. 基于空間相鄰的地類轉換傾向性模型的構建及應用[J]. 生態(tài)學報，2009,29(1):337-343.

[10] 劉紀遠. 中國資源環(huán)境遙感宏觀調查與動態(tài)研究[M]. 北京：中國科學技術出版社，1995.

[11] 張金屯，邱 揚，鄭鳳英. 景觀格局的數量研究方法[J]. 山地學報，2000,19(4):346-352.

[12] 傅伯杰. 黃土區(qū)農業(yè)景觀空間格局分析[J]. 生態(tài)學報，1995,15(2):336-341.

[13] Tumer M G. Spatial and temporal analysis of landscape patterns[J]. Landscape Ecology，1990,4(1):21-30.

[14] 鐘來元，粵西低山丘陵景觀演變及其驅動力研究[J]. 湖南師范大學自然科學學報，2005,28(2):85-86.

[15] 何承剛，馮 彥，楊燕平. 西雙版納林地景觀演變過程及其驅動力分析[J]. 云南地理環(huán)境研究，2008,20(5):12-14.

[16] 王永軍，李團勝，劉 康，等. 榆林地區(qū)景觀格局及其破碎化評價[J] . 資源科學，2005,27(2):161-166.

[17] 馬克明，傅伯杰. 北京東靈山地區(qū)景觀格局及破碎化評價[J]. 植物生態(tài)學報，2000,24(3):320-326.

Analysis of change and patch landscape characteristics of forestland

ZENG Chun-yang1, ZENG Guang-yu2, CHENG Li-hua3, LIU Zhi-bin1, LIAO Yuan4, LI Bo4
(1.Guangxi Forest Inventory & Planning Institute, Nanning 530011, Guangxi, China; 2. Guangxi Forest Research Institute, Nanning 530002, Guangxi, China; 3. Guangxi Land Surveying & Mapping Institute, Nanning 530023, Guangxi, China;4. Forestry Bureau of Cenxi City, Cenxi 543200, Guangxi, China)

By using the data of forest inventory planning and design in Cenxi city in 1999 & 2009, the forestland and landscape patch changes were analyzed with the support of GIS and landscape ecology principle.The results show that in 1999～2009 the total area of the forestland increased slightly, the increasing of forestland, especially greatly increasing of woodland led to the improvement of the whole forestland quality. The value of shape index of the landscape was decreasing, shape trending to be simplified. The heterogeneity of the landscape was decreasing, shape trending to be law. The density index of the landscape was increasing, the degree of interference was increasing. The forestland landscape presented diversification in study area, the fragmentation of forestland landscape was increasing,the mainly reason of change was human activities.

forestland; patch characteristics; landscape ecology; fragmentation; Cenxi City

S718.557；S757.2

A

1673-923X(2012)08-0033-04

2012-03-10

廣西區(qū)科技廳攻關項目(林科攻0992007-5)

曾春陽(1981—)，男，江西于都人，工程師，碩士，主要從事林業(yè)調查規(guī)劃設計、營造林等工作; E-mail:zchy12000@163.com

[本文編校：吳 彬]