1. 引言
低频天气图是一类新的预报的方法,它与常规的统计学预报、数值预报等方法有所差异,但又存在关联的预报10~30天延伸期的预报方法。通常将用滤波器滤波后得到的低频流场图、低频高度场图、低频降水量图等一些低频要素场图称为低频天气图。
近年来,气象灾害频繁发生的高影响天气事件越来越严重,尤其是城市内涝、泥石流、洪灾等极端强降水造成的灾害不容忽视,这对社会经济、农业领域、环境领域以及人民生活等方面产生严重影响,例如在2010年受持续强降水影响,长江流域多处发生超警戒洪水,中下游干流不同地点水位上涨,长江流域面临20年来的最高水位。2009年7月我国共有17个省份超过60个站点逐日降水量达到强降水的标准,其中四川较强降雨导致一些乡镇出现洪涝、山体滑坡、泥石流造成通讯、交通中断,部分居民住房严重损毁。2010年海南出现自1961年以来少见的强降水,20日的平均降水量超过1060.1毫米,平均达到暴雨级别的天数为6.6天,是1961年以来同期最多,导致部分江河水位超过警戒水位,出现严重内涝,部分中小学停课。随着社会经济,科学技术和国家建设的快速发展,政府和公众需要越来越高的气象服务,这对我们气象行业的要求也越来越高。我国气象部门的预报主要为常规预报和月尺度的气候预测,传统观念认为:日天气预报具有周时间尺度的上限,并且气候预测的下限为一个月时间尺度,并且它们之间存在不可预报的时间差。所以用低频天气图的方法来预报“天气预报”和“气候预测”两者10到30天“时差”的这种预报是具有研究性的。而低频天气图的建立,是基于季节内振荡的特性,通过滤波的方式从中选取低频信号,最后建立10到30天的天气预报,做好这时段的预报对人们的生活、社会的经济都有重要的意义。本文根据低频天气图方法,结合四川的天气情况,创建四川盆地夏季10~30 d强降水过程的预报方法,为政府部门的防灾减灾部署工作预留充足的救灾时间,那么将使防灾减灾工作更具有针对性,这会使公民的生命安全和财产得到保障,减少不必要的损失。
20世纪70年代初热带大气季节内震荡由Madden and Julian (1971; 1972)发现,后来以Madden和Julian的名字被称为Madden Julian Oscillation,是目前全球发现的最强低频信号 [1] [2] 。目前的研究成果是关于大气低频的研究主要集中在它的活动及其对热带的异常现象以及对全球气候系统的产生的影响,它的数值模拟(预报)研究,以及与厄尔尼诺与南方涛动事件的相互联系这三个方面 [3] 。由于热带MJO的准周期性、与热带对流的密切相关性以及其时间尺度与延伸期预报时效的匹配性,使得MJO成为了目前延伸期天气预报的重要桥梁 [4] 。当前针对热带MJO的研究,是用于延伸期预报,针对热带MJO的延伸期预报运用的是统计模型和动力模式两种方法,统计模型主要利用滞后回归模型、自回归模型、组合相似法和经验位相传播 [5] [6] [7] [8] 等方法来实时预报MJO,也有在多种统计方法相结合的基础上对MJO进行集合预报,如韩国利用三种统计学的理论即小波的周期分析、多元线性回归和奇异谱分析相结合的方式来预报对MJO,再与统计预报相结合,预报时效可以达到24天 [9] 。
基于统计预测模型,信飞根据矢量场的正交分解方法,自回归和多元回归的统计学方法建立了一个低频统计模型(30~50天),并创建了未来的5候的低频流场图 [10] 。丁一汇与梁萍基于大气低频信号和低频流场的统计模型,制作出梅雨季节的降水区的延伸期预报流程图 [11] 。陈丽娟等根据大气环流模式的特点对高程场进行了延伸期的扩展预报 [12] 。覃志年等通过延伸期预报方法对汛期的暴雨过程做了试验 [10] 。牛法宝等根据大气低频振荡的特征变化对不同季节的降雨做了延伸期预报的应用,并获得了很好的结果 [13] 。孙国武提出了运用低频天气图对延伸期降水的预报方法,利用低频天气图上的低频天气系统造成的南北气流辐合进行延伸期降水预报,后来进一步研究出了关键区大气低频预报方法 [14] 。20世纪80年代末,90年代初孙国武和章基嘉等就已经有了使用大气低频震荡来制作低频天气图预报中长期天气预报的设想 [15] ,并且和何金海教授等人共同将低频振荡预报方法应用于气象业务部门 [16] 。直到2009年左右,低频天气系统预报方法才被中国气象局作为一种新的方法进行试验研究,并将这种预报方法应用于气象业务预报中,在最近几年取得了突破性进展,其中上海、江苏、辽宁等地区应用了该方法在汛期强降水预报过程中 [17] [18] 。通过在不同区域2008年~2012年的应用之后,发现强降水过程10~45天的预报的效果是较好的。目前,低频天气图方法在各个气象部门进行推广应用,低频天气图方法作为延伸期天气预报的关键已经成为重点方法之一。
2. 资料和方法
2.1. 资料概况
本文利用资料有:
1) 四川2013~2017年5~8月逐日NCEP再分析数据,水平分辨率为2.5˚ × 2.5˚,垂直方向为,海平面,500 hPa,要素场主要包括U风场、V风场,通过处理将nc数据转换为dat格式,区域范围为0˚~70˚N,60˚~150˚E。
2) 2013年~2017年四川省155个站的逐日降水量。因为数据原因,最终选取132个站点,四川盆地132个站点分布图如图1所示。
2.2. 方法
本文通过对比不同滤波器的性能,最后选取Butterworth带通滤波器,它可以在序列内提取出物理量的周期,所以一直被广泛运用到气象问题的研究工作中,因为季节内振荡在30~50天有时间的周期性、空间的持续性和系统发展的准定常性 [19] ,所以最终选取滤波带通为30~50天。首先,利用低频天气图方法对降水过程做预报,用滤波后的数据建立低频流场图。其次,根据历史各例的统计分析及四川所处地理位置划分出关键区,分析低频气旋、低频反气旋等低频天气系统的活动周期和变化路径。最后分析低频天气系统与强降水过程的配置即强降水过程发生时的低频系统生成的区域、分布及连续出现的频数,

Figure 1. 132 stations in the Sichuan Basin
图1. 四川盆地132个站点分布图
找出最有利于四川强降水过程的低频系统配置,创建四川地区的预测模型,根据外推30~50天的方法,来预报自预报日之后10~30天中四川地区强降水的发生时段。
3. 低频天气图特征
3.1. 低频天气图的建立
四川省位于97˚21'E~108˚33'E和26˚03'N~34˚19'N之间,在长江的上游位置,居于中国西南,是青藏高原和长江中下游平原的过渡带之间,所有位势高低悬殊,导致川西地势高,川东地势低的特点十分明显,盆地西面海拔高度在3000 m以上,多为高原、山地等;盆地东面海拔较西部低,海拔在500~2000米之间,多为盆地、丘陵等。由于四川盆地地形的因素,低层的大气变化很复杂,所以本文选取500 hPa等压面对其进行分析研究。通过所选取的带通滤波器,滤波所用序列为U,滤波后的序列Un中第n个数据为:
其中:
式中,N为时间,
为滤波的开始时间和结束时间。
将2013~2017年5月1日到8月31日的500 hPa逐日u,v合成风场数据通过滤波器滤波,建立出逐日的低频天气图(如图2)。根据低频天气图中流场的分布特征,对低频天气系统的活动周期、低频天气系统的变化路径等参数进行研究,所谓低频系统及为低频气旋、低频反气旋等天气系统。
3.2. 低频天气系统活动周期
通过历史各例的统计分析,划分关键区的前提是不阻扰低频天气系统的区域性和完整性,为了更方便的分析低频天气系统的变化特征,在低频天气图的区域内划分出若干个区域。由于四川位于亚热带地区,按经纬度划将其划分为8个关键区(图3),1区:亚热带高压及台风活跃区(0˚~25˚N,120˚~150˚E):对水汽的输送起平衡的作用;2区:南海海域区域附近(0˚~25˚N,100˚~120˚E):水汽的来源地;3区:阿拉伯海、孟加拉湾附近(0˚~25˚N,60˚~100˚E):水汽的来源地;4区:东海至日本海附近(25˚~50˚N,120~150˚E):夏季风对降水的影响;5区:蒙古至河套附近(25˚~50˚N,100˚~120˚E):经向风、纬向风的作用;6区:青藏高原及以北地区(25˚~50˚N,60˚~100˚E):高原系统对降水的影响;7区:贝湖及以东地区(50˚~70˚N,100˚~140˚):低槽的影响;8区:贝湖以西地区(50˚~70˚N,60˚~100˚E):脊的强弱。
通过分析低频天气图,发现低频系统的移动很慢,但其移动的周期在每个关键区中有较大的差异,然而,低频天气图方法的关键又是低频气旋与低频反气旋在各关键区中的周期变化来决定。
本文通过普查2013~2017年5月~8月逐日低频天气图,归纳整理各关键区低频系统的中心位置、移动路径、活动周期等。从而找出低频系统在各关键区内的30~50天内的振荡规律,有时,先生成低频气旋在生成低频反气旋,又生成低频气旋或只生成低频气旋或低频反气旋,或同时生成低频气旋与低频反

Figure 2. 500 hPa low-frequency weather map (July 4, 2017) (A: low-frequency anticyclone, C: low-frequency cyclone)
图2. 500 hPa低频天气图(2017年7月4日) (A:低频反气旋,C:低频气旋)

Figure 3. Low frequency weather system critical area
图3. 低频天气系统关键区
气旋等……,各关键区内的低频系统并不同步。从而整理了各关键区低频气旋与低频反气旋的周期变化(表1),由表可以看出关键区1、3、6的振荡周期较短,而关键区4、5、7、8的振荡周期较长,其中关键区1、3属于低纬关键区,关键区4、5、7、8属于中高纬关键区。
3.3. 低频天气系统变化路径
低频天气系统的变化方式大致分为了两种及移动型(图4)和振荡型(图5),所谓移动型为低频活动周期内的低频天气系统在某一关键区产生,产生后向东西方向移动,或者向南北方向移动,最后移出产生时的关键区。振荡型路径则是指低频活动周期内的低频天气系统在关键区产生后就在原地振荡后自己消亡,始终没有移动出生成时的关键区。通过整理得到2013~2017年5~8月500 hPa低频天气图上各关键区低频天气系统的变化路径为:关键区1、关键区2、关键区3、关键区5、关键区6、关键区7生成的低频天气系统为移动型,关键区4、关键区8生成的低频系统为振荡型,其中关键区1、关键区2生成的低频系统为向西或东北方向移动;关键区3、关键区6生成的低频天气系统为向东北方向移动;关键区5、关键区7生成的低频天气系统为向东移动;关键区4、关键区8生成的低频天气系统大部分为原地振荡型。
3.4. 低频天气系统与强降水的配置
由于四川盆地地域较复杂,在该地区的降水分布以及降水强度有明显的区域性特征 [20] ,四川盆地降水主要集中在中部和东部地区,取2013年~2017年5月~8月逐日降水资料,设含有n个值的气象要素序列(这里为降水量)按升序排列后为X1,X2,…,Xm,…,Xn,则某个随机量值 ≤ Xm的几率可用下式来计算,即:P = (m − 0.31)/(n + 0.38),其中:m为Xm的序号;n为气象要素序列的个数;常值为经验

Table 1. Low-frequency synoptic oscillation period in key areas (unit: days)
表1. 各关键区低频天气图振荡周期(单位:天)
系数。第95个百分位值是P = 95%相应的Xm的值,定义四川强降水事件的阈值为第95 百分位的降水量。当某观测站的日降水量超出阈值时,规定该观测站产生了1次极端降水事件。通过计算,得到四川盆地第95 百分位的降水极端事件的阈值平均值为32.1 mm。盆地绝大部分的降水量都能达到20 mm,所以在实际业务中为了与天气预报更好的对应,将强降水阈值定为20 mm,此降水量为大雨级别。根据这个等级,制定四川地区夏季强降水过程的标准,当盆地地区内超过40个台站的日降水量 ≥ 20 mm,为一个强降水日。通过普查得到2013~2017年4月~8月四川盆地共出现强降水过程42次,其中2017年7月7日为一次强降水日,降水分布如图6。
分析500 hpa低频天气图,统计各关键区内低频气旋与反气旋生成的频数,得到当低频气旋在关键区1中活动时,强降水发生的次数为18次(18/42),当低频反气旋在关键区1中活动时,强降水发生的次数为12次(12/42);当低频气旋在关键区2中活动时,强降水发生的次数为13次(13/42),当低频反气旋在关键区2中活动时,发生强降水的次数为16次(16/42);当低频气旋在关键区3中活动时,强降水发生的次数为21次(21/42),当低频反气旋在关键区3中活动时,强降水发生的次数为14次(14/42;当低频气旋在关键区4中活动时,强降水发生的次数为14次(14/42),当低频反气旋在关键区4中活动时,发生强降水的次数为13次(13/42);当低频气旋在关键区5中活动时,强降水发生的次数为11次(11/42),当低频反气旋在关键区5中活动时,强降水发生的次数为6次(6/42);当低频气旋在关键区6中活动时,强降水发生的次数为18次(18/42),当低频反气旋在关键区6中活动时,强降水发生的次数为19次(19/42);当低频气旋在关键区7中活动时,强降水发生的次数为13次(13/42),当低频反气旋在关键区7中活动时,强降水发生的次数为19次(19/42);当低频气旋在关键区8中活动时,强降水发生的次数为23次(23/42),当低频反气旋在关键区8中活动时,强降水发生的次数为16次(16/42)。
4. 低频天气图在四川盆地的应用
4.1. 建立四川盆地夏季强降水低频预报模型
本文整理统计42个强降水日在8个关键区中低频系统出现的频数、活动周期、移动路径等特征,分析低频系统配置对四川盆地强降水的作用,制作出符合四川盆地强降水的低频预报模型(图7)。图中的A与C分别表示低频气旋和低频反气旋。由图可以得到,在关键区1、关键区3、关键区5、关键区8中有低频气旋活动时,关键区2、关键区7中会出现低频反气旋活动,而关键区4、关键区6中则有低频天气系统配合另外6个关键区活动,此时,来自北方的冷空气(箭头(1))与来自南方的暖湿气流(箭头(2))在四川盆地内交汇,这种形势很有利于强降水过程的发生,这也是低频天气图预报方法所阐述的重要的符合形势的模型。这种预报模型不仅考虑了南北气流相遇,还考虑了地形因素及局地小地形和川北的大巴山的抬升作用、西太平洋亚热带高压环流西侧的西南暖湿气流、西南低涡的发展、南海的东南气流、孟加拉湾西南气流、向东移动的高原系统、印度季风槽等实际的大气环流中关键的天气系统对盆地降雨的影响。这样使得关键区内的低频天气系统与实际强降雨天气过程有相应的配合关系,证明了使用低频天气图进行10~30天的天气预报是可以应用降雨过程的。
4.2. 强降水过程的预报流程
本文通过四川盆地强降水低频预报模型,整理出强降水的预报流程(图8),第一步为资料的选取,主要取区域流场等格点资料,通过滤波,绘制出低频天气图,然后根据四川盆地的地理条件及位置确定关键区,分析关键区内低频系统及低频系统与强降水之间的关系,最后建立强降水的低频预报模型,确定关键区中低频系统生成及周期,做出相关预报。
4.3. 2017年四川夏季强降水过程预报应用
受青藏高原系统的移动、西南涡的发展、南海东南气流和西南季风等因素的影响,四川盆地汛期主要发生在5月~9月,本文根据低频天气图方法,建立四川地区强降水低频预报模型,主要将2017年夏季6月~7月的强降水过程应用于此预报模型中,对其进行预报试验(表2)。由于季节内振荡有30~50天的时间周期性、空间持续性和系统生成的准定常性,所以在每次预报时需要分析前30天的逐日低频天气图,整理记载各个关键区中低频系统的生成时间、变化路径,切合预报模型中与低频系统变化一致的日

Figure 6. Sichuan Province, July 7, 2017 heavy precipitation daily precipitation distribution map
图6. 四川省2017年7月7日强降水日降水分布图

Figure 7. Low-frequency prediction model of heavy precipitation in Sichuan Basin (C: low-frequency cyclone, A: low-frequency anticyclone, arrow (1): cold air in the north, arrow (2): warm south air in the south)
图7. 四川盆地强降水低频预报模型(C:低频气旋,A:低频反气旋,箭头(1):北方冷空气,箭头(2):南方暖湿空气)
期,依据统计出低频系统在各关键区的周期,来外推未来10~30天中有与低频预报模型中低频系统配置符合相对应的时段,这个时段则是发生强降水可能性最大的时间。
本次试验预报降水过程13次,预报过程中正确有8次,特别是7月5日~7日预报十分准确,达州站特大暴雨(如图9)带来很大损失,提前25天预报出来将减少一些损失,漏报2次,空报3次,发现正确率为61.53%,同时发现空报次数较多,这是因为低频天气图的周期变化的原因,3次空报的时间的实际低频系统的变化比预报时间推迟10~20天,所以实际应用时对低频系统的变化的分析不能忽视,2次漏报分别是6月下旬与8月上旬的两段时间内,其中6月30日至7月2日为漏报,这时间段甘孜内有暴雨的产生,若漏报,将会产生一些气象灾害,给公众的生活带来一些损害,所以要对模型前期预报的分析结果进行修订,这能大大提高预报的准确率。总的来说,低频天气图方法用于四川盆地的强降水过程

Figure 9. Precipitation in Rida Prefecture with heavy precipitation
图9. 强降水日达州降水量分布图

Table 2. Prediction of heavy precipitation process in Sichuan Basin from June to July 2017
表2. 2017年6~7月四川盆地强降水过程预报
的预报结果是比较成功的。通过分析还发现预报试验的预报结果和实况来对比,发现运用低频天气图方法对预报6月份的强降水过程效果要更好一些。
4.4. 低频天气图的天气学意义
低频天气系统与常规天气图上的天气系统存在较大的差别,低频天气系统所反映的是某些天气过程中天气系统的发生、维持、变化的过程。通过分析2017年某一强降水日的低频流场图(图10(a))和实况流场图(图10(b))。发现,低频气旋和低频反气旋的生成的位置在四川盆地所在的30˚N附近汇合,相应的低频偏北降水过程的气流与低频偏南气流在四川盆地的分布与其低频系统分布相对应,同时期实况流场图上的变化也是如此,从实况流场图上也可以得到偏北气流与偏南气流在四川盆地所在的30˚N附近汇合。低频天气图上偏南气流可以反映暖而湿的气流北上,两支气流分别来源于东亚季风、印度季风;而偏北气流则反映出中高纬区域的冷空气南下。冷空气与暖湿空气在该地区汇合,是造成降水过程的重要因素。虽然低频天气系统与实况天气系统有所差异,但两者的环流系统是相互关联的。低频气旋对应的环流系统为实况天气图上的气旋,低频反气旋所对应的环流系统是实况天气图上的反气旋。然而低频天气系统有30~50天低频振荡,其特点是有明显的时间周期性、空间持续性和系统生成的准定常性,在实况的天气图上的大气环流并不没有这些时空特点,所以在30天以前做不出相关预报。这也说明了低频天气系统对于大气环流系统的变化也起到相应的作用。
以四川地区为例,来分析低频气旋、低频反气旋与实况环流系统的联系,2017年四川大部分地区暴雨时间在6月下旬至7月上旬,如图10(a)、(b)可得:在蒙古至河套地区有一低频反气旋,而相应的在四川地区有相对应的低频的偏南气流相对应。随着蒙古至河套的低频反气旋向南移动,在阿拉伯海、孟加拉湾地区生成低频气旋。低频气旋与低频反气旋在四川盆地28˚N附近汇合并停留,导致暴雨的形成。
4.5. 完善技术
任何新预报技术方法都不是固定的,需要通过实践检验,在实践中改进、充实和完善,甚至借鉴吸吸取其他方法的经验,才能为之所用,下面从三个方面进行探索与思考:
1) 低频天气图方法有很明显的地域性,适合在区省、市的气象单位进行应用推广。但应用时,需要着重考虑地域的大气环流、天气系统、地形地理等方面的特点,把这些特点加在低频天气图方法中,就会得到比较有地域性的低频天气图方法。
2) 低频天气图方法要以天气动力学原理为基础,把实际天气图(未滤波)上的高频变化的天气系统滤去,改为具有缓变特征的低频天气图上的低频系统。
3) 低频天气图方法预报流程中,有几个预报流程(关键区的选择、低频图的分析和预测模型的建立) 非常重要。如分析低频天气图,作为预报人员要对天气学原理和技术规则十分了解,能够读懂低频天气图,了解低频系统的演变特征,这是低频天气图方法的基本功。
5. 结论
本文建立逐日低频天气图,普查低频系统的周期,分析低频系统的移动路径来阐述低频天气图的特征,在联系与强降水的配置关系,得到低频天气图上低频系统生成位置,与实际大气环流中有利于产生降水的天气系统位置相对应。这说明通过低频天气图方法来预测夏季10~30天强降水的产生,在天气学方面是有意义的。
确定符合四川盆地夏季强降雨时段的低频系统的关键区,分别为亚热带高压及台风活跃区(1区)、南海海域区域(2区)、阿拉伯海、孟加拉湾附近(3区)、东海至日本海附近(4区)、蒙古至河套附近(5区)、青藏高原及以北地区(6区)、贝湖及以东地区(7区)、贝湖以西地区(8区),其中关键区1、3、6的振荡周期较短,而关键区4、5、7、8的振荡周期较长,其中关键区1、3属于低纬关键区地区,关键区4、5、7、8属于中高纬关键区。

Figure 10. June 23, 2017 Flow Field Map ((a) Low Frequency Flow Field Diagram, (b) Live Flow Field Diagram)
图10. 2017年6月23日流场图((a)低频流场图,(b)实况流场图)
分析四川盆地强降水与低频系统配置的关系,制作出四川盆地低频天气图的预报模型。关键区1、关键区3、关键区5、关键区8中有低频气旋活动时,关键区2、关键区7中会出现低频反气旋活动,而关键区4、关键区6中则有低频系统配合另外6个关键区活动,此时最有利于强降水过程的发生。这不仅考虑了南北气流相遇的原因,还考虑了不仅考虑了南北气流相遇,还考虑了地形因素及局地小地形和川北的大巴山的抬升作用、西太平洋亚热带高压环流西侧的西南暖湿气流、西南低涡的发展、南海的东南气流、孟加拉湾西南气流、向东移动的高原系统、印度季风槽等实际的大气环流中关键的天气系统对盆地降水的影响。
最后利用低频天气图方法,对2017年夏季6月~7月强降水过程进行预报应用,得到的结论是,预报降水的正确率为61.53%,同时发现空报次数较多,这是因为低频天气图的周期变化的原因,3次空报的时间的实际低频系统的比变化预报时间推迟10~20天,所以实际应用中对低频系统的变化的分析不能忽视,并且要对模型前期预报的分析结果进行修订,这能大大提高预报的准确率。总的来说,低频天气图方法用于四川盆地的强降水过程的预报结果是比较成功的。通过分析还发现预报试验的预报结果和实况来对比,发现运用低频天气图方法对预报6月份的强降水过程效果要更好一些。
致 谢
本文是在肖天贵老师的指导下完成的,不管是开始定题目还是在查资料准备以及撰写论文的全过程中,肖老师一直都耐心地给予我指导和意见,使我在思考问题及撰写论文方面都有了很大提高;在每次设计遇到问题时老师不辞辛苦的讲解才使得我的论文顺利的进行。在此肖老师表示衷心地感谢!还要感谢和我同组的同学们,在撰写论文的过程中,遇到的一些困难,正是你们的帮助,才得以顺利解决,在此表示深深的谢意。
基金项目
国家自然科学基金(91337215,41575066)、国家科技支撑计划(2015BAC03B05)、全国大学生创新训练计划项目(201610621033)支持。
参考文献