劉斌濤，宋春風(fēng)，史 展，陶和平

(1.中國科學(xué)院 水利部成都山地災(zāi)害與環(huán)境研究所，成都 610041；2.西南大學(xué) 地理科學(xué)學(xué)院，重慶 400715)

?

蘆山地震災(zāi)區(qū)土壤流失方程研究

劉斌濤1，宋春風(fēng)1，史 展2，陶和平1

(1.中國科學(xué)院 水利部成都山地災(zāi)害與環(huán)境研究所，成都 610041；2.西南大學(xué) 地理科學(xué)學(xué)院，重慶 400715)

蘆山地震誘發(fā)大量崩塌、滑坡、泥石流等次生地質(zhì)災(zāi)害，植被破壞嚴(yán)重，生態(tài)環(huán)境受到嚴(yán)重威脅。建立土壤流失方程是進(jìn)行區(qū)域土壤侵蝕定量評估的最有效方法。通過對蘆山地震水土流失調(diào)查與評價工作的總結(jié)與提煉，介紹了蘆山地震災(zāi)區(qū)土壤流失方程的基本形式，給出了降雨侵蝕力因子(R值)、土壤可蝕性因子(K值)、坡長坡度因子(LS值)、地表覆蓋與管理因子(C值)、水土保持措施因子(P值)的具體計算方法。本方程對蘆山地震災(zāi)區(qū)水土流失調(diào)查與動態(tài)監(jiān)測、水土保持規(guī)劃、災(zāi)后生態(tài)恢復(fù)重建具有參考意義。

水土保持；蘆山地震災(zāi)區(qū)；土壤流失；水利普查；降雨侵蝕力

1 研究背景

2013年4月20日，四川省雅安市蘆山縣境內(nèi)發(fā)生Ms7.0級強烈地震，這是四川省自2008年汶川特大地震以來發(fā)生的第2次嚴(yán)重地震災(zāi)害，強震波及四川雅安、成都、樂山、眉山、甘孜、涼山、德陽等市州的32個縣(市、區(qū))，受災(zāi)人口約218.4萬人，遇難196人，失蹤2人，受傷14 785人。蘆山地震發(fā)生后，筆者負(fù)責(zé)了“蘆山地震災(zāi)區(qū)水土流失調(diào)查與評價”應(yīng)急項目。經(jīng)過1 a的調(diào)查和分析，獲取了蘆山地震災(zāi)區(qū)1∶50 000、重災(zāi)區(qū)1∶10 000的土壤侵蝕強度數(shù)據(jù)，掌握了蘆山地震災(zāi)區(qū)震前、震后土壤侵蝕的面積、強度及空間分布特征，編制了蘆山地震災(zāi)區(qū)分縣(市、區(qū))和重災(zāi)區(qū)分鄉(xiāng)鎮(zhèn)的土壤侵蝕強度圖，支撐了《“4·20”蘆山強烈地震水土保持專題公報》的及時發(fā)布?！疤J山地震災(zāi)區(qū)水土流失調(diào)查與評價”工作是在第一次全國水利普查水土保持情況普查、四川省水土保持情況普查工作基礎(chǔ)上開展的。蘆山地震災(zāi)區(qū)水土流失調(diào)查綜合使用了遙感三參量調(diào)查法[1]、土壤流失方程計算法和空間抽樣調(diào)查法[2-3]，多種方法相互比較，取得優(yōu)良的調(diào)查精度。本文是對蘆山地震災(zāi)區(qū)水土流失調(diào)查工作中土壤流失方程計算方法的總結(jié)，可以為我國區(qū)域水土流失調(diào)查，尤其是西南土石山區(qū)水土流失調(diào)查工作提供借鑒和參考。

2 方程基本形式

Wischmeier等[4-5]于20世紀(jì)60年代提出了通用土壤流失方程(Universal Soil Loss Equation，USLE)模型。經(jīng)過50多年的發(fā)展，USLE模型已經(jīng)成為全世界廣泛應(yīng)用的土壤流失預(yù)報方程[6-10]。USLE模型在我國應(yīng)用十分廣泛，并出現(xiàn)了眾多國產(chǎn)化方程，如中國土壤流失方程(CSLE)[2, 11-12]、中國坡面水蝕預(yù)報模型[13-14]、云南省金沙江流域土壤流失方程[15-16]、黑龍江省土壤流失方程[17-18]、北京土壤流失方程[19-20]等。蘆山地震災(zāi)區(qū)土壤流失方程(Soil Loss Equation for the Area Affected by 4·20 Lushan Earthquake, LSSLE)仍采用USLE模型的基本形式，其公式為

A=R×K×L×S×C×P。

(1)

式中：A為土壤侵蝕模數(shù)[t/(km2·a)]；R為降雨侵蝕力因子[MJ·mm/(hm2·h·a)]；K為土壤可蝕性因子[t·km2·h/(hm2·MJ·mm)]；L為坡長因子(無量綱，0～1)；S為坡度因子(無量綱，0～1)；C為地表覆蓋與管理因子(無量綱，0～1)；P為水土保持措施因子(無量綱，0～1)。

3 影響因子計算方法

3.1 降雨侵蝕力因子(R值)

降雨侵蝕力(rainfall erosivity)是指降雨引起土壤侵蝕的潛在能力[21]。本研究采用第一次全國水利普查水土保持情況普查使用的降雨侵蝕力算法計算蘆山地震災(zāi)區(qū)降雨侵蝕力數(shù)據(jù)，其形式[2, 22-24]為

(2)

其中：

(3)

(4)

(5)

(6)

3.2 土壤可蝕性因子(K值)

環(huán)境政策綜合氣候(Environmental Policy-Integrated Climate，EPIC)模型是最常用的計算土壤可蝕性的方程之一，在我國應(yīng)用非常廣泛[25-30]。在蘆山地震災(zāi)區(qū)土壤流失方程中用EPIC模型計算土壤可蝕性K值。EPIC模型的公式為：

(7)

(8)

式中：Sa為砂粒(2～0.05 mm)含量(%)；Si為粉砂(0.05～0.002 mm)含量(%)；Cl為黏粒(<0.002 mm)含量(%)；C為有機碳含量(%)。

張科利等[27, 31]研究指出，EPIC模型計算的K值與中國各地區(qū)實測的K值相差較大，并提出K值修正公式為

K=-0.013 83+0.515 75KEPIC。

(9)

3.3 坡度坡長因子(LS值)

地形因子(LS)是決定土壤侵蝕量的關(guān)鍵性因子。USLE，RUSLE模型給出了LS因子的計算公式[4-5, 32]。我國學(xué)者針對LS因子的修正算法開展了大量研究，代表性的有劉寶元算法[33-34]、江忠善算法[13-14]、楊子生算法[15-16, 35]等。

楊子生[15-16, 35]在金沙江流域進(jìn)行了大量觀測，建立一個適用于金沙江流域陡坡環(huán)境的坡長(L)因子算法。蘆山地震災(zāi)區(qū)地形特征與金沙江流域較為接近，因此本文使用楊子生算法計算蘆山地震災(zāi)區(qū)L因子，其公式為

(10)

筆者整理分析了西南土石山區(qū)的徑流小區(qū)資料，提出了10°～25°及≥25°兩個坡度范圍的S因子計算算法，修正了現(xiàn)有的劉寶元算法[33, 36]。修正的S因子計算公式為

(11)

3.4 地表覆蓋與管理因子(C值)

本文將地表覆蓋類型分為農(nóng)地、林地、草地3類。農(nóng)地C值與農(nóng)作物類型、種植制度、耕地管理方式等有關(guān)。根據(jù)我國土地利用分類制度，本文將耕地分為水田、旱地2類分別確定C值。第一次全國水利普查水力侵蝕野外調(diào)查確定了每一個調(diào)查地塊的農(nóng)作物類型、植被覆蓋度、農(nóng)作物種植制度等信息。根據(jù)楊子生等[37-43]的研究成果，確定該區(qū)域水稻、小麥、玉米、黃豆、甘薯、馬鈴薯、油菜、蔬菜等作物生長期內(nèi)的C值分別為0.010 0，0.434 5，0.351 7，0.355 1，0.386 4，0.372 1，0.427 3和0.500 0。然后通過野外調(diào)查獲得的種植制度和種植結(jié)構(gòu)信息，進(jìn)行空間統(tǒng)計，確定各縣域耕地的C值。圖1展示了蘆山地震災(zāi)區(qū)21個縣(市、區(qū))水田和旱地的C值。

林地C值與林地郁閉度和林下地表覆蓋度密切相關(guān)。在USLE模型手冊[5]中，將林地劃分為：①未擾動林地；②放牧、火燒或采伐過的林地；③采伐后進(jìn)行了跡地整理準(zhǔn)備再造林的林地。我國現(xiàn)行的土地利用分類系統(tǒng)無法支撐這樣的處理方式。本文將林地分為茂密林地、有林地、灌木林地、其他林地等4類，分別確定它們的C值。根據(jù)USLE模型手冊[5]以及王萬忠等[44]、楊子生等[16]的研究，本文分別確定4種林地的C值賦值標(biāo)準(zhǔn)。根據(jù)第一次全國水利普查水力侵蝕野外調(diào)查單元中的林地郁閉度、地表覆蓋度調(diào)查參數(shù)，通過空間統(tǒng)計獲取了蘆山地震災(zāi)區(qū)縣域林地的C值。圖2展示了蘆山地震災(zāi)區(qū)21個縣(市、區(qū))茂密林地、有林地、灌木林地和其它林地的C值。

圖1 蘆山地震災(zāi)區(qū)耕地C值Fig.1 Value C of cultivated land in Lushan earthquake area

圖2 蘆山地震災(zāi)區(qū)林地C值Fig.2 Value of C of forest in Lushan earthquake area

草地C值主要取決于植被覆蓋度。江忠善等[13-14]提出了草地C值與植被覆蓋度的定量關(guān)系，其公式為

(12)

式中：C為草地C值；V為草地植被覆蓋度(%)。

本文使用HJ-1A/1B多光譜影像和蘆山地震災(zāi)區(qū)1∶50 000土地利用圖，借助像元二分法[45]計算了蘆山地震災(zāi)區(qū)2010年的植被覆蓋度數(shù)據(jù)，然后利用式(12)計算了全區(qū)草地的C值。

3.5 水土保持措施因子(P值)

蘆山地震災(zāi)區(qū)主要的水土保持措施有水平梯田(地)、坡式梯田(地)等類型。其中，在我國南方地區(qū)種植水稻的耕地稱為田，對應(yīng)的水土保持措施為水平梯田、坡式梯田；種植旱作作物的耕地稱為地，相應(yīng)的水土保持措施為水平梯地、坡式梯地。我國學(xué)者對水平梯田(地)、坡式梯田(地)的減蝕效益開展了大量研究，通過前人研究成果，本文確定該區(qū)域水田、水旱輪作田(地)、水平梯地、園地梯地、坡式梯田(地)的P值分別為0.010 0，0.052 8，0.136 2，0.103 5，0.187 9。利用第一次全國水利普查耕地野外調(diào)查資料，通過空間統(tǒng)計獲得每一個縣域耕地的P值。圖3給出了蘆山地震災(zāi)區(qū)21個縣(市、區(qū))旱地的P值。從圖3中可以看出，地處四川盆地西南部臺地、淺丘區(qū)的夾江、東坡、丹棱、雨城、名山等縣(區(qū))的P值較低，這與當(dāng)?shù)卮罅糠植嫉乃递喿魈?地)以及農(nóng)林復(fù)合經(jīng)營方式有關(guān)；盆周山區(qū)各縣P值普遍較高，這主要是因為這些地區(qū)坡耕地較多的緣故。

圖3 蘆山地震災(zāi)區(qū)旱地P值Fig.3 Value of P of dry land in Lushan earthquake area

4 結(jié) 論

通過對“蘆山地震災(zāi)區(qū)水土流失調(diào)查與評價”工作的總結(jié)與提煉，筆者建立了蘆山地震災(zāi)區(qū)土壤流失方程(LSSLE)，給出了降雨侵蝕力因子(R值)、土壤可蝕性因子(K值)、坡長坡度因子(LS值)、地表覆蓋與管理因子(C值)、水土保持措施因子(P值)的具體計算方法。該方程對于指導(dǎo)蘆山地震災(zāi)區(qū)水土流失調(diào)查與監(jiān)測、水土保持規(guī)劃和災(zāi)后生態(tài)恢復(fù)重建具有重要參考意義。

筆者利用蘆山地震災(zāi)區(qū)土壤流失方程(LSSLE)已經(jīng)計算獲得了災(zāi)區(qū)1∶50 000和重災(zāi)區(qū)1∶10 000的土壤侵蝕空間數(shù)據(jù)，編制了災(zāi)區(qū)分縣(市、區(qū))和重災(zāi)區(qū)分鄉(xiāng)鎮(zhèn)的土壤侵蝕強度圖，統(tǒng)計獲取災(zāi)區(qū)分鄉(xiāng)鎮(zhèn)的土壤侵蝕模數(shù)、年土壤侵蝕量、土壤侵蝕面積等數(shù)據(jù)，充分展示方程的應(yīng)用前景，特別是在大比例尺土壤侵蝕制圖方面的應(yīng)用能力。

筆者在后期工作中，對蘆山地震災(zāi)區(qū)土壤流失方程(LSSLE)進(jìn)行了發(fā)展和完善，將模型拓展成四川省土壤流失方程(Sichuan Soil Loss Equation，SCSLE)，實現(xiàn)對四川全省1∶50 000空間精度的土壤侵蝕調(diào)查與制圖。限于篇幅，本文只給出了蘆山地震災(zāi)區(qū)土壤流失方程(LSSLE)的基本形式和各個因子的計算方法，后續(xù)將進(jìn)一步對方程的精度驗證、蘆山地震災(zāi)區(qū)土壤侵蝕空間分布特征、蘆山地震誘發(fā)的土壤侵蝕以及基于無人機遙感和LSSLE的大比例尺(1∶10 000，1∶1 000)的土壤侵蝕制圖技術(shù)進(jìn)行專門分析。

[1] SL190—2007，土壤侵蝕分類分級標(biāo)準(zhǔn)[S]. 北京：中國水利水電出版社, 2008.

[2] 李智廣，符素華，劉寶元. 我國水力侵蝕抽樣調(diào)查方法[J]. 中國水土保持科學(xué)，2012, 10(1): 77-81.

[3] 劉寶元，郭索彥，李智廣，等. 中國水力侵蝕抽樣調(diào)查[J]. 中國水土保持，2013,(10): 26-34. [4] WISCHMEIER W H, SMITH D D. Predicting Rainfall-Erosion Losses from Cropland East of the Rocky Mountains[R]. Washington, USA:United States Department of Agriculture, 1965.

[5] WISCHMEIER W H, SMITH D D. Predicting Rainfall Erosion Losses—A Guide to Conservation Planning[R]. Washington, USA:United States Department of Agriculture, 1978.

[6] GRIMM M, JONES R, MONTANARELLA L. Soil Erosion Risk in Europe[R].Ottawa,Canada:European Soil Bureau and Institute for Environment & Sustainability, 2002.

[7] WALL G J, COOTE D R, PRINGLE E A,etal. Revised Universal Soil Loss Equation for Application in Canada——A Handbook for Estimating Soil Loss from Water Erosion in Canada[R]. Ispra, Italia:Agriculture and Agri-Food Canada, 1997.

[8] BAGARELLO V, FERRO V, PAMPALONE V. A New Expression of the Slope Length Factor to Apply USLE-MM at Sparacia Experimental Area (Southern Italy)[J]. CATENA. 2013, 102: 21-26.

[9] KINNELL P I A. Event Soil Loss, Runoff and the Universal Soil Loss Equation Family of models: A Review[J]. Journal of Hydrology, 2010, 385(1/2/3/4): 384-397.

[10]MOEHANSYAH H, MAHESHWARI B L, ARMSTRONG J. Field Evaluation of Selected Soil Erosion Models for Catchment Management in Indonesia[J]. Biosystems Engineering，2004, 88(4): 491-506.

[11]章文波，劉寶元. 基于GIS的中國土壤侵蝕預(yù)報信息系統(tǒng)[J]. 水土保持學(xué)報，2003, 17(2): 89-92.

[12]祝 贏，章文波，劉素紅，等. 第一次全國水利普查侵蝕模數(shù)的批量計算方法——基于CSLE和GIS的土壤水蝕模數(shù)計算器設(shè)計與應(yīng)用[J]. 水土保持通報，2012, 32(5): 291-295.

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

[14]江忠善，鄭粉莉，武 敏. 中國坡面水蝕預(yù)報模型研究[J]. 泥沙研究，2005，(4): 1-6.

[15]楊子生. 滇東北山區(qū)坡耕地土壤流失方程研究[J]. 水土保持通報，1999, 19(1): 1-9.

[16]楊子生. 云南省金沙江流域土壤流失方程研究[J]. 山地學(xué)報,2002, 20(增1): 1-9.

[17]張憲奎. 黑龍江省土壤流失方程應(yīng)用簡介[J]. 國土與自然資源研究，1995，(1): 30-32.

[18]張憲奎，許靖華. 黑龍江省土壤流失方程的研究[J]. 水土保持通報，1992, 12(4): 1-9.

[19]畢小剛，段淑懷，李永貴，等. 北京山區(qū)土壤流失方程探討[J]. 中國水土保持科學(xué)，2006, 4(4): 6-13.

[20]劉寶元，畢小剛，付素華. 北京土壤流失方程[M]. 北京: 科學(xué)出版社, 2010.

[21]章文波，謝 云，劉寶元. 降雨侵蝕力研究進(jìn)展[J]. 水土保持學(xué)報, 2002, 16(5): 43-46.

[22]劉斌濤，陶和平，宋春風(fēng)，等. 1960—2009年中國降雨侵蝕力的時空變化趨勢[J]. 地理研究，2013, 32(2): 245-256.

[23]劉斌濤，陶和平，宋春風(fēng)，等. 我國西南山區(qū)降雨侵蝕力時空變化趨勢研究[J]. 地球科學(xué)進(jìn)展，2012, 27(5): 499-509.

[24]劉斌濤，陶和平，宋春風(fēng)，等. 基于重心模型的西南山區(qū)降雨侵蝕力年內(nèi)變化分析[J]. 農(nóng)業(yè)工程學(xué)報，2012, 28(21): 113-120.

[25]ZHANG K L, SHU A P, XU X L,etal. Soil erodibility and its Estimation for Agricultural Soils in China[J]. Journal of Arid Environments，2008, 72(6): 1002-1011.

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

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

[28]SUN W, SHAO Q, LIU J,etal. Assessing the Effects of Land Use and Topography on Soil Erosion on the Loess Plateau in China[J]. Catena，2014, 121: 151-163.

[29]ZHOU W, WU B. Assessment of Soil Erosion and Sediment Delivery Ratio Using Remote Sensing and GIS: A Case Study of Upstream Chaobaihe River Catchment, North China[J]. International Journal of Sediment Research，2008, 23(2): 167-173.

[30]RAO E, OUYANG Z, YU X,etal. Spatial patterns and impacts of Soil Conservation Service in China[J]. Geomorphology，2014, 207: 64-70.

[31]ZHANG K L, SHU A P, XU X L,etal. Soil Erodibility and Its Estimation for Agricultural Soils in China[J]. Journal of Arid Environments，2008, 72(6): 1002-1011.

[32]FOSTER G R, YODER D C, WEESIES G A,etal. Revised Universal Soil Loss Equation User’s Guide[R]. Washington, USA:United States Department of Agriculture, 2003.

[33]LIU B Y, NEARING M A, RISSE L M. Slope Gradient Effects on Soil Loss for Steep Slopes[J]. Transactions of the Asae，1994, 37(6): 1835-1840.

[34]LIU B Y, NEARING M A, SHI P J,etal. Slope Length Effects on Soil Loss for Steep Slopes[J]. Soil Science Society of America Journal，2000, 64(5): 1759-1763.

[35]楊子生. 滇東北山區(qū)坡耕地土壤侵蝕的地形因子[J]. 山地學(xué)報，1999, 17(增1): 17-19.

[36]SL 592—2012,水土保持遙感監(jiān)測技術(shù)規(guī)范[S]. 北京: 中國水利水電出版社, 2012.

[37]楊子生. 滇東北山區(qū)坡耕地土壤侵蝕的作物經(jīng)營因子[J]. 山地學(xué)報，1999, 17(增1): 20-22.

[38]楊子生. 云南省金沙江流域土壤流失方程研究[J]. 山地學(xué)報. 2002, 20(增1): 1-9.

[39]唐 寅，代 數(shù)，蔣光毅，等. 重慶市坡耕地植被覆蓋與管理因子C值計算與分析[J]. 水土保持學(xué)報，2010, 24(6): 53-59.

[40]陳正發(fā). 基于RUSLE模型的重慶市土壤流失方程研究[D]. 重慶: 西南大學(xué), 2011.

[41]史志華，蔡崇法，丁樹文，等. 基于GIS和RUSLE的小流域農(nóng)地水土保持規(guī)劃研究[J]. 農(nóng)業(yè)工程學(xué)報，2002, 18(4): 172-175.

[42]劉培娟. 魯中南中低山丘陵典型侵蝕區(qū)水土流失規(guī)律研究[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2008.

[43]劉培娟，馬文貴，楊吉華，等. 魯中南山區(qū)徑流小區(qū)不同坡度條件下四種植被的水土流失規(guī)律研究[J]. 水土保持研究. 2007, 14(6): 357-359.

[44]王萬忠，焦菊英. 中國的土壤侵蝕因子定量評價研究[J]. 水土保持通報, 1996, 16(5): 1-20.

[45]穆少杰，李建龍，陳奕兆，等. 2001—2010年內(nèi)蒙古植被覆蓋度時空變化特征[J]. 地理學(xué)報. 2012, 67(9): 1255-1268.

(編輯：黃 玲)

Soil Loss Equation of Lushan Earthquake Area

LIU Bin-tao1，SONG Chun-feng1，SHI Zhan2，TAO He-ping1

(1.Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu 610041, China； 2. School of Geography Science, Southwest University, Chongqing 400715, China)

Many secondary disasters such as collapse, landslide and debris flow were induced by Lushan earthquake. The vegetation was seriously destroyed and ecological environment was under serious threat. The most effective method for quantificationally evaluating the regional soil erosion is to establish the soil loss equation. Through summarizing the soil loss investigation and assessment results in Lushan earthquake area, we described the soil loss equation which had been used in Lushan earthquake area. The specific calculation methods of rainfall erosivity factor, soil erodibility factor, slope length and slope gradient factor, surface cover and management factor and the soil and water conservation measures factor were given in this paper. This equation can be referenced for the investigation and dynamic monitoring of soil and water loss, the soil and water conservation planning and the ecological reconstruction of Lushan earthquake area.

soil and water conservation；Lushan earthquake area；soil loss；water conservancy survey；rainfall erosivity

2014-08-08；

2014-09-12

國家重點基礎(chǔ)研究計劃(973)項目(2015CB452706)；國家自然科學(xué)基金項目(41201457)；地理空間信息工程國家測繪地理信息局重點實驗室經(jīng)費資助項目(201117)

劉斌濤(1984-)，男，山東德州人，助理研究員，碩士，主要從事冰緣地貌過程與生態(tài)過程、土壤侵蝕、山地垂直帶譜等領(lǐng)域的研究工作，(電話)15982060905(電子信箱)lbt609@163.com。

10.11988/ckyyb.20140674

2016,33(01):15-19

S157

A

1001-5485(2016)01-0015-05