基于GIS和RUSLE的玉溪市土壤侵蝕研究

孫詠琦，李建華，劉士鑫，付景

（云南農業大學，a.水利學院；b.資源與環境學院，昆明 650201）

土壤作為關乎人類生存的物質基礎與自然資源，具有不可再生的屬性，土壤侵蝕不僅因有機質流失而引起土地退化、土壤肥力下降，還會因泥沙攜帶氮、磷、鉀等營養物質和其他污染物進入河道系統而造成水環境惡化［1］。土壤侵蝕直接影響區域生態環境和社會經濟的可持續發展，成為各國普遍關注的重要環境問題之一。近些年來，結合地理信息系統和遙感進行空間分析的手段日益成熟，結合RS/GIS技術和修正通用土壤流失方程（Universal soil loss equation，RUSLE）的區域土壤侵蝕定量評估，顯著提高了監測及預測效率和空間可視化能力［2，3］，在世界中小流域或區域的土壤流失研究中得到了廣泛應用［4］。

通過閱讀文獻了解到，國內外學者自19世紀就開始對土壤侵蝕進行研究［5-9］。隨著研究的深入，基于經驗的土壤侵蝕模型RUSLE逐漸被眾人認可。該模型因操作簡單計算方便，加之RS和GIS的迅猛發展，已被廣泛地應用在一些區域的土壤侵蝕的相關研究，成為主流的土壤侵蝕研究方法［10］。Roberto等［11］利用RUSLE模型計算了位于越南北部地區的Lo River流域的懸移質；Tran等［12］利用基于模糊邏輯的模型改良了RUSLE模型；怡凱等［13］綜合利用RUSLE模型研究遼寧省朝陽市2001—2010年的土壤侵蝕時空分異特征。

2020年1月，習近平總書記再次來到云南，心系云南生態環境保護。而玉溪市正處云貴高原西部的滇中高原區，是滇中城市經濟圈南大門；玉溪市水系發達，且土壤侵蝕主要為水力侵蝕（水土流失等），玉溪市部分地區由于自然因素以及人為因素的影響，導致自然災害頻發（主要為滑坡、泥石流等），因此，基于玉溪市優質的地理環境條件以及土壤和水資源條件，根據重點開發的政策導向和整體的土壤侵蝕情況開展土壤侵蝕研究工作，綜合玉溪市土壤侵蝕分布情況以及造成土壤侵蝕主要因素，可以對該地區土地利用規劃與生態環境工程管理政策的制定提供一定的科學依據，為區域的可持續發展提供決策支持。

1 研究區域與方法

1.1 研究區概況

玉溪市位于云南省滇中地區，地勢西北高、東南低，哀牢山脈為境內主要山峰，其斜貫新平縣與元江縣，呈西北向東南走向。全市最高點大雪鍋山海拔為3 137 m，最低點海拔約為327 m，平均海拔在1 500?1 800 m。

1.2 研究過程與方法

考慮到玉溪市的基本情況，本研究采用修正的通用土壤侵蝕方程計算，公式：

式中，A為現實土壤侵蝕量［t/（hm2·a）］；R為降雨侵蝕力因子［（MJ·mm）/（hm2·h·a）］；K為土壤可蝕性因子［（t·hm2·h）/（hm2·MJ·mm）］；LS為地形坡度因子；C為地表植被覆蓋與管理因子；P為水土保持措施因子。

1.2.1 降雨侵蝕力因子降水作為引起水土流失的主要氣候因素，是造成土壤侵蝕的主要因素。考慮到模型的適用性以及考慮更容易獲得氣象因子數據，本研究采用Arnoldus在1977年提出的簡易計算方法，該公式只考慮了年降雨量以及月降水量，具有一定代表性而且應用范圍較廣，公式：

式中，R為年降雨侵蝕力［（MJ·mm）/（hm2·h·a）］，P為年平均降水量（mm），Pj為第j月降水量（mm）。降水量數據在中國氣象數據網獲得，計算得到玉溪市及周邊各氣象站點的年降雨侵蝕力，通過在Arc-GIS軟件中進行克里金插值，得到流域年降雨侵蝕力柵格圖，其中部分月份結果如圖1所示。

圖1 玉溪市單數月平均降水量

1.2.2 地形坡度因子地形坡度因子測算主要應用DEM數據。閱讀文獻可知，目前對于地形坡度因子已經有大量學者進行研究，本研究選用具有一定代表性且具有一定普適性的公式進行計算，利用DEM數據通過坡度提取、柵格計算器計算得到坡度因子（S），計算公式如下。

坡長因子表達式：

式中，L代表標準化坡長上的土壤侵蝕量；m代表坡長因子指數；θ為坡度（°）。地形坡度因子計算結果如圖2所示。

圖2 地形坡度因子

1.2.3 土壤可侵蝕因子土壤可侵蝕因子計算模型一般用EPIC模型，該模型主要參數為土壤的粉粒、沙粒、黏粒含量與土壤有機碳含量。在世界土壤數據庫（HWSD）中下載全球土壤數據，利用柵格計算器根據已有公式計算出土壤可侵蝕因子（K）。計算公式：

式中，S1為砂粒含量（%）；S2為粉粒含量（%）；n為黏粒含量（%）；C為有機碳含量（%）。土壤可侵蝕因子計算結果如圖3所示。

圖3 玉溪市土壤可侵蝕因子分布

1.2.4 植被覆蓋管理因子主要利用歸一化植被指數（NDVI）來計算植被覆蓋因子（C），再通過投影、重采樣和掩膜提取這些準備工作之后，利用公式計算植被覆蓋度和植被覆蓋因子。植被覆蓋管理因子如圖4所示。

圖4 玉溪市植被覆蓋管理因子分布

式中，ρNIR為近紅外波段，ρR為紅波段，NDVIMIN為NDVI的最小值，NDVIMAX為NDVI的最大值，f為植被覆蓋度，C為植被覆蓋管理因子。

1.2.5 水土保持措施因子水土保持措施因子（P因子）是采取水保措施后的土壤流失量與未實施水保措施的順坡耕作地塊的土壤流失量之比。提取玉溪市內耕地、林地、草地、裸地和建設用地、水體和濕地共6種類別，結合土地利用數據由DEM生成坡度信息，參照方廣玲等［14］的研究成果對P因子賦值，取值范圍為0?1，0代表不會發生土壤侵蝕的地區，1代表未采取任何控制措施的地區。賦值如表1所示，水土保持措施因子計算結果如圖5所示。

表1 不同土地利用類型P因子值

圖5 水土保持措施因子分布圖

1.3 數據收集

本研究主要采用玉溪市的降雨、土壤、DEM、植被和土地利用數據。在國家氣象科學數據中心下載玉溪市及周邊縣區月均降雨量數據，在世界土壤數據庫（HWSD）中獲取土壤數據；在地理空間數據云中下 載90 m分辨率DEM數據，并且 以Landsat 8的玉溪市遙感影像計算植被指數，土地利用數據則采用中國科學院資源環境科學與數據中心中提供的數據進行分析。根據上述方法代入公式，利用ArcGIS軟件進行處理，并計算得到RUSLE模型所需的因子空間分布的30 m×30 m柵格數據，最終得到研究區土壤侵蝕強度空間分布情況。

2 結果與分析

將各土壤侵蝕因子的柵格數據的分辨率統一定為30 m。根據修正的通用土壤侵蝕方程計算土壤侵蝕程度，并乘以美制轉公制常數100，將單位t/（hm2·a）轉換為t/（km2·a），得出研究區土壤侵蝕量柵格數據。玉溪市地處西南地區，并且以紅壤為主，因此參照水利部頒布的土壤侵蝕分類分級標準（SL 190—2007），該市的土壤侵蝕強度分為微度、輕度、中度、強烈、極強烈和劇烈6個級別，分級標準如表2所示。

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

通過計算最終得出研究區土壤侵蝕的分布情況（圖6），研究區內平均土壤侵蝕量為1 043.03 t/（km2·a），其中級別為微度的面積為863 497.62 hm2，占該市總面積的57.844%，級別為輕度的面積為499 884.03 hm2，占該市總面積的33.486%，級別為中度的面積為112 005.72 hm2，占該市總面積的7.503%，級別為強烈的面積為14 646.96 hm2，占該市總面積的0.981%，土壤侵蝕強度級別為極強烈及以上的區域面積為2 772.63 hm2，占該市總面積的0.186%。因此通過分析可知，玉溪市的土壤侵蝕主要表現為大部分區域輕度、微度侵蝕，局部區域強烈和極強烈等特征。

圖6 玉溪市土壤侵蝕強度分布特征

2.1 不同行政區內的土壤侵蝕特征分析

根據表3分析玉溪市各縣的土壤侵蝕量以及不同等級土壤侵蝕在各行政單元的分布情況，全市內主要土壤侵蝕級別為微度和輕度，玉溪市整體土壤侵蝕嚴重區域集中分布在西南部以及東南部，主要在元江、綠汁江以及曲江流域附近，由于哀牢山國家自然保護區位于新平縣，且該地正處于東亞季風熱帶、南亞季風熱帶和青藏橫斷山系三大自然地理區域的交匯處，當氣流運動到這里后，氣溫降低，有利于凝結降水，因此容易發生水土流失，新平縣的平均土壤侵蝕量為799.079 8 t/（km2·a），其土壤侵蝕中度及以上占全縣的6.720 6%，而元江縣的平均土壤侵蝕量為969.908 7 t/（km2·a），全縣內土壤侵蝕級別大多數為微度和輕度。玉溪市水系發達，縱觀整體可知該市的土壤侵蝕較為嚴重的區域集中分布在各大河流附近，例如新平縣的土壤侵蝕中度以上區域主要分布在漠沙江附近，易門縣、峨山縣土壤侵蝕中度以上區域主要集中在綠汁江流域，根據《截至2015年底云南省礦產資源儲量簡表》可知，玉溪市礦山主要分布在易門縣以及華寧縣，而礦山的開采破壞土壤結構，因此，易門縣和華寧縣的土壤侵蝕量相對較高，分別為571.142 1 t/（km2·a）和777.671 2 t/（km2·a）。

表3 各行政單元土壤侵蝕情況

2.2 不同土地利用類型的土壤侵蝕特征分析

結合遙感影像以及實地調查，選取玉溪市內耕地、林地、草地、工礦用地、水域及水利設施用地進行研究。由表4可知，玉溪市土壤侵蝕較為嚴重的地區主要地類為耕地、林地、草地，中度及以上侵蝕面積比例超過35%，玉溪市耕地面積263 804.44 hm2，占玉溪市總面積的17.66%，由于耕地受人類活動干擾頻繁，常年耕作導致土壤結構遭到破壞，作物生長周期的原因導致植被覆蓋具有季節變化特征，因此土壤侵蝕也相對嚴重。由于玉溪市局部地區存在植被退化的現象，嚴重影響了林地以及草地的植被覆蓋，導致了這些區域的侵蝕強度較高，其中全市林草地中度以上侵蝕面積近30%。玉溪市礦產資源較為豐富，已知礦產有鐵、錳、鉻、銅、鎳、磷、大理石、煤和硅石等23種，其中鐵、銅、磷、鎳4種優勢礦產是該區重要的礦產資源，玉溪市工礦用地占地11 060.78 hm2，其中輕度以上土壤侵蝕占29.24%。由于地區工業礦床多且分布廣泛，大量長期的煤礦開采造成表層物質松散，土地植被以及地表、地下含水巖層遭到破壞，與其他形式的土地破壞相比，具有形式特殊、數量大、難補救的特點［13］。

表4 不同土地利用類型土壤侵蝕情況

2.3 不同植被覆蓋度土壤侵蝕狀況

現將植被覆蓋程度分為5個等級，以便于分析玉溪市不同植被覆蓋條件下的土壤侵蝕狀況。由表5可知，植被覆蓋度小于10%的區域，土壤侵蝕等級主要為微度和輕度，土壤侵蝕程度微度的面積占59.30%，輕度土壤侵蝕的面積占36.34%；植被覆蓋度為10%?30%區域，土壤侵蝕程度微度的面積占41.45%，輕度土壤侵蝕的面積占48.53%；植被覆蓋度為45%?60%區域，受土壤侵蝕面積較大的影響，其土壤侵蝕量也不容忽視。對不同植被覆蓋條件下的區域進行分析（圖7），隨著植被覆蓋度增加，輕度以下侵蝕面積增加，絕大部分中度以上侵蝕發生在植被覆蓋度低于60%的區域。

表5 不同植被覆蓋程度土壤侵蝕情況（單位：%）

圖7 不同植被覆蓋度各土壤侵蝕強度等級面積比例

3 小結與討論

利用GIS與RUSLE模型對玉溪市的土壤侵蝕程度進行研究，結果表明，研究區的平均土壤侵蝕量為1 043.03 t/（km2·a），該市土壤侵蝕主要以微度和輕度為主，其所占面積分別為57.844%和33.486%，土壤侵蝕強烈及以上部分主要分布在易門縣、華寧縣以及元江縣。根據土地利用類型與土壤侵蝕程度結合分析，玉溪市土壤侵蝕較為嚴重的地類基本分布在耕地、林地以及草地，侵蝕強度為中度及以上的區域侵蝕面積比例超過35%，雖然林地和草地的植被覆蓋較高也具有一定水土保持作用，但是植被退化導致了局部區域侵蝕強度較高，全市林草地中度以上侵蝕面積近30%。而玉溪市礦產資源豐富，采礦對土壤的傷害仍不容忽視。現有的土壤侵蝕預測模型大多是地方性的，不能用于較大區域，也很難直接推廣到其他地區，同時對各侵蝕影響因子尚無系統全面的研究，缺乏足夠的試驗觀測資料，影響了模型的實用性［15］，而本研究缺乏對相關參數的實地驗證，在后續研究中予以完善；另外模型因子計算方法應進一步完善，提高精度。