梁占岐, 李錦榮, 郭建英, 邢恩德, 劉艷萍, 珊 丹, 侯雪彤
(1.水利部 牧區(qū)水利科學(xué)研究所, 呼和浩特 010020; 2.磴口縣水務(wù)局, 內(nèi)蒙古 磴口 015200)
?
退耕還林工程建設(shè)對(duì)土壤侵蝕的影響
——以安塞縣為例
梁占岐1, 李錦榮1, 郭建英1, 邢恩德1, 劉艷萍1, 珊 丹1, 侯雪彤2
(1.水利部 牧區(qū)水利科學(xué)研究所, 呼和浩特 010020; 2.磴口縣水務(wù)局, 內(nèi)蒙古 磴口 015200)
以黃土丘陵區(qū)退耕還林典型縣——安塞縣為研究對(duì)象,基于3S技術(shù)與RUSLE土壤侵蝕模型,從縣域尺度上分析評(píng)價(jià)了黃土高原區(qū)退耕還林前后土壤侵蝕變化。研究結(jié)果表明:與退耕前1999年相比較,2010年土壤侵蝕強(qiáng)度在空間上發(fā)生明顯改變,總體上,土壤侵蝕強(qiáng)度有明顯減小的趨勢(shì)。極強(qiáng)烈侵蝕下降幅度最大(下降13.73%),主要向強(qiáng)烈轉(zhuǎn)移,占10.45%;中度侵蝕由35.92%增加到59.98%,主要由強(qiáng)烈侵蝕轉(zhuǎn)移而來(lái),占27.08%;微度和輕度變化較小。土壤侵蝕既有增強(qiáng)區(qū)域也有減弱區(qū)域,總體趨勢(shì)以減弱為主。侵蝕加強(qiáng)區(qū)域主要是由中度和強(qiáng)烈向強(qiáng)烈、極強(qiáng)烈和劇烈轉(zhuǎn)移,轉(zhuǎn)移面積較小;發(fā)生增強(qiáng)的區(qū)域主要是以草地覆蓋為主;其次是原耕地,退耕還林后轉(zhuǎn)化為林地和草地,部分地區(qū)出現(xiàn)裸露斑塊,呈現(xiàn)部分小面積土壤侵蝕加劇。土壤侵蝕減弱區(qū)域主要是極強(qiáng)烈、劇烈向中度、強(qiáng)烈等轉(zhuǎn)移;轉(zhuǎn)移面積主要發(fā)生區(qū)域溝壑區(qū)低蓋度草地、溝道以及河道邊緣裸地。隨著退耕還林還草工程的實(shí)施,荒山林草覆蓋度增加,林草植被對(duì)降水進(jìn)行截留下滲,緩減洪峰流量,減少降雨洪峰對(duì)溝道和河道側(cè)沖刷,降低土壤侵蝕強(qiáng)度。安塞縣實(shí)施退耕還林11年后,土壤侵蝕以強(qiáng)度侵蝕為主(46.47%)轉(zhuǎn)中度侵蝕為主(59.98%),全縣平均土壤侵蝕由1998年的9 780 t/(km2·a)轉(zhuǎn)為2010年的5 460 t/(km2·a),每年約減少土壤侵蝕量1 274萬(wàn)t。退耕還林工程建設(shè)對(duì)于安塞縣控制水土流失和改善生態(tài)環(huán)境有著重要作用。研究結(jié)果將對(duì)該區(qū)域水土流失治理提供參考依據(jù)。
退耕還林; 土壤侵蝕; 安塞縣; RUSLE
黃土高原的生態(tài)環(huán)境惡化已成為阻礙當(dāng)?shù)乜沙掷m(xù)發(fā)展的重大問(wèn)題,水土流失、土地沙化問(wèn)題日趨嚴(yán)重。水土流失造成土壤生產(chǎn)力下降、土壤貧瘠、植被覆蓋度降低。特別是該地區(qū)的陡坡耕地,人為加劇水土流失。國(guó)家在1999年提出了退耕還林(草)政策,在黃土高原區(qū)主要是將大于25度的耕地,種植當(dāng)?shù)剜l(xiāng)土樹(shù)種,進(jìn)行生態(tài)修復(fù)。改變土地利用方式,增加地表覆蓋度,降低水土流失。土壤侵蝕的發(fā)生發(fā)展特性與地表植被、土地利用方式息息相關(guān)[1-2]。土地利用方式的改變、植被蓋度的增加必將引起一系列的生態(tài)過(guò)程變化。生態(tài)過(guò)程對(duì)生態(tài)環(huán)境及其小氣候產(chǎn)生積極影響,更進(jìn)一步促進(jìn)植被生長(zhǎng)。所以,有必要對(duì)退耕還林建設(shè)工程實(shí)施若干年后,黃土高原的水土流失情況展開(kāi)研究。
本文以全國(guó)退耕還林典型縣——安塞縣為例,利用TM遙感影像解譯土地利用狀況、植被覆蓋度和DEM數(shù)據(jù),在ArcGIS地理信息系統(tǒng)的支持下,利用ULSE模型,對(duì)黃土高原區(qū)安塞縣退耕還林前后的土地利用、植被覆蓋度、坡耕地面積引起的土壤侵蝕強(qiáng)度的變化、安塞縣自退耕還林生態(tài)建設(shè)工程實(shí)施以來(lái)土壤侵蝕強(qiáng)度變化做定量研究。
1.1研究區(qū)概況
安塞縣位于陜西省延安市以北38 km處(108°51′36″—109°26′6″E,36°30′40″—37°19′25″N),北靠榆林市靖邊縣,東及東南與寶塔區(qū)、子長(zhǎng)縣相連,南接甘泉縣,西依志丹縣,海拔高度為997~1 731 m,全縣總土地面積2 950 km2。地貌類型主要為黃土梁澗、梁峁?fàn)铧S土丘陵和溝谷階地[3],溝壑密度為4.7萬(wàn)條/km2。境內(nèi)水土流失嚴(yán)重,全縣水土流失面積2 832 km2,占全縣總面積的96%[4]。安塞縣地處西北內(nèi)陸黃土高原,屬中溫帶大陸性半干旱季風(fēng)氣候,夏秋多雨,冬季嚴(yán)寒干燥,年平均氣溫為8.8℃,最高值為9.3℃,最低為8℃,年平均日照時(shí)數(shù)為2 395.6 h,無(wú)霜期157 d,降雨主要集中在夏季的7月、8月、9月,平均降水量505.3 mm。土壤以黃綿土為主,土壤發(fā)育層厚25~60 cm,pH值為7~8.5,呈弱堿性或中堿性,保水保肥性能差,有機(jī)質(zhì)含量不足1%,土質(zhì)綿軟,結(jié)構(gòu)疏松,遇水后極易崩解。境內(nèi)植被較差,天然植被破壞殆盡,南部分布有少量天然次生林,森林覆蓋率為17.7%[5-6]。安塞縣1999年開(kāi)始實(shí)施退耕還林還草工程,截止2010年底,安塞縣累計(jì)完成退耕還林77 906.67 hm2,其中退耕地造林41 040 hm2(農(nóng)民人均0.31 hm2),荒山地造林35 600 hm2,封山育林1 266.67 hm2。安塞農(nóng)民人均退耕面積0.31 hm2,是全國(guó)退耕農(nóng)民人均退耕面積的2.7倍。
1.2ULSE模型
在土壤水蝕預(yù)報(bào)方面,ULSE是比較成熟的數(shù)學(xué)模型,它具有建模資料易于獲取、預(yù)報(bào)精度較高、應(yīng)用方便等優(yōu)點(diǎn)。通用土壤流失方程(USLE)由Wischemier等[7]于1958年提出,為提高各因子計(jì)算的通用性,美國(guó)農(nóng)業(yè)部又提出了修正方程RUSLE,模型的數(shù)學(xué)表達(dá)式為:
A=R·K·L·S·C·P
(1)
式中:A——某地塊多年平均降雨徑流年土壤流失量(t/hm2);R——年降雨侵蝕力因子[MJ·mm/(hm2·h)];K——土壤可蝕性因子;L——坡長(zhǎng)因子,無(wú)量綱,反映坡長(zhǎng)對(duì)土壤流失的影響;S——坡度因子,無(wú)量綱,反映坡度對(duì)土壤流失的影響;C——植被因子,無(wú)量綱,反映農(nóng)作物、牧草或樹(shù)木覆蓋對(duì)土壤流失的影響;P——土壤保持措施因子,無(wú)量綱,反映等高耕作等土壤保持措施與順坡耕作農(nóng)地小區(qū)的土壤流失率的比值。
1.3R值計(jì)算
1.4K值計(jì)算
K因子反映了土壤對(duì)侵蝕的敏感性及降水所產(chǎn)生的徑流量與徑流速率的大小。影響K因子的因素很多,但一般來(lái)說(shuō)主要與土壤的粒徑和有機(jī)質(zhì)有關(guān)。依據(jù)前人的研究,K值的大小與土壤質(zhì)地有較高的相關(guān)性。Sharpley等[9]在EPIC棋型中,把土壤可蝕性因子K的計(jì)算公式發(fā)展為:
(2)
式中:Sa——砂粒(0.1~2 mm)含量百分比;Si——粉粒(0.002~0.1 mm)含量百分比;Cl——黏粒(<0.002 mm)含量百分比;C——有機(jī)碳含量百分含量;Sn——粉粒和黏粒含量,Sn=1-Sa/100。利用1∶100萬(wàn)土壤分類圖(數(shù)據(jù)來(lái)源于國(guó)家自然科學(xué)基金委員會(huì)“中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心”http:∥westdc.westgis.ac.cn),提取所需要的信息計(jì)算K值。數(shù)據(jù)庫(kù)中提供的T_SAND砂粒(0.0625~2 mm),黏粒T_CLAY(<0.002 mm)和粉粒T_SILY(0.002~0.062 5 mm)的分類標(biāo)準(zhǔn)不一致,利用3次樣條插值法求得0.1 mm處的值[10],利用公式(2)計(jì)算K值,結(jié)果見(jiàn)附圖8。
1.5LS因子計(jì)算
1.5.1L坡長(zhǎng)因子利用DEM柵格數(shù)據(jù)和ArcGIS,根據(jù)公式(3)提取坡長(zhǎng)λ。利用L和坡長(zhǎng)的關(guān)系求取坡長(zhǎng)因子。
(3)
式中:λ——坡長(zhǎng)(m);m——坡長(zhǎng)效應(yīng)指數(shù),當(dāng)坡度
1.5.2S坡度因子緩坡上選用McCool等[11]研究的坡度公式,在陡坡上采用Liu等[12]的坡度公式,見(jiàn)公式(4)。
(4)
式中:S——坡度因子;θ——坡度。
按照上述方法和公式(4),計(jì)算得到安塞縣坡長(zhǎng)坡度因子LS(附圖9)。
1.6C作物管理及植被蓋度因子的確定
植被具有截留降雨、減緩徑流、保土固土等生態(tài)功能,對(duì)控制土壤侵蝕起著決定性的作用,山地植被覆蓋率的大小直接影響土壤侵蝕程度的強(qiáng)弱。因而,不管是早期的土壤侵蝕人工普查,還是近來(lái)各種理論下建立的土壤侵蝕量預(yù)測(cè)模型,都將植被覆蓋度作為衡量地表水土流失程度的一個(gè)重要因子。
通用土壤流失量方程USLE中,植被因子C與植被類型、植被覆蓋度有關(guān),是根據(jù)地表植物覆蓋狀況不同而反映植被對(duì)土壤流失影響的因素,它是方程諸因子中相對(duì)變化范圍較大的因子。
1.6.1植被蓋度的提取植被蓋度的計(jì)算現(xiàn)在大多數(shù)是利用遙感提取植被指數(shù),利用植被指數(shù)與外業(yè)調(diào)查的植被蓋度建立相關(guān)關(guān)系,最終得到植被蓋度圖。而最常用的就是利用植被NDVI指數(shù),來(lái)估算植被蓋度。
根據(jù)安培浚等[13]在西北干旱區(qū)利用NDVI所建立的植被指數(shù)與植被覆蓋度的關(guān)系來(lái)確定植被覆蓋度c。利用NDVI計(jì)算植被覆蓋度c的公式為:
(5) 式中:NDVI——所求像元的歸一化植被指數(shù);NDVImin,NDVImax——研究區(qū)NDVI的最小值和最大值。
1.6.2C因子的確定C值主要受植被覆蓋度和土地利用現(xiàn)狀的制約,本次研究將參考國(guó)內(nèi)外的研究成果[13-16],結(jié)合安塞縣的實(shí)際情況利用公式(6)計(jì)算得各地類年均C因子。利用ArcGIS中的Raster Calculator,利用下面公式計(jì)算得到C因子(附圖10)。
(6)
1.7P保護(hù)措施因子的確定
P因子為采用專門措施后的土壤流失量與采用順坡種植時(shí)的土壤流失量的比值。國(guó)內(nèi)P值的獲得基本上是根據(jù)區(qū)域特點(diǎn)對(duì)土地的不同利用方式賦值。最終參照劉得俊等[17]的研究成果,結(jié)合我們分類得到的土地利用現(xiàn)狀圖進(jìn)行P因子賦值。坡耕地在實(shí)施退耕還林工程后,農(nóng)業(yè)用地基本上屬于等高耕作和梯田,侵蝕量比較大而且沒(méi)有任何保護(hù)措施的坡耕地基本都種植喬木和灌木,部分轉(zhuǎn)化為草地。所以可以據(jù)此對(duì)不同土地利用下的P因子賦值:退耕前林地、高覆蓋度草地和低覆蓋度草地的均賦值為1,建筑用地和水域的賦值為0,而農(nóng)田則賦值為0.30;退耕后考慮到安塞縣從1999年的退耕還林工程的實(shí)施,>25°的坡耕地全部退耕地為林草地,轉(zhuǎn)化之后的有林地和灌木林地均有侵蝕防治措施,有林地和灌木林地進(jìn)行反坡整地、魚(yú)鱗坑等。所以在2010年的土地利用新增有林地和灌木林地的P值取0.8,未變化有林地和灌木林地取1。退耕后耕地基本是較為平坦的旱地和梯田,P值取0.35,其他土地利用類型P值退耕還林前后無(wú)變化(表1)。
表1 退耕還林前后P值變化
注:退耕的有林地、灌木林地P=0.8以及剩余的地勢(shì)較為平坦的耕地和梯田P=0.35。
2.1土壤侵蝕等級(jí)變化
在得到ULSE通用方程的因子圖后,利用ArcGIS柵格計(jì)算器,把各因子相乘得到土壤流失數(shù)據(jù),按照水利部LS190—2007標(biāo)準(zhǔn)將土壤流失進(jìn)行分級(jí)(表2),得到土壤侵蝕強(qiáng)度。由于其單位為英制,需要進(jìn)行單位轉(zhuǎn)換,乘以系數(shù)224.2,即可轉(zhuǎn)換為t/(km2·a)的公制單位,得到各像元的年土壤流失量。
從表2可以明顯看到,與退耕前1999年相比較,2010年土壤侵蝕強(qiáng)度在空間上發(fā)生明顯改變,總體上,土壤侵蝕強(qiáng)度有明顯減小的趨勢(shì)。與1999年相比較,在實(shí)施退耕還林的11年后,土壤侵蝕的輕度明顯降低,水土流失得到有效改善。劇烈侵蝕面積由1999年的2.75%下降到0.17%,極強(qiáng)烈侵蝕面積由14.55%下降到0.82%,下降了13.73%,下降幅度最大。強(qiáng)烈侵蝕面積由46.74%下降到39.02%。退耕還林工程實(shí)施1年后,侵蝕強(qiáng)度有了顯著降低。從表2只能看到單純面積的變化,想進(jìn)一步對(duì)侵蝕強(qiáng)度動(dòng)態(tài)變化趨勢(shì),分析土壤侵蝕強(qiáng)度是如何轉(zhuǎn)移變化的,還需通過(guò)轉(zhuǎn)移矩陣方法獲得(表3)。
表2 安塞縣土壤侵蝕強(qiáng)度
2.2土壤侵蝕空間轉(zhuǎn)移變化
從表3看,微度、輕度變化不大,大部分還是保持原有侵蝕水平,只有部分侵蝕加劇向更高一級(jí)侵蝕轉(zhuǎn)移;中度侵蝕28.82%保持不變,其他分別向強(qiáng)烈、極強(qiáng)烈和劇烈轉(zhuǎn)移,轉(zhuǎn)移面積分別占6.89%,0.15%,0.05%;從土地利用類型來(lái)看,轉(zhuǎn)移面積主要發(fā)生在耕地、林地、中蓋度草地向中蓋度草地和低蓋度草地轉(zhuǎn)移。
表3 安塞縣土壤侵蝕強(qiáng)度轉(zhuǎn)移矩陣
強(qiáng)烈有27.08%向中度轉(zhuǎn)移,19.32%保持不變,0.25%向極強(qiáng)烈等級(jí)轉(zhuǎn)移,約0.08%向劇烈轉(zhuǎn)移。強(qiáng)烈侵蝕主要發(fā)生在耕地,大部分是坡度較的坡耕地,其中還有部分中低蓋度草地,2010年后19.32%保持不變,主要土地利用類型為草地;27.08%向中度轉(zhuǎn)移的主要是耕地轉(zhuǎn)變?yōu)榱值睾筒莸?,?jīng)過(guò)十多年恢復(fù),植被蓋度有顯著增加,侵蝕等級(jí)呈現(xiàn)下降趨勢(shì)。向極強(qiáng)烈和劇烈轉(zhuǎn)移的面積主要是分布在低蓋度草地。可能是退耕還林還草中個(gè)別地區(qū)在林地和草地中出現(xiàn)裸露地表斑塊。
極強(qiáng)烈大部分向強(qiáng)烈轉(zhuǎn)移,占10.45%;只有0.25%保持不變,有3.8%向中度轉(zhuǎn)移,向微度和輕度轉(zhuǎn)移的共約0.05%,保持侵蝕等級(jí)降低趨勢(shì)。1999年極強(qiáng)烈侵蝕區(qū)主要分布在低蓋度草地和部分裸地,從地理位置上分布溝壑深處,隨著退耕還林還草工程的實(shí)施,封山對(duì)植被恢復(fù)起到積極作用。
劇烈侵蝕區(qū)主要分布在溝道位置以及河流邊緣,其中有0.27%向中度侵蝕轉(zhuǎn)移,2.32%向強(qiáng)烈侵蝕轉(zhuǎn)移,有0.16%向極強(qiáng)烈轉(zhuǎn)移。這主要是由于退耕還林工程實(shí)施,林地和草地對(duì)部分降水進(jìn)行截留下滲,緩減洪峰流量,減少降雨洪峰對(duì)溝道和河道側(cè)沖刷,降低土壤侵蝕強(qiáng)度。
安塞縣2010年的全縣平均土壤侵蝕為5 460 t/(km2·a),屬于強(qiáng)度侵蝕,其研究結(jié)果與張巖等[18]利用中國(guó)土壤侵蝕模型(CSLE)對(duì)吳起縣2009年縣域平均土壤侵蝕強(qiáng)度的研究結(jié)果基本相似,安塞縣和吳起縣同屬于黃土高原區(qū),土壤類型和水土流失現(xiàn)狀較為相似,其水土流失模數(shù)具有參考價(jià)值。
總體來(lái)看,安塞縣退耕還林前以強(qiáng)度侵蝕為主,占安塞縣土地總面積的46.47%,退耕還林后的2010年以中度侵蝕為主,占安塞縣土地總面的59.98%。由此可見(jiàn),安塞縣退耕還林后水土流失得到了有效遏制,生態(tài)建設(shè)成果突出,但土壤侵蝕狀況仍然較為嚴(yán)重,需要繼續(xù)加強(qiáng)鞏固退耕還林工程建設(shè)成果,充分發(fā)揮退耕還林工程建設(shè)的水土保育效應(yīng)。
安塞縣實(shí)施退耕還林后11年來(lái),水土流失控制效果明顯,較高強(qiáng)度土壤侵蝕的土地面積在大幅度減少。實(shí)施退耕還林還草工程11年后,土壤侵蝕強(qiáng)度發(fā)生顯著變化。土壤侵蝕既有增強(qiáng)區(qū)域也有減弱區(qū)域,總體趨勢(shì)以減弱為主。侵蝕加強(qiáng)區(qū)域主要是中度和強(qiáng)烈向強(qiáng)烈、極強(qiáng)烈和劇烈轉(zhuǎn)移,轉(zhuǎn)移面積較??;發(fā)生增強(qiáng)的區(qū)域主要是以草地覆蓋為主;其次是原耕地,退耕還林后轉(zhuǎn)化為林地和草地,部分地區(qū)出現(xiàn)裸露斑塊,呈現(xiàn)部分小面積土壤侵蝕加劇。土壤侵蝕發(fā)生減弱區(qū)域主要是極強(qiáng)烈、劇烈向中度、強(qiáng)烈等轉(zhuǎn)移;轉(zhuǎn)移面積主要發(fā)生區(qū)域溝壑區(qū)低蓋度草地、溝道以及河道邊緣裸地,隨著退耕還林還草工程的實(shí)施,在荒山造林,封山育林使原有荒山林草覆蓋度增加,林草植被對(duì)降水進(jìn)行截留下滲,緩減洪峰流量,減少降雨洪峰對(duì)溝道和河道側(cè)沖刷,降低土壤侵蝕強(qiáng)度。
研究區(qū)退耕還林11年后土壤侵蝕分布以中度侵蝕最多,其面積占安塞縣國(guó)土總面積的59.98%,而退耕還林前安塞縣土壤侵蝕以強(qiáng)度侵蝕為主,占安塞縣國(guó)土總面積的46.47%;安塞縣退耕還林11年后的全縣平均土壤侵蝕為5 460 t/(km2·a),屬于強(qiáng)度侵蝕,其退耕還林前1998年的平均土壤侵蝕為9 780 t/(km2·a),為退耕還林后的1.80倍,每年約減少土壤侵蝕量1 274萬(wàn)t。由此可見(jiàn),退耕還林工程建設(shè)對(duì)于安塞縣控制水土流失和改善生態(tài)環(huán)境有著重要作用。
[1]張魯,周躍,張麗彤.國(guó)內(nèi)外土地利用與土壤侵蝕關(guān)系的研究現(xiàn)狀與展望[J].水土保持研究,2008,15(3):43-48.
[2]郭碧云,王光謙,張正峰.內(nèi)蒙古農(nóng)牧交錯(cuò)區(qū)土地利用/覆被變化及土壤侵蝕:以太仆寺旗為例[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2011,19(5):9-16.
[3]焦峰,溫仲明,石輝,等.黃土高原安塞縣土地結(jié)構(gòu)[J].山地學(xué)報(bào),2004,22(4):406-410.
[4]杜英,楊改河,劉志超等.黃土丘陵區(qū)退耕還林(草)生態(tài)經(jīng)濟(jì)系統(tǒng)的能值分析:以安塞縣為例[J].干旱地區(qū)農(nóng)業(yè)研究,2008,26(5):189-196.
[5]陳源泉,高旺盛.基于農(nóng)業(yè)生態(tài)服務(wù)價(jià)值的農(nóng)業(yè)綠色GDP核算:以安塞縣為例[J].生態(tài)學(xué)報(bào),2007,27(1):250-259.
[6]董孝斌,高旺盛,嚴(yán)茂超.基于能值理論的農(nóng)牧交錯(cuò)帶兩個(gè)典型縣域生態(tài)經(jīng)濟(jì)系統(tǒng)的耦合效應(yīng)分析[J].農(nóng)業(yè)工程學(xué)報(bào),2005,21(11):1-6.
[7]Wischmeier W H, Smith D D. Predicting rainfall erosion losses: A guide to conservation planning [Z]. Washington D C: US Department of Agricultural Science and Education Administration,1978.
[8]孫保平,趙廷寧,齊實(shí).USLE在西吉縣黃土丘陵溝壑區(qū)的應(yīng)用[J].中國(guó)科學(xué)院水利部西北水土保持研究所集刊(黃土高原試驗(yàn)區(qū)土壤侵蝕和綜合治理減沙效益研究專集),1990(2):50-58.
[9]Sharpley A N, Williams J R. EPIC-erosion/productivity impact calculator: 1. Model documentation.[J]. Technical Bulletin-United States Department of Agriculture,1990,4(4):206-207.
[10]蔡永明,張科利,李雙才.不同粒徑制間土壤質(zhì)地資料的轉(zhuǎn)換問(wèn)題研究[J].土壤學(xué)報(bào),2003,40(4):511-517
[11]McCool D K, Foster G R, Mutchler C K, et al. Revised slope length factor for the universal soil loss equation[J]. Transactions of the Asae,1987,30(5):1387-1396.
[12]Liu B Y, Nearing M A, Shi P J, et al. Slope length effects on soil loss for steep slopes[J]. Soil Science Society of America Journal, 2000,64(5):1759-1763.
[13]安培浚,劉樹(shù)林,領(lǐng)耀文,等.植被指數(shù)遙感定量研究:以民勤綠洲為例[J].遙感技術(shù)與應(yīng)用,2005,20(6):576-581.
[14]Wischmeier W H, Smith D D. Predicting rainfall erosion losses: A guide to conservation planning. U. S.[M]. Washington: Science and Education Administration,1987.
[15]Wischmeier W H, Smith D D. Predicting rainfall erosion losses from cropland east of the Rocky Mountains: A guide for selection of practices for soil and water conservation[M]. Washington: Science and Education Administration,1965.
[16]江忠善,王志強(qiáng),劉志.黃土丘陵區(qū)小流域土壤侵蝕空間變化定量研究[J].土壤侵蝕與水土保持學(xué)報(bào),1996,1(2):1-9.
[17]劉得俊,李潤(rùn)杰,王文卿,等.基于地理信息系統(tǒng)的西寧市土壤侵蝕監(jiān)測(cè)的實(shí)現(xiàn)[J].水土保持研究,2006,13(5):111-113.
[18]張巖,劉憲春,李智廣,等.利用侵蝕模型普查黃土高原土壤侵蝕狀況[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(10):165-171.
Effects of Returning Farmland to Forest Construction on Soil Erosion—A Case Study of Ansai County
LIANG Zhanqi1, LI Jinrong1, GUO Jianying1, XING Ende1,LIU Yanping1, SHAN Dan1, HOU Xuetong2
(1.InstituteofWaterResourcesforPastoralArea,MinistryofWaterResources,Hohhot010020,China; 2.WaterworksAuthorityofDengkouCounty,Dengkou,InnerMongolia015200,China)
Selecting the typical loess hilly-gully region returning farmland to forest County as the research example, based on 3S technology and the soil erosion RUSLE model, we analyzed and evaluated loess soil erosion change before and after returning farmland to forest loess plateau zones on the county scale. The results showed that compared with the former farmland in 1999, soil erosion intensity obviously changed on the spacial scale in 2010, on the whole, the soil erosion intensity decreased significantly. Extremely strong erosion decreased most (13.73% decline), shift of very strong erosion mainly to strongly accounted for 10.45%. Moderate erosion increased from 35.92% to 59.98%, which was transferred from the strong erosion, accounting for 27.08%. Slight and mild erosion changes were small. The enhanced area and weakened area of soil erosion were observed, and the overall trend of erosion density was declining. The transfer of moderate erosion to strong, and very strong erosion occurred in strengthened erosion area, but transferred area was small; enhanced area was mainly dominated by grassland coverage; the second one was the farmland converted to forestland and grassland, and bare patch appeared in some parts with the increased intensity of soil erosion. The transfer of strong erosion to moderate, and very strong erosion to strong erosion occurred in weakened erosion areas. The areas was grassland with low coverage, channel, and bare land of channel gully edge. With the implementation of Grain for Green Project, grass coverage increased in the former barren land, rainfall is trapped by vegetation and infiltrates in soil. Plants play the roles in mitigating peak flow and reducing channel erosion,which makes the soil erosion intensity decrease. After the implementation of 11-year Ansai County Forest, soil erosion intensity changed from the strong erosion (46.47%) to moderate erosion (59.98%), the average soil erosion in the county changed form 9 780 t/(km2·a) in 1998 into 5 460 t/(km2·a) in 2010, the annual reducing soil erosion was 12.74 million t. Returning farmland to forest construction played an important role in the control of soil erosion and improvement of ecological environment. These results will provide reference for soil erosion control in this region.
returning farmland to forest; soil erosion; Ansai County; RUSLE
2015-05-14
2015-06-23
國(guó)家林業(yè)公益性行業(yè)科研專項(xiàng)“中國(guó)南方退耕還林工程建設(shè)效益評(píng)價(jià)研究”(201004018)
梁占岐(1962—),男,內(nèi)蒙古鄂爾多斯人,碩士,高級(jí)工程師,主要從事草地水土保持與生態(tài)環(huán)境保護(hù)及草地節(jié)水灌溉研究。E-mail:mkslzq@126.com
李錦榮(1980—),男,內(nèi)蒙古呼和浩特人,博士,高級(jí)工程師,主要從事荒漠化防治及生態(tài)遙感研究。E-mail:lijinrong918@126.com
S157
A
1005-3409(2016)03-0077-05