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

    田間原位試驗(yàn)分析長(zhǎng)期機(jī)械作業(yè)下稻麥輪作地塊土壤入滲性能

    2023-06-12 03:25:14孫浩田徐高明李毅念何瑞銀丁啟朔
    關(guān)鍵詞:實(shí)度耕作層底層

    梁 磊,孫浩田,徐高明,李毅念,何瑞銀,丁啟朔

    ?農(nóng)業(yè)水土工程?

    田間原位試驗(yàn)分析長(zhǎng)期機(jī)械作業(yè)下稻麥輪作地塊土壤入滲性能

    梁 磊,孫浩田,徐高明,李毅念,何瑞銀,丁啟朔※

    (1. 南京農(nóng)業(yè)大學(xué)工學(xué)院,南京 10031; 2. 江蘇省智能化農(nóng)業(yè)裝備重點(diǎn)實(shí)驗(yàn)室,南京 10031)

    田間原位不同深度入滲試驗(yàn)是表達(dá)土壤分層狀態(tài)、展示土層物理分異以及定量土壤剖面水功能變化的關(guān)鍵。為了探究不同深度水稻土的入滲能力及保水作用,該研究以華東稻麥輪作區(qū)小農(nóng)戶長(zhǎng)期機(jī)械化耕整模式下代表地塊的土層分異為目標(biāo),設(shè)計(jì)田間原位不同深度入滲試驗(yàn)。在試驗(yàn)地塊內(nèi)開挖7個(gè)不同深度的入滲坑并在坑底進(jìn)行入滲試驗(yàn),然后滲透48 h分層測(cè)取土壤含水率,研究不同坑底深度(坑深)土壤的入滲能力和入滲后各土層含水率的變化。結(jié)果表明,不同深度入滲試驗(yàn)準(zhǔn)確表達(dá)了不同坑深土壤的水分入滲及土層持水分異,同時(shí)也能清晰地鑒別出犁底層所在位置和厚度,犁底層始于15 cm深,且耕作層與犁底層分異明顯,耕作層平均緊實(shí)度為1 005.79 kPa,犁底層平均緊實(shí)度為1 910.73 kPa;土壤剖面分析表明,耕作層土壤形態(tài)疏松,根系分布稠密,犁底層土壤容重大,孔隙度小,透水性差,心土層土壤鐵錳斑點(diǎn)較多,結(jié)構(gòu)性差;土壤入滲參數(shù)隨坑深的增加而減少,其中0~15 cm坑深范圍內(nèi)平均的平均入滲速率和累計(jì)入滲量分別為>20~30 cm的17.04倍和18.06倍;通過對(duì)比初始含水率和滲透48 h后含水率,得到坑深在15 cm范圍內(nèi)的水分入滲深度均為距表土20 cm,而大于15 cm坑深的水分入滲深度均為坑深以下10 cm;利用Horton、Kostiakov和Philip 3個(gè)模型擬合不同坑深入滲,結(jié)果表明,Kostiakov模型的擬合參數(shù)與實(shí)際相符且2最高(0.98~0.99),均方根誤差最小(0.01~0.77 mm/min)。入滲參數(shù)與土壤容重、含水率、總孔隙度和田間持水量有極顯著的相關(guān)性(<0.01),與土壤緊實(shí)度相關(guān)性不顯著(>0.05)。該研究通過原位不同深度入滲試驗(yàn),說明華東稻麥輪作區(qū)小農(nóng)戶生產(chǎn)模式下長(zhǎng)期機(jī)械化耕作導(dǎo)致水稻土明顯的分層和土壤水功能垂直分異,進(jìn)而導(dǎo)致耕作層與耕作層以下土層表現(xiàn)出顯著的入滲能力差異。研究可為稻麥輪作地區(qū)的機(jī)械耕作、灌溉等提供借鑒。

    土壤;入滲;水分;稻麥輪作區(qū);耕作層;犁底層

    0 引 言

    長(zhǎng)江中下游稻麥輪作區(qū)土壤周年干濕交替,土壤結(jié)構(gòu)保護(hù)及耕作層水管理也隨水-旱作物的輪茬交替而周年切換。集約化生產(chǎn)模式及機(jī)械化的發(fā)展導(dǎo)致小農(nóng)戶生產(chǎn)方式下普遍以旋代耕,加之周年多次田間機(jī)械作業(yè)造成土壤壓實(shí),多重因素導(dǎo)致耕層淺薄、土壤緊實(shí)、孔隙破壞,進(jìn)而增加作物產(chǎn)量降低的風(fēng)險(xiǎn)[1-2]。集約化稻麥輪作系統(tǒng)的可持續(xù)性必然要求增加耕作層厚度、下移犁底層,而深松是目前普遍應(yīng)用的保護(hù)性耕作措施之一[3]。

    為打破犁底層、增加土壤保水量,北方旱作制深松耕深通常達(dá)到30 cm[3-4]。對(duì)稻麥輪作制而言,如果在旱作的麥季實(shí)施深耕深松創(chuàng)建厚度達(dá)30 cm的耕層,雖可打破犁底層提高土壤蓄水,但此做法在生產(chǎn)中常面臨諸多質(zhì)疑,人們擔(dān)心麥季破壞的犁底層失去保水阻滲的效果,輪作到稻季時(shí)水分下滲加劇、形成漏斗田、增加稻季用水。稻麥輪作周年混亂的農(nóng)藝現(xiàn)狀及尚存爭(zhēng)議的耕作措施要求基于生產(chǎn)實(shí)際闡明田間土壤的物理性質(zhì)及耕作深度對(duì)土壤入滲能力的影響[5-6]。同時(shí)基于土壤精細(xì)分層并逐層定量土壤入滲能力及土層間差異性的研究尚未見諸報(bào)道。

    近年來耕作與土壤入滲關(guān)系的研究越發(fā)受到關(guān)注,侯賢清等[7]研究表明,2 a免耕+1 a深松和2 a深松+1 a免耕能顯著提高土壤水分入滲,降低土壤體積質(zhì)量和增大土壤孔隙度。張西平等[8]研究表明耕作對(duì)土壤入滲參數(shù)的影響比種植方式大。姚毓香等[9]針對(duì)深松效應(yīng)研究時(shí)發(fā)現(xiàn)適當(dāng)增加鏟間距有利于加速土壤水分入滲。安晶等[10]研究表明深松有利于創(chuàng)造上虛下實(shí)的土壤結(jié)構(gòu),提高上層土壤的保水和蓄水能力。長(zhǎng)期機(jī)械化生產(chǎn)導(dǎo)致農(nóng)田土壤分層,土壤剖面從上至下分異成為耕作層、犁底層和心土層[4]?;趯?duì)土層分異性的關(guān)注,孫蓓等[11]對(duì)田間耕作層和犁底層土壤分別取樣并開展室內(nèi)模擬試驗(yàn),據(jù)此提出可連續(xù)測(cè)量耕作層與犁底層入滲過程的方法。劉淙琮等[12]獲取田間耕作層和犁底層土壤進(jìn)行室內(nèi)土柱模擬入滲試驗(yàn),表明犁底層顯著提高上層土壤的含水量,并且犁底層深度的增加可有效促進(jìn)水分在上層土壤的儲(chǔ)存。

    已報(bào)道的耕作與土壤入滲關(guān)系的研究主要集中在耕后土壤檢測(cè),展示耕后土壤剖面的入滲能力變化,難以指導(dǎo)稻麥輪作地區(qū)的耕作深度,同時(shí)對(duì)不同深度的入滲過程鮮有研究,更未見到田間原位開挖不同深度的坑底入滲試驗(yàn)及相關(guān)報(bào)道。對(duì)長(zhǎng)江中下游區(qū)域的稻麥輪作制而言,尤其缺乏針對(duì)田間水稻土不同開挖深度坑底入滲試驗(yàn)所得的不同土層入滲能力的數(shù)據(jù)。因此有必要設(shè)計(jì)一種基于不同深度處理的田間原位入滲試驗(yàn),準(zhǔn)確定量土壤剖面精細(xì)分層的各土層入滲及蓄水能力。本文設(shè)計(jì)了一種精準(zhǔn)控制田間開挖深度的坑底入滲試驗(yàn),用于準(zhǔn)確測(cè)定不同深度土層的入滲能力及入滲后各土層的含水率變化,闡明長(zhǎng)期集約化機(jī)械化稻麥輪作對(duì)土壤分層及土壤水行為的影響,以期為稻麥輪作地區(qū)機(jī)械耕作、灌溉等提供借鑒。

    1 材料與方法

    1.1 研究地概況

    田間試驗(yàn)位于南京市六合區(qū)八百橋試驗(yàn)田(118°55′E,32°25′N),平均海拔22.7 m,地勢(shì)平坦,土壤質(zhì)地為壤質(zhì)黏土,土壤是長(zhǎng)期(10 a以上)稻麥輪作而形成的水稻土。該區(qū)處于亞熱帶季風(fēng)氣候帶,年平均溫度為15.7 ℃,年平均降水量為1 000 mm左右,該地區(qū)生產(chǎn)周年全程機(jī)械化。試驗(yàn)地為小農(nóng)戶土地經(jīng)營(yíng)模式,地塊分散、細(xì)碎化嚴(yán)重,周年(10 a以上)機(jī)械化作業(yè)環(huán)節(jié)繁多且長(zhǎng)期以旋代耕,耕層淺薄、土層分異明顯。土壤入滲試驗(yàn)進(jìn)行于2021年11月水稻收獲后,試驗(yàn)前人工清除地表留茬。

    1.2 田間試驗(yàn)方法

    1.2.1 土壤剖面入滲試驗(yàn)設(shè)計(jì)

    耕作后受擾動(dòng)的土壤產(chǎn)生較大的孔隙,大的孔隙不具有保水阻滲作用,為了探究深松耕深以下未受擾動(dòng)土壤的入滲能力及保水作用,本文參照常規(guī)深松的最大耕深(圖1)設(shè)置5、10、15、20、25和30 cm共6個(gè)不同開挖深度,挖坑截面為60 cm×60 cm,移除坑內(nèi)土壤并以坑底作為入滲面進(jìn)行入滲試驗(yàn)。不同深度的入滲坑及坑底入滲能夠直接準(zhǔn)確測(cè)定各土層的入滲能力并區(qū)分土層間的物理分異。試驗(yàn)采用(0 cm)表土入滲作為對(duì)照,因此共計(jì)7個(gè)不同的入滲試驗(yàn),各處理3次重復(fù)。

    圖1 土壤分層以及深松耕深示意圖

    1.2.2 土壤剖面入滲試驗(yàn)過程

    試驗(yàn)采用鐵框法[13],鐵框截面為正方形(圖2a),分為大框(邊長(zhǎng)50 cm、高20 cm)和小框(邊長(zhǎng)25 cm、高15 cm)。將土壤表面輕輕鏟平,將大小框按照?qǐng)D2b所示切入土中3 cm,且相鄰大框的距離大于50 cm。為避免水流過大沖擊表層土壤,用水桶沿著框壁四周對(duì)大小框同時(shí)注水,并使大小框內(nèi)水位均為5 cm,以便大小框內(nèi)的水垂直入滲。注水后,立即用秒表計(jì)時(shí),然后用量杯持續(xù)向框內(nèi)外加水維持水位5 cm。因入滲初期土壤基質(zhì)吸力和重力的作用,入滲速率先快后慢,所以記錄時(shí)間間隔由快變慢,分別在0、0.5、1、2、3、4、5、6、7、8、9、10 min及10 min后每間隔10 min記錄,1 h后間隔30 min記錄[14],直到穩(wěn)滲狀態(tài)停止試驗(yàn)。入滲總歷時(shí)120 min,初始入滲速率為0~1 min的入滲速率[15-16],穩(wěn)定入滲速率為單位時(shí)間內(nèi)入滲速率趨于同一入滲速率,本研究在90 min時(shí)達(dá)到穩(wěn)定入滲,平均入滲速率表示到達(dá)穩(wěn)定入滲狀態(tài)后,累計(jì)入滲量與到達(dá)穩(wěn)滲時(shí)間的比值,取前90 min內(nèi)的入滲總量與時(shí)間的比值。因本文入滲總時(shí)間為120 min,故取120 min的滲透總量作為累計(jì)入滲量進(jìn)行分析。入滲速率為

    式中V為第次測(cè)量時(shí)的入滲速率,mm/min;Q為第次測(cè)量時(shí)量筒加入的水量,mL;為小鐵框的橫截面積,cm2;T為第次測(cè)量的時(shí)間間隔,min。

    第次測(cè)量時(shí)的累計(jì)入滲量為

    式中為總加水次數(shù)。

    圖2 鐵框法測(cè)定田間土壤水分入滲過程

    1.3 指標(biāo)測(cè)定與計(jì)算方法

    1.3.1 土壤物理參數(shù)

    入滲試驗(yàn)前,在鐵框旁隨機(jī)選取21個(gè)點(diǎn)進(jìn)行緊實(shí)度測(cè)試,使用土壤緊實(shí)度儀(TJSD-750-II型,浙江托普儀器有限公司)測(cè)取,每隔2.5 cm深度記錄一次,直到深度40 cm。在緊實(shí)度測(cè)試點(diǎn)旁取土壤容重,用環(huán)刀(高5 cm,體積100 cm3)采集0~5、>5~10、>10~15、>15~20、>20~25、>25~30、>30~35、>35~40、>40~45、>45~50 cm 10個(gè)土層的土樣,每層3次重復(fù),取出的土樣立即稱取質(zhì)量獲取濕土質(zhì)量,用來測(cè)定各土層的初始含水率。同時(shí)在較為平坦的地面上挖出深50 cm、寬50 cm的地坑用以觀察水稻土剖面。

    入滲試驗(yàn)結(jié)束后,為了闡明不同深度坑底入滲的滲透深度,依據(jù)《麥秸稈機(jī)械化全量還田與水稻機(jī)插秧集成技術(shù)》[17],繼續(xù)維持入滲框內(nèi)水位滲透48 h,環(huán)刀法分層取出入滲面下部土樣,每層3次重復(fù),測(cè)定各土層滲泡后的含水率,用來檢驗(yàn)初始含水率和滲泡后含水率的差異顯著性并基于此判斷入滲到達(dá)的土層深度。

    將土樣帶回實(shí)驗(yàn)室置于105 ℃條件下烘干48 h,稱取干土質(zhì)量,然后其與環(huán)刀體積的比值即為土壤容重,相應(yīng)的土壤總孔隙度計(jì)算式為

    土壤含水率為

    式中為含水率,%;為濕土質(zhì)量,g;m為干土質(zhì)量,g。

    田間持水量是指排水良好的土地上充分灌溉或者降水后,水分充分下滲,經(jīng)過一段時(shí)間,多余的重力水排走后,土壤所能維持的最大含水量,是大多數(shù)植物可利用水的上限[20]。將土樣放置在充滿水的平底金屬托盤中,靜置24 h,土樣飽和后,將其靜置在干燥的土壤上8 h排除重力水,然后稱質(zhì)量,按式(5)計(jì)算田間持水量:

    式中w為田間持水量,%;m為飽和濕土質(zhì)量,g。

    1.3.2 土壤入滲模型

    土壤入滲模型可以反映土壤入滲性能,利用相關(guān)數(shù)學(xué)模型對(duì)入滲速率擬合可以更直觀描述不同土層的入滲規(guī)律。本研究分別選取Horton、Koistakov、Philip 入滲模型[21-24]對(duì)土壤入滲速率與入滲時(shí)間進(jìn)行擬合,各模型計(jì)算式及參數(shù)見表1。

    表1 土壤入滲模型及參數(shù)說明

    1.4 數(shù)據(jù)處理

    試驗(yàn)數(shù)據(jù)使用Excel 2019進(jìn)行基本處理,利用SPSS 26.0對(duì)土壤物理和入滲參數(shù)進(jìn)行方差分析(one-way ANOVA,顯著性水平為0.05)。圖和表中的數(shù)據(jù)均為平均值,圖中的誤差線為標(biāo)準(zhǔn)差。使用Pearson相關(guān)性分析表示土壤物理和入滲參數(shù)之間的相關(guān)關(guān)系,使用Origin 2020軟件繪圖,利用Matlab R2020軟件進(jìn)行土壤入滲模型擬合。

    2 結(jié)果與分析

    2.1 土壤物理性質(zhì)

    土壤容重變化是反映土層物理分異的基礎(chǔ)指標(biāo)之一。表2顯示0~50 cm深度范圍內(nèi)土壤容重以15 cm為分界,呈現(xiàn)上下土層間的顯著性差異(<0.05)。除此之外,土壤含水率也反映出各土層的土壤物理分異性,土壤含水率在0~15 cm范圍內(nèi)隨著土層深度增加呈現(xiàn)遞減變化趨勢(shì),同時(shí)顯示15 cm以上與15 cm以下顯著變化。田間持水量在0~20 cm內(nèi)每5 cm顯著降低,之后不再顯著變化,總孔隙度在0~15 cm內(nèi)每5 cm變化不顯著,而15 cm以上與15 cm以下顯著變化,表現(xiàn)為以15 cm深度為分界,呈上下顯著性差異(<0.05)。可見,土壤分層取樣分析及土壤層物理分異性表達(dá)涉及的土壤容重、含水率、田間持水量、總孔隙度等多個(gè)土壤物理基礎(chǔ)指標(biāo)反映的土層物理分異特征各不相同。

    表2 0~50 cm土壤物理參數(shù)

    注:同列不同的小寫字母表示在0.05水平上土層間存在顯著性差異。下同。

    Note: Different lowercase letters in same column indicate significant differences among soil layers at the 0.05 level. Same as below.

    土壤緊實(shí)度是反映土壤抵抗外力壓實(shí)以及破碎的一種能力,直接影響耕作阻力以及作物根系的下扎[25],通過測(cè)量土壤緊實(shí)度能夠確定犁底層在耕地中的深度[26]。由圖3可知,土壤緊實(shí)度隨土壤深度呈“幾”型的變化趨勢(shì)。圖3緊實(shí)度結(jié)合表2土壤物理參數(shù)可知,耕作層位于0~15 cm,平均緊實(shí)度為1 005.79 kPa,犁底層處于15~20 cm處,平均緊實(shí)度為1 910.73 kPa,心土層處于20 cm以下,平均緊實(shí)度為1 412.22 kPa,犁底層緊實(shí)度約為耕作層緊實(shí)度的2倍。犁底層的致密性影響植物根系的延伸生長(zhǎng)和對(duì)下層養(yǎng)分的吸收,還影響氣體和熱量的上下交換,更是影響農(nóng)田的水循環(huán)[26-27]。同時(shí)犁底層具有避免其以下土壤壓實(shí)的風(fēng)險(xiǎn),有研究表明在收割機(jī)軋后15 cm內(nèi)土壤緊實(shí)度顯著變化,但不造成15 cm以下土壤顯著破壞[19]。

    圖3 0~40 cm土壤緊實(shí)度

    2.2 土壤剖面結(jié)構(gòu)特征

    土壤剖面結(jié)構(gòu)特征對(duì)土壤水分交換、降雨入滲、地下水補(bǔ)給及物質(zhì)交換等有重要影響[28]。圖4可以看出,耕作層土壤較為疏松,分布著較多根系,土壤顆粒間空隙較大,土壤以團(tuán)聚體形式為主,總體結(jié)構(gòu)性較好,利于水分入滲,是機(jī)械耕作的主要影響區(qū)域。犁底層土壤較為緊實(shí),根系下扎困難,土壤孔隙度小,土壤以大塊狀分布,透水性差,結(jié)合表2可知,犁底層對(duì)比耕作層土壤容重增大,含水率減小。心土層處于犁底層之下,與上層土壤物質(zhì)交換較少,土壤堅(jiān)實(shí),鐵錳斑點(diǎn)較多,結(jié)構(gòu)性差,養(yǎng)分含量低。結(jié)合表2和圖3可知,心土層緊實(shí)度比犁底層小,容重、含水率、田間持水量和總孔隙度并無顯著性差異。

    圖4 水稻土剖面結(jié)構(gòu)特征

    2.3 不同坑深的土壤入滲過程

    由圖5a可知,不同坑深的土壤在0~5 min內(nèi)水分入滲較快,處于瞬變階段。在6~90 min處于滲漏階段,此階段隨著入滲時(shí)間的增加土壤入滲速率逐漸減小且變化趨于平緩。在90 min后,水的入滲達(dá)到穩(wěn)滲階段,此階段土壤入滲速率變化不明顯,接近穩(wěn)定狀態(tài),并且隨著坑深的增加,穩(wěn)滲速率減小。從圖中可以看出大于15 cm坑深的土壤水分入滲速率相差較小。由圖5b可知,在入滲前90 min,由于前期入滲較快,累計(jì)入滲量呈對(duì)數(shù)型增長(zhǎng),到達(dá)90 min后,累計(jì)入滲量不再快速增長(zhǎng),保持穩(wěn)定。與入滲速率相同,隨著坑深的增加,累計(jì)入滲量在減小,且大于15 cm坑深的累計(jì)入滲量相差較小,這可能是由于15 cm以下土壤的土壤緊實(shí)度均大于表層土層,此處土壤常年未受到擾動(dòng),土壤較為濕黏,孔隙差。

    圖5 不同坑深的土壤入滲過程

    土壤初始入滲速率、穩(wěn)定入滲速率和平均入滲速率是反映土壤入滲能力的重要指標(biāo),比較其大小關(guān)系可以了解不同土層入滲能力的強(qiáng)弱[29]。從表3可知,不同坑深的土壤初始入滲速率和穩(wěn)定入滲速率呈現(xiàn)顯著性差異(<0.05)。各個(gè)坑深的土壤初始入滲速率、穩(wěn)定入滲速率、平均入滲速率和累計(jì)入滲量均表現(xiàn)為隨坑深遞增,數(shù)值遞減,其中0~15 cm坑深范圍內(nèi)平均的平均入滲速率和累計(jì)入滲量分別為>20~30 cm的17.04倍和18.06倍,耕作層與耕作層以下的土壤入滲能力差異較大。

    表3 不同坑深的土壤入滲參數(shù)

    注:不同小寫字母表示不同坑深間差異顯著(<0.05)。

    Note: Different lowercase letters indicate significant differences among pit depths at the 0.05 level.

    2.4 不同坑深滲后土壤含水率

    同一坑深條件下的滲后不同土層含水率變化的顯著性是反映土層間含水率差異的關(guān)鍵指標(biāo),以表4的縱向數(shù)據(jù)間差異顯著性檢驗(yàn)為判斷依據(jù)。不同坑深下滲后同一土層間的含水率變化顯著性是反映水分從開始入滲到“截止”入滲的深度的關(guān)鍵指標(biāo),以表4的橫向數(shù)據(jù)間差異顯著性檢驗(yàn)為判斷依據(jù)。由表4可知,田間持水量和初始含水率均以>15~20 cm為分界呈現(xiàn)上層與下層土壤間的顯著性差異(<0.05),這緊密對(duì)應(yīng)前文所述的犁底層鑒定結(jié)果。經(jīng)過48 h滲透后,對(duì)于0、5和10 cm坑深的滲后各土層而言,以>20~25 cm為分界呈現(xiàn)上層與下層含水率的顯著性變化(<0.05)。同樣的>25~30 cm為分界也發(fā)生在15 cm坑深的入滲處理。對(duì)于20 cm坑深則表現(xiàn)出以>30~35 cm土層為分界。對(duì)于25 cm坑深的入滲處理,>35~40 cm土層是其分界。30 cm的坑深入滲處理則以>40~45 cm為分界。

    表4 不同坑深滲后田間持水量與含水率變化

    注:大寫字母表示不同土層差異顯著(<0.05),小寫字母表示同一土層不同坑深差異顯著(<0.05)。

    Note: Uppercase letters indicate significance at different soil layers (<0.05), lowercase letters indicate significance at different pit depths at the same layer (<0.05).

    由表4可知,對(duì)于0、5和10 cm坑底入滲處理,滲后土層含水率與初始含水率在>20~25 cm土層間的差異性不顯著,表明0、5和10 cm坑深的水分入滲達(dá)20 cm。同理,對(duì)于15 cm坑深,>25~30 cm土層間的差異性不顯著(>0.05)。對(duì)于20 cm坑深,>30~35 cm土層間的差異性不顯著(>0.05)。對(duì)于25 cm坑深,>35~40 cm土層間的差異性不顯著(>0.05)。對(duì)于30 cm坑深,>40~45 cm土層間的差異性不顯著(>0.05)。

    將表4不同坑深滲后土層含水率與對(duì)應(yīng)土層的田間持水量進(jìn)行橫向比較發(fā)現(xiàn),滲透48 h后,坑深以下5 cm土壤均含有重力水,處于過飽和狀態(tài)。通過縱向與橫向?qū)Ρ瓤芍?,滲透48 h后,坑深小于15 cm(等效于淺旋耕未打破犁底層),水分入滲深度可達(dá)20 cm,而坑深大于或等于15 cm時(shí),水分入滲深度為坑深以下10 cm。鑒于常規(guī)旋耕和水田打漿深度在10~15 cm[30-31]并未打破犁底層,這印證了生產(chǎn)中人們對(duì)水田犁底層具有保水保肥功能的基本認(rèn)知的正確性。

    2.5 土壤入滲模型模擬結(jié)果分析

    為了進(jìn)一步探究不同坑深土壤水分入滲速率與時(shí)間的變化關(guān)系以及入滲模型在水稻土的適用性,采用Horton、Kostiakov和Philip 3個(gè)常用模型擬合入滲速率,擬合結(jié)果如表5所示。決定系數(shù)2的取值范圍為[0,1],2越接近1,表明擬合程度越佳。均方根誤差(root mean squared error,RMSE)越趨近于0表明模型的擬合度越好。Horton 模型的2為0.91~0.97,RMSE平均值為0.02~1.63 mm/min,表明該模型對(duì)不同坑深土壤入滲的過程擬合程度較好,但與表3相比,模型擬合的c過高,不符合真實(shí)穩(wěn)定入滲速率,所以該模型不適用擬合本研究土壤水分入滲速率。Kostiakov 模型的2為0.98~0.99,RMSE為0.01~0.77 mm/min,表明該模型對(duì)于土壤入滲速率擬合度良好。Philip 模型的2為0.95~0.99,RMSE為0.03~0.80 mm/min,表明該模型對(duì)不同坑深的入滲過程擬合程度較好,表征土壤入滲能力的吸滲率隨著坑深增加土壤入滲能力越弱,這與實(shí)際土壤入滲能力相一致,但模型擬合的c為負(fù)值,與實(shí)際不符,說明該模型不適用于本研究。綜上,Kostiakov 模型最適合于分析水稻土不同坑深的水分入滲過程。在Kostiakov模型中,參數(shù)表示土壤入滲的第1時(shí)段平均入滲速率,第1時(shí)段平均入滲速率越大土壤入滲性能越好,參數(shù)為衰減系數(shù),參數(shù)越大表示土壤入滲速率隨時(shí)間增加遞減越快[32]。從表5可知,參數(shù)隨著坑深的增加而減小,與表3的初始入滲速率規(guī)律相同,參數(shù)隨著坑深的增加并無明顯差異。

    表5 不同坑深的土壤入滲模型擬合參數(shù)

    注:2為決定系數(shù);RMSE為均方根誤差。

    Note:2is dertermination coefficient;RMSEis root mean squared error.

    2.6 入滲參數(shù)與土壤物理相關(guān)性分析

    由表4可知,當(dāng)坑深小于15 cm入滲時(shí),水分入滲可達(dá)20 cm土壤層,當(dāng)坑深大于或等于15 cm開始入滲時(shí),水分入滲深度為坑深以下10 cm。因此,為了探究入滲參數(shù)與土壤初始物理參數(shù)之間的關(guān)系,對(duì)0、5、10 cm坑深的土壤物理參數(shù)取20 cm內(nèi)的平均土壤初始物理參數(shù),對(duì)15、20、25、30 cm坑深的土壤物理參數(shù)取其坑深以下10 cm內(nèi)的平均土壤初始物理參數(shù),對(duì)土壤初始物理參數(shù)與入滲參數(shù)進(jìn)行Pearson相關(guān)性分析如表6所示。土壤入滲參數(shù)與容重呈極顯著負(fù)相關(guān)(<0.01),與含水率和田間持水量呈極顯著正相關(guān)(<0.01),與土壤緊實(shí)度不顯著(>0.05)。這表明田間土壤容重越大,入滲性能越差。劉姍姍等[33]研究表明入滲參數(shù)與土壤初始含水率、土壤緊實(shí)度和黏粒含量均呈負(fù)相關(guān),曾辰等[34]研究表明累積入滲量隨初始含水率的增大而減小,然而本文含水率越高,入滲性能反而越好,這可能是由于表層土壤含水率相對(duì)于底層含水率高的原因,突顯了田間原位土壤入滲的特殊性。對(duì)于田間耕作,耕作層越疏松,越有利于水分的入滲,然而卻不利于作物根系的下扎,如何平衡好土壤物理參數(shù),對(duì)土壤水循環(huán)及耕作具有重要意義。

    表6 入滲參數(shù)與土壤初始物理參數(shù)相關(guān)關(guān)系

    注:**表示極顯著(<0.01)。

    Note: ** mean extreme significance (<0.01).

    3 討 論

    農(nóng)田耕作層與犁底層土壤性質(zhì)不盡相同,上層耕作層土壤結(jié)構(gòu)良好、結(jié)構(gòu)孔隙豐富,而下層犁底層土壤致密,多表現(xiàn)為基粒材質(zhì)孔隙,因此不同土層的水分入滲有著不同表現(xiàn)。本研究表明,試驗(yàn)地水稻土耕作層位于15 cm以上,土壤較為疏松,分布著較多根系,犁底層處于15 cm以下處,土壤較為緊實(shí),根系較少,這些基于精細(xì)土層控制的坑深入滲所得結(jié)論與HUO等[19,35]研究結(jié)果一致,從土壤水的角度進(jìn)一步印證以往區(qū)域調(diào)查研究所得土層結(jié)論,同時(shí)本研究進(jìn)一步發(fā)現(xiàn)犁底層平均緊實(shí)度約為耕作層的2倍。

    土壤容重、緊實(shí)度、含水率等指標(biāo)的土層分異也印證土壤剖面視覺評(píng)價(jià)方法的可行性,疏松多孔的土壤層有利于作物根系的生長(zhǎng)發(fā)育,緊實(shí)少孔的土壤不利于作物根系下扎,土壤剖面配合土壤結(jié)構(gòu)體視覺評(píng)價(jià)方法是生產(chǎn)中簡(jiǎn)單、可靠、實(shí)用的田間原位觀測(cè)評(píng)價(jià)方法[36]。耕作層厚度與本地長(zhǎng)期的耕作模式有很大的關(guān)系,本文地區(qū)長(zhǎng)期10~15 cm深度的旋耕和水田打漿[30-31]決定了耕作層的厚度為15 cm,而對(duì)于北方旱作土壤,犁底層出現(xiàn)的位置為15~30 cm[37-38],水稻土相比北方旱作土壤,耕作層偏薄,犁底層出現(xiàn)位置偏上(>15~20 cm),不利于作物根系下扎,影響作物最終產(chǎn)量的提高。因而,針對(duì)稻麥輪作區(qū)的水稻土,耕作層淺薄,土壤板結(jié),孔隙度低等問題[1-2],有必要進(jìn)行深松耕作,使犁底層下移,消除土壤板結(jié),改善土壤結(jié)構(gòu),提高土壤入滲能力[39]。

    土壤入滲能力除了受水頭、溫度等外在因素的影響,也受土壤自身物理特性(尤其是結(jié)構(gòu)孔隙)的影響,如質(zhì)地、容重、含水率、緊實(shí)度、孔隙度等[24]。田間土壤隨著土層深度的增加,各層土壤物理性質(zhì)表現(xiàn)出差異性[40],土壤物理性質(zhì)的不同就會(huì)引起土壤入滲能力的不同。從試驗(yàn)結(jié)果可知,初始入滲速率、穩(wěn)定入滲速率和平均入滲速率隨著土壤深度的增加而減小,同余林等[41]研究結(jié)論一致。當(dāng)滲透48 h后,小于15 cm坑深入滲時(shí),水稻土入滲深度均為距表土20 cm,而包含等[42]通過室內(nèi)模擬降雨試驗(yàn)表明砂壤土的最佳保水性土層深度為20~35 cm。王曉彤等[43]研究表明層狀土壤黏土夾層能夠明顯改變水分的入滲特性。肖慶禮等[44]表明土壤下部黏土夾層不僅對(duì)水分入滲起到顯著的阻礙作用,而且具有較強(qiáng)的持水保水能力。這些研究表明層狀土壤中存在土層保水阻滲,水稻土的犁底層具有土壤孔隙小、緊實(shí)和黏重的特點(diǎn)[26-27],這些特點(diǎn)起著儲(chǔ)存水分和保水阻滲的作用,同時(shí)YI等[45]研究表明降低犁底層的飽和導(dǎo)水率可以阻礙水滲透,且隨著水稻栽培年限的增加,犁底層的厚度和容重在增加。3個(gè)常用模型對(duì)不同坑深土壤水分入滲過程擬合優(yōu)度從大到小為:Kostiakov模型、Philip模型、Horton模型,結(jié)果表明Kostiakov模型對(duì)于水稻土的適用性比較高,這與羅瑩麗等[22-23,46]在非水稻土水分入滲模型的研究中得出結(jié)果相似。

    4 結(jié) 論

    為了獲取不同深度水稻土的入滲能力及保水作用,用以指導(dǎo)深松耕深,本文在長(zhǎng)期機(jī)械作業(yè)的稻麥輪作地塊內(nèi)開挖7個(gè)不同深度的入滲坑并在坑底進(jìn)行入滲試驗(yàn),研究了不同坑深土壤的入滲能力和入滲后各土層含水率的變化,得出如下結(jié)論:

    1)通過分析田間水稻土的物理參數(shù)及剖面得出,耕作層位于0~15 cm土深,分布較多根系,平均緊實(shí)度為1 005.79 kPa,犁底層處于>15~20 cm土深,容重大,孔隙度小,平均緊實(shí)度為1 910.73 kPa,犁底層緊實(shí)度約為耕作層的2倍。

    2)7個(gè)控制土層深度的初始入滲速率、穩(wěn)定入滲速率、平均入滲速和累計(jì)入滲量均表現(xiàn)為隨著入滲面深度遞增,數(shù)值遞減,其中0~15 cm坑深范圍內(nèi)平均的平均入滲速率和累計(jì)入滲量分別為>20~30 cm的17.04倍和18.06倍,耕作層與耕作層以下的土壤入滲能力差異較大。

    3)通過對(duì)比初始含水率和滲透48 h后含水率的顯著性得出,坑深小于15 cm入滲時(shí),水分入滲深度均為距表土20 cm,而坑深大于等于15 cm入滲時(shí),水分入滲深度均為坑底以下10 cm,這表明犁底層具有保水阻滲效果。

    4)Horton、Kostiakov和Philip 3個(gè)模型中,Kostiakov模型的擬合參數(shù)與實(shí)際相符,可較好地模擬水稻土入滲過程。土壤物理和入滲參數(shù)的相關(guān)性表明,入滲參數(shù)與容重、含水率、總孔隙度和田間持水量有著較強(qiáng)的相關(guān)性,與土壤緊實(shí)度相關(guān)性不顯著。

    [1] 李春林,丁啟朔,陳青春. 水稻土的先期固結(jié)壓力測(cè)定與分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(8):141-144.

    LI Chunlin, DING Qishuo, CHEN Qingchun. Measurement and analysis of precompression stress of soil in rice field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(8): 141-144. (in Chinese with English abstract)

    [2] KUMAR V, BELLINDER R R, GUPTA R K, et al. Role of herbicide-resistant rice in promoting resource conservation technologies in rice-wheat cropping systems of India: a review[J]. Crop Protection, 2008, 27(3/4/5): 290-301.

    [3] 劉俊安,王曉燕,李洪文,等. 基于土壤擾動(dòng)與牽引阻力的深松鏟結(jié)構(gòu)參數(shù)優(yōu)化[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(2):60-67.

    LIU Jun’an, WANG Xiaoyan, LI Hongwen, et al. Optimization of structural parameters of subsoiler based on soil disturbance and traction resistance[J]. Transaction of the Chinese Society for Agricultural Machinery, 2017, 48(2): 60-67. (in Chinese with English abstract)

    [4] 鄭侃,何進(jìn),李洪文,等. 基于離散元深松土壤模型的折線破土刃深松鏟研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(9):62-72.

    ZHENG Kan, HE Jin, LI Hongwen, et al. Research on polyline soil-breaking blade subsoiler based on subsoiling soil model using discrete element method[J]. Transaction of the Chinese Society for Agricultural Machinery, 2016, 47(9): 62-72. (in Chinese with English abstract)

    [5] WANG C, MAO X, HATANO R. Modeling ponded infiltration in fine textured soils with coarse interlayer[J]. Soil Science Society of America Journal, 2014, 78(3): 745-753.

    [6] ZETTL J, LEE B S, HUANG M, et al. Influence of textural layering on field capacity of coarse soils[J]. Canadian Journal of Soil Science, 2011, 91(2): 133-147.

    [7] 侯賢清,賈志寬,韓清芳,等. 不同輪耕模式對(duì)旱地土壤結(jié)構(gòu)及入滲蓄水特性的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(5):85-94.

    HOU Xianqing, JIA Zhikuan, HAN Qingfang, et al. Effects of different rotational tillage patterns on soil structure, infiltration and water storage characteristics in dryland[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(5): 85-94. (in Chinese with English abstract)

    [8] 張西平,程伍群,繩莉麗,等. 耕作及種植方式對(duì)土壤入滲參數(shù)和畦灌水流運(yùn)動(dòng)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(12):89-97.

    ZHANG Xiping, CHENG Wuqun, SHENG Lili, et al. Effects of tillage and planting patterns on soil infiltration parameters and water flow of border irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(12): 89-97. (in Chinese with English abstract)

    [9] 姚毓香,高喜杰,高鵬洋,等. 深松鏟橫向間距對(duì)土壤水分入滲的影響[J]. 農(nóng)機(jī)化研究,2020,42(9):18-23.

    YAO Yuxiang, GAO Xijie, GAO Pengyang, et al. Influence of tine spacing on soil-water infiltration[J]. Journal of Agricultural Mechanization Research, 2020, 42(9): 18-23. (in Chinese with English abstract)

    [10] 安晶,張玉龍,鄒洪濤,等. 深松和壓實(shí)對(duì)旱地棕壤硬度和水分入滲性能的影響[J]. 土壤通報(bào),2016,47(5):1107-1112.

    AN Jing, ZHANG Yulong, ZOU Hongtao, et al. Quantifying the influence of deep ripping and compaction on soil strength and water infiltration of brown earth[J]. Chinese Journal of Soil Science, 2016, 47(5): 1107-1112. (in Chinese with English abstract)

    [11] 孫蓓,馬玉瑩,雷廷武,等. 農(nóng)地耕層與犁底層土壤入滲性能的連續(xù)測(cè)量方法[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(4):118-124.

    SUN Bei, MA Yuying, LEI Tingwu, et al. Method for continuously measuring soil infiltrability of plow-cultivated layer and hardpan in farmland[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(4): 118-124. (in Chinese with English abstract)

    [12] 劉淙琮,王利書,楊寶斌. 不同農(nóng)田土壤壓實(shí)層深度對(duì)水分入滲和分布的影響[J]. 北方園藝,2021(19):107-113.

    LIU Congcong, WANG Lishu, YANG Baobin. Effects of soil compaction depth on water in filtration and distribution in different farmland[J]. Northern Horticulture, 2021(19): 107-113. (in Chinese with English abstract)

    [13] 依艷麗. 土壤物理研究法[M]. 北京:北京大學(xué)出版社,2009.

    [14] 孟凡旭,王樹森,馬迎梅,等. 不同果農(nóng)復(fù)合種植模式土壤入滲能力及其影響因素[J]. 干旱區(qū)研究,2020,37(6):1469-1477.

    MENG Fanxu, WANG Shusen, MA Yingmei, et al. Soil infiltration ability and influencing factors under different intercropping patterns in orchards in the gully region of loess plateau[J]. Arid Zone Research, 2020, 37(6): 1469-1477. (in Chinese with English abstract)

    [15] 何苗苗,劉芝芹,王克勤,等. 滇池流域不同植被覆蓋土壤的入滲特征及其影響因素[J]. 水土保持學(xué)報(bào),2022,36(3):181-187.

    HE Miaomiao, LIU Zhiqin, WANG Keqin, et al. Soil infiltration characteristics and influencing factors under different vegetation cover types in Dianchi Lake Basin[J]. Journal of Soil and Water Conservation, 2022, 36(3): 181-187. (in Chinese with English abstract)

    [16] 吳佩瑤,秦富倉,董曉宇,等. 坡度影響下的砒砂巖區(qū)裸露坡面水分入滲特征及模擬[J]. 水土保持學(xué)報(bào),2022,36(2):64-69.

    WU Peiyao, QIN Fucang, DONG Xiaoyu, et al. Water infiltration characteristics and simulation on exposed slope in pisha rock area under the influence of slope[J]. Journal of Soil and Water Conservation, 2022, 36(2): 64-69. (in Chinese with English abstract)

    [17] 江蘇省農(nóng)業(yè)農(nóng)村廳. 麥秸稈機(jī)械化全量還田與水稻機(jī)插秧集成技術(shù)[EB/OL]. (2014-06-04)[2022-06-28]. http://nynct.jiangsu.gov.cn/art/2014/6/4/art_13467_5056407.html

    [18] 丁啟朔,葛雙洋,任駿,等. 水稻土深松阻力與土壤擾動(dòng)效果研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2017,48(1):47-56,63.

    DING Qishuo, GE Shuangyang, REN Jun, et al. Characteristics of subsoiler traction and soil disturbance in paddy soil[J]. Transaction of the Chinese Society for Agricultural Machinery, 2017, 48(1): 47-56, 63. (in Chinese with English abstract)

    [19] HUO L F, LIANG L, ABBAS A, et al. Soil disturbance under small harvester traffic in paddy-based smallholder farms in China[J]. Agronomy Journal, 2020, 112(2): 1441-1451.

    [20] ZHAO W X, LI J S, YANG R M, et al. Crop yield and water productivity responses in management zones for variable-rate irrigation based on available soil water holding capacity[J]. Transactions of the ASABE, 2017, 60(5): 1659 - 1667.

    [21] 陳磊,張科利,李業(yè)桐,等. 喀斯特坡耕地塊石出露對(duì)土壤水分入滲的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(10):111-118.

    CHEN Lei, ZHANG Keli, LI Yetong, et al. Effects of rock outcropping on soil water infiltration in karst cropland slope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(10): 111-118. (in Chinese with English abstract)

    [22] 羅瑩麗,韋杰,劉春紅. 紫色土坡耕地埂坎裂隙發(fā)育對(duì)土壤入滲的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(21):116-123.

    LUO Yingli, WEI Jie, LIU Chunhong. Water seepage in soil bunds under different crack development degrees on the purple-soil sloping farmlands[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(21): 116-123. (in Chinese with English abstract)

    [23] 魏永霞,王鶴,劉慧,等. 生物炭對(duì)黑土區(qū)土壤水分及其入滲性能的影響[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(9):290-299.

    WEI Yongxia, WANG He, LIU Hui, et al. Effect of Biochar on Soil Moisture and Its Infiltration Performance in Black Soil Area[J]. Transaction of the Chinese Society for Agricultural Machinery, 2019, 50(9): 290-299. (in Chinese with English abstract)

    [24] 劉月梅,張興昌. EN-1對(duì)黃土性固化土水分垂直入滲特征的影響[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2012,43(11):65-73.

    LIU Yuemei, ZHANG Xingchang. Effects of EN-1 soil stabilizer on water vertical infiltration characteristics for loess stabilized soil[J]. Transaction of the Chinese Society for Agricultural Machinery, 2012, 43(11): 65-73. (in Chinese with English abstract)

    [25] 王慧杰,郝建平,馮瑞云,等. 微孔深松耕降低土壤緊實(shí)度提高棉花產(chǎn)量與種籽品質(zhì)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(8):7-14.

    WANG Huijie, HAO Jianping, FENG Ruiyun, et al. Microhole subsoiling decreasing soil compaction, and improving yield and seed quality of cotton[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(8): 7-14. (in Chinese with English abstract)

    [26] JE?áBEK J, ZUMR D, DOSTáL T. Identifying the plough pan position on cultivated soils by measurements of electrical resistivity and penetration resistance[J]. Soil and Tillage Research, 2017, 174: 231-240.

    [27] 王立梅. “犁底層”對(duì)北方半干旱地區(qū)農(nóng)田水循環(huán)的影響[J]. 鄉(xiāng)村科技,2019(15):119,121.

    [28] 吉恒瑩,邵明安,賈小旭. 土壤剖面結(jié)構(gòu)特征對(duì)坡面產(chǎn)流產(chǎn)沙過程的影響[J]. 土壤通報(bào),2018,49(2):441-446.

    JI Hengying, SHAO Min’an, JIA Xiaoxu. Impact of soil profile structure on infiltration and erosion processes[J]. Chinese Journal of Soil Science, 2018, 49(2): 441-446. (in Chinese with English abstract)

    [29] 云慧雅,畢華興,王珊珊,等. 不同林分類型土壤理化特征及其對(duì)土壤入滲過程的影響[J]. 水土保持學(xué)報(bào),2021,35(6):183-189.

    YUN Huiya, BI Huaxing, WANG Shanshan, et al. Soil physical and chemical characteristics of different forest types and their effects on soil infiltration process[J]. Journal of Soil and Water Conservation, 2021, 35(6): 183-189. (in Chinese with English abstract)

    [30] 方會(huì)敏,姬長(zhǎng)英,F(xiàn)arman Ali Chandio,等. 基于離散元法的旋耕過程土壤運(yùn)動(dòng)行為分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(3):22-28.

    FANG Huimin, JI Changying, FARMAN Ali Chandio, et al. Analysis of soil dynamic behavior during rotary tillage based on distinct element method[J]. Transaction of the Chinese Society for Agricultural Machinery, 2016, 47(3): 22-28. (in Chinese with English abstract)

    [31] 丁啟朔,呂緒敏,孫克潤(rùn),等. 基于逆向工程的攪漿刀作業(yè)性能界定與關(guān)鍵設(shè)計(jì)參數(shù)研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2021,52(3):68-74.

    DING Qishuo, LYU Xumin, SUN Kerun, et al. Performance and design parameter definition on puddling knife using reverse engineering[J]. Transaction of the Chinese Society for Agricultural Machinery, 2021, 52(3): 68-74. (in Chinese with English abstract)

    [32] 張志華,桑玉強(qiáng),孔玉華,等. 太行山低山丘陵區(qū)坡位對(duì)土壤滲透性能影響研究[J]. 灌溉排水學(xué)報(bào),2018,37(9):43-47.

    ZHANG Zhihua, SANG Yuqiang, KONG Yuhua, et al. The variation of water infiltration rate over various slopes in taihang mountain[J]. Journal of Irrigation and Drainage, 2018, 37(9): 43-47. (in Chinese with English abstract)

    [33] 劉姍姍,白美健,許迪,等. Green-Ampt模型參數(shù)簡(jiǎn)化及與土壤物理參數(shù)的關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(1):106-110.

    LIU Shanshan, BAI Meijian, XU Di, et al. Parameters simplification of Green-Ampt infiltration models and relationships between infiltration and soil physical parameters[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(1): 106-110. (in Chinese with English abstract)

    [34] 曾辰,王全九,樊軍. 初始含水率對(duì)土壤垂直線源入滲特征的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2010,26(1):24-30.

    ZENG Chen, WANG Quanjiu, FAN Jun. Effect of initial water content on vertical line-source infiltration characteristics of soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(1): 24-30. (in Chinese with English abstract)

    [35] 丁啟朔,李楊,BELAL Eisa Adam,等. 基于田間攝像的多參數(shù)水稻土深松擾動(dòng)行為與效應(yīng)研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(10):44-55.

    DING Qishuo, LI Yang, BELAL Eisa Adam, et al. Subsoiler-induced paddy soil disturbance and effects based on video-assisted multi-index quantification[J]. Transaction of the Chinese Society for Agricultural Machinery, 2019, 50(10): 44-55. (in Chinese with English abstract)

    [36] MURPHY B W, CRAWFORD M H, DUNCAN D A, et al. The use of visual soil assessment schemes to evaluate surface structure in a soil monitoring program[J]. Soil and Tillage Research, 2013, 127: 3-12.

    [37] 宋強(qiáng),夏可,楊斌,等. 旱作區(qū)典型土類穿透阻力分布特征及耕層厚度確定[J]. 水土保持學(xué)報(bào),2021,35(3):369-377.

    SONG Qiang, XIA Ke, YANG Bin, et al. Distribution characteristic of penetration resistance and determination of plough layer thickness of typical soils in dry farming area[J]. Journal of Soil and Water Conservation, 2021, 35(3): 369-377. (in Chinese with English abstract)

    [38] 喬磊,黃明鏡,張吳平,等. 山西省典型縣域農(nóng)田耕作層厚度空間格局及影響因素[J]. 中國(guó)土壤與肥料,2020(2):75-82.

    QIAO Lei, HUANG Mingjing, ZHANG Wuping, et al. Spatial pattern of cultivated layer thickness and its influencing factors in typical area of Shanxi province[J]. Soil and Fertilizer Sciences in China, 2020(2): 75-82. (in Chinese with English abstract)

    [39] GAO L, BECKER E, LIANG G, et al. Effect of different tillage systems on aggregate structure and inner distribution of organic carbon[J]. Geoderma, 2017, 288: 97-104.

    [40] 李勇,趙云澤,勾宇軒,等. 黃淮海旱作區(qū)土壤壓實(shí)度空間分布特征及其影響因素[J]. 農(nóng)業(yè)工程學(xué)報(bào),2021,37(13):83-91.

    LI Yong, ZHAO Yunze, GOU Yuxuan, et al. Spatial distribution characteristics and influence factors of degree of compaction in dry-farming Huang-Huai-Hai Plain of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(13): 83-91. (in Chinese with English abstract)

    [41] 余林,肖復(fù)明,徐海寧,等. 江西安福不同類型毛竹林土壤入滲與貯水特征[J]. 林業(yè)科技開發(fā),2013,27(2):37-41.

    YU Lin, XIAO Fuming, XU Haining, et al. Soil water holding capacities and infiltration characteristics in different types of Phyllostachys edulis plantations at Anfu, Jiangxi Province[J]. Journal of Forestry Engineering, 2013, 27(2): 37-41. (in Chinese with English abstract)

    [42] 包含,侯立柱,劉江濤,等. 室內(nèi)模擬降雨條件下土壤水分入滲及再分布試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2011,27(7):70-75.

    BAO Han, HOU Lizhu, LIU Jiangtao, et al. Experiment on process of soil water infiltration and redistribution under simulated rainfall[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 70-75. (in Chinese with English abstract)

    [43] 王曉彤,胡振琪,賴小君,等. 黏土夾層位置對(duì)黃河泥沙充填復(fù)墾土壤水分入滲的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019, 35(18):87-94.

    WANG Xiaotong, HU Zhenqi, LAI Xiaojun, et al. Influence of clay interlayer position on infiltration of reclaimed soil filled with Yellow River sediment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(18): 86-93. (in Chinese with English abstract)

    [44] 肖慶禮,黃明斌,邵明安,等. 黑河中游綠洲不同質(zhì)地土壤水分的入滲與再分布[J]. 農(nóng)業(yè)工程學(xué)報(bào),2014,30(2):124-131.

    XIAO Qingli, HUANG Mingbin, SHAO Ming’an, et al. Infiltration and drainage processes of different textural soil moisture in middle reaches of Heihe River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(2): 124-131. (in Chinese with English abstract)

    [45] YI J, QIU W W, HU W, et al. Effects of cultivation history in paddy rice on vertical water flows and related soil properties[J]. Soil and Tillage Research, 2020, 200: 104613.

    [46] 黨宏宇,陳洪松,邵明安. 喀斯特地區(qū)不同層次土石混合介質(zhì)對(duì)土壤水分入滲過程的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(8):38-43.

    DANG Hongyu, CHEN Hongsong, SHAO Ming'an. Effects of laminated rock fragments on soil infiltration processes in Karst regions[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(8): 38-43. (in Chinese with English abstract)

    Soil infiltration of rice-wheat rotation field under long-term mechanical treatment based on field in situ experiments

    LIANG Lei, SUN Haotian, XU Gaoming, LI Yinian, HE Ruiyin, DING Qishuo※

    (1.210031,;2.210031,)

    Field in situ infiltration experiments at different depths can be used to express the state of soil stratification, and the physical differentiation of soil layers, in order to quantify the changes in the water of the soil profile. This study aims to obtain the infiltration capacity and water retention of paddy soils. The field in situ soil infiltration experiments were also carried out at different depths to identify the soil stratification in the representative plot under long-term mechanized tillage of smallholder farmers in the rice-wheat rotation region of eastern China. Seven infiltration pits of different depths were excavated in the experimental plots. After that, the infiltration experiments were conducted at the bottom of the pits. Soil water content was then measured in the layers for 48 hours of infiltration. A systematic investigation was made to explore the soil infiltration capacity at the bottom of the pits and the soil water content of each layer. The results showed that an accurate description was achieved in the water infiltration and water holding capacity of the soil at different pits after the infiltration experiments at different depths. The location and thickness of the plow pan were also clearly identified during this time. The plow pan was started at a depth of about 15 cm, indicating the outstanding difference between the cultivated layer and the plow pan. The average soil penetration resistances of the cultivated layer and the plow pan were 1 005.79 kPa, and 1 910.73 kPa, respectively. The soil profile showed that the soil in the cultivated layer shared a loose morphology and dense root distribution, whereas, the soil in the plow pan presented a high bulk density, low porosity, and poor permeability, while the soil in the subsoil layer was in more iron-manganese spots and poor structure. Furthermore, the soil infiltration parameters decreased with the increasing pit depth. The average infiltration rate and cumulative infiltration in the 0-15 cm pit depth range were 17.04 and 18.06 times higher, respectively, than those in the >20-30 cm pit depth range. The three infiltration models were fitted using Horton, Kostiakov and Philip. Specifically, the Kostiakov model had the highest2(0.98-0.99) and the smallest root mean square error (0.01-0.77 mm/min), indicating the consistency of the fitted parameters. Infiltration parameters were extremely significant correlated with the soil bulk density, water content, total porosity, and field capacity (<0.01), but not with the soil penetration resistance (>0.05). Field in situ infiltration experiment at different depths was an important tool to identify the soil stratification and quantify the water function differentiation in the soil profiles, according to the soil profile, cone penetration and sampling. Long-term mechanized tillage under the smallholder production model in the rice-wheat rotation region of eastern China can be expected to result in the apparent stratification and vertical differentiation of soil water functions in paddy soils. In turn, there are also some significant differences in the infiltration capacity between the cultivated layer and soil layers below the cultivated layer. This study can provide a reference for mechanical cultivation and irrigation in rice-wheat rotation regions.

    soils; infiltration; water; rice-wheat rotation; cultivated layer; plow pan

    2022-09-26

    2023-03-01

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2022YFD2300304);江蘇省研究生科研創(chuàng)新計(jì)劃項(xiàng)目(KYCX21_0573)

    梁磊,博士生,研究方向?yàn)橥寥老到y(tǒng)科學(xué)。Email:liangupright@163.com

    丁啟朔,博士,教授,研究方向?yàn)橥寥?植物-機(jī)器系統(tǒng)科學(xué)。Email:qsding@njau.edu.cn

    10.11975/j.issn.1002-6819.202209207

    S152.7+2; S152.2+3; S220

    A

    1002-6819(2023)-07-0110-09

    梁磊,孫浩田,徐高明,等. 田間原位試驗(yàn)分析長(zhǎng)期機(jī)械作業(yè)下稻麥輪作地塊土壤入滲性能[J]. 農(nóng)業(yè)工程學(xué)報(bào),2023,39(7):110-118. doi:10.11975/j.issn.1002-6819.202209207 http://www.tcsae.org

    LIANG Lei, SUN Haotian, XU Gaoming, et al. Soil infiltration of rice-wheat rotation field under long-term mechanical treatment based on field in situ experiments[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2023, 39(7): 110-118. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.202209207 http://www.tcsae.org

    猜你喜歡
    實(shí)度耕作層底層
    耕地耕作層資源保護(hù)利用的意義、問題和對(duì)策
    航天企業(yè)提升采購能力的底層邏輯
    土壤緊實(shí)度對(duì)毛竹冬筍和春筍產(chǎn)量的影響
    土壤緊實(shí)度對(duì)裸土侵蝕強(qiáng)度影響的實(shí)驗(yàn)與分析
    耕作層剝離再利用實(shí)踐與政策差距淺談
    耕作層剝離再利用實(shí)踐與政策差距分析
    基于加速度補(bǔ)償?shù)耐寥谰o實(shí)度測(cè)量方法與傳感器設(shè)計(jì)
    內(nèi)陸沙區(qū)丘間地土壤緊實(shí)度特征及其與植被的關(guān)系
    回到現(xiàn)實(shí)底層與悲憫情懷
    小說林(2014年5期)2014-02-28 19:51:47
    略論“底層”
    雜文選刊(2013年7期)2013-02-11 10:41:11
    舔av片在线| 99riav亚洲国产免费| 免费在线观看亚洲国产| 美女被艹到高潮喷水动态| 日韩精品中文字幕看吧| 欧美又色又爽又黄视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧美一区二区三区黑人| 午夜福利成人在线免费观看| 在线十欧美十亚洲十日本专区| 日韩欧美 国产精品| 国产私拍福利视频在线观看| 搡老妇女老女人老熟妇| 最新在线观看一区二区三区| 欧美高清成人免费视频www| 日韩欧美 国产精品| 99精品欧美一区二区三区四区| 91久久精品国产一区二区成人 | 亚洲色图 男人天堂 中文字幕| svipshipincom国产片| 在线十欧美十亚洲十日本专区| 女人高潮潮喷娇喘18禁视频| 国产精品女同一区二区软件 | 一a级毛片在线观看| 制服丝袜大香蕉在线| 免费在线观看亚洲国产| 在线国产一区二区在线| 伊人久久大香线蕉亚洲五| 国产成人福利小说| 午夜激情欧美在线| 国产伦一二天堂av在线观看| 亚洲真实伦在线观看| 中文字幕人成人乱码亚洲影| 亚洲片人在线观看| 午夜日韩欧美国产| 熟女电影av网| 欧美日韩瑟瑟在线播放| 欧美激情久久久久久爽电影| 国产爱豆传媒在线观看| 97人妻精品一区二区三区麻豆| 男女床上黄色一级片免费看| 叶爱在线成人免费视频播放| 国产亚洲av嫩草精品影院| a在线观看视频网站| 三级毛片av免费| 一卡2卡三卡四卡精品乱码亚洲| 国产成+人综合+亚洲专区| 婷婷精品国产亚洲av在线| 亚洲av日韩精品久久久久久密| 成人av在线播放网站| 一级毛片女人18水好多| 国内精品美女久久久久久| 黄片小视频在线播放| 高清毛片免费观看视频网站| 国产探花在线观看一区二区| 色哟哟哟哟哟哟| 欧美黄色片欧美黄色片| 男人的好看免费观看在线视频| 亚洲成a人片在线一区二区| 国产 一区 欧美 日韩| 日日摸夜夜添夜夜添小说| 一二三四在线观看免费中文在| 免费在线观看影片大全网站| 日韩成人在线观看一区二区三区| 久久久成人免费电影| 午夜福利欧美成人| 色噜噜av男人的天堂激情| www.999成人在线观看| 免费观看的影片在线观看| 最近最新中文字幕大全电影3| 日本成人三级电影网站| 91麻豆精品激情在线观看国产| 成人无遮挡网站| 亚洲第一电影网av| 久久精品91蜜桃| 欧美黄色片欧美黄色片| 在线免费观看不下载黄p国产 | 一二三四社区在线视频社区8| 欧美成狂野欧美在线观看| 日本一本二区三区精品| 久久久久久九九精品二区国产| 哪里可以看免费的av片| 狂野欧美白嫩少妇大欣赏| 亚洲国产欧美网| 村上凉子中文字幕在线| 国产黄色小视频在线观看| 国产97色在线日韩免费| 成年版毛片免费区| 亚洲激情在线av| 青草久久国产| 亚洲国产欧美网| 久久久国产欧美日韩av| 亚洲av熟女| 精品久久久久久成人av| 国产一级毛片七仙女欲春2| 一二三四社区在线视频社区8| 在线观看免费午夜福利视频| 免费看光身美女| www日本黄色视频网| 亚洲av熟女| 在线播放国产精品三级| 国内少妇人妻偷人精品xxx网站 | 久久久久久久精品吃奶| 国内精品久久久久精免费| ponron亚洲| 国产av不卡久久| 国产精品 国内视频| 女警被强在线播放| 真人做人爱边吃奶动态| 国产精品久久久久久久电影 | 又紧又爽又黄一区二区| av欧美777| 麻豆av在线久日| 男女那种视频在线观看| 99国产极品粉嫩在线观看| 在线观看美女被高潮喷水网站 | 母亲3免费完整高清在线观看| 国产久久久一区二区三区| 在线观看免费视频日本深夜| 哪里可以看免费的av片| 特大巨黑吊av在线直播| 欧美日韩黄片免| 国产一区二区在线av高清观看| 精品无人区乱码1区二区| 国产伦在线观看视频一区| 美女午夜性视频免费| 又大又爽又粗| 国产精品电影一区二区三区| 巨乳人妻的诱惑在线观看| 久久亚洲真实| 国产av麻豆久久久久久久| 嫩草影视91久久| 性色av乱码一区二区三区2| 欧美日韩亚洲国产一区二区在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲精品色激情综合| 国产美女午夜福利| 中文字幕久久专区| 白带黄色成豆腐渣| 国产精品av视频在线免费观看| 俺也久久电影网| 窝窝影院91人妻| 女人高潮潮喷娇喘18禁视频| 男人的好看免费观看在线视频| 日日摸夜夜添夜夜添小说| 日韩成人在线观看一区二区三区| 亚洲精品久久国产高清桃花| 国产成人精品久久二区二区免费| 久久精品国产亚洲av香蕉五月| 亚洲欧美日韩高清专用| 搞女人的毛片| 一二三四社区在线视频社区8| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品合色在线| 久久久久久国产a免费观看| 夜夜夜夜夜久久久久| 欧美日韩精品网址| 露出奶头的视频| 亚洲欧美精品综合久久99| 可以在线观看毛片的网站| 免费看光身美女| 亚洲国产精品999在线| 免费人成视频x8x8入口观看| 久久久精品大字幕| 最好的美女福利视频网| 久久九九热精品免费| 午夜日韩欧美国产| 欧美黄色片欧美黄色片| 精品无人区乱码1区二区| 国产毛片a区久久久久| 日韩欧美三级三区| 91麻豆av在线| 日本免费一区二区三区高清不卡| 欧美成人性av电影在线观看| 怎么达到女性高潮| 亚洲欧美日韩无卡精品| 麻豆国产av国片精品| 免费在线观看视频国产中文字幕亚洲| 免费高清视频大片| 亚洲男人的天堂狠狠| 国产一区二区三区视频了| 99国产综合亚洲精品| 日韩欧美在线二视频| 久久香蕉精品热| 久久久久久国产a免费观看| 欧美成人免费av一区二区三区| 中文亚洲av片在线观看爽| 九九热线精品视视频播放| 亚洲avbb在线观看| 成人av一区二区三区在线看| 亚洲黑人精品在线| 国产亚洲精品综合一区在线观看| 久久精品国产清高在天天线| 桃色一区二区三区在线观看| 12—13女人毛片做爰片一| 51午夜福利影视在线观看| 男女那种视频在线观看| 日本 av在线| 亚洲国产色片| 国产精品1区2区在线观看.| 成年女人永久免费观看视频| 特大巨黑吊av在线直播| 欧美又色又爽又黄视频| tocl精华| 日本成人三级电影网站| а√天堂www在线а√下载| 男女那种视频在线观看| 精品国产三级普通话版| 亚洲一区高清亚洲精品| 欧美三级亚洲精品| 69av精品久久久久久| 亚洲在线自拍视频| 最新中文字幕久久久久 | 又紧又爽又黄一区二区| 亚洲精品中文字幕一二三四区| 夜夜躁狠狠躁天天躁| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 国产毛片a区久久久久| 国产成人福利小说| 两个人看的免费小视频| 亚洲国产看品久久| av片东京热男人的天堂| 国产爱豆传媒在线观看| 一本久久中文字幕| 网址你懂的国产日韩在线| 琪琪午夜伦伦电影理论片6080| 亚洲精品乱码久久久v下载方式 | 51午夜福利影视在线观看| 看免费av毛片| 久久中文看片网| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 嫩草影院精品99| 99久久精品国产亚洲精品| 色综合婷婷激情| 男女做爰动态图高潮gif福利片| 亚洲人成伊人成综合网2020| 国产高清视频在线播放一区| 嫁个100分男人电影在线观看| 国产熟女xx| 欧美黑人欧美精品刺激| 校园春色视频在线观看| 精品久久久久久,| 久久午夜亚洲精品久久| 99视频精品全部免费 在线 | 午夜视频精品福利| 久久久久久久精品吃奶| av黄色大香蕉| 亚洲成人免费电影在线观看| 男女床上黄色一级片免费看| 欧美zozozo另类| 级片在线观看| 特级一级黄色大片| 日本精品一区二区三区蜜桃| 成熟少妇高潮喷水视频| 好男人电影高清在线观看| 99riav亚洲国产免费| 99国产精品一区二区三区| 哪里可以看免费的av片| 黑人欧美特级aaaaaa片| 黄色女人牲交| 国产成人欧美在线观看| 日本五十路高清| 欧美成狂野欧美在线观看| 老熟妇乱子伦视频在线观看| 精品国产乱码久久久久久男人| 免费观看精品视频网站| 97超视频在线观看视频| 99久久成人亚洲精品观看| 精品久久久久久,| 久99久视频精品免费| 在线a可以看的网站| 偷拍熟女少妇极品色| 亚洲 欧美一区二区三区| 国产精品久久久久久精品电影| 美女高潮的动态| 日本黄大片高清| 亚洲国产欧美一区二区综合| 成人高潮视频无遮挡免费网站| 九九热线精品视视频播放| 夜夜躁狠狠躁天天躁| 黄色女人牲交| 亚洲精品乱码久久久v下载方式 | 亚洲国产欧美一区二区综合| 成人国产综合亚洲| 一夜夜www| 中文字幕高清在线视频| 亚洲精品美女久久久久99蜜臀| 欧美国产日韩亚洲一区| 两个人的视频大全免费| 久久性视频一级片| 99视频精品全部免费 在线 | 男人舔女人的私密视频| 一级作爱视频免费观看| 国产一级毛片七仙女欲春2| 国产精品 国内视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美性猛交╳xxx乱大交人| 亚洲人成网站高清观看| 黄片小视频在线播放| 婷婷丁香在线五月| 欧美三级亚洲精品| 琪琪午夜伦伦电影理论片6080| 性欧美人与动物交配| 黑人巨大精品欧美一区二区mp4| 中文字幕人妻丝袜一区二区| 高清毛片免费观看视频网站| 亚洲欧美日韩高清在线视频| 亚洲中文字幕一区二区三区有码在线看 | 在线永久观看黄色视频| 国产97色在线日韩免费| 色在线成人网| 国产单亲对白刺激| 少妇人妻一区二区三区视频| 欧美中文日本在线观看视频| 亚洲欧美日韩无卡精品| 狂野欧美激情性xxxx| 午夜福利在线在线| 国内精品久久久久久久电影| 成人国产一区最新在线观看| 国产成人av教育| 悠悠久久av| 午夜免费激情av| 老司机在亚洲福利影院| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 亚洲国产欧美人成| 国产单亲对白刺激| 曰老女人黄片| 性色av乱码一区二区三区2| 亚洲 欧美一区二区三区| www日本黄色视频网| 老鸭窝网址在线观看| 精品久久久久久久久久久久久| 全区人妻精品视频| 麻豆国产97在线/欧美| e午夜精品久久久久久久| 动漫黄色视频在线观看| 性色avwww在线观看| 亚洲专区国产一区二区| 丝袜人妻中文字幕| 成人特级黄色片久久久久久久| 久久精品91无色码中文字幕| 亚洲人成网站在线播放欧美日韩| 亚洲五月婷婷丁香| 亚洲一区二区三区色噜噜| 午夜福利视频1000在线观看| 国内精品美女久久久久久| 国产成人精品久久二区二区91| 国产伦精品一区二区三区四那| 精品不卡国产一区二区三区| www.www免费av| 亚洲欧美精品综合一区二区三区| 中亚洲国语对白在线视频| 网址你懂的国产日韩在线| 国产成人精品无人区| 欧美绝顶高潮抽搐喷水| 亚洲aⅴ乱码一区二区在线播放| www日本在线高清视频| 国产伦人伦偷精品视频| 日韩精品中文字幕看吧| 免费在线观看成人毛片| 黄色 视频免费看| 十八禁人妻一区二区| 麻豆成人av在线观看| 嫩草影院入口| 女同久久另类99精品国产91| 俺也久久电影网| 日本黄色视频三级网站网址| 国产淫片久久久久久久久 | 国产久久久一区二区三区| 国产高清激情床上av| 精品99又大又爽又粗少妇毛片 | 99精品在免费线老司机午夜| 欧美一区二区国产精品久久精品| 99久久99久久久精品蜜桃| 嫩草影院精品99| 香蕉久久夜色| 欧美激情在线99| 曰老女人黄片| 99riav亚洲国产免费| 亚洲成人免费电影在线观看| 麻豆国产av国片精品| 国产高清视频在线播放一区| x7x7x7水蜜桃| 亚洲成人中文字幕在线播放| 成人亚洲精品av一区二区| 色噜噜av男人的天堂激情| 亚洲国产中文字幕在线视频| 欧美午夜高清在线| 啦啦啦观看免费观看视频高清| 久久精品夜夜夜夜夜久久蜜豆| 18禁黄网站禁片午夜丰满| 欧美日韩国产亚洲二区| 露出奶头的视频| 国产熟女xx| 一个人免费在线观看电影 | 黄片大片在线免费观看| 精品国产亚洲在线| 国产精品美女特级片免费视频播放器 | 中亚洲国语对白在线视频| 日本成人三级电影网站| 1024手机看黄色片| 国产精品影院久久| 特大巨黑吊av在线直播| 久久天堂一区二区三区四区| 中文字幕久久专区| 成在线人永久免费视频| 色老头精品视频在线观看| 免费人成视频x8x8入口观看| 亚洲中文字幕日韩| 国产淫片久久久久久久久 | 日本一二三区视频观看| 亚洲国产精品合色在线| 欧美xxxx黑人xx丫x性爽| 国产亚洲av高清不卡| 亚洲黑人精品在线| 一个人看视频在线观看www免费 | 亚洲精品美女久久av网站| 精品久久久久久久久久久久久| 亚洲成人精品中文字幕电影| 国产亚洲精品久久久com| 亚洲精品在线美女| 无人区码免费观看不卡| 亚洲乱码一区二区免费版| 国产成人系列免费观看| 91麻豆av在线| 免费电影在线观看免费观看| 中文字幕久久专区| av在线天堂中文字幕| 成人高潮视频无遮挡免费网站| 日本 欧美在线| 中文亚洲av片在线观看爽| 欧美色视频一区免费| 亚洲 国产 在线| 国产精品亚洲av一区麻豆| 国产又色又爽无遮挡免费看| 亚洲精品乱码久久久v下载方式 | 久99久视频精品免费| 国产av在哪里看| 久久草成人影院| 在线观看一区二区三区| 两个人看的免费小视频| 午夜福利免费观看在线| 国产精品国产高清国产av| 成人无遮挡网站| 人妻夜夜爽99麻豆av| 亚洲人成伊人成综合网2020| 在线观看一区二区三区| 淫妇啪啪啪对白视频| 午夜福利免费观看在线| 国产毛片a区久久久久| 欧美黄色淫秽网站| 亚洲欧美日韩无卡精品| 级片在线观看| 久久久国产精品麻豆| 免费看光身美女| 亚洲精品久久国产高清桃花| 亚洲专区中文字幕在线| 国产视频内射| 人人妻人人澡欧美一区二区| 99久久国产精品久久久| 国内少妇人妻偷人精品xxx网站 | 国产三级在线视频| 蜜桃久久精品国产亚洲av| 男女床上黄色一级片免费看| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av| 精品国产亚洲在线| 久久久久性生活片| 九九久久精品国产亚洲av麻豆 | 亚洲美女视频黄频| av天堂在线播放| 桃色一区二区三区在线观看| 国产蜜桃级精品一区二区三区| 一级毛片女人18水好多| 观看美女的网站| 一级黄色大片毛片| 在线观看免费午夜福利视频| 美女大奶头视频| 久久久久精品国产欧美久久久| 久久久水蜜桃国产精品网| 久久久久免费精品人妻一区二区| 久久天躁狠狠躁夜夜2o2o| 我的老师免费观看完整版| 非洲黑人性xxxx精品又粗又长| 国产又色又爽无遮挡免费看| 中文字幕久久专区| 亚洲欧美一区二区三区黑人| 国产极品精品免费视频能看的| 老司机午夜十八禁免费视频| 久久99热这里只有精品18| 免费观看精品视频网站| 精品免费久久久久久久清纯| 国产极品精品免费视频能看的| xxxwww97欧美| 午夜免费观看网址| 亚洲成av人片在线播放无| 天天躁日日操中文字幕| 精品电影一区二区在线| 成人亚洲精品av一区二区| 波多野结衣高清作品| 国产精品99久久久久久久久| 十八禁人妻一区二区| 欧美又色又爽又黄视频| 黄色丝袜av网址大全| 欧美日韩中文字幕国产精品一区二区三区| 婷婷亚洲欧美| 操出白浆在线播放| 国内揄拍国产精品人妻在线| 啦啦啦韩国在线观看视频| 男女之事视频高清在线观看| 国产99白浆流出| 成人鲁丝片一二三区免费| 色吧在线观看| 国产伦一二天堂av在线观看| 一本久久中文字幕| 午夜福利18| 午夜成年电影在线免费观看| 国产成人福利小说| 成年版毛片免费区| 黄色片一级片一级黄色片| 中文字幕人成人乱码亚洲影| 国产精品电影一区二区三区| 亚洲av日韩精品久久久久久密| 91九色精品人成在线观看| 国产伦人伦偷精品视频| 老熟妇乱子伦视频在线观看| 国产精品久久视频播放| 小说图片视频综合网站| 亚洲精品久久国产高清桃花| 欧美日韩综合久久久久久 | 一本一本综合久久| 人妻夜夜爽99麻豆av| 在线观看免费午夜福利视频| 国产av在哪里看| www国产在线视频色| 久久精品夜夜夜夜夜久久蜜豆| 国产av一区在线观看免费| 亚洲 欧美一区二区三区| 哪里可以看免费的av片| 精品久久久久久久人妻蜜臀av| 九色成人免费人妻av| 黑人操中国人逼视频| 国内毛片毛片毛片毛片毛片| 超碰成人久久| 久久久久久久精品吃奶| 变态另类成人亚洲欧美熟女| a在线观看视频网站| 国产99白浆流出| 在线观看免费午夜福利视频| www.自偷自拍.com| 国产69精品久久久久777片 | 香蕉国产在线看| ponron亚洲| 少妇的丰满在线观看| 免费高清视频大片| 精品福利观看| 亚洲欧美激情综合另类| 亚洲av免费在线观看| 不卡av一区二区三区| 一级作爱视频免费观看| 国产精品av视频在线免费观看| 亚洲国产欧美网| 美女cb高潮喷水在线观看 | 淫妇啪啪啪对白视频| 岛国在线免费视频观看| ponron亚洲| 99久久成人亚洲精品观看| 久久中文字幕一级| 人妻夜夜爽99麻豆av| 亚洲男人的天堂狠狠| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品久久久久久毛片| 国产成人欧美在线观看| 19禁男女啪啪无遮挡网站| 国语自产精品视频在线第100页| aaaaa片日本免费| 国产成年人精品一区二区| 欧美激情久久久久久爽电影| 中文字幕精品亚洲无线码一区| 久久中文看片网| 在线十欧美十亚洲十日本专区| 国产美女午夜福利| 国产精品免费一区二区三区在线| 日韩国内少妇激情av| 一进一出抽搐动态| 精品国产美女av久久久久小说| 国产精品99久久久久久久久| 一级作爱视频免费观看| 他把我摸到了高潮在线观看| 老司机深夜福利视频在线观看| 亚洲人成电影免费在线| 香蕉av资源在线| 天天一区二区日本电影三级| 欧美日韩福利视频一区二区| 午夜免费观看网址| 欧美精品啪啪一区二区三区| 在线免费观看不下载黄p国产 | 欧美日韩国产亚洲二区| av女优亚洲男人天堂 | 亚洲av中文字字幕乱码综合| 中国美女看黄片| 国产精品久久久av美女十八| 国产激情偷乱视频一区二区| 国模一区二区三区四区视频 | 国产又黄又爽又无遮挡在线| 岛国在线免费视频观看| 久久久久亚洲av毛片大全| 精品久久久久久久久久免费视频| 啪啪无遮挡十八禁网站| 村上凉子中文字幕在线| 99久久精品一区二区三区| 亚洲五月婷婷丁香|