• <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)的敏感性分析
    亚洲精品国产一区二区精华液| 国产欧美亚洲国产| 精品国产超薄肉色丝袜足j| 激情视频va一区二区三区| 午夜精品国产一区二区电影| 建设人人有责人人尽责人人享有的| 午夜激情av网站| 咕卡用的链子| 一级a爱视频在线免费观看| 欧美日韩av久久| 日韩av在线免费看完整版不卡| 精品国产乱码久久久久久小说| 一级毛片 在线播放| 又粗又硬又长又爽又黄的视频| √禁漫天堂资源中文www| 久久人妻熟女aⅴ| 亚洲,欧美精品.| 99精国产麻豆久久婷婷| 亚洲欧美精品自产自拍| 久久精品亚洲av国产电影网| 亚洲欧美日韩高清在线视频 | 中文欧美无线码| 亚洲色图 男人天堂 中文字幕| 国产成人免费观看mmmm| 久久综合国产亚洲精品| 美女国产高潮福利片在线看| 亚洲图色成人| 91精品三级在线观看| 免费看十八禁软件| 精品欧美一区二区三区在线| 少妇猛男粗大的猛烈进出视频| 精品国产超薄肉色丝袜足j| 丁香六月天网| 另类亚洲欧美激情| 亚洲一区二区三区欧美精品| 18禁裸乳无遮挡动漫免费视频| 国产视频首页在线观看| 大片免费播放器 马上看| 妹子高潮喷水视频| 另类亚洲欧美激情| 狠狠婷婷综合久久久久久88av| 日韩av不卡免费在线播放| 亚洲国产欧美在线一区| 久久精品国产综合久久久| 丝袜脚勾引网站| 亚洲av在线观看美女高潮| 黄片小视频在线播放| 91麻豆精品激情在线观看国产 | 美国免费a级毛片| 超碰成人久久| 亚洲精品一区蜜桃| 国产亚洲精品第一综合不卡| 亚洲精品久久成人aⅴ小说| 香蕉丝袜av| 国产成人av教育| 一区二区三区四区激情视频| 飞空精品影院首页| 免费在线观看完整版高清| 丝袜喷水一区| 国产免费又黄又爽又色| 国产成人精品在线电影| 这个男人来自地球电影免费观看| 国产男女内射视频| 美国免费a级毛片| 亚洲精品日本国产第一区| 欧美乱码精品一区二区三区| 久久久久久久国产电影| 三上悠亚av全集在线观看| 欧美黄色片欧美黄色片| 欧美日韩视频精品一区| 少妇的丰满在线观看| svipshipincom国产片| 亚洲成人手机| 久久久久久久大尺度免费视频| 日本91视频免费播放| 国产精品av久久久久免费| 久久99精品国语久久久| 日本猛色少妇xxxxx猛交久久| 国产成人a∨麻豆精品| 97精品久久久久久久久久精品| 亚洲av在线观看美女高潮| 男人添女人高潮全过程视频| 夫妻性生交免费视频一级片| 一区在线观看完整版| 国产亚洲欧美精品永久| 高清欧美精品videossex| 亚洲精品国产区一区二| 在线观看免费视频网站a站| 欧美精品人与动牲交sv欧美| 我的亚洲天堂| 97在线人人人人妻| cao死你这个sao货| av网站免费在线观看视频| 成人黄色视频免费在线看| 久久久久国产一级毛片高清牌| 久久青草综合色| 五月天丁香电影| 亚洲专区中文字幕在线| 日韩,欧美,国产一区二区三区| 久久天躁狠狠躁夜夜2o2o | h视频一区二区三区| 久久久国产精品麻豆| 国产黄色免费在线视频| 欧美中文综合在线视频| 欧美乱码精品一区二区三区| 蜜桃国产av成人99| 亚洲精品一二三| 欧美激情高清一区二区三区| 51午夜福利影视在线观看| 国产一区二区三区av在线| 精品欧美一区二区三区在线| 国产成人啪精品午夜网站| 国精品久久久久久国模美| 老司机亚洲免费影院| 中文字幕制服av| kizo精华| 国产精品 国内视频| 最新在线观看一区二区三区 | 国产一级毛片在线| a级片在线免费高清观看视频| 男女之事视频高清在线观看 | 搡老乐熟女国产| 国产99久久九九免费精品| 精品福利观看| 免费少妇av软件| 午夜福利在线免费观看网站| bbb黄色大片| 亚洲国产中文字幕在线视频| 免费不卡黄色视频| 亚洲三区欧美一区| 女人精品久久久久毛片| 你懂的网址亚洲精品在线观看| 中文字幕亚洲精品专区| 啦啦啦中文免费视频观看日本| 久久久精品区二区三区| 9色porny在线观看| 一级毛片电影观看| 免费看十八禁软件| 交换朋友夫妻互换小说| 婷婷色麻豆天堂久久| 欧美在线一区亚洲| 久久精品久久久久久噜噜老黄| 精品一区二区三区av网在线观看 | 一边摸一边做爽爽视频免费| 人人妻,人人澡人人爽秒播 | 韩国精品一区二区三区| 好男人电影高清在线观看| 伊人久久大香线蕉亚洲五| 成人国产一区最新在线观看 | 人体艺术视频欧美日本| 成人亚洲欧美一区二区av| 成人亚洲欧美一区二区av| 免费观看人在逋| 人成视频在线观看免费观看| 中国美女看黄片| 亚洲av欧美aⅴ国产| 欧美日韩成人在线一区二区| 后天国语完整版免费观看| 国产精品三级大全| 日韩免费高清中文字幕av| 午夜视频精品福利| 国产91精品成人一区二区三区 | 日韩中文字幕欧美一区二区 | 黑丝袜美女国产一区| 国产一区有黄有色的免费视频| 高清欧美精品videossex| 99国产精品99久久久久| 午夜福利视频精品| 美女福利国产在线| 久久久久久久精品精品| 丝袜美足系列| 老司机影院成人| 九色亚洲精品在线播放| 美女福利国产在线| 啦啦啦在线观看免费高清www| 十分钟在线观看高清视频www| 操美女的视频在线观看| 国产国语露脸激情在线看| 国产极品粉嫩免费观看在线| 免费观看a级毛片全部| 国产一区二区三区av在线| 操出白浆在线播放| 亚洲av男天堂| 欧美成狂野欧美在线观看| 亚洲一区二区三区欧美精品| 欧美精品一区二区免费开放| 日本vs欧美在线观看视频| 国产成人91sexporn| 国产精品 欧美亚洲| 青春草亚洲视频在线观看| 久久综合国产亚洲精品| netflix在线观看网站| 国产在视频线精品| 一边摸一边抽搐一进一出视频| 91精品伊人久久大香线蕉| 伊人亚洲综合成人网| 国产亚洲精品久久久久5区| 高清黄色对白视频在线免费看| 久久久亚洲精品成人影院| 视频在线观看一区二区三区| 亚洲av日韩精品久久久久久密 | 亚洲色图 男人天堂 中文字幕| 最近中文字幕2019免费版| 日本黄色日本黄色录像| 久久久欧美国产精品| 黄片小视频在线播放| 免费看十八禁软件| 久久热在线av| 成年av动漫网址| 国产免费现黄频在线看| 亚洲人成电影免费在线| 热re99久久国产66热| 乱人伦中国视频| 久久国产精品大桥未久av| 国产成人av教育| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| 国产一区二区激情短视频 | 亚洲第一av免费看| 色精品久久人妻99蜜桃| 亚洲情色 制服丝袜| 欧美成人午夜精品| kizo精华| 亚洲成av片中文字幕在线观看| 老司机影院成人| 99久久99久久久精品蜜桃| 秋霞在线观看毛片| 少妇被粗大的猛进出69影院| 高清黄色对白视频在线免费看| 人成视频在线观看免费观看| xxx大片免费视频| 国产精品久久久久成人av| 好男人电影高清在线观看| 精品人妻在线不人妻| 亚洲视频免费观看视频| 黑丝袜美女国产一区| 婷婷色综合大香蕉| 少妇的丰满在线观看| 男女边吃奶边做爰视频| 国产福利在线免费观看视频| 国产三级黄色录像| 人成视频在线观看免费观看| 国产精品久久久久久人妻精品电影 | 亚洲专区中文字幕在线| 一本综合久久免费| 秋霞在线观看毛片| 老司机靠b影院| 日韩一区二区三区影片| 日日夜夜操网爽| 啦啦啦中文免费视频观看日本| 一区二区三区乱码不卡18| 亚洲av男天堂| 欧美久久黑人一区二区| 精品一区二区三区av网在线观看 | 亚洲国产最新在线播放| 国产有黄有色有爽视频| 性高湖久久久久久久久免费观看| 欧美日本中文国产一区发布| 亚洲中文av在线| 亚洲av电影在线进入| 国产成人精品久久二区二区91| 久久精品亚洲av国产电影网| 免费不卡黄色视频| 爱豆传媒免费全集在线观看| 欧美另类一区| 大片电影免费在线观看免费| 一边亲一边摸免费视频| 在线观看免费高清a一片| 咕卡用的链子| 久久青草综合色| 9热在线视频观看99| 久久天堂一区二区三区四区| 国产一区二区在线观看av| 午夜福利视频精品| 美女中出高潮动态图| 母亲3免费完整高清在线观看| 久久综合国产亚洲精品| 韩国精品一区二区三区| 久久天躁狠狠躁夜夜2o2o | 秋霞在线观看毛片| 国产伦人伦偷精品视频| 久久国产精品人妻蜜桃| 亚洲成人免费av在线播放| 电影成人av| 黑人猛操日本美女一级片| 在线 av 中文字幕| 亚洲国产日韩一区二区| av网站在线播放免费| 亚洲成人免费电影在线观看 | 国产成人免费观看mmmm| 在现免费观看毛片| 亚洲中文av在线| 国产精品人妻久久久影院| av不卡在线播放| 飞空精品影院首页| 青草久久国产| 最近手机中文字幕大全| 亚洲欧美色中文字幕在线| 久久天堂一区二区三区四区| 成人手机av| 免费高清在线观看视频在线观看| 久久99精品国语久久久| 久久99一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 91麻豆av在线| 国产高清国产精品国产三级| 黄片播放在线免费| 超碰成人久久| 亚洲少妇的诱惑av| 久久精品人人爽人人爽视色| 日本黄色日本黄色录像| 电影成人av| 七月丁香在线播放| 亚洲 欧美一区二区三区| 国产精品99久久99久久久不卡| 国产爽快片一区二区三区| 亚洲人成网站在线观看播放| 国产熟女欧美一区二区| 男女边摸边吃奶| 黄色毛片三级朝国网站| 国产精品麻豆人妻色哟哟久久| 亚洲国产欧美一区二区综合| 亚洲欧美日韩高清在线视频 | 国产成人欧美| 久久99一区二区三区| 亚洲欧美日韩高清在线视频 | 制服人妻中文乱码| 999精品在线视频| 两个人免费观看高清视频| 人人澡人人妻人| 三上悠亚av全集在线观看| 成年人免费黄色播放视频| 久久天躁狠狠躁夜夜2o2o | 中文字幕人妻熟女乱码| 桃花免费在线播放| 欧美日韩亚洲国产一区二区在线观看 | av在线app专区| av在线app专区| 久久鲁丝午夜福利片| 在线看a的网站| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品久久午夜乱码| 亚洲午夜精品一区,二区,三区| 亚洲欧美色中文字幕在线| 人体艺术视频欧美日本| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 纵有疾风起免费观看全集完整版| 欧美另类一区| 日韩 亚洲 欧美在线| 亚洲精品美女久久久久99蜜臀 | 中国美女看黄片| 亚洲国产欧美网| 国产不卡av网站在线观看| 亚洲国产av新网站| 亚洲五月色婷婷综合| 人人妻人人澡人人爽人人夜夜| 国产精品三级大全| 中文字幕制服av| 国产成人精品久久二区二区91| 高潮久久久久久久久久久不卡| 水蜜桃什么品种好| 亚洲成色77777| 欧美另类一区| 亚洲av成人不卡在线观看播放网 | 久久中文字幕一级| 一区二区三区精品91| 国产三级黄色录像| av在线app专区| 久热这里只有精品99| 久久精品人人爽人人爽视色| 国产精品三级大全| 欧美+亚洲+日韩+国产| 亚洲激情五月婷婷啪啪| 女人爽到高潮嗷嗷叫在线视频| 免费看不卡的av| 在线av久久热| 国产精品免费视频内射| 国产高清不卡午夜福利| 亚洲av电影在线观看一区二区三区| 国产一区有黄有色的免费视频| 日韩av在线免费看完整版不卡| 午夜激情久久久久久久| 在线 av 中文字幕| 国产精品麻豆人妻色哟哟久久| 免费观看av网站的网址| 女人精品久久久久毛片| 高清不卡的av网站| 欧美在线黄色| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 国产日韩欧美亚洲二区| 亚洲五月婷婷丁香| 亚洲午夜精品一区,二区,三区| 成人亚洲精品一区在线观看| 久久久精品国产亚洲av高清涩受| 欧美国产精品va在线观看不卡| 在线精品无人区一区二区三| 少妇的丰满在线观看| 亚洲国产欧美日韩在线播放| 亚洲av欧美aⅴ国产| 99九九在线精品视频| 亚洲,欧美精品.| 中文字幕色久视频| 久久久久久人人人人人| 久久久亚洲精品成人影院| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 欧美日韩综合久久久久久| 日本猛色少妇xxxxx猛交久久| 中文字幕亚洲精品专区| 亚洲av国产av综合av卡| www日本在线高清视频| 日韩,欧美,国产一区二区三区| 久久99精品国语久久久| 九色亚洲精品在线播放| 少妇裸体淫交视频免费看高清 | 欧美日韩黄片免| 男女边摸边吃奶| 我要看黄色一级片免费的| kizo精华| 汤姆久久久久久久影院中文字幕| 欧美人与性动交α欧美软件| 波野结衣二区三区在线| 久久99热这里只频精品6学生| 精品人妻在线不人妻| 91麻豆av在线| 日本91视频免费播放| 精品国产乱码久久久久久小说| a级毛片黄视频| 久久鲁丝午夜福利片| 久久国产精品人妻蜜桃| 国产午夜精品一二区理论片| 大码成人一级视频| 男人舔女人的私密视频| 一本—道久久a久久精品蜜桃钙片| 丝袜在线中文字幕| 亚洲一区二区三区欧美精品| 老汉色∧v一级毛片| 国产精品 国内视频| 丰满少妇做爰视频| 男女床上黄色一级片免费看| 成人影院久久| 久久狼人影院| 自拍欧美九色日韩亚洲蝌蚪91| 女人被躁到高潮嗷嗷叫费观| 国产色视频综合| 亚洲国产精品成人久久小说| 精品国产乱码久久久久久男人| 脱女人内裤的视频| 国产欧美日韩一区二区三 | 国产在线免费精品| 欧美久久黑人一区二区| 91麻豆精品激情在线观看国产 | 国产免费又黄又爽又色| 桃花免费在线播放| 乱人伦中国视频| 国产极品粉嫩免费观看在线| 欧美成人精品欧美一级黄| 国产真人三级小视频在线观看| 精品国产乱码久久久久久小说| 国产成人精品久久久久久| 午夜免费男女啪啪视频观看| videosex国产| 99久久精品国产亚洲精品| 婷婷色综合www| 两个人免费观看高清视频| 十八禁人妻一区二区| 国产精品一区二区免费欧美 | 美女脱内裤让男人舔精品视频| 亚洲av男天堂| 亚洲精品成人av观看孕妇| 黄频高清免费视频| netflix在线观看网站| 国产伦理片在线播放av一区| 亚洲精品国产一区二区精华液| 精品人妻一区二区三区麻豆| 国产97色在线日韩免费| 亚洲欧美清纯卡通| 国产精品一区二区在线观看99| 久久久久久久精品精品| 欧美+亚洲+日韩+国产| 日韩人妻精品一区2区三区| 成年人午夜在线观看视频| 精品免费久久久久久久清纯 | 少妇人妻 视频| 免费看av在线观看网站| 亚洲一码二码三码区别大吗| 超碰成人久久| 国产伦人伦偷精品视频| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 女人高潮潮喷娇喘18禁视频| 精品国产一区二区三区四区第35| 国产精品国产三级专区第一集| 亚洲精品日本国产第一区| 亚洲男人天堂网一区| 欧美国产精品一级二级三级| 久久精品aⅴ一区二区三区四区| 一边亲一边摸免费视频| 久久狼人影院| 日韩av免费高清视频| 新久久久久国产一级毛片| 极品人妻少妇av视频| 自线自在国产av| 亚洲欧美成人综合另类久久久| 欧美 日韩 精品 国产| 69精品国产乱码久久久| 久久久久国产精品人妻一区二区| 黄频高清免费视频| 亚洲精品日本国产第一区| 1024视频免费在线观看| 九色亚洲精品在线播放| 亚洲人成电影观看| 啦啦啦啦在线视频资源| 在线观看免费日韩欧美大片| 国产视频一区二区在线看| 国产女主播在线喷水免费视频网站| 99国产精品一区二区三区| 男女免费视频国产| av不卡在线播放| 高清黄色对白视频在线免费看| 丝袜脚勾引网站| 色精品久久人妻99蜜桃| 一区二区三区激情视频| 十八禁人妻一区二区| 丰满饥渴人妻一区二区三| 香蕉国产在线看| 久久久久久人人人人人| 少妇人妻 视频| 国产亚洲精品第一综合不卡| 久久久久视频综合| 国产精品麻豆人妻色哟哟久久| 狂野欧美激情性xxxx| 久久狼人影院| 国产视频首页在线观看| 亚洲av美国av| xxxhd国产人妻xxx| 午夜福利免费观看在线| 免费在线观看视频国产中文字幕亚洲 | 国产成人精品在线电影| 亚洲国产欧美网| 免费观看av网站的网址| 99国产精品一区二区三区| 欧美日韩综合久久久久久| 亚洲欧美清纯卡通| 亚洲一码二码三码区别大吗| 水蜜桃什么品种好| a级毛片在线看网站| 国产黄色视频一区二区在线观看| 天天躁夜夜躁狠狠久久av| 另类精品久久| 美女主播在线视频| 亚洲伊人色综图| 一本色道久久久久久精品综合| 精品视频人人做人人爽| 2018国产大陆天天弄谢| 宅男免费午夜| 亚洲精品久久午夜乱码| 夫妻性生交免费视频一级片| 亚洲av片天天在线观看| 亚洲国产日韩一区二区| 国产成人a∨麻豆精品| 丁香六月天网| 丰满人妻熟妇乱又伦精品不卡| 热99久久久久精品小说推荐| 久久久久精品国产欧美久久久 | 精品少妇一区二区三区视频日本电影| 最新的欧美精品一区二区| 久久国产精品影院| 91国产中文字幕| 欧美亚洲日本最大视频资源| 又大又黄又爽视频免费| 男女午夜视频在线观看| 啦啦啦在线观看免费高清www| 宅男免费午夜| 亚洲精品av麻豆狂野| 国产精品九九99| 91九色精品人成在线观看| 老鸭窝网址在线观看| 我要看黄色一级片免费的| 精品福利永久在线观看| 97人妻天天添夜夜摸| 色精品久久人妻99蜜桃| 性色av一级| 国产麻豆69| 中文字幕制服av| 一区二区av电影网| 精品人妻1区二区| 成人三级做爰电影| 男人添女人高潮全过程视频| 男女国产视频网站| 国产成人一区二区在线| av国产精品久久久久影院| 人体艺术视频欧美日本| 亚洲一区中文字幕在线| 汤姆久久久久久久影院中文字幕| 久久久久久人人人人人| bbb黄色大片| 日本欧美视频一区| 国产高清国产精品国产三级| 亚洲久久久国产精品| 中国国产av一级| 久久精品国产亚洲av涩爱| 久久国产精品人妻蜜桃| 少妇猛男粗大的猛烈进出视频| 亚洲欧美一区二区三区黑人| 777久久人妻少妇嫩草av网站| 欧美精品一区二区免费开放| 亚洲av成人不卡在线观看播放网 | 国产成人免费观看mmmm| 欧美av亚洲av综合av国产av| 欧美大码av| 性高湖久久久久久久久免费观看| 99国产精品99久久久久|