天津百年气温序列的均一性检验订正与极端气温变化趋势.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * 天津百年气温序列的均一性检验订正与极端气温变化趋势 郭军 2010.9.29沈阳 常用的均一性检验方法 气候资料不连续的产生有很多原因： 从台站的迁移和仪器的变更到计算时间平均方法的变化等 。和断点不连续相比，趋势不均一通常更难以检测，因为它们可以叠加在气候变化之上 。? 常用的方法：T检验、Potter方法、标准正态检验（SNHT）方法、二相回归（E-P）、序列均一性的多元分析（MASH）等等 PMF检验法 (Penalized Maximal F test ) —惩罚最大F检验 在没有参考序列可用的条件下，采用PMF检测对天津1891～2009年逐月平均最高和最低气温序列的非均一性进行检验。 PMF检验允许待检测序列在其整个时段上存在一个线性趋势（即趋势成分没有漂移），通过迭代算法，对待检测序列的年循环、线性趋势、落后1自相关进行协调估计。 通过应用经验的惩罚函数，使得误报警率和检验能力的非均匀分布问题也大大减少，利用该方法进行突变点检验，可以不使用参考序列。 序列{Xt}的原假设是通过TPR3（Common-trend two-phase regression model based maximal F test, Wang, 2003）评估回归模型来检测的，即 同时备择假设采用了两相回归模型的形式，即 其中，εt是均值为0、方差为未知常数σ2的随机误差，β为线性趋势成分。认为μ1≠μ2并接受备择假设。T=k时的点称作突变点。 Δ=|μ1-μ2|称为平均突变大小。位于时刻k ∈{2,…,n-1}的变点的统计显著性是通过F统计量来评估的： 其中P(k)是建立的经验性的惩罚因子，其中 RHtestV3 输入数据格式 annual series MM=00 and DD=00. 百年气温序列的检验与订正 时段 站址 经纬度 高度 观测次数 台站环境 1891-1920 不详 39o07’N 117o12’E 不详 3次 不详 1921-1931 自由道22号 39o08’N 117o11’E 16 6次 市区 1932-1954 自由道22号 39o08’N 117o11’E 16 24次 市区 1955-2000 气象台路100号 39o06’N 117o10’E 3.3 4次 市区 2001-2009 友谊路 38o59’N 117o43’E 5.6 4次 市区 天津城市气候观测站气温器测沿革 断点日期 检验统计量 95%置信区间 是否显著 1897年3月 270.2484 22.5434~27.4394 是 1907年6月 658.0518 23.2365~28.3145 是 1926年8月 136.5347 22.8391~27.8010 是 1927年11月 82.5084 23.8504~29.1270 是 1993年12月 97.4607 23.9850~29.3040 是 检测出变点的位置和相关统计信息 逐月平均最低气温PMF检验（95%置信度） 年平均最低气温PMF检验（95%置信度） 断点日期 检验统计量 95%置信区间 是否显著 1907年 48.1053 23.0814-38.7180 是 1925年 12.7305 27.9654- 55.0880 是 1993年 35.0688 27.9057- 54.8958 否 订正后的年平均最低气温 断点日期 检验统计量 95%置信区间 是否显著 1952年10月 109.9827 19.4628- 23.8157 是 断点日期 检验统计量 95%置信区间 是否显著 1952 47.7822 13.7343-27.2992 是 订正后的年平均最高气温 订正前后的年平均气温 订正合理性检验 序列的选择（尽量多的考虑，周边的长序列） 北京----原始序列 北京----谢庄 区域----北京、保定、塘沽、大连、青岛算数平均 辽东半岛平均序列---任国玉等 英国气象局的 Hadley 中心和 EAST ANGLIA大学气候研究中心(CRU)合作制作的全球气温数据集- HadCRUT3, 陆地部分为 CRU 研制, 即 CRUTEM3该数据集中中国及邻近地区长序列站点共计186站(含国内 41站） 各参考序列变化情况 华北区域各站气温情况 天津订正后的序列与个参考序列的比值，1951年以后的比值变化较小，表明1951年以后的天津气温序列的变化与北京、华北、辽东半岛是基本一致的。 天津年平均气温变化 天津百年来极端气温变化 使用的天津市1910～2009年逐日最低气温和最高气温观测资料 依照STARDEX（STAtistical and Regional