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

    環(huán)境因子對(duì)沙地人工楊樹林樹干液流的驅(qū)動(dòng)影響*

    2021-05-22 06:43:36武鵬飛劉云強(qiáng)李冬梅陳志軍馬長(zhǎng)明
    中國(guó)農(nóng)業(yè)氣象 2021年5期
    關(guān)鍵詞:液流時(shí)滯樹干

    武鵬飛,劉云強(qiáng),李冬梅,陳志軍,馬長(zhǎng)明**

    環(huán)境因子對(duì)沙地人工楊樹林樹干液流的驅(qū)動(dòng)影響*

    武鵬飛1,劉云強(qiáng)2,李冬梅3,陳志軍3,馬長(zhǎng)明1**

    (1.河北農(nóng)業(yè)大學(xué)林學(xué)院,保定 071000;2.遼寧農(nóng)業(yè)職業(yè)技術(shù)學(xué)院,營(yíng)口 115009;3.河北省唐山市灤南縣國(guó)營(yíng)林場(chǎng),唐山 063502)

    以唐山市灤南縣林場(chǎng)107速生楊(cv.‘74/76’)為研究對(duì)象,使用滴灌法將土壤體積含水量分別設(shè)置為15%~20%、10%~15%和5%~10%,利用熱擴(kuò)散式探針法(TDP)測(cè)定楊樹邊材液流速率,同時(shí)監(jiān)測(cè)太陽輻射、氣溫、大氣濕度、大氣壓、風(fēng)速等氣象因子,探尋不同環(huán)境條件下107楊樹樹干日間和夜間的液流變化規(guī)律,分析不同環(huán)境下林木對(duì)水分的利用,為速生豐產(chǎn)林制定合理的水分灌溉策略提供指導(dǎo)。結(jié)果表明,晴天樹干液流較陰雨天啟動(dòng)時(shí)間早、峰值高、持續(xù)時(shí)間長(zhǎng)、夜間波動(dòng)大;土壤缺水條件下,樹干液流啟動(dòng)時(shí)間晚、峰值低、持續(xù)時(shí)間短,夜間液流變化平緩且增高;樹干液流與氣象因子存在“時(shí)滯效應(yīng)”,樹干液流滯后于太陽輻射40min,而提前于氣溫、相對(duì)濕度、VPD(水汽壓虧缺)60min;樹干液流與VPD呈現(xiàn)“遲滯回環(huán)”關(guān)系,時(shí)滯錯(cuò)位后氣象因子對(duì)樹干液流均具有顯著驅(qū)動(dòng)作用,日間樹干液流與氣溫相關(guān)性最大,驅(qū)動(dòng)力較強(qiáng),而夜間與VPD相關(guān)性最大,驅(qū)動(dòng)力較強(qiáng);經(jīng)過時(shí)滯錯(cuò)位后樹干液流與氣象因子擬合效果優(yōu)于時(shí)滯錯(cuò)位前。可見,環(huán)境因子對(duì)樹干液流具有明顯的驅(qū)動(dòng)作用,樹干液流可以較好地反映沙地土壤的水分情況和林木的水分消耗特征。

    樹干液流;熱擴(kuò)散探針;冀東沙地;107速生楊;環(huán)境因子

    樹干液流是土壤液態(tài)水進(jìn)入根系,由莖輸導(dǎo)組織運(yùn)送至冠層的過程中在樹體內(nèi)部產(chǎn)生的液體流動(dòng)。研究表明,樹干液流中蒸騰作用占到了90%以上,因此,研究樹干液流的動(dòng)態(tài)特征以及環(huán)境因子對(duì)其的驅(qū)動(dòng)作用,對(duì)探究植物自身體內(nèi)生理特性和水分動(dòng)態(tài)變化具有重要意義[1]。目前熱擴(kuò)散探針法(TDP)因操作簡(jiǎn)單、結(jié)果精準(zhǔn)等因素被廣泛應(yīng)用于林木耗水研究中[2],國(guó)內(nèi)外學(xué)者利用該方法針對(duì)不同樹種樹干液流進(jìn)行了研究,研究焦點(diǎn)主要集中于樹干液流的時(shí)間變化規(guī)律、空間變化規(guī)律以及影響因素等,如韓永貴等[3]指出樹干液流呈明顯季節(jié)性變化規(guī)律,黨宏忠等[4]研究指出沙地樟子松(var)樹干液流具有顯著的方位差異,Bhusal等[5]研究指出土壤水分條件對(duì)植物樹干液流起到一定限制作用。環(huán)境因子以及土壤條件是影響樹干液流的主要因素,然而,目前針對(duì)沙地土壤條件下環(huán)境因子以及土壤水分對(duì)樹干液流的驅(qū)動(dòng)影響等研究相對(duì)薄弱。

    107楊(cv.‘74/76’)具有生長(zhǎng)快、干形好、用途廣泛等特點(diǎn),是栽培廣泛的主要速生樹種之一,對(duì)改善生態(tài)環(huán)境、促進(jìn)經(jīng)濟(jì)建設(shè)、防風(fēng)固沙具有重要作用[6],同時(shí),107楊是水肥敏感型樹種,水分的多寡直接影響其生長(zhǎng)[7]。冀東沙地是遼河黃泛平原立地亞區(qū)的主要立地類型之一,沙壤土具有通透性好但保水能力差的特點(diǎn),更加凸顯了土壤水分的重要性。因此,本研究選擇冀東沙壤土立地類型,以107楊為研究對(duì)象,采用熱擴(kuò)散探針法,研究環(huán)境因子對(duì)樹干液流的驅(qū)動(dòng)影響,以及107楊潛在耗水規(guī)律,為沙地人工速生豐產(chǎn)林制定合理的水分灌溉策略提供指導(dǎo)。

    1 材料與方法

    1.1 研究區(qū)概況

    試驗(yàn)在河北省唐山市灤南縣(N38°48′23″,E115°24′58″,海拔20m)進(jìn)行,位于河北省東北部,地處燕山南麓灤河西岸,屬于冀東沙地地區(qū),灤南縣位于東部季風(fēng)區(qū),屬于暖溫帶半濕潤(rùn)季風(fēng)氣候。冬季受西伯利亞和蒙古冷空氣的影響,盛行偏北風(fēng);夏季受海洋氣團(tuán)和太平洋高壓的影響,盛行偏南風(fēng)。具有春季干燥多風(fēng)、夏季悶熱多雨、秋季晝暖夜寒、冬季寒冷少雪的特點(diǎn)。全年平均氣溫10.60℃,年平均日照2853h,年平均降水量658mm,無霜期186d。

    1.2 土壤情況

    利用Berrersize2000激光粒度分析儀對(duì)研究區(qū)土壤進(jìn)行分析,結(jié)果見表1,根據(jù)國(guó)際制土壤分級(jí)標(biāo)準(zhǔn),確定土壤類型為壤質(zhì)砂土,土壤基本特征見表2,該類型土壤基本呈無結(jié)構(gòu)狀態(tài),土壤大孔隙較多,透氣性和下滲性較好,但吸附作用較弱,保水與保肥能力較差。

    表1 研究區(qū)土壤粒徑分布

    表2 研究區(qū)土壤基本特性

    1.3 樣樹情況

    在唐山市灤南縣國(guó)有林場(chǎng)所轄的“107速生楊”人工林中,選取生長(zhǎng)狀況良好、林齡為4a的試驗(yàn)樹9棵,分別編號(hào),測(cè)量其胸徑、樹高、冠幅,結(jié)果見表3。

    表3 試驗(yàn)樹的胸徑、樹高和冠幅

    1.4 樹干液流測(cè)定

    采用熱擴(kuò)散原理測(cè)算樹干液流速率,利用熱擴(kuò)散式莖流計(jì)進(jìn)行測(cè)定,傳感器采用AV-3665R(北京產(chǎn))。探針規(guī)格為20mm,數(shù)據(jù)采集器型號(hào)為RR-1016。通過上下兩根探針的溫度差來計(jì)算液流速率,溫差越大,液流速率越小,反之溫差越小,液流速率越大。

    樹干液流測(cè)定時(shí)間為2020年4月初?10月底。在林木的南方向,距樹基部1.3m處安裝探針,在探針安裝處刮去粗糙的樹皮,然后采用特定規(guī)格的鉆頭,依據(jù)儀器自帶鉆模沿樹干縱向垂直鉆取直徑為2mm的孔洞。將探針插入樹木邊材,探針插好后,用1/4泡沫球?qū)⑻结樄潭ǎ饷姘^緣、防輻射材料,最后用膠帶密封,防止雨水進(jìn)入。將探針反饋線與數(shù)據(jù)采集器連接,用以自動(dòng)采集和記錄液流數(shù)據(jù),采樣頻率設(shè)定為60s,數(shù)據(jù)記錄間隔期為600s?;贕ranier原理計(jì)算液流速率的原始公式(Granier, 1985)為

    V=0.0119·K1.231(1)

    將液流速率轉(zhuǎn)化成莖流量,即

    FS=AS·V(3)

    式中,F(xiàn)s為液流通量(cm3·s?1);As為邊材面積(cm2)。

    1.5 邊材面積測(cè)定

    采用生長(zhǎng)錐法進(jìn)行測(cè)定。為避免對(duì)樣樹造成傷害從而影響樹干液流,在樣木附近選擇不同胸徑的107速生楊15棵,測(cè)量胸徑和樹皮厚度,利用生長(zhǎng)錐在測(cè)定樣木相同高度區(qū)域鉆取樹干木芯,將鉆得的木芯組織進(jìn)行染色處理,觀察其邊材和心材分布,測(cè)量并記錄統(tǒng)計(jì)邊材厚度,以此確定邊材面積。利用樹木胸徑與邊材面積的顯著相關(guān)性建立方程,計(jì)算得出每棵樣樹的邊材面積,即

    As=14.523D?65.857

    (R2=0.9919,P<0.01)(4)

    式中,D為胸徑(cm);As為邊材面積(cm?2);R2為方程的決定系數(shù)。

    1.6 土壤水分控制

    利用Spectrum TDR300土壤水分分析儀(北京產(chǎn)),測(cè)得土壤最大體積含水量即土壤中水分占有的體積與土壤總體積的比值為20.8%。

    采用滴灌方式灌水,即在樹干根部按照樹行方向鋪設(shè)滴灌管,每棵樹根部?jī)蓚?cè)設(shè)有2個(gè)滴頭。灌水期間(8月11?24日)每日進(jìn)行土壤含水量監(jiān)測(cè),通過控制灌水時(shí)間,將試驗(yàn)分為3個(gè)水分梯度,分別為15%~20%(用H表示)、10%~15%(用M表示)和5%~10%土壤體積含水量(用L表示)。

    1.7 氣象因子測(cè)定

    在林場(chǎng)內(nèi)安裝HOBO自動(dòng)氣象監(jiān)測(cè)站,自動(dòng)監(jiān)測(cè)并記錄氣溫(Ta,℃)、空氣相對(duì)濕度(RH,%)、太陽輻射強(qiáng)度(Rs,%)、大氣壓(kPa)以及風(fēng)速(m·s?1)等數(shù)據(jù)。數(shù)據(jù)采樣間隔為10min??諝馑麎禾澣保╒PD)計(jì)算式為

    VPD=0.611e[17.502Ta/(Ta+240.97)](1?RH)(5)

    式中,VPD為空氣水汽壓虧缺(kPa);Ta為氣溫(℃);RH為空氣相對(duì)濕度(%)。

    1.8 數(shù)據(jù)處理與統(tǒng)計(jì)

    選擇樹干液流主要驅(qū)動(dòng)因素VPD(空氣水汽壓虧缺)較高的8月[8],以全月1?31日數(shù)據(jù)進(jìn)行分析。采用Microsoft Excel對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行整理與計(jì)算,運(yùn)用SPSS 21.0對(duì)液流速率和氣象因子進(jìn)行相關(guān)分析,Origin 2018繪圖。

    2 結(jié)果與分析

    2.1 天氣類型對(duì)樹干液流的驅(qū)動(dòng)影響

    選取3個(gè)典型天氣,分別為晴天、陰天和雨天,分析不同天氣條件下107楊樹干液流的日變化差異,結(jié)果見圖1。由圖1a可以看出,在3種典型天氣條件下,日間(6:00?19:00)107楊樹干液流速率變化趨勢(shì)一致,均為先上升后下降的“幾”字型變化曲線,但不同天氣也存在一定的差異。晴天,樹干液流在6:00左右啟動(dòng),然后迅速上升,10:00?16: 30液流速率一直保持在較高水平(峰值為17.99cm·h?1),隨后開始迅速下降,持續(xù)時(shí)間約為6.5h;陰天,樹干液流啟動(dòng)時(shí)間在7:00左右,較晴天晚約1h,11:00?16:00樹干液流速率一直保持較高水平(峰值為9.85cm·h?1),持續(xù)約5h,隨后開始迅速下降,液流速率峰值較晴天時(shí)降低;雨天,樹干液流啟動(dòng)時(shí)間在7:00左右,隨后液流速率上升,但是維持在較低水平(峰值為5.76cm·h?1),持續(xù)時(shí)間約為5h??梢?,天氣類型對(duì)107楊日間樹干液流啟動(dòng)時(shí)間和液流速率峰值均有明顯影響,晴天液流啟動(dòng)時(shí)間比陰天和雨天早約1h,液流峰值和持續(xù)時(shí)間也大于陰天和雨天。由圖1b可見,3種典型天氣條件下,夜間(19:00?次日6:00)樹干液流速率均較低,夜間液流占全天液流比例多為5%~20%。夜間液流速率整體變化趨勢(shì)基本一致,呈現(xiàn)先下降后上升的“U”字型變化曲線,但從數(shù)值上看也存在一定的差異。晴天夜間樹干液流速率下降幅度大于陰天和雨天,雨天夜間液流速率變化幅度較為平緩;晴天夜間樹干液流速率的谷值(1.94cm·h?1)小于陰天(2.23cm·h?1)和雨天(2.24cm·h?1)。

    可見,天氣類型對(duì)樹干液流具有明顯的驅(qū)動(dòng)作用,晴天啟動(dòng)早、峰值高、持續(xù)時(shí)間長(zhǎng),夜間波動(dòng)大,驅(qū)動(dòng)力較強(qiáng),而陰天和雨天啟動(dòng)稍晚、峰值低、持續(xù)時(shí)間短,夜間波動(dòng)平緩,驅(qū)動(dòng)力較弱。

    2.2 土壤水分對(duì)樹干液流的驅(qū)動(dòng)影響

    土壤水分對(duì)樹干液流的影響較大,林木的液流啟動(dòng)時(shí)間、液流量以及峰值等時(shí)刻都受土壤水分的影響。為規(guī)避天氣類型的影響,選擇3d晴天且不同灌水處理下的液流數(shù)據(jù)進(jìn)行分析,結(jié)果見圖2。由圖2a可以看出,不同土壤水分條件下,日間(6:00?19:00)樹干液流速率存在一定差異。H(15%~20%土壤體積含水量)條件下,樹干液流啟動(dòng)時(shí)間在6:30左右,M(10%~15%土壤體積含水量)條件下,樹干液流啟動(dòng)時(shí)間在6:40左右,較H延遲10min啟動(dòng),而L(5%~10%土壤體積含水量)條件下,樹干液流啟動(dòng)最晚,在近7:00時(shí)。從樹干液流峰值來看,L條件下樹干液流速率峰值(8.00cm·h?1)低于H(15.49cm·h?1)和M(11.84cm·h?1),但其到達(dá)峰值的時(shí)間較前兩者要提前,呈現(xiàn)出樹干液流啟動(dòng)晚,到達(dá)峰值時(shí)間早,且峰值偏低的特征??梢?,土壤水分對(duì)日間樹干液流啟動(dòng)時(shí)間、液流量和液流速率峰值均有影響,15%~20%土壤體積含水量條件下樹干液流啟動(dòng)時(shí)間比10%~15%土壤體積含水量和5%~10%土壤體積含水量條件下均有所提前,樹干液流速率峰值也均比后兩者大。由圖2b可見,夜間(19:00?6:00)液流速率在不同土壤水分條件下也存在一定差異。H和M條件下夜間液流速率變化起伏程度要大于L條件;L條件下夜間液流速率的谷值(0.63cm·h?1)大于M(0.33cm·h?1)和H(0.12cm·h?1)條件,與日間液流速率峰值變化相反。

    由此可見,土壤水分對(duì)樹干液流具有明顯的驅(qū)動(dòng)作用,低水分條件下,樹干液流啟動(dòng)時(shí)間晚(7:00左右)、峰值低(8.00cm·h?1)、夜間波動(dòng)小,驅(qū)動(dòng)力較弱,而土壤水分充足時(shí),樹干液流啟動(dòng)早(6:30左右)、峰值高(15.49cm·h?1),夜間波動(dòng)大,驅(qū)動(dòng)力較強(qiáng)。

    注:H、M和L分別為15%~20%、10%~15%和5%~10%土壤體積含水量處理。

    Note: H, M and L were 15%?20%, 10%?15%, 5%?10% soil volume water content, respectively.

    2.3氣象因子對(duì)樹干液流的驅(qū)動(dòng)影響

    2.3.1 時(shí)滯分析

    研究表明,樹干液流隨著氣象因子的變化而變化,但是兩者存在一個(gè)時(shí)間差,這種現(xiàn)象就是時(shí)滯效應(yīng)。為規(guī)避天氣類型和土壤水分的影響,選擇晴天且供水充足條件下4d(8月7、14、21和22日)的107楊樹干液流數(shù)據(jù),按照觀測(cè)時(shí)間順序,將107楊樹干液流速率與氣溫、相對(duì)濕度、VPD和太陽輻射4項(xiàng)氣象因子數(shù)據(jù)對(duì)應(yīng)排列,按照10min的步長(zhǎng)開始錯(cuò)位移動(dòng),進(jìn)行相關(guān)性分析,直到相關(guān)系數(shù)的絕對(duì)值達(dá)到最大值并出現(xiàn)拐點(diǎn),結(jié)果見圖3。

    相關(guān)系數(shù)絕對(duì)值達(dá)到最大的時(shí)間,即樹干液流與氣象因子之間的實(shí)際時(shí)滯時(shí)間。從圖3可以看出,樹干液流速率與40min前太陽輻射的相關(guān)性最高,而與60min后氣溫、相對(duì)濕度、VPD的相關(guān)性最高。

    2.3.2 氣象因子驅(qū)動(dòng)

    選取氣象因子中主要的驅(qū)動(dòng)因子,即太陽輻射和水汽壓虧缺(VPD),經(jīng)時(shí)滯位移后與樹干液流進(jìn)行對(duì)應(yīng)分析,結(jié)果見圖4。由圖4a樹干液流速率與VPD響應(yīng)曲線中可以看出,一旦液流速率達(dá)到最大值,其變化趨勢(shì)將趨于穩(wěn)定,隨著驅(qū)動(dòng)因子進(jìn)一步增強(qiáng),樹干液流也不會(huì)再增加;而隨著驅(qū)動(dòng)因子的強(qiáng)度逐步減弱,樹干液流并非按比例立即減小而是逐漸衰減,并沿著不同的路徑后退,這種關(guān)系即為“遲滯回環(huán)”(hysteresis loops)關(guān)系,即上升曲線和下降曲線不重合。圖4b顯示,樹干液流對(duì)VPD的響應(yīng)回環(huán)過程呈順時(shí)針方向,而與太陽輻射的響應(yīng)過程并沒有表現(xiàn)出這種關(guān)系。

    注:氣溫(Ta)、空氣相對(duì)濕度(RH)、水汽壓虧缺(VPD)、太陽輻射強(qiáng)度(Rs),虛線表示相關(guān)系數(shù)達(dá)到最大值時(shí)所對(duì)應(yīng)的時(shí)滯時(shí)間。下同。

    Note: Ta is the air temperature, RH is the air relative humidity, Rs is the solar radiation intensity, VPD is the saturated vapor pressure difference. The dashed line indicate the time when correlation coefficient reached the maximum value. The same as below.

    將樹干液流速率與氣象因子進(jìn)行時(shí)滯錯(cuò)位,日間(6:00?19:00)和夜間(19:00?6:00)分別進(jìn)行相關(guān)性分析,結(jié)果如表4。由表可以看出,在日間,樹干液流與各個(gè)氣象因子之間均表現(xiàn)為極顯著相關(guān)性,尤其與空氣溫度和相對(duì)濕度相關(guān)性較高,相關(guān)系數(shù)絕對(duì)值達(dá)到0.940以上,其中與相對(duì)濕度呈負(fù)相關(guān)關(guān)系,而與太陽輻射和VPD之間的相關(guān)系數(shù)分別為0.901和0.923。在夜間,樹干液流與氣溫、相對(duì)濕度、太陽輻射以及VPD均表現(xiàn)為極顯著相關(guān)性,其中與相對(duì)濕度呈負(fù)相關(guān)關(guān)系,但是與氣溫、相對(duì)濕度、太陽輻射之間的相關(guān)系數(shù)絕對(duì)值均比日間減小,而與VPD的相關(guān)系數(shù)比日間增大。可見,氣象因子對(duì)樹干液流具有明顯的驅(qū)動(dòng)作用,日間氣溫、相對(duì)濕度和太陽輻射對(duì)樹干液流的影響比夜間大,而夜間VPD對(duì)樹干液流的影響比日間大;相對(duì)濕度與樹干液流之間呈負(fù)相關(guān)關(guān)系。

    采用逐步回歸法建立樹干液流速率與上述環(huán)境因子的多元線性模型,以5%和10%的可靠性作為因變量入選和剔除臨界值,最終的回歸模型參數(shù)見表5。由表可以看出,時(shí)滯錯(cuò)位前和時(shí)滯錯(cuò)位后的回歸方程有一定的差異,擬合效果也不一樣,經(jīng)過時(shí)滯錯(cuò)位后的擬合方程相關(guān)系數(shù)大,擬合效果好。因此,在分析氣象因子對(duì)樹干液流的驅(qū)動(dòng)影響時(shí),應(yīng)將各驅(qū)動(dòng)因子進(jìn)行時(shí)滯錯(cuò)位后再進(jìn)行樹干液流回歸模型的構(gòu)建,可以有效減小誤差。

    表4 時(shí)滯錯(cuò)位后樹干液流與氣象因子的相關(guān)性

    注:**表示P<0.01。

    Note:**is P<0.01.

    表5 樹干液流與氣象因子的多元回歸模型估計(jì)

    注:SF為樹干液流,R2為方程的決定系數(shù),P為顯著性檢驗(yàn)。

    Note: SF is the tree trunk sap flow, R2is the coefficient, and P is the significance test level.

    3 討論與結(jié)論

    3.1 討論

    大量研究表明,環(huán)境因子對(duì)樹木的生長(zhǎng)發(fā)育存在重要的影響,尤其對(duì)樹干液流的驅(qū)動(dòng)影響巨大[9]。對(duì)冀東沙地不同天氣類型條件下107楊人工林樹干液流變化規(guī)律的研究表明,晴天樹干液流速率較陰雨天啟動(dòng)時(shí)間早、峰值高、持續(xù)時(shí)間長(zhǎng)、夜間波動(dòng)大,這與賈天宇等[10?11]研究指出在整個(gè)生長(zhǎng)季內(nèi)樹干液流速率晴天明顯大于陰雨天,趙春彥等[12]研究指出胡楊()陰雨天樹干液流啟動(dòng)時(shí)間較晴天推遲,樹干液流停止時(shí)間提前,液流歷時(shí)縮短,且峰值顯著縮小等研究結(jié)果一致。出現(xiàn)這種結(jié)果的主要原因是晴天太陽輻射強(qiáng),溫度高,空氣濕度低,導(dǎo)致植物蒸騰量增大,從而使樹干液流增大,而陰雨天太陽輻射弱,溫度低,空氣濕度高,植物蒸騰減弱,液流也隨之減小。有學(xué)者指出,樹干液流在中午左右會(huì)出現(xiàn)“午休”現(xiàn)象,即中午左右樹體為了保存水分而降低蒸騰,使樹干液流減小[13]。然而本研究并未發(fā)現(xiàn)此現(xiàn)象,這可能與107楊為了適應(yīng)沙地土壤環(huán)境做出的自我水分調(diào)控策略有關(guān)。

    夜間液流對(duì)樹木彌補(bǔ)水分虧缺和維持自身水分平衡有重要意義[14]。從研究結(jié)果中可以看出,樹干夜間液流非常微弱,但并不為0,這與劉瀟瀟等[15]研究結(jié)果一致。對(duì)于夜間液流的存在,有的學(xué)者認(rèn)為是樹木夜間補(bǔ)水所致,也有學(xué)者認(rèn)為是樹木夜間蒸騰所致,而鄢春華等[16]研究認(rèn)為夜間樹木莖干補(bǔ)水和夜間蒸騰同時(shí)發(fā)生。針對(duì)夜間液流的組成和精確測(cè)量等,仍然有較大的研究意義。

    土壤水分主要從潛在水平影響樹干液流從而影響林木的生長(zhǎng)。本研究通過設(shè)置不同水分梯度,研究土壤水分對(duì)樹干液流的驅(qū)動(dòng)影響,結(jié)果表明土壤在缺水條件下,樹干液流啟動(dòng)時(shí)間晚,減弱時(shí)間早,峰值降低,而夜間樹干液流速率谷值卻增高,并且變化幅度較為平緩,這與王媛等[17]研究表明不同土壤含水量下白樺氣孔蒸騰作用對(duì)土壤含水量變化具有較高的敏感性,吳鵬等[18]研究表明土壤含水量的變化是導(dǎo)致雨天夜間液流大于晴天和陰天的主要因素等結(jié)果一致。然而劉春鵬等[19]研究認(rèn)為土壤含水量是影響雨天的主要環(huán)境因子,對(duì)晴天和陰天的樹干液流影響不大,這與張曉艷等[20]研究指出土壤含水量主要影響晴天樹干液流,而與陰雨天相關(guān)性不顯著的結(jié)論相反。出現(xiàn)這種現(xiàn)象的原因可能是不同樹種對(duì)不同土壤和水分條件產(chǎn)生了自我水分利用機(jī)制。

    氣象因子主要是瞬態(tài)特征影響樹干液流。在擬合樹干液流與氣象因子模型中,考慮兩者之間的時(shí)滯效應(yīng)對(duì)提高擬合精度至關(guān)重要。本研究指出樹干與氣象因子存在明顯時(shí)滯效應(yīng),即樹干液流滯后于太陽輻射40min,而提前于氣溫、相對(duì)濕度、VPD60min,這與李少宇等[21]研究指出樹干液流與太陽輻射、氣溫、相對(duì)濕度存在明顯滯后效應(yīng)結(jié)果一致。產(chǎn)生這一現(xiàn)象的原因是蒸騰散失水分主要從樹冠開始,樹干基部水分需要等樹冠存儲(chǔ)的水分消耗到一定程度產(chǎn)生梯度差才開始向上運(yùn)輸,兩者之間產(chǎn)生的時(shí)間差導(dǎo)致了時(shí)滯效應(yīng)。針對(duì)時(shí)滯效應(yīng)不同學(xué)者在研究中存有不同的觀點(diǎn),楊潔等[22]研究指出樹干液流與環(huán)境因子的時(shí)滯效應(yīng)與樹干大小有關(guān),與樹種無關(guān),而姚增旺等[23]研究指出樹干液流與環(huán)境因子的時(shí)滯效應(yīng)與樹高及樹冠的暴露程度無關(guān)。

    時(shí)滯錯(cuò)位后,無論日間還是夜間,樹干液流與氣象因子之間均存在極顯著相關(guān)性。日間樹干液流與空氣溫度相關(guān)性最高,而夜間則與水汽壓虧缺相關(guān)系數(shù)最高,這與張婕等[14,24?25]研究結(jié)果一致。氣象因子對(duì)樹干液流的驅(qū)動(dòng)主要在于太陽輻射,太陽輻射增高,空氣溫度升高,大氣濕度降低,飽和水汽壓也隨之降低,導(dǎo)致樹木蒸騰增強(qiáng),從而使樹干液流增大。

    樹木的生長(zhǎng)年限不同,樹干液流的變化規(guī)律也有可能不同;樹木不同徑向上樹干液流也存在一定的差異,這會(huì)導(dǎo)致對(duì)樹木耗水以及與氣象因子之間的關(guān)系研究上存在一定誤差[26],然而,近些年對(duì)這些方向的研究仍有欠缺。因此,對(duì)樹木整個(gè)生長(zhǎng)季進(jìn)行研究,并分析不同生長(zhǎng)年限、不同徑向上的樹干液流變化規(guī)律以及對(duì)環(huán)境因子的響應(yīng)都具有重大的研究意義。

    3.2 結(jié)論

    環(huán)境因子對(duì)沙地土壤條件下樹干液流存在重要的驅(qū)動(dòng)影響。晴天天氣下,樹干液流速率較陰雨天啟動(dòng)早、峰值高、持續(xù)時(shí)間長(zhǎng)、夜間波動(dòng)大,驅(qū)動(dòng)力較強(qiáng);土壤低水分條件下,樹干液流啟動(dòng)晚、峰值低、夜間波動(dòng)小,驅(qū)動(dòng)力較弱;樹干液流與氣象因子之間存在“時(shí)滯效應(yīng)”,即樹干液流速率與40min前太陽輻射的相關(guān)性最高,而與60min后氣溫、相對(duì)濕度、水汽壓虧缺的相關(guān)性最高;將氣象因子時(shí)滯錯(cuò)位后與樹干液流進(jìn)行回歸模型擬合效果好,該模型更能說明氣象因子與樹干液流之間的數(shù)學(xué)關(guān)系,為楊樹速生豐產(chǎn)林的精準(zhǔn)灌溉策略提供有力依據(jù)。

    [1]楊瑞,喻理飛,戴全厚.喀斯特區(qū)云南鼠刺樹干液流及土壤水分動(dòng)態(tài)[J].水土保持通報(bào),2015,35(2):89-93.

    Yang R,Yu L F,Dai Q H.Stem sap flow and soil moisture dynamics ofin Karst area[J].Bulletin of Soil and Water Conservation,2015,35(2):89-93.(in Chinese)

    [2]任啟文,畢君,李聯(lián)地,等.冀北山地油松樹干液流特征及其對(duì)環(huán)境因子的響應(yīng)[J].水土保持學(xué)報(bào),2017,31(5): 183-188.

    Ren Q W,Bi J,Li L D,et al.The characteristics of trunk sap flow ofin the mountainous area of northern Hebei and its response to environmental factors[J].Journal of Soil and Water Conservation,2017, 31(5):183-188.(in Chinese)

    [3]韓永貴,高陽,韓磊,等.基于模糊規(guī)則的人工神經(jīng)網(wǎng)絡(luò)模擬新疆楊蒸騰耗水[J].應(yīng)用生態(tài)學(xué)報(bào),2020,31(5): 1525-1534.

    Han Y G,Gao Y,Han L,et al.Artificial neural network based on fuzzy rules to simulate transpiration of Xinjiang poplar[J].Chinese Journal of Applied Ecology,2020,31(5): 1525-1534.(in Chinese)

    [4]黨宏忠,馮金超,韓輝.沙地樟子松邊材液流速率的方位差異特征[J].林業(yè)科學(xué),2020,56(1):29-37.

    Dang H Z,Feng J C,Han H.The azimuth difference characteristics of the sapwood sap flow rate ofon sandy land[J].Scientia Silvae Sinicae,2020, 56(1):29-37.(in Chinese)

    [5]Bhusal N,Han S G,Yoon T M.Impact of drought stress on photosynthetic response,leaf water potential,and stem sap flow in two cultivars of bi-leader apple trees(×Borkh.)[J].Scientia Horticulturae,2019,246: 535-543.

    [6]秘洪雷,蘭再平,孫尚偉,等.滴灌栽培楊樹人工林細(xì)根空間分布特征[J].林業(yè)科學(xué)研究,2017,30(6):946-953.

    Mi H L,Lan Z P,Sun S W,et al.Spatial distribution characteristics of fine roots in poplar plantations cultivated by drip irrigation[J].Forest Research,2017,30(6):946-953. (in Chinese)

    [7]Yong H,Zaiping L,Shangwei S,et al.Effects of drip irrigation on the growth and use efficiency of water and fertilizer of Young'107' poplar plantation[J].Journal of Northeast Forestry University,2015,43(11):37-41.

    [8]Pasqualotto G,Carraro V,Menardi R,et al.Calibration of Granier-Type(TDP) sap flow probes by a high precision electronic potometer[J].Sensors,2019,19(10):2419.

    [9]崔鴻俠,唐萬鵬,胡文杰,等.神農(nóng)架2樹種樹干液流特征及與環(huán)境因子關(guān)系[J].西北林學(xué)院學(xué)報(bào),2020,35(1):37-44.

    Cui H X,Tang W P,Hu W J,et al.The characteristics of trunk sap flow of two tree species in Shennongjia and its relationship with environmental factors[J].Journal of Northwest Forestry University,2020,35(1):37-44.(in Chinese)

    [10]賈天宇,劉廷璽,段利民,等.半干旱沙丘草甸過渡帶人工楊樹蒸騰耗水規(guī)律[J].生態(tài)學(xué)雜志,2020,39(10): 3255-3264.

    Jia T Y,Liu Y X,Duan L M,et al.Transpiration and water consumption of artificial poplars in the semi-arid dune-meadow transition zone[J].Chinese Journal of Ecology,2020,39(10):3255-3264.(in Chinese)

    [11]高峻,孟平,張勁松,等.華北石質(zhì)山區(qū)杏樹蒸騰規(guī)律和蒸降差分析[J].中國(guó)農(nóng)業(yè)氣象,2009,30(4):538-542.

    Gao J,Meng P,Zhang J S,et al.Analysis of transpiration law and evapotranspiration difference of apricot trees in rocky mountain area of North China[J].Chinese Journal of Agrometeorology,2009,30(4):538-542.(in Chinese)

    [12]趙春彥,司建華,馮起,等.胡楊()樹干液流特征及其與環(huán)境因子的關(guān)系[J].中國(guó)沙漠,2014, 34(3):718-724.

    Zhao C Y,Si J H,Feng Q,et al.trunk sap flow characteristics and its relationship with environmental factors[J].Journal of Desert Research,2014,34(3):718-724. (in Chinese)

    [13]凌錦霞,史彥江,羅達(dá),等.平歐雜種榛樹干莖流對(duì)環(huán)境因子的響應(yīng)[J].新疆農(nóng)業(yè)科學(xué),2020,57(5):785-797.

    Ling J X,Shi Y J,Luo D,et al.Response of stem flow of Ping'ou hybrid hazel to environmental factors[J].Xinjiang Agricultural Sciences,2020,57(5):785-797.(in Chinese)

    [14]張婕,蔡永茂,陳立欣,等.北京山區(qū)元寶楓夜間液流活動(dòng)特征及影響因素[J].生態(tài)學(xué)報(bào),2019,39(9):3210-3223.

    Zhang J,Cai Y M,Chen L X,et al.Characteristics and influencing factors of night sap flow ofin Beijing mountainous area[J].Acta Ecologica Sinica,2019, 39(9):3210-3223.(in Chinese)

    [15]劉瀟瀟,李國(guó)慶,閆美杰,等.黃土高原主要樹種樹干液流研究進(jìn)展[J].水土保持研究,2017,24(3):369-373.

    Liu X X,Li G Q,Yan M J,et al.Research progress on trunk sap flow of main tree species in the Loess Plateau[J]. Research of Soil and Water Conservation,2017,24(3): 369-373.(in Chinese)

    [16]鄢春華,王蓓,鄒振東,等.九寨溝針闊混交林的夜間液流及其分配特征研究[J].北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020,56(4):732-738.

    Yan C H,Wang B,Zou Z D,et al.Study on the night sap flow and its distribution characteristics of Jiuzhaigou conifer and broadleaf mixed forest[J].Acta Scientiarum Naturalium Universitatis Pekinensis,2020,56(4):732-738.(in Chinese)

    [17]王媛,魏江生,周梅,等.大興安嶺南段白樺樹干液流對(duì)土壤水分的響應(yīng)[J].水土保持研究,2020,27(4):128-133.

    Wang Y,Wei J S,Zhou M,et al.Response of birch sap flow to soil moisture in the southern part of Daxinganling[J]. Research of Soil and Water Conservation,2020,27(4): 128-133.(in Chinese)

    [18]吳鵬,楊文斌,崔迎春,等.喀斯特區(qū)天峨槭()樹干液流特征及其與環(huán)境因子的相關(guān)分析[J].生態(tài)學(xué)報(bào),2017,37(22):7552-7567.

    Wu P,Yang W B,Cui Y C,et al.The characteristics oftrunk sap flow and its correlation with environmental factors in Karst area[J].Acta Ecologica Sinica,2017,37(22):7552-7567.(in Chinese)

    [19]劉春鵬,朱美秋,薛帥,等.干旱條件下荊條樹干液流對(duì)不同天氣的響應(yīng)[J].福建林業(yè)科技,2014,41(1):34-39.

    Liu C P,Zhu M Q,Xue S,et al.Responses of vitex twig sap flow to different weathers under drought conditions[J]. Journal of Fujian Forestry Science and Technology,2014, 41(1):34-39.(in Chinese)

    [20]張曉艷,褚建民,孟平,等.民勤綠洲荒漠過渡帶梭梭[C.A.Mey)Bunge]樹干液流特征及其對(duì)環(huán)境因子的響應(yīng)[J].生態(tài)學(xué)報(bào),2017,37(5): 1525-1536.

    Zhang X Y,Zhu J M,Meng P,et al.The characteristics ofC.A.Mey Bunge trunk sap flow and its response to environmental factors in the Minqin oasis-desert transition zone[J].Acta Ecologica Sinica,2017, 37(5):1525-1536.(in Chinese)

    [21]李少寧,魯紹偉,趙云閣,等.北京典型天氣下的4種闊葉樹種液流特征及其影響因素[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2019,35(2):189-196.

    Li S Y,Lu S W,Zhao Y G,et al.Sap flow characteristics and influencing factors of four broad-leaved tree species under typical weather in Beijing[J].Journal of Ecology and Rural Environment,2019,35(2):189-196.(in Chinese)

    [22]楊潔,呂金林,何秋月,等.黃土丘陵區(qū)遼東櫟和刺槐樹干液流時(shí)滯效應(yīng)與蒸騰特征的關(guān)聯(lián)性[J].應(yīng)用生態(tài)學(xué)報(bào),2019,30(8):2607-2613.

    Yang J,Lv J L,He Q Y,et al.Correlation between time lag effects of trunk sap flow and transpiration characteristics ofandin the hilly Loess Plateau[J].Chinese Journal of Applied Ecology, 2019,30(8):2607-2613.(in Chinese)

    [23]姚增旺,褚建民,吳利祿,等.民勤綠洲荒漠過渡帶梭梭樹干液流的時(shí)滯特征[J].應(yīng)用生態(tài)學(xué)報(bào),2018,29(7): 2339-2346.

    Yao Z W,Zhu J M,Wu L L,et al.The time lag characteristics of the stem sap flow ofin the Minqin oasis-desert transition zone[J].Chinese Journal of Applied Ecology,2018,29(7):2339-2346.(in Chinese)

    [24]孔喆,陳勝楠,律江,等.歐美楊單株液流晝夜組成及其影響因素分析[J].林業(yè)科學(xué),2020,56(3):8-20.

    Kong Z,Chen S N,Lv J,et al.Analysis of the circadian and influencing factors of the sap flow of individual poplar trees in Europe and America[J].Scientia Silvae Sinicae, 2020,56(3):8-20.(in Chinese)

    [25]王志超,許宇星,竹萬寬,等.雷州半島尾葉桉人工林夜間耗水特征及驅(qū)動(dòng)因素[J].浙江農(nóng)林大學(xué)學(xué)報(bào),2020, 37(4):646-653.

    Wang Z C,Xu Y X,Zhu W K,et al.Characteristics and driving factors of night water consumption ofplantation in Leizhou Peninsula[J].Journal of Zhejiang A&F University,2020,37(4):646-653.(in Chinese)

    [26]劉超,李春友,張勁松,等.生長(zhǎng)季大葉女貞樹干液流的徑向分布特征[J].中國(guó)農(nóng)業(yè)氣象,2011,32(2):196-202.

    Liu C,Li C Y,Zhang J S,et al.The radial distribution characteristics of the stem sap flow ofin the growing season[J].Chinese Journal Agricultural Meteorology,2011,32(2):196-202.(in Chinese)

    Driving Influence of Environmental Factors on the Sap Flow of the Artificial Poplar Forest on Sandy Land

    WU Peng-fei1, LIU Yun-qiang2, LI Dong-mei3, CHEN Zhi-jun3, MA Chang-ming1

    (1.College of Forestry, Agricultural University of Hebei, Baoding 071000, China;2.Liaoning Agricultural Technical College, Yingkou 115009; 3.State-owned Forest Farm of Luannan County, Tangshan City, Hebei Province,Tangshan 063502)

    Stem sap flow is an important basis to characterize the transpiration of trees. Studying the dynamic characteristics of sap flow and the driving effects of environmental factors on it is of great significance to explore the physiological characteristics and water dynamic changes of trees.cv.‘74/76’ is one of the main fast-growing tree species widely cultivated and also a water-fertilizer-sensitive tree species. The amount of water directly affects its growth.cv.‘74/76’ in Luannan county forest farm in Tangshan city was taken as research object, the drip irrigation method was used to set the soil volumetric water content at 15%?20%, 10%?15% and 5%?10% respectively by controlling the irrigation time, the soil volumetric water content was respectively set as 15%?20%, 10%?15% and 5%?10% by drip irrigation, the sapwood flow rate ofcv.‘74/76’ was measured by the thermal dissipation probe(TDP), and the solar radiation, air temperature, atmospheric humidity, atmospheric pressure, wind speed and other meteorological factors were simultaneously monitored by using the HOBO automatic weather monitoring station to explore the change rules of sap flow during day and night under different environmental conditions, establish a regression model between meteorological factors and the sap flow, analyze the water utilization of trees in different environments in order to provide guidance for the formulation of reasonable water irrigation strategies for fast-growing and high-yield plantation. Tharted earlier, peaked higher, lasted longer, and fluctuated more at night than cloudy and rainy de results showed that the sap flow was driven by different weather conditions, it on sunny days stays; The sap flow was also driven by soil moisture, under the condition of soil water shortage, the start time was late, the peak was low, the duration was short, and it changed steadily and increasingly at night. After time lag analysis, it is found that there is a "time lag effect" between the sap flow and meteorological factors, the sap flow lagged behind the solar radiation for 40 min, and was ahead of temperature, relative humidity and VPD (vapor pressure deficit) for 60min. There was a "hysteresis loop" relationship between sap flow and VPD, meteorological factors all had a significantly driving effect on the sap flow after the time lag dislocation, it had the greatest correlation with air temperature and stronger driving force during the day, while had the greatest correlation with VPD and stronger driving force at night. After time lag dislocation, the imitative effect of sap flow and meteorological factors was better than before. It can be seen that the environmental factors had an obvious driving effect on the sap flow, especially weather types, meteorological factors, and soil moisture, the regression model established between the weather factors and the sap flow after the time lag can better explain the relationship between the two mathematical relationship,meanwhile, sap flow can better reflect the water condition of sandy soil and the water consumption characteristics of trees, which provided a strong basis for the precise irrigation strategy of poplar fast-growing and high-yielding forests.

    Sap flow; Thermal dissipation probe; Sandy land in eastern Hebei;cv.‘74/76’; Environmental factors

    10.3969/j.issn.1000-6362.2021.05.005

    武鵬飛,劉云強(qiáng),李冬梅,等. 環(huán)境因子對(duì)沙地人工楊樹林樹干液流的驅(qū)動(dòng)影響[J].中國(guó)農(nóng)業(yè)氣象,2021,42(5):402-411

    2020?10?19

    “十三五”國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFD0600401)

    馬長(zhǎng)明,副教授,從事森林培育研究,E-mail:machangming@126.com

    武鵬飛,E-mail:2542790171@qq.com

    猜你喜歡
    液流時(shí)滯樹干
    為什么樹干不是方的?
    軍事文摘(2021年16期)2021-11-05 08:49:16
    帶有時(shí)滯項(xiàng)的復(fù)Ginzburg-Landau方程的拉回吸引子
    一種全釩液流電池并網(wǎng)控制系統(tǒng)的仿真設(shè)計(jì)
    為什么要在樹干上刷白漿
    為什么要在樹干上刷白漿
    為什么要在樹干上刷一層白漿
    H2-Fe3+/Fe2+氧化還原液流電池初步研究
    一階非線性時(shí)滯微分方程正周期解的存在性
    一類時(shí)滯Duffing微分方程同宿解的存在性
    我國(guó)—項(xiàng)液流電池國(guó)際標(biāo)準(zhǔn)提案獲準(zhǔn)通過
    色综合婷婷激情| 国产69精品久久久久777片 | 91九色精品人成在线观看| 一个人看的www免费观看视频| 12—13女人毛片做爰片一| 午夜a级毛片| 国产亚洲欧美在线一区二区| 国产欧美日韩一区二区精品| 五月玫瑰六月丁香| 九九热线精品视视频播放| 丝袜人妻中文字幕| www.精华液| 麻豆一二三区av精品| 一个人看的www免费观看视频| 国产成人系列免费观看| 日本a在线网址| 三级毛片av免费| 真实男女啪啪啪动态图| 国产精品1区2区在线观看.| 女人高潮潮喷娇喘18禁视频| 午夜福利免费观看在线| 亚洲人成电影免费在线| 免费观看精品视频网站| 久久精品91蜜桃| 亚洲欧美日韩无卡精品| 国产三级黄色录像| 国产精品美女特级片免费视频播放器 | 淫妇啪啪啪对白视频| 亚洲国产精品合色在线| 国产一区在线观看成人免费| а√天堂www在线а√下载| 黄色片一级片一级黄色片| 亚洲欧洲精品一区二区精品久久久| 国产亚洲精品综合一区在线观看| 国产三级中文精品| 听说在线观看完整版免费高清| www.www免费av| 亚洲avbb在线观看| 国产精品国产高清国产av| 精品人妻1区二区| 国产高清视频在线观看网站| 亚洲成av人片免费观看| 高潮久久久久久久久久久不卡| 国产视频一区二区在线看| 午夜福利欧美成人| 麻豆av在线久日| 又爽又黄无遮挡网站| 人人妻人人看人人澡| 免费在线观看影片大全网站| 国产精品一区二区三区四区久久| 欧美色视频一区免费| 国产又黄又爽又无遮挡在线| 全区人妻精品视频| 亚洲av成人av| 男插女下体视频免费在线播放| 美女扒开内裤让男人捅视频| 麻豆成人av在线观看| 99re在线观看精品视频| 免费人成视频x8x8入口观看| 高清在线国产一区| 欧美三级亚洲精品| 美女高潮喷水抽搐中文字幕| 中国美女看黄片| 法律面前人人平等表现在哪些方面| 老司机深夜福利视频在线观看| 免费看a级黄色片| 女生性感内裤真人,穿戴方法视频| 毛片女人毛片| 级片在线观看| 国产黄片美女视频| 在线免费观看的www视频| 中国美女看黄片| 级片在线观看| 亚洲成人免费电影在线观看| 亚洲熟妇中文字幕五十中出| 99久久精品国产亚洲精品| 一级毛片高清免费大全| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 国产成人av教育| 一级毛片精品| 人妻夜夜爽99麻豆av| 免费av毛片视频| 免费人成视频x8x8入口观看| 久久人妻av系列| 在线a可以看的网站| 在线观看一区二区三区| 国产成人影院久久av| 欧美国产日韩亚洲一区| 在线观看66精品国产| 亚洲美女视频黄频| 美女cb高潮喷水在线观看 | 国产一区二区激情短视频| 麻豆国产av国片精品| 亚洲在线观看片| 国产精品98久久久久久宅男小说| 99精品欧美一区二区三区四区| 桃红色精品国产亚洲av| 丁香六月欧美| 人人妻,人人澡人人爽秒播| 淫秽高清视频在线观看| 亚洲欧美精品综合久久99| 不卡av一区二区三区| av天堂在线播放| 精品久久久久久久毛片微露脸| 精品无人区乱码1区二区| 亚洲午夜理论影院| 啦啦啦免费观看视频1| 69av精品久久久久久| 欧美中文综合在线视频| 国产视频一区二区在线看| 18禁国产床啪视频网站| 欧美日韩一级在线毛片| 香蕉久久夜色| 99re在线观看精品视频| 99riav亚洲国产免费| 亚洲人成网站高清观看| 国产精品久久久久久亚洲av鲁大| 欧美黑人巨大hd| 日本在线视频免费播放| av在线蜜桃| 亚洲熟妇熟女久久| 黄色片一级片一级黄色片| 香蕉久久夜色| 一a级毛片在线观看| 亚洲精华国产精华精| 午夜久久久久精精品| 日本黄色视频三级网站网址| 视频区欧美日本亚洲| 亚洲欧美激情综合另类| 国产黄片美女视频| 亚洲国产欧美人成| 午夜福利成人在线免费观看| 全区人妻精品视频| 成人特级av手机在线观看| 桃色一区二区三区在线观看| h日本视频在线播放| 一区福利在线观看| 香蕉丝袜av| 男女床上黄色一级片免费看| 欧美激情在线99| 成人性生交大片免费视频hd| 九色国产91popny在线| 国产精品久久久人人做人人爽| 国产激情久久老熟女| 老司机深夜福利视频在线观看| 美女免费视频网站| 色精品久久人妻99蜜桃| 国产亚洲av高清不卡| 亚洲人成伊人成综合网2020| 亚洲国产看品久久| 欧美另类亚洲清纯唯美| 亚洲av第一区精品v没综合| 国产成人系列免费观看| av黄色大香蕉| 久久久久久大精品| 亚洲国产欧美网| 99国产综合亚洲精品| 高潮久久久久久久久久久不卡| 日本撒尿小便嘘嘘汇集6| 精品日产1卡2卡| 国产综合懂色| 99久久国产精品久久久| www日本在线高清视频| 偷拍熟女少妇极品色| 国产精品久久久久久亚洲av鲁大| 午夜日韩欧美国产| 久久九九热精品免费| 成年女人看的毛片在线观看| 丰满人妻熟妇乱又伦精品不卡| 色精品久久人妻99蜜桃| 男女下面进入的视频免费午夜| 久久久久免费精品人妻一区二区| 欧美中文日本在线观看视频| 黄色视频,在线免费观看| 欧美绝顶高潮抽搐喷水| 国产精品爽爽va在线观看网站| 久久久水蜜桃国产精品网| 午夜精品在线福利| 日韩大尺度精品在线看网址| 亚洲av成人不卡在线观看播放网| 日韩欧美免费精品| 成人三级做爰电影| av天堂中文字幕网| 国内揄拍国产精品人妻在线| 国产精品久久视频播放| 九色国产91popny在线| 九九在线视频观看精品| 色吧在线观看| 精品欧美国产一区二区三| 亚洲av成人精品一区久久| 久久香蕉精品热| 五月玫瑰六月丁香| 日韩精品中文字幕看吧| 在线观看一区二区三区| 精品久久久久久,| 午夜福利成人在线免费观看| svipshipincom国产片| 欧美日韩福利视频一区二区| 午夜福利欧美成人| 亚洲国产精品合色在线| 老司机福利观看| 国产伦精品一区二区三区视频9 | 一级毛片高清免费大全| 国产高清videossex| 中文字幕精品亚洲无线码一区| 久久久久九九精品影院| 18禁美女被吸乳视频| www.999成人在线观看| 欧美日韩国产亚洲二区| 亚洲国产欧美人成| 又黄又爽又免费观看的视频| 51午夜福利影视在线观看| 看黄色毛片网站| 我要搜黄色片| 亚洲欧美精品综合一区二区三区| 99国产综合亚洲精品| 欧美高清成人免费视频www| 深夜精品福利| 亚洲人成网站高清观看| 国产三级在线视频| 久久久久国产一级毛片高清牌| 999久久久国产精品视频| 搡老妇女老女人老熟妇| 久久香蕉精品热| aaaaa片日本免费| 国产乱人视频| 一本一本综合久久| 在线a可以看的网站| 黄片大片在线免费观看| 亚洲最大成人中文| 亚洲五月婷婷丁香| 一a级毛片在线观看| 欧美乱码精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 99久久成人亚洲精品观看| av欧美777| 99久久综合精品五月天人人| 九九久久精品国产亚洲av麻豆 | 精品久久久久久久末码| 一进一出好大好爽视频| 不卡一级毛片| 黄色成人免费大全| 国产精品一区二区免费欧美| 国产人伦9x9x在线观看| 国产成人精品久久二区二区免费| 国产亚洲欧美98| 成在线人永久免费视频| 久久久久久久久久黄片| 日韩欧美三级三区| 国产欧美日韩精品亚洲av| 琪琪午夜伦伦电影理论片6080| 在线观看美女被高潮喷水网站 | 成年人黄色毛片网站| 青草久久国产| 一进一出抽搐动态| 女的被弄到高潮叫床怎么办| 中文天堂在线官网| 99国产精品一区二区蜜桃av| 少妇的逼水好多| 亚洲欧美成人综合另类久久久 | 国产真实伦视频高清在线观看| 18禁裸乳无遮挡免费网站照片| 91精品一卡2卡3卡4卡| 亚洲怡红院男人天堂| 好男人在线观看高清免费视频| 亚洲国产日韩欧美精品在线观看| 日本猛色少妇xxxxx猛交久久| 18禁动态无遮挡网站| 国产免费一级a男人的天堂| 久久综合国产亚洲精品| av在线老鸭窝| 男人舔奶头视频| 色播亚洲综合网| 可以在线观看毛片的网站| 中文乱码字字幕精品一区二区三区 | 精品无人区乱码1区二区| 日韩,欧美,国产一区二区三区 | 精品欧美国产一区二区三| av在线老鸭窝| 天堂中文最新版在线下载 | 欧美另类亚洲清纯唯美| 国产精品人妻久久久久久| 欧美一区二区精品小视频在线| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品成人久久小说| 久久久久久久久中文| av在线蜜桃| 在线免费观看的www视频| 成人无遮挡网站| 高清毛片免费看| 少妇高潮的动态图| 日韩大片免费观看网站 | av线在线观看网站| 在线a可以看的网站| 久久亚洲精品不卡| 亚洲精品国产av成人精品| 欧美另类亚洲清纯唯美| 最近中文字幕2019免费版| 久久这里只有精品中国| 丰满少妇做爰视频| 天堂网av新在线| 中文亚洲av片在线观看爽| 国产精品一及| 国产在视频线精品| 成人综合一区亚洲| 日本与韩国留学比较| 日本一二三区视频观看| 国产精品.久久久| 中文字幕熟女人妻在线| 久久这里有精品视频免费| 亚洲伊人久久精品综合 | 国产女主播在线喷水免费视频网站 | 天天躁日日操中文字幕| 色5月婷婷丁香| 一边亲一边摸免费视频| 国产乱人偷精品视频| 亚洲乱码一区二区免费版| 春色校园在线视频观看| 在线天堂最新版资源| 最近中文字幕高清免费大全6| 久久这里有精品视频免费| 美女国产视频在线观看| 日本三级黄在线观看| 黄片无遮挡物在线观看| 麻豆精品久久久久久蜜桃| 国产av码专区亚洲av| 久久6这里有精品| 少妇的逼水好多| 国产成年人精品一区二区| 免费看av在线观看网站| 2022亚洲国产成人精品| 午夜精品在线福利| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区免费毛片| 国产精品一区二区三区四区免费观看| 视频中文字幕在线观看| 看片在线看免费视频| 高清日韩中文字幕在线| 亚洲精品一区蜜桃| 热99在线观看视频| 午夜精品国产一区二区电影 | 国产精品野战在线观看| 免费电影在线观看免费观看| 亚洲精品乱码久久久v下载方式| 精品人妻视频免费看| 国产乱来视频区| 日韩人妻高清精品专区| 一本一本综合久久| 麻豆久久精品国产亚洲av| 成人无遮挡网站| 日本熟妇午夜| 丰满乱子伦码专区| 成人一区二区视频在线观看| 国产亚洲91精品色在线| 99热精品在线国产| 国产欧美另类精品又又久久亚洲欧美| 国产成人a区在线观看| 国产成人91sexporn| 内地一区二区视频在线| 国产精品一区二区性色av| 久久久久久伊人网av| 国产一级毛片七仙女欲春2| 天天躁夜夜躁狠狠久久av| 丰满乱子伦码专区| 国产人妻一区二区三区在| 欧美一级a爱片免费观看看| 国产精品一区二区性色av| 男女国产视频网站| 国产视频内射| 成人亚洲欧美一区二区av| 床上黄色一级片| 如何舔出高潮| av在线蜜桃| 一夜夜www| 免费观看精品视频网站| 日韩中字成人| 1024手机看黄色片| 精品少妇黑人巨大在线播放 | 最近2019中文字幕mv第一页| 看十八女毛片水多多多| 我的女老师完整版在线观看| 韩国av在线不卡| 午夜福利在线观看免费完整高清在| 午夜福利在线观看吧| .国产精品久久| 两个人视频免费观看高清| 欧美激情国产日韩精品一区| 一个人免费在线观看电影| 高清视频免费观看一区二区 | 少妇熟女欧美另类| 国产免费一级a男人的天堂| 1024手机看黄色片| 永久免费av网站大全| 中文乱码字字幕精品一区二区三区 | 91久久精品电影网| 热99re8久久精品国产| 日本免费一区二区三区高清不卡| 午夜激情欧美在线| 亚洲精品,欧美精品| 成年免费大片在线观看| 国产午夜精品久久久久久一区二区三区| 久久精品国产亚洲网站| 成人三级黄色视频| 久久精品国产99精品国产亚洲性色| 久久99蜜桃精品久久| 欧美丝袜亚洲另类| 2022亚洲国产成人精品| 少妇熟女欧美另类| 久久国产乱子免费精品| 九九爱精品视频在线观看| 久久99热这里只有精品18| 99热这里只有精品一区| 菩萨蛮人人尽说江南好唐韦庄 | 国产熟女欧美一区二区| 亚洲国产最新在线播放| 久久久久国产网址| 中国国产av一级| 亚洲av中文字字幕乱码综合| 乱人视频在线观看| 青春草视频在线免费观看| 国产成人福利小说| 97超碰精品成人国产| 亚洲av电影在线观看一区二区三区 | 美女内射精品一级片tv| 人人妻人人澡人人爽人人夜夜 | 最近手机中文字幕大全| 国产成人精品久久久久久| 久久久久网色| 日产精品乱码卡一卡2卡三| 干丝袜人妻中文字幕| 亚洲成人中文字幕在线播放| 亚洲美女搞黄在线观看| 国产欧美另类精品又又久久亚洲欧美| 天堂av国产一区二区熟女人妻| 丰满乱子伦码专区| 国产探花极品一区二区| 午夜福利在线观看吧| 精品熟女少妇av免费看| 国产成人aa在线观看| 国产三级中文精品| 国产一级毛片在线| 亚洲av男天堂| av专区在线播放| 久久久久久大精品| 人妻制服诱惑在线中文字幕| 亚洲va在线va天堂va国产| 成人毛片60女人毛片免费| 欧美xxxx性猛交bbbb| 2022亚洲国产成人精品| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| av黄色大香蕉| 国产一区二区在线av高清观看| 国产亚洲精品av在线| 国产高清有码在线观看视频| 欧美日本亚洲视频在线播放| 国产乱来视频区| 精品久久久久久久久亚洲| 欧美高清性xxxxhd video| 99热网站在线观看| 国产淫语在线视频| 久久久精品欧美日韩精品| 不卡视频在线观看欧美| 一本久久精品| 久久久久久九九精品二区国产| 成人性生交大片免费视频hd| 亚洲欧美精品专区久久| 成年av动漫网址| av免费在线看不卡| 日韩一区二区视频免费看| 青青草视频在线视频观看| 国产精品不卡视频一区二区| 亚洲精品国产成人久久av| 国产精品野战在线观看| 国产成人精品一,二区| 国产 一区精品| 国产精华一区二区三区| 欧美不卡视频在线免费观看| 亚洲内射少妇av| av在线老鸭窝| 美女黄网站色视频| 男女那种视频在线观看| 日韩大片免费观看网站 | 一个人看视频在线观看www免费| 国产精品女同一区二区软件| 亚洲精华国产精华液的使用体验| 久久久久久国产a免费观看| 青春草视频在线免费观看| 97超视频在线观看视频| 99久久成人亚洲精品观看| 久久婷婷人人爽人人干人人爱| 色尼玛亚洲综合影院| 午夜a级毛片| 久热久热在线精品观看| 边亲边吃奶的免费视频| 精品一区二区三区视频在线| 亚洲国产高清在线一区二区三| 亚洲国产最新在线播放| 亚洲av二区三区四区| 欧美成人a在线观看| 最近中文字幕2019免费版| 日本一二三区视频观看| 日本免费在线观看一区| 丝袜美腿在线中文| 免费看光身美女| 免费看美女性在线毛片视频| 97超碰精品成人国产| av在线天堂中文字幕| 特级一级黄色大片| av在线播放精品| 国产精品乱码一区二三区的特点| 免费大片18禁| 日日摸夜夜添夜夜添av毛片| 搡老妇女老女人老熟妇| 免费无遮挡裸体视频| 少妇猛男粗大的猛烈进出视频 | 中文乱码字字幕精品一区二区三区 | 中文字幕av在线有码专区| 成人特级av手机在线观看| 天堂av国产一区二区熟女人妻| 久久久久久大精品| 国产午夜精品论理片| 乱人视频在线观看| 日韩欧美在线乱码| 激情 狠狠 欧美| 国产精品麻豆人妻色哟哟久久 | 七月丁香在线播放| 国产精品永久免费网站| 搡女人真爽免费视频火全软件| 国产精品久久久久久久久免| 国产一区二区在线观看日韩| 国产av一区在线观看免费| 久久久久久伊人网av| 黄色欧美视频在线观看| 国产黄a三级三级三级人| 国产在线男女| 成人漫画全彩无遮挡| 在线播放无遮挡| 视频中文字幕在线观看| 99久久九九国产精品国产免费| 天天躁日日操中文字幕| av福利片在线观看| 一区二区三区免费毛片| 全区人妻精品视频| 国产精品嫩草影院av在线观看| 国产爱豆传媒在线观看| 久久人人爽人人片av| 亚洲自拍偷在线| 国产91av在线免费观看| 婷婷六月久久综合丁香| 真实男女啪啪啪动态图| 国产一区二区亚洲精品在线观看| 97超碰精品成人国产| 男女国产视频网站| 中国美白少妇内射xxxbb| 99九九线精品视频在线观看视频| 人妻制服诱惑在线中文字幕| 男女视频在线观看网站免费| 网址你懂的国产日韩在线| videossex国产| 白带黄色成豆腐渣| 一区二区三区免费毛片| 国产熟女欧美一区二区| 成人av在线播放网站| 只有这里有精品99| 超碰97精品在线观看| 黄片wwwwww| 又爽又黄a免费视频| 亚州av有码| 亚洲在线自拍视频| 国产成年人精品一区二区| 亚洲三级黄色毛片| a级毛片免费高清观看在线播放| 高清视频免费观看一区二区 | 欧美性猛交╳xxx乱大交人| 久久久久久久久中文| 美女xxoo啪啪120秒动态图| 精品少妇黑人巨大在线播放 | 欧美+日韩+精品| 综合色丁香网| 在线a可以看的网站| 少妇猛男粗大的猛烈进出视频 | 又爽又黄a免费视频| 九色成人免费人妻av| 在线播放国产精品三级| 美女国产视频在线观看| av在线老鸭窝| 久久精品综合一区二区三区| 最近最新中文字幕免费大全7| 日本三级黄在线观看| 精品一区二区免费观看| 国产日韩欧美在线精品| 九色成人免费人妻av| 国产伦理片在线播放av一区| 国产免费视频播放在线视频 | 国产乱人视频| 久久精品国产亚洲网站| 精品人妻视频免费看| 亚洲国产精品久久男人天堂| 少妇丰满av| 久久久久久久亚洲中文字幕| 综合色av麻豆| 国产毛片a区久久久久| 男人和女人高潮做爰伦理| av在线观看视频网站免费| 亚洲av电影不卡..在线观看| 免费播放大片免费观看视频在线观看 | 又爽又黄无遮挡网站| 99久久精品热视频| 天堂网av新在线| 久久综合国产亚洲精品| 免费av不卡在线播放| 日韩亚洲欧美综合| 自拍偷自拍亚洲精品老妇|