樹輪δ13C記錄的阿勒泰地區(qū)近160 a夏季氣溫變化

張瑞波，尚華明，魏文壽，袁玉江*，喻樹龍，張同文，范子昂，陳 峰，秦 莉

（1.中國氣象局烏魯木齊沙漠氣象研究所；新疆維吾爾自治區(qū)樹木年輪生態(tài)重點(diǎn)實(shí)驗(yàn)室；中國氣象局樹木年輪理化研究重點(diǎn)實(shí)驗(yàn)室，新疆 烏魯木齊830002；2.中國科學(xué)院寒區(qū)旱區(qū)環(huán)境與工程研究所冰凍圈科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，甘肅 蘭州730000）

阿爾泰山全長(zhǎng)2 100 km，位于中國、哈薩克斯坦、俄羅斯和蒙古交界段，呈西北—東南走向橫亙于亞洲中部[1]。阿爾泰山森林資源豐富，在海拔1 400～2 400 m的山區(qū)最大降水帶分布有西伯利亞落葉松（Larix sibirica Ledeb.）原始森林，俄羅斯學(xué)者等較早在阿爾泰山北坡開展了樹輪氣候?qū)W研究[2]。20世紀(jì)80年代，李江風(fēng)、袁玉江等就在中國境內(nèi)阿爾泰山南坡東部開展了大量的樹輪氣候和水文學(xué)研究，建立了該地區(qū)的樹輪年表序列，重建了溫度、降水和額爾齊斯河徑流量等氣候水文長(zhǎng)序列[3]，但由于采樣條件的限制，早期的采樣點(diǎn)大多位于海拔較低的森林下限，這里樹輪寬度的主要限制因子多為降水。近幾年來，通過對(duì)該區(qū)域森林上限區(qū)域的樹輪研究發(fā)現(xiàn)，位于阿爾泰山南坡森林上限的西伯利亞落葉松寬度主要受生長(zhǎng)季溫度的影響[4]。隨著樹木年輪研究技術(shù)的發(fā)展和國內(nèi)實(shí)驗(yàn)室條件的改善，樹輪研究已不再僅僅局限于樹輪寬度，在樹輪密度[5]，樹輪同位素[6-7]，樹輪圖像[8]等多個(gè)方面獲得很大的進(jìn)步。在阿爾泰山，有大量的樹輪寬度[9-13]、密度研究結(jié)果[14-15]，利用樹輪穩(wěn)定同位素研究重建該地區(qū)過去氣候變化研究較少，植物中碳穩(wěn)定同位素的研究已被廣泛應(yīng)用于植物生理、生態(tài)環(huán)境及全球環(huán)境變化等領(lǐng)域[16]。樹輪穩(wěn)定碳同位素以能夠提供分辨率到年甚至季節(jié)的優(yōu)勢(shì)，在研究過去環(huán)境變化及全球碳循環(huán)方面[17-18]具有重要意義。張瑞波等（2012）分析了阿爾泰山東段樹輪穩(wěn)定碳同位素對(duì)氣候的響應(yīng)[19]，本研究利用采自阿爾泰山南坡的樹輪樣本，提取樹δ13C數(shù)據(jù)，利用樹輪δ13C重建了阿勒泰地區(qū)過去160 a夏季平均氣溫，并分析了氣候變化特征，為揭示阿爾泰山過去氣候變化提供新資料。

1 資料和方法

1.1 研究區(qū)概況及樣本采集

阿爾泰山脈是亞洲宏偉山系之一，該區(qū)域?qū)俦睖貛Т箨懶詺夂?，具有冬寒漫長(zhǎng)，春旱多風(fēng)，夏短炎熱，秋涼氣爽的氣候特點(diǎn)。在海拔1 400～2 400 m的山區(qū)最大降水帶分布有西伯利亞落葉松（Larix sibirica Ledeb.），該樹種耐干旱、嚴(yán)寒，一般5月發(fā)芽，6—7月為速生期，9月開始落葉進(jìn)入休眠期。

樹木年輪樣本分別采自位于阿爾泰山西部哈巴河林場(chǎng)的齊巴阿克什 QBA （47°25′N，86°30′E；1 660～1 720 m），大致接近森林下線，選取25棵樹，分別采集了50根粗芯（用于樹輪穩(wěn)定同位素分析）和50根細(xì)芯（用于樹輪寬度分析），位于阿爾泰山中部的阿勒泰林場(chǎng)的科勒迭難布拉克KLD（47°59′N，88°22′E；2 300～2400 m），大致接近森林上線，選取31棵樹，分別采集了30根粗芯和61根細(xì)芯，位于阿爾泰山東部的富蘊(yùn)林場(chǎng)卡依爾特后山的正格ZGE（47°42′N，89°53′E；2 090～2250 m），選取 25 株樹，分別采集了50根粗芯和50根細(xì)芯。

1.2 樹輪寬度和δ13C序列的建立

3個(gè)采樣點(diǎn)樹木年輪寬度年表的研制經(jīng)過以下過程：（1）按照樣本的基本處理程序，對(duì)樹輪細(xì)芯樣本進(jìn)行干燥、固定、打磨、交叉定年，用精度為0.001 mm的輪寬測(cè)量?jī)x和MeasureJ2X程序進(jìn)行輪寬測(cè)量；（2）用國際年輪庫的COFECHA定年質(zhì)量控制程序進(jìn)行交叉定年的檢驗(yàn)[20]；（3）采用ARSTAN年表研制程序完成QBA、KLD和ZGE3個(gè)采點(diǎn)樹輪寬度年表的建立[21]。

不少研究者發(fā)現(xiàn)，盡管全木和纖維素同位素分析在絕對(duì)值明顯不同，但是二者的變化非常相似[22]，本研究提取樹輪全木的δ13C。Leavitt等[23]發(fā)現(xiàn)當(dāng)樹木樣品數(shù)目從一個(gè)增加到4個(gè)時(shí)，Pearson相關(guān)系數(shù)和累積平均系列與真實(shí)系列差的平方和分別從0.89和1.54改善到0.96和0.47，但超過4個(gè)樣本后，上述兩個(gè)參數(shù)幾乎沒有改變。考慮到樣本剝離的難度及復(fù)本量原則，分別選取QBA采點(diǎn)9棵、KLD采點(diǎn)10棵和ZGE采點(diǎn)8棵年輪較寬、缺輪較少，敏感度較高，沒有異常的樹輪粗樣芯，對(duì)比寬度年表交叉定年，確定每棵粗芯準(zhǔn)確日歷年并做兩面標(biāo)記。從樹皮到樹心依次在顯微鏡下用手術(shù)刀對(duì)樹木樣芯逐輪剝離，為防止碳源污染，樣品在玻璃墊板上切割?？紤]到可能存在的幼齡效應(yīng)[24]的影響，舍去靠近樹心的20 a；將樣芯中相同年份的樣品混合，置于已按照年代編號(hào)的離心管中，在60℃下低溫干燥24 h后，用混合研磨儀將樣品磨細(xì)至100目，稱取100 μg左右的全樣樣品，用小錫杯包好，同時(shí)稱量實(shí)驗(yàn)室標(biāo)準(zhǔn)樣，每10 a加一個(gè)實(shí)驗(yàn)室標(biāo)準(zhǔn)樣以確定系統(tǒng)誤差，利用FLASH HT在950℃條件下通氧燃燒，生成的氣體經(jīng)純化后經(jīng)ConfloIII進(jìn)入Delta V Advantage穩(wěn)定氣體同位素質(zhì)譜儀分析其δ13C組成，δ13C分析的系統(tǒng)誤差小于0.2‰。樹輪δ13C分析結(jié)果以δ13CPDB表示[25]：δ13CPDB（‰）=[（13C/12C）樣品/（13C/12C）標(biāo)準(zhǔn)－1]×1 000。

自工業(yè)革命以來，大氣CO2濃度持續(xù)升高，大氣 CO2中 δ13C下降會(huì)影響植物組織δ13C的變化[26-27]。為了真實(shí)地反映樹輪δ13C所記錄的氣候變化信息，在探討樹輪δ13C與氣候要素關(guān)系前有必要去除大氣CO2濃度持續(xù)升高所帶來的影響；按照校正工業(yè)革命前的水平（-6.4‰）進(jìn)行校正，δ13Cr=δ13Cplant－（δ13Catm+6.4），得到校正后的穩(wěn)定碳同位素序列，其中δ13Cplant和δ13Catm分別為樹輪實(shí)測(cè)δ13C值和大氣背景δ13C值。上述校正過程中的1850—2003年的大氣背景δ13C值來源于McCarroll等[28]，2004-2008年的數(shù)據(jù)根據(jù)1973—2003年的數(shù)據(jù)進(jìn)行線性趨勢(shì)外推得到，由此可以得到阿爾泰山南坡東中西部三條樹輪δ13C序列（圖1）。

1.3 氣象資料

圖1 阿爾泰山南坡δ13C趨勢(shì)序列

本研究所用的氣象資料來源于中國氣象局國家氣候中心，選取具有代表性、時(shí)間序列連續(xù)、資料較為完整、可信度較高的阿勒泰地區(qū)阿勒泰（47°44′N，88°05′E，海拔 737.7 m）、布爾津（47°42′N，86°52′E，海拔 475.5 m）、福海（47°04′N，87°28′E，海拔502.8 m）、富蘊(yùn)（46°59′N，89°31′E，海拔 810.5 m）、哈巴河（48°03′N，86°24′E，海拔 534.5 m）、吉木乃（47°26′N，85°52′E，海拔 983.9 m）和青河（46°40′N，90°23′E，海拔 1 220.0m）等 7 個(gè)氣象站 1962—2012年平均氣溫、平均最高氣溫、平均最低氣溫和降水量逐月資料，采用Mann-Kendell方法和Double-mass analysis方法對(duì)該站的溫度和降水資料進(jìn)行了均一性檢驗(yàn)[29]，結(jié)果表明溫度和降水記錄無隨機(jī)突變和明顯不均勻分布情況，這7個(gè)氣象站自建站以來到現(xiàn)在也沒有遷站記錄。觀測(cè)資料顯示（1962—2012），中國境內(nèi)的阿爾泰山南坡阿勒泰地區(qū)年均溫3.76℃，年降水量175 mm，近50 a來，阿勒泰地區(qū)降水和溫度具有不同程度的上升趨勢(shì)，降水量呈雙峰型，但仍然以夏季為主（圖2）。

空間代表性分析降水資料來自英國East Anglia大學(xué)的Climatic Research Unit（簡(jiǎn)稱CRU）1901年1月到2009年12月的高分辨率全球逐月格點(diǎn)數(shù)據(jù)集（CRU TS 3.1），其空間分辨率為 0.5°×0.5°。

1.4 研究方法

采用Dendro Clim 2002相關(guān)分析和響應(yīng)分析解釋樹輪徑向生長(zhǎng)對(duì)氣候的響應(yīng)[30]；采用區(qū)域平均序列的逐步回歸方法建立歷史溫度變化序列；采用國際年輪研究中常用的“逐一剔除法”進(jìn)行交叉檢驗(yàn)，從誤差縮減值RE、相關(guān)系數(shù)r、一階差相關(guān)系數(shù)rd、符號(hào)檢驗(yàn)值z(mì)、一階差符號(hào)檢驗(yàn)值z(mì)d和乘積平均數(shù)t等幾個(gè)方面對(duì)重建方程進(jìn)行交叉檢驗(yàn)。如果這些檢驗(yàn)統(tǒng)計(jì)量中的誤差縮減值或其它某幾個(gè)能通過檢驗(yàn)，則說明，該重建方程是穩(wěn)定的，由其重建出的溫度是可靠的。采用對(duì)稱延伸法消除子波變換效應(yīng)的Morlet子波變換和功變譜分析提取時(shí)間序列的變化準(zhǔn)周期。

2 結(jié)果與分析

2.1 樹輪δ13C序列對(duì)氣候的響應(yīng)及樹木生理學(xué)意義

研究表明，樹木年輪碳同位素序列主要反映區(qū)域性氣候變化而非局地氣候因子的影響[31-32]。相關(guān)分析表明，阿勒泰KLD樹輪δ13C序列與阿勒泰地區(qū)區(qū)域平均氣溫有較好的相關(guān)性，其中能夠代表夏季氣溫的7月平均氣溫與樹輪δ13C序列顯著正相關(guān)，相關(guān)系數(shù)高達(dá)0.608。阿爾泰山南坡樹輪δ13C序列與7月平均氣溫呈明顯的正相關(guān)關(guān)系，從生理角度講，阿爾泰山山區(qū)（采樣區(qū)）主要位于亞寒帶濕潤、半濕潤區(qū)，充沛的降水使得生長(zhǎng)在上樹線的樹木在生長(zhǎng)季時(shí)不會(huì)存在缺水的情況，相反，樹木的各種生命活動(dòng)都必須在一定的溫度范圍內(nèi)進(jìn)行，光合速率通過溫度影響光合酶的活性[33]，7月溫度偏高，較高的光合速率導(dǎo)致進(jìn)入葉片細(xì)胞內(nèi)的CO2濃度的降低，反映在δ13C值偏大。在生長(zhǎng)期（尤其是7月），阿爾泰山西伯利亞落葉松樹木碳同位素分餾的主要控制因子為光合速率（溫度）。也就是說，夏季平均氣溫是位于上樹線樹木δ13C變化的主要限制因子，由此可見，阿爾泰山西伯利亞落葉松樹輪δ13C序列與7月平均氣溫的關(guān)系是有其樹木生理學(xué)意義的。阿爾泰山北坡樹輪對(duì)氣候的響應(yīng)研究也表明，δ13C與7月平均氣溫顯著相關(guān)[34]。另外，Alexander等[35]發(fā)現(xiàn)位于西伯利亞東部的落葉松樹輪穩(wěn)定碳同位素與該地區(qū)6—7月的平均最高氣溫顯著正相關(guān)（r=0.46，p<0.001），Valerie等[36]研究發(fā)現(xiàn)，阿拉斯加白云杉樹輪δ13C與該地區(qū)夏季（5—8月）平均氣溫相關(guān)顯著，Polona等[37]研究表明，斯洛文尼亞挪威云杉樹輪δ13C值與生長(zhǎng)季內(nèi)溫度顯著正相關(guān)，Sheu等研究表明，臺(tái)灣冷杉樹輪穩(wěn)定碳同位素對(duì)5—10月溫度響應(yīng)較好[38]；Liu等[39]研究表明，賀蘭山油松樹輪δ13C序列與夏季（6—8月）平均氣溫相關(guān)顯著，這些研究結(jié)果均與阿爾泰山西伯利亞落葉松樹輪δ13C序列與7月平均氣溫顯著正相關(guān)的結(jié)果較為一致。

圖2 阿勒泰地區(qū)近50 a氣溫和降水的年變化和月變化（1962—2012）

2.2 阿勒泰地區(qū)過去160 a平均氣溫的重建及檢驗(yàn)

將采自阿勒泰科勒迭難布拉克的樹輪δ13C序列與阿勒泰地區(qū)7月平均氣溫（1962—2009年）進(jìn)行逐步回歸分析，得到轉(zhuǎn)換函數(shù)：

轉(zhuǎn)換函數(shù)的相關(guān)系數(shù)是0.608，方差解釋量為36.9%，調(diào)整自由度后的解釋方差26.929，遠(yuǎn)遠(yuǎn)超過0.000 1的極顯著水平，由該方程可以重建出1850—2009年阿爾泰7月平均氣溫（圖3）。

圖3 阿勒泰地區(qū)7月平均氣溫重建序列

方程檢驗(yàn)值見表1，7月平均氣溫重建序列的一階差符號(hào)檢驗(yàn)（zd)，符號(hào)檢驗(yàn)（z）均通過99%的信度檢驗(yàn)，說明該重建序列無論是高頻變化還是低頻變化，與實(shí)測(cè)序列一致性均較好；t值通過0.001的顯著性水平，RE值大于0.3，可見重建方程是較為穩(wěn)定可靠的。因此，利用該方程可較好的重建阿勒泰過去160年7月平均氣溫（圖4）。

表1 重建方程的統(tǒng)計(jì)檢驗(yàn)

圖4 重建溫度值（虛線）和實(shí)測(cè)值（實(shí)線）比較

2.3 近160 a氣溫變化特征分析

利用功率譜分析方法對(duì)阿勒泰地區(qū)夏季平均氣溫重建序列進(jìn)行周期分析，取最大滯后M=53，相當(dāng)于序列長(zhǎng)度的1/3，當(dāng)功率譜大于一定置信水平所對(duì)應(yīng)的臨界值，則認(rèn)為相應(yīng)周期顯著。分析發(fā)現(xiàn)，阿勒泰地區(qū)過去160 a的夏季氣溫以2 a的周期最為顯著，與“準(zhǔn)兩年脈動(dòng)（QBO）”十分接近，超過了0.05的顯著性水平。已有研究表明，QBO的影響存在于較大的范圍，其可能與海氣間相互耦合振蕩有關(guān)[40]?；趯?duì)稱延伸法的Morlet子波變換表明（圖5），阿爾泰地區(qū)過去160 a夏季平均氣溫周期是隨時(shí)間變化的，1900年以前，30 a的變化準(zhǔn)周期較為顯著，而1900年以后，20 a左右的變化準(zhǔn)周期明顯，1900年左右和1980年左右還存在著10 a的變化準(zhǔn)周期.10 a和20 a的變化準(zhǔn)周期可能對(duì)應(yīng)太陽黑子活動(dòng)11 a和22 a準(zhǔn)周期，說明太陽活動(dòng)對(duì)這一區(qū)域的溫度變化也有影響[41]。這些周期暗示了阿勒泰地區(qū)過去溫度變化除受當(dāng)?shù)貧夂蜃兓绊懲?，還受到更大范圍氣候變化的影響。

圖5 夏季溫度重建序列的MORLET小波變換圖

近160 a來，阿勒泰地區(qū)夏季溫度19世紀(jì)末有緩慢上升趨勢(shì)，20世紀(jì)初，夏季溫度整體偏高，40年代以后，夏季溫度明顯下降，直到60年代，60年代成為近160 a夏季溫度最低的時(shí)期，60年代至今，阿爾泰山南坡夏季氣溫經(jīng)歷了160 a來最長(zhǎng)和最為強(qiáng)烈的升溫時(shí)期。雖然有氣象觀測(cè)記錄以來，阿勒泰地區(qū)夏季氣溫升溫明顯，但是20世紀(jì)初夏季溫度要高于近30 a平均氣溫（圖6）。有研究表明，新疆的冰川在20世紀(jì)初退縮相對(duì)加快。這種退縮一直延續(xù)至本世紀(jì)70年代。這種退縮標(biāo)志著新疆偏暖，這種變暖從l9世紀(jì)90年代中期開始呈上升趨勢(shì)，至20世紀(jì)40年代又逐漸下降，其間出現(xiàn)一些小的峰谷[42]。本研究所獲得的過去160 a阿勒泰地區(qū)夏季溫度變化與新疆冰川記錄完全一致。

圖6 阿勒泰地區(qū)夏季溫度的年代際變化

將重建的阿勒泰地區(qū)夏季溫度序列與大范圍CRU7月氣溫空間相關(guān)分析表明（圖7），重建的阿勒泰地區(qū)夏季平均氣溫序列不僅能代表整個(gè)阿爾泰山夏季溫度變化，也能夠較好的代表包括北疆、中亞部分地區(qū)和俄羅斯南部大范圍的夏季溫度變化。

圖7 阿勒泰地區(qū)7月平均氣溫與CRU 3.10的7月平均氣溫空間相關(guān)（1960—2009）

3 結(jié)論

阿爾泰山南坡西伯利亞落葉松森林上線的樹輪δ13C序列對(duì)阿勒泰地區(qū)夏季氣溫響應(yīng)較好，與7月平均氣溫明顯正相關(guān)，相關(guān)系數(shù)達(dá)到0.608（p<0.000 1）；利用樹輪δ13C序列可較好的重建阿勒泰地區(qū)過去160 a的夏季氣溫。

近160 a來，阿勒泰地區(qū)夏季溫度19世紀(jì)末有緩慢上升趨勢(shì)，20世紀(jì)初，夏季溫度整體偏高，40年代以后，夏季溫度明顯下降，直到60年代，60年代成為近160年夏季溫度最低的時(shí)期，60年代至今，阿爾泰山南坡夏季氣溫經(jīng)歷了160 a來最長(zhǎng)和最為強(qiáng)烈的升溫時(shí)期。雖然有氣象觀測(cè)記錄以來，阿勒泰地區(qū)夏季氣溫升溫明顯，但是20世紀(jì)初夏季溫度要高于近30 a平均氣溫。

阿勒泰地區(qū)過去160 a的夏季氣溫以2 a的周期最為顯著，1900年以前，30 a的變化準(zhǔn)周期較為顯著，而1900年以后，20 a左右的變化準(zhǔn)周期明顯，1900年左右和1980年左右還存在著10 a的變化準(zhǔn)周期，這些周期暗示了阿勒泰地區(qū)過去溫度變化除受當(dāng)?shù)貧夂蜃兓绊懲猓€受到更大范圍氣候變化的影響，進(jìn)一步證明了氣候變化的全球性，重建的阿勒泰地區(qū)夏季平均氣溫序列不僅能代表整個(gè)阿爾泰山夏季溫度變化，也能夠較好的代表包括北疆、中亞部分地區(qū)和俄羅斯南部大范圍的夏季溫度變化。

[1]Schwikowski M.，Eichler A.，Kalugin I.，et al.Past climate variability in the Altai[J].PAGES news 2009，17（1）：44-45.

[2]O vtchinnikov D.，Adamenko M.,Panushkina I..A 1105-yeartree-ring chronology in Altairegion and its application for reconstruction of summer temperatures[J].Geolines，2000，11：121-122.

[3]李江風(fēng).新疆年輪氣候年輪水文研究[M].北京：氣象出版社，1989.

[4]張同文，魏文壽，袁玉江，等.阿勒泰西部樹輪年表特征分析[J].中國沙漠，2007，27（6）：1040-1047.

[5]Wang，L.Duan，J.P.，Chen,J.，et al.Temperature reconstruction from tree-ring maximum density of Balfour spruce in eastern Tibet,China[J].International Journal of Climatology.2010，30：972-979.

[6]尚華明，張瑞波，魏文壽，等.胡楊樹輪δ13C記錄的艾比湖地區(qū)夏季最高溫度變化 [J].沙漠與綠洲氣象，2013，7（5）：7-13.

[7]Liu Y.，Wang R.Y.,Leavitt S.W.，et al.Individual and pooled tree-ring stable-carbon isotope series in Chinese pine from the Nan Wutai region,China:Common signal and climate relationships[J],Chemical Geology，2012，330-331:17-26.

[8]張瑞波，袁玉江，魏文壽，等.利用樹輪灰度重建烏孫山被批 4—5 月平均最低溫度[J].中國沙漠，2008，28（5）：848-854.

[9]Panyushkina I.P.，Ovtchinnikov D.V.,Adamenko M.F..Mixed response of decadal variability in Larch tree-ring chronologies from upper tree-lines of the Russian Altai[J].Tree-ring Research，2005，61（1）：33-42.

[10]胡義成，袁玉江，魏文壽，等.用樹木年輪重建阿爾泰東部 6—7月平均溫度序列 [J].中國沙漠，2012，32（4）：1003-1009.

[11]胡義成，魏文壽，袁玉江，等.新疆阿勒泰中東部樹輪寬度年表特征及其氣候響應(yīng)[J].沙漠與綠洲氣象，2011，5（1）：6-12.

[12]李漠巖，袁玉江，魏文壽，等.阿勒泰西部樹輪寬度年表特征及其氣候響應(yīng)[J].沙漠與綠洲氣象，2011，5（6）：16-21.

[13]張同文，袁玉江，魏文壽，等.使用樹輪資料重建阿勒泰西部年積雪深度≥0cm日數(shù)的長(zhǎng)期變化[J].沙漠與綠洲氣象，2010，4（3）：6-11.

[14]Sidorova O.V.，Saurer M.,Myglan V.S.，et al.A multiproxy approach for revealing recent climatic changes in the Russian Altai[J].Clim Dyn，2012，38:175-188.

[15]Chen F.,Yuan Y.J.，Wei W.S.，et al.Climatic response of ring width and maximum latewood density of Larix sibirica in the Altay Mountains,receals recent warming trends[J].Annals Forest Science，2012，69:723-733.

[16]Michele E M，Kingston J D，Marino B D.Carbon isotope evidence for emergence of C4 plants in the Neogene from Pakistan and Kenya[J].Nature，1994，367，162-165.

[17]Wilson A T,Grinsted M J.13C/12C in cellulose and lignin paleathermometer.Nature，1977，265:133-135.

[18]Leavitt S W，Long A.Seasonal stable-carbon isotope variability in tree ringspossible paleoenvironmental signals [J].Chemical Geology （Isotope Geoscience Science），1991，87:59-70.

[19]張瑞波，袁玉江，魏文壽，等.西伯利亞落葉松樹輪穩(wěn)定碳同位素對(duì)氣候的響應(yīng) [J].干旱區(qū)研究，2012，29（2）：328-334.

[20]Holmes R L.Computer-assisted quality control in treering dating and measurement[J].Tree-Ring Bulletin，1983，43:69-75.

[21]Cook E R.A time series analysis approach to tree-ring atandardization[D].Tucson:The University of Arizona，1985.

[22]Borella S.，Leuenberger M.，Saurer M.,et al.Reducing uncertainties in δ13C analysis of tree rings:Pooling,milling，and cellulose extraction[J].Journal of Geophysical Research，1998，103：19519-19526.

[23]Leavitt S.W.，Long A.Sampling strategy for stable carbon isotope analysis of tree rings in pine [J].Nature，1984，311：145-147.

[24]侯愛敏，彭少麟.周國逸.樹木年輪δ13C含量幼齡效應(yīng)的定量化探討[J].生態(tài)學(xué)報(bào)，2001，21（3）：430-433.

[25]Coplen T B.Discontinuance of SMOW and PDB[J].Nature，1995，375：285.

[26]Stuiver M.Atmospheric carbon dioxide and carbon reservoir Change[J].Science，1978，199：253-258.

[27]Stuiver M,Braziunas T.Tree cellulose 13C/12C isotopic ratios and climatic change[J].Nature，1987，328：58-60.

[28]McCarroll D，Loader N J.Stable isotopes in tree rings[J].Quaternary Science Reviews,2004，23（7）：771-801.

[29]Mann H.B.Non-parametric testagainsttrend.Econometrika，1945，13：245-259

[30]Biondi，F(xiàn).andWaikul,K.2004.DENDROCLIM 2002：a C++program for statistical calibration of climate signals in tree-ring chronologies[J].Comput.Geosci.-UK 30，303-311.

[31]Anderu L.，PlanellsO.，GutierrezE.，etal.Climatic significance of tree-ring width and δ13C in a Spanish pine foreat network[J].Tellus 2008，60B，771-781.

[32]Mazany T.，Lerman J.C.，Long A..Carbon-13 in tree-ring cellulose as an indicator of past climates[J].Nature 1980，287，432-435.

[33]Beerling D J.Predicting leaf gas exchange and δ13C response to the past 30 000 years of global environmental change[J].New Phytologist，1994，128：425-433.

[34]Sidorova O.V.，Saurer M.,Myglan V.S.，et al.A multiproxy approach for revealing recent climatic changes in Russian Altai[J].Climate Dynamics，2012，38：175-188.

[35]Alexander V K,Kerstin S T,Anatolii N,et al.Climate signals in tree-ring width,density and δ13C from larches in Eastern Siberia （Russia）[J].Chemical Geology，2008，252：31-41.

[36]Valerie A B,Glenn P J,Bruce P F,et al.Reconstruction of Summer temperatures in interior Alaska from tree-ring proxies:Evidence for changing synoptic climate regimes[J].Climatic Change，2004，63：91-120.

[37]Polona H，Tom L.Stable carbon isotopes in Norway Spruce（Picea abies（L.）Karst.） tree rings at two sites in Slovenia[J].Zbornik Gozdarstva in Lesarstva，2009，88：43-52.

[38]Sheu D D,Kou P,Chiu C H,et al.Variability of tree ring δ13C in Taiwan fir:Growth effect and response to May-October temperatures[J].Geochimica et Cosmochimica Acta，1996，60（1）：171-177.

[39]Liu Y.，Ma L.M,Cai Q.F，et al.Reconstruction of summer temperature （June-August） at Mt.Helan,China,from tree-ring stable carbon isotope values since AD 1890[J].Science in China:Series D，2002，45（12）：1127-1136.

[40]錢維宏，朱亞芬，葉謙.赤道東太平洋海溫異常的年際和年代際變率[J].科學(xué)通報(bào)，1998，43（10）：1098-l102.

[41]Rasnusson E.M.，Wang,X.L.，Ropelewski，C.F.The biennial component of ENSO variability[J].Journal of Marine Systems 1990，1：71-96.

[42]胡文康.20世紀(jì)塔克拉瑪干沙漠環(huán)境及其變遷[J].干旱區(qū)研究，1992，9（4）：1-9.