黃碩文，李健*，張欣佳，鄧聯(lián)文，張金萍

（1.鄭州大學(xué)地球科學(xué)與技術(shù)學(xué)院，鄭州 450001；2.鄭州大學(xué)水利科學(xué)與工程學(xué)院，鄭州 450001）

2019年9月，黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展被列為重大國家戰(zhàn)略，水土流失相關(guān)研究是生態(tài)保護(hù)的重要組成部分。河南省處于黃河中下游，是農(nóng)業(yè)大省和糧食主產(chǎn)區(qū)，但部分地區(qū)在氣候、地形等因素的綜合作用下極易發(fā)生土壤侵蝕，對耕地產(chǎn)生破壞，直接影響糧食作物的產(chǎn)量。為了更好地保護(hù)有限的耕地資源，實現(xiàn)區(qū)域高質(zhì)量發(fā)展，有必要對河南省土壤侵蝕狀況進(jìn)行定量評價，分析其時空變化特征，探究引發(fā)土壤侵蝕的主導(dǎo)因素，并制定相關(guān)的保護(hù)措施。

對土壤侵蝕進(jìn)行評價主要有實地觀測和建立土壤侵蝕模型兩種方法。在實際應(yīng)用中，由于某些地區(qū)的客觀條件限制難以進(jìn)行實地觀測，建立土壤侵蝕模型進(jìn)行侵蝕量預(yù)報是常用的方法。根據(jù)建立的方法用途和所模擬的侵蝕過程可以將土壤侵蝕模型分為經(jīng)驗?zāi)Ｐ?、物理模型和分布模型[1]。代表性的物理模型有水蝕預(yù)報模型WEPP[2]、歐洲土壤侵蝕預(yù)報模型EUROSEM[3]、非點(diǎn)源地區(qū)流域環(huán)境反應(yīng)模型ANSWERS[4]等，分布模型以SHE 模型[5]最為典型。物理模型和分布模型由于面向過程建模，所需的資料和參數(shù)設(shè)置過多，難以在特定區(qū)域外進(jìn)行推廣。土壤侵蝕經(jīng)驗?zāi)Ｐ突诖罅繉崪y數(shù)據(jù)進(jìn)行統(tǒng)計分析，結(jié)合土壤侵蝕過程的機(jī)理建立，國際上應(yīng)用較廣的主要有通用水土流失方程USLE 和其修正形式RUSLE。RUSLE 對USLE 中各因子含義及計算方法進(jìn)行了必要的修正，增加了復(fù)雜坡度下土壤流失量的計算方法，具有結(jié)構(gòu)簡明、計算因子含義明顯的特點(diǎn)，且容易與GIS 結(jié)合進(jìn)行分析，在國內(nèi)外土壤侵蝕定量評價中應(yīng)用廣泛[6-9]。土壤侵蝕過程中同時受到多種因素的影響，因此分析其內(nèi)在驅(qū)動力對于水土保持政策的制定至關(guān)重要。地理探測器可用于分析不同自變量對因變量的解釋能力，在土壤侵蝕領(lǐng)域的應(yīng)用也得到了學(xué)者們的認(rèn)可[10-11]。

本研究利用RUSLE模型與GIS技術(shù)相結(jié)合，對河南省全境土壤侵蝕強(qiáng)度進(jìn)行定量評價，分析其時空變化情況，并利用地理探測器分析土壤侵蝕影響因素，為河南省水土流失防治重點(diǎn)保護(hù)區(qū)的劃定和生態(tài)環(huán)境保護(hù)工作提供參考。

1 材料與方法

1.1 研究區(qū)概況

河南?。?1°23′～36°22′N，110°21′～116°39′E）地處中原地區(qū)，是國家重要的糧食生產(chǎn)基地（圖1），由于季風(fēng)氣候的影響以及跨海河、黃河、淮河、長江四大水系的地理位置，研究區(qū)內(nèi)常常出現(xiàn)局部土壤強(qiáng)侵蝕，對生態(tài)環(huán)境造成破壞，但近幾年河南加強(qiáng)了水土保持措施的制定與實施，水土流失情況得到了一定改善。

圖1 河南省行政區(qū)劃及數(shù)字高程模型Figure 1 Administrative division and digital elevation model of Henan Province

1.2 數(shù)據(jù)來源

本研究所選用的數(shù)據(jù)為河南省行政區(qū)劃圖、降雨數(shù)據(jù)（2008、2013、2018年）、ETM+衛(wèi)星遙感數(shù)據(jù)、數(shù)字高程模型（DEM）數(shù)據(jù)、地貌類型數(shù)據(jù)和土壤數(shù)據(jù)。其中降雨數(shù)據(jù)為中國氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn）的中國地面氣候資料日值數(shù)據(jù)集；土壤數(shù)據(jù)（質(zhì)地、有機(jī)質(zhì)等）來源于黑河計劃數(shù)據(jù)管理中心（http：//westdc.west?gis.ac.cn）的“基于世界土壤數(shù)據(jù)庫的中國土壤數(shù)據(jù)集（V1.1）”，從中截取了河南省地區(qū)數(shù)據(jù)；DEM 90 m分辨率原始高程數(shù)據(jù)和ETM+衛(wèi)星遙感影像來源于地理空間數(shù)據(jù)云平臺（http：//www.gscloud.cn）；地貌形態(tài)類型數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境數(shù)據(jù)中心（http：//www.resdc.cn）。將各個因子數(shù)據(jù)坐標(biāo)系轉(zhuǎn)化為高斯-克里格投影坐標(biāo)系，柵格數(shù)據(jù)分辨率統(tǒng)一為90 m×90 m。

1.3 RUSLE模型

RUSLE模型充分考慮了影響土壤侵蝕的多個因子，對復(fù)雜、分布廣泛的土壤侵蝕及模型參數(shù)變異性等問題均提出較好的解決方案，其公式為：

式中：A為模型計算得到的土壤侵蝕模數(shù)，t·hm-2·a-1，本研究乘以100，將單位換算為t·km-2·a-1；R為降雨侵蝕力因子；K為土壤可侵蝕性因子；L為坡長因子；S為坡度因子；C為作物覆蓋和管理因子；P為水土保持措施因子[12]。RUSLE 模型通過計算各個因子的指標(biāo)值對土壤侵蝕強(qiáng)度進(jìn)行定量評價。

降雨侵蝕力因子（R）：降雨通過對土壤進(jìn)行沖擊引發(fā)水力侵蝕，是進(jìn)行土壤侵蝕預(yù)報的重要因子之一。本研究采用基于月平均降雨量和年平均降雨量數(shù)據(jù)的Wischmeier公式[13]進(jìn)行降雨侵蝕力的計算，公式為：

式中：R為年降雨侵蝕力，MJ·mm·hm-2·h-1；Pi、P分別為月平均、年平均降雨量，mm。

土壤可蝕性因子（K）：土壤可蝕性因子反映了土壤自身理化性質(zhì)對于土壤侵蝕的抵抗能力。同等條件下，K值越大，土壤被沖蝕的可能性就越大。本研究根據(jù)已有參數(shù)采用Williams 等[14]在EPIC 模型中的算法，利用土壤機(jī)械組成和土壤有機(jī)質(zhì)進(jìn)行K值計算。具體計算公式為：

式中：Sn=1-Sd/100；Sd為砂粒含量，%；Si為粉粒含量，%；Cl為黏粒含量，%；C為有機(jī)質(zhì)含量，%。

地形因子（LS）：地形因子直接影響坡面侵蝕的速率，是土壤侵蝕的動力因子。在RUSLE 模型中通常將坡長L因子與坡度S因子合并為LS 地形因子進(jìn)行計算。國內(nèi)學(xué)者基于不同地貌類型區(qū)和徑流小區(qū)的實測數(shù)據(jù)對模型中地形因子計算參數(shù)進(jìn)行了修正[15]，本研究采用第四次土壤侵蝕普查算法[16]計算LS因子，公式如下：

表1 不同土地利用類型的水土保持因子（P值）Table 1 Water and soil conservation factors（P values）of different land use types

式中：λ為坡長；θ為坡度；m為坡長系數(shù)。

植被覆蓋和管理因子（C）：植被覆蓋與管理因子（0≤C≤1）反映了植被覆蓋對土壤侵蝕的抑制作用，為無量綱數(shù)，其值越大則植被發(fā)揮減蝕的作用越小。C因子的值通常通過植被覆蓋度計算，本研究使用歸一化植被指數(shù)NDVI 來計算地表植被覆蓋度c，計算公式如下：

C因子的計算采用蔡崇法等[17]算法的修正形式，避免了當(dāng)c值在0～0.1之間時異常值（C＞1）的出現(xiàn)，如式（6）所示：

水土保持因子（P）：水土保持措施因子P（0≤P≤1），其值越小代表水土保持措施越好，土壤侵蝕越不易發(fā)生。目前P的賦值主要由實地考察資料結(jié)合經(jīng)驗決定，本研究參照了美國農(nóng)業(yè)部手冊537 號[18]和覃杰香等[19]的相關(guān)研究，對不同土地利用類型區(qū)域和不同坡度的耕地分別按表1和表2進(jìn)行賦值。

1.4 地理探測器

地理探測器由因子探測器、風(fēng)險探測器、生態(tài)探測器和交互作用探測器四部分組成。它通過探測要素的空間分層異質(zhì)性來揭示其背后驅(qū)動力[20]。

因子探測器可用于探測某個因素X對因變量Y的解釋程度，其大小可用q值來度量（0≤q≤1），其值越大表示自變量X對因變量Y的解釋力越強(qiáng)[20]。以此可以識別出影響土壤侵蝕的主導(dǎo)因素。

表2 不同坡度耕地的水土保持因子（P值）Table 2 Water and soil conservation factors（P values）of cultivated land on different slopes

交互探測器用于識別不同自變量因子之間的交互作用[20]，其交互方式見表3。

表3 自變量對因變量的交互作用方式Table 3 Types of interaction between two covariates

生態(tài)探測器可以對兩個因子進(jìn)行比較，判斷它們的空間分布影響是否有顯著差異[20]。風(fēng)險探測器判斷單個影響因子不同層間屬性是否有顯著差異，可用于識別土壤侵蝕高風(fēng)險區(qū)。

本研究選取土壤侵蝕強(qiáng)度作為因變量，坡度、海拔、年平均降雨量、植被覆蓋度、土地利用類型、地形地貌作為土壤侵蝕的影響因子，并對6 類因子進(jìn)行離散化處理。多年平均降雨量等間距分成9 類，植被覆蓋度分為8 類（＜0.3、0.3～0.4、0.4～0.5、0.5～0.6、0.6～0.7、0.7～0.8、0.8～0.9、0.9～1），坡度分為8 類（＜5°、5°～10°、10°～15°、15°～20°、20°～25°、25°～30°、30°～35°、＞35°），海拔分為4 類（＜500 m、500～1 000 m、1 000～1 500 m、＞1 500 m），土地利用類型分為林地、草地、水域、居民用地、建筑用地、未利用地6 個類別，地形地貌數(shù)據(jù)使用《中華人民共和國地貌圖集（1∶100 萬）》中的類別編號。然后以河南省為范圍進(jìn)行隨機(jī)采樣，共選取10 000個樣本點(diǎn)，每個樣本點(diǎn)賦予以上各類屬性值作為地理探測器的輸入數(shù)據(jù)。

圖2 河南省2008—2018年降雨侵蝕力因子分布Figure 2 Distribution of rainfall erosivity factors in Henan Province from 2008 to 2018

2 結(jié)果與討論

2.1 RUSLE模型土壤侵蝕分布計算結(jié)果

RUSLE 模型中降雨侵蝕力因子、土壤可侵蝕性因子、地形因子、植被覆蓋與管理因子、水土保持因子計算結(jié)果如圖2～圖6所示。

將以上各因子結(jié)果相乘，得到河南省2008—2018年的土壤侵蝕分布情況（圖7）。

圖3 河南省土壤可侵蝕性因子分布Figure 3 Distribution of soil erosibility factors in Henan Province

圖4 河南省地形因子分布Figure 4 Distribution of topographic factors in Henan Province

2.2 土壤侵蝕時空變化特征

根據(jù)《土壤侵蝕分類分級標(biāo)準(zhǔn)》（SL 190—2007），并結(jié)合河南省土壤侵蝕量實際分布情況對研究區(qū)土壤侵蝕強(qiáng)度進(jìn)行分級：平均侵蝕模數(shù)0～500 t·km-2·a-1劃分為微度侵蝕區(qū)；平均侵蝕模數(shù)500～2 500 t·km-2·a-1劃分為輕度侵蝕區(qū)；平均侵蝕模數(shù)2 500 t·km-2·a-1以上區(qū)域統(tǒng)一劃分為中度及以上侵蝕區(qū)。河南省各類侵蝕區(qū)的面積占比見表4。

圖7 和表4 顯示，河南省2008—2018 年土壤侵蝕分布廣、程度低：整體以微度侵蝕為主，其面積占比超過95%，且2008 年以來逐年遞增，分布在東部及中部沖積平原區(qū)；輕度侵蝕地區(qū)主要分布在西部、南部山區(qū)及中部風(fēng)沙平原區(qū)，其面積占比呈逐年下降趨勢，從圖7 可以看出豫西部地區(qū)下降尤為顯著；中度及以上侵蝕區(qū)域集中在西北及北部黃河流域，多年來總體變化不大。

圖5 河南省2008—2018年植被覆蓋與管理因子分布Figure 5 Distribution of vegetation cover and management factors in Henan Province from 2008 to 2018

圖6 河南省2008—2018年水土保持因子分布Figure 6 Distribution of water and soil conservation factors in Henan Province from 2008 to 2018

圖7 河南省2008—2018年土壤侵蝕模數(shù)分布Figure 7 Distribution of soil erosion modulus in Henan Province from 2008 to 2018

表4 2008—2018年河南省不同土壤侵蝕等級面積比例（%）Table 4 Proportion of area of different soil erosion levels in Henan Province from 2008 to 2018（%）

對各年份土壤侵蝕分級圖進(jìn)行疊加分析得到2008—2013、2013—2018 年兩個階段的土壤侵蝕強(qiáng)度轉(zhuǎn)移矩陣（表5、表6）和轉(zhuǎn)移圖（圖8）。

根據(jù)圖表顯示，2008—2013年河南省共有2 738.01 km2地區(qū)土壤侵蝕強(qiáng)度降低，243.88 km2侵蝕強(qiáng)度增加。侵蝕強(qiáng)度降低主要集中在輕度轉(zhuǎn)化為微度，轉(zhuǎn)化面積達(dá)2 324.77 km2。從前文分析可知侵蝕下降主要集中在西部地區(qū)，這與河南省近年來在西部山地生態(tài)區(qū)實施的一系列保護(hù)性措施有關(guān)。這些措施增加了植被面積，有效抑制了土壤侵蝕。侵蝕強(qiáng)度增加主要集中在微度侵蝕增強(qiáng)為輕度，主要是因為豫西北地區(qū)以山地為主，同時受河流侵蝕及耕地分布不合理等自然及人為因素的共同影響，治理難度較大。

2013—2018 年河南省共有946.57 km2地區(qū)土壤侵蝕強(qiáng)度減小，843.89 km2地區(qū)侵蝕強(qiáng)度加大。與2008—2013 年相比，侵蝕強(qiáng)度明顯加強(qiáng)。這可能是由于2018年河南省降雨量明顯增加，水力侵蝕增強(qiáng)。

表5 河南省2008—2013年土壤侵蝕強(qiáng)度轉(zhuǎn)移矩陣（km2）Table 5 Transfer matrix of soil erosion intensity in Henan Province from 2008 to 2013（km2）

表6 河南省2013—2018年土壤侵蝕強(qiáng)度轉(zhuǎn)移矩陣（km2）Table 6 Transfer matrix of soil erosion intensity in Henan Province from 2013 to 2018（km2）

2.3 土壤侵蝕影響因素分析

前述研究結(jié)果揭示了河南省土壤侵蝕的時空演變特征，借助地理探測器進(jìn)一步探討其驅(qū)動因子。

圖8 土壤侵蝕動態(tài)轉(zhuǎn)移分布圖Figure 8 Dynamic transfer map of soil erosion

表7 土壤侵蝕影響因子的解釋力（q值）Table 7 The q value of influencing factors of soil erosion

地理探測器中因子探測器結(jié)果如表7 所示，不同影響因子對土壤侵蝕強(qiáng)度的解釋力有明顯差異，排序為坡度＞土地利用類型＞地貌類型＞年平均降雨量＞海拔＞植被覆蓋度。其中坡度對土壤侵蝕強(qiáng)度的解釋力最強(qiáng)，高達(dá)39.63%，可以看作是影響河南省土壤侵蝕的主導(dǎo)因子；其次是土地利用類型，解釋力為11.27%；植被覆蓋度的q值最小，解釋力不足1%，單因子作用下對土壤侵蝕的影響并不顯著。

根據(jù)交互探測器結(jié)果（表8），不同影響因子通過交互作用對土壤侵蝕空間分布的解釋力遠(yuǎn)大于單因子。坡度與其他因子的交互作用明顯大于剩余因子間的交互作用，其中與土地利用類型的協(xié)同作用是顯著控制因子，解釋力高達(dá)57.14%，明顯高于坡度的單因子作用。此外，植被覆蓋度與其他因子的交互作用大大增強(qiáng)了單因子對土壤侵蝕的解釋力，特別是與坡度因子的交互作用，q值增幅極大，說明了在地形起伏較大的地區(qū)有必要加強(qiáng)植被建設(shè)，通過因子間交互作用抑制土壤侵蝕的發(fā)生。

通過風(fēng)險探測器對不同影響因子各分類的土壤侵蝕平均強(qiáng)度進(jìn)行統(tǒng)計，識別出高風(fēng)險區(qū)，結(jié)果（表9）表明坡度＞35°區(qū)域發(fā)生土壤侵蝕的風(fēng)險最高。河南省土壤侵蝕強(qiáng)度隨著坡度等級的增加而增大，坡度范圍為0～25°時增長緩慢，大于25°時土壤侵蝕強(qiáng)度急劇增加，這部分區(qū)域主要分布在河南省西部的三門峽以及南陽市，土地利用類型以耕地為主，植被覆蓋度較低，極易引發(fā)水土流失。降雨量和地貌類型因子不同分類間平均土壤侵蝕強(qiáng)度無顯著差異，植被覆蓋度的高風(fēng)險區(qū)分布在0～0.3 區(qū)間，且隨著植被覆蓋度的升高，土壤侵蝕強(qiáng)度明顯降低，說明植被對研究區(qū)土壤侵蝕有重要影響。當(dāng)海拔大于1 500 m 時，侵蝕強(qiáng)度比前3 個分類有明顯提升，通過與降雨量圖層的疊置分析，發(fā)現(xiàn)這些區(qū)域也對應(yīng)著強(qiáng)降雨區(qū)，兩因子間協(xié)同作用造成了強(qiáng)烈的侵蝕。

表8 土壤侵蝕影響因子交互作用下的解釋力（q值）Table 8 The q values of dominant interactions between soil erosion influencing factors

表9 各影響因子侵蝕高風(fēng)險區(qū)域及平均侵蝕模數(shù)Table 9 High risk areas of soil erosion and its mean value of erosion modulus

生態(tài)探測器結(jié)果表明，降雨和坡度、降雨和土地利用類型、植被覆蓋度和土地利用類型這三對因子對土壤侵蝕的空間分布影響具有顯著差異，與前面分析相結(jié)合，可以推斷出坡度、土地利用類型、植被覆蓋度和降雨量因子對河南省土壤侵蝕具有主導(dǎo)作用，影響著土壤侵蝕空間分布格局。

總體來說，影響河南省土壤侵蝕的因素多種多樣。河南地形復(fù)雜，雖大部分屬于平原地區(qū)，但北部、西部、南部分別與太行山脈、秦嶺余脈和大別山脈相接，這些區(qū)域多為丘陵和山地，坡度變化幅度大，降雨頻繁且集中，兩者間的協(xié)同作用加速了土壤侵蝕的過程。此外，河南作為農(nóng)業(yè)大省，土地類型多為旱地，林地和草地占比較少，植被可以有效降低雨滴的動能，減少雨水的沖刷作用，避免對土壤的破壞，如果不注重耕作區(qū)的造林防護(hù)，勢必會造成生態(tài)系統(tǒng)的破壞。為了推動黃河流域高質(zhì)量發(fā)展，有必要轉(zhuǎn)變土地利用方式：對于坡耕地，應(yīng)加強(qiáng)退耕還林還草工程建設(shè)，提高植被覆蓋率；對于暴雨易發(fā)地區(qū)，應(yīng)設(shè)立監(jiān)測站進(jìn)行預(yù)警，避免由局部強(qiáng)降雨而引發(fā)的土壤侵蝕惡性循環(huán)；對于西部黃河流經(jīng)的黃土區(qū)，可考慮建設(shè)大型的攔泥蓄水工程，控制泥沙下泄，同時又為耕作區(qū)提供水源保證。

3 結(jié)論

（1）河南省土壤侵蝕以微度侵蝕為主，且呈逐年遞增趨勢，輕度侵蝕區(qū)域面積逐年遞減，中度及以上侵蝕區(qū)域隨年降雨量小幅變化。整體土壤侵蝕狀況呈改善趨勢，其原因是河南省近年來積極實施一系列水土整治措施，取得了良好的成效。

（2）2008—2018 年河南省土壤侵蝕分布主要由輕度侵蝕轉(zhuǎn)變?yōu)槲⒍惹治g，地區(qū)主要集中在西部地區(qū)，山地區(qū)域土壤侵蝕得到抑制，但仍需注意由于局部強(qiáng)降雨而引發(fā)的水土流失。

（3）各因子對土壤侵蝕的解釋力大小依次為坡度＞土地利用類型＞地貌類型＞年平均降雨量＞海拔＞植被覆蓋度。其中坡度的解釋力最大，與其他因子的交互作用也最明顯，是影響土壤侵蝕空間格局分布的主要因素，因此有必要加強(qiáng)地形復(fù)雜地區(qū)的水土治理工作，積極擴(kuò)大森林面積，恢復(fù)生態(tài)平衡，從而實現(xiàn)高質(zhì)量發(fā)展。