劉文輝， 璩向?qū)帲?李曉宇， 曹園園

(寧夏大學西北土地退化與生態(tài)恢復省部共建國家重點實驗室培育基地/西北退化生態(tài)系統(tǒng)恢復與重建教育部重點實驗室，寧夏銀川 750021)

土壤侵蝕是土壤或其地面組成物質(zhì)在自然營力作用下或在自然營力與人為活動的綜合作用下被剝蝕、破壞、分離、搬運和沉積的過程[1]。土壤侵蝕包括水蝕、風蝕、重力侵蝕、人為侵蝕等不同類型。其影響是破壞土壤結構、減少土壤養(yǎng)分，不僅使土地生產(chǎn)力下降、水土資源流失，而且造成河湖庫塘的泥沙淤積，加劇洪澇、干旱等災害發(fā)生，已成為危及人類生存與發(fā)展的重要環(huán)境問題之一[2]。土壤侵蝕在自然狀態(tài)下也一直存在，但劇烈的人類活動加速了土壤侵蝕的發(fā)生[3]。因此，關注大型生態(tài)工程驅(qū)動下的生態(tài)環(huán)境要素與土壤侵蝕過程的變化[2]很有必要。

修正通用土壤流失方程(the revised universal soil loss equation，RUSLE)是在美國農(nóng)業(yè)部開發(fā)的通用土壤流失方程(the universal soil loss equation，USLE)的基礎上于1992年改進完善而來，是目前世界上應用最廣泛的基于坡面的水蝕預報經(jīng)驗模型。隨著地理信息系統(tǒng)(GIS)技術不斷盛行，基于GIS運用USLE/RUSLE定量評估土壤侵蝕，已成為合理規(guī)劃、利用區(qū)域水土資源的有效手段，在世界許多國家得到廣泛研究和應用。我國學者已將USLE/RUSLE模型應用于小流域、流域、區(qū)域尺度的土壤侵蝕預測預報。秦偉等使用RUSLE模型并改進坡長算法，評估了黃土高原四面窯溝小流域的土壤侵蝕狀況[4]；龐國偉通過對人為因素的定量表征在RUSLE模型基礎上研究了黃土高原延河流域的土壤侵蝕與人類活動之間的關系[5]；郭建英等研究了退耕還林前后陜西吳起縣土壤侵蝕變化[6]；白曉松對北方山區(qū)土壤侵蝕進行了評價[7]。在GIS支持下，通過建立土壤侵蝕因子數(shù)據(jù)庫，用土壤侵蝕模型進行土壤侵蝕定量評價，是區(qū)域和大中流域土壤侵蝕研究的基本趨勢[8]。

本試驗采用RUSLE模型，在GIS與RS支持下對寧夏彭陽縣2000、2008、2014年3年的土壤侵蝕進行研究，考察退耕還林(草)工程前中后3期的土壤侵蝕變化情況，分析不同流域、不同地質(zhì)區(qū)域的土壤侵蝕變化狀況，引入地形位指數(shù)與灰度關聯(lián)分析，研究坡度與人類活動對彭陽縣土壤侵蝕的影響，以期為區(qū)域生態(tài)環(huán)境建設、退耕還林(草)工程的下一步開展提供科學參考。

1 材料與方法

1.1 研究區(qū)概況

彭陽縣位于寧夏回族自治區(qū)南部，屬固原市，縣域位于106°32′～107°47′E、35°41′～36°17′N之間，地處黃土丘陵溝壑區(qū)第二副區(qū)，海拔1 248～2 418 m。地形分為北部黃土丘陵區(qū)、中部河谷殘源區(qū)、西南部土石質(zhì)山區(qū)3個自然類型區(qū)，地勢由西北向東南呈波狀傾斜，境內(nèi)地形支離破碎，山高坡陡，溝壑縱橫，是寧夏水土流失最嚴重的縣之一，也是全國重點水土保持區(qū)。彭陽縣境內(nèi)河流水系屬于涇河流域，主要有茹河、紅河、安家川河。紅河是涇河一級支流，茹河、安家川河是涇河二級支流。土壤類型主要包括黃綿土、黑壚土、山地灰褐土等。多年平均降水量為433.6 mm，降水主要集中于6—10月，呈南多北少、西豐東欠的趨勢。多年平均水面蒸發(fā)量為960 mm，全年無霜期為120～160 d，全縣平均溫度為 7.0 ℃。

彭陽縣自20紀70年代開始水土保持生態(tài)恢復工程，經(jīng)過近40年的不懈努力，縣域內(nèi)生態(tài)治理卓有成效，2003年被國家環(huán)保部列為第8批全國生態(tài)示范區(qū)，“彭陽模式”成為全國生態(tài)治理的典范。自2000年退耕還林(草)工程開始后，彭陽縣大力開展退耕還林還草、坡改梯、小流域治理等生態(tài)工程。在實施系列生態(tài)工程，特別是對生態(tài)環(huán)境影響巨大的退耕還林(草)工程后，彭陽縣土壤侵蝕狀況對退耕還林(草)工程生態(tài)效益的評價具有重要意義。

1.2 研究方法

1.2.1 數(shù)據(jù)來源 本研究所用數(shù)據(jù)來源：(1)彭陽縣2000、2008、2014年3年境內(nèi)月降水量數(shù)據(jù)，由彭陽縣水務局提供；(2)參考彭陽縣1 ∶30萬土壤類型圖[9]并結合實地考察，在ArcGIS 10.2軟件中生成30 m分辨率的矢量土壤類型圖，根據(jù)《寧夏土壤》[10]與中國土壤數(shù)據(jù)庫(http：//vdb3.soil.csdb.cn/)為土壤類型圖添加土壤質(zhì)地屬性數(shù)據(jù)；(3)彭陽縣30 m分辨率的ASTER GDEM數(shù)據(jù)來源于中國科學院計算機網(wǎng)絡信息中心地理空間數(shù)據(jù)云平臺(http：//www.gscloud.cn)；(4)基于Landsat/TM5(2000年8月30日)、Landsat/TM7(2008年9月5日、2014年9月11日)三景影像，來源于中國科學院計算機網(wǎng)絡信息中心地理空間數(shù)據(jù)云平臺(http：//www.gscloud.cn)。在ENVI5.0軟件中經(jīng)幾何校正、研究區(qū)裁切、輻射定標、大氣校正后，應用監(jiān)督分類且經(jīng)實地驗證后解譯獲得3期土地利用，生成30 m分辨率的矢量土地利用圖，地類劃分為耕地、林草地、未利用土地、水域灘涂、建設用地5類，同時計算生成植被覆蓋度柵格圖。以上所有圖采用UTM投影，地理坐標系均為WGS-1984。

1.2.2 數(shù)據(jù)處理 本研究采用的修正通用土壤流失方程(RUSLE)如下[11]：

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

(1)

式中：A為單位面積上的土壤流失量，t/(hm2·年)；R為降雨侵蝕力因子，(MJ·mm)/(hm2·h·年)；K為土壤可蝕性因子，[t/(hm2·年)]·[hm2/(MJ·mm)]；L為坡長因子；S為坡度因子；C為植被覆蓋與管理因子；P為水土保持措施因子。L、S、C、P均為無量綱因子。不同因子計算公式方法見表1、表2。

表1 RUSLE各因子計算方法

表2 彭陽縣水土保持措施因子(P)賦值

根據(jù)公式(1)計算2000、2008、2014年代表退耕還林工程的前中后3期年均土壤侵蝕量/土壤侵蝕模數(shù)，按照我國水利部頒布的SL190—2007《土壤侵蝕分級分類標準》[17]進行土壤侵蝕強度分級。

2 結果與分析

經(jīng)過計算，彭陽縣土地利用以旱地與林草地為主，少量的建設用地、未利用土地、水域灘涂。因此，彭陽縣土壤侵蝕主要來源于耕地與林草地。經(jīng)過退耕還林建設從2000—2014年耕地面積減少，林草地面積增加，耕地減少的幅度略高于林草地增加的幅度。

彭陽縣25°以下區(qū)域占全縣總面積超過80%，且五大地類均集中分布于25°以下區(qū)域。經(jīng)過退耕，至2014年旱地主要分布在25°以下地形區(qū)域，其中8～15°最多。林草地面積則逐年大幅增加，主要分布于8°～35°之間。林草地在15°以上區(qū)域分布面積遠遠大于旱地，且15°～25°區(qū)域分布最多，向兩邊依次遞減。2008—2014年是退耕還林還草的加速階段，其耕地減少面積與林草地增加面積變化最多。表明彭陽縣退耕還林(草)工程對旱地、林草地變化影響很大，土壤侵蝕年季變化的主要影響因子來自于地類的變化及與之相關的退耕還林等生態(tài)工程的進行。

2.1 土壤侵蝕模數(shù)的時空變化

從全縣侵蝕模數(shù)來看，2000年彭陽縣為65.47 t/(hm2·年)，屬于強烈侵蝕；2008年為119.04 t/(hm2·年)，屬于極強烈侵蝕；2014年為84.99 t/(hm2·年)，屬于極強烈侵蝕，但低于2008年的侵蝕模數(shù)，土壤侵蝕得到一定的控制。從流域來看，3個流域3期的土壤侵蝕模數(shù)均呈先上升后下降的趨勢，整體較2000年都有所惡化。其中，安家川流域在3期均處于較高的侵蝕強度與侵蝕模數(shù)，土壤侵蝕狀況最為嚴重(表3)。

彭陽縣12個鄉(xiāng)鎮(zhèn)分布在2個地貌地形單元中，西南部的古城鄉(xiāng)、新集鄉(xiāng)大部位于六盤山土石質(zhì)山區(qū)，其他鄉(xiāng)鎮(zhèn)位于黃土丘陵溝壑區(qū)。經(jīng)計算發(fā)現(xiàn)，土石質(zhì)山區(qū)鄉(xiāng)鎮(zhèn)整體土壤侵蝕模數(shù)均低于黃土丘陵區(qū)，2000、2008、2014年土石質(zhì)山區(qū)土壤侵蝕模數(shù)分別為47.11、76.87、55.63 t/(hm2·年)，而同期黃土丘陵區(qū)為71.06、133.11、94.95 t/(hm2·年)，這與六盤山土石質(zhì)山區(qū)植被覆蓋度高于黃土丘陵溝壑區(qū)，同時人為開發(fā)活動也較小有關。

表3 分流域侵蝕模數(shù)及侵蝕強度

2.2 土壤侵蝕強度轉移矩陣分析

2.2.1 不同年份不同侵蝕強度比例 2000、2008、2014年3年的6個侵蝕強度級別占縣域總面積比例見圖1，微度侵蝕占比最多，面積均超過縣域面積的35%。其中，2008年微度侵蝕占比最多，2014年最少。侵蝕程度最高的劇烈侵蝕以2008年面積占比最大，具體表現(xiàn)為2008年>2014年>2000年。從2000—2014年強烈及以上侵蝕強度占比之和逐年擴大，2000—2008年階段擴大最多。表明從2000—2014年研究區(qū)彭陽縣土壤侵蝕強度強烈以上等級貢獻逐年增加，2000—2008年是一個急速變化的階段。2008年微度侵蝕及劇烈侵蝕面積占比均是3年中最大的，2008年水土流失狀況表現(xiàn)了向2個侵蝕強度極端發(fā)展。水土流失狀況相對于2000、2008、2014年均有所惡化。

通過ArcGIS字段計算器功能，對2000—2008、2008—2014年2個階段土壤侵蝕強度等級轉移情況進行計算。2000—2008、2008—2014年沒有發(fā)生轉移的面積占縣域面積基本在60%左右，2000—2008年侵蝕強度轉移主導方向為向上轉移，向上轉移占縣域面積的32.09%，向上轉移1級最多；而2008—2014年轉移為主導方向為向下轉移，向下轉移1級最多。表明2008年侵蝕情況狀況惡化嚴重，2014年有了較大的改觀。這與退耕還林深入推進，地類在2008—2014年變化較大有關。

2.2.2 2000—2008年土壤侵蝕強度轉移矩陣 2000—2008年土壤侵蝕強度轉移矩陣見表4。

表4 2000—2008年土壤侵蝕強度轉移矩陣

2000年有40.44%面積在2008年發(fā)生了轉移。從2000年轉出情況來看，轉出面積大小依次為輕度侵蝕>極強烈侵蝕>中度侵蝕>微度侵蝕>強烈侵蝕>劇烈侵蝕。輕度侵蝕主要轉向中度侵蝕及以上級別。從微度侵蝕到極強烈侵蝕均表現(xiàn)為向上轉移的趨勢，向上轉移總面積遠遠超過向下轉移總面積。

從轉入情況來看，2000年各侵蝕強度向2008年劇烈侵蝕轉移占研究區(qū)總面積最多，其次為向2008年中度侵蝕轉。從轉入構成來看，除微度侵蝕，輕度及以上侵蝕均是以低等級侵蝕強度向高等級侵蝕強度轉移的趨勢。這與退耕還林還草等生態(tài)工程后，農(nóng)地轉為林草地，初期林草地植被覆蓋度較低、植被不穩(wěn)定有關。

2.2.3 2008—2014年土壤侵蝕強度轉移矩陣 2008—2014年土壤侵蝕強度轉移矩陣見表5。

表5 2008—2014年土壤侵蝕強度轉移矩陣

2008年有39.46%的面積在2014年發(fā)生轉移，但是轉移幅度相對于2000年有所變小(圖2)。從2008年轉出來看，劇烈侵蝕轉出最多，在2014年轉向劇烈侵蝕以下強度，轉出面積大小情況為劇烈侵蝕>中度侵蝕>微度侵蝕>輕度侵蝕>極強烈侵蝕>強烈侵蝕。2008—2014年階段總體呈各侵蝕強度向下轉移趨勢，向下轉移總面積大于向上轉移總面積，其中劇烈侵蝕和極強烈侵蝕向下轉移最多。表明2008年土壤侵蝕情況在2014年得到了一定控制，尤其是劇烈侵蝕和中度侵蝕情況有所減少。

2.2.4 兩階段轉移變化 從轉入看2014年極強烈侵蝕接受了最多的轉入，其中由2008年的劇烈侵蝕向下轉移而來，具體轉入面積大小情況為極強烈侵蝕>輕度侵蝕>中度侵蝕>劇烈侵蝕>強烈侵蝕>微度侵蝕。

整體來看，2000—2008年土壤侵蝕狀況急速惡化，低級別侵蝕強度向高級別侵蝕強度轉移，2008年出現(xiàn)了極強侵蝕與弱侵蝕面積擴大的兩極分化現(xiàn)象，而2014年土壤侵蝕狀況相比2008年有所控制，特別是在極強烈侵蝕上。這與退耕還林、坡改梯等生態(tài)工程逐步完善穩(wěn)定之后，有效地控制了土壤侵蝕狀況。

2.3 土壤侵蝕的地形位分布指數(shù)分析

對不同侵蝕強度下的時空分布進行統(tǒng)計計算，結合地類、地形位分布指數(shù)分析侵蝕強度的變化。地形位分布指數(shù)采用如下公式[18]：

P=(Sie/Si)/(Se/S)。

(2)

式中：P是地形位分布指數(shù)；Sie為e地形位下第i級土壤侵蝕的面積；Si為整個研究區(qū)內(nèi)第i級土壤侵蝕的總面積；Se表示整個工作區(qū)內(nèi)e種地形位的總面積；S為整個研究區(qū)域的面積。P值越大，說明某土壤侵蝕類型出現(xiàn)頻率越高。參照水利部《土壤侵蝕分類分級標準》將坡度分為0～8、8～15、15～25、25～35、>35[17]。

從圖3可以看到，由于耕地、林草地集中分布于25°以下區(qū)域，所以6個侵蝕強度均在25°以下區(qū)域有較高的侵蝕地形位指數(shù)。而劇烈侵蝕強度下，35°以上區(qū)間的土壤侵蝕地形位指數(shù)3年均達到極高值，同時從中度侵蝕開始至2014年在35°以上區(qū)域地形位指數(shù)逐漸升高，說明經(jīng)過退耕還林后，中度以上土壤侵蝕強度在較高坡度區(qū)域主要是坡度提高導致了侵蝕的發(fā)生。集中分布于較低坡度的耕地主要產(chǎn)生了微度、輕度、中度侵蝕，林草地對強烈、極強烈、劇烈侵蝕產(chǎn)生了主要貢獻，耕地與林草地在微度、中度侵蝕產(chǎn)生了交互作用。

在今后水土保持工作與退耕還林工程中，8°以下區(qū)域要杜絕輕度與中度侵蝕的發(fā)生，關鍵是做好耕地的水土流失防治。在8°～15°區(qū)域做好坡改梯工作，在田埂上合理安排喬灌草的護坡工程，對于冬小麥地塊要做到桔梗還田減少暴雨的直接沖刷擊濺侵蝕與地表徑流的形成。中度侵蝕、強烈侵蝕、極強烈侵蝕、35°以下的劇烈侵蝕主要由于退耕還林后開始后，植被覆蓋度降低、植被存活率低、工程后管護不到位、植被不穩(wěn)定造成的，所以要相應地合理配置喬灌草，做好之后的管護工作。

2.4 土壤侵蝕灰度關聯(lián)分析

為了更好地了解彭陽縣土壤侵蝕主要影響因子對土壤侵蝕的作用，使用灰度關聯(lián)分析對人為作用顯著的水土保持措施因子(P)、植被覆蓋與管理因子(C)進行計算分析?；叶汝P聯(lián)分析公式如下：

(2)

式中：ξ(K)為關聯(lián)系數(shù)，ρ為分辨系數(shù)。分辨系數(shù)不宜取值太大，取值范圍為0.1～0.5，本研究選取0.1。

在ArcGIS 10.1軟件中，以象元為單位分別計算P、C因子相對于土壤侵蝕模數(shù)的關聯(lián)系數(shù)，最后使用P、C因子關聯(lián)系數(shù)的平均值作為關聯(lián)度。關聯(lián)度取值范圍為0～1，越接近于1，因子對土壤侵蝕的作用影響越大，分析結果見表6。

表6 彭陽縣P、C因子與土壤侵蝕模數(shù)關聯(lián)度

從表6可以看出，3年中水土保持措施因子均對土壤侵蝕的影響最大，說明地類的變化是引起彭陽縣土壤侵蝕變化的主要原因。因此，在以后生態(tài)建設與退耕還林工作中，要注意合理地規(guī)劃調(diào)整地類。

3 結論

運用GIS和RS技術，基于RUSLE模型定量計算寧夏退耕還林還草、生態(tài)恢復的典范縣——彭陽縣的土壤侵蝕模數(shù)，在退耕還林(草)工程的背景下考察彭陽縣2000、2008、2014年3期土壤侵蝕強度的時空變化，同時引入地形位指數(shù)與灰度關聯(lián)分析，研究彭陽縣坡度變化、植被覆蓋度、地類變化對土壤侵蝕的影響。結果顯示：(1)2000年彭陽縣土壤侵蝕模數(shù)為65.47 t/(hm2·年)，屬強烈侵蝕；2008年為 119.041 t/(hm2·年)，屬于極強烈侵蝕；2014年為 84.99 t/(hm2·年)，屬于極強烈侵蝕，但低于2008年的侵蝕模數(shù)，土壤侵蝕得到一定的控制。(2)從空間分布上看，彭陽縣的3個流域在3期土壤侵蝕模數(shù)均是呈先上升后下降的趨勢，整體較2000年都有所惡化。其中，安家川流域在3期均處于較高的侵蝕強度與侵蝕模數(shù)，土壤侵蝕狀況最為嚴重。土石質(zhì)山區(qū)的土壤侵蝕模數(shù)均低于黃土丘陵溝壑區(qū)。(3)從侵蝕強度轉移情況來看，2000—2008年間向上轉移為主要趨勢，2008—2014年為向下轉移趨勢，2008年為侵蝕情況最為嚴重的1年。(4)從土壤侵蝕強度的坡度分布來看，35°以上區(qū)域在中度以上侵蝕強度主要是坡度提高導致侵蝕的發(fā)生。25°以下區(qū)域主要是耕地與林草第交互作用產(chǎn)生了土壤侵蝕。(5)通過灰度關聯(lián)分析，彭陽縣3期土壤侵蝕人為影響中，地類的變化對土壤侵蝕的發(fā)生作用最大。

從2000、2008、2014年3年的土壤侵蝕及其強度變化可以看出，在退耕還林(草)工程開展的背景下，彭陽縣土壤侵蝕受地類的變化出現(xiàn)了變化，2008年土壤侵蝕嚴重至2014年土壤侵蝕受到抑制，反映了退耕還林還草后生態(tài)恢復治理的成效。本試驗僅研究了間隔3年的土壤侵蝕及其變化，今后還須對連續(xù)時間尺度上土壤侵蝕的變化進行完整的研究。

：

[1]吳發(fā)啟，張洪江. 土壤侵蝕學[M]. 北京：科學出版社，2012：241.

[2]冷疏影，馮仁國，李 銳，等. 土壤侵蝕與水土保持科學重點研究領域與問題[J]. 水土保持學報，2004，18(1)：1-6.

[3]怡 凱，王詩陽，王 雪，等. 基于RUSLE模型的土壤侵蝕時空分異特征分析——以遼寧省朝陽市為例[J]. 地理科學，2015，35(3)：365-372.

[4]秦 偉，朱清科，張 巖. 基于GIS和RUSLE的黃土高原小流域土壤侵蝕評估[J]. 農(nóng)業(yè)工程學報，2009，25(8)：157-163.

[5]龐國偉. 人為作用對土壤侵蝕環(huán)境影響的定量表征[D]. 北京：中國科學院研究生院，2012.

[6]郭建英，李錦榮，劉鐵軍，等. 退耕還林工程建設對吳起縣土地利用/覆被變化及其土壤侵蝕的影響[J]. 水土保持研究，2013，20(5)：6-11.

[7]白曉松. 基于RUSLE的北方山區(qū)土壤侵蝕定量研究與生態(tài)適宜性評價[D]. 保定：河北農(nóng)業(yè)大學，2010.

[8]謝紅霞. 延河流域土壤侵蝕時空變化及水土保持環(huán)境效應評價研究[D]. 西安：陜西師范大學，2008.

[9]中國科學院寧夏回族自治區(qū)固原縣考察隊.黃土高原典型地區(qū)寧夏固原縣綜合農(nóng)業(yè)區(qū)劃與應用[M]. 銀川：寧夏人民出版社，1988：331.

[10]王吉智. 寧夏農(nóng)業(yè)勘查設計院.寧夏土壤[M]. 銀川：寧夏人民出版社，1990：544.

[11]劉寶元，謝 云，張科利. 土壤侵蝕預報模型[M]. 北京：中國科學技術出版社，2001：251.

[12]章文波，付金生. 不同類型雨量資料估算降雨侵蝕力[J]. 資源科學，2003，25(1)：37-43.

[13]Wang B，Zheng F，Guan Y. Improved USLE-K，factor prediction：a case study on water erosion areas in China[J]. International Soil & Water Conservation Research，2016，4(3)：168-176.

[14]崔 晨. 基于DEM的土壤侵蝕模型中地形因子的研究[D]. 西安：西北大學，2010.

[15]潘建平. RUSLE及其影響因子的快速計算分析[J]. 地質(zhì)災害與環(huán)境保護，2008，19(1)：88-92.

[16]蔡崇法，丁樹文. 應用USLE模型與地理信息系統(tǒng)IDRISI預測小流域土壤侵蝕量的研究[J]. 水土保持學報，2000，14(2)：19-24.

[17]中華人民共和國水利部. 土壤侵蝕分類分級標準：SL190—2007[S]. 北京：中國標準出版社，2008.

[18]高鳳杰，張 柏，王宗明，等. 基于GIS與USLE的牡丹江市退耕還林前后水土流失變化研究[J]. 農(nóng)業(yè)現(xiàn)代化研究，2010，31(5)：612-616.