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

    喀斯特區(qū)天峨槭(Acer wangchii)樹干液流特征及其與環(huán)境因子的相關(guān)分析

    2018-01-03 08:29:47楊文斌崔迎春趙文君舒德遠(yuǎn)侯貽菊丁訪軍
    生態(tài)學(xué)報(bào) 2017年22期
    關(guān)鍵詞:陰天液流晴天

    吳 鵬, 楊文斌, 崔迎春, 趙文君, 舒德遠(yuǎn), 侯貽菊,丁訪軍,*

    1 中國(guó)林業(yè)科學(xué)研究院荒漠化研究所,北京 100091 2 貴州省林業(yè)科學(xué)研究院,貴陽(yáng) 550005

    喀斯特區(qū)天峨槭(Acerwangchii)樹干液流特征及其與環(huán)境因子的相關(guān)分析

    吳 鵬1, 2, 楊文斌1, 崔迎春2, 趙文君2, 舒德遠(yuǎn)2, 侯貽菊2,丁訪軍2,*

    1 中國(guó)林業(yè)科學(xué)研究院荒漠化研究所,北京 100091 2 貴州省林業(yè)科學(xué)研究院,貴陽(yáng) 550005

    采用TDP(Thermal Dissipation Probe)熱擴(kuò)散探針法,于2012年2月—2014年3月,對(duì)喀斯特地區(qū)天峨槭(Acerwangchii)樹干液流速率進(jìn)行了長(zhǎng)期連續(xù)的監(jiān)測(cè),并同步監(jiān)測(cè)了林分空氣溫度(Ta)、相對(duì)濕度(RH)、太陽(yáng)輻射(Solar)、風(fēng)速(WS)、降雨量(Rainfall)、土壤含水量(SWC)等環(huán)境因子,探討了不同時(shí)間尺度下天峨槭樹干液流特征及其與環(huán)境因子的關(guān)系。結(jié)果表明:1)不同天氣條件下的樹干液流為晴天>陰天>雨天,且均呈明顯的“晝高夜低”變化規(guī)律;2)天峨槭樹干的日平均液流量為5.08 kg/d,不同季節(jié)表現(xiàn)為夏季((8.38±5.32) kg/d)>秋季((5.16±3.99) kg/d)>春季((4.86±3.77) kg/d)>冬季((1.94 ±1.40) kg/d);3)月平均液流量為153.64 kg,年平均液流量達(dá)1838.40 kg;4)小時(shí)尺度下,影響晴天、陰天、雨天全天樹干液流的主要環(huán)境因子都是Solar,但各環(huán)境因子對(duì)樹干液流的影響程度又因晝夜、季節(jié)及降雨量的不同而存在差異;5)從年度范圍來看:小時(shí)尺度下,Solar、Ta、RH、SWC10cm和WS可以共同解釋樹干液流的63.50%;日尺度下,Solar、Ta、SWC10cm和Rainfall可共同解釋其68.50%;月尺度下,單個(gè)因子Ta就能解釋其74.80%;且隨著時(shí)間尺度的縮小,各環(huán)境因子入選回歸方程的個(gè)數(shù)有增加的趨勢(shì),而對(duì)樹干液流的解釋程度(R2)則有降低的趨勢(shì)。對(duì)比其他地區(qū)研究結(jié)果,其環(huán)境因子對(duì)樹干液流的影響差異都很大;但總的來說,無論在何種時(shí)間尺度上,Solar(或光合有效輻射(PAR))和Ta基本上都是影響樹干液流的主要環(huán)境因子,且各環(huán)境因子對(duì)雨天樹干液流的解釋程度都不高,本研究亦如此;因此,以環(huán)境因子對(duì)雨天樹干液流進(jìn)行預(yù)測(cè)的時(shí)候可能會(huì)存在誤差,應(yīng)特別注意土壤含水量對(duì)其的影響。

    喀斯特;天峨槭;樹干液流;環(huán)境因子;不同時(shí)間尺度

    植物蒸騰是土壤-植物-大氣連續(xù)體水熱傳輸過程中一個(gè)極為重要的環(huán)節(jié)[1],也是植物本身生長(zhǎng)過程中極為重要的生理生態(tài)過程之一,植物體內(nèi)超過90%的樹干液流(忽略了樹木新陳代謝過程中對(duì)極小部分水分的利用)都用以蒸騰[2];而水分是植物生長(zhǎng)的主要限制性生態(tài)因子之一[3],樹木對(duì)水分利用的主要?jiǎng)恿κ钦趄v需求,只有在蒸騰拉力的作用下[4- 5],這些水分才能被用于植物不同的組織器官,以保證其健康生長(zhǎng),因此可以用樹干液流來表征植物的蒸騰量[6];通過對(duì)植物樹干液流變化規(guī)律及其與相關(guān)環(huán)境因子關(guān)系的研究,對(duì)于理解植物水分需耗、估算不同時(shí)間尺度下植被蒸散具有重要意義[7- 8]。有研究表明,植物的蒸騰耗水主要受植物個(gè)體的生物學(xué)結(jié)構(gòu)(樹干直徑、冠幅、邊材面積等)[9- 10]、土壤水分[6,10-13]和氣象因子[1- 2,6,8,10,12-14]等的影響;生物學(xué)結(jié)構(gòu)主要決定了樹干液流的潛力,土壤供水決定了樹干液流的總體水平,而氣象因子則決定了樹干液流的瞬間變動(dòng)[10,15]。

    目前,已有的研究多是針對(duì)常態(tài)地貌上植物的樹干液流來進(jìn)行,而對(duì)喀斯特地貌上的研究甚少,僅張中峰[16]、黃玉清[17]、劉延惠[18]等人分別對(duì)巖溶地區(qū)的青岡櫟(Cyclobalanopsisglauca)和短萼海桐(Pittosporumbrevicalyx)樹干液流特征及其與環(huán)境因子的關(guān)系進(jìn)行了研究??λ固氐孛采系纳种脖?其生境與相同氣候條件下常態(tài)地貌上的明顯不同[19],特殊的地貌構(gòu)造使其生態(tài)水文過程更加復(fù)雜,對(duì)其樹木蒸騰、生理生態(tài)過程以及對(duì)土壤間歇性干旱的適應(yīng)機(jī)制尚缺乏足夠認(rèn)識(shí)。為此,本文以喀斯特森林頂極群落主要優(yōu)勢(shì)樹種天峨槭(Acerwangchii)為研究對(duì)象,利用熱擴(kuò)散技術(shù)對(duì)其樹干液流進(jìn)行連續(xù)不間斷的監(jiān)測(cè),并結(jié)合自動(dòng)氣象站對(duì)環(huán)境氣象因子的同步測(cè)定數(shù)據(jù),旨在探討:1)不同時(shí)間尺度下天峨槭樹干液流特征及其與環(huán)境因子的內(nèi)在關(guān)系是否存在差異;2)不同時(shí)間尺度下環(huán)境因子對(duì)樹干液流的擬合對(duì)簡(jiǎn)化分析樹干液流是否具有意義;3)喀斯特地貌上樹干液流特征及其對(duì)環(huán)境因子的響應(yīng)與其他地區(qū)是否存在差異。以期能為喀斯特森林生態(tài)需水估算以及退化生態(tài)系統(tǒng)植被恢復(fù)與重建過程中樹種的篩選提供基礎(chǔ)數(shù)據(jù)。

    1 研究區(qū)概況

    研究區(qū)位于貴州喀斯特森林生態(tài)系統(tǒng)定位觀測(cè)研究站,地處茂蘭國(guó)家級(jí)自然保護(hù)區(qū)內(nèi)(107°52.2′—108°5.7′ E,25°9.3′—25°20.8′ N)。保護(hù)區(qū)最低海拔430.0 m,最高1078.6 m,平均在550.0—850.0 m之間;屬中亞熱帶季風(fēng)濕潤(rùn)氣候,具有春秋溫暖、冬無嚴(yán)寒、夏無酷暑、雨量充沛的特點(diǎn);年平均氣溫15.3℃,≥10℃積溫4598.6℃;年平均降雨量1752.5 mm,集中分布于4—10月;年平均相對(duì)濕度83.0%;無霜期315 d;全年日照時(shí)數(shù)1272.8 h,日照百分率29.0%;太陽(yáng)年輻射總量為63289.80 kW/m2。保護(hù)區(qū)內(nèi)主要是由純質(zhì)石灰?guī)r及白云巖構(gòu)成的喀斯特地貌,裸巖率在80%以上,屬于典型的喀斯特生境;土壤稀少,土層淺薄(20 cm左右)且不連片,多存于巖石縫隙之中,土壤富鈣和富鹽基化,有機(jī)質(zhì)含量高。該保護(hù)區(qū)是目前世界上同緯度地區(qū)殘存下來僅有的、原生性強(qiáng)且相對(duì)穩(wěn)定的喀斯特森林生態(tài)系統(tǒng),同時(shí)也是喀斯特區(qū)原生性森林(常綠落葉闊葉混交林)分布面積最大的區(qū)域[20]。

    2 研究方法

    2.1 樣地調(diào)查與樣木選擇

    2012年,對(duì)站區(qū)原有的4個(gè)固定樣地(30 m×30 m)進(jìn)行了調(diào)查[21],樣地基本概況(表1)。

    表1 樣地基本概況

    根據(jù)調(diào)查結(jié)果,在樣地內(nèi)選取徑級(jí)適中(接近樣地內(nèi)的平均胸徑)且生長(zhǎng)良好、樹干通直、無病蟲害、冠幅相似的3株天鵝槭作為監(jiān)測(cè)樣樹(表2)。

    表2 天峨槭樣樹基本特征

    2.2 樹干液流測(cè)定

    利用熱擴(kuò)散液流探針(TDP- 30,Dynamax Inc., Houston, USA)產(chǎn)生的溫差來計(jì)算天峨槭樹干液流速率,TDP探針纜線與CR1000數(shù)據(jù)采集器連接,以10 min時(shí)間間隔進(jìn)行數(shù)據(jù)采集;采集時(shí)間:2012年2月25日—2014年3月22日;根據(jù)Granier經(jīng)驗(yàn)公式[22]來?yè)Q算樹干邊材瞬時(shí)液流速率(Js):

    Js=0.0119×[(ΔTmax-ΔT)/ΔT]1.231

    (1)

    式中,Js是液流密度(瞬時(shí)液流速率)(g H2O/(m2/s)),即單位時(shí)間通過單位邊材面積的液流量[23];ΔTmax為上下探針間的最大晝夜溫差(℃),ΔT瞬時(shí)溫差(℃)。

    Fs=Js×SA×3600

    (2)

    式中,Fs為液流通量(每小時(shí)的液流速率)(g H2O/h);SA為胸徑處邊材面積;3600是時(shí)間換算系數(shù)[23]。

    邊材面積(SA)的確定:在樣地外選擇12株不同胸徑的天峨槭,在樹干胸徑處用生長(zhǎng)錐鉆取樣芯。將樣芯帶回實(shí)驗(yàn)后,用游標(biāo)卡尺分別測(cè)量樣芯的心材寬、邊材寬、半徑和皮厚等值,計(jì)算樣木邊材面積和胸徑,建立了天峨槭邊材面積與胸徑的回歸關(guān)系式:

    SA=0.07426×D1.8077(R2=0.7898)

    (3)

    式中,D為樣木的胸徑,R為相關(guān)系數(shù)。

    日間或夜間的樹干液流對(duì)日總液流量的貢獻(xiàn)率,計(jì)算公式如下[12]:

    (4)

    式中,CR為日間或夜間的樹干液流對(duì)日累積量的貢獻(xiàn)率(%);Fsi為i時(shí)間的液流速率(g H2O/h);i為時(shí)間段(1、2、…、t);t為日間(以太陽(yáng)輻射大于零定義為“日間”)或夜間(太陽(yáng)輻射等于零為“夜間”)的時(shí)長(zhǎng);24為全天時(shí)長(zhǎng)。

    基于小時(shí)尺度下的液流速率(Fs)可以計(jì)算日尺度下的液流量(Fd)和月尺度下的液流量(Fm):

    (5)

    (6)

    式中,Fd為日液流量(kg H2O/d);0.001是單位換算系數(shù);Fm為月液流量(kg H2O);n為每月的天數(shù);Fdj為該月第j天的日液流量。

    2.3 環(huán)境因子測(cè)定

    環(huán)境因子由安裝在觀測(cè)塔上的Campbell自動(dòng)氣象站數(shù)據(jù)采集器同步采集,主要包括:太陽(yáng)輻射(Solar)、氣溫(Ta)、相對(duì)濕度(RH)、風(fēng)速(WS)、降雨量(Rainfall)、不同層次的土壤含水量(SWC10cm、SWC20cm、SWC40cm)等。

    2.4 數(shù)據(jù)處理方法

    1)異常值的處理:首先建立液流速率正常值與空氣溫度和太陽(yáng)輻射的擬合模型,然后利用異常值(缺失值)所對(duì)應(yīng)的空氣溫度和太陽(yáng)輻射,反過來計(jì)算其液流速率,用該液流速率來代替其異常值(缺失值);異常值的判別方法參考文獻(xiàn)[4,24]。

    2)小時(shí)尺度下的時(shí)滯效應(yīng):采用錯(cuò)位對(duì)比法來確定[25]。

    3)數(shù)據(jù)分析:采用Excel 2010進(jìn)行繪圖,SPSS 22.0進(jìn)行方差分析(Repeated-measures ANOVA)、多重比較(LSD最小顯著差法)、相關(guān)性分析、逐步回歸分析及通徑分析等;通徑系數(shù)計(jì)算方法參見文獻(xiàn)[26]。

    3 結(jié)果與分析

    3.1 小時(shí)尺度下樹干液流特征及其與環(huán)境因子的相關(guān)性

    3.1.1典型晴天和陰天樹干液流特征

    選取典型晴天和陰天天峨槭樹干液流的平均值,按月、分不同季節(jié)(春季:3—5月,夏季:6—8月,秋季:9—11月,冬季:12—次年2月)分析樹干液流的日變化(圖1)。小時(shí)尺度下,樹干液流的日變化表現(xiàn)為晴天(388.60g/h)>陰天(216.30g/h);不同季節(jié)表現(xiàn)為夏季(686.35g/h晴和359.06g/h陰)>春季(393.68g/h晴和219.23g/h陰)>秋季(372.19g/h晴和198.79g/h陰)>冬季(102.19g/h晴和88.00g/h陰),且不同天氣條件下冬季與其他季節(jié)差異均顯著(P<0.05),夏季與秋季差異均顯著(P<0.05),而其他季節(jié)間差異則均不顯著(P>0.05);相同天氣條件下,同一季節(jié)不同月份樹干液流的日變化亦存在差異。從日變化的曲線波動(dòng)來看,除冬季以外,晴天樹干液流的日變化多呈平滑的“單峰”曲線;陰天則多為“單峰”“雙峰”和“多峰”曲線并存。

    圖1 不同季節(jié)、不同天氣條件下(典型晴天和陰天)天峨槭樹干液流的日變化Fig.1 Diurnal variation of stems sap flow in Acer wangchii under different weather conditions of different seasons

    由表3可知,不同季節(jié)晴天樹干液流從“‘啟動(dòng)’→達(dá)到‘峰值’→回落到‘低值’”的整個(gè)過程較陰天均呈現(xiàn)出:?jiǎn)?dòng)時(shí)間較早、到達(dá)峰值和回落到低值時(shí)間較晚、整個(gè)過程持續(xù)時(shí)間較長(zhǎng)的變化規(guī)律;由“夏季→冬季→夏季”,不同天氣條件下樹干液流的啟動(dòng)和到達(dá)峰值的時(shí)間均表現(xiàn)為先“推后”后“提前”,回落到低值的時(shí)間均為先“提前”后“推后”,整個(gè)過程持續(xù)時(shí)間均為先“縮短”后“延長(zhǎng)”。

    3.1.2典型雨天樹干液流特征

    按照降雨量的大小(中雨:10.0—24.9mm、大雨:25.0—49.9mm),選取連續(xù)降雨2—3d的數(shù)據(jù)作為典型雨天樹干液流的日變化(圖2)。

    雨天樹干液流的日變化多呈“多峰”或不規(guī)則上下波動(dòng),平均值為64.27g/h,且中雨天(70.52g/h)要高于大雨天(58.03g/h);與當(dāng)年同月晴天和陰天樹干液流的日變化相比,雨天的要顯著小于晴天和陰天(P<0.05)。

    表3 不同季節(jié)、不同天氣條件下樹干液流日變化持續(xù)過程時(shí)間表

    圖2 典型雨天天峨槭樹干液流和降雨量的日變化Fig.2 Diurnal variation of rainfall and stems sap flow in Acer wangchiiof rainy day

    3.1.3小時(shí)尺度下樹干液流的晝夜差異

    不同天氣條件下的樹干液流均呈明顯的“晝高夜低”變化規(guī)律,且差異均顯著(P<0.05);晴天、陰天、雨天,日間樹干液流的平均值分別是夜間的12.37倍、6.59倍和2.49倍。另外,不同天氣條件下日間和夜間樹干液流對(duì)日總液流量的貢獻(xiàn)率也不同;晴天,晝、夜貢獻(xiàn)率的平均值分別為90.53%和9.47%,陰天分別為84.08%和15.92%,雨天分別為77.26%和22.74%,夜間貢獻(xiàn)率以雨天為最高,陰天次之,晴天最低。晴天和陰天不同季節(jié)夜間的貢獻(xiàn)率則均以冬季為最高,春季次之,夏季最低(表4)。

    表4 小時(shí)尺度下樹干液流的晝夜差異/(g/h)

    同列之間凡標(biāo)有不同小寫字母的表示差異顯著(P<0.05);x為對(duì)日總液流量的貢獻(xiàn)率

    3.1.4小時(shí)尺度下樹干液流與環(huán)境因子的相關(guān)性

    綜合圖1、圖2可知,不同天氣條件下樹干液流均受環(huán)境因子的影響,但其影響方式、作用大小并不清楚。為此,我們進(jìn)一步分析了各環(huán)境因子(Solar、Ta、RH、WS、SWC10cm、SWC20cm、SWC40cm和Rainfall)與樹干液流的相關(guān)系數(shù)、直接及間接通徑系數(shù)(表5)。

    (1)樹干液流晝夜變化與環(huán)境因子的相關(guān)性

    由表5可知(根據(jù)相關(guān)系數(shù)、直接和間接通徑系數(shù)絕對(duì)值的大小進(jìn)行判斷),不同天氣條件下,Solar都是影響全天樹干液流的主要環(huán)境因子;Solar、Ta是影響晴天和陰天日間的主要環(huán)境因子,而Solar、SWC10cm是影響雨天日間的主要環(huán)境因子;Ta、RH是影響晴天和陰天夜間的主要環(huán)境因子,而SWC10cm是影響雨天夜間的主要環(huán)境因子。另外,不同天氣條件下樹干液流對(duì)土壤水分的利用和響應(yīng)也存在差異,晴天和陰天受日間較強(qiáng)的Solar和Ta等環(huán)境因子的影響,致使土壤表層含水量減小,較深層的土壤水分就成了夜間樹干液流主要的水分來源,所以SWC20cm對(duì)晴天和陰天夜間的樹干液流有較明顯的促進(jìn)作用;雨天全天及晝夜,由于日間太陽(yáng)輻射小、溫度較低,樹木蒸騰潛力小,且水分供應(yīng)充足,極高的土壤表層含水量(SWC10cm)反而對(duì)樹干液流產(chǎn)生了抑制作用。

    表5 小時(shí)尺度下樹干液流晝夜變化與環(huán)境因子的通徑分析結(jié)果

    **P<0.01,*P<0.05;Ta、RH、WS、Solar、SWC分別代表空氣溫度Air temperature、相對(duì)濕度Relative humidity、風(fēng)速Wind speed、太陽(yáng)輻射Solar radiation和土壤含水量Soil water content

    表6 小時(shí)尺度下樹干液流晝夜變化與環(huán)境因子的逐步回歸分析結(jié)果

    (2)樹干液流季節(jié)變化與環(huán)境因子的相關(guān)性

    不同天氣條件下各環(huán)境因子對(duì)樹干液流的響應(yīng)及與其擬合的回歸方程因季節(jié)的不同也存在差異(表7,因篇幅所限,此表省略了間接通徑系數(shù)的計(jì)算)。晴天,影響春季、秋季和冬季樹干液流的環(huán)境因子主要是Solar和Ta,直接通徑系數(shù)都在0.5左右;影響夏季的主要是Solar(0.602),Ta對(duì)其直接影響較小(0.149)。陰天,不同季節(jié)樹干液流受環(huán)境因子的影響差異較大,春季主要受Solar、Ta和RH的影響,夏季的是RH,秋季的是Solar,冬季則受Solar和Ta的影響較大。大雨天,主要受Solar、SWC10cm和Ta的影響,相較于Ta對(duì)晴天、陰天樹干液流明顯的促進(jìn)作用,Ta對(duì)大雨天樹干液流抑制作用明顯;中雨天,影響樹干液流的主要是Solar。不同天氣條件,RH對(duì)樹干液流多是抑制作用;而WS對(duì)樹干液流的作用因天氣和季節(jié)的不同而存在差異。逐步回歸分析結(jié)果表明(表7):入選方程的各環(huán)境因子與樹干液流的相關(guān)系數(shù)以及各擬合方程的回歸系數(shù)(R2)均達(dá)到極顯著水平(P<0.01);晴天,不同季節(jié)各環(huán)境因子對(duì)樹干液流的解釋程度在82.00%—93.40%之間;陰天,各環(huán)境因子對(duì)樹干液流的解釋程度在81.60%—93.90%之間;雨天,各環(huán)境因子對(duì)樹干液流的解釋程度較低,大雨天和中雨天分別為39.80%和47.70%。

    3.2 日尺度下樹干液流特征及其與環(huán)境因子的相關(guān)性

    由圖3可知,天峨槭的日液流量在0.01—24.00 kg/d之間,平均值為5.08 kg/d。不同季節(jié)日平均液流量分別為:春季((4.86±3.77) kg/d)、夏季((8.38±5.32) kg/d)、秋季((5.16±3.99) kg/d)和冬季的((1.94 ±1.40) kg/d),除春季和秋季間差異不顯著以外(P>0.05),其余各季節(jié)間均達(dá)到極顯著差異(P<0.01);冬季的日液流量最低,低溫(7.56℃)和少的降雨量(104.75 mm)可能是其主要限制因子;夏季的日液流量最高,主要得益于高溫(23.33℃)和強(qiáng)太陽(yáng)輻射(134.15 W/s)致使其蒸騰潛力很大,且水分供應(yīng)充足(降雨量1121.00 mm),致使其液流量很高;秋季的日液流量顯著低于夏季,可能是由于該季節(jié)降雨量偏少(412.70 mm)、土壤含水量過低等原因造成的;春季的降雨量雖然充足(1002.40 mm),但其空氣溫度(16.48℃)和太陽(yáng)輻射(110.70 W/s)較夏季要低,致使其蒸騰潛力不高,使得其日液流量也顯著低于夏季。

    逐步回歸分析結(jié)果表明(表8):入選方程的各環(huán)境因子對(duì)樹干日液流量的直接作用以及方程的回歸系數(shù)均達(dá)到顯著(P<0.05)或極顯著水平(P<0.01),說明所擬合方程均能很好的揭示不同季節(jié)及全年樹干日液流量與環(huán)境因子的變化規(guī)律;除冬季各環(huán)境因子對(duì)日液流量的解釋程度(33.90%)較低以外,其余各季節(jié)的解釋程度均超過了70.00%。從全年來看,各環(huán)境因子對(duì)其解釋程度也高達(dá)68.50%,而Solar和Ta是影響全年日液流量的主要環(huán)境因子,可解釋日液流量的65.80%,占總解釋程度的96.06%。

    表7 小時(shí)尺度下樹干液流季節(jié)變化與環(huán)境因子的逐步回歸分析結(jié)果

    注:**P<0.01,*P<0.05

    圖3 天峨槭樹干液流與空氣溫度、相對(duì)濕度、太陽(yáng)輻射、土壤含水量、風(fēng)速及降雨量的日均值變化Fig.3 Changes of daily of stems sap flow in Acer wangchii Tair, RH, Solar radiation, SWC, WS and Rainfall

    尺度Scale季節(jié)Season環(huán)境因子Environmentalfactors相關(guān)系數(shù)Correlationcoefficient直接通徑系數(shù)Directpathcoefficient回歸方程Regressionequation樣本數(shù)NR2FP日尺度春季Solar0.768**0.620**Fd=3.043+0.029Solar+0.273Ta-2060.723175.8240.000DailyscaleTa0.566**0.363**12.291SWC10cmSWC10cm-0.217**-0.192**夏季RH-0.818**-0.533**Fd=25.188-0.366RH+0.029So-lar+1840.794172.0820.000Solar0.810**0.359**0.572Ta-1.605WSTa0.368**0.159**WS-0.081-0.105**秋季Solar0.816**0.742**Fd=0.457+0.043Solar+0.176Ta-1820.724155.9020.000Ta0.477**0.172**6.298SWC10cmSWC10cm-0.272**-0.136**冬季SWC20cm0.452**0.400**Fd=-6.573+12.852SWC20cm-1850.33923.0390.000Rainfall-0.105-0.220**0.093Rainfall+0.009Solar+0.031RHSolar0.175**0.335**RH0.225**0.339**年度Solar0.753**0.581**Fd=1.119+0.036Solar+0.246Ta-7570.685407.9760.000Ta0.599**0.374**18.133SWC10cm-0.025RainfalSWC10cm-0.114**-0.141**Rainfall-0.139**-0.049*月尺度Monthlyscale年度Ta0.865**0.865**Fm=-58.540+13.207Ta250.74868.1330.000

    3.3 月尺度下樹干液流特征及其與環(huán)境因子的相關(guān)性

    由圖4可知,在月尺度下(2012年2月份只有5天數(shù)據(jù),不能代表該月的情況,故舍去;2014年3月份有22天的數(shù)據(jù),基本上可以代表該月的情況,故對(duì)其保留,為便于比較,該月的月液流量=該月日液流量的平均值×該月的天數(shù)),天峨槭的月平均液流量為153.64 kg,其中以2012年8月份的417.61 kg最高,2014年2月份的37.48 kg最低。不同年際間同一月份或同一季節(jié)的液流量因環(huán)境因子的變化也存在差異;以3月份為例,2012年、2013年和2014年的月液流量分別為86.68、131.19 kg和37.68 kg,2013年分別是2012年和2014年的1.51倍和3.48倍;2012年春、夏、秋、冬各季的樹干液流總量分別為:502.36、959.59、490.56、218.57 kg,2013年的則分別為472.36、582.46、448.64、128.78 kg,2012年夏季和冬季的樹干液流總量分別為2013年同季的1.65倍和1.70倍。2012年3月1日—2014年3月22日間,天峨槭樹干液流總量為3830.05 kg,年平均液流量可達(dá)1838.40 kg。逐步回歸分析結(jié)果顯示(表8),在月尺度下,Ta是影響其液流量的主要環(huán)境因子,與液流量呈極顯著的正相關(guān)關(guān)系(P<0.01),可解釋液流量的74.80%,Ta每升高1℃,液流量可增加13.21 kg。

    4 結(jié)論與討論

    小時(shí)尺度下,天峨槭樹干液流的日變化表現(xiàn)為晴天>陰天>雨天(圖1、圖2),這主要是由于晴天Solar強(qiáng)、Ta高、RH小,植物蒸騰量大,致使其液流也大,而陰雨天Solar弱、RH大,植物蒸騰量小,致使其液流也較小,另外,雨天極高的RH使葉片內(nèi)外的蒸汽壓梯度降低,促使葉片氣孔關(guān)閉,也制約了樹干液流的上升[27]。晴天樹干液流主要受Solar和Ta的影響,而這兩個(gè)環(huán)境因子在晴天多為“單峰”變化,且波動(dòng)較小,使得其樹干液流也多呈平滑“單峰”曲線變化;陰天樹干液流主要受Solar的影響,而Solar在陰天上下波動(dòng)較大,致使其樹干液流出現(xiàn)了“單峰”“雙峰”和“多峰”并存的現(xiàn)象;雨天除了各環(huán)境因子對(duì)其影響以外,還有“根壓”[27]以及降雨和光照的間斷性等對(duì)其的影響,致使其日變化呈“多峰”或不規(guī)則上下波動(dòng);此結(jié)果與張璇[12]、張中峰[16]、劉華[28]、凡超[29]等人的研究結(jié)果基本一致。有研究認(rèn)為,在中午前后,植物為了保存體內(nèi)水分,會(huì)短暫關(guān)閉或減小葉片氣孔開合程度,使其蒸騰量降低,導(dǎo)致樹干液流出現(xiàn)明顯或較明顯的“午休”現(xiàn)象[30- 31];而在本研究中,除了天氣原因?qū)е聵涓梢毫鞒识鄠€(gè)峰值或上下波動(dòng)以外,并未出現(xiàn)明顯的“午休”現(xiàn)象;此結(jié)果與張中峰等人[16]的研究結(jié)果一致;其原因可能是植物為了適應(yīng)喀斯特區(qū)高溫高濕的環(huán)境以及生境的高度異質(zhì)性而做出的水分調(diào)控策略,也可能是由于該區(qū)域太陽(yáng)輻射的峰值不足以使植物的葉片氣孔關(guān)閉而形成“光合午休”現(xiàn)象[12]。

    與大多數(shù)研究結(jié)果一致[1,10,12- 13,29,31- 32],不同天氣條件下天峨槭樹干液流的日變化均呈明顯的“晝高夜低”變化規(guī)律,晝、夜樹干液流因季節(jié)的不同及降雨量的大小而存在差異。有學(xué)者認(rèn)為,樹干的夜間液流是為了彌補(bǔ)日間蒸騰導(dǎo)致的水分虧缺而在樹體內(nèi)儲(chǔ)存水分[30,33],白天強(qiáng)烈的蒸騰作用,使樹體與土壤之間形成了較大水勢(shì)差,夜間在蒸騰停止后,在根壓作用下根系吸水方式由日間的被動(dòng)變?yōu)橐归g的主動(dòng),持續(xù)保持吸水能力,導(dǎo)致了夜間液流活動(dòng)的存在[30- 31];本研究中,晴天和陰天樹干夜間液流的平均值分別為53.32 g/h和56.02 g/h,均要高于雨天的35.79 g/h,其原因主要是因?yàn)橛晏熘参锏恼趄v量小,且由于降雨的存在,相對(duì)濕度高,水分供應(yīng)充足,植物體內(nèi)的水分虧缺不明顯,致使雨天的夜間液流較晴天和晴天要低;但是從夜間液流對(duì)日總液流量的貢獻(xiàn)率來看,由于雨天夜間液流速率與日間的差距較晴天、陰天夜間和日間的差距要小的多,導(dǎo)致其夜間液流的貢獻(xiàn)率反而最高。

    圖4 天峨槭樹干液流與空氣溫度、相對(duì)濕度、太陽(yáng)輻射、土壤含水量、風(fēng)速及降雨量的月均值變化Fig.4 Changes of monthly total of stems sap flow in Acer wangchii, Tair, RH, Solar radiation, SWC, WS and Rainfall

    天峨槭樹干的日、月平均液流量分別為5.08 kg/d和153.64 kg,年液流量達(dá)1838.40 kg,是湖南會(huì)同杉木(Cunninghamialanceolata)樹干液流的1.5倍以上[34],比華北落葉松(Larixprincipis-rupprechtii)生長(zhǎng)季的日液流量(3.17 kg/d)還要高[35];這可能與本研究所選樣樹邊材面積較大有關(guān),但相較其生長(zhǎng)的喀斯特生境(土壤稀少,土層淺薄,地表水缺乏),也可能是得益于其植物根系能從更深的巖溶內(nèi)部含水層中吸取水分的能力[17];但要低于程靜等[32]對(duì)鼎湖山針闊混交林主要優(yōu)勢(shì)樹種樹干液流的研究結(jié)果;這可能是因?yàn)闃涓梢毫鞯拇笮∈芏喾N因子的影響,立地條件、氣象因子、土壤溫度、水分以及土壤肥力等因子的不同,都能導(dǎo)致植物對(duì)水分的利用產(chǎn)生很大差異[32]。還有研究表明,同一樹種即使胸徑相同或相近,其液流差異也很大[34,36],這可能是由于樹種之間生物學(xué)結(jié)構(gòu)的不同而導(dǎo)致的。

    小時(shí)尺度下,影響晴天、陰天、雨天全天樹干液流的主要環(huán)境因子都是Solar,但各環(huán)境因子對(duì)樹干液流的影響又因晝夜差異、季節(jié)不同以及降雨量的大小而存在差異(表5、表7)。有研究認(rèn)為[37],晴天日間WS有利于植物的蒸騰作用,而夜間或陰雨天WS對(duì)液流影響不大,WS過大還會(huì)降低液流水平;在本研究中,WS對(duì)夜間的樹干液流都有明顯的促進(jìn)作用,反而對(duì)晴天日間的作用不明顯,對(duì)陰雨天日間的樹干液流還存在抑制作用,這個(gè)可能與不同的生境條件有關(guān)。

    土壤含水量(SWC)對(duì)樹干液流影響較為復(fù)雜。一般來說,在小尺度上,SWC對(duì)樹干液流無顯著影響,適宜在月、季節(jié)等大尺度上對(duì)其變化進(jìn)行分析[38];有學(xué)者認(rèn)為[12],SWC對(duì)雨天夜間的樹干液流影響顯著,它的變化是導(dǎo)致雨天夜間液流大于晴、陰天的主要因素;而王艷兵等人[33]通過對(duì)華北落葉松樹干液流的研究得出:SWC0—40cm是限制夜間樹干液流的主要因子,且月均土壤含水量對(duì)月蒸騰貢獻(xiàn)率影響顯著,但對(duì)雨天夜間樹干液流影響不顯著(P>0.05);還有研究認(rèn)為,土壤水分虧缺在一定程度上能限制生長(zhǎng)盛期刺槐(Robiniapseudoacacia)的蒸騰作用[10],但對(duì)展葉期刺槐樹干液流的影響不大[14];曹文強(qiáng)等[27]則認(rèn)為,SWC是解釋雨天遼東櫟(Quercusliaotungensis)樹干液流的主要環(huán)境因子,但對(duì)晴天和陰天樹干液流的影響不大;但張曉艷等[39]通過對(duì)民勤梭梭(Haloxylonammodendron(C. A. Mey)Bunge)樹干液流的研究,得出了與其正好相反的結(jié)論:土壤含水量與晴天樹干液流的相關(guān)性顯著,但與陰雨天相關(guān)性則不顯著;其原因可能與樹種的不同以及樹種根系對(duì)當(dāng)?shù)赝寥浪值倪m應(yīng)機(jī)制有關(guān)。本研究中,小時(shí)尺度下,SWC10cm是影響雨天全天和日間、夜間的主要環(huán)境因子,且對(duì)樹干液流有較強(qiáng)的抑制作用;SWC20cm對(duì)晴天和陰天夜間的樹干液流有明顯的促進(jìn)作用,但對(duì)晴天日間和全天的樹干液流存在抑制作用。在日尺度下,SWC10cm對(duì)春季、秋季和全年的樹干液流抑制作用較明顯;而SWC20cm是影響冬季樹干液流最主要的環(huán)境因子(表5、表7、表8);此結(jié)果與劉延惠等[18]對(duì)該地區(qū)短萼海桐樹干液流的研究亦存在一定差異,其原因可能是由于兩樹種在同一群落中所處的層次不同,造成了對(duì)水分利用機(jī)制的差異。

    從年度范圍來看,不同時(shí)間尺度下,各環(huán)境因子對(duì)樹干液流的影響也存在一定差異(表7、表8);小時(shí)尺度下,Solar、Ta、RH、SWC10cm和WS可以共同解釋樹干液流的63.50%,其中以Solar的直接影響最大,而RH和SWC10cm對(duì)其則有明顯的抑制作用;日尺度下,Solar、Ta、SWC10cm和Rainfall可以共同解釋樹干液流的68.50%,其中以Solar的直接影響最大,Ta次之,Rainfall最??;月尺度下,單個(gè)因子Ta就能解釋樹干液流74.80%;我們發(fā)現(xiàn),隨著時(shí)間尺度增大,與樹干液流顯著相關(guān)的環(huán)境因子個(gè)數(shù)有減少的趨勢(shì),而對(duì)其解釋程度則有增加的趨勢(shì),此結(jié)果與王文杰等[4]對(duì)興安落葉松(Larixgmelinii)樹干液流的研究結(jié)果一致,他們還認(rèn)為,影響樹干液流的主要環(huán)境因子隨著時(shí)間尺度的增大,有從地上氣象因子向地下土壤環(huán)境因子(土壤溫度、濕度等)轉(zhuǎn)變的趨勢(shì),本研究因沒有對(duì)土壤溫度的監(jiān)測(cè),故無法論證此結(jié)論。

    對(duì)比其他地區(qū)研究結(jié)果:王文杰等[4]認(rèn)為,在小時(shí)尺度上,光合有效輻射(PAR)、5 cm土壤溫度(Ts5 cm)、RH和Ta是影響興安落葉松樹干液流的主要環(huán)境因子;在日尺度上,Ts5cm、PAR和Ta是其主要影響因子;在月尺度上,Ts5cm和SWC5cm對(duì)其影響顯著。趙仲輝等[34]的研究則得出,在小時(shí)尺度上,Solar、Ta和Ts20cm是影響杉木樹干液流的主要環(huán)境因子;在日尺度上,Ta、SWC30cm、Solar和SWC20cm是其主要影響因子;在月尺度上,Ta和RH對(duì)其影響顯著。而通過對(duì)晉西黃土區(qū)遼東櫟樹干液流的研究,張建國(guó)等[40]認(rèn)為空氣水汽壓虧缺(VPD)對(duì)樹干液流能進(jìn)行很好的擬合,隋旭紅等[30]則認(rèn)為Ta和Solar是影響其樹干液流的主要環(huán)境因子,而陳寶強(qiáng)等[41]的研究則得出:不同月份影響樹干液流的環(huán)境因子也不同,5、6月份主要是VPD和PAR對(duì)其產(chǎn)生影響,而7、8月份則主要是VPD和Ta對(duì)其影響顯著;曹文強(qiáng)等[27]通過對(duì)太岳山遼東櫟樹干液流的研究則得出:晴天影響其樹干液流的主要環(huán)境因子是Solar、Ta和RH,陰天主要是Solar和RH,而雨天則是SWC5—10cm、Rainfall和Ta。對(duì)比同處喀斯特地區(qū)的樹干液流與環(huán)境因子的研究結(jié)果,其差異也很大[16- 18]。

    總的來說,由于地域、物種、生境、所監(jiān)測(cè)的環(huán)境因子、時(shí)間尺度、天氣條件以及季節(jié)等的不同,還有數(shù)據(jù)分析方法的差異等等,即使處于同一生境中的同一物種,各環(huán)境因子對(duì)其樹干液流的影響差異也很大。但是,我們也發(fā)現(xiàn),無論在何種時(shí)間尺度上,Solar(或PAR)和Ta基本上都是影響各樹種樹干液流的主要環(huán)境因子[1,2,10,13- 14,16,18,24,27- 31,34- 35,39],究其原因,Solar(或PAR)作為光合作用的主要驅(qū)動(dòng)力,決定了葉片氣孔的開張程度,從而主要影響著植物蒸騰的強(qiáng)弱[34]。另外,我們還發(fā)現(xiàn),雨天各環(huán)境因子對(duì)樹干液流的解釋程度較晴天和陰天要低的多[13,27,42],本研究亦是如此;這可能是由于雨天Solar弱,RH高,再加上降雨事件的發(fā)生(降雨的歷時(shí)、強(qiáng)度、雨量的大小、間斷性等因素)又對(duì)Solar、PAR、Ta、SWC、WS、Ts等環(huán)境因子產(chǎn)生影響,使得影響樹干液流的環(huán)境因子眾多,另外還有“根壓”的存在[27],使其雨天樹干液流的影響因子較晴天和陰天要復(fù)雜的多;當(dāng)然,這也可能與不同天氣條件下物種自身的生理特性有關(guān)[36],具體原因還有待于進(jìn)一步深入研究;所以,以環(huán)境因子對(duì)雨天的樹干液流進(jìn)行預(yù)測(cè)的時(shí)候可能會(huì)存在誤差,尤其應(yīng)注意不同土層深度的土壤含水量對(duì)其樹干液流的影響。

    [1] 黃德衛(wèi), 張德強(qiáng), 周國(guó)逸, 劉世忠, Dennis O, 李躍林. 鼎湖山針闊葉混交林優(yōu)勢(shì)種樹干液流特征及其與環(huán)境因子的關(guān)系. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(5): 1159- 1166.

    [2] 徐世琴, 吉喜斌, 金博文. 典型荒漠植物沙拐棗莖干液流密度動(dòng)態(tài)及其對(duì)環(huán)境因子的響應(yīng). 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(2): 345- 353.

    [3] Hansen S. Daisy, a flexible Soil-Plant-Atmosphere system Model. 2002. [2017-02- 13].http://www.dina.kvl.dk/daisy/ftp/DaisyDescription.pdf.

    [4] 王文杰, 孫偉, 邱嶺, 祖元?jiǎng)? 劉偉. 不同時(shí)間尺度下興安落葉松樹干液流密度與環(huán)境因子的關(guān)系. 林業(yè)科學(xué), 2012, 48(1): 77- 85.

    [5] Ford C R, Hubbard R M, Kloeppel B D, Vose J M. A comparison of sap flux-based evapotranspiration estimates with catchment-scale water balance. Agricultural and Forest Meteorology, 2007, 145(3/4): 176- 185.

    [6] 張雷, 孫鵬森, 劉世榮. 樹干液流對(duì)環(huán)境變化響應(yīng)研究進(jìn)展. 生態(tài)學(xué)報(bào), 2009, 29(10): 5600- 5610.

    [7] Bleby T M, Mcelrone A J, Jackson R B. Water uptake and hydraulic redistribution across large woody root systems to 20m depth. Plant, Cell and Environment, 2010, 33(12): 2132- 2148.

    [8] 于貴瑞, 王秋鳳. 植物光合、蒸騰與水分利用的生理生態(tài)學(xué). 北京: 科學(xué)出版社, 2010.

    [9] Braun P, Schmid J. Sap flow measurements in grapevines (VitisviniferaL.) 1. stem morphology and use of the heat balance method. Plant and Soil, 1999, 215(1): 39- 45.

    [10] 吳芳, 陳云明, 于占輝. 黃土高原半干旱區(qū)刺槐生長(zhǎng)盛期樹干液流動(dòng)態(tài). 植物生態(tài)學(xué)報(bào), 2010, 34(4): 469- 476.

    [11] Ford C R, Goranson C E, Mitchell R J, Will R E, Teskey R O. Diurnal and seasonal variability in the radial distribution of sap flow: predicting total stem flow in Pinus taeda trees. Tree Physiology, 2004, 24(9): 941- 950.

    [12] 張璇, 張會(huì)蘭, 王玉杰, 王云琦, 劉春霞, 楊坪坪, 潘聲雷. 縉云山典型樹種樹干液流日際變化特征及與氣象因子關(guān)系. 北京林業(yè)大學(xué)學(xué)報(bào), 2016, 38(3): 11- 20.

    [13] 王小菲, 孫永玉, 李昆, 張春華, 李彬. 山合歡樹干液流的季節(jié)變化. 生態(tài)學(xué)雜志, 2013, 32(3): 597- 603.

    [14] 于占輝, 陳云明, 杜盛. 黃土高原半干旱區(qū)人工林刺槐展葉期樹干液流動(dòng)態(tài)分析. 林業(yè)科學(xué), 2009, 45(4): 53- 59.

    [15] 孫慧珍, 周曉峰, 康紹忠. 應(yīng)用熱技術(shù)研究樹干液流進(jìn)展. 應(yīng)用生態(tài)學(xué)報(bào), 2004, 15(6): 1074- 1078.

    [16] 張中峰, 黃玉清, 李先琨, 莫凌, 焦繼飛, 尤業(yè)明. 巖溶區(qū)青岡櫟樹干液流特征及其與環(huán)境因子關(guān)系. 中國(guó)巖溶, 2008, 27(3): 228- 234.

    [17] 黃玉清, 張中峰, 何成新, 趙平, 袁維園, 焦繼飛, 尤業(yè)明. 巖溶區(qū)青岡櫟整樹蒸騰的季節(jié)變化. 應(yīng)用生態(tài)學(xué)報(bào), 2009, 20(2): 256- 264.

    [18] 劉延惠, 舒德遠(yuǎn), 崔迎春, 侯貽菊, 趙文君, 丁訪軍. 茂蘭喀斯特森林亞優(yōu)勢(shì)種短萼海桐樹干液流特征及其環(huán)境因子響應(yīng). 水土保持學(xué)報(bào), 2016, 30(5): 205- 211.

    [19] 朱守謙. 喀斯特森林生態(tài)研究(Ⅱ). 貴陽(yáng): 貴州科技出版社, 1997.

    [20] 俞國(guó)松, 王世杰, 容麗, 冉景丞. 茂蘭喀斯特森林主要演替群落的凋落物動(dòng)態(tài). 植物生態(tài)學(xué)報(bào), 2011, 35(10): 1019- 1028.

    [21] 胡蕖, 魏魯明, 楊成華, 崔迎春, 丁訪軍. 茂蘭喀斯特峰叢中部常綠落葉闊葉混交林群落特征及多樣性研究. 福建林業(yè)科技, 2009, 36(1): 43- 48.

    [22] Granier A. A new method of sap flow measurement in tree stems. Annales Des Sciences Forestieres, 1985, 42(2): 193- 200.

    [23] 陳立欣, 張志強(qiáng), 李湛東, 張文娟, 張曉放, 董克宇, 王國(guó)玉. 大連4種城市綠化喬木樹種夜間液流活動(dòng)特征. 植物生態(tài)學(xué)報(bào), 2010, 34(5):535- 546.

    [24] 李海濤, 向樂, 夏軍, 林耀明, 梁濤. 應(yīng)用熱擴(kuò)散技術(shù)對(duì)亞熱帶紅壤區(qū)濕地松人工林樹干邊材液流的研究. 林業(yè)科學(xué), 2006, 42(10): 31- 38.

    [25] 王華, 歐陽(yáng)志云, 鄭華, 王效科, 倪永明, 任玉芬. 北京城區(qū)常見樹種生長(zhǎng)季樹干液流的時(shí)滯特征. 應(yīng)用生態(tài)學(xué)報(bào), 2009, 20(9): 2111- 2117.

    [26] 杜家菊, 陳志偉. 使用SPSS線性回歸實(shí)現(xiàn)通徑分析的方法. 生物學(xué)通報(bào), 2010, 45(2): 4- 6.

    [27] 曹文強(qiáng), 韓海榮, 馬欽彥, 康峰峰, 藺琛. 山西太岳山遼東櫟夏季樹干液流通量研究. 林業(yè)科學(xué), 2004, 40(2): 174- 177.

    [28] 劉華, 佘春燕, 白志強(qiáng), 李倩, 劉端, 韓燕梁. 不同徑級(jí)的西伯利亞紅松樹干液流及蒸騰耗水特征的差異. 西北植物學(xué)報(bào), 2016, 36(2): 390- 397.

    [29] 凡超, 邱燕萍, 李志強(qiáng), 李建光, 張邦躍, 袁沛元. 荔枝樹干液流速率與氣象因子的關(guān)系. 生態(tài)學(xué)報(bào), 2014, 34(9): 2401- 2410.

    [30] 隋旭紅, 張建軍, 文萬榮. 晉西黃土區(qū)遼東櫟、山楊樹干液流比較研究. 生態(tài)學(xué)報(bào), 2011, 31(16): 4791- 4798.

    [31] 徐先英, 孫保平, 丁國(guó)棟, 郭樹江, 柴成武. 干旱荒漠區(qū)典型固沙灌木液流動(dòng)態(tài)變化及其對(duì)環(huán)境因子的響應(yīng). 生態(tài)學(xué)報(bào), 2008, 28(3): 895- 905.

    [32] 程靜, 歐陽(yáng)旭, 黃德衛(wèi), 劉世忠, 張德強(qiáng), 李躍林. 鼎湖山針闊葉混交林4種優(yōu)勢(shì)樹種樹干液流特征. 生態(tài)學(xué)報(bào), 2015, 35(12): 4097- 4104.

    [33] 王艷兵, 德永軍, 熊偉, 王彥輝, 李振華, 劉千. 華北落葉松夜間樹干液流特征及生長(zhǎng)季補(bǔ)水格局. 生態(tài)學(xué)報(bào), 2013, 33(5): 1375- 1385.

    [34] 趙仲輝, 康文星, 田大倫, 項(xiàng)文化, 閆文德. 湖南會(huì)同杉木液流變化及其與環(huán)境因子的關(guān)系. 林業(yè)科學(xué), 2009, 45(7): 127- 132.

    [35] 熊偉, 王彥輝, 徐德應(yīng). 寧南山區(qū)華北落葉松人工林蒸騰耗水規(guī)律及其對(duì)環(huán)境因子的響應(yīng). 林業(yè)科學(xué), 2003, 29(2): 1- 7.

    [36] 孫鵬森, 馬履一, 王小平, 翟明普. 油松樹干液流的時(shí)空變異性研究. 北京林業(yè)大學(xué)學(xué)報(bào), 2000, 22(5): 1- 6.

    [37] 涂潔, 劉琪璟, 李海濤, 林耀明. 江西千煙洲濕地松生長(zhǎng)旺季樹干液流動(dòng)態(tài)及影響因素分析. 林業(yè)科學(xué), 2008, 44(1): 46- 51.

    [38] 馬玲, 趙平, 饒興權(quán), 蔡錫安, 曾小平, 陸平. 馬占相思樹干液流特征及其與環(huán)境因子的關(guān)系. 生態(tài)學(xué)報(bào), 2005, 25(9): 2145- 2151.

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

    [40] 張建國(guó), 久米朋宣, 大規(guī)恭一, 山中典和, 杜盛. 黃土高原半干旱區(qū)遼東櫟的樹干液流動(dòng)態(tài). 林業(yè)科學(xué), 2011, 47(4): 63- 69.

    [41] 陳寶強(qiáng), 張建軍, 張艷婷, 田寧寧. 晉西黃土區(qū)遼東櫟和山楊樹干液流對(duì)環(huán)境因子的響應(yīng). 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(3): 746- 754.

    [42] 劉春鵬, 朱美秋, 薛帥, 馬長(zhǎng)明, 支恩波, 李向軍. 干旱條件下荊條樹干液流對(duì)不同天氣的響應(yīng). 福建林業(yè)科技, 2014, 41(1): 34- 39.

    CharacteristicsofsapflowandcorrelationanalysiswithenvironmentalfactorsofAcerwangchiiinthekarstarea

    WU Peng1,2, YANG Wenbin1, CUI Yingchun2, ZHAO Wenjun2, SHU Deyuan2, HOU Yiju2, DING Fangjun2,*

    1InstituteofDesertificationStudies,ChineseAcademyofForestry,Beijing100091,China2GuizhouAcademyofForestry,Guiyang550005,China

    The Thermal Dissipation Probe (TDP) method was used to conduct long-term continuous observations from February 2012 to March 2014 on the sap flow velocity ofAcerwangchiiin karst regions. In addition, environmental factors including air temperature (Ta), relative humidity (RH), solar radiation (Solar), wind speed (WS), rainfall, soil water content (SWC), and the characteristics ofA.wangchiisap flow were investigated, and correlations between environmental factors and different temporal scales were determined. The findings revealed the following: 1) sap flow velocity in relation to different weather conditions, in descending order, was sunny, cloudy, and rainy days, and the sap flow was generally high in the daytime and low at night; 2) the average amount of sap flow was 5.08 kg/d, with the sap flow quantity varying in different seasons, in descending order of (8.38 ± 5.32) kg/d in summer, (5.16 ± 3.99) kg/d in autumn, (4.86 ± 3.77) kg/d in spring, and (1.94 ± 1.40) kg/d in winter; 3) the monthly sap flow was 153.64 kg and the annual sap flow was 1838.40 kg, on average; 4) on the hourly scale, solar radiation over 1 day was the main environmental factor influencing sap flow regardless of the weather conditions, while the influence of the other environmental factors on sap flow varied in the daytime, at night, and in different seasons and rainfalls; and 5) over 1 year: on the hourly scale, 63.50% of the sap flow could be interpreted jointly by the environmental factors including solar radiation,Ta, RH, SWC10cmand WS. On the daily scale, 68.50% of the sap flow could be interpreted jointly by the environmental factors including solar radiation,Ta, SWC10cmand rainfall. On the monthly scale, the sap flow quantity could be interpreted with the single environmental factor,Ta, which contributed 74.80% to the total. With the time scale reducing, the more environmental factors were involved in the regression equation, the lesser interpretation of sap flow quantity (R2) would be. In comparison to similar findings from other regions, the impact of environment factors on sap flow diverged greatly. However, on temporal scales, solar radiation or PAR andTawere the two major factors that affected sap flow, and no environmental factor contributed in the interpretation of sap flow on rainy days, which is consistent with the present study. Therefore, estimates of the impact of environmental factors on sap flow on rainy days deviate, and particular attention should be paid to the role of soil water content.

    karst;Acerwangchii; sap flow; environmental factors; different temporal scale

    黔科合人字(2011)12號(hào);黔科合人才團(tuán)隊(duì)[2014]4004號(hào);黔林科合J字〔2015〕14號(hào);黔科合SY字〔2014〕3039號(hào);黔科合院士站〔2014〕4006;黔林科合J字〔2014〕09號(hào)

    2016- 09- 25; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 07- 12

    *通訊作者Corresponding author.E-mail: ding3920034@163.com

    10.5846/stxb201609251934

    吳鵬, 楊文斌, 崔迎春, 趙文君, 舒德遠(yuǎn), 侯貽菊, 丁訪軍.喀斯特區(qū)天峨槭(Acerwangchii)樹干液流特征及其與環(huán)境因子的相關(guān)分析.生態(tài)學(xué)報(bào),2017,37(22):7552- 7567.

    Wu P, Yang W B, Cui Y C, Zhao W J, Shu D Y, Hou Y J, Ding F J.Characteristics of sap flow and correlation analysis with environmental factors ofAcerwangchiiin the karst area.Acta Ecologica Sinica,2017,37(22):7552- 7567.

    猜你喜歡
    陰天液流晴天
    它若安好,便是晴天
    一起來寫“雪”
    太陽(yáng)躲在哪里
    啟蒙(3-7歲)(2020年9期)2020-09-30 05:45:24
    給陰天染上色彩
    晴天有時(shí)下豬
    陰天
    一種全釩液流電池并網(wǎng)控制系統(tǒng)的仿真設(shè)計(jì)
    H2-Fe3+/Fe2+氧化還原液流電池初步研究
    哲理漫畫
    健康必讀(2014年1期)2014-06-23 14:58:22
    我國(guó)—項(xiàng)液流電池國(guó)際標(biāo)準(zhǔn)提案獲準(zhǔn)通過
    观看av在线不卡| 亚洲av.av天堂| 日韩av不卡免费在线播放| 国产精品.久久久| 99视频精品全部免费 在线| 亚洲欧美精品自产自拍| 久久99蜜桃精品久久| 三级国产精品片| 热re99久久精品国产66热6| 亚洲国产精品国产精品| 精品人妻熟女毛片av久久网站| 久久这里有精品视频免费| 免费大片黄手机在线观看| 国产成人精品久久久久久| 国产精品成人在线| 日本猛色少妇xxxxx猛交久久| 在线观看美女被高潮喷水网站| 亚洲综合精品二区| 亚洲婷婷狠狠爱综合网| 嫩草影院新地址| 99热网站在线观看| √禁漫天堂资源中文www| 久久97久久精品| 亚洲精品日本国产第一区| 国产黄色视频一区二区在线观看| 妹子高潮喷水视频| 精品人妻一区二区三区麻豆| 国产欧美亚洲国产| 亚洲精品,欧美精品| 色吧在线观看| 国产成人免费无遮挡视频| 黄色视频在线播放观看不卡| 丝袜喷水一区| 欧美日本中文国产一区发布| 2022亚洲国产成人精品| 一边亲一边摸免费视频| 亚洲精品亚洲一区二区| 亚洲情色 制服丝袜| 亚洲精品国产av成人精品| 日韩一区二区三区影片| 一级毛片aaaaaa免费看小| 性高湖久久久久久久久免费观看| 欧美亚洲 丝袜 人妻 在线| 美女国产视频在线观看| 国产精品偷伦视频观看了| 色吧在线观看| 99热这里只有是精品50| 交换朋友夫妻互换小说| 日日撸夜夜添| 午夜免费鲁丝| 国产91av在线免费观看| 麻豆成人av视频| 免费高清在线观看视频在线观看| 少妇的逼好多水| 18禁在线无遮挡免费观看视频| 丝袜在线中文字幕| 另类亚洲欧美激情| 永久免费av网站大全| 26uuu在线亚洲综合色| 夫妻午夜视频| 99九九在线精品视频 | 欧美成人午夜免费资源| av天堂中文字幕网| 国产精品福利在线免费观看| a级毛片免费高清观看在线播放| 国模一区二区三区四区视频| 日本黄色片子视频| 不卡视频在线观看欧美| 男女边摸边吃奶| 日本午夜av视频| 免费观看a级毛片全部| 最后的刺客免费高清国语| 日韩大片免费观看网站| 国产午夜精品久久久久久一区二区三区| 日韩人妻高清精品专区| 国产女主播在线喷水免费视频网站| 午夜福利网站1000一区二区三区| 一级二级三级毛片免费看| 成人影院久久| 一级a做视频免费观看| 国产黄片美女视频| 夜夜爽夜夜爽视频| 精品亚洲成a人片在线观看| 亚洲欧美精品专区久久| 欧美精品一区二区免费开放| 天天躁夜夜躁狠狠久久av| 亚洲av.av天堂| 啦啦啦在线观看免费高清www| 久久精品久久精品一区二区三区| 春色校园在线视频观看| 91成人精品电影| 亚洲丝袜综合中文字幕| 九九在线视频观看精品| 国产一区二区三区av在线| 美女国产视频在线观看| 91久久精品国产一区二区成人| 99热6这里只有精品| 极品少妇高潮喷水抽搐| 欧美日韩综合久久久久久| 国产 一区精品| 午夜激情久久久久久久| 插阴视频在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 日韩,欧美,国产一区二区三区| 日韩av在线免费看完整版不卡| 中国三级夫妇交换| 18禁在线无遮挡免费观看视频| 性色av一级| 免费黄色在线免费观看| 一级毛片aaaaaa免费看小| 韩国高清视频一区二区三区| 精品亚洲成a人片在线观看| 精品酒店卫生间| 中文天堂在线官网| 免费不卡的大黄色大毛片视频在线观看| 久久午夜福利片| 香蕉精品网在线| 精品一区二区三卡| 国产高清不卡午夜福利| 哪个播放器可以免费观看大片| 国产一区亚洲一区在线观看| 在现免费观看毛片| 亚洲国产毛片av蜜桃av| 三级国产精品欧美在线观看| 久久精品国产自在天天线| 国产亚洲欧美精品永久| 国产精品熟女久久久久浪| 久久狼人影院| 亚洲国产精品一区二区三区在线| 少妇被粗大猛烈的视频| 久久国产亚洲av麻豆专区| 国产美女午夜福利| 国产在线男女| 精品国产露脸久久av麻豆| 久久精品久久久久久噜噜老黄| 一级毛片久久久久久久久女| 欧美 日韩 精品 国产| 国产精品人妻久久久久久| 永久免费av网站大全| 777米奇影视久久| 精品久久久精品久久久| 伊人久久国产一区二区| 色94色欧美一区二区| 亚洲精品,欧美精品| 亚洲内射少妇av| 久久99热这里只频精品6学生| 少妇 在线观看| 下体分泌物呈黄色| 伦精品一区二区三区| 一区在线观看完整版| 国产精品.久久久| 欧美日韩一区二区视频在线观看视频在线| 亚洲av二区三区四区| 免费观看a级毛片全部| 亚洲国产欧美在线一区| 777米奇影视久久| 亚洲欧洲精品一区二区精品久久久 | 国产成人一区二区在线| 国产精品99久久久久久久久| 国产欧美亚洲国产| 国产男女超爽视频在线观看| 成人毛片a级毛片在线播放| 国内少妇人妻偷人精品xxx网站| 国产极品天堂在线| 免费观看a级毛片全部| 日本wwww免费看| xxx大片免费视频| 天美传媒精品一区二区| 王馨瑶露胸无遮挡在线观看| 99热6这里只有精品| 亚洲不卡免费看| 三上悠亚av全集在线观看 | 国产日韩一区二区三区精品不卡 | 成人二区视频| a 毛片基地| 国产精品国产三级国产专区5o| 亚洲欧美日韩另类电影网站| a级毛色黄片| 99热这里只有是精品50| 久久国产精品男人的天堂亚洲 | 国产亚洲av片在线观看秒播厂| 亚洲av成人精品一区久久| 欧美97在线视频| 另类亚洲欧美激情| 黄色一级大片看看| 在线观看人妻少妇| 久久婷婷青草| 交换朋友夫妻互换小说| 亚洲国产欧美在线一区| 一级毛片 在线播放| 能在线免费看毛片的网站| 中文字幕制服av| av有码第一页| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩卡通动漫| 国产爽快片一区二区三区| 黄色配什么色好看| 少妇人妻久久综合中文| 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区蜜桃 | 成人国产av品久久久| 欧美丝袜亚洲另类| 亚洲欧美日韩东京热| 精品久久久久久久久亚洲| 久久97久久精品| 黑人巨大精品欧美一区二区蜜桃 | av在线播放精品| 国产一级毛片在线| 看免费成人av毛片| 日韩制服骚丝袜av| 中文字幕精品免费在线观看视频 | 久久精品国产自在天天线| a级毛片免费高清观看在线播放| 久久国内精品自在自线图片| 日韩一区二区三区影片| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区视频9| 少妇裸体淫交视频免费看高清| 国产美女午夜福利| 欧美+日韩+精品| av播播在线观看一区| 国产黄片视频在线免费观看| 一级av片app| 精品少妇黑人巨大在线播放| 免费黄频网站在线观看国产| 一级a做视频免费观看| 美女内射精品一级片tv| 大香蕉97超碰在线| 亚洲伊人久久精品综合| 2018国产大陆天天弄谢| 久久国产亚洲av麻豆专区| 亚洲精品第二区| 日本91视频免费播放| 女性生殖器流出的白浆| 欧美日韩视频精品一区| 免费观看性生交大片5| 亚洲精品456在线播放app| 日本黄色日本黄色录像| 69精品国产乱码久久久| 久久6这里有精品| 中文字幕免费在线视频6| av网站免费在线观看视频| 一二三四中文在线观看免费高清| 亚洲欧美成人精品一区二区| 国产精品三级大全| 免费少妇av软件| 26uuu在线亚洲综合色| 最近2019中文字幕mv第一页| 天堂中文最新版在线下载| 一区二区三区四区激情视频| 欧美高清成人免费视频www| 中国三级夫妇交换| 亚洲av在线观看美女高潮| 婷婷色麻豆天堂久久| 欧美日韩视频精品一区| 午夜老司机福利剧场| 另类精品久久| 欧美日韩亚洲高清精品| a级毛片免费高清观看在线播放| 啦啦啦中文免费视频观看日本| 国产精品久久久久久av不卡| 精品人妻熟女av久视频| av在线app专区| 日韩精品免费视频一区二区三区 | 另类亚洲欧美激情| 内射极品少妇av片p| 国产一区二区三区综合在线观看 | 两个人免费观看高清视频 | 国产一区二区在线观看日韩| 欧美日本中文国产一区发布| 日韩欧美一区视频在线观看 | 国产又色又爽无遮挡免| 色哟哟·www| 亚洲国产日韩一区二区| 人妻制服诱惑在线中文字幕| 国产精品不卡视频一区二区| 国产高清三级在线| 高清在线视频一区二区三区| 一级爰片在线观看| 久久久a久久爽久久v久久| 精品亚洲成a人片在线观看| 午夜91福利影院| 韩国高清视频一区二区三区| 男人狂女人下面高潮的视频| 久久国产精品男人的天堂亚洲 | 日韩av免费高清视频| 91在线精品国自产拍蜜月| 午夜激情福利司机影院| 高清在线视频一区二区三区| 内射极品少妇av片p| 最新的欧美精品一区二区| 9色porny在线观看| av在线老鸭窝| 在线观看av片永久免费下载| 精品少妇久久久久久888优播| 亚洲国产欧美日韩在线播放 | 我要看日韩黄色一级片| 亚洲国产欧美日韩在线播放 | 亚洲怡红院男人天堂| 国产视频首页在线观看| 日本vs欧美在线观看视频 | 99精国产麻豆久久婷婷| 自拍偷自拍亚洲精品老妇| 十分钟在线观看高清视频www | 桃花免费在线播放| 免费观看性生交大片5| 国产无遮挡羞羞视频在线观看| 一级毛片aaaaaa免费看小| 精品熟女少妇av免费看| 日韩av在线免费看完整版不卡| 2018国产大陆天天弄谢| 国产片特级美女逼逼视频| 少妇的逼水好多| 中文欧美无线码| av国产久精品久网站免费入址| 成人综合一区亚洲| 亚洲美女黄色视频免费看| 性色avwww在线观看| 18+在线观看网站| 五月伊人婷婷丁香| 亚洲精品日本国产第一区| 啦啦啦中文免费视频观看日本| 街头女战士在线观看网站| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 九九在线视频观看精品| 久久影院123| 在线观看www视频免费| 黑丝袜美女国产一区| 日本猛色少妇xxxxx猛交久久| 亚洲综合精品二区| 亚洲人成网站在线观看播放| 国产欧美日韩综合在线一区二区 | 亚洲av免费高清在线观看| av在线app专区| 看免费成人av毛片| 老司机影院毛片| 少妇丰满av| 日日摸夜夜添夜夜爱| 亚洲av不卡在线观看| 欧美精品高潮呻吟av久久| 久久久久久久久久久免费av| 黄色视频在线播放观看不卡| 免费大片黄手机在线观看| 国产男女内射视频| 亚洲国产精品一区三区| 亚洲人与动物交配视频| 丁香六月天网| 亚洲真实伦在线观看| 久久99热这里只频精品6学生| 国产乱人偷精品视频| 高清黄色对白视频在线免费看 | 大香蕉97超碰在线| tube8黄色片| 精品国产乱码久久久久久小说| 午夜91福利影院| 亚洲av在线观看美女高潮| 美女中出高潮动态图| 99热这里只有精品一区| 人人妻人人爽人人添夜夜欢视频 | 国产片特级美女逼逼视频| 国产亚洲一区二区精品| 国产在线男女| 久久久欧美国产精品| 这个男人来自地球电影免费观看 | 麻豆成人av视频| 一级毛片黄色毛片免费观看视频| 国产成人精品无人区| 国产亚洲91精品色在线| 国产片特级美女逼逼视频| 又黄又爽又刺激的免费视频.| 日韩,欧美,国产一区二区三区| 2022亚洲国产成人精品| 91久久精品国产一区二区成人| 99热这里只有精品一区| 亚洲成人av在线免费| 成人黄色视频免费在线看| 国产成人精品久久久久久| 2022亚洲国产成人精品| 亚洲高清免费不卡视频| 美女视频免费永久观看网站| 三上悠亚av全集在线观看 | 18禁裸乳无遮挡动漫免费视频| 精品久久久精品久久久| 黄色毛片三级朝国网站 | 久久久久精品性色| 精品一区二区免费观看| 男的添女的下面高潮视频| 18禁裸乳无遮挡动漫免费视频| 亚洲,一卡二卡三卡| 精品99又大又爽又粗少妇毛片| 午夜免费男女啪啪视频观看| 国产精品一区二区性色av| 777米奇影视久久| 亚洲国产最新在线播放| 中文字幕人妻熟人妻熟丝袜美| av天堂中文字幕网| 亚洲欧美精品专区久久| 久久韩国三级中文字幕| 久久精品熟女亚洲av麻豆精品| 国内揄拍国产精品人妻在线| 久久精品国产亚洲av天美| 国产伦精品一区二区三区视频9| av在线app专区| 日韩三级伦理在线观看| 80岁老熟妇乱子伦牲交| 自线自在国产av| 亚洲激情五月婷婷啪啪| 乱系列少妇在线播放| 亚洲成色77777| 国产一区二区在线观看av| 久久久久久久久大av| 高清av免费在线| 在线亚洲精品国产二区图片欧美 | 欧美老熟妇乱子伦牲交| 日韩大片免费观看网站| 国产一级毛片在线| 日韩av在线免费看完整版不卡| 春色校园在线视频观看| 久久久久国产网址| 黑人猛操日本美女一级片| 在线观看人妻少妇| 3wmmmm亚洲av在线观看| 免费在线观看成人毛片| 九草在线视频观看| 男的添女的下面高潮视频| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 亚洲av.av天堂| videos熟女内射| 亚洲天堂av无毛| 亚洲人与动物交配视频| 国产免费福利视频在线观看| 亚洲高清免费不卡视频| 大香蕉久久网| 国产亚洲欧美精品永久| 乱人伦中国视频| 欧美bdsm另类| 特大巨黑吊av在线直播| 最新中文字幕久久久久| 狠狠精品人妻久久久久久综合| 久久热精品热| 一本久久精品| 亚洲成色77777| 国产成人a∨麻豆精品| 欧美一级a爱片免费观看看| 国产成人免费观看mmmm| 亚洲av综合色区一区| 成人毛片a级毛片在线播放| 国模一区二区三区四区视频| 黄色毛片三级朝国网站 | 一本—道久久a久久精品蜜桃钙片| 精品视频人人做人人爽| 青青草视频在线视频观看| 免费看光身美女| 在线观看美女被高潮喷水网站| 国产伦精品一区二区三区四那| 国产午夜精品久久久久久一区二区三区| 国产午夜精品久久久久久一区二区三区| 王馨瑶露胸无遮挡在线观看| 伊人久久精品亚洲午夜| 97在线视频观看| 精品亚洲成a人片在线观看| 人人妻人人澡人人看| 免费看日本二区| 男人舔奶头视频| 色婷婷久久久亚洲欧美| 在线天堂最新版资源| 亚洲欧洲精品一区二区精品久久久 | 午夜91福利影院| 精品酒店卫生间| 免费观看性生交大片5| 色5月婷婷丁香| 我要看日韩黄色一级片| 国产欧美日韩精品一区二区| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕| 一区二区av电影网| 国产爽快片一区二区三区| 亚洲美女黄色视频免费看| 国产精品国产三级国产专区5o| 日韩欧美 国产精品| 亚洲国产精品成人久久小说| 国产男女超爽视频在线观看| 在线观看美女被高潮喷水网站| 91精品一卡2卡3卡4卡| 91久久精品国产一区二区三区| 日韩一区二区视频免费看| 少妇的逼好多水| 大又大粗又爽又黄少妇毛片口| 在线观看三级黄色| 一个人免费看片子| 看非洲黑人一级黄片| 国产精品久久久久久精品古装| 国产极品天堂在线| 久久97久久精品| 国语对白做爰xxxⅹ性视频网站| 黄色配什么色好看| 在线观看美女被高潮喷水网站| av视频免费观看在线观看| 国产男女内射视频| 日韩精品免费视频一区二区三区 | a 毛片基地| 久久免费观看电影| 成年人午夜在线观看视频| 国产免费福利视频在线观看| 日本欧美国产在线视频| 在线看a的网站| 久久久久精品性色| 男女边摸边吃奶| 秋霞在线观看毛片| 只有这里有精品99| av女优亚洲男人天堂| 国产精品久久久久久久久免| 日本与韩国留学比较| 中文字幕亚洲精品专区| 精品国产乱码久久久久久小说| 午夜老司机福利剧场| 夜夜看夜夜爽夜夜摸| 黄色毛片三级朝国网站 | 国产亚洲欧美精品永久| 国产又色又爽无遮挡免| 在线观看国产h片| 不卡视频在线观看欧美| 91精品国产国语对白视频| 久久鲁丝午夜福利片| 国产高清有码在线观看视频| 国产成人一区二区在线| 久久99热这里只频精品6学生| 欧美人与善性xxx| 美女福利国产在线| 日韩av在线免费看完整版不卡| 九九爱精品视频在线观看| √禁漫天堂资源中文www| 人妻制服诱惑在线中文字幕| 下体分泌物呈黄色| 搡女人真爽免费视频火全软件| 国产日韩欧美在线精品| 久久青草综合色| 永久免费av网站大全| 香蕉精品网在线| 五月玫瑰六月丁香| 九草在线视频观看| 欧美精品亚洲一区二区| 菩萨蛮人人尽说江南好唐韦庄| 久久久欧美国产精品| 99久久综合免费| 最新的欧美精品一区二区| 精品亚洲成a人片在线观看| 欧美性感艳星| 国产免费福利视频在线观看| 亚洲av不卡在线观看| 色哟哟·www| 久久女婷五月综合色啪小说| 人妻系列 视频| 最近最新中文字幕免费大全7| 91精品一卡2卡3卡4卡| 夜夜骑夜夜射夜夜干| 亚洲精品日韩av片在线观看| 久久婷婷青草| 王馨瑶露胸无遮挡在线观看| 在线观看免费日韩欧美大片 | 天堂中文最新版在线下载| 色视频www国产| 国产男人的电影天堂91| 女性生殖器流出的白浆| 日韩 亚洲 欧美在线| 成人国产麻豆网| 欧美高清成人免费视频www| 婷婷色综合www| 一区二区三区精品91| 一级爰片在线观看| 久久 成人 亚洲| 伊人亚洲综合成人网| 狠狠精品人妻久久久久久综合| 一级毛片电影观看| 国产有黄有色有爽视频| 嘟嘟电影网在线观看| 亚洲欧洲日产国产| 久久久国产精品麻豆| 三级国产精品片| 在线观看三级黄色| 久久久久久久久久久免费av| 国产乱来视频区| 七月丁香在线播放| 日本黄色日本黄色录像| 最近2019中文字幕mv第一页| www.色视频.com| 欧美激情国产日韩精品一区| av国产精品久久久久影院| 人妻夜夜爽99麻豆av| 高清毛片免费看| 美女xxoo啪啪120秒动态图| 在线看a的网站| 老熟女久久久| 观看美女的网站| 午夜福利在线观看免费完整高清在| 国产一级毛片在线| 久久久精品94久久精品| 99视频精品全部免费 在线| 亚洲国产欧美日韩在线播放 | 不卡视频在线观看欧美| 亚洲人与动物交配视频| 一本色道久久久久久精品综合| 伊人久久精品亚洲午夜| 51国产日韩欧美| 女人精品久久久久毛片| 亚洲国产精品999| 亚洲欧美中文字幕日韩二区| 国产精品秋霞免费鲁丝片| 高清视频免费观看一区二区| a 毛片基地| 丰满人妻一区二区三区视频av| 精品人妻熟女毛片av久久网站| 久久久久精品性色|