摘要：以遙感（RS）和地理信息系統(tǒng)（GIS）為平臺，以四川省茂縣為研究區(qū)域，依據(jù)土壤侵蝕程度的評價原則選取對研究區(qū)域土壤侵蝕影響較大的5個評價因子，包括坡度、植被蓋度、土壤可蝕性因子、降雨侵蝕力和水土保持措施因子，運用層次分析法（AHP）計算各評價因子的權(quán)重，然后根據(jù)多因子分析將土壤侵蝕程度分成5個等級——微度侵蝕、輕度侵蝕、中度侵蝕、強度侵蝕、極強度侵蝕，得到茂縣土壤侵蝕的等級分布。結(jié)果表明，茂縣土壤侵蝕還在可控制范圍內(nèi)，沒有強度侵蝕和極強度侵蝕的地方。

關(guān)鍵詞：土壤侵蝕；時空變化；預測；RS；GIS；四川省茂縣

中圖分類號：S157.1 文獻標識碼：A 文章編號：0439-8114（2016）13-3310-06

DOI：10.14088/j.cnki.issn0439-8114.2016.13.011

土壤侵蝕，即風和水對地表母質(zhì)及其他地面組成物質(zhì)的剝蝕、搬運和沉積，是一個隨時隨地都可能發(fā)生的自然地表過程，受水力、風力、凍融、重力等外營力的作用，使土地資源遭到破壞和流失[1]。它是當今世界最嚴重的環(huán)境問題之一，嚴重威脅著人類的生產(chǎn)和發(fā)展，一直受到各國研究者的普遍關(guān)注和重視[2]。而中國是世界上水土流失最為嚴重的國家之一，近幾年來各主管部門雖然在土壤綜合治理方面做出了很大努力，但由于中國人口眾多，人均水土資源相對匱乏，土壤侵蝕依然嚴重。因而對不同條件的土壤侵蝕規(guī)律還需要進行大量深入的研究，以便為新形勢下的區(qū)域水土保持的治理和決策工作提供科學參考。

近年來，國外在土壤侵蝕機理、土壤流失方程、土壤侵蝕與生產(chǎn)力關(guān)系的評價方法、小流域水土保持規(guī)劃等方面都取得了新的進展[3]。目前，GIS與計算機技術(shù)進行統(tǒng)一規(guī)劃是國外水土保持規(guī)劃中采用的最新手段之一，該法采用數(shù)據(jù)庫與地理信息相結(jié)合的方式來描述土壤侵蝕危害的空間分布特性，目前該方法在水土保持方面已取得了很大的進展[4]。而中國對于土壤侵蝕的研究還停留在理論研究的層面上，很少對土壤的各項功能和土壤質(zhì)量的評價指標體系進行綜合質(zhì)量評價，很少把土壤質(zhì)量評價的結(jié)果在生產(chǎn)實踐中進行應(yīng)用以及推廣，而且很少對這些研究的結(jié)果進行驗證[5]。因此未來應(yīng)加強對土壤質(zhì)量評價指標體系與評價方法的研究，針對不同的土壤功能，選擇不同的土壤質(zhì)量評價指標體系，對其閾值、定量化、評價方法模型等進行深入的研究。

茂縣作為川西北地區(qū)干旱河谷的腹心地帶，經(jīng)過2008年5月12日汶川大地震以后，區(qū)域內(nèi)水土流失不斷加劇、地形地貌發(fā)生突變、生態(tài)系統(tǒng)遭到損壞，此外，其植被覆蓋率大大降低，山體出現(xiàn)大面積的裸露，動物的棲息地也遭到了不同程度的毀壞[6]。這一系列自然災害的發(fā)生，給當?shù)氐纳鷳B(tài)環(huán)境帶來了嚴重破壞，同時也不可避免地加重了土壤侵蝕的程度。因此，本研究提取茂縣地區(qū)在2000、2005、2010、2015年影響土壤侵蝕的5個評價因子，充分發(fā)揮RS的數(shù)據(jù)處理能力和GIS強大的空間分析功能，通過5個評價因子對茂縣的土壤侵蝕進行時間和空間上的對比和分析，并做出科學的評價以及趨勢變化，為區(qū)域內(nèi)水土保持工作提供可靠的科學依據(jù)。

1 研究區(qū)概況

茂縣位于四川省阿壩藏族羌族自治州東南部，與北川、安縣、綿竹、什邡、彭州、汶川、理縣、松潘縣相鄰[7]。南北寬94.8 km，東西長116.5 km，地理坐標為北緯31°24′-32°17′，東經(jīng)102°56′-104°10′，如圖1所示。茂縣作為青藏高原東麓自然地理狀況的典型代表，其地質(zhì)構(gòu)造比較復雜，區(qū)域內(nèi)斷裂十分發(fā)育，新構(gòu)造運動強烈[8]。整個地勢表現(xiàn)為由西北向東南方向逐漸降低，地表切割則由北向南加劇。茂縣是岷江上游干旱河谷區(qū)的腹心地段，海拔1 700～3 100 m，年均氣溫12 ℃，年降雨量410 mm，蒸發(fā)量1 500 mm，土壤以養(yǎng)分極低的燥褐土、石灰性褐土為主，屬高寒貧瘠干旱地帶[9]。由于區(qū)域內(nèi)干旱河谷獨特的氣候特點——干旱少雨，使區(qū)域內(nèi)植被的生長面臨極大的威脅，生長在這里的植被多是低矮、多刺的旱生型灌叢。如今，茂縣地區(qū)干旱河谷已經(jīng)呈現(xiàn)出一種次生的、脆弱的生態(tài)系統(tǒng)狀況[10]。區(qū)域內(nèi)森林資源的演變順序表現(xiàn)為常綠闊葉林→落葉闊葉林→灌叢→干旱灌叢→稀樹草坡（荒草地）→半荒漠→荒漠。

2 研究方法

2.1 數(shù)據(jù)的收集與預處理

2.1.1 數(shù)據(jù)來源 針對5個評價因子，所收集的相關(guān)數(shù)據(jù)主要包括來自中國氣象數(shù)據(jù)共享網(wǎng)的日降雨量數(shù)據(jù)（2000、2005、2010、2015年）、四川省1∶1 000 000土壤類型數(shù)據(jù)、四川省行政區(qū)劃圖、茂縣地區(qū)30 m分辨率的數(shù)字高程模型（DEM）數(shù)據(jù)以及4年的Landsat4-5 TM的遙感影像數(shù)據(jù)。

2.1.2 數(shù)據(jù)預處理 在遙感影像的成像過程中，電磁波會穿過厚度超過1 000 km的大氣層，在此過程中電磁波會受到大氣層反射、吸收、散射等影響作用，從而導致遙感影像分辨率的下降，甚至影像的失真[11]。因此，為了消除影像的失真，獲取真實的地面影像，必須對原始圖像進行圖像預處理工作，從而更充分地利用原始觀測的遙感數(shù)據(jù)。這是使用遙感影像的基礎(chǔ)性工作，它包括波段合并、幾何校正、圖像鑲嵌，圖像裁剪以及圖像增強等環(huán)節(jié)[12]。

2.2 研究過程

2.2.1 提取評價因子

1）地形因子。一般隨著坡度的增加土壤侵蝕強度會增強，但地形在土壤侵蝕中的變化還是遲緩和相對穩(wěn)定的，在描述地形的各個指標中，以坡度和坡長因子對土壤侵蝕強度的影響最大，其中又以坡度尤甚[13]。一般情況下，隨著坡度的增加，土壤侵蝕強度會加劇，因此直接采用DEM在地形分析中計算出的坡度來衡量地形因子對于土壤侵蝕強度的影響。依據(jù)土壤侵蝕遙感調(diào)查工作技術(shù)細則（表1）對茂縣坡度數(shù)據(jù)進行重分類，得到坡度因子等級分布，如圖2所示。

2）植被因子。植被是侵蝕動力的抑制因子，植被蓋度和植被結(jié)構(gòu)直接影響著土壤侵蝕強度。大量的野外實測表明，植被蓋度越高，地表的抗蝕能力就越強[14]。植被蓋度是最常用來表達植被對土壤侵蝕程度的指標[15]。遙感技術(shù)是目前植被蓋度指標獲取最為有效的手段之一，已被廣泛地利用。本研究使用TM數(shù)據(jù)提取研究區(qū)的植被指數(shù)。

f=（NDVI-NDVImin）/（NDVImax-NDVImin） （1）

式中，f為植被指數(shù)；NDVI為所求像元的植被指數(shù)，NDVImin表示研究區(qū)內(nèi)NDVI的最小值，NDVImax表示研究區(qū)內(nèi)NDVI的最大值。將得到的植被蓋度因子按照分級標準（表2）進行分級，得到4年的植被蓋度等級分布（圖3）。

3）土壤質(zhì)地因子。由于受到植物、氣候垂直自然帶的制約，茂縣的土壤類型也具有明顯的垂直地帶性[16]，從低到高依次為黃壤、黃棕壤、暗棕壤、寒棕壤、寒氈壤、寒凍氈壤和高山寒漠土，如圖4所示。

土壤質(zhì)地是土壤侵蝕的抗侵蝕因子，土壤侵蝕與成土母質(zhì)直接相關(guān)，主要表現(xiàn)在土壤的透水性、抗蝕性和抗沖性，這些特征又決定于土壤顆粒組成和結(jié)構(gòu)[17]。在通用土壤侵蝕方程中，用土壤可蝕性因子K值來定量反映土壤對侵蝕的影響，但直接測定K值要求較高，在實際應(yīng)用中，常根據(jù)成土母質(zhì)的不同來確定土壤侵蝕的影響強度[18]，采用查表法（表3）得到，其侵蝕強度分級標準見表4。

4）降雨侵蝕力因子。降雨是產(chǎn)生土壤侵蝕的基本動力和觸發(fā)因子，降雨強度越大，土壤侵蝕強度也越大，降雨強度、歷時及其在年內(nèi)的分配狀況都直接影響到土壤侵蝕的發(fā)生和侵蝕量的大小。降雨具有年際變化大、年內(nèi)分配不均和地域差異大的特點。在中國大部分地區(qū)，由于干濕交替強烈，表土疏松，汛前長期干旱，植被萌發(fā)力差，因此土壤的透水性和抗蝕能力差，很難抵御第一次暴雨對土壤的侵蝕[19]。因此，降雨侵蝕力是反映大氣降雨到地表產(chǎn)生沖擊的一種動量表示，這一因子必須采用一定時間系列的降雨量和降雨強度數(shù)據(jù)來計算。根據(jù)王萬忠等[20]對中國主要氣象站點的年降雨侵蝕力計算公式進行計算，具體為：

R=0.207（P×I60/100）×1.205 （2）

式中，R表示降雨侵蝕力；P表示年降雨量；I60表示年內(nèi)最大60 min的降雨量，通過查閱資料，得到西南地區(qū)的I60為11.2 mm，由于資料不足，文中直接以西南地區(qū)的研究資料代替茂縣。通過計算，得出2000、2005、2010、2015年的降雨量分別為4 268、5 091、4 506、5 018 mm，代入式中計算得到降雨侵蝕力分別為350.5、433.5、374.2、426.0。降雨侵蝕力因子強度分級標準見表5。

5）人為因素。土壤侵蝕時常發(fā)生，這是地球表面不斷變化的一種自然規(guī)律，通常稱為正常的自然地質(zhì)侵蝕。然而隨著人類活動的逐漸頻繁，農(nóng)作物輪作周期短、森林砍伐和超負荷放牧等使土壤侵蝕和土壤生成之間的自然平衡遭到破壞，又加速了土壤侵蝕的發(fā)生[21]。人為因素雖然對土壤侵蝕產(chǎn)生了十分重要的影響，但是在進行定量的計算中一般不作為計算因子，因此本研究中將人為因素賦為常數(shù)1。

2.2.2 土壤侵蝕綜合評價

1）確定各評價因子權(quán)重。按照各指標的重要程度，在歸納前人的研究成果[22]和綜合專家知識后得出如下準則：地形（B1）>降雨侵蝕力（B4）>植被覆蓋度（B2）>土壤可蝕性（B3）>水土保持措施（B5），其矩陣及權(quán)重值（W）見表6。

2）土壤侵蝕影響因子分級及監(jiān)測。前面通過層次分析法確定了影響土壤侵蝕的各因子權(quán)重，結(jié)合前人研究成果[23]和關(guān)于土壤侵蝕類型及強度分級標準的規(guī)定把土壤侵蝕強度分為5級進行標準化處理，侵蝕強度為微度的等級為一，賦值1；輕度等級為二，賦值3；中度等級為三，賦值5；強度等級為四，賦值7；極強度等級為五，賦值9。

采用加權(quán)和法作為評價土壤侵蝕程度的模型，評價模型的表達式為：

式中，S為土壤侵蝕程度綜合評價的總分值；n為評價因子個數(shù)；Wk為第k個評價因子的權(quán)重；Vk為第k個評價因子的量化值。加權(quán)指數(shù)和法即將各個評價因子的權(quán)重和量化值的乘積求和，得到綜合評價分值。

根據(jù)層次分析法得到的各因子權(quán)重以及上述評價模型[24]對坡度、植被蓋度、降雨侵蝕力、土壤可蝕性等4個主要評價因子的等級圖進行疊加，然后通過疊加的綜合評價分值對土壤侵蝕程度進行分區(qū)，按綜合分的高低劃分：分值越高，表示侵蝕程度越強；分值越低，表示越弱。侵蝕等級劃分為微度侵蝕、輕度侵蝕、中度侵蝕、強度侵蝕、極強度侵蝕5個等級，如表7所示。

3 結(jié)果與分析

根據(jù)坡度、植被蓋度、降雨侵蝕力、土壤可蝕性因子的綜合評價分，將綜合評價分按微度侵蝕、輕度侵蝕、中度侵蝕、強度侵蝕、極強度侵蝕5個等級來劃分[25]，得到土壤侵蝕等級分布圖（圖5）。

通過統(tǒng)計對比得到各等級土壤侵蝕面積的單元數(shù)，具體見表8。由表8可知，2000年微度侵蝕面積占總面積的0.02%，輕度侵蝕面積占39.21%，中度侵蝕面積占60.77%；2005年微度侵蝕面積占總面積的0.03%，輕度侵蝕面積占45.37%，中度侵蝕面積占54.60%。由此可知，與2000年相比，2005年微度侵蝕和輕度侵蝕面積增加，中度侵蝕面積有所減少，說明相關(guān)部門采取了一定的有效措施對其進行了強有力的控制。2010-2015年，土壤侵蝕出現(xiàn)了逆轉(zhuǎn)情況，輕度侵蝕有所減少，中度侵蝕有所增加，主要是由于2008年汶川地震后，該區(qū)域水土急劇流失，植被遭到了很大的破壞，因此土壤侵蝕變得更加嚴重。從總體上來說，茂縣地區(qū)土壤侵蝕還在可控制范圍內(nèi)，沒有強度侵蝕和極強度侵蝕的地方，基本上是輕度侵蝕和中度侵蝕，微度侵蝕面積雖然有所變化，但變化很小，基本上可以忽略。

4 結(jié)論

本研究結(jié)果表明，在影響土壤侵蝕程度的眾多因子中，有變化的是植被蓋度因子和降雨侵蝕力因子，但是在所選取的4年中，降雨量相近，所以降雨侵蝕力的影響作用并不大，而植被卻是主要因子，因此可以預測當研究區(qū)域內(nèi)植被狀況有所好轉(zhuǎn)的情況下，土壤侵蝕程度會有不同程度的降低，因此相關(guān)部門應(yīng)該采取強有力的措施，做好區(qū)域內(nèi)植被保護和生態(tài)建設(shè)工作。

近年來，研究人員對岷江上游土壤侵蝕的問題做了大量研究，其研究都是大區(qū)域范圍內(nèi)，對周邊具體縣鎮(zhèn)研究并不多見[26]。因此，此次選擇茂縣為研究區(qū)域，縮小了研究范圍，對其做出了更加具體的分析。結(jié)合此次研究結(jié)果，在此提出以下幾點具有針對性的意見：①對坡耕地的開發(fā)和工程建設(shè)進行改革，完成緩坡耕地保土耕作措施的推廣與應(yīng)用；②茂縣水土保持工作的重點是封山育林育草，退耕還林，退牧還草，合理調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，實施宜林則林、宜草則草、宜農(nóng)則農(nóng)，禁止天然林砍伐和亂墾土地；③加速水土保持林（草）建設(shè)，興建小型水利水保工程，保護林草植被，涵養(yǎng)水源，重點防止對生態(tài)環(huán)境造成新的破壞；④應(yīng)加大森林及灌叢的保護力度、加快對道路建設(shè)和工礦開發(fā)引發(fā)的水土流失區(qū)的治理[27]。

參考文獻：

[1] 李國瑞，王貴平，封九梁，等.土壤侵蝕模型研究的現(xiàn)狀與發(fā)展趨勢[J].太原理工大學學報，2003，34（1）：100-103.

[2] 伏耀龍.岷江上游干旱河谷區(qū)土壤質(zhì)量評價及侵蝕特征研究[D].陜西楊凌：西北農(nóng)林科技大學，2012.

[3] 楊勤科，李 銳，曹明明.區(qū)域土壤侵蝕定量研究的國內(nèi)外進展[J].地球科學進展，2006，21（8）：850-853.

[4] 楊艷生，史德明.長江三峽區(qū)土壤侵蝕研究[M].南京：東南大學出版社，1994.

[5] 趙曉麗，張增祥，周全斌，等.中國土壤侵蝕現(xiàn)狀及綜合防治對策研究[J].水土保持學報，2002，16（1）：40-43.

[6] 趙英時.遙感應(yīng)用分析原理與方法[M].北京：科學出版社，2003.

[7] 鄭 睿.基于RS和GIS的漓江上游土壤侵蝕研究[D].昆明：昆明理工大學，2011.

[8] 劉武林.岷江流域土壤侵蝕演變與治理對策研究[D].四川雅安：四川農(nóng)業(yè)大學，2008.

[9] 章文波，付金生.不同類型雨量資料估算降雨侵蝕力[J].資源科學，2003，25（1）：36-41.

[10] 馮軍輝.基于GIS平臺汶川震后新生水土流失分布特征研究[D].重慶：西南大學，2012.

[11] 鄧良基，侯大斌，王昌全，等.四川自然土壤和旱耕地土壤可蝕性特征研究[J].中國水土保持，2003（7）：23-25.

[12] 房世波，楊武年，張新時.基于GIS的農(nóng)業(yè)土壤侵蝕分級和退耕還林決策——以川西干旱河谷地區(qū)為例[J].青島大學學報（工程技術(shù)版），2005，20（1）：72-77.

[13] 郭麗霞.基于GIS的土壤侵蝕潛在危險度評價——以共和盆地塔拉灘為例[D].西寧：青海師范大學，2013.

[14] 王 茜，王 衛(wèi)，呂昌河.基于GIS和RS的土地利用與土壤侵蝕關(guān)系研究——以冀北地區(qū)為例[J].中國水土保持科學，2006， 4（6）：37-41.

[15] 張愛國，張平倉，楊勤科.區(qū)域水土流失土壤因子研究[M].北京：地質(zhì)出版社，2003.

[16] 羅細芳，姚小華.水土流失機理與模型研究進展[J].江西農(nóng)業(yè)大學學報：自然科學，2004，26（5）：813-817.

[17] KNOESEN D M，SMITHERS J C. The development and assessment of a regionalised daily rainfall disaggregation model for South Africa[J].2009，54（3）：217-233.

[18] 徐澤水.關(guān)于層次分析中幾種標度的模擬評估[J].系統(tǒng)工程理論與實踐，2000，20（7）：58-62.

[19] 張建平，樊 宏，葉延瓊.岷江上游土壤侵蝕及其防治對策[J].水土保持學報，2002，16（S1）：19-22.

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

[21] 蔡強國，劉紀根.關(guān)于我國土壤侵蝕模型研究進展[J].地理科學進展，2003，22（3）：242-250.

[22] 孫 凡，游 翔，劉伯云，等.四川德陽地區(qū)土壤侵蝕因子與敏感性分析[J].西南師范大學學報（自然科學版），2007，32（4）：47-52.

[23] JASROTIA A S， DHIMAN S D， AGGARWAL S P. Rainfall-runoff and soil erosion modeling using remote sensing and GIS technique-a case study of tons watershed[J]. Journal of the Indian Society of Remote Sensing，2002，30（3）：167-180.

[24] KINIRY J R，WILLIAMS J R，SRINIVASAN R. Soil and water assessment tool user’s manual[J]. Nature Clinical Practice Rheumatology，2007，3（3）：119.

[25] 郭廷輔.我國水土保持工作現(xiàn)狀、問題和對策[J].地理研究，1995（4）：1-7.

[26] 周忠軒，吳 鋼，邵國凡.遙感探測土地植被覆蓋指數(shù)的準確度評估[J].應(yīng)用生態(tài)學報，2004，15（1）：36-38.

[27] 黃光忠，劉向東，何 飛.岷江上游地震災后受損植被狀況及其恢復重建對策[J].四川林業(yè)科技，2009，30（3）：95-99.