于德水 盧杰，2，3*

（1.西藏農(nóng)牧學(xué)院高原生態(tài)研究所，西藏 林芝860000；2.西藏高原森林生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，西藏 林芝860000；3.西藏林芝高山森林生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站，西藏林芝860000）

樹木年輪氣候?qū)W作為氣象學(xué)的一類運(yùn)用樹木年輪分析的溫度變化的分支學(xué)科，在19世紀(jì)至今的發(fā)展過(guò)程中日漸成熟，當(dāng)前已被廣泛運(yùn)用在各項(xiàng)氣象數(shù)據(jù)的重建，包括在溫度[1，2]、降水[3，4]、帕爾默干旱（PDSI）[5，6]、干燥指數(shù)（DI）[7，8]、Walter指數(shù)（WI）等[9]。通過(guò)樹木年輪的變化，運(yùn)用統(tǒng)計(jì)學(xué)原理進(jìn)行分析，以擴(kuò)充或延伸區(qū)域內(nèi)的氣象資料[10]。樹木作為區(qū)域環(huán)境變化的常見指示性植物，如全球范圍內(nèi)的ENSO現(xiàn)象[11]、La Nina現(xiàn)象、火山爆發(fā)[12]等，對(duì)森林生態(tài)系統(tǒng)產(chǎn)生影響的氣候變化，以及推測(cè)森林對(duì)區(qū)域內(nèi)氣候的反饋?zhàn)饔?。而近年利用樹木年輪分析森林火?zāi)[13]、地震[14]、山體滑坡[15]等災(zāi)害性環(huán)境因素，仍然作為樹木年輪學(xué)的關(guān)注熱點(diǎn)。樹木年輪氣候?qū)W作為從單一區(qū)域內(nèi)的氣候重建到全球樹木年輪網(wǎng)絡(luò)的重建統(tǒng)一[16]，當(dāng)前在全球針對(duì)樹木年輪氣候標(biāo)準(zhǔn)年表的重建正在逐步與全球網(wǎng)絡(luò)進(jìn)行相互補(bǔ)充并不斷的發(fā)展。

落葉松屬為樹木年輪氣候?qū)W的傳統(tǒng)研究對(duì)象，由于有著對(duì)溫度變化十分敏感的特性，在一段時(shí)間內(nèi)受到全球樹木年輪氣候?qū)W家的廣泛關(guān)注，本文選用當(dāng)前國(guó)內(nèi)研究較為深入、研究資料較為全面的兩種樹種西伯利亞落葉松（Larix sibirica）及太白紅杉（Larix potaniniivar.chinensis）進(jìn)行討論。

1 樹木年輪標(biāo)準(zhǔn)年表的建立

從Douglass[17]及Fritts[18]提出建立標(biāo)準(zhǔn)化年表，并確定標(biāo)準(zhǔn)化的主要內(nèi)容是去除樹木自身的生長(zhǎng)優(yōu)勢(shì)、樹木群落內(nèi)的干擾因素及環(huán)境中所存在的多種噪音后而建立的年表，無(wú)數(shù)的樹木氣候?qū)W家提出了大量的去除生長(zhǎng)勢(shì)數(shù)學(xué)模型，直至今日我們所運(yùn)用的常規(guī)趨勢(shì)標(biāo)準(zhǔn)仍然是上個(gè)世紀(jì)樹木年輪氣象學(xué)家Cook[19]所建立的“三步去勢(shì)法”。

1.1 生長(zhǎng)趨勢(shì)擬合

1985年Cook建立了標(biāo)準(zhǔn)化回歸模型，并開發(fā)出ARSTAN[20]軟件，使其囊括了Cook所建立的標(biāo)準(zhǔn)去勢(shì)方法。通過(guò)實(shí)用公式得出樹木年輪寬度指數(shù)（Ii）：

（1）式中Ti表示樹木年輪在LINTAB樹輪測(cè)序儀下的實(shí)測(cè)數(shù)值，Gi表示樹輪實(shí)測(cè)數(shù)值通過(guò)SPSS等統(tǒng)計(jì)軟件所計(jì)算出的期望值。Cook在實(shí)際中針對(duì)年輪較短的樹種發(fā)現(xiàn)，實(shí)用公式存在與數(shù)據(jù)偏離的結(jié)果。于是重新針對(duì)實(shí)用公式提出了自然對(duì)數(shù)公式

（2）式中InTi代表樹木實(shí)測(cè)樹輪的自然對(duì)數(shù)，InGi為樹輪期望值的自然對(duì)數(shù)，InIi則表示樹木寬度指數(shù)的自然對(duì)數(shù)，通過(guò)對(duì)公式對(duì)數(shù)化得到的數(shù)值將不會(huì)受時(shí)間空間的影響而存在數(shù)據(jù)偏移，但存在有兩種問(wèn)題，首先該對(duì)數(shù)公式極易受到年輪缺失影響，其次是得到的指數(shù)序列方差并不趨近于0.5的中心值導(dǎo)致得出的結(jié)論較為分散。

1.2 均值化

Cook在進(jìn)行去勢(shì)的過(guò)程中采用的雙權(quán)重平均法，將得到的年表指數(shù)（Iy）納入到標(biāo)準(zhǔn)年表（STD）中。

（3）式中Iy為年表指數(shù)，Ii為年輪寬度指數(shù)，k為常數(shù)，S為|Ii-Iy|的中位數(shù)。Cook在使用雙權(quán)重平均法建立STD年表時(shí)發(fā)現(xiàn)，對(duì)于少量的小樣本其精度得以保證，但當(dāng)樣本數(shù)過(guò)多（i＞6）時(shí)將會(huì)導(dǎo)致誤差過(guò)大，而此時(shí)選用第二種方式——算數(shù)平均法，將更加穩(wěn)健，同時(shí)在針對(duì)大量的樣本時(shí)，誤差將小于雙權(quán)重平均法。

當(dāng)前樹木氣候?qū)W家在進(jìn)行STD年表的制作過(guò)程中，常包含有大量的實(shí)驗(yàn)樣本，因此多選用算數(shù)平均法進(jìn)行數(shù)據(jù)處理。

1.3 STD年表的轉(zhuǎn)化

由于STD年表在建立完成后仍包含有大量的環(huán)境噪音，因此還需將STD年表作進(jìn)一步的處理。將樹木年輪中所包含的自相關(guān)指數(shù)重新返回至STD年表中進(jìn)行運(yùn)算，將得到差值年表RES，此時(shí)的RES年表在包含有外界氣候信號(hào)的同時(shí)，相比于STD年表去除了樹木自身的生長(zhǎng)優(yōu)勢(shì)及種內(nèi)競(jìng)爭(zhēng)等噪音信號(hào)的影響。實(shí)驗(yàn)經(jīng)驗(yàn)總結(jié)，此時(shí)進(jìn)行分析樹木年輪對(duì)氣候變化的響應(yīng)最為適宜。但RES年表在包含有氣候的影響的信號(hào)同時(shí)，仍然含有部分的森林群落大環(huán)境中的各項(xiàng)低頻的種間干擾因素，針對(duì)這個(gè)問(wèn)題樹木氣候?qū)W家開始思考如何解決種間干擾對(duì)樹木年輪所產(chǎn)生的影響。此時(shí)有學(xué)者提出建立自回歸標(biāo)準(zhǔn)化年表ARS，該年表是將樹種自身的信號(hào)響應(yīng)單獨(dú)提取，將信號(hào)回歸至RES年表，在進(jìn)行統(tǒng)計(jì)學(xué)處理后所得的ARS年表，不僅去除了自身的生長(zhǎng)因素的影響，而且還去除了森林小區(qū)域內(nèi)的干擾因素。但在去除干擾因素的同時(shí)也去除了一部分樹木對(duì)氣候變化響應(yīng)微弱的信號(hào)。因此Cook總結(jié)，在使用三種年表的過(guò)程中，應(yīng)先進(jìn)行初步的相關(guān)分析，通常相關(guān)度小于0.85的年表不適用于當(dāng)前區(qū)域或樹種。吳祥定指出針對(duì)中國(guó)的不同區(qū)域，在受地形的影響諸如在青藏高原，部分地區(qū)樹木年表的相關(guān)度可能不足0.75，但經(jīng)擬合后發(fā)現(xiàn)樹木年輪對(duì)氣候變化的相應(yīng)，卻有著極高的相關(guān)度。因此吳祥定提出，針對(duì)不同區(qū)域在建立樹輪年表時(shí)，相關(guān)度可以適時(shí)調(diào)整。[21]

2 西伯利亞落葉樹木年輪對(duì)氣候變化響應(yīng)的概況

西伯利亞落葉松（Larix sibirica）主要分布在新疆地區(qū)，同時(shí)也作為干旱半干旱地區(qū)樹木年輪氣候?qū)W的熱門樹種，研究表明，西伯利亞落葉松對(duì)區(qū)域內(nèi)的溫度變化存在著顯著的相關(guān)性[22]。而在深入的研究后西伯利亞落葉松被證明受到海拔的限制，表現(xiàn)出與氣溫及降水呈不同梯度的響應(yīng)變化[23]。

尚華明等人在21世紀(jì)初在新疆阿爾泰地區(qū)利用STD年表分析西伯利亞落葉松對(duì)6月溫度的變化，并提出與當(dāng)月的氣溫呈顯著正相關(guān)為開端。[24]國(guó)內(nèi)諸多學(xué)者相繼在新疆地區(qū)展開對(duì)西伯利亞落葉松的研究，從樹輪密度到同位素13C[25-26]，從生長(zhǎng)量到徑向生長(zhǎng)[27-28]。諸多學(xué)者分析西伯利亞落葉松生長(zhǎng)發(fā)育的過(guò)程中表現(xiàn)出對(duì)環(huán)境適應(yīng)性的生長(zhǎng)特征。直至崔宇等人運(yùn)用樹木年輪氣候?qū)W的研究方法，對(duì)西伯利亞落葉松樹木年輪過(guò)去的359a的初夏氣溫的重建，使用STD年表進(jìn)行分析，得出林線上樹木年輪寬度主要受氣溫的影響，通過(guò)分析STD年表后發(fā)現(xiàn)研究區(qū)域內(nèi)氣候變化的冷暖期及其變化規(guī)律。[29]張晴等人在天山東部采用RES年表，詳細(xì)地分析了在天山地區(qū)不同海拔梯度下的西伯利亞落葉松與氣候變化的關(guān)系，并發(fā)現(xiàn)氣候變暖使樹木徑向生長(zhǎng)不斷減小的特征，在低海拔地區(qū)特征表現(xiàn)更為顯著。在海拔上中下限西伯利亞落葉松對(duì)氣溫升高的敏感性逐步降低，高海拔地區(qū)西伯利亞落葉松的徑向生長(zhǎng)主要受溫度的影響，中低海拔地區(qū)主要受降水與溫度的共同影響。[30]隨后的幾年內(nèi)學(xué)者們運(yùn)用不同的研究方式，不同的研究手段，試圖將新疆地區(qū)整體氣候變化進(jìn)行統(tǒng)一，尋找西伯利亞落葉松樹木年輪與氣溫關(guān)系最顯著的月份，并希望將這種變化的響應(yīng)信號(hào)歸入到全球氣溫變化的研究范圍中。直至焦亮等人在阿爾泰山中段分析從1798年以來(lái)6～7月的氣溫變化，證明了6～7月平均最低氣溫是研究區(qū)內(nèi)樹木年輪徑向生長(zhǎng)的主控氣候因子引證了前人的研究觀點(diǎn)，并證明研究區(qū)內(nèi)的氣候變化主要受厄爾尼諾-南方驅(qū)動(dòng)（ENSO）、太陽(yáng)黑子活動(dòng)及全球溫度變化（SST）的驅(qū)動(dòng)，在國(guó)內(nèi)率先將西伯利亞落葉松樹木年輪徑向生長(zhǎng)與氣候變化的關(guān)系擴(kuò)大至全球。[31]

西伯利亞落葉松對(duì)氣候變化的響應(yīng)在新疆地區(qū)主要是受到當(dāng)年6月氣溫的影響并呈顯著的正相關(guān)性，且與當(dāng)年6月的降水呈顯著的負(fù)相關(guān)性。[32]但也有學(xué)者研究發(fā)現(xiàn)，西伯利亞落葉松的徑向生長(zhǎng)對(duì)與上一年8月和當(dāng)年8月的平均氣溫、最高氣溫以及當(dāng)年8月的最低氣溫呈顯著負(fù)相關(guān)。[33]在不同海拔梯度的研究中發(fā)現(xiàn)，在高緯度地區(qū)西伯利亞落葉松主要受到溫度的限制，而在低海拔地區(qū)則主要受到降水的限制。[34]在天山地區(qū)也有學(xué)者對(duì)西伯利亞落葉松的徑向生長(zhǎng)進(jìn)行復(fù)本研究，并采用RES年表得出了與新疆地區(qū)的研究學(xué)者相似的研究結(jié)果。更有學(xué)者運(yùn)用ArcGIS程序?qū)SO模型引入到樹木年輪徑向生長(zhǎng)對(duì)氣候響應(yīng)關(guān)系的模擬中，且模擬效果同樣印證了阿爾泰山與天山兩地區(qū)的研究成果。[35]

通過(guò)對(duì)西伯利亞落葉松在阿爾泰山與天山兩地區(qū)的研究對(duì)比發(fā)現(xiàn)，在研究區(qū)域范圍內(nèi)，樹木年輪與氣候變化的研究實(shí)驗(yàn)中最先使用STD年表在大范圍內(nèi)開展概況研究。因?yàn)镾TD年表包含有大量的響應(yīng)因素，所以針對(duì)區(qū)域內(nèi)的氣候變化有著更加廣泛的包容性。同時(shí)由于其中還包含有樹木與區(qū)域范圍內(nèi)的相關(guān)信號(hào)，因此在研究的過(guò)程中，所關(guān)聯(lián)的響應(yīng)信號(hào)可以擴(kuò)大至區(qū)域內(nèi)的所有氣候及干擾因素，包括自然干擾在內(nèi)的所有信號(hào)，均可采用STD年表進(jìn)行討論。在針對(duì)氣候變化的某單一環(huán)境因素時(shí)，可以通過(guò)采用RES年表對(duì)氣候變化的因素進(jìn)行相關(guān)性分析，此時(shí)采用STD年表將會(huì)受到區(qū)域內(nèi)其他因素的干擾而造成信號(hào)的衰減而使敏感度降低。而在采用ARS年表進(jìn)行研究分析的過(guò)程中要將研究的信號(hào)單獨(dú)進(jìn)行分析，由于其去除了大量的外界干擾，最大程度的表現(xiàn)出單一氣象因子的影響，但在建立ARS年表的過(guò)程中也去除了部分樹木對(duì)氣候變化的信號(hào)，而使得ARS年表表現(xiàn)出對(duì)氣候信號(hào)的相關(guān)性降低，因此在使用三種年表的過(guò)程中要根據(jù)研究?jī)?nèi)容進(jìn)行預(yù)先評(píng)估，并在其最后采用COFECHA程序進(jìn)行信號(hào)敏感度的信任關(guān)系檢驗(yàn)。

3 太白紅杉樹木年輪對(duì)氣候變化響應(yīng)的概況

太白紅杉（Larix potaninii var.chinensis）作為陜西、山東等我國(guó)中部地區(qū)的優(yōu)勢(shì)樹種，主要分布在秦嶺太白山周邊地區(qū)，作為國(guó)內(nèi)中部地區(qū)研究樹木年輪與氣候變化的主要樹種。當(dāng)前研究證明太白紅杉對(duì)氣候變化的響應(yīng)十分復(fù)雜，與當(dāng)年溫度、降水和上一年溫度、降水均存在有相關(guān)性，與全球氣溫變化及太平洋濤動(dòng)（ODP）也存在有顯著的相關(guān)性。本世紀(jì)初，吳建國(guó)在分析太白山地區(qū)7種植物的分布性研究時(shí)發(fā)現(xiàn)太白紅杉在徑向生長(zhǎng)的過(guò)程中受到當(dāng)年及上一年氣溫及降水的影響。[36]但當(dāng)時(shí)還并不清楚具體時(shí)期及相應(yīng)的響應(yīng)規(guī)律，但此時(shí)已有學(xué)者提出隨著全球氣溫的不斷變化，在未來(lái)太白紅杉的適宜區(qū)將可能逐步衰減直至消失。由康永祥等開展太白紅杉樹木年輪對(duì)氣候變化響應(yīng)的研究，但此時(shí)并未運(yùn)用樹木氣象學(xué)的研究方法，而是采用函數(shù)分析的方式對(duì)太白紅杉樹木徑向生長(zhǎng)與氣候變化的相關(guān)性進(jìn)行函數(shù)分析，并提出太白紅杉對(duì)當(dāng)年4月的平均氣溫有顯著的正相關(guān)性，與當(dāng)年8月的降水呈顯著的負(fù)相關(guān)性，而與當(dāng)年10月的降水呈顯著的正相關(guān)性。[37]由于未采用樹木標(biāo)準(zhǔn)年表的研究形式，因此所提出的研究結(jié)論存在有一定的誤差，在之后的研究過(guò)程中，其觀點(diǎn)作為研究基礎(chǔ)，被學(xué)者們不斷地更新。2015年由秦進(jìn)等在太白山南北兩坡開展實(shí)驗(yàn)，并在分析樹木年輪的三種年表所包含的信息后，采用RES年表與太白紅杉的樹木年輪對(duì)氣候變化的響應(yīng)進(jìn)行分析，并得出在研究區(qū)內(nèi)太白紅衫的徑向生長(zhǎng)主要受前一年6月降水影響的結(jié)論。針對(duì)不同的坡向，證明：①在南坡樹木年輪主要與當(dāng)年5～6月的平均氣溫呈顯著的正相關(guān)，與當(dāng)年1～4月的降水呈顯著負(fù)相關(guān)；②在北坡樹木年輪與當(dāng)年及前一年1～6月的平均氣溫呈顯著正相關(guān)，與當(dāng)年8月的降水呈顯著負(fù)相關(guān)。同時(shí)研究發(fā)現(xiàn)北坡樹木年輪對(duì)氣溫的變化更加敏感。[38]這也是在研究區(qū)域內(nèi)最先采用樹木年輪建立樹木年表的研究方式，隨后大量學(xué)者開展研究實(shí)驗(yàn)所得的結(jié)論也印證了秦進(jìn)等人的觀點(diǎn)，奠定了研究當(dāng)?shù)貐^(qū)域內(nèi)樹木徑向生長(zhǎng)與氣候變化的基礎(chǔ)。而針對(duì)太白山地區(qū)樹木徑向生長(zhǎng)與氣候變化開展了深入的研究，有研究表明在太白山區(qū)域內(nèi)的氣溫變化與全球氣溫的變化存在著顯著的相關(guān)性。而在太白山地區(qū)的樹木徑向生長(zhǎng)也受到全球氣溫變化的影響。[39]印證了在未來(lái)太白紅杉的生存環(huán)境必將遭受嚴(yán)酷的考驗(yàn)。

在太白山進(jìn)行的氣候變化的研究同樣在秦嶺鰲山地區(qū)被重復(fù)，并在研究的過(guò)程中采用了STD年表，同樣得出相似的研究結(jié)論[40]，在不同研究區(qū)域中同樣印證了全球氣溫變化后，太白紅杉的適宜區(qū)將逐步減小直至消失[41]。

對(duì)太白紅杉的研究結(jié)果及其研究方法簡(jiǎn)單分析后，可發(fā)現(xiàn)在探究區(qū)域氣溫與當(dāng)?shù)貥浞N的生長(zhǎng)變化的響應(yīng)性研究中，最先采用的是函數(shù)公式的模擬研究，但研究結(jié)論較為廣泛，還需深入分析，而采用年表分析的過(guò)程中，尤其是采用ARS年表，由于在去除大量的環(huán)境噪聲的同時(shí)破壞了環(huán)境中的部分氣候信號(hào)，致使對(duì)氣溫的變化無(wú)法通過(guò)信號(hào)的檢驗(yàn)。而在研究的程中采用ArcGIS與STD年表的結(jié)合同樣可以得到相當(dāng)?shù)臏?zhǔn)確結(jié)論。[42，43]

對(duì)太白紅杉的研究過(guò)程進(jìn)行分析，發(fā)現(xiàn)太白山地區(qū)與太白紅杉相關(guān)的研究，不單證明區(qū)域內(nèi)樹木徑向生長(zhǎng)對(duì)氣候變化的關(guān)系，同時(shí)也提供了豐富的研究手段及思路，擴(kuò)展了樹木年輪氣候?qū)W的研究方法。

4 結(jié)語(yǔ)

樹木年輪氣候?qū)W的提出最先是用于氣象數(shù)據(jù)的印證與擴(kuò)展，在隨著技術(shù)手段的不斷發(fā)展，Cook的研究模式被廣泛地運(yùn)用在樹木徑向生長(zhǎng)對(duì)氣候變化的各個(gè)階段，而隨后的學(xué)者諸如Blasing[44]、Robinson[45]、Holmes[46]等不斷地優(yōu)化分析公式，將樹木年表不斷地深化，并對(duì)Cook等人所開發(fā)ARSTAN程序不斷地更新，將上世紀(jì)樹木生態(tài)學(xué)家所提出的研究公式及方法進(jìn)行整合，隨后創(chuàng)建的COFECHA程序，使樹木年輪學(xué)家在建立樹木年輪STD年表、RES年表、ARS年表等三種標(biāo)準(zhǔn)年表時(shí)更為精準(zhǔn)、更為系統(tǒng)，不斷地豐富學(xué)者們的研究方法，并不斷地提供新的研究思路。通過(guò)對(duì)兩種樹木徑向生長(zhǎng)與氣候變化的文獻(xiàn)分析，當(dāng)前的研究進(jìn)展仍然是對(duì)研究區(qū)域內(nèi)的特有樹種的徑向生長(zhǎng)與氣候變化的響應(yīng)進(jìn)行分析，并運(yùn)用不同的分析手段，探究研究區(qū)域內(nèi)的氣候變化，并將區(qū)域內(nèi)樹木年輪對(duì)氣候的響應(yīng)結(jié)論延伸至全球氣候變化的大范圍內(nèi)。