胡悅琴，馬燕飛，張偉科

（邯鄲學(xué)院地理系，河北邯鄲 056005）

0 引言

土地利用/覆被變化研究始終是全球變化研究的熱點之一，其研究框架是在1990年由美國全球變化委員會首次提出的，并第一次對該類研究給予研究名稱、提出研究方向，IGBP（全球地圈與生物圈計劃組織）和IHDP（全球變化中的人文領(lǐng)域計劃）組織于1995年共同對LUCC（土地利用/覆被變化）的可行性進(jìn)行分析，制定核心項目計劃，將LUCC研究列入核心計劃[1]。持續(xù)10余年的LUCC計劃后，2005年又共同推出了全球土地計劃(GLP)，目前，GLP已成為土地系統(tǒng)變化研究的又一核心研究計劃[2]。隨著3S 技術(shù)的發(fā)展，土地利用/覆被變化的研究也取得了很多研究成果，通過對比分析不同時期的LUCC 信息，可以獲取土地利用變化的空間位置、數(shù)量和類型，研究土地利用類型變化的過程、特點和驅(qū)動因素等。如Brondizio 等[3]基于遙感數(shù)據(jù)，探究了巴西亞馬遜地區(qū)1971—2001年間熱帶雨林的時空變化規(guī)律及其區(qū)域分異的動因；Dang 等[4]基于系統(tǒng)動力學(xué)和地理信息系統(tǒng)(GIS)，并綜合考慮社會經(jīng)濟(jì)和生物物理因素，構(gòu)建了越南達(dá)喀爾省LUCC 驅(qū)動模型，結(jié)果顯示該模型有助于規(guī)劃者和決策者制定可持續(xù)的規(guī)劃和管理戰(zhàn)略。

隨著國際上對土地利用/覆被變化的研究，中國學(xué)者也緊密結(jié)合土地利用與全球變化進(jìn)行了研究，并結(jié)合3S 理論和技術(shù)，在土地利用的數(shù)量變化研究、驅(qū)動力研究、土地利用變化模擬研究方面取得了豐碩的研究成果。土地利用變化包括土地利用類型的面積變化、空間變化和質(zhì)量變化。目前，土地利用的數(shù)量變化的指標(biāo)有土地利用變化的幅度、土地利用變化的速度和土地利用變化的區(qū)域差異，通過各地類之間的轉(zhuǎn)化，得出土地利用的轉(zhuǎn)化矩陣。一些學(xué)者提出了一些方法和模式，包括土地利用動態(tài)度的模型、土地利用程度的模型和土地利用變化方向的模型、土地利用空間變化的模型等[5]。如劉紀(jì)遠(yuǎn)等[6]基于3S 技術(shù)綜合運用動態(tài)度、轉(zhuǎn)移矩陣等方法，探究了20世紀(jì)80年代至21世紀(jì)初期中國土地利用變化的時空特征；張海燕等[7]在3S的基礎(chǔ)上運用轉(zhuǎn)移矩陣、動態(tài)度的方法，定量分析了全國退牧還草工程區(qū)的土地利用變化情況等。

有關(guān)LUCC驅(qū)動力的研究在土地利用變化機制研究中一直占主導(dǎo)地位。驅(qū)動力一般分為自然驅(qū)動力和人為驅(qū)動力，其中自然驅(qū)動力相對較為穩(wěn)定，發(fā)揮著累積性效應(yīng)，人為驅(qū)動力則相較活躍[8]。所以，目前的研究更多集中在對土地利用變化人為驅(qū)動力的探討方面。

中國是轉(zhuǎn)型的發(fā)展中國家，人地關(guān)系高度緊張，環(huán)境問題突出[9]。因此，在流域尺度上，利用土地利用轉(zhuǎn)移矩陣、土地利用動態(tài)度模型等研究中國陸地區(qū)域LUCC 時空格局變化，并定性分析土地利用類型變化的驅(qū)動因素，這樣可以更好地了解其近年的現(xiàn)狀，以期為各流域的生態(tài)環(huán)境健康和資源的可持續(xù)發(fā)展提供科學(xué)依據(jù)，為政府制定相關(guān)管理政策提供理論依據(jù)[10]。

1 資料與方法

1.1 研究區(qū)概況

中國陸地區(qū)位于亞洲東部，太平洋西岸，背陸面海，海陸兼?zhèn)洌I(lǐng)土面積約960萬km2，地勢西高東低呈現(xiàn)三級階梯分布，大部分地區(qū)屬于溫帶、副熱帶季風(fēng)氣候，表層在季風(fēng)與青藏高原共同影響下形成了三大區(qū)，即東部季風(fēng)區(qū)、西北干旱區(qū)和青藏高寒區(qū)，早期的土地利用遵照這種規(guī)律。但隨著國內(nèi)社會經(jīng)濟(jì)的持續(xù)上漲和工業(yè)化、城鎮(zhèn)化進(jìn)程的加速，人地矛盾問題越來越尖銳，資源環(huán)境與生態(tài)問題也越來越嚴(yán)重。

1.2 數(shù)據(jù)來源與處理

本研究數(shù)據(jù)基礎(chǔ)來自于全球30 m 地表覆蓋遙感數(shù)據(jù)產(chǎn)品GlobeLand30 (http://www.globallandcover.com/GLC30Download/index.aspx)的2000 年和2010 年2 期的地表覆蓋數(shù)據(jù)，在流域尺度的研究上利用了全國流域矢量圖，參照了全國流域分區(qū)圖。利用ENVI軟件對下載的GlobeLand30 的2 期數(shù)據(jù)進(jìn)行拼接，用ArcGIS 根據(jù)研究區(qū)邊界進(jìn)行裁剪，得到2 期中國土地利用圖，統(tǒng)計各地類的類型并對統(tǒng)計結(jié)果進(jìn)行分析（如圖1所示）。全球30 m地表覆蓋遙感數(shù)據(jù)產(chǎn)品將土地分成了耕地、森林、草地、灌木地、濕地、水體、苔原、人造地表、裸地、冰川和永久積雪10 類土地利用類型[7]（具體見表1）。由于研究區(qū)域為中國陸地區(qū)域，研究海域的變化沒有意義，再者中國陸地區(qū)域沒有苔原類型，所以研究區(qū)有9 個土地利用類型（耕地、森林、草地、灌木地、濕地、水體、人造地表、裸地、冰川和永久積雪）。

1.3 研究方法

本研究采用定性分析與定量分析相結(jié)合的方法，并結(jié)合ENVI、ArcGIS、Excel 等軟件對數(shù)據(jù)進(jìn)行分析，利用土地利用轉(zhuǎn)移矩陣、土地利用動態(tài)度模型分析研究區(qū)2000—2010年土地利用時間與空間的變化特征，并定性的分析土地利用的驅(qū)動因素。用GIS軟件流域矢量圖和土地利用圖進(jìn)行疊加統(tǒng)計，用Excel軟件對.dbf數(shù)據(jù)處理得到各流域土地利用類型的變化量，用ArcGIS 軟件的工具，選擇相應(yīng)的數(shù)據(jù)區(qū)域，即可形成土地利用轉(zhuǎn)移矩陣。而單一土地利用動態(tài)度模型，又稱為地類變化率指數(shù)，該模型能夠反映某一特定時期不同地類的轉(zhuǎn)移速率和變動幅度[12-15]。計算見公式(1)。

表1 GlobeLand 30 m數(shù)據(jù)產(chǎn)品土地利用類型及定義[11]

式中，K為研究時段內(nèi)某一土地利用類型的動態(tài)度指數(shù)，La、Lb分別為研究期初和研究期末某土地利用類型的面積，T為研究時段。

2 結(jié)果與分析

2.1 中國陸地區(qū)土地利用時空格局變化

2.1.1 土地利用類型面積變化 表2顯示了中國陸地區(qū)域2 個時期各土地利用類型的面積和變化情況數(shù)據(jù)，分析中國陸地區(qū)域2000—2010 年土地利用類型變化可知，2000年研究區(qū)內(nèi)各土地利用類型面積從高到低依次為草地、森林、耕地、裸地、人造地表、水體、冰川和永久積雪、灌木地、濕地。2010 年與2000 年土地類型的面積由高到低的次序一樣，只是所占的比例發(fā)生了變化。從變化量（見圖2）來看，草地、耕地和濕地的面積呈下降趨勢，分別下降140413、20480、692 km2，其他土地利用類型的面積在同一時期呈增長趨勢，裸地和人造地表的變化達(dá)99645、26302 km2，其次為灌木地、冰川和永久積雪、森林、水體，變化面積分別為15241、8232、6879、5396 km2。

表2 2個時期土地利用類型面積和變化情況

2.1.2 土地利用類型空間變化 土地利用轉(zhuǎn)換矩陣能有效地說明研究區(qū)內(nèi)土地利用類型轉(zhuǎn)變的過程和流向[16]。如表3 所示，研究區(qū)內(nèi)主要是耕地、林地、草地和裸地之間的轉(zhuǎn)換。其中，草地的轉(zhuǎn)出量與轉(zhuǎn)入量都是最多的，并且，轉(zhuǎn)出量大于轉(zhuǎn)入量，這與近10年來草地面積大幅度減少的現(xiàn)象是一致的。從去向看，草地轉(zhuǎn)為裸地的面積最大，為206439 km2，其次是轉(zhuǎn)化為森林和耕地，分別是166135、103560 km2。而森林、耕地與裸地是草地增加的主要來源，面積轉(zhuǎn)化量分別是157559、105058、98646 km2。森林變化中，森林轉(zhuǎn)化為草地和耕地的量是最多的，但其轉(zhuǎn)入來源也主要是草地和耕地，分別為157559、140057 km2?？偟膩砜?，森林面積總量略有增加，這說明森林的變化主要是空間位置的變化。耕地變化中，轉(zhuǎn)出量與轉(zhuǎn)入量最多的都是森林和草地，且面積相當(dāng)，究其耕地變化的主要原因就是耕地轉(zhuǎn)化為人造地表的面積為68770 km2，但人造地表轉(zhuǎn)化為耕地的面積為47137 km2，這與研究期內(nèi)耕地面積減少的趨勢大致相同。綜合來看，研究區(qū)內(nèi)耕地、森林和草地主導(dǎo)了這一時期內(nèi)的交換變化。

2.1.3 土地利用類型動態(tài)度分析 土地利用的動態(tài)度研究是國內(nèi)外學(xué)者進(jìn)行LUCC 研究的重要手段，能客觀地表示土地利用變化的強度與深度，該模型能夠反映某一特定時期不同地類的轉(zhuǎn)移速率和變動幅度[17]。根據(jù)公式得出研究區(qū)土地利用動態(tài)度圖（如圖3 所示），2000—2010 年灌木地和人造地表表現(xiàn)出明顯的變化速率，變化速率都達(dá)到了1.80%，森林和耕地的變化速率較小，分別為0.03%、0.10%。另外，除了耕地、草地和濕地動態(tài)變化為負(fù)(K<0)外，其余地類動態(tài)變化均為正(K>0)?？偟膩碚f，土地利用類型的動態(tài)度變化緩慢。

2.2 流域尺度土地利用時空格局分析

中國陸地區(qū)域分成了10個流域區(qū)，分別為松花江區(qū)(A)、遼河區(qū)(B)、海河區(qū)(C)、黃河區(qū)(D)、淮河區(qū)(E)、長江區(qū)(F)、東南諸河(G)、珠江區(qū)(H)、西南諸河(J)、西北諸河(K)。各流域土地利用類型圖如圖4所示。

表3 2000、2010年土地利用類型轉(zhuǎn)換矩陣 km2

分析各流域區(qū)土地利用類型的面積及變化（見表4）可知，2010 年淮河區(qū)土地利用類型以耕地、人造地表和森林為主，西南諸河流域主要的3 種土地利用類型由高到低依次是草地、森林和裸地，西北諸河流域則以裸地、草地和耕地類型為主，其他流域都以草地、森林和耕地類型為主，但土地利用類型的組成具有顯著差異。此外，西北諸河流域的水體、冰川和永久積雪、灌木地的面積在所有流域中是最大的，松花江流域濕地類型的面積是最多的。

對比這10 個流域區(qū)的土地利用類型面積的變化可知，西北諸河和西南諸河流域的耕地類型呈增長趨勢，分別增加了2874、1561 km2，其他流域均呈減少的趨勢，但淮河流域的耕地類型減少的最多，減少了4905 km2；黃河區(qū)森林類型的增加與草地類型的減少較為明顯，分別變化了7124、9614 km2；西北諸河和西南諸河的草地減少得最為明顯，分別減少了111451、17632 km2，但同時裸地增長也明顯，分別增加了92562、7224 km2；各流域人造地表的面積都呈增長的趨勢，淮河區(qū)與長江區(qū)增長的最多，分別增長了6848、5497 km2。西南諸河冰川和永久積雪的面積增長的最多，增長了9707 km2，其他位于研究區(qū)內(nèi)東南部的流域區(qū)沒有此土地利用類型。

2.3 典型流域土地利用時空格局變化——以黃河流域為例

黃河是中國的母親河，黃河流域西起青藏高原，北鄰戈壁沙漠，流經(jīng)黃淮海平原，跨越三大階梯，流域面積遼闊，且黃河流域的各土地利用類型變化很明顯，在此把黃河流域單獨拿出來作為典型流域分析其土地利用時空格局變化[18]。土地利用類型如圖5所示。

表4 流域尺度土地利用類型面積與變化情況 km2

2.3.1 黃河流域土地利用類型面積變化 分析黃河流域2000年和2010年土地利用類型的面積和變化（見圖6）可以得出，研究區(qū)內(nèi)的土地類型以草地、耕地和森林為主，其次是裸地和人造地表，濕地、水體和灌木地類型的面積相對較小，冰川和永久積雪的面積最少。分析研究區(qū)內(nèi)各地類的變化可知，草地、耕地和裸地的面積在減少，分別減少了9551、3756、276 km2，其他類型的面積呈增長趨勢，其中草地的面積減少的最多，森林的面積增加的最多。

觀察組1例感染，發(fā)生率為1.89%，而對照組3例感染、2例壓瘡、1例肺炎，發(fā)生率為11.32%，組間比較差異有統(tǒng)計學(xué)意義（P<0.05）。同時，與對照組比較，觀察組的護(hù)理滿意度評分高，且切口愈合和術(shù)后血糖控制時間均較短，組間對比差異有統(tǒng)計學(xué)意義（P<0.05），見表1。

2.3.2 黃河流域土地利用類型空間變化 黃河流域土地利用類型轉(zhuǎn)移矩陣如表5 所示，轉(zhuǎn)出量和轉(zhuǎn)入量最大的是草地類型，其次是耕地和森林，其中草地、耕地和裸地的轉(zhuǎn)出量大于轉(zhuǎn)入量，這與黃河流域各土地利用類型的變化是一致的。草地主要向耕地、森林、裸地轉(zhuǎn)化，分別轉(zhuǎn)化了22996、13918、6587 km2，但草地的轉(zhuǎn)入來源也是耕地、森林和裸地，其次除耕地向草地轉(zhuǎn)移外，耕地向森林和人造地表的轉(zhuǎn)化量也很大?？偟膩砜?，草地、耕地和森林主導(dǎo)了黃河流域這一時期的交換變化。

2.3.3 黃河流域土地利用動態(tài)度分析 根據(jù)土地利用動態(tài)度公式，得出黃河流域土地利用動態(tài)度圖（見圖7），分析可知，灌木地與冰川和永久積雪的變化率最為顯著，分別是9.15%、10.91%，冰川和永久積雪類型因為基數(shù)小，主要是草地向其轉(zhuǎn)移的面積，故而轉(zhuǎn)移面積雖不多，但變化率顯著，耕地、草地和濕地的變化速率較小，都小于0.5%，變化幅度較小。除了耕地、草地和裸地的動態(tài)度為負(fù)數(shù)外，其余各地類動態(tài)度都為正?？傊魍恋乩妙愋偷淖兓什町愶@著。

3 驅(qū)動因素分析

表5 黃河流域土地利用轉(zhuǎn)移矩陣 km2

影響區(qū)域土地利用類型時空格局變化的驅(qū)動因素通?？煞譃樽匀灰蛩睾腿藶橐蛩豙19]。這2種因素都促進(jìn)了土地利用類型的變化，自然因素和人為因素在不同的空間和時間尺度上對土地等級的變化有不同的影響。

在長序列的時間維度內(nèi)，自然因素主導(dǎo)了土地使用的空間格局，在歷史條件下，決定和影響土地使用變化的背景條件，但人為因素是區(qū)域土地使用變化的內(nèi)在驅(qū)動因素，在很短的時間內(nèi)和在較小的區(qū)域范圍內(nèi)都是如此[20-21]。此外，在短期內(nèi)，自然因素相對穩(wěn)定，變化緩慢，對土地利用的影響相對有限，但對區(qū)域土地利用格局的影響是持續(xù)和穩(wěn)定的。

3.1 自然因素

中國幅員遼闊，地勢西高東低，呈階梯狀分布，地形復(fù)雜多樣，有高原、山嶺、平原、丘陵和盆地等地形類型，河流大都自西向東流。但山區(qū)面積廣大，平原較少，使耕地資源不足。地形影響土地利用類型及其變化方向，對土地利用空間格局的形成具有重要影響[22]。地形因素在土地利用類型變化中表現(xiàn)在地表物質(zhì)和能量再分配，并間接參與土壤和生態(tài)環(huán)境之間的物質(zhì)和能量交換。海拔、坡度、坡向和地形對土地利用時空格局變化的影響表現(xiàn)為海拔、坡度和坡向的差異，這導(dǎo)致太陽輻射的吸收和地面輻射的地理分化的差異[23]。另一方面，由于化石燃料的大量使用，使全球氣候變暖，還有旱澇、火災(zāi)、病蟲害、荒漠化等災(zāi)害，使耕地、草地類型的面積減少，裸地面積增加[24]。不同尺度的流域區(qū)自然因素的影響力也是不同的，但是在研究期內(nèi)，自然因素對土地利用的影響作用不太明顯。

3.2 人為因素

3.2.1 政策因素 土地利用是在一定的政策背景下進(jìn)行的，土地利用的形式和方式是由政策的約束作用決定的，政府的政策通常促進(jìn)或限制土地使用的變化，并在土地利用變化中發(fā)揮宏觀調(diào)控作用[25]。土地有償使用制度和產(chǎn)權(quán)制度對規(guī)范土地流轉(zhuǎn)具有重要意義。在具體區(qū)域，土地利用的類型也是在具體的經(jīng)濟(jì)水平和政策背景下制定的。例如，減免農(nóng)業(yè)稅調(diào)動了農(nóng)民生產(chǎn)積極性，促進(jìn)了耕地的保護(hù)，使耕地面積增長。退耕還林、還草、還濕等生態(tài)環(huán)保政策的實施，提高了土地利用率，森林和灌木地不斷擴(kuò)大。

由于政策因素，各流域內(nèi)土地利用類型時空格局的變化是不一樣的。比如黃河流域，由于黃土高原地區(qū)植被破壞和水土流失嚴(yán)重，生態(tài)環(huán)境脆弱，政府實施了三北防護(hù)林工程、天然林保護(hù)工程、退耕還林等政策，采取積極的土地利用和生態(tài)恢復(fù)政策，黃河流域耕地和裸地類型面積減少，森林和灌木地類型增加。

研究區(qū)內(nèi)，東部經(jīng)濟(jì)發(fā)展水平高，人造地表的面積增長的較多，而西北部經(jīng)濟(jì)發(fā)展水平比較低，耕地面積增長，但荒漠化比較嚴(yán)重，裸地面積增長，草地面積增加。

3.2.3 人口因素 人口因素與土地利用類型的變化有密切的關(guān)系，中國是人口大國，2000 年祖國大陸約12.7 億人，到了2010 年達(dá)到了13.7 億。人口不斷增長，土地是人類生產(chǎn)生活的場所，人們往往通過改變土地利用的類型和結(jié)構(gòu)來滿足生產(chǎn)生活的需要[28]，從而增加對土地的干預(yù)程度，不斷促進(jìn)土地利用方式的更替和演變[29-30]。因此人口密集的流域，耕地和草地的面積減少明顯，人造地表增長較多。

4 土地利用管理措施建議

4.1 土地利用存在的問題

研究區(qū)內(nèi)，草地類型的面積不斷減少，10 年內(nèi)面積減少了140413 km2，草地主要變?yōu)榱寺愕?，裸地面積10年內(nèi)增加了99645 km2，各子流域內(nèi)，草地類型的面積雖然有增長的，但幅度較小，主要還是以減少為主，這說明研究區(qū)內(nèi)草地荒漠化較嚴(yán)重，生態(tài)環(huán)境被破壞，造成生物棲息地也被破壞，生物多樣性減少。

人造地表的面積不斷增加，10年內(nèi)增加了26302 km2，究其來源，主要來自于耕地，10 年來耕地的面積減少了20480 km2，各子流域內(nèi)，大部分利于人造地表的面積是增加的，只有極少流域是減少的。人口增長，技術(shù)進(jìn)步，耕地的單位產(chǎn)量增加，人類對于非農(nóng)的需求增大，故而人造地表的面積增長。

4.2 土地利用管理措施建議

根據(jù)國家相關(guān)法律法規(guī)，保護(hù)草地、濕地等土地利用類型，嚴(yán)查非法擠占，保護(hù)生態(tài)環(huán)境，保護(hù)生物棲息地、魚類遷徙路線和重要濕地等生態(tài)空間[31]。加強森林資源保護(hù)和管理，落實林地利用條例，確保林地面積只增不減，擴(kuò)大森林面積和草地面積。在河流沿線修建森林和草地來保護(hù)河岸，并加強土壤采掘、挖掘砂土和采石等活動的管理?？刂迫嗽斓乇砻娣e用量，提高經(jīng)濟(jì)技術(shù)，進(jìn)而提高利用率，減少對其他土地利用類型的占用。各流域要因地制宜，制定適合本地生態(tài)發(fā)展的政策法規(guī)，合理利用土地。

5 結(jié)論

研究區(qū)內(nèi)土地利用類型主要以草地、森林和耕地為主，其中草地、耕地和濕地面積不斷減少，分別減少了140413、20480、692 km2，裸地和人造地表的面積增長的比較明顯，增長了99645、26302 km2；草地—裸地、草地—森林之間的轉(zhuǎn)化較為劇烈，草地轉(zhuǎn)化為裸地的面積為206439 km2，草地轉(zhuǎn)化為森林面積是166135 km2，而森林轉(zhuǎn)化成草地的面積是157559 km2；灌木地和人造地表的變化速率最大，達(dá)到了1.80%，但總體土地利用類型的動態(tài)度變化緩慢。

各子流域內(nèi)，土地利用類型時空格局變化不同，除淮河區(qū)、西南諸河與西北諸河外，大部分流域以草地、耕地和森林為主，西北諸河草地類型的減少與裸地類型的增加十分明顯，另外，淮河區(qū)與長江區(qū)更低的減少與人造地表的增加變化也很大。

黃河流域以草地、耕地和森林為主，草地的面積減少的最多，減少了9551 km2，主要轉(zhuǎn)化為了耕地、森林和裸地；森林和灌木地增長的最多，分別增加了7123、4007 km2，主要來源是草地；灌木地與冰川和永久積雪的變化率最為顯著，分別是9.15%和10.91%，耕地、草地和濕地的變化速率較小，都小于0.5%，各土地利用類型的變化率差異顯著。

土地利用時空格局發(fā)生變化，離不開自然因素和人為因素的影響，自然因素的影響作用是持續(xù)且穩(wěn)定的，不斷累積才能產(chǎn)生結(jié)果；人為因素的影響是活躍且多方面的，不斷地對土地進(jìn)行干預(yù)，不斷促進(jìn)土地利用方式的更替和演變。