ZDCF+在气候时间序列上应用及MPT+成因一些探讨.pdfVIP

ZDCF 在气候时间序列上的应用及MPT 成因的一些探讨 沈宏 邹宛彤 陈星 （南京大学大气科学学院，南京，210093） 摘 要 互相关函数（Cross-Correlation Function，CCF）作为一种统计工具，已在大气科学的各个研究领 域取得了广泛应用。首次尝试将Tal Alexander （1997）在天文上提出的ZDCF（Z-transformed Discrete Correlation Function）方法应用在气候时间序列CCF的计算上。利用时间序列模型对ZDCF和传统的插值方 法进行了比较，并使用ZDCF研究了0.9Ma前后的古气候转型问题。用具有一定复杂气候特征的第四纪的气候 数据作为时间序列模型对ZDCF与插值法进行对比。结果显示，与原始模型数据CCF比较（1）ZDCF可以显示 插值方法由于数据稀疏而无法体现出来的小周期相关性， （2）ZDCF和插值法的相关系数极值均有一定的 相位差异，但一般情况下ZDCF显示的相位更可靠，（3）关于相关系数的大小，ZDCF普遍大于插值法而接近 模型的真实CCF值，（4）对于稀疏数据，插值法得出的相关系数对所插入值的数量不敏感，（5）基于插值法 的CCF的绝对相关性判断很难实现，而ZDCF可以实现。另外，ZDCF作为一种傅里叶变换，可以避免插值法由 于插值方式不同而对结果产生主观性影响，并且相对于天文学上的DCF （Discrete Correlation Function） 方法（R. A. Edelson、J. H. Crolik，1988）有着一系列的优点，但同时ZDCF也有一定的缺点，包括对于 分辨率尺度内有自相关的序列产生较大误差和对某些数据给出的CCF的时间延迟范围较小，但多数数据符合 ZDCF的应用条件，并可以基于ZDCF得到可靠的结果。第四纪和第三纪晚期到大约9OO ka冰期变化周期由40ka 转变为100ka的中更新世转变(MPT，mid-Pleistocene transition）问题一直是古气候研究中的一个有较多 争议的问题。低纬度地区的碳同位素和高纬度氧同位素的ZDCF表明，气候转型前气候系统的内部驱动基本 不构成气候周期的驱动因素，而转型后呈现一定的相关性。太阳辐射和高纬度氧同位素的ZDCF显示，气候 转型之后相关性加强，太阳辐射驱动的影响加强。从而支持了Scott Rutherford等（2000）关于0.9Ma气候 转型机制的理论，即造成气候转型的原因是高低纬度之间的气候联系或南北半球之间的热量传输，同时太 阳辐射也是主要的驱动机制。 关键词：ZDCF,，MPT，相关性，太阳辐射，轨道，低纬度，热量传输 1 引言 CCF 已在大气科学的各个研究领域取得了广泛应用。计算不同时间序列的 CCF 有助于 发现物理因子的相关性，从而进行深入机制的研究，也可以检验理论与数据的相符程度。对 于对时间连续分布的时间序列可以直接计算 CCF 。但对于多数时离散形式的古气候数据， 传统的做法是先用内插值的方法使时间序列对时间连续分布，再代入计算 CCF 。但这样得 到的结果通常因插值方法的不同而有所差异，从而带有一定的主观性。EdelsonKrolik （1988）在天文上提出的 DCF 方法为离散时间序列相关性的处理提供了另一种途径，但 DCF 由于统计学上的缺陷，如装箱方法（binning method ）的选取问题、理论本身的有偏性 问题，通常无法得到优于插值法的结果。ZDCF 作为 DCF 方法的改进，不仅克服了 DCF 方 法的诸多缺点，还加入了误差估计，从而更完整有效地反映离散时间序列的相关性。 MPT 的物理成因问题一直是古气候研究领域的一个有争议的问题。基于 ZDCF 方法针 对 MPT 问题对一些古气候数据进行了分析，包括高低纬度之间的辐射梯度驱动理论 （RaymoNisancioglu ，2003 ），结合年平均辐射的理论（Loutre 等，2004 ）和高低纬度或跨 越赤道的能量传输（Scott 等，2000 ），但始终没有定论，并且有些理论缺乏支持的证据。本 文基于 ZDCF 方法对 MPT 可能的太阳辐射和能量传输影响因子进行了研究。 2 ZDCF算法介绍 2.1 基于 Z 变换的相关函数估计 气候学上的要素之间的相互作用可能需要一段时间，因此两已知时间序列的相关性不应 该只考虑零时间延迟点