杜曉晴,牛健植,侯琨,張英虎,孟晨,駱紫藤
(北京林業(yè)大學水土保持學院,水土保持與荒漠化防治教育部重點實驗室,100083,北京)
華北土石山區(qū)優(yōu)先流區(qū)與基質(zhì)流區(qū)土壤特性分析
杜曉晴,牛健植?,侯琨,張英虎,孟晨,駱紫藤
(北京林業(yè)大學水土保持學院,水土保持與荒漠化防治教育部重點實驗室,100083,北京)
在非飽和土壤中,優(yōu)先流是一種較為常見的現(xiàn)象。針對華北方土石山區(qū)的優(yōu)先流現(xiàn)象,在鷲峰森林公園選取6塊樣地,利用染色示蹤法獲取優(yōu)先流路徑,分析林地優(yōu)先流區(qū)和基質(zhì)流區(qū)土壤孔隙度與水分變化的差異,以及土壤初始含水率對優(yōu)先流(最大染色深度、50%染色深度、染色比例)的影響。結(jié)果表明:1)土壤中優(yōu)先流區(qū)非毛管孔隙度和總孔隙度比基質(zhì)流區(qū)的大,毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)無明顯差異;2)試驗后優(yōu)先流區(qū)土壤含水率變少,含水率減少平均比例是21%,基質(zhì)流區(qū)試驗前后含水率幾乎無變化,含水率變化比例是0;3)初始含水率與染色比例呈正相關(R2=0.679),與50%染色深度呈正相關(R2=0.527),與最大染色深度相關性小(R2=0.236)。
優(yōu)先流; 土壤孔隙度; 土壤含水率變化比例; 染色示蹤法; 初始含水率
優(yōu)先流是一種常見的土壤水分運動現(xiàn)象[1]。土壤學家早在19世紀中后期就已發(fā)現(xiàn)優(yōu)先流[2],但長期以來并沒有得到重視。直到20世紀80年代,人們在美國長島多處飲用水井中發(fā)現(xiàn)用于殺死科羅拉多番茄瓢蟲的殺蟲劑,才開始關注這種對人類健康有潛在威脅的水分快速運移現(xiàn)象[3]。優(yōu)先流是土壤水分運動的重要組成部分,對地表和地下水質(zhì)有重要影響。獲得優(yōu)先流路徑的方法有染色示蹤法等[4]。在野外獲得優(yōu)先流后,一般認為,染色區(qū)是優(yōu)先流區(qū),未染色區(qū)是基質(zhì)流區(qū)[5]。導致優(yōu)先流發(fā)生的因素很多,例如土壤中自然或人為產(chǎn)生的裂隙、根或動物活動形成的孔隙[6]等,并且土壤結(jié)構(gòu)特性、植物根系分布、植被生長狀況[7]及灌溉方法差異[8]等也影響優(yōu)先流區(qū)和基質(zhì)流區(qū)的分布。其中,土壤結(jié)構(gòu)是最本質(zhì)的影響因素[6]。土壤結(jié)構(gòu)指標包括土壤密度、孔隙度、粒度、含水率等,這些差異都會對優(yōu)先流的發(fā)生產(chǎn)生影響。華北土石山區(qū)土少石多,是典型的優(yōu)先流發(fā)生區(qū)域[9],但在該地區(qū)有關土壤特性在優(yōu)先流區(qū)和基質(zhì)流區(qū)的差異很少有人研究。筆者選取位于華北土石山區(qū)的鷲峰森林公園為試驗地,利用染色示蹤法,研究土壤孔隙度和含水率在優(yōu)先流區(qū)和基質(zhì)流區(qū)的差異,以及初始含水率對優(yōu)先流的影響,以期為華北土石山區(qū)優(yōu)先流研究提供科學依據(jù)。
研究區(qū)位于華北土石山區(qū)的北京鷲峰國家森林公園(E116°28′,N39°54′)。鷲峰國家森林公園位于北京城西30 km的群山懷抱之中,南連太行山,北接燕山山脈,總面積832.04 hm2,林地面積811.173 hm2,最高海拔1 153 m。該區(qū)氣候為華北暖溫帶半濕潤半干旱大陸性氣候,年均氣溫12.5 ℃,年最高氣溫為41.6 ℃,最低氣溫為-19.6 ℃,無霜期為193 d。研究區(qū)多年平均降雨量為630 mm,夏季雨水多,春秋干旱,冬季寒冷,年蒸發(fā)量1 800~2 000 mm。土壤類型在海拔70~900 m之間為山地淋溶褐土,900 m以上為棕壤。該區(qū)植物群系屬暖溫帶落葉闊葉林,主要為人工林營造植物群落,主要樹種包括國槐(Sophorajaponica)、側(cè)柏(Platycladusorientalis)、栓皮櫟(Quercusvariabilis)、刺槐(Robiniapseudoacacia)、油松(Pinustabulaeformis)、五角楓(AcermonoMaxim),夾雜有槲櫟(Quercusaliena)、鵝耳櫟(Carpinusturczaninowii)等。
2.1 試驗設計
2013年7月,在鷲峰森林公園選取選取國槐、側(cè)柏、槲櫟3塊主林地,每塊主林地在較平坦地段各選2塊試驗樣地作為試驗地點,共設置6塊樣地(表1)。
在每塊樣地劃定1.2 m×1.2 m區(qū)域,在中心1 m×1 m區(qū)域的周圍插上鐵板,防止施加染色劑后外滲。每塊樣地染色后做垂直剖面,均分成8份,共9層,用數(shù)碼相機拍照染色剖面(圖1)。每個樣地60 cm深,每層10 cm,分6層取土樣,即每個樣地取6個土樣。
2.2 野外示蹤試驗
用噴壺施加染色劑,施加染色劑時噴頭盡量挨近地面,一次性完成,以免噴灑時對地面產(chǎn)生沖擊力。染色24 h后,選擇地面染色區(qū)作為研究范圍,消除邊界影響。根據(jù)M.Bundt等[5]的試驗,利用染色示蹤法表示優(yōu)先流路徑(圖2),染色區(qū)域(黑色)是優(yōu)先流區(qū),未染色區(qū)域(空白)是基質(zhì)流區(qū)。
圖1 田間試驗設計圖Fig.1 Field experimental setup
圖2 優(yōu)先流路徑示意圖Fig.2 Path of preferential flow
2.3 優(yōu)先流評價指標
通過最大染色深度、50%染色深度和染色比例等優(yōu)先流評價指標[10]反映優(yōu)先流運移特性。其中:最大染色深度是入滲后在土壤剖面上達到的最大入滲深度,能夠描述優(yōu)先流的快速運動距離特征;50%染色深度是指在染色面積分數(shù)曲線上染色面積為50%對應的深度,與最大染色深度一樣是描述優(yōu)先流快速運動的一項指標;染色比例指剖面被染色的區(qū)域與總剖面面積之比,表明在該區(qū)域優(yōu)先流發(fā)生面積占整個剖面的比例。
2.4 土壤結(jié)構(gòu)指標測定方法
1)土壤孔隙度。用環(huán)刀(直徑5 cm,高5 cm)分別在優(yōu)先流區(qū)和基質(zhì)流區(qū)取土樣,測定土壤非毛管孔隙度、毛管孔隙度、總孔隙度,并用基質(zhì)流區(qū)孔隙度減優(yōu)先流區(qū)孔隙度得出孔隙度基質(zhì)流區(qū)與優(yōu)先流區(qū)差值:差值大于0,說明基質(zhì)流區(qū)孔隙度大于優(yōu)先流區(qū)孔隙度;差值小于0,說明優(yōu)先流區(qū)孔隙度大于基質(zhì)流區(qū)孔隙度。
2)土壤密度。利用環(huán)刀法測定土壤密度。土壤密度是單位體積土壤(不含孔隙)的烘干質(zhì)量,可間接體現(xiàn)土壤孔隙特征。
3)土壤質(zhì)量含水率(簡稱含水率)。用烘干法分別測定染色前含水率和染色后在優(yōu)先流區(qū)和基質(zhì)流區(qū)含水率,并計算含水率變化比例。
W=(W1-W2)/W2。
式中:W為含水率變化比例,%;W1為試驗后土壤含水率,%;W2為試驗前土壤含水率,%。W>0說明試驗后含水率增加,W<0說明試驗后含水率減少。
4)優(yōu)先流路徑。利用染色示蹤法[4],使用濃度5 g/L的亮藍溶液染色劑進行染色獲取染色路徑,施加染色劑體積見表1。24 h后染色結(jié)束,用直尺測量優(yōu)先流路徑最大深度為最大染色深度;通過繪制染色面積分數(shù)曲線[11]查看50%染色深度,最后計算染色比例。
3.1 優(yōu)先流區(qū)和基質(zhì)流區(qū)土壤孔隙度分析
3.1.1 土壤非毛管孔隙度分析 基質(zhì)流區(qū)與優(yōu)先流區(qū)非毛管孔隙度差值見圖3??梢钥闯?,試驗1到試驗4共4個試驗樣地24個土樣中優(yōu)先流區(qū)大于基質(zhì)流區(qū)的非毛管孔隙度有16個,占66.7%。其中,在20~30 cm和30~40 cm土層上優(yōu)先流區(qū)非毛管孔隙度全都大于基質(zhì)流區(qū)非毛管孔隙度,可知,在20~40 cm土層,非毛管孔隙度越大,越容易發(fā)生優(yōu)先流。這是因為,非毛管孔隙是大孔隙,大孔隙是影響優(yōu)先流發(fā)生的重要因素,水分能在土壤中形成快速運動的大孔隙流。此外,土壤非毛管孔隙度對土壤的通透性和滲透性均具有重要影響,非毛管孔隙越大,土壤通透性和滲透性越好[11]。何東寧等[12]研究得出,非毛管孔隙度增大可以增加林地入滲率和入滲量。該試驗與本試驗結(jié)論相似,在0~40 cm土層上,越往下,非毛管孔隙度在優(yōu)先流區(qū)大于基質(zhì)流區(qū)的現(xiàn)象越明顯。
圖3 基質(zhì)流區(qū)與優(yōu)先流區(qū)非毛管孔隙度差值Fig.3 Difference of non-capillary porosity between matrix flow area and preferential flow area
3.1.2 土壤毛管孔隙度分析 基質(zhì)流區(qū)與優(yōu)先流區(qū)毛管孔隙度差值見圖4??梢钥闯觯?個試驗樣地24個土樣中土壤優(yōu)先流區(qū)與基質(zhì)流區(qū)的毛管孔隙度差異不是很明顯,其中優(yōu)先流區(qū)大于基質(zhì)流區(qū)的毛管孔隙度為13個,占55%;但在不同土層上毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)的關系不同,在0~10 cm土層和30~40 cm土層基質(zhì)流區(qū)毛管孔隙度均大于優(yōu)先流區(qū)毛管孔隙度,在其余土層基本上是優(yōu)先流區(qū)毛管孔隙度大于基質(zhì)流區(qū)。這是因為,毛管孔隙度的直徑很小,對水流的影響差異性很大,對優(yōu)先流的發(fā)生影響較小,但在不同土層上毛管孔隙對優(yōu)先流的影響不同。
圖4 基質(zhì)流區(qū)與優(yōu)先流區(qū)毛管孔隙度差值Fig.4 Difference of capillary porosity between matrix flow area and preferential flow area
3.1.3 土壤總孔隙度分析 基質(zhì)流區(qū)與優(yōu)先流區(qū)總孔隙度差值見圖5。可以看出,4個試驗樣地24個土樣中,優(yōu)先流區(qū)大于基質(zhì)流區(qū)的總管孔隙度為15個,占62.5%,即發(fā)生優(yōu)先流路徑的區(qū)域總孔隙度較大,且10~20 cm層、20~30 cm層、30~40 cm層比較明顯。這是因為,總孔隙度是毛管孔隙度和非毛管孔隙度的總和,直接影響土壤的蓄水性和通氣性[13],總孔隙度的存在能促進優(yōu)先流發(fā)生。
由量化孔隙度可知,非毛管孔隙度和總孔隙度大的土壤更容易發(fā)生優(yōu)先流,毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)的差異較小。非毛管孔隙度和總孔隙度越大越容易發(fā)生優(yōu)先流現(xiàn)象,且在土壤中部比較明顯(比如非毛管孔隙度在20~40 cm、總孔隙度在10~50 cm土層)。盛豐等[10]通過試驗測得總孔隙度隨土層加深逐漸增大,但沒有得出孔隙度與優(yōu)先流的關系。本試驗更詳細地說明了土壤孔隙度對優(yōu)先流的影響。
圖5 基質(zhì)流區(qū)與優(yōu)先流區(qū)總孔隙度差值Fig.5 Difference of total soil porosity between matrix flow area and preferential flow area
3.2 優(yōu)先流區(qū)與基質(zhì)流區(qū)試驗前后含水率變化
3.2.1 優(yōu)先流區(qū)含水率變化 優(yōu)先流區(qū)土壤含水率變化比例見表2??梢钥闯?,優(yōu)先流區(qū)土壤含水率變化比例平均值為-20.14%,經(jīng)計算,平均增加比例是9.70%,平均減少比例是51.16%。試驗1和試驗3比試驗2和試驗4施加染色劑多,施加染色劑多時導致試驗樣地的平均含水率變化比例增大,對優(yōu)先流區(qū)試驗前后含水率的變化影響更大。優(yōu)先流區(qū)含水率有減少趨勢,原因有多種:最大可能是施加染色劑后改變土壤結(jié)構(gòu),增加土壤孔隙的通透性,使孔隙的連通性增強,空氣流動活躍,蒸發(fā)量增大,使含水率減小;另一種可能是,由于試驗地點在土石山區(qū),有可能水流流經(jīng)孔隙時攜帶土壤水入滲到土壤深部,使土壤含水率減小;也可能是水流流動時發(fā)生水分駐留導致有的地方含水率減小。本試驗也可以說明優(yōu)先流的發(fā)生區(qū)域含水率有減小的趨勢,主要原因在于森林表層蓄水減少,地下儲水量增加。
表2 優(yōu)先流區(qū)土壤含水率變化比例
3.2.2 基質(zhì)流區(qū)含水率變化 基質(zhì)流區(qū)土壤含水率變化比例見表3??梢钥闯觯|(zhì)流區(qū)樣地平均含水率變化比例平均值是0。經(jīng)計算,試驗前后含水率平均增加比例是8.25%,平均減少比例是7.58%?;|(zhì)流區(qū)定義為沒有染色的區(qū)域,即施加染色劑后,溶液滲入土壤中時沒有經(jīng)過此區(qū)域,因此可證明試驗結(jié)論是合理的。但是,通過觀察,發(fā)現(xiàn)每塊樣地平均含水率變化比例是有變化的,有增有減,基質(zhì)流區(qū)雖然沒有染色劑流過;而優(yōu)先流區(qū)有水分經(jīng)過,并且水分會揮發(fā),對基質(zhì)流區(qū)土壤含水率存在影響。試驗1和試驗3比試驗2和試驗4施加染色劑多,平均含水率變化比例更大,說明施加染色劑量不同對基質(zhì)流區(qū)含水率變化也是有影響的。可知,基質(zhì)流區(qū)試驗前后含水率變化不大,幾乎為0。
表3 基質(zhì)流區(qū)土壤含水率變化比例
3.3 初始含水率對優(yōu)先流路徑的影響
初始含水率與染色深度以及染色比例的關系見圖6。可知,初始含水率與染色比例呈正相關關系(R2=0.679),與50%染色深度呈正相關關系(R2=0.527),與最大染色深度相關性較小(R2=0.236)。表明初始含水率對染色比例和50%染色深度影響較大,但對最大染色深度的影響較小。這是因為,初始含水率高,土壤孔隙水分含量大,吸附水、毛管水就多,在滲潤和滲漏階段會形成較少的吸附水和毛管水,由于土壤中吸附力和毛管力的作用,使得染色劑發(fā)生水平運動,在重力作用下染色劑會向下運動,水分在水平和垂直運動過程中會形成較大的染色面積,致使對應的50%染色深度和染色比例會增大。
圖6 初始含水率與染色深度以及染色比例的關系Fig.6 Relationship between the antecedent soil moisture and the dyeing depth and the dyeing percentage
不同試驗初始含水率與優(yōu)先流評價指標見表4。可以看出,施加染色劑量不同對優(yōu)先流路徑影響不同,施加染色劑50 L的樣地,最大染色深度、50%染色深度和染色比例均比施加染色劑量30 L的要大??芍┘尤旧珓?yōu)先流路徑的影響明顯。
研究初始含水率對水分運移有一定影響,盛豐等[10]通過試驗得出,土壤初始含水率與土壤入滲深度有相關關系,施加染色劑量相同(20 L)的條件下,初始含水率與土壤入滲深度成正比。這一結(jié)果與本研究結(jié)果相同,但盛豐等[10]在試驗中還得出施加少于20 L染色劑時,初始含水率與土壤入滲深度成反比的結(jié)論。高朝俠等[14]在系統(tǒng)分析初始含水率對優(yōu)先流的影響時也得出,隨著初始含水率的增大濕潤鋒運移距離也增大。土壤含水率對優(yōu)先流的影響也存在爭議,D.B.Jaynes等[15]和K.J.Kung等[16]發(fā)現(xiàn),灌溉能快速增加土壤含水率,同時可增加殺蟲劑運移速度;但有人也得出相反的結(jié)論,如I.Hincapié等[13]研究發(fā)現(xiàn),在土壤較干的時候,施水會促進殺蟲劑運動,在潮濕土壤中會減緩溶質(zhì)運移運動,因此,在不同地區(qū)研究初始含水率對優(yōu)先流路徑的影響很有意義。此外,土壤的優(yōu)先流路徑不僅與初始含水率有關,也與外在因素施加染色劑量有關。在野外自然狀態(tài)下,施加染色劑量等同于降雨量,降雨量也會影響優(yōu)先流的發(fā)生[17]。
表4 不同試驗初始含水率與優(yōu)先流評價指標
1)在試驗區(qū)發(fā)現(xiàn)優(yōu)先流區(qū)和基質(zhì)流區(qū)土壤孔隙度特征有明顯差異。土壤中非毛管孔隙度和總孔隙度大的更容易發(fā)生優(yōu)先流,非毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)差異不是很明顯。
2)試驗前后平均含水率變化比例在優(yōu)先流區(qū)和基質(zhì)流區(qū)有差異。在優(yōu)先流區(qū),試驗后平均含水率變小,平均含水率變化比例為20.14%,基質(zhì)流區(qū)試驗后平均含水率幾乎無變化;但施加染色劑量的多少對每塊樣地的平均含水率變化有影響,平均含水率在施加染色劑50 L的樣地比施加染色劑30L的樣地變化比例大。
3)初始含水率與優(yōu)先流染色比例和50%染色深度呈正相關關系,對最大染色深度的影響較小。施加染色劑越多的樣地,染色深度、50%染色深度和染色比例越大。
土壤孔隙度、土壤密度、初期含水率等對優(yōu)先流有重要影響。在相同質(zhì)地中,土壤密度越大,孔隙率越低,土壤密度越小,孔隙率越大[18],孔隙率與孔隙度密切相關。相同質(zhì)地中,非毛管孔隙越大,初始含水率越低,毛管孔隙越大,初始含水率越大。土壤孔隙對優(yōu)先流路徑有重要影響。本文雖然沒有試驗分析土壤密度對優(yōu)先流的影響,但土壤密度在土壤中的作用是不可忽視的。綜上分析可知,總孔隙越大,越容易發(fā)生優(yōu)先流;因此,可以推測土壤密度越小越容易發(fā)生優(yōu)先流現(xiàn)象,但還須后續(xù)試驗進一步驗證。
由于施加染色劑體積的差異導致試驗重復性小,在今后的研究中,應增加試驗數(shù)量,以便能建立起適合當?shù)氐膬?yōu)先流模型。優(yōu)先流對水文水資源、泥石流滑坡等有重要影響,探索土壤特性對優(yōu)先流的影響,利于人類與自然的和諧相處。
[1] Liu Henry H, Zhang Rui, Bodvarsson G S. An active region model for capturing fractal flow patterns in unsaturated soils: Model development[J]. Journal of contaminant hydrology, 2005, 80(1): 18-30
[2] Gerke H H, Dusek J, Vogel T, et al. Two-dimensional dual-permeability analyses of a bromide tracer experiment on a tile-drained field[J]. Vadose Zone Journal, 2007, 6(3): 651-667
[3] 牛健植, 余新曉. 優(yōu)先流問題研究及其科學意義[J]. 中國水土保持科學, 2005, 3(3): 110-116
[4] Flury M, Wai Na Na. Dyes as tracers for vadose zone hydrology[J]. Reviews of Geophysics, 2003, 41(1):1-37
[5] Bundt M, J?ggi M, Blaser P, et al. Carbon and nitrogen dynamics in preferential flow paths and matrix of a forest soil[J]. Soil Science Society of America Journal, 2001, 65(5): 1529-1538
[6] Jarvis N J. A review of non-equilibrium water flow and solute transport in soil macropores: Principles, controlling factors and consequences for water quality[J]. European Journal of Soil Science, 2007, 58(3): 523-546
[7] 張英虎,牛健植,朱蔚利,等.森林生態(tài)系統(tǒng)林木根系對優(yōu)先流的影響[J].生態(tài)學報,2015,35(6):1-13
[8] Klaus J, Zehe E, Elsner M, et al. Macropore flow of old water revisited: experimental insights from a tile-drained hillslope[J]. Hydrology and Earth System Sciences, 2013, 17(1): 103-118
[9] 張英虎, 牛健植, 杜曉晴, 等. 鷲峰國家森林公園土壤優(yōu)先流現(xiàn)象分析[J]. 水土保持學報, 2013 ,27(1): 41-45
[10] 盛豐, 方妍,張仁鐸. 土壤水非均勻流動的碘- 淀粉染色示蹤研究[J]. 土壤, 2012, 44(1): 144-148
[11] Schaik N Van.Spatial variability of infiltration patterns related to site characteristic in a semi-arid watershed[J]. Catena,2009,78:36-47
[12] 何東寧, 王占林, 張洪勛. 青海樂都地區(qū)森林涵養(yǎng)水源效能研究[J]. 植物生態(tài)學報與地植物學學報, 1991, 15(1): 71-78
[13] Hincapié I, Germann P F. Impact of initial and boundary conditions on preferential flow[J]. Journal of contaminant hydrology, 2009, 104(1): 67-73
[14] 高朝俠, 徐學選, 趙傳普, 等. 土壤初始含水率對優(yōu)先流的影響[J]. 中國水土保持科學, 2014, 12(1): 46-54
[15] Jaynes D B, Ahmed S I, Kung K J S, et al. Temporal dynamics of preferential flow to a subsurface drain[J]. Soil Science Society of America Journal, 2001, 65(5): 1368-1376
[16] Kung K J ,Steenhuis T S, Kladivko E J, et al. Impact of preferential flow on the transport of adsorbing and non-adsorbing tracers[J]. Soil Science Society of America Journal, 2000, 64(4): 1290-1296
[17] 張洪江, 程金花, 史玉虎, 等. 三峽庫區(qū)花崗巖林地坡面優(yōu)先流對降雨的響應[J]. 北京林業(yè)大學學報, 2004, 26(5): 6-9
[18] 郭漢清, 張勇, 李特, 等. 關帝山典型針葉林土壤理化性狀研究[J]. 山西農(nóng)業(yè)大學學報:自然科學版, 2014, 34(2):162-167
(責任編輯:宋如華)
Soil characteristics of preferential flow area and matrix flow area in Rocky Mountain Area of Northern China
Du Xiaoqing, Niu Jianzhi, Hou Kun, Zhang Yinghu, Meng Chen, Luo Ziteng
(School of Soil and Water Conservation, Beijing Forestry University, Key Lab. of Soil & Water Conservation and Desertification Combating of the Ministry of Education, 100083,Beijing,China)
Preferential flow is a major infiltration process in unsaturated soils. We conducted a dye-staining experiment in Beijing Jiufeng National Forest Park, in order to analyze the differences of soil porosity and soil moisture in preferential flow area and matrix flow area, as well as the effect of antecedent soil moisture (ASM) on preferential flow characteristics (i.e., maximum dyed soil depth, 50% dyed soil depth, dyed soil percentage). Results indicated that: 1) the non-capillary porosity and the total porosity of the soil in preferential flow area were statistically higher than those in the matrix area, and no statistical difference in capillary porosity was found between the two areas; 2) soil moisture in the preferential flow area decreased by 21% after the experiment, while there was no difference in soil moisture in the matrix flow area;3) ASM was positively correlated with the dyed soil percentage (R2=0.679) and 50% dyed soil depth (R2=0.527); in contrast, no apparent correlations were found between ASM and the maximum dyed soil depth (R2=0.236).
preferential flow; soil porosity; soil moisture variation; dye staining; antecedent soil moisture
2014-04-02
2015-03-28
項目名稱: 國家自然科學基金“林木根系對溶質(zhì)優(yōu)先運移的影響機制研究”(41271044)
杜曉晴(1988—),女,碩士研究生。主要研究方向:林業(yè)生態(tài)工程與森林水文。E-mail:446378984@qq.com
?通信作者簡介: 牛健植(1974—),女,博士,副教授。主要研究方向:森林生態(tài)水文與水土保持。E-mail: nexk@bjfu.edu.com
S152.4
A
1672-3007(2015)03-0051-07