基于多場景模型的沙漠-綠洲交錯帶林草生態網絡模擬

蘇 凱 于 強 YANG Di 張啟斌 楊 斕 孫小婷

(1.北京林業大學精準林業北京市重點實驗室， 北京 100083； 2.佛羅里達大學地理系， 蓋恩斯維爾 FL 32611)

0 引言

生態系統是人類生存和發展的物質基礎，也是人類的生命支持系統[1]。人類的生存和發展取決于生態系統的穩定性和健康。在干旱地區，沙漠-綠洲生態交錯帶作為缺水沙漠和富水農田系統之間的過渡區，在確保區域生態安全和維護區域內部穩定方面發揮著關鍵作用。然而，沙漠-綠洲生態交錯帶也是脆弱的，其狀態極易受發展策略的影響[2]。當發展策略以經濟優先時，景觀實現了經濟的發展，但導致了地下水位下降，森林遭受破壞，荒漠化加速，濕地退化等生態問題日益嚴重[2-3]，使區域生態系統不平衡且難以恢復[4]。如何有效管理荒漠化和遏制沙漠擴張一直是一個世界性的問題。

生態網絡通過帶狀植被、河流等形成的線性廊道或者是潛在線性生態廊道將荒漠綠洲交錯區相互孤立的生態斑塊連接，形成“點、線、面”相結合的空間生態安全格局，提升景觀的自我調節能力，維持區域生態環境穩定，抵御沙漠化風險[5]。從生態網絡的視角來看，土地沙漠化及其造成的一系列影響的過程實際是生態網絡中生態節點與生態廊道不斷被破壞與消失的過程[6]。在這一過程中，網絡結構的完整性不斷降低，生態功能逐步下降，網絡中生態、物質與能量的交換逐漸受阻。

生態網絡的構建已經形成基本的研究范式：生態源地的識別與潛在生態廊道的提取。識別生態源地，主要是通過評估生態適宜性、生態重要性或生態連通性。其中基于生態系統服務的生態重要性評估最為常見[7]。提取生態廊道的方法通常基于土地覆蓋的價值分配構建阻力面，阻力面描述了物種穿越不同生境斑塊的難度，反映了物種對生態過程的水平阻力。阻力面表征了景觀異質性對生態過程流量的影響。最小成本分析通常用于提取生態走廊。但這種方法忽略了生態源地之間的引力，沒有考慮生態源地之間的相互作用。引力模型早期被應用于城市間的空間結構和零售市場研究[8]，后來也被廣泛用于區域經濟聯系[9]、城市間相互作用結構，以及城市間貿易、物流聯系等研究[10]。目前引力模型廣泛應用于度量區域經濟聯系量或空間相互作用量。

然而，沙漠-綠洲生態交錯帶的森林-草地生態網絡(FG生態網絡)是一個動態網絡，而且是一個脆弱的系統，受自然條件和經濟活動的影響。因此，本文以中國西北地區典型的沙漠-綠洲交錯帶為研究對象，提出林草生態網絡(FG生態網絡)構建的綜合模型，該模型包括兩個子模型，利用改進的生態阻力模型獲取生態流流動的最小累計阻力面，利用改進的生態引力模型計算生態源地之間的引力。通過逐一比較各生態源地之間的合力，利用最小成本分析提取生態廊道，構建連接生態源地的生態網絡。同時，制定多種發展戰略，探討發展戰略中生態保護與經濟發展的比例對生態網絡演變的影響。采用復雜網絡分析方法，對研究區林草生態網絡進行拓撲分析和空間結構分析。旨在為干旱區經濟發展和生態環境保護研究提供新的視角，為干旱區可持續發展提供理論支持。

1 材料與方法

1.1 研究區概況

研究區域(106°7′～107°23′E，39°56′～40°56′N)位于烏蘭布和沙漠東北緣，面積5 923.9 km2，海拔968～1 338 m。研究區東部是中國的母親河——黃河；東北部的平原是中國重要的糧食產區——河套平原；西側是狼山；西南部是烏蘭布和沙漠的腹地。根據1983—2016年磴口沙漠生態站的氣象資料，該地區年平均降水量為144.5 mm，年平均蒸發量為2 398 mm。研究區域的位置如圖1所示。

圖1 研究區位置Fig.1 Location of research area

1.2 數據來源與處理

使用2018年夏季的Sentinel 2遙感影像，分辨率為20 m。數據集由歐洲航天局提供(http:∥www.esa.int/ESA)。通過輻射定標、大氣校正、圖像增強和幾何校正等進行圖像的預處理。利用最大似然監督分類方法對遙感圖像進行分類，并提取研究區域內的土地利用類型信息。從土地利用數據中提取水網和道路網，通過多緩沖分析獲取距離因子數據。從土地利用數據中提取住宅區、道路網和水網，通過密度分析獲取密度因子數據。從處理后的圖像中提取歸一化差異植被指數(NDVI)、修正歸一化差異水體指數(MNDWI)和溫度植被干旱指數(TVDI)。DEM由中國科學院計算機網絡信息中心地理空間數據云平臺(http:∥www.gscloud.cn)提供。地下水埋深的空間分布數據來源文獻[11]。

1.3 方法

1.3.1生態網絡多場景仿真模型

FG生態網絡多場景仿真模型由兩個子模型組成。其中利用生態阻力模型模擬不同發展策略中生態源地之間的最小累計耗費生態阻力，利用生態引力模型模擬生態源地之間的相互作用。通過比較不同發展方案時生態源的綜合受力情況，分析FG生態網絡的變化。

在干旱地區，盡管水資源稀缺，但仍然可以自然形成相對穩定的生態系統[12]。一些研究認為，生態資源地之間存在相互依賴的關系，通過生態網絡，在干旱地區建立了較為穩定的生態格局，維護了區域生態安全。FG生態網絡由生態源和連接生態源的生態廊道組成[13]。在生態網絡中，生態源是提供各種物質和能量的空間單元或生態系統，生態廊道是能量和物質流動的重要通道。在FG生態網絡中，生態源地之間的相互作用會產生生態引力，但同時也會受到來自基質的阻力，生態源地之間互作用的示意圖如圖2所示，圖中F12=F21，為生態引力，fmin表示生態源地之間的最小累積生態阻力。

圖2 生態源地間作用力示意圖Fig.2 Schematic of interaction forces between ecological sources

生態網絡多場景仿真模型的詳細算法步驟如下：

(1)計算生態源地之間的最大生態引力和最小累積阻力。當生態引力大于生態阻力時，源地之間將形成生態廊道，生態流將通過廊道在源地之間傳遞。當生態引力小于生態阻力時，源地之間就沒有生態廊道。當生態引力等于生態阻力時，源地之間的生態廊道將非常脆弱，基質的微小變化可能導致生態廊道斷裂。

(2)重復計算，直至遍歷完成整個FG生態網絡的源地，最終獲得FG生態網絡空間分布網格模擬結果。然后模擬下一個開發計劃中FG生態網絡的空間分布，直到所有發展策略模擬結束。

由于政策的實施是統一的，因此往往不適合當地條件。本文通過調整發展策略中經濟發展和生態保護的比例，設定了11個方案。其中2個是極端發展計劃：沒有經濟發展的純生態保護模式和沒有生態保護的純經濟發展模式。實際上，研究區的發展既是生態保護，也是經濟發展。但是，該地區的生態環境極其脆弱，是控制荒漠化的重要區域。我國政府頒布了更為嚴格的生態保護政策，嚴格限制了該地區的經濟發展，重點是改善和保護生態。因此，研究區目前的發展模式可以看作是一種純粹的生態保護模式。不同情景的發展模式將導致生態環境條件的差異。每個發展計劃的發展策略定義為1，生態保護的比例用m表示，而n代表經濟發展的比例，則m+n=1。當m=1.0時，n=0，(1.0，0)代表純生態保護模式，即以生態保護為重點的研究區發展策略。當m=0，n=1.0時，(0,1.0)代表純經濟發展模型，即研究區域優先考慮經濟發展。最后，獲得了11個開發方案：(0,1.0)、(0.1,0.9)、(0.2,0.8)、(0.3,0.7)、(0.4,0.6)、(0.5,0.5)、(0.6,0.4)、(0.7,0.3)、(0.8,0.2)、(0.9,0.1)、(1.0,0)。

1.3.2生態阻力模型

生態資源具有向外擴展的能力[14]。在擴張過程中，由于土地利用和景觀的差異，會遇到生態障礙，受到生態阻力的影響。荷蘭生態學家KNAAPEN提出了最小累積阻力(MCR)[15]，本研究根據研究區性質對MCR模型進行改進。本文認為不同源地具有不同的生態能量，它們輻射和傳遞能量的能力不僅與其面積有關，而且與源的類型和形狀有關。目前，一些學者對此進行了研究和討論[16]。生態源的大小將對物質交換或能量流產生強烈影響，其邊緣也會影響物質、能量和未來的物種分布。因此本文增加了形狀指數(Li)和NDVI、MNDWI來描述生態資源的特征。由于研究區位于干旱地區，藍源(河流、湖泊)水平高于綠源(森林、草地)水平，本文提出了生態源地質量的概念，并改進了MCR模型。生態源地質量的定義是任何生態源地都具有向外輻射，傳遞能量或接收能量的特性，其生態源地質量與生態源地的面積、形狀、屬性有關。

(1)

式中Ei——生態源地i的周長

Ai——生態源地i的面積

(2)

式中mi——生態源地i的生態源地質量

本研究將基質阻力分為兩類：①建設用地集中區和農田等生態屏障。②地形坡度、植被覆蓋度、土地利用類型、地下水埋深、土壤干旱程度等生態障礙。一般來說，當地表變得陡峭時，植被覆蓋度和土壤干旱指數會降低，地下水埋藏得更深，阻力系數會越大，這將阻止生態資源的擴大。建設用地集中區和耕地是人類活動的主要集中區，土地利用方式難以改變。通過對研究區域的分析，改進的MCR模型更適合于研究區域的生態阻力模擬，改進的公式為

(3)

式中VMCR-m——生態源地擴展的最小累積生態阻力

Dij——生態源地i和j的空間距離

Ri——生態源地i運動過程中的阻力系數

1.3.3生態引力模型

生態引力模型是研究兩個物體或空間相互作用的模型，源于牛頓萬有引力定律[17]。19世紀，學者研究城市系統的相互作用時將引力模型引入到地理空間的研究中。后來，地理學家對引力模型進行了廣泛的研究，并用于分析城市之間貿易的相互作用，以及國家之間貿易交流形式的研究。2010年，KEUM[18]證明了引力模型在國際貿易和旅游業的可用性。

本文嘗試將引力模型引入景觀生態領域，計算生態源地之間的生態引力。生態引力與最小累積耗費距離的平方成反比，并且與生態源地質量成正比。生態引力模型為

(4)

式中Fij——生態源地i和生態源地j之間的生態引力

αi、αj——生態源地i、j的類型

G——常數，通常為1

mj——生態源地j的生態源地質量

1.3.4生態網絡性能分析

本文選擇復雜網絡分析中衡量網絡拓撲結構和節點重要性的6個圖論指標揭示各發展策略下FG生態網絡的連接特征，為FG生態網絡的優化提供重要參考[19]。從網絡整體性與節點的重要性角度對生態網絡的性能進行分析，其中連通性c(G)[20]、核數k(G)[21]、網絡密度d(G)[22]、平均節點連接數a(G)4個指標評價網絡的整體特征，而介數、PageRank值則評價節點的的重要性程度[23]。

表1 生態網絡性能分析指標Tab.1 Indicators for ecological network performance analysis

1.3.5骨架廊道識別

克魯斯卡爾(Kruskal)算法是骨架廊道識別的常見算法[24]，其主要思想是將初始最小生成樹邊為0，在迭代中，每次選擇最小代價邊(即本文中的生態廊道)加入最小生成樹集合，最終提取出骨架廊道[25]。記原網絡Graph中有v個節點，e條邊，將Graph中的頂點提取，形成新的網絡Graph(new)，此時該網絡中沒有邊，將原圖中的e條邊按照邊權從小到大排序，以邊權最小的邊開始，遍歷原網絡中的所有邊，直到v個節點都在同一個連通分量中。循環的條件為：若某條邊兩側的節點在Graph(new)中不在同一個連通分量中，將該邊增加到Graph(new)中。Kruskal算法的圖例描述如圖3所示，圖3a為原始的加權網絡，圖3b～3i為骨架廊道的提取過程。

圖3 Kruskal算法圖例描述Fig.3 Legend description of Kruskal algorithm

2 結果與分析

2.1 生態源地識別與提取

生態源地提取是生態網絡提取的第一步。根據景觀生態學中的“源-匯”理論:生態源地是多種生態服務的源頭，具備景觀連續性與完整性，對于維持生態系統穩定具有重要作用。選取對烏蘭布和沙漠東北緣具有重要生態意義的林地、草地、水體3種景觀類型斑塊，通過各斑塊面積、NDVI平均值、MNDWI平均值及斑塊形狀指數篩選生態源地。利用ArcMap軟件對各斑塊面積進行統計，利用Zonal Statics工具提取NDVI平均值、MNDWI平均值，采用Fragstats軟件計算斑塊形狀指數。通過熵值法確定斑塊面積、斑塊NDVI平均值、MNDWI平均值及斑塊形狀指數的權重，綜合評價各生態景觀斑塊的生態功能重要性程度，根據各斑塊生態功能重要性程度的統計分布特征，選取重要性排序中從高到低占據重要性總和前60%的生態用地斑塊作為生態源地。最終獲得研究區生態源地共481個，細碎及孤立的斑塊在篩選過程中被大量去除，景觀單元中的核心斑塊得到了保留，提取結果如圖4所示。

圖4 研究區生態源地分布Fig.4 Distribution of ecological sources in study area

其中水體型源地共37個，面積占所有生態用地斑塊的25.47%，水體型源地主要分布在研究區東部，西部區域也有零星水體分布。黃河是研究區最大面積的水域，黃河、奈倫湖等面積較大的水域均由黃河水引流形成。林草地型源地與水體型源地相比更加破碎，生態源地總數為439個，占全部生態源地的91.26%，而面積僅占全部源地的74.52%，林草地分散分布在研究區中部農田集中區、沙漠交錯地帶、狼山南麓以及黃河東岸，其他區域也有少量分布。

2.2 生態阻力模型構建

景觀生態流從源向外流動需要克服阻力，受多種因素的影響。研究區的生態環境特征由12個因素評估，包括地形因子(海拔、坡度)、植被覆蓋因子、水文因子(MNDWI、TVDI、地下水埋深)、景觀因子和密度因子(居民點密度、路網密度、水網密度)等。每個因子的阻力根據評估系統設定，如表2所示。

此外，發展策略將影響某些因素的阻力，并最終影響研究區域中基質的最小累積消耗阻力。例如，植被覆蓋度越高，水網越近，生態保護模式中的生態阻力越小，而經濟發展模式則相反。純生態保護模式((1.0,0)模式)和純經濟發展模式((0,1.0)模式)的生態阻力如圖5所示。在(1.0,0)模式下，基質的生態阻力為17～32，從東北向西南逐漸增大。阻力較低的地區主要位于東北部，由于黃河附近水資源豐富；而阻力較大的地區位于烏蘭布和沙漠、狼山和黃河的東岸，這里是庫布其沙漠的西緣。在這種發展模式下，嚴格限制經濟發展，優先考慮生態保護的發展策略，使得水資源和林草地得到較好的保護。研究區的生態源地及其周邊區域的生態阻力都較低。在(0,1.0)模式中，水資源主要用于農業生產和工業制造，而很少用于生態建設和保護，這使得在東北地區引起較高的基質阻力。在中部的烏蘭布和沙漠地區以及狼山的前沖積扇中，基質阻力相對較低。以上情況表明，在現狀的基礎上大力發展經濟將嚴重阻礙生態流在源地之間的流動，影響生態的發展。

最小累積生態阻力是生態流從一個生態源地到另一個生態源地，克服生態阻力面的最小累積成本。在改進的MCR模型的基礎上，采用成本距離模型計算不同生態源地之間的累積生態阻力。(1.0,0)模式和(0,1.0)模式的最小累積生態阻力如圖6所示。由圖6a可知，在(0,1.0)模式中，累計生態阻力最大值是380 727，最小值是28 459。累計阻力較大的區域主要位于南部的烏蘭布和沙漠。一方面是因為該區域生態環境條件較為惡劣，地下水埋深、水體密度、NDVI等多因子的阻力均較高。另一方面，該區域生態景觀斑塊較少，且與現存的生態用地距離較遠，生態阻力的累積作用較為明顯。研究區北部同樣出現了累積阻力較大的區域，盡管該位置生態條件較好，但這個區域是耕地，為生態屏障，阻力較大。且在該發展模式中，農業生產中大量從地下水、地表水中使用水資源，導致生態源地的供水水量較少，這也導致了在(0,1.0)發展模式中基質有比較大的生態累積阻力。這表明，在研究區內大力發展經濟將嚴重阻礙生態流在源地之間的流動，嚴重影響生態源的發展。由圖6b可知，在(1.0,0)發展模式中，累積生態阻力最大值是346 774，最小值是17 467。由于生態條件惡劣，累積阻力較大的區域同樣主要位于南部的烏蘭布和沙漠。不同于(0，1.0)模式，在該模式中生態阻力較小的區域顯著增加。

表2 生態阻力評價指標體系Tab.2 Ecological resistance evaluation index system

圖5 (0,1.0)模式和(1.0,0)模式的生態阻力Fig.5 Ecological resistance in (0, 1.0) mode and (1.0, 0) mode

圖6 (0,1.0)模式和(1.0,0)模式的最小累積生態阻力Fig.6 Minimum cumulative ecological resistance in (0, 1.0) mode and (1.0, 0) mode

2.3 生態引力模型

通過將生態源地的質量和其間的最小累積耗費距離代入生態引力模型，可以計算出任何生態源地之間的生態引力。雖然生態源地會受到基質的阻礙作用，但它們之間仍然可能通過生態廊道形成聯系。生態引力越大，生態源地之間的聯系越緊密，它們之間的生態流傳遞能量、物質就越頻繁。而且當生態引力足夠大時，生態源地會通過不斷吸收和合并外部資源的方式，擴大其自身規模。但如果生態引力很小，在生態源地之間形成生態廊道的可能性不大，除非基質發生變化，否則它們之間將不會形成聯系。

以1號生態源地為例，表3顯示了1號生態源地與其他生態源地之間的生態引力。由于生態引力在1號生態源地和342號生態源地或233號生態源地之間引力最大，因此在它們之間形成生態廊道的可能性最高。相比之下，由于在1號生態源地和35號生態源地或66號生態源地之間生態引力最小，形成生態廊道最困難。

表3 生態源地之間的生態引力(源地1)Tab.3 Ecological gravitation between ecological sources (source 1)

2.4 生態網絡多場景模擬

圖7 (0.1,0.9)模式到(0.9,0.1)模式的最小累積生態阻力Fig.7 Minimum cumulative ecological resistance from (0.1, 0.9) mode to (0.9, 0.1) mode

生態引力由生態源地的性質和其間的引力半徑決定，并且不隨著基質的變化而波動。相比之下，源地之間的最小累積生態阻力與基質密切相關。通過使用柵格計算器模塊獲得本研究中其他9個發展策略的累積生態阻力。如圖7所示，9種模式的累積生態阻力范圍為27 743～379 056、27 064～377 915、26 132～375 922、25 875～373 169、23 509～270 376、22 205～368 743、21 413～363 205、19 976～359 492和18 603～352 301。累積阻力大的區域位于南部的烏蘭布和沙漠以及北部河套平原附近的農業開發區。隨著經濟發展的重要性增加，生態保護的重要性逐漸降低，生態積累阻力較小的灰色地帶逐漸減少，這將使生態資源地之間能源和物質交換的可能性變小。

通過模擬不同開發方案中基質的最小累積阻力，比較生態源地之間的生態阻力和引力，確定其是否能形成生態廊道，然后利用Python腳本語言模擬生態引力與生態源阻力之間的相互作用，構建FG生態網絡。模擬結果如圖8所示。

圖8中FG生態網絡分布不均勻。例如，北方生態源地之間的生態廊道數量較多，生態網絡較為復雜，而南方的FG生態網絡則具有相反的特征。(1.0,0)模式的FG生態網絡比(0,1.0)模式更復雜，(1.0,0)模式下生態源地之間的聯系更緊密。在(1.0,0)模式向(0,1.0)模式的轉變中，發展策略中經濟發展的比例增加，一些地區的生態阻力逐漸增加。生態源地之間的最小累積生態阻力大于生態引力，這導致生態廊道破裂并導致FG生態網絡的破壞。

2.5 生態網絡性能分析

采用Matlab計算不同情景下FG生態網絡的連通性，核心數和密度、平均節點連接數以及介數和PageRank值。計算結果如表4所示。(1.0,0)模式是目前研究區域的發展模式，其網絡連通性最大，為329，占網絡生態資源總量的68.24%。結果表明，去除329個生態源地可以將(1.0,0)模型的FG生態網絡轉化為平凡圖，從而完全喪失網絡的生態功能。

圖8 (0,1.0)模式和(1.0,0)模式中的FG生態網絡Fig.8 Forest-grass ecological network in (0, 1.0) mode and (1.0, 0) mode

表4 11個模擬場景下FG生態網絡的結構特征Tab.4 Structure characteristics of forest-grass ecological network in 11 simulated scenarios

從網絡的整體特征來看，隨著發展策略中經濟發展的比重增加、生態保護的比重減小，FG生態網絡的連通性逐漸降低，并且FG生態網絡的穩定性隨著連通性的降低而惡化。在純粹經濟發展的(0,1.0)模式下，FG生態網絡的連通性最低。僅需要去除265個生態資源可以將(0,1.0)模式的FG生態網絡轉化為平凡圖。核數的變化趨勢與連通性的相同。隨著經濟發展的比重逐漸增大，FG生態網絡的核數減少，同時網絡的連通性和穩健性也在逐漸減弱。網絡平均節點連接數的變化趨勢與核數和連通性相同。就網絡密度而言，網絡密度在從(1.0,0)到(0,1.0)模式的演變過程中先增加后減小。FG生態網絡的密度在(0.9,0.1)模式中最大。從節點的重要性程度來看，節點的最大介數在(0.9,0.1)模式中最大，隨著發展策略中經濟發展的比重增加而持續減小。隨著發展策略中經濟發展比重的增加、生態保護的比重減小，節點的最大PageRank值逐漸降低。

2.6 關鍵節點及骨架廊道識別

提取研究區生態節點的介數、PageRank值，統計分布中前10%的生態節點作為關鍵生態節點，利用最小生成樹算法，提取研究區潛在生態廊道中的骨架廊道，結果如圖9所示。

圖9 研究區關鍵生態節點及骨架廊道Fig.9 Key ecological nodes and skeleton corridors of study area

研究區識別出36個關鍵生態節點。關鍵生態節點是研究區生態網絡的核心組成部分，其生態功能的完好，有利于區域生態物質、能量和信息的交換，對于維持區域物種多樣性、保持生態環境穩定、防止沙漠擴張具有至關重要的作用，應著重加強對關鍵生態節點的生態建設與保護，實行較嚴格的生態管控措施，防止關鍵生態節點被破壞。絕大部分關鍵生態節點位于研究區中部，如扎木溫多爾、嘎咕塔拉、阿茨嘎查以及納林湖附近，少量位于東部如奈倫湖國家濕地公園及黃河附近，南部僅有個別點位分布，位于哈拉特爾瑞輝附近。

經提取的骨架廊道基本保留了研究區生態廊道的骨架部分，細碎部分在提取過程中被舍棄。對研究區骨架廊道的長度分布進行統計可知，骨架廊道總長度為1 236.89 km。因此應對骨架廊道進行重點建設與保護，尤其是廊道110、160、235、297、304、311、457、487。

3 結論

(1)考慮了生態網絡中的相互作用，并使用引力模型來量化生態源地之間的相互作用。基于引力模型和MCR模型，構建了沙漠-綠洲交錯帶FG生態網絡多場景模擬模型，用于提取不同情景下的生態網絡。

(2)從(1.0,0)模式到(0,1.0)模式，隨著經濟發展比例的增加和生態保護比例的下降，FG生態網絡逐漸被破壞。在(0,1.0)模型中，FG生態網絡的結構被破壞，FG生態網絡的連通性、穩定性和核數達到最低水平。這表明，在目前的自然條件下，干旱地區的大規模經濟發展將對FG生態網絡造成巨大破壞。研究發現(0.9,0.1)模式的FG生態網絡更復雜，網絡的穩定性和連通性指標更大。然而，隨著經濟的發展，FG生態網絡的結構遭到越來越嚴重的破壞。這表明研究區可能有經濟發展空間，但空間不大。在干旱地區，經濟發展與生態破壞之間存在平衡。應根據當地條件發展經濟，避免過度開采造成的生態惡化，加劇荒漠化。

(3)根據潛在生態網絡拓撲結構特征，識別出36個關鍵生態節點；利用最小生成樹算法，識別出研究區骨架廊道總長度為1 236.89 km。對于關鍵生態節點的生態建設，應實行較嚴格的生態管控措施，防止因關鍵生態節點被破壞導致的區域生態惡化。骨架廊道是研究區生態廊道的主體部分，是網絡其他部分連通的基本保障，因此應對骨架廊道進行重點建設與保護。