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

    水分虧缺程度對(duì)變量灌溉水分傳感器埋設(shè)位置預(yù)判的影響

    2018-11-24 01:37:08李秀梅趙偉霞李久生栗巖峰
    關(guān)鍵詞:黏粒灌水土層

    李秀梅,趙偉霞,李久生,栗巖峰

    ?

    水分虧缺程度對(duì)變量灌溉水分傳感器埋設(shè)位置預(yù)判的影響

    李秀梅,趙偉霞,李久生※,栗巖峰

    (中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100048)

    基于土壤黏粒含量遴選可代表平均土壤含水率的點(diǎn)位,是確定變量灌溉系統(tǒng)土壤水分傳感器網(wǎng)絡(luò)埋設(shè)位置的重要方法。為了評(píng)估水分虧缺程度對(duì)傳感器埋設(shè)位置的影響,該研究通過(guò)在不同土壤可利用水量管理區(qū)內(nèi)設(shè)置不同的水分虧缺灌溉處理,基于土壤含水率時(shí)間穩(wěn)定性原理,分析了水分虧缺程度和土壤性質(zhì)對(duì)冬小麥主要根系層內(nèi)土壤含水率空間分布結(jié)構(gòu)相似性和直接代表平均土壤含水率點(diǎn)位的影響。結(jié)果表明,受土壤性質(zhì)和水分虧缺程度的影響,在土壤砂粒含量隨土層深度增加而增大的1區(qū),土壤含水率空間分布結(jié)構(gòu)相似性?xún)H在2016年的雨養(yǎng)處理達(dá)到了顯著水平(0.05),而在土壤剖面質(zhì)地均勻的2區(qū),2016年的雨養(yǎng)和高、中、低水分虧缺處理和2017年的中、低水分虧缺處理均達(dá)到了顯著水平(0.05)。在冬小麥主要根系層內(nèi),不同土壤性質(zhì)管理區(qū)直接代表平均土壤含水率點(diǎn)位占總測(cè)點(diǎn)數(shù)的比例基本相等。受水分虧缺程度的影響,1區(qū)直接代表平均土壤含水率的測(cè)點(diǎn)比例隨水分虧缺程度減小而增加,2區(qū)則呈先減小后略有增加趨勢(shì)。除2016年高水分虧缺處理外,0~0.2、>0.2~0.4、>0.4~0.6 m土層代表平均土壤含水率點(diǎn)位的黏粒含量與該土層平均黏粒含量之間均存在顯著的線(xiàn)性關(guān)系(0.05),2 a擬合系數(shù)變化范圍為0.66~1.03,且2017年擬合系數(shù)隨土壤水分虧缺程度增加呈增大趨勢(shì)。因此,在砂壤土地塊變量灌溉管理區(qū)內(nèi),基于管理區(qū)平均土壤黏粒含量進(jìn)行土壤水分傳感器埋設(shè)位置遴選時(shí),需根據(jù)擬采用的水分虧缺管理模式對(duì)擬合系數(shù)進(jìn)行修正。

    土壤;水分;傳感器;網(wǎng)絡(luò);時(shí)間穩(wěn)定性;埋設(shè)位置;冬小麥;變量灌溉

    0 引 言

    變量灌溉處方圖的生成依賴(lài)于各種傳感器網(wǎng)絡(luò)實(shí)時(shí)監(jiān)測(cè)和反饋不同管理區(qū)內(nèi)土壤水分、作物生長(zhǎng)和氣象條件等信息[1-2]。土壤水分傳感器用于監(jiān)測(cè)土壤水分對(duì)作物的滿(mǎn)足程度,是生成變量灌溉處方圖常用的工具之一[3-5]。因?yàn)槊糠N類(lèi)型的土壤水分傳感器均有一定的測(cè)量范圍,因此確保土壤水分傳感器埋設(shè)位置處所測(cè)定的土壤含水率對(duì)田塊內(nèi)平均值的代表性是變量灌溉技術(shù)應(yīng)用時(shí)面臨的主要任務(wù)[4]。

    自1985年Vachaud等[6]介紹了土壤含水率時(shí)間穩(wěn)定性概念后,國(guó)內(nèi)外學(xué)者先后應(yīng)用土壤含水率時(shí)間穩(wěn)定性原理在不同土壤物理性質(zhì)[7-10]、海拔[11-12]、研究尺 度[13-14]和植被覆蓋[15]條件下開(kāi)展了直接代表平均土壤含水率(以下簡(jiǎn)稱(chēng)為“代表含水率”)的樣點(diǎn)位置研究。2017年,Zhao等[16]在變量灌溉管理中開(kāi)展了基于時(shí)間穩(wěn)定性原理的土壤水分傳感器埋設(shè)位置研究,通過(guò)在不同土壤可利用水量管理區(qū)設(shè)置相同的灌水下限,以均一灌溉管理為對(duì)照,研究了變量灌溉管理下主要根系層的土壤含水率空間結(jié)構(gòu)相似性和時(shí)間穩(wěn)定性,證明了當(dāng)作物耗水主要來(lái)源于灌溉時(shí),變量灌溉管理雖會(huì)改變農(nóng)田內(nèi)土壤水分空間格局的整體相似性,但每個(gè)管理區(qū)內(nèi)的土壤含水率時(shí)間穩(wěn)定性依然存在。研究還發(fā)現(xiàn),每個(gè)管理區(qū)內(nèi)代表含水率點(diǎn)位的黏粒含量與對(duì)應(yīng)管理區(qū)的平均黏粒含量之間存在的線(xiàn)性關(guān)系,可用于土壤性質(zhì)變異較大的農(nóng)田中土壤水分傳感器埋設(shè)位置的預(yù)判。

    因?yàn)闈駶?rùn)土壤中具有時(shí)間穩(wěn)定性的測(cè)點(diǎn)在土壤干燥過(guò)程中,其穩(wěn)定性可能消失[17-18],因此在不同的土壤水分狀態(tài)下,代表平均土壤含水率的點(diǎn)位并不完全相同[19]。Zhang等[20]也發(fā)現(xiàn)不同土壤水分條件下的土壤含水率大小在空間上的排序并不相同。因此,尋找不同水分虧缺程度下平均土壤含水率點(diǎn)位的“平均屬性”對(duì)確定平均土壤含水率點(diǎn)位更具指導(dǎo)意義[21]。

    為了充分發(fā)揮變量灌溉技術(shù)在干旱半干旱地區(qū)的節(jié)水潛力,將非充分灌溉原理與變量灌溉技術(shù)相結(jié)合是發(fā)展干旱半干旱條件下作物變量灌溉管理的關(guān)鍵[3-4, 22],因此本文基于不同水分虧缺灌溉管理?xiàng)l件下的土壤含水率時(shí)間穩(wěn)定性分析,研究不同濕潤(rùn)程度下土壤含水率空間分布格局的相似性,分析土壤性質(zhì)和土壤水分虧缺程度對(duì)遴選代表含水率點(diǎn)位的影響,并驗(yàn)證非充分灌溉制度下,以管理區(qū)內(nèi)平均土壤黏粒含量作為代表含水率點(diǎn)位的預(yù)判方法的可行性,研究結(jié)果將為非充分變量灌溉管理?xiàng)l件下土壤水分傳感器網(wǎng)絡(luò)布設(shè)方法的制定提供依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)地概況

    試驗(yàn)在河北省涿州市東城坊鎮(zhèn)進(jìn)行(39°27¢N,115°51¢E,海拔42 m)。該試驗(yàn)場(chǎng)位于河北省中部,坐落于太行山山洪沖積扇平原上,屬暖溫帶半濕潤(rùn)季風(fēng)區(qū)。多年平均降水量563.3 mm,年平均溫度11.6 ℃。土壤質(zhì)地以壤土和砂壤土為主。

    試驗(yàn)區(qū)為圓形噴灌機(jī)變量灌溉系統(tǒng)[23]控制區(qū)域的1/4,面積1.64 hm2。基于土壤可利用水量(available soil water holding capacity, AWC)將試驗(yàn)區(qū)劃分為4個(gè)管理區(qū),1~4區(qū)0~0.6 m田間持水量分別為0.21、0.22、0.23和0.25 cm3/cm3;0~1.0 m深度AWC變化范圍依次為152~161、>161~171、>171~185、>185~205 mm[24]。因?yàn)?區(qū)0.4 m以下土層主要為南水北調(diào)管道工程回填的砂礫石,本文分析未包含4區(qū)。在1~3區(qū),土壤性質(zhì)沿垂直剖面表現(xiàn)出不同的分層特征,1區(qū)0~0.2、>0.2~0.4、>0.4~0.6 m土層土壤砂粒體積分?jǐn)?shù)隨深度增加呈明顯增加趨勢(shì)(57.8%~89.5%);2區(qū)土壤砂粒和黏粒體積分?jǐn)?shù)隨深度變化較小,分別為59.8%~65.8%和4.6%~5.7%;3區(qū)>0.2~0.4 m土層存在明顯的黏土夾層,該土層砂粒和黏粒體積分?jǐn)?shù)分別為47.9%和8.6%[25]。

    1.2 試驗(yàn)設(shè)計(jì)

    3個(gè)管理區(qū)內(nèi),分別按照面積相等的原則設(shè)置灌水處理。1區(qū)和3區(qū)分別劃分為4個(gè)子區(qū),2區(qū)劃分為5個(gè)子區(qū),每個(gè)子區(qū)為1個(gè)灌水處理。2015—2016年冬小麥生育期內(nèi)(以下簡(jiǎn)稱(chēng)2016年冬小麥),設(shè)置為雨養(yǎng)和不同灌水下限(為田間持水量百分比)的處理(表1),灌水定額為20 mm。因?yàn)槊總€(gè)子區(qū)內(nèi)到達(dá)灌水下限的時(shí)間不同,所以采用部分噴頭開(kāi)啟式變量灌溉方法進(jìn)行每個(gè)子區(qū)的灌溉管理[23]。試驗(yàn)結(jié)果表明,由于黏土夾層對(duì)入滲水流的阻水作用[26],3區(qū)0.2~0.4 m土層土壤含水率較低,灌溉次數(shù)較多,難以實(shí)現(xiàn)基于灌水下限的水分虧缺管理。因此,2016—2017年冬小麥生育期內(nèi)(以下簡(jiǎn)稱(chēng)2017年冬小麥),除雨養(yǎng)處理外,各灌水處理灌水量設(shè)置為不同灌水定額(為設(shè)定灌水定額的百分比)(表1),采用脈沖式變量灌溉方法通過(guò)調(diào)節(jié)電磁閥占空比[23]進(jìn)行每個(gè)子區(qū)的灌溉管理。所有子區(qū)均采用相同的灌水時(shí)間,灌水下限值設(shè)計(jì)為AWC最小的1區(qū)田間持水量的70%。揚(yáng)花期之前灌水定額設(shè)定為20 mm,之后為30 mm。2 a雨養(yǎng)處理除了苗期灌水(20 mm)和施肥期間灌水(20 mm)外,其他時(shí)間均不灌溉。

    表1 田間試驗(yàn)處理設(shè)置

    注:水分虧缺程度設(shè)置為雨養(yǎng)、高、中和低4個(gè)水平。雨養(yǎng)處理僅在苗期和施肥期時(shí)進(jìn)行灌水,其他時(shí)間均不灌溉;2區(qū)低水分虧缺設(shè)置了2個(gè)處理,分別記為低1和低2。

    Note: A total of 4 water deficit levels include rain-fed, severe, medium, and low deficit levels in each management zone. There is no additional irrigation except for seedling emergence and fertilization stages in rain-fed treatments. In zone 2, 2 low deficit levels (Low 1 and Low 2) are included.

    將整個(gè)田塊劃分成12 m×12 m的網(wǎng)格,Trime-T3管式土壤水分探管(IMKO,德國(guó))埋設(shè)在每個(gè)網(wǎng)格的中心,1區(qū)、2區(qū)和3區(qū)的每個(gè)灌水處理內(nèi)分別布置6個(gè)、10個(gè)和3個(gè)探管,共86個(gè)。判斷灌水時(shí)間時(shí),2016年每日測(cè)量所有灌水處理0.2~0.4 m土層的土壤含水率,2017年每日測(cè)量R1-4處理0~0.4 m土層的平均土壤含水率,當(dāng)土壤含水率下降到設(shè)定下限時(shí)進(jìn)行灌水。為評(píng)估土壤含水率時(shí)間穩(wěn)定性,每7~10 d分層測(cè)量1次所有處理的土壤含水率,每0.2 m 1層,測(cè)量至1.2 m深度。2016年和2017年冬小麥生育期內(nèi)分別測(cè)量了10次和13次,共23次。因?yàn)?區(qū)每個(gè)子區(qū)內(nèi)布置的測(cè)量點(diǎn)數(shù)較少(3個(gè)),因此在土壤含水率時(shí)間穩(wěn)定性分析時(shí)未考慮3區(qū)。

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

    根據(jù)Spearman秩相關(guān)系數(shù)(r)原理[6]評(píng)價(jià)土壤含水率在不同時(shí)間的空間分布格局相似性:

    式中R是任一測(cè)點(diǎn)在第次觀(guān)測(cè)時(shí)土壤含水率測(cè)量值的秩,R'為測(cè)量值在相同測(cè)點(diǎn)、不同觀(guān)測(cè)次數(shù)時(shí)的秩。r越趨于1表示土壤含水率在2個(gè)測(cè)量時(shí)間下空間分布格局相似性越強(qiáng)。為測(cè)點(diǎn)數(shù),1區(qū)4個(gè)子區(qū)測(cè)點(diǎn)數(shù)量均為6個(gè),2區(qū)5個(gè)子區(qū)測(cè)點(diǎn)數(shù)量均為10個(gè)。

    計(jì)算每個(gè)測(cè)點(diǎn)的相對(duì)差分值[6],以判斷特定測(cè)點(diǎn)土壤含水率在時(shí)間尺度的穩(wěn)定性:

    式中為土壤含水率觀(guān)測(cè)次數(shù)。2016年和2017年冬小麥生育期內(nèi)分別監(jiān)測(cè)10次和13次。

    為了準(zhǔn)確捕捉不同采樣時(shí)間的土壤含水率均值,并有效地減少系統(tǒng)偏差,應(yīng)綜合考慮均值和方差,因此利用均方根誤差(root mean square error,RMSE)確定最佳的測(cè)點(diǎn)位置。

    RMSEi值最趨于0的測(cè)點(diǎn)可設(shè)置直接代表平均土壤含水率的傳感器。

    2 結(jié)果與分析

    2.1 不同灌溉制度下土壤含水率的空間結(jié)構(gòu)相似性

    2 a冬小麥生育期內(nèi),不同灌水處理的累積灌水量和降水量如圖1所示。2016年和2017年總降水量分別為80.1和80.3 mm,遠(yuǎn)小于華北地區(qū)冬小麥400 mm的平均需水量[27]。由于初始土壤含水率的空間變異性和管理區(qū)之間不同的土壤分層特征,2016年不同灌水處理累積灌水量并不總是隨灌水下限降低而減少,1區(qū)最大累積灌水量(380 mm)發(fā)生在F1-3處理,分別比雨養(yǎng)、F1-2和F1-4處理多334、80和40 mm;2區(qū)最大累積灌水量發(fā)生在F2-5處理(416 mm),分別比雨養(yǎng)、F2-3、F2-2和F2-4處理多370、190、110和60 mm。3區(qū)最大累積灌水量發(fā)生在F3-3處理(396 mm),分別比雨養(yǎng)、F3-1、F3-2處理多350、16和16 mm。與1區(qū)相比,2區(qū)(不包括F2-5)累積灌水量減少132 mm。2017年,1區(qū)、2區(qū)和3區(qū)相同處理的累積灌水量相同,雨養(yǎng)處理和3個(gè)灌水處理的累積灌水量分別為40、153、263和370 mm;R2-5處理累積灌水量為436 mm。

    圖1 冬小麥生育期內(nèi)各區(qū)灌水量和降水量

    表2 不同水分虧缺處理土壤含水率秩相關(guān)系數(shù)

    注:*和**分別表示在0.05和0.01水平上顯著。下同。

    Note: Symobol*and**mean significance at probability levels of 0.05 and 0.01, respectively. Same as below.

    2.2 不同測(cè)點(diǎn)土壤含水率的時(shí)間穩(wěn)定性

    每個(gè)測(cè)點(diǎn)的δ在時(shí)間尺度上的標(biāo)準(zhǔn)差σ如圖2所示。根據(jù)Starks等[29]的劃分標(biāo)準(zhǔn),當(dāng)σ<0.05時(shí),測(cè)點(diǎn)具有時(shí)間穩(wěn)定性。2016年1區(qū)所有處理0~0.2、>0.2~0.4和>0.4~0.6 m土層σ<0.05的測(cè)點(diǎn)數(shù)量占總測(cè)點(diǎn)的比例分別為83%,91%和65%;2區(qū)所有處理相應(yīng)土層該比例分別為85%,92%和80%。2017年1區(qū)相應(yīng)土層該比例分別為91%、91%和95%;2區(qū)分別為96%、100%和93%。這說(shuō)明大部分測(cè)點(diǎn)都具有較強(qiáng)的時(shí)間穩(wěn)定性,尤其是0~0.4 m土層,且除了2017年的>0.4~0.6 m土層外,2區(qū)任一土層具有時(shí)間穩(wěn)定性的測(cè)點(diǎn)數(shù)量均高于1區(qū)。

    圖2 各測(cè)點(diǎn)不同土層土壤含水率標(biāo)準(zhǔn)差

    每個(gè)測(cè)點(diǎn)均方根誤差RMSE值如圖3所示。在任一土層,所有測(cè)點(diǎn)的RMSE均在不同土層深度表現(xiàn)出明顯的波動(dòng),且>0.4~0.6 m土層RMSE波動(dòng)幅度最大,這表明各測(cè)點(diǎn)在剖面方向偏離平均土壤含水率的程度與深度有關(guān),代表平均土壤含水率的測(cè)點(diǎn)和深度具有不一致性。2 a同一測(cè)點(diǎn)的RMSE值不完全相等,2016年RMSE值整體上大于2017年,但各測(cè)點(diǎn)的變化趨勢(shì)基本保持一致,說(shuō)明各測(cè)點(diǎn)與平均土壤含水率的相對(duì)偏離程度具有年際間的時(shí)間穩(wěn)定性。

    將RMSE<0.1[30]的點(diǎn)作為代表含水率的點(diǎn)位,代表含水率的測(cè)點(diǎn)個(gè)數(shù)占總測(cè)點(diǎn)數(shù)的比例如表3所示。代表含水率的測(cè)點(diǎn)個(gè)數(shù)占總測(cè)點(diǎn)個(gè)數(shù)的比例與水分虧缺程度密切相關(guān)。2016年1區(qū)雨養(yǎng)和高、中、低水分虧缺處理0~0.6 m土層代表含水率的測(cè)點(diǎn)數(shù)量占總測(cè)點(diǎn)數(shù)量的比例分別為17%、22%、28%和33%;2區(qū)雨養(yǎng)和高、中、低水分虧缺處理0~0.6 m土層代表含水率的測(cè)點(diǎn)數(shù)量占總測(cè)點(diǎn)數(shù)量的比例分別為31%、33%、13%和38%(F2-4和F2-5處理均值)。2017年1區(qū)相應(yīng)水分虧缺處理該比例分別為39%、34%、44%和53%;2區(qū)則分別為52%、43%、43%和49%(R2-4和R2-5處理均值)。2 a試驗(yàn),除中水分虧缺處理外,2區(qū)任一水分虧缺處理代表含水率的測(cè)點(diǎn)比例均高于1區(qū)。2016年1區(qū)0~0.2、>0.2~0.4和>0.4~0.6 m土層代表含水率的測(cè)點(diǎn)比例分別為25%、46%和4%;2區(qū)0~0.2、>0.2~0.4和>0.4~0.6 m土層代表含水率的測(cè)點(diǎn)比例分別為36%、40%和17%。2017年1區(qū)和2區(qū)相應(yīng)土層代表含水率的測(cè)點(diǎn)分別為57%、44%、27%和72%、48%、22%,除1區(qū)>0.2~0.4 m土層外,該比例明顯高于2016年。這與2 a水分管理方法一致時(shí),代表含水率的測(cè)點(diǎn)個(gè)數(shù)占總測(cè)點(diǎn)個(gè)數(shù)的比例在年際間變化較小[23]的規(guī)律不同。另外,2 a試驗(yàn)1區(qū)和2區(qū)0~0.2、>0.2~0.4、>0.4~0.6 m土層代表含水率的測(cè)點(diǎn)比例平均值分別為41%、45%、16%和54%、44%、19%,具有不同土壤剖面特征的管理區(qū)同一土層代表含水率的測(cè)點(diǎn)比例差別較小的結(jié)果表明,代表含水率的測(cè)點(diǎn)個(gè)數(shù)占總測(cè)點(diǎn)個(gè)數(shù)的比例受土壤剖面的分層特征影響相對(duì)較小。

    圖3 各測(cè)點(diǎn)土壤含水率均方根誤差

    表3 不同土層代表含水率的測(cè)點(diǎn)與總測(cè)點(diǎn)數(shù)的比值

    2.3 代表平均土壤含水率點(diǎn)位的確定方法

    為研究土壤水分虧缺程度對(duì)代表含水率點(diǎn)位“平均屬性”的影響,分別將所有管理區(qū)雨養(yǎng)和高、中、低水分虧缺處理內(nèi)不同土層RMSE<0.1的測(cè)量點(diǎn)的平均黏粒含量作為自變量,將對(duì)應(yīng)土層內(nèi)代表含水率點(diǎn)位的黏粒含量當(dāng)作因變量,對(duì)二者進(jìn)行線(xiàn)性擬合(=,為擬合系數(shù)),結(jié)果如表4所示。

    表4 0~0.6 m土層含水率代表測(cè)點(diǎn)的黏粒體積分?jǐn)?shù)(y)與管理區(qū)平均黏粒體積分?jǐn)?shù)(x)的擬合關(guān)系

    注 Note:*,<0.05;**,<0.01。

    2016年雨養(yǎng)、中、低水分虧缺處理和2017年所有水分虧缺處理和之間均呈顯著的線(xiàn)性關(guān)系,2016年雨養(yǎng)和高水分虧缺處理值低于中、低水分虧缺處理,2017年值隨土壤水分虧缺程度增加呈增大趨勢(shì)。值隨水分虧缺程度的變化趨勢(shì)在年際間存在差異,這主要是因?yàn)橥寥婪謱犹卣鞯挠绊懯沟?016年各處理灌水量并未隨灌水下限的降低而減少(圖1)。2 a擬合系數(shù)的變化范圍為0.66~1.03,小于前期在該地塊不同管理區(qū)灌水下限均設(shè)置為0.45 AWC時(shí)得出的值(1.15)[11],這再次表明代表含水率的點(diǎn)位主要受土壤黏粒含量的影響,但與土壤黏粒含量的值在2017年隨土壤水分虧缺程度增加呈增大趨勢(shì)。

    3 結(jié) 論

    在華北平原土壤可利用水量變異較大的地塊,分別基于雨養(yǎng)和3種非充分灌溉制度開(kāi)展了為期2 a的冬小麥變量灌溉試驗(yàn),利用冬小麥生育期內(nèi)布設(shè)的86個(gè)測(cè)點(diǎn),共23次觀(guān)測(cè)值,分析了水分虧缺程度和土壤質(zhì)地對(duì)土壤含水率空間分布時(shí)間穩(wěn)定性的影響。主要結(jié)論如下:

    1)受水分虧缺程度和土壤剖面分層特征的影響,在砂粒含量隨土層深度增加的1區(qū),0~0.6 m土層的平均土壤含水率空間分布格局相似性隨土壤水分虧缺程度減小呈先減小后增加趨勢(shì),且相似性?xún)H在2016年的雨養(yǎng)處理達(dá)到了顯著水平。在沿垂直剖面方向砂粒含量分布均勻的2區(qū),0~0.6 m土層平均土壤含水率空間分布格局相似性雖然隨水分虧缺程度減小呈增加趨勢(shì),但2016年所有處理和2017年中、低水分虧缺處理相似性均達(dá)到了顯著水平。

    2)受水分虧缺程度的影響,2 a試驗(yàn)不同土層代表含水率的測(cè)點(diǎn)數(shù)量占總測(cè)點(diǎn)的比例差異較大。但受土壤剖面分層特征的影響較小,1區(qū)和2區(qū)同一土層代表含水率的測(cè)點(diǎn)比例平均值并無(wú)明顯差別,0~0.2、>0.2~0.4、>0.4~0.6 m土層該值分別為41%、45%、16%和54%、44%、19%。

    3)2 a試驗(yàn)中,除2016年高水分虧缺處理外,其他所有處理0~0.2、>0.2~0.4、>0.4~0.6 m土層所有測(cè)量點(diǎn)的平均黏粒含量與對(duì)應(yīng)土層內(nèi)代表含水率點(diǎn)位的黏粒含量之間均存在顯著的線(xiàn)性關(guān)系(<0.05),線(xiàn)性方程擬合系數(shù)為0.66~1.03,且2017年隨土壤水分虧缺程度增加呈增大趨勢(shì)。

    [1] Howell T A, Evett S R, O'Shaughnessy S A, et al. Advanced irrigation engineering: Precision and precise[J]. Journal of Agricultural Science and Technology, 2012, 29(1): 1-9.

    [2] O’Shaughnessy S A, Evett S R, Colaizzi P D. Dynamic prescription maps for site-specific variable rate irrigation of cotton[J]. Agricultural Water Management, 2015, 159: 123-138.

    [3] O'Shaughness S A, Rush C. Encyclopedia of Agriculture and Food Systems[M]. New York: Academic Press, 2014.

    [4] Evans R G, LaRue J, Stone K C, et al. Adoption of site-specific variable rate sprinkler irrigation systems[J]. Irrigation Science, 2013, 31(4): 871-887.

    [5] Sui Ruixiu, Baggard J. Wireless sensor network for monitoring soil moisture and weather conditions[J]. Applied Engineering in Agriculture, 2015, 31(2): 193-200.

    [6] Vachaud G, Passerat D S A, Balabanis P, et al. Temporal stability of spatially measured soil water probability density function[J]. Soil Science Society of American Journal, 1985, 49(4): 822-828.

    [7] Mohanty B P, Skaggs T H. Spatio-temporal evolution and time-stable characteristics of soil moisture within remote sensing footprints with varying soil, slope, and vegetation[J]. Advances in Water Resources, 2001, 24(9): 1051-1067.

    [8] Western A W, Zhou S L, Grayson R B, et al. Spatial correlation of soil moisture in small catchments and its relationship to dominant spatial hydrological processes[J]. Journal of Hydrology, 2004, 286(1): 113-134.

    [9] Ran Youhua, Li Xin, Jin Rui, et al. Strengths and weaknesses of temporal stability analysis for monitoring and estimating grid-mean soil moisture in a high-intensity irrigated agricultural landscape[J]. Water Resources Research, 2017, 53(1): 283-301.

    [10] Hu Wei, Shao Mingan, Han Fengpeng, et al. Watershed scale temporal stability of soil water content[J]. Geoderma, 2010, 158(3): 181-198.

    [11] 胡偉. 黃土高原小流域土壤含水量與飽和導(dǎo)水率的時(shí)空變異[D]. 北京: 中國(guó)科學(xué)院研究生院, 2009.

    Hu Wei. The Mean Clay Percentile was Used for a Priori Identification for Pplacement of Soil Moisture Sensors[D]. Beijing: Graduate University of Chinese Academy of Sciences, 2009.

    [12] Jacobs J M, Mohanty B P, Hsu E C, et al. SMEX02: Field scale variability, time stability and similarity of soil moisture [J]. Remote sensing of Environment, 2004, 92(4): 436-446.

    [13] Go'mez-Plaza A, Alvarez-Rogel J, Albaladejo J, et al. Spatial patterns and temporal stability of soil moisture across a range of scales in a semi-arid environment[J]. Hydrological Processes, 2000, 14(7): 1261-1277.

    [14] Li B, Rodell M. Spatial variability and its scale dependency of observed and modeled soil moisture over different climate regions[J]. Hydrology Earth System Sciences, 2013, 17(3): 1177-1188.

    [15] López-Vicente M, álvarez S. Stability and patterns of topsoil water content in rainfed vineyards, olive groves, and cereal fields under different soil and tillage conditions[J]. Agricultural Water Management, 2018, 201: 167-176.

    [16] Zhao Weixia, Li Jiusheng, Yang Rumiao, et al. Determinging placement criteria of moisture sensors through temporal stability analysis of soil water contents for a variable rate irrigation system[J]. Precision Agriculture, 2018, 19(4): 648-665.

    [17] Schneider K, Huisman J A, Breuer L, et al. Temporal stability of soil moisture in various semi-arid steppe ecosystems and its application in remote sensing[J]. Journal of Hydrology, 2008, 359(1): 16-29.

    [18] Zhao Long, Yang Kun, Qin Jun, et al. Spatio temporal analysis of soil moisture observations within a Tibetan mesoscale area and its implication to regional soil moisture measurements[J]. Journal of Hydrology, 2013, 482(5): 92-104.

    [19] Gao Lei, Shao Mingan. Temporal stability of soil water storage in diverse soil layers[J]. Catena, 2012, 95(95): 24-32.

    [20] Zhang Pingping, Shao Mingan, Zhang Xingchang. Temporal stability of soil moisture on two transects in a desert area of northwestern China[J]. Environmental Earth Sciences, 2016, 75(2): 161-175.

    [21] Grayson R B, Western A W. Towards areal estimation of soil water content from point measurements: Time and space stability of mean response[J]. Journal of Hydrology, 1998, 207(1): 68-82.

    [22] Evans R G, Sadler E J. Methods and technologies to improve efficiency of water use[J]. Water Resources Research, 2008, 44(7): 1-15.

    [23] 趙偉霞, 李久生, 楊汝苗, 等. 田間試驗(yàn)評(píng)估圓形噴灌機(jī)變量灌溉系統(tǒng)水量分布特性[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2014, 30(22): 53-62. Zhao Weixia, Li Jiusheng, Yang Rumiao, et al. Field evaluation of water distribution characteristics center pivot irrigation system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(22): 53-62. (in Chinese with English abstract)

    [24] 趙偉霞, 李久生, 楊汝苗, 等. 基于土壤水分空間變異的變量灌溉作物產(chǎn)量及節(jié)水效果[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2017, 33(2): 1-7. Zhao Weixia, Li Jiusheng, Yang Rumiao, et al. Yields and water-saving effects of crops as affected by variable rate irrigation management based on soil water spatial variation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(2): 1-7. (in Chinese with English abstract)

    [25] Zhao Weixia, Li Jiusheng, Yang Rumiao, 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.

    [26] 李久生, 楊風(fēng)艷, 栗巖峰. 層狀土壤質(zhì)地對(duì)地下滴灌水氮分布的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2009, 25(7): 25-31. Li Jiusheng, Yang Fengyan, Li Yanfeng. Water and nitrogen distribution under subsurface drip fertigation as affected by layered-textural soils[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(7): 25-31. (in Chinese with English abstract)

    [27] 胡瑋, 嚴(yán)昌榮, 李迎春, 等. 冀京津冬小麥灌溉需水量時(shí)空變化特征[J]. 中國(guó)農(nóng)業(yè)氣象, 2013, 34(6): 648-654. Hu Wei, Yan Changrong, Li Yingchun, et al. Spatial and temporal variation of irrigation water requirement for winter wheat in Jijingjin region[J]. Chinese Journal of Agrometeorology, 2013, 34(6): 648-654. (in Chinese with English abstract)

    [28] Van W B, Mulligan M, Poesen J. Spatial patterns of soil water balance on intensively cultivated hillslopes in a semiaridenvironment: The impact of rock fragments and soil thickness[J]. Hydrological Process, 2000, 14(10): 1811-1828.

    [29] Starks P J, Heathman G C, Jackson T J, et al. Temporal stability of soil moisture profile[J]. Journal of Hydrology, 2006, 324(1): 400-411.

    [30] Dettori M, Cesaraccio C, Motroni A, et al. Using CERES- Wheat to simulate durum wheat production and phenology in Southern Sardinia, Italy[J]. Field Crops Research, 2011, 120(1): 179-188.

    Influence of water stress level on determination of soil moisture sensor position under variable rate irrigation

    Li Xiumei, Zhao Weixia, Li Jiusheng※, Li Yanfeng

    (100048,)

    Determining positions to represent mean soil water content based on soil clay contents is an alternative method for positioning soil water content sensors in wireless sensor networks for a variable rate irrigation system. The field was divided into 4 management zones according to available soil water holding capacity (AWC). Two of the 4 management zones were selected to arrange the rain-fed treatment and 3 irrigation treatments representing different water stress levels to assess the effect of the levels on the placement of soil water content sensors under variable rate irrigation system. In zone 1, sand fraction largely increased with depth with AWC within 1-m soil profile ranging from 152 to 161 mm. In zone 2, a relatively uniform profile was observed along the profile with AWC within 1-m soil profile ranging from 161 to 171 mm. Based on the temporal stability of soil water content, the effects of soil water status and soil properties on the similarity of soil water spatial pattern and the positions directly representing the plot-mean soil water content were studied. Results showed that both of soil texture and water stress had effect on the structure similarity of soil water content distribution. In zone 1, the average Spearman’s rank correlation coefficient of 0-0.6 m was significant at the probability level of 0.05 only in the rain-fed treatment in the 2016 season. In zone 2, the Spearman’s rank correlation coefficient was significant at the probability level of 0.05 in all treatments in the 2016 season and in the medium and low water stress treatments in the 2017 season. The percentages of positions directly representing the mean soil water content were almost the same in zones 1 and 2. Affected by soil water status, the percentages increased as the level of soil water stress decreased in zone 1. While in zone 2, as the severity of water stress decreased, the percentages decreased and then had a slight increase. In general, significant linear regressions (0.05) between the mean clay content and the clay content representing the mean soil water content sites were found in layers 0-0.2, 0.2-0.4, and 0.4-0.6 m for all the treatments in 2016 and 2017, except for that in the severe water stress treatment in 2016. The fitted equation coefficients ranged from 0.66 to 1.03 in the 2 seasons, demonstrating a clearly increasing trend as the severity of water stress increased in 2017. When the mean clay content was used for a priori identification for positioning soil water content sensors in the management zones under variable rate irrigation system in a field with sandy loam soil, the strategies of water stress management should be considered in determining a fitted equation coefficients.

    soils; water content; sensors; networks; time stability; placement; winter wheat; variable rate irrigation

    李秀梅,趙偉霞,李久生,栗巖峰. 水分虧缺程度對(duì)變量灌溉水分傳感器埋設(shè)位置預(yù)判的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(23):94-100.doi:10.11975/j.issn.1002-6819.2018.23.011 http://www.tcsae.org

    Li Xiumei, Zhao Weixia, Li Jiusheng, Li Yanfeng. Influence of water stress level on determination of soil moisture sensor position under variable rate irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(23): 94-100. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.23.011 http://www.tcsae.org

    2018-06-24

    2018-08-10

    “十三五”國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0400104);中國(guó)水科院科研專(zhuān)項(xiàng)(2016TS05)共同資助

    李秀梅,河北邯鄲人,博士生,主要從事灌溉原理與技術(shù)方面的研究。Email:lixiumei11126@126.com

    李久生,河北邢臺(tái)人,研究員,博士生導(dǎo)師,主要從事灌溉原理與技術(shù)方面的研究。Email:lijs@iwhr.com

    10.11975/j.issn.1002-6819.2018.23.011

    S274.3

    A

    1002-6819(2018)-23-0094-07

    猜你喜歡
    黏粒灌水土層
    黏粒對(duì)紅黏土微觀(guān)結(jié)構(gòu)及力學(xué)性質(zhì)的影響
    土釘噴錨在不同土層的支護(hù)應(yīng)用及效果分析
    灌水取球
    番茄灌水掌握技巧
    冬季棚菜灌水四關(guān)鍵
    土層 村與人 下
    土層——伊當(dāng)灣志
    土層 沙與土 上
    不同黏粒含量黃土的人工切坡穩(wěn)定性探討
    黏粒含量對(duì)黃土物理力學(xué)性質(zhì)的影響
    国产 一区精品| 如何舔出高潮| 两个人看的免费小视频| av网站免费在线观看视频| 80岁老熟妇乱子伦牲交| 欧美日韩亚洲高清精品| 黄色毛片三级朝国网站| 菩萨蛮人人尽说江南好唐韦庄| 男的添女的下面高潮视频| 成年女人在线观看亚洲视频| 亚洲欧美中文字幕日韩二区| 男女国产视频网站| 在线看a的网站| 男女无遮挡免费网站观看| 午夜激情av网站| 日韩中字成人| 欧美av亚洲av综合av国产av | 两性夫妻黄色片| 国产免费又黄又爽又色| 777久久人妻少妇嫩草av网站| 在线亚洲精品国产二区图片欧美| 国产色婷婷99| 国产成人精品福利久久| 免费高清在线观看日韩| 99re6热这里在线精品视频| 久久亚洲国产成人精品v| 桃花免费在线播放| 69精品国产乱码久久久| 高清黄色对白视频在线免费看| kizo精华| 日韩精品免费视频一区二区三区| 国产深夜福利视频在线观看| 999久久久国产精品视频| 热re99久久精品国产66热6| 精品少妇一区二区三区视频日本电影 | 国产精品国产三级国产专区5o| 亚洲欧美中文字幕日韩二区| 午夜91福利影院| 在线观看免费日韩欧美大片| 在线免费观看不下载黄p国产| xxxhd国产人妻xxx| 欧美精品国产亚洲| 伦理电影大哥的女人| 欧美成人午夜免费资源| 9191精品国产免费久久| 午夜久久久在线观看| 免费播放大片免费观看视频在线观看| 如日韩欧美国产精品一区二区三区| 大话2 男鬼变身卡| 巨乳人妻的诱惑在线观看| 日本欧美视频一区| 日韩精品免费视频一区二区三区| 老汉色av国产亚洲站长工具| 美女大奶头黄色视频| 亚洲精品日韩在线中文字幕| 性少妇av在线| 欧美国产精品一级二级三级| 一二三四中文在线观看免费高清| 一二三四中文在线观看免费高清| 亚洲一区中文字幕在线| 韩国精品一区二区三区| 国产成人午夜福利电影在线观看| 美女高潮到喷水免费观看| av.在线天堂| 久久99一区二区三区| 婷婷色综合www| 午夜激情久久久久久久| 久久久久久伊人网av| 亚洲av福利一区| 欧美少妇被猛烈插入视频| 国产精品成人在线| 精品视频人人做人人爽| 国产xxxxx性猛交| videos熟女内射| 纯流量卡能插随身wifi吗| 久久国内精品自在自线图片| 国产综合精华液| 黄片播放在线免费| 伦理电影免费视频| 国产亚洲最大av| 亚洲精品国产一区二区精华液| 爱豆传媒免费全集在线观看| 99热网站在线观看| 爱豆传媒免费全集在线观看| 精品国产乱码久久久久久小说| 大话2 男鬼变身卡| 菩萨蛮人人尽说江南好唐韦庄| 免费少妇av软件| 精品少妇久久久久久888优播| 免费在线观看完整版高清| 18在线观看网站| 国产日韩欧美在线精品| 亚洲av综合色区一区| 一级毛片我不卡| av国产精品久久久久影院| 男女边吃奶边做爰视频| 一级片'在线观看视频| 欧美精品一区二区免费开放| 亚洲av电影在线进入| 丁香六月天网| 精品卡一卡二卡四卡免费| 男人爽女人下面视频在线观看| 亚洲综合色网址| 国产男人的电影天堂91| 女人精品久久久久毛片| 亚洲一级一片aⅴ在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品一区蜜桃| 高清欧美精品videossex| 王馨瑶露胸无遮挡在线观看| 91在线精品国自产拍蜜月| 免费观看无遮挡的男女| av在线播放精品| 王馨瑶露胸无遮挡在线观看| 日韩欧美精品免费久久| 不卡av一区二区三区| 男男h啪啪无遮挡| 亚洲成人手机| 日韩免费高清中文字幕av| 丝袜脚勾引网站| 春色校园在线视频观看| 丰满少妇做爰视频| 亚洲国产欧美网| 亚洲伊人久久精品综合| 亚洲av国产av综合av卡| 高清黄色对白视频在线免费看| 亚洲精品国产av蜜桃| 久久久久精品性色| 天天躁日日躁夜夜躁夜夜| 我的亚洲天堂| 伦理电影免费视频| 亚洲内射少妇av| 各种免费的搞黄视频| 一本—道久久a久久精品蜜桃钙片| 制服诱惑二区| 日韩一本色道免费dvd| 日本-黄色视频高清免费观看| 色94色欧美一区二区| h视频一区二区三区| 欧美av亚洲av综合av国产av | 国产免费现黄频在线看| 日韩电影二区| 国产在线视频一区二区| 国产高清不卡午夜福利| 人成视频在线观看免费观看| 最近最新中文字幕免费大全7| 成人18禁高潮啪啪吃奶动态图| 亚洲四区av| 国产精品久久久av美女十八| 久久久久国产网址| 亚洲熟女精品中文字幕| av天堂久久9| 最近手机中文字幕大全| 欧美精品亚洲一区二区| av在线老鸭窝| 黄频高清免费视频| 日本色播在线视频| 日韩成人av中文字幕在线观看| 精品人妻熟女毛片av久久网站| 成人漫画全彩无遮挡| 狂野欧美激情性bbbbbb| 在线观看一区二区三区激情| 亚洲精品美女久久久久99蜜臀 | 国产一区二区 视频在线| 国产白丝娇喘喷水9色精品| 天天躁夜夜躁狠狠久久av| 精品一区二区免费观看| 天美传媒精品一区二区| 99久久精品国产国产毛片| 国产成人a∨麻豆精品| 啦啦啦在线免费观看视频4| 天美传媒精品一区二区| 大片免费播放器 马上看| 亚洲四区av| 女人被躁到高潮嗷嗷叫费观| 999久久久国产精品视频| 欧美成人精品欧美一级黄| 女的被弄到高潮叫床怎么办| 啦啦啦啦在线视频资源| 丰满饥渴人妻一区二区三| 蜜桃国产av成人99| 国产深夜福利视频在线观看| 在线亚洲精品国产二区图片欧美| 亚洲av综合色区一区| 婷婷色av中文字幕| 欧美日韩亚洲高清精品| 久久久久久久大尺度免费视频| 色婷婷久久久亚洲欧美| 婷婷色综合www| 欧美人与性动交α欧美软件| 免费日韩欧美在线观看| 女性被躁到高潮视频| 美国免费a级毛片| 一区二区三区乱码不卡18| 最近2019中文字幕mv第一页| 亚洲三区欧美一区| 久久久久久久国产电影| 国产精品久久久av美女十八| 久久这里有精品视频免费| 成年美女黄网站色视频大全免费| 亚洲精品美女久久av网站| 男的添女的下面高潮视频| 久久精品亚洲av国产电影网| 女人精品久久久久毛片| 中国国产av一级| 国产在视频线精品| 黄色 视频免费看| 国产又爽黄色视频| 成年人免费黄色播放视频| 叶爱在线成人免费视频播放| 街头女战士在线观看网站| 欧美精品一区二区免费开放| 免费在线观看黄色视频的| 国产精品人妻久久久影院| 少妇猛男粗大的猛烈进出视频| 亚洲精品第二区| 日本91视频免费播放| 亚洲精品美女久久av网站| a级毛片在线看网站| 少妇人妻久久综合中文| 久久影院123| 久久久久久久精品精品| 久久午夜福利片| av福利片在线| 人妻人人澡人人爽人人| 可以免费在线观看a视频的电影网站 | 成人国产av品久久久| 国产精品蜜桃在线观看| 亚洲av免费高清在线观看| 一二三四在线观看免费中文在| 91成人精品电影| 精品少妇黑人巨大在线播放| 天天躁夜夜躁狠狠久久av| 国产精品亚洲av一区麻豆 | 国产男女内射视频| 99久久综合免费| 久热久热在线精品观看| 赤兔流量卡办理| 日本免费在线观看一区| 欧美日韩视频精品一区| 亚洲精品国产色婷婷电影| 不卡av一区二区三区| 午夜福利视频在线观看免费| 叶爱在线成人免费视频播放| 老司机影院毛片| 日韩伦理黄色片| 国产成人免费无遮挡视频| 免费观看性生交大片5| 日本爱情动作片www.在线观看| 久久人人爽av亚洲精品天堂| 亚洲欧美一区二区三区久久| 久久国产亚洲av麻豆专区| 在线免费观看不下载黄p国产| 久久久久精品性色| 久久人人97超碰香蕉20202| 老熟女久久久| 黄色一级大片看看| 精品福利永久在线观看| 99久久中文字幕三级久久日本| 人人妻人人添人人爽欧美一区卜| a 毛片基地| av在线老鸭窝| 亚洲av电影在线进入| 看非洲黑人一级黄片| 久久热在线av| 久久精品夜色国产| 日本欧美视频一区| 免费高清在线观看日韩| av免费观看日本| 七月丁香在线播放| 精品少妇久久久久久888优播| 久久午夜福利片| 免费日韩欧美在线观看| 久久久久久免费高清国产稀缺| 免费少妇av软件| 黑丝袜美女国产一区| 男人操女人黄网站| 国产一区亚洲一区在线观看| 亚洲精品久久午夜乱码| 亚洲精品日本国产第一区| 晚上一个人看的免费电影| 国产淫语在线视频| 伦理电影大哥的女人| 18在线观看网站| 大码成人一级视频| 美女福利国产在线| 丝袜脚勾引网站| 精品亚洲乱码少妇综合久久| 考比视频在线观看| 老司机影院毛片| 久久免费观看电影| 国产精品亚洲av一区麻豆 | 国产精品国产三级专区第一集| 免费不卡的大黄色大毛片视频在线观看| 精品卡一卡二卡四卡免费| 九九爱精品视频在线观看| 99re6热这里在线精品视频| 1024视频免费在线观看| 男女边吃奶边做爰视频| 最近2019中文字幕mv第一页| 亚洲精品久久成人aⅴ小说| 在线看a的网站| 这个男人来自地球电影免费观看 | 国产有黄有色有爽视频| 99热国产这里只有精品6| 日本爱情动作片www.在线观看| 久久久精品区二区三区| av视频免费观看在线观看| 叶爱在线成人免费视频播放| 亚洲美女视频黄频| 亚洲婷婷狠狠爱综合网| 欧美日韩亚洲高清精品| 国产欧美亚洲国产| 少妇人妻 视频| 久久韩国三级中文字幕| 久久久久久人妻| 国产一区二区激情短视频 | 精品国产乱码久久久久久小说| 中文字幕av电影在线播放| 中文字幕亚洲精品专区| 99re6热这里在线精品视频| 人人妻人人澡人人爽人人夜夜| 女人精品久久久久毛片| 有码 亚洲区| 欧美日韩av久久| 丝袜在线中文字幕| 亚洲国产精品999| 亚洲成色77777| 看非洲黑人一级黄片| 久久99热这里只频精品6学生| 久久久国产欧美日韩av| 丝袜美腿诱惑在线| 亚洲精品国产一区二区精华液| 啦啦啦啦在线视频资源| 国产精品.久久久| 一本—道久久a久久精品蜜桃钙片| 国产精品嫩草影院av在线观看| 男人添女人高潮全过程视频| 国产精品 国内视频| 国产日韩一区二区三区精品不卡| 观看av在线不卡| 成年美女黄网站色视频大全免费| 日本色播在线视频| 一二三四中文在线观看免费高清| 韩国精品一区二区三区| 国产伦理片在线播放av一区| 日本午夜av视频| 91精品伊人久久大香线蕉| 三级国产精品片| 看免费成人av毛片| 精品一区二区免费观看| av电影中文网址| 免费少妇av软件| 日韩av不卡免费在线播放| 在线观看国产h片| 欧美 日韩 精品 国产| 亚洲成色77777| 成年女人毛片免费观看观看9 | 国产精品国产三级专区第一集| 爱豆传媒免费全集在线观看| 我要看黄色一级片免费的| 国产精品久久久av美女十八| 人人妻人人澡人人爽人人夜夜| 国产精品亚洲av一区麻豆 | 亚洲欧洲精品一区二区精品久久久 | 欧美黄色片欧美黄色片| 26uuu在线亚洲综合色| 一本大道久久a久久精品| h视频一区二区三区| 亚洲av中文av极速乱| 热re99久久国产66热| 国产精品香港三级国产av潘金莲 | 国产人伦9x9x在线观看 | 欧美精品一区二区大全| 有码 亚洲区| 亚洲欧美一区二区三区国产| 国产伦理片在线播放av一区| 两性夫妻黄色片| 日韩成人av中文字幕在线观看| 国产黄色视频一区二区在线观看| 精品亚洲乱码少妇综合久久| 国产高清不卡午夜福利| 日韩大片免费观看网站| 老司机影院成人| 成年动漫av网址| 免费人妻精品一区二区三区视频| 女人精品久久久久毛片| 美女国产高潮福利片在线看| a级片在线免费高清观看视频| 观看av在线不卡| 午夜91福利影院| 免费看不卡的av| 日本午夜av视频| 国产成人a∨麻豆精品| 九九爱精品视频在线观看| 女人久久www免费人成看片| 日本vs欧美在线观看视频| av卡一久久| 亚洲伊人色综图| 久久精品国产亚洲av高清一级| 精品一区二区三区四区五区乱码 | 国产精品一二三区在线看| 国产在线免费精品| av免费在线看不卡| 午夜激情av网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品亚洲av一区麻豆 | av电影中文网址| 1024香蕉在线观看| 国产麻豆69| 男人爽女人下面视频在线观看| 老熟女久久久| 国产黄频视频在线观看| 2021少妇久久久久久久久久久| 久久这里有精品视频免费| 欧美人与性动交α欧美精品济南到 | 丰满乱子伦码专区| 国产精品蜜桃在线观看| 亚洲第一区二区三区不卡| 国产日韩一区二区三区精品不卡| 国产爽快片一区二区三区| 最近最新中文字幕大全免费视频 | 成年av动漫网址| 亚洲精品美女久久久久99蜜臀 | 久久久久网色| 亚洲欧美色中文字幕在线| 亚洲一区二区三区欧美精品| 国产亚洲欧美精品永久| 国产乱来视频区| 午夜免费男女啪啪视频观看| 精品人妻偷拍中文字幕| 国产亚洲欧美精品永久| 如何舔出高潮| 色94色欧美一区二区| 久久久国产一区二区| 久久精品夜色国产| 90打野战视频偷拍视频| 日韩av免费高清视频| 自线自在国产av| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 精品国产一区二区三区四区第35| 亚洲色图 男人天堂 中文字幕| 日韩欧美精品免费久久| 18+在线观看网站| 国产精品欧美亚洲77777| 国产老妇伦熟女老妇高清| 欧美变态另类bdsm刘玥| av福利片在线| 一级a爱视频在线免费观看| 深夜精品福利| 精品少妇一区二区三区视频日本电影 | 人妻一区二区av| 久久精品国产自在天天线| 高清欧美精品videossex| 黄色一级大片看看| av网站在线播放免费| 女的被弄到高潮叫床怎么办| 满18在线观看网站| 黑人猛操日本美女一级片| 亚洲图色成人| 欧美日韩一区二区视频在线观看视频在线| 99香蕉大伊视频| 亚洲精品av麻豆狂野| 亚洲美女搞黄在线观看| 青春草国产在线视频| 国产又色又爽无遮挡免| 观看av在线不卡| 亚洲欧美一区二区三区国产| 亚洲国产精品999| 超色免费av| 毛片一级片免费看久久久久| 考比视频在线观看| 一本色道久久久久久精品综合| 亚洲成人一二三区av| 国产av码专区亚洲av| 日本av手机在线免费观看| 日韩视频在线欧美| 免费播放大片免费观看视频在线观看| a级毛片黄视频| 亚洲经典国产精华液单| 国产日韩一区二区三区精品不卡| 亚洲伊人色综图| 国产黄色免费在线视频| 国产精品熟女久久久久浪| 99re6热这里在线精品视频| 一本—道久久a久久精品蜜桃钙片| 国产亚洲一区二区精品| 黑丝袜美女国产一区| 亚洲情色 制服丝袜| 亚洲精品一二三| 国产精品 欧美亚洲| 极品少妇高潮喷水抽搐| 少妇人妻 视频| 男人舔女人的私密视频| 午夜日本视频在线| 色播在线永久视频| 日韩一区二区视频免费看| 成人国产av品久久久| 亚洲精品第二区| 90打野战视频偷拍视频| 国产成人精品婷婷| 亚洲国产成人一精品久久久| 熟女电影av网| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品美女久久久久99蜜臀 | 在线亚洲精品国产二区图片欧美| 男女免费视频国产| 日日摸夜夜添夜夜爱| 亚洲第一av免费看| 国产一区二区在线观看av| 国产精品久久久久久久久免| 亚洲精品视频女| 热99国产精品久久久久久7| 欧美最新免费一区二区三区| 一区二区三区四区激情视频| h视频一区二区三区| 可以免费在线观看a视频的电影网站 | 男女免费视频国产| 国产精品欧美亚洲77777| 欧美成人午夜精品| 两个人免费观看高清视频| 母亲3免费完整高清在线观看 | 色哟哟·www| 免费黄网站久久成人精品| 最近最新中文字幕大全免费视频 | 欧美成人精品欧美一级黄| 精品国产露脸久久av麻豆| 久久精品久久精品一区二区三区| 午夜免费观看性视频| 国产精品亚洲av一区麻豆 | 99热全是精品| 亚洲精品aⅴ在线观看| 丰满少妇做爰视频| 黄色配什么色好看| 国产一区二区 视频在线| 乱人伦中国视频| 一级,二级,三级黄色视频| 久久国内精品自在自线图片| 久久久欧美国产精品| 国产男人的电影天堂91| 日本91视频免费播放| 日本wwww免费看| 女的被弄到高潮叫床怎么办| 91aial.com中文字幕在线观看| 一级a爱视频在线免费观看| 边亲边吃奶的免费视频| 99久国产av精品国产电影| 熟女电影av网| 精品国产国语对白av| 色94色欧美一区二区| 日韩熟女老妇一区二区性免费视频| 欧美国产精品一级二级三级| 午夜免费观看性视频| 夫妻午夜视频| 久久久久久久亚洲中文字幕| 亚洲国产av影院在线观看| 人妻系列 视频| 交换朋友夫妻互换小说| 五月天丁香电影| 在线看a的网站| 女性被躁到高潮视频| 制服丝袜香蕉在线| 欧美国产精品va在线观看不卡| 人体艺术视频欧美日本| 亚洲成人手机| 国产男女超爽视频在线观看| 嫩草影院入口| 国产成人精品福利久久| 久久青草综合色| 捣出白浆h1v1| 综合色丁香网| 久久久久精品人妻al黑| 亚洲精品久久午夜乱码| 久久鲁丝午夜福利片| 免费观看无遮挡的男女| 黄色一级大片看看| 国产精品久久久久久av不卡| 久久久国产精品麻豆| 啦啦啦视频在线资源免费观看| 欧美日本中文国产一区发布| 久久青草综合色| 黄色配什么色好看| 一二三四中文在线观看免费高清| 久久久久国产网址| 久久韩国三级中文字幕| 久久久久久久久久久免费av| 大片免费播放器 马上看| 国产精品av久久久久免费| 春色校园在线视频观看| 日本免费在线观看一区| 久久精品国产综合久久久| 在线天堂中文资源库| 午夜免费鲁丝| 中文乱码字字幕精品一区二区三区| 精品亚洲成国产av| av在线app专区| 在线精品无人区一区二区三| 国产精品偷伦视频观看了| 三级国产精品片| 只有这里有精品99| 女人被躁到高潮嗷嗷叫费观| av免费观看日本| 欧美 日韩 精品 国产| 91aial.com中文字幕在线观看| 精品国产超薄肉色丝袜足j| 超色免费av| 国产精品一二三区在线看| 国产极品天堂在线| 少妇熟女欧美另类| 国产精品国产三级专区第一集| 99热全是精品| 制服诱惑二区| 如何舔出高潮| 中文字幕人妻熟女乱码|