安徽省气温序列的均一性检验供参习.doc

安徽省气温序列的均一性检验 　 1，2???? 吴必文2???? 惠? 军1 （1.合肥工业大学理学院?? 合肥?? 230009；2.安徽省气候中心） 　??? 要? 利用ｔ检验法和SNHT法对安徽省部分台站1961~2005年月平均气温序列进行均一性检验。当两种方法检验结果相同时，认为该间断点是客观存在的，其后通过相应台站的历史沿革对检验结果进行了分析，并以合肥站为例给出检验和分析的具体过程。检验结果表明，安徽省大多数台站的温度资料质量基本可靠，而台站迁移则是造成气温序列非均一性的主要原因。 关键词 ??均一性检验?? ｔ检验法??? SNHT?? 气温?? 安徽 引 言 随着全球变暖的趋势越来越明显，气候变化已成为一个备受关注的问题。气候变化研究的基础是建立在能够反映气候真实变化的均一性长序列数据。然而由于观测台站的迁移，周围环境的变化，观测方法和观测仪器的变更等常会使序列数据出现非均一性[1]。 近年来国内外很多专家致力于气候资料的非均一性检验的研究[2~8]，Alexandresson[2]发展了广泛应用的SNHT（标准正态检验）法，Easterling[3]等发展了E - P技术，全球历史气候网采用该技术对全球气温和气压进行均一性检验，Potter[4]对19个台站100年降水序列进行均一性检验。国内气候资料的均一性研究起步较晚。刘小宁等[5]，宋超辉等[6]，鞠晓慧等[7]分别针对降水，平均气温，太阳辐射的均一性方法做了总结，利用这些方法对我国部分台站的年降水量，温度序列，太阳总辐射月总量进行了实用性检验，得出了仪器变化和站址的迁移是引起气象资料非均一性的主要原因。翟盘茂[8]采用了Easterling等的均一性检验手段对我国探空和湿度资料进行检验，得出了20世纪60年代我国探空序列存在明显不均一的结论，指出这种不均一性是和探空仪器的变换和辐射的订正相联系。 本文简要介绍了ｔ检验法和SNHT方法以及检验过程，利用两种方法分别对安徽省20个代表台站1961－2005年月平均气温资料进行了均一性检验，并对合肥站作进一步的分析，给出了相关结论。 1? 检验方法简介 在统计学上，所谓均一性序列，就是序列样本{χi}（i=1，2，L，n）取自同一总体因而具有相同的概率分布函数以及相同的数字特征，如数学期望u。而非均一性序列，是指不同时段的序列{χi}取自不同的样本总体，即存在间断点，间断点前后的序列有不同的数学期望u，u+Δ。而检验这些间断点的核心思路基本上都是选取临近台站，构建参考序列，与检测台站形成比值序列Ｚ进行检验。 ????????????????????????????????????????? 。 本文利用ｔ检验法和SNHT法进行检验，当两种检验方法均检验出间断点存在时，认为该点是合理不连续点，是客观存在的。下面简要介绍这两种方法。 1.1? ｔ检验法 ｔ检验法可以检验两组数据的均值是否有显著性差异，其检验统计量为 ??????????????????????（1） 假设检验统计量?????????????? 　　　　　 （2） 若ｗ值大于临界值[9]就认为序列是非均一的。且间断点位置在ｋ处。 ?????? ?????????? S22分别为序列ｚ前ｋ个，后ｎ-ｋ个样本的方差。 1.2?? SNHT 法 ????? ??????????????????????????????????????????????? ???????? 假设检验统计量为：?????????????????????????4） 　[2]时认为序列存在间断点，且间断点位置在ｋ处。其中 ???? 同上。 2?? 检验过程 2.1?? 参照台站的选取 从1961年起温度资料较为完整的站点中，选择距离待测站点较近、地理环境相似、没有明显不均一点的站点 计算他们之间相关系数。相关系数通过显著性检验（α=0.05），本文选取1961～2005年的资料，序列长度为45，即当相关系数大于0.294时可以通过检验，然后将相关系数最大的5个台站作为参考站点。 2.2?? 序列资料的预处理 2.2.1??? 对于待测台站气温序列?Xi?作变换：????????????? ??。 2.2.2??? ?g(yi)? ，由选出来的5个参照站点加权平均获得： ， 其中?? ， 表示第i个参照站点月平均气温序列的均值， Pi表示第i个参照站与待检台站的相关系数， 　　表示第??i?个参照站加权因子。 2.2.3??? 形成待测序列Z ，??进一步对q?作标准化处理： ??????????????（5） 其中为?的算术平均值，??Sq 为该序列的标准差，?Zi? 就是要进行非均一性检验的序列。此时序列的均值为0，方差为 1。??????? 2.3?? 正态性检验 两