,

(中國(guó)科學(xué)院 東北地理與農(nóng)業(yè)生態(tài)研究所 黑土區(qū)農(nóng)業(yè)生態(tài)重點(diǎn)實(shí)驗(yàn)室，黑龍江 哈爾濱 150081)

0 引 言

溝道侵蝕是短期內(nèi)積累的徑流水在狹窄的渠道經(jīng)常反復(fù)沖刷，把土壤剝離成具有一定深度的狹窄區(qū)域的侵蝕過(guò)程[1]。其發(fā)生發(fā)展過(guò)程對(duì)現(xiàn)代地貌發(fā)育及演化過(guò)程有著重要的影響，從高地向低谷產(chǎn)生的徑流和產(chǎn)沙與溝道發(fā)生發(fā)展過(guò)程密切相關(guān)，且加劇了水蝕產(chǎn)生的非點(diǎn)源污染程度。

溝道侵蝕是人類(lèi)活動(dòng)誘發(fā)的土地變化和極端降雨共同作用的結(jié)果，是長(zhǎng)期存在的一個(gè)最受關(guān)注的具有普遍性、緊迫性和復(fù)雜性的世界性環(huán)境問(wèn)題[2]。它是產(chǎn)沙的主要來(lái)源[3-4]，也是土地退化的主要過(guò)程及最嚴(yán)重的表現(xiàn)形式[2,5]。

溝道侵蝕研究，如果從Charles Lyell的溝道描述算起至今170多年[6]，從Rubey[7]的工作算起90年，從Ireland的工作算79年[6]。目前，對(duì)溝道侵蝕的研究已經(jīng)進(jìn)入成熟期。

除了極地氣候外，溝道侵蝕在其它氣候條件下均有發(fā)生，多數(shù)情況下與不合理的土地管理方式直接相關(guān)[8]。過(guò)去的1 500年，溝道侵蝕在世界范圍內(nèi)都有報(bào)道。如，美國(guó)西南部水沖溝的形成就是由于過(guò)度放牧和氣候變化而致[9]；英國(guó)第九世紀(jì)和第十世紀(jì)溝道的形成是由人類(lèi)活動(dòng)誘發(fā)植被變化所導(dǎo)致的匯水區(qū)水文過(guò)程的變化引起的[10]；第十四世紀(jì)德國(guó)的溝道形成是因土地的高強(qiáng)度利用和極端降雨而致[11]。

我國(guó)東北黑土區(qū)，在過(guò)去的百余年間土地被大規(guī)模、高強(qiáng)度的快速開(kāi)發(fā)利用，導(dǎo)致水土流失現(xiàn)象嚴(yán)重[12-13]?，F(xiàn)有長(zhǎng)度大于100 m的侵蝕溝數(shù)量為29.6萬(wàn)條，損毀耕地面積約占該區(qū)域的0.5%，且侵蝕面積和強(qiáng)度仍在不斷增加，侵蝕溝數(shù)量和密度也隨之增加，抵御自然災(zāi)害能力明顯下降[14]。東北黑土區(qū)已成為我國(guó)水土流失的3個(gè)急需治理的區(qū)域之一。

溝道形成過(guò)程具有三維本質(zhì)，受到各種因素和過(guò)程的影響，不僅研究難度大，而且很難預(yù)測(cè)。盡管多數(shù)的溝道侵蝕是由土地利用變化或是極端氣候所誘導(dǎo)或加速的，但其結(jié)果是長(zhǎng)期的歷史過(guò)程所形成的[15]。溝道侵蝕一旦發(fā)生，發(fā)展的維度和景觀(guān)的連通性迅速增加，導(dǎo)致沉積物向下坡沉積的危險(xiǎn)性增加，影響土壤的其它退化過(guò)程[16]。農(nóng)業(yè)用地的水道和土壤特性的許多沉積及化學(xué)的危害都與溝道有關(guān)。溝道形成后，很難實(shí)施有效的防治技術(shù)，而且投入成本很高。因此，從全球的角度而言，溝道侵蝕作為土地退化的主要驅(qū)動(dòng)力不容忽略。

本文旨在從國(guó)內(nèi)外溝道侵蝕的分類(lèi)、多樣性和發(fā)生過(guò)程、影響因素及其臨界條件方面進(jìn)行評(píng)述，并對(duì)今后需要探討的問(wèn)題予以展望。

1 溝道侵蝕的多樣性

溝道一詞，在西班牙稱(chēng)為“carvaka”，法國(guó)“ravine”，馬達(dá)加斯加“l(fā)avaka”，阿拉伯“wadi”，南非“donga”，巴西“vocoroca”，阿根廷“barranco”。根據(jù)最古老的西班牙詞典，“carvaka”一詞源于‘肉’‘car’和‘洞穴’‘vava’的組合。其詞源與戰(zhàn)爭(zhēng)之后埋葬殉職人員的普通深坑有關(guān)。字典同時(shí)定義溝道為由于極端降雨對(duì)易侵蝕土地作用形成的大切口。

不同學(xué)科定義溝道的標(biāo)準(zhǔn)不同。形態(tài)學(xué)和地貌學(xué)的標(biāo)準(zhǔn)是：景觀(guān)中匯水區(qū)面積為10 km2或小于10 km2，有一定深度，溝壁陡峭，且少有植被的切溝[17]。水文學(xué)的標(biāo)準(zhǔn)是：暴雨過(guò)程中濺蝕洪水形成的臨時(shí)水線(xiàn)[18]。農(nóng)業(yè)措施標(biāo)準(zhǔn)是：水流溝的寬度和深度不能進(jìn)行正常的土壤耕作[19]。不穩(wěn)定性標(biāo)準(zhǔn)是：流域內(nèi)以前不存在確定的溝渠，而近期形成的切溝[20]。

Davis[21]則將溝道定義為平坦廣袤的平原上那些原有的正常河道處于儲(chǔ)水能力持續(xù)降低階段的溝渠。Ireland等[6]首次提出了溝道形成理論，強(qiáng)調(diào)了溝道作為一個(gè)重要現(xiàn)象的可能性，并第一次提出Bt-層抗性在控制溝道深度和溝頭形狀中所起的重要作用。

由于溝道形態(tài)是溝道過(guò)程的產(chǎn)物，而且溝道的形態(tài)代表著溝道過(guò)程的過(guò)去、現(xiàn)在和未來(lái)，因此，明確溝道的形態(tài)學(xué)是評(píng)估溝道過(guò)程的第一步。Ireland等[6]基于在美國(guó)南加州的研究結(jié)果，首次對(duì)溝道進(jìn)行了分類(lèi)，并對(duì)溝道的形態(tài)和樣式進(jìn)行了簡(jiǎn)單的描述，他們提出了溝道的6種形態(tài)：線(xiàn)性的、球莖的、樹(shù)枝狀的、格狀的、并行的和復(fù)合狀的。

因此，溝道有多種形式，但多數(shù)是指丘陵坡溝道(Hillslope gullies)。一般而言，溝道比細(xì)溝大，但比河道小[22]。典型的丘陵坡溝道長(zhǎng)度遠(yuǎn)大于寬度，具有一個(gè)臨界閾值，即坡度和匯水面積的函數(shù)。當(dāng)剪切應(yīng)力超過(guò)特定值時(shí)即可發(fā)生[23]。由于最大的滲透速率不會(huì)超過(guò)最大的降雨強(qiáng)度，由此主要的徑流過(guò)程在匯水區(qū)和暴雨之間差異較大[24]。

國(guó)外侵蝕溝分類(lèi)標(biāo)準(zhǔn)存在一定差異，F(xiàn)oster[25]在1986年根據(jù)侵蝕溝的持久性將其分為臨時(shí)性溝道(Ephemeral gully)和永久性溝道(Permanent gully)，認(rèn)為臨時(shí)溝道每次侵蝕的寬度和深度大于細(xì)溝，但小于永久溝道，不妨礙耕作，但不能消除其痕跡。隨著年復(fù)一年耕作與侵蝕的交替，臨時(shí)切溝可發(fā)展為永久溝道。

Poesen等[26]將侵蝕溝道分為3種類(lèi)型：永久溝道、岸坡溝道和臨時(shí)溝道。

永久性溝道被定義為：那些農(nóng)業(yè)土地中存在的溝渠深度較深(典型的深度為0.5 m以上，有的深度達(dá)20～30 m)，以至于不能進(jìn)行常規(guī)的農(nóng)機(jī)具耕作，且不能去除的渠道[1]。換句話(huà)說(shuō)，永久性溝道就是具有橫截面且寬度和深度正常耕作農(nóng)機(jī)具不能實(shí)現(xiàn)橫跨的渠道。該定義強(qiáng)制性提出用最低深度0.5 m來(lái)區(qū)分于細(xì)溝。永久溝道通常是大雨期間或雨后，由于地表徑流和壤中流形成的沖積、崩積沉積物的切力而形成的[27-28]?？臻g上，永久溝道可以是不連續(xù)的。Morgan[18]根據(jù)永久性溝道所在的空間位置將其分為谷底侵蝕溝和谷壁侵蝕溝。

岸坡溝道是當(dāng)集中徑流橫跨土堤時(shí)形成的。由于陡峭的坡度梯度，岸坡溝道可以通過(guò)水力侵蝕、管涌侵蝕和最終的土體移動(dòng)快速形成。一旦形成，通過(guò)溝頭挺進(jìn)遷移到中度坡土壤表面，并進(jìn)入江河或農(nóng)業(yè)梯田[27]。岸坡溝道的形成不僅需要很長(zhǎng)的時(shí)間，還很難控制[29]。

由年際間的正常耕作所填平的擴(kuò)展了的細(xì)溝稱(chēng)之為臨時(shí)溝道[30]。即臨時(shí)溝道是比細(xì)溝大，但比永久溝道小的渠道切口，常規(guī)耕作能夠?qū)⑵淙コ?，但是額外的徑流發(fā)生在同一地點(diǎn)又會(huì)形成溝道[31-33]。

臨時(shí)溝道是主要的產(chǎn)沙來(lái)源，占土壤侵蝕產(chǎn)沙量的20%～100%[33]，它是農(nóng)業(yè)耕地中最為嚴(yán)重的侵蝕類(lèi)型，在土壤未被植被覆蓋的種床準(zhǔn)備階段、作物器官形成和種植過(guò)程中均可發(fā)生。臨時(shí)溝道會(huì)引起嚴(yán)重的土壤流失現(xiàn)象，導(dǎo)致土壤質(zhì)量逐漸退化和土壤潛在生產(chǎn)力下降[34-38]。

臨時(shí)溝道形成區(qū)作物被沖刷走，并被沉積物掩埋。填埋運(yùn)作會(huì)引起土壤退化，降低農(nóng)田的長(zhǎng)期生產(chǎn)力[39-40]。而且，臨時(shí)溝道成為高地向谷底運(yùn)輸徑流和沉積物的有效渠道，并很快演化為永久溝道[41-42]。

基于室外觀(guān)測(cè)，Casalí等[31]描述了臨時(shí)溝道的3種主要類(lèi)型：經(jīng)典的、排水的和不連續(xù)的臨時(shí)溝道。經(jīng)典的臨時(shí)溝道是由相同地塊的徑流發(fā)生區(qū)的集中徑流形成，而排水的溝道來(lái)自于地塊上游區(qū)的集中排水流形成，排水流抵達(dá)田塊的上游區(qū)，侵蝕下游區(qū)的農(nóng)田。不連續(xù)的溝道常常發(fā)生在由管理措施所形成的坡度突變的區(qū)域，如與道路相鄰的農(nóng)田。臨時(shí)溝道的溝頭或許就是由這些坡度的不連續(xù)性引起的，之后向上游移動(dòng)。

臨時(shí)性溝道和永久性溝道侵蝕的主要差異是前者深度淺(通常深度只有0.5 m)，且存在時(shí)間短，每年的耕作能夠?qū)⑵涮盥馵26]。

國(guó)內(nèi)最早將侵蝕溝分為淺溝、切溝、沖溝和干溝等[43]，后來(lái)又認(rèn)為完整的溝蝕序列一般包括細(xì)溝、淺溝、切溝和沖溝，認(rèn)為淺溝由主細(xì)溝演變而來(lái)，并能發(fā)展為切溝。劉元保等[44]認(rèn)為，淺溝侵蝕是暴雨發(fā)生時(shí)匯集于集流槽底部，由徑流沖刷形成的新的侵蝕溝槽過(guò)程中所造成的土壤侵蝕，即發(fā)生在坡面順坡集流槽的底部(寬度一般在2 m 以?xún)?nèi))。淺溝一般不阻礙普通耕作，但犁耕不能消除淺溝形態(tài)痕跡[45]。鄭粉莉等[46]認(rèn)為細(xì)溝屬于面蝕的范疇，而淺溝、切溝和沖溝是侵蝕溝發(fā)育的3個(gè)不同階段。

中國(guó)農(nóng)業(yè)百科全書(shū)土壤卷編委會(huì)在《土壤侵蝕與水土保持分支條目》[47]中指出，淺溝由坡耕地上主細(xì)溝發(fā)展形成，其橫斷面因不斷的再侵蝕和再耕作呈弧形擴(kuò)展，無(wú)明顯的溝緣。在黃土丘陵區(qū)淺溝的深度均大于耕層厚度，一般為20～30 cm，也有超過(guò)50 cm的；而在南方花崗巖風(fēng)化殼丘陵斜坡上，也可發(fā)生淺溝，其寬度約1 m 左右，深度可超過(guò)0.5 m。

鄭粉莉等[48]在參編《Encyclopedia of Soil Science》溝蝕(Gully)條目時(shí)，將我國(guó)定義的淺溝侵蝕與美國(guó)定義的臨時(shí)性侵蝕(Ephemeral gully)歸為同一類(lèi)溝蝕類(lèi)型。

由于定義的多樣性和溝道侵蝕評(píng)估的復(fù)雜性，人們不斷地對(duì)此問(wèn)題進(jìn)行重復(fù)探索，認(rèn)知也在不斷深入[49-51]。

2 溝道侵蝕發(fā)生過(guò)程

Ireland等[22]認(rèn)為，溝道是由下坡的沖刷、溝頭的切割及其快速擴(kuò)展、自我愈合和趨于穩(wěn)定的過(guò)程。Horton[52]關(guān)于溝道形成的經(jīng)典論文指出，當(dāng)降雨強(qiáng)度超過(guò)入滲能力時(shí)，水就會(huì)在洼地積累，溢出的水以不規(guī)則湍流向下坡流動(dòng)，對(duì)土壤施加剪切應(yīng)力。分水線(xiàn)周邊的剪切應(yīng)力小于夾帶的閾值，而下坡的水流水位增高，超過(guò)閾值，便形成了“非侵蝕帶”，這個(gè)“非侵蝕帶”是降雨強(qiáng)度與入滲能力的差異、水力粗糙度、表面臨界牽引力和所在坡度的函數(shù)?！胺乔治g帶”以外，片蝕與剪切應(yīng)力相稱(chēng)，由此形成流水的深度和水表面的梯度，即侵蝕溝。

FAO提出，溝道的形成發(fā)育是幾個(gè)過(guò)程單一或者同時(shí)作用的結(jié)果，包括(1)流水和磨料(土壤或瓦礫)對(duì)底部或溝道壁的沖刷；(2)溝頭(跌水效應(yīng))瀑布式的侵蝕。溝道發(fā)育的速率和程度與徑流量和徑流速度緊密相關(guān)。溝道形成需要相對(duì)大量的水來(lái)提供用于分離和運(yùn)送土壤的能量。水量與匯水區(qū)面積的大小和徑流特性緊密相連，暴雨徑流是搬運(yùn)溝道瓦礫的主要外力[53]。因此，溝道主要通過(guò)沖積而形成，即沖積過(guò)程是引起溝道形成的初始因素。

實(shí)際上，溝道形成是一個(gè)復(fù)雜的地貌系統(tǒng)過(guò)程，而不是簡(jiǎn)單的沖積過(guò)程，它的形成涉及到溝頭、溝岸、溝底、溝頭上面的坡地及匯水面積等。Bryan和Yair[54]指出，形成溝道的主要過(guò)程包括雨濺、表面徑流、面蝕、集中徑流和管道徑流。Hudson[28]認(rèn)為，溝道的形成是水道過(guò)程和形態(tài)平衡被破壞的結(jié)果，(1)要么是渠道承載的洪水徑流量的增加；(2)要么是渠道承載水流的能力下降。土地利用的變化或匯水面積的增加都會(huì)導(dǎo)致水流的增加；渠道的粗燥度增加或是沉積后水流量降低都會(huì)降低渠道承載水的能力。由此，Hudson[28]指出溝道形成的過(guò)程是不能自我修復(fù)的。他的這種說(shuō)法，與溝道可能隨著時(shí)間的推移會(huì)逐漸穩(wěn)定的提法形成鮮明對(duì)比。

景可[55]認(rèn)為切溝侵蝕發(fā)展的方式主要有：(1)溝頭的溯源侵蝕，溝頭前進(jìn)是切溝發(fā)展的主要方式之一；(2)溝坡橫向侵蝕的發(fā)展，主要方式是瀉溜、崩塌和滑坡；(3)垂向的下切侵蝕，在其它條件相同的情況下，切溝下切侵蝕的強(qiáng)度取決于水流侵蝕強(qiáng)度。

Bocco[56]認(rèn)為潛蝕是切溝形成和發(fā)展的最重要形式之一。Roberson 和Hanson[57]認(rèn)為溯源侵蝕是切溝發(fā)展的主要方式之一。Majid[58]通過(guò)對(duì)切溝侵蝕過(guò)程的監(jiān)測(cè)結(jié)果表明，切溝的發(fā)展是一系列過(guò)程作用的結(jié)果，其中包括滲漏、土壤蠕動(dòng)、塌方和沖刷。沿切溝發(fā)生的侵蝕形態(tài)和空間分布規(guī)律表明，滲漏和地表徑流這2 種水文過(guò)程在切溝發(fā)展過(guò)程中起著非常重要的作用。而且，溝頭下切侵蝕的主要過(guò)程是滲漏，地表徑流則是維持溝頭溯源侵蝕的必要過(guò)程。

于國(guó)強(qiáng)等[59]對(duì)黃土高原坡溝系統(tǒng)的研究表明，坡溝系統(tǒng)上部位移以“沉降”模式為主，下部位移以“剪切”模式為主，溝頭溯源區(qū)是坡溝系統(tǒng)重力侵蝕最為強(qiáng)烈的部位；剪切塑性區(qū)域主要分布于坡面和溝坡大部分區(qū)域，張拉塑性區(qū)域主要分布于梁峁頂和梁峁坡上部。

鄭粉莉等[46]結(jié)合她們團(tuán)隊(duì)研究的成果，提出溝道侵蝕包括溝頭溯源侵蝕、溝壁崩塌和溝底下切3 個(gè)子過(guò)程。溯源侵蝕會(huì)影響匯水面積，改變坡面坡度。因此,在當(dāng)前溝頭處測(cè)量的坡度及匯水面積并不能代表溝道形成之初的情況。

3 溝道侵蝕影響因素

總體上，影響溝道侵蝕的因素有侵蝕動(dòng)力因子(降雨、上方匯流及地下潛流)、地形因子(坡度、坡長(zhǎng)及坡形等)、地表植被、土壤、土地類(lèi)型和人類(lèi)活動(dòng)等[46]。Watson等[30]曾指出，參與溝道形成的因素包括3個(gè)方面：(1)徑流的體積、速度和類(lèi)型；(2)土壤內(nèi)含物對(duì)侵蝕的敏感性或溝道的侵蝕性；(3)土地利用和保護(hù)措施引起的植被覆蓋的變化。因此，影響溝道侵蝕的因素可以簡(jiǎn)單歸結(jié)為：能量、抗性和保護(hù)措施。

眾所周知，雨滴擊濺和徑流沖刷是水土流失的動(dòng)力[60-61]。而極端的降雨，尤其當(dāng)降雨的能量與地表徑流的合力超過(guò)土壤顆粒耐分離的能力時(shí)，就會(huì)發(fā)生溝蝕[62]。

降雨對(duì)溝道侵蝕的影響主要體現(xiàn)在降雨強(qiáng)度、降雨量等方面。淺溝侵蝕主要受降雨強(qiáng)度影響，基本與降雨量無(wú)關(guān)。當(dāng)降雨強(qiáng)度由1.31 mm·min-1增大到3.52 mm·min-1時(shí)，在降雨量大體相當(dāng)?shù)那闆r下，每毫米降雨引起的侵蝕量由0.543 kg增加到2.626 kg。當(dāng)降雨強(qiáng)度增大2.687倍時(shí)，淺溝侵蝕量增加4.836倍[63]。意大利西西里島典型小流域內(nèi)的研究結(jié)果表明，當(dāng)3 d 最大降雨量達(dá)到51 mm 時(shí)，流域內(nèi)即開(kāi)始發(fā)生淺溝侵蝕；雖然一年內(nèi)平均監(jiān)測(cè)到7 次侵蝕性降雨，但淺溝的發(fā)生和發(fā)育均是由單次侵蝕性暴雨所導(dǎo)致的[64]。此外，降雨在年內(nèi)的分布特征、降雨雨型等對(duì)淺溝發(fā)育過(guò)程也有重要影響。

次降雨條件下，坡上方來(lái)水使淺溝侵蝕帶的侵蝕量增加12%～84%，年平均增加38%～66%[65]。而且，上方來(lái)水引起的侵蝕產(chǎn)沙量隨降雨強(qiáng)度的增加而增大，其增加幅度為4.4%～83.5%；1L·min-1的上方匯流可引起單位長(zhǎng)度淺溝剝離率達(dá)0.1～0.3 kg·(m·min)-1；上方來(lái)水使淺溝溝坡徑流流速較無(wú)上方來(lái)水時(shí)增加12%～24%，尤其是上方來(lái)水使淺溝溝槽徑流流速增加45.6%～58.4%。上方來(lái)水時(shí)的淺溝土槽的徑流含沙量較無(wú)上方來(lái)水時(shí)的含沙量增大5.4%～287%[66]。由此提出，溝道上方來(lái)水來(lái)沙對(duì)淺溝侵蝕帶產(chǎn)沙量具有重要影響。

地下徑流的潛蝕影響著坡面淺溝侵蝕的發(fā)生和發(fā)展過(guò)程，但通過(guò)土壤大孔隙或土壤管道運(yùn)輸?shù)牡叵聺摿骱茈y被觀(guān)測(cè)到，常常被研究人員所忽視[67]。Wilson等[35]研究表明，地下潛流的存在加快了淺溝的發(fā)生過(guò)程，隨著徑流的沖刷，土壤孔隙慢慢變大，表層土壤發(fā)生崩塌，產(chǎn)生的侵蝕量也變大，崩塌的土壤顆粒被徑流搬運(yùn)后便形成了淺溝，這一過(guò)程受土壤性質(zhì)和土地利用等多因素的影響。由此認(rèn)為，針對(duì)地表徑流的一般水保措施，很難對(duì)地下潛流產(chǎn)生的侵蝕進(jìn)行防治。

存在于地表的黃土、泥灰?guī)r以及其他沉積物，由于其相對(duì)豐度高且抗侵蝕能力低，是最普遍的巖性因素。風(fēng)成沉積是世界所有黃土帶最易形成溝道侵蝕的原因，在美國(guó)、以色列、阿根廷、波蘭和比利時(shí)均有發(fā)生[68-70]。

黃土區(qū)由溝道侵蝕所引起的生態(tài)環(huán)境脆弱性，主要源于其對(duì)土壤顆粒分離的抗性低，對(duì)瓦解倒塌的敏感性高以及不能滲透但促進(jìn)滲漏的地質(zhì)層次配置等因素。泥灰?guī)r對(duì)侵蝕的敏感性，主要源于它們的黏性結(jié)構(gòu)、柔軟度和地球化學(xué)特性，如具有較高的pH值和鈉離子含量[71-74]。韓魯艷等[75]發(fā)現(xiàn)，黃土侵蝕量和水穩(wěn)性團(tuán)聚類(lèi)因子及有機(jī)質(zhì)含量之間分別存在極顯著與顯著相關(guān)關(guān)系，鑒于土壤分析的易獲性，由此提出可選擇>0.5 mm水穩(wěn)性團(tuán)聚體與有機(jī)質(zhì)含量作為反映土壤侵蝕程度的指標(biāo)。

疏松的砂巖、泥巖和頁(yè)巖經(jīng)常在溝道中出現(xiàn)[76-78]。除了特定的基巖屬性外，有些沉積層次配置，如黃土耕作層的配置或頁(yè)巖沙石交替層[79]，能夠促進(jìn)形成有利于滲漏的條件，易導(dǎo)致壁體塌陷。然而，溝道侵蝕在火成巖和變質(zhì)巖巖性區(qū)也同樣發(fā)生，但比例較小[80]。這些巖石堅(jiān)硬固相的本質(zhì)解釋了為什么嚴(yán)重溝道侵蝕經(jīng)常發(fā)生在深層土、風(fēng)化層或腐泥土層環(huán)境中[81-82]。

研究發(fā)現(xiàn)地下巖石的非連續(xù)性控制著溝道的形成[83-86]，基巖對(duì)切口的極限抗性、溝道橫切面積的維度以及景觀(guān)巖石學(xué)意義的溝道分布均受其控制[87-88]。

地形因子對(duì)淺溝侵蝕的影響主要體現(xiàn)在坡度、坡長(zhǎng)和坡形等特征方面。

坡度影響坡面的受雨面積及受雨量，從而影響坡面徑流、入滲能力和徑流動(dòng)能的大小。淺溝發(fā)生的臨界坡度是臨界動(dòng)能的具體體現(xiàn)之一。黃土丘陵區(qū)淺溝發(fā)生的臨界坡度為18.2°，坡度越大，坡面物質(zhì)重力在沿坡面向下方向的分量也就越大，穩(wěn)定性就越差、越容易發(fā)生淺溝侵蝕[88]。劉元保等[44]基于野外調(diào)查發(fā)現(xiàn)，淺溝頂端到分水嶺的距離以及坡面順坡集流槽的間距與坡度呈線(xiàn)性相關(guān)關(guān)系。

坡長(zhǎng)影響徑流的匯集過(guò)程及徑流量的多少，決定著淺溝侵蝕的發(fā)生。就某一條淺溝而言，在一定坡長(zhǎng)范圍內(nèi)侵蝕量隨坡長(zhǎng)的增加而增加，但超過(guò)一定的坡長(zhǎng)范圍后，由于泥沙負(fù)荷的增加，徑流挾沙力減小，淺溝侵蝕量減小。黃土丘陵區(qū)淺溝發(fā)生臨界坡長(zhǎng)特征值范圍為20～60 m，平均為40 m 左右。由于淺溝已有了固定形態(tài)，淺溝侵蝕發(fā)生的臨界匯水面積就等于發(fā)生淺溝侵蝕的臨界坡長(zhǎng)與間距的乘積，該特征值為300～1 200 m2，以400～800 m2之間居多，平均為657 m2[63]。

坡形是坡度和坡長(zhǎng)的組合形態(tài)，其決定著徑流的匯集方式和過(guò)程，進(jìn)而影響淺溝的分布形式和分布密度，由此影響著侵蝕量。唐克麗等[89]研究表明，凸凹形坡與凹形坡上的土壤流失量均高于直形坡的原因是凸凹形坡與凹形坡易發(fā)生淺溝侵蝕，而直形坡一般尚無(wú)淺溝形成，在其它條件相同的情況下，凸凹形地塊的侵蝕量較凹形坡地塊的侵蝕量增加了一倍，其中淺溝侵蝕量為總侵蝕量的86.7%。Deoliverira[90]研究了坡度和坡形對(duì)切溝的影響，結(jié)果表明切溝常常在坡面中段坡度最陡處出現(xiàn)，有的甚至在中間出現(xiàn)，而下部卻是淺溝侵蝕。這是坡度大導(dǎo)致徑流下切能力突然增大的原因。

隨著GIS和遙感成像技術(shù)的出現(xiàn)，人們更加關(guān)注地表過(guò)程。秦偉等[91]基于Quickbird高分辨率遙感影像和數(shù)字高程模型，提取了坡面淺溝分布的地形參數(shù)，發(fā)現(xiàn)在黃土丘陵溝壑區(qū)，坡面坡度、長(zhǎng)度、坡向以及匯水坡長(zhǎng)是影響坡面淺溝數(shù)量的主要地形要素，發(fā)生淺溝侵蝕臨界坡度的上限與下限分別介于26°～27°和15°～20°，而臨界坡長(zhǎng)介于50～80 m。

植被覆蓋通過(guò)削弱雨滴擊濺，增加入滲率和增強(qiáng)土壤抗沖性等方面的作用而有效地減少淺溝侵蝕。野外調(diào)查發(fā)現(xiàn)，已經(jīng)發(fā)育成淺溝的坡面一經(jīng)撂荒后，因淺溝底部生草的根系從而抑制淺溝侵蝕的發(fā)生[75]。而且，隨著植被演替的恢復(fù)，土壤侵蝕量呈現(xiàn)明顯的下降趨勢(shì)。土壤抗蝕性綜合指數(shù)強(qiáng)弱為灌木群落階段>多年生草本和蒿類(lèi)群落階段>一二年生草本群落階段。與一二年生草本群落階段相比，灌木群落階段與多年生草本和蒿類(lèi)群落階段的土壤抗蝕性分別增加了3.6～6.7 倍和5.7～9.3倍[75]。就優(yōu)勢(shì)植物群落而言，狼牙刺群落、刺槐和沙棘群落抗侵蝕性能最好，豬毛蒿群落抗侵蝕性能最低[92]。于國(guó)強(qiáng)等[59]研究發(fā)現(xiàn)，有植被覆蓋時(shí)，根系加固作用改善了坡面淺層的土體應(yīng)力，有效地降低了坡面土體的集中應(yīng)力；減小了坡溝系統(tǒng)坡面淺層土體位移，使水平位移減少15%；鉛垂位移減少2.5%；并減少了塑性屈服區(qū)體積，使剪切塑性區(qū)減少11.5%，張拉塑性區(qū)減少84.0%；減輕了重力侵蝕的破壞程度，但并未改變坡溝系統(tǒng)以剪切破壞力為主的屈服模式。該研究為更好理解植被根系調(diào)控重力侵蝕研究提供了新的技術(shù)方法。

溝道侵蝕與其它形式土壤退化的區(qū)別是地下深層過(guò)程不同。Rubey[7]是首位研究溝道侵蝕的地下過(guò)程重要性的科學(xué)家。地下過(guò)程因素主要涉及管路、隧道、巖石學(xué)因子和土壤剖面過(guò)程。這些過(guò)程的研究主要涉及土壤的可分散性或黏粒含量、活躍的地下水流和坡度梯度等。但是，與表層動(dòng)態(tài)比較，由于亞表層現(xiàn)象中存在更大的差異性，因此評(píng)估起來(lái)難度更大，導(dǎo)致有關(guān)土壤開(kāi)裂、地下水和滲漏方面的研究工作十分有限。

經(jīng)驗(yàn)分析表明，永久溝道與臨時(shí)溝道相比：深層形態(tài)趨于產(chǎn)生更高的土壤侵蝕模數(shù)。這似乎與已有的結(jié)論相違：即當(dāng)臨時(shí)溝道保持不變時(shí)，溝道演化呈降低趨勢(shì)[93]。今后的研究應(yīng)該側(cè)重什么條件(地表限制、土地管理、植被重建的程度)會(huì)導(dǎo)致不同的結(jié)果，這有利于解釋現(xiàn)有的觀(guān)測(cè)結(jié)果，也有利于制定及實(shí)施防控和恢復(fù)策略。

4 溝道侵蝕發(fā)生的臨界條件

有關(guān)溝道形成的地形閾值的評(píng)估是溝道侵蝕最熱門(mén)的研究課題。Schumm和 Hadley[94]認(rèn)識(shí)到流域的地形特征對(duì)溝道發(fā)育的重要性，并認(rèn)為半干旱山谷區(qū)土壤侵蝕的發(fā)生有一個(gè)臨界坡度閾值。Schumm[95]進(jìn)一步強(qiáng)調(diào)了地形閾值在侵蝕沉積土地演化中的意義。Patton 和Schumm[96]有關(guān)美國(guó)科羅拉多州西北部溝道的研究發(fā)現(xiàn)，溝道上方的局地坡面臨界坡度(S)和上坡匯水面積(A)之間一般存在反向關(guān)系，數(shù)據(jù)的底線(xiàn)值能夠建立坡面與流域面積的臨界關(guān)系(圖1)，即對(duì)一定排水面積的區(qū)域而言，高于直線(xiàn)(臨界值)以上的區(qū)域就會(huì)產(chǎn)生溝道。此關(guān)系可用于辨識(shí)潛在的不穩(wěn)定的溝谷，有助于土地經(jīng)營(yíng)者確定不穩(wěn)定的區(qū)域，制定有效的經(jīng)濟(jì)措施。他們同時(shí)強(qiáng)調(diào)，此特定的數(shù)量關(guān)系僅適用于他們研究的區(qū)域，即科羅拉多州西北部的半干旱Piceance彎和黃彎流域。對(duì)于具有其他特性的區(qū)域，還要考慮額外的因素給予分析，但流域穩(wěn)定性的基本理論具有應(yīng)用價(jià)值。由此開(kāi)創(chuàng)了溝道侵蝕地形臨界條件研究的新紀(jì)元。此項(xiàng)研究工作的開(kāi)展，源于當(dāng)時(shí)有關(guān)溝道的臨界坡度可能與徑流程度相關(guān)的認(rèn)知，但缺少測(cè)定排水的方法。更受到Burkham[97]在美國(guó)新墨西哥州半干旱區(qū)域有關(guān)排水和流域面積高度相關(guān)的研究結(jié)果的啟發(fā)(圖1)。

Begin和Schumm[98]及Vandaele等[99]提出水力侵蝕導(dǎo)致溝道發(fā)育的臨界關(guān)系可以用冪函數(shù)S=aAb表示，即由溝頭排水面積和溝頭坡度的對(duì)數(shù)坐標(biāo)軸說(shuō)明，或者說(shuō)溝道侵蝕是坡度和匯水面積的函數(shù)[23]。其中，S單位為m·m-1；A單位為hm2；a值代表著溝蝕發(fā)生所需的臨界值；b為一無(wú)量綱量。b值在理論上則代表著匯水面積的相對(duì)重要性，它受到氣候、土壤、主導(dǎo)徑流過(guò)程、降雨特征和植被等土地利用狀況的影響。

Desmet等[100]認(rèn)為溝道侵蝕是坡面集中徑流侵蝕的結(jié)果，其形成和發(fā)展取決于一定的徑流量和徑流動(dòng)能，而上方匯水面積及坡度決定了徑流量和徑流動(dòng)能的大小。Poesen等[27]認(rèn)為，只有在徑流產(chǎn)生的剪切力大于土壤的抗蝕力時(shí)才會(huì)有溝道的出現(xiàn)。Sidorchuk[24]也認(rèn)為當(dāng)徑流特性，如剪切應(yīng)力超過(guò)一個(gè)關(guān)鍵值時(shí)，才有可能發(fā)生溝道。而地表徑流強(qiáng)度主要是由匯水面積和局地坡度等地貌參數(shù)決定的[101]。

在坡度一定的前提下，一般從濕潤(rùn)氣候到干旱氣候，產(chǎn)生同樣徑流能量所需的匯水面積逐漸變大；地下徑流過(guò)程的加入使得產(chǎn)生等同徑流能量的匯水面積變小，如果是以塊體運(yùn)動(dòng)和崩塌等地貌為主，b值甚至可以變?yōu)樨?fù)值，即意味著S-A 關(guān)系呈正相關(guān)關(guān)系。

低頻高強(qiáng)度降水使較小的匯水面積就能產(chǎn)生同樣的徑流能量；良好的植被覆蓋度一方面可以緩沖降水對(duì)地表的沖擊，另一方面還可以增加土壤抗蝕性及降水入滲，從而減少地表徑流，因此，產(chǎn)生同樣的徑流能量就需要更大的匯水面積。這就是為什么淺溝侵蝕常發(fā)生于坡度較陡的坡面上，而且一般多發(fā)生于具有一定匯水面積的坡面中部和中下部的主要原因。

圖1 美國(guó)科羅拉多州Piceance流域溝道和非溝道區(qū)坡度與排水面積關(guān)系[96]Fig.1 Slope-area relation for gullied and ungullied valley in Piceance basin of Northwestern Colorado[96]

顯然，在坡度一定的情況下，切溝下切形成時(shí)所需的臨界面積要比淺溝大。淺溝的發(fā)育規(guī)模形態(tài)要比切溝小，而侵蝕體積大小與剪切力有關(guān)[102]。另一方面，剪切力主要由水流量和坡度決定，在以表面漫流為主的流域中，流量可以用匯水面積來(lái)代替，也就是說(shuō)剪切力主要由匯水面積和坡度決定。因此，在假定其他因素相同的情況下，小侵蝕形態(tài)只需較小的匯水面積就可以形成，而發(fā)育成切溝則需要相對(duì)較大的匯水面積。

一般而言，對(duì)于不同的地區(qū)，都會(huì)有不同的溝道起始臨界值a和b。b值正值是與表面漫流產(chǎn)生的侵蝕有關(guān)，而負(fù)值則與滲蝕(Seepage erosion)和發(fā)育崩塌的地下過(guò)程等有關(guān)[103]。

臨界模型中的參數(shù)一般通過(guò)實(shí)測(cè)溝道上方的匯水面積(A)和局地坡度(S)數(shù)據(jù)進(jìn)行推求，方法有兩種：一是分散數(shù)據(jù)的下限所對(duì)應(yīng)的參數(shù)；二是通過(guò)數(shù)據(jù)的回歸得出的參數(shù)。這兩套參數(shù)用于不同的研究目的：第一種主要進(jìn)行溝蝕的預(yù)測(cè)，第二種主要用于對(duì)主導(dǎo)溝蝕過(guò)程的分析。只要確定了臨界模型中的系數(shù)a和b值，就可以利用其預(yù)測(cè)將發(fā)生溝蝕的區(qū)域。

Begin和Schumm[100]通過(guò)理論分析得到的b值在0.2～0.4之間。Morgan 和Mngomezulu[103]認(rèn)為b值小于0.2時(shí)將會(huì)有地下過(guò)程參與。

不同學(xué)者基于各自研究區(qū)域的地形及侵蝕特點(diǎn)，建立了不同地區(qū)的淺溝侵蝕臨界模型。

胡剛等[104]分析東北漫川漫崗黑土區(qū)淺溝和切溝侵蝕與地貌參數(shù)發(fā)現(xiàn)：(1)無(wú)論切溝還是淺溝的上坡匯水面積和局地坡度都呈較高的負(fù)冪函數(shù)相關(guān)關(guān)系；(2)坡度一定時(shí),切溝開(kāi)始下切所需的臨界面積要比淺溝大；同樣,匯水條件一定時(shí),發(fā)育切溝需要的坡度比淺溝更大；(3)切溝發(fā)生的地貌臨界公式為：S =0.1161 A-0.4457，淺溝公式為：S =0.0631 A-0.4643。此外，他們還校驗(yàn)了Moore等(1988)[105]的溝蝕發(fā)生公式，通過(guò)與野外實(shí)測(cè)淺溝和切溝發(fā)生位置對(duì)比，臨界模型預(yù)測(cè)的溝蝕位置較好地反映了野外實(shí)際狀況。由此提出，地下徑流過(guò)程的參與可能是最終導(dǎo)致切溝的S-A 回歸趨勢(shì)線(xiàn)斜率即b值較淺溝小的原因，只有在SAb>a時(shí)，才會(huì)有溝道侵蝕的發(fā)生。

對(duì)于淺溝和切溝局地坡度和匯水面積的測(cè)量及量算有不同標(biāo)準(zhǔn)，一旦切溝或淺溝形成就會(huì)發(fā)生溯源侵蝕，溯源侵蝕會(huì)影響匯水面積，改變坡面坡度。一般而言,在臨界坡度下坡面侵蝕量最大，片蝕最容易過(guò)渡到細(xì)溝侵蝕，進(jìn)而發(fā)展到溝蝕[106]。

李斌兵等[107]通過(guò)GPS 實(shí)測(cè)數(shù)據(jù)并結(jié)合GIS空間分析與統(tǒng)計(jì)回歸方法，建立了適用于黃土高原丘陵區(qū)的切溝侵蝕臨界模型(SA0.1351> 1.948)，發(fā)現(xiàn)在黃土高原丘陵溝壑區(qū)，隨著坡度的增大，發(fā)生切溝侵蝕的臨界值a 增大，高強(qiáng)度降雨致使判定式中匯水面積的指數(shù)b值減小，從而減少了匯水面積的影響作用。同時(shí)臨界模型的驗(yàn)證結(jié)果表明，基于所建臨界模型提取的淺溝侵蝕分布區(qū)與野外實(shí)際相吻合。

Torri和Poesen[108]對(duì)比了世界各地63 個(gè)已經(jīng)建立的形成淺溝的臨界地形模型，并對(duì)不同環(huán)境狀況下的參數(shù)值變化規(guī)律進(jìn)行了總結(jié)，認(rèn)為需要建立更完善的，包括土地利用、氣候變化和自然災(zāi)害因素在內(nèi)的基于物理過(guò)程的模型來(lái)判定淺溝發(fā)生的臨界條件。因?yàn)闇系赖漠a(chǎn)生是由作用在溝頭上的動(dòng)力過(guò)程所控制，這些過(guò)程包括地表漫流、地下水導(dǎo)致的滲流和潛蝕以及塊體塌陷或崩塌等因素。

有關(guān)淺溝形成、臨界地形條件及其影響因素的研究，鄭粉莉等[46]做了非常詳細(xì)的評(píng)述，感興趣的讀者可以閱讀原文。

5 溝道侵蝕研究展望

作為土壤侵蝕科學(xué)中一個(gè)獨(dú)立的研究領(lǐng)域，溝道侵蝕研究發(fā)表的論文只占土壤侵蝕方面的10%。分析溝道侵蝕的研究?jī)?nèi)容，大約有2/3以上的研究集中在溝道的觀(guān)測(cè)、退化速率的定量化、退化過(guò)程的評(píng)估及其影響因素的分析等，長(zhǎng)期和短期評(píng)價(jià)主要基于對(duì)正射投影的解釋?zhuān)渌悄Ｐ秃蜏系李A(yù)測(cè)指數(shù)的研究，很少一部分是實(shí)驗(yàn)室研究、定性描述和防治研究。溝道預(yù)測(cè)指數(shù)的研究主要基于環(huán)境指標(biāo)的分析(如地形等因素)來(lái)預(yù)測(cè)溝道的地點(diǎn)和嚴(yán)重程度[16]。

上世紀(jì)90年代初，學(xué)者們主要開(kāi)展了有關(guān)溝道測(cè)定和定性方面的研究，之后實(shí)驗(yàn)室模擬研究受到重視。1999-2005年間，開(kāi)始開(kāi)發(fā)和利用侵蝕模型開(kāi)展研究。2006年以來(lái)，溝道侵蝕預(yù)測(cè)指數(shù)的研究受到重視。同期，溝道侵蝕對(duì)土壤質(zhì)量、生產(chǎn)力和生產(chǎn)成本的研究也逐步受到關(guān)注。目前的研究趨向于更復(fù)雜的表面過(guò)程和地下過(guò)程分析。這些研究加深了人們對(duì)溝道侵蝕發(fā)生發(fā)展過(guò)程及其防治的認(rèn)識(shí)，也為溝道侵蝕危害性評(píng)價(jià)提供了一定的依據(jù)。深入開(kāi)展對(duì)溝道的監(jiān)測(cè)、研究和模擬，是預(yù)測(cè)環(huán)境變化對(duì)溝道侵蝕速率影響的基礎(chǔ)，有助于人類(lèi)深入認(rèn)識(shí)溝道侵蝕發(fā)生的機(jī)理，從而采取阻控措施以減緩溝道侵蝕的持續(xù)發(fā)生。

但是，從時(shí)間和強(qiáng)度角度看，溝道侵蝕并不是一個(gè)均質(zhì)同類(lèi)的過(guò)程，具有多要素和非一致性的形成特點(diǎn)，依賴(lài)于尺度的行為，且年內(nèi)存在著差異性(氣候促使和深層剖面的依賴(lài)性)[16]。而且，長(zhǎng)期研究表明全球范圍內(nèi)溝道侵蝕的驅(qū)動(dòng)力并不相同。例如，二十世紀(jì)中葉機(jī)械的引進(jìn)是歐洲近期溝道侵蝕的主要因素[109-110]。由于溝道侵蝕過(guò)程的復(fù)雜性、非線(xiàn)性、突發(fā)性及區(qū)域差異性等特點(diǎn)，加之研究方法的局限和受人為耕作活動(dòng)的特殊性影響，溝道侵蝕的研究尚有許多亟需深入探討的問(wèn)題。

鄭粉莉等[47]提出，目前溝道侵蝕過(guò)程的定量研究、溝道發(fā)生的臨界地形和動(dòng)力條件、發(fā)育過(guò)程的定量描述、水流的剝離方程及泥沙搬運(yùn)能力，包括切溝侵蝕的流域侵蝕預(yù)報(bào)模型、動(dòng)態(tài)監(jiān)測(cè)切溝侵蝕過(guò)程的新方法等研究還相對(duì)薄弱，有待加強(qiáng)。且至今尚未有公認(rèn)的溝蝕過(guò)程定量表達(dá)式，致使現(xiàn)有侵蝕預(yù)報(bào)模型尚不能預(yù)報(bào)溝蝕量，并對(duì)控制溝道侵蝕動(dòng)態(tài)的基本過(guò)程分析關(guān)注的較少，迫切需要加強(qiáng)溝蝕過(guò)程的研究。

土壤侵蝕發(fā)生規(guī)律具有很強(qiáng)的區(qū)域特征，與土壤性質(zhì)和氣候緊密相關(guān)。我國(guó)黑土區(qū)具有典型的漫川漫崗地形和相對(duì)集中的降雨特征，土壤黏重，雖然研究者對(duì)黑土區(qū)土壤侵蝕機(jī)理的研究已經(jīng)開(kāi)展了一些工作，但是相對(duì)其他區(qū)域的研究還比較少。因此，今后要深入研究東北黑土區(qū)溝道侵蝕空間分布特征，定量化標(biāo)識(shí)東北黑土區(qū)細(xì)溝及細(xì)溝間侵蝕，界定侵蝕溝生成的臨界條件，明晰溝蝕發(fā)生發(fā)展過(guò)程及其多種侵蝕外營(yíng)力交互、耦合與疊加作用機(jī)制，揭示溝道侵蝕機(jī)理，提出東北黑土區(qū)溝道侵蝕防治、修復(fù)措施，為選擇與布設(shè)提供科學(xué)依據(jù)。

根據(jù)國(guó)內(nèi)外溝道侵蝕研究的進(jìn)展和成果，以下方面仍需要進(jìn)一步加強(qiáng)和明確：(1)從水力學(xué)、降雨、地形、土壤和土地利用角度，明確特定區(qū)域溝道發(fā)生、發(fā)展和填充的水力閾值、降雨閾值和地形閾值；(2)從土壤抗切割的能力、土壤剖面中滲透條件或坡度的不穩(wěn)定性方面，定量評(píng)估地下過(guò)程對(duì)溝道侵蝕的動(dòng)態(tài)學(xué)和敏感性影響；(3)凍融過(guò)程、土體開(kāi)裂、地下水水文學(xué)和滲漏力對(duì)溝道形成的作用，溝道侵蝕與水文和其他土壤退化過(guò)程的互作作用；(4)既然細(xì)溝和片蝕有一個(gè)可接受的土壤損失值，那么溝道侵蝕是否需要確定一個(gè)最小的或者是一個(gè)可以接受的土壤損失值？以此確定什么條件下，需要采取控制措施；(5)管理措施的空間模式是如何影響溝道侵蝕空間變化和沉積速率的？不同土地利用、不同氣候環(huán)境下，不同時(shí)空尺度溝道侵蝕對(duì)整個(gè)土壤流失和沉積的貢獻(xiàn)率如何？(6)對(duì)易發(fā)生溝道侵蝕的區(qū)域，在溝道侵蝕的初期，土地管理、植被重建對(duì)溝道侵蝕發(fā)生發(fā)展的影響程度，什么樣的農(nóng)業(yè)或工程措施可以盡早有效地控制溝頭的挺進(jìn)；(7)有效阻止和控制溝道的技術(shù)方法以及溝道治理計(jì)劃成敗的教訓(xùn)等，由此制定防控和恢復(fù)策略；(8)研發(fā)預(yù)測(cè)各種時(shí)空條件下溝道侵蝕速率、溝道對(duì)水文學(xué)、產(chǎn)沙量和景觀(guān)演化的適宜模型；(9)溝道普查的方法也需要改進(jìn)，除了對(duì)溝道形態(tài)的時(shí)間和空間測(cè)定，也要對(duì)排水量和產(chǎn)沙量等進(jìn)行測(cè)定。陸地?cái)z影測(cè)量方法和三維攝影重建技術(shù)已經(jīng)用于測(cè)定和模擬細(xì)溝和溝道，今后需要探討不同時(shí)空尺度下，監(jiān)測(cè)和研究適宜溝道類(lèi)型的測(cè)定技術(shù)；(10)全球變化對(duì)溝道侵蝕的直接(降雨量、降雨強(qiáng)度和降雨時(shí)空分布變化)和間接影響(全球增溫影響植被覆蓋、土壤含水量變化以及人類(lèi)社會(huì)經(jīng)濟(jì)因素的變化)也需要關(guān)注。

致 謝

感謝“國(guó)家科技基礎(chǔ)條件平臺(tái)-國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)共享服務(wù)平臺(tái)-東北黑土科學(xué)數(shù)據(jù)中心(http://northeast.geodata.cn)”提供數(shù)據(jù)支撐。

參考文獻(xiàn)(References)：

[1] MADISON W I.Soil science society of America,glossary of soil science terms[R].Soil Science Society of America,2001.

[2] BILLI P,DRAMIS F.Geomorphological investigation on gully erosion in the Rift Valley and the northern highlands of Ethiopia[J].Catena,2003,50(2-4): 353-368.

[3] KRAUSE A K,FRANKS S W,KALMA J D,et al.Multi parameter finger printing of sediment deposition in a small gullied catchment in SE Australia[J].Catena,2003,53(4): 327-348.

[4] DE VENTE J,POESEN J,VERSTRAETEN G.The application of semi-quantitative methods and reservoir sedimentation rates for the prediction of basin sediment yield in Spain[J].Journal of Hydrology,2005,305(1-4): 63-86.

[5] PIMENTEL D,BURGESS M.Soil erosion threatens food production[J].Agriculture,2013,3(3): 443-463.

[6] IRELAND H A.“Lyell”gully,a record of a century of erosion[J].Journal of Geology,1939,47(1): 47-63.

[7] RUBEY W W.Gullies in the Great Plains formed by sinking of the ground[J].American Journal of Science,1928,15(89): 417-422.

[8] DOTTERWEICH M,RODZIK J,ZGLOBICKI W,et al.High resolution gully erosion and sedimentation processes,and land use changes since the Bronze Age and future trajectories in the Kazimierz Dolny area(Naleczow Plateau,SE-Poland)[J].Catena,2012,95: 50-62.

[9] WEBB R H,HEREFORD R.Floods and geomorphic change in the southwestern United States:An historical perspective[C].Proceedings of Seventh Federal Interagency Sedimentation Conference,Mar.25-29,2001.Reno,NV,USA:Federal Interagency Sedimentation Committee,2001:IV30-IV37.

[10] HARVEY A M.Holocene hillslope gully systems in the How gill Fells,Cumbria[M]//In:ANDERSON M G，eds.Advances in Hillslope Processes.New York:John Wiley and Sons Ltd,1996,2: 731-752.

[11] BORK H R,LI Y,ZHAO Y T,et al.Land use changes and gully development in the upper Yangtze River Basin,SW-China[J].Journal of Mountain Science,2001,19(2): 97-103.

[12] ZHANG X Y,CRUSE R M,SUI Y Y,et al.Soil compaction induced by small tractor traffic in Northeast China[J].Soil Science Society of America Journal,2006,70(2): 613-619.

[13] 劉興土,閻百興.東北黑土區(qū)水土流失與糧食安全[J].中國(guó)水土保持,2009(1): 17-19.

LIU X T,YAN B X.Soil erosion anf grain safety in the black soil area of Northeast in China[J].Soil and Water Conservation in China,2009(1): 17-19.

[14] LIU X B,ZHANG X Y,WANG Y X,et al.Soil degradation: A problem threatening the sustainable development of agriculture in Northeast China[J].Plant Soil and Environment,2010,56(2): 87-97.

[15] CHAPLOT V,BROWN J,DLAMINI P,et al.Rainfall simulation to identify the storm-scale mechanisms of gully bank retreat[J].Agricultural Water Management,2011,98(11):1704-1710.

[17] EUSTACE A H,PRINGLE M J,DENHAM R J.A risk map for gully locations in central Queensland,Australia[J].European Journal of Soil Science,2011,62:(3) 431-441.

[18] MORGAN R P C.Soil erosion and conservation[M].Cornwall: Blackwell Publishing,2005.

[19] Soil Science Society of America,Glossary of Soil Science Terms[M/OL].[S.l.]:[s.n.]，2016.[2018-06-01].https://www.soils.org/publications/soils-glossary#.

[20] BETTIS E A,THOMPSON D M.Gully erosion[J].Rangelands,1985,7(2): 70-72.

[21] DAVIS W M.Channels,valleys and intermont detrital plains[J].Science,1927,66(1708): 272-274.

[22] KNIGHTON D.Fluvial forms and processes a new perspective[M].London: Routledge,1998.

[23] SIDORCHUK A.Stages in gully evolution and self-organized criticality[J].Earth Surface Processes and Landforms,2006,31(11): 1329-1344.

[24] DUNNE T.Field studies of hillslope flow processes[M]//In:KIRKBY M J,eds.Hill slope hydrology.Chichester: Wiley,1978,227-93.

[25] FOSTER G R.Understanding ephemeral gully erosion[J].Soil Conservation,1986,2:90-125.

[26] POESEN J,NACHTERGAELE J,VERSTRAETEN G,et al.Gully erosion and environmental change:Importance and research needs[J].Catena,2003,50(2-4): 91-133.

[27] RUSTOMJI P.Analysis of gully dimensions and sediment texture from southeast Australia for catchment sediment budgeting[J].Catena,2006,67(2): 119-127.

[28] HUDSON N W.Soil conservation[M].London: Batsford,1985.

[29] WU Y Q,CHENG H.Monitoring of gully erosion on the Loess Plateau of China using a global positioning system[J].Catena,2005,63(2-3): 154-166.

[30] WSTSON D A,LAFLEN J M,FRANTI T G.Estimating ephemeral gully erosion[J].Journal of Municipal Problems,1986,61: 86-2020.

[32] CAPRA A,SCICOLONE B.Ephemeral gully erosion in a wheat-cultivated area in Sicily(Italy)[J].Biosystems Engineering,2002,83(1): 119-126.

[33] CAPRA A.Ephemeral gully and gully erosion in cultivated land:A review[M]//In:EVELYN C L,eds.Drainage basins and catchment management:Classification,modelling,and environmental assessment.New York: Nova Science Publishers,Inc,2013: 109-141.

[34] AUZET DR A V,BOIFFIN DR J,PAPY DR F,et al.Rill erosion as a function of the characteristics of cultivated catchments in the north of France[J].Catena,1993,20(1-2): 41-62.

[35] WILSON G V,CULLUM R F,R?MKENS M J M.Ephemeral gully erosion by preferential flow through a discontinuous soil-pipe[J].Catena,2008,73(1): 98-106.

[36] TAGUAS E V,YUAN Y,BINGNER R L,et al.Modelling the contribution of ephemeral gully erosion under different soil managements: A case study in an olive orchard microcatchment using the AnnAGNPS model[J].Catena,2012,98: 1-16.

[37] CHAPLOT V,BROZEC E C L,SILVERA N,et al.Spatial and temporal assessment of linear erosion in catchments under sloping lands of northern Laos[J].Catena,2005,63(2-3): 167-184.

[38] LIU H H,ZHANG T Y,LIU B Y,et al.Effects of gully erosion and gully filling on soil depth and crop production in the black soil region,northeast China[J].Environment Earth Science,2013,68(6): 1723-1732.

[39] WOODWARD D E.Method to predict cropl and ephemeral gully erosion[J].Catena,1999,37(3-4): 393-399.

[40] TANG J L,CHENG X Q,ZHU B,et al.Rainfall and tillage impacts on soil erosion of sloping cropland with subtropical monsoon climate:A case study in hilly purple soil area,China[J].Journal of Mountain Science,2015,12(1): 134-144.

[41] VALENTIN C,POESEN J,YONG L.Gully erosion:Impacts,factors and control[J].Catena,2005,63(2-3): 132-153.

[42] VANWALLAGHEM T,BORK H R,POESEN J,et al.Rapid development and infilling of a buried gully under cropland,central Belgium[J].Catena,2005,63(2-3):221-243.

[43] 朱顯謨.黃土區(qū)土壤侵蝕的分類(lèi)[J].土壤學(xué)報(bào),1956,4(2): 99-115.

ZHU X M.Classification on the soil erosion in the loess region[J].Acta Pedologica Sinica,1956,4(2): 99-115.

[44] 劉元保,朱顯謨,周佩華,等.黃土高原坡面溝蝕的類(lèi)型及其發(fā)生發(fā)展規(guī)律[J].水土保持研究,1988,1: 9-18.

LIU Y B,ZHU X M,ZHOU P H,et al.The laws of hillslope channel erosion occurrence and development on Loess Plateau[J].Research of Soil and Water Conservation,1988,1: 9-18.

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

ZHENG F L,WU M,ZHANG Y B,et al.Ephemeral gully development process at Loess Steep Hillslope[J].Scientia Geographica Sinica,2006,26(4): 438-442.

[46] 鄭粉莉,徐錫蒙,覃 超.溝蝕過(guò)程研究進(jìn)展[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(8): 48-59,116.

ZHENG F L,XU X M,QIN C.A review of gully erosion on process research[J].Transactions of the Chinese Society for Agricultural Machinery,2016,47(8): 48-59,116.

[47] 中國(guó)農(nóng)業(yè)百科全書(shū)土壤卷編委會(huì)．土壤侵蝕與水土保持分支條目[M]．北京：農(nóng)業(yè)出版社,1996．

Soil Volume Editorial Board of China Agricultural Encyclopedia.Entry of soil erosion and water conservation[M]．Beijing:China Agriculture Press,1996．

[48] ZHENG F L,HUANG C H.Gully erosion[M]∥In:RATTAN L,eds.Encyclopedia of soil science.New York: Marcel Dekker,Inc,2002: 630-634.

[50] POESEN J.Challenges in gully erosion research[J].Landform Analysis,2011,17: 5-9.

[51] ST?CKER C,ELTNER A,KARRASCH P.Measuring gullies by synergetic application of UAV and close range photogrammetry-A case study from Andalusia,Spain[J].Catena,2015,132: 1-11.

[52] HORTON R E.Erosional development of streams and their drainage basins; hydrophysical approach to quantitative morphology[J].Bulletin of the American Geological Society,1945,56(3): 275-370.

[53] PIEST R F,SPOMER R G.Sheet and gully erosion in the Missouri valley loessial region[J].Transactions American Society of Agricultural Engineers,1968,11: 850-853.

[54] BRYAN R B,YAIR A.Perspectives on studies of badlands geomorphology[M]//In:BRYAN R B and YAIR A,eds.Badlands geomorphologyand piping.Norwich: GeoBooks,1982: 1-12.

[55] 景 可.黃土高原溝谷侵蝕研究[J].地理科學(xué),1986,6(4): 340-347.

JING K.A study on gully erosion on the Loess Plateau[J].Scientia Geographica Sinica,1986,6(4): 340-347.

[56] BOCCO G.Gully erosion—processes and models[J].Progress in Physical Geography,1991,15(4): 392-406.

[57] ROBERSON K M,HANSON G J.Gully headcut advance[J].Transactions of ASAE,1996,39(1): 33-38.

[58] SOUFI M.澳大利亞一個(gè)林地環(huán)境下的切溝發(fā)展過(guò)程和侵蝕趨勢(shì)研究[J].中國(guó)水土保持,2002,23(7): 26-27.

SOUFI M.Research on gully development process and erosion trend in a woodland environment in Australia[J].Soil and Water Conservation in China,2002,23(7): 26-27.

[59] 于國(guó)強(qiáng),張 霞,張茂省,等.植被對(duì)黃土高原坡溝系統(tǒng)重力侵蝕調(diào)控機(jī)理研究[J].自然資源學(xué)報(bào),2012,27(6): 922-932.

YU G Q,ZHANG X,ZHANG M S,et al.Mechanism of vegetation regulating on gravitational erosion in the slope-gully system on Loess Plateau[J].Journal of Natural Resources,2012,27(6): 922-932.

[60] PARKNER T,PAGE M J,MARUTANI T,et al.Development and controlling factors of gullies and gully complexes,East Coast,New Zealand[J].Earth Surface Processes and Landforms,2006,31(2),187-199.

[61] NYSSEN J,POESEN J,VEYRET-PICOT M,et al.Assessment of gully erosion rates through interviews and measurements:A case study from northern Ethiopia[J].Earth Surface Processes and Landforms,2006,31(2): 167-185.

[62] HADLEY R F,LAL R,ONSTAD C A,et al.Recent developments in erosion and sediment yield studies[M].Technical Documents in Hydrology,UNESCO,1985.

[63] 張科利.陜北黃土丘陵溝壑區(qū)坡耕地淺溝及其防治途徑[D].楊陵: 中國(guó)科學(xué)院水利部水土保持研究所,1988.

ZHANG K L.Shallow furrow of sloping farmland and its control in Loess Hilly and gully region of Northern Shaanxi Province[D].Yangling: Institute of Soil and Water Conservation,CAS and MWR,1988.

[64] CAPRA A,PORTO P,SCICOLONE B.Relationships between rainfall characteristics and ephemeral gully erosion in a cultivated catchment in Sicily(Italy)[J].Soil and Tillage Research,2009,105(1): 77-87.

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

ZHENG F L,KANG S Z.Erosion and sediment yield in different zones of loess slopes[J].Acta Geographica Sinica,1998,53(5): 422-428.

[66] 鄭粉莉,高學(xué)田.黃土坡面土壤侵蝕過(guò)程與模擬[M].西安:陜西人民出版社,2000: 96-119.

ZHENG F L,GAO X T.Soil ersion progresses and modeling at loessial hillslope[M].Xi′an:Shaanxi People′s Publishing House,2000:96-119.

[67] WILSON G.Understanding soil pipe flow and its role in ephemeral gully erosion[J].Hydrological Processes,2011,25(15): 2354-2364.

[68] AVNI Y.Gully incision as a key factor in desertification in an environment,the Negev highlands,Israel[J].Catena,2005,63(2-3): 185-220.

[69] ZAIMES G N,SCHULTZ R C.Assessing riparian conservation land management practice impacts on gully erosion in Iowa[J].Environmental Management,2012,49(5): 1009-1021.

[70] MAUGNARD A,CORDONNIER H,DEGRE A,et al.Uncertainty assessment of ephemeral gully identification,characteristics and topographic threshold when using aerial photographs in agricultural settings[J].Earth Surface Processes and Landforms,2014,39(10): 1319-1330.

[71] COLLISON A J C.The cycle of instability:Stress release and fissure flow as controls on gully head retreat[J].Hydrological Processes,2001,15(1): 3-12.

[72] LESSCHEN J P,KOK K,VERBURG P H,et al.Identification of vulnerable areas for gully erosion under different scenarios of land abandonment in Southeast Spain[J].Catena,2007,71(1): 110-121.

[73] FAUST D,SCHMIDT M.Soil erosion processes and sediment fluxes in a Mediterranean landscape of marls(Campina de Cadiz,SW Spain)[J].Zeitschrift Fur Geomorphologie,2009,53(2):247-265.

[74] SAMANI A N,AHMADI H,MOHAMMADI A,et al.Factors controlling gully advancement and models evaluation(Hableh Rood Basin,Iran)[J].Water Resources Management,2010,24(8): 1531-1549.

[75] 韓魯艷,賈燕鋒,王 寧,等.黃土丘陵溝壑區(qū)植被恢復(fù)過(guò)程中的土壤抗蝕與細(xì)溝侵蝕演變[J].土壤,2009,41(3): 483-489.

HAN L Y,JIA Y F,WANG N,et al.Soil anti-erodibility and soil erosion evolution during process of vegetation recovering in loess hilly-gully region[J].Soils,2009,41(3): 483-489.

[76] BETTS H D,TRUSTRUM N A,DE ROSE R C.Geomorphic changes in a complex gully system measured from sequential digital elevation models,and implications for management[J].Earth Surface Processes and Landforms,2003,28(10): 1043-1058.

[77] GRELLIER S,KEMP J,JANEAU J L,et al.The indirect impact of encroaching trees on gully extension: A 64 year study in a sub-humid grassland of South Africa[J].Catena,2012,98: 110-119.

[78] SEUTLOALI K E,BECKEDAHL H R,DUBE T,et al.An assessment of gully erosion along major armored roads in south-eastern region of South Africa:A remote sensing and GIS approach[J].Geocarto International,2016,31(2): 225-239.

[79] OKOYEH E I,AKPAN A E,EGBOKA B C E,et al.An assessment of the influences of surface and subsurface water level dynamics in the development of gullies in Anambra State,Southeastern Nigeria[J].Earth Interactions,2014,18(4): 1-24.

[80] HAILE G W,FETENE M.Assessment of soil erosion hazard in Kilie catchment,East Shoa,Ethiopia[J].Land Degradation and Development,2012,23(3): 293-306.

[81] TEBEBU T Y,ABIY A Z,DAHLKE H E,et al.Surface and subsurface flow effect permanent gully formation and upland erosion near Lake Tana in the northern highlands of Ethiopia[J].Hydrology and Earth System Sciences Discussions,2010,14(11): 2207-2217.

[82] VOARINTSOA N R G,COX R,RAZANATSEHENO M O M,et al.Relation between bedrock geology,topography and lavaka distribution in Madagascar[J].South African Journal of Geology,2012,115(2): 225-250.

[83] BEAVIS S G.Structural controls on the orientation of erosion gullies inmid-western New South Wales,Australia[J].Geomorphology,2000,33(1-2): 59-72.

[84] PARKNER T,PAGE M,MARDEN M,et al.Gully systems under undisturbed indigenous forest,East Coast Region,New Zealand[J].Geomorphology,2007,84(3-4): 241-253.

[85] BOARDMAN J,PARSONS A J,HOLLAND R,et al.Development of badlands and gullies in the Sneeuberg,Great Karoo,South Africa[J].Catena,2003,50(2-4): 165-184.

[86] FRANKL A,POESEN J,HAILE M,et al.Quantifying long-term changes in gully networks and volumes in dryland environments:The case of Northern Ethiopia[J].Geomorphology,2013,201: 254-263.

[87] MARARAKANYE N,LE ROUX J J.Gully location mapping at a national scale for South Africa[J].South African Geographical Journal,2012,94(2): 208-218.

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

ZHANG K L.A study on the distribution of erosion and sediment yeild on loess slope and the relationship between the distribution and rain characteristics[J].Jounal of Sediment Research,1991(4): 39-47.

[89] 唐克麗,鄭世清,席道勤,等.杏子河流域坡耕地的水土流失及其防治[J].水土保持通報(bào),1983,3(5): 43-48.

TANG K L,ZHENG S Q,XI D Q,et al.Water and soil lossess and their control in Hillslope land of Xingzihe Valley[J].Bulletin of Soil and Water Conservation,1983,3(5): 43-48.

[90] DEOLIVERIRA M A T.Slope geometry and gully erosion development: Bananal,Sao Paulo,Brazil[J].Zeitschrift Fur Geomorphologie,1990,34(4): 423-434.

[91] 秦 偉,朱清科,趙磊磊,等.基于RS和GIS 的黃土丘陵溝壑區(qū)淺溝侵蝕地形特征研究[J].農(nóng)業(yè)工程學(xué)報(bào),2010,26(6): 58-64.

QIN W,ZHU Q K,ZHAO L L,et al.Topographic characteristics of ephemeral gully erosion in loess hilly and gully region based on RS and GIS[J].Transactions of the Chinese Society of Agricultural Engineering,2010,26(6): 58-64.

[92] 張振國(guó),黃建成,焦菊英,等.安塞黃土丘陵溝壑區(qū)退耕地植物群落土壤抗蝕性分析[J].水土保持研究,2008,15(1): 28-31.

ZHANG Z G,HUANG J C,JIAO J Y,et al.Analysis on erosion resistance of different vegetation communities in abandoned lands in Ansai hilly gully loess region[J].Research of Soil and Water Conservation,2008,15(1): 28-31.

[93] CAPRA A,SPADA L C.Medium-term evolution of some ephemeral gullies in Sicily(Italy)[J].Soil and Tillage Research,2015,154: 34-43.

[94] SCHUMM S A,HADLEY R F.Arroyos and the semiarid cycle of erosion[J].American Journal of Science,1957,255(3): 164-174.

[95] SCHUMM S A.Geomorphic thresholds and complex response of drainage systems[J].Fluvial Geomorphology,1973: 299-310.

[96] PATTON P,SCHUMM S.Gully erosion,Northwestern Colorado-threshold phenomenon[J].Geology,1975,3(2): 88-90.

[97] BURKHAM D E.Hydrology of cornfield wash area and the effects of land-treatment practices,Sandoval County,New Mexico,1951-1960[J].Pediatric Research,1967,41(3): 201-201.

[98] BEGIN Z B,SCHUMM S A.Instability of alluvial valley floors: A method for its assessment[J].Transactions of the ASAE,1979,22: 347-350.

[99] VANDAELE K,POESEN J,GOVERS G,et al.Geomorphic threshold conditions for ephemeral gully incision[J].Geomorphology,1996,16(1): 161-173.

[100] DESMET P J J,POESEN J,GOVERS G,et al.Importance of slope gradient and contributing area for optimal prediction of the initiation and trajectory of ephemeral gullies[J].Catena,1999,37(3-4): 377-392.

[101] VANDEKERCKHOVE L,POESEN J,WIJDENES D O,et al.Topographical thresholds for ephemeral gully initiation in intensively cultivated areas of the Mediterranean[J].Catena,1998,33(3-4): 271-292.

[102] GOVERS G.Rill erosion on arable land in central Belgium: Rates,controls and predictability[J].Catena,1991,18(2):133-155.

[103] MORGEN R P C,Mngomezulu D.Threshold conditions for initiation of valley-side gullies in the middle veld of Swaziand[J].Catena,2003,50(2-4): 401-414.

[104] 胡 剛,伍永秋,劉寶元,等.東北漫川漫崗黑土區(qū)淺溝和切溝發(fā)生的地貌臨界模型探討[J].地理科學(xué),2006,26(4) :449-454.

HU G,WU Y Q,LIU B Y,et al.Geomorphic threshold model for ephemeral gully incision in rolling hills with black soil in Northeast China[J].Scientia Geographica Sinica,2006,26(4) :449-454.

[105] MOORE I D,BURCH G J,MACKENZIE D H.Topographic effects on the distribution of surface soil water and the location of ephemeral gullies[ J].Transactions of the ASEA,1988,31(4): 1098 -1107.

[106] 靳長(zhǎng)興.坡度在坡面侵蝕中的作用[J].地理研究,1996,15(3): 57-63.

JIN C X.The role of slope gradient on slope ersion[J].Geographical Research,1996,15(3): 57-63.

[107] 李斌兵,鄭粉莉,張 鵬,等.黃土高原丘陵溝壑區(qū)小流域淺溝和切溝侵蝕區(qū)的界定[J].水土保持通報(bào),2008,28(5): 16-20.

LI B B,ZHENG F L,ZHANG P,et al.Geomorphic threshold determination for ephemeral gully and gully erosion areas in the Loess Hilly Gully Region[J].Bulletin of Soil and Water Conservation,2008,28(5):16-20.

[108] TORRI D,POESEN J.A review of topographic threshold conditions for gully head development in different environments[J].Earth-Science Reviews,2014,130: 73-85.

[110] GUTIéRREZ A G,SCHNABEL S,CONTADOR F L.Gully erosion,land use and topographical thresholds during the last 60 years in a small rangeland catchment in SW Spain[J].Land Degradation and Development,2010,20(5): 535-550.