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

    天山中西段不同地區(qū)雪嶺云杉徑向生長對氣候變暖的響應差異

    2017-11-10 10:43:58張艷靜鄭宏偉于瑞德楊美琳石冰冰
    植物研究 2017年3期
    關鍵詞:烏蘇年表氣候因子

    張艷靜 鄭宏偉 于瑞德 楊美琳 石冰冰 王 蕾

    (1.中國科學院新疆生態(tài)與地理研究所荒漠環(huán)境研究室,烏魯木齊 830011; 2.中國科學院大學,北京 100049; 3.新疆林業(yè)科學院現代林業(yè)研究所,烏魯木齊 830011)

    天山中西段不同地區(qū)雪嶺云杉徑向生長對氣候變暖的響應差異

    張艷靜1,2鄭宏偉1*于瑞德1楊美琳1石冰冰1,2王 蕾3

    (1.中國科學院新疆生態(tài)與地理研究所荒漠環(huán)境研究室,烏魯木齊 830011;2.中國科學院大學,北京 100049;3.新疆林業(yè)科學院現代林業(yè)研究所,烏魯木齊 830011)

    為研究烏蘇和石河子地區(qū)雪嶺云杉(Piceaschrenkiana)徑向生長對氣候變暖的響應差異狀況,利用響應函數及滑動相關分析等樹輪氣候學方法,研究了2個地區(qū)森林下限雪嶺云杉徑向生長與氣候因子的響應關系。結果表明:烏蘇和石河子地區(qū)雪嶺云杉生長均與生長季的氣候因子呈顯著相關關系。但烏蘇和石河子地區(qū)雪嶺云杉生長對當地氣候因子的響應也存在差異,烏蘇地區(qū)雪嶺云杉徑向生長與上一年7、9月平均氣溫呈顯著負相關,與當年8月降水、上一年9月及當年8月平均相對濕度和上一年7~10月的scPDSI呈顯著正相關。石河子地區(qū)雪嶺云杉徑向生長與當年1月平均濕度和當年1、2月的scPDSI呈顯著正相關。此外,烏蘇地區(qū)樹輪寬度指數具有與溫度升高而下降的“分離現象”,而石河子地區(qū)樹輪寬度指數具有與溫度升高相一致的趨勢。升溫導致的水分脅迫是造成烏蘇地區(qū)雪嶺云杉徑向生長與溫度變化趨勢相反的重要因素。生長季的溫度和降水的增加對石河子地區(qū)雪嶺云杉徑向生長有促進作用。烏蘇地區(qū)雪嶺云杉徑向生長對5~7月平均氣溫和降水量變化敏感性上升;石河子地區(qū)雪嶺云杉徑向生長對5~7月降水量變化敏感性上升,而對5~7月平均氣溫變化敏感性下降。結果顯示,氣候變化的區(qū)域差異是造成烏蘇和石河子地區(qū)雪嶺云杉徑向生長趨勢不同以及各年表序列對溫度和降水變化敏感性變化的主要氣候因素。

    樹木年輪;徑向生長;氣候變暖;敏感性;天山;響應差異

    樹木年輪以其定年準確、分辨率高、連續(xù)性強以及重建精度好等特點,成為研究過去全球變化(PAGES)的重要技術途徑之一[1]。大多數樹木年輪的寬窄變化能夠記錄樹木生長時的氣候變化信息,探討樹木徑向生長的主要限制因子不僅有助于解釋植被對氣候變化響應的規(guī)律性,也可以為預測未來氣候變化及模擬未來氣候變化條件下植被生長的潛在變化[2~4]。然而,年輪寬度與氣候要素之間存在著復雜的關系,不同地域不同氣候因子條件下,同一樹種樹木生長存在差異[5],氣候條件的變化也將會影響樹木生長與氣候要素的關系,并進一步影響樹木徑向生長的過程[6]。這一結果已得到國內外大量研究的證實,如自20世紀80年代以來,北半球高緯度高海拔地區(qū)樹木生長對氣溫升高表現出敏感度降低的現象[7]。此外,在阿拉斯加地區(qū),D’Arrigo等[8]報道了白云杉(Piceaglauca)生長與氣溫的關系在1970年前后出現差別,隨著氣溫升高,生長與氣溫的關系減弱。在意大利地區(qū),Carrer等[9]對意大利北部林線處歐洲落葉松(Larixdecidua)樹木生長進行了研究,發(fā)現樹木生長在19~20世紀對氣候變化的敏感性發(fā)生了轉變。在長白山地區(qū),于健等[10]研究發(fā)現長白山落葉松(Larixolgensis)生長與氣溫的關系在1989年前后出現了不同的響應特征,發(fā)生了響應模式的改變。王淼等[11]報道了大氣增溫有利于長白山地區(qū)闊葉紅松(Pinuskoraiensis)主要樹種的徑向生長。同時,一些研究表明不同坡向、樹種、海拔和緯度樹木生長會對氣候變暖存在明顯地響應差異,如Leonelli等[12]研究表明,高海拔地區(qū)南坡和西坡瑞士石松(Pinuscembra)生長對氣溫的敏感性在升溫后下降,而北坡瑞士石松生長對氣溫敏感性卻呈增加趨勢。李廣起等[13]研究發(fā)現長白山紅松年輪寬度具有與溫度升高相一致的趨勢,而魚鱗云杉(Piceajezoensisvar.komarovii)年輪寬度則出現隨溫度升高而下降的“分離現象”。郭濱德等[14]研究了川西不同坡向紫果云杉(Piceapurpurea)和岷江冷杉(Abiesfaxoniana)樹木生長對氣候變暖的響應差異,發(fā)現隨著氣溫升高云杉生長顯著增加,而冷杉生長降低,但不顯著。此外,張文濤等[15]研究指出氣溫的升高不僅影響蘆芽山不同海拔的華北落葉松(Larixprincipis-rupprechtii)與氣溫因子關系的海拔差異,而且對其徑向生長的海拔差異也有影響。劉敏等[16]研究了不同緯度紅松徑向生長趨勢及其與氣候因子的關系,掲示了不同緯度紅松徑向生長對氣候因子響應的差異及其氣候變化影響下紅松的動態(tài)特征、適應性及敏感性。而雪嶺云杉(Piceaschrenkiana)作為天山森林的主要建群樹種之一,不同海拔雪嶺云杉對氣候因子的響應也存在差異[17~19]。但是,氣候變暖對不同地區(qū)雪嶺云杉生長的影響及樹木生長對氣候變化的時間動態(tài)響應特征是否均存在差異?目前,天山地區(qū)雪嶺云杉有關這些問題的研究較少,因此,開展烏蘇和石河子地區(qū)的不同區(qū)域樹木生長對氣候變暖響應差異的研究將是對不同區(qū)位氣候條件下森林植被對氣候因子響應研究的有益補充。

    本文將利用樹木年輪資料分析樹木徑向生長與氣候因子的關系,掲示烏蘇和石河子地區(qū)雪嶺云杉徑向生長的主要限制因子,并對比分析不同地區(qū)(烏蘇和石河子)樹木生長對氣候因子的不同響應。在氣候變暖背景下,我們探究天山中西部不同地區(qū)森林下限雪嶺云杉樹木生長對氣候變化的響應差異,確定氣候變化前后樹木生長與氣候因子相關關系的變化,以期了解氣候變化對天山中西部不同地區(qū)雪嶺云杉徑向生長的影響,從而對天山地區(qū)氣候變化背景下的森林植被管理提供理論依據。

    1 研究區(qū)概況

    烏蘇與石河子地區(qū)位于天山北坡中西部,準噶爾盆地南緣,地理坐標為44°09′~43°47′N,84°21′~85°55′E。天山北坡中西部屬溫帶大陸性半干旱氣候,冬季寒冷,夏季炎熱,春秋季較短,冬夏溫差大,日照時間長,降水量少,蒸發(fā)量大。垂直地帶性比較明顯,海拔1 800~2 700 m分布著大量天然森林,以雪嶺云杉為主,林內植被有灌木和草類,土壤為黑褐色。這兩個地區(qū)的氣候總體特征相似,但也存在一些不同,烏蘇地區(qū)溫度高于石河子地區(qū),而其降水低于石河子地區(qū)。總體而言,石河子地區(qū)濕度高于烏蘇地區(qū)。2個樣地的雪嶺云杉均有成片的天然林,林內沒有人為的擾動和自然破壞。雪嶺云杉是研究區(qū)的優(yōu)勢種,常組成純林。該樹種為地帶性植被,屬于耐蔭樹種,壽命長,喜歡生長在氣候濕潤的中山帶陰坡、半陰坡及山中河谷[19]。雪嶺云杉的年輪可形成明顯界限,每一年的年輪較易辨認,對氣候變化有很好的敏感性,已被證實是樹輪氣候學分析的理想樹種[20]。

    2 材料和方法

    2.1 樣本采集

    2015年8月進行了樣本的采集工作,采樣點分布圖如圖1所示。按照樹木年輪氣候學的基本要求、原理和方法,在樹木生長過程對氣候變化更為敏感的森林下限采集樣本[20],兩個采樣點都分別隨機采集了23(46根樹芯)、22棵樹(48根樹芯),每棵樹用生長錐在胸徑位置從南北、東西方向取兩個完好的樹芯作為樣本,樣本樹都是健康樹,采樣點的基本情況見表1。

    圖1 樹輪采樣點位置示意圖Fig.1 Location of the tree-ring sampling site

    采樣點Samplingsites樹種Treespecies緯度Latitude(N)經度Longitude(E)海拔Elevation(m)樣本數(株)Samplenumber坡向Aspect坡度Slope(°)烏蘇Wusu雪嶺云杉Piceaschrenkiana44°09′84°21′190442(23)NE30.9石河子Shihezi雪嶺云杉Piceaschrenkiana43°47′85°55′228038(22)NE?NW35.9

    2.2 建立年表

    樣本帶回實驗室后,進行加工處理:自然干燥、固定、打磨,直至在顯微鏡下能夠清晰的辨認出早材、晚材細胞及輪界分明為止[21]。在顯微鏡下用畫骨架法對樣本進行目視交叉定年,然后使用測量精度為0.01 mm的LINTAB寬度測量儀進行輪寬測量[22]。為保證定年和測量的準確性,用國際年輪庫的COFECHA程序[23]進行交叉定年質量檢驗,確保每一年輪具有準確的日歷年齡,本研究保留了年輪序列之間高度相關的樹芯數據,即42根樹芯(烏蘇)以及38根樹芯(石河子)參與制定兩地區(qū)各自最終的年表。

    樹輪寬度年表的建立是由國際年輪庫ARSTAN程序[24]完成,所有樣芯的去趨勢方法均選取相對保守的負指數函數,以去除樹木本身生長趨勢與非氣候因素對樹木生長的影響,然后用雙權重平均法對測量的輪寬曲線進行標準化,并得出標準化年表(STD)、差值年表(RES)和自回歸標準化年表(ARS)。由于差值年表中保存了較高質量的高頻信息[25],最終選取差值年表作為分析對象。

    2.3 氣候資料

    氣候資料來自烏蘇氣象站(44°25′N,84°40′E,海拔478.7 m)和石河子氣象站(44°19′N,86°03′E,海拔442.9 m)1954~2012年的氣象數據,此數據來自中國氣象科學數據共享網服務網(http://cdc.cma.gov.cn)。采用Mann-Kandell方法[26]和double-mass analysis方法[27]對2個氣象站的氣溫和降水進行了均一性檢驗。檢驗結果表明:2個氣象站的氣候數據無隨機突變和明顯不均勻分布情況,代表了自然界的主要氣候變化特征。所選指標包含2個氣象站點的月平均氣溫、月均降水量和月平均相對濕度(圖2),另外,本文還采用了更能精確評估干旱變化的基于改進帕爾默干旱指數(scPDSI),數據來源于英國東英格利亞大學氣候研究中心(climatic research unit,CRU),可從英國氣象數據中心下載(http://www.badc.rl.ac.uk)[28]。scPDSI數據時間段為1961~2012年,空間分辨率為0.5°×0.5°,選取離烏蘇和石河子采樣點最近格點(44°15′N,84°45′E;44°15′N,86°15′E)的scPDSI數據。

    圖2 烏蘇(A、B)和石河子地區(qū)(C、D)的月均溫度、月均降水量、平均相對濕度和scPDSI指數Fig.2 Monthly mean precipitation,monthly mean temperature,relative humidity and scPDSI of Wusu(A,B) and Shihezi(C,D)

    氣象資料顯示烏蘇地區(qū)年平均溫度以1969年為分界點劃分為平穩(wěn)期(stable)和升溫期(warmer),年平均溫度從5.91℃(1969年)增至8.15℃(2012年),升溫速度大于2.5℃·50 a-1,而石河子地區(qū)自1954年以來處于升溫期,年平均溫度從5.75℃(1954年)增至7.17℃(2012年),升溫速度大于1.2℃·50 a-1(圖3)。烏蘇和石河子地區(qū)年降水量表現出緩慢增加的變化趨勢,兩研究區(qū)年降水量最小和最大分別為84.7 mm(1965年)、124.9 mm(1978年)和339.1 mm(2010年)、351.8 mm(2010)。對比分析發(fā)現,烏蘇和石河子地區(qū)自1969年以來都屬于明顯的升溫期,但烏蘇地區(qū)的年平均溫度均高于石河子地區(qū),升溫速度也大于石河子地區(qū);而石河子地區(qū)年降水量高于烏蘇地區(qū),相對烏蘇地區(qū)而言,石河子地區(qū)氣候條件較濕潤。

    圖3 烏蘇和石河子地區(qū)1954~2012年平均溫度和降水量的變化Fig.3 Variation of annual average temperature and total precipitation of Wusu and Shihezi in 1954-2012

    2.4 數據分析

    年輪寬度指數與氣候要素之間的相關關系和2個年表與氣候因子的滑動相關系數的計算,本研究分別選用了Dendroclim 2002軟件和SPSS軟件[29],旨在揭示不同地區(qū)樹木生長與氣候要素的關系,評估烏蘇和石河子地區(qū)雪嶺云杉生長節(jié)律的地區(qū)差異及其對氣候變暖的響應差異。同樣計算2個年表與月氣候因子之間的相關系數和滑動相關系數,從而得到不同地區(qū)樹木生長與氣候因子的關系,以及關系的穩(wěn)定性。氣候因子選自距采樣點最近的烏蘇、石河子氣象站自建站以來(1968~2012年)的前一年6月至當年10月的月平均氣溫、月降水量、月平均相對濕度和scPDSI指數。

    3 結果與分析

    3.1 年表的基本統計特征

    表2給出了兩研究區(qū)雪嶺云杉年表的統計參數,其中平均敏感度介于0.236~0.334,在年表所要求的0.15~0.80,所分析年表屬于敏感系列范圍[14]。標準差介于0.248~0.512,第一主成分的解釋量達到59.6%~62.9%,年表序列間、樹間和樹內相關系數達到0.533~0.613,表明2個采樣點的樣本均對群體的代表性較好即序列的共性較強,年表中的環(huán)境信息含量較大[30]。然而,2個采樣點的信噪比達到23.416~40.148,均大于3,EPS達到0.979~0.985,均超過了0.85的最低閾值[30],說明本次研究采樣較成功,年表中包含的氣候信息較多,樹木生長受氣候因子的影響較強,兩研究區(qū)的雪嶺云杉對氣候變化的響應敏感,從表2可知,烏蘇地區(qū)的年表各項統計參數均較高,其樹木生長對氣候變化的響應更敏感[31]。

    表2 天山中西部采樣點的差值年表(RES)的統計特征值

    注:R1.不同序列間的平均相關系數;R2.不同樹間的平均相關系數;R3.同一棵樹不同樣本之間的平均相關系數;M.S..平均敏感度;S.D..標準差;SNR.信噪比;EPS.樣本總體代表性;PCA1.第一主成分解釋量

    Note:R1.Inter-series mean correlation; R2.Mean correlations between trees; R3.Mean correlation within the same tree; M.S..Mean sensitivity; S.D..Standarad deviation; SNR.Signal to noise ratio; EPS.Expressed population signal; PCA1.The variance expressed by the first principal component

    3.2 雪嶺云杉徑向生長與氣候因子的相關關系

    為了建立可靠的年表,本研究中采納了樣本總體代表性EPS(Expressed population signal)[32]大于0.85的年份作為年表的起始點,這樣高于該起點的年表序列才是可靠的,烏蘇地區(qū)有效年表的起始年為1968年,而石河子地區(qū)有效年表的起始年為1941年,其2個采樣點的差值年表如圖4所示。

    圖4 研究區(qū)雪嶺云杉年輪寬度(RES)年表和樣本量實線代表輪寬指數;虛線為樣本量;粗豎線條是樣本總體代表性EPS>0.85的樣本量起點,表示年表在該起點之后是可靠的Fig.4 P.schrenkiana tree-ring width chronology(RES)and the amount of samples Solid line stands for the ring-width index,and dash line for numbers of core samples;Thick line for the start of expressed population signal(EPS>0.85),the chronology is reliable after the start.

    從圖4可以看出,烏蘇與石河子地區(qū)雪嶺云杉年表的樹輪指數變化比較大,說明其受氣候變化的影響比較大[33]。為了探討烏蘇和石河子地區(qū)樹木生長與氣候因子的關系,本文選擇了2個地區(qū)1968~2012年的月平均氣溫、月降水量、月平均相對濕度和scPDSI等氣候因子與采樣點的年表進行了相關分析,分析結果如圖5所示。

    圖5 樹輪度年表與各月氣象要素之間的相關統計 P6~P12.上一年4~12月;1~10.當年1~10月;RH.平均相對濕度;scPDSI.基于改進帕爾默干旱指數 *P<0.05Fig.5 The correlation coefficients between the tree-ring chronology and the climate factors P6~P12. Previous June-December; 1~10. Current year January-October; RH. Average relative humidity; scPDSI. Self-calibrating palmer drought severity index *P<0.05

    統計結果表明(圖5),烏蘇和石河子地區(qū)雪嶺云杉樹木生長對當地氣候因子的響應有一致的地方,但也存在響應差異。2個采樣點雪嶺云杉徑向生長均對生長季的氣候因子敏感,相對于石河子地區(qū)而言,烏蘇地區(qū)對水分和溫度更敏感。具體表現為:烏蘇地區(qū)雪嶺云杉徑向生長與上一年7、9月和當年5~7月平均氣溫呈顯著負相關,與當年5~8月降水量、上一年9月及當年5~8月平均相對濕度和上一年7~10月及當年3~10月scPDSI呈顯著正相關;石河子地區(qū)雪嶺云杉樹木生長與當年5~7月平均氣溫呈顯著負相關,與當年5~7月降水量、當年1、5、7月平均相對濕度和當年1~10月scPDSI呈顯著正相關??梢?,烏蘇和石河子地區(qū)雪嶺云杉森林下限樹木生長并不是主要受降水的影響,這兩個地區(qū)森林下限的雪嶺云杉徑向生長對水熱條件均呈現明顯的相關關系,樹木生長均受溫度、降水綜合的影響(圖5)。

    3.3 樹木生長與重要氣候因子的動態(tài)關系

    相關分析的結果顯示兩地區(qū)雪嶺云杉樹木生長對生長季5~7溫度和降水具有較高的敏感性。我們進一步對烏蘇和石河子地區(qū)雪嶺云杉徑向生長與生長季5~7月平均氣溫和降水量關系的時間動態(tài)進行了分析(圖6)。結果表明,自1969年隨著溫度升高,這兩個地區(qū)的雪嶺云杉徑向生長與氣候因子的關系均不穩(wěn)定。詳細地講,烏蘇地區(qū)雪嶺云杉年表與生長季5~7月平均氣溫和降水量的關系均逐漸增強,該年表與生長季溫度呈顯著負相關,與生長季降水呈顯著正相關,滑動相關系數都逐漸增大;與烏蘇地區(qū)不同,石河子地區(qū)雪嶺云杉年表與生長季5~7月平均氣溫的關系逐漸減弱,而與生長季5~7月降水量的關系逐漸增強,該年表與生長季溫度呈負相關,但該滑動相關系數逐漸減小,由顯著性變?yōu)椴伙@著性的氣候因子,與生長季降水呈顯著正相關,滑動相關系數也逐漸增大。這些分析結果表明在氣候變暖背景下烏蘇地區(qū)和石河子地區(qū)雪嶺云杉徑向生長對氣候因子的響應呈現不同,說明氣候條件的改變可以引起樹木生長與氣候因子關系的變化,以及徑向生長的區(qū)域差異。

    圖6 烏蘇和石河子地區(qū)年表序列與氣象數據的滑動相關分析,滑動窗口為20年Fig.6 Moving correlation between ring-width chronology and meteorological inWushu and Shihezi,the moving window is 20 years.

    3.4 增溫對雪嶺云杉徑向生長的影響

    年表和溫度序列的趨勢對比(圖7)顯示:烏蘇地區(qū)森林下限雪嶺云杉年表序列具有隨溫度升高而降低的趨勢,烏蘇地區(qū)雪嶺云杉年表中明顯地體現了溫度上升的“分離現象”,升溫導致樹木的水分脅迫加劇是出現分離現象的主要原因;但石河子地區(qū)森林下限的雪嶺云杉與溫度之間表現出一定的“吻合現象”,即隨著溫度的升高,該地區(qū)的雪嶺云杉徑向生長呈上升趨勢,這與石河子地區(qū)的氣候條件相對于烏蘇地區(qū)較濕潤有關,烏蘇地區(qū)年平均溫度均高于石河子地區(qū),烏蘇地區(qū)近50年來升溫速度大于2.2℃·50 a-1,石河子地區(qū)升溫速度大于1.2℃·50 a-1,而石河子近50年來年降水量均高于烏蘇地區(qū)(圖3)。

    圖7 1954~2012年烏蘇和石河子地區(qū)雪嶺云杉年表序列與年平均溫度一致性分析 黑實線代表年表序列;灰實線為年平均溫度;短劃線為相應時間序列的擬合趨勢線Fig.7 Consistency analysis between ring-width chronologies of Wusu and Shihezi and annual temperature in 1954-2012 Black solid line stands for the ring-width chronology,and the gray solid line for the annual average temperature.Dash line is the trend line.

    4 討論

    4.1 雪嶺云杉生長與氣候因子關系的區(qū)域差異

    目前普遍認為森林上限的樹木生長主要受溫度的影響而森林下限的樹木生長主要受降水的影響[19],而侯愛敏、彭劍鋒等[34~35]認為森林下限樹木生長受溫度、降水等氣候因子綜合影響,森林下限樹輪寬度與氣候因子有著復雜的相關關系,這種關系受氣候因子之間的相互制衡和因物種而異的樹木生長節(jié)律的影響。通過烏蘇和石河子地區(qū)的森林下限年表與氣候因子的相關分析,發(fā)現烏蘇和石河子地區(qū)樹木徑向生長限制因子并不是單一受降水的影響,而是受溫度、降水綜合影響的。但這兩個地區(qū)雪嶺云杉樹木生長對氣候因子的響應存在很多的異同點,可能與這兩個地區(qū)的生態(tài)環(huán)境的不同有關。詳細分析如下:

    共同之處表現為2個方面:①烏蘇和石河子地區(qū)雪嶺云杉徑向生長都與當年5~7月平均氣溫呈顯著負相關,與當年5~7月降水量、5、7月平均相對濕度和當年3~10月scPDSI呈顯著正相關,而雪嶺云杉生長期集中于5~8月,表明生長季的氣候因子對雪嶺云杉的徑向生長影響很大。從樹木生長的生理學角度來理解,森林下限土壤濕度對樹木生長起到關鍵性的控制作用,在干旱半干旱地區(qū),溫度和降水共同調控引起的土壤濕度變化在樹木徑向生長中起著很重要的作用[29]。烏蘇和石河子地區(qū)5~7月的降水量相對樹木生長需要來說還不充足,生長季5~7月的高溫加重了土壤水分的匱乏[16]。土壤干旱脅迫會降低樹木葉水勢,使葉肉細胞對水分的吸附力增大,限制了葉內水分的運動,從而使樹木的光合作用速率減弱,影響到樹木的生長[36]。此外,雪嶺云杉是喜濕樹種,分布在土壤濕度和空氣濕度較大的陰溝與陰坡,自然更新與生長發(fā)育均較好[19],雪嶺云杉對土壤濕度非常敏感。②烏蘇和石河子地區(qū)雪嶺云杉徑向生長與scPDSI指數的相關性均為正相關(圖5),這一結果與天山北坡中部[18]、內蒙古東部[37]、賀蘭山北部[38]樹木徑向生長對PDSI的相關相同,scPDSI指數使用彭曼公式計算蒸散發(fā),更適合代表干旱半干旱地區(qū)的自然干旱程度,更好地指示土壤水分變化,表征土壤對樹木生長的可供水量[39]。

    不同之處表現為:烏蘇地區(qū)雪嶺云杉徑向生長對溫度和水分的響應更敏感,除了2個采樣點均對生長季5~7的氣候因子存在顯著關系外,烏蘇地區(qū)樹木生長與上一年7、9月平均氣溫呈顯著負相關,與當年8月降水、上一年9月及當年8月平均相對濕度和上一年7~10月的scPDSI呈顯著正相關。石河子地區(qū)樹木生長與當年1月平均濕度和1~2月的scPDSI呈顯著正相關。烏蘇地區(qū),溫度高、降水量少使土壤蒸發(fā)量增大,其scPDSI較低,說明土壤比石河子地區(qū)干旱。溫度對烏蘇地區(qū)樹木生長的影響具有明顯的“滯后效應”,上一年7、9月氣溫過高,將會導致該地區(qū)生長季前的降水大量蒸發(fā),減少土壤含水量的積累,易形成水分脅迫,不利于樹體在生長季后期營養(yǎng)物質的積累,從而影響來年生長季初期的營養(yǎng)物質的供給,繼而影響樹木生長[31]。7~10月是雪嶺云杉樹木生長季中末期,此時濕潤的條件有利于積累更多的光合產物及木質化的過程,繼而影響影響來年樹木生長的營養(yǎng)物質的供給及當年樹木生長[35]。正如結果顯示烏蘇地區(qū)樹木生長與當年8月降水量、平均相對濕度和上一年7~10月scPDSI呈顯著正相關。與烏蘇地區(qū)不同,石河子地區(qū)溫度較低,降水量較多,平均相對濕度和scPDSI高,說明土壤比烏蘇地區(qū)濕潤。1~2月是雪嶺云杉樹木生長季前期,是其早材形成的關鍵時期,該時段濕潤的條件對快速生長的早材有促進作用[40],故當年1月平均相對濕度和1~2月的scPDSI均與樹木生長呈顯著正相關。

    4.2 氣候變化對雪嶺云杉生長的影響

    分布于不同地區(qū)森林下限的雪嶺云杉生境條件隨著氣候變化而發(fā)生變化,即樹木分布范圍內的水熱條件的增加,生境條件的變化對樹木徑向生長有顯著的影響作用[13]?;瑒酉嚓P(圖6)結果表明烏蘇和石河子地區(qū)森林下限雪嶺云杉徑向生長與生長季氣候因子的關系并不穩(wěn)定,而是隨著時間推移表現出一定的變化,這一滑動相關的結果與蘆芽山華北落葉松、長白山紅松和魚鱗云杉的滑動相關研究結果一致[6,13,15]?;瑒酉嚓P(圖6)的結果顯示烏蘇和石河子地區(qū)生長季的水熱條件的增加,使得氣候因子對樹木生長的限制作用也產生了明顯的變化。值得注意的是,隨著溫度升高,烏蘇地區(qū)年表與生長季溫度呈顯著負相關,與生長季降水呈顯著正相關,相關系數都逐漸增大(圖6),樹木生長受水分限制作用逐漸增強。隨著氣溫升高,植物體的蒸騰作用以及土壤水分的蒸散都會加快,加快了水分散失和養(yǎng)分消耗,從而影響樹木徑向生長[13],升溫造成的水分缺失是烏蘇地區(qū)樹木生長與溫度上升趨勢相背離的最主要原因。石河子地區(qū)年表與生長季溫度呈負相關,但該相關系數逐漸減小,樹木生長與溫度的響應由顯著逐漸變?yōu)椴伙@著;與生長季降水呈顯著正相關,相關系數也逐漸增大(圖6),而氣候變化促進了該地區(qū)樹木生長。

    烏蘇和石河子地區(qū)溫度、降水量增加引起了當地雪嶺云杉徑向生長與氣候因子關系的變化,進一步引起了生長節(jié)律的區(qū)域差異。烏蘇地區(qū)近50年溫度上升幅度顯著,升溫速度大于2.5℃·50 a-1,而該地區(qū)雪嶺云杉年輪指數呈下降趨勢,表明溫度上升阻礙了該地區(qū)雪嶺云杉的徑向生長。烏蘇地區(qū)樹木生長與生長季溫度顯著負相關以及受水分限制逐漸增強結果與長白山的魚鱗云杉、岳樺和落葉松[5,13,41]的研究結果一致。隨著溫度過高,將會導致該地區(qū)的降水大量蒸發(fā),減少土壤含水量的積累,而同期的降水比較豐富,則可以起到補償作用[30];且烏蘇地區(qū)降水量增加不顯著(圖3),溫度升高土壤含水量加劇蒸發(fā),降水量增加的補償作用不足,樹木呼吸作用和代謝作用增強,生長所需的養(yǎng)料開始增加,從而導致樹木生長緩慢[14];而石河子地區(qū)升溫速度大于1.2℃·50 a-1,石河子地區(qū)降水量明顯高于烏蘇地區(qū),而溫度明顯低于烏蘇地區(qū)(圖3),該地區(qū)雪嶺云杉年輪指數呈上升趨勢,桑衛(wèi)國等[31]也發(fā)現,隨著溫度的升高,天山昭蘇地區(qū)和小渠子地區(qū)也呈上升趨勢。石河子地區(qū)降水量增加對溫度升高起到了補償作用,石河子地區(qū)雪嶺云杉的生長受溫度的限制作用由顯著轉變?yōu)椴伙@著,隨著氣候變化,徑向生長與溫度的關系減弱。石河子地區(qū)樹木生長對溫度的敏感度降低的現象與北半球高緯度高海拔地區(qū)、阿拉斯加地區(qū)[7~8]樹木生長的研究結果一致。結果表明,溫度和降水量的改變引起了烏蘇和石河子地區(qū)森林下限雪嶺云杉樹木生長與氣候因子關系的不同變化。

    5 結論

    不同地區(qū)的雪嶺云杉徑向生長對氣候因子的響應具有明顯差異。烏蘇地區(qū)樹木生長與上一年7、9月平均氣溫呈顯著負相關,與當年8月降水、上一年9月及當年8月平均相對濕度和上一年7~10月的scPDSI呈顯著正相關;而石河子地區(qū)樹木生長與當年1月平均濕度和1~2月的scPDSI呈顯著正相關。與石河子地區(qū)對比分析發(fā)現,烏蘇地區(qū)雪嶺云杉徑向生長受水熱條件的限制更強,該響應差異的存在與兩個地區(qū)生態(tài)環(huán)境不同有關。

    本研究采用樹木氣候學方法分析了氣候變暖對烏蘇和石河子地區(qū)森林下限雪嶺云杉徑向生長的影響,結果表明溫度升高抑制了烏蘇地區(qū)雪嶺云杉的生長,而促進了石河子地區(qū)雪嶺云杉的生長,氣候條件的不同引起了烏蘇和石河子地區(qū)樹木生長與氣候因子的區(qū)域差異和雪嶺云杉徑向生長節(jié)律的區(qū)域差異。氣候變暖造成了烏蘇地區(qū)雪嶺云杉年表序列對溫度和降水變化敏感性上升,造成了石河子地區(qū)雪嶺云杉年表了序列對降水變化敏感性上升,對溫度變化敏感性下降。

    1.邵雪梅.樹輪年代學的若干進展[J].第四紀研究,1997,17(3):265-271.

    Shao X M.Advancements in dendrochronology[J].Quaternary Sciences,1997,17(3):265-271.

    2.楊青霄,朱良軍,王曉春.涼水自然保護區(qū)紅松樹輪年表建立及特征年分析[J].植物研究,2015,35(3):418-424.

    Yang Q X,Zhu L J,Wang X C.Development ofPinuskoraiensisTree-ring chronology and master year analysis in Liangshui National Natural Reserve,China[J].Bulletin of Botanical Research,2015,35(3):418-424.

    3.范敏杰,袁玉江,魏文壽,等.用樹木年輪重建伊犁南天山北坡西部的降水量序列[J].干旱區(qū)地理,2007,30(2):268-273.

    Fan M J,Yuan Y J,WEI W S,et al.Reconstruction of precipitation series from tree-rings at the northern slopes of south Tianshan Mountains of Yili prefecture[J].Arid Land Geography,2007,30(2):268-273.

    4.秦莉,袁玉江,喻樹龍,等.賽里木湖流域雪嶺云杉(Piceaschrenkiana)樹木徑向生長對氣候變化的響應[J].中國沙漠,2015,35(1):113-119.

    Qin L,Yuan Y J,Yu S L,et al.Response of tree-ring growth ofPiceaschrenkianato climate change in the Sayram Lake Basin,Xinjiang,China[J].Journal of Desert Research,2015,35(1):113-119.

    5.于大炮,王順忠,唐立娜,等.長白山北坡落葉松年輪年表及其與氣候變化的關系[J].應用生態(tài)學報,2005,16(1):14-20.

    Yu D P,Wang S Z,Tang L N,et al.Relationship between tree-ring chronology ofLarixolgensisin Changbai Mountains and the climate change[J].Chinese Journal of Applied Ecology,2005,16(1):14-20.

    6.張文濤,江源,王明昌,等.蘆芽山陽坡不同海拔白杄徑向生長對氣候變暖的響應[J].植物生態(tài)學報,2013,37(12):1142-1152.

    Zhang W T,Jiang Y,Wang M C,et al.Responses of radial growth to climate warming inPiceameyeritrees growing at different elevations on the southern slope of Luya Mountain[J].Chinese Journal of Plant Ecology,2013,37(12):1142-1152.

    7.Briffa K R,Schweingruber F H,Jones P D,et al.Reduced sensitivity of recent tree-growth to temperature at high northern latitudes[J].Nature,1998,391(6668):678-682.

    8.D’arrigo R,Mashig E,Frank D,et al.Temperature variability over the past millennium inferred from Northwestern Alaska tree rings[J].Climate Dynamics,2005,24(2-3):227-236.

    9.Carrer M,Urbinati C.Long-term change in the sensitivity of tree-ring growth to climate forcing inLarixdecidua[J].New Phytologist,2006,170(4):861-872.

    10.于健,徐倩倩,何秀,等.長白山東坡落葉松樹輪寬度對氣候響應的分離效應[J].中南林業(yè)科技大學學報,2013,33(3):89-97.

    Yu J,Xu Q Q,He X,et al.Response divergence ofLarixolgensistree-ring widths to climate variation in eastern slope of Changbai Mountain,northeast China[J].Journal of Central South University of Forestry & Technology,2013,33(3):89-97.

    11.王淼,白淑菊,陶大立,等.大氣增溫對長白山林木直徑生長的影響[J].應用生態(tài)學報,1995,6(2):128-132.

    Wang M,Bai S J,Tao D L,et al.Effect of rise in air-temperature on tree ring growth of forest on Changbai Mountain[J].Chinese Journal of Applied Ecology,1995,6(2):128-132.

    12.Leonelli G,Pelfini M,Battipaglia G,et al.Site-aspect influence on climate sensitivity over time of a high-altitudePinuscembratree-ring network[J].Climatic Change,2009,96(1-2):185-201.

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

    Li G Q,Bai F,Sang W G.Different responses of radial growth to climate warming inPinuskoraiensisandPiceajezoensisvar.Komaroviiat their upper elevational limits in Changbai,China[J].Chinese Journal of Plant Ecology,2011,35(5):500-511.

    14.郭濱德,張遠東,王曉春.川西高原不同坡向云、冷杉樹輪對快速升溫的響應差異[J].應用生態(tài)學報,2016,27(2):354-364.

    Guo B D,Zhang Y D,Wang X C.Response of Picea purpurea andAbiesfaxonianatree rings at different slope aspects to rapid warming in western Sichuan,China[J].Chinese Journal of Applied Ecology,2016,27(2):354-364.

    15.張文濤,江源,王明昌,等.蘆芽山陽坡不同海拔華北落葉松徑向生長對氣候變化的響應[J].生態(tài)學報,2010,35(19):6481-6488.

    Zhang W T,Jiang Y,Wang M C,et al.Responses of radial growth inLarixprincipis-rupprechtiito climate change along an elevation gradient on the southern slope of Luya Mountain[J].Acta Ecologica Sinica,2010,35(19):6481-6488.

    16.劉敏,毛子軍,厲悅,等.不同緯度闊葉紅松林紅松徑向生長對氣候因子的響應[J].應用生態(tài)學報,2016,27(5):1341-1352.

    Liu M,Mao Z J,Li Y,et al.Response of radial growth ofPinuskoraiensisin broad-leaved Korean pine forests with different latitudes to climatical factors[J].Chinese Journal of Applied Ecology,2016,27(5):1341-1352.

    17.郭允允,劉鴻雁,任佶,等.天山中段樹木生長對氣候垂直梯度的響應[J].第四紀研究,2007,27(3):322-331.

    Guo Y Y,Liu H Y,Ren J,et al.Responses of tree growth to vertical climate gradient in the middle section of the Tianshan mountains[J].Quaternary Sciences,2007,27(3):322-331.

    18.喻樹龍,袁玉江,秦莉,等.天山北坡中部不同海拔高度雪嶺云杉樹輪寬度氣候響應對比分析[J].沙漠與綠洲氣象,2016,10(3):30-38.

    Yu S L,Yuan Y J,Qin L,et al.Tree-ring-width growth responses ofPiceaschrenkianato climate change for different elevations in the central Tianshan Mountains[J].Desert and Oasis Meteorology,2016,10(3):30-38.

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

    Wang T.Ecological study onPiceaschrenkianaforest along an altitudinal gradient in the central Tianshan Mountains[D].Wuhan:Wuhan University,2004.

    20.吳祥定.樹木年輪與氣候變化[M].北京:氣象出版社,1990:83-85.

    Wu X D.Tree rings and climate change[M].Beijing:China Meteorological Press,1990:83-85.

    21.Stokes M A,Smiley T L.An introduction to tree-ring dating[M].Tucson,Ariz,USA:The University of Arizona Press,1996.

    22.Holmes R L,Adams R K,Fritts H C.Tree-Ring chronologies of western north America:California,eastern Oregon and Northern Great Basin,with procedures used in the chronology development work,including users manuals for computer programs COFECHA and ARSTAN.Chronology Series Ⅵ[R].Tucson:Laboratory of Tree-Ring Research,University of Arizona,1986.

    23.Holmes R L.Computer-assisted quality control in tree-ring dating and measurement[J].Tree-Ring Bulletin,1983,44(3):69-75.

    24.Cook E R,Holmes R L.Program ARSTAN user manual[R].Tucson,Arizona:Laboratory of Tree-Ring Research,University of Arizona,1984.

    25.邵雪梅,吳祥定.華山樹木年輪年表的建立[J].地理學報,1994,49(2):174-181.

    Shao X M,Wu X D.Tree-ring chronologies forPinusarmandifranch from Huashan,China[J].Acta Geographica Sinica,1994,49(2):174-181.

    26.Mann H B.Nonparametric tests against trend[J].Econometrica,1945,13(3):245-259.

    27.Kohler M A.On the use of double-mass analysis for testing the consistency of meteorological records and for making required adjustments[J].Bulletin of American Meteorological Society,1949,30:188-189.

    28.Van Der Schrier G,Barichivich J,Briffa K R,et al.A scPDSI-based global data set of dry and wet spells for 1901-2009[J].Journal of Geophysical Research,2013,118(10):4025-4048.

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

    30.Wigley T W L,Briffa K R,Jones P D.On the average value of correlated time series,with applications in dendroclimatology and hydrometeorology[J].Journal of Climate and Applied Meteorology,1984,23(2):201-213.

    31.桑衛(wèi)國,王云霞,蘇宏新,等.天山云杉樹輪寬度對梯度水分因子的響應[J].科學通報,2007,52(19):2292-2298.

    Sang W G,Wang Y X,Su H X,et al.Response of tree-ring width to rainfall gradient along the Tianshan Mountains of northwestern China[J].Chinese Science Bulletin,2007,52(21):2954-2962.

    32.張同文,袁玉江,魏文壽,等.開都河流域天山樺樹輪寬度年表的建立及其氣候響應[J].生態(tài)學報,2015,35(9):3034-3042.

    Zhang T W,Yuan Y J,Wei W S,et al.Development of the tree-ring width chronology ofBetulatianschanicafor the Kaiduhe River Watershed and its climate response[J].Acta Ecologica Sinica,2015,35(9):3034-3042.

    33.伊建平,張海州,李留振,等.天山中部雪嶺云杉樹木年輪寬度與氣候因子關系[J].林業(yè)建設,2008,(2):46-49.

    Yi J P,Zhang H Z,Li L Z,et al.Relationship between tree-ring width ofPiceaschrenkianaand climatic factors in the middle area of Tianshan mountain[J].Forestry Construction,2008,(2):46-49.

    34.侯愛敏,彭少麟,周國逸.樹木年輪對氣候變化的響應研究及其應用[J].生態(tài)科學,1999,18(3):16-23.

    Hou A M,Peng S L,Zhou G Y.The study of the reactions of tree rings to the climate change and its applications[J].Ecologic Science,1999,18(3):16-23.

    35.彭劍鋒,勾曉華,陳發(fā)虎,等.天山東部西伯利亞落葉松樹輪生長對氣候要素的響應分析[J].生態(tài)學報,2006,26(8):2723-2731.

    Peng J F,Gou X H,Chen F H,et al.The responses of growth ring width variations ofLarixsibiricaLedb to climatic change in eastern Tianshan Mountains[J].Acta Ecologica Sinica,2006,26(8):2723-2731.

    36.郭建平,高素華.CO2濃度和輻射強度變化對沙柳光合作用速率影響的模擬研究[J].生態(tài)學報,2004,24(2):181-185.

    Guo J P,Gao S H.Growth chamber study on the impact of CO2concentration and radiant intensity onSalixpsammophilaphotosynthetic rate[J].Acta Ecologica Sinica,2004,24(2):181-185.

    37.梁爾源,邵雪梅,劉鴻雁,等.樹輪所記錄的公元1842年以來內蒙古東部渾善達克沙地PDSI的變化[J].科學通報,2007,52(14):1694-1699.

    Liang E Y,Shao X M,Liu H Y,et al.Tree-ring based PDSI reconstruction since AD 1842 in the Ortindag Sand Land,east Inner Mongolia[J].Chinese Science Bulletin,2007,52(19):2715-2721.

    38.陳峰,袁玉江,魏文壽,等.樹輪記錄的賀蘭山北部5~7月PDSI變化[J].氣候變化研究進展,2010,6(5):344-348.

    Chen F,Yuan Y J,Wei W S,et al.Long-term drought severity variations recorded in tree rings in the northern Helan Mountains[J].Advances in Climate Change Research,2010,6(5):344-348.

    39.王兆禮,李軍,黃澤勤,等.基于改進帕默爾干旱指數的中國氣象干旱時空演變分析[J].農業(yè)工程學報,2016,32(2):161-168.

    Wang Z L,Li J,Huang Z Q,et al.Spatiotemporal variations analysis of meteorological drought in China based on scPDSI[J].Transactions of the Chinese Society of Agricultural Engineering,2016,32(2):161-168.

    40.陳峰,袁玉江,魏文壽,等.樹輪記錄的伊犁地區(qū)近354年帕爾默干旱指數變化[J].高原氣象,2011,30(2):355-362.

    Chen F,Yuan Y J,Wei W S,et al.Variations of long-term palmer drought index in recent 354 years in Yili based on tree-ring record[J].Plateau Meteorology,2011,30(2):355-362.

    41.Yu D P,Gu H Y,Wang J D,et al.Relationships of climate change and tree ring ofBetulaermaniitree line forest in Changbai Mountain[J].Journal of Forestry Research,2005,16(3),187-192.

    “One thousand plan” from Chinese academy of sciences(Y474161;Y672141);“Young scholars of western” from Chinese academy of sciences(2015-XBQN-B-22);International(regional) cooperation and the exchange program(41361140361)

    introduction:ZHANG Yan-Jing(1991—),female,postgraduate,engaged in the research of dendroclimatology.

    date:2016-12-20

    ResponseDifferencesofRadialGrowthofPiceaschrenkianatoClimateWarminginMidwesternTianshanMountains

    ZHANG Yan-Jing1,2ZHENG Hong-Wei1*YU Rui-De1YANG Mei-Lin1SHI Bing-Bing1,2WANG Lei3

    (1.Laboratory of Environment Change in Arid Lands,Xinjiang Institute of Ecology and Geography Chinese Academy of Sciences,Urumqi 830011;2.University of Academy of Sciences,Beijing 100049;3.Insititute of Modern Forestry,Xinjiang Academy of Forestry Science,Urumqi 830011)

    The experiment was conducted to study the response differences of radial growth ofPiceaschrenkianato climate warming in Wusu and Shihezi regions. Dendrochronological methods including response function and evolutionary intervals analysis method were applied to study the relationships between radial growth and climatic factors ofP.schrenkianaat the lower forest border of the two regions. The radial growth ofP.schrenkianaat two regions was significantly correlated to the local climatic factors in the growing season, but the differences in response were statistically significant. In Wusu, the radial growth ofP.schrenkianawas significantly negatively correlated to the average temperature in previous July and September but significantly positively correlated to precipitation in current August, the average relative humidity in previous September and current August, and scPDSI from July to October last year. In Shihezi, the radial growth ofP.schrenkianawas significantly positively correlated to the average relative humidity in current January, and scPDSI in current January and February. The tree ring width index in Wusu showed a decreasing trend namely “separation phenomenon” with the increasing of temperature, while that in Shihezi was increased with the increasing temperature. ForP.schrenkianain Wusu, water stress caused by temperature rising was the main factor causing the opposite tendency between radial growth and temperature change. The increase of temperature and precipitation in the growing season mutually enhanced theP.schrenkianagrowth in Shihezi. The sensitivity of radial growth ofP.schrenkianato average temperature and precipitation from May to July was increasing in Wusu. The sensitivity of radial growth ofP.schrenkianato precipitation from May to July was increased in Shihezi, while the sensitivity to average temperature was decreased from May to July. Therefore, the regional differences in climate change were the main factor in the radial growth difference at two regions and the sensitivity change to temperature and precipitation in chronological series.

    tree rings;radial growth;climate warming;sensitivity;Tianshan Mountains;response difference

    中國科學院千人計劃科研項目(Y474161;Y672141);中國科學院“西部青年學者”項目(2015-XBQN-B-22);國際(地區(qū))合作與交流項目(41361140361)

    張艷靜(1991—),女,碩士研究生,主要從事樹木年輪氣候學研究。

    * 通信作者:E-mail:hzheng@ms.xjb.ac.cn

    2016-12-20

    * Corresponding author:E-mail:hzheng@ms.xjb.ac.cn

    Q948

    A

    10.7525/j.issn.1673-5102.2017.03.004

    猜你喜歡
    烏蘇年表氣候因子
    青海省果洛州高原鼠兔宜生區(qū)分布及其對氣候變化的響應
    草地學報(2024年6期)2024-01-01 00:00:00
    年表
    世界建筑(2022年10期)2022-10-25 08:26:48
    王錫良年表
    石魯年表
    藝術品鑒(2019年12期)2020-01-18 08:47:30
    妙手仁心的白衣天使
    唇槍舌劍下的赤子之心
    四川水稻產量及其構成要素對不同生育期氣候因子的響應分析
    天山北坡NDVI對氣候因子響應的敏感性分析
    沙漠地區(qū)微波地表發(fā)射率年內變化規(guī)律與氣候因子的關系分析
    梁玉繩年表
    天一閣文叢(2013年1期)2013-11-06 10:01:20
    9191精品国产免费久久| 亚洲美女黄色视频免费看| 亚洲成人国产一区在线观看 | 午夜老司机福利片| 90打野战视频偷拍视频| 国产亚洲av高清不卡| 十八禁网站网址无遮挡| 国产老妇伦熟女老妇高清| 午夜激情久久久久久久| 欧美日韩成人在线一区二区| 91精品三级在线观看| 免费少妇av软件| 女人被躁到高潮嗷嗷叫费观| 老司机在亚洲福利影院| 成人漫画全彩无遮挡| 国产精品成人在线| 成人亚洲精品一区在线观看| 巨乳人妻的诱惑在线观看| 免费在线观看视频国产中文字幕亚洲 | 精品酒店卫生间| 亚洲欧美色中文字幕在线| 久久精品国产综合久久久| 欧美久久黑人一区二区| 91aial.com中文字幕在线观看| 在线天堂中文资源库| 国产精品久久久久久久久免| 日日爽夜夜爽网站| 日韩精品有码人妻一区| 国产男女超爽视频在线观看| 亚洲精品久久成人aⅴ小说| 久热爱精品视频在线9| 日日爽夜夜爽网站| 欧美黄色片欧美黄色片| 日韩免费高清中文字幕av| 国产一区二区 视频在线| 女人精品久久久久毛片| 五月开心婷婷网| 一二三四在线观看免费中文在| 超色免费av| 在线天堂中文资源库| 丁香六月天网| 亚洲欧美清纯卡通| 亚洲精品国产一区二区精华液| 欧美精品一区二区免费开放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲第一青青草原| 搡老岳熟女国产| 老司机在亚洲福利影院| 嫩草影院入口| 日本猛色少妇xxxxx猛交久久| 色吧在线观看| 香蕉国产在线看| 午夜福利网站1000一区二区三区| 久久久久精品人妻al黑| 精品亚洲成国产av| 国产爽快片一区二区三区| 亚洲精品中文字幕在线视频| 久久久久精品性色| av电影中文网址| 我的亚洲天堂| av卡一久久| 看免费av毛片| 国产精品人妻久久久影院| 亚洲熟女精品中文字幕| 交换朋友夫妻互换小说| 90打野战视频偷拍视频| 亚洲成人av在线免费| 老司机影院毛片| 看免费成人av毛片| 欧美日本中文国产一区发布| 亚洲第一青青草原| 亚洲少妇的诱惑av| 精品亚洲乱码少妇综合久久| 欧美日韩亚洲高清精品| 国产精品国产三级国产专区5o| 国产片特级美女逼逼视频| 色播在线永久视频| 我的亚洲天堂| 色视频在线一区二区三区| 久久影院123| 久久久精品94久久精品| 爱豆传媒免费全集在线观看| 各种免费的搞黄视频| 各种免费的搞黄视频| 十分钟在线观看高清视频www| 亚洲av成人精品一二三区| 午夜福利网站1000一区二区三区| 叶爱在线成人免费视频播放| 黄片小视频在线播放| 18禁国产床啪视频网站| 飞空精品影院首页| 亚洲精品一区蜜桃| 亚洲精品一二三| 欧美黑人欧美精品刺激| 免费在线观看完整版高清| 欧美精品高潮呻吟av久久| 免费少妇av软件| 一区福利在线观看| 欧美在线黄色| 九九爱精品视频在线观看| 丝瓜视频免费看黄片| 国产精品蜜桃在线观看| 久久久久久久国产电影| a 毛片基地| 超碰成人久久| 欧美人与性动交α欧美软件| 女性被躁到高潮视频| 婷婷色综合www| av国产久精品久网站免费入址| 亚洲精品久久午夜乱码| 久久国产亚洲av麻豆专区| 多毛熟女@视频| 成年女人毛片免费观看观看9 | 国产亚洲最大av| bbb黄色大片| 国产福利在线免费观看视频| 国产极品天堂在线| 国产欧美日韩综合在线一区二区| 国产av国产精品国产| 人人妻人人添人人爽欧美一区卜| 亚洲 欧美一区二区三区| 飞空精品影院首页| 日韩大片免费观看网站| 大码成人一级视频| 成人18禁高潮啪啪吃奶动态图| 成人国语在线视频| 久久久精品国产亚洲av高清涩受| 日韩大片免费观看网站| 高清av免费在线| 日本vs欧美在线观看视频| bbb黄色大片| 精品国产乱码久久久久久小说| av免费观看日本| 国产 一区精品| 大香蕉久久成人网| 国产黄色免费在线视频| 如日韩欧美国产精品一区二区三区| 哪个播放器可以免费观看大片| 精品国产乱码久久久久久小说| 性高湖久久久久久久久免费观看| 考比视频在线观看| 在线观看免费视频网站a站| 天天躁狠狠躁夜夜躁狠狠躁| 久久国产精品男人的天堂亚洲| 精品人妻在线不人妻| 亚洲国产日韩一区二区| 1024香蕉在线观看| 久久久久精品人妻al黑| 啦啦啦 在线观看视频| 亚洲中文av在线| 精品视频人人做人人爽| h视频一区二区三区| 亚洲图色成人| 激情视频va一区二区三区| 国产乱来视频区| 在线免费观看不下载黄p国产| 建设人人有责人人尽责人人享有的| 在线 av 中文字幕| 免费日韩欧美在线观看| 嫩草影院入口| 日日摸夜夜添夜夜爱| 国产成人av激情在线播放| 丁香六月欧美| 中文字幕精品免费在线观看视频| 老司机深夜福利视频在线观看 | 国产成人欧美| 天堂俺去俺来也www色官网| 久久人人爽人人片av| 黄色视频在线播放观看不卡| 亚洲三区欧美一区| 中文字幕亚洲精品专区| 男女无遮挡免费网站观看| 久久鲁丝午夜福利片| 51午夜福利影视在线观看| 天堂中文最新版在线下载| 纵有疾风起免费观看全集完整版| 亚洲熟女精品中文字幕| 99久久精品国产亚洲精品| 一本—道久久a久久精品蜜桃钙片| 欧美精品av麻豆av| 热re99久久国产66热| 国产精品熟女久久久久浪| 九色亚洲精品在线播放| 亚洲精品国产区一区二| 91成人精品电影| 日韩一区二区三区影片| 建设人人有责人人尽责人人享有的| 自线自在国产av| 精品一品国产午夜福利视频| 亚洲久久久国产精品| 欧美精品一区二区大全| 中文欧美无线码| 亚洲国产av影院在线观看| 欧美国产精品一级二级三级| 久久人人爽av亚洲精品天堂| 亚洲自偷自拍图片 自拍| 青春草国产在线视频| 嫩草影视91久久| 亚洲精品久久成人aⅴ小说| 90打野战视频偷拍视频| 亚洲激情五月婷婷啪啪| 久久性视频一级片| 麻豆av在线久日| 日韩av不卡免费在线播放| 丝袜美腿诱惑在线| 久久久久久久久久久久大奶| 国产熟女午夜一区二区三区| 精品视频人人做人人爽| 精品国产一区二区三区四区第35| 久久亚洲国产成人精品v| 天天添夜夜摸| 亚洲激情五月婷婷啪啪| 亚洲欧洲国产日韩| 大码成人一级视频| 日本av手机在线免费观看| 精品国产国语对白av| 97精品久久久久久久久久精品| 欧美日韩精品网址| 纯流量卡能插随身wifi吗| 国产一区有黄有色的免费视频| 欧美激情极品国产一区二区三区| 国产精品二区激情视频| av.在线天堂| 免费在线观看黄色视频的| 久久99热这里只频精品6学生| 亚洲 欧美一区二区三区| 亚洲欧美中文字幕日韩二区| 狠狠精品人妻久久久久久综合| 日韩精品免费视频一区二区三区| 国产一区二区三区综合在线观看| 免费观看人在逋| 国产精品成人在线| 日韩一卡2卡3卡4卡2021年| 黑人欧美特级aaaaaa片| 久久久久网色| 一个人免费看片子| 你懂的网址亚洲精品在线观看| 日本猛色少妇xxxxx猛交久久| 久久综合国产亚洲精品| 亚洲欧美清纯卡通| 精品一区二区三卡| 国产黄频视频在线观看| 9色porny在线观看| 性色av一级| 国产午夜精品一二区理论片| 国产精品 国内视频| 一级毛片电影观看| 日韩中文字幕欧美一区二区 | 亚洲四区av| 欧美 日韩 精品 国产| 久久精品熟女亚洲av麻豆精品| 黄色视频不卡| 制服诱惑二区| 国产精品欧美亚洲77777| 国产精品女同一区二区软件| 女性被躁到高潮视频| av线在线观看网站| 亚洲国产日韩一区二区| 满18在线观看网站| 日韩制服骚丝袜av| 青青草视频在线视频观看| 不卡av一区二区三区| 午夜福利视频在线观看免费| 国产精品久久久人人做人人爽| 波野结衣二区三区在线| 亚洲精品久久久久久婷婷小说| 一本色道久久久久久精品综合| 久久精品国产综合久久久| 美女视频免费永久观看网站| 国产精品久久久人人做人人爽| 免费女性裸体啪啪无遮挡网站| 欧美日韩成人在线一区二区| 亚洲av成人精品一二三区| 黑人猛操日本美女一级片| netflix在线观看网站| 99热全是精品| 精品人妻一区二区三区麻豆| 国产探花极品一区二区| 天天躁夜夜躁狠狠久久av| 亚洲欧美成人综合另类久久久| 日韩成人av中文字幕在线观看| 水蜜桃什么品种好| 午夜影院在线不卡| 你懂的网址亚洲精品在线观看| 日本av免费视频播放| 国产精品三级大全| 五月天丁香电影| 亚洲欧美色中文字幕在线| 美女中出高潮动态图| 国产免费现黄频在线看| 欧美日韩福利视频一区二区| 两个人看的免费小视频| 日本欧美国产在线视频| av线在线观看网站| 久久鲁丝午夜福利片| 亚洲精品久久久久久婷婷小说| 大片电影免费在线观看免费| 国产精品麻豆人妻色哟哟久久| 国产成人欧美| 国产探花极品一区二区| 亚洲精品,欧美精品| 成人影院久久| 黄色视频在线播放观看不卡| 搡老岳熟女国产| 国产精品av久久久久免费| 婷婷色综合大香蕉| 成人漫画全彩无遮挡| 久久久久国产精品人妻一区二区| 国产又爽黄色视频| 五月开心婷婷网| 哪个播放器可以免费观看大片| 夫妻午夜视频| 一本久久精品| 成年av动漫网址| 亚洲久久久国产精品| 人成视频在线观看免费观看| 啦啦啦在线观看免费高清www| 最近中文字幕高清免费大全6| 欧美 亚洲 国产 日韩一| 亚洲av成人不卡在线观看播放网 | 久久国产精品男人的天堂亚洲| 免费观看人在逋| 午夜老司机福利片| 久久精品久久精品一区二区三区| 婷婷色av中文字幕| 亚洲,欧美精品.| 欧美亚洲 丝袜 人妻 在线| 伦理电影免费视频| 免费黄色在线免费观看| 两个人免费观看高清视频| 纵有疾风起免费观看全集完整版| 色综合欧美亚洲国产小说| 国产精品嫩草影院av在线观看| 天天影视国产精品| 狠狠婷婷综合久久久久久88av| 亚洲综合色网址| 欧美乱码精品一区二区三区| 亚洲伊人色综图| 美女脱内裤让男人舔精品视频| 最新在线观看一区二区三区 | 一区在线观看完整版| 夜夜骑夜夜射夜夜干| 久久精品aⅴ一区二区三区四区| 亚洲精品日韩在线中文字幕| 男女国产视频网站| 国产精品国产三级国产专区5o| 久久毛片免费看一区二区三区| 国产精品偷伦视频观看了| 国产精品一区二区在线不卡| 国产无遮挡羞羞视频在线观看| a级片在线免费高清观看视频| 视频在线观看一区二区三区| 欧美久久黑人一区二区| 日韩中文字幕欧美一区二区 | 一区二区日韩欧美中文字幕| 大话2 男鬼变身卡| 欧美中文综合在线视频| 亚洲四区av| 国产精品香港三级国产av潘金莲 | 日韩熟女老妇一区二区性免费视频| 国产熟女午夜一区二区三区| 人人澡人人妻人| 亚洲国产成人一精品久久久| 亚洲av电影在线进入| av福利片在线| 国产成人av激情在线播放| 另类精品久久| 五月开心婷婷网| 日韩制服骚丝袜av| netflix在线观看网站| 亚洲少妇的诱惑av| 久久人妻熟女aⅴ| 国产亚洲午夜精品一区二区久久| 欧美日韩精品网址| 久久精品aⅴ一区二区三区四区| 国产女主播在线喷水免费视频网站| av网站在线播放免费| 午夜激情av网站| 青春草视频在线免费观看| 成人毛片60女人毛片免费| www.熟女人妻精品国产| 下体分泌物呈黄色| 80岁老熟妇乱子伦牲交| 亚洲精品日韩在线中文字幕| 我的亚洲天堂| 国产精品久久久久久精品古装| 国产一区二区三区综合在线观看| 精品国产乱码久久久久久小说| 欧美日韩国产mv在线观看视频| 一级爰片在线观看| 久久免费观看电影| 精品第一国产精品| 另类精品久久| 国产av精品麻豆| 久久久国产一区二区| 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| 日韩精品有码人妻一区| 毛片一级片免费看久久久久| 男女午夜视频在线观看| 狠狠婷婷综合久久久久久88av| 桃花免费在线播放| 国产亚洲av高清不卡| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 咕卡用的链子| 一本一本久久a久久精品综合妖精| 狂野欧美激情性xxxx| 曰老女人黄片| 亚洲欧美成人精品一区二区| 少妇被粗大猛烈的视频| 女人爽到高潮嗷嗷叫在线视频| 女人高潮潮喷娇喘18禁视频| 成人毛片60女人毛片免费| 色综合欧美亚洲国产小说| 韩国av在线不卡| 91精品三级在线观看| 亚洲国产av影院在线观看| 日本欧美国产在线视频| 在线观看免费视频网站a站| 最新的欧美精品一区二区| 国产视频首页在线观看| 999久久久国产精品视频| 精品国产一区二区三区久久久樱花| 大码成人一级视频| 丰满饥渴人妻一区二区三| 人体艺术视频欧美日本| √禁漫天堂资源中文www| 国产精品久久久久久久久免| 黑丝袜美女国产一区| 成人三级做爰电影| 亚洲三区欧美一区| 国产成人精品在线电影| 午夜激情久久久久久久| 国产精品蜜桃在线观看| 欧美xxⅹ黑人| 一级毛片电影观看| 国产精品久久久久久久久免| 黄频高清免费视频| 精品一区在线观看国产| 久久久久久久大尺度免费视频| 丝袜美腿诱惑在线| 搡老乐熟女国产| av线在线观看网站| 男的添女的下面高潮视频| 精品一区二区三区av网在线观看 | 夫妻午夜视频| 国产精品蜜桃在线观看| 制服丝袜香蕉在线| 久久久久网色| avwww免费| 一级毛片黄色毛片免费观看视频| 中国三级夫妇交换| 欧美国产精品va在线观看不卡| 80岁老熟妇乱子伦牲交| 在线观看免费午夜福利视频| 精品卡一卡二卡四卡免费| 午夜影院在线不卡| 天天操日日干夜夜撸| 成人国产av品久久久| 国精品久久久久久国模美| 亚洲国产精品成人久久小说| 男女床上黄色一级片免费看| 日本wwww免费看| 精品少妇一区二区三区视频日本电影 | 男人舔女人的私密视频| 日韩人妻精品一区2区三区| 王馨瑶露胸无遮挡在线观看| 国产av一区二区精品久久| 国产一区二区在线观看av| 亚洲欧美日韩另类电影网站| 女人久久www免费人成看片| 国产日韩一区二区三区精品不卡| 国产精品麻豆人妻色哟哟久久| 亚洲天堂av无毛| av女优亚洲男人天堂| 女人久久www免费人成看片| 欧美另类一区| 人人妻,人人澡人人爽秒播 | avwww免费| 日韩不卡一区二区三区视频在线| 亚洲欧美激情在线| 男的添女的下面高潮视频| 亚洲国产最新在线播放| 国产在线免费精品| 成人18禁高潮啪啪吃奶动态图| 亚洲第一av免费看| 免费观看a级毛片全部| 国产精品免费视频内射| 国产精品av久久久久免费| 午夜免费鲁丝| 午夜福利乱码中文字幕| 亚洲av男天堂| 飞空精品影院首页| 精品福利永久在线观看| 咕卡用的链子| 十八禁网站网址无遮挡| 丰满乱子伦码专区| 波多野结衣av一区二区av| 一区二区av电影网| 性色av一级| 欧美在线一区亚洲| 91精品伊人久久大香线蕉| 久久毛片免费看一区二区三区| 日本vs欧美在线观看视频| 亚洲三区欧美一区| 久久99一区二区三区| 99九九在线精品视频| 在线天堂中文资源库| 亚洲精品国产色婷婷电影| 亚洲国产毛片av蜜桃av| 成人国产麻豆网| 在线亚洲精品国产二区图片欧美| 亚洲激情五月婷婷啪啪| a级毛片黄视频| 精品国产露脸久久av麻豆| 午夜免费鲁丝| 欧美黑人精品巨大| 午夜老司机福利片| 欧美 日韩 精品 国产| 欧美日韩成人在线一区二区| a级毛片黄视频| 大香蕉久久网| 人人澡人人妻人| 精品一品国产午夜福利视频| 女人精品久久久久毛片| 18禁观看日本| 亚洲天堂av无毛| 亚洲欧美日韩另类电影网站| 你懂的网址亚洲精品在线观看| 精品卡一卡二卡四卡免费| 99re6热这里在线精品视频| 欧美少妇被猛烈插入视频| 波多野结衣av一区二区av| 日本一区二区免费在线视频| 黄色 视频免费看| 国产 一区精品| 亚洲av综合色区一区| 国产极品天堂在线| 亚洲欧美中文字幕日韩二区| 纵有疾风起免费观看全集完整版| 精品国产国语对白av| 欧美激情极品国产一区二区三区| 中文字幕最新亚洲高清| 考比视频在线观看| 国产男女内射视频| 国产成人系列免费观看| 精品一区二区免费观看| 午夜精品国产一区二区电影| 如日韩欧美国产精品一区二区三区| 麻豆乱淫一区二区| 亚洲激情五月婷婷啪啪| 曰老女人黄片| 亚洲七黄色美女视频| av电影中文网址| 欧美日韩综合久久久久久| 91精品国产国语对白视频| 国产在线视频一区二区| 两个人看的免费小视频| 蜜桃在线观看..| 久久久久人妻精品一区果冻| 一区二区日韩欧美中文字幕| 国产激情久久老熟女| 成人午夜精彩视频在线观看| 亚洲综合色网址| 午夜影院在线不卡| 一本色道久久久久久精品综合| 亚洲av在线观看美女高潮| 久久久久精品久久久久真实原创| 国产精品久久久久久久久免| 男人操女人黄网站| 蜜桃国产av成人99| 久久99精品国语久久久| 精品少妇一区二区三区视频日本电影 | 国产淫语在线视频| 精品午夜福利在线看| 在线观看人妻少妇| 女人精品久久久久毛片| 十八禁高潮呻吟视频| bbb黄色大片| 老司机靠b影院| 青草久久国产| 国产亚洲一区二区精品| 久久性视频一级片| 久久久精品区二区三区| 亚洲色图 男人天堂 中文字幕| 韩国高清视频一区二区三区| 在线观看免费日韩欧美大片| 国产精品一区二区精品视频观看| 男女下面插进去视频免费观看| 国精品久久久久久国模美| 亚洲欧美色中文字幕在线| 少妇的丰满在线观看| 亚洲成人免费av在线播放| 18禁裸乳无遮挡动漫免费视频| 18禁动态无遮挡网站| 看非洲黑人一级黄片| 亚洲精品国产av蜜桃| 视频区图区小说| 男女午夜视频在线观看| 桃花免费在线播放| 欧美精品高潮呻吟av久久| 一边亲一边摸免费视频| 考比视频在线观看| 两个人看的免费小视频| 亚洲男人天堂网一区| 我的亚洲天堂| 九草在线视频观看| 久久韩国三级中文字幕| 国产成人精品在线电影| 国产免费一区二区三区四区乱码| 精品国产一区二区三区久久久樱花| 美国免费a级毛片| 狠狠婷婷综合久久久久久88av| 一本大道久久a久久精品|