基于地形梯度的淇河流域土地利用時(shí)空變異分析

馮小燕，朱文博，張靜靜，李理，徐帥博，朱連奇

（河南大學(xué)環(huán)境與規(guī)劃學(xué)院，河南 開封 475004）

山區(qū)是特殊的國土空間單元，具有較強(qiáng)的資源梯度性、災(zāi)害多發(fā)性和生態(tài)脆弱性等特點(diǎn)；不合理的山區(qū)土地利用將導(dǎo)致植被破壞和嚴(yán)重的水土流失[1-2]。而地形通過立地條件影響陸地生態(tài)系統(tǒng)物質(zhì)與能量的再分配，并通過干擾發(fā)生頻率及強(qiáng)度成為影響山區(qū)土地利用的主導(dǎo)因子[3]。因山區(qū)土地和環(huán)境因子的多樣及復(fù)雜性、生態(tài)系統(tǒng)對人類擾動(dòng)的敏感性，二者之間的耦合更成為學(xué)界關(guān)注的熱點(diǎn)[4-11]。地形是決定土地利用的重要環(huán)境因子，目前多見于單一地形因子對土地利用/覆被類型的影響，比如從地形海拔高程角度結(jié)合土地利用程度綜合指數(shù)對不同垂直帶土地利用類型分異程度進(jìn)行分析；基于高程和坡度，利用景觀研究的分布指數(shù)分析不同土地利用類型時(shí)空變化的規(guī)律等[12]?；贒EM和遙感技術(shù)、地理信息系統(tǒng)方法把地形對土地利用格局的影響進(jìn)行分析，揭示了山區(qū)土地利用的多維變化規(guī)律，指出地形構(gòu)成了土地利用格局的基本框架，分析了山區(qū)土地利用變化的總體趨勢[13-26]。但已有研究對不同地形梯度土地利用類型轉(zhuǎn)移的定量化研究相對薄弱[27-28]。而土地利用與覆被變化的研究對于山地生態(tài)系統(tǒng)的穩(wěn)定及服務(wù)有重要影響，如何確定不同垂直梯度的土地利用結(jié)構(gòu)，合理管控山地不同海拔高程土地利用類型的轉(zhuǎn)化是穩(wěn)定山區(qū)土地資源、保護(hù)其利用的可持續(xù)性，使山地生態(tài)系統(tǒng)服務(wù)價(jià)值最大化的理論基礎(chǔ)。因此開展基于地形梯度的山區(qū)土地利用時(shí)空變異規(guī)律研究就具有理論意義和現(xiàn)實(shí)意義。

本文選擇具有明顯地形過渡特征的淇河流域，綜合選取多種地形因子，基于1990、2000、2010和2015年四期土地利用/覆被變化（LUCC）數(shù)據(jù)，運(yùn)用轉(zhuǎn)移矩陣、區(qū)域統(tǒng)計(jì)分析和地形分布指數(shù)等方法，系統(tǒng)分析淇河流域土地利用時(shí)空變異特征，探討地形因子對土地利用格局的分異作用以及土地利用動(dòng)態(tài)變化的地形梯度差異特征，以期為該區(qū)土地利用規(guī)劃與管理以及生態(tài)環(huán)境建設(shè)與保護(hù)提供決策支持。因此研究國土空間格局的變化規(guī)律，對促進(jìn)國土空間有序、適度、可持續(xù)開發(fā)利用，強(qiáng)化國土空間治理極為重要。

1 研究區(qū)概況

淇河是海河水系衛(wèi)河的支流，發(fā)源于山西省陵川縣，經(jīng)山西壺關(guān)縣、河南省的輝縣市、林州市、鶴壁市、淇縣、浚縣，匯入衛(wèi)河，全長約165 km。淇河流域（113°17′～114°23′E，35°32′～36°04′N）往東與華北平原相接（圖1），流域面積2227 km2，90%的區(qū)域?qū)儆谏絽^(qū)。海拔自西向東遞減，范圍為31～1862 m，平均海拔為781 m。地貌類型復(fù)雜，主要有山地、低山丘陵和平原等。屬暖溫帶半濕潤季風(fēng)氣候，四季分明，平均降水量為600 mm左右，年均氣溫為12 ℃左右。太行山作為京津和華北平原的生態(tài)屏障，使得位于其南段的淇河流域生態(tài)環(huán)境也比較脆弱，因此該區(qū)急需加強(qiáng)土地利用規(guī)劃與管理。

2 研究方法與數(shù)據(jù)來源

2.1 研究方法

2.1.1 轉(zhuǎn)移矩陣 引入轉(zhuǎn)移矩陣來分析各土地利用類型之間的相互轉(zhuǎn)化情況，進(jìn)而揭示土地利用變化的內(nèi)部過程及趨勢。公式如下[29]：

式中：Sij為地類面積；n為地類的類型數(shù)；i和j分別為研究時(shí)段初期和末期的地類序號(hào)。根據(jù)公式可計(jì)算不同地形因子梯度上各地類之間的相互轉(zhuǎn)換特征。

基于1990、2000、2010和2015年四期土地利用數(shù)據(jù)，分析淇河流域土地利用多維變化，構(gòu)建1990—2000，2000—2010、2010—2015 和 1990—2015年土地利用轉(zhuǎn)移矩陣，分析1990—2015年淇河流域土地利用的變化方向。并從不同海拔、坡度和坡向的視角出發(fā)，提取不同土地利用現(xiàn)狀及變化在不同海拔、坡度和坡向的分布，得到淇河流域土地利用結(jié)構(gòu)及變化的地形梯度特征。

2.1.2 地形分布指數(shù) 為消除不同地形梯度上面積差異帶來的量綱影響，本研究引入地形分布指數(shù)來分析土地利用在各地形因子梯度上的分布特征，分布指數(shù)[30]是一個(gè)標(biāo)準(zhǔn)化、無量綱指數(shù)，公式如下：

式中：P為分布指數(shù)；Sie為e地形梯度上i種地類面積；Si為i種地類的總面積；Se為e地形因子的總面積；S為研究區(qū)的面積。P值的大小代表某地類出現(xiàn)的頻率，當(dāng)P＞1時(shí)，表示i地類在e地形梯度上的分布為優(yōu)勢分布，值越大，代表優(yōu)勢度越高。根據(jù)公式可計(jì)算土地利用類型以及轉(zhuǎn)移矩陣得到的地類之間轉(zhuǎn)換情況在各地形因子分布的優(yōu)勢區(qū)間。

2.2 數(shù)據(jù)來源

本研究所用1990、2000和2010年土地利用/覆被變化（LUCC）數(shù)據(jù)來源于國家地球系統(tǒng)數(shù)據(jù)共享平臺(tái)-黃河中下游科學(xué)數(shù)據(jù)中心（http://www.geodata.cn/），2015年的土地利用/覆被變化（LUCC）是基于2015年的Landsat 8遙感影像 (來源于地理空間數(shù)據(jù)云，http://www.gscloud.cn/)，參考該區(qū)往年土地利用類型圖，進(jìn)行人工目視解譯所得。土地利用數(shù)據(jù)均為1∶10萬的矢量圖，包括25個(gè)二級(jí)類，根據(jù)需要，合并為6個(gè)一級(jí)類，分別為耕地、林地、草地、水域、建設(shè)用地及未利用地。

2.3 數(shù)據(jù)處理

地形數(shù)據(jù)包括海拔、坡度及坡向，均由DEM計(jì)算得到，DEM數(shù)據(jù)來源于ASTER GDEM數(shù)據(jù)，分辨率為30 m。依據(jù)研究區(qū)地形特點(diǎn)和土地利用分布情況，將生成的高程圖、坡度圖和坡向圖重新進(jìn)行分級(jí)[28]，利用ArcGIS中的重分類工具對海拔、坡度和坡向進(jìn)行分級(jí)，將海拔劃分為＜200 m、200～400 m、400～600 m、600～800 m、800～1 000 m、1 000～1200 m、1200～1400 m、1400～1600 m、＞1600 m共九級(jí) ；將坡度劃分為 0°～5°、5°～8°、8°～15°、15°～25°、25°～35°、＞35°共六級(jí)，坡度分級(jí)時(shí)結(jié)合了研究區(qū)的地形地類特點(diǎn)，采用這樣分級(jí)范圍能夠體現(xiàn)出研究區(qū)坡度梯度明顯的差異性；將坡向劃分為無坡向（-1°）、北坡（0°～22.5°，337.5°～360°）、東北坡（22.5°～67.5°）、東坡（67.5°～112.5°）、東南坡（112.5°～157.5°）、南坡（157.5°～202.5°）、西南坡（202.5°～247.5°）、西 坡（247.5°～292.5°） 和 西 北 坡（292.5°～337.5°），研究區(qū)地形因子分級(jí)及面積比例見圖2和表1。

圖2 研究區(qū)海拔、坡度和坡向的空間分布Fig. 2 Spatial distributions of elevation, slope and slope direction in the study area

表1 海拔、坡度、坡向分級(jí)及面積比例Table 1 Classi fi cation and area proportions of elevation,slope and slope direction

3 結(jié)果與分析

3.1 土地利用整體變化

由表2和圖3可以看出1990—2000年淇河流域土地利用變化幅度不大，2000—2010年和2010—2015年變化幅度較劇烈；25年來變化最大的為耕地和草地，耕地面積凈減少了137.19 km2，占6.19%，草地面積凈增加了88.39 km2，占3.97%。建設(shè)用地面積凈增加了49.21 km2，占2.21%。水域面積凈增加了20.16 km2，占0.91%。林地面積和未利用地面積變化不大，分別減少了0.96%和增加了0.03%。近年來淇河流域退耕還林還草、生態(tài)保護(hù)區(qū)建設(shè)、濕地保護(hù)等政策的實(shí)施使得草地和水域面積大幅增加。

表2 淇河流域不同年份土地利用變化Table 2 Land use change in Qi River Basin among different years

圖31990、2000、2010和2015年土地利用類型圖Fig. 3 Land use types in 1990, 2000, 2010 and 2015

由圖4可知，1990—2000年間耕地面積減少了7.49 km2，主要轉(zhuǎn)為草地、建設(shè)用地和林地，少部分轉(zhuǎn)為了水域，同時(shí)也有部分耕地和林地轉(zhuǎn)變?yōu)椴莸?，草地變化不大；林地增加量略大于減少量，未利用地和水域基本不變；建設(shè)用地增加的主要來源是耕地。2000—2010年土地利用變化幅度相比前十年明顯增大，耕地變化面積最大，減少72.30 km2；林地增加量和減少量相當(dāng)，主要是與耕地和草地相互轉(zhuǎn)化；水域面積呈增加趨勢；草地增加明顯，主要來源是耕地和林地；建設(shè)用地主要來源是耕地，其次是草地和林地。2010—2015年各用地類型相互轉(zhuǎn)化最明顯，主要是耕、林、草之間的轉(zhuǎn)化。草地面積動(dòng)態(tài)變化最大，呈現(xiàn)增長態(tài)勢；其次為林地，減少量多于增加量；耕地面積主要轉(zhuǎn)化為草地和林地，也有部分草地轉(zhuǎn)化為耕地，草地和林地與耕地之間相互轉(zhuǎn)化，但耕地面積減少量明顯大于增加量；建設(shè)用地增加的主要來源是草地，小部分來源于耕地。1990—2015年耕地變化面積最大，25年來減少137.19 km2，主要轉(zhuǎn)化為草地、林地和建設(shè)用地；草地增加幅度較大，25年來增加88.39 km2；林地呈減少趨勢，水域呈增加趨勢；建設(shè)用地主要來源是耕地，其次為草地。

3.2 基于地形梯度的土地利用結(jié)構(gòu)變化

圖41990—2015年淇河流域土地利用轉(zhuǎn)移方向Fig. 4 Transformation of land uses in Qi River Basin from 1990 to 2015

圖5 土地利用結(jié)構(gòu)的海拔梯度分布Fig. 5 Elevation gradient distribution of land use structure

3.2.1 土地利用結(jié)構(gòu)的海拔梯度差異 圖5為淇河流域1990年和2015年不同海拔分級(jí)土地利用結(jié)構(gòu)：耕地主要分布在400 m海拔以下地區(qū)，1990年和2015年面積變化不大，分別占51.98%和50.49%；在1200～1400 m等級(jí)再次出現(xiàn)峰值，分別占比16.11%和15.57%。林地主要分布在1200～1600 m海拔區(qū)，1990年和2015年略有差別，分別占46.13%和49.51%，在低于400 m和大于1600 m海拔區(qū)占比重均較小。草地主要分布在200～600 m海拔區(qū)，1990年和2015年分別占54.34%和52.12%，在1400—1600 m等級(jí)再次出現(xiàn)峰值，分別占10.34%和11.92%；水域和建設(shè)用地分布相似，均在低于400 m海拔區(qū)域所占面積最大；在低于400 m海拔區(qū)，1990年和2015年水域分布差別較大，分別占水域總面積的63.44%和50.69%；在低于400 m海拔區(qū)，1990年和2015年建設(shè)用地分布差別較明顯，分別占建設(shè)用地總面積的75.18%和66.13%。

從圖6淇河流域1990年和2015年不同坡度分布可知：耕地在0°～5°等級(jí)的分布比重最大，1990年和2015年的占比分別為31.43%和30.45%；在8°～15°等級(jí)達(dá)到第二個(gè)峰值，分別占27.52%和26.87%。林地和草地在坡度分布上類似，隨坡度增大比重呈先增后減，均在15°～25°等級(jí)達(dá)到峰值，1990年和2015年比重變化不明顯，比重均在31%左右；水域和建設(shè)用地分布較為類似，隨坡度增大比重不斷降低，均在0°～5°等級(jí)的面積比重最大，水域在1990年和2015年的比重差別較大，分別為43.78%和32.41%；在0°～5°的坡度等級(jí)中，1990年和2015年建設(shè)用地的占比差距較大，分別占比50.02%和43.23%。

從圖7可以得出：耕地在東南坡、南坡、東坡、東北坡分布比例較大，即主要分布在陽坡、半陽坡；林地在東坡、東北坡、北坡和西北坡、東南坡分布較多，即半陽坡、半陰坡和陰坡；草地在東坡、東南坡、東北坡和西北坡分布較多，即多分布于半陽坡和半陰坡；水域明顯朝向東北坡、東坡、東南坡、南坡分布。

3.2.2 土地利用變化的海拔梯度差異 由圖8可知各土地利用類型在不同海拔梯度的變化情況：耕地在海拔200 m以下、200～400 m和1200～1600 m區(qū)域面積呈減少趨勢，其中在200～400 m和1400～1600 m區(qū)域內(nèi)的減少達(dá)到兩個(gè)峰值，減少量最大，分別減少63.21 km2和35.13 km2；在200 m以下平坦寬闊的平原區(qū)，耕地減少的原因是人為開發(fā)建設(shè)占用，導(dǎo)致耕地轉(zhuǎn)變?yōu)榻ㄔO(shè)用地；在200～400 m和1200～1600 m高程區(qū)域內(nèi)，耕地減少是因?yàn)樯絽^(qū)耕地分散零碎，交通不便，耕作成本高，經(jīng)濟(jì)效益低下和國家的退耕還林還草政策；耕地在400～1200 m的區(qū)域面積內(nèi)呈穩(wěn)定狀態(tài)。林地面積增加主要在海拔1200～1600 m區(qū)域，增加了12.41 km2；在400～1200 m面積呈減少趨勢，減少了34.20 km2；在200～400 m區(qū)域內(nèi)呈穩(wěn)定狀態(tài)。草地在海拔400～600 m面積減少，減少面積為16.00 km2；在200～400 m、600～800 m 和 1 000～1600 m 高程級(jí)別上面積增加，其中在200～400 m區(qū)域內(nèi)增加的最多，達(dá)到峰值狀態(tài)，增加了46.76 km2。水域和建設(shè)用地增加集中在低海拔地區(qū)，隨海拔增大比重呈減少趨勢。

圖6 土地利用結(jié)構(gòu)的坡度梯度分布Fig. 6 Slope gradient distribution of land use structure

圖7 土地利用結(jié)構(gòu)的坡向分布Fig. 7 Slope direction distribution of land use structure

圖81990—2015年土地利用變化的海拔梯度分布Fig. 8 Elevation gradient distribution of land use change from 1990 to 2015

由圖9可知各土地利用類型在不同坡度梯度的變化情況：耕地在0°～25°區(qū)域內(nèi)呈現(xiàn)減少趨勢，其中在 0°～5°和 8°～15°區(qū)域出現(xiàn)峰值，分別減少了49.00 km2和41.53 km2。林地增加面積主要在低于15°區(qū)域，在8°～15°范圍增加最多，增加了16.32 km2；減少面積主要在高于15°區(qū)域，在15°～25°坡度級(jí)別上減少的最多，減少面積為5.23 km2。草地在0°～35°區(qū)域內(nèi)呈現(xiàn)增加趨勢，其中在 15°～25°區(qū)域出現(xiàn)變化的峰值，增加最多，增加面積為8.72 km2，增加面積始終多于減少面積；水域在0°～35°區(qū)域內(nèi)呈現(xiàn)增加狀態(tài)，增加面積主要集中在8°～25°區(qū)域；建設(shè)用地在0°～35°區(qū)域內(nèi)始終呈現(xiàn)增加狀態(tài)，且增加集中在 0°～15°地區(qū)，在 0°～5°區(qū)域增加值最大，增加面積為16.34 km2。

圖91990—2015年土地利用變化的坡度梯度分布Fig. 9 Slope gradient distribution of land use change from 1990 to 2015

3.2.3 土地利用類型分布指數(shù)變化特征 由圖10可知土地利用類型在海拔梯度上的優(yōu)勢位分布：1990年和2015年兩個(gè)時(shí)期的耕地分布指數(shù)曲線形態(tài)基本相似，其分布指數(shù)基本保持穩(wěn)定。優(yōu)勢分布區(qū)間沒有變化，均為＜400 m和1200～1400 m高程等級(jí)。通過前述分析可知，25年間，耕地面積除在＜400 m、1200～1600 m和＞1600 m等級(jí)降幅明顯，在其余的高程帶呈現(xiàn)增長狀態(tài)，而分布指數(shù)消除了研究區(qū)耕地面積減少帶來的差異，說明在高程梯度上耕地面積的增減幅度和分布指數(shù)變化基本一致。林地在不同的高度梯度上具有很強(qiáng)的適應(yīng)性，其空間分布格局的顯著變化是優(yōu)勢分布范圍逐漸向高地形梯度縮小，1990年其優(yōu)勢位分布在600～1600 m和＞1600 m高程梯度上，2015年其優(yōu)勢位分布在600～1600 m區(qū)域，與1990年相比，2015年的優(yōu)勢位減少了1個(gè)高程梯度。草地的優(yōu)勢分布區(qū)間沒有明顯變化，1990和2015年的優(yōu)勢區(qū)間均在200～1 000 m高程級(jí)別上，集中分布于中高地形梯度區(qū)；在1 000～1600 m地形梯度區(qū)，草地分布指數(shù)有所上升。在1990年和2015年兩個(gè)時(shí)期中，水域的優(yōu)勢區(qū)間上升了兩個(gè)地形梯度。研究期間建設(shè)用地優(yōu)勢位保持不變，均為＜200 m和200～400 m高程級(jí)別；在9個(gè)地形梯度中，2個(gè)時(shí)期的分布指數(shù)曲線相似，除＜200 m和200～400 m地形梯度間，2015年的分布指數(shù)均大于1990年，說明地形因素對人類開發(fā)建設(shè)活動(dòng)的影響在減少。

圖10 1990—2015年土地利用類型在海拔梯度上的分布指數(shù)Fig. 10 Elevation gradient distribution index of land use types from 1990 to 2015

圖111990—2015年土地利用類型在坡度梯度上的分布Fig. 11 Slope gradient distribution index of land use types from 1990 to 2015

從圖11可知土地利用類型在坡度梯度上的優(yōu)勢位分布：1990年和2015年兩個(gè)時(shí)期耕地的分布指數(shù)曲線形態(tài)基本一致，其分布指數(shù)保持穩(wěn)定；優(yōu)勢分布區(qū)間沒有變化，均為 0°～5°、5°～8°和 8°～15°區(qū)域，耕地面積在0°～15°降幅明顯，而分布指數(shù)消除了研究區(qū)耕地面積減少帶來的差異，說明在坡度梯度上耕地面積的增減幅度和分布指數(shù)變化基本一致。林地在不同的坡度梯度上具有很強(qiáng)的適應(yīng)性，其優(yōu)勢位分布于 15°～25°、25°～35°、＞35°坡度區(qū)間；在0°～15°坡度區(qū)間，2015年分布指數(shù)高于1990年；在 15°～35°、＞35°坡度區(qū)間，2015 年的分布指數(shù)明顯低于1990年。草地的優(yōu)勢分布區(qū)間沒有明顯變化，1990年優(yōu)勢區(qū)間為15°～35°坡度級(jí)別，2015年的優(yōu)勢區(qū)間為15°～35°和＞35°區(qū)域，集中分布于中高地形梯度區(qū)；1990年到2015年，在8°～25°地形梯度區(qū)，草地分布指數(shù)有所降低。在1990年和2015年間，水域的優(yōu)勢位均為 0°～5°和5°～8°坡度梯度；在8°～35°、＞35°坡度范圍內(nèi)，2015年的分布指數(shù)始終高于1990年的分布指數(shù)。研究期間建設(shè)用地優(yōu)勢位保持不變，均為0°～5°和5°～8°之間；在 6個(gè)地形梯度中，2個(gè)時(shí)期的分布指數(shù)曲線相似，在0°～8°坡度分布區(qū)間中，1990年分布指數(shù)高于2015年，在8°～35°和＞35°坡度區(qū)間中2015年的分布指數(shù)均大于1990年。

4 結(jié)論與討論

4.1 結(jié)論

1）淇河流域土地利用總體變化特征以耕地減少為主。2000—2010年和2010—2015年變化幅度較劇烈；25年來變化最大的為耕地和草地，耕地面積凈減少最多，草地凈增加最多，建設(shè)用地面積凈增加次之，水域面積凈增加最少。林地面積和未利用地面積變化不大。近年來淇河流域退耕還林還草、生態(tài)保護(hù)區(qū)建設(shè)、濕地保護(hù)等政策的實(shí)施使得草地和水域面積大幅增加。而建設(shè)用地?cái)U(kuò)張明顯，淇河流域人地矛盾依然嚴(yán)峻。

2）淇河流域土地利用結(jié)構(gòu)具有顯著的地形梯度差異。就土地利用結(jié)構(gòu)的地形梯度分布而言，耕地、建設(shè)用地、草地、水域多位于低海拔和小坡度地區(qū)，林地處于高海拔地區(qū)。土地利用變化在地形梯度上的分布：在低海拔區(qū)域，耕地面積減少，草地、水域和建設(shè)用地面積增加。在中海拔區(qū)域，草地和林地面積減少。林地面積增加主要在高海拔區(qū)域。土地利用變化在坡度梯度上的分布：在低坡度區(qū)域，耕地面積減少，建設(shè)用地和林地面積增加；草地和水域始終處于增加狀態(tài)。

3）土地利用類型優(yōu)勢位在海拔梯度和坡度梯度保持穩(wěn)定，變化微小。耕地、草地和建設(shè)用地在海拔梯度和坡度梯度的優(yōu)勢位變化微小。在海拔梯度上，林地的空間分布格局發(fā)生了顯著變化，優(yōu)勢分布范圍逐漸向高地形梯度縮小。

4.2 討論

淇河流域處于山地與平原的過渡地帶，區(qū)域內(nèi)地形分異較大，通過研究區(qū)域內(nèi)土地利用格局的地形梯度特征，可以進(jìn)一步分析不同地形條件下土地利用方式的分布選擇。無量綱的分布指數(shù)消除了不同地形梯度面積分布比例不同帶來的影響，更能反映出土地利用類型的優(yōu)勢分布情況。分析發(fā)現(xiàn)土地利用方式有著明顯的地形梯度特征分異，低海拔區(qū)間和低坡度區(qū)間高程低，坡度平緩，氣候條件適宜，人為活動(dòng)對土地利用方式干擾性大，為建設(shè)用地、耕地和水域的優(yōu)勢分布區(qū)間，高海拔區(qū)間和高坡度區(qū)間人類對土地利用方式干擾度低，為林地的優(yōu)勢分布區(qū)間，符合山地平原過渡帶的自然規(guī)律和社會(huì)經(jīng)濟(jì)規(guī)律，地形因子很大程度決定了土地利用格局。研究中沒有考慮地形起伏度及其他地形因子的影響，并只做了定性的驅(qū)動(dòng)力分析，探究綜合更多的地形因子對土地利用格局及其變化的影響并定量分析變化的驅(qū)動(dòng)力是下一步工作的重點(diǎn)。