陳峰，袁玉江，張瑞波，喻樹(shù)龍

（1.中國(guó)氣象局烏魯木齊沙漠氣象研究所，新疆樹(shù)木年輪生態(tài)實(shí)驗(yàn)室，中國(guó)氣象局樹(shù)輪年輪理化研究重點(diǎn)開(kāi)放實(shí)驗(yàn)室，新疆烏魯木齊， 830002；2.蘭州大學(xué)西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，西部環(huán)境與氣候變化研究院，甘肅蘭州 730000）

基于樹(shù)輪寬度的巴丹吉林沙漠南緣山地NDVI重建與分析

陳峰1，2，袁玉江1，張瑞波1，喻樹(shù)龍1

（1.中國(guó)氣象局烏魯木齊沙漠氣象研究所，新疆樹(shù)木年輪生態(tài)實(shí)驗(yàn)室，中國(guó)氣象局樹(shù)輪年輪理化研究重點(diǎn)開(kāi)放實(shí)驗(yàn)室，新疆烏魯木齊， 830002；2.蘭州大學(xué)西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，西部環(huán)境與氣候變化研究院，甘肅蘭州 730000）

利用巴丹吉林沙漠南緣山地3個(gè)采樣點(diǎn)的青海云杉樹(shù)輪樣本研制出區(qū)域樹(shù)輪寬度年表。相關(guān)分析發(fā)現(xiàn)區(qū)域樹(shù)輪寬度年表與5—7月NDVI變化有較好相關(guān)性，相關(guān)系數(shù)為0.670（P<0.001）。利用線性回歸模型重建了巴丹吉林沙漠南緣山地1765—2010年5—7月NDVI變化，方差解釋量為44.9%。NDVI重建序列揭示在1765—2010年期間巴丹吉林沙漠南緣山地有10個(gè)植被生長(zhǎng)良好時(shí)段和10個(gè)植被生長(zhǎng)較差時(shí)段。同時(shí)，NDVI重建序列的低值與河西走廊極端干旱歷史事件有著良好的一致性。多窗譜周期分析發(fā)現(xiàn)，巴丹吉林沙漠南緣山地5—7月NDVI重建序列具有11.4，8.1，4.9，3.4，2.4 a的準(zhǔn)周期變化。交叉小波分析發(fā)現(xiàn)太陽(yáng)黑子活動(dòng)是區(qū)域NDVI變化的重要驅(qū)動(dòng)力之一。基于NDVI與區(qū)域干濕變化良好相關(guān)，分析了極端低值年份和極端高值年份的矢量風(fēng)場(chǎng)距平變化發(fā)現(xiàn)區(qū)域NDVI變化與大范圍氣候場(chǎng)變化有顯著關(guān)聯(lián)，發(fā)現(xiàn)當(dāng)西風(fēng)增強(qiáng)時(shí)，研究區(qū)氣候偏濕，有利于樹(shù)木生長(zhǎng)，NDVI值偏高，形成較寬的樹(shù)輪。當(dāng)研究區(qū)被來(lái)自北面沙漠地區(qū)氣流控制時(shí)候，研究區(qū)氣流偏干，不利于樹(shù)木生長(zhǎng)。

巴丹吉林沙漠；樹(shù)木年輪；NDVI重建；西風(fēng)

河西走廊南山——祁連山分布了大片原始針葉林，因其氣候及冰川水資源的變化對(duì)于下游河西走廊綠洲農(nóng)業(yè)區(qū)的重要影響，一直是中國(guó)樹(shù)輪氣候水文重建研究最為關(guān)注的地區(qū)之一[1-4]。而在河西走廊北部的巴丹吉林沙漠和騰格里沙漠南緣山地，由于其森林稀疏，人煙稀少，年降水不足150 mm，難以形成穩(wěn)定的地表徑流，在以往的樹(shù)輪氣候研究中較少受到關(guān)注[5]。巴丹吉林沙漠和騰格里沙漠南緣山地雖然屬于溫帶荒漠干旱氣候，但其北坡地形相對(duì)平緩，在陰坡和迎風(fēng)坡分布成片天然針葉林，最近幾年成為區(qū)域樹(shù)輪氣候環(huán)境研究熱點(diǎn)地區(qū)之一，研究人員在這一地區(qū)已建立了多條降水和干旱指數(shù)樹(shù)輪序列[5-7]，但針對(duì)該地區(qū)植被生長(zhǎng)狀況的歷史變化研究還相對(duì)較少。由于靠近巴丹吉林沙漠和騰格里沙漠，該地區(qū)的針葉林生長(zhǎng)變化對(duì)于改善區(qū)域生態(tài)環(huán)境和遏制土地沙漠化有著極為重要的作用。因此，我們有必要對(duì)其植被生長(zhǎng)狀況歷史變化進(jìn)行深入研究。

由GIMMS（global inventory monitoring and modeling studies）提供的歸一化植被指數(shù)（normalized difference vegetation index，NDVI）是目前最常用的植被生長(zhǎng)監(jiān)測(cè)資料之一[8]。由于其消除了火山爆發(fā)、太陽(yáng)高度角、傳感器靈敏度隨時(shí)間變化等影響，能夠有效指示區(qū)域植被生長(zhǎng)變化[9]。但由于該NDVI數(shù)據(jù)觀測(cè)時(shí)間較短，無(wú)法滿足植被生長(zhǎng)長(zhǎng)期變化的研究需要，而許多研究已經(jīng)證實(shí)通過(guò)建立樹(shù)輪指數(shù)與NDVI的關(guān)系，能夠有效利用樹(shù)輪資料解決這一問(wèn)題。目前，樹(shù)輪指數(shù)和植被指數(shù)之間關(guān)系的研究多集中在北半球中高緯度地區(qū)。Malmstr?m等[10]在美國(guó)阿拉斯加的研究結(jié)果表明樹(shù)木年輪寬度與利用NDVI估算得到森林凈初級(jí)生產(chǎn)力之間有顯著的正相關(guān)；Berner等[11]發(fā)現(xiàn)俄羅斯和加拿大高緯度地區(qū)1982—2008年夏季森林NDVI值與樹(shù)木年輪寬度正相關(guān)；在北美和歐亞大陸中高緯地區(qū)針葉林帶、美國(guó)中部大平原等均發(fā)現(xiàn)樹(shù)木年輪寬度與植被NDVI有較強(qiáng)的相關(guān)性[12-15]。近年來(lái)，在中國(guó)北方和青藏高原相關(guān)地區(qū)利用樹(shù)輪資料和NDVI數(shù)據(jù)進(jìn)行區(qū)域植被生長(zhǎng)變化重建研究取得了許多重要的進(jìn)展。例如He和Shao（2006）利用德令哈地區(qū)樹(shù)輪寬度指數(shù)重建了該地區(qū)過(guò)去1 000 a來(lái)草地8月的NDVI變化曲線[16]；Liang等（2009）的研究揭示了中國(guó)內(nèi)蒙古東部半干旱草原上樹(shù)輪寬度與草原生長(zhǎng)季NDVI之間的顯著正相關(guān)關(guān)系，表明可以利用樹(shù)輪寬度指示過(guò)去草原地上生物量的變化[17]。本文基于三個(gè)采樣點(diǎn)青海云杉（Picea crassifolia）樹(shù)輪資料，利用線性回歸模型重建巴丹吉林沙漠南緣山地NDVI變化，并以此為基礎(chǔ)分析巴丹吉林沙漠南緣山地NDVI的歷史變化特征。巴丹吉林沙漠南緣山地因其特殊的地理位置，是河西走廊對(duì)氣候變化的最為敏感區(qū)和脆弱區(qū)之一。該研究工作對(duì)于進(jìn)一步了解中國(guó)北方沙漠邊緣區(qū)植被變化提供數(shù)據(jù)支持，對(duì)認(rèn)識(shí)氣候變化背景下中國(guó)北方沙漠邊緣區(qū)的生態(tài)環(huán)境演變具有重要意義。

1 資料和方法

1.1 研究區(qū)概況及樹(shù)輪資料獲取

受溫帶大陸性氣候和西風(fēng)的綜合影響，巴丹吉林沙漠南緣山地以溫帶荒漠干旱氣候?yàn)橹?。研究區(qū)的主要?dú)夂蛱卣鞅憩F(xiàn)為降水稀少，溫度變化幅度大，雨熱同期。根據(jù)張掖氣象站1951—2010年的氣候資料，該地區(qū)年平均降水為130.6 mm，年平均溫度為7.3℃，7月和8月是溫度和降水最高的兩個(gè)月。隨著海拔高度變化，研究區(qū)降水和植被覆蓋有著明顯的變化，由下至上依次為荒漠草原帶，針葉林帶以及灌叢帶。林下土壤為灰褐森林土。山地針葉林的建群種為青海云杉（Picea crassifolia），青海云杉已經(jīng)被證明適于利用樹(shù)輪寬度進(jìn)行氣候變化研究[5]。

樹(shù)輪采樣工作分別于2009年7月和2011年6月進(jìn)行。在祁連山國(guó)家級(jí)自然保護(hù)區(qū)的東大山和龍首山兩個(gè)保護(hù)站的管轄區(qū)內(nèi)布設(shè)了3個(gè)采樣點(diǎn)（LLS，DDS1和DDS2），采樣樹(shù)種為青海云杉（圖1）。采樣點(diǎn)分布于海拔2 800～3 185 m之間，坡向?yàn)楸逼潞蜄|北坡，坡度為10～30°，森林郁閉度為0.1～0.3（表1）。在上述3個(gè)采樣點(diǎn)共計(jì)采集了100棵樹(shù)的193個(gè)樣芯。依照樹(shù)輪樣本處理的基本程序，首先對(duì)所采樣芯進(jìn)行晾干、粘貼、磨平。用砂紙將樣芯打磨至年輪輪廓清楚，細(xì)胞清晰可見(jiàn)之后，在顯微鏡下在樹(shù)芯表面進(jìn)行記號(hào)標(biāo)識(shí)。然后用精度為0.001 mm的Velmex輪寬測(cè)量?jī)x進(jìn)行輪寬測(cè)量。利用COFECHA程序?qū)徊娑杲Y(jié)果進(jìn)行檢驗(yàn)[18]。樹(shù)輪寬度年表的建立是利用ARSTAN程序完成的[19]。我們將3個(gè)采樣點(diǎn)的所有樹(shù)輪寬度序列利用hugershoff生長(zhǎng)函數(shù)進(jìn)行去除生長(zhǎng)趨勢(shì)處理[19]。由于樹(shù)輪序列之間存在較高相關(guān)系數(shù)（r=0.58），以雙權(quán)重平均法將上述去除生長(zhǎng)趨勢(shì)后的序列合并成均值為1、最小值為0的無(wú)量綱的區(qū)域樹(shù)輪寬度標(biāo)準(zhǔn)化年表。最后再使用自回歸模型消除森林內(nèi)部樹(shù)木之間競(jìng)爭(zhēng)所導(dǎo)致的低頻變化,得到區(qū)域樹(shù)輪寬度差值年表（圖2）。因?yàn)闃?shù)輪差值年表在重建巴丹吉林沙漠南緣山地NDVI效果更為理想，故在以下的研究中我們只展示樹(shù)輪寬度差值年表重建結(jié)果。

1.2 氣候和NDVI資料的選取

離研究區(qū)最近的，且觀測(cè)時(shí)間最長(zhǎng)的氣象站為張掖氣象站（38°56′N，100°26′E，海拔1 483.7 m，記錄時(shí)段1951—2008年）。經(jīng)初步檢驗(yàn)確認(rèn)，張掖氣象站的記錄可靠，無(wú)明顯的系統(tǒng)誤差?？紤]到前一年的氣候條件可能對(duì)當(dāng)年生長(zhǎng)季的樹(shù)木生長(zhǎng)產(chǎn)生一定影響，結(jié)合青海云杉的主要生長(zhǎng)季節(jié)（5—8月）[20]，本研究所采用的氣候要素為上一年7月至當(dāng)年9月的月平均溫度和月降水量。所使用的NDVI數(shù)據(jù)為GIMMS（global inventory monitoring and modeling studies）工作組提供的逐月平均NDVI數(shù)據(jù)集，空間分辨率8 km×8 km[8]。一般認(rèn)為NDVI值超過(guò)0.1表示有植被覆蓋，NDVI值越大則植被覆蓋率越高[21]。本研究從該數(shù)據(jù)集中提取出覆蓋研究區(qū)（38.90°～ 39.10°N，100.70°～101.00°E）的逐月平均NDVI值，選取的研究時(shí)段為1982年1月至2006年12月。在這一過(guò)程中，我們剔除了非植被地物對(duì)NDVI的干擾。從NDVI值的月變化來(lái)看，從5月到7月呈明顯上升趨勢(shì)，并7月達(dá)到峰值，之后逐漸下降。

利用相關(guān)函數(shù)（Pearson）分析樹(shù)輪寬度年表與氣候和NDVI的關(guān)系，并建立樹(shù)輪寬度與NDVI的回歸方程。為評(píng)價(jià)氣候重建效果，采用逐一剔除檢驗(yàn)方法進(jìn)行驗(yàn)證，計(jì)算的統(tǒng)計(jì)量有方差解釋量，符號(hào)檢驗(yàn)，誤差縮減值，效率系數(shù)等[22]。利用多窗譜分析對(duì)重建序列進(jìn)行周期分析[23]。

2 結(jié)果

2.1 樹(shù)輪寬度指數(shù)與氣候、NDVI的關(guān)系

對(duì)樹(shù)輪寬度差值年表與張掖氣象站的月平均溫度、月降水量和NDVI進(jìn)行相關(guān)分析的結(jié)果表明,差值年表與上年7—9月降水量和當(dāng)年5—7月降水量的相關(guān)關(guān)系十分顯著，超過(guò)了95%置信水平;而差值年表與月平均溫度的相關(guān)系數(shù)較低，沒(méi)有通過(guò)95%置信水平的月份（圖3）。從氣候相關(guān)分析結(jié)果看，巴丹吉林沙漠南緣山地降水對(duì)于青海云杉生長(zhǎng)的影響大于對(duì)溫度的影響，這符合研究區(qū)處于干旱區(qū)的實(shí)際情況。從對(duì)降水相關(guān)看,研究區(qū)青海云杉生長(zhǎng)與上年7月至當(dāng)年6月的降水量最為相關(guān)（r= 0.728），這種相關(guān)關(guān)系是由樹(shù)木生理過(guò)程所決定的。上年7月到當(dāng)年6月包括了上年生長(zhǎng)季末期、冬季樹(shù)木休眠期和當(dāng)年生長(zhǎng)季，該時(shí)段降水偏多不僅有利于當(dāng)年生長(zhǎng)季中青海云杉樹(shù)輪的生長(zhǎng)，還將對(duì)次年樹(shù)輪生長(zhǎng)產(chǎn)生影響[24]。以上分析揭示降水是研究區(qū)青海云杉生長(zhǎng)最為主要的氣候限制因子。從圖3中還可以看出，巴丹吉林沙漠南緣山地青海云杉當(dāng)年生長(zhǎng)季的NDVI與樹(shù)輪差值年表相關(guān)性較好。從5月開(kāi)始，青海云杉生長(zhǎng)逐漸旺盛，占樹(shù)輪寬度最主要部分--早材快速生長(zhǎng)，NDVI與青海云杉樹(shù)輪寬度相關(guān)系數(shù)也開(kāi)始上升，并于6—7月達(dá)到最高，這說(shuō)明樹(shù)輪寬度變化與生長(zhǎng)季內(nèi)的NDVI變化存在顯著關(guān)聯(lián)。其中，樹(shù)輪寬度差值年表與5—7月NDVI變化的聯(lián)系最為緊密，相關(guān)系數(shù)為0.670，通過(guò)了99%顯著性檢驗(yàn)。同時(shí)，我們還發(fā)現(xiàn)5—7月NDVI變化與上年7月至當(dāng)年6月的降水量也存在顯著相關(guān)（r=0.457，p<0.05）。在降水量相對(duì)充足，土壤水分的脅迫作用減弱的情況下，良好的水分條件會(huì)促進(jìn)青海云杉葉片的光合作用,導(dǎo)致NDVI值上升，增加青海云杉光合作用所產(chǎn)生的積累量，從而使得青海云杉徑向生長(zhǎng)增加，形成較寬的樹(shù)輪。反之，則會(huì)有較低的NDVI值和較窄樹(shù)輪。因此，研究區(qū)青海云杉樹(shù)輪寬度變化在可以反映巴丹吉林沙漠南緣山地降水變化的同時(shí)，也能夠揭示逐年生長(zhǎng)季內(nèi)NDVI的變化。

2.2 NDVI重建

基于巴丹吉林沙漠南緣山地區(qū)域樹(shù)輪寬度差值年表與生長(zhǎng)季的NDVI具有很好的相關(guān)關(guān)系，我們利用巴丹吉林沙漠南緣山地區(qū)域樹(shù)輪寬度差值年表重建生長(zhǎng)季的5—7月NDVI值，得到重建方程：

NDVI=0.12×RC+0.02，（1）

其中，NDVI為生長(zhǎng)季5—7月的平均值，RC為巴丹吉林沙漠南緣山地區(qū)域樹(shù)輪寬度差值年表。方程方差解釋量達(dá)到44.9%，調(diào)整后的解釋方差為42.5%。根據(jù)重建方程（1），并基于子樣本信號(hào)強(qiáng)度（SSS）>0.85的閥值[25]，重建了巴丹吉林沙漠南緣山地1765—2010年間的5—7月NDVI序列變化。圖4為校準(zhǔn)期1982—2006年間由區(qū)域樹(shù)輪差值年表重建的5—7月NDVI與實(shí)測(cè)值的對(duì)比，二者表現(xiàn)出很好的一致性。利用逐一剔除法對(duì)重建方程進(jìn)行穩(wěn)定性檢驗(yàn)，結(jié)果顯示其誤差縮減值和效率系數(shù)分別為0.358和0.333，一階差相關(guān)系數(shù)為0.48，達(dá)到了99%置信水平（表2），符號(hào)檢驗(yàn)也達(dá)到95%置信水平。以上檢驗(yàn)結(jié)果說(shuō)明重建方程穩(wěn)定可靠。3.3NDVI重建序列特征

對(duì)巴丹吉林沙漠南緣山地NDVI序列進(jìn)行15年低通濾波處理，顯示在1765—2010年間，NDVI經(jīng)歷了10個(gè)高值期和10個(gè)低值期，植被生長(zhǎng)較好的時(shí)期有：1772—1787年，1798—1804年，1823—1829年，1839—1848年，1866—1875年，1885—1911年，1935—1946年，1951—1958年，1969— 1984年和2005—2010年。植被生長(zhǎng)較差的時(shí)段有：1765—1771年，1788—1797年，1805—1822年，1830—1838年，1849—1865年，1876—1884年，1912—1934年，1947—1950年，1959—1968年和1985—2004年,其中1912—1934年這23 a植被生長(zhǎng)最差（圖5）。20世紀(jì)初以來(lái)的NDVI年際波動(dòng)幅度和時(shí)間間隔明顯高于之前的時(shí)段。從20世紀(jì)80年代末期以來(lái)，區(qū)域NDVI有顯著上升趨勢(shì)。在這246 a間植被生長(zhǎng)最好的10 a為:1786年，1790年，1855年，1888年，1889年，1899年，1955年，1980年，1982年，2010年。最差的10年為:1789年，1895年，1925年，1926年，1928年，1950年，1957年，1962年，1981年，1986年。另外,NDVI重建序列存在著顯著的周期性變化，主要的振蕩準(zhǔn)周期有11.4年，8.1年，4.9年，3.4年和2.4年。上述準(zhǔn)周期變化與騰格里沙漠東緣賀蘭山北部的3—7月帕爾默干旱指數(shù)重建序列周期變化具有良好的一致性[26]，這說(shuō)明同處亞洲戈壁沙漠南緣的兩者受類似的氣候波動(dòng)影響。

4 討論

將NDVI重建序列介于Mean（平均值）±σ（標(biāo)準(zhǔn)差）值之間的范圍視為巴丹吉林沙漠南緣山地植被正常生長(zhǎng)變化，把NDVI

NDVI重建序列2.4～4.9 a的變化周期可能與ENSO有關(guān)[30]。在本NDVI重建序列中發(fā)現(xiàn)的11.4年變化周期說(shuō)明了區(qū)域植被生長(zhǎng)變化不僅受到干旱脅迫的影響，同時(shí)還受到太陽(yáng)活動(dòng)影響[27]。交叉小波分析也進(jìn)一步證實(shí)區(qū)域NDVI變化與太陽(yáng)黑子活動(dòng)的11年周期波動(dòng)有著顯著關(guān)聯(lián)（圖6）。

現(xiàn)代氣候觀測(cè)資料證實(shí)研究區(qū)氣候主要受到西風(fēng)活動(dòng)控制。研究區(qū)降水重建結(jié)果也發(fā)現(xiàn)研究區(qū)降水與太平洋海溫、大西洋海溫有著緊密關(guān)聯(lián)，并深受西風(fēng)和東亞夏季風(fēng)交互作用影響[5，31]。我們從NDVI重建序列中挑選有氣候資料記錄（1948—2010年）的5個(gè)最高值和6個(gè)最低值年份，繪制出5—7月500 hPa矢量風(fēng)場(chǎng)距平合成圖（圖7）。在NDVI值偏高年份，中緯度西風(fēng)輸送增強(qiáng)，亞洲戈壁沙漠地區(qū)受到西風(fēng)氣流控制，同時(shí)來(lái)自低緯地區(qū)的濕潤(rùn)氣流能夠更加深入內(nèi)陸，導(dǎo)致巴丹吉林沙漠南緣山地降水偏多，從而緩解干旱壓力，促進(jìn)研究區(qū)青海云杉生長(zhǎng)。相反，在NDVI值偏低年份，中緯度西風(fēng)氣流偏南，亞洲夏季風(fēng)偏弱，研究區(qū)主要受到亞洲戈壁沙漠腹地干燥氣流控制，干旱氣候加重了干旱壓力，不利于青海云杉生長(zhǎng)。以上分析也佐證了研究區(qū)氣候和植被生長(zhǎng)變化受到了西風(fēng)和亞洲夏季風(fēng)共同作用，其中西風(fēng)驅(qū)動(dòng)的氣候變化在研究區(qū)氣候環(huán)境變化中起到了主導(dǎo)作用[32-36]。最近30年來(lái)NDVI上升趨勢(shì)的發(fā)生與中國(guó)西北西部暖濕化相同步[32-35]，這說(shuō)明兩個(gè)地區(qū)很可能具有類似的水汽來(lái)源，并與西風(fēng)活動(dòng)有著緊密聯(lián)系[5，26]。這一暖濕化趨勢(shì)為研究區(qū)脆弱的生態(tài)環(huán)境改善提供良好的氣候條件，這不僅直接反映在降水和濕度得到增加方面，同時(shí)在促進(jìn)區(qū)域植被生長(zhǎng)，改善生態(tài)效應(yīng)方面得到了體現(xiàn)。

5 結(jié)論

（1）巴丹吉林沙漠南緣山地近246 a的5—7月的NDVI經(jīng)歷了10個(gè)植被生長(zhǎng)良好階段和10個(gè)植被生長(zhǎng)較差階段，其中NDVI重建序列中的低值與河西走廊地區(qū)極端干旱事件有對(duì)應(yīng)關(guān)系。

（2）巴丹吉林沙漠南緣山地NDVI重建序列存在11.4，8.1，4.9，3.4和2.4 a的變化準(zhǔn)周期，與同處亞洲戈壁沙漠南緣的賀蘭山北部地區(qū)樹(shù)輪氣候序列的周期具有良好的一致性。同時(shí)，交叉小波分析證實(shí)太陽(yáng)黑子活動(dòng)對(duì)于區(qū)域NDVI值變化有著重要影響大。

（3）5—7月500 hPa矢量風(fēng)場(chǎng)距平表明巴丹吉林沙漠南緣山地NDVI與區(qū)域氣候變化有著顯著作用。通過(guò)矢量風(fēng)場(chǎng)距平分析發(fā)現(xiàn)，巴丹吉林沙漠南緣山地NDVI主要受到西風(fēng)變化的影響，當(dāng)西風(fēng)增強(qiáng)時(shí)，降水有可能增加，有利于青海云杉生長(zhǎng)，同時(shí)也與季風(fēng)氣流活動(dòng)有著顯著相關(guān)。

[1]陳峰，魏文壽，袁玉江，等.基于多點(diǎn)樹(shù)輪序列的1768—2006年甘肅降水量變化[J].中國(guó)沙漠，2013，33（5）：1520-1526.

[2]Yang B，Qin C，Shi F，et al.Tree ring-based annual streamflow reconstruction for the heihe river in arid northwestern china from AD 575 and its implications for water resource management[J].The Holocene，2012，22：773-784.

[3]ChenF,YuanYJ,WeiWS,etal.Temperature reconstruction from tree-ring maximum latewood density of QinghaispruceinmiddleHexiCorridor,China[J]. Theoretical and Applied Climatology，2011，107：633-643.

[4]李穎俊，勾曉華，方克艷，等.祁連山東部188 a上年8月至當(dāng)年6月降水量的樹(shù)輪重建[J].中國(guó)沙漠，2012，32（5）：1393-1401.

[5]Chen F，Yuan YJ，Wei WS，et al.Tree-ring-based annual precipitation reconstruction for the Hexi Corridor,NW China:consequences for climate history on and beyond the mid-latitude Asian continent[J].Boreas，2013，42：1008-1021.

[6]Kang SY，Yang B，Qin C.Recent tree-growth reduction in north central China as a combined result of a weakened monsoon and atmospheric oscillations[J].Climatic Change，2012，115（3-4）：519-536.

[7]Liu Y，Lei Y，Sun B，et al.Annual precipitation variability inferredfromtree-ringwidthchronologiesinthe Changling-Shoulu region,China,during AD 1853-2007[J]. Dendrochronologia，2013，31（4）：290-296.

[8]Tucker CJ，Slayback DA，Pinzon JE，et al.Higher northern latitude normalized difference vegetation index andgrowingseasontrendsfrom1982to1999. International Journal of Biometeorology，2001，45：184-190.

[9]Holben BN.Characteristics of maximum-value composite images for temporal AVHRR data.International Journal of Remote Sensing，1986，7：1435-1445.

[10]Malmstr?mCM，ThompsonMV，JudayGP，etal. Interannualvariationinglobalscalenetprimary production：testingmodelestimates.Global Biogeochemical Cycles，1997，11（3）：367-392.

[11]Berner LT，Beck PSA，Bunn AG，et al.High-latitude tree growth and satellite vegetation indices:correlations and trends in Russia and Canada（1982-2008）.Journal of Geophysical Research，2011，116：G01015.

[12]Kaufmann RK，D’Arrigo RD，Laskowski C，et al.The effect of growing season andsummer greenness on northern forests.Geophysical Research Letters，2004，31：L09205.

[13]D’Arrigo RD，Malmstr?m CM，Jacoby GC，et al. Correlation between maximum latewood density of annual treeringsandNDVIbasedestimatesofforest productivity.Journal of Remote Sensing，2000，21（11）：2329-2336.

[14]Wang J，Rich PM，Price KP，et al.Relations between NDVI and tree productivity in the central Great Plains. International Journal of Remote Sensing，2004，25（16）：3127-3138.

[15]Lopatin E，Kolstr?m T，Spiecker H.Determination of forest growth trends in Komi Republic（northwestern Russia）：combination of tree-ring analysis and remote sensingdata.BorealEnvironmentResearch，2006，11（5）：341-353.

[16]He JC，Shao XM.Relationships between tree-ring width index and NDVI of grassland in Delingha.Chinese Science Bulletin，2006，51：1106-1114.

[17]Liang EY，Eckstein D，Liu HY.Assessing the recent grassland greening trend in a long-term context based on tree-ringanalysis:acasestudyinNorthChina. Ecological Indicators，2009，9：1280-1283.

[18]Holmes RL.Computer-assisted quality control in treering dating and measurement[J].Tree-ring Bulletin, 1983，43：69-75.

[19]Cook ER,Kairiukstis LA.Methods of Dendrochronology [M].Dordrecht：Kluwer Academic Press，1990.

[20]李雁，梁爾源，邵雪梅.柴達(dá)木盆地東緣青海云杉樹(shù)輪細(xì)胞結(jié)構(gòu)變化特征及其對(duì)氣候的指示[J].應(yīng)用生態(tài)學(xué)報(bào)，19（3）：524-532.

[21]王文志，劉曉宏，陳拓，等.基于祁連山樹(shù)輪寬度指數(shù)的區(qū)域NDVI重建.植物生態(tài)學(xué)報(bào)，2010，34（9）：1033-1044.

[22]李江風(fēng)，袁玉江，由希堯.樹(shù)木年輪水文學(xué)研究與應(yīng)用[M].北京：科學(xué)出版社，2000.

[23]Mann ME，LEEJM.Robustestimateestimationof background noise and signal detection in climatic time series[J].Climatic Change，1996，33：409-445.

[24]陳峰，袁玉江，魏文壽，等.樹(shù)輪記錄的甘肅山丹近224年降水變化[J].地理與地理信息科學(xué)，2010，26（5）：82-86.

[25]Wigley TML，Briffa KR，Jones PD.On the average value ofcorrelatedtimeseries,withapplicationsin dendroclimatology and hydrometeorology[J].Journal of Applied Meteorology and Climatology，1984，23：201-213.

[26]Li J,Chen F,Cook E,et al.Drought reconstruction for north central China from tree rings：The value of the Palmer drought severity index[J].International Journal of Climatology，2007，27：903-909.

[27]袁林.西北災(zāi)荒史[M].蘭州：甘肅人民出版社，1994：479-618.

[28]成愛(ài)芳，趙景波，馮起.甘肅隴中地區(qū)清代干旱災(zāi)害與氣候變化研究[J].中國(guó)沙漠，2011，31（6）：1611-1616.

[29]Wang HJ.The Weakening oftheAsianMonsoon Circulation after the End of 1970’s[J].Advances in Atmospheric Sciences，2001，18：376-386.

[30]朱炳瑗，李棟梁.ENSO現(xiàn)象與甘肅省夏季降水[J].高原氣象，1989，8（1）：64-69.

[31]Yang B，Qin C,Br?uning A，et al.Rainfall history for the Hexi Corridor in the arid northwest China during the past 620 years derived from tree rings[J].International Journal of Climatology，2011，31（8）：1166-1176.

[32]陳峰，袁玉江，魏文壽，等.呼圖壁河流域樹(shù)輪年表的建立及其特征分析[J].沙漠與綠洲氣象，2008，2（2）：7-10.

[33]許崇海，徐影，羅勇.新疆地區(qū)21世紀(jì)氣候變化分析[J].沙漠與綠洲氣象，2008，2（3）：1-7.

[34]陳峰，王慧琴，袁玉江，等.樹(shù)輪最大密度記錄的吉爾吉斯斯坦天山山區(qū)公元1650年以來(lái)的7—8月溫度變化[J].沙漠與綠洲氣象，2014，8（4）：1-7.

[35]施雅風(fēng)，沈永平，李棟梁，等.中國(guó)西北氣候由暖干向暖濕轉(zhuǎn)型的特征和趨勢(shì)探討[J].第四紀(jì)研究，2003，23（2）：152-164.

[36]黃偉，陳建徽，張肖劍，等.現(xiàn)代氣候條件下降水變化的“西風(fēng)模態(tài)”空間范圍及其影響因子初探.中國(guó)科學(xué)：地球科學(xué)，2015，45：379-388.

Reconstruction of Regional NDVI Using Tree-ring Width Chronology in the Southern Edge of the Badain Jaran Desert,Northwestern China

CHEN Feng1，2，YUAN Yujiang1，ZHANG Ruibo1，YU Shulong1
（1.Xinjiang Laboratory of Tree-ring Ecology/Key Laboratory of Tree-ring Physical and Chemical Research of China Meteorological Administration,Institute of Desert Meteorology，China Meteorological Administration，Urumqi 830002，China；2.Key Laboratory of Western China’s Environmental Systems（MOE），Lanzhou University，Lanzhou 730000，China）

A regional tree-ring chronology of Picea crassifolia spanning 1765-2010 was developed from the southern edge of the Badain Jaran Desert,northwestern China.The climate/NDVI response analysis shows the regional tree-ring chronology have a good correlation（r=0.670）with May-July NDVI.Based on the regional tree-ring chronology,we developed a May-July NDVI reconstruction of the southern edge of the Badain Jaran Desert for the period AD 1765-2010.The NDVI/treegrowth model accounts for 44.9%of the observed NDVI variance during the period 1982-2006.Ten high NDVI periods（NDVI values higher than the average of reconstruction）and ten low NDVI periods（NDVI values lower than the average）was found during 1765-2010 respectively,and reflected the strong-weak changes of the Westerlies.There is a reasonable agreement with climatic extremes previously estimated from documentary records in northwest China.The results reveal common climatic extremes over much of Gansu.Many of these events have had profound impacts on the peoples of the Hexi Corridor over the past several centuries.The multitaper method（MTM）spectral analysis indicates the existence of 11.4 years，8.1 years，4.9 years，3.4 years and 2.4 years cycles,which may potentially be the fingerprints of some proposed climate change forcings. Preliminary analysis of links between large-scale climatic variation and the NDVI reconstructions shows that there is a relationship between extremes in NDVI reconstruction and anomalous atmospheric circulation in the region，and revealed that the westerlies have the great effects on the NDVI changes of the study area.

Badain Jaran Desert；tree-ring；NDVI reconstruction；westerlies

P467

B

1002-0799（2015）04-0001-07

陳峰，袁玉江，張瑞波，等.基于樹(shù)輪寬度的巴丹吉林沙漠南緣山地NDVI重建與分析[J].沙漠與綠洲氣象，2015，9（4）：1-7.

10.3969/j.issn.1002-0799.2015.04.001

2014-11-07；

2015-05-30

科技部公益性行業(yè)（氣象)）科研專項(xiàng)“沙漠地區(qū)歷史氣候變化序列重建研究”（GYHY201206014）；中國(guó)沙漠氣象科學(xué)研究基金“樹(shù)輪記錄的新疆東部近300年干濕變化及其與周邊地區(qū)之關(guān)聯(lián)”（Sqj2013015）；中央級(jí)基本科研業(yè)務(wù)費(fèi)項(xiàng)目“基于樹(shù)輪密度的秦嶺中西部溫度的重建和分析”（IDM201105）共同資助。

陳峰（1982-），男，副研究員，主要從事樹(shù)木年輪與氣候變化研究。E-mail：feng653@163.com