喻樹龍，李淑娟，姜盛夏，張同文，張瑞波，秦 莉，尚華明，張合理

（中國氣象局烏魯木齊沙漠氣象研究所，中國氣象樹木年輪理化研究重點(diǎn)開放實(shí)驗(yàn)室，新疆樹木年輪生態(tài)實(shí)驗(yàn)室，新疆 烏魯木齊830002）

樹木生長(zhǎng)不僅受自身遺傳因子的影響， 還受外界環(huán)境因子的支配， 其中海拔高度是影響樹木生長(zhǎng)的重要環(huán)境因子，隨海拔高度升高，水分和熱量條件產(chǎn)生變化，使樹木的生長(zhǎng)環(huán)境復(fù)雜多樣，進(jìn)而造成海拔梯度上樹輪徑向生長(zhǎng)的差異。 基于海拔高度變化對(duì)樹木年輪生長(zhǎng)影響的研究， 開始于20 世紀(jì)90 年代[1]。 國內(nèi)樹木氣候?qū)W學(xué)者曾在祁連山[2，3]、阿尼瑪卿山[4，5]、呂梁山[6]、川西高原[7]、蘆芽山[8]等區(qū)域開展過不同海拔高度樹木年輪氣候響應(yīng)研究， 發(fā)現(xiàn)同一坡面的不同海拔高度樹輪寬度指數(shù)年表特征呈一定的變化規(guī)律， 這是由于不同海拔高度水熱條件差異及高海拔樹木生長(zhǎng)季較短造成的。 在天山山區(qū)的西部[9-10]、中部[11]和東部[12]也分別開展過樹輪寬度梯度研究，朱海峰等[9]發(fā)現(xiàn)在天山北坡西部的伊犁地區(qū)， 地形對(duì)雪嶺云杉與氣候要素之間的關(guān)系影響較大。在較大空間范圍上，對(duì)比不同區(qū)域間樹木徑向生長(zhǎng)的梯度變化特征的相關(guān)研究開展較為有限， 本研究利用在中國境內(nèi)天山北坡5 個(gè)區(qū)域采集的樹輪梯度樣本， 探討天山北坡不同區(qū)域樹輪寬度年表特征在海拔梯度上的變化規(guī)律。

1 資料和方法

1.1 樹木年輪采樣和年表研制

天山北坡屬溫帶大陸性干旱半干旱氣候， 以雪嶺云杉（Picea schrenkiana）為主的原始森林主要集中在海拔1250~2900 m 的中山帶迎風(fēng)坡[13]，垂直地帶性比較明顯。分別在天山北坡的東部、中部和西部進(jìn)行樹木年輪梯度樣本采集（圖1），西部采樣區(qū)域在伊犁地區(qū)的特克斯林場(chǎng)的喬拉克鐵熱克（QL）和新源林場(chǎng)的卡勒克瑪（KL），中部在沙灣林場(chǎng)的大鹿角灣（DL）和石頭溝（ST），東部位于奇臺(tái)林場(chǎng)的寬溝（KG），采樣時(shí)間在2012—2015 年，各區(qū)域均選擇同一小流域，在森林分布帶內(nèi)間隔海拔高度100 m 左右的布設(shè)采樣點(diǎn)（表1）。 由于山區(qū)地形和森林分布差異，采樣區(qū)域的上下林緣海拔高度各不相同，各區(qū)域采樣點(diǎn)的個(gè)數(shù)也不相同， 其中采樣點(diǎn)數(shù)較少的新源林場(chǎng)的KL 和沙灣林場(chǎng)的DL 位于前山帶林區(qū)，KL 的上林緣在2100 m， DL 的下林緣在2155 m 左右。 5 個(gè)區(qū)域共選取35 個(gè)采樣點(diǎn)進(jìn)行采樣，每個(gè)采樣點(diǎn)均選擇20 棵以上樹木， 共采集了935 棵樹，每棵樹鉆取2 個(gè)樣芯。 5 個(gè)區(qū)域中，除了西部的QL 下林緣采樣點(diǎn)和中部的DL 2500 m 以下采樣點(diǎn)的坡度較緩，在5°~15°，其余采樣點(diǎn)均在20°以上（圖2），采樣區(qū)域林內(nèi)植被有草類和灌木， 土壤多為山地灰褐色森林土，各采樣點(diǎn)的郁閉度均＞0.3。

圖1 天山北坡樹輪采樣點(diǎn)和氣象站示意圖

樹木年輪樣本按照樹木年代學(xué)的基本原理和研究步驟進(jìn)行預(yù)處理，待樣本干燥后，進(jìn)行樣本粘貼固定、打磨、樹輪標(biāo)識(shí)，使用樹木年輪寬度測(cè)量系統(tǒng)進(jìn)行樹輪寬度測(cè)量， 利用折線圖法和國際年輪庫的COFECHA 程序進(jìn)行測(cè)量結(jié)果的交叉定年和質(zhì)量控制， 最后利用ARSTAN 程序建立35 個(gè)采樣點(diǎn)樹輪寬度標(biāo)準(zhǔn)化年表，取所有年表SSS 值均在0.85 以上的公共區(qū)間1932—2012 年進(jìn)行梯度分析。

1.2 研究方法

采用相關(guān)分析、 去趨勢(shì)對(duì)應(yīng)分析和主成分分析等方法分析不同海拔高度樹輪寬度年表間的一致性和差異性。對(duì)應(yīng)分析[14]也稱關(guān)聯(lián)分析、R-Q 型因子分析，是一種多元相關(guān)變量統(tǒng)計(jì)分析技術(shù)，通過分析由定性樣本構(gòu)成的交互匯總表來揭示樣本間的聯(lián)系，可以揭示各個(gè)樣本之間的差異。 由于在對(duì)應(yīng)分析中變量第二軸（縱軸） 的位置很大程度上依賴第一軸（橫軸）而產(chǎn)生的“拱形效應(yīng)”，采用區(qū)間去趨勢(shì)法進(jìn)行對(duì)應(yīng)分析[15]。

2 結(jié)果與分析

2.1 年表特征參數(shù)

天山北坡東部和中部采樣點(diǎn)的平均樹輪寬度均表現(xiàn)為隨海拔高度的升高而減少（圖2），平均樹齡呈增加的趨勢(shì)，與朱海峰等[9]研究一致，在西部的南天山北坡森林下樹線樹木生長(zhǎng)強(qiáng)于上樹線， 而QL和KL 沒有呈現(xiàn)相同規(guī)律， 可能是由于伊犁地區(qū)水分和熱量條件為天山北坡最好區(qū)域。 在這一區(qū)域沿海拔梯度采集樹齡較大的樹木時(shí)， 多選擇生長(zhǎng)在地形復(fù)雜、巖石多、樹木生長(zhǎng)處的土層較薄或坡度較大的區(qū)域，局地環(huán)境變化較大，因而未能較好地表現(xiàn)出輪寬和樹齡沿海拔高度變化的規(guī)律性。

表1 天山北坡樹木年輪采樣點(diǎn)概況

圖2 5 個(gè)區(qū)域不同海拔高度采樣點(diǎn)的平均坡度（a）、平均樹輪寬度（b）、平均樹齡（c）和年表長(zhǎng)度和樣本量（d，紅線藍(lán)線為年表長(zhǎng)度、黑線為樣本量）

2.1.1 平均敏感度

平均敏感度是衡量樹木生長(zhǎng)對(duì)氣候敏感性的重要指標(biāo)。 從圖3 可以看出，天山北坡中部的ST區(qū)域、東部KG 區(qū)域1800 m 以上區(qū)域均表現(xiàn)為隨海拔高度升高敏感度降低的趨勢(shì)， 與祁連山中部[3]和川西米亞羅林區(qū)[7]相同，在天山北坡東部的巴里坤林區(qū)[11]低海拔采樣點(diǎn)的平均敏感度要高于高海拔區(qū)域， 平均敏感度在0.2 以上的采樣點(diǎn)均出現(xiàn)在海拔2100 m 以下，說明在這些區(qū)域低海拔采樣點(diǎn)對(duì)氣候變化的敏感性要高于高海拔采樣點(diǎn)。 西部的QL 區(qū)域從2200~2300 m 處分為2 個(gè)部分， 呈現(xiàn)隨海拔升高呈降低—升高—降低的波狀趨勢(shì)， 這種趨勢(shì)與阿尼瑪卿山青海云杉寬度年表的平均敏感度[4]隨海拔變化的趨勢(shì)和呂梁山出現(xiàn)的森林中部的平均敏感度高于上下樹線的規(guī)律[6]基本一致。 天山北坡已有研究樹輪年表的特征參數(shù)[16~33]也表現(xiàn)為在2000 m 以上平均敏感度隨海拔高度的升高而降低。

2.1.2 標(biāo)準(zhǔn)差和信噪比

標(biāo)準(zhǔn)差為偏離平均水平的生長(zhǎng)偏差， 信噪比為氣候信號(hào)與氣候噪音的比值，二者的值越大，表明氣候?qū)淠镜挠绊懺斤@著， 樹輪年表中的氣候信息越多[34]。 從圖3 中可以看出，中部的DL 和ST 不同海拔采樣點(diǎn)的標(biāo)準(zhǔn)差均在0.25 以上， 而西部＞0.25 的值出現(xiàn)在森林中部，東部在低海拔區(qū)域。除了西部的QL，標(biāo)準(zhǔn)差最大值均出現(xiàn)在低海拔區(qū)域，在祁連山中部[3]也出現(xiàn)過同樣現(xiàn)象。 信噪比沒有出現(xiàn)與祁連山中部相同的梯度變化規(guī)律， 西部的KL 區(qū)域和東部KG 區(qū)域信噪比和標(biāo)準(zhǔn)差隨海拔高度變化與平均敏感度的變化基本相同。

2.1.3 一階自相關(guān)

反映上年氣候?qū)淠灸贻喩L(zhǎng)的持續(xù)性影響。5個(gè)區(qū)域自相關(guān)系數(shù)的顯著性水平在0.01 以上，其中東部的KG 和西部的KL 自相關(guān)系數(shù)最低， 中部的DL 和ST 大部分采樣點(diǎn)的自相關(guān)系數(shù)在0.7 以上，2400 m 以上的自相關(guān)系數(shù)在0.9 以上，2 個(gè)區(qū)域高海拔區(qū)域樹輪生長(zhǎng)受上年的氣候影響更顯著。 西部的QL 最大自相關(guān)系數(shù)出現(xiàn)在2400 m 的森林中部，上下林緣的自相關(guān)低于中部。 在天山北坡當(dāng)年樹輪生長(zhǎng)對(duì)上年氣候均有顯著響應(yīng)， 在2400 m 左右高海拔區(qū)域的森林中上部， 氣候?qū)漭喌某掷m(xù)性影響要強(qiáng)于其他海拔高度。

圖3 年表特征參數(shù)與海拔高度的關(guān)系

2.1.4 樹間相關(guān)系數(shù)

反映同一采樣點(diǎn)林木生長(zhǎng)受到氣候影響所導(dǎo)致的同步性變化大小， 相關(guān)系數(shù)越大表明對(duì)氣候的響應(yīng)越顯著，采樣點(diǎn)同步性越好，包含越多的氣候信息[34]。 從圖3e 中可以看出，天山北坡不同區(qū)域樹間相關(guān)系數(shù)沒有隨海拔高度變化的規(guī)律性， 未出現(xiàn)阿尼瑪卿山隨海拔升高而遞減[4]的趨勢(shì)。 5 個(gè)區(qū)域中，除了西部QL 的上樹線采樣點(diǎn)、 東部KG 的2400 m和1800~1900 m 幾個(gè)采樣點(diǎn)外，樹間相關(guān)系數(shù)都在0.7 以上，說明在天山山區(qū)高度差較小的區(qū)域內(nèi)取樣樹木徑向生長(zhǎng)的同步性更好。

2.1.5 樣本對(duì)總體的解釋信號(hào)

天山北坡EPS 區(qū)域平均值在0.9 以上，5 個(gè)區(qū)域中中部的DL、ST 和西部的KL 的EPS 在0.95 以上，而西部QL 的2500 m 以上和2100 m，東部KG的2200~2400 m 處采樣點(diǎn)低于0.9， 未表現(xiàn)出海拔梯度的規(guī)律性， 說明海拔高度變化不是影響樣本對(duì)總體的解釋信號(hào)的主要因素。

2.1.6 第一主成分方差解釋量

第一主成分方差解釋量是考量氣候因子對(duì)輪寬生長(zhǎng)的限制作用的重要指標(biāo)之一，百分比越大，限制作用越強(qiáng)，年表包含的氣候信息越多。區(qū)域平均最大值出現(xiàn)在2000~2200 m，除了西部的KL 是在1900 m出現(xiàn)最大值外， 其他3 個(gè)區(qū)域也在這個(gè)高度范圍內(nèi)出現(xiàn)極值，說明天山北坡森林中部2000~2200 m 的樹輪寬度年表氣候因子的限制最強(qiáng)，氣候信息最多。

2.2 樹輪年表與海拔梯度的關(guān)系

為了描述年表間的同步性關(guān)系， 計(jì)算5 個(gè)區(qū)域不同海拔高度間的互相關(guān)系數(shù)。 從互相關(guān)結(jié)果可以看出（圖4），5 個(gè)區(qū)域的高海拔采樣點(diǎn)中，西部的QL 2100 m 以上、中部ST 2300 m 以上和DL 2400 m 以上的高海拔采樣點(diǎn)相關(guān)較好， 而KG 并沒有表現(xiàn)出高海拔區(qū)域間互相關(guān)更好的規(guī)律性[11]，KG 的高海拔采樣點(diǎn)雖然坡度均在30°以上，但坡向有較大扭轉(zhuǎn)，坡向從NE—W—NW 變化，而其余區(qū)域的高海拔采樣點(diǎn)的坡向卻沒有顯著扭轉(zhuǎn)， 說明坡向?qū)Ω吆０螀^(qū)域樹輪寬度生長(zhǎng)的影響較大。在低海拔采樣點(diǎn)中，西部KL 和東部KG 的2000 m 以下低海拔采樣點(diǎn)間的互相關(guān)更好，在川西米亞羅林區(qū)[7]也出現(xiàn)同樣規(guī)律。 西部QL 的1890 m、 中部DL 2600 m 和ST 2206 m 等采樣點(diǎn)的坡度與相近采樣點(diǎn)有較大差別，這些采樣點(diǎn)與其他采樣點(diǎn)的互相關(guān)系數(shù)也較低，坡度也是樹輪生長(zhǎng)的影響因子。

圖4 不同區(qū)域樹輪寬度海拔梯度年表互相關(guān)分析

在相同海拔高度上（圖5），2500 m 以上的上樹線和高海拔區(qū)域采樣點(diǎn)間的樹輪寬度年表間均有較好的相關(guān)性， 說明在整個(gè)天山北坡高海拔和上樹線間的樹輪生長(zhǎng)一致性較好， 受大范圍的氣候變化影響。 而2200~2300 m 的森林中部采樣點(diǎn)間的相關(guān)較低，正、負(fù)相關(guān)均有出現(xiàn)，這是由于不同區(qū)域降水的局地性差異引起的森林中部樹輪徑向生長(zhǎng)的差異。2100 m 以下低海拔區(qū)域中，中部的ST 和東部的KG下林緣間的相關(guān)較高，而其他采樣點(diǎn)一致性較差，但均表現(xiàn)為正相關(guān)， 說明天山北坡東部和中部低海拔區(qū)域采樣點(diǎn)的一致性要好于西部。

為了進(jìn)一步了解海拔梯度上采樣點(diǎn)和不同距離的采樣點(diǎn)之間的差異， 利用去趨勢(shì)對(duì)應(yīng)分析分別對(duì)5 個(gè)區(qū)域樹輪采樣點(diǎn)進(jìn)行分析，并通過年表得分前2項(xiàng)（橫軸為第一項(xiàng)，縱軸為第二項(xiàng)）的正負(fù)值對(duì)不同海拔年表進(jìn)行分類（圖6、7）。從圖6 可以看出，西部QL 除了2760 m 的上樹線采樣點(diǎn)， 其他采樣點(diǎn)基于第一項(xiàng)年表得分的0 值線在2200~2300 m 分為兩類，而第二項(xiàng)得分的分類則與海拔高度無關(guān)。 QL 海拔高度在2600~2760 m 和2100~2200 m 采樣點(diǎn)的歐式距離較短，4 個(gè)采樣點(diǎn)間的差異小于2300~2500 m 3 個(gè)采樣點(diǎn)，這可能是由于前4 個(gè)采樣點(diǎn)的坡向?yàn)镹W，而后3 個(gè)采樣點(diǎn)坡向?yàn)閃NW，體現(xiàn)為較小坡向扭轉(zhuǎn)的影響。 其中2300 m 采樣點(diǎn)與其他采樣點(diǎn)差異最大，歐式距離均較長(zhǎng)，由于這個(gè)采樣點(diǎn)的平均敏感度和標(biāo)準(zhǔn)差是QL 區(qū)域內(nèi)最高的， 這可能是除了坡向和坡度外，包括土層厚度、石崖等其他小生境的影響造成。 KL 的第一項(xiàng)得分在1900~2000 m 將采樣點(diǎn)分為2 類， 其中1900 m 也表現(xiàn)為與區(qū)域內(nèi)其他采樣點(diǎn)的歐式距離均較長(zhǎng)， 這一采樣點(diǎn)的平均敏感度等年表特征均為區(qū)域最高值，坡度、坡向與相鄰采樣點(diǎn)無顯著差異， 平均年輪寬度和樹齡也是區(qū)域內(nèi)最大， 進(jìn)一步說明其他小生境條件也會(huì)對(duì)樹輪徑向生長(zhǎng)具有顯著的影響。中部DL 第一項(xiàng)得分的0 值線在2300~2400 m，ST 的正負(fù)值分界線在2095~2206 m，區(qū)域內(nèi)2300 m 以上的4 個(gè)采樣點(diǎn)和2095 m 以下3 個(gè)采樣點(diǎn)間歐式距離很短，而2206 m處的采樣點(diǎn)與其他采樣點(diǎn)距離均較長(zhǎng)， 這是由于該采樣點(diǎn)的坡度較相鄰采樣點(diǎn)的坡度小，平均敏感度和標(biāo)準(zhǔn)差也是區(qū)域內(nèi)最小。 天山北坡東部的KG 在2002~2098 m 采樣點(diǎn)分為兩類， 而第二項(xiàng)得分與其他區(qū)域相同也未表現(xiàn)出隨海拔高度變化的規(guī)律。 從上述分析可以看出， 如果不考慮未達(dá)到中山帶上下限樹線的DL 和KL 區(qū)域，其余3 個(gè)區(qū)域基于第一項(xiàng)得分的分界線表現(xiàn)出從西向東降低的規(guī)律性。 從采樣點(diǎn)間的歐式距離來看， 采樣點(diǎn)小生境的差異是采樣點(diǎn)間相似程度的重要影響因素。

圖5 天山北坡不同區(qū)域相同海拔高度年表間互相關(guān)分析

圖6 5 個(gè)采樣區(qū)域去趨勢(shì)對(duì)應(yīng)分析

利用5 個(gè)區(qū)域的所有采樣點(diǎn)進(jìn)行去趨勢(shì)對(duì)應(yīng)分析（圖7）。 從圖7 可以看出天山北坡西部的QL、KL和東部KG 區(qū)域內(nèi)各采樣點(diǎn)分布在上樹線采樣點(diǎn)的兩側(cè)且區(qū)域內(nèi)采樣點(diǎn)間的歐式距離較短， 而中部的DL 和ST 的采樣點(diǎn)分布較為分散，說明東西部樹輪寬度生長(zhǎng)有顯著的差異。 從海拔高度來看， 各區(qū)域2400 m 以上的高海拔采樣點(diǎn)多集中在第一項(xiàng)正值區(qū)內(nèi)且采樣點(diǎn)間的歐式距離較短，2000 m 以下采樣點(diǎn)則在0 值線和負(fù)值區(qū)內(nèi)但距離較長(zhǎng)，2300~2100 m的采樣點(diǎn)分布分散，說明在天山北坡，樹輪寬度生長(zhǎng)的一致性高海拔區(qū)域最好，低海拔區(qū)域次之，而森林中部一致性最差。

主成分分析結(jié)果表明（圖8）， 結(jié)果5 個(gè)區(qū)域中，西部KL 的第一主分量（PC1）的貢獻(xiàn)率最高，達(dá)到78%，西部的QL 次之，中部和東部的PC1 貢獻(xiàn)率在43%~48%，由于KL 采樣點(diǎn)均位于較低海拔，KL 年表間的互相關(guān)均較為顯著（圖4），對(duì)應(yīng)分析（圖7）也顯示KL 年表間的歐式距離較短，說明5個(gè)區(qū)域中KL 的年表間的一致性最好，西部要好于中部和東部。 5 個(gè)區(qū)域的載荷向量均為負(fù)值，西部的QL、中部的DL 和ST 均表現(xiàn)為高海拔采樣點(diǎn)載荷向量絕對(duì)值大于低海拔采樣點(diǎn)，而東部的KG 則相反，表明在同一區(qū)域存在同時(shí)影響不同海拔高度樹輪寬度形成的環(huán)境要素，西部和中部高海拔樹木對(duì)該要素的響應(yīng)更顯著，而東部的KG 因坡向扭轉(zhuǎn)等生境變化的影響顯示為低海拔區(qū)域響應(yīng)顯著。西部2 個(gè)區(qū)域的PC2 貢獻(xiàn)率在10%~12%， 中部和東部在25%~27%， 說明除了同時(shí)影響整個(gè)區(qū)域的環(huán)境要素外，中部和東部還存在其他影響較大的環(huán)境因子，而西部其他環(huán)境要素的影響則較小。 中部和東部的PC2 載荷向量沿海拔梯度的變化與PC1 相同， 中部和東部的高低海拔載荷向量呈反向變化，但載荷向量絕對(duì)值中部為低海拔大于高海拔，東部為高海拔大于低海拔，說明與PC1 相反，對(duì)次一級(jí)影響因子中部是低海拔采樣點(diǎn)響應(yīng)顯著，而東部是高海拔采樣點(diǎn)。 中部和東部PC3 貢獻(xiàn)率在10%以上，與西部的PC2 相同，說明在更為干旱的天山北坡中部和東部，對(duì)樹木徑向生長(zhǎng)的影響較大的環(huán)境因子要多于西部。

圖7 區(qū)域間和海拔高度去趨勢(shì)對(duì)應(yīng)分析

圖8 天山北坡不同區(qū)域載荷向量與海拔高度的關(guān)系

3 討論

3.1 海拔高度變化對(duì)樹輪徑向生長(zhǎng)的影響

在天山北坡隨著海拔高度升高氣溫隨之降低，降水量則存在最大降水帶， 一般認(rèn)為最多的高度在海拔1700～3200 m 附近，而各個(gè)季節(jié)降水量最大高度帶的高度是隨季節(jié)變化而變化， 一般為冷季低暖季高[35]。天山北坡的森林帶分布在最大降水帶，因而造成在海拔梯度上不同季節(jié)的影響樹輪生長(zhǎng)的氣候要素都各不相同， 這也是樹輪年表產(chǎn)生海拔梯度變化的主要原因。 天山北坡樹輪寬度、樹齡、敏感度等樹輪特征均存在隨海拔高度變化的現(xiàn)象，平均輪寬、平均敏感度與海拔高度呈反向變化， 低海拔采樣點(diǎn)對(duì)氣候變化響應(yīng)更敏感， 平均樹齡則與海拔高度為正相關(guān)，這是由于高海拔區(qū)域氣溫降低，樹木生長(zhǎng)期變短，樹輪較低海拔更不易生病和腐朽，因而每年的生長(zhǎng)量少且高樹齡樹多。 反映氣候因子對(duì)樹輪生長(zhǎng)影響的年表特征參數(shù)中， 標(biāo)準(zhǔn)差最大值均出現(xiàn)在低海拔區(qū)域，存在海拔高度上升標(biāo)準(zhǔn)差降低的規(guī)律。不同區(qū)域的信噪比、一階自相關(guān)、樹間相關(guān)和EPS 則未表現(xiàn)出顯著的海拔梯度變化規(guī)律。 第一主分量在2000~2100 m 左右的森林中部出現(xiàn)最大值， 這一高度也是春季和秋季最大降水帶[35]，氣候因子限制作用最強(qiáng)。

綜合互相關(guān)分析和對(duì)應(yīng)分析發(fā)現(xiàn)， 在同一梯度采樣區(qū)域， 高海拔采樣點(diǎn)間的樹輪生長(zhǎng)變化的同步性最好，但起始的海拔高度不同，例如西部在2100 m，中部在2300~2400 m，2000 m 以下的低海拔采樣點(diǎn)次之，2000~2400 m 左右的森林中部一致性最差。此外， 在對(duì)應(yīng)分析中還發(fā)現(xiàn)了一個(gè)由西到東降低的分界線，從2300 m 下降到了2000 m 左右，這一分界線可能與天山北坡最大降水帶東西部的分布差異有關(guān)。雖然天山北坡東西部采樣區(qū)域的距離在800 km以上，不同區(qū)域間在高海拔區(qū)域仍然有較好的同步性，特別是接近森林上樹限的2400 m 以上采樣點(diǎn)，說明高海拔區(qū)域的樹輪寬度年表和密度一樣，都可以反映較大范圍內(nèi)的氣候變化，由于高海拔采點(diǎn)樹輪寬度生長(zhǎng)多受氣溫的影響，溫度在空間的均一性遠(yuǎn)好于降水，進(jìn)而造成遠(yuǎn)距離樹輪采點(diǎn)間較好的一致性。 在低海拔區(qū)域，東部和中部下林緣采樣點(diǎn)間也有較好的一致性，而西部與中部和東部沒有顯著相關(guān)性，這可能與西部采樣點(diǎn)均位于伊犁河谷，低海拔區(qū)域降水量要遠(yuǎn)多于東部和中部， 低海拔區(qū)域樹輪徑向生長(zhǎng)多為降水影響， 因而造成區(qū)域間的差異。

3.2 坡向和坡度對(duì)樹輪的影響

坡向扭轉(zhuǎn)是天山北坡樹輪生長(zhǎng)同步性差異的重要因素， 東部KG 的高海拔采樣點(diǎn)坡向發(fā)生較大扭轉(zhuǎn)， 年表互相關(guān)系數(shù)顯著小于其他區(qū)域高海拔采樣點(diǎn)， 西部QL 海拔高度在2600~2760 m 和2100~2200 m 采樣點(diǎn)間有較小坡向扭轉(zhuǎn)，年表間的歐式距離則表現(xiàn)為坡向相近的更短， 在青藏高原阿尼瑪卿山東部山地也發(fā)現(xiàn)坡向扭轉(zhuǎn)是海拔梯度上影響祁連圓柏生長(zhǎng)變化的重要因子[5]。 在低海拔區(qū)域，西部KL 1660~1800 m 的坡向從NNW 轉(zhuǎn)向NE， 中部ST 1742~2095 m 的3 個(gè)采樣點(diǎn)坡向變化是ENE—NNW—N， 但年表間的互相關(guān)和歐式距離卻沒有顯著差異， 說明坡向扭轉(zhuǎn)可能在高海拔對(duì)樹輪生長(zhǎng)同步性影響更大。

坡度是平均敏感度的重要影響因子， 在天山北坡西部、 中部不同海拔高度均有出現(xiàn)坡度較大采樣點(diǎn)的平均敏感度大于相鄰采樣點(diǎn)的現(xiàn)象， 而坡度變化較大相鄰采樣點(diǎn)間互相關(guān)系數(shù)也較低， 例如西部QL 的1890 m、 中部DL 的2600 m 和ST 的2206 m等采樣點(diǎn)的坡向與相鄰采樣點(diǎn)接近， 但坡度相差在15°以上，與其他采樣點(diǎn)的互相關(guān)系數(shù)較低，因而坡度也是樹輪生長(zhǎng)的重要影響因子。 在天山北坡的不同海拔高度上，坡向、坡度等地形地貌是樹輪寬度生長(zhǎng)的重要影響因子， 會(huì)干擾海拔高度變化引起樹輪生長(zhǎng)垂直變化規(guī)律。

3.3 其他小生境的影響

在西部伊犁地區(qū)QL 的2300 m 采樣點(diǎn)和KL 的1900 m 采樣點(diǎn)與相鄰采樣點(diǎn)相關(guān)較為顯著，對(duì)應(yīng)分析中表現(xiàn)出與其他采樣點(diǎn)有較大差異， 歐式距離均較長(zhǎng)， 且2 個(gè)采樣點(diǎn)的平均敏感度和標(biāo)準(zhǔn)差是區(qū)域內(nèi)最高的，采樣點(diǎn)的坡向和坡度與相鄰采樣點(diǎn)相近，這可能是郁閉度、土壤、土層厚度、石崖等其他小生境對(duì)樹輪徑向生長(zhǎng)造成的影響。

4 結(jié)論

天山北坡不同區(qū)域的樹輪寬度、樹齡、敏感度等樹輪特征存在隨海拔高度變化的現(xiàn)象， 其中平均樹輪寬度、平均敏感度與海拔高度呈反向變化，平均樹齡則與海拔高度為正相關(guān)， 而其他特征參數(shù)無顯著的梯度規(guī)律。在同一區(qū)域內(nèi)和不同區(qū)域間，樹輪寬度生長(zhǎng)的一致性均表現(xiàn)為高海拔區(qū)域最好， 低海拔區(qū)域次之，中海拔區(qū)域最差。整個(gè)天山北坡高海拔和上樹線間的樹輪生長(zhǎng)一致性較好， 受大范圍的氣溫變化影響， 可以用于大空間尺度的樹木年輪氣候?qū)W研究中。 東部和中部低海拔區(qū)域采樣點(diǎn)的一致性要好于西部， 在天山區(qū)域東西部歷史氣候差異研究中可以在低海拔區(qū)域采集年代較長(zhǎng)的樹輪樣本。坡向、坡度等采樣點(diǎn)小生境的差異是影響沿海拔梯度的樹輪生長(zhǎng)垂直變化規(guī)律變異的重要要素。 在天山北坡同一區(qū)域存在同時(shí)影響不同海拔高度樹輪寬度形成的環(huán)境要素， 西部和中部高海拔樹木對(duì)該要素的響應(yīng)更顯著， 而東部因坡向扭轉(zhuǎn)等生境變化的影響顯示為低海拔區(qū)域響應(yīng)顯著， 在樹木年輪樣本采集時(shí)必須要考慮坡度和坡向的影響。 此外更為干旱的天山北坡中部和東部， 對(duì)樹木徑向生長(zhǎng)的影響較大的環(huán)境因子要多于西部伊犁地區(qū)。 本文中僅分析了海拔高度造成的樹輪寬度年表特征差異， 樹輪梯度年表的氣候響應(yīng)、 基于海拔梯度變化的樹輪氣候分異和空間氣候重建等方面的研究將在下一步進(jìn)行。