張艷靈 張 紅

(山西大學(xué)環(huán)境與資源學(xué)院)

1 引言

土壤侵蝕伴隨著人類社會的發(fā)展，產(chǎn)生的環(huán)境問題越來越嚴(yán)重，對全球經(jīng)濟(jì)社會的可持續(xù)發(fā)展產(chǎn)生了很大的負(fù)面影響。全球土壤侵蝕面積約1.643 ×107km2，占地球陸地總面積的10.95%。我國是世界上遭受土壤侵蝕頗為嚴(yán)重的國家之一，據(jù)估算，每年因土壤流失而造成的間接經(jīng)濟(jì)損失在100 億元以上［1］。土壤侵蝕引發(fā)土地退化，使土地的生產(chǎn)率降低，進(jìn)而影響到產(chǎn)量和品質(zhì)。由于土壤流失，會將污染物質(zhì)帶入地表徑流中，直接或間接地對異地(侵蝕區(qū)的相鄰地區(qū)，包括位于侵蝕流域的下游地區(qū)、湖泊和近海地區(qū))的環(huán)境及經(jīng)濟(jì)的可持續(xù)發(fā)展帶來威脅，直接影響到人類的生存發(fā)展條件［2］。因此，應(yīng)用科學(xué)合理的土壤侵蝕估算模型，對區(qū)域土壤流失量進(jìn)行估算具有重要的現(xiàn)實(shí)意義。

美國學(xué)者在20 世紀(jì)70年代，提出了土壤流失通用方程USLE(Universal Soil Loss Equation)，之后又出現(xiàn)了歐洲學(xué)者的EUROSEM 和LISEM 模型，澳大利亞學(xué)者的GUEST 模型等。這幾種土壤流失量估算理論模型，幾乎都包含了與土壤流失相關(guān)的因子。歐洲學(xué)者建立的LISEM 模型，首次將土壤侵蝕模型與GIS 技術(shù)有效結(jié)合［3］，但USLE 仍是目前世界范圍內(nèi)應(yīng)用于土壤流失預(yù)報最為普遍的通用方程。我國引入土壤流失通用方程USLE 后，對其各因子進(jìn)行了深入探索，取得了顯著成效。

2 USLE 土壤流失方程的發(fā)展

通用土壤流失方程USLE 是美國學(xué)者研究建立的一種經(jīng)驗(yàn)性模型，立足研究由降雨引起的水動力土壤侵蝕，在美國最初是專門為農(nóng)田土壤侵蝕預(yù)報而創(chuàng)立的，后來逐漸延伸到其他類型土地的土壤侵蝕預(yù)報［4］。通用土壤流失方程USLE，其數(shù)學(xué)表達(dá)是影響土壤流失的6 個變量因子相乘的方程形式:

這個公式的形成，經(jīng)歷了幾十年的探索與研究。早在1917年，美國密蘇里大學(xué)M.F.Miller 就創(chuàng)立了“侵蝕小區(qū)”，為后來的土壤侵蝕研究打下了堅(jiān)實(shí)的基礎(chǔ)。1936年，H.L.Cook 通過對大量的徑流小區(qū)試驗(yàn)數(shù)據(jù)及觀測資料進(jìn)行分析，得出制約土壤侵蝕的3 個因子，即土壤可蝕性因子、降雨侵蝕力因子和植被覆蓋度因子［5］。水土保持學(xué)者辛格(zingg)，在1940年提出了能夠全面反映土壤流失因素的一般公式，根據(jù)已有的實(shí)驗(yàn)數(shù)據(jù)及觀測資料，建立了單位面積土壤流失量與坡度和坡長的指數(shù)關(guān)系:

式中，A 是單位面積土地上的平均土壤流失量;C 是一個回歸常數(shù);S 是坡度;L 是坡長［6］。

從(2)式中可以看出，辛格當(dāng)時考慮到影響土壤流失量的只是地面坡度與坡長兩個因素。

水土保持學(xué)者史密斯(Smith)于1941 -1946年發(fā)展了辛格的公式:

式中，C 是作物因子，它與(2)式中C 的含義不同;P 是水土保持措施因子。史密斯在方程中增加了植物保護(hù)土壤的作用因子C 和水土保持措施因子P。

布勞因在1947年仔細(xì)研究了不同土壤類型、輪作和坡長情況，進(jìn)一步充實(shí)了土壤侵蝕預(yù)報方程。后來，美國地域性的水土保持協(xié)會在方程中又增補(bǔ)了降雨因素，使方程基本完善起來。

1971年，Wischmeier 在對美國8 000 多個試驗(yàn)小區(qū)資料進(jìn)行研究后，改良以前的方程式，提出了現(xiàn)在流行的土壤流失通用方程式:

式中，A 為年平均土壤流失量，t/(km2·a);R為降雨徑流侵蝕因子;K 為土壤可蝕性因子;LS 為地形地貌因子(坡長、坡度因子);C 為耕作管理(植被覆蓋)因子;P 為水土保持因子。

該方程式的特點(diǎn)是考慮到的侵蝕因子更為全面，把影響土壤流失量的各因子都定量化，然后將各因子相乘，就可以估算出該地區(qū)的土壤侵蝕量，以制定符合該地條件的農(nóng)田管理或水土保持措施。

3 通用土壤流失方程的研究及應(yīng)用

3.1 降雨徑流侵蝕因子

自Wischmeier 提出以次降雨總動能E 與30 min 最大雨強(qiáng)I30之積作為降雨侵蝕力的指標(biāo)，并應(yīng)用于USLE 方程，從此，降雨侵蝕力在世界范圍內(nèi)受到學(xué)者們的熱忱關(guān)注與研究。由于每個地域的降水過程及強(qiáng)度資料獲得并不容易，因此國內(nèi)外學(xué)者研究提出了不同地區(qū)降雨侵蝕力的另外表示方法［7］。Richardson 建立的冪函數(shù)結(jié)構(gòu)日降雨量估算公式，是日降雨量估算方法中的經(jīng)典，但模型中的參數(shù)具有明顯的區(qū)域差異和季節(jié)變化，不便直接應(yīng)用。伍育鵬等［8］選擇分布在我國不同氣候帶的10 個氣象站降雨數(shù)據(jù)，對比國內(nèi)外幾種年降雨侵蝕力計(jì)算方法，得出分區(qū)域進(jìn)行擬合使用年降雨量的指數(shù)函數(shù)，用于估算我國多年平均年降雨侵蝕力，這種方法不僅便捷，而且較精確。根據(jù)降雨的時間尺度，分年、月及日降雨資料的估算方法，其中以日降雨量方法計(jì)算精度最高。謝云等［9］借助我國8 個氣象站的次降雨和日降雨資料，通過研究次、日降雨關(guān)系，簡化了EI30指標(biāo)用于估算降雨侵蝕力的方法，而且便于更廣泛的使用。但是，日最大雨強(qiáng)不易獲得，仍然限制著在區(qū)域土壤侵蝕研究中的應(yīng)用。章文波等［10］利用我國71 個代表性氣象站的資料，創(chuàng)建并利用日雨量估算降雨侵蝕力的模型，并提出了模型參數(shù)的估計(jì)方法。結(jié)果表明，此方法能夠用來估算多年平均降雨侵蝕力及其季節(jié)分布，并且估算的降雨侵蝕力精度與所在區(qū)域的降雨特征有關(guān)。

3.2 土壤可蝕性因子

土壤可蝕性因子是研究土壤侵蝕的基礎(chǔ)，是土壤侵蝕的一個非常重要的指標(biāo)，它能反映出土壤本身抵抗侵蝕的特性。國外學(xué)者對土壤可蝕性K 值進(jìn)行了很多的研究，Wischmeier 等的諾模公式以及Williams 等發(fā)展的EPIC 公式是應(yīng)用最多的經(jīng)驗(yàn)公式。我國從20 世紀(jì)50年代開始進(jìn)入土壤可蝕性研究領(lǐng)域，80年代后國外較為成熟的經(jīng)驗(yàn)公式才逐漸被引用來研究土壤可蝕性K 值，其中對這兩種經(jīng)驗(yàn)公式進(jìn)行對比研究的有鄭海金等［11］，他們創(chuàng)建了15 個標(biāo)準(zhǔn)徑流小區(qū)并經(jīng)過連續(xù)5年的定點(diǎn)監(jiān)測，運(yùn)用這兩種公式計(jì)算并分析紅壤土壤的可蝕性K 值，并以此來對比諾模公式與EPIC 公式在當(dāng)?shù)氐倪m宜性。國內(nèi)學(xué)者也有直接利用經(jīng)驗(yàn)公式來研究土壤侵蝕的，張金池等［12］通過結(jié)合傳統(tǒng)統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)，利用EPIC 公式的土壤可蝕性K 值算法，研究了鄧下小流域土壤可蝕性K 值空間分布狀況及不同植被覆蓋對其的影響。張科利等［12，13］對黃土高原地區(qū)土壤可蝕性強(qiáng)弱和其變化規(guī)律進(jìn)行了研究分析，通過回顧已經(jīng)搞過的土壤侵蝕研究成果，提出了我國土壤可蝕性標(biāo)準(zhǔn)小區(qū)定義。張科利等［14］通過野外觀測，計(jì)算出一組土壤可蝕性實(shí)測值，并通過此值分析探討了我國土壤可蝕性的估算問題，結(jié)果表明，國外現(xiàn)有的可蝕性估算模型中估算的K 值大于我國土壤可蝕性實(shí)測值。史學(xué)正等［15］用人工模擬降雨來分析我國亞熱帶土壤的可蝕性，還用諾漠公式估算土壤可蝕性K 值，由于與人工模擬降雨測得的K 值相差大，因此，他認(rèn)為該地大多數(shù)土壤不能直接用諾漠方法來估算K 值。

3.3 地形地貌因子

地形因子分為坡度因子(S)和坡長因子(L)，一般把它們放在一起作為一個因子考慮，國內(nèi)外的專家們對此進(jìn)行了大量的研究。坡度因子對侵蝕作用有著很大的影響。李鵬等［16］運(yùn)用動力原理開展了侵蝕坡度是否存在臨界值的研究，在室內(nèi)模擬實(shí)驗(yàn)基礎(chǔ)上研究了黃土高原地區(qū)坡面水動力與侵蝕產(chǎn)沙的關(guān)系，對揭示坡面侵蝕動力過程也有很大的參考價值。鄭良勇等［17］通過室內(nèi)模擬實(shí)驗(yàn)對黃土高原陡坡土壤侵蝕特性進(jìn)行了研究，結(jié)果顯示出陡坡徑流平均流速隨著徑流量和坡度的增大而呈波動增加趨勢。江忠善等［18］以我國坡面水蝕預(yù)報模型研究成果為基礎(chǔ)，加之考慮坡面侵蝕特征，提出了我國坡面水蝕預(yù)報模型的基本形式。坡長因子也對侵蝕作用有著很大的影響，我國學(xué)者們對坡長L 因子的研究多為其值的計(jì)算。zingg［19］發(fā)現(xiàn)土壤流失量與坡長L 的1.6 次方極為相關(guān)。楊子生［20］通過在滇東北山區(qū)設(shè)置試驗(yàn)小區(qū)并對其坡長進(jìn)行了連續(xù)3年的觀測，得到了相應(yīng)的坡長值。坡長因子與坡度因子的測算，除了應(yīng)用實(shí)測法外，還可依據(jù)相應(yīng)的坡長因子與坡度因子相關(guān)公式，然后使用AcrGis 中的柵格工具自動生成坡度、坡長因子值。

3.4 耕作管理(植被覆蓋)因子

早在1882年，Wollny 就通過小區(qū)實(shí)驗(yàn)來分析植被覆蓋對土壤侵蝕的影響。自20 世紀(jì)30年代開展土壤流失預(yù)報研究以來，關(guān)于水土保持措施與土壤侵蝕關(guān)系的研究一直很受重視［21］。蔡崇法等［22］根據(jù)實(shí)地調(diào)查的第一手資料，建立了植被管理因子與該指數(shù)的相關(guān)經(jīng)驗(yàn)方程，以此確定了定量計(jì)算USLE 方程中C 因子指標(biāo)的方法。楊學(xué)明等［23］利用土壤流失公式，模擬計(jì)算出了吉林省兩市黑土區(qū)坡耕地在不同土地管理方式下的土壤流失強(qiáng)度，據(jù)此得出保護(hù)性耕作能夠有效地減小吉林黑土的流失強(qiáng)度。宋現(xiàn)鋒等［24］采用多時相Landsat TM/ETM 遙感圖像和天氣模擬相結(jié)合的方法，來估算植物管理因子值。馬超飛等［25］對影像數(shù)據(jù)中的植被指數(shù)進(jìn)行歸一化處理，得到了植被管理因子，并能表現(xiàn)出植被的分布。從我國的研究中不難發(fā)現(xiàn)，對于植被管理因子的研究方法不盡一致，而且不同研究得到的植被管理因子值也有一定的差異，也不存在可比性。所以，現(xiàn)有研究得出的植被管理因子值，也只能在相應(yīng)的研究區(qū)域使用，不能推廣到更大的范圍。

3.5 水土保持因子

水土保持因子表征了水土保持措施對土壤流失的作用。國內(nèi)外學(xué)者在對水土保持因子的長期研究中，取得了許多成果。1941年，Smith D.D 首次將水土保持措施因子引用到土壤流失估算方程中對土壤侵 蝕 進(jìn) 行 估 算［26］。1965年，由 Wischmeier 與Smith［27］創(chuàng)建的土壤流失方程USLE 中，也把水土保持措施放到對土壤侵蝕影響的因子中。水土保持措施因子P 值，通常要采用實(shí)地測定法獲取，因子取值為0 -1 之間。國內(nèi)學(xué)者也對水土保持因子進(jìn)行了大量研究。水建國等［28］在研究水土保持因子時，限制農(nóng)作物的種植標(biāo)準(zhǔn)，結(jié)果得出不同土地利用方式的年水土流失情況不同，作物的覆蓋度對泥沙流失的影響有差別。符素華等［29］利用北京密云水庫20 個坡面徑流試驗(yàn)小區(qū)的實(shí)測數(shù)據(jù)，分析得出流域內(nèi)不同水土保持措施下的水土保持效益值，為水土流失治理提供了科學(xué)依據(jù)。范建榮等［30］利用不同流域總計(jì)343 次產(chǎn)流產(chǎn)沙的實(shí)測數(shù)據(jù)，計(jì)算出了7種主要水土保持措施因子值，也能為水土流失治理提供良好的科學(xué)依據(jù)。焦菊英等［31］根據(jù)黃土高原幾個地區(qū)的草地徑流小區(qū)數(shù)據(jù)，研究了降雨指標(biāo)與坡耕地侵蝕性標(biāo)準(zhǔn)之間的關(guān)系，得出結(jié)論為當(dāng)降雨指標(biāo)小于坡耕地侵蝕性標(biāo)準(zhǔn)時，草地并不會發(fā)揮其水土保持作用，反之才會產(chǎn)生水土保持的效益。

4 結(jié)語

在幾十年的時間里，通用土壤流失方程經(jīng)過眾多學(xué)者與專家對其不斷探索、實(shí)踐、發(fā)展和完善，現(xiàn)已成為計(jì)算土壤侵蝕量的基本工具，越來越受到人們的關(guān)注。通過以上的介紹，我們可以看到國外的土壤侵蝕研究已經(jīng)取得了眾多的成果，并且在土壤侵蝕模型方面占主導(dǎo)地位。不難發(fā)現(xiàn)，對于土壤侵蝕的計(jì)算是要有大量的試驗(yàn)和監(jiān)測數(shù)據(jù)做基礎(chǔ)的，并且要依據(jù)科學(xué)建模，然后因地制宜地結(jié)合當(dāng)?shù)赝寥狼治g特征來估算土壤侵蝕量［32］。我國許多學(xué)者對土壤侵蝕的研究也在不斷發(fā)展和深化，同時積累了大量的試驗(yàn)和監(jiān)測數(shù)據(jù)。在我國目前尚未有自己測算土壤侵蝕的工具及計(jì)算公式情況下，土壤通用流失方程(USLE)不失為一種好的估算公式與方法。通過USLE 公式的運(yùn)用，我們可以估算某一地區(qū)的土壤流失量，了解其水土流失狀況，從而為土地管理和防治水土流失提供一定的科學(xué)依據(jù)［33］。土壤通用流失方程經(jīng)過幾十年的發(fā)展，雖然得到了長足發(fā)展，但它仍然是一種經(jīng)驗(yàn)性的估算方法，在一些方面還是具有一定的缺陷與不足。例如，土壤流失量的估算要依據(jù)大量的統(tǒng)計(jì)數(shù)據(jù)，在數(shù)據(jù)收集的過程中難免會遇到難以收集的情況，這樣就給方程應(yīng)用造成了困難。再者，方程本身是描述定量因子的，不會顯示出數(shù)據(jù)的變化。實(shí)踐證明，USLE 不太適用于壟作、等高耕作以及那些使泥沙就地沉積的帶狀耕作措施等［34］。

［1］鄭粉莉，王占禮，楊勤科.我國土壤侵蝕科學(xué)研究回顧和展望［J］.自然雜志，2008，30(1):12-16.

［2］李占斌，朱冰冰，李 鵬.土壤侵蝕與水土保持研究進(jìn)展［J］.土壤學(xué)報，2008，45(5):802-809.

［3］丁 華，張 勇，夏 進(jìn)，等.大連市水土流失控制模擬研究［J］.環(huán)境科學(xué)研究，2007，20(3):168-172.

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

［5］汪東川，盧玉東.國外土壤侵蝕模型發(fā)展概述［J］.中國水土保持科學(xué)，2004，2(2):35-40.

［6］謝 云，林 燕，張 巖.通用土壤流失方程的發(fā)展與應(yīng)用［J］.地理科學(xué)進(jìn)展，2003，22(3):279-287.

［7］章文波，付金生.不同類型雨量資料估算降雨侵蝕力［J］.資源科學(xué)，2003，25(1):35-41.

［8］伍育鵬，謝 云，章文波.國內(nèi)外降雨侵蝕力簡易計(jì)算方法的比較［J］.水土保持學(xué)報，2001，15(3):31-34.

［9］謝 云，章文波，劉寶元.用日雨量和雨強(qiáng)計(jì)算降雨侵蝕力［J］.水土保持通報，2001，21(6):53-56.

［10］章文波，謝 云，劉寶元.利用日雨量計(jì)算降雨侵蝕力的方法研究［J］.地理科學(xué)，2002，22(6):705-711.

［11］鄭海金，楊 潔，喻榮崗，等.紅壤坡地土壤可蝕性K值研究［J］.土壤通報，2010，41(2):425-428.

［12］張金池，李海東，林 杰，等.基于小流域尺度的土壤可蝕性K 值空間變異［J］.生態(tài)學(xué)報，2008，28(5):2 199-2 206.

［13］張科利，蔡永明，劉寶元，等.黃土高原地區(qū)土壤可蝕性及其應(yīng)用研究［J］.生態(tài)學(xué)報，2001，21(10):1 687-1 695.

［14］張科利，彭文英，楊紅麗.中國土壤可蝕性值及其估算［J］.土壤學(xué)報，2007，44(1):7-13.

［15］史學(xué)正，于東升，呂喜璽.用人工模擬降雨儀研究我國亞熱帶土壤的可蝕性［J］.水土保持學(xué)報，1995，9(3):38-42.

［16］李 鵬，李占斌，鄭良勇.黃土陡坡土壤侵蝕臨界動力機(jī)制試驗(yàn)研究［J］.泥沙研究，2008(1):17-20.

［17］鄭良勇，李占斌，李 鵬.黃土高原陡坡土壤侵蝕特性試驗(yàn)研究［J］.水土保持研究，2003，10(2):47-49.

［18］江忠善，鄭粉莉.坡面水蝕預(yù)報模型研究［J］.水土保持學(xué)報，2004，18(1):66-69.

［19］Zingg A W.Degree and length of land slope as it affeets soil loss in runoff［J］.Agricultural Engineering，1940，21(2):59-64.

［20］楊子生.滇東北山區(qū)坡耕地土壤侵蝕的地形因子［J］.山地學(xué)報，1999，17(增刊):16-18.

［21］范建榮，王念忠，陳 光，等.東北地區(qū)水土保持措施因子研究［J］.中國水土保持科學(xué)，2011，9(3):75-78.

［22］蔡崇法，丁樹文，史志華，等.應(yīng)用USLE 模型與地理信息系統(tǒng)IDRISI 預(yù)測小流域土壤侵蝕量的研究［J］.水土保持學(xué)報，2000，14(2):19-24.

［23］楊學(xué)明，張曉平，方華軍.不同管理方式下吉林省農(nóng)田黑土流失量［J］.土壤通報，2003，34(5):389-393.

［24］宋現(xiàn)鋒，段 崢，牛海山，等.土壤侵蝕模型中植被管理因子的遙感估算［J］.北京林業(yè)大學(xué)學(xué)報，2009，31(3):58-63.

［25］馬超飛，馬建文，布和敖斯?fàn)?USLE 模型中植被覆蓋因子的遙感數(shù)據(jù)定量估算［J］.水土保持通報，2001，21(4):6-9.

［26］Smith D.D.Interpretation of Soil Conservation Data for Field Use［M］.Agricultural Engineering，1941:173-175.

［27］Wischmeier W H，Smith D.D.Predicting Rainfallerosion Losses from Cropland East of the Rocky Mountains［M］.Washinhton:USDA Agricultural Handbook，1965:282.

［28］水建國，孔繁根，鄭俊臣.紅壤坡地不同耕作影響水土流失的試驗(yàn)［J］.水土保持學(xué)報，1989，3(1):84-90.

［29］符素華，吳敬東，段淑懷，等.北京密云石匣小流域水土保持措施對土壤侵蝕的影響研究［J］.水土保持學(xué)報，2001，15(2):21-24.

［30］范建榮，王念忠，陳 光，等.東北地區(qū)水土保持措施因子研究［J］.中國水土保持科學(xué)，2011，9(3):75-78.

［31］焦菊英，王萬忠.人工草地在黃土高原水土保持中的減水減沙效益與有效蓋度［J］.草地學(xué)報，2001，9(3):177-182.

［32］謝 云，林 燕，張 巖.通用土壤流失方程的發(fā)展與應(yīng)用［J］.地理科學(xué)進(jìn)展，2003，22(3):279-287.

［33］楊學(xué)明，張曉平，方華軍.不同管理方式下吉林省農(nóng)田黑土流失量［J］.土壤通報，2003，34(5):389-393.

［34］汪東川，盧玉東.國外土壤侵蝕模型發(fā)展概述［J］.中國水土保持科學(xué)，，2004，2(2):35-40.