周子淵,王 立,王小平,馬海龍,李夢茹
(1.甘肅農(nóng)業(yè)大學(xué) 林學(xué)院,甘肅 蘭州730070;2.甘肅省定西市水土保持科學(xué)研究所,甘肅 定西743000)
水土資源是人類生存最基本的條件,水土流失直接影響到水、土資源的開發(fā)、利用和保護問題。黃土高原西部黃土丘陵溝壑區(qū)屬旱作農(nóng)業(yè)區(qū),降雨分布不均且水土流失嚴重。保蓄土壤水分,減少地表徑流,提高降水利用率,是發(fā)展該區(qū)旱作農(nóng)業(yè)的關(guān)鍵[1]。由定西安家溝流域黃土丘陵溝壑區(qū)的土壤理化性質(zhì)及其氣候條件決定,水力侵蝕是造成該區(qū)域水土流失的主要侵蝕方式,而降雨侵蝕是該區(qū)域水力侵蝕的主要發(fā)生形式。降雨對土壤侵蝕的影響已有較多研究,包括降雨侵蝕力、侵蝕性降雨標準以及降雨對不同地區(qū)土壤侵蝕的影響等[2-7]。魏天興[8]對黃土殘塬溝壑區(qū)降雨侵蝕的研究表明,場降雨量、降雨強度與小流域產(chǎn)沙模數(shù)呈線性相關(guān),防止暴雨侵蝕是防治黃土殘塬溝壑區(qū)水土流失的關(guān)鍵。本文通過對天然降雨在不同土地利用條件下產(chǎn)流產(chǎn)沙特征及其規(guī)律的研究,以期探索出能夠提高坡耕地土壤保水能力、減少水土流失的土地利用方式,為推動甘肅黃土高原西部旱作農(nóng)業(yè)的可持續(xù)發(fā)展和進一步研究區(qū)域水土保持綜合治理提供科學(xué)的理論依據(jù)。
試驗區(qū)域處于黃土高原西部黃土丘陵溝壑區(qū)的定西市水土保持科學(xué)研究所實驗基地,海拔約2 000m,年平均氣溫6.3℃,≥0℃積溫2 933.5℃,≥10℃積溫2 239.1℃,極端最高溫度34.3℃,最低溫度-27.1℃,空氣相對濕度 65.8%,年日照時數(shù)2 408.6h,無霜期141d。多年平均降水量427mm,年蒸發(fā)量1 510mm,為典型的雨養(yǎng)農(nóng)業(yè)區(qū)。土壤類型為典型的黃綿土,土質(zhì)綿軟,土層深厚,質(zhì)地均勻,儲水性能良好;0—20cm土壤容重平均為1.17g/cm3,凋萎含水率7.3%,飽和含水率21.9%。
本研究共設(shè)兩組試驗,分別在10°,15°的坡地上進行,每組試驗選擇4個不同處理,試驗小區(qū)進行隨機區(qū)組排列。試驗小區(qū)均為標準徑流小區(qū),各小區(qū)配備自記雨量計,記錄每次降雨的降雨量及其過程。小區(qū)下部布設(shè)V型集流槽連接徑流收集池,池內(nèi)放置容量為13.5L的集流桶,用以收集徑流和泥沙。如表1所示,兩組試驗按照生長的植物種類各設(shè)4個標準徑流小區(qū),分別為:對照小區(qū)荒地(W)、春小麥(SW)、灌木林沙棘(H)、喬木林油松(P)。8個標準徑流小區(qū)布設(shè)在同一小氣候區(qū)域內(nèi),其中荒地、春小麥的面積為5m×10m,灌木林沙棘、喬木林油松的面積為10m×10m。
表1 各試驗小區(qū)基本立地條件
降雨量及其過程采用自記雨量計自動統(tǒng)計。不同處理的降雨產(chǎn)沙量及其產(chǎn)沙過程與徑流過程同步,人工測定,在開始產(chǎn)流后,每5min測量1次徑流,在集流桶中用體積法求得渾水總量(ml),同時采集混合水樣500ml,將徑流樣靜置24h,除去上層清水,余下的用濾紙濾得泥沙,烘干8h稱重,測定徑流樣泥沙含量。試驗準備階段盡量減少對地表破壞,保持地表原狀。
2.2.1 徑流量及侵蝕量測定 徑流產(chǎn)生后,每5min在集流桶中用體積法求得渾水總量(ml),同時采集混合水樣500ml,過濾后烘干稱重,計算水樣的泥沙含量,進一步計算侵蝕量。侵蝕量=
2.2.2 統(tǒng)計分析 采用Excel 2010和SPSS 19.0統(tǒng)計分析軟件對試驗數(shù)據(jù)進行統(tǒng)計分析。
在多次降雨中選取1次典型天然降雨進行分析(2011年8月20日),選取標準按照國家氣象部門規(guī)定的大雨為降雨強度12h內(nèi)降雨量達到15~29.9mm,或24h內(nèi)降雨量達到25~49.9mm的降雨,及試驗地降雨土壤侵蝕發(fā)生程度來確定。該次降雨其降雨量為28mm,降雨歷時52,30min平均降雨強度為32.3mm/h,屬于典型的大雨并在試驗地造成明顯的降雨土壤侵蝕,研究本次降雨土壤侵蝕過程具有較好的代表性。不同處理的地表徑流過程如圖1所示,無壤中流產(chǎn)生。
如圖1所示,在降雨初期5min,降雨強度不足10mm/h,各處理都沒有產(chǎn)流。持續(xù)5min的低強度降雨,已逐步濕潤土層,土壤表層的入滲率減小。而在降雨強度突然增至16.4mm/h時,來不及入滲的雨水形成地表徑流。坡度較大時不同處理的產(chǎn)流時間均較早,不同處理產(chǎn)流時間經(jīng)LSD法檢驗差異顯著(p<0.05),15°荒地產(chǎn)流時間最早,為6.2min,10°灌木林沙棘產(chǎn)流時間最晚,為8.93min。相同植被覆蓋下坡度為10°處理的產(chǎn)流時間較15°處理晚,10°荒地的產(chǎn)流時間比15°荒地晚1.43min,10°春小麥的產(chǎn)流時間比15°春小麥晚2.31min,10°灌木林沙棘的產(chǎn)流時間比15°灌木林沙棘晚2.08min,10°喬木林油松的產(chǎn)流時間比15°喬木林油松晚1.25min。不同坡度下10°處理的產(chǎn)流時間早晚順序為:10°灌木林沙棘(10°H)>10°春小麥(10°SW)>10°喬木林油松(10°P)>10°荒地(10°W),15°處理產(chǎn)流時間早晚順序為:15°灌木林沙棘(15°H)>15°喬木林油松(15°P)>15°春小麥(15°SW)>15°荒地(15°W)。
圖1 2011年8月20日不同處理的地表徑流過程
不同處理徑流量差異顯著,15°春小麥徑流量最大,為23.22m3/hm2;10°灌木林沙棘徑流量最小,為4.07m3/hm2。不同處理的徑流過程都經(jīng)歷了產(chǎn)流、峰值、穩(wěn)定、消減4個表征時刻。相同植物覆蓋下10°處理的徑流量較15°處理的小。10°處理在產(chǎn)流25min后徑流量達到峰值,15°處理在20min后徑流量達到峰值。這是由于地表面高低不平等微地形的影響,坡地會滯蓄部分水量,但隨著坡度的增大,坡面變陡,地面的這種滯蓄能力下降,部分被積蓄在洼地的水量隨著坡度的增加匯流速度加快,匯集形成徑流。流域坡度越大,地面的積水、滯水作用越弱,徑流量越多[9]。不同坡度下10°處理的徑流量排序為:10°春小麥>10°喬木林油松>10°荒地>10°灌木林沙棘,15°處理的徑流量排序為:15°春小麥>15°荒地>15°灌木林沙棘>15°喬木林油松。
由圖2可以看出,不同處理的侵蝕量差異顯著,15°春小麥侵蝕量最大,為1 132.29kg/hm2,10°灌木林沙棘侵蝕量最小,為19.51kg/hm2。在次降雨條件下不同處理的土壤侵蝕過程都經(jīng)歷了發(fā)生、峰值、消減3個表征時刻。相同植被覆蓋下10°處理的侵蝕量較15°處理的小。10°處理在產(chǎn)流20min后侵蝕量達到峰值,15°處理在產(chǎn)流15min后侵蝕量達到峰值。侵蝕量峰值出現(xiàn)的時間和徑流量峰值出現(xiàn)的時間并不一致,這是因為在降雨初期,地表含水量極低,雨滴降落到地面發(fā)生擊濺侵蝕,濺起表層干燥土壤,雨滴的濺蝕作用明顯,所以在徑流峰值出現(xiàn)之前侵蝕量達到了峰值。隨著降雨時間的推移,地表形成結(jié)皮,侵蝕量出現(xiàn)消減[10]。盡管降雨量、降雨侵蝕力對泥沙流失量產(chǎn)生影響,但泥沙最終通過徑流流失,所以徑流對泥沙流失量影響最大[11]。當10°,15°處理徑流量達到峰值時,10°,15°處理的侵蝕量達到侵蝕量消減階段的最大值。不同坡度下10°處理的侵蝕量排序為:10°春小麥>10°喬木林油松>10°荒地>10°灌木林沙棘,15°處理的侵蝕量排序為:15°春小麥>15°荒地>15°灌木林沙棘>15°喬木林油松。
圖2 2011年8月20日不同處理的產(chǎn)沙過程
由圖1和圖2可以看出,該次降雨不同處理其徑流總量、侵蝕總量排序均為:15°春小麥>10°春小麥>10°喬木林油松>15°荒地>15°灌木林沙棘>10°荒地>15°喬木林油松>10°灌木林沙棘。15°處理的徑流總量、侵蝕總量分別是10°處理的1.11,1.36倍。10°荒地與15°荒地相比,徑流總量、侵蝕總量分別減少了36.1%,17.54%;10°春小麥與15°春小麥相比,徑流總量、侵蝕總量分別減少了16.35%,12.25%;10°灌木林沙棘與15°灌木林沙棘相比,徑流總量、侵蝕總量分別減少了39.31%,77.04%;10°喬木林油松與15°喬木林油松相比,徑流總量、侵蝕總量分別增加了67.63%,85.96%。
層次聚類分析法是目前國內(nèi)外使用最多的一種方法。這種方法的基本思想是:先將n個樣品各自看成一類,然后規(guī)定樣品之間的距離和類與類之間的距離[12]。對2011年8月20日該次降雨各處理的徑流量、侵蝕量進行層次聚類分析,取類間距離d=5時,不同處理的水土流失狀況分為4個層次,10°荒地、10°喬木林油松、15°荒地為1個層次,10°灌木林沙棘、15°喬木林油松為1個層次,10°春小麥、15°春小麥為1個層次,15°灌木林沙棘為1個層次。當d=10時,不同處理的水土流失狀況分為3個層次,10°荒地、10°喬木林油松、15°荒地、10°灌木林沙棘、15°喬木林油松為1個層次,10°春小麥、15°春小麥為1個層次,15°灌木林沙棘為1個層次。由此可得,在10°的坡度下,灌木林沙棘的水土保持效應(yīng)優(yōu)于其他3種植物措施,在15°的坡度下,喬木林油松的水土保持效應(yīng)最佳。
2011年試驗區(qū)年降雨量為372.1mm,侵蝕性降雨量為121.6mm,占年降雨量的32.68%(表2)。侵蝕性降雨集中在6—9月出現(xiàn),近乎造成了全年土壤侵蝕的發(fā)生。10°荒地處理的年徑流模數(shù)、年侵蝕總量分別比15°荒地處理減少了36.13%,17.53%,10°春小麥處理的年徑流模數(shù)、年侵蝕總量分別比15°春小麥處理減少了15.64%,21.54%,10°灌木林沙棘處理的年徑流模數(shù)、年侵蝕總量分別比15°灌木林沙棘處理減少了60.73%,85.95%,10°喬木林油松處理的年徑流模數(shù)、年侵蝕總量分別比15°喬木林油松處理增加了67.65%,76.49%。不同處理年侵蝕模數(shù)排序為:15°春小麥>10°春小麥>10°喬木林油松>15°灌木林沙棘>15°荒地>10°荒地>15°喬木林油松>10°灌木林沙棘。這與次降雨條件下各處理徑流總量、侵蝕總量排序基本一致。
表2 不同處理年降雨徑流及土壤侵蝕特征
(1)在次降雨條件下不同處理徑流過程分為產(chǎn)流、峰值、穩(wěn)定、消減4個過程,土壤侵蝕分為發(fā)生、峰值、消減3個過程。
(2)在10°,15°的坡度下,灌木林沙棘、喬木林油松分別具有明顯的水土保持效應(yīng)。次降雨條件下,不同處理的徑流總量、侵蝕總量排序均為:15°春小麥>10°春小麥>10°喬木林油松>15°荒地>15°灌木林沙棘>10°荒地>15°喬木林油松>10°灌木林沙棘。
(3)對該次降雨各處理的徑流量、侵蝕量進行層次聚類分析,取類間距離d=5時,不同處理的水土流失狀況分為4個層次;當d=10時,不同處理的水土流失狀況分為3個層次,其中在10°的坡度下,灌木林沙棘的水土保持效應(yīng)優(yōu)于其他3種植物措施,在15°的坡度下,喬木林油松的水土保持效應(yīng)最佳。
(4)通過對2011年年降雨徑流及土壤侵蝕特征的分析,得出不同處理年侵蝕模數(shù)排序為:15°春小麥>10°春小麥>10°喬木林油松>15°灌木林沙棘>15°荒地>10°荒地>15°喬木林油松>10°灌木林沙棘。這與次降雨條件下各處理徑流總量、侵蝕總量排序基本一致。再次驗證了在10°,15°的坡度下,H,P分別具有明顯的水土保持效應(yīng)。所以結(jié)合試驗地的氣候、環(huán)境特點,10°坡面適宜種植灌木林沙棘,15°坡面適宜退耕還林。
[1] 張曉艷,王立,黃高寶,等.道地藥材保護性耕作對坡耕地土壤侵蝕的影響[J].水土保持學(xué)報,2008,22(2):58-61.
[2] 章文波,謝云,劉寶元.用雨量和雨強計算次降雨侵蝕力[J].地理研究,2002,21(3):384-390.
[3] Liu Baoyuan,Nearing M A.Practical thresholds for separating erosive and non-erosive stroms[J].Transaction of the ASAE,2002,45(6):1843-1847.
[4] 王萬忠.黃土地區(qū)降雨特性與土壤流失關(guān)系的研究(Ⅲ):關(guān)于侵蝕性降雨的標準問題[J].水土保持通報,1984,3(2):58-63.
[5] 程慶杏,呂萬民,吳百林.土壤侵蝕的雨量標準研究初報[J].中國水土保持科學(xué),2004,2(3):90-92.
[6] 江東,盧喜平,蔣光毅,等.降雨因素對紫色土坡地土壤侵蝕影響的試驗研究[J].西南大學(xué)學(xué)報,2009,31(1):140-144.
[7] 左長清,馬良.天然降雨對紅壤坡地侵蝕的影響[J].水土保持學(xué)報,2005,19(2):1-4.
[8] 魏天興.黃土殘塬溝壑區(qū)降雨侵蝕分析[J].水土保持學(xué)報,2001,15(4):47-50.
[9] 范世香,韓紹文.地面坡度對地表徑流影響的實驗研究[J].水土保持通報,1991,11(4):7-10.
[10] 王生鑫,王立,黃高寶,等.糧草豆隔帶種植保護性耕作對坡耕地土壤水蝕的影響[J].水土保持學(xué)報,2010,24(4):40-43.
[11] 袁東海,王兆騫,陳欣,等.不同農(nóng)作措施紅壤坡耕地水土流失特征研究[J].水土保持學(xué)報,2001,15(4):66-69.
[12] 鮑艷,胡振琪,柏玉,等.主成分聚類分析在土地利用生態(tài)安全評價中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報,2006,22(8):87-90.