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

    張廣才嶺北部三大硬闊樹木生長-氣候關(guān)系的時空變異

    2017-05-25 00:37:30蘇金娟王曉春
    生態(tài)學(xué)報 2017年5期
    關(guān)鍵詞:水曲柳最低溫度年表

    蘇金娟, 王曉春

    東北林業(yè)大學(xué)生態(tài)研究中心, 哈爾濱 150040

    張廣才嶺北部三大硬闊樹木生長-氣候關(guān)系的時空變異

    蘇金娟, 王曉春*

    東北林業(yè)大學(xué)生態(tài)研究中心, 哈爾濱 150040

    樹木年輪在時空尺度上的比較可以更好地反映環(huán)境變化對樹木生長的影響,在認(rèn)識氣候變化對森林生態(tài)系統(tǒng)的影響上具有重要意義。采用樹木年代學(xué)方法分析了張廣才嶺北部地區(qū)闊葉紅松林中主要闊葉樹種-水曲柳(Fraxinusmandshurica)、黃菠蘿(Phellodendronamurense)和胡桃楸(Juglansmandshurica)徑向生長與氣候關(guān)系的時空變異。結(jié)果表明,在同一地點樹種間氣候響應(yīng)差異明顯,胡桃楸受降水和最低溫度共同作用,而黃菠蘿和水曲柳則主要受最低溫度限制,這表明樹木生長與氣候因子的關(guān)系具有一定的物種特異性。隨著溫度和降水格局的改變,三大硬闊年輪與氣候關(guān)系在空間水平上存在差異,方正和西大圈樣點的胡桃楸與6—8月最低溫度呈顯著正相關(guān)(P<0.05),而鳳凰山樣點與5、6月最低溫度和降水呈正相關(guān)(P<0.05);黃菠蘿和水曲柳徑向生長隨著降水空間格局的變化,其生長季末期相關(guān)性程度逐漸減弱。1980年后張廣才嶺北部出現(xiàn)氣溫顯著升高,在升溫前三大硬闊的生長趨勢相對一致,而在升溫后黃菠蘿和水曲柳樹輪寬度隨溫度升高呈上升趨勢,但胡桃楸卻出現(xiàn)隨溫度升高而生長下降的“分異現(xiàn)象”。如果未來增溫趨勢持續(xù)發(fā)生或者加重,可以推斷在張廣才嶺北部胡桃楸可能受干旱脅迫加劇,其可能出現(xiàn)生長衰退,但增溫可能更有利于黃菠蘿和水曲柳的生長。

    樹木年輪;闊葉紅松林;干旱脅迫;全球變暖;徑向生長

    近百年來,地球氣候正經(jīng)歷著一次以變暖為主要特征的顯著變化。全球變暖的影響在北半球中高緯度地區(qū)最為突出,尤其是在內(nèi)陸生態(tài)系統(tǒng)[1- 4]。越來越多的數(shù)據(jù)表明中國東北地區(qū)自1970年以來增溫幅度顯著提高,這種明顯的增溫導(dǎo)致從1980年開始進(jìn)入一種暖干化趨勢。這種變暖會在不同尺度上影響森林生態(tài)系統(tǒng)的生產(chǎn)力,已有研究表明氣候變暖后森林生產(chǎn)力下降[5- 7],也有研究顯示氣候變暖會改變森林樹木的生長趨勢[8]、物候和分布地區(qū)[9]。張廣才嶺北部的闊葉紅松林,也是氣候變化最為顯著的地區(qū)之一。三大硬闊-水曲柳、黃菠蘿和胡桃楸是闊葉紅松林中的重要樹種,也是紅松(Pinuskoraiensis)的主要伴生樹種,它們在森林環(huán)境形成、紅松更新及闊葉紅松林的結(jié)構(gòu)與功能動態(tài)中起著非常重要的作用。此外,由于樹木生理特性的限制,樹輪的清晰程度有所不同,闊葉樹種樹輪研究較針葉樹種更為困難,闊葉樹種樹輪研究相對較少,但是闊葉樹種年輪可能含有與針葉樹種不一樣的信息[10-11]。因此,分析三大硬闊徑向生長與氣候關(guān)系,是探索整個闊葉紅松林生態(tài)系統(tǒng)對氣候變化響應(yīng)的重要組成部分。

    樹木年代學(xué)在時間和空間尺度上的研究能夠為樹木徑向生長時空變化提供可靠的信息[12]。目前,國外對闊葉樹種樹輪-氣候響應(yīng)時空變化研究已有一些報道,例如,Huang等[13]發(fā)現(xiàn)加拿大美洲山楊(Populustremuloides)和北美白樺(Betulapapyrifera)在南部主要受夏季濕度影響,而在北部主要受溫度影響。Pierre等[9]發(fā)現(xiàn)隨著經(jīng)度變化,夏季降水減少,溫度增加,無?;?Quercuspetraea)對5、6月份溫度和降水的相關(guān)性也減弱。另外,在中歐和法國西部對櫟樹以及其它闊葉樹種的研究均表明海拔、經(jīng)度、緯度變化對樹木的生長有不同的影響[14- 16]。在國內(nèi)張寒松等[17]用水曲柳年表重建了長白山過去240年降水變化;而王蔚蔚等[18]發(fā)現(xiàn)長白山水曲柳對溫度信號更敏感,重建了長白山地區(qū)秋季溫度;李牧等[19]在敦化地區(qū)用紅松和三大硬闊年表重建了生長季最低溫度。這些研究都集中在空間分析和重建上,然而過去幾十年研究表明,在我國青藏高原和東北地區(qū)樹木生長對氣候響應(yīng)是不穩(wěn)定的,在氣候變暖后出現(xiàn)了明顯的轉(zhuǎn)型[20- 22],即樹木生長分異現(xiàn)象[23]。那么在全球變暖的背景下在張廣才嶺北部三大硬闊的樹木生長分異現(xiàn)象時空格局會發(fā)生變化呢?

    本文以張廣才嶺北部三大硬闊樹輪為研究對象,探討“三大硬闊樹輪與氣候關(guān)系在空間格局上的響應(yīng)變化以及氣候變暖后在時間上是否具有穩(wěn)定性”這一科學(xué)問題,研究結(jié)果旨在揭示三大硬闊樹輪生長對氣候變化的不同響應(yīng)特點,以期了解未來氣候變暖對于森林生態(tài)系統(tǒng)的影響。

    1 材料和方法

    1.1 采樣點概況

    在黑龍江省張廣才嶺北部地區(qū)沿經(jīng)度選取3個人為干擾較少的地區(qū)(圖1):方正縣雙子山國家森林公園、勃利縣西大圈森林公園和雞東縣鳳凰山國家級自然保護(hù)區(qū),作為本研究的采樣點(表1)。研究區(qū)域?qū)儆陂L白山系支脈張廣才嶺北部,該區(qū)域?qū)儆跍貛Т箨懶詺夂?四季較為分明,雨熱同期,年均溫2.8—5℃,無霜期100—160 d,年總降水量400—600 mm(圖2,表2),土壤多為暗棕色森林土。張廣才嶺地貌形態(tài)復(fù)雜,地形以山地為主,自然條件豐富,植物種類眾多,植被區(qū)系類型以溫帶分布為主。主要森林植被類型為針闊葉混交林,主要包括白樺(Betulaplatyphylla)林、落葉闊葉林、雜木林、闊葉紅松林以及蒙古櫟(Quercusmongolica)林等多種類型。喬木樹種以紅松、黃菠蘿、水曲柳、胡桃楸、白樺、紫椴(Tiliaamurensis)和蒙古櫟等為主;灌木植物以金銀忍冬(Loniceramaackii)、毛榛子(Corylusmandshurica)、胡枝子(Lespedezabicolor)等為優(yōu)勢種;草本植物尖齒蹄蓋蕨(Athyriumdentilobum)、小葉芹(Aegopodiumalpestre)、木賊(Equisetumhyemale)等為主。張廣才嶺北部是東亞針闊葉混交林分布中心,三大硬闊——胡桃楸、黃菠蘿、水曲柳廣泛分布于張廣才嶺,是張廣才嶺北部山區(qū)頂極群落(闊葉紅松林)的主要闊葉樹種,并且采樣地區(qū)受人為干擾較少,是研究樹輪-氣候關(guān)系的理想?yún)^(qū)域。

    圖1 三大硬闊樹木年輪采樣點和氣象站分布Fig.1 Map of sampled sites for three hardwood tree rings and weather station near its sampled site

    1.2 樣本采集及年表建立

    圖2 研究地區(qū)氣象站月平均氣溫和月總降水量變化 Fig.2 Monthly mean temperature and total precipitation at each study area

    2014年7月在方正縣雙子山國家森林公園(FZ),勃利縣西大圈森林公園(XDQ)和雞東縣鳳凰山國家級自然保護(hù)區(qū)(FHS)分別在北坡下坡位進(jìn)行樹輪采樣(表1,圖1)。選擇生長良好的樹木,在胸高處(1.3 m)用生長錐鉆取樹輪樣本。將取到的樣本放入塑料管內(nèi),進(jìn)行編號,并記錄每棵樹的生境。每個樹種采集15—20棵樹,每棵樹采集2個樹輪樣芯。

    將樣品帶回實驗室進(jìn)行固定、干燥、打磨后,用骨架示意圖方法進(jìn)行交叉定年[24],之后用Velmex年輪分析儀進(jìn)行測量,精確到0.001 mm。測量后用COFECHA程序[25]對交叉定年和測量結(jié)果進(jìn)行檢驗,以確保其結(jié)果準(zhǔn)確性。經(jīng)交叉定年后的樹輪序列,利用ARSTAN程序進(jìn)行去趨勢和標(biāo)準(zhǔn)化[26],去趨勢方法采用負(fù)指數(shù)或線性函數(shù)進(jìn)行,降低因年齡生長速率不同而帶來的影響。用雙重平均法進(jìn)行年輪曲線標(biāo)準(zhǔn)化,最終建立3個采樣點黃菠蘿、胡桃楸和水曲柳9個標(biāo)準(zhǔn)年表(圖3)用以進(jìn)行年輪-氣候響應(yīng)分析。

    1.3 統(tǒng)計分析

    選擇離采樣較近的通河和雞西氣象站分別對應(yīng)方正(FZ)和鳳凰山(FHS)采樣點(表1),由于西大圈與依蘭和雞西氣象站距離相等,本研究將依蘭站和雞西站氣象站平均作為西大圈氣象數(shù)據(jù)。氣象數(shù)據(jù)由中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc.cma.gov.cn/)提供(表2)。氣候因子主要包括:月均最低溫(Tmin)、月均最高溫(Tmax)、月均溫(T)、月總降水量(P)、相對濕度(RH),其中相對濕度是由空氣中的絕對濕度與同溫度下的飽和絕對濕度的比值,得出一個百分比。由于上一年氣候因子可能影響當(dāng)年樹木徑向生長[24],因此分析了上一年10月到當(dāng)年9月共12個月的氣候因子與黃菠蘿、胡桃楸和水曲柳標(biāo)準(zhǔn)年表的相關(guān)性。由于氣候變化對樹木徑向生長有累積和長期的影響[22],本研究分析了年際、季節(jié)與樹木徑向生長的關(guān)系。結(jié)合生理生態(tài)意義將季節(jié)劃分為:冬季(前一年12月至當(dāng)年2月)、春季(3—5月)、夏季(6—8月)、秋季(9—11月)。

    表2 不同采樣點對應(yīng)氣象站點信息

    根據(jù)相關(guān)分析結(jié)果分別將同一地點水曲柳、胡桃楸和黃菠蘿提取同一地點主成分年表和不同地點同一樹種的主成分年表用于分析最低溫度對不同地點間樹種的影響。采用SPSS19.0完成主成分提取及相關(guān)統(tǒng)計分析;采用Sigmapolt 12.5等軟件完成圖表制作。

    圖3 不同地點三大硬闊標(biāo)準(zhǔn)年表年輪指數(shù)變化Fig.3 Variations of the standard chronologies for three hardwoods tree species at different sitesFZ、XDQ、FHS分別代表方正、西大圈、鳳凰山

    2 結(jié)果與分析

    2.1 年表特征

    在9個三大硬闊標(biāo)準(zhǔn)年表中,不同樣點間三大硬闊年齡差異較大,其中方正水曲柳的樹齡最大(208a),鳳凰山胡桃楸的樹齡最小(64a)。平均敏感度反映了年輪寬度對氣候變化的敏感程度,胡桃楸的平均敏感度普遍高于黃菠蘿和水曲柳,這說明了胡桃楸對氣候信號更敏感。一階自相關(guān)反映了上一年氣候狀況對當(dāng)年年輪寬度生長的影響,一階自相關(guān)都較高,說明了氣候變化對三大硬闊生長影響的“滯后效應(yīng)”較為明顯。9個標(biāo)準(zhǔn)年表的樣本總體代表性都高于0.85(表3),表明這些樣本可以代表種群的總體水平??傮w來說,9個標(biāo)準(zhǔn)年表的平均敏感度、信噪比和樣本代表性都很高,年表質(zhì)量較好,這說明這些年表適合用于樹木年輪學(xué)氣候研究。

    表3 三大硬闊標(biāo)準(zhǔn)年表特征的主要統(tǒng)計特征

    9個年表間在整體上具有較高的一致性,但局部也存在差異(表4)。方正水曲柳和西大圈黃菠蘿年表間相關(guān)性最高(r=0.64),鳳凰山水曲柳與方正胡桃楸和方正水曲柳年表間相關(guān)系數(shù)均最低(r=-0.03, 0.03),這可能與采樣點空間距離、樹種有關(guān)。另外,同一地點不同樹種間的相關(guān)系數(shù)也存在差異,鳳凰山胡桃楸與黃菠蘿年表間的相關(guān)系數(shù)是0.39,而與水曲柳年表的相關(guān)為0.15,這表明樹木的生理特征差異也會在一定程度影響年表對氣候的響應(yīng)。

    表4 不同地點、樹種標(biāo)準(zhǔn)年表間相關(guān)系數(shù)(公共區(qū)間1949—2013年)

    2.2 三大硬闊樹輪與氣候關(guān)系的空間變異

    從9個三大硬闊標(biāo)準(zhǔn)年表與氣候因子相關(guān)分析結(jié)果來看(圖4),不同樣點三大硬闊與氣候因子響應(yīng)關(guān)系存在差異,方正胡桃楸徑向生長與生長季(6—9月)最低溫度和降水均呈明顯正相關(guān),其中年表序列與6、7月最低溫度及6、7、9月降水均呈顯著正相關(guān);西大圈胡桃楸與生長季最低溫度和降水相關(guān)性較弱,僅與9月份降水顯著正相關(guān);鳳凰山胡桃楸在生長季前(3—5月)及生長季與最低溫度均呈正相關(guān),其中5、6月最低溫度呈顯著正相關(guān),而8、9月降水則呈負(fù)相關(guān)。黃菠蘿和水曲柳對降水相關(guān)性不明顯,在各個樣點均表現(xiàn)較弱程度的正相關(guān),各樣點黃菠蘿和水曲柳與最低溫度呈明顯的正相關(guān),且鳳凰山樣點年表序列與最低溫度的相關(guān)程度明顯高于西大圈和方正樣點。

    圖4 三大硬闊標(biāo)準(zhǔn)年表與月最低溫度、降水的相關(guān)系數(shù)Fig.4 Correlation coefficients between standard tree-ring chronologies and monthly climate variables (monthly mean minimum temperatures and monthly total precipitation) for different sites and tree species-表示前一年

    圖5 三大硬闊標(biāo)準(zhǔn)年表與高溫、均溫、相對濕度的相關(guān)性Fig.5 Correlation coefficients between the standard tree-ring chronologies and monthly climate variables (monthly mean maximum temperature, mean temperature, and monthly mean relative humidity) for three hardwoods at different sites

    從圖5可以看出,方正胡桃楸徑向生長與4、7月最高溫負(fù)相關(guān),與生長季相對濕度呈正相關(guān);西大圈胡桃楸徑向生長與溫度、相對濕度相關(guān)性明顯相反,與溫度呈現(xiàn)一致的負(fù)相關(guān),而相對濕度則呈現(xiàn)正相關(guān),其中7—9月相對濕度呈顯著正相關(guān);鳳凰山胡桃楸與最高溫、均溫以及相對濕度相關(guān)性較弱,僅與7月最高溫、相對濕度達(dá)到顯著水平。不同樣點黃菠蘿徑向生長與最高溫、均溫、相對濕度差異明顯,方正、西大圈樣點相對濕度均呈正相關(guān),而鳳凰山樣點則表現(xiàn)明顯的負(fù)相關(guān),另外鳳凰山樣點與均溫相關(guān)程度明顯高于方正、西大圈。3個樣點水曲柳徑向生長與生長季相對濕度均呈正相關(guān),最高溫表現(xiàn)明顯不同,方正樣點在生長季末期最高溫明顯正相關(guān),而鳳凰山和西大圈則表現(xiàn)較弱的負(fù)相關(guān)。

    表5 三大硬闊標(biāo)準(zhǔn)年表與季節(jié)氣候因子相關(guān)性

    *P<0.05,**P<0.01

    從表6中可以看出,方正年際、季節(jié)溫度和降水對胡桃楸徑向生長均呈正相關(guān),其中夏季的最低溫度和降水對胡桃楸的徑向生長具有極顯著的促進(jìn)作用;西大圈胡桃楸對最低溫度和降水則表現(xiàn)較弱程度的正相關(guān),僅秋季的降水達(dá)到了顯著水平(P<0.05);鳳凰山胡桃楸較方正和西大圈對最低溫度、降水的相關(guān)性程度有所增強(qiáng),與春季和秋季的最低溫度都達(dá)到了顯著正相關(guān),與春季降水達(dá)到極顯著水平。3個樣點的黃菠蘿和水曲柳年際、季節(jié)相關(guān)結(jié)果表明(表6),年最低溫度與黃菠蘿和水曲柳都達(dá)到了極顯著正相關(guān)(P<0.01),而降水相關(guān)性均呈較弱的正相關(guān)關(guān)系。這說明了最低溫度是限制黃菠蘿和水曲柳徑向生長的主要限制因子。此外3個采樣點從東到西降水在減少,其黃菠蘿和水曲柳與降水相關(guān)性程度隨樣點變化逐漸減弱。

    2.3 三大硬闊樹輪-氣候關(guān)系的時間變異

    3個采樣點溫度和降水變化表明(圖6),3個樣點在1980年左右發(fā)現(xiàn)溫度顯著升高,1980年前平均最低溫度緩慢上升甚至下降,但是1980年后平均最低溫度均以0.03℃/a速率持續(xù)上升。降水變化則不明顯,3個樣點1980年前后一致的表現(xiàn)為下降趨勢。

    圖6 1950—2013年研究區(qū)域年最低溫度和總降水量變化Fig.6 Variation of annual mean minimum temperature and total precipitation from 1950 to 2013

    三大硬闊不同樹種升溫前后年輪指數(shù)變化較大(圖7),升溫后3個樣點的胡桃楸生長速率較升溫前均呈下降趨勢;黃菠蘿和水曲柳年輪指數(shù)升溫前后變化波動較小,其變化規(guī)律基本一致,生長速率持續(xù)上升。提取同一地點不同樹種的第一主成分年表以及不同地點同一樹種第一主成分年表與最低溫度關(guān)系比較(圖8)結(jié)果顯示,3個樣點間年輪指數(shù)變化差異不明顯,1980年氣候變暖后生長速率均是上升趨勢,其上升幅度較升溫前均有所下降;但是不同樹種間差異較大,胡桃楸隨溫度升高年輪指數(shù)下降,而黃菠蘿和水曲柳具有與溫度升高相一致的趨勢。

    圖7 升溫(1980年)前后不同地點三大硬闊年輪指數(shù)隨最低溫度變化趨勢Fig.7 Changes of tree-ring index for three hardwoods with annual minimum temperature at different sites before and after rapid warming (1980)

    圖8 升溫(1980年)前后同一地點不同樹種以及不同地點相同樹種第一主成分隨最低溫度變化趨勢Fig.8 Changes of the first components of chronologies for the same tree species and sites with annual mean minimum temperature before and after rapid warming (1980)

    3 討論

    3.1 氣候因子和空間環(huán)境變化對三大硬闊徑向生長的影響

    三大硬闊徑向生長主要受溫度和降水限制,最低溫度是水曲柳和黃菠蘿徑向生長的主要限制因子,胡桃楸徑向生長受最低溫度和降水共同限制,這與敦化地區(qū)三大硬闊研究結(jié)果一致[19]。不同樹種徑向生長響應(yīng)差異可能與其遺傳特性有關(guān),水曲柳具有一定的耐陰性,喜濕潤環(huán)境并且光合速率最強(qiáng)[27-28],而溫度對樹木生理生態(tài)的最直接影響是改變光合速率,故水曲柳受最低溫度的限制作用最強(qiáng)。從對水分需求的生理特性來看,胡桃楸耗水量和萎蔫系數(shù)均大于黃菠蘿和水曲柳[28-29],并且水曲柳和黃菠蘿能夠適應(yīng)一定的干旱脅迫,而胡桃楸在干旱脅迫中不能生長[30-31]。因此,胡桃楸對于水分的限制作用大于黃菠蘿和水曲柳。彭劍鋒等[32]在天山地區(qū)對天山云杉(Piceaschrenkiana)和西伯利亞落葉松(Larixsibirica)的研究發(fā)現(xiàn)遺傳特性是造成樹木對同一生境響應(yīng)不同的原因。

    張廣才嶺北部屬于濕潤、半濕潤地區(qū),但由于地理空間位置的差異,3個樣點的環(huán)境條件并不完全一致,其不同樣點溫度、降水與樹木徑向生長響應(yīng)有所差異。方正樣點胡桃楸與生長季6、7、9月最低溫度、降水顯著正相關(guān)。方正樣點位于最西端,氣候較濕潤,溫度和降水共同調(diào)控胡桃楸徑向生長,6—9月是樹木生長最旺盛時期,良好的水熱條件能夠促進(jìn)胡桃楸快速生長,體現(xiàn)了胡桃楸對水熱條件的同步正響應(yīng)。西大圈樣點胡桃楸與生長季最低溫度、降水相關(guān)程度減弱,僅與9月最低溫度、降水呈顯著正相關(guān),主要是由于西大圈位于中間,其氣候條件較溫暖,加劇了水分的蒸發(fā),此外9月份是生長季末期,過多的降水不利于營養(yǎng)物質(zhì)的積累以及木質(zhì)化的進(jìn)程,導(dǎo)致水分對胡桃楸的促進(jìn)作用下降,兩者對胡桃楸徑向生長的正響應(yīng)均降低。鳳凰山胡桃楸與最低溫度和降水顯著相關(guān)性則出現(xiàn)在5、6月,這可能與鳳凰山地理位置更靠東部有關(guān),氣候條件為暖干,5月處于生長期早期,樹木對水分的需求量大,較高的溫度加速了冰雪融化和土壤溫度升溫,促進(jìn)了根系的活動和地上部分的萌發(fā),從而有利于生長季早期樹木生長,故鳳凰山胡桃楸提前表現(xiàn)顯著的正相關(guān)。黃菠蘿和水曲柳與降水相關(guān)性程度在生長季隨樣點變化逐漸減弱,這是因為在半濕潤地區(qū), 高溫可以加快植物蒸騰和土壤蒸發(fā), 加劇水分脅迫作用, 從而很大程度上抑制樹體的生長,當(dāng)降水量低于某一值時, 溫度與樹木生長呈現(xiàn)明顯負(fù)相關(guān), 降水量越少, 這種關(guān)系越明顯與桑衛(wèi)國[33]等對天山云杉研究結(jié)果一致。這似乎可以說明三大硬闊徑向生長隨著樣點(經(jīng)度)的變化,不同的水熱條件對其徑向生長影響有一定變化,但具體機(jī)制有待于進(jìn)一步研究。

    3.2 氣候變暖對三大硬闊年輪生長的影響

    20世紀(jì)80年代東北地區(qū)溫度顯著升高[34- 36],其中最低氣溫增高更明顯[37],最低溫升溫對植物的生長影響是復(fù)雜的,有研究證明最低溫度升高對植物生長有促進(jìn)作用[38- 40],也有人認(rèn)為最低溫度升高對植物生長有抑作用[7]。最低溫度升高對三大硬闊徑向生長在不同樹種間的差異明顯。1980年前三大硬闊生長趨勢基本一致(圖7),1980年后黃菠蘿和水曲柳徑向生長隨著溫度升高而上升,但胡桃楸在升溫后并未隨溫度升高而呈現(xiàn)上升趨勢,反而出現(xiàn)下降趨勢,胡桃楸的主成分年表結(jié)果也呈現(xiàn)這一趨勢,胡桃楸徑向生長與最低溫度升高出現(xiàn)了“分異現(xiàn)象”。響應(yīng)分異現(xiàn)象廣泛存在于北半球中高緯度地區(qū)[20-21,41]。胡桃楸徑向生長受降水和最低溫度共同作用,最低溫度升高加快植物蒸騰作用和土壤蒸發(fā),從而導(dǎo)致干旱脅迫抑制胡桃楸徑向生長。從胡桃楸年輪與氣候相關(guān)分析可以推測,升溫造成的干旱脅迫可能是胡桃楸出現(xiàn)分異現(xiàn)象主要原因,而且隨著自西向東降水增多,其響應(yīng)分異程度也隨著經(jīng)度梯度逐漸減弱,這似乎證明了這一結(jié)論。這與長白山魚鱗云杉(Piceajezoensisvar.komarovii)上限樹輪與溫度出現(xiàn)的分異現(xiàn)象一致[39]。李宗善等[42]對川西地區(qū)岷江冷杉(Abiesfaxoniana)研究發(fā)現(xiàn)日照時數(shù)的變化是引起分異現(xiàn)象的原因,川西地區(qū)屬于青藏高原氣候,海拔高、空氣稀薄、日照時數(shù)長對植物的影響較大。黃菠蘿和水曲柳主要受最低溫度限制,降水對其影響較小,升溫導(dǎo)致的干旱加劇可能沒有成為主要的限制因素,即使夏季最低溫度有一定程度上升,可能目前并沒有超過這兩種樹種的溫度界限,因此并沒有發(fā)生顯著的變化,尚華明等[43]也證明了這一結(jié)論。另外,最低溫度升高沒有影響到黃菠蘿和水曲柳的同時,也可能降低了其對環(huán)境變化的敏感性,即增強(qiáng)了其對高溫脅迫的抵抗力,使其在顯著升溫后沒有出現(xiàn)響應(yīng)分異現(xiàn)象。雖然1980年后的顯著升溫,黃菠蘿和水曲柳徑向生長沒有顯著變化,但是如果持續(xù)變暖,尤其是夏季最低溫度升高,超出其耐受溫度或者引起干旱脅迫,也會對其生長趨勢有一定影響。

    4 結(jié)論

    張廣才嶺北部闊葉紅松林是氣候變化最為顯著的地區(qū)之一,本研究通過對闊葉紅松林中水曲柳、黃菠蘿和胡桃楸徑向生長與氣候關(guān)系分析表明:在空間尺度上,方正生長季同步的水熱條件促進(jìn)了胡桃楸徑向生長,西大圈胡桃楸在生長季與最低溫度、降水則表現(xiàn)較弱的正相關(guān),鳳凰山胡桃楸與5、6月最低溫度、降水呈現(xiàn)正相關(guān)。黃菠蘿和水曲柳徑向生長與經(jīng)度變化吻合,隨著從西向東降水增多,其年表序列與降水的正相關(guān)作用逐漸減弱。在時間尺度上,20世紀(jì)80年代溫度顯著升高后,三大硬闊的徑向生長在不同樹種間出現(xiàn)明顯差異,黃菠蘿和水曲柳在顯著升溫后生長趨勢沒有明顯變化,但胡桃楸生長趨勢出現(xiàn)與升溫相反的“分異現(xiàn)象”。未來在持續(xù)變暖的環(huán)境下胡桃楸可能出現(xiàn)衰退,而升溫對黃菠蘿和水曲柳的影響不大。三大硬闊與氣候關(guān)系的時空變異機(jī)理尚不完全清楚,將來可通過三大硬闊的樹輪解剖學(xué)或生理生態(tài)響應(yīng)實驗進(jìn)一步分析其時空變異機(jī)制。

    [1] Serreze M C, Walsh J E, Chapin III F S, Osterkamp T, Dyurgerov M, Romanovsky V, Oechel W C, Morison J, Zhang T, Barry R G. Observational evidence of recent change in the northern high-latitude environment. Climatic Change, 2000, 46(1/2): 159- 207.

    [2] Andreu L, Gutiérrez E, Macias M, Ribas M, Bosch O, Camarero J J. Climate increases regional tree-growth variability in Iberian pine forests. Global Change Biology, 2007, 13(4): 804- 815.

    [3] Diego Galván J, Büntgen U, Ginzler C, Grudd H, Gutiérrez E, Labuhn I, Camarero J J. Drought-induced weakening of growth-temperature associations in high-elevation Iberian pines. Global and Planetary Change, 2015, 124: 95- 106.

    [4] IPCC. Summary for policymakers of climate change 2013: the physical science basis//Stocker T F, Qin D, Plattner G K, Tignor M, Allen S K, Boschung J, Nauels A, Xia Y, Bex V, Midgley P M. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press, 2014.

    [5] 朱建華, 侯振宏, 張治軍, 羅云建, 張小全. 氣候變化與森林生態(tài)系統(tǒng): 影響、脆弱性與適應(yīng)性. 林業(yè)科學(xué), 2007, 43(11): 138- 145.

    [6] 劉國華, 傅伯杰. 全球氣候變化對森林生態(tài)系統(tǒng)的影響. 自然資源學(xué)報, 2001, 16(1): 71- 78.

    [7] Clark D A, Piper S C, Keeling C D, Clark D B. Tropical rain forest tree growth and atmospheric carbon dynamics linked to interannual temperature variation during 1984- 2000. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(10): 5852- 5857.

    [8] D′Arrigo R D, Kaufmann R K, Davi N, Jacoby G C, Laskowski C, Myneni R B, Cherubini P. Thresholds for warming-induced growth decline at elevational tree line in the Yukon Territory, Canada. Global Biogeochemical Cycles, 2004, 18(3): GB3021, doi: 10.1029/2004GB002249.

    [9] Mérian P, Bontemps J D, Bergès L, Lebourgeois F. Spatial variation and temporal instability in climate-growth relationships of sessile oak (Quercuspetraea[Matt.] Liebl.) under temperate conditions. Plant Ecology, 2011, 212(11): 1855- 1871.

    [10] 高露雙, 王曉明, 趙秀海. 長白山闊葉紅松林共存樹種徑向生長對氣候變化的響應(yīng). 北京林業(yè)大學(xué)學(xué)報, 2013, 35(3): 24- 31.

    [11] Takahashi K, Okuhara I. Comparison of climatic effects on radial growth of evergreen broad-leaved trees at their northern distribution limit and co-dominating deciduous broad-leaved trees and evergreen conifers. Ecological Research, 2012, 27(1): 125- 132.

    [12] Tardif J, Camarero J J, Ribas M, Gutiérrez E. Spatiotemporal variability in tree growth in the Central Pyrenees: climatic and site influences[J]. Ecological Monographs, 2003, 73(2): 241- 257.

    [13] Huang J G, Tardif J C, Bergeron Y, Denneler B, Berninger F, Girardin M P. Radial growth response of four dominant boreal tree species to climate along a latitudinal gradient in the eastern Canadian boreal forest. Global Change Biology, 2010, 16(2): 711- 731.

    [14] Lara A, Villalba R, Wolodarsky-Franke A, Aravena J C, Luckman B H, Cuq E. Spatial and temporal variation inNothofaguspumiliogrowth at tree line along its latitudinal range (35°40′- 55°S) in the Chilean Andes. Journal of Biogeography, 2005, 32(5): 879- 893.

    [15] Friedrichs D A, Trouet V, Büntgen U, Frank D C, Esper J, Neuwirth B, L?ffler J. Species-specific climate sensitivity of tree growth in Central-West Germany. Trees, 2009, 23(4): 729- 739.

    [17] 張寒松, 韓士杰, 李玉文, 張軍輝. 利用樹木年輪寬度資料重建長白山地區(qū)240年來降水量的變化. 生態(tài)學(xué)雜志, 2007, 26(12): 1924- 1929.

    [18] 王蔚蔚, 張軍輝, 戴冠華, 王秀秀, 韓士杰, 張寒松, 王云. 利用樹木年輪寬度資料重建長白山地區(qū)過去240年秋季氣溫的變化. 生態(tài)學(xué)雜志, 2012, 31(4): 787- 793.

    [19] 李牧, 王曉春. 敦化三大硬闊、紅松年輪-氣候關(guān)系及生長季低溫重建. 南京林業(yè)大學(xué)學(xué)報:自然科學(xué)版, 2013, 37(3): 29- 34.

    [20] Zhang Y X, Wilmking M. Divergent growth responses and increasing temperature limitation of Qinghai spruce growth along an elevation gradient at the northeast Tibet Plateau. Forest Ecology and Management, 2010, 260(6): 1076- 1082.

    [21] Zhang Y X, Wilmking M, Gou X H. Changing relationships between tree growth and climate in Northwest China. Plant Ecology, 2009, 201(1): 39- 50.

    [22] Yu D P, Liu J Q, Benard J L, Zhou L, Zhou W M, Fang X M, Wei Y W, Jiang S W, Dai L M. Spatial variation and temporal instability in the climate-growth relationship of Korean pine in the Changbai Mountain region of Northeast China. Forest Ecology and Management, 2013, 300: 96- 105.

    [23] Briffa K R, Schweingruber F H, Jones P D, Osborn T J, Shiyatov, S G, Vaganov E A. Reduced sensitivity of recent tree-growth to temperature at high northern latitudes. Nature, 1998, 391(6668): 678- 682.

    [24] Fritts H C. Relationships of ring widths in arid-site conifers to variations in monthly temperature and precipitation. Ecological Monographs, 1974, 44(4): 411- 440.

    [25] Holmes R L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin, 1983, 43: 69- 78.

    [26] Cook E, Holmes R L. User Manual for Program ARSTAN, Laboratory of Tree-Ring Research. Tucson: University of Arizona, 1986.

    [27] 馬建路. 水曲柳、胡桃楸、黃波羅和紫椴適生立地條件的對比研究. 林業(yè)資源管理, 1996, (4): 52- 56.

    [28] 王晶英, 王展, 趙雨森. 水曲柳等樹種光合、耗水及萎蔫系數(shù)比較研究. 林業(yè)科技, 2010, 35(3): 10- 13.

    [29] 孫慧珍, 孫龍, 王傳寬, 周曉峰. 東北東部山區(qū)主要樹種樹干液流研究. 林業(yè)科學(xué), 2005, 41(3): 36- 42.

    [30] 梁淑娟, 潘攀, 孫志虎, 王慶成. 坡位對水曲柳及胡桃楸生長的影響. 東北林業(yè)大學(xué)學(xué)報, 2005, 33(3): 18- 19.

    [31] 肖冬梅, 王淼, 姬蘭柱. 水分脅迫對長白山闊葉紅松林主要樹種生長及生物量分配的影響. 生態(tài)學(xué)雜志, 2004, 23(5): 93- 97.

    [32] 彭劍峰, 勾曉華, 陳發(fā)虎, 李金豹. 天山云杉和西伯利亞落葉松的樹輪氣候記錄. 生態(tài)環(huán)境, 2005, 14(4): 460- 465.

    [33] 桑衛(wèi)國, 王云霞, 蘇宏新, 陸兆華. 天山云杉樹輪寬度對梯度水分因子的響應(yīng). 科學(xué)通報, 2007, 52(19): 2292- 2298.

    [34] 王紹武, 趙宗慈, 唐國利. 中國的氣候變暖. 國際政治研究, 2009, (4): 1- 11, 194- 194.

    [35] 周秀杰, 王鳳玲, 吳玉影, 那濟(jì)海, 潘華盛, 王勇. 近60年來黑龍江省與東北及全國氣溫變化特點分析. 自然災(zāi)害學(xué)報, 2013, 22(2): 124- 129.

    [36] 李騰, 何興元, 陳振舉. 東北南部蒙古櫟徑向生長對氣候變化的響應(yīng)——以千山為例. 應(yīng)用生態(tài)學(xué)報, 2014, 25(7): 1841- 1848.

    [37] 孫鳳華, 楊修群, 路爽, 楊素英. 東北地區(qū)平均、最高、最低氣溫時空變化特征及對比分析. 氣象科學(xué), 2006, 26(2): 157- 163.

    [38] 于大炮, 王順忠, 唐立娜, 代力民, 王慶禮, 王紹先. 長白山北坡落葉松年輪年表及其與氣候變化的關(guān)系. 應(yīng)用生態(tài)學(xué)報, 2005, 16(1): 14- 20.

    [39] 李廣起, 白帆, 桑衛(wèi)國. 長白山紅松和魚鱗云杉在分布上限的徑向生長對氣候變暖的不同響應(yīng). 植物生態(tài)學(xué)報, 2011, 35(5): 500- 511.

    [40] 王淼, 白淑菊, 陶大立, 單建平. 大氣增溫對長白山林木直徑生長的影響. 應(yīng)用生態(tài)學(xué)報, 1995, 6(2): 128- 132.

    [41] Driscoll W W, Wiles G C, D′Arrigo R D, Wilmking M. Divergent tree growth response to recent climatic warming, Lake Clark National Park and Preserve, Alaska. Geophysical Research Letters, 2005, 32(20): L20703, doi: 10.1029/2005GL024258.

    [42] 李宗善, 劉國華, 傅伯杰, 張齊兵, 胡嬋娟, 羅淑政. 川西臥龍國家級自然保護(hù)區(qū)樹木生長對氣候響應(yīng)的時間穩(wěn)定性評估. 植物生態(tài)學(xué)報, 2010, 34(9): 1045- 1057.

    [43] 尚華明, 魏文壽, 袁玉江, 喻樹龍, 張同文, 瓦合提·艾則買提, 李新建. 阿爾泰山南坡樹輪寬度對氣候變暖的響應(yīng). 生態(tài)學(xué)報, 2010, 30(9): 2246- 2253.

    Spatio-temporal variations in climate-growth relationships of three hardwood tree species across the north Zhangguangcai Mountains, northeast China

    SU Jinjuan, WANG Xiaochun*

    CenterforEcologicalResearch,NortheastForestryUniversity,Harbin150040,China

    Forests play a critical role in the global ecological balance. One of the major uncertainties in global climate change predictions is how the spatio-temporal dynamics of forests will be affected by global warming. Tree rings at different spatial and temporal scales can reflect the impacts of environmental change on tree growth, and can be critical for understanding the influence of climate change on forest ecosystems. Here, we investigated the spatio-temporal patterns of variation in climate-growth relationships of three common broadleaved tree species (Fraxinusmandshurica,Phellodendronamurense, andJuglansmandshurica) in a broadleaved Korean pine mixed forest across the north Zhangguangcai Mountains, northeast China. Tree-ring chronologies from three sample sites were established, transformed into principal component analysis (PCA), and analyzed through correlation analysis. With these analyses, climate factors limiting the radial growth were identified, as well as the detailed spatial radial growth-climate associations. The responses of the tree species at similar sites to climate factors differed significantly. Radial growth ofJ.mandshuricawas limited by monthly total precipitation and minimum temperature, whereasP.amurenseandF.mandshuricawere mainly affected by monthly mean minimum temperature. This suggests that climate-growth relationships are affected by species-specific characteristics. With changes in temperature and precipitation patterns, the three hardwoods differed in their spatial growth response to climate factors. The radial growth ofJ.mandshuricain Fangzheng and Xidaquan was positively correlated with minimum June to August temperatures (P<0.05). In contrast, it was positively correlated in Fenghuangshan with minimum temperature and precipitation in May and June (P<0.05). The radial growth ofP.amurenseandF.mandshuricadecreased from west to east with precipitation changes, since the correlation coefficients between tree-ring index and precipitation decreased from west to east. Temperatures in the north Zhangguangcai Mountains have increased since 1980, and the radial growth of the three hardwoods was consistent before 1980. Because climate factors affect growth differently, these boreal tree species showed inconsistent spatio-temporal responses to the recent warming, withJ.mandshuricabeing probably the most sensitive species. Radial growth ofP.amurenseandF.mandshuricaincreased with rapid warming since 1980, whereas tree growth ofJ.mandshuricadecreased after 1980. Therefore,J.mandshuricamay be subjected to drought stress and a subsequent decline in radial growth in the north Zhangguangcai Mountains if global warming continues or exacerbates. In contrast,P.amurenseandF.mandshuricamay benefit from increased temperatures through enhanced radial growth. For future research, the spatio-temporal variation of growth-climate relationship of the three hardwoods could be explored further by investigating the tree ring anatomical and physiological responses.

    tree rings; broadleaved Korean pine mixed forest; drought stress; global warming; radial growth

    國家自然科學(xué)基金面上項目(41471168, 31370463);教育部新世紀(jì)優(yōu)秀人才支持計劃(NCET- 12-0810);黑龍江省留學(xué)歸國科學(xué)基金(LC2012C09)共同資助

    2015- 09- 25;

    日期:2016- 07- 13

    10.5846/stxb201509251973

    *通訊作者Corresponding author.E-mail: wangxc-cf@nefu.edu.cn

    蘇金娟, 王曉春.張廣才嶺北部三大硬闊樹木生長-氣候關(guān)系的時空變異.生態(tài)學(xué)報,2017,37(5):1484- 1495.

    Su J J, Wang X C.Spatio-temporal variations in climate-growth relationships of three hardwood tree species across the north Zhangguangcai Mountains, northeast China.Acta Ecologica Sinica,2017,37(5):1484- 1495.

    猜你喜歡
    水曲柳最低溫度年表
    年表
    世界建筑(2022年10期)2022-10-25 08:26:48
    京藏高速民和至西寧段冬季路面與地面最低溫度變化特征及相關(guān)性研究
    錦州地區(qū)各預(yù)報產(chǎn)品溫度預(yù)報質(zhì)量分析
    錦州地區(qū)各預(yù)報產(chǎn)品溫度預(yù)報質(zhì)量分析
    寶龍店水曲柳無性系生長早期評價
    王錫良年表
    石魯年表
    淺談水曲柳營林的速產(chǎn)豐產(chǎn)技術(shù)要點
    水曲柳和落葉松人工純林與混交林的碳儲量
    森林工程(2018年4期)2018-08-04 03:23:10
    林區(qū)水曲柳營林的速生豐產(chǎn)技術(shù)要點
    午夜影院日韩av| 噜噜噜噜噜久久久久久91| 亚洲精品一区av在线观看| 欧美绝顶高潮抽搐喷水| 国产aⅴ精品一区二区三区波| 午夜福利成人在线免费观看| 久久久国产成人免费| 国产一区二区在线观看日韩| 日韩成人伦理影院| 午夜精品一区二区三区免费看| 色视频www国产| 韩国av在线不卡| 桃色一区二区三区在线观看| 老司机午夜福利在线观看视频| 日本在线视频免费播放| 国产精品亚洲美女久久久| 黄色日韩在线| 日日摸夜夜添夜夜添小说| 精品无人区乱码1区二区| 久久韩国三级中文字幕| 综合色丁香网| 欧美性猛交╳xxx乱大交人| 国产精品日韩av在线免费观看| 中文字幕熟女人妻在线| 狠狠狠狠99中文字幕| 六月丁香七月| 久久久精品欧美日韩精品| 亚洲成人av在线免费| 最近最新中文字幕大全电影3| 日本在线视频免费播放| 一级av片app| 午夜久久久久精精品| 欧美绝顶高潮抽搐喷水| 国产午夜福利久久久久久| 国产视频一区二区在线看| 亚洲在线观看片| 一级黄片播放器| 可以在线观看的亚洲视频| 成人性生交大片免费视频hd| 久久久久性生活片| 亚洲四区av| 国产黄色视频一区二区在线观看 | 国产一区二区亚洲精品在线观看| 久久精品国产自在天天线| 欧美日韩在线观看h| 蜜臀久久99精品久久宅男| 国产av在哪里看| 欧美最新免费一区二区三区| 联通29元200g的流量卡| 久久精品国产亚洲av天美| 亚洲不卡免费看| av黄色大香蕉| 人妻夜夜爽99麻豆av| 免费高清视频大片| 三级毛片av免费| 啦啦啦观看免费观看视频高清| 久久亚洲国产成人精品v| 成人鲁丝片一二三区免费| 久久午夜福利片| 噜噜噜噜噜久久久久久91| 精品午夜福利在线看| 亚洲性夜色夜夜综合| 51国产日韩欧美| 又黄又爽又刺激的免费视频.| 国产亚洲精品综合一区在线观看| 真实男女啪啪啪动态图| 又爽又黄无遮挡网站| 男人和女人高潮做爰伦理| av女优亚洲男人天堂| 高清毛片免费观看视频网站| 俺也久久电影网| 最近在线观看免费完整版| 欧美国产日韩亚洲一区| 国产在线男女| 色视频www国产| 禁无遮挡网站| 日韩成人av中文字幕在线观看 | 亚洲国产精品成人久久小说 | 午夜日韩欧美国产| 97人妻精品一区二区三区麻豆| 一个人免费在线观看电影| 欧美中文日本在线观看视频| 99久久九九国产精品国产免费| 亚洲精品在线观看二区| 日本色播在线视频| 一进一出抽搐gif免费好疼| 午夜影院日韩av| 亚洲av五月六月丁香网| 少妇高潮的动态图| 国内精品美女久久久久久| 亚洲精品日韩av片在线观看| 欧美激情久久久久久爽电影| 国产综合懂色| 最近2019中文字幕mv第一页| 亚洲精品国产av成人精品 | 色综合色国产| 国语自产精品视频在线第100页| 看片在线看免费视频| 啦啦啦韩国在线观看视频| 在线播放国产精品三级| 亚洲国产精品国产精品| 久久久久久久久久久丰满| 国产精品99久久久久久久久| 麻豆久久精品国产亚洲av| 欧美三级亚洲精品| 校园春色视频在线观看| 最后的刺客免费高清国语| 深夜精品福利| 少妇裸体淫交视频免费看高清| 国产精品久久久久久av不卡| 国产亚洲精品久久久com| 午夜亚洲福利在线播放| 亚洲国产欧美人成| 欧美一区二区国产精品久久精品| 亚洲无线在线观看| 一边摸一边抽搐一进一小说| 美女被艹到高潮喷水动态| 欧美又色又爽又黄视频| 亚洲高清免费不卡视频| 十八禁国产超污无遮挡网站| 97碰自拍视频| av视频在线观看入口| 国产在视频线在精品| 日本黄大片高清| 亚洲国产欧美人成| 欧美高清性xxxxhd video| 国产亚洲精品久久久久久毛片| 免费不卡的大黄色大毛片视频在线观看 | 国产美女午夜福利| 成人永久免费在线观看视频| 亚洲一级一片aⅴ在线观看| 99在线视频只有这里精品首页| 赤兔流量卡办理| 国产三级中文精品| 国模一区二区三区四区视频| 麻豆一二三区av精品| 天天躁夜夜躁狠狠久久av| www.色视频.com| 久久久久久伊人网av| 精品不卡国产一区二区三区| 欧美激情久久久久久爽电影| 全区人妻精品视频| 人人妻人人澡人人爽人人夜夜 | 十八禁网站免费在线| 久久午夜亚洲精品久久| 久久99热6这里只有精品| 午夜爱爱视频在线播放| 一本久久中文字幕| 久久韩国三级中文字幕| 精品欧美国产一区二区三| 一级黄片播放器| 真人做人爱边吃奶动态| 亚洲欧美精品自产自拍| 精华霜和精华液先用哪个| 我的老师免费观看完整版| .国产精品久久| 精品熟女少妇av免费看| 日韩人妻高清精品专区| 亚洲av.av天堂| 国模一区二区三区四区视频| 亚洲性久久影院| 黄色日韩在线| 小蜜桃在线观看免费完整版高清| 亚洲图色成人| 午夜福利成人在线免费观看| 非洲黑人性xxxx精品又粗又长| 春色校园在线视频观看| 在线免费十八禁| 国产激情偷乱视频一区二区| 麻豆久久精品国产亚洲av| 精品乱码久久久久久99久播| 在线观看美女被高潮喷水网站| 一本一本综合久久| 干丝袜人妻中文字幕| 精品人妻一区二区三区麻豆 | 精品久久久久久成人av| 精品日产1卡2卡| 国产色婷婷99| 中文字幕av在线有码专区| 亚洲国产高清在线一区二区三| 国产一区二区亚洲精品在线观看| 尤物成人国产欧美一区二区三区| 国产午夜精品久久久久久一区二区三区 | 99热全是精品| .国产精品久久| 97在线视频观看| 午夜精品国产一区二区电影 | 欧美激情在线99| 亚洲成人中文字幕在线播放| 此物有八面人人有两片| 九九久久精品国产亚洲av麻豆| 超碰av人人做人人爽久久| 在现免费观看毛片| 成人精品一区二区免费| 高清日韩中文字幕在线| 成年女人毛片免费观看观看9| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av成人精品一区久久| 亚洲国产精品国产精品| 国产伦在线观看视频一区| 成人av一区二区三区在线看| 精品一区二区免费观看| avwww免费| 国产精品免费一区二区三区在线| 国产欧美日韩精品亚洲av| 深爱激情五月婷婷| 性插视频无遮挡在线免费观看| 欧美另类亚洲清纯唯美| 久久精品91蜜桃| 美女xxoo啪啪120秒动态图| 秋霞在线观看毛片| 级片在线观看| 内地一区二区视频在线| 尤物成人国产欧美一区二区三区| 欧美最黄视频在线播放免费| 国产午夜精品论理片| 熟女电影av网| 搡女人真爽免费视频火全软件 | 日韩欧美一区二区三区在线观看| 久久亚洲精品不卡| 亚洲国产欧美人成| 性色avwww在线观看| 亚洲不卡免费看| av.在线天堂| 国产国拍精品亚洲av在线观看| 最新在线观看一区二区三区| 国产一级毛片七仙女欲春2| 听说在线观看完整版免费高清| 99久久成人亚洲精品观看| 神马国产精品三级电影在线观看| 男人和女人高潮做爰伦理| 小蜜桃在线观看免费完整版高清| 国产午夜福利久久久久久| 身体一侧抽搐| 国产一区二区三区av在线 | 国产亚洲精品av在线| 嫩草影视91久久| 国产成年人精品一区二区| 亚洲欧美中文字幕日韩二区| 亚洲人成网站高清观看| 亚洲人成网站在线播放欧美日韩| avwww免费| 日产精品乱码卡一卡2卡三| 国产精品无大码| 亚洲自偷自拍三级| 免费搜索国产男女视频| 校园春色视频在线观看| 看十八女毛片水多多多| 国产精品女同一区二区软件| 国产欧美日韩精品一区二区| 久久久久久久久久黄片| 国产精品永久免费网站| 天天躁夜夜躁狠狠久久av| 国产精品精品国产色婷婷| 非洲黑人性xxxx精品又粗又长| 免费黄网站久久成人精品| 女的被弄到高潮叫床怎么办| av黄色大香蕉| 不卡视频在线观看欧美| 最好的美女福利视频网| 丰满乱子伦码专区| 免费人成视频x8x8入口观看| 99热6这里只有精品| 九九爱精品视频在线观看| 亚洲在线观看片| 亚洲人成网站在线播放欧美日韩| avwww免费| 日日干狠狠操夜夜爽| 亚洲国产精品国产精品| eeuss影院久久| 狂野欧美白嫩少妇大欣赏| 国产熟女欧美一区二区| 久久人人爽人人爽人人片va| 老司机影院成人| 欧洲精品卡2卡3卡4卡5卡区| 看片在线看免费视频| 国产老妇女一区| 国产精品一及| 18禁在线无遮挡免费观看视频 | 日韩欧美精品v在线| 麻豆国产av国片精品| 国产精品1区2区在线观看.| 超碰av人人做人人爽久久| 免费人成在线观看视频色| 神马国产精品三级电影在线观看| 69人妻影院| 午夜福利成人在线免费观看| 女同久久另类99精品国产91| 免费高清视频大片| 99热全是精品| 国产精品野战在线观看| 色哟哟哟哟哟哟| 久久亚洲精品不卡| 激情 狠狠 欧美| 亚洲最大成人手机在线| 国内久久婷婷六月综合欲色啪| 欧美不卡视频在线免费观看| www.色视频.com| 亚洲va在线va天堂va国产| 日韩欧美在线乱码| 哪里可以看免费的av片| 中文字幕久久专区| 亚洲国产欧洲综合997久久,| 日本一本二区三区精品| 高清日韩中文字幕在线| 日韩欧美国产在线观看| 国产伦精品一区二区三区视频9| 午夜精品在线福利| 国产色婷婷99| 你懂的网址亚洲精品在线观看 | 国产亚洲91精品色在线| 热99re8久久精品国产| 伊人久久精品亚洲午夜| 久久这里只有精品中国| 国产精品久久久久久av不卡| 欧洲精品卡2卡3卡4卡5卡区| 寂寞人妻少妇视频99o| 国产成人a∨麻豆精品| 国产成人freesex在线 | av福利片在线观看| 精品久久久久久久末码| 国产精品1区2区在线观看.| 一区二区三区四区激情视频 | 久久国产乱子免费精品| 国产单亲对白刺激| 精华霜和精华液先用哪个| 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 精品人妻熟女av久视频| 大型黄色视频在线免费观看| 三级毛片av免费| 成人av一区二区三区在线看| 人妻少妇偷人精品九色| 大型黄色视频在线免费观看| 少妇人妻精品综合一区二区 | 五月伊人婷婷丁香| 午夜福利成人在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 如何舔出高潮| 深夜a级毛片| 美女高潮的动态| 亚洲精品亚洲一区二区| 美女 人体艺术 gogo| 日韩欧美 国产精品| 久久九九热精品免费| 国产精品亚洲一级av第二区| 一级毛片我不卡| 三级经典国产精品| 国模一区二区三区四区视频| 亚洲人成网站在线播放欧美日韩| 免费观看的影片在线观看| 国产高清不卡午夜福利| 亚洲精品国产av成人精品 | 午夜激情欧美在线| av天堂在线播放| 中文字幕精品亚洲无线码一区| 国产高潮美女av| 99久国产av精品国产电影| 欧美三级亚洲精品| 99在线视频只有这里精品首页| 亚洲一区二区三区色噜噜| 99久国产av精品国产电影| 你懂的网址亚洲精品在线观看 | 舔av片在线| 啦啦啦观看免费观看视频高清| 婷婷色综合大香蕉| 国模一区二区三区四区视频| 成人欧美大片| 免费观看在线日韩| 精品国产三级普通话版| 俺也久久电影网| 精品久久久噜噜| 蜜臀久久99精品久久宅男| 免费看美女性在线毛片视频| 国内精品久久久久精免费| 男女下面进入的视频免费午夜| 如何舔出高潮| 国产精品av视频在线免费观看| 精品99又大又爽又粗少妇毛片| 成人三级黄色视频| 久99久视频精品免费| 婷婷亚洲欧美| 日韩欧美精品v在线| 免费高清视频大片| 日日摸夜夜添夜夜爱| 99热这里只有精品一区| 日韩在线高清观看一区二区三区| 成人性生交大片免费视频hd| 91午夜精品亚洲一区二区三区| 日韩 亚洲 欧美在线| 国产三级中文精品| 干丝袜人妻中文字幕| 97人妻精品一区二区三区麻豆| 特大巨黑吊av在线直播| 男女啪啪激烈高潮av片| 午夜福利视频1000在线观看| 国产精品99久久久久久久久| 高清日韩中文字幕在线| 欧美xxxx性猛交bbbb| 免费av观看视频| 男女做爰动态图高潮gif福利片| 日本一二三区视频观看| 欧美潮喷喷水| 日本 av在线| 美女高潮的动态| 我的老师免费观看完整版| 午夜免费激情av| 久久精品91蜜桃| 亚洲aⅴ乱码一区二区在线播放| 国产伦精品一区二区三区四那| 亚洲乱码一区二区免费版| 99国产精品一区二区蜜桃av| 91精品国产九色| 一个人看的www免费观看视频| 国产一区亚洲一区在线观看| 国产视频内射| 欧美日韩综合久久久久久| 18禁在线播放成人免费| 国产成人a区在线观看| 国内揄拍国产精品人妻在线| 国产在视频线在精品| 五月玫瑰六月丁香| 高清毛片免费看| 欧美日韩在线观看h| 国内精品宾馆在线| 久久久精品大字幕| 亚洲欧美精品自产自拍| 最近视频中文字幕2019在线8| 男人狂女人下面高潮的视频| av在线蜜桃| 成人午夜高清在线视频| 精品久久久噜噜| 精品一区二区免费观看| 夜夜夜夜夜久久久久| 日本与韩国留学比较| 亚洲欧美精品综合久久99| 小说图片视频综合网站| 国产精品国产高清国产av| 丰满乱子伦码专区| 日本一本二区三区精品| 成人国产麻豆网| 成人鲁丝片一二三区免费| 神马国产精品三级电影在线观看| 99久国产av精品| 午夜a级毛片| 免费看日本二区| 国产av在哪里看| 国产欧美日韩一区二区精品| 欧美激情久久久久久爽电影| 国产精品久久久久久精品电影| 久久久久性生活片| 欧美最新免费一区二区三区| a级毛片免费高清观看在线播放| 成年女人毛片免费观看观看9| 在线观看av片永久免费下载| 在现免费观看毛片| 国产白丝娇喘喷水9色精品| 欧美性猛交╳xxx乱大交人| 淫妇啪啪啪对白视频| 国产精品一及| 亚洲人成网站在线播| 自拍偷自拍亚洲精品老妇| 俺也久久电影网| 日韩,欧美,国产一区二区三区 | 美女高潮的动态| 永久网站在线| 天天躁日日操中文字幕| 国产探花极品一区二区| 久久久久久大精品| 免费一级毛片在线播放高清视频| 精品国产三级普通话版| 热99re8久久精品国产| 波野结衣二区三区在线| 99九九线精品视频在线观看视频| 女的被弄到高潮叫床怎么办| 天天躁夜夜躁狠狠久久av| 毛片一级片免费看久久久久| 看黄色毛片网站| 亚洲熟妇中文字幕五十中出| 搡老熟女国产l中国老女人| 丝袜美腿在线中文| 色哟哟哟哟哟哟| 亚洲成a人片在线一区二区| 亚洲欧美成人综合另类久久久 | 免费看av在线观看网站| 日本-黄色视频高清免费观看| 九色成人免费人妻av| 干丝袜人妻中文字幕| 黄色日韩在线| 国产蜜桃级精品一区二区三区| 久久99热6这里只有精品| 亚洲中文字幕一区二区三区有码在线看| 毛片一级片免费看久久久久| 亚洲不卡免费看| 久久精品国产自在天天线| 精品午夜福利视频在线观看一区| 高清日韩中文字幕在线| 午夜日韩欧美国产| 六月丁香七月| 久久精品人妻少妇| 国产91av在线免费观看| 亚洲熟妇熟女久久| 久久久久国产网址| 一本一本综合久久| 亚洲国产精品sss在线观看| 精品人妻一区二区三区麻豆 | 欧美国产日韩亚洲一区| 久久久久久国产a免费观看| 国内精品一区二区在线观看| 亚洲欧美日韩高清专用| 日本色播在线视频| 亚洲在线观看片| 我要搜黄色片| 久久精品91蜜桃| 国产亚洲欧美98| av免费在线看不卡| 人妻丰满熟妇av一区二区三区| 深夜精品福利| 国产av一区在线观看免费| 欧美日本亚洲视频在线播放| 国产成人影院久久av| 看十八女毛片水多多多| 男女之事视频高清在线观看| 午夜日韩欧美国产| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 97超视频在线观看视频| 天堂网av新在线| 国产亚洲91精品色在线| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久av不卡| 日韩av不卡免费在线播放| 亚洲美女搞黄在线观看 | 亚洲精华国产精华液的使用体验 | 天堂影院成人在线观看| 床上黄色一级片| 日韩一区二区视频免费看| 亚洲自偷自拍三级| 一本精品99久久精品77| 草草在线视频免费看| 女人被狂操c到高潮| 99国产极品粉嫩在线观看| 老司机影院成人| 欧美三级亚洲精品| 日本欧美国产在线视频| 欧美一区二区精品小视频在线| 给我免费播放毛片高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产高清有码在线观看视频| 人妻制服诱惑在线中文字幕| 在线观看午夜福利视频| 自拍偷自拍亚洲精品老妇| 久久人妻av系列| 免费观看人在逋| 亚洲一区二区三区色噜噜| 国产精品美女特级片免费视频播放器| 亚洲精品在线观看二区| 日本黄色片子视频| 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 插阴视频在线观看视频| 亚洲经典国产精华液单| 国产黄色小视频在线观看| 亚洲精华国产精华液的使用体验 | av中文乱码字幕在线| 国产片特级美女逼逼视频| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 成熟少妇高潮喷水视频| 一级毛片aaaaaa免费看小| 我要搜黄色片| 波野结衣二区三区在线| 国产亚洲av嫩草精品影院| 欧美一区二区精品小视频在线| 国产成人aa在线观看| av中文乱码字幕在线| 亚洲久久久久久中文字幕| 精品久久久久久久久亚洲| 亚洲天堂国产精品一区在线| 国产一区二区在线av高清观看| 国产成人一区二区在线| 欧美zozozo另类| 国产午夜福利久久久久久| 国产中年淑女户外野战色| 精品人妻视频免费看| 欧美高清性xxxxhd video| 深夜精品福利| 少妇的逼水好多| 一夜夜www| 一级黄色大片毛片| 在现免费观看毛片| 人人妻人人澡欧美一区二区| 九九爱精品视频在线观看| 久久久精品欧美日韩精品| 午夜福利在线观看免费完整高清在 | 国产精品久久视频播放| 国产精品福利在线免费观看| 人妻丰满熟妇av一区二区三区| 毛片一级片免费看久久久久| 成人av一区二区三区在线看| 国产v大片淫在线免费观看| 国产精品日韩av在线免费观看| 国产熟女欧美一区二区| 老司机福利观看| 亚洲最大成人av| 三级国产精品欧美在线观看| 日本在线视频免费播放| 少妇裸体淫交视频免费看高清| 久久久久久九九精品二区国产| 中文在线观看免费www的网站| а√天堂www在线а√下载| 久久热精品热| 国产在线精品亚洲第一网站| 一个人观看的视频www高清免费观看|