華北土石山區(qū)優(yōu)先流區(qū)與基質(zhì)流區(qū)土壤特性分析

杜曉晴，牛健植，侯琨，張英虎，孟晨，駱紫藤

(北京林業(yè)大學(xué)水土保持學(xué)院，水土保持與荒漠化防治教育部重點(diǎn)實(shí)驗(yàn)室，100083，北京)

華北土石山區(qū)優(yōu)先流區(qū)與基質(zhì)流區(qū)土壤特性分析

杜曉晴，牛健植?，侯琨，張英虎，孟晨，駱紫藤

(北京林業(yè)大學(xué)水土保持學(xué)院，水土保持與荒漠化防治教育部重點(diǎn)實(shí)驗(yàn)室，100083，北京)

在非飽和土壤中，優(yōu)先流是一種較為常見的現(xiàn)象。針對(duì)華北方土石山區(qū)的優(yōu)先流現(xiàn)象，在鷲峰森林公園選取6塊樣地，利用染色示蹤法獲取優(yōu)先流路徑，分析林地優(yōu)先流區(qū)和基質(zhì)流區(qū)土壤孔隙度與水分變化的差異，以及土壤初始含水率對(duì)優(yōu)先流(最大染色深度、50%染色深度、染色比例)的影響。結(jié)果表明：1)土壤中優(yōu)先流區(qū)非毛管孔隙度和總孔隙度比基質(zhì)流區(qū)的大，毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)無(wú)明顯差異；2)試驗(yàn)后優(yōu)先流區(qū)土壤含水率變少，含水率減少平均比例是21%，基質(zhì)流區(qū)試驗(yàn)前后含水率幾乎無(wú)變化，含水率變化比例是0；3)初始含水率與染色比例呈正相關(guān)(R2=0.679)，與50%染色深度呈正相關(guān)(R2=0.527)，與最大染色深度相關(guān)性小(R2=0.236)。

優(yōu)先流； 土壤孔隙度； 土壤含水率變化比例； 染色示蹤法； 初始含水率

優(yōu)先流是一種常見的土壤水分運(yùn)動(dòng)現(xiàn)象[1]。土壤學(xué)家早在19世紀(jì)中后期就已發(fā)現(xiàn)優(yōu)先流[2]，但長(zhǎng)期以來(lái)并沒有得到重視。直到20世紀(jì)80年代，人們?cè)诿绹?guó)長(zhǎng)島多處飲用水井中發(fā)現(xiàn)用于殺死科羅拉多番茄瓢蟲的殺蟲劑，才開始關(guān)注這種對(duì)人類健康有潛在威脅的水分快速運(yùn)移現(xiàn)象[3]。優(yōu)先流是土壤水分運(yùn)動(dòng)的重要組成部分，對(duì)地表和地下水質(zhì)有重要影響。獲得優(yōu)先流路徑的方法有染色示蹤法等[4]。在野外獲得優(yōu)先流后，一般認(rèn)為，染色區(qū)是優(yōu)先流區(qū)，未染色區(qū)是基質(zhì)流區(qū)[5]。導(dǎo)致優(yōu)先流發(fā)生的因素很多，例如土壤中自然或人為產(chǎn)生的裂隙、根或動(dòng)物活動(dòng)形成的孔隙[6]等，并且土壤結(jié)構(gòu)特性、植物根系分布、植被生長(zhǎng)狀況[7]及灌溉方法差異[8]等也影響優(yōu)先流區(qū)和基質(zhì)流區(qū)的分布。其中，土壤結(jié)構(gòu)是最本質(zhì)的影響因素[6]。土壤結(jié)構(gòu)指標(biāo)包括土壤密度、孔隙度、粒度、含水率等，這些差異都會(huì)對(duì)優(yōu)先流的發(fā)生產(chǎn)生影響。華北土石山區(qū)土少石多，是典型的優(yōu)先流發(fā)生區(qū)域[9]，但在該地區(qū)有關(guān)土壤特性在優(yōu)先流區(qū)和基質(zhì)流區(qū)的差異很少有人研究。筆者選取位于華北土石山區(qū)的鷲峰森林公園為試驗(yàn)地，利用染色示蹤法，研究土壤孔隙度和含水率在優(yōu)先流區(qū)和基質(zhì)流區(qū)的差異，以及初始含水率對(duì)優(yōu)先流的影響，以期為華北土石山區(qū)優(yōu)先流研究提供科學(xué)依據(jù)。

1 研究區(qū)概況

研究區(qū)位于華北土石山區(qū)的北京鷲峰國(guó)家森林公園(E116°28′，N39°54′)。鷲峰國(guó)家森林公園位于北京城西30 km的群山懷抱之中，南連太行山，北接燕山山脈，總面積832.04 hm2，林地面積811.173 hm2，最高海拔1 153 m。該區(qū)氣候?yàn)槿A北暖溫帶半濕潤(rùn)半干旱大陸性氣候，年均氣溫12.5 ℃，年最高氣溫為41.6 ℃，最低氣溫為-19.6 ℃，無(wú)霜期為193 d。研究區(qū)多年平均降雨量為630 mm，夏季雨水多，春秋干旱，冬季寒冷，年蒸發(fā)量1 800～2 000 mm。土壤類型在海拔70～900 m之間為山地淋溶褐土，900 m以上為棕壤。該區(qū)植物群系屬暖溫帶落葉闊葉林，主要為人工林營(yíng)造植物群落，主要樹種包括國(guó)槐(Sophorajaponica)、側(cè)柏(Platycladusorientalis)、栓皮櫟(Quercusvariabilis)、刺槐(Robiniapseudoacacia)、油松(Pinustabulaeformis)、五角楓(AcermonoMaxim)，夾雜有槲櫟(Quercusaliena)、鵝耳櫟(Carpinusturczaninowii)等。

2 研究方法

2.1 試驗(yàn)設(shè)計(jì)

2013年7月，在鷲峰森林公園選取選取國(guó)槐、側(cè)柏、槲櫟3塊主林地，每塊主林地在較平坦地段各選2塊試驗(yàn)樣地作為試驗(yàn)地點(diǎn)，共設(shè)置6塊樣地(表1)。

在每塊樣地劃定1.2 m×1.2 m區(qū)域，在中心1 m×1 m區(qū)域的周圍插上鐵板，防止施加染色劑后外滲。每塊樣地染色后做垂直剖面，均分成8份，共9層，用數(shù)碼相機(jī)拍照染色剖面(圖1)。每個(gè)樣地60 cm深，每層10 cm，分6層取土樣，即每個(gè)樣地取6個(gè)土樣。

2.2 野外示蹤試驗(yàn)

用噴壺施加染色劑，施加染色劑時(shí)噴頭盡量挨近地面，一次性完成，以免噴灑時(shí)對(duì)地面產(chǎn)生沖擊力。染色24 h后，選擇地面染色區(qū)作為研究范圍，消除邊界影響。根據(jù)M.Bundt等[5]的試驗(yàn)，利用染色示蹤法表示優(yōu)先流路徑(圖2)，染色區(qū)域(黑色)是優(yōu)先流區(qū)，未染色區(qū)域(空白)是基質(zhì)流區(qū)。

圖1 田間試驗(yàn)設(shè)計(jì)圖Fig.1 Field experimental setup

圖2 優(yōu)先流路徑示意圖Fig.2 Path of preferential flow

2.3 優(yōu)先流評(píng)價(jià)指標(biāo)

通過最大染色深度、50%染色深度和染色比例等優(yōu)先流評(píng)價(jià)指標(biāo)[10]反映優(yōu)先流運(yùn)移特性。其中：最大染色深度是入滲后在土壤剖面上達(dá)到的最大入滲深度，能夠描述優(yōu)先流的快速運(yùn)動(dòng)距離特征；50%染色深度是指在染色面積分?jǐn)?shù)曲線上染色面積為50%對(duì)應(yīng)的深度，與最大染色深度一樣是描述優(yōu)先流快速運(yùn)動(dòng)的一項(xiàng)指標(biāo)；染色比例指剖面被染色的區(qū)域與總剖面面積之比，表明在該區(qū)域優(yōu)先流發(fā)生面積占整個(gè)剖面的比例。

2.4 土壤結(jié)構(gòu)指標(biāo)測(cè)定方法

1)土壤孔隙度。用環(huán)刀(直徑5 cm，高5 cm)分別在優(yōu)先流區(qū)和基質(zhì)流區(qū)取土樣，測(cè)定土壤非毛管孔隙度、毛管孔隙度、總孔隙度，并用基質(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)刀法測(cè)定土壤密度。土壤密度是單位體積土壤(不含孔隙)的烘干質(zhì)量，可間接體現(xiàn)土壤孔隙特征。

3)土壤質(zhì)量含水率(簡(jiǎn)稱含水率)。用烘干法分別測(cè)定染色前含水率和染色后在優(yōu)先流區(qū)和基質(zhì)流區(qū)含水率，并計(jì)算含水率變化比例。

W=(W1-W2)/W2。

式中：W為含水率變化比例，%；W1為試驗(yàn)后土壤含水率，%；W2為試驗(yàn)前土壤含水率，%。W>0說明試驗(yàn)后含水率增加，W<0說明試驗(yàn)后含水率減少。

4)優(yōu)先流路徑。利用染色示蹤法[4]，使用濃度5 g/L的亮藍(lán)溶液染色劑進(jìn)行染色獲取染色路徑，施加染色劑體積見表1。24 h后染色結(jié)束，用直尺測(cè)量?jī)?yōu)先流路徑最大深度為最大染色深度；通過繪制染色面積分?jǐn)?shù)曲線[11]查看50%染色深度，最后計(jì)算染色比例。

3 結(jié)果與分析

3.1 優(yōu)先流區(qū)和基質(zhì)流區(qū)土壤孔隙度分析

3.1.1 土壤非毛管孔隙度分析 基質(zhì)流區(qū)與優(yōu)先流區(qū)非毛管孔隙度差值見圖3?？梢钥闯?，試驗(yàn)1到試驗(yàn)4共4個(gè)試驗(yàn)樣地24個(gè)土樣中優(yōu)先流區(qū)大于基質(zhì)流區(qū)的非毛管孔隙度有16個(gè)，占66.7%。其中，在20～30 cm和30～40 cm土層上優(yōu)先流區(qū)非毛管孔隙度全都大于基質(zhì)流區(qū)非毛管孔隙度，可知，在20～40 cm土層，非毛管孔隙度越大，越容易發(fā)生優(yōu)先流。這是因?yàn)?，非毛管孔隙是大孔隙，大孔隙是影響?yōu)先流發(fā)生的重要因素，水分能在土壤中形成快速運(yùn)動(dòng)的大孔隙流。此外，土壤非毛管孔隙度對(duì)土壤的通透性和滲透性均具有重要影響，非毛管孔隙越大，土壤通透性和滲透性越好[11]。何東寧等[12]研究得出，非毛管孔隙度增大可以增加林地入滲率和入滲量。該試驗(yàn)與本試驗(yàn)結(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?？梢钥闯觯?個(gè)試驗(yàn)樣地24個(gè)土樣中土壤優(yōu)先流區(qū)與基質(zhì)流區(qū)的毛管孔隙度差異不是很明顯，其中優(yōu)先流區(qū)大于基質(zhì)流區(qū)的毛管孔隙度為13個(gè)，占55%；但在不同土層上毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)的關(guān)系不同，在0～10 cm土層和30～40 cm土層基質(zhì)流區(qū)毛管孔隙度均大于優(yōu)先流區(qū)毛管孔隙度，在其余土層基本上是優(yōu)先流區(qū)毛管孔隙度大于基質(zhì)流區(qū)。這是因?yàn)?，毛管孔隙度的直徑很小，?duì)水流的影響差異性很大，對(duì)優(yōu)先流的發(fā)生影響較小，但在不同土層上毛管孔隙對(duì)優(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個(gè)試驗(yàn)樣地24個(gè)土樣中，優(yōu)先流區(qū)大于基質(zhì)流區(qū)的總管孔隙度為15個(gè)，占62.5%，即發(fā)生優(yōu)先流路徑的區(qū)域總孔隙度較大，且10～20 cm層、20～30 cm層、30～40 cm層比較明顯。這是因?yàn)?，總孔隙度是毛管孔隙度和非毛管孔隙度的總和，直接影響土壤的蓄水性和通氣性[13]，總孔隙度的存在能促進(jìn)優(yōu)先流發(fā)生。

由量化孔隙度可知，非毛管孔隙度和總孔隙度大的土壤更容易發(fā)生優(yōu)先流，毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)的差異較小。非毛管孔隙度和總孔隙度越大越容易發(fā)生優(yōu)先流現(xiàn)象，且在土壤中部比較明顯(比如非毛管孔隙度在20～40 cm、總孔隙度在10～50 cm土層)。盛豐等[10]通過試驗(yàn)測(cè)得總孔隙度隨土層加深逐漸增大，但沒有得出孔隙度與優(yōu)先流的關(guān)系。本試驗(yàn)更詳細(xì)地說明了土壤孔隙度對(duì)優(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ū)試驗(yàn)前后含水率變化

3.2.1 優(yōu)先流區(qū)含水率變化 優(yōu)先流區(qū)土壤含水率變化比例見表2?？梢钥闯?，優(yōu)先流區(qū)土壤含水率變化比例平均值為-20.14%，經(jīng)計(jì)算，平均增加比例是9.70%，平均減少比例是51.16%。試驗(yàn)1和試驗(yàn)3比試驗(yàn)2和試驗(yàn)4施加染色劑多，施加染色劑多時(shí)導(dǎo)致試驗(yàn)樣地的平均含水率變化比例增大，對(duì)優(yōu)先流區(qū)試驗(yàn)前后含水率的變化影響更大。優(yōu)先流區(qū)含水率有減少趨勢(shì)，原因有多種：最大可能是施加染色劑后改變土壤結(jié)構(gòu)，增加土壤孔隙的通透性，使孔隙的連通性增強(qiáng)，空氣流動(dòng)活躍，蒸發(fā)量增大，使含水率減??；另一種可能是，由于試驗(yàn)地點(diǎn)在土石山區(qū)，有可能水流流經(jīng)孔隙時(shí)攜帶土壤水入滲到土壤深部，使土壤含水率減??；也可能是水流流動(dòng)時(shí)發(fā)生水分駐留導(dǎo)致有的地方含水率減小。本試驗(yàn)也可以說明優(yōu)先流的發(fā)生區(qū)域含水率有減小的趨勢(shì)，主要原因在于森林表層蓄水減少，地下儲(chǔ)水量增加。

表2 優(yōu)先流區(qū)土壤含水率變化比例

3.2.2 基質(zhì)流區(qū)含水率變化 基質(zhì)流區(qū)土壤含水率變化比例見表3。可以看出，基質(zhì)流區(qū)樣地平均含水率變化比例平均值是0。經(jīng)計(jì)算，試驗(yàn)前后含水率平均增加比例是8.25%，平均減少比例是7.58%。基質(zhì)流區(qū)定義為沒有染色的區(qū)域，即施加染色劑后，溶液滲入土壤中時(shí)沒有經(jīng)過此區(qū)域，因此可證明試驗(yàn)結(jié)論是合理的。但是，通過觀察，發(fā)現(xiàn)每塊樣地平均含水率變化比例是有變化的，有增有減，基質(zhì)流區(qū)雖然沒有染色劑流過；而優(yōu)先流區(qū)有水分經(jīng)過，并且水分會(huì)揮發(fā)，對(duì)基質(zhì)流區(qū)土壤含水率存在影響。試驗(yàn)1和試驗(yàn)3比試驗(yàn)2和試驗(yàn)4施加染色劑多，平均含水率變化比例更大，說明施加染色劑量不同對(duì)基質(zhì)流區(qū)含水率變化也是有影響的?？芍|(zhì)流區(qū)試驗(yàn)前后含水率變化不大，幾乎為0。

表3 基質(zhì)流區(qū)土壤含水率變化比例

3.3 初始含水率對(duì)優(yōu)先流路徑的影響

初始含水率與染色深度以及染色比例的關(guān)系見圖6?？芍?，初始含水率與染色比例呈正相關(guān)關(guān)系(R2=0.679)，與50%染色深度呈正相關(guān)關(guān)系(R2=0.527)，與最大染色深度相關(guān)性較小(R2=0.236)。表明初始含水率對(duì)染色比例和50%染色深度影響較大，但對(duì)最大染色深度的影響較小。這是因?yàn)?，初始含水率高，土壤孔隙水分含量大，吸附水、毛管水就多，在滲潤(rùn)和滲漏階段會(huì)形成較少的吸附水和毛管水，由于土壤中吸附力和毛管力的作用，使得染色劑發(fā)生水平運(yùn)動(dòng)，在重力作用下染色劑會(huì)向下運(yùn)動(dòng)，水分在水平和垂直運(yùn)動(dòng)過程中會(huì)形成較大的染色面積，致使對(duì)應(yīng)的50%染色深度和染色比例會(huì)增大。

圖6 初始含水率與染色深度以及染色比例的關(guān)系Fig.6 Relationship between the antecedent soil moisture and the dyeing depth and the dyeing percentage

不同試驗(yàn)初始含水率與優(yōu)先流評(píng)價(jià)指標(biāo)見表4。可以看出，施加染色劑量不同對(duì)優(yōu)先流路徑影響不同，施加染色劑50 L的樣地，最大染色深度、50%染色深度和染色比例均比施加染色劑量30 L的要大。可知，施加染色劑對(duì)優(yōu)先流路徑的影響明顯。

研究初始含水率對(duì)水分運(yùn)移有一定影響，盛豐等[10]通過試驗(yàn)得出，土壤初始含水率與土壤入滲深度有相關(guān)關(guān)系，施加染色劑量相同(20 L)的條件下，初始含水率與土壤入滲深度成正比。這一結(jié)果與本研究結(jié)果相同，但盛豐等[10]在試驗(yàn)中還得出施加少于20 L染色劑時(shí)，初始含水率與土壤入滲深度成反比的結(jié)論。高朝俠等[14]在系統(tǒng)分析初始含水率對(duì)優(yōu)先流的影響時(shí)也得出，隨著初始含水率的增大濕潤(rùn)鋒運(yùn)移距離也增大。土壤含水率對(duì)優(yōu)先流的影響也存在爭(zhēng)議，D.B.Jaynes等[15]和K.J.Kung等[16]發(fā)現(xiàn)，灌溉能快速增加土壤含水率，同時(shí)可增加殺蟲劑運(yùn)移速度；但有人也得出相反的結(jié)論，如I.Hincapié等[13]研究發(fā)現(xiàn)，在土壤較干的時(shí)候，施水會(huì)促進(jìn)殺蟲劑運(yùn)動(dòng)，在潮濕土壤中會(huì)減緩溶質(zhì)運(yùn)移運(yùn)動(dòng)，因此，在不同地區(qū)研究初始含水率對(duì)優(yōu)先流路徑的影響很有意義。此外，土壤的優(yōu)先流路徑不僅與初始含水率有關(guān)，也與外在因素施加染色劑量有關(guān)。在野外自然狀態(tài)下，施加染色劑量等同于降雨量，降雨量也會(huì)影響優(yōu)先流的發(fā)生[17]。

表4 不同試驗(yàn)初始含水率與優(yōu)先流評(píng)價(jià)指標(biāo)

4 結(jié)論與討論

1)在試驗(yàn)區(qū)發(fā)現(xiàn)優(yōu)先流區(qū)和基質(zhì)流區(qū)土壤孔隙度特征有明顯差異。土壤中非毛管孔隙度和總孔隙度大的更容易發(fā)生優(yōu)先流，非毛管孔隙度在優(yōu)先流區(qū)與基質(zhì)流區(qū)差異不是很明顯。

2)試驗(yàn)前后平均含水率變化比例在優(yōu)先流區(qū)和基質(zhì)流區(qū)有差異。在優(yōu)先流區(qū)，試驗(yàn)后平均含水率變小，平均含水率變化比例為20.14%，基質(zhì)流區(qū)試驗(yàn)后平均含水率幾乎無(wú)變化；但施加染色劑量的多少對(duì)每塊樣地的平均含水率變化有影響，平均含水率在施加染色劑50 L的樣地比施加染色劑30L的樣地變化比例大。

3)初始含水率與優(yōu)先流染色比例和50%染色深度呈正相關(guān)關(guān)系，對(duì)最大染色深度的影響較小。施加染色劑越多的樣地，染色深度、50%染色深度和染色比例越大。

土壤孔隙度、土壤密度、初期含水率等對(duì)優(yōu)先流有重要影響。在相同質(zhì)地中，土壤密度越大，孔隙率越低，土壤密度越小，孔隙率越大[18]，孔隙率與孔隙度密切相關(guān)。相同質(zhì)地中，非毛管孔隙越大，初始含水率越低，毛管孔隙越大，初始含水率越大。土壤孔隙對(duì)優(yōu)先流路徑有重要影響。本文雖然沒有試驗(yàn)分析土壤密度對(duì)優(yōu)先流的影響，但土壤密度在土壤中的作用是不可忽視的。綜上分析可知，總孔隙越大，越容易發(fā)生優(yōu)先流；因此，可以推測(cè)土壤密度越小越容易發(fā)生優(yōu)先流現(xiàn)象，但還須后續(xù)試驗(yàn)進(jìn)一步驗(yàn)證。

由于施加染色劑體積的差異導(dǎo)致試驗(yàn)重復(fù)性小，在今后的研究中，應(yīng)增加試驗(yàn)數(shù)量，以便能建立起適合當(dāng)?shù)氐膬?yōu)先流模型。優(yōu)先流對(duì)水文水資源、泥石流滑坡等有重要影響，探索土壤特性對(duì)優(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)先流問題研究及其科學(xué)意義[J]. 中國(guó)水土保持科學(xué), 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)林木根系對(duì)優(yōu)先流的影響[J].生態(tài)學(xué)報(bào),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] 張英虎, 牛健植, 杜曉晴, 等. 鷲峰國(guó)家森林公園土壤優(yōu)先流現(xiàn)象分析[J]. 水土保持學(xué)報(bào), 2013 ，27(1): 41-45

[10] 盛豐, 方妍，張仁鐸. 土壤水非均勻流動(dòng)的碘- 淀粉染色示蹤研究[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)學(xué)報(bào)與地植物學(xué)學(xué)報(bào), 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] 高朝俠, 徐學(xué)選, 趙傳普, 等. 土壤初始含水率對(duì)優(yōu)先流的影響[J]. 中國(guó)水土保持科學(xué), 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] 張洪江, 程金花, 史玉虎, 等. 三峽庫(kù)區(qū)花崗巖林地坡面優(yōu)先流對(duì)降雨的響應(yīng)[J]. 北京林業(yè)大學(xué)學(xué)報(bào), 2004, 26(5): 6-9

[18] 郭漢清, 張勇, 李特, 等. 關(guān)帝山典型針葉林土壤理化性狀研究[J]. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版, 2014, 34(2)：162-167

(責(zé)任編輯：宋如華)

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

項(xiàng)目名稱： 國(guó)家自然科學(xué)基金“林木根系對(duì)溶質(zhì)優(yōu)先運(yùn)移的影響機(jī)制研究”(41271044)

杜曉晴(1988—)，女，碩士研究生。主要研究方向：林業(yè)生態(tài)工程與森林水文。E-mail：446378984@qq.com

?通信作者簡(jiǎn)介： 牛健植(1974—)，女，博士，副教授。主要研究方向：森林生態(tài)水文與水土保持。E-mail: nexk@bjfu.edu.com

S152.4

A

1672-3007(2015)03-0051-07