1、引言

台风是造成浙江地区暴雨灾害的重要天气系统。台风源自海洋，是复杂海气相互作用的产物[1]。按台风中心和暴雨落区的相对位置，台风暴雨可分为台风本体暴雨和台风远距离暴雨[2,3]。台风远距离暴雨具有降雨量大、突发性强和降水时段集中等特点[4]，在实际预报服务中难度很大。近几十年来，台风远距离暴雨研究取得不少新成果，一些研究表明，台风远距离暴雨常与西风槽以及北方冷空气南下等有关[5,6]。在有利的大尺度环流背景下，中小尺度系统活动尤其是中尺度对流云团在台风倒槽内发展是造成远距离暴雨的一个重要原因[7]。刘文婷[8]通过研究1410号台风“麦德姆“，指出在对流不稳定以及弱的条件对称不稳定环境下，对流单体可能会沿着辐合线爆发造成暴雨；由于秋季冷空气逐渐开始活跃，高安宁等[9]认为秋台风在暖湿气流和入侵的冷空气影响下形成辐合，导致温湿度梯度增大，促进了中尺度对流的发生和发展，对暴雨起到增强作用；黄莉等[10]认为秋季适当冷空气的入侵对台风暴雨有明显的增幅作用。另外，蒋尚城等[11,12]认为，台风东侧低空急流在中低纬度系统相互作用引发的暴雨中起重要作用；当台风外围出现东南急流，同时地面有波动发展时，应警惕暴雨的发生[13,14]。当然水汽输送对台风暴雨也很重要，源源不断的水汽输送有利于台风环流的维持及雨带中对流活动的加强，从而促进暴雨的发展[15,16,17]。以上工作从不同角度给出了远距离台风暴雨，特别是秋季远距离台风暴雨的形成原因及相关分析方法。

1949—2017年登陆我国广东的秋台风（9—11月）共有70个，平均每年1个左右。期间引起浙江东北部较大范围出现50 mm以上降雨量的共有10个台风，占秋季登陆广东台风的14%，相当于每20 a出现3个，属于小概率事件。其中有7个造成了区域性大暴雨，包括最近10 a出现的3个，表明与此类台风相关的浙江大暴雨事件有增多趋势。由于此类过程累积雨量往往较大，再加上当时已进入汛期尾声，当地水库蓄水比较丰沛，防汛压力大，对雨量预报要求高，所以预报服务的难度非常高。

本文利用自动站常规资料、卫星云图以及美国国家环境预报中心（National Centers for Environmental Prediction,NECP）的1°×1°再分析资料等，对1822号南海台风“山竹”登陆广东减弱期间浙江东北部出现大暴雨的成因进行诊断分析，以期了解此次大暴雨过程的触发机制和物理过程，并对以后类似暴雨过程的预报有所帮助和启发。

2、“山竹”概况和浙江雨情

2.1“山竹”概况

2018年第22号台风“山竹”于9月16日17时（北京时，下同）前后以强台风在广东省江门市台山沿海登陆，登陆时中心附近最大风力14级（45 m/s），登陆后往西北偏西方向移动，强度逐渐减弱，17日20时停止编报。“山竹”是当年登陆我国的最强台风，具有范围大、强度强和台风倒槽明显等特征。

2.2雨情

浙江东北部部分地区于16日下午开始出现阵雨，17日下午雨渐止，强降水集中时段出现在16日夜里—17日中午前后。16日16时—17日16时浙江东北部普遍出现了暴雨（见图1a），其中嘉兴和宁波普降大暴雨。嘉兴面雨量122.2 mm，最大为海宁223.9 mm，雨量超过200 mm有3个站，超过100 mm有74个站；宁波面雨量113.9 mm，最大为象山麦地山294.6 mm，雨量超过250 mm有2个站，超过200 mm有16个站，超过100 mm有193个站。本次过程有降水集中、雨强强和雨量大等特点。

小时雨量分布图表明16日23时起海宁站降水开始增强（见图1b），一直到17日04时前后雨势减弱，小时雨强在20 mm以上的强降雨持续了4个多小时，其中17日01—02时小时雨强达94.2 mm。强对流云团东南移后，17日白天宁波中北部普遍出现1 h雨强超过20 mm的短时强降水，象山麦地山站10—13时3 h累计降水279.8 mm（见图1c）。可见雨强特别强。

图1 2018年9月16日16时—9月17日16时浙江省雨情与地区雨量分布

2.3降水性质分析

本次强降水过程伴有雷电（见图2a）。9月16日16时—17日16时，嘉兴共发生雷电1 372次，其中正地闪116次，负地闪1 256次，闪电主要发生在17日00—03时，3 h共发生闪电872次；宁波共发生雷电6 733次，其中正地闪1 213次，负地闪5 520次，从17日02时起开始增多，13时以后明显减少，其中17日09—13时，4 h共发生闪电4 379次。强雷电的发生区域、集中活动期和强降水出现的时间以及落区基本一致，表明此次强降水过程和对流系统的活动相关性较大。

9月17日02时海宁站1 h雨强达94.2 mm,12时象山麦地山站1 h雨强达180.8 mm。分别对海宁站17日02时（见图2b）及象山麦地山站17日12时（见图2c)0.5°仰角雷达回波作垂直剖面，由图可见，海宁站附近最强回波达52 dBz，象山麦地山站附近最强回波达60 dBz以上，两者的回波顶高均超过0℃线甚至超过了-20℃线，达到了强对流的回波高度，两个站强回波的出现时间也和闪电频数的时间分布相对应，进一步表明本次过程是对流性质的降水。

3、暴雨成因分析

浙江东北部的暴雨出现在距离1822号台风“山竹”1 500 km以外的地方，那么这个暴雨是否与台风活动有关？我们先从相对直观的云图进行分析。

3.1对流云团

从风云2E卫星红外1通道的黑体亮度温度（Black Body Temperature,TBB）可见（见图3），“山竹”范围大，南北长约14个纬距，东西宽约12个经距。16日下午起宁波西南侧“山竹”倒槽东北部有分散性对流云团产生（图略），16日20时浙江北部有一明显强对流云团（见图3a），最强云顶亮温达到-70℃以下，随后强对流云团缓慢南压；17日02时强对流云团位于嘉兴南部、杭州湾地区以及宁波东北部（见图3b);17日08时强对流云团位于宁波、舟山及其以东海面上（见图3c），强度略有减弱，最强云顶亮温在-60℃以下，宁波慈溪、镇海、北仑等地在08时前后均出现了1 h超过30 mm的短时强降水，镇海站小时雨强超过10 mm的时间达5 h;17日14时强对流云团主体已南压至台州附近（见图3d），该对流云团最强云顶亮温达到了-70℃，台州北侧的宁海和象山在10—13时受该强对流云团影响，出现了强降水和强雷电，象山站2 h最大雨强达到了175.4 mm;17日16时之后，浙江东北部强降水过程基本结束。云图分析可见，对流云团先是在“山竹”外围东北侧生成，随着“山竹”的西移，对流云团逐渐脱离台风外围云系，并向东南方向移动，先加强后减弱。

图2强降水时段内降水性质分析

从以上分析可见对流云团的发生与“山竹”外围云系有一定的关联，随着“山竹”的西移，其他系统又造成了云团进一步发展和往东南方向移动，那么影响系统具体有哪些呢？

3.2影响系统

3.2.1副高

李春虎等[18]研究认为西行台风和副热带高压（以下简称“副高”）的相互作用会使副高的中心位置西伸，强度增强。9月13日起我国中东部大部分地区都受副高控制，16日08时500 hPa（图略）副高呈块状，脊线在30°N附近，最西伸至108°E附近，一直到17日（见图4），浙江东北部都处于副高控制下。持续副高控制下，浙江东北部以晴热天气为主，气温较高，能量积聚明显。

3.2.2台风倒槽

从地面图上分析，“山竹”北侧有台风倒槽发展，特别是东北侧倒槽发展旺盛，9月16日14时倒槽北伸至江苏南部（图略），浙江东北部位于倒槽东北侧，有明显的东南急流输送和风向辐合。从台风倒槽的位置分析，暴雨区位于倒槽东北侧，与台风有一定的关联。随着“山竹”登陆并西移，17日倒槽减弱。

3.2.3冷空气

9月16日河套以东有高空槽东移，16日20时浙江东北部500 hPa以西南风为主（见图4a),17日08时（见图4b），浙江东北部500 hPa转为西偏北风向，江苏大多为西北风向，说明有弱冷空气渗透。700hPa环流形势（图略）和500 hPa相似。850 hPa上，17日08时浙江与江苏交界处有东西向的弱辐合带，是西北风和西偏南风的辐合，浙江东北部的强降水区域正好位于该辐合线的南侧。地面上16日夜里弱冷空气扩散南下，地面辐合生成并发展，本次强降水过程也与该系统的触发有关。

图3不同时刻1 000 hPa气压场（线条，单位：dapgm）及亮温（阴影，单位：℃）

分析实况地面流场图，16日20时嘉兴海宁附近有一偏西气流和东南气流形成的辐合中心（见图5a）。随着弱冷空气的渗透，23时该处有偏北气流和东南气流的辐合（见图5b),17日05时在宁波的镇海和北仑附近出现了偏北气流和偏南气流形成的辐合中心（见图5c),08时辐合中心南压至象山港附近，且辐合有所加强，12时以后影响浙江东北部的中尺度辐合明显减弱。地面辐合的发展减弱比雨势的增大和减小提前了1～3 h左右。

由影响系统分析可见，在前期副高控制的浙江东北部，连续晴热积聚了能量；“山竹”倒槽东北侧的对流云团在沿海高湿高能区域有所发展；弱冷空气的侵入使得低层辐合加强，造成了对流云团的进一步发展，强降水是由这些系统先后影响共同造成的。

4、物理量诊断

4.1稳定度和能量条件

4.1.1探空物理量

对流云团的发生、发展还可以从本次降水过程的环境物理量特征中得到印证（见表1）。9月16日08时—9月17日08时嘉兴和宁波抬升指数（Lifted Index,LI）在-3～-4℃，有较明显的条件性不稳定度；而K指数普遍在32℃以上，尤其是16日20时和17日08时均在35℃以上，超过了本地强对流发生发展的阈值，表明降水过程中大气层结有明显的不稳定性。16日08时嘉兴和宁波的对流有效位能（Convective Available Potential Energy,CAPE）值均在1 500 J/kg左右，16日20时嘉兴增至1 809 J/kg，宁波增至1 687 J/kg，表明16日白天CAPE积聚，为对流性天气的发生、发展提供了有利条件，该指标对暴雨的预报有12 h左右的提前量。

图4不同时刻形势场图（红色虚线：500 hPa温度场，单位:℃；紫色实线：500 hPa高度场，单位：dagpm；棕色实线：图a为500 hPa槽线、图b为850 hPa切变线；风向杆：500 hPa风场）

图5不同时刻地面流场图

表1 2018年9月16日08时—9月17日08时探空物理量

4.1.2假相当位温θse分布

“山竹”登陆前，江南大部地区由副高控制，气温较高，能量得以积聚。从925 hPa上假相当位温（θse）的水平分布来看，16日14时江苏到苏浙交界处有θse密集带（见图6a），基本呈东北-西南向分布，在密集带西北侧θse较低，有干冷空气；在其东南侧θse较高，有暖湿气流，θse密集带的存在表示有弱的干冷空气向台风倒槽侵入；同时，θse密集带也代表了能量锋的形成，由于能量锋有利于对流系统的发生、发展，因此沿着θse密集带有出现强降水的可能[19]。浙江东北部此时位于该θse密集带前方高值区中心，这阶段有不稳定能量积聚，之后该θse密集带逐渐往南移，促使浙江东北部低层的暖湿空气抬升，导致降水过程增强；浙江东北部地区一直到17日02时θse都普遍在350 K以上，能量较高，此后，θse逐渐减小，925 hPa上θse的增大和减小，比雨势的增大和减小提前了12h左右。

沿121°E，做28°～35°N之间的θse垂直剖面图。由图6b可见，16日14时850 hPa以下的长江以南地区有明显的θse密集带，浙江东北部高温、高湿，淮河流域有干冷空气向南推进，促进了对流的生成和发展，并引起强降水；20时（图略）该θse密集带有所减弱，17日02时（图略）θse密集带南移并继续减弱。垂直剖面上θse密集带的强弱变化比雨势增减提前了12 h左右。

4.1.3湿位涡分析

大暴雨的发生、发展常伴有气旋性垂直涡度的迅速增大[20]，因此研究气旋性垂直涡度的变化可以更好地了解暴雨成因。不考虑ω的水平变化，湿位涡（Moist Potential Vorticity,MPV）在P坐标系中的表达式为[21]:

公式1

式中：ζ为绝对涡度矢量；f为科氏参数；θe为相当位温。第一项（MPV1）包含了惯性稳定性(ζ+f)和对流稳定性的贡献，在北半球绝对涡度一般为正，当大气处于对流不稳定时，MPV1<0，反之则MPV1>0；第二项（MPV2）包括水平风垂直切变和湿斜压性∇pθe。倾斜涡度发展理论[21]表明：在湿位涡守恒的制约下，水平风垂直切变、湿斜压变化和大气垂直稳定度都可导致垂直涡度的发展。林小红等[22]研究提出MPV2可对台风中尺度对流系统（Mesoscale Convective System,MCS）的发展起加强作用；在MPV1<0，即对流不稳定条件下，与暖湿气流相关的MPV2的发展可触发气旋性垂直涡度的增大，并引起对流性强降水的发生和加强[23]。

图6 2018年9月16日14时θse分布（单位：K)

沿着121°E做湿位涡垂直剖面，在台风“山竹”登陆之前，600 hPa以下24°～33°N的我国东南沿海MPV1均为负值，登陆前后大的负值中心在26°N附近，最小值可达-1.0 PVU（图略），随后MPV1负值中心向北移动，强度略有减弱。17日02时30°N附近中低层MPV1有弱的负值（见图7a），表明有对流不稳定能量，若有低层辐合，将促进上升运动的发展，在台风环流水汽输送的配合下，有利于对流性降水的发生、发展；700 hPa以下的低层，MPV2值为正（见图7b），气旋性垂直涡度得以增长，此时正是嘉兴海宁降水最强时段。17日08时30°N附近MPV1的负值区范围有所扩大（见图7c），高度有所增加，MPV2仍为正值，因此强降水维持。17日14时，MPV1负值中心移至30°N附近800 hPa高度上，且值有所增大，主要是由于副高西伸北抬，此时降水逐渐减弱停止，气温上升，大气层结不稳定，而在这一时段内MPV2大值中心略有南移（见图7d），在29°N附近的低层有弱的正值中心，宁波象山地区（约29.5°N附近）刚经历了2～3 h的持续强降水，与上文所提强降水过程有东南移的特点一致。由上文可知，湿位涡能较好的诊断和预报暴雨过程。

4.2水汽及来源分析

为了分析本次大暴雨过程的水汽来源，应用美国国家海洋和大气管理局（National Oceanic Atmospheric Administration,NOAA）等开发的混合单粒子拉格朗日综合轨迹模式（HYbrid SingleParticle Lagrangian Integrated Trajectory,HYSPLIT-4）和全球资料同化系统（Global Data Assimilation System,GDAS）资料[24,25]，对2018年9月16日20时和9月17日08时宁波逐小时的48 h后向水汽源轨迹进行分析（见图8）。由图可见，本次大暴雨发生前和发生时的低层水汽源均来自于台风中心位置的东侧，且随着台风的西移，水汽源位置也逐渐西移，说明本次大暴雨过程中水汽来源和台风有关。

图7 2018年9月17日不同时刻湿位涡沿121°E经向垂直剖面图

图8不同时刻以宁波为中心不同高度48 h水汽后向轨迹

持续不断的水汽供应是暴雨形成的重要条件，取本次强降水区域范围：28°～31°N,119°～123°E，计算该时段该区域4个方向截面（1 000 hPa和100 hPa）的水汽通量收支。计算公式如下（东、南、西、北4个方向的水汽通量收支分别用Fe、Fs、Fw、Fn表示，总的用F表示）[26,27]:

公式2

公式3

公式4

公式5

式中：P0为100 hPa,Ps为表面气压，λ为经度，φ为纬度；总的水汽通量收支公式为：

公式6

图9所选区域4个边界面水汽水平通量及其总和

正值表示该区域水汽流入，负值表示该区域水汽流出。图9中各线是所选区域4个方向水汽水平通量及其总和。从图9可以看到水汽流入总量比较充沛，16日08—20时水汽呈递增趋势，20时水汽流入达最大值，约40×104g/（cm·hPa·s）左右，此后水汽流入总量逐渐下降，17日02—08时有一个快速减小过程，此后呈缓慢波动减小趋势；17日0时起降水逐渐增强，而浙江东北部的强降水在17日13时之后减弱，16时基本结束。16日“山竹”靠近并登陆我国，由于其位于整个关注区域的南侧，南侧水汽流入逐渐增多，16日20时之后成为水汽来源最多的地方，最大值超过50×104g/（cm·hPa·s）；东侧水汽在“山竹”登陆后迅速减少，从早期的60×104g/（cm·hPa·s）以上减少至10×104g/（cm·hPa·s）以下；北侧和西侧主要是水汽输出，较稳定的输出在10～20×104g/（cm·hPa·s）左右。

图1 0不同时刻850 hPa水汽通量（彩色阴影）、水汽通量散度（彩色阴影）及风场（风向杆）

通过分析可知，水汽收入主要来自关注区域的东面和南面，水汽来源主要来自台风环流，偏东和偏南气流的输入对浙江东北部的强降水起了很大作用，水汽的增减比雨势的增减有12 h左右的提前。

通常情况下，水汽通量和水汽通量散度能较好地反映水汽大小、输送方向及辐合情况，强降水的持续也和急流有较好的相关性[28]。16日14时，“山竹”及其倒槽所在区域水汽通量值较大（图略），浙江大部分地区在台风外围环流中，浙江和福建交界处有一水汽通量大值区，最大值可达28 g/（cm2·hPa·s）。图10a是9月16日20时850 hPa的水汽通量分布，浙江处于东南风急流中，闽浙交界处的水汽通量大值区已有所北抬，强度略有减弱，最大值达到20 g/（cm2·hPa·s），输送到浙江东北部的水汽通量在8～16 g/（cm2·hPa·s）左右；17日00时前后强降水过程开始，17日02时（图略）该处水汽通量仍维持较大值，东南气流的持续输入是造成强降水的原因之一，17日08时（图略）水汽通量值开始迅速减小，14时浙江东北部水汽通量基本在4 g/（cm2·hPa·s）以下（见图10b），强降水过程已经减弱并趋于结束。

9月16日20时，浙江东北部850 hPa的水汽通量散度也出现了一个弱的负值中心（见图10c），极小值约为-4×10-5g/（cm2·hPa·s），且还有东南风的风速辐合（蓝色椭圆形内），这都有利于水汽的辐合；17日08时宁波地区水汽通量散度极值达到了-6×10-5g/（cm2·hPa·s），此后6 h宁波大部地区都出现了强降水，随后东南急流减弱，水汽通量散度减小（见图10d），降水也随之减弱并结束。

由此可见，此次强降水过程850 hPa水汽通量和水汽通量散度场的变化对强降水的大小和落区有一定的指示效应，其增减相比雨量的增减有6 h左右的提前。

5、总结与讨论

统计分析表明登陆广东的台风虽然离浙江较远，但台风倒槽可能给浙江带来暴雨，特别是秋台风，如有冷空气结合，会使降雨出现增幅。通过分析2018年9月16日16时—17日16时强降水过程以及台风“山竹”所起的作用，得到如下结论：

(1）本次强降水过程具有降水时间集中、雨强强、累积雨量大等特点，并伴有强雷电；强降水发生在台风“山竹”登陆广东后减弱过程期间，“山竹”台风倒槽为对流云团在浙江东北部产生提供了条件，副高提供了能量，弱冷空气的侵入使得低层辐合加强，使对流云团进一步发展，强降水是这些系统先后影响共同造成的。

(2）通过对湿位涡的分析，发现在对流不稳定条件下，即MPV1<0时，与暖湿气流相联系的MPV2的发展，触发了垂直涡度的增长，促进了本次大暴雨过程的发生和发展。强降水发生在能量锋区附近，且浙江东北部925 hPaθse值上升至350 K以上；能量场的位置和梯度大小对强降水预报有12 h左右的提前量。

(3）本次大暴雨过程的水汽主要来自台风东侧，充沛的偏东和偏南气流为大暴雨的发生、发展提供了必要的水汽条件；低层水汽辐合，高层水汽辐散，使上升运动发展，并促进强降水的形成，此次大暴雨过程中低层水汽通量分布和水汽通量散度的增减相比雨量的增减有6 h左右的提前。

(4）在预报实践过程中，虽然主要侧重台风倒槽、冷空气、低空急流等位置和相互作用的分析，但是对当地层结的稳定性判断，以及水汽、能量等的输入情况及变化也需要格外关注。

参考文献：

[1]赖巧珍,马雷鸣,黄伟,等.台湾岛附近海洋对0908号台风“莫拉克”的响应特征[J].海洋学报,2013,35(3):65-77.

[2]陈联寿,丁一汇.西太平洋台风概论[M].北京:科学出版社,1979:8.

[3]陈联寿.热带气旋研究和业务预报技术的发展[J].应用气象学报,2006,17(6):672-681.

[4]单磊,谭桂容,姚叶青,等.一次台风远距离暴雨水汽条件及输送过程研究[J].热带气象学报,2014,30(2):353-360.

[5]杨金锡,洪吉.能量锋生与台风倒槽暴雨[J].南京气象学院学报,1986(1):47-54.

[6]蒋尚城,林楠.85年9号台风与辽宁特大暴雨的卫星云图分析[J].北京大学学报(自然科学版),1988,24(3):351-361.

[7]陈德花,张伟,尹烈,等.“海贝思”残留螺旋云带的中尺度暴雨成因分析[J].暴雨灾害,2015,34(2):133-142.

[8]刘文婷.登陆台风Matmo(2014)北侧对流雨带的发生及发展机制研究[D].杭州:浙江大学,2018.

[9]高安宁,张瑞波.2011年秋季强台风“纳沙”导致广西灾害成因分析[J].灾害学,2013,28(1):54-58.

[10]黄莉,白龙,李紫甜.适当冷空气对秋季台风暴雨增幅作用研究[J].灾害学,2018,33(2):38-44.

[11]蒋尚城,谢安.西太平洋台风与北方暴雨[C]//北方灾害性天气文集.北京:气象出版社,1981:6-10.

[12]蒋尚城.远距离台风影响西风带特大暴雨的过程模式[J].气象学报,1983,41(2):147-158.

[13]徐威,周顺武,葛旭阳,等.西北太平洋热带气旋快速增强阶段的风速分布特征[J].热带气象学报,2017,33(2):259-266.

[14]陈有利,钱燕珍,潘灵杰,等.一次与台风相关联的浙江东北部暴雨成因及预报难点分析[J].干旱气象,2018,36(2):272-281.

[15]李英,陈联寿,徐祥德.水汽输送影响登陆热带气旋维持和降水的数值试验[J].大气科学,2005,29(1):91-98.

[16]冀春晓,陈联寿,赵放.登陆台风Matsa维持机理的数值研究[J].气象学报,2007,65(6):888-895.

[17]周玲丽,翟国庆,王东海,等.0713号“韦帕”台风暴雨的中尺度数值研究和非对称性结构分析[J].大气科学,2011,35(6):1046-1056.

[18]李春虎,黄福均,罗哲贤.台风活动对副热带高压位置和强度的影响[J].高原气象,2002,21(6):576-582.

[19]曹宗元,陈淑琴,刘飞,等.相似路径台风“天鹅”(1515)和“珊珊”(0613)降水差异分析[J].海洋预报,2018,35(3):48-56.

[20]陶诗言.中国之暴雨[M].北京:科学出版社,1980:1-10.

[21]吴国雄,蔡雅萍,唐晓菁.湿位涡和倾斜涡度发展[J].气象学报,1995,53(4):387-405.

[22]林小红,吴建成,刘通易,等.三个典型登闽空心结构台风强降水分布差异分析[C]//第35届中国气象学会年会S1灾害天气监测､分析与预报.福州:福建省气象台,2018.

[23]郭宇光,钱燕珍,方艳莹,等.“杜鹃”登陆减弱后所致的宁波大暴雨成因分析[J].暴雨灾害,2018,37(4):356-363.

[26]方艳莹,钱燕珍,郭宇光,等.“灿鸿”台风大暴雨的水汽诊断分析[J].浙江气象,2016,37(3):7-10,26.

[27]刘家峻,聂新旺,梅娜,等.超强台风“威马逊”强度变化条件的诊断分析[J].海洋预报,2016,33(6):51-56.

[28]高晓梅,王令军,闫凤君,等.台风“达维”影响山东的暴雨落区及路径特征分析[J].海洋预报,2015,32(2):72-78.

方艳莹,钱燕珍,申华羽,吕劲文,薛国强.1822号台风“山竹”引起浙江东北部大暴雨成因分析[J].海洋预报,2020,37(04):86-96.

基金:宁波市“科技创新2025”重大专项(2019B10025);浙江省宁波市科技局项目(2014C50024)