基于USLE的北黃嶺河小流域土壤侵蝕量估算

高立新

(岫巖滿族自治縣水利事務(wù)服務(wù)中心,遼寧 岫巖 114300)

土壤侵蝕是土壤顆粒在水力、風(fēng)力等綜合作用下持續(xù)分離、運移的過程,土壤侵蝕發(fā)生后必將對土地資源造成破壞,引發(fā)土地資源退化、干旱等次生災(zāi)害。流域內(nèi)土壤侵蝕則會加大流域輸沙量,引起水量變化及水質(zhì)惡化,水生態(tài)功能衰減。小流域?qū)儆谝粋€完整的自然匯水單元,也是水土流失和土壤侵蝕發(fā)生、發(fā)展的基本單元,研究小流域土壤侵蝕過程、水土保持過程是水土保持研究的重要組成部分。USLE(Universal Soil Loss Equation,通用土壤流失方程)是美國農(nóng)業(yè)部最早提出的土壤侵蝕量經(jīng)驗預(yù)報公式,對于因降雨所引發(fā)的水動力土壤侵蝕預(yù)測較為適用,與其他方法相比,該模型結(jié)構(gòu)簡單,數(shù)據(jù)量少,結(jié)果真實可靠,故該模型自20世紀(jì)30年代提出后便先后應(yīng)用于土壤侵蝕試驗及區(qū)域土壤侵蝕量定量研究。我國于20世紀(jì)80年代引入該模型后,逐漸在不同區(qū)域土壤侵蝕預(yù)報中進(jìn)行了應(yīng)用。本文依托岫巖縣蘇子溝鎮(zhèn)北黃嶺河小流域,對流域土壤侵蝕量展開估算分析,為構(gòu)建該流域水生態(tài)功能分區(qū)以及流域水資源保護(hù)、合理利用提供依據(jù)。

1 流域概況

北黃嶺河小流域位于岫巖縣蘇子溝鎮(zhèn),涉及蘇子溝鎮(zhèn)的北黃嶺1個行政村,流域總面積為1935hm2,流域土地利用情況見表1。北黃嶺河小流域整體呈北高南低地勢,多丘陵,河川交錯,地形主要以山地和丘陵為主,夾雜小面積沖積平原和盆地,低山區(qū)在全縣總面積中約占78%以上。流域平均海拔為211.64m。小流域位于中溫帶氣候區(qū),夏季炎熱多雨,冬季寒冷干燥,多年平均年降水量為820mm,且從上游向下游呈遞增趨勢,降雨主要集中在每年的6—9月,期間降水量在全年降水量中的占比超70%,暴雨大多發(fā)生在7—8月。流域內(nèi)土壤類型以棕壤和山地草甸土為主:棕壤土多見于低山丘陵,土層厚度為0.5～1.5m;草甸土則多見于河流兩岸,土層較厚,為1.0～3.0m,土壤肥力高。

北黃嶺河小流域主要存在輕度水力侵蝕,水土流失面積為461hm2,在流域總面積中占比23.82%,土壤侵蝕較為嚴(yán)重。因受制于治理資金及管理模式,流域內(nèi)水土流失始終未得到系統(tǒng)治理,并呈日益嚴(yán)峻之態(tài)勢。流域內(nèi)山洪災(zāi)害頻繁發(fā)生,造成耕地沖毀、道路沖斷,嚴(yán)重影響農(nóng)業(yè)生產(chǎn)。洪水沖刷河道兩岸導(dǎo)致水土流失加劇,河道內(nèi)泥沙量持續(xù)增大,抬高河床,河道泄洪能力持續(xù)減弱;土壤養(yǎng)分流失,土壤肥力持續(xù)減弱;徑流沖刷帶走表土養(yǎng)分,土層變薄,土壤肥力下降,糧食產(chǎn)量降低。

2 USLE模型及因子

USLE是美國農(nóng)業(yè)部于20世紀(jì)60年代在大量觀測數(shù)據(jù)及人工模擬試驗的基礎(chǔ)上,所構(gòu)建的特定區(qū)域土壤侵蝕量估算的方程模型,該模型涵蓋影響特定區(qū)域水土流失的主要影響因子[1],公式如下:

A=RKLSCP

(1)

式中:A為所研究區(qū)域土壤侵蝕量,t/(hm2·a);R為年降雨侵蝕力因子,MJ·mm/(hm2·h·a);K為研究區(qū)土壤可侵蝕因子,t·hm2·h/(hm2·MJ·mm);L為坡長因子,無量綱;S為坡度因子,無量綱;C為管理因子,無量綱;P為水土保持措施因子,無量綱。

2.1 年降雨侵蝕力因子

降雨徑流及雨水擊濺作用對土壤侵蝕影響較大,所謂降雨侵蝕力就是降雨過程可能引起土壤侵蝕的潛力。根據(jù)降雨過程資料所得出的降雨侵蝕力顯然最為可靠,但相關(guān)資料較難獲取,故本文在應(yīng)用USLE模型預(yù)測土壤侵蝕量時,主要根據(jù)月降水量進(jìn)行降雨侵蝕力的簡便估算[2],公式如下:

(2)

式中:Pi為項目區(qū)月降水量,mm;其余參數(shù)意義同前。

將所獲得的北黃嶺河小流域1966—2021年的氣象資料代入式(2),便可得出小流域多年平均降雨侵蝕力因子為290.94MJ·mm/(hm2·h·a),因流域面積不大,可視為流域內(nèi)降雨均一分布,故整個流域內(nèi)年降雨侵蝕力因子取值一致。

2.2 土壤可侵蝕因子

該因子體現(xiàn)的是在降雨徑流沖刷、雨水擊濺作用下,土壤被分散、運移的可能程度[3],也是進(jìn)行項目區(qū)土壤侵蝕量估算的主要指標(biāo)之一。該流域土類以棕壤和山地草甸土為主,可以參考當(dāng)前學(xué)者對各類土壤可侵蝕因子的研究成果,得出北黃嶺河小流域棕壤和山地草甸土可侵蝕因子K,分別取0.27t·hm2·h/(hm2·MJ·mm)和0.268t·hm2·h/(hm2·MJ·mm)。

2.3 坡長和坡度因子

小流域坡長和坡度是其余條件不變的情況下,某一特定坡長和坡度下土壤侵蝕量與標(biāo)準(zhǔn)小區(qū)土壤侵蝕量之比,對土壤水力侵蝕影響較大。在實際運用中,兩個因子通常結(jié)合起來,統(tǒng)一作為項目區(qū)地形因子進(jìn)行考慮。

USLE模型中的坡長因子按照下式確定:

L=(λ/22.1)m

(3)

式中:λ為實際坡長,m;m為坡長指數(shù),應(yīng)根據(jù)允許坡度θ的大小及以下范圍取值:

(4)

依托北黃嶺河小流域1∶10000地形圖,采用ArcGIS軟件進(jìn)行等高線校正和矢量化,并構(gòu)建小流域柵格數(shù)字高程模型,再借助ArcGIS軟件中的Spatial Analysis模塊在柵格數(shù)字高程模型中提取坡度、坡長,進(jìn)而計算坡度和坡長因子[4]。其中,土地利用現(xiàn)狀圖主要根據(jù)實地考察繪制而成;土壤類型圖則依托全國土壤普查成果圖,并以土屬為成圖單元繪制而成;流域數(shù)字高程模型則基于流域等高線圖制成。

2.4 管理及水土保持措施因子

管理因子即項目區(qū)植被覆蓋后土壤流失總量與實際土壤流失量之比,一般在0～1之間取值,該因子大小主要體現(xiàn)的是植被覆蓋、輪作制度、管理措施的綜合效果。此處應(yīng)用Erdas軟件對北黃嶺河小流域遙感影像實施人工目視解譯,并根據(jù)目視解譯結(jié)果進(jìn)行判讀和校驗,最終生成小流域土地利用圖層,對流域土地利用類型進(jìn)行劃分,并根據(jù)土地利用現(xiàn)狀對不同利用類型管理因子賦值。具體見表2。

水土保持措施因子是采取修筑梯田、等高耕作等特定水保措施后區(qū)域土壤侵蝕量與順坡種植時土壤侵蝕量之比,是水土保持措施采取后對土壤流失量控制程度的體現(xiàn)[5],取值在0～1之間,取值為0意味著水土保持措施效果十分顯著且不存在土壤侵蝕,而取值為1則意味著未采取任何水土保持措施[6]。

根據(jù)流域管理局所提供的北黃嶺河小流域水土保持措施圖及所收集到的2010—2020年項目區(qū)水土流失綜合治理資料,結(jié)合相關(guān)文獻(xiàn),對流域不同水土保持措施因子取值進(jìn)行了合理估計,結(jié)果見表3。

表3 北黃嶺河小流域水土保持措施因子取值

3 土壤侵蝕量估算過程及結(jié)果

根據(jù)統(tǒng)計結(jié)果,北黃嶺河小流域年土壤侵蝕總量為12.1萬t,年土壤侵蝕模數(shù)為900t/(km2·a),基本上為輕度侵蝕。流域內(nèi)57.17%的面積土壤侵蝕模數(shù)在5.0t/(km2·a)以下,為輕度侵蝕;劇烈侵蝕、極強(qiáng)烈侵蝕、強(qiáng)烈侵蝕及中度侵蝕的面積占比分別為5.22%、10.0%、10.0%、17.61%。流域水土流失現(xiàn)狀結(jié)果見表4。值得注意的是,即便是再完善的經(jīng)驗性土壤侵蝕預(yù)報模型,仍面臨基礎(chǔ)數(shù)據(jù)較難獲取的難題[7],本文則根據(jù)水利行業(yè)《土壤侵蝕分類分級標(biāo)準(zhǔn)》(SL 190—2007),將小流域土壤侵蝕量估算基礎(chǔ)數(shù)據(jù)簡化成土地利用類型、坡度等簡便且能準(zhǔn)確獲取的指標(biāo)[8],從而快速得出流域土壤侵蝕量測算結(jié)果。該流域內(nèi)土壤侵蝕表現(xiàn)出較大的空間差異性,通過進(jìn)一步分析看出,造成這種空間差異的原因主要在于土地利用類型以及地形地貌起伏變化。

表4 北黃嶺河小流域水土流失現(xiàn)狀

3.1 土地利用類型

土地利用類型差異主要體現(xiàn)在北黃嶺河小流域植被類型及植被覆蓋度差異方面,通過全面統(tǒng)計流域內(nèi)不同土地利用類型所對應(yīng)的土壤侵蝕強(qiáng)度,得到具體土地類型下的土壤侵蝕模數(shù)均值及侵蝕總量[9],結(jié)果見表5。其中林地年侵蝕總量取值最大,在總侵蝕量中占比為95.87%,說明林地是北黃嶺河小流域現(xiàn)狀土地利用類型下最容易發(fā)生水土流失的類型,必須加強(qiáng)重點治理。本研究依托項目區(qū)2021年全年遙感影像數(shù)據(jù),流域內(nèi)耕地年侵蝕模數(shù)均值為159.44t/(km2·a),年侵蝕總量均值為1421100t,在總侵蝕量中占比為8.69%。疏林地、草地、工況用地、水域、其他用地土壤侵蝕程度均較小。

表5 北黃嶺河小流域不同土地利用類型土壤侵蝕結(jié)果

基于以上對不同土地利用類型土壤侵蝕結(jié)果的分析,有助于北黃嶺河小流域管理部門優(yōu)化土地利用結(jié)構(gòu),為流域土壤侵蝕及水土流失治理提供依據(jù)。

3.2 地形地貌坡度

地形地貌起伏變化對流域土壤侵蝕空間差異的影響特征在北黃嶺河小流域下游城鎮(zhèn)密集區(qū)域及中游低山丘陵區(qū)表現(xiàn)得尤為突出。流域內(nèi)不同坡度下土壤侵蝕模數(shù)統(tǒng)計分析結(jié)果見表6。該小流域內(nèi)不同坡度下土壤侵蝕結(jié)果較為接近,且整體上隨著坡度的增大,土壤侵蝕程度加劇[10]。坡度在5°以內(nèi)的區(qū)域在流域總面積中占1.19%,占比最小,且對應(yīng)的土壤侵蝕量占比也為最小;坡度在25°以上的區(qū)域在流域內(nèi)占比為41.86%,土壤侵蝕量占比為31.35%。

表6 北黃嶺河小流域不同坡度下土壤侵蝕結(jié)果

4 結(jié) 論

綜上所述,本文依托遙感及地理信息系統(tǒng)技術(shù),應(yīng)用USLE土壤流失量估算方程模型,對北黃嶺河小流域土壤侵蝕量進(jìn)行分析估算后得出以下結(jié)論:流域內(nèi)年土壤侵蝕模數(shù)為900t/(km2·a),年土壤侵蝕總量平均為163589740t,總土壤侵蝕面積為461hm2,在流域總面積中占比為23.82%。所得到的流域土壤侵蝕量預(yù)測結(jié)果與相關(guān)研究基本吻合,結(jié)果可靠度高。該流域土壤侵蝕表現(xiàn)出十分明顯的空間差異性,土壤侵蝕受土地利用類型及地形坡度等的影響較大,其中林地、疏林地和耕地土壤侵蝕量占比最大;坡度在25°以上的區(qū)域土壤侵蝕占比最大。為此,流域管理部門應(yīng)當(dāng)將以上土壤侵蝕嚴(yán)重區(qū)域作為下一步水土流失治理的重點。