• <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)通過
    国产精品一区二区免费欧美| 男女下面进入的视频免费午夜 | 黑人巨大精品欧美一区二区mp4| 99久久人妻综合| 丰满的人妻完整版| 久久精品成人免费网站| 亚洲国产中文字幕在线视频| 国产精品久久久久成人av| 丝袜在线中文字幕| 久久精品aⅴ一区二区三区四区| 一个人免费在线观看的高清视频| 99热国产这里只有精品6| 黑人欧美特级aaaaaa片| 久久久久久免费高清国产稀缺| 久久精品aⅴ一区二区三区四区| 亚洲欧美一区二区三区黑人| 国产成人av教育| 午夜免费鲁丝| 热re99久久精品国产66热6| 国产av一区二区精品久久| 亚洲欧美一区二区三区黑人| 免费久久久久久久精品成人欧美视频| 日韩国内少妇激情av| 久久人人97超碰香蕉20202| 一二三四社区在线视频社区8| 久久久久国内视频| 国产高清国产精品国产三级| 丰满饥渴人妻一区二区三| 久久精品aⅴ一区二区三区四区| 在线看a的网站| 欧美丝袜亚洲另类 | av中文乱码字幕在线| 欧美激情极品国产一区二区三区| av福利片在线| 热re99久久国产66热| 亚洲成人久久性| 免费av毛片视频| 一级片免费观看大全| 大码成人一级视频| www日本在线高清视频| 91字幕亚洲| av欧美777| 久久精品国产清高在天天线| 亚洲av日韩精品久久久久久密| 18禁美女被吸乳视频| 国产成人精品久久二区二区91| 国产精品偷伦视频观看了| 久久这里只有精品19| 夜夜躁狠狠躁天天躁| 激情视频va一区二区三区| 欧美午夜高清在线| 一区二区三区精品91| 伦理电影免费视频| 国产无遮挡羞羞视频在线观看| 黄片小视频在线播放| av超薄肉色丝袜交足视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品免费视频内射| 亚洲,欧美精品.| 一级片免费观看大全| 多毛熟女@视频| 国产乱人伦免费视频| av欧美777| 一级片免费观看大全| 多毛熟女@视频| 国产精品免费视频内射| 男人舔女人下体高潮全视频| 母亲3免费完整高清在线观看| 精品欧美一区二区三区在线| 亚洲av五月六月丁香网| 1024香蕉在线观看| 青草久久国产| 久久人人爽av亚洲精品天堂| 一进一出抽搐gif免费好疼 | 又黄又粗又硬又大视频| 亚洲第一av免费看| 又大又爽又粗| 亚洲精品久久午夜乱码| 成人av一区二区三区在线看| 国产亚洲欧美在线一区二区| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 91字幕亚洲| 亚洲av电影在线进入| 国产精品秋霞免费鲁丝片| 成人精品一区二区免费| 国产在线观看jvid| 亚洲五月天丁香| 亚洲精品国产色婷婷电影| 亚洲五月色婷婷综合| 在线永久观看黄色视频| 成人手机av| 91老司机精品| 国产精品影院久久| 欧美在线一区亚洲| 两个人看的免费小视频| 91字幕亚洲| 久久久精品欧美日韩精品| 亚洲人成伊人成综合网2020| 亚洲av五月六月丁香网| 亚洲精品美女久久av网站| 18禁黄网站禁片午夜丰满| 黑人操中国人逼视频| 亚洲成人免费av在线播放| 亚洲国产精品一区二区三区在线| 91大片在线观看| 热re99久久精品国产66热6| 99在线人妻在线中文字幕| 日韩有码中文字幕| 高清欧美精品videossex| 欧美av亚洲av综合av国产av| 一区二区三区精品91| 精品国内亚洲2022精品成人| 色在线成人网| 久久久久国内视频| 久久久久久久久中文| 亚洲少妇的诱惑av| 亚洲五月天丁香| 精品久久久久久久久久免费视频 | 男人的好看免费观看在线视频 | 国产免费现黄频在线看| 香蕉国产在线看| 精品一区二区三区av网在线观看| 天天躁夜夜躁狠狠躁躁| 亚洲精品中文字幕一二三四区| 久9热在线精品视频| 多毛熟女@视频| 777久久人妻少妇嫩草av网站| 久久人人97超碰香蕉20202| 国产精品自产拍在线观看55亚洲| 国产成人影院久久av| 99久久国产精品久久久| 成人免费观看视频高清| 麻豆av在线久日| 99精品久久久久人妻精品| 亚洲欧洲精品一区二区精品久久久| av电影中文网址| 亚洲精品一区av在线观看| 亚洲av电影在线进入| 久久久国产成人免费| 久久影院123| 精品福利观看| 丁香欧美五月| 男女之事视频高清在线观看| 十分钟在线观看高清视频www| 最新美女视频免费是黄的| 九色亚洲精品在线播放| 在线视频色国产色| 久久久国产精品麻豆| 欧美人与性动交α欧美精品济南到| av天堂久久9| 99国产精品一区二区三区| 新久久久久国产一级毛片| 久久香蕉国产精品| 久久久久九九精品影院| 丰满人妻熟妇乱又伦精品不卡| 不卡av一区二区三区| www.自偷自拍.com| 精品福利永久在线观看| 亚洲中文日韩欧美视频| 亚洲成人国产一区在线观看| 在线av久久热| 国产无遮挡羞羞视频在线观看| 露出奶头的视频| 婷婷丁香在线五月| 大型黄色视频在线免费观看| 中亚洲国语对白在线视频| 50天的宝宝边吃奶边哭怎么回事| 国产又爽黄色视频| 两性夫妻黄色片| 国产熟女xx| 国产av一区二区精品久久| 99精品久久久久人妻精品| 老司机深夜福利视频在线观看| 91麻豆精品激情在线观看国产 | 一本综合久久免费| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看| 国产深夜福利视频在线观看| 十八禁人妻一区二区| 91精品国产国语对白视频| 国产又爽黄色视频| 如日韩欧美国产精品一区二区三区| 欧美日韩国产mv在线观看视频| 在线观看一区二区三区| 在线观看免费午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 国产真人三级小视频在线观看| 高清毛片免费观看视频网站 | 亚洲,欧美精品.| 国产麻豆69| 亚洲午夜精品一区,二区,三区| 亚洲成av片中文字幕在线观看| 国产97色在线日韩免费| 久久香蕉精品热| 一个人免费在线观看的高清视频| 久久人妻av系列| 午夜精品国产一区二区电影| videosex国产| 久久这里只有精品19| 国产精品电影一区二区三区| 999久久久国产精品视频| 久久久久国产精品人妻aⅴ院| 亚洲中文字幕日韩| 久久中文看片网| 亚洲欧美一区二区三区黑人| 91麻豆av在线| 一级片'在线观看视频| 精品欧美一区二区三区在线| 丝袜人妻中文字幕| 女人高潮潮喷娇喘18禁视频| 老司机靠b影院| 亚洲少妇的诱惑av| 亚洲国产毛片av蜜桃av| 国产av一区二区精品久久| 久久精品影院6| 中文字幕人妻丝袜制服| 亚洲九九香蕉| 久久久久久久久久久久大奶| 亚洲aⅴ乱码一区二区在线播放 | 日本 av在线| 欧美在线一区亚洲| 国产精品一区二区在线不卡| 美女国产高潮福利片在线看| 男人舔女人下体高潮全视频| 久久精品成人免费网站| 50天的宝宝边吃奶边哭怎么回事| 老熟妇仑乱视频hdxx| 视频区欧美日本亚洲| 老司机午夜福利在线观看视频| 在线国产一区二区在线| 日韩欧美国产一区二区入口| 亚洲久久久国产精品| 国产成人影院久久av| 免费av中文字幕在线| 女性被躁到高潮视频| av中文乱码字幕在线| 国产三级黄色录像| 久9热在线精品视频| 精品国内亚洲2022精品成人| 国产无遮挡羞羞视频在线观看| 久久天堂一区二区三区四区| 欧美在线一区亚洲| 操美女的视频在线观看| 亚洲av片天天在线观看| 久久人人97超碰香蕉20202| 久久久久久久午夜电影 | 精品福利观看| 亚洲片人在线观看| 久久久久久久精品吃奶| 欧美丝袜亚洲另类 | 亚洲av片天天在线观看| 九色亚洲精品在线播放| 在线观看一区二区三区| 最近最新中文字幕大全免费视频| xxxhd国产人妻xxx| 18禁国产床啪视频网站| 成熟少妇高潮喷水视频| 91字幕亚洲| 老鸭窝网址在线观看| 日日爽夜夜爽网站| 日韩欧美三级三区| 免费久久久久久久精品成人欧美视频| 日韩成人在线观看一区二区三区| 99久久国产精品久久久| 免费在线观看影片大全网站| 亚洲黑人精品在线| 成人精品一区二区免费| 99国产精品免费福利视频| 久久草成人影院| 不卡一级毛片| 18禁裸乳无遮挡免费网站照片 | 美女 人体艺术 gogo| 757午夜福利合集在线观看| 免费高清视频大片| svipshipincom国产片| 欧美+亚洲+日韩+国产| 色精品久久人妻99蜜桃| 午夜福利影视在线免费观看| 亚洲熟妇中文字幕五十中出 | 国产精品免费一区二区三区在线| 欧美成人免费av一区二区三区| 身体一侧抽搐| 午夜免费激情av| 亚洲一区二区三区不卡视频| 啦啦啦 在线观看视频| 日本撒尿小便嘘嘘汇集6| 高潮久久久久久久久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品中文字幕在线视频| 天堂俺去俺来也www色官网| 一边摸一边抽搐一进一出视频| 欧美国产精品va在线观看不卡| 免费看a级黄色片| 国内毛片毛片毛片毛片毛片| 亚洲avbb在线观看| 久久国产精品影院| 精品人妻在线不人妻| 电影成人av| 精品久久久久久成人av| 欧美成人午夜精品| 91精品三级在线观看| 别揉我奶头~嗯~啊~动态视频| 可以在线观看毛片的网站| 亚洲专区中文字幕在线| 欧美成人性av电影在线观看| 久久人妻av系列| 12—13女人毛片做爰片一| 亚洲中文av在线| 色婷婷av一区二区三区视频| 亚洲色图av天堂| 他把我摸到了高潮在线观看| 视频区图区小说| 黑人操中国人逼视频| 亚洲第一青青草原| 日韩欧美在线二视频| 欧美日本亚洲视频在线播放| 夜夜夜夜夜久久久久| 国产一区二区三区视频了| 交换朋友夫妻互换小说| 国产精品久久久久成人av| 午夜福利免费观看在线| 窝窝影院91人妻| 欧美黑人精品巨大| 99久久综合精品五月天人人| 久久久国产欧美日韩av| 日本vs欧美在线观看视频| 一进一出抽搐gif免费好疼 | 淫秽高清视频在线观看| 免费久久久久久久精品成人欧美视频| 欧美在线黄色| 很黄的视频免费| 国产无遮挡羞羞视频在线观看| 欧美成人性av电影在线观看| 欧美日韩国产mv在线观看视频| 19禁男女啪啪无遮挡网站| 天堂影院成人在线观看| 一区福利在线观看| 精品午夜福利视频在线观看一区| 亚洲国产精品999在线| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜两性在线视频| 18美女黄网站色大片免费观看| 免费在线观看视频国产中文字幕亚洲| 亚洲七黄色美女视频| 亚洲午夜理论影院| 国产精品日韩av在线免费观看 | 久久精品亚洲精品国产色婷小说| 欧美在线黄色| 99精品欧美一区二区三区四区| 欧美日本中文国产一区发布| 99久久国产精品久久久| 高清av免费在线| 日韩中文字幕欧美一区二区| 88av欧美| 国产成人一区二区三区免费视频网站| 窝窝影院91人妻| www.熟女人妻精品国产| 亚洲欧美日韩另类电影网站| 国产黄a三级三级三级人| svipshipincom国产片| 777久久人妻少妇嫩草av网站| 国产成年人精品一区二区 | 女性被躁到高潮视频| 免费在线观看完整版高清| e午夜精品久久久久久久| 美女国产高潮福利片在线看| 国产成年人精品一区二区 | 国产男靠女视频免费网站| 亚洲av熟女| 欧美久久黑人一区二区| 神马国产精品三级电影在线观看 | 人妻久久中文字幕网| 久久婷婷成人综合色麻豆| 免费高清在线观看日韩| 伦理电影免费视频| 97碰自拍视频| 久久精品91无色码中文字幕| 看免费av毛片| 一a级毛片在线观看| 男人操女人黄网站| 人妻丰满熟妇av一区二区三区| 亚洲伊人色综图| 国产成人精品久久二区二区91| 免费av中文字幕在线| 亚洲成人久久性| 日本免费a在线| 99久久综合精品五月天人人| 亚洲av五月六月丁香网| 亚洲精品一二三| 国产精品成人在线| 9191精品国产免费久久| 亚洲精品在线美女| 长腿黑丝高跟| 亚洲一卡2卡3卡4卡5卡精品中文| 精品高清国产在线一区| 亚洲自拍偷在线| 黑丝袜美女国产一区| 久久精品国产99精品国产亚洲性色 | 男女做爰动态图高潮gif福利片 | 最新在线观看一区二区三区| 波多野结衣av一区二区av| 日本撒尿小便嘘嘘汇集6| 久久久久亚洲av毛片大全| 国产精品成人在线| 热99国产精品久久久久久7| 黄片小视频在线播放| 少妇的丰满在线观看| 亚洲成人精品中文字幕电影 | 久久婷婷成人综合色麻豆| 国产色视频综合| 在线观看免费视频日本深夜| 超色免费av| 夜夜看夜夜爽夜夜摸 | 国产精品国产高清国产av| 日本a在线网址| 日韩国内少妇激情av| 亚洲aⅴ乱码一区二区在线播放 | 亚洲色图 男人天堂 中文字幕| 亚洲中文字幕日韩| 欧美成狂野欧美在线观看| 999久久久精品免费观看国产| 欧美精品一区二区免费开放| 黑人巨大精品欧美一区二区蜜桃| 操美女的视频在线观看| 亚洲国产精品sss在线观看 | 美女大奶头视频| www国产在线视频色| 一级a爱片免费观看的视频| 老司机福利观看| av在线播放免费不卡| 美女扒开内裤让男人捅视频| 国产蜜桃级精品一区二区三区| 黑人欧美特级aaaaaa片| av有码第一页| 91在线观看av| 亚洲一码二码三码区别大吗| 一级a爱视频在线免费观看| 18美女黄网站色大片免费观看| 麻豆成人av在线观看| 最近最新免费中文字幕在线| 欧美性长视频在线观看| 欧美日韩亚洲综合一区二区三区_| 午夜福利影视在线免费观看| 91av网站免费观看| 欧美不卡视频在线免费观看 | 日日摸夜夜添夜夜添小说| 亚洲精品粉嫩美女一区| 免费观看精品视频网站| 亚洲激情在线av| 最新美女视频免费是黄的| 一级,二级,三级黄色视频| 国产片内射在线| 久久久国产欧美日韩av| 最近最新免费中文字幕在线| 亚洲午夜理论影院| 99国产精品免费福利视频| 国产精品永久免费网站| 精品一区二区三卡| 女警被强在线播放| 激情在线观看视频在线高清| 啦啦啦 在线观看视频| 亚洲avbb在线观看| 男人舔女人的私密视频| 欧美色视频一区免费| 777久久人妻少妇嫩草av网站| xxxhd国产人妻xxx| 国产不卡一卡二| x7x7x7水蜜桃| 国产一区二区三区视频了| 亚洲aⅴ乱码一区二区在线播放 | 天堂俺去俺来也www色官网| 叶爱在线成人免费视频播放| 人妻丰满熟妇av一区二区三区| 18禁国产床啪视频网站| 精品国产乱子伦一区二区三区| 久久久久久久午夜电影 | 18美女黄网站色大片免费观看| 久久久久久免费高清国产稀缺| 18禁黄网站禁片午夜丰满| 久久国产乱子伦精品免费另类| 久久伊人香网站| 亚洲激情在线av| 精品国产乱子伦一区二区三区| 大陆偷拍与自拍| 久久精品亚洲熟妇少妇任你| 1024视频免费在线观看| 久久久久久人人人人人| 亚洲成国产人片在线观看| 国产不卡一卡二| 国产精品永久免费网站| 国产亚洲欧美在线一区二区| 热99国产精品久久久久久7| 亚洲五月天丁香| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 日韩 欧美 亚洲 中文字幕| www.熟女人妻精品国产| videosex国产| 日韩精品中文字幕看吧| 久久性视频一级片| 婷婷精品国产亚洲av在线| 校园春色视频在线观看| 国产一卡二卡三卡精品| 色播在线永久视频| 电影成人av| 曰老女人黄片| 国产亚洲精品久久久久5区| 国产高清国产精品国产三级| 黄色毛片三级朝国网站| 精品国产美女av久久久久小说| 成年女人毛片免费观看观看9| 亚洲一区二区三区色噜噜 | 天堂√8在线中文| 久久影院123| 国产真人三级小视频在线观看| 女警被强在线播放| 精品国产超薄肉色丝袜足j| 交换朋友夫妻互换小说| 成人永久免费在线观看视频| 美女高潮到喷水免费观看| 人人妻人人添人人爽欧美一区卜| 国产主播在线观看一区二区| 亚洲免费av在线视频| ponron亚洲| 香蕉国产在线看| 动漫黄色视频在线观看| 99在线视频只有这里精品首页| 夜夜看夜夜爽夜夜摸 | 夫妻午夜视频| 国产又色又爽无遮挡免费看| 成人国语在线视频| 精品久久久久久,| 性欧美人与动物交配| 最好的美女福利视频网| 热99re8久久精品国产| 黄色 视频免费看| 后天国语完整版免费观看| 又紧又爽又黄一区二区| 国产成人精品无人区| 免费人成视频x8x8入口观看| 91老司机精品| 国产免费av片在线观看野外av| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品第一综合不卡| 日本一区二区免费在线视频| 欧美日韩黄片免| 1024视频免费在线观看| 丝袜美足系列| 国产精品九九99| 中文字幕av电影在线播放| 成年人免费黄色播放视频| 人妻丰满熟妇av一区二区三区| 亚洲久久久国产精品| 9色porny在线观看| 久久中文看片网| 国产av一区在线观看免费| av在线播放免费不卡| 色综合欧美亚洲国产小说| 国产精品久久久av美女十八| 老司机午夜十八禁免费视频| 大型黄色视频在线免费观看| 亚洲色图 男人天堂 中文字幕| 亚洲成人久久性| 免费在线观看亚洲国产| 人人妻,人人澡人人爽秒播| x7x7x7水蜜桃| 亚洲少妇的诱惑av| 日韩一卡2卡3卡4卡2021年| 日本 av在线| 日日爽夜夜爽网站| 母亲3免费完整高清在线观看| 日韩免费av在线播放| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区久久 | 最近最新免费中文字幕在线| 国产精品综合久久久久久久免费 | 国产99久久九九免费精品| 亚洲一码二码三码区别大吗| 欧美亚洲日本最大视频资源| 欧美成人免费av一区二区三区| 国产一区二区激情短视频| 中文字幕精品免费在线观看视频| 香蕉国产在线看| 亚洲第一av免费看| 久久影院123| 免费女性裸体啪啪无遮挡网站| 女人被躁到高潮嗷嗷叫费观| 一级片'在线观看视频| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 97超级碰碰碰精品色视频在线观看| 99精品欧美一区二区三区四区| 在线十欧美十亚洲十日本专区| 波多野结衣一区麻豆| 男女下面进入的视频免费午夜 | 69精品国产乱码久久久| 在线观看免费视频网站a站| 国产精品电影一区二区三区| 18美女黄网站色大片免费观看| 婷婷精品国产亚洲av在线| 两个人看的免费小视频| 欧美在线一区亚洲| 天天躁夜夜躁狠狠躁躁| 国产高清videossex| 一夜夜www| 女同久久另类99精品国产91| 久久精品国产99精品国产亚洲性色 | 午夜91福利影院| 国产免费现黄频在线看| 老司机午夜十八禁免费视频| 水蜜桃什么品种好| 久久性视频一级片| 亚洲成人免费av在线播放| 伦理电影免费视频|