上方來(lái)水來(lái)沙對(duì)淺溝侵蝕產(chǎn)沙及水動(dòng)力參數(shù)的影響

車小力，王文龍，2?，郭軍權(quán)，韓芳芳

(1.西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院，712100，陜西楊凌;2.中國(guó)科學(xué)院水利部水土保持研究所土壤侵蝕與旱地農(nóng)業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，712100，陜西楊凌;3.延安職業(yè)技術(shù)學(xué)院，745000，陜西延安;4.長(zhǎng)安大學(xué)環(huán)境科學(xué)與工程學(xué)院，710064，西安)

土壤侵蝕在黃土高原地區(qū)非常嚴(yán)重，制約著當(dāng)?shù)亟?jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展，而淺溝侵蝕作為坡面土壤侵蝕的一個(gè)重要過(guò)程，在黃土地區(qū)坡耕地上的分布十分廣泛，是一種人為耕墾與強(qiáng)烈徑流侵蝕相互作用的特殊侵蝕形態(tài)，對(duì)黃土區(qū)溝蝕的發(fā)生發(fā)展有著重要的影響，歷來(lái)為學(xué)者所重視。淺溝是在徑流沖刷和人為耕作的反復(fù)作用下形成的瓦背狀起伏的地形［1］。有關(guān)土壤侵蝕過(guò)程的研究，在片蝕和細(xì)溝侵蝕研究方面已取得了重要的成果，但對(duì)淺溝侵蝕過(guò)程的研究仍較薄弱［2］。有關(guān)資料表明，淺溝占坡面范圍大，侵蝕量可占坡面侵蝕量的17% ～85%［3-6］。國(guó)內(nèi)外學(xué)者對(duì)淺溝侵蝕進(jìn)行了一定的研究，主要涉及淺溝侵蝕分布特征［7］、淺溝侵蝕過(guò)程及其影響因素［8-9］、淺溝水流水動(dòng)力參數(shù)特征［10］等，并建立了淺溝侵蝕預(yù)報(bào)模型［11-12］。鄭粉莉等［13］利用室內(nèi)人工建筑淺溝發(fā)育初期的雛形模型，研究淺溝發(fā)育不同階段的侵蝕過(guò)程及淺溝侵蝕對(duì)坡面侵蝕產(chǎn)沙的貢獻(xiàn)，指出淺溝發(fā)育不同階段對(duì)應(yīng)于不同的淺溝侵蝕過(guò)程，淺溝侵蝕量在發(fā)育初期和中期階段占總坡面侵蝕產(chǎn)沙量的58%，而在后期為26% ～59%。王文龍等［14］通過(guò)模擬降雨試驗(yàn)的方法對(duì)比了淺溝與細(xì)溝水動(dòng)力參數(shù)的關(guān)系，指出紊動(dòng)和急流沖刷是淺溝區(qū)別于細(xì)溝侵蝕的最重要的特征。鄭粉莉等［15］利用室內(nèi)雙土槽徑流小區(qū)，通過(guò)模擬降雨試驗(yàn)研究了坡度為20°時(shí)上方來(lái)水來(lái)沙對(duì)坡下方淺溝侵蝕產(chǎn)沙的影響，結(jié)果表明，上方來(lái)水引起的侵蝕產(chǎn)沙量隨降雨強(qiáng)度的增加可增加4.4% ～83.5%，而在相同上方來(lái)水條件下，上方來(lái)沙強(qiáng)度的減少，可使下方淺溝侵蝕產(chǎn)沙量增加。由于淺溝侵蝕過(guò)程的復(fù)雜性，現(xiàn)有研究成果在淺溝水流水動(dòng)力參數(shù)方面的試驗(yàn)研究在黃土高原地區(qū)較少，且以室內(nèi)人工模擬試驗(yàn)居多［8-14］，而野外徑流小區(qū)觀測(cè)資料表明，上方來(lái)水來(lái)沙對(duì)坡面侵蝕具有重要的影響［15-16］。因此，筆者通過(guò)野外放水沖刷試驗(yàn)，定量研究上方來(lái)水來(lái)沙對(duì)坡面淺溝侵蝕產(chǎn)沙及其水動(dòng)力參數(shù)的影響，以期為淺溝侵蝕過(guò)程定量化研究和坡面土壤侵蝕預(yù)報(bào)模型的建立提供科學(xué)依據(jù)。

1 研究區(qū)概況

研究區(qū)位于陜西省延安市燕兒溝小流域內(nèi)，地處延安市南3 km處，地勢(shì)西北高，東南低，地形以梁、峁為主。土壤為典型的黃綿土，粉粒含量多，有機(jī)質(zhì)含量為0.9%，耕作層土壤密度為1.15 g/m3左右。該區(qū)年平均降雨量480～630 mm，降雨年際變化大，多集中于夏季，雨型以短歷時(shí)暴雨為主。

2 試驗(yàn)設(shè)計(jì)與研究方法

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

試驗(yàn)地布設(shè)在燕兒溝流域廟河村對(duì)面山上的蘋果園里，盡可能選取原來(lái)就發(fā)育有淺溝的完整坡面，或依照具體地形修整成兩肩高中間低的瓦背狀淺溝。設(shè)寬1 m、長(zhǎng)10 m、坡度為26°淺溝小區(qū)。試驗(yàn)利用水窖供水，采用水泵和閥門控制流量，用標(biāo)準(zhǔn)徑流桶對(duì)放水流量進(jìn)行率定。將10 m小區(qū)分為上下2個(gè)，通過(guò)小區(qū)相連的方式來(lái)模擬上方來(lái)水對(duì)下方小區(qū)水動(dòng)力學(xué)參數(shù)及產(chǎn)沙的影響。在小區(qū)上方放置一個(gè)與小區(qū)同寬的溢流槽，保證試驗(yàn)中水是均勻地從溢流槽流出，在小區(qū)中部，將有無(wú)數(shù)小孔的塑料水管，用紗布包裹，并使其形狀與坡面淺溝凹狀相符，將水管與水箱相連，保證試驗(yàn)中水流均勻地流向下方坡面。順小區(qū)縱向2、4、6和8 m處設(shè)定4個(gè)觀測(cè)斷面，以測(cè)定流速、流寬和流深。小區(qū)下方設(shè)有集流槽，根據(jù)集流桶內(nèi)的水沙體積可計(jì)算出徑流量與泥沙量。放水流量分別采用上方放水流量5 L/min、下方放水流量5 L/min;上方放水流量10 L/min、下方放水流量10 L/min;上方放水流量15 L/min、下方放水流量10 L/min;另外還對(duì)下方5 m小區(qū)進(jìn)行放水流量為5和10 L/min的無(wú)上方來(lái)水來(lái)沙試驗(yàn)。單場(chǎng)試驗(yàn)重復(fù)次數(shù)為2。野外放水沖刷試驗(yàn)系統(tǒng)示意圖見(jiàn)圖1。

圖1 野外放水沖刷試驗(yàn)系統(tǒng)示意圖Fig.1 Sketch map of outdoor scouring erosion experimental apparatus

2.2 試驗(yàn)過(guò)程及數(shù)據(jù)獲取

試驗(yàn)前測(cè)量土壤密度和前期含水量，試驗(yàn)開(kāi)始，記錄產(chǎn)流時(shí)間，產(chǎn)流后前3 min內(nèi)每1 min接1次水樣，3 min后每3 min接1次水樣。用體積法量測(cè)流量，用取樣烘干法測(cè)含沙量，用高錳酸鉀溶液示蹤法每隔5 min測(cè)量各斷面處流速，用測(cè)針?lè)繙y(cè)斷面水深，流寬用卡尺量測(cè)，試驗(yàn)歷時(shí)40 min。

通過(guò)野外放水沖刷試驗(yàn)，用斷面平均流速、雷諾數(shù)、弗勞德數(shù)、阻力系數(shù)、剪切力、水流功率表示淺溝水流的水力學(xué)參數(shù)。主要計(jì)算公式如下。

斷面平均流速v為各時(shí)段流速大小的平均值，單位是m/s。

式中:Re為雷諾數(shù)，是水流慣性力與黏性力之比，量綱為1，Re＜500，水流為層流，Re＞500，水流處于紊流狀態(tài);v為平均流速，m/s;R為水力半徑，m;ν為運(yùn)動(dòng)黏性系數(shù)，m2/s。

式中:Fr為弗勞德數(shù)，是水流的慣性力與重力的比值?？梢杂肍r反映水流流態(tài)，當(dāng)Fr＜1時(shí)水流為緩流，當(dāng)Fr＞1時(shí)水流為急流;g為重力加速度，m/s2。

式中:f為阻力系數(shù)，表征了水流在流經(jīng)坡面時(shí)所受的阻力情況;J為水力能坡，m/m;v1為斷面平均流速，m/s。

式中:S為水流剪切力，Pa，是反應(yīng)水流對(duì)淺溝溝床土壤剝蝕力大小的參數(shù);Q為水流平均含沙量時(shí)的密度，kg/m3，取 1.8 kg/m3。

式中:P為水流功率，N/(m·s)，表征水流流動(dòng)時(shí)的挾沙能力的參數(shù)。

3 結(jié)果與分析

3.1 上方來(lái)水對(duì)淺溝水動(dòng)力參數(shù)的影響

3.1.1 上方來(lái)水對(duì)淺溝水流流速的影響 不同放水流量下，上方含沙水流匯入淺溝侵蝕槽后，導(dǎo)致淺溝水流流速明顯增大。從圖2可以看出，上方匯水匯沙導(dǎo)致淺溝下方流速明顯增大，且隨著上方匯水流量的增大而增大。由于受微地貌的影響，淺溝侵蝕過(guò)程流速的變化是不均勻的，隨著侵蝕的發(fā)生，水流流速出現(xiàn)整體遞減的趨勢(shì)，說(shuō)明坡面淺溝侵蝕過(guò)程以侵蝕-搬運(yùn)為主，隨著上方泥沙的運(yùn)移，下方水流進(jìn)行挾沙需要消耗較多的能量，其水流沖擊能力相對(duì)減少，因此，淺溝水流流速變小。

圖2 26°時(shí)不同放水流量有無(wú)上方來(lái)水流速隨時(shí)間的變化Fig.2 Flow velocity over time with or without up-slope runoff in different water discharge at slope of 26°

3.1.2 上方來(lái)水對(duì)淺溝水流雷諾數(shù)的影響 有無(wú)上方匯流對(duì)淺溝水流雷諾數(shù)具有明顯的影響。雷諾數(shù)隨著放水流量的增大而增大，上方來(lái)水來(lái)沙的匯入，使下坡淺溝水流的紊動(dòng)性顯著增加，雷諾數(shù)變大，從而增大了淺溝侵蝕產(chǎn)沙量。在坡度為26°時(shí)，無(wú)上方來(lái)水時(shí)淺溝水流雷諾數(shù)變化于1 165～8 689之間，大于500，淺溝水流屬于紊流狀態(tài);與無(wú)上方來(lái)水相比，有上方來(lái)水時(shí)的淺溝水流雷諾數(shù)增大33% ～76%。

從圖3可以看出，雷諾數(shù)隨著放水流量的增大而增大，上方匯水匯沙使下方水流雷諾數(shù)呈現(xiàn)遞增趨勢(shì)，當(dāng)下坡放水流量相同時(shí)，雷諾數(shù)隨著上坡放水流量的增大而增大;因此，在坡度為26°時(shí)，淺溝水流屬于紊流，在上方來(lái)水來(lái)沙的作用下，淺溝水流雷諾數(shù)迅速增大，淺溝水流搬運(yùn)泥沙的能力增強(qiáng)，侵蝕產(chǎn)沙量增大。

圖3 26°時(shí)不同放水流量有無(wú)上方來(lái)水雷諾數(shù)隨時(shí)間的變化Fig.3 Reynolds number over time with or without up-slope runoff in different water discharge at slope of 26°

3.1.3 上方來(lái)水對(duì)淺溝水流弗勞德數(shù)的影響 弗勞德數(shù)隨著水流紊動(dòng)強(qiáng)度的增大而增大，上方來(lái)水來(lái)沙的匯入對(duì)淺溝水流弗勞德數(shù)有顯著的影響。無(wú)上方來(lái)水淺溝水流的弗勞德數(shù)介于1.21～4.35之間，其流態(tài)屬于急流。與無(wú)上方來(lái)水相比，上方來(lái)水使淺溝水流的弗勞德數(shù)增加21% ～47%，當(dāng)下坡放水量相同時(shí)，上方來(lái)水為15 L/min的弗勞德數(shù)比10 L/min的大26%。

從圖4可以看出，淺溝水流弗勞德數(shù)隨著上方來(lái)水來(lái)沙的匯入顯著增大，上方來(lái)水來(lái)沙的增大亦使弗勞德數(shù)變大;弗勞德數(shù)最大值出現(xiàn)在放水初期，最小值出現(xiàn)在放水末期，隨著放水時(shí)間的持續(xù)，淺溝形態(tài)趨于穩(wěn)定，水流紊動(dòng)強(qiáng)度減小，故弗勞德數(shù)隨放水時(shí)間呈遞減趨勢(shì)。

圖4 26°時(shí)不同放水流量有無(wú)上方來(lái)水弗勞德數(shù)隨時(shí)間的變化Fig.4 Froude number over time with or without up-slope runoff in different water discharge at slope of 26°

3.1.4 上方來(lái)水對(duì)淺溝水流阻力系數(shù)的影響 水流進(jìn)入淺溝后，必然受到摩擦阻力的影響，其阻力主要來(lái)自水流中沙粒本身的摩擦力對(duì)水流的阻礙作用。在坡度為26°時(shí)，淺溝水流的阻力系數(shù)介于0.69～1.98之間，上方來(lái)水的匯入使淺溝水流阻力系數(shù)減小，與無(wú)上方來(lái)水相比其值減小11% ～13%。

從圖5可以看出:放水流量與阻力系數(shù)成反比，上方來(lái)水匯入淺溝侵蝕槽后，阻力系數(shù)明顯減小，且隨著上方來(lái)水的增大而減小;整個(gè)放水過(guò)程中，由于上方來(lái)沙匯入下方淺溝侵蝕槽，水流摩擦力增大，故阻力系數(shù)隨著放水時(shí)間的持續(xù)逐漸變大。

3.1.5 上方來(lái)水對(duì)淺溝水流剪切力的影響 剪切力反應(yīng)水流對(duì)土壤剝蝕力的大小，上方來(lái)水進(jìn)入淺溝侵蝕槽后，增加了下方匯水流量，淺溝水流剪切力相應(yīng)增大。在坡度為26°時(shí)，不同放水流量淺溝水流剪切力介于1.763～4.521之間，與無(wú)上方來(lái)水相比，上方來(lái)水使淺溝水流剪切力增加18%～42%。

圖5 26°時(shí)不同放水流量有無(wú)上方來(lái)水阻力系數(shù)隨時(shí)間的變化Fig.5 Daycy-weisbach coefficient over time with or without up-slope runoff in different water discharge at slope of 26°

從圖6可以看出:淺溝侵蝕槽在接受上方來(lái)水后，其水流剪切力顯著增大，同時(shí)，放水流量的增加也使淺溝水流剪切力增大。淺溝水流剪切力隨著放水時(shí)間的持續(xù)，呈現(xiàn)出遞增趨勢(shì)，只有在放水流量為“上5 L/min，下5 L/min”時(shí)，其值在整個(gè)放水階段比較平穩(wěn)。

圖6 26°時(shí)不同放水流量有無(wú)上方來(lái)水剪切力隨時(shí)間的變化Fig.6 Shearing force over time with or without up-slope runoff in different water discharge at slope of 26°

3.1.6 上方來(lái)水對(duì)淺溝水流功率的影響 水流功率表征水流流動(dòng)時(shí)的挾沙能力，流動(dòng)水體將勢(shì)能轉(zhuǎn)化為動(dòng)能，對(duì)淺溝侵蝕槽造成破壞。上方來(lái)水使淺溝水流的水流功率增加29% ～72%，同時(shí)，放水流量的增加，也使淺溝水流的水流功率增大。

從圖7可以看出:淺溝水流功率隨上方來(lái)水的匯入而顯著增大，下坡放水流量相同時(shí)，水流功率隨上方來(lái)水量的增大而增大;水流功率隨放水時(shí)間的持續(xù)而遞增，唯有在“上5 L/min，下5 L/min”時(shí)，其值波動(dòng)比較平緩。

圖7 26°時(shí)不同放水流量有無(wú)上方來(lái)水水流功率隨時(shí)間的變化Fig.7 Stream power over time with or without up-slope runoff in different water discharge at slope of 26°

3.2 上方來(lái)水對(duì)淺溝侵蝕產(chǎn)沙量的影響

不同放水流量對(duì)坡面淺溝侵蝕產(chǎn)沙過(guò)程具有顯著影響，而上方來(lái)水的匯入加劇了這一過(guò)程。根據(jù)實(shí)驗(yàn)設(shè)計(jì)分析上方來(lái)水對(duì)淺溝侵蝕產(chǎn)沙的影響，點(diǎn)繪出淺溝水流的產(chǎn)沙過(guò)程線。產(chǎn)沙過(guò)程大致為，隨著放水時(shí)間的延長(zhǎng)，產(chǎn)沙量在出現(xiàn)峰值后逐漸減少，最后趨于穩(wěn)定。

從圖8可以看出，上方來(lái)水對(duì)淺溝侵蝕的產(chǎn)沙量貢獻(xiàn)很大，上方來(lái)水量大的，產(chǎn)沙增量大。與無(wú)上方來(lái)水相比，有上方來(lái)水的淺溝侵蝕產(chǎn)沙過(guò)程較為復(fù)雜，有多峰多谷的特點(diǎn)。其最大值出現(xiàn)在產(chǎn)流初期，此時(shí)由于土粒松散，易被剝蝕，再加上上方來(lái)水匯入淺溝侵蝕槽，產(chǎn)沙過(guò)程線便出現(xiàn)峰值。隨著放水時(shí)間的延長(zhǎng)，待坡面土壤水分飽和時(shí)，即坡面產(chǎn)流達(dá)到峰值后，產(chǎn)沙過(guò)程又出現(xiàn)峰值。當(dāng)產(chǎn)流漸漸趨于穩(wěn)定，淺溝侵蝕產(chǎn)沙最終趨于穩(wěn)定。

圖8 26°時(shí)不同放水流量有無(wú)上方來(lái)水淺溝侵蝕產(chǎn)沙量隨時(shí)間的變化Fig.8 Sediment yield over time with or without up-slope runoff in different water discharge at slope of 26°

3.3 淺溝侵蝕產(chǎn)沙量與水動(dòng)力參數(shù)的關(guān)系

上方來(lái)水引起下方淺溝侵蝕產(chǎn)沙相對(duì)增量可用下式計(jì)算:

式中:ΔM為產(chǎn)沙相對(duì)增量，g/min;Msx為有上方來(lái)水時(shí)淺溝侵蝕產(chǎn)沙量，g/min;Mx為相對(duì)應(yīng)的無(wú)上方來(lái)水時(shí)下方淺溝侵蝕產(chǎn)沙量，g/min。剪切力和水流功率計(jì)算式同上，ΔP為水流功率相對(duì)增量，N/(m·s);ΔS為剪切力相對(duì)增量，Pa。

上方來(lái)水對(duì)坡下方淺溝侵蝕產(chǎn)沙量的影響與淺溝水力學(xué)參數(shù)有密切的關(guān)系，用相關(guān)分析法對(duì)實(shí)驗(yàn)資料進(jìn)行分析，得知上方來(lái)水引起的產(chǎn)沙量相對(duì)增量與淺溝水流的剪切力和水流功率相對(duì)增量的相關(guān)性較好，淺溝侵蝕產(chǎn)沙量相對(duì)增量與剪切力相對(duì)增量和水流功率相對(duì)增量均呈冪函數(shù)關(guān)系(表1)。

從圖9可以看出，上方來(lái)水引起的產(chǎn)沙相對(duì)增量與水流剪切力和水流功率相對(duì)增量呈冪函數(shù)關(guān)系，即隨著淺溝水流剪切力和水流功率的增大，侵蝕產(chǎn)沙量迅速增大。

表1 水流剪切力、水流功率與產(chǎn)沙量的關(guān)系Tab.1 Relationship between sediment yield and flow shearing force and stream power

4 結(jié)論與討論

1)上方來(lái)水的匯入使淺溝水流流速明顯增大，雷諾數(shù)、弗勞德數(shù)、水流功率和剪切力分別增大了33% ～76%、21% ～47%、29% ～72%、18% ～42%，阻力系數(shù)減少了11%～13%，導(dǎo)致淺溝侵蝕產(chǎn)沙量明顯增大。除流速和弗勞德數(shù)外，其余水動(dòng)力參數(shù)隨放水時(shí)間的延長(zhǎng)呈遞增趨勢(shì)。

2)上方來(lái)水對(duì)淺溝侵蝕產(chǎn)沙量貢獻(xiàn)很大，上方來(lái)水量大的，產(chǎn)沙增量亦愈大。與無(wú)上方來(lái)水情況相比，有上方來(lái)水的淺溝侵蝕產(chǎn)沙過(guò)程較為復(fù)雜，有多峰多谷的特點(diǎn)，其產(chǎn)沙過(guò)程大致為，隨著放水時(shí)間的延長(zhǎng)，產(chǎn)沙量在出現(xiàn)峰值后逐漸減少，最后趨于穩(wěn)定。

3)淺溝水流剪切力和功率是影響淺溝侵蝕產(chǎn)沙量的重要水動(dòng)力學(xué)參數(shù)。上方來(lái)水引起坡下方淺溝侵蝕產(chǎn)沙量相對(duì)增量與淺溝水流的剪切力和水流功率相對(duì)增量均呈冪函數(shù)關(guān)系。

上方來(lái)水來(lái)沙的匯入，加速了坡下方淺溝侵蝕的發(fā)展，并使淺溝侵蝕過(guò)程更為復(fù)雜。筆者在現(xiàn)有條件下做了一些有益的探索，而諸如不同含沙量的上方來(lái)水、不同下墊面影響因子以及復(fù)雜自然環(huán)境因子的影響在今后的工作中還需進(jìn)一步研究。

圖9 不同放水流量水流剪切力相對(duì)增量和水流功率相對(duì)增量與產(chǎn)沙量相對(duì)增量的關(guān)系Fig.9 Relationship between stream power，shearing force increment and sediment yield increment relatively in differentwater discharge

［1］許建民.黃土高原淺溝發(fā)育主要影響因素及其防治措施研究［J］.水土保持學(xué)報(bào)，2008，22(4):39-41

［2］劉志，江忠善.降雨因素和坡度對(duì)片蝕影響的研究［J］.水土保持通報(bào)，1994，14(2):19-22

［3］江忠善，王志強(qiáng)，劉志.黃土丘陵區(qū)站流域土壤侵蝕空間變化定量研究［J］.土壤侵蝕與水土保持學(xué)報(bào)，1996，2(1):1-9

［4］鄭粉莉，康紹忠.黃土坡面不同侵蝕帶侵蝕產(chǎn)沙關(guān)系及其機(jī)理［J］.地理學(xué)報(bào)，1998，53(5):422-428

［5］張科利.黃土坡面侵蝕產(chǎn)沙分配及其與降雨特征關(guān)系的研究［J］.泥沙研究，1991(4):39-47

［6］Zheng Fenli，Huang Chihua.Gully erosion［A］.In:Lal Rattan edited.Encyclopedia of Soil Science［M］.New York:Marcel Dekker，Inc，2002:630-634

［7］張科利，唐克麗，王斌科.黃土高原坡面淺溝侵蝕特征值的研究［J］.水土保持學(xué)報(bào)，1991，5(2):8-13

［8］Woodward D E.Method to predict crop and ephemeral gully erosion［J］.Catena，1999，37(3-4):393-399

［9］Nachtergaele J，Poesen J，Vandekerckhove L，et al.Testing the ephemeral gully erosion model for two mediterranean environments［J］.Earth Surface Processes and Landforms，2001，26(1):17-30

［10］鄭粉莉，唐克麗.降雨動(dòng)能對(duì)坡耕地細(xì)溝侵蝕影響的研究［J］.人民黃河，1995，17(7):22-24

［11］肖培青，鄭粉莉.上方來(lái)水來(lái)沙對(duì)細(xì)溝侵蝕產(chǎn)沙過(guò)程的影響［J］.水土保持通報(bào)，2001，21(1):23-25

［12］肖培青，鄭粉莉.細(xì)溝侵蝕過(guò)程與細(xì)溝水流水力學(xué)參數(shù)的關(guān)系研究［J］.水土保持學(xué)報(bào)，2001，15(1):54-57

［13］鄭粉莉，武敏，張玉斌，等.黃土陡坡裸露坡耕地淺溝發(fā)育過(guò)程研究［J］.地理科學(xué)，2006，26(4):438-442

［14］王文龍，雷阿林，李占斌，等.土壤侵蝕鏈內(nèi)細(xì)溝淺溝切溝動(dòng)力機(jī)制研究［J］.水科學(xué)進(jìn)展，2003，14(4):471-475

［15］鄭粉莉，高學(xué)田.黃土坡面土壤侵蝕過(guò)程與模擬［M］.西安:陜西人民出版社，2000:95-124

［16］陳浩.降雨特征和上方來(lái)水對(duì)產(chǎn)沙的綜合影響［J］.水土保持學(xué)報(bào)，1992，6(2):17-23