巴丹吉林沙漠南緣191 a來年降水量重建及初步分析

張永香，於琍，尹紅

（國家氣候中心，北京 100081）

巴丹吉林沙漠南緣191 a來年降水量重建及初步分析

張永香，於琍，尹紅

（國家氣候中心，北京 100081）

沙漠及其邊緣地區生態環境脆弱，對氣候變化敏感。但沙漠地區有限的森林資源限制了區域百年到千年尺度上的歷史氣候變化研究。利用采自巴丹吉林南緣的青海云杉年輪寬度資料，重建了區域近191 a（1815—2005年）來的年降水（前一年7月至當年8月的總降水量）變化序列。重建的相關系數是0.636，方差解釋量為40.4%，調整自由度后的解釋方差R2adj為0.392。重建結果穩定可靠。分析區域過去年降水變化結果可見，19世紀該區域干濕變化頻繁，20世紀前半段主要以干旱為主，干濕轉變較少。20世紀20年代的干旱事件在巴丹吉林南緣的干旱持續時間更長。周期分析的結果表明，區域年降水量變化有2 a、4 a、64 a等周期。

巴丹吉林沙漠南緣；樹木年輪；年降水量重建

全球變暖已成為不爭的事實。近期發布的IPCC第五次評估報告明確指出[1]，近一個世紀以來的全球變暖一半以上由人類活動造成。這種由于人類活動引起的氣候系統變化是史無前例的，同時這一變化對生態系統的影響也極為深刻。最新的模式預估結果表明未來全球變暖仍將持續，而隨著變暖的加劇，對不同區域的影響也不盡相同。未來的中緯度地區，干濕兩極分化將更為嚴重，濕潤的區域將更為濕潤而干旱的區域將更加干旱[2]。

巴丹吉林沙漠是我國四大沙漠之一，地處對氣候變化響應敏感的季風邊緣區，對氣候變化極為敏感[3]。已有的研究表明近半個世紀以來，隨著變暖的加劇巴丹吉林沙漠周邊地區諸如民勤綠洲、額濟納綠洲等土地荒漠化不斷加劇，沙塵天氣頻發，生態環境不斷惡化，沙漠腹地的湖泊也出現不同程度的縮小[4，5]。在此區域展開氣候變化歷史研究有助于更好地理解周邊地區土地沙漠化和腹地湖泊退縮的氣候背景，為沙漠氣候資源利用和周邊地區沙漠化治理提供理論依據，對當地生態環境保護和社會經濟發展具有重要的現實意義。

目前關于該區域的研究主要集中于地質歷史時期的氣候變遷[6，7]。樹木年輪資料是研究過去千年和百年尺度氣候變化的可靠代用資料，但沙漠內有限的森林分布限制了基于此手段對區域百年尺度的氣候變化研究。也有學者借助沙漠綠洲中的胡楊做相關研究，但結果表明，該樹種的徑向生長更多反映的是地下水的變化[8]。沙漠以南的祁連山區，森林資源豐富。已有的研究表明，利用該地的針葉樹種可以重建當地及鄰近區域歷史時期的氣候變化[9-11]。因此，本研究嘗試利用采自沙漠南緣的青海云杉重建區域的氣候變化，以期能為研究該地區氣候變化，環境保護以及社會經濟發展規劃提供百年尺度的氣候背景。

1 樹輪樣本的采集和年表建立

本研究所用資料是2006年采自巴丹吉林南緣東大山樣點（DDS）的青海云杉的寬度資料（圖1）。樣品的采集選擇了受人類活動干擾較少的青海云杉生長的最上限處，平均海拔約2 826 m，林地相對較開闊。所有樣芯均是借助生長錐在健康的、沒有明顯機械傷害的成年青海云杉采得。考慮到樹木與氣候的關系，在采樣中盡量選擇樹齡較大的樹木，坡面、海拔差異較小的范圍內采集樣芯。共采集28棵樹34個樣芯。

在樣品的前處理過程中，依據國際樹輪年代學的基本處理方法對所有樣芯依次進行固定和晾干。然后將固定好的樣芯依次用由粗到細不同大小粒度的干砂紙打磨。所用最細的砂紙粒度為600目，打磨到在顯微鏡下能夠清晰分辨細胞大小為止。再用骨架示意圖法對所有采集的樹芯進行了初步輔助交叉定年[12]，接著采用精度為0.01 mm的量測儀對所有樣芯進行輪寬測量，之后用COFECHA程序[13]對定年和量測結果進行檢驗。最后，在確定定年無誤后，所有的序列用負指數函數去除其生長趨勢，對負指數函數不能去趨勢的個別樣芯采用2/3步長的樣條函數去除生長趨勢。所有去趨勢序列通過雙權重平均方法借助ARSTAN程序建立3個樹輪年表（標準年表，差值年表和自回歸年表）。其中差值年表在1815年達到穩定通過分析，最終采用差值年表重建區域近200 a來的氣候變化歷史（圖2）。

2 氣象資料和區域氣候分析

本文所用的氣象資料均來自中國氣象科學數據集。在選擇氣象站點時，主要考慮巴丹吉林南緣各氣象站點的資料情況，并遵從與采樣點的距離近、環境差異較小，站點受城市化影響較小[14]資料可靠性好的選點原則，最終選取了6個氣象站點（鼎新（DX）、阿拉善右旗（ALSYQ）、張掖（ZY）、山丹（SD）、永昌（YC）、民勤（MQ））的資料作為主要的氣象資料來源。借助Person相關分析，計算了樣點附近各站點記錄的溫度降水的空間一致性。相關分析的結果顯示，除距離相對較遠的DX站外，各站的記錄具有較強的空間一致性。就前一年7月至當年8月的降水而言，位于巴丹吉林沙漠內部的ALSYQ站的年降水與SD站相關可達0.742，意味著利用SD站資料重建的年降水可以反映沙漠南緣的水分條件波動（表1）。本文最終選取了與采樣點氣候相對一致的SD站資料與樹輪寬度資料進行氣候重建。從SD站記錄的多年平均氣候狀況（圖3）來看，該區域氣候具有明顯的大陸性氣候特征，冬季嚴寒干燥，雨熱同期，年平均氣溫為6.4℃，7月平均最高氣溫可達20.6℃，1月平均最低氣溫為-10.1℃。年平均降水量為200.2 mm，5—9月的降水占全年的83.5%。

首先，利用Person相關分析計算了差值年表與SD站的單月氣溫與降水的關系。相關計算結果顯示，該樣點樹輪徑向生長與當年6月氣溫顯著負相關（-0.303），與當年1月和前一年7月8月的降水顯著正相關（0.273，0.274，0.289），其他月份均不顯著。這說明水分條件是限制當地樹木徑向生長的主要限制因子。利用前一年7月至當年8月（P7C8）的降水總和與樣點差值年表計算結果表明，該樣點的樹輪寬度資料能夠很好地記錄區域年降水的變化。其中與SD站的年降水相關可達0.636，與ALSYQ相關可達0.576。這表明DDS樣點的樹輪寬度資料可以用于重建并能反映巴丹吉林南緣歷史時期區域年降水的變化。

3 巴丹吉林南緣年降水量重建

利用1953年以來SD站年降水量資料（P7C8）和DDS樣點的寬度年表進行了回歸分析，并建立了重建方程，最終得到轉換函數：

式中：Pp7c8為該地區上年7月到當年8月的總年降水量重建值，DDSres為樹輪差值年表。該轉換函數的單相關系數是0.636，方差解釋量為40.4%，調整自由度后的解釋方差R2adj=0.392，F2，97=33.89，超過0.000 1的顯著水平，由該方程可以重建出1815—2005年巴丹吉林沙漠南緣上年7月到當年8月的年降水量（圖4）。對比重建值和實測值可見（圖5），重建值在高低頻上均能很好地反映實測降水的變化。

4 分析與討論

祁連山區的樹輪寬度資料已被廣泛用于歷史時期氣候變化研究。陳峰等[11]利用研究區附近樣點重建的過去兩百多年的前一年8月到當年6月的降水。對比二者可見，兩條重建序列在典型干旱事件如20世紀20年代的干旱事件上趨勢相近（圖5）。分析結果表明，本研究重建序列均值為238.1，方差為605.97，而陳峰等重建序列的均值為153.2，方差為389.15。本研究的重建結果在降水多寡年的變幅和時間上均更為清晰。除了樹輪記錄本身的敏感性差異之外，本文重建的降水總量包含了兩年的降水最大月的降水量，因此其變化更明顯。

為了更好理解過去區域降水的趨勢變化，利用t檢驗[15]，采取10 a的步長檢測了區域過去191 a以來的年降水量的變化趨勢（圖4b）。通過無信息先驗分布的貝葉斯方法計算了各突變點的顯著性[16]。突變點的正負表征降水的增加或減少。趨勢分析的結果表明，19世紀區域整體為偏濕，干濕轉換相對頻繁，而20世紀則整體呈干旱趨勢，干旱時期較長，干濕轉換緩慢。19世紀存在四次顯著的干濕轉換（1839年變濕，1849年變干，1867年變濕，1876年變干，1887年變濕）（圖4）。主要干旱時段與《西北災荒史》[17]中記錄的典型期基本一致，如1821—1847年，1857—1871年，1877—1884年。20世紀區域的干旱時期較長，且屬于長期漸變的干旱趨勢，其中1954年有變濕趨勢，1971年變干，2002年變濕，但這幾個干濕轉換并不顯著。整個20世紀中偏干的時期從世紀初一直持續到世紀中葉。其中西北地區最大持續時間最長的干旱事件，20世紀20年代在本重建序列中反映的最為清晰。相比于柴達木[18]和祁連山中東段[19]重建的降水結果，20世紀20年代的干旱在巴丹吉林沙漠南緣持續時間更長。

利用小波分析對區域年降水進行了周期分析（圖6）。周期分析的結果表明，區域年降水量存在2 a、4 a、64 a和百年尺度上的震蕩。其中，4 a期的周期與袁林等[17]研究的甘青寧地區存在4～6 a的干旱周期相一致。長周期信號較強主要在20世紀，短周期在近期也有所體現。

[1]Stocker T F,Qin,D,Plattner G-K,et al.Climate change 2013：The physical science basis.Intergovernmental Panel on Climate Change,Working Group I Contribution to the IPCC Fifth Assessment Report（AR5）2013.

[2]Seneviratne S I,Nicholls N,Easterling D,et al.Changes in climate extremes and their impacts on the natural physical environment.Managing the risks of extreme events and disasters to advance climate change adaptation，2012：109-230.

[3]魏文壽.現代沙漠對氣候變化的響應與反饋：以古爾班通古特沙漠為例[J].科學通報，2000，45（6）：636-641.

[4]白美蘭，沈建國，裴浩，等.氣候變化對沙漠化影響的評估[J].氣候與環境研究，2002，7（4）：457-464.

[5]尚可政，董光榮，王式功，等.我國北方沙區氣候變化對全球變暖的響應[J].中國沙漠，2001，21（4）：387-392.

[6]董光榮，高全洲，鄒學勇，等.晚更新世以來巴丹吉林沙漠南緣氣候變化[J].科學通報，199，40（13）：423-427.

[7]楊小平.近3萬年來巴丹吉林沙漠的景觀發育與雨量變化[J].科學通報，2000，45（4）：428-434.

[8]劉普幸，陳發虎，勾曉華，等.額濟納旗近100 a來胡楊年表的建立與響應分析[J].中國沙漠，2005，764-768.

[9]王亞軍，陳發虎，勾曉華.利用樹木年輪資料重建祁連山中段春季降水的變化[J].地理科學，2001，21（4）：373-376.

[10]康興成，程國棟，陳發虎，等.祁連山中部公元904年以來樹木年輪記錄的旱澇變化[J].冰川凍土，2003，25（5）：37-44.

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

[12]Stokes M A,Smiley T L.1968.An Introduction to Tree Ring Dating.Chicago:the University of Chicago Press.

[13]Holmes R L，1983.Computer-assisted quality control in tree ring dating and measurement.Tree Ring Bull 43：69-78.

[14]白虎志，任國玉，張愛英，等.城市熱島效應對甘肅省溫度序列的影響[J].高原氣象，2006，25（1）：90-94.

[15]Rodionov,S.A sequential algorithm for testing climate regime shifts.GEOPHYSICAL RESEARCH LETTERS, 2004：31,L09204,doi:10.1029/2004GL019448.

[16]Carlin B P,Louis T A,Bayes and empirical Bayes methods for data analysis.Statistics and Computing 1997，7：153-154.

[17]袁林.西北災荒史[M].蘭州：甘肅人民出版社，1994.

[18]邵雪梅，黃磊，劉洪濱，等.樹輪記錄的青海德令哈地區千年降水變化[J].中國科學：D輯，2004，32（2）：145-153.

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

Annual Precipitation Reconstruction over Last 191 Years at the South Edge of Badain Jaran Desert Based on Tree Ring Width Data

ZHANG Yongxiang，YU Li，YIN Hong
（National Climate Center，Beijing 100081，China）

Desert and its surrounding area are vulnerable and sensitive to climate change.Due to the limited resource of forest,few paleo-climate researches has been conducted at Badain Jaran Desert at century and millenary scale.Based on tree ring width data sampled at south edge of Badain Jaran,a 191 annual precipitation has been reconstructed.The annual precipitation is the total precipitation from previous July to current year August.The correlation between tree ring chronology and observed annual precipitation is 0.636.The reconstruction is reliable and stable and could explain variance 40.4%,adjusted 0.392.During last 200 years,there were frequent dry and wet switches during the 19th century,but less abrupt change in the 20th century.Drought is the domain condition in the first half part of the 20th century.The 1920s drought event is a long lasting and extreme drought period at the south edge of Badain Jaran Desert.

Badain Jaran Desert；tree ring width；annual precipitation reconstruction

P467

B

1002-0799（2015）01-0012-05

張永香，於琍，尹紅.巴丹吉林沙漠南緣191 a來年降水量重建及初步分析[J].沙漠與綠洲氣象，2015，9（1）：12-16.

10.3969/j.issn.1002-0799.2015.01.003

2014-10-28；

2014-11-22

氣象行業專項（GYHY201206014）和沙漠氣象科學研究基金共同資助。

張永香（1979-），女，副研究員，主要從事歷史氣候與氣候變化研究。E-mail：zhangyx@cma.gov.cn