近四百年北半球气温变化的分析emd方法的应用-中国气象局广州.pdf

第20 卷 第1 期 热 带 气 象 学 报 Vol.20, No. 1 2004 年02月 JOURNAL OF TROPICAL METEOROLOGY Feb.,2004 文章编号：1004-4965 （2004 ）01-0090-07 近四百年北半球气温变化的分析： EMD方法的应用 1 2 林振山 ， 汪曙光 （1.南京师范大学地理科学学院，江苏 南京 210097；2.南京大学大气科学系，江苏 南京 210093） 摘 要：介绍了时间序列的经验模态分解(EMD)方法，将EMD方法应用于近400年北半 球重建平均气温资料的分析，从气候时间序列中提取气候信号中各个尺度的变化，进而分 析气候的周期性振荡在20世纪的20年代变暖、40年代变冷及70年代变暖的可能作用。 关 键 词：经验模态分解； 气候变率； 全球变暖 中图分类号：P423.3.4 文献标识码：A 1 引 言 传统的气候时间序列的谱分析方法多建立在傅立叶谱分析的基础上，傅立叶谱分析的 优点是可以在频域内获得非常高的分辨率，但在时域内毫无分辨率可言。小波分析尽管在 保留信号时域特征的同时也具有一定的频域分辨率，但是海森堡测不准原理限制了其时域 [1,2] 和频域精度的提高，同时会产生很多虚假的谐波 。希尔伯特变换与小波变换相比，可以 获得更精细的时频局部特征，却不适合宽波段信号。同时，气候系统具有复杂的非线性特 征，需要分析的气候时间序列也多是非线性、非平稳的数据，各种经典的时间序列分析方 法遇到非线性、非平稳的数据就会显得力不从心，或者得到的结果没有清晰的物理含义。 [1] Huang N E于1998年提出了希尔伯特－黄变换 （Hilbert-Huang Transform, HHT），这 种变换非常适合具有非线性和非平稳特征的气候时间序列。HHT方法将信号作经验模态分 解（empirical mode decomposition, EMD ），即逐级进行平稳化处理，把不同周期的波动从 原信号中分离出来，并且该波动是平稳的，最后得到的是趋势分量，不同尺度的波动被定 义成为本征模态函数(Intrinsic Mode Function, IMF)，IMF波动分量具有显著的缓变波包的特 性，而趋势分量则是单调函数或者均值函数。不同的IMF分量是平稳信号，具有非线性特征， 收稿日期：2002-08-25；修订日期：2003-03- 10 基金项目：国家“十·五”“211”重点学科建设项目“不同时空尺度环境破坏和建设”资助 作者简介：林振山（1955-），男，福建人，教授，主要从事气候、环境、生态领域的研究。 1 期 林振山等：近四百年北半球气温变化的分析：EMD方法的应用 91 其缓变波包特征意味着不同特征尺度波动的波幅随时间变化，因而也具有时间上的局域化 特征。对IMF分量进行希尔伯特变换，得到的希尔伯特谱和小波谱具有近似的特征，但却提 供了更加清晰详细的局部细节特征。从时频分析的角度，希尔伯特－黄变换方法的最终结 果是得到在时域和频域都有较高分辨率特征的谱图。目前，EMD方法已经成功应用于湍流、 [3,4] 地震及大气科学等研究领域 。 本文将EMD 方法应用于北半球四百年平均气温资料的分析，得到一系列IMF 分量分 别表示各种不同尺度上的变化，比较分析大尺度的变化特征，得到一些有意义的结论。 2 经验正交分解方法和希尔伯特－黄变换方法 经验模态分解（EMD ）方法是一种全新的处理非平稳数据序列的方法，EMD 方法将时 间信号X (t) 分解成一系列本征模态函数IMF ，每个IMF 分量具有如