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

    干濕交替次數(shù)對土壤干密度和飽和導(dǎo)水率的影響

    2022-07-13 03:15:52孫召軍李金山賈艷輝
    灌溉排水學(xué)報 2022年6期
    關(guān)鍵詞:導(dǎo)水率原狀變異性

    孫召軍,李金山,賈艷輝,李 浩

    干濕交替次數(shù)對土壤干密度和飽和導(dǎo)水率的影響

    孫召軍1,2,李金山1*,賈艷輝1,李 浩1

    (1.中國農(nóng)業(yè)科學(xué)院 農(nóng)田灌溉研究所 河南省節(jié)水農(nóng)業(yè)重點試驗室,河南 新鄉(xiāng) 453002;2.中國農(nóng)業(yè)科學(xué)院 研究生院,北京 100081)

    研究干濕交替對原狀土和擾動土土壤干密度()和飽和導(dǎo)水率()的影響特征。通過室內(nèi)試驗對不同干濕交替條件下原狀土和擾動土的干密度和飽和導(dǎo)水率測定,分析了不同干濕交替過程對土壤干密度和飽和導(dǎo)水率的變化特征及規(guī)律的影響。原狀土干密度在1.41~1.89 g/cm3之間波動,擾動土干密度在1.47~1.85 g/cm3之間波動;原狀土在第1次干濕交替后的干密度較原始干密度有大幅上升,隨著干濕交替次數(shù)的增加逐漸趨于小范圍波動,總體變異系數(shù)為8%;擾動土總體處于上升趨勢,總體變異系數(shù)為5%;原狀土干密度變化幅度顯著大于擾動土(<0.05)。原狀土飽和導(dǎo)水率在0.88~3.22 cm/h之間波動,擾動土飽和導(dǎo)水率在0.21~3.54 cm/h之間波動;原狀土在經(jīng)歷不同干濕交替次數(shù)后呈上升的趨勢,總體變異系數(shù)為39%;擾動土在經(jīng)歷不同干濕交替次數(shù)后飽和導(dǎo)水率呈現(xiàn)先下降后上升的趨勢,總體變異系數(shù)為94%;原狀土飽和導(dǎo)水率變化幅度顯著小于擾動土(<0.05)。本試驗條件下,土壤經(jīng)歷干濕交替后形成了較多水道,土壤透水性增強;大田翻耕后,表層土壤處于擾動松散狀態(tài),干密度較小,飽和導(dǎo)水率較大,而第1次灌水以后,干密度變大,飽和導(dǎo)水率降低,且干密度和飽和導(dǎo)水率在后續(xù)的生產(chǎn)過程中趨于穩(wěn)定狀態(tài)。

    干濕交替;土壤干密度;飽和導(dǎo)水率;土壤水分特性

    0 引 言

    土壤干密度(Dry density, DD)是土壤的基本物理性質(zhì)之一,主要反映土壤內(nèi)部總孔隙體積的變化[1],受母質(zhì)組成成分、氣候變化、人類活動、有機質(zhì)量、生物作用等的綜合影響[2]。土壤飽和導(dǎo)水率()是土壤重要的物理性質(zhì)之一,研究土壤在飽和狀態(tài)下的土壤入滲速率可以反映土壤水分運移規(guī)律,主要受土壤質(zhì)地、孔隙度、結(jié)構(gòu)等因素的影響[3-8]。土壤粒徑組成(Particle size distribution, PSD)是反應(yīng)土壤結(jié)構(gòu)特征的重要因素之一,對土壤物理力學(xué)性質(zhì)(孔隙比、液塑限、強度和變形等)和水理性質(zhì)(滲透性、濕陷性和保水性等)有非常大的影響[9-11]。

    【研究意義】在降雨、灌溉和蒸發(fā)的影響下,土壤經(jīng)歷不斷吸水和失水的過程。土壤在失水時收縮,進而產(chǎn)生表面裂縫和體積收縮;土壤吸水時膨脹,土壤體積增大且裂縫消失,含水率越大膨脹越大,嚴重影響和改變了土壤結(jié)構(gòu)和孔隙分布,進而改變了土壤水分運動規(guī)律[12-14]。因此,研究不同干濕交替次數(shù)對土壤干密度和飽和導(dǎo)水率的影響機理,對提高水資源利用效率、發(fā)展節(jié)水灌溉以及農(nóng)業(yè)可持續(xù)發(fā)展具有重要意義。

    【研究進展】目前已有較多針對土壤水分運動規(guī)律在空間上變異的研究。自然條件下表層土壤因受到干濕交替、土壤裂隙、蟲洞和根孔等因素的影響,大孔隙數(shù)量較多,土壤入滲能力相對較高;甘淼等[15]和祝景彬等[16]研究發(fā)現(xiàn)淺層土壤較為緊實,土壤干密度較大,但隨著土層深度的增加,土壤干密度有減小的趨勢。許振欣等[17]研究發(fā)現(xiàn)隨著土壤深度的增加,飽和導(dǎo)水率呈先增加后減少的趨勢,且表層土壤(0~10 cm)飽和導(dǎo)水率明顯高于底層土壤(10~50 cm);毛娜等[18]研究發(fā)現(xiàn)20 cm以上的飽和導(dǎo)水率顯著高于20~200 cm深度的土壤,干密度與飽和導(dǎo)水率顯著負相關(guān);李濤等[19]研究發(fā)現(xiàn)飽和導(dǎo)水率具有很強的空間變異性,與黏粒、粉粒、砂粒和有機質(zhì)量具有一定相關(guān)性,而與干密度幾乎沒有空間相關(guān)性。胡偉等[4]認為盡管原狀土和擾動土在干密度上保持一致,但由于原狀土含有非均勻的孔隙,而擾動土的孔隙分布均勻,即二者大小孔隙的分布不一致,造成二者的飽和導(dǎo)水率測定值不一樣;Sandin等[20]利用X射線層分析技術(shù)研究飽和導(dǎo)水率動態(tài)變化和土壤結(jié)構(gòu)性孔隙相關(guān)關(guān)系的過程中發(fā)現(xiàn),土壤飽和導(dǎo)水率與孔隙連通率相關(guān)。王力等[21]發(fā)現(xiàn)土壤滲透性受粒度組成的影響最為直接,二者聯(lián)系也最為密切。土壤干密度具有明顯的空間變異性,且隨著土壤深度的增加呈減小趨勢,這是由于機械作用和人為活動對表層土壤影響較大,對深層土壤影響較小;飽和導(dǎo)水率的大小受土壤孔隙的影響,干密度越小,土壤松散,孔隙就越多,故飽和導(dǎo)水率與干密度呈顯著的負相關(guān)關(guān)系,但這些研究大多針對土壤水分運動規(guī)律在空間上的變異,較少關(guān)注在時間上的變異。因此,為了探求土壤水分運動規(guī)律在時間和空間上的變異,有必要對不同干濕交替次數(shù)下土壤干密度和飽和導(dǎo)水率的變化規(guī)律進行研究。

    農(nóng)業(yè)生產(chǎn)過程中人類活動、自然環(huán)境等影響因素顯著改變了土壤孔隙結(jié)構(gòu)的形成和重組,進而影響土壤水分入滲運動參數(shù)的變化態(tài)勢?!厩腥朦c】目前,關(guān)于土壤水分物理參數(shù)在時間上變異的研究較少,使得大田中翻耕、灌溉、降雨、蒸發(fā)等對土壤入滲速率、水運動特征、養(yǎng)分運移規(guī)律等認識不夠深刻?!緮M解決的關(guān)鍵問題】因此,本研究利用原狀土和擾動土對比試驗,通過室內(nèi)短時間干濕交替的方法探求土壤水分物理參數(shù)在時間上的變異特征,得到土壤入滲特性在作物生育期內(nèi)經(jīng)歷翻耕、降雨、灌溉和蒸發(fā)等過程中的動態(tài)變化及規(guī)律,以期為田間灌溉的智能化管理和節(jié)水模型構(gòu)建提供技術(shù)支撐。

    1 材料與方法

    1.1 研究區(qū)域

    研究區(qū)域位于河南省新鄉(xiāng)縣七里營鎮(zhèn)中國農(nóng)業(yè)科學(xué)院新鄉(xiāng)綜合試驗基地(E113.5°,N35.9°)。新鄉(xiāng)縣地處古黃河沖積半原的北翼和太行山前沖洪積扇的南緣地帶,海拔70~82 m,屬于暖溫帶大陸性季風(fēng)氣候,年均降水量580 mm,多集中在7—8月,年均蒸發(fā)量2 000 mm,年均氣溫14.1 ℃,1月最冷,平均氣溫0.7 ℃,7月最熱,平均氣溫27.1 ℃,年均日照時間2 407.7 h,年均無霜期200.5 d。以冬小麥-夏玉米連作為主。試驗區(qū)土壤主要為粉砂質(zhì)黏壤土,土壤干密度在0.7 m土層深度內(nèi)平均值為1.47 g/cm3,0~70 cm土層平均田間持水率為26.04%,土壤物理性質(zhì)見表1[22]。

    表1 試驗田土壤物理性質(zhì)

    1.2 試驗設(shè)計和研究方法

    1.2.1 供試材料

    為保證土樣的代表性,選取畦田中間部位作為取樣點,按照1.5 m×1.5 m的樣方進行取土。原狀土取土?xí)r,剝離表層10 cm后取樣,10~30 cm為第1層,30~50 cm為第2層,50~70 cm為第3層,測定飽和導(dǎo)水率土樣利用環(huán)刀(體積120 cm3,內(nèi)徑61.8 cm,高4 cm)每層取5個樣共15個,測定干密度土樣利用環(huán)刀每層取5個樣共15個;擾動土取土?xí)r,每層土壤取1袋,風(fēng)干后碾碎過0.9 mm孔篩,根據(jù)每層土壤的干密度和風(fēng)干后土壤的含水率進行分層裝填到環(huán)刀內(nèi),每層取5個樣,共15個。

    1.2.2 試驗設(shè)計

    1)干濕交替共做4次交替,即每組5個中第1個土樣不做干濕處理,第2個土樣做1次干濕交替,以此類推,原狀土和擾動土干濕交替方法相同。均先將土壤在試驗室級高純水中浸水至飽和狀態(tài)后用烘箱105 ℃烘至恒質(zhì)量,再將烘干的土泡飽和后烘干,此步驟重復(fù)。此試驗中所用試驗室級高純水是利用易普易達EPED-X2-10T制備,制備過程中在設(shè)備沖洗完成后開始制水,根據(jù)操作手冊要求純度值在18 MΩ·cm以上的水可用。具體處理方法見表2。

    表2 干濕交替處理方法

    2)擾動土處理方法。將原狀土室內(nèi)風(fēng)干,碾碎過0.9 mm孔篩,裝入密封袋中。裝填環(huán)刀之前測土壤含水率和原狀土初始干密度,擾動土取樣部位的初始干密度第1、第2、第3層分別為1.53、1.56、1.47 g/cm3,裝填前的含水率第1、第2、第3層分別為1.74%、1.92%、1.93%。每層擾動土的干密度與原狀土對應(yīng),裝土質(zhì)量增加擾動土含水率。在裝入120 cm3環(huán)刀過程中,逐層填土均勻夯實,每層1 cm厚,接觸面打毛,防止出現(xiàn)明顯分層現(xiàn)象。

    1.2.3 觀測項目與方法

    1)土壤干密度的測定。初始干密度采用環(huán)刀法測定,即土樣烘箱中105 ℃烘至恒質(zhì)量,得到土塊質(zhì)量,除以環(huán)刀體積得到土壤初始干密度;干濕交替后的干密度采用蠟封法測定,取土塊的一部分在105 ℃的烘箱中烘干至恒質(zhì)量后,放入干燥皿中冷卻,冷卻后在土塊系上細繩并稱量其質(zhì)量,將土塊放入剛過熔點的蠟液中,土塊浸沒后立即取出,檢查土塊周圍的蠟?zāi)?,若有起泡?yīng)用針刺破,再用蠟液補平,再次冷卻后稱蠟封土塊的質(zhì)量,最后將土塊完全浸入到試驗室級高純水中,利用排水法測定土塊體積。排水法測定土塊體積是將蠟封后的土塊冷卻后放入試驗室級高純水中,測得其在水中的質(zhì)量,空氣中質(zhì)量與水中質(zhì)量的差值,即為常溫下蠟封后土塊的體積,減去蠟的體積即得到土塊的體積。

    環(huán)刀法測干密度計算式:

    式中:d為土的干密度(g/cm3);0為干土質(zhì)量(g);為環(huán)刀的體積(cm3)。

    蠟封法測干密度公式:

    式中:0為干土質(zhì)量(g);n為試樣加蠟質(zhì)量(g);nw為試樣加蠟在水中的質(zhì)量(g);為土的濕密度(g/cm3);d為土的干密度(g/cm3);wT為純水在℃時的密度(g/cm3),準確至0.01 g/cm3;n為蠟的密度(g/cm3),準確至0.01 g/cm3;為含水率(%)。

    2)飽和導(dǎo)水率的測定。土樣在試驗室級高純水中充分飽和后,用定水頭法測定飽和導(dǎo)水率,使用馬氏瓶供水,利用滲透儀(TST 55型)進行測定,水頭高度維持在15 cm[23]。每個土壤樣品重復(fù)6次,每次重復(fù)間隔3 min(即0.05 h),取其平均值。

    式中:為土壤飽和導(dǎo)水率(cm/h);為流量,滲透過一定截面積的水量(cm3);為飽和土層厚度,滲透經(jīng)過的距離(cm);為環(huán)刀橫截面積(cm2);為水頭差(cm);為滲透過水量時所需時間(h)。

    3)土壤粒度的測定。將原始土樣中的大顆粒進行分散細碎,放置風(fēng)干后,過0.9 mm的孔篩,將土樣放入百特公司BT-401型超聲波循環(huán)分散器循環(huán)至百特公司BT-9300HT型激光粒度分布儀中測定土壤粒度,測量范圍0.1~1 000 μm,循環(huán)介質(zhì)采用試驗室級高純水。為減小試驗誤差,每個土樣測定重復(fù)5次,取均值。按照美國制分類,黏粒粒徑小于2 μm,粉粒粒徑在2~50 μm之間,砂粒粒徑在50~1 000 μm之間。

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

    使用Excel 2019處理數(shù)據(jù),利用SPSS 23.0進行變異系數(shù)、方差分析和相關(guān)性分析,使用Origin 2019畫圖。

    2 結(jié)果與分析

    2.1 粒徑組成

    原狀土和擾動土室內(nèi)干濕交替后粒徑組成的統(tǒng)計結(jié)果(表3)顯示,粒徑組成在經(jīng)歷干濕交替后變異系數(shù)非常小,總體呈弱變異性,故本文試驗暫不考慮干濕交替對粒徑組成產(chǎn)生的影響,在后續(xù)的試驗中會針對該內(nèi)容進行試驗分析論證。其中有部分數(shù)據(jù)的變異系數(shù)變化幅度較大,如原狀土第2層砂粒占比在經(jīng)歷干濕交替后變異系數(shù)達到18.34%,從總體數(shù)據(jù)特征分析是因為測驗誤差引起的。

    2.2 干濕交替后土壤干密度的變化

    圖1為原狀土和擾動土在經(jīng)歷室內(nèi)干濕交替后的干密度變化。原狀土和擾動土室內(nèi)干濕交替后干密度的統(tǒng)計結(jié)果(表4)顯示,原狀土第1、第2、第3層干密度均值分別為1.80、1.67、1.58 g/cm3,根據(jù)數(shù)據(jù)分析顯示每層均為弱變異性,但是最小值與最大值差異較大;擾動土第1、第2、第3層干密度均值分別為1.70、1.66、1.60 g/cm3,根據(jù)數(shù)據(jù)分析顯示每層均為弱變異性,但是最小值與最大值差異較大。試驗表明,原狀土干濕交替后的干密度大幅增加,因此認為干濕交替對土壤的孔隙結(jié)構(gòu)影響非常大。

    表3 室內(nèi)干濕交替后粒徑組成的統(tǒng)計特征

    注 變異系數(shù)≤10%時為弱變異性,10%<變異系數(shù)<100%時為中變異性,變異系數(shù)≥100%時為強變異性。

    試驗表明,擾動土與原狀土的變化趨勢相似,土壤干濕交替后,干密度大幅增加。原狀土和擾動土的干密度在T1處理后有了大幅增加,這是因為干濕交替后土壤孔隙大幅減小;原狀土和擾動土在T2處理或T3處理以后,呈輕微減小趨勢,這是因為隨著干濕交替次數(shù)的增加,部分土壤孔隙遭到破壞,并形成了新的孔隙結(jié)構(gòu);在T3、T4處理過程中原狀土和擾動土的每層土壤在干濕交替后處在小范圍的波動中,干濕交替次數(shù)對土壤孔隙的增加影響逐漸減小,形成了穩(wěn)定的孔隙結(jié)構(gòu),但第1層擾動土在T3、T4處理過程中干密度變化較大,說明干濕交替過程中孔隙收縮力產(chǎn)生的作用效果不穩(wěn)定。原狀土的干密度在T0處理時第1層最大,第2層次之,第3層最小,在T4處理時干密度也呈現(xiàn)此規(guī)律,說明原狀土在干濕交替過程中孔隙收縮力產(chǎn)生的作用效果比較穩(wěn)定;擾動土的干密度在T0處理時第2層最大、第1層次之、第3層最小,在T4處理時干密度第1層最大、第2層次之、第3層最小,這說明擾動土由于初始狀態(tài)中孔隙結(jié)構(gòu)遭到破壞,在干濕交替過程中孔隙收縮力產(chǎn)生的作用效果不穩(wěn)定,隨著干濕交替次數(shù)的增加,新形成的孔隙結(jié)構(gòu)逐漸趨于穩(wěn)定。

    圖1 干濕交替后原狀土和擾動土干密度變化趨勢

    表4 室內(nèi)干濕交替后干密度的統(tǒng)計特征

    注 變異系數(shù)≤10%時為弱變異性,10%<變異系數(shù)<100%時為中變異性,變異系數(shù)≥100%時為強變異性。

    2.3 干濕交替后土壤飽和導(dǎo)水率的變化

    圖2為原狀土和擾動土在經(jīng)歷室內(nèi)干濕交替后的飽和導(dǎo)水率變化。

    通過原狀土和擾動土室內(nèi)干濕交替后飽和導(dǎo)水率的統(tǒng)計結(jié)果(表5)顯示,原狀土第1、第2、第3層飽和導(dǎo)水率均值分別為1.94、1.18、1.38 cm/h,飽和導(dǎo)水率在每層呈中變異性,說明在經(jīng)歷干濕交替以后土壤飽和導(dǎo)水率均產(chǎn)生了較大的變化,第1層土樣中摻雜了作物秸稈等異物且孔隙結(jié)構(gòu)受外界影響較大,一定程度上影響了水道的形成,第2層受外界因素影響較第1層小,較第3層大且土樣中雜物較少密實程度介于第1層和第3層,第3層土樣受外界因素影響較小且無雜物密實程度高,所以第1層和第3層的飽和導(dǎo)水率相比第2層要低一些;擾動土第1、第2、第3層飽和導(dǎo)水率均值分別為0.87、0.31、2.17 cm/h,飽和導(dǎo)水率在每層呈中變異性,說明土壤擾動后在經(jīng)歷泡水、烘干過程重組了土壤孔隙結(jié)構(gòu)分布格局,填充時第2層干密度最大、第1層次之、第3層最小,第2層干密度最大且雜物較第1層少,土壤密實程度較高,第1層干密度較第2層小但雜物較多對孔隙結(jié)構(gòu)的重組影響較大,第3層干密度最小且雜物最少故變異系數(shù)介于第1、第2層之間,說明擾動狀態(tài)下的土壤密實度越高土壤孔隙結(jié)構(gòu)重組影響就越小,由于擾動土的結(jié)構(gòu)遭到徹底破壞,裝填時形成了新的結(jié)構(gòu),重組后的孔隙結(jié)構(gòu)在經(jīng)歷干濕交替后的變化較原狀土小。趨勢圖顯示,原狀土隨著干濕交替次數(shù)的增加,第1、第2、第3層土壤的飽和導(dǎo)水率均有上升的趨勢,這是因為在干濕交替過程中土壤孔隙率發(fā)生了變化,造成部分孔隙產(chǎn)生破壞,但上升速度有所差異,這與干密度和孔隙結(jié)構(gòu)有關(guān)系;第1層土壤的飽和導(dǎo)水率在T3處理時急速上升,T4處理后飽和導(dǎo)水率又急速下降;第2層土壤的飽和導(dǎo)水率在T3處理時有上升,T4處理后飽和導(dǎo)水率緩慢下降;第3層土壤的飽和導(dǎo)水率處于波動狀態(tài),總體有上升趨勢,T4后飽和導(dǎo)水率最大。

    擾動土第1層和第3層的試驗數(shù)據(jù)顯示,T0處理飽和導(dǎo)水率均大于干濕交替后的土壤;第2層土壤中,T0處理飽和導(dǎo)水率只高于T1、T2處理的土壤,但是差值不大。所有的數(shù)據(jù)顯示,T1處理后土壤的飽和導(dǎo)水率急速下降,隨著干濕交替次數(shù)的增加,飽和導(dǎo)水率逐漸增大,這是因為T1處理后土壤孔隙率大幅減小對土壤孔隙結(jié)構(gòu)產(chǎn)生了較大的影響,干濕交替對土壤結(jié)構(gòu)的影響是先收縮后逐漸膨脹。不同干濕交替次數(shù)后,泡水后土體呈不同程度的膨脹。不同干濕交替后飽和導(dǎo)水率的變化趨勢非常明顯,呈先下降后上升的趨勢,在第2層中由于干密度最大,飽和導(dǎo)水率變化很小,規(guī)律依然明顯。

    圖2 干濕交替后原狀土和擾動土飽和導(dǎo)水率變化趨勢

    表5 室內(nèi)干濕交替后飽和導(dǎo)水率的統(tǒng)計特征

    注 變異系數(shù)≤10%時為弱變異性,10%<變異系數(shù)<100%時為中變異性,變異系數(shù)≥100%時為強變異性;由于飽和導(dǎo)水率較小,為了更加清晰的顯示變化規(guī)律,本文以厘米/小時(cm/h)為計算單位。

    3 討 論

    3.1 干濕交替對土壤干密度的影響

    原狀土和擾動土的干密度均值逐層遞減且第1層最高,這是由于農(nóng)業(yè)生產(chǎn)活動中,機械作用和人類活動使土壤壓實,導(dǎo)致干密度較大,鞠忻倪等[24]的研究成果與之相似。干密度均值之所以逐層降低,還因為在干濕交替過程中土壤表現(xiàn)出吸水土體膨脹、失水土壤收縮的現(xiàn)象,且干濕交替在一定程度上加劇了土體的變形各向異性,土壤的裂隙會隨干濕交替次數(shù)的增加而逐漸增大[25]。在干濕交替過程中,干密度呈快速增加和緩慢增加2個階段,最終趨于穩(wěn)定狀態(tài)[26-28],這可能是由于土壤吸水到失水的過程中,水分子對土粒產(chǎn)生壓力,使土粒變得緊實,減少了土壤孔隙[29-31]。本文結(jié)合以往研究中土壤干密度在時間和空間上的變異性,以原狀土和擾動土對比,在時空尺度上開展研究,原狀土和擾動土室內(nèi)干濕交替后,原狀土第1、第2、第3層干密度分別在1.56~1.89、1.45~1.80、1.41~1.66 g/cm3之間波動,變異系數(shù)分別是7%、7%、6%,這是由于受機械作用和人類活動的作用,土壤表層在一定程度上被壓實,隨著土層深度的增加受影響的程度越來越小;擾動土第1、第2、第3層干密度分別在1.53~1.85、1.56~1.72、1.47~1.68 g/cm3之間波動,變異系數(shù)分別是6%、3%、5%,這是由于表層土壤含有植物根莖以及土壤顆粒成分組成的影響,擾動土的干密度均值也是逐層遞減趨勢。原狀土干密度變化幅度顯著大于擾動土(<0.05),說明土壤經(jīng)過擾動破壞后,土壤干密度結(jié)構(gòu)產(chǎn)生重組,并且翻耕后機械作用和人類活動對土壤的干密度影響較大。

    3.2 干濕交替對土壤飽和導(dǎo)水率的影響

    原狀土第1層的飽和導(dǎo)水率均值最高,這是由于表層經(jīng)歷翻耕、機械和人為作用,對土壤結(jié)構(gòu)產(chǎn)生了一定的破壞,連通孔隙較多;擾動土第1層的飽和導(dǎo)水率均值低于第3層,這是由于擾動破壞后重新裝填重組了土壤的孔隙結(jié)構(gòu)。原狀土第1、第2、第3層飽和導(dǎo)水率均值呈先下降后上升趨勢,這是因為干濕交替后土壤中失去抗拉強度的裂隙會有所閉合,加之土體內(nèi)部密實度已經(jīng)增大,所以水分滲流阻礙增強,滲透性質(zhì)相應(yīng)減弱;隨著干濕交替次數(shù)的增加,土體中存在的縱橫向裂隙在干濕交替作用下產(chǎn)生了連通效應(yīng),成為新的水分滲流通道,飽和導(dǎo)水率呈增加趨勢[32]。針對水分滲漏通道問題,已有研究中表明土壤中有機質(zhì)量和礦物質(zhì)量對非毛細滲流通道的形成存在一定的影響[33],這是因為土壤水力傳導(dǎo)度與土壤質(zhì)地、土壤孔隙結(jié)構(gòu)、含水率和基質(zhì)勢存在相關(guān)關(guān)系[34],土壤中有機質(zhì)量越多團粒結(jié)構(gòu)越多且穩(wěn)定性越好[35],所以有機質(zhì)量高的土壤入滲能力也大,本文試驗沒有對有機質(zhì)、礦物指標等測定,是因為本研究是針對土壤水分物理參數(shù)作為重點研究對象,在后續(xù)的試驗中會補充有機質(zhì)、礦物指標等的測定內(nèi)容。擾動土第1、第2、第3層飽和導(dǎo)水率均值也呈先下降后上升的趨勢,但第3層飽和導(dǎo)水率最大,這是由于干濕交替會影響?zhàn)ち!⒎哿?、砂粒的懸浮和遷移,而擾動條件下的顆粒遷移較原狀土更活躍,巫尚蔚等[36]認為粉粒量的增加弱化了原有的強力鏈。從團聚體形成角度上來說,原狀土和擾動土在經(jīng)歷了干濕交替后飽和導(dǎo)水率呈現(xiàn)相似的規(guī)律,是因為干濕交替對團聚體形成初期產(chǎn)生非常大的影響,土體的密實程度大幅提高,隨著干濕交替次數(shù)的增加而逐漸平穩(wěn)[37]。研究表明,在干濕交替1~2次促進大團聚體的形成[38],土體在失水時,團聚體粒徑減??;在吸水時,團聚體粒徑增大[39-40]。所以在首次干濕交替時,土體形成了新的團粒結(jié)構(gòu),導(dǎo)致土壤密實度迅速增加,飽和導(dǎo)水率迅速下降,隨后逐漸上升并趨于穩(wěn)定(圖2)。原狀土第1、第2、第3層飽和導(dǎo)水率分別在1.24~3.22、0.88~1.75、0.99~1.75 cm/h之間變化,變異系數(shù)分別是36%、30%、21%;擾動土第1、第2、第3層飽和導(dǎo)水率分別在0.42~1.70、0.21~0.39、0.41~3.54 cm/h之間變化,變異系數(shù)分別是55%、22%、52%。原狀土飽和導(dǎo)水率變化幅度顯著小于擾動土(<0.05),這是因為干濕交替對土壤孔隙結(jié)構(gòu)進行了破壞重組,加之土壤顆粒組成的影響,形成了先下降后上升的趨勢;擾動土飽和導(dǎo)水率低于原狀土,是與土壤團聚體的形成有關(guān),擾動土在經(jīng)歷破壞后形成了新的團粒結(jié)構(gòu),土壤密實程度大于原狀土。

    4 結(jié) 論

    1)不同干濕交替次數(shù)下,原狀土和擾動土的干密度整體變化規(guī)律基本一致,均上升后趨于穩(wěn)定。原狀土和擾動土干密度總體變異系數(shù)分別為8%、5%,原狀土干密度變化幅度顯著大于擾動土(<0.05)。

    2)不同干濕交替次數(shù)下,原狀土的孔隙結(jié)構(gòu)產(chǎn)生了變化,使飽和導(dǎo)水率增大,干濕交替過程中連通孔隙增多。原狀土和擾動土飽和導(dǎo)水率總體變異系數(shù)分別為39%、94%,原狀土飽和導(dǎo)水率變化幅度顯著小于擾動土(<0.05)。

    3)擾動土在干濕交替后重組的孔隙結(jié)構(gòu),由均勻分布變成了定向分布,并且借助土壤收縮形成了大孔隙,擾動后土壤的密實度降低,水道較多。

    [1] ALAOUI A, LIPIEC J, GERKE H H. A review of the changes in the soil pore system due to soil deformation: A hydrodynamic perspective[J]. Soil and Tillage Research, 2011, 115/116: 1-15.

    [2] HU W, SHAO M A, SI B C. Seasonal changes in surface bulk density and saturated hydraulic conductivity of natural landscapes[J]. European Journal of Soil Science, 2012, 63(6): 820-830.

    [3] SAUER T J, CLOTHIER B E, DANIEL T C. Surface measurements of the hydraulic properties of a tilled and untilled soil[J]. Soil and Tillage Research, 1990, 15(4): 359-369.

    [4] 胡偉, 邵明安, 王全九, 等. 取樣尺寸對土壤飽和導(dǎo)水率測定結(jié)果的影響[J]. 土壤學(xué)報, 2005, 42(6): 1 040-1 043.

    HU Wei, SHAO Ming’an, WANG Quanjiu, et al. Effects of sampling size on measurements of soil saturated hydraulic conductivity[J]. Acta Pedologica Sinica, 2005, 42(6): 1 040-1 043.

    [5] 王子龍, 趙勇鋼, 趙世偉, 等. 退耕典型草地土壤飽和導(dǎo)水率及其影響因素研究[J]. 草地學(xué)報, 2016, 24(6): 1 254-1 262.

    WANG Zilong, ZHAO Yonggang, ZHAO Shiwei, et al. Study on soil saturated hydraulic conductivity and its influencing factors in typical grassland of farmland conversion[J]. Acta Agrestia Sinica, 2016, 24(6): 1 254-1 262.

    [6] 黃婉霞, 鄧羽松, 謝福倩, 等. 花崗巖崩崗不同部位土壤飽和導(dǎo)水率特征及其影響因素[J]. 應(yīng)用生態(tài)學(xué)報, 2020, 31(7): 2 431-2 440.

    HUANG Wanxia, DENG Yusong, XIE Fuqian, et al. Characteristics of soil saturated hydraulic conductivity on different positions and their controlling factors of granite collapsing gullies[J]. Chinese Journal of Applied Ecology, 2020, 31(7): 2 431-2 440.

    [7] 張瑜, 馮紹元, 蘇童, 等. 不同耕作方式下新復(fù)墾區(qū)春玉米試驗研究[J]. 灌溉排水學(xué)報, 2018, 37(8): 64-70.

    ZHANG Yu, FENG Shaoyuan, SU Tong, et al. Impact of different tillage on dynamics of soil moisture and nutrient and the growth of spring maize in newly reclaimed land[J]. Journal of Irrigation and Drainage, 2018, 37(8): 64-70.

    [8] 陳姣, 吳鳳平, 王輝, 等. 生物炭對南方紅壤和水稻土水力學(xué)特性的影響分析[J]. 灌溉排水學(xué)報, 2020, 39(9): 73-80.

    CHEN Jiao, WU Fengping, WANG Hui, et al. The effects of biochar amendment on hydraulic properties of red and paddy soils in Southern China[J]. Journal of Irrigation and Drainage, 2020, 39(9): 73-80.

    [9] 趙茜. 凍融循環(huán)與干濕交替對黃土滲透各向異性及空間分異性的影響研究[D]. 西安: 西安建筑科技大學(xué), 2019.

    ZHAO Qian. Investigation on the influence of freezing-thawing cycle and drying-wetting alternation on the permeability anisotropy and spatial variability of loess[D]. Xi’an: Xi’an University of Architecture and Technology, 2019.

    [10] 馮連昌, 鄭晏武. 中國濕陷性黃土[M]. 北京: 中國鐵道出版社, 1982.

    FENG Lianchang, ZHENG Yanwu. Collapsible loess in China[M]. Beijing: China Railway Publishing House, 1982.

    [11] 唐大雄. 工程巖土學(xué)[M]. 北京: 地質(zhì)出版社, 1987.

    TANG Daxiong. Ock and soil engineering[M]. Beijing: Geological Publishing House, 1987.

    [12] CABIDOCHE Y M, RUY S. Field shrinkage curves of a swelling clay soil: Analysis of multiple structural swelling and shrinkage phases in the prisms of a Vertisol[J]. Soil Research, 2001, 39(1): 143.

    [13] 黃傳琴. 干濕交替過程中土壤脹縮特征的試驗研究[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2008.

    HUANG Chuanqin. Experimental study on soil shrinking and swelling characteristics during the alternative drying and wetting processes[D]. Yangling: Northwest A&F University, 2008.

    [14] CRESCIMANNO G, PROVENZANO G. Soil shrinkage characteristic curve in clay soils: Measurement and prediction[J]. Soil Science Society of America Journal, 1999, 63(1): 25-32.

    [15] 甘淼, 賈玉華, 李同川, 等. 黃土區(qū)坡溝系統(tǒng)容重、飽和導(dǎo)水率和土壤含水量變化分析[J]. 干旱區(qū)研究, 2018, 35(2): 315-324.

    GAN Miao, JIA Yuhua, LI Tongchuan, et al. Variation of soil bulk density, saturated hydraulic conductivity and soil moisture content in a slope-gully unit on the northern loess plateau[J]. Arid Zone Research, 2018, 35(2): 315-324.

    [16] 祝景彬, 賀慧丹, 李紅琴, 等. 牧壓梯度下高寒草甸土壤容重及持水能力的變化特征[J]. 水土保持研究, 2018, 25(5): 66-71.

    ZHU Jingbin, HE Huidan, LI Hongqin, et al. Characteristics of soil bulk density and soil water-holding capacity in alpine meadow under grazing gradients[J]. Research of Soil and Water Conservation, 2018, 25(5): 66-71.

    [17] 許振欣, 鄧羽松, 林立文, 等. 南亞熱帶典型人工林土壤飽和導(dǎo)水率特征及其影響因素研究[J]. 北京林業(yè)大學(xué)學(xué)報, 2021, 43(4): 100-107.

    XU Zhenxin, DENG Yusong, LIN Liwen, et al. Characteristics of soil saturated hydraulic conductivity and its influencing factors of typical plantations in South Subtropical Zone[J]. Journal of Beijing Forestry University, 2021, 43(4): 100-107.

    [18] 毛娜, 黃來明, 邵明安. 黃土區(qū)坡面尺度不同植被類型土壤飽和導(dǎo)水率剖面分布及影響因素[J]. 土壤, 2019, 51(2): 381-389.

    MAO Na, HUANG Laiming, SHAO Ming’an. Profile distribution of soil saturated hydraulic conductivity and controlling factors under different vegetations on slope in loess region[J]. Soils, 2019, 51(2): 381-389.

    [19] 李濤, 李曉琳, 郝新梅, 等. 應(yīng)用土壤質(zhì)地預(yù)測干旱區(qū)葡萄園土壤飽和導(dǎo)水率空間分布[J]. 農(nóng)業(yè)工程學(xué)報, 2014, 30(9): 84-92.

    LI Tao, LI Xiaolin, HAO Xinmei, et al. Predicting spatial distribution of soil saturated hydraulic conductivity by soil texture on vineyard in arid region[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(9): 84-92.

    [20] SANDIN M, KOESTEL J, JARVIS N, et al. Post-tillage evolution of structural pore space and saturated and near-saturated hydraulic conductivity in a clay loam soil[J]. Soil and Tillage Research, 2017, 165: 161-168.

    [21] 王力, 李喜安, 何軍, 等. 不同黏粒含量對黃土滲透系數(shù)影響的試驗研究[J]. 水土保持通報, 2018, 38(2): 95-101.

    WANG Li, LI Xi’an, HE Jun, et al. Experimental study on permeability coefficient of loess with different clay content[J]. Bulletin of Soil and Water Conservation, 2018, 38(2): 95-101.

    [22] 龐婕, 韓其晟, 周爽, 等. 水氣互作對溫室番茄生長、產(chǎn)量和水分利用效率的影響[J]. 灌溉排水學(xué)報, 2022, 41(1): 87-94.

    PANG Jie, HAN Qisheng, ZHOU Shuang, et al. The integrative effects of irrigation and aeration on growth and water use efficiency of greenhouse tomato[J]. Journal of Irrigation and Drainage, 2022, 41(1): 87-94.

    [23] 王翔翔, 程志龍, 蘇光辰, 等. 測定時間對淮北平原砂姜黑土飽和導(dǎo)水率的影響[J]. 灌溉排水學(xué)報, 2021, 40(12): 136-141.

    WANG Xiangxiang, CHENG Zhilong, SU Guangchen, et al. Saturated soil hydraulic conductivity is not constant but varies with the time it is measured [J]. Journal of Irrigation and Drainage, 2021, 40(12): 136-141.

    [24] 鞠忻倪, 賈玉華, 甘淼, 等. 黃土溝壑區(qū)不同地形部位土壤大孔隙特征研究[J]. 土壤學(xué)報, 2018, 55(5): 1 098-1 107.

    JU Xinni, JIA Yuhua, GAN Miao, et al. Characteristics of soil macropores in the gully area of loess plateau as affected by terrain[J]. Acta Pedologica Sinica, 2018, 55(5): 1 098-1 107.

    [25] 蘇立君, 趙茜, 劉華, 等. 干濕交替過程中原狀黃土的脹縮變形特性及裂隙形態(tài)演化規(guī)律[J]. 天津大學(xué)學(xué)報(自然科學(xué)與工程技術(shù)版), 2021, 54(3): 255-267.

    SU Lijun, ZHAO Qian, LIU Hua, et al. Swelling and shrinkage behaviors and evolution law of crack morphology of undisturbed loess during wetting-drying cycles[J]. Journal of Tianjin University (Science and Technology), 2021, 54(3): 255-267.

    [26] 張猛. 干濕交替過程中土壤容重、水分特征曲線和熱特性的動態(tài)變化特征[D]. 北京: 中國農(nóng)業(yè)大學(xué), 2017.

    ZHANG Meng. The dynamics of soil bulk density, water retention curve and thermal properties under wetting and drying cycles[D]. Beijing: China Agricultural University, 2017.

    [27] 馬璠, 王健, 張鵬輝, 等. 干濕交替處理對土壤抗沖性的影響試驗研究[J].中國農(nóng)學(xué)通報, 2014, 30(35): 178-184.

    MA Fan, WANG Jian, ZHANG Penghui, et al. Effects of dry-wet alternation treatment on anti-scourability of soil[J]. Chinese Agricultural Science Bulletin, 2014, 30(35): 178-184.

    [28] 張鵬輝. 土壤干濕交替中抗侵蝕能力的變化及坡面細溝發(fā)育敏感性的影響研究[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2014.

    ZHANG Penghui. The study on the effect of the soil erosion resistance in soil alternate drying-wetting and to the sensitivity of the rill evolution[D]. Yangling: Northwest A & F University, 2014.

    [29] 張素, 熊東紅, 校亮, 等. 干濕交替對土壤性質(zhì)影響的研究[J]. 土壤通報, 2017, 48(3): 762-768.

    ZHANG Su, XIONG Donghong, XIAO Liang, et al. Influence of dry-wet cycling on soil properties[J]. Chinese Journal of Soil Science, 2017, 48(3): 762-768.

    [30] 茹豪, 張建軍, 李玉婷, 等. 黃土高原土壤粒徑分形特征及其對土壤侵蝕的影響[J]. 農(nóng)業(yè)機械學(xué)報, 2015, 46(4): 176-182.

    RU Hao, ZHANG Jianjun, LI Yuting, et al. Fractal features of soil particle size distributions and its effect on soil erosion of loess plateau[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(4): 176-182.

    [31] 秦偉, 左長清, 晏清洪, 等. 紅壤裸露坡地次降雨土壤侵蝕規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報, 2015, 31(2): 124-132.

    QIN Wei, ZUO Changqing, YAN Qinghong, et al. Regularity of individual rainfall soil erosion in bare slope land of red soil[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(2): 124-132.

    [32] 趙茜, 楊金熹, 趙晉萍. 凍融和干濕循環(huán)對原狀黃土滲透系數(shù)的影響[J].中國地質(zhì)災(zāi)害與防治學(xué)報, 2020, 31(2): 119-126.

    ZHAO Qian, YANG Jinxi, ZHAO Jinping. Experimental study on permeability of undisturbed loess under freezing-thawing and drying-wetting cycling conditions[J]. The Chinese Journal of Geological Hazard and Control, 2020, 31(2): 119-126.

    [33] 喬照華. 土壤有機質(zhì)含量與土壤物理性能參數(shù)的相關(guān)性分析[J]. 中國農(nóng)村水利水電, 2008(2): 3-4.

    QIAO Zhaohua. Analysis of interrelation between organic matter content and soil physical property parameters[J]. China Rural Water and Hydropower, 2008(2): 3-4.

    [34] 雷志棟, 楊詩秀, 謝森傳. 土壤水動力學(xué)[M]. 北京: 清華大學(xué)出版社, 1988.

    LEI Zhidong, YANG Shixiu, XIE Senchuan. Soil water dynamics[M]. Beijing: Tsinghua University Press, 1988.

    [35] 黃昌勇. 土壤學(xué)[M]. 北京: 中國農(nóng)業(yè)出版社, 2000.

    HUANG Changyong. Soil science[M]. Beijing: Chinese Agriculture Press, 2000.

    [36] 巫尚蔚, 楊春和, 張超, 等. 粉粒含量對尾礦力學(xué)特性的影響[J]. 巖石力學(xué)與工程學(xué)報, 2017, 36(8): 2 007-2 017.

    WU Shangwei, YANG Chunhe, ZHANG Chao, et al. The effects of silt content on the mechanical properties of tailings[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(8): 2 007-2 017.

    [37] BIRCH H F. The effect of soil drying on humus decomposition and nitrogen availability[J]. Plant and Soil, 1958, 10(1): 9-31.

    [38] DENEF K, SIX J, BOSSUYT H, et al. Influence of dry-wet cycles on the interrelationship between aggregate, particulate organic matter, and microbial community dynamics[J]. Soil Biology and Biochemistry, 2001, 33(12/13): 1 599-1 611.

    [39] 劉艷, 馬茂華, 吳勝軍, 等. 干濕交替下土壤團聚體穩(wěn)定性研究進展與展望[J]. 土壤, 2018, 50(5): 853-865.

    LIU Yan, MA Maohua, WU Shengjun, et al. Soil aggregates as affected by wetting-drying cycle: A review[J]. Soils, 2018, 50(5): 853-865.

    [40] 劉艷. 干濕交替條件下土壤團聚特性及其與植物地上生產(chǎn)力的關(guān)系[D]. 重慶: 中國科學(xué)院大學(xué)(中國科學(xué)院重慶綠色智能技術(shù)研究院), 2018.

    LIU Yan. Under the environmental condition of wetting and drying cycles: Edaphic aggregate characteristics and their relationships with floral above-ground productivity[D]. Chongqing: Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, 2018.

    Effects of Wetting-drying Cycles on Bulk Density and Saturated Hydraulic Conductivity of Soils

    SUN Zhaojun1,2, LI Jinshan1*, JIA Yanhui1, LI Hao1

    (1. Farmland Irrigation Research Institute of Chinese Academy of Agricultural Sciences, Henan Key Laboratory of Water-saving Agriculture, Xinxiang 453002, China; 2.Graduate School of Chinese Academy of Agricultural Sciences, Beijing 100081, China)

    Changes in density and hydraulic conductivity of soil not only modulate subsurface hydrological processes but also affect root growth and root uptake of water and nutrients from soil. Understanding the mechanisms underlying their changes is thus critical to mitigating flooding and improving bioavailability of soil water. In this paper we studied how wetting-drying cycles often seen in the field change bulk density and saturated hydraulic conductivity of the soil.Intact (mimicking zero-tillage) and disturbed (mimicking tillage) soil samples were taken from a cropped field. They were then subject to wetting-drying at different frequencies. During and at the end of the wetting-drying cycle experiment, we measured bulk density and saturated hydraulic conductivity of each soil sample.Bulk density of the intact and disturbed soil samples varied in the range of 1.41~1.89 g/cm3and 1.47~1.85 g/cm3respectively. The bulk density of intact soils after the first wetting-drying cycle increased significantly followed by slight fluctuations in other wetting-drying cycles, with an average coefficient of variation 8%. In contrast, bulk density of the disturbed soils increased with the increase in wetting-drying cycles, with an average coefficient of variation 5%. Overall, the wetting-drying cycles affected the bulk density of intact soils more significantly than the disturbed soils (<0.05). The saturated hydraulic conductivity of the intact and disturbed soils varied in the range of 0.88~3.22 cm/h and 0.21~3.54 cm/h, respectively. The saturated hydraulic conductivity of the intact soils increased with the increase in wetting-drying cycles, with an average coefficient of variation 39%. In contrast, the hydraulic conductivity of the disturbed soil decreased first followed by a decline with the increase in wetting-drying cycles, with an average coefficient of variation 94%. Wetting-drying cycles affected the hydraulic conductivity of the disturbed soils more significantly than the intact soils (<0.05).Wetting-drying cycles resulted in cracks thereby reducing the bulk density and increasing saturated hydraulic conductivity of the soils. The effects, however, depend on tillage. Overall, wetting-drying cycles affect the bulk density of zero-tillage soil and saturated hydraulic conductivity of tilled soil more significantly.

    wetting-drying cycles; bulk density; saturated hydraulic conductivity; soil moisture characteristics

    S275

    A

    10.13522/j.cnki.ggps.2021635

    孫召軍, 李金山, 賈艷輝, 等. 干濕交替次數(shù)對土壤干密度和飽和導(dǎo)水率的影響[J]. 灌溉排水學(xué)報, 2022, 41(6): 89-96.

    SUN Zhaojun, LI Jinshan, JIA Yanhui, et al. Effects of Wetting-drying Cycles on Bulk Density and Saturated Hydraulic Conductivity of Soils[J]. Journal of Irrigation and Drainage, 2022, 41(6): 89-96.

    1672 - 3317(2022)06 - 0089 – 08

    2021-12-22

    中國農(nóng)業(yè)科學(xué)院基本科研業(yè)務(wù)費專項(Y2021YJ07);農(nóng)村應(yīng)急性防洪排澇技術(shù)和裝備研發(fā);河南省重點研發(fā)與推廣專項(科技攻關(guān))項目(202102110279)

    孫召軍(1989-),男。碩士研究生,主要從事節(jié)水灌溉技術(shù)研究。E-mail: sunzj186@163.com

    李金山(1972-),男。研究員,主要從事節(jié)水灌溉技術(shù)研究。E-mail: lijinshan72@126.com

    責(zé)任編輯:白芳芳

    猜你喜歡
    導(dǎo)水率原狀變異性
    陜北山地蘋果園土壤飽和導(dǎo)水率和植物導(dǎo)水率特征
    基于線性源法與圖像處理的土壤飽和導(dǎo)水率快速測量方法
    基于細觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    毓慶宮惇本殿明間原狀陳列的復(fù)原
    紫禁城(2020年8期)2020-09-09 09:38:04
    重現(xiàn)「皇太子之宮」 毓慶宮原狀陳設(shè)復(fù)原記
    紫禁城(2020年8期)2020-09-09 09:37:56
    咳嗽變異性哮喘的預(yù)防和治療
    海城市土壤飽和導(dǎo)水率空間分布特征
    原狀黃土各向異性及卸載變形特征試驗研究
    咳嗽變異性哮喘的中醫(yī)治療近況
    清肺止咳湯治療咳嗽變異性哮喘40例
    or卡值多少钱| 夜夜躁狠狠躁天天躁| 亚洲午夜理论影院| 一区二区三区激情视频| 欧美午夜高清在线| 性插视频无遮挡在线免费观看| 欧美xxxx性猛交bbbb| 精品久久久久久久久久免费视频| 国产午夜福利久久久久久| 中文字幕人成人乱码亚洲影| 天美传媒精品一区二区| 女人被狂操c到高潮| 啦啦啦韩国在线观看视频| 一级黄色大片毛片| 国产欧美日韩一区二区三| 久久久久久久精品吃奶| 免费看日本二区| 99久久99久久久精品蜜桃| 一进一出抽搐gif免费好疼| 亚洲av免费在线观看| 午夜老司机福利剧场| 蜜桃亚洲精品一区二区三区| 免费av观看视频| 日本免费一区二区三区高清不卡| 国产三级黄色录像| 在线免费观看不下载黄p国产 | av福利片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 深爱激情五月婷婷| 成人鲁丝片一二三区免费| 免费av观看视频| 欧美+亚洲+日韩+国产| 九色成人免费人妻av| 欧美成狂野欧美在线观看| 欧美国产日韩亚洲一区| 亚洲成人免费电影在线观看| www.999成人在线观看| 男人的好看免费观看在线视频| 丰满乱子伦码专区| 99久久精品国产亚洲精品| 国产亚洲精品综合一区在线观看| 国产三级中文精品| or卡值多少钱| 特大巨黑吊av在线直播| 国产精品1区2区在线观看.| 午夜福利欧美成人| 午夜福利在线观看吧| 色尼玛亚洲综合影院| 国产白丝娇喘喷水9色精品| 婷婷色综合大香蕉| 99国产精品一区二区蜜桃av| 欧美日韩亚洲国产一区二区在线观看| 欧美激情久久久久久爽电影| 一a级毛片在线观看| a级毛片a级免费在线| 欧美成人性av电影在线观看| 亚洲五月婷婷丁香| 国产欧美日韩一区二区精品| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品乱码久久久v下载方式| 精品一区二区三区av网在线观看| 国产伦在线观看视频一区| 小说图片视频综合网站| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区精品小视频在线| 12—13女人毛片做爰片一| 国产真实伦视频高清在线观看 | 91久久精品电影网| 久久久久久久久久黄片| 一个人免费在线观看的高清视频| 一本一本综合久久| 美女 人体艺术 gogo| 好男人在线观看高清免费视频| 精品午夜福利视频在线观看一区| 看片在线看免费视频| 国产午夜福利久久久久久| 国产野战对白在线观看| 在线免费观看的www视频| 深爱激情五月婷婷| 亚洲成人久久爱视频| 国产av一区在线观看免费| 9191精品国产免费久久| 色哟哟哟哟哟哟| 91麻豆精品激情在线观看国产| 国产一区二区亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器| 国产视频内射| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品合色在线| 亚洲乱码一区二区免费版| 欧美性猛交黑人性爽| 欧美xxxx黑人xx丫x性爽| 欧美色视频一区免费| 直男gayav资源| 韩国av一区二区三区四区| 一进一出抽搐gif免费好疼| 国产国拍精品亚洲av在线观看| 嫁个100分男人电影在线观看| 国产亚洲欧美98| 欧美一级a爱片免费观看看| 午夜两性在线视频| 国产私拍福利视频在线观看| av视频在线观看入口| 简卡轻食公司| 国产综合懂色| 久久香蕉精品热| 日本 av在线| 欧美成人性av电影在线观看| 91久久精品国产一区二区成人| 国产精品久久久久久久久免 | 久久香蕉精品热| 人妻夜夜爽99麻豆av| 九九久久精品国产亚洲av麻豆| av在线天堂中文字幕| www.www免费av| 久久久国产成人精品二区| 色综合婷婷激情| 日本成人三级电影网站| 男女下面进入的视频免费午夜| 亚洲av.av天堂| 亚洲18禁久久av| 丰满乱子伦码专区| 一夜夜www| 精品欧美国产一区二区三| 老司机午夜十八禁免费视频| 国产综合懂色| 亚洲最大成人中文| 亚洲第一区二区三区不卡| 怎么达到女性高潮| 人人妻人人看人人澡| 99国产精品一区二区三区| 欧美三级亚洲精品| 搡老岳熟女国产| 最近视频中文字幕2019在线8| 如何舔出高潮| 亚洲成人久久性| 欧美性感艳星| 欧美另类亚洲清纯唯美| 成人特级黄色片久久久久久久| 亚洲熟妇熟女久久| 午夜福利成人在线免费观看| 日韩欧美 国产精品| 俄罗斯特黄特色一大片| 午夜久久久久精精品| 每晚都被弄得嗷嗷叫到高潮| 一区二区三区四区激情视频 | 天堂√8在线中文| 99热6这里只有精品| 99国产精品一区二区三区| 久久精品国产清高在天天线| 99精品久久久久人妻精品| 好看av亚洲va欧美ⅴa在| 国内揄拍国产精品人妻在线| 一个人免费在线观看电影| 91麻豆av在线| 久久久久久九九精品二区国产| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 久久热精品热| 一边摸一边抽搐一进一小说| 一本综合久久免费| 十八禁网站免费在线| 日本撒尿小便嘘嘘汇集6| 伦理电影大哥的女人| 久久久久久久久久黄片| 成熟少妇高潮喷水视频| 国内少妇人妻偷人精品xxx网站| 搡老妇女老女人老熟妇| av专区在线播放| 婷婷色综合大香蕉| 久久九九热精品免费| netflix在线观看网站| 色噜噜av男人的天堂激情| 成人无遮挡网站| 美女被艹到高潮喷水动态| 亚洲成人中文字幕在线播放| 日本一二三区视频观看| 18+在线观看网站| 午夜a级毛片| 亚洲最大成人手机在线| 国产毛片a区久久久久| 国产精品一区二区三区四区久久| 中出人妻视频一区二区| 国产精品综合久久久久久久免费| 天天躁日日操中文字幕| 国产精品不卡视频一区二区 | 免费人成在线观看视频色| 深夜a级毛片| 国产成人影院久久av| 色综合婷婷激情| 99久久精品热视频| 两人在一起打扑克的视频| 变态另类成人亚洲欧美熟女| 亚洲最大成人av| 高清在线国产一区| 亚洲精品一区av在线观看| 三级毛片av免费| 伦理电影大哥的女人| 亚洲精品粉嫩美女一区| 亚洲精品乱码久久久v下载方式| 91久久精品电影网| 18禁裸乳无遮挡免费网站照片| 看黄色毛片网站| 亚洲片人在线观看| 欧美3d第一页| 动漫黄色视频在线观看| 男女视频在线观看网站免费| 91在线观看av| 久久中文看片网| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 亚洲aⅴ乱码一区二区在线播放| 夜夜躁狠狠躁天天躁| 搡老岳熟女国产| 久9热在线精品视频| 人妻制服诱惑在线中文字幕| 欧美性猛交黑人性爽| 直男gayav资源| 人人妻,人人澡人人爽秒播| 色精品久久人妻99蜜桃| 国产免费男女视频| 国产精品一区二区三区四区久久| 色尼玛亚洲综合影院| 午夜福利在线在线| 麻豆成人午夜福利视频| 欧美精品啪啪一区二区三区| 露出奶头的视频| 久久久久精品国产欧美久久久| 波多野结衣高清无吗| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 亚洲熟妇熟女久久| 真实男女啪啪啪动态图| 日韩成人在线观看一区二区三区| 亚洲成人久久爱视频| 亚洲专区国产一区二区| 欧美性猛交╳xxx乱大交人| 国产又黄又爽又无遮挡在线| 国产午夜精品久久久久久一区二区三区 | 欧美三级亚洲精品| 精品久久久久久久末码| 黄色一级大片看看| 超碰av人人做人人爽久久| 淫妇啪啪啪对白视频| 18禁黄网站禁片免费观看直播| av女优亚洲男人天堂| 成人国产综合亚洲| 免费大片18禁| 色综合婷婷激情| 99热只有精品国产| 国产成人啪精品午夜网站| 国产三级中文精品| 欧美日韩福利视频一区二区| 久久国产精品影院| 欧美zozozo另类| 国产在线男女| 麻豆国产97在线/欧美| 国产精品永久免费网站| 日本与韩国留学比较| 99在线视频只有这里精品首页| 99久久九九国产精品国产免费| 国产精品爽爽va在线观看网站| 精品久久久久久久久久久久久| 亚洲自拍偷在线| 久久久色成人| 中文亚洲av片在线观看爽| 日本a在线网址| 亚洲av中文字字幕乱码综合| 亚洲av成人av| 男女那种视频在线观看| 国产91精品成人一区二区三区| 99久国产av精品| 国产乱人视频| 简卡轻食公司| 成人av一区二区三区在线看| 久久性视频一级片| 一级黄色大片毛片| 亚洲精华国产精华精| 亚洲一区二区三区不卡视频| 欧美日韩中文字幕国产精品一区二区三区| 免费观看精品视频网站| 国产成年人精品一区二区| 欧美日韩亚洲国产一区二区在线观看| 搡老岳熟女国产| x7x7x7水蜜桃| 亚洲专区中文字幕在线| 99精品久久久久人妻精品| 九九在线视频观看精品| 久久久久国内视频| 亚洲最大成人中文| 日本a在线网址| 在线播放无遮挡| 亚洲国产欧洲综合997久久,| 天堂影院成人在线观看| 男女做爰动态图高潮gif福利片| 女人被狂操c到高潮| 大型黄色视频在线免费观看| 国产真实乱freesex| 亚洲 欧美 日韩 在线 免费| 99riav亚洲国产免费| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站 | 亚洲av免费在线观看| 别揉我奶头 嗯啊视频| 亚洲精品一区av在线观看| 亚洲人成电影免费在线| 久久久久久九九精品二区国产| 青草久久国产| 久久精品91蜜桃| 欧美黑人巨大hd| 亚洲欧美日韩东京热| av在线天堂中文字幕| 波多野结衣巨乳人妻| 深夜精品福利| 日本一本二区三区精品| 久久婷婷人人爽人人干人人爱| 亚洲国产欧美人成| 国产又黄又爽又无遮挡在线| 免费电影在线观看免费观看| 午夜久久久久精精品| 精品国产三级普通话版| 久久精品国产亚洲av天美| 亚洲国产精品999在线| 日本撒尿小便嘘嘘汇集6| 国产乱人视频| 一级黄片播放器| 久久热精品热| 中文字幕精品亚洲无线码一区| 悠悠久久av| 午夜精品在线福利| 一区二区三区激情视频| 给我免费播放毛片高清在线观看| 国产精品98久久久久久宅男小说| 夜夜爽天天搞| 国语自产精品视频在线第100页| 69人妻影院| 国产伦在线观看视频一区| 亚洲狠狠婷婷综合久久图片| 9191精品国产免费久久| 偷拍熟女少妇极品色| netflix在线观看网站| 九九热线精品视视频播放| 丁香六月欧美| 麻豆国产97在线/欧美| 看免费av毛片| 国产国拍精品亚洲av在线观看| 成人美女网站在线观看视频| 亚洲熟妇熟女久久| 欧美潮喷喷水| 国产亚洲欧美在线一区二区| 色综合婷婷激情| a级毛片免费高清观看在线播放| 国产精品电影一区二区三区| 久久精品综合一区二区三区| 人人妻人人看人人澡| 久久精品综合一区二区三区| 亚洲午夜理论影院| 欧美日韩乱码在线| 婷婷亚洲欧美| 国模一区二区三区四区视频| 亚洲国产高清在线一区二区三| 亚洲精品久久国产高清桃花| 亚洲国产高清在线一区二区三| 俺也久久电影网| 国产精品久久久久久久电影| 午夜精品在线福利| ponron亚洲| 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| 国内精品久久久久久久电影| 国产男靠女视频免费网站| 亚洲人成伊人成综合网2020| 日韩欧美精品v在线| 国产精品伦人一区二区| 男人舔女人下体高潮全视频| 丰满人妻熟妇乱又伦精品不卡| 久久久成人免费电影| 久久伊人香网站| 色综合亚洲欧美另类图片| 亚洲国产精品久久男人天堂| 丰满人妻一区二区三区视频av| av在线老鸭窝| 精品久久国产蜜桃| 久久久久免费精品人妻一区二区| 91在线精品国自产拍蜜月| 成人高潮视频无遮挡免费网站| 男女之事视频高清在线观看| 亚洲av五月六月丁香网| www.www免费av| 动漫黄色视频在线观看| 能在线免费观看的黄片| 精品人妻1区二区| 成年免费大片在线观看| 日韩中字成人| 精品人妻一区二区三区麻豆 | 亚洲成人久久性| 熟女人妻精品中文字幕| 亚洲精品一区av在线观看| 性欧美人与动物交配| 五月伊人婷婷丁香| 尤物成人国产欧美一区二区三区| 亚洲av美国av| 国产av麻豆久久久久久久| 久久久久久大精品| 啦啦啦韩国在线观看视频| 搡女人真爽免费视频火全软件 | 久久久久久久精品吃奶| 变态另类成人亚洲欧美熟女| 两个人视频免费观看高清| 国产伦精品一区二区三区四那| 有码 亚洲区| 国产在线精品亚洲第一网站| 久久久成人免费电影| 成人亚洲精品av一区二区| 久久精品国产亚洲av香蕉五月| 日韩欧美国产在线观看| 九色成人免费人妻av| 中文字幕av在线有码专区| 俄罗斯特黄特色一大片| 老熟妇乱子伦视频在线观看| 午夜福利视频1000在线观看| 成人毛片a级毛片在线播放| 在线免费观看不下载黄p国产 | 嫩草影院新地址| 国产精品国产高清国产av| 天天躁日日操中文字幕| 日本黄大片高清| 精品人妻熟女av久视频| 无遮挡黄片免费观看| 99久国产av精品| 久久婷婷人人爽人人干人人爱| 九色国产91popny在线| 变态另类丝袜制服| 91九色精品人成在线观看| 51午夜福利影视在线观看| 白带黄色成豆腐渣| 午夜精品在线福利| 日本五十路高清| 最近在线观看免费完整版| 悠悠久久av| 免费在线观看亚洲国产| 五月伊人婷婷丁香| 国产aⅴ精品一区二区三区波| 亚洲国产欧洲综合997久久,| netflix在线观看网站| 欧美成人性av电影在线观看| 午夜两性在线视频| 欧美xxxx性猛交bbbb| 久久国产乱子伦精品免费另类| 99热这里只有精品一区| 国产成年人精品一区二区| 波多野结衣巨乳人妻| 国产伦人伦偷精品视频| 国产亚洲精品久久久久久毛片| 在现免费观看毛片| 亚洲乱码一区二区免费版| 99热这里只有精品一区| 亚洲熟妇中文字幕五十中出| 国产熟女xx| 夜夜爽天天搞| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人精品二区| 久久精品国产亚洲av香蕉五月| 精品久久久久久久久久免费视频| 国产一区二区三区在线臀色熟女| 成年版毛片免费区| 真人一进一出gif抽搐免费| 男人舔奶头视频| 精品日产1卡2卡| 亚洲av.av天堂| 国产激情偷乱视频一区二区| 国产伦人伦偷精品视频| 全区人妻精品视频| 免费在线观看成人毛片| 最新中文字幕久久久久| 久久久精品大字幕| 亚洲av一区综合| 亚洲精品一卡2卡三卡4卡5卡| 欧美色视频一区免费| 人人妻人人澡欧美一区二区| 欧美高清性xxxxhd video| 男女床上黄色一级片免费看| 亚洲午夜理论影院| 午夜福利在线观看吧| 成人毛片a级毛片在线播放| 国产单亲对白刺激| 精品乱码久久久久久99久播| 国产精品不卡视频一区二区 | 国产高清视频在线播放一区| 一级av片app| 色吧在线观看| 色综合婷婷激情| 免费搜索国产男女视频| 欧美一区二区精品小视频在线| 男人狂女人下面高潮的视频| 在线免费观看的www视频| 午夜日韩欧美国产| 最近最新免费中文字幕在线| 国产一级毛片七仙女欲春2| 亚洲成a人片在线一区二区| 亚洲美女视频黄频| 看片在线看免费视频| 欧美zozozo另类| 麻豆久久精品国产亚洲av| 99热这里只有是精品在线观看 | 午夜视频国产福利| 欧美zozozo另类| 国产免费一级a男人的天堂| 色视频www国产| 91麻豆精品激情在线观看国产| 中文字幕高清在线视频| 欧美在线一区亚洲| 人妻夜夜爽99麻豆av| 男人狂女人下面高潮的视频| 亚洲自偷自拍三级| 国产伦在线观看视频一区| 最后的刺客免费高清国语| 国产精品影院久久| 国产在视频线在精品| 欧美黄色片欧美黄色片| 91九色精品人成在线观看| 亚洲精品久久国产高清桃花| 亚洲精品在线观看二区| 黄色丝袜av网址大全| 久久久久久大精品| 麻豆成人av在线观看| 在线免费观看不下载黄p国产 | a在线观看视频网站| 神马国产精品三级电影在线观看| 色尼玛亚洲综合影院| 亚洲av日韩精品久久久久久密| 国产白丝娇喘喷水9色精品| 熟女人妻精品中文字幕| 久久这里只有精品中国| 亚洲中文日韩欧美视频| 97热精品久久久久久| 亚洲 国产 在线| 日韩欧美在线二视频| 欧美日韩中文字幕国产精品一区二区三区| 一个人看的www免费观看视频| 一进一出抽搐动态| 国产免费av片在线观看野外av| 人人妻人人看人人澡| 午夜福利在线观看吧| 人人妻,人人澡人人爽秒播| 国产在视频线在精品| 极品教师在线免费播放| 国产精品亚洲一级av第二区| 国产69精品久久久久777片| 草草在线视频免费看| 午夜视频国产福利| 久久久久国内视频| 精品久久久久久,| 男女那种视频在线观看| 精品免费久久久久久久清纯| 麻豆久久精品国产亚洲av| 成人亚洲精品av一区二区| 国产一区二区亚洲精品在线观看| 黄色日韩在线| 国产视频内射| 国产高清视频在线观看网站| 国产探花极品一区二区| 99精品久久久久人妻精品| 国产高清激情床上av| 午夜日韩欧美国产| 男人的好看免费观看在线视频| 精品人妻视频免费看| 精品国内亚洲2022精品成人| 欧美zozozo另类| 欧美+亚洲+日韩+国产| 午夜久久久久精精品| 欧美日韩乱码在线| 国产高清视频在线播放一区| 别揉我奶头 嗯啊视频| 极品教师在线免费播放| 国产一区二区三区视频了| 色综合欧美亚洲国产小说| 色综合站精品国产| 全区人妻精品视频| 如何舔出高潮| 欧美性猛交黑人性爽| 全区人妻精品视频| 欧美成人免费av一区二区三区| 有码 亚洲区| 又紧又爽又黄一区二区| 女人被狂操c到高潮| 哪里可以看免费的av片| 成年版毛片免费区| 在线看三级毛片| 99热只有精品国产| 婷婷亚洲欧美| 欧美色欧美亚洲另类二区| 亚洲自偷自拍三级| 国产白丝娇喘喷水9色精品| 国产日本99.免费观看| 五月玫瑰六月丁香| 亚洲国产精品合色在线| x7x7x7水蜜桃| www.999成人在线观看| 亚洲av二区三区四区| 国产日本99.免费观看| 亚洲美女视频黄频| 国产精品电影一区二区三区| 我要搜黄色片| 中文字幕av成人在线电影| 观看免费一级毛片| 丝袜美腿在线中文| 2021天堂中文幕一二区在线观| 亚洲欧美日韩高清在线视频| 嫩草影院入口| 亚洲精品日韩av片在线观看| 亚洲av中文字字幕乱码综合| 午夜福利视频1000在线观看|