基于WEPP的螞蟻河流域土壤侵蝕狀況研究

周 歡，陳圣波， 程 雪

(吉林大學地球探測科學與技術學院，吉林長春 130026)

?

基于WEPP的螞蟻河流域土壤侵蝕狀況研究

周 歡，陳圣波*， 程 雪

(吉林大學地球探測科學與技術學院，吉林長春 130026)

摘要[目的]探尋準確、有效的方法研究螞蟻河流域土壤侵蝕狀況。[方法]利用WEPP模型對螞蟻河流域的侵蝕狀況進行定量研究，得到近30年來螞蟻河流域的土壤侵蝕分布狀況，并結合模型參數進行分析。[結果] 該流域土壤侵蝕較嚴重的區域位于東南部和北部，侵蝕量總體由中間向四周先逐漸降低后又升高。對所得到的侵蝕量進行換算，結果顯示，該研究區的平均侵蝕模數范圍為0～1 020 t/(km2·a)，屬于微度和輕度侵蝕。[結論]WEPP模型的應用為東北黑土區土壤侵蝕預報提供一種可行方案，也為水土保持規劃提供依據。

關鍵詞WEPP；土壤侵蝕；螞蟻河流域

土壤侵蝕可以導致土地退化、干旱、洪澇等災害發生，引起生態環境惡化。近年來由于人類活動的增加，土壤侵蝕日益加劇。因此，研究土壤侵蝕狀況，分析影響土壤侵蝕的因素，有利于了解不同地區土壤侵蝕的空間分布，為土壤侵蝕預報及治理提供依據。

針對土壤侵蝕問題，國內外有不少學者以不同的方式進行了研究。范昊明等[1]從自然地理環境和人類活動角度出發，分析影響東北黑土區土壤侵蝕的主要因素。由于該研究主要是定性研究，較不精確。湯立群等[2]按自然流域劃分子流域，將子流域概化分區，根據不同區域土壤侵蝕和產沙特點，應用水文學、泥沙運動學基本原理建立了具有物理成因的概念性產沙數學模型。該模型在黃土丘陵區具有一定精度，在其他地區的適用性有待驗證。1965年，Wischmeier等[3]提出通用土壤流失方程USLE(Universal soil loss equation)，該模型是美國研制的用于定量預報農地或草地坡面多年平均土壤流失的一個經驗性土壤侵蝕預報模型，主要研究由降雨引起的水動力土壤侵蝕[4]。由于在應用USLE時，必須對每一個因子的應用范圍有清楚的了解，因此對于不同地區各個因子存在差異性，模型的推廣受到限制。WEPP (Water erosion prediction project) 模型是用來代替USLE的新一代侵蝕預報模型，它可以模擬和預測不同時間尺度、不同土地利用類型的徑流量和土壤侵蝕量，同時還具備模擬和預測土壤水分的滲入、蒸發、農作物生長的功能[5]。由于該模型是一個過程模型，所以可外推到更廣泛的區域，甚至是沒有試驗數據的區域[6]。基于此，筆者以黑龍江省螞蟻河流域為例，應用WEPP模型模擬流域1982～2012年的總侵蝕量，得到螞蟻河流域的侵蝕狀況，分析出現該侵蝕量分布狀況的主要原因，以期為東北黑土區土壤的侵蝕預報提供一種可行方案，也為水土保持規劃提供依據。

1材料與方法

1.1研究區概況螞蟻河流域位于黑龍江省，長白山西坡，地理坐標為127°～129° E，44°～45° N，流域面積為10 184.6 km2(圖1)。流域內的螞蟻河是松花江右岸一級支流，發源于長白山張廣才嶺地區，向西流至黑龍江省尚志市尚志縣，干流長285 km。

圖1　研究區位置Fig.1　The position of study area

1.2數據收集運行WEPP模型坡面版，分別需要氣候數據、作物管理數據、土壤數據和坡面數據，而每個數據庫的生成都需要多個參數來支持。因此，需要收集氣候、作物類型、土壤類型等參數。氣候數據采用從“中國氣象科學數據共享服務網”上下載的尚志氣象站1982～2012年的每日最高溫、最低溫及降雨量數據；作物類型由“全球土地覆蓋數據GlobCover2009”影像確定；土壤參數主要通過“中國土壤數據庫”中尚志地區土壤信息得到；坡面數據通過30 m分辨率ASTER GDEM得到的螞蟻河流域的DEM經處理得以確定。

1.3WEPP模型侵蝕量模擬

1.3.1WEPP模型坡面版基本原理。WEPP模型坡面版將坡面侵蝕分為細溝侵蝕和細溝間侵蝕。細溝侵蝕是指雨滴擊濺和坡面水流對土壤進行剝蝕和搬運的過程；細溝間侵蝕是指在細溝內土壤所發生剝蝕、搬運和沉積的過程。WEPP模型用穩態泥沙連續方程來描述泥沙運動：

(1)

式中，x表示某點沿下坡方向的距離，單位為m；G表示輸沙量，單位為kg/(s·m)；Dr表示細溝侵蝕速率，單位為kg/(s·m2)；Di表示細溝間泥沙輸移到細溝的速率，單位為kg/(s·m2)。

若水流剪切力>臨界土壤剪切力，輸沙率<泥沙輸移能力時，細溝內以搬運為主：

(2)

Dc=Kr(τf-τc)

(3)

式中，Dc表示細溝水流的剝蝕能力，單位為kg/(s·m2)；Tc表示細溝間泥沙輸移能力，單位為kg/(s·m)；Kr表示細溝可蝕性參數，單位為s/m；τf為水流剪切壓力，單位為Pa；τc為臨界剪切力，單位為Pa。

若輸沙率大于泥沙輸移能力，則以沉積為主：

(4)

式中，Vf表示有效沉積速率，單位為m/s；q為單寬水流流量，單位為m2/s；β為雨滴擾動系數[7]。

1.3.2WEPP模型運用。由于WEPP模型主要適用于坡度較緩和坡長較短坡面的模擬，因此，選取平均坡度<35°、平均坡長≤30m的位置，運行WEPP模型，得到螞蟻河流域不同坡面位置的侵蝕量，其流程如圖2所示。

2結果與分析

2.1侵蝕量模擬結果運行WEPP模型，得到193個坡面單元的30a侵蝕量。在ArcGIS中，根據這些值進行插值，得到整個流域的土壤侵蝕量總體分布圖(圖3)。由圖3可知，該流域土壤侵蝕較嚴重的區域位于東南部和北部，侵蝕量總體由中間向四周先逐漸降低后又升高。依據中華人民共和國水利部1997年頒布的東北黑土區土壤侵蝕強度分級標準(表1)[8]，對所得到的侵蝕量進行換算，結果顯示，該研究區的平均侵蝕模數范圍為0～1 020 t/(km2·a)，屬于微度和輕度侵蝕。

2.2自然因素對土壤侵蝕的影響分析螞蟻河流域的土壤侵蝕狀況，在氣候方面，螞蟻河流域屬于溫帶大陸性季風氣候。氣候差異明顯，春季風大干旱，夏季溫暖濕潤，秋季氣候涼爽，冬季干燥寒冷。氣溫、降水、季風和霜凍等隨季節變化顯著，對流域內的土壤侵蝕有一定影響。

表1東北黑土區土壤侵蝕強度分級標準

Table 1Classification standard of soil erosion intensity in black soil region of Northeast China

級別Grade平均侵蝕模數Averageerrosionmodulust/(km2·a)級別Grade平均侵蝕模數Averageerrosionmodulust/(km2·a)微度Microdegree<200強度Intensity5000～8000輕度Mild200～2500極強度Extremelyintensity8000～15000中度Moderate2500～5000劇烈Severe>15000

在植被方面，該流域的植被覆蓋度高達90%以上，由于植被覆蓋度直接影響著地面徑流和土壤侵蝕模數，植被覆蓋度越高，地表徑流量越少，土壤侵蝕模數越小[9]，因此該流域

圖3　螞蟻河流域土壤侵蝕分布情況Fig.3　Soil errosion of the Ant River basin

的土壤侵蝕量較少，侵蝕程度較輕。

土壤的機械組成及理化性質都會不同程度的影響土壤侵蝕。螞蟻河流域土壤以砂粒為主，占60%～70%，上層土壤有機質含量達8.12%。但隨侵蝕年限的延長，粗粒增多，土壤結構粗粒化，腐殖質含量降低，土壤理化性質退化，抗侵蝕能力會逐漸降低。

在地形方面，侵蝕量受坡度和坡長的影響，其中受坡度的影響更大。由于坡度影響水流速度、滲透量和徑流量，即影響地表侵蝕方式、強度和過程，在一定坡度范圍內，土壤侵蝕量隨坡度的增大而增大[10]。因此，對比螞蟻河流域坡度與土壤侵蝕量可發現，侵蝕量較大的地區其對應的坡度也比較陡，說明在坡度較陡的地區水土流失比較嚴重。

2.3人為因素對土壤侵蝕的影響人為因素包括過度放牧、森林采伐等生產建設活動。螞蟻河流域主要包含黑龍江省延壽縣、尚志市局部和方正縣局部，受這些城鎮農業開墾和工業用地的影響，土壤質量不斷下降，致使該流域的水土流失趨重[11-12]。

3結論

該研究利用WEPP模型計算出螞蟻河流域緩坡度、短坡長坡面的侵蝕量，并通過ArcGIS插值得到了整個流域的侵蝕量分布狀況。從建立WEPP模型所應用的數據出發，分析不同因素對土壤侵蝕的影響，為土壤侵蝕預報提供了一種可行方案，同時也為水土流失防治提供依據。

參考文獻

[1] 范昊明，蔡國強，王紅閃.中國東北黑土區土壤侵蝕環境[J].水土保持學報，2004,18(2)：66-70.

[2] 湯立群，陳國祥，蔡名揚.黃土丘陵區小流域產沙數學模型[J].河海大學學報，1990,18(6)：10-16.

[3] WISCHMEIER W H，SMITH D D.Predicting rainfall-erosion losses from crop land east of the Rocky Mountains.Agriculture Handbook No.282[M].Washington DC：US Department of Agriculture，1965.

[4] 潘美慧，伍永秋，任斐鵬，等.基于USLE的東江流域土壤侵蝕量估算[J].自然資源學報，2010，25(12)：2154-2164.

[5] 王建勛，鄭粉莉，江忠善，等.WEPP 模型(坡面版)在黃土高原丘陵溝壑區的適用性評價：以坡度因子為例[J].泥沙研究，2008(6):52-60.

[6] 何建林.WEPP模型預測參數在紫色土區的研究[D].重慶：西南大學，2010.

[7] 項甜甜.大興安嶺東坡侵蝕-沉積過程及其指示意義研究[D].長春：吉林大學，2014.

[8] 土壤侵蝕分級分類標準：SL 190—96[S].中華人民共和國水利部，1997：9-12.

[9] 羅志軍，劉耀林.基于RS與GIS的植被覆蓋度與水土流失關系研究：以三峽庫區秭歸縣為例[J].國土資源科技管理，2008，25(3)：6-10.

[10] 查軒，黃少燕，陳世發.退化紅壤地土壤侵蝕與坡度坡向的關系：基于GIS的研究[J].自然災害學報，2010,19(2)：322-329.

[11] 崔巍，高青峰，白音包力皋，等.中小河流治理水生態功能兩級分區研究：以哈爾濱市螞蟻河流域為例[J].中國水利水電科學研究學報，2014,(4)：394-12401.

[12] 李立新，陳英智，馬景勝.螞蟻河流域水土流失與水土保持生態建設[J].水土保持研究，2004,11(1)：97-100.

Study on Soil Erosion Status of the Ant River Basin Based on WEPP model

ZHOU Huan,CHEN Sheng-bo*,CHENG Xue

(College of Geo Exploration Science and Technology,Jilin University,Changchun,Jilin 130026)

Abstract[Objective] The aim was to find an accurate and effective method for study soil erosion status in the Ant River basin.[Method] WEPP model was used for quantitative research on the erosion of the Ant River region to get soil erosion distribution results in the nearly 30 years and combined with the model parameters for analysis.[Result] The soil erosion serious area in the basin was located in the southeast and the north,and the total amount of erosion was gradually decreased and then increased from the central to the surrounding area.The obtained erosion amount was converted.The results showed that the average erosion modulus of study area was 0-1 020 t/(km2·a),which belongs to the micro degree and slight erosion.[Conclusion] The application of WEPP model provides a feasible scheme for the northeast black soil area to predict soil erosion and the basis for soil and water conservation.

Key wordsWEPP; Soil erosion; Ant River basin

作者簡介周歡(1990- )，女,滿族，遼寧本溪人，碩士研究生，研究方向：遙感地質。*通訊作者，教授，博士生導師，從事遙感地質研究。

收稿日期2016-02-23

中圖分類號S 181

文獻標識碼A

文章編號0517-6611(2016)10-074-02