張永

（中國(guó)科學(xué)院地理科學(xué)與資源研究所，陸地表層格局與模擬重點(diǎn)實(shí)驗(yàn)室，北京 100101）

經(jīng)驗(yàn)?zāi)B(tài)分解方法對(duì)祁連圓柏年輪寬度序列標(biāo)準(zhǔn)化初探

張永

（中國(guó)科學(xué)院地理科學(xué)與資源研究所，陸地表層格局與模擬重點(diǎn)實(shí)驗(yàn)室，北京 100101）

利用經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）方法對(duì)祁連圓柏樹輪寬度序列進(jìn)行標(biāo)準(zhǔn)化處理的初步探討，分析該方法在圓柏寬度序列生長(zhǎng)量訂正中的應(yīng)用潛力。利用EMD方法對(duì)祁連山地區(qū)的65個(gè)圓柏樣芯寬度序列的逐級(jí)分解，較好地展示了各樣芯由高頻至低頻的自然變化特征，其中最低頻的趨勢(shì)項(xiàng)雖包含了線性、負(fù)指數(shù)、類樣條函數(shù)等多種形式，但末端基本趨于平穩(wěn)，與樹木的生物學(xué)遺傳特性較為吻合，基本反映了與樹木遺傳相關(guān)的生物學(xué)趨勢(shì)。與傳統(tǒng)線性或負(fù)指數(shù)函數(shù)方法對(duì)祁連圓柏序列的標(biāo)準(zhǔn)化過程相比，傳統(tǒng)處理中常常出現(xiàn)髓心端上翹或樹皮端異常增大的問題，通過新的EMD方法的訂正得到了改進(jìn)，進(jìn)一步分析了EMD方法在祁連圓柏輪寬序列生長(zhǎng)量訂正中的優(yōu)勢(shì)及不足。兩種方法得到的祁連圓柏年表在整個(gè)時(shí)段上高頻變化極為一致，僅在接近髓心端的100a內(nèi)傳統(tǒng)方法處理得到的指數(shù)值偏低，方差變化較大；而EMD處理得到的指數(shù)年表整個(gè)時(shí)段方差相對(duì)穩(wěn)定。

樹輪寬度序列；標(biāo)準(zhǔn)化方法；經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）

氣候變化是影響世界環(huán)境、人類健康與福利和全球經(jīng)濟(jì)持續(xù)性發(fā)展的重要因素，了解氣候變化的規(guī)律及提高對(duì)未來氣候變化的預(yù)估水平成為現(xiàn)今地球科學(xué)最重要的研究課題之一。作為過去氣候環(huán)境變化的天然信息庫(kù)，樹木年輪以其定年準(zhǔn)確、分辨率高、連續(xù)性強(qiáng)、空間分布廣、環(huán)境變化指示意義明確且可精確定量化等優(yōu)勢(shì)成為研究過去氣候變化的一種重要代用指標(biāo)[1-10]。

然而，樹木年輪的形成與變異除了取決于氣候和其它環(huán)境因子以外，還受樹木本身的遺傳因子控制。樹輪的生長(zhǎng)可以用線性聚類模式[11]表示如下：

Rt=At+Ct+△D1t+△D2t+Et（1）

其中Rt表示年輪寬度；At表示生長(zhǎng)趨勢(shì)；Ct表示氣候信息；△D1t表示內(nèi)部干擾；△D2t表示外部干擾；Et表示其它。

而我們主要研究環(huán)境因子中的氣候因子對(duì)樹木生長(zhǎng)的影響，因此，對(duì)實(shí)際測(cè)量得到的年輪寬度值，在分析其能否反映出某些氣候要素的變化之前，應(yīng)設(shè)法消除樹齡對(duì)年輪寬度的這種緩慢影響，即樹木本身支配其生長(zhǎng)的因素（遺傳因子）應(yīng)被剔除，從而得到僅由氣候因子影響所決定的年輪序列變化量，該過程被稱為生長(zhǎng)量訂正。然而，樹齡的影響不像氣候因子那樣表現(xiàn)在逐年的輪寬相對(duì)增量上，而主要表現(xiàn)為緩慢的趨勢(shì)性變化。一般來說，多數(shù)樹木的年輪寬度均有下述變化趨勢(shì)：幼齡期年輪寬度比較窄，后隨樹齡的進(jìn)一步增加和樹干直徑增大的幾何效應(yīng)，到達(dá)樹木生長(zhǎng)最旺盛的時(shí)期時(shí)，平均的年輪寬度常為整個(gè)生長(zhǎng)期內(nèi)的極大值，此后，樹木的年輪寬度隨樹齡的增加而趨于減小，減小的速度開始較快，后來越來越慢，到一定時(shí)期就大體趨于平穩(wěn)，被稱為樹木生長(zhǎng)趨勢(shì)[12-13]。

如何去除樹木生長(zhǎng)趨勢(shì)的影響，對(duì)樹輪資料進(jìn)行處理時(shí)存在的去趨勢(shì)偏差、低頻信號(hào)的提取和保留等也是樹木年輪氣候?qū)W研究中的一個(gè)重要部分[12]，這限制和影響著樹輪資料在研究氣候變化時(shí)的準(zhǔn)確性。

1 傳統(tǒng)的樹輪生長(zhǎng)量訂正方法

精確地確定出樹木的生長(zhǎng)趨勢(shì)并作出準(zhǔn)確的生長(zhǎng)量訂正是一項(xiàng)十分困難的事情，這是由于隨著樹種和環(huán)境的不同，其樹木的群體生長(zhǎng)趨勢(shì)也有較大差異。此外，生長(zhǎng)量訂正的另一目的是要把樹輪寬度資料整理成可比較的形式。長(zhǎng)期以來，為了達(dá)到上述兩個(gè)目的，樹木年輪氣候?qū)W者提出了許多生長(zhǎng)量訂正的方法和手段[4]，設(shè)法消除這種影響，被稱為年輪寬度序列標(biāo)準(zhǔn)化[13-14]。根據(jù)生長(zhǎng)趨勢(shì)曲線估算的思路和方法大體上可歸為3類：樹輪解析估計(jì)法、大樣本估算法、多種函數(shù)擬合法[13-14]。

隨著樹輪學(xué)者對(duì)樹輪去趨勢(shì)理論認(rèn)識(shí)的提高及經(jīng)驗(yàn)的積累，多種函數(shù)逐步被嘗試用來表達(dá)樹木年輪隨樹齡變化的生長(zhǎng)趨勢(shì)曲線。其中主要包括線性函數(shù)、雙曲線函數(shù)、指數(shù)函數(shù)、多項(xiàng)式函數(shù)及樣條函數(shù)等。但每種函數(shù)也都有自己的局限性，其中目前以線性或負(fù)指數(shù)函數(shù)應(yīng)用最為廣泛，通常也被認(rèn)為是比較適用于干旱、半干旱區(qū)樹種生長(zhǎng)趨勢(shì)擬合的方法[15-16]。特殊情況會(huì)考慮樣條函數(shù)，樣條函數(shù)[15-16]一般用來擬合較濕潤(rùn)環(huán)境下樹木生長(zhǎng)的持續(xù)性以及由樹木相互競(jìng)爭(zhēng)產(chǎn)生的非同步擾動(dòng)。應(yīng)用中發(fā)現(xiàn)負(fù)指數(shù)和直線擬合在樹輪曲線兩端的擬合常出現(xiàn)較大偏差，樣條函數(shù)一定程度上可以適應(yīng)環(huán)境干擾信號(hào)的擬合，但當(dāng)樣條函數(shù)的步長(zhǎng)選擇不當(dāng)時(shí)，雖然有可能看起來與原序列的擬合效果甚佳，但此時(shí)可能不僅消除了樹木自身的生長(zhǎng)趨勢(shì)，也將受長(zhǎng)期氣候影響所造成的樹木年輪寬度緩慢變化濾掉，造成“信號(hào)”或“噪聲”重疊的低頻部分的“過度擬合”，很難判定最終得出的氣候信息是否完整可靠，所以在選用其進(jìn)行生長(zhǎng)量訂正的過程中有一定爭(zhēng)議。

近年來，研究中常有報(bào)道，即使采樣點(diǎn)的樹木生長(zhǎng)對(duì)溫度敏感，卻依然很難重建出中世紀(jì)暖期（MWP）、小冰期（LIA）等近千年中典型氣候事件，其主要原因就是去趨勢(shì)的問題，如Mann等[8]利用樹輪重建的北半球過去近千年來的溫度距平序列，因中世紀(jì)暖期不明顯而受到較大爭(zhēng)議。就此問題，國(guó)外學(xué)者Briffa[18-19]，Esper[20]開展了關(guān)于樹輪低頻信息提取方面的研究，提出了區(qū)域曲線標(biāo)準(zhǔn)化（Regional Curve Standardization,RCS）方法，但局限性在于該方法必須保證對(duì)樣本總體生長(zhǎng)曲線（RC）有較好估計(jì)的情況下方可使用，否則會(huì)有較大偏差，另外對(duì)于單個(gè)樣芯真實(shí)生理年齡的估計(jì)方法也存在較多爭(zhēng)議。Melvinand Briffa[21]建立了“signal-free”方法用于解決十年至百年尺度的異常變化對(duì)傳統(tǒng)的樹木生長(zhǎng)擬合曲線的影響，但也同樣存在低頻變化的缺失以及使用不同擬合方法導(dǎo)致的不確定性的問題；國(guó)內(nèi)的樹輪學(xué)者還提出了總體生長(zhǎng)趨勢(shì)曲線擬合法[22]，嘗試進(jìn)一步完善去趨勢(shì)過程中低頻信息的保留以及樹木生長(zhǎng)總體趨勢(shì)的估計(jì)，但目前仍處于初步研究階段，未推廣應(yīng)用。

2 經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）方法與圓柏生長(zhǎng)量訂正

2.1 EMD方法介紹

經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical mode decomposition, EMD)[14]，是20世紀(jì)末由Huang提出來的一個(gè)信號(hào)分析方法，EMD分解的主要思想是利用波動(dòng)上、下包絡(luò)的平均值去確定“瞬時(shí)平衡位置”，進(jìn)而分解IMF分量，每次提取都將高頻信號(hào)分離出，剩下低頻信號(hào)，本質(zhì)上是一個(gè)過濾器。按照它所確定的程序，通過逐步篩選，把一個(gè)混合信號(hào)序列（或其導(dǎo)數(shù)，視所需的分解精度而定）分解成有限個(gè)更為簡(jiǎn)單的反映某一個(gè)物理信號(hào)在窄波段上隨時(shí)間變化的IMF分量。這些分量保持了原來信號(hào)的自然振蕩特性，并且是相互獨(dú)立的。當(dāng)所有的IMF分量被提取完畢之后，剩余項(xiàng)被稱為原始信號(hào)的趨勢(shì)項(xiàng)。對(duì)于一個(gè)給定的時(shí)間序列x（t），可以通過下面的基本步驟完成它的分解：

（1）找出x（t）所有的極大值點(diǎn)并將其用三次樣條函數(shù)擬合成原始數(shù)據(jù)的上包絡(luò)線；類似的找出極小值點(diǎn)擬合成下包絡(luò)線；并分別記為u（t）、v（t）和m（t），其中m（t）滿足：

（2）計(jì)算x（t）和m（t）之差，并記為：h1（t）=x（t）-m（t）。

如果m（t）等于零，或者滿足設(shè)定的精度要求，那么，h1（t）就給出了第一個(gè)固有模態(tài)函數(shù)，并記為IMF1；否則，重復(fù)步驟（1）和（2）。

（3）用殘差r1（t）=x（t）-h1（t）代替原始信號(hào)x（t），并重復(fù)上面的篩選過程，可以得到第二個(gè)IMF2。

（4）重復(fù)步驟（3），依次得到ri（t）和hi+1（t）（i= 2，3，…，n）。

如果rn-1（t）的極值點(diǎn)不超過2個(gè)，則整個(gè)分解過程結(jié)束。EMD分解最終把原序列分解成有限個(gè)（n個(gè)）IMF分量，不同層次得IMF分量都可能對(duì)應(yīng)某一物理背景。因此，對(duì)于一個(gè)混沌序列，找出其主要特征量即對(duì)應(yīng)主要物理背景的特征層次，對(duì)混沌序列的預(yù)測(cè)具有重要意義。而原始序列就可以由這些IMF分量及均值之和表示：

分解后所得到每個(gè)IMF分量具有如下特征：（1）從全局特性看，極值點(diǎn)數(shù)必須與過零點(diǎn)數(shù)一致或者至多相差一個(gè)；（2）在某一個(gè)局部點(diǎn)，極大值包絡(luò)和極小值包絡(luò)在該點(diǎn)的算術(shù)平均和為零。為了分辨分解結(jié)果中可能的虛假分量，可以通過計(jì)算IMF和細(xì)節(jié)分量與原序列的相關(guān)系數(shù)，得到與原序列顯著相關(guān)的IMF和細(xì)節(jié)分量。

2.2 EMD對(duì)祁連圓柏輪寬序列的分解

由上文樹輪寬度序列的聚類模式表達(dá)式可以看出，樹輪寬度序列正是一個(gè)多信號(hào)混合的混沌序列，為了分析其中最主要的氣候因子對(duì)輪寬變化的影響，借助于EMD分解的方法，將各信號(hào)逐級(jí)進(jìn)行分離，去除與年齡相關(guān)的生物趨勢(shì)及其他非氣候噪音是一種可行的思路。利用EMD運(yùn)算程序，對(duì)祁連山孔崗木祁連圓柏[24]（Kgm02）02a和09a樣芯寬度測(cè)量序列進(jìn)行逐級(jí)分解（圖1和圖2）?？梢钥闯鰞蓚€(gè)樹芯的寬度序列均包含了7個(gè)本征函數(shù)分量，每個(gè)IMF分量對(duì)應(yīng)于一個(gè)窄波段的信號(hào)，反映樹輪序列記錄的該尺度信號(hào)的自然波動(dòng)過程，IMF7均為原始變量的最低頻趨勢(shì)分量，各分量之和可還原為原始變量序列（藍(lán)線）。

圖1所示Kgm0202a芯的7條IMF分量反映了該樣芯由高頻至低頻不同準(zhǔn)周期尺度上的瞬時(shí)波動(dòng)特征，而本研究中主要關(guān)注EMD分解所得趨勢(shì)項(xiàng)對(duì)樹輪寬度序列生長(zhǎng)趨勢(shì)的擬合效果，趨勢(shì)項(xiàng)IMF7分量1692—2000年間由高至低，逐漸趨于平穩(wěn)，反映的生長(zhǎng)趨勢(shì)近似于負(fù)指數(shù)函數(shù)擬合曲線。

圖2展示的Kgm0209a樣芯的IMF7趨勢(shì)分量變化過程與前者有較大差異，該序列的IMF7分量1651—1757年逐步加速至峰值，1757年之后開始逐步減慢，至19世紀(jì)末趨于平穩(wěn)，與樹木自幼齡期到達(dá)壯年的自然生長(zhǎng)過程極為吻合。該曲線近似于樣條函數(shù)擬合曲線，但不同之處在于樣條函數(shù)的步長(zhǎng)選取較為主觀，研究中如果選取不當(dāng)則會(huì)導(dǎo)致過度擬合的情形，而EMD分解則是遵循序列本身的自然震蕩進(jìn)行逐級(jí)分解，不存在人為假定和預(yù)先設(shè)定函數(shù)形式的問題。

Kgm0202a和Kgm0209a雖然均是按照同一原理進(jìn)行分解，可因?yàn)樾蛄斜旧戆恼鹗庍^程不同，也就是說兩條序列受起始年靠近髓心位置不同及局地小生境的不同，得到的序列生長(zhǎng)趨勢(shì)項(xiàng)也表現(xiàn)了較大差異。因此，靠一個(gè)模式或一個(gè)一成不變的公式是很難滿足完全剔除掉所有樣芯的生物學(xué)效應(yīng)的。在實(shí)際應(yīng)用中根據(jù)所采集樣芯的生態(tài)環(huán)境及樹種，參照多序列的綜合分析，才是保證較合理剔除各樣芯生長(zhǎng)趨勢(shì)的關(guān)鍵。

進(jìn)一步對(duì)Kgm02樣點(diǎn)的65個(gè)樣芯分別進(jìn)行了EMD分解，各序列的最低頻分量趨勢(shì)項(xiàng)見圖3。由圖3可以看出，EMD分解能較好地適應(yīng)樹輪序列趨勢(shì)項(xiàng)多樣不一的情況，針對(duì)各序列不同的生境及生長(zhǎng)情況逐步分離各序列所包含的多尺度震蕩分量，同時(shí)避免了選擇樣條函數(shù)擬合步長(zhǎng)選取的主觀性，有助于較準(zhǔn)確的獲取趨勢(shì)項(xiàng)。

圖1 Kgm0202a寬度序列EMD分解各分量

圖2 Kgm0209a寬度序列EMD分解各分量

圖3 Kgm02各樣芯寬度序列EMD分解的趨勢(shì)項(xiàng)分量

趨勢(shì)項(xiàng)包含了傳統(tǒng)使用的線性函數(shù)、負(fù)指數(shù)函數(shù)、類樣條函數(shù)等各種形式，但總體而言，趨勢(shì)項(xiàng)在末端均表現(xiàn)出趨于平穩(wěn)的特點(diǎn)，該特點(diǎn)與樹木生長(zhǎng)晚期生理學(xué)效應(yīng)減弱并趨于平穩(wěn)的特征較為吻合，進(jìn)一步證實(shí)了EMD在趨勢(shì)項(xiàng)獲取中的可行性。2.3EMD與傳統(tǒng)方法對(duì)圓柏生長(zhǎng)量訂正的比較

為了檢驗(yàn)EMD方法在祁連圓柏生長(zhǎng)量訂正中的實(shí)際效果，遵循常規(guī)的樹輪曲線標(biāo)準(zhǔn)化訂正原理（原始曲線除以趨勢(shì)項(xiàng)），分別按照傳統(tǒng)擬合法和EMD方法對(duì)各序列進(jìn)行趨勢(shì)項(xiàng)分離，并對(duì)傳統(tǒng)方法和EMD方法獲得的65個(gè)輪寬序列去除生長(zhǎng)趨勢(shì)后得到的指數(shù)序列進(jìn)行逐一對(duì)比。圖4展示了四組具有代表性的序列細(xì)節(jié)對(duì)比圖，上下兩圖分別為一組，共4組（a，b；c，d；e，f；g，h）。每組上圖為輪寬曲線及傳統(tǒng)線擬合（線性或負(fù)指數(shù)）及EMD分解的趨勢(shì)項(xiàng)，下圖為對(duì)應(yīng)的去趨勢(shì)后指數(shù)序列。

其中最為突出的差異之一表現(xiàn)為兩種方法得到的指數(shù)序列在接近髓心端的差異。以祁連圓柏代表性樣芯Kgm0211a（圖4a和4b）為例可以看出，傳統(tǒng)方法擬合的趨勢(shì)項(xiàng)是斜率為負(fù)的直線（圖4a紅線），EMD分解得到的趨勢(shì)項(xiàng)（圖4a藍(lán)線）自1600年開始先慢后快，至1850年左右趨于平穩(wěn)。相對(duì)應(yīng)兩種方法得到的指數(shù)序列在髓心端（最后50a）差異較大，傳統(tǒng)方法得到的序列值（圖4b紅線）高于EMD方法得到的指數(shù)值（圖4b藍(lán)線）大約0.2個(gè)標(biāo)準(zhǔn)差，導(dǎo)致髓心端異常上翹。

相對(duì)于髓心端的異常，兩種方法對(duì)部分序列樹皮端的訂正也存在一定的差異。祁連圓柏代表性序列為Kgm0219b（圖4c黑線），兩種方法對(duì)原始寬度序列的擬合線（圖4c紅線和藍(lán)線）較為相近，僅在末端直線擬合表現(xiàn)了持續(xù)性的下降，而EMD趨勢(shì)項(xiàng)趨于水平值。相對(duì)應(yīng)得到的指數(shù)序列前者在后5a出現(xiàn)了整個(gè)時(shí)段的異常極大值，而后者則較好地避免了因趨勢(shì)線的誤差導(dǎo)致的偏差。

結(jié)合以上兩點(diǎn)不同，還有一類序列在兩種方法的對(duì)比中，整個(gè)時(shí)段會(huì)表現(xiàn)為波動(dòng)性偏差。以圓柏Kgm0222b（圖4e黑線）為例可以看出，傳統(tǒng)方法的擬合線仍為直線（圖4e紅線），而EMD得到的趨勢(shì)項(xiàng)自序列始端（1787年）生長(zhǎng)量的變化過程為由慢至快加速，1865年左右達(dá)到峰值，后逐漸減速至1990年趨于水平（圖4e藍(lán)線）。相對(duì)應(yīng)的兩條指數(shù)序列偏差可分為三段：1785—1830年（傳統(tǒng)方法的指數(shù)序列偏低）、1830—1935年（傳統(tǒng)方法的指數(shù)序列偏高）和1935—2000年（傳統(tǒng)方法的指數(shù)序列偏低）。而這些偏差分析可知，正是由于傳統(tǒng)擬合線對(duì)樹木生物學(xué)效應(yīng)低估的結(jié)果，EMD方法得到的指數(shù)序列則避免了這個(gè)偏差。

祁連圓柏Kgm0202a（圖4g，4h）則代表了兩種方法得到的較為一致的指數(shù)序列的例子。由圖4g可以看出，傳統(tǒng)方法擬合得到的趨勢(shì)線與EMD方法分解得到的趨勢(shì)項(xiàng)一致，相對(duì)應(yīng)訂正后的指數(shù)序列也幾乎完全重合，高低頻變化完全一致。對(duì)本文所選樣本祁連圓柏的統(tǒng)計(jì)分析表明，當(dāng)輪寬序列本身所包含生長(zhǎng)趨勢(shì)項(xiàng)較為單調(diào)，且生長(zhǎng)過程中不存在異常小生境擾動(dòng)時(shí)，兩種方法得到的結(jié)果也會(huì)相對(duì)比較一致。

雖然整體上EMD對(duì)大部分樣芯的訂正都表現(xiàn)出了較大優(yōu)勢(shì)，但對(duì)祁連圓柏65個(gè)樣芯的分析中也發(fā)現(xiàn)個(gè)別樣芯由EMD分解得到的趨勢(shì)項(xiàng)也會(huì)出現(xiàn)失效的情況。本文分析發(fā)現(xiàn)Kgm0215a和Kgm0218c兩個(gè)樣芯的EMD分解趨勢(shì)項(xiàng)在生長(zhǎng)量訂正中需要特殊處理。為了更清楚地看出EMD方法與傳統(tǒng)的負(fù)指數(shù)和線性函數(shù)擬合方法對(duì)多復(fù)本量樹輪曲線訂正的差異，本文還對(duì)比了兩種方法最終建立的Kgm02研究點(diǎn)（65芯）的樹輪年表，結(jié)果如圖5。

由圖5可以看出，兩條序列整體所反映的高頻信息極為一致，雖然由上面的分析可知少部分樣芯在傳統(tǒng)方法處理中會(huì)出現(xiàn)樹皮端的變形，但所有序列平均之后樹皮端變形被削弱，故兩種方法得到的該樣點(diǎn)年表在樹皮端基本一致；然而，因?yàn)閭鹘y(tǒng)方法對(duì)大部分樣芯接近髓心端的部分生物學(xué)效應(yīng)剔除效果較差，因此兩種方法得到的所有序列平均指數(shù)在髓心端的100a內(nèi)出現(xiàn)了較大差異。EMD方法的年表在1430—1500年略高于傳統(tǒng)方法得到的指數(shù)值，可能反映了該時(shí)期樹木生長(zhǎng)環(huán)境較為適宜的氣候條件，然而，具體的氣候指示意義需進(jìn)一步分析。

綜上所述，EMD方法與傳統(tǒng)方法對(duì)祁連圓柏生長(zhǎng)量訂正的改進(jìn)主要體現(xiàn)在序列的前后兩端，EMD方法分解得到的趨勢(shì)項(xiàng)與樹木生長(zhǎng)自幼齡期開始慢—快—慢的生長(zhǎng)特征較為吻合，因此對(duì)于序列兩端生長(zhǎng)量的訂正極為適應(yīng)，效果遠(yuǎn)高于線性或負(fù)指數(shù)函數(shù)擬合。在實(shí)際應(yīng)用中可重點(diǎn)分析去趨勢(shì)序列前后兩端方差的異常變化，進(jìn)而評(píng)估樹輪生長(zhǎng)量的訂正效果。然而，對(duì)于一種新方法的應(yīng)用，仍然需要更大量樣本實(shí)驗(yàn)的對(duì)比驗(yàn)證，尤其要根據(jù)不同的樹種及不同的采樣環(huán)境進(jìn)行深入分析，確定最適宜的生長(zhǎng)量訂正方法和手段。

3 結(jié)論

圖4 Kgm02各樣芯寬度序列傳統(tǒng)方法擬合（紅線）和EMD分解趨勢(shì)分量擬合（藍(lán)線）對(duì)比

圖5 Kgm02樣點(diǎn)傳統(tǒng)年表與EMD年表對(duì)比

EMD是對(duì)一個(gè)時(shí)間序列多尺度變化特征全面分析的新方法，可以將時(shí)間序列包含的不同尺度的自然振蕩和趨勢(shì)逐級(jí)分離，有利于更好地認(rèn)識(shí)目標(biāo)序列的多尺度變化特征和總體趨勢(shì)。本文針對(duì)祁連圓柏寬度序列生長(zhǎng)量的訂正進(jìn)行了傳統(tǒng)方法和EMD分析手段的全面對(duì)比分析，初步探討EMD分解趨勢(shì)項(xiàng)在氣候變化研究中對(duì)生長(zhǎng)量訂正的改進(jìn)效果及可應(yīng)用性。分析初步表明，與傳統(tǒng)方法相比，EMD方法在生長(zhǎng)量訂正的應(yīng)用在序列髓心端和樹皮端表現(xiàn)出最為明顯的提高。傳統(tǒng)線性或負(fù)指數(shù)的擬合在樹皮端會(huì)引起少數(shù)樣芯序列5～10a的異常極大值；而在髓心端對(duì)大部分樣芯先慢后快再慢至平穩(wěn)的生長(zhǎng)趨勢(shì)的反映較差，因此在髓心端大部分樣芯生長(zhǎng)量訂正失效，出現(xiàn)接近髓心端100a內(nèi)異常上翹或降低的現(xiàn)象，因此在研究中樹輪年表的后50～100a常常出現(xiàn)較大的方差而被遺棄，EMD方法在一定程度上緩解了這一問題。

總體而言，EMD方法在祁連圓柏寬度序列生長(zhǎng)量訂正中的適用性得到了初步證實(shí)，為生長(zhǎng)量的訂正提供了一種新的可選的方法。與傳統(tǒng)方法相結(jié)合，對(duì)不同樣芯靈活選取最合適的擬合方法將會(huì)得到更加準(zhǔn)確的氣候代用序列。該方法在改進(jìn)樹輪寬度生長(zhǎng)量訂正及標(biāo)準(zhǔn)化的應(yīng)用中有較大優(yōu)勢(shì)和潛力。然而，針對(duì)不同的樹種、不同的研究環(huán)境，對(duì)該方法的應(yīng)用仍需謹(jǐn)慎，需要更多的實(shí)驗(yàn)進(jìn)行驗(yàn)證和分析，尤其是對(duì)每一個(gè)樣芯的情況均需要仔細(xì)研究，綜合分析，才可能得到真正能反映氣候信息的可靠序列，得到對(duì)氣候變化的合理解釋和分析結(jié)果。

[1]ICSU.The Initial core projects[M].IGBP“Global Change”Report No.12,Berne,Switzerland.1990.330.

[2]張永香.巴丹吉林南緣近兩百年來年降水重建及初步分析[J].沙漠與綠洲氣象，2015，9（1）：12-16.

[3]尚華明，魏文壽，袁玉江，等.帕米爾東北部昆侖圓柏850年樹輪寬度年表的建立及其氣候意義[J].沙漠與綠洲氣象，2015，9（1）：6-11.

[4]Zhu h F,Zheng Y h,Shao X m,et al.Millennial

temperature reconstruction based on tree-ring widths of Qilian juniper from Wulan,Qinghai Province,China[J]. Chinese Science Bulletin，2008，53：3914-3920.

[5]袁玉江，喻樹龍，穆桂金，等.天山北坡瑪納斯河355a來年徑流量的重建與分析[J].冰川凍土，2005，27：411-417.

[6]Liang E,Liu X,Yuan Y,et al.The 1920s drought recorded by tree rings and historical documents in the semi-arid areas of Northern China[J].Climatic Change，2006，79：403-432.

[7]Yang B,Qin C,Wang J,Hem,et al.A 3，500-year treering record of annual precipitation on the northeastern Tibetan Plateau[J].PNAS，2014，111：2903-2908.

[8]Liu Y,Sun J,Song H,et al.Tree-ring hydrologic reconstructions for the Heihe River watershed,western China since AD 1430[J].Water Research，2010，44：2781-2792.

[9]Shao X,Huang L,Liu H,et al.Reconstruction of precipitation variation from tree rings in recent 1000 years in Delingha,Qinghai[J].Science in China Seriesd：Earth Sciences，2005，48：939-949.

[10]Zhang,Q B,Michael N E,Lyu L X.Moisture dipole over the Tibetan Plateau during the past five and a half centuries[J].Nature Communication.2015，6：8062，doi：10.1038/ncomms9062.

[11]Cook E R,Kairiukstis La.Methods of Dendrochronology[M].Netherlands:Kluwer Academic，1990

[12]Fritts h C.Tree Ringsand Climate[M].Academic Press, London，1976.

[13]吳祥定.樹木年輪與氣候變化[M].氣象出版社，1990.

[14]Phipps R L.Annual growth of suppressed chestnut oak and red maple,a basis for hydrologic inference[M],U.S.geological Survey Professional Paper 485-c，Washington D.C.1967.

[15]Cook E R,Briffa K R,Shiyatov S，et al.Tree-ring standardization and growth-trend estimation[M].In:Cook ER,Kairiukstis L A（eds.）.Methods of Dendrochronology.Kluwer Academic Publishers，1990，104-123.

[16]Cook E R.A time series analysis approach to tree-ring standardization[M].Disscrtation for the Doctoral Degree. Tucson:University of Arizona，1985，1-171.

[17]Mann M E,Bradley R S,Hughes M K.North hemisphere temperature during the last millennium:inferences, uncertainties and limitations[J].Geophysical Research Letters，1999，26：759-762.

[18]Briffa K R,Jones Pd,Schweingruber F H.Tree-ring reconstructions of summer temperature patterns across western North America since 1600[J].Journal of Climate，1992，5：735-754.

[19]Briffa K R,Jones Pd,Schweingruber F H，et al.Tree-ring variables as proxy-climate indicators:Problems with low frequency signals[J].Springer-Verlag，Berlin.、1996.

[20]Esper J,Cook E R,Schweingruber F H.Low-frequency signals in long tree-ring chronologies for reconstructing past temperature variability[J].Science，2002，295：2250-2253.

[21]Melvin，T M and Briffa,K R.A"Signal-Free"approach to Dendroclimatic Standardisation[J].Dendrochronologia 26，2008，71-86.doi:10.1016/j.dendro.2007.12.001.

[22]徐巖，邵雪梅.柴達(dá)木盆地東緣祁連圓柏輪寬序列標(biāo)準(zhǔn)化的方法研究[J].地理學(xué)報(bào)，2006，9：919-928.

[23]Huang N E,Shen Z,Long R S，et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary series analysis[J].Proc.R. Soc.Lond.A.1998，454：899-995.

[24]田沁花，周秀驥，勾曉華，等.祁連山中部近500年來降水重建序列分析[J].中國(guó)科學(xué)：地球科學(xué)，2012，42（4）：536-544.

Study of Standardization of the Qilian Juniper Ring-width Series Based on Empirical Mode Decomposition Method

ZHANG Yong
（Key Laboratory of Land Surface Pattern and Simulation,Institute of geographic Sciences and Natural Resources Research,Chinese academy of Sciences，Beijing 100101，China）

Standardization of tree-ring measurements is one of the most important procedures in dendroclimatology.In this paper,applying Empirical Mode decomposition（EMD）method to the Qilian Juniper（Sabina przewalskii kom）ring-width samples,the potential and limitations of the EMD method are investigated.Tree-ring series are firstly decomposed by EMD to produce a number of Intrinsic mode Functions（IMFs）in different frequency band with different physical meanings. The lowest-frequency signals mainly reflect the age-dependent growth trend,which are employed to correct the raw ring-width series to get the standardization chronology.Both the chronologies from the traditional and EMD method preserve the same high-frequency,but the different low-frequency variations in the 50-100yrs near to the pith.In this period,the variations of tree-ring index based on traditional method are large in general,while those based on EMD method are more stable.EMD method is a suitable tool for preserving more reliable climate information than traditional method, especially for the pith periods and young periods of tree-ring series,which can be applied to fit the growth trend of Qilian Juniper.

tree ring index;standardized method;Empirical Mode Decomposition（EMD）

P468

A

1002-0799（2015）05-0009-07

張永.經(jīng)驗(yàn)?zāi)B(tài)分解方法對(duì)祁連圓柏年輪寬度序列標(biāo)準(zhǔn)化初探[J].沙漠與綠洲氣象，2015，9（5）：9-15.

10.3969/j.issn.1002-0799.2015.05.002

2015-04-02；

2015-06-30

國(guó)家自然科學(xué)面上基金（41471087和41175066）資助。

張永（1979-），男，助理研究員，主要從事古氣候重建與氣候變化研究。E-mail：zhangyong@igsnrr.ac.cn