2000-2018 年甘肅省土地利用變化分析

郭麗麗

（甘肅省自然資源規劃研究院，甘肅 蘭州 730000）

土地利用/土地覆被變化(LUCC)研究一直以來都是地球科學研究的熱點方向[1]，土地利用變化研究不僅反映出土地資源的自然屬性，也揭示了地區的經濟發展狀況。土地利用變化關會影響到氣候變化、生物多樣性變化以及人與自然環境的相互關系[2]。土地利用變化研究通過分析土地利用類型的時空變化規律、土壤屬性的變化、土地功能的轉變等，對生態環境和生物多樣性的改善有重要意義。國內外學者針對土地利用變化主題做了大量研究，主要體現在土地利用類型面積的變化[3-7]、土地利用類型的重心變化規律[8-12]、土地利用變化的景觀格局變化[13-18]等方面。

文章基于已有的土地利用數據產品，采用土地利用轉移矩陣、土地利用類型重心遷移、景觀格局指數等方法分析甘肅省2000-2010 和2010-2018年兩個時期內的土地利用格局變化，揭示了甘肅省近20 年的土地利用類型的變化規律，為今后的土地利用規劃、土地資源管理以及土地利用監測等方面提供數據支持。

1 研究區與數據

1.1 研究區域

甘肅省位于中國西部地區，地處黃河中上游，地域遼闊，位于北緯32°11′～42°57′，東經92°13′～108°46′之間，總面積42.58 萬km2（如圖1 所示）。甘肅省下轄12 個地級市、2 個自治州，17 個市轄區、5 個縣級市、57 個縣、7 個自治縣。甘肅地形呈狹長狀，地貌復雜多樣，山地多，平地少，全省山地和丘陵占總土地面積的78.2%。全省土地利用率為56.93%，尚未利用的土地有28681.4 萬畝，占全省總土地面積的42.05%。

圖1 研究區位置圖

1.2 數據

土地利用數據來源于資源環境數據云平臺（http://www.resdc.cn/）2000 年、2010 年和2018 年的中國土地利用現狀遙感監測數據，分辨率為1km，基于Landsat 遙感影像，通過人工目視解譯生成土地利用數據，土地利用類型包括耕地、林地、草地、水域、居民地和未利用土地6 個一級類型以及25個二級類型，實驗所用數據為一級土地利用類型。所需矢量數據來源于中國自然資源部標準地圖服務中心（http://bzdt.ch.mnr.gov.cn/）。

2 研究方法

2.1 土地利用變化轉移矩陣

土地利用變化轉移矩陣是土地利用變化的類型轉移面積按矩陣的形式展示，是表示在某一段時期內，以定量地表明不同土地利用類型之間的轉化情況，反映出土地利用的用地結構與變化方向，以及各地類之間的相互轉化關系。

2.2 重心遷移

土地利用類型的重心遷移是反映出不同類型的土地利用類型的空間格局變化，能夠直觀反映出不同土地利用類型在某一時期的重心變化的過程，重心計算公式如下：

其中，Xi和Yi表示第i 個土地利用類型的經度和緯度坐標，LTij表示第i 類土地利用類型中的第j個斑塊的土地利用類型面積，Xij和Yij表示第i 類土地利用類型中的第j 個斑塊的經度和緯度坐標。

2.3 景觀格局指數

景觀格局指數主要包括兩方面，一是指景觀的空間結構特征以及人為或自然的干擾在內的各種生態過程在不同尺度上作用的結果，另外也指自然的和人為的多種因素相互作用所產生的一定區域生態環境體系的綜合反映。描述景觀格局的指數比較多[19-24]，本實驗主要用到：斑塊密度（PD）、景觀破碎化指數（SPLIT）、蔓延度指數（CONTAG）、香農均度指數（SHEI）、香農多樣性指數（SHDI）、散布與并列指數（IJI）和斑塊聚合度（AI）。指數表示的意義如下：

其中，（1）PD 表示單位面積上的斑塊數量，也是景觀破碎化的指標，PD 值越大，地塊破碎化越嚴重；（2）SPLIT 呈現出景觀被分割的破碎程度，反映出土地景觀格局的復雜性，也反映出人類對生態環境系統的影響，SPLIT 值越大，說明人類對生態系統的影響越大；（3）CONTAG 反映了不同類型的斑塊的團聚程度和蔓延趨勢，CONTAG 值越大，說明斑塊之間有很好的連通性；（4）SHEI 反映了景觀中的不同斑塊類型的分配均勻度和優勢度，SHEI 值越大，表明斑塊類型分布均勻，多樣性越大；（5）SHDI 反映出景觀中不同斑塊類型的復雜性和變異性，SHDI 值越大，表明斑塊類型增加斑塊類型呈現均衡化趨勢分布；（6）IJI 表征不同斑塊類型之間的總體散布與并列狀況，IJI 值越大，表明各斑塊類型之間的相鄰的邊長越接近，各斑塊之間的比鄰概率均等；（7）AI是度量景觀中斑塊聚集程度的指數，AI 值越大，說明斑塊高度聚集。

3 研究結果與分析

3.1 土地利用類型轉移矩陣分析

通過2000、2010 和2018 年甘肅省的土地利用數據（如圖2 所示），建立2000-2010 年和2010-2018 年的土地利用變化轉移矩陣，分析10 年間甘肅省的土地利用變化情況。2000-2010 年的土地利用變化見表1，6 類土地利用類型轉化面積最大為耕地和未利用地，其次是草地，水體和建設用地轉化面積較少。其中，耕地轉化面積為1284km2，主要轉化為草地；未利用的轉化面積為1224km2，主要轉化為耕地和草地；草地主要轉化為耕地和林地。通過統計土地類型變化比例，發現耕地和水體的轉化比例比較高，分別為1.964%和1.792%。

圖2 2000-2018 年的土地利用分布圖

2010-2018 年的土地利用變化見表2，和2000-2010 年之間的轉化比較類似，也是耕地和未利用地轉化面積最大，其次是草地，水體和建設用地轉化面積較少，但是轉移的土地利用類型不同。其中，耕地轉化面積為787km2，主要轉化為建設用地；未利用的轉化面積為742km2，主要轉化為耕地和建設用地；草地主要轉化為耕地、建設用地和未利用地，耕地和水體的轉化比例較高，分別為1.204%和1.391%。

表1 2000-2010 年之間的土地利用變化轉移矩陣（單位：km2）

表2 2010-2018 年之間的土地利用變化轉移矩陣（單位：km2）

通過對比2 個時期的土地利用變化矩陣，發現2000-2010 年之間耕地轉化為草地，符合“退耕還林、還草”的政策，同時為了確保耕地的基數，部分未利用地也被轉化為耕地；而2010-2018 年期間是經濟快速發展的階段，建設用地的面積呈現增長，促進城市化的快速發展。

3.2 土地利用類型重心變化分析

土地利用類型的重心變化反映出不同類型土地利用的空間格局變化，圖3 是2000-2018 年間甘肅省的6 類土地利用類型的重心遷移變化結果。從圖可知：2010-2018 年的各類土地類型重心遷移量均高于2000-2010 時期，說明2010-2018 年之間的土地利用格局變化較大。在整個時期內，耕地的重心依次向西北方向遷移；林地的重心開始向東北方向遷移，之后向東南方向遷移；草地的重心開始向南方向遷移，之后向東南方向遷移；水體的重心先開始向東南方向遷移，之后向西北方向遷移；建設用地的重心開始向東方向遷移，之后向西北方向遷移；未利用土地的重心開始向西方向遷移，之后向西北方向遷移。

圖3 反映出整個甘肅省的總體土地利用類型重心遷移情況，為了進一步分析局部地區的土地利用類型的遷移情況，實驗展示出以市級為單元的土地利用類型的重心分布格局（如圖4 所示）。由圖4可知：2000-2018 年之間的土地利用類型重心變化較大的地區為張掖市、慶陽市和平涼市地區，其他地區的土地利用類型的重心變化不大。其中，張掖市在2010-2018 年間的草地的重心向西方向遷移；平涼市在2000-2010 年間的耕地重心向東北方向遷移；慶陽市在2010-2018 年間的耕地重心向東南方向遷移。

圖3 2000-2018 年間的土地利用類型的重心遷移結果

圖4 2000-2018 年間市級單元的土地利用類型重心空間格局

3.3 土地利用類型景觀格局變化結果

利用Fragstats4.2 軟件分別計算2000、2010 和2018 年的土地利用類型的7 個景觀格局指數，結果見表3。結果表明：2018 年的PD 值較高，表明在2010-2018 年間的斑塊破碎化較為嚴重；2018 年的CONTAG 值較低，表明在2010-2018 年間的斑塊破碎，連通性較低；2018 年的IJI 值較高，表明各斑塊之間的比鄰概率較高；2018 年的SPLIT 值有所下降，表明與2000-2010 年相比，2010-2018 年之間的人為影響略小；兩個時期內的SPLIT、SHDI 和AI 指數基本一致，表明2000-2018 年之間的斑塊變異性很小，聚合度變化比較小。

表3 土地利用景觀格局指數

4 小結

本研究分析了甘肅省近20 年間的土地利用變化規律，建立土地利用變化轉移矩陣，分析不同土地類型的轉化量和轉化率，分析2 個時期（2000-2010 年和2010-2018 年）內各類土地利用類型的重心遷移規律，并分區域統計每個時期的各類土地利用類型的重心位置，分析重心坐標變化趨勢。結合景觀格局指數分析不同年份的土地利用空間結構，對比分析隨時間的變化規律。為甘肅省今后的土地利用規劃研究以及土地資源管理宏觀調控方面提供技術支持和數據基礎。