• <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
    很黄的视频免费| 精品久久久久久成人av| 国产国拍精品亚洲av在线观看| 嫩草影视91久久| 99久久成人亚洲精品观看| 久久99热6这里只有精品| 久久久国产成人免费| 久久久久国内视频| 久9热在线精品视频| 国产麻豆成人av免费视频| 久9热在线精品视频| 99在线视频只有这里精品首页| 亚洲成a人片在线一区二区| 啦啦啦啦在线视频资源| 亚洲精品456在线播放app | 国产精品久久久久久久电影| 熟妇人妻久久中文字幕3abv| 免费一级毛片在线播放高清视频| 亚洲久久久久久中文字幕| 国产在线精品亚洲第一网站| 啦啦啦观看免费观看视频高清| 淫妇啪啪啪对白视频| 国产精品电影一区二区三区| 韩国av一区二区三区四区| 久久久久久久久久成人| 精品福利观看| 不卡视频在线观看欧美| 久久人人精品亚洲av| 国产一级毛片七仙女欲春2| 男人舔奶头视频| eeuss影院久久| 欧美日韩中文字幕国产精品一区二区三区| 久久久久国内视频| 国产真实乱freesex| 简卡轻食公司| 国产中年淑女户外野战色| 亚洲国产色片| 可以在线观看的亚洲视频| 日本撒尿小便嘘嘘汇集6| 又爽又黄无遮挡网站| 久久草成人影院| 如何舔出高潮| 蜜桃亚洲精品一区二区三区| 干丝袜人妻中文字幕| 日本撒尿小便嘘嘘汇集6| 久久久久性生活片| 国产亚洲欧美98| 日本欧美国产在线视频| 国产精品久久久久久av不卡| 少妇猛男粗大的猛烈进出视频 | 少妇裸体淫交视频免费看高清| 最近视频中文字幕2019在线8| 极品教师在线免费播放| 色吧在线观看| 成人国产麻豆网| 国产精品av视频在线免费观看| 在线免费观看的www视频| avwww免费| 五月伊人婷婷丁香| 久久精品国产亚洲av香蕉五月| 婷婷六月久久综合丁香| 九九久久精品国产亚洲av麻豆| 久久久久久久久大av| 草草在线视频免费看| 日日撸夜夜添| 久久久久久伊人网av| 最新中文字幕久久久久| 午夜免费成人在线视频| 亚洲自拍偷在线| 亚洲美女黄片视频| 韩国av一区二区三区四区| 尤物成人国产欧美一区二区三区| 久久精品国产99精品国产亚洲性色| 99在线视频只有这里精品首页| 高清毛片免费观看视频网站| 最后的刺客免费高清国语| 夜夜爽天天搞| 国产成人影院久久av| 久久久精品欧美日韩精品| 老熟妇仑乱视频hdxx| 中国美白少妇内射xxxbb| 99久久成人亚洲精品观看| 成人二区视频| 最近在线观看免费完整版| www.www免费av| 好男人在线观看高清免费视频| 久久午夜亚洲精品久久| 国产精品一区www在线观看 | 欧美最黄视频在线播放免费| 亚洲成人久久爱视频| 国产精华一区二区三区| 精品一区二区三区视频在线| 99在线视频只有这里精品首页| 麻豆国产av国片精品| 国产真实乱freesex| 国产美女午夜福利| 美女高潮的动态| 久久精品国产99精品国产亚洲性色| 精品无人区乱码1区二区| 非洲黑人性xxxx精品又粗又长| 亚洲一区高清亚洲精品| 精品久久久久久久久久久久久| 男女那种视频在线观看| 成人精品一区二区免费| 少妇的逼水好多| 在线看三级毛片| 天美传媒精品一区二区| 啪啪无遮挡十八禁网站| 国内毛片毛片毛片毛片毛片| 亚洲无线在线观看| 99久久中文字幕三级久久日本| 搡女人真爽免费视频火全软件 | 人妻制服诱惑在线中文字幕| 白带黄色成豆腐渣| 天堂av国产一区二区熟女人妻| 别揉我奶头 嗯啊视频| 亚洲avbb在线观看| 国产在视频线在精品| av在线老鸭窝| 哪里可以看免费的av片| 我的女老师完整版在线观看| 色综合站精品国产| 69av精品久久久久久| 99久久无色码亚洲精品果冻| 成人午夜高清在线视频| 亚洲aⅴ乱码一区二区在线播放| 一本久久中文字幕| 1024手机看黄色片| 国产精华一区二区三区| 一进一出抽搐动态| 国产精品永久免费网站| 午夜日韩欧美国产| 三级男女做爰猛烈吃奶摸视频| 国产乱人视频| 久久99热这里只有精品18| 久久精品国产清高在天天线| 亚洲美女视频黄频| 欧美又色又爽又黄视频| 美女免费视频网站| 99在线人妻在线中文字幕| 99国产精品一区二区蜜桃av| 国产精品嫩草影院av在线观看 | 九九爱精品视频在线观看| 免费在线观看成人毛片| 成人性生交大片免费视频hd| 好男人在线观看高清免费视频| 男人的好看免费观看在线视频| 久久久久国内视频| 色精品久久人妻99蜜桃| 狠狠狠狠99中文字幕| 日本a在线网址| 舔av片在线| 男女视频在线观看网站免费| 色精品久久人妻99蜜桃| 听说在线观看完整版免费高清| 九色成人免费人妻av| 久久中文看片网| 日韩欧美一区二区三区在线观看| 3wmmmm亚洲av在线观看| 欧美日本亚洲视频在线播放| 日韩精品青青久久久久久| 人妻夜夜爽99麻豆av| 欧美日韩瑟瑟在线播放| 一个人看视频在线观看www免费| 国产亚洲精品av在线| 色综合婷婷激情| 赤兔流量卡办理| 成人精品一区二区免费| 看片在线看免费视频| 啦啦啦韩国在线观看视频| 天堂网av新在线| 日本一本二区三区精品| 国产精品一区二区免费欧美| 午夜久久久久精精品| 99久久成人亚洲精品观看| 久久久久精品国产欧美久久久| 久久久久精品国产欧美久久久| 一级av片app| 白带黄色成豆腐渣| 精品一区二区三区视频在线观看免费| 大型黄色视频在线免费观看| 久久人人精品亚洲av| 国产高清视频在线播放一区| 国内少妇人妻偷人精品xxx网站| 一进一出抽搐gif免费好疼| 18禁黄网站禁片免费观看直播| 国产男人的电影天堂91| 乱人视频在线观看| 国产精品女同一区二区软件 | 久久精品国产99精品国产亚洲性色| 久久人妻av系列| 久久久久免费精品人妻一区二区| 免费av不卡在线播放| 久久久午夜欧美精品| 99久久精品热视频| 乱人视频在线观看| 国产真实乱freesex| 国产v大片淫在线免费观看| 久久精品国产亚洲av天美| 少妇裸体淫交视频免费看高清| 91麻豆精品激情在线观看国产| 亚洲av免费在线观看| 丝袜美腿在线中文| 亚洲精品乱码久久久v下载方式| 在线播放无遮挡| 欧美又色又爽又黄视频| 国产成年人精品一区二区| 亚洲五月天丁香| 男插女下体视频免费在线播放| 两个人视频免费观看高清| 最近最新免费中文字幕在线| 桃色一区二区三区在线观看| 国产成人a区在线观看| 老熟妇乱子伦视频在线观看| 波多野结衣高清无吗| 国产高清激情床上av| 香蕉av资源在线| 变态另类丝袜制服| 国产精品自产拍在线观看55亚洲| 精品免费久久久久久久清纯| 18禁裸乳无遮挡免费网站照片| 可以在线观看的亚洲视频| aaaaa片日本免费| 夜夜夜夜夜久久久久| 亚洲美女搞黄在线观看 | 在线看三级毛片| 国产视频内射| 精品午夜福利视频在线观看一区| 亚洲不卡免费看| 成人亚洲精品av一区二区| 在线播放国产精品三级| 午夜精品一区二区三区免费看| 看片在线看免费视频| 特级一级黄色大片| 色综合色国产| 精品一区二区三区人妻视频| 亚洲精品亚洲一区二区| 亚洲国产精品sss在线观看| 久久99热这里只有精品18| 国产av不卡久久| 超碰av人人做人人爽久久| 午夜福利在线观看吧| 国产精品嫩草影院av在线观看 | 亚洲人成网站在线播| 99久久成人亚洲精品观看| 精品国内亚洲2022精品成人| 老师上课跳d突然被开到最大视频| a级一级毛片免费在线观看| 少妇熟女aⅴ在线视频| 亚洲精品日韩av片在线观看| 色哟哟·www| 一进一出抽搐动态| 亚洲avbb在线观看| 天堂影院成人在线观看| 国产成人a区在线观看| 亚洲av中文字字幕乱码综合| 免费大片18禁| 欧美xxxx性猛交bbbb| 真实男女啪啪啪动态图| 小蜜桃在线观看免费完整版高清| 国产乱人伦免费视频| 3wmmmm亚洲av在线观看| 69人妻影院| 蜜桃亚洲精品一区二区三区| 内射极品少妇av片p| 日韩一本色道免费dvd| 国产蜜桃级精品一区二区三区| 国产av不卡久久| 成人特级av手机在线观看| 国产黄色小视频在线观看| 俺也久久电影网| 在线播放国产精品三级| 免费无遮挡裸体视频| 国产av在哪里看| a级毛片a级免费在线| 成人av在线播放网站| 精品日产1卡2卡| 美女 人体艺术 gogo| 婷婷丁香在线五月| 天堂av国产一区二区熟女人妻| 一级黄色大片毛片| 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 一本一本综合久久| 一本精品99久久精品77| 亚洲电影在线观看av| 久久中文看片网| 免费观看的影片在线观看| 亚洲av不卡在线观看| 精品人妻一区二区三区麻豆 | 午夜精品久久久久久毛片777| 人妻久久中文字幕网| 天堂网av新在线| 一区二区三区高清视频在线| 色哟哟·www| 99热6这里只有精品| 搡老熟女国产l中国老女人| 美女大奶头视频| 校园人妻丝袜中文字幕| 老司机午夜福利在线观看视频| 99热这里只有是精品50| 亚洲欧美日韩高清在线视频| 国产亚洲精品综合一区在线观看| 欧美最新免费一区二区三区| 久久亚洲精品不卡| 国模一区二区三区四区视频| 九九爱精品视频在线观看| 伦理电影大哥的女人| 久久精品国产亚洲av涩爱 | 成人特级av手机在线观看| 给我免费播放毛片高清在线观看| 人妻少妇偷人精品九色| 国产久久久一区二区三区| 日本熟妇午夜| 神马国产精品三级电影在线观看| 露出奶头的视频| 国产视频内射| 欧美人与善性xxx| 国产极品精品免费视频能看的| av中文乱码字幕在线| 亚洲av美国av| 国产精品久久久久久av不卡| 麻豆av噜噜一区二区三区| 欧美一级a爱片免费观看看| 白带黄色成豆腐渣| 午夜影院日韩av| 久久久久久国产a免费观看| 最近最新免费中文字幕在线| aaaaa片日本免费| 成人一区二区视频在线观看| h日本视频在线播放| 少妇被粗大猛烈的视频| 欧美一区二区国产精品久久精品| 一本久久中文字幕| 国产色爽女视频免费观看| 男人舔奶头视频| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 99久久久亚洲精品蜜臀av| 嫩草影院精品99| 国产精品日韩av在线免费观看| 最后的刺客免费高清国语| 国产精品伦人一区二区| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 丝袜美腿在线中文| 国产 一区 欧美 日韩| 国产爱豆传媒在线观看| av视频在线观看入口| 尤物成人国产欧美一区二区三区| 亚洲经典国产精华液单| 少妇人妻一区二区三区视频| 亚洲av熟女| 一本精品99久久精品77| 久久精品人妻少妇| 两个人视频免费观看高清| 亚洲在线自拍视频| 日韩精品青青久久久久久| 欧美性感艳星| 99久国产av精品| 看黄色毛片网站| 亚洲电影在线观看av| 日本 av在线| 免费大片18禁| 蜜桃久久精品国产亚洲av| 欧美日韩黄片免| 99九九线精品视频在线观看视频| 干丝袜人妻中文字幕| h日本视频在线播放| 国产精品国产三级国产av玫瑰| 性插视频无遮挡在线免费观看| 亚洲一区高清亚洲精品| 国产精品99久久久久久久久| 亚洲精品在线观看二区| 精品久久久久久久久av| 在现免费观看毛片| 精品久久久久久久久久久久久| 国内精品美女久久久久久| or卡值多少钱| 亚洲av电影不卡..在线观看| 免费观看在线日韩| 亚洲狠狠婷婷综合久久图片| 免费看av在线观看网站| 黄色一级大片看看| 日韩欧美三级三区| 成年免费大片在线观看| 在线观看免费视频日本深夜| 床上黄色一级片| 三级毛片av免费| 波多野结衣巨乳人妻| 日韩欧美精品v在线| 男插女下体视频免费在线播放| 色5月婷婷丁香| 国产蜜桃级精品一区二区三区| 国产精品美女特级片免费视频播放器| 国产精品99久久久久久久久| 美女黄网站色视频| 色视频www国产| 日本熟妇午夜| 91av网一区二区| 18禁在线播放成人免费| 精品一区二区三区视频在线| 日日夜夜操网爽| 91久久精品电影网| 亚洲av一区综合| 亚洲在线观看片| 婷婷色综合大香蕉| 久久久国产成人免费| 久久久久国产精品人妻aⅴ院| 亚洲自偷自拍三级| 99精品久久久久人妻精品| 欧美最黄视频在线播放免费| 亚洲精品亚洲一区二区| 制服丝袜大香蕉在线| 日韩欧美在线二视频| 成年女人看的毛片在线观看| 国产一区二区激情短视频| 欧美丝袜亚洲另类 | 午夜福利欧美成人| 精品久久久久久久末码| 窝窝影院91人妻| 人人妻人人澡欧美一区二区| 高清毛片免费观看视频网站| netflix在线观看网站| 国内揄拍国产精品人妻在线| 桃红色精品国产亚洲av| 国产探花在线观看一区二区| 久久精品人妻少妇| 欧美日韩乱码在线| 欧洲精品卡2卡3卡4卡5卡区| 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| 亚洲欧美日韩高清专用| 人人妻人人看人人澡| 美女大奶头视频| 久久午夜福利片| 联通29元200g的流量卡| 人人妻,人人澡人人爽秒播| 天美传媒精品一区二区| 色尼玛亚洲综合影院| 亚洲一区二区三区色噜噜| 很黄的视频免费| 欧美最新免费一区二区三区| 精品免费久久久久久久清纯| 麻豆一二三区av精品| 美女高潮喷水抽搐中文字幕| 午夜福利成人在线免费观看| 人人妻,人人澡人人爽秒播| 久9热在线精品视频| 在线观看免费视频日本深夜| 国国产精品蜜臀av免费| 成人美女网站在线观看视频| 在线观看一区二区三区| 又爽又黄无遮挡网站| 大型黄色视频在线免费观看| 欧美成人免费av一区二区三区| 国内精品美女久久久久久| 51国产日韩欧美| 自拍偷自拍亚洲精品老妇| 亚洲av电影不卡..在线观看| 97超视频在线观看视频| 午夜日韩欧美国产| 男女之事视频高清在线观看| 又紧又爽又黄一区二区| 久99久视频精品免费| 亚洲乱码一区二区免费版| 特大巨黑吊av在线直播| 亚洲av免费在线观看| 赤兔流量卡办理| 观看美女的网站| 免费无遮挡裸体视频| 国产精华一区二区三区| 一区二区三区高清视频在线| 成人二区视频| 色综合亚洲欧美另类图片| 99热精品在线国产| 偷拍熟女少妇极品色| 精品福利观看| 久久久国产成人免费| 亚洲精品一区av在线观看| 婷婷丁香在线五月| 亚州av有码| 免费av观看视频| 日韩av在线大香蕉| 精品一区二区三区人妻视频| 亚洲一区高清亚洲精品| 久久久成人免费电影| 精品国内亚洲2022精品成人| 国产伦人伦偷精品视频| 一级黄色大片毛片| 国产麻豆成人av免费视频| 亚洲成人中文字幕在线播放| 十八禁国产超污无遮挡网站| 亚洲av日韩精品久久久久久密| 老司机深夜福利视频在线观看| 成人无遮挡网站| 12—13女人毛片做爰片一| 亚洲一区高清亚洲精品| 国产精品,欧美在线| 真人做人爱边吃奶动态| 国产高清有码在线观看视频| 国产午夜福利久久久久久| aaaaa片日本免费| 热99在线观看视频| 国产av麻豆久久久久久久| 嫩草影院新地址| 亚洲国产高清在线一区二区三| 热99re8久久精品国产| 成人三级黄色视频| 又黄又爽又刺激的免费视频.| 狂野欧美激情性xxxx在线观看| 色综合色国产| 成人性生交大片免费视频hd| 婷婷色综合大香蕉| 人妻丰满熟妇av一区二区三区| 偷拍熟女少妇极品色| 麻豆成人av在线观看| 欧美xxxx性猛交bbbb| 男人舔女人下体高潮全视频| 免费在线观看日本一区| 天美传媒精品一区二区| 永久网站在线| 久久草成人影院| 久久久久久久久久黄片| 色5月婷婷丁香| www.www免费av| 日日啪夜夜撸| 国产在视频线在精品| 噜噜噜噜噜久久久久久91| 黄色配什么色好看| 成人av在线播放网站| 欧美日韩精品成人综合77777| 亚洲在线观看片| 搡老妇女老女人老熟妇| 国产一区二区三区视频了| 国产亚洲欧美98| 久久久成人免费电影| 小蜜桃在线观看免费完整版高清| 性插视频无遮挡在线免费观看| 国产爱豆传媒在线观看| 九九爱精品视频在线观看| 麻豆国产97在线/欧美| 不卡视频在线观看欧美| 最近在线观看免费完整版| 美女免费视频网站| 欧美最黄视频在线播放免费| 亚洲aⅴ乱码一区二区在线播放| 两个人的视频大全免费| 欧美另类亚洲清纯唯美| 男女边吃奶边做爰视频| 亚洲av二区三区四区| 三级毛片av免费| 久久久成人免费电影| 国产精品人妻久久久久久| 一边摸一边抽搐一进一小说| 人人妻,人人澡人人爽秒播| 国产在视频线在精品| av在线天堂中文字幕| 中文亚洲av片在线观看爽| 欧美丝袜亚洲另类 | 日日夜夜操网爽| 一夜夜www| 人妻夜夜爽99麻豆av| 成熟少妇高潮喷水视频| 久久久久久国产a免费观看| 亚洲经典国产精华液单| bbb黄色大片| 亚洲美女视频黄频| 亚洲午夜理论影院| www.www免费av| 国产午夜精品久久久久久一区二区三区 | 久久99热这里只有精品18| 人人妻人人看人人澡| 99riav亚洲国产免费| 国产午夜精品久久久久久一区二区三区 | bbb黄色大片| 久久国产乱子免费精品| 国产精品一区二区免费欧美| 69人妻影院| 久久热精品热| 亚洲图色成人| 欧美性感艳星| 三级国产精品欧美在线观看| 亚洲av电影不卡..在线观看| a级毛片免费高清观看在线播放| 一级黄色大片毛片| 欧美3d第一页| 哪里可以看免费的av片| 99久久九九国产精品国产免费| 热99re8久久精品国产| 我的女老师完整版在线观看| 老司机深夜福利视频在线观看| 别揉我奶头~嗯~啊~动态视频| 草草在线视频免费看| 中国美白少妇内射xxxbb| 日本欧美国产在线视频| 最新中文字幕久久久久| 18+在线观看网站| 日韩欧美精品免费久久| 亚洲美女黄片视频| 免费在线观看日本一区| 色精品久久人妻99蜜桃| 亚洲美女黄片视频| 成人一区二区视频在线观看| 精品久久久久久久末码| 亚洲自拍偷在线| 美女被艹到高潮喷水动态| videossex国产| 亚洲美女黄片视频| 国产精品电影一区二区三区| 中国美女看黄片| 日本色播在线视频| 小蜜桃在线观看免费完整版高清| 久久天躁狠狠躁夜夜2o2o| 国产一区二区三区av在线 |