張 啟,方歐婭

1 山西師范大學(xué)生命科學(xué)學(xué)院,臨汾 041004 2 中國科學(xué)院植物研究所植被與環(huán)境變化國家重點(diǎn)實驗室,北京 100093

森林健康是森林生態(tài)系統(tǒng)健康與恢復(fù)的簡稱,是指森林在發(fā)揮必要的生態(tài)服務(wù)功能的同時,保持其自身良性存在和更新的狀態(tài)[1]。近年來由于全球氣候變化所帶來的高溫和干旱等極端氣候事件的頻發(fā),使得世界上大多數(shù)國家的森林生態(tài)系統(tǒng)都存在不同程度的退化[2-4],主要表現(xiàn)為樹木在生長發(fā)育過程中的生理機(jī)能下降、生產(chǎn)力降低、生長發(fā)育滯緩或死亡等現(xiàn)象。在我國,森林覆蓋率已由建國初期的12.5%上升至2013年的21.63%[5]。然而與快速恢復(fù)性增長的森林覆蓋率呈鮮明反差的是,許多現(xiàn)存的森林處于破碎化或次生演替狀態(tài)。中國森林退化問題突出,森林健康狀況堪憂,影響了森林發(fā)揮必要的生態(tài)服務(wù)功能[1],使得森林健康的研究十分迫切。

目前國內(nèi)外學(xué)者對森林健康的研究提出了許多評價方法[6-8],但是由于對森林群落缺乏長期觀測,人們難以得知其健康歷史。樹木年輪記錄了樹木個體的生長歷史,具有精確定年、分辨率高、連續(xù)性強(qiáng)和復(fù)本量好等特點(diǎn)[9],在反映過去森林衰退方面有著突出的貢獻(xiàn)[10-12]。研究表明,樹木在生長過程中會受到不同范圍生物和非生物多元化組合的影響,可能出現(xiàn)樹冠枯死、低生長或死亡等不健康生長的現(xiàn)象[4,13-14]。由于限制因子發(fā)生的順序、時間和強(qiáng)度的不同,其影響可能是持續(xù)或者短暫的,森林健康在空間上的表現(xiàn)較為復(fù)雜[2,15-17]。因此,過去的森林不健康事件在發(fā)生程度和持續(xù)時間上有何特征是一個亟待解決的生態(tài)學(xué)問題。對森林健康歷史時空特征以及其驅(qū)動因子的深入研究很大程度上豐富了歷史時期森林生態(tài)學(xué)特征的研究。

忻州地處山西高原中北部,位于呂梁山生態(tài)脆弱區(qū)。針對上述問題,我們該選擇地區(qū)保存較好的4個森林為研究對象,試圖利用樹木年輪學(xué)評價該地區(qū)森林健康狀況,旨在研究以下兩個具體問題：1)忻州地區(qū)森林健康狀況在過去100年中發(fā)生了怎樣的變化？2)研究時間段內(nèi)多次森林不健康事件在空間上有無同步性？研究結(jié)果所得森林健康時空信息也可為林業(yè)部門制定合理的經(jīng)營管理方案提供數(shù)據(jù)依據(jù)。

1 研究材料與方法

1.1 研究區(qū)域概況

本研究區(qū)位于山西省忻州市地區(qū)(38°08′—39°40′N,110°53′—113°58′E),西臨黃河,東止太行山。該區(qū)地形西高東低,逐步傾斜,北、西、南三面環(huán)山,東部開闊平坦,為忻定盆地的主體部分。境內(nèi)流徑滹沱河、云中河及牧馬河等重要河流。由于受溫帶大陸性季風(fēng)氣候的影響,該地區(qū)夏季多東南風(fēng),冬季多西北風(fēng),春溫高于秋溫。年平均氣溫為8.5℃左右,年平均降水量為405 mm,降水集中于7—9月份[18]。

1.2 樹輪樣品采集與氣象資料

忻州地區(qū)自然植被資源豐富,全區(qū)森林面積4.07×105 hm2,森林覆蓋率為16.4%。本研究選取忻州地區(qū)植被覆蓋率較高的4個縣作為采樣點(diǎn)。其中寧武梅洞(38°71′N,111°96′E；海拔1772 m；青扦Piceawilsonii)和大石洞林場(38°92′N,112°08′E；海拔1752 m；青扦)屬于呂梁山脈,坡度5°—10°。繁峙沙河寬灘林場(39°06′N,113°47′E；海拔2184 m；臭冷衫Abiesnephrolepis)屬于五臺山區(qū)域,坡度為0°—5°。岢嵐縣中寨林場(38°58′N,111°43′E；海拔1745 m；油松Pinustabuliformis)屬于黑茶山區(qū)域,坡度為15°—20°(圖1)。4個采樣點(diǎn)區(qū)域天然次生林植被豐富,主要有白扦林、青扦林、華北落葉松林、白樺林、山楊林、遼東櫟林、油松林等。野外樣品采集于2012年10月,我們選擇區(qū)域內(nèi)年齡較老的樹木,使用生長錐在胸高處(1.3 m)采集樹芯樣本。取樣方向為沿山坡的等高線方向,每株立木采集一根樹芯樣本,共采集了86棵樹的樹芯樣本。

圖1 采樣點(diǎn)示意圖Fig.1 The location map of the sampling sites

本研究選用距離采樣點(diǎn)較近位置的原平氣象站(38°44′N,112°43′E；海拔828.2 m,數(shù)據(jù)自1954年起測定)和五寨氣象站(38°55′N,111°49′E；海拔1401 m,數(shù)據(jù)自1957年起測定)的氣候數(shù)據(jù)(圖1)。氣候因子包括月總降水量(Monthly total precipitation,Pm)以及月平均氣溫(Monthly mean temperature,Tm)。由于兩地同屬暖溫帶季風(fēng)氣候區(qū)域,氣象站數(shù)據(jù)差異性較小,因此將兩個氣象站的氣候數(shù)據(jù)平均值作為該區(qū)域的氣候數(shù)據(jù)[19],分別分析不同樣點(diǎn)與區(qū)域氣候的響應(yīng)關(guān)系。干旱指數(shù)Palmer drought severity index(PDSI)數(shù)據(jù)來自荷蘭皇家氣象研究所的數(shù)據(jù)共享網(wǎng)(http://climexp.knmi.nl/),區(qū)域范圍為38.59°—39.06°N,111.43°—113.47°E。

1.3 樹輪數(shù)據(jù)的分析方法

樣本在自然狀態(tài)下干燥后,按照基本程序進(jìn)行固定、磨光處理和交叉定年[20]。在測量精度為0.001 mm的LINTAB樹輪寬度測量儀(RINNTECH,Heidelberg,German)上測量每一年輪的寬度,并利用國際樹木年輪庫的軟件程序COFECHA[21]進(jìn)行定年質(zhì)量控制。更正因缺輪、偽輪和測量誤差等帶來的錯誤,最終使樹芯樣本上的每一個樹輪都對應(yīng)其實際生長年份。應(yīng)用ARSTAN[22]程序,用負(fù)指數(shù)函數(shù)或直線擬合每個測量序列的生長趨勢,用加權(quán)平均法將去趨勢序列合并成樹輪寬度指數(shù)序列。研究計算了寧武縣兩個采樣點(diǎn)年表間的滑動相關(guān)系數(shù)矩陣,將梅洞和大石洞的樣本合并。同時設(shè)定了分析健康歷史的起始年,采用子樣本信號強(qiáng)度SSS(Sub-sample Signal Strength)>85%的樣本量作為分析的起始點(diǎn),高于該樣本量的健康序列認(rèn)為是可靠的。

森林中樹木個體生長的健康狀況由于樹木自身或局部生長小環(huán)境的原因會產(chǎn)生差別,在不利環(huán)境條件下,一些樹木具有較強(qiáng)的恢復(fù)力而能夠及時響應(yīng)環(huán)境的變化,而一些樹木受到外界影響后不能及時恢復(fù)到健康的狀況,從而導(dǎo)致了樹木的健康與不健康的分異。不健康樹木較健康樹木年輪生長較慢,而且不健康樹木的低生長是一個持續(xù)的過程。結(jié)合樹木年輪數(shù)據(jù),定義樹木的不健康時段表現(xiàn)為樹輪指數(shù)連續(xù)5 a低于正常值?；诖硕x,從下面3個步驟挑選出不健康時段的樹木：1)挑出樹木至少持續(xù)5 a低生長片段(樹輪指數(shù)小于0.9)；2)挑出同時期樹輪指數(shù)平均值低于0.7,且期間指數(shù)最小值低于0.5的樹為非健康樹；3)統(tǒng)計各點(diǎn)每年中出現(xiàn)非健康樹的樣本量所占百分比。將挑出的非健康樹年份中的樹輪值平均作為非健康樹年表,同時將非健康樹輪值剔除平均剩余值作為健康樹年表。通過對比各采樣點(diǎn)結(jié)果,找出該區(qū)域時間和空間上非健康歷史的特征。

為了研究該地區(qū)森林生長在對當(dāng)?shù)貧夂蜃兓捻憫?yīng)特征,研究利用DendroClim2002程序[23]分別分析各點(diǎn)的樹輪寬度標(biāo)準(zhǔn)年表與逐月氣候因子的相關(guān)關(guān)系,逐月氣候因子包括前1年10月到當(dāng)年9月的氣候指標(biāo)(月總降水量和月平均氣溫、PDSI)。

2 研究結(jié)果

2.1 樹輪樣本及交叉定年結(jié)果

經(jīng)過對樹芯樣本測量和交叉定年后,得到了共86個年齡準(zhǔn)確可用于分析的樹輪數(shù)據(jù)。各年表的統(tǒng)計特征見表1,其中該地區(qū)年齡最大的樹為223 a,年齡最小的樹為25 a。各采樣點(diǎn)序列間平均相關(guān)系數(shù)都達(dá)到了0.5以上,表明該區(qū)域內(nèi)各點(diǎn)樹木個體間的輪寬變化均較為一致,對外界環(huán)境變化有著一致的響應(yīng)。一階自相關(guān)系數(shù)也都達(dá)到了0.6以上,表明上一年氣候狀況對當(dāng)年樹輪生長有著持續(xù)性影響。樣本總體代表性都高于0.85,各樣點(diǎn)含著大量的共同信號。

表1 各采樣點(diǎn)標(biāo)準(zhǔn)年表的統(tǒng)計特征Table 1 Statistic characteristics of standard tree-ring chronologies in the sampling sites

2.2 樹輪和氣候因子的響應(yīng)關(guān)系

各樣點(diǎn)標(biāo)準(zhǔn)年表與氣候響應(yīng)分析結(jié)果如圖2所示。從圖中可以看出,研究區(qū)域樹輪寬度指數(shù)與月平均氣溫和月總降水量的相關(guān)結(jié)果有較好的一致性。各樣點(diǎn)樹輪主要響應(yīng)于生長季中期的溫度與生長季末期的降水,其中,以與當(dāng)年5—6月月平均溫度的負(fù)相關(guān)關(guān)系最為顯著(r=-0.51—-0.26,P<0.05)；在樹木徑向生長與降水因子的關(guān)系中,主要與7月份的降水顯著的正相關(guān)(r=0.29—0.43,P<0.05)。從樹輪指數(shù)與PDSI相關(guān)的分布圖可以看出,4個采樣點(diǎn)樹輪與該區(qū)域PDSI呈顯著的正相關(guān)關(guān)系(r=0.22—0.55,P<0.05)。

圖2 樹輪標(biāo)準(zhǔn)年表和氣候因子的相關(guān)關(guān)系Fig.2 Correlation between standard tree-ring chronologies and climate factors左圖為各采樣點(diǎn)標(biāo)準(zhǔn)年表和PDSI相關(guān)關(guān)系,下圖為采樣點(diǎn)和溫度降水的相關(guān)關(guān)系; P為上一年月份,C為當(dāng)年月份,虛線和*表示該月份相關(guān)達(dá)到95%顯著水平

2.3 健康樹與非健康樹的徑向生長歷史

各采樣點(diǎn)健康和非健康樹輪指數(shù)序列如圖3所示。結(jié)果表明,岢嵐地區(qū)健康樹1900—1930年間,生長量呈現(xiàn)明顯的下降趨勢,此時非健康樹樹輪指數(shù)也較低,其中1910—1930年期間均值僅為0.43；在1960年之后健康樹和非健康樹也表現(xiàn)出低生長的趨勢,一直持續(xù)到2005年。繁峙地區(qū)非健康樹在1920—1940年間樹輪指數(shù)均值為0.55；健康樹和非健康樹在1990—2010年期間表現(xiàn)出明顯的持久下降趨勢。寧武地區(qū)非健康樹在1910—1940年期間,樹輪指數(shù)均值為0.51,且在1926年極低；在1970—1980年和2000年以后樹輪指數(shù)較低,均值分別為0.58和0.54存在穩(wěn)定下降趨勢。

圖3 各采樣點(diǎn)健康樹年表和非健康樹年表Fig.3 Ring-width index of healthy and unhealthy trees in the four sampling sites黑色和灰色折線分別為健康樹和非健康樹樹輪指數(shù)值及其對應(yīng)的11年滑動平均值,虛線為SSS>85%的起始年

研究統(tǒng)計了各采樣點(diǎn)非健康樹所占百分比情況(圖4),發(fā)現(xiàn)樣本量百分比峰值年較為一致。在過去1個世紀(jì)的公共時期內(nèi),岢嵐地區(qū)最為顯著的不健康事件發(fā)生在1913—1936年期間,持續(xù)了25年,通過定義判斷的非健康樹木量達(dá)總體的40%以上。而寧武地區(qū)的非健康樹樣本量從1915年就開始持續(xù)增多,1923年達(dá)到峰值,期間非健康樹所占比例達(dá)到30%以上。超過30%樣本量的時期持續(xù)了18 a。繁峙地區(qū)為1925—1942年,持續(xù)了18 a。70年代的非健康事件僅在岢嵐和寧武地區(qū)有所體現(xiàn),都是持續(xù)了5 a,寧武從1973年開始發(fā)生,岢嵐從1983年開始。20世紀(jì)末,該區(qū)域的不健康事件在3個研究區(qū)域都有體現(xiàn),都是從1999年開始。其中沙河地區(qū)持續(xù)時間最長為14 a,一直持續(xù)到采樣時間；岢嵐地區(qū)11 a；寧武地區(qū)持續(xù)時間較短,為6 a。

圖4 各采樣點(diǎn)非健康樹樣本所占百分比情況Fig.4 Percentage of unhealthy trees in the four sampling sites灰色陰影區(qū)為非健康樹樣本量百分比大于30%以上的時期,虛線為總樣本量

3 討論

3.1 樹木徑向生長與氣候的關(guān)系

忻州市地處中緯度地帶內(nèi)陸,受大陸性季風(fēng)氣候以及北部地區(qū)的風(fēng)蝕,使得該地區(qū)的年蒸發(fā)量遠(yuǎn)超過于降水量。在干旱及半干旱地區(qū),降水量是樹木生長的主要限制因子之一[24]。本文中樹輪寬度主要與生長季中期5—6月的月平均溫度以及7月總降水量呈顯著的相關(guān)關(guān)系。樹輪寬度與降水和溫度的這種相關(guān)關(guān)系與戴君虎等[25]在山西忻州五臺山地區(qū)、李穎俊等[26]在蘆芽山地區(qū)、蔡秋芳等[27]在呂梁山地區(qū)的研究結(jié)果相似。忻州地區(qū)5—6月份降水量不足,較高的溫度必然會加速樹木的蒸騰量,使樹木體內(nèi)的水分散失過快,同時也會加速土壤水分的蒸發(fā)量,形成水分脅迫限制了樹木的生理代謝活動,從而降低了樹木的生長量。7月份降水充足,同時較高的氣溫和地溫也使該地區(qū)凍土融化,有效的提供樹木徑向細(xì)胞分裂和伸長所需的必要水分,促進(jìn)了樹木的生長。樹輪與PDSI的相關(guān)系數(shù)結(jié)果要高于降水或溫度(圖4),也表明該地區(qū)的樹木主要受到水分作用的限制,更多地體現(xiàn)在土壤水分的有效供應(yīng)方面。Cai等[28]在山西太行山中部地區(qū)和Sun等[29]在汾河流域的研究也得到了類似的結(jié)果。

3.2 森林健康的時空特征

由于微環(huán)境、限制因子、樹木個體耐受力、群落競爭等狀況的影響,樹木生長個體間存在著差異[30]。本文利用生長指數(shù)的差異性挑選出特定時期健康樹和非健康樹,得到了該區(qū)域在時間和空間尺度上的森林健康信息,并且非健康樹所占百分比的統(tǒng)計進(jìn)一步補(bǔ)充了不健康事件發(fā)生的程度。各采樣點(diǎn)非健康樹樣本情況(圖4)顯示該區(qū)域出現(xiàn)高于30%閾值的不健康事件在空間上表現(xiàn)出明顯的同步性,主要集中在1910—1940、1970—1987和1990—2012年期間。各地區(qū)不同時期的不健康事件持續(xù)時間和強(qiáng)度不盡相同,如岢嵐地區(qū)在20世紀(jì)30年代不健康事件持續(xù)時間最長強(qiáng)度也最大。20世紀(jì)70年代時期,岢嵐和寧武非健康樹比例都較高,而該時期繁峙地區(qū)非健康僅占少數(shù)比例。20世紀(jì)末,繁峙地區(qū)不健康事件持續(xù)時間最長。梁爾源等[31]對中國中西部地區(qū)樹輪研究發(fā)現(xiàn),極端干旱事件導(dǎo)致20世紀(jì)20年代和20世紀(jì)30年代初期大范圍的樹木生長下降現(xiàn)象。Liu等[10]探究了亞洲中部半干旱地區(qū)森林衰退現(xiàn)象,也發(fā)現(xiàn)20世紀(jì)70年代末80年代初和21世紀(jì)初期樹輪表現(xiàn)出顯著的生長下降。我們認(rèn)為,大尺度的氣候事件導(dǎo)致森林不健康事件發(fā)生的空間格局具有很強(qiáng)的一致性。不同林分特征的森林對氣候變化的抗干擾能力和恢復(fù)能力有顯著的差異[32],并且不同樹種對氣候變化的敏感度不同[30],這可能是影響不健康事件發(fā)生強(qiáng)度和持續(xù)時間不一致的重要因素。

3.3 干旱事件對森林健康的指示

Liu等[10]認(rèn)為氣候變暖和干旱加速了亞洲中部半干旱地區(qū)森林樹木生長的下降。Liang等[11]探究了青藏高原東北部半干旱地區(qū)云杉林在全球變暖背景下的森林衰退現(xiàn)象,也認(rèn)為干旱事件導(dǎo)致了森林衰退,快速的氣候變暖對森林樹木的生長有著非常重要的影響。本研究發(fā)現(xiàn)該區(qū)域三次較為嚴(yán)重的非健康事件爆發(fā)的高峰主要集中在1930s、1970s和2000s左右。根據(jù)研究區(qū)的氣象記錄,忻州地區(qū)5—9月平均氣溫與同一時期降水量顯著負(fù)相關(guān),即研究區(qū)氣候具有高溫少雨、低溫多雨的組合模式。在中國北方地區(qū), 高溫少雨總是和干旱聯(lián)系在一起。20世紀(jì)20—30年代前后由于東亞夏季風(fēng)減弱而導(dǎo)致中國北方降水量明顯減少,而這次著名干旱事件在許多樹輪研究中也得到印證[31,33]。李強(qiáng)等[34]對山西寧武地區(qū)降水重建資料顯示,1966—1972年和1988—2003年為降水較少的時期(低于多年平均降水量413 mm),其中1972年干旱尤為嚴(yán)重。蔡秋芳等[27]對山西呂梁山氣溫變化的研究中發(fā)現(xiàn),1994—2002年的溫度是過去近170年中溫度最高的時段,這一持續(xù)高溫事件對此次干旱事件的爆發(fā)無疑貢獻(xiàn)巨大。Williams等[35]認(rèn)為干旱會導(dǎo)致氣孔關(guān)閉,光合作用降低,從而限制了形成層活動,使得樹木在生長季受到嚴(yán)重的水分脅迫。本研究發(fā)現(xiàn)的3次較大的不健康事件與歷史時期嚴(yán)重干旱事件時間相吻合,表明該地區(qū)降水減少或高溫所帶來的極端干旱是影響該區(qū)域森林健康生長的主要原因。

20世紀(jì)中葉以來的氣候變暖是不容置疑的[36]。忻州地區(qū)地形崎嶇,多為山地丘陵,水分條件較差,干旱事件發(fā)生的頻率和干旱程度在未來氣候變化背景下都可能會增加。高溫對育苗造成的嚴(yán)重危害和干旱對造林成活率和保存率的嚴(yán)重影響等都是該地區(qū)林業(yè)發(fā)展的不利因素。要針對這些不利的限制性氣候因子,選擇合理的造林季節(jié),避開高溫和干旱時期因地適宜地選擇耐旱樹種進(jìn)行栽培和經(jīng)營,實現(xiàn)林業(yè)的可持續(xù)發(fā)展。

4 結(jié)論

本文以忻州地區(qū)4個采樣點(diǎn)森林為研究對象,利用樹木年輪生態(tài)學(xué)方法重建了自1900年以來忻州森林的健康歷史,分析了森林健康狀況的時空變化特征,并對產(chǎn)生不健康事件的原因進(jìn)行深入探討。研究結(jié)果表明：過去100年中該地區(qū)森林不健康事件主要集中發(fā)生在1910—1940、1970—1987和1990—2012年期間,4個采樣點(diǎn)在20世紀(jì)30年代和世紀(jì)末不健康事件在空間上較為一致,20世紀(jì)70年代的不健康事件在寧武和岢嵐地區(qū)表現(xiàn)較強(qiáng),而在繁峙地區(qū)僅有少量樹木表現(xiàn)出非健康生長。在經(jīng)歷了3次較嚴(yán)重的不健康事件的影響,樹木能夠度過環(huán)境不利時期依然存活,表明該地區(qū)的森林具有較強(qiáng)的抵抗力或生態(tài)彈性。升溫或降水減少所引起的極端干旱是造成該區(qū)域森林非健康生長的主要原因,也是不健事件空間一致性的主要驅(qū)動力。研究結(jié)果有助于進(jìn)一步認(rèn)識區(qū)域森林健康在全球氣候背景下時間和空間上的變化特征以及其驅(qū)動因子,對科學(xué)地制定森林保護(hù)措施具有積極的現(xiàn)實意義。在以后的研究中仍需加強(qiáng)空間同步信號的鑒別能力,在更廣泛的區(qū)域和氣候背景下開展森林健康的研究,來提高我們對森林衰退的預(yù)防和治理能力。

致謝：中國科學(xué)院植物研究所樹木年輪實驗室提供樹輪樣本,張齊兵研究員在論文寫作過程中給予建議,邱紅巖老師幫助實驗工作,特此致謝。

[1] 高均凱. 森林健康基本理論及評價方法研究[D]. 北京: 北京林業(yè)大學(xué), 2007.

[2] Amoroso M M, Daniels L D, Larson B C. Temporal patterns of radial growth in declining Austrocedrus chilensisforests in Northern Patagonia: the use of tree-rings as an indicator of forest decline. Forest Ecology and Management, 2012, 265(1): 62-70.

[3] Brouwers N C, Mercer J, Lyons T, Poot, Veneklaas E, Hardy G. Climate and landscape drivers of tree decline in a Mediterranean ecoregion. Ecology and Evolution, 2013, 3(1): 67-79.

[4] Rodríguez-Catón M, Villalba R, Srur A M, Luckman B. Long-term trends in radial growth associated withNothofaguspumilioforest decline in Patagonia: integrating local-into regional-scale patterns. Forest Ecology and Management, 2015, 339: 44-56.

[5] 劉世榮, 代力民, 溫遠(yuǎn)光, 王暉. 面向生態(tài)系統(tǒng)服務(wù)的森林生態(tài)系統(tǒng)經(jīng)營: 現(xiàn)狀、挑戰(zhàn)與展望. 生態(tài)學(xué)報, 2015, 35(1): 1-9.

[6] 高廣磊, 信忠保, 丁國棟, 李叢叢, 張佳音, 梁文俊, 安云, 賀宇, 肖萌, 李文葉. 基于遙感技術(shù)的森林健康研究綜述. 生態(tài)學(xué)報, 2013, 33(6): 1675-1689.

[7] 王彥輝, 肖文發(fā), 張星耀. 森林健康監(jiān)測與評價的國內(nèi)外現(xiàn)狀和發(fā)展趨勢. 林業(yè)科學(xué), 2007, 43(7): 78-85.

[8] Innes J L. Methods to estimate forest health. Silva Fennica, 1993, 27(2): 145-157.

[9] 邵雪梅. 樹輪年代學(xué)的若干進(jìn)展. 第四紀(jì)研究, 1997, 17(3): 265-271.

[10] Liu H Y, Williams A P, Allen C D, Guo D L, Wu X C, Anenkhonov O A, Liang E Y, Sandanov D V, Yin Y, Qi Z H, Badmaeva N K. Rapid warming accelerates tree growth decline in semi-arid forests of Inner Asia. Global Change Biology, 2013, 19(8): 2500-2510.

[11] Liang E Y, Leuschner C, Dulamsuren C, Wagner B, Hauck M. Global warming-related tree growth decline and mortality on the north-eastern Tibetan plateau. Climatic Change, 2016, 134(1/2): 163-176.

[12] 姚啟超, 王曉春, 肖興威. 小興安嶺紅皮云杉年輪-氣候關(guān)系及其衰退原因. 應(yīng)用生態(tài)學(xué)報, 2015, 26(7): 1935-1944.

[13] Millar C I, Stephenson N L. Temperate forest health in an era of emerging megadisturbance. Science, 2015, 349(6250): 823-826.

[14] Whyte G, Howard K, Hardy G E S J, Burgess T I. The tree decline recovery seesaw; a conceptual model of the decline and recovery of drought stressed plantation trees. Forest Ecology and Management, 2016, 370: 102-113.

[15] Mitchell P J, Battaglia M, Pinkard E A. Counting the costs of multiple stressors: is the whole greater than the sum of the parts? Tree Physiology, 2013, 33(5): 447-500.

[16] Zhang Q B, Alfaro R I. Spatial synchrony of the two-year cycle budworm outbreaks in central British Columbia, Canada. Oikos, 2003, 102(1): 146-154.

[17] 胡海清, 趙致奎, 王曉春, 張遠(yuǎn)東. 基于樹輪火疤塔河蒙克山樟子松林火災(zāi)的頻度分析. 生態(tài)學(xué)報, 2010, 30(23): 6372-6379.

[18] 佟杰, 張永國. 忻州市城區(qū)空氣污染特征及防治對策. 忻州師范學(xué)院學(xué)報, 2007, 23(5): 63-65.

[19] Kimball K D, Keifer M. Climatic comparisons with tree-ring data from montane forests: are the climatic data appropriate? Canadian Journal of Forest Research, 1988, 18(4): 385-390.

[20] Cook E R, Kairiukstis L A. Methods of Dendrochronology: Applications in the Environmental Sciences. Dordrecht: Kluwer Academic Publishers, 1990: 51-55.

[21] Holmes R L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin, 1983, 44(3): 69-75.

[22] Cook E R, Holmes R L. Users Manual for ARSTAN. Tucson: Laboratory of Tree-Ring Research, University of Arizona, 1986.

[23] Biondi F, Waikul K. DENDROCLIM2002: a C++ program for statistical calibration of climate signals in tree-ring chronologies. Computers & Geosciences, 2004, 30(3): 303-311.

[24] Liang E Y, Shao X M, Hu Y X, Lin J X. Dendroclimatic evaluation of climate-growth relationships of Meyer spruce (Piceameyeri) on a sandy substrate in semi-arid grassland, north China. Trees, 2001, 15(4): 230-235.

[25] 戴君虎, 潘嫄, 崔海亭, 唐志堯, 劉鴻雁, 曹燕麗. 五臺山高山帶植被對氣候變化的響應(yīng). 第四紀(jì)研究, 2005, 25(2): 216-223.

[26] 李穎俊, 王尚義, ?？〗? 方克艷, 李曉嵐, 栗燕, 布文麗, 李玉晗. 蘆芽山華北落葉松(Larixprincipis-rupprechtii)樹輪寬度年表對氣候因子的響應(yīng). 生態(tài)學(xué)報, 2016, 36(6): 1608-1618.

[27] 蔡秋芳, 劉禹, 包光, 雷鶯, 孫鉑. 樹輪記錄的呂梁山地區(qū)公元1836年以來5—7月平均氣溫變化. 科學(xué)通報, 2010, 55(20): 2033-2039.

[28] Cai Q F, Liu Y. Climatic response of Chinese pine and PDSI variability in the middle Taihang Mountains, north China since 1873. Trees, 2013, 27(2): 419-427.

[29] Sun J Y, Liu Y, Sun B, Wang R Y. Tree-ring based PDSI reconstruction since 1853 AD in the source of the Fenhe river basin, Shanxi province, China. Science China Earth Sciences, 2012, 55(11): 1847-1854.

[30] 吳祥定. 樹木年輪與氣候變化. 北京: 氣象出版社, 1990: 77-234.

[31] 梁爾源, 邵雪梅, 黃磊, 王麗麗. 中國中西部地區(qū)樹木年輪對20世紀(jì)20年代干旱災(zāi)害的指示. 自然科學(xué)進(jìn)展, 2004, 14(4): 469-474.

[32] 周丹卉, 賀紅士, 李秀珍, 周春華, 王緒高, 陳宏偉. 小興安嶺不同年齡林分對氣候變化的潛在響應(yīng). 北京林業(yè)大學(xué)學(xué)報, 2007, 29(4): 110-117.

[33] Liang E Y, Liu X H, Yuan Y J, Qin N S, Fang X Q, Huang L, Zhu H F, Wang L L, Shao X M. The 1920s drought recorded by tree rings and historical documents in the semi-arid and arid areas of northern China. Climatic Change, 2006, 79(3/4): 403-432.

[34] 李強(qiáng), 劉禹, 蔡秋芳, 孫軍艷, 易亮, 宋惠明, 王雷. 山西寧武地區(qū)1686年以來年降水重建. 第四紀(jì)研究, 2006, 26(6): 999-1006.

[35] Williams A P, Allen C D, Millar C I, Swetnam T W, Michaelsen J, Still C J, Leavitt S W. Forest responses to increasing aridity and warmth in the southwestern United States. Proceedings of the National Academy of Sciences, 2010, 107(50): 21289-21294.

[36] IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK and New York, NY: Cambridge University Press, 1-1535.