賀 丹，陳彩虹，佘濟云，陸 禹，董雅雯

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

基于生態(tài)阻力面綜合構(gòu)建法的三亞市景觀格局分析

賀 丹，陳彩虹，佘濟云，陸 禹，董雅雯

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

為提高景觀格局穩(wěn)定性，增強生態(tài)網(wǎng)絡(luò)的空間連通性，以MCR模型為基礎(chǔ)，采用生態(tài)阻力面綜合構(gòu)建法、熱點分析法、核密度分析法等方法，分析構(gòu)建了海南省三亞市景觀的生態(tài)源地、生態(tài)廊道、生態(tài)節(jié)點等組分，在此基礎(chǔ)上對三亞市的景觀格局進行了綜合分析。結(jié)果表明：三亞市生態(tài)節(jié)點密度分布并沒有出現(xiàn)特別松散或者特別集中的區(qū)域，其中密度較大的區(qū)域主要分布在吉陽區(qū)東北部和天涯區(qū)西南部，與熱點分析的高值聚類區(qū)域和低值聚類區(qū)域基本相符。其中，三亞市一級生態(tài)節(jié)點共10個，全部分布在吉陽區(qū)，屬非常重要生態(tài)節(jié)點，應(yīng)優(yōu)先建設(shè)；二級生態(tài)節(jié)點共14個，絕大部分分布在天涯區(qū)，屬重要生態(tài)節(jié)點，應(yīng)重點建設(shè)；三級生態(tài)節(jié)點共33個，隨機分布在三亞市的四個區(qū)，屬一般生態(tài)節(jié)點，應(yīng)加強建設(shè)。研究成果可為三亞市生態(tài)環(huán)境建設(shè)提供科學(xué)的參考依據(jù)。

MCR模型；生態(tài)阻力面；熱點分析；核密度分析；景觀格局；三亞市

三亞市地處海南島南端，是中國最南部的熱帶濱海旅游城市，居中國四大一線旅游城市之首。旅游業(yè)的蓬勃發(fā)展為當(dāng)?shù)貛砹司薮蟮慕?jīng)濟效益，同時也給三亞市的生態(tài)環(huán)境帶來了巨大的威脅[1]。因此，改善三亞市生態(tài)環(huán)境，實現(xiàn)生態(tài)過程的平衡發(fā)展，提出景觀格局優(yōu)化措施，將為三亞市生態(tài)環(huán)境改善、旅游業(yè)可持續(xù)發(fā)展提供科學(xué)的依據(jù)，對三亞市生態(tài)文明建設(shè)意義重大。目前，國內(nèi)外對于景觀格局的研究主要是以Fomran[2]提出的景觀格局優(yōu)化理論和俞孔堅[3]提出的景觀生態(tài)學(xué)理論為基礎(chǔ)，以Knaapen[4]提出的最小耗費距離模型為方法，對不同區(qū)域進行了景觀格局優(yōu)化的研究，如張海鷗[5]以成都市為研究對象，采用最小耗費距離模型建立阻力面，對成都市不同時期景觀生態(tài)安全格局動態(tài)變化進行了初步研究；劉杰等[6]以滇池流域為研究對象，以GIS空間技術(shù)為基礎(chǔ)，采用最小耗費距離模型構(gòu)建了生態(tài)源地、生態(tài)廊道、生態(tài)節(jié)點等景觀組分，對區(qū)域景觀格局進行了優(yōu)化。雖然景觀格局研究目前取得了較多的成就，但仍有很多問題值得探討，如事物在空間上的相互影響對生態(tài)阻力面的影響[7]，所選取的生態(tài)節(jié)點相互之間的影響及建設(shè)順序等。本文以三亞市為研究對象，采用綜合阻力面構(gòu)建法構(gòu)建生態(tài)阻力面、采用核密度分析和熱點分析對生態(tài)節(jié)點相互之間的影響及建設(shè)順序進行研究，探討三亞市景觀格局優(yōu)化措施，為三亞市可持續(xù)發(fā)展提供科學(xué)依據(jù)。

1 研究區(qū)概況

海南省三亞市位于海南島最南端，地處北緯 18°09′34″～ 18°37′27″、 東 經(jīng) 108°56′30″～109°48′28″，東西長約 91.6 km，南北寬約 51 km，總面積1 919.58 hm2。三亞市東鄰陵水縣，西毗樂東縣，南臨南海，北依保亭縣，地勢北高南低，北部為山地和丘陵，南部為平原和臺地，自北向南逐漸傾斜。三亞屬于熱帶海洋性季風(fēng)氣候區(qū)，受南海海洋氣候影響較大，長夏無冬，全年高溫多雨，氣溫年較差小，有“天然溫室”之稱。三亞市范圍內(nèi)流域面積大于100 hm2的河流有4條，南部水庫眾多，水資源豐富。三亞市林地面積14.43萬hm2[8-10]，占全市土地面積的75.5%。

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源與處理

以三亞市1∶10 000地形圖和2010年二類調(diào)查數(shù)據(jù)為基礎(chǔ)，結(jié)合2015年遙感影像進行解譯，合并地類，得到三亞市土地利用類型分布圖。通過給不同地類賦阻力值和克里金插值法分別得到顯性與隱性阻力面，運用空間主成分分析得到綜合阻力面，再通過計算距離成本得到累計耗費表面。結(jié)合三亞市生態(tài)源地和累計耗費距離表面進行空間分析和水文分析，構(gòu)建三亞市生態(tài)廊道和生態(tài)節(jié)點，并對生態(tài)節(jié)點進行熱點分析和核密度分析，得出三亞市景觀格局現(xiàn)狀。

2.2 研究方法

2.2.1 生態(tài)源地選取

篩選生態(tài)服務(wù)功能價值高的景觀類型定義為生態(tài)斑塊，從中挑選規(guī)模較大、結(jié)構(gòu)較優(yōu)、對生態(tài)過程影響有決定性作用的生態(tài)斑塊集合體作為生態(tài)源地[7]。

2.2.2 阻力面構(gòu)建

(1)顯性阻力面構(gòu)建

顯性生態(tài)阻力面是從土地利用類型考慮，對不同的景觀類型賦上相應(yīng)的阻力值[11-12]（表1）。

表1 各景觀類型對應(yīng)阻力值Table 1 The resistance value of different landscape patterns

(2)隱性阻力面構(gòu)建

隱性生態(tài)阻力值是考慮事物在空間上的相互影響，采用克里金插值法構(gòu)建，公式為：

式中：Z(x0)為未知樣點值；Z(xi)為未知樣點周圍的已知樣點值，λi為第i個已知樣本點的權(quán)重，n為已知樣本點的個數(shù)[13]。

(3)綜合阻力面構(gòu)建

綜合阻力面同時考慮前面兩者，采用生態(tài)阻力面綜合構(gòu)建法，對顯性生態(tài)阻力面和隱性生態(tài)阻力面進行加權(quán)計算，公式為：

式中，Cr為綜合阻力面，Dr為顯性阻力面，Rr為隱性阻力面，a、b分別為顯性阻力面和隱性阻力面所占權(quán)重（a+b=1），顯性阻力面權(quán)值為0.7，隱性阻力面的權(quán)值為0.3[7,12]。

2.2.3 生態(tài)廊道構(gòu)建

生態(tài)廊道是運用Arcgis中空間分析的最短路徑分析得出，公式如下：

式中：MCR為最小累積阻力值，Dij為任意景觀單元i與生態(tài)源地j在空間上的距離，Ri是指在這一運動中景觀單元i的生態(tài)阻力值阻力值，m是指景觀類型數(shù)量，n是指基本的單元總數(shù)[14-15]。

2.2.4 生態(tài)節(jié)點分析

生態(tài)節(jié)點是運用Arcgis中空間分析的山體陰影工具得到[16-17]。

(1)熱點分析

其中xj是要素j的屬性值，wi,j是要素i和j之間的空間權(quán)重，n為要素總數(shù)，且：

(2)核密度分析

核密度分析是通過離散點數(shù)據(jù)或者線數(shù)據(jù)進行內(nèi)插的方法，對落入搜索區(qū)的點賦予不同的權(quán)重，離中心越近的點權(quán)重越大，越遠則越小。

3 結(jié)果與分析

3.1 生態(tài)源地和生態(tài)阻力面分析

3.1.1 生態(tài)源地分析

選取寧遠河、三亞河、藤橋河等67個面積大于10公頃的連續(xù)水域[18]作為“生態(tài)源地”（圖1），總體分布特征是南部數(shù)量少規(guī)模大，北部數(shù)量多規(guī)模小。分析原因為三亞市南部沿海地區(qū)為三亞市城區(qū)，景觀布局經(jīng)過人為設(shè)計與重點建設(shè)，景觀破碎度低，生態(tài)源地規(guī)模大；北部耕地和森林廣泛，很大程度上保持了原生態(tài)景觀類型，因而景觀破碎度高，生態(tài)源地規(guī)模小。

圖1 三亞市景觀格局分布Fig.1 The landscape pattern distribution of Sanya

3.1.2 生態(tài)阻力面分析

顯性生態(tài)阻力面僅從景觀類型阻力值考慮，景觀類型分布復(fù)雜，所以顯性阻力面的阻力值相互之間差異大；隱性生態(tài)阻力面考慮了不同景觀類型之間的相互影響，所以阻力值差異小。綜合生態(tài)阻力面（圖2）是顯性生態(tài)阻力面和隱性生態(tài)阻力面的疊加，綜合考慮了不同景觀類型阻力值的差異及其相互之間的影響，比前面兩者更真實、全面地反映現(xiàn)實生態(tài)系統(tǒng)的情況。三亞市西南部沿海區(qū)域的生態(tài)阻力值最大，因該區(qū)域是三亞市的城市中心地區(qū)，建設(shè)用地較多，不利于生態(tài)流通過；而東南部和中部區(qū)域森林、水域分布眾多，建設(shè)用地稀少，生態(tài)阻力值較小，生態(tài)流運行暢通。

3.2 生態(tài)廊道分析

根據(jù)生態(tài)廊道穿過的景觀類型以及自身的功能特征可以分為建設(shè)用地生態(tài)廊道、耕地生態(tài)廊道和森林生態(tài)廊道三種類型，三亞市生態(tài)廊道主要為森林生態(tài)廊道（表2）。

圖2 三亞市綜合阻力面Fig.2 Comprehensive resistance surface of Sanya

表2 三亞市各類型生態(tài)廊道Table 2 The types of ecological corridor in Sanya

三亞市建設(shè)用地生態(tài)廊道是連通人類社會經(jīng)濟環(huán)境與自然生態(tài)系統(tǒng)環(huán)境的重要渠道，主要分布地是三亞市城區(qū)，主要土地利用類型為建設(shè)用地，生態(tài)流易受到干擾與阻礙，應(yīng)加強保護并維持建設(shè)用地生態(tài)廊道的現(xiàn)有規(guī)模，維護三亞市建設(shè)用地生態(tài)廊道的穩(wěn)定性；耕地生態(tài)廊道是阻隔農(nóng)業(yè)生產(chǎn)帶來的有毒有害物質(zhì)進入生態(tài)源地的天然屏障，分布地主要土地利用類型為耕地，生態(tài)流的流通性較低，應(yīng)保持耕地生態(tài)廊道的現(xiàn)有寬度，設(shè)置農(nóng)業(yè)生產(chǎn)隔離區(qū)域，保證農(nóng)業(yè)生產(chǎn)和生態(tài)環(huán)境的可持續(xù)發(fā)展；森林生態(tài)廊道是連接生態(tài)源地和森林景觀的橋梁，主要分布地是三亞市森林覆蓋較高的丘陵地帶，主要土地利用類型為林地，生態(tài)功能與生態(tài)源地相似，是三種生態(tài)廊道中要優(yōu)先重點建設(shè)的一類生態(tài)廊道。

3.3 生態(tài)節(jié)點分析

生態(tài)節(jié)點是生態(tài)廊道與生態(tài)阻力面最大阻力值的交匯處，也是生態(tài)廊道的最薄弱處，應(yīng)加強建設(shè)。由熱點分析可知（圖3），三亞市中部偏東地區(qū)的紅色區(qū)域為阻力值的熱點區(qū)域，Z＞1.65（P＜0.10），表示這個區(qū)域生態(tài)節(jié)點阻力值高且被高值包圍，空間自相關(guān)性強且為正相關(guān)；中部偏西地區(qū)的藍色區(qū)域為阻力值的冷點區(qū)域，Z＜-1.65（P＜0.10），表示這個區(qū)域的生態(tài)節(jié)點阻力值低，且其周邊的生態(tài)節(jié)點阻力值也低，空間自相關(guān)性強且為正相關(guān)；其余節(jié)點Z得分在-1.65到1.65之間，表示這些生態(tài)節(jié)點空間自相關(guān)性弱，為隨機分布。從核密度角度分析，三亞市生態(tài)節(jié)點密度分布并沒有出現(xiàn)特別松散或者特別集中的區(qū)域，其中密度較大的區(qū)域主要分布在吉陽區(qū)東北部和天涯區(qū)西南部，與熱點分析的高值聚類區(qū)域和低值聚類區(qū)域基本相符。

生態(tài)節(jié)點分為森林生態(tài)節(jié)點、耕地生態(tài)節(jié)點和建設(shè)用地生態(tài)節(jié)點三種類型。為了對生態(tài)節(jié)點的建設(shè)順序進行綜合評價，根據(jù)熱點分析的Z得分和P值，進一步將各類型生態(tài)節(jié)點劃分為一級生態(tài)節(jié)點（Z＞1.65，P＜0.10）、二級生態(tài)節(jié)點（Z＜-1.65，P＜0.10）和三級生態(tài)節(jié)點（-1.65≤Z≤1.65，P≥0.10）。

三亞市共57個生態(tài)節(jié)點，森林生態(tài)節(jié)點46個、耕地生態(tài)節(jié)點8個、建設(shè)用地生態(tài)節(jié)點3個。森林生態(tài)節(jié)點數(shù)量最多，主要是由于三亞市森林資源豐富，林地面積大，森林生態(tài)節(jié)點在整個生態(tài)系統(tǒng)中處于戰(zhàn)略地位，對于生態(tài)流在生態(tài)源地和森林之間的暢行有重大作用；耕地生態(tài)節(jié)點與建設(shè)用地生態(tài)節(jié)點分別分布在三亞市東部農(nóng)田區(qū)域與東南沿海城區(qū)，二者皆受人為干擾的影響，耕地生態(tài)節(jié)點的作用是保護耕地生態(tài)廊道，以防農(nóng)業(yè)活動破壞生態(tài)廊道；而建設(shè)用地生態(tài)節(jié)點對提高城區(qū)人民生活質(zhì)量和優(yōu)化城區(qū)生態(tài)安全有直接作用，是城區(qū)生態(tài)建設(shè)的重點。三亞市一級生態(tài)節(jié)點共10個，全部分布在吉陽區(qū)，屬非常重要生態(tài)節(jié)點，包括5個一級森林生態(tài)節(jié)點、4個一級耕地生態(tài)節(jié)點和1個一級建設(shè)用地生態(tài)節(jié)點，應(yīng)優(yōu)先重點建設(shè)；二級生態(tài)節(jié)點共14個，絕大部分分布在天涯區(qū)，屬重要生態(tài)節(jié)點，包括12個二級森林生態(tài)節(jié)點、1個三級耕地生態(tài)節(jié)點、1個二級建設(shè)用地生態(tài)節(jié)點，應(yīng)重點建設(shè)；三級生態(tài)節(jié)點共33個，隨機分布在三亞市的四個區(qū)，屬一般生態(tài)節(jié)點，包括29個三級森林生態(tài)節(jié)點、3個三級耕地生態(tài)節(jié)點和1個三級建設(shè)用地生態(tài)節(jié)點，應(yīng)加強建設(shè)。

圖3 三亞市生態(tài)節(jié)點核密度分布圖Fig.3 Kernel density of ecological nodes in Sanya

圖4 三亞市生態(tài)節(jié)點等級分布圖Fig.4 Grade distribution of ecological nodes in Sanya

4 結(jié)論與討論

（1）生態(tài)源地南部數(shù)量少規(guī)模大，北部數(shù)量多規(guī)模小，應(yīng)當(dāng)增加南部生態(tài)源地數(shù)量，擴大北部生態(tài)源地規(guī)模，提高其生態(tài)服務(wù)價值；中部地區(qū)的生態(tài)廊道較多，呈環(huán)狀分布，是生態(tài)流高速運行的絕佳通道，南部沿海城區(qū)和北部地區(qū)生態(tài)廊道較少，生態(tài)流流通阻力大，應(yīng)增加該區(qū)域生態(tài)廊道寬度和長度。

表3 三亞市生態(tài)節(jié)點重要性等級劃分Table 3 The classification of ecological nodes’ importance in Sanya

（2）三亞市生態(tài)節(jié)點總體分布比較均勻，局部出現(xiàn)較為集中，呈現(xiàn)“大分散小集中”的格局。生態(tài)節(jié)點按類型、等級進行保護，對于一級森林生態(tài)節(jié)點考慮建設(shè)森林公園，并因地制宜加強對二級、三級森林生態(tài)節(jié)點的保護力度；對于各級耕地生態(tài)節(jié)點可采用退耕還林的方式進行生態(tài)文明建設(shè)；對于各級建設(shè)用地生態(tài)節(jié)點，可考慮增加綠化面積或建設(shè)成小型公園。

本研究將熱點分析和核密度分析等方法引入景觀格局的研究，為三亞市生態(tài)文明建設(shè)提供理論依據(jù)和數(shù)據(jù)支撐。然而，本次研究僅對生態(tài)節(jié)點進行等級劃分，未考慮其形狀特征。因此，如何確定生態(tài)節(jié)點的具體形狀及面積等問題還有待進一步探究。

[1]李 娟, 彭鎮(zhèn)華, 高 健. 三亞市土地利用變化及影響因素驅(qū)動力分析[J]. 安徽農(nóng)業(yè)科學(xué), 2012, 40(36): 17794-17797.

[2]O Neill R V, Kummel J R, Gardner R H, et al. I ndices of landscape Pattern[J]. Landscape Ecology, 1988, 1(3): 153-162.

[3]俞孔堅. 生物保護的景觀生態(tài)安全格局[J]. 生態(tài)學(xué)報, 1999,19(1): 8-15.

[4]Knaapen J P, Scheffer M, Harms B. Estimating habitat isolation in landscape planning[J]. Landscape and Urban Planning, 1992,23:1-16.

[5]張海鷗. 成都市不同時期景觀生態(tài)安全格局動態(tài)變化初步研究[D]. 雅安: 四川農(nóng)業(yè)大學(xué), 2011.

[6]劉 杰, 葉 晶, 楊 婉, 等. 基于GIS的滇池流域景觀格局優(yōu)化[J]. 自然資源學(xué)報, 2012, 27(5): 801-808.

[7]陸 禹, 佘濟云, 陳彩虹, 等. 基于粒度反推法的景觀生態(tài)安全格局優(yōu)化——以海口市秀英區(qū)為例[J]. 生態(tài)學(xué)報, 2015,35(19): 1-10.

[8]季小梅, 張永戰(zhàn), 朱大奎. 三亞海岸演變與人工海灘設(shè)計研究[J]. 第四紀(jì)研究, 2007, 27(5): 853-860.

[9]楊 帆, 楊傳金, 孫 寧,等. 三亞紅樹林景觀特點及保護利用對策[J]. 中南林業(yè)調(diào)查規(guī)劃, 2012, 31(2): 31-34.

[10]周亞東, 周兆德. 基于GIS與Fragstats的海南島森林景觀格局研究[J]. 中南林業(yè)科技大學(xué)學(xué)報, 2015, 35(5): 78-83

[11]Costanza R, Arge R, Groot R, et al. The value of the world’s ecosystem services and natural capital[J]. Nature, 1987, 387:253-260.

[12]YE Yuyao, SU Yongxian, et al. Construction of an ecological resistance surface model and its application in urban expansion simulations[J]. Journal of Geographical Sciences, 2015, 25(2):211-224.

[13]孔紅梅, 韋琳林, 劉國華. Kriging插值法在植物物種地理分布空間格局研究中的應(yīng)用[J]. 生態(tài)環(huán)境學(xué)報, 2010, 19(5):1165-1169.

[14]劉孝富, 舒儉民, 張林波. 最小累積阻力模型在城市土地生態(tài)適宜性評價中的應(yīng)用——以廈門市為例[J]. 生態(tài)學(xué)報, 2010,30(2): 421-428.

[15]周 銳, 王新軍, 蘇海龍, 等. 平頂山新區(qū)生態(tài)用地的識別與安全格局構(gòu)建[J]. 生態(tài)學(xué)報, 2015, 35(6): 2003-2012.

[16]烏尼圖, 岳德鵬, 張亦超, 等. 基于3S技術(shù)的城鄉(xiāng)區(qū)域生態(tài)節(jié)點的識別與分析[J]. 環(huán)境科學(xué)與技術(shù), 2014,37(2):184-189.[17]馬 力, 顧 敏, 吳建明, 等. 京滬高鐵蘇州相城段生態(tài)廊道景觀節(jié)點規(guī)劃研究[J]. 中國園藝文摘, 2014(1): 121-122

[18]任慧君. 區(qū)域生態(tài)安全格局評價與構(gòu)建研究[D]. 北京: 北京林業(yè)大學(xué), 2011.

Landscape pattern analysis based on comprehensive construction method of ecological resistance surface in Sanya

HE Dan, CHEN Caihong, SHE Jiyun, LU Yu, DONG Yawen
(Central South University of Forestry and Technology, Changsha 410004, Hunan, China)

In this paper, landscape pattern optimization of Sanya was studied. It utilizes comprehensive construction method of ecological resistance surface, hot spot analysis and kernel density constructed the ecological source, ecological corridors, ecological nodes,explored the way of optimization of landscape ecological pattern in Sanya, Hainan province.The results showed that: density distribution of ecological nodes in Sanya substantially matches the hot spot analysis.It is not particularly loose or concentrated, and density regions mainly in the northeastern of Jiyang district and southwestern of Tianya district. There are 10 level 1 ecological nodes distributing in Jiyang district,and they are so important that should be given priority to the construction.There are 14 level 2 ecological nodes distributing in Tianya district,and they are so important that should be given more attention to the construction.There are 33 level 3 ecological nodes distributing in Sanya, and they are common ecological nodes with strengthening construction.The study of landscape pattern optimization will provide scienti fi c theoretical references perfect for ecological environment construction and ecosystems in Sanya.

MCR model; ecological resistance surface; hot spot analysis; kernel density; landscape pattern; Sanya

S727.5

A

1673-923X(2017)06-0083-06

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

2016-04-20

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

賀 丹，碩士研究生

陳彩虹，博士，教授；E-mail：shejiyun@126.com

賀 丹，陳彩虹，佘濟云，等. 基于生態(tài)阻力面綜合構(gòu)建法的三亞市景觀格局分析[J].中南林業(yè)科技大學(xué)學(xué)報，2017,37(6): 83-88.

[本文編校：吳 彬]