高璐，王學(xué)鳳，曲本亮，曹永強

（1.遼寧師范大學(xué)城市與環(huán)境學(xué)院，遼寧大連 116029;2.中國水利水電科學(xué)研究院，北京 100038; 3.大連市環(huán)境科學(xué)設(shè)計研究院，遼寧大連 116023）

基于RUSLE模型的遼寧省土壤侵蝕定量研究

高璐1，王學(xué)鳳2，曲本亮3，曹永強1

（1.遼寧師范大學(xué)城市與環(huán)境學(xué)院，遼寧大連 116029;2.中國水利水電科學(xué)研究院，北京 100038; 3.大連市環(huán)境科學(xué)設(shè)計研究院，遼寧大連 116023）

以降雨量、土壤、遙感影像及土壤侵蝕普查數(shù)據(jù)等為基礎(chǔ)數(shù)據(jù)，運用GIS技術(shù)，結(jié)合修正的通用土壤流失方程（RUSLE）模型對遼寧省土壤侵蝕狀況進行研究，分析了土壤侵蝕的空間分布特征以及與不同土地利用類型的關(guān)系。研究結(jié)果表明:2011年遼寧省土壤侵蝕面積為459.36萬hm2，年均土壤流失量為1.64億t，年平均侵蝕模數(shù)為3 637.8 t/（km2·a），屬于中度侵蝕，土壤侵蝕面積呈現(xiàn)上升趨勢;遼西低山丘陵區(qū)的土壤侵蝕最為嚴重，侵蝕面積為191.99萬hm2，城市受人為不合理的生產(chǎn)活動影響嚴重，導(dǎo)致土壤侵蝕面積增加;林地依舊是遼寧省土壤侵蝕發(fā)生的最主要的土地類型，占侵蝕總面積的57.66%，占土地利用類型的38.20%。

土壤侵蝕;RUSLE模型;定量分析;遼寧省;土地利用類型

1 研究背景

土壤侵蝕是全球性的重大環(huán)境問題之一［1］，是限制當(dāng)今人類生存與發(fā)展的全球性環(huán)境災(zāi)害。土壤侵蝕會造成土壤養(yǎng)分流失、水質(zhì)污染、生產(chǎn)力下降、惡化生態(tài)和人類生存環(huán)境［2］，在國內(nèi)外受到長期重視，因此對土壤侵蝕進行調(diào)查、研究、監(jiān)測和評估具有重要的現(xiàn)實意義。

對于這一重大的環(huán)境問題，19世紀(jì)國外學(xué)者就開始對土壤侵蝕進行研究［3-5］，并且相繼建立了一系列的基于經(jīng)驗的土壤侵蝕模型［6］。RUSLE模型［7-10］是代表模型USLE的修正模型，修正后的公式結(jié)構(gòu)合理、參數(shù)代表性普遍;結(jié)合GIS和RS，具有很強的實用性和綜合能力，將成為土壤侵蝕研究工具和水土保持規(guī)劃的主要研究手段。O.Terranova［11］運用RUSLE模型和GIS對地中海區(qū)域進行土壤侵蝕風(fēng)險的模擬實驗;V.Prasannakumar［12］運用RUSLE模型對印度的喀拉拉邦區(qū)域進行土壤侵蝕風(fēng)險評估;徐麗芬［13］以中國渤海為例，采用RUSLE模型和信息傳播模型在不同降雨情況下對土壤侵蝕進行了風(fēng)險評估;陳思旭［14］、張旭群［15］、李天宏［16］等學(xué)者也在不同方面對土壤侵蝕進行了研究。本文采用高分辨率遙感數(shù)據(jù)以及空間技術(shù)等手段，結(jié)合RUSLE模型以及土壤侵蝕強度分級標(biāo)準(zhǔn)，對遼寧省2011年土壤侵蝕狀況進行定量分析。

2 區(qū)域概況及研究方法

2.1 研究區(qū)域概況

研究區(qū)位于我國東北地區(qū)的南部，處于118°53'E—125°46'E，38°43'N—43°26'N之間，全省總面積為1 463.162萬hm2。地勢大致為:自北向南、自東西兩側(cè)向中部傾斜，山地丘陵分列東西兩廂。屬于溫帶大陸性季風(fēng)氣候區(qū)，是東北地區(qū)降水量最多的省份，年降水量在600～1 100 mm之間。遼寧省森林覆蓋率為37%，面積約為5萬km2，中部為平原棕黃土地帶，西部為低山丘陵生草棕色森林土黃白土地帶，北部土質(zhì)結(jié)構(gòu)差，水土流失最嚴重。在2005年第4次土壤侵蝕遙感普查中，土壤侵蝕面積占全省總面積比例的28.91%［17］，土壤侵蝕情況不容樂觀，所造成的生態(tài)環(huán)境退化，會加劇貧困、影響社會穩(wěn)定、阻礙社會進步，并將在未來成為影響遼寧省生態(tài)和經(jīng)濟可持續(xù)發(fā)展的重大問題。

2.2 所用數(shù)據(jù)來源

所用數(shù)據(jù)為2011年遼寧省土壤侵蝕面積（分市）統(tǒng)計表、遼寧省第4次土壤侵蝕遙感普查數(shù)據(jù)、遼寧省高程圖、遼寧省各市1980—2011年逐日降雨數(shù)據(jù)、遼寧省1∶150萬土壤類型圖、遼寧省地形圖及2011年遼寧省土地利用資料。

2.3 研究方法

本文選用目前應(yīng)用最廣、適用范圍最大的修正的通用土壤流失方程（RUSLE）進行土壤侵蝕的量化，其計算公式為

式中:A為年均土壤侵蝕模數(shù)（t/（km2·a））;R為降雨侵蝕力因子（（MJ·mm）/（km2·h·a））;K為土壤可蝕性因子（（t·km2·h）/（km2·MJ·mm））; LS為坡長坡度因子，無量綱;C為地表植被覆蓋與管理因子，無量綱;P為水土保持措施因子，無量綱。將遙感數(shù)據(jù)和非遙感數(shù)據(jù)相結(jié)合，經(jīng)過各因子的算式，計算出各因子的值，最終將各因子逐年相乘獲得研究區(qū)內(nèi)土壤侵蝕模數(shù)。本文只對遼寧省水蝕狀況進行研究分析。

3 結(jié)果分析

3.1 遼寧省土壤侵蝕現(xiàn)狀分析

根據(jù)RUSLE模型的式（1）以及《遼寧省水利基礎(chǔ)信息集》［18］中所獲得的數(shù)據(jù)，遼寧省現(xiàn)有土壤侵蝕面積為459.36萬hm2，占全省總面積的31.39%，年均土壤流失量為1.64億t，年平均侵蝕模數(shù)為3 637.8 t/（km2·a），屬于中度侵蝕。與2005年第4次遙感普查結(jié)果（土壤侵蝕面積為422.95萬hm2，年均土壤流失量為1.18億t，平均土壤侵蝕模數(shù)為2 801 t/（km2·a））相比，侵蝕面積增加了36.41萬hm2，由原來的28.9%上升到31.39%，增加了2.49%，平均每年增加0.28%，說明遼寧省土壤侵蝕呈上升趨勢。

3.2 遼寧省侵蝕強度等級分析

在研究區(qū)內(nèi)，侵蝕強度等級以輕度侵蝕、中度侵蝕以及強烈侵蝕為主，伴有少量的極強烈侵蝕和劇烈侵蝕，詳見表1、表2。以2011年為例，輕度侵蝕面積為237.69萬hm2，占全省總面積的16.25%;中度侵蝕面積為121.23萬hm2，占全省總面積的8.29%;強烈侵蝕面積為64.57萬hm2，占全省總面積的4.41%;極強烈侵蝕面積為27.94萬hm2，占全省總面積的1.91%;劇烈侵蝕面積為7.93萬hm2，占全省總面積的0.54%。與2005年土壤侵蝕強度等級評價結(jié)果比較分析，遼寧省的侵蝕面積較大，由平均侵蝕模數(shù)評判土壤侵蝕強度，研究區(qū)內(nèi)依舊為中度侵蝕，但其受降水變化以及人為活動等因素的影響較大，導(dǎo)致了平均侵蝕模數(shù)和侵蝕面積的不斷增長。

表1 土壤侵蝕強度分級標(biāo)準(zhǔn)Table 1Classification of level of soil erosion

表2 遼寧省土壤侵蝕等級評價結(jié)果Table 2Evaluation results of level of soil erosion in Liaoning province

3.3 遼寧省土壤侵蝕空間差異性分析

3.3.1 各市土壤侵蝕現(xiàn)狀分析

遼寧省2005年和2011年各市土壤侵蝕面積對比如圖1。

圖1 2005年和2011年遼寧省各市土壤侵蝕面積對比Fig.1Comparison among areas of soil erosion of all cities in Liaoning province in 2005 and 2011

對比各城市土壤侵蝕分布情況，侵蝕面積最大的是朝陽市，面積為80.85萬hm2，盤錦市的侵蝕面積最小，為0.67萬hm2。強烈以下的侵蝕最嚴重的是朝陽市、丹東市和大連市，面積分別為77.53，49.27，48.3萬hm2;極強烈和劇烈侵蝕最嚴重的城市為丹東市，面積為6.12萬hm2。土壤侵蝕占各市總面積40%以上的有營口市、葫蘆島市、大連市和朝陽市;在30%～40%的有阜新市、丹東市和鞍山市;在20%～30%的有鐵嶺市、錦州市和本溪市;其余各市在20%以下。

遼寧省各市在近6 a內(nèi)，土壤侵蝕面積均出現(xiàn)增加的現(xiàn)象。侵蝕面積的增加，主要受自然因素和人為因素的影響，而人為不合理的生產(chǎn)活動是造成土壤加速侵蝕的主導(dǎo)因素。近幾年遼寧省經(jīng)濟飛速發(fā)展，開礦、改擴建電廠、新建高速公路以及城市基礎(chǔ)建設(shè)等各類開發(fā)建設(shè)項目大量增加，造成了由開發(fā)建設(shè)活動帶來的短時間內(nèi)侵蝕面積的增加。盲目地追求經(jīng)濟效益，從而破壞了生態(tài)平衡;同時由于管理上的疏忽與不當(dāng)，放任這些不合理的土地利用的存在，更加加劇了水土流失的嚴重性從而促使土壤侵蝕面積增加。

3.3.2 各區(qū)域土壤侵蝕現(xiàn)狀分析

遼寧省2005年和2011年各區(qū)域土壤侵蝕面積對比如圖2。

圖2 2005年和2011年各區(qū)域土壤侵蝕面積對比Fig.2Comparison among areas of soil erosion of regions in Liaoning province in 2005 and 2011

土壤侵蝕面積區(qū)域分布仍以遼西低山丘陵區(qū)的錦州、朝陽、阜新和葫蘆島4市的土壤侵蝕最為嚴重。該地區(qū)土壤侵蝕面積為191.99萬hm2，占全省總侵蝕面積的41.80%。其中輕度侵蝕面積為110.92萬hm2，中度侵蝕面積為49.49萬hm2，強烈侵蝕面積為22.22萬hm2，極強烈侵蝕面積為7.96萬hm2，劇烈侵蝕為1.40萬hm2。遼寧東部山區(qū)侵蝕面積為158.04萬hm2，占全省總侵蝕面積的34.40%，其中中度以上侵蝕面積為40.73萬hm2。遼寧中部丘陵漫崗區(qū)侵蝕面積為109.29萬hm2，占全省總侵蝕面積的23.80%，其中中度以上侵蝕面積為28.13萬hm2。

地形環(huán)境、地貌環(huán)境、土壤環(huán)境以及社會人文環(huán)境都是影響遼寧各區(qū)域土壤侵蝕面積增加的因素。遼西地區(qū)的輕度、中度侵蝕主要分布在旱地上，強烈、極強烈、劇烈侵蝕主要分布在中覆蓋度的草地、有林地和疏林地上。從圖2和表3中2005年侵蝕面積發(fā)現(xiàn)遼西研究區(qū)內(nèi)輕度侵蝕面積呈現(xiàn)減少的趨勢，而中度以上侵蝕面積呈現(xiàn)增長的趨勢，而且增長的幅度較大，所以更應(yīng)該對中度以上侵蝕進行重點治理，減少耕地的分布是減少遼西地區(qū)土壤侵蝕的主要措施。自然因素所導(dǎo)致的土壤侵蝕現(xiàn)狀很難改變，但是，人為因素導(dǎo)致的可以通過合理、適當(dāng)?shù)恼{(diào)整土地利用類型從而達到抑制土壤侵蝕的發(fā)生。

表3 遼西地區(qū)不同研究時間侵蝕面積對比Table 3Comparison of erosion areas of western Liaoning province in 2005 and 2011

4 土壤侵蝕和不同土地利用類型的關(guān)系

就土壤侵蝕而言，土地利用類型是造成水土流失的重要的因子，是土壤侵蝕的促進因素。從遼寧省2011年和2005年土壤侵蝕與土地利用的關(guān)系發(fā)現(xiàn)（表4），林地是侵蝕發(fā)生的最主要的土地類型，林地侵蝕占整個遼寧省侵蝕面積的57.66%（2011年）和27.10%（2005年），增長了30.56%。其次，耕地也是侵蝕發(fā)生的另一個主要的土地類型，占遼寧省侵蝕面積的30.62%（2011年）和36.40% （2005年），但相比較侵蝕面積有所減少。將2005年和2011年比較，林地的侵蝕面積增加得最多，為150.12萬hm2，同時草地、水域以及城鄉(xiāng)、工礦等建設(shè)用地的面積也有著不同程度的增加。充分表明在2005—2011年間，由于人口的激增、對糧食迫切的需求以及人類不同程度的開發(fā)活動等所造成的不同土地利用類型的土壤侵蝕面積，正在不斷增加。

表4 2005年和2011年遼寧省不同土地利用類型的土壤侵蝕面積Table 4Areas of soil erosion with different types of land use in Liaoning province in 2005 and 2011

由于經(jīng)濟建設(shè)的迅猛發(fā)展，各類開發(fā)項目對林地的需求也大幅度增加，特別是國家對基礎(chǔ)設(shè)施建設(shè)投入的進一步加大，對林地的需求也迅速膨脹，同時退耕還林和封山育林的一些措施也使得林地的面積增加，因此林地的侵蝕面積呈現(xiàn)相對增長的趨勢。研究范圍內(nèi)林地多質(zhì)量較差，結(jié)構(gòu)單一，均為水土保持較為脆弱的林地。同時人口的增加對蠶場生態(tài)環(huán)境的超負荷壓力而產(chǎn)生的蠶場沙化現(xiàn)象;各類建設(shè)項目的開礦、采石等對生態(tài)影響較大的林地經(jīng)濟項目，極易造成對環(huán)境以及生態(tài)系統(tǒng)的破壞，都是造成林地土壤侵蝕最根本的原因。因此要在現(xiàn)有林地的基礎(chǔ)上實施工程措施，加快水保林建設(shè)，對單一林地實施混交，合理利用開發(fā)林地，因地制宜地使用林地，科學(xué)有效地管理和保護林地。

5 結(jié)論

（1）遼寧省2011年土壤侵蝕面積為459.36萬hm2，年均土壤流失量為1.64億t，年平均侵蝕模數(shù)為3 637.8 t/（km2·a），屬于中度侵蝕，與2005年第4次遙感普查結(jié)果相比，侵蝕面積增加了36.41萬hm2，占比例由原來的28.9%上升到31.39%，增加了2.49%，平均每年增加0.28%，土壤侵蝕面積土壤侵蝕呈現(xiàn)上升趨勢。

（2）侵蝕面積最大的城市是朝陽市，面積為80.85萬hm2，而盤錦市是侵蝕面積最小的城市，為0.67萬hm2。城市受人為不合理的生產(chǎn)活動影響嚴重，從而造成土壤加速侵蝕。土壤侵蝕面積區(qū)域分布仍以遼西地區(qū)的土壤侵蝕最為嚴重，為191.99萬hm2;遼寧東部山區(qū)侵蝕面積為158.04萬hm2;遼寧中部丘陵漫崗區(qū)侵蝕面積為109.29萬hm2。地形環(huán)境、地貌環(huán)境、土壤環(huán)境以及社會人文環(huán)境都是影響遼寧省各區(qū)域土壤侵蝕面積增加的因素。

（3）林地是侵蝕發(fā)生的最主要的土地類型，全省有一半以上侵蝕面積發(fā)生在林地。研究范圍內(nèi)林地多質(zhì)量較差、頻現(xiàn)蠶場沙化現(xiàn)象和各類對生態(tài)環(huán)境有影響的建設(shè)項目，因此導(dǎo)致土壤侵蝕格外嚴重。因此要加強林地的基礎(chǔ)設(shè)施建設(shè)，合理利用開發(fā)林地，因地制宜地使用林地，科學(xué)有效地管理和保護林地。

［1］劉娜，王克林，張偉，等.土壤侵蝕及其評價、校驗方法研究進展［J］.中國農(nóng)學(xué)通報，2011，27（18）:1-6，（LIU Na，WANG Ke-lin，ZHANG Wei，et al.Progress on the Study of Soil Erosion，Evaluation and Validation［J］. Chinese Agricultural Science Bulletin，2011，27（18）:1-6.（in Chinese））

［2］張燕，張洪，彭補拙，等.不同土地利用方式下農(nóng)地土壤侵蝕與養(yǎng)分流失［J］.水土保持通報，2003，23（1）: 23-31.（ZHANG Yan，ZHANG Hong，PENG Bu-zhuo，et al.Soil Erosion and Nutrient Loss of Various Land Use Patterns［J］.Bulletin of Soil and Water Conservation，2003，23（1）:23-31.（in Chinese））

［3］吳發(fā)啟，張洪江.土壤侵蝕學(xué)［M］.北京:科學(xué)出版社，2012.（WU Fa-qi，ZHANG Hong-jiang.Soil Erosion［M］. Beijing:Science Press，2012.（in Chinese））

［4］楊勤科，李銳，曹明明.區(qū)域土壤侵蝕定量研究的國內(nèi)外進展［J］.地球科學(xué)進展，2006，21（8）:849-856. （YANG Qin-ke，LI Rui，CAO Ming-ming.Advances of Quantitative Assessment on Regional Soil Erosion［J］. Advantage in Earth Sciences，2006，21（8）:849-856. （in Chinese））

［5］YU Z W.Surface Interpolation from Irregularly Distributed Points Using Surface Splines with Fortran Program［J］. Computers＆Geosciences，2001，27（7）:877-882.

［6］張光輝.土壤侵蝕模型研究現(xiàn)狀及展望［J］.水科學(xué)進展，2002，13（3）:389-396.（ZHANG Guang-hui.Research Situation and Prospect of the Soil Erosion Model［J］.Advances in Water Science，2002，13（3）:389-396.（in Chinese））

［7］WISCHMEIER W H，SMITH D D.Predicting RainfallerosionLossesfromCroplandEastoftheRocky Mountains:A Guide for Selection of Practices for Soil and Water Conservation［R］.Washington:U.S.Department of Agriculture，1965.

［8］WISCHMEIER W H，SMITH D D.Predicting Rainfall E-rosion Losses:A Guide to Conservation Planning［R］. Washington:U.S.Department of Agriculture，1978.

［9］RENARD K G，F(xiàn)OSTER G R，WEEIES G A.Predicting Soil Erosion by Water:A Guide to Conservation Planning with the Revised Universal Soil Loss Equation（RUSLE）［R］.Washington:U.S.Department of Agriculture，1997.

［10］陳云明，劉國彬，鄭粉莉，等.RUSLE侵蝕模型的應(yīng)用及進展［J］.水土保持研究，2004，11（4）:80-83.（CHEN Yun-ming，LIU Guo-bin，ZHENG Fen-li，et al.Progress and Application on Soil Erosion Model of RUSLE［J］.Research of Soil and Water Conservation，2004，11（4）:80-83.（in Chinese））

［11］TERRANOVA O.Soil Erosion Risk Scenarios in the Mediterranean Environment Using RUSLE and GIS:An Application Model for Calabria（Southern Italy）［J］.Geomorphology，2009，112（3）:258-245.

［12］PRASANNAKUMAR V.Estimation of Soil Erosion Risk within A Small Mountainous Sub-watershed in Kerala，India，Using Revised Universal Soil Loss Equation（RUSLE） and Geo-information Technology［J］.Geoscience Frontiers，2012，3（2）:209-215.

［13］XU L F.Risk Assessment of Soil Erosion in Different Rainfall Scenarios by RUSLE Model Coupled with Information Diffusion Model:A Case Study of Bohai Rim，China［J］. Catena，2012（100）:74-82.

［14］陳思旭，楊小喚，肖林林，等.基于RUSLE模型的南方丘陵山區(qū)土壤侵蝕研究［J］.資源科學(xué)，2014，36（6）:1288-1297.（CHEN Si-xu，YANG Xiao-huan，XIAO Lin-lin，et al.Study of Soil Erosion in the Southern Hillside Area of China Based on RUSLE Model［J］.Resources Science，2014，36（6）:1288-1297.（in Chinese））

［15］張旭群，陳耀強，陳浩昆，等.基于GIS和RUSLE的粵東黃岡河流域土壤侵蝕評估［J］.中國水土保持，2013，（2）:34-37.（ZHANG Xu-qun，CHEN Yao-qiang，CHEN Hao-kun，et al.Evaluation on Soil Erosion of Huanggang River Basin in Eastern Guangdong Based on GIS and RUSLE［J］.Soil and Water Conservation in China，2013，（2）:34-37.（in Chinese））

［16］李天宏，鄭麗娜.基于RUSLE模型的延河流域2001-2010年土壤侵蝕動態(tài)變化［J］.自然資源學(xué)報，2012，27（7）:1164-1175.（LI Tian-qi，ZHENG Li-na.Soil E-rosion Changes in the Yanhe Watershed from 2001 to 2010 Based on RUSLE Model［J］.Journal of Natural Resources，2012，27（7）:1164-1175.（in Chinese））

［17］曹忠杰，蔡景平，何建明，等.遼寧省第四次土壤侵蝕遙感普查成果分析［J］.水土保持應(yīng)用技術(shù)，2007，（5）: 21-22.（CAO Zhong-jie，CAI Jing-ping，HE Jian-ming，et al.Investigation Result Analysis of the Fourth Remote Sensing of Soil Erosion in Liaoning Province［J］.Technology of Soil and Water Conservation，2007，（5）:21-22. （in Chinese））

［18］遼寧省水利廳.遼寧省水利基礎(chǔ)信息集［M］.北京:中國水利水電出版社，2013.（Department of Water Resources of Liaoning Province.Basic Information Collection of Liaoning Provincial Water Conservancy［M］.Beijing: China Water and Power Press，2013.（in Chinese））

（編輯:姜小蘭）

Quantitative Study on Soil Erosion in Liaoning Province Based on RUSLE Model

GAO Lu1，WANG Xue-feng2，QU Ben-liang3，CAO Yong-qiang1
（1.School of Urban and Environmental Sciences，Liaoning Normal University，Dalian 116029，China; 2.China Institute of Water Resources and Hydropower Research，Beijing 100038，China; 3.Dalian Municipal Design and Research Institute of Environmental Science，Dalian 116023，China）

Based on the rainfall，soil，remote sensing image and census data of soil erosion，in this paper，we used GIS and 3S technology in association with revised universal soil loss equation（RUSLE）model to analyze the spatial distribution characteristics of the soil erosion and relationship between the distribution and land type in Liaoning province.The results of this research are as follows:1）area of soil erosion of Liaoning province in 2011 is 4 593 600 hm2，with annual mean soil loss of 164 million tons and the average annual erosion modulus of 3 637.8 t/（km2·a）and it belongs to moderate erosion，moreover，the area of soil erosion shows the tendency of increase;2）soil erosion in hilly region of western Liaoning is the most severe and area of soil erosion is 1.919 9 million hm2，and cities are seriously affected by unreasonable human activities that lead to the increase of area of soil erosion;3）forest land is still the main land type with soil erosion in Liaoning province，accounting for about 57.66%of the total area of erosion and 38.20%of the land type.

soil erosion;RUSLE model;quantitative analysis;Liaoning province;type of land utilization

S157

A

1001-5485（2015）12-0025-05

10.11988/ckyyb.20140966

2015，32（12）:25-29

2014-11-18;

2015-02-09

國家自然科學(xué)基金項目（51579126）;流域水循環(huán)模擬與調(diào)控國家重點實驗室開放基金項目（IWHR-SKL-201501）

高璐（1991-），女，遼寧大連人，碩士研究生，從事水環(huán)境方面的研究，（電話）13352418116（電子信箱）gl1991sara@163.com。

曹永強（1972-），男，內(nèi)蒙古豐鎮(zhèn)人，教授，博士，從事水文水資源方面的研究，（電話）0411-82158364（電子信箱）caoyongqiang@lnnu.edu.cn。