2. 中央气象台,北京 100081
2. National Meteorological Centre, Beijing 100081
近年来,利用经向风指数 (或南风分量) 研究、预测季风雨带和旱涝位置变化,取得不少研究成果。例如,王亚非等[1]利用区域20°~30°N、110°~140°E及区域30°~40°N、110°~140°E的850 hPa南风分量之差,定义了一个反映梅雨现象的东亚夏季风指数,对长江中下游梅雨预报预测提供了分析和参考依据;姚才等[2]利用历年观测资料研究了6月同期东亚—西北太平洋环流与华南降水的关系,根据区域20°~26°N、105°~117°E 850 hPa南风分量,定义了一个包含了西南季风、副热带高压以及中高纬度西风槽等各影响系统信息的华南夏季风指数,表明季风指数强 (弱) 的年代与华南降水偏多 (少) 的年代有很好的对应关系;李向红等[3]利用1980—2001年NCEP和TBB资料,采用逐例与合成分析相结合的方法,分析了夏季风影响期间广西大范围暴雨发生前各层天气系统和孟加拉湾强对流云团的配置演变,发现暴雨前孟加拉湾强对流发展和经向风加强等特征,建立了广西季风暴雨中期预报平台;尹东屏等[4]利用平均场方法,研究暴雨时梅雨锋结构,发现经向锋生函数的正值区呈直柱状,与南风等风速线的密集区相重叠,纬向分布的锋生函数指示了暴雨区的范围;林爱兰等[5-6]利用经向风指数分析了1958—2000年南海对流层经向风转向 (从北风转为南风) 的早晚与华南春夏季逐月平均降水量时空分布特征,研究了广东省开汛迟早、华南前汛期锋面降水的多寡与南海夏季风爆发异常之间的联系,得到南海低层转南风有利于气流在华南地区辐合、形成降水的若干事实。以上成果表明,经向风指数在阶段性降水分析和预测方法上,有显著的应用价值。而气候学上阶段性降水是多个短期过程累加起来的,因此,从短期天气预报 (1~3 d) 角度考虑, 探讨利用经向风指数,预报短期暴雨过程开始时间、暴雨落区和暴雨强度是可行的,对提高短期暴雨预报准确率是有益的。基于这一思路,笔者在总结二十多年预报经验的基础上,根据华南锋面活动特点,设置广西暴雨预报关键区南北风指数,利用1970年以来的EC 850 hPa风场资料和广西降水资料,建立南北风指数计算系统,该系统在每日自动处理08和20时EC资料后,分别形成一个折线图和一个平面分布图,其中南北风指数折线图用于预报暴雨过程开始时间,南北风指数平面分布图用于预报暴雨落区和强度,在广西暴雨过程开始时间、暴雨落区、暴雨强度的预报中,取得较好的预报效果,对提高广西暴雨预报准确率有重要的参考价值,并已应用到日常预报业务中。
1 资料、方法和系统环境 1.1 资料(1) 暴雨资料取自广西气象台雨量数据库,从1970—2009年,时间长度共40年。
(2) 数值预报产品资料为EC 850 hPa 1970—2009年v分量实况和预报资料,空间分辨率2.5°×2.5°,时间分辨率为12 h (08和20 h),时间长度共40年。
1.2 计算方法所用方法包括统计分析、合成分析、离散格点坐标分析等。
1.3 系统环境研究目的是用于日常天气预报业务,重点是暴雨预报,因此从方便计算、推广、检验、运行稳定原则出发,研制、建立了“南北风指数计算系统”。
系统环境和设置:
(1) 硬件环境:联想电脑,Intel P4 1.7G,40G硬盘,125M内存,21′显示器,设置显示分辨率为1280×1024。
(2) 软件环境:Windows2000系统,VisualBasic 6.0高级编程语言。
(3) 具备ECMWF数值预报产品资料,内容包括:气压、位势高度、气温、风及湿度等。通过相应的计算增加u、v分量等产品。
层次:1000、850、700、500和200 hPa。
时效:0、24、48、72、96、120、144、168、192、216和240 h。
范围:0°~75°N、50°~140°E。
分辨率:2.5°×2.5°。
(4) 系统研制:系统主要是由资料路径设置,欧洲中心资料处理,计算u和v分量等模块组成。根据人机交互工作平台,计算关键区南北风指数,输出图像主要有一个折线图和一个平面分布图。
2 南北风指数计算方法广西地处华南西部,是我国多暴雨省区之一,一年四季均有暴雨发生,其中汛期暴雨发生次数明显多于非汛期。尽管目前针对华南暴雨的研究取得很多成果和进展[7-20],但是由于华南暴雨发生机理复杂,预报难度大,并且有效预报方法、预报工具缺少,预报准确率仍偏低。华南暴雨预报难点在于,降水机理与中高纬的江淮地区有明显差别,准地转理论基本上不适用,预报业务中的难点,主要体现在夏季气象要素分布接近均匀,风、压、温、湿等气象要素的水平梯度比温带地区小得多,常出现均压区,除热带气旋外,气压场和天气区分布的关系不明显,常常造成一些暴雨过程漏报或降雨量预报偏小。因此,在理论不能解决实际问题的情况下,总结经验归纳预报方法和工具,对提高暴雨预报准确率尤为重要。近年来数值预报能力显著提高,形势预报时效延长,在做低纬地区天气分析时,依靠低空风场 (流线和等风速线) 的变化更容易确定中小尺度天气系统,这正是暴雨预报的关键,因此设置暴雨预报关键区,量化分析低涡、切变线活动正是在这一思路下形成的。
2.1 关键区的选取关键区的选取基于冷、暖空气活动路径。为此选取20°~27.5°N、105°~110°E范围为广西暴雨预报关键区 (图 1),网格距2.5°×2.5°,分北风区、南风区两部分。北风区为27.5°N,105°~110°E的3个格点,西路、中路和东路冷空气影响时,出现正变压、负变温,统计为一次冷空气活动过程。南风区为20°~25°N、105°~110°E的9个格点,当偏南风活跃,尤其出现季风急流时,v分量迅速增大,为降水提供不稳定能量和水汽条件。
在27.5°N、105°~110°E范围,对850 hPa所有格点的v分量求和,定义为北风指数 (IN)。
${I_N} = \sum\limits_{i = 1}^n {v(i)} n = 3$ | (1) |
其中v(i) 是经向风分量,n是格点数。
在20°~25°N、105°~110°E范围,对850 hPa所有格点的v分量求总和,定义为南风指数 (ISW)。
${I_{sw}}= \sum\limits_{i = 1}^n {} \sum\limits_{i = 1}^n {v(i,j)} n = 3$ | (2) |
其中v(i, j) 是经向风分量,n是区域格点数。
当南风加强甚至出现急流时,与冷空气在关键区内交汇,便出现降水天气过程,根据南北风分量的大小,可以判断出有无日暴雨10站以上的暴雨过程。
3 预报指标的确定天气预报中,预报指标的设定主要基于天气出现概率的分布和预报员的分析经验、习惯等,本文中,概率分布通过以下方式计算:
(1) 采用南风指数和北风指数构建一个101×81的离散格点坐标 (图 2),采用式 (3) 计算每个离散格点下,南风指数和北风指数对应的暴雨事件发生的概率
${p_{{I_{sw}},{I_N}}} = \frac{{N + }}{{{N_ + } + {N_0}}}$ | (3) |
式中N+和No是落在圆形区域:{R(x, y), (x-ISW)2+(y-IN)2≤d2}的锋面过程发生暴雨事件次数和不发生暴雨事件次数。如果N+和No都为0,则该格点上记为缺省值。
(2) 将经过第一步计算之后缺省的网格用插值的方法获得取值,再对整个格点场进行平滑,图中实线为概率相同点连线,每隔0.1(即10%) 画一根。另外,南风指数和北风指数都相同的多次过程中,暴雨事件发生与否的记录 (+或o) 可能叠加显示在同一个位置。
将广西日暴雨达到或超过10个站点称为一次全区性暴雨事件。根据1970—2009年1115次锋面活动过程中,全区性暴雨事件发生时对应的南风指数和北风指数、全区性暴雨事件不发生时对应的南风指数和北风指数,计算得到全区性暴雨事件发生的概率 (等值线) 随南风指数和北风指数的分布。
图 2中显示了在1970—2009年40年间 (1—12月),广西受锋面过程影响,出现了440次 (天) 全区性暴雨事件 (+),675次 (天) 无全区性暴雨事件 (o),因此得出全年锋面活动发生全区性锋面暴雨事件的过程比例为39.5%。
图 2还显示了锋面活动有、无全区性暴雨事件发生时南、北风指数概率分布的全貌,根据等值线0.5覆盖区间分析,表明全区性暴雨事件发生主要集中在南风指数ISW在30~75 m·s-1、北风指数IN在-40~10 m·s-1区间。有两个大值中心,分别位于 (IN,ISW)=(9,58),(-25,50) 附近,两个中心之间的区域,暴雨事件落点密集。等值线0.5区间以外,南风指数和北风指数或偏大或偏小,表明全区性暴雨事件不容易发生。从单个样本分析,主要为非汛期 (10月至次年3月) 锋面过程,说明秋冬季北风偏大,降水系统停留时间短,降水总量偏小。而南风偏大区间为后汛期的8—9月,冷空气难到达华南,缺少动力条件,不容易出现暴雨。
虽然图 2反映年全区性暴雨事件发生的概率分布特征,但是华南锋面暴雨事件发生主要集中在4—7月,因此有必要再分析和了解4—7月全区性暴雨事件发生的概率分布特征。
图 3是4—7月全区性暴雨事件发生概率随南风指数和北风指数的分布。发生暴雨事件次数总共为359次,不发生暴雨事件次数总共为175次。等值线0.5覆盖区间,全区性暴雨事件发生同样集中在南风指数30~75 m·s-1、北风指数-40~10 m·s-1区间。而发生全区性锋面暴雨事件的比例上升为67.2%,证明4—7月是广西全区性暴雨事件发生的集中时段。
因此,全区性暴雨事件发生时,基于出现概率分布特征,南风指数大于30 m·s-1,北风指数小于10 m·s-1,是判断全区性暴雨发生的一个初步指标。
从预报员分析经验、习惯考虑,图 2和图 3反映的概率分布中,北风偏弱时,北风指数在10 m·s-1的附近,暴雨事件出现次数少反而概率高,是因为多数无冷空气影响事件并未真实反映在其中 (无:o),因此属于偶然性事件,代表性不强。同样,北风偏强时,北风指数≤-20 m·s-1区间,暴雨事件出现次数少也反而概率高,也属于偶然性事件,代表性不强。因此进一步综合考虑了不同月份关键区K指数、θse分布,热力条件和水汽条件,不损伤样本以及考虑了预报员分析习惯等因素后,通过经验订正,确定广西全区性暴雨发生的指标为:南风指数大于32 m·s-1,同时北风指数小于5 m·s-1。
3.1 暴雨过程开始时间预报指标暴雨预报涉及时间、空间变化,需要准确判断影响系统相遇的时间,才能减少空报、漏报。通常情况下,广西锋面暴雨发生前,暴雨预报关键区内,盛行偏南风。当冷锋过境时,850 hPa北风指数迅速减小到5 m·s-1以下,是暴雨过程开始的信号。
经统计,在1970—2009年440次锋面暴雨过程中,暴雨发生前24 h,南风通常迅速增大形成急流,冷空气进入广西北部关键区域时,与西南急流构成切变辐合区。
按照南风指数大于32 m·s-1,同时北风指数小于5 m·s-1的起报条件,在南北风指数折线图中,对440次出现日暴雨10站以上锋面暴雨 (328次) 过程进行大样本检验,符合这一条件的过程比例为74.5%,这一比例充分利用了样本总数,也属于一个比较合理的范围,因此,南风指数大于32 m·s-1,同时北风指数小于5 m·s-1,作为暴雨过程开始时间的预报指标。
3.2 暴雨落区预报指标初步确定了暴雨开始时间后,暴雨落区的判断非常关键。以往对暴雨落区的预报,通常考虑在切变线附近的急流轴左前方,主观判断成分大。通过对440次锋面暴雨过程南北风指数分析表明,暴雨通常出现在v分量0等值线南侧的南风急流核附近区域,南风急流核的闭合等值线为8 m·s-1或以上,大部分情况下,暴雨出现在8 m·s-1闭合等值线范围内。需要指出的是,暴雨预报需要考虑冷空气的路径,根据路径对暴雨落区作适当调整,比如当冷空气路径偏东,暴雨落区需要向左侧扩大至4 m·s-1等值线的部分范围。
利用南北风指数平面分布图,对440次锋面暴雨过程分析统计显示,在v分量南风8 m·s-1闭合等值线内出现暴雨的过程334次比例为75.9%。这一比例符合预报员的期望值,同时也属于一个比较合理的范围,因此,南风指数8 m·s-1闭合等值线可作为暴雨落区预报指标。
3.3 暴雨强度预报指标暴雨强度预报基于影响系统的强度及滞留时间,天气预报中往往根据高空槽的深浅、切变线、锋面的强弱来综合判断,主观性也较大。通过对南北风指数分析显示,当v分量南风区、北风区同时出现闭合等值线时,表明系统强度较大,容易出现大暴雨或特大暴雨。如果北风区不出现闭合等值线,则表明北风势力偏小,系统强度弱。
利用南北风指数平面分布图,对440次锋面暴雨过程的南北风指数统计分析显示,同时出现v分量南风区、北风区闭合等值线时,在南风12 m·s-1闭合等值线范围内,出现大暴雨或特大暴雨的过程 (295次) 比例为67%。其中v分量北风闭合等值线在-4 m·s-1或以下,大部分情况下,v分量北风比南风小,因此北风闭合等值线阈值定为-4 m·s-1,南风闭合等值线在12 m·s-1或以上。因此,同时出现v分量南风区、北风区闭合等值线,南风区出现12 m·s-1或以上闭合等值线,可作为暴雨强度预报指标。
4 预报应用实例2010年6月20日,广西出现一次超大范围暴雨过程 (图 4),日暴雨36站 (其中大暴雨9站),这次暴雨过程造成138万人受灾,因灾死亡1人,直接经济损失12.9亿元,达到特大型气象灾害等级。
经分析,这次大范围暴雨过程是在副热带高压减弱南落,高空槽东移,与南移进入广西的地面冷空气和850 hPa切变线相遇发生的。在这次暴雨过程中,利用基于v分量的南北风指数准确预报了暴雨开始时间、暴雨落区和暴雨强度,显示了较好指标意义。
4.1 暴雨开始时间的预报暴雨开始时间的预报,关键在南风大于32 m·s-1的同时,北风小于5 m·s-1的时间点。图 5显示,6月18日20时,广西关键区南风指数为72 m·s-1,北风指数为21 m·s-1,关键区盛行偏南风,没有达到大范围锋面暴雨发生的起报指标条件。
6月19日20时,关键区南风指数为69 m·s-1,北风指数为-1 m·s-1,同时达到起报条件,因此6月19日20时,是暴雨开始的时间点。此时500 hPa有高原槽东南移影响广西,地面冷锋南移进入广西,满足了大范围锋面暴雨发生的起报条件。实况暴雨开始时间是从6月19日20时开始,预报与实况吻合。
4.2 暴雨落区的预报暴雨落区预报的关键在分析南风闭合等值线变化范围。在这次暴雨过程中,冷空气自湖南和湖北进入广西,影响该地区,因此分析暴雨落区时,应该考虑向左侧扩大一些。图 6是6月19—20日08和20时v分量分布图,图中8 m·s-1的南风闭合等值线持续3个时次,南风急流核从广西北部向东南方向移动并减弱,最强时段出现在19日20时至20日08时,这也是暴雨发生的集中时段。北风-4 m·s-1闭合等值线持续4个时次,最强时段也出现在19日20时至20日08时。因此暴雨落区在南风8 m·s-1闭合等值线向西扩大到4 m·s-1等值线的范围内,实况暴雨落区与预报基本吻合。
暴雨强度预报的关键在分析v分量南风区、北风区是否同时出现强闭合中心。图 6中,前3个时次都同时出现南风区、北风区闭合中心,表明前3个时次系统强度大。6月19日08和20时,南风区域12 m·s-1闭合等值线内,最大值达到16 m·s-1。6月20日08时,v分量北风区最大值达到-15 m·s-1,表明从19日20时至20日08时,对流发展旺盛,有利于强降水的发生,达到大暴雨或特大暴雨降水量级。实况是出现9站100 mm以上的大暴雨。20日20时后,广西境内南风闭合中心 (急流核) 消失,北风控制,强降雨结束。
5 结论采用离散格点坐标方法,计算基于南北风指数对应全区性暴雨事件的概率分布,通过经验订正,结合南北风指数平面分布图,归纳出广西暴雨过程开始时间、暴雨落区和暴雨强度预报指标。
(1) 暴雨过程开始时间的预报,当南风指数大于32 m·s-1,同时北风指数小于5 m·s-1时,达到起报条件,利用EC 08和20时850 hPa风场预报资料,计算24~240 h的南北风指数,可提前作出对应时次的预报。
(2) 暴雨落区预报,根据EC 08和20时850 hPa风场预报资料,计算出v分量南风闭合等值线,暴雨主要出现在8 m·s-1闭合等值线区域。
(3) 暴雨强度预报,根据v分量南风区、北风区是否同时出现强闭合中心,在南风区是否有12 m·s-1闭合等值线,判断暴雨的强度。
王亚非, 张雁, 陈菊英, 2001. 一个能反映梅雨现象的东亚夏季风指数[J]. 气候与环境研究, 6(2): 146-152. |
姚才, 钱维宏, 林振敏, 等, 2011. 华南季风降水对应的环流指数[J]. 应用气象学报, 22(3): 302-311. DOI:10.11898/1001-7313.20110306 |
李向红, 蒋丽娟, 薛荣康, 等, 2006. 广西季风暴雨中期预报平台研制思路与应用[J]. 气象, 32(5): 110-115. DOI:10.7519/j.issn.1000-0526.2006.05.019 |
尹东屏, 张备, 孙燕, 等, 2010. 2003年和2006年梅汛期暴雨的梅雨锋特征分析[J]. 气象, 36(6): 1-6. DOI:10.7519/j.issn.1000-0526.2010.06.001 |
林爱兰, 李春晖, 郑彬, 等, 2010. 南海夏季风爆发与华南前汛期锋面降水气候平均的联系[J]. 中山大学学报 (自然科学版), 49(4): 134-138. |
林爱兰, 谷德军, 李春晖, 等, 2010. 南海夏季风爆发与华南前汛期锋面降水异常变化的联系[J]. 中山大学学报 (自然科学版), 49(5): 128-132. |
林良勋, 吴乃庚, 黄忠, 等, 2009. 广东2008年罕见龙舟水特点及成因诊断分析[J]. 气象, 35(4): 43-50. DOI:10.7519/j.issn.1000-0526.2009.04.006 |
包澄谰, 王两铭, 李真光, 1979. 华南前汛期暴雨的研究[J]. 气象, 5(10): 6-10. |
吴志伟, 江志红, 何金海, 2006. 近50年华南前汛期降水、江淮梅雨和华北雨季旱涝特征对比分析[J]. 大气科学, 30(3): 391-401. |
陈见, 钱俊, 高安宁, 等, 2006. 广西暴雨卫星云图与热带流场变化综合预报方法[J]. 气象科技, 34(2): 138-142. |
杨贵名, 毛冬艳, 姚秀萍, 2006. 梅雨期一次黄淮气旋发展的干侵入特征分析[J]. 热带气象学报, 22(2): 176-183. |
陈见, 钱俊, 黄明策, 2006. 广西暴雨过程的Ⅴ-3θ结构特征检验及预报应用[J]. 气象科学, 26(3): 334-340. |
高安宁, 陈见, 李生艳, 等, 2009. 两次西南低涡造成广西暴雨差异的对比分析[J]. 气象科学, 29(4): 557-563. |
黄忠, 吴乃庚, 冯业荣, 等, 2008. 2007年6月粤东持续性暴雨的成因分析[J]. 气象, 34(4): 53-60. DOI:10.7519/j.issn.1000-0526.2008.04.007 |
郑彬, 梁建茵, 林爱兰, 等, 2006. 华南前汛期的锋面降水和夏季风降水划分日期的确定[J]. 大气科学, 30(6): 1207-1216. |
伍红雨, 杜尧东, 秦鹏, 2011. 华南暴雨的气候特征及变化[J]. 气象, 34(10): 1262-1269. DOI:10.7519/j.issn.1000-0526.2011.10.009 |
伍志方, 曾沁, 吴乃庚, 等, 2011. 广州"5.7"高空槽后和"5.14"槽前大暴雨过程对比分析[J]. 气象, 37(7): 838-846. DOI:10.7519/j.issn.1000-0526.2011.07.008 |
黄明策, 2006. 广西暴雨时空分布特征[J]. 广西气象, 27(3): 9-13. |
吴兴国, 2001. 广西汛期暴雨若干特征分析[J]. 广西气象, 22(3): 9-12. |
吴兴国, 2001. 广西非汛期暴雨事件的特征分析[J]. 广西气象, 22(2): 7-9. |