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

    生物結(jié)皮坡面不同降雨歷時的產(chǎn)流特征

    2019-02-21 00:32:10趙允格張子輝谷康民郭雅麗王閃閃
    農(nóng)業(yè)工程學(xué)報 2019年23期
    關(guān)鍵詞:結(jié)皮土坡產(chǎn)流

    楊 凱,趙 軍,趙允格,,張子輝,孫 會,谷康民,郭雅麗,王閃閃

    生物結(jié)皮坡面不同降雨歷時的產(chǎn)流特征

    楊 凱1,2,趙 軍2,趙允格2,4※,張子輝1,孫 會3,谷康民1,郭雅麗1,王閃閃4

    (1. 西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院,楊凌 712100; 2. 西北農(nóng)林科技大學(xué)水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室,楊凌 712100; 3. 西北農(nóng)林科技大學(xué)林學(xué)院,楊凌 712100;4. 中國科學(xué)院大學(xué),北京 100049)

    黃土高原退耕還林(草)工程實施后,生物結(jié)皮廣泛發(fā)育,顯著影響坡面產(chǎn)流。已有大量研究探索了生物結(jié)皮對徑流的影響,但相關(guān)結(jié)論存在較大分歧。該研究以黃土高原典型生物結(jié)皮坡面為研究對象,通過人工模擬降雨試驗,研究了生物結(jié)皮坡面產(chǎn)流過程。結(jié)果表明:生物結(jié)皮坡面較翻耕后的裸土坡面顯著降低了初始產(chǎn)流時間,裸土坡面初始產(chǎn)流時間是生物結(jié)皮坡面的1.59~3.04倍。生物結(jié)皮蓋度與初始產(chǎn)流時間之間呈顯著的負相關(guān)關(guān)系;降雨15 min時與60 min時生物結(jié)皮對坡面徑流的影響發(fā)生逆轉(zhuǎn),90 mm/h的雨強下,當降雨歷時為15 min時,生物結(jié)皮坡面較裸土坡面增加75.42%的徑流;當降雨歷時為60 min時,生物結(jié)皮坡面徑流量較裸土坡面降低52.42%;生物結(jié)皮影響了土壤水分入滲速率,導(dǎo)致生物結(jié)皮坡面與裸土坡面隨降雨歷時變化的產(chǎn)流特征出現(xiàn)差異,裸土坡面降雨60 min時的入滲率較15 min時降低了34.30%,高蓋度生物結(jié)皮坡面降低了6.38%;生物結(jié)皮對坡面入滲產(chǎn)流的影響與降雨歷時有極大的關(guān)系,降雨歷時不同,很可能得到截然相反的結(jié)論,考慮生物結(jié)皮因素的野外降雨試驗,降雨歷時應(yīng)不少于45 min。研究結(jié)果為解釋生物結(jié)皮影響坡面入滲產(chǎn)流方面存在的分歧提供了科學(xué)依據(jù),進一步明確了干旱半干旱地區(qū)生物結(jié)皮的水文效應(yīng)。

    徑流;水動力學(xué);生物土壤結(jié)皮;初始產(chǎn)流時間;徑流系數(shù)

    0 引 言

    生物土壤結(jié)皮(生物結(jié)皮)是由生長于土壤表層和近地表數(shù)毫米內(nèi)的藍綠藻、苔蘚、地衣、真菌,以及許多景觀中常見的其它非維管束植物成分與土壤復(fù)合而形成的具有生命活性的復(fù)雜復(fù)合體,是干旱半干旱地區(qū)普遍存在的生物地被物,對土壤水分、養(yǎng)分、理化屬性及抗侵蝕性等具有重要的影響[1]。生物結(jié)皮由于微生物分泌物粘結(jié)、假根系/菌絲體的捆綁束縛作用,會形成水平方向極其穩(wěn)定的層狀結(jié)構(gòu)[2-3]。同時,生物結(jié)皮的發(fā)育顯著改變了土壤表面特性,如粗糙度、持水性、斥水性等[4-5](王媛, 2014 #82)。表層土壤在水分循環(huán)中起著重要的作用,如降水入滲、地表徑流等過程都是以表層土壤為介質(zhì)發(fā)生和轉(zhuǎn)化的,因此生物結(jié)皮的發(fā)育影響了土壤水分的入滲和產(chǎn)流。就此,國內(nèi)已進行了大量的研究,但目前所得結(jié)論存在較大分歧。一些研究認為生物結(jié)皮延長了水分在地表的停滯時間[6],從而促進了入滲。Galun等[7]在以色列沙漠地區(qū)的研究發(fā)現(xiàn),藍藻地衣幾分鐘內(nèi)即可吸收其本身干重或體積的3~13倍的水分,進而降低產(chǎn)流。澳大利亞的研究認為生物結(jié)皮改善了下層土壤結(jié)構(gòu),增加了水分入滲,降低了坡面產(chǎn)流[8-9]。肖波等[10]使用砂質(zhì)壤土的裝填土槽,接種生物結(jié)皮后通過1 h的人工降雨試驗發(fā)現(xiàn),生物結(jié)皮坡面較無生物結(jié)皮坡面減少了49%~64%的徑流量。

    另有研究認為生物結(jié)皮堵塞了地表土壤孔隙,導(dǎo)致入滲降低和徑流增加[11]。Rodríguez-Caballero[12]在西班牙南部研究了自然降雨下坡面尺度生物結(jié)皮對徑流的影響,結(jié)果表明生物結(jié)皮的存在增加了坡面徑流量。李新榮等[13]在沙坡頭人工植被固沙區(qū)比較了自然降水后生物結(jié)皮的濕潤峰,結(jié)果表明,生物結(jié)皮降低了水分入滲率。Zhao等[14]在黃土高原地區(qū)采用15 min放水試驗研究藻結(jié)皮和蘚結(jié)皮對坡面產(chǎn)流的影響,發(fā)現(xiàn)生物結(jié)皮較無生物結(jié)皮坡面增加了10%~15%的徑流量。

    除此,還有研究認為生物結(jié)皮對入滲無明顯作用。Booth[15]使用高強度耐壓水管噴灑藻結(jié)皮和無藻結(jié)皮的樣方,試驗結(jié)果表明,2類樣方上的入滲率無明顯差異。Eldridge等[16]認為在生物結(jié)皮發(fā)育良好的地區(qū),下層土壤的大空隙對水分入滲的影響占絕對優(yōu)勢,生物結(jié)皮的影響微不足道。而在生物結(jié)皮缺失的地區(qū),地表土壤缺乏大孔隙,水分入滲率本來就低。因此,生物結(jié)皮對入滲的影響應(yīng)歸因于土壤物理性質(zhì)、水分進入土壤通道的不同以及地表侵蝕史的差異,而非生物結(jié)皮。綜上所述,國內(nèi)外專家學(xué)者關(guān)于生物結(jié)皮對水分入滲及徑流的影響已經(jīng)取得不少成果,但結(jié)論仍存在較大分歧。因此,有關(guān)生物結(jié)皮對坡面產(chǎn)流的影響,仍需要進一步研究。

    黃土高原是中國乃至全球土壤侵蝕最為嚴重的區(qū)域之一。在20世紀末,由于退耕還林(草)工程的實施,生物結(jié)皮在退耕坡地上廣泛發(fā)育,平均蓋度達到60%~70%,甚至更高[17],是退耕坡地表面比例最大的覆蓋因子[18],顯著影響該區(qū)坡面的水土流失。已有的研究證實,生物結(jié)皮顯著影響坡面入滲及產(chǎn)流過程[19-20]。在黃土丘陵區(qū),暴雨多集中在6—9月,且降雨歷時不固定,亦有可能影響到生物結(jié)皮坡面的產(chǎn)流特征。因此,本文借助野外人工模擬降雨手段,研究了生物結(jié)皮坡面隨降雨歷時的入滲產(chǎn)流特征,以揭示降雨歷時的差異對生物結(jié)皮坡面入滲產(chǎn)流的影響,為黃土高原的生物結(jié)皮坡面水文過程研究及水土保持工作提供理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    試驗于2018年8月在陜西省榆林市定邊縣楊井鎮(zhèn)的退耕撂荒地進行。該區(qū)屬于陜北黃土高原北部與毛烏素沙漠南緣的過渡地帶,屬溫帶半干旱內(nèi)陸性氣候,多年平均氣溫為8.7 ℃,多年平均降水量316.9 mm,年蒸發(fā)量2 490 mm,地形屬于黃土高原的丘陵溝壑區(qū),土壤以石灰性黃綿土為主。試驗樣地主要植被為達烏里胡枝子()、早熟禾()、百里香()、茵陳蒿()、賴草()、長芒草()等。樣地生物結(jié)皮平均蓋度79.2%,主要由蘚結(jié)皮、藻結(jié)皮及二者混合構(gòu)成。其中,蘚結(jié)皮優(yōu)勢種主要有短葉對齒蘚()、土生對齒蘚()及銀葉真蘚()等。

    1.2 試驗設(shè)計

    樣地選取及小區(qū)處理:選取蘚結(jié)皮及蘚藻混合結(jié)皮為主的撂荒坡地作為研究樣地,樣地基本情況見表1。

    表1 樣地概況

    樣地面積約為80 m×20 m,坡度15°左右。用薄鋼板在該樣地圈建10 m×2.1 m的試驗小區(qū),用剪刀去除小區(qū)內(nèi)的高等植物冠層。根據(jù)生物結(jié)皮總蓋度差異設(shè)置2個坡面處理,將無干擾的結(jié)皮坡面作為高蓋度生物結(jié)皮組(平均蓋度79.2%),用鏟子移除一部分結(jié)皮后的結(jié)皮坡面作為低蓋度生物結(jié)皮組(平均蓋度43.6%),同時以翻耕平整過的裸土坡面作為對照,每個試驗處理設(shè)3個重復(fù),共計9個徑流小區(qū)。

    模擬降雨試驗:采用槽式擺噴頭下噴式人工模擬降雨系統(tǒng)(EL-RS3/5,北京易科立德生態(tài)環(huán)境科技有限責(zé)任公司)。降雨強度設(shè)為90 mm/h,降雨歷時1 h。

    1.3 試驗步驟

    降雨開始前率定雨強,保證降雨的均勻性與精確度[21]。降雨開始后,當坡面整體開始產(chǎn)流,且出水口形成連續(xù)水流時,記下此時的產(chǎn)流時間并接取第一個泥沙徑流樣,此后每隔3 min收集徑流泥沙樣。降雨過程中采用高錳酸鉀測取徑流流過1 m坡段的時間,計算流速。泥沙測定采用烘干法。

    1.4 數(shù)據(jù)分析

    產(chǎn)流率、入滲速率、徑流系數(shù)、佛汝德數(shù)(值)及水深計算公式如下

    式中為產(chǎn)流率,mm/min;R為第次取樣的產(chǎn)流量,L;為入滲速率,mm/min;為降雨強度,mm/min;為徑流系數(shù);為產(chǎn)流時間,min;為坡面面積,m2;為坡度,(°);為佛汝德數(shù);為坡面平均流速,m/s;為重力加速度,m/s2;為水深,m;為單寬流量,m3/(m·min)。

    在小區(qū)上中下部各選取5個點布設(shè)25 cm×25 cm樣方,記錄樣方中藻、蘚、裸土、植物根基、枯落物出現(xiàn)的頻次,以各類物種占調(diào)查總點數(shù)的百分數(shù)作為生物結(jié)皮蓋度。用針狀糙度計測定粗糙度。糙度計由52根金屬針組成,針隨測定地面的凸凹狀況而自由升降,針上部各點相對與參照基準面的高度變化反映地面的起伏程度。用數(shù)碼照相機將針的起伏狀況拍攝下,用Profile meter程序處理照片,計算每根針相對于參照基準面的高度。粗糙度用各測點高度的標準差表示[22]。每個試驗小區(qū)分別測定平行于等高線和垂直于等高線方向的地面粗糙度各9次,取平均值。有機質(zhì)含量采用重鉻酸鉀外加熱法測定。顆粒組成采用馬爾文MS2000激光粒度分析儀進行測定。

    本文數(shù)據(jù)采用excel2010和spss19.0進行處理分析,分析前對數(shù)據(jù)進行正態(tài)分布檢測和方差齊性檢驗,對初始產(chǎn)流時間、產(chǎn)流率、入滲速率進行單因素方差分析、LSD多重比較。對初始產(chǎn)流時間、產(chǎn)流率、以及生物結(jié)皮蓋度、粗糙度進行Pearson雙尾相關(guān)性分析。

    2 結(jié)果與分析

    2.1 不同處理對坡面初始產(chǎn)流時間的影響

    圖1為90 mm/h降雨條件下3個處理坡面產(chǎn)流量的時間動態(tài)。隨著生物結(jié)皮蓋度的降低,初始產(chǎn)流時間增加,平均初始產(chǎn)流時間由大到小依次為裸土7.0 min、低蓋度生物結(jié)皮4.4 min、高蓋度生物結(jié)皮2.3 min,3種處理之間差異顯著(<0.05)。裸土坡面初始產(chǎn)流時間是生物結(jié)皮坡面的1.59~3.04倍。初始產(chǎn)流時間與生物結(jié)皮蓋度間的Pearson雙尾相關(guān)分析表明,初始產(chǎn)流時間與生物結(jié)皮蓋度之間存在極顯著的負相關(guān)關(guān)系(=?0.925,<0.01)。

    注:不同小寫字母表示差異顯著(P<0.05),下同。

    2.2 坡面產(chǎn)流隨降雨歷時變化過程

    不同處理坡面的瞬時產(chǎn)流量隨降雨歷時變化過程如圖2所示??梢娐阃疗旅娉跏籍a(chǎn)流時間高于不同蓋度生物結(jié)皮坡面,產(chǎn)流后徑流量快速增加。高蓋度生物結(jié)皮坡面產(chǎn)流早于裸土和低蓋度生物結(jié)皮坡面,降雨15 min后趨于平穩(wěn)。低蓋度生物結(jié)皮坡面徑流變化處于二者之間,降雨45 min后趨于穩(wěn)定。

    圖2 不同處理坡面產(chǎn)流量隨降雨歷時的變化

    各處理坡面不同降雨歷時的產(chǎn)流率如圖3所示。各處理坡面降雨歷時15 min時產(chǎn)流率由高至低依次為高蓋度生物結(jié)皮0.20 mm/min、裸土0.13 mm/min、低蓋度生物結(jié)皮0.10 mm/min,高蓋度生物結(jié)皮坡面產(chǎn)流率顯著高于低蓋度生物結(jié)皮(<0.05)。降雨歷時30 min時各處理坡面產(chǎn)流率差異不顯著(>0.05)。降雨歷時為45 min時裸土坡面產(chǎn)流率顯著高于高蓋度生物結(jié)皮(<0.05)。降雨歷時60 min時裸土坡面產(chǎn)流率達到最大,為0.56 mm/min,顯著高于高蓋度生物結(jié)皮坡面(0.28 mm/min)。生物結(jié)皮坡面與裸土坡面在不同降雨時段的產(chǎn)流量有逆轉(zhuǎn)表現(xiàn),當降雨歷時為15 min時,高蓋度生物結(jié)皮坡面較裸土增加75.42%的徑流;當降雨歷時為60 min時,高蓋度生物結(jié)皮坡面徑流較裸土降低52.42%。

    圖3 不同降雨歷時各坡面產(chǎn)流率

    為解析降雨歷時的影響,對15、30、45、60 min時的產(chǎn)流率與生物結(jié)皮蓋度、容重及粗糙度進行了Pearson雙尾相關(guān)分析(表2)。結(jié)果表明,15、30 min產(chǎn)流率與粗糙度、容重及生物結(jié)皮蓋度均無相關(guān)關(guān)系;而隨著降雨時間的增加,45 min產(chǎn)流率與這3項指標顯著相關(guān),60 min降雨率與容重顯著相關(guān),與粗糙度、生物結(jié)皮蓋度極顯著相關(guān)。相關(guān)系數(shù)絕對值隨著降雨歷時的增加而變大。

    表2 不同降雨歷時產(chǎn)流率與影響因素相關(guān)性

    注:*和**分別表示在0.05水平和0.01水平(雙側(cè))上顯著相關(guān)。

    Note: * and ** indicated significant correlation at< 0.05 and<0.01 (two-tail), respectively.

    2.3 徑流系數(shù)及入滲速率隨降雨歷時的變化特征

    圖4為不同處理坡面徑流系數(shù)隨降雨歷時的變化過程。各處理坡面在產(chǎn)流初始時的徑流系數(shù)由大到小依次為高蓋度生物結(jié)皮0.08,裸土0.07,低蓋度生物結(jié)皮0.04。生物結(jié)皮坡面徑流系數(shù)隨降雨歷時的變化趨勢均為先增加,最終趨于穩(wěn)定。裸土坡面徑流系數(shù)變化趨勢表現(xiàn)為先增加,在降雨歷時為60 min時有下降趨勢。高蓋度生物結(jié)皮坡面徑流系數(shù)隨降雨歷時變化幅度最??;裸土坡面的徑流系數(shù)在產(chǎn)流后隨降雨歷時增加迅速增大,在后續(xù)的降雨過程中始終高于結(jié)皮坡面;低蓋度結(jié)皮坡面徑流系數(shù)介于在二者之間。

    圖4 不同處理坡面徑流系數(shù)隨降雨歷時的變化

    圖5為各處理坡面不同降雨歷時的入滲速率。裸土坡面的入滲速率在前30 min時與生物結(jié)皮坡面差異不顯著。45 min時的入滲速率由大到小依次為高蓋度生物結(jié)皮1.18 mm/min、低蓋度生物結(jié)皮1.13 mm/min、裸土0.98 mm/min,高蓋度生物結(jié)皮坡面入滲速率顯著大于裸土坡面(<0.05)。裸土坡面入滲速率在60 min時降低至0.89 mm/min,顯著低于高蓋度生物結(jié)皮坡面(1.17 mm/min,<0.05)。

    圖5 不同降雨歷時各處理坡面入滲速率變化

    2.4 坡面流速與流態(tài)隨降雨歷時的變化特征

    不同處理降雨過程中流速的時間動態(tài)見圖6a。其中裸土坡面的平均流速為0.23 m/s,顯著高于有結(jié)皮覆蓋的坡面(<0.05),隨著時間的增加呈波浪狀增長趨勢,在降雨末期(50 min后)有下降趨勢。高蓋度生物結(jié)皮覆蓋的坡面平均流速最低,為0.05 m/s,沒有明顯的變化趨勢。低蓋度生物結(jié)皮坡面流速為0.06 m/s,略大于高蓋度生物結(jié)皮坡面,在降雨中期流速波動幅度較大。

    佛汝德數(shù)(值)是流體內(nèi)慣性力與重力的比值,是判別水流狀態(tài)的重要參數(shù)。研究表明當值>1時,水流為急流,值<1時,水流為緩流[23]。由圖6b可見,3種坡面處理的流態(tài)變化特征與流速相似。裸土坡面的Fr值顯著高于生物結(jié)皮坡面,最大值達到了5.09,平均值為3.41,是生物結(jié)皮坡面的6.8~7.6倍。表明生物結(jié)皮在坡面上的覆蓋改變了降雨過程中坡面徑流的流態(tài)。

    3 討 論

    不同降雨歷時的產(chǎn)流率及生物結(jié)皮蓋度Pearson相關(guān)分析結(jié)果表明,產(chǎn)流率在降雨前30 min的與生物結(jié)皮蓋度無顯著關(guān)系,隨著降雨歷時的增加,45、60 min的產(chǎn)流率與生物結(jié)皮蓋度顯著相關(guān),且相關(guān)系數(shù)絕對值隨降雨歷時增加而變大。表明隨著降雨歷時的增加,生物結(jié)皮對坡面產(chǎn)流的影響也隨之增大。這可能是由于裸土坡面和生物結(jié)皮坡面的結(jié)構(gòu)隨降雨歷時變化的差異造成的。裸土坡面土壤結(jié)構(gòu)未破壞時孔隙度較大,水分入滲在降雨前期為大孔隙流,而生物結(jié)皮的飽和導(dǎo)水率低于無結(jié)皮土壤[24],此時裸土坡面產(chǎn)流量低于結(jié)皮坡面。隨著降雨歷時的增加,裸土坡面的土壤結(jié)構(gòu)容易被雨滴擊濺破壞,濺蝕產(chǎn)生的細土粒堵塞土壤孔隙,形成大面積的物理結(jié)皮,其致密的表層結(jié)構(gòu)會顯著降低表層土壤的入滲速率,導(dǎo)致坡面的產(chǎn)流量迅速增大[25]。已有研究表明,生物結(jié)皮的存在可以削減雨滴動能,減少雨滴對下層土壤結(jié)構(gòu)的破壞,增加土壤水穩(wěn)定性[3,26]。因此,生物結(jié)皮坡面的表層土壤在降雨過程中較裸土坡面更加穩(wěn)定,水分入滲率隨降雨歷時變化不大?;谝陨显?,生物結(jié)皮坡面的產(chǎn)流率隨降雨歷時變化平穩(wěn),裸土坡面由于表層土壤結(jié)構(gòu)的不穩(wěn)定性導(dǎo)致產(chǎn)流隨降雨歷時增加而迅速增大。

    本研究團隊在該地區(qū)采用線源入流入滲法測定了高蓋度生物結(jié)皮小區(qū)與裸土小區(qū)的水分入滲時間動態(tài)。分析發(fā)現(xiàn),裸土小區(qū)前6 min的初始入滲速率高于高蓋度生物結(jié)皮小區(qū),隨著入滲時間的增加,裸土小區(qū)的入滲速率迅速下降,在7~20 min時段低于高蓋度生物結(jié)皮小區(qū),隨后緩慢下降至穩(wěn)定數(shù)值。高蓋度結(jié)皮小區(qū)初始入滲速率較穩(wěn)定,8 min后開始緩慢下降,最后下降至穩(wěn)定的入滲速率。該結(jié)果證明了生物結(jié)皮影響了水分入滲的過程,從入滲的角度解釋了裸土坡面與生物結(jié)皮坡面在不同降雨歷時表現(xiàn)出產(chǎn)流特征差異的原因。

    生物結(jié)皮的覆蓋顯著地降低了暴雨降雨條件下坡面的流速(圖6a),從而改變了坡面徑流的佛汝德數(shù)(值),裸土坡面的值顯著高于結(jié)皮坡面,平均值是生物結(jié)皮坡面的6.8~7.6倍。Savat和Ploey[27]的試驗表明,當坡面徑流的佛汝德數(shù)為2~3時,細溝發(fā)育的概率很大。張科利等[28]認為佛汝德數(shù)≥1是細溝發(fā)生的水動力臨界值。這意味著生物結(jié)皮的存在能夠?qū)⑵旅嫠鞯牧鲬B(tài)從急流變?yōu)榫徚?,減少侵蝕發(fā)生的可能性。本試驗中裸土坡面的值隨降雨歷時增加呈波浪式上升的趨勢,而生物結(jié)皮坡面的變化趨勢則相對平穩(wěn),這表明裸土坡面表層土壤狀態(tài)在降雨過程中一直不斷變化,物理結(jié)皮發(fā)育經(jīng)歷形成-破壞-再形成的過程[29],因此,裸土坡面的流態(tài)變化趨勢起伏較大,而生物結(jié)皮坡面由于結(jié)皮的結(jié)構(gòu)穩(wěn)定,徑流的流態(tài)變化相對較平緩。該結(jié)果也表明,生物結(jié)皮改變了暴雨條件下坡面徑流的流態(tài),裸土坡面徑流流態(tài)隨降雨歷時變化劇烈,而生物結(jié)皮坡面流態(tài)隨降雨歷時變化平緩,生物結(jié)皮減少了徑流侵蝕力,避免坡面出現(xiàn)溝蝕,這也是降雨中后期生物結(jié)皮坡面較裸土坡面產(chǎn)流更少的原因之一。

    在本次試驗中,降雨歷時不同,生物結(jié)皮對產(chǎn)流的影響截然不同,裸土坡面在降雨前30 min的產(chǎn)流率與生物結(jié)皮坡面并無顯著差異,降雨45 min后顯著高于高蓋度結(jié)皮坡面。試驗結(jié)果表明,考慮生物結(jié)皮因素的野外降雨試驗,降雨歷時不少于45 min得到的結(jié)果是比較可靠的。除降雨歷時外,坡面產(chǎn)流還受降雨類型、土壤性質(zhì)、結(jié)皮類型等多種因素的影響[30-32],當以上影響因素發(fā)生改變時,本試驗揭示的產(chǎn)流規(guī)律是否仍然適用,是下一步亟待解決的問題,還需要做更多的試驗研究。

    4 結(jié) 論

    1)生物結(jié)皮坡面較翻耕后的裸土坡面顯著降低了初始產(chǎn)流時間,裸土初始產(chǎn)流時間是生物結(jié)皮坡面的1.59~3.04倍。生物結(jié)皮蓋度與初始產(chǎn)流時間之間呈顯著的負相關(guān)關(guān)系。

    2)降雨15 min時與60 min時生物結(jié)皮對坡面徑流的影響發(fā)生逆轉(zhuǎn),90 mm/h的雨強下,當降雨歷時為15 min時,生物結(jié)皮坡面較裸土坡面增加75.42%的徑流;當降雨歷時為60 min時,生物結(jié)皮坡面徑流較裸土坡面降低52.42%。

    3)生物結(jié)皮影響了土壤水分入滲速率,導(dǎo)致生物結(jié)皮坡面與裸土坡面隨降雨歷時變化的產(chǎn)流特征出現(xiàn)差異,裸土坡面60 min時的入滲率較15 min時降低了34.30%,高蓋度生物結(jié)皮坡面只降低了6.38%。

    4)生物結(jié)皮對坡面入滲產(chǎn)流的影響與降雨歷時有極大的關(guān)系,降雨歷時不同,很可能得到截然相反的結(jié)論,因此考慮生物結(jié)皮因素的野外降雨試驗,降雨歷時應(yīng)不少于45 min。這為解釋生物結(jié)皮影響坡面入滲產(chǎn)流方面存在的分歧提供了數(shù)據(jù)參考。

    [1]Belnap J. Biological soil crusts in deserts: A short review of their role in soil fertility, stabilization, and water relations[J]. Archiv Fur Hydrobiologie Algological Studies, 2003, 148(4): 113-126.

    [2]張元明. 荒漠地表生物土壤結(jié)皮的微結(jié)構(gòu)及其早期發(fā)育特征[J]. 科學(xué)通報,2005,50(1):42-47.

    [3]楊凱,趙允格,馬昕昕. 黃土丘陵區(qū)生物土壤結(jié)皮層水穩(wěn)性[J]. 應(yīng)用生態(tài)學(xué)報,2012,23(1):173-177.

    Yang Kai, Zhao Yunge, Ma Xinxin. Water stability of biological soil crusts in hilly regions of Loess Plateau, Northwest China[J]. Chinese Journal of Applied Ecology, 2012, 23(1): 173-177. (in Chinese with English abstract)

    [4]張培培,趙允格,王媛,等. 黃土高原丘陵區(qū)生物結(jié)皮土壤的斥水性[J]. 應(yīng)用生態(tài)學(xué)報,2014,25(3):657-663.

    Zhang Peipei, Zhao Yunge, Wang Yuan, et al. Impact of biological soil crusts on soil water repellence in the hilly Loess Plateau region, China[J]. Chinese Journal of Applied Ecology, 2014, 25(3): 657-663. (in Chinese with English abstract)

    [5]王媛,趙允格,姚春竹,等. 黃土丘陵區(qū)生物土壤結(jié)皮表面糙度特征及影響因素[J]. 應(yīng)用生態(tài)學(xué)報,2014,25(3):647-656.

    Wang Yuan, Zhao Yunge, Yao Chunzhu, et al. Surface roughness characteristics of biological soil crusts and its influencing factors in the hilly Loess Plateau region, China[J]. Chinese Journal of Applied Ecology, 2014, 25(3): 647-656. (in Chinese with English abstract)

    [6]Loope W L, Gifford G F. Influence of a soil microfloral crust on select properties of soils under pinyon-juniper in southeastern Utah[J]. Journal of Soil & Water Conservation, 1972, 27(4): 164-167.

    [7]Galun M, Bubrick P, Garty J. Structural and metabolic diversity of two desert lichen populations[J]. Journal- Hattori Botanical Laboratory, 1982, 53: 321-324.

    [8]Greene R, Chartres C J, Hodgkinson K C. The effects of fire on the soil in a degraded semiarid woodland. I. Cryptogam cover and physical and micromorphological properties[J]. Australian Journal of Soil Research, 1990, 28(5): 755-777.

    [9]Eldridge D J. Cryptogams, vascular plants, and soil hydrological relations: some preliminary results from the semiarid woodlands of eastern Australia[J]. Great Basin Naturalist, 1993, 53(1): 48-58.

    [10]肖波,趙允格,邵明安. 黃土高原侵蝕區(qū)生物結(jié)皮的人工培育及其水土保持效應(yīng) [J]. 草地學(xué)報,2008,16(1):28-33.

    Xiao Bo, Zhao Yunge, Shao Mingan. Artificial cultivation of biological soil crust and its effects on soil and water conservation in water-wind erosion crisscross region of Loess Plateau, China[J]. Acta Agrestia Sinica, 2008, 16(1): 28-33. (in Chinese with English abstract)

    [11]Al-Qinna M I, Abu-Awwad A M. Infiltration rate measurements in arid soils with surface crust[J]. Irrigation Science, 1998, 18(2): 83-89.

    [12]Rodríguez-Caballero E, Cantón Y, Jetten V. Biological soil crust effects must be included to accurately model infiltration and erosion in drylands: An example from Tabernas Badlands[J]. Geomorphology, 2015, 241: 331-342.

    [13]李新榮,馬鳳云,龍立群,等. 沙坡頭地區(qū)固沙植被土壤水分動態(tài)研究[J]. 中國沙漠,2001,21(3):217-222.

    Li Xinrong, Ma Fengyun, Long Liqun, et al. Soil water dynamics under sand-fixing vegetation in Shapotou area[J]. Journal of Desert Research, 2001, 21(3): 217-222. (in Chinese with English abstract)

    [14]Zhao Yunge, Xu Mingxiang. Runoff and soil loss from revegetated grasslands in thehilly Loess Plateau region, China: Influence of biocrust patches and plant canopies[J]. Journal of Hydrologic Engineering, 2013, 18(4): 387-393.

    [15]Booth W E. Algaeas pioneers in plant succession and their importance in erosion control[J]. Ecology, 1941, 22(1): 38-46.

    [16]Eldridge D J, Tozer M E, Slangen S. Soil hydrology is independent of microphytic crust cover: Further evidence from a wooded semiarid Australian rangeland[J]. Arid Soil Research & Rehabilitation, 1997, 11(2): 113-126.

    [17]Chamizo S, Cantón Y, Rodríguez-Caballero E, et al. Biocrusts positively affect the soil water balance in semiarid ecosystems[J]. Ecohydrology, 2016, 9(7): 1208-1221.

    [18]王一賀,趙允格,李林,等. 黃土高原不同降雨量帶退耕地植被-生物結(jié)皮的分布格局 [J]. 生態(tài)學(xué)報,2016,36(2):377-386.

    Wang Yihe, Zhao Yunge, Li Lin, et al. Distribution patterns and spatial variability of vegetation and biocrusts in revegetated lands in different rainfall zones of the Loess Plateau region, China[J]. Acta Ecologica Sinica, 2016, 36(2): 377-386. (in Chinese with English abstract)

    [19]葉菁,卜崇峰,楊永勝,等. 翻耙干擾下生物結(jié)皮對水分入滲及土壤侵蝕的影響[J]. 水土保持學(xué)報,2015,29(3):22-26.

    Ye Jing, Bu Chongfeng, Yang Yongsheng, et al. Influences of biological soil crusts on water infiltration and soil erosion under harrowing disturbance[J]. Journal of Soil and Water Conservation, 2015, 29(3): 22-26. (in Chinese with English abstract)

    [20]謝申琦,高麗倩,趙允格,等. 模擬降雨條件下生物結(jié)皮坡面產(chǎn)流產(chǎn)沙對雨強的響應(yīng)[J]. 應(yīng)用生態(tài)學(xué)報,2019,30(2):35-41.

    Xie Shenqi, Gao Liqian, Zhao Yunge, et al. Responses of runoff and soil loss from biological soil crustal slope to rainfall intensity under simulated rainfall[J]. Chinese Journal of Applied Ecology, 2019, 30(2): 35-41. (in Chinese with English abstract)

    [21]謝云,林小鵑,劉英娜,等. 槽式擺噴頭下噴式人工模擬降雨機的雨強及其空間分布率定[J]. 水土保持通報,2008,28(4):1-6.

    Xie Yun, Lin Xiaojuan, Liu Yingna, et al. Calibration of simulated rainfall intensity and its spatial distribution for trough rainfall simulator[J]. Bulletin of Soil and Water Conservation, 2008, 28(4): 1-6. (in Chinese with English abstract)

    [22]許明祥,劉國彬,溫仲明,等. 黃土丘陵區(qū)小流域土壤特性時空動態(tài)變化研究[J]. 水土保持通報,2000,20(1):20-23.

    Xu Mingxiang, Liu Guobin, Wen Zhongming, et al. Temporal and spatial variation of soil characters in small catchment of Loess Hilly areas[J]. Bulletin of Soil and Water Conservation, 2000, 20(1): 20-23. (in Chinese with English abstract)

    [23]張寬地,王光謙,孫曉敏,等. 模擬植被覆蓋條件下坡面流水動力學(xué)特性[J]. 水科學(xué)進展,2014,25(6):825-834.

    Zhang Kuandi, Wang Guangqian, Sun Xiaomin, et al. Hydraulic characteristic of overland flow under different vegetation coverage[J].Advances in Water Science, 2014, 25(6): 825-834. (in Chinese with English abstract)

    [24]肖波,趙允格,邵明安. 陜北水蝕風(fēng)蝕交錯區(qū)兩種生物結(jié)皮對土壤飽和導(dǎo)水率的影響[J]. 農(nóng)業(yè)工程學(xué)報,2007,23(12):35-40.

    Xiao Bo, Zhao Yunge, Shao Mingan. Effects of biological soil crust on saturated hydraulic conductivity in water-wind erosion crisscross region, North of Shaanxi Province, China[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2007, 23(12): 35-40. (in Chinese with English abstract)

    [25]劉俊娥,王占禮,高素娟. 黃土坡面片蝕過程試驗研究[J]. 水土保持學(xué)報,2011,25(3):35-39.

    Liu June, Wang Zhanli, GaoSujuan. Experimental study of sheet erosion processes on Loess Hill slope[J]. Journal of Soil and Water Conservation, 2011, 25(3): 35-39. (in Chinese with English abstract)

    [26]秦寧強,趙允格. 生物土壤結(jié)皮對雨滴動能的響應(yīng)及削減作用[J]. 應(yīng)用生態(tài)學(xué)報,2011,22(9):2259-2264.

    Qin Ningqiang, Zhao Yunge. Responses of biological soil crust to and its relief effect on raindrop kinetic energy[J]. Chinese Journal of Applied Ecology, 2011, 22(9): 2259-2264. (in Chinese with English abstract)

    [27]Savat J, Ploey J D. Sheetwash and rill development by surface flow[C]//Badland Geomorphology and Piping, Norwich: Geobooks, 1998: 113-126.

    [28]張科利,鐘德鈺. 黃土坡面溝蝕發(fā)生機理的水動力學(xué)試驗研究[J]. 泥沙研究,1998(3):76-82.

    Zhang Keli, Zhong Deyu. The study of gully erosion mechanism of the hydrodynamic experimental on loess slope surface[J]. Journal of Sediment Research, 1998(3): 76-82. (in Chinese with English abstract)

    [29]蔡強國,陸兆熊. 黃土發(fā)育表土結(jié)皮過程和微結(jié)構(gòu)分析的試驗研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,1996(4):363-370.

    Cai Qiangguo, Lu Zhaoxiong. Experimental study of surface crusts process and microfabric analysis on loess soil[J]. Journal of Basic Science and Engineering, 1996(4): 363-370. (in Chinese with English abstract)

    [30]鄔鈴莉,王云琦,王晨灃,等. 降雨類型對北方土石山區(qū)坡面土壤侵蝕的影響[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(24):157-164.

    Wu Lingli, Wang Yunqi, Wang Chenfeng, et al. Effect of rainfall patterns on hillslope soil erosion in rocky mountain area of North China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(24): 157-164. (in Chinese with English abstract)

    [31]吳秋菊,吳發(fā)啟,王林華. 土壤結(jié)皮坡面流水動力學(xué)特征[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(1):73-80.

    Wu Qiuju, Wu Faqi, Wang Linhua. Hydrodynamic characteristics of overland flow under soil crusts condition[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(1): 73-80. (in Chinese with English abstract)

    [32]彭浩,李忠武,劉春,等. 湘中低山丘陵區(qū)坡面產(chǎn)流輸沙對降雨、土壤類型及水保措施的綜合響應(yīng)特征[J]. 水土保持學(xué)報,2019,33(2):60-67.

    Peng Hao, Li Zhongwu, Liu Chun, et al. Comprehensive response characteristics of runoff and sediment yield on rainfall, soil type and water conservation measures in hilly area of central Hunan Province[J]. Journal of Soil and Water Conservation, 2019, 33(2): 60-67. (in Chinese with English abstract)

    Characteristics of runoff on biological soil crust slope in different rainfall durations

    Yang Kai1,2, Zhao Jun2, Zhao Yunge2,4※, Zhang Zihui1, Sun Hui3, Gu Kangmin1, Guo Yali1, Wang Shanshan4

    (1.,,712100,;2,,,712100,; 3.,,712100,; 4.,100049,)

    The Loess Plateau in China is one of the most severely eroded regions of the world. Since the implementation of “Grain for Green” ecological restoration project, biological soil crusts (biocrusts) were widely distributed in this region, which significantly affected surface runoff. Numerous studies have explored the effect of biocrusts on runoff. However, the conclusions were still widely different. In the Loess Plateau region, rainfall is mostly concentrated in June to September, and the rainfall duration is not fixed, which may affect the runoff characteristics of biocrust slopes. This study investigated characteristics of runoff from biocrust slope in different rainfall durations in the Loess Plateau region by using artificial simulated rainfall experiment. The experiment was conducted in the revegetated grassland of northern Shaanxi Provence, China. The experiment site was about 80 m×20 m, and the slope gradient was approximately 15°. The biocrust types were mainly moss crust and moss cyanobacteria mixed crust in this site and their average coverage was 79.2%. The dimensions of the experimental plots were 10 m×2.1 m (length×width).Canopy of higher plants in the plots was removed with scissors. According to the range of local biocrust coverage, two treatments were set: 1) slopes with undisturbed biocrust as a high coverage biocrust (the average biocrust coverage were 79.2%); 2) the slopes with removal of a part of the biocrusts by shovels, which simulated the low biocrust cover situation (the average biocrust coverage were 43.6%). Meanwhile, ploughing plots were set as the control group. The rainfall intensity was set as 90mm/h and the duration was1 hour. The results showed that the initial runoff time of biocrust slope was significantly reduced compared to the bare soil slope. The initial runoff yield time of bare soil was 1.59-3.04 times that of the biocrust slopes. There was a significant negative correlation between biocrust coverage and initial runoff generation time; Conclusion of the influence of biocrusts on runoff yield was contradictory during the first 15 min and 60 min. For 90 mm/h rainfall intensity, runoff from biocrust slope increased by 75.42% compared to bare soil when the rainfall duration was the 15 min. While, runoff from biocrust slope was decreased by 52.42% compared to the bare soil when the rainfall lasted to 60 min; the infiltration rate of soil moisture was affected by biocrusts. The infiltration rate of bare soil slope with 60 min rainfall was 34.30% lower than that with 15 min. The infiltration rate of high coverage biocrust slope with 60 min rainfall was only 6.38% lower than that with 15 min, which may cause the difference of runoff yield between bare soil slope and biocrust slope; the effect of biocrust on slope infiltration and runoff is closely related to rainfall duration. Different periods of rainfall are likely to lead to inconsistent conclusions. Therefore, the duration of rainfall experiment considering the factors of biocrust should be no less than 45 min. The study provides scientific evidences for explaining the differences in the effect of biocrusts on infiltration and runoff, and further clarifies the hydrological effect of biocrusts in arid and semi-arid areas.

    runoff; hydrodynamics; biological soil crust; initial runoff time; runoff coefficient

    楊 凱,趙 軍,趙允格,張子輝,孫 會,谷康民,郭雅麗,王閃閃. 生物結(jié)皮坡面不同降雨歷時的產(chǎn)流特征[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(23):135-141.doi:10.11975/j.issn.1002-6819.2019.23.017 http://www.tcsae.org

    Yang Kai, Zhao Jun, Zhao Yunge, Zhang Zihui, Sun Hui, Gu Kangmin, Guo Yali, Wang Shanshan. Characteristics of runoff on biological soil crust slope in different rainfall durations[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(23): 135-141. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.23.017 http://www.tcsae.org

    2019-07-12

    2019-11-12

    國家自然科學(xué)基金重點項目(41830758);國家自然科學(xué)基金面上項目(41571268)

    楊 凱,博士生,主要從事生物結(jié)皮抗侵蝕方面的研究。Email:yangkai0409@163.com。

    趙允格,博士,研究員,博士生導(dǎo)師,主要從事生物結(jié)皮生態(tài)功能及土壤系統(tǒng)中水、溶質(zhì)運移研究。Email:zyunge@ms.iswc.ac.cn

    10.11975/j.issn.1002-6819.2019.23.017

    S157.1; S181

    A

    1002-6819(2019)-23-0135-07

    猜你喜歡
    結(jié)皮土坡產(chǎn)流
    產(chǎn)流及其研究進展
    NaCl與Na2SO4風(fēng)沙土鹽結(jié)皮力學(xué)特性比較分析?
    不同坡面單元人工降雨產(chǎn)流試驗與分析
    黃土高原不同發(fā)育階段生物結(jié)皮的導(dǎo)水和持水特征
    黃土高原典型生物結(jié)皮對坡面產(chǎn)流產(chǎn)沙過程的影響
    北京山區(qū)側(cè)柏林地坡面初始產(chǎn)流時間影響因素
    地表粗糙度對黃土坡面產(chǎn)流機制的影響
    土壤風(fēng)蝕環(huán)境下鹽結(jié)皮研究及展望?
    上海SMP公園土坡場
    戶外探險(2015年2期)2015-02-26 23:30:21
    一个人看视频在线观看www免费 | 欧美精品啪啪一区二区三区| 成人亚洲精品av一区二区| 黄片大片在线免费观看| 亚洲七黄色美女视频| 999久久久精品免费观看国产| 国产乱人伦免费视频| 在线a可以看的网站| 在线十欧美十亚洲十日本专区| 香蕉久久夜色| 两性午夜刺激爽爽歪歪视频在线观看| 国产亚洲精品av在线| 国内久久婷婷六月综合欲色啪| 午夜福利成人在线免费观看| 国产麻豆成人av免费视频| 欧美成狂野欧美在线观看| 国产精华一区二区三区| 欧美激情在线99| 嫩草影院入口| 免费大片18禁| 日韩欧美免费精品| 久久久国产精品麻豆| 亚洲精品国产精品久久久不卡| 久久精品国产清高在天天线| 欧美日本亚洲视频在线播放| 精品久久久久久,| 精品久久久久久,| 国产三级在线视频| 日本与韩国留学比较| 久久国产乱子伦精品免费另类| 日韩欧美国产一区二区入口| 老熟妇仑乱视频hdxx| 亚洲成人免费电影在线观看| 国语自产精品视频在线第100页| 亚洲成av人片在线播放无| 欧美3d第一页| 99国产精品一区二区蜜桃av| 中文字幕人成人乱码亚洲影| 国产精品乱码一区二三区的特点| 色视频www国产| 乱人视频在线观看| 久久欧美精品欧美久久欧美| av专区在线播放| 久99久视频精品免费| 十八禁网站免费在线| 麻豆久久精品国产亚洲av| 国产激情欧美一区二区| 亚洲成人久久爱视频| 欧美日韩乱码在线| av天堂在线播放| 村上凉子中文字幕在线| 亚洲精品日韩av片在线观看 | 国产一区二区亚洲精品在线观看| 久久香蕉国产精品| 天美传媒精品一区二区| 午夜免费激情av| 久久久久性生活片| 欧美3d第一页| 尤物成人国产欧美一区二区三区| 精品久久久久久久久久免费视频| 成人欧美大片| 国产av一区在线观看免费| www.熟女人妻精品国产| 亚洲欧美日韩高清在线视频| 夜夜躁狠狠躁天天躁| 天堂网av新在线| 久久久久免费精品人妻一区二区| 99国产精品一区二区三区| 岛国在线观看网站| 亚洲av五月六月丁香网| 午夜两性在线视频| 国产视频一区二区在线看| 国内精品美女久久久久久| 国产精品综合久久久久久久免费| ponron亚洲| 亚洲成人精品中文字幕电影| 男人舔女人下体高潮全视频| 色吧在线观看| 国产精品三级大全| 叶爱在线成人免费视频播放| 制服丝袜大香蕉在线| 国产精品久久久久久人妻精品电影| 一边摸一边抽搐一进一小说| 真人一进一出gif抽搐免费| 狂野欧美白嫩少妇大欣赏| 国产中年淑女户外野战色| 99久国产av精品| 欧洲精品卡2卡3卡4卡5卡区| 日日夜夜操网爽| 小蜜桃在线观看免费完整版高清| 国产精品野战在线观看| 好看av亚洲va欧美ⅴa在| 国产高清videossex| 中文字幕精品亚洲无线码一区| 久久久久久久精品吃奶| 宅男免费午夜| 精品人妻偷拍中文字幕| 丁香六月欧美| 蜜桃亚洲精品一区二区三区| 99热这里只有精品一区| 亚洲七黄色美女视频| 久久这里只有精品中国| 亚洲成人久久性| 无限看片的www在线观看| av片东京热男人的天堂| 午夜福利免费观看在线| 舔av片在线| 午夜福利视频1000在线观看| e午夜精品久久久久久久| 色尼玛亚洲综合影院| 一区二区三区国产精品乱码| 亚洲专区国产一区二区| 成人国产一区最新在线观看| 99久久精品热视频| 欧美一区二区国产精品久久精品| 男女视频在线观看网站免费| 真人一进一出gif抽搐免费| 久久久久性生活片| 天堂影院成人在线观看| 免费看光身美女| av在线蜜桃| www日本黄色视频网| 波多野结衣高清作品| 国产伦人伦偷精品视频| a在线观看视频网站| 亚洲人成伊人成综合网2020| 亚洲第一电影网av| 中文字幕av在线有码专区| av天堂中文字幕网| 人人妻,人人澡人人爽秒播| 久久亚洲精品不卡| 乱人视频在线观看| 青草久久国产| www国产在线视频色| 久久人妻av系列| 成人永久免费在线观看视频| 色哟哟哟哟哟哟| 91av网一区二区| 女生性感内裤真人,穿戴方法视频| 91久久精品国产一区二区成人 | 黄色丝袜av网址大全| 欧美日本视频| 国产亚洲精品久久久久久毛片| 色在线成人网| 国产成年人精品一区二区| 观看美女的网站| 亚洲av免费高清在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲av成人精品一区久久| 国产探花极品一区二区| 夜夜看夜夜爽夜夜摸| 日本与韩国留学比较| 久久中文看片网| 国产伦人伦偷精品视频| 国产亚洲av嫩草精品影院| 91麻豆av在线| 国产激情偷乱视频一区二区| 小说图片视频综合网站| 亚洲欧美日韩无卡精品| 深夜精品福利| 日韩av在线大香蕉| 久久草成人影院| 国产97色在线日韩免费| 观看免费一级毛片| 成人永久免费在线观看视频| 国产精品综合久久久久久久免费| 国产午夜精品久久久久久一区二区三区 | 国产三级中文精品| 精品电影一区二区在线| 久久这里只有精品中国| 免费大片18禁| 久久久国产成人免费| 高清在线国产一区| 一个人看视频在线观看www免费 | 免费看十八禁软件| 欧美色欧美亚洲另类二区| 最新美女视频免费是黄的| 国内精品一区二区在线观看| 国产爱豆传媒在线观看| 久久亚洲精品不卡| 欧美一区二区亚洲| 在线观看免费午夜福利视频| 国产在视频线在精品| 天天一区二区日本电影三级| 又粗又爽又猛毛片免费看| 88av欧美| 亚洲狠狠婷婷综合久久图片| 久久久久九九精品影院| 国产v大片淫在线免费观看| 少妇的逼水好多| 白带黄色成豆腐渣| 精品人妻1区二区| 亚洲18禁久久av| 成人无遮挡网站| 最新在线观看一区二区三区| 最近最新中文字幕大全电影3| 一个人看视频在线观看www免费 | 久9热在线精品视频| 国产久久久一区二区三区| 99精品久久久久人妻精品| 女人被狂操c到高潮| 啦啦啦观看免费观看视频高清| 757午夜福利合集在线观看| 村上凉子中文字幕在线| 亚洲内射少妇av| 十八禁人妻一区二区| 婷婷丁香在线五月| 亚洲国产欧美网| 动漫黄色视频在线观看| 一区福利在线观看| 亚洲欧美日韩卡通动漫| 国产日本99.免费观看| 一级毛片女人18水好多| 一级黄色大片毛片| 久久精品国产自在天天线| 国产精品三级大全| 欧美色视频一区免费| 最近最新中文字幕大全电影3| 手机成人av网站| av中文乱码字幕在线| 变态另类丝袜制服| tocl精华| 久久久久精品国产欧美久久久| 欧美绝顶高潮抽搐喷水| 亚洲真实伦在线观看| 亚洲av成人av| 男插女下体视频免费在线播放| 色av中文字幕| 韩国av一区二区三区四区| 日本黄大片高清| 搡女人真爽免费视频火全软件 | 日韩中文字幕欧美一区二区| 精品一区二区三区人妻视频| 日韩中文字幕欧美一区二区| 欧美区成人在线视频| 国产黄片美女视频| 国产精品1区2区在线观看.| 波野结衣二区三区在线 | 亚洲中文字幕日韩| 全区人妻精品视频| 亚洲国产高清在线一区二区三| 中文亚洲av片在线观看爽| 欧美日韩综合久久久久久 | 午夜日韩欧美国产| 久久久久久久久中文| 国产成+人综合+亚洲专区| 无限看片的www在线观看| 在线十欧美十亚洲十日本专区| 特级一级黄色大片| 国产精品一区二区三区四区免费观看 | 日日摸夜夜添夜夜添小说| 亚洲电影在线观看av| 在线视频色国产色| 极品教师在线免费播放| 久久久国产成人免费| 99热这里只有精品一区| 91在线观看av| 久久久久久九九精品二区国产| 亚洲18禁久久av| 1024手机看黄色片| 国产aⅴ精品一区二区三区波| 婷婷亚洲欧美| 国产精品久久久久久亚洲av鲁大| www.999成人在线观看| 久久午夜亚洲精品久久| 亚洲不卡免费看| 成人高潮视频无遮挡免费网站| 亚洲自拍偷在线| 国产精品自产拍在线观看55亚洲| 美女 人体艺术 gogo| a级一级毛片免费在线观看| 亚洲成人久久性| АⅤ资源中文在线天堂| 亚洲成av人片免费观看| 亚洲在线自拍视频| 欧美日韩黄片免| 亚洲av成人av| 国产淫片久久久久久久久 | 亚洲成av人片免费观看| 亚洲人成网站高清观看| 又黄又粗又硬又大视频| 又黄又粗又硬又大视频| 天堂av国产一区二区熟女人妻| 成人av在线播放网站| 18禁国产床啪视频网站| 免费av观看视频| 狂野欧美激情性xxxx| 香蕉av资源在线| 一个人免费在线观看的高清视频| 日韩欧美精品免费久久 | 尤物成人国产欧美一区二区三区| 亚洲av五月六月丁香网| 亚洲精品一区av在线观看| 90打野战视频偷拍视频| 日韩高清综合在线| av女优亚洲男人天堂| 午夜激情福利司机影院| 国产v大片淫在线免费观看| 黄色丝袜av网址大全| a级一级毛片免费在线观看| 99国产精品一区二区三区| 精品久久久久久久毛片微露脸| 亚洲国产精品sss在线观看| 五月伊人婷婷丁香| 亚洲,欧美精品.| 老汉色av国产亚洲站长工具| 亚洲av日韩精品久久久久久密| 午夜福利免费观看在线| 最近最新中文字幕大全免费视频| 青草久久国产| 欧美日韩福利视频一区二区| 91麻豆av在线| 尤物成人国产欧美一区二区三区| 久久精品国产综合久久久| 亚洲精品在线观看二区| 88av欧美| 日韩有码中文字幕| 麻豆一二三区av精品| 国产午夜精品久久久久久一区二区三区 | 国内精品久久久久精免费| 成人av一区二区三区在线看| 天天添夜夜摸| 国产精品亚洲av一区麻豆| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| 一边摸一边抽搐一进一小说| 国产色爽女视频免费观看| 91字幕亚洲| or卡值多少钱| 精品无人区乱码1区二区| 老司机深夜福利视频在线观看| 国产不卡一卡二| 禁无遮挡网站| 国产精品电影一区二区三区| 嫩草影视91久久| 搡老妇女老女人老熟妇| 久久伊人香网站| 亚洲成人久久性| 18禁在线播放成人免费| 久久久国产成人精品二区| 伊人久久精品亚洲午夜| 国产蜜桃级精品一区二区三区| 亚洲成av人片免费观看| 成人欧美大片| 99国产极品粉嫩在线观看| 99国产精品一区二区蜜桃av| 天堂影院成人在线观看| 日韩人妻高清精品专区| 欧美日韩乱码在线| 观看免费一级毛片| 久久欧美精品欧美久久欧美| 嫩草影院入口| 成人一区二区视频在线观看| 国产午夜精品久久久久久一区二区三区 | 夜夜看夜夜爽夜夜摸| 色综合站精品国产| 性欧美人与动物交配| 久久国产精品影院| 亚洲熟妇熟女久久| 两人在一起打扑克的视频| 99久久九九国产精品国产免费| 久久精品国产亚洲av香蕉五月| 欧美乱妇无乱码| 久久久久免费精品人妻一区二区| 久久亚洲精品不卡| 亚洲人成电影免费在线| 又黄又粗又硬又大视频| 欧美黑人欧美精品刺激| 成年免费大片在线观看| 日本成人三级电影网站| 国产成人福利小说| 一卡2卡三卡四卡精品乱码亚洲| 国产极品精品免费视频能看的| 免费大片18禁| 九色国产91popny在线| 日韩欧美精品v在线| 俺也久久电影网| 级片在线观看| 久久久久久九九精品二区国产| 搡老妇女老女人老熟妇| 制服人妻中文乱码| 久久久久久人人人人人| 欧美zozozo另类| 欧美一级a爱片免费观看看| 一进一出抽搐gif免费好疼| 观看美女的网站| 黄色丝袜av网址大全| 久久欧美精品欧美久久欧美| 亚洲精品久久国产高清桃花| 无限看片的www在线观看| 波多野结衣高清作品| 精品99又大又爽又粗少妇毛片 | 可以在线观看的亚洲视频| 99国产精品一区二区三区| 少妇的逼好多水| 狂野欧美激情性xxxx| 国产午夜福利久久久久久| 国产黄色小视频在线观看| 国产 一区 欧美 日韩| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 国产伦在线观看视频一区| 99精品欧美一区二区三区四区| 女人被狂操c到高潮| 亚洲 欧美 日韩 在线 免费| 我要搜黄色片| 久久人妻av系列| 亚洲乱码一区二区免费版| 青草久久国产| 久久精品亚洲精品国产色婷小说| 亚洲av免费在线观看| 亚洲一区高清亚洲精品| 久久精品91蜜桃| 偷拍熟女少妇极品色| 天堂网av新在线| 九九久久精品国产亚洲av麻豆| 两人在一起打扑克的视频| 岛国在线观看网站| 亚洲av美国av| 人妻久久中文字幕网| 国产爱豆传媒在线观看| 欧美成人性av电影在线观看| 黑人欧美特级aaaaaa片| 亚洲国产精品sss在线观看| 麻豆国产av国片精品| 特级一级黄色大片| 长腿黑丝高跟| 精品久久久久久成人av| 国产三级中文精品| 国产成人影院久久av| 亚洲国产日韩欧美精品在线观看 | 亚洲男人的天堂狠狠| 男人舔奶头视频| 国产亚洲av嫩草精品影院| 国产精品国产高清国产av| 丰满乱子伦码专区| 国产乱人伦免费视频| 有码 亚洲区| 丰满人妻熟妇乱又伦精品不卡| 国产乱人视频| 国产成人影院久久av| 99国产精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 国产av一区在线观看免费| 精品人妻1区二区| 一本一本综合久久| www日本在线高清视频| 90打野战视频偷拍视频| 九色国产91popny在线| 欧美av亚洲av综合av国产av| 在线观看日韩欧美| 欧美区成人在线视频| 亚洲av第一区精品v没综合| 日日夜夜操网爽| 久久婷婷人人爽人人干人人爱| 最新在线观看一区二区三区| 国产高潮美女av| 十八禁人妻一区二区| 国产精品国产高清国产av| 亚洲精品一卡2卡三卡4卡5卡| 婷婷精品国产亚洲av| 我要搜黄色片| 一级毛片高清免费大全| 欧美黄色片欧美黄色片| 99热这里只有精品一区| 男女午夜视频在线观看| 十八禁人妻一区二区| 成年版毛片免费区| 欧美性感艳星| 两个人看的免费小视频| 999久久久精品免费观看国产| 亚洲无线在线观看| 蜜桃亚洲精品一区二区三区| 老司机在亚洲福利影院| 欧美午夜高清在线| 久久九九热精品免费| 窝窝影院91人妻| 亚洲精品一区av在线观看| 亚洲七黄色美女视频| 特大巨黑吊av在线直播| 免费av不卡在线播放| 哪里可以看免费的av片| 国产日本99.免费观看| 亚洲色图av天堂| 亚洲,欧美精品.| 在线免费观看不下载黄p国产 | 午夜福利18| 免费人成视频x8x8入口观看| 中文字幕人妻熟人妻熟丝袜美 | 国产一区二区亚洲精品在线观看| 国产97色在线日韩免费| 香蕉久久夜色| 高清日韩中文字幕在线| 国产精品永久免费网站| 五月伊人婷婷丁香| 亚洲一区高清亚洲精品| 日韩欧美国产一区二区入口| 国产淫片久久久久久久久 | 日韩大尺度精品在线看网址| 久久精品国产亚洲av香蕉五月| 老司机午夜福利在线观看视频| 国产精品综合久久久久久久免费| 桃红色精品国产亚洲av| 欧美激情久久久久久爽电影| 国产黄色小视频在线观看| netflix在线观看网站| 90打野战视频偷拍视频| 国产伦精品一区二区三区四那| 床上黄色一级片| 人妻久久中文字幕网| 国产成人系列免费观看| 少妇人妻精品综合一区二区 | 欧美xxxx黑人xx丫x性爽| 91九色精品人成在线观看| 噜噜噜噜噜久久久久久91| 一级毛片高清免费大全| 亚洲精品国产精品久久久不卡| 看免费av毛片| 亚洲精品色激情综合| 国产高清videossex| av欧美777| 成人亚洲精品av一区二区| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 亚洲精品亚洲一区二区| 国产精品美女特级片免费视频播放器| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 波野结衣二区三区在线 | 校园春色视频在线观看| 3wmmmm亚洲av在线观看| 亚洲欧美精品综合久久99| 在线视频色国产色| 成人鲁丝片一二三区免费| 99国产极品粉嫩在线观看| 一级黄片播放器| 午夜福利高清视频| 在线观看免费午夜福利视频| 亚洲专区中文字幕在线| 欧美日韩福利视频一区二区| 1000部很黄的大片| 久久久色成人| 久久久久久九九精品二区国产| 久99久视频精品免费| 日韩精品青青久久久久久| 青草久久国产| 色哟哟哟哟哟哟| 99在线视频只有这里精品首页| 色哟哟哟哟哟哟| 好男人在线观看高清免费视频| 国产免费一级a男人的天堂| av女优亚洲男人天堂| 日韩人妻高清精品专区| 99热6这里只有精品| 神马国产精品三级电影在线观看| 老司机午夜十八禁免费视频| 久久欧美精品欧美久久欧美| 欧美日韩福利视频一区二区| 中文亚洲av片在线观看爽| 亚洲性夜色夜夜综合| 精品久久久久久久久久免费视频| 亚洲国产精品999在线| 国产淫片久久久久久久久 | 国产精品女同一区二区软件 | 变态另类成人亚洲欧美熟女| 亚洲精品国产精品久久久不卡| 欧美国产日韩亚洲一区| 欧美+亚洲+日韩+国产| 亚洲国产精品成人综合色| 小蜜桃在线观看免费完整版高清| 蜜桃亚洲精品一区二区三区| 成年女人看的毛片在线观看| а√天堂www在线а√下载| 国内精品久久久久精免费| 午夜a级毛片| 亚洲欧美激情综合另类| 久久香蕉精品热| 岛国在线免费视频观看| 小蜜桃在线观看免费完整版高清| 香蕉av资源在线| 狂野欧美激情性xxxx| 色哟哟哟哟哟哟| 中国美女看黄片| 亚洲国产精品合色在线| 国产一区二区三区在线臀色熟女| 观看免费一级毛片| 亚洲国产精品sss在线观看| 亚洲国产色片| 日本成人三级电影网站| 日本a在线网址| 亚洲国产欧美网| 变态另类丝袜制服| 窝窝影院91人妻| 听说在线观看完整版免费高清| 伊人久久精品亚洲午夜| 国产精品 国内视频| 久久精品91蜜桃| 亚洲国产色片| 人人妻,人人澡人人爽秒播| 成年女人毛片免费观看观看9| 男人舔女人下体高潮全视频| 亚洲国产精品久久男人天堂| 美女黄网站色视频| 日日夜夜操网爽| 欧美性感艳星| 国产黄片美女视频| 毛片女人毛片| 波多野结衣高清无吗| 美女高潮的动态| 搡女人真爽免费视频火全软件 | 日本黄色片子视频| 精品国产三级普通话版| 欧美大码av| 午夜福利免费观看在线| www日本在线高清视频|