降雨侵蝕過程中黑土團聚體流失特征

盧　嘉，鄭粉莉,2,*，安　娟，李桂芳

1 西北農(nóng)林科技大學水土保持研究所，黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室，楊凌　712100 2 中國科學院水利部水土保持研究所，楊凌　712100 3 臨沂大學水土保持與環(huán)境保育研究所，臨沂　276005 4 中國科學院大學，北京　100049

?

降雨侵蝕過程中黑土團聚體流失特征

盧嘉1，鄭粉莉1,2,*，安娟3，李桂芳2,4

1 西北農(nóng)林科技大學水土保持研究所，黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室，楊凌712100 2 中國科學院水利部水土保持研究所，楊凌712100 3 臨沂大學水土保持與環(huán)境保育研究所，臨沂276005 4 中國科學院大學，北京100049

摘要:坡面侵蝕過程中土壤團聚體流失可直接反映土壤團聚體的破碎程度以及雨滴打擊和徑流搬運之間的相互作用?；谀M降雨試驗，研究降雨侵蝕過程中黑土坡面土壤團聚體流失特征。試驗處理包括黑土區(qū)常見的2個侵蝕性降雨(50、100mm/h降雨強度)和2個坡度(5°和7.5°)以及2種地面處理(裸露休閑和秸稈覆蓋)。結(jié)果表明：裸露休閑處理下徑流含沙量是秸稈覆蓋處理的27.5—141.3倍，且不同處理下含沙量最大值均出現(xiàn)在降雨初期。覆蓋秸稈與無覆蓋試驗處理下土壤團聚體流失均以<0.25 mm 微團聚體為主。秸稈覆蓋試驗處理下泥沙中<0.25 mm微團聚體流失量占團聚體流失總量的34.5%—56.8%，而在裸露休閑處理下其值達到82%以上。秸稈覆蓋處理下的各粒級團聚體流失量較裸露休閑試驗處理減少了33.3%以上，其中差異最明顯的是≥1 mm粒級團聚體與<0.25 mm微團聚體，二者分別較裸露休閑處理減少了43.1%—96.4%和99.0%以上。秸稈覆蓋處理下的0.25—2 mm粒級團聚體流失比例較裸露休閑處理明顯增加。秸稈覆蓋試驗處理下流失團聚體的平均重量直徑(MWD)和幾何平均直徑(GMD)分別是裸露休閑處理的1.5—2.9和1.7—2.0倍；而秸稈覆蓋處理下的平均重量比表面積(MWSSA)和分形維數(shù)(D)較裸露休閑處理分別減少了26.2%—32.9%和5.1%—6.7%。

關(guān)鍵詞：模擬降雨；降雨侵蝕過程；團聚體流失；黑土區(qū)

坡面土壤侵蝕主要來源于土壤團聚體的破壞，而雨滴打擊和坡面薄層水流為破碎團聚體的流失提供了動力條件[1- 3]。降雨條件下土壤團聚體破碎機制主要包括消散作用(快速濕潤條件下土壤團聚體中空氣壓縮產(chǎn)生的應力而引起的破碎)、機械作用(雨滴打擊和徑流搬運等)、黏粒膨脹作用(土壤礦物濕潤后非均勻膨脹引起的破碎)和土壤物理-化學分散破壞作用[1, 4]。有關(guān)研究結(jié)果表明，降雨條件下土壤團聚體破碎機制主要為消散作用和機械作用[4- 5]。消散作用主要發(fā)生在降雨初期，而降雨過程中土壤團聚體的破碎主要受雨滴打擊作用，即雨滴對大團聚體的拆分和徑流搬運過程對團聚體的破壞作用[6- 8]。陳曉燕等[9]研究表明，降雨強度越大土壤團聚體分散作用越明顯。周一楊等[10]基于濺蝕試驗研究發(fā)現(xiàn)，>1 mm粒級團聚體不易發(fā)生遷移，雨滴擊濺作用主要體現(xiàn)為對>1 mm粒級團聚體進行拆分。郭進[11]等通過模擬降雨發(fā)現(xiàn)降雨初期侵蝕泥沙中微團聚體含量在65%以上，隨降雨進行團聚體的破壞方式由最初的雨滴分離濺散轉(zhuǎn)變?yōu)樗骷羟?、搬運。因此，降雨條件下增加地表覆蓋可降低消散和機械作用對土壤團聚體的破壞，減少土壤團聚體流失。秸稈覆蓋作為水土保持耕作措施之一，不僅可改善土壤結(jié)構(gòu)[12- 13]，抑蒸保墑、培肥地力，同時還可增加土壤入滲和土壤含水量，減少地表徑流、土壤團聚體流失和土壤侵蝕的發(fā)生[14- 15]。

黑土具有良好的團粒結(jié)構(gòu)和較高的有機質(zhì)含量，黑土區(qū)是我國主要的商品糧產(chǎn)區(qū)。但近幾十年來由于過度墾殖和開發(fā)，致使該區(qū)域土壤侵蝕加劇，水土流失日趨劇烈，黑土功能嚴重退化。土壤團聚體作為土壤結(jié)構(gòu)的基本組成單位，對土壤侵蝕有重要的影響作用。因此開展降雨侵蝕過程中黑土團聚體流失特征及其破碎機制的研究十分必要?，F(xiàn)有研究結(jié)果表明，黑土團聚體的破碎機制主要是消散作用和粘粒膨脹作用[16- 17]，而對于降雨侵蝕過程中團聚體破碎機制及其流失特征的研究鮮見報道。據(jù)此，本研究采用模擬降雨試驗，對比分析不同降雨強度和坡度條件下，秸稈覆蓋對黑土區(qū)土壤團聚體流失的影響，探討降雨侵蝕過程中土壤團聚體的破碎機制，以期為黑土寶貴資源的保護提供科學依據(jù)。

1材料與方法

1.1試驗裝置與材料

降雨試驗在黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室人工模擬降雨大廳進行。降雨設(shè)備為側(cè)噴式人工降雨裝置，采用2組單噴頭對噴，降雨高度為16 m，所有雨滴均可達到終點速度，降雨特性與天然降雨特性相似[18]，降雨均勻度大于85%。試驗所用土槽為8.0 m(長)×1.5 m(寬)×0.6 m(深)的坡度可調(diào)式鋼槽，坡度調(diào)節(jié)范圍為0—35°。在試驗土槽下端設(shè)徑流收集裝置，用于采集試驗過程中的徑流泥沙樣。

供試土壤取自吉林省榆樹市劉家鎮(zhèn)(北緯44°43′28″，東經(jīng)126°11′47″)的0—20 cm的耕作層土壤。土壤質(zhì)地為粉壤土，其顆粒組成(美國農(nóng)業(yè)部土壤質(zhì)地分級制)為：砂粒(2000—50 μm)含量3.3%，粉粒(50—2 μm)含量76.4%，粘粒(< 2 μm)含量20.3%；耕作層土壤容重為1.20 g/cm3，有機質(zhì)(重鉻酸鉀氧化-外加熱法)含量為23.81 g/kg，pH值(水浸提法，水土比2.5∶1)為5.92。

1.2試驗設(shè)計

試驗地表處理包括：裸露休閑和秸稈覆蓋。秸稈覆蓋是在土槽表面均勻鋪設(shè)長約5 cm的玉米秸稈，覆蓋密度為1 kg/m2。試驗設(shè)計的降雨強度和坡度根據(jù)當?shù)貙嶋H情況。詹敏等[19]和張憲奎等[20]研究表明，導致黑土區(qū)土壤流失的主要降雨類型為短歷時、高強度降雨，持續(xù)時間多在1 h左右，且該區(qū)中度侵蝕的瞬時雨強標準為I10≥0.71 mm/min。所以結(jié)合此標準設(shè)計2個降雨強度(50、100 mm/h，即0.83、1.67 mm/min)，降雨歷時為60 min。東北黑土區(qū)的坡度為1—8°，很少達到10°[21]。因此本試驗設(shè)計5°和7.5°，2個坡度。每個試驗處理重復2次。

表1　試驗設(shè)計

1.3試驗步驟

為保證模擬降雨過程試驗土槽有良好的透水性，試驗土槽底部均勻打孔，并在底部裝10 cm細沙，沙層之上分為犁底層和耕作層。根據(jù)課題組2007—2013年野外調(diào)查發(fā)現(xiàn)，東北典型黑土區(qū)農(nóng)耕地耕作層深度大多為20 cm左右，黑土層以下為亞粘狀的黃土層。因此根據(jù)試驗土槽的實際情況，設(shè)計20 cm的粘黃土，以模擬犁底層；在其上填充20 cm的黑土，以模擬耕作層。黑土耕作層土壤容重變化范圍為1.15—1.25 g/cm3，犁底層土壤容重變化范圍為1.28—1.48 g/cm3[22]。結(jié)合野外實際調(diào)查資料，本研究中耕層土壤容重設(shè)計為1.20 g/cm3，犁底層土壤容重為1.35 g/cm3。為了保證填土的均勻性，采用分層填土法，每5 cm為一層。填裝上層土之前，打毛土壤表面防止發(fā)生分層現(xiàn)象。試驗土槽準備好后，利用25 mm/h降雨強度對坡面進行預降雨至坡面產(chǎn)流為止。預降雨24 h后進行正式降雨。

每次正式降雨前對降雨強度進行率定，以保證降雨的均勻度達到試驗要求。降雨試驗過程中當坡面產(chǎn)流后，記錄產(chǎn)流時間，同時接取第一個徑流泥沙樣，隨后每隔6 min連續(xù)采集兩次徑流泥沙樣，將接取的徑流泥沙樣先稱重，一部分用于分析坡面侵蝕過程，一部分進行濕篩分析，將徑流泥沙樣依次通過孔徑為5、2、1、0.5、0.25 mm的套篩，按濕篩法測定團聚體粒級分布和流失狀況。降雨結(jié)束后，將侵蝕泥沙樣和篩分后的各粒級團聚體(≥5、2—5、1—2、0.5—1、0.25—0.5、<0.25 mm，6個粒級)風干稱重，編號裝入自封袋中。

1.4土壤團聚體特征指標

采用團聚體平均重量直徑(MWD)、幾何平均直徑(GMD)、團聚體平均重量比表面積(MWSSA)和分形維數(shù)(D)4個指標來評價不同地表條件下降雨侵蝕過程中土壤團聚體流失特征。其計算式為：

(1)

(2)

(3)

式中，xi為土壤粒級的平均直徑(mm)；wi為不同粒級土壤團聚體占總土壤團聚體流失量的比例；ρ1=……ρn=2.65 g/cm3；i為所分級團聚體組別；n為所分級團聚體粒級總組數(shù)。

D計算公式(楊培嶺[23])如下：

(4)

式中，W為直徑小于Ri累計質(zhì)量；Wt為總質(zhì)量；Ri為相鄰兩粒級土壤團聚體平均直徑；Rmax為最大粒級土壤團聚體平均直徑。

2結(jié)果與分析

2.1秸稈覆蓋對坡面侵蝕的影響

2.1.1坡面徑流量和侵蝕量

在50和100 mm/h降雨強度處理下，秸稈覆蓋可明顯減少地表徑流16.1%—37.1%(表2)。這是因為秸稈覆蓋有保護地表的作用，減弱了雨滴動能，同時又有對雨水的吸收和對徑流的分散作用，增大了土壤入滲。坡面侵蝕方面，覆蓋秸稈后坡面侵蝕量較裸露休閑對照處理下減少了97.7%以上?？梢?，秸稈覆蓋能有效的減少徑流量和侵蝕泥沙量，且與削減徑流作用相比，秸稈覆蓋的減沙效應更為明顯。覆蓋秸稈后地表糙率增加，同時削減雨滴打擊，消耗徑流沖刷力，從而削弱了徑流搬運泥沙的能力。Rajan和Khera[24]研究表明，0.6 kg/m2的秸稈覆蓋程度同裸露休閑處理相比徑流量和產(chǎn)沙量分別減少66%和87%。這與本研究中坡面徑流和侵蝕量的減少率有所差異。這是因為兩組試驗土壤質(zhì)地不同，Rajan和Khera所研究的土壤為砂壤土，孔隙度較大，入滲量要大于黑土；其次，本試驗為了保證試驗前期土壤條件基本一致，在正式降雨試驗前進行了前期預降雨，使得前期土壤含水量較高(土壤含水量為(21.11±0.89)%)，導致土壤入滲作用較少，坡面徑流量增加，進而導致秸稈覆蓋下徑流量的減少率較小；而兩組試驗秸稈覆蓋量的不同的是造成坡面侵蝕量減少率有所差異的主要原因。由此可見不同質(zhì)地土壤在覆蓋秸稈后均可減少地表徑流和坡面土壤流失。因此在黑土區(qū)采取秸稈還田能有效保護耕地質(zhì)量。

2.1.2坡面侵蝕過程與團聚體破碎

坡面侵蝕是降雨與土壤界面之間的響應過程，受降雨過程和下墊面變化的影響。在持續(xù)降雨條件下，由于土壤團聚體破碎方式的不同，其侵蝕過程也有所差異。對比秸稈覆蓋和裸露休閑處理下含沙量變化發(fā)現(xiàn)，秸稈覆蓋可以延緩產(chǎn)流時間6—10 min(圖1)。覆蓋秸稈消除雨滴動能后，土壤團聚體免遭雨滴打擊破壞，土壤空隙破壞程度降低，有利于入滲，從而減緩坡面徑流和侵蝕的發(fā)生。坡面產(chǎn)流時間隨降雨強度的增加也相應提前。50 mm/h降雨強度下裸露休閑處理的含沙量表現(xiàn)為降雨初期最大，隨降雨歷時的增加呈快速下降后又趨于穩(wěn)定的變化趨勢。而100 mm/h降雨強度下裸露休閑處理的含沙量隨降雨歷時的增加表現(xiàn)為在從坡面產(chǎn)流到降雨歷時10 min期間呈快速增加的趨勢，而后坡面徑流含沙量隨降雨歷時的增加呈急速下降后又趨于穩(wěn)定的變化趨勢。裸露休閑處理下含沙量均表現(xiàn)為降雨初期出現(xiàn)最大值，這是因為導致黑土團聚體破碎的主要原因為快速濕潤處理下的消散作用[17]。降雨過程的初始階段就是一個模擬土壤快速濕潤的過程，土壤團聚體發(fā)生孔隙內(nèi)部封閉的空氣壓力破碎，地表徑流將破碎團聚體沖走，使徑流含沙量在降雨初期達到最大，可以認為孔隙內(nèi)部封閉的空氣壓力破碎團聚體(即消散作用)是產(chǎn)流前期主要的泥沙來源[25]。降雨初期不同降雨強度下含沙量變化有所差異的原因是，100 mm/h降雨強度下雨滴打擊力較大，可進一步引起表層土壤團聚體的破碎，增大了泥沙量；而50 mm/h降雨強度下的雨滴打擊力相對較弱，徑流以搬運地表松散物質(zhì)為主。隨著降雨歷時的增加，侵蝕物質(zhì)來源減少，含沙量呈下降趨勢；最后坡面徑流量與泥沙來源趨于穩(wěn)定，徑流含沙量也趨于穩(wěn)定。

表2　裸露休閑和秸稈覆蓋試驗處理下坡面徑流量和侵蝕量

*減少率(%)=(秸稈覆蓋-裸露休閑)/裸露休閑×100%; **大寫字母表示同一行不同處理下，在P<0.05水平上差異顯著；小寫字母表示同一列不同處理下，在P<0.05水平上差異顯著

不同降雨強度和坡度下秸稈覆蓋處理下含沙量遠小于裸露休閑處理。秸稈覆蓋處理下的含沙量均小于0.3 g/L(圖1)。裸露休閑處理下的含沙量是秸稈覆蓋處理的27.5—141.3倍。秸稈覆蓋處理下的含沙量均表現(xiàn)為在降雨初期較高，隨著降雨的進行呈下降趨勢，最后呈現(xiàn)平緩趨勢。分析原因是因為除快速濕潤作用引起的土壤團聚體破碎以外，覆蓋秸稈后消除雨滴動能，從而削弱了由雨滴打擊土壤團聚體所造成的團聚體破碎，導致侵蝕泥沙來源的減少，含沙量呈下降趨勢?？梢娊涤赀^程中，除氣爆所引起的消散作用，機械破壞作用對土壤團聚體破碎的影響也較大。進一步說明雨滴打擊是致使黑土團聚體分散的主要動力之一，土壤團聚體的破碎為坡面侵蝕提供物質(zhì)來源。

圖1　裸露休閑與秸稈覆蓋試驗處理的含沙量隨降雨歷時的變化Fig.1　Sediment concentration versus rainfall duration in bare land treatment and straw mulch cover treatment

2.2坡面土壤團聚體流失2.2.1侵蝕泥沙團聚體粒級分布

土壤團聚體可劃分為大團聚體(≥0.25 mm)和微團聚體(<0.25 mm)[26]。土壤團聚體破碎機制的不同，必然導致土壤團聚體流失的差異。表3描述了裸露休閑和秸稈覆蓋處理下侵蝕泥沙中各粒級團聚體流失量，圖2為流失團聚體各粒級的分布狀況。土壤水穩(wěn)性團聚體數(shù)量和組成對土壤結(jié)構(gòu)的穩(wěn)定性具有重要作用，影響土壤通透性、抗蝕性，是土壤肥力的主要指標之一。從圖2可以看出，供試黑土≥0.25 mm水穩(wěn)性團聚體含量達53.8%，可見供試土壤具有較好的團聚結(jié)構(gòu)。不同處理條件下各粒級土壤團聚體流失量均表現(xiàn)為秸稈覆蓋<裸露休閑(圖2)。不同處理下侵蝕泥沙中≥0.25 mm (≥2 mm、1—2、0.5—1 和0.25—0.5 mm)各粒級團聚體流失量與<0.25 mm粒級團聚體之間存在顯著性差異，侵蝕泥沙中團聚體的流失以<0.25 mm微團聚體為主。秸稈覆蓋處理下，泥沙中<0.25 mm微團聚體流失量占團聚體流失總量的34.5%—56.8%，裸露休閑處理下達到82%以上。不同處理下侵蝕泥沙團聚體的流失以<0.25 mm微團聚體為主，這與申艷等[27]和黃滿湘等[28]研究結(jié)果一致，說明侵蝕過程中，細顆粒泥沙容易被選擇搬運，且搬運距離較遠。對于本研究設(shè)計的5°坡面和50 mm/h降雨強度試驗處理，在有秸稈覆蓋條件下，流失團聚體的粒級均為<0.25 mm微團聚體，這是因為這個處理下雨強和坡度較小，加之覆蓋秸稈后，雨滴對土壤表面的打擊作用減弱，同時削弱了坡面徑流的搬運能力。此處理下坡面徑流量侵蝕量均較小，侵蝕量僅為0.9 g/m2。秸稈覆蓋處理下各個粒級團聚體流失量較裸露休閑下減少了33.3%以上，其中差異最明顯的是≥1 mm粒級團聚體與<0.25 mm微團聚體。秸稈覆蓋處理下≥1 mm粒級團聚體較裸露休閑處理下的減少率達到43.1%—96.4%，<0.25 mm微團聚體流失量減少率更是達到99%以上。裸露休閑與秸稈覆蓋之間各粒級的流失比例及流失量的差異表明，秸稈覆蓋在消除雨滴打擊對團聚體的破壞作用的同時，也削弱了徑流搬運破壞團聚體的能力?？梢?，裸露休閑處理下黑土團聚體的破碎主要來自于雨滴打擊對團聚體的破壞，雨水濕潤和徑流攜帶搬運過程也可引起其團聚體的破碎。而覆蓋秸稈后消除了雨滴打擊對土壤團聚體的破壞作用，坡面徑流減小，團聚體的破碎主要來自于雨水的濕潤作用。

表3　裸露休閑和秸稈覆蓋試驗處理下不同粒級土壤團聚體的流失量

*裸露休閑The bare land；秸稈覆蓋：The straw mulch cover; **同一行不同小寫字母表示不同粒徑團聚體在P<0.05水平上差異顯著

圖2　裸露休閑和秸稈覆蓋試驗處理下流失團聚體各粒級流失量占流失總量比例Fig.2　Soil aggregate size distribution in bare land and straw mulch cover treatment

不同試驗處理下侵蝕泥沙中流失團聚體粒級分布存在差異(圖2)。除5°坡面和50 mm/h降雨強度試驗處理外，其它幾種試驗處理下覆蓋秸稈0.25—2 mm (1—2、0.5—1 和0.25—0.5 mm)各粒級團聚體流失比例與裸露休閑處理相比明顯增加，分別增加了0.8—28.7、2.4—28.2、4.3—27.7倍。與供試土壤相比，侵蝕泥沙中≥0.25 mm粒級團聚體的含量則均有所減小，其中裸露休閑處理下≥0.25 mm粒級團聚體含量較供試土壤中減少了67.8%—94.4%。這是因為在降雨侵蝕過程中，大團聚體在雨滴打擊和雨水濕潤的破壞作用下分散成小團聚體，團聚體穩(wěn)定性下降(表4)，從而使得團聚體組成結(jié)構(gòu)也變?yōu)橐孕×綀F聚體(<0.25 mm 微團聚體)為主。而秸稈覆蓋消除了雨滴對地表團聚體的打擊作用，使得大團聚破碎程度降低，由大團聚體破碎所供給的微團聚體流失量減少，從而導致泥沙中大團聚體流失比例增大，也說明雨滴打擊對大團聚體的破壞作用很大。王春燕等[29]研究表明，紅壤侵蝕程度越高，>0.25 mm水穩(wěn)性團聚體的含量越低，團聚體的分散度越大，團聚體越容易遭到破壞，穩(wěn)定性減小。Bernard和Eric[30]認為>0.25 mm水穩(wěn)性大團聚體含量與土壤侵蝕程度和土壤團聚體穩(wěn)定性二者之間均呈負相關(guān)。本研究中>1 mm粒級團聚體流失比例低于17%?？梢姶髨F聚體不易發(fā)生遷移，雨滴打擊主要對其進行拆分[10]。在100 mm/h降雨強度下，≥0.25 mm各粒級團聚體流失比例均隨著坡度的增加而增大，增加幅度為0.7—6.8倍。說明坡度的增加使得坡面徑流侵蝕力增強，從而有利于大團聚體的流失。侵蝕泥沙中團聚體粒級的分布受到雨滴打擊、徑流搬運和團聚體穩(wěn)定性共同作用。大粒級團聚體由微團聚體或土壤單粒通過膠結(jié)物質(zhì)膠結(jié)而成，內(nèi)聚力相對較小[31]，在較強雨滴擊濺分散下易破碎為較小粒級團聚體，該結(jié)果通過對比供試土壤與各不同處理下的團聚體粒級分布及穩(wěn)定性特征可以證實(圖2和表4)；而<0.25 mm微團聚體的形成主要取決于粘粒和三二氧化物的粘結(jié)作用，微團聚體孔隙小，彎曲程度大，其內(nèi)聚力也較大，在侵蝕過程中不易破碎，團聚體較為穩(wěn)定[32]，這與閆峰陵[33]等研究結(jié)果相同。綜上可見，降雨擊濺分散和徑流搬運能力對黑土區(qū)土壤不同粒級團聚體流失均有所影響。

2.2.2侵蝕泥沙中微團聚體流失過程

以上分析了侵蝕過程中團聚體流失特征，發(fā)現(xiàn)降雨侵蝕過程中黑土團聚體的流失以<0.25 mm微團聚體為主。因此，這里著重分析這部分微團聚體在降雨過程中的流失狀況(圖3)。

裸露休閑處理下<0.25 mm微團聚體流失量明顯大于秸稈覆蓋處理下的。裸露休閑處理下坡面<0.25 mm微團聚體流失量隨降雨過程呈先快速增加后緩慢上升趨勢，其中降雨初期100 mm/h降雨強度下微團聚體流失量的增漲趨勢明顯大于50 mm/h降雨強度下。秸稈覆蓋處理下坡面<0.25 mm微團聚體流失變化過程可分為3個階段。第1階段從坡面產(chǎn)流到產(chǎn)流后5 min內(nèi)，此時坡面微團聚體流失量隨降雨歷時的增加呈快速增加的趨勢，這與裸露休閑處理的第一階段趨勢一致。這一方面主要是因為降雨初期徑流搬運能力較弱，因此優(yōu)先遷移微團聚體；另一方面造成裸露休閑處理下微團聚體流失量遠大于秸稈覆蓋處理下的原因為，雨滴打擊引起團聚體的破碎一般發(fā)生在降雨初期的幾分鐘，而覆蓋秸稈后消除了雨滴對地表團聚體的打擊作用。第2階段秸稈覆蓋處理下微團聚體流失量急劇下降，而裸露休閑處理下此階段微團聚體流失量隨降雨歷時的增加呈緩慢增加的趨勢。這是因為裸露休閑處理下侵蝕物質(zhì)的來源不僅來自于松散物質(zhì)，還有來自破碎的大團聚體；在降雨過程中，隨著雨滴對團聚體連續(xù)的打擊作用，機械外力累積能量逐漸克服團聚體中膠結(jié)物質(zhì)膠結(jié)力的影響，大團聚體逐漸破碎為較小的團聚體，為微團聚體的流失提供了物質(zhì)來源，從而增加了坡面微團聚體的流失量，導致裸露處理下微團聚體流失量隨降雨歷時的增加而呈現(xiàn)增加的趨勢，表現(xiàn)出雨滴對大團聚體的拆分能力。但是覆蓋秸稈后，消除了雨滴打擊力，減小了雨滴對大團聚體的破碎供給微團聚體，導致微團聚體流失量減少。第3階段為降雨后期，秸稈覆蓋處理下微團聚體流失量呈緩慢下降趨勢，這是因為雨滴打擊對大團聚體破碎作用較弱，加上侵蝕物質(zhì)來源的減少，徑流對土壤搬運能力也較弱，因此微團聚體流失量減小。

圖3　裸露休閑和秸稈覆蓋試驗處理下微團聚體流失量隨降雨歷時變化Fig.3　Micro-aggregate loss versus rainfall duration in bare land treatment and straw mulch cover treatment

2.3表征土壤團聚體流失指標的分析

土壤團聚體平均重量直徑(MWD)和幾何平均直徑(GMD)常作為評價土壤團聚體狀況的指標，其值越大表示土壤的團聚度越高，團聚體穩(wěn)定性就越強。土壤團聚體分形維數(shù)(D )反映了土壤水穩(wěn)定性團聚體及水穩(wěn)性大團聚體含量對土壤結(jié)構(gòu)與穩(wěn)定性的影響趨勢，其值越大，則團聚體的分散度越大[34]。本文選取平均重量直徑(MWD)、幾何平均直徑(GMD)、團聚體平均重量比表面積(MWSSA)和分形維數(shù)(D)4個指標進一步分析裸露休閑處理和秸稈覆蓋處理下黑土團聚體流失特征。

表4　裸露休閑和秸稈覆蓋試驗處理下流失團聚體特征指標

*裸露休閑The bare land；秸稈覆蓋The straw mulch cover；供試土壤The tested soil; **MWD：平均重量直徑 Mean Weight Diameter；GMD：幾何平均直徑 Geometric Mean Diameter；MWSSA：平均重量比表面積 Mean Weight Soil Specific；D：分形維數(shù) Fraction Dimension

從表4中看出，除5°坡面和50 mm/h降雨強度試驗處理外，各團聚體特征指標在裸露休閑和秸稈覆蓋處理下表現(xiàn)出顯著性差異。秸稈覆蓋處理下流失團聚體的MWD和GMD表現(xiàn)為最大。秸稈覆蓋處理下流失團聚體的MWD和GMD分別是裸露休閑處理下的1.5—4.3和1.7—2.0倍。造成這種現(xiàn)象的主要原因是覆蓋秸稈后雨滴打擊拆分土壤團聚體的作用減弱，雨滴打擊和徑流搬運破壞的能量不能破碎這些穩(wěn)定的團聚體，從而使大粒徑團聚體含量明顯增加，泥沙中大粒級團聚體含量較高，從而導致流失團聚體的MWD和GMD有所增加。同時MWD和GMD兩個指標隨著坡度的增加而增大，進一步說明坡度有利于大團聚體的流失。而秸稈覆蓋處理下MWSSA和D較裸露休閑處理下分別減少了26.2%—32.9%和5.1%—6.9%。流失團聚體的MWSSA和D值越大，說明侵蝕泥沙中小粒級團聚體含量增多，團聚體破壞性變大，相應的穩(wěn)定性變差。Castrignano等[34]發(fā)現(xiàn)，當 D 值接近2時，團聚體主要由數(shù)量很少的大團聚體組成，而隨著 D值的增大，土壤中小粒級團聚體的數(shù)量隨之增加。吳承禎等[35]研究指出，>0.25 mm粒級團聚體含量越低，其團聚體粒徑分布的分形維數(shù)越高。本研究秸稈覆蓋處理下≥0.25 mm粒級團聚體流失比例遠大于裸露休閑處理下(圖2)，這也是秸稈覆蓋處理下D值小于裸露休閑處理下的原因。可見即使在團聚體含量較高的黑土，雨滴打擊對土壤團聚體的分散作用仍較為明顯，秸稈覆蓋可有效防止土壤團聚體的破碎和流失。

3結(jié)論

本研究基于模擬降雨試驗和試驗現(xiàn)場篩分泥沙中的團聚體，研究了降雨侵蝕過程中的東北黑土區(qū)坡面土壤團聚體流失特征和破碎機制，得到如下研究結(jié)論：

(1)秸稈覆蓋試驗處理下的坡面徑流量和侵蝕量分別較裸露休閑試驗處理減少了16.3%—37.1%和97.7%—99.5%?？梢娕c削減徑流作用相比，秸稈覆蓋的減沙效應更為明顯。

(2)不同降雨強度和坡度下裸露休閑處理的徑流含沙量是秸稈覆蓋處理的27.5—141.3倍。含沙量最大值出現(xiàn)在降雨初期。說明降雨過程中黑土團聚體破碎機制主要為消散作用和機械作用，且降雨初期以消散作用為主。

(3)秸稈覆蓋試驗處理各粒級團聚體流失量較裸露休閑試驗處理減少了33.3%以上，其中差異最明顯的是≥1 mm粒級團聚體與<0.25 mm微團聚體，分別較裸露休閑處理下減少了43.1%—96.4%和99.0%以上。覆蓋與無覆蓋秸稈試驗處理下土壤團聚體流失均以<0.25 mm 微團聚體為主。覆蓋秸稈后0.25—2 mm各粒級團聚體流失比例較裸露休閑條件下明顯增加?？梢?，秸稈覆蓋消除了雨滴對坡面團聚體的打擊作用和徑流搬運破壞團聚體的能力。秸稈覆蓋處理下土壤團聚體的破碎方式以消散作用和雨水濕潤作用為主。

(4)秸稈覆蓋處理下，<0.25 mm微團聚體流失量隨降雨歷時的增加呈快速增加—急劇遞減—緩慢下降的趨勢;而裸露休閑處理下，其變化表現(xiàn)為快速增加—緩慢上升趨勢，表明雨滴打擊能夠影響土壤團聚體的破碎。

(5)秸稈覆蓋試驗處理下流失團聚體的MWD和GMD最大，分別是裸露休閑處理下的1.5—4.3和1.7—2.0倍，而MWSSA和D與裸露休閑試驗處理相比分別減少了26.2%—32.9%和5.1%—6.9%。說明秸稈覆蓋可有效防治土壤團聚體的破碎和流失。

參考文獻(References):

[1]Le Bissonnais Y. Aggregate stability and assessment of soil crustability and erodibility: I. Theory and methodology. European Journal of Soil Science, 1996, 47(4): 425- 437.

[2]Zheng F L, Merrill S D, Huang C H, Tanaka D L, Darboux F, Liebig M A, Halvorson A D. Runoff, soil erosion, and erodibility of conservation reserve program land under crop and hay production. Soil Science Society of America Journal, 2004, 68(4): 133- 139.

[3]Roth C H, Eggert T. Mechanisms of aggregate breakdown involved in surface sealing, runoff generation and sediment concentration on loess soils. Soil and Tillage Research, 1994, 32(2/3): 253- 268.

[4]Ramos M C, Nacci S, Pla I. Effect of raindrop impact and its relationship with aggregate stability to different disaggregation forces. Catena, 2003, 53(4): 365- 376.

[5]史志華, 閆峰陵, 李朝霞, 王天巍, 蔡崇法. 紅壤表土團聚體破碎方式對坡面產(chǎn)流過程的影響. 自然科學進展, 2007, 17(2): 217- 224.

[6]Legout C, Leguedois S, Le Bissonnais Y, Malam I O. Splash distance and size distributions for various soils. Geoderma, 2005, 124(3/4): 279- 292.

[7]Young R A, Wiersma J L. The role of rainfall impact in soil detachment and transport. Water Resources Research, 1973, 9(6): 1629- 1636.

[8]Kinnell P I A. The mechanics of raindrop induced flow transport. Australian Journal of Soil Research, 1990, 28(4): 497- 516.

[9]陳曉燕, 牛青霞, 周繼, 魏朝富, 謝德體, 何丙輝. 人工模擬降雨條件下紫色土陡坡地土壤顆粒分布空間變異特征. 水土保持學報, 2010, 24(5): 163- 168.

[10]周一楊, 王恩姮, 陳祥偉. 模擬降雨條件下黑土濺蝕與團聚體分選特征. 水土保持學報, 2008, 22(6): 176- 179.

[11]郭進, 文安邦, 嚴冬春, 史忠林, 巨莉. 三峽庫區(qū)紫色土坡地土壤顆粒流失特征. 水土保持學報, 2012, 26(3): 18- 21.

[12]張仁陟, 黃高寶, 蔡立群, 羅珠珠, 李玲玲, 謝軍紅. 幾種保護性耕作措施在黃土高原旱作農(nóng)田的實踐. 中國生態(tài)農(nóng)業(yè)學報, 2013, 21(1): 61- 69.

[13]劉婷, 賈志寬, 張睿, 鄭甲成, 任世春, 楊寶平, 聶俊峰, 劉艷紅, 王海霞. 秸稈覆蓋對旱地土壤水分及冬小麥水分利用效率的影響. 西北農(nóng)林科技大學學報: 自然科學版, 2010, 38(7): 69- 76.

[14]Myers J L, Wagger M G. Runoff and sediment loss from three tillage systems under simulated rainfall. Soil and Tillage Research, 1996, 39(1/2): 115- 129.

[15]張少良, 張興義, 劉曉冰, 劉爽, 于同艷. 典型黑土侵蝕區(qū)不同耕作措施的水土保持功效研究. 水土保持學報, 2009, 23(3): 11- 15.

[16]王彬. 東北典型薄層黑土區(qū)土壤可蝕性關(guān)鍵因子分析與土壤可蝕性計算[D]. 楊凌: 西北農(nóng)林科技大學, 2009.

[17]張孝存, 鄭粉莉. 基于Le Bissonnais法的東北黑土區(qū)土壤團聚體穩(wěn)定性研究. 陜西師范大學學報: 自然科學版, 2009, 37(5): 82- 86.

[18]周佩華, 張學棟, 唐克麗. 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室土壤侵蝕模擬實驗大廳降雨裝置. 水土保持通報, 2000, 20(4): 27- 30.

[19]詹敏, 厲占才, 信玉林. 黑土侵蝕區(qū)降雨參數(shù)與土壤流失關(guān)系. 黑龍江水專學報, 1998, (1): 40- 43.

[20]張憲奎, 許靖華, 盧秀琴, 鄧育江, 高德武. 黑龍江省土壤流失方程的研究. 水土保持通報, 1992, 12(4): 9- 18.

[21]張曉平, 梁愛珍, 申艷, 李文鳳, 張學林, 王玉璽, 解運杰, 劉鳳飛, 楊學明. 東北黑土水土流失特點. 地理科學, 2006, 26(6): 687- 692.

[22]水利部, 中國科學院, 中國工程院. 中國水土流失防治與生態(tài)安全: 東北黑土區(qū)卷. 北京: 科學出版社, 2010: 64- 66.

[23]楊培嶺, 羅遠培, 石元春. 用粒徑的重量分布表征的土壤分形特征. 科學通報, 1993, 38(20): 1896- 1899.

[24]Rajan B, Khera K L. Effect of tillage and mode of straw mulch application on soil erosion in the submontaneous tract of Punjab, India. Soil and Tillage Research, 2006, 88(1/2): 107- 115.

[25]陳正發(fā), 史東梅, 謝均強, 張兵. 紫色土旱坡地土壤團聚體穩(wěn)定性特征對侵蝕過程的影響. 中國農(nóng)業(yè)科學, 2011, 44(13): 2721- 2729.

[26]Edwards A P, Bremner J M. Micro-aggregates in soils. Journal of Soil Science, 1967, 18(1): 64- 73.

[27]申艷, 張曉平, 梁愛珍, 李文鳳, 楊學明. 黑土坡耕地土壤流失形態(tài)分析——以一次降雨為例. 干旱地區(qū)農(nóng)業(yè)研究, 2008, 26(6): 224- 229.

[28]黃滿湘, 章申, 晏維金. 農(nóng)田暴雨徑流侵蝕泥沙對氮磷的富集機理. 土壤學報, 2003, 40(2): 306- 310.

[29]王春燕, 黃麗, 譚文峰, 胡紅青. 幾種侵蝕紅壤中有機質(zhì)和團聚體的關(guān)系. 水土保持學報, 2007, 21(3): 52- 56.

[30]Bernard B, Eric R. Aggregate stability as an indicator of soil susceptibility to runoff and erosion: validation at several levels. Catena, 2002, 47(2): 133- 149.

[31]Six J, Bossuyt H, Degryze S, Denef K. A history of research on the link between (micro) aggregates, soil biota, and soil organic matter dynamics. Soil Tillage Research, 2004, 79(1): 7- 31.

[32]彭新華, 張斌, 趙其國. 紅壤侵蝕裸地植被恢復及土壤有機碳對團聚體穩(wěn)定性的影響. 生態(tài)學報, 2003, 23(10): 2176- 2183.

[33]閆峰陵, 史志華, 蔡崇法, 李朝霞. 紅壤表土團聚體穩(wěn)定性對坡面侵蝕的影響. 土壤學報, 2007, 44(4): 577- 583.

[34]Castrignano A, Stelluti M. Fractal geometry and geostatistics for describing the field variability of soil aggregation. Journal of Agriculture Engineering Research, 1999, 73(1): 13- 18.

[35]吳承禎, 洪偉. 不同經(jīng)營模式土壤團粒結(jié)構(gòu)的分形特征研究. 土壤學報, 1999, 36(2): 162- 166.

An experimental study of Mollisol aggregate loss characteristics during rainfall erosion processes

LU Jia1, ZHENG Fenli1,2,*, AN Juan3, LI Guifang2,4

1InstituteofSoilandWaterConservation,StateKeyLaboratoryofSoilErosionandDrylandFarmingonLoessPlateau,NorthwestA&FUniversity,Yangling712100,China2InstituteofSoilandWaterConservation,ChineseAcademyofSciencesandMinistryofWaterResources,Yangling712100,China3InstituteofWaterandSoilConservationandEnvironmentalProtection,LinyiUniversity,Linyi276005,China4UniversityofChineseAcademyofSciences,Beijing100049,China

Abstract：Soil aggregate loss during hillslope erosion processes directly reflects the interaction among the degree of soil aggregate breakdown, raindrop impact and runoff transportation. Different soil aggregate breakdown mechanisms result in the size distribution of soil aggregate is different. Furthermore, soil aggregate loss has an important influence on hillslope soil loss. However, existing studies have mostly focused on aggregate stability. Rainfall simulation experiments were conducted to investigate Mollisol aggregate loss during rainfall erosion processes. The experimental treatments included two rainfall intensities (50 mm/h and 100 mm/h), representative of erosive rainfalls in the black soil region of northeast China; two common slope gradients (5° and 7.5°); and two surface conditions (bare land and straw mulch cover). Each experimental treatment had two replications. The tested soil was the Mollisol, which was collected from the upper 20 cm of the plow layer in a maize field in Yushu City, Jilin Province, which is typical of the black soil of the region of northeast China. This study was conducted in the rainfall simulation laboratory of the State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Yangling City, China. A rainfall simulator system with a side-sprinkler was used to apply rainfall. The soil pan was 8 m long, 1.5 m wide and 0.6 m deep. During the rainfall simulations, runoff samples were collected in 15-L buckets. Half of each runoff sample was immediately processed through a set of sieves with apertures of 5 mm, 2 mm, 1 mm, 0.5 mm, and 0.25 mm. After sieving, the lost aggregate samples and the remaining half of each runoff samples were oven-dried to calculate the soil and aggregate losses. The results showed that sediment concentration in the bare land treatments was 27.5—141.3 times greater than that in the straw mulch cover treatments. The maximum values for sediment concentration were observed at the beginning of the rainfall simulation, in all treatments. This indicated that the dominant Mollisol aggregate breakdown mechanisms were slaking and mechanical breakdown, and that the effect of slaking was mainly exerted in the initial stages of rainfall. The <0.25 mm micro-aggregates were the main aggregate size fraction lost in both the bare land and the straw mulch cover treatments. The <0.25 mm micro-aggregate loss comprised 34.5—56.8% of the total aggregate loss in the straw mulch cover treatments, whereas the loss was >82% in the bare land treatments. Compared with the bare land treatments, the loss of every size aggregate decreased >33.3% in the straw mulch cover treatments; thus, the most significant differences in soil aggregate loss were observed in the ≥1 mm and <0.25 mm size fractions between the bare land and the straw mulch cover treatments. Compared with the bare land treatments, the loss of the ≥1 mm size fraction and the <0.25 mm micro-aggregates decreased 43.1%—96.4% and 99.0%, respectively, in the straw mulch cover treatments. The proportion of 0.25—2 mm aggregates in the sediment from the straw mulch cover treatments increased and was greater than that of the bare land treatments. These results indicated that raindrop impact was the main driving force for aggregate breakdown. The straw mulch cover eliminated the effect of raindrop impact on soil aggregate breakdown, and also decreased runoff transport capacity. The mean weight diameter (MWD) and geometric mean diameter (GMD) of aggregates in the sediments for the straw mulch cover treatments were 1.5—2.9 and 1.7—2.0 times greater than those for the bare land treatments. Compared with the bare land treatments, the mean weight soil specific (MWSSA) and fraction dimension (D) decreased 26.2%—32.9% and 5.1%—6.7%, respectively. The above four indicators reflect the aggregate loss characteristics of Mollisols.

Key Words:rainfall simulation; rainfall soil erosion process; aggregates loss; black soil region

基金項目:國家自然科學基金資助項目(41271299)

收稿日期:2014- 10- 28; 網(wǎng)絡(luò)出版日期:2015- 08- 21

*通訊作者

Corresponding author.E-mail: flzh@ms.iswc.ac.cn

DOI:10.5846/stxb201410282108

盧嘉，鄭粉莉，安娟，李桂芳.降雨侵蝕過程中黑土團聚體流失特征.生態(tài)學報,2016,36(8):2264- 2273.

Lu J, Zheng F L, An J, Li G F.An experimental study of Mollisol aggregate loss characteristics during rainfall erosion processes.Acta Ecologica Sinica,2016,36(8):2264- 2273.