基于GIS和RS的太谷縣土壤侵蝕動態變化分析

王國芳

（山西農業大學資源環境學院，山西太谷030801）

基于GIS和RS的太谷縣土壤侵蝕動態變化分析

王國芳

（山西農業大學資源環境學院，山西太谷030801）

土壤侵蝕的發生使得土壤資源受到破壞，土壤肥力下降，生態環境惡化，洪澇災害加劇，已嚴重威脅到人類的生活。以山西省太谷縣為研究區，在ENVI，ArcGIS軟件的支持下，以1990，2007，2015年的TM遙感影像為數據源，提取了影響土壤侵蝕強度的坡度因子、土地利用類型因子和植被覆蓋度因子，對太谷縣3個時期的土壤侵蝕強度進行了分級，并對近25 a間太谷縣的土壤侵蝕動態變化進行了分析。結果表明，研究區土壤侵蝕強度在25 a間主要呈現加重趨勢，微度侵蝕和輕度侵蝕呈先增加后減少的趨勢，中度侵蝕呈現逐年增加趨勢，強烈侵蝕、極強烈侵蝕和劇烈侵蝕呈現先減少后增加的趨勢。太谷縣南部區域土壤侵蝕長期比較嚴重，需采取有效措施加以防治。

土壤侵蝕；坡度；植被覆蓋度；土地利用

在自然營力與人類活動的綜合作用下，土壤或其他地面組成物質被剝蝕、破壞、分離、搬運和沉積，形成了土壤侵蝕[1-3]。土壤侵蝕的成因復雜，危害嚴重，其造成大量土壤資源被蠶食和破壞，耕地面積大量減少，土壤肥力和質量隨之下降，地表植被遭到嚴重破壞，生態環境隨之失調惡化，引發了各種地質災害，造成了巨大的經濟損失[4-5]。因此，土壤侵蝕問題已經成為環境、土壤、農業科學各界共同關注的熱點問題[6]。人們也積極地采取水利、生物、農業工程技術相結合的方式進行有效的治理。RS和GIS技術以及土壤侵蝕模型的發展，為土壤侵蝕研究提供了更加全面詳實可靠的數據，可以及時地掌握土壤侵蝕情況，為區域合理進行生態建設提供依據，從而減少不必要的經濟損失。國內外學者早在20世紀80年代就利用影像數據通過目視解譯進行土壤侵蝕的調查研究[7]，有學者通過對影響土壤侵蝕的因子加權進行土壤侵蝕確定[8]；之后有學者通過影像分類、經驗統計模型、物理過程模型、分布式模型、數字高程模型等方法對土壤進行了研究[9-13]。

筆者以山西省晉中市太谷縣為研究區，獲取研究區1990，2007，2015年3期遙感影像數據和DEM數據，進一步分析提取研究區域的坡度、植被覆蓋度和土地利用類型因子，然后對坡度、植被覆蓋度和土地利用類型3個因子基于ArcGIS進行了疊置分析，依據土壤侵蝕強度分級標準，最終確定研究區各年土壤侵蝕強度及25 a間土壤侵蝕變化情況。

1 材料和方法

1.1 研究區概況

太谷縣位于山西省晉中市，地處晉中盆地東北部，東經112°28′～113°12′，北緯37°12′～37°32′。屬溫帶季風氣候，境內四季分明，年平均氣溫小于10℃為冬季，大于22℃為夏季；冬季干冷少雪，夏季炎熱多雨，秋季秋高氣爽，春季干旱多風。地貌形態有山地、丘陵、平原3種，由東南向西北傾斜，東山主脈即為汾、漳分水嶺東段和中段。

1.2 數據來源與研究方法

1.2.1 數據來源研究收集到的數據資料主要包括遙感影像、DEM數據、土地利用現狀圖、土壤類型圖和相關的說明文檔。研究獲取的遙感數據是1990年9月16日的LANDSAT_5 TM影像，2007年8月14日的LANDSAT_5 TM影像和2015年8月4日的LANDSAT_8影像。這3期遙感影像云覆蓋率為0，影像質量較好，能滿足本研究的需要。DEM數據采用的是GDEMDEM30M分辨率數字高程數據。本研究中的遙感數據和DEM數據均來源于中國科學院計算機網絡信息中心國際科學數據鏡像網（http://www.gscloud.cn）。研究中的數據處理基于ArcGIS10.0和ENVI 5.1軟件。

1.2.2 數據預處理鑒于本研究主要提取的信息是植被指數、植被覆蓋度和土地利用類型，選擇Landsat TM/ETM＋遙感影像數據的有利于反演植被反射率的TM2波段、有利于鑒別植被種類的TM3波段和有利于調查生物量的TM4波段作為組合波段[14]。

本研究區域覆蓋了2幅DEM影像，所以，在處理影像之前首先進行了拼接，之后對DEM影像和組合之后的TM432影像根據太谷縣行政區劃邊界進行了影像的裁剪。為了突出主要信息、改善圖像的顯示質量、便于信息的提取和識別，對遙感圖像進行了輻射校正、大氣校正和增強處理。最后獲得了太谷縣1990年TM數據、2007年TM數據、2015年TM數據和DEM數據。

1.2.3 坡度因子提取坡度是對地面傾斜程度的定量描述，是影響土地利用、土壤侵蝕的主要地形因子之一[15]。本研究基于已有的DEM數據（圖1），應用ArcGIS中的坡度計算模塊提取出太谷縣的坡度因子，并根據土壤侵蝕強度分級標準中對坡度的分級標準，將研究區的坡度分為6類[16]：坡度≤5°，5°＜坡度≤8°，8°＜坡度≤15°，15°＜坡度≤25°，25°＜坡度≤35°和坡度＞35°，并生成分級坡度圖（圖2）。

1.2.4 植被覆蓋度提取采用歸一化植被指數NDVI來反映土地覆蓋的植被狀況。NDVI＝（NIRR）/（NIR＋R）。對于Landsat ETM＋數據，NIR波段對應的TM4，R波段對應的TM3波段；對于Landsat 8數據，NIR波段對應的TM5，R波段對應的TM4波段。為了減少誤差、提高精度，本研究采用NDVI＝（float（TM4）－float（TM3））/（float（TM4）＋float（TM3））計算NDVI值，2015年遙感數據采用NDVI＝（float（TM5）－float（TM4））/（float（TM5）＋float（TM4））計算NDVI值，并對計算得到的NDVI數據基于ENVI的波段運算進行二值化處理。

依據像元二分模型原理，利用fvc＝（NDVINDVIsoil）/（NDVIveg－NDVIsoil）公式計算出研究區3期的植被覆蓋度[17-19]。其中，NDVIsoil和NDVIveg分別取二值化NDVI概率累積表上5%和95%處的NDVI值[13-15]。從求得的植被覆蓋圖（圖3，4，5）上可以看到，顏色越淺的地方其植被覆蓋度值越大，顏色越深的地方植被覆蓋度值越小。

根據土壤侵蝕強度分級標準對植被覆蓋的分類要求，將植被覆蓋度分為5類：植被覆蓋度≤30%，30%＜植被覆蓋度≤45%，45%＜植被覆蓋度≤60%，60%＜植被覆蓋度≤75%，植被覆蓋度＞75%，通過決策樹分類對太谷縣植被覆蓋度進行劃分，得到1990，2007，2015年3個時期的太谷縣植被覆蓋度。

1.2.5 土地利用類型提取根據土壤侵蝕強度分級標準對土地利用類型的分類要求，將土地利用類型劃分為居民地及工礦用地、耕地、林地、草地、水域5種類型。首先，利用非監督分類方法對研究區遙感數據進行初步的分類；然后，根據5種類型的解譯標志建立訓練樣區，對非監督分類結果再進行監督分類，并且對分類結果進行聚類處理，避免出現小圖斑；最后，獲得了該研究區1990，2007，2015年3個時期的土地利用分布圖。

1.2.6 土壤侵蝕提取依據土壤侵蝕分類分級標準，參照研究區的土地利用類型、植被狀況和地形數據制訂了土壤侵蝕強度分級標準表[20]（表1）。

表1 土壤侵蝕強度分級標準

首先，將提取出的坡度數據、植被覆蓋度數據和土地利用類型數據基于ArcGIS的柵格轉面工具轉換成矢量面文件；然后，利用疊加分析功能進行疊加分析；最后，按照土壤侵蝕強度分級標準表中的組合對疊加的結果進行分類，根據分類結果形成了太谷縣1990，2007，2015年的土壤侵蝕強度數據。

2 結果與分析

2.1 土壤侵蝕面積變化分析

通過3期土壤侵蝕強度數據分析可以得到（圖6），1990年微度侵蝕占研究區總面積的35.65%，2007年微度侵蝕占到40.22%，而2015年微度侵蝕僅占22.9%，有明顯的下降趨勢。輕度侵蝕1990年占29.92%，2007年占30.82%，2015年占28.22%，變化不是特別明顯。對于中度侵蝕逐年增加，從1990年占22.33%到2015年增加到了30.67%；而對于強烈侵蝕、極強烈侵蝕和劇烈侵蝕，都表現為先減少后增加的趨勢。其中，強烈侵蝕由1990年占9.01%減少到2007年的4.25%又增加到2015年的13.13%，變化幅度較大；極強烈侵蝕和劇烈侵蝕有所增加，但是增加比例不大。

整體上看，1990，2007年從微度侵蝕到劇烈侵蝕所占面積比例逐漸減少，以微度侵蝕為主。2015年中度侵蝕所占比例最大，其次為輕度和微度侵蝕。

2.2 土壤侵蝕強度與植被覆蓋度關系分析

選取3期土壤侵蝕數據與對應的植被覆蓋度數據進行疊加分析，統計出3期不同植被覆蓋度下各類土壤侵蝕強度分級的面積，對3期中同一類別的面積求其平均值獲得不同植被覆蓋度下各類土壤侵蝕強度分級的3期平均面積（圖7）。

微度侵蝕主要發生在植被覆蓋度小于30%和植被覆蓋度大于75%的區域，可能是由于植被覆蓋度小于30%的區域主要分布在坡度小的平原城市、村莊集中的地方，所以侵蝕強度小；而植被覆蓋度的增加可以緩減土壤侵蝕的發生，所以植被覆蓋度大的區域，侵蝕強度相對較小。輕度侵蝕隨著植被覆蓋度的增加而逐漸增加，而中度侵蝕、強烈侵蝕、極強烈侵蝕和劇烈侵蝕隨著植被覆蓋度的減少，其面積在增加；強烈侵蝕和極強烈侵蝕主要發生在植被覆蓋度小于45%的區域。

2.3 土壤侵蝕強度與土地利用類型關系分析

通過土壤侵蝕強度數據與土地利用類型數據的疊加分析，研究區內居民建設用地和水域土壤侵蝕強度以微度侵蝕為主；耕地以輕度侵蝕為主；草地以中度侵蝕為主。在各種土地利用類型中，3期平均侵蝕面積由大到小依次為耕地、草地、林地、居民建設用地及水域，所占比例分別為32.18%，29.92%，21.17%,16.38%，0.35%。在各類侵蝕中，耕地類型的輕度侵蝕占的比例最大，達到了21.19%；其次為草地類型的中度侵蝕比較大，達到了14.69%；林地類型的微度侵蝕和居民地的微度侵蝕達到了14%。而水域對于侵蝕的貢獻不大。1990—2015年間，耕地類型的輕度侵蝕由24.02%降到了19.33%，中度侵蝕由9.73%降到了5.62%，草地的中度侵蝕由1990年占的比例9.06%增加到2007年的14.04%又增加到2015年的20.96%。草地類型的強烈侵蝕、極強烈侵蝕和劇烈侵蝕都呈現逐漸增加的趨勢。這可能是由于不合理的耕作方法，甚至棄耕、退耕還林不理想都會造成水土流失、土壤侵蝕加重。

3 結論與討論

本研究利用RS和GIS技術，結合地面調查資料，利用1990，2007，2015年3期遙感影像數據，提取了地面坡度、植被覆蓋度2個土壤侵蝕因子和利用非監督和監督分類相結合的方法進行了土地利用類型分類，依據土壤侵蝕強度分級標準對太谷縣3個時期的土壤侵蝕強度及近25 a間太谷縣的土壤侵蝕動態變化進行了分析，結果表明，1990—2007年，太谷縣土壤侵蝕總體呈現出逐漸減輕的趨勢，其中，以微度侵蝕面積增加最為明顯，增加了47.93 km2，增幅為4.7%；強烈侵蝕、極強烈侵蝕和劇烈侵蝕的面積均有大幅度減少，減幅分別為4.9%，1.8%和0.4%。2007—2015年太谷縣土壤侵蝕總體呈現出逐漸加重的趨勢，其中，劇烈侵蝕在這8a間增加了0.6%，極強烈侵蝕增加了37.09 km2，增幅為3.5%；強烈侵蝕增加了9%；微度侵蝕、輕度侵蝕、中度侵蝕都有所減少。

從空間變化方面看，太谷土壤侵蝕出現加劇與減弱并存的狀況，總體上土壤侵蝕呈加重趨勢，分析其原因主要是由于自然原因，如植被覆蓋度較低，坡度較高，流失面積擴散較快和人為原因破壞所造成；通過調查，山區本來退耕還林的耕地效果不佳，占用各種土地、破壞植被、開荒種地等因素造成了水土流失。

通過對太谷縣的土壤侵蝕程度分析，強烈侵蝕、極強烈侵蝕和劇烈侵蝕主要分布在太谷縣的南部，其中，大塔山植被稀疏，西北部童山禿嶺，千胡坪土層厚，地面平坦，但垣的邊沿侵蝕嚴重，植被很差，水土流失嚴重。這些侵蝕嚴重的區域多數發生在坡度大于25°以上、植被覆蓋度小于30%的區域。

［1］王茜.基于GIS和RS的冀北地區土地利用與土壤侵蝕關系研究[D].石家莊：河北師范大學，2005.

［2］張俊飚.中國土壤侵蝕影響因素及其危害分析[J].云南環境科學，2001，20（2）：4-7，29.

［3］李海燕.土壤侵蝕危害及其防治措施研究現狀[J].寧夏農林科技，2011，52（1）：71-72，77.

［4］李占斌，朱冰冰，李鵬.土壤侵蝕與水土保持研究進展[J].土壤學報，2008（5）：802-809.

［5］彭坷珊.中國土壤侵蝕影響因素及其危害分析[J].首都師范大學學報，2000，21（2）：88-94.

［6］孔澤祥.基于GIS和RS的沂河上游重點地區土壤侵蝕監測方法研究[D].濟南：山東師范大學，2012.

［7］頡耀文，陳懷錄.數字遙感影像判讀法在土壤侵蝕調查中的應用[J].蘭州大學學報：自然科學版，2002，38（2）：157-162.

［8］張增樣，趙曉麗.基于遙感和地理信息系統（GIS）的山區土壤侵蝕強度數值分析[J].農業工程學報，1998，14（3）：77-83.

［9］劉愛霞，王靜，劉正軍.三峽庫區土壤侵蝕遙感定量監測：基于和修正通用土壤流失方程的研究[J].自然災害學報，2009，18（4）：25-30.

［10］南秋菊，華洛.國內外土壤侵蝕研究進展[J].首都師范大學學報：自然科學版，2003，24（2）：86-95.

［11］鄭粉莉.土壤侵蝕預報模型研究進展[J].水土保持通報，2001，21（6）：16-18，32.

［12］蔡強國，劉紀根.關于我國土壤侵蝕模型研究進展[J].地理科學進展，2003，22（3）：142-150.

［13］張巖，劉憲春，李智廣，等.利用侵蝕模型普查黃土高原土壤侵蝕狀況[J].農業工程學報，2012，28（10）：165-171.

［14］王國芳.基于遙感技術的植被覆蓋度的動態監測研究[J].山西農業科學，2015，43（5）：592-595.

［15］蔣榮.地形因子對貴州喀斯特地區坡面土壤侵蝕的影響[D].南京：南京大學，2013.

［16］馬曉勇，柳青艷.基于RS與GIS的山西省主要煤田區土壤侵蝕研究[J].山西農業科學，2014，42（9）：987-990．

［17］徐定，彭道黎.基于像元二分模型的森林郁閉度估測方法[J].東北林業大學學報，2013（2）：119-122.

［18］馬娜，胡云鋒.基于遙感和像元二分模型的內蒙古正藍旗植被覆蓋度格局和動態變化[J].地理科學，2012（2）：251-256.

［19］郭芬芬，范建容.基于像元二分模型的昌都縣植被蓋度遙感估算[J].中國水土保持，2010（5）：65-67.

［20］費燕明.汾河流域土壤侵蝕動態變化分析[D].太原：太原理工大學，2012.

Dynamic Analysis of Soil Erosion in Taigu County Based on GIS and RS

WANGGuofang
（College ofResources and Environment，Shanxi Agricultural University，Taigu 030801，China）

The occurrence of soil erosion causes the destruction of soil resources,the decline of soil fertility,the deterioration of ecological environment and the aggravation of flood and waterlogging,which have seriously threatened the life of human beings.In this paper,taking Taigu county of Shanxi province as the study area,based on the ENVI and ArcGIS,using remote sensing image data of 1990,2007,2015 TM,three soil erosion factors type of the land slope,vegetation cover and land use were extracted.Based on the classification standard of soil erosion intensity,the soil erosion intensity during of the three periods and the dynamic changes of soil erosion in Taigu county during the past 25 years were analyzed.The results showed that soil erosion intensity in the study area during the first 25 years showed a trend of increasing.Micro erosion and slight erosion showed a trend of first increasing and then decreasing.The moderate erosion showed a trend of increasing year by year.The trend of strong erosion,extremely strong erosion and severe erosion first decreased and then increased.Based on the dynamic analysis ofsoil erosion in Taigu countyduring25 years,it shows that the soil erosion in the south of the Taigu county is more serious in the long term,and it is necessary to take effective measures to prevent and control the soil erosion.

soil erosion;slope;vegetation coverage;land use

S157

A

1002-2481（2016）12-1813-05

10.3969/j.issn.1002-2481.2016.12.19

2016-09-22

山西農業大學科技創新基金項目（201322）

王國芳（1980-），女，山西左權人，講師，主要從事3S技術應用研究工作。