商洛市景觀格局粒度效應響應及景觀格局分析

房舒 趙培

摘 要：基于景觀格局粒度響應能更科學準確地研究區域景觀格局特征，為促進地區可持續發展、制定環境修復政策提供理論依據。以秦嶺南麓商洛市為研究對象，基于2020年商洛市土地利用狀況，探討商洛市景觀粒度效應響應，確定最適宜商洛市景觀格局研究的粒度。并在此粒度下，研究商洛市景觀格局特征。結果表明，商洛市景觀級別和斑塊類型級別的最適宜研究粒度都為90 m。商洛市整體生態環境較好，但仍存在斑塊破碎化程度高，有一定脆弱性的問題。未來發展中應更加注意提升林地和水域斑塊的穩定性，實現經濟發展與生態保護的雙贏。

關鍵詞：景觀格局指數；粒度效應；最適宜粒度；景觀格局特征

中圖分類號：P901? ? 文獻標識碼：Ａ? ? ?文章編號：1674-0033（2023）02-0022-09

引用格式：房舒，趙培.商洛市景觀格局粒度效應響應及景觀格局分析[J].商洛學院學報，2023，37（2）：22-30.

Abstract： Landscape pattern based on grain size effect can more scientifically and accurately help study the characteristics of regional landscape pattern， and provide a theoretical basis for promoting regional sustainable development and formulating environmental restoration policies. Taking Shangluo City in the south of Qinling Mountains as the research object， based on the land use status of Shangluo City in 2020， the response of landscape grain size effect in Shangluo City was discussed to determine the most suitable granularity for the study of landscape pattern in Shangluo City. And under this granularity， the characteristics of landscape pattern in Shangluo City were studied. The results showed that the most suitable research grain size of landscape level and patch type level in Shangluo City was 90 m. The overall ecological environment of Shangluo City was good， but there were still problems of high degree of patch fragmentation and certain vulnerability. In the future， we should pay attention to enhance the stability of forest land and water patches， and achieve a win-win between economic development and environmental protection.

Key words： landscape pattern index; grain size effect; most suitable grain size; landscape pattern characteristics

景觀格局指數能定量化地探究景觀格局特征，為區域的生態安全和可持續發展提供基礎數據，在土地利用方式的提升、景觀格局的優化等多個領域廣泛應用[1-3]。尺度是決定景觀格局研究科學與否的關鍵，景觀格局指數對空間尺度的變化非常敏感，景觀格局的尺度效應是景觀格局特征推繹的基礎[4-5]。尺度包括時間和空間尺度，空間尺度包含空間幅度和粒度[6]。在研究過程中，研究區域和研究年限確定后，時間尺度和空間幅度也隨之確定，尺度效應只需關注粒度[7]。景觀格局粒度效應響應，是指隨空間數據柵格單元大小改變，景觀格局指數發生改變的過程，目的是通過臨界點或轉折點確定研究景觀格局的最佳粒度[8]。已有研究者針對不同的行政區、經濟區、流域、山地、河谷等進行了景觀格局效應響應研究[1-3，9]。基于景觀格局粒度響應，在最適宜粒度下研究區域景觀格局特征可為促進地區可持續發展、制定環境修復政策提供理論依據和技術支持[10]。

秦嶺是我國南北氣候的分界線和重要的生態屏障[11]。2019年9月陜西省修訂的《秦嶺生態環境保護條例》中將保護秦嶺生態環境，實現經濟社會可持續發展提上新的高度[12]。習近平總書記來陜考察時多次強調秦嶺生態環境問題。然而多年來秦嶺生態環境保護和修復研究主要集中在秦嶺北麓[13]。秦嶺南麓資源豐富、環境優美，經濟發展和生態環境建設之間的矛盾突出[14]。研究發現秦嶺東南部生態風險較高，且有繼續增高的趨勢[15]。雖然早有學者在2008年就提出要及早關注秦嶺南麓生態環境的保護，但目前為止針對秦嶺南麓生態環境評價及修復的研究相對較少。

商洛市地處陜西省東南部，是全國生態環境建設重點試點示范區。隨著人口增長、經濟發展，商洛市生態風險問題日益突出[16]。保護水源地生態環境，確保受水區水質水量安全，實現人與自然和諧共處，一直是該區的熱點問題，但是針對商洛市開展的相關研究較少。因此，本文以秦嶺南麓的商洛市為研究對象，探討商洛市景觀粒度效應響應，在最適宜的粒度下，研究商洛市景觀級別和斑塊類型級別的景觀格局特征，為秦嶺南麓的景觀格局優化和可持續發展提供依據。

1? 研究區概況

商洛市（108°34′～111°01′E，33°02′～34°24′N）地處陜西省東南部，是處于秦嶺南麓的城市，也是長江水系漢江支流丹江的源頭。商洛市是南水北調中線工程的重要水源地，是典型的山區水源地[17]。

商洛市的地形有盆地、低山和中山三大類，是以中低山為主的土石山區，土壤主要以褐土、棕壤、黃棕壤為主。氣候主要為南暖溫帶向北亞熱帶過渡的大陸性山地季風氣候，四季分明，雨熱同期，是陜西省林業大市、重點林區。年均降水量為710～930 mm，年均氣溫為7.8～13.9 ℃，日照約為1 860～2 130 h，無霜期為210 d。

2? 數據來源與研究方法

2.1 數據來源

下載商洛市2020年夏季Landsat 8遙感影像，進行預處理和增強處理后，在ENVI5.3中進行非監督ISODATA分類，將商洛市土地利用類型劃分為耕地、林地、草地、水域、居民用地、未利用土地六類，分類后總精度大于90%。

2.2 研究方法

研究路線如圖1所示。首先從15～450 m選取了22個不同的粒度，其中15～210 m粒度之間的間隔設為15 m，210～450 m粒度之間的間隔設為30 m。然后依據文獻[9]及研究區的土地利用特征，選擇了景觀和斑塊類型兩個級別中較為常用且能代表不同指標類別的景觀指數。在景觀級別上選擇了8個景觀格局指數（見表1），在斑塊類型級別上選擇了6個景觀格局指數（見表2）。最后運用景觀格局指數計算軟件FRAGSTATS 4.2，計算了不同粒度兩個級別的景觀格局指數，并制作粒度響應曲線，結合SPSS回歸分析曲線估計結果，得到適宜商洛市景觀格局研究的最適宜粒度，并在此粒度下計算和分析商洛市景觀格局特征。

3? 結果與分析

3.1 景觀級別的景觀指數粒度選擇

2020年商洛市8個景觀級別的景觀指數的粒度響應曲線如圖2所示。由圖2可以看出，8個景觀級別的景觀指數的響應曲線有三種趨勢。第一種，曲線無明顯規律，包括面積—邊緣指標的LPI（平穩—波動上升）、聚散性指標的PD（下降—上升—下降）和IJI（平穩上升—波動下降）及多樣性指標的SHEI（不斷波動）。第二種，景觀指數隨粒度增大呈上升趨勢，包括形狀指標的PAFRAC，其粒度響應曲線在粒度為45 m處出現拐點，在360 m后出現小幅度上下波動。第三種，景觀指數隨粒度增大呈下降趨勢，包括SPLIT、AI、LSI，均屬于聚散性指標，其中SPLIT較為特別，屬于波動下降。AI隨粒度增大呈平穩下降，下降趨勢無明顯變化。LSI在粒度為210 m處出現拐點，下降趨勢變大。

將景觀級別的景觀指數與粒度的關系進行模擬，結果如表3。大部分景觀級別的景觀指數的尺度效應可以較好地用函數模擬。其中多樣性指標的SHEI沒有合適的函數模擬，PAFRAC適合用冪函數模擬，而其他指數均適合用三次函數模擬。聚散性指標IJI與面積—邊緣指標的LPI的擬合R2<0.9。而形狀指標的PAFRAC與聚散性指標AI、LSI的擬合均有R2>0.97。

結合景觀格局指數粒度效應曲線和曲線擬合的結果，得到商洛市景觀級別的景觀指數第一尺度域為75～105 m，最適宜研究粒度為90 m（見表4）。

3.2 斑塊類型級別的景觀指數粒度選擇

2020年商洛市6個斑塊級別的景觀指數的粒度響應曲線如圖3所示。

由圖3可以看出，6個斑塊級別的景觀指數的響應曲線有四種趨勢。第一種，景觀指數曲線無明顯規律，包括面積—邊緣指標的LPI和聚散性指標的PD、NP。LPI的粒度響應圖中，居民地、水域、未利用地的LPI幾乎無變化；草地的LPI在15～210 m呈上升趨勢，但之后有小幅下降，LPI在210 m之后出現較大波動；林地的LPI在粒度為15～105 m小幅波動，在105～150 m波動變大，之后為波動上升；耕地的LPI總體呈現波動下降的趨勢。PD的粒度響應圖中，未利用地的PD幾乎無變化；水域的PD在15～195 m呈略微上升趨勢，之后呈現略微下降趨勢；居民地、草地、林地的PD均呈現隨粒度的增大而下降的趨勢；耕地的PD在粒度為15～45 m呈現略微下降趨勢，粒度在45～210 m呈現上升趨勢且上升幅度先小后大，在210 m之后呈下降趨勢。NP的粒度響應圖與PD變化趨勢相似。第二種，景觀指數隨粒度增大而增大，為形狀指標的PAFRAC。居民地、草地、林地、耕地的PAFRAC均有隨粒度的增大而上升的趨勢，其中居民地的上升趨勢較為明顯；但對于未利用地來說，沒有明顯的上升或下降趨勢，存在較大波動；而對于水域來說，則存在整體下降趨勢。第三種，景觀指數隨粒度增大呈下降趨勢，為聚散性指標AI。在其粒度響應圖中，所有土地利用類型均存在明顯的景觀指數隨粒度增大而下降的趨勢，其中草地、林地的下降趨勢較小，水域的下降趨勢較大，未利用地在粒度為195 m之后呈現出波動下降。第四種，景觀指數基本不隨粒度增大變化，為面積—邊緣指標的PLAND。其中草地、林地的景觀指數均在粒度為420 m處出現一不明顯的拐點；耕地的景觀指數整體上存在較小的波動。

將斑塊類型級別景觀指數與粒度的關系進行模擬，結果如表5所示。由表5可以看出，大部分斑塊級別的景觀指數的尺度效應可以較好地用函數模擬，其中LPI的未利用地沒有合適的函數模擬，林地雖然可以用指數函數進行模擬但擬合效果差；草地、林地、耕地、水域都可以用三次函數進行較好地模擬。PLAND景觀指數所有斑塊類型均沒有合適的函數模擬，粒度效應不明顯。PAFRAC景觀指數中，居民地、林地、水域均可以用三次函數進行模擬，但水域的擬合程度較低；草地可以用對數函數進行較好地模擬；耕地適宜用S模型模擬；未利用地沒有合適的函數模擬。PD的景觀指數除林地適合用二次函數進行模擬外，其余斑塊類型均適合用三次函數進行模擬。NP的函數擬合結果與PD相似。AI的六種斑塊類型均可用三次函數進行模擬，且擬合效果好，R2均在0.99以上。

結合景觀格局指數粒度效應曲線和曲線擬合的結果，得到商洛市斑塊類型級別景觀指數第一尺度域為75～105 m，最適宜研究粒度為90 m（見表6）。

3.3 商洛市景觀格局特征

商洛市的主要土地利用類型是草地，占總面積的40.83%，接著是林地和耕地，分別占總面積的35.75%和22.45%。居民用地、水域和未利用地的比例都不足1%。最適宜商洛市景觀格局研究的粒度為90 m。90 m粒度下，商洛市景觀級別、斑塊類型級別的景觀指數如表7、表8。

景觀級別的景觀格局指數能說明商洛市景觀的整體水平，見表7。LPI的計算結果小于10，說明商洛市景觀基質不突出，斑塊較為破碎。PAFRAC的計算結果為1.57，商洛市景觀形狀復雜程度不強。PD的計算結果為0.52，說明商洛市景觀較統一。SPLIT的計算結果為96.27，說明商洛市景觀細化程度較高。IJI的計算結果小于50，說明商洛市斑塊分布聚集程度不明顯。AI的計算結果大于80，說明商洛市景觀斑塊間聯系緊密。LSI的計算結果接近150，說明商洛市景觀同類斑塊聚集程度較好。SHEI的計算結果大于0.6，說明商洛市景觀類型較為豐富。

斑塊級別的景觀格局指數用于比較商洛市不同土地利用類型的景觀格局特征，見表8。從LPI來看，草地、耕地、林地斑塊中含有優勢度大的斑塊，居民地、水域、未利用地斑塊相對更破碎。從PLAND來看，草地是商洛市的主要景觀類型，占地面積比例大于40%，林地、耕地景觀占地面積比例均大于20%，居民地、水域、未利用地斑塊所占面積比例都較小，均不足1%。從PAFRAC來看，水域的斑塊計算值最大，形狀最沒有規律，居民地的斑塊計算值最小，形狀最為規則，其他斑塊形狀不規則。PD和NP的變化趨勢與LPI保持一致。從AI來看，商洛市AI最高的為林地，聚集性最強；最低的為水域，較為分散。

4? 討論與結論

本文選取了22個粒度，8個景觀級別、6個斑塊類型級別的景觀指數，對商洛市景觀格局指數粒度效應進行研究，發現在景觀級別和斑塊類型級別最適宜研究區的景觀粒度都為90 m。這與陳永林等[18]在長沙市的景觀粒度效應研究中得出的結果一致。同時，翟俊等[19]認為流域景觀格局分析的最佳空間粒度為90 m。楊馗等[20]認為90 m的粒度能更好地體現空間異質性。

整體來看，商洛一區六縣均位于中國南水北調中線水源地生態環境保護規劃區內[17]，而商洛四面環山，人口少、耕地少、企業少、經濟增長緩慢，城市土地集約利用、土地投資、經濟水平在陜西省10個城市中都處于最低水平[21]。商洛的土地利用主要是草地和林地。雖然商洛整體環境較好，植被碳匯能力超過500 g/（cm2·a）[22]，但是由于商洛處于山區[23]，仍然存在景觀斑塊細碎化、破碎化，脆弱化程度高的問題，具有一定的景觀生態風險。

商洛市森林覆蓋率高，但天然次生林的比例較大，進行工程管護的森林面積達160多萬公頃[24-25]。近年來商洛建成區植被覆蓋度有所改善[25]，林地景觀雖然存在一定優勢，但穩定性不足。商洛市山陽縣礦產資源豐富，2005年前因露天開采森林資源遭到破壞，景觀脆弱性增加；但2005年后，礦區環境有所改善，生態風險略有降低[15]。2015年到2019年，商洛的資源利用效率呈上升趨勢[26]。商洛市雖然為南水北調中線主要供水城市，但其水域斑塊仍存在分散、破碎的問題。為保證供水質量，南水北調中線工程實施以后商洛市開展了水污染治理、水環境保護工作來提高丹江庫區及上游水質質量[27]，但對由此帶來的土地結構改變、區域環境響應缺乏系統性的認識和研究，同時，急需治理的水土流失面積仍有232萬公頃[28]。可以看出，提高商洛市景觀格局的穩定性，主要還是要從林地斑塊和水域斑塊的改善入手。

依托得天獨厚的氣候資源、優美的自然環境、無公害農產品等，商洛市未來的發展目標是打造“中國康養之都”。未來自然空間景觀結構的穩定和提升將是商洛城市發展的重要基礎，同時商洛經濟發展也離不開生態環境良好帶來的優勢[27]。因此商洛市未來仍需繼續優化景觀結構，尋求經濟發展與生態保護相協調的發展模式。

參考文獻：

[1]? 華琳，黃志霖，馬良，等.三峽庫區低山丘陵區多尺度景觀指數響應及適宜粒度[J].生態學報，2022，42（11）：4703-4717.

[2]? 周海菊，劉小英，胡靚達，等.基于最佳分析粒度的廣西北部灣經濟區景觀格局動態變化分析[J].生態與農村環境學報，2022，38（5）：545-555.

[3]? 任梅，王志杰，王志泰，等.黔中喀斯特山地城市景觀格局指數粒度效應——以安順市為例[J].生態學雜志，2018，37（10）：3137-3145.

[4]? 陳雅如，肖文發，滕明君，等.三峽庫區景觀格局粒度效應及其對土地利用變化過程的響應[J].自然資源學報，2018，33（4）：588-599.

[5]? 葉均艷，王亞娟，劉小鵬，等.基于最佳粒度的生態移民安置區土地利用景觀格局分析——以寧夏紅寺堡區為例[J].水土保持研究，2018，25（2）：300-306，314.

[6]? 劉學錄.我國景觀生態學尺度效應研究進展[J].甘肅農業大學學報，2021，56（6）：1-9.

[7]? 頓耀龍，王軍，白中科，等.松嫩平原西部土地整理區景觀格局指數的粒度效應研究[J].水土保持研究，2014，21（5）：66-71.

[8]? 陸禹，佘濟云，羅改改，等.基于粒度反推法和GIS空間分析的景觀格局優化[J].生態學雜志，2018，37（2）：534-545.

[9]? FANG S， ZHAO Y H， HAN L， et al. Analysis of landscape patterns of arid valleys in China， based on grain size effect[J].Sustainability，2017，9（12）：1-20.

[10] LEUVEN R S， POUDEVIGNE I. Riverine landscape dynamics and ecological risk assessment[J].Freshwater Biology，2002，47（4）：845-865.

[11] 張曉鵬，郭永峰.秦嶺北麓臨渭段生態環境保護分析[J].防護林科技，2019（9）：65-66，86.

[12] 梁宏賢.全面貫徹習近平生態文明思想 以立法保障秦嶺寧靜和諧美麗——修訂《陜西省秦嶺生態環境保護條例》的幾點體會[J].中國人大，2019（24）：51-52.

[13] 李紹文.淺析秦嶺北麓生態系統保護與修復[J].陜西水利，2019（2）：103-105，110.

[14] 李懷印，徐剛，肖穎.及早重視保護秦嶺南麓生態環境[J].陜西林業，2008（2）：28-30.

[15] LIU D， CHEN H， ZHANG H， et al. Spatiotemporal evolution of landscape ecological risk based on geomorphological regionalization during 1980—2017： A case study of Shaanxi Province， China[J]. Sustainability，2020，12（3）：1-16.

[16] 楊維鴿，趙培，馬垚垚，等.山區城市空間擴張特征及影響因素分析[J].甘肅科學學報，2019，31（6）：129-135.

[17] 孫丕苓，楊海娟，劉慶果.南水北調重要水源地的土地生態安全動態研究——以陜西省商洛市為例[J].自然資源學報，2012，27（9）：1520-1530.

[18] 陳永林，謝炳庚，李曉青.長沙市土地利用格局變化的空間粒度效應[J].地理科學，2016，36（4）：564-570.

[19] 翟俊，侯鵬，趙志平，等.青海湖流域景觀格局空間粒度效應分析[J].國土資源遙感，2018，30（3）：159-166.

[20] 楊馗，信桂新，蔣好雨，等.基于最佳尺度的景觀生態風險時空變化研究：以重慶市江津區為例[J].生態與農村環境學報，2021，37（5）：576-586.

[21] ZHOU J， CAO X S. What is the policy improvement of China's land consolidation？ Evidence from completed land consolidation projects in Shaanxi Province [J]. Land Use Policy，2020，99：1-13.

[22] WEI X D， YANG J， LUO P P， et al. Assessment of the variation and influencing factors of vegetation NPP and carbon sink capacity under different natural conditions[J].Ecological Indicators，2022，138：1-15.

[23] BIN L. Spatial development explorations of middle-and small-sized mountainous cities： A case study of Shangluo City， China[J].Journal of Landscape Research，2012，4（12）：22-24.

[24] 張春玲.陜西省漢江、丹江非點源污染及控制對策[J].西北水資源與水工程，2002（1）：18-25.

[25] CUI L， ZHAO Y H， LIU J C， et al. Vegetation coverage prediction for the Qinling Mountains using the CA–Markov model[J]. ISPRS International Journal of Geo-Information，2021，10（10）：1-17.

[26] YU X H， LIU Y， ZHANG Z Y， et al. Urban spatial structure features in Qinling mountain area based on ecological network analysis-case study of Shangluo City[J]. Alexandria Engineering Journal，2022，61（12）：12829-12845.

[27] 劉建林，張雁，高英，等.南水北調（中線）工程水源地商洛水資源可持續利用研究[J].陜西農業科學，2010，56（1）：50-52.

[28] 張雁.山區水源地建設的生態環境效應評價與保護研究[D].西安：西安理工大學，2018：1-50.