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

    模擬氣候變暖對東北興安落葉松徑向生長的影響

    2017-06-26 11:28:08劉玉佳顧卓欣王曉春
    生態(tài)學(xué)報(bào) 2017年8期
    關(guān)鍵詞:松嶺小枝興安

    劉玉佳,顧卓欣,王曉春

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

    ?

    模擬氣候變暖對東北興安落葉松徑向生長的影響

    劉玉佳,顧卓欣,王曉春*

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

    2004年在帶嶺、孫吳、松嶺和塔河構(gòu)成的緯度梯度上分別選取幼齡興安落葉松(Larixgmelinii),將其統(tǒng)一移栽到緯度較低的帽兒山地區(qū),以模擬溫度升高對興安落葉松徑向生長的影響。在2013年生長季(6—9月)分別在4個(gè)緯度移栽原地和移栽地(帽兒山)的樹木主干和小枝上取樣,主干用微生長錐鉆取樣芯,小枝直接截取一段,主干和小枝的樣品用石蠟切片法得到擴(kuò)大細(xì)胞和成熟細(xì)胞的大小及數(shù)量的動(dòng)態(tài)變化。結(jié)果表明:溫度升高導(dǎo)致興安落葉松主干擴(kuò)大細(xì)胞大小和數(shù)量均增加,但是主干徑向生長與溫度升高不呈正比。生長季前期(6月),帶嶺小枝成熟細(xì)胞數(shù)量明顯高于其他兩組,但在8月晚材細(xì)胞數(shù)量偏少,表明其還在繼續(xù)生長。生長季后期(8月),松嶺興安落葉松小枝擴(kuò)大細(xì)胞平均大小為5 μm,分化數(shù)為4個(gè),而松嶺移栽到帽兒山后擴(kuò)大細(xì)胞平均大小為9.9 μm,分化數(shù)為7個(gè),增長了近1倍。松嶺成熟細(xì)胞大小較其他兩地偏大,可能與其水分狀況有關(guān)。溫度升高不能直接改變植物生長周期一類由基因控制的生理特性,但是對樹木不同構(gòu)件的影響力是不同的,小枝受到的影響要大于主干。當(dāng)全球變暖尤其是伴隨降水量減少時(shí),將減少興安落葉松徑向生長,包括管胞直徑、數(shù)量和生長周期。通過緯度變化模擬氣候變暖,深入分析不同構(gòu)件樹木生長對氣候變化的響應(yīng),為預(yù)測氣候變化對興安落葉松生長及分布的影響提供依據(jù)。

    興安落葉松;年輪;暖化;形成層;切片

    近幾十年以來,人類活動(dòng)引起溫室氣體排量增加,致使地表溫度升高,導(dǎo)致全球持續(xù)變暖[1],模型預(yù)測未來升溫可能更為明顯[2]。預(yù)計(jì)到2100年,全球地表溫度將相對1980—1999年的水平升高1.1—6.4℃[3],而我國溫度將升高3.9—6.0℃[4]。北部高緯度預(yù)計(jì)將有更大幅度的變暖[5-6],全球升溫對高緯度地區(qū)影響顯著[1,7],這可能導(dǎo)致高緯度寒溫帶森林面積大幅度下降[7-8]。興安落葉松(Larixgmelinii)林是我國北方森林的優(yōu)勢植被類型[9],在大興安嶺地區(qū)占到70%以上[10],其碳儲(chǔ)量約占全國森林總碳儲(chǔ)量的8.5%[11]。據(jù)張新時(shí)[12]、程肖俠等[13]預(yù)測,溫度、降水量增加到一定程度時(shí),興安落葉松林可能向北大幅移動(dòng)。因此,探索興安落葉松對全球變暖的響應(yīng)具有重要意義。

    在樹木徑向生長與氣候變化的響應(yīng)關(guān)系上,溫度一直是一個(gè)主要研究因子。有些樹木年輪寬度與生長季前后的溫度呈負(fù)相關(guān),即溫度越高徑向生長量越少[14],還有一些樹木剛好與之相反[15]。這樣當(dāng)溫度變化時(shí),不同樹木的響應(yīng)就會(huì)出現(xiàn)差異[16-17]。對不同時(shí)間段的溫度變化,徑向生長也有不同響應(yīng),5—8月的溫度升高可以促進(jìn)山楊(Populusdavidiana)的徑向生長[18]。在樹木徑向生長過程中,管胞的一些形態(tài)特征同年輪寬度一樣也傳遞著重要?dú)夂蛐畔19]。管腔尺寸和氣候波動(dòng)有緊密的聯(lián)系[20],2—3月的平均溫度升高將有利于管胞直徑的變大[21]。管胞直徑、數(shù)量、管胞壁厚等一些形態(tài)特征存在著較大的差異性[22-23],在氣候變化后,管胞的數(shù)量、大小、直徑、壁厚等特征不可能同時(shí)增加[24]。在每年早材形成的5—7月,溫度對其影響較高。由于早材細(xì)胞形成于春季,春季溫度的適度提升有利于形成層活動(dòng),可形成較大、較多的管胞。入秋后,溫度下降,形成層活動(dòng)趨于停止,形成的管胞較小、較少[25-26]。在對北美云杉(Piceasitchensis)的研究中發(fā)現(xiàn),管胞數(shù)量和直徑與日平均溫度有著明顯的相關(guān)關(guān)系[27]。在徑向生長的形成層活動(dòng)中,形成層活動(dòng)對于溫度的響應(yīng)存在一定的差異,在植物休眠階段通過局部加熱發(fā)現(xiàn),形成層活化區(qū)僅為加熱區(qū)域[28],形成層開始活動(dòng)后逐漸分化出早材細(xì)胞,如果降低溫度,對形成層影響較小,主要是晚材比例下降[29]。

    本研究通過將不同緯度梯度的興安落葉松幼苗,移植到其分布區(qū)南緣的帽兒山森林生態(tài)系統(tǒng)國家野外觀測研究站,模擬升溫對興安落葉松的徑向生長(包括年輪寬度、莖和枝的木質(zhì)部管胞特征及形成層分化情況)的影響,以探討:(1)不同程度的溫度升高對興安落葉松徑向生長有什么樣的影響,影響是否有差別?(2)不同程度的升溫給興安落葉松徑向生長帶來的影響在不同部位是否有差別?(3)興安落葉松移栽前后的同一時(shí)期,徑向生長是否存在差別?

    1 材料與方法

    1.1 研究區(qū)域概況

    本研究將實(shí)驗(yàn)地分為兩個(gè)部分,一部分為移栽原地,另一部分為移栽地。在興安落葉松自然分布區(qū)內(nèi)由南向北選擇4個(gè)實(shí)驗(yàn)點(diǎn):帶嶺、孫吳、松嶺和塔河,構(gòu)建1條主要由溫度驅(qū)動(dòng)的緯度梯度樣帶(表1)。于2004年秋季,根據(jù)其起源、樹齡和經(jīng)營方式等選取上述4個(gè)地點(diǎn)的興安落葉松幼苗、地被物和土壤A、B層,移植到其分布區(qū)南緣的帽兒山森林生態(tài)系統(tǒng)國家野外觀測研究站實(shí)驗(yàn)地,構(gòu)成帶嶺、孫吳、松嶺和塔河4個(gè)實(shí)驗(yàn)處理[29]。

    實(shí)驗(yàn)地采用隨機(jī)設(shè)計(jì),每個(gè)處理設(shè)置4塊重復(fù)樣地,每塊樣地面積為10 m×10 m,移栽木行距為1.5 m[29]。

    1.2 樣品采集

    2013年生長季(6—9月),在移栽原地(帶嶺、孫吳、松嶺和塔河)和移栽地(帽兒山)分別隨機(jī)選取3棵樣樹,用微生長錐在每棵樣樹上取3個(gè)樣芯并放入福爾馬林-乙酸-乙醇(FAA)固定液中固定。每棵樣樹上再分別取3個(gè)5 cm長小枝同樣放入FAA固定液中固定。固定好的樣品后帶回實(shí)驗(yàn)室后將樣芯放入70%乙醇溶液中,隨后再以5%乙醇濃度遞增直至100%的乙醇溶液脫水浸蠟包埋,隨后對包埋好樹芯進(jìn)行石蠟切片,樣芯切片后對木質(zhì)化和非木質(zhì)化部分主要用染色劑來區(qū)分,番紅可以將木質(zhì)化部分染成紅色,固綠將非木質(zhì)化部分染成藍(lán)綠色。用光學(xué)顯微鏡觀察切片封好的片子,通過測量細(xì)胞形成層擴(kuò)大細(xì)胞的尺寸,結(jié)合樣地處理后的氣候數(shù)據(jù)進(jìn)行分析。在時(shí)間尺度上通過觀察分化細(xì)胞的變化,了解木質(zhì)部形成的季節(jié)動(dòng)態(tài)變化。在生長季結(jié)束時(shí),用生長錐在塔河、松嶺、孫吳和帶嶺的移栽地(帽兒山)取樣,將樹芯帶回室內(nèi)實(shí)驗(yàn)室后打磨,在光學(xué)顯微鏡下觀察成熟管胞的徑向長度(以下簡稱管胞直徑)及數(shù)量。

    表1 興安落葉松移栽原地及移栽地的生境基本特征[29]

    1.3 統(tǒng)計(jì)分析

    采用單因素方差分析統(tǒng)計(jì)了2013年移栽原地和移栽地樹干和小枝之間擴(kuò)大細(xì)胞和成熟細(xì)胞大小、數(shù)量以及管胞直徑、數(shù)量,用多重比較來分析每個(gè)移栽原地之間的差別和移栽后的差別。

    2 結(jié)果與分析

    2.1 移栽原地自然溫度梯度下興安落葉松徑向生長的差異

    溫度升高對興安落葉松徑向生長有一定促進(jìn)作用(圖1)。4個(gè)移栽原地構(gòu)成的緯度梯度中,緯度較低的帶嶺和孫吳的平均年輪寬度顯著高于緯度較高的松嶺和塔河。從時(shí)間序列上看,緯度最高的塔河每年徑向生長量波動(dòng)不大,每年徑向生長量基本上在4個(gè)移栽原地中都是最低的。緯度最低的帶嶺每年的徑向生長量都相對較高,而孫吳在個(gè)別年份徑向生長量突出。

    在2013年生長季末,將當(dāng)年的年輪樹芯打磨光滑后在顯微鏡下測得管胞直徑并且按著由早材向晚材方向排列(圖2)。在生長季初期,4個(gè)移栽原地中塔河的管胞直徑明顯低于其它3個(gè)移栽原地,帶嶺、孫吳和松嶺3個(gè)移栽原地之間沒有明顯差別。在生長季后期的晚材管胞直徑上,可以看出帶嶺管胞的直徑要高于其它移栽原地的管胞直徑,塔河移栽原地依然是各原地中最低的。由此,管胞的直徑與緯度相反,緯度越高,興安落葉松管胞直徑越小。從管胞數(shù)量看,帶嶺在生長季階段數(shù)量最多,然后隨緯度的升高而遞減,松嶺的管胞數(shù)量比孫吳略低,緯度最高的塔河管胞數(shù)量是4個(gè)移栽原地中最少的,生長期也是最短的。溫度差異導(dǎo)致帶嶺興安落葉松管胞直徑和數(shù)量均最高,塔河的均最低。

    圖1 帶嶺、孫吳、松嶺和塔河移栽原地年輪寬度比較 Fig.1 Comparison of tree-ring width of Larix gmelinii in Dailing, Sunwu, Songling and Tahe

    圖2 帶嶺、孫吳、松嶺和塔河移栽原地徑向管胞直徑從早材到晚材變化比較Fig.2 Comparison of radial tracheid length from the first earlywood cell to the last latewood cell at Dailing, Sunwu, Songling and Tahe

    2.2 模擬升溫對興安落葉松主干形成層擴(kuò)大細(xì)胞的影響

    圖3 在生長季不同階段移栽地和原地微生長芯擴(kuò)大細(xì)胞數(shù)量及大小比較Fig.3 Comparison of the expanded cell number and size at different time of growing season for original and transplanting sites

    帽兒山地區(qū)移栽木從分化細(xì)胞的數(shù)量上看(圖3),6月初,孫吳擴(kuò)大細(xì)胞數(shù)量相對較少,擴(kuò)大細(xì)胞大小略大。整體看,3個(gè)地點(diǎn)移栽木在擴(kuò)大細(xì)胞大小上沒有顯著差別,而擴(kuò)大細(xì)胞數(shù)量則是松嶺移栽木顯著高于孫吳移栽木。6月中旬,3組移栽木的擴(kuò)大細(xì)胞數(shù)量沒有顯著差別,擴(kuò)大細(xì)胞大小則是孫吳顯著高于松嶺,并且達(dá)到擴(kuò)大細(xì)胞大小的峰值。7月中旬,3組移栽木的擴(kuò)大細(xì)胞數(shù)量均達(dá)到峰值,帶嶺顯著高于其它兩組,孫吳擴(kuò)大細(xì)胞數(shù)量最少。帶嶺和松嶺移栽木的擴(kuò)大細(xì)胞大小也在7月達(dá)到峰值,均顯著高于孫吳擴(kuò)大細(xì)胞。8月下旬,移栽木擴(kuò)大細(xì)胞數(shù)量與大小的趨勢相同,帶嶺移栽木擴(kuò)大細(xì)胞數(shù)量和大小均為最高,顯著高于其它地區(qū)的移栽木(圖3)。8月份的這一趨勢也與移栽木之前原地的緯度相一致,即原地緯度越高停止生長越早。綜合移栽木擴(kuò)大細(xì)胞大小和數(shù)量分析,帶嶺移栽木在7月出現(xiàn)峰值,且均高于其它兩地的移栽木;孫吳移栽木變化曲線幅度不大;松嶺移栽木則是在生長季階段細(xì)胞分化最早的移栽木,但是在8月和孫吳移栽木開始一起逐漸停止生長,時(shí)間上要早于帶嶺的移栽木(圖3)。

    由移栽原地與移栽地興安落葉松形成層擴(kuò)大細(xì)胞活動(dòng)可知,松嶺地區(qū)興安落葉松與移栽地一樣在6月初期也開始分化細(xì)胞開始生長,而在8月份的時(shí)候逐漸停止分化;雖然這一習(xí)性整體上并沒有因?yàn)闇囟鹊纳叨淖?但是在細(xì)節(jié)上二者存在差異較大。松嶺原地在8月份時(shí),分化的擴(kuò)大細(xì)胞平均大小為5 μm,分化數(shù)量為4。而移栽地在8月份時(shí)擴(kuò)大細(xì)胞平均大小為9.9 μm,分化數(shù)量為7,增長了近一倍;說明溫度升高雖然不能改變松嶺地區(qū)興安落葉松的生活習(xí)性,但是可以影響各個(gè)發(fā)育階段的生長幅度。8月份原地取樣,擴(kuò)大細(xì)胞均值最小的為松嶺,主要原因是切片中松嶺幾乎全部都是接近停止生長的狀態(tài),而其次的孫吳出現(xiàn)了部分接近停止分化的樣本。擴(kuò)大細(xì)胞在數(shù)量上松嶺也是移栽原地樣樹中最少的(表2)。

    表2 原地與移栽地樹干擴(kuò)大細(xì)胞數(shù)量和大小比較

    *表示移栽地和原地差異顯著(P< 0.05);M:移栽地Transplanting site,Y:原地Original sites

    2.3 模擬升溫對興安落葉松小枝形成層擴(kuò)大與成熟細(xì)胞的影響

    從不同緯度移栽至帽兒山的興安落葉松小枝(2—3年生)細(xì)胞大小可以看出(圖4),6月初,松嶺移栽木的擴(kuò)大細(xì)胞大小略大于其它兩地的移栽木,但不顯著。6月中旬,松嶺移栽木的擴(kuò)大細(xì)胞大小顯著高于其它兩地的移栽木。7月份,3組移栽木的擴(kuò)大細(xì)胞大小差異不顯著,孫吳移栽木擴(kuò)大細(xì)胞的大小達(dá)到峰值。8月份,孫吳移栽木的擴(kuò)大細(xì)胞大小顯著低于其它兩地的移栽木。帶嶺和松嶺達(dá)到擴(kuò)大細(xì)胞大小的峰值。9月份3地的移栽木差異不顯著,但是孫吳還是略低的。孫吳移栽木的生長時(shí)期相對較短,帶嶺的增長速率則是最高的。

    從擴(kuò)大細(xì)胞數(shù)量上看(圖4),6月份,孫吳移栽木顯著小于其它地點(diǎn)的移栽木。7月份,3個(gè)地點(diǎn)的移栽木細(xì)胞數(shù)量沒有顯著差異,均達(dá)到峰值。8月份孫吳移栽木的細(xì)胞數(shù)量顯著低于其它兩地的移栽木。9月份,3個(gè)地點(diǎn)的移栽木在小枝擴(kuò)大細(xì)胞數(shù)量上沒有顯著差異。綜合不同地點(diǎn)移栽木小枝擴(kuò)大細(xì)胞大小和數(shù)量比較,可以得出孫吳小枝的生長時(shí)間最短,松嶺略低于帶嶺。

    從移栽木小枝成熟細(xì)胞大小上看(圖4),松嶺移栽木的成熟細(xì)胞明顯大于其它兩地,6月份的成熟細(xì)胞中全部都是早材細(xì)胞,早材多與水分相關(guān),水分充足時(shí)管胞直徑較大,因?yàn)樗蓭X相比其它兩地濕度差最大,導(dǎo)致早材形成的管胞較大。8月份,松嶺移栽木成熟細(xì)胞大小明顯低于其它兩地。8月份已經(jīng)開始出現(xiàn)晚材細(xì)胞,這樣成熟細(xì)胞大小的均值被拉低,晚材細(xì)胞數(shù)量要高于其它兩地,接近停止生長,導(dǎo)致細(xì)胞大小均值上顯著低于其它地區(qū)。

    從移栽木小枝成熟細(xì)胞數(shù)量上看(圖4),帶嶺地區(qū)要高于其它兩組,說明在小枝活動(dòng)中帶嶺移栽木優(yōu)先開始生長活動(dòng)。在隨后月份中,數(shù)量都是逐漸增加的,但是增加速度最快的是松嶺移栽木,它在6月份活動(dòng)相對較晚,但是從7月份開始就在數(shù)量上領(lǐng)先,并在8月份顯著高于其它組。結(jié)合成熟細(xì)胞大小,可以發(fā)現(xiàn)松嶺移栽木在8月份時(shí)數(shù)量明顯超過帶嶺,但是在平均大小上卻略小。所以在8月份時(shí)孫吳和松嶺移栽木都出現(xiàn)大量晚材細(xì)胞,相比較而言,帶嶺移栽木晚材細(xì)胞最少,生長時(shí)期更長。

    圖4 移栽地興安落葉松小枝擴(kuò)大、成熟細(xì)胞大小及數(shù)量比較Fig.4 Comparisons of the enlarged and mature cell size and number at braches of Larix gmelinii from different transplanting sites

    6月份的小枝在移栽原地中,帶嶺的擴(kuò)大細(xì)胞大小明顯高于其它移栽原地(圖5),數(shù)量上帶嶺和松嶺明顯高于塔河(圖5)。而6月份的各個(gè)移栽原地中成熟細(xì)胞的大小無顯著差異(圖5),數(shù)量上看,帶嶺最高,塔河最低(圖5)。在生長季初期,興安落葉松的小枝徑向生長量與緯度分布相吻合,即溫度在生長季初期小枝的徑向生長中是關(guān)鍵因子。在不同緯度的移栽原地中,小枝成熟細(xì)胞的大小都沒有顯著差異,緯度較低的帶嶺地區(qū),成熟細(xì)胞的數(shù)量在生長季初期、中期和晚期都顯著高于其它移栽原地。此外,帶嶺地區(qū)的擴(kuò)大細(xì)胞大小和數(shù)量在生長季中期明顯高于其它移栽原地。

    從移栽地和原地小枝細(xì)胞數(shù)量和大小比較可以看出(表3),溫度升高有利于興安落葉松小枝的徑向生長。但是移栽地中帶嶺移栽木的小枝并沒有因?yàn)闇囟壬叨黾蛹?xì)胞數(shù)量和大小,相反在生長季初期,帶嶺成熟細(xì)胞大小上甚至明顯低于原地和其它地區(qū)的移栽木。在生長季末,帶嶺移栽木小枝成熟細(xì)胞數(shù)量比原地減少了近1/3。在松嶺和帶嶺地區(qū)在生長季初的6月,其擴(kuò)大細(xì)胞的數(shù)量和大小以及成熟細(xì)胞的數(shù)量和大小沒有差異,表明與移栽后的小枝在帽兒山有共同的生長軌跡,以表示對移栽地的適應(yīng)。

    圖5 原地興安落葉松小枝擴(kuò)大、成熟細(xì)胞大小及數(shù)量比較Fig.5 Comparisons of the enlarged and mature cell size and number at braches of Larix gmelinii from original sites

    表3 原地與移栽地小枝擴(kuò)大、成熟細(xì)胞數(shù)量和大小比較

    *表示移栽地和原地差異顯著(P< 0.05);M:移栽地Transplanting site,Y:原地Original sites

    3 討論

    樹木年輪是研究樹木生長與氣候關(guān)系的重要途徑,具有定年準(zhǔn)確、連續(xù)性強(qiáng)、分辨率高、易于獲取和復(fù)本量大等優(yōu)點(diǎn),已被廣泛用于全球變化以及古氣候研究[15- 26,30- 33]。樹輪解剖特征是反映樹木生長的一種更高分辨率代用指標(biāo),也是樹輪在細(xì)胞尺度上的分量。已被廣泛用于全球變化研究[19- 24,34-35]。近年來,東北地區(qū)溫度顯著升高[30,31]。在寒冷環(huán)境中生長的樹木,木質(zhì)部是由溫度主導(dǎo)代謝活動(dòng)[32],因此溫度變化影響和改變樹木的年輪特征[33]。隨緯度梯度變化的主要?dú)夂蛞蜃邮菧囟萚35],不同緯度移栽模擬了溫度升高對興安落葉松徑向生長的影響。移栽后,孫吳的年輪寬度超過低緯度的帶嶺,表明增溫幅度不同對興安落葉松徑向生長的影響也存在差異。樹輪對環(huán)境變化的記錄上除了寬度外,樹輪解剖結(jié)構(gòu)也記錄不少的生態(tài)環(huán)境信息[34, 36]。樹輪解剖結(jié)構(gòu)相對于樹輪寬度長期對氣候的響應(yīng),更能展現(xiàn)短期內(nèi)氣候變化對樹木徑向生長的影響。2013年興安落葉松的管胞大小及數(shù)量在不同緯度地區(qū)有明顯的區(qū)分,帶嶺管胞直徑相對較大,管胞數(shù)量明顯多于其它地區(qū),這導(dǎo)致帶嶺的徑向生長量較大。塔河年均溫度比帶嶺低約6℃,年均降水量也是4個(gè)移栽原地中最少的,其管胞大小和樹輪明顯低于其它地區(qū),導(dǎo)致塔河地區(qū)的徑向增量最小。

    樹木生長發(fā)育過程受到外部因素和內(nèi)部因素的控制,外部因素主要有溫度、降水和日照等。物候觀測、遙感數(shù)據(jù)等證實(shí)了溫度升高提早樹木展葉期、延遲落葉期,進(jìn)而延長了植物的生長期[37-38],而生長期的延長會(huì)顯著影響樹木的徑向生長[39]。高緯度地區(qū)樹木形成層的活動(dòng)主要依賴夏季溫度,在低溫時(shí)期提升溫度,可以促進(jìn)形成層開始活動(dòng)[28]。移栽后,松嶺興安落葉松樹干生長時(shí)期延長,導(dǎo)致生長期內(nèi)形成層分化的細(xì)胞數(shù)量和大小增加,進(jìn)而使徑向生長量增加,Vaganov等[40]也發(fā)現(xiàn)北方樹線附近的落葉松管胞直徑會(huì)因長期的夏季升溫而顯著增加。移栽原地中,帶嶺緯度最低,在6月份,移栽后樹干擴(kuò)大細(xì)胞大小并未顯著增加,細(xì)胞數(shù)量也有減少趨勢。這表明溫度在一定閾值內(nèi)增加有利于樹木的徑向生長,超過這個(gè)閾值繼續(xù)增加可能減少樹木徑向[41]。不同緯度樣地落葉松南移(升溫)后,徑向生長上表現(xiàn)一定程度上增加趨勢,但模擬升溫的強(qiáng)度與相應(yīng)的生長量或形成層活動(dòng)并不成正比,這可能是移栽(升溫)后控制落葉松徑向生長的限制因子發(fā)生轉(zhuǎn)變的結(jié)果[19,42];不同緯度興安落葉松南移后,溫度對落葉松的限制作用相對減弱,其他的環(huán)境因子如升溫所引起的暖干旱可能會(huì)使得水分取代溫度成為新的主控因子。已有研究表明,東北升溫在一定程度上促進(jìn)了樹木徑向生長,但是近幾十年徑向增長的速度有所減緩,這可能與溫度的閾值效應(yīng)以及升溫帶來的干旱脅迫有關(guān)[30,43]。高緯度地區(qū),樹木生長(或新細(xì)胞產(chǎn)生)約于8月中旬停止,生長季長度主要由形成層活動(dòng)的開始日期決定[44];相反,本研究則發(fā)現(xiàn)移栽地(升溫)與原地興安落葉松莖干形成層擴(kuò)大細(xì)胞的數(shù)量及大小,在生長季初期差異不顯著,而生長季末期差異顯著,表明升溫可能會(huì)以推遲生長季結(jié)束時(shí)間的方式延長生長季。

    研究不同功能型植物的小枝生長特征,可以更好地理解植物功能性狀隨環(huán)境變化發(fā)生的一些改變[45]。在生理上許多復(fù)雜的器官有單獨(dú)的代謝機(jī)能,小枝可以看作是樹木生長的獨(dú)立器官,具有一定的獨(dú)立性[46],是樹木最有活力的部分之一[47]。小枝的次級組織較少,在內(nèi)部能很快的進(jìn)行營養(yǎng)物質(zhì)和碳水化合物流動(dòng)和轉(zhuǎn)移,相比整株樹木對環(huán)境的響應(yīng)更顯著或者具有與主干不同的響應(yīng)特征[48]。移栽前,在生長季初期,原地小枝擴(kuò)大細(xì)胞大小與溫度呈正比,即隨緯度升高,溫度降低,擴(kuò)大細(xì)胞大小降低;在生長季末期,小枝成熟細(xì)胞數(shù)量也與溫度成正比,帶嶺成熟細(xì)胞數(shù)量是最多的。移栽到低緯度地區(qū)后,松嶺小枝成熟細(xì)胞大小高于其它地點(diǎn)的移栽木,但在成熟細(xì)胞數(shù)量上的優(yōu)勢不明顯。孫吳小枝成熟和擴(kuò)大細(xì)胞都略低于來自松嶺的移栽木,可能是由于移栽后帽兒山水分高于原地所致。在生長季中,小枝擴(kuò)大細(xì)胞大小和數(shù)量差異都不顯著,說明小枝的徑向生長與溫度的升高不成正比。但在生長季初期的成熟細(xì)胞數(shù)量可以看出移栽木形成層活動(dòng)的起始時(shí)間還是略有不同。相對樹木的莖干而言,移栽后興安落葉松的小枝木質(zhì)部解剖特征相比主莖變化更為顯著,小枝對氣候變化(升溫)響應(yīng)更為敏感,樹木小枝的解剖特征相比主干更適合被用于研究樹木生長與氣候關(guān)系的變化。升溫后小枝擴(kuò)大細(xì)胞數(shù)量明顯減少(表3),表明升溫后,興安落葉松的生長策略已由強(qiáng)調(diào)多細(xì)胞同時(shí)分化的數(shù)量速度取勝模式向強(qiáng)調(diào)分化持續(xù)時(shí)間的時(shí)間取勝模式轉(zhuǎn)變[42]。

    年輪中記載了大量的信息,是眾多環(huán)境因子在不同時(shí)間段內(nèi)綜合作用結(jié)果,很難確定某一特定環(huán)境因子在特定時(shí)期內(nèi)對樹木生長的單獨(dú)影響[49]。一般認(rèn)為,溫度是限制高緯度地區(qū)樹木生長的關(guān)鍵因子,也是影響寒帶或溫帶森林樹木形成層活動(dòng)或解剖特征重要因素[42]。雖然移栽原地落葉松樹輪寬度并未隨溫度梯度而呈現(xiàn)出線性的上升或下降的趨勢,但整體表現(xiàn)出高緯度年輪較窄,低緯度年輪較寬的特征(圖1),這也從側(cè)面印證了緯度梯度上溫度的限制作用。孫吳地區(qū)落葉松輪寬整體上高于帶嶺地區(qū);一方面,可能是由于該地區(qū)土層相對較深(A+B層厚為(20.7+79.3) cm),能為落葉松生長提供更多的養(yǎng)分;另一方面,異常充沛的降水也為落葉松生長提供必要的水分條件;如2004、2005年孫吳地區(qū)生長季初期(4—5月)的降水量分別為104.7、158 mm遠(yuǎn)高于近幾十年來均值(69 mm),充沛的降水提早落葉松的形成層活動(dòng)并加速細(xì)胞的分裂與分化,進(jìn)而形成寬輪;2009年充沛的生長季(5—8月,530 mm近似年平均540 mm)降水也有可能是導(dǎo)致該年生長量較高的原因。此外,種間或種內(nèi)關(guān)系等也可能是產(chǎn)生以上差異的部分原因。

    4 結(jié)論

    本研究通過將帶嶺、孫吳、松嶺和塔河4個(gè)不同緯度的興安落葉松幼苗南移至帽兒山地區(qū),來模擬氣候變暖對興安落葉松年輪寬度以及莖、枝木質(zhì)部解剖特征的影響。溫度升高在一定程度上有利于興安落葉松的徑向生長,但是溫度升高與生長量的增加不成正比;降水、種間或種內(nèi)關(guān)系等溫度以外干擾因素的影響可能是產(chǎn)生上述現(xiàn)象的原因。不同緯度興安落葉松南移(升溫)后,溫度對落葉松的限制作用相對減弱,水分對樹木生長的限制作用增強(qiáng)。南移(升溫)后,興安落葉松的生長策略已由強(qiáng)調(diào)多細(xì)胞同時(shí)分化的數(shù)量取勝模式向強(qiáng)調(diào)分化持續(xù)時(shí)間的時(shí)間取勝模式轉(zhuǎn)變。升溫使興安落葉松的生長季結(jié)束延遲,間接延長了生長季長度。興安落葉松小枝解剖特征相比主干解剖特征對升溫響應(yīng)更敏感,小枝相比主干而言更適合被用于氣候變化研究。

    [1] IPCC. Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK: Cambridge University Press, 2007.

    [2] Marotzke J, Forster P M. Forcing, feedback and internal variability in global temperature trends. Nature, 2015, 517(7536): 565- 570.

    [3] Hansen J, Ruedy R, Sato M, Reynolds R. Global surface air temperature in 1995: return to pre-Pinatubo level. Geophysical Research Letters, 1996, 23(13): 1665- 1668.

    [4] 《氣候變化國家評估報(bào)告》編寫委員會(huì). 氣候變化國家評估報(bào)告. 北京: 科學(xué)出版社, 2007.

    [5] Bonsal B R, Zhang X, Vincent L A, Hogg W D. Characteristics of daily and extreme temperatures over Canada. Journal of Climate, 2001, 14(9): 1959- 1976.

    [6] Vincent L A, Mekis E. Changes in daily and extreme temperature and precipitation indices for Canada over the twentieth century. Atmosphere-Ocean, 2006, 44(2): 177- 193.

    [7] Krankina O N, Dixon R K. Forest management options to conserve and sequester terrestrial carbon in the Russian Federation. World Resource Review, 1994, 6(1): 88- 101.

    [8] 周廣勝, 王玉輝, 張新時(shí). 中國植被及生態(tài)系統(tǒng)對全球變化反應(yīng)的研究與展望. 中國科學(xué)院院刊, 1999, 14(1): 28- 32.

    [9] 周以良, 董世林, 聶紹荃. 黑龍江樹木志. 哈爾濱: 黑龍江科學(xué)出版社, 1986.

    [10] 徐化成. 中國大興安嶺森林. 北京: 科學(xué)出版社, 1998.

    [11] 周玉榮, 于振良, 趙士洞. 我國主要森林生態(tài)系統(tǒng)碳貯量和碳平衡. 植物生態(tài)學(xué)報(bào), 2000, 24(5): 518- 522.

    [12] 張新時(shí). 研究全球變化的植被-氣候分類系統(tǒng). 第四紀(jì)研究, 1993, 13(2): 157- 169.

    [13] 程肖俠. 氣候變化背景下中國東北森林的演替動(dòng)態(tài)[D]. 北京: 中國科學(xué)院研究生院(大氣物理研究所), 2007.

    [14] 何維明, 董鳴. 升高氣溫對旱柳光合和生長的影響. 林業(yè)科學(xué), 2003, 39(1): 160- 164.

    [15] 劉春延, 谷建才, 李吉躍, 陳平, 陸貴巧, 田國恒. 塞罕壩華北落葉松生長與氣候因子的相關(guān)分析. 北京林業(yè)大學(xué)學(xué)報(bào), 2009, 31(4): 102- 105.

    [16] 熊偉, 王彥輝, 于澎濤, 劉海龍, 時(shí)忠杰, 管偉. 六盤山南坡華北落葉松(Larixprincipis-rupprechtii) 樹干直徑生長及其對氣象因子的響應(yīng). 生態(tài)學(xué)報(bào), 2007, 27(2): 432- 441.

    [17] 張先亮, 何興元, 陳振舉, 崔明星, 黎娜. 大興安嶺山地樟子松徑向生長對氣候變暖的響應(yīng)——以滿歸地區(qū)為例. 應(yīng)用生態(tài)學(xué)報(bào), 2011, 22(12): 3101- 3108.

    [18] 王雨茜, 高露雙, 趙秀海. 長白山山楊年輪年表及其與氣候變化的關(guān)系. 東北林業(yè)大學(xué)學(xué)報(bào), 2013, 41(1): 10- 13.

    [19] Olano J M, Eugenio M, García-Cervigón A I, Folch M, Rozas V. Quantitative tracheid anatomy reveals a complex environmental control of wood structure in continental Mediterranean climate. International Journal of Plant Sciences, 2012, 173(2): 137- 149.

    [20] Carrer M, von Arx G, Castagneri D, Petit G. Distilling allometric and environmental information from time series of conduit size: the standardization issue and its relationship to tree hydraulic architecture. Tree Physiology, 2015, 35(1): 27- 33.

    [21] DeSoto L, De la Cruz M, Fonti P. Intra-annual patterns of tracheid size in the Mediterranean treeJuniperusthuriferaas an indicator of seasonal water stress. Canadian Journal of Forest Research, 2011, 41(6): 1280- 1294.

    [22] 趙西平, 郭明輝. 落葉松管胞形態(tài)與氣候變化關(guān)系的實(shí)證分析. 東北林業(yè)大學(xué)學(xué)報(bào), 2009, 37(7): 45- 48.

    [23] Deslauriers A, Morin H. Intra-annual tracheid production in balsam fir stems and the effect of meteorological variables. Trees, 2005, 19(4): 402- 408.

    [24] 趙西平, 郭明輝. 人工林紅松管胞形態(tài)對氣候變化的響應(yīng). 東北林業(yè)大學(xué)學(xué)報(bào), 2009, 37(8): 1- 4.

    [25] 高洪娜, 高瑞馨. 氣象因子對樹木生長量影響研究綜述. 森林工程, 2014, 30(2): 6- 9.

    [26] M?kinen H, Njd P, Saranp?? P. Seasonal changes in stem radius and production of new tracheids in Norway spruce. Tree Physiology, 2003, 23(14): 959- 968.

    [27] Barnett J R, Miller H. The effect of applied heat on graft union formation in dormantPiceasitchensis(Bong.) Carr. Journal of Experimental Botany, 1994, 45(1): 135- 143.

    [29] 霍宏. 氣候暖化對興安落葉松光合和生長影響的研究[D]. 哈爾濱: 東北林業(yè)大學(xué), 2007.

    [30] 朱良軍, 楊婧雯, 朱辰, 王曉春. 林隙干擾和升溫對小興安嶺紅松和臭冷杉徑向生長的影響. 生態(tài)學(xué)雜志, 2015, 34(8): 2085- 2095.

    [31] 劉實(shí), 閆敏華, 隋波. 東北三省冬季氣溫變化的有關(guān)研究進(jìn)展. 氣候變化研究進(jìn)展, 2009, 5(6): 357- 361.

    [32] Rossi S, Morin H, Deslauriers A, Plourde P Y. Predicting xylem phenology in black spruce under climate warming. Global Change Biology, 2011, 17(1): 614- 625.

    [33] 封曉輝, 程瑞梅, 肖文發(fā), 王瑞麗, 王曉榮, 劉澤彬. 基于日均溫度的華山松徑向生長敏感溫度研究. 生態(tài)學(xué)報(bào), 2012, 32(5): 1450- 1457.

    [34] Fonti P, Von Arx G, García-González I, Eilmann B, Sass-Klaassen U, G?rtner H, Eckstein D. Studying global change through investigation of the plastic responses of xylem anatomy in tree rings. New Phytologist, 2010, 185(1): 42- 53.

    [35] 方精云. 地理要素對我國溫度分布影響的數(shù)量評價(jià). 生態(tài)學(xué)報(bào), 1992, 12(2): 97- 104.

    [36] Chave J, Coomes D, Jansen S, Lewis S L, Swenson N G, Zanne A E. Towards a worldwide wood economics spectrum. Ecology letters, 2009, 12(4): 351- 366.

    [37] Zheng J Y, Ge Q S, Hao Z X, Wang W C. Spring phenophases in recent decades over eastern china and its possible link to climate changes. Climatic Change, 2006, 77(3): 449- 462.

    [38] Chen X Q, Hu B, Yu R. Spatial and temporal variation of phenological growing season and climate change impacts in temperate eastern China. Global Change Biology, 2005, 11(7): 1118- 1130.

    [39] Hasenauer H, Nemani R R, Schadauer K, Running S W. Forest growth response to changing climate between 1961 and 1990 in Austria. Forest Ecology and Management, 1999, 122(3): 209- 219.

    [40] Vaganov E A, Hughes M K, Shashkin A V. Growth Dynamics of Conifer Tree Rings. Berlin Heidelberg: Springer, 2006.

    [41] Drake J E, Aspinwall M J, Pfautsch S, Rymer P D, Reich P B, Smith R A, Crous K Y, Tissue D T, Ghannoum O, Tjoelker M G. The capacity to cope with climate warming declines from temperate to tropical latitudes in two widely distributedEucalyptusspecies. Global Change Biology, 2015, 21(1): 459- 472.

    [42] Vaganov E A, Anchukaitis K J, Evans M N. How well understood are the processes that create dendroclimatic records? A mechanistic model of the climatic control on conifer tree-ring growth dynamics // Hughes M K, Swetnam T W, Diaz H F, eds. Dendroclimatology. Netherlands: Springer, 2011: 11: 37- 75.

    [43] 李牧. 東北三大硬闊年輪與氣候關(guān)系[D]. 哈爾濱: 東北林業(yè)大學(xué), 2013.

    [44] Mikola P. Temperature and tree growth near the northern timberline // Kozlowski T T, ed. Tree Growth. New York: Ronald Press, 1962: 265- 274.

    [45] Sun S C, Jin D M, Shi P L. The leaf size-twig size spectrum of temperate woody species along an altitudinal gradient: an invariant allometric scaling relationship. Annals of Botany, 2006, 97(1): 97- 107.

    [46] Westoby M, Wright I J. The leaf size-twig size spectrum and its relationship to other important spectra of variation among species. Oecologia, 2003, 135(4): 621- 628.

    [47] 楊冬梅. 貢嘎山植物小枝功能性狀在不同生活型和不同海拔間的比較研究[D]. 成都: 中國科學(xué)院成都生物研究所, 2009.

    [48] 郭慶學(xué), 柴捷, 錢鳳, 張騰達(dá), 劉文丹, 陶建平. 不同木本植物功能型當(dāng)年生小枝功能性狀差異. 生態(tài)學(xué)雜志, 2013, 32(6): 1465- 1470.

    Impact of simulated climate warming on the radial growth ofLarixgmeliniiin northeast China

    LIU Yujia,GU Zhuoxin,WANG Xiaochun*

    CenterforEcologicalResearch,NortheastForestryUniversity,Harbin150040,China

    Climate warming has been widely reported across the globe, especially in high-latitude regions. There is an increasing body of evidence that a temperature increase would greatly affect tree growth and forest ecosystem. To simulate the effects of climate warming on the radial growth ofLarixgmelinii, eight-year-old young trees in Dailing, Sunwu, Songling, and Tahe (from low to high latitude) were transplanted to Maoershan (the lowest latitude) in 2004. We collected the micro-growing cores from the main trunks and segments of the branchlets from trees growing at the four transplanting origins and the one transplanting destination in the growing season of 2013. The enlarged and mature cells in main trunk and branchlet in each year were identified, and their sizes and quantities were measured by paraffin section methods. The results showed that the size and quantity of enlarged and mature cells in theLarixgmeliniitrunk increased with the simulated temperature increase, but the radial growth increment of the main trunk did not respond to the temperature increase. In the early part of the growing season (June), the quantity of mature cells in Dailing was significantly higher than that of the other two groups, whereas the cell number in August was relatively lower, indicating a pattern of continuous growth. In the late part of the growing season (August), the size and differentiation number of enlarged cells in the transplanting origin of Songling was 5 μm and 4, respectively, whereas these values were 9 μm and 7 in the transplanting destination (Maoershan), representing a nearly 2-fold increase. The size of mature cells in the transplanting origin of Songling was larger than that in the other two groups, which may be related to the change in water status. Temperature increase didn′t directly change the physiological characteristics, such as plant growth period, which might be controlled by genetic factors. We found that the branchlets were more affected by the temperature increase than the main trunk. The radial growth ofLarixgmelinii, including the diameter, number, and growth period of the tracheids, is expected to decline when global temperature increases, especially during epochs with the contemporary precipitation decrease. Our provenance experiments provide important original data to benchmark the predictions on the climate warming effects onLarixgmeliniigrowth and forest ecosystem changes in northeast China.

    Larixgmelinii; tree rings; global warming; cambium; slicing technology

    國家“十二五”科技支撐項(xiàng)目(2011BAD37B01);教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃(NCET- 12-0810);長江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃(IRT- 15R09)共同資助

    2015- 12- 17; 網(wǎng)絡(luò)出版日期:2016- 10- 29

    10.5846/stxb201512172515

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

    劉玉佳,顧卓欣,王曉春.模擬氣候變暖對東北興安落葉松徑向生長的影響.生態(tài)學(xué)報(bào),2017,37(8):2684- 2693.

    Liu Y J,Gu Z X,Wang X C.Impact of simulated climate warming on the radial growth ofLarixgmeliniiin northeast China.Acta Ecologica Sinica,2017,37(8):2684- 2693.

    猜你喜歡
    松嶺小枝興安
    江門市老促會(huì)
    源流(2024年2期)2024-05-04 21:42:24
    白鷺和大白鷺羽毛顯微結(jié)構(gòu)觀察
    高松嶺:一副干凈的手銬
    北風(fēng)的圍巾
    這個(gè)條件好
    上海故事(2020年10期)2020-12-14 04:03:16
    興安行吟(組詩)
    北極光(2019年12期)2020-01-18 06:22:16
    興安四月樹
    文苑(2019年22期)2019-12-07 05:29:32
    吉林松嶺老禿頂子春季大雪 極致美景“驚艷了世界”
    環(huán)境(2019年5期)2019-05-24 02:32:19
    追夢興安
    草原歌聲(2018年2期)2018-12-03 08:14:12
    回望興安
    鴨綠江(2013年12期)2013-03-11 19:42:10
    久久久久久免费高清国产稀缺| 精华霜和精华液先用哪个| 国产精品免费一区二区三区在线| 一边摸一边抽搐一进一小说| 国产高清videossex| 一本大道久久a久久精品| 麻豆成人av在线观看| 亚洲美女黄片视频| 亚洲七黄色美女视频| 午夜激情福利司机影院| 窝窝影院91人妻| 精品一区二区三区视频在线观看免费| 88av欧美| 久久性视频一级片| 中文字幕最新亚洲高清| 亚洲中文字幕日韩| av欧美777| 18禁美女被吸乳视频| 国产日本99.免费观看| 制服人妻中文乱码| 国产午夜福利久久久久久| 日韩欧美一区二区三区在线观看| 99riav亚洲国产免费| 日本三级黄在线观看| 国产熟女xx| 亚洲av成人不卡在线观看播放网| 国产精品自产拍在线观看55亚洲| 久久精品国产99精品国产亚洲性色| 亚洲午夜理论影院| 少妇人妻一区二区三区视频| 大型av网站在线播放| 欧美3d第一页| 欧美日韩瑟瑟在线播放| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 国产v大片淫在线免费观看| netflix在线观看网站| 午夜福利高清视频| 国产精品精品国产色婷婷| 色播亚洲综合网| 精品欧美一区二区三区在线| 欧美不卡视频在线免费观看 | 人妻久久中文字幕网| 国产精品亚洲美女久久久| av中文乱码字幕在线| 午夜福利免费观看在线| 欧美3d第一页| 国产乱人伦免费视频| 巨乳人妻的诱惑在线观看| 免费看美女性在线毛片视频| 欧美极品一区二区三区四区| 久热爱精品视频在线9| avwww免费| 亚洲aⅴ乱码一区二区在线播放 | 亚洲中文av在线| 18禁美女被吸乳视频| 国产区一区二久久| 欧美日韩瑟瑟在线播放| 国产精品一及| 一级作爱视频免费观看| 亚洲国产日韩欧美精品在线观看 | 久久精品91蜜桃| 精品国产美女av久久久久小说| 亚洲国产欧美一区二区综合| 两个人看的免费小视频| 日韩免费av在线播放| 国产一区二区在线av高清观看| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 久久香蕉精品热| 在线免费观看的www视频| 午夜免费激情av| 最近最新中文字幕大全电影3| 国产亚洲av高清不卡| 少妇被粗大的猛进出69影院| 1024手机看黄色片| 人妻丰满熟妇av一区二区三区| 精品一区二区三区视频在线观看免费| 精品久久久久久,| 黄色视频不卡| 91国产中文字幕| 日韩欧美国产在线观看| 亚洲色图av天堂| 神马国产精品三级电影在线观看 | 久久人妻av系列| 日韩欧美在线乱码| 手机成人av网站| 国产av在哪里看| 黄片小视频在线播放| 久久中文字幕一级| 国产精华一区二区三区| av片东京热男人的天堂| 黑人欧美特级aaaaaa片| www日本在线高清视频| 12—13女人毛片做爰片一| 国产一区二区在线观看日韩 | 久久人妻av系列| 国产精品国产高清国产av| 大型黄色视频在线免费观看| avwww免费| 午夜福利免费观看在线| 一级片免费观看大全| 97超级碰碰碰精品色视频在线观看| 正在播放国产对白刺激| 国产免费男女视频| 午夜视频精品福利| 九九热线精品视视频播放| 99久久精品国产亚洲精品| 久久精品亚洲精品国产色婷小说| 久久婷婷人人爽人人干人人爱| 亚洲精品久久国产高清桃花| 99热这里只有是精品50| 午夜福利在线在线| 亚洲精华国产精华精| 久久国产精品人妻蜜桃| 亚洲真实伦在线观看| 欧美成人性av电影在线观看| 久久人妻av系列| 亚洲中文字幕一区二区三区有码在线看 | 在线a可以看的网站| 亚洲精品色激情综合| 成人av在线播放网站| 亚洲成av人片在线播放无| 欧美日本亚洲视频在线播放| 18禁裸乳无遮挡免费网站照片| 最新美女视频免费是黄的| 12—13女人毛片做爰片一| 天天添夜夜摸| 正在播放国产对白刺激| 日本五十路高清| 欧美成人性av电影在线观看| 精品国内亚洲2022精品成人| 国产成人系列免费观看| 真人一进一出gif抽搐免费| 真人做人爱边吃奶动态| 亚洲国产看品久久| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美日韩东京热| 国产蜜桃级精品一区二区三区| 久久久久九九精品影院| 免费一级毛片在线播放高清视频| 国产亚洲精品第一综合不卡| 亚洲五月婷婷丁香| 亚洲专区字幕在线| 亚洲精品色激情综合| 久久这里只有精品19| 天天一区二区日本电影三级| 亚洲精品美女久久久久99蜜臀| 中文字幕精品亚洲无线码一区| 日韩国内少妇激情av| 国产蜜桃级精品一区二区三区| 亚洲色图 男人天堂 中文字幕| 色综合婷婷激情| 亚洲国产看品久久| 99久久国产精品久久久| 欧美 亚洲 国产 日韩一| 国产免费av片在线观看野外av| 久久久国产成人免费| 一本综合久久免费| 88av欧美| 亚洲avbb在线观看| 老熟妇仑乱视频hdxx| 久久人妻福利社区极品人妻图片| 国产在线精品亚洲第一网站| 亚洲成人中文字幕在线播放| 在线观看日韩欧美| 亚洲国产精品合色在线| 欧美日韩中文字幕国产精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲 欧美 日韩 在线 免费| 少妇裸体淫交视频免费看高清 | 亚洲国产精品sss在线观看| 久久天躁狠狠躁夜夜2o2o| 国产探花在线观看一区二区| av欧美777| 欧美人与性动交α欧美精品济南到| 在线观看舔阴道视频| av欧美777| 午夜精品久久久久久毛片777| 女人爽到高潮嗷嗷叫在线视频| 国产高清激情床上av| 制服诱惑二区| 亚洲自偷自拍图片 自拍| 日韩欧美一区二区三区在线观看| 桃色一区二区三区在线观看| 亚洲精品在线美女| 亚洲七黄色美女视频| 99热这里只有是精品50| 在线观看一区二区三区| 午夜激情福利司机影院| 亚洲 欧美 日韩 在线 免费| 黄色成人免费大全| 麻豆国产97在线/欧美 | 久久国产乱子伦精品免费另类| 黄频高清免费视频| 国产单亲对白刺激| 老司机深夜福利视频在线观看| 99国产极品粉嫩在线观看| 中文字幕人妻丝袜一区二区| 日本 欧美在线| 老司机福利观看| 一本大道久久a久久精品| 亚洲精品一区av在线观看| 亚洲av电影在线进入| 操出白浆在线播放| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久电影 | 超碰成人久久| 亚洲欧美日韩无卡精品| 舔av片在线| 精品无人区乱码1区二区| 99久久99久久久精品蜜桃| 神马国产精品三级电影在线观看 | 欧美日本视频| 97碰自拍视频| av有码第一页| 亚洲欧美激情综合另类| 亚洲欧美日韩东京热| 久久亚洲精品不卡| 日本 欧美在线| 日韩欧美国产在线观看| 亚洲精品中文字幕在线视频| 欧美性猛交黑人性爽| 桃色一区二区三区在线观看| 在线观看舔阴道视频| 日韩 欧美 亚洲 中文字幕| 日本 欧美在线| 99热这里只有是精品50| 午夜福利高清视频| 国产黄a三级三级三级人| 亚洲欧美日韩高清专用| 国产精品一区二区免费欧美| 久久久久久免费高清国产稀缺| 精品乱码久久久久久99久播| 真人一进一出gif抽搐免费| 高潮久久久久久久久久久不卡| 黄色成人免费大全| 久久中文看片网| 成年人黄色毛片网站| 欧美国产日韩亚洲一区| 日本成人三级电影网站| 色综合亚洲欧美另类图片| 亚洲五月天丁香| 久久人妻av系列| 国产精品综合久久久久久久免费| 一个人观看的视频www高清免费观看 | 特大巨黑吊av在线直播| 婷婷精品国产亚洲av在线| a在线观看视频网站| 成人永久免费在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 国产成人av激情在线播放| 巨乳人妻的诱惑在线观看| 一进一出抽搐动态| 久久婷婷成人综合色麻豆| 国产成人精品无人区| 日日摸夜夜添夜夜添小说| 成人国产一区最新在线观看| 黄频高清免费视频| 久久午夜综合久久蜜桃| 午夜福利欧美成人| 欧美3d第一页| 婷婷精品国产亚洲av在线| 91国产中文字幕| 婷婷亚洲欧美| 听说在线观看完整版免费高清| 99热这里只有是精品50| 国产亚洲精品一区二区www| 国产精品久久久人人做人人爽| 国产日本99.免费观看| 色哟哟哟哟哟哟| 99久久无色码亚洲精品果冻| 免费在线观看视频国产中文字幕亚洲| 中文字幕熟女人妻在线| 窝窝影院91人妻| 亚洲激情在线av| 久久草成人影院| 男人舔女人的私密视频| 一级作爱视频免费观看| 身体一侧抽搐| 禁无遮挡网站| 夜夜看夜夜爽夜夜摸| 岛国在线免费视频观看| 亚洲自偷自拍图片 自拍| av有码第一页| 一级毛片精品| АⅤ资源中文在线天堂| 久久久久久久久免费视频了| 美女黄网站色视频| 美女扒开内裤让男人捅视频| 18禁黄网站禁片免费观看直播| 99久久国产精品久久久| 18禁观看日本| 国产真实乱freesex| 亚洲熟女毛片儿| 大型av网站在线播放| 国产亚洲精品一区二区www| 国产一区二区在线av高清观看| 制服人妻中文乱码| 国产91精品成人一区二区三区| 国产v大片淫在线免费观看| 国产精品久久久久久亚洲av鲁大| 三级毛片av免费| 精品国产亚洲在线| 青草久久国产| 亚洲国产高清在线一区二区三| 亚洲国产欧美一区二区综合| 亚洲自偷自拍图片 自拍| 90打野战视频偷拍视频| 99国产精品99久久久久| 亚洲精品粉嫩美女一区| 国产99久久九九免费精品| 亚洲成人国产一区在线观看| 伦理电影免费视频| 一进一出好大好爽视频| 天堂影院成人在线观看| 在线永久观看黄色视频| 亚洲第一电影网av| 色综合亚洲欧美另类图片| 亚洲av电影不卡..在线观看| av国产免费在线观看| 国内精品久久久久久久电影| 又粗又爽又猛毛片免费看| 国产熟女午夜一区二区三区| 久久婷婷成人综合色麻豆| 俺也久久电影网| 国语自产精品视频在线第100页| 性色av乱码一区二区三区2| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 久久国产精品人妻蜜桃| 88av欧美| 亚洲欧美日韩高清在线视频| 亚洲熟妇熟女久久| 成人欧美大片| 国产亚洲欧美在线一区二区| 亚洲中文日韩欧美视频| or卡值多少钱| 国产亚洲精品av在线| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 国产精品九九99| 啦啦啦观看免费观看视频高清| 久久久久久久久久黄片| 精品国产乱子伦一区二区三区| 国产乱人伦免费视频| 制服丝袜大香蕉在线| 韩国av一区二区三区四区| 久久久久久大精品| 亚洲成av人片在线播放无| 亚洲无线在线观看| 韩国av一区二区三区四区| 90打野战视频偷拍视频| 久久久久免费精品人妻一区二区| 在线观看美女被高潮喷水网站 | 亚洲成人国产一区在线观看| 亚洲第一欧美日韩一区二区三区| 丁香欧美五月| 亚洲av电影不卡..在线观看| 国产又色又爽无遮挡免费看| 特大巨黑吊av在线直播| 国产高清激情床上av| 国产激情久久老熟女| 美女午夜性视频免费| 久久精品国产99精品国产亚洲性色| 国产三级在线视频| 欧美中文日本在线观看视频| 亚洲国产看品久久| www.熟女人妻精品国产| 久久精品国产清高在天天线| 亚洲在线自拍视频| 国产精品一及| 欧美久久黑人一区二区| 一二三四社区在线视频社区8| 夜夜躁狠狠躁天天躁| 亚洲av美国av| 一级毛片精品| 国产精品自产拍在线观看55亚洲| 老司机深夜福利视频在线观看| 男女之事视频高清在线观看| 成年版毛片免费区| 亚洲精华国产精华精| 变态另类成人亚洲欧美熟女| 俺也久久电影网| 日本黄色视频三级网站网址| 99在线视频只有这里精品首页| 日韩欧美 国产精品| 香蕉国产在线看| 亚洲成人久久性| 麻豆成人午夜福利视频| 亚洲 国产 在线| 国产成人一区二区三区免费视频网站| 久久人人精品亚洲av| 欧美丝袜亚洲另类 | 亚洲无线在线观看| 欧美一级a爱片免费观看看 | 国产成人欧美在线观看| 12—13女人毛片做爰片一| 啪啪无遮挡十八禁网站| netflix在线观看网站| 19禁男女啪啪无遮挡网站| 欧美大码av| av在线播放免费不卡| 欧美黄色淫秽网站| 亚洲成人久久性| 日韩高清综合在线| 国产熟女xx| 男女床上黄色一级片免费看| 老司机福利观看| 天天躁夜夜躁狠狠躁躁| 一边摸一边抽搐一进一小说| netflix在线观看网站| 久久99热这里只有精品18| 日本撒尿小便嘘嘘汇集6| 国产午夜福利久久久久久| 99久久无色码亚洲精品果冻| 99国产精品99久久久久| 亚洲成人国产一区在线观看| 狂野欧美激情性xxxx| 欧美一区二区精品小视频在线| 俺也久久电影网| 老熟妇乱子伦视频在线观看| 亚洲av电影在线进入| 夜夜躁狠狠躁天天躁| 99久久精品热视频| 色综合欧美亚洲国产小说| 国产男靠女视频免费网站| 日本一区二区免费在线视频| 妹子高潮喷水视频| 国模一区二区三区四区视频 | 精品不卡国产一区二区三区| 精品国产乱码久久久久久男人| 男男h啪啪无遮挡| 变态另类丝袜制服| 欧美一级a爱片免费观看看 | 欧美大码av| 国产又色又爽无遮挡免费看| 免费在线观看成人毛片| ponron亚洲| 国产av又大| 久久国产精品影院| 国产三级黄色录像| 女人被狂操c到高潮| 国产三级中文精品| 成人手机av| 777久久人妻少妇嫩草av网站| 黄频高清免费视频| 男女视频在线观看网站免费 | 欧美日韩中文字幕国产精品一区二区三区| 欧美大码av| 欧美午夜高清在线| 中出人妻视频一区二区| 午夜视频精品福利| 久久久久国内视频| 超碰成人久久| or卡值多少钱| 成人欧美大片| 亚洲,欧美精品.| 国产激情久久老熟女| 久久国产乱子伦精品免费另类| av福利片在线| 亚洲中文字幕日韩| 少妇裸体淫交视频免费看高清 | 午夜福利在线在线| 后天国语完整版免费观看| 精品少妇一区二区三区视频日本电影| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 久久久久久大精品| 欧美激情久久久久久爽电影| 91国产中文字幕| 此物有八面人人有两片| 精品久久久久久成人av| 国产欧美日韩精品亚洲av| 韩国av一区二区三区四区| 99久久无色码亚洲精品果冻| 两个人的视频大全免费| 免费在线观看成人毛片| 男女做爰动态图高潮gif福利片| 欧美一级a爱片免费观看看 | 国产精品亚洲一级av第二区| 亚洲人成电影免费在线| 日韩 欧美 亚洲 中文字幕| 欧美久久黑人一区二区| 亚洲性夜色夜夜综合| 欧美精品亚洲一区二区| 久99久视频精品免费| 97人妻精品一区二区三区麻豆| 精品久久久久久久久久久久久| 九色成人免费人妻av| 老司机深夜福利视频在线观看| 这个男人来自地球电影免费观看| www.自偷自拍.com| 欧美另类亚洲清纯唯美| 哪里可以看免费的av片| 少妇熟女aⅴ在线视频| 禁无遮挡网站| 色在线成人网| 国产又黄又爽又无遮挡在线| 露出奶头的视频| 欧美日韩精品网址| 一级作爱视频免费观看| 中文字幕人成人乱码亚洲影| 啦啦啦观看免费观看视频高清| 国内毛片毛片毛片毛片毛片| 中文亚洲av片在线观看爽| 精品无人区乱码1区二区| 亚洲精品久久成人aⅴ小说| 国内精品久久久久精免费| 精品国产乱码久久久久久男人| 国产视频一区二区在线看| x7x7x7水蜜桃| 婷婷六月久久综合丁香| 天堂av国产一区二区熟女人妻 | 在线观看舔阴道视频| 中国美女看黄片| 国产av又大| 亚洲真实伦在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产精品日韩av在线免费观看| 日韩精品中文字幕看吧| 别揉我奶头~嗯~啊~动态视频| 搡老岳熟女国产| 久久久久久久久免费视频了| 免费观看人在逋| 亚洲成人免费电影在线观看| 成年免费大片在线观看| 精品乱码久久久久久99久播| 成人永久免费在线观看视频| 禁无遮挡网站| 美女大奶头视频| 黄色a级毛片大全视频| 黄频高清免费视频| 老汉色∧v一级毛片| 亚洲第一欧美日韩一区二区三区| 免费高清视频大片| 亚洲最大成人中文| 别揉我奶头~嗯~啊~动态视频| 99re在线观看精品视频| 一级毛片精品| 亚洲国产欧美一区二区综合| 亚洲av成人不卡在线观看播放网| 精品久久久久久久久久久久久| 嫩草影院精品99| 亚洲专区国产一区二区| 狂野欧美白嫩少妇大欣赏| 久久亚洲精品不卡| 在线国产一区二区在线| 丝袜人妻中文字幕| 亚洲美女黄片视频| 夜夜躁狠狠躁天天躁| 国产高清有码在线观看视频 | 日韩大码丰满熟妇| 村上凉子中文字幕在线| 99热这里只有是精品50| 日日干狠狠操夜夜爽| 两个人的视频大全免费| 精品久久久久久成人av| 午夜福利免费观看在线| 熟女电影av网| 久久精品国产99精品国产亚洲性色| 五月玫瑰六月丁香| 色老头精品视频在线观看| 日本一区二区免费在线视频| 色综合亚洲欧美另类图片| 搡老熟女国产l中国老女人| 黄片大片在线免费观看| 国产精品久久久av美女十八| 国产激情欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 国产av在哪里看| 久久久久久久久免费视频了| 岛国在线免费视频观看| 成人国产一区最新在线观看| 午夜久久久久精精品| 久久中文字幕人妻熟女| 欧美日韩黄片免| 校园春色视频在线观看| 亚洲无线在线观看| 熟女电影av网| 精品欧美一区二区三区在线| 亚洲熟妇熟女久久| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美激情综合另类| 男女床上黄色一级片免费看| 婷婷亚洲欧美| 久久久久性生活片| 波多野结衣高清作品| 亚洲国产日韩欧美精品在线观看 | 97人妻精品一区二区三区麻豆| √禁漫天堂资源中文www| 美女黄网站色视频| 一级毛片高清免费大全| 18禁美女被吸乳视频| 在线a可以看的网站| 成年人黄色毛片网站| 成人国产一区最新在线观看| 成年女人毛片免费观看观看9| 亚洲美女视频黄频| 久久久精品大字幕| 极品教师在线免费播放| 午夜激情av网站| 国产av又大| 日韩成人在线观看一区二区三区| 国产片内射在线| 19禁男女啪啪无遮挡网站| 一本大道久久a久久精品| 色精品久久人妻99蜜桃| 悠悠久久av| 性欧美人与动物交配| 在线永久观看黄色视频| 成人欧美大片| 日日干狠狠操夜夜爽| 欧美日韩黄片免| 欧美3d第一页|