气象学报  2013, Vol. 71 Issue (5): 825-838   PDF    
http://dx.doi.org/10.11676/qxxb2013.069
中国气象学会主办。
0

文章信息

徐洪雄, 徐祥德, 陈斌, 陈联寿, 朱复成. 2013.
XU Hongxiong, XU Xiangde, CHEN Bin, CHEN Lianshou, ZHU Fucheng. 2013.
双台风生消过程涡旋能量、水汽输送相互影响的三维物理图像
The structure change and energy moisture transport physical image in the development and decay processes of binary typhoon vortices
气象学报, 71(5): 825-838
Acta Meteorologica Sinica, 71(5): 825-838.
http://dx.doi.org/10.11676/qxxb2013.069

文章历史

收稿日期:2012-12-10
改回日期:2013-03-29
双台风生消过程涡旋能量、水汽输送相互影响的三维物理图像
徐洪雄1, 徐祥德2 , 陈斌2, 陈联寿2, 朱复成1    
1. 南京信息工程大学, 南京, 210044;
2. 中国气象科学研究院灾害天气国家重点实验室, 北京, 100081
摘要:重点研究“莫拉克”台风发展并登陆台湾,以及“天鹅”台风消亡阶段两者相互作用的问题。通过诊断分析发现“莫拉克”与“天鹅”移动过程存在双台风涡旋互旋、吸引与合并现象。采用双台风中心连线的垂直剖面移动坐标分析法可揭示出双台风涡度、风场三维结构,演变过程中双台风的涡度、动能强度呈反向变化关系,在双台风生消过程中,动能、位涡场分布存在显著“连体”通道特征。并揭示出双台风涡旋各自生、消过程水汽、动能可能存在的相互影响及其涡旋结构变化的内在关联。对“天鹅”消亡、“莫拉克”引发暴雨过程,采用Flexpart-WRF耦合模式模拟“质点群”轨迹,模拟结果再现了双台风生消阶段“天鹅”台风水汽“粒子群”向“莫拉克”低层气旋式输入通道,且在“莫拉克”涡旋高层反气旋式卷出的三维立体动态图像。通过剔除“天鹅”台风涡旋数值模拟试验进一步印证了“天鹅”台风趋于消弱过程其水汽、动能输送为“莫拉克”台风发展与维持做出了一定贡献。基于以上合成分析、轨迹和数值模拟技术综合分析提出了能揭示“天鹅”消亡、“莫拉克”发展过程能量、水汽输送相互影响的三维物理图像。
关键词双台风相互作用     水汽通道     轨迹     三维结构    
The structure change and energy moisture transport physical image in the development and decay processes of binary typhoon vortices
XU Hongxiong1, XU Xiangde2 , CHEN Bin2, CHEN Lianshou2, ZHU Fucheng1    
1. Nanjing University of Information Science and Technology, Nanjing 210044, China;
2. State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China
Abstract:A possible mechanism responsible for the exceptional rainfall over Taiwan by Typhoon Morakot (2009) (named Morakot simply here after) was explored by the means of the investigation of the binary typhoon vortices (BTV) interaction. The analysis of the three-dimensional thermal-dynamic structure revealed that there is substantial moisture, energy transport from Typhoon Goni (2009) to Morakot. The inherent correlation between the two vortex respective development and decay is also obvious, indicating aoccurrence of mutually rotating. The dynamic vertical structure of the two vortexes presents an opposite tendency, namely, the intensity of Typhoon Goni vortex becomes weaker and asymmetric with the time elapsed, which is in contrast to the continuous development of "Morakot" during the period of the BTV interaction. However, the horizontal structure of energy and potential vorticity distribution of the two vortexes is characterized by obvious connecting body transportation. A significant moisture transport channel from Typhoon Goni to Morakot can be deduced from the further analysis of the correlation between the satellite blackbody temperature (TBB) over the Taiwan areas and the moisture column flux around its surrounding areas. Using the trajectories model Flexpart coupling with a mesoscale weather forecasting model WRF, we reproduced this moisture transport pathway by the backward tracking of majority of air parcels initialized over the cores of Goni. In addition, the numerical simulation experiments with the Goni vortex being removed or not further confirmed that the moisture and energy transportation from Goni is one of the key factors for the development and maintenance of Morakot. Based on the above data diagnosis, trajectory analysis and numerical simulation, a three-dimensional image of BTV interaction is provided to describe the energy, moisture transport of BTV, and the development and decay of BTV.
Key words: Interaction between binary typhoon vortices     Moisture transport channel     Trajectory     Three-dimensional structure    
1 引 言

双热带气旋活动是造成复杂热带气旋路径的重要原因之一,双热带气旋的相互作用可导致两者的气旋性互旋、相互吸引、排斥或者热带气旋的停滞打转等现象,从而使路径变得复杂,这往往给预报工作带来很大困难(陈联寿等,1979)。因此,双热带气旋的相互作用一直是大气科学界研究的一个重要课题,也是台风研究领域十分关注的理论难点。

早在20世纪20年代初,Fujiwhara(19211923)就对一般的双涡旋问题进行了理论分析和实验考察,发现了双涡旋相互作用的一些现象,这些现象后来被称作藤原效应。半个多世纪以来,许多学者对双台风的相互作用曾作了不少研究。例如,Br and (1970)根据22个个例研究了北太平洋的双热带气旋,指出双热带气旋的旋转依赖于其中心距离,发现相互作用明显的距离在700海里以内。陈联寿等(20012002)指出,一般情况下当两个台风之间的距离足够近时,两个强的低值系统会产生吸附作用,促使两个台风有合并的趋势,合并成一个环流系统,从天气学的角度提出了台风强度变化的一种新的可能机制。相邻中尺度涡旋与台风涡旋的相互作用可以使台风加强。此外,学者们通过数值模拟(罗哲贤等,20012002Wu et al,2003Prieto et al,2003; Yang et al,2008)探讨了双台风的相互作用对于台风路径、强度以及台风暴雨的影响。

本研究选取引发异常暴雨的“莫拉克”(Morakot)台风作为研究双台风相互影响的重点目标。2009年第8号热带气旋“莫拉克”给中国台湾省带来一场历史罕见的灾害。台湾南部遭遇了50多年以来最强的一次暴雨天气过程,全岛过程雨量普遍在500—1500 mm,而台湾中南部许多县市的过程雨量在1500—2500 mm。如此之强的罕见特大暴雨造成台湾中南部地区洪水暴发,泥石流等次生地质灾害泛滥,是台湾气象史上伤亡最惨重的台风。受“莫拉克”正面袭击,福建、浙江等地的沿海地区还普降暴雨到特大暴雨,累积雨量100—500 mm(台湾环境保护局(Environment Protection Administration(Taiwan)),2010;Chanson,2010)。“莫拉克”台风造成台湾与大陆沿海异常灾害的成因是学术界十分关注的热点话题,“莫拉克”台风为何登陆台湾后涡旋持续增强?而且台风涡旋为何存在如此巨大的能量?其涡旋系统维持的强对流及其强水汽输送的机制究竟是什么?更值得注意的是“莫拉克”台风持续增强过程中位于其西侧的“天鹅”(Goni)台风却在与“莫拉克”台风互旋过程中消亡。因此,“天鹅”与“莫拉克”台风生消过程是否存在相互关联?即双台风生、消过程是否存在着能量、水汽相互输送机制?关于“莫拉克”与“天鹅”双台风相互作用问题,已有研究(Xu et al,2011)表明,“天鹅”台风对于“莫拉克”台风强度及其产生的暴雨有重要贡献,“莫拉克”与“天鹅”双台风相互影响问题,有关文献从双台风相互作用不同阶段的特征提出了具有学术价值的观点,如Wu等(2012)认为,双台风相互作用前期夏季西南水汽流与台风涡旋相互作用亦是重要的影响因子,其模拟结果表明,前期“莫拉克”对“天鹅”台风发展具有影响。但在“天鹅”台风消亡过程中,“莫拉克”台风涡旋却发展并引发异常暴雨灾害,后期双台风相互作用仍是值得研究的问题。本研究重点在于后期,即“天鹅”台风消亡,莫拉克台风引发灾害的阶段。

针对上述双台风相互影响的理论问题,本研究以双台风的三维结构为切入点,采用资料合成分析、数值模拟以及拉格朗日轨迹计算相结合的技术途径剖析以下有关双台风相互作用问题:(1)探讨“天鹅”与“莫拉克”双台风互旋、吸引过程台风涡旋三维结构变化趋势特征;(2)揭示双台风生、消过程涡旋间的能量、水汽输送途径以及其相互影响;(3)研究双台风相互影响过程中能量、水汽交换与输送三维结构及其综合物理图像,深化认识双台风一个台风消亡和另一个台风发展并引发灾害过程的关联机制。2 资料与模式2.1 资 料

采用中国气象局上海台风研究所6 h一次的最佳路径资料,包括台风中心位置(经度、纬度)和强度(最大风速和台风中心气压)等信息研究双台风互旋过程。采用美国国家环境预报中心(NCEP)/美国国家大气研究中心(NCAR)的再分析资料(时间分辨率为6 h,水平分辨率为1°×1°,垂直方向1000—10 hPa共26层)来认识双台风间的能量、水汽交换、转化过程,欧洲中期天气预报中心(ECMWF)的高分辨率(25 km)YOTC(Year of Tropical Convection)全球分析资料来分析移动过程中双台风涡旋动力、热力内部结构的时空演变,台湾省的降水观测资料则来自其台风数据库(2009)。

2.2 整层水汽

众所周知,单层的水汽通量无法清楚地表达空间的水汽输送情况,为了追踪水汽源,本研究采用NCAR/NCEP再分析1°×1°资料,分析陆面或海面到300 hPa高度整层大气的水汽及其输送特征,其中,整层纬向水汽输送Qλ,经向水汽输送Qφ的计算方法如下:

式中,u、v分别为纬向风和经向风分量,q为比湿,ps为地面气压。

2.3 整层水汽相关矢量

为进一步认识“莫拉克”台风引发暴雨过程的水汽源及其水汽输送结构,采用相关矢量方法,以追踪台湾地区水汽源,其相关矢量的数学模型为

式中,Rqu和qv分量的相关合成矢量,Rx(x,y)、Ry(x,y)分别为整层水汽qu、qv分量相关场。

2.4 数值模式

为了验证“天鹅”台风与“莫拉克”台风的相互作用,采用新一代中尺度数值天气预报模式(WRF)进行数值模拟。模拟时间为2009年8月7日00时—9日00时(世界时,下同),使用NCEP 1°×1°再分析资料作为初始场及侧边界驱动场,网格选取两重嵌套区域,区域中心为(25°N,125°E),粗网格水平方向为163 ×124个格点,水平格距为30 km,细网格水平方向为160 ×127个格点,水平格距为10 km,垂直分为28 层,模式顶高50 hPa,时间步长180 s,输出时间间隔1 h。微物理过程内外层均采用Lin方案(Lin et al,1983),长波辐射采用RRTM方案(Mlawer et al,1997),短波辐射采用Dudhia方案(Dudhia,1989),积云对流方案采用Kain-Fritsch对流方案(Kain,2004)。

2.5 粒子扩散模式

粒子扩散模式Flexpart(Lagrangian Particle Dispersion Model)是由挪威学院大气系空气研究所和气候研究所开发(Stohl et al,2005),它是一个拉格朗日粒子扩散模式,可以用来追踪大气成分轨迹。Flexpart模式能够计算大量的粒子群轨迹,用这个轨迹来描述示踪物在大气中的传输和扩散过程。而Flexpart-WRF耦合了WRF和Flexpart两个模式,使其可以适用于中小尺度系统。在此采用Flexpart-WRF粒子扩散模式来计算“天鹅”台风粒子群轨迹,采用WRF模拟逐时的输出作为驱动场,粒子群初始时刻放置于“天鹅”台风区域,其范围为(18.1—18.9°N,108.1—110°E),高度为100—3000 m,模拟时间为2009年8月7日00时至9日00时。3 双台风互旋过程能量、水汽相互影响通道特征

“天鹅”与“莫拉克”台风路径(图 1)表明,两者路径轨迹走向均为东南—西北向,登陆地点分别在台湾至福建省和海南岛;2009年8月7—10日“莫拉克”台风由原偏西路径折向北,与此同时“天鹅”台风6—9日亦出现路径“打转”异常现象,即由原向西北方向行进,6—7日出现折向西南然后又转向东南的路径。分析“天鹅”与“莫拉克”6—9日相对位置发现两台风涡旋相互“靠拢”,且呈现出双台风反时针“互旋”现象。

图 1 “莫拉克”(蓝色)和“天鹅”(红色)台风最佳路径及“互旋”现象

(绿色线为同时刻两台风间的连线)
Fig. 1 Best tracks of Typhoon Morakot(blue line) and Typhoon Goni(red line)with the line(green line)connected the two TC centers

从“莫拉克”与“天鹅”台风区域中层(500 hPa)涡度及中高层(300 hPa)台风“暖心”温度变化(图 2)可见,2009年8月7—8日,“莫拉克”处于维持加强阶段,其中层涡度加强,中高层温度上升,表明暖心的发展,8日起“莫拉克”涡旋“暖心”持续加强,而“天鹅”则处于消弱趋势,且7—9日其中层涡度明显减弱,中高层暖心温度明显下降,两者涡旋强度呈反向变化。

图 2 NCEP再分析资料“莫拉克”与“天鹅”台风中心的500 hPa涡度(a)及300 hPa温度(b)变化 Fig. 2 Temporal variations of the central vorticity at 500 hPa(a) and the central temperature at 300 hPa(b)for Morakot and Goni

为进一步研究双台风水平流场及其水汽输送的相互影响问题,从台风风速、涡度等物理量和水汽相关矢量的水平结构,探讨双台风相互作用过程中的输送“通道”,并揭示双台风“互旋”运动规律。考虑双台风为互旋移动过程,计算双台风涡旋三维动力结构时采用移动坐标,即以双台风中心间连线作为风场、涡度场垂直剖面图的坐标横轴。结果表明(图 3a、b),在“天鹅”、“莫拉克”双台风中心连线的垂直剖面上,7—8日双台风中心两侧风速随时间变化呈准对称垂直结构特征。“天鹅”台风在7日台风中心两侧存在明显的风速轴,8日台风中心两侧风速轴接近消失,且风速减小,而“莫拉克”台风两侧前后两时刻均维持显著的准对称风速轴,且7—8日西侧最大风速从30 m/s增大至35 m/s,东侧则从40 m/s增大至50 m/s,以上分析表明双台风涡旋风场垂直结构变化亦呈相反变化趋势。

图 3 “天鹅”和“莫拉克”两台风中心连线“斜线”风速(由YOTC资料计算,单位:m/s)垂直剖面

(a.2009年8月7日00时,b. 8日00时;阴影区为大于14 m/s的区域,虚框区为台风区域)
Fig. 3 Wind speed cross section along the line connected the two TC centers calculated from the YOTC data(dashed line denotes the TC area)(a. 00:00 UTC 7 August and b. 00:00 UTC 8 August; the shaded is for the area in which the speed >14 m/s; the dashed frame is for the typhoon area)

分析双台风强涡度柱结构变化特征(图 4)可见,双台风强涡度区的垂直结构变化表现为,7—8日“天鹅”台风由较强涡旋柱(涡度极值为50×10-5s-1),转变为弱涡旋柱(涡度极值小于30×10-5s-1),“莫拉克”则相反转变为强涡旋柱(涡度极值超过80×10-5s-1)。从以上7—8日双台风强涡度区垂直剖面图分析,可发现双台风涡旋强涡度区呈反向变化趋势,即7—8日“天鹅”台风涡旋垂直动力结构发生剧变,呈消弱特征,而“莫拉克”台风涡旋垂直结构则显著加强。

图 4图 3但为涡度剖面(阴影区涡度大于10×10-5s-1) Fig. 4 As in Fig. 3 but for vorticity with its values greater than 10×10-5s-1 shaded

图 5给出“天鹅”、“莫拉克”双台风2009年8月5日12时—9日12时850 hPa水平风场结构变化过程,5日12时,“天鹅”与“莫拉克”台风涡旋为相对“独立”单体;6日12时—7日12时,双台风风场涡旋结构间气流有“连体”通道(图 5,黑色流线与阴影区所示),值得注意的是,此阶段双台风在“互旋”过程中两者距离缩小,出现“相互吸引”现象(图 5,虚线所示);7日12时,“天鹅”台风涡旋削弱,涡旋准对称性已破坏;至8日12时,“天鹅”台风涡旋处于消亡过程的残涡状态,但原“天鹅”向“莫拉克”动能通道还持续维持(图 5 d、e),而“莫拉克”台风涡旋呈加强趋势;9日12时,“莫拉克”台风引发特大暴雨后在中国东南沿海地区登陆。以上台风涡旋的生、消过程亦反映了“连体”双台风涡旋能量相互影响的特征,亦可揭示出“莫拉克”台风涡旋发展,强度增强可能与“天鹅”台风消亡过程中的能量输送、转化有关。

图 5 2009年8月5日12时—9日12时(a—e)间隔24 h的NCEP 1°× 1°再分析资料500 hPa风场(A、B分别代表“天鹅”和“莫拉克”台风涡旋;阴影区动能大于100 m2/s2,间隔为200 m2/s2;实线为动能通道,虚线为双台风中心连线) Fig. 5 Wind fields(vectors)with kinetic energy greater than 100(m2/s2)(color-shaded)at 12:00 UTC from 5 to 9 August 2009(a-e),respectively,calculated from the NCEP/NCAR 1°× 1°reanalysis data(Solid line is kinetic energy channel,dashed line denotes the line link Morakot(B)with Goni(A))

分析“天鹅”与“莫拉克”双台风涡旋位涡场(图 6)的分布可知,2009年8月5日12时—9日12时“天鹅”(A)、“莫拉克”(B)台风区域位涡高值中心亦呈“互旋”现象(图 6b、c、d箭头所示),在5日12时(图 6a)、6日12时(图 6b)双台风涡旋各自呈“分离”的位涡环形高值中心,7日12时(图 6c)、8日12时(图 6d)“天鹅”与“莫拉克”涡旋间出现带状双台风位涡“连体”通道。另外,随时间演变双台风位涡高值区距离趋于缩小亦呈“相互吸引”状态,其中,8日12时“莫拉克”台风位涡高值环形圈范围扩大,极值增强超过1 PUV(图 6 虚线椭圆长轴形变及其高值区变化),而天鹅位涡高值环形圈则呈收缩状态。9日12时“天鹅”台风位涡高值区显著消弱,且其位涡高值环形圈消失。

图 6 2009年8月5日12时—9日12时(a—e)间隔24 h的500 hPa位涡(单位:PUV)分布

(阴影间隔0.4 PUV; 椭圆形(虚线)为双台风位涡高值中心间覆盖范围,粗箭头为台风位涡移动方向)
Fig. 6 Potential vorticity greater than 0.4(PUV)with the interval of 0.4 at(a)12:00 UTC 5,(b)12:00 UTC 6,(c)12:00 UTC 7,(d)12:00 UTC 8 and (e)12:00 UTC 9 August 2009

“莫拉克”造成的台湾特大暴雨过程的对流云水汽输送结构亦可描述出上述双台风间水汽输送、交换及其相互影响结构。2009年8月7—9日“莫拉克”台风登陆造成了台湾岛历史罕见的暴雨灾害,为了分析“莫拉克”台风引发异常暴雨的水汽源及其水汽输送结构,图 7给出发生异常暴雨的台湾地区云顶黑体辐射亮温(TBB)和周边整层水汽通量的相关矢量场,可见在台湾区域异常暴雨前后时段(2009年8月4—11日)对流云水汽源主要来自“天鹅”至“莫拉克”双台风涡旋相互影响“连体”南侧的强水汽通量相关矢量通道以及偏南水汽通道。

图 7 台湾区域2009年8月4—11日云顶辐射亮温与周边区域整层水汽通量相关矢量场(实线为水汽相关通道) Fig. 7 Correlation vector field between the averaged satellite bright temperature over the Taiwan region and its peripheral vertically integrated moisture fluxes during the period from 4 to 11 August 2009(solid line denotes the moisture channel)
4 双台风“粒子群”输送三维轨迹及双台风相互影响模拟试验

为进一步揭示双台风相互影响过程水汽输送通道及其三维“粒子群”轨迹特征,采用Flexpart-WRF轨迹模式,计算分析了“天鹅”台风趋于消亡过程(2009年8月7—9日)涡旋质点“粒子群”输入“莫拉克”台风涡旋的动态过程,计算结果表明(图 8),该初始“粒子群”由“天鹅”台风水汽(前向轨迹)低层输入台风“莫拉克”,并描述出台风“天鹅”水汽“粒子群”向台风“莫拉克”输送的“连体”轨迹通道,这表明“天鹅”台风趋于消亡阶段其水汽、能量向“莫拉克”台风输送的可能轨迹,此外,试验研究揭示了“天鹅”水汽“粒子群”向“莫拉克”涡旋低层气旋式输入通道,在“莫拉克”涡旋高层则反气旋式卷出的三维立体动态物理图像。

图 8 Flexpart-WRF模拟的2009年8月7日23时(a),8日10时(b),8日21时(c),9日00时(d)水汽三维立体轨迹 Fig. 8 Flexpart-WRF-simulated 3D-trajectories at 23 :00 UTC 7(a)10:00 UTC 8(b)21:00 UTC 8(c),and 00:00 UTC 9(d)August 2009

Xu等(2011)利用剔除“天鹅”台风的敏感试验方案,验证了“天鹅”台风对“莫拉克”造成台湾地区异常暴雨的重要贡献。但双台风互旋、生消与吸引复杂过程的物理机制尚需进一步剖析,本研究重点从双台风动力、热力“三维”结构特征的演变及其相互影响的视角,探讨了“天鹅”台风消亡过程“莫拉克”台风发展并登陆引发异常暴雨的双台风相互作用。针对“天鹅”消亡,“莫拉克”发展维持过程,采用WRF-TC台风Bogus模型中的剔除台风模块(Low-Nam et al,2001)进行敏感试验,以分析“天鹅”台风的能量与水汽输送对于“莫拉克”台风增强、涡旋发展的影响。

图 9为控制试验与“莫拉克”和“天鹅”台风观测路径,在初始的18小时,模式较好地模拟出“莫拉克”的路径,在第18—24小时,其路径与实况有较大偏差,导致登陆台湾岛过程路径偏南。台风进入台湾海峡后与观测比较一致,但移速稍快于观测。“天鹅”台风移动路径与观测相比,很好地模拟出了其打转路径,但台风移动速度比观测慢。另外,对比观测与控制试验的降水量(图 10),结果表明观测与控制试验最大降水都超过1200 mm,且强降水位于台湾南部,模拟分布与观测分布较一致,但是,控制试验的最强降水区域比观测偏南,这是由于“莫拉克”登陆过程模拟路径偏南引起的降水位置偏移。总体上,控制试验能够较好地模拟出此次异常暴雨过程,且模拟出由于双台风相互作用而引起的“天鹅”打转路径,因此,本研究利用此次模拟结果结合剔除“天鹅”的敏感性试验研究“天鹅”台风对“莫拉克”的影响是合适的。

图 9 2009年8月7日00时—9日00时每6 h的莫拉克和天鹅台风最佳路径(黑色)与控制试验模拟路径(蓝色)的对比 Fig. 9 Best tracks of the typhoons(Morakot and Goni)(black) and the simulated tracks(blue)from 00:00 UTC 7 to 00:00 UTC 9 August 2009
图 10 2009年8月7日18时—8日12时累积降水

(a.观测,b. 控制试验)
Fig. 10 Accumulated rainfall from the observation(a) and the control experiment(b)during 18:00 UTC 7 to 12:00 UTC 8 August 2009

图 11为850 hPa控制试验与敏感性试验风场对比,控制试验中能够清晰地发现两个台风涡旋“天鹅”和“莫拉克”,而在敏感性试验中“天鹅”已经消失,仅剩下“莫拉克”涡旋。对比剔除“天鹅”台风涡旋的敏感试验与初始双台风涡旋的控制试验的模拟结果可发现,剔除“天鹅”台风后,“莫拉克”台风动能平均减少达15%(图 12a),平均水汽辐合量减少甚至超过23%(图 12b),说明“天鹅”台风涡旋的存在对于“莫拉克”台风发展与维持过程中动能、水汽的增加有显著的贡献。

图 11 2009年8月7日00时850 hPa风场(a.控制试验,b.敏感性试验) Fig. 11 850 wind field at 00:00UTC 7 August 2009 from the control(a) and the sensitivity experiment(b)
图 12 控制试验与敏感性试验“莫拉克”台风中心区域(以台风为中心相距4°的方形区域)内700 hPa总动能(a)和850 hPa总水汽散度(b)的变化 Fig. 12 Temporal variations of the simulated 700 hPa energy(a) and 800 hPa moisture flux divergence(b)within the Morakot center region(4°×4°)

对比图 13控制试验和敏感性试验模拟的“莫拉克”台风区域(图中方框所示)内风场、整层水汽输送通量与涡度等物理量可知,图 13a、c、e方框中“莫拉克”台风涡旋区风场、整层水汽与涡度高值范围与强度均显著大于图 13b、d、f方框中的相应物理量。另外,“天鹅”的存在对于“莫拉克”垂直结构亦有较大贡献(图 14),剔除“天鹅”试验与控制试验相比,右侧风速轴显著减弱,最大风速从45 m/s降低至27 m/s,相对而言,左风速轴变化较小。此模拟结果进一步印证了“天鹅”台风与“莫拉克”台风能量、涡度以及水汽输送相互作用及其影响效应,且“天鹅”在消亡阶段,对于“莫拉克”结构亦有显著的影响。

图 13 2009年8月8日00时控制试验(a、c、e)与敏感性试验(b、d、f)对比(a,b. 500 hPa 风场(矢量)和动能(阴影,m2/s2),c、d. 整层水汽通量(矢量)和通量值(阴影,kg/(m·s),e、f. 500 hPa涡度(阴影,10-5s-1)) Fig. 13 Comparison between control(a,c,e) and sensitive(b,d,f)experiments in terms of 500 hPa velocity(vector) and kinetic energy(shaded,m2/s2)(a,b),column vertically integrated moisture flux(vectors)with their magnitudes(color-shaded,kg/(m·s))(c,d) and vorticity(color-shaded,10-5s-1)(e,f)at 00:00 UTC 8 August 2009
图 14 模式模拟2009年8月7日18时“莫拉克”台风中心经向水平风速垂直剖面

(a. 控制试验,b. 剔除“天鹅”台风的敏感试验)
Fig. 14 Longitude-vertical cross sections the horizontal velocity along the Morakot typhoon center at 18:00 UTC 7 August 2009

(a. the control experiment,b. the sensitive experiment)

通过以上资料合成分析、数值模拟以及拉格朗日轨迹计算的综合分析,可以给出双台风生、消过程中能量、水汽输送的三维物理图像(图 15),从图 15可以看出,“莫拉克”(B)与“天鹅”(A)的台风涡旋两者互旋,相互吸引,并产生能量、水汽相互转换的情况。两台风在互旋过程中距离缩小,在环境场水汽的作用下,“天鹅”台风与“莫拉克”台风涡旋间出现能量、水汽输送的“连体”通道。“天鹅”台风通过此“连体”通道向“莫拉克”输送水汽、能量,“天鹅”台风粒子群三维动态轨迹亦表明“天鹅”台风水汽流中低层气旋式卷入“莫拉克”台风,且“莫拉克”向“天鹅”反馈水汽流相对较弱。两个台风呈现出反向变化趋势,即“天鹅”涡旋减弱,涡旋风场结构趋于非对称,而“莫拉克”涡旋增强,且维持显著准轴对称结构,双台风结构与强度反向变化特征亦印证了上述双台风相互作用特征。

图 15 双台风(A和B)互旋与相互作用的物理图像

(黑色实线箭头表示高低层流线,白色箭头为水汽输送,虚线箭头表示两者互旋方向)
Fig. 15 Schematic diagram for the mutual rotation and interaction of BTVs(A and B)

(The black solid arrows denote airflow; the white arrows denote the moisture flux channel; the dashed arrows illustrate the cyclonic rotation)
5 结 论

双台风涡旋的相互作用问题是台风理论研究的难点之一(陈联寿等,1979),本研究通过诊断结合数值模拟以及轨迹模拟对“莫拉克”与“天鹅”双台风相互作用进行研究,结论归纳如下:

(1)针对“莫拉克”与“天鹅”分别处于发展与消亡过程,综合分析表明双台风涡旋仍存在互旋,并相互吸引。从双台风涡旋三维结构变化的视角揭示出双台风生消过程其两个涡旋垂直结构强度呈反向演变特征。

(2)位涡、动能合成场分析可揭示出“天鹅”台风趋于消亡过程其与“莫拉克”台风涡旋间位涡、动能的“连体”通道。

(3)针对“莫拉克”台风引发特大暴雨过程,台湾区域云顶辐射亮温与周边整层水汽输送通量相关矢量场亦可描述出“天鹅”与“莫拉克”双台风整层水汽存在显著的“输送通道”结构,该双台风水汽输送通道亦是“莫拉克”台风引发特大暴雨的重要水汽源。

(4)采用Flexpart-WRF模式轨迹模型模拟了“天鹅”台风涡旋消弱过程其涡旋“粒子群”输入“莫拉克”台风涡旋过程的三维立体动态特征,并发现该阶段初始粒子群由“天鹅”台风中低层气旋式卷入“莫拉克”台风涡旋,且随着时间延长,粒子群从“莫拉克”台风涡旋高层顺时针卷出的轨迹特征。

(5)剔除天鹅台风涡旋的敏感性模拟试验结果进一步印证了双台风生消过程“天鹅”台风涡旋对“莫拉克”台风动能、水汽输送的重要贡献,揭示了双台风生成、消亡过程中涡旋间能量、水汽输送及其相互影响的机制。

(6)综合以上合成分析、轨迹模拟和数值模拟结论,提出了双台风生、消过程中能量、水汽输送相互影响的三维物理图像。

参考文献
陈联寿, 丁一汇. 1979. 西太平洋台风概论. 北京: 科学出版社,491pp
陈联寿, 孟智勇. 2001. 我国热带气旋研究十年进展. 大气科学, 25(3): 420-431
陈联寿, 徐祥德, 罗哲贤等. 2002. 热带气旋动力学引论. 北京: 气象出版社,317pp
罗哲贤, 马镜娴. 2001. 副热带高压南侧双台风相互作用的数值研究. 气象学报, 59(4): 450-458
罗哲贤, 徐祥德, 陈联寿. 2002. 涡旋合并过程的数值研究. 大气科学, 26(6): 807-816
台湾中央气象局台风数据库. 2009. Typhoon Morakot. http://rdc28.cwb.gov.tw/data.php?num=2009080804&year=2009&c_name=%B2%F6%A9%D4%A7J&e_name=MORAKOT
Brand S. 1970. Interaction of binary tropical cyclones of the western North Pacific Ocean. J Appl Meteor, 9(3): 433-441
Chanson H. 2010. The impact of Typhoon Morakot on the southern Taiwan coast. Shore Beach, 78(2): 33-37
Dudhia J. 1989. Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model. J Atmos Sci, 46(20): 3077-3107
Environment Protection Administration (Taiwan). 2010. Extreme events and disasters are the biggest threat to Taiwan-Typhoon Morakot. http://www.epa.gov.tw (Accessed 15 August 2010)
Fujiwhara S. 1921. The natural tendency towards symmetry of motion and its application as a principle in meteorology. Quart J Roy Meteor Soc, 47(200): 287-293
Fujiwhara S. 1923. On the growth and decay of vertical systems. Quart J Roy Meteor Soc, 49(206): 75-104
Kain J S. 2004. The Kain-Fritsch convective parameterization: An update. J Appl Meteor, 43(1): 170-181
Lin Y L, Farley R D, Orville H D. 1983. Bulk parameterization of the snow field in a cloud model. J Climate Appl Meteor, 22(6): 1065-1092
Low-Nam S, Davis C. 2001. Development of a Tropical Cyclone Bogussing Scheme for the MM5 System. Preprints of the Eleventh PSU/NCAR Mesoscale Model Users Workshop. June 25-27, Boulder, Colorado: 130-134
Mlawer E J, Taubman S J, Brown P D, et al. 1997. Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. J Geophys Res, 102(D14): 16663-16682
Prieto R, McNoldy B D, Fulton S R, et al. 2003. A classification of binary tropical cyclone-like vortex interactions. Mon Wea Rev, 131(11): 2656-2666
Stohl A, Forster C, Frank A, et al. 2005. Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2. Atmos Chem Phys, 5(9): 2461-2474
Wu C C, Huang T S, Huang W P, et al. 2003. A new look at the binary interaction: Potential vorticity diagnosis of the unusual southward movement of tropical storm Bopha (2000) and its interaction with supertyphoon Saomai (2000). Mon Wea Rev, 131(7): 1289-1300
Wu X, Fei J F, Huang X G, et al. 2012. A numerical study of the interaction between two simultaneous storms: Goni and Morakot in September 2009. Adv Atmos Sci, 29(3): 561-574, doi:10.1007/s00376-011-1014-7
Xu X D, Lu C, Xu H X, et al. 2011. A possible mechanism responsible for exceptional rainfall over Taiwan from Typhoon Morakot. Atmos Sci Lett, 12(3): 294-299
Yang C C, Wu C C, Chou K H, et al. 2008. Binary interaction between Typhoons Fengshen (2002) and Fungwong (2002) based on the potential vorticity diagnosis. Mon Wea Rev, 136(12): 4593-4611