翟榮飛++趙翠薇++楊莉

摘 要 以烏江北源2000年、2010年兩期土地利用和DEM數(shù)據(jù)為基礎(chǔ)，運(yùn)用GIS空間分析軟件，探究喀斯特流域土地利用變化的空間規(guī)律.結(jié)果表明： 1）2000—2010年期間，耕地大幅度減少（-2.29%），林地、草地和水域用地卻在增加（+1.53%、+0.21%、+0.47%）.2）土地利用類型面積變化差異顯著，在高程1 200～1 600 m、坡度15°～25°、與河流主干距離<10 km區(qū)域內(nèi)變化最為劇烈.3）主要地類轉(zhuǎn)化包括：以耕地轉(zhuǎn)向林地、草地、水域用地為主，多發(fā)生在高程1 200～1 600 m、坡度15°～25°的與河流主干距離<10 km區(qū)域內(nèi)；林地、草地轉(zhuǎn)向水域用地，多發(fā)生在高程<1 200 m、坡度8°～15°的與河流主干距離<10 km區(qū)域內(nèi)；林地向居住用地轉(zhuǎn)化，多發(fā)生在高程1 200～1 600 m、坡度8°～15°與河流主干距離<10 km區(qū)域內(nèi).4）土地利用變化程度與高程、坡度和河流主干距離有高度的相關(guān)性，且隨著高程、坡度的增加先增加再減小，隨著與河流主干距離增加明顯急劇較小的趨勢(shì).

關(guān)鍵詞 喀斯特流域；土地利用變化；高程；坡度；與河流主干距離

中圖分類號(hào) K903 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 10002537（2016）01000907

Natural Factor Analysis of Land Use Change in Karst Mountain Area

—North Source of Wujiang River

ZHAI Rongfei， ZHAO Cuiwei*， YANG Li

（Geographic and Environmental Sciences， Guizhou Normal University， Guiyang 550001， China）

Abstract

The natural environment as a basic factor that affects land use change， which directly has impact on the way and degree of land use. Especially on the karst mountain areas where the ecological environment is fragile， the uniqueness of the natural environment leads to the typical form of land use. On the basis of two issues DEM data （ 2000， 2010 year） on land use of north source in Wujiang river， GIS spatial analysis software was employed to explore the space law of land use change in karst basin. Our results showed that： （1） Cultivated land was reduced by a large margin （2.29%） during 2000—2010， while forest land， grass land and waters were increased by+1.53%，+0.21%，+0.47%， respectively. Significant differences were found in different land use types with the most dramatic change occurred in the area with the elevation of 1 200～1 600 m， gradient of 15°～25°， and river trunk distance <10 km. （2） Main transformation includes that cultivated land was mainly turned to forest land， grassland， water land， and they occurred in the area with the elevation of 1 200～1 600 m， slope of 15°～25° and river trunk distance < 10 km area. Forest land， grass land were converted to water land in the area where elevation<1 200 m， gradient of 8°～15° and river trunk distance <10 km. Forest land tumed to residential land elevation of 1 200～1 600 m， gradient of 8°～15° and river trunk distance<10 km area. （3） Degree of land use change a high correlation with elevation， slope， and rivers backbone distance. With an increased elevation and slope change increases first and then decreases. With an increase in distance of river trunk the change was found to decrease sharply.

Key words Karst basin； land use change； elevation； slope； distance to backbone

土地利用是人類活動(dòng)對(duì)自然環(huán)境的綜合表現(xiàn)，也是影響全球環(huán)境變化的重要組成部分.而土地利用類型的變化將對(duì)能量交換、水分循環(huán)和土壤侵蝕等陸地主要生態(tài)過程的結(jié)果和功能產(chǎn)生深遠(yuǎn)影響，同時(shí)也能直接反映一個(gè)區(qū)域生態(tài)環(huán)境和社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展程度 [14].

近年來喀斯特地區(qū)土地利用研究已成為國(guó)內(nèi)外廣泛關(guān)注的焦點(diǎn)，其研究主要包括：喀斯特土地利用結(jié)構(gòu)、類型特征[510]；土地利用變化的時(shí)空差異[1115]以及土地利用類型、結(jié)構(gòu)和變化受地形因子的影響程度[1620]等.生態(tài)環(huán)境脆弱地區(qū)，社會(huì)經(jīng)濟(jì)發(fā)展落后、人地矛盾突出、土地利用變化劇烈，從而反過來限制區(qū)域的可持續(xù)發(fā)展.地形因子作為自然環(huán)境的重要因素，特別是高程、坡度，對(duì)土地利用形式和時(shí)空分布產(chǎn)生直接影響；在流域區(qū)尺度上，水資源狀況影響人類對(duì)土地利用程度，尤其在自然條件獨(dú)特的喀斯特山地，其地形崎嶇、地表嚴(yán)重缺水，導(dǎo)致土壤能力低下、養(yǎng)分及水土流失狀況嚴(yán)重，因而地形、水資源狀況對(duì)喀斯特流域區(qū)的土地利用形態(tài)有著不可忽略的影響[5].

本文以烏江北源流域?yàn)槔?000年與2010年兩期土地利用矢量圖，運(yùn)用GIS空間分析軟件，在對(duì)高程、坡度及與河流主干距離分級(jí)的基礎(chǔ)上，分析烏江北源流域區(qū)2000—2010年土地利用變化在不同高程、坡度及與河流主干距離

的變化情況，探究喀斯特流域土地利用變化空間規(guī)律，為喀斯特流域區(qū)的土地資源利用、規(guī)劃以及政府決策提供依據(jù).

1 研究區(qū)概況

圖1 烏江北源區(qū)位圖

Fig.1 The geographic location map of the Wujiang north source

烏江，長(zhǎng)江上游右岸的最大支流，有南北兩源，其中南源三岔河，北源六沖河.六沖河又名烏江北源，地理位置介于東經(jīng)104°19′46.46″～106°7′18.86″、北緯26°31′～28.96″之間，流域區(qū)總面積為10 874 km2，小部分流經(jīng)云南省高雄鎮(zhèn)，大部分流經(jīng)貴州畢節(jié)、大方、赫章、納雍、織金、黔西以及威寧縣，貴州省境內(nèi)流域面積為9 934 km2（圖1）.

流域地勢(shì)呈西北東南走向，海拔913～3 034 m，碳酸鹽巖溶面積占流域區(qū)總面積的71.22%，典型喀斯特地貌特征，包括峰叢洼地、峰叢溝谷、峰叢谷地和溶丘洼地等.大部分地方屬亞熱帶濕潤(rùn)氣候，年均日照1 231 h、氣溫13.09 ℃、降雨量1 129.04 mm、無霜期266天，冬無嚴(yán)寒，夏無酷暑，四季分明，氣候宜人，適合于多種農(nóng)作物和動(dòng)物的生長(zhǎng)、繁衍.其經(jīng)濟(jì)情況表現(xiàn)為，2013年初步核算全區(qū)生產(chǎn)總值為1 041.93億元，按可比價(jià)格計(jì)算，比上年增長(zhǎng)15.1%，第一、二、三產(chǎn)業(yè)結(jié)構(gòu)比為18.8∶43.2∶38.0.按常住人口計(jì)算，人均生產(chǎn)總值為15 953元.

2 數(shù)據(jù)來源及方法

2.1 數(shù)據(jù)來源及處理

數(shù)據(jù)來源：（1）地球系統(tǒng)科學(xué)數(shù)據(jù)共享平臺(tái)獲取30×30 m分辨率的DEM；（2）遙感數(shù)據(jù)解譯貴州省2000年和2010年兩期土地利用.

數(shù)據(jù)處理：（1） 利用ArcGIS空間分析，對(duì)30×30 m DEM進(jìn)行水文分析，確定河網(wǎng)分布及研究區(qū)范圍（見圖1）；（2）利用ArcGIS中分析工具，提取研究區(qū)高程、坡度以及與河流主干距離分級(jí)專題圖和屬性數(shù)據(jù)庫；（3）疊加分析，獲取研究區(qū)2000～2010年土地利用變化圖和屬性數(shù)據(jù)；（4）疊加研究區(qū)土地利用變化圖與高程、坡度以及與河流主干距離分級(jí)圖，統(tǒng)計(jì)分析研究區(qū)土地利用變化在不同高程、坡度、與河流主干距離的分布規(guī)律.

2.2 研究方法

在GIS與RS技術(shù)支持下，利用DEM數(shù)據(jù)提取高程、坡度與坡向等地形因子信息，根據(jù)研究區(qū)區(qū)地形特點(diǎn)以及《貴州農(nóng)業(yè)地貌區(qū)劃》進(jìn)行高程、坡度和坡向與河流主干距離分級(jí).將研究區(qū)土地利用變化圖和高程、坡度和河流主干距離圖進(jìn)行疊加，統(tǒng)計(jì)分析在不同高程、坡度和河流主干距離帶上土地利用變化特征.

2.2.1 土地利用分類 參照全國(guó)土地利用現(xiàn)狀分類，結(jié)合研究目的，將土地利用類型劃分為林地、耕地、草地、水域、居住用地、工礦用地、交通用地、未利用地8類.

2.2.2 自然條件因子分級(jí) 研究區(qū)海拔介于913～3 034 m之間，屬云貴高原山地.高程決定局地的溫度和光照狀況，影響人類利用土地的難易程度，從而直接對(duì)土地利用的空間分布產(chǎn)生重要影響.又因研究區(qū)高程<1 200 m的占總面積的2.17%、>2 000 m占總面積的11.59%.因此，將研究區(qū)高程<1 200 m和>2 000 m界定為第一、四高程帶，1 200～2 000 m之間的400 m為間隔平均分為2個(gè)高程帶.

坡度是表征局部地表坡面在空間的傾斜程度的地形因子，坡度的大小直接影響地表物質(zhì)流動(dòng)和能量轉(zhuǎn)換的規(guī)模和強(qiáng)度，是制約土地利用的空間布局.按照《貴州農(nóng)業(yè)地貌區(qū)劃》坡度分級(jí)指標(biāo)，將研究區(qū)坡度分為微坡（0°～5°）、較緩坡（5°～8°）、緩坡（8°～15°）、較陡坡（15°～20°）、陡坡（20°～25°）和急陡坡（>25°）共6個(gè)等級(jí).

由于研究區(qū)為流域尺度下的喀斯特山地，地表水較為缺乏，水源因素對(duì)土地利用變化的影響不可忽略，所以將研究區(qū)與河流主干距離每隔10 km劃一級(jí)，將研究區(qū)與河流主干距離劃分為：近距離（<10 km）、較近距離（10～20 km）和遠(yuǎn)距離（>20 km）3個(gè)等級(jí)（圖2）.

其中高程1 600～2 000 m、坡度8°～15°、與河流主干距離<10 km所占面積比例最大（見表1）.

3 結(jié)果

3.1 地類面積變化總體特征

2000～2010年烏江北源流域區(qū)以林地、耕地和草地為主要地類，其中地類變化主要表現(xiàn)為：耕地在大幅度減少，林地、水域用地、草地、居住用地以及未利用地在增加，工礦用地和交通用地基本無明顯變化（見表2）.

由圖3可得，烏江北源2000～2010年以林地、耕地、草地為主要地類，林地在地勢(shì)較高的西北分布居多，耕地在地勢(shì)較低的東南分布較廣，而草地卻零星分布在整個(gè)研究區(qū).其中，耕地在明顯減少，林地、草地和水域用地卻在增加.這是由于研究區(qū)地勢(shì)陡峭，地形破碎，碳酸鹽巖發(fā)育較好，屬水土流失嚴(yán)重的西南喀斯特區(qū)，退耕還林政策的實(shí)施，將易造成水土流失的耕地向林草地轉(zhuǎn)型.研究區(qū)水能資源豐富，10 MW以上的河流有8條（包括干流），其中洪家渡電站是我國(guó)“西電東送”重要電站，從2000年11月8日起我國(guó)正式進(jìn)入西部大開發(fā)，這就是導(dǎo)致研究區(qū)水域用地在2000～2010年間面積急劇增加的主要原因.由于烏江北源區(qū)地形特點(diǎn)導(dǎo)致社會(huì)經(jīng)濟(jì)發(fā)展滯后，使之交通用地、工礦用地面積較少且基本無變化.

3.2 不同自然條件的土地利用變化分析

3.2.1 不同高程帶的土地利用變化 由表3可得，2000～2010年期間烏江北源流域區(qū)在1 200～1 600 m高程范圍內(nèi)，地類變化面積最大，占總轉(zhuǎn)化面積的43.58%，其次為1 600～2 000 m高程范圍內(nèi)占33.81%，轉(zhuǎn)化面積相對(duì)較小的為>2 000 m高程，只占總轉(zhuǎn)化面積的2.35%.表明研究區(qū)土地轉(zhuǎn)化主要在1 200～1 600 m高程范圍內(nèi)發(fā)生，隨高程的增加，土地轉(zhuǎn)化面積先增加再減小.總體來講，以耕地向林地、草地、水域地轉(zhuǎn)出為主，且在1 200～1 600 m高程范圍內(nèi)耕地向林地轉(zhuǎn)化面積最大，在研究期間交通用地、工礦用地與其他地類幾乎無轉(zhuǎn)化.主要地類轉(zhuǎn)化具體表現(xiàn)為：在<1 200 m高程范圍內(nèi)，耕地轉(zhuǎn)化為林地、草地、水域、未利用地；林地、草地轉(zhuǎn)化為水域用地.在1 200～1 600 m高程范圍內(nèi)，耕地向林地、草地、水域、居住用地、未利用地轉(zhuǎn)化；林地向水域、居住用地轉(zhuǎn)化；草地向水域用地轉(zhuǎn)化.在1 200～1 600 m高程范圍內(nèi)，耕地向林地、草地、居住用地、未利用地轉(zhuǎn)化，林地向水域用地轉(zhuǎn)化.在>2 000 m高程范圍內(nèi)以耕地向林地、草地轉(zhuǎn)化為主.

3.2.2 不同坡度帶的土地利用變化 由表4可得，2000～2010年期間烏江北源流域區(qū)在坡度15°～25°范圍內(nèi)，地類變化最為頻繁，占總轉(zhuǎn)化面積的31.68%，其次為坡度8°～15°范圍內(nèi)占28.07%，而坡度>35°范圍內(nèi)地類轉(zhuǎn)化最小，占總轉(zhuǎn)化的6.41%.表明研究區(qū)土地轉(zhuǎn)化主要在坡度為8°～25°范圍內(nèi)發(fā)生，隨坡度的增加，土地轉(zhuǎn)化面積表現(xiàn)為先增加再減小.總體表現(xiàn)為：以耕地向林地、草地、居住用地、水域、未利用地轉(zhuǎn)化為主，且在坡度為15°～25°范圍內(nèi)耕地向林地轉(zhuǎn)化面積最大，交通用地和工礦用地在研究期間幾乎無轉(zhuǎn)化.具體表現(xiàn)為：在坡度為0°～5°、5°～8°、8°～15°、15°～25°、25°～35°范圍內(nèi)，地類轉(zhuǎn)化方向相近，主要轉(zhuǎn)化地類包括：耕地轉(zhuǎn)化為林地、草地、水域、居住用地、未利用地；林地轉(zhuǎn)化為水域、居住用地；草地轉(zhuǎn)化為水域用地.在坡度>35°范圍內(nèi)主要以耕地轉(zhuǎn)化為林地、草地、水域、未利用地；林地、草地轉(zhuǎn)化為水域用地為主.

3.3 不同河流干道距離的土地利用變化

由表5可得，2000～2010年期間烏江北源流域區(qū)在與河流主干距離<10 km范圍內(nèi)地類變化最為頻繁，占總轉(zhuǎn)化面積的68.40%，其次為與河流主干距離10～20 km范圍內(nèi)占26.72%，轉(zhuǎn)化面積最少的為與河流主干距離>20 km范圍內(nèi)占4.88%.表明研究區(qū)土地轉(zhuǎn)化主要在與河流主干距離<10 km范圍內(nèi)發(fā)生，表現(xiàn)為隨著與河流主干距離的增加土地轉(zhuǎn)化面積急劇減小的趨勢(shì).總體來講，以耕地向林地、草地、水域用地、居住用地轉(zhuǎn)出為主；在與河流主干距離<10 km范圍內(nèi)耕地向林地轉(zhuǎn)化面積最大；交通用地和工礦用地在研究期間與其他地類幾乎無轉(zhuǎn)化.具體表現(xiàn)為，在與河流主干距離<10 km范圍內(nèi)主要地類轉(zhuǎn)化包括，耕地轉(zhuǎn)化為林地、草地、水域、居住用地、未利用地；林地轉(zhuǎn)化為水域、居住用地，草地轉(zhuǎn)化為水域用地；在與河流主干距離10～20 km范圍內(nèi)，耕地轉(zhuǎn)化為林地、草地、居住用地、未利用地，林地轉(zhuǎn)化為居住用地；在與河流主干距離>20 km范圍內(nèi)只有耕地向林地、草地、居住用地、未利用地轉(zhuǎn)出.

4 結(jié)論與討論

綜合考慮高程、坡度和與河流主干距離因素，表明2000～2010年在烏江北源流域區(qū)：（1）土地利用變化程度隨著高程、坡度的增加先增加再減小，隨與河流主干距離增加顯急劇減小的趨勢(shì).（2）土地利用變化以耕地向林地、草地和水域用地轉(zhuǎn)化為主，且發(fā)生在高程1 200～1 600 m、坡度15°～25°的與河流主干距離<10 km區(qū)域內(nèi)最為劇烈.（3）主要地類轉(zhuǎn)化包括：以耕地轉(zhuǎn)向林地、草地、水域用地為主，多發(fā)生在高程1 200～1 600 m、坡度15°～25°的與河流主干距離<10 km區(qū)域內(nèi)；林地、草地轉(zhuǎn)向水域用地，多發(fā)生在高程<1 200 m、坡度8°～15°的與河流主干距離<10 km區(qū)域內(nèi)；林地向居住用地轉(zhuǎn)化，多發(fā)生在高程1 200～1 600 m、坡度8°～15°的與河流主干距離<10 km區(qū)域內(nèi).（4）基本無轉(zhuǎn)化地類包括：交通用地、工礦用地與其他地類的轉(zhuǎn)化.

由于國(guó)家退耕還林和“西電東送”政策的實(shí)施，使得烏江北源流域區(qū)土地利用變化劇烈，以發(fā)生在與河流主干近距離的中海拔、較陡坡范圍內(nèi)耕地向林地、草地、水域用地的轉(zhuǎn)出和發(fā)生在與河流主干近距離的低海拔、緩坡區(qū)域內(nèi)耕地、草地和林地向水域用地的轉(zhuǎn)入最為明顯.退耕還林政策的不斷進(jìn)行，使得易造成水土流失的耕地向林草地轉(zhuǎn)型，生態(tài)環(huán)境得以改善；但由于喀斯特流域地形特點(diǎn)導(dǎo)致經(jīng)濟(jì)發(fā)展滯后，大量青壯年外出打工，主要?jiǎng)趧?dòng)力減少，使得部分耕地也向未利用地轉(zhuǎn)化；同時(shí)社會(huì)經(jīng)濟(jì)發(fā)展遲緩，使之交通用地、工礦用地面積較少且基本無變化.

烏江北源流域區(qū)土地利用變化受高程、坡度以及與河流距離影響較大，這種影響符合貴州喀斯特地貌下的自然環(huán)境和社會(huì)經(jīng)濟(jì)規(guī)律，遵循地形特征，充分利用高程、坡度以及水資源特點(diǎn)，結(jié)合區(qū)域社會(huì)經(jīng)濟(jì)狀況，合理、科學(xué)制定烏江北源流域區(qū)的土地利用方式和類型，以達(dá)到生態(tài)、社會(huì)和經(jīng)濟(jì)一體化的發(fā)展.

參考文獻(xiàn)：

[1] 郭洪峰，許月卿.基于地形因子的大都市邊緣區(qū)土地利用變化分析—以北京市平谷區(qū)為例[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2013，18（1）：178187.

[2] 梁發(fā)超，劉黎明.基于地形梯度的土地利用類型分布特征分析——以湖南省瀏陽市為例[J].資源科學(xué)， 2010，32（11）：21382144.

[3] 郭洪峰，許月卿，吳艷芳.基于地形梯度的土地利用格局與時(shí)空變化分析——以北京市平谷區(qū)為例[J].經(jīng)濟(jì)地理， 2013，33（1）：160166.

[4] 劉智勇，張 鑫，李小冰.禿尾河流域不同地形條件下土地利用動(dòng)態(tài)變化[J].西北林學(xué)院學(xué)報(bào)， 2012，27（3）：167173.

[5] 李 正，王 軍，白中科，等.喀斯特地區(qū)土地利用變化研究——以貴州省為例[J]. 地域研究與開發(fā)， 2011，30（2）：143148.

[6] 馮朝陽，于 勇，高吉喜，等.地形因子對(duì)京西門頭溝區(qū)土地利用/覆蓋變化的影響[J].山地學(xué)報(bào)， 2007，25（3）：274279.

[7] 楊格格，楊艷昭.南方紅壤丘陵地區(qū)土地利用變化特征[J].地理科學(xué)進(jìn)展， 2010，29（4）：483488.

[8] 馬廷剛，秦占飛，常慶瑞，等.縣域土地利用的地形特征分析——以陜西省戶縣為例[J].干旱地區(qū)農(nóng)業(yè)研究， 2011，29（6）：191221.

[9] 崔衛(wèi)國(guó)，文 倩，劉艷艷，等.基于DEM的醴陵市土地利用空間格局分析[J].資源科學(xué)， 2008，30（2），228234.

[10] 孫 麗，陳煥偉，潘家文，等.運(yùn)用DEM剖析土地利用類型的分布及時(shí)空變化——以北京延慶縣為例[J].山地學(xué)報(bào)， 2004，22（6）：762766.

[11] 寧艷梅，周文元，曾曉妹，等.湖南省衡南縣土地利用結(jié)構(gòu)的區(qū)域差異分析[J].熱帶地理， 2008，28（2）：172176.

[12] 熊康寧，白利妮，彭賢偉，等.不同尺度喀斯特地區(qū)土地利用變化研究[J].中國(guó)巖溶， 2005，24（1）：4147.

[13] 鄒 敏，吳泉源，逢杰武，等.基于DEM的龍口市土地利用空間格局與時(shí)空變化研究[J]. 測(cè)繪科學(xué)， 2007，32（6）：173175.

[14] 李 松，張小雷，李壽山，等.新疆綠洲城市土地利用景觀格局時(shí)空演化及驅(qū)動(dòng)力研究[J].經(jīng)濟(jì)地理， 2013，33（12）：161168.

[15] 趙銳鋒，陳亞寧，李紅衛(wèi)，等.塔里木河干流區(qū)土地覆被變化與景觀格局分析[J].地理學(xué)報(bào)， 2009，64（1）：95106.

[16] 余新曉，張曉明，牛麗麗，等.黃土高原流域土地利用/覆被動(dòng)態(tài)演變及驅(qū)動(dòng)力分析[J].農(nóng)業(yè)工程學(xué)報(bào)， 25（7）：219225.

[17] 鐘德燕，常慶瑞，宋豐驥，等.黃土丘陵溝壑區(qū)土地利用空間分布與地形因子關(guān)系研究[J].干旱區(qū)資源與環(huán)境， 2012，26（6）：102107.

[18] 黃海洋，楊慶媛，王成媛，等.基于DEM的土地利用類型與地形因子關(guān)系研究——以重慶市酉陽縣麻旺鎮(zhèn)為例[J].西南大學(xué)學(xué)報(bào)：自然科學(xué)版， 2009，31（4）：159164.

[19] 周元媛，何騰兵.喀斯特地區(qū)縣域土地利用變化及其原因分析：以貴州省黔西縣為例[J].貴州農(nóng)業(yè)科學(xué)， 2013，41（5）：168172.

[20] 韓建平，賈寧鳳.土地利用與地形因子關(guān)系研究——以磚窯溝流域?yàn)槔齕J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)， 2010，18（5）：10711075.