張東旭，程金花，王 偉，張洪江，劉 濤

?

基于O-ring統(tǒng)計的公路邊坡土壤優(yōu)先流路徑分布分析

張東旭1，程金花1※，王 偉2，張洪江1，劉 濤2

（1. 北京林業(yè)大學(xué)水土保持學(xué)院，北京 100083； 2. 交通運輸部科學(xué)研究院，北京 100029）

引入優(yōu)先流理論分析氯鹽污染物在邊坡土壤中的遷移機制，對揭示公路工業(yè)鹽污染物遷移規(guī)律和保護路域環(huán)境具有重要意義。以鶴大高速低路基公路邊坡土壤為例，采用染色示蹤試驗對土壤優(yōu)先流路徑的位置進行確定，并應(yīng)用景觀生態(tài)學(xué)中的空間點格局分析方法O-ring統(tǒng)計，對不同影響半徑土壤優(yōu)先流路徑水平分布特征及其空間關(guān)聯(lián)性進行分析。結(jié)果表明：邊坡土壤優(yōu)先流路徑發(fā)育深度可達40 cm；相同影響半徑優(yōu)先流路徑的數(shù)量隨著土層深度的增加而減少，同一土層優(yōu)先流路徑的數(shù)量隨影響半徑的增大逐漸減少；隨著影響半徑的增加，聚集分布趨于顯著；不同影響半徑的優(yōu)先流路徑之間的關(guān)聯(lián)性相對獨立。O-ring統(tǒng)計，能從統(tǒng)計學(xué)角度量化分析土壤優(yōu)先流路徑的分布特征，豐富了優(yōu)先流的研究方法。

統(tǒng)計；土壤；滲透；空間點格局分析；公路邊坡；優(yōu)先流；路徑；分布特征

0 引 言

中國東北地區(qū)冬季降雪量較大，受經(jīng)濟條件限制，目前多數(shù)公路長期使用工業(yè)鹽（氯鹽）進行除雪。東北平原微丘區(qū)公路多采用低路基設(shè)計方式，其路基高度一般在2 m以內(nèi)，采用路面散排結(jié)合放緩邊坡的設(shè)計方式[1]。融雪后，工業(yè)鹽通過公路排水系統(tǒng)進入邊坡表層土壤，導(dǎo)致土壤有機質(zhì)含量下降和土壤板結(jié)，綠化植物發(fā)生嚴(yán)重的鹽害，繼而危害路基的穩(wěn)定[2-3]。邊坡土壤層結(jié)構(gòu)在公路修筑過程中受到擾動，表現(xiàn)出非均質(zhì)性質(zhì)，水分和溶質(zhì)在通過時以非平衡流方式運動[4]，即為優(yōu)先流。引入優(yōu)先流分析工業(yè)鹽污染物在邊坡土壤中的遷移機制，對揭示公路工業(yè)鹽污染物遷移規(guī)律和保護路域環(huán)境具有重要意義[5-6]。

目前關(guān)于優(yōu)先流的研究主要集中在林地、農(nóng)地等自然立地條件的土壤[7-8]。近年來，隨著開發(fā)建設(shè)項目日漸增多，對地下水的污染和對環(huán)境的破壞日益加重，一些學(xué)者開始研究垃圾填埋場和礦區(qū)的優(yōu)先流[9-10]，但公路邊坡中的優(yōu)先流還鮮有研究[11]。優(yōu)先流路徑可定義為水分和溶質(zhì)在土壤中發(fā)生優(yōu)先傳導(dǎo)過程的通道[12]。對于優(yōu)先流路徑的研究，最直觀的方法就是通過染色示蹤試驗對優(yōu)先流路徑進行識別和標(biāo)記。近年來，有少數(shù)學(xué)者研究了優(yōu)先流路徑的數(shù)量形態(tài)特征和空間位置分布[13]。魏虎偉等[14]采用水分穿透曲線和Poiseulle方程分析了土壤大孔隙路徑的數(shù)量、形態(tài)特征，在一定程度上揭示了優(yōu)先流路徑的數(shù)量特征，但沒有對其分布狀態(tài)進行深入研究，存在一定的局限性。呂文星等[15]采用Ripley’s K()函數(shù)的分析方法，對重慶四面山林地10～20 cm土壤優(yōu)先流路徑分布特征進行研究，得出優(yōu)先流路徑大多為均勻和聚集2種分布狀態(tài)，但該方法在分析過程中無法消除中小尺度累積效應(yīng)對大尺度的影響，存在一定的偏差，不能準(zhǔn)確反映出優(yōu)先流路徑的分布狀態(tài)[16]。

本文以鶴大高速（G11）低路基公路邊坡土壤染色圖像數(shù)據(jù)為基礎(chǔ)，應(yīng)用景觀生態(tài)學(xué)中空間點格局分析方法O-ring統(tǒng)計，利用()分布函數(shù)展現(xiàn)不同土層不同影響半徑優(yōu)先流路徑的分布，來揭示邊坡土壤中優(yōu)先流路徑水平分布特征及其相互之間的空間關(guān)聯(lián)性，同時豐富土壤優(yōu)先流的研究方法，使景觀生態(tài)學(xué)中O-ring統(tǒng)計這種基于統(tǒng)計學(xué)的方法更好地應(yīng)用在小尺度的土壤優(yōu)先流路徑空間分布的研究上。

1 研究區(qū)概況

本研究區(qū)位于鶴崗至大連高速公路（G11）丹東至通化段寬甸互通處（124°43'39.70"E，40°43'58.07"N），此處路段為低路基（路基高度≤2 m）公路，且路邊無排水溝。該區(qū)地處遼東斷塊山地丘陵區(qū)，地勢自北向南、自西北向東南降低，呈階梯狀，平均海拔為500 m，屬東亞季風(fēng)區(qū)，溫帶濕潤大陸性季風(fēng)氣候。年平均氣溫6.5 ℃。年平均降水量在1 000～1 200 mm之間，降雨多集中在6～8月，占全年降雨量的2/3，無霜期129 d。土壤發(fā)育較好，多為褐土。在研究區(qū)的公路邊坡上，選擇3個相鄰的點作為試驗樣地，樣地基本情況介紹見表1。

表1 寬甸互通邊坡土壤性質(zhì)

2 研究方法

2.1 染色示蹤試驗

試驗樣地布設(shè)在低路基公路邊坡上，邊坡坡度較緩，為15°，坡面上種植紫穗槐和高羊茅等灌草植物。在坡面上選取較為平緩的區(qū)域作為觀測區(qū)，并選取2塊相鄰的區(qū)域進行平行試驗。首先清除觀測區(qū)域內(nèi)的枯落物，平整土壤表面，將長70 cm、寬70 cm、高30 cm、厚0.5 cm的自制鐵框鉛直埋入邊坡土壤中，埋入深度為20 cm。將距鐵框內(nèi)壁5 cm范圍內(nèi)的土壤壓實，防止染料隨縫隙下滲影響試驗結(jié)果。根據(jù)當(dāng)?shù)厮鸭慕笛┵Y料和路面產(chǎn)流特性，估算出在大雪（降雪量>5 mm）時產(chǎn)生的融雪水徑流量約為10 L。配制10 L濃度為4 g/L的亮藍溶液作為示蹤溶液，通過積水滲透方式[9]，使用自制恒流泵以150 mL/min的流量，將預(yù)制的亮藍溶液自邊坡頂部均勻的輸入鐵框內(nèi)部，鐵框內(nèi)溶液入滲較為均勻，沒有產(chǎn)生積水，然后用塑料布遮蓋鐵框上部，防止外部條件變化對試驗產(chǎn)生影響。染色示蹤試驗完成24 h后除去覆蓋的塑料布和埋置的鐵框挖掘邊坡剖面，用標(biāo)尺標(biāo)注剖面的長度與寬度，使用高像素數(shù)碼相機采集垂直染色剖面圖像，各垂直剖面拍攝照片2張，剖面旁放置標(biāo)準(zhǔn)灰階比色卡，用以校正圖像色彩，減少明暗不均所帶來的誤差。在染色范圍內(nèi)取其中50 cm×50 cm面積，以10 cm深度為標(biāo)準(zhǔn)，分層挖掘水平剖面并拍攝照片，根據(jù)染色劑分布狀況將剖面分為染色區(qū)和未染色區(qū)，圖1為染色試驗示意圖。

2.2 圖像處理

使用ARCMAP 10.2軟件對采集到的原始染色圖像進行平面幾何校正，水平剖面染色圖像大小為50 cm× 50 cm（500×500像素）。使用Adobe Photoshop CS3軟件中圖像調(diào)整中顏色替換功能，先參照原始照片，用吸管工具對照片中染色土壤的顏色采樣，確定染色區(qū)域。通過設(shè)置色相、飽和度將染色土壤顏色替換為黑色，未染色區(qū)土壤為白色。圖像數(shù)值化轉(zhuǎn)換可采用Image-Pro Plus 6.0軟件中的Bitmap功能，黑色區(qū)域值為255，白色區(qū)域為0，并以Excel格式輸出。染色剖面圖像處理步驟見圖2。

2.3 優(yōu)先流路徑影響半徑及其位置信息提取

水分通過土壤中的優(yōu)先流路徑快速下滲到深層，由于受土壤含水率變化、土壤顆粒組成等因素的影響，其染色的范圍要大于土壤中實際的優(yōu)先流路徑，可將染色所在區(qū)域的半徑看作是土壤優(yōu)先流路徑的“影響半徑”[17]。圖3為優(yōu)先流路徑信息的提取過程，具體操作如下：將邊坡土壤各層水平剖面視為一個整體，根據(jù)染色影響區(qū)域的空間拓?fù)潢P(guān)系，利用形態(tài)學(xué)計算方法，將染色影響區(qū)域劃分成多個獨立的斑塊，染色影響區(qū)的面積即為斑塊所占像素點的個數(shù)。把這些斑塊近似看作圓，根據(jù)染色面積計算其影響半徑。可將影響半徑劃分為1.0～2.5、>2.5～5.0、>5.0～10.0和≥10.0 mm 4個等級。用Image ProPlus 6.0軟件進行提取，得到優(yōu)先流路徑的影響半徑及其重心位置坐標(biāo)。根據(jù)坐標(biāo)便可進行優(yōu)先流路徑空間分布格局和關(guān)聯(lián)性的分析。

2.4 O-ring統(tǒng)計

O-ring統(tǒng)計分析方法是在Ripley’s K函數(shù)和Mark函數(shù)的基礎(chǔ)上演變而來的，其方法是使用圓環(huán)代替圓作為統(tǒng)計工具，可以有效消除Ripley’s K分析中小尺度累積效應(yīng)對大尺度造成的影響，理論上提高了空間格局分析的精確性[18]。目前，該方法通常應(yīng)用在動植物群落這種大尺度的空間分布格局及其關(guān)聯(lián)性的研究上[19]。

O-ring統(tǒng)計分為單變量統(tǒng)計和雙變量統(tǒng)計[20]。式（1）為單變量統(tǒng)計，用于分析某一影響半徑的優(yōu)先流路徑的空間分布格局，式（2）為雙變量統(tǒng)計，用于分析2種不同影響半徑優(yōu)先流路徑之間的空間關(guān)聯(lián)性。

式中為半徑，mm；為圓環(huán)寬度，mm；11()是水平分布格局函數(shù)；為觀測區(qū)域內(nèi)某一影響半徑的優(yōu)先流路徑的數(shù)量；R()表示該影響半徑的優(yōu)先流路徑中的第點為圓心，為半徑，為寬的圓環(huán)；Points[R()]表示圓環(huán)內(nèi)該影響半徑的優(yōu)先流路徑的數(shù)量；Area[R()]表示圓環(huán)的面積，mm2。

式中O12()是基于個體與個體間距離的關(guān)聯(lián)性函數(shù)；1為觀測區(qū)域內(nèi)較小影響半徑的優(yōu)先流路徑的數(shù)量；R1,i()表示該影響半徑優(yōu)先流路徑中的第點為圓心；Points2[R1,i()]表示圓環(huán)內(nèi)較大影響半徑的優(yōu)先流路徑的數(shù)量；Area[R1,i()]表示圓環(huán)的面積，mm2。

使用生態(tài)學(xué)軟件Programita 2010對空間格局?jǐn)?shù)據(jù)進行分析。以樣地邊長的1/5(100 mm)作為研究尺度，1 mm作為步長，通過19次Monte Carlo模擬計算獲取95%的置信區(qū)間，即上下包跡線。在某一空間尺度上，對于單變量統(tǒng)計，若11()值大于上包跡線，呈聚集分布；落入上下包跡線內(nèi)，呈隨機分布；小于下包跡線，呈均勻分布；對于雙變量統(tǒng)計，若12()值大于上包跡線，二者呈空間上正關(guān)聯(lián)；在上下包跡線之內(nèi)，二者相互獨立；小于下包跡線，二者呈空間上負(fù)關(guān)聯(lián)。

3 結(jié)果與分析

表層0～10 cm土壤染色區(qū)域為80%以上（圖4），且染色較為均勻，這主要是由于表層土壤中基質(zhì)流作用要遠(yuǎn)遠(yuǎn)大于優(yōu)先流作用，染色區(qū)域并非完全因優(yōu)先流作用而染色[21]，故把該層看作是基質(zhì)流作用區(qū)，不對表層0～10 cm土壤進行優(yōu)先流的分析。圖5為公路邊坡豎直剖面的染色圖像，3個平行試驗樣地的土壤染色深度均為40 cm，選擇1#樣地為例展示其土壤染色剖面優(yōu)先流路徑的分布特征及其空間關(guān)聯(lián)性。

3.1 邊坡土壤優(yōu)先流路徑數(shù)量及位置

通過圖像提取技術(shù)從處理后的染色圖像中提取出公路邊坡土壤優(yōu)先流路徑的位置信息。經(jīng)統(tǒng)計，10～20、>20～30和>30～40 cm土層不同影響半徑優(yōu)先流路徑數(shù)量見表2。10～20、>20～30、>30～40 cm土層優(yōu)先流路徑總量分別為713、510、341個，優(yōu)先流路徑總數(shù)和相同影響半徑優(yōu)先流路徑的數(shù)量均隨著土層深度的增加而減少。受到人為擾動和土壤顆粒組成的影響，表層土壤容重較小，土質(zhì)疏松，大孔隙較多[22]，因此優(yōu)先流路徑數(shù)量較多。而深層土壤容重較大，土壤緊實，受到的擾動也少，故優(yōu)先流路徑數(shù)量較少。在同一土層中，1.0～10.0 mm影響半徑的優(yōu)先流路徑分布較多，≥10.0 mm的分布極少，小影響半徑的優(yōu)先流路徑數(shù)量要遠(yuǎn)遠(yuǎn)多于大影響半徑優(yōu)先流路徑的數(shù)量，這是由于邊坡土壤較為緊實，很難形成較大的孔隙。這與王彬儼等[8]采用土壤水分穿透曲線法和Zhang等[23]采用CT掃描測定的農(nóng)地土壤優(yōu)先流路徑即大孔隙數(shù)量變化趨勢的結(jié)果相一致。

表2 各土層優(yōu)先流路徑數(shù)量

3.2 邊坡土壤優(yōu)先流路徑分布格局

表3為各層土壤水平剖面上，各影響半徑的優(yōu)先流路徑在不同尺度上的分布狀態(tài)。圖6為公路邊坡土壤4種影響半徑的優(yōu)先流路徑空間分布格局，11()函數(shù)曲線的上下波動表示優(yōu)先流路徑空間分布的變化狀態(tài)。由于>30～40 cm土層≥10.0 mm的優(yōu)先流路徑數(shù)量較少，曲線不能真實地反映出分布狀態(tài)。結(jié)合表3和圖6可以看出，相同土層不同影響半徑優(yōu)先流路徑分布格局曲線的最大值（峰值）隨著影響半徑的增大，其對應(yīng)的距離尺度也會相應(yīng)增加，即波峰的位置逐漸向較大尺度的方向移動。較小影響半徑（1.0≤＜5.0 mm）的優(yōu)先流路徑以隨機分布為主，較大影響半徑（≥5.0 mm）以聚集分布為主，即影響半徑越大，聚集分布越明顯。隨著距離尺度的增大，優(yōu)先流路徑分布格局曲線的呈現(xiàn)出先增大后減小最后逐漸趨于穩(wěn)定的趨勢，即優(yōu)先路徑的分布由聚集分布逐漸向隨機分布變化，且聚集狀態(tài)的峰值一般出現(xiàn)在較小尺度上（0～40 mm之間）。>20～30 cm土層優(yōu)先路徑變化最為復(fù)雜，表明該土壤層是優(yōu)先流形態(tài)變化的過渡區(qū)域。

表3 各影響半徑優(yōu)先流路徑的分布狀態(tài)

圖6中，從不同土層相同影響半徑()的優(yōu)先流路徑分布特征來看，1.0≤＜2.5 mm和2.5≤＜5.0 mm區(qū)間的優(yōu)先流路徑分布格局曲線差異較大，在0～40 mm尺度上有較大幅度的波動，呈現(xiàn)出先升高后降低的趨勢，且主要表現(xiàn)為隨機分布；5.0≤＜10.0 mm區(qū)間的優(yōu)先流路徑分布格局曲線趨勢較為一致，在0～20 mm尺度上呈現(xiàn)均勻分布，隨后由隨機分布過渡到聚集分布，最后趨于隨機分布的變化規(guī)律；≥10.0 mm的優(yōu)先流路徑分布格局曲線表現(xiàn)為先迅速增大而后下降的趨勢，并伴有小幅波動，在20～100 mm尺度上均呈現(xiàn)出聚集分布。在相同影響半徑下，隨土層深度的增加，其分布格局曲線的波峰出現(xiàn)的尺度位置呈現(xiàn)不規(guī)則的變化，優(yōu)先流路徑的分布狀態(tài)也有較大的變化，說明優(yōu)先流路徑在土壤垂直剖面上的連通性較差。陳曉冰等[24]在三峽庫區(qū)紫色砂巖土農(nóng)地的研究中認(rèn)為優(yōu)先流路徑在土壤豎直剖面上呈現(xiàn)出良好的連通性，與本研究的結(jié)果存在一定的差異。這是由于公路邊坡含有大量的礫石，水流在向深層土壤流動的過程中，會經(jīng)過石塊的切割阻擋，導(dǎo)致水流的運動方向發(fā)生了偏移，使實際流動軌跡的彎曲度加大，造成優(yōu)先流路徑出現(xiàn)先分離后合并的現(xiàn)象，即所謂的延遲效應(yīng)[25]。由于本試驗各土層之間的距離較大，為10 cm，也會導(dǎo)致各水平剖面間表現(xiàn)出較差的連通性。植物根系和生物活動在邊坡土壤中較少，也很難形成連通性較好的優(yōu)先流路徑。土壤組成的不均會導(dǎo)致土壤分層，水分隨優(yōu)先流路徑下滲的過程中也會出現(xiàn)水平運動的現(xiàn)象，這說明公路邊坡中優(yōu)先流路徑的空間異質(zhì)性規(guī)律[26]。

3.3 邊坡土壤不同影響半徑優(yōu)先流路徑的空間關(guān)聯(lián)性

以往對農(nóng)林地的研究表明，相同影響半徑優(yōu)先流路徑的分布不僅與土壤的風(fēng)化、崩解、淋溶等過程以及機械組成、土壤養(yǎng)分、生物活動[27]等內(nèi)外因素有關(guān)，還受到其他影響半徑優(yōu)先流路徑的影響[24]，故本研究對優(yōu)先流路徑的空間格局關(guān)聯(lián)性進行了分析。各土層得到的結(jié)果較為一致，以20～30 cm土層為例，不同優(yōu)先流路徑之間的空間關(guān)聯(lián)性見圖7?？臻g關(guān)聯(lián)性函數(shù)12()曲線均在上下包跡線之間，表明各土層不同影響半徑優(yōu)先流路徑之間大多呈相互獨立的關(guān)系，即空間關(guān)聯(lián)性較差，與農(nóng)林地等自然立地的研究存在一定差異。其原因是由公路邊坡自身的結(jié)構(gòu)特性所決定的。在農(nóng)地中，土壤大孔隙會影響土壤結(jié)構(gòu)的穩(wěn)定性，大量小的優(yōu)先流路徑會塌陷成為大的優(yōu)先流路徑[28]。而公路邊坡土壤結(jié)構(gòu)較為緊實，且礫石含量較高，礫石之間的大孔隙穩(wěn)定性較好[29]，不易發(fā)生崩塌，較大優(yōu)先流路徑的形成很難受到較小優(yōu)先流路徑的影響。公路邊坡是由松散填料刷坡形成的邊坡層，以及上面覆蓋種植土組成，邊坡土壤中大孔隙主要是由于碎石土細(xì)粒填充不足而形成的，大小不一的礫石和土粒分布不均，導(dǎo)致形成的優(yōu)先流路徑分布不均，隨機性較強，不同影響半徑的優(yōu)先流路徑關(guān)聯(lián)性也相對獨立。植物根系和蟲洞是形成聯(lián)通性較好優(yōu)先流路徑的另一重要因子[30]。根據(jù)現(xiàn)場土壤剖面的調(diào)查，邊坡中植物根系多集中在0～10 cm表層土壤中，深層土壤中植物根系和生物活動的痕跡幾乎沒有，故本文不再對生物因子進行討論。

4 結(jié) 論

本文對鶴大高速寬甸互通處公路邊坡土壤進行染色示蹤試驗，并采用O-ring統(tǒng)計方法對優(yōu)先流路徑的空間分布格局和空間關(guān)聯(lián)性進行了分析，得到結(jié)論如下

1）通過分層挖掘土壤剖面，發(fā)現(xiàn)邊坡土壤優(yōu)先流路徑發(fā)育深度可達40 cm；相同影響半徑優(yōu)先流路徑的數(shù)量隨著土層深度的增加而減少，同一土層優(yōu)先流路徑的數(shù)量隨影響半徑的增大逐漸減少。

2）在所研究尺度（0～100 mm）下，公路邊坡土壤優(yōu)先流路徑的空間分布狀態(tài)表現(xiàn)為：10～40 cm土層中，較小影響半徑（1.0≤＜5.0 mm）的優(yōu)先流路徑以隨機分布為主，較大影響半徑（≥5.0 mm）以聚集分布為主。隨著影響半徑的增加，聚集分布趨于顯著。不同土層優(yōu)先流路徑的連通性較差。

3）對于優(yōu)先流路徑的空間關(guān)聯(lián)性分析顯示，各土層不同影響半徑優(yōu)先流路徑相互之間的關(guān)聯(lián)性相對獨立，這與公路邊坡土壤的形成過程和自身結(jié)構(gòu)有關(guān)。邊坡土壤中大孔隙主要是由于碎石土細(xì)粒填充不足而形成的，有著較強的隨機性和較好的穩(wěn)定性。

[1] 王偉，陳宗偉，崔慧珊，等. 平微區(qū)高速公路低路基設(shè)計及其水土保持效果初探[J]. 交通建設(shè)與管理，2011(9)： 88－89.

[2] 陳宗偉，張立塔，劉濤，等. 氯鹽型融雪劑的環(huán)境影響及其替代材料研發(fā)進展[J]. 生態(tài)科學(xué)，2011，30(5)：556－561.

Chen Zongwei, Zhang Lita, Liu Tao, et al. Environmental impacts of chloride-based deicing salt and research progress of its alternative materials[J]. Ecological Science, 2011, 30(5): 556－561. (in Chinese with English abstract)

[3] Fay L, Shi X M. Environmental impacts of chemicals for snow and ice control: state of the knowledge[J]. Water, Air, & Soil Pollution, 2012, 223(5): 2751－2770.

[4] Jarvis N, Etana A, Stagnitti F. Water repellency, near-saturated infiltration and preferential solute transport in a macroporous clay soil[J]. Geoderma, 2008, 143(3/4): 223－230.

[5] 張千千，李向全，王效科，等. 城市路面降雨徑流污染特征及源解析的研究進展[J]. 生態(tài)環(huán)境學(xué)報，2014，23(2)：352－358.

Zhang Qianqian, Li Xiangquan, Wang Xiaoke, et al. Research advance in the characterization and Source apportionment of pollutants in urban roadway runoff [J]. Ecology and Environmental Sciences, 2014, 23(2): 352－358. (in Chinese with English abstract)

[6] Chang P G. Chloride characterization from pavement runoff using automated samplers and specific conductivity sensors at three Eastern Massachusetts locations[D]. Masters Projects, 2009.

[7] 張財寶，劉目興，王文德，等. 三峽庫區(qū)森林土壤優(yōu)先流染色特征[J]. 水土保持學(xué)報，2012，26(2)：80－84.

Zhang Caibao, Liu Muxing, Wang Wende, et al. Dyeing characteristic of preferential flow under forestland in the Three Gorges reservoir area[J]. Journal of Soil and Water Conservation, 2012, 26(2): 80－84. (in Chinese with English abstract)

[8] 王彬儼，程金花，張洪江，等. 北京昌平區(qū)農(nóng)地土壤大孔隙特征[J]. 水土保持學(xué)報，2012，26(3)：189－193.

Wang Binyan, Cheng Jinhua, Zhang Hongjiang, et al. Macropores properties of agricultural land a Changping district in Beijing [J]. Journal of Soil and Water Conservation, 2012, 26(3): 189－193. (in Chinese with English abstract)

[9] 邱戰(zhàn)洪，何春木，朱兵見，等. 不同降雨模式下山谷型垃圾填埋場水分運移及其穩(wěn)定性研究[J]. 巖土力學(xué)，2012，33(10)：3151－3155.

Qiu Zhanhong, He Chunmu, Zhu Bingjian, et al. Investigations of water transport in valley-type MSW landfills and their stabilities subjected to various rainfall patterns[J]. Rock and Soil Mechanics, 2012, 33(10): 3151－3155. (in Chinese with English abstract)

[10] Zhen Qing, Ma Wenmei, Li Mingming, et al. Effects of vegetation and physicochemical properties on solute transport in reclaimed soil at an opencast coalmine site on the Loess Plateau, China[J]. Catena, 2015, 133(3): 403－411.

[11] 官琦，徐則民，田林.植被發(fā)育玄武巖斜坡土體大孔隙尺寸及其主要影響因素[J]. 應(yīng)用生態(tài)學(xué)報，2013，24(10)：2888－2896. Guan Qi, Xu Zemin, Tian Lin. Sizes of soil macropores and related main affecting factors on a vegetated basalt slope[J]. Chinese Journal of Applied Ecology, 2013, 24(10): 2888－2896. (in Chinese with English abstract)

[12] Isensee A R, Helling C S, Gish T J, et al. Groundwater residues of atrazine, alachlor and cyanazine under no tillage practices[J]. Chemosphere, 1988, 17(1): 165－174.

[13] Weiler M, Flühler H. Inferring flow types from dye patterns in macroporous soils[J]. Geoderma, 2004, 120(1/2): 137－153.

[14] 魏虎偉，程金花，張洪江，等. 四面山2種林地大孔隙特征與優(yōu)先流關(guān)系研究[J]. 水土保持學(xué)報，2012，28(4)：263－268.

Wei Huwei, Cheng Jinhua, Zhang Hongjiang, et al. Relationships between soil macropores properties and preferential flow about two kind of woodlands in simian mountains[J]. Journal of Soil and Water Conservation, 2012, 28(4): 263－268. (in Chinese with English abstract)

[15] 呂文星，張洪江，吳煜禾，等. 基于點格局分析的林地表層土壤優(yōu)先路徑水平分布特征[J]. 水土保持學(xué)報，2012，26(6)：68－74.

Lü Wenxing, Zhang Hongjiang, Wu Yuhe, et al. Horizontal distribution characteristics of the preferential flow paths in surface soil of the forest land based on the spatial point pattern analysis[J]. Journal of Soil and Water Conservation, 2012, 26(6): 68－74. (in Chinese with English abstract)

[16] 董靈波，劉兆剛，張博，等. 基于Ripley L和O-ring函數(shù)的森林景觀空間分布格局及其關(guān)聯(lián)性[J]. 應(yīng)用生態(tài)學(xué)報，2014，25(12)：3429－3436.

Dong Lingbo, Liu Zhaogang, Zhang Bo, et al. Forest landscapes’spatial point patterns and associations based on Ripley L and O-ring functions[J]. Chinese Journal of Applied Ecology, 2014, 25(12): 3429－3436. (in Chinese with English abstract)

[17] Wiegand T, He F, Hubbell SP. A systematic comparison of summary characteristics for quantifying point patterns in ecology[J]. Ecography, 2013, 36(1): 92－103.

[18] Weiler M, Naef F. An experimental tracer study of the role of macropores in infiltration in grassland soils[J]. Hydrological Processes, 2003, 17(2): 477－493.

[19] Wiegand T, Moloney KA. Rings, circles, and null-models for point pattern analysis in ecology[J]. Oikos, 2004, 104(2): 209－229.

[20] 王鑫廳，梁存柱，王煒. 尺度與密度：測定不同尺度下的種群密度[J]. 植物生態(tài)學(xué)報，2013，37(2)：104－110.

Wang Xinting, Liang Cunzhu, Wang Wei. Scale and density: measuring local neighborhood density at different spatial scales[J]. Chinese Journal of Plant Ecology, 2013, 37(2): 104－110. (in Chinese with English abstract)

[21] 王康，張仁鐸，王富慶，等. 土壤水分運動空間變異性尺度效應(yīng)的染色示蹤入滲試驗研究[J]. 水科學(xué)進展，2007，18(2)：158－163.

Wang Kang, Zhang Renduo, Wang Fuqing, et al. Characterizing scale dependent heterogeneity of soil water movement through dying infiltration experiments[J]. Advances in Water Science, 2007, 18(2): 158－163. (in Chinese with English abstract)

[22] 張英虎，牛健植，杜曉晴，等. 鷲峰國家森林公園土壤優(yōu)先流現(xiàn)象分析[J]. 水土保持學(xué)報，2013，27(1)：41－45.

Zhang Yinghu, Niu Jianzhi, Du Xiaoqing, et al. Analysis of preferential flow phenomenon in Jiufeng Natiional Forest Park[J]. Journal of Soil and Water Conservation, 2013, 27(1): 41－45. (in Chinese with English abstract)

[23] Zhang Zhongbin, Peng Xinhua, Zhou Hu, et al. Characterizing preferential flow in cracked paddy soils using computed tomography and breakthrough curve[J]. Soil & Tillage Research, 2015, 146: 53－65.

[24] 陳曉冰，張洪江，程金花，等. 基于數(shù)量生態(tài)學(xué)的農(nóng)地優(yōu)先路徑水平分布特征研究[J]. 農(nóng)業(yè)機械學(xué)報，2015，46(7)：130－138.

Chen Xiaobing, Zhang Hongjiang, Cheng Jinhua, et al. Horizontal distribution characteristics of preferential flow paths in farmland based on quantitative ecology[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 130－138. (in Chinese with English abstract)

[25] Allaire S E, Gupta S C, Nieber J, et al. Role of macropore continuity and tortuosity on solute transport in soils: 1.Effects of initial and boundary conditions[J]. Journal of Contaminant Hydrology, 2002, 58(3/4): 299－321.

[26] 徐宗恒，徐則民，官琦，等. 不同植被發(fā)育斜坡土體優(yōu)先流特征[J]. 山地學(xué)報，2012，30(5)：521－527.

Xu Zongheng, Xu Zemin, Guan Qi, et al. The characteristic of preferential flow in different vegetated slope soils[J]. Journal of Mountain Science, 2012, 30(5): 521－527. (in Chinese with English abstract)

[27] 王偉，張洪江，程金花，等. 四面山闊葉林土壤大孔隙特征與優(yōu)先流的關(guān)系[J]. 應(yīng)用生態(tài)學(xué)報，2010，21(5)：1217－1223.

Wang Wei, Zhang Hongjiang, Cheng Jinhua, et al. Macropore characteristics and its relationships with the preferential flow in broadleaved forest soils of Simian Mountains[J]. Chinese Journal of App lied Ecology, 2010, 21(5): 1217－1223. (in Chinese with English abstract)

[28] Cey E E, Rudolph D L. Field study of macropore flow processes using tension infiltration of a dye tracer in partially saturated soils[J]. Hydrological Processes, 2009, 23(12): 1768－1779.

[29] 李顯平. 碎石土的顆粒特征及其優(yōu)先流研究進展[J]. 中國科技信息，2012，20(6)：58－59.

Li Xianping. The Characteristic of Debris’ Grain Composition and Study on Debris Preferential Flow[J]. China Science and Technology Information, 2012, 20(6): 58－59. (in Chinese with English abstract)

[30] 王大力，尹澄清. 植物根孔在土壤生態(tài)系統(tǒng)中功能[J]. 生態(tài)學(xué)報，2000，20(5)：869－874.

Wang Dali, Yi Chengqing. Functions of root channels in the soil system[J]. Acta Ecologica Sinica, 2000, 20(5): 869－874. (in Chinese with English abstract)

Distribution characteristic analysis of preferential flow paths in highway slope soils based on O-ring statistic

Zhang Dongxu1, Cheng Jinhua1※, Wang Wei2, Zhang Hongjiang1, Liu Tao2

(1.,,100083,;2.,100029,)

The soil layer structure of the slope would be disturbed during highway construction, leading to heterogeneity of soil properties. Water and solute are always moving in a non-equilibrium flow. Introducing the theory of preferential flow to analyze the migration mechanism of chloride pollutants in slope soil is of great significance to reveal the law of the migration of industrial salt in the protection of road environment. In order to observe the distribution of preferential flow paths in slope soil directly, dye infiltration experiment was used in soils of low subgrade slope of Heda highway. Image processing techniques and morphological computation procedures were used to determine the influence radius and coordinates of preferential flow paths. Spatial point pattern analysis method (O-ring statistic) of landscape ecology with 19 times simulation tests by using Monte Carlo stochastic was used to analyze the spatial distribution characteristics and the spatial correlation of the preferential flow paths with different influence radius. The results showed that the development depth of preferential flow paths in slope soil could reach to 40 cm by dye infiltration experiment. Plant roots were concentrated in topsoil (0-10 cm), and there was no trace of plant roots and biological activities in deeper soil. The number of preferential flow paths with same influence radius decreased as the soil depth increased, while the number of preferential flow paths decreased with the increasing of influence radius in the same soil layer. In the scale of this research, preferential flow paths showed a random distribution state at the range of 1.0≤＜5.0 mm and an aggregation distribution state when≥5.0 mm. With the gradual increase of influence radius, aggregated distribution pattern was more evident. At the range of 5≤＜10.0 mm, preferential flow paths showed an uniform distribution state in small scale (0-20 mm). The peak of aggregation state usually appeared at a smaller scale (0-40 mm). A large number of gravels lead to a change in the direction of actual movement of water flow, and preferential flow paths will be separated and then merged. The difference of soil particle composition can also lead to lateral water flow, which shows poor correlation between adjacent soil layers. The spatial correlation analysis of preferential flow path shows that the correlation between different preferential flow paths is relatively independent which related to the structure of soil. Large pores in slope soil are mainly formed by the filling of the fine particles of gravel soil, which has strong randomness, good stability, and not easy to collapse. This kind of non-random distribution condition shows that it is very necessary to study the distribution of preferential flow paths in slope soil. O-ring statistic is much more simple and intuitive than the traditional water breakthrough curve method, and it can enrich soil preferential flow research methods. It can also quantitatively analyze the distribution characteristics and spatial correlation of soil preferential flow paths from a statistical point of view. In future, landscape ecology spatial point pattern analysis method will be used well in study on spatial distribution of soil preferential paths.

statistics; soils; infiltration; spatial point pattern analysis; highway slope; preferential flow; paths; distribution characteristics

10.11975/j.issn.1002-6819.2017.04.023

S152.5

A

1002-6819(2017)-04-0161-08

2016-05-30

2016-07-18

國家自然科學(xué)基金資助項目（41271300、30900866）；中央級公益性科研院所基本科研業(yè)務(wù)費項目（20140601）。

張東旭，男，山東濟寧人，博士生，研究方向為水土保持。北京 北京林業(yè)大學(xué)水土保持學(xué)院，100083。Email：dongxu_zhang1989@126.com

程金花，女，山東聊城人，教授，主要從事土壤侵蝕與流域管理研究。北京 北京林業(yè)大學(xué)水土保持學(xué)院，100083。 Email：Jinhua_cheng@126.com

張東旭，程金花，王 偉，張洪江，劉 濤. 基于O-ring統(tǒng)計的公路邊坡土壤優(yōu)先流路徑分布分析[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(4)：161－168. doi：10.11975/j.issn.1002-6819.2017.04.023 http://www.tcsae.org

Zhang Dongxu, Cheng Jinhua, Wang Wei, Zhang Hongjiang, Liu Tao. Distribution characteristic analysis of preferential flow paths in highway slope soils based on O-ring statistic[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(4): 161－168. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.04.023 http://www.tcsae.org