毛志远 付丹红 黄彦彬* 李光伟 敖 杰 蔡杏富
1)(海南省人工影响天气中心， 海口 570203) 2)(海南省南海气象防灾减灾重点实验室， 海口 570203) 3)(中国科学院大气物理研究所， 北京 100029) 4)(海南省三沙市气象局， 三沙 573100)
台风是影响我国的重要气象灾害之一[1]，登陆台风除了直接带来暴雨、强风外，还会导致由极端降水诱发的洪涝和山体滑坡等次生灾害，造成重大人员伤亡与经济财产损失。2014年超强台风威马逊(1409)在登陆的3 d时间里，造成广东、广西和海南等地区共4.68×105hm2农作物受损、1100万人受灾，经济损失约265.5亿元[2-3]。2019年在浙江登陆的超强台风利奇马(1909)引发暴雨、泥石流等灾害，受灾人数达到1400万人，直接经济损失515亿元[4]。
近年随着探测技术的发展，卫星、飞机、雷达和地面雨滴谱仪等逐渐成为观测台风云系三维结构特征和降水过程的主要手段[5-17]。Chen等[18]对登陆我国福建的台风莫拉克(0908)滴谱特征进行分析，发现眼墙和外围雨带的微物理特征存在显著差异，眼墙降水滴谱分布比外围雨带更宽，对比不同降水阶段以及整个降水过程，滴谱分布均呈单峰型，峰值为0.4～0.7 mm。超强台风苏迪罗(1513)从福建登陆后北移，对2015年8月10日影响安徽产生降水的雨滴谱特征参量分析显示，降水云系主要为层状云，降水粒子以小粒子为主，平均粒径为0.6～0.85 mm，直径为1～1.25 mm的粒子对雨强贡献最大[19]。林文等[20]分析了台风麦德姆(1410)在福建省福清市高山镇登陆后不同部位的降水强度和雨滴谱特征，发现台风右侧云系雨滴平均谱谱宽由宽变窄，小滴数浓度先增后减，大滴逐渐减少，移动路径中轴后侧残留云系出现大滴数浓度和谱宽突增，雨强小于10 mm·h-1时，以小粒子贡献为主，雨强大于10 mm·h-1时，大滴贡献随雨强增大而增大。申高航等[21]选取4个强降水中心雨滴谱数据，对比超强台风利奇马(1909)影响山东降水的雨滴谱特征，结果表明：不同降水中心雨滴谱演变特征差异明显，但雨滴尺度谱和速度谱主要为双峰型，且尺度谱双峰中不同降水中心在大尺度峰值对应一致，小尺度峰值不一致。李欣等[22]对比台风利奇马(1909)和台风巴威(2008)北上期间影响青岛产生强降水的形成机制和滴谱特征，发现两次降水中雨滴的碰并增长对降水贡献最显著，随着雨强增大，小雨滴数浓度下降，而中等以上尺度的雨滴谱占比增加。
平均每年有6～7个台风影响海南岛，其中登陆台风为2.6个[23]，带来强风、暴雨灾害的同时，台风降水占海南岛年降水量的30%，对缓解干旱有积极作用[24]。目前，针对影响和登陆海南岛的台风研究主要集中在天气形势诊断、数值模拟及气象灾害统计[25-31]。2018年16号台风贝碧嘉(1816)对海南岛影响时间长达8 d，移动缓慢，长时间维持较高强度，造成海南岛大范围受灾和巨额经济损失。本文利用S波段单偏振多普勒天气雷达CINRAD/SA和OTT-Parsivel二代雨滴谱仪分析2018年8月14—16日台风贝碧嘉(1816)降水云系结构与降水形成特征，分析台风不同位置和不同云系结构下微物理量演变和雨滴谱分布特征，获取台风外围云系微物理特征，这对台风外围云系降水估测以及对台风带来的灾害天气预报预警有重要作用。
1.1 台风贝碧嘉(1816)概况
2018年8月9日08:00(北京时，下同)南海海域有热带云团加强为热带低压，并于10日09:00在海南省琼海市潭门镇沿海地区登陆，先后经过琼海、定安、海口和文昌等市县，10日下午从文昌市东北部出海。11日10:00在广东省阳江市海陵岛再次登陆，12日14:00在广东西部近海加强成为2018年第16号台风贝碧嘉(热带风暴级，编号1816)。14日台风贝碧嘉(1816)在广东沿海徘徊时，此时南亚高压逐渐减弱，西西伯利亚冷涡逐渐东移，副热带高压脊减弱东退至120°E附近，17:00台风贝碧嘉(1816)加强为强热带风暴，14日夜间向西偏南方向移动。15日随着南亚高压的持续减弱，副热带高压(简称副高)加强西伸，台风贝碧嘉(1816)位于副高南侧，台风贝碧嘉(1816)北侧偏东风加强，在北侧偏东风引导下，台风贝碧嘉(1816)向雷州半岛方向移动，21:40在广东省雷州市沿海登陆，登陆时由强热带风暴级减弱为热带风暴级，16日02:00移入北部湾，随后在北部湾东部海面再次加强为强热带风暴级，之后稳定西行。
1.2 台风贝碧嘉(1816)降水过程演变特征
8月14—16日受台风贝碧嘉(1816)影响，海南岛出现降水过程，强降水主要出现在海南岛北部地区。14日00:00—15日00:00海南岛24 h累积雨量最大值约为50 mm，南部地区降水相对较弱，24 h累积雨量约为10 mm(图1a)。随着台风发展加强，降水过程逐渐增强，15日00:00—16日00:00北部地区24 h累积雨量超过250 mm，其中海口和屯昌地区24 h累积雨量分别为120.9 mm和39.1 mm(图1b)。由逐小时雨强(图2)可以看到，海口站和屯昌站降水过程主要分8月14日00:00—16:00和14日20:00—16日00:00两个阶段，其中14日00:00—16:00的降水过程相对较弱，两站最大雨强分别为4.2 mm·h-1和1.8 mm·h-1，14日20:00—16日00:00降水过程相对较强，海口和屯昌最大雨强可达19.3 mm·h-1和5.7 mm·h-1，此时台风贝碧嘉(1816)已由热带风暴级加强为强热带风暴级，移入广东省湛江市东南方的南海海面上。15日11:00台风中心位于20.9°N，111.9°E，近中心最大风力为25 m·s-1，中心气压为925 hPa，稳定西行向雷州半岛靠近，对海南岛降水的影响最为显著。
雷达观测数据来源于海南省海口市永庄气象观测站(20.0°N，110.3°E，简称海口站)的S波段单偏振多普勒天气雷达CINRAD/SA。两台OTT-Parsivel二代雨滴谱仪分别位于海口站和屯昌县气象观测站(19.4°N，110.1°E，简称屯昌站)。雨滴谱仪采样面积为18×3 cm2，等效体积直径(D，单位：mm)测量范围为0～26 mm，按非等间距间隔将其分为32个通道，速度(V，单位：m·s-1)测量范围为0～22.4 m·s-1，同样按非等间距分为32个速度通道，每次采样时间间隔60 s，一次采样间隔内的粒子谱数据为32×32=1024个。2018年8月14日17:00 台风贝碧嘉(1816)加强为强热带风暴级后向西偏南移动，从广东省雷州市登陆直至进入北部湾期间，位于其移动路径左侧的海南岛受螺旋雨带影响产生的强降水过程的雨滴谱特征进行观测，其中海口站位于螺旋雨带内侧，相比屯昌站更靠近台风移动路径。
3.1 台风贝碧嘉(1816)降水云系结构特征
图3为2018年8月14日22:09—16日00:35海口站雷达组合反射率因子。由图3可见，14日22:09台风外围雨带的前段已经影响海南岛东北部，整体以层状云为主，雨带北部嵌有回波较强的对流性云块。至23:35，海口站和屯昌站影响云系依然表现为台风外围云系的层状云，雷达组合反射率因子约为30 dBZ。15日01:38影响海口站和屯昌站的降水云均以层状云为主，伴随弱对流过程，海口站组合反射率因子为35 dBZ，屯昌站为25 dBZ。随着台风中心逐渐向西移动，15日03:35外围雨带影响范围向西南部扩大，影响海口站和屯昌站的降水云系仍以层状云为主；
07:34雨带影响范围覆盖整个海南岛，开始有对流性云块影响海口地区，最大组合反射率因子达60 dBZ；
此时影响屯昌站的以层状云为主，组合反射率因子稳定在30 dBZ。15日11:22 外围雨带整体结构变得松散，逐渐向西北移出海南岛，随着台风主体向雷州市靠近，19:34外围雨带主体对海口站和屯昌站的影响减弱，此后两个站点的降水主要来自组合反射率因子低于25 dBZ的尾部残留云团，15日21:40贝碧嘉(1816)在雷州市登陆，云团强度进一步减弱。台风贝碧嘉(1816)对海南岛降水的影响主要为台风登陆前的外围云系，台风在雷州市登陆后对海南岛的影响减弱。
3.2 台风贝碧嘉(1816)降水云系雨滴谱微物理特征
受设备自身限制和环境影响，OTT-Parsivel雨滴谱仪在观测大雨时，无法解决雨滴重叠问题，样本数量超过计数器上限将出现饱和溢出现象，因此对Parsivel雨滴谱仪观测数据质量控制尤为重要[32-36]，分别剔除原始数据中D>8 mm的雨滴粒子，因为在自然降水中不会出现大于该粒径的大雨滴。根据图2站点雨量数据演变，较强降水过程主要集中在8月15日03:35—8月16日00:36，利用该时段观测采集的雨滴谱数据计算并分析该台风影响下降水的雨滴谱特征，计算公式见文献[18]。
3.2.1 降水过程雨滴谱微物理特征演变特征
图4为8月15日海口站和屯昌站雨滴总数浓度、雨强(I,单位：mm·h-1)、雨滴最大直径、雨滴平均直径和雨滴谱逐时间演变。结合图3可知，8月15日海口站03:35—07:35主要受层状云影响，03:35海口站以北有显著的对流回波，但海口站至屯昌站以层状云为主，最强回波出现在5 km高度，组合反射率因子约为35 dBZ。03:35—07:35，I≤30 mm·h-1，雨滴平均直径约为1.0 mm，雨滴谱谱型为单峰，峰值出现在D≤1 mm的小雨滴端，雨滴最大直径达到7.5 mm，仅出现在04:11，04:35和04:36，且数浓度小于0.5 m-3·mm-1，该时段雨滴总数浓度约为400 m-3。07:35后高回波强度的对流性云团进入海口站观测范围，海口地区上空附近云系的最大组合反射率因子达60 dBZ，有45 dBZ左右的强回波悬垂。对流云团影响期间，海口站雨强和雨滴总数浓度显著增加，D>7 mm的大粒子频繁出现且数浓度增加至1 m-3·mm-1以上，雨滴总数浓度和雨强在11:07达到最大值，分别为1106.02 m-3和115.58 mm·h-1。此时雨滴谱仍保持单峰型，D≤1 mm的小雨滴浓度迅速增大，10:48D≤1 mm小雨滴数浓度瞬时值超过4500 m-3·mm-1，整体谱型拓宽，D为4～8 mm的大粒径雨滴数浓度增加。随着台风贝碧嘉(1816)向雷州市靠近，11:03外围雨带整体结构已出现松散，海口站附近的云系回波组合反射率因子约为50 dBZ。随着云系主体向海南岛西北部移动，15日11:37—16日00:36海口站附近平均组合反射率因子低于35 dBZ，对流强度减弱，组合反射率因子40 dBZ的弱对流和层状云交替影响，雨滴总数浓度为400～800 m-3，滴谱明显变窄，D>2 mm的中大雨滴数浓度减小，以小雨滴为主，降水减弱，雨强降低至30 mm·h-1以下。
影响屯昌站的降水云系均为组合反射率因子约为35 dBZ的层状云，由图4可见，屯昌站的雨强小于海口站，平均直径约为0.8 mm，雨滴最大直径不超过5 mm，各微物理特征参量演变曲线的峰谷值差距不大，降水分为3个时段：第1段为15日03:35—11:35，第2段为15日11:35—19:35，这两个时段降水雨滴总数浓度和雨强曲线峰谷值变化对应较好，雨滴最大直径为4～5 mm，最大雨强5.81 mm·h-1出现在06:58，两个时段的雨滴谱型以单峰为主，多峰出现次数较少，峰值分布在D<1.5 mm的区间，谱型较海口站明显偏窄，谱宽变化不明显，雨滴浓度集中在D<1.5 mm的小雨滴端；
第3段降水出现在19:35后，此时屯昌站上空雷达回波小于25 dBZ，雨滴平均直径无变化，但雨滴最大直径平均值降至1.35 mm，D>2 mm的雨滴数浓度显著减少，谱宽进一步缩小。23:35后总数浓度出现大幅度提升，峰值出现在D≤1 mm的小雨滴端，雨滴总数浓度突增至2884.48 m-3，最大雨滴D<2 mm，但雨强并未显著增加，平均雨强为1.1 mm·h-1。
3.2.2 雨滴平均分布谱特征
利用雨滴谱数据，分别计算两个站点的平均谱(图5)，由图5可见，海口站和屯昌站平均谱呈单峰型，海口站谱宽更宽，雨滴最大直径达到7.5 mm，屯昌站雨滴最大直径仅为4.25 mm。两个站点的平均谱峰值均出现在D≤1 mm，海口站的峰值出现在0.812 mm，屯昌站的峰值出现在0.437 mm，说明两个站点降水均以小雨滴为主，但在D<0.5 mm的范围内，屯昌站相比于海口站雨滴数浓度更高；
随着雨滴直径增大，屯昌站雨滴数浓度迅速减小，但海口站仍存在较低浓度的大雨滴。综合分析，此次台风影响下的降水过程，大雨滴在下落时破碎生成更多小雨滴，导致小雨滴在总雨滴的数浓度中占主要地位。
为进一步分析降水平均谱在不同雨强下的分布特征，将整个台风降水过程的雨滴谱数据按照雨强I分为6个等级：0 由图6可见，随着雨强的增大，两个站点不同雨强条件下的谱宽逐渐增大，海口站和屯昌站不同等级雨强的平均谱谱型基本保持单峰。平均谱同一个雨强范围内海口站的谱宽比屯昌站更宽，海口站的最大谱宽出现在10 mm·h-1 3.2.3 不同粒径雨滴对数浓度和雨强的贡献 将雨滴直径D分为4个等级：0 由图7可见，两个站点雨滴数浓度在总数浓度中的占比随直径的增大而减小，0 Ulbrich[37]提出的Gamma分布引入形状因子反映雨滴变形对雨滴谱的影响，相较于Marshall等[38]提出的M-P分布更能客观描述雨滴谱谱分布，对各类降水均有较好适应性，被广泛应用于云物理研究中[39-42]。利用Gamma分布对本次降水过程的雨滴谱进行拟合，表达形式为 N(D)=N0DμeλD。 (1) 式(1)中，N(D)(单位：m-3·mm-1)为空间上单位粒径间隔和体积内的雨滴数浓度，N0(单位：m-3·mm-1-μ)是与离子浓度相关的参数，λ(单位：mm-1)是雨滴谱分布的斜率因子，μ是形状因子，量纲为1。 阶矩法广泛应用于估算Gamma分布的3个参量N0，λ和μ，郑娇恒等[43]认为阶矩法估算分布参量时不同矩量的选择对拟合效果有影响，使用低阶或高阶矩量改善某个量的拟合精度同时会造成稳定性降低，中间阶量的拟合效果更佳。因此，本文采用二、四、六阶矩对Gamma分布的3个参量进行估算[44]，并分别对海口站和屯昌站的雨滴谱进行拟合对比： 海口站：N(D)=517.198D-0.025e-1.43D， (2) 屯昌站：N(D)=7518.575D1.487e-3.847D。 (3) 式(2)和式(3)给出了海口站和屯昌站的拟合结果。同时，对比图8拟合谱与观测谱分布可以发现，不论是海口站还是屯昌站，D<0.5 mm处的拟合谱对雨滴数浓度有所高估，其中海口站在直径约为1 mm的观测值略高于拟合值； Gamma函数中，N0，λ和μ并不是相互独立的，其中λ和μ存在二项式函数关系，N0与μ直接相关，且不同地形条件和降水系统间的函数关系也不尽相同[41,45]。因此，获取本次降水λ和μ的关系，不仅能提高对海南岛台风降水滴谱分布特征的认识，也可为数值模式中降水滴谱分布公式的本地化提供参考。 分别对海口站和屯昌站的λ和μ进行拟合，图9为两个站点λ和μ的散点分布和拟合关系。由图9可知，两站点的λ均随μ的增大而增大，且随着μ的增大,屯昌站的离散程度大于海口站； 研究表明： 1) 台风贝碧嘉(1816)降水在时间和空间上有显著不对称性：海口站受对流云块影响，特征参量曲线起伏明显，时间分布不均，表现为阵性强降水； 2) 整个降水过程海口站和屯昌站的平均谱呈单峰型，峰值出现在D≤1 mm的小雨滴端。海口站具有更宽的谱分布，屯昌站在D≤1 mm范围内存在更多雨滴。不同等级雨强下，海口站和屯昌站的平均谱基本保持单峰，随着雨强增大，两个站点的谱宽拓宽，各尺度雨滴的数浓度增加，其中海口站D>3 mm中大粒径雨滴端的增幅更明显，说明中大粒径雨滴出现和数浓度显著增多，且伴随其他粒径雨滴数浓度的增加是导致雨强增大的主要原因。 3) 海口站和屯昌站D≤1 mm的雨滴数浓度在总数浓度中占比超过50%，但对总雨强的贡献小。屯昌站1 mm 4) 海口站和屯昌站的Gamma函数拟合结果较好反映台风贝碧嘉(1816)降水的雨滴谱特征，斜率因子λ随形状因子μ增加而增加，符合二项式分布。 本文基于台风贝碧嘉(1816)降水过程的观测数据，得到海南岛台风降水不同影响云系下的滴谱分布特征，得到一些有意义的结论。但对降水云中的微物理演变过程还需要借助卫星观测和数值模式等方法进一步探讨。
而海口站在该粒径范围的总数浓度为3610 m-3，占比为3.08%，对总雨强的贡献却占56.61%。
在D>1 mm处，海口站和屯昌站的拟合结果均与观测相近，对屯昌站雨滴谱的拟合效果优于海口站。整体而言，在D≤1 mm的小雨滴端，拟合值与观测值略有差异，但拟合谱曲线的分布趋势和峰值变化与观测基本一致，拟合值与观测值相关系数分别为0.79和0.92(达到0.05的显著性水平)，相关较显著，能较好反映本次台风降水的雨滴谱分布。
海口站具有更宽的μ值范围，相同μ值条件下屯昌站λ值更大。两个站点的λ和μ均存在较好的二项式关系，拟合结果的相关系数分别为0.98和0.93(达到0.05显著性水平)。
屯昌站影响云系稳定为层状云，雨滴谱各特征参量曲线波动起伏较海口站小且均匀和连续。在雨滴谱演变上，海口站雨滴谱保持单峰，谱宽较宽；
屯昌站雨滴谱以单峰为主，双峰较少，谱宽较窄。两个站点的雨滴谱随雨强增大均在1 mm的直径范围出现数浓度增加，同时谱宽迅速拓宽，大粒径雨滴数浓度增加。