安徽省1980—2020年土壤侵蝕時空變化特征

徐少杰, 鄧良, 趙明松*, 陳宣強, 邱士其

(1.安徽理工大學空間信息與測繪工程學院，淮南 232001；2.礦山采動災害空天地協同監測與預警安徽省教育廳重點實驗室，淮南 232001； 3.礦區環境與災害協同監測煤炭行業工程研究中心， 淮南 232001； 4.安徽省地質測繪技術院，合肥 230022)

土壤侵蝕嚴重威脅著人類的生存和發展，在水力、風力或凍融等外營力的作用下，土壤侵蝕會導致土壤及其母質被破壞、搬運以及堆積，造成土壤的退化、肥力下降以及山體滑坡等自然災害的加劇，嚴重制約著社會經濟的發展[1]。隨著對土壤侵蝕的不斷關注，研究土壤侵蝕的方法也在不斷被完善。從19世紀40～50年代開發并得到實踐應用的通用土壤流失方程(universal soil loss equation，USLE)[2]、修正的通用土壤流失方程(revised universal soil loss equation，RUSLE)、水蝕預測項目(water erosion prediction project，WEPP)等土壤侵蝕模型[3]，到如今經過各國學者研究改正的中國土壤流失方程(Chinese soil loss equation，CSLE)、生態系統服務和交易的綜合評估(integrated valuation of ecosystem services and trade-offs，InVEST)[4]、生態系統服務的人工智能(artificial intelligence for ecosystem services，ARIES)等模型。

應用不同的模型，探討不同區域的土壤侵蝕特征為當前土壤侵蝕研究的主要熱點之一。如Wang等[5]運用InVEST模型，分析了粵北地區在2000—2008年內水土流失的時空演變特征。李娜等[6]基于CSLE模型，利用采樣調查與地理信息系統(geographic information system，GIS)相結合，探討了天山北坡西白楊溝流域的土壤侵蝕強度變化的驅動機制及其時空變化。影響土壤侵蝕的因素較多，應根據當地研究區的實際情況，采取不同的影響因子以及適宜的模型進行土壤侵蝕的研究。陳起偉等[7]基于土地利用類型、巖性、基巖裸露率等因子，針對貴州省喀斯特地區進行了土壤侵蝕及土壤抗侵蝕年限的研究，為喀斯特區土壤侵蝕的防治提供了參考。曾成等[8]基于RUSLE模型對喀斯特峰叢洼地的土壤侵蝕展開研究，并探討了土壤侵蝕與土壤養分之間的關系，結果表明：人類活動是造成該地區土壤侵蝕及養分流失的重要因素，為喀斯特地區生態環境治理提供了重要的數據支撐。RUSLE模型因其普適性，被廣泛應用于土壤侵蝕的研究之中，為合理規劃黃河流域在蘭州段的水土資源，武翠翠等[9]基于RUSLE模型對黃河流域蘭州段的土壤侵蝕敏感性進行定量評價，該地區土壤侵蝕敏感性面積占比較大。肖武等[10]利用USLE模型分析了安徽省巢湖流域2000—2010年的土壤侵蝕空間變化趨勢，并識別出土壤侵蝕的關鍵區域。陳朝良等[11]基于RUSLE模型在黃河上游湟水流域開展土壤侵蝕研究，結果表明：2000—2015年該流域土壤侵蝕強度整體呈減弱趨勢。研究表明，RUSLE等模型在平原以及河谷平原地區具有較高的適用性。

在安徽省的土壤侵蝕的研究之中，大多數均為小區域或短時間段內的土壤侵蝕時空演變分析。如管禮松等[12]基于CSLE模型精確評估了安徽省皖南區域2010—2018年的土壤侵蝕狀況。趙明松等[13]定量分析了安徽省大別山區域的土壤侵蝕空間分布特征，并探討了地形因素對土壤侵蝕的影響。但當下卻很少有人對安徽省近年來的土壤侵蝕時空變化進行完整的研究和分析。因此，現以安徽省為研究區，利用1980年、2000年、2005年、2010年、2015年、2020年6個時期的降雨量、土地利用現狀圖、地形數據、土壤類型等數據，借助GIS分析探討安徽省40年以來土壤侵蝕的時空演變，并著重探討土壤侵蝕和土地利用類型之間的關系，為今后安徽省土壤流失的治理保護提供重要依據。

1 材料與方法

1.1 研究區概括

安徽省位于中國華東地區，地處長江、淮河中下游以及長江三角洲腹地。介于114°54′E～119°37′E，29°41′N～34°38′N，土地面積約為13.94萬km2。屬暖溫帶與亞熱帶的過渡地區，年平均氣溫為14～17 ℃，年平均降水量為773～1 670 mm[14]。全省分成淮北與沿淮平原、沿江平原、江淮丘陵崗地、皖西大別山區、皖南丘陵山地區5個地貌區，全省地勢西南高、東北低，平均海拔119.3 m。土壤類型主要為潮濕雛形土、水耕人為土、濕潤淋溶土、正常新成土、正常有機土等。安徽省數字高程模型及地理分區如圖1所示。

圖1 安徽省數字高程模型及地理分區Fig.1 Digital elevation model and geographical division of Anhui Province

1.2 數據來源

(1)安徽省1980年、2000年、2005年、2010年、2015年、2020年的月平均降雨量數據，來源于安徽省氣象信息中心及安徽省統計年鑒。

(2)安徽省土壤圖以及《安徽土種》，比例尺為1∶50萬，來源于安徽省土壤普查辦公室。

(3)90 m空間分辨率的數字高程模型SRTM DEM(https://srtm.csi.cgiar.org/SELECTION/inputCoord.asp)。

(4)安徽省1980年、2000年、2005年、2010年、2015年、2020年土地利用圖，由Landsat TM影像圖監督分類而來，空間分辨率為30 m，分類精度均達到了88.9%。

(5)2010年安徽省典型土壤類型理化性質，來源于國家科技基礎性工作專項“我國土系調查與《中國土系志》編制”中的安徽省土系調查的206個典型樣點數據集。

1.3 RUSLE模型因子確定

基于RUSLE模型以及GIS技術評價安徽省40年來土壤侵蝕的時空變化特征，其公式為

A=RKLSCP

(1)

式(1)中：A為單位面積上的年平均土壤侵蝕模數，t/(hm2·a);R為降雨侵蝕力因子， MJ·mm/(hm2·h·a),采用周伏建等[15]的降雨侵蝕力因子計算公式，利用安徽省各站點的月平均降雨量分別計算每個站點的R因子，插值生成各時期R值分布圖；K為土壤可蝕性因子，t·h/(MJ·mm);基于Williams[16]在EPIC模型中的估算方法計算K值，并根據張科利等[17]的研究對K值進行修正；L、S分別為坡長、坡度因子，L因子采用Wischmeier等[18]提出的公式進行計算，S因子在緩坡(坡度≤14°)采用Liu等[19]的算法進行計算，而在陡坡(坡度>14°)則采用McCool等[20]提出的公式；C為植被覆蓋與管理因子[21]。根據前人的研究[22-23]，為不同的土地利用類型賦以不同的C值:旱地為0.31，建設用地為0.2，水田為0.18，草地為0.06，疏林地為0.017，其他林地為0.006,水域和裸地均為0；P為水土保持措施因子。根據以往學者的相關研究[24-26]，結合安徽省的土地利用圖，完成水土保持措施因子的賦值:旱地為0.35，水田為0.15，草地、疏林地及林地均為1，水域、裸地和建設用地均為0。

按照水利部《土壤侵蝕分類分級標準》(SL190—2007)，根據土壤侵蝕模數將安徽省土壤侵蝕等級劃分為6類：微度侵蝕[500 t/(km2·a)以下]、輕度侵蝕[500～2 500 t/(km2·a)]、中度侵蝕[2 500～5 000 t/(km2·a)]、強度侵蝕[5 000～8 000 t/(km2·a)]、極強度侵蝕[8 000～15 000 t/(km2·a)]以及劇烈侵蝕[15 000 t/(km2·a)以上]。

2 結果與分析

2.1 安徽省土壤侵蝕量變化

圖2是安徽省1980、2000、2005、2010、2015、2020年的土壤侵蝕強度空間分布圖。表1為安徽省土壤侵蝕量的分類及統計情況。可以發現，安徽省土壤侵蝕整體呈現南高北低的分布特征。大部分區域為微度及輕度侵蝕，中度及以上侵蝕主要分布在安徽省西部和南部，北部以微度、輕度侵蝕為主。

圖2 安徽省1980—2020年土壤侵蝕強度空間分布Fig.2 Spatial distribution of soil erosion intensity in Anhui Province from 1980 to 2020

表1 安徽省1980—2020年各土壤侵蝕強度的土壤侵蝕量

圖3 安徽省土壤侵蝕強度變化空間分布Fig.3 Spatialdistribution of soil erosion intensity change in Anhui Province

自1980年，安徽省的土壤侵蝕總量及平均土壤侵蝕模數總體呈現先減少后增加的趨勢。1980年安徽省土壤侵蝕總量達到6 510.81 萬t/a，2010年土壤侵蝕總量下降至3 450.88 萬t/a，之后土壤侵蝕有所加劇，2020年土壤侵蝕總量上升至4 202.89 萬t/a。與1980年相比，2000年、2005年、2010年、2015年、2020年土壤侵蝕總量分別下降了23.82%、51.32%、47.00%、41.70%、35.45%，但相比于2010年，2015年和2020年的土壤侵蝕總量有所增加。

2.2 安徽省土壤侵蝕面積變化

圖3為安徽省1980—2020年間6個時間段內的土壤侵蝕強度空間變化圖。在1980—2000年、2000—2005年各地區土壤侵蝕狀況呈現減弱的趨勢，侵蝕強度變化的區域主要集中在皖西大別山區、皖南丘陵山地區。相比于2005年的土壤侵蝕狀況，2010年、2015年、2020年皖西大別山區、皖南丘陵山地區土壤侵蝕呈現加強的趨勢。

圖4為1980—2020年各等級侵蝕強度的面積轉移流向圖。中度以上侵蝕等級占比較少，因此，將強度侵蝕、極強度侵蝕以及劇烈侵蝕歸為一類，即中度以上侵蝕。

圖4 安徽省1980—2020年侵蝕強度面積轉移桑基圖Fig.4 Sankey plot of erosion intensity area transfer from 1980 to 2020 in Anhui Province

安徽省40年來土壤侵蝕強度有所減弱，在2005年及以前，土壤治理措施非常有成效，微度侵蝕面積呈現不斷增加的趨勢，有較大部分的輕度侵蝕向微度侵蝕進行轉移，較大面積的土壤侵蝕得到治理及改善。在2005年以后有所回升，但回升趨勢較低，但相比于1980年來說，總體上還是處于降低趨勢。從圖4中可以發現，相比于1980年，2020年安徽省微度侵蝕面積有所增加，而輕度、中度及中度以上的侵蝕面積均得到了一定程度的降低。

3 討論

3.1 降雨量變化對土壤侵蝕的影響

在安徽省近些年里，特別是在2020年，安徽省降雨頻繁，為探究在降雨“驟增”的情況下，安徽省土壤侵蝕變化情況，對比分析安徽省1980—2020年土壤侵蝕總量與年平均降雨量，圖5。

圖5 安徽省1980—2020年侵蝕量與年平均降雨對比Fig.5 Comparison of the amount of erosion and annual average rainfall in Anhui Province from 1980 to 2020

在2005年，安徽省總侵蝕量達到最低，在降雨量相近的情況下，安徽省一系列的土壤治理措施使得安徽省土壤侵蝕得到了較大的改善，總土壤侵蝕量得到了較大幅度的降低；相比于2005年降雨，2010年、2015年、2020年年平均降雨量有了較大幅度的增長。特別是2020年，相比于2015年，年平均降雨曲線出現了明顯的“陡坡”，雖然一定的水土保持措施，減緩了水土流失，但降雨量的暴增使得全省2010年、2015年、2020年的土壤侵蝕總量有所增加。

盡管2020年強降雨以及洪澇災害頻發，但安徽省40年來的土壤治理措施使得土壤抗侵蝕能力得到了大大的提升，侵蝕量并沒有因為年降雨平均量的“驟增”而迅速提升。安徽省近年來開展的一系列水土保持與治理工程，例如：自2010年，安徽省每年平均審批生產建設項目水土保持方案近500個。2019年達到了837個，2019年內治理小流域24條，治理水土流失493.8 km2，使水土保持監測與保護效率進一步提高[27-28]。

3.2 不同土地利用類型下土壤侵蝕強度差異

根據安徽省土壤侵蝕變化特征，以2005年為節點，統計1980—2005年、2005—2020年、1980—2020年間的土地利用情況，對不同土地利用類型的土壤侵蝕情況進行分析。

3.2.1 各土地利用類型土壤侵蝕強度變化

表2反映了各土地利用類型(剔除水體)在1980年、2005年、2020年的土壤侵蝕平均侵蝕模數及總侵蝕量信息。可以發現，整體上草地的侵蝕模數和侵蝕量相對較大；建設用地以及裸地的平均侵蝕模數及侵蝕量最低。同時對比各土地利用類型的侵蝕模數標準差，草地平均侵蝕模數的標準差加大，這表示在草地區域內土壤侵蝕模數浮動較大，同時與平均侵蝕模數有較大的差異；建筑用地及裸地在2005年以前標準差均維持在1以下，但在2020年標準差有所提升，參照兩類區域內總侵蝕量占比較少，可以得知，建筑用地及裸地區域的土壤侵蝕僅存在于個別區域。這與程先富等[22]的研究結果相符，草地區域的侵蝕最為嚴重，林地次之，而建筑用地和裸地區域的侵蝕情況較弱。

表2 各類土地利用類型土壤侵蝕概況

3.2.2 土地利用變化對侵蝕強度的影響

進一步分析土地利用的改變對土壤侵蝕強度的影響，由于建筑用地及裸地的侵蝕量較少，僅統計1980—2020年水田、旱地、林地及草地區域的土地利用面積轉移以及土壤侵蝕強度的變化情況，見表3。

表3 土地利用改變區域土壤侵蝕強度變化面積

土地利用類型發生改變的區域占據了較少的一部分，較大面積的土地利用未發生變化，未發生土地利用轉移的區域的土壤侵蝕強度仍以減弱為主，由此可知，安徽省的一系列治理措施取得了較高的成效。在土地利用發生改變的區域，除了林地，其他3類用地的侵蝕強度減弱的面積均比侵蝕強度增加的面積多。在水田區域，侵蝕強度得到降低的區域集中在水田-林地。旱地區域同樣主要集中在旱地-林地的區域。同時草地區域也以侵蝕強度減弱為主，在草地-林地的區域有175.36 km2的面積侵蝕強度得到減弱；而林地則有較大面積的土壤侵蝕強度得到加劇，主要集中在林地-草地區域。當地表裸露的旱地、草地及水田轉移為能有效阻擋降雨動能的林地后，有效地遏制了土壤侵蝕，而當林地向其他土地利用類型轉移時，旱地、草地等區域無法抵擋雨水的沖刷，從而導致土壤侵蝕強度加劇。

4 結論

(1)安徽省土壤侵蝕強度呈現南高北低的趨勢，整體以微度侵蝕為主。

(2)1980—2005年，安徽省整體土壤侵蝕呈現減弱的趨勢，在2015—2020年土壤侵蝕因降雨略有回升。侵蝕強度變化的區域主要集中在皖西大別山區、皖南丘陵山地區。

(3)旱地、水田和草地轉移為林地能有效地降低土壤的侵蝕強度，特別是在侵蝕較為嚴重的草地區域，應根據實際情況加強林地的建設。

安徽省一系列的水土管理與治理措施使土壤侵蝕得到了一定的改善。但在安徽南方，由于降雨量較多導致南方土壤侵蝕均較為嚴重，仍存在著較多的強度及以上的侵蝕，今后應重點加強皖西大別山區、皖南丘陵山地區的水土防護和治理。