陸 斌，張勝利，* ，李 侃，馬國(guó)棟

(1.西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院，楊凌 712100;2.陜西秦嶺森林生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站，楊凌 712100;3.西北農(nóng)林科技大學(xué)水土保持研究所，楊凌 712100)

壤中流的產(chǎn)生是森林涵養(yǎng)水源和調(diào)節(jié)徑流的主要方式，大孔隙在其中起著重要的作用。大孔隙可促進(jìn)側(cè)向表層流和快速旁路流，被認(rèn)為是為坡面淺土層暴雨洪峰提供重要優(yōu)先流的路徑［1-2］;Wilson［3］、Shipitalo［4］、Weiler［5］等研究發(fā)現(xiàn)，雖然土壤中大孔隙數(shù)量很少，但對(duì)傳導(dǎo)水流有著重要作用，可明顯地增加入滲，減少地表徑流。Carneira等［6］研究發(fā)現(xiàn)，大孔隙的存在使土壤的飽和導(dǎo)水率增加了幾倍，且對(duì)總?cè)霛B量貢獻(xiàn)率可達(dá)到 85%;Wilson、Noguchi［6］等認(rèn)為直徑大于l mm的孔隙對(duì)傳導(dǎo)水流有著重要的因素;而Luxmoore［7］發(fā)現(xiàn)大中孔隙的不同組合有著不同的導(dǎo)水速率。土壤中大孔隙的孔徑及數(shù)量特征直接反映了林地優(yōu)先路徑的組成情況，研究林地土壤大孔隙分布規(guī)律有助于理解林區(qū)土壤的導(dǎo)水性能和壤中流的發(fā)展過(guò)程。

目前有關(guān)大孔隙的研究并未給出嚴(yán)格的定義和孔徑范圍，一般將孔徑大于0.03 mm，0.06 mm，0.085 mm，0.1 mm或1 mm的孔隙定義為大孔隙;孔徑范圍大致為 0.03—3 mm［1，6，8-9］。受動(dòng)物、植被根系或土壤自身等因素的影響，不同地區(qū)不同植被類型的大孔隙分布具有不同特點(diǎn)，這也導(dǎo)致了土壤水分運(yùn)動(dòng)空間上的差異［10-11］。所以，要進(jìn)一步完善林地土壤水分運(yùn)動(dòng)及流域生態(tài)水文過(guò)程理論，必須深入開(kāi)展不同地區(qū)多重小尺度上有關(guān)大孔隙的研究。

本研究選擇在秦嶺南坡火地塘林場(chǎng)進(jìn)行。秦嶺南坡林區(qū)的水源供給集中了陜西省73%的水資源，其涵養(yǎng)水源狀況對(duì)周邊及其下游地區(qū)社會(huì)和經(jīng)濟(jì)發(fā)展影響重大。在林區(qū)，徑流多以壤中流的形式匯集，較大空隙的優(yōu)先流運(yùn)動(dòng)起到了至關(guān)重要的作用［12-14］。因此，系統(tǒng)的分析研究該區(qū)林地土壤大孔隙組成與數(shù)量分布規(guī)律及其對(duì)土壤水分運(yùn)移過(guò)程的影響，有助于進(jìn)一步理解流域水分運(yùn)動(dòng)規(guī)律以及林區(qū)壤中流的產(chǎn)匯流特征。同時(shí)，該區(qū)也是我國(guó)南水北調(diào)中線工程的主要水源地，故該研究還可為調(diào)水區(qū)水資源管理提供科學(xué)依據(jù)。

1 研究區(qū)概況及研究方法

1.1 研究區(qū)概況

火地塘林區(qū)位于秦嶺南坡中山地帶中心地段，選擇在火地塘林區(qū)開(kāi)展森林大孔隙的相關(guān)研究極具典型性和代表性［13］。該區(qū)地處北緯 33°18'—33°28'，東經(jīng) 108°20'—108°39'，屬長(zhǎng)江水系漢江流域。林區(qū)海拔范圍900—2450 m，年平均氣溫8—10℃，年平均降水900—1200 mm，年水面蒸發(fā)量800—950 mm，年總?cè)照諘r(shí)數(shù)為1100—1300 h。林區(qū)腐殖質(zhì)層厚約7—10 cm，森林覆蓋率達(dá)91.8%。主要成林樹(shù)種有銳齒櫟(Quercus aliena var.acuteserrata)、油松(Pinus tabulaeformis)、華山松(Pinus armandii)、紅樺(Betula albo-sinensis)、光皮樺(Butula luminifera);伴生樹(shù)種有:漆樹(shù)(Toxicodendron)、鵝耳櫟(Carpinus turczanoinwii)、木姜子 (Litsea phngens)、青榨械(Acerdavidii);主要灌木有:榛子(Corylus heterophylla)、山柳(Salix)、刺懸鉤子(Rubus phngens)、胡枝子(Lespedeza bicolor Thunb)、山馬蝗(Desmodium sp.)、衛(wèi)矛(Euonymus alata)等。

試驗(yàn)區(qū)土壤以山地棕壤為主，平均厚度為50 cm。成土母質(zhì)以片麻巖、花崗巖、結(jié)晶片巖及似斑狀花崗閃長(zhǎng)巖風(fēng)化的坡積和殘積物為主。土壤中石礫分布較多且形狀差異較大，油松林地石礫多介于2—30 mm，平均體積含量為21.08%;針闊混交林和闊葉林地，≤30 mm的礫石含量在4%—8%之間，且多鑲嵌有較大直徑的半分化巖石。林地土壤中根系較豐富，多分布于0—30 cm土層，直徑≤1 mm的根系生物量最大。其中針闊混交林和闊葉林地表層的腐根較多。

1.2 研究方法

1.2.1 樣地設(shè)置

實(shí)驗(yàn)選在林區(qū)最大的自然集水區(qū)火地溝流域進(jìn)行。該區(qū)森林植被、地形地貌、土壤等具有秦嶺南坡中山地帶的典型特征。區(qū)內(nèi)植被為松櫟林和松樺林為主的落葉闊葉和針闊葉混交林，主要樹(shù)種有銳齒櫟、山楊、白樺、油松和華山松等［13］。根據(jù)林區(qū)的主要成林樹(shù)種分布，本研究設(shè)置了4種主要林分類型(表 1)。

表1 研究地土壤(0—60 cm)的基本特征Table 1 Selected soil physical and chemical properties of the experimental site(mean±SD)

1.2.2 大孔隙的計(jì)算方法

測(cè)量土壤大孔隙的方法有很多種(如示蹤技術(shù)，土壤切片法，穿透曲線法，張力入滲儀測(cè)量法，X射線CT掃描攝像以及雷達(dá)探測(cè)技術(shù)等)［15-19］。從森林土壤大孔隙水文效應(yīng)的角度出發(fā)，采用土壤水分穿透曲線和Poiseulle方程相結(jié)合的方法較合適［7］。該法測(cè)定的是“瓶頸”大小孔隙，而在微小尺度上，一個(gè)連續(xù)水流路徑的出流速率受最小尺寸空隙(亦稱為pore necks)控制［1］。與其它測(cè)定方法相比，該法所得參數(shù)更適合于描述大孔隙與土壤水分運(yùn)動(dòng)的關(guān)系，并適于林地多根系多石礫土壤的研究［8］。

根據(jù)時(shí)忠杰［8］、石輝［15］、Radulovich［17］等的研究，大孔隙定義為田間持水量到飽和含水量之間的土壤孔隙。田間持水量到飽和含水量之間的排水過(guò)程首先是大孔隙排水，然后排水孔隙的孔徑逐漸減小。對(duì)某一固定土樣，最大孔隙半徑就是第一次出水時(shí)的半徑;而土樣的最小半徑就是當(dāng)水流開(kāi)始達(dá)到穩(wěn)定時(shí)的半徑［9，15］。試驗(yàn)中土壤水分運(yùn)動(dòng)速率較慢，處于層流范圍，假設(shè)土壤孔隙為圓形的情況下，可根據(jù)土壤水分穿透曲線理論，將流量方程(1)與Poiseulle方程(2)聯(lián)合求解大孔隙的當(dāng)量孔徑(3):

式中，Q為單位流量(cm3/s);A為孔隙面積(cm2);v為流速(cm/s);r為孔隙孔徑(cm);τ為水流實(shí)際路徑的彎曲系數(shù)，本研究取1.2;L為土柱長(zhǎng)度(cm);η為水的粘滯系數(shù)(g/cm/s);ΔP為壓力水頭(cm);t為從首次加水開(kāi)始記時(shí)的時(shí)間(s)。對(duì)任意時(shí)間水分出流流量進(jìn)行觀測(cè)，并利用式(3)計(jì)算出相應(yīng)的孔隙孔徑。將兩孔隙孔徑之間按照一定間隔劃分出一個(gè)孔隙范圍，以其平均值作為計(jì)算值，當(dāng)間隔排水量為Qe時(shí)，計(jì)算相應(yīng)大孔隙數(shù)量(n):

1.2.3 土樣采集與出流速率的測(cè)定。

2012年5—8月，用環(huán)刀(直徑7 cm、高5.2 cm)分層(A層約0—20 cm、B層約20—40 cm、BC層約40—60 cm)采集土壤樣品，平放并盡量減少顛簸以保持土樣的原狀，帶回實(shí)驗(yàn)室測(cè)定水分穿透曲線。每種植被類型設(shè)3—5個(gè)樣點(diǎn)。將采集的原狀土壤樣品置于水中12 h以充分飽和，然后放置于粗砂上放置12 h，使其達(dá)到田間持水量，并在環(huán)刀上邊緣與土壤樣品接觸的部分涂上凡士林，以保證水分不能通過(guò)土壤與環(huán)刀壁之間的孔隙直接流出。在土柱上部放置一張濾紙，以防止水分沖刷破壞表層土壤結(jié)構(gòu)。土柱下部過(guò)濾板的出流孔密度為5個(gè)/cm2。用馬氏瓶控制水頭2 cm，從加水開(kāi)始計(jì)時(shí)，當(dāng)土柱下部有水流出時(shí)，每隔5s或10s收集出流水量，測(cè)定其體積，直到水流通量達(dá)到穩(wěn)定為止。然后根據(jù)測(cè)定的出流水量作水分穿透曲線，確定其水流通量達(dá)到穩(wěn)定的時(shí)間，再參照1.2.2計(jì)算大孔隙的相關(guān)指標(biāo)。根據(jù)水流穩(wěn)定后的流量計(jì)算飽和導(dǎo)水率(Ks)，統(tǒng)一換算為10℃的飽和導(dǎo)水率(K10)。

土壤有機(jī)質(zhì)含量測(cè)定采用重鉻酸鉀法。

2 結(jié)果與分析

2.1 不同植被類型大孔隙分布特征

研究區(qū)土壤大孔隙的當(dāng)量孔徑范圍為0.3—3.8 mm(表2)，該計(jì)算結(jié)果與王偉［18］和 Germann 等［1］所測(cè)土壤大孔隙孔徑范圍基本一致。最小孔徑0.3 mm，與土壤學(xué)中土壤通氣孔中的粗孔隙孔徑下限(孔徑大于0.3或0.2 mm)一致［19］，說(shuō)明研究地森林土壤大孔隙主要為孔徑大于0.3 mm的通氣孔隙。

表2 不同當(dāng)量孔徑的大孔隙密度Table 2 Number of soil macropores in each radius interval

大孔隙密度為土壤單位斷面積上的平均大孔隙數(shù)量。火地塘林區(qū)，0—60 cm土壤大孔隙密度在0.14×104—7.18×104個(gè)/m2之間，平均為 1.98×104個(gè)/m2，變異系數(shù)為0.42。不同樣地表層(0—20 cm)土壤大孔隙平均密度大小順序?yàn)?針闊混交林＞落葉闊葉林＞華山松林＞油松林;0—60 cm土層，大小順序則為針闊混交林＞華山松林＞落葉闊葉林＞油松林。針闊混交林地土壤大孔隙密度最大。0—60 cm土層的大孔隙平均密度的變異系數(shù)大小順序?yàn)?針闊混交林＞華山松林＞落葉闊葉林＞油松林。土壤大孔隙斷面面積與土柱總過(guò)水?dāng)嗝婷娣e的百分比(簡(jiǎn)稱大孔隙面積比)能夠綜合反映大孔隙半徑及密度等因素對(duì)土壤水分運(yùn)動(dòng)的影響。各樣地大孔隙面積比平均為4.59%，標(biāo)準(zhǔn)差為2.81%，變異系數(shù)為0.61。大孔隙面積比最高的為油松林地的A層，較小的為油松林的C層和華山松林的C層。0—60 cm土層大孔隙平均面積比的順序?yàn)?針闊混交林＞油松林＞落葉闊葉林＞華山松林。油松林地的大孔隙面積比較大可能是土壤中礫石分布影響較大半徑大孔隙密度的緣故;華山松林地石礫含量較低且質(zhì)地較黏重，可能對(duì)大孔隙的發(fā)育有一定限制。針闊混交林和落葉闊葉林的土壤大孔隙平均密度和面積比的變異系數(shù)均較大，這可能與該林分類型較為豐富的腐爛根系有一定關(guān)系，也可能是由于所處林地土壤表層鑲嵌的石塊以及底部散布的礫石所致。不同類型林地、不同土壤層中大孔隙密度變化具有較大差異，說(shuō)明植物組成對(duì)土壤大孔隙特征具有較為顯著的影響。

2.2 大孔隙垂直分布規(guī)律

隨著林地土壤深度的增加，地表植物根系能夠作用到的范圍逐漸縮小，植物對(duì)土壤結(jié)構(gòu)的影響程度也逐漸降低［20］。缺少了植物根系的擠壓與切割以及腐根對(duì)土體的疏松，林地底層的土壤一般會(huì)較其上層變得更加緊實(shí)，從而影響著土壤大孔隙的形成。各林地土壤層中，較大孔徑孔隙的密度較小，而較小孔徑孔隙的密度較大，且各當(dāng)量孔徑的大孔隙密度隨土層分布基本表現(xiàn)為上層大、下層小的特點(diǎn)(表2)。采用新復(fù)極差法的多重比較統(tǒng)計(jì)(P＜0.05)對(duì)不同林分不同土層深度中各當(dāng)量孔徑的大孔隙密度進(jìn)行分析，表明隨林地土壤深度增加，各級(jí)別大孔隙的數(shù)量整體呈顯著下降趨勢(shì)(表3)。林地土壤表層，當(dāng)量孔徑0.5—3.8 mm范圍內(nèi)的大孔隙密度明顯高于底層。林地表層土壤中，當(dāng)量孔徑為2.0—3.8 mm的平均大孔隙密度在135—768個(gè)/m2之間，而底層僅為45—121個(gè)/m2。該孔徑級(jí)別土壤大孔隙的形成可能與植物根系腐爛后形成的孔徑較大的植物根孔有關(guān)［21］。當(dāng)量孔徑為0.5—2.0 mm的大孔隙，在表層土壤中的平均數(shù)量可達(dá)到4522個(gè)/m2。地被物中灌木、草本植物的淺層根系的腐爛與生長(zhǎng)是植物根孔更新形成的主要來(lái)源，因此當(dāng)量孔徑較大的大孔隙多在土壤表層范圍內(nèi)聚集［21］。當(dāng)量孔徑為0.3—0.5 mm的大孔隙密度在相鄰?fù)翆娱g變化不顯著，說(shuō)明數(shù)量最多的低孔徑大孔隙在林地土壤中的空間分布狀態(tài)是較為均勻的。對(duì)各當(dāng)量孔徑的大孔隙密度與對(duì)應(yīng)有機(jī)質(zhì)含量進(jìn)行相關(guān)分析，二者具有極顯著的正相關(guān)性(P＜0.01)(表4)，表明有機(jī)質(zhì)對(duì)大孔隙的形成有重要作用，對(duì)較大孔徑孔隙的作用尤為明顯。土壤有機(jī)質(zhì)具有一定的粘結(jié)力，能夠使相對(duì)松散的土壤顆粒通過(guò)有機(jī)質(zhì)自身的膠結(jié)作用形成團(tuán)粒結(jié)構(gòu)。森林枯枝落葉分解和根系腐爛導(dǎo)致表層土壤有機(jī)質(zhì)含量增加，土壤團(tuán)聚體增多，相應(yīng)地土壤孔隙率增大、孔隙尺度增大、連通性變好。Tsukamota［1］研究發(fā)現(xiàn)單個(gè)大孔隙通常較短，但通過(guò)土壤松散有機(jī)物質(zhì)的節(jié)點(diǎn)連接，可形成大孔隙網(wǎng)絡(luò)。Aubertin［1］研究認(rèn)為此類大孔隙可能至少占到森林土壤體積的35%，但會(huì)隨土層的加深而迅速減小。

表3 大孔隙密度隨土壤深度變化的差異顯著性Table 3 Significance of difference of macropore density in each soil depth

表4 有機(jī)質(zhì)含量與不同當(dāng)量孔徑大孔隙密度的相關(guān)關(guān)系Table 4 Correlation analysis between the number of macropores and organic matter for each pore size

2.3 大孔隙對(duì)飽和導(dǎo)水率的影響

土壤層的飽和導(dǎo)水率是由不同類型不同性質(zhì)的水分運(yùn)移路徑所決定的，而土壤大孔隙是其中一種重要的路徑表現(xiàn)形式［1，22］。對(duì)各土層不同當(dāng)量孔徑的大孔隙密度和土壤飽和導(dǎo)水率進(jìn)行相關(guān)分析(表5)，發(fā)現(xiàn)當(dāng)量孔徑為0.3—0.5 mm的大孔隙密度對(duì)飽和導(dǎo)水率具有顯著的影響(P＜0.05)。由于不同林地不同土層內(nèi)該級(jí)別大孔隙的數(shù)量最大且變化差異較小，因此整體上其對(duì)飽和導(dǎo)水率空間變化的影響程度較低。當(dāng)量孔徑為0.5—3.8 mm的大孔隙密度對(duì)飽和導(dǎo)水率的影響極顯著(P＜0.01)，加之其在不同林分、不同土層中數(shù)量變化較大(表2)，因此其對(duì)飽和導(dǎo)水率的空間變化影響程度較高。大孔隙面積比與飽和導(dǎo)水率顯著相關(guān)(P＜0.05)也反映了大孔隙孔徑與密度分布共同影響土壤水分運(yùn)移。隨孔徑降低，大孔隙密度與飽和導(dǎo)水率的相關(guān)程度逐漸降低，說(shuō)明土壤中水分傳導(dǎo)特性與大孔隙孔徑特征具有密切聯(lián)系。油松林地土壤表層(0—20 cm)豐富的根系石礫分布所形成的較大孔徑的孔隙使其具有較強(qiáng)的導(dǎo)水性能;針闊混交林地底層(40—60 cm)的破碎基巖是其具有較好深層導(dǎo)水能力的重要原因。大孔隙的水分傳導(dǎo)作用也是各土層不同徑級(jí)大孔隙在土壤水分、質(zhì)地和寄生生物等因素影響下相互作用、相互貫通的過(guò)程。

表5 飽和導(dǎo)水率與不同當(dāng)量孔徑大孔隙密度的相關(guān)關(guān)系Table 5 Correlation analysis between the number of macropores and saturated water conductivity for each pore size

進(jìn)一步研究大孔隙對(duì)水分運(yùn)動(dòng)影響，對(duì)當(dāng)量孔徑＞1.5mm大孔隙密度與飽和導(dǎo)水率進(jìn)行回歸分析，發(fā)現(xiàn)孔徑＞1.5mm的大孔隙密度與飽和導(dǎo)水率呈冪函數(shù)關(guān)系，且擬合度較高(圖1)。孔徑＞1.5mm的大孔隙密度僅占總大孔隙密度的5.37%，但決定了飽和導(dǎo)水率84.4%的變異。Beven等［1］認(rèn)為較大孔徑的孔隙可有效排走土壤中的空氣，減小了氣壓對(duì)入滲的抑制;大孔隙也增加了側(cè)向滲透和垂直下滲，從而影響導(dǎo)水率的變化。在考慮大孔隙存在時(shí)，土壤中水分運(yùn)移主要由通過(guò)大孔隙中的水流決定，特別是在土壤接近飽和的情況下，大孔隙中的水流速度遠(yuǎn)大于土壤基質(zhì)流，是非達(dá)西流［1，16］。Peterson［1］等測(cè)量孔徑＞3 mm的某單一孔隙出流時(shí)，發(fā)現(xiàn)其雷諾數(shù)達(dá)到3000，為湍流;有的大孔隙可增加水流通量，從1.7到2.8×10-5m/s，孔隙率卻僅增加了0.002%。林地土壤較大孔徑的孔隙主要分布在表層(表2)，所以更有利于土壤對(duì)降水的吸收與存蓄。

圖1 孔徑＞1.5mm的大孔隙密度與飽和導(dǎo)水率的關(guān)系Fig.1 The relationship between density of soil macropore＞1.5mm and saturated water conductivity

林區(qū)土壤平均大孔隙率在1.6%—13.3%之間。飽和導(dǎo)水率隨著大孔隙率的增加呈冪函數(shù)增大;進(jìn)一步研究也發(fā)現(xiàn)飽和導(dǎo)水率先隨大孔隙率的增加而增大，到一定程度反而下降，呈拋物線型變化，且擬合度也達(dá)到了0.81。但在大孔隙體積所占比例小于5%時(shí)，均表現(xiàn)為飽和導(dǎo)水率隨大孔隙體積比增大而升高的趨勢(shì)(圖2)。分析其原因，可能是大孔隙體積比增大到一定程度時(shí)，決定其排水能力的主要因素轉(zhuǎn)為大孔隙聯(lián)通性和水流型態(tài)，而非其數(shù)量的多少［15］。大孔隙率的微小變化會(huì)引起水流通量的突變，其流態(tài)可能是層流，也可能是界于層流與紊流之間的過(guò)渡流。水流在不規(guī)則大孔隙網(wǎng)絡(luò)中傾向于頸縮效應(yīng)［1］。Beven和Germann研究發(fā)現(xiàn)土壤樣品中連續(xù)且直徑為0.3 mm的大孔隙要比周圍不連續(xù)的直徑為100 mm的孔隙傳導(dǎo)更多的水分。土壤中的垂直大孔隙受到植物根系、動(dòng)物活動(dòng)和土壤結(jié)構(gòu)等因素的影響，其間往往不是聯(lián)通的，存在著斷孔和死孔。雖存在土壤有機(jī)物質(zhì)的節(jié)點(diǎn)連接［8］，但其作用有限，并且隨土層加深大孔隙密度也顯著減小(表2)，從而影響著對(duì)水分的傳導(dǎo)。

圖2 大孔隙率與飽和導(dǎo)水率的關(guān)系Fig.2 The relationship between macroporosity and saturated water conductivity

3 結(jié)論與討論

本研究得出秦嶺南坡中山帶水源林區(qū)土壤大孔隙當(dāng)量孔徑集中于0.3—3.8 mm，與以往定義范圍基本一致，但孔徑范圍上限較大，這可能與該區(qū)石礫含量較大并廣泛分布著基巖碎屑和破碎巖石有關(guān)，且林區(qū)植被的腐根對(duì)較大孔徑的大孔隙形成也起重要作用［23］。根系腐爛不僅有助于提高土壤有機(jī)質(zhì)，促進(jìn)團(tuán)聚體的形成;也為新根的伸延提供了有效通道。多數(shù)情況，石礫與植物根系交織在一起，共同影響大孔隙的發(fā)育以及土壤水分運(yùn)動(dòng)。土壤大孔隙的形成因素較為復(fù)雜，若需獲得更為全面的了解，則需要定量區(qū)分各個(gè)因子的作用。

＞0.3 mm的各當(dāng)量孔徑大孔隙密度與飽和導(dǎo)水率均顯著相關(guān)，這與王偉等［18］研究結(jié)果基本一致。石輝等［15］在研究岷江上游森林土壤大孔隙特征時(shí)，發(fā)現(xiàn)孔徑大于1.4 mm的孔隙數(shù)量對(duì)飽和導(dǎo)水率影響最大，可以決定飽和導(dǎo)水率69%的變異。本研究中，當(dāng)量孔徑＞1.5 mm的大孔隙密度較小，但決定了飽和導(dǎo)水率84.4%的變異?？讖捷^大(當(dāng)量孔徑為1.5—3.8 mm)的大孔隙，在實(shí)驗(yàn)林區(qū)不同類型林地土壤中廣泛分布，加之該類型大孔隙的水力傳導(dǎo)性能較強(qiáng)，對(duì)飽和導(dǎo)水率的變化影響程度也較高。一定條件下大孔隙的連通性對(duì)水分的入滲起主導(dǎo)作用［24］。但因技術(shù)原因，本研究未考慮孔隙聯(lián)通性，在今后研究中需加強(qiáng)技術(shù)研發(fā)，予以綜合考慮。

總體上，研究區(qū)針闊混交林和落葉闊葉林土壤大孔隙平均密度和面積比較大，因此，針闊混交林和落葉闊葉林較利于改良土壤滲透性，具有較好的涵養(yǎng)水源和調(diào)節(jié)徑流的功能。影響流域水分運(yùn)動(dòng)的因素很多，除植被類型和土壤外，也包括降雨、地質(zhì)和地形等［25］。以往研究中，很少討論當(dāng)降雨強(qiáng)度處于那個(gè)范圍時(shí)，更有利于大孔隙傳導(dǎo)水分;是否大孔隙度越大，導(dǎo)水性越好，還是存在上下限。此外，未來(lái)還需進(jìn)一步討論大孔隙的變異性和有效性，為水源林區(qū)水文生態(tài)理論的完善和模型的構(gòu)建提供參考。

［1］Beven K，Germann P.Macropores and water flow in soils.Water Resources Research，1982，18(5):1311-1325.

［2］Noguchi S，Tsuboyama Y，Sidle R C，Hosoda I.Morphological characteristics of macropores and the distribution of preferential flow pathways in a forested slope segment.Soil Science Society of America Journal，1999，63(5):1413-1423.

［3］Wilson G V，Jardine P M，Luzmoore R J，Jones J R.Hydrology of a forested hillslope during storm events.Geoderma，1990，46(1/3):119-138.

［4］Shipitalo M J，Dick W A，Edwards W M.Conservation tillage and macropore factors that affect water movement and the fate of chemicals.Soil and Tillage Research，2000，53(3/4):167-183.

［5］Weiler M，Naef F.Simulating surface and subsurface initiation of macropore flow.Journal of Hydrology， 2003， 273(1/4):139-154.

［6］Cameira M R，F(xiàn)ernando R M，Pereira L S.Soil macropore dynamics affected by tillage and irrigation for a silty loam alluvial soil in southern Portugal.Soil and Tillage Research，2003，70(2):131-140.

［7］Luxmoore R J.Micro-，meso-，and macroporosity of soil.Soil Science Society of America Journal，1981，45(3):671-672.

［8］Shi Z J，Wang Y H，Xu L H，Yu P T，Xiong W，Xu D P.Soil macropore characteristics under typical vegetations in Liupan Mountains.Chinese Journal of Applied Ecology，2007，18(12):2675-2680.

［9］Liu W，Qu Z Q，Ying P F.Soil macropore and its studying methodology.Chinese Journal of Applied Ecology，2001，12(6):465-468.

［10］Bouma J.Influence of soil macroporosity on environmental quality.Advances in Agronomy，1991，46:1-39.

［11］Zehe E，F(xiàn)lühler H.Preferential transport of isoproturon at a plot scale and a field scale tile-drained site.Journal of Hydrology，2001，247(1/2):100-115.

［12］Li J Z，Pei T F，Niu L H，Li X Y.Simulation and model of interflow on hillslope of forest catchment.Scientia Silvae Sinicae，1999，35(4):2-8.

［13］Zhang S L，Lei R D，Lü Y L，Ma Y.Water-balance of forest ecosystem in Huoditang area of Qinling mountain.Bulletin of Soil and Water Conservation，2000，20(6):18-22.

［14］Bronsterta A，Plate E J.Modelling of runoff generation and soil moisture dynamics for hillslopes and micro-catchments.Journal of Hydrology，1997，198(1/4):177-195.

［15］Shi H，Chen F Q，Liu S R.Macropores properties of forest soil and its influence on water effluent in the upper reaches of Minjiang Rive.Acta Ecologica Sinica，2005，25(3):507-512.

［16］Li W L，Jin C J，Wang A Z，Pei T F，Guan D X.Amount and vertical distribution of macropores in forest soils in Changbai Mountains.Chinese Journal of Applied Ecology，2007，18(6):1213-1218.

［17］Radulovich R，Solorzano E，Sollins P.Soil macropore size distribution from water breakthrough curves.Soil Science Society of America Journal，1989，53(2):556-559.

［18］Wang W，Zhang H J，Cheng J H，Wu Y H，Du T C，Wang Y.Macropore characteristics and its relationships with the preferential flow in broadleaved forest soils of simian mountains.Chinese Journal of Applied Ecology，2010，21(5):1217-1223.

［19］Zhu Z X.Pedology.Beijing:Agriculture Press，1983.

［20］Wang D L，Yin C Q.Functions of root channels in the soil system.Acta Ecologica Sinica，2000，20(5):869-874.

［21］Gish T J，Gimenez D，Rawls W J.Impact of roots on ground water quality.Plant and Soil，1998，200(1):47-54.

［22］Cheng J H.Research on the model of Preferential Flow in the Woodland of Granite Region Three Gorges of Yangtze River［D］.Beijing:Beijing Forestry University Doctoral Dissertation，2005.

［23］Shi H，Liu S R.The Macroporosity properties of forest soil and its eco-hydrological significance.Journal of Mountain Science，2005，23(5):533-539.

［24］Wu J Q，Zhang J F，Gao R.Physical simulation experiments of effects of macropores on soil water infiltration characteristics.Transactions of the CSAE，2009，25(10):13-18.

［25］Zhang Z Q，Wang L X，Yu X X，Klhgaoefr E.Impacts of forest vegetation on runoff generation mechanisms:a review.Journal of Natural Resources，2001，16(1):79-84.

參考文獻(xiàn):

［8］時(shí)忠杰，王彥輝，徐麗宏，于澎濤，熊偉，徐大平.六盤山典型植被下土壤大孔隙特征.應(yīng)用生態(tài)學(xué)報(bào)，2007，18(12):2675-2680.

［9］劉偉，區(qū)自清，應(yīng)佩峰.土壤大孔隙及其研究方法.應(yīng)用生態(tài)學(xué)報(bào)，2001，12(6):465-468.

［12］李金中，裴鐵璠，牛麗華，李曉宴.森林流域坡地壤中流模型與模擬研究.林業(yè)科學(xué)，1999，35(4):2-8.

［13］張勝利，雷瑞德，呂瑜良，馬玉.秦嶺火地塘林區(qū)森林生態(tài)系統(tǒng)水量平衡研究.水土保持通報(bào)，2000，20(6):18-22.

［15］石輝，陳鳳琴，劉世榮.岷江上游森林土壤大孔隙特征及其對(duì)水分出流速率的影響.生態(tài)學(xué)報(bào)，2005，25(3):507-512.

［16］李偉莉，金昌杰，王安志，裴鐵璠，關(guān)德新.長(zhǎng)白山主要類型森林土壤大孔隙數(shù)量與垂直分布規(guī)律.應(yīng)用生態(tài)學(xué)報(bào)，2007，18(10):2179-2154.

［18］王偉，張洪江，程金花，吳煜禾，杜士才，王冉.四面山闊葉林土壤大孔隙特征與優(yōu)先流的關(guān)系.應(yīng)用生態(tài)學(xué)報(bào)，2010，21(5):1217-1223.

［19］朱祖祥.土壤學(xué)(上冊(cè)).北京:農(nóng)業(yè)出版社，1983.

［20］王大力，尹澄清.植物根孔在土壤生態(tài)系統(tǒng)中的功能.生態(tài)學(xué)報(bào)，2000，20(5):869-874.

［22］程金花.長(zhǎng)江三峽花崗巖區(qū)林地坡面優(yōu)先流模型研究［D］.北京:北京林業(yè)大學(xué)，2005.

［23］石輝，劉世榮.森林土壤大孔隙特征及其生態(tài)水文學(xué)意義.山地學(xué)報(bào)，2005，23(5):533-539.

［24］吳繼強(qiáng)，張建豐，高瑞.大孔隙對(duì)土壤水分入滲特性影響的物理模擬試驗(yàn).農(nóng)業(yè)工程學(xué)報(bào)，2009，25(10):15-17.

［25］張志強(qiáng)，王禮先，余新曉，Klhgaoefr E.森林植被影響徑流形成機(jī)制研究進(jìn)展.自然資源學(xué)報(bào)，2001，16(1):79-84.