王衛(wèi)林 葉燎原 楊昆 王桂林 蘇飛

摘要：土壤侵蝕空間分布特征是進行土壤侵蝕防治規(guī)劃、實踐的重要基礎與依據?；诘乩硇畔⑾到y(tǒng)（GIS）和遙感技術（RS），提取了洱海流域降雨侵蝕力、土壤可蝕性、坡度坡長、地表覆蓋、水土保持措施5個主要影響水土流失因子，采用修正的通用土壤流失模型（RUSLE）估算土壤侵蝕量，并運用GIS空間統(tǒng)計分析功能，分析了土壤侵蝕在海拔、坡度與土地利用類型等方面的空間分布特征。結果表明，洱海流域年平均土壤侵蝕模數為64.68 t/（hm2·a），屬于強度侵蝕，土壤侵蝕面積1 404.279 km2，占總流域面積的61.12%，2 000～2 500 m高程帶、15°～45°坡度帶，以及未利用土地區(qū)域是研究區(qū)土壤侵蝕防治的重點地區(qū)。

關鍵詞：洱海流域；修正的通用土壤侵蝕方程；地理信息系統(tǒng)；水土流失；空間分布特征

中圖分類號：S157 文獻標識碼：A 文章編號：0439-8114（2015）13-3108-06

DOI：10.14088/j.cnki.issn0439-8114.2015.13.010

Analysis on Spatial Characteristics of Soil Erosionin Erhai River Basin

Based on GIS and RUSLE

WANG Wei-lin1a，1b，YE Liao-yuan2，YANG Kun1a，2，WANG Gui-lin1b，2，SU Fei1b，2

（1a.School of Information Science and Technology； 1b.School of Tourism and Geographic Sciences， Yunnan Normal University， Kunming 650500， China；2.GIS Technology Research Center of Resource and Environment in Western China， Ministry of Education， Kunming 650500， China）

Abstract：The spatial distribution characteristics are important basis for the planning and application of soil erosion prevention. Based on geographic information systems（GIS） and remote sensing（RS）， five primary factors affecting soil erosion in the ErHai Lake basin were extracted，i.e.，rainfall erosivity， soil erodibility， terrain， vegetation coverage and soil and water conservation measures. Meanwhile， the amount of soil erosion was estimated using the revised universal soil loss equation （RUSLE） model. Furthermore， the spatial distribution characteristics were also analyzed with the application of GIS in the view of elevation， slope and land use types. The results showed that the soil erosion modulus were 64.68t/（hm2·a），which belonged to the grade of intensity erosion， and a erodible area of 1404.279km2（61.12% of the whole area）.The results also revealed that soil erosion ，especially soil erosion with high grade， mostly distributed in the following areas， that is elevation area from 2000 to 2500 meters， slope area from 15 to 45 degree， and unused land. Those areas were the key zone for soil erosion.

Key words： Erhai river basin； revised universal soil loss equation （RUSLE）； geographic information systems （GIS）； soil erosion； spatial distribution characteristics

土壤侵蝕是世界上的主要生態(tài)破壞現象之一[1]，它是危及人類生存和發(fā)展的重要環(huán)境問題之一。

土壤侵蝕的研究始于19世紀晚期，1887年德國土壤學家Bwald Wollny進行相關研究時，只限于表面現象的觀察和定性描述[2]。從20世紀20年代開始，有不少學者開始進行相關的定量研究，但直到60年代，此類研究基本上仍是經驗所得。隨著相鄰學科的發(fā)展、測驗方法的改進以及計算機技術的興起，逐步建立了一些模擬侵蝕物理過程的數學模型。進入80年代，隨著科學技術的迅猛發(fā)展和計算機技術的日益完善，侵蝕產沙的研究迅速發(fā)展，大量研究成果問世。特別到了90年代后，RS和GIS的發(fā)展和廣泛應用為人們進一步研究提供了有力的工具，模型的研究也開始由集總式向分布式過渡，國內外學者提出了大量的土壤侵蝕模型[3]。

中國是世界上水土流失較強烈的國家之一，定期調查并公告全國土壤侵蝕及其治理狀況，是客觀反映水土保持治理成效、評估治理效益編制治理規(guī)劃的基礎[4]。國內外對于區(qū)域土壤侵蝕綜合評價和制圖方法，可以歸納為4種方法，包括常規(guī)地理制圖法、坡面土壤侵蝕模型與GIS集成法、區(qū)域土壤侵蝕模型方法和抽樣調查方法[5，6]。中國比較常用的方法是依據水利部頒布的土壤侵蝕分級分類標準進行評價[7]。

1 研究區(qū)域與數據來源

洱海流域位于瀾滄江、金沙江和元江三大水系分水嶺地帶，屬瀾滄江—湄公河水系，流域面積2 565 km2，地理坐標在東經100°05′-100°17′、北緯25°36′-25°58′之間。洱海地處云南省大理白族自治州境內，是云南省第二高原淡水湖泊，風光明媚，素有“高原明珠之稱”。洱海流域屬低緯度高原亞熱帶季風氣候，干濕分明，氣候溫和，日照充足。流域內的地帶性土壤為紅壤，隨著海拔的變化，由低到高依次為紅壤、黃紅壤、黃棕壤、暗棕壤、亞高山草甸及高山草甸土等。

本研究基礎數據主要包括遙感、專題圖件和實測數據。遙感影像數據為2012年洱海流域30分辨率的Landsat TM多光譜影像。專題圖件數據包括矢量流域邊界圖、數字高程模型（DEM）以及土壤類型數據，DEM數據來源于國家1∶5 000數字高程模型，分辨率為25，土壤矢量類型數據來自國家自然科學基金委員會“中國西部環(huán)境與生態(tài)科學數據中心”。實測數據為洱海流域6個氣象站點1997～2010年日降雨數據，來源于云南省水文水資源局和中國氣象科學數據共享服務網。

2 模型與方法

2.1 修正通用土壤侵蝕模型

Wischmeier在20世紀60年代提出的通用土壤流失方程（universalsoillossequation，USLE）是當前應用最廣泛的經驗模型[8]。修正通用土壤侵蝕方程（reviseduniversalsoillossequation，RUSLE）從技術性和確定因子的算法上對USLE進行了改進，因而應用范圍更廣，預測精度更高。此方程形式簡單，所需的參數較易獲取，且各因子具有物理意義，被廣泛應用于土壤侵蝕風險的評估分析[9，10]。在參照國內外專家土壤侵蝕研究成果的基礎上，結合洱海流域的自然環(huán)境特征，決定采用修正的通用土壤流失方程對洱海流域的土壤進行分析，其模式為：

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

式中，A為年單位面積的土壤侵蝕量（t/hm2·a）；R為降雨侵蝕力因子（MJ·mm/hm2·h·a）；K為土壤可蝕性因子（t·h/MJ·mm）；L為坡長因子（m）；S為坡度因子（%）；C為植被覆蓋因子；P為水土保持措施因子。

2.1.1 降雨侵蝕力因子（R） 降雨侵蝕力主要衡量雨滴沖擊的動能，反映降雨這一氣候因素對土壤侵蝕的潛在作用能力。降雨侵蝕力計算利用日雨量采用下列計算公式計算降雨侵蝕力[11，12]。

R=0.444 88P0.969 82

式中，R為降雨侵蝕力系數，單位為（MJ·mm）/（hm2·h·a）；P為年平均降雨量，單位為mm。

本研究利用公式（2），根據洱海流域各氣象站點1997～2010年觀測資料，逐年計算各站點的年平均降雨量，利用ArcGIS中的IDW插值得到洱海流域降雨侵蝕力表面模型（圖1）。

2.1.2 土壤可蝕性因子（K） 土壤可蝕性因子K反映土壤被降雨侵蝕力分離、沖蝕和搬運的難易程度。本研究結合土壤類型數據，在相關學者的研究基礎上[13，14]，近似確定出研究區(qū)不同類型土壤可蝕性因子K，并將該值輸入到土壤類型矢量圖中，通過矢量—柵格數據轉換，獲取洱海流域可蝕性因子K值得分布范圍（圖2）。

2.1.3 地形因子（LS） 地形地貌特征對土壤侵蝕的影響集中表現在坡長與坡度兩方面，因此，一般用坡長（L）與坡度（S）因子估算地形因素對土壤侵蝕的影響。目前，在DEM上進行坡長計算的方法按原理有非累積流量的直接計算法、累積流量的單位匯水面積計算法和水流強度指數的間接計算法[15]，本研究坡度坡長因子LS采用彭建等[12]的計算方式獲得。

本研究以洱海流域25 m分辨率的DEM數據，獲取洱海流域可蝕性因子LS得分布范圍（圖3）。

2.1.4 植被與經營管理因子（C） 植被與經營管理因子C反映的是有關植被覆蓋和管理對土壤侵蝕的綜合作用，其值介于0～1之間[16]。本研究以2012年的TM遙感影像為數據源，在ENVI4.8遙感圖像處理軟件支持下，結合實地調研，并參考GoogleEarth影像、地形地貌以及土壤類型等資料，采用監(jiān)督分類的方法，進行人工交互解譯，得到洱海流域土地利用類型數據。根據研究需要，土地利用類型主要分為林地、草地、水田、園地、濕地、居民點、水域、未利用地八大類。根據洱海流域土地利用以及實際經營情況，并參考美國農業(yè)手冊537號[17]以及相關研究[18，19]的基礎上，對各種土地利用類型賦予不同的C值（表1）。

2.1.5 水土保持措施因子（P） P因子也稱為侵蝕防治措施因子，是指特定保護措施下的土壤流失量與未實施保護措施之前相應地地順坡耕作時的土壤流失量之比值。其值介于0～1之間[20]。通常，包含于這一因子中的侵蝕控制措施有調整水體形態(tài)、斜坡坡度等，通過改變小地形、攔蓄地表徑流、增加降雨入滲等手段減少土壤侵蝕[21]。本研究區(qū)主要以林地、水田為主，根據研究區(qū)特殊的地形環(huán)境，參照美國農業(yè)部手冊537號[17]以及相關研究的基礎上[22，23]，對各種土地利用類型賦予不同的P值（表1），生成CP因子的分布圖（圖4）。

2.2 研究方法

本研究以遙感影像為信息源，地理信息系統(tǒng)為分析工具，依據遙感、降水、土壤、數字高程模型等資料，運用修正通用土壤侵蝕模型（RUSLE），生成土壤侵蝕主要因子圖層，然后利用ArcGIS軟件，將R、K、L、S、C和P各因子圖層轉化到統(tǒng)一坐標系，像元大小為30 m×30 m進行疊加相乘處理，得到流域土壤侵蝕模數，形成洱海流域土壤侵蝕強度柵格分布圖[t/（hm2·a）]，并依據水利部水土保持司的《土壤侵蝕分類分級標準（SL190-1996）》，并考慮到研究區(qū)土壤侵蝕量的實際分布特征，得到洱海流域土壤侵蝕強度分布圖（圖5）。

3 結果與分析

3.1 流域土壤侵蝕量與土壤侵蝕強度

通過計算得到洱海流域土壤侵蝕面積1 404.279 km2，占總面積的61.12%，土壤侵蝕總量為1 685.3×105 t/a，平均土壤侵蝕模數64.68 t/（hm2·a），屬于強度侵蝕。并對土壤侵蝕強度分布圖進行分類統(tǒng)計分析，得到洱海流域各土壤侵蝕強度等級平均土壤侵蝕模數與土壤侵蝕量（表2）。

從表2可以看出，盡管極劇烈侵蝕面積比例不到6%，但其土壤侵蝕量比重卻達到80%以上，其他土壤侵蝕強度等級占的比重較低，因此，極劇烈侵蝕是研究區(qū)土壤侵蝕防治的重點。

3.2 洱海流域土壤侵蝕空間特征分析

3.2.1 土壤侵蝕的海拔分布特征 從洱海流域不同海拔帶的平均土壤侵蝕模數（表3）可以看出，盡管土壤侵蝕與海拔具有較大的相關性，但土壤侵蝕強度與海拔無直接的線性關系。除<2 000 m、>3 500 m高程帶其值較低外，其余各高程帶上的土壤侵蝕模數均較大，且差異顯著，其中，2 000～2 500 m高程帶的平均土壤侵蝕模數較高，屬極強度侵蝕，其他高程帶的土壤侵蝕分屬于強度、中度侵蝕。平均土壤侵蝕模數之所以在不同高程帶呈現不同的分布特征主要緣于研究區(qū)高程的空間分布以及植被覆蓋率的差異，即海拔在2 000 m以下高程帶，大多數屬于平原地帶，坡度較低，土地類型多為田地、濕地以及少量的草地，對水土的保持能力相對較強；而2 000～2 500 m高程帶占研究區(qū)近50%的土地面積，大多山高坡陡，土地類型多為水田、未利用土地、居民點以及少量的園地和草地，人為活動強烈，破壞性相對較大，因而，對水土保持能力相對較低；而2 500～3 000 m高程帶是洱海流域林地、其他林地等高植被覆蓋地類的集中分布區(qū)，因而，對水土流失有一定的保持能力；而3 000～3 500 m高程帶主要以林地、草地為主，但大部分區(qū)域坡度較低，其平均土壤侵蝕模數相對較低；而大于3 500 m高程帶主要以低矮的草甸為主，雖然坡度較高，但常年以低溫為主，加之凍土的作用，對水土流失起到了積極的保持作用。

由洱海流域不同高程帶的土壤侵蝕量（表3）可以看出，2 000～2 500 m高程帶其土壤侵蝕量占了土地總面積的70%左右，而2 500～3 000 m高程帶的土壤侵蝕量比例也遠高于其他高程帶比例，因此，2 000 m～2 500 m高程帶是研究區(qū)土壤侵蝕治理的重點地區(qū)，2 500 m～3 000 m高程帶上的土壤侵蝕也應引起重視。

由不同高程帶上土壤侵蝕強度等級的面積比例構成（表4）可以看出，洱海流域不同海拔帶上的土壤侵蝕強度特征各異，各海拔帶上以微度、輕度侵蝕侵蝕強度等級占絕對比例（67%～97%），強度、極強度侵蝕以及劇烈侵蝕占的面積比例普遍偏低。此外，中度侵蝕在2 500～2 500 m以及2500～300 m高程帶上占一定比例，2 000～2 500 m以及2 500～ 3 000 m高程帶上極劇烈侵蝕的面積比例也較高。

3.2.2 土壤侵蝕的坡度分布特征 由洱海流域不同坡度帶的平均土壤侵蝕模數與土壤侵蝕量（表5）可以看出：①研究區(qū)平均土壤侵蝕模數與坡度呈顯著正相關。其中，0°～6°坡度帶平均土壤侵蝕模數最低，屬微度侵蝕，6°～15°、15°～25°坡度帶屬強度侵蝕，>25°坡度帶屬于極強度侵蝕；②盡管15°～45°坡度帶占全縣面積的53%，其土壤侵蝕量卻占了總侵蝕量的78%以上。因此，15°～45°坡度帶是洱海流域土壤侵蝕防治的重點區(qū)域，其中，以15°～35°坡度帶上的土壤侵蝕治理最為關鍵。

由不同坡度帶土壤侵蝕強度等級面積比例構成（表6）可以看出，不同土壤侵蝕強度等級在各坡度帶上均有分布，但其面積構成有較大的差異：①0°～6°坡度帶上，微度侵蝕占絕對比重（87%），其他侵蝕強度等級比例較低；②6°～15°、15°～25°坡度帶仍以微度、輕度侵蝕為主（72%～80%），中度侵蝕、極劇烈侵蝕占了一定比例，其他侵蝕強度等級較低；③25°～35°、35°～45°、>45°坡度帶上的侵蝕強度等級比例分布特征基本相似，都以微度、輕度、中度侵蝕為主（76%～80%），強度、極劇烈侵蝕次之（13%～15%），其他侵蝕強度等級比例較低。

3.2.3 土壤侵蝕的土地利用類型分布特征 由洱海流域不同土地利用類型的平均土壤侵蝕模數（表7）可以看出，未利用土地的平均土壤侵蝕模數是其他土地利用類型的30倍以上（0除外）。同時，盡管未利用土地不足研究區(qū)土地面積的10%，其侵蝕量卻占到全縣土壤侵蝕總量的83%左右，林地、草地、水田與園地的土壤侵蝕量比例遠小于其土地總面積比重，因此，未利用土地是研究區(qū)土壤侵蝕防治的重點土地利用類型。

由洱海流域不同土地利用類型的土壤侵蝕強度等級面積比例構成（表8）可以看出：①林地、草地及園地以微度、輕度侵蝕為主（80%～83%），中度侵蝕也占一定比例（12%～14%），以及總面積比例不超過5%～6%的強度、極強度侵蝕，極少劇烈、極劇烈侵蝕（0.1%～3%）；②水田以微度侵蝕為主（58%），輕度（16%）、中度侵蝕（10%）次之，強度、極強度侵蝕也占一定比例（12%），其他侵蝕強度等級面積比例不足2%；③未利用土地分布在各個侵蝕強度等級，但以極劇烈侵蝕為主（54%），微度侵蝕次之（28%），其他侵蝕強度等級面積比例相差不大，均在1%～5%之間。

4 結論

采用RUSLE模型估算洱海流域年平均土壤侵蝕量，結果表明，洱海流域年平均土壤侵蝕模數為64.68 t/（hm2·a），屬強度侵蝕；而土壤侵蝕的空間特征表明，2 000～2 500 m高程帶、15°～45°坡度帶，以及未利用土地區(qū)域是研究區(qū)土壤侵蝕防治的重點地區(qū)，相關部門應采取相關措施以減少土壤流失量。其中2 000～2 500 m高程帶的防治重點在于維持原生植被的自然演替、退耕還林還草，15°～45°坡度帶宜開展坡改梯、等高耕作等生物、工程措施以保持水土，未利用土地防治的基本方向則是積極引導其向有林地、草地等環(huán)境友好型土地利用/土地覆被類型的轉型[12]。

參考文獻：

[1] 肖 寒，歐陽志云，趙景柱，等.海南島生態(tài)系統(tǒng)土壤保持空間分布特征及生態(tài)經濟價值評估[J].生態(tài)學報，2000，20（4）：552-558.

[2] MEYERL D.Evaluation of the universal soil loss equation[J].Journal of Soil and Water Conservation，1984，39：99-104.

[3] 張有全，宮輝力，趙文吉.基于遙感與GIS的土壤侵蝕強度快速估測方法[J].生態(tài)環(huán)境，2007，16（1）：102-107.

[4] 程 琳，楊勤科，謝紅霞，等.基于GIS和CSLE的陜西省土壤侵蝕定量評價方法研究[J].水土保持學報，2009，23（5）：61-66.

[5] 楊勤科，李 銳，曹明明.區(qū)域土壤侵蝕定量研究的國內外進展[J]. 地球科學進展，2006，21（9）：31-38.

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

[7] SL190-1996，土壤侵蝕分類分級標準[S].

[8] WISCHMEIER W H，SMITH D D.Predicting Rainfall Erosion Losses from Cropland in East of the Rocky Mountains[M].Washington，DC：US Department of Agriculture，1965.

[9] 姚 妤，張 沛，嚴力蛟，等.基于RUSLE和景觀安全格局的土壤侵蝕風險格局研究：以甘肅省甘南藏族自治州迭部縣為例[J]. 水土保持通報，2011，31（3）：161-167.

[10] 江青龍，謝永生，張應龍，等.京津水源區(qū)小流域土壤侵蝕及其空間分異[J].水土保持通報，2011，31（1）：249—255.

[11] 楊子生.云南省金沙江流域土壤流失方程研究[J].山地學報， 2002，20（增刊）：1-9.

[12] 彭 建，李丹丹，張玉清.基于GIS和RUSLE的滇西北山區(qū)土壤侵蝕空間特征分析——以云南省麗江縣為例[J].山地學報，2007，25（5）：584-556.

[13] 楊子生.滇東北山區(qū)坡耕地土壤可蝕性因子[J].山地學報，1999，17（增刊）：10-15.

[14] 張科利，彭文英，楊紅麗.中國土壤可蝕性值及其估算[J].土壤學報，2007，44（1）：7-13.

[15] 晉 蓓，劉學軍，甄 艷，等.ArcGIS環(huán)境下DEM的坡長計算與誤差分析[J]. 地球信息科學學報，2010，12（5）：700-705.

[16] 湯麗潔，舒 暢.基于GIS和RS的巢湖流域水土流失評估[J].水土保持通報，2013，33（1）：307-308.

[17] WISCHMEIER W H，SMITH D D. Predicting rainfall erosion losses： A guide to conservation planning[A]. USDA. Agriculture Handbook No.537[C]. Washington， DC： United States Department of Agriculture，1978.

[18] 楊子生.滇東北山區(qū)坡耕地土壤侵蝕的作物經營因子[J].山地學報，1999，17（增刊）：19-21.

[19] 張旭群，陳耀強，陳浩昆，等.基于GIS和RUSLE的粵東黃岡河流域土壤侵蝕評估[J].中國水土保持，2013（2）：34—36.

[20] 柯克比M J，摩根R P C.土壤侵蝕[M].王禮先，吳 斌，洪惜英，譯.北京：水利水電出版社，1987.

[21] 秦 偉，朱清科，張 巖.基于GIS和RUSLE的黃土高原小流域土壤侵蝕評估[J].農業(yè)工程學報，2009，25（8）：157-163.

[22] 楊子生.滇東北山區(qū)坡耕地土壤侵蝕的水土保持措施因子[J].山地學報，1999，17（增刊）：22-24.

[23] 王萬忠，焦菊英.中國的土壤侵蝕因子定量評價研究[J].水土保持通報，1996，16（5）：1-20.