楊曉+黎武+羊秀娟

摘 要：土壤侵蝕預(yù)報(bào)模型是由美國(guó)引進(jìn)的，經(jīng)過(guò)不斷地消化與改進(jìn)，現(xiàn)已成為我國(guó)進(jìn)行水土流失計(jì)算與預(yù)測(cè)的重要工具。該文簡(jiǎn)單介紹了土壤侵蝕預(yù)報(bào)模型的發(fā)展歷程、公式及其應(yīng)用。

關(guān)鍵詞：土壤侵蝕預(yù)報(bào)模型；土壤流失方程；研究進(jìn)展

中圖分類號(hào) S159 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731（2017）13-0072-03

Review of Research Progress in Soil Erosion Prediction Model

Yang Xiao et al.

（School of Land and Resources，China West Normal University，Nanchong 637009，China）

Abstract：The soil erosion predication model was imported from America. Though researchers continuous improve and absorb that has became an important tool in calculation and prediction of soil erosion. This paper briefly introduces the course of its development，formula and application.

Key words：Soil erosion prediction model；USLE；Research progress

1 歷史背景

RUSLE（Revised Universal Soil Loss Equation）模型為修正的通用土壤流失模型，是1997年由美國(guó)農(nóng)業(yè)部自然資源保護(hù)局（NRCS）正式?jīng)Q定實(shí)施的一個(gè)模型[1]，其前身為USLE（Universal Soil Loss Equation），即通用土壤流失方程，它的提出在美國(guó)及其他國(guó)家的土壤侵蝕方面都得到了廣泛的響應(yīng)。20世紀(jì)50年代開始，美國(guó)以在落基山脈以東地區(qū)8250個(gè)徑流試驗(yàn)小區(qū)收集的資料為依據(jù)，在Wischmeier的領(lǐng)導(dǎo)下得到通用土壤流失方程USLE，并作為定量坡面年平均土壤流失量的經(jīng)驗(yàn)性土壤侵蝕預(yù)報(bào)模型[2-3]，1956年美國(guó)農(nóng)業(yè)部頒布的農(nóng)業(yè)手冊(cè)282號(hào)標(biāo)志著土壤流失方程的正式推廣及應(yīng)用，但其在全球的通用性具有一定的局限性，并不完全適用于各國(guó)的實(shí)際情況。我國(guó)于20世紀(jì)40年代建立第一個(gè)水土保持試驗(yàn)站——天水水土保持試驗(yàn)區(qū)[4]，在大量研究和觀測(cè)的基礎(chǔ)上，出現(xiàn)了許多具有區(qū)域性的土壤侵蝕定量模型。20世紀(jì)80年代開始，引進(jìn)通用土壤流失方程與中國(guó)土壤侵蝕特點(diǎn)相結(jié)合，進(jìn)行了一系列的訂正研究[5]，開始建立具有中國(guó)特色的土壤侵蝕預(yù)報(bào)模型。隨著研究的深入，研究者們開創(chuàng)了將GIS、RS與USLE進(jìn)行結(jié)合的新思路，并得到了推廣[6-7]，GIS與USLE二者的結(jié)合不僅能夠體現(xiàn)土壤侵蝕的空間異質(zhì)性，還能夠有效地提高水土保持措施空間布置的合理性[2]，其較為全面的考慮了影響土壤侵蝕的各個(gè)因素，且具有一定的精度[8]，并將其應(yīng)用到了對(duì)地震危險(xiǎn)性與地震風(fēng)險(xiǎn)的評(píng)估[9]。1992年美國(guó)農(nóng)業(yè)部農(nóng)業(yè)研究局（USDA-ARS）提出了USLE的新一代模型RUSLE[1]，該模型較USLE具有很大的改善，其數(shù)據(jù)范圍更為廣泛，侵蝕因子的計(jì)算也更加精確，適用范圍更廣[8]，目前已發(fā)展成為流域管理規(guī)劃及對(duì)土壤侵蝕風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)的有效手段[10]，之后美國(guó)農(nóng)業(yè)部又先后頒布了RUSLE2、RUSLE3等相關(guān)模型。由于各研究區(qū)實(shí)地情況不同，直接采用美國(guó)的因子取值并不合適，故在我國(guó)各個(gè)地區(qū)具有屬于自己的RUSLE，便于方便地區(qū)水土流失的計(jì)算與預(yù)測(cè)。

2 基本公式及改進(jìn)

RUSLE模型與USLE模型的結(jié)構(gòu)基本相同，主要有降雨侵蝕力R，標(biāo)準(zhǔn)小區(qū)條件下土壤可侵蝕量K，坡長(zhǎng)與坡度因子LS，覆蓋與管理因子C，水土保持措施因子P，其基本公式為：

[A=R×K×LS×C×P] （1）

式中，A即為預(yù)測(cè)土壤侵蝕量[1]。

劉寶元通過(guò)借鑒美國(guó)土壤侵蝕預(yù)報(bào)模型建立了具有中國(guó)特色的中國(guó)土壤流失方程CSLE（Chinese Soil Loss Equation），簡(jiǎn)單實(shí)用，其公式[11，12]為：

[A=R×K×L×S×B×E×T] （2）

其中，A為多年平均土壤流失量，R為降雨侵蝕力，K為土壤可蝕性，L、S分別為坡長(zhǎng)、坡度因子，B、E、T分別為水土保持的生物措施、工程措施、耕作措施因子。

馮泮極通過(guò)計(jì)算及單位換算，得到適用于福山流域的通用土壤流失方程，其公式[13]為：

[A=0.224RKLSCP] （3）

其中，0.224為美制單位換算成kg/m2的換算系數(shù)。

孫文義運(yùn)用修正的土壤流失方程對(duì)黃土高原不同生態(tài)系統(tǒng)的土壤保持量進(jìn)行估算，其公式[14]為：

[Ac=Ap-Ar=RKLS（1-CP）] （4）

其中，Ac是單位面積上土壤保持量，Ap是單位面積上潛在的土壤侵蝕量，Ar是單位面積上實(shí)際土壤侵蝕量；因未考慮植被覆蓋和水土保持措施情況下的土壤侵蝕量，故C=1，P=1。

此外，還有金沙江流域土壤流失方程[15]，與SDR相結(jié)合的流域土壤流失脆弱性評(píng)價(jià)模型[16]，與GIS技術(shù)相結(jié)合的土壤侵蝕預(yù)報(bào)模型[17，18]等。

3 研究進(jìn)展

經(jīng)過(guò)不斷地探索、實(shí)踐、發(fā)展與完善，土壤侵蝕預(yù)報(bào)模型已經(jīng)成為計(jì)算土壤侵蝕量的重要工具，取得了顯著的成果，并在各地區(qū)出現(xiàn)了因地制宜的土壤侵蝕預(yù)報(bào)模型，對(duì)各地區(qū)水土流失治理與防治工作奠定了堅(jiān)定的基礎(chǔ)[19]。

降雨侵蝕力是降雨引起的土壤侵蝕潛在的能力，是降雨的物理性質(zhì)的函數(shù)[20]。美國(guó)學(xué)者Wischmeier認(rèn)為判斷土壤流失量最好的指標(biāo)是暴雨動(dòng)能與30min最大降雨強(qiáng)度的乘積，其表達(dá)式為[R=E·Iao][21]，但其算法數(shù)據(jù)要求較高，獲取較為困難；我國(guó)降雨侵蝕力基本沿用了其形式，但在不同的地區(qū)由于降雨動(dòng)能和時(shí)段雨強(qiáng)的不同出現(xiàn)了不同的組合形式[22]。學(xué)者們積極地改進(jìn)土壤侵蝕模型，以便本土降雨侵蝕力因子的計(jì)算。如張憲奎[23]的黑龍江省R因子的最佳指標(biāo)：[R=E60I30]；黃炎和[24]閩東南地區(qū)的[R=1120.019P1.5682i]；曲麗英[25]福建省的[R=i=112（-2.6398+0.3046Pi）]，等等。

土壤可蝕性是20世紀(jì)30年代由Middeton提出的評(píng)價(jià)土壤受侵蝕侵蝕營(yíng)力破壞性能大小的一個(gè)指標(biāo)，也是土壤對(duì)分離和搬運(yùn)的侵蝕作用的敏感性[26]，是指在雨滴擊濺和徑流沖刷等外營(yíng)力的作用下，土壤被分散與搬移的難易程度，而土壤可蝕性指標(biāo)是指試驗(yàn)中利用土壤流失量對(duì)土壤中某些性質(zhì)進(jìn)行測(cè)定來(lái)對(duì)土壤可蝕性進(jìn)行評(píng)價(jià)[27]；受空間變化、土壤性質(zhì)以及人類活動(dòng)等因素的影響，土壤可蝕性測(cè)定的方法也存在差異[26]。當(dāng)前研究在相關(guān)概念的定義、指標(biāo)的選取、分析時(shí)空分異特征等方面都具有一定的成果[28]，如吳昌廣[29]對(duì)三峽庫(kù)區(qū)采用修正后的幾何平均粒徑模型中K值計(jì)算方法；劉斌濤[30]對(duì)青藏高原地區(qū)采用模型計(jì)算和面積加權(quán)分析方法來(lái)探索K值的空間分布特征，等等。

坡長(zhǎng)坡度復(fù)合因子是地表徑流輸沙能力的量度[31]，坡長(zhǎng)因子在很大程度上決定模型計(jì)算的精度，坡長(zhǎng)指數(shù)隨著坡度的變化而不同[32]，將坡地進(jìn)行分段計(jì)算每段的LS值，進(jìn)而求出整個(gè)坡面LS值，能夠有效地反映不同坡形下土壤侵蝕的狀況[33]。對(duì)LS的計(jì)算一般都是以DEM為基礎(chǔ)，而在DEM中對(duì)LS值進(jìn)行提取會(huì)因分辨率不同導(dǎo)致獲取的LS值不固定，故選擇合適的DEM對(duì)研究區(qū)域LS值的合理獲取十分必要。例如，王邦隱[34]研究表明對(duì)陜北縣南溝流域LS值的較準(zhǔn)確提取應(yīng)采用分辨率為5m的DEM數(shù)據(jù)；吳東亮[35]以View GIS為平臺(tái)，開發(fā)模塊對(duì)LS值在研究流域坡面上的分布進(jìn)行計(jì)算，并將其結(jié)果存放在關(guān)系型數(shù)據(jù)庫(kù)；楊勤科[36]結(jié)合Van Remortel代碼開發(fā)了LS-Reg工具，方便對(duì)較大流域地區(qū)進(jìn)行LS值的自動(dòng)計(jì)算。

覆蓋與管理因子是指地表在一定的覆蓋和管理措施下土壤流失量與在同等條件下的經(jīng)過(guò)適時(shí)翻耕、連續(xù)休閑的對(duì)照地的土壤流失量之比，不同作物生長(zhǎng)期的降雨量決定著C值的大小，降低C值能使治理土壤侵蝕的投入資金最小化[37]。人工賦值法、小區(qū)實(shí)驗(yàn)法和基于植被蓋度的遙感數(shù)據(jù)定量估算法是當(dāng)前確定C值的3個(gè)主要方法，植被覆蓋度越高，C值越小[38]；不同的土地利用方式對(duì)C值的計(jì)算有不同的公式[39]，各小區(qū)標(biāo)準(zhǔn)不同，C值的計(jì)算公式也不相同[40，41]。

水土保持措施因子是用來(lái)評(píng)價(jià)最佳管理措施的一個(gè)重要因子，是指在同等條件下，實(shí)施等高耕作、梯田等水土保持措施之后的土壤流失量與在標(biāo)準(zhǔn)小區(qū)上的比值[42]。西南土石山區(qū)水平梯田在不同坡度下具有不同的P值[43]；以遼寧土地利用類型為統(tǒng)計(jì)單元，通過(guò)研究計(jì)算各水土保持因子值為水土保持提供可靠的參考[44]；劉寶元將水土保持措施細(xì)分為生物措施、工程措施、耕作措施用來(lái)規(guī)范各因子的計(jì)算[42]。

隨著研究者對(duì)土壤侵蝕預(yù)報(bào)模型研究的不斷深入，其應(yīng)用相當(dāng)廣泛，與GIS、RS等技術(shù)的結(jié)合使得土壤侵蝕預(yù)報(bào)模型成功的在農(nóng)業(yè)、礦業(yè)、林業(yè)等各方面的應(yīng)用中占據(jù)了一定的地位，為各研究區(qū)水土流失的預(yù)防與治理提供了一定的手段。在未來(lái)的發(fā)展中，土壤侵蝕預(yù)報(bào)模型應(yīng)進(jìn)一步完善，在操作計(jì)算方面更加簡(jiǎn)化，實(shí)現(xiàn)應(yīng)用的普遍性，使之成為進(jìn)行水土流失治理與防治的強(qiáng)有力的手段。

參考文獻(xiàn)

[1]陳云明，劉國(guó)彬，鄭粉莉，等.RUSLE侵蝕模型的應(yīng)用及進(jìn)展[J].水土保持研究，2004，11（4）：80-83.

[2]傅世鋒，查軒.基于GIS和USLE的東圳庫(kù)區(qū)土壤侵蝕量預(yù)測(cè)研究[J].地球信息科學(xué)，2008，10（3）：390-395.

[3]K.G.RENARD，G.R.FOSTER，G.A.WEESIES，et al. Predicting soil erosion by water ：A guide to conservation planning with the Revised Universal Soil Loss Equation[J].Handbook Department of Agriculture Washington DC，1993，（703）：

[4]許國(guó)華.羅德民博士與中國(guó)的水土保持事業(yè)[J].中國(guó)水土保持，1984（1）：41-44.

[5]謝云，林燕，張巖. 通用土壤流失方程的發(fā)展與應(yīng)用[J].地理科學(xué)進(jìn)展，2003，22（3）：179-187.

[6]卜兆宏，唐萬(wàn)龍，楊林章，等.水土流失定量遙感方法新進(jìn)展及其在太湖流域的應(yīng)用[J].土壤學(xué)報(bào)，2003，40（1）：1-9.

[7]E.Z.NYAKATAWA，K.C.REDDY，J.L.LEMUNYON. Predicting soil erosion in conservation tillage cotton production systems using the revised universal soil loss wquation（RUSLE）[J].Soil Tillage Research，2001，（57）：213-224.

[8] 李宏偉，鄭鈞瀠，彭慶衛(wèi)，等.國(guó)外土壤侵蝕預(yù)報(bào)模型研究進(jìn)展[J].中國(guó)人口·資源與環(huán)境，2016，26（S1）：183-185.

[9]奈科拉索娃，柯索波科夫，陶夏新，等.借助USLE識(shí)別中國(guó)中部地震活動(dòng)性空間模式[J].世界地震工程，2016，32（2）：170-175.

[10]周璟，張旭東，何丹，等.基于GIS與RUSLE的武陵山區(qū)小流域土壤侵蝕評(píng)價(jià)研究[J].長(zhǎng)江流域資源與環(huán)境，2011，20（4）：468-474.

[11]蔡強(qiáng)國(guó)，劉紀(jì)根.關(guān)于我國(guó)土壤侵蝕模型研究進(jìn)展[J].地理科學(xué)進(jìn)展，2003，22（3）：142-150.

[12]鄭粉莉，劉峰，楊勤科，等.土壤侵蝕預(yù)報(bào)模型研究進(jìn)展[J].水土保持通報(bào)，2001，21（6）：16-18+32.

[13]馮泮極，遲明春，曲延恩.通用土壤流失方程在福山流域的應(yīng)用[J].山東水利科技，1997（4）：14-16.

[14]孫文義，邵全琴，劉紀(jì)遠(yuǎn).黃土高原不同生態(tài)系統(tǒng)水土保持服務(wù)功能評(píng)價(jià)[J].自然資源學(xué)報(bào)，2014，29（3）：365-376.

[15]楊子生.云南省金沙江流域土壤流失方程研究[J].山地學(xué)報(bào)，2002，20（S1）：1-9.

[16]陳煉鋼，錢新，施勇，等.基于RUSLE和SDR的香溪河流域土壤流失脆弱區(qū)識(shí)別[J].中國(guó)科技論文，2012，7（5）：395-402.

[17]范麗麗，沈珍瑤，劉瑞民.基于GIS的大寧河流域土壤侵蝕評(píng)價(jià)及其空間特征研究[J].北京師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2007，43（5）：563-566.

[18]龍?zhí)煊?，喬敦，安?qiáng)，等.基于GIS和RULSE的三峽庫(kù)區(qū)土壤侵蝕量估算分析[J].灌溉排水學(xué)報(bào)，2012，31（2）：33-37.

[19]張艷靈，張紅.通用土壤流失方程研究進(jìn)展[J].山西水土保持科技，2013（2）：12-15.

[20]王萬(wàn)忠.黃土地區(qū)降雨侵蝕力R指標(biāo)的研究[J].中國(guó)水土保持，1987（12）：36-40+67.

[21]王萬(wàn)忠.黃土地區(qū)降雨特性與土壤流失關(guān)系的研究Ⅱ——降雨侵蝕力指標(biāo)R值的探討[J].水土保持通報(bào)，1983，（5）：62-64，26.

[22]張文源. 貴州喀斯特黃壤坡面降雨侵蝕特征研究[D].北京：北京林業(yè)大學(xué)，2015.

[23]張憲奎，盧秀琴，詹敏，等.土壤流失預(yù)報(bào)方程中R指標(biāo)的研究[J].水土保持科技情報(bào)，1991（4）：49-48.

[24]黃炎和，盧程隆，鄭添發(fā)，等.閩東南降雨侵蝕力指標(biāo)R值的研究[J].水土保持學(xué)報(bào)，1992 6（4）：1-5.

[25]曲麗英，劉廉海. 廈門市降雨侵蝕力指標(biāo)R值研究[J].中國(guó)水利，2005（16）：60-61.

[26]宋陽(yáng)，劉連友，嚴(yán)平，等.土壤可蝕性研究述評(píng)[J].干旱區(qū)地理，2006 29（1）：124-131.

[27]王彬，鄭粉莉，M R M M J.水蝕過(guò)程的土壤可蝕性研究述評(píng)[J].水土保持研究，2013，20（1）：277-286.

[28]魏慧，趙文武，王晶. 土壤可蝕性研究述評(píng)[J].應(yīng)用生態(tài)學(xué)報(bào)，2017（8）：1-13.

[29]吳昌廣，曾毅，周志翔，等.三峽庫(kù)區(qū)土壤可蝕性K值研究[J].中國(guó)水土保持科學(xué)，2010，8（3）：8-12.

[30]劉斌濤，陶和平，史展，等.青藏高原土壤可蝕性K值的空間分布特征[J].水土保持通報(bào)，2014，34（4）：11-16.

[31]LAND.MOORE，JOHNP.WILSON，陳奇伯. 坡長(zhǎng)、坡度因素的簡(jiǎn)單計(jì)算方法（Ⅰ）[J].水土保持科技情報(bào)，1995，（2）：30-33.

[32]秦偉，朱清科，張巖. 通用土壤流失方程中的坡長(zhǎng)因子研究進(jìn)展[J].中國(guó)水土保持科學(xué)，2010，8（2）：117-124.

[33]張照錄，崔繼紅. 通用土壤流失方程最新研究改進(jìn)分析[J].地球信息科學(xué)，2004，6（4）：51-55.

[34]汪邦穩(wěn)，楊勤科，劉志紅，等.基于DEM和GIS的修正通用土壤流失方程地形因子值的提取[J].中國(guó)水土保持科學(xué)，2007，5（2）：18-23.

[35]吳東亮，劉鵬舉，唐小明，等.基于GIS的柵格化坡面徑流路徑模擬與LS值計(jì)算[J].北京林業(yè)大學(xué)學(xué)報(bào)，2001，23（5）：10-14.

[36]楊勤科，郭偉玲，張宏鳴，等.基于DEM的流域坡度坡長(zhǎng)因子計(jì)算方法研究初報(bào)[J].水土保持通報(bào)，2010，30（2）：203-206+211.

[37]馮強(qiáng)，趙文武.USLE/RUSLE中植被覆蓋與管理因子研究進(jìn)展[J].生態(tài)學(xué)報(bào)，2014，34（16）：4461-4472.

[38]鄂立思.烏裕爾河流域土壤侵蝕時(shí)空分布特征及影響因素研究[D].哈爾濱：哈爾濱師范大學(xué)，2015.

[39]許晨純.小流域尺度土壤侵蝕評(píng)價(jià)指數(shù)模型的修正及應(yīng)用[D].南京：南京信息工程大學(xué)，2016.

[40]江忠善，王志強(qiáng)，劉志.黃土丘陵區(qū)小流域土壤侵蝕空間變化定量研究[J].土壤侵蝕與水土保持學(xué)報(bào)，1996，2（1）：1-9.

[41]劉秉正，劉世海，鄭隨定.作物植被的保土作用及作用系數(shù)[J].水土保持研究，1999，6（2）：33-37，114.

[42]范建榮，王念忠，陳光，等.東北地區(qū)水土保持措施因子研究[J].中國(guó)水土保持科學(xué)，2011，9（3）：75-78，92.

[43]劉斌濤，宋春風(fēng)，史展，等.西南土石山區(qū)水平梯田的水土保持措施因子[J].中國(guó)水土保持，2015（4）：36-39.

[44]劉振宇.基于水土保持措施因子及其分布特征探析[J].水資源開發(fā)與管理，2017（2）：37-40，44. （責(zé)編：張宏民）