梁凤霞等

摘要 利用吉林市旅游局提供的雾凇实况资料，人工观测站吉林和北大壶的相对湿度资料以及中国气象局下发的T639数值预报产品中的物理量资料，即1 000 hPa与925 hPa的相对湿度、2 m温度、1 000 hPa与850 hPa的温度和10 m风速，对吉林市雾凇预报进行研究，得到T639数值预报产品中各物理量对雾凇预报的预报指标。

关键词 雾凇；T639数值预报；吉林吉林

中图分类号 P456 文献标识码 A 文章编号 1007-5739（2015）14-0243-02

1 资料使用及雾凇标准

1.1 雾凇实况资料

由吉林市旅游局提供的2009年1月至2013年2月雾凇实况资料，该资料由旅游局相关人员在每年冬季雾凇多发期实地观测得到。观测资料中把雾凇分为“好”“一般”“无”3个等级，该等级是由旅游局制定的。2009年1月至2013年2月，旅游局共提供309 d的雾凇实况数据，详细分级情况如图1所示。

1.2 气象要素实况资料

相对湿度资料用人工观测站吉林（站号：54172）和北大壶（站号：54175）资料，时间段是从2009年1月至2013年2月。

1.3 T639数值预报产品

使用由中国气象局下发的T639数值预报产品中的以下4个要素：1 000 hPa与925 hPa的相对湿度、2 m温度、1 000 hPa与850 hPa的温度和10 m风速。其中，T639数值预报资料的时间是2009年1月至2013年2月。

2 数值预报物理量要素的选取及数值预报数据的读取

2.1 数值预报中物理量的选取

由于雾凇是发生在近地面层的现象，因此根据雾凇的形成机理和吉林市气象台多年积累的预报经验以及参考大量文献，最终确定4个近地面层的要素作为影响雾凇是否出现的因素，这4个条件分别为水汽条件、温度条件、风速条件、垂直方向温度差条件。对应数值预报产品中的物理量，水汽条件用1 000 hPa（吉林）和925 hPa（北大壶）相对湿度表示，气温用2 m温度代替，风速用10 m风速代替，垂直方向温度差用1 000 hPa与850 hPa温度差表示。

2.2 数值预报产品中各物理量数据的读取方法

T639数值预报产品中，1 000 hPa与925 hPa的相对湿度、2 m温度和1 000 hPa与850 hPa的温度都是micaps第4类数据，因此要将格点数据转化为站点数据，采用双线性插值法读取数据，而10 m风速是micaps第3类数据，采用反距离平方差值法读取数据。

3 统计数值预报产品中各物理量要素阈值时需要说明的问题

3.1 统计阈值

需要统计阈值的物理量包括4个，分别是1 000 hPa与925 hPa相对湿度、2 m温度、T1000-850和10 m风速。统计物理量阈值时，对某一物理量来说，按吉林（54172）和北大壶（54175）2站分别统计；对某一站点来说，分24 h和48 h逐一统计。阈值分“有”和“无”之间的阈值，在有雾凇出现时，分“好”和“一般”之间的阈值。

3.2 不同站点相对湿度预报

吉林（54172）位于北纬43°46′48″，东经126°31′12″，台站海拔高度198.8 m，北大壶（54175）位于北纬43°15′00″，东经126°22′12″，海拔高度537.0 m。由于2站点海拔高度不同，因此2个站点采用不同等压面上的相对湿度，吉林站用1 000 hPa等压面上的相对湿度预报值，北大壶站点用925 hPa等压面上的相对湿度预报值。

3.3 相对湿度预报值订正

吉林市地理环境比较特殊，其位于长白山向松辽平原的过渡地带，城区坐落于开阔平坦的松花江河谷盆地，周围被峰峦起伏连绵不断的低山丘陵所环绕，松花江穿城而过，江水冬季不封冻，从吉林市区溯松花江而上15 km是丰满水电站，冬季江水通过水轮机组出口流出的水仍有4 ℃，松花江从丰满电站大坝到吉林市区的下游出现了70 km的不封冻的松花江江面，不封冻的江面源源不断地蒸发水汽，在这局部地区形成了独特的小气候，是形成吉林雾凇的关键条件[1]。并且根据平时的工作经验和简单的对比发现：T639产品中的相对湿度与实况相比差距较大。因此，在分析数值预报中水汽条件即相对湿度时，要根据相对湿度实况对相对湿度预报值进行订正。

订正相对湿度预报值时，采用最小二乘法。

设经验公式为Y=AX+B，X为数值预报产品的预报值，Y为相对应的实况值。X取X1，X2，X3……Xn时，Y等于Y1，Y2，Y3……Yn，n为样本数。

利用最小二乘法得到：M=[Yi-（AXi+B）]2

求A和B的值，使M值最小，即将问题转化为求函数M=M（A，B）在哪些点取得最小值问题。

X1，X2，X3……Xn以及Y1，Y2，Y3……Yn是已知的，样本数n也已知，所以A和B可求[2-3]。

4 统计T639产品中各物理量阈值

4.1 相对湿度阈值分析

（1）考虑到雾凇是在夜间形成，而且都是天亮以前就已经形成，并且旅游局每天都是6：30左右给我们提供雾凇实况，因此对于各因子，我们均选用数值预报物理量产品20：00—5：00这一段的预报场资料。对于相对湿度，我们把夜间的相对湿度平均值作为预报因子，先求出20：00—5：00的预报场的平均值，把这个平均值作为数值预报的预报值。下文提到的相对湿度预报值就是指这个相对湿度的平均值。

（2）考虑到相对湿度这个物理量的月变化，采取分月份订正的方法。于是我们将相对湿度预报值与实况值按月分类，分别订正。吉林和北大壶2站的相对湿度经验公式分别如表1、2所示。endprint