1. 引言
暴雨是宿迁夏季主要灾害性天气,其引发的内涝、洪涝灾害对区域社会经济与民生安全威胁显著[1]。宿迁地处江淮北部气候过渡带,夏季暴雨受东亚夏季风、副热带高压及中纬度低值系统共同影响,降水系统复杂且微物理过程具独特性[2],但当前研究多聚焦宏观环流,对雨滴谱(DSD)特征及分云型(层状云、积雨云、层积混合云) Z-R关系分析较少,制约暴雨微物理认知与雷达定量降水估测(QPE)精度提升[3]。
雨滴谱直接反映降水微物理机制,地域与云型差异性显著。郑州“7・20”极端暴雨中,积雨云降水最强时段粒子数浓度与质量加权平均直径(Dm)极高,Dm达2.82 mm,直接致极端雨强[1];淮北地区观测显示,层状云降水以直径 < 2 mm的雨滴为主,积雨云中雨滴对雨强贡献超75% [4];南京地区证实,层状云雨滴谱参量稳定、谱宽最窄,积雨云谱宽最宽,三类云Gamma分布参数差异显著[5];安宁河谷对流云粒子数浓度是层状云3.7倍,90%降水强度来自直径 > 1 mm粒子[6];广东东莞不同云型降水参数遵循“积状云 > 混合云 > 层状云”规律,层状云以单峰谱为主[7];海南屯昌三类云降水均以小雨滴为主,对流云中大雨滴对雨强贡献显著[8],这些成果凸显宿迁分云型研究必要性。
Z-R关系(
)的云型依赖性亦被多地证实。淮北地区层状云与积雨云Z-R关系系数分异明显,层状云呈“小a大b”特征[4];南京地区积雨云Z-R关系系数a值最大[5];庐山观测发现,对流云a值大、b值小[9];张欢等研究表明,对流云降水Z-R拟合效果优于层状云[10];短时强降水与持续性强降水的积雨云Z-R关系对比显示,大雨强下积雨云Z-R关系b值更接近经典大陆性对流云[11]。此外,雨滴下落速度具地域差异,郑州“7・20”暴雨中Atlas经验公式低估大雨滴落速[1];贵州乌蒙山区随雨强增大,2~3 mm雨滴贡献升至70%以上[12];青藏高原东北部雨滴落速标准曲线拟合存在低估[13],而宿迁不同云型雨滴落速规律尚不明确。
Gamma函数拟合雨滴谱的参数差异显著,南京地区三类云Gamma参数不同[5];淮北地区虽以层状云降水为主,但对流云对总降水贡献达75% [4];柳州“龙舟水”降水Gamma函数对小雨滴拟合优于大雨滴[14];昭苏夏季降水Gamma拟合也具云型差异[15],曾庆伟等研究为雨滴谱参数化提供方法支撑[16]。
鉴于此,本研究利用2022~2024年宿迁夏季DSG1雨滴谱仪观测数据,结合多源资料筛选9次暴雨样本,分析不同云型雨滴谱特征、Gamma函数分布和Z-R关系,旨在揭示宿迁暴雨微物理差异,为雷达分析提供支撑,补充沿淮过渡带数据[17]-[19]。
2. 资料和方法
2.1. 资料总述
观测使用DSG1降水现象观测仪,仪器采样面积0.0054 m2,粒径测量范围0~25 mm,粒子落速测量范围0~22.4 m/s,粒子测量等级分为32种直径和32种速度共计1024种,取样间隔为60 s。
本文选用2022~2024年宿迁夏季(6~8月)雨滴谱资料,结合降水量、雷达和云观测,将24小时降水量50 mm以上且小时雨强不低于0.1 mm/h视为一次暴雨样本,按照云状分为积雨云、层状云和层积混合云暴雨。根据筛选,得到9次暴雨过程,其中2次层状云暴雨(S),4次积层混合云暴雨(T)以及3次积雨云暴雨(C)。
2.2. 处理方法
将仪器每分钟输出的雨滴总数不小于10视为一个降水分钟,剔除直径
和
的粒子。剔除粒子下落速度与Atlas速度经验公式
(1)
计算的速度相差大于5 m/s的数据。
利用Battaglia等的轴比订正法进行轴比订正。
(2)
其中,AR是轴比,Deq是等效粒子直径。
首先计算出单位体积、单位尺度间隔内的雨滴数浓度,
(3)
其中,
是第i个尺度通道的雨滴直径,
是第i个尺度通道、第j个速度通道下的雨滴个数,A (单位:m2)是雨滴谱仪的采样面积,Δt是采样间隔60 s,
是第j个速度通道的下落速度,
是第i个尺度通道的尺度间隔。
(单位:mm−1∙m−3)表示直径介于
与
的单位体积单位尺度间隔内的雨滴数。结合表1中的计算公式,由此可以计算出数浓度(
),雨强(R)、雨水含量(W)和反射率因子(Z);同时计算出粒子微物理结构的不同特征尺度参数:算术平均直径(
)、质量平均直径(
)、体积中值直径(
)、峰值直径(
)、粒子谱宽(
)。其中
是Atlas速度经验公式(详见式(1)),
为水的密度,取1 g/cm3。
选用Gamma函数分布
(4)
对数浓度(
)和直径(D)的谱分布进行拟合,式中,
(单位:
)是截距参数,Λ (单位:mm−1)是斜率参数,μ是形状因子。利用阶矩法估计谱分布的参数。
Table 1. Meanings and calculation formulas of DSD microphysical parameters
表1. 雨滴谱微物理特征量的含义和计算公式
特征量 |
符号 |
单位 |
定义 |
公式 |
数浓度 |
|
m−3 |
单位空间的雨滴总量 |
|
雨强 |
|
mm/h |
单位时间的降水量 |
|
雨水含量 |
|
g/m3 |
单位空间雨滴总质量 |
|
反射率因子 |
|
mm6∙mm−3 |
降水回波强度 |
|
算术平均直径 |
|
mm |
全部雨滴的直径总和除以雨滴的总数 |
|
质量平均直径 |
|
mm |
单位体积内所有粒子直径加权质量相对于粒子总质量的平均直径 |
|
体积中值直径 |
|
mm |
含水量的一半是由半径大于此值的大雨滴所组成的 |
|
峰值直径 |
|
mm |
最大频率直径 |
最大值所对应的直径 |
粒子谱宽 |
|
mm |
粒子最大直径和最小直径差 |
|
3. 结果分析
3.1. 降水微物理特征分析
3.1.1. 雨滴下落速度
图1为不同类型暴雨的雨滴谱尺度–速度散点分布及拟合曲线:利用速度直径公式
得出的黑色拟合曲线,基本处于各尺度箱线图的四分位数之间,且与平均加权下降速度一致性较高,表明拟合效果良好;而Atlas经验速度公式参考曲线位于拟合曲线下方,且随粒子直径增大,两曲线间距扩大,尤其是当D > 1 mm时,经验曲线在四分位数之外,这意味着按Atlas经验速度公式计算的雨强会低估宿迁降水强度,中、大雨滴数浓度贡献率越高,低估程度越大。此外,积雨云(
)、层状云(
)、层积混合云(
)暴雨的雨滴尺度-速度拟合函数存在差异,想更精准的计算雨强,需考虑采用对应的尺度–速度函数。
3.1.2. 微物理特征量的平均特征分析
雨滴谱的微物理特征量可反映降水基本特征。表2是不同类型暴雨雨滴谱微物理特征量的平均值:数浓度(
)范围为381~613 m−3,其中积雨云浓度最大,层积混合云浓度最小;平均雨强(1.96~5.72 mm/h)与含水量(0.12~0.29 g/m3)特征一致,均为积雨云最大、层状云最小;特征直径(算数平均直径
、质量平均直径
、体积中值直径
、峰值直径
)中,
、
、
均为积雨云最大、层状云最小;积雨云和层积混合云的
相同,且均大于层状云的
。
Figure 1. Droplet size-velocity spectrum distribution and boxplots of each size bin for different types of rainstorms ((a) Cumulonimbus cloud, (b) Stratiform cloud, (c) Stratocumulus mixed cloud) (The color scale represents the number of raindrops on a logarithmic scale, the red * represents the measured average weighted velocity, the blue solid line represents the Atlas empirical velocity fitting curve, and the black solid line is the falling velocity fitting curve)
图1. 不同类型暴雨的雨滴尺度–速度谱分布和各尺度档的箱线图((a) 积雨云,(b) 层状云,(c) 层积混合云) (色阶代表对数尺度的雨滴数,红色*代表实测平均加权速度,蓝色实线表示Atlas经验速度拟合曲线,黑色实线是下落速度的拟合曲线)
Table 2. Mean values of DSD microphysical parameters for rainstorms of different cloud types
表2. 不同云类型暴雨雨滴谱微物理特征量的平均值
类型 |
|
|
|
特征直径(mm) |
|
|
|
|
积雨云 |
613.01 |
5.72 |
0.29 |
0.72 |
1.26 |
1.13 |
0.54 |
层状云 |
417.79 |
1.96 |
0.12 |
0.64 |
0.93 |
0.84 |
0.52 |
层积混合云 |
380.96 |
2.91 |
0.16 |
0.70 |
1.08 |
0.98 |
0.54 |
3.1.3. 不同尺寸雨滴对降水的贡献
为明确各尺度雨滴对降雨微物理参量的贡献,将雨滴按直径分为3档:小雨滴(
)、中雨滴(
)、大雨滴(
),分别计算其对总数浓度(
)、总雨强(R)和总回波强度(Z)的贡献率和不同降水档下各尺度粒子对雨强(R)的贡献率。
图2为积雨云、层状云、层积混合云暴雨中不同档位雨滴对
、R、Z的贡献率:宿迁夏季暴雨中,小雨滴数浓度占比最高,层状云达87.3%,积雨云和层积混合云均为80.8%;随粒子尺度增大,积雨云和层积混合云的中、大雨滴数浓度占比高于层状云。尽管小雨滴数占比约80%,但对雨强的主要贡献来自中雨滴,不同类型暴雨中雨滴贡献率为60.4%~67.9%。雷达反射率因子(Z)的贡献中,积雨云以大雨滴为主(65%),层状云和层积混合云以中雨滴为主(分别为64.5%、51.6%);值得注意的是,层积混合云中大雨滴对Z的贡献率达45.7%。因Z与直径(D)的6次方及数浓度成正比,故即使大雨滴数占比不足1%,对Z的贡献仍远大于小雨滴。综上,宿迁暴雨以小雨滴为主,中雨滴是雨强的主要贡献者,中、大雨滴则主要影响雷达反射率因子。
图3为不同云类型暴雨各档位雨滴对不同强度雨强的贡献率:不同暴雨类型中小雨滴对雨强的贡献率随降水强度的增大呈现减小的趋势,大雨滴的贡献率呈增大趋势;当雨强小于1 mm/h时,对降水率起最大贡献的为小雨滴,层状云中最大为92.3%;在雨强大于5 mm/h的降水中,中雨滴对降水率贡献始终最大,其中当雨强为10~20 mm/h时贡献率达到最大。当雨强大于20 mm/h,对降水率的贡献依次是中雨滴、大雨滴、小雨滴,大雨滴贡献率在30%左右;层状云和层积混合云中大雨滴在雨强大于1 mm/h时才开始对降水有贡献;积雨云中,当雨强大于10 mm/h时,大雨滴对降水的贡献就超过小雨滴;而层状云和层积混合云雨强大于20 mm/h时,大雨滴对降水的贡献就超过小雨滴。整体来看,中雨滴在不同雨强下对降水率的贡献起主要作用,大雨滴对雨强的贡献率逐渐增大,小雨滴的贡献率逐渐减小,当雨强大于20mm/h时,对降水起主要作用的为中、大雨滴。
Figure 2. Contribution rates of logarithmic concentration, rain intensity, and reflectivity factor of raindrops in different size bins ((a) Cumulonimbus cloud, (b) Stratiform cloud, (c) Stratocumulus mixed cloud)
图2. 不同档位雨滴对数浓度、雨强、反射率因子的贡献率((a) 积雨云,(b) 层状云,(c) 层积混合云)
Figure 3. Contribution rates of raindrops in different size bins to precipitation intensities of various magnitudes ((a) Cumulonimbus cloud, (b) Stratiform cloud, (c) Stratocumulus mixed cloud)
图3. 不同档位雨滴对不同强度雨强的贡献率((a) 积雨云,(b) 层状云,(c) 层积混合云)
3.2. 典型个例微物理量的演变特征
选取三个典型个例分析三类暴雨微物理量(粒子数浓度
、反射率因子Z、降雨强度R、含水量W、质量加权平均直径
、粒子谱宽
)随时间(世界时)的演变。
图4(a)为2022年6月26日04时至11时积雨云暴雨:各微物理参量随时间起伏剧烈,峰谷差显著。数浓度、雨强和含水量在两个数量级内波动,
在0~4 mm间变化,粒子谱宽0~6 mm,各参量趋势一致,且起伏次数多。
图4(b)为2022年6月23日04时至10时层状云暴雨:各参量变化趋势一致,但起伏幅度小。粒子数密度在一个数量级内波动,Z、R、W数值较小,Dm变化范围0~2 mm,粒子谱宽0~4 mm,参量间相关性良好。
图4(c)为2024年7月8日13时至9日02时层积混合云暴雨:各参量起伏程度介于前两者之间。2024年7月8日16:30前为层状云降雨,参量波动小(粒子数浓度102~103 m−3),质量加权平均直径
和谱宽
符合层状云特征;16:30~17:00转为积雨云降雨,参量骤升后于17:30骤降,谱宽接近6 mm (符合积雨云特征);17:00后重回层状云;7月9日01:00~02:00再次出现积雨云,变化特征相似。从
曲线看,积雨云时段
比层状云时段高1~2 mm,可初步以
作为积雨云降雨的判断依据。
3.3. 雨滴谱特征
3.3.1. 粒子谱Gamma函数拟合
图5是不同云类型粒子谱Gamma函数拟合图:由图a到c可以看出,Gamma函数拟合的小雨滴的拟合情况较好,峰值直径
均集中在0.45~0.65 mm,但对于粒子数浓度峰值的拟合偏小,大雨滴的拟合值总体偏小,且随着雨滴直径的的增加,拟合偏差越来越大。由图d可知,不同类型暴雨拟合曲线,积雨云暴雨的粒子谱宽最大,其次是层积混合云,层积云最小,三种暴雨的数浓度峰值和对应的峰值直径较为接近,且均为小雨滴,与前文分析对暴雨数浓度贡献率最大的为小雨滴的结论一致。此外,平均粒子谱Gamma函数拟合与层积混合云的Gamma函数拟合最为接近,介于层状云和层积混合云之间。在实际的天气过程中,多为以层状云或积雨云为主导的层积混合云降水,因此平均Gamma参数拟合可以作为分析暴雨过程提供,再判断暴雨主导云类型,适当向主导云类型倾斜,可以更加精确地分析尺度–数浓度关系。
(a)
(b)
(c)
Figure 4. (a) Temporal variation curves of microphysical parameters for cumulonimbus rainstorms; (b) Temporal variation curves of microphysical parameters for stratiform cloud rainstorms; (c) Temporal variation curves of microphysical parameters for stratocumulus mixed cloud rainstorms
图4. (a) 积雨云暴雨微物理参量随时间变化曲线;(b) 层状云暴雨微物理参量随时间变化曲线;(c) 层积混合云暴雨微物理参量随时间变化曲线
Figure 5. Scatter plots and corresponding gamma function fitting curves of different types of rainstorms ((a) Cumulonimbus cloud, (b) Stratiform cloud, (c) Stratocumulus mixed cloud) and comparison diagram of gamma function fitting curves among different types ((d) Gamma function fitting curves of rainstorms for different cloud types)
图5. 不同类型暴雨散点图及其Gamma函数拟合曲线((a) 积雨云,(b) 层状云,(c) 层积混合云)与不同类型Gamma函数拟合曲线对比图((d) 不同云类型暴雨Gamma函数拟合曲线)
3.3.2. Z-R关系
Z-R关系是雷达定量估测降水的基础,目前雷达定量估测降水是通过
来完成的。因雨滴谱分布差异引起的Z-R关系也随不同地区而变化,
的系数(a)和指数(b)与雨滴谱分布密切相关,按云类型分布拟合关系,能更好提高降水估测精度。
图6对不同云类型暴雨使用最小二乘法对Z-R关系拟合的曲线;可以看出,积雨云(
)、层状云(
)、层积混合云中(
)暴雨Z-R关系与标准关系
均存在差异,整体系数与标准关系式相比偏小,系数a变化范围为160.7~236.1;指数较为相似,指数b变化范围为1.33~1.42。对宿迁暴雨总样本质量控制后拟合得的平均Z-R关系为
,不同云类型暴雨和平均Z-R拟合曲线均在
下方,目前雷达估测降水所采用的标准关系式会高估淮北地区夏季暴雨降水,尤其对层状云和层积混合云暴雨误差相对较大。
Figure 6. Comparison diagram of fitting functions for Z-R relationships of rainstorms of different cloud types and the average Z-R relationship
图6. 不同云类型暴雨和平均Z-R关系拟合函数对比图
4. 结论
利用2022~024年宿迁夏季(6~8月)暴雨进行分类和统计处理,得到了积雨云、层状云、层级混合云等3类不同类型暴雨的微物理参量、Gamma函数分布以及Z-R关系等,并进行了相应的分析,结果如下:
1. 雨滴下落速度与经验公式差异显著:宿迁不同云型暴雨雨滴落速拟合曲线(
)拟合效果良好,而Atlas经验公式曲线位于拟合曲线下方,且随雨滴直径增大(尤其D > 1 mm))间距扩大,按该经验公式计算会低估降水强度,且中、大雨滴浓度贡献率越高,低估程度越大,不同云型暴雨需用对应落速函数计算雨强。
2. 微物理参量具云型差异性:数浓度(
)、平均雨强(R)、含水量(W)及多数特征直径(
、
、
)均呈“积雨云 > 层积混合云 > 层状云”规律,仅峰值直径(
)为积雨云与层积混合云相同且大于层状云,数浓度(
)范围381~613 m−3,平均雨强(R)范围1.96~5.72 mm/h,雨水含量(W)范围0.12~0.29 g/m3。
3. 不同尺寸雨滴贡献各有侧重:数浓度上,小雨滴(
)占比最高(层状云87.3%,积雨云、层积混合云均80.8%);雨强上,中雨滴(
)是主要贡献者(贡献率60.4%~67.9%);雷达反射率因子(Z)上,积雨云以大雨滴(
)为主(65%),层状云、层积混合云以中雨滴为主(分别64.5%、51.6%),且大雨滴因Z与D6成正比,虽数量占比低但对Z贡献大。
4. 微物理参量起伏相似且可用
初判云型:暴雨类型决定了其微物理参量的时间演变特征,积雨云暴雨参量起伏最剧烈、峰谷差最显著且波动次数多,层状云暴雨参量起伏幅度最小且参量间相关性良好,层积混合云暴雨参量起伏程度介于两者之间,且会随云型(层状云与积雨云)转换呈现“低波动–骤升骤降–低波动”的阶段性变化,同时质量加权平均直径
可作为初步判断云型的依据。
5. 雨滴谱Gamma函数拟合存在特点:Gamma函数对小雨滴拟合较好,峰值直径
集中在0.45~0.65 mm,但对粒子数浓度峰值及大雨滴拟合值偏小,且雨滴直径越大拟合偏差越大;谱宽呈“积雨云 > 层积混合云 > 层状云”,三者数浓度峰值及对应
接近且均为小雨滴;平均Gamma参数拟合与层积混合云最接近,可作为暴雨分析参考,实际需结合主导云类型调整。
6. Z-R关系与标准式差异明显且具云型依赖性:宿迁暴雨Z-R关系整体系数a在160.7到236.1之间,小于标准式
,指数b在1.33到1.42之间与标准式相近;总样本拟合关系为
,标准式会高估宿迁暴雨降水,尤其层状云和层积混合云误差较大,分云型拟合Z-R关系可提高降水估测精度。
本论文建立的宿迁分云型Z-R关系能优化传统标准式,更精准反演不同云型暴雨雨强,提升本地雷达定量降水估测精度,同时明确的不同云型暴雨微物理特征可助力通过雷达回波快速判断暴雨云型、预判雨强变化,为暴雨灾害预警与防控提供支撑,一定程度上补充江淮北部气候过渡带宿迁的暴雨微物理基础数据;但本论文也存在一定局限性,在仅依托2022~2024年夏季9次暴雨样本,样本量少、时间跨度短,代表性受限,无法充分反映宿迁夏季暴雨微物理特征的年际变化规律,可能导致研究得出的雨滴谱特征结论与分云型Z-R关系的普适性受限,无法完全代表宿迁长期夏季暴雨的真实微物理状况。