甘肅省榆中縣地形對(duì)耕地分布影響

甘肅省榆中縣地形對(duì)耕地分布影響摘要：針對(duì)現(xiàn)有研究在地形對(duì)耕地分布影響中較少考慮地形因子間交互作用，且忽視了溝壑密度等因子影響的問(wèn)題,以甘肅省榆中縣為例，在地貌分區(qū)基礎(chǔ)上利用分形指數(shù)分析了耕地分布特征，并通過(guò)地理探測(cè)器分析高程、坡度、坡向和溝壑密度對(duì)耕地分布的影響。結(jié)果表明，河谷平原區(qū)與黃土塬區(qū)耕地連片性較好，河谷平原區(qū)、黃土丘陵區(qū)和北部石質(zhì)低山丘陵區(qū)耕地空間分布較穩(wěn)定;各地貌區(qū)影響耕地分布的主導(dǎo)因子不同，南部石質(zhì)山地和黃土丘陵區(qū)是高程，黃土塬和河谷平原區(qū)是坡度，北部石質(zhì)低山丘陵區(qū)是溝壑密度;地形因子交互作用均表現(xiàn)為增強(qiáng),其中溝壑密度影響作用不可忽略。關(guān)鍵詞：耕地分布；地形；地理探測(cè)器；溝壑密度；榆中縣EffectofTopographyonCultivatedLandDistributionin

YuzhongCountyofGansuProvinceAbstract:Aimatexistingstudiesontheeffectsoftopographyincultivatedlanddistributionseldomanalyzeeffectofinteractionofterrainfactors,andignoreeffectofgullydensity,thispapertakesYuzhongcountyofGansuprovinceasanexample,analyzesfeatureofcultivatedlanddistributiononthebasisofgeomorphiczoningwiththehelpoffractalindex,andanalyzestheeffectsofelevation,slope,aspectandgullydensityoncultivatedlanddistributionthroughgeographicaldetector.Theresultsshowthat：①Theintegrityoffarmlandintherivervalleyplainandtheloesstablelandisgood;thespataldistributionoffarmlandisstableintherivervalleyplan,loesshillregonsandnorthstonyhillyarea.②Thedomnantfactorsinfuencingthedistributionofcutivated6andarevariousindifferentgeomorphicareas.Thedominantfactortosouthstonymountainsand6oesshi6regions,the6oesstabeandandriverva6eypains,andthenorthstonyhi6yareaiseevation,s6ope,andgullydensity,respectively.③Theinteractionofterrainfactorsisenhanced,andtheeffectofgullydensitycannotbeignored.Keywords'distributionofcultivatedland；terrain；geographcdetector；gullydensty；Yuzhongcounty0引言耕地是土地利用的重要類型之一，也是人類賴以生存的重要場(chǎng)所和生產(chǎn)資料，不僅擔(dān)負(fù)著為人類提供食物保障的重任，從人類生存與可持續(xù)發(fā)展的角度來(lái)看，也是國(guó)家生態(tài)安全的重要保障耕地

的形成與存在受氣候、地形、土壤等條件的影響，同時(shí)高程、坡度等在很大程度上決定著局部環(huán)境中的溫度、降雨和水分等條件，進(jìn)而影響耕地的開發(fā)難度、質(zhì)量及其分布⑵。因此，分析地形對(duì)耕地分布的影響具有重要的理論和現(xiàn)實(shí)意義（現(xiàn)有地形對(duì)耕地分布影響研究主要在地貌類型分區(qū)、地形因子［5=7］或地形綜合指數(shù)分級(jí)基礎(chǔ)上#，運(yùn)用分形理論分布指數(shù)［11=12］和耕地動(dòng)態(tài)度⑵等方法，討論不同分級(jí)區(qū)域中耕地分布的時(shí)空特征。研究主要發(fā)現(xiàn)，耕地隨地形因子等級(jí)變化分布具有一定的梯度性，在地勢(shì)平坦的地貌類型區(qū)和地形因子低級(jí)別區(qū)分布面積較大，且在地勢(shì)起伏度小的地區(qū)連片性好、穩(wěn)定性強(qiáng)。也有學(xué)者在定性分析的基礎(chǔ)上，采用相關(guān)系數(shù)#;!$、回歸方程和典范對(duì)應(yīng)分析#項(xiàng)等方法，定量探究地形因子與耕地分布的相關(guān)性。一般認(rèn)為，耕地分布與高程、坡度和地形起伏度等地形因子具有較顯著的相關(guān)性，與坡向相關(guān)性不顯著?？v觀現(xiàn)有研究，大多突出單個(gè)或多個(gè)地形因子對(duì)耕地分布的影響程度及其差異，不能很好地揭示不同地形因子之間的交互作用。地理探測(cè)器是近年來(lái)出現(xiàn)的一種用于探測(cè)空間分異性，以及揭示背后驅(qū)動(dòng)力的一組統(tǒng)計(jì)學(xué)方法口句，尤為突出的是對(duì)因子交互作用的探測(cè)功能，在地形因子影響耕地分布的研究中具有較大的應(yīng)用潛力。此外，現(xiàn)有研究中一般考慮高程、坡度、坡向和地形起伏度等地形因子較少考慮區(qū)域特征進(jìn)行地形因子選擇。對(duì)于千溝萬(wàn)壑、土壤侵蝕嚴(yán)重的黃土地區(qū)，溝壑密度是應(yīng)當(dāng)考慮的地形因子。該因子在一定程度上能夠表征溝谷侵蝕狀況，是地形發(fā)育階段、降水量、地勢(shì)高差、土壤滲透能力和地表抗侵蝕能力的綜合標(biāo)示值口一宵，反映著地表破碎程度、水土流失情況，進(jìn)而影響耕地分布數(shù)量和質(zhì)量⑵」（甘肅省榆中縣位于黃土高原，區(qū)內(nèi)溝壑縱橫，南北縱跨5個(gè)地貌類型區(qū)域，地形具有明顯過(guò)渡特征。以榆中縣為研究區(qū)，在地貌分區(qū)基礎(chǔ)上，利用分形指數(shù)和地理探測(cè)器等方法定量分析高程、坡度、坡向和溝壑密度對(duì)耕地分布影響，以及不同因子之間的交互作用,可為同類研究提供方法上的借鑒，也為研究區(qū)優(yōu)化耕地資源、改善耕地質(zhì)量提供一定的科學(xué)依據(jù)。1研究區(qū)榆中縣位于甘肅省中部，省會(huì)蘭州東郊,103°4&'E?104°35'E,35°33'N?36°27'N，面積為3301km$（如圖1所示，該地區(qū)在黃土高原與青藏高原交界處,地形由南向北依次可劃分為南部石質(zhì)山地、黃土塬區(qū)、黃土丘陵區(qū)、北部石質(zhì)低山丘陵區(qū)和河谷平原區(qū)南部石質(zhì)山地是祁連山東延余脈，有興隆山和馬銜山2條平行山嶺，呈西北）東南走向，雨水豐沛，植被茂盛，是全縣的生態(tài)屏障。黃土塬區(qū)由興隆山山前沖積、洪積平原和宛川河谷地組成，地勢(shì)平坦，交通便利，是全縣工農(nóng)業(yè)最發(fā)達(dá)地區(qū)。黃土丘陵區(qū)為中世紀(jì)馬蘭黃土覆蓋泥質(zhì)砂巖形成，干旱少雨,植被覆蓋度低，地表抗侵蝕能力差。北部石質(zhì)低山丘陵為基巖區(qū)，上覆厚度不等的黃土層。河谷平原區(qū)位于最北端，沿黃河呈西南）東北走向，河谷兩岸發(fā)育五級(jí)階地，居民點(diǎn)與農(nóng)業(yè)種植區(qū)位于二級(jí)階地上。區(qū)丘地山山區(qū)低區(qū)質(zhì)區(qū)陵質(zhì)原￥南黃黃北河注：該圖基于蘭州市自然資源局網(wǎng)站下載的審圖號(hào)為甘8（區(qū)丘地山山區(qū)低區(qū)質(zhì)區(qū)陵質(zhì)原￥南黃黃北河注：該圖基于蘭州市自然資源局網(wǎng)站下載的審圖號(hào)為甘8（2011）14號(hào)的標(biāo)準(zhǔn)地圖制作，底圖無(wú)修改區(qū)陵2數(shù)據(jù)與方法2.1數(shù)據(jù)來(lái)源研究采用30m分辨率的ASTERGDEMV2全球數(shù)字高程數(shù)據(jù)，下載自地理空間數(shù)據(jù)云。耕地?cái)?shù)據(jù)、地貌區(qū)劃數(shù)據(jù)和行政邊界數(shù)據(jù)來(lái)自寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心。其中，耕地?cái)?shù)據(jù)為2.5mSPOT-5遙感影像人機(jī)交互解譯成果，包括平原旱地、坡旱地、溝谷地、梯田和條田5個(gè)二級(jí)類。2.2數(shù)據(jù)處理地形因子的獲取在ArcGIS軟件中實(shí)現(xiàn)。坡度和坡向直接利用現(xiàn)有工具計(jì)算生成。溝壑密度定義為單位面積溝谷總長(zhǎng)度，采用樣方分析思想，通過(guò)水文分析工具和克里金模型插值計(jì)算生成溝壑密度。其提取過(guò)程是：首先，利用水文分析工具提取出溝谷網(wǎng)絡(luò)［17'25］;然后，生成方形網(wǎng)格并計(jì)算每個(gè)網(wǎng)格中的溝壑密度值，其中網(wǎng)格太小無(wú)法獲得想要的插值效果，網(wǎng)格太大又會(huì)降低數(shù)據(jù)精度，綜合考慮插值效果和精度問(wèn)題，經(jīng)多次實(shí)驗(yàn)將樣方面積控制在10km2左右，最終確定網(wǎng)格邊長(zhǎng)為3240m，即30m分辨率數(shù)據(jù)像元邊長(zhǎng)的10&倍；最后，分析插值點(diǎn)數(shù)據(jù)特征，數(shù)據(jù)服從正態(tài)分布并存在趨勢(shì)，采用泛克里金插

值法輸出30m分辨率溝壑密度數(shù)據(jù)，插值點(diǎn)數(shù)據(jù)統(tǒng)計(jì)特征如表1所示。表1插值點(diǎn)數(shù)據(jù)統(tǒng)計(jì)特征統(tǒng)計(jì)特征項(xiàng)數(shù)值統(tǒng)計(jì)特征項(xiàng)數(shù)值總數(shù)354.00000偏度=0.07008最小值0.10441峰度3.06570最大值2.091101/4分位數(shù)0.72080平均值0.96814中位數(shù)0.97710標(biāo)準(zhǔn)差0.331723/4分位數(shù)1.19940地理探測(cè)器模型要求自變量為類型量，對(duì)地形因子連續(xù)數(shù)據(jù)進(jìn)行離散化處理。如圖2所示，用300m間距將高程劃分為8個(gè)級(jí)別：1級(jí)(1358?TOC\o"1-5"\h\z600 m)、2 級(jí)(1 600—1 900 m)、3 級(jí)(1 900?200 m)、4 級(jí)(2 200-2 500 m)、5 級(jí)(2 500?2800 m)、6 級(jí)(2 800?3 100 m)、7 級(jí)(3 100?3400m)和8級(jí)(3400m)以上(1")。將坡度劃分為1級(jí)(0°?3°)、2級(jí)(3°?7°)、3級(jí)(7°?15°)、4級(jí)(15°?25°)、5級(jí)(25°?35°)和6級(jí)(〉35°)6個(gè)級(jí)別，即平坡、緩坡、斜坡、緩陡坡、陡坡和急陡坡[4'21]0坡向劃分為平地(一1)、北(0°?22.5°,337.5°?360°)、東北(22.5°?67.5°)、東(67.5°?112.5°)、東南(112.5°?157.5°)、南(157.5°?202.5°)、西南(202.5°?247.5°)、西(247.5°?292.5°)、西北(292.5°?337.5°)共9個(gè)等級(jí)。溝壑密度采取自然間斷點(diǎn)法，以(km/km2)為單位，劃分為1級(jí)(0?0.25)、2級(jí)(0.25?5)、3級(jí)(0.5?0.65)、4級(jí)(0.65?0.8)、5級(jí)(0.8?0.95)、6級(jí)(0.95?1.1)、7級(jí)(1.1?1.25)和8級(jí)(1.25?2.1)8個(gè)級(jí)別。使用ArcGIS漁網(wǎng)工具建立30mX30m網(wǎng)格，通過(guò)漁網(wǎng)點(diǎn)工具將各柵格像元值提取到點(diǎn)，以300m間隔進(jìn)行采樣獲得地理探測(cè)器運(yùn)行數(shù)據(jù)。高程分級(jí)/m1358-16001600-1900?1900-220002200-25002500-28002800-31003100-34003400?3671I坡度分級(jí)/(°)|?0~3|?3~7■7-15M15-25?25~35|?>35■東北溝壑密度分級(jí)/(km-km2)■0-0.25■0.25?0.500.50-0.650.65-0.800.80-0.950.95-1.10」1.10?1.25■1.25?2.10(a)高程(b)坡度 (c)坡向 (d)溝壑密度注：該圖基于蘭州市自然資源局網(wǎng)站下載的審圖號(hào)為甘S(2011)14號(hào)的標(biāo)準(zhǔn)地圖制作，底圖無(wú)修改。圖2榆中縣高程、坡度、坡向和溝壑密度分布圖2.3分維數(shù)和穩(wěn)定性指數(shù)分維數(shù)可以比較準(zhǔn)確地測(cè)定地物的復(fù)雜度和破碎程度(10)。董連科蛔根據(jù)Mandelbrot^27)的理論提出了適用于"維歐式空間關(guān)系的分形公式，如式(1)所示。TOC\o"1-5"\h\z-(*)1/C?—1=kr("-1-D，-1>/D，-1V(r)1/" (1)式中：，=2時(shí)，就是二維歐式空間的面積與周長(zhǎng)的分形公式，令S*)和P(r)分別代表以*為測(cè)量尺度的地物圖斑面積和周長(zhǎng)，得出表達(dá)式如式(2)所示。P(r)1/D=*〈1-D)/DS(r)1/2 (2)式中：兩側(cè)同時(shí)取以10為底的對(duì)數(shù)，得出式(3)。0A(r)=(2/D)lgP(r)+C (3)式中：D為所要求的分維值*A(r)為某地類圖斑面積*P(r)為該地類圖斑周長(zhǎng)C為待定常數(shù)。式(3)表示各地類圖斑的周長(zhǎng)面積之間的相互關(guān)系。對(duì)某一土地利用類型而言，將各斑塊的周長(zhǎng)一面積數(shù)據(jù)在0A(r)和0P(r)坐標(biāo)系中作回歸分析擬合，可得到該地類斑塊的周長(zhǎng)一面積擬合直線的斜率值(2/D)，進(jìn)而求得各地類的分維數(shù)D。D值介于1?2之間，它的值越大，表明地類圖斑的破碎程度越高。當(dāng)D=1.5時(shí)，空間結(jié)構(gòu)最活躍，即最不穩(wěn)定狀態(tài)(3'10]o由此，徐建華等(28)定義了穩(wěn)定性指數(shù)-如式(4)所示。-值越大，說(shuō)明耕地的空間形態(tài)結(jié)構(gòu)越穩(wěn)定。-=|D—1.5| (4)2.4地理探測(cè)器地理探測(cè)器是由王勁峰等(16)提出的，是用于探測(cè)空間分異性并挖掘背后驅(qū)動(dòng)力的一組統(tǒng)計(jì)學(xué)方法。其核心思想是基于這樣的假設(shè)：如果某個(gè)自變量對(duì)某個(gè)因變量有重要影響，那么自變量和因變量

的空間分布應(yīng)該具有相似性[29=30]'地理探測(cè)器分因子探測(cè)器、風(fēng)險(xiǎn)探測(cè)器、交互作用探測(cè)器、生態(tài)探測(cè)器4個(gè)模塊，本研究主要通過(guò)因子探測(cè)和交互作用探測(cè)分析地形因子對(duì)耕地分布的影響。因子探測(cè)：用于探測(cè)自變量X多大程度上解釋因變量丫的空間分異性，用g值度量[31],如式(5)、式(6)所示。(5)="Nz"=sswgx— N" —SST式中：Z=1，2)；L為自變量X的分層;Nz和N分別為層Z內(nèi)和區(qū)域內(nèi)的單元數(shù)"Z和"分別為層Z和全區(qū)的丫值的方差；SSW為層內(nèi)方差之和；SST為區(qū)域總方差gx為自變量對(duì)因變量的解釋力，&范圍0，1$，q值越接近于1，(5)交互作用探測(cè)：通過(guò)比較單因子g值與雙因子q值大小，判斷因子交互作用的方向及方式，是地理探測(cè)器相對(duì)于其他統(tǒng)計(jì)方法的最大優(yōu)勢(shì)。2個(gè)自變量對(duì)因變量交互作用類型如表2所示#16$。SSW—"Nz"z$，SST—N"z—1(6)表2自變量對(duì)因變量交互作用類型判據(jù)交互作用q(X1nX2)<Min(q(X1),q(X2))非線性減弱Min(q(X1)，q(X2))Vq(X1+X2)VMax(q(X1),q(X2))單因子非線性減弱q(X1+X2)>Max(q(X1),q(X2))雙因子增強(qiáng)q(X1nX2)=q(X1)+q(X2)獨(dú)立q(X1+X2)>q(X1)+q(X2)非線性增強(qiáng)地貌類型區(qū)河谷平原面積/km2墾殖率地貌類型區(qū)河谷平原面積/km2墾殖率/(%)面積/km2墾殖率/(%)面積/km2墾殖率/(%)面積/km2墾殖率/(%)面積/km2墾殖率/(%)8.141.290.000.0074.5311.790.000.0073.8611.69224.8663.060.000.0033.729.460.050.015.981.6825.811.712.060.14629.9841.670.880.0613.880.923.180.411.870.24144.9818.620.000.0057.367.3713.1959.820.000.000.000.000.000.001.466.63表3不同地貌類型區(qū)耕地面積統(tǒng)計(jì)耕地類型南部石質(zhì)山地黃土塬區(qū)黃土丘陵區(qū)北部石質(zhì)低山丘陵區(qū)平原旱地 溝谷地 梯田 條田 坡旱地3結(jié)果與分析3.1耕地面積分布情況研究區(qū)耕地總面積為1315.79km2,平原旱地面積為275.18km2,溝谷地面積為3.93km2,坡旱注：該圖基于蘭州市自然資源局網(wǎng)站下載的審圖號(hào)為甘S(2011)14號(hào)的標(biāo)準(zhǔn)地圖制作，底圖無(wú)修改。圖3榆中縣不同地貌區(qū)耕地分布情況山區(qū)地山區(qū)低區(qū)地質(zhì)區(qū)陵質(zhì)原地地旱石墟丘石平田田谷旱原部土土部谷條梯溝坡平南黃黃北河區(qū)陵丘3.2耕地分形特征1)不同地貌類型區(qū)耕地分形特征。不同地貌區(qū)耕地分維數(shù)與穩(wěn)定性指數(shù)具有顯著特征(表4)。分維數(shù)表示耕地斑塊復(fù)雜度與破碎程度，穩(wěn)定性指數(shù)值反映耕地在自然狀態(tài)下抵御外界干擾維持自身形態(tài)的能力從耕地的分維數(shù)來(lái)看，地勢(shì)平坦的黃土塬區(qū)與河谷平原區(qū)明顯小于其他地貌區(qū)，分別為1.470和1.230,表明這2個(gè)地貌區(qū)內(nèi)耕地斑塊形狀規(guī)則，連片地面積為152.54km2,梯田面積為883.21注：該圖基于蘭州市自然資源局網(wǎng)站下載的審圖號(hào)為甘S(2011)14號(hào)的標(biāo)準(zhǔn)地圖制作，底圖無(wú)修改。圖3榆中縣不同地貌區(qū)耕地分布情況山區(qū)地山區(qū)低區(qū)地質(zhì)區(qū)陵質(zhì)原地地旱石墟丘石平田田谷旱原部土土部谷條梯溝坡平南黃黃北河區(qū)陵丘3.2耕地分形特征1)不同地貌類型區(qū)耕地分形特征。不同地貌區(qū)耕地分維數(shù)與穩(wěn)定性指數(shù)具有顯著特征(表4)。分維數(shù)表示耕地斑塊復(fù)雜度與破碎程度，穩(wěn)定性指數(shù)值反映耕地在自然狀態(tài)下抵御外界干擾維持自身形態(tài)的能力從耕地的分維數(shù)來(lái)看，地勢(shì)平坦的黃土塬區(qū)與河谷平原區(qū)明顯小于其他地貌區(qū)，分別為1.470和1.230,表明這2個(gè)地貌區(qū)內(nèi)耕地斑塊形狀規(guī)則，連片表4不同地貌類型區(qū)耕地分維數(shù)和穩(wěn)定性指數(shù)地貌分區(qū)分維數(shù)（D）穩(wěn)定性指數(shù)（S）南部石質(zhì)山地1.5280.028黃土塬區(qū)1.4700.030黃土丘陵區(qū)1.5810.081北部石質(zhì)低山丘陵區(qū)1.6340.134河谷平原區(qū)1.2300.270石質(zhì)山地（1.528）,表明這些分區(qū)的耕地分布相對(duì)破碎，黃土丘陵與北部石質(zhì)低山丘陵溝谷分布相對(duì)密集，高程與坡度范圍跨度較大，南部石質(zhì)山地地勢(shì)最高，坡度以陡坡為主，這些相對(duì)復(fù)雜的地表形態(tài)使得耕地分維數(shù)較高"32#。從耕地的穩(wěn)定性指數(shù)來(lái)看,河谷平原最大（0.270）,空間分布最穩(wěn)定；黃土丘陵與北部石質(zhì)低山丘陵區(qū)次之（分別為0.081和0.134）,空間分布穩(wěn)定性居中；黃土塬與南部石質(zhì)山地取值最?。ǚ謩e為0.030和0.028）,穩(wěn)定性最差。河谷平原區(qū)、黃土丘陵與北部石質(zhì)低山丘陵區(qū)人類活動(dòng)主要以農(nóng)業(yè)種植為主,耕地類型不容易發(fā)生改變，穩(wěn)定良好黃土塬地區(qū)工業(yè)生產(chǎn)較發(fā)達(dá)，耕地穩(wěn)定性低，受人類干擾程度較大南部石質(zhì)山地耕地分維數(shù)較大，穩(wěn)定性低「32#。2）不同耕地在各個(gè)地貌類型區(qū)的分形特征。耕地開墾順應(yīng)地形特征，平原旱地分布在地勢(shì)平坦開闊處，坡旱地分布在具有一定坡度的地表上，為保持水土，提高農(nóng)作物產(chǎn)量，對(duì)自然坡面進(jìn)行改造修整,開墾梯田、條田，溝谷地一般分布在黃土地區(qū)丘陵的溝谷間。表5顯示了不同地貌區(qū)各耕地分位數(shù)與穩(wěn)定性指數(shù)。可以看出，各個(gè)地貌區(qū)平原旱地分維數(shù)一般小于其他類耕地，分布于坡面上的梯田、坡旱地分維數(shù)較大，說(shuō)明平地處的耕地比坡面上的耕地破碎度低，連片性好。原因?yàn)辄S土地區(qū)坡面上細(xì)溝、淺溝侵蝕強(qiáng)烈，斜坡下部被沖溝、切溝所分割，坡面破碎，使得耕地分布破碎度較大。平原旱地在河谷平原區(qū)穩(wěn)定性指數(shù)最大為0.305；梯田在黃土丘陵區(qū)與北部石質(zhì)低山丘陵具有顯著分布優(yōu)勢(shì)，穩(wěn)定性指數(shù)較大,分別為0.218和0.185；坡旱地在南部石質(zhì)山地相比其他地貌區(qū)穩(wěn)定性指數(shù)較大，為0.145；所有類型耕地在黃土塬區(qū)穩(wěn)定性指數(shù)都小于0.1。由此看出，各類耕地穩(wěn)定性指數(shù)受地形影響存在差異。榆中縣溝谷地與條田分布極少，所得出分維數(shù)與穩(wěn)定性指數(shù)代表性不大。表5不同耕地在各個(gè)地貌類型區(qū)分維數(shù)和穩(wěn)定性指數(shù)地貌類型區(qū)耕地種類分維數(shù)（C）穩(wěn)定性指數(shù)―）平原旱地1.3520.148溝谷地——南部石質(zhì)山地梯田1.5160.016條田—坡旱地1.6450.145平原旱地1.4330.067溝谷地—黃土塬區(qū)梯田1.5940.094條田——坡旱地1.4990.001平原旱地1.4530.047溝谷地1.3560.144黃土丘陵區(qū)梯田1.7180.218條田1.1250.375坡旱地1.4490.051平原旱地1.3760.124溝谷地1.5120.012北部石質(zhì)低山丘陵區(qū)梯田1.6850.185條田——坡旱地1.6090.109平原旱地1.1950.305溝谷地——河谷平原區(qū)梯田——條田—坡旱地1.6580.158注:-表示地貌區(qū)內(nèi)不存在此類耕地。3.3地形對(duì)耕地分布的影響由分維數(shù)與穩(wěn)定性指數(shù)可以看出，不同地貌區(qū)耕地分布特征不同，地形因子從不同角度定量表達(dá)了地貌形態(tài)特征，選擇能夠表征榆中縣地貌形態(tài)特征的高程、坡度、坡向和溝壑密度，通過(guò)地理探測(cè)器定量分析這4個(gè)地形因子對(duì)耕地分布的影響大小。1）不同地形因子對(duì)耕地分布的影響。表6顯示了因子探測(cè)結(jié)果，可以看出，不同地貌區(qū)各地形因子對(duì)耕地分布的影響大小具有顯著差異。高程對(duì)耕地分布影響力在南部石質(zhì)山地和黃土丘陵區(qū)均排第一,分別為0.177和0.132；在北部石質(zhì)低山丘陵區(qū)影響力高達(dá)0.255,但略低于溝壑密度，排第二。耕地分布隨高程變化具有一定的梯度性，這3個(gè)地貌區(qū)域內(nèi)高程范圍跨度較大，所以高程對(duì)耕地分布影響力顯著坡度在地勢(shì)平坦的河谷平原區(qū)與黃土塬區(qū)對(duì)耕地分布影響力顯著,分別為0.101和0.127,在其他地貌區(qū)影響力較小。王歡等曾指出坡度在地形相對(duì)平坦的地區(qū)對(duì)土壤侵蝕影響作用顯著，大于地形相對(duì)陡峭地區(qū)，這與本文坡度在地勢(shì)平坦地區(qū)對(duì)耕地分布影響力顯著具有相似性?？赡茉蚴窍鄬?duì)陡峭地區(qū)地形復(fù)雜，其他地形因子影響作用更大，因此，坡度對(duì)耕地分布影響力減弱。坡向在各個(gè)地貌區(qū)對(duì)耕地分布影響力都小于其他地形因子。溝壑密度對(duì)耕地分布影響力在北部石質(zhì)低山丘陵區(qū)排第一，為0.27&；在黃土丘陵區(qū)排第三，為0.027；其他地貌區(qū)均排在第二。研究區(qū)溝谷縱橫，耕地分布受溝壑密度影響總體較大。高程在河谷平原區(qū)只有1級(jí)（1358-1600m）,模型無(wú)法計(jì)算河谷平原區(qū)高程對(duì)耕地分布影響大小。表6不同地貌類型區(qū)地形因子對(duì)耕地分布影響q值統(tǒng)計(jì)地貌類型區(qū)高程坡度坡向溝壑密度南部石質(zhì)山地0.1770.0390.0290.0H8黃土塬區(qū)0.0410.12H0.0220.093黃土丘陵區(qū)0.1320.0450.00H0.02H北部石質(zhì)低山丘陵區(qū)0.2550.0380.00H0.2H8河谷平原—0.1010.0530.102注:-表示g值未經(jīng)模型計(jì)算。2）地形因子交互作用對(duì)耕地分布影響。耕地分布同時(shí)受多個(gè)地形因子的影響。分析得知，地形因子兩兩交互作用對(duì)耕地分布影響作用均增強(qiáng)。將地形因子交互作用影響力排在前三的因子組合及影響力g值整理得到表7。南部石質(zhì)山地、黃土塬區(qū)對(duì)耕地分布影響最大的地形因子分別為高程與坡度,與溝壑密度交互作用后影響力顯著增強(qiáng)，分別為0.294和0.207。黃土丘陵區(qū)影響力最顯著的地形因子組合為高程和坡度，高程與溝壑密度、坡度與溝壑密度交互影響力略低，分別為0.183和0.153。北部石質(zhì)低山丘陵區(qū)對(duì)耕地分布影響作用排在前三的地形因子組合