佛山市地形與土地利用空間關系研究

劉亮亮

（滄州市地產開發服務中心，河北 滄州 061000）

0 引言

地形因子中的高程、坡向、坡度是土地利用的重要影響因素，同時地形因子構成了土地利用結構的重要框架，所以，土地利用（Land Use）/土地覆蓋變化（Land Cover Change）（簡稱LUCC）研究的重要內容之一就是土地利用的時空變化與地形因子的關系。地形因子決定著區域的光照情況、土壤質地和水源涵養能力等，進而影響人們利用土地的方式和程度。因此，地形是土地利用的重要影響因素，并且地形差異是土地利用結構和空間分布格局分異的重要影響因素。綜上所述，分析地形對土地利用時空分布變化的影響，可以了解不同級別地形因子下各區域土地利用類型的空間分布變化。

我國學者對土地利用時空變化與地形因子的關系研究主要集中于高程、坡度、坡向與土地利用的關系，如分析地形因子對土地利用類型、產值及其經濟發展的影響，分析土地利用變化規律和驅動力，分析土地利用的適宜性等，較少綜合研究土地利用類型分布特征與高程、坡度、坡向、地形位指數和分布指數的關系。鑒于此，筆者分析廣東省佛山市土地利用類型分布特征，將其與高程、坡度、坡向、地形位指數和分布指數之間的關系進行空間連接，旨在梳理土地利用的空間變化與地形因子的關系，為區域土地資源的科學、可持續利用提供可靠的理論依據。

1 數據來源與研究方法

1.1 數據來源與處理

利用佛山市2011年的統計數據和佛山市的數字高程模型（Digital Elevation Model，DEM）數據，結合研究區域的土地利用現狀，將土地類型劃分為耕地、園地、林地、建設用地、水域和水利設施用地（簡稱水域）及其他用地6類。

佛山市海拔為0～740 m。為便于研究土地利用與各個地形因子之間的關系，筆者依據一定的社會經濟規律和自然規律，將佛山市的高程按照其在各個級別的分布密度分為7級，分別是0～20 m、20～60 m、 60～120 m、120～200 m、200～300 m、300～400 m和400～740 m。

坡度劃分依據《土地利用更新調查技術規定》，根據研究區域的土地利用類型和地形特征，將坡度分為平坡（0°～5°）、緩坡（5°～15°）、斜坡（15°～25°）、陡坡（25°～35°）和險坡（35°～90°）。

坡向共劃分為4個等級，以正北方向為0°，順時針方向度量，依次是陰坡（0°～45°、315°～360°）、半 陽 坡（45°～135°、225°～315°）、陽坡（135°～225°）、平緩坡（平坦無坡向地區）。

1.2 研究方法

在結果分析中，為方便分析土地利用空間分布與地形梯度之間的關系，筆者引入地形位指數來度量地形梯度，其計算公式如下：

根據式（1）推導可知，海拔越高、坡度越陡的地點，地形位指數越大；海拔越低、坡度越緩的地點，地形位指數越小；海拔較高但坡度較緩或海拔較低但坡度較陡的地點，地形位指數適中。

利用地形位指數分析土地利用空間分布與地形梯度之間的關系，受到坡度分級的客觀合理性、土地利用面積比例差異的影響，需要使用地形位分布指數消除這些影響，進而科學地描述不同土地利用類型在地形梯度上的分布特征。地形位分布指數計算公式如下：

式（2）中：P代表第種土地利用類型在地形位上的分布指數，代表地形因子；S為第種地形位下第種土地利用類型的面積，單位為km；S為整個研究區域第種土地利用類型的總面積，單位為km；為整個研究區域第種地形位的總面積，單位為km；為研究區域總面積，單位為km。地形位分布指數是一個標準化、無量綱的指數。如果P大于1，表示第種地類在第級地形位上的分布屬于優勢分布，P值越大，優勢度越高，該級別地形適宜該類土地發育和分布的程度較高；反之，如果P小于1，表示第種地類在第級地形位上的分布屬于劣勢分布，P值越小，優勢度越低，該級別地形適宜該類土地發育和分布的程度較低；如果P值等于1，則表示第種地類在第級地形位上的比重與研究區內該地類的比重相等。

2 結果與分析

2.1 土地利用類型分布特征隨高程梯度的變化

由圖1可知，隨著海拔的升高，除了林地外的其他各種土地利用類型的分布指數逐漸減小。海拔60 m以下耕地分布指數最大，其中海拔20 m以下其分布指數均大于1.00；隨著海拔的升高，各高程梯度上的分布指數逐漸減小，海拔300 m處其分布指數基本趨近于0。這表明在海拔0～60 m區域耕地出現頻率最大，是耕地的優勢分布區，且低海拔區域較適宜耕作。海拔60 m以下園地分布指數比較穩定，均在1.00以上，而后隨海拔升高迅速降低，海拔大于300 m時其分布指數趨向于0，表明高程梯度對園地的限制作用顯著。海拔20 m以下林地分布指數都不超過0.50，海拔20～60 m區域林地分布指數呈現顯著增大趨勢，海拔大于60 m的區域林地分布指數緩慢增大。建設用地和水域分布指數類似，海拔小于20 m的區域其分布指數均大于1.00，海拔20～60 m區域其分布指數迅速下降至0.30左右，海拔200 m以上區域其分布指數基本接近于0，說明建設用地和水域的分布指數有緊密的聯系，有高度的空間自相關性。隨著海拔的升高，建設用地分布指數下降迅速，說明其分布指數與高程存在極大的關聯性，同時表明建設用地的選址應盡可能在低海拔區。海拔60 m以下其他用地分布指數均大于1.00，海拔120 m以上其分布指數迅速減小，并趨向于0。這說明佛山市未利用土地在海拔60 m以下分布密集，便于開發利用，可開發利用的土地資源有沙地、裸地、草地、設施農用地。

圖1 各高程區土地利用類型分布指數

2.2 土地利用類型分布特征隨坡度梯度的變化

坡度是土地利用的重要影響因素，尤其對建設用地、水域、耕地的分布影響顯著，陡坡和險坡對人類的各種活動都有較大的限制性。由表1可知，平坡區建設用地、水域面積最大，占該區土地利用總面積的76.60%，是該區域的主導土地利用方式。隨著坡度的增大，其他用地在緩坡區所占比例達到峰值，這主要是因空閑地、設施農用地、草地增加所致。同時，緩坡區林地所占比例迅速增大，林地、建設用地和水域共同成為該區域的主導土地利用方式，占該區土地利用總面積的83.93%，其中建設用地、水域下降幅度最大，主要是由于隨著坡度的增大，地形對人類活動的限制性增強。在斜坡以上坡度區間，林地所占比例在70%以上，呈絕對優勢分布。

表1 不同坡度分級區域土地利用類型分布特征

由圖2可知，在坡度小于5°的區域，耕地面積分布最大，是該坡度梯度上的主導土地利用方式。耕地分布指數隨坡度增大而逐漸減小的原因主要是坡度越大，人類從事農業生產活動的難度越大，投入與產出比越小，所以耕地主要分布在比較平坦的區域。在坡度小于35°的區域，耕地分布指數均大于0.60；在坡度大于35°的區域，其分布指數為0.34。這說明當地耕地資源相對缺乏，人地矛盾比較突出，為了求生存，人們會在坡度大于25°的區域種植糧食作物。在坡度小于15°的區域，園地分布指數都大于1.00；當坡度大于15°時，隨著坡度的增加，園地的分布指數逐漸減小，坡度大于35°時園地分布指數僅為0.17，說明園地對坡度的適宜性差。在坡度小于5°的區域，林地分布指數小于1.00；當坡度大于5°時，其分布指數迅速增大且大于1.00，是該坡度梯度上的主導土地利用方式，且隨著坡度的增加，其分布比例逐漸增大，表明林地對坡度的限制條件有較好的適宜性。在坡度小于5°的區域，水域分布指數大于1.00，而后迅速遞減。在坡度小于15°的區域，建設用地分布指數大于0.83，而后遞減趨勢明顯增強，主要是由于隨著坡度的增加，人類生產建設的難度增大。其他用地主要分布在25°以下的區域，其中15°以下區域其分布指數大于1.00，表明其他用地在15°以下的區域集中比較分布。

圖2 各坡度區土地利用類型分布指數

2.3 土地利用類型分布特征隨坡向梯度的變化

坡向作為地形因子的重要因素，對土地利用中的光照條件影響顯著。由圖3可知，耕地主要分布在平緩坡、陽坡、半陽坡，分布指數大于0.55，其中平緩坡其分布指數都大于1.00。這3個坡向梯度區域為耕地的坡向優勢區，占該區域土地總面積的96.5%，表明耕地對光照、熱量、水分條件的要求較高。園地這與耕地相似，也主要分布在平緩坡、陽坡、半陽坡，其分布指數都大于0.70，這3個坡向梯度區域同樣為園地的坡向優勢區。林地和其他用地的分布幾乎不受坡向限制，在各個坡向上其分布指數變化都不大，也就是說林地和其他用地的分布基本上沒有方向性，說明林地和其他用地對各個坡向具有較強的適宜性，與耕地、園地的適宜條件相反，林地和其他用地受光熱水的影響較小。建設用地主要分布在平緩坡，其分布指數大于1.00，在陽坡、半陽坡的分布指數大于陰坡，但差別不大，說明建設用地對坡向的適宜性較強。水域主要分布在平地，在各坡向上變化幅度不大。

圖3 各坡向區土地利用類型分布指數

2.4 土地利用類型分布特征隨地形位梯度的變化

由圖4可知，佛山市地形位梯度在0～1.60，當地形位梯度為0～0.60時，耕地的分布指數均大于0.60，其中0～0.20區間其分布指數大于1.00，耕地分布優勢明顯；當地形位梯度大于0.60時，耕地分布指數迅速減小，最終其分布指數在地形位指數為1.20～1.40時趨向于0。這表明在低地形位梯度上，耕地存在分布優勢，而在較高地形位梯度上，耕地對高程和坡度的適宜性差，耕地分布指數小。當地形位指數為0～0.20時，林地分布指數為0.27；地形位指數大于0.21時，林地分布指數均大于2.00。這說明林地在低地形位梯度上受耕地、建設用地、水域等分布的影響較大，而在較高地形位梯度上林地分布存在明顯優勢，表明林地比耕地、園地、建設用地對高程和坡度的適宜性強。當地形位梯度為0～0.40時，建設用地分布指數大于0.70，其中地形位梯度為0～0.20時建設用地和水域分布指數大于1.00。這表明人類的生產建設活動對高程和坡度的依賴性強，在建設用地分配使用中需重點考慮地形位梯度和水域分布。當地形位梯度為0～0.60時，其他用地分布指數大于0.58，其中地形位梯度為0～0.40時其分布指數均大于1.00，尤其是地形位梯度為0.21～0.40時，除了林地外，其分布指數最大，表明在較低地形位梯度上其他用地優勢區位明顯。在以后的生產建設中，應注重科學合理開發其他用地，注重開發與生態建設并行。

圖4 各地形位梯度土地利用類型分布指數

3 結論與討論

在一定程度上，地形是土地利用空間分布差異的決定因素。土地利用類型分布特征隨地形梯度的變化能反映該地域的地形條件，反映各個土地利用類型的分布優勢區位，反映各級地形梯度上的土地主導利用方式，也能闡述各種土地利用類型的適宜條件、限制因素及未來的發展方向。

在大尺寸地形位梯度上，佛山市土地利用類型分布特征存在明顯差異：耕地、園地和其他用地主要分布在第一、二、三級地形位梯度上，表明耕地、園地、其他用地對高程和坡度具有較差的土地適應性。其原因主要是低海拔區域光、熱、水資源豐富，比較適宜耕作；坡度越大，人類從事農業生產活動和生產建設的難度越大，投入與產出比越小，所以耕地、園地、其他用地主要分布在比較平坦的區域。同種土地類型在高程、坡度、坡向、地形位梯度上的分布指數具有一定的相似性。林地在較低高程梯度等級區域上分布指數較低，而在中高高程梯度等級區域上具有明顯分布優勢；林地在較低坡度梯度上分布指數較低，而在其他坡度梯度上分布優勢明顯；林地在坡向梯度上的分布特征不存在顯著變化；林地在地形位梯度上的分布指數與在高程梯度、坡度梯度上的分布指數相似，研究區域林地分布指數的變化不大。這表明林地對高程、坡度、坡向、地形位梯度具有較好的適應性，在較低高程梯度、較小坡度梯度、陽坡上的分布受耕地、園地分布擠壓明顯。

從分布指數的變化角度來看，耕地、園地、建設用地、水域、其他用地的分布指數隨著地形位梯度的增加而逐漸減少，并趨近于0，表明其對高程和坡度具有較差的土地適應性；隨著地形位梯度的升高，林地分布指數逐漸增大，并趨向于各級地形位梯度上的優勢分布，除了表明林地對地形位因子有較好的適應性外，還表明當地水土保持良好、生態效益明顯，各個景觀組分在低、中、高地形位梯度上呈現出良好的立體組合。

該研究采用高程梯度、坡度梯度、坡向梯度、地形位梯度綜合分析土地利用類型分布特征，大體上各地形梯度分級區域的土地利用類型分布指數具有一致性。在同一地形位梯度上，坡度等級上的分布指數和高程等級上的分布指數存在一定程度的差異。單個地形因子雖能反映土地利用類型的分布，但是地形位因子更能反映土地利用類型的實際分布趨勢，避免土地利用類型選擇時既要考慮高程梯度也要考慮坡度梯度兩方面因素的煩瑣，不僅減少了工程量，也提高了規劃質量與精度。該研究基于地形位梯度分析土地利用類型分布特征，定量分析了土地利用變化規律，可為佛山市土地利用規劃提供依據。