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

    祁連山青海云杉樹干液流密度的優(yōu)勢度差異

    2017-06-22 14:04:27萬艷芳于澎濤劉賢德王順利王彥輝
    生態(tài)學(xué)報(bào) 2017年9期
    關(guān)鍵詞:液流云杉晴天

    萬艷芳, 于澎濤,劉賢德,王順利,王彥輝,熊 偉

    1 甘肅農(nóng)業(yè)大學(xué)林學(xué)院, 蘭州 730070 2 中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所,國家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室, 北京 100091 3 甘肅省祁連山水源涵養(yǎng)林研究院, 張掖 734000

    ?

    祁連山青海云杉樹干液流密度的優(yōu)勢度差異

    萬艷芳1,2, 于澎濤2,*,劉賢德1,3,王順利3,王彥輝2,熊 偉2

    1 甘肅農(nóng)業(yè)大學(xué)林學(xué)院, 蘭州 730070 2 中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所,國家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室, 北京 100091 3 甘肅省祁連山水源涵養(yǎng)林研究院, 張掖 734000

    以祁連山排露溝小流域青海云杉林為研究對象,選取有代表性的優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木各3—5株,2015年6月16日至10月14日應(yīng)用熱擴(kuò)散技術(shù)對不同優(yōu)勢度青海云杉樹干液流密度進(jìn)行測定,并同步測定相關(guān)的林外氣象因子。結(jié)果表明:(1)青海云杉液流密度呈晝高夜低趨勢,晴天液流密度變化幅度較大,而陰雨天變化幅度較小。(2)晴天樹木優(yōu)勢度越大,其液流在日內(nèi)的啟動(dòng)越早,結(jié)束越晚,峰值也越大;優(yōu)勢木的平均液流密度為(0.0758±0.0475) mL cm-2min-1,是亞優(yōu)勢木的1.5倍,是中等木和被壓木的1.68倍。(3)青海云杉平均液流密度基本呈現(xiàn)6月份最大,其次是8月份,9、10月份明顯減小,且優(yōu)勢木>亞優(yōu)勢木>中等木>被壓木。(4)相關(guān)性分析和逐步回歸表明,青海云杉日均液流密度與太陽輻射強(qiáng)度、飽和水氣壓差和空氣溫度呈正相關(guān)關(guān)系,與空氣相對濕度和降雨量呈負(fù)相關(guān)關(guān)系。影響優(yōu)勢木、亞優(yōu)勢木和中等木液流密度的主要?dú)庀笠蜃邮翘栞椛鋸?qiáng)度,被壓木液流密度主要受空氣相對濕度的影響。

    樹干液流;優(yōu)勢度;氣象因子;青海云杉林;祁連山

    對干旱地區(qū)的森林而言,蒸散是樹木、林下植被和地面向大氣輸送的水汽總量,是森林生態(tài)系統(tǒng)的最大水分輸出分量。其中,林木蒸騰是蒸散的最重要組分,在土壤-植被-大氣水分傳輸過程中起著重要作用,也對流域產(chǎn)水量變化有著重要的影響。

    樹干液流測定是一種準(zhǔn)確估計(jì)樹木單株蒸騰量的方法,并結(jié)合邊材面積、胸徑或葉面積等空間純量,可以實(shí)現(xiàn)由單株液流向林分蒸騰的尺度擴(kuò)展[1- 3]。目前測定樹干液流的方法有很多,包括莖部熱量平衡法、樹干熱平衡法、熱脈沖法、熱擴(kuò)散方法等[4],其中熱擴(kuò)散方法具有能夠長期地連續(xù)監(jiān)測、時(shí)間分辨率高,同時(shí)具有對樣木損傷較少、數(shù)據(jù)采集自動(dòng)化等優(yōu)點(diǎn),是測定樹干液流的常用方法之一[3]。

    國內(nèi)外研究者應(yīng)用熱擴(kuò)散方法測定樹干液流,研究單株樹木的樹干液流時(shí)空動(dòng)態(tài)[5- 7]和樹干液流與影響因子之間的關(guān)系,如樹木形態(tài)特征[7- 8]、氣象因子[7- 8]和土壤供水水平[9- 11]等因子對樹干液流的影響,這些研究主要通過測定一株或幾株樣樹推算林分液流,并假設(shè)相同徑階樹木的樹干傳輸水分的能力是一致的,再根據(jù)樣地內(nèi)樹木的邊材面積或胸徑來推算樣地的總蒸騰量[1- 3],但此法忽略了其它樹形因子的影響。事實(shí)上,在天然林內(nèi),樹木個(gè)體之間在形態(tài)和立地環(huán)境上存在很大的差異,如優(yōu)勢度的不同[12- 14],導(dǎo)致樹木接受的光照不同,可能會影響單株樹木的液流,進(jìn)而會影響單株樹木以及整個(gè)林分的蒸騰[9]。因此,有必要對樹木進(jìn)行優(yōu)勢度劃分,進(jìn)行樹干液流密度的研究。

    青海云杉林(Piceacrassifolia)是祁連山的主要建群種,呈斑塊狀分布在海拔2500—3300 m的陰坡、半陰坡上,對涵養(yǎng)水源、凈化水質(zhì)、保持水土具有重要的作用[15]。但目前針對青海云杉的蒸騰研究較少,僅見Chang和魏瀟等[11,16- 17]分析了單株青海云杉樹干液流徑向變化以及液流時(shí)空動(dòng)態(tài)及其影響因素,而未考慮樹木優(yōu)勢度差異對液流的影響。因此,本研究選取位于祁連山排露溝小流域的青海云杉,研究樹干液流密度的優(yōu)勢度差異,為準(zhǔn)確估算青海云杉林分蒸騰量奠定基礎(chǔ)。

    1 試驗(yàn)區(qū)概況

    試驗(yàn)區(qū)位于祁連山排露溝小流域(38°33′58.7″N,100°17′18.0″E,海拔范圍2600—3800 m),屬于溫帶氣候,在海拔2600 m處,年均氣溫0.5℃,年均降水量435.5 mm,且集中在5—9月份,年均蒸發(fā)量1051.7 mm,年日照時(shí)數(shù)1892.6 h,年均相對濕度60%。試驗(yàn)區(qū)森林類型單一,主要是以青海云杉(Piceacrassifolia)林為主的寒溫性針葉林,呈斑塊狀分布在陰坡、半陰坡;陽坡以草地為主,零星分布有祁連圓柏(Sabinaprezewalskii)林和灌叢;灌叢優(yōu)勢種有鬼箭錦雞兒(Caraganajubata)、高山柳(Salixcupularis)、銀露梅(Dasiphoradavurica)等。草本主要有苔草(Koeleriacrisrata)、馬藺(Irisensata)、馬先蒿(Pedicularisspp.)等。土壤主要有山地森林灰褐土、山地栗鈣土、草甸土、亞高山灌叢草甸土和高山寒漠土等類型。

    2 研究材料與方法

    2.1 樣地概況

    在海拔2700 m處選取了一塊大小為25 m×50 m的蘚類青海云杉林樣地,其坡向350°,位于坡面的中下部,平均坡度為23°;該林分密度為1128株/hm2,林齡86年,郁閉度0.55,平均樹高(10.6±5.6) m,平均胸徑(15.3±9.1) cm,平均冠幅(3.6±1.1) m。林下灌木層主要由銀露梅(Dasiphoradavurica)組成,蓋度在4%左右;草本層有披針苔草(Koeleriacrisrata)、蘚生馬先蒿(Pedicularisspp.)、棘豆(Oxytropiskansuensis)等,蓋度在37%左右;苔蘚層平均厚度為(8.4±3.1) cm,蓋度在30%左右。

    2.2 樣樹選取

    依據(jù)樣地內(nèi)樹木在林冠內(nèi)的相對位置,按照Kraft樹冠優(yōu)勢度分級標(biāo)準(zhǔn)[18],將樣地內(nèi)青海云杉按林木個(gè)體優(yōu)勢度分為5級:I代表優(yōu)勢木,Ⅱ代表亞優(yōu)勢木,Ⅲ代表中等木,Ⅳ代表被壓木,V代表瀕死木。對樣地內(nèi)的樹木進(jìn)行每木檢尺,并記錄其優(yōu)勢度等級。

    根據(jù)樣地內(nèi)青海云杉的樹木優(yōu)勢度分布,在每個(gè)優(yōu)勢度級別中選擇3—5棵樹作為樣樹(枯死木除外),共選取了13株樣樹(表1)。樣樹的胸徑介于10—38 cm,樹高介于6—21 m,冠幅介于2—7 m。其中,優(yōu)勢木的平均胸徑、樹高和冠幅都最大,分別為(28.18±6.93) cm、(19.4±1.0) m和(5.2±0.8) m,其次是亞優(yōu)勢木和中等木,被壓木都最小,分別為10.52 cm、6.1 m和3.6 m。

    基于該實(shí)驗(yàn)地建立的青海云杉邊材寬度與胸徑之間的函數(shù)關(guān)系[19](式1)計(jì)算出每株樣樹的邊材寬度:

    y=22.19×e0.025x(R2=0.85,n=17)

    (1)

    式中,y為青海云杉林木邊材寬度(mm),x為青海云杉林木的胸徑(cm)。

    本研究選取的13株樣樹2的邊材寬度介于28—57 mm(表1),其中,優(yōu)勢木的邊材面積最大,平均值為(45.45±8.15) mm,介于39.39—56.88 mm,其次是亞優(yōu)勢木((34.07±1.08) mm)和中等木((30.81±1.61) mm),被壓木最小,為28.87 mm。

    表1 樣樹特征參數(shù)

    2.3 樹干液流密度測定

    樹干液流密度用德國Ecomatik公司生產(chǎn)的SF-L型熱擴(kuò)散液流計(jì)進(jìn)行測定。SF-L樹干液流測定儀是基于熱擴(kuò)散原理設(shè)計(jì)的,它由4個(gè)探針(記為S0、S1、S2和S3)組成,探針長度為20 mm,所測出的液流密度是線平均值(有別于熱脈沖測定系統(tǒng)的點(diǎn)平均值)[4,20],探針基本上能整合青海云杉液流密度隨邊材厚度變化而可能出現(xiàn)的變異(即液流密度在邊材內(nèi)的徑向變化),測定值反映的是平均液流密度。測定探針之間的溫差,再由溫差推算出液流密度(Js,mL cm-2min-1)(式2和3):

    (2)

    式中,dtmax是液流密度為0時(shí)的dtact,dtact的計(jì)算公式如下(式3):

    (3)

    式中,T1-0、T1-2和T1-3分別為探針S1和S0、S1和S2、S1和S3之間的溫差(℃)。

    數(shù)據(jù)采集時(shí)間間隔為10 min,測定時(shí)間為2015年6月16日至10月14日。

    2.4 環(huán)境因子測定

    在試驗(yàn)樣地內(nèi)設(shè)有美國Campbell公司生產(chǎn)的CR3000自動(dòng)氣象站,用于連續(xù)采集森林冠層上方降水(P,mm)、太陽輻射強(qiáng)度(Rs,w/m2)、空氣溫度(T,℃)、空氣相對濕度(Rh,%)和林內(nèi)土壤濕度(Ms,%)和土壤溫度(Ts,℃)等數(shù)據(jù)。其中冠層上方環(huán)境因子的觀測高度為30 m,土壤體積含水率分別在0—10、10—20、20—40、40—60 cm和60—80 cm5個(gè)土層進(jìn)行測定,土壤溫度分別在10、20、40、60 cm和80 cm5個(gè)土層進(jìn)行測定。數(shù)據(jù)采集時(shí)間間隔為10 min。

    根據(jù)空氣溫度和空氣相對濕度,計(jì)算飽和水汽壓差(VPD,kPa)如下(式4和式5):

    E=0.611×e

    (4)

    (5)

    式中,E為飽和水汽壓(kPa);T為空氣溫度(℃);Rh為空氣相對濕度(%)。

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

    利用Excel 2007和SPSS 19.0對數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析和作圖。

    3 結(jié)果與分析

    3.1 典型天氣下林外氣象因子日變化

    太陽輻射強(qiáng)度呈晝高夜低趨勢,白天太陽輻射強(qiáng)度變化曲線呈單峰型或多峰型,夜間太陽輻射強(qiáng)度非常微弱,一般保持在5 W/m2以下,太陽輻射強(qiáng)度晴天>陰天>雨天(圖1)。如8月4日晴天,太陽輻射強(qiáng)度日變化曲線為典型的拋物線,太陽輻射從6:10開始啟動(dòng),之后快速上升,于13:50到達(dá)峰值,峰值為1340.17 W/m2,之后快速降低,于20:50到達(dá)谷值。8月19日陰天和8月7日雨天,太陽輻射強(qiáng)度呈多峰形;陰雨天太陽輻射啟動(dòng)時(shí)間和結(jié)束時(shí)間與晴天基本相同,但出現(xiàn)峰值的次數(shù)與峰值間差異較大;8月19日陰天出現(xiàn)3個(gè)持續(xù)時(shí)間較短且數(shù)值較小的峰值,于11:30到達(dá)第1個(gè)峰值,峰值為1048.59 W/m2,比晴天小291.58 W/m2;而8月7日雨天出現(xiàn)2個(gè)峰值,于13:40到達(dá)第1個(gè)峰值,峰值為522.59 W/m2,明顯小于晴天和陰天的峰值。

    在不同天氣條件下,飽和水氣壓差在一日之內(nèi)的變化幅度并不相同(圖1)。0:00—8:00期間飽和水氣壓差呈現(xiàn)出雨天>晴天>陰天,8:00—24:00期間飽和水氣壓差呈現(xiàn)出晴天>陰天>雨天。晴天與陰天飽和水氣壓差變化趨勢基本相同,呈晝高夜低的變化趨勢,但陰天的飽和水氣壓差略小于晴天的;雨天的飽和水氣壓差曲線與晴天和陰天差別較大,飽和水氣壓差在0:00—8:10期間緩慢降低,8:10—24:00期間變化幅度較小,基本維持在0.23 kPa的較低水平,明顯低于該時(shí)段晴天與陰天的飽和水氣壓差。

    圖1 典型天氣下太陽輻射強(qiáng)度和飽和水氣壓差日變化Fig.1 Daily courses of solar radiation intensity and vapor pressure deficit in different weather condition

    3.2 典型天氣下樹干液流密度日變化

    液流密度日變化呈現(xiàn)明顯的晝夜變化規(guī)律(圖2),即白天液流密度明顯高于清晨和夜晚,變幅較大;而清晨和晚上的液流密度緩慢,但并不為0,且變幅相對較小;這與太陽輻射強(qiáng)度日變化曲線相似(圖1)。晴天(8月4日)液流密度呈寬峰型曲線,太陽輻射強(qiáng)度從6:10啟動(dòng),液流密度則從7:00啟動(dòng),滯后于太陽輻射50 min,之后液流快速上升,于12:20到達(dá)峰值,峰值為(0.1107±0.0511) mL cm-2min-1,之后液流緩慢減小,19:00左右急劇下降,于22:20左右降低到夜間水平,而太陽輻射強(qiáng)度于20:50降低到5 W/m2以下,液流結(jié)束時(shí)間明顯比太陽輻射滯后1.5 h,這表明液流的啟動(dòng)和結(jié)束與太陽輻射強(qiáng)度緊密相關(guān)。

    天氣條件影響液流密度,陰天(8月19日)和雨天(8月7日)的液流曲線明顯不同于晴天(圖2),陰天和雨天液流呈窄峰型曲線,陰天液流曲線變化比較平緩,雨天液流曲線變化幅度相對較小。而且,陰天的啟動(dòng)時(shí)間比晴天滯后30 min,雨天明顯比晴天推遲6 h;陰天((0.0824±0.0369) mL cm-2min-1)和雨天((0.0524±0.0229) mL cm-2min-1)的峰值明顯小于晴天的.統(tǒng)計(jì)表明,晴天、陰天和雨天的平均液流密度分別為(0.0590±0.0320) mL cm-2min-1、(0.0363±0.0281) mL cm-2min-1和(0.0250±0.0086) mL cm-2min-1,平均夜間液流密度分別為(0.0291±0.0136) mL cm-2min-1、(0.0124±0.0083) mL cm-2min-1和(0.0207±0.0113) mL cm-2min-1。

    圖2 典型天氣下樹干液流密度日變化Fig.2 Daily courses of sap flow density in different weather condition

    3.3 不同優(yōu)勢度青海云杉樹干液流密度日變化

    8月4日晴天,日最高氣溫20.25℃,最低氣溫11.93℃,平均相對濕度38.14%,7月31日至8月3日期間降雨較多,故8月4日不同優(yōu)勢度樹木所處林地的土壤水分供應(yīng)充足(21.4±2.6)%,可以忽略土壤水分條件對不同優(yōu)勢度青海云杉液流日變化造成的差異。

    晴天優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木液流密度呈單峰型曲線。其中,優(yōu)勢木、亞優(yōu)勢木和中等木液流密度變化幅度較大,而被壓木晝夜變化不明顯(圖3)。液流從清晨7:00—8:00啟動(dòng),12:00—13:00達(dá)到峰值,20:00—23:00降到夜間水平,但不為0。

    優(yōu)勢度越大,液流啟動(dòng)越早,結(jié)束時(shí)間越晚,液流時(shí)間跨度越長。即優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木的液流啟動(dòng)時(shí)間分別在7:00、7:10、7:30和7:50,結(jié)束時(shí)間分別在22:40、22:10、21:10和20:00。而且,優(yōu)勢木、亞優(yōu)勢木和中等木出現(xiàn)峰值的時(shí)間相近,于12:00—13:00到達(dá)峰值,但峰值差異較大,優(yōu)勢木峰值為(0.1487±0.0294) mL cm-2min-1,是亞優(yōu)勢木的1.3倍,是中等木的1.9倍;而被壓木則表現(xiàn)為一個(gè)略有波動(dòng)的寬平曲線,峰值不明顯,從10:40開始一直維持到17:00左右,液流密度基本保持在0.055 mL cm-2min-1左右。統(tǒng)計(jì)表明,優(yōu)勢木的平均液流密度最大,為(0.0758±0.0475) mL cm-2min-1,是亞優(yōu)勢木的1.5倍、是中等木和被壓木的1.68倍。但被壓木平均夜間液流密度最大,為(0.0419±0.0029) mL cm-2min-1,是中等木和優(yōu)勢木的1.5倍,是亞優(yōu)勢木的1.8倍。

    圖3 8月4日(晴天)不同優(yōu)勢度青海云杉樹干液流密度日變化Fig.3 Daily courses of sap flow density for each crown class of Picea crassifolia on August 4, 2015

    3.4 不同優(yōu)勢度青海云杉樹干液流密度季節(jié)變化

    優(yōu)勢木平均液流密度6月份最大,其次是8月份,9、10月份明顯減小,且較接近,二者相差1.81 mL cm-2d-1(表2)。亞優(yōu)勢木、中等木和被壓木平均液流密度與優(yōu)勢木的變化趨勢基本一致,但被壓木8月份稍有反彈,比6月份的平均液流密度大2.80 mL cm-2d-1。另外,隨著優(yōu)勢度的減小,液流密度月均值呈現(xiàn)減小的趨勢,即優(yōu)勢木>亞優(yōu)勢木>中等木>被壓木。優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木的液流密度月均值之間的差異也隨著生長季的進(jìn)程而減小。

    3.5 青海云杉日均液流密度與林外氣象因子的關(guān)系

    青海云杉日均液流密度與太陽輻射強(qiáng)度(Rs)、飽和水氣壓差(VPD)和空氣溫度(T)呈正相關(guān)關(guān)系,而與空氣相對濕度(Rh)和降雨量(P)呈負(fù)相關(guān)關(guān)系(表3)。

    林外氣象因子對日均液流密度的影響程度隨著優(yōu)勢度的減小而減小(即相對系數(shù)R2減小)。其中,優(yōu)勢木、亞優(yōu)勢木和中等木日均液流密度與Rs、VPD和T呈極顯著的正相關(guān)關(guān)系,而被壓木與VPD呈顯著的正相關(guān),與Rs和T僅呈正相關(guān),且相關(guān)系數(shù)較??;青海云杉日均液流密度與Rh均呈極顯著的負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)的絕對值隨優(yōu)勢度的減小而減小;優(yōu)勢木和亞優(yōu)勢木日均液流密度與降雨過程中的降雨量呈極顯著的負(fù)相關(guān),而中等木與P僅呈顯著的負(fù)相關(guān),被壓木與P呈負(fù)相關(guān)關(guān)系??傮w來看,優(yōu)勢木日均液流密度受林外氣象因子的影響最大,其次是亞優(yōu)勢木和中等木,被壓木最小。

    表2 不同月份日均液流密度變化(mL cm-2 d-1)

    (1)數(shù)據(jù)表示為:平均值±S.D.;(2) 7月份為29日—31日的平均值.

    表3 液流密度與氣象因子的Pearson相關(guān)系數(shù)

    **, *分別表示該值在0.01 和 0.05 水平上差異顯著

    多元逐步回歸分析表明,影響青海云杉液流密度的主要?dú)庀笠蜃右驑淠緝?yōu)勢度等級而異(表4)。

    影響優(yōu)勢木、亞優(yōu)勢木和中等木液流密度的主要?dú)庀笠蜃邮翘栞椛鋸?qiáng)度,被壓木液流密度主要受空氣相對濕度的影響。

    表4 液流密度與氣象因子的多元回歸模型

    4 討論

    本研究中青海云杉林在個(gè)體形態(tài)(胸徑10—38 cm、樹高6—21 m)和功能上分化明顯,13株被測樣樹的平均液流密度介于12.84—43.35 mL cm-2d-1.而且,對于胸徑差別較小,優(yōu)勢度等級為中等木(DBH 10.43 cm)和被壓木(DBH 10.52 cm)的2株樣樹,液流密度的平均值分別為(21.80±18.68) mL cm-2d-1和(12.24±8.84) mL cm-2d-1,分別介于2.64—92.48 mL cm-2d-1和0.68—36.20 mL cm-2d-1,說明樹木優(yōu)勢度較大有利于液流產(chǎn)生,其原因可能與樹木優(yōu)勢度大時(shí)接受的太陽輻射能量較多有關(guān)。

    青海云杉樹干液流日動(dòng)態(tài)因樹木優(yōu)勢度等級而異。晴天液流密度呈單峰型曲線,液流從7:00—8:00啟動(dòng),12:00—13:00出現(xiàn)最大值,20:00—23:00降到夜間水平,但不為0。這與銀白楊(PopulusalbaL.×P.talassica)[7]、紅松(Pinuskoraiensis)[13]、木荷(Schimasuperba)[14]和冷杉(Abiesfabri)[21]的晴天液流曲線變化趨勢一致,均呈單峰曲線。而相比于華北落葉松(Larixprincipis-rupprechtii)[22]液流的液流曲線,青海云杉未出現(xiàn)“午休”現(xiàn)象。這是由于液流密度的日變化與大多數(shù)樹種葉片光合速率的日進(jìn)程基本同步[23],而青海云杉在午間氣孔開度可能處于變小狀態(tài),但飽和水氣壓差增大,從而促進(jìn)了青海云杉的蒸騰作用,使得由于氣孔關(guān)閉而引起的“午休”現(xiàn)象未出現(xiàn)。優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木的液流密度日變化趨勢基本一致,但樹木優(yōu)勢度越大,樹干液流啟動(dòng)時(shí)間越早,結(jié)束時(shí)間越晚,這是由于樹木優(yōu)勢度越大,樹木在林分中占據(jù)的空間位置越有利,能獲取太陽輻射也越多,對液流的驅(qū)動(dòng)作用也越強(qiáng)[23];另外,樹木優(yōu)勢度越大,峰值越大,反映了優(yōu)勢度大的樹木樹體水容和根系吸水能力較強(qiáng)[24]。

    青海云杉樹干液流季節(jié)動(dòng)態(tài)因樹木優(yōu)勢度等級而異。平均液流密度基本呈現(xiàn)6月份最大,其次是8月份,9、10月份明顯減小。這是由于青海云杉液流密度季節(jié)變化與葉面積指數(shù)的季節(jié)變化有很大關(guān)系[25];8月份青海云杉的葉面積指數(shù)較大且穩(wěn)定,介于3.0—4.0,液流密度較高;而9、10月份葉面積指數(shù)減小到2.5左右,青海云杉已進(jìn)入生長末期,液流密度明顯較低。優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木的季節(jié)變化趨勢基本一致,但被壓木8月份出現(xiàn)反彈。這可能是由于被壓木多處于主林層之下,林內(nèi)氣象因子對液流的影響較大,當(dāng)8月份林內(nèi)溫度較高時(shí),會促進(jìn)青海云杉體內(nèi)的代謝加快,從而促進(jìn)液流密度加快。而且,還發(fā)現(xiàn)液流密度月均值呈現(xiàn)出優(yōu)勢木>亞優(yōu)勢木>中等木>被壓木,表明樹木優(yōu)勢度越大其液流密度也越高,其原因可能與優(yōu)勢度大的樹木獲得的水、熱、營養(yǎng)等資源較多有關(guān)。

    此外,青海云杉日均液流密度還受林外氣象因子的影響,日均液流密度與太陽輻射強(qiáng)度、飽和水氣壓差和空氣溫度呈正相關(guān)關(guān)系,而與空氣相對濕度和降雨量呈負(fù)相關(guān)關(guān)系。這與銀白楊[7]、紅松[13]、冷杉[21]和華北落葉松[22]等研究結(jié)論基本一致。優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木日均液流密度均受這些林外氣象因子的影響,且影響程度依次減小。而且,影響的主導(dǎo)因子因樹木優(yōu)勢度等級也不同,優(yōu)勢木、亞優(yōu)勢木和中等木液流密度的主要受太陽輻射強(qiáng)度的影響,而被壓木液流密度主要受空氣相對濕度的影響。這是由于優(yōu)勢木、亞優(yōu)勢木和中等木獲取太陽輻射較為容易,且根系吸水范圍也大[14,24],而被壓木處于獲取光源不利位置,且根系吸水范圍有限,空氣相對濕度的適度增大,能有效地補(bǔ)充植物體內(nèi)的水分虧損,當(dāng)溫度升高時(shí),會促進(jìn)植物體內(nèi)的代謝加快,從而促進(jìn)液流密度加快。孫龍等[13]也發(fā)現(xiàn)影響紅松優(yōu)勢木、亞優(yōu)勢木和被壓木液流密度的主要環(huán)境因子為光合有效輻射和蒸汽壓虧缺,吳永波等[21]也發(fā)現(xiàn)影響冷杉液流密度的主要環(huán)境因子是光合有效輻射。

    5 結(jié)論

    本文通過對祁連山排露溝小流域青海云杉樹干液流密度研究,得到以下結(jié)論:

    (1)青海云杉樹干液流密度呈晝高夜低趨勢,晴天液流密度呈寬峰型曲線,變化幅度較大,而陰天和雨天液流呈窄峰型曲線,陰天液流曲線變化比較平緩,雨天液流曲線變化幅度相對較小。

    (2)晴天優(yōu)勢木、亞優(yōu)勢木、中等木和被壓木液流密度呈單峰型曲線,液流從7:00—8:00啟動(dòng),12:00—13:00達(dá)到峰值,20:00—23:00降到夜間水平。但樹木優(yōu)勢度越大,液流密度啟動(dòng)時(shí)間越早,結(jié)束時(shí)間越晚,液流時(shí)間跨度越長,峰值也越大。而且,優(yōu)勢木的平均液流密度最大,為(0.0758±0.0475) mL cm-2min-1,是亞優(yōu)勢木的1.5倍、是中等木和被壓木的1.68倍。但被壓木平均夜間液流密度最大,為(0.0419±0.0029) mL cm-2min-1,是中等木和優(yōu)勢木的1.5倍,是亞優(yōu)勢木的1.8倍。

    (3)青海云杉平均液流密度基本呈現(xiàn)6月份最大,其次是8月份,9、10月份明顯減小,但被壓木8月份稍有反彈。而且,液流密度月均值隨著優(yōu)勢度的減小而減小,液流密度月均值之間的差異隨著生長季的進(jìn)程而減小。

    (4)青海云杉日均液流密度與太陽輻射強(qiáng)度、飽和水氣壓差和空氣溫度呈正相關(guān)關(guān)系,而與空氣相對濕度和降雨量呈負(fù)相關(guān)關(guān)系。太陽輻射強(qiáng)度是影響優(yōu)勢木、亞優(yōu)勢木和中等木液流密度的主要?dú)庀笠蜃?而被壓木液流密度主要受空氣相對濕度的影響。

    [1] Grier C G, Running S W. Leaf area of mature Northwestern coniferous forests: relation to site water balance. Ecology, 1977, 58(4): 893- 899.

    [2] Hatton T J, Vertessy R A. Transpiration of plantationPinusradiataestimated by the heat pulse method and the Bowen Ratio. Hydrological Processes, 1990, 4(3): 289- 298.

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

    [4] 熊偉, 王彥輝, 于澎濤, 劉海龍, 徐麗宏, 時(shí)忠杰, 莫菲. 華北落葉松樹干液流的個(gè)體差異和林分蒸騰估計(jì)的尺度上推. 林業(yè)科學(xué), 2008, 44(1): 34- 40.

    [6] K?stner B, Granier A, Cermák J. Sapflow measurements in forest stands: methods and uncertainties. Annales Des Sciences Forestières, 1998, 55(1/2): 13- 27.

    [7] 張俊, 李曉飛, 李建貴, 王華, 黃春堂, 閔首軍, 李剛, 張福海, 田旭, 孔軍. 干旱荒漠區(qū)銀白楊樹干液流動(dòng)態(tài). 生態(tài)學(xué)報(bào), 2013, 33(18): 5655- 5660.

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

    [9] 李振華, 王彥輝, 于澎濤, 王艷兵, 韓新生, 曹恭祥. 華北落葉松液流速率的優(yōu)勢度差異及其對林分蒸騰估計(jì)的影響. 林業(yè)科學(xué)研究, 2015, 28(1): 8- 16.

    [10] 馮永建, 馬長明, 王彥輝, 杜阿朋. 華北落葉松人工林蒸騰特征及其與土壤水勢的關(guān)系. 中國水土保持科學(xué), 2010, 8(1): 93- 98.

    [11] Chang X X, Zhao W Z, He Z B. Radial pattern of sap flow and response to microclimate and soil moisture in Qinghai spruce (Piceacrassifolia) in the upper Heihe River Basin of arid northwestern China. Agricultural and Forest Meteorology, 2014, 187: 14- 21.

    [12] Granier A. Evaluation of transpiration in a Douglas-fir stand by means of sap flow measurements. Tree Physiology, 1987, 3(4): 309- 20.

    [13] 孫龍, 王傳寬, 楊國亭, 張全智, 周曉峰. 應(yīng)用熱擴(kuò)散技術(shù)對紅松人工林樹干液流通量的研究. 林業(yè)科學(xué), 2007, 43(11): 8- 14.

    [14] 梅婷婷, 王傳寬, 趙平, 蔡錫安, 劉曉靜, 張全智. 木荷樹干液流的密度特征. 林業(yè)科學(xué), 2010, 46(1): 40- 47.

    [15] 車克鈞, 傅輝恩, 王金葉. 祁連山水源林生態(tài)系統(tǒng)結(jié)構(gòu)與功能的研究. 林業(yè)科學(xué), 1998, 34(5): 29- 37.

    [16] Chang X X, Zhao W Z, Liu H, Wei X, Liu B, He Z B. Qinghai spruce (Piceacrassifolia) forest transpiration and canopy conductance in the upper Heihe River Basin of arid northwestern China. Agricultural and Forest Meteorology, 2014, 198- 199: 209- 220.

    [17] 魏瀟, 常學(xué)向, 楊淇越, 劉冰, 王國華. 祁連山青海云杉(Piceacrassifolia)夜間樹干液流特征及影響因素. 冰川凍土, 2015, 37(1): 87- 94.

    [18] 德國林業(yè)和林產(chǎn)品研究中心, ICP項(xiàng)目協(xié)調(diào)中心編輯. 空氣污染對森林影響的統(tǒng)一采樣、評價(jià)、監(jiān)測和分析的方法與標(biāo)準(zhǔn)手冊. 北京: 中國科學(xué)技術(shù)出版社, 2002.

    [19] 常學(xué)向, 趙文智, 何志斌, 王金鳳. 青海云杉(Piceacrassifolia)邊材心材邊界的確定及樹干傳輸水分的空間格局. 冰川凍土, 2013, 35(2): 483- 489.

    [20] 馬玲, 趙平, 饒興權(quán),蔡錫安, 曾小平. 喬木蒸騰作用的主要測定方法. 生態(tài)學(xué)雜志, 2005, 24(1): 88- 96.

    [21] 吳永波, 薛建輝. 岷江流域冷杉樹干液流的動(dòng)態(tài)變化規(guī)律. 南京林業(yè)大學(xué)學(xué)報(bào): 自然科學(xué)版, 2005, 29(6): 61- 64.

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

    [23] 趙平, 饒興權(quán), 馬玲, 蔡錫安, 曾小平. 馬占相思(Acaciamangium)樹干液流密度和整樹蒸騰的個(gè)體差異. 生態(tài)學(xué)報(bào), 2006, 26(12): 4050- 4058.

    [24] 馬履一, 王華田, 林平. 北京地區(qū)幾個(gè)造林樹種耗水性比較研究. 北京林業(yè)大學(xué)學(xué)報(bào), 2003, 25(2): 1- 7.

    [25] 常學(xué)向, 趙文智. 荒漠綠洲農(nóng)田防護(hù)樹種二白楊生長季節(jié)樹干液流的變化. 生態(tài)學(xué)報(bào), 2004, 24(7): 1436- 1441.

    Variation in sap flow density among levels of tree dominance inPiceacrassifoliain the Qilian Mountains

    WAN Yanfang1, 2, YU Pengtao2, *, LIU Xiande1, 3, WANG Shunli3, WANG Yanhui2, XIONG Wei2

    1ForestryofGansuAgricultureUniversity,Lanzhou730070,China2ResearchInstituteofForestEcology,EnvironmentandProtection,ChineseAcademyofForestry,Beijing100091,China3AcademyofWaterResourcesConservationForestsinQilianMountainofGansuProvince,Zhangye734000,China

    Tree transpiration plays a determining role in water balance for forest stands and in water yield from forested catchments. In the present study, an experiment was conducted in the Pailugou watershed in the Qilian Mountains, in the arid region of Northwest China. In a 86-year-oldPiceacrassifoliaforest stand, 3—5 trees from each of dominant, subdominant, intermediate and suppressed trees were chosen as sample trees. The sap flux density for these trees was measured using the thermal dissipation probe (TDP) method from June 16 to October 14, 2015. Furthermore, the related meteorological factors, including solar radiation (Rs), air temperature (T), air relative humidity (Rh), and precipitation (P) were simultaneously monitored by an automatic meteorological station in the outer forest. The results showed that: (1) Daily sap flow density was higher during the day than during the night. On sunny days, daily variation in sap flow density exhibited large amplitudes, whereas it exhibited little variation on cloudy and rainy days. (2) On sunny days, daily sap flow density began earlier in the morning, ended later at night, and exhibited a larger maximum of sap flow density as tree dominance class increased. The value of the mean sap flow density of dominant trees was (0.0758 ± 0.0475) mL cm-2min-1, which was 1.5 times of that of subdominant trees and 1.68 times that of intermediate and suppressed trees. (3) Mean daily sap flow density for all dominant trees showed that the maximum values were achieved in June, with the next highest values in August, and decreasing values in September and October. In addition, sap flow density exhibited the following order: dominant > subdominant > intermediate > suppressed trees. (4) Correlation and stepwise regression analyses indicated that mean daily sap flow density was significantly correlated with solar radiation, vapor pressure deficit, and air temperature, and was negatively correlated with air relative humidity and precipitation. Solar radiation was the main meteorological factor influencing mean daily sap flow density of dominant, subdominant, and intermediate trees, whereas air relative humidity primarily influenced suppressed trees.

    sap flow; dominance; meteorological factors;Piceacrassifoliaforest; Qilian Mountains

    國家自然科學(xué)基金項(xiàng)目(91225302,91425301,31360201,41461004);甘肅省基礎(chǔ)研究創(chuàng)新群體(145RJIG337)

    2016- 01- 28; 網(wǎng)絡(luò)出版日期:2016- 12- 19

    10.5846/stxb201601280203

    *通訊作者Corresponding author.E-mail: yupt@caf.ac.cn

    萬艷芳, 于澎濤,劉賢德,王順利,王彥輝,熊偉.祁連山青海云杉樹干液流密度的優(yōu)勢度差異.生態(tài)學(xué)報(bào),2017,37(9):3106- 3114.

    Wan Y F, Yu P T, Liu X D, Wang S L, Wang Y H, Xiong W.Variation in sap flow density among levels of tree dominance inPiceacrassifoliain the Qilian Mountains.Acta Ecologica Sinica,2017,37(9):3106- 3114.

    猜你喜歡
    液流云杉晴天
    想和白云握手的樹
    它若安好,便是晴天
    一起來寫“雪”
    晴天有時(shí)下豬
    云 杉
    不同云杉品種幼林苗高的遺傳變異研究
    一種全釩液流電池并網(wǎng)控制系統(tǒng)的仿真設(shè)計(jì)
    西秦嶺山地云杉育苗技術(shù)
    H2-Fe3+/Fe2+氧化還原液流電池初步研究
    我國—項(xiàng)液流電池國際標(biāo)準(zhǔn)提案獲準(zhǔn)通過
    热99re8久久精品国产| 一区二区三区高清视频在线| 国产精品,欧美在线| 少妇裸体淫交视频免费看高清 | 极品教师在线免费播放| 欧美色视频一区免费| 午夜福利成人在线免费观看| 搞女人的毛片| 国产成人系列免费观看| 亚洲七黄色美女视频| 麻豆国产av国片精品| 高清在线国产一区| 一级黄色大片毛片| 午夜福利一区二区在线看| 国产亚洲精品久久久久5区| 可以在线观看的亚洲视频| 午夜久久久在线观看| 美国免费a级毛片| 午夜福利高清视频| 国产欧美日韩精品亚洲av| 一边摸一边抽搐一进一小说| 久久草成人影院| 国产在线精品亚洲第一网站| 欧美三级亚洲精品| 波多野结衣高清作品| 男女视频在线观看网站免费 | 真人做人爱边吃奶动态| 国产成人欧美| 国产精品二区激情视频| 啪啪无遮挡十八禁网站| 国内久久婷婷六月综合欲色啪| 亚洲自偷自拍图片 自拍| 久久久久久人人人人人| 欧美黑人欧美精品刺激| 午夜免费激情av| 校园春色视频在线观看| 成年版毛片免费区| 黄色女人牲交| 这个男人来自地球电影免费观看| 色综合欧美亚洲国产小说| 午夜日韩欧美国产| 欧美+亚洲+日韩+国产| 天堂√8在线中文| 一区二区三区精品91| www.精华液| 在线十欧美十亚洲十日本专区| 欧美又色又爽又黄视频| 国产一区二区在线av高清观看| 久久性视频一级片| 一级毛片女人18水好多| 999精品在线视频| 国产一卡二卡三卡精品| 女生性感内裤真人,穿戴方法视频| or卡值多少钱| 精品久久久久久久毛片微露脸| 色综合亚洲欧美另类图片| 在线观看一区二区三区| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 欧洲精品卡2卡3卡4卡5卡区| 一边摸一边做爽爽视频免费| 99精品欧美一区二区三区四区| 1024视频免费在线观看| 淫妇啪啪啪对白视频| 久久精品91蜜桃| 宅男免费午夜| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 波多野结衣高清无吗| 变态另类丝袜制服| av福利片在线| 91九色精品人成在线观看| 婷婷精品国产亚洲av在线| videosex国产| 国产精品亚洲av一区麻豆| 国产熟女xx| 日韩精品青青久久久久久| 啦啦啦免费观看视频1| 国产又色又爽无遮挡免费看| 久久婷婷人人爽人人干人人爱| 国产高清激情床上av| 美女免费视频网站| 亚洲激情在线av| 国产男靠女视频免费网站| 日本成人三级电影网站| 一二三四在线观看免费中文在| 757午夜福利合集在线观看| 欧美日韩瑟瑟在线播放| 黑人巨大精品欧美一区二区mp4| 日韩精品中文字幕看吧| 成人国产综合亚洲| 侵犯人妻中文字幕一二三四区| 亚洲成人免费电影在线观看| 午夜福利高清视频| 色在线成人网| 国产一卡二卡三卡精品| or卡值多少钱| 在线观看午夜福利视频| 欧美丝袜亚洲另类 | 国产真实乱freesex| 精品久久久久久久久久免费视频| 亚洲中文字幕一区二区三区有码在线看 | 国产精品一区二区免费欧美| 精品国产超薄肉色丝袜足j| 亚洲精品色激情综合| 黄片大片在线免费观看| 色综合婷婷激情| 午夜免费成人在线视频| 国产精品久久久久久人妻精品电影| 99热这里只有精品一区 | 成熟少妇高潮喷水视频| 亚洲熟女毛片儿| 国产成年人精品一区二区| 欧美色视频一区免费| 国产激情偷乱视频一区二区| 国产精品久久久av美女十八| 欧美zozozo另类| 亚洲男人的天堂狠狠| 久热爱精品视频在线9| 国产精品久久久人人做人人爽| 国产精品 国内视频| 欧美国产日韩亚洲一区| 12—13女人毛片做爰片一| 亚洲精品国产精品久久久不卡| 黄色片一级片一级黄色片| 亚洲av中文字字幕乱码综合 | 美女扒开内裤让男人捅视频| 精品国产超薄肉色丝袜足j| 久久久水蜜桃国产精品网| 在线观看日韩欧美| 超碰成人久久| 久久中文看片网| 午夜精品在线福利| av电影中文网址| 国产av在哪里看| 免费女性裸体啪啪无遮挡网站| 欧美成人免费av一区二区三区| 精品福利观看| 亚洲国产欧洲综合997久久, | www日本黄色视频网| 久久亚洲真实| 日日摸夜夜添夜夜添小说| www.自偷自拍.com| 久久天躁狠狠躁夜夜2o2o| xxx96com| 天天躁夜夜躁狠狠躁躁| 亚洲国产精品成人综合色| 99热这里只有精品一区 | 国产真人三级小视频在线观看| 日韩精品免费视频一区二区三区| 亚洲熟妇熟女久久| 免费看美女性在线毛片视频| 久久中文字幕一级| 国产一区二区在线av高清观看| 90打野战视频偷拍视频| 国产精品自产拍在线观看55亚洲| 欧美成人一区二区免费高清观看 | 国产私拍福利视频在线观看| 亚洲欧美一区二区三区黑人| 丝袜人妻中文字幕| 日韩有码中文字幕| 日本撒尿小便嘘嘘汇集6| 精品国产乱子伦一区二区三区| 超碰成人久久| 午夜激情av网站| 老司机靠b影院| 在线观看www视频免费| 亚洲中文字幕日韩| 国产精品综合久久久久久久免费| 日韩免费av在线播放| 欧美+亚洲+日韩+国产| 午夜福利18| 午夜免费观看网址| 国产av又大| 久久 成人 亚洲| 动漫黄色视频在线观看| 久久狼人影院| 国产精品爽爽va在线观看网站 | 国产伦一二天堂av在线观看| 岛国视频午夜一区免费看| 无人区码免费观看不卡| 国产精品免费视频内射| 欧美性猛交╳xxx乱大交人| 国产激情久久老熟女| 国产午夜福利久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久电影中文字幕| 亚洲全国av大片| 日日干狠狠操夜夜爽| 一本精品99久久精品77| 国产三级在线视频| 国产色视频综合| 91在线观看av| 中文字幕人妻丝袜一区二区| 精品久久蜜臀av无| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 亚洲国产精品久久男人天堂| 天天躁夜夜躁狠狠躁躁| 精品久久久久久久末码| 中文资源天堂在线| 免费在线观看完整版高清| 午夜老司机福利片| 又紧又爽又黄一区二区| 欧美性猛交黑人性爽| 国产精品久久电影中文字幕| 日韩欧美国产一区二区入口| 视频在线观看一区二区三区| 亚洲一区二区三区不卡视频| 亚洲精品美女久久久久99蜜臀| 久久久久久亚洲精品国产蜜桃av| 999精品在线视频| 黄频高清免费视频| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av香蕉五月| 最近在线观看免费完整版| 一本综合久久免费| 亚洲av成人av| www国产在线视频色| 亚洲国产高清在线一区二区三 | 国产成人av激情在线播放| 欧美日本亚洲视频在线播放| 一本精品99久久精品77| 免费在线观看影片大全网站| 啪啪无遮挡十八禁网站| 香蕉久久夜色| 日韩三级视频一区二区三区| 亚洲男人的天堂狠狠| ponron亚洲| 母亲3免费完整高清在线观看| 一级毛片女人18水好多| 国产蜜桃级精品一区二区三区| 麻豆国产av国片精品| 亚洲最大成人中文| 欧美日韩中文字幕国产精品一区二区三区| 久久人妻av系列| av天堂在线播放| 18禁裸乳无遮挡免费网站照片 | 长腿黑丝高跟| 丰满人妻熟妇乱又伦精品不卡| 久久精品人妻少妇| 久久午夜亚洲精品久久| 国产精品免费视频内射| 久久中文字幕人妻熟女| 日日夜夜操网爽| 一区福利在线观看| 美女 人体艺术 gogo| 国语自产精品视频在线第100页| 成人三级黄色视频| 国产精品一区二区三区四区久久 | 99久久综合精品五月天人人| 搞女人的毛片| 婷婷精品国产亚洲av在线| 国产精品,欧美在线| 88av欧美| 俄罗斯特黄特色一大片| 制服诱惑二区| 久久热在线av| 国产精品亚洲美女久久久| 97碰自拍视频| 免费高清在线观看日韩| 国产极品粉嫩免费观看在线| 国产区一区二久久| 欧美大码av| 亚洲一区中文字幕在线| 日韩 欧美 亚洲 中文字幕| 欧美在线一区亚洲| 日韩欧美在线二视频| 日本成人三级电影网站| 美女高潮喷水抽搐中文字幕| 一进一出抽搐动态| 色播在线永久视频| 人人妻人人澡欧美一区二区| 亚洲七黄色美女视频| 午夜成年电影在线免费观看| 欧美成人午夜精品| 亚洲精品美女久久久久99蜜臀| 一进一出抽搐动态| 在线观看舔阴道视频| 美女免费视频网站| 国产精品自产拍在线观看55亚洲| 亚洲成av人片免费观看| 50天的宝宝边吃奶边哭怎么回事| 日本成人三级电影网站| 亚洲久久久国产精品| 日韩大码丰满熟妇| 一本精品99久久精品77| 在线观看免费视频日本深夜| 香蕉av资源在线| 欧美成狂野欧美在线观看| 国产成人欧美在线观看| 欧美日本视频| 亚洲熟女毛片儿| 91成人精品电影| 欧美黄色片欧美黄色片| 久久中文字幕人妻熟女| 久久精品国产99精品国产亚洲性色| 成人精品一区二区免费| 国内少妇人妻偷人精品xxx网站 | 久久草成人影院| 一区二区三区国产精品乱码| 自线自在国产av| 美女高潮到喷水免费观看| 搡老妇女老女人老熟妇| 可以在线观看的亚洲视频| 国产1区2区3区精品| 亚洲精品国产区一区二| 日韩欧美三级三区| cao死你这个sao货| 精品久久久久久久久久久久久 | 国产精品,欧美在线| 国产av又大| 国产亚洲精品综合一区在线观看 | 日韩欧美一区二区三区在线观看| 成年人黄色毛片网站| 波多野结衣巨乳人妻| 丁香欧美五月| a在线观看视频网站| 亚洲中文日韩欧美视频| 中文字幕人成人乱码亚洲影| 亚洲五月婷婷丁香| 黄色毛片三级朝国网站| 欧美日本亚洲视频在线播放| 制服丝袜大香蕉在线| 两个人看的免费小视频| 久久天堂一区二区三区四区| 日本在线视频免费播放| 国内揄拍国产精品人妻在线 | 国产成人欧美| 欧美色视频一区免费| 中文亚洲av片在线观看爽| 一个人观看的视频www高清免费观看 | av在线天堂中文字幕| 国产精品野战在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲欧美日韩高清在线视频| 亚洲久久久国产精品| 中文字幕最新亚洲高清| 在线永久观看黄色视频| 亚洲avbb在线观看| 亚洲片人在线观看| 久久伊人香网站| 搡老岳熟女国产| 精品一区二区三区四区五区乱码| 韩国精品一区二区三区| 亚洲成人久久爱视频| 国产乱人伦免费视频| 99在线人妻在线中文字幕| 九色国产91popny在线| 免费av毛片视频| 亚洲国产中文字幕在线视频| 国产av一区二区精品久久| 欧美另类亚洲清纯唯美| 一区二区日韩欧美中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 久久99热这里只有精品18| 美女扒开内裤让男人捅视频| 国产主播在线观看一区二区| 国产av一区在线观看免费| 久热这里只有精品99| 精品久久久久久成人av| 国产蜜桃级精品一区二区三区| 国产精品久久电影中文字幕| 久久久国产成人免费| 精品一区二区三区av网在线观看| 久久久久久人人人人人| 人成视频在线观看免费观看| 怎么达到女性高潮| 日韩精品青青久久久久久| 99久久无色码亚洲精品果冻| 中文在线观看免费www的网站 | 久久久久亚洲av毛片大全| 最近最新免费中文字幕在线| 热re99久久国产66热| 日韩欧美一区二区三区在线观看| 久久久国产成人免费| 在线十欧美十亚洲十日本专区| 99热只有精品国产| 首页视频小说图片口味搜索| 久久精品91蜜桃| 男女之事视频高清在线观看| 50天的宝宝边吃奶边哭怎么回事| 99re在线观看精品视频| 日韩有码中文字幕| 黄色成人免费大全| 亚洲人成77777在线视频| 欧美国产日韩亚洲一区| 欧美日本亚洲视频在线播放| 国产亚洲欧美98| 午夜福利18| www国产在线视频色| 97人妻精品一区二区三区麻豆 | 好男人电影高清在线观看| 女警被强在线播放| 国产亚洲精品第一综合不卡| 精品日产1卡2卡| 美女免费视频网站| 亚洲国产精品久久男人天堂| 免费电影在线观看免费观看| 精品国内亚洲2022精品成人| 90打野战视频偷拍视频| 中文字幕精品免费在线观看视频| 亚洲黑人精品在线| 亚洲中文字幕日韩| 精品熟女少妇八av免费久了| 亚洲精品一区av在线观看| 日韩一卡2卡3卡4卡2021年| 欧美又色又爽又黄视频| 免费看日本二区| 国产不卡一卡二| 午夜成年电影在线免费观看| 黄片小视频在线播放| 我的亚洲天堂| 亚洲无线在线观看| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 黄网站色视频无遮挡免费观看| 午夜精品在线福利| 亚洲av成人不卡在线观看播放网| 首页视频小说图片口味搜索| 国产高清激情床上av| 日韩欧美国产在线观看| 久久久水蜜桃国产精品网| 精品久久久久久久人妻蜜臀av| 亚洲人成77777在线视频| 国产伦在线观看视频一区| 狠狠狠狠99中文字幕| 老司机在亚洲福利影院| 一区二区三区精品91| 老司机靠b影院| 黄色片一级片一级黄色片| 久久欧美精品欧美久久欧美| 成年免费大片在线观看| 美女午夜性视频免费| 国产激情久久老熟女| 男人舔女人下体高潮全视频| 成人18禁在线播放| 岛国视频午夜一区免费看| 日本一本二区三区精品| 亚洲精品国产精品久久久不卡| 国产亚洲精品综合一区在线观看 | 国产亚洲欧美98| www日本黄色视频网| 久久伊人香网站| 男人舔女人的私密视频| 非洲黑人性xxxx精品又粗又长| 久久久国产欧美日韩av| 亚洲av电影不卡..在线观看| 宅男免费午夜| 一级毛片高清免费大全| 亚洲七黄色美女视频| 一二三四在线观看免费中文在| 亚洲黑人精品在线| 亚洲va日本ⅴa欧美va伊人久久| 国内毛片毛片毛片毛片毛片| 99久久99久久久精品蜜桃| 久久久久久国产a免费观看| 欧美最黄视频在线播放免费| 91麻豆精品激情在线观看国产| 国产成人啪精品午夜网站| 免费电影在线观看免费观看| 一级毛片高清免费大全| 两个人视频免费观看高清| 又紧又爽又黄一区二区| 操出白浆在线播放| 亚洲国产精品sss在线观看| 久久伊人香网站| 中文字幕人妻丝袜一区二区| 在线观看舔阴道视频| 欧美在线一区亚洲| 国产一区二区三区视频了| 久久这里只有精品19| 欧美在线黄色| 久久热在线av| 亚洲狠狠婷婷综合久久图片| 亚洲色图av天堂| 欧美成狂野欧美在线观看| 波多野结衣高清无吗| 国产成人欧美在线观看| 宅男免费午夜| 99久久99久久久精品蜜桃| 亚洲一码二码三码区别大吗| 中文亚洲av片在线观看爽| 一进一出抽搐动态| 美女免费视频网站| 两个人看的免费小视频| 日韩视频一区二区在线观看| 免费av毛片视频| 看免费av毛片| 精品免费久久久久久久清纯| 女同久久另类99精品国产91| 无人区码免费观看不卡| av欧美777| 日韩大尺度精品在线看网址| 久久中文字幕人妻熟女| 久热爱精品视频在线9| 亚洲成av人片免费观看| 亚洲美女黄片视频| 成人亚洲精品av一区二区| 岛国视频午夜一区免费看| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 亚洲 国产 在线| 人妻久久中文字幕网| 中文字幕人妻丝袜一区二区| 日日夜夜操网爽| 国产一区二区激情短视频| 中文字幕高清在线视频| 99国产综合亚洲精品| 一本综合久久免费| 成人午夜高清在线视频 | 少妇裸体淫交视频免费看高清 | 伦理电影免费视频| 国产免费av片在线观看野外av| 亚洲avbb在线观看| 久久国产精品男人的天堂亚洲| 国产区一区二久久| 此物有八面人人有两片| 中文字幕人妻熟女乱码| 12—13女人毛片做爰片一| 午夜精品在线福利| 亚洲精品在线美女| 中国美女看黄片| 侵犯人妻中文字幕一二三四区| 男女视频在线观看网站免费 | 中文在线观看免费www的网站 | 亚洲第一青青草原| 两人在一起打扑克的视频| 欧美成人午夜精品| 中国美女看黄片| 成年版毛片免费区| or卡值多少钱| 久久久久久久久中文| 母亲3免费完整高清在线观看| 亚洲三区欧美一区| 欧美乱妇无乱码| 亚洲精品中文字幕在线视频| 国产成人av教育| 国产男靠女视频免费网站| 国产精品亚洲av一区麻豆| 精华霜和精华液先用哪个| 成人18禁在线播放| 亚洲第一青青草原| 成人一区二区视频在线观看| 国产一区二区激情短视频| 色尼玛亚洲综合影院| 国产蜜桃级精品一区二区三区| 成人18禁在线播放| 免费av毛片视频| 国产久久久一区二区三区| av在线播放免费不卡| 成人亚洲精品av一区二区| 男人舔女人下体高潮全视频| 一卡2卡三卡四卡精品乱码亚洲| 精品国产一区二区三区四区第35| 亚洲第一电影网av| 亚洲精品美女久久久久99蜜臀| 窝窝影院91人妻| 久久久久久久午夜电影| 脱女人内裤的视频| 亚洲中文av在线| 中文字幕高清在线视频| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品综合一区在线观看 | 亚洲欧美一区二区三区黑人| 天天一区二区日本电影三级| 欧美zozozo另类| 国产精品免费视频内射| 久久精品成人免费网站| 国产视频一区二区在线看| 最近在线观看免费完整版| 午夜福利欧美成人| 久久 成人 亚洲| 久久天躁狠狠躁夜夜2o2o| 日本五十路高清| 国产aⅴ精品一区二区三区波| 免费观看人在逋| 老司机靠b影院| 美女扒开内裤让男人捅视频| 久久久国产精品麻豆| 国产精品亚洲一级av第二区| 女人被狂操c到高潮| 免费高清视频大片| 又黄又爽又免费观看的视频| 麻豆成人午夜福利视频| 久久精品国产亚洲av香蕉五月| 国产亚洲精品av在线| 给我免费播放毛片高清在线观看| 欧美在线黄色| 精品高清国产在线一区| a在线观看视频网站| 国产成人系列免费观看| 窝窝影院91人妻| 午夜福利在线观看吧| 午夜福利18| 亚洲 欧美一区二区三区| 男女视频在线观看网站免费 | 日韩视频一区二区在线观看| 精品国产亚洲在线| 欧美日韩亚洲国产一区二区在线观看| videosex国产| 日本黄色视频三级网站网址| 人妻久久中文字幕网| 国产精品爽爽va在线观看网站 | 免费人成视频x8x8入口观看| 此物有八面人人有两片| 久久热在线av| cao死你这个sao货| av天堂在线播放| 怎么达到女性高潮| 99久久无色码亚洲精品果冻| 日本五十路高清| 一级毛片女人18水好多| 国产精品久久久久久人妻精品电影| 午夜福利在线在线|