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

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

(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