“碧利斯”暴雨增幅高分辨率数值模拟及诊断分析

王黎娟，任晨平，3，崔晓鹏，冉令坤

(1．气象灾害教育部重点实验室(南京信息工程大学)，江苏南京210044;2．南京信息工程大学大气科学学院，江苏南京210044;3．中国人民解放军94754部队气象台，浙江嘉兴314013;4．中国科学院大气物理研究所云降水物理与强风暴实验室，北京100029)

0 引言

我国东邻西北太平洋，是世界上热带气旋(台风)登陆最多、灾害最严重的国家之一(陈联寿和孟智勇，2001)。台风是自然界中最强烈的造雨系统，台风Gloria(6312号)给台湾百新带来的最大24 h雨量达1 248 mm;台风Carla(6718号)的24 h降水达1 672 mm;台风Nina(7503号)在河南林庄的24 h累积降水达到1 062 mm，为中国大陆特大暴雨之最(陶诗言，1980)。

统计分析(董美莹等，2009)显示，我国引发特大暴雨的热带气旋年频数为1.8个，占我国登陆热带气旋总数的20%，发生暴雨突然增幅的热带气旋年频数为1.9个，占登陆总数的21%。针对登陆台风暴雨，科学家已经开展了大量的研究工作，并取得了很多有意义的研究成果。Atallah et al．(2007)应用准地转理论和位涡理论对降水明显分布在路径左侧和路径右侧的两类台风开展了对比分析，指出由于西风槽携带的冷空气侵入导致台风变性，其引发的降水分布在路径左侧，而降水发生在路径右侧则是台风环流和下游脊相互作用的结果。垂直风切变是制约台风发展的重要条件之一，同时也是影响对流及降水分布的重要因素之一。石顺吉等(2009)研究发现，环境风垂直切变在“碧利斯”(0604)一波非对称性降水中起到了重要作用。Rogers et al．(2003)对Bonnie(1998)飓风的数值研究表明，垂直风切变的显著变化对降水分布有重要影响，热带气旋内最强对流发生在强切变下风方向左侧。数值模拟发现，当热带气旋位于两环副高之间时，可出现暴雨增幅，在带状副高南侧时则很少出现增幅(董美莹等，2009)。孙建华等(2006)在研究Herb(9608号)台风造成北方特大暴雨过程时发现，强降水主要由两个中尺度对流系统(mesoscale convective system，MCS)造成;下垫面对台风暴雨增幅有着重要作用，陆地上的湖泊、大型水库、江河湖海可以向登陆热带气旋输送水汽和能量，有利于降水的维持，不同的干湿下垫面对台风暴雨的水汽输送、能量反馈有不同作用(董美莹等，2009)。郑庆林和吴军(1996)通过数值模拟研究发现，陡峭地形抬升辐合作用是暴雨增幅的重要因子，山地地形对大气的拖曳效应往往使得热带气旋降水出现明显增幅。Lin et al．(2002)研究认为，地形在“碧利斯”(0604)引发暴雨的过程中起到了重要作用。蔡则怡和宇如聪(1997)利用 LASG η坐标有限区域数值模式对“75.8”大暴雨开展了模拟，指出山脉相对高度和陡峭程度，以及山脉与暴雨系统的相对位置对暴雨的强度十分敏感，地形稍有改变，暴雨强度将减弱3/4以上。2006年第4号强热带风暴“碧利斯”登陆后，与西南季风相互作用(顾春利等，2010;Wang et al．，2010)，得到充足的水汽供应，气旋性环流在陆地上维持较长时间，并引发暴雨增幅，造成湖南、广东、江西等地区发生严重洪涝灾害，即季风水汽输送的突然增强对暴雨增幅有重要贡献。

由以上简要回顾不难发现，针对登陆台风导致的强降水，包括暴雨增幅，科学家已经开展了较多研究工作，取得了很多有意义的研究成果，但是由于台风暴雨过程不仅与登陆台风本身有关，还与台风同周边环流系统的相互作用过程以及下垫面的特性等有重要关系，其物理机制十分复杂，目前的研究进展还远远不能满足实际业务需求。台风暴雨，尤其是暴雨突然增幅，仍然是预报难点。2006年第4号强热带风暴“碧利斯”在登陆后出现明显的暴雨增幅过程，为这方面的研究提供了较典型的个例。虽然前期基于该个例已经开展了大量的研究工作(康志明等，2008;毛丽娜和潘益农，2009;石顺吉等，2009;Gao et al．，2009;郭英莲和徐海明，2010;Wang et al．，2010;罗玲等，2012;闫菲等，2012)，但是针对其登陆和登陆后的暴雨增幅过程开展高分辨数值模拟和动力诊断的研究还很少。本文首先利用美国NCEP/NCAR 1°×1°全球再分析资料、中国气象局台风年鉴资料以及地面观测资料等，借助ARPS(Advanced Regional Prediction System)中尺度非静力数值模式，对2006年第4号强热带风暴“碧利斯”的登陆过程以及登陆后的暴雨增幅过程开展观测分析、高分辨率数值模拟和模拟结果验证;进一步利用高分辨率模拟资料，选取散度垂直通量(周冠博等，2012)和湿位涡(吴国雄等，1995)两个动力因子对“碧利斯”引发的强降水进行初步的动力诊断分析;最后给出主要结论。

1 强热带风暴“碧利斯”概况和环流背景

2006年第4号强热带风暴“碧利斯”是该年度给我国造成严重灾害的热带气旋之一。“碧利斯”于7月8日下午在雅蒲岛东北约300 km处的西北太平洋洋面上生成，生成后向西北方向移动;9日下午加强为热带风暴，11日凌晨增强为强热带风暴，并逐渐向台湾省东北部一带沿海靠近，13日22时20分(世界时，下同)在我国台湾省宜兰附近沿海登陆，登陆时中心附近最大风力11级(30 m/s)。登陆后，“碧利斯”继续向西移动并穿过台湾省北部;随后进入台湾海峡并继续向西北方向移动;“碧利斯”于14日12时50分在我国福建省霞浦附近沿海再次登陆，登陆时中心附近最大风力11级(30 m/s)，中心气压975 hPa。二次登陆后“碧利斯”继续向偏西方向移动，强度减弱为热带风暴，15日凌晨进入江西省，下午风暴中心减弱为热带低压，并向西偏南方向移动;17日晚在越南北部地区减弱消失。“碧利斯”最强仅达到强热带风暴等级，但与之伴随的降水却十分集中，强度大，持续时间长，灾情严重。受“碧利斯”影响，13—17日福建、浙江、江西、湖南、广东、广西六省普遍出现大暴雨和特大暴雨(图1)，引发洪涝、山体滑坡、泥石流等自然灾害，给人民生命财产带来极大损失，据统计，受灾人员2 962.2万人，因灾死亡612人，失踪208人，直接经济损失高达266亿元人民币(康志明等，2008)。

1.1 “碧利斯”降水特点

“碧利斯”强降水主要出现在7月13—17日(图1)，浙江、福建、湖南南部、广东和云南大部分地区先后出现暴雨和特大暴雨，部分地区累计降水量达300～500 mm，其深入内陆时间之长、降水强度之大、影响范围之广，在历史上极为少见(卢珊等，2012)。登陆后强降水主要分布在“碧利斯”路径左侧(图1)，降水分布极不均匀，存在若干降水极值中心(图1、2)。7月14日12—18时，6 h累积降水小于25 mm(图2b)，14日18时—15日06时的逐6 h累积降水量均超过95 mm，12 h累积降水量达到特大暴雨(＞120 mm)级别(图 2c、d)，即在广东、湖南、江西三省交界处出现急剧暴雨增幅。“碧利斯”降水时空分布与众不同并出现明显暴雨增幅，造成了严重的灾害，应对其成因进行深入研究。

1.2 “碧利斯”环流特征

图3给出了2006年7月14日00时(暴雨增幅发生前)和18时(暴雨增幅开始时刻)的高、低空环流场。可见，在“碧利斯”整个登陆期间，东亚中高纬度地区500 hPa上两槽一脊环流形势稳定少变(图3)，副热带高压加强西伸，与大陆高压(脊)打通，在中纬度地区形成一道高压坝，导致“碧利斯”难以向北移动，登陆后在周围高压环流的包围下继续西进，且移动速度减缓，有利于其在陆上继续维持(康志明等，2008)。

图1 2006年7月13日00时—17日00时(世界时，下同)“碧利斯”过程雨量(单位:mm;图中黑点表示台风中心位置，时间间隔为6 h，数字标明当前时刻)Fig．1 Accumulated rainfall associated with Bilis from 0000 UTC 13 to 0000 UTC 17 July 2006(units:mm;the black dots in the figure indicate TC centers during the period with 6 h interval，and the digits mark the current time)

与此同时，低纬地区从阿拉伯海、孟加拉湾、中南半岛、南海到我国华南地区为一支宽广的水汽输送通量大值带(图4)，随着季风涌的北传源源不断地向风暴提供充足水汽(卢珊等，2012)。值得注意的是，暴雨增幅发生前，中南半岛附近地区的水汽输送强度相对较弱，水汽输送通道存在明显的断裂带(图3a、4a)，而暴雨增幅发生时刻，随着 80 ～90°E附近越赤道气流的增强，西南季风水汽输送明显加强(图3b、4b)，原来断裂的水汽输送通道完全打通，携带大量水汽进入已深入大陆内部的“碧利斯”环流中。这对“碧利斯”暴雨增幅的发生有重要贡献。

图2 “碧利斯”登陆期间6 h累计降水分布(单位:mm) a．14日06—12时;b．14日12—18时;c．14日18时—15日00时;d．15日00—06时;e．15日06—12时;f．15日12—18时Fig．2 The 6 h accumulated rainfall associated with Bilis in July 2006(units:mm) a．0600—1200 UTC 14 July;b．1200—1800 UTC 14 July;c．1800 UTC 14—0000 UTC 15 July;d．0000—0600 UTC 15 July;e．0600—1200 UTC 15 July;f．1200—1800 UTC 15 July

图3 2006年7月14日00时(a)和18时(b)高低空环流场(等值线为500 hPa位势高度，单位:dagpm;风向杆表示200 hPa风速大于30 m/s的风速矢;矢量箭头表示850 hPa风速大于12 m/s的风速矢)Fig．3 The circulation fields at(a)0000 UTC and(b)1800 UTC 14 July 2006(contour:500 hPa geopotential height with the unit of dagpm;barb:wind vector with speed larger than 30 m/s at 200 hPa;arrow:wind vector with speed larger than 12 m/s at 850 hPa)

2 数值模式与模拟方案设计

2.1 数值模式简介

ARPS(Advanced Regional Prediction System)是美国Oklahoma大学风暴分析和预报中心研发的高分辨率、非地转、非静力平衡的三维有限区域中尺度数值模式。模式采用全弹性非静力平衡方程组，包括动量、热量、质量、水物质和能量等预报方程和状态方程。模式垂直坐标采用广义高度地形追随坐标，水平为Arakawa-C交错网格，时间积分方案采用二阶蛙跳方案，主要物理过程包括大气辐射方案、可分辨尺度云微物理过程、次网格尺度积云对流参数化方案等。该数值模式可以模拟再现降水云系的微物理结构特征，可用于降水云系的微物理过程和降水机制的研究。

图42006年7月14日00时(a)和18时(b)850 hPa水汽通量(单位:g/(s·hPa·cm);阴影区表示大于等于 15 g/(s·hPa·cm))Fig．4 Water vapor flux vector(units:g/(s·hPa·cm);shadings:≥15 g/(s·hPa·cm))at 850 hPa at(a)0000 UTC and(b)1800 UTC 14 July 2006

2.2 模拟方案设计

利用ARPS模式，采用双层嵌套，粗、细网格水平分辨率分别为27和3 km。由于在本文前面分析(图3、4)中以及他人研究(卢珊等，2012)中已经发现，越赤道气流和西南季风涌爆发造成的水汽输送对“碧利斯”强降水有重要影响，因此粗网格模拟范围不仅覆盖了台风影响范围，同时也包括了西至索马里越赤道气流附近地区以及赤道以南的部分地区，以充分描述南半球和西南季风等对此次过程的影响。积云对流参数化方案选取Kain-Fritch方案，云微物理方案选择Lin-Tao冰相微物理方案，模式利用NCEP/NCAR 1°×1°全球再分析资料形成初始场和侧边界条件，从2006年7月14日00时开始运行，到2006年7月16日00时结束，共计48 h;细网格区域未使用积云对流参数化方案，模拟时段为2006年7月14日12时—15日12时，共计24 h，重点关注“碧利斯”暴雨增幅时段和区域(图5)。

图5 模拟区域设置Fig．5 The area of numerical simulation

3 模拟结果验证

3.1 路径及强度对比

模拟路径与日本气象厅区域气象中心(the RegionalSpecialized Meteorological Center，简 称RSMC)发布的最佳路径相比(图6)可知，模式较好地模拟再现了强热带风暴“碧利斯”的移动路径，登陆后模拟与观测的“碧利斯”均向偏西方向移动，在细节上存在一定差别。模拟时间开始于7月14日00时，模拟的登陆地点比实况略偏南0.5个纬度左右，这是由于模拟初始时刻，初始场中台风中心的位置与最优观测之间已经存在一定的误差。登陆后，模拟的“碧利斯”先向西北方向移动，14日18时起由西北转向西南，12 h后再次转向西北移动。6 h后，转为向偏西方向移动，与观测略有差别。但从总体上看，模式还是较好地再现了“碧利斯”的移动路径。

图6 “碧利斯”观测路径与模拟路径对比Fig．6 Observed and simulated tracks of Bilis

强度方面，模拟时段内，模拟的“碧利斯”中心最低海平面气压的变化趋势与RSMC最优观测资料吻合度较高(图7)，由于初始场存在误差，两者在模拟起始阶段存在较大差距。随着时间的推移，两者演变趋势逐渐趋于一致。15日00时以后，模拟的台风强度有所加强，而观测强度此时也停止减弱，这可能是由于西南季风涌的明显北推，输送大量水汽进入“碧利斯”环流，导致“碧利斯”强度得以维持;此后，“碧利斯”强度进一步缓慢减弱，观测与模拟具有较好一致性。

图7 观测与模拟的“碧利斯”中心最低气压对比Fig．7 Observed and simulated minimum sea surface pressure of Bilis

3.2 高低空环流形势对比

图8为观测与模拟的高低空环流形势。可见，模拟的高低空环流形势及其演变与实况吻合得非常好，模式较好地模拟再现了中高纬度两槽一脊的环流形势及其维持特征(图8a—j)、副高的形态特征和演变及其与陆地高压(脊)打通形成中纬度高压坝的特征、以及低纬西南季风水汽输送(图8a—j)、季风输送明显加强(图8c、d)和越赤道气流明显加强过程等(图8c、d)。这些关键因子的成功模拟再现，保证了模式能够较好地模拟再现“碧利斯”的移动路径、强度演变特征以及强降水分布等。

3.3 降水对比

图9为7月14日18时—15日06时的逐6 h累积降水对比。可见，观测降水主要分布在湖南南部、广东东部、江西和福建南部，其中湖南南部最强降水强度达96 mm/(6 h)。模式模拟出了这一强降水中心，位于广东、湖南和江西三省交界处，较实况略偏东南，强度也略偏强。广东东南部较弱的雨带也得到较好的再现，极值中心略向东北侧偏移(图9a、b)。15日 00—06时(图 9c、d)，强降水主要位于广东中南部、江西东南部和湖南南部，其中湖南南部和广东中南部的降水最强。此时段降水的模拟效果比前6 h略差，模拟的湖南南部和江西东南部的降水中心与观测相比位置偏南，没有模拟出位于广东中南部的降水中心。尽管与观测降水分布和强度均存在一定差距，但考虑到模拟过程的复杂性和时空尺度较小等因素，可以说该模拟结果是比较令人满意的。

综合上述模拟验证分析可知，模式较好地模拟再现了强热带风暴“碧利斯”在我国大陆的登陆过程以及登陆之后的强降水过程，发生在湖南、江西和广东三省交界附近地区的暴雨增幅过程的模拟也取得了令人满意的结果。下面将利用高分辨率的模拟资料，结合动力诊断量，对暴雨增幅过程的强降水落区开展诊断分析。

4 暴雨增幅诊断分析

4.1 散度垂直通量诊断分析

暴雨系统的发生、发展与垂直运动以及辐合辐散有着密切的联系，因此本文首先选取了散度垂直通量对“碧利斯”的暴雨增幅过程开展诊断分析。冉令坤和楚艳丽(2009)在研究华北暴雨落区和雨带移动方向时，为了更加准确地描述暴雨系统的发展和演变，将水平散度与垂直运动这两个与暴雨发生直接相关的关键因子结合起来，构造出一个新的动力诊断量—散度垂直通量:

图 8 观测(a，c，e，g，i)和模拟(b，d，f，h，j)的高低空环流形势对比(等值线表示 500 hPa 位势高度，单位:dagpm;风向杆表示250 hPa大于25 m·s－1的风场;矢量箭头表示850 hPa大于12 m·s－1的风场) a，b．7 月14 日00 时;c，d．7 月14 日18时;e，f．7月15 日00时;g，h．7月15 日06 时;i，j．7月16 日00时Fig．8 (a，c，e，g，i)Observed and(b，d，f，h，j)simulated circulations(contour:500 hPa geopotential height with the unit of dagpm;barb:wind vector with speed larger than 25 m/s at 250 hPa;arrow:wind vector with speed larger than 12 m/s at 850 hPa) a，b．0000 UTC 14 July;c，d．1800 UTC 14 July;e，f．0000 UTC 15 July;g，h．0600 UTC 15 July;i，j．0000 UTC 16 July

对Q的绝对值取850～500 hPa的质量权重垂直积分，得到:

利用上述ARPS模式模拟得到的“碧利斯”暴雨增幅过程高分辨率资料，计算了暴雨增幅期间6 h平均的散度垂直通量(图10)。7月14日18时—15日00时，降水主要分布于江西和湖南南部、广东中北部及福建西南和东北部，25 mm以上的暴雨主要发生在广东、湖南、江西三省交界附近，广东西部和东部以及福建东北部地区降水极大值达到130 mm以上。Qsum对强降水有着非常明显的指示意义，其分布与暴雨大值区十分吻合(图10a)。图10b为沿着6 h累积降水极大值中心所做的纬向垂直剖面，降水极大值区对应散度垂直通量的大值区，正值区从低层一直延伸到500 hPa附近，表明在500 hPa以下的对流层低层存在很强的辐合上升运动，最强中心在750 hPa附近;负值区从500 hPa一直到200 hPa以上，表明在对流层中高层以辐散为主，极值中心位于300 hPa附近。15日00—06时(图10c、d)，降水主要落区稍向东南移动，但仍主要分布在广东、湖南和江西三省交界附近地区(图10c)，Qsum的分布与强降水落区仍保持非常好的一致性。高度—经度垂直剖面(图10d)显示，Q的正值区仍主要位于500 hPa以下区域，维持低层辐合抬升、高层辐散下沉的垂直动力结构，因此在此时段内强降水得以持续。

图9 模拟(a，c)与观测(b，d)的6 h累积降水量(单位:mm) a，b．2006年7月14日18时—15日00时;c，d．2006年7月15日00—06时Fig．9 (a，c)Simulated and(b，d)observed 6 h accumulated rainfall(units:mm)a，b．from 1800 UTC 14 to 0000 UTC 15 July 2006;c，d．from 0000 UTC to 0600 UTC 15 July 2006

4.2 湿位涡诊断分析

位势涡度(potential vorticity，PV)综合反映了大气的热力和动力性质，有着广泛的应用(寿绍文，2010)。发生降水特别是强降水(暴雨)时，必须考虑水汽的作用，从而引入了湿位涡的概念。对于湿空气，用相当位温θe代替位温θ，可得到湿位涡的表达式(吴国雄等，1995):

湿位涡不仅表征了大气的动力和热力性质，还考虑了水汽的影响，运用湿位涡开展降水诊断，可以探索各个动力、热力和水汽条件与降水之间的关系，进而揭示降水，特别是强降水发生、发展的内在物理机制(安洁和张立凤，2004;寿绍文，2010)。

在p坐标系下，忽略垂直速度ω的水平变化，方程(3)可以展开(留小强等，1994)为:

其中:ζp表示相对涡度(单位:s－1);f是科氏参数(单位:s－1);θe为相当位温(单位:K)。湿位涡的单位(PVU)为 10－6m2·s－1·K·kg－1。与散度垂直通量类似，对湿位涡的绝对值同样取850～500 hPa的质量权重垂直积分，得到:

图10 6 h累计降水和6 h平均散度垂直通量 a．2006年7月14日18时—15日00时6 h累计降水(阴影;单位:mm)和垂直积分散度垂直通量(等值线;单位:kg·Pa2·m－2·s－3);b．2006年7月14日18时—15日00沿图10a中6 h降水量最大值区的散度垂直通量气压—经度剖面(单位:10－4Pa·m·s－3;直方图代表6 h累计降水，单位:mm);c．2006年7月15日00—06时6 h累计降水(阴影;单位:mm)和垂直积分散度垂直通量(等值线;单位:kg·Pa2·m－2·s－3);d．2006年7月15日00—06时沿图10c中6 h降水量最大值区的散度垂直通量气压—经度剖面(单位:10－4Pa·m·s－3;直方图代表6 h累计降水，单位:mm)Fig．10 The 6 h accumulated rainfall and averaged Q a．the 6 h accumulated rainfall(shadings;units:mm)and averaged Qsum(contour;units:kg·Pa2·m －2·s－3)from 1800 UTC 14 to 0000 UTC 15 July 2006;b．the pressure-longitude vertical cross section of Q along the maximum rainfall in Fig．10a from 1800 UTC 14 to 0000 UTC 15 July 2006(units:10 －4Pa·m·s－3;bars indicate the 6 h accumulated precipitation);c．the 6 h accumulated rainfall(shadings;units:mm)and averaged Qsum(contour，units:kg·Pa2·m －2·s－3)from 0000 UTC to 0600 UTC 15 July 2006;d．the pressure-longitude vertical cross section of Q along the maximum rainfall in Fig．10c from 0000 UTC to 0600 UTC 15 July 2006(units:10－4Pa·m·s－3;bars indicate the 6 h accumulated precipitation)

利用上述“碧利斯”暴雨增幅过程的高分辨率模拟资料，计算850～500 hPa垂直积分的湿位涡，分析其分布与暴雨增幅强降水落区的关系。由图11可见，Msum与6 h累计强降水也具有较好的对应关系，强降水中心附近存在明显的Msum大值区，该大值区东西方向水平尺度约为150 km，南北方向约为100 km，中尺度特征十分显著。但强降水中心与Msum大值中心并不完全对应，Msum大值中心主要位于强降水中心边缘附近，由Msum的定义可知，这种对应可能与强降水中心周围湿等熵面较为倾斜(有较大的相当位温水平梯度)有关。值得注意的是，Msum在降水区东北侧还存在一个范围较大、强度较弱的闭合等值线区，该区域的存在可能与“碧利斯”本身的环流分布有关(该时刻“碧利斯”中心位于强降水中心的东北侧)，而湖南、广东和江西三省交界附近的强降水则是由“碧利斯”台风环流中的中尺度系统引发的。

图11 模拟的6 h累积地面降水(阴影;单位:mm)和6 h平均的850～500 hPa垂直积分的湿位涡(等值线;单位:10－4Pa·K·m－1·s－1)分布 a．2006年 7月14日18时—15日 00时;b．2006年7月15日00—06时Fig．11 Simulated 6 h accumulated rainfall(shadings;units:mm)and averaged Msum(contour;units:10 －4Pa·K·m －1·s－1) a．from 1800 UTC 14 to 0000 UTC 15 July 2006;b．from 0000 UTC to 0600 UTC 15 July 2006

5 结论

利用中尺度非静力数值模式ARPS对2006年第4号强热带风暴“碧利斯”的登陆过程及登陆后发生在广东、湖南和江西三省交界的暴雨增幅过程开展了高分辨率数值模拟，并利用观测资料进行了较为细致的模拟验证;随后，利用经过验证、认为较为可靠的高分辨率资料，针对暴雨增幅开展了动力诊断分析，主要结论如下:

1)模式较好地模拟再现了“碧利斯”的发展、演变和登陆过程，模拟的台风路径和强度变化、高低空环流特征和演变与观测均较为一致。尤其值得指出的是，“碧利斯”登陆我国大陆后引起的广东、江西、湖南三省交界处的暴雨增幅过程也取得了较为满意的模拟效果。

2)利用高分辨率模拟资料，结合散度垂直通量Q和湿位涡M对强降水开展的诊断分析表明，暴雨增幅期间，散度垂直通量Q和湿位涡M的分布与强降水均有良好的对应关系，对降水落区的指示意义较强。其中散度垂直通量与强降水落区的对应关系较好，湿位涡不仅很好地反映了与强降水(系统)相关的湿位涡分布，同时也很好地反映了“碧利斯”环流本身所对应的湿位涡分布特征。

安洁，张立凤．2004．暴雨过程中湿位涡的中尺度时空特征［J］．气象科学，24(1):72-80．

蔡则怡，宇如聪．1997．LASG η坐标有限区域数值预报模式与一次登陆台风特大暴雨的数值试验［J］．大气科学，21(4):459-471．

陈联寿，孟智勇．2001．我国热带气旋研究十年进展［J］．大气科学，25(3):420-432．

董美莹，陈联寿，郑沛群，等．2009．登陆热带气旋暴雨突然增幅和特大暴雨之研究进展［J］．热带气象学报，25(4):495-502．

顾春利，倪允琪，徐又平，等．2010．南海地区准单、双周震荡对“碧利斯”降水维持影响分析［J］．暴雨灾害，29(1):7-13．

郭英莲，徐海明．2010．对流层中上层干空气对“碧利斯”台风暴雨的影响［J］．大气科学学报，33(1):99-109．

康志明，陈涛，钱传海，等．2008.0604号强热带风暴“碧利斯”特大暴雨的诊断研究［J］．高原气象，27(3):596-607．

留小强，王田民，吴宝俊，等．1994．湿位涡方程及其应用［J］．大气科学，18(5):569-575．

卢珊，王黎娟，管兆勇，等．2012．低纬季风涌影响登陆台风“榴莲”(0103)和“碧利斯”(0604)暴雨增幅的比较［J］．大气科学学报，35(2):175-185．

罗玲，娄小芬，胡亮，等．2012.0604号强热带风暴“碧利斯”结构特征及其暴雨成因分析［J］．气象科学，doi:10.3969/2012jms．0097．

毛丽娜，潘益农．2009．环境风垂直风切变对热带气旋“碧利斯”的影响［J］．气象科学，29(4):507-512．

冉令坤，楚艳丽．2009．强降水过程中垂直螺旋度和散度通量及其拓展形式的诊断分析［J］．物理学报，58(11):8094-8106．

石顺吉，余锦华，张大林．2009．热带风暴Bilis(2006)登陆期间一波非对称降水分布成因的探讨［J］．热带海洋学报，28(1):34-42．

寿绍文．2010．位涡理论及其应用［J］．气象，36(3):9-18．

孙建华，齐琳琳，赵思雄．2006．“9608”号台风登陆北上引发北方特大暴雨的中尺度对流系统研究［J］．气象学报，64(1):57-71．

陶诗言．1980．中国之暴雨［M］．北京:科学出版社:225．

吴国雄，蔡雅萍，唐晓箐．1995．湿位涡和倾斜涡度发展［J］．气象学报，53(4):387-405．

闫菲，李艳伟，周毓荃，等．2012．受“碧利斯”影响的一次中尺度对流系统模拟研究［J］．大气科学学报，35(6):737-745．

郑庆林，吴军．1996．地形对9216号台风暴雨增幅影响的数值研究［J］．南京气象学院学报，19(1):8-17．

周冠博，崔晓鹏，高守亭．2012．台风“凤凰”登陆过程的高分辨数值模拟及其降水的诊断分析［J］．大气科学，36(1):23-34．

Atallah E，Bosart L F，Aiyyer A R．2007．Precipitation distribution associated with landfalling tropical cyclones over the eastern United States［J］．Mon Wea Rev，135:2185-2205．

Gao S Z，Meng Z Y，Zhang F Q，et al．2009．Observational analysis of heavy rainfall mechanisms associated with severe tropical Storm Bilis(2006)after its landfall［J］．Mon Wea Rev，137:1881-1897．

Lin L L，Ensley D B，Chiao S，et al．2002．Orographic influences on rainfall and track deflection associated with the passage of a tropical cyclone［J］．Mon Wea Rev，130(12):2929-2950．

Rogers R，Chen S，Tenerelli J，et al．2003．A numerical study of the impact of vertical shear on the distribution of rainfall in Hurricane Bonnie(1998)［J］．Mon Wea Rev，131:1577-1599．

Wang L J，Lu S，Guan Z Y，et al．2010．Effects of low-latitude monsoon surge on the increase in downpour from tropical storm Bilis［J］．J Trop Meteor，16(2):101-108．