杨宁，张其林

(南京信息工程大学1.气象灾害省部共建教育部重点实验室;2.大气物理学院，江苏南京210044)

0 引言

台风是地球上最强烈的风暴，在北太平洋西部(国际日期线以西，100°E以东)称为台风，在北大西洋及东太平洋称之为飓风。一般始发于洋面温度超过26℃的热带或副热带海洋上，从热带低压发展成为热带风暴，进一步加强成为强热带风暴后最终形成台风。整个台风的生命周期可维持几天到几周不等，其大部分时间都位于洋面上，登陆后由于下垫面的变化，台风系统中的水汽供应量减少，能量无法得到补充，难以维持整个台风系统，致使其逐步削弱直至消亡。最近的统计表明最近的几十年中台风的总体强度呈现上升趋(Emanuel，2005;Webster et al.，2005)。

多年来，闪电活动一致被认为与对流活动的微物理和动力学过程有着密切的联系。当上升气流速度大于10 m/s时，就能促使电荷分离从而建立起足够强大的电场。有研究表明，当台风中上升气流速度突破20 m/s时，台风强度将快速增强，同时台风中闪电活动也会出现明显的增长(Black et al.，1994)。这表明台风强度与闪电频数之间存在着一定的联系。国外的一些学者使用LLDN(long range lightning detection network)资料分析了2005年发生在大西洋上的Rita与Katrina台风，结果表明，在台风强度迅速增强期间眼壁上有闪电爆发(Squires and Businger，2008)。随后Shao et al.(2005)利用LASA(Los Alamos Sferic Array)资料同样分析了Rita与Katrina台风，发现在台风增强期间闪电爆发，数量达到600次/h。在这之后Thomas et al.(2010)研究了Rita、Katrina、Emily台风，发现在台风进入减弱阶段之前，系统内部正地闪数量会增加。国内曾利用LIS(Lightning Imaging Sensor)资料分析了1998—2005年太平洋区域热带气旋中的闪电特征，结果表明在西北太平洋热带气旋的眼壁附近、内雨带和外雨带内普遍存在闪电活动(雷小途等，2009)。随后潘伦湘等(2010)研究了西北太平洋地区2005—2008年的7个台风的闪电空间分布及时间演变特征，发现西北太平洋地区强台风成熟时平均闪电密度呈现出三圈结构。

为了进一步了解台风过程中闪电活动特征，本文选取了2005—2010年西太平洋洋面共计55次台风过程，并对台风系统内的闪电活动及其与最大风速之间的关系进行了统计和分析。

1 资料及方法介绍

本文闪电定位数据由美国华盛顿大学地球与空间科学中心建设的WWLLN(World Wide Lightning Location Network)网络提供。该网络采用闪电发生时产生的3～30 kHz甚低频电磁辐射信号对闪电进行定位。由于甚低频信号频率比较低，波长较长，并可以通过电离层反射。这些特性使其可以传输至几千千米之外并被接收到，使得该网络可以远距离对闪电进行探测。WWLLN的优势在于可以实时不间断地进行探测，提供一个完整台风过程的闪电定位资料(Lay et al.，2007)。该探测网利用分布于全球范围内的各子站点，采用改进的时间到达法对闪电进行定位。每个站点将实时接收到的信号传输至网络中央数据处理中心进行处理、计算，最终得到闪电发生位置。WWLLN网络分布在全球的站点从2003年的11个，已经扩展到2010年的50个。

对2004年澳大利亚国家闪电定位网与WWLLN系统的对比后发现，在这期间WWLLN系统的总体探测效率达到了13%，平均探测精度为3.4 km(Rodger et al.，2005)。近期有学者使用WWLLN资料与美国国家闪电探测网资料作比较，发现WWLLN的总体探测效率达到10.30%，探测精度达到4.03～4.98 km。对于电流强度大于35 kA的闪电，WWLLN探测效率开始持续增长，对于电流强度大于55 kA的闪电，探测效率达到了20%，这说明WWLLN对强闪电有着良好的探测效率(Sergio et al.，2010)。虽然WWLLN系统的总体探测效率约为10%左右，但这并不影响其捕捉台风中各时刻闪电发生的总体趋势，其在大尺度范围内实时提供闪电数据的能力是其他系统无法比拟的。

本文台风资料来源于中国气象局上海台风所。台风强度的定义取自《热带气旋等级》国家标准(中华人民共和国国家质量监督检验检疫总局和中国国家标准化管理委员会，2006)。定义台风的最大风速大于32.7 m/s，强台风的最大风速大于41.5 m/s，超强台风的最大风速大于51 m/s。因西太平洋地区属于台风高发区域，所以从该区域中得出的结论可具有一定代表性。图1是本文选用的2005—2010年55次台风的移动路径及发生区域。该55次台风分别为:2005年的Damrcy、Talim、Sonca、Nesat、Haitang、Matsa、Mawar、Nabi、Khanun、Saola、Longwang、Kirogi;2006年的Chanchu、Ewiniar、Saomai、Shanshan、Yagi、Xangsane、Cimaron、Utor;2007年的Kong-rey、Yutu、Usagi、Man-yi、Sepat、Nari、wipha、Krosa、Kajiki;2008年的Rammasun、Nakri、Fengshen、Fung-wong、Hagupit、Jangmi、Dolphin;2009年的Kugira、Chan-hom、Molave、Morakot、Vamco、Choi-wan、Parma、Melor、Lupit、Nida;2010年的Conson、Chanthu、Dianmu、Kompasu、Malou、Fanapi、Malakas、Megi、Chaba。

图1 2005—2010西北太平洋地区55次台风过程分布及路径汇总Fig.1 The 55 typhoons'distribution and paths over northwest Pacific during 2005—2010

台风在其生命期间不断地移动，考虑到台风移动时速比较缓慢，最高时速约在25 km/h。故本文利用每个整点时刻至下一整点时刻，以台风中心点为原点的780 km半径范围内发生的闪电来代表该时刻的闪电数，并视台风在这一时段内做匀速运动。由于台风的水平尺度一般为1 000 km左右，因此将半径r选为780 km是比较合理的，基本可以将整个台风系统包含在内。台风路径数据中每隔6 h给定一次台风中心经纬度，考虑到6 h的时间步长会对精度产生影响，故采用线性插值方法，对原始台风定位数据进行插值处理，得出逐时台风中心位置。

2 台风过程中闪电时空分布特征

为了详细分析台风系统内闪电的普遍规律，本文对55次台风过程中的闪电活动进行了统计分析，发现台风过程中的闪电活动具有一定的相似性。因此，下面选取5次台风过程进行介绍，如表1所示。从表中可见，台风过程中伴随着明显的闪电活动，但由于个例之间的差异，各台风系统内的总闪数变化幅度较大，最高日平均闪电数达到1 842次/d，最低仅为342次/d。考虑WWLLN系统对地闪与云闪的探测效率分别为10.3%和4.82%(Sergio et al.，2010)，修正后与表2陆地各地区不同对流系统比较后发现，台风系统中闪电频数低于陆地对流系统，但系统内的闪电活动总体来说仍然比较频繁。

表1 2005—2010年间5次台风过程闪电活动特征Table 1 The characteristics of five typhoons'lightning activity during 2005—2010

表2 陆地各地区对流系统闪电频数Table 2 The lightning frequency of convection system on land

台风系统一般可根据径向距离划分为3个区域，分别为眼区、云墙区以及外层区(包括外云带和内云带)。图2给出了这5次台风过程中的闪电空间分布特征。可以看出，台风过程中闪电活动具有3个高密集区。首先，在位于台风中心20～40 km处的云墙区处存在着一个显著的闪电高发带，闪电频数随径向距离的增加迅速增大。其次，在距离中心280～340 km范围内出现第二个密集区，闪电频数较内雨带中有明显增加。随后，闪电数随径向距离振荡波动并在距中心440～580 km附近区域激增，出现第三个密集区。

图3进一步给出了5次台风过程中闪电活动的时变特征。可以看出，各台风系统中闪电数量随时间的变化并不连续，短时间内闪电的活动变得异常活跃，且在闪电频数到达峰值后闪电活动迅速减弱，呈现振荡变化的状态。由图3还可以发现，在台风强度增强(中心气压降低)阶段，系统内闪电数激增，闪电活动在这一阶段相当活跃呈现爆发状态。而在台风强度减弱阶段，系统内闪电活动较少，变化较为平稳。

图2 各台风系统中闪电随径向距离分布特征a.韦帕;b.海豚;c.环高;d.彩云;e.妮妲Fig.2 The lightning characteristics of typhoon systems with the radial distance distributiona.Wipha;b.Dolphin;c.Vamco;d.Choi-wan;e.Nida

图3 台风系统内闪电随时间变化特征(柱状图代表台风系统内每小时闪电数;虚线代表台风中心最低气压)a.韦帕;b.海豚;c.环高;d.彩云;e.妮妲Fig.3 The lightning characteristics of typhoon system with time variation(The bar chart denotes the number of lightning per hour in typhoon system;the dotted line denotes the minimum pressure in typhoon center)a.Wipha;b.Dolphin;c.Vamco;d.Choi-wan;e.Nida

台风强度增加时系统内的闪电数发生增长，可能是由于此时系统内变得更有效的起电机制导致的。通过飞机穿越大西洋台风所得到的观测资料发现，在台风中存在着过冷水滴、软雹和冰晶粒子(Black and Hallett，1986)。当过冷水滴凝结在软雹粒子上时，粒子质量增加并下坠，并与质量较轻的冰晶粒子碰撞，两者之间的碰撞摩擦使电荷发生转移。由于冰的热电效应，温度较高的雹粒子带负电荷，而温度较低的冰粒子带正电荷。再经过重力分离作用，带正电荷的冰粒子随上升气流升至上部，带负电荷的雹粒子因重力作用沉降至下部，在系统中形成正负电荷区。对于冰晶与雹粒子单离散性谱分布，荷电的速率为

其中:Ng、Ni分别为雹粒与冰晶的浓度;Dg、Di分别为这两类粒子的直径;Egi是碰撞反弹系数;ΔV=Vg－Vi是碰撞雹粒子和冰晶粒子的速度差。

当台风强度开始增强时，中心气压开始降低，促使气流不断向系统中心辐合并作上升运动。上升过程中水汽凝结释放出潜热补给台风能量，使得系统内气流的上升运动越来越强。当上升气流增强时，根据上述起电机制可知此时电荷分离过程将更加迅速。通过Black and Hallett(1999)对约20 a穿越台风的飞行观测资料总结后发现，闪电活动与速度大于10 m/s的上升气流存在着密切的联系。荷电速率的增加导致系统中电荷总量增加，这需要更多次数的放电使系统内的电荷被中和，在系统内就表现为闪电活动的增强。每次闪电发生后，系统中的电荷量便会相应减少，而此时系统内的起电机制会不断作用补充消耗的电荷，继续维持闪电的发生。

3 台风过程中闪电频数与最大风速的特征

为了研究台风中闪电频数与台风最大风速之间的相关性，本文选取了2005—2010年间发生在西太平洋地区的55次台风过程并进行了统计。下面以200906号台风莫拉菲Molave为例。由2009年7月15—19日间莫拉菲的移动路径(图4)可见，它始发于一个位于菲律宾以东海域的热带低气压，7月16日20时(北京时)，升级为热带风暴。得益于西太平洋高海温、强劲的辐散辐合及风切变的环境，其强度快速增强，于17日下午增强成强热带风暴。因减速及洋面温度较高的原因，台风强度进一步增强，在18日清晨升级成台风。在进入南海北部后，因得到西南季风带来的水汽支援，莫拉菲的强度继续缓慢增强，直至登陆前中心风力仍保持增长。19日凌晨，台风莫拉菲在深圳附近一带登陆，于05时(北京时)减弱为强热带风暴，随后急剧减弱，逐渐消散。

图5为整个莫拉菲台风过程中逐时闪电频数变化与最大风速的时变特征。从图中可以看出，闪电频数随着台风风速的增强而增加，闪电活动于18日07时(北京时)开始呈现活跃状态，随后系统内闪电爆发并于10时到达峰值，达到峰值后系统内闪电数目迅速减少。而在闪电数达到了峰值后系统内风速极值于12时左右出现。两者变化的趋势存在一致性。值得注意的是，闪电频数最大值出现时刻先于最大风速出现的时刻。取18日00—18时闪电数和风速两者的数据做相关性分析，得出两者的相关系数为0.89。

图6给出了55次台风过程中最大风速滞后于闪电峰值的时间分布曲线，横坐标表示时滞时间，纵坐标为出现该时滞时间的概率。由图6可见，最大风速滞后于闪电峰值的概率明显大于超前的概率，经过统计，83.6%的台风过程存在风速峰值滞后于闪电峰值这一现象，41.8%的滞后时间为0～24 h，滞后时间峰值约为6 h，平均滞后时间为26.6 h。

图4 2009年7月15—19日莫拉菲台风移动路径示意Fig.4 Typhoon Molave path during July 15 and 19 in 2009

图5 莫拉菲台风过程中逐时闪电数与近中心风速特征(柱状图表示台风系统内逐时的闪电数;虚线表示台风近中心风速)a.完整过程;b.局部放大Fig.5 The number of lightning and characteristics of wind speed in the center of Molave typhoon(The bar chart denotes the number of lightning in typhoon system per hour;the dotted line denotes the wind speed in typhoon center)a.the whole process;b.the detailed view

图6 2005—2010年台风最大风速与闪电峰值时滞分布(横坐标中正值表示闪电峰值提前于风速极值出现的天数，负值表示风速极值先于闪电峰值出现的天数;纵坐标为出现这一时滞的概率)Fig.6 Lag time between the maximum lightning activity and the maximum winds of typhoons during 2005—2010(Positive values imply that the lightning activity peaked before the maximum winds，and negative values imply the opposite;the vertical ordinate indicates the occurrence probability of lag time)

4 结论与讨论

本文利用全球闪电定位网(WWLLN)的闪电定位数据与中国气象局上海台风所提供的台风定位资料，对2005—2010年发生在西太平洋地区的55次台风过程中闪电活动与最大风速两者间的关系进行了研究。得到以下结论

1)西太平洋台风系统中发生的闪电大致有3个高发区，分别位于距离中心20～40 km的云墙区，280～340 km及440～580 km的外云带中。台风系统中内云带中的闪电活动发生较少。

2)在台风过程中，闪电往往在台风强度增强阶段呈现爆发状态。在此期间系统内闪电活动频繁，闪电数目在短时间内剧增，在达到峰值后急剧减少并回落至较少发生的状态。在台风强度减弱阶段，系统内闪电活动较少，变化较为平稳。登陆后，台风强度迅速减弱。系统内的闪电也随之相应减少。

3)大约83.6%的台风存在风速峰值滞后于闪电峰值这一现象。最大风速的出现相对于闪电数的爆发存在一定的滞后性，滞后时间峰值约为6 h，一般滞后时间在0～24 h内。

本文仅利用了WWLLN闪电定位资料和上海台风研究所提供的台风路径数据，对西太平洋上空台风过程中的闪电特征进行了初步分析和统计。由于WWLLN探测效率的局限，下一步工作的重点是将雷达、卫星等资料融合至台风闪电时空分布特征分析中，可能会得到更好的结果。随着对台风中闪电活动与系统内微物理过程两者间关系的进一步研究，还可将闪电活动作为参考因子加入到台风模型中，更好地对台风过程进行预测。

雷小途，张义军，马明.2009.西北太平洋热带气旋的闪电特征及其与强度关系的初步分析［J］.海洋学报，31(4):29-38.

潘伦湘，邱秀书，刘东霞，等.2010.西北太平洋地区强台风的闪电活动特征［J］.中国科学D辑，40(2):252-260.

郄秀书，张其林，冯桂力，等.2007.中尺度对流系统中雷电的人工引发及特征放电参量的测量和估算［J］.中国科学D辑，37(6):564-572.

中华人民共和国国家质量监督检验检疫总局，中国国家标准化管理委员会.2006.GB/T 19201-2006热带气旋等级［S］.北京:中国标准出版社.

Black R A，Hallett J.1986.Observations of the distribution of ice in hurricanes［J］.J Atmos Sci，43:802-822.

Black R A，Hallett J.1999.Electrification of the hurricane［J］.J Atmos Sci，56:2004-2028.

Black R A，Bluestein H B，Black M L.1994.Unusually strong vertical motions in a Caribbean hurricane［J］.Mon Wea Rev，122:2722-2739.

Buechler D E，Driscoll K T，Goodman S J，et al.2000.Lightning activity within a Tornadic thunderstorm observed by the optical transient detector(OTD)［J］.Geophys Res Lett，27(D15):2253-2256.

Carey L D，Rutledge S A.1998.Electrical and multi parameter radar observations of a severe hailstorm［J］.J Geophys Res，103(D12):13979-14000.

Emanuel K.2005.Increasing destructiveness of tropical cyclones over the past 30 years［J］.Nature，436:686-688.

Lang R，Maurelli A N，Van der Zande W J，et al.2002.Forward modeling and retrieval of water vapor from the Global Ozone Monitoring Experiment:Treatment of narrow band absorption spectra［J］.J Geophys Res，107，4300.doi:10.1029/2001JD001453.

Lay E H，Jacobson A R，Holzworth R H，et al.2007.Local time variation in land/ocean lightning flash density as measured by the World Wide Lightning Location Network［J］.J Geophys Res，112，D13111.doi:10.1029/2006JD007944.

Rodger C J，Brundell J B，Dowden R L.2005.Location accuracy of VLF Worldwide Lightning Location(WWLL)network:Post-algorithm upgrade［J］.Ann Geophys，23:277-290.

Sergio F A，Kristen L C，Thomas J G Jr.2010.An evaluation of the Worldwide Lightning Location Network(WWLLN)using the National Lightning Detection Network(NLDN)as ground truth［J］.J Geophys Res，115，D18206.doi:10.1029/2009JD013411.

Shao X M，Harlin J，Stock M，et al.2005.Katrina and Rita were lit up with lightning［J］.EOS，Trans Amer Geophys Union，86:398-399.

Squires K，Businger S.2008.The morphology of eyewall lightning outbreaks in two category 5 hurricanes［J］.Mon Wea Rev，136:1706-1726.

Thomas J N，Solorzano N，Cummer S A，et al.2010.Polarity and energetics of inner core lightning in three intense North Atlantic hurricanes［J］.JGeophysRes，115，A00E15.doi:10.1029/2009JA014777.

Webster P J，Holland G J，Curry J A，et al.2005.Changes in tropical cyclone number，duration，and intensity in a warming environment［J］.Science，309:1844-1846.