基于WEPP模型的多級(jí)梯田暴雨侵蝕模擬初步研究

賈立志, 高建恩, 許秀泉, 白先發(fā), 張?jiān)?/p>

(1.中國科學(xué)院 水利部 水土保持研究所， 陜西 楊凌 712100; 2.中國科學(xué)院大學(xué)， 北京 100049; 3.西北農(nóng)林科技大學(xué) 資源與環(huán)境學(xué)院， 陜西 楊凌 712100)

基于WEPP模型的多級(jí)梯田暴雨侵蝕模擬初步研究

賈立志1,2, 高建恩1,3, 許秀泉1,2, 白先發(fā)1,2, 張?jiān)?

(1.中國科學(xué)院 水利部 水土保持研究所， 陜西 楊凌 712100; 2.中國科學(xué)院大學(xué)， 北京 100049; 3.西北農(nóng)林科技大學(xué) 資源與環(huán)境學(xué)院， 陜西 楊凌 712100)

針對(duì)2013年6—9月份安塞縣馬家溝多級(jí)梯田遭受暴雨侵蝕嚴(yán)重等問題，運(yùn)用新一代水蝕預(yù)報(bào)模型WEPP模型對(duì)馬家溝多級(jí)梯田在暴雨條件下的侵蝕情況進(jìn)行模擬，并將WEPP模型模擬的結(jié)果與實(shí)測(cè)資料進(jìn)行了對(duì)比。結(jié)果表明：WEPP 模型能夠較好地模擬黃土丘陵區(qū)多級(jí)梯田在暴雨條件下的徑流量以及土壤侵蝕量，隨著梯田級(jí)數(shù)的增加，WEPP模型對(duì)徑流量與侵蝕量的模擬精度逐漸降低。由于受該區(qū)地下水的影響，模擬徑流量總體小于采用徑流系數(shù)法算得的徑流量，由于梯田田埂的防蝕作用，模擬土壤侵蝕量總體大于實(shí)測(cè)侵蝕量。研究為黃土丘陵區(qū)梯田防蝕以及梯田修筑提供科學(xué)依據(jù)及有效方法，為黃土丘陵區(qū)多級(jí)梯田的暴雨侵蝕預(yù)測(cè)提供了合理化建議。

多級(jí)梯田； 侵蝕； 暴雨； WEPP模型； 黃土丘陵區(qū)

水土資源流失與浪費(fèi)嚴(yán)重已成為阻礙黃土高原地區(qū)經(jīng)濟(jì)發(fā)展的一個(gè)重要的因素，梯田作為黃土高原地區(qū)廣泛推廣應(yīng)用的一種水土保持措施，為該區(qū)水土保持做出了巨大的貢獻(xiàn)[1]。梯田是指在坡耕地上沿著等高線修筑的田面平整的臺(tái)階狀耕作田塊，或按照等高線修建的田面[2]。然而由于在修筑梯田時(shí)對(duì)一些自然災(zāi)害考慮不足，當(dāng)遇到暴雨等自然災(zāi)害時(shí)，下方梯田面匯集大量徑流，梯田田埂就會(huì)遭到嚴(yán)重的沖毀[3]，不僅會(huì)造成該地水土流失、養(yǎng)分流失等問題，甚至造成農(nóng)作物減產(chǎn)甚至絕產(chǎn)，因此對(duì)于梯田侵蝕的模擬與預(yù)報(bào)對(duì)黃土丘陵區(qū)的水土保持工作具有重要的現(xiàn)實(shí)意義。本文運(yùn)用新一代水蝕預(yù)報(bào)模型WEPP模型對(duì)暴雨條件下安塞縣馬家溝多級(jí)梯田斷面侵蝕情況進(jìn)行模擬，將模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，并且探求該區(qū)多級(jí)梯田在暴雨沖刷下的侵蝕規(guī)律，以及梯田田埂與坡面林草措施在防蝕方面的作用與修筑標(biāo)準(zhǔn)。為黃土丘陵區(qū)梯田防蝕以及梯田修筑提供科學(xué)依據(jù)及有效方法，并且為黃土丘陵區(qū)多級(jí)梯田侵蝕的預(yù)測(cè)提供了合理化建議。

1 數(shù)據(jù)監(jiān)測(cè)與整理

1.1 研究區(qū)概況

安塞縣馬家溝流域位于延安市西北35 km，屬于典型的黃土丘陵溝壑區(qū)，區(qū)內(nèi)地形破碎，植被稀少，地形陡峻，水土流失嚴(yán)重[4]。其中溝蝕、面蝕嚴(yán)重，土壤侵蝕模數(shù)約為14 000 t/(km2·a)，多年平均降雨約為500 mm，80%集中在6—9月，并且多以短歷時(shí)大強(qiáng)度的暴雨為主[5]。曹新莊梯田所處位置經(jīng)緯度為：36.8983905N,109.2286867E，土壤以黃綿土為主，為典型的粉壤土。土壤容重在1.08 g/cm3左右，其中粉粒含量67.06%，黏粒含量5.53%，砂粒含量26.27%，有機(jī)質(zhì)含量0.9%[6-7]。梯田斷面總共有8級(jí)，按高程從上到下依次為1—8號(hào)梯田斷面。

1.2 監(jiān)測(cè)數(shù)據(jù)與降雨資料整理

運(yùn)用徠卡全站儀與ArcGIS軟件結(jié)合對(duì)梯田侵蝕量進(jìn)行測(cè)定與計(jì)算，具體步驟如下：(1) 在5月份降雨前運(yùn)用全站儀對(duì)多級(jí)梯田地形數(shù)據(jù)進(jìn)行測(cè)量包括經(jīng)度、緯度、高程；(2) 將測(cè)量數(shù)據(jù)導(dǎo)入ArcGIS，構(gòu)建矢量點(diǎn)格式圖層并且“內(nèi)插構(gòu)建TIN”；(3) 在9月份降雨過后再運(yùn)用全站儀對(duì)多級(jí)梯田地形數(shù)據(jù)測(cè)量并按照第2步操作，生成降雨后TIN圖；(4) 將降雨前后的兩個(gè)TIN圖層做差便可得到各級(jí)梯田土壤侵蝕量以及匯水面積。

參照多級(jí)梯田實(shí)際地形，每一級(jí)梯田所產(chǎn)生的地表徑流量都是上方梯田所匯集的降雨量。每一級(jí)上方來水量都需要根據(jù)匯水面積計(jì)算。所以，由于下級(jí)梯田的匯水面積大，所以上方來水量也比上級(jí)梯田要多，對(duì)梯田造成的侵蝕也越嚴(yán)重。多級(jí)梯田的各級(jí)產(chǎn)流量可用徑流系數(shù)法根據(jù)地表徑流量公式計(jì)算得出：

(1)

式中：W——地表徑流量(m3)；F——匯水面積(m2)；R——該區(qū)6—9月降雨量(mm)；K——徑流系數(shù)。

根據(jù)安塞試驗(yàn)站所觀測(cè)統(tǒng)計(jì)的氣象資料，馬家溝流域2013年6—9月累計(jì)降雨量為1 000.6 mm。由于各級(jí)梯田的海拔高度不同，地下水出露情況也不同，受地下水與各級(jí)梯田不同的管理措施影響下的土壤含水量也不同。根據(jù)實(shí)地測(cè)量可獲得各級(jí)梯田土壤含水量、土壤容重以及土壤入滲率數(shù)據(jù)，同時(shí)參照《延安地區(qū)實(shí)用水文手冊(cè)》[8-9]，各個(gè)梯田斷面的5月份降雨前土壤含水量、土壤容重、土壤入滲率以及徑流系數(shù)選取如表1所示。

表1 各梯田斷面土壤數(shù)據(jù)以及徑流系數(shù)選取

在監(jiān)測(cè)的2013年6—9月的暴雨中，7月持續(xù)強(qiáng)降雨累計(jì)達(dá)524 m，是該區(qū)域年平均降水量的1.04倍[5]，共計(jì)有66次單場(chǎng)日降雨，日降雨量如圖1所示，平均日降雨量為15.2 mm，最大日降雨量為100 mm，最小日降雨量為0.2 mm。其中6月20日、7月4日、7月8日、8月24日1 h最大降雨量分別為13.6，24，11.4，22.4 mm。

2 WEPP模型準(zhǔn)備

2.1 WEPP模型簡(jiǎn)介

WEPP模型是美國農(nóng)業(yè)部開發(fā)的新一代水蝕預(yù)報(bào)模型，主要包括氣候、土壤、植物生長(zhǎng)、殘留物分解、冬季過程、灌溉、水文過程、地表徑流、侵蝕等9個(gè)模塊，輸入?yún)?shù)有氣候、地形、土壤和管理措施等4個(gè)數(shù)據(jù)庫參數(shù)模塊[10]。該模型克服了傳統(tǒng)的經(jīng)驗(yàn)?zāi)Ｐ蚒SLE的缺點(diǎn)，在調(diào)整參數(shù)之后對(duì)產(chǎn)沙、產(chǎn)流、土壤水分的模擬具有較高的精度[10-12]。

圖1 馬家溝流域日降雨過程

2.2 WEPP模型數(shù)據(jù)庫建立

WEPP模型坡面版是WEPP 模型中最基本、最簡(jiǎn)單的模型版本，運(yùn)行所需要建立的四個(gè)數(shù)據(jù)庫分別是氣候、坡面、土壤、作物數(shù)據(jù)庫。

根據(jù)安塞縣試驗(yàn)站統(tǒng)計(jì)的2013年6—9月的日序列的氣象數(shù)據(jù)，包括日降雨量、太陽輻射量、風(fēng)速、風(fēng)向、最高氣溫、最低氣溫、露點(diǎn)溫度等數(shù)據(jù)。按照隨機(jī)氣候生成器CLIGEN格式建立WEPP模型運(yùn)行所需要的氣象數(shù)據(jù)庫[13-14]。在坡面數(shù)據(jù)庫的建立方面，根據(jù)WEPP模型中自帶的分段功能模塊，可將該梯田分為多個(gè)斷面[15]，通過對(duì)每個(gè)斷面的模擬可以實(shí)現(xiàn)對(duì)整個(gè)梯田的模擬，由于地勢(shì)原因，上級(jí)梯田所產(chǎn)生的徑流與泥沙會(huì)通過流失進(jìn)入下一級(jí)，并且在該級(jí)梯田發(fā)生蒸發(fā)、入滲、侵蝕等一系列過程。根據(jù)實(shí)地測(cè)量可將各級(jí)梯田分為水平田面段和梯田田坎段兩部分，梯田坡度、坡長(zhǎng)等地形數(shù)據(jù)如表2所示。

表2 各級(jí)梯田斷面坡度、坡長(zhǎng)概況

土壤數(shù)據(jù)庫總共包括6個(gè)參數(shù)，其中土壤反照率、粉粒含量、初始飽和導(dǎo)水率、細(xì)溝間土壤可蝕性在輸入土壤的砂粒含量、黏粒含量、有機(jī)質(zhì)含量、巖屑含量及陽離子交換變量之后可通過WEPP模型自帶的公式計(jì)算得到。土壤臨界剪切力、細(xì)溝土壤可蝕性、有效水力傳導(dǎo)系數(shù)則是在已有的觀測(cè)資料基礎(chǔ)之上，結(jié)合前人的研究而取得[16-17]。本文選取休閑地土壤參數(shù)分別為有效水力傳導(dǎo)系數(shù)為19.7 mm/h，細(xì)溝土壤可蝕性為0.026 s/m，土壤臨界剪切力為3.51 Pa[18]。

作物數(shù)據(jù)庫是WEPP數(shù)據(jù)庫中涉及參數(shù)最多、最為復(fù)雜的一個(gè)，其包含8個(gè)子數(shù)據(jù)庫，70多個(gè)參數(shù)，其中玉米(Zeamays)、裸地、向日葵(Helianthusannuus)、谷子(Olyzasativa)等采用梁改革[19]與郝連安[20]所建立的安塞地區(qū)的作物數(shù)據(jù)庫，由于WEPP模型沒有關(guān)于土豆(Solanumtuberosum)的數(shù)據(jù)庫，因此選用苜蓿代替。各梯田斷面作物與土壤類型如表3所示。

表3 各級(jí)梯田斷面作物與土壤概況

3 結(jié)果與分析

將WEPP模型模擬的8個(gè)梯田斷面的土壤侵蝕量、徑流量與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比分析，并且運(yùn)用Nash-Sutcliffe模型有效性對(duì)模型進(jìn)行評(píng)價(jià)[21]。在8個(gè)梯田斷面的不同坡長(zhǎng)、坡度、管理措施條件下,WEPP模型對(duì)土壤侵蝕量與徑流量的模擬精度也表現(xiàn)出不同。

3.1 徑流量模擬效果分析

將WEPP模型模擬數(shù)據(jù)與徑流系數(shù)法算得徑流量進(jìn)行對(duì)比，計(jì)算得到模擬值與徑流系數(shù)法算得值之間絕對(duì)誤差與相對(duì)誤差如表4示。模擬值與徑流系數(shù)法算得值基本吻合并且R2達(dá)到0.87為顯著相關(guān)(如圖2所示)，大于一般認(rèn)為的模型可用R2=0.5，因此WEPP模型能較好地模擬黃土丘陵區(qū)多級(jí)梯田坡面在暴雨條件下的地表徑流量。

表4為地表徑流量模擬值和徑流系數(shù)法計(jì)算得值以及誤差對(duì)比?？梢钥闯?，第1級(jí)梯田徑流量模擬值與徑流系數(shù)法計(jì)算得值吻合度最高，相對(duì)誤差僅為8.55%，主要原因?yàn)榈?級(jí)梯田產(chǎn)生徑流不受其他梯田影響，并且海拔較高，受地下水影響小，因此WEPP模型在模擬時(shí)具有較高的精度。第3，4，5級(jí)梯田模擬徑流量比徑流系數(shù)法計(jì)算得值分別高27.29%，14.38%，16.56%，因?yàn)榈?，2級(jí)梯田田埂在暴雨沖刷過程中有效地?cái)r擋了徑流，并且由于匯流面積較小，徑流量沒有沖毀梯田田埂，所以大部分徑流被攔蓄在上級(jí)梯田，沒有對(duì)下級(jí)梯田造成顯著影響。從第6級(jí)梯田開始徑流系數(shù)法計(jì)算得值比模擬徑流量大，并且隨著梯田級(jí)數(shù)增加相對(duì)誤差顯著增加，因?yàn)殡S著匯流面積的增加，土壤含水量增加明顯，梯田徑流逐漸沖毀梯田田埂并且對(duì)下級(jí)梯田徑流產(chǎn)生具有貢獻(xiàn)。第6級(jí)梯田徑流系數(shù)法計(jì)算得值與模擬徑流量間相對(duì)誤差為16.30%，由于第7，8級(jí)梯田田埂被完全沖毀，所以徑流系數(shù)法計(jì)算得值顯著高于模擬值，其相對(duì)誤差分別為30.87%，32.23%，說明WEPP模型在模擬多級(jí)梯田暴雨沖刷時(shí)，隨著梯田級(jí)數(shù)增加，精度有降低趨勢(shì)。

圖2 徑流量模擬值與徑流系數(shù)法計(jì)算得值對(duì)比

通過對(duì)比可以看出，模型模擬的徑流量總體上明顯小于徑流系數(shù)法計(jì)算得值，分析主要原因是該區(qū)地下水豐富，可使土壤保持在相對(duì)濕潤狀態(tài)，使得土壤含水量較高，這樣降水不容易下滲，更容易使土壤含水量達(dá)到飽和狀態(tài)，對(duì)于地表徑流的補(bǔ)給起到了重要的作用[22]，WEPP模型在對(duì)該區(qū)地表徑流量模擬時(shí)對(duì)地下水的含量考慮不夠充分，所以對(duì)于地表徑流量的模擬值要小于徑流系數(shù)法計(jì)算得值。所以在運(yùn)用WEPP模型對(duì)該地區(qū)徑流量模擬時(shí)應(yīng)當(dāng)充分考慮地下水對(duì)于地表徑流的補(bǔ)給以及土壤含水量的影響。

表4 地表徑流量模擬值和徑流系數(shù)法計(jì)算得值對(duì)比

注：AE 為絕對(duì)誤差的絕對(duì)值，AER 為相對(duì)誤差的絕對(duì)值。

3.2 土壤侵蝕量模擬效果分析

在土壤侵蝕量方面，WEPP模型模擬數(shù)據(jù)與實(shí)測(cè)土壤侵蝕量進(jìn)行對(duì)比，結(jié)果如圖3所示。

圖3 土壤侵蝕量模擬值與實(shí)測(cè)值對(duì)比

模擬值與實(shí)測(cè)值之間的R2達(dá)到0.79為顯著相關(guān)，因此WEPP模型能較好地模擬黃土丘陵區(qū)多級(jí)梯田坡面在暴雨條件下的土壤侵蝕。

通過對(duì)比可以看出(表5)，第1，2級(jí)梯田模擬值與實(shí)測(cè)值之間的相對(duì)誤差比較大，分別為127.05%，104.72%，主要原因?yàn)椋诘?，2級(jí)梯田田埂分別為27，26 cm，梯田田埂是分布于梯田外圍的主要用于減水減沙的埂坎，其具有攔蓄徑流、防止侵蝕、保水保肥等作用[23]，當(dāng)面臨2013年6—9月暴雨時(shí)，梯田田埂能夠有效地?cái)r蓄徑流防止土壤侵蝕的發(fā)生，梯田田埂對(duì)于緩減徑流對(duì)下級(jí)梯田沖刷起到了重要的作用，并且在一定范圍內(nèi)，梯田田埂高度越高，其對(duì)泥沙攔蓄效果越好。然而在模型坡面版中沒有相對(duì)應(yīng)的攔擋措施模塊，沒有充分考慮梯田田埂的作用，導(dǎo)致模擬值遠(yuǎn)遠(yuǎn)大于實(shí)測(cè)值。

表5 土壤侵蝕量模擬值與實(shí)測(cè)值對(duì)比

注：AE 為絕對(duì)誤差的絕對(duì)值，AER 為相對(duì)誤差的絕對(duì)值

從數(shù)據(jù)分析來看，第8級(jí)梯田模擬侵蝕量與實(shí)測(cè)侵蝕量之間的相對(duì)誤差為169.01%，該級(jí)梯田田埂高度最高為28 cm(表6為梯田田埂暴雨沖蝕情況)，所以起到攔蓄徑流以及防蝕效果最好，WEPP模型的模擬效果也就最差。根據(jù)蔣定生[24]在中國科學(xué)院安塞實(shí)驗(yàn)站的觀測(cè)表明，梯田田埂高度超過25 cm時(shí)能夠很好地保持水土，在整個(gè)訊期都不發(fā)生水土流失。然而在面臨極端暴雨時(shí)，尤其是僅2013年7月12日一天降雨量達(dá)到100 mm，多級(jí)梯田特別是下級(jí)梯田土壤含水量接近飽和，降雨無法入滲，形成大量的徑流，對(duì)梯田田埂沖刷嚴(yán)重，造成梯田田埂被嚴(yán)重沖毀，無法起到蓄水保土、防止侵蝕的目的。因此對(duì)于多級(jí)梯田的下級(jí)梯田應(yīng)當(dāng)提高防蝕措施的設(shè)計(jì)標(biāo)準(zhǔn)，可將梯田田埂的設(shè)計(jì)年限由原來的20 a一遇的暴雨提高至50 a一遇的暴雨，才能切實(shí)地起到防止土壤侵蝕的目的。

在第 3，6 級(jí)梯田模擬值與實(shí)測(cè)值誤差不大，相對(duì)誤差僅為29.23%與32.87%，因?yàn)榈?，6級(jí)梯田為裸地，人為干預(yù)較少，并且梯田田埂高度分別為21 cm和22 cm，沒有能夠有效地?cái)r蓄地表徑流來防止土壤侵蝕的發(fā)生，WEPP模型在對(duì)其進(jìn)行模擬時(shí)能夠較好地反映實(shí)地情況。

表6 梯田田埂暴雨沖蝕情況

4 結(jié) 論

(1) WEPP模型能夠較好地模擬黃土丘陵區(qū)多級(jí)梯田在暴雨條件下的產(chǎn)流量與土壤侵蝕量，隨著梯田級(jí)數(shù)的增加WEPP模型的模擬精度有降低趨勢(shì)，并且對(duì)產(chǎn)流量的模擬精度要高于土壤侵蝕量的精度。

(2) WEPP模型對(duì)多級(jí)梯田暴雨沖刷的徑流量模擬效果較好，由于該區(qū)地下水豐富的原因，徑流系數(shù)法計(jì)算得值總體高于模擬值。WEPP模型對(duì)下級(jí)梯田的模擬精度低于上級(jí)梯田。

(3) 由于WEPP模型缺少攔擋措施模塊，在模擬土壤侵蝕量時(shí)，模擬值總體上要大于實(shí)測(cè)值。隨著梯田級(jí)數(shù)增加，模擬值與實(shí)測(cè)值之間的相對(duì)誤差值增加顯著。梯田田埂在應(yīng)對(duì)暴雨沖刷時(shí)對(duì)多級(jí)梯田的防蝕起到了重要的作用。

[1] 吳普特,高建恩.黃土高原水土保持新論[M].鄭州:黃河水利出版社,2006.

[2] 辛樹幟,蔣德麟.中國水土保持概論[M].北京:農(nóng)業(yè)出版社,1982.

[3] 吳發(fā)啟,張玉斌,宋娟麗,等.水平梯田環(huán)境效應(yīng)的研究現(xiàn)狀及其發(fā)展趨勢(shì)[J].水土保持學(xué)報(bào),2003,17(5):28-31.

[4] 趙誠信.安塞縣水土保持實(shí)驗(yàn)區(qū)的土壤侵蝕情況[J].泥沙研究,1986(1):81-87.

[5] 陳鳳娟,張曉萍,付艷玲,等.陜北安塞縣馬家溝流域近20a來土地利用變化分析[J].水土保持通報(bào),2011,31(1):244-248.

[6] 唐科明,汪邦穩(wěn),曹穎.黃土高原土壤分離速率的季節(jié)變化研究[J].水土保持學(xué)報(bào),2010,24(2):57-60,65.

[7] Zhang Guanghui, Tang Mingke, Zhang Xingchang. Temporal variation in soil detachment under different land uses in the Loess Plateau of China [J].Earth Surface Process and Landforms, 2009,34(9):1302-1309.

[8] 延安地區(qū)水電局.《延安地區(qū)實(shí)用水文手冊(cè)》[M].延安地區(qū)水電局,1971.

[9] 盧龍彬,付強(qiáng),黃金柏.黃土高原北部水蝕風(fēng)蝕交錯(cuò)區(qū)產(chǎn)流條件及徑流系數(shù)[J].水土保持研究,2013,20(4):17-23.

[10] 張玉斌,鄭粉莉,賈媛媛.WEPP模型概述[J].水土保持研究,2004,11(4):146-149.

[11] 梁改革,高建恩,韓浩,孟巖,幸定武,趙文君.基于作物需水與降雨徑流調(diào)控的隔坡梯田結(jié)構(gòu)優(yōu)化[J].中國水土保持科學(xué),2011,9(1):24-32.

[12] 王建勛.WEPP模型(坡面版)在黃土高原丘陵溝壑區(qū)的適用性評(píng)價(jià)[D].陜西楊凌：西北農(nóng)林科技大學(xué),2006.

[13] Flanagan D C, Meyer C R, Yu B, et al. Evaluation and enhancement of the CLIGEN weather generator[C]∥International Symposium of ASAE and Soil Erosion Research for the 21stcentury. 2001:107-110.

[14] Elliot W J, Arnold C D. Validation of the weather generator CLIGEN with precipition data from UGANDA[J]. Trans.ASAE,2001,44(1):53-58.

[15] Frankenberger D C, Flanagan J R. Water erosion prediction project (WEPP) windows interface tutorial[M]. Indiana: USDA-Agricultural Research Service & Purdue University National Soil Erosion Research Laboratory West Lafayette, 2004.

[16] 田均良,良一民,劉普靈.黃土高原丘陵溝壑區(qū)中尺度生態(tài)農(nóng)業(yè)建設(shè)探索[M].鄭州:黃河水利出版社,2003.

[17] 王建勛,鄭粉莉,江中善.WEPP模型坡面版在黃土丘陵溝壑區(qū)的適用性評(píng)價(jià)[J].水土保持通報(bào),2007,27(2):50-55.

[18] 郝連安,高建恩,張?jiān)?等.基于土壤水分動(dòng)態(tài)的梯田蘋果園水窖配置[J].中國水土保持科學(xué),2012,10(3):57-63.

[19] 梁改革.基于降雨徑流調(diào)控的坡地梯田結(jié)構(gòu)優(yōu)化初步研究[D].北京:中國科學(xué)院研究生院,2010.

[20] 郝連安.黃土丘陵溝壑區(qū)隔坡梯田結(jié)構(gòu)及集蓄措施配置優(yōu)化[D].陜西楊凌:西北農(nóng)林科技大學(xué),2012.

[21] Zhang X C. Calibration, refinement, and application of the WEPP model for simulation climatic impact on wheat production[J]. American Society of Agricultural Engineer, 2004,47(4):1075-1085.

[22] 陸垂裕,孫青言,李慧,等.基于水循環(huán)模擬的干旱半干旱地區(qū)地下水補(bǔ)給規(guī)律評(píng)價(jià)[J].水利學(xué)報(bào),2014,45(6):701-711.

[23] 周孚明.山地梯田埂坎設(shè)計(jì)探討[J].水土保持研究,2004,11(3):281-282.

[24] 蔣定生.三川河和無定河流域的基本農(nóng)田建設(shè)述評(píng)[J].水土保持學(xué)報(bào),1989,3(4):61-66.

StudyonRainstormErosionSimulationofMulti-LevelTerracesBasedontheWEPPModel

JIA Li-zhi1,2, GAO Jian-en1,3, XU Xiu-quan1,2, BAI Xian-fa1,2, ZHANG Yuan-xing3

(1.InstituteofSoilandWaterConservation,CAS&MWR,Yangling,Shaanxi712100,China; 2.UniversityofChineseAcademyofSciences,Beijing100049,China; 3.CollegeofNaturalResourcesandEnvironment,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

With respect to the prob lem of serious rainstorm erosion of multi-level terraces in Majiagou river basin from June to September in 2013, this article used the WEPP model to simulate the soil erosion of rainstorm of multi-level terraces, and then compared the result of simulation with the observed data. The results show that WEPP model can well simulate the runoff and soil loss of multi-level terraces under rainstorm in Loess Hilly and Gully Region. With the increase of terrace level the precision of WEPP model of the simulation on runoff and soil loss reduced; simulation of runoff is less than the result calculated by runoff coefficient method due to the influence of the groundwater. Simulation of soil erosion is greater than the measured erosion as a whole because of the prevention erosion of the terrace ridge. This research result provided scientific basis and effective method for the terrace anticorrosive and terraces buildup and provided scientific basis for the forecast of multi-level terrace erosion under rainstorm in the loess hilly region.

multi-level terraces; soil erosion; rainstorm; WEPP model; loess hilly region

2014-04-24

：2014-06-20

國家“十二五”科技支撐計(jì)劃項(xiàng)目課題“黃土丘陵溝壑區(qū)水土保持與高效農(nóng)業(yè)關(guān)鍵技術(shù)集成示范”(2011BAD31B05);國家自然科學(xué)基金“復(fù)雜下墊面暴雨徑流侵蝕相似性模擬實(shí)驗(yàn)研究”(41371276);中國科學(xué)院水利部水土保持研究所知識(shí)創(chuàng)新工程專項(xiàng)(A315021304)。

賈立志(1988—),男,山東濱州人,碩士研究生,主要研究方向?yàn)樗临Y源高效利用。E-mail:jializhi10000@126.com

高建恩(1962—),男,山西運(yùn)城人,博士,博士生導(dǎo)師,研究員,主要研究方向:徑流調(diào)控河流治理。E-mail:gaojianen@126.com

S157

：A

：1005-3409(2014)06-0001-05