1资料与方法本
本研究所采用的气象资料是国家气象信息中心整编的月值数据,经过质量控制和检查验证数据无误。
研究将不同时间段的平均气温作为研究对象,采用Mann-Kendall非参数秩次统计检验法和最小二乘法等方法来分析气温变化的总趋势,并采用距平法和UFk统计量等方法分析气温变化的阶段性,采用滑动t检验和M-K非参数检验法来进行气温的突变分析。
在使用M-K检验进行趋势分析时,采用假设H0表示时间序列(x1、x2,…,xn)是数据样本独立同分布,不存在趋势;备择假设H1是双边检验:对于所有的i,j≤n,且i≠j,xi和xj的分布是不同的。检验的统计变量S的计算如下式:
其中,t为任意给定结点的时间范围,n为样本数,Zc为经过修正的Kendall符号检验统计量,收敛于标准正态分布,并通过下式计算:
在双边的趋势检验中,给定置信水平α,如果经过修正的Kendall符号检验统计量Zc≥标准正态分布的(1-α/2)分位数±Z1-α/2,就可以拒绝原假设H0,即在置信水平α下,时间序列数据存在明显的上升或下降趋势。
其中,变化趋势的大小可用Kendall倾斜度β来表示,其计算公式如下:
其中,1<j<i<n。当β>0时,表示上升的趋势,反之则表示下降的趋势。
除用于趋势分析外,M-K方法还可用于突变检验。对于具有n个样本量的时间序列x_1,x_2,…,x_n,可以构建一秩序列: