基于粒度反推法的昌江縣景觀格局優化

陳 怡，陸 禹

(國家林業和草原局中南調查規劃設計院，長沙 410014)

生態環境是人與自然共同作用的結果[1]，不同的景觀斑塊組成了異質的景觀格局[2]。景觀生態學認為，不同的景觀類型對生態過程有不同的影響機制[3-4]，合理的景觀格局是實現可持續發展的基礎[5]。景觀格局通過關鍵局部控制生態過程，從而促進區域協調可持續發展[6]。利用最小累積阻力(minimal cumulative resistance，簡稱MCR)模型[7-8]通過分析生態關鍵局部的空間結構[9-10]，進而構建生態網絡，是實現景觀格局優化的主要方法之一。目前對生態網絡的研究集中在整體資源整合優化上[11]，對于組成元素的特征結構研究較少。優化的準確性基于生態源地選擇的合理性，粒度反推法[8，12]從景觀格局現狀出發，反推最佳生態源地空間結構的思路具有啟發意義，為科學選擇生態源地提供了基礎。同時生態節點對于提升生態結構和強化生態廊道有關鍵作用，是生態網絡構建的重要組成內容，在優化中應加以重視。

昌江縣生態旅游資源豐富，是海南西部旅游發展的中心，目前正以“山海黎鄉大花園”為目標建設森林城市，生態環境質量對其發展至關重要[13]。在開展建設的同時不可避免對生態環境產生破壞，加之旅游人數的與日俱增，當地出現了水土流失、植被退化、局部環境質量低的現象[14]，如何取得經濟社會和生態雙贏是昌江縣急需解決的關鍵問題。本研究采用空間網絡分析、粒度反推法和主成分分析等方法對昌江縣進行景觀格局優化研究，旨在探索主要的生態結構，明確重要的生態節點，提出科學的優化措施，以促進昌江縣生態和經濟協調發展。

1 研究區概況

昌江縣是海南省西部沿?？h城， 108°38′—109°17E，18°53′—19°30′N，總面積161 952 hm2，北、東、南面分別與儋州市、白沙縣、樂東縣、東方市接壤，土地利用類型以耕地和林地為主，分別占總面積的26.19%和65.27%。西北部多為平原，東南部以丘陵為主，海拔由沿海向內陸抬升；熱帶季風氣候，全年無冬，年均氣溫24℃，年均降雨量1 600 mm，昌化江自東南向西北沿南緣入海，水源充足，土地肥沃，具有發展熱帶高效農業的優良條件。

2 數據與方法

2.1 數據來源及處理

以昌江縣2019年多規合一數據、林地變更數據、DEM數字高程模型，以及2018年遙感影像為研究基礎，在大地2000坐標系中配準各類數據，將景觀類型分為林地、園地、耕地、草地、水域、未利用地和建設用地7種類型[11-12，15]，通過實地調查校對調整景觀類型數據。

2.2 研究方法

2.2.1 景觀格局優化方法

基于MCR模型的景觀格局優化通過構建生態網絡實現，主要包括四個步驟：

1)生態源地的提取。生態源地由生態服務功能較高的景觀類型組成[8，15]，本研究選擇林地、水域和草地。利用粒度反推法和主成分分析法確定生態源地的空間位置，首先從分析的連續性出發，以50 m為間隔尺度，形成16個50～800 m不同研究尺度的昌江縣生態景觀組分結構；然后從整體性和連通性兩方面選取組分數、最大組分斑塊數、連接度、分離度、聚合度、鄰接比率、平均臨近距離、內聚力、連接度增加百分率等9個測定指標來分析景觀組分結構，其中組分數用于判斷生態源地的數量，其他指標判斷景觀組分結構的穩定性，利用主成分分析計算各景觀組分結構的整體得分，確定最優景觀組分結構；進而根據最優景觀組分的空間位置選擇相連或相隔很近的林地、水域、草地組成生態源地。

2)生態阻力面的構建。參考相關研究[10]，以景觀類型的單位面積生態服務價值為標準，利用空間分析和地統計分析，從顯性和隱性生態阻力兩方面構建生態阻力面[10]，景觀類型阻力值[8]參照見表1。同時以空間主成分分析法，基于顯性和隱性生態阻力面構建綜合生態阻力面，根據DEM計算坡度，對生態阻力值進行修正。

表1 景觀類型的生態阻力值景觀類型單位面積生態服務價值(元/hm2·a)生態阻力值林地19334.020園地12870.345耕地6114.355草地6406.570水域40676.41未利用地371.480建設用地-8852.1100

3)生態廊道的計算。利用MCR模型分析各生態源地之間的最小累積阻力路徑，即為生態廊道的空間位置。

4)生態節點的確定。以修正后的生態阻力為高程模型，利用ArcGIS空間分析模塊計算生態阻力的最大路徑，與生態廊道的交叉點為生態節點。網絡分析能基于空間拓撲關系、跨度以及網絡內在聯系[16]等進行空間狀態分析，在充分考慮生態系統連接度和生態廊道空間位置的條件下，利用網絡分析產生每個生態節點的服務范圍，以得到生態節點的面積和形狀。

2.2.2 粒度反推法

粒度反推法是基于反推思想，從規模方面利用不同尺度的柵格模擬生態系統中的景觀類型，通過景觀測定指標確定目標特征的最優空間結構，進而選取興趣區域的一種方法[10]。其思想基于生態學，認為結構連接性良好、規模較大的生態系統更穩定。目前該方法主要應用在生態源地分析提取、自然保護地整合優化、動植物保護等方面。在生態源地提取中，粒度反推法利用不同尺度的柵格，反映生態系統主要組成結構的變化，從生態系統規模和穩定性的角度分析確定生態源地潛在位置，配合主成分分析等方法，實現生態源地的客觀提取。

3 結果與分析

3.1 生態源地

提取昌江縣的林地、水域和草地景觀，生成50～800 m 16個不同粒度的景觀組分結構，并在Fragstats中計算景觀組分結構的景觀格局指數，結果見表2。

先對測定指標進行標準化處理。分離度和平均臨近距離與其他景觀格局指數的指向結果相反，取相反數進行分析。以累積貢獻率大于80%為標準，共提取2個主成分(表3)。主成分1中連接度、最大組分斑塊數、聚合度、內聚力、鄰接比率和平均臨近距離荷載較高，反映了景觀組分結構的集聚性，為集聚性指標；主成分2中連通性增加百分率和分離度指標荷載較高，反映了景觀組分結構隨粒度變化表現出的整體分離狀態，為分離性指標。

利用荷載系數除以特征值的平方根，計算主成分的特征向量，其與標準化數據矩陣相乘得到主成分表達式：

Z1=0.350X1-0.384X2+0.245X3+0.446X4+

0.348X5-0.153X6+0.448X7-0.354X8

(1)

Z2=0.051X1+0.223X2+0.512X3-0.061X4+

0.431X5+0.687X6-0.054X7+0.148X8

(2)

將數據代入(1)式和(2)式，得到測定指標在主成分上的得分，與相應主成分的方差貢獻率相乘計算各粒度的景觀組分結構的整體連通性得分(圖1)。粒度的增加逐漸簡化生態系統的主要組成部分，面積較小的零星斑塊逐漸融入較大組分或被剔除，增加了組分之間的距離，因此隨粒度的增加景觀組分的整體連通性得分呈現總體下降趨勢。得分在下降的過程中存在趨勢突然改變的情況，表明景觀組分達到適宜的結構，即達到了生態系統的極簡形態，又保留了必要的關鍵組成部分。昌江縣的景觀組分結構在600 m粒度時得分明顯高于兩側粒度得分，說明該粒度的景觀組分結構最優。

參照600 m粒度的景觀組分結構，從景觀格局現狀提取生態景觀斑塊形成景觀組分，并根據景觀組分的空間連接特征確定生態斑塊的增加位置，轉換地類以形成生態源地。600 m粒度的景觀組分結構存在組分數23，表明生態系統存在23個生態源地，其中9個在現狀中為彼此相連的昌化江景觀組分，合并形成1個生態源地，最終得到14個生態源地(圖2)。根據景觀組分的空間連接特征，共21塊非生態景觀斑塊需轉換成生態斑塊，以保證生態源地的完整性和生態功能發揮，其中17個耕地斑塊，面積12.90 hm2；未利用地斑塊1個，面積0.04 hm2；建設用地斑塊3個，面積1.65 hm2。

圖2 生態源地空間分布

昌江縣西北部和東南部的生態源地規模較大，中部規模較小，這與地形地貌基本吻合。中部平原居多，海拔在200 m以內，平坦的地形有利于農業發展，耕地分布范圍大。耕地的經濟效益較高，是糧食生產的安全保障，人為干擾強烈，該區域的生態斑塊遭到耕地蠶食并在空間上被隔斷，規模小、數量多，難以形成大規模的生態源地。昌江縣東南部為霸王嶺國家森林公園，中北部為保梅嶺自然保護區，地形起伏大，用地開發受到限制，生態系統保持完好，森林景觀相互連接，組成了最大的生態源地。昌江縣西北沿海主要為低矮丘陵，海拔400 m左右，包括昌化林場和大面積的國有防護林、后備土地資源，植被以混交林、灌木林和桉樹為主，多連接成片形成生態源地。

昌江縣生態源地總面積為104 759.57 hm2，占陸地面積的64.69%，生態景觀數量多，可作為海南西部的生態核心發揮生態服務功能。最大生態源地的面積為79 954.34 hm2，以天然混交林為主，生態服務功能價值高，應重點保護；最小生態源地的面積為162.44 hm2，位于老洋地村和唐村交界處，植被為橡膠林，生態源地的邊緣較規整，說明受規劃影響較深，人為干擾強烈。生態源地規模的差異反映出生態演替機制的不同，在生態建設中應差別對待。

3.2 生態阻力面

計算任意空間位置到生態源地的顯性生態阻力和隱性生態阻力，利用空間主成分分析進行疊加，對比預測模型的擬合程度，選擇J-Bessel模型進行模擬預測，形成綜合生態阻力面，擬合優度為0.7642。利用坡度修正系數對綜合生態阻力面進行修正(圖3)。通過對比顯性生態阻力和隱性生態阻力的差異局部，得出生態陷阱的位置，包括叉河鎮的大阪田、烏烈鎮的峨溝村、石碌鎮的山竹溝村和牙營村、七叉鎮的燕窩嶺等，這些區域多為被耕地包圍的林地或礦區，林地斑塊面積小、破碎化程度高，難以維系自身生態系統的良性發展，是生態建設中需要重點建設但又容易被忽略的區域。

圖3 阻力面的構建

3.3 生態廊道

昌江縣有生態廊道32條(圖4)，主要分布在中部平原地區，總長度51 644 m。中部的生態廊道較長，四周的較短，說明中部地區受人為干擾強烈，生態景觀組分相隔較遠。根據景觀基質將生態廊道分為381段，其中：耕地生態廊道25 235 m、林地生態廊道19 476 m、水域生態廊道6 228 m、建設用地生態廊道614 m、園地生態廊道88 m。耕地生態廊道和林地生態廊道是將昌江縣生態系統連接成整體的主體，分別占生態廊道長度的48.86%和37.71%，耕地是阻隔生態過程的主要因素。

耕地是農業生產的主要載體，農藥化肥使用較多，耕地生態廊道的主要作用是防止農業污染影響生態過程，應加寬生態廊道的寬度，形成緩沖帶，并以抗農藥化肥污染強的植物進行建設。林地生態服務功能較高，與生態源地功能相似，應從林地的垂直和水平空間結構出發加強建設，提升生態穩定性。

圖4 生態廊道和生態節點的空間分布

3.4 生態節點

昌江縣有生態節點29個(圖4)，其中11個位于林地、15個位于耕地、3個位于水域。從生態節點的布局來看，可將中部地區分為三條東西向帶狀區域，耕地生態節點位于北部和南部帶狀區域，林地和水域生態節點位于中部帶狀區域。為滿足生產生活的需要，中部帶狀區域建有白馬山水庫等大面積水域，以及石碌河等河流流經區域，沿河分布有保平林場，林地和水域生態節點主要分布于此。

生態節點能有效提高生態系統的整體性和生態廊道的穩定性，采用粒度反推法分析存在生態節點時生態系統的連接度，進而確定生態節點的合適尺度(圖5)，得到350 m粒度時景觀組分連接度最強，可參考該尺度建設生態節點。

利用網絡分析根據生態節點的合適尺度、空間位置、生態廊道的走勢計算生態節點的規模和形狀(圖6)，共形成生態節點范圍29個，總面積201.49 hm2，其中最大面積20.71hm2，最小面積5.21 hm2。

圖5 連接度分析

圖6 生態節點的空間形態

4 結論與建議

4.1 結論

1)利用粒度反推法和主成分分析，從整體連通性的角度計算得出600 m粒度的生態景觀組分是昌江縣的最優景觀組分結構，昌江縣有生態源地14個，中部地區數量多、規模小，東南部和西北沿海數量少、規模大，源地規模的差異反映出生態演替機制的不同。生態源地的分布特征與人為干擾存在明顯的相關性，地形地貌是影響其空間分布的主要因素。昌江縣需將21塊非生態景觀斑塊轉換為生態景觀，總面積14.59 hm2，以連接生態景觀組分形成生態源地。

2)結合顯性和隱性生態阻力，及地形地貌構建了綜合生態阻力面，顯性和隱性生態阻力的差異指明了生態系統中存在的脆弱點，能有效指導生態建設。

3)昌江縣有生態廊道32條，長度51 644 m，以耕地和林地生態廊道為主，不同基質的生態廊道功能和建設措施不同；有生態節點29個，空間位置主要受景觀類型和人為干擾的影響。研究從生態系統連接度和生態廊道走勢的角度，采用網絡分析計算了生態節點的空間形態，既能提高生態系統的連接度，又能保護生態廊道的穩定性。昌江縣生態節點的合適建設尺度為350 m，共形成29個生態節點區域，總面積201.49 hm2。

4.2 建議

1)昌江縣生態源地在中部地區數量多、規模小，呈現破碎化，同時生態廊道和生態節點也集中在此區域，這與中部地區的農業發展建設息息相關。為實現區域協調可持續發展，在發展生產的同時應重視生態環境建設，著力加強中部區域生態源地的保護和管理，做好生態廊道和生態節點的建設，實施區別化空間管制，采取限制人類活動、設置緩沖帶等措施保護生態網絡。

2)根據生態阻力面顯示的生態脆弱點，主要為叉河鎮的大阪田、烏烈鎮的峨溝村、石碌鎮的山竹溝村和牙營村、七叉鎮的燕窩嶺等被耕地包圍的林地或礦區，這些區域是生態系統中的生態陷阱，應加強生態建設防治生態退化。一是要加強該區域生態監測和管理，確保植被不退化、不遭到人為破壞；二是要降低周邊區域人為干擾，通過生態建設擴大區域內生態景觀斑塊的規模；三是嚴格管理礦區生產，及時做好污染防范措施和實施生態恢復。

3)根據生態廊道和生態節點的類型進行生態網絡建設，建設采取的方式與措施應和生態廊道、生態節點相適應。生態廊道中水域生態廊道和耕地生態廊道較多，水域生態廊道可結合美麗鄉村做好生態建設，兼顧景觀價值，促進昌江縣全域旅游發展；耕地生態廊道周邊應盡量減少農藥、化肥的使用，降低污染風險。生態節點對強化生態系統意義重大，同時對區域協調發展價值明顯，應切實做好規劃，與城鄉空間規劃對接，在發揮生態節點生態保護價值的同時，將其規劃成各類城市、鄉村生態公園，發揮生態和社會效益。