基于CA—Markov模型的巴里坤湖流域土地利用動態變化分析與預測

呂娜+李政海+鮑雅靜+李卓玲+張靖+劉麗+劉翀+徐媛+譚嫣辭+宮曉雙

摘 要：隨著經濟和農業的快速發展，人類活動對流域生態環境的影響日益顯著。在巴里坤哈薩克自治縣巴里坤湖流域環境變化已經得到了當地政府的重視。本文利用CA-Markov模型對2006—2016年該流域的土地利用動態變化進行分析，并模擬預測2020年的土地利用變化。結論如下：（1）2006—2016年間，該地區林地、水域面積小幅增加；草地、未利用土地面積減少；居民點、工礦用地、鹽堿地、耕地面積增加較明顯；（2）模擬結果表明，巴里坤湖及周邊生態系統將呈現出持續惡化的趨勢。到2020年，草地、湖泊和冰川、沼澤地面積減少，而灘地、水庫塘坑、鹽堿地耕地面積增加較多，將導致入湖及生態用水減少，會導致湖濱濕地及周邊生態系統退化、消失，使得其生態屏障的功能也隨之下降。

關鍵詞：預測；土地利用變化；元胞自動機；巴里坤流域

中圖分類號：F301.2 文獻標識碼：A DOI 編碼：10.3969/j.issn.1006-6500.2017.07.017

Analysis and Prediction of Land Use Dynamic Change in Barkol Lake Watershed Based on CA-Markov Model

LYU Na， LI Zhenghai， BAO Yajing， LI Zhuoling， ZHANG Jing， LIU Li， LIU Chong， XU Yuan， TAN Yanci， GONG Xiaoshuang

（College of Environment and Resources， Dalian Minzu University， Dalian， Liaoning 116600， China）

Abstract：With the rapid development of economy and agriculture， the impact of human activities on the ecological environment of the basin is becoming increasingly significant. The environmental changes in the Barkol Lake Basin in the Barry Kun Kazakh Autonomous County have got the attention of the local government. In this paper， the CA-Markov model was used to analyze the dynamic changes of land use in this basin from 2006 to 2016， and simulate the forecast of land use change in 2020. Conclusion as below：（1） from 2006 to 2016， the area of forest land， water area increased slightly； grassland， unused land area decreased； residential， industrial land， saline land and cultivated area increased more obvious； （2） the simulation results show that the Barkol Lake and the surrounding ecosystem will show a trend of continued deterioration. By 2020， grassland， lake and glacier， marsh area decreased， and beach， reservoir pond， saline land to increase the area， will lead to the lake and ecological water reduction， will lead to a large number of lake wetlands disappear， making its ecological barrier function also will decline.

Key words： prediction； land use change； cellular automata； Barkol basin

土地是人類生存的基礎，其動態變化是自然、社會、經濟生態復合系統演變的直接表現[1]。自20世紀以來，由于人口急劇增長，可利用的土地資源相對越來越少，土地利用問題逐漸引起世界各國的重視；同時隨著經濟的迅速增長，土地利用同生態環境之間的問題日益突出，并得到了人們的關注；1995年“國際地圈與生物圈計劃”（IGBP）和“全球環境變化人文計劃”（IHDP）聯合提出了“土地利用和土地覆蓋變化”（LUC）研究計劃，隨后土地利用變化成為全球環境變化研究的前沿和熱點課題[2]。土地利用是人類根據其生活和社會發展的需要，依據土地的性質特點與其經濟區位而對土地進行利用、開發與改造的社會經濟行為[3]，這一系列行為對生態環境與生態服務功能產生重要影響，并以累積的形式影響全球環境變化；目前人類面臨的許多環境問題，諸如從局地到全球尺度的氣候變化、生物多樣性的喪失與生態系統服務功能下降和土地退化等都與土地利用/覆蓋變化有關[4]；因而對土地利用結構變化特征及其規律的研究，是開展土地資源優化配置和合理調控的理論基礎和關鍵內容，對于區域產業結構轉型升級、土地資源可持續利用與生態環境保護具有重要的理論和指導意義[5]。

巴里坤湖位于新疆哈密市東天山區域，是巴里坤縣主要的湖泊，面積由歷史最大的800多km2萎縮到112 km2，湖水水質下降、礦化度高達280 g·L-1，流域內的植被退化嚴重，草場沙化、鹽堿化嚴重，流域生態服務功能大大減弱[6]。對巴里坤縣的自然生態系統和當地農牧民的生產與生活產生不利的影響。因此，通過對巴里坤2006、2010、2016年三期土地利用進行分析與預測，得出該區土地利用變化速度、變化程度及未來變化趨勢，為分析巴里坤湖流域土地資源開發利用強度評估流域生態系統安全狀況，為規避不合理的資源開發方式和經濟增長方式，實現該流域的社會經濟可持續發展提供依據。

2 研究區概況與研究方法

2.1 研究區概況

巴里坤湖流域位于新疆哈密市巴里坤哈薩克自治縣的巴里坤盆地。坐標范圍：東經92°05′35″～93°48′01″，北緯43°22′37″～44°14′33″。位于東天山的兩支山脈之間，流域面積為5.91×105 hm2。屬溫帶大陸性氣候，暖季涼爽、冷季嚴寒、光照不足，無霜期短，降水偏少，蒸發量大，年均降水量202 mm，蒸發量1 638 mm。四季不明顯，只有冷暖兩季，光照充足，無霜期98～104 d，年均氣溫1 ℃，≥10 ℃年積溫為1 500～1 900 ℃。

2.2 數據來源與數據處理

使用2006年、2010年陸地衛星（Landsat 5）TM遙感圖像，以及2016年Landsat 8 OLI遙感圖像，空間分辨率為30 m，云量均<2%。使用2000年已校正的TM遙感圖像對數據進行幾何校正，校正方法為二次多項式，重采樣方法為最近鄰法，均方根誤差（RMS Error）小于2。使用ArcGIS 10.3軟件對遙感圖像進行目視解譯，清除碎屑多邊形（最小圖斑為4 hm2），最后獲得2006年、2010年和2016年三期土地利用圖。

依據已有的土地利用分類原則和分類系統和參考一些學者工作的基礎上，并結合當地土地利用的實際情況，最終確定6個1級類和20個2級類反映巴里坤流域土地利用。土地利用類型為：林地（灌木林地、疏林地、其他林地）、草地（高覆蓋度草地、中覆蓋度草地、低覆蓋度草地）、水域（湖泊、水庫、塘坑、冰川積雪、灘地）、城鄉居民點和工礦用地（城鎮用地、農村居民點）、未利用土地（沙地、戈壁、鹽堿地、沼澤地、裸土地、裸巖石礫地）、耕地（耕地）。

2.3 CA-Markov模型預測

土地利用預測模型使用IDRISI TerrSet軟件的CA-Markov模型，該模型集成元胞自動機模型（Cellular Automaton，CA）和馬爾可夫鏈模型（Markov），極大地提高了傳統Markov模型的模擬精度，并通過設定約束條件的形式，又擁有約束性細胞自動機（Constrained CA）的優點[7]。模型建立步驟包括：首先，采用該軟件的Markov模塊，計算2010—2016年巴里坤湖泊流域土地利用面積轉移矩陣；其次，以各土地利用類型為因變量（1為是，0為否），分別選取DEM，到城鎮、公路和河流湖泊的距離為自變量，使用Logistic模塊，建立各土地利用類型的Logistic回歸方程，作為類型變化的約束條件；最后，在CA-Markov模型中模擬2020年的土地利用變化。

3 結果與分析

3.1 土地利用變化分析

對巴里坤湖流域土地利用一級類面積[8-9]（2006—2016年均值）按從大到小排列：草地>未利用土地>耕地>林地>水域>居民點、工礦用地。其中，草地為該流域的最大土地利用類型，面積約占該流域總面積的62.9%；未利用土地所占比例次之，面積占20.8%；其后是耕地，面積占7.8%；林地、水域和居民點、工礦用地在該流域中比例較小，不到10%，分別為4.2%，3.8%及0.4%。

2006—2016年該地區的林地變化不大，僅從24 760.5 hm2增加到24 896.2 hm2，主要以灌木林增加為主；草地在10年間減少了4 324 hm2，其中以中覆蓋度草地減少為主（1 161 hm2）；水域面積有所增加（1 208 hm2），湖泊、水庫塘坑及冰川積雪均表現相對各自面積的不小幅度的增加（湖泊：16.8%，水庫塘坑：50.7%，冰川積雪：79.2%，波動幅度較大），但灘地面積減少13.4%；居民點、工礦用地增加近346 hm2，增加比例為16.4%；未利用土地面積總體減少了1 190 hm2，其中沼澤地減少147 hm2，裸巖石礫地減少1 160 hm2，戈壁少量減少，但鹽堿地面積增加了431 hm2；耕地面積增加較大，為3 824 hm2。各土地利用二級類變化如圖1所示。

3.2 土地利用變化預測

為了驗證該模型的有效性，采用模擬相近年份的土地利用變化，來驗證該模型的有效性，通過建立2006—2010年的土地利用面積轉移矩陣，向2016年進行模擬，然后采用交叉驗證的方式進行精度對比。本次檢驗的卡帕系數（Kappa）為0.986≥0.75，表明模型模擬效果較好。通過CA-Markov模型模擬，獲得了2020年巴里坤湖流域的土地利用變化情況（一級類和二級類變化見圖2）。到2020年林地有小幅增加（132 hm2），主要以其他林地增加為主；草地減少1 139 hm2，中度、低度覆蓋度的草地均有小幅度的減少；水域總體面積將減少171 hm2，其中湖泊減少比例較大，約13.5%（1 445 hm2），冰川積雪減少比例為10.2%（190 hm2），而水庫塘坑面積小幅增加，灘地面積增幅較大為13.9%（1 445 hm2）；城鄉居民點和工礦用地較2016年略有增加：154 hm2；未利用土地面積減少722 hm2，總體減少幅度不大；耕地面積增加較多為1 747 hm2。由預測結果可以看出，該地區水域面積、草地面積與未利用土地面積均將減少，以草地面積減少量最大，而耕地面積增加幅度較大，湖泊面積轉為灘地面積量將增加。

4 結 論

通過對巴里坤2006、2010、2016年三期土地利用進行分析與預測，得出該區土地利用變化速度、變化程度及未來變化趨勢。結論如下：（1）在2006—2016年該地區的林地面積小幅增加，草地面積減少（以中覆蓋度草地減少為主），水域面積有所增加（湖泊、水庫塘坑及冰川積雪都表現出小幅增加），居民點、工礦用地面積增加，耕地面積增幅較大；（2）生態用水減少、生態屏障受損，巴里坤湖及周邊生態系統將呈現出持續惡化的趨勢。在土地利用模擬中，到2020年，中度、低度覆蓋度的草地面積均有小幅度的減少；湖泊和冰川面積呈減少趨勢，灘地面積增加，水庫塘坑面積增加，沼澤地面積減少，鹽堿地面積增加；耕地面積增加較多，將導致農業灌溉用水量增加。這些土地利用變化將會導致進入巴里坤湖泊及濕地生態用水減少，會導致大量湖塘及湖濱濕地的消失，使得其凈化水質的功能也隨之下降。為此在未來的發展中，該地區行政管理者應采取有效的手段控制耕地面積、優化農業結構、推廣農業節水灌溉、節制放牧等措施，是改善湖區周邊生態環境、提高生態服務功能。

參考文獻：

[1]倪紹祥，譚少華.近年來我國土地利用/覆被變化研究的進展[C]//全國土地覆被變化及其環境效應學術研討會，2002.

[2]潘佩佩，楊桂山，蘇偉忠.土地利用變化對土地生產力的影響研究進展[J].地理科學進展，2012，31（5）：539-550.

[3]劉平輝，郝晉珉.土地利用分類系統的新模式——依據土地利用的產業結構而進行劃分的探討[J].中國土地科學，2003，17（1）：16-26.

[4]張靖，同麗嘎，董建軍，等.土地利用/覆蓋變化研究方法[J].內蒙古大學學報（自然科學版），2013，44（3）：329-336.

[5]余德貴，吳群.基于Logistic-Markov的土地利用結構變化多因素驅動預測模型研究與應用[J].水土保持通報，2017，37（1）：104-110.

[6]劉紀遠，匡文慧，張增祥，等.20世紀80年代末以來中國土地利用變化的基本特征與空間格局[J].地理學報，2014，69（1）：3-14.

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

[8]張亞麗，黃珺嫦，蔡運龍，等.基于規劃協調的鄉鎮土地利用統一分類研究[J].地域研究與開發，2011，30（5）：150-155.

[9]文銳，吳宇哲.《土地利用現狀分類》實施對中國土地管理影響之管見[J].資源科學，2010，32（4）：731-736.