黃土丘陵區土地利用變化與地形梯度關系研究

莊逐舟, 黃秋昊, 石 云

(1.南京大學 地理與海洋科學學院， 南京 210023; 2.寧夏大學 資源環境學院， 銀川 750021)

?

黃土丘陵區土地利用變化與地形梯度關系研究

莊逐舟1, 黃秋昊1, 石 云2

(1.南京大學 地理與海洋科學學院， 南京 210023; 2.寧夏大學 資源環境學院， 銀川 750021)

土地利用/土地覆被變化(LUCC)是人類活動與自然環境相互作用最直接的表現形式，是反映人類活動程度的重要因子。以彭陽縣為例，應用1995年和2014年25 m×25 m的DEM數據和土地利用現狀數據，運用GIS空間分析功能，采用地形位指數以及分布指數定量分析了1995—2014年彭陽縣土地利用基于地形梯度的變化特征。結果表明：(1) 1995—2014年，彭陽縣耕地和未利用地面積大量減少，草地面積大量增加，林地、水域和建筑用地面積均有不同程度的增加；(2) 隨著地形位的升高，耕地、水域、建筑用地的分布指數逐漸降低，林地和草地的分布指數逐漸升高，未利用地的分布指數先升后降。草地的優勢分布區間增加了2個地形位，未利用地的優勢分布區間降低了2個地形位；(3) 水域的優勢度區間集中在低地形位區間，耕地和建筑用地的優勢度區間集中于低中地形位區間，林地、草地和未利用地的優勢度區間集中于中高地形位區間。

土地利用變化; 地形位指數; 分布指數; 彭陽縣

土地利用變化是土地利用/覆被變化研究的重要內容[1-2]，也是全球環境變化和可持續發展研究的熱點領域[3-8]。地形是自然環境中的重要因素，地形特征在一定程度上決定土地利用的方向和方式[9-13]，其空間特征也影響是土地利用的演變過程[14]。近年來，李增加等[15]分析了西雙版納1976—2007年不同地形下的土地利用變化；崔步禮[16]等以青海湖流域為例，從地形角度(坡度、坡向、高程)分析了流域內土地利用的空間分布格局；邱揚等[17]從土地利用類型與地形分布格局的變化分析著手，研究黃土丘陵溝壑區小流域土地利用的時空分布格局及其影響因子；陳楠等[18]基于地形位指數與分布指數，探究了巴中市不同地形梯度下不同土地利用結構的分布特征。上述研究均揭示了土地利用在地形梯度上呈現的變化規律，但研究區域較為集中于中國西南山區以及南方低山丘陵區，涉及黃土丘陵溝壑區土地利用空間分布與變化規律的研究相對較少。彭陽縣位于寧夏回族自治區南部，地處黃土丘陵溝壑區，也是全國的重點水土保持區，以彭陽縣為研究區，利用GIS的空間分析功能，結合25 m×25 m的DEM數據以及土地利用現狀數據，采用地形位指數以及分布指數，定量分析彭陽縣土地利用變化與地形梯度的關系，為黃土丘陵溝壑區縣域土地利空間格局研究、土地利用布局的優化及土地資源的可持續利用提供理論依據。

1　研究區概況

彭陽縣位于寧夏回族自治區東南部，地處黃土高原中部丘陵溝壑區，地貌由梁峁、殘垣、川地、和溝壑組成，其中以梁峁狀丘陵為主，屬黃土丘陵地貌，其間有河谷川臺及殘塬分布，地理坐標介于東經106°32′—106°58′，北緯35°41′—36°17′。土地總面積2 529 km2，海拔1 286～2 416 m，地形西北高，東南低，由西北向東南呈波狀傾斜。境內有紅河、茹河、安家川河縱橫東西(均流入涇河)，三河流域形成以6塊殘塬、2條長川、1條碎陵狹長河谷為主體，5條大溝、100多條小溝組成的溝壑縱橫的復雜地貌，是寧夏地區水土流失最為嚴重的縣域之一，同時也是全國重點水土保持區。

2　數據準備及研究方法

2.1數據準備

研究選用數據源包括彭陽縣25 m×25 m，1995年TM影像和2014年資源三號高分遙感影像，目視判讀解譯的1995年和2014年兩期土地利用數據，并結合黃土丘陵溝壑區自身特點，將研究區土地利用/覆被類型劃分為耕地、林地、草地、水域、建筑用地和未利用地六類。其他數據包括統計及社會經濟數據。

2.2研究方法

2.2.1數字地形分析結合彭陽縣的地形地貌特征，應用25×25 m的DEM數據以及地形分析，將由DEM中提取的高程、坡度信息進行分級(圖1—2)。高程以200 m為間隔，按<1 300，1 300～1 500，1 500～1 700，1 700～1 900，>1 900 m分為五級；坡度按0°～5°，5°～8°，8°～15°，15°～25°，25°～35°和>35°分為六級[19-21]。

圖2　彭陽縣坡度分級

2.2.2地形位指數結合GIS地形分析，將坡度和高程組合成一個地形位指數。該指數能夠定量分析土地利用變化與地形梯度的關系，并能夠綜合描述坡度與高程變化對土地利用變化的影響[22]，其計算公式如式(1)所示：

(1)

2.2.3地形位分布指數分布指數是用來描述不同土地利用類型在地形位梯度上的分布情況[22],分布指數可以消除地形位梯度分段和土地利用面積差異的影響。其公式如(2) 所示：

P=(Sie/Si)×(S/Se)

(2)

式中：P代表分布指數；e地形因子，分別為高程、坡度、地形位；Sie代表e地形因子特定等級下的i地類的面積；Si代表i地類的面積；Se代表整個區域e地形因子特定等級下的總面積；S代表整個區域的面積。

分布指數是一個標準化、無量綱的指數[22]。分布指數曲線越平緩，表明某種土地利用類型的分布與標準分布的偏離越小，其對地形差異的適宜性越大；反之，表明某種土地利用類型對地形具有較強的選擇性，在其優勢地形位上該組分發育較多。若P>1，表示地形位e是該土地利用類型分布的優勢地形位，P越大，優勢越明顯。反之若P<1，表示地形位e的該土地利用類型分布為非優勢地形位。在強烈的人為干擾作用下，不同時期的P值變化體現了土地利用類型在相應地形位上發生的調整。

根據公式(1)，在ArcGIS軟件平臺下，計算得到彭陽縣地形位指數，并采用等數量重分類(quantile reclassify)法[23]將地形位(0.276～1.030)劃分為13級(等級間隔為0.058)，分屬低、中、高地形位3個梯度，如表1所示。

表1　彭陽縣地形位分級

3　結果與分析

3.1基于高程分級的土地利用變化特征分析

根據高程分級圖與土地利用圖疊加得到彭陽縣1995年和2014年各土地利用類型高程分級分布，并選擇面積變化較大的耕地、草地和未利用地，各土地利用類型高程分級面積統計表，見表2。彭陽縣水域和建筑用地只分布于高程小于1 300 m的區域內，耕地面積隨高程增加而減少，這可能是因為高程高的區域不適合人們居住以及勞作。1995—2014年，耕地和未利用地面積大量減少，草地面積大量增加。高程大于1 300 m的區域，土地利用主要由草地和耕地構成。變化較為明顯的為耕地和未利用地，前者面積共減少12 778 hm2；高程大于1 700 m的區域基本不存在耕地，后者面積共減少66 916 hm2，1995年分布于各高程分級，至2014年只分布于高程小于1 300 m的高程區域內。而水域和建筑用地雖只分布于高程小于1 300 m的區域內，但面積也都有所增加。

以上彭陽縣1995—2014年土地利用類型分布在高程分級的變化體現區域生態恢復略有成效，從20世紀90年代開始，彭陽縣大力施行封山禁牧、退耕還林還草、荒山造林等生態恢復工程，導致未利用地和耕地面積大量減少，生態環境得到很好的改善，草地面積大量增加，林地和水域面積也有所增加。此外，彭陽縣經濟發展勢頭迅猛，人口的增加以及生態移民的實施，使得建筑用地面積也呈現快速增長趨勢。

3.2基于坡度分級的土地利用變化特征分析

根據坡度分級圖與土地利用圖疊加得到彭陽縣1995年和2014年各土地利用類型坡度分級分布圖，并選擇面積變化較大的耕地、草地和未利用地，各土地利用類型坡度分級面積統計表，見表3。彭陽縣水域和建筑用地的面積隨著坡度的增大而減小，而其他土地利用類型面積隨著坡度變化則沒有明顯的連續性變化。水域和建筑用地面積的增加主要發生在0°～5°的坡度分級內，占研究期間相應地類增加面積的94.75%和66.72%。建筑用地面積的增加反映20 a來彭陽縣的經濟發展；林地面積的增加主要集中在坡度小于25°的范圍內，占研究期間林地增加面積的95.06%。耕地和未利用地面積的減少體現在所有的坡度分級，大部分耕地和未利用地轉化為草地也導致了草地面積在所有坡度分級內的增加。

從1995年以來，退耕還林還草、封山禁牧、荒山造林等生態恢復工程的實施，是彭陽縣1995—2014年土地利用變化呈現坡度分級特征的最主要原因。

3.3基于地形梯度的土地利用變化特征分析

結合彭陽縣地形位指數，并分析彭陽縣1995年和2014年耕地、林地、草地、水域、建筑用地和未利用地在各地形位上的分布指數，見圖3所示。彭陽縣土地利用類型在地形位梯度上的優勢分布區間見表4。

(1) 耕地。耕地的分布指數在地形位上表現為逐漸降低的趨勢，1995年和2014年的曲線形態基本相似，其優勢度區間為1～5，而在6～13區間耕地分布指數均小于1，但在地形位2～5的區間上2014年比1995年占優勢，這主要是坡耕地面積大量減少。同時，隨著地形位的增加，耕地的分布指數開始趨向于0，表明該地形位區間內耕地的面積逐漸在減少。

表2　彭陽縣1995年與2014年各高程分級土地利用類型面積統計　hm2

表3　彭陽縣1995年與2014年各坡度分級土地利用類型面積統計　hm2

表4　彭陽縣土地利用類型在地形位梯度上的優勢分布區間

圖3　基于地形位分級的彭陽縣1995年和2014年各土地利用類型分布指數對比

(2) 林地。林地的分布指數在地形位上表現為逐漸升高的趨勢，其優勢度區間為7～13。2014年與1995年相比，在10～13地形位區間上，2014年的分布指數小于1995年的分布指數，表明林地分布趨向較高地形位區間。由于研究區域實行退耕還林還草并且調整了農業結構，減少了玉米、小麥等作物在坡耕地的種植，增加了桃、李、杏等經果林的種植。

(3) 草地。草地的分布指數在地形位上表現為逐漸升高的趨勢，且優勢分布區間最廣，地形適應彈性最大，1995年的優勢度區間為7～13而2014年的優勢度區間為5～13。相比而言，1995年草地的優勢度區間范圍較小，分布區間范圍也較小，2014年草地的優勢度區間范圍已擴展到第5區間，這與20 a來彭陽縣實施的退耕還林還草、荒山造林工程有密切關系。同時，草地優勢度區間的變化與草地面積隨坡度和高程梯度的變化相關。

(4) 水域。水域的分布指數隨著地形位的增大而減小，第1和第2區間是水域的優勢度區間。2014年與1995年相比，水域分布指數的變化較小，變化趨勢較為相似，優勢度區間沒有發生變化。而在3～13區間上分布指數接近于0，表明水域的分布對定形的選擇有嚴格的限制，同時也表明這20 a來水域所分布的地形保持穩定，沒有發生較大變化，其主要分布在低地形位區域。

(5) 建筑用地。建筑用地的分布指數隨著地形位的增大而減小，其優勢度區間為1～4。2014年與1995年相比，在6～10地形位區間，2014年的分布指數均小于1995年的分布指數。說明隨著社會經濟的發展，人口的增加，在低地形位區域，居民用地、道路、采礦等建筑用地規模在不斷擴張。此外，地形位指數從第11開始，分布指數減小至0，這主要是高地形位區間受地形限制較大，同時彭陽縣這20 a來實行了生態移民政策，將處于山區，交通不便，環境惡劣的居民遷移出來，導致較高地形位上沒有建筑用地分布。

(6) 未利用地。未利用地的分布指數曲線1995年和2 014存在較為明顯的不同。1995年的分布指數曲線隨著地形位的增大，分布指數先增大后減小，其優勢度區間為6～11，而2014年的分布指數曲線變化相對劇烈，特別是從第9地形位區間開始，分布指數分布指數直接降到1以下，其優勢度區間范圍也縮小為6～9。可以看出，彭陽縣這20 a來土地利用類型的變化，最主要便是集中在對未利用地的處理，通過荒山造林等生態恢復工程，將高地形位區間的未利用地大量轉化為草地，使得2014年的未利用地的優勢度區間范圍較小，分布區間范圍也較小。

除未利用地外，其他土地利用類型1995年和2014年的分布指數在地形位分布上較為一致。水域的優勢度區間集中在低地形位區間，耕地和建筑用地的優勢度區間集中于低中地形位區間，林地、草地和未利用地的優勢度區間集中于中高地形位區間。1995—2014年，耕地、水域和建筑用地的地形位分布指數趨向減少，其分布趨向于較低的地形位區間，而林地、草地和未利用地則趨向于較高的地形位區間。

4　討論與結論

4.1討 論

(1) 本文從高程，坡度單個因子以及地形位梯度兩個方面，對彭陽縣1995—2014年土地利用類型在地形梯度上的分布及變化特征進行研究分析，不難發現運用高程、坡度以及地形梯度都能夠反映研究區域20 a來土地利用的變化，并且總體上各土地利用類型在高程、坡度上的變化與在地形梯度上的變化一致，但是地形位更能反映土地利用在空間格局上的變化規律，例如近20 a來草地面積在高程帶和坡度帶上均有大量增加但總體變化較為復雜，不過在地形位分布上呈現一定的規律性。

(2) 黃土丘陵溝壑區是生態環境極端脆弱區，經濟貧困落后，該區域土地資源的可持續利用關系到區域生態安全，因此需要加強退耕還林草工程、荒山造林、水土保持生態修復、生態移民等工程引起的土地利用變化的影響研究，對評估區域土地資源可持續利用和生態環境重建具有重要指導意義。彭陽縣20 a的土地利用變化反映該區域生態恢復的重點與變化規律，彭陽縣土地利用/覆被景觀格局變化的時間節點為2000年，從土地利用數據充分反映了國家的退耕還林草工程、水土保持生態修復、生態移民等土地可持續利用政策成為土地利用變化的重要影響因素。2007年國家對退耕還林還草工程結構性調整，和水土保持工作的支持，實現水土資源可持續利用和生態環境可持續發展，保障和改善民生。

(3) 本文只是利用高程、坡度、地形位指數以及土地利用程度綜合指數定量分析了在地形梯度上的土地利用變化特征，并沒有考慮坡向以及其他地形因子的影響，綜合更多的地形因子對土地利用變化特征進行探究，并定量分析其變化是接下來可以繼續研究的方向。

4.2結 論

(1) 1995—2014年，彭陽縣耕地和未利用地面積大量減少，分別減少12 778，66 916 hm2，草地面積大量增加，增加74 973 hm2，林地、水域和建筑用地面積均有不同程度的增加，分別為3 873,438，1 208 hm2。

(2) 土地利用變化與高程的關系。1995年到2014年，彭陽縣水域和建筑用地只分布于高程小于1 300 m的區域內，耕地面積隨高程增加而減少；高程大于1 300 m的區域，土地利用類型主要為草地和耕地；高程大于1 700 m的區域基本無耕地；1995年耕地分布于各高程分級，2014年主要分布于高程小于1 300 m的高程區域內。(3) 土地利用變化與坡度的關系。1995—2014年，彭陽縣水域和建筑用地的面積隨著坡度的增大而減小，其他土地利用類型面積隨著坡度變化則沒有明顯的連續性變化；由于坡度較大的區域并不適宜人們居住，因此1995年到2014年彭陽縣建筑用地面積的增加集中在0°～5°的坡度分級內；林地面積的增加主要集中在坡度小于25°的范圍內；耕地和未利用地面積的減少貫穿了所有的坡度分級；大部分耕地和未利用地轉化為草地也導致了草地面積在所有坡度分級內的增加。

(4) 土地利用變化與地形位指數的關系。1995—2014年,隨著地形位的升高，耕地、水域、建筑用地的分布指數逐漸降低；林地和草地的分布指數逐漸升高；未利用地的分布指數先升后降；草地的優勢分布區間增加了2個地形位，未利用地的優勢分布區間降低了2個地形位。整體而言，水域的優勢度區間集中在低地形位區間，耕地和建筑用地的優勢度區間集中于低中地形位區間，林地、草地和未利用地的優勢度區間集中于中高地形位區間。研究期內，耕地、水域和建筑用地的地形位分布指數趨向減少，其分布趨向于較低的地形位區間，而林地、草地和未利用地則趨向于較高的地形位區間。

[1]Lambin E F, Baulies X, Bockstael N, et al. Land-use and land-cover change (LUCC): implementation strategy[R]. IGBP Rep. 48, IHDP Rep. 10, Int, 1999.

[2]Liu Jiyuan, Liu Mingliang, Zhuang Dafang, et al. Study on spatial pattern of land use change in China during 1995—2000[J]. Science in China：Series D, 2003,46(4):373-384.

[3]Erickson D L. Rural land use and land cover change implications for local planning in the fiver Raisin watershed[J]. Land Use Police, 1995,12(3):223-236.

[4]劉紀遠,布和敖斯爾.中國土地利用變化現代過程時空特征的研究：基于衛星遙感數據[J].第四紀研究,2000,20(3):229-239.

[5]劉紀遠,張增祥,莊大方,等.20世紀90年代中國土地利用變化時空特征及其成因分析[J].地理研究,2003,22(1):1-12.

[6]Liu J Y, Liu M L, Zhuang D F, et a1. Study on spatial pattern of land use change in China during 1995—2000[J]. Science in China,2003,48(4):373-384.

[7]趙堅,曹春香,賈慧聰,等.干旱區土地利用/覆蓋時空動態變化分析:以準格爾旗為例[J].干旱區資源與環境,2012,26(8):93-96.

[8]呂靜,常慶瑞.黃土高原土石山區土地利用空間分布研究[J].干旱區資源與環境,2012,26(8):93-96.

[9]崔衛國,文倩,劉艷艷,等.基于DEM的醴陵市土地利用空格局分析[J].資源科學,2008,30(2):228-234.

[10]崔步禮,李小雁,姜廣輝,等.基于DEM的山地丘陵區土地利用/覆被研究:以青海湖流域為例[J].自然資源學報,2011,26(5):871-880.

[11]鄒敏,吳泉源,逄杰武.基于DEM的龍口市土地利用空間格局與時空變化研究[J].測繪科學,2007,32(6):174-175.

[12]湯國安,宋佳.基于DEM坡度圖制圖中坡度分級方法的比較研究[J].水土保持學報,2006,20(2):157-162.

[13]馬士彬,張勇榮,安裕倫.山區城市土地利用動態空間分布特征:以貴州省六盤水市為例[J].自然資源學報,2012,27(3):489-496.

[14]吳學偉,龔文峰,袁力.基于DEM黑龍江省土地利用定量分析[J].東北林業大學學報,2009,37(9):74-75.

[15]李增加,馬友鑫,李紅梅,等.西雙版納土地利用/覆被變化與地形的關系[J].植物生態學報,2008,32(5):1091-1103.

[16]崔步禮,李小雁,姜廣輝,等.基于DEM的山地丘陵區土地利用覆被研究:以青海湖流域為例[J].自然資源學報,2011,26(5):871-880.

[17]邱揚,傅伯杰,王軍,等.黃土丘陵小流域土地利用的時空分布及其地形因子的關系[J].自然資源學報,2003,18(1):20-29.

[18]陳楠,楊武年,李娟,等.巴中市丘陵山區土地利用在地形梯度上的分布特征[J].水土保持通報,2012,32(1):185-188.

[19]李壁成.小流域水土流失與綜合治理遙感監測[M].北京:科學出版社,1995.

[20]范紅霞.紙坊溝流域景觀格局動態及地形分異研究[D].陜西楊凌：西北農林科技大學,2011.

[21]井學輝.新疆額爾齊斯河流域植被景觀格局與生物多樣性空間變化規律研究[D].北京：中國林業科學研究院,2008.

[22]喻紅,曾輝,江子瀛,等.快速城市化地區景觀組分在地形梯度上的分布特征研究[J].地理科學,2001,21(2):64-69.

[23]李升發,董玉祥.廣州市土地利用垂直分層模式及其演變過程[J].自然資源學報,2013,28(6):1019-1029.

Relationship Between the Land Use Change and the Terrain Gradient in Counties in the Loess Hilly and Gully Region

ZHUANG Zhuzhou1, HUANG Qiuhao1, SHI Yun2

(1.CollegeofGeographicandOceanographicSciences,NanjingUniversity,Nanjing210023,China;2.CollegeofResourcesandEnvironment,NingxiaUniversity,Yinchuan750021,China)

Land-use/cover change (LUCC) is the most direct manifestation of the interaction between human activities and natural environment, and is an important factor to reflect the extent of human activities as well. Based on the DEM data in 1995 and 2014 and the land use status, we quantitatively analyzed the characteristics of land use change based on terrain gradient in Pengyang County with the aid of GIS space analysis and terrain niche index, dispersion index, and the comprehensive index. Firstly, from 1995 to 2014, the area of cultivated land and unused land decreased dramatically. The grassland area, on the other hand, increased significantly, and the area of woodland, water bodies and construction land also showed the increasing trend. Secondly, with the elevation of the terrain niche, the distribution index of cultivated land, water bodies and construction land declined, while areas of woodland and grassland were rising, and the distribution index of unused land first rose then declined. The dominance distribution section of grassland increased two terrain niches, while that of unused land decreased two. Overall, water bodies′ dominance distribution section concentrates at the low position; cultivated land and construction land dominance distribution sections are mainly at the low-to-mid position and woodland, grassland and unused land dominance distribution sections are mainly at high position. Thirdly， compared with the comprehensive index of the extent of land use in 1995, that increased in all terrain niches in 2014. This research can provide the theoretical basis for the optimization of land use pattern and sustainable utilization of land resources in Pengyang County.

land use change; terrain niche; distribution index; Pengyang County

2015-07-30

2015-08-25

國家自然科學基金“基于景觀格局變化的黃土高原縣域退耕還林還草生態效益評價研究”(41161081)

莊逐舟(1992—),男,福建泉州人,研究生,從事GIS應用方面研究。E-mail:zhuangzhuzhou@163.com

石云(1977—),女,寧夏銀川人,博士,副教授,長期從事GIS應用教學及研究。E-mail:shiyun@nxu.edu.cn

F301

A

1005-3409(2016)04-0331-07