在分析河南省冬春季沙尘天气气候特征的基础上,筛选出T213数值预报产品中与风速和能见度相关性较好的预报因子,采用卡尔曼滤波方法,分站建立了河南省冬、春季的风速和能见度短期预报方程并利用相应的沙尘天气分级量化标准,在一定程度上实现了河南省冬、春季沙尘天气的短期分站、分级预报。业务运行时,又增加了两条判别指标进行后处理有效地消除了液滴对能见度预报的影响,提高了沙尘天气预报的准确率。
Based on the analysis of climatic characteristics of dust weather in winter spring in Henan some variables, which have close relationship with visibility and wind speed are selected from the T213 Numerical Product as predicting factors (predictor), and then, with Kalman filtering method, the prediction equation of visibility and wind speed in spring and winter in Henan is developed. In terms of discriminative standard of dust weather the quantitative dust weather forecast in spring and winter for the stations in Henan is realized. In the forecast operation another two discriminative standard were added to eliminate the effect of liquid droplets on the visibility therefore the prediction accuracy is increased.
沙尘是我国西北和华北地区冬、春季经常出现的灾害性天气,它能加剧土地沙漠化,对大气环境造成严重污染,给交通和供电线路产生重要影响。因此,对沙尘天气做出较为准确的预报对防灾减灾具有十分重要的实际意义。我国气象工作者针对沙尘暴已做了很多研究工作,如:对当地沙尘天气进行天气气候学分析和诊断[
卡尔曼滤波方法在气象领域中的应用己经得到认可,尤其是在对连续性变量(如温度、湿度等)的预报方面己得到广泛的应用[
《地面气象观测规范》中将沙尘天气分为浮尘扬沙、沙尘暴三类:①浮尘——尘土、细沙均匀地浮游在空中,使水平能见度<10.0km; ②扬沙——由于大风将地面尘沙吹起,使空气相当混浊,水平能见度在1.0km至10.0km以内。③沙尘暴——由于强风将地面大量尘沙吹起,使空气很混浊,水平能见度<1.0km。这种分级法,虽然在一定程度上反映了沙尘现象的强弱,但没有给出确切的量值,给沙尘天气的深入研究带来了诸多不便。因此,我们参照已有的研究成果[
沙尘天气的分级量化标准
|
利用郑州、安阳等11个国家基准(基本)站1961~2000年40年的气象资料,对河南沙尘暴、扬沙、浮尘天气发生的气候特点、变化趋势等进行了较详尽的分析。
沙尘暴主要发生在我省黄淮平原,尤以沿黄一带及黄河故道、黄泛区最为严重,发生最多的是开封、郑州,平均每年5.4次和3.2次; 相对于东部平原来说,西部丘陵山区及豫南沙尘暴发生较少,一般在10次以下,如卢氏近40年来没有发生过沙尘暴天气。
扬沙天气在各地均有发生,切相差悬殊。发生最多的是豫西的卢氏、三门峡及豫东的西华,平均每年30次、28次和27.5次。最少的是孟津平均每年2.4次。
浮尘天气也是一种全省性的天气现象,其区域分布不同于扬沙,而与沙尘暴相似,主要分布在沿黄一带及豫东黄泛区,全省以郑州最多,平均每年28.3次,较少的是固始、栾川、西峡,平均每年不到1次。
不同年代和年份,河南各地沙尘暴、扬沙、浮尘的发生有显著的不同,表现为随着时间的推移,其发生次数有明显的起伏变化,尤其从20世纪70年代后,呈现出显著的递减趋势。
沙尘暴:沙尘暴的发生具有明显的年代际变化。以郑州、开封为代表的全省大多数地区,40年中以60年代最多,之后逐渐减少,至90年代最少,其中60年代的沙尘暴总数占40年总数的50%~75%。
扬沙:河南扬沙的发生次数和发生频率要远远高于沙尘暴。全省绝大多数地区以60~80年代较多,90年代扬沙天气最少。
浮尘:河南浮尘的发生次数和发生频率也远远高于沙尘暴。其年代际变化与扬沙相似,全省绝大多数地区以60~80年代较多,90年代浮尘天气最少。
河南沙尘暴、扬沙和浮尘的发生有明显的季节变化,冬、春季较多,夏、秋季较少。
沙尘暴主要发生在12月~次年6月份的7个月里,一般占全年总数的85%~100%, 又集中在春季的3~5月,占50%~100%。
扬沙的季节变化与沙尘暴相似,但主要发生在11月~次年5月份的7个月里,一般占全年总数的88%~93%, 又集中在春季的3~5月,一般占40%~60%。
浮尘的季节变化与沙尘暴、扬沙相似,但主要发生在11月~次年6月份的8个月里,一般占全年总数的90%以上,又集中在冬、春的1~4月,一般占40%~60%。
卡尔曼滤波方法是目前国际上较为先进的一种新的递推滤波方法。卡尔曼滤波方法通过利用前一时刻预报误差反馈到原来的预报方程,并及时修正预报方程系数,以此提高下一时刻的预报精度。
若滤波对象为离散时间线性动态系统,并认为天气预报对象是具有这种特征的的动态系统,可用下列两组方程来描述:
式(1)为预报方程,
式(2)为状态方程。可以简化为
式(1)和式(3)用来描述离散时间的线性动态系统,用最小二乘法,可推导出一组递推滤波系统,在实际运算中,通常是利用递推滤波系统来进行预报量的计算。
由于卡尔曼滤波对象是离散时间的线性动态系统,因此,将沙尘这个不连续量的预报转化为对风綱随这两个连续量的预报,通过对它们的预报达到对沙尘天气的预报。
本方法采用T213的每6h一次的预报场(昨日20点发布的30、36、42、48小时预报场)作为预报资料,分站建立每6h一次的风速和能见度的预报方程。建立方程时,使用了2003年11~12月各代表站的02、08、14、20时四个时次的风速和能见度实况资料和T213的预报资料。虽然用小样本建立的预报方程统计特性差,但卡尔曼滤波方法自身可以不断通过误差反馈订正预报方程系数来适应预报对象的统计特征[
为了在实际预报工作中能够实现预报的自动化我们选取了能够直接获取风速和能见度资料的安阳、新乡、三门峡等11个站进行站点沙尘天气预报。
针对每一站点,将风速和能见度作为预报量,选取预报因子,采用卡尔曼滤波方法建立预报方程,并做逐日滚动预报。
由于卡尔曼滤波初始方程是建立在少量历史资料基础上的,为了减少干扰,预报因子一般不超过4个[
风速和能见度往往与特定的地面及高空天气形势相联系。通过普查历史天气图得出,导致沙尘天气发生的天气系统一般为:高空有一较深的西风槽,槽后有强冷中心; 地面有一庞大的冷高压,高压前有气旋强烈发展。另外,天气形势的变化与动力场和热力场的变化也有较大关联。
因此,在进行风速和能见度预报因子的选取时,首先根据上述形势分析和物理量分析结论选取了一些初级因子,并增加一些复合因子,然后经过反复分析和相关性计算确定最终入选因子。不同区域的站点,预报因子并不完全相同。
下面以郑州站为例,给出风速和能见度的预报因子,因子中用到的各站点物理量均通过T213产品每6h一次的格点预报资料进行插值到站点计算后获得。
郑州站风速的预报选取了地面气压差、地面变压差、100hPa全风速及850hPa温度梯度四个因子。
(1) 地面气压差
(2) 地面24小时变压差
(3) 1000hPa全风速
(4) 温度梯度
郑州站能见度的预报选取了500hPa高度差、稳定度、湿度项、垂直速度四个因子。
(1) 500hPa高度差
(2) 近地面层稳定度
(3) 近地面层湿度项
(4) 上升运动项
利用T213数值预报产品的6h物理量预报场,计算出所需的预报因子,并读取当日各站的风速和能见度实况,用卡尔曼滤波方法对预报方程系数进行修正后,即可实现风速和能见度的滚动预报。
利用
利用2004年1~5月的资料,对建立的沙尘预报方法进行了业务试验。从11个站的分站预报结果(
用卡尔曼滤波方法做全省沙尘天气预报的试验结果(2004年1~5月)
|
分析空报原因主要是因为沙尘天气产生因子比较复杂,除了与地形有关外,与本站的前期气温、降水、地表状况和干旱情况等气候背景特征也有很大关系。另外,对于能见度的预报方程而言,因子无法排除液滴作用(降雨或大雾)所带来的能见度下降,虽然有时预报能见度较低,但该站当时有降水天气出现时,沙尘天气一般不会发生。
根据上述分析,将预报方法进行了后处理,増加了两条判别指标。
(1) 方程输出沙尘预报等级后,再利用该站前15天的降水实况与历史平均值进行对比,如果总降水量比历史同期偏多20%以上时,该站沙尘天气预报结果降低一个等级(如方程预报结果为沙尘暴,则实际发布预报结论为扬沙,以此类推)。
(2) 如果预报站点前一天实况己经出现降水或T213数值产品预报6h内有降水发生时,该站沙尘天气预报结果相应降低一个等级。
建立的河南省冬春季沙尘天气短期预报业务流程见
河南省冬春季沙尘天气短期预报业务流程图
2004年11月-2005年5月,河南省冬分析空报原因主要是因为沙尘天气的春季沙尘短期预报方法正式投入了业务运行,并表现出了较高的预报准确率。特别是对2005年河南省最强的两次沙尘天气(3月10日和4月7日),均做出了准确的预报。从全省11个站的沙尘业务预报结果(
正式投入业务运行的沙尘预报方法的预报结果(2004年12月~2005年5月)
|
(1) 试验表明,将对沙尘这个不连续量的预报转换为对连续的风速和能见度的预报是可行的,设计的预报方法是成功的。用卡尔曼滤波方法做沙尘天气预报是可行的,但必须利用前期降水量、地表状况等气候特征指标进行后处理,以减少预报方法的空报率。
(2) 此方法在一定程度上实现了沙尘天气的分站、分级预报,但目前只适用于基本站,而对于非基本站,由于无法及时得到6h一次的风速、能见度和降水量等实况资料,因此无法利用卡尔曼滤波方法及时订正方程系数,进行沙尘预报。
本文得到河南省气象局李平高级工程师的指导和帮助在此表示衷心感谢。
彭维耿, 赵光平, 陈豫英.宁夏春季沙尘暴与气象要素及环流指数的关系.气象, 2005, 31(3):17-21.
张仁健, 徐永福, 韩志伟.北京春季沙尘暴的近地面特征.气象, 2005, 31(2):8-11.
申红喜, 李秀连, 石步鸠.北京地区两次沙尘(暴)天气过程对比分析.气象, 2004, 30(2):12-16.
赵翠光.人工神经元网络方法在沙尘暴短期预报中的应用.气象, 2004, 30(4):39-41.
黄嘉佑, 谢庄.卡尔曼滤波在天气预报中的应用.气象, 1993, 19(4):3-7.
陆如华, 徐传玉, 张玲等.卡尔曼滤波在天气预报中的应用技术.数值预报产品评价公报, 1996, 5-6双月刊:28-36.
陈雷, 刘开富, 李英.卡尔曼滤波在短期气候预测中的应用.气象, 2001, 27(10):42-45.
魏文秀, 任彪, 杨海龙等.卡尔曼滤波技术在暴雨中期预报中的应用.气象, 1998, 24(3):46-49.
张冬斌, 尚可政等.沙尘天气的定量化指数及其应用.干旱气象, 2003, 21(3):58-61.