基于轉移矩陣的準格爾旗土地利用變化分析

張建國, 李晶晶,2, 殷寶庫, 高云飛, 劉曉燕

(1.黃河上中游管理局, 陜西 西安710021; 2.西北農(nóng)林科技大學 水土保持研究所,陜西 楊凌 712100; 3.黃河水利委員會, 河南 鄭州 450003)

土地利用變化類因子是影響水土流失變化的重要類型因子，了解土地利用變化是分析區(qū)域水土流失強度和面積變化的基礎。目前比較常用的土地利用遙感動態(tài)監(jiān)測信息的分析方法主要包括傳統(tǒng)的數(shù)理統(tǒng)計法、構造參數(shù)法、土地利用轉移矩陣法以及景觀格局特征分析法[1]。土地利用轉移矩陣在土地利用變化分析中具有重要作用，并得到較為廣泛的應用[2-4]。相對于其他分析方法，其科學性及優(yōu)勢性在于矩陣的生成靈活多變，分析的結果不受地物類別、數(shù)量限制，同時又可以對多時期的數(shù)據(jù)進行分析處理[5]。土地利用轉移矩陣是將土地利用變化的類型轉移面積以矩陣的形式列出，作為用地結構與變化方向分析的基礎，可以直觀地反映出某一區(qū)域土地利用面積的變化，也可細致反映各地類之間的相互轉化關系及來源概率，進而了解轉移前后各地類的結構特征，在土地利用變化分析和地類模擬方面具有重要作用。

準格爾旗位于毛烏素沙漠和西北黃土高原過渡地帶的沙化區(qū)，屬于陜晉蒙農(nóng)牧交錯區(qū)，生態(tài)環(huán)境脆弱，水土流失嚴重。研究該區(qū)土地利用類型變化規(guī)律，對于了解該土地利用變化驅動因素和水土流失變化研究具有十分重要的意義。本文在遙感與GIS技術的支持下，通過引入土地利用轉換矩陣，對準格爾旗兩不變年度的土地利用格局進行分析，獲取土地利用變化特征，并結合土壤侵蝕數(shù)據(jù)，以期分析該區(qū)域水土流失變化情況，為該地區(qū)土地可持續(xù)利用與管理提供基礎數(shù)據(jù)與決策支持。

1 研究區(qū)概況及信息源

準格爾旗位于內(nèi)蒙古西南部、鄂爾多斯高原東部，地處黃土高原邊緣，是晉陜蒙金三角地帶，位于東經(jīng)110°05′—111°27′，北緯39°16′—40°21′，南北長116.5 km，東西寬115.2 km，總面積7 692 km2。

研究數(shù)據(jù)主要來源于2013—2017年全國水土流失動態(tài)監(jiān)測與公告項目黃河流域監(jiān)測成果：準格爾旗2014和2016年遙感影像及土地利用和土壤侵蝕解譯成果。土地利用主要包括耕地、園地、林地、草地、居民點及工礦交通用地、水域及水利設施用地及其他土地這7個1級地類，其中其他土地包括沙地和裸巖裸地等地類。信息源以天繪1號遙感影像為主，輔以高分1號影像和資源3號衛(wèi)星影像。融合影像包括紅、綠、藍3個波段，影像分辨率達到2 m。準格爾旗水蝕區(qū)坡度信息由1∶1萬DEM提取，風蝕區(qū)坡度信息由1∶5萬DEM提取。高斯—克呂格(Gauss-Kruger)投影，CGCS 2000國家大地坐標系，1985國家高程基準，D_China_2 000基準面。

2 土地利用轉移矩陣

2.1 土地利用轉移矩陣定義

土地利用轉移矩陣，就是根據(jù)同一地區(qū)不同時相的土地覆蓋現(xiàn)狀的變化關系，求得一個二維矩陣。通過對得到的轉移矩陣進行分析，能夠得到2個時相，不同的地類之間相互轉化的情況，它描述了不同的土地利用的類型在不同年份發(fā)生變化的土地類別以及發(fā)生變化的位置和變化面積。不僅能夠反映上述靜態(tài)的固定區(qū)域固定時間的各個地類面積數(shù)據(jù)，還能夠反映更加豐富的初期各個地類的面積轉出以及末期各個地類面積的轉入情況。從面積入手，反應區(qū)域土地利用變化。面積變化首先反映在不同土地利用類型的總量變化上，通過分析土地利用類型的總量變化，可了解土地利用變化總的變化趨勢和土地利用結構的變化。

轉移矩陣的應用主要表現(xiàn)在：采用GIS技術對兩期土地利用圖進行疊置分析從而獲得土地利用轉移矩陣，研究兩期之間各土地利用類型的轉移變化情況，了解研究期初各類型土地的流失去向以及研究期末各類型土地的來源與構成，通常結合土地利用數(shù)量變化、程度變化、空間變化等指數(shù)進行分析[6]。具體操作可描述為:在ArcMap中同時打開兩期土地利用現(xiàn)狀數(shù)據(jù)，根據(jù)土地利用類型字段對現(xiàn)有圖進行融合，利用疊加分析模塊下的交叉分析功能，對融合后的兩期土地利用數(shù)據(jù)進行疊加交叉操作，之后對獲得的交集圖層屬性表賦值并計算其面積，導出dbf屬性表，并在Excel中打開，執(zhí)行數(shù)據(jù)透視表和數(shù)據(jù)透視圖操作，即可生成土地利用轉移矩陣。

根據(jù)以上描述，運用ArcGIS的空間分析功能對準格爾旗2014，2016年2期土地利用現(xiàn)狀數(shù)據(jù)分析，獲取了準格爾旗兩期的土地利用轉移矩陣(表1)。

2.2 土地利用轉移矩陣動態(tài)變化指標

土地利用轉移矩陣動態(tài)變化指標主要包括土地利用類型轉換概率(P)和土地利用類型動態(tài)度(D)。

在獲取的土地利用轉移矩陣的基礎上，為了更好地分析土地利用動態(tài)變化的主導過程，建立土地利用轉移模型，其模型公式為：

(1)

式中：Pij——2014—2016年土地利用類型i轉換為土地利用類型j的轉換概率；Sij——土地利用類型i轉變?yōu)橥恋乩妙愋蚸的面積(km2)。根據(jù)式(1)中計算，獲得土地利用轉換概率矩陣(表2)。

土地利用動態(tài)度可以定量描述區(qū)域土地利用變化的速度，它對比較土地利用變化的區(qū)域差異和預測未來土地利用變化趨勢具有積極的作用。單一土地利用類型動態(tài)度表達的是研究區(qū)域一定時段范圍內(nèi)某土地利用類型的數(shù)量變化情況，其表達式為：

D=(Ab-Aa)/Aa×1/T×100%

(2)

式中：D——研究時段某一土地利用類型動態(tài)度；Aa——研究時期初某一土地利用類型的數(shù)量；Ab——研究期末某一土地利用類型的數(shù)量；T——研究時段長，當T的時段設定為年時，D值就是該研究區(qū)某種土地利用類型的年變化率。

表1 準格爾旗土地利用轉移矩陣

表2 準格爾旗土地利用轉換概率矩陣 %

根據(jù)公式(1)和(2)，獲得準格爾旗某一土地利用類型的年變化率如表4所示。

表3 準格爾旗土地利用年變化率

3 結果與分析

3.1 土地利用變化的主導過程

根據(jù)表1—3可以看出，準格爾旗2014—2016年3 a內(nèi)的各土地利用類型面積總體變化不顯著。由準格爾旗土地利用年變化率可以看出，該區(qū)內(nèi)園地、其他土地、居民點及工礦交通用地和耕地的年變化率為正，表明這4種土地利用類型的面積呈增加趨勢；林地、草地和水域及水利設施用地年變化率為負，表明這3種土地利用類型的面積呈減少趨勢。準格爾旗3 a的土地利用變化的主導主要集中于其他土地和耕地面積的增加、農(nóng)村居民點的擴張以及工礦交通用地的擴建，林地和草地大量減少，轉變?yōu)槠渌恋乩妙愋汀?/p>

3.2 土地利用變化率分析

林地和草地年變化率分別為-0.15%和-0.17%，林地面積由2014年3 456.94 km2減少到2016年3 440.92 km2，減少了16.02 km2；水域及水利設施用地年轉換率為-0.31%，面積由2014年的269.49 km2，減少到2016年266.6 km2。盡管林地和草地的自身轉換率比較高，分別為98.34%的林地和97.79%的草地，但自身轉換之間可能導致該區(qū)域的林地和草地的覆蓋度發(fā)生變化，結合減少的水域及水利設施面積可以推斷，該區(qū)域的植被覆蓋度出現(xiàn)降低趨勢。

居民點及工礦交通用地年變化率為1.31%，面積由2014年的474.68 km2增加到2016年的493.31 km2，凈增加18.63 km2，除92.53%來源于自身外，4.71%的草地和1.26%的林地轉變?yōu)榫用顸c及工礦交通用地。調(diào)查研究發(fā)現(xiàn)，導致居民點及工礦交通用地面積增加的主要原因為準旗境內(nèi)擁有豐富礦業(yè)資源，新建工礦用地、擴建工礦用地較多，破壞了林地和草地面積；除此之外，正在新修的呼準鐵路和大路鎮(zhèn)周圍新建的公路也是導致林地草地面積減少，居民點及工礦交通用地增加的原因。

其他土地面積年轉化率2.31%，增加了10.53 km2，除86.7%為原有的其他土地外，其余主要由草地、居民點及工礦交通用地和林地轉移而來，轉移比例分別占3.65%，3.38%和4.48%。其他土地面積增加主要原因為人為擾動因素的影響，在建設前期場地平整階段的土地面積以及工礦用地棄土棄渣場植被恢復前表土覆蓋增加較多，在影像上顯示為裸地。

3.3 土地利用變化綜合分析與水土流失變化響應

準格爾旗2014和2016年不同土壤侵蝕強度面積如表4所示。該區(qū)域水土流失面積為除微度侵蝕以外的其他侵蝕的面積之和，可以看出準格爾水土流失面積增加了70.80 km2，結合土地利用變化可以看出，增加的水土流失面積與林地、草地和水域及水利設施用地減少有密切關系。除直接減少的林地和草地的面積導致了水土流失面積增加外，水域面積的減少以及98.34%的林地和97.79%的草地自身轉換引起的林地和草地覆蓋度等級的變化，也是導致該區(qū)域水土流失面積增加得另一個原因。

表4 準格爾旗土壤侵蝕統(tǒng)計

4 結 論

(1) 2014—2016年，準格爾旗各土地利用類型均發(fā)生了變化，其中園地、其他土地、居民點及工礦交通用地和耕地面積表現(xiàn)出增加趨勢；林地、草地和水域及水利設施用地面積表現(xiàn)出呈減少趨勢。林地和草地面積的減少與李忠峰等[7]研究的研究結果不一致，導致結果不一致的原因主要是2014—2016年，該區(qū)礦區(qū)的開采及鐵路的建設等人為擾動因素增加導致了該區(qū)的林地和草地面積減少。通過土地利用轉移矩陣的變化可以看出減少的林地和草地主要向居民點及工礦交通用地和其他土地轉移，轉移概率表現(xiàn)為：林地轉移工礦交通用地和其他土地的分別為0.79%和0.42%，是林地轉移耕地的4.6倍和2.5倍，草地轉移工礦交通用地和其他土地分別為0.68%和0.34%，是草地轉移耕地的4倍和2倍。

(2) 2014—2016年該區(qū)水土流失面積增加了70.80 km2。水土流失面積的增加主要與林地和草地面積減少以及林地和草地植被覆蓋度的變化有關。陳桂榮等[8]研究了2015年準格爾旗水土流失動態(tài)變化，表明該區(qū)水土流失增加與呼準鐵路和其它公路的修建、大路鎮(zhèn)周圍施工帶來了人為擾動，從而植被面積減少，工礦、裸地面積增加，此外，準格爾旗降雨量逐年減少，導致植被覆蓋度整體降低，也是導致該區(qū)水土流失增加的原因之一。

[1] 鄭艷慧,靲海亮,賈露,等.基于馬爾柯夫模型平頂山市土地動態(tài)變化預測[J].北京測繪,2011(4):1-4.

[2] 張偉,彭淑貞,程鵬,等.泰安市土地利用動態(tài)變化研究[J].國土與自然資源研究,2011(4):17-19.

[3] 岳東霞,杜軍,劉俊艷,等.基于RS和轉移矩陣的涇河流域生態(tài)承載力時空動態(tài)研究[J].生態(tài)學報,2011,31(9):2550-2558.

[4] 劉瑞,朱道林.基于轉移矩陣的土地利用變化信息挖掘方法探討[J].資源科學,2010,32(8):1544-1550.

[5] 陸平.基于轉移矩陣的城市土地利用變化分析[J].北京測繪,2011(1):13-16.

[6] 喬偉峰,盛業(yè)華,方斌,等.基于轉移矩陣的高度城市化區(qū)域土地利用演變信息挖掘:以江蘇省蘇州市為例[J].地理研究,2013,32(8):1497-1507.

[7] 李忠鋒,王一謀,馬安青.準格爾旗土地利用變化與影響因素分析[J].水土保持通報,2003,23(3):42-44.

[8] 陳桂榮,杜婧.基于遙感監(jiān)測分析準格爾旗水土流失動態(tài)變化[J].內(nèi)蒙古水利,2016(9):41-42.