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

    川西亞高山不同年齡紫果云杉徑向生長對氣候因子的響應(yīng)

    2016-03-10 06:00:39趙志江康東偉李俊清
    生態(tài)學(xué)報 2016年1期
    關(guān)鍵詞:氣候因子年齡

    趙志江,康東偉,李俊清,*

    1 福建中咨工程咨詢有限公司,福州 350003

    2 青島農(nóng)業(yè)大學(xué)園林與林學(xué)院,青島 266109

    3 北京林業(yè)大學(xué)省部共建森林培育與保護教育部重點實驗室, 北京 100083

    川西亞高山不同年齡紫果云杉徑向生長對氣候因子的響應(yīng)

    趙志江1,2,康東偉3,李俊清3,*

    1 福建中咨工程咨詢有限公司,福州350003

    2 青島農(nóng)業(yè)大學(xué)園林與林學(xué)院,青島266109

    3 北京林業(yè)大學(xué)省部共建森林培育與保護教育部重點實驗室, 北京100083

    摘要:運用樹木年輪氣候?qū)W的基本方法,建立王朗自然保護區(qū)紫果云杉在集中分布上限區(qū)域的年輪寬度年表,選取差值年表分析不同年齡云杉的徑向生長同逐月氣候因子的相關(guān)及響應(yīng)關(guān)系,結(jié)果顯示:幼齡組云杉年表的敏感度高于中齡組和老齡組云杉,幼齡組云杉對生長季前及生長季的氣溫狀況顯著正相關(guān);中齡組云杉年表僅與當年4月份和7月份的月平均最低氣溫顯著正相關(guān);老齡組云杉的年輪寬度指數(shù)同上年生長季(上年8月份)的月平均氣溫和月平均最低溫顯著負相關(guān),上年生長季高溫的“滯后效應(yīng)”在老齡組云杉體現(xiàn)的更為突出;幼齡組與中齡組云杉對當年6月份降水持續(xù)增加顯示出明顯的負相關(guān)關(guān)系,上年12月份的降水會對幼齡組和老齡組云杉徑向生長不利。研究表明幼齡組云杉包含的氣候信息要優(yōu)于中齡組和老齡組云杉,在該區(qū)域進行相關(guān)研究時應(yīng)根據(jù)研究需要選取不同年齡跨度的云杉年表。

    關(guān)鍵詞:紫果云杉;年齡;徑向生長;氣候因子

    樹木在自然環(huán)境下的生長主要依賴于生長期的氣候條件,但是諸如微氣候、上年生長、樹齡等附加因子同樣會對樹木的生長造成一定的影響[1]。通常情況下,依據(jù)樹木年輪氣候?qū)W的技術(shù)操作流程,通過移除非氣候“噪音”來放大樹輪資料中的氣候信號,這個過程稱為“標準化”,其前提假設(shè)為:經(jīng)過“標準化”處理之后,樹木生長與氣候的關(guān)系不會依賴于樹木年齡而有所不同,即在給定的時間序列內(nèi)不同年齡的樹木對于氣候的響應(yīng)一致[1- 2]。但是,生理生態(tài)學(xué)的相關(guān)研究表明,與生長相關(guān)的環(huán)境信號有可能會依賴于年齡因素,例如,多數(shù)針葉樹種的老樹光合速率較低[3]。因此,年齡對于樹木生長-氣候響應(yīng)的不確定性影響逐漸受到樹木年輪氣候?qū)W研究的重視[4]。

    近些年,相關(guān)研究主要側(cè)重于檢驗樹木生長對于氣候的響應(yīng)關(guān)系在不同齡級或徑級間是否一致,研究結(jié)果會因樹種和地點的不同而有所差異。部分研究表明年齡不會明顯影響樹木年表的氣候響應(yīng),例如,加拿大南部山區(qū)的幼齡與老齡落葉松(LarixlyaliiParl.)的年表對于氣候響應(yīng)基本一致[5- 6],歐洲中部阿爾卑斯山的五針松(Pinuscembra)不同年齡組年表所含的氣候信號沒有明顯變化[7],渾善達克沙地不同徑級油松的年表對氣候因子的響應(yīng)存在明顯的一致性[8]。

    然而,多數(shù)研究表明幼齡樹與老齡樹的徑向生長與氣候因子間的關(guān)系存在明顯差異。例如,在中部斯堪的納維亞山脈,老齡的(大于250a)歐洲赤松(PinussylvestrisL.)比中齡的對氣候響應(yīng)更為敏感[9];祁連山低于200a的祁連圓柏(SabinaprzewalskiiKom.)與其他4組高齡樹相比,在氣候響應(yīng)方面存在明顯的不一致[10];葡萄牙西北部的海岸松(Pinuspinaster)幼齡樹的早材對于氣候響應(yīng)更為敏感,而老齡樹的晚材更易受到更易受到氣候的影響[11];我國長白山北坡低齡油松(Pinuskoraiensis)年表與高齡年表亦存在明顯差異[12]。

    川西亞高山位于我國青藏高原東緣,屬于典型的氣候敏感區(qū)域。紫果云杉(Piceapurpurea)作為該區(qū)域分布的典型優(yōu)勢樹種,年齡因素是否會對其年表以及氣候響應(yīng)產(chǎn)生影響亟待揭示,其結(jié)果會直接影響今后利用紫果云杉樹輪資料進行氣候重建的可靠性。為避免出現(xiàn)因取樣年齡不同造成掩蓋部分氣候信息的現(xiàn)象,本文旨在闡明不同齡級的紫果云杉年表對于氣候因子的響應(yīng)規(guī)律,以期為在本區(qū)域深入開展樹輪生態(tài)學(xué)的針對性研究提供合理、科學(xué)的依據(jù)。

    1材料與方法

    1.1野外取樣

    本文研究地點位于四川省王朗國家級自然保護區(qū),保護區(qū)面積為322.97km2,海拔跨度2300—4980m,屬于丹巴-松潘半濕潤季風(fēng)氣候,冬季寒冷干燥、日照強烈、降水較少,夏季溫暖濕潤、日照較少、降水集中。年均溫為2.3℃,年降水量約為1100mm[13]。區(qū)內(nèi)紫果云杉主要分布于陽坡海拔2600—3000m的范圍內(nèi)。在紫果云杉典型集中分布區(qū)的上限布設(shè)一條20m(坡向方向)×100m(水平方向)的代表性樣帶,樣帶中心海拔為3003m(32°54.2′ N,104°03.2′ E),用生長錐對樣帶內(nèi)胸徑大于10cm云杉在胸徑處(離地高度為130cm)進行取樣,胸徑小于20cm的取1個樣芯,胸徑大于20cm的取1—2個樣芯。由于該區(qū)域大徑級(胸徑大于40cm)紫果云杉分布較為分散,因此在樣帶外選取一定數(shù)量的大徑級樹進行補充取樣,各采樣點與樣帶處于同一海拔范圍,以滿足不同年齡樣芯的需求。

    1.2年齡分組及年表建立

    將樣芯處理完成并進行測量,根據(jù)Duncan[14]與Norton等[15]的方法判斷接近和未到髓心的樣芯的年齡,即樹木在胸徑處的年齡,然后按照100a為標準[2]將樣芯分為3組:幼齡組(Young Age Classes,YAC,小于100a的樣芯)、中齡組(Middle Age Classes,MAC,年齡介于100—200a的樣芯)、老齡組(Old Age Classes,OAC,大于200a的樣芯)[2]。分組完成后,利用COFECHA程序?qū)Ω髂挲g組的云杉年輪序列進行質(zhì)量檢驗并完成交叉定年[16],其中不包括心腐、破損、鉆偏嚴重的樣本,3個齡級交叉定年成功的樣芯數(shù)量分別為:幼齡組取自33棵樹的33個樣芯,中齡組取自24棵樹的27個樣芯,老齡組取自18棵樹的27個樣芯。然后利用ARSTAN軟件移除生長趨勢來建立年表[17],所采用的去趨勢方法為負指數(shù)函數(shù)擬合或任意斜率的線性回歸擬合[2],所建立的3個齡級云杉的標準年表見圖1。

    圖1 王朗自然保護區(qū)不同年齡紫果云杉標準年表與樣本數(shù)量Fig.1 Standard chronologies and sample numbers of P. purpurea for different age classes in the Wanglang Natural Reserve

    1.3敏感度與胸徑大小的關(guān)系

    姜慶彪等發(fā)現(xiàn)輪寬樣芯的敏感度同胸徑大小之間存在一定的關(guān)系[8],這在一定程度上也會反映樹木大小對于樹木敏感性的影響。敏感度的計算公式如下:

    (1)

    式中,xi是第i個年輪寬度值或指數(shù);xi+1是第i+1個年輪寬度值或指數(shù);n是該樣本年輪總數(shù)。

    本文利用SimaPlot 10.0軟件提供的線性回歸方程計算了本研究中各交叉定年成功且到達髓心的樣芯的敏感度同胸徑處年齡的關(guān)系。

    1.4不同年表同氣候因子的關(guān)系分析

    分析不同齡級云杉的差值年表同氣候因子的關(guān)系,包括相關(guān)分析與響應(yīng)函數(shù)分析,分析過程通過DendroCLIM 2002軟件完成[18],氣象數(shù)據(jù)通過山地小氣候模型(MTCLIM)對采樣點進行模擬,其中初始氣象數(shù)據(jù)選取與采樣點處于同一氣候區(qū)且直線距離最近的松潘氣象站(32°39′ N,103°34′ E,海拔2850m)的數(shù)據(jù),時間跨度為1951—2009年,選取上年6月至當年9月共計16個月份的月平均氣溫(Tmean)、月平均最高氣溫(Tmax)、月平均最低氣溫(Tmin)與月降水(P)。

    2結(jié)果與分析2.1不同齡級年表的質(zhì)量與統(tǒng)計分析

    不同齡級紫果云杉標準年表與差值年表的統(tǒng)計參數(shù)見表1。公共區(qū)間分析中的各參數(shù)均反映出差值年表的質(zhì)量要優(yōu)于標準年表。標準年表與差值年表的平均敏感度(MS)、標準差(SD)、樣芯間平均相關(guān)系數(shù)(RBAR)、樣本總體代表性(EPS)、信噪比(SNR)、第一特征根解釋量(PC1)均體現(xiàn)如下趨勢:即幼齡組與老齡組年表高于中齡組年表,尤以幼齡組年表最高。樣芯敏感度與對應(yīng)樣樹年齡的回歸分析表明(圖2),敏感度與胸徑大小之間不存在線性相關(guān)關(guān)系;從散點分布上看,除部分年齡較小的樣芯敏感度較大外,樣芯的敏感度與年齡之間并不無明顯規(guī)律。

    表1王朗自然保護區(qū)不同齡級紫果云杉標準年表與差值年表的特征統(tǒng)計

    Table 1Dendrochronological statistics for standard and residual chronologies of theP.purpureaat different age classes in the the Wanglang Nature Reserve

    統(tǒng)計特征Statisticscharacters幼齡組YACSTDRES中齡組MACSTDRES老齡組OACSTDRES樣芯數(shù)量/樹No.ofcores/trees33/3327/2427/18時間序列Recordperiod1911—20091852—20091667—2009平均敏感度MS0.1230.1480.1310.1230.1060.132標準差SD0.1720.1210.2570.1070.1940.114一階自相關(guān)系數(shù)AC10.591-0.1920.879-0.0020.777-0.093公共區(qū)間分析Commonintervalanaly-sis1950—20091900—20091800—2009樣芯間平均相關(guān)系數(shù)RBAR0.2810.3780.2630.2990.2630.352樣本總體代表性EPS0.9100.9400.8650.8850.8870.923信噪比SNR10.15515.7896.4137.6937.83111.970第一特征根解釋量PC1/%33.441.235.234.232.438.7

    YAC:幼齡組young age classes;MAC:中齡組middle age classes;OAC:老齡組old age classes;STD:標準年表standard chronologies;RES:差值年表residual chronologies;MS:平均敏感度 mean sensitivity;SD:標準差standard deviation;AC1:一階自相關(guān)系數(shù) autocorrelation order 1;RBAR:樣芯間平均相關(guān)系數(shù) mean correlation among radii;EPS:樣本總體代表性 express population signal;SNR:信噪比 signal-to-noise ratio;PC1:第一特征根解釋量 variance in first eigenvector

    圖2 王朗自然保護區(qū)紫果云杉樣芯敏感度與年齡大小的相關(guān)關(guān)系 Fig.2 Correlation between sensitivity of P. purpurea samplings and ages in the Wanglang Nature Reseve

    2.2不同齡級年表與氣候因子的關(guān)系

    由于各齡級云杉差值年表的信噪比等參數(shù)均優(yōu)于標準年表,因此本文選取差值年表同各氣候因子進行相關(guān)及響應(yīng)分析[19],結(jié)果見圖3。對于幼齡組云杉年表,相關(guān)分析的結(jié)果表明,年輪指數(shù)同下列氣候因子顯著正相關(guān):上年11月份、當年2、6、9月份的月平均氣溫,上年11月份、當年2、9月份的月平均最高氣溫,當年2、3、4、6、7、8、9月份的月平均最低氣溫;年輪指數(shù)還同上年12月與當年6月的降水顯著負相關(guān)。響應(yīng)函數(shù)分析的結(jié)果表明,年輪指數(shù)與當年2月份的月平均最高氣溫、月平均最低氣溫以及當年4月份的月平均最低氣溫顯著正相關(guān)。

    對于中齡組云杉,同年輪指數(shù)具有明顯相關(guān)關(guān)系的氣候因子相對較少,僅在相關(guān)分析中有所體現(xiàn)。中齡云杉的徑向生長同當年2月份、7月份的月平均最低氣溫顯著正相關(guān),同當年6月份的降水顯著負相關(guān)。

    對于老齡組云杉,相關(guān)分析顯示,上年8月份的月平均氣溫、月平均最低氣溫以及當年9月份的月平均最低氣溫與年輪指數(shù)顯著負相關(guān);相關(guān)與響應(yīng)函數(shù)分析的結(jié)果均表明上年12月份的降水對徑向生長起明顯的負作用。

    圖3 不同齡組紫果云杉差值年表與月氣候因子的相關(guān)及響應(yīng)關(guān)系(*P<0.05)Fig.3 Correlations and responses between P. purpurea residual chronology with different age classes and monthly climatic factors in the Wanglang Nature Reserve (* P<0.05)

    3討論

    3.1樹齡對年表參數(shù)的影響

    通過3個齡級年表特征對比(表1),發(fā)現(xiàn)幼齡云杉年表的質(zhì)量明顯高于中齡與老齡云杉年表,中齡云杉的年表質(zhì)量最差,說明幼齡云杉對于周邊環(huán)境的變化更為敏感,樣芯敏感性與年齡大小關(guān)系也印證了這一點(圖2)。姜慶彪等在對渾善達克沙地不同徑級油松(Pinustubulaeformis)研究時也發(fā)現(xiàn),小徑級油松具有更高的敏感度[8]。王婷分析了天山云杉(Piceaschrenkiana)隨年齡變化的年表特征,結(jié)果顯示中齡組的年表質(zhì)量要優(yōu)于幼齡組與老齡組[20]。一般認為大徑級[21]或大齡級[22]樹木對于氣候響應(yīng)更加的敏感,與本研究的結(jié)果有所差別,主要原因可能是不同地點不同樹種隨年齡(徑級)變化的生理過程不一致[8]。在王朗自然保護區(qū)紫果云杉的主要分布區(qū),幼齡或老齡云杉更適合進行氣候響應(yīng)分析研究。

    3.2不同樹齡年表對氣候因子的響應(yīng)差異

    不同樹齡年表對于氣溫的敏感性存在很大的差異性。幼齡組云杉對生長季前(上年11月份、當年2月份)及生長季(6月份、9月份)的氣溫狀況均顯示出明顯的相關(guān)性,生長季前及生長季升溫均有利于其徑向生長,特別是當年2、3、4月份及6、7、8、9月份的月平均最低氣溫。生長季前的升溫有利于幼齡云杉度過寒冷冬季,減少凍害的發(fā)生,同時通過增強細胞活性、促進營養(yǎng)物質(zhì)運輸?shù)确绞綖闃淠旧L做好充足準備[8]。生長季初(當年6月份)的升溫能夠提高云杉體內(nèi)與生長相關(guān)的酶活性、提高光合速率、促進營養(yǎng)物質(zhì)積累。生長季末(當年9月份)較高的氣溫有利于延長生長期時間,促進寬輪的產(chǎn)生。在這一區(qū)域最低氣溫有可能是影響幼齡云杉形成層活性的關(guān)鍵因素,由此限制木質(zhì)部的管胞分裂與擴展[23]。生長季較低的氣溫與土壤溫度會通過影響水分有效性限制形成層活性,進而不利于樹木徑向生長[24]。而中齡組云杉差值年表僅與當年4月份和7月份的月平均最低氣溫顯著正相關(guān),中齡組云杉標準年表同氣候的關(guān)系還顯示當年6月份較高的月平均氣溫和月平均最低氣溫有利于云杉徑向生長,說明生長季前及生長季低溫是限制中齡云杉徑向生長的主要氣溫因子,即中齡組云杉生長對當年生長季前及生長季最低氣溫較為敏感。與幼齡組和中齡組云杉不同,老齡組云杉的差值年表年輪寬度指數(shù)同上年生長季(上年8月份)的月平均氣溫和月平均最低溫顯示出明顯的負相關(guān)關(guān)系,另,標準年表的氣候相關(guān)分析也顯示上年7月份、8月份和當年9月份的月平均氣溫和月平均最高氣溫與輪寬指數(shù)顯著負相關(guān)??梢娚L旺季和生長季末的高溫會過度消耗云杉體內(nèi)儲存的光合產(chǎn)物,不利于云杉徑向生長。

    不同齡級云杉生長對于降水的響應(yīng)也不完全一致。幼齡組與中齡組云杉對當年6月份降水持續(xù)增加顯示出明顯的負相關(guān)關(guān)系,可能的原因是生長季充沛降雨造成云杉根部的無氧呼吸,從而消耗較多的有機物限制年輪寬度增加[25]。而幼齡組與老齡組云杉年表同冬季(上年12月份)的降水顯著負相關(guān),12月份過多的降水會加劇云杉的凍害或機械損傷,從而影響到下年的生長[26]。

    目前,還不能從機理上深入揭示不同齡級云杉生長對于氣候因子的響應(yīng)差異,部分學(xué)者認為遺傳基因[2]和林分動態(tài)[27]等造成的樹木生理機能的不同是其主要原因。樹齡不同,則生理機能不同,隨之樹木對氣候的響應(yīng)也會不同。例如,幼齡云杉在生長季末(當年9月份)會受到低溫脅迫,過早的遭受凍害,從而降低活性影響生長;而對于老齡云杉,由于長期的適應(yīng)過程,會對早期凍害產(chǎn)生一定的抗性不易受損,但是生長季末(當年9月份)的高溫會增強呼吸速率,不利于有機物質(zhì)的積累。幼齡與老齡云杉對于冬季(上年12月份)的降水更加敏感,而對于位于生長壯年的中齡云杉具有更強的抵御冬季嚴寒的能力。幼齡和中齡云杉的生長會受到當年生長季(當年6月份)充足降水的抑制,而對于老齡云杉而言,由于枝干粗壯、冠幅較大對于水分的需求量大,會降低根部無氧呼吸的強度。

    4結(jié)論

    在王朗自然保護區(qū)紫果云杉的集中分布上限區(qū)域,幼齡組云杉年表的統(tǒng)計參數(shù)要優(yōu)于中齡組與老齡組云杉年表,與之相對應(yīng),幼齡組云杉對于環(huán)境的敏感性也高于其他兩個齡組。在溫度方面,幼齡組云杉生長同當年的氣溫狀況表現(xiàn)出良好的相關(guān)性,老齡組云杉生長主要受到上年生長季高溫的顯著抑制,中齡組云杉生長主要得益于當年4月份和7月份較高的低溫。在降水方面,當年6月份持續(xù)降水將會明顯不利于幼齡組與中齡組云杉的生長,上年12月份的降水亦會對幼齡組和老齡組云杉產(chǎn)生抑制作用。鑒于此,本文認為今后在進行樹輪氣候?qū)W研究時,應(yīng)根據(jù)研究需要恰當選取建立年表的樣芯的年齡跨度,以免損失或放大氣候信息。

    致謝:感謝王朗自然保護區(qū)的大力支持。

    參考文獻(References):

    [1]吳祥定. 樹木年輪與氣候變化. 北京: 氣象出版社, 1990: 44-145.

    [2]Szeicz J M, MacDonald G M. Age-dependent tree-ring growth responses of subarctic white spruce to climate. Canadian Journal of Forest Research, 1994, 24(1): 120-132.

    [3]Bond B J. Age-related changes in photosynthesis of woody plants. Trends in Plant Science, 2000, 5(8): 349-353.

    [4]彭劍峰, 劉玉振, 王婷. 神農(nóng)山白皮松不同齡組年輪-氣候關(guān)系及PDSI 重建. 生態(tài)學(xué)報, 2014, 34(13): 3509-3518.

    [5]Colenutt M E, Luckman B H. The dendrochronological characteristics of alpine larch. Canadian Journal of Forest Research, 1995, 25(5): 777-789.

    [6]Colenutt M E, Luckman B H. Dendrochronological investigation ofLarixlyalliiat Larch Valley, Alberta. Canadian Journal of Forest Research, 1991, 21(8): 1222-1233.

    [7] Esper J, Niederer R, Bebi P, Frank D. Climate signal age effects-Evidence from young and old trees in the Swiss Engadin. Forest Ecology and Management, 2008, 255(11): 3783-3789.

    [8]姜慶彪, 高露雙, 王曉明, 趙秀海. 渾善達克沙地油松樹輪寬度與氣候因子的關(guān)系. 應(yīng)用與環(huán)境生物學(xué)報, 2012, 18(3): 405-410.

    [9]Linderholm H W, Linderholm K. Age-dependent climate sensitivity ofPinussylvestrisL. in the central Scandinavian Mountains. Boreal Environment Research, 2004, 9(4): 307-317.

    [10]Yu G R, Liu Y B, Wang X C, Ma K P. Age-dependent tree-ring growth responses to climate in Qilian juniper (SabinaprzewalskiiKom.). Trees, 2008, 22(2): 197-204.

    [11]Vieira J, Campelo F, Nabais C. Age-dependent responses of tree-ring growth and intra-annual density fluctuations ofPinuspinasterto Mediterranean climate. Trees, 2009, 23(2): 257-265.

    [12]王曉明, 趙秀海, 高露雙, 姜慶彪. 長白山北坡不同年齡紅松年表及其對氣候的響應(yīng). 生態(tài)學(xué)報, 2011, 31(21): 6378-6387.

    [13] Taylor A H, Jang S W, Zhao L J, Liang C P, Miao C J, Huang J Y. Regeneration patterns and tree species coexistence in old-growthAbies-Piceaforests in southwestern China. Forest Ecology and Management, 2006, 223(1-3): 303-317.

    [14]Duncan R P. An evaluation of errors in tree age estimates based on increment cores in Kahikatea (Dacrycarpusdacrydioides). New Zealand Natural Sciences, 1989, 16: 31-37.

    [15]Norton D A, Palmer J G, Ogden J. Dendroecological studies in New Zealand 1. An evaluation of tree age estimates based on increment cores. New Zealand Journal of Botany,1987, 25(3): 373-383.

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

    [17]Cook E R. A Time Series Analysis Approach to Tree Ring Standardization (Dendrochronology, Forestry, Dendroclimatology, Autoregressive Process)[D]. Tucson: The University of Arizona, 1985.

    [18]Biondi F, Waikul K. DENDROCLIM2002: A C++ program for statistical calibration of climate signals in tree-ring chronologies. Computers & Geosciences, 2004, 30(3): 303-311.

    [19]王曉明, 趙秀海, 高露雙, 姜慶彪. 長白山北坡林線處岳樺年輪年表及其與氣候的關(guān)系. 應(yīng)用與環(huán)境生物學(xué)報, 2012, 18(1): 9-16.

    [20]王婷. 天山中部不同海拔高度天山云杉林的生態(tài)學(xué)研究[D]. 武漢: 武漢大學(xué), 2004.

    [21]Chhin S, Hogg E H, Lieffers V J, Huang R M. Potential effects of climate change on the growth of lodge pole pine across diameter size classes and ecological regions. Forest Ecology and Management, 2008, 256(10): 1692-1703.

    [22]Carrer M, Urbinati C. Age-dependent tree-ring growth responses to climate inLarixdeciduaandPinuscembra. Ecology, 2004, 85(3): 730-740.

    [23]Deslauriers A, Morin H, Begin Y. Cellular phenology of annual ring formation ofAbiesbalsameain the Quebec boreal forest (Canada). Canadian Journal of Forest Research, 2003, 33(2): 190-200.

    [24]Liang E Y, Wang Y F, Xu Y, Liu B, Shao X M. Growth variation inAbiesgeorgeivar.smithiialong altitudinal gradients in the Sygera Mountains, southeastern Tibetan Plateau. Trees, 2010, 24(2): 363-373.

    [25]Lloyd A H, Graumlich L J. Holocene dynamics of treeline forests in the Sierra Nevada. Ecology, 1997, 78(4): 1199-1210.

    [26]徐倩倩. 長白山東坡長白落葉松生長和氣候因子關(guān)系的研究[D]. 北京: 北京林業(yè)大學(xué), 2011.

    [27] Szeicz J M, Macdonald G M. Recent white spruce dynamics at the subarctic alpine treeline of north-western Canada. Journal of Ecology, 1995, 83(5): 873-885.

    Age-dependent radial growth responses ofPiceapurpureato climatic factors in the subalpine region of Western Sichuan Province, China

    ZHAO Zhijiang1,2, KANG Dongwei3, LI Junqing3,*

    1FujianCIECCEngineeringConsultingCo.Ltd,Fuzhou350003,China

    2CollegeofLandscapeArchitectureandForestry,QingdaoAgriculturalUniversity,Qingdao266109,China

    3TheKeyLaboratoryforSilvicultureandConservationofMinistryofEducation,BeijingForestryUniversity,Beijing100083,China

    Abstract:Although tree growth is mainly dependent on natural climate conditions, factors such as microclimate, growth of the previous year, tree age, and others may also affect tree growth. In dendroclimatological studies, it is generally assumed that the relationships between tree growth and climate factors are age-independent after removing the biological growth trends related to tree age. However, in some multi-aged forest stands, climate-growth relationships may be biased, because the trees′ chronologies respond differently to climate factors depending on their age or size. To assess this age-dependent effect, in this study, we tested the consistency of climate-growth responses in tree-ring series from Picea purpurea trees of different age classes in the Wanglang Nature Reserve, Sichuan province. In the concentrated distribution area of the upper limit, spruces were grouped into three age classes: trees younger than 100 years (young age class, YAC), trees 100—200 years of age (middle age class, MAC), and trees older than 200 years (old age class, OAC). Residual chronologies of the three age classes were built to analyze the climate-growth relationships using correlation and response functions. Correlation analysis indicated that the radial growth of the YAC was significantly and positively correlated with mean monthly temperature of the previous November and current February, June, and September; mean monthly maximum temperature of the previous November and current February and September; and mean monthly minimum temperature of the current February, March, April, June, July, August, and September, and it was significantly and negatively affected by monthly precipitation of the previous December and current June. The ring width of the MAC was positively correlated with mean monthly minimum temperature of the current April and July, and was negatively correlated with monthly precipitation of the current June. The residual chronology of the OAC was significantly and negatively affected by mean monthly temperature of the previous August, mean monthly minimum temperature of the previous August and current September, and monthly precipitation of the previous December. Response analysis showed that radial growth of the YAC was significantly and positively correlated with mean monthly maximum temperature of the current February, and mean monthly minimum temperature of the current February and April. There was no significant correlation between radial growth of the MAC and climate factors. The ring width of the OAC was significantly and negatively affected by monthly precipitation of the previous December.Overall, our results revealed that with an increase in tree age, the sensitivity of spruce to climate factors was reduced. The responses of young spruce were significantly correlated with temperature before the growing season and in the growing season. The chronology of the middle-aged spruce showed a significant and positive correlation with monthly mean minimum temperatures in the current April and July. The ring-width index of old spruce was significantly negatively correlated with monthly mean temperature and monthly mean minimum temperature of the previous August. The “l(fā)ag effect” of high temperatures in the previous growing season was prominent in the old spruce. Spruce within the young and middle-aged groups showed a significant negative correlation with current June precipitation. Adequate precipitation in December was not beneficial for the radial growth of young and old spruce. Our study may demonstrate that future studies should consider the adaptive chronology of age span for spruce in this area to avoid losing or magnifying the climate signals.

    Key Words:Picea purpurea; age; radial growth; climatic factors

    DOI:10.5846/stxb201409121815

    *通訊作者Corresponding author.E-mail: lijq@bjfu.edu.cn

    收稿日期:2014- 09- 12;

    修訂日期:2015- 06- 08

    基金項目:國家林業(yè)公益性行業(yè)科研專項(201404422);青島農(nóng)業(yè)大學(xué)高層次人才科研基金項目(631403)

    趙志江,康東偉,李俊清.川西亞高山不同年齡紫果云杉徑向生長對氣候因子的響應(yīng).生態(tài)學(xué)報,2016,36(1):173- 179.

    Zhao Z J, Kang D W, Li J Q.Age-dependent radial growth responses ofPiceapurpureato climatic factors in the subalpine region of Western Sichuan Province, China.Acta Ecologica Sinica,2016,36(1):173- 179.

    猜你喜歡
    氣候因子年齡
    內(nèi)蒙古不同類型草原植物群落葉片氮、磷含量變化研究
    變小的年齡
    忘記年齡
    TOO YOUNG TO LOCK UP?
    我們都在年齡里生死
    海峽姐妹(2017年3期)2017-04-16 03:06:38
    氣候因子對烤煙質(zhì)量風(fēng)格特色的影響
    年齡歧視
    算年齡
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    天山北坡NDVI對氣候因子響應(yīng)的敏感性分析
    啦啦啦观看免费观看视频高清| 在线免费观看不下载黄p国产 | 国产av一区在线观看免费| 免费人成视频x8x8入口观看| 亚洲黑人精品在线| 精品欧美国产一区二区三| 少妇猛男粗大的猛烈进出视频 | 国产精品三级大全| 国产白丝娇喘喷水9色精品| 久久精品91蜜桃| 国产真实伦视频高清在线观看 | 又爽又黄a免费视频| 亚洲第一电影网av| 别揉我奶头 嗯啊视频| 又黄又爽又刺激的免费视频.| 国产色婷婷99| 亚洲在线观看片| 又黄又爽又刺激的免费视频.| 国产精品一区二区免费欧美| 极品教师在线免费播放| 亚洲va日本ⅴa欧美va伊人久久| 18+在线观看网站| 美女被艹到高潮喷水动态| 窝窝影院91人妻| 成人无遮挡网站| 国产精品一区二区三区四区免费观看 | 最近在线观看免费完整版| 69人妻影院| 久久欧美精品欧美久久欧美| 一级黄片播放器| 亚洲va在线va天堂va国产| 97热精品久久久久久| 一边摸一边抽搐一进一小说| 精品久久久久久久末码| 日本一二三区视频观看| 国产亚洲精品av在线| 欧美xxxx性猛交bbbb| 别揉我奶头~嗯~啊~动态视频| 美女高潮的动态| 久久久国产成人精品二区| 97超级碰碰碰精品色视频在线观看| 桃色一区二区三区在线观看| 国产单亲对白刺激| 99热6这里只有精品| 内地一区二区视频在线| 男女那种视频在线观看| 色5月婷婷丁香| 国产老妇女一区| 午夜免费成人在线视频| 国产成年人精品一区二区| 天堂√8在线中文| 国内毛片毛片毛片毛片毛片| 日本 欧美在线| 亚洲欧美清纯卡通| 欧美xxxx性猛交bbbb| 欧美成人a在线观看| 毛片女人毛片| 一本一本综合久久| 熟女人妻精品中文字幕| 国产精品久久久久久久久免| 联通29元200g的流量卡| 日韩大尺度精品在线看网址| 国产在视频线在精品| 精品国内亚洲2022精品成人| 国产午夜精品论理片| 看十八女毛片水多多多| 亚洲一区二区三区色噜噜| 久久人人精品亚洲av| 黄色女人牲交| 日韩欧美免费精品| 在现免费观看毛片| 99久国产av精品| 韩国av在线不卡| 国内揄拍国产精品人妻在线| 中国美白少妇内射xxxbb| 一本久久中文字幕| 99riav亚洲国产免费| 99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| 免费av不卡在线播放| 亚洲成人久久性| 美女大奶头视频| 国内揄拍国产精品人妻在线| 国产主播在线观看一区二区| 十八禁网站免费在线| 日韩欧美 国产精品| 久久九九热精品免费| 成人无遮挡网站| 国产精品美女特级片免费视频播放器| 精品日产1卡2卡| 亚洲图色成人| 精品久久久久久久久久免费视频| 1024手机看黄色片| 午夜影院日韩av| 自拍偷自拍亚洲精品老妇| 3wmmmm亚洲av在线观看| 日韩欧美一区二区三区在线观看| 国产成人福利小说| 色精品久久人妻99蜜桃| 成人av在线播放网站| 色尼玛亚洲综合影院| 国产探花在线观看一区二区| 在线观看66精品国产| 欧美三级亚洲精品| 成人特级黄色片久久久久久久| 亚洲欧美日韩东京热| 成人综合一区亚洲| 色av中文字幕| 午夜日韩欧美国产| 一区二区三区高清视频在线| 亚洲 国产 在线| 国产aⅴ精品一区二区三区波| 99热这里只有是精品在线观看| 亚洲精品色激情综合| 国产成年人精品一区二区| 日韩av在线大香蕉| 男女视频在线观看网站免费| 最新在线观看一区二区三区| 免费看日本二区| 99热精品在线国产| 亚洲色图av天堂| h日本视频在线播放| 深夜a级毛片| a级毛片a级免费在线| 制服丝袜大香蕉在线| 国产成人影院久久av| 国产69精品久久久久777片| 精品一区二区三区av网在线观看| 中文字幕av在线有码专区| 全区人妻精品视频| 国产91精品成人一区二区三区| 天美传媒精品一区二区| 99在线视频只有这里精品首页| 欧美3d第一页| 精品久久久久久久久久久久久| 色尼玛亚洲综合影院| 亚洲av日韩精品久久久久久密| 日本爱情动作片www.在线观看 | 久久久久久久久久黄片| 国产一区二区激情短视频| 欧美一区二区亚洲| 欧美色视频一区免费| 蜜桃亚洲精品一区二区三区| 99热精品在线国产| 国产精品99久久久久久久久| av在线亚洲专区| 国产精品久久久久久久电影| 国产三级在线视频| 久久人人爽人人爽人人片va| 午夜免费激情av| 九九爱精品视频在线观看| 麻豆av噜噜一区二区三区| 日本 av在线| 亚洲国产欧洲综合997久久,| 色播亚洲综合网| 在现免费观看毛片| 欧美高清性xxxxhd video| 欧美高清成人免费视频www| 老司机深夜福利视频在线观看| 国产白丝娇喘喷水9色精品| av在线观看视频网站免费| 人人妻人人澡欧美一区二区| 国产av在哪里看| 亚洲成人中文字幕在线播放| 成年女人永久免费观看视频| 男插女下体视频免费在线播放| 乱码一卡2卡4卡精品| 看黄色毛片网站| 蜜桃亚洲精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 国产精品精品国产色婷婷| www.www免费av| 精品人妻熟女av久视频| 国产在线精品亚洲第一网站| 精品不卡国产一区二区三区| 国产不卡一卡二| 国产亚洲91精品色在线| 97热精品久久久久久| av国产免费在线观看| 久久久久久久久久黄片| 国产精品久久视频播放| 精品日产1卡2卡| 91在线观看av| 国产精品一区二区三区四区免费观看 | 岛国在线免费视频观看| 久99久视频精品免费| 免费在线观看日本一区| 在线a可以看的网站| 男女啪啪激烈高潮av片| 91av网一区二区| 欧美一区二区国产精品久久精品| 亚洲最大成人手机在线| 高清在线国产一区| 国产精品乱码一区二三区的特点| 永久网站在线| 久久久久久久久久久丰满 | 天堂av国产一区二区熟女人妻| 成人永久免费在线观看视频| 97碰自拍视频| 色综合站精品国产| 99热6这里只有精品| 中文资源天堂在线| 婷婷精品国产亚洲av在线| 国产主播在线观看一区二区| 99九九线精品视频在线观看视频| 亚洲av不卡在线观看| 国产亚洲精品综合一区在线观看| 国产精品综合久久久久久久免费| 国产成人福利小说| 国产精品美女特级片免费视频播放器| 中文字幕免费在线视频6| 中文资源天堂在线| 国产av麻豆久久久久久久| 天堂av国产一区二区熟女人妻| 欧美日韩瑟瑟在线播放| 黄片wwwwww| 观看美女的网站| 人妻久久中文字幕网| 黄色配什么色好看| 99精品久久久久人妻精品| 亚洲欧美精品综合久久99| 国产伦在线观看视频一区| 久久欧美精品欧美久久欧美| av国产免费在线观看| 日韩欧美国产一区二区入口| 搡女人真爽免费视频火全软件 | 成人国产综合亚洲| 色哟哟哟哟哟哟| 搡女人真爽免费视频火全软件 | 男女做爰动态图高潮gif福利片| 亚洲精品久久国产高清桃花| 啦啦啦韩国在线观看视频| 日日啪夜夜撸| 99国产精品一区二区蜜桃av| 久久久久久久久久久丰满 | 久久精品夜夜夜夜夜久久蜜豆| 国产黄色小视频在线观看| 国产69精品久久久久777片| 日韩,欧美,国产一区二区三区 | 国产高清激情床上av| 极品教师在线免费播放| xxxwww97欧美| 亚洲va日本ⅴa欧美va伊人久久| 国产蜜桃级精品一区二区三区| 毛片一级片免费看久久久久 | 婷婷精品国产亚洲av| 一a级毛片在线观看| 亚洲av熟女| 女的被弄到高潮叫床怎么办 | 日本一二三区视频观看| 国产精品电影一区二区三区| 亚洲美女黄片视频| 狠狠狠狠99中文字幕| 最近最新免费中文字幕在线| 国产亚洲精品久久久com| 精品一区二区三区av网在线观看| 老熟妇仑乱视频hdxx| 狠狠狠狠99中文字幕| 免费观看精品视频网站| 色5月婷婷丁香| 亚洲av.av天堂| 高清毛片免费观看视频网站| 免费在线观看日本一区| 精品人妻视频免费看| 毛片女人毛片| 美女被艹到高潮喷水动态| 亚洲精品成人久久久久久| 成人午夜高清在线视频| 又爽又黄a免费视频| 久久热精品热| 狂野欧美激情性xxxx在线观看| 18禁在线播放成人免费| 日日撸夜夜添| 亚洲精品影视一区二区三区av| 亚洲第一电影网av| 国产亚洲精品久久久久久毛片| 国产探花在线观看一区二区| 精品日产1卡2卡| 黄色一级大片看看| 搡女人真爽免费视频火全软件 | 人妻夜夜爽99麻豆av| 亚洲精品粉嫩美女一区| 亚洲欧美激情综合另类| 国产探花极品一区二区| 99热网站在线观看| 久9热在线精品视频| 国产亚洲精品综合一区在线观看| 色噜噜av男人的天堂激情| 日本一二三区视频观看| 亚洲av成人精品一区久久| 国产老妇女一区| 亚洲一级一片aⅴ在线观看| 国产成人福利小说| 日本 欧美在线| 欧美高清性xxxxhd video| 欧美绝顶高潮抽搐喷水| 亚洲精品亚洲一区二区| 美女高潮喷水抽搐中文字幕| 亚洲精品影视一区二区三区av| 国产精品永久免费网站| 我的老师免费观看完整版| 99热只有精品国产| 亚洲成a人片在线一区二区| 色哟哟哟哟哟哟| 久久久久国内视频| 波野结衣二区三区在线| 精品免费久久久久久久清纯| 欧美日韩国产亚洲二区| 乱人视频在线观看| 美女被艹到高潮喷水动态| 精品乱码久久久久久99久播| 在线播放国产精品三级| 久久久久久伊人网av| 蜜桃久久精品国产亚洲av| 亚洲美女黄片视频| 麻豆精品久久久久久蜜桃| 一进一出抽搐gif免费好疼| 亚洲精品久久国产高清桃花| 精品福利观看| 一夜夜www| 免费av观看视频| 国产国拍精品亚洲av在线观看| 真实男女啪啪啪动态图| 亚州av有码| 日韩欧美国产在线观看| 亚洲av.av天堂| 亚洲不卡免费看| 久久国产乱子免费精品| 蜜桃亚洲精品一区二区三区| 日韩欧美精品v在线| 99热6这里只有精品| 两个人的视频大全免费| 欧美色视频一区免费| av.在线天堂| 久久午夜福利片| 亚洲色图av天堂| 日韩欧美免费精品| 欧美日韩国产亚洲二区| 校园人妻丝袜中文字幕| 老司机深夜福利视频在线观看| a级毛片免费高清观看在线播放| 99热精品在线国产| 国产亚洲91精品色在线| 日韩一区二区视频免费看| 久久精品影院6| 欧美一级a爱片免费观看看| 97超视频在线观看视频| 成人美女网站在线观看视频| 国产精品人妻久久久影院| 久久人妻av系列| 久久精品国产清高在天天线| 国产成人a区在线观看| 国产精品,欧美在线| 蜜桃久久精品国产亚洲av| 久久精品久久久久久噜噜老黄 | 黄色女人牲交| 观看免费一级毛片| 国产成年人精品一区二区| 国产高清三级在线| 亚洲国产欧洲综合997久久,| 成人欧美大片| 淫秽高清视频在线观看| 色5月婷婷丁香| 欧美性感艳星| 日韩国内少妇激情av| 一个人免费在线观看电影| 中国美白少妇内射xxxbb| 国产精品精品国产色婷婷| 搡老熟女国产l中国老女人| 俄罗斯特黄特色一大片| 特级一级黄色大片| 狠狠狠狠99中文字幕| 亚洲最大成人中文| 色播亚洲综合网| 国产乱人伦免费视频| 毛片一级片免费看久久久久 | 国内毛片毛片毛片毛片毛片| 少妇熟女aⅴ在线视频| 国产高清有码在线观看视频| 婷婷精品国产亚洲av| 国产精品久久视频播放| 精品一区二区三区视频在线观看免费| 在现免费观看毛片| 大型黄色视频在线免费观看| 麻豆av噜噜一区二区三区| 亚洲成人精品中文字幕电影| 中文字幕免费在线视频6| 久久亚洲精品不卡| 又爽又黄a免费视频| 久久中文看片网| 观看免费一级毛片| 变态另类丝袜制服| 欧美高清成人免费视频www| 最近视频中文字幕2019在线8| 中文资源天堂在线| 成人午夜高清在线视频| 欧美日韩综合久久久久久 | 欧美不卡视频在线免费观看| 午夜老司机福利剧场| 女同久久另类99精品国产91| 日本欧美国产在线视频| 午夜福利高清视频| 两人在一起打扑克的视频| 2021天堂中文幕一二区在线观| 国产成人aa在线观看| 小说图片视频综合网站| 香蕉av资源在线| 大型黄色视频在线免费观看| 久久久久久大精品| 老师上课跳d突然被开到最大视频| 日本色播在线视频| 中文字幕精品亚洲无线码一区| av福利片在线观看| 中文字幕熟女人妻在线| 国内精品久久久久精免费| 中文字幕av成人在线电影| 免费观看的影片在线观看| 不卡视频在线观看欧美| АⅤ资源中文在线天堂| 亚洲成a人片在线一区二区| 国产女主播在线喷水免费视频网站 | a级一级毛片免费在线观看| 丰满的人妻完整版| 久久精品国产99精品国产亚洲性色| 此物有八面人人有两片| 亚洲国产色片| 女同久久另类99精品国产91| 日本a在线网址| 在线观看舔阴道视频| 免费高清视频大片| 亚洲av成人精品一区久久| 久久这里只有精品中国| 蜜桃亚洲精品一区二区三区| 日韩精品青青久久久久久| 亚洲最大成人手机在线| 波野结衣二区三区在线| 岛国在线免费视频观看| 欧美激情国产日韩精品一区| 波多野结衣高清作品| 别揉我奶头 嗯啊视频| 日韩欧美精品v在线| 97超级碰碰碰精品色视频在线观看| 一区二区三区高清视频在线| 春色校园在线视频观看| 国产精品1区2区在线观看.| 欧美成人性av电影在线观看| 日日摸夜夜添夜夜添小说| 免费高清视频大片| 色5月婷婷丁香| 色播亚洲综合网| 亚洲电影在线观看av| 成人性生交大片免费视频hd| 免费看a级黄色片| 精品久久久久久成人av| 欧美黑人欧美精品刺激| 日日摸夜夜添夜夜添av毛片 | 999久久久精品免费观看国产| 99在线视频只有这里精品首页| av天堂在线播放| 极品教师在线视频| 国国产精品蜜臀av免费| 欧美另类亚洲清纯唯美| 国内揄拍国产精品人妻在线| 亚洲在线自拍视频| 欧美性猛交╳xxx乱大交人| 99久久九九国产精品国产免费| 一边摸一边抽搐一进一小说| 大型黄色视频在线免费观看| 午夜福利成人在线免费观看| 深夜精品福利| 尾随美女入室| 精品久久久久久久末码| 白带黄色成豆腐渣| 男人和女人高潮做爰伦理| 婷婷精品国产亚洲av在线| 亚洲人成网站在线播| 亚洲av免费在线观看| 麻豆一二三区av精品| 嫩草影院新地址| 小说图片视频综合网站| 国产成人av教育| 久久国产乱子免费精品| 男女那种视频在线观看| 俄罗斯特黄特色一大片| 女同久久另类99精品国产91| 午夜激情欧美在线| 国产主播在线观看一区二区| 日日摸夜夜添夜夜添小说| 91在线精品国自产拍蜜月| 中文亚洲av片在线观看爽| 狂野欧美白嫩少妇大欣赏| 少妇被粗大猛烈的视频| 高清在线国产一区| 变态另类丝袜制服| 亚洲国产精品合色在线| 好男人在线观看高清免费视频| 国产精品福利在线免费观看| 久久精品国产鲁丝片午夜精品 | 久久久久久久久久成人| 嫁个100分男人电影在线观看| 精品人妻偷拍中文字幕| 亚洲男人的天堂狠狠| www.色视频.com| 成人精品一区二区免费| 啪啪无遮挡十八禁网站| 国产精品综合久久久久久久免费| 人妻久久中文字幕网| 亚洲熟妇中文字幕五十中出| 国产成人一区二区在线| 男人的好看免费观看在线视频| 亚洲av中文av极速乱 | 色哟哟·www| 久久精品国产亚洲av涩爱 | 噜噜噜噜噜久久久久久91| 伦精品一区二区三区| www.www免费av| 日日啪夜夜撸| 亚洲天堂国产精品一区在线| av在线天堂中文字幕| av天堂中文字幕网| x7x7x7水蜜桃| 免费av毛片视频| 中文字幕av在线有码专区| 日本 欧美在线| 在线天堂最新版资源| 久久久久性生活片| 我的女老师完整版在线观看| 91狼人影院| 春色校园在线视频观看| 国产精品无大码| 99热这里只有是精品50| 天美传媒精品一区二区| 三级国产精品欧美在线观看| 老司机午夜福利在线观看视频| 国产精品免费一区二区三区在线| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 国产精品一区二区性色av| 精品无人区乱码1区二区| 日韩欧美一区二区三区在线观看| 亚洲成av人片在线播放无| 久久天躁狠狠躁夜夜2o2o| 国产一区二区三区视频了| 如何舔出高潮| 国产精品三级大全| 国产精品98久久久久久宅男小说| 很黄的视频免费| 久久久久久伊人网av| 搡老熟女国产l中国老女人| 亚洲欧美清纯卡通| 波多野结衣高清无吗| 亚洲av电影不卡..在线观看| 久久精品影院6| 老熟妇乱子伦视频在线观看| 午夜激情欧美在线| 久久婷婷人人爽人人干人人爱| 欧美高清成人免费视频www| 99riav亚洲国产免费| 国产三级在线视频| 九九爱精品视频在线观看| 露出奶头的视频| 亚洲人成网站在线播放欧美日韩| xxxwww97欧美| 亚洲18禁久久av| 免费在线观看日本一区| 国产91精品成人一区二区三区| 欧美+日韩+精品| 麻豆国产av国片精品| 国产不卡一卡二| 中文字幕高清在线视频| 亚洲va在线va天堂va国产| 国产真实伦视频高清在线观看 | 免费看美女性在线毛片视频| 午夜福利在线观看免费完整高清在 | 99久久久亚洲精品蜜臀av| 精品午夜福利视频在线观看一区| 国产精品国产三级国产av玫瑰| 一个人看视频在线观看www免费| 在线免费十八禁| 一区福利在线观看| 日韩欧美在线乱码| 日本 欧美在线| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 少妇丰满av| 乱码一卡2卡4卡精品| 色在线成人网| 色av中文字幕| 男女下面进入的视频免费午夜| 国产 一区 欧美 日韩| 国产精品三级大全| 韩国av一区二区三区四区| 嫩草影院新地址| 国产亚洲精品久久久久久毛片| 国产成人aa在线观看| av国产免费在线观看| 国产 一区 欧美 日韩| 亚洲av日韩精品久久久久久密| 日韩,欧美,国产一区二区三区 | 亚洲av成人精品一区久久| 露出奶头的视频| 色综合色国产| 国产精品嫩草影院av在线观看 | 亚洲国产精品成人综合色| 永久网站在线| 久久99热6这里只有精品| 久久精品国产亚洲av香蕉五月| av专区在线播放| 国产高清有码在线观看视频| 午夜福利在线观看免费完整高清在 | 亚洲成人免费电影在线观看| 少妇人妻精品综合一区二区 | 老司机福利观看| 十八禁网站免费在线| 熟妇人妻久久中文字幕3abv|