陸清平, 趙翠薇, 王 杰

(1.貴州師范大學 地理與環(huán)境科學學院, 貴陽 550025; 2.喀斯特山地生態(tài)環(huán)境保護與資源利用協(xié)同創(chuàng)新中心, 貴陽 550025)

土地利用/覆被與地球大氣、土壤、水文和生物及物質(zhì)循環(huán)等息息相關[1-2]，對人類的生存環(huán)境和發(fā)展影響深刻持久[3-4]，已成為全球生態(tài)環(huán)境變化的熱點研究[5-7]。地形是土地利用格局形成的基礎，在一定程度上決定了土地利用的結構和方向[8-9]，使得土地利用在地形上呈現(xiàn)明顯梯度特征[10]。

近年來，土地利用結合地形的研究較為豐富，Liyew Birhanu等探究了埃塞俄比亞土地利用與海拔和坡度的關系[11]，李京京等運用綜合高程和坡度的地形位指數(shù)分析了汾河流域土地利用的地形梯度演變規(guī)律[12]。從研究內(nèi)容看，主要包括地形的土地利用圖譜變化[13-14]和景觀生態(tài)空間風險及效應分析等[15-17]；研究區(qū)域主要集中在生態(tài)脆弱區(qū)、流域及高原山地等[18-21]，且國外更加偏重對人為因素的探究。對喀斯特山區(qū)研究有林地的地形梯度變化[22]，城市建設用地擴展梯度特征[23]和流域的土地利用地形梯度演化規(guī)律等[24]，然都是將研究區(qū)整體作為純喀斯特或非喀斯特地貌分析其地類梯度特征，而Cvijic[25]、安裕倫[26]、賀中華[27]等發(fā)現(xiàn)，純喀斯特與非喀斯特間還存在與兩者景觀差異較大的“半喀斯特”或“過渡喀斯特”地貌，即亞喀斯特地貌[28]。

當前，對純喀區(qū)、非喀區(qū)及亞喀地貌區(qū)的土地利用地形梯度格局研究較匱乏，3種地貌區(qū)土地利用地形梯度特征是否存在差異，若有差異，差異體現(xiàn)在哪些方面尚不清楚?；诖耍再F州黔中地區(qū)為研究區(qū)，根據(jù)巖性劃分喀斯特區(qū)、非喀區(qū)和亞喀區(qū)，運用地形位指數(shù)和優(yōu)勢分布指數(shù)，揭示土地利用類型在不同地貌區(qū)的地形梯度規(guī)律，從而為統(tǒng)籌國土資源的合理空間配置提供基礎數(shù)據(jù)支撐，也為區(qū)域因地制宜開展資源精細和差異化管理及可持續(xù)發(fā)展提供參考。

1 研究區(qū)概況

研究區(qū)位于貴州中偏西和南部，包括安順市全境，外加六枝特區(qū)、織金縣、清鎮(zhèn)市和花溪區(qū)，共1個縣級市、4個區(qū)和7個縣，總面積為20 283.80 km2。地勢為西北高東南低，海拔范圍314～2 225 m，地形高差較大，其中西北部織金縣和西部六枝特區(qū)海拔最高，處于貴州西部第一階梯向中部第二階梯的過渡帶(貴州地勢自西向東分為三級階梯[29])，河流縱橫且河谷深切；中部西秀區(qū)、東部花溪區(qū)和北部清鎮(zhèn)市為貴州高原面，東南部和南部區(qū)域海拔最低，地處第三階梯地勢帶范圍。巖性以連續(xù)型石灰?guī)r為主，還伴隨著大量的非碳酸鹽巖等；地質(zhì)成分以灰?guī)r碎屑巖互層和灰?guī)r夾碎屑巖組合為主。地處亞熱帶季風區(qū)，年均溫、年均降雨量分別在14～16℃和1 300～1 400 mm，雨熱同季利于農(nóng)業(yè)生產(chǎn)。土地利用以林地、耕地和草地為主，其中林地面積最大，占研究區(qū)的38.67%，耕地次之，占36.52%，草地占10.32%，建設用地達4.47%。2017年GDP總量2 098.58億元，占全省的14.17%，常住人口達480.60萬人，城鎮(zhèn)化率53.32%略高于貴州省(46%)，是貴州城鎮(zhèn)化及經(jīng)濟發(fā)展水平較高的區(qū)域。

2 材料與方法

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

研究基礎數(shù)據(jù)包括土地利用和DEM(數(shù)字高程模型)及巖性數(shù)據(jù)。土地利用數(shù)據(jù)由貴州省2017年土地調(diào)查年度變更成果數(shù)據(jù)庫裁剪獲得；巖性數(shù)據(jù)由1∶20萬圖件掃描、精校正、矢量化等處理后得到；坡度和高程數(shù)據(jù)由日本地球觀測衛(wèi)星ALOS(分辨率為12.5 m，重采樣為30 m)網(wǎng)站(https:∥search.asf.alaska.edu/)下載DEM，在ArcGIS空間分析工具中提取得到。

根據(jù)2017年全國土地利用現(xiàn)狀分類標準，結合研究區(qū)實際需要，分為草地、旱地、有林地、灌木林地、水域(水庫和坑塘及河流等水面)、未利用地(沼澤地、沙地、內(nèi)陸灘涂和裸地)、農(nóng)村居民點、其他林地、水田、園地(果園和其他園地)、道路、基礎設施用地(港口碼頭、溝渠、管道運輸、設施農(nóng)用地、機場和水工建筑用地)、工礦用地、城鎮(zhèn)建設用地和水澆地共15種地類。利用ArcGIS柵格計算器計算研究區(qū)地形位指數(shù)并采用自然斷點法進行分級；最后將土地利用與地形位指數(shù)和地貌圖疊加，獲得不同地貌區(qū)各地類的地形位指數(shù)。

2.2 地貌類型劃分

安裕倫等認為巖性是區(qū)分純喀斯特、非喀斯特和亞喀斯特地貌的主要因素[28]，蘇琪嬌等總結出3種類型分布特征，即純喀斯特在純碳酸鹽巖間發(fā)育強烈，亞喀斯特在不純碳酸鹽巖區(qū)集中，非喀斯特與非碳酸鹽巖密切相關[30]。因此本次研究以1∶20萬巖性數(shù)據(jù)為基礎，將巖性進行合并歸類(表1)，最終獲得研究區(qū)3種喀斯特地貌類型分布范圍(圖1)。

表1 地貌類型劃分標準及面積

圖1 地貌類型分布

2.3 研究方法

2.3.1 地形位指數(shù) 單一因子難以全面表征地形特征[31]。土地利用類型分布與高程和坡度有關，故采用坡度和高程組成的地形位指數(shù)綜合表征地形特征。計算公式如下：

T=lg[(E/E0+1)×(S/S0+1)]

(1)

式中:E,E0分別表示高程及平均高程值；S及S0表示坡度和坡度均值，坡度和高程增加時,T隨之增大,相反則會減小,而當坡度增大高程減小,或者高程增加坡度減小時,T則處于前兩者的中間。

2.3.2 分布指數(shù) 地形對地類分布產(chǎn)生的影響，實質(zhì)可認為是各梯度上地類的分布頻率，故用分布指數(shù)來表達景觀分布特征，還可消除地形梯度下各個土地利用類型面積比例差異而造成的量綱差異。公式如下：

P=(Sie/Si)/(Se/S)

(2)

式中:P為分布指數(shù),分布優(yōu)勢和p值兩者具有正相關性;Sie表示i地類在e地形內(nèi)的面積;Si為i地類的總面積,Se表示e地形的總面積;當p=1時,表明e地形內(nèi)i地類的面積占比與研究區(qū)內(nèi)i的面積比例相等;p<1時,該地類在e梯度內(nèi)呈劣勢分布,反之為優(yōu)勢分布,且P越大,分布優(yōu)勢越明顯。

3 結果與分析

3.1 地形位指數(shù)梯度分區(qū)

根據(jù)公式(1)獲得研究區(qū)地形位指數(shù)分布圖，采用自然斷點法分為10個等級，并歸為低地形梯度、一般地形梯度、較高地形梯度和高地形梯度4個類型，統(tǒng)計出各地形位等級的面積(表2)。統(tǒng)計結果顯示，研究區(qū)以3～7地形梯度為主，3～7級各地形梯度面積比例均大于10%，合計面積占比達68.26%。低地形梯度面積占15.22%，一般地形梯度區(qū)占44.58%，較高地形梯度區(qū)占27.70%，高地形位面積比例最小，僅2.50%。高地形位分布在西北和西部及南部區(qū)域，以織金、六枝和紫云最顯著，織金縣、六枝特區(qū)海拔高、坡度大，紫云縣處于貴州高原向廣西丘陵過渡地帶，地形破碎，坡度較大，地形位指數(shù)大；平壩、花溪、西秀區(qū)和普定縣南部及惠水縣一帶，為貴州高原面，地勢較為平坦，為喀斯特溶蝕盆地、谷地，雖然海拔較高但坡度小，成為低地形位集中區(qū)域(圖2)。

圖2 研究區(qū)地形位指數(shù)

表2 地形梯度區(qū)間及面積統(tǒng)計

3.2 不同地貌區(qū)土地利用的地形分布

3.2.1 純喀斯特區(qū) 純喀斯特地貌類型在研究區(qū)域內(nèi)分布范圍較廣，占研究區(qū)總面積的55.57%。根據(jù)公式(2)和Origin 2017軟件獲得研究區(qū)不同地形梯度的分布指數(shù)(圖3)。

由圖3可知，水域分布指數(shù)變化最大，在第1梯度內(nèi)指數(shù)達4.57，在1～3梯度內(nèi)為優(yōu)勢分布，隨著地形梯度的增加分布指數(shù)逐漸降低，說明水域集中分布在低地形梯度區(qū)，在地形位高的區(qū)域水域面積較少。水田在1～3梯度為優(yōu)勢分布，水澆地和園地是1～4梯度，最大分布指數(shù)分別是3.60，2.66及1.99，表明3種地類對坡度海拔低的區(qū)域選擇性較強，主要分布在西秀、平壩、普定和惠水等地。旱地在第5梯度有最大值1.16，優(yōu)勢區(qū)間為2～6，說明其對地形選擇性較小，能在較高梯度有較好的適應性[19]。有林地和其他林地在4～8梯度為優(yōu)勢分布區(qū)，與地形呈現(xiàn)先增后減趨勢，多分布在地形位一般區(qū)域(4～6梯度)，有林地在第8梯度時分布指數(shù)為1，表明第8梯度內(nèi)的有林地與在研究區(qū)內(nèi)的比例相同，是其優(yōu)勢向劣勢的分界點，而其他林地分界點在8～9地形梯度間。灌木林地在1～5為劣勢分布，6～10梯度為優(yōu)勢分布，與梯度變化呈正相關關系，在較高地形位適宜性較好[32]。草地地形梯度特征呈倒“U”型變化，第4梯度前為劣勢而5～10呈優(yōu)勢分布，第7梯度達到優(yōu)勢分布峰值1.64，后又逐漸降低。農(nóng)村居民點和城鎮(zhèn)建設用地在1～3梯度為優(yōu)勢分布，在第1梯度有最大指數(shù)2.37，4.05，均隨著梯度增加而減小，農(nóng)村居民點由優(yōu)勢向劣勢轉(zhuǎn)變發(fā)生在4～5地形位，而城鎮(zhèn)建設用地過渡區(qū)間為3～4，說明城鎮(zhèn)用地對地形的選擇性比農(nóng)村居民點強。未利用地分布指數(shù)與梯度表現(xiàn)為正相關性，6～10梯度為優(yōu)勢分布，第10梯度有最大分布值2.06，主要集中在紫云西北部、東部和關嶺中部。

圖3 純喀斯特區(qū)域土地利用類型在各地形梯度上的分布指數(shù)

3.2.2 非喀斯特區(qū) 非喀斯特地貌類型在研究區(qū)分布范圍較小，主要集中在鎮(zhèn)寧和紫云，在織金、清鎮(zhèn)和六枝及花溪有少許分布(圖4)，占研究區(qū)總面積的21.48%。

由圖4可知，有林地優(yōu)勢分布為5～9地形梯度階段，草地為5～10，灌木林地和其他林地為6～10，除有林地外，其余均與地形梯度呈正相關，在第10梯度有最大分布值，分別是1.35，1.67，2.12，有林地在第7梯度時就達到最大分布值1.15，但幾類生態(tài)用地均分布在較高地形位區(qū)域。農(nóng)村居民點、水田、園地、道路、基礎設施用地、城鎮(zhèn)建設用地、水澆地及水域均隨著地形梯度增加而減小，在第1梯度均有最大分布值，城鎮(zhèn)建設用地和水域高達7.36，6.18，其他幾類均在2～2.5之間。旱地在3～6梯度為優(yōu)勢分布，最大分布值1.05位于第4梯度，在地形梯度上表現(xiàn)為先增后減特征。未利用地與地形位具有正相關性，7～10地形位為優(yōu)勢分布，其余階段為劣勢分別，第10梯度有最大值3.24。

圖4 非喀斯特區(qū)域土地利用類型在各地形梯度上的分布指數(shù)

3.2.3 亞喀斯特區(qū) 亞喀斯特地貌類型在研究區(qū)內(nèi)分布范圍較小，主要集中在惠水、長順、花溪、西秀、普定和織金等地方，其他區(qū)域也有少許分布(圖5)，占研究區(qū)總面積的22.95%。

圖5 亞喀斯特區(qū)域土地利用類型在各地形梯度上的分布指數(shù)

由圖5可知，工礦用地、道路、城鎮(zhèn)建設用地和水澆地在1～3地形梯度為優(yōu)勢分布，工礦和城鎮(zhèn)建設用地在第2梯度有最大值3.02，3.31，道路在第1梯度時有最大值2.41，3種地類優(yōu)勢分布主要集中在六枝、普定、西秀等縣城中心和清鎮(zhèn)市東南部及花溪東北部區(qū)域；農(nóng)村居民點、水田和基礎設施用地在1～4梯度為優(yōu)勢分布，水田在第1梯度有最大值2.47，而農(nóng)村居民點和基礎設施用地最大值在第2梯度，分別是2.03，1.85；園地優(yōu)勢分布區(qū)域為1～4梯度，在第2梯度有最大值1.80；水域優(yōu)勢區(qū)間為1～2梯度，最大優(yōu)勢分布值在第1梯度，達9.49，兩區(qū)間分布值相差8.16；上述9地類均與地形梯度呈現(xiàn)負相關關系，水域變化最為明顯，表明在亞喀斯特區(qū)域，這幾種土地利用類型對低地形的選擇性最強，如水田在平壩、西秀和惠水等地分布較廣。灌木林地和未利用地在6～10梯度為優(yōu)勢分布，與梯度成單調(diào)遞增關系，在第10梯度有最大值2.39，1.66，而旱地、有林地、草地和其他林地均是先增加后減小，旱地在3～7梯度為優(yōu)勢分布，草地為4～9，有林地是4～7，而其他林地則為5～10，分布指數(shù)均沒有較大差異，但旱地分布指數(shù)更接近于1，說明旱地在3～7梯度的比重接近于整個亞喀斯特內(nèi)旱地的比重。

3.3 研究區(qū)土地利用分布指數(shù)差異

純喀斯特與非喀斯特及亞喀斯特存在著土地利用地形的差異，主要表現(xiàn)在：(1) 土地利用從優(yōu)勢向劣勢以及從劣勢向優(yōu)勢分布過渡的地形梯度帶不同，如有林地在純喀斯特地貌過渡地形梯度為4和8，非喀斯特為5和9，亞喀斯特則是4和7，其對非喀斯特地形適應更強。(2) 土地利用類型變化趨勢和速率不同，變化趨勢上如農(nóng)村居民點在純喀與非喀區(qū)是單調(diào)遞減，而亞喀區(qū)則是先增后減趨勢；在變化速率上如旱地在純喀區(qū)由0.83增加至1.16后降至0.63，非喀則是從0.64增加到1.08在降至0.94，而亞喀區(qū)則是由0.66上升至1.05再降至0.92，旱地在非喀區(qū)變化速率低于亞喀區(qū)，更遠低于純喀區(qū)。(3) 耕地的優(yōu)勢分布地形位面積占比順序是：純喀區(qū)(26.26%)<亞喀區(qū)(26.61%)<非喀區(qū)(26.64%)，水田在純喀區(qū)優(yōu)勢區(qū)間為1～3，亞喀和非喀區(qū)則為1～4，同時亞喀區(qū)在優(yōu)勢向劣勢的過渡值(1.07)比非喀區(qū)(1.12)小，旱地在純喀區(qū)優(yōu)勢分布最小梯度為第2級，亞喀和非喀區(qū)則為第3等級且由劣勢向優(yōu)勢過渡值分別是1.03，1.05，非喀區(qū)優(yōu)勢值大于亞喀區(qū)。綜上表明耕地在非喀區(qū)地形選擇性小于亞喀區(qū)，更低于純喀區(qū)。

4 討 論

4.1 基于地形梯度的土地利用一般規(guī)律

通過對不同梯度上地類及其分布指數(shù)進行統(tǒng)計分析，發(fā)現(xiàn)存在區(qū)域差異性的同時，也有共同的分布特征。

(1) 無論是純喀、非喀還是亞喀區(qū)，農(nóng)村居民點、水域水田、園地、水澆地和基礎設施及城鎮(zhèn)建設用地等7類景觀，均隨梯度增加逐漸從優(yōu)勢向劣勢分布過渡，都在1～2梯度內(nèi)呈現(xiàn)優(yōu)勢分布，水域在第1梯度有最大優(yōu)勢分布值，其余6類在第2梯度分布范圍最廣，且水田為1～2梯度內(nèi)生活生產(chǎn)的主要利用類型，原因是貴州省平地本底稀少[33]，水田對地形和水資源條件要求較高。旱地在3～6梯度為優(yōu)勢分布，但其優(yōu)勢分布指數(shù)均趨近于1，表明旱地在此區(qū)間的分布比重與整個研究區(qū)旱地的比重相近，與梯度呈現(xiàn)先增后減的變化，同時旱地在2～7梯度面積分布較多，與貴州1991年開始的坡改梯工程有密切關系[34]。(2) 草地、灌木林地、未利用地、有林地和其他林地優(yōu)勢分布均在較高地形位，除有林地外其他景觀類型均隨梯度增加而增加，有林地表現(xiàn)為先增加后減小趨勢，表明有林地適宜性范圍小于這幾種地類，這與植被的垂直地帶性分異有關，隨著地形位的增加，坡度和高程逐漸增加，溫度逐漸下降，有林地生長環(huán)境受限，被灌木和草地及未利用地(裸地)等所替代，第8梯度后氣溫持續(xù)降低，草地優(yōu)勢范圍逐漸縮小，發(fā)生向灌草叢和未利用地的演替。(3) 總體而言，地形位較小區(qū)域，生活條件較好，適宜人類從事生產(chǎn)生活，隨著地形位增加，環(huán)境逐漸變差，人類發(fā)展受限，適宜生態(tài)用地的發(fā)展，具體表現(xiàn)就是生產(chǎn)生活用地的分布指數(shù)逐漸趨近于0，而生態(tài)用地分布指數(shù)與之相反。

4.2 土地利用類型空間分布影響因素

土地利用類型空間分布是由社會經(jīng)濟和政府政策及自然因素綜合作用的結果[35]，坡度高程等決定其空間分布格局[36]。研究區(qū)土地利用具有顯著的地形梯度特征，從自然角度出發(fā)，純喀區(qū)經(jīng)過千萬年形成復雜地表地下二元結構，土層淺薄、保水性差，而亞喀區(qū)則擁有較厚的土層，土壤保肥性好覆被度高，宜農(nóng)耕地多，且非喀區(qū)最為理想[28]。研究表明非喀區(qū)最大持水率均值>純喀區(qū)>亞喀區(qū)，且同等條件下亞喀區(qū)土壤侵蝕量最高[37]，這就是為什么在純喀區(qū)和非喀區(qū)水域優(yōu)勢地形位區(qū)為1～3，而亞喀區(qū)為1～2的原因。人類生產(chǎn)活動主要在地形位較小的區(qū)域，對高海拔坡度大區(qū)域的涉及較少，所以林草地及未利用地得以在地形位較大區(qū)域占據(jù)優(yōu)勢。社會經(jīng)濟涵蓋城鎮(zhèn)化、經(jīng)濟、交通建設和對生活高質(zhì)量追求等，是土地利用分布的內(nèi)在動力[38]，由于工業(yè)和城鎮(zhèn)化進程的加快，經(jīng)濟發(fā)展集聚效應誘使人們?yōu)樵黾邮杖脒h離農(nóng)村，此外，生活質(zhì)量的不斷上升，土地屬性不斷變化，城市草地、林地和公園綠地面積增長，使得城鎮(zhèn)建設用地在低地形位不斷向外擴張，并逐漸占據(jù)優(yōu)勢。政策方面受退耕還林和精準扶貧的影響，大于25°的耕地嚴格實行退耕，造成高海拔坡度大的區(qū)域出現(xiàn)向森林和草地及未利用地轉(zhuǎn)型現(xiàn)象[39]；同時，異地扶貧搬遷，人口向低地形位聚集，高地形位大量耕地撂荒，發(fā)生自然向草地、灌木和其他林地(疏林地)演替，使得林草地的分布優(yōu)勢增加。

4.3 研究特色及不足

景觀類型格局研究成果豐碩，然而已有研究，只是簡單指出其空間分布，較少進行對比不同地貌的空間格局研究，尤其是不同喀斯特地貌景觀類型地形梯度的對比探究。本文特色之處在于將景觀與不同地貌類型和地形梯度三者聯(lián)系起來，指出地類在各地貌的不同地形位上的分布特征。不足之處在于，不同地貌上景觀類型的地形梯度影響因素，只是簡單從定性方面進行闡述，沒有定量探究，同時也未開展長時間序列的演化特征分析，對地形梯度上不同地貌的景觀類型演化規(guī)律和影響因素定量研究是下一步工作的重點。

5 結 論

研究區(qū)土地利用具有顯著的地形梯度特征，農(nóng)村居民點、水田、園地、道路、基礎設施用地、城鎮(zhèn)建設用地、水澆地和水域等在1～3地形梯度有顯著分布優(yōu)勢，有林地、灌木和草地及未利用地在7～10地形梯度優(yōu)勢分布突出，而旱地優(yōu)勢主要集中于低和一般地形位。研究區(qū)不同地貌的城鎮(zhèn)建設用地、基礎設施用地和灌木林地具有相同的優(yōu)勢分布地形位區(qū)間，其他土地利用類型雖在整體上表現(xiàn)出較強相似的地形梯度特征，但不同地貌內(nèi)的相同景觀類型地形梯度特征具有差異性，在優(yōu)勢分布范圍變化和趨勢及速率有所不同。有林地對非喀區(qū)地形適應范圍最廣(地形位范圍5～9)、純喀區(qū)次之(4～8)，亞喀區(qū)較窄(4～7)；草地在純喀區(qū)和非喀區(qū)都集中在地形位5～10范圍，亞喀區(qū)最差(4～9)；純喀區(qū)水田對地形選擇性最強(1～3)，亞喀和非喀區(qū)選擇性較弱(1～4)；農(nóng)村居民點在純喀區(qū)與非喀區(qū)均呈現(xiàn)單調(diào)遞減趨勢，在亞喀區(qū)則表現(xiàn)為先增加再單調(diào)減??；旱地和未利用地在非喀區(qū)的優(yōu)勢分布變化速率低于亞喀區(qū)，更遠低于純喀區(qū)，而水域卻恰恰與之相反。