气溶胶对雷暴云微物理过程和起电影响的数值模拟

将云滴冻结方案植入已有的二维雷暴云起、放电模式,结合一次山地雷暴个例,探讨了气溶胶浓度对雷暴云微物理过程、起电以及空间电荷结构的影响。结果表明:气溶胶浓度增加,云滴数目增多,尺度降低,雨滴含量减少;云滴冻结导致冰晶在低温区快速生长,冰晶数浓度增加,尺度减小,当气溶胶浓度高于1000 cm-3后小冰晶难以增长成大尺度的霰粒子,因此霰粒子数浓度先增加后急剧减少。此外,气溶胶浓度的大小不会影响雷暴云的电荷结构特征,但会对云内的起电强度产生明显的作用:当气溶胶浓度较低时,增加气溶胶浓度,更多的冰晶和霰粒子发生碰撞使得云内起电过程增强,空间电荷密度增加;当气溶胶浓度高于1000 cm-3后,少量的霰粒子和小冰晶的出现抑制了非感应起电过程,导致电荷密度降低。     

重庆地区一次雷暴云电过程及其对初始云滴数浓度响应的数值模拟

采用耦合了Saunders和Takahashi两种非感应起电参数化方案的RAMS(Regional Atmospheric Modeling System)模式,对重庆地区一次雷暴过程进行模拟,对比分析了两种起电参数化方案下,电荷开始分离时和雷暴云发展到成熟阶段时的水成物粒子的分布、所带电荷密度以及雷暴云的电荷结构分布。模拟结果表明,在Saunders起电参数化方案下,雷暴云的电荷结构从起电到放电都呈现偶极性特征,而在Takahashi参数化方案下,雷暴云的电荷结构则由反偶极性发展成正偶极性。为研究CCN(cloud condensation nuclei)对雷暴云的影响,本文进行了两组敏感性试验,随着云滴初始数浓度增加,雷暴云的电荷结构没有发生极性翻转,但雷暴云中电荷量增加,电荷分布区域变大,有利于闪电发生。在Saunders起电参数化方案下,当云滴初始数浓度大于2 000 cm-3时,电荷量变小。通过分析微物理量场和微物理过程发现,随着云滴初始数浓度增加,冰相粒子质量混合比增加,在Saunders起电参数化方案下,当云滴初始数浓度大于2 000 cm-3时,霰粒子质量混合比减小。验证了CCN的变化能影响云的微物理过程,从而影响雷暴云的电荷分布以及闪电的发生,尤其是冰相物质的变化显著影响了雷暴云的起电过程。     

冰晶繁生对雷暴云非感应起电影响的数值模拟

为了研究冰晶繁生在雷暴云发展过程中对非感应起电过程的影响,利用三维雷暴云模式在理想层结环境下,对雷暴云内各种水成物粒子、电荷以及电场分布情况进行数值模拟。模拟结果表明:在雷暴云发展和成熟阶段,有繁生过程参与的雷暴云中下部存在一个冰晶聚集区域,从而使得云内冰晶的数量较无繁生过程增大约1 5%～18%,且聚集的区域范围更大;同时,繁生过程的加入也使得霰粒子数量也比无繁生过程时增大约20%;霰冰非感应电荷转移的正区一般位于霰粒子浓度高值区附近,而负区位于冰晶和霰粒子浓度高值区相重合的区域;冰晶繁生过程通过影响雷暴云中冰晶和霰粒子浓度和分布位置,使得雷暴云非感应起电的强度和位置发生改变,导致云内起电过程提前约5～6 min。     

大气冰核对雷暴云电过程影响的数值模拟

利用已有的二维雷暴云起、放电模式模拟了一次雷暴天气,并通过敏感性试验研究了冰核浓度变化对雷暴云动力、微物理及电过程的影响。结果表明：随着大气冰核浓度的增加,雷暴云发展提前,上升气流速度和下沉气流速度均呈现降低的趋势。大气冰核浓度提升有利于异质核化过程增强,冰晶在高温区大量生成,而同质核化过程被抑制,因此冰晶整体含量降低,引起低温区中霰粒含量降低和高温区中霰粒尺度降低。在非感应起电过程中,正极性非感应起电率逐渐减小,负极性非感应起电率逐渐增大。由于液态水含量随大气冰核浓度的增加逐渐降低,高温度冰晶携带电荷的极性由负转变为正的时间有所提前。在感应起电过程中,由于霰粒尺度减小及云滴的快速消耗,感应起电率极值逐渐降低。冰晶优先在高温区生成而带负电,不同大气冰核浓度下的雷暴云空间电荷结构在雷暴云发展初期均呈现负的偶极性电荷结构。在雷暴云旺盛期,随着冰核浓度增加,空间电荷结构由三极性转变为复杂四极性。在雷暴云消散阶段不同个例均呈现偶极性电荷结构,且随着冰核浓度的增加电荷密度值逐渐减小。     

湖北地区一次雷暴云电荷结构和放电特征的数值模拟研究

在三维强风暴动力—电耦合数值模式中引入非感应起电参数化方案、感应起电参数化方案以及放电参数化方案,对湖北宜昌2014年6月19日一次闪电过程中雷暴云电荷结构和放电特征进行了模拟分析。模拟结果表明,当云内粒子增多、增大,大部分霰粒子逐渐降落到中低层,上部正电荷区减小,底部正电荷堆范围开始扩大,中部负电荷区和底部正电荷区成为主要的起电区域,这种底部正电荷区较厚的三极性电荷结构不利于地闪的产生。在粒子带电分析中,霰与冰晶粒子携带的电荷量均大于云滴,说明霰与冰晶之间非感应碰撞是云中主要的起电过程。虽然云滴的电荷量较小,但霰与云滴之间感应碰撞的作用不可忽视。结合电荷结构的分布,发现底部正电荷堆的垂直分布高度与霰粒子、云滴的电荷浓度的分布有关,且霰与云滴电荷浓度的累积区与底部正电荷堆相一致。     

气溶胶对雷暴云起电以及闪电发生率影响的数值模拟

本文利用二维耦合气溶胶模块的雷暴云起电模式,结合一次南京雷暴个例,进行250 m分辨率雷暴云起电模拟实验,探讨了气溶胶浓度对雷暴云空间电荷分布以及闪电发生率的影响。在这个气溶胶模块中,假定一个三模态的气溶胶对数分布,考虑了气溶胶活化过程。结果显示:(1)随着气溶胶浓度增大,雷暴云电荷结构保持为三极型。(2)当气溶胶浓度从50 cm-3增加至1000 cm-3时,水成物粒子浓度上升,雷暴云电荷量和闪电发生率增加明显。(3)气溶胶浓度在1000~3000 cm-3范围时,云水竞争限制了冰晶的增长,导致雷暴云上部主正电荷堆电荷量降低。云滴和霰粒子浓度缓慢上升促进中部主负电荷堆和底部次正电荷堆电荷量继续增大。闪电发生率保持稳定。(4)当气溶胶浓度大于3000 cm-3时,水成物粒子浓度稳定,云内的电荷量以及闪电发生率保持为一定量级。     

冰晶异质核化对雷暴云电过程影响的数值模拟

引入一种新型冰晶异质核化方案, 基于二维雷暴云模式, 探讨雷暴云电过程对三种异质核化的响应。结果表明: 浸润核化是冰晶生成的最重要异质核化过程, 较高数浓度的冰晶消耗雷暴云内液态水含量, 抑制淞附过程, 导致霰粒子比含水量低, 表现为较强的负极性非感应起电率; 接触核化生成的冰晶量最少, 仅对雷暴云中下层3~5 km处的冰晶有贡献, 同时霰粒子数浓度较低, 导致该方案下的起电过程最弱; 沉积核化主要影响云砧处的冰晶, 有利于提高霰收集云滴的效率, 表现为极高的霰比含水量, 促进低温区非感应起电过程的发生。总体上来看, 三个方案下的电荷结构均由较复杂的多极性发展为偶极性。其中浸润方案中主正电荷区的抬升最明显, 而接触方案过低的冰晶分布高度与沉积方案过高的冰晶分布高度, 都直接导致了次正电荷区更快消散。      

广东一次飑线过程中一个雷暴单体成熟阶段的电荷结构演变特征的数值模拟

为了进一步印证以往观测反推得到的广东地区雷暴云多偶极性电荷结构的结论,利用加入了起放电参数化方案的WRF模式,模拟了广东在2017年5月8日发生的一次飑线过程,并对这次飑线过程中一个雷暴单体成熟期的电荷结构演变特征进行分析,通过分析动力、云水含量、各水成物粒子混合比及携带电荷情况,讨论了电荷结构的形成及演变机制。结果表明,成熟阶段的单体,电荷结构从三极性逐渐演变为偶极性。这是因为在成熟初期,霰粒子在有效液态水含量适中且温度较高的地方与冰晶/雪花粒子发生了非感应碰撞,因此底部霰粒子携带正电,雷暴云底部形成次正电荷区,电荷结构为三极性。而在成熟后期,由于丰富的云水含量,使冰粒子的凇附过程增强,霰不断增加,冰晶和雪花不断被消耗,温度较暖区域与霰共存的冰晶和雪花急剧减少,使得该区域大小冰粒子的非感应碰撞起电急剧减少,此处霰粒子不能再通过非感应碰撞获得正电荷,底部次正电荷区随之消失,雷暴云的电荷结构转变为偶极性。此结果和以往观测反推得到的结论不同,这表明,对南方雷暴电荷结构还需继续深入认识。      

大气冰核谱对雷暴云微物理过程及起电影响的数值模拟

利用已有的二维对流云模式,讨论了三种不同的冰核谱对雷暴云微物理、起电及电荷结构的影响。模拟结果表明:(1)不同的冰核谱环境对雷暴云中冰相粒子的含量及分布具有明显作用。冰核谱的垂直温区越大,产生的冰相粒子分布越广。在冰核浓度较大的个例中,冰晶和霰粒子的含量高,更多的小冰相粒子出现在海拔更高的区域;(2)高温区冰核的数量会对上升气流速度产生显著影响。高温区的冰核越多,冰相粒子在微物理发展过程中释放的潜热越多,上升气流强,对流发展越旺盛;(3)在低温区冰晶浓度高的谱环境个例中,雷暴云中的非感应起电率和感应起电率高,导致起电量增加。高温区冰核多的谱环境,大量冰晶和霰获得正电荷形成次正电荷区,电荷结构呈现三极性;而高温区冰核少的谱环境,参与起电的水成物粒子少,易形成偶极性电荷结构。     

不同水汽条件下气溶胶对雷暴云电过程影响的数值模拟研究

为全面了解水汽在气溶胶影响雷暴云电过程中的作用,本研究在已有的二维雷暴云起、放电模式基础上,通过改变相对湿度和气溶胶初始浓度（文中气溶胶浓度均指气溶胶数浓度）进行敏感性数值模拟试验。结果表明:（1）随着气溶胶浓度升高,雷暴云产生更多的小云滴,降水过程受到抑制。而当水汽含量升高时,云滴数浓度的增长速度更快,雨滴数浓度升高,缓解了降水变弱的趋势。（2）水汽含量较低时,随着气溶胶浓度升高,更多小云滴被带入冻结层形成大量小冰晶,霰粒含量升高,雷暴云起电过程增强。气溶胶浓度升高至一定的量级（3000 cm?3）时,冰晶尺度减小和雨滴浓度降低抑制霰粒生长,雷暴云起电过程受到削弱。感应起电和非感应起电过程随气溶胶浓度升高呈先增强后减弱的趋势。水汽含量的升高促进了冰相粒子的增长,起电过程呈现持续增强的趋势,气溶胶浓度为3000 cm?3时起电率达到极值,电荷密度的增幅扩大。（3）水汽含量较低时,雷暴云难以发展成深厚的系统,气溶胶浓度变化对其影响不明显,电荷结构由三极性发展,在消散期演变为偶极性电荷结构;水汽含量较高时,雷暴云迅速发展成深厚的系统,随着气溶胶浓度升高,在雷暴发展旺盛阶段电荷分布表现为多层复杂结构。研究显示水汽含量在气溶胶浓度变化对雷暴云微物理、起电过程及电荷结构的作用中扮演重要角色。      

对流强度对雷暴云微物理发展和电荷结构影响的数值模拟

为了探讨对流强度大小对雷暴云内微物理发展和起电过程的影响,基于已有的二维积云起、放电模式,改变其扰动温度进行敏感性试验。试验结果表明:对流强度对雷暴云内微物理过程、起电率及后续电荷结构的产生均有一定程度的影响:1）对流强度较小时,冰晶粒子极大值在高温区（高于-13.8℃）出现,对流强度较大时,上升风明显增强,将更多的水汽带入高空,气溶胶活化过程明显增强,使得云滴粒子明显增多,冰晶粒子较早产生,冰晶粒子极大值在低温区（低于-13.8℃）出现,发展过程更为剧烈;同时,较高的对流强度也使得降雨量增多,霰粒子数目也在对流发展旺盛时期显著增多。2）非感应起电率主要和冰晶-霰的碰并分离过程有关,对流强度较大时,非感应起电率较大,电荷结构持续时间较长,过程明显,感应起电率也较强。3）对流强度较大时,电荷结构更为复杂,雷暴云发展初期基本呈现为三极性,发展旺盛时期底部正电荷区域嵌入一个较小的负电荷区,呈现四极性电荷结构,雷暴云发展末期基本呈现偶极性电荷结构;对流强度较小时,发展初期、旺盛时期均呈现三极性电荷结构,发展末期呈现偶极性电荷结构。     

一次暖云降水主导的暴雨过程中电荷结构特征及其成因的模拟研究

为了认识以暖云强降水为主导的对流单体中的电荷结构特征及其形成原因, 利用加入了起放电参数化方案的WRF模式, 模拟了2017年5月7日广州局地突发的以暖云降水为主导的特大暴雨过程, 分析讨论了此次过程中一个单体成熟发展阶段的电荷结构的特征及其成因。结果表明, 此次以暖云降水为主导的特大暴雨过程中的单体对流强度较弱, 云顶高度低于同地区典型对流过程, 强回波区由大雨滴形成, 范围较小, 顶较低, 对流运动向0℃层以上输送的过冷水较少, 不利于冰相粒子形成, 导致大小冰相粒子含量均较少, 其中含量最多的冰相粒子为雪花, 其次依次为霰、冰晶、冰雹。云内起电较弱, 以非感应起电为主。非感应起电主要以对流区中-15℃层以下正的起电率为主, 感应起电率以对流区中的负极性为主。对流区中空间净电荷呈三极性结构, 其中中部负电荷区和底部正电荷区中心电荷密度及电荷区范围相当, 上部正电荷区相对较弱, 范围较小。对流区外围仅有弱的中部负电荷区和底部正电荷区。中部负电荷区由带负电荷的冰晶和雪花共同主导, 上部正电荷区由带正电荷的雪花主导, 底部正电荷区主要是由带正电荷的霰粒子及带正电荷的雨滴主导。强起电区和放电区重合, 主要集中在回波中心上部35~50 dBZ的对流区。      

CCN数浓度及参数化方案对一次层云降水过程影响的数值研究

利用中尺度WRF模式(V3.7),采用WDM6双参数微物理方案,对2014年7月26日12时—28日06时发生在华东地区的一次层状云降水天气进行数值模拟。通过改变模式中初始云凝结核(CCN)数浓度及参数化方案,进行敏感试验,对模拟结果进行对比分析。改变CCN数浓度的结果表明,CCN数浓度对降水的影响复杂、非线性,随着CCN数浓度的增大,降水量减小。云水、霰混合比始终增加,雨水混合比表现为先增加后减小再增加的趋势,冰晶混合比则与之相反,呈现先减小再增加再减小的趋势,雪晶混合比呈现先减小后增加的趋势;改变CCN参数化方案的结果表明,两者模拟降水落区有差别,三参数方案更接近实际;降水产生后,三参数方案的CCN浓度一直高于双参数方案,且数值变化不大;双参数方案的结果显示暖云降水加强,冷云降水略弱,三参数方案则显示暖云降水较弱,冷云降水较强。     

Idealized simulation of the Colorado hailstorm case: comparison of bulk and detailed microphysics

One of the purposes of the Fourth Cloud Modeling Workshop was to compare different microphysical treatments. In this paper, the results of a widely used bulk treatment and five versions of a detailed microphysical model are presented. Sensitivity analysis was made to investigate the effect of bulk parametrization, ice initiation technique, CCN concentration and collision efficiency of rimed ice crystal–drop collision. The results show that: (i) The mixing ratios of different species of hydrometeors calculated by bulk and one of the detailed models show some similarity. However, the processes of hail/graupel formation are different in the bulk and the detailed models. (ii) Using different ice initiation in the detailed models' different processes became important in the hail and graupel formation. (iii) In the case of higher CCN concentration, the mixing ratio of liquid water, hail and graupel were more sensitive to the value of collision efficiency of rimed ice crystal–drop collision. (iv) The Bergeron–Findeisen process does not work in the updraft core of a convective cloud. The vapor content was always over water saturation; moreover, the supersaturation gradually increased after the appearance of precipitation ice particles.