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

格陵兰近千年气候变化的多尺度分析1 杜建丽，林振山，俞鸣同，张真真，刘会玉 南京师范大学地理科学学院，南京 210046 E-mail: linzhenshan@ 18 摘 要：用经验模态分解（EMD ）方法对格陵兰冰盖GISP2 冰芯古气候代用指标δ O 序列 进行分析，结果表明：北半球高纬地区气候波动具有 3、6.5、12、24 、49 、96 、213 、468 年准周期性波动，既有NESO 的影响，也受制于太阳活动的周期。近1000 年来气候变化总 趋势出现历时约490 年的中世纪暖期和历时约570 年的小冰期，其间还存在次级的冷暖期变 化；EMD 第7 ，第8 本征模函数IMF7 和IMF8 波动振幅以及EMD 趋势项在 1350 年出现明 显的转型，表明 1350 年为中世纪暖期和小冰期的分界。研究结果对北半球高纬地区近千年 气候变化的波动规律和变化分期作了较详细的划分。 关键词：格陵兰冰芯；EMD；波动周期；小冰期 0. 前言 过去 2ka 的气候变化是“过去全球变化” （ 两大国际研究计划的重要研究内容[1,2,3]，其中 1 ka 来的气候变化是现代器测资料与过去代用 指标序列衔接的关键时段,也是人类活动影响不断加剧的时期,气候冷暖变化阶段和周期以及 现代气温的快速增暖成为各国科学家关注的热题。 目前，对2ka 来气候研究主要关注多种代用指标所揭示的中世纪暖期、小冰期在不同地 区的反应，中世纪暖期、小冰期是否具有全球意义，暖期和小冰期内的次级气候波动及其峰 谷值点，中世纪暖期向小冰期转型的分界， 以及现代增暖期的变幅等热点问题。中世纪暖 期以及随后的小冰期在欧洲、北美中高纬地区的多种代用指标里有所反映。Briffa K R, Lipp J, 等研究西欧多处树轮发现，大约 950～1200AD 气候暖期对树木生长和稳定同位素变化具 有明显的影响[4] 18 ；Dansgaard W，Grove J M 和 Stuiver M,等研究冰岛和格陵兰冰芯 δ O 序列 ， ， 时指出其中具有中世纪暖期的信号[5 6 7] ；Bodri 和 Cermak 从捷克深钻孔反演的温度变化也 发现了中世纪暖期的存在[8] 。中国东部和西部对中世纪暖期的反应差别较大。王绍武等根据 史料、树木年轮、冰芯等代用资料建立了 1000～2000 AD 中国东部、西部及全国平均气温 距平序列, 认为中国东部在 1040～1130AD；1170～1260AD 存在中世纪暖期,而西部无明显 的中世纪暖期[9] 。北半球小冰期研究已有较多的报道。捷克钻孔对古地温研究表明小冰期出 现在 1400～1500AD，1600～1700AD，最冷时期为 1650AD；Stuiver M,对格陵兰冰芯 δ18O 研究表明小冰期出现在 1300～1900AD；王绍武研究中国近千年气候变化指出，中国东部小 冰期出现 1600～1690AD，1780～1870AD，西部的反应可提前至 11 世纪。 这些研究表明千年来气候变化出现中世纪暖期和小冰期现象，但是在各地区的反应有所 差异，表现在暖期和冰期的起止年代不同，气候期内的次级波动次数和峰谷极值出现的年代 也不同；中世纪暖期向小冰期气候转型的分界年代也不一致，以及千年来对气候变化影响显 著的小尺度高频波动的周期规律尚不清楚。这些研究多数从时间序列数据的直观变化中总结 规律，发现事件的突变，或作简单的数据平滑处理，缺少严密的科学方法推导非平稳序列的 数学模式，蕴涵在非平稳序列中的多尺度变化规律难以发现，阈值和临界点的确定带有较大 的主观臆断，因此众说纷纭，难以达成共识。 1本课题得到教育部博士点基金项目 座机电话号码010 “1ka-10ka 全新世东亚季风的多尺度诊断和预测”的资助。 - 1 - 经验模态分解（EMD ）方法是目前最好、最有效地进行非平稳序列多尺度分析的方法[10] 。 1998年Huang[11]提出了经验模态分解，1999年又做了一些改进[12] 。EMD可以对一个信号序 列进行不同尺度 （频率）的波动或趋势逐级分解，产生一系列具有不同特征尺度的数据序列， 所表达的物理意义明确。EMD 方法具有自适应的特性,适宜于非平稳信号的分解，而被美国 NASA宇航中心采用[10] 。 1. 研究方法 1.1 经验模态分解 EMD 简介 自然界的波谱信号往往是多种不同尺度信号的迭代结果。被检测的信号综合了不同尺度 信号的因素，使得信号波谱杂乱无章，难以分辨信号中隐含的不同尺度信号的物理意义。过 去对第四纪样品代用指标检测数据的分析，多数以目视方法、概率统计方法做简单解释；或 用简单的数据平滑处理