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

    日光溫室番茄不同空間尺度蒸散量變化及主控因子分析

    2017-05-25 08:00:37龔雪文孫景生馬筱建王萬寧
    農(nóng)業(yè)工程學(xué)報 2017年8期

    龔雪文,劉 浩,孫景生※,張 昊,馬筱建,王萬寧

    ?

    日光溫室番茄不同空間尺度蒸散量變化及主控因子分析

    龔雪文1,2,劉 浩1,孫景生1※,張 昊1,2,馬筱建1,2,王萬寧1,2

    (1. 中國農(nóng)業(yè)科學(xué)院農(nóng)田灌溉研究所,農(nóng)業(yè)部作物需水與調(diào)控重點開放實驗室,新鄉(xiāng) 453003;2. 中國農(nóng)業(yè)科學(xué)院研究生院,北京 100081)

    明確日光溫室作物不同空間尺度蒸散量及變化規(guī)律是提高水分利用效率、實現(xiàn)農(nóng)業(yè)水資源合理配置的關(guān)鍵。該文針對華北地區(qū)典型日光溫室,于2015—2016年在中國農(nóng)業(yè)科學(xué)院新鄉(xiāng)綜合試驗基地,以滴灌番茄為研究對象,參考20 cm標準蒸發(fā)皿的累積蒸發(fā)量,設(shè)計充分灌溉和虧缺灌溉2種水平,研究不同水平下番茄葉片蒸騰、單株耗水(用莖流速率表征)和群體蒸散量的日變化和生育期變化,并采用通徑分析法確定影響不同空間尺度蒸散量的主控因子。結(jié)果表明:葉片蒸騰和氣孔導(dǎo)度隨太陽輻射變化,峰值出現(xiàn)在10:00—14:00之間,移栽54~58 d后充分和虧缺處理的葉片蒸騰和氣孔導(dǎo)度開始出現(xiàn)差異;充分和虧缺處理的單株耗水在晴天差異最大,陰雨天最小,且滯后太陽輻射約1 h;全生育期充分和虧缺處理的日群體蒸散量分別在0.32~6.65和0.15~5.91 mm/d之間變化,群體蒸散量在盛果期最大,占總耗水量的31.7%~34.7%。凈輻射對葉片、單株和群體尺度的蒸騰量影響均顯著,而水汽壓差僅對單株和群體尺度蒸散量影響顯著,估算日光溫室番茄單株耗水和群體蒸散量時需考慮風(fēng)速影響。水分脅迫條件下,考慮葉溫變量可顯著提高單株耗水和群體蒸散量的估算精度。研究可為不同空間尺度蒸散量轉(zhuǎn)換方法的選擇以及尺度提升理論模型的構(gòu)建提供借鑒。

    蒸散;蒸騰;通徑分析;氣孔導(dǎo)度;莖流速率

    0 引 言

    近年來,多尺度蒸散量之間如何轉(zhuǎn)換與提升越來越受到關(guān)注[1-3]。葉片、單株和群體尺度是3個不同空間尺度上表觀作物水分利用效率的重要指標,且不同尺度蒸散量所受到的環(huán)境影響因子也不相同[4]。探索多尺度蒸散量變化規(guī)律,明確各尺度環(huán)境主控因子是進一步評價農(nóng)業(yè)用水管理的基礎(chǔ)。

    日光溫室是一個封閉半封閉獨立結(jié)構(gòu)體,內(nèi)部能量平衡和水汽平衡自成體系,該系統(tǒng)與多尺度蒸散量變化息息相關(guān)。目前,基于溫室環(huán)境下的相關(guān)報道多集中在葉片蒸騰[5-6]、單株耗水[7-8]和蒸散量[9-10]的單尺度方面,由于不同空間尺度蒸散量之間既有差異又存在一定關(guān)聯(lián)性,是一個復(fù)雜多變的過程,尤其在水分脅迫條件下,各尺度環(huán)境影響因子差異更大,而有關(guān)這方面的研究報道還比較少見。在多尺度問題上,蔡甲冰等[11-12]對冬小麥返青后蒸騰量的尺度效應(yīng)及轉(zhuǎn)換關(guān)系進行了研究,并采用多元線性回歸方法進行了尺度提升;Zhang等[4]針對夏玉米多尺度蒸散量變化以及環(huán)境控制因子進行了討論,并利用通徑分析法確定了各尺度蒸散量的主控因子;徐璐等[13]確定了典型荒漠區(qū)2種植物多尺度蒸騰的主控因子分別為太陽輻射和水汽壓差。然而,針對日光溫室內(nèi)部高溫高濕和低風(fēng)速特點,以上結(jié)論可能有所不同,且采用線性回歸模型進行尺度提升受當?shù)丨h(huán)境條件的限制。

    研究表明,氣象因子和土壤水分狀況是影響作物耗水的主要因素,其中太陽輻射是驅(qū)動作物蒸騰的主要驅(qū)動力[4,14],而葉片周圍的水汽壓差和葉溫變化也會間接影響蒸騰量[15]。Wilson等[16]認為植株莖稈液流與冠層內(nèi)部的水汽壓差一致性較高,且無水分脅迫條件下的相關(guān)性最好。土壤水分影響農(nóng)田耗水主要通過改變作物根區(qū)土壤水勢實現(xiàn)的,水分脅迫條件下土壤水勢低,根系吸水阻力大,蒸騰速率下降。陳倩倩等[17]認為輕度和中度土壤水分脅迫可降低全生育期單株玉米蒸騰量的16.3%和36.9%;姚勇哲等[18]認為土壤水分狀況是影響番茄蒸騰量的主要決策因子,其決策系數(shù)為27.4%;大量研究表明,水分脅迫對作物蒸騰和農(nóng)田耗水具有抑制作用[19-20]。然而,充分和虧缺供水條件下,對不同空間尺度蒸散量變化規(guī)律以及不同尺度主控環(huán)境因子的研究報道還比較少見,這些關(guān)鍵性指標是用于制定溫室環(huán)境閾值,指導(dǎo)農(nóng)業(yè)用水的理論基礎(chǔ)。

    綜上,為研究日光溫室作物不同空間尺度蒸散量及變化規(guī)律,以華北地區(qū)典型日光溫室番茄為研究對象,于2015—2016年對充分和虧缺2種水分條件下的葉片蒸騰、單株耗水和群體蒸散量進行系統(tǒng)研究,利用通徑分析法確定影響不同尺度蒸散量的主控環(huán)境因子,旨在為進一步研究多尺度蒸散量轉(zhuǎn)換方法以及構(gòu)建尺度提升理論模型提供借鑒。

    1 材料與方法

    1.1 試驗區(qū)概況

    于2015和2016年的3—6月在中國農(nóng)業(yè)科學(xué)院新鄉(xiāng)綜合試驗基地的日光溫室內(nèi)進行,該試驗站位于新鄉(xiāng)縣七里營鎮(zhèn)(35°9¢N、113°5¢E,海拔78.7 m),屬暖溫帶大陸性季風(fēng)氣候區(qū),多年平均降雨量為540 mm,蒸發(fā)量在1 910 mm左右,全年平均氣溫為14.5 ℃,年日照時數(shù)為2 395 h,無霜期200 d。試驗所用日光溫室坐北朝南,占地510 m2(60 m × 8.5 m),下沉0.5 m,為鋼架結(jié)構(gòu),覆蓋聚乙烯薄膜。日光溫室后墻和山墻內(nèi)鑲嵌有60 cm厚的保溫材料,頂部用5 cm厚的保溫棉被覆蓋,以確保室內(nèi)溫度控制在10~35 ℃范圍內(nèi)。試驗區(qū)0~80 cm土層的土壤質(zhì)地為壤土,每20 cm為1層,剖面土壤容重依次為1.47、1.44、1.52和1.54 g/cm3,田間持水量(體積分數(shù),下同)依次為0.31、0.28、0.32和0.38 cm3/cm3,80~100 cm為砂壤土,容重和田間持水量分別為1.46 g/cm3和0.33 cm3/cm3。

    1.2 試驗設(shè)計及過程

    2 a番茄品種均為“金頂新星”,分別于2015年1月5日和2016年1月6日育苗,2015年3月8日和2016年3月9日移栽,采用滴灌方式灌水(滴頭間距為33 cm,滴頭流量為1.1 L/h),確保滴頭與植株一一對應(yīng),采用精度為0.001 m3的水表控制各試驗區(qū)灌水量。采用20 cm標準蒸發(fā)皿累計蒸發(fā)量(E)控制每次灌水量,當E達到(20 ± 2)mm時進行灌水。參考劉浩[21]對日光溫室滴灌番茄多年試驗研究成果,試驗設(shè)計0.9E為充分供水(T0.9)和0.5E為虧缺供水(T0.5)2種灌溉水平,試驗小區(qū)畦長8.0 m,寬1.1 m,采用雙行種植模式,行距45 cm,株距33 cm。每個處理4個重復(fù),每個重復(fù)50株,試驗小區(qū)采用完全隨機區(qū)組排列。各小區(qū)之間埋設(shè)60 cm深的塑料薄膜,以防止水分側(cè)滲。幼苗移栽后,為確保幼苗成活率并加強其長勢,所有小區(qū)均以滴灌方式灌水 20 mm。2 a不同水分處理的灌水時間和灌水量見表1。苗期不做水分處理,且灌水時間和灌水量一致,試驗起止時間分別為2015年4月24日—2015年7月15日和2016年4月19日—2016年7月10日。番茄各生育期的時間段分別為:苗期(2015年3月10日—4月15日、2016年3月9日—4月10日)、開花坐果期(2015年4月16日—5月15日、2016年4月11日—5月10日)、盛果期(2015年5月16日—6月15日、2016年5月11日—6月10日)和采摘期(2015年6月16日—7月15日、2016年6月11日—7月10日)。

    表1 2015和2016年不同水分處理灌水時間和灌水量

    注:T0.9和T0.5分別表示充分供水(0.9E)和虧缺供水(0.5E);E為20 cm標準蒸發(fā)皿累計蒸發(fā)量,下同。

    Note: T0.9and T0.5were full irrigation (0.9E) and deficit irrigation (0.5E), respectively;Eis accumulated evaporation in 20-cm standard evaporation pan; Same as below.

    1.3 試驗觀測項目

    1.3.1 土壤水分

    采用TRIME-IPH時域反射儀(Micromodultechnik GmbH,Germany)測定0~100 cm土壤層含水率,為確保儀器測定的準確性,定期采用取土烘干法對儀器進行矯正。滴灌為局部灌溉,采用劉浩[21]確定的最能反應(yīng)根區(qū)水分狀況的位置進行測量,即同1條滴灌帶的2個滴頭中間位置。由于幼苗期不做水分處理,這期間,含水率每7~10 d測1次,之后每隔7 d測1次。

    1.3.2 葉片蒸騰速率

    采用LI-6400XT便攜式光合作用測量系統(tǒng)(LI-COR,USA),在番茄每個生育期選擇晴朗無云日測定葉片蒸騰速率(T)和氣孔導(dǎo)度(g)的日變化過程,測量時間范圍在7:00—18:00之間,每隔1~1.5 h測定1次。對于生育期變化,每隔7~10 d,選擇11:00—11:30之間測定葉片的Tg等指標。在測量時,隨機選取形態(tài)均勻無病蟲害的番茄植株,選取葉片生長點以下第3片功能葉測定其中部位置,并保持所測葉片與太陽光垂直以便葉片得到充分光誘導(dǎo),每個處理測定12片。為避免因環(huán)境條件改變而引起的測量誤差,每次測量過程在30 min以內(nèi)完成。

    1.3.3 單株蒸騰速率

    作物蒸騰量與莖稈液流量密切相關(guān),測定植株莖稈液流是獲取作物蒸騰量的有效方法之一[7-8]。采用包裹式莖流計(Flow32-1k system,Dynamax,Houston,TX,USA)于2015年5月15日—7月15日和2016年4月18日—7月10日測定番茄植株莖桿液流速率。每個處理隨機挑選6株形態(tài)均勻無病蟲害的植株,探頭安裝在第3枝節(jié)與第4枝節(jié)之間,并與地表保持20 cm高度,以避免土壤熱反射的干擾。為確保探頭與莖稈緊密接觸,所選探頭尺寸規(guī)格需滿足番茄莖稈直徑要求。采用CR1000數(shù)據(jù)采集器(Campbell,USA)每15 min自動記錄1次數(shù)據(jù)。其標準化處理方法見文獻[7]。

    1.3.4 群體蒸散量

    于2016年采用2臺稱重式蒸滲儀(長×寬×高=1 m× 1 m × 1.2 m)分別測定T0.9和T0.5的群體蒸散量(evapotranspiration,ET),蒸滲儀安裝在小區(qū)中部,且1.0 m深度內(nèi)的土壤質(zhì)地與田間相同,移栽時挑選6株長勢均勻無病蟲害的幼苗定植于蒸滲儀內(nèi),其間距和行距與大田相同,待植株長至40 cm高時做搭架處理。蒸滲儀內(nèi)番茄耗水量采用水量平衡法計算,即

    ·=(-1-)/+(1)

    式中為時間段內(nèi)農(nóng)田耗水量,mm;為蒸滲儀表面積,mm2;-1和分別為-1時刻和時刻蒸滲儀內(nèi)土體質(zhì)量,g;為水的密度,1.0 g/cm3;I為時段內(nèi)進入到蒸滲儀內(nèi)的水量,mm3。

    1.3.5 氣象數(shù)據(jù)

    日光溫室中部安裝有1套氣象自動監(jiān)測系統(tǒng),該系統(tǒng)包括凈輻射(NR LITE2, Kipp & Zonen, Delft, Netherlands)、總輻射(LI200X, Campbell Scientific, Inc, USA)、溫濕度(CS215, Campbell Scientific, Inc, USA)、風(fēng)速傳感器(Wind Sonic, Gill, UK)。冠層溫度用位于冠層上方30 cm處的紅外測溫儀測定(SI-111, Campbell Scientific, Inc, USA),探頭朝北,精度為±0.2 ℃。所有氣象數(shù)據(jù)每隔30 min由CR1000數(shù)據(jù)采集系統(tǒng)(Campbell,USA)自動記錄。表2列出了2015和2016年試驗期間日光溫室番茄不同生育階段的環(huán)境特征,其中參考作物蒸散量(reference evapotranspiration,ET0)選用陳新明等[22]推薦的方法獲得。采用20 cm標準蒸發(fā)皿(直徑20 cm,深11 cm)測定水面蒸發(fā)量,于每日7:30— 8:00之間進行測量,蒸發(fā)皿至于冠層上方30 cm處,并隨著冠層高度進行調(diào)整,每次測量完畢后更換蒸發(fā)皿中的蒸餾水,并重新添加到20 mm水位處。2 a日光溫室氣象因子平均值如表2所示。

    1.3.6 通徑分析法

    通徑分析法是基于結(jié)構(gòu)方程建模,包括統(tǒng)計技術(shù),如影響因子分析和線性回歸分析,并將其集成到多元統(tǒng)計分析中。這種方法是用于識別、估計和驗證所有類型的因果模型,以及確定變量之間的相互關(guān)系、行動模式以及從獨立到從屬變量的影響程度。該方法闡明了所有與因變量有關(guān)的直接影響和間接影響,從獨立因變量與多重自變量之間確定多元線性方程[23]。本文中用于通徑分析的數(shù)據(jù)均為小時尺度數(shù)據(jù),葉片蒸騰速率T為2015年和2016年日變化數(shù)據(jù),為明確不同水分條件下單株莖流速率S和群體蒸散量ET主控環(huán)境因子的差異性,特選用具有明顯水分差異且連續(xù)觀測的試驗數(shù)據(jù)(2016年6月1日—7月10日)進行分析,以上所選環(huán)境參數(shù)均為測量TS和ET時所對應(yīng)的小時尺度數(shù)據(jù)。

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

    采用Excel 2010分析數(shù)據(jù)并繪圖,采用SPSS 9.0進行通徑分析、顯著性分析和模型構(gòu)建。

    表2 2015和2016年日光溫室番茄不同生育階段氣象因子平均值

    2 結(jié)果與分析

    2.1 不同水分處理根系層土壤水分分布特征

    研究表明,番茄80%以上根系分布在0~60 cm土壤層[24],該土壤層水分分布特征直接影響著根系生長以及對土壤水分的利用,進而影響葉片、植株乃至群體尺度的耗水特性。圖1為2種不同水分條件下根系層平均含水率隨生育期的變化過程??梢钥闯?,移栽54 d后不同水分處理的土壤含水率開始出現(xiàn)明顯差異,T0.9處理的根系層含水率高于T0.5處理。從整個生育期的變化過程來看,2015年番茄移栽54 d后,2016年移栽58 d后,T0.9與T0.5處理的差異開始逐漸明顯,但這種差異隨著生育期的推進并無明顯增大跡象,說明0.9E灌溉水平在滿足作物耗水的同時,能確保根系層土壤含水率維持在一個相對較高水平,而0.5E灌溉水平,由于根系層土壤含水率相對較低(平均含水率體積分數(shù)在0.19~0.20 cm3/cm3之間),導(dǎo)致根系吸水困難,抑制了番茄的耗水過程。從土壤水分出現(xiàn)差異到試驗結(jié)束期間,2015和2016年T0.9的平均含水率體積分數(shù)分別為0.25和0.24 cm3/cm3,T0.5分別為0.20和0.19 cm3/cm3,2 a T0.9的平均含水率分別較T0.5高20.5%和21.7%(<0.05)。

    圖1 2015和2016年不同處理日光溫室番茄根系層含水率變化

    2.2 不同水分處理番茄葉片蒸騰速率和氣孔導(dǎo)度變化

    2.2.1 日變化規(guī)律

    圖2給出了2015和2016年T0.9和T0.5處理在不同生育階段葉片蒸騰速率(T)和氣孔導(dǎo)度(g)的日變化規(guī)律。可以看出,Tg的日變化均表現(xiàn)為開口向下的單峰曲線,最大值出現(xiàn)在10:00—14:00之間。如上所述,2015、2016年T0.9和T0.5根層含水率分別在移栽54 d和58 d后開始出現(xiàn)差異,從圖2可看出,2015年4月30日和2016年5月3日T0.9與T0.5的Tg并無顯著差異(>0.05)。一般情況下,高水處理的Tg高于低水處理[21],本研究得出類似結(jié)論,如2015年5月24日和6月22日T0.9處理的日平均Tg分別較T0.5高15.9%和30.9%、14.9%和19.9% (<0.05)。

    注:所有單位以H2O計。

    2.2.2 生育期動態(tài)變化

    圖3給出了不同水分處理番茄葉片Tg在整個生育期內(nèi)的動態(tài)變化,可以看出,Tg在苗期小,盛果期增大而采摘期又減小的變化趨勢,2016年T0.9處理分別在移栽58 d后逐漸高于T0.5處理。2 a T0.9和T0.5的T分別在3.76~13.68和3.89~11.98 mmol/(m2·s)之間變化,g分別在0.25~0.89 mol/(m2·s)和0.17~0.92 mol/(m2·s)之間變化。2015和2016年全生育期內(nèi)T0.9的平均T分別比T0.5處理高10.8%和14.7%,T0.9的平均g分別比T0.5處理高21.8%和16.3%(<0.05)。

    圖3 2015和2016年不同水分處理日光溫室番茄葉片蒸騰速率和氣孔導(dǎo)度生育期變化

    2.3 不同水分處理番茄植株莖稈液流變化規(guī)律

    2.3.1 生育期動態(tài)變化

    圖4給出了2015和2016年日莖流速率的逐日變化過程,監(jiān)測時間分別從移栽后69 d和41 d開始??梢钥闯?,2015和2016年T0.9與T0.5處理的日莖流速率差異分別在移栽后第69天和58天后開始出現(xiàn)差異,且隨著生育期的推進,二者的差異有增大趨勢,其中晴天差異最大,而陰雨天和灌水后的差異最小。2 a T0.9的日莖流速率分別較T0.5處理高54.8%和41.2%(<0.05)。不同水分處理之間,日光溫室番茄的日最大莖流速率差異較大,2015年T0.9和T0.5處理分別為6.74和4.61 mm/d,2016年T0.9和T0.5處理分別為6.98和4.51 mm/d,T0.9的日最大莖流速率分別較T0.5處理高31.6%和35.4%(<0.05)。

    圖4 2 a不同處理溫室番茄日莖流速率生育期動態(tài)變化

    2.3.2 不同典型天氣日變化規(guī)律

    選擇水分處理出現(xiàn)差異后的晴天(2015年6月11日和2016年6月16日)、多云天(2015年6月20日和2016年6月21日)和陰雨天(2015年6月24日和2016年6月23日)3種典型天氣狀況進行對比分析,2 a 3種典型天氣不同水分處理單株番茄莖稈液流和太陽輻射(R)的日變化曲線相似。以2016年為例(圖5),不同水分處理的莖流速率呈單峰曲線變化,莖流從8:00左右開始啟動并隨著R的增大而迅速增大,最大R和莖流速率均出現(xiàn)在12:00—13:00之間,之后莖流速率隨R減弱而降低,20:00之后以及整個晚上的莖流速率都處于較低水平。莖流速率變化相對R具有滯后性,平均滯后R約1 h,這是因為番茄植株的生理活動需要適應(yīng)R的變化[4]。水分脅迫抑制了番茄植株的蒸騰作用,在晴天尤為明顯,2016年晴天T0.9的日最大莖流速率分別比T0.5高85.3%(< 0.05),而陰雨天T0.9和T0.5的日莖流速率無明顯差異。

    圖5 2016年不同處理日光溫室番茄植株莖流速率典型天氣日變化

    2.4 不同水分處理番茄群體蒸散量變化規(guī)律

    2.4.1 日變化規(guī)律

    稱重式蒸滲儀可精確測量群體尺度蒸散量(evapotranspiration,ET)的日變化過程,掌握作物1 d內(nèi)的耗水規(guī)律。選擇番茄耗水強度最大的盛果期,選擇6月20—23日蒸滲儀的日變化數(shù)據(jù)進行分析(圖6)??梢钥闯觯琓0.9和T0.5處理的ET最大值均出現(xiàn)在13:00附近,T0.9處理的峰值明顯高于T0.5處理,最大差值為0.42 mm/h。群體尺度ET同樣隨R變化,但T0.9和T0.5處理的ET隨R的變化波動幅度不同,如6月20日R從0增加到 658.2 W/m2時,T0.9的ET從0.06 mm增加到0.69 mm,而T0.5僅從0.01 mm增加到0.31 mm,T0.9的增長速率較T0.5高52.4%。陰雨天(6月23日),由于R較小,T0.9與T0.5的ET之間無明顯差異。

    圖6 2016年不同處理日光溫室番茄群體尺度蒸散量日變化

    2.4.2 生育期動態(tài)變化

    全生育期內(nèi),番茄日ET呈先增大后減小的變化趨勢,ET在番茄盛果期達到最大(圖7),即移栽后第96~115天,T0.9和T0.5分別占全生育期總ET的31.7%和34.7%,日平均蒸散量分別為4.45和2.21 mm/d。全生育期T0.9和T0.5的日蒸散量分別在0.32~6.65和0.15~ 5.91 mm/d之間變化,生育期T0.9和T0.5的總ET分別為348.96和268.69 mm,T0.9較T0.5高23%。苗期、開花坐果期、盛果期和采摘期T0.9的ET分別較T0.5高0.9%、1.2%、15.6%和53.8%。可以看出,定植60 d后,T0.9與T0.5處理差異逐漸增大,尤其到了盛果期之后,其中采摘期的差異達到最大,最大可達3.47 mm/d??梢姡谕饨绛h(huán)境相同條件下(輻射水汽壓差和溫度等),根系層土壤含水率狀況是影響作物耗水過程的重要因素。

    圖7 2016年不同水分處理日光溫室番茄群體蒸散量的生育期動態(tài)變化

    2.5 影響日光溫室不同空間尺度蒸散量的主控環(huán)境因子分析

    選擇凈輻射(R)、水汽壓差()、空氣溫度(T)、葉片溫度(T)和室內(nèi)風(fēng)速()作為主控環(huán)境因子進行通徑分析。表3為T0.9和T0.5處理T、S和ET與主控環(huán)境因子的通徑分析結(jié)果,通過篩選具有顯著性指標的變量因子構(gòu)成多元線性回歸方程參數(shù),并對其進行顯著性分析可知,不同水分條件下各線性方程的變量因子不同,其中RTS和ET尺度上的影響均達到顯著性水平。從不同尺度分析結(jié)果來看,對于葉片尺度T,僅有1個顯著因子R,T0.9和T0.5處理的直接通徑系數(shù)分別為0.88和0.82。對單株尺度S,變量R、T和均達到顯著水平,2種水分處理的決定系數(shù)和直接通徑系數(shù)均為R最大,變量T通過其他環(huán)境變量的間接作用高于其他環(huán)境變量,所不同的是在T0.5水分處理條件下,S的估算需考慮T的直接作用。對于群體尺度ET,T0.9和T0.5的決定系數(shù)同樣是R最大,分別為0.95和0.91。對T0.9處理,變量T通過R和的間接作用均高于其他環(huán)境變量,而T0.5處理的間接作用則表現(xiàn)為T最高??梢?,R對葉片、植株和群體3個尺度的耗水作用均顯著高于其他環(huán)境變量,說明R是影響蒸散量的主要環(huán)境因子。

    表3 不同水分條件下葉片蒸騰速率、單株莖流速率和群體蒸散量與主控環(huán)境因子的通徑分析結(jié)果

    注:T為葉片蒸騰速率;S為莖流速率;ET為蒸散量;R為凈輻射;T為空氣溫度;為風(fēng)速;T為葉片溫度。

    Note:Tis leaf transpiration rate,Sis sap flow rate, FET is evapotranspiration,Ris net radiation,Tis air temperature,is wind speed,Tis leaf temperature.

    3 討 論

    3.1 日光溫室番茄不同作物空間尺度蒸散量變化及其影響因素分析

    水是植物進行一切生理活動并合成干物質(zhì)的主要成分,其中95%~99%的水分用于蒸騰作用。植物的蒸騰作用受自身生理特性和周圍環(huán)境因子的共同影響,其中環(huán)境因子包括太陽輻射(R)、溫度(T)、相對濕度(RH)、風(fēng)速()和根系層土壤水分狀況等[4]。在外界氣象條件一致的情況下,土壤水分含量是影響葉片蒸騰(T)、氣孔導(dǎo)度(g)、單株莖流速率(S)和群體蒸散量的主要因子,土壤水分通過影響植物氣孔開度間接影響植物體內(nèi)水分通過氣孔進入大氣的過程。當土壤水分不受限制時,根區(qū)土壤含水率對作物的耗水過程影響不明顯,當植物處于水分脅迫狀態(tài)時,作物的g和蒸騰作用都會隨根區(qū)土壤水分的減小而下降[25]。這是因為當作物經(jīng)受水分脅迫時,會不斷在根區(qū)產(chǎn)生脫落酸并經(jīng)由木質(zhì)部傳輸?shù)饺~片減小氣孔開度,進而抑制葉片的蒸騰作用[26]。

    單葉片尺度上,Tg密切相關(guān),共同受土壤水分的影響,本研究認為T0.9的Tg分別較T0.5高(圖2)。研究表明:葉齡越大,葉片內(nèi)的葉綠素含量就越低,葉片中的超氧化物歧化酶和過氧化氫酶活性降低,從而導(dǎo)致葉片的蒸騰速率降低[27]。本研究中,番茄打頂后,所測功能葉(葉片生長點下第3片葉)相對固定,并隨著生育期的推進而逐漸衰老,最終導(dǎo)致Tg的下降(圖3)。R是影響Tg的關(guān)鍵因子,并隨R的增大而逐漸增大,Tg的最大值出現(xiàn)在10:00—14:00之間,之后隨R的減小而降低,這與Zhang[4]等的研究結(jié)論一致。另外在晴朗無云的下午時段,隨著和T的不斷升高,植物為避免過多水分流失而減小氣孔開度,從而引起Tg的迅速下降[4]。

    單株尺度方面,通常用莖流速率來表征作物的蒸騰速率,莖流速率同樣隨R的變化而波動,R和莖流速率的峰值均出現(xiàn)在12:00—13:00之間,但莖流速率滯后R約1 h(圖4),這是由于番茄進行蒸騰作用時所需水分從根系經(jīng)由莖稈木質(zhì)部最終到達葉片需要一定時間,該過程需適應(yīng)R的變化[4]。本研究認為,當R超過650 W/m時(圖4),由于作物蒸騰作用強烈,根系吸水速率不能完全適應(yīng)作物蒸騰耗水的需求,會引起番茄植株莖流速率的間歇性變化,這也是導(dǎo)致中午12:00—13:00期間莖流速率變化不穩(wěn)定的主要因素,本結(jié)論與Zhang等[4,7,28]一致。根區(qū)土壤水分含量在一定程度上影響了單株番茄的耗水過程,植株長期處于水分脅迫條件下,加速了葉片的衰老過程,葉細胞活性不斷降低,其生理機制逐漸弱化,導(dǎo)致植株蒸騰能力降低[29],劉浩[7]對日光溫室番茄和Yamane等[30]對桃樹的研究結(jié)論同樣表明水分虧缺狀況下的莖流量顯著降低,本研究中T0.9的日最大莖流速率分別較T0.5處理高31.6%和35.4%,T0.9的日平均莖流速率較T0.5高41.2%~54.8%(圖4)。

    群體尺度方面,ET的日變化過程與S相似,同樣受R的影響,并在12:00—13:00期間有輕微的下凹趨勢(圖6),其原因與S相同。考慮ET與植株蒸騰、氣孔導(dǎo)度、光合作用等生理過程密切相關(guān),因此,可在R超過650 W/m時采用遮陽等措施降低室內(nèi)輻射量,以減少作物的奢侈耗水[31]。另外,20:00之后,ET逐漸降低,但并不為0,這部分能量是白天儲存于土壤和保溫墻內(nèi)的熱能,由于夜晚這部分能量較小,導(dǎo)致夜間ET的變化相對穩(wěn)定。日光溫室番茄的ET在盛果期達到最大(圖7),盛果期ET占全生育期總ET的31.7 %~34.7 %。移栽58 d后T0.9與T0.5的ET開始出現(xiàn)差異,并隨著生育期推進逐步增大,于采摘期達到最大,為53.8%,這與Liu等[20]等研究結(jié)論一致。牛勇等[10]利用大型稱重式蒸滲儀研究了日光溫室黃瓜ET的逐日變化規(guī)律,認為滴灌條件下黃瓜的ET在盛果期最大,占總ET的69.6%;Wang等[32]對日光溫室地下滴灌黃瓜ET的研究同樣表明,盛果期ET最大,在60.9~100.42 mm之間變化。類似研究結(jié)果在溫室甜瓜[33]和辣椒[14]等作物上均得到驗證。這是因為在結(jié)果盛期,輻射強度大以及高溫低濕的環(huán)境特征導(dǎo)致葉片蒸騰速率增加,加之果實膨大期對水分的需求比較旺盛,作物根系需源源不斷從土壤中吸收水分以滿足自身的蒸騰需求[10,14,32-33]。

    3.2 日光溫室不同空間尺度蒸散量的主控因子分析

    影響不同空間尺度蒸散量的環(huán)境因子包括R、T和等,其中R是多尺度蒸散量的主控因子。具體而言,R既能誘導(dǎo)氣孔開啟與關(guān)閉,又決定空氣溫度與相對濕度的變化,太陽光照可以提高葉面溫度,使葉內(nèi)外的水汽壓差增大,從而加強蒸騰速率。類似研究在溫室辣椒[14]和香蕉[15]]等作物上得到了驗證。第2個主控環(huán)境因子為,在單株尺度和群體尺度上的影響均達到顯著水平,第3和第4個環(huán)境主控因子依次為T和,而T僅對水分脅迫條件下的S和ET有效。Zhang等[4]認為對于充分供水條件的玉米,R的影響作用最大,第2個和第3個影響因子分別為和T,而的影響作用最小;徐璐等[13]認為在干旱荒漠區(qū)影響梭梭和檉柳不同尺度蒸騰作用的影響因子T大于,這可能是由于不同植物間自身特性差異以及對外界環(huán)境因子的響應(yīng)程度不同造成的;而對于熱帶雨林,影響莖桿速率的主要因子為,R次之[34],這可能是由于環(huán)境特征不同造成的。通過以上分析,可得到2個結(jié)論:1)日光溫室環(huán)境下,估算單株尺度和群體尺度ET時,是不可忽略的主要環(huán)境因子之一;2)水分脅迫條件下,考慮葉溫變量可顯著提高S和ET的估算精度,而對于充分供水作物則無需考慮。張大龍等[35]認為氣象環(huán)境因子對甜瓜蒸騰量的影響力很大程度上取決于土壤水分含量,當供水虧缺或過飽和時,甜瓜蒸騰與環(huán)境因子的相關(guān)性較差。由于土壤性質(zhì)和含水率本身存在著較大的空間變異性,以及植物蒸騰受多方面環(huán)境因子的影響,這對解釋水分脅迫條件下多尺度蒸散量主控因子造成一定難度,因此有關(guān)水分脅迫條件下,多尺度蒸散量主控因子的研究報道較少。而充分供水條件下,Rana等[36]通過分析環(huán)境參數(shù)對ET的敏感程度發(fā)現(xiàn),變量R、和對ET均具有較高的敏感性??梢姡ㄟ^考慮水分和環(huán)境因素的空間變異性來確定多尺度蒸散量的主控因子是進一步構(gòu)建尺度提升模型以及多尺度間轉(zhuǎn)換方法的關(guān)鍵。

    本研究未提出從單葉到植被冠層到群體尺度ET的蒸散量提升模型以及轉(zhuǎn)換方法,主要限制條件如下:1)從單葉尺度外推到單株尺度時需考慮植株不同層次處葉片蒸騰速率的綜合響應(yīng)機制,這就需要測定不同層次處不同遮蔭狀況葉片的蒸騰速率,并考慮冠層內(nèi)部微環(huán)境變化對葉片氣孔開度的影響,從而構(gòu)建基于單葉片的單株蒸騰綜合模型,實現(xiàn)單葉到植被冠層的提升。而采用葉片蒸騰和同期環(huán)境參數(shù)上推到植被冠層尺度的經(jīng)驗性模型(通常為多元線性回歸方程),可能具有一定的局限性;2)從單株外推到群體尺度ET時需考慮作物生育期、土壤水分和環(huán)境條件的空間變異性,本文雖考慮了土壤水分的空間變異性,但缺少對多尺度環(huán)境變量的監(jiān)測資料,而以往外推模型同樣是依靠經(jīng)驗性回歸方法估算群體蒸散量。因此,建議針對以上主要限制因素進行系統(tǒng)試驗研究并驗證,最終實現(xiàn)從經(jīng)驗性模型向機理性模型的轉(zhuǎn)變。

    4 結(jié) 論

    本文以日光溫室番茄為研究對象,通過2 a田間試驗測定了不同尺度蒸散量的變化過程,分析了不同水分條件下多尺度蒸散量的變化規(guī)律以及主控環(huán)境因子。結(jié)論如下:

    1)不同水分處理葉片蒸騰速率和氣孔導(dǎo)度在移栽54~58 d后出現(xiàn)差異,2 a充分供水(T0.9)的生育期內(nèi)葉片蒸騰速率和氣孔導(dǎo)度分別比虧缺供水(T0.5)處理高10.8 %~14.7 %和16.3 %~21.8 % (<0.05),全生育期葉片蒸騰速率和氣孔導(dǎo)度的變化規(guī)律表現(xiàn)為苗期小,盛果期大。

    2)T0.9和T0.5處理的莖流速率在晴天差異最大,陰雨天最小,莖流速率隨太陽輻射的波動而變化,但滯后太陽輻射約1 h。T0.9的日莖流速率較T0.5高41.2%~54.8% (<0.05)。

    3)群體ET表現(xiàn)為先增大后減小的變化趨勢,在盛果期達到最大,全生育期T0.9和T0.5處理的日ET分別在0.32~6.65 mm/d和0.15~5.91 mm/d之間變化。

    4)日光溫室環(huán)境下,影響不同尺度蒸散量的環(huán)境因子依次為凈輻射>水汽壓差>氣溫>風(fēng)速>葉片溫度。凈輻射在葉片、單株以及群體尺度上的影響均達到顯著水平,估算單株尺度蒸騰量和群體尺度ET時,風(fēng)速是不可忽略的環(huán)境因子之一;水分脅迫條件下,考慮變量葉溫可顯著提高單株蒸騰和群體ET的估算精度。

    [1] 崔遠來,董斌,李遠華,等. 農(nóng)業(yè)灌溉節(jié)水評價指標與尺度問題[J]. 農(nóng)業(yè)工程學(xué)報,2007,23(7):1-7. Cui Yuanlai, Dong Bin, Li Yuanhua, et al. Assessment indicators and scales of water saving in agricultural irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(7):1-7. (in Chinese with English abstract)

    [2] 許迪. 灌溉水文學(xué)尺度轉(zhuǎn)換問題研究綜述[J]. 水利學(xué)報,2006,37(2):141-149. Xu Di. Review on scaling study in irrigation hydrology [J]. Journal of Hydraulic Engineering, 2006, 37(2): 141-149. (in Chinese with English abstract)

    [3] 張寶忠,許迪,劉鈺,等. 多尺度蒸散發(fā)估測與時空尺度拓展方法研究進展[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(6):8-16. Zhang Baozhong, Xu Di, Liu Yu, et al. Review of multi-scale evapotranspiration estimation and spatio-temporal scale expansion[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(6): 8-16. (in Chinese with English abstract)

    [4] Zhang Baozhong, Xu Di, Liu Yu, et al. Multi-scale evapotranspiration of summer maize and the controlling meteorological factors in north China [J]. Agricultural and Forest Meteorology, 2016, 216(5): 1-12.

    [5] 李銀坤,武雪萍,吳會軍,等. 水氮條件對溫室黃瓜光合日變化及產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報,2010,26(增刊1): 122-129. Li Yinkun, Wu Xueping, Wu Huijun, et al. Effects of water and nitrogen conditions on the diurnal variation of photosynthesis and yield of cucumber in greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(Supp.1): 122-129. (in Chinese with English abstract)

    [6] Tartachnyk I I, Blanke M M. Photosynthesis and transpiration of tomato and CO2fluxes in a greenhouse under changing environmental conditions in winter[J]. Annals of Applied Biology, 2007, 150(2): 149-156.

    [7] 劉浩,孫景生,段愛旺,等. 溫室滴灌條件下番茄植株莖流變化規(guī)律試驗[J]. 農(nóng)業(yè)工程學(xué)報,2010,26(10):77-82. Liu Hao, Sun Jingsheng, Duan Aiwang, et al. Experiments on variation of tomato sap flow under drip irrigation conditions in greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010,26(10): 77-82. (in Chinese with English abstract)

    [8] Qiu Rangjian, Kang Shaozhong, Du Taisheng, et al. Effect of convection on the Penman-Monteith model estimates of transpiration of hot pepper grown in solar greenhouse[J]. Scientia Horticulturae, 2013, 160(3): 163-171.

    [9] 劉海軍,黃冠華,Josef Tanny, 等. 溫室內(nèi)香蕉樹蒸騰量預(yù)報的試驗研究[J]. 農(nóng)業(yè)工程學(xué)報,2007,23(3):174-179. Liu Haijun, Huang Guanhua, Josef Tanny, et al. Estimation of evapotranspiration of banana plants in a greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(3): 174-179. (in Chinese with English abstract)

    [10] 牛勇,劉洪祿,吳文勇,等. 基于大型稱重式蒸滲儀的日光溫室黃瓜蒸騰規(guī)律研究[J]. 農(nóng)業(yè)工程學(xué)報,2011,27(1):52-56. Niu Yong, Liu Honglu, Wu Wenyong, et al. Cucumber transpiration by large–scale weighting lysimeter in solar greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(1): 52-56. (in Chinese with English abstract)

    [11] 蔡甲冰,許迪,劉鈺,等. 冬小麥返青后騰發(fā)量時空尺度效應(yīng)的通徑分析[J]. 農(nóng)業(yè)工程學(xué)報,2011,27(8):69-76.Cai Jiabing, Xu Di, Liu Yu, et al. Path analysis on spatio-temporal scaling effect of crop evapotranspiration in growing seasons for winter wheat after reviving[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(8): 69-76. (in Chinese with English abstract)

    [12] 蔡甲冰,許迪,劉鈺,等. 冬小麥返青后作物騰發(fā)量的尺度效應(yīng)及其轉(zhuǎn)換研究[J]. 水利學(xué)報,2010,41(7):862-869. Cai Jiabing, Xu Di, Liu Yu, et al. Scaling effects and transformation of crop evapotranspiration for winter wheat after reviving[J]. Journal of Hydraulic Engineering, 2010, 41(7): 862-869. (in Chinese with English abstract)

    [13] 徐璐,李平衡,趙鑫,等. 2種典型荒漠植物多尺度蒸騰日變化特征及其尺度轉(zhuǎn)換研究[J]. 石河子大學(xué)學(xué)報,2014,32(1):79-85. Xu Lu, Li Pingheng, Zhao Xin, et al. Study on diurnal change characteristics of transpiration in different scales and up- scaling of two typical desert plants[J]. Journal of Shihezi University, 2014, 32(1): 79-85. (in Chinese with English abstract)

    [14] Qiu Rangjian, Kang Shaozhong, Li Fusheng, et al. Energy partitioning and evapotranspiration of hot pepper grown in greenhouse with furrow and drip irrigation methods[J]. Scientia Horticulturae, 2011, 129(4): 790-797.

    [15] Demrati H, Boulard T, Fatnassi H, et al. Microclimate and transpiration of a greenhouse banana crop [J]. Biosystems Engineering, 2007, 98(1): 66-78.

    [16] Wilson K B, Hanson P J, Mulholland P J, et al. A comparison of methods for determining forest evapotranspiration and its components: Sap-flow, soil water budget, eddy covariance and catchment water balance[J]. Agricultural and Forest Meteorology, 2001, 106(2): 153-168.

    [17] 陳倩倩,范陽陽,郝影賓,等. 不同土壤水分含量對玉米氣孔發(fā)育過程和蒸騰耗水量的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究,2011,29(3):75-79. Chen Qianqian, Fan Yangyang, Hao Yingbin, et al. Effects of different soil water content on stomata development and water consumption of maize[J]. Agricultural Research in the Arid Areas, 2011, 29(3): 75-79. (in Chinese with English abstract)

    [18] 姚勇哲,李建明,張榮,等. 溫室番茄蒸騰量與其影響因子的相關(guān)分析及模型模擬[J]. 應(yīng)用生態(tài)學(xué)報,2012,23(7):1869-1874. Yao Yongzhe, Li Jianming, Zhang Rong, et al. Greenhouse tomato transpiration and its affecting factors: Correlation analysis and model simulation[J]. Chinese Journal of Applied Ecology, 2012, 23(7): 1869-1874. (in Chinese with English abstract)

    [19] 龔雪文,孫景生,劉浩,等. 基于20 cm蒸發(fā)皿蒸發(fā)量制定的華北地區(qū)溫室黃瓜滴灌灌水制度[J]. 應(yīng)用生態(tài)學(xué)報,2015,26(11):3381-3388. Gong Xuewen, Sun Jingsheng, Liu Hao, et al. Irrigation scheduling with a 20 cm standard pan for drip-irrigated cucumber grown in solar greenhouse in the North China Plain[J]. Chinese Journal of Applied Ecology, 2015, 26(11): 3381-3388. (in Chinese with English abstract)

    [20] Liu Hao, Duan Aiwang, Li Fusheng, et al. Drip irrigation scheduling for tomato grown in solar greenhouse based on pan evaporation in North China plain[J]. Journal of Integrative Agriculture, 2013, 12 (3): 520-531.

    [21] 劉浩. 溫室番茄需水規(guī)律與優(yōu)質(zhì)高效灌溉指標研究[D]. 北京:中國農(nóng)業(yè)科學(xué)院,2010. Liu Hao. Water Requirement and Optimal Irrigation Index for Effective Water Use and High-Quality of Tomato in Greenhouse[D]. Beijing: Chinese Academy of Agricultural Sciences, 2010. (in Chinese with English abstract)

    [22] 陳新明,蔡煥杰,李紅星,等. 溫室內(nèi)作物騰發(fā)量計算與驗證[J]. 水科學(xué)進展,2007,18(6):812-815. Chen Xinming, Cai Huanjie, Li Hongxing, et al. Calculation and verification of crop evapotranspiration in greenhouse[J]. Advances in Water Science, 2007, 18(6): 812-815. (in Chinese with English abstract)

    [23] Seker H, Serin Y. Explanation of the relationships between seed yield and some morphological traits in smooth bromegrass() by path analysis[J]. European Journal of Agronomy, 2004, 21(1): 1-6.

    [24] 李波,任樹梅,楊培嶺,等. 供水條件對溫室番茄根系分布及產(chǎn)量影響[J]. 農(nóng)業(yè)工程學(xué)報,2007,23(9):39-44. Li Bo, Ren Shumei, Yang Peiling, et al. Impacts of water supply on root distribution and yield of tomato in greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(9): 39-44. (in Chinese with English abstract)

    [25] 陳仁升,康爾泗,趙文智,等. 中國西北干旱區(qū)樹木蒸 騰對氣象因子的響應(yīng)[J]. 生態(tài)學(xué)報,2004,24(3):477-485. Chen Rensheng, Kang Ersi, Zhao Wenzhi, et al. Trees transpiration response to meteorological variables in arid regions of Northwest China[J]. Acta Ecologica Sinica, 2004, 24(3): 477-485. (in Chinese with English abstract)

    [26] 胡笑濤,康紹忠,張建華,等.番茄垂向分根區(qū)交替控制滴灌室內(nèi)試驗及節(jié)水機理[J]. 農(nóng)業(yè)工程學(xué)報,2005,21(7):1-5. Hu Xiaotao, Kang Shaozhong, Zhang Jianhua, et al. Water-saving mechanism and efficiency of vertical partial- root zone alternative controlled drip irrigation of tomato ()[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2005, 21(7): 1-5. (in Chinese with English abstract)

    [27] 許守民,石連旋,胡波,等. 大豆葉片衰老過程中光合作用與有關(guān)酶活性變化的研究[J]. 東北師大學(xué)報,1997,8(4):75-78. Xu Shoumin, Shi Lianxuan, Hu Bo, et al. Changes in photosynthetic capacity and related enzyme activities during soybean() leaf senescence[J]. Journal of Northeast Normal University, 1997, 8(4): 75-78. (in Chinese with English abstract)

    [28] 彭致功,楊培嶺,段愛旺,等. 日光溫室條件下番茄植株蒸騰規(guī)律研究[J]. 干旱地區(qū)農(nóng)業(yè)研究,2004,22(1):62-65.Peng Zhigong, Yang Peiling, Duan Aiwang, et al. Experiment on transpiration of tomato in solar-heated greenhouse [J]. Agricultural Research in the Arid Areas, 2004, 22(1): 62-65. (in Chinese with English abstract)

    [29] Tardieu F. Plant tolerance to water deficit: Physical limits and possibilities for progress[J]. Comptes Rendus Geoscience, 2005, 337(1): 57-67.

    [30] Yamane T, Hamana Y, Nakano M. Detection of water-deficit stress from daily sap flow profiles in peach [J]. Journal of the Japanese Society for Horticultural Science, 2011, 80(4): 383-389.

    [31] 丁小濤,金海軍,張紅梅,等. 遮蔭處理對溫室四種蔬菜生長及光合作用日變化的影響[J]. 浙江農(nóng)業(yè)學(xué)報,2010,22(1):51-56. Ding Xiaotao, Jin Haijun, Zhang Hongmei, et al. Effect of shading on growth and diurnal photosynthetic changes of four vegetables in glasshouse[J]. Acta Agriculturae Zhejiangensis, 2010, 22(1): 51-56. (in Chinese with English abstract)

    [32] Wang Zhenying, Liu Zuoxin, Zhang Zikun, et al. Subsurface drip irrigation scheduling for cucumber (L.) grown in solar greenhouse based on 20 cm standard pan evaporation in Northeast China[J]. Scientia Horticulturae, 2009, 123(1):51-57.

    [33] Sensoy S, Ertek A, Gedik I, et al. Irrigation frequency and amount affect yield and quality of field-grown melon (L.)[J]. Agricultural Water Management, 2007, 88(1/2/3): 269-274.

    [34] Granier A, Huc R, Barigah S T. Transpiration of natural rain forest and its dependence on climatic factors[J]. Agricultural and Forest Meteorology, 1996, 78(78): 19-29.

    [35] 張大龍,常毅博,李建明,等. 大棚甜瓜蒸騰規(guī)律及其影響因子[J]. 生態(tài)學(xué)報,2014,34(4):953-962. Zhang Dalong, Chang Yibo, Li Jianming, et al. The critical factors of transpiration on muskmelon in plastic greenhouse[J]. Acta Ecologica Sinica, 2014, 34(4): 953-962. (in Chinese with English abstract)

    [36] Rana G, Katerji N. A Measurement based sensitivity analysis of the Penman-Monteith actual evapotranspiration model for crops of different height and in contrasting water status[J]. Theoretical & Applied Climatology, 1998, 60(1): 141-149.

    Variation of evapotranspiration in different spatial scales for solar greenhouse tomato and its controlling meteorological factors

    Gong Xuewen1,2, Liu Hao1, Sun Jingsheng1※, Zhang Hao1,2, Ma Xiaojian1,2, Wang Wanning1,2

    (1453003,2.100081,)

    Knowledge of multi-scale evapotranspiration (ET) is the basis of increasing agricultural water use efficiency and realizing optimal water resources allocation. However, there is little information about the characteristics of multi-scale ET and the difference in the quantitative response of leaf and plant transpiration and farm-scale ET to the meteorological factors, especially for the solar greenhouse, which limits our understanding of agricultural water and energy cycles. In this study, the ET at different scales was investigated in solar greenhouse tomato and their main influential factors were analyzed. Experiments were carried out in a typical solar greenhouse to investigate variation in leaf and plant transpiration and farmland ET of drip-irrigated tomato under different water conditions in 2 continuous seasons. Field trials were conducted during the growing season from March to June, 2015 and 2016 at Xinxiang Comprehensive Experimental Station, Chinese Academy of Agricultural Sciences (35°9′N, 113°47′E, and an altitude of 78.7 m), located in Xinxiang, Henan province. Tomato was used in the experiment on a plot of 8.8 m2(8.0 m long by 1.1 m wide) with 2 rows and 50 plants. There are 4 replications for each treatment. Tomato seedlings were transplanted on March 8th, 2015 and March 9th, 2016 respectively and irrigated with a drip irrigation system. Irrigation amount was determined based on the accumulated evaporation in 20 cm pan, and 2 treatments were designed with full irrigation and deficit irrigation. Leaf transpiration was measured with a photosynthesis system, plant transpiration (expressed as sap flow rate) measured by a sap flow meter and farm-scale ET measured by a weighing lysimeter. Besides, solar radiation, net radiation, relative humidity and air temperature inside the solar greenhouse were constantly monitored with an automatic weather station at 2.0 m height above ground level. And then the difference of quantitative response of leaf transpiration, plant transpiration and farm-scale ET to the meteorological factors was performed using the path analysis method. Results showed that leaf transpiration rate and stomatal conductance gradually increased with the increase in solar radiation, and reached their maximum values at 10:00–14:00, and then rapidly decreased with solar radiation. Difference in leaf transpiration rate between full irrigation and deficit irrigation appeared at 54-58 days after transplanting. Daily average of leaf transpiration rate in the full irrigation was 10.8% and 14.7% higher than that in deficit irrigation, respectively. Tomato sap flow rate changed with weather conditions, and the difference in plant transpiration between the full irrigation and deficit irrigation was the largest in the sunny day, and the smallest in rainy days, and the plant transpiration was related to solar radiation but lagged behind 1 hour. Daily average of leaf transpiration rate in the full irrigation was 54.8% and 41.2% higher than that in deficit irrigation, respectively in the 2 years. The daily farm-scale ET in the full irrigation and deficit irrigation ranged between 0.32-6.65 mm/d and 0.15-5.91 mm/d, respectively, over the whole growth stage, and the farm-scale evapotranspiration reached the maximum in the full fruit period, accounting for 31.7%-34.7% of the total ET. The path analysis of evapotranspiration at different scales and its controlling meteorological factors indicated that leaf transpiration rate, plant transpiration and farm-scale evapotranspiration could be characterized by the net solar radiation, and the vapor pressure deficit also had significant influence on the scale of individual plant and farmland scale. The influence of wind speed should be considered in estimating plant transpiration rate and farm-scale evapotranspiration. Besides, adding leaf temperature could significantly improve the estimation accuracy of plant transpiration and farm-scale evapotranspiration under the condition of deficit irrigation.This paper provides valuable information for constructing conversion methods and theoretic model at various spatio-temporal scales.

    evapotranspiration; transpiration; path analysis; stomatal conductance; sap flow rate

    10.11975/j.issn.1002-6819.2017.08.023

    S161.4;S626.5;S641.2

    A

    1002-6819(2017)-08-0166-10

    2016-08-01

    2016-12-10

    中央級科研院所基本科研業(yè)務(wù)費專項(中國農(nóng)業(yè)科學(xué)院農(nóng)田灌溉研究所)資助項目

    龔雪文,男,河南安陽人,博士生,主要從事作物水分生理與高效利用方面的研究。新鄉(xiāng) 中國農(nóng)業(yè)科學(xué)院農(nóng)田灌溉研究所,453003。 Email:gxw068@126.com

    孫景生,男,遼寧建平人,研究員,博士生導(dǎo)師,主要從事節(jié)水灌溉和作物高效用水技術(shù)方面的研究。新鄉(xiāng) 中國農(nóng)業(yè)科學(xué)院農(nóng)田灌溉研究所,453003。Email:jshsun623@163.com

    龔雪文,劉浩,孫景生,張昊,馬筱建,王萬寧.日光溫室番茄不同空間尺度蒸散量變化及主控因子分析[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(8):166-175. doi:10.11975/j.issn.1002-6819.2017.08.023 http://www.tcsae.org

    Gong Xuewen, Liu Hao, Sun Jingsheng, Zhang Hao, Ma Xiaojian, Wang Wanning. Variation of evapotranspiration in different spatial scales for solar greenhouse tomato and its controlling meteorological factors[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(8): 166-175. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.08.023 http://www.tcsae.org

    成人高潮视频无遮挡免费网站| 一区福利在线观看| 长腿黑丝高跟| 最后的刺客免费高清国语| 亚洲人与动物交配视频| 黄色欧美视频在线观看| 哪个播放器可以免费观看大片| av又黄又爽大尺度在线免费看 | 波多野结衣高清作品| 国国产精品蜜臀av免费| 一区二区三区免费毛片| 高清毛片免费看| 国产亚洲精品av在线| 日本黄大片高清| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人av| 亚洲五月天丁香| 三级毛片av免费| 国内精品久久久久精免费| 变态另类成人亚洲欧美熟女| 丰满人妻一区二区三区视频av| 超碰av人人做人人爽久久| 国产黄色视频一区二区在线观看 | 在线a可以看的网站| 亚洲国产精品sss在线观看| 精品久久久久久久久av| 国产精品嫩草影院av在线观看| 校园春色视频在线观看| 亚洲最大成人中文| av在线观看视频网站免费| 国产黄a三级三级三级人| 欧美xxxx黑人xx丫x性爽| 亚洲一区高清亚洲精品| kizo精华| 只有这里有精品99| 深爱激情五月婷婷| 久久久国产成人精品二区| 波多野结衣巨乳人妻| 午夜福利在线观看免费完整高清在 | 亚洲国产精品成人综合色| 美女高潮的动态| 成人欧美大片| 两个人的视频大全免费| 亚洲精品色激情综合| 国产精品久久久久久av不卡| 日本三级黄在线观看| 久久久久国产网址| 精华霜和精华液先用哪个| 我要看日韩黄色一级片| 一级黄片播放器| 国产精品电影一区二区三区| 在线观看av片永久免费下载| 成人特级av手机在线观看| 午夜老司机福利剧场| 亚洲av电影不卡..在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲一区二区三区色噜噜| 国产蜜桃级精品一区二区三区| 国产高清有码在线观看视频| 欧美极品一区二区三区四区| 久久久久久久久久久免费av| 国产成人91sexporn| 亚洲国产精品成人综合色| 中文字幕免费在线视频6| 精品少妇黑人巨大在线播放 | 欧美性猛交╳xxx乱大交人| 精品人妻熟女av久视频| 久久久欧美国产精品| 人妻夜夜爽99麻豆av| 国产男人的电影天堂91| 熟女人妻精品中文字幕| 日本一本二区三区精品| 悠悠久久av| 国产一区亚洲一区在线观看| 成人av在线播放网站| 日韩视频在线欧美| 1000部很黄的大片| 国产黄a三级三级三级人| 国产精品久久久久久精品电影小说 | 亚洲三级黄色毛片| 日韩一本色道免费dvd| 国产精品一及| videossex国产| 极品教师在线视频| 国产精品国产高清国产av| eeuss影院久久| 在线观看美女被高潮喷水网站| 国产日本99.免费观看| 免费av观看视频| 大型黄色视频在线免费观看| 久久久久久国产a免费观看| 日韩 亚洲 欧美在线| 联通29元200g的流量卡| 国产高潮美女av| 亚洲av.av天堂| 国产精品久久视频播放| 中文字幕制服av| 日韩欧美在线乱码| 黄片wwwwww| 91久久精品国产一区二区三区| 国产成人aa在线观看| 国产精品野战在线观看| 2022亚洲国产成人精品| 综合色丁香网| 精品久久久久久久久久久久久| 国产精品乱码一区二三区的特点| 婷婷色综合大香蕉| 日韩高清综合在线| 国产精品野战在线观看| 91av网一区二区| 少妇猛男粗大的猛烈进出视频 | 深夜精品福利| 最近的中文字幕免费完整| 欧美变态另类bdsm刘玥| 22中文网久久字幕| 免费大片18禁| 高清在线视频一区二区三区 | 久久精品91蜜桃| 深夜精品福利| 美女高潮的动态| 国内精品一区二区在线观看| av免费在线看不卡| 99在线视频只有这里精品首页| 欧美潮喷喷水| 亚洲成人精品中文字幕电影| 亚洲五月天丁香| 日韩欧美精品v在线| 噜噜噜噜噜久久久久久91| 亚洲一区高清亚洲精品| 人人妻人人澡人人爽人人夜夜 | 成人特级黄色片久久久久久久| 久久精品国产清高在天天线| 美女cb高潮喷水在线观看| 老女人水多毛片| 精品熟女少妇av免费看| 青春草亚洲视频在线观看| 久久精品久久久久久久性| 日韩国内少妇激情av| 婷婷色综合大香蕉| 好男人在线观看高清免费视频| 欧美潮喷喷水| 少妇被粗大猛烈的视频| 99在线视频只有这里精品首页| 又黄又爽又刺激的免费视频.| 欧美三级亚洲精品| 成年女人永久免费观看视频| 少妇熟女欧美另类| h日本视频在线播放| 菩萨蛮人人尽说江南好唐韦庄 | 国产中年淑女户外野战色| 日韩成人伦理影院| www.色视频.com| 久久国内精品自在自线图片| 欧美+日韩+精品| www日本黄色视频网| 亚洲五月天丁香| 大又大粗又爽又黄少妇毛片口| av.在线天堂| 麻豆国产97在线/欧美| 久久综合国产亚洲精品| 精品久久久久久久久久久久久| 寂寞人妻少妇视频99o| 九九久久精品国产亚洲av麻豆| 夜夜爽天天搞| 国产日韩欧美在线精品| 久久99精品国语久久久| 久久国内精品自在自线图片| 久久亚洲国产成人精品v| 哪个播放器可以免费观看大片| videossex国产| av专区在线播放| 日韩一本色道免费dvd| 午夜亚洲福利在线播放| 亚洲精品乱码久久久v下载方式| 你懂的网址亚洲精品在线观看 | 麻豆国产av国片精品| 日韩成人伦理影院| 日韩三级伦理在线观看| 午夜激情福利司机影院| 九九在线视频观看精品| 国产精品精品国产色婷婷| 精品99又大又爽又粗少妇毛片| 成人亚洲欧美一区二区av| 亚洲熟妇中文字幕五十中出| 午夜福利成人在线免费观看| 一区二区三区高清视频在线| 夜夜爽天天搞| 国产在线男女| 菩萨蛮人人尽说江南好唐韦庄 | 一级av片app| 网址你懂的国产日韩在线| 久久精品国产清高在天天线| 日本撒尿小便嘘嘘汇集6| 九草在线视频观看| 少妇的逼好多水| 91精品国产九色| 成人亚洲欧美一区二区av| 男女边吃奶边做爰视频| 99久久成人亚洲精品观看| 国产女主播在线喷水免费视频网站 | 欧美又色又爽又黄视频| 男人的好看免费观看在线视频| 成人欧美大片| 久久久久久伊人网av| 热99re8久久精品国产| 日本熟妇午夜| 国产精品人妻久久久影院| 成年免费大片在线观看| 日本黄大片高清| 青青草视频在线视频观看| 精品久久久久久久久久久久久| 日韩一本色道免费dvd| 天天躁日日操中文字幕| 国产探花在线观看一区二区| 成人三级黄色视频| 高清毛片免费观看视频网站| 91精品一卡2卡3卡4卡| 国产私拍福利视频在线观看| 99国产精品一区二区蜜桃av| 好男人在线观看高清免费视频| 成年女人永久免费观看视频| 亚洲欧美日韩高清专用| 一卡2卡三卡四卡精品乱码亚洲| 亚洲三级黄色毛片| 国产女主播在线喷水免费视频网站 | 欧美精品国产亚洲| 亚洲高清免费不卡视频| 别揉我奶头 嗯啊视频| 亚洲人成网站在线播放欧美日韩| av又黄又爽大尺度在线免费看 | 国产蜜桃级精品一区二区三区| 日韩成人av中文字幕在线观看| 99热6这里只有精品| 欧美潮喷喷水| 22中文网久久字幕| 国产精品日韩av在线免费观看| 国产麻豆成人av免费视频| 亚洲欧美成人精品一区二区| 黄片无遮挡物在线观看| 观看美女的网站| 欧美极品一区二区三区四区| 国产精品.久久久| 草草在线视频免费看| 免费搜索国产男女视频| 午夜精品国产一区二区电影 | 少妇人妻一区二区三区视频| 国产麻豆成人av免费视频| 在线观看66精品国产| 亚洲无线在线观看| 成人美女网站在线观看视频| 亚洲精品成人久久久久久| 久久久国产成人免费| 欧美激情久久久久久爽电影| 成人漫画全彩无遮挡| 久久精品国产鲁丝片午夜精品| 亚洲激情五月婷婷啪啪| 最近手机中文字幕大全| 三级男女做爰猛烈吃奶摸视频| 国产白丝娇喘喷水9色精品| 五月玫瑰六月丁香| 国产精品久久电影中文字幕| 欧美日本亚洲视频在线播放| 麻豆成人午夜福利视频| 欧美日韩综合久久久久久| 91久久精品国产一区二区三区| 午夜久久久久精精品| 此物有八面人人有两片| 在线观看午夜福利视频| 亚洲欧美日韩东京热| 国产成人精品婷婷| 成人亚洲精品av一区二区| 国产免费一级a男人的天堂| 一夜夜www| 啦啦啦观看免费观看视频高清| 男女做爰动态图高潮gif福利片| 久久久久网色| 黄色日韩在线| 一夜夜www| 国产伦精品一区二区三区四那| 久久精品国产亚洲av香蕉五月| 秋霞在线观看毛片| 国产精品一区二区性色av| av视频在线观看入口| 午夜免费男女啪啪视频观看| 国产男人的电影天堂91| 国产精品麻豆人妻色哟哟久久 | 一进一出抽搐gif免费好疼| 深爱激情五月婷婷| 欧美高清成人免费视频www| 亚州av有码| 亚洲av不卡在线观看| 国产乱人视频| 亚洲欧美中文字幕日韩二区| 天堂中文最新版在线下载 | 国产在线男女| 变态另类丝袜制服| 精品久久久噜噜| 少妇丰满av| 欧美丝袜亚洲另类| 22中文网久久字幕| 亚洲国产精品成人久久小说 | 热99re8久久精品国产| 成人特级黄色片久久久久久久| 久久6这里有精品| 欧美变态另类bdsm刘玥| 欧美性猛交黑人性爽| 中文在线观看免费www的网站| 国产精品无大码| av女优亚洲男人天堂| 人体艺术视频欧美日本| 精品久久久久久久末码| 国产精品.久久久| 我要搜黄色片| 亚洲成av人片在线播放无| 亚洲精品粉嫩美女一区| 联通29元200g的流量卡| 亚洲国产欧洲综合997久久,| 国产又黄又爽又无遮挡在线| 久久久成人免费电影| 啦啦啦啦在线视频资源| 日韩精品青青久久久久久| 给我免费播放毛片高清在线观看| 在线观看66精品国产| 91在线精品国自产拍蜜月| 亚洲精华国产精华液的使用体验 | 欧美日韩国产亚洲二区| 床上黄色一级片| 又黄又爽又刺激的免费视频.| 男人舔女人下体高潮全视频| 韩国av在线不卡| 国产高潮美女av| 禁无遮挡网站| 黄色日韩在线| 好男人视频免费观看在线| 国产精品一区二区三区四区免费观看| 国产高清有码在线观看视频| 国产成人freesex在线| 欧美性猛交黑人性爽| av专区在线播放| 国产成年人精品一区二区| 在线播放国产精品三级| 哪个播放器可以免费观看大片| 日韩一区二区三区影片| 国产极品天堂在线| 最近的中文字幕免费完整| 成人毛片60女人毛片免费| 能在线免费看毛片的网站| 国产成人freesex在线| 午夜亚洲福利在线播放| 不卡视频在线观看欧美| 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 国产成人福利小说| 美女被艹到高潮喷水动态| 丰满人妻一区二区三区视频av| 久久精品国产亚洲av天美| 亚洲欧美日韩无卡精品| 久久草成人影院| 村上凉子中文字幕在线| 99热这里只有是精品在线观看| 久久久午夜欧美精品| 一进一出抽搐动态| 一进一出抽搐gif免费好疼| 日本欧美国产在线视频| 成人毛片a级毛片在线播放| 最近手机中文字幕大全| 欧美日韩一区二区视频在线观看视频在线 | 啦啦啦啦在线视频资源| 国产av在哪里看| 国产成人aa在线观看| 悠悠久久av| 午夜爱爱视频在线播放| 亚洲精品乱码久久久久久按摩| 毛片女人毛片| 全区人妻精品视频| 六月丁香七月| 波野结衣二区三区在线| 亚洲高清免费不卡视频| 国产熟女欧美一区二区| 日本欧美国产在线视频| 一区二区三区高清视频在线| 日本欧美国产在线视频| 一进一出抽搐gif免费好疼| 久久久国产成人免费| 成人毛片60女人毛片免费| 国产精品嫩草影院av在线观看| 亚洲第一区二区三区不卡| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 久久国内精品自在自线图片| 亚洲成人久久性| 亚洲av免费高清在线观看| 天堂中文最新版在线下载 | 国产黄片美女视频| 亚洲精品久久国产高清桃花| 在线a可以看的网站| 中国美女看黄片| h日本视频在线播放| 欧美成人a在线观看| 黄色一级大片看看| 亚洲五月天丁香| 国产高清不卡午夜福利| 一个人观看的视频www高清免费观看| 伦理电影大哥的女人| 成人午夜高清在线视频| 色综合色国产| 日韩欧美精品免费久久| 国产熟女欧美一区二区| 久久久久网色| 亚洲精品自拍成人| 欧美3d第一页| 日韩精品有码人妻一区| 亚洲天堂国产精品一区在线| 久久久久久久久久黄片| 国产色婷婷99| 色哟哟·www| 麻豆一二三区av精品| 麻豆国产97在线/欧美| 波多野结衣巨乳人妻| 舔av片在线| 久久亚洲国产成人精品v| 亚洲av一区综合| 成人永久免费在线观看视频| 亚洲欧美日韩卡通动漫| 精品不卡国产一区二区三区| 最近手机中文字幕大全| 欧美精品一区二区大全| 可以在线观看的亚洲视频| ponron亚洲| 51国产日韩欧美| 欧美+亚洲+日韩+国产| 国产精品久久久久久久久免| 热99在线观看视频| 国产综合懂色| 又粗又硬又长又爽又黄的视频 | 国产一区二区在线观看日韩| 九九爱精品视频在线观看| 久久久久久伊人网av| 国内久久婷婷六月综合欲色啪| 久久久成人免费电影| av在线老鸭窝| 一边亲一边摸免费视频| 少妇的逼水好多| 久久久精品欧美日韩精品| 特级一级黄色大片| 成人午夜高清在线视频| 99视频精品全部免费 在线| 最近手机中文字幕大全| 午夜福利高清视频| 免费看av在线观看网站| av国产免费在线观看| 性欧美人与动物交配| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 免费电影在线观看免费观看| 日韩一区二区视频免费看| 国产一区二区亚洲精品在线观看| 亚洲五月天丁香| 女同久久另类99精品国产91| 亚洲国产精品成人久久小说 | 欧美另类亚洲清纯唯美| 国产精品一及| 天堂av国产一区二区熟女人妻| 一级毛片电影观看 | 亚洲色图av天堂| 国产精品久久视频播放| 悠悠久久av| 久久久久性生活片| 精华霜和精华液先用哪个| 久久久久久久久久久丰满| 久久精品夜夜夜夜夜久久蜜豆| 乱码一卡2卡4卡精品| 亚洲18禁久久av| 亚洲在线自拍视频| 激情 狠狠 欧美| 麻豆成人av视频| 夜夜爽天天搞| 我的老师免费观看完整版| 日韩 亚洲 欧美在线| 神马国产精品三级电影在线观看| 亚洲精品456在线播放app| 久久久成人免费电影| 国内精品久久久久精免费| 国产视频首页在线观看| 国产精品国产高清国产av| 久久久久九九精品影院| 成人毛片a级毛片在线播放| 久久99热6这里只有精品| avwww免费| 精品久久久久久久久久久久久| 久久精品影院6| 日韩成人伦理影院| 美女内射精品一级片tv| av.在线天堂| 狂野欧美白嫩少妇大欣赏| 亚洲欧美精品专区久久| 蜜桃久久精品国产亚洲av| 免费不卡的大黄色大毛片视频在线观看 | 亚洲va在线va天堂va国产| 久久久久久九九精品二区国产| 免费看av在线观看网站| 国产私拍福利视频在线观看| 亚洲成av人片在线播放无| 欧美在线一区亚洲| 亚洲电影在线观看av| 麻豆乱淫一区二区| avwww免费| 国产日韩欧美在线精品| 美女xxoo啪啪120秒动态图| 国产高清不卡午夜福利| 亚洲人与动物交配视频| 欧美bdsm另类| 黄色欧美视频在线观看| 观看美女的网站| 插逼视频在线观看| 99热这里只有是精品在线观看| 亚洲美女视频黄频| 日韩一区二区视频免费看| 亚洲av熟女| 欧美日韩一区二区视频在线观看视频在线 | 欧美极品一区二区三区四区| 国产男人的电影天堂91| 天天一区二区日本电影三级| 一区二区三区免费毛片| 日韩一本色道免费dvd| 国产精品久久电影中文字幕| 哪里可以看免费的av片| 日本熟妇午夜| 一本一本综合久久| 亚洲精品成人久久久久久| 亚洲人成网站高清观看| 日本av手机在线免费观看| av福利片在线观看| 久久精品夜色国产| 免费观看a级毛片全部| 国产色爽女视频免费观看| 亚洲av成人av| 少妇的逼水好多| 激情 狠狠 欧美| 日日摸夜夜添夜夜添av毛片| 亚洲一区二区三区色噜噜| 舔av片在线| 不卡视频在线观看欧美| 少妇猛男粗大的猛烈进出视频 | 国产老妇女一区| 男女视频在线观看网站免费| 九草在线视频观看| 两个人视频免费观看高清| 国产精品不卡视频一区二区| 深夜精品福利| 啦啦啦观看免费观看视频高清| 少妇的逼好多水| 久久这里有精品视频免费| av又黄又爽大尺度在线免费看 | 久久草成人影院| 少妇人妻一区二区三区视频| 一卡2卡三卡四卡精品乱码亚洲| 看非洲黑人一级黄片| 中出人妻视频一区二区| av国产免费在线观看| 99久国产av精品| 成人国产麻豆网| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 在线播放国产精品三级| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品亚洲一区二区| 国产淫片久久久久久久久| 国产高潮美女av| 少妇熟女欧美另类| 亚洲精品自拍成人| av黄色大香蕉| or卡值多少钱| 九九在线视频观看精品| 日本撒尿小便嘘嘘汇集6| 久久久久久久久大av| 乱码一卡2卡4卡精品| 综合色av麻豆| 九草在线视频观看| 久99久视频精品免费| 精品久久久久久久久亚洲| 十八禁国产超污无遮挡网站| 日日干狠狠操夜夜爽| 国内精品美女久久久久久| 成年av动漫网址| 舔av片在线| 日本黄大片高清| 免费黄网站久久成人精品| 我要看日韩黄色一级片| 亚洲中文字幕一区二区三区有码在线看| 人妻久久中文字幕网| 一级毛片电影观看 | 天天躁夜夜躁狠狠久久av| 精品人妻一区二区三区麻豆| 搡老妇女老女人老熟妇| 欧美精品一区二区大全| av黄色大香蕉| av在线播放精品| 日韩一本色道免费dvd| 最近最新中文字幕大全电影3| 少妇裸体淫交视频免费看高清| 我的老师免费观看完整版| 免费黄网站久久成人精品| 国产精品精品国产色婷婷| 午夜福利在线观看免费完整高清在 | 中文字幕av成人在线电影| 在线播放国产精品三级| 欧美一区二区亚洲| 国产爱豆传媒在线观看| 尤物成人国产欧美一区二区三区| 日韩成人av中文字幕在线观看| 一区二区三区高清视频在线| 日本在线视频免费播放| 天美传媒精品一区二区| 国产爱豆传媒在线观看|