鄭加興佘冬立徐翠蘭劉冬冬沈 暉

（1.河海大學(xué)南方地區(qū)高效灌排與農(nóng)業(yè)水土環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，江蘇南京 210098；2.河海大學(xué)水資源高效利用與工程安全國(guó)家工程研究中心，江蘇南京 210098；3.江蘇省土地開發(fā)整理中心，江蘇南京 210017）

不同雨強(qiáng)和坡度條件海涂鹽土邊坡侵蝕細(xì)溝發(fā)育過程

鄭加興1，佘冬立1，2，徐翠蘭3，劉冬冬1，沈 暉1

（1.河海大學(xué)南方地區(qū)高效灌排與農(nóng)業(yè)水土環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，江蘇南京 210098；2.河海大學(xué)水資源高效利用與工程安全國(guó)家工程研究中心，江蘇南京 210098；3.江蘇省土地開發(fā)整理中心，江蘇南京 210017）

采用室內(nèi)人工模擬降雨試驗(yàn)，在不同坡度（11.3°、21.8°、35.0°）和降雨強(qiáng)度（85 mm／h、110mm／h、125mm／h）條件下，模擬海涂鹽土邊坡細(xì)溝發(fā)育過程。坡面侵蝕細(xì)溝發(fā)育動(dòng)態(tài)采用數(shù)碼攝像進(jìn)行監(jiān)測(cè)，通過Image?Pro Plus（IPP）6圖像分析軟件對(duì)圖像和數(shù)據(jù)進(jìn)行整理。結(jié)果表明：（a）不同坡度和降雨強(qiáng)度條件下，徑流產(chǎn)生后，很快在坡面形成細(xì)溝且細(xì)溝溝網(wǎng)密度變化較快；細(xì)溝數(shù)量、細(xì)溝深度、細(xì)溝平均寬度和細(xì)溝密度均隨降雨時(shí)間延長(zhǎng)而增大。（b）緩坡度條件下，土壤跌坑發(fā)育慢，細(xì)溝侵蝕率低；隨著坡度增大，坡面跌坑形成后水流聚集迅速，水流侵蝕沖刷力強(qiáng)，細(xì)溝發(fā)育快。（c）細(xì)溝發(fā)育越深，坡面細(xì)溝數(shù)量也越多；坡面侵蝕率與細(xì)溝平均深度及水流雷諾數(shù)分別在p＜0.05和p＜0.01水平上呈顯著的正線性相關(guān)關(guān)系。

海涂鹽土邊坡；邊坡坡度；降雨強(qiáng)度；細(xì)溝侵蝕；相關(guān)性分析

我國(guó)是世界土壤侵蝕最嚴(yán)重的國(guó)家之一，其土壤侵蝕范圍遍及全國(guó)各地，威脅社會(huì)和農(nóng)業(yè)可持續(xù)發(fā)展［1?2］。細(xì)溝侵蝕為土壤侵蝕的重要類型之一，是細(xì)溝溝岸和溝底土壤被細(xì)溝中的股流所分散、剝離和輸移的過程。侵蝕細(xì)溝的產(chǎn)生將加劇坡面侵蝕程度，細(xì)溝侵蝕量占坡面總侵蝕量的絕大部分［3?4］。細(xì)溝侵蝕發(fā)育過程及侵蝕產(chǎn)沙特征的研究對(duì)水土保持規(guī)劃、土壤侵蝕模型模擬等方面具有重要意義。

細(xì)溝侵蝕的發(fā)生取決于坡面水流水力學(xué)特性和坡面下墊面條件，包括降雨徑流條件、土壤性質(zhì)、地形因素、坡度、降雨雨強(qiáng)等。鄭粉莉等［5?7］、李君蘭等［8］研究了土壤侵蝕過程中包括細(xì)溝侵蝕在內(nèi)的不同侵蝕方式發(fā)生的臨界判定式，所構(gòu)建的侵蝕形態(tài)演變過程動(dòng)態(tài)監(jiān)測(cè)方法與模擬模型填補(bǔ)了該領(lǐng)域研究空白，并為流域侵蝕防治和預(yù)報(bào)模型建立提供了重要的科學(xué)依據(jù)。Gomez等［9］研究了土壤表面粗糙度對(duì)細(xì)溝發(fā)育過程的影響。作為降雨、徑流水動(dòng)力作用界面，土壤是決定侵蝕過程最重要的內(nèi)在因素。不同類型的土壤，其結(jié)構(gòu)穩(wěn)定性與土壤抵抗降雨濺蝕能力、降雨入滲以及土壤結(jié)皮形成、徑流的沖刷和泥沙的顆粒特性等侵蝕特性方面都存在著很大的差異，因此雨滴擊濺?徑流沖刷下坡面侵蝕細(xì)溝的發(fā)育過程及其水力學(xué)機(jī)理也有所差異。目前我國(guó)在這方面的研究主要集中于西北黃土、東北黑土、南方紅壤等，對(duì)于灘涂鹽堿地粉砂土的研究資料較少，極大限制了對(duì)該類型土壤細(xì)溝侵蝕過程的研究［10］。筆者在模擬降雨試驗(yàn)的基礎(chǔ)上，分析海涂圍墾區(qū)鹽堿土在坡度和降雨雨強(qiáng)雙因素條件下細(xì)溝的發(fā)育過程，以期為防治該區(qū)域水利工程邊坡侵蝕提供科學(xué)支持。

1 試驗(yàn)材料與方法

1.1 土壤材料

試驗(yàn)土樣取自江蘇省南通市如東縣東陵墾區(qū)（121°22′E，32°36′N）。土壤砂粒（粒徑0.25～0.05 mm）、粉粒（粒徑0.05～0.002mm）、和黏粒（粒徑小于0.002mm）的質(zhì)量分?jǐn)?shù)分別為38.6%、52.9%和8.5%，有機(jī)質(zhì)質(zhì)量比為3.26g／kg，鈉離子質(zhì)量比為1.6g／kg，交換性鈉的質(zhì)量分?jǐn)?shù)為68.9%，1∶5土樣和水經(jīng)過振蕩、離心后的電導(dǎo)率為5.98mS／cm。試驗(yàn)土樣經(jīng)自然風(fēng)干后過4mm篩，充分混合均勻。

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

模擬降雨試驗(yàn)采用組合側(cè)噴式單噴頭降雨器［11］，降雨噴頭組成包括噴頭體、碎流擋板和出流孔板部件。降雨器支架高度為4m，雨滴上噴高度為0.5m，雨滴到達(dá)水平地面的高度為4.5 m，雨滴降落速度可達(dá)到自然雨滴速度的95%以上，降雨均勻度為87%。降雨區(qū)采用2個(gè)噴頭互噴式，相隔距離5 m，形成疊加降雨區(qū)。每個(gè)噴頭裝置連接1根帶有控制閥門的供水管，供水采用中型泵取水（揚(yáng)程45m，流量6m3／s）。通過調(diào)節(jié)壓力，可獲得不同的降雨強(qiáng)度（本試驗(yàn)降雨強(qiáng)度分別為85 mm／h、110 mm／h、125 mm／h）。每次試驗(yàn)降雨持續(xù)時(shí)間為72min。試驗(yàn)用土槽采用Meyer［12］類似的設(shè)計(jì)。土槽長(zhǎng)50cm、寬30 cm、高15cm。土槽底部鋪設(shè)直徑2～10 mm的石子，厚度約為4.5 cm，裝土前再鋪1層紗布，保持平整，裝填土壤密度控制在1.3 g／cm3左右。土槽坡度設(shè)置為11.3°、21.8°、35.0°共3個(gè)梯度。每個(gè)處理進(jìn)行2次重復(fù)。

坡面侵蝕過程中細(xì)溝發(fā)育動(dòng)態(tài)采用數(shù)碼攝像記錄方式進(jìn)行監(jiān)測(cè)，通過Image?Pro Plus（IPP）6圖像分析軟件對(duì)圖像和數(shù)據(jù)進(jìn)行整理以得出細(xì)溝發(fā)育的動(dòng)態(tài)變化過程。本文主要選取各場(chǎng)降雨試驗(yàn)下坡面產(chǎn)生徑流后3min、6min和9min細(xì)溝發(fā)育過程進(jìn)行繪圖對(duì)比研究。同時(shí)，坡面開始產(chǎn)生徑流后每隔2 min取徑流、泥沙樣，觀察徑流、侵蝕動(dòng)態(tài)變化。采用高錳酸鉀染色劑法結(jié)合數(shù)碼攝像機(jī)錄像與IPP圖像處理分析方法，每隔2min分別測(cè)定細(xì)溝數(shù)量和細(xì)溝間的水流速度。試驗(yàn)過程中降雨水溫18～22℃。根據(jù)實(shí)測(cè)降雨強(qiáng)度資料與徑流資料，計(jì)算坡面水流動(dòng)力學(xué)參數(shù)［13］，包括水流雷諾數(shù)Re、弗勞德數(shù)Fr和曼寧糙率系數(shù)n。試驗(yàn)數(shù)據(jù)采用Excell和Spss12.0進(jìn)行分析處理。

2 結(jié)果分析

2.1 降雨強(qiáng)度對(duì)坡面細(xì)溝發(fā)育過程的影響

降雨強(qiáng)度對(duì)坡面細(xì)溝的形成和發(fā)育有著重要影響［14］。圖1呈現(xiàn)了坡度為35.0°情況下，降雨強(qiáng)度分別為85mm／h、110mm／h、125mm／h時(shí)各次降雨下鹽土坡面細(xì)溝形態(tài)發(fā)育過程。海涂圍墾區(qū)鹽土粉砂性強(qiáng)，降雨過程中坡面開始產(chǎn)流，細(xì)溝即開始發(fā)育，隨著產(chǎn)流時(shí)間的推進(jìn)，坡面土壤細(xì)溝數(shù)量增加，細(xì)溝寬度和深度不斷增大，之后趨于穩(wěn)定。對(duì)各次降雨試驗(yàn)下坡面侵蝕細(xì)溝數(shù)、寬度和深度發(fā)育的統(tǒng)計(jì)結(jié)果見表1。不同降雨強(qiáng)度下，海涂鹽土坡面細(xì)溝的發(fā)育主要體現(xiàn)在細(xì)溝數(shù)量和深度的發(fā)育。各徑流小區(qū)坡面經(jīng)濺蝕、表土入滲飽和等過程后，從坡面中下部開始產(chǎn)生徑流。在降雨水流的淘蝕作用下，坡面中下部逐漸出現(xiàn)跌坑。跌坑處由于積水形成一股水流，增加了水流的侵蝕力，同時(shí)土壤表面形態(tài)的變化遭到破環(huán)，使土壤抗侵蝕能力下降，直至徑流在該處突然下泄，跌坑不斷向坡面上部發(fā)展，坡面侵蝕細(xì)溝形成。隨著降雨過程的進(jìn)行，水流對(duì)細(xì)溝的溝底切割加深，侵蝕量增大。

坡度為35.0°、降雨強(qiáng)度分別為85 mm／h、110mm／h、125mm／h的坡面，開始產(chǎn)流的時(shí)間略有差異，分別為33min、29min和27 min。從各坡面產(chǎn)流后3 min細(xì)溝發(fā)育的圖像可以看出海涂_鹽土坡面降雨產(chǎn)流后不久就能形成細(xì)溝，此時(shí)細(xì)溝數(shù)量分別為121條、90條、63條，細(xì)溝平均寬度分別為0.31 cm、0.27 cm、0.34 cm。隨著降雨侵蝕的繼續(xù)進(jìn)行，坡面細(xì)溝溝壁穩(wěn)定性降低，加上此時(shí)溝壁土層含水量不斷增大，土塊逐漸失去平衡，坍塌至溝底，導(dǎo)致細(xì)溝寬度相應(yīng)地增加。徑流產(chǎn)生9 min的時(shí)候，細(xì)溝網(wǎng)已經(jīng)發(fā)育完全，此時(shí)各降雨強(qiáng)度下坡面細(xì)溝數(shù)量分別為137條、125條、81條，細(xì)溝平均寬度分別為0.46 cm、0.34 cm、0.51 cm。在降雨強(qiáng)度為 85 mm／h、110mm／h、125mm／h的情況下，細(xì)溝侵蝕密度分別為19.15%、27.03%、27.01%，較細(xì)溝發(fā)育初期（產(chǎn)流后3min）分別增加了1.6倍、2.7倍、3.0倍，細(xì)溝密度隨時(shí)間的變化也不斷增加。

2.2 坡度對(duì)細(xì)溝動(dòng)態(tài)發(fā)育過程的影響

如圖2所示，降雨強(qiáng)度為110 mm／h，坡度分別為11.3°、21.8°和35.0°時(shí)，細(xì)溝發(fā)育過程明顯。坡度為11.3°時(shí)，徑流產(chǎn)生初期，坡面濺蝕不明顯；跌坑產(chǎn)生后，水流需較長(zhǎng)時(shí)間在跌坑內(nèi)聚集，水流侵蝕力弱，跌坑向上發(fā)育變慢，延緩了細(xì)溝的形成；隨著坡度增大，坡面跌坑形成后水流聚集迅速，水流侵蝕沖刷力增強(qiáng)，細(xì)溝發(fā)育加快。

表2是降雨強(qiáng)度為110mm／h，坡度分別為11.3°、21.8°和35.0°情況下，徑流產(chǎn)生后細(xì)溝發(fā)育參數(shù)。隨著降雨時(shí)間的持續(xù)，細(xì)溝數(shù)量不斷增加；坡度越大，細(xì)溝數(shù)量越多，在徑流開始3 min時(shí)，坡度為35.0°的細(xì)溝數(shù)量是坡度為21.8°和11.3°的1.43倍和3.75倍；坡面產(chǎn)流后9 min，坡度為35.0°的細(xì)溝數(shù)量是坡度為21.8°和11.3°的1.14倍和1.79倍；平均細(xì)溝深度和寬度隨著時(shí)間的持續(xù)而增大，坡度為21.8°時(shí)，細(xì)溝最大平均深度分別是坡度在11.3°和35°下最大平均深度的13.1倍和1.16倍。

2.3 侵蝕參數(shù)與水力學(xué)參數(shù)的相關(guān)性分析

坡面細(xì)溝的發(fā)育是坡面水流沖刷與土壤抗蝕作用的綜合結(jié)果［15］。表3給出了部分試驗(yàn)參數(shù)。不同降雨強(qiáng)度和坡度條件下，徑流產(chǎn)生的時(shí)間并沒有顯著的規(guī)律性，這與土壤初始狀態(tài)有關(guān)。降雨強(qiáng)度相同，坡度變化的情況下，總產(chǎn)沙量隨時(shí)間變化明顯。在降雨強(qiáng)度為110 mm／h、坡度為35.0°時(shí)的總產(chǎn)沙量是坡度為21.8°和11.3°時(shí)總產(chǎn)沙量的2.31倍和9.97倍。坡度為11.3°和21.8°，降雨強(qiáng)度變化的情況下，總產(chǎn)沙量不斷增大。坡度為35.0°、降雨強(qiáng)度變化的情況下，總產(chǎn)沙量先增大、后減小。說明存在一個(gè)臨界坡度，在降雨強(qiáng)度變化的情況下總產(chǎn)沙量先增加到最大值，然后逐漸減小并達(dá)到穩(wěn)定。

為了研究水力參數(shù)與細(xì)溝發(fā)育過程的關(guān)系，試驗(yàn)取用坡面產(chǎn)生徑流后0～9 min相應(yīng)的試驗(yàn)數(shù)據(jù)，利用Spss軟件進(jìn)行了Pearson相關(guān)性檢驗(yàn)（雙側(cè)檢驗(yàn)），結(jié)果如表4所示。研究結(jié)果表明，細(xì)溝平均深度與細(xì)溝數(shù)量在0.01水平下兩兩顯著相關(guān)，說明細(xì)溝平均深度越深，坡面形成的細(xì)溝數(shù)量也就越多；細(xì)溝平均深度與坡面溝網(wǎng)密度、坡面侵蝕率在0.05水平下兩兩顯著相關(guān)；坡面侵蝕率與雷諾數(shù)在0.01水平下顯著相關(guān)且隨著坡面侵蝕率與水流雷諾數(shù)成線性相關(guān)（Rs＝1.9576Re-4.6432；R2＝0.71），這樣的結(jié)果可能緣于雷諾數(shù)的增加，水流沖刷形成的慣性力作用增強(qiáng)，產(chǎn)流量增大，產(chǎn)流挾帶的泥沙量增加，坡面侵蝕率也隨之增加；流速與弗勞德數(shù)在0.01水平下顯著相關(guān)。然而坡面侵蝕率與流速、弗勞德數(shù)、曼寧糙率系數(shù)等水力學(xué)參數(shù)關(guān)系不明顯，需要進(jìn)一步研究。

3 結(jié) 語

細(xì)態(tài)發(fā)育過程的研究對(duì)土壤侵蝕過程演變研究具有重要意義。細(xì)溝發(fā)育受土壤質(zhì)地、坡度、降雨強(qiáng)度等因素的影響。筆者采用人工模擬降雨，在不同坡度和降雨強(qiáng)度條件下對(duì)細(xì)溝發(fā)育過程和水力參數(shù)進(jìn)行分析，主要結(jié)論如下：

a.降雨強(qiáng)度對(duì)海涂鹽土坡面細(xì)溝動(dòng)態(tài)發(fā)育過程的影響明顯；不同降雨強(qiáng)度下，海涂鹽土坡面細(xì)溝的發(fā)育主要體現(xiàn)在細(xì)溝數(shù)量和深度的發(fā)育。

b.在不同坡度和降雨強(qiáng)度條件下，海涂鹽土坡面細(xì)溝數(shù)量、細(xì)溝深度、細(xì)溝平均寬度和細(xì)溝密度隨著降雨歷時(shí)的持續(xù)而不斷增加。

c.緩坡條件下，海涂鹽土土壤跌坑發(fā)育慢，細(xì)溝侵蝕率低；隨著坡度增大，坡面跌坑形成后水流聚集迅速，水流侵蝕沖刷力強(qiáng)，細(xì)溝發(fā)育快。

d.細(xì)溝發(fā)育越深，坡面細(xì)溝數(shù)量也越多；海涂鹽土坡面侵蝕率與細(xì)溝平均深度及水流雷諾數(shù)分別在p＜0.05和p＜0.01水平上呈顯著的正線性相關(guān)關(guān)系。

［1］SHE Dongli，F(xiàn)EI Yuanhang，LIU Zhipeng.Soil erosion characteristics on saline soil hillslope in coastal region of China：field investigation and rainfall simulation experiment［J］.Catena，2014，121：176?185.

［2］劉曉冰，周克勤，苗淑杰，等.土壤侵蝕影響作物產(chǎn)量及其因素分析［J］.土壤與作物，2012，1（4）：205?211.（LIU Xiaobing，ZHOU Keqin，MIAO Shujie，et al.Crop yield and relevant factors as affected by soil erosion［J］.Soil and Crop，2012，1（4）：205?211.（in Chinese）

［3］KIMARO D N，POESEN J，MSANYA B M，et al.Magnitude of soil erosion on the northern slope of the Uluguru Mountains，Tanzania：interrill and rill erosion［J］.Catena，2008，75：38?44.

［4］張晴雯，雷廷武，高佩玲，等.黃土區(qū)細(xì)溝侵蝕過程中輸沙能力確定的解析法［J］.中國(guó)農(nóng)業(yè)科學(xué)，2004，37（5）：699?703.（ZHANG Qingwen，LEI Tingwu，GAO Peiling，et al.An analytic method for determining transport capacity of concentrated flow in eroding rills［J］.Scientia Agricultura Sinica，2004，37（5）：699?703.（in Chinese））

［5］鄭粉莉，唐克麗，陳文亮.細(xì)溝侵蝕過程的研究方法［J］.中國(guó)科學(xué)院水利部西北水土保持研究所，1993，17（1）：107?111.（ZHENG Fenli，TANG Keli，CHEN Wenliang.Research methods of rill erosion process［J］.Memoir of NISWC，Academic Snica and Ministry of Water Resources，1993，17（1）：107?111.（in Chinese））

［6］鄭粉莉.發(fā)生細(xì)溝侵蝕的臨界坡長(zhǎng)與坡度［J］.中國(guó)水土保持，1989（8）：23?24.（Zheng Fenli.The critical slope length and slope when rill erosion happens［J］.Soil and Water Conservation in China，1989（8）：19?24.（in Chinese））

［7］鄭粉莉，唐克麗，周佩華.坡耕地細(xì)溝侵蝕的發(fā)生發(fā)展和防治途徑的探討［J］.水土保持學(xué)報(bào)，1987，1（1）：36?48.（Zheng Fenli，Tang Keli，Zhou Peihua.Approach to the genesis and development of roll erosion on slope land and the way to control［J］.Journal of Soil and Water Conservation，1987，1（1）：36?48.（in Chinese））

［8］李君蘭，蔡強(qiáng)國(guó)，孫莉英，等.降雨強(qiáng)度、坡度及坡長(zhǎng)對(duì)細(xì)溝侵蝕的交互效應(yīng)分析［J］.中國(guó)水土保持科學(xué)，2011，9（6）：8?13.（LI Junlan，CAI Qiangguo，SUN Liying，et al.Analysis of interaction effects of rainfall intensity，slope degree and slope length on rill erosion［J］.Science of Soil and Water Conservation，2011，9（6）：8?13.（in Chinese））

［9］GOMEZ J A，DARBOUX F，NEARING M A.Development and evolution of rill network under simulated rainfall［J］.Water Resource Research，2003，1148：1029?2002.

［10］LIU Dongdong，SHE Dongli，YU Shuang'en，et al.Rainfall intensity and slope gradient effects on sediment losses and splash from a saline?sodic soil under coastal reclamation［J］.Catena，2015，128：54?62.

［11］陳文亮.組合側(cè)噴式野外人工模擬降雨裝置［J］.水土保持通報(bào)，1984（5）：43?48.（CHEN Wenliang.Artificial simulate rainfall device of combination of side jet［J］.Bulletin of Soil and Water Conservation，1984（5）：43?48.（in Chinese））

［12］GUO Tailong，WANG Quanjiu，LI Dingqiang，et al.Flow hydraulic characteristic effect on sediment and solute transport on slope erosion［J］.Catena，2013，107：145?153.

［13］佘冬立，劉冬冬，彭世彰，等.海涂圍墾區(qū)排灌工程邊坡土壤侵蝕過程的水動(dòng)力學(xué)特征［J］.水土保持學(xué)報(bào)，2014，28（1）：1?5.（She Dongli，Liu Dongdong，Peng Shizhang，et al.Hydrodynamic characters of runoff in soil erosion processes on drainage?irrigation engineering slope of farmland in coastal region［J］.Journal of Soil and Water Conservation，2014，28（1）：1?5.（in Chinese））

［14］BERGER C，SCHULZE M，RIEKE D，et al.Rill development and soil erosion：a laboratory study of slope and rainfall intensity［J］.Earth Surface Processes and Landforms，2010，35：1456?1467.

［15］張永東，吳淑芳，馮浩，等.黃土陡坡細(xì)溝侵蝕動(dòng)態(tài)發(fā)育過程及其發(fā)生臨界動(dòng)力條件試驗(yàn)研究［J］.泥沙研究，2013（2）：25?32.（ZHANG Yongdong，WU Shufang，F(xiàn)ENG Hao，et al.Experimental study of rill dynamic development process and its critical dynamic conditions on loess slope［J］.Journal of Sediment Research，2013（2）：25?32.（in Chinese））

Evolution process of erosion rill on saline soil slope in coastal reclamation with different rainfall intensities and slope gradients

ZHENG Jiaxing1，SHE Dongli1，2，XU Cuilan3，LIU Dongdong1，SHEN Hui1
（1.Key Laboratory of Efficient Irrigation?Drainage and Agricultural Soil?Water Environment in Southern China，Ministry of Education，Hohai University，Nanjing 210098，China；2.National Engineering Research Center of Water Resources Efficient Utilization and Engineering Safety，Hohai University，Nanjing 210098，China；3.Jiangsu Land Consolidation and Rehabilitation Center，Nanjing 210017，China）

Simulated rainfall experiments were conducted with three slope gradients（11.3°，21.8°，and 35.0°）and three rainfall intensities（85 mm／h，110 mm／h，and 125 mm／h）to investigate the rill evolution process on saline soil slopes in coastal reclamation.The evolution process of rill erosion was recorded with a digital camera，and then the image data were analyzed with Image?Pro Plus（IPP）6 software.The results are as follows：（a）Rills formed soon after the runoff occurred and the rill density changed rapidly with different slope gradients and rainfall intensities.The rill number，rill depth，average rill width，and rill density increased with the rainfall duration.（b）At low slope gradients the drop pits developed slowly，and the rill erosion rate was low.With the increase of the slope gradient，water concentrated in the drop pits rapidly，and the erosion rate increased，leading to rapid evolution of rills.（c）The number of rills increased with the rill depth.There were significant linear positive correlations of the slope erosion rate with the average depth of rills andReatp＜0.05 andp＜0.01，respectively.

saline soil slope in coastal reclamation；slope gradient；rainfall intensity；rill erosion；correlation analysis

S157

A

1000-1980（2015）04-0313-06

10.3876／j.issn.1000-1980.2015.04.006

2014-0930

國(guó)家自然科學(xué)基金（51109063，41471180）

鄭加興（1989—），男，安徽池州人，碩士研究生，主要從事水資源規(guī)劃研究。E?mail：846766720＠qq.com

佘冬立，教授。E?mail：shedongli＠hhu.edu.cn