沈晨++萬(wàn)大娟++尹輝++楊龍++周慧++熊玲

摘要：對(duì)大通河流域的土壤侵蝕進(jìn)行估算和定量評(píng)價(jià)，以期為該流域的水土保持和生態(tài)環(huán)境保護(hù)工作提供科學(xué)合理的依據(jù)?；谛拚耐ㄓ猛寥懒魇Х匠蹋≧USLE），運(yùn)用遙感和地理信息技術(shù)，對(duì)大通河流域土壤侵蝕狀況進(jìn)行估算，并探討了不同土地利用類(lèi)型和不同坡度下土壤侵蝕的差異性。結(jié)果表明，大通河流域土壤侵蝕總面積為12 683.66 km2，占流域總面積的83.83%，年均侵蝕模數(shù)為6 274.76 t/km2。侵蝕強(qiáng)度整體表現(xiàn)為微度侵蝕和輕度侵蝕。土壤侵蝕主要集中在草地，侵蝕面積占比高達(dá)55.50%，灌木林次之。在0～25°坡度范圍內(nèi)，侵蝕面積隨著坡度的增大呈先增加后逐漸減少的趨勢(shì)，但在>25°時(shí)達(dá)到最大侵蝕面積。微度侵蝕主要發(fā)生在0～10°坡度范圍內(nèi)，中度以上侵蝕等級(jí)的侵蝕面積均在>25°坡度范圍時(shí)達(dá)到最大值。林草過(guò)牧濫伐、城鎮(zhèn)無(wú)序擴(kuò)張等不合理的人類(lèi)活動(dòng)是引起土壤侵蝕的主要原因，適度控制人類(lèi)活動(dòng)、提高植被覆蓋率，是流域內(nèi)今后土壤侵蝕治理及水土流失防治的重點(diǎn)。

關(guān)鍵詞：土壤侵蝕；修正的通用土壤流失方程；估算；大通河流域

中圖分類(lèi)號(hào)：S157.1 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：0439-8114（2017）08-1453-05

DOI：10.14088/j.cnki.issn0439-8114.2017.08.014

Study on Estimating Soil Erosion in Datong River Basin Based on RUSLE

SHEN Chen1，WAN Da-juan1，YIN Hui2，YANG Long3，ZHOU Hui2，XIONG Ling1

（1.Collge of Resources and Environment Science， Hunan Normal University， Changsha 410081，China；2.Department of Tourism，Huizhou University，Huizhou 516007，Guangdong，China；3.Qinghai Province Key Laboratory of Phycical Geography and Environment Process，Qinghai Normal University，Xining 810008，China）

Abstract： The soil erosion in Datong River basin was estimated and quantitatively evaluated in order to provide a scientific and rational basis for soil and water conservation and environmental protection of the basin. Based on the revised universal soil loss equation（RUSLE）， supported by remote sensing and geographical information system techniques， the situation of soil erosion in Datong River basin was estimated， and the differences in soil erosion under different land use types and different slope was discussed. The results showed that the erosion area in Datong River basin was 12 683.66 km2， accounting for 83.83% of the total basin area， and the annual mean soil erosion modulus was 6 274.76 t/km2. The erosion intensity was overall dominated by micro and slight degree erosion. Soil erosion was mainly concentrated in the grassland， whose erosion area accounting was up to 55.50%， followed by shrubbery. In the range of 0～25° slope， the soil erosion area increase with the increase of slope gradient until to a certain value， after which the erosion area decreased， but reaching the maximum value when the slope gradient greater than 25°. The micro degree erosion occurs mainly in the range of 0～10° slope， and the area of moderate erosion and above all reached the maximum value when the slope gradient greater than 25°. Irrational human activities， such as excessive deforestation， overgrazing and urban sprawl were the main cause of soil erosion. And to moderately control human activities and to improve vegetation coverage were the focus of future water loss and soil erosion prevention in Datong River basin.

Key words： soil erosion； revised universal soil loss equation； estimation； Datong River basin

土壤侵蝕是土壤及其母質(zhì)在水力、風(fēng)力、凍融、重力等外力作用下被破壞、剝蝕、搬運(yùn)和沉積的過(guò)程，是當(dāng)今人類(lèi)面臨的一種普遍、持續(xù)性強(qiáng)的地質(zhì)災(zāi)害[1]。流域內(nèi)的土壤侵蝕會(huì)帶來(lái)一系列嚴(yán)重后果，極易造成土地資源的破壞，增加河流輸沙量，造成河流水量和水質(zhì)的變化，最終引起生態(tài)環(huán)境惡化。大通河是黃河的二級(jí)支流，湟水河的一級(jí)支流，發(fā)揮著涵養(yǎng)水源及調(diào)控黃河源生態(tài)系統(tǒng)的重要作用。近年來(lái)，由于人類(lèi)活動(dòng)的強(qiáng)烈干擾，加之該流域獨(dú)特的土壤地理環(huán)境，導(dǎo)致其土壤侵蝕程度不斷加劇，迫切需要對(duì)流域內(nèi)土壤侵蝕進(jìn)行科學(xué)的估算和評(píng)價(jià)，為水土流失防治提供依據(jù)。

定量評(píng)估區(qū)域土壤侵蝕、分析土壤侵蝕空間分布特征是防治區(qū)域水土流失的基礎(chǔ)[2]。傳統(tǒng)的土壤侵蝕分析及預(yù)測(cè)方法需要大量的時(shí)間及人力物力投入，為有效開(kāi)展土壤侵蝕預(yù)測(cè)及水土保持防治工作，國(guó)內(nèi)外學(xué)者結(jié)合各自所在地區(qū)的實(shí)際情況提出了不同的預(yù)報(bào)模型，如美國(guó)的WEEP[3]（Water erosion prediction project）、歐洲的EUROSEM[4]（European soil erosion model）、荷蘭的LISEM[5]（Limburg soil erosion model）等。修正的通用土壤流失方程[6]（Revised universal soil loss equation，RUSLE）是基于通用土壤流失方程[7]（USLE）所提出，針對(duì)USLE模型中的不足加以重新評(píng)估，使其具有更高的準(zhǔn)確度和更強(qiáng)的實(shí)用性，是目前國(guó)內(nèi)外應(yīng)用最廣泛的土壤侵蝕模型。近年來(lái)，國(guó)內(nèi)學(xué)者將該模型廣泛地應(yīng)用在中國(guó)土壤侵蝕和水土流失的預(yù)測(cè)及治理工作中，取得了豐碩的成果[8，9]。在借鑒前人經(jīng)驗(yàn)的基礎(chǔ)上，基于修正的通用土壤流失方程（RUSLE），依據(jù)大通河流域的現(xiàn)狀選擇合理的方法計(jì)算RUSLE中各個(gè)因子值，運(yùn)用遙感和地理信息技術(shù)，對(duì)大通河流域的土壤侵蝕進(jìn)行估算和定量評(píng)價(jià)，以期為該流域的水土保持和生態(tài)環(huán)境保護(hù)工作提供科學(xué)合理的依據(jù)。

1 研究區(qū)概況與數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

大通河位于青藏高原東北部，發(fā)源于青海省天峻縣托勒南山，流域范圍介于98°30′-103°15′E， 36°30′-38°25′N(xiāo)之間，流域面積15 130 km2，干流全長(zhǎng)560.7 km，自西北向東南依次流經(jīng)青海省的天俊、祁連等11個(gè)縣（區(qū)），與湟水匯合后最終注入黃河，是華熱藏族地區(qū)最大的河流之一。大通河流域深居西北內(nèi)陸，冬長(zhǎng)夏短，氣候嚴(yán)寒，冰凍期較長(zhǎng)，屬典型的大陸高原性氣候，年平均降水量300～600 mm，年內(nèi)降水主要集中在6～9月，占年降水量的70%左右，10月至翌年5月為旱季[10]。地勢(shì)總體西北高、東南低，海拔在2 122～5 232 m之間；流域內(nèi)的地貌類(lèi)型主要有冰蝕構(gòu)造高山、侵蝕構(gòu)造中山、構(gòu)造剝蝕低山丘陵、冰川冰水堆積臺(tái)地及堆積平原[11]；主要土地利用類(lèi)型有建筑用地、林地、灌木林、草地、旱地、裸地、水域等，主要土壤類(lèi)型有高山石質(zhì)土、山地棕褐土、黑鈣土、栗鈣土、草甸土、沼澤土等。

1.2 數(shù)據(jù)來(lái)源

采用的數(shù)據(jù)有1993～2012年大通河流域及其周邊共5個(gè)氣象站點(diǎn)的月降雨量，來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)氣候資料月值數(shù)據(jù)集；大通河流域30 m×30 m分辨率數(shù)字高程模型（DEM）；2011年美國(guó)陸地衛(wèi)星Landsat/TM影像，空間分辨率為30 m×30 m；土壤數(shù)據(jù)選自青海省土壤圖，來(lái)自第二次土壤普查；土地利用數(shù)據(jù)選自青海省2010年土地利用現(xiàn)狀圖。

2 研究方法

對(duì)大通河流域土壤侵蝕量估算研究采用修正的通用土壤流失方程（RUSLE），該方程全面地考慮了自然要素對(duì)土壤侵蝕的作用機(jī)制，并綜合體現(xiàn)在降雨侵蝕力、土壤可蝕性等多個(gè)不同的因子上，具有較強(qiáng)的科學(xué)性和廣泛的實(shí)用性。其一般表達(dá)式為：

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

式中，A為年土壤侵蝕量（t/km2）；R為年降雨侵蝕力因子（MJ·mm/hm2·h）；K為土壤可蝕性因子[t·h/（MJ·mm）]；L為坡長(zhǎng)因子；S為坡度因子；C為植被覆蓋與管理因子；P為水土保持措施因子（L、S、C、P均無(wú)量綱）。

2.1 降雨侵蝕力因子

降雨侵蝕力是降水及其產(chǎn)生的徑流所具有的引起土壤侵蝕的潛在能力，是土壤的水蝕動(dòng)力[12]，其因子（R）與降雨量、降雨強(qiáng)度、降雨歷時(shí)、雨滴的大小以及雨滴下降速度有關(guān)，但由于降雨過(guò)程等資料實(shí)地監(jiān)測(cè)和數(shù)據(jù)獲取難度較大，因而國(guó)內(nèi)外學(xué)者提出了各種基于降雨量指標(biāo)計(jì)算R的簡(jiǎn)易算法[13，14]。采用FAO通過(guò)改進(jìn)Fournier[15]提出的流域降雨侵蝕力指數(shù)的方法來(lái)計(jì)算R，該方法綜合考慮了年降水量和逐月降水量，且數(shù)據(jù)獲取較為容易。其表達(dá)式如下：

F=■■

式中，pi為第i月降雨量，p為年降雨量。建立R與該指數(shù)的線性關(guān)系式：

R=a×F+b （1）

式中，R為降雨侵蝕力因子，a、b為常數(shù)，其數(shù)值取決于氣候條件。根據(jù)大通河流域的氣候條件以及國(guó)內(nèi)外其他相似地區(qū)的類(lèi)比分析，確定a與b的取值分別為4.17和-152。選取大通河流域及其周邊5個(gè)站點(diǎn)1993～2012年的降雨數(shù)據(jù)，通過(guò)各個(gè)站點(diǎn)的年降雨量及各月降雨量，利用關(guān)系式（1）計(jì)算R。

根據(jù)多年平均R計(jì)算結(jié)果，通過(guò)對(duì)不同插值法的比較分析，在ArcGIS中選取精度較高的反距離權(quán)重法進(jìn)行插值，得到流域降雨侵蝕力因子的空間分布圖。

3 結(jié)果與分析

3.1 大通河流域土壤侵蝕現(xiàn)狀分析

由表3可知，大通河流域土壤侵蝕總面積為12 683.66 km2，占流域總面積的83.83%，年均侵蝕模數(shù)為6 274.76 t/km2。從輕度侵蝕到劇烈侵蝕， 其面積分別為4 914.48、3 229.24、1 450.57、872.82、952.11、1 264.44 km2，占流域侵蝕總面積比分別為38.75%、25.46%、11.44%、6.87%、7.51%、9.97%，侵蝕面積先逐漸減少，后有小幅度增加。流域內(nèi)侵蝕整體表現(xiàn)為以微度侵蝕和輕度侵蝕為主，二者侵蝕面積占比高達(dá)64.21%；中度等級(jí)以上侵蝕面積共4 539.94 km2，占侵蝕總面積比例為35.79%，其中劇烈侵蝕主要集中在流域下游及環(huán)流域邊緣的裸地分布區(qū)，侵蝕面積所占比例為9.97%。

3.2 大通河流域土壤侵蝕隨土地利用類(lèi)型的分布特征

不同土地利用類(lèi)型對(duì)土壤侵蝕的影響主要表現(xiàn)在植被覆蓋度、植被類(lèi)型及水土保持措施等因子的差異上。根據(jù)大通河流域不同土地利用類(lèi)型下土壤侵蝕面積統(tǒng)計(jì)結(jié)果（表4），土壤侵蝕主要集中在草地，侵蝕面積為7 039.44 km2，占流域侵蝕面積比高達(dá)55.50%，從大通河流域土地利用構(gòu)成上看，草地是流域內(nèi)主要的土地利用類(lèi)型，加之盲目挖金采藥、過(guò)度放牧導(dǎo)致的草場(chǎng)退化，可能是造成流域內(nèi)土壤侵蝕面積集中在草地的原因。灌木林侵蝕面積為2 512.34 km2，僅次于草地，占侵蝕總面積的19.81%。雖然侵蝕面積大，但由于灌木林和草地植被覆蓋度高，因此土壤侵蝕總量并不大。建筑用地和水域的侵蝕面積極少，僅分別占流域侵蝕總面積的0.45%、0.76%。裸地侵蝕面積占比為12.51%，由于植被覆蓋度低，缺乏完善的水土保持措施，該區(qū)域也是土壤侵蝕的主要地帶。

3.3 大通河流域土壤侵蝕隨坡度的分布特征

坡度是影響土壤侵蝕的主要因子之一，不同地形起伏下的土壤侵蝕面積和強(qiáng)度呈現(xiàn)出一定的差異性。根據(jù)大通河流域不同坡度下土壤侵蝕面積的統(tǒng)計(jì)（表5），在0～25°坡度范圍內(nèi)，隨著坡度的增大，各坡度范圍的侵蝕面積分別為1 704.90、2 652.56、 2 051.97、1 625.87、1 274.37 km2，侵蝕面積呈先增加后逐漸減少的趨勢(shì)，在>25°坡度范圍內(nèi)的侵蝕面積達(dá)到最大值3 369.75 km2。在0～10°坡度范圍內(nèi)，微度侵蝕到劇烈侵蝕的面積逐漸減小，在>15°坡度范圍內(nèi)，微度侵蝕到劇烈侵蝕的面積先逐漸減小，在強(qiáng)烈侵蝕處達(dá)到侵蝕面積最小值，后逐漸增大。從坡度與侵蝕等級(jí)的相關(guān)性上來(lái)看，微度侵蝕主要發(fā)生在0～10°坡度范圍內(nèi)，其侵蝕面積占微度侵蝕總面積的52.7%；中度以上侵蝕等級(jí)的侵蝕面積均在>25°坡度范圍內(nèi)達(dá)到最大值；隨著坡度的增大，劇烈侵蝕的面積隨之增加，在坡度>25°時(shí)達(dá)到劇烈侵蝕最大面積689.53 km2，占劇烈侵蝕總面積的53.5%。

4 小結(jié)與討論

在地理信息系統(tǒng)技術(shù)的支持下，基于修正的通用土壤流失方程（RUSLE），通過(guò)合理選取模型中各因子的計(jì)算方法，估算了大通河流域土壤侵蝕狀況，并探討了不同土地利用類(lèi)型和不同坡度下土壤侵蝕的差異性。主要結(jié)論如下：

1）總體上看，大通河流域土壤侵蝕總面積為12 683.66 km2，占流域總面積的83.83%，年均侵蝕模數(shù)為6 274.76 t/km2。侵蝕強(qiáng)度整體表現(xiàn)為以微度侵蝕和輕度侵蝕為主，劇烈侵蝕所占比例為9.97%。從大通河流域土壤侵蝕強(qiáng)度分級(jí)圖中可以看出，劇烈侵蝕主要集中在流域中下游及環(huán)流域邊緣的裸地分布區(qū)，其中不合理的人類(lèi)活動(dòng)諸如林草過(guò)牧濫伐、城鎮(zhèn)無(wú)序擴(kuò)張等，加之裸地植被覆蓋度低，缺乏有效的水土保持措施，是造成流域劇烈侵蝕的主要原因。賈俊姝等[22]關(guān)于大通縣1995、2005年土壤侵蝕量估算研究中所得年侵蝕模數(shù)分別為6 855、5 183 t/km2，相比較之下，本研究估算的年均侵蝕模數(shù)與其相差并不大，目前流域內(nèi)的侵蝕強(qiáng)度呈增強(qiáng)的趨勢(shì)。

2）從大通河流域土地利用結(jié)構(gòu)上看，草地是流域內(nèi)主要的土地類(lèi)型，土壤侵蝕也集中發(fā)生在草地，其侵蝕面積占比高達(dá)55.50%。雖然草地植被覆蓋度較高，但由于存在退化，加之過(guò)度放牧、濫挖亂采等現(xiàn)象，草地的土壤保持和涵養(yǎng)水源的能力也會(huì)逐漸降低，出現(xiàn)大面積的土壤侵蝕和水土流失，侵蝕等級(jí)也會(huì)逐漸增加。因此，在流域內(nèi)通過(guò)封禁修復(fù)、植樹(shù)種草等措施提高植被覆蓋度，減輕人類(lèi)活動(dòng)的不利影響，是區(qū)域土壤侵蝕防治的重點(diǎn)。

3）從不同坡度上的土壤侵蝕分異特征來(lái)看，微度侵蝕和輕度侵蝕主要發(fā)生在0～10°坡度范圍內(nèi)，各坡度范圍的侵蝕總面積在>25°時(shí)達(dá)到最大值，中度以上侵蝕等級(jí)的侵蝕面積也均在>25°時(shí)達(dá)到最大值，結(jié)合其他地區(qū)的研究結(jié)果，坡度在25°以上的區(qū)域是土壤侵蝕易發(fā)生的地區(qū)，故針對(duì)>25°坡度范圍的土壤侵蝕防治也是流域內(nèi)土壤侵蝕預(yù)防及治理工作的重點(diǎn)。

參考文獻(xiàn)：

[1] 史德明.如何正確理解有關(guān)水土保持術(shù)語(yǔ)的討論[J].土壤侵蝕與水土保持學(xué)報(bào)，1998，4（4）：89-91.

[2] 李占斌，朱冰冰，李 鵬.土壤侵蝕與水土保持研究進(jìn)展[J].土壤學(xué)報(bào)，2008，45（5）：802-809.

[3] NEARING M A，F(xiàn)OSTER G R，LANE L J. A process-based soil erosion model for USDA-Water Erosion Prediction Project Technology[J].Trans ASAE，1989（32）：1587-1593.

[4] MORGAN R. The European Soil Erosion Model： An Update on Its Structure and Research Base， In： RICKSON R（ED.） Conserving Soil Resources： European Perspectives[M].Cambridge：CAB International，1994.

[5] DE ROO A，WESSELING C G，RITSMA C G. LISEM： A single-event， physically based hydrological and soils erosion model for drainage basin.Ⅰ：Theory， input and output[J].Hydrological Processes，1996（10）：1107-1117.

[6] RENARD K G，F(xiàn)OSTER G R，WEESIES G A，et al. Predicting Soil Erosion by Water： A Guide to Conservation Planning with the Revised Universal Soil Loss Equation（RUSLE）[M].Washington DC：United States Department of Agriculture，1997.

[7] WISCHMEIER W H，SMITH D D. Predicting Rainfall Erosion Losses from Cropland East of the Rocky Mountains：Guide for Selection of Practices for Soil and Water Conservation[M].Washington DC：US Government Printing Office，1965.

[8] 方廣玲，香 寶，趙 衛(wèi)，等.基于GIS和RUSLE的拉薩河流域土壤侵蝕研究[J].水土保持學(xué)報(bào)，2015，29（3）：6-12.

[9] 趙明松，李德成，張甘霖.基于RUSLE模型的安徽省土壤侵蝕及其養(yǎng)分流失評(píng)估[J].土壤學(xué)報(bào)，2016，53（1）：28-38.

[10] 劉賽艷，黃 強(qiáng)，劉登峰，等.青藏高原大通河流域氣候要素時(shí)空演變特性分析[J].水資源與水工程學(xué)報(bào)，2015，26（3）：24-29，34.

[11] 吳恒卿，劉賽艷，黃 強(qiáng)，等.基于SWAT模型的大通河流域徑流模擬[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，43（9）：210-216.

[12] 劉寶元，謝 云，張科利.土壤侵蝕預(yù)報(bào)模型[M].北京：中國(guó)科學(xué)技術(shù)出版社，2001.

[13] RICHARDSON C W，F(xiàn)OSTER G R，WRIGHT D A. Estimation of erosion index from daily rainfall amount[J].Transactions of the ASAE，1983，26（1）：153-156.

[14] 伍育鵬，謝 云，章文波.國(guó)內(nèi)外降雨侵蝕力簡(jiǎn)易計(jì)算方法的比較[J].水土保持學(xué)報(bào)，2001，15（3）：31-34.

[15] FOURNIER F. Climat et Erosion[M].Paris：Presses Universitarires de France，1960.

[16] 劉吉峰，李世杰，秦寧生，等.青海湖流域土壤可蝕性K值研究[J].干旱區(qū)地理，2006，29（3）：321-326.

[17] MCCOOL D K，F(xiàn)OSTER G R，MUTCHLER C K，et al. Revised slope length factor for the universal soil loss equation[J].Transactions of the American Society of Agricultural Engineers，1989，32：1571-1576.

[18] LIU B Y，NEARING M A，RISSE L M. Slope gradient effects on soil loss for steep slopes[J].Transactions of the American Society of Agricultural Engineers，1994，37：1835-1840.

[19] WISCHMEIER W H，SMITH D D. Predicting Rainfall Erosion Losses：A Guide to Conservation Planting[M].United States：Science and Education Administration，1978.

[20] 蔡崇法，丁樹(shù)文，史志華，等.應(yīng)用USLE模型與地理信息系統(tǒng)IDRISI預(yù)測(cè)小流域土壤侵蝕量的研究[J].水土保持學(xué)報(bào)，2000， 14（2）：19-24.

[21] 楊 娟，葛劍平，李慶斌.基于GIS和USLE的臥龍地區(qū)小流域土壤侵蝕預(yù)報(bào)[J].清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，46（9）：1526-1529.

[22] 賈俊姝，陳金蓮，高 欣，等.基于RUSLE的大通縣土壤侵蝕量估算[J].西北林學(xué)院學(xué)報(bào)，2012，27（2）：56-61.