格陵兰近千年气候变化的多尺度分析.doc

士 j ■ M nJ ■ a M ■ jT 亠山格陵兰近千年气候变化的多尺度分析本课题得到教育部博士点基金项目 （20060139010） “ 1ka-10ka全新世东亚季风的多尺度诊断和预测” 的资助杜建丽，林振山，俞鸣同，张真真，刘会玉南京师范大学地理科学学院，南京 （210046）E-mail: linzhenshan@摘 要：用经验模态分解（EMD ）方法对格陵兰冰盖 GISP2冰芯古气候代用指标 S8O序列进行分析，结果表明：北半球高纬地区气候波动具有 3、6.5、12、24、49、96、213、468年准周期性波动，既有 NESO的影响，也受制于太阳活动的周期近 1000年来气候变化总趋势出现历时约 490年的中世纪暖期和历时约 570年的小冰期，其间还存在次级的冷暖期变化；EMD第乙第8本征模函数IMF7和IMF8波动振幅以及 EMD趋势项在1350年出现明 显的转型，表明1350年为中世纪暖期和小冰期的分界研究结果对北半球高纬地区近千年 气候变化的波动规律和变化分期作了较详细的划分关键词：格陵兰冰芯；EMD ;波动周期；小冰期0.刖言过去2ka的气候变化是 过去全球变化” （PAGES）和 气候变率与可预测性”（CLIVAR ） 两大国际研究计划的重要研究内容 ［1,2,3］，其中1 ka来的气候变化是现代器测资料与过去代用指标序列衔接的关键时段，也是人类活动影响不断加剧的时期，气候冷暖变化阶段和周期以及 现代气温的快速增暖成为各国科学家关注的热题。

目前，对2ka来气候研究主要关注多种代用指标所揭示的中世纪暖期、 小冰期在不同地区的反应，中世纪暖期、小冰期是否具有全球意义，暖期和小冰期内的次级气候波动及其峰 谷值点，中世纪暖期向小冰期转型的分界， 以及现代增暖期的变幅等热点问题中世纪暖期以及随后的小冰期在欧洲、北美中高纬地区的多种代用指标里有所反映 Briffa K R, LippJ,等研究西欧多处树轮发现，大约 950〜1200AD气候暖期对树木生长和稳定同位素变化具有明显的影响 ⑷；Dansgaard W，Grove J M 和Stuiver M,等研究冰岛和格陵兰冰芯 S8O序列时指出其中具有中世纪暖期的信号 ［5,6, 7］ ； Bodri和Cermak从捷克深钻孔反演的温度变化也发现了中世纪暖期的存在［8］中国东部和西部对中世纪暖期的反应差别较大王绍武等根据 史料、树木年轮、冰芯等代用资料建立了 1000〜2000 AD中国东部、西部及全国平均气温距平序列，认为中国东部在 1040〜1130AD ； 1170〜1260AD存在中世纪暖期，而西部无明显 的中世纪暖期［9］北半球小冰期研究已有较多的报道捷克钻孔对古地温研究表明小冰期出 现在1400〜1500AD，1600〜1700AD，最冷时期为 1650AD ； Stuiver M,对格陵兰冰芯 护0 研究表明小冰期出现在 1300〜1900AD ;王绍武研究中国近千年气候变化指出，中国东部小 冰期出现1600〜1690AD，1780〜1870AD，西部的反应可提前至 11世纪。

这些研究表明千年来气候变化出现中世纪暖期和小冰期现象， 但是在各地区的反应有所差异，表现在暖期和冰期的起止年代不同， 气候期内的次级波动次数和峰谷极值出现的年代也不同；中世纪暖期向小冰期气候转型的分界年代也不一致， 以及千年来对气候变化影响显著的小尺度高频波动的周期规律尚不清楚 这些研究多数从时间序列数据的直观变化中总结规律，发现事件的突变，或作简单的数据平滑处理， 缺少严密的科学方法推导非平稳序列的数学模式，蕴涵在非平稳序列中的多尺度变化规律难以发现， 阈值和临界点的确定带有较大的主观臆断，因此众说纷纭，难以达成共识HI田魂QP衣幵t±«■« I ■ -K经验模态分解(EMD )方法是目前最好、最有效地进行非平稳序列多尺度分析的方法 ［10］1998年Huang［11］提出了经验模态分解，1999年又做了一些改进［12］EMD可以对一个信号序 列进行不同尺度(频率)的波动或趋势逐级分解，产生一系列具有不同特征尺度的数据序列， 所表达的物理意义明确 EMD方法具有自适应的特性，适宜于非平稳信号的分解， 而被美国NASA宇航中心采用［10］1. 研究方法1.1经验模态分解(EMD)简介自然界的波谱信号往往是多种不同尺度信号的迭代结果。

被检测的信号综合了不同尺度信号的因素，使得信号波谱杂乱无章， 难以分辨信号中隐含的不同尺度信号的物理意义 过去对第四纪样品代用指标检测数据的分析， 多数以目视方法、概率统计方法做简单解释；或用简单的数据平滑处理后，判读数据的地学意义这些处理没有把数据的原本真实的物理意 义反映出来，或扭曲了数据的真实物理意义，使解释意义失真经验模态分解(EMD )方法是一种处理非平稳数据序列的全新方法，其本质是对一个 信号进行平稳化处理，将信号中不同尺度的波动或趋势逐级分解开来， 形成一系列具有不同尺度的数据序列，每一个序列称为一个本征模函数( intrinsic mode function, IMF )分量，其中不同尺度数据的物理意义得以保留， 并很好地被提取和表达 最低频率的分量(IMF , res)代表原始信号的总趋势或均值的时间序列EMD方法的处理过程相对简单，适合于对大容量数据的处理具体方法 ［10］是：找出原始数据序列X(t)所有的极大值点，其中t是时刻，t=1,…,T用三次样条函数拟合，形成原 数据序列的上包络线；相应地找出 X(t)所有的极小值点，同样用三次样条函数拟合，形成原数据序列的下包络线；上下包络线的均值为原数据序列的平均包络线 m1(t);将原始数据序列X(t)减去平均包络，得到一个去掉低频的新数据序列 h1(t)。

h1 (t)=X(t) — m1(t) (1)通常情况下，一次处理后h1(t)仍然不是一个平稳数据序列， 因此，将其重复上述处理过程k次直到所得到的平均包络趋于零为止hk(t) = h (k-1)(t) — mk(t) (2)由此，得到第一本征模函数(IMF1 )分量C1 (t):C1(t)= hk(t) ⑶IMF1分量代表原始数据序列中最高频的组分，保留了原始数据中最高频信号的物理特 征原始数据中其他尺度的信号可以用同样的方法继续提取 将原始数据序列X(t)减去第一个IMF分量C1(t)，可以得到一个去掉高频组分的差值数据序列 r1(t)对r1(t)进行上述平稳化处理过程可以得到第 2个IMF分量C2(t),如此重复下去直到最后一个差值序列 rn(t)不可再分解为止此时，rn(t)代表原始数据序列的趋势或均值的时间序列：r2(t) =r 1(t)- C2(t), - n(t) = r p-n (t)- C n(t) (4)Hua ng将这样的处理过程形象地比喻为 筛”过程最后，原始的数据序列即可由这些IMF分量以及一个均值或趋势项表示：nX(t)= ”Cj(t) + rn(t). ⑸-# -士 i ■ M nJ ■ iB M ■ *'亠山«■« K ■ * 1.2代用指标的含义和来源研究古气候变化规律和机制最主要的参数是温度变化。

在古气候变化多种代用指标中， 冰盖降水护0变化对气温变化的敏感性和时标的测定准确性相对较高 其基本原理是水汽海陆循环中，即0是地表温度的函数温度升高时，降水中 厝0较高；温度降低时，降水中 尹0 较低因此，用 即0变化的时间序列直接表达气温的变化规律是切实可行的（请引用前人文 献）格陵兰冰盖钻孔GISP2冰芯e?80变化代表了北半球高纬地区近千年气候变化的规律 对这一规律的深入研究和拓展，有助于对北半球未来气候变化的预测90年代，格陵兰冰盖研究项目取得了大量有关第四纪气候变化代用指标的数据 美国华盛顿大学第四纪同位素实验室完成的格陵兰冰盖 GISP2冰芯 护0数据序列，时间跨度从公元818〜1987年，具有每年间隔，信号分辨率高，数据连续性好，间隔均匀，数据量大，内含气 候变化信息丰富等优点，适合于做多尺度变化研究（注明数据来源：如世界数据中心 WDC的古环境数据库，或原文出处）2. 格陵兰冰盖GISP2冰芯S18O序列的EMD分析从数据的直角坐标显示（图1），几乎找不到什么规律，既看不出衣8o的多时间尺度特征 性,也看不出周期变化在时间域中的分布情况，因此，必须对数据做必要的处理31-32-33-34-35-36-37-38-39800 1000 1200 1400 1600 1800 2000图1格陵兰冰芯 目80近1000年来的趋势禾U用EMD运算软件，对格陵兰冰盖 GISP2冰芯公元818a〜1987a近千个数据的 护0序列进行分解，得到的 8个本征模函数（intrinsic mode function, IMF ）分量及其数据趋势分量res （图2,表1）。

每个IMF分量图形表示不同尺度或一个窄波段的 即0变化特征信号 -士 j ■ M nJ ■ a M ■ jT 亠山HI田魂QP衣幵t±«■« I ■ -K— . 1 ■ 1 1 1 t1 'a ... 月… 1 L a J?.1 ■ L U1 L. Ln .- «ul- . 'j J 41A.-! !. ■iJ ! !.■ -Jji 1 -'l Ji . . U . . 1,L" .'..1'：.2~""T~| ■— p, - ■"-h I = ■1000 110012001300 1400 1500 1600 1700 1800 19001000 1100120013001400 1500 16001700 18001900o-12LT Ml2O-2900ii ■ i r ・ i ・ i l ... 1 ■ :.. .... . : T M ■- V ' - w 、 ■ - •• J > "■ v f、' ■ F ■ " ■ ；V 、 ■匚1 j j -900 10001100 1200 13001400 15001600 1700 18001900o-13r Ml' ■• J* r 严、1广 Ji '•.”「一 r丄一一A 4 LA JlLS~uIII111』1 1 1』U1LJ U11*900~1000110012001300140015001600170018001900"T7~1 ~十~1 ~~1^~1 ~n-■1IILl11Ji ”一]11\ i—9001000110012001300140015001600170018001900■1"-I1I..1T1J*I1-■11||” 1ilIi1。