詹 松，王文龍，黃鵬飛，李宏偉，李建明，王 貞，羅 婷

（1.西北農(nóng)林科技大學(xué) 資源與環(huán)境學(xué)院，陜西 楊凌712100；2.中國科學(xué)院 水利部 水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點(diǎn)試驗(yàn)室，陜西 楊凌712100；3.西北農(nóng)林科技大學(xué) 水土保持研究所，陜西 楊凌712100；4.江西省水土保持科學(xué)研究院 土壤侵蝕與防治重點(diǎn)試驗(yàn)室，江西 南昌330029；5.浙江省水利水電勘測設(shè)計(jì)院，浙江 杭州；6.中國水電顧問集團(tuán) 華東勘測設(shè)計(jì)研究院，浙江 杭州310014；7.杭州大地科技有限公司，浙江 杭州310000）

道路是土地利用方式之一，道路的形成實(shí)際是對原始地貌環(huán)境進(jìn)行改造的過程。道路使原地面地形地貌、地表植被等發(fā)生改變，相關(guān)物理性質(zhì)也發(fā)生較大變化，自然降雨作用下，其土壤侵蝕過程和方式與農(nóng)耕地有所不同［1］。國外學(xué)者［2］較早已開始了道路侵蝕研究，認(rèn)為道路侵蝕是河流泥沙來源之一。道路形成后，經(jīng)過不斷碾壓，入滲率降低，徑流系數(shù)增大，排水系統(tǒng)密度往往較其它下墊面大，道路因易于匯水而常常在其周圍沖刷成溝頭，據(jù)調(diào)查，高達(dá)58%公路通常直接或間接與河流相接，而公路周圍形成的溝道是主要連接方式之一［3］。隨著道路侵蝕研究的不斷深入，針對不同的公路、鐵路、林區(qū)和礦區(qū)等非硬化和硬化道路，學(xué)者們［4］從不同類型道路侵蝕過程、防護(hù)措施、防護(hù)效益展開了一系列研究，普遍認(rèn)為道路侵蝕可以使流域產(chǎn)沙量增加，而采取適當(dāng)?shù)姆雷o(hù)措施可以有效減少道路產(chǎn)生的泥沙量。

由于道路類型復(fù)雜多樣，因而道路侵蝕預(yù)測研究一直以來進(jìn)展緩慢，還無法對道路侵蝕進(jìn)行有效地估算，其原因主要是對相關(guān)機(jī)理的認(rèn)識還不夠深入，國內(nèi)研究主要對產(chǎn)流產(chǎn)沙和防護(hù)措施效益進(jìn)行相關(guān)的探討。王貞［5］采用野外放水沖刷試驗(yàn)分析了非硬化路面產(chǎn)流產(chǎn)沙特征。崔歡虎［6］探討了不同下墊面產(chǎn)流特征，認(rèn)為道路徑流損失可達(dá)80%以上。田鳳霞、劉剛、鄭世清等［7－9］分別通過室內(nèi)人工降雨或放水沖刷試驗(yàn)探討了不同雨強(qiáng)，不同坡度和不同植被覆蓋度下植物路產(chǎn)流產(chǎn)沙過程和水動力學(xué)特征，認(rèn)為水力學(xué)特征受植物路蓋度影響較大。曹世雄等［10］通過在農(nóng)田田間路面多年種植牧草并進(jìn)行定點(diǎn)監(jiān)測后認(rèn)為其侵蝕量可減少54%～78%，并有助于當(dāng)?shù)鼐坝^生態(tài)的恢復(fù)。

綜上所述，有關(guān)道路的研究主要集中在探討路面的產(chǎn)流產(chǎn)沙特征，而對水力學(xué)特征及其相關(guān)響應(yīng)關(guān)系研究較少，特別是在開發(fā)建設(shè)項(xiàng)目中，煤礦區(qū)非硬化路面的研究顯得更為薄弱。為加強(qiáng)礦區(qū)非硬化路面機(jī)理的深入探討，通過野外人工模擬降雨試驗(yàn)對比分析非硬化路面與原生地面的產(chǎn)流產(chǎn)沙、水動力參數(shù)變化特征及其響應(yīng)關(guān)系，為礦區(qū)非硬化道路水土流失的預(yù)測與治理提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

本文研究地點(diǎn)為神府東勝礦區(qū)。該區(qū)位于陜西府谷縣西部、神木縣北部和內(nèi)蒙古自治區(qū)鄂爾多斯（伊克昭盟）的南部，內(nèi)有烏蘭木倫河和窟野河流經(jīng)（108°36′—110°45′E，37°20′—39°30′N），屬于水蝕風(fēng)蝕交錯(cuò)區(qū)。該區(qū)屬于干旱、半干旱大陸季風(fēng)氣候區(qū)，年降雨量約在194.7～530.5mm，年際變化較大，且年內(nèi)分布不均，多暴雨，多集中在夏、秋季，蒸發(fā)量約為降雨量的5倍。土壤以黃土和風(fēng)沙土為主，結(jié)構(gòu)松散，抗蝕性差，植被覆蓋率低。受大風(fēng)和暴雨影響，水土流失嚴(yán)重，年平均侵蝕模數(shù)達(dá)1.5×104t／（km2·a），易發(fā)生滑坡、崩塌等自然災(zāi)害。

1.2 下墊面條件

試驗(yàn)于2010年8—9月在神東煤田區(qū)神木縣西溝鄉(xiāng)六道溝村一塊撂荒地進(jìn)行，該地地表植被覆蓋度相近，原生地面小區(qū)保持同一坡度內(nèi)地面平整并不做任何處理。非硬化路面是在原生地面上經(jīng)過不斷地碾壓而形成的，地表植被不復(fù)存在，地面硬度和強(qiáng)度都大于原生地面，采用控制土壤密度的方法人工構(gòu)筑非硬化路面，先測算多條實(shí)際非硬化路面土壤密度范圍（1.60～1.70g／cm3），然后用鐵鍬翻動原生地面20cm深并除去植被，整平后灑水夯實(shí)，并分3次分別覆土3～5cm厚度，前兩次灑水夯實(shí)，最后一次夯實(shí)不灑水并用環(huán)刀法測土壤密度（控制在1.70g／cm3左右），TDR測土壤含水量，達(dá)到要求后即可準(zhǔn)備試驗(yàn)。通過過篩分析，試驗(yàn)地土壤機(jī)械組成詳見表1［5］。

表1 試驗(yàn)地土體物理特性

1.3 試驗(yàn)小區(qū)布設(shè)

試驗(yàn)小區(qū)長寬為3m×1m，四周用40cm寬、1mm厚的鋼板插入地下25cm，地表出露15cm，為防止邊際效應(yīng)，小區(qū)周圍10cm內(nèi)處理與小區(qū)內(nèi)相同，小區(qū)下端設(shè)置集流槽。降雨設(shè)施采用12根6m，20根3m的6mm鋼管搭設(shè)3m高的降雨架，并用防風(fēng)布包圍四周以防止風(fēng)對降雨的影響。為保證降雨雨滴達(dá)到自然降雨終點(diǎn)速度，試驗(yàn)采用有壓供水，即利用潛水泵供水從小區(qū)旁8m高處的蓄水箱內(nèi)直接給降雨噴頭供水，通過壓力表和閥門控制降雨強(qiáng)度。由于試驗(yàn)地形條件限制，原生地面選擇坡度為5°，10°和18°，非硬化路面坡度為3°，6°，9°和12°，初始降雨強(qiáng)度都設(shè)置為1.0，1.5，2.0，2.5和3.0mm／min（數(shù)據(jù)處理時(shí)以實(shí)際量測為準(zhǔn)）。

試驗(yàn)開始前用塑料布蓋住試驗(yàn)小區(qū)，用量筒呈梅花形測定1min內(nèi)降雨量，即降雨強(qiáng)度，并采用克里斯琴森系數(shù)CU計(jì)算降雨均勻度。雨強(qiáng)達(dá)到要求后，揭開小區(qū)遮蓋的塑料布并開始記錄產(chǎn)流時(shí)間，開始產(chǎn)流時(shí)記下產(chǎn)流時(shí)間并接取徑流泥沙樣，測定水深、徑流寬和流速等數(shù)據(jù)。產(chǎn)流后前3min內(nèi)每分鐘用取樣桶接取一次泥沙樣，高猛酸鉀示蹤法在同一斷面不同水流位置測定3次流速（觀測段為2m），測針法測同一斷面水深3次和尺子測徑流寬3次。產(chǎn)流3min后，每3min重復(fù)一次上述過程。單次降雨時(shí)間總計(jì)45min（產(chǎn)流時(shí)間除外），并采用比重瓶法獲取徑流含沙量。

1.4 數(shù)據(jù)計(jì)算

降雨形成的徑流是引起土壤侵蝕的主要動力，對徑流的相關(guān)水力參數(shù)進(jìn)行計(jì)算有助于深入揭示非硬化路面產(chǎn)沙的本質(zhì)。運(yùn)用水力學(xué)相關(guān)理論和方法對徑流流速（V）、雷諾數(shù)（Re）、弗勞德數(shù)（Fr）、曼寧系數(shù)（n）、Darcy—Weisbach阻力系數(shù)、水流剪切力（τ）、水流功率（ω）等水力指標(biāo)進(jìn)行了計(jì)算。

雷諾數(shù)（Re）是水流慣性力與黏粘滯力的比值，是判定水流流態(tài)的重要依據(jù)。Re＜500時(shí)，水流為層流；Re＞2 000時(shí)，流水為紊流，計(jì)算公式為：

式中：V——徑流流速（m／s）；R——水力半徑（m）；Vm——含沙水流運(yùn)動黏滯性系數(shù)

式中：Sv——體積含沙量百分?jǐn)?shù)；d50——懸沙（試驗(yàn)上用）的中值粒徑（mm）；V0——運(yùn)動黏滯系數(shù)（m2／s），其值根據(jù)所測水溫查水力學(xué)表得到。

弗勞德數(shù)Fr是水流慣性力與重力的比值。如果Fr＜1，水流為緩流；Fr≥1為急流。計(jì)算公式為：

式中：g——重力加速度，取值9.8m／s2。

曼寧系數(shù)n和達(dá)西阻力系數(shù)f反映了下墊面對水流阻力作用，計(jì)算式為：

式中：J——水流能坡，近似為坡度正切值；h——水流深度（m）。

水流剪切力τ反映徑流對土壤的剝蝕能力，計(jì)算公式為：

式中：g——含沙水的密度（g／cm3）。

水流功率（ω）反映了一定高度的水體所具有的勢能，為單位面積水體勢能隨時(shí)間的變化率，計(jì)算公式為：

2 結(jié)果與分析

2.1 原生地面與非硬化路面水力參數(shù)比較

2.1.1 流速比較 坡面徑流流速是反映下墊面物理性質(zhì)的重要水力參數(shù)，也是計(jì)算其它水力參數(shù)的基礎(chǔ)。在不同坡度、雨強(qiáng)和下墊面等組合處理中，以測得的流速作為各時(shí)段平均流速，并對相對穩(wěn)定后的流速取平均值作為全時(shí)段的平均流速，根據(jù)雨強(qiáng)、坡度繪制流速變化關(guān)系圖（圖1）。如圖1所示，非硬化路面在不同坡度不同流量下的流速均明顯大于原生地面流速，主要是原生地面有植被覆蓋，增大了地表阻力系數(shù)，徑流隨坡面流動過程中受植被阻擋，流速相對變緩，而非硬化路面缺乏植被阻礙，流速較大，流速范圍為0.13～0.25m／s，最大流速值較劉剛、田鳳霞等［7－8］裸露土質(zhì)路面流速偏大，但都在同一個(gè)數(shù)量級，主要受土壤密度差異影響。兩下墊面相近坡度（3°與5°，5°與6°，10°與9°，12°與10°）進(jìn)行流速比較發(fā)現(xiàn)，相同降雨強(qiáng)度下，非硬化路面流速是原生地面的2.2～4.2倍，且較穩(wěn)定，可認(rèn)為相同條件下，兩者的流速比值為一常數(shù)，該值不隨坡度、雨強(qiáng)而發(fā)生較大變化，表明現(xiàn)有采用流量和坡度而不考慮下墊面影響推求流速的方法是可行的。非硬化路面和原生地面流速隨雨強(qiáng)增大而增大，但變化趨勢有所不同。原生地面變化平緩，可用二次函數(shù)進(jìn)行模擬：

式中：V——流速；I——雨強(qiáng)；a，b，c——常數(shù)。

對2mm／min雨強(qiáng)、10°坡度條件下流速修正后模擬表明，各坡度下，決定系數(shù)達(dá)到0.91以上。非硬化路面流速隨雨強(qiáng)變化波動較大，隨坡度增大，波動幅度越大，主要是在坡面下部可見細(xì)溝、切溝等，由此產(chǎn)生跌水耗能而影響徑流流速。而坡度對流速的變化不顯著，這與Nearing［11］結(jié)論相似。

2.1.2 水力參數(shù)隨坡度、雨強(qiáng)變化特點(diǎn)對比 根據(jù)穩(wěn)定后平均流速、水深等相關(guān)參數(shù)和公式計(jì)算獲取水流剪切力（t）、功率（ω）、弗勞德數(shù)（Fr）、雷諾數(shù)（Re）等水力參數(shù)，結(jié)果如圖2所示。

水流剪切力和水流功率是反映徑流對土壤剝蝕能力的一組水力參數(shù)，由圖2可知，同一下墊面水流剪切力和水流功率隨坡度和雨強(qiáng)的增大而增大。在坡度較大時(shí)，非硬化路面水流剪切力隨雨強(qiáng)先增大后減小，但減小幅度較小，原因是在小雨強(qiáng)下，坡面流以薄層水流運(yùn)動，隨著雨強(qiáng)的增大，水流形態(tài)發(fā)生變化，出現(xiàn)滾波流［12］，水流水深不斷變化而不穩(wěn)定，而水面寬度則有所增加，但坡面流不同于明渠水流，水深和水寬數(shù)值都較小，在一定范圍內(nèi)，水力半徑會有所減小，從而使剪切力偏小。在相近坡度下，非硬化路面水流剪切力和水流功率分別比原生地面增大了0.4～9倍，0.8～20倍，雨強(qiáng)越大，增大倍數(shù)越小。

圖1 流速隨坡度和雨強(qiáng)的變化

試驗(yàn)中，F(xiàn)r均小于1，表明水流為緩流，F(xiàn)r越小，水流流動越緩慢。在同一雨強(qiáng)相近坡度試驗(yàn)中，非硬化路面Fr均大于原生地面（包括18°坡面），平均增大了66%。非硬化路面Re均大于500，甚至大于2 000，即多數(shù)情況下流態(tài)屬于過渡流，少數(shù)為紊流，并隨雨強(qiáng)增大有增大趨勢。原生地面Re絕大部分情況下在34～400之間，小雨強(qiáng)時(shí)低于100，僅在大雨強(qiáng)大坡度時(shí)超過1 000，因此原生地面流態(tài)基本為緩流。

非硬化路面和原始路面達(dá)西阻力系數(shù)、曼寧系數(shù)均隨坡度增大而增大，但非硬化路面達(dá)西阻力系數(shù)變化幅度比原生地面小。非硬化路面達(dá)西阻力系數(shù)隨雨強(qiáng)變化不明顯，而原生地面隨著雨強(qiáng)的增大先增大后減小。非硬化路面曼寧系數(shù)是原生地面的10.9%～53.5%，雨強(qiáng)較小時(shí)，非硬化路面曼寧系數(shù)大于原生地面，隨著雨強(qiáng)的增大而發(fā)生波動，雨強(qiáng)較大時(shí)非硬化路面曼寧系數(shù)低于原生地面，即小雨強(qiáng)時(shí)與田鳳霞植物路曼寧系數(shù)比土質(zhì)路大的結(jié)論有出入［7］，可能是小雨強(qiáng)時(shí)原生地面受植被影響徑流水深比非硬化路面小。

圖2 原生地面（5°，10°，18°）與非硬化路面（3°，6°，9°，12°）水力參數(shù)變化

2.2 非硬化路面產(chǎn)沙分析

2.2.1 產(chǎn)沙量和徑流量關(guān)系對比分析 非硬化路面屬于土質(zhì)路面，是原生地面經(jīng)過人畜和車輛碾壓后形成的。其物理性質(zhì)與原生地面有較大區(qū)別，路面的產(chǎn)沙量也隨著物理特性的改變而相應(yīng)發(fā)生改變。產(chǎn)沙量與徑流量關(guān)系如圖3所示，非硬化路面產(chǎn)沙量隨徑流增大明顯，坡度和徑流量越大，產(chǎn)沙量也越大，產(chǎn)沙量上升趨勢明顯，產(chǎn)沙量和徑流量關(guān)系式為：

式中：M——產(chǎn)沙量；Q——徑流量；a，n——常數(shù)。

原生地面產(chǎn)沙量（M原生）明顯比非硬化路面產(chǎn)沙量M非硬少，最大值不超過5kg，相近坡度不同流量下僅為非硬化路面的0.8%～14%，產(chǎn)沙量隨徑流的變化關(guān)系同樣可用冪函數(shù)進(jìn)行模擬，決定系數(shù)在0.94以上，但參數(shù)a比非硬化路面低3個(gè)數(shù)量級。

圖3 原生地面與非硬化路面產(chǎn)沙量與徑流量關(guān)系

2.2.2 侵蝕模數(shù)與平均剪切力關(guān)系對比分析 侵蝕模數(shù)是侵蝕過程單位面積上土壤流失量。分析可知，非硬化路面侵蝕模數(shù)與剪切力之間關(guān)系不明顯，各點(diǎn)離散度較大，坐標(biāo)軸上呈“扇形”較均勻分布，而原生地面分布比較集中，呈“箭”形狀分布，多數(shù)點(diǎn)集中箭頭部分，即剪切力0～10N／m2，產(chǎn)沙量0～1.5kg范圍內(nèi)，對兩下墊面分別用指數(shù)函數(shù)和二次函數(shù)進(jìn)行擬合：

非硬化路面擬合式?jīng)Q定系數(shù)R2僅為0.362，而原生地面二次函數(shù)的決定系數(shù)達(dá)0.849。由公式（9）—（10）可知，兩者發(fā)生侵蝕的臨界剪切力不同，原生地面的臨界剪切力為1.1N／m2，而非硬化路面結(jié)果與實(shí)際不符，非硬化路面擬合公式不符合該地面特征。（1）在坡度較大時(shí)，非硬化路面水流隨雨強(qiáng)增大，水流以滾動方式流動，水流對地表作用力不均勻，地表易形成跌坑，跌坑形成后對地面的作用力不僅有切應(yīng)力，還有跌水重力作用，增大了地表產(chǎn)沙量，從而導(dǎo)致擬合結(jié)果出現(xiàn)偏差；（2）剪切力采用的是平均剪切力，而非瞬時(shí)剪切力。因此原生地面侵蝕模數(shù)與剪切力關(guān)系密切，非硬化路面受其它因素影響關(guān)系較復(fù)雜。

2.2.3 侵蝕模數(shù)與平均水流功率關(guān)系對比分析 分析可知，非硬化路面和原生地面水流侵蝕模數(shù)隨水流功率增大而增大，非硬化路面可用二次函數(shù)描述，原生地面可用對數(shù)函數(shù)描述：

從公式（11）—（12）可知，非硬化路面發(fā)生侵蝕時(shí)不存在臨界功率，而原生地面臨界值為0.031N／（m·s），結(jié)果似乎與部分放水沖刷試驗(yàn)結(jié)論相悖［13］，但從試驗(yàn)過程分析，結(jié)論并不矛盾，在非硬化路面上，降雨與沖刷不同，在徑流形成之前就會對地表形成擊濺侵蝕，在徑流形成后，在剪切力未達(dá)到臨界剪切力之前徑流不會引起新的侵蝕，那么此時(shí)就只對擊濺侵蝕土壤進(jìn)行搬運(yùn)，就造成了未達(dá)到臨界侵蝕力就造成侵蝕的假象，而原生地面有植被保護(hù)，雨滴直接作用不大，主要以形成的徑流侵蝕為主。

2.2.4 侵蝕模數(shù)與動能關(guān)系對比分析 坡面侵蝕量大小還與能量有關(guān)，坡面侵蝕的過程其實(shí)是能量變化過程，而含沙水流的動能大小直接反映水流對侵蝕面的剝蝕能力大小。通過將侵蝕模數(shù)與坡面過程總動能建立關(guān)系，有助于對侵蝕量進(jìn)行模擬和預(yù)測。動能采用動能定理公式獲取，即：

式中：Q——坡面全過程徑流總量（m3）；V——平均流速（m／s）；E——總動能（J）。

分析可知，非硬化路面和原生地面侵蝕模數(shù)隨動能增大而增大，非硬化路面侵蝕總動能明顯大于原生地面，而兩者降雨總動能相同，造成該結(jié)果主要原因有：（1）入滲因素。非硬化路面經(jīng)過碾壓，土壤密度大于原生地面，土壤孔隙度低于原生地面，不易形成壤中通路，雨滴降落到地面后不易入滲，主要以坡面徑流形式存在；（2）植被因素。原生地面有植被存在，雨滴在與地面接觸之前首先接觸植被，植被莖葉對雨滴進(jìn)行攔截，降低了雨滴動能，使最終降落到地面的雨滴能量低于試驗(yàn)設(shè)計(jì)總動能。通過二次函數(shù)和線性擬合，非硬化路面和原始路面侵蝕模數(shù)與動能關(guān)系分別為表達(dá)式（14）—（15）。

通過表達(dá)式（14）—（15）分析可知，非硬化路面和原生地面都不存在臨界動能，令E＝0，發(fā)現(xiàn)M原生＝0.272，M非硬＝6.188，即在徑流有初始動能之前就發(fā)生了侵蝕，可以推測是有雨滴直接作用于地面而引起侵蝕，即雨滴擊濺侵蝕，此過程雨滴在降落到地面時(shí)具有一定的速度沖量，在與地表作用的瞬間使受打擊的土壤發(fā)生位移。從大小看，M非硬＞M原生，說明雨滴對非硬化地面作用遠(yuǎn)大于原生地面，結(jié)合上述總動能影響因素推斷，結(jié)果符合降雨侵蝕實(shí)際過程，雨滴對非硬化路面侵蝕作用約是原生地面的30倍。

有關(guān)坡面侵蝕量定量描述方面，多數(shù)學(xué)者主要從力和功能兩個(gè)方向?qū)搅骷羟辛?、功率和動能進(jìn)行了探討，具體哪種方式對坡面侵蝕描述結(jié)論不一，本研究探討了侵蝕模數(shù)對徑流平均水流剪切力、平均水流功率和總動能響應(yīng)關(guān)系，認(rèn)為徑流從能量角度描述坡面侵蝕較好，即動能可以較好描述侵蝕模數(shù)與水動力參數(shù)的關(guān)系（徑流動能＞平均水流功率＞平均水流剪切力），結(jié)果與劉俊娥［14］坡面片流降雨試驗(yàn)相似。

3 結(jié)論

（1）通過降雨試驗(yàn)比較了非硬化路面和原生地面的水動力參數(shù)變化，結(jié)果表明非硬化路面與原生地面差異較大。由雷諾數(shù)判斷，在試驗(yàn)條件下，非硬化路面流態(tài)主要為過渡流，而原生地面主要為層流，但都表現(xiàn)為緩流。非硬化路面水流流速、剪切力、佛汝德數(shù)、雷諾數(shù)和功率均大于原生地面，曼寧系數(shù)、達(dá)西阻力系數(shù)小于原生地面。

（2）非硬化路面產(chǎn)沙量隨徑流量增大，而原生地面受地面植被對降雨的影響，產(chǎn)沙量隨徑流量變化緩慢，產(chǎn)沙量僅為非硬化路面的0.8%～14%，因此，減少對原生地面的干擾可以減少水土流失。

（3）非硬化路面侵蝕模數(shù)與水動力參數(shù)響應(yīng)關(guān)系與原生地面不一致。非硬化路面坡面侵蝕模數(shù)與平均剪切力、平均水流功率相關(guān)性不大，而原生地面與平均剪切力相關(guān)性明顯，與平均水流功率不明顯。但總體而言，非硬化路面侵蝕模數(shù)隨動能響應(yīng)明顯，決定系數(shù)較剪切力和水流功率高，能很好地描述降雨侵蝕過程，即非硬化路面侵蝕模數(shù)對水動力參數(shù)響應(yīng)的強(qiáng)弱關(guān)系為：徑流動能＞平均水流功率＞平均水流剪切力。

［1］ 張科利，徐憲利，羅麗芳.國內(nèi)外道路侵蝕研究回顧與展望［J］.地理科學(xué)，2008，28（1）：119-123.

［2］ Fransen P J B，Phillips C J，F(xiàn)ahey B D.Forest road erosion in New Zealand：Overview［J］.Earth Surface Processes and Landforms，2001，26（2）：165-174.

［3］ Nyssen J，Poesen J，Moeyersons J，et al.Impact of road building on gully erosion risk：A case study from the northern Ethiopian Highland ［J］.Earth Surface Processes and Landforms，2002，27（12）：1267-1283.

［4］ Carr W W，Ballard T M.Hydroseeding forest roadsides in British Columbia for erosion control［J］.Journal of Soil and Water Conservation，1980，35（1）：33-35.

［5］ 王貞，王文龍，羅婷，等.非硬化路面侵蝕產(chǎn)沙規(guī)律野外模擬試驗(yàn)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2010，26（12）：79-83.

［6］ 崔歡虎，閆翠萍，靖華，等.黃土丘陵區(qū)不同下墊面降雨產(chǎn)流與蓄水窖工藝技術(shù)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2004，20（2）：101-104.

［7］ 田風(fēng)霞，劉剛，鄭世清，等.草本植物對土質(zhì)路面徑流水動力學(xué)特征及水沙過程的影響［J］.農(nóng)業(yè)工程學(xué)報(bào)，2009，25（10）：25-29.

［8］ 劉剛，鄭世清，田風(fēng)霞，等.坡度與降雨強(qiáng)度對植物路路面侵蝕動力參數(shù)的影響［J］.水土保持學(xué)報(bào)，2009，3（4）：44-49.

［9］ 鄭世清，田風(fēng)霞，王占禮，等.植物路與土質(zhì)路產(chǎn)流產(chǎn)沙過程的比較試驗(yàn)［J］.泥沙研究，2009（4）：1-6.

［10］ 曹世雄，陳莉，高旺盛.在黃土丘陵區(qū)土質(zhì)路面種草［J］.生態(tài)學(xué)報(bào)，2005，25（7）：1754-1763.

［11］ Nearing M A，Simanton J R，Norton L D，et al.Soil erosion by surface water flow on a stony，semiarid hillslope［J］.Earth Surface Processes and Landforms，1999，24（8）：677-686.

［12］ 張寬地，王光謙，呂宏興，等.模擬降雨條件下坡面流水動力學(xué)特性研究［J］.水科學(xué)進(jìn)展，2012，23（2）：229-235.

［13］ 羅婷，王文龍，王貞，等.非硬化土路土壤剝蝕率與水動力學(xué)參數(shù)分析［J］.人民黃河，2011，3（4）：96-102.

［14］ 劉俊娥，王占禮，高素娟，等.黃土坡面片蝕過程動力學(xué)機(jī)理試驗(yàn)研究［J］.農(nóng)業(yè)工程學(xué)報(bào)，2012，28（7）：144-149.