王莉娜，李文龍，王素芳，陳 迪,許 靜

(1.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，甘肅 蘭州 730020； 2.甘肅省環(huán)境監(jiān)測(cè)中心站，甘肅 蘭州 730020；

3.甘南州國土資源局，甘肅 甘南 747000； 4.蘭州財(cái)經(jīng)大學(xué)農(nóng)林經(jīng)濟(jì)管理學(xué)院，甘肅 蘭州 730020)

?

基于遙感和USLE模型的2000-2010年甘肅省土壤侵蝕變化評(píng)價(jià)

王莉娜1,2，李文龍1，王素芳3，陳 迪1,許 靜4

(1.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，甘肅 蘭州 730020； 2.甘肅省環(huán)境監(jiān)測(cè)中心站，甘肅 蘭州 730020；

3.甘南州國土資源局，甘肅 甘南 747000； 4.蘭州財(cái)經(jīng)大學(xué)農(nóng)林經(jīng)濟(jì)管理學(xué)院，甘肅 蘭州 730020)

摘要：本研究將通用土壤流失方程(USLE)和遙感、地理信息技術(shù)相結(jié)合，定量評(píng)價(jià)甘肅省2000-2010年土壤侵蝕狀況。研究結(jié)果表明，甘肅省主要侵蝕強(qiáng)度為輕度和中度侵蝕，2000、2005和2010年甘肅省土壤侵蝕面積分別為24 872.9、27 503.6和28 000.8 km2。主要侵蝕地區(qū)分布在甘肅省南部地區(qū)；林地的土壤侵蝕面積最大，草地次之，這兩類土地利用類型為主要侵蝕地帶。隴南山地、甘南高原和隴中黃土高原為土壤侵蝕的主要地貌類型。2000-2010年土壤侵蝕變化分析表明，造成土壤侵蝕的原因除了降水量分布不均、土地利用類型和地貌類型的變化外，人類對(duì)生態(tài)環(huán)境的不合理利用也是另一重要原因。

關(guān)鍵詞：土壤侵蝕；土壤流失方程(USLE)；遙感

土壤侵蝕是地球表面的一種自然現(xiàn)象，是指地球表面的土壤及其母質(zhì)在水力、風(fēng)力、凍融、重力等外力的作用及人為因素影響下發(fā)生的各種破壞、分離、搬運(yùn)和沉積的現(xiàn)象[1]。地表物質(zhì)在外在驅(qū)動(dòng)力的作用下，其自身產(chǎn)生位移與再分配，其內(nèi)在的營養(yǎng)成分也會(huì)流失于水體之中[2]，這不僅導(dǎo)致土地的生產(chǎn)能力下降，也會(huì)使水體產(chǎn)生富營養(yǎng)化。

土壤侵蝕預(yù)報(bào)是監(jiān)測(cè)水土流失和評(píng)價(jià)水保措施效益的有效手段，侵蝕模型則是進(jìn)行土壤流失監(jiān)測(cè)和預(yù)報(bào)的重要工具[3-4]。通用土壤流失方程(Universal Soil Loss Equation,USLE)是美國研制的用于定量預(yù)報(bào)農(nóng)田或草地坡面多年平均土壤流失量的一個(gè)經(jīng)驗(yàn)性土壤侵蝕預(yù)報(bào)模型，主要研究由降水引起的水動(dòng)力土壤侵蝕[5]。從20世紀(jì)30年代開始，通用土壤流失方程經(jīng)過大量的土壤侵蝕試驗(yàn)和定量研究得到了不斷的發(fā)展完善。1940年Zingg[6]通過研究坡長、坡度對(duì)土壤侵蝕的影響，土壤侵蝕預(yù)報(bào)模型被正式應(yīng)用于土壤侵蝕的研究。較之其它基于物理過程的模型(如水利侵蝕預(yù)報(bào)模型WEPP和動(dòng)力侵蝕模型KINEROS)，USLE結(jié)構(gòu)簡單，所需數(shù)據(jù)量少，結(jié)果可靠性高，在世界范圍內(nèi)得到了廣泛應(yīng)用[7-8]，在定量評(píng)價(jià)土壤侵蝕強(qiáng)度、土地資源合理利用和水土保持規(guī)劃等方面起到了重要的作用[9-10]。我國從20世紀(jì)80年代以來，開始引進(jìn)通用土壤流失方程，在我國土壤侵蝕預(yù)報(bào)工作中起到了積極的作用，并取得了重要成果[11-14]。

本研究采用通用土壤流失方程，結(jié)合遙感、地理信息技術(shù)，定量評(píng)價(jià)甘肅省2000-2010年土壤侵蝕狀況，分析近10年甘肅省的土壤侵蝕變化，進(jìn)而為甘肅省生態(tài)環(huán)境質(zhì)量狀況評(píng)價(jià)提供科學(xué)依據(jù)。

1材料與方法

1.1研究區(qū)概況

甘肅位于中國三大高原交匯地帶，地貌復(fù)雜多樣，河谷、高原、山地、戈壁與沙漠交錯(cuò)分布，是典型的山地型高原地貌，全省呈東西長，南北短的空間格局，而黃土高原以丘陵溝壑和高原溝壑地貌為主。不同地區(qū)的氣候類型差異顯著，主要包括干旱區(qū)、高寒區(qū)與亞熱帶濕潤區(qū)。氣候干燥少雨，氣溫日差大，年平均氣溫在0～14 ℃，年均降水量200~700 mm，降水主要集中在7-9月。

1.2研究數(shù)據(jù)

本研究利用的主要數(shù)據(jù)包括：1)甘肅省空間分辨率為90 m的數(shù)字高程模型(DEM)；2)基于Landsat5 TM影像解譯的甘肅省2000、2005和2010年3期土地利用類型/土地覆蓋圖；3)甘肅省土壤類型圖；4)甘肅省2000、2005和2010年81個(gè)氣象站點(diǎn)的降水量觀測(cè)數(shù)據(jù)；5)2000、2005和2010年MODIS月植被指數(shù)產(chǎn)品；6)甘肅省土壤質(zhì)地、有機(jī)質(zhì)含量等土壤屬性數(shù)據(jù)。

1.3USLE模型及其因子的確定方法

本研究利用RS和GIS技術(shù)，并結(jié)合USLE模型來對(duì)土壤侵蝕量進(jìn)行預(yù)測(cè)。USLE模型如下：

A=R×K×L×S×C×P

(1)

式中，A為年平均土壤流失量(t·hm-2·a-1)；R為降雨侵蝕力因子(MJ·mm·hm-2·h-1·a-1)；S為坡度因子，無量綱；P為水土保持措施因子，無量綱；L為坡長因子，無量綱；K為土壤可蝕性因子(t·h·MJ-1·mm-1)；C為植被覆蓋度因子，無量綱。

1.3.1降水侵蝕力因子的估算降水侵蝕力(R)是指降水過程對(duì)土壤流失存在的潛在影響力。甘肅省降水空間分布差異顯著，因此，降水對(duì)全省各地水土流失的影響也不同。國內(nèi)外學(xué)者對(duì)降水侵蝕力因子的計(jì)算方法做了大量的研究，本研究采用Wischmeier和Smith[5]提出的基于月平均降水量和年降水量的經(jīng)驗(yàn)公式：

(2)

式中，Pi為月降水量，P為年降水量，單位均為mm。

降水量數(shù)據(jù)來自甘肅省81個(gè)縣市氣象站點(diǎn)的觀測(cè)資料。通過計(jì)算月平均降水與年降水，然后采用Wischmeier和Smith[5]提出的經(jīng)驗(yàn)公式來計(jì)算2000、2005和2010年各縣市的降雨侵蝕力。

1.3.2土壤可蝕性因子的估算土壤可侵蝕性因子(K)大小是指標(biāo)準(zhǔn)條件下，土壤與土壤剖面對(duì)各種侵蝕與水動(dòng)力過程的平均敏感程度。本研究采用Wischmeier等[15]的侵蝕力影響模型 (ErosionProductivityImpactCalculator，EPIC)估算土壤可蝕性因子K值。依據(jù)土壤顆粒組成與有機(jī)質(zhì)含量等資料推算研究區(qū)內(nèi)不同土地利用類型的可蝕性因子，計(jì)算公式如下：

K={0.2+0.3e[-0.025 6S1 (1.0-S2/100)]}×[S2/(n+S2)]0.3×[1.0-0.25C/(C+e(3.72-2.95C)]×[1.0-0.7S3/(S3+e(-5.51+22.9S3 )]

(3)

式中，S1為土壤砂粒含量(%)，S2為土壤粉粒含量(%)，n為土壤粘粒含量(%)，C為土壤有機(jī)碳含量(%)，S3=1-S1/100。

利用在全省采集的2 900個(gè)樣品，分析土壤砂粒、粉粒、粘粒和有機(jī)質(zhì)含量，并采用VanBemmelen因數(shù)1.724[16]，計(jì)算有機(jī)碳含量，利用土壤可蝕性K值計(jì)算公式，近似確定出研究區(qū)不同土壤類型的可蝕性因子K值。

1.3.3植被覆蓋因子的估算植被蓋度(C)是植被在垂直投影面占總面積的百分比，本研究利用像元二分模型計(jì)算研究區(qū)的植被覆蓋度作為植被因子度量指標(biāo)，計(jì)算公式如下：

Fg=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)

(4)

式中，F(xiàn)g為植被蓋度，NDVIsoil為純土壤像元的NDVI值，即在NDVI最小值概率分布表中最接近0.5%位置處的NDVI值；NDVIveg為純植被像元的NDVI值，即在NDVI最大值概率分布表中最接近0.5%位置處的NDVI值[17]。利用MODIS植被歸一化差值植被指數(shù)(NDVI)計(jì)算甘肅省2000、2005和2010年植被覆蓋度。

1.3.4坡長和坡度因子的估算在地形因素中，反映土壤坡面侵蝕的坡長(L因子)和坡度(S因子)是兩個(gè)主要因子。坡長增加，坡面來水面積增加，向坡下方流動(dòng)的徑流量、流速、水力半徑也隨著增大，侵蝕量也隨之增大。坡長L因子計(jì)算公式是由Liu等[18]提出的，具體公式為：

L=(λ/22.13)α

(5)

α=β/(β+1)

(6)

β=(sinθ/0.089 6)/[3.0(sinθ)0.8+0.56]

(7)

式中，L為坡長因子，λ為水平坡長，α為坡長指數(shù)，θ為坡度。

坡度S因子采用的是McCool等人在1987年提出并由國內(nèi)學(xué)者Liu等[19]于1994年改進(jìn)的計(jì)算公式：

(8)

式中，S為坡度因子，θ為坡度。

1.3.5水土保持措施因子的估算水土保持措施因子[20](P)指為抵制水土侵蝕，而采取的保持水土的措施。根據(jù)土地利用類型不同，P的取值范圍在0-1之間，采取水土保持措施越少則P值越大。本研究參照劉寶元等[21]相關(guān)研究中P值的確定方法，得出了甘肅省林地和未利用地的水土保持措施因子為1.0，其次草地和濕地為0.7，耕地為0.3，城鄉(xiāng)居民點(diǎn)和工礦用地及水域?yàn)?。利用甘肅省2000、2005、2010年土地利用/覆蓋圖進(jìn)行全省水土保持因子(P)的確定。

1.3.6土壤侵蝕量的估算利用ArcGIS軟件對(duì)侵蝕因子圖層進(jìn)行疊加分析，以此得到各區(qū)域的年土壤侵蝕量。通過參考水利部頒發(fā)的土壤侵蝕強(qiáng)度分類分級(jí)標(biāo)準(zhǔn)(SL-190-96)，并且結(jié)合甘肅省的實(shí)際情況，如表1標(biāo)準(zhǔn)進(jìn)行全省土壤侵蝕量的等級(jí)劃分，得到甘肅省2000、2005和2010年土壤侵蝕強(qiáng)度分布圖(圖1)。

2結(jié)果與分析

2.1甘肅省2000-2010年土壤侵蝕及其變化

通過分析甘肅省土壤侵蝕強(qiáng)度分布狀況得出：甘肅省土壤侵蝕主要發(fā)生在甘南、隴南和隴中等地區(qū)，另外在甘肅省中部的天祝藏族自治縣也有分布，強(qiáng)度主要為輕度和中度侵蝕(圖1)。根據(jù)甘肅省土壤侵蝕強(qiáng)度統(tǒng)計(jì)結(jié)果(表2)可知，2000、2005和2010年甘肅省土壤侵蝕(不含無明顯侵蝕類型)面積分別為24 873.0、27 503.6和28 000.7 km2，分別占土地總面積的6.16%、6.81%和6.52%。從輕度侵蝕到重度侵蝕，侵蝕面積隨著強(qiáng)度等級(jí)的增加而減少。在2000-2010年期間重度侵蝕面積分別占侵蝕總面積的20.32%、23.22%和20.86%，中度侵蝕面積分別占侵蝕總面積的24.81%、25.94%和25.74%，輕度侵蝕面積分別占侵蝕總面積的54.87%、50.84%和53.39%。說明甘肅省的土壤侵蝕強(qiáng)度以輕度和中度侵蝕為主。

表1　甘肅省土壤侵蝕強(qiáng)度分級(jí)表

圖1　2000、2005和2010年甘肅省土壤侵蝕強(qiáng)度分級(jí)圖

表2　2000, 2005和2010年甘肅省不同土壤侵蝕強(qiáng)度面積及比例

2.2甘肅省土壤侵蝕與土地利用類型的關(guān)系

分析不同土地利用類型對(duì)土壤侵蝕的影響，將土壤侵蝕圖層與2000、2005和2010年3期土地利用/覆蓋圖層進(jìn)行疊加統(tǒng)計(jì)分析，得出甘肅省各土地利用類型的不同侵蝕強(qiáng)度的土壤侵蝕面積。統(tǒng)計(jì)結(jié)果表明，2000、2005和2010年不同土地利用類型的侵蝕面積存在很大差異性(表3)，其中林地侵蝕總面積最大，分別為15 244.6、16 111.0和16 324.9 km2，分別占該土地利用類型總面積的26.83%、28.15%和29.86%，10年間侵蝕面積呈增加趨勢(shì)。總體來看，土壤侵蝕面積從大到小依次為林地>草地>耕地>未利用地。在3個(gè)時(shí)期中，草地土壤侵蝕總面積分別為7 670.9、8 830.5和9 042.4 km2，占該土地利用類型總面積的7.06%、8.08%和7.96%；耕地土壤侵蝕總面積分別為1 567.1、2 199.6和2 215.7 km2，占該土地利用類型總面積的2.42%、3.32%和3.18%。在3個(gè)時(shí)期中，草地和耕地的土壤侵蝕面積逐漸增加。未利用地的土壤侵蝕面積在3個(gè)時(shí)期中分別為327.9、294.1和417.7 km2，占該土地利用類型總面積的0.19%、0.19%和0.23%，2010年面積大幅度增加，未利用地土壤侵蝕面積所占比例先減小后增加，這是由于2010年對(duì)甘肅省行政邊界進(jìn)行了修正，各類土地總面積均有所增加。

林地土壤侵蝕最嚴(yán)重，一方面由于林地中疏林地的植被較稀疏，本身水土保持能力比較差；另一方面林地多分布在甘肅省南部，該地區(qū)降水量多，加上長期的濫砍濫伐等不合理利用，導(dǎo)致林地的侵蝕加劇。草地土壤侵蝕多分布在甘南地區(qū)，該地區(qū)農(nóng)牧民人口日益增加，為了滿足生活需求，牲畜數(shù)量也大幅增加，超載過牧引起草地退化日益加劇，導(dǎo)致草地的水土保持功能減弱。甘肅省的耕地以旱地為主，發(fā)生土壤侵蝕的耕地大多分布在隴東黃土高原地區(qū)，該地區(qū)耕地地表比較破碎，降水多集中在夏季汛期，容易發(fā)生水土流失，其它大部分耕地多分布在降水量較少的地區(qū)，故侵蝕面積所占比例較小。需要關(guān)注的是，隨著人們對(duì)農(nóng)牧產(chǎn)品需求的不斷增加，濫砍濫伐、毀林開荒、超載過牧、濫采濫挖等行為將會(huì)導(dǎo)致土壤侵蝕狀況進(jìn)一步惡化。要科學(xué)合理地利用林地和草地，并對(duì)其進(jìn)行保護(hù)，以防其發(fā)生退化；而對(duì)未利用地、疏林地等植被覆蓋度較低的地類進(jìn)行生態(tài)恢復(fù)，增強(qiáng)其水土保持功能。

2.3甘肅不同地貌類型與土壤侵蝕的關(guān)系

甘肅省共有7類主要的地貌類型，分別為北方山地、河西走廊、阿爾金山和阿拉善、祁連山地、隴中黃土高原、甘南高原和隴南山地(圖2)。分布在最北部的北方山地海拔1 000～3 600 m，靠近騰格里沙漠和巴丹吉林沙漠，多為荒漠，氣候干旱少雨。河西走廊海拔1 400～2 900 m，地勢(shì)平坦，光熱充足，水資源較豐富。阿爾金山橫跨青海、新疆與甘肅三省，平均海拔4 000m以上，由系列山脈和谷地組成。阿拉善海拔處于800～1 600 m，地勢(shì)向北傾斜，區(qū)域內(nèi)分布眾多丘陵、剝蝕低山，其中山間為大小不一的盆地，地表主要由流沙與沙礫覆蓋。祁連山區(qū)域大部分海拔在4 000 m以上，內(nèi)部分布著多個(gè)平行山脈，且山脈間由一系列谷地與盆地組成，植被垂直地帶性分布明顯。隴中黃土高原海拔1 400～2 200 m，受流水侵蝕作用，地表破碎。甘南高原位于甘肅省西南部，是青藏高原的一部分，海拔處于3 000～4 200 m，區(qū)域內(nèi)呈波狀起伏狀，高原與山地相間分布，主要土地利用類型為草地。隴南山地位于甘肅東南部，為秦嶺山脈的西延部分，地勢(shì)東低西高，海拔從東部的1 500 m上升到西部的4 500 m以上。

表3　甘肅省不同土地利用類型土壤侵蝕面積及比例

圖2　甘肅省主要地貌類型圖

分析不同地貌類型對(duì)土壤侵蝕的影響，將土壤侵蝕圖層與甘肅省地貌類型圖層做疊置分析，統(tǒng)計(jì)甘肅省各類地貌類型的不同侵蝕強(qiáng)度的土壤侵蝕面積(表4)?？梢钥闯?，北方山地以及阿爾金山和阿拉善幾乎沒有土壤侵蝕發(fā)生，分析其原因主要是該地區(qū)降水量稀少，地貌類型決定了其大部分地區(qū)的坡度均小于5°，另一方面是USLE模式存在一定的缺陷性，在評(píng)價(jià)過程中沒有考慮風(fēng)蝕的影響；在河西走廊南部地區(qū)有少量輕度土壤侵蝕發(fā)生，這些地區(qū)除降水量較多外，坡度一般大于15°。祁連山地的土壤侵蝕主要發(fā)生在南部地區(qū)(天祝藏族自治縣境內(nèi))，2000、2005和2010年土壤侵蝕面積分別為1 152.5、661.3和682.8 km2(表4)。發(fā)生土壤侵蝕的地區(qū)一般雨水充沛，坡度多大于25°。甘肅省土壤侵蝕面積大的地貌類型為隴南山地和甘南高原，隴南山地東部地區(qū)土壤侵蝕強(qiáng)度較大，甘南高原東北部和東南部土壤侵蝕強(qiáng)度較大。2000、2005和2010年隴南山地土壤侵蝕面積分別為10 504.4、10 899.5和10 033.8 km2；甘南高原土壤侵蝕面積分別為10 212.4、11 432.9和11 659.8 km2(表4)。這兩種地貌類型侵蝕面積相差不大，侵蝕比例均在25%以上。分析隴南山地引起土壤侵蝕的主要原因發(fā)現(xiàn)，除了降水量較多、山地坡度大等自然因素外，人類對(duì)森林的大量砍伐，過度利用也是導(dǎo)致土壤侵蝕的加劇的重要因素，甘南高原也存在與隴南山地類似的狀況。隴中黃土高原的土壤侵蝕主要發(fā)生在東部、中部和西部的小部分地區(qū)，2000、2005和2010年土壤侵蝕面積分別為2 667.2、4 311.2和5 392.6 km2，土壤侵蝕比例在5.49%以下。黃土高原雖然地表比較破碎，但是大部分地區(qū)降水量很少，不容易造成土壤流失，而降水量較大的地區(qū)，土壤侵蝕較嚴(yán)重。同時(shí)發(fā)現(xiàn)，同期的降水量呈持續(xù)增加趨勢(shì)，是導(dǎo)致該區(qū)土壤侵蝕逐年增加的主要因素。

表4　甘肅省不同地貌類型土壤侵蝕面積及比例

3結(jié)論

目前，甘肅省土壤侵蝕調(diào)查方法以野外現(xiàn)地實(shí)測(cè)與訪問為主，本研究通過引入遙感和USLE模型對(duì)估算甘肅省的土壤侵蝕方法具有一定的改進(jìn)作用。利用空間分析技術(shù)可以很好地克服傳統(tǒng)方法的局限性，傳統(tǒng)的野外調(diào)查及長期觀測(cè)不僅費(fèi)時(shí)費(fèi)力，而且觀測(cè)區(qū)域小，周期長，但隨著RS和GIS技術(shù)的不斷發(fā)展，不僅可以實(shí)時(shí)掌握水土流失狀況，還能對(duì)區(qū)域水土流失做出預(yù)測(cè)和合理規(guī)劃。本研究利用RS和GIS技術(shù)，結(jié)合USLE模型，對(duì)甘肅省土壤侵蝕量的時(shí)空變化進(jìn)行了估算和分析。主要結(jié)論如下：

1)甘肅省土壤侵蝕狀況嚴(yán)重，2000、2005和2010年甘肅省土壤侵蝕面積分別為24 872.9、27 503.6和28 000.8 km2，侵蝕強(qiáng)度主要以輕度和中度侵蝕為主。

2)甘肅省不同地區(qū)土壤侵蝕空間差異性大，主要侵蝕地區(qū)分布在甘肅省南部，降水量的分布不均、土地利用類型和地貌類型的變化是造成這種差異的主要原因。林地的土壤侵蝕面積最大，草地次之，這兩種土地利用類型為甘肅省主要侵蝕地帶。隴南山地、甘南高原和隴中黃土高原為土壤侵蝕的主要地貌類型。

3)造成土壤侵蝕的原因除了降水量分布不均、土地利用類型和地貌類型的變化外，人類活動(dòng)對(duì)生態(tài)環(huán)境的不合理利用也可能是另一重要原因。

4)USLE模式存在一定的缺陷性，在評(píng)價(jià)過程中沒有考慮風(fēng)蝕的影響，評(píng)價(jià)結(jié)果與實(shí)際侵蝕情況有一定的誤差，在今后的研究過程中應(yīng)加以改進(jìn)。

參考文獻(xiàn)(References)

[1]楊勤科,李銳,劉詠梅.區(qū)域土壤侵蝕普查方法的初步討論.中國水土保持科學(xué),2008,6(3):1-7.

Yang Q K,Li R,Liu Y M.Discussions on methods of national regional soil erosion survey.Science of Soil and Water Conservation,2008,6(3):1-7.(in Chinese)

[2]徐炳成,山侖.半干旱黃土丘陵區(qū)白羊草人工草地生產(chǎn)力和土壤水分特征研究.草業(yè)科學(xué),2004,21(6):6-10.

Xu B C,Shan L.Aboveground productivity and soil water characteristics of cultivated Old World bluestem grassland in semi-arid loess hilly-gully areas on the Loess Plateau.Pratacultural Science,2004,21(6):6-10.(in Chinese)

[3]張巖,袁建平,劉元寶.土壤侵蝕預(yù)報(bào)模型中的植被覆蓋與管理因子研究進(jìn)展.應(yīng)用生態(tài)學(xué)報(bào),2002,13(8):1033-1036.

Zhang Y,Yuan J P,Liu B Y.Advance in researches on vegetation cover and management factor in the soil erosion prediction model.Chinese Journal of Applied Ecology,2002,13(8):1033-1036.(in Chinese)

[4]張科利,劉寶元,蔡永明.土壤侵蝕預(yù)報(bào)研究中的標(biāo)準(zhǔn)小區(qū)問題論證.地理研究,2000,19(3):297-302.

Zhang K L,Liu B Y,Cai Y M.The standard of unit plot in soil loss prediction of China.Geographical Research,2000,19(3):297-302.(in Chinese)

[5]Wischmeier W H,Smith D D.Predicting Rainfall Erosion Losses From Cropland East of the Rocky Mountains.Agriculture Handbook,1965,No.282.

[6]Zingg A W.Degree and length of land slope as it affects soil loss in run-off.Agrucyktyrak Engineering,1940,21:59-64.

[7]黃金良,洪華生,張珞平,杜鵬飛.基于GIS和USLE的九龍江流域土壤侵蝕量預(yù)測(cè)研究.水土保持學(xué)報(bào),2005,18(5):75-79.

Huang J L,Hong H S,Zhang L P,Du P F.Study on predicting soil erosion in Jiulong River watershed based on GIS and USLE.Journal of Soil and Water Conservation,2005,18(5):75-79.(in Chinese)

[8]鄭淑華,王堃,趙萌莉,韓國棟,馮雨峰.北方農(nóng)牧交錯(cuò)區(qū)草地生態(tài)系統(tǒng)服務(wù)間接價(jià)值的初步評(píng)估——以太仆寺旗和沽源縣境內(nèi)為例.草業(yè)科學(xué),2009,26(9):18-23.

Zheng S H,Wang K,Zhao M L,Han G D,Feng Y F.Primary evaluation of the indirect value on rangeland ecosystem services in Northern agro pastoral ecotone-A case study in Taipusi banner and Guyuan league.Pratacultural Science,2009,26(9):18-23.(in Chinese)

[9]曹麗花,劉合滿,趙世偉.退化高寒草甸土壤有機(jī)碳分布特征及與土壤理化性質(zhì)的關(guān)系.草業(yè)科學(xué),2011,28(8):1411-1415.

Cao L H,Liu H M,Zhao S W.Istribution of soil organic carbon and its relationship with soil physical and chemical properties on degraded alpine meadows.Pratacultural Science,2011,28(8):1411-1415.(in Chinese)

[10]李洪勛.草帶在防治坡耕地土壤侵蝕中的作用.草業(yè)科學(xué),2005,22(1):94-97.

Li H X.The role of grass strip in prevention of water and soil erosion on sloping land.Pratacultural Science,2005,22(1):94-97.(in Chinese)

[11]牟金澤,孟慶枚.陜北部分中小流域輸沙量計(jì)算.人民黃河,1983(4):35-37.

Mou J Z,Meng Q M.Part of Shaanxi small watershed sediment discharge calculation.Yellow River,1983(4):35-37.(in Chinese)

[12]鄭粉莉,劉峰,楊勤科,江忠善.土壤侵蝕預(yù)報(bào)模型研究進(jìn)展.水土保持通報(bào),2001,21(6):17-18.

Zheng F L,Liu F,Yang Q K,Jiang Z S.Review of research progress in soil erosion prediction model.Bulletin of Soil and Water Conservation,2001,21(6):17-18.(in Chinese)

[13]王文娟,張樹文,李穎,卜坤.基于GIS和USLE的三江平原土壤侵蝕定量評(píng)價(jià).干旱區(qū)資源與環(huán)境，2008,22(9):112-117.

Wang W J,Zhang S W,Li Y,Bu K.Quantitative sssessment of soil erosionin Sanjiang Plain based on GIS and USLE.Journal of Arid Land Resources and Environment,2008,22(9):112-117.(in Chinese)

[14]馬琦,王琦.幾種草被植物的水土保持效應(yīng)研究.草業(yè)科學(xué),2006,22(10):72-74.

Ma Q,Wang Q.The effect of several grassy species on soil and water conservation.Pratacultural Science,2006,22(10):72-74.(in Chinese)

[15]Wischmeier W H,Johnson C B,Cross B V.A soil erodibility nomograph farmland and construction sites.Journal of Soil and Water Conservation,1971,26(5):189-193.

[16]鮑士旦.土壤農(nóng)化分析.北京:中國農(nóng)業(yè)出版社,2000.

Bao S D.Soil Agro-chemistrical Analysis.Beijing:China Agriculture Press,2000.(in Chinese)

[17]朱曉麗,李文龍,薛中正,許靜,王浩.基于3S技術(shù)的甘南州生態(tài)健康與生態(tài)承載力耦合.草業(yè)科學(xué),2011,28(6):65-71.

Zhu X L,Li W L,Xue Z Z,Xu J,Wang H.Coupling of ecological health and ecological carrying capacity in the Gannan pastoral area based on the “3S” technology.Pratacultural Science,2011，28(6):65-71.(in Chinese)

[18]Liu B Y,Nearing M A,Shi P J,Jia Z W.Slope length effects on soil loss for steep slopes.Soil Science Society of America Journal,2000,64(5):1759-1763.

[19]Liu B Y,Nearing M A,Risse L M.Slope gradient effects on soil loss for steep slopes.Transactions of American Society of Agriculture Engineers,1994,37(6):1835-1840.

[20]吳禮福.黃土高原土壤侵蝕模型及應(yīng)用.水土保持通報(bào),1996,16(5):29-35.

Wu L F.Soil erosion model and its application in Loess Plateau.Bulletin of Soil and Water Conservation,1996,16(5):29-35.(in Chinese)

[21]劉寶元,張科利,焦菊英.土壤可蝕性及其在侵蝕預(yù)報(bào)中的應(yīng)用.自然資源學(xué)報(bào),1999,4(4):345-350.

Liu B Y,Zhang K L,Jiao J Y.Soil erodibility and its use in soil erosion prediction model.Journal of Natural Resources,1999,4(4):345-350.(in Chinese)

(責(zé)任編輯武艷培)

Variation in soil erosion based on USLE model and remote sensing technology during 2000-2010 in Gansu Province

Wang Li-na1,2, Li Wen-long1, Wang Su-fang3, Chen Di1,Xu Jing4

(1.College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China;2.Gansu Province Environmental Monitoring Center, Lanzhou 730020, China;3.Gannan Land and Resources Bureau of Gansu Province, Gannan 747000, China;4.Lanzhou University of Finance and Economics, Lanzhou 730020, China)

Abstract:Combined USLE (Universal Soil Loss Equation) with remote sensing, geographic information technologies, the distribution area of soil erosion was evaluated quantificationally in Gansu Province from 2000 to 2010. The results showed that the mainly erosion intensity was mild and moderate erosion in Gansu Province.The soil area with significant erosion was 24 872.9, 27 503.6 and 28 000.8 km2 in 2000, 2005 and 2010, respectively. The region of soil erosion mainly distributed in the southern part of Gansu Province, and the main soil erosion area distribute in forest and grassland. The largest soil erosion area distribute in forest,followed by grassland. The main landforms types of erosion are mountain land in Longnan, Gannan and loess plateau. The variation in soil erosion from 2000 to 2010 was also analyzed, it indicated that the reasons of soil erosion were not only by uneven distribution of precipitation, land-use types and landform types, but also irrational useage of ecological environment.

Key words:soil erosion; universal soil loss equation (USLE); remote sensing

Corresponding author:Li Wen-longE-mail: wllee@lzu.edu.cn

中圖分類號(hào)：S157;S127

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1001-0629(2016)2-0176-08

通信作者：李文龍(1977-)，男，甘肅蘭州人，教授，博士，主要從事草地遙感與地理信息系統(tǒng)研究。E-mail:wllee@lzu.edu.cn

基金項(xiàng)目：國家自然科學(xué)基金(31170430、41471450)；國家社會(huì)科學(xué)基金(14CJY010)第一作者：王莉娜(1982-)，女，甘肅臨洮人，碩士，主要從事草地遙感與地理信息系統(tǒng)研究。E-mail:dq954@126.com

*收稿日期：2015-04-13接受日期：2015-04-24

DOI:10.11829/j.issn.1001-0629.2015-0203

王莉娜,李文龍,王素芳,陳迪,許靜.基于遙感和USLE模型的2000-2010年甘肅省土壤侵蝕變化評(píng)價(jià).草業(yè)科學(xué),2016,33(2)：176-183.

Wang L N,Li W L,Wang S F,Chen D,Xu J.Variation in soil erosion based on USLE model and remote sensing technology during 2000-2010 in Gansu Province.Pratacultural Science，2016,33(2)：176-183.

前植物生產(chǎn)層