• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    雨強和坡度對黃土陡坡地淺溝形態(tài)特征影響的定量研究

    2017-07-12 18:45:38徐錫蒙鄭粉莉
    農(nóng)業(yè)工程學報 2017年11期
    關鍵詞:淺溝細溝格網(wǎng)

    徐錫蒙,鄭粉莉,2※,武 敏

    (1. 西北農(nóng)林科技大學水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室,楊凌 712100;2. 中國科學院水利部水土保持研究所,楊凌 712100;3 沈陽農(nóng)業(yè)大學水利學院,沈陽 110866)

    雨強和坡度對黃土陡坡地淺溝形態(tài)特征影響的定量研究

    徐錫蒙1,鄭粉莉1,2※,武 敏3

    (1. 西北農(nóng)林科技大學水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室,楊凌 712100;2. 中國科學院水利部水土保持研究所,楊凌 712100;3 沈陽農(nóng)業(yè)大學水利學院,沈陽 110866)

    淺溝形態(tài)特征是建立陡坡地坡面淺溝侵蝕預報模型的基礎。為了定量研究黃土陡坡地淺溝形態(tài)特征,在長8 m、寬2 m、深0.6 cm的試驗土槽上制作了雛形淺溝,設計了2個降雨強度(50、100 mm/h)和3個淺溝發(fā)生的典型坡度(15°、20°、25°),利用模擬降雨和徑流沖刷(10 L/min)相結(jié)合的試驗方法定量分析了黃土陡坡地的淺溝形態(tài)特征。結(jié)果表明:降雨強度和坡度的增加均加快了坡面淺溝侵蝕過程并使淺溝溝槽寬度和深度不斷增加,25°和100 mm/h降雨強度下的淺溝溝槽平均寬度和深度比15°和50 mm/h降雨強度下的分別增加1.40和0.61倍。根據(jù)測針板法得到的3 cm×10 cm精度的地表高程值數(shù)據(jù),在Surfer軟件中生成不同試驗處理下的地面數(shù)字高程模型(DEM, digital elevation model)及水流流路圖等,發(fā)現(xiàn)坡度的增加使兩側(cè)坡面細溝匯入淺溝溝槽的坡長增大,而降雨強度的增加則導致淺溝溝槽兩側(cè)坡面細溝匯入淺溝溝槽的坡長縮短,同時,溝道密度、地面割裂度和淺溝復雜度均隨著降雨強度和坡度的增加而呈現(xiàn)增大的趨勢,三者分別變化于0.74~1.48 m/m2、0.13~0.29和1.64~2.84之間,而不同降雨強度和坡度條件下淺溝溝槽寬深比變化于0.65~1.27之間。基于不同試驗處理下的DEM,根據(jù)相鄰格網(wǎng)關系在水平方向上計算方向?qū)?shù)后發(fā)現(xiàn),方向?qū)?shù)格網(wǎng)等值線圖可以有效地反映坡面淺溝和細溝的長度、表面積及侵蝕最嚴重的淺溝溝底位置。

    土壤;侵蝕;徑流;淺溝;形態(tài)特征;黃土陡坡地;模擬試驗;方向?qū)?shù)

    0 引 言

    淺溝侵蝕是在集中水流沖刷和人類耕作的共同作用下由坡面主細溝發(fā)育而來的典型線狀集中水流侵蝕,是坡面侵蝕過程中重要的侵蝕方式之一[1-3]。淺溝能夠被正常的耕作措施消除,但會在下一次徑流沖刷時發(fā)生在坡面的相同部位[4-5],隨著耕作和侵蝕的交替進行,淺溝可以發(fā)展成為固定的切溝[6-7]。淺溝侵蝕是細溝侵蝕向切溝侵蝕演變過程中的過渡侵蝕類型,在坡面形態(tài)演變和坡溝侵蝕系統(tǒng)中起著承上啟下的作用[8-10]。淺溝在黃土高原丘陵溝壑區(qū)第一副區(qū)和第二副區(qū)的分布面積占溝間地總面積的78%,其平均分布密度為每15~20 m寬坡面分布一條淺溝[11]。野外調(diào)查結(jié)果表明,淺溝侵蝕量占梁坡侵蝕產(chǎn)沙量的35%~46%以上。因此,淺溝侵蝕在黃土陡坡土壤侵蝕中占有重要地位,需要對坡面淺溝侵蝕進行深入研究并建立能夠預測淺溝侵蝕的坡面侵蝕預報模型。

    淺溝的形成和發(fā)展是坡面集中股流侵蝕的結(jié)果,受地形臨界條件的制約[12-14]。因此,淺溝常分布于坡度較陡的坡面上,而且一般多發(fā)生于具有一定匯水面積的坡面中部和中下部[14-15]。在淺溝發(fā)育過程中,影響淺溝侵蝕的因素主要有侵蝕動力因子(降雨[4]、上方匯流[16]、地下潛流[8])、地形因子(坡度、坡長、坡形等)[11-12]、地表植被、土壤、土地利用類型和人類活動[17]等。這些研究結(jié)果對認識淺溝侵蝕的發(fā)生和發(fā)展過程有重要作用。而在黃土丘陵區(qū)獨特的地形地貌特征條件下,淺溝的發(fā)生和發(fā)展使得原來平整的坡面成為獨特的瓦背狀地形,并形成了一系列新的傾斜面,增大了地面起伏度并增大了雨滴和徑流與坡面的接觸面積從而加強了侵蝕作用,同時瓦背狀地形使坡面徑流泥沙橫溢進入淺溝溝槽,最終增大了坡面輸沙量[18]。因此,研究淺溝形態(tài)特征對認識黃土陡坡地侵蝕產(chǎn)沙規(guī)律和建立相應的淺溝侵蝕預報模型具有重要意義。

    在美國農(nóng)業(yè)部等部門聯(lián)合研發(fā)的淺溝侵蝕預報模型(ephemeral gully erosion model,EGEM模型)[19]中,淺溝形態(tài)特征例如淺溝長度、深度、分叉及橫斷面形態(tài)在模型輸入和模型運算過程中被限定,使得其地區(qū)適應性較差[20-21],而通過借鑒AnnAGNPS和RSULE2等模型的基本原理預報淺溝侵蝕也需要以分析淺溝形態(tài)為基礎,提取水流流路并估算淺溝侵蝕量[22-23]。因此,研究淺溝形態(tài)特征可為建立適合陡坡地形的坡面淺溝侵蝕預報模型提供重要的理論支持。然而,目前關于淺溝形態(tài)特征的定量研究在野外測量方面多集中于對淺溝長度及相鄰淺溝間隔的描述[14-15],在室內(nèi)模擬方面主要集中于對淺溝溝槽的寬度和深度的定點測量[16-17],且測量的范圍和精度不足以將研究成果應用于侵蝕預報模型中。為此,本研究利用模擬降雨和徑流沖刷結(jié)合的試驗方法,用測針板法結(jié)合Surfer軟件生成了雨季徑流沖刷作用后的淺溝形態(tài)的地面數(shù)字高程模型(digital elevation model ,DEM),量化了黃土陡坡地的淺溝形態(tài)特征,以期為黃土陡坡地淺溝侵蝕預報模型的建立提供科學依據(jù)和理論基礎。

    1 材料與方法

    1.1 試驗材料

    試驗在黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室人工模擬降雨大廳進行。試驗選用側(cè)噴式降雨系統(tǒng),降雨高度16 m,能保證雨滴達到終點速度,降雨均勻度大于80%,雨滴直徑分布與天然降雨相似[24]。供試土槽為8 m長、2 m寬的固定式液壓升降鋼槽,其深度為0.6 m,坡度調(diào)節(jié)范圍0~30°,鋼槽底部每1 m 長排列4個孔徑為2 cm的排水孔以保證試驗過程中排水良好。鋼槽上端設置了放水水箱,能夠模擬自然條件下的上方匯水情況,鋼槽下端設置了集流裝置,用來收集徑流泥沙樣品。

    安塞縣(109°19′23″E、36°51′30″N)位于黃土高原丘陵溝壑區(qū)腹地,陡峭的地形條件、易于侵蝕的土質(zhì)以及不合理的人為開墾等造成了淺溝地貌在此區(qū)域內(nèi)的廣泛分布。因此本研究選取的供試土壤即為此地區(qū)的黃綿土,其顆粒組成為:砂粒(>50 μm)質(zhì)量分數(shù)占28.3%,粉砂粒(2~50 μm)質(zhì)量分數(shù)占58.1%,黏粒(<2 μm)質(zhì)量分數(shù)占13.6%,屬粉壤土,試驗土壤的采集樣地為當?shù)氐湫偷霓r(nóng)耕地,有明顯的犁底層,所以試驗土壤樣品采集時分為耕作層和犁底層2層,分別進行采集和保存,在裝填土槽時填入相應的土層。用重鉻酸鉀-外加熱法測定有機質(zhì)質(zhì)量分數(shù)為5.9 g/kg;水浸提法測定土壤pH值為7.9(土水比為1∶2.5)。

    1.2 試驗設計

    黃土坡面的淺溝侵蝕過程主要包括當年雨季前人為橫向犁耕及當年雨季徑流沖刷侵蝕過程。本研究通過在室內(nèi)人工建造淺溝發(fā)育初期的雛形,利用模擬降雨結(jié)合徑流沖刷研究了春季犁耕后當年雨季徑流沖刷作用后的淺溝形態(tài)特征。根據(jù)野外坡耕地調(diào)查和實際測量的春季人為橫向犁耕后的淺溝形態(tài)特征參數(shù),設計了本文中的淺溝雛形,以保證淺溝形態(tài)在研究區(qū)的典型性和代表性。在土槽距頂部2~8 m處用刮板制作了淺溝雛形,淺溝溝槽位于土槽中間,淺溝溝底與兩側(cè)溝坡高差12 cm,淺溝雛形的橫斷面為弧形(圖1)。制作淺溝雛形的刮板為2 m長的木板(與試驗土槽寬度相同),木板兩端與最低點的高程差為12 cm,并與淺溝雛形的弧形橫斷面形態(tài)一致。在每次試驗處理前用相同的刮板制作淺溝雛形,以保證雛形淺溝形態(tài)的一致。

    圖1 淺溝雛形形態(tài)Fig.1 Initial ephemeral gully morphology

    根據(jù)黃土高原短歷時、高強度侵蝕性降雨標準[25](即I5=1.52 mm/min,5 min瞬時雨量為7.6 mm),設計降雨強度為50和100 mm/h?;邳S土高原淺溝分布較廣且發(fā)育活躍的典型坡度(淺溝侵蝕一般發(fā)育在18°~35°的坡面上,平均坡度為23°)[11,18],并結(jié)合15°是淺溝發(fā)生的臨界坡度的研究結(jié)果[26],本研究設計了15°、20°和25°共3個陡坡地坡度[3,27],每次試驗歷時70 min。為了模擬野外實際情況下淺溝溝槽發(fā)育過程中的侵蝕動力情況[17],根據(jù)野外實際情況下的上方匯水面積及產(chǎn)流系數(shù),設計上方匯水流量為10 L/min。各試驗處理重復2次,所有試驗數(shù)據(jù)取平均值并計算標準偏差(SD)。

    1.3 試驗步驟

    為保持土壤原有結(jié)構(gòu)免遭破壞,不對試驗土壤采取過篩和研磨處理。填土時,首先用紗布填充試驗土槽底部的排水孔,在土槽底部填10 cm細沙,以保證良好的透水性。然后在沙層上覆蓋紗布,并裝填40 cm的黃綿土。填土時將40 cm的土層分為耕層和犁底層,其中耕層深度為20 cm,容重為1.06~1.10 g/cm3,犁底層深度為20 cm,容重為1.25 g/cm3。填土時采用分層裝土,每次裝土深度為5 cm,以保證試驗土槽裝土的均勻性。每填完一層后,用齒耙將土層表面耙松,再填裝下一層土壤,以保證2個土層能夠很好地接觸。為了減小邊界效應的影響,在填土時將試驗土槽的四周邊壁壓實。在填裝耕層土壤后,用傳統(tǒng)的犁耕方式對土槽進行橫向犁耕,并用刮板將距土槽頂部2~8 m的表層土壤刮去,制作淺溝雛形。

    為了保證試驗前坡面土壤含水量的一致性,選用30 mm/h降雨強度進行預降雨至坡面產(chǎn)流,然后靜置12 h后開始正式降雨。正式降雨開始前,對降雨強度和匯水流量進行率定,當實測降雨強度和匯水流量與目標降雨強度和匯水流量的差值小于5%方可進行正式降雨和沖刷試驗。

    由于測針板法精度可靠,操作簡單便且實用性強,在溝蝕監(jiān)測中得到了廣泛的應用[3,28]。因此,在每場降雨結(jié)束后,采用測針板法測量坡面淺溝及細溝形態(tài)。本研究中采用的測針板寬度1.6 m,由53根測針構(gòu)成,測針間距為3 cm。試驗土槽寬2 m,當完成某一橫斷面1 m寬的測量后,將測針板移動測量同一斷面另1 m寬的坡面形態(tài)。測量完一個斷面的高程值后,將測針板沿坡長方向移動10 cm測量下一個斷面的高程值。這樣,長8 m、寬2 m的試驗土槽共被分成80個橫斷面,在每個橫斷面上每隔3 cm讀取一個數(shù)值,共讀取67個,一次試驗中共測量8 m×2 m試驗土槽的5 360個三維坐標數(shù)據(jù),并將數(shù)據(jù)導入Surfer軟件,生成DEM。

    1.4 淺溝形態(tài)指標

    溝道密度(ρ)是指單位研究區(qū)域內(nèi)所有溝道的總長度,能夠反映坡面的破碎程度,計算式為

    式中Li為第i條溝道及其分叉的總長度,m;A為研究坡面的表面積,m2。

    地面割裂度(D)是指單位研究區(qū)域內(nèi)所有溝道的平面面積之和,為無量綱單位,能夠更為全面地反映坡面的破碎程度及淺溝侵蝕強度,計算式為

    式中Ai為第i條溝道的表面積,m2。

    淺溝復雜度(c)是指淺溝溝槽及其分叉的總長度與對應的垂直有效長度的比值,能夠反映向淺溝溝槽匯水的細溝水流流路的豐富度,計算式為

    式中Le為淺溝溝槽及其分叉的總長度,m;Lev為淺溝溝槽的垂直有效長度,m。

    淺溝溝槽寬深比(RWD)是指淺溝溝槽寬度與對應深度的比值,該參數(shù)是無量綱參數(shù),可以反映淺溝溝槽形狀的變化,計算式為

    式中Wj為第j個監(jiān)測點處的淺溝溝槽寬度,m;Dj為第j個監(jiān)測點處的淺溝溝槽深度,m。

    1.5 方向?qū)?shù)計算

    將高程點數(shù)據(jù)導入Surfer軟件,生成DEM,并根據(jù)相鄰格網(wǎng)的關系進行方向?qū)?shù)計算,獲得DEM格網(wǎng)在水平X軸正方向上的坡度及坡度變化率[29]。本文選取的3種方向?qū)?shù)為:

    一階導數(shù):計算DEM格網(wǎng)表面沿水平方向的坡度,等于坡度矢量和指定方向上的單位矢量的乘積。在某一特定格網(wǎng)上,沿水平軸方向,格網(wǎng)處坡度上升即為正值,坡度下降為負值,可用于判定坡面溝道發(fā)生的位置并量化地形起伏的劇烈程度。

    二階導數(shù):計算DEM格網(wǎng)表面沿水平方向的坡度的變化率,為一階導數(shù)的方向?qū)?shù)。經(jīng)過“二階導數(shù)”分析生成的格網(wǎng)文件可以顯示在水平方向上坡度變化率的等值線。在某一個特定的格網(wǎng)上,水平方向上的二階導數(shù)值越大,說明格網(wǎng)處于地表形態(tài)轉(zhuǎn)折處,例如淺溝溝槽底部和溝壁等位置。

    方向曲率Ks:計算沿水平方向的剖面的切線傾角的變化率。方向曲率是對于一個表面f(x,y)在某一方向上的切面傾角變化率的絕對值,均為正值。它與二次導數(shù)有些類似,其極值能展示水平方向上地表形態(tài)的轉(zhuǎn)折位置。

    為了比較不同試驗處理的地表形態(tài)特征及侵蝕空間分布,對降雨后坡面DEM格網(wǎng)進行方向?qū)?shù)計算得到一階導數(shù)、二階導數(shù)和方向曲率格網(wǎng)文件,并創(chuàng)建了各方向?qū)?shù)文件在水平方向上的剖面線并統(tǒng)計方向?qū)?shù)值,對比分析不同試驗處理下坡面侵蝕強度的空間分布。

    2 結(jié)果與分析

    2.1 不同降雨強度和坡度條件下淺溝形態(tài)對比

    根據(jù)測針板法得到的雨后地表形態(tài)的三維坐標數(shù)據(jù),在Surfer軟件中分別制作了各試驗處理的坡面等高線圖、水流流路圖、格網(wǎng)矢量地圖和3D曲面圖(圖2)。

    在降雨強度為50 mm/h的條件下,坡度為15°時(圖2a),坡面土壤侵蝕主要發(fā)生在淺溝溝槽,經(jīng)過70 min,溝頭溯源侵蝕活躍,淺溝溝頭雖未達到坡頂位置,但其上方也出現(xiàn)了斷續(xù)的細溝下切溝頭,淺溝兩側(cè)的坡面僅有很淺的非固定水流流路出現(xiàn),發(fā)生了輕微的細溝侵蝕。而坡度為20°時(圖2b),淺溝溝頭已溯源侵蝕至坡頂位置,淺溝溝槽寬度和深度略有增加,淺溝溝槽兩側(cè)的細溝侵蝕發(fā)育明顯,細溝深度和寬度均增大,已形成了固定的水流流路。當坡度增加到25°時(圖2c),淺溝溝槽與兩側(cè)細溝進一步加寬和加深,淺溝溝槽深度已接近犁底層位置(深度為(18.8±0.6) cm),其中部分淺溝溝槽深度已經(jīng)超過犁底層(深度>20 cm)。對比3個坡度下的淺溝形態(tài)可知,在相同降雨強度和降雨歷時條件下,隨著坡度的增大,淺溝溝槽的溝頭溯源侵蝕、溝壁崩塌和溝底下切侵蝕速率增加,同時淺溝溝槽兩側(cè)坡面上的細溝發(fā)育也更加完善,土壤侵蝕加劇。

    圖2 各試驗處理下坡面水流流路和淺溝形態(tài)Fig.2 Flow path and ephemeral gully morphology in different treatments

    在降雨強度為100 mm/h,坡度為15°時(圖2d),坡面侵蝕主要為淺溝溝槽侵蝕和兩側(cè)坡面上固定流路的細溝侵蝕,與相同坡度50 mm/h降雨強度條件下的坡面侵蝕狀況相比(圖2a),淺溝溝槽寬度略有增加,淺溝溝槽深度由于受到犁底層的限制變化不大,而兩側(cè)坡面上的細溝侵蝕明顯加劇,細溝寬度和細溝深度均明顯增大,匯入淺溝溝槽的細溝數(shù)量也增多。坡度為20°時(圖2e),淺溝溝槽深度已超過耕層深度到達犁底層(平均深度為(24.2±2.8) cm),兩側(cè)坡面細溝網(wǎng)發(fā)育也更為完善。坡度為25°時(圖2f),淺溝溝槽深度較20°條件下略有增加,但溝壁崩塌侵蝕明顯加快,淺溝溝槽寬度增加明顯(平均寬度為(35.0±2.4) cm),幾乎為20°條件下淺溝溝槽寬 (平均寬度為(16.9±0.8) cm)的2倍,而由于分叉細溝的合并,兩側(cè)細溝數(shù)量減少但細溝寬度和深度明顯增加,同時,坡度的增加使得兩側(cè)的細溝水流無法在8 m的坡長內(nèi)匯入淺溝溝槽而直接匯入出水口。

    不同處理下淺溝溝槽的形態(tài)差別是由淺溝發(fā)育過程決定的。淺溝溝槽侵蝕過程包括溝頭溯源侵蝕、溝底下切侵蝕和溝壁擴張侵蝕3種過程(圖3),隨著降雨強度和坡度的變化,淺溝溝槽的溝頭溯源侵蝕、溝底下切侵蝕和溝壁崩塌侵蝕速率發(fā)生了變化[17],在70 min的降雨后停止在了不同的發(fā)育階段,并決定了淺溝溝槽的最終形態(tài)。在降雨強度為50 mm/h的試驗條件下,坡度為15°的坡面在70 min降雨歷時內(nèi)主要以淺溝溝頭溯源侵蝕為主;坡度為20°的坡面已由溝頭溯源侵蝕為主的階段轉(zhuǎn)入了以溝底下切侵蝕為主的階段;而坡度為25°的坡面處于以溝底下切侵蝕為主的階段,同時存在少量的溝壁崩塌。在100 mm/h降雨強度下,淺溝侵蝕過程明顯加快。坡度為15°的坡面在70 min降雨結(jié)束后就進入了溝底下切侵蝕和溝壁崩塌侵蝕并存的階段,坡度為20°的坡面則處于淺溝溝槽溝底下切侵蝕的后期,淺溝溝槽深度明顯加大,而坡度為25°的坡面則在70 min的降雨后結(jié)束了溝底下切侵蝕階段,進入了以溝壁崩塌擴張侵蝕為主階段,淺溝溝槽寬度明顯增大。圖3展示了100 mm/h和25°條件下的淺溝侵蝕過程,并隨著降雨歷時的延長呈現(xiàn)出分別以溝頭溯源侵蝕、溝底下切侵蝕和溝壁擴張侵蝕為主的3個淺溝發(fā)育階段。

    圖3 降雨強度100 mm/h和坡度25°條件下的淺溝侵蝕過程Fig.3 Ephemeral gully erosion processes under condition of 100 mm/h and 25°

    對比各試驗處理的水流流路圖可知(圖2),在淺溝溝槽發(fā)育的同時,淺溝溝槽兩側(cè)坡面上的細溝侵蝕也明顯受坡度和雨強的影響。在相同的降雨強度下,隨著坡度的增加,原來的斷續(xù)細溝逐漸形成連續(xù)細溝并連接成細溝網(wǎng)絡,細溝寬度和深度不斷增加[30]。此外,細溝匯入淺溝溝槽所需坡長逐漸增大,匯入淺溝溝槽時與淺溝溝槽所夾得銳角角度逐漸變小。例如,在100 mm/h的降雨強度下,15°坡面左側(cè)的細溝每隔2~3 m坡長匯入淺溝溝槽,20°坡面左側(cè)的細溝每隔4~5 m坡長匯入淺溝溝槽,而25°坡面在8 m的坡長內(nèi)左側(cè)只有1條細溝匯入淺溝溝槽;類似的情況也發(fā)生在右側(cè)坡面,15°坡面右側(cè)的細溝在坡中部匯入淺溝溝槽,同時在坡面中部也形成了一條沒有匯入淺溝溝槽的細溝,20°坡面右側(cè)的細溝在坡面下部匯入淺溝溝槽并在坡面下部形成了一條沒有匯入淺溝槽的細溝,而25°坡面右側(cè)的細溝直接匯入了出水口,沒有匯入淺溝溝槽。而對比相同坡度不同降雨強度下的細溝形態(tài)可知,降雨強度的增大促進了兩側(cè)細溝網(wǎng)的發(fā)育,使細溝變寬變深,而且縮短了細溝向淺溝溝槽匯水所需的坡長,即降雨強度的增大使得坡面兩側(cè)細溝網(wǎng)絡豐富度和復雜度增加。

    土壤侵蝕直接作用的結(jié)果是使地面破碎化。降雨強度和地面坡度的增加,導致淺溝侵蝕和細溝侵蝕加劇,地表更加破碎,但降雨強度和坡度對坡面侵蝕影響的方式和程度不同,各試驗處理的最終地表形態(tài)也不盡同[31-32],因此,需要對坡面淺溝形態(tài)特征進行量化,從坡面地形特征參數(shù)的變化方面探究坡面侵蝕的強弱。

    2.2 不同降雨強度和坡度下淺溝形態(tài)特征參數(shù)

    坡面淺溝形態(tài)特征參數(shù)可以量化反映坡面土壤侵蝕方式及侵蝕程度的強弱?;诟髟囼炋幚淼腄EM數(shù)據(jù),利用式(1)~(4)計算了各試驗處理下的淺溝形態(tài)特征參數(shù)(表1),量化了降雨強度和坡度對淺溝形態(tài)特征參數(shù)的影響??傮w來看,降雨強度和坡度的增大均加快了坡面淺溝侵蝕過程,25°和100 mm/h條件下的淺溝平均寬度和深度分別較15°和50 mm/h條件下的平均寬度和深度增大了1.40和0.61倍。

    表1 不同降雨強度和坡度下淺溝形態(tài)特征參數(shù)Table 1 Ephemeral gully morphological indicators under different rainfall intensities and slope gradients

    50 mm/h降雨強度下,25°坡面的淺溝溝槽平均寬度與15°和20°坡面相比有顯著差異,為15°坡面淺溝溝槽平均寬度的1.23倍;而當降雨強度為100 mm/h時,15°坡面的淺溝溝槽平均寬度和平均深度與50mm/h雨強下的平均寬度和平均深度均無顯著差異, 20°坡面的淺溝溝槽平均深度則顯著增加了35.2%,25°坡面的淺溝溝槽平均寬度和平均深度與50 mm/h雨強下的平均寬度和平均深度均存在顯著性差異,分別增加了94.4%和57.3%。這說明在極端降雨和陡峭地形條件下坡面淺溝溝槽發(fā)育速度明顯加快,這也是造成黃土陡坡地土壤侵蝕嚴重的重要原因[4]。

    隨著降雨強度和坡度的增大,溝道密度、地面割裂度和淺溝復雜度均呈現(xiàn)增大趨勢,三者分別變化于0.74~1.48 m/m2、0.13~0.29和1.64~2.84之間。溝道密度和地面割裂度不僅由淺溝溝槽的形態(tài)所決定,也與坡面兩側(cè)細溝形態(tài)及細溝網(wǎng)發(fā)育情況有關[30,33];而淺溝復雜度主要是由匯入淺溝溝槽的細溝數(shù)量及形態(tài)決定。50 mm/h降雨強度下,隨著坡度的增加,雖然淺溝溝槽的平均寬度和平均深度變化不大,但坡面兩側(cè)細溝網(wǎng)發(fā)育卻十分明顯,細溝總長度、寬度和深度都明顯增大,因此溝道密度、地面割裂度和淺溝復雜度均不斷增大。而100 mm/h降雨強度下,坡度由15°增加到20°時,溝道密度、地面割裂度和淺溝復雜度均顯著增大,而當坡度從20°增加到25°時,三者的數(shù)值變化則相對較小,說明在坡度為20°時坡面兩側(cè)細溝已經(jīng)發(fā)育完善,坡度增加并沒有使細溝形態(tài)特征發(fā)生明顯變化。

    淺溝溝槽寬深比可以反映淺溝溝槽形狀的變化。不同降雨強度和坡度下,淺溝溝槽寬深比為0.65~1.27,其隨降雨強度和坡度的增加而增大,這主要與淺溝溝槽的發(fā)育階段和主導侵蝕過程有關(圖3)。在100 mm/h的降雨強度下,20°坡面的淺溝溝槽寬深比小于15°和25°坡面的淺溝溝槽寬深比,主要原因是20°坡面條件下的淺溝溝槽在降雨結(jié)束時正處于下切侵蝕的活躍期,淺溝溝槽深度大,而溝槽寬度則與15°坡面的淺溝溝槽寬度差異較小,最終使得20°坡面的淺溝溝槽寬深比最小。在黃土陡坡地上,淺溝寬深比隨著降雨侵蝕的進行會逐漸穩(wěn)定在1左右[34],而不同的降雨強度和坡度條件使得淺溝溝槽在70 min的降雨結(jié)束后停留在以不同侵蝕過程為主導的各個階段,造成淺溝溝槽形態(tài)的不同。

    2.3 坡面侵蝕空間分布與淺溝橫斷面形態(tài)

    為了進一步說明淺溝溝槽及坡面侵蝕形態(tài)的變化,對比分析坡面侵蝕強度的空間分布,將各場次降雨后利用測針板測量獲得的高程數(shù)據(jù)導入Surfer軟件,生成了3個試驗處理(a:降雨強度為50 mm/h和坡度為20°;b:降雨強度為50 mm/h和坡度為25°;c:降雨強度為100 mm/h和坡度為25°)雨后的坡面DEM,并利用式(5)~(7)在水平方向上根據(jù)DEM相鄰格網(wǎng)數(shù)值關系進行方向?qū)?shù)計算得到了一階導數(shù),二階導數(shù)和方向曲率格網(wǎng)表面的等值線圖(圖4)。

    雨后坡面DEM等值線圖展示了不同試驗處理下坡面地形狀況(圖4);一階導數(shù)等值線圖可以反映水平方向上坡度分布狀況,在水平方向上坡度上升該格網(wǎng)點值為正值,坡度下降格網(wǎng)點的值為負值,而等值線密集的地方即為坡度變化劇烈的地方。對比坡面水流流路圖(圖2)可知,等值線密集部位與水流流路即溝底位置重合,可用于估算細溝及淺溝長度。二階導數(shù)等值線圖(圖4)反映了水平方向上坡度的變化率,在特定的格網(wǎng)處,坡度上升,坡度變化率值為正,坡度下降,坡度變化率值為負,圖中等值線主要密集分布在坡度值變化較快的地方,也就是坡面淺溝溝槽和細溝溝槽位置,因此等值線密集的區(qū)域可用于展示坡面細溝及淺溝的平面分布位置及所占表面積。方向曲率等值線圖(圖4)展示了在水平方向上切面傾角變化率絕對值的變化情況,所有的格網(wǎng)值均為正值,圖中等值線密集的地方為溝槽在水平方向上的曲率最大值,能夠詳盡展示溝底所在的位置,同時也反映了坡面侵蝕最嚴重的部位。

    對比3個不同試驗處理的雨后坡面地形和各方向?qū)?shù)等值線圖可以發(fā)現(xiàn),隨著坡度和雨強的增大,淺溝溝槽處的等值線彎曲度增大,且淺溝溝槽兩側(cè)的細溝也明顯增多加長,細溝處等高線彎曲度也增大(圖4)。對比一階導數(shù)等值線圖可知(圖4),雨強和坡度的增加促進了坡面水流流路的發(fā)育,特別是淺溝溝槽兩側(cè)細溝水流流路逐漸加長并匯入了淺溝溝槽,試驗處理A中淺溝溝槽兩側(cè)的斷續(xù)細溝沒有完全連接,而B處理中細溝水流流路連接形成連續(xù)細溝,C處理中細溝水流流路則更加明顯。二階導數(shù)等值線圖則顯示,坡面上細溝和淺溝表面積隨著坡度和降雨強度的增大也明顯增大。方向曲率等值線圖展示了侵蝕最嚴重的溝底位置,對比不同處理的方向曲率等值線圖可知,隨著坡度從20°增加到25°,試驗處理B中不僅淺溝溝槽溝底下切侵蝕強度明顯增大,淺溝溝槽兩側(cè)的細溝溝底面積也明顯增多,而且坡面右側(cè)還出現(xiàn)了一條連續(xù)的細溝;而隨著降雨強度從50增大到100 mm/h,試驗處理C中淺溝溝槽兩側(cè)的細溝發(fā)育明顯加強,侵蝕嚴重的溝底分布面積明顯加大。

    圖4 雨后坡面DEM與水平方向上各方向?qū)?shù)等值線圖Fig.4 Contour map of slope surface after rain and its directional derivatives grids in horizon direction

    為展示坡面淺溝和細溝橫斷面形態(tài),基于圖3的等值線文件創(chuàng)建了各試驗處理3 m坡長處的剖面線文件,并統(tǒng)計了該剖面線上一階導數(shù)、二階導數(shù)、方向曲率和降雨后坡面DEM的格網(wǎng)表面值(圖5)。

    圖5 雨后坡面DEM與各方向?qū)?shù)格網(wǎng)在3 m坡長處的剖面圖Fig.5 Profile map of digital elevation model slope surface after rain and directional derivatives grids at 3 m of slope length

    對比不同處理雨后坡面DEM的剖面圖可知(圖5a),在3 m坡長處,隨著降雨強度與坡度的增大,淺溝溝槽深度變大,淺溝溝槽形態(tài)逐漸由寬淺型向?qū)捝畋葹?變化,淺溝溝槽兩側(cè)細溝也增多,細溝侵蝕強度變大。50 mm/h降雨強度下,20°坡面淺溝溝槽寬80 cm、深24 cm,25°坡面的淺溝溝槽寬75 cm、深25 cm,均呈現(xiàn)寬淺型,而在100 mm/h和25°條件下,淺溝溝槽寬度為54 cm,深為38 cm,同時兩側(cè)坡面上也有3條細溝發(fā)育,左側(cè)槽邊有一條寬為6 cm深為5 cm的細溝,其旁有一條寬為25 cm,最大深度為8 cm細溝,而右側(cè)坡面細溝寬度為13 cm,深度為5 cm(圖5a)??梢姡涤陱姸群推露鹊脑龃蟠龠M了坡面淺溝溝槽和兩側(cè)細溝的發(fā)育,而淺溝溝槽的深度和寬度明顯大于兩側(cè)細溝侵蝕,是坡面侵蝕最嚴重的部位。對比各處理的不同方向?qū)?shù)的格網(wǎng)值可知,一階導數(shù)表示的是雨后坡面DEM格網(wǎng)在水平方向上的坡度,所以在溝底位置一階導數(shù)的格網(wǎng)表面值為0(圖5b);二階導數(shù)值表示的是坡度的變化率,因此,原始坡面溝底位置即一階導數(shù)值為0的地方,二階導數(shù)值出現(xiàn)峰值(圖5c);類似的,方向曲率格網(wǎng)表面值在溝底位置最大,且所有格網(wǎng)表面值均為正值(圖5d)。由各方向?qū)?shù)剖面圖可知,各方向?qū)?shù)格網(wǎng)表面值可反映坡面的侵蝕分布狀況,展示坡面上細溝和淺溝溝底的位置。

    3 結(jié) 論

    在長8 m、寬2 m、深0.6 m的試驗土槽上制作了淺溝發(fā)育初期的雛形,利用模擬降雨和徑流沖刷(10 L/min)相結(jié)合的試驗方法,定量研究了不同降雨強度和坡度條件下坡面淺溝形態(tài)特征,結(jié)論如下:

    1)降雨強度和坡度的增大均加快了坡面淺溝侵蝕過程,25°和100 mm/h條件下的淺溝平均寬度和深度分別較15°和50 mm/h條件下的平均寬度和深度增大1.40和0.61倍。坡度的增加使兩側(cè)細溝匯入淺溝溝槽所需的坡長變大,細溝匯入淺溝溝槽時與淺溝溝槽所夾得銳角角度也逐漸變??;而降雨強度的增大則縮短了細溝向淺溝溝槽匯水所需的坡長。

    2)受淺溝溝槽形態(tài)和坡面兩側(cè)細溝網(wǎng)發(fā)育情況的影響,溝道密度、地面割裂度和淺溝復雜度均隨著降雨強度和坡度的增大呈現(xiàn)增大趨勢,三者分別變化于0.74~1.48 m/m2、0.13~0.29和1.64~2.84之間。而淺溝寬深比受淺溝發(fā)育過程中不同主導侵蝕過程的影響,在不同降雨強度和坡度的試驗條件下變化于0.65~1.27之間。

    3)方向?qū)?shù)格網(wǎng)等值線圖可以反映坡面淺溝和細溝的長度、表面積及侵蝕最嚴重的溝底位置,剖面線也可以展示坡面上細溝和淺溝溝槽底部的位置。

    目前,有關淺溝形態(tài)特征及建立基于淺溝形態(tài)的坡面侵蝕預報模型研究有待加強。

    [1] 朱顯謨. 黃土區(qū)土壤侵蝕的分類[J]. 土壤學報,1956,4(2):99-115. Zhu Xianmo. Soil erosion classification in the loessial region[J]. Acta Pedologica Sinica, 1956, 4(2): 99-115. (in Chinese with English abstract)

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

    [3] 鄭粉莉,徐錫蒙,覃超. 溝蝕過程研究進展[J]. 農(nóng)業(yè)機械學報,2016,47(8):48-59,116. Zheng Fenli, Xu Ximeng, Qin Chao. A review of gully erosion processes research[J]. Transaction of the Chinese Society for Agricultural Machinery, 2016, 47(8): 48-59, 116. (in Chinese with English abstract)

    [4] Gong J G, Jia Y W, Zhou Z H, et al. An experimental study on dynamic processes of ephemeral gully erosion in loess landscapes[J]. Geomorphology, 2011, 125(1): 203-213.

    [5] 崔明,蔡強國,張永光,等. 漫崗黑土區(qū)坡耕地中雨季淺溝發(fā)育機制[J]. 農(nóng)業(yè)工程學報,2007,23(8):59-65. Cui Ming, Cai Qiangguo, Zhang Yongguang, et al. Development of ephemeral gully during rainy season in the slope land in rolling-hill black-soil region of Northeast China[J]. Transactions of the Chinese Society of AgriculturalEngineering (Transactions of the CSAE), 2007, 23(8): 59-65. (in Chinese with English abstract)

    [6] Di Stefano C, Ferro V, Pampalone V, et al. Field investigation of rill and ephemeral gully erosion in the Sparacia experimental area, South Italy[J]. Catena, 2013, 101: 226-234.

    [7] 武敏,鄭粉莉. 淺溝侵蝕過程及預報模型研究進展[J]. 水土保持研究,2004,11(4):113-116. Wu Min, Zheng Fenli. Research progress on ephemeral gully erosion process and its prediction model[J]. Research of Soil and Water Conservation, 2004, 11(4): 113-116. (in Chinese with English abstract)

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

    [9] Poesen J, Nachtergaele J, Verstraeten G, et al. Gully erosion and environmental change: Importance and research needs[J]. Catena, 2003, 50(2): 91-133.

    [10] 胡剛,伍永秋,劉寶元,等. 東北漫崗黑土區(qū)淺溝侵蝕發(fā)育特征[J]. 地理科學,2009,29(4):545-549. Hu Gang, Wu Yongqiu, Liu Baoyuan, et al. Growth characteristic of ephemeral gully in rolling hills of black soils in northeast China[J]. Scientia Geographica Sinica, 2009, 29(4): 545-549. (in Chinese with English abstract)

    [11] 張科利,唐克麗,王斌科. 黃土高原坡面淺溝侵蝕特征值的研究[J]. 水土保持學報,1991,5(2):8-13. Zhang Keli, Tang Keli, Wang Binke. A study on characteristic value of shallow gully erosion on slope farmland in the Loess Plateau[J]. Journal of Soil and Water Conservation, 1991, 5(2): 8-13. (in Chinese with English abstract)

    [12] 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. [13] Maugnard A, Van Dyck S, Bielders C L. Assessing the regional and temporal variability of the topographic threshold for ephemeral gully initiation using quantile regression in Wallonia (Belgium)[J]. Geomorphology, 2014, 206: 165-177.

    [14] Cheng H, Zou X, Wu Y, et al. Morphology parameters of ephemeral gully in characteristics hillslopes on the Loess Plateau of China[J]. Soil and Tillage Research, 2007, 94(1): 4-14.

    [15] 秦偉,朱清科,趙磊磊,等. 基于RS 和GIS 的黃土丘陵溝壑區(qū)淺溝侵蝕地形特征研究[J]. 農(nóng)業(yè)工程學報,2010,26(6):58-64. Qin Wei, Zhu Qingke, Zhao Leilei, 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 (Transactions of the CSAE), 2010, 26(6): 58-64. (in Chinese with English abstract)

    [16] 鄭粉莉,武敏,張玉斌,等. 黃土陡坡裸露坡耕地淺溝發(fā)育過程研究[J]. 地理科學,2006,26(4):438-442. Zheng Fenli, Wu Min, Zhang Yubin, et al. Ephemeral gully development process at loess steep hillslope[J]. Scientia Geographica Sinica, 2006, 26(4): 438-442. (in Chinese with English abstract)

    [17] 郭明明,王文龍,李建明,等. 黃土區(qū)坡耕地耕作對淺溝徑流產(chǎn)沙及其形態(tài)發(fā)育特征的影響[J]. 農(nóng)業(yè)工程學報,2015,31(15):114-123. Guo Mingming, Wang Wenlong, Li Jianming, et al. Effect of tillage on runoff and sediment yields and morphology development characteristic of ephemeral gully in loessial region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(15): 114-123. (in Chinese with English abstract)

    [18] 姜永清,王占禮,胡光榮,等. 瓦背狀淺溝分布特征分析[J].水土保持研究,1999,6(2):181-184. Jiang Yongqing, Wang Zhanli, Hu Guangrong, et al. Distribution features of shallow gully[J]. Research of Soil and Water Conservation, 1999, 6(2): 181-184. (in Chinese with English abstract)

    [19] Woodward D E. Method to predict cropland ephemeral gully erosion[J]. Catena, 1999, 37(3): 393-399.

    [20] Capra A, Mazzara L M, Scicolone B. Application of the EGEM model to predict ephemeral gully erosion in Sicily, Italy[J]. Catena, 2005, 59(2): 133-146.

    [21] Nachtergaele J, Poesen J, Sidorchuk A, et al. Prediction of concentrated flow width in ephemeral gully channels[J]. Hydrological Processes, 2002, 16(10): 1935-1953.

    [22] Taguas E V, Yuan Y, Bingner R L, et al. Modeling 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.

    [23] Daggupati P, Sheshukov A Y, Douglas-Mankin K R. Evaluating ephemeral gullies with a process-based topographic index model[J]. Catena, 2014, 113: 177-186.

    [24] 鄭粉莉,趙軍. 人工模擬降雨大廳及模擬降雨設備簡介[J].水土保持研究,2004,11(4):177-178.

    [25] 周佩華,王占禮. 黃土高原土壤侵蝕暴雨標準[J]. 水土保持通報,1987,7(1):38-44. Zhou Peihua, Wang Zhanli. Soil erosion storm rainfall standard in the Loess Plateau[J]. Bulletin of Soil and Water Conservation, 1987, 7(1): 38-44. (in Chinese with English abstract)

    [26] 江忠善,鄭粉莉,武敏. 中國坡面水蝕預報模型研究[J]. 泥沙研究,2005(4):1-6. Jiang Zhongshan, Zheng Fenli, Wu Min. Prediction model of water erosion on hillslopes[J]. Journal of Sediment Research, 2005( 4): 1-6. (in Chinese with English abstract)

    [27] 王紅雷,王秀茹,王希,等. 采用 SCS-CN 水文模型和GIS 確定雨水集蓄工程的位置[J]. 農(nóng)業(yè)工程學報,2012,28(22):108-114. Wang Honglei, Wang Xiuru, Wang Xi, et al. Identification of potential rainwater harvesting sites using SCS-CN and GIS[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(22): 108-114. (in Chinese with English abstract)

    [28] Vinci A, Brigante R, Todisco F, et al. Measuring rill erosion by laser scanning[J]. Catena, 2015, 124: 97-108.

    [29] 白世彪,王建,常直楊. Surfer 10地學計算機制圖[M]. 北京:科學出版社,2012:205-211.

    [30] Shen H, Zheng F, Wen L, et al. An experimental study of rill erosion and morphology[J]. Geomorphology, 2015, 231: 193-201.

    [31] 沈海鷗,鄭粉莉,溫磊磊,等. 降雨強度和坡度對細溝形態(tài)特征的綜合影響[J]. 農(nóng)業(yè)機械學報,2015,46(7):162-170. Shen Haiou, Zheng Fenli, Wen Leilei, et al. Effects of rainfall intensity and slope gradient on rill morphological characteristics[J]. Transactions of the Chinese Society of Agricultural Machinery (Transactions of the CSAM), 2015, 46(7): 162-170. (in Chinese with English abstract)

    [32] Berger C, Schulze M, Rieke-Zapp D, et al. Rill development and soil erosion: A laboratory study of slope and rainfall intensity[J]. Earth Surface Processes and Landforms, 2010, 35(12): 1456-1467.

    [33] Brunton D A, Bryan R B. Rill network development and sediment budgets[J]. Earth Surface Processes and Landforms, 2000, 25(7): 783-800.

    [34] 龔家國,周祖昊,賈仰文,等. 黃土區(qū)淺溝侵蝕溝槽發(fā)育及其水流水力學基本特性模擬實驗研究[J]. 水土保持學報,2010,24(5):92-96. Gong Jiaguo, Zhou Zuhao, Jia Yangwen, et al. Simulation experiment of ephemeral gully erosion on basis hydraulics parameters of concentrated flow and erosion morphology in loess area[J]. Journal of Soil and Water Conservation, 2010, 24(5): 92-96. (in Chinese with English abstract)

    Quantification study of rainfall intensity and slope gradient impacts on ephemeral gully morphological characteristic on steep loessial hillslope

    Xu Ximeng1, Zheng Fenli1,2※, Wu Min3
    (1. State Key Laboratory of Soil Erosion and Dryland Farming on Loess Plateau, Institute of Soil and Water Conservation, Northwest A&F University, Yangling 712100, China; 2. Institute of Soil and Water Conservation CAS & MWR, Yangling 712100, China; 3.College of Water Conservancy, Shenyang Agriculture University, Shenyang 110866, China)

    Ephemeral gully is formed in the cycle of concentrated flow and tillage practices, and contributed a lot to the sediment yield on the loessial hillslope as well as in the watershed. Ephemera gully morphology is the base of establishing ephemeral gully erosion prediction model on the steep loessial hillslope. Thus, to quantify the ephemeral gully morphological characteristics on steep loessial hillslope in different rainfall intensities and slope gradients, an 8 m long, 2 m wide and 0.6 m deep slope adjustable soil pan was used to make an initial ephemeral gully channel on the soil bed according to the topographic characteristics of natural ephemeral gully after tillage and before rainy season. The initial ephemeral gully was placed at the middle of the soil bed with a depth of 12 cm. The soil used in this study was loess soil (fine-silty and mixed, with 28.3% sand, 58.1% silt, 13.6% clay), classified as Calcic Cambisols (USDA Soil Taxonomy). Soil materials were collected from 0 to 20 cm in the Ap horizon of a well-drained site in Ansai County, Shaanxi Province, and packed according to natural soil structure of the farmland on Loess Plateau. Two rainfall intensities (50 and 100 mm/h) and three typical slopes on which ephemeral gully occurred and developed (15°, 20°, and 25°) were designed in this experiment. Simulated rainfall and runoff scouring experiments were carried out at rainfall simulation laboratory of the State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau. Needle board method was used to measure the topography after 70 min experiment and the DEMs (digital elevation models) were generated in Surfer 10 with a resolution of 3 cm ×10 cm. The results showed that, the increasing of rainfall intensity and slope gradient accelerated the ephemeral gully erosion processes on the steep loessial hillslope. Average ephemeral gully width and depth in 25oand 100 mm/h condition were 1.40 and 0.61 times larger than those in 15° and 50 mm/h condition. Based on DEMs after rain and flow path figures, it could be concluded that the increasing of slope gradient increased the slope length required for the converging of rills and ephemeral gully channel, and the angle of rill and ephemeral gully channel at converging point was decreased; while the increasing of rainfall intensity decreased the slope length required for the converging of rills and ephemeral gully channel. Gully density, surface dissected degree and tortuosity complexity of ephemeral gully increased with the increase of rainfall intensity and slope gradient, varying from 0.74 to 1.48 m/m2, from 0.13 to 0.29, and from 1.64 to 2.84, respectively, while ephemeral gully channel width to depth ratio ranged from 0.65 to 1.27 and was the smallest when slope gradient was 20o. Directional derivative distribution was generated from original DEMs after rain according to the relationship between the neighbor grids, and it could be concluded that contour map of directional derivative grids reflected the length, surface area and gully bottom position of ephemeral gully and rills. More studies should be done on the ephemeral gully morphology and hillslope ephemeral gully erosion prediction model.

    soils; erosion; runoff; ephemeral gully; morphological characteristics; steep loessial hillslope; simulated experiment; directional derivatives

    10.11975/j.issn.1002-6819.2017.11.016

    S157.1

    A

    1002-6819(2017)-11-0124-09

    徐錫蒙,鄭粉莉,武 敏. 雨強和坡度對黃土陡坡地淺溝形態(tài)特征影響的定量研究[J]. 農(nóng)業(yè)工程學報,2017,33(11):124-132.

    10.11975/j.issn.1002-6819.2017.11.016 http://www.tcsae.org

    Xu Ximeng, Zheng Fenli, Wu Min. Quantification study of rainfall intensity and slope gradient impacts on ephemeral gully morphological characteristic on steep loessial hillslope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(11): 124-132. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.11.016 http://www.tcsae.org

    2016-11-18

    2017-05-02

    國家自然科學基金資助項目(41271299);水利部黃土高原水土流失過程與控制重點實驗室開放課題基金項目(2017001);黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室重要方向創(chuàng)新項目(A314021403-C2)

    徐錫蒙,男,河南衛(wèi)輝人,博士生,主要從事土壤侵蝕過程與機理研究。楊凌 西北農(nóng)林科技大學水土保持研究所,712100。

    Email:xuxm@nwsuaf.edu.cn

    ※通信作者:鄭粉莉,女,陜西藍田人,博士,教授,博士生導師,主要從事土壤侵蝕過程及預報和侵蝕環(huán)境效應評價研究。楊凌 西北農(nóng)林科技大學水土保持研究所,712100。Email:flzh@ms.iswc.ac.cn

    猜你喜歡
    淺溝細溝格網(wǎng)
    黃土高原淺溝空間分異特征
    黑土坡面細溝形態(tài)及剖面特征試驗研究
    實時電離層格網(wǎng)數(shù)據(jù)精度評估
    陜北子洲“7?26”暴雨后坡耕地細溝侵蝕及其影響因素分析
    南方紅壤丘陵區(qū)馬尾松人工林地淺溝形態(tài)特征
    黃土丘陵溝壑區(qū)淺溝發(fā)育動態(tài)監(jiān)測與形態(tài)定量研究
    細溝發(fā)育與形態(tài)特征研究進展
    基于QuickBird影像的黃土丘陵區(qū)坡面淺溝侵蝕速率研究
    基于空間信息格網(wǎng)與BP神經(jīng)網(wǎng)絡的災損快速評估系統(tǒng)
    尖山河小流域坡耕地的細溝侵蝕研究
    亚洲欧美一区二区三区国产| 国产熟女欧美一区二区| 老司机亚洲免费影院| 亚洲伊人久久精品综合| 永久免费av网站大全| 亚洲中文av在线| 欧美老熟妇乱子伦牲交| 亚洲美女搞黄在线观看| 99热全是精品| 老司机亚洲免费影院| 好男人视频免费观看在线| 爱豆传媒免费全集在线观看| 亚洲熟女精品中文字幕| 精品国产乱码久久久久久男人| 成年女人在线观看亚洲视频| 黑人猛操日本美女一级片| 看十八女毛片水多多多| 飞空精品影院首页| 在线观看免费高清a一片| 麻豆av在线久日| 午夜精品国产一区二区电影| 国产日韩一区二区三区精品不卡| 成人毛片a级毛片在线播放| 考比视频在线观看| www.自偷自拍.com| 卡戴珊不雅视频在线播放| 制服人妻中文乱码| 大码成人一级视频| 少妇的丰满在线观看| 日本wwww免费看| 国产黄频视频在线观看| 日韩人妻精品一区2区三区| 熟妇人妻不卡中文字幕| 国产一区二区 视频在线| 精品国产乱码久久久久久男人| 婷婷成人精品国产| 爱豆传媒免费全集在线观看| 精品人妻偷拍中文字幕| 久久久久国产网址| 精品国产一区二区三区四区第35| 国产精品蜜桃在线观看| 秋霞在线观看毛片| 国产精品久久久久成人av| 国产精品一区二区在线不卡| 久久精品国产自在天天线| 亚洲经典国产精华液单| 久久这里有精品视频免费| 国产在线一区二区三区精| 黑人巨大精品欧美一区二区蜜桃| 如何舔出高潮| 飞空精品影院首页| 另类精品久久| 一本—道久久a久久精品蜜桃钙片| 国产白丝娇喘喷水9色精品| 秋霞伦理黄片| 久久久a久久爽久久v久久| 国产无遮挡羞羞视频在线观看| 国产 精品1| 啦啦啦在线免费观看视频4| 精品一区在线观看国产| 久久影院123| 97在线视频观看| 亚洲欧洲日产国产| 一边亲一边摸免费视频| 亚洲av福利一区| 亚洲国产欧美在线一区| 日韩大片免费观看网站| 亚洲一级一片aⅴ在线观看| 亚洲国产欧美在线一区| 黄色一级大片看看| 一边摸一边做爽爽视频免费| 超碰成人久久| 99久久中文字幕三级久久日本| 亚洲经典国产精华液单| 亚洲第一av免费看| 亚洲图色成人| 日日摸夜夜添夜夜爱| 26uuu在线亚洲综合色| 久久免费观看电影| 久久精品人人爽人人爽视色| 秋霞伦理黄片| 另类亚洲欧美激情| 人妻 亚洲 视频| 亚洲精品久久成人aⅴ小说| 日日啪夜夜爽| 国产xxxxx性猛交| tube8黄色片| 最近2019中文字幕mv第一页| 色视频在线一区二区三区| 成人黄色视频免费在线看| 欧美+日韩+精品| 久久精品国产亚洲av天美| 九色亚洲精品在线播放| 精品国产国语对白av| 午夜免费鲁丝| 97在线视频观看| 日韩中文字幕视频在线看片| 精品一区二区免费观看| 国产一区二区在线观看av| 99久久中文字幕三级久久日本| 美女国产视频在线观看| 久久精品亚洲av国产电影网| 一级片免费观看大全| 亚洲精品aⅴ在线观看| 国产av码专区亚洲av| 亚洲欧美精品自产自拍| 美女脱内裤让男人舔精品视频| 久久毛片免费看一区二区三区| 丁香六月天网| 少妇 在线观看| 国产精品一二三区在线看| 激情视频va一区二区三区| 纯流量卡能插随身wifi吗| 午夜精品国产一区二区电影| 涩涩av久久男人的天堂| 久久精品久久久久久久性| 欧美国产精品va在线观看不卡| 亚洲成人手机| 一二三四中文在线观看免费高清| 综合色丁香网| 男男h啪啪无遮挡| 高清不卡的av网站| 啦啦啦中文免费视频观看日本| 黄片播放在线免费| 国产一区二区 视频在线| 两个人看的免费小视频| 两个人看的免费小视频| 91在线精品国自产拍蜜月| 黄色毛片三级朝国网站| 大片免费播放器 马上看| 99久久精品国产国产毛片| 波多野结衣一区麻豆| 看免费成人av毛片| 国产av一区二区精品久久| 亚洲国产欧美日韩在线播放| 久久青草综合色| 亚洲欧洲国产日韩| 在线观看免费日韩欧美大片| 国产精品一区二区在线观看99| 一个人免费看片子| 80岁老熟妇乱子伦牲交| 亚洲伊人色综图| 精品亚洲成a人片在线观看| h视频一区二区三区| 精品人妻在线不人妻| 97精品久久久久久久久久精品| 久久久久网色| 亚洲欧美精品综合一区二区三区 | 国产精品三级大全| 69精品国产乱码久久久| 国产av精品麻豆| 国产不卡av网站在线观看| 久久精品国产亚洲av天美| 亚洲国产成人一精品久久久| 侵犯人妻中文字幕一二三四区| 久久国内精品自在自线图片| 我要看黄色一级片免费的| 午夜91福利影院| 精品久久久精品久久久| 亚洲四区av| 97人妻天天添夜夜摸| 国产乱来视频区| 美女视频免费永久观看网站| 天天躁夜夜躁狠狠躁躁| 女的被弄到高潮叫床怎么办| 999久久久国产精品视频| 亚洲精品国产av蜜桃| 欧美成人精品欧美一级黄| 啦啦啦啦在线视频资源| 久久久久久久国产电影| 女的被弄到高潮叫床怎么办| 中文字幕制服av| 午夜影院在线不卡| av视频免费观看在线观看| 看十八女毛片水多多多| 永久免费av网站大全| 天堂中文最新版在线下载| 十八禁网站网址无遮挡| 岛国毛片在线播放| 亚洲精品美女久久av网站| 亚洲人成77777在线视频| 亚洲经典国产精华液单| 久久鲁丝午夜福利片| 久久女婷五月综合色啪小说| 高清黄色对白视频在线免费看| 国产精品国产三级专区第一集| 亚洲综合色惰| 亚洲国产欧美在线一区| 国产一区二区 视频在线| 国产精品久久久久久av不卡| 女性被躁到高潮视频| 丝袜在线中文字幕| 一区二区日韩欧美中文字幕| 国产黄色视频一区二区在线观看| 久久午夜福利片| 卡戴珊不雅视频在线播放| 国产成人aa在线观看| 婷婷成人精品国产| videossex国产| 国产精品免费大片| 18禁裸乳无遮挡动漫免费视频| 免费观看a级毛片全部| 亚洲一码二码三码区别大吗| 性高湖久久久久久久久免费观看| 亚洲精品美女久久久久99蜜臀 | 男人操女人黄网站| 天天操日日干夜夜撸| 老鸭窝网址在线观看| 精品国产露脸久久av麻豆| 国产精品 国内视频| 精品国产乱码久久久久久小说| 午夜激情久久久久久久| 啦啦啦中文免费视频观看日本| 日韩成人av中文字幕在线观看| 亚洲欧洲日产国产| 日韩大片免费观看网站| 色婷婷av一区二区三区视频| 中文字幕人妻熟女乱码| 亚洲精品乱久久久久久| 亚洲av福利一区| 男女国产视频网站| 日韩成人av中文字幕在线观看| 丝袜美足系列| 成年人午夜在线观看视频| 国产精品免费视频内射| 天天躁夜夜躁狠狠躁躁| 秋霞伦理黄片| 五月开心婷婷网| 老鸭窝网址在线观看| 26uuu在线亚洲综合色| 黑人猛操日本美女一级片| 亚洲精品久久午夜乱码| av电影中文网址| 国产无遮挡羞羞视频在线观看| 十八禁高潮呻吟视频| 一区二区日韩欧美中文字幕| 国产不卡av网站在线观看| 香蕉丝袜av| 国产成人精品无人区| 五月伊人婷婷丁香| 国产亚洲欧美精品永久| 亚洲,一卡二卡三卡| 777久久人妻少妇嫩草av网站| 91成人精品电影| 熟妇人妻不卡中文字幕| 久久久欧美国产精品| 欧美变态另类bdsm刘玥| 看免费成人av毛片| 男女啪啪激烈高潮av片| 日本wwww免费看| 国产成人av激情在线播放| 性色avwww在线观看| 日韩一区二区三区影片| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av涩爱| 国产女主播在线喷水免费视频网站| 精品国产一区二区久久| 在线观看国产h片| 亚洲一区中文字幕在线| 啦啦啦啦在线视频资源| h视频一区二区三区| 丝袜美腿诱惑在线| 日本vs欧美在线观看视频| 午夜福利网站1000一区二区三区| 日本黄色日本黄色录像| 99国产精品免费福利视频| 亚洲欧美成人精品一区二区| 巨乳人妻的诱惑在线观看| 在线观看人妻少妇| 国产欧美日韩一区二区三区在线| 女的被弄到高潮叫床怎么办| 久久久久久免费高清国产稀缺| 久久精品夜色国产| 婷婷成人精品国产| 精品亚洲成国产av| 青春草视频在线免费观看| 欧美老熟妇乱子伦牲交| 丝袜美腿诱惑在线| 国产精品一国产av| 亚洲欧洲精品一区二区精品久久久 | 爱豆传媒免费全集在线观看| 99久久中文字幕三级久久日本| 国产精品秋霞免费鲁丝片| 性少妇av在线| 亚洲伊人久久精品综合| 久久久精品国产亚洲av高清涩受| 中文精品一卡2卡3卡4更新| 男女高潮啪啪啪动态图| 我的亚洲天堂| 国产免费一区二区三区四区乱码| 欧美日韩av久久| 黄色 视频免费看| 亚洲欧美色中文字幕在线| 欧美人与善性xxx| 欧美日韩一区二区视频在线观看视频在线| 国产黄频视频在线观看| 中国三级夫妇交换| 人人妻人人添人人爽欧美一区卜| av不卡在线播放| 香蕉精品网在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美精品高潮呻吟av久久| 少妇被粗大的猛进出69影院| 亚洲男人天堂网一区| 久久精品国产亚洲av高清一级| 男人操女人黄网站| 亚洲欧洲日产国产| 亚洲 欧美一区二区三区| 亚洲五月色婷婷综合| 国产免费又黄又爽又色| 国产在线免费精品| 两性夫妻黄色片| 国产亚洲最大av| 人妻人人澡人人爽人人| 亚洲精品国产av成人精品| 一区二区三区乱码不卡18| 国产无遮挡羞羞视频在线观看| 美女中出高潮动态图| 熟女电影av网| 日韩制服骚丝袜av| 综合色丁香网| 侵犯人妻中文字幕一二三四区| 精品人妻在线不人妻| 免费观看av网站的网址| 人人妻人人澡人人看| 久久国产亚洲av麻豆专区| 国产亚洲一区二区精品| 亚洲,欧美精品.| 久久99精品国语久久久| 黑人猛操日本美女一级片| 最近的中文字幕免费完整| 亚洲精品久久成人aⅴ小说| 日本-黄色视频高清免费观看| 中文天堂在线官网| 在线观看三级黄色| 欧美精品av麻豆av| 交换朋友夫妻互换小说| 精品福利永久在线观看| 中文精品一卡2卡3卡4更新| 久久精品久久精品一区二区三区| 亚洲欧美日韩另类电影网站| 欧美成人精品欧美一级黄| 老汉色av国产亚洲站长工具| 亚洲欧美精品自产自拍| 亚洲av.av天堂| 国产免费视频播放在线视频| 久久久久国产精品人妻一区二区| 国产xxxxx性猛交| 男女啪啪激烈高潮av片| 成人漫画全彩无遮挡| 人妻系列 视频| 日本黄色日本黄色录像| av福利片在线| 一二三四中文在线观看免费高清| 女的被弄到高潮叫床怎么办| 波野结衣二区三区在线| 99热国产这里只有精品6| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美亚洲二区| 国产在线视频一区二区| 18在线观看网站| 高清av免费在线| 欧美 亚洲 国产 日韩一| 久久久久久久亚洲中文字幕| 国产av精品麻豆| 美女国产视频在线观看| 午夜福利影视在线免费观看| 亚洲国产毛片av蜜桃av| 国产成人精品久久二区二区91 | 丝袜美腿诱惑在线| 成人国语在线视频| 春色校园在线视频观看| 亚洲国产毛片av蜜桃av| 日本-黄色视频高清免费观看| 免费久久久久久久精品成人欧美视频| 超碰成人久久| 丁香六月天网| 最近最新中文字幕大全免费视频 | 中文字幕亚洲精品专区| 看免费成人av毛片| 香蕉国产在线看| 久久 成人 亚洲| 99久久人妻综合| 亚洲精品aⅴ在线观看| 伦理电影免费视频| 亚洲av综合色区一区| 欧美激情 高清一区二区三区| 日韩一本色道免费dvd| 伦理电影免费视频| 黄网站色视频无遮挡免费观看| 一区福利在线观看| 街头女战士在线观看网站| 亚洲,一卡二卡三卡| 久久久久久久亚洲中文字幕| 精品国产乱码久久久久久男人| 一本—道久久a久久精品蜜桃钙片| 午夜日韩欧美国产| 巨乳人妻的诱惑在线观看| 电影成人av| 久久精品久久精品一区二区三区| 免费黄频网站在线观看国产| 亚洲精品aⅴ在线观看| 国产一区二区 视频在线| 久久精品国产亚洲av天美| 高清黄色对白视频在线免费看| av一本久久久久| 欧美 日韩 精品 国产| 日本免费在线观看一区| av一本久久久久| 久热久热在线精品观看| 久久久精品免费免费高清| 2022亚洲国产成人精品| 国产男人的电影天堂91| 欧美人与性动交α欧美软件| 天天影视国产精品| 在线天堂最新版资源| 亚洲av.av天堂| 99久久精品国产国产毛片| 伦精品一区二区三区| 免费播放大片免费观看视频在线观看| 国产av国产精品国产| 欧美成人午夜精品| 国产不卡av网站在线观看| 欧美日韩成人在线一区二区| 伦理电影免费视频| 国产成人欧美| 日韩,欧美,国产一区二区三区| 精品国产一区二区久久| 免费在线观看视频国产中文字幕亚洲 | 搡老乐熟女国产| 国产欧美日韩一区二区三区在线| 日本-黄色视频高清免费观看| 国产综合精华液| 国产一级毛片在线| 日本免费在线观看一区| 国产人伦9x9x在线观看 | 欧美亚洲 丝袜 人妻 在线| 国产1区2区3区精品| 久久狼人影院| 成人黄色视频免费在线看| 欧美 亚洲 国产 日韩一| 亚洲人成网站在线观看播放| 亚洲成色77777| 一区二区三区精品91| 国产成人免费无遮挡视频| 青春草视频在线免费观看| 只有这里有精品99| 亚洲精品国产av成人精品| 国产成人午夜福利电影在线观看| 精品国产乱码久久久久久小说| av不卡在线播放| 国产午夜精品一二区理论片| 免费在线观看黄色视频的| 亚洲欧美一区二区三区国产| 岛国毛片在线播放| www.精华液| 亚洲av在线观看美女高潮| 最新的欧美精品一区二区| 夜夜骑夜夜射夜夜干| 国产探花极品一区二区| 亚洲综合色网址| 国产爽快片一区二区三区| 久久久久久久久久久久大奶| 91精品国产国语对白视频| 国产精品国产av在线观看| 日韩制服骚丝袜av| 韩国av在线不卡| 久久久国产精品麻豆| 久久久精品94久久精品| 成人二区视频| av网站免费在线观看视频| 热re99久久国产66热| kizo精华| 最近最新中文字幕免费大全7| 激情视频va一区二区三区| 亚洲综合精品二区| 人妻一区二区av| 免费久久久久久久精品成人欧美视频| 九九爱精品视频在线观看| 亚洲av欧美aⅴ国产| 久久精品国产综合久久久| 欧美黄色片欧美黄色片| 爱豆传媒免费全集在线观看| 久热久热在线精品观看| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 亚洲第一青青草原| 亚洲综合色网址| 国产在线一区二区三区精| 国产极品天堂在线| 五月天丁香电影| 人体艺术视频欧美日本| 亚洲精品久久午夜乱码| 国产亚洲精品第一综合不卡| 一级,二级,三级黄色视频| 久久 成人 亚洲| 欧美97在线视频| 亚洲国产av影院在线观看| 下体分泌物呈黄色| 欧美日本中文国产一区发布| 国产成人精品婷婷| 妹子高潮喷水视频| 天天操日日干夜夜撸| 免费观看性生交大片5| 免费黄频网站在线观看国产| 免费在线观看黄色视频的| av不卡在线播放| 欧美+日韩+精品| 免费不卡的大黄色大毛片视频在线观看| 精品少妇一区二区三区视频日本电影 | 一级毛片 在线播放| 女人被躁到高潮嗷嗷叫费观| 天天躁夜夜躁狠狠久久av| 2021少妇久久久久久久久久久| 亚洲成色77777| 国产精品女同一区二区软件| 日韩欧美一区视频在线观看| 国产精品人妻久久久影院| 国产在线一区二区三区精| 亚洲天堂av无毛| 三级国产精品片| 国产精品 国内视频| 欧美 亚洲 国产 日韩一| 亚洲av电影在线进入| 日日爽夜夜爽网站| 黄片小视频在线播放| 国产片内射在线| 欧美在线黄色| av又黄又爽大尺度在线免费看| 久久精品国产a三级三级三级| 大香蕉久久成人网| 丝袜脚勾引网站| 免费观看av网站的网址| 少妇猛男粗大的猛烈进出视频| 亚洲精华国产精华液的使用体验| 啦啦啦在线观看免费高清www| 久久精品国产a三级三级三级| 男人操女人黄网站| 国产精品免费视频内射| 久久99一区二区三区| 婷婷色综合www| 久久女婷五月综合色啪小说| 色视频在线一区二区三区| 国产av一区二区精品久久| 国产一区亚洲一区在线观看| 极品少妇高潮喷水抽搐| 搡女人真爽免费视频火全软件| 久久久久久久大尺度免费视频| 国产精品无大码| 青春草国产在线视频| 日韩av在线免费看完整版不卡| 亚洲国产精品成人久久小说| 永久网站在线| 午夜91福利影院| 国产男女超爽视频在线观看| 黄色毛片三级朝国网站| 亚洲久久久国产精品| av一本久久久久| 亚洲精品久久成人aⅴ小说| 欧美精品国产亚洲| 久久久a久久爽久久v久久| 国产一区亚洲一区在线观看| 久久狼人影院| 熟女少妇亚洲综合色aaa.| 精品人妻在线不人妻| 有码 亚洲区| 久久国产精品男人的天堂亚洲| 26uuu在线亚洲综合色| 如何舔出高潮| 国产精品无大码| 精品国产超薄肉色丝袜足j| 亚洲一区二区三区欧美精品| 黄色视频在线播放观看不卡| 交换朋友夫妻互换小说| 国产精品久久久久久精品电影小说| 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| 寂寞人妻少妇视频99o| 久久国内精品自在自线图片| 免费观看av网站的网址| 欧美 亚洲 国产 日韩一| 亚洲第一区二区三区不卡| 免费在线观看视频国产中文字幕亚洲 | 一边亲一边摸免费视频| 亚洲国产精品国产精品| 99热国产这里只有精品6| 丝袜脚勾引网站| 午夜福利在线免费观看网站| 热re99久久国产66热| 免费高清在线观看日韩| 性色avwww在线观看| 在线天堂最新版资源| videossex国产| 大陆偷拍与自拍| 日韩免费高清中文字幕av| 国产有黄有色有爽视频| 欧美人与性动交α欧美软件| 免费观看在线日韩| 国产精品嫩草影院av在线观看| 国产成人精品久久二区二区91 | 亚洲精品乱久久久久久| 亚洲第一av免费看| 国产精品三级大全| 宅男免费午夜| 日韩成人av中文字幕在线观看| 免费大片黄手机在线观看| 午夜福利视频精品| 最近中文字幕2019免费版| 免费在线观看视频国产中文字幕亚洲 | 少妇熟女欧美另类| 这个男人来自地球电影免费观看 | 午夜福利影视在线免费观看| 精品国产一区二区三区久久久樱花| 亚洲国产最新在线播放| 亚洲av电影在线进入|