基于移動窗口法的豫西山地丘陵地區景觀異質性分析

李棟科，丁圣彥,*，梁國付，趙清賀，湯 茜，孔令華

(1.教育部黃河中下游數字地理技術重點實驗室，開封 475004； 2.河南大學環境與規劃學院， 開封 475001)

基于移動窗口法的豫西山地丘陵地區景觀異質性分析

李棟科1,2，丁圣彥1,2,*，梁國付1,2，趙清賀1,2，湯 茜1,2，孔令華2

(1.教育部黃河中下游數字地理技術重點實驗室，開封 475004； 2.河南大學環境與規劃學院， 開封 475001)

研究典型區域景觀異質性過程中，對特征尺度的判斷尤為重要。基于3S技術，以豫西山地丘陵地區鞏義市為研究區域，在分析其景觀組分與結構的基礎上，采用半變異函數和移動窗口相結合的方法，確定研究景觀異質性的特征尺度，并在特征尺度下選取西北—東南方向和東北—西南方向兩條樣線，分析研究區不同干擾背景下景觀異質性的空間特征。結果表明：豫西鞏義地區總體上以旱地、建設用地、水澆地、草地和有林地景觀類型為主，其中旱地所占比例最大(30.9%)，溝渠所占比例最小(0.1%)；根據特征尺度分析確定研究區鞏義市景觀異質性的合適尺度為1000 m；在特征尺度下，研究區景觀破碎化以河谷平原和丘陵為主的鞏義東部偏南區域最高，以人為干擾較少的南部低山有林地區最低，景觀多樣性分布特征則呈現4個聚集中心，分別為北部河谷平原農業區、城市建成區、南部低山區和東南丘陵區；不同干擾背景下，鞏義地區景觀異質性根據城市化程度呈現以市區為中心向外輻射發展的特征，同時受海拔高度因素影響較大，即在海拔200 m以下區域景觀受人為干擾強烈，海拔200—500 m丘陵區域受海拔因子和人為干擾因子共同作用，海拔500 m以上區域海拔因子起主導作用。研究結果可為豫西山地丘陵地區及同類地區景觀異質性的研究尺度、景觀格局優化和土地的有效管理提供依據。

移動窗口法；特征尺度；景觀異質性；山地丘陵地區；河南西部

景觀生態學注重異質性的研究[1]，其研究核心是強調空間異質性與生態過程、尺度之間的相互作用關系[2- 5]。景觀異質性是空間斑塊性和空間梯度的綜合反映[6- 7]，是景觀格局的重要特征。分析景觀異質性是研究景觀功能和動態的基礎，對探索景觀結構與自然生態過程和社會經濟活動的關系[8- 9]、研究區土地資源合理利用和可持續發展均具有重要的理論和現實意義[10]。目前分析景觀異質性多是在深入理解各個景觀指數對應公式及生態學意義基礎上[6, 11- 12]，本著代表性、統一性和便于解釋原則，直接用景觀指數來表征景觀異質性。但是，景觀指數具有尺度效應、方向性及其他行為學特征[2, 13- 15]，目前的研究較少綜合以上特征對區域景觀異質性進行分析。

對于特定研究區域，選擇合適的研究尺度對研究結果尤為重要。有研究表明，利用地統計學中的半變異函數和移動窗口法都可以較好地分析景觀指數與空間變量的變異關系，以及判定研究區域的特征尺度。其中，移動窗口法最早由Whittaker提出用于分析植被沿水分梯度的變化[16]，隨后被應用于城鄉交錯帶景觀異質性等方面的研究[17]，其原理是在研究區內選取一定大小的窗口進行移動，形成新的可在地理信息系統(GIS)中運算的數據圖形[18]，實現在區域或地區尺度上對景觀指標的量化[19- 20]，從空間尺度上明確展示景觀異質性特征；半方差分析是把統計相關系數的大小作為一個距離的函數來反映一個采樣點與其相鄰采樣點的空間關系的一種地統計學方法，可用于分析特征尺度，目前已被廣泛應用于多個領域。但是，綜合半變異函數和移動窗口法來判定區域景觀特征尺度的研究較少。

鞏義市位于中岳嵩山北麓，地貌類型多樣，地勢南高北低，山地、丘陵、平原類型齊全。近年來，隨著區域社會經濟發展，該地區的景觀發生了巨大變化，引起較多關注。但是，大部分學者多從村域和農戶角度研究農田生態經濟系統投入產出特征、農田系統能量投入產出特征、農戶遷居模式等[21- 23]，以及從恢復生態學角度研究森林景觀恢復[24]，較少從景觀生態學角度研究不同干擾背景下景觀異質性變異特征。

因此，本研究基于地理信息系統，采用移動窗口法和半變異函數相結合的方法，以鞏義市為研究區域，探討該區域景觀異質性變化的特征尺度和景觀空間變異特征，以期為研究區域土地的有效管理及可持續利用提供依據。

1 研究區概況

作為豫西地區經濟活力較高的三縣之一的鞏義市[25]，位于黃河中游地區，處于國家第2臺階向第3臺階過渡邊緣，地理位置介于北緯34°31′—34°52′，東經112°49′—113°17′，在河南省省會鄭州與古都洛陽之間，南依嵩山，北臨黃河，地勢東南高西北低(圖1)。鞏義市土地總面積為1041 km2，其中山地、丘陵面積共919 km2，河谷、平原面積117 km2，屬典型的山地丘陵地區。山地分布著保護較完整的大面積林地，植物資源豐富，以溫帶落葉闊葉林為主；低山丘陵地區以荒山灌草叢、經濟林和人工林為主；河谷平川地區以精耕細作的農業為主，或者是農林間作。鞏義地處暖溫帶半濕潤偏干旱區，屬于溫帶季風氣候，四季分明，年平均氣溫14.4 ℃,7月最熱，平均氣溫27.3 ℃，1月最冷，平均氣溫0.2 ℃，平均的降雨量為640.9 mm，無霜期為220 d。

圖1 鞏義市地理位置圖及海拔梯度Fig.1 Location and DEM of Gongyi City

2 研究方法

2.1 數據處理

2.1.1 數據選擇

本研究選用的數據有：2011年鞏義市土地利用現狀圖(當地政府部門提供，比例尺較大，數據涉密)，比例尺為1∶1萬；2012年鞏義實地調查得出的景觀類型組成、地貌特征和高程等相關數據。

2.1.2 數據處理平臺

本研究所采用的數據處理軟件主要為：Fragstats 3.3、ArcGIS 9.3以及SPSS 17.0、GS+ 7.0數據分析等軟件。

2.1.3 景觀分類

在對研究區景觀功能認識的基礎上，根據土地利用現狀，將研究區景觀類型劃分為旱地、建設用地、水澆地、草地、有林地、灌木林、河流、內陸灘涂、道路、果園、水庫坑塘、裸地、溝渠共13種景觀類型，并由矢量格式轉換為柵格格式，柵格大小為30 m×30 m，進而進行格局指數的計算。

2.2 景觀指數的選取

景觀破碎化程度和多樣性指數是衡量景觀異質性的重要指標。根據前人研究，本文采用最大斑塊指數(LPI)、斑塊密度(PD)代表景觀破碎化程度，破碎化程度越高，空間異質性程度就越大，反之破碎化程度越低，空間異質性就越小；采用香農多樣性指數(SHDI)、聚集度指數(CONTAG)、景觀均勻度指數(SHEI)代表研究區的景觀多樣性，其中，SHDI代表景觀要素的多少和各景觀要素面積占比的變化，CONTAG描述景觀中不同生態系統的團聚程度或延展趨勢，SHEI則表示景觀中不同生態系統的分配均勻程度。

2.2.1 最大斑塊指數(LPI)

LPI用于揭示最大斑塊對整個景觀(或類型)的影響程度:

(1)

式中，aij是景觀中斑塊的面積，A為景觀總面積。

2.2.2 斑塊密度指數(PD)

PD反映斑塊破碎化程度，同時也反映景觀空間異質性程度。斑塊密度越大，破碎化程度就越高，空間異質性程度就越大，反之，破碎化程度越低，空間異質性就越小:

(2)

式中，PD整個景觀的斑塊密度，Ni斑塊數量，A景觀的總面積。

2.2.3 香農多樣性指數(SHDI)

SHDI值大小反映景觀要素的多少和各景觀要素面積占比的變化:

(3)

式中，Pi表示景觀類型i所占面積的比例，m表示景觀類型數。

2.2.4 聚集度指數(CONTAG)

CONTAG用于描述景觀中不同生態系統的團聚程度或延展趨勢：

(4)

式中，Pi為第i種景觀類型所占景觀總面積的比例，m為景觀類型種類數。gik為隨機選擇的兩個相鄰的柵格屬于類型i和k的概率。蔓延度指數通常是用來度量景觀斑塊形狀的復雜程度。取值范圍：0

2.2.5 景觀均勻度指數(SHEI)

SHEI用于揭示景觀中不同生態系統的分配均勻程度。

(5)

式中，Pi為第i種景觀類型占景觀總面積的比例，m為景觀類型的種類數。

2.3 移動窗口法

首先，將矢量格式的土地利用圖轉換為柵格格式(30 m×30 m)；其次，將柵格數據導入Fragstats3.3軟件中，采用標準法計算研究區整體景觀水平上最大斑塊指數(LPI)、斑塊密度(PD)、香農多樣性指數(SHDI)、聚集度指數(CONTAG)、景觀均勻度指數(SHEI)等景觀指數；利用移動窗口法，移動窗口半徑分別設為500，750，1000，1250，1500，1750，2000 m，從研究區左上角開始移動，每次移動1個柵格，計算每個窗口內上述景觀指數值(根據Fragstats3.3使用說明書，對于窗口邊界上不能被整分的柵格由Fragstats3.3軟件自動判定取舍)，并賦予各個窗口的中心柵格，最終獲得各個景觀指數的空間分布柵格圖。

2.4 半變異函數

地統計學是一系列檢測、估算變量在空間上的相關關系和格局的統計方法[26]，半方差分析(又稱變異函數分析)是地統計學中的一個重要組成部分[27- 30]，目前已在生態學、土壤學等領域應用。本研究在移動窗口法的基礎上，采用地統計學中的半變異函數方法判定景觀特征尺度。半變異函數是把統計相關系數的大小作為一個距離的函數。區域化變量Z(x)在點x和x+h處的值Z(x)與Z(x+h)差的方差的一半稱為區域化變量Z(x)的半變異函數[31]，記為r(h)。定義為：

r(x,h)= 1/2Var[Z(x)-Z(x+h)]

半變異函數曲線可以反映一個采樣點與其相鄰采樣點的空間關系。在半變異函數曲線上有兩個非常重要的點，即間隔為0時的點和半變異函數趨勢平穩時的拐點，由這兩個點產生4個相應的參數：塊金值(Nugget)、變程(Range)、基臺值(Still)、偏基臺值(Partial Still)。理論上，空間變異是永恒存在的，樣點任意一方或兩者共同作用產生了塊金值(Nugget)，記為C0，當采樣點間距增大時，r(h)從初始的塊金值達到一個相對穩定的常數，該常數就為基臺值，記為C+C0。基臺值與塊金值間的差值為偏基臺值，也就是結構方差C。塊金值與基臺值比值即為塊基比，記為C0/(C+C0)，其值可以估算隨機因素在空間總變量中重要性，進而反映變量在空間上的變異程度。C0/(C+C0)的比值大小代表空間變異程度的高低，其值越小，說明空間自相關越明顯，也就越穩定[23,32- 33]。本研究通過地統計軟件GS + 模擬各個移動窗口下景觀指數半變異函數模型，根據各個景觀指數的C0/(C+C0)值的變化規律，分析不同指數的空間變異特征對尺度變換的響應。結合前人研究，當比值達到穩定狀態時，即說明景觀指數在空間的變異趨于穩定，可判定為研究區域景觀的特征尺度[32- 33]。

2.5 樣線設置

鞏義市地勢是自東南向西北呈階梯狀急劇降低，分別為中山、低山、丘陵、河谷平川,樣線選取盡量覆蓋所有景觀類型。本文選取西北—東南方向樣線，長為28.3 km；東北—西南方向樣線，長為30.5 km；然后基于特征尺度，利用ArcGIS軟件中3D Analyst功能模塊將兩條樣線上的信息提取并進行梯度變化分析。

圖2 鞏義市2011年景觀分類圖及樣線設置Fig.2 The classified map of Gongyi city and transects setting in 2011

3 結果與分析

3.1 研究區景觀的基本特征

該地區占較大面積的景觀類型是旱地，面積為322.0 km2，占總面積的30.9%，為該地區的基質；其次為建設用地，水澆地，草地和有林地的面積，其面積分別為190.0，142.4、135.2 km2和125.8 km2，占總面積的比例分別為18.2%、13.7%、13.0%、12.1%；其他景觀類型面積比較小，均未超過總面積的10.0%(表1)。

3.2 特征尺度

根據景觀格局指數空間變異特征值趨勢圖(圖3)，在多個連續尺度上對區域景觀斑塊組成與構型的空間變異特征進行探討，且在此基礎上確定景觀斑塊組成與構型的分析尺度(特征尺度)。

圖3中橫坐標(X軸)代表不同的窗口半徑，縱坐標(Y軸)代表的是塊基比值(C0/(C0+C))。塊基比值越低，空間變異程度越低，空間自相關性越明顯，也就越穩定；越高，空間變異程度越高，越不穩定。3個景觀格局指數的穩定性隨著窗口半徑的增加呈現先增加后降低，到一定數值后趨于穩定，因此可以進行特征尺度的判定。結果表明，粒度750 m和粒度為1000 m時，趨勢圖上出現明顯的拐點。窗口為500 m和750 m時呈現上升趨勢，變化不穩定，因此750 m并不能反映研究區域的特征尺度。3個景觀指數的空間變異特征值在1000 m左右時開始趨于穩定，表明該尺度可以反映研究區景觀格局變化空間變異特征的內在尺度。另外，尺度過大(1250、1500、1750、2000 m)會致使較多的空間信息規律損失。因此，本研究選擇1000 m作為分析鞏義市景觀格局的理想特征分析尺度。

表1 各景觀類型的面積及所占比例

圖3 景觀格局指數空間變異特征值趨勢圖Fig.3 The trend of characteristic values of spatial heterogeneity within different extents

3.3 特征尺度下景觀異質性特征

3.3.1 景觀破碎化空間分析

2011年研究區景觀破碎化空間分布狀況如圖4所示，最大斑塊指數有3個高值區，分別為鞏義市區、西部河谷平原城鎮建成地區和南部低山有林地區，其中最大高值斑塊為鞏義市區。鞏義市區和西部河谷平原城鎮建成區林地和綠地面積較少，景觀類型以建設用地為主，人為干擾劇烈；而南部低山有林地區，由于封山育林，人工油松林已初成規模，林地面積較大，景觀類型以有林地為主，且海拔在800 m以上，地勢較高，人為干擾因素較弱。斑塊密度在以上3個地區較小，是由于該地區優勢類型斑塊較大且較完整，破碎化程度較低，空間異質性較小。斑塊密度較高地區位于鞏義東部偏南區域，大致由鞏義城區東南部近郊向鞏義東部延展，地形以河谷平原和丘陵為主，地勢較平坦，景觀類型多樣，受人類活動影響較大，所以該地區破碎化程度較高，并且圖中顯示兩指數的相關性較高。

圖4 研究區景觀破碎化空間分布Fig.4 Spatial distribution of landscape fragmentation of the study area

3.3.2 景觀多樣性空間分析

景觀多樣性空間分布如圖5所示，香農多樣性指數低值區出現在東南丘陵地區，南部低山地區和中部建成區，尤其在有林地占絕對優勢東南部地區，景觀類型較為單一，景觀幾乎是均質的，香農多樣性值較低；香農多樣性指數高值出現在市區北部，該地區為平原地區，開發較早，擁有水澆地、旱地、林地、水域、建設用地等景觀類型。景觀均勻度指數表現出與其相似分布狀況。聚集度指數與景觀多樣性指數有較高負相關性，鞏義市區和中部建成區，以建設用地、工礦用地為主，受人類活動影響較大，大斑塊團聚程度較高，呈規則形狀分布；東南丘陵地區和南部低山丘陵地區以荒山灌叢、有林地為主，不同生態系統團聚程度低。

圖5 研究區景觀多樣性空間分布Fig.5 Spatial distribution of landscape diversity of the study area

3.4 不同干擾背景下景觀異質性特征分析

3.4.1 人為干擾背景下景觀異質性分析

在西北—東南方向上，鞏義市景觀異質性特征如圖6所示。從遠郊到近郊，最大斑塊指數總體呈下降趨勢，在離中心城區約5 km處達到拐點，隨后急劇上升，在市中心達到頂點，隨后在城市中心東南近郊最大斑塊指數迅速下降。這是由于城市建設用地景觀類型占絕對優勢，人為干擾劇烈所致。最大斑塊指數總體呈現“W”形狀；聚集度指數呈現類似結果。景觀均勻度指數從遠郊到近郊不斷增加，在近郊離城市中心約5 km處達到最大，然后不斷減少，在市中心達到最小值，隨著離城市中心越來越遠，指數逐漸增加，在遠郊又呈下降趨勢，并向東南方向不斷延伸，總體呈現“M”形狀。香農多樣性和斑塊密度也呈現類似結果。

圖6 景觀水平上的景觀指數沿西北—東南方向樣線變化Fig.6 Change in landscape-level metrics along the NW-SE transect

在東北—西南方向上，鞏義市景觀異質性特征如圖7所示。遠郊到市中心距離約13 km，最大斑塊指數和聚集度指數的峰值出現在市中心，在近郊變化較為劇烈，出現多次的波動，在市中心西南又有回升趨勢，均勻度指數，多樣性指數和斑塊密度出現類似的特點，在市中心出現谷值，在近郊出現峰值，在遠郊指數相對較低。東北—西南方向的景觀格局指數在城市郊區波動具有相類似的特點，變化明顯的地區出現在離市中心7 km范圍，有多次頻繁的波動，景觀類型復雜多樣。

圖7 景觀水平上的景觀指數沿東北—西南方向樣線變化Fig.7 Change in landscape-level metrics along the NE-SW transect

通過對比兩條樣線可以發現，景觀破碎化程度以市區為中心向外輻射波動。在城市中心斑塊密度達到最低值，是由于斑塊類型較為單一，主要以建設用地為主，數量較少，但每個斑塊的面積較大，從而造成斑塊密度較小。而在城市郊區，斑塊類型多樣，但每個斑塊的面積較小，從而造成斑塊密度較大，甚至出現峰值。所以由城市中心到郊區景觀由單一、均質和連續的整體趨向于復雜、異質和不連續的斑塊鑲嵌體。但是，由于選取的兩條樣線方向不同，景觀指數的變化趨勢也有所不同。東北-西南方向景觀格局指數變化比較劇烈，這可能與研究區地勢有關，如鞏義市中心向西南遠郊方向地勢較平坦，城市建設用地和水澆地占較大優勢，受人為活動影響較大，而鞏義市中心向東南遠郊方向以低山為主，有林地斑塊占較大優勢，建設用地和農業用地相對較少，受人為活動影響較小。因此，推斷鞏義市景觀異質性變化受海拔因素影響較大。

3.4.2 景觀異質性隨海拔高度變化規律

豫西鞏義市地處在嵩山北麓，最高峰玉柱峰海拔高度為1487 m，最低點鞏義市沙魚溝鄉洛口灘，海拔為103.8 m，上下海拔落差為1383.2 m，其地勢是自南部向西北呈階梯狀急劇下降。南部和東南部分地區海拔在1000 m以上，屬中山區域；由中山區域向西北方向是一系列海拔在500—1000 m的低山，構型成低山區；由低山區域再向西北方向為廣闊的黃土丘陵區，其海拔為200—500 m；丘陵區西北向即為其境內河谷平原區。

為了進一步分析景觀異質性隨海拔高度的變化規律，結合圖2，圖6和表2發現，在研究區西北-東南方向上，各景觀指數曲線圖的極值均出現在低海拔受人為干擾嚴重的城市建成區中心附近。對比最大斑塊指數(LPI)曲線、香農多樣性指數(SHDI)曲線和景觀均勻度指數(SHEI)曲線發現，自西北向東南方向過城市中心后，海拔高度升至200 m出現了較為明顯的較小波幅的連續波動，說明隨海拔高度攀升，自海拔200 m后人為因素減少而海拔因子開始顯著作用于該研究區景觀，同時也表明海拔200—500 m黃土丘陵區是海拔因子和人為干擾共同作用的區域。當海拔進入500—1000 m低山區向高海拔區延展的過程中，景觀異質性呈均衡平穩波動的態勢，證明隨著海拔梯度的快速躍升，景觀受海拔因子的影響強烈，人類干擾因素顯著降低。以灌木林和有林地為主的低山區景觀由于海拔的變化被不同程度的自然分割，適度破碎化，斑塊數量略有增多，進而增加了生境類型，更有利于生物多樣性的保護。另外最大斑塊指數(LPI)曲線上峰弧度平滑，略有跨度，下峰急變斗轉無弧度，表明當地海拔落差變化較快，確實對景觀產生了明顯影響。對比斑塊密度(PD)曲線和聚集度指數(CONTAG)曲線發現，進一步印證當海拔從500—1000 m低山區到高海拔區延展的過程中出現較為明顯的均衡平穩波動并維持在一個適度水平上，表明隨著海拔的攀升，各類型斑塊面積占比趨于穩定，景觀中不同生態系統的分配均勻度適中，人為干擾因素明顯降低，海拔因子對500 m以上高海拔區域的景觀變化起主導作用。

綜上所述，豫西鞏義地區海拔200 m以下區域景觀受人為干擾強烈；海拔200—500 m丘陵區域是海拔因子和人為干擾因子共同作用的過度區域，且隨著海拔上升人為干擾因子對景觀作用力逐步減弱；海拔因子對500 m以上高海拔區域的景觀變化起主導作用。

表2 景觀異質性與海拔的相關性

4 結論和討論

本文基于鞏義市2011年30 m×30 m柵格圖像數據，運用景觀生態學、空間統計學的理論與方法[34- 35]，在地理信息系統技術支持下，對豫西山地丘陵典型地區的鞏義市的景觀格局進行探討。研究表明，基于GIS的梯度和景觀指數分析，利用移動窗口的方法定量地測定景觀空間格局及其變化規律是可行的，且可實現景觀指標的空間化與可視化，有利于直觀地了解區域景觀格局和對干擾進行定量化研究[36- 37]。

采用移動窗口計算景觀格局指數時，移動窗口半徑大小的選擇對于分析景觀格局的變化有重要影響。在鞏義市山地丘陵地區，當移動窗口半徑較小時，景觀格局指數變化較大，當移動窗口半徑到1000 m或更大時，景觀格局指數變化相對比較平緩。本研究結果表明，采用1000 m移動窗口半徑來分析研究區景觀格局時，既能保留梯度特征，又能通過景觀指標的變化真實地反映空間格局的變化。在不同的區域，隨著景觀組成和構型的變化，景觀異質性也在不斷變化，分析景觀格局時采用的移動窗口半徑也不一樣，但是本研究提出的方法，可以為研究其它地區時選擇合適的研究尺度提供參考。

移動窗口形成的景觀指數空間分布圖，集中形成了4個中心，分別為鞏義市區、東南丘陵地區、南部低山地區和北部平原河川農業區，這與鞏義地區空間格局的分布特點密切相關。中心地區周圍約5 km處是景觀格局指數變化最強烈的地區，這4個中心恰恰代表鞏義4種不同程度的人為干擾區。市區是人為干擾最大的區域，山地區人為干擾最小，景觀表現為均質化的城市景觀，景觀要素類型少，最大斑塊指數大，多樣性低；河川平原和丘陵地區人為干擾程度居中，多樣性高，聚集度低。人類活動對鞏義景觀的影響取決于離城市中心的距離，沿兩樣線能清楚識別出城市化中心，鞏義城市空間的擴展受其特殊地形結構的影響，主要從市區沿東部和北部進行拓進和擴張。從另一個方面說明，所選用表示景觀組成和構型方面的指數，如多樣性指數、蔓延度指數等可以很好地表征研究區景觀格局特征，對于這些指數深入地分析，可以很好理解和區域景觀異質性特征。此外，除人為干擾外，海拔高度對研究區景觀異質性也起到相當重要的作用，如在地勢較為平坦的低海拔地區(200 m以下)受人為影響顯著，在中海拔丘陵地區(200—500 m)景觀異質性受海拔高度和人為干擾共同作用，而在500 m以上的高海拔區域，海拔因子對景觀變化起主導作用。

景觀格局指數是有尺度相關性的，景觀格局指數本身的尺度相關性是怎樣的？在表征景觀異質性時，隨尺度的變化又是怎樣的？這些問題在理解景觀異質性特征、特別是借用景觀格局指數來分析時，是急需要回答的問題，在下一步研究中還需要深入考慮和分析。此外，分析景觀異質性要特別要從景觀不同的粒度和幅度來研究，研究中主要探討了在30 m粒度下，不同干擾程度地區的景觀異質性的幅度相關特征，而沒有分析別的粒度，以及不同粒度和幅度組合下景觀異質性的特征，這對全面理解景觀異質性特征有很重要意義，這也是在以后研究中需要深入分析的地方。

[1] Risser P G.Champaign: Illinois Natural History Survey // Landscape Ecology: Directions and Approaches, Special Publication Number 2.Champaign, IL: Illinois Natural History Survey, 1984.

[2] Forman R T T, Godron M.Landscape Ecology.New York: John Wiley and Sons, 1986: 1- 47.

[3] Naveh Z.Some remarks on recent developments in landscape ecology as a transdisciplinary ecological and geographical science.Landscape Ecology, 1991, 5(2): 65- 73.

[4] Wu J G.Landscape Ecology, Cross-disciplinarity, and sustainability science.Landscape Ecology, 2006, 21(1): 1- 4.

[5] Turner M G.Landscape ecology: the effect of pattern on process.Annual Review of Ecology and Systematics, 1989, 20: 171- 197.

[6] Wu J G.Landscape Ecology—Pattern, Process, Scale and Hierarchy.Beijing: Higher Education Press, 2007: 19- 20.

[7] Pickett S T A, Cadenasso M L.Landscape ecology: spatial heterogeneity in ecological systems.Science, 1995, 269(5222): 331- 334.

[8] Turner M G.Landscape ecology in North America: past, present, and future.Ecology, 2005, 86(8): 1967- 1974.

[9] Forman R T T.Some general principles of landscape and regional ecology.Landscape Ecology, 1995, 10(3): 133- 142.

[10] Gao J B, Cai Y L.Spatial heterogeneity of landscape fragmentation at multi-scales.Scientia Geographica Sinica, 2010, 30(5): 742- 748.

[11] Ye Y Q, Chen G J.An analysis on the landscape patterns based on the GIS technology in the upper Min River basin.Resources and Environment in the Yangtze Basin, 2006, 15(1): 112- 115.

[12] Sui X, Tong F, Yao C Q, Yang Z F.Landuses and landscape patterns of Qinghai Province in the Yellow River Basin.Bulletin of Soil and Water Conservation, 2007, 27(2): 123- 127.

[13] Wu J G.Landscape Ecology―Concepts and Theorieshinese.Journal of Ecology, 2000, 19(1): 42- 52.

[14] Bu R C, Hu M Y, Chang Y, Li X Z, He H S.A correlation analysis on landscape metrics.Acta Ecologica Sinica, 2005, 25(10): 2765- 2775.

[15] Shen W J, Wu J G, Lin Y B, Ren H, Li Q F.Effects of changing grain size on landscape pattern analysis.Acta Ecologica Sinica, 2003, 23(12): 2506- 2519.

[16] Whittaker R H.Vegetation of the Siskiyou Mountains, Oregon and California.New York: Ecological Society of America, 1960: 30- 49.

[17] McDonnell M J, Piekett S T A.Ecosystem structure and function along urban-rural gradients: An unexpected opportunity for ecology.Ecology, 1990, 71(4): 1231- 123.

[18] Baker W L, Cai Y M.The r.le programs for multiscale analysis of landscape structure using the GRASS geographical information system.Landscape Ecology, 1992, 7(4): 291- 302.

[19] Yin H W, Kong F H.Spatio-temporal gradient analysis is of urban green space in Ji′nan City.Acta Ecologica Sinica, 2005, 25(11): 3010- 3018.

[20] Pickett S T A, Cadenasso M L, Grove J M, Nilon, C H, Pouyat R V, Zipperer W C, Costanza R.Urban ecological systems: Linking terrestrial ecological, physical, and socioeconomic components of metropolitan areas.Annual Review of Ecology and Systematics, 2001, 32: 127- 157.

[21] Li X J, Qiao J J.Impact of landform on input-output of man-land system in farmland of mountainous region: a micro-study of a small village in Wugou of Henan Province.Geographical Research, 2004, 23(6): 717- 727.

[22] Qiao J J, Ding P F.Comparative study of energy input and output in villages farmland systems: the case of three types of villages in Gongyi county, Henan Province.Resources Science, 2004, 26(5): 139- 146.

[23] Qiao J J, Li X J.Impact of LULC on human-land systems: a micro-study of three villages in Gongyi City, Henan Province.Research of Soil and Water Conservation, 2006, 13(2): 91- 94.

[24] Chen J, Liang G F, Ding S Y.Landscape connectivity analysis for the forest landscape restoration: a case study of Gongyi City.Beijing: Acta Ecologica Sinica, 2012, 32(12): 3773- 3781.

[25] Qiao J J.Rural Territorial Economy of China.Beijing: Science Press, 2008: 128- 129.

[26] Journal A, Huijbregis C.Mining Geostatistics.New York: Academic Press, 1978.

[27] Matheron G.Principles of geostatistics.Economic Geology, 1963, 58(8): 1246- 1266.

[28] Griffith J A, Martinko E A, Price K P.Landscape structure analysis of Kansas at three scales.Landscape Urban Plan, 2000, 52(1): 45- 61.

[29] Fortin M J, Boots B, Csillag F, Remmel T K.On the role of spatial stochastic models in understanding landscape indices in ecology.Oikos, 2003, 102(1): 203- 212.

[30] Smith A C, Koper N, Francis C M, Fahrig L.Confronting collinearity: comparing methods for disentangling the effects of habitat loss and fragmentation.Landscape Ecology, 2009, 24(10): 1271- 1285.

[31] Cai B F, Yu R.Comparison on spatial scale analysis methods in landscape ecology.Acta Ecologica Sinica, 2008, 28(5): 2280- 2282.

[32] Xue D D, She G H, Wen X R, Sun H H, Li Q W, Zhou Y C.Scale effect of landscape pattern of Nanjing Zhongshan Scenic Spot based on statistic analysis.Journal of Southwest Forestry University, 2012, 32(1): 30- 35.

[33] Yue W Z, Xu J H, Xu L H, Tan W Q, Mei A X.Spatial variance characters of urban synthesis pattern indices at different scales.Chinese Journal of Applied Ecology, 2005, 16(11): 2053- 2059.

[34] Xiao D N.Landscape Ecology.Beijing: Science Press, 2003: 89- 124.

[35] Yu Z R.Landscape Ecology.Beijing: Chemical Industry Press, 2008: 44- 78.

[36] Liu X, Guo Q X.Landscape pattern in Northeast China based on moving window method.Chinese Journal of Applied Ecology, 2009, 20(6): 1415- 1422.

[37] Wang J, Sun J, Yang X J, Liu W Z, Zaccarelli N, Zhang X L, Wang X Y.An analysis of disturbance on social-ecological system a tmultiple scales based on NDVI, case study in Yuzhong county of Gansu Province.Acta Ecologica Sinica, 2009, 29(3): 1622- 1628.

參考文獻:

[6] 鄔建國.景觀生態學——格局、過程、尺度與等級.北京: 高等教育出版社, 2007: 19- 20.

[10] 高江波, 蔡運龍.區域景觀破碎化的多尺度空間變異研究.地理學報, 2010, 30(5): 742- 748.

[11] 葉延瓊, 陳國階.GIS支持下的岷江上游流域景觀格局分析.長江流域資源與環境, 2006, 15(1): 112- 115.

[12] 隋欣, 童飛, 姚長青, 楊志峰.黃河流域青海片土地利用景觀空間格局分析.水土保持通報, 2007, 27(2): 123- 127.

[13] 鄔建國.景觀生態學―概念與理論.生態學雜志, 2000, 19(1): 42- 52.

[14] 布仁倉, 胡遠滿, 常禹, 李秀珍, 賀紅士.景觀指數之間的相關分析.生態學報, 2005, 25(10): 2765- 2775.

[15] 申衛軍, 鄔建國, 林永標, 任海, 李勤奮.空間粒度變化對景觀格局分析的影響.生態學報, 2003, 23(12): 2506- 2519.

[19] 尹海偉, 孔繁花.濟南市城市綠地時空梯度分析.生態學報, 2005, 25(11): 3010- 3018.

[21] 李小建, 喬家君.地形對山區農田人地系統投入產出影響的微觀分析——河南省鞏義市吳溝村的實證研究.地理研究, 2004, 23(6): 717- 727.

[22] 喬家君, 丁鵬飛.村域農田系統能量投入產出特征比較研究——以河南省鞏義市3個不同類型村為例.資源科學, 2004, 26(5): 139- 146.

[23] 喬家君, 李小建.土地利用/土地覆被對人地系統的影響分析——以河南省鞏義市3個不同類型村為例.水土保持研究, 2006, 13(2): 91- 94.

[24] 陳杰, 梁國付, 丁圣彥.基于景觀連接度的森林景觀恢復研究——以鞏義市為例.生態學報, 2012, 32(12): 3773- 3781.

[25] 喬家君.中國鄉村地域經濟論.北京: 科學出版社, 2008: 128- 129.

[31] 蔡博峰, 于嶸.景觀生態學中的尺度分析方法.生態學報, 2008, 28(5): 2280- 2282.

[32] 薛冬冬, 佘光輝, 溫小榮, 孫浩晗, 李勤文, 周易春.基于地統計分析的南京鐘山風景區景觀格局尺度效應分析.西南林業大學學報, 2012, 32(1): 30- 35.

[33] 岳文澤, 徐建華, 徐麗華, 談文琦, 梅安新.不同尺度下城市景觀綜合指數的空間變異特征研究.應用生態學報, 2005, 16(11): 2053- 2059.

[34] 肖篤寧.景觀生態學.北京: 科學出版社, 2003: 89- 124.

[35] 宇振榮.景觀生態學.北京: 化學工業出版社, 2008: 44- 78.

[36] 劉昕, 國慶喜.基于移動窗口法的中國東北地區景觀格局.應用生態學報, 2009, 20(6): 1415- 1422.

[37] 王俊, 孫晶, 楊新軍, 劉文兆, Zaccarelli N, 張向龍, 汪興玉.基于NDVI的社會- 生態系統多尺度干擾分析——以甘肅省榆中縣為例.生態學報, 2009, 29(3): 1622- 1628.

LandscapeheterogeneityofmountainousandhillyareainthewesternHenanProvincebasedonmovingwindowmethod

LI Dongke1,2，DING Shengyan1,2,*，LIANG Guofu1,2，ZHAO Qinghe1,2，TANG Qian1,2，KONG Linghua2

1KeyLaboratoryofGeospatialTechnologyFortheMiddleandLowerYellowRiverRegions,MinistryofEducation，Kaifeng475004，China2CollegeofEnvironment&Planning,HenanUniversity,Kaifeng475001,China

In analyzing landscape heterogeneity in a typical area, appropriate selection of characteristic scale is very important.In this paper, Gongyi City, located in mountainous and hilly area of western Henan Province, has been chosen as our study area.Based on 3S technology, the landscape composition and structure of this area are analyzed.Meanwhile, the characteristic scale of landscape heterogeneity is determined by using the semi-variable function and moving window method.According to the characteristic scale, we selected two sample lines, which extend from the direction of NW to SE and the direction of NE to SW.These two sampled lines are set to analyze the spatial characteristics of landscape heterogeneity under different disturbances.Finally, the results show that landscape of the study area mainly consists of dry land, construction land, irrigated land, grassland and woodland.Among these landscape types, the dry lands have the largest proportion, accounting for 30.9% of the whole area, while the irrigation canals have the least area proportion, accounting for 0.1% of the whole area.Based on the semi-variable function and the moving window method, the characteristic scale for analyzing landscape heterogeneity of the study area is 1000 m.Under this characteristic scale, the landscape fragmentation of the study area is mainly concentrated in the valley plain and hilly regions of the southeastern part of the Gongyi City.The lowest fragmentation regions are in the forestland of southern hilly area of Gongyi City.Landscape diversity of the study area is characterized by four gathering centers, which are the northern valley plain agricultural area, constructed urban areas, the southeast hills and the low mountains in southern area.The first center exhibits the highest landscape diversity, while the landscape diversity of the other three is of the lowest.Under the influence of different disturbances, the landscape heterogeneity of the study area show features of radiating from the city center to surrounding rural area.In addition, landscape heterogeneity is also related to altitude in hilly regions.We find that the landscape heterogeneity has dramatic changes below 200 m because of strong human disturbance.However, at the altitude of 200—500 m, human disturbance and altitude are combined to affect the landscape heterogeneity.In addition, at 500 meters above sea level, altitude is the dominant factor in influencing landscape heterogeneity.The results of this study can provide valuable guidance for further study on landscape pattern optimization and efficient management of land use in the mountainous and hilly area of western Henan Province.

moving window method; characteristic scale; landscape heterogeneity; mountainous and hilly area; Western Henan Province

國家自然科學基金資助項目(41371195,41071118)

2013- 10- 28;

2014- 04- 01

10.5846/stxb201310282595

*通訊作者Corresponding author.E-mail: syding@henu.edu.cn

李棟科，丁圣彥，梁國付，趙清賀，湯茜，孔令華.基于移動窗口法的豫西山地丘陵地區景觀異質性分析.生態學報,2014,34(12):3414- 3424.

Li D K，Ding S Y，Liang G F，Zhao Q H，Tang Q，Kong L H.Landscape heterogeneity of mountainous and hilly area in the western Henan Province based on moving window method.Acta Ecologica Sinica,2014,34(12):3414- 3424.