大庆市近50年来的气温变化及突变分析

大庆市近50年来的气温变化及突变分析
高洪程1 ，周明丽，娄德君2
（1. 大庆市气象局，163000，2. 齐齐哈尔市气象局，161006）
摘 要： 利用Morlet 小波变换方法和3种突变检测方法对大庆市近50年气温变化的周期特
征及突变特征进行了综合分析，结果表明：近50年来大庆市的气温变化具有1.5a 、4a 左右
的年际尺度的周期变化和12～15a 左右的年代际尺度的周期变化，气温变化的阶段性特征明
显。

气温的冷暖交替具有突变特征，突变点为1987年。

关键词： 气温变化 Morlet 小波变换 突变
1引言
气候变化是一个非常复杂的问题，在存在线性变化的同时也存在着一定的非线性变化。

线性变化包括周期变化和趋势变化等，它是气候变化的重要组成部分，非线性变化指的是气候的相对稳定性和突变现象等。

近50年来大庆市气温总体上呈上升趋势，近50年约上升了
2.8℃，与全球和全国的气温变化趋势一致[1]，但其变化又有明显的局地特征。

深入研究大庆市的气温变化规律，得出其变化的周期性规律和突变特征，对工农业生产、能源利用和保护环境等都具有十分重要的意义。

2方法和资料
本文中周期分析采用小波变换方法，研究表明，将Morlet 小波应用于气候变化分析，小波系数的变化趋势与气候信号的起伏是基本一致的[2]，适用于对气候资料序列进行时频局部化分析。

这里将本文应用的Morlet 小波及其计算公式作一说明。

Morlet 小波母函数形式为
()e e x x icx 2412--=ψπ
其子小波为 )(1)(,a b x a
x b a -ψ=ψ 其中，a 为尺度参数，b 为平移参数。

伸缩尺度a 与Fourier 分析中周期T 的关系为
224c c a
T ++=π
小波变换所用资料为大庆市1957~2004年月平均气温距平序列，资料总样本数n=576。

在进行小波变换之前，为消除边界效应，对资料序列向两端做了对称性延拓[2]。

由于在气候序列的突变检测中累积距平法（CA ）、滑动的t 检验法（MTT ）和Mann-Kendall 法（M-K ）各有利弊，本文结合3种方法对大庆市48年的年平均气温序列可能的突变点进行综合检测。

CA 法、MTT 法和M-K 法的应用已有过很多介绍[3]，本文不再详述。

气温突变所用资料为1957~2004年的年平均气温序列。

3结果分析
3.1 大庆市气温变化的基本特征
图1为1957年至2004年大庆市年平均气温距平随时间的变化。

由图可见，近50年来大庆市气温总体上呈上升趋势，近50年约上升了2.8℃，比全球和全国的气温上升幅度要大
得多。

统计资料得出，10年气温距平平均值（1961～2000年）分别为-1.1、–0.6、0.1、0.7℃,由10年气温距平平均值变化可以看出，每10年气温距平平均值差值在0.6～0.7℃左右，即大庆市的年代气温呈稳定上升趋势。

由图还可以看出，1970年代以前为大庆的较冷时期，温度均为负距平，而1988年以后为较暖时期，气温均为正距平，这与全国的气温增暖时间基本一致。

图1 1957～2004年大庆市年平均气温距平曲线
3.2 大庆市气温变化的多时间尺度分析
图2 1957.1~2004.12大庆市月平均气温距平Morlet小波变换系数的实部图2为Morlet小波系数的实部，由图可以清楚地看出大庆市的气温变化存在多时间尺度结构。

首先，大庆市的气温变化存在着明显的12～15a左右的年代际尺度的周期性变化。

12～15a左右的周期存在于2000年以前的时域中，气温变化表现为三个偏冷期和三个偏暖期，即1960年以前的偏冷期；60年代初至60年代中期的偏暖期，1962年前后最暖；60年代后期至70年代中期的偏冷期，其中1970年前后最冷；以及70年代中期至80年代初期的相对偏暖期、80年代中期至末期的相对偏冷期和90年代的相对偏暖期。

将图2中10a尺度以下的部分放大（图略），可以看出近50年来大庆市的气温变化还存在1.5a、4a、8a左右的年际尺度的周期性变化。

不同尺度的周期信号都比较显著，而且在时域上呈现出阶段性特征，尺度越小对应的冷暖交替越频繁。

4a左右的周期出现在1957~1986年期间以及1995年后，在1957年至70年代末期表现得尤为显著，在1995年后尺度有减小的趋势。

而8a左右的周期在1975年以后至今表现得较为突出。

1.5a左右的周期在时域上表现得比较零乱复杂，这里不作详细分析。

对小波功率进行分析，可以看出不同尺度周期信号的强度，由小波功率谱（图3）可以看出，4a左右的周期振荡最为显著，强度最大，其次为1.5a和12～15a左右的周期较显著，8a左右的周期振荡信号也较强，表现为小波功率的闭合中心，20年以上尺度的周期信号最弱。

经检验，1.5a、4a、12～15a的周期振荡都通过了信度为0.05的F检验，说明这几个尺度的周期振荡对近50年来大庆市的气温变化贡献较大，在气温变化中起主要作用。

各不同
尺度相互作用的综合结果是使得大庆市的气温在70年代以前处于偏冷期，在70年代至80年代中期气温在正常值附近波动变化，而80年代末期至今，气温处于偏暖期。

图3 1957.1~2004.12大庆市月平均气温距平Morlet 小波变换的功率谱
3.3 大庆市气温变化的突变分析
气候突变是普遍存在于气候系统中的一个重要现象，它表现为气候在时空上从一个统计特性到另一个统计特性的急剧变化。

由于均值突变能够较好地反映气候基本状况的变化，本文着重分析气温的均值突变。

-30-25
-20
-15
-10
-5
图4 1957～2004年大庆市年平均气温累积距平曲线
由累积距平曲线（图4）可以看出，1987年以前，大庆市年平均气温呈波动式变化趋势，87年后气温呈直线上升趋势，1987年为突变点。

MTT 法中取n =5年，选取信度01.0=α，查表得αt 为3.355，如αt T j ≥0，则认为在第j 年前、后n 年两时段可能会有突变。

图5为n =5时MTT 法计算出的0T 随时间的变化曲线，可以看出，有两个点通过信度检验，1987年为突变点，即1987年前后的5年大庆市的年平均气温可能存在突变。

由M-K 法计算出的)(k d u 随时间的变化曲线（图略）可以清楚地看到近50年来大庆市的气温具有显著的增暖趋势，曲线1c 和2c 的交叉点位于信度线之外的1982年，结合CA 法和MTT 法结果可以看出该点不是突变点。

综上分析得出近50年大庆市气温的冷暖交替具有突变特征，突变点为1987年。

T的变化曲线
图5 MTT法中
4 小结
通过上述分析，可以得出近50年来大庆市气温变化的主要特征：
近50年来大庆市的气温变化具有多层次的时间尺度结构，时域中存在着1.5a、4a左右的年际尺度的周期变化和12～15a左右的年代际尺度的周期变化，气温变化的阶段性特征明显。

近50年来大庆市气温的冷暖交替具有突变特征，突变点为1987年。

参考文献：
1王绍武. 近百年气候变化与变率的诊断研究[J]. 气象学报，1994，52（3）：261～273.
2林振山，邓自旺. 子波气候诊断技术的研究[M]. 北京：气象出版社，1999.
3符淙斌，王强. 气候突变的定义和检测方法[J]. 大气科学，1992，16（4）：482 ～493.
Period and Jump Analysis of Temperature Variation in Daqing in Recengt 50 years
GAO Hong_cheng1，LOU De_jun2，LI Zhi_min2
(1. Meteorological Bureau of Daqing City，Daqing 163000，2. Meteorological Bureau of Qiqihar City，Qiqihar 161006)
Abstract：The periods and jump characters of temperature variation in Daqing in Recent 50 years are analyzed by using Morlet wavelet transforming and three kinds of jump test methods. The results show that there are 1.5a and 4a annual periodic change and 12～15a decadal periodic change of Daqing’s temperature variation in recent 50 years, and there are jump characters in Daqing’s temperature variation, the jump point lies in 1987.
Keywords: temperature variation Morlet wavelet jump。