数值预报在航空气象中的应用.ppt

,数值预报在航空气象预报中的应用,一、WAFS的应用,提 要,WAFS系统及其模式简介 WAFS的GRIB数据介绍 WAFS数值预报产品的应用,WAFS定义,WAFS（英文全称：World Area Forecast System）是国际民航组织和世界气象组织专为航空气象设计的一套航空气象服务系统，1982年建立。 目的：为航空气象用户提供图象和数字格式的气象预报产品。,世界区域预报中心-WAFC,WAFS有两个中心，即伦敦世界区域预报中心（英国国家气象局）和华盛顿世界预报中心（国家环境预报中心 NCEP）。 负责模式的运行与航空气象预报产品的制作,WAFS卫星广播系统组成,WAFS的产品通过三个卫星组成的卫星广播系统向外广播,民航的WAFS接收系统,两颗卫星属于华盛顿区域预报中心，覆盖了太平洋区域（ POR）和大西洋区域（AOR)，叫ISCS（International Satellite Communication system）；另一个卫星属于伦敦区域预报中心，覆盖了印度洋区域（IOR），叫SADIS（Satellite distribution system）。 目前中国民航有4套WAFS接收系统,其中两套SADIS，北京和广州各一套；两套ISCS，北京和上海各一套。 1996年此系统安装调试，1997年在北京民航气象中心开始业务运行。 SADIS和ISCS的所有数据已入民航气象数据库。,美国WAFC模式的简介,The NWS Global Aviation Model,水平坐标： 254谱截断，高斯格点是768384，相当于经纬度格点的0.50.5； 垂直区域：垂直区域从地面（sigma=1）到大气顶层（sigma=0），从地面到顶层一共64层，地面和顶层附近增加其分辨率。最低层的大气压约为997.3hpa，最顶层为0.27hpa； 垂直坐标：sigma坐标，Lorenz格点，Arakawa and Mintz (1974)二次守恒差分方案。 垂直分辨率：64层不等距的sigma层，地面气压1000hpa，低于800hpa以下是15层，100hpa以上是24层。 计算机的运行：T254模式24小时预报需要12分钟。,Improvements in Data Assimilation Systems Increased use of automated aircraft reports,英国模式的介绍,Unified Model,全球大气模式一天运行4次（以00Z、06Z、12Z、18Z为初始场），00Z的资料模式在0150Z开始运行，模式最长预报有效时间168小时（6天），48小时预报结果机器运行需18分钟，144小时预报结果机器运行需50分钟。 唯一的全球格点模式。,Unified Model,基本方程：完全的非近似的原始方程，非静力有限差分模式； 水平格点：球形的经纬度坐标。南北极点分别是90S和90N。升级后的全球模式其水平分辨率是经度上0.8333度（432行），纬度上0.5555度（325列），在中纬度地区水平分辨率近似60公里。Arakawa C的变量交错的格点； 垂直格点：垂直方向38层，Charney_Philips的变量交错的格点。,Insert figures hi_res_model_fig1.ps hi_res_model_fig3.ps,WAFS产品种类,T4传真图 BUFR(Binary Universal Form for the Representation of meteorological data）码重要天气预告 业务气象资料（OPMET资料） GRIB 数据（binary grid data） 未来的产品,GRIB资料的应用,飞行文件 风温度预告图 OPMET情报 飞行计划 制作航线计划 优化航路 飞机配载状况 航路重要天气预报 积冰 颠簸 制作航线计划 任意航线剖面 航站天气预报,WAFS是根据航空气象的特点而设计的天气预报系统，产品丰富。 国家气象局和日本发布的数值预报产品着重点是一般的天气形式预报，没有航空气象需要的一些重要天气要素（如对流层顶高度、最大风的高度、积冰等）的预报，而这些产品都是制作飞行文件所必需的资料。 另外国家气象局和日本发布的数值预报产品对航空气象的应用来说预报时效太长，最短的为24小时，而航空气象要求的预报时效一般为9小时。WAFS的产品则不同，它提供6、12、18、24、30、36小时的预报，特别适合于航空气象业务的应用。,WAFS数值预报产品的优点,数值预报可视化系统,系统流程图,气象数据库中定时提取GRIB资料,数据插值计算,数据分类、统一计算处理成为二进制、无格式数据,GRADS基本图形制作,图形产品二次处理,根据不同天气类型图形分类打包,输出显示,常用资料,降水资料,大风资料,积冰指数,对流指数,对流层顶资料,最大风速资料,业务运行资料,美国数值预报产品,地面气压、总降水资料、对流降水资料。 最大风速层所在高度、风向、风速。 对流层顶所在气压、温度。 2米高度上的温度和湿度。 10米高度上的风向和风速。 标准气压层：1000hPa70hPa九个层次的风向、风速、位势高度、温度、相对湿度和垂直上升速度等。 资料的预报时次为：6小时、12小时、18小时、24小时、30小时、36小时6个预报时次。,最大风速层资料：包括最大风速层所在高度、风向、风速。 对流层顶资料：包括对流层顶所在气压、温度。 标准气压层资料：其中包括850hPa70hPa八个层次的风向、风速、温度、相对湿度等。 资料的预报时次为：6小时、12小时、18小时、24小时、30小时、36小时6个预报时次。,英国数值预报产品,GRIB数据诊断分析,为了使数值预报产品在天气预报中发挥更好的作用，在原基本量的基础上计算出诊断量。如风场引起的低层质量辐合、高层辐散；表示湿度情况的温度露点差、水汽通量、水汽通量散度。表示冷空气强度的温度平流、全风速、地面气压场、温度场高度场的6小时变化；表示能量的如相当位温以及位温平流等。,2002年8月22日1200UTC 500hPa高度和涡度预报,2002年8月22日1200UTC 500hPa温度平流预报,2002年8月22日1200UTC 850hPa高度和涡度预报,2002年8月22日1200UTC 1000hPa流线预报,对流指数： SWEAT指数 TT指数 850hPa与500hPa的温差 e平流 K指数 对流降水量,诊断分析产品对流指数介绍,对流指数,强天气威胁指数简称“SWEAT 指标”，记作“I” I = 12 Td850 + 20（T-49）+28 8+ 5+125 （S+0.2），其中S=sin(a500-a850)，式中各项不为负数。 强天气威胁指数（SWEAT）目前在许多国家和地区得到应用。它是根据328次龙卷风资料和日常预报经验得出的一个预报指数。,上图是根据WAFS美国华盛顿27日1200UTC预报资料制作的SWEAT指数18小时(28日0600UTC)预报场。由图可以看出：北京处于SWEAT指数的梯度区，且指数的值在250以上。 据统计，雷暴出现前SWEAT指数多在150以上。同时，统计显示，在雷暴持续时间较长或伴随强对流天气的雷暴中，SWEAT指数多在300以上，可见SWEAT指数越大，出现雷暴可能性越大，强度也可能越大。,SWEAT指数预报图,K指数（气团指标） K=（T850-T500）+TD850-（T-TD）700,温度直减率,中层饱和程度,低层水汽条件,对流指数,K指数能够反映大气的层结稳定情况, K指数越大,层结越不稳定，但它不能明显表示出整个大气的层结不稳定情况。K值大小与可能出现雷暴活动的关系为： K20，无雷暴 20K25，孤立雷暴 25K30，零星雷暴 30K35，分散雷暴 K35，成片雷暴,对流指数,TT指数 TT = T850 + Td850 - 2T500 全总指数（TT）是20世纪70年代由Miller引入的。TT越大越容易发生对流天气。,位温及位温平流,位温e的定义为空气沿干绝热过程变化到P=1000HPA时的温度。 e在假绝热过程中守恒。大气层结稳定度可用e的垂直分布表征，即当 0时，大气层结呈对流不稳定，即大气中e随高度升高而降低时大气为对流不稳定。位温平流同时代表了温度和湿度的平流，直接反映了相应高度层的e变化趋势，使我们可以更方便的判断各层e的垂直分布变化情况从而更准确地判断大气层结的未来趋势。,对流指数,位温平流,e平流的表达式为：,当中低层的e平流为正，中高层的e平流为负时，说明低层有正的能量输入，高层有负的能量输入，易产生对流性天气且进一步维持。当中低层的e平流为正，中高层的e平流也为正，但低层平流量远大于高层平流量时，也要特别注意不稳定能量的潜在威胁。 因此选取了低层850hPa及中层500hPa的e平流来表征大气层结的对流不稳定变化趋势。,图 为2003年5月28日0600UTC 700hPae值，黑点为北京（根据WAFS美国华盛顿据27日1200UTC预报资料制作）。单位:,为了判断上下层能量平流的差别，我们统计500hPa与850hPa的e平流差值。64次雷暴日中，28次为负值，即500hPa相对于850hPa为小的能量平流，占总数的43.8%。但在伴随10mm以上降水量的9次较强雷暴中，有7次为上负下正，可见中低层e平流差对于强雷雨和持续雷雨的预报具有一定指示意义，但对于一般性雷暴，则指示意义不大。此结果也说明了强雷雨需要较持续的不稳定能量的输入。,2002年8月22日1200UTC 500hPa假相当位温平流预报,2002年8月22日1200UTC 850hPa假相当位温平流预报,垂直速度,垂直速度由积分连续方程得到： p=p0+( )(p0-p) 这里垂直速度仅指大尺度运动的垂直速度，不包括对流性上升与下沉运动的速度。正值为下沉运动，负值为上升运动。,图 为2003年5月28日0600UTC 700hPa垂直速度,阴影区域为上升运动区（根据WAFS美国华盛顿据27日1200UTC预报资料制作）,统计63次雷暴日中，700hPa及850hPa均为上升气流的次数为28次，占总数的44.4%，而有一层以上为上升气流的次数为48次，占总数的76.2%，说明雷暴出现时并非中低层均为上升气流，而是在中低层存在上升气流即有可能出现雷暴。 同时，在20次强雷暴中，两层均为上升气流的次数为12次，占总数的60%；850hPa为上升气流的16次，占总数的80%；700hPa为上升气流的14次，占总数的70%。而且中低层上升速度越大越有利于雷暴的产生。,航路积冰预报,积冰指数(RH-50)*2)*(T*(T+14)/-49) 指数为正，代表孤立的潜在积冰区，而最可能出现积冰的区域，积冰指数接近100。,2002年9月13日0600UTC 500hPa积冰预报,2002年9月13日0600UTC对流层顶高度,2002年9月13日0600UTC最大风速及高度预报,1、 颠簸指数 LP指数法： 当L0时预报有颠簸 ，计算出概率P,诊断分析产品介绍,判断： 当70%P50%，预报有轻度颠簸； 当85%P70%，预报有中度颠簸； 当P85%，预报有严重颠簸；,特点： 同时考虑了风场水平切变、垂直切变和温度场水平切变对颠簸的影响,颠簸指数 Ellrod 指数 垂直风切变 VWS= 散度 DIV= 总形变 DEF= 其中切变项DSH= 拉伸项 DST=,诊断分析产品介绍,根据EI判断： 当8EI4时，预报有轻度颠簸； 当12EI8时，预报有中度颠簸； 当EI12以上时，预报有严重颠簸；,特点： 主要考虑风场垂直和水平切变以及形变对颠簸的影响,2004年11月9日雷雨天气的分析,个例1：,天气过程简述,9日北京时0740-1042本场及周围阵雨，0746-0915雷雨，下午至夜间小雨，整个过程降水量9.8mm。,各气象要素的分析,湿度 稳定度指数（sweat index） 700hpa的风场(低空急流） 能量值 降水量,要素选择说明,主要以美国WAFS的2004年11月8日12UTC的数值预报资料来分析。 本报告分析的各要素从预报雷雨有指示的角度出发，没有明显指示的要素没有分析。,地面的相对湿度的分布表明地面湿度大,低层1000HPA的T-Td，能看出来有明显的水汽输送，温度场表现不明显,700hPaT-Td与风场，700hPa存在着明显的低空急流，形成了底层很强的暖湿平流，加强了层结不稳定,Thetae分布：高能舌区与雷雨的对应区,雷雨出现在700hPa高能舌区,雷雨出现在700hPa高能舌区,对这次雷雨过程，SWEAT指数大于300的区域与雷雨区对应较好。,对流降水量与地面雷雨区对应好，而且从降水量看这次过程的预报与观测吻合好。,对流降水量与地面雷雨区对应好，而且从降水量看这次过程的预报与观测吻合好。,700hPa的涡度区,500hPa的涡度区,小结,总的说来，WAFS的数值预报产品对此次雷雨是有很好的指示意义的。 以上只是对此次雷雨过程一些参数的描述，但不一定适合其它的雷雨过程。 对不同天气形势下出现的雷雨，其气象参数及物理量的指示不一定相同。要针对不同的形势总结参数的使用，有的参数在这次过程有用，但在另一次过程中不一定有明显的指示作用。 任何资料的使用都应积累经验，总结其优势。,个例2:2005年5月30日冷涡雷雨,5月30日至31日雷雨,形势介绍 5月30日08时从500hpa至850有冷涡，地面锋面，随着500hpa的冷空气东移，有对流云生成东移，移动过程加强，中午之后管区开始有对流回波，回波一直维持，本场1420UTC至1908UTC雷雨； 500hpa冷涡后部的高压脊向东伸，预示着冷涡后部有冷空气下渗，31日受低层辐合气流和冷涡的影响，0610UTC至0715UTC本场雷雨，市区冰雹；1052UTC至1218UTCB本场雷雨，其中1110UTC-1130UT冰雹。,5月30日至31日雷雨WAFS,500HPA的系统变化,5月30日至31日雷雨WAFS,近地层辐合带的维持,5月30日至31日雷雨WAFS,稳定度指数,5月30日至31日雷雨WAFS,总结,WAFS数值预报产品在航空气象预报及服务中发挥了很大的作用； 进一步做好释用,使WAFS数值预报产品在航空气象预报中起到更大的作用； 以此为初始场制作中尺度数值预报模式，为民航气象预报及服务提供更精细的产品。,二、日本数值预报,一天之中我们接收四个时次初始场的日本传真图，但业务工作使用的主要是00和12UTC两个时次的传真图资料。,日本传真图的种类 ：,500hPa高度、涡度预报,地面降水和风场预报,500hPa温度和 700hPaT-Td度预报,850hPa 风、温和700hPa垂直速度预报,48小时地面气压、降水预报,48小时500hPa高度、涡度预报,日本数值预报图上的主要内容所用的符号、单位、等值线的间隔、形状