松嫩平原西部景觀自然性評價及動態變化

林奕伶, 馬長迪, 劉吉平

(吉林師范大學 資源保護與利用研究所, 吉林 四平 136000)

目前，對“自然性”這一概念還沒有明確的定義，可將其理解為未受人類干擾的生態系統狀態，也可理解為當前狀態與原始狀態的差異[1]。國內外針對“自然性”這一特性已開展了部分研究，如馮澤深[2]、張棟[3]、張敏[4]、占昕[5]等利用層次分析法，分別對雁棲河、潮河、河北山區森林、閩江河口濕地保護區等地的自然性進行定量研究。前人的研究方法主要為層次分析法，研究對象多數為河流和森林等生態系統，而對景觀自然性的研究較少。景觀的自然性及其評價是景觀管理和自然保護的重要課題[6]，也可作為生態系統恢復和景觀規劃的參考[7-10]。如何選擇科學的研究方法、構建合理的評價指標體系對景觀自然性進行評價是目前需要解決的問題。

層次分析法是一種定性與定量結合的賦權方法，但存在主觀性過強的缺點。為了兼顧主觀賦權與客觀賦權，本文在層次分析法的基礎上引用了信息熵模型，利用綜合評價法進行權重賦值，可以在一定程度上減少因主觀成分過強或數據差異過大而產生的誤差。松嫩平原地區生態環境脆弱，抗干擾能力差。近年來人類對自然環境的影響逐漸加強，導致松嫩平原生態環境的進一步惡化，主要表現為土地沙化和濕地鹽漬化[11]。本文以松嫩平原西部為研究區域，對松嫩平原西部的景觀自然性及其動態變化進行研究，探究松嫩平原的景觀自然性可以作為人類開發、改造自然環境的依據和診斷，從而引導人類更加科學、有效地保護自然。

1 研究區域與方法

1.1 研究區概況

松嫩平原西部坐標為43°59′—48°33′E，121°37′—126°39′N(圖1)，總面積為1.01×105km2，氣候為中溫帶大陸性半濕潤半干旱季風氣候，地貌類型以山前傾斜平原、低平原以及部分河谷平原為主[12]。主要土地利用類型為森林、草地、濕地、農田，優勢景觀為農田景觀。土壤類型主要為黑土、黑鈣土、栗鈣土、草甸土、鹽堿土、風沙土以及沼澤土[13]。主要植被類型為虎尾草群落、堿茅群落、蘆葦群落以及羊草群落[14]，植被覆蓋率較高。主要城市包括白城市、松原市、大慶市、齊齊哈爾市。人口總數為1 731萬人(2015年)，人口密度較大。經濟發展水平較高，GDP總量為8 155億元(2015年)。產業以農業、制造業為主，是我國重要的商品糧基地。1981年至今依次設立了莫莫格、向海、扎龍、大布蘇狼牙壩、查干湖5個國家級自然保護區。

1.2 數據來源及處理過程

本文以探究松嫩平原西部不同區域景觀自然性狀態為目的，采用網格分析法，考慮到網格劃分過小時，會導致數據量過大，而網格劃分較大時會導致計算結果粗略。景觀尺度一般劃分為4個尺度等級，即1～106 km2為小尺度，106～1 010 km2為中尺度，1 010～1 012 km2為大尺度，大于1 012 km2為巨尺度。為兼顧科學性和可操作性，參考相關文獻[15]，利用ArcGIS 10.2軟件將研究區劃分為155個網格單元。網格大小為625 km2，研究區邊界地區的網格面積略大或略小于625 km2。本文的景觀尺度處于中尺度，能表示具體區域的變化，是松嫩平原西部較為適宜的景觀自然性評價尺度。

圖1 研究區位置

土地利用數據來源于1985年、2000年、2015年的Landsat TM/OLI遙感數據。首先利用ENVI軟件對遙感影像的5，4，3波段進行標準假彩色合成，并對合成的遙感影像進行前期處理。然后用1∶10萬地形圖對其進行糾正。最后對其進行目視解譯并提取耕地、林地、草地、濕地等土地利用類型。采用抽樣統計精度驗證法，在野外設立驗證點，并結合高分辨率遙感影像對解譯結果進行檢驗，最終解譯精度達到90.7%以上，基本可以滿足本文的研究需求。

道路數據來源于1∶50萬黑龍江省、吉林省公路基礎數據庫的道路要素圖。人口、GDP數據來源于各地區統計年鑒。利用SPSS軟件對各縣、市的人口、GDP數據進行統計，在ArcGIS 10.2中將其與劃分的155個網格進行疊加、計算，最終得到1985—2015年松嫩平原西部各網格的人口密度與GDP數據。

1.3 研究方法

1.3.1 構建評價指標體系 在參考國內外自然性評價的基礎上，以科學性、整體性、代表性、可量化性、層次性為主要原則，構建了以景觀指數、土地利用類型、人類干擾為主體的綜合評價體系。景觀指數是對景觀格局進行定量分析的重要工具，研究景觀變化必然要用到相應的景觀指數。土地利用類型也是衡量景觀自然狀態的重要依據，土地利用和土地覆蓋變化是影響景觀結構、功能及動態的最普遍的主導因素之一，同時也是景觀生態學和全球生態學中極重要和頗具挑戰性的研究領域之一[16-17]。由于人類活動對景觀影響的普遍性與深刻性，在作為人類生存環境的各種景觀中，人類活動對景觀演化起著主導作用[18]，因此在評價指標中加入了人類干擾相關因素，如人口、GDP等。

評價體系包括目標層(A)、準則層(B)、指標層(C)3個層次。

其中目標層為景觀自然性評價(A)，準則層包括壓力指標(B1)、狀態指標(B2)、響應指標(B3)。在PSR模型的基礎上，向下劃分出13個指標。指標層包括人口(C1)、GDP(C2)、道路密度(C3)、人為干擾度(C4)、景觀多樣性(C5)、景觀均勻度(C6)、破碎化程度(C7)、景觀優勢度(C8)、景觀形狀指數(C9)、景觀結構穩定性[19](C10)、連通性(C11)、NDVI指數(C12)、自然景觀比例(C13)。

1.3.2 評價指標權重的確定及評價模型的構建 本文采用主觀賦權法與客觀賦權法結合的綜合評價法。

在主觀權重的確定上，使用的方法是層次分析法(AHP)，本文利用yaahp 0.5.3軟件進行賦值，首先，構建層次結構模型。結構模型確立后構建判斷矩陣，采用專家打分法，對各指標的相對重要程度進行對比賦值，賦值方法為1～9標度法。通過層次單排序與層次總排序進行一致性檢驗，4個判斷矩陣CR均小于0.1，通過一致性檢驗。

同時，本文利用信息熵模型計算客觀權重。1948年，香農在熱力學研究的基礎上提出了信息熵公式，他認為，任何信息都存在冗余，冗余大小與信息中每個符號(數字、字母或單詞)的出現概率或者說不確定性有關。信息熵模型的基本思想為：評價因子差異程度的大小能夠反映該因子在整個綜合評價體系中的地位，而評價因子差異程度的大小可以利用信息熵進行度量[20]。

分別利用層次分析法和信息熵模型計算出評價指標的主觀權重、客觀權重。最后利用綜合評價法將二者進行結合，公式為

W=αw主觀權重+(1-α)w客觀權重

式中：W為綜合權重值。其中；0≤α≤1，本文中取α=0.7。

計算后得到的綜合權重見表1，其中權重值較高的指標有破碎化程度、景觀優勢度、自然景觀比例、GDP。

確定綜合權重后，利用多因子加權評價模型對松嫩平原西部景觀自然性進行分析，計算公式為：

式中：N為景觀自然性值；Wi為因子i的綜合權重值；Pi為標準化處理后第i項指標的評價值。

表1 綜合評價法所得權重值

2 結果與分析

2.1 松嫩平原西部景觀自然性評價

從整體上看，1985年、2000年、2015年松嫩平原西部的景觀自然性分別為0.313 7，0.278 1，0.269 8，呈不斷降低趨勢，1985—2000年下降率為11%，2000—2015年下降率為3%(圖2)。下降的主要原因是壓力指標的變化，1985—2000年、2000—2015年兩個時期的壓力指標分別上升了47%，9%(圖2)。而壓力指標中增長最快的是GDP指標，兩個時期的GDP增長率分別為749%和368%。其次為響應指標的變化，1985—2000年、2000—2015年兩個時期的響應指標分別下降了7.77%，1.30%。

圖2 1985-2015年各指標及景觀自然性值

從地區看，以生態保護為目的的自然保護區的景觀自然性整體上高于城市景觀自然性。1985—2015年，4個城市中僅有大慶市的景觀自然性呈上升狀態，其余3個城市的景觀自然性均有不同程度的下降(表2)。在5個保護區中，扎龍國家級自然保護區始終是景觀自然性值最高的地區，并且該保護區的景觀自然性呈上升狀態，1985—2015年的上升率為0.46%。莫莫格國家級保護區、查干湖國家級自然保護區的景觀自然性呈先上升、后下降狀態(表2)。莫莫格國家級自然保護區的景觀自然性在1985—2000年的上升率為9.86%，在2000—2015年的下降率為27.84%。查干湖國家級自然保護區的景觀自然性在1985—2000年的上升率為4.17%，在2000—2015年的下降率為16.21%。向海國家級自然保護區的景觀自然性值下降最快，下降率為34.32%。大布蘇狼牙壩國家級自然保護區的景觀自然性在1985—2015年的下降率為14.73%。

為探究松嫩平原西部景觀自然性的空間分異規律，利用ArcGIS 10.2軟件，將松嫩平原西部三年各個網格的景觀自然性值進行克里金插值，得到松嫩平原西部景觀自然性空間分布圖(圖3)。

表2 松嫩平原西部主要城市及自然保護區景觀自然性值

圖3 松嫩平原西部景觀自然性空間分布

由圖3可知，松嫩平原西部景觀自然性值較高的區域集中在研究區南部，以扎龍國家級自然保護區為中心向四周遞減。

1985年景觀自然性值最高的地區為扎龍國家級自然保護區及其北部，為0.197 5。這一區域自然景觀比例較高，為60.35%，植被覆蓋率較高，NDVI值為0.40，濕地是這一區域的優勢景觀，濕地率為28.29%。其次為向海國家級自然保護區及其南部，為0.188 8，這一區域人口密度為15.6人/km2，人口稀少，經濟發展水平較低，NDVI值為0.29，自然景觀比例為66.12%。1985年景觀自然性最低值出現在松原市扶余縣，為0.117 0，這一區域的優勢景觀為耕作景觀，耕地率為76.00%。

2000年景觀自然性值最高的地區為莫莫格國家級自然保護區中部及其北部，為0.211 6。這一區域自然景觀比例為58.27%，其中濕地率為20.80%、草地率為13.25%。2000年松嫩平原西部整體人為干擾值為0.591 7，而該區域人為干擾值為0.439 4，說明該區域受人為干擾程度較低。其次為扎龍國家級自然保護區中部及其北部，為0.171 6，這一區域自然景觀比例為65.79%，連通性為2.55，但該區域人口密度和GDP都高于松嫩平原西部整體平均值，導致自然性下降。2000年景觀自然性最低值出現在松原市長嶺縣附近，為0.109 4。這一區域自然景觀比例為11.64%，NDVI指數較小，為0.101 3，耕地率為76.67%。

2015年景觀自然性值最高的地區為扎龍國家級自然保護區東部，為0.222 5。該區域自然景觀比例為76.60%，其中濕地率為44.73%，草地率為28.69%，NDVI值為0.336 5，高于平均值0.263 7。2015年景觀自然性最低值出現在白城市、齊齊哈爾西部及白城南部，分別為0.090 0，0.100 0，0.103 0，這些區域均為人類活動強度較高的地區，優勢景觀為城市景觀、農業景觀，耕地率分別為79.28%，84.71%，77.00%。

2.2 松嫩平原西部景觀自然性的動態變化

利用ArcGIS 10.2軟件，將1985—2000年、2000—2015年、1985—2015年3個時段的景觀自然性變化值進行克里金插值，得到松嫩平原西部景觀自然性動態變化圖(圖4)。

圖4 1985-2015年松嫩平原西部景觀自然性動態變化

研究區大部分區域景觀自然性呈下降狀態，僅少數區域景觀自然性呈上升狀態，這些區域集中在自然保護區及保護區附近地區。1985—2000年，景觀自然性上升的地區主要集中在研究區中部的莫莫格國家級自然保護區，而景觀自然性下降的地區主要有向海國家級自然保護區和齊齊哈爾中部地區。2000—2015年，景觀自然性上升的地區主要為扎龍國家級自然保護區、莫莫格國家級自然保護區和松原部分地區，景觀自然性下降的地區主要為齊齊哈爾南部地區。

1985—2015年，景觀自然性上升的地區僅有扎龍國家級保護區及保護區以東、莫莫格國家級保護區、大慶南部地區、松原少數地區，其余地區均呈下降趨勢,尤以齊齊哈爾中部及南部、白城西部及南部下降最快。

2.3 松嫩平原西部景觀自然性的空間聚集性分析

利用Geoda軟件，距離閾值設為26 km，對1985年、2000年、2015年的景觀自然性值進行全局Moran′sI和局部空間自相關分析。3個時期的全局Moran′sI指數分別為0.036 9，0.052 1，0.011 4，呈先上升后下降趨勢。從空間分布上看，松嫩平原西部空間自相關結果的南北差距較大(圖5)。高—高自相關類型主要集中在研究區北部，特別是扎龍國家級保護區及其周邊區域，表明該保護區自然性值與周圍區域自然性值都較高，空間差異性小，呈空間正相關。低—低自相關類型主要集中在研究區南部，表明研究區南部自然性值較低，空間差異性小，呈空間正相關。

圖5 松嫩平原西部景觀自然性空間分布

從數量上看，高—高自相關類型、低—低自相關類型均呈先增加，后減少的變化趨勢(表3)。結合3個時期的克里金插值圖可知，1985年景觀自然性最高值出現在扎龍國家級自然保護區，隨著莫莫格國家級自然保護區自然性的上升，2000年景觀自然性最高值出現在莫莫格國家級自然保護區和扎龍國家級自然保護區兩個區域，故這一時期高—高自相關類型的數量達到最高值，2015年景觀自然性最高值僅出現在扎龍國家級自然保護區，故高—高自相關類型的數量相應降低。1985—2000年期間，研究區北部大部分地區景觀自然性均呈下降狀態，低—低自相關類型的數量隨之上升，而2000—2015年期間，隨著松原市長嶺縣等地景觀自然性值的上升，低—低自關類型數量也隨之減少。低—高自相關類型呈逐年增加趨勢，高—低自相關類型呈逐年減少趨勢。無空間自相關類型呈先減少、后增加趨勢。

表3 1985-2015年研究區局部空間自相關類型統計

3 結 論

(1) 景觀自然性受景觀破碎度、景觀優勢度、自然景觀比例、GDP這4個因子的影響較大，其中自然景觀比例所占的權重值最高，增加自然景觀的占比是提高景觀自然性的關鍵。

(2) 松嫩平原西部景觀自然性呈逐年下降趨勢，其中1985—2000年下降速率最快，2000—2015年下降率減慢，說明21世紀后推出的鹽堿地生態改良、禁牧休牧等措施產生了一系列成效。景觀自然性下降的原因主要是壓力指標的上升和響應指標的下降，日后的工作重點應以降低人類干擾為主。

(3) 景觀自然性值較高且呈上升狀態的地區主要為自然保護區及保護區周邊地區，如扎龍國家級保護區及保護區以東，莫莫格國家級保護區等。而人類活動密集的城市景觀自然性較低。造成這種差異的原因主要是自然保護區景觀的響應指標較高，而城市景觀的壓力指標較高。

(4) 松嫩平原西部景觀自然性的空間自相關結果南北差異較大，高—高自相關類型、低—低自相關類型均呈先增加，后減少的變化趨勢。

本文從景觀指數、土地利用類型、人類干擾三方面探討松嫩平原西部景觀自然性變化。影響一個地區景觀自然性的因素還有很多，如生物因素、致災因子等，這些因素的數據目前難以獲取而不能準確定量化，如何將這些因素定量化、探究這些因素對景觀自然性的影響可作為未來研究的重點。同時，本文僅對松嫩平原西部景觀自然性及其動態變化進行評價，未對評價結果進行驗證，如何科學、有效地驗證評價結果也是未來研究的重點之一。