王雯燕,张颖梅,唐文哲,魏俊涛,李晓冬
(西安市气象局,西安 710016)
王雯燕,张颖梅,唐文哲,等.西安区域自动站气候界限值订正及质量控制 [J].陕西气象,2017(5):38-40.
1006-4354(2017)05-0038-03
2017-03-06
王雯燕(1972—),女,陕西大荔人,汉族,学士,高工,从事综合观测数据分析应用。
西安市气象局科研创新基金“自动气象站数据质量控制方法研究”(2014M-7)
西安区域自动站气候界限值订正及质量控制
王雯燕,张颖梅,唐文哲,魏俊涛,李晓冬
(西安市气象局,西安 710016)
基于天气学原理,探索不同高度区域自动气象站气候界限值的订正方法,在西安区域自动站质量控制软件中试运行,并对运行结果进行分析。结果表明:用相关系数法计算区域站的参考站符合大气变化规律,区域站要素界限值高度订正方法具有代表性,方法简单,易于操作;区域站温度传感器出现问题少,气压传感器故障较多。此方法为其他地方区域站实时质量控制提供参考。
区域自动站;相关系数;质量控制;气候界限值
近几年城市经济快速发展,对短时临近预报、气象预警、决策服务和重大气象服务保障要求越来越高,气象资料的质量直接影响预报、服务效果。北京、上海、广西、湖北、江西等地相继开展了气象资料实时质量控制研究工作[1-5],陕西白水城[6]从数据生产、传输、检查等几个方面设计了自动气象站数据质量控制体系,理念先进,然而在实际建站过程中,选择较好的设备和符合规范的观测场所是较难解决的问题。在设备和测站环境暂时无法改变的情况下,气候背景资料检查是一种好的统计检验方法[7],在气象数据资料质量审核中得到普遍应用。为了提高区域自动气象站(下简称“区域站”)观测的数据质量,西安也基于气候极值开展了相应的质量控制方法。目前,列入中国气象局、陕西省气象局考核的站点122个,其中只有41站建站时间较长,历史记录时长等于或接近10年,还有近半数站点为近3年所建,统计数据不适合作为极值指标使用;加之西安地形多样,有平原、丘陵、沟、塬、山地等,导致区域站海拔高度高差较大,最低352 m,最高1 576 m,复杂的地形特点和较短的历史观测资料为区域站用气候极值法进行实时质量控制增加了难度。本文在传统质量控制理论[8-10]的基础上,根据气象要素空间一致性特点,提出西安地区获取区域站气候界限值的方法,投入业务应用并对运行结果进行分析,探索观测数据误差产生的原因,为区域站实时质量控制提供借鉴。
西安国家基本气象站7个,其中,泾河国家基本气象站2005年由西安国家基准站迁站建成,原西安国家基准站建于1951年,迁站后原址保留温度、气压、湿度等六要素观测,其历史资料连续性较好,城区天气预报预测通常以该站观测资料为基准,因此本研究也将原址观测的数据作为国家基本气象站资料使用,各站建站时间如表1所示。
表1 西安国家基本气象站及建站时间
西安区域站2006年7月开始建设,随后站点不断增加,截止2016年底,投入业务使用并参加中国气象局、陕西省气象局考核的站点有122个。
空间一致性是指气象参数具有一定的空间分布特点,其有效性取决于观测站网的密度和被检参数与空间的相关程度[9]。利用空间一致性确定区域站气候界限值的主要思路是选择与区域站相关系数较高的自动站作为参考站,再根据要素随高度变化规律进行订正。
2.1 区域站气候界限值计算方法
区域站气候界限值计算具体步骤如下。
(1)逐月统计7个国家基本气象站近30年温度、气压等要素的最高和最低气候界限值,建立各站气候界限值序列。
(2)以区域站全月逐日小时数据作为一个序列,逐对计算与基本气象站数据之间的相关系数,采用T分布检验法进行显著性检验[11]。在通过显著性检验的台站中,指定相关系数最大的基本气象站为参考站。
互联网 20世纪60年代,美国和苏联卷入了古巴导弹危机,核战争似乎一触即发,每个国家都在考虑核攻击后的情景。兰德研究员保罗·巴兰(Paul Baran)试图解决核攻击后的一个问题:在核袭击造成毁灭性后果后,如何保持官方沟通渠道的畅通。他的答案为互联网提供了一个框架。
(3)根据区域站与参考站的海拔高度差,计算区域站的气候界限值。
2.2 区域站气候界限值订正公式
2.2.1 温度 温度界限值订正方法是依据气象学基本原理,高度每升高(降低)100 m,温度降低(升高)0.65 ℃,订正计算公式如下:
月最高温度界限值
tmax=tmax1-(h1-h)×0.006 5 ,
(1)
月最低温度界限值
tmin=tmin1-(h1-h)×0.006 5。
(2)
式中,tmax1为参考站月最高温度界限值;tmin1为参考站月最低温度界限值;h1为区域站海拔高度;h为参考站海拔高度。
2.2.2 气压 按照天气学原理,近地面大气中气压随高度的增加按指数规律递减,气压界限值订正方法采用拉普拉斯压高公式。拉普拉斯压高公式是在大气静力学方程的基础上推导出来适用于垂直范围内气压随高度变化的公式,具体公式如下:
月最高气压界限值
pmax=pmax1/ 10((h1-h)/18 400(1+Tmax/273)),
(3)
月最低气压界限值
pmin=pmin1/ 10((h1-h)/18 400(1+Tmin/273))。
(4)
式中,pmax1为参考站最高气压界限值;pmin1为参考站最低气压界限值;h1为区域站海拔高度;h为参考站海拔高度。Tmax、Tmin为区域站月最高、最低温度的平均值。
2015年9月至2016年5月,对西安两要素以上区域站的10 min数据进行实时质量控制,气压、温度显示有问题的数据数量分别为3 795条、308条,分类统计如表2所示。由表2可知,西安区域站温度存在问题少,气压存在问题较多,下面主要针对气压的疑误数据进行分析。
表2 2015年9月—2016年5月西安区域站气压、温度气候界限值10 min实时质量控制结果 条
3.1 超界限值
气压超界限值有两种情况:超最高值和超最低值。
超最高值主要出现在两个时间时段,2015年11月29日—2015年12月30日和2016年1月23日—24日。第一个时段期间疑误数据1 024条、持续时间较长,查看天气形势并无任何过程发生,检查后发现,疑误数据主要发生在V8852、V8853、V8864三个区域站,观测值比附近同海拔高度的区域站数值明显偏大,由此断定这个阶段发生气压超界限值的区域站气压传感器出现故障;2016年1月22—24日出现的气压超界限值,经分析是西安地区出现寒潮降温天气,近地面有较强的冷高压系统影响,属于正常现象。
超最低值有两种情况:一是固定气压,数值为800.4 hPa,由于区域站配置的气压最低值参数为800 hPa,当传感器没有连接时,传输回来的值通常为800.4 hPa;二是出现超低异常连续值也会提醒疑误数据超界限值,2016年2月24日,V8851站气压连续数小时在883.9 hPa和885 hPa之间跳变,经查是该站的气压传感器出现故障导致。
3.2 无数据
无数据指要素栏为/////的情况,出现原因有两个,区域站建设时中心站需定义宏、配置参数,如果将要素字段设定为空时,传输回来的要素栏显示/////;另一方面,当两要素站使用六要素主板,温度、降水正常传输,其它要素因为传感器没有连接,数据传输回来时要素栏显示/////。
3.3 传输异常
传输异常主要指要素栏出现 9999.9 (温度为99.9)的情况,当区域站传感器没有连接或设备出现故障时,传输回来的要素栏显示9999.9(温度为99.9)。
根据空间一致性特点确定西安地区各区域气象站的参考自动站,讨论了不同海拔高度区域站以100 m为间隔订正计算气候界限值的方法,对传输数据业务试运行结果进行统计分析,得到如下结论。
(1)相关系数法计算的区域站参考站比较符合大气变化规律;区域站要素界限值高度订正方法具有代表性,计算方法简单,易于操作。
(2)区域站温度传感器出现问题少,气压传感器发生故障较多。
(3)理论上,风速也可以依据气象学原理利用高度差进行订正,但在与实况资料对比验证时却发现,订正结果与实测值相差较大,而且由于大部分区域站建站时间较短,没有合适的值可供参考,因此没有对风速实施质量控制;相对湿度受温度、水汽含量等影响,变化较为复杂,在不同的地形位置随高度变化不同,用高度订正不适宜,所以在质量控制参数中只以10%为界设置了相对湿度极端最低值,对低于10%的的值提出报警。
[1] 窦以文,屈玉贵,陶士伟,等.北京自动气象站实时数据质量控制应用[J].气象,2008,34(8):77-81.
[2] 孙娟,胡平.上海自动气象站实时质量控制系统研究[J].气象水文海洋仪器,2009,26(4):38-41.
[3] 王超球,黄理,程爱珍,等.无人值守自动气象站实时数据质量控制技术的探讨[J].气象研究与应用,2007,28(z3):64-66.
[4] 王海军,杨志彪,杨代才,等.自动气象站实时资料自动质量控制方法及其应用[J].气象,2007,33(10):102-109.
[5] 李志鹏,张玮,黄少平,等.自动气象站数据实时质量控制业务软件设计与实现[J].气象,2012,38(3):371-376.
[6] 白水成,李社宏,周林.自动气象站数据质量控制体系设计[J].陕西气象, 2016(3):42-46.
[7] 任芝花,熊安元.地面自动站观测资料三级质量控制业务系统的研制[J].气象, 2007,33(1):19-24.
[8] 王新华,罗四维,刘小宁,等.国家级地面自动站A文件质量控制方法及软件开发[J].气象,2006,32(3):107-112.
[9] 刘小宁,任芝花.地面气象资料质量控制方法研究概述[J].气象科技,2005,33(3):199-203.
[10] 熊安元.北欧气象观测资料的质量控制[J].气象科技,2003,31(5):314-320.
[11] 魏凤英.现代气候统计诊断与预测技术[M].2版.气象出版社,2009:18-19.
P412.1
B