资料质量控制和分析方法.docVIP

第二章 资料、质量控制和分析方法 2.1 研究区域和资料 本文研究范围为石家庄市所辖全部区县。区域边界及其所用的17个气象观测站点位置分布见图2.1。区域经度范围为114.13°～115.21°，纬度范围为37.59°～38.45°，海拔高度在35.9m～255.5m之间。各气象站基本信息见表1。 （表1 应增加观测记录年份和迁站时间等信息） 图2.1 石家庄地区17个气象站点分布 表1 石家庄地区17个观测站基本信息 站点 经度 纬度 海拔高度(m) 石家庄市区 114.42° 38.03° 81 新乐 114.68° 38.35° 70.8 辛集 115.21° 37.91° 35.9 赞皇 114.36° 37.65° 137.4 平山 114.2° 38.25° 131 井陉 114.13° 38.03° 255.5 正定 114.56° 38.15° 71 行唐 114.55° 38.45° 96.2 灵寿 114.38° 38.3° 108.9 藁城 114.81° 38.01° 53.5 晋州 115.06° 38.01° 42.4 无极 114.96° 38.2° 45.4 深泽 115.18° 38.18° 38.1 栾城 114.63° 37.88° 52.9 赵县 114.73° 37.76° 38.5 元氏 114.53° 37.75° 66.4 高邑 114.61° 37.59° 48.9 本文后面气温和降水两章所用资料为石家庄市气象局业务科提供的石家庄地区17个观测站的逐日平均气温、最低、最高气温和降水资料。由于各站资料的时间序列长度不同，为了对比分析，因此在温度、降水变化特征中统一采用1972-2008年37年的资料。周期分析中所用资料为石家庄市区站建站以来即1955年至2008年的逐日平均气温、降水资料。 城市气候一章城市化对气温和降水量的影响分析所用资料为石家庄市区站和16个郊区站1972-2008年逐日平均气温和降水资料。其中5.1和5.2小节所用资料为石家庄市区站以及藁城、元氏、平山和新乐四个距离市区站较近的郊区站1962～2009年近48年的逐日平均气温、最低气温和最高气温。四个郊区站分布在市区站东、南、西、北四个方向，海拔高度与市区站接近，因此不必做高度订正，减少了分析误差；且探测环境受城市化影响较小。分析日变化时采用2009年3月至2010年2月加密的逐时气温资料。 所用资料均由石家庄气象局相关业务单位进行了初步的质量控制，订正了由于各种人为因素造成的错误值。 2.2资料均一化订正 本文借鉴李庆祥等[65]检验和订正国家基本、基准站气温数据时所采用的E-P方法，对所用气温资料进行了均一性检验和订正。周雅清等[32]也曾采用这种方法检验和订正华北地区的气温资料。降水资料的时间序列对于台站位置变动等影响没有温度那么敏感，所以本文对降水资料没有进行严格的均一性检验和订正。 首先对所有台站建立气温变化序列dT/dt，计算待检台站dT/dt 序列和它附近台站的dT/dt序列的相关系数，选取5个与之正相关最大的台站作为参考台站。利用参考台站序列和待检序列相关系数的平方作为权重建立参考dT/dt 序列，再反算得到参考序列。然后采用滑动t检验法检测待检序列和参考序列的差异序列的不连续点，并通过台站历史沿革资料确定其合理性。如果不连续点附近存在站址迁移，则认为该不连续点是合理的，并进一步进行订正；否则保留原始记录。对于那些被确认的有合理不连续点的待检序列，计算其与参考序列的差异序列在不连续点前后5年的差值，以此作为补偿值来订正不连续点。为了与以后的资料合并使用，订正是以最近的资料为基础从后向前进行的。由于上述方法要求序列子段必须多于5年，因此，如果不连续点是在开始年后或结束年前的5年内，检验的结果也只能出现在开始后或结束前的第5年。这就需要参考台站的沿革资料进行鉴别，确定其真实不连续点的位置和补偿值。通过订正，使得包括台站迁移、仪器更换、观测方法改变等对资料均一性的影响尽可能减少到最小。 在对各站年平均气温进行均一性检验过程中（图2.2），发现17个站均在1993年出现断点，排除了迁站和仪器更换的原因，因此各站的年平均气温在1993年是一次正常的突变。各站年平均最低气温（图略）与年平均气温类似，断点出现在1987年和1993年，同样排除了迁站和仪器更换的原因。各站年平均最高气温则比较均一，没有出现一致的断点现象。 图2.2 石家庄地区各站年平均气温均一性检验 （黑实线为=0.05的显著性水平） 通过对各站年平均气温、年平均最高、最低气温进行检验，发现有1个站（5.9%）的平均气温资料存在非均一性引起的偏差，有1个站（5.9%）的平均最高气温资料存在非均一性引起的偏