佘宇晨，陳彩虹，陸 禹，陳冬洋

（中南林業(yè)科技大學，湖南 長沙 410004）

海南省五指山市林地景觀異質性研究

佘宇晨，陳彩虹，陸 禹，陳冬洋

（中南林業(yè)科技大學，湖南 長沙 410004）

運用海南省五指山市1998年、2003年、2007年3期遙感影像數(shù)據，2010年森林資源二類調查數(shù)據和1993年、1998年、2003年、2008年4期森林資源一類清查數(shù)據，應用景觀生態(tài)學、森林可持續(xù)經營和地統(tǒng)計學理論和方法，選取聚合度（AI）、景觀形狀指數(shù)（LSI）、邊緣密度（ED）、破碎度指數(shù)（FN）、連接度指數(shù)（CONNECT）、斑塊密度（PD）等6個景觀格局指數(shù)，利用景觀fragstats4.0軟件計算五指山市不同林地景觀類型的景觀格局指數(shù),在此基礎上對五指山市林地景觀的異質性進行綜合分析與研究。研究結果表明，近20年內五指山市的林地景觀類型分布一直較為均勻，復雜度基本較少改變; 各林地景觀類型的斑塊密度，均呈明顯的增加趨勢，形狀指數(shù)與邊緣密度變化方向一致，各斑塊間相互影響及轉化較小，景觀斑塊間類型的改變主要是通過林地與非林地間的轉化完成，景觀聚合度的年變化不大。研究成果對林地景觀動態(tài)評價、生態(tài)脆弱性分析和林地綜合管理具有較大的理論與實踐指導意義。

林地景觀異質性；遙感影像；景觀格局指數(shù)；海南省五指山市

對于“景觀異質性”這個概念國內外有較多定義,相對一致的觀點認為，“景觀異質性是指在景觀中對一個物種或高級生物組織的內在決定作用的資源（或某種形狀）在空間（時間）上的變異程度（強度）”[1]。許多學者在景觀異質產生的機理上的觀點有所分歧[2-3]。產生景觀異質性的機制正是基于熱力學原理，即在開放系統(tǒng)中，能量由當前狀態(tài)轉化為另一種狀態(tài)，伴隨著新結構的建立而使此系統(tǒng)增加了異質性[4-6]。景觀異質性起源于系統(tǒng)與系統(tǒng)要素之間的原生差異，同時也來源于外界干擾（如人類錯誤生態(tài)行為的影響）與系統(tǒng)運動的不平衡，是隨著其中某一項景觀要素出現(xiàn)的相對頻率變化而產生變化的[7]。當景觀中存在≤1種景觀要素時，可以認為該景觀是均質的，由于人為干擾或是自然演變，產生另外一種景觀要素并占據有一定比例的時候，可以認為景觀的同質性受到了破壞，產生了異質性。景觀異質性會伴隨著該要素出現(xiàn)頻率的增加而不斷增加，直到達到某一臨界閾值占據主導地位的時候，此要素的相對頻率繼續(xù)增加而景觀異質性開始下降，景觀則趨于均質化[8]。

景觀異質性使得各個景觀個體之間產生內部與邊緣的差異，這種格局會對有機體的運動、物質流、能量流等產生影響，而且對生態(tài)系統(tǒng)功能的維持有著非常重要的意義[9]。景觀異質性的研究主要有空間異質性（Spatial heterogeneity）、時間異質性（Temporal heterogeneity）、功能異質性（Functional heterogeneity）等3個方面。近年來，由于空間統(tǒng)計學的蓬勃發(fā)展，出現(xiàn)了一系列用來統(tǒng)計、度量變量在空間上的分布特征和影響程度的方法[10]。目前，景觀異質性的分析方法主要是采用景觀格局指數(shù)法。本文的研究成果對林地景觀動態(tài)評價、生態(tài)脆弱性分析和林地綜合管理具有較好的理論與實踐指導意義，對景觀格局的評價分析具有一定的現(xiàn)實意義。

1 材料與方法

1.1 研究區(qū)概況

五指山市地處海南島中南部腹地，經緯度介于 109°19′～ 109°44′E、18°38′～ 19°02′N， 東 鄰瓊中縣、西接樂東縣、南臨保亭縣、北至白沙縣。全境東西寬43.1 km，南北長45.0 km，總面積1 128.87 km2，占海南島面積的3.5%。境內地勢東、南、北三面高、西面低，地貌類型分為：中山、低山、高丘、低丘、臺地和河流谷地，主要山脈為五指山山脈，自東經阿陀嶺向西延伸，把本市分成南北兩片。大小河流120多條，年均徑流量1.3×1010m3。氣候類型屬熱帶雨林季風氣候。土壤類型多樣，既有地帶性土壤，也有非地帶性土壤，其垂直分布規(guī)律為：海拔750 m以上為山地黃壤；海拔400～750 m為赤紅壤；海拔150～400 m為磚紅壤。非地帶性土壤有水稻土，分為潴育型水稻土、淹育型水稻土、沼澤型水稻土、潛育型水稻土、滲育型水稻土等5個亞類。全市土地總面積112 886.67 hm2，其中：林地102 015.35 hm2，占土地總面積的90.37%；非林地10 871.31 hm2，占土地總面積的9.63%。全市活立木總蓄積9 694 617 m3，其中：林分蓄積9 686 463 m3，占總蓄積的99.92%；散生木蓄積2 391 m3，占總蓄積的0.02%；四旁樹蓄積5 763 m3，占總蓄積的0.06%。全市森林覆蓋率為86.44%。但是近年來隨著社會經濟發(fā)展和人為干擾程度加深，五指山市林地景觀破碎化程度增大，生態(tài)系統(tǒng)呈現(xiàn)脆弱化趨勢，生態(tài)環(huán)境保護迫在眉睫。

1.2 數(shù)據來源與處理

采用1998年、2003年和2007年3期海南省五指山市TM遙感影像數(shù)據為研究基礎，分別為1998年、2003年兩期SPOT數(shù)據（空間分辨率為10 m），以及2007年的ALOS數(shù)據（空間分辨率為2.5 m）；另外由于SPOT數(shù)據為全色波段數(shù)據，在部分地區(qū)判讀中存在困難，輔以1998年、2003年2期的TM數(shù)據（假彩色，空間分辨率為30 m）加以判讀；海南省五指山市在1993年、1998年、2003年和2008年四期森林資源連續(xù)清查資料，數(shù)據格式為 shape file 文件，包括屬性數(shù)據和空間數(shù)據； 2010年森林資源二類調查數(shù)據，以小班為單元，包括地類、林種、優(yōu)勢樹種、齡組、起源、林齡、郁閉度、可及度、自然度、蓄積量等詳細信息；2000～2013年五指山市統(tǒng)計年鑒資料。同時采用北京54坐標系配準各類數(shù)據，對3期TM遙感影像數(shù)據進行實地調查解譯，并獲得林地景觀類型的矢量數(shù)據。

圖1 固定樣地分布Fig.1 Distribution of fixation plots

1.3 景觀指標的計算

（1）聚合度（AI）

式中：gij為基于單倍法的斑塊類型像元之間的節(jié)點數(shù)；maxgij為基于單倍法的斑塊類型i像元之間的最大節(jié)點數(shù)。等于gij的實際值除以該類型最大限度聚集在一起時的gij最大值。

（2）景觀形狀指數(shù)LSI

式中：ei指類型i的邊緣長度或周長（用柵格表面數(shù)目表示），包括涉及斑塊類型i的所有景觀邊界線和背景邊緣；min為ei的最小可能值。LSI就等于相關斑塊類型的總邊緣長度除以總邊緣長度最小可能值[11]。

（3）邊緣密度ED

式中：eik為景觀中相應斑塊類型的總邊緣長度，包括涉及該斑塊類型的景觀邊界線和背景部分；A為景觀總面積。ED就等于總邊緣長度除以景觀總面積，乘以10 000，單位為：米/公頃。

（4）破碎度指數(shù)（FN）

破碎度表示景觀中斑塊被分割的破碎程度，反映了景觀空間結構的復雜性[11]。

公式如下：

式中：Ci為景觀i的破碎度，Ni為景觀i的斑塊數(shù)，Ai為景觀i的總面積。

（5）連接度指數(shù)CCONNECT

式中：cijk為與用戶指定臨界距離范圍內的，和斑塊類型i有關的斑塊j與k的連接情況；ni為景觀中斑塊類型i的斑塊數(shù)量。連通度等于某類型斑塊中所有斑塊之間的節(jié)點數(shù)目再除以所有可能的節(jié)點數(shù)目乘以100轉化為百分比[11]。

（6）斑塊密度PD

式中：N等于景觀斑塊中斑塊的總數(shù)；A等于總景觀面積（m2）；PD等于景觀斑塊中斑塊數(shù)除以總景觀面積，在乘以10 000，再乘以100（轉化100 hm2）。注：總景觀面積包括任何內部背景。單位：斑塊數(shù)/100 hm2[11]。

2 結果與分析

通過Fragstats4.0計算可以得出五指山市1998年、2003年、2007年的林地景觀指標值,結果見表1。

由表1可知，各景觀類型中闊葉混交林的形狀指數(shù)最大為57.901 7，而非林地的形狀指數(shù)最小為2.707 7。因形狀指數(shù)LSI值越大，斑塊周邊越發(fā)達，表明1998年各景觀類型中，闊葉混交林周邊最發(fā)達，非林地周邊最不發(fā)達。闊葉混交林面積最大，從而導致混交林分布較為廣泛。非林地的破碎度FN最小為0，而針葉類破碎度FN最大為14 959.64，表明非林地斑塊在給定距離范圍內沒有與它類型相同的斑塊，而針葉類則較多。連接度指數(shù)中，非林地連接度最大為30，竹林類的連接度最小為0.418 1，說明非林地類與周邊斑塊的連通性最大，因為本文研究的主要是林地景觀，因該數(shù)據處理時將非林地（水域、建設用地、未利用地等）合并成為非林地，因而非林地斑塊較為分散，分布比較均勻，與周邊的連通性較好。在9類景觀類型中，聚合度指數(shù)差別不大，且聚合度均大于60，表明五指山市景觀較為聚集，各景觀類型之間排列較為緊密。闊葉混交林斑塊密度PD最大為1.055 5，表明闊葉林出現(xiàn)頻數(shù)較大，為五指山市主要的景觀基質。針葉林的為邊緣密度ED最大為31.082 1，其他灌木類的邊緣密度ED最小為0.4404。表明，針葉林斑塊與其他類型斑塊接觸面較大，不易受到危害。而其他灌木林則相反，易遭到破壞，穩(wěn)定性較差。景觀分離度指標中，它的單位是比率，取值范圍為[0,1)，當景觀只有一個斑塊時，DIVISION=0；當景觀最大限度的細化時，DIVISION取得最大值。說明：DIVISION的基礎是累計斑塊面積的分布，可以理解為在景觀中隨機選擇的兩個像元，這兩個像元不屬于同一個斑塊類型的概率。同理：求出2003年、2007年五指山市景觀類型指標值。

由圖2可知，五指山各景觀類型的邊緣密度ED多數(shù)呈明顯的增加趨勢，硬闊類與軟闊類景觀的增加較為明顯。闊葉混交林、竹林類、喬經濟林類的邊緣密度ED則呈減小的趨勢。各景觀類型的邊緣密度ED大小排序為：針葉類＞軟闊類＞闊葉混交林＞硬闊類＞其他灌木類＞灌經濟林類＞竹林類＞喬經濟林類＞非林地。

表1 五指山市1998、2003、2007年各景觀類型指數(shù)Table 1 All kinds of landscape indexes of Wuzhishan City in 1998, 2003 and 2007

圖2 1998,2003,2007年邊緣密度動態(tài)變化Fig.2 Dynamic changes of edge density in 1998, 2003, 2007

由圖3可知，五指山市各景觀類型的斑塊密度PD，除竹林類、喬經濟林類均呈明顯的增加趨勢，其中軟闊類變化幅度最為明顯。同年中，各林地景觀類型斑塊密度PD大小排序為：闊葉混交林＞軟闊類＞針葉類＞硬闊類＞其他灌木類＞竹林類＞灌經濟林類＞喬經濟林類＞非林地。

由圖4可知，五指山市各景觀類型的聚合度AI指數(shù)大多呈減小趨勢，只有硬闊類、喬經濟林類、灌經濟林類、其他灌木類呈減小趨勢。表明這四種景觀類型斑塊的發(fā)展和其他景觀類型有所不同，抵抗外界干擾的能力與其他景觀類型有所區(qū)別。

由圖5可知，除非林地之外的各景觀類型的連接度均較小，說明所計算的景觀類型存在3種可能：1）僅由1個或少數(shù)幾個斑塊組成；2）該景觀中所有斑塊類型只含有1個或少數(shù)幾個斑塊數(shù)量；3）斑塊連通性較小。同時各景觀類型的連接度CONNECT動態(tài)變化幅度較小，不存在明顯變化。

圖3 1998,2003,2007年斑塊密度動態(tài)化Fig.3 Dynamic changes of patch density in 1998, 2003 and 2007

圖4 1998,2003,2007年聚合度動態(tài)變化Fig.4 Dynamic changes of polymerization degree in 1998, 2003 and 2007

圖5 1998,2003,2007年連接度動態(tài)變化Fig.5 Dynamic changes of connectivity degree in of 1998, 2003 and 2007

由圖6可知，五指山市各景觀類型的形狀指數(shù)LSI多呈增加的趨勢，竹林類、喬經濟林類呈下降的趨勢。在1998～2007年年間，軟闊類的形狀指數(shù)LSI上升幅度最大，竹林類下降幅度最大。

圖6 1998,2003,2007年形狀指數(shù)動態(tài)變化Fig.6 Dynamic changes of landscape shape index in 1998, 2003 and 2007

3 結論與討論

3.1 結 論

（1）通過分析1998、2003、2007年連續(xù)3期的五指山市林地景觀，由圖1，得1998年闊葉混交林的形狀指數(shù)LSI最高，非林地類最低；針葉林類的破碎度FN最大，非林地類最?。挥蓤D4，非林地類的連接度CCONNECT最大，闊葉混交林最小；針葉林的聚合度AI最大，喬木經濟林的最小，闊葉混交林的斑塊密度PD最大，非林地的最??；針葉林的邊緣密度ED最大，其他灌木林類最小。2003年闊葉混交林的形狀指數(shù)LSI最高，非林地類最低；針葉林類的破碎度FN最大，非林地類最?。环橇值仡惖倪B接度CCONNECT最大，竹林類最小；針葉林的聚合度AI最大，喬木經濟林的最小，闊葉混交林的斑塊密度PD最大，非林地的最?。会樔~林的邊緣密度ED最大，非林地類最小。2007年，闊葉混交林的形狀指數(shù)LSI最高，非林地類最低；針葉林類的破碎度FN最大，非林地類最?。环橇值仡惖倪B接度CCONNECT最大，闊葉混交林最小；針葉林的聚合度AI最大,竹林類的最小，闊葉混交林的斑塊密度PD最大，非林地的最小；針葉林的邊緣密度ED最大，非林地類最小。

（2）近10年內五指山市林地景觀類型分布較為均勻，復雜度基本沒變。各景觀類型的斑塊密度，均呈明顯的增加趨勢，形狀指數(shù)與邊緣密度變化方向一致，各斑塊間相互影響及轉化較小，景觀斑塊間類型的改變主要是通過林地與非林地間的轉化完成，景觀聚合度的年變化不大。

3.2 討 論

（1）運用目視解譯進行判讀，人為主觀影響較大，將五指山市的景觀類型分為闊葉混交類、喬木經濟林類、硬闊類、針葉類、軟闊類、竹林類、灌木經濟林類、其他灌木類、非林地類等九大景觀類而非小班類型的研究尺度，研究尺度相對較大，因此將會給所預測結果造成一定的誤差[12]。

（2）景觀格局的變化是受人為因素、自然條件因素、社會發(fā)展因素等眾多驅動力所影響的。本文主要是對景觀格局的變化進行分析及預測，鮮有討論對引起景觀因素變化的驅動力的分析，在今后的研究過程中，可嘗試從此方面著手進行分析[13-14]。

[1]李哈濱,王政權,王慶成. 空間異質性研究理論與方法[J].應用生態(tài)學報, 1998,9(6):651-657.

[2]John AB. Wildlife and Landscape Ecology [M]. New York;Springer-verlag, 1997:881-908.

[3]Michael RM. Landscape Ecology and Management.Management[M]. Montreal: Polyscience Publications Inc., 1987:1320-1336.

[4]Gordon M, Baudry J, Formand RTT. Trans. By Zhao-Y.Thermodynamics foundation and information theory of landscape heterogeneity[J]. Chinese Journal of Ecology, 1995, 14(2):27-36.

[5]肖篤寧,布仁倉,李秀珍. 生態(tài)空間理論和景觀異質性[J].生態(tài)學報, 1997,17(5):453-461.

[6]郭晉平,陽漢熙,薛俊杰.關帝山森林景觀異質性及其動態(tài)的研究[J].應用生態(tài)學報,1999,10(2):167-168.

[7]趙 羿,李月輝. 實用景觀生態(tài)學[M].北京:科學出版社,2001:256-257.

[8]郭晉平. 森林景觀生態(tài)研究[M].北京:北京大學出版社,2001: 18-20.

[9]Galbraith L M, Burns C W. Linking land-use, water body type and water quality in southern New Zealand[J]. Landscape Ecology, 2007, 22(2):231-241.

[10]馬文俊. 西洞庭湖區(qū)森林景觀斑塊劃分及空間格局研究[D].長沙：中南林業(yè)科技大學,2011.

[11]鄭新奇,付梅臣. 景觀格局空間分析技術及其應用[M].北京:科學出版社,2010.

[12]劉立武,佘濟云,程玉娜.五指山市喬木林動態(tài)變化特征分析[J].中南林業(yè)科技大學學報,2014,34(4):1-5.

[13]張煜星,嚴恩萍,夏朝宗,等.基于多期遙感的三峽庫區(qū)森林景觀破碎化演變研究[J].中南林業(yè)科技大學學報,2013,33(7):1-7.

[14]李瑞霞.杭州市景觀格局及動態(tài)變化研究[D].金華:浙江師范大學,2013.

Study on heterogeneity of forest landscape in Wuzhishan city, Hainan province

SHE Yu-chen, CHEN Cai-hong, LU Yu, CHEN Dong-yang
(Central South University of Forestry and Technology, Changsha 410004, Hunan, China)

Based on the data of Wuzhishan city’s 1998, 2003, and 2007 three period remote sensing images, 2010 forest resources II class investigation data and 1993, 1998, 2003, and 2008 four period category inventory data, by adopting applying landscape ecology,sustainable forest management and geostatistical statistics theory and method, six landscape pattern indexes, including degree of polymerization (AI), landscape shape index(LSI), Edge density(ED), Fragmentation index(FN), Connectivity index(CONNECT) and patch density index (PD) were selected, the landscape pattern indexes of different woodland landscape types of Wuzhishan city were computed by using software of landscape fragstats 4.0. On the basis of above results, Wuzhishan forest landscape heterogeneity was comprehensive analyzed and studied. The results show that in recently 20 years, the distribution of Wuzhishan city forest landscape types have been relatively uniform, and less complexity basic changed; Various types of forest landscape patch density all showed a trend of obvious increase and shape index and edge density changed in the same direction, the mutual influence and transformation between the plaques and were smaller, the landscape type changes were accomplished mainly through the transformation of forest-land and nonforestland, the annual change of landscape polymerization degree didn’t changed much. The research findings have great theory and practical guiding significance in dynamic evaluation of forest landscape, ecological vulnerability analysis andfForest land comprehensive management.

woodland landscape heterogeneity; remote sensing image; forest landscape pattern index; Wuzhishan city of Hainan province

S718.55+1.1

A

1673-923X(2015)11-0136-06

10.14067/j.cnki.1673-923x.2015.11.024

2015-01-10

湖南省軟科學重點項目（2013ZK2031）；國家林業(yè)公益性行業(yè)項目（201004032）；海南省林業(yè)廳重點科研項目（201108478）

佘宇晨，碩士研究生

陳彩虹，教授； E-mail：chencaihong056@163.com

佘宇晨，陳彩虹，陸 禹，海南省五指山市林地景觀異質性研究[J].中南林業(yè)科技大學學報, 2015, 35(11): 136-141.

[本文編校：吳 毅]