曹丹妮, 邸涵悅, 郭忠錄, 馬美景2,
(1.湖北方源東力電力科學(xué)研究有限公司, 武漢 430062;2.武漢市水務(wù)科學(xué)研究院, 武漢 430014; 3.華中農(nóng)業(yè)大學(xué) 水土保持研究中心, 武漢 430070)
土壤優(yōu)先流是指土壤水分快速穿過大孔隙系統(tǒng),并繞過土壤基質(zhì)沿著優(yōu)先途徑在土體剖面內(nèi)運動的現(xiàn)象[1-2]。優(yōu)先流是土壤中普遍存在的一種水分運動,它與地表地下水質(zhì)、土壤養(yǎng)分吸收息息相關(guān)[3],土壤中優(yōu)先流的發(fā)生既可以增加入滲、減少地表徑流,但也會引起土壤養(yǎng)分流失、地下水污染等[4-5],因此開展土壤優(yōu)先流特征研究對揭示土壤入滲機制及水土保持工作具有重要意義[6]。
目前,國內(nèi)外學(xué)者采用不同方法(染色示蹤法、樹脂填充法、張力滲透儀法、離子示蹤法等)分析了土地利用方式、土壤初始含水量、根系特性、植被覆蓋、耕作方式、降雨強度等因素對土壤孔隙特征及土壤優(yōu)先流形態(tài)發(fā)育特征的影響[7-9]。但是這些研究主要集中在農(nóng)地和林地內(nèi)且研究對象較為單一,針對不同演替下紫色土區(qū)土壤性質(zhì)和根系性質(zhì)等如何影響土壤優(yōu)先流形態(tài)發(fā)育特征研究較少。土壤質(zhì)地、土壤中孔隙的連通性和曲折性、土壤有機質(zhì)含量等因素對土壤中優(yōu)先流的運動和發(fā)展有重要的影響。Clark等[10]通對美國東部露天煤礦近地表區(qū)土壤優(yōu)先流特征研究發(fā)現(xiàn),重新造林的區(qū)域的土壤剖面染色面積和入滲速率明顯高于草地;Jiang等[11]通過對不同耕作方式下土壤優(yōu)先流特征研究發(fā)現(xiàn),耕作區(qū)與未耕作區(qū)的土壤優(yōu)先流特征間有明顯差異,且不同耕作方式下土壤中優(yōu)先流特征也不同;陳曉冰等[12]對重慶市四面山針闊混交林、竹林和草地進行染色示蹤試驗發(fā)現(xiàn),3種植被條件下土壤中砂粒含量增加、粉粒和黏粒含量的減少在一定程度上均有助于土壤中優(yōu)先流的發(fā)生;呂文星[13]通過對三峽庫區(qū)坡面荒地、玉米地、柑橘地進行亮藍(lán)染色示蹤試驗發(fā)現(xiàn),土壤總孔隙度、毛管孔隙度和土壤有機質(zhì)含量越高,越有助于土壤中水分的快速運移。根系是土壤層的重要組成部分,草本植物、灌木、喬木等具有不同的根系系統(tǒng),而這些植被的根系系統(tǒng)組合會構(gòu)成較為復(fù)雜的空間網(wǎng)絡(luò)結(jié)構(gòu),從而形成復(fù)雜的根系通道,加快優(yōu)先流過程[3,7]。以胡家山小流域內(nèi)荒地、草地、灌木、針葉林、喬—灌混合林為研究樣地,運用原位染色示蹤技術(shù),對不同演替階段的土壤優(yōu)先流形態(tài)發(fā)育特征及其影響因素進行對比分析,為該區(qū)域不同演替階段土壤入滲機制提供一定的理論基礎(chǔ)和實踐參考。
試驗地位于湖北省丹江口市胡家山小流域,地理坐標(biāo)為111°12′22.0″—111°15′20.5″E,32°44′17.8″—32°49′15.6″N,多年平均氣溫16.1℃,多年平均降雨量797.6 mm,屬于北亞熱帶半濕潤季風(fēng)氣候。流域內(nèi)土壤類型主要為紫色土、石灰土和黃棕壤。
選取研究區(qū)內(nèi)荒地、草地、灌木、針葉林、喬—灌混合林為研究樣地,進行染色示蹤試驗。研究點土壤類型均為紫色土。草地植被覆蓋度約為96%;灌木群落植被覆蓋度約為60%;針葉林和喬—灌混合林的郁閉度約為88%,92%,地表枯枝落葉層厚度為3~4 cm。樣地詳細(xì)情況見表1。
表1 樣地參數(shù)
選取相對平坦的荒地、草地、灌木、針葉林、喬—灌混合林為試驗樣地。試驗前,先去除試驗樣地內(nèi)的枯枝落葉、樹枝和礫石等,并平整土壤表面。利用噴壺將10 L濃度為4 g/L亮藍(lán)溶液均勻噴灑在邊長為50 cm的正方形樣地中,整個試驗過程中無地表徑流產(chǎn)生,噴灑完畢后再用薄膜和綠色植被覆蓋在樣方表面上,防止蒸發(fā)和靜置過程中降雨等外部因素的影響。染色24 h后,開挖3個垂直剖面(第一個剖面距離邊框15 cm、第2個剖面距離邊框25 cm、第3個剖面距離邊框35 cm),總計60個剖面,并用數(shù)碼相機沿坡面垂直方向拍攝圖像。首先將拍攝的剖面染色照片經(jīng)Photoshop 2017軟件處理,之后將其導(dǎo)入Matlab編寫的程序中計算并提取土壤染色剖面各特征參數(shù)(未染色區(qū)域顏色為白色,染色區(qū)域調(diào)為黑色)[6]。
分4層、每層10 cm采集土樣,測定土壤機械組成、質(zhì)地、有機質(zhì)等基本理化性質(zhì)。其中,容重、孔隙度采用環(huán)刀法測定,土壤機械組成采用沉降法測定[14],土壤有機質(zhì)含量采用重鉻酸鉀外加熱法測定。
分4層、每層10 cm采用體積為500 cm3的環(huán)刀采集根系樣品(3個重復(fù)),測定土層中根系含量。將樣品置于0.25 mm土樣篩上,用流水沖去根系上的土樣,利用LA-S系列植物圖像分析系統(tǒng)獲取根長密度,方法為:將根系無重疊的放入根盤(根系呈分散狀態(tài)),倒入干凈自來水至根系完全呈淹沒狀態(tài),點擊EPSON Scan圖標(biāo),開始掃描根系并獲取根長密度。
5種演替梯度下,土壤優(yōu)先流垂直分布形態(tài)特征見圖1,5種演替梯度下的土壤中均產(chǎn)生較為明顯的優(yōu)先流,隨演替的進行和土層深度的增加,土壤中基質(zhì)流發(fā)育程度越不明顯,但優(yōu)先流的數(shù)量和運移路徑均呈增大趨勢。
荒地和草地群落中,0—5 cm的表層土壤均被亮藍(lán)溶液大面積染色,且水平方向上染色分布均勻,說明該深度范圍內(nèi)優(yōu)先流發(fā)育不明顯,土壤水分以基質(zhì)流形式運移為主?;牡貎?nèi),土壤剖面中20 cm以下染色區(qū)域較小,開挖土壤剖面發(fā)現(xiàn),20 cm左右存在粘滯層,土壤的滲透性較低,土壤中優(yōu)先路徑的連通性與延伸性較弱[6]。草地內(nèi),剖面染色區(qū)域位于0—30 cm的土壤中,且呈聚集分布,優(yōu)先路徑的連通性較強,一方面可能與土壤斥水性有關(guān),因為土壤斥水性會阻滯或延遲水分入滲,從而迫使水流借助大孔隙、裂隙、根孔等優(yōu)先通道運動[15];另一方面與草本植物的根系特性有關(guān),與群落4相比,群落1、群落2和群落3的染色深度較大,原因是群落1、群落2和群落3的優(yōu)勢種分別是白茅、小蓬草和牛筋草,群落4的優(yōu)勢種是藎草,與藎草相比、白茅根系和牛筋草根系較粗,小蓬草為直根系草本植物,平均根系直徑也較大,形成較為明顯的根系通道,加速水分運移速率。灌木中,群落1內(nèi)土壤基質(zhì)流發(fā)生層次較深且水平方向上染色分布均勻,主要分布在0—10 cm深度范圍內(nèi),優(yōu)先流運移較淺,主要分布在10—25 cm深度范圍內(nèi),且土壤剖面中僅有1條明顯的優(yōu)先流路徑;其余3個群落的土壤剖面中均存在2條染色路徑較深的優(yōu)先流路徑,其中,群落4的2條染色路徑最深,其鋒部可達(dá)36 cm。針葉林內(nèi),0—3 cm的表層土壤表層染色較為均勻,剖面中亮藍(lán)溶液繞過基質(zhì)區(qū)后呈網(wǎng)狀擴散式蔓延下滲,形成明顯的優(yōu)先路徑。喬—灌混合林內(nèi)基質(zhì)流發(fā)生層次也比較淺,土壤剖面3 cm以下土壤水分分化較嚴(yán)重,土壤剖面染色區(qū)域呈不規(guī)則狀,優(yōu)先流路徑大部分由短小的分支組成。林地土壤染色最大深度可達(dá)50 cm,遠(yuǎn)大于其余3種演替梯度下土壤染色最大深度,一是因為與其余演替階段相比,林地土壤中根系直徑較大且分布較深,導(dǎo)致土壤垂直連通性較好,優(yōu)先流現(xiàn)象更明顯;二是因為林地群落地表具有3~4 cm的枯枝落葉層,枯枝落葉層內(nèi)部與根系之間形成適合厭氧動物與微生物生存的環(huán)境,同時也會加速土壤中死根系的腐爛,根系腐爛后產(chǎn)生孔隙的幾率較高,進一步影響土壤水分入滲[16-17]。
注:Ⅰ為荒地,Ⅱ為草地,Ⅲ為灌木,Ⅳ為針葉林,Ⅴ為喬—灌混合林,下圖同。
5種演替梯度下,土壤剖面染色面積比見圖2,土壤染色面積比在同一演替階段不同深度處以及同一深度處不同演替階段間均表現(xiàn)出明顯的異質(zhì)性。5種演替梯度下針葉林土壤染色面積比最大、其次是喬—灌混合林、灌木、草地,荒地的土壤染色面積比最小。
圖2 土壤剖面染色面積比變化
荒地和草地群落中,0—5 cm的表層土壤均被染色,土壤染色面積比可達(dá)80%以上,隨土層深度增加,染色面積比均減小,但減小趨勢較均勻。草地內(nèi),10 cm以下,群落1的土壤染色面積比較大,群落4的土壤染色面積比較小,原因是群落1的優(yōu)勢種是白茅和小蓬草,群落4的優(yōu)勢種是藎草和小薹草,可能是因為藎草和小薹草均為須根草本植物、其平均根系直徑較小,而白茅根系較粗且延伸深度較大,小蓬草為直根系草本植物、其平均根系直徑較大,根系直徑增大一定程度上會導(dǎo)致根系與土壤接觸面之間的縫隙增大,進而提高土壤中水分運移的速率。灌木、喬—灌混合林內(nèi),隨土層深度增加,染色面積比也呈減小趨勢,但并非單調(diào)遞減,其中,灌木中群3和群落4在16,12 cm左右、喬—灌混合林中群落2和群落4在15,25 cm左右均出現(xiàn)局部反彈現(xiàn)象,說明在此處優(yōu)先流均發(fā)生側(cè)滲現(xiàn)象。灌木中群落2的染色面積比隨土層深度增加很少出現(xiàn)較大反彈上升的情況,說明以白刺花和小薹草為優(yōu)勢種的群落的土壤中發(fā)生側(cè)向流的幾率較小。針葉林內(nèi),由于存在水分的側(cè)向入滲,隨土層深度增加,土壤染色面積比整體上呈“S”型趨勢減小,先在0—10 cm土層中迅速減小,在10—20 cm土層中略有增大趨勢,20 cm以下又逐漸減小,土壤染色最大深度可達(dá)50 cm;群落3在10—30 cm染色較為均勻,染色面積比為36%~42%,說明在此深度范圍內(nèi),土壤中水分沿著穩(wěn)定的滲流通道向深層土壤運移。
土壤染色路徑數(shù)量的空間變化可以反映出優(yōu)先路徑的垂直連通狀況和分支性[6,9]。染色路徑寬度(Stained Path Width,SPW)是指在某一指定深度處,被染色的優(yōu)先流區(qū)域所占的寬度。演替梯度下土壤染色路徑數(shù)分布特征見圖3,各個群落內(nèi),土壤剖面的染色路徑總數(shù)和SPW<20 mm的數(shù)量隨土層深度的增加其變化規(guī)律較一致,均隨土層深度增加表現(xiàn)出多峰狀態(tài),但峰值隨土層深度增加呈減小趨勢;SPW>200 mm的數(shù)量隨土層深度增加呈減小趨勢。5種演替梯度下,每個剖面土壤染色路徑總數(shù)出現(xiàn)峰值的位置均不同,但都集中在0—20 cm的表層土壤中,說明水分運動在土壤表層分化較嚴(yán)重。土壤水分運移以優(yōu)先流為主時,SPW<20 mm的數(shù)量越多,土壤水分運移以基質(zhì)流為主時,SPW>200 mm的數(shù)量越多[18],由圖3可知,隨演替的進行,染色路徑總數(shù)、SPW<20 mm的數(shù)量和SPW>200 mm的數(shù)量波動趨勢逐漸增大。土壤中染色路徑總數(shù)和SPW<20 mm的數(shù)量隨演替的進行呈增大趨勢,而SPW>200 mm的數(shù)量隨演替的進行呈減小趨勢,一定程度上說明隨演替的進行,土體優(yōu)先流現(xiàn)象越明顯。
圖3 土壤剖面染色路徑數(shù)分布特征
荒地內(nèi),土壤中染色路徑總數(shù)為0~33個,在2 cm左右,染色路徑總數(shù)值最小,說明此深度處,水分運動以基質(zhì)流為主,無明顯優(yōu)先流;與其余演替梯度相比,該演替梯度下SPW<20 mm的數(shù)量較少,但在10 cm左右,SPW<20 mm的數(shù)量最大,說明該位置水分分化較嚴(yán)重。草地內(nèi),土壤染色路徑數(shù)量最大值出現(xiàn)在7—15 cm的土層中,群落1的染色路徑總數(shù)和SPW<20 mm的最大值出現(xiàn)在15 cm左右,其深度略大于其余3個群落,與該群落內(nèi)優(yōu)勢種白茅的根系特性有關(guān);在10 cm左右土層中,群落4中染色路徑總數(shù)值最大,可達(dá)54個,10 cm以下,染色路徑總數(shù)迅速減小,說明此處土壤優(yōu)先流分化最嚴(yán)重。灌木中,群落4內(nèi)染色路徑數(shù)隨土層深度增加變動幅度相對較小,說明該群落的土壤優(yōu)先路徑在垂直方向上的連通性較好。林地內(nèi),土壤染色路徑數(shù)量隨土層深度增加變化波動較頻繁,說明林地內(nèi)土壤優(yōu)先流路徑隨土層深度變化極不穩(wěn)定。針葉林內(nèi),在20 cm左右,群落2,群落3和群落4內(nèi)染色路徑總數(shù)和SPW<20 mm的數(shù)量均增大,說明此處優(yōu)先流分化較嚴(yán)重。喬—灌混合林中,土層深度在25 cm左右時,群落4的染色路徑數(shù)量值較大,開挖土壤剖面發(fā)現(xiàn),此深度處,剖面內(nèi)存在一個直徑約4 cm的蛇洞,可能是因為蛇洞的存在會使本來緊實的土壤變的松動,促進了土壤水分和溶質(zhì)的運移。
土壤理化性質(zhì)與土壤染色面積比之間相關(guān)性見表2,土壤容重、毛管孔隙度、土壤有機質(zhì)含量與土壤染色面積比之間的相關(guān)性較強?;牡睾筒莸刂?,土壤容重與土壤染色面積比之間存在極顯著負(fù)相關(guān)性(p<0.01),相關(guān)值分別-0.67,-0.71,是因為土壤容重越大,土層越緊實,土壤中大孔隙數(shù)量越少,土壤中水分入滲的優(yōu)先流路徑減少[6]?;牡?、草地和灌木群落中,土壤有機質(zhì)含量與土壤染色面積比之間存在極顯著正相關(guān)性(p<0.01),相關(guān)值分別0.68,0.64,0.63;喬—灌混合林中,土壤有機質(zhì)含量與土壤染色面積比之間也存在顯著正相關(guān)性(p<0.05),相關(guān)值為0.57,是因為有機質(zhì)的存在會促進土壤團粒結(jié)構(gòu)的形成,團粒和團粒之間的非毛管孔隙越多,越有利于水分入滲[19]。
表2 土壤理化性質(zhì)與土壤染色面積比之間相關(guān)性
土壤質(zhì)地對土壤中優(yōu)先路徑的形成與發(fā)展有一定的影響。鄭欣等[20]通過對北京地區(qū)2種類型土壤(淋溶褐土和潮土)研究發(fā)現(xiàn),2種類型土壤的砂粒含量與土壤染色面積比之間存在極顯著正相關(guān)關(guān)系,粉粒含量和黏粒含量與土壤染色面積比之間存在極顯著負(fù)相關(guān)關(guān)系(p<0.01)。孫龍等[21]對不同林齡柑橘林地土壤特性與優(yōu)先路徑分布特征進行相關(guān)性分析發(fā)現(xiàn),砂粒含量、粉粒含量和黏粒含量并沒有對優(yōu)先路徑構(gòu)成顯著影響。研究發(fā)現(xiàn),灌木和針葉林中,粉粒含量與土壤染色面積比之間存在顯著負(fù)相關(guān)關(guān)系(p<0.05);5種演替梯度下,砂粒含量、黏粒含量與土壤染色面積比之間無顯著相關(guān)關(guān)系(p>0.05),說明5種演替梯度下,砂粒含量和黏粒含量對土壤中優(yōu)先流染色面積影響不顯著,這可能與各演替梯度下土壤基質(zhì)差異不大有關(guān)。
根長密度是指單位土壤體積所含根系的總長度,一定程度上影響著土壤中水分的運移和傳輸。將5個演替梯度下土壤剖面中的植物根系按照直徑D≤0.5 mm,0.5 圖4 演替梯度下根長密度分布特征 根長密度與土壤染色面積比之間相關(guān)性見表3,根長密度與土壤染色面積比之間總體上呈正相關(guān)關(guān)系。土壤染色面積比與荒地、草地、灌木和喬—灌混合林的根長密度間具有極顯著相關(guān)性(p<0.01),相關(guān)值分別為0.80,0.66,0.6,0.77,而與針葉林相關(guān)性系數(shù)為0.43,說明植物根系極大地影響著荒地、草地、灌木和喬—灌混合林草地優(yōu)先流的發(fā)生發(fā)展,對針葉林的影響較草地小。 表3 根長密度與土壤染色面積比之間相關(guān)性 不同根系徑級的根長密度對土壤優(yōu)先流有著不同的影響。草地和灌木內(nèi),土壤染色面積比與根系直徑為D≤0.5 mm和0.5 王趙男等[22]通過對長白山系榛子灌木林根系進行染色示蹤試驗發(fā)現(xiàn),D<1 mm的根系對土體優(yōu)先流的貢獻度最大。研究發(fā)現(xiàn),演替初期,D≤2 mm的細(xì)根與土壤染色面積比之間的相關(guān)性較大;演替種后期,2 (1) 隨演替的進行,土壤中優(yōu)先流數(shù)量增多且運移路徑增大;5種演替梯度下針葉林土壤染色面積比最大,其次是喬—灌混合林、灌木、草地,荒地的土壤染色面積比最??;隨演替的進行,SPW<20 mm染色路徑數(shù)量呈增大趨勢,SPW>200 mm染色路徑數(shù)量呈減小趨勢。 (2) 土壤容重、毛管孔隙度、土壤有機質(zhì)含量與土壤染色面積比之間的相關(guān)性較強。其中,荒地和草地中,土壤容重與土壤染色面積比之間存在極顯著負(fù)相關(guān)性(p<0.01);荒地、草地和灌木群落中,土壤有機質(zhì)含量與土壤染色面積比之間存在極顯著正相關(guān)性(p<0.01);喬—灌混合林中,土壤有機質(zhì)含量與土壤染色面積比之間也存在顯著正相關(guān)性(p<0.05)。砂粒含量、黏粒含量與土壤染色面積比之間無顯著相關(guān)關(guān)系。 (3) 植物根系主要分布在0—20 cm表層土壤,各徑級根系的根長密度隨土層深度增加均減小。不同徑級根系對土壤優(yōu)先流路徑影響不同,5種演替梯度下,D≤0.5 mm的根長密度與土壤染色面積比之間相關(guān)性最顯著,其次是0.55 結(jié) 論