长江中下游5-8月降水量的变化特征分析

长江中下游5－8月降水量的变化特征分析

翟景秋 孔玉寿 李艳兵

解放军理工大学气象学院研究生三队（211101）

E-mail:lansedetiankong_916@

摘 要：本文利用线性倾向估计和Morlet 连续小波变换分析了长江中下游地区四个台站1952－1999年5－8月降水量的时间变化特征．结果表明：（1）从1952－1999年近50年长江中下游5－8月降水量呈减少趋势；（2）降水序列存在19～20年和9～14年的显著周期。在周期分析的基础上本文还进行了降水序列的拟合试验，可作为未来降水预测的一个依据。

关键词： 降水量， 特征分析，小波变换，线性倾向估计

1. 引 言

长江中下游地区是我国经济比较发达的地区，旱涝灾害所造成的损失也就更为严重，如1998年长江流域的洪涝灾害给我国造成了巨大的经济损失，对工农业生产和人民的日常生活产生严重的影响。而影响该地区旱涝的降水量主要集中在5－8月份，因此，着重对5－8月降水量的研究具有很重要的现实意义。

本文的目的是用统计诊断技术对长江中下游地区5－8月降水量的时间变化特征进行研究分析：包括趋势分析、周期分析和拟合。以便为气候预测提供科学依据和技术支持。

2. 资料和方法

本文使用1952－1999年长江中下游地区沿江4个测站（武汉、九江、南京、上海）5－8月降水量资料。用线性倾向估计分析各个测站5－8月降水量的变化趋势；利用小波分析来研究长江中下游地区5－8月降水量变化的多时间尺度演变特征，并在此基础上进行拟合试验。

2.1 线性倾向估计方法简介[1]

线性倾向估计是用在气候变化趋势分析中的一种分离气候趋势的统计方法，它的主要思想就是用一条合理的直线表示变量与时间之间的关系。

用i x 表示样本容量为n 的某一气候变量，用表示i t i x 所对应的时间，建立i x 与之间的一元线性回归方程：

i t ˆ(1,2,,i i )x

a bt i n =+= 此式可以看作一种特殊的、最简单的线性回归形式，方程中的回归系数和回归常数可以用最小二乘进行估计。

2.2 小波分析原理[2,3,4]

小波分析最初由法国地球物理学家Morlet 分析地震时间序列时引入的，后来Daubechies - 1 -

与Mallat 等在小波分析的理论与应用方面做了大量工作，使小波分析得以成功应用到各个领域。其优于传统的谐波分析之处在于具有多分辨率（multi-resolution ）也叫多尺度（multi-scale ）的特点，可以由粗及细地逐步观察信号。

本文取复数形式的Morlet 小波作母小波函数，Morlet 小波表达式为

()20/2i t t t e

e ωϕ−= （1） Morlet 小波的时间尺度与周期T 有如下关系：

a T a = （2）

当02ωπ=时T ，尺度a 和周期T 在数量上视为同等。因此，Morlet 小波可用于周期分析

a ≈[4]。

由于Morlet 小波的实部和虚部存在一个2π

的位相差，因此，当观察小波系数的模时可

以消除小波本身的振荡，这样作为结果的复值小波系数的模平方（能量密度）的取值可较好地用于鉴别各种特征尺度及其在轴上的位置。由此，定义小波方差：

t ()()2var ,f a W a b ∞−∞=∫db （3）

在一定尺度下，表示时间序列中该尺度周期波动的强弱（能量大小）．小波方差随尺度变化过程绘制的图，称为小波方差变化图，该图能反映时间序列中所包含的各种尺度（周期）的波动及其强弱（能量大小）随尺度变化的特性．因此，通过此图可以非常方便地查找一个时间序列中起主要作用的尺度（周期）。

()var a 2.3 拟合原理

曲线拟合问题是指：已知一组（二维）数据，即平面上 n 个点(),1i i x y i n =…，寻求

一个函数（曲线），

使它在某种准则下与所有数据点最为接近。拟合有各种方法，最基本的做法是提供一组拟合所用函数来作为输入参数，经计算可以得到该组函数的基本表达形式，这就是拟合结果。本文资料为降水序列，即时间与降水量数据对，如果能拟合得到函数()y f x =()'f t ，则可以直观反映降水量随时间的变化规律，并可作为预测未来降水的一个依据。

3. 趋势分析

利用线性倾向估计研究长江中下游地区5－8月降水量的变化趋势，从表1可以看出：这几个测站的回归系数b 和相关系数r 都是小于零的，表明长江中下游地区5－8月降水量是随着时间的增加而呈下降趋势的。对相关系数进行显著性检验，取显著性水平a ＝0.01，可得=0.3613。从表1可以看出，九江和上海两站的相关系数0.01r r 都是大于，这表明九江和上海两站5－8月降水量的这种下降趋势在a ＝0.01显著性水平上是显著的，而武汉和南京两站均未达到一定的显著性水平，虽然也呈下降趋势，但趋势不明显。由此可以得出0.01r - 2 -

这样的结论：从1952－1999年近50年，无论从长江中下游整体还是从各个测站的角度分析，5－8月降水量均呈下降趋势。

表1 线性倾向估计

回归常数a 回归系数b 相关系数r

南京789.513 -0.10542 -0.144929

上海1040.49 -0.241258 -0.388549

武汉877.324 -0.095842 -0.100879

九江1356.57 -0.333705 -0.393702

4. 周期分析

图1 南京降水序列Morlet复小波方差图

图2 南京降水序列Morlet复小波变换实部图(灰度表示值的大小，白色最大，黑色最小)

以南京降水为例。由图１可以看出，小波方差在周期为2、9和19年时出现极值，说明降水序列存在相应的周期。再从图2来观察水平方向小波系数大小的周期性起伏，可以发现，19年和9年周期在整个时间跨度上都比较明显，2年周期主要出现在1975年之前。以19年时间尺度的降水变化为例，主要经历了5个时期的交替变换：1959年以前的偏多期，1960

- 3 -