全新世氣候?qū)η瓿叨忍柣顒拥目赡茼憫?yīng)

尹志強 韓延本 王博

（中國科學(xué)院國家天文臺，北京 100012）

全新世氣候?qū)η瓿叨忍柣顒拥目赡茼憫?yīng)

尹志強 韓延本 王博

（中國科學(xué)院國家天文臺，北京 100012）

利用熟知的山洞石筍資料及重建的太陽黑子數(shù)獲取的時間序列，討論了長時間尺度上氣候?qū)μ柣顒幼兓目赡茼憫?yīng)。小波分析顯示氣候變化與太陽活動變化在約1000及約2000年周期上具有一致性，這可能反映了全新世時期太陽活動的調(diào)制作用。另外，小波功率譜顯示在公元前2000年附近缺失千年周期，這種突變及氣候條件的長期變化可能給許多地區(qū)帶來社會性影響。

古氣候，太陽黑子數(shù)，全新世

1 引言

眾所周知，氣候系統(tǒng)受許多因素的影響，例如地球軌道變化、太陽活動變化、大型熱帶火山噴發(fā)、地面覆蓋差異及溫室氣體等。其中，太陽活動的變化是關(guān)鍵因素之一。目前有許多有關(guān)太陽活動變化與各種氣候記錄相關(guān)性的研究成果，但涉及千年時間尺度的太陽活動驅(qū)動氣候變化的機制與該驅(qū)動下氣候系統(tǒng)響應(yīng)之間關(guān)系的研究相對較少[1-4]。

由于缺少更長時間的直接觀測或儀器記錄數(shù)據(jù)，我們只能依靠代用指標(biāo)來研究全新世時期的古氣候。例如，樹木年輪、珊瑚、鐘乳石及冰芯等代用指標(biāo)通常能很好地反映太陽活動、氣溫、濕度及季風(fēng)強度等[5-7]。也有許多學(xué)者曾嘗試從陸地及海洋代用指標(biāo)中獲取十年至百萬年時間尺度的氣候變化，研究太陽活動對氣候變化的驅(qū)動[3,7]。

對于千年周期尺度的氣候變化，格陵蘭冰芯及北大西洋的巖芯很好地記錄了全新世期間氣候變化[8-9]。許多學(xué)者也在古氣候資料，例如在Big Lake沉積物記錄[10]、加拿大[11]、北美花粉記錄[12]中找到了這些變化周期。Clemens等[13]認(rèn)為千年尺度的氣候變化是由太陽活動世紀(jì)周期外差得來的。Bütikofer等[14]研究表明，全新世氣候記錄的典型特征是寬波段的準(zhǔn)周期變化，而不具有某一顯著的譜峰。

本文對收集到的一些數(shù)據(jù)做了小波分析，以進一步研究全新世千年尺度的氣候變化，并討論存在的周期振蕩及氣候?qū)μ柣顒涌赡茼憫?yīng)的時變特征。

2 數(shù)據(jù)及分析方法

2.1 重建的太陽黑子數(shù)

我們可以使用宇宙成因核素（如14C，10Be等）含量變化的高頻分量作為太陽活動變化的代用指標(biāo)。宇宙成因核素是由宇宙線輻射與大氣分子相互作用形成的，其含量受太陽活動強度變化的調(diào)制。Solanki等[15]根據(jù)中緯度地區(qū)樹木年輪10年采樣間隔的高精度14C資料分析獲取的δ14C記錄，給出了覆蓋全新世的太陽黑字?jǐn)?shù)（SSN）的重建結(jié)果。該數(shù)據(jù)在公元1600年之后的結(jié)果與實際觀測有很好的一致性。雖然10年采樣間

隔的數(shù)據(jù)不能用來分析太陽活動的高頻波動，但足以分析百年及更長的周期性振蕩。

2.2δ18O記錄

從冰芯、有孔蟲殼或石筍中獲取的穩(wěn)定的氧同位素比值（18O/16O）記錄被廣泛用于古氣候的研究。千年尺度變化的分析對數(shù)據(jù)的要求更為嚴(yán)格，我們需要更長時間跨度及更高分辨率的信息。根據(jù)Fleitmann等[16]的分析報告，選取了兩個高分辨率及較長時間跨度的低緯度δ18O記錄。其中，董哥（Dongge）洞（25°17′N，108°5′E）的數(shù)據(jù)被認(rèn)為記錄了過去9千年亞洲季風(fēng)的信息，數(shù)據(jù)采樣間隔從1年至10年不等，平均采樣間隔為4年。該記錄與大氣14C記錄的十年至百年尺度的互相關(guān)分析顯示，該頻率的一些季風(fēng)變化被認(rèn)為是由太陽活動的變化引起的[1,17]。為了便于進一步分析，我們利用分段3次Hermite方法對該數(shù)據(jù)做了內(nèi)插處理，得到5年間隔的數(shù)據(jù)序列。另一個從阿曼北部與南部及也門的4個山洞的全新世石筍中獲取的高分辨率δ18O數(shù)據(jù)，其提供了沿12°—23°N方向的沉積作用的詳細波動信息[16,18-19]。我們同樣利用分段3次Hermite方法以10年分辨率對該數(shù)據(jù)做了內(nèi)插處理。下文中這兩個代用指標(biāo)分別簡稱為Dongge及Qunf。圖1繪出了所有數(shù)據(jù)的曲線，以展示強或弱的突變事件。

2.3 小波分析

為了估計一些數(shù)據(jù)頻率的非平穩(wěn)性，Morlet提出了小波變換。小波是一類用于對一給定的函數(shù)在空間及時間域上做局部化分析的函數(shù)。Morlet小波將三角函數(shù)與高斯函數(shù)相結(jié)合，通過時頻域上的小波濾波，將時間序列的信號與噪聲分離開。這樣可以有效消除噪聲的影響[20-22]。我們在研究中利用小波變換的方法對上述代用指標(biāo)的時變信息在時間及頻率域上分別進行了估計。

圖1 SSN（紅）、Dongge（藍）及Qunf（黑）時間序列的變化曲線

3 結(jié)果及討論

圖2給出了SSN，Dongge及Qunf數(shù)據(jù)序列的小波變換結(jié)果。由于本文主要對百年及千年尺度長期波動的變化特征做研究，故把小波分析的周期范圍選為64～4096年。小波圖（圖2a，2c，2e）中，背景置為白色，由藍色至紅色表示能量由低到高的變化。利用時間序列的協(xié)方差對小波功率做了歸一化處理，以給出功率相對白噪聲的大小。圖2b，2d，2f分別給出了3個序列的全局小波譜。

圖2a顯示全新世期間SSN時間序列具有明顯的長周期波動（包括約220，517，1000及2000年等）。雙世紀(jì)周期的變化在190～220年之間，千年周期長度的變化在900～1100年之間，2000年周期長度的變化幅度為2040～2400年。圖2b的全局小波變換顯示所有的周期都在置信水平以上。圖2c—d顯示Dongge時間數(shù)據(jù)具有類似的特征。Dongge序列的主要周期約為210，567，1000，2320年。Qunf序列的分析顯示其幾個顯著的周期約為：1000，1820及2400年（如圖2e—f）。幾乎所有的δ18O的主要變化周期與SSN的變化周期接近。

全新世的十年與百年尺度上的10Be及14C中的約90年的Gleisberg周期及約200年的Suess周期以及其他太陽活動的記錄，同樣很好地在環(huán)境代用指標(biāo)記錄中展現(xiàn)出來[7,23-28]。在百年及千年尺度上的小波分析顯示SSN的約200年周期幅度較大的時段主要為公元前約8000—6000年間、公元前約4000—2700年間、公元前約700—50年間，以及公元1500年前后的幾個時段內(nèi)。

Dongge記錄同樣顯示出約200年的周期（圖2c）。雖然該記錄的譜分析顯示出顯著的208年及86年的太陽周期，且認(rèn)為其受太陽驅(qū)動活動的影響[17]，但局部小波功率譜分析顯示該周期在Dongge與SSN數(shù)據(jù)中時間不一致（圖2a與2c）。圖2e顯示了Qunf序列的時變特征，我們發(fā)現(xiàn)在置信水平之上并不存在約200年周期（圖2f）。

上述的分析表明，約1000年及2000年的周期波動是SSN，Dongge及Qunf序列的更長期的變化特征，其都在置信水平之上。與其他波動相比，這些振蕩似乎更加穩(wěn)定。尤其是3個序列的約1000年周期幾乎在整個時間跨度上一直存在。另外，Dongge的周期長度在公元前600年附近較短，在公元前2000年之后稍微變長。對Qunf記錄，在公元前2000年之前，該周期長度更為穩(wěn)定，隨后稍微變短。這3個時間序列的局部小波功率譜顯示，在公元前2000年附近，約1000年周期有明顯的下降（圖2）。Debret等[29]給出了14C含量[10]

及10Be[30]的類似分析結(jié)果，主要的約1000年周期同樣在5000-0 B.P.呈現(xiàn)衰減。

為了展示這些序列的千年尺度波動，我們對每個時間序列做了歸一化尺度平均小波功率譜，圖3繪出了相應(yīng)的結(jié)果（最大值為1，最小值為0）。千年尺度變化的小波尺度約為900～1100年。

由圖3可以看到，Dongge及Qunf在千年尺度上有明顯的相關(guān)性。公元前7000—1000年間，Dongge序列的波動幅度與Qunf非常一致。公元前8500—7000年間，SSN的千年準(zhǔn)周期幅度相對較大，隨后變小。雖然Dongge及Qunf的小波譜的平均能量似乎與SSN的相關(guān)性較差，它們在公元前2000年附近的最小值可能與上述的SSN千年尺度缺失相關(guān)。

連續(xù)小波變換結(jié)果表明，在過去的11000年期間，連續(xù)的千年尺度振蕩周期在1000～2000年左右。它們與在北大西洋發(fā)現(xiàn)的全新世約1500年周期不同[9-10]。另外，除Dongge與Qunf記錄的夏季日照的趨勢性變化外，兩數(shù)據(jù)記錄的許多突變事件也與14C序列相關(guān)[1,16,31]，我們的分析認(rèn)為顯著的約1000年波動以及由這兩個記錄獲取的全新世約2000年波動可能來自于太陽活動的驅(qū)動。Debret等[29]同樣給出了北大西洋記錄的類似結(jié)果。

圖2 SSN、Dongge及Qunf數(shù)據(jù)序列的小波變換。（a）、（c）、（e）分別為SSN、Dongge、Qunf時間序列的局部小波功率譜；黑色虛線為Morlet小波變換端部效應(yīng)的邊界，該邊界以下的小波譜的可靠性較差；黑色實線為95%置信水平的置信線。（b）、（d）、（f）分別為SSN、Dongge、Qunf時間序列的小波全局功率譜，虛線為以白噪聲為背景的95%置信水平的置信線

圖3 SSN、Dongge及Qunf時間序列的歸一化尺度平均小波功率譜

4 結(jié)論

通過再分析的中低緯度山洞石筍的δ18O數(shù)據(jù)序列，小波分析顯示全新世氣候的千年尺度的變化至少存在約一千年和約兩千年的兩個周期波動。我們的結(jié)果表明這兩個周期可能由太陽活動驅(qū)動。雖然譜分析發(fā)現(xiàn)δ18O記錄呈現(xiàn)出顯著的太陽活動百年峰值[17]，然而小波分析顯示它們在同一時間內(nèi)與太陽活動周期并沒有明顯的相關(guān)性。從這一結(jié)果來看，全新世期間在千年尺度上，反應(yīng)亞洲季風(fēng)的δ18O代用指標(biāo)，ITCZ或印度夏季季風(fēng)的平均緯度位置是由太陽活動驅(qū)動的。與百年尺度的變化相比，千年尺度可能是氣候系統(tǒng)變

化的更基本的周期，這種波動很少受如復(fù)雜地形，地面覆蓋，火山驅(qū)動及溫室氣體等其他因素的影響[4]?？紤]到在中緯度[32]及高緯度[29]各種氣候代用指標(biāo)發(fā)現(xiàn)的千年尺度太陽驅(qū)動[33]，我們認(rèn)為這種影響可能不是一種局部特征，而是廣泛地存在于整個氣候系統(tǒng)。

致謝:感謝Solanki教授提供重建的太陽黑子數(shù)序列，同時感謝NOAA古氣候項目提供替代數(shù)據(jù)。小波變換程序由C. Torrence及G. Comp開發(fā)，可通過如下鏈接獲?。篽ttp://paos.colorado.edu/research/wavelets/。

[1]Wang Y, Cheng H, Edwards R L, et al. The Holocene Asian monsoon: links to solar changes and North Atlantic climate. Science, 2005, 308(5723): 854-857.

[2] Rind D. The Sun’s role in climate variations. Science, 2002, 296(5568): 673-667.

[3]Bard E, Frank M. Climate change and solar variability: what’s new under the Sun. Earth and Planetary Science Letters, 2006, 248(1-2): 1-14.

[4]Wanner H, Beer J, Bütikofer J, et al. Mid- to Late Holocene climate change: an overview. Quat Sci Rev, 2008, 27(19-20): 1791-1828.

[5]Beer J, Mende W, Stellmacher R. The role of the sun in climate forcing. Quat Sci Rev, 2000, 19(1-5): 403-415.

[6]Mayewski P A, Rohlingb E E, Stager J C, et al. Holocene climate variability. Quaternary Research, 2004, 62(3): 243-255.

[7]Versteegh G J M. Solar forcing of climate 2: evidence from the past. Space Science Reviews, 2005, 120(3-4): 243-286.

[8] O’Brien S R, Mayewski P A, Meeker L D, et al. Complexity of Holocene climate as reconstructed from a Greenland ice core. Science, 1995, 270(5244): 1962-1964.

[9]Bond G, Showers W, Cheseby M, et al. A Pervasive millennialscale cycle in North Atlantic Holocene and glacial climates. Science, 1997, 278(5341): 1257-1266.

[10]Bond G, Kromer B, Beer J, et al. Persistent solar in fl uence on North Atlantic climate during the Holocene. Science, 2001, 294(5549): 2130-2136.

[11]Cumming B F, Laird K R, Bennett J R, et al. Persistent millennialscale shifts in moisture regimes in western Canada during the past six millennia. Proceedings of the National Academy of Sciences, 2002, 99(25): 16117-16121.

[12]Viau A E, Gajewski K, Fines P, et al. Widespread evidence of 1500 yr climate variability in North America during the past 14000 yr. Geology, 2002, 30: 455-458.

[13]Clemens S C. Millennial-band climate spectrum resolved and linked to centennial-scale solar cycles. Quat Sci Rev, 2005, 24(5-6): 521-531.

[14] Bütikofer J. [Online] Available:http://www.giub.unibe.ch/klimet/ docs/diplom_jbuetikofer.pdf (2014/1/1).

[15]Solanki S K, Usoskin I G, Kromer B, et al. Unusual activity of the Sun during recent decades compared to the previous 11,000 years. Nature, 2004, 431: 1084-1087.

[16]Fleitmann D, Burns S. J, Mangini A, et al. Holocene ITCZ and Indian monsoon dynamics recorded in stalagmites from Oman and Yemen (Socotra). Quat Sci Rev, 2007, 26(1-2): 170-188.

[17]Dykoski C A, Edwards R L, Cheng H, et al. A high-resolution, absolute-dated Holocene and deglacial Asian monsoon record from Dongge Cave, China. Earth and Planetary Science Letters, 2005, 233(1-2): 71-86.

[18]Neff U, Burns S J, Mangini A, et al. Strong coherence between solar variability and the monsoon in Oman between 9 and 6 kyr ago. Nature, 2001, 411: 290-293.

[19]Fleitmann D, Burns S J, Mudelsee M, et al. Holocene Forcing of the Indian Monsoon Recorded in a Stalagmite from Southern Oman. Science, 2003, 300(5626): 1737-1739.

[20]Kumar P, Foufoula-Georgiou E. Wavelet analysis for geophysical applications. Rev Geophys, 1997, 35(4): 385-412.

[21]Torrence C, Compo G. A Practical guide to wavelet analysis. Bull Amer Meteor Soc, 1998, 79(1): 61-78.

[22]Yin Z Q, Ma L H, Han Y B, et al. Long-term variations of solar activity. Chin Sci Bull, 2007, 52(20): 2737-2741.

[23]Mergentaler J. On the long-period changes of solar activity. Acta Astronomica, 1984, 34(2): 263-267.

[24]Ogurtsov M G, Nagovitsyn Y A, Kocharov G E, et al. Long-period cycles of the Sun’s activity recorded in direct solar data and proxies. Sol Phys, 2002, 211(1-2), 371-394.

[25]Ogurtsov M G, Jungner H, Kocharov G E, et al. On the link between northern fennoscandian climate and length of the quasieleven-year cycle in galactic cosmic-ray flux. Sol Phys, 2003, 218(1-2): 345-357.

[26]Raspopov O M, DergachevV A, Kolstr?m T. Periodicity of climate conditions and solar variability derived from dendrochronological and other palaeoclimatic data in high latitudes. Palaeogeogr Palaeoclimatol Palaeoecol, 2004, 209(1-4): 127-139.

[27]Velasco V M, Mendoza B. Assessing the relationship between solar activity and some large scale climatic phenomena. Advances in Space Research, 2008, 42(5): 866-878.

[28]Ma L H, Vaquero J M. Is the Suess cycle present in historical naked-eye observations of sunspots? New Astronomy, 2009, 14(3): 307-310.

[29]Debret M, Bout-Roumazeilles V, Grousset F, et al. The origin of the 1500-year climate cycles in Holocene North-Atlantic records. Clim Past, 2007, 3(4): 569-575.

[30]Vonmoos M, Beer J, Muscheler R J. Large variations in Holocene solar activity: constraints from 10Be in the Greenland Ice Core Project ice core. Journal of Geophysical Research: Space Physics, 2006, 111(A10): A10105.

[31]Wang Y, Cheng H, Edwards R L, et al. Millennial- and orbitalscale changes in the East Asian monsoon over the past 224,000 years. Nature, 2008, 451: 1090-1093.

[32]Nederbragt A J, Thurow J. Geographic coherence of millennialscale climate cycles during the Holocene. Palaeogeography Palaeoclimatology Palaeoecology, 2005, 221(3-4): 313-324.

[33]Usoskin I G, Solanki S K, Kovaltsov G A. Grand minima and maxima of solar activity: new observational constraints. A&A, 2007, 471(1): 301-309.

Possible Response of Climate to Solar Activity on Millennial Scale during the Holocene

Yin Zhiqiang, Han Yanben, Wang bo
(National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012)

Based on well-known series from land-based caves and reconstructed sunspot numbers, the nature of climatic response to solar variability is discussed over a long-time scale. Wavelet analysis reveals coherence between the climatic and solar oscillations with periods of ～1000- and ～2000-yr that may ref l ect a modulation of solar activity during the Holocene. Moreover, the wavelet power spectra present that ～1000-yr cycle of these series is absent around ～2000 BC and this abrupt change and prolonged changes in climatic conditions might have brought about major societal challenges to many regions.

paleoclimate, Sun sunspots, Holocene

10.3969/j.issn.2095-1973.2014.04.002

2014年1月1日；

2014年3月14日

尹志強（1979—），E-mail: yinzhq@bao.ac.cn

資助信息：國家重大科學(xué)研究計劃項目（2012CB957801）