1、数值预报在航空气象数值预报在航空气象预报中的应用预报中的应用数值预报在航空气象预报中的应用一、WA F S 的应用WAFS系统及其模式简介WAFS的GRIB数据介绍WAFS数值预报产品的应用提 要WA F S 系统及其模式简介WAFS（英文全称：World Area Forecast System）是国际民航组织和世界气象组织专为航空气象设计的一套航空气象服务系统，1982年建立。目的：为航空气象用户提供图象和数字格式的气象预报产品。WA F S 定义WA F S（英文全称：Wo r l d A r e a F o rWAFS有两个中心，即伦敦世界区域预报中心（英国国家气象局）和华盛顿世界预报
2、中心（国家环境预报中心 NCEP）。负责模式的运行与航空气象预报产品的制作世界区域预报中心-WA F C WA F S 有两个中心，即伦敦世界区WA F S 卫星广播系统组成WA F S 的产品通过三个卫星组成的卫星两颗卫星属于华盛顿区域预报中心，覆盖了太平洋区域（POR）和大西洋区域（AOR)，叫ISCS（International Satellite Communication system）；另一个卫星属于伦敦区域预报中心，覆盖了印度洋区域（IOR），叫SADIS（Satellite distribution system）。目前中国民航有4套WAFS接收系统,其中两套SADIS，北京和
3、广州各一套；两套ISCS，北京和上海各一套。1996年此系统安装调试，1997年在北京民航气象中心开始业务运行。SADIS和ISCS的所有数据已入民航气象数据库。民航的WA F S 接收系统两颗卫星属于华盛顿区域预报中心，覆盖 美国美国WAFC模式的简介模式的简介水平坐标：254谱截断，高斯格点是768384，相当于经纬度格点的0.50.5；垂直区域：垂直区域从地面（sigma=1）到大气顶层（sigma=0），从地面到顶层一共64层，地面和顶层附近增加其分辨率。最低层的大气压约为997.3hpa，最顶层为0.27hpa；垂直坐标：sigma坐标，Lorenz格点，Arakawa and Mi
4、ntz(1974)二次守恒差分方案。垂直分辨率：64层不等距的sigma层，地面气压1000hpa，低于800hpa以下是15层，100hpa以上是24层。计算机的运行：T254模式24小时预报需要12分钟。T h e N WS G l o b a l A v i a t i o n Mo d e l 水I m p r o v e m e n t s i n D a t a A s s i m i l a t 英国模式的介绍英国模式的介绍 英国模式的介绍全球大气模式一天运行4次（以00Z、06Z、12Z、18Z为初始场），00Z的资料模式在0150Z开始运行，模式最长预报有效时间168小时（6
5、天），48小时预报结果机器运行需18分钟，144小时预报结果机器运行需50分钟。唯一的全球格点模式。U n i f i e d Mo d e l 全球大气模式一天运行4 次（以0 0 Z基本方程：完全的非近似的原始方程，非静力有限差分模式；水平格点：球形的经纬度坐标。南北极点分别是90S和90N。升级后的全球模式其水平分辨率是经度上0.8333度（432行），纬度上0.5555度（325列），在中纬度地区水平分辨率近似60公里。Arakawa C的变量交错的格点；垂直格点：垂直方向38层，Charney_Philips的变量交错的格点。U n i f i e d Mo d e l 基本方程：完
6、全的非近似的原始方程，Insert figureshi_res_model_fig1.pshi_res_model_fig3.psI n s e r t f i g u r e sT4传真图BUFR(Binary Universal Form for the Representation of meteorological data）码重要天气预告业务气象资料（OPMET资料）GRIB GRIB 数据（数据（binary grid databinary grid data）未来的产品 WA F S 产品种类T 4 传真图飞行文件风温度预告图OPMET情报飞行计划制作航线计划优化航路飞机配载状况
7、航路重要天气预报积冰颠簸制作航线计划任意航线剖面航站天气预报G R I B 资料的应用飞行文件WAFS是根据航空气象的特点而设计的天气预报系统，产品丰富。国家气象局和日本发布的数值预报产品着重点是一般的天气形式预报，没有航空气象需要的一些重要天气要素（如对流层顶高度、最大风的高度、积冰等）的预报，而这些产品都是制作飞行文件所必需的资料。另外国家气象局和日本发布的数值预报产品对航空气象的应用来说预报时效太长，最短的为24小时，而航空气象要求的预报时效一般为9小时。WAFS的产品则不同，它提供6、12、18、24、30、36小时的预报，特别适合于航空气象业务的应用。WA F S 是根据航空气象的特
8、点而设计的天气预报系统，产品丰富。数值预报可视化系统气象数据库气象数据库中定时提取中定时提取GRIB资料资料 数据插值计算数据插值计算数据分类、统一计数据分类、统一计算处理成为二进制算处理成为二进制、无格式数据、无格式数据GRADS基本基本图形制作图形制作图 形 产 品图 形 产 品二次处理二次处理根据不同天根据不同天气类型图形气类型图形分类打包分类打包输出显示输出显示常 用常 用资料资料降水资料降水资料大风资料大风资料积冰指数积冰指数对流指数对流指数对流层顶资料对流层顶资料最大风速资料最大风速资料业务运业务运行资料行资料系统流程图 气象数据库中定时提取G R I B 资料 数据插值计算地面气
9、压、总降水资料、对流降水资料。地面气压、总降水资料、对流降水资料。最大风速层所在高度、风向、风速。最大风速层所在高度、风向、风速。对流层顶所在气压、温度。对流层顶所在气压、温度。2米高度上的温度和湿度。米高度上的温度和湿度。10米高度上的风向和风速。米高度上的风向和风速。标准气压层：标准气压层：1000hPa70hPa九个层次的风向、风速、位势九个层次的风向、风速、位势高度、温度、相对湿度和垂直上升速度等。高度、温度、相对湿度和垂直上升速度等。资料的预报时次为：资料的预报时次为：6小时、小时、12小时、小时、18小时、小时、24小时、小时、30小小时、时、36小时小时6个预报时次。个预报时次。
10、美国数值预报产品最大风速层资料：包括最大风速层所在高度、风最大风速层资料：包括最大风速层所在高度、风向、风速。向、风速。对流层顶资料：包括对流层顶所在气压、温度。对流层顶资料：包括对流层顶所在气压、温度。标准气压层资料：其中包括标准气压层资料：其中包括850hPa70hPa八个八个层次的风向、风速、温度、相对湿度等。层次的风向、风速、温度、相对湿度等。资料的预报时次为：资料的预报时次为：6小时、小时、12小时、小时、18小时、小时、24小时、小时、30小时、小时、36小时小时6个预报时次。个预报时次。英国数值预报产品英国数值预报产品 英国数值预报产品 为了使数值预报产品在天气预报中发挥为了使数
11、值预报产品在天气预报中发挥更好的作用，在更好的作用，在原基本量原基本量的基础上计算出的基础上计算出诊断量诊断量。如风场引起的。如风场引起的低层质量辐合、高低层质量辐合、高层辐散层辐散；表示湿度情况的；表示湿度情况的温度露点差、水温度露点差、水汽通量、水汽通量散度汽通量、水汽通量散度。表示冷空气强度。表示冷空气强度的的温度平流、全风速、地面气压场、温度温度平流、全风速、地面气压场、温度场高度场的场高度场的6小时变化小时变化；表示能量的如；表示能量的如相相当位温以及位温平流当位温以及位温平流等。等。G R I B 数据诊断分析 为了使数值预报产品在天气预报中2 0 0 2 年8 月2 2 日1 2
12、 0 0 U T C 5 0 0 h P a 高度和涡度预2 0 0 2 年8 月2 2 日1 2 0 0 U T C 5 0 0 h P a 温度平流预报2 0 0 2 年8 月2 2 日1 2 0 0 U T C 8 5 0 h P a 高度和涡度预2 0 0 2 年8 月2 2 日1 2 0 0 U T C 1 0 0 0 h P a 流线预报对流指数：SWEAT指数TT指数850hPa与500hPa的温差e平流K指数对流降水量对流指数：诊断分析产品对流指数介绍强天气威胁指数简称“SWEAT 指标”，记作“I”I=12 Td850+20（T-49）+28 8+5+125 （S+0.2），
13、其中S=sin(a500-a850)，式中各项不为负数。强天气威胁指数（SWEAT）目前在许多国家和地区得到应用。它是根据328次龙卷风资料和日常预报经验得出的一个预报指数。对流指数 强天气威胁指数简称“S WE A T 指标”，记作“I”上图是根据WAFS美国华盛顿27日1200UTC预报资料制作的SWEAT指数18小时(28日0600UTC)预报场。由图可以看出：北京处于SWEAT指数的梯度区，且指数的值在250以上。据统计，雷暴出现前SWEAT指数多在150以上。同时，统计显示，在雷暴持续时间较长或伴随强对流天气的雷暴中，SWEAT指数多在300以上，可见SWEAT指数越大，出现雷暴可能
14、性越大，强度也可能越大。SWEAT指数预报图指数预报图 上图是根据WA F S 美国华盛顿2 7 日1 2 0K指数（气团指标）（气团指标）K=（T850-T500）+TD850-（T-TD）700 温度直减率中层饱和程度低层水汽条件对流指数对流指数 K 指数（气团指标）温度直减率中层饱和程度低层水汽条件对K指数能够反映大气的层结稳定情况,K指数越大,层结越不稳定，但它不能明显表示出整个大气的层结不稳定情况。K值大小与可能出现雷暴活动的关系为：K20，无雷暴 20K25，孤立雷暴 25K30，零星雷暴 30K35，分散雷暴 K35，成片雷暴，成片雷暴K 指数能够反映大气的层结稳定情况,K 指数
15、越大,层结越不稳定 TT指数指数TT=T850+Td850-2T500全总指数（TT）是20世纪70年代由Miller引入的。TT越大越容易发生对流天气。对流指数 T T 指数 位温e的定义为空气沿干绝热过程变化到P=1000HPA时的温度。e在假绝热过程中守恒。大气层结稳定度可用e的垂直分布表征，即当 0时，大气层结呈对流不稳定，即大气中e随高度升高而降低时大气为对流不稳定。位温平流同时代表了温度和湿度的平流，直接反映了相应高度层的e变化趋势，使我们可以更方便的判断各层e的垂直分布变化情况从而更准确地判断大气层结的未来趋势。Pe对流指数对流指数 位温及位温平流 位温 e 的定义为空气沿干绝热
16、过e平流的表达式为：dydvdxduee 当中低层的e平流为正，中高层的e平流为负时，说明低层有正的能量输入，高层有负的能量输入，易产生对流性天气且进一步维持。当中低层的e平流为正，中高层的e平流也为正，但低层平流量远大于高层平流量时，也要特别注意不稳定能量的潜在威胁。因此选取了低层850hPa及中层500hPa的e平流来表征大气层结的对流不稳定变化趋势。位温平流 e 平流的表达式为：图 为2003年5月28日0600UTC 700hPae值，黑点为北京（根据WAFS美国华盛顿据27日1200UTC预报资料制作）。单位:图 为2 0 0 3 年5 月2 8 日0 6 0 0 U T C 7 0
17、 0 h P a e 值 为了判断上下层能量平流的差别，我们统计500hPa与850hPa的e平流差值。64次雷暴日中，28次为负值，即500hPa相对于850hPa为小的能量平流，占总数的43.8%。但在伴随10mm以上降水量的9次较强雷暴中，有7次为上负下正，可见中低层e平流差对于强雷雨和持续雷雨的预报具有一定指示意义，但对于一般性雷暴，则指示意义不大。此结果也说明了强雷雨需要较持续的不稳定能量的输入。为了判断上下层能量平流的差别，我们统计5 0 0 h2 0 0 2 年8 月2 2 日1 2 0 0 U T C 5 0 0 h P a 假相当位温平2 0 0 2 年8 月2 2 日1 2
18、 0 0 U T C 8 5 0 h P a 假相当位温平垂直速度由积分连续方程得到：p=p0+()(p0-p)这里垂直速度仅指大尺度运动的垂直速度，不包括对流性上升与下沉运动的速度。正值为下沉运动，负值为上升运动。yvxu垂直速度 垂直速度由积分连续方程得到：图 为2003年5月28日0600UTC 700hPa垂直速度,阴影区域为上升运动区（根据WAFS美国华盛顿据27日1200UTC预报资料制作）图 为2 0 0 3 年5 月2 8 日0 6 0 0 U T C 7 0 统计63次雷暴日中，700hPa及850hPa均为上升气流的次数为28次，占总数的44.4%，而有一层以上为上升气流的
19、次数为48次，占总数的76.2%，说明雷暴出现时并非中低层均为上升气流，而是在中低层存在上升气流即有可能出现雷暴。同时，在20次强雷暴中，两层均为上升气流的次数为12次，占总数的60%；850hPa为上升气流的16次，占总数的80%；700hPa为上升气流的14次，占总数的70%。而且中低层上升速度越大越有利于雷暴的产生。统计6 3 次雷暴日中，7 0 0 h P a 及8 5 0 h P a积冰指数积冰指数(RH-50)*2)*(T*(T+14)/-49)指数为正，代表孤立的潜在积冰区，而指数为正，代表孤立的潜在积冰区，而最可能最可能出现积冰的区域，积冰指数接近出现积冰的区域，积冰指数接近1
20、00。航路积冰预报 积冰指数(R H-5 0)*2)*(T*(T+2 0 0 2 年9 月1 3 日0 6 0 0 U T C 5 0 0 h P a 积冰预报2 0 0 2 年9 月1 3 日0 6 0 0 U T C 对流层顶高度2 0 0 2 年9 月1 3 日0 6 0 0 U T C 最大风速及高度预报1、颠簸指数LP指数法：当L0时预报有颠簸，计算出概率P7.2680.7180.3182.52uTuLznn0.5 911Lpe判断：当70%P50%，预报有轻度颠簸；当85%P70%，预报有中度颠簸；当P85%，预报有严重颠簸；特点：同时考虑了风场水平切变、垂直切变和温度场水平切变对
21、颠簸的影响1、颠簸指数诊断分析产品介绍判断：特点：颠簸指数Ellrod 指数 垂直风切变 VWS=散度 DIV=总形变 DEF=其中切变项DSH=拉伸项 DST=EIVWSDEFDIVVZ 根据EI判断：当8EI4时，预报有轻度颠簸；当12EI8时，预报有中度颠簸；当EI12以上时，预报有严重颠簸；()dudvdxdy22DSHDST()d vd ud xd y()dudvdxdy特点：主要考虑风场垂直和水平切变以及形变对颠簸的影响颠簸指数诊断分析产品介绍 根据E I 判断：特 20042004年年1111月月9 9日雷雨天气的分析日雷雨天气的分析个例1：2 0 0 4 年1 1 月9 日雷雨
22、天气的分析个例1：9日北京时0740-1042本场及周围阵雨，0746-0915雷雨，下午至夜间小雨，整个过程降水量9.8mm。天气过程简述9 日北京时0 7 4 0-1 0 4 2 本场及周围阵雨，0 7湿度稳定度指数（sweat index）700hpa的风场(低空急流）能量值降水量各气象要素的分析湿度主要以美国WAFS的2004年11月8日12UTC的数值预报资料来分析。本报告分析的各要素从预报雷雨有指示的角度出发，没有明显指示的要素没有分析。要素选择说明主要以美国WA F S 的2 0 0 4 年1 1 月8 日1 2 U T地面的相对湿度的分布表明地面湿度大低层1000HPA的T-T
23、d，能看出来有明显的水汽输送，温度场表现不明显低层1 0 0 0 H P A 的T-T d，能看出来有明显的水汽输送，温度7 0 0 h P a T-T d 与风场，7 0 0 h P a 存在着明显的低空急流T h e t a e 分布：高能舌区与雷雨的对应区雷雨出现在7 0 0 h P a 高能舌区雷雨出现在7 0 0 h P a 高能舌区对这次雷雨过程，S WE A T 指数大于3 0 0 的区域与雷雨区对应对流降水量与地面雷雨区对应好，而且从降水量看这次过程的预报与对流降水量与地面雷雨区对应好，而且从降水量看这次过程的预报与观测吻合好。对流降水量与地面雷雨区对应好，而且从降水量看这次过
24、程的预报7 0 0 h P a 的涡度区5 0 0 h P a 的涡度区总的说来，WAFS的数值预报产品对此次雷雨是有很好的指示意义的。以上只是对此次雷雨过程一些参数的描述，但不一定适合其它的雷雨过程。对不同天气形势下出现的雷雨，其气象参数及物理量的指示不一定相同。要针对不同的形势总结参数的使用，有的参数在这次过程有用，但在另一次过程中不一定有明显的指示作用。任何资料的使用都应积累经验，总结其优势。小结总的说来，WA F S 的数值预报产品对此次雷雨是有很好的指示日期时段降水性质系统降水量 5月月30日日 1420-1928UTC -TSRA冷涡冷涡 5月月31日日 0616-0715UTC1
25、052-1218UTC -TSRA，市区冰雹市区冰雹 TSRA 1110-1130B本场冰本场冰雹雹冷涡冷涡90日期S指数K指数雷雨时段5月月30日日08时时-2.223.95月月30日日20时时-7.016.91420-1928UTC，-TSRA5月月31日日08时时-2.825.10616-0715UTC，-TSRA，市区冰雹市区冰雹5月月31日日20时时-1.237.11052-1218UTC，TSRA 1110-1130B本场冰雹本场冰雹个例2:2 0 0 5 年5 月3 0 日冷涡雷雨日期时段降水性质系统降水形势介绍 5月30日08时从500hpa至850有冷涡，地面锋面，随着500
26、hpa的冷空气东移，有对流云生成东移，移动过程加强，中午之后管区开始有对流回波，回波一直维持，本场1420UTC至1908UTC雷雨；500hpa冷涡后部的高压脊向东伸，预示着冷涡后部有冷空气下渗，31日受低层辐合气流和冷涡的影响，0610UTC至0715UTC本场雷雨，市区冰雹；1052UTC至1218UTCB本场雷雨，其中1110UTC-1130UT冰雹。5 月3 0 日至3 1 日雷雨形势介绍500HPA的系统变化5 月3 0 日至3 1 日雷雨WA F S 5 0 0 H P A 的系统变化近地层辐合带的维持5 月3 0 日至3 1 日雷雨WA F S 近地层辐合带的维持稳定度指数5
27、月3 0 日至3 1 日雷雨WA F S 稳定度指数5 月3 0 日至3 1 日雷雨WA F SWAFS数值预报产品在航空气象预报及服务中发挥了很大的作用；进一步做好释用,使WAFS数值预报产品在航空气象预报中起到更大的作用；以此为初始场制作中尺度数值预报模式，为民航气象预报及服务提供更精细的产品。总结WA F S 数值预报产品在航空气象预报及服务中发挥了很大的作二、日本数值预报一天之中我们接收四个时次初始场的日本传真图，但业务工作使用的主要是00和12UTC两个时次的传真图资料。一天之中我们接收四个时次初始场的日本传真图，但业务工作使用的日本传真图的种类：5 0 0 h P a 高度、涡度预
28、报地面降水和风场预报5 0 0 h P a 温度和 7 0 0 h P a T-T d 度预报8 5 0 h P a 风、温和7 0 0 h P a 垂直速度预报4 8 小时地面气压、降水预报4 8 小时5 0 0 h P a 高度、涡度预报日本数值预报图上的主要内容所用的符号、单位、等值线的间隔、形状和有关说明 日本数值预报图上的主要内容所用的符号、单位、等值线的间隔、形符号符号内容内容单位单位等值线等值线说明说明间隔间隔形状形状HEIGHT(M)位势高度位势高度gpm60粗实线粗实线高低中心分别标注高低中心分别标注H、L，300hPa以上以上层次间隔层次间隔120gpmISOTRCH(KT
29、)等风速线等风速线n mile/h20虚线虚线在在300hPa以上层次用以上层次用PRECIP(MM)2436降水量降水量mm5虚线虚线标有中心数值，方括号内表示出预报时标有中心数值，方括号内表示出预报时段段P-VEL(hPa/H)P坐标垂直速度坐标垂直速度hPa/h10虚线虚线标有正负中心值，标有正负中心值，0值线用细实线，上值线用细实线，上升运动区（负值区）用阴影区表示升运动区（负值区）用阴影区表示SURFACE PRESS(hPa)地面气压地面气压hPa4实线实线高低中心分别标注高低中心分别标注H、L，北半球图上，北半球图上等直线间隔为等直线间隔为10hPaTEMP(C)温度温度3粗实线
30、粗实线冷暖中心分别标注冷暖中心分别标注C、W，300hPa以上以上层次直接标出各点温度值层次直接标出各点温度值T-TD(C)温度露点差温度露点差6细实线细实线以阴影区表示以阴影区表示T-TD3区区VORT(10*-6/SEC)涡度涡度10-6/s20虚线虚线标有正负中心数值，标有正负中心数值，0值线用细实线，值线用细实线，正涡度区用阴影区表示正涡度区用阴影区表示WIND ARROW风矢风矢直接用羽矢填在各计算格点上，规定同直接用羽矢填在各计算格点上，规定同常规天气图常规天气图符号内容单位等值线说明间隔形状H E I G H T(M)位势高度g p三、中尺度数值模式MM5非静力平衡中尺度模式MM
31、5是由美国国家大气研究中心和宾夕法尼亚州立大学联合开发的第5代中尺度天气预报模式，可以广泛用于大气科学研究，特别是在对中小尺度强对流系统、锋面、海陆风、山地环流和城市热岛等的理论和业务预报中有其独特的优势。非静力平衡中尺度模式MM5 是由美国国家大气研究中心和宾夕法尼民航北京气象中心从2005年12月开始引入MM5V3.7版本，主要目的是为民航气象预报及服务提供气象数值预报产品。民航北京气象中心从2 0 0 5 年1 2 月开始引入MM5 V 3.7 版本民航气象MM5模式目前预测范围为中国区域，中心位置中心位置设在首都机场（40.080111N，116.584556 E），水平四重嵌套，参数
32、如下：水平分辨率（KM）格点数 覆盖区域 D01 81 64*88 中国D02 27 82*73 华北D03 9 70*79 北京及其周边 D04 3 73*82 首都机场及 北京各区县B MC MM5 模式概况民航气象MM5 模式目前预测范围为中国区81km9km27km3km9km8 1 k m 9 k m 2 7 k m 3 k m 9 k m 模式垂直分为不等距的23个层。模式概况 模式垂直分为不等距的2 3 个 层。模式初始场初始场来自NCEP的1 1 全球分析场，侧边界条件侧边界条件使用NCEP预报场，每6小时换一次。每日0000、1200 UTC各作一次预报，积分48小时。模式概
33、况 模式初始场来自N C E P 的1 1 全球分析场INTERPF 模式运行 预报结果后处理 RIP4 绘图系统 系统 Vis5d 绘图系统 Grads 绘图系统 服务器 REGRIDDER TERRAIN PREGRID NCEP 资料 模式运行流程框架目前输出的产品有：高空：高度和涡度、风温、流线（层次：1000、925、850、700、500hPa）地面：海平面气压和风场、总降水、对流降水所有产品均每三小时输出一次！所有产品均每三小时输出一次！模式产品简介目前输出的产品有：0.1 0.2 0.4 0.8 1.6 3.2 6.4 12.8 25.6 51.2 106.4100110120130140oE203040No502006年年8月月1日日0000UTC6h 累积降水量累积降水量 36h预报预报实况实况MM5降水预报检验2 0 0 6 年8 月1 日0 0 0 0 U T C 3 6 h 预报实况MM5 降水预报8 5 0 h P a 风、温预告5 0 0 h P a 流线预报5 0 0 h P a 高度、涡度预报地面气压场、风场预报地面降水预报谢谢大家！
