雷暴云起电机制及其数值模拟的回顾与进展

对国内外有关雷暴云起电机制和起/放电模式的研究结果和进展进行了回顾,总结了这两方面研究的主要进展及所涉及的重要问题。详细介绍了离子、粒子及其他起电机制的发展历程,并从非云模式和云模式这两方面综述了国内外起/放电模式的一些结果和进展。目前对雷暴云内电荷的产生已有较为系统地了解,但仍需对雷暴内动力、微物理和起电过程及它们之间的相关性做更深入的研究。进一步完善雷暴云起/放电模式,并应用模拟结果和观测资料来研究动力、微物理过程及电结构的相互作用机制。     

基于高分辨率中尺度气象模式的实际雷暴过程的数值模拟试验

雷暴数值预报的实际应用离我们还有多远？本文对此进行了尝试, 即利用一个复杂的高分辨率中尺度气象模式驱动一个三维雷电模式, 在只采用常规气象观测资料的条件下, 对北京的一次实际雷暴过程进行模拟试验, 分析了雷暴云的宏观动力、 微物理过程及电结构的时空变化特征以及其可能的相互作用机制。结果表明: 利用高分辨中尺度模式预报出的三维气象场作为雷电模式的初始场, 完全可以不需添加虚假的扰动来触发雷暴云的发展, 高分辨中尺度模式的预报场本身所包含的水平非均匀、 垂直强非静力性及较强的对流不稳定信息足够促发雷暴云的剧烈发展; 用较为真实的三维气象场作为初始场模拟产生的电场分布特征与云微物理分布特征及环境气象要素的分布结构非常协调, 得到的雷暴云的电荷结构特征以及电结构随时间的演变特征更为复杂, 更真实的体现了实际雷暴云本身发展的复杂性, 同时, 模式能够模拟出合理的云闪及正负云地闪, 且模拟的闪电频数随时间发展演变趋势基本与观测实况基本吻合, 从而表现了对雷电天气潜在的预报能力。本次模拟的北京雷暴云在发展过程中, 水物质霰的最大质量比、 最大正电场强度及闪电频数随模拟时间的演变发展趋势非常相似。     

播撒金属丝对雷暴云电结构影响模拟研究

利用积云动力一起电二维轴对称模式,模拟分析了雷暴不同发展阶段在云中上升气流速度最大区播撒金属丝对雷暴云电结构的总体影响。通过与未播撒模拟结果相比较,发现各种降水粒子的荷电率及电场强度都有所减小,这种减小在雷暴发展的成熟阶段尤为明显,雷暴云内的扩散、电导和感应起电机制作用减弱。通过不同播撒位置和播撒率的敏感性试验得到,在云中强上升所泫区播撒对垂直电场的抑制效果比在云底和云顶更好,不同的播撒率引起的电     

雷暴云不同空间电荷结构数值模拟研究

利用一个三维时变双参数动力电耦合模式对决定雷暴云空间电荷结构的三个初始场：反转温度、中心最大扰动位温和扰动区域进行模拟对比研究。结果表明,选取不同的参数值,雷暴云电结构有很大差异,从而解释了在不同地区、季节、强度的雷暴云中所观测到得不同的空间电荷结构的原因。     

雷暴微物理过程对电活动影响的数值模拟研究

在一维时变雷暴云起电放电模式的基础上,分别利用对流降水试验合作计划(CCOPE)和夏季雷暴放电与降水研究(STEPS)一次雷暴的实测探空资料进行了初步模拟,并通过敏感性试验研究了软雹垂直通量、地面温度和反转温度对雷暴云起电放电过程的影响.模式考虑的冰晶生成机制为Fletcher 和H-M冰晶生成机制.对比CCOPE和STEPS雷暴个例的模拟与观测结果,发现垂直电场正负极大值出现的位置与实测结果较接近,但电场极值差异较大,模拟结果明显大于实测结果.敏感性试验表明,当软雹垂直通量增加时,雷暴云首次放电时间降低,而单位时间内发生闪电的次数(闪频)随之增大.当地面温度较高时,雷暴云中负电荷的垂直厚度较小,雷暴云的正、负电荷最大密度较大.当反转温度较低时,冰晶所带负电荷的垂直范围较大.在敏感性试验中,当模式中地面温度为32℃,反转温度＜-26℃时,雷暴云出现反极性电荷结构.该模式能较直观地模拟雷暴发展过程中的电荷结构和电场时空分布等.     

江淮地区副高后部强对流的分析和数值模拟研究

选用1989年7月14日20时南京站的实测探空为初始资料,应用许焕斌的二维可压缩大气云尺度模式,对副热带高压（以下简称副高）后部的强雷暴产生及发展进行了数值模拟和分析。模拟结果与实际雷暴的出现情况较一致,并通过不同初始扰动方案模拟结果的对比分析,表明扰动方案1对这次过程的模拟效果最好。     

海气耦合随机动力模式的Fokker—Planck方程

将描写海温和气温交互作用的随机动力模式化为一个Ｆｏｃｋｋｅｒ－Ｐｌａｎｃｋ方程（ＦＰＥ），然后用矩阵连分法进行求解，并对二氧化碳增温效应进行了计算。当ＣＯ２从３３０ｐｐｍ增加到６６０ｐｐｍ时增温１．２℃，模式中存在的优势周期为３￣４年，当ＣＯ２倍增时各周期有延长的趋势。     

中低层水平风速对闪电和降水影响的数值模拟

利用三维雷暴云动力电耦合模式,对无风、弱风、中等强度风和较强风这四种不同中低层水平风速情况下雷暴的发展特征进行了敏感性试验。模拟结果表明,随着中低层水平风速增大,雷暴云内粒子生长受到抑制,云内起电过程将变得缓慢,起电区域发生转移,更长的起电时间使云内产生更大的电场、发生更剧烈的放电过程,但若风速过大,云内电场将被严重抑制而不发生放电,中等强度风最有利于闪电发生。中低层水平风速在一定范围内增长将使降水增多,但若风速过大,降水将减少。中低层水平风对闪电和降水影响明显,在对雷暴云闪电和降水的研究中应当综合考虑其对云内动力、微物理过程的影响。     

CISK、蒸发-风反馈机制和热带大气30-60天振荡

建立一个包含ＣＩＳＫ机制和蒸发－风反馈机制（Ｅｖａｐｏｒａｔｉｏｎ－ＷｉｎｄＦｅｅｄｂａｃｋＭｅｃｈａｎｉｓｍ，简写为：ＥＷＦＭ）的热带大气运动的无耗散、线性斜压半地转模式，引入反映ＣＩＳＫ的无量纲凝结潜热参数η和反映ＥＷＦＭ的无量纲参数α，并求得模式系统的解析解。动力学分析和模式大气的计算结果进一步说明了ＣＩＳＫ和ＥＷＦＭ的共同作用是激发和驱动热带大气低频振荡的重要物理过程。     

雷暴电过程对动力发展的影响研究

利用三维时变积云动力—电过程耦合模式对雷暴中的电过程对动力发展的影响进行了模拟计算。结果表明,随着雷暴云体发展,云内形成三极性电荷结构。电场力对雷暴发展的直接作用与动力场相比很小,电涡度10-6·s-1比动力涡度小两个量级。但电场力影响了雷暴中的微物理过程,特别是水成物粒子的降落速度,并通过对三相水成物之间转换过程的调制,改变了三相水成物的时空结构,使云中水汽总量增加41%,释放潜热增加19.4%,从而改变了流场结构,电涡度达到10-4·s-1量级,和动力涡度相当。模拟结果指出,雷暴内电过程在云体成熟阶段对雷暴发展有较大影响。     

Investigations of the development of thunderstorm with hail. Part 3. Numerical simulation of cloud evolution

The three-dimensional nonstationary model of a convective cloud is used for investigating a thunderstorm with hail which developed over Pyatigorsk on May 29, 2012 and produced a severe hailstorm. The values of cloud characteristics (liquid water content, ice content, vertical velocity, etc.) are obtained. The importance ofconsidering wind shear is noted. The simulation results are used to analyze the transformation of precipitation field and the electric charge structure of the analyzed cloud during its development.     

Contribution of the ground corona ions to the convective charging mechanism

The convective charging mechanism of thunderclouds is based on the vertical transport of space charge generated by corona from ground irregularities under the influence of the surface electric field. The present work estimates the amount of charge which is expected to reach cloud base by conduction and convection processes during the lifetime of a thunderstorm. This estimate is made using the numerical model PICASSO, previously designed to characterize the evolution of this corona space charge between ground level and cloud base. Experimentally determined values of surface electric field are introduced into the model in order to initiate the computation. These values are based on six events documented during four different field experiments carried out in Florida and in France. As an upper limit, the convective transport is uniformly applied as a linearly increasing vertical air speed profile and competes with the conductive transport. The fraction of the positive charge generated at the surface by corona which finally reaches the upper limit of the layer varies between 26 and 86%, essentially depending on the electric field evolution at altitude. Assuming that the vertical transport conditions remain the same over an area of 10 km×10 km, the overall charge amount can be roughly estimated. It ranges between about 63 and over 300 C. Because the present assumptions probably lead to an overestimate of these amounts, such a range suggests that the convective charging mechanism is unlikely to be able to account for the major electrification process of the thundercloud. However, it could be considered as a relevant mechanism contributing to the lower positive charge center of the thundercloud, often observed close to cloud base.     

初值中云变量对T213预报性能的影响

采用云变量自由变化的方案，在国家气象中心全球业务模式T213L31的初始场中增加有关云变量的信息，通过2005年6—8月和2015年12月—2006年2月各3个月的连续滚动对比试验的统计分析和个例预报分析，研究探讨了全球模式初值中增加云变量对模式预报性能的影响。初步研究结果表明:采用自由变化的方案在初值场中增加云的信息，使模式能够较为合理地描述出模式预报初期与云相关变量分布和变化特征，降低了spin-up现象对模式前期降水预报能力的影响，同时对500 hPa形势场预报也有一定程度提高。     

基于LAPS资料的一次冰雹过程数值催化模拟研究

基于3维冰雹云模式,采用局地分析预报系统高时空分辨率分析场作为初始场,对鄂西北地区2009年4月15日发生的一次降雹天气过程进行数值模拟,分析了强对流发展过程的流场、雷达回波、水成物粒子演变特征以及催化机制。研究表明:基于高时空分辨率资料的3维冰雹云模式,能够较好地模拟冰雹云的形成发展演变及人工催化后微物理过程响应。模式输出的雷达回波强度及回波顶高与雷达实测资料相近,最大反射率因子分别为70 dBz、65 dBz,强回波中心基本上位于6~7 km处。冰雹云发展过程中呈现典型的低层辐合、高层辐散特征,消亡过程流场相反;冰雹胚胎形成主要源于冻滴自动转化,碰并云水和雨水促进冰雹增长;不同时间、不同剂量的催化对于防雹效果差异显著,催化时间越早,催化剂量越大,效果越好。     

Fine-resolution simulation of cloud-to-ground lightning and thundercloud charge transfer

2-D 12.5 m-resolution simulations of cloud-to-ground (CG) lightning discharge processes have been performed using an improved stochastic lightning model for different types of cloud charge distributions, such as dipole, tripole, bi-dipole and multi-layer charge structures produced from the numerical simulation of thundercloud electrification. The modelling produced the fine branched channel structure of CG lightning and the results illustrate the relations between CG lightning channel propagation and cloud charge distribution. The simulated features of CG lightning are associated with the observed results. The simulation studies are essential in our understanding of complex charge transfer processes caused by CG lightning discharges in thunderclouds. The induced charges of opposite polarity are deposited or embedded in the local volumes where the bidirectional leaders passed during a CG lightning discharge. Although the embedding affects charge structure only in a pair of significant positive and negative charge regions closest to the ground, the electric field strength acutely weakens and electrostatic energy in thunderclouds is significantly consumed when the discharge terminates. In addition to simulating the upward and downward breakdown of the initial leader to ground and the ensuing return stroke (RS), the simulation assumes that current continues to flow in the channel to ground and determines the upward breakdown until the end of the discharge. For the subsequent discharge sub-process, the upward leader channel tends to transfer the charges with the same polarity as the RS, while the downward leader channel favors transfer of opposite charge to ground. In the sub-processes of a few CG flash simulations, the magnitude of the opposite charges from the downward leader exceeds that of charges with the same polarity from the upward leader so that the net charges transferred to the ground have a reversed polarity to the RS. The simulation presents similar features of CG lightning as those observed in realistic bipolar CG lightning.     

Numerical simulation of the relationship between electrification and microphysics in the prelightning stage of thunderstorms

To further investigate the influence of cloud base temperature, updraft velocity and precipitation particle constitution on cloud electrification, five thunderstorms in various regions of China were simulated by using the three-dimensional compressible hailstorm numerical model including inductive and non-inductive charging mechanisms. The results indicate that changes of cloud base temperature have an influence on the initial electrification. Comparison of the above cases shows that in the case of warm cloud base and moderate updraft velocity (< 20 m s^{−  1}), active electrification occurred below the − 10 °C level before moving upward to the − 20 °C level. In contrast, when cloud base is cold and updraft velocity is intensive, the main charging region is at the − 20 °C or even higher level. In that case, the vertical extent of the main negative charge region becomes larger with the increase of cloud base temperature. Apart from the main dipolar or tripolar charge structure, some smaller charge regions with relatively high values of charge density may also appear. Frozen drops, originating mainly from supercooled raindrops, mainly get electrified through charging interactions with snow at or below the − 20 °C level. They are responsible for the negative charge region near the melting level at the initial stage of precipitation if there is a large supercooled raindrop content. Non-inductive charging during hail-snow collisions is rather weak, resulting in the charge density on hail of no more than − 0.01 nC m^− 3.     

Aircraft observations of electrical conductivity in warm clouds

Aircraft observations of electrical conductivity and cloud microphsical, dynamical and other electrical parameters were made in warm stratocumulus and cumulus clouds forming during the summer monsoon seasons (June-September) of 1983 and 1985 in the Deccan Plateau region, India. A Gerdien type cylindrical condenser was used for the measurement of electrical conductivity. The variations in the electrical conductivity are observed to be closely associated with the updrafts and downdrafts in the cloud, liquid water content, cloud droplet charge and coro-na discharge current. The value of electrical conductivity in warm clouds is found to be in the order of 10-12 ohm-1 m-1 which is two orders higher than that observed in clear-air at cloud-base levels in some regions by other investigators.Classical static electricity concepts predict reduced conductivity values inside clouds. Cloud electrical conductivi-ty measurements, particularly in warm clouds are few and the results are contradictory. The recently identified mech-anism of vertical mixing in clouds lends support to coovective charge separation mechanism with inherent larger than clear-air values for cloud electrical conductivity and therefore consistent with the measurements reported herein.     

Shapley-模糊神经网络方法在华南台风卫星云图的长时效滚动预测中的应用

为了更好地利用大量的卫星云图观测资料来提高台风暴雨的预报能力,解决并提高对台风强降水云系变化的预报精度,延长对未来云系变化的预报时效,构建基于合作对策Shapley-模糊神经网络的华南区域台风卫星云图非线性智能计算滚动集合预测模型,对增强卫星云图资料在台风暴雨天气预报中的实用性和及时性具有重要意义.依据2013—201...     

雷暴云宏观特性对其电过程影响的数值模拟研究

利用一维雷暴云起电放电模式,初步模拟了STEPS(夏季雷暴降水与闪电研究计划)一次雷暴个例的一些基本电学特性。模拟的电荷结构为:雷暴云上部为正电荷,下部为负电荷,即符合典型雷暴云的偶极性电场分布特征,放电层电场随时间变化规律呈锯齿状分布。在此基础上,对云底高度和温度垂直递减率进行了敏感性试验,研究这两种因子对雷暴电过程的影响。结果表明:当云底高度降低时,放电高度升高,闪电频数增加,首次放电时间随之提前,即放电过程变强。温度垂直递减率增大,闪电频数降低,首次放电时间随之延后,正负电荷的分布范围减小。     

闪电起始过程时空特征的宽带干涉仪三维观测

闪电的起始位置和起始阶段发展速度是闪电研究中的重要问题。2010年夏季,使用架设在广州市从化区的两套甚高频 (VHF) 宽带干涉仪对闪电的起始阶段放电过程进行三维定位观测。对观测数据给出的地闪和云闪的起始高度分布特征以及起始阶段击穿过程的时空发展特征进行统计和对比分析,结果表明:闪电的起始高度分布呈双峰值特征,分别在5.0 km和8.8 km有两个明显的分布峰值,符合雷暴云三极性总体电荷结构的描述。对起始阶段闪电放电发展速度的计算表明,云闪和地闪在起始阶段的前15 ms内的平均发展速度均在104~105 m·s-1量级之间;多数云闪、地闪起始阶段前15 ms内的平均发展速度表现出减速趋势,但云闪个例中起始阶段前10 ms存在减速趋势的比例更高,且其中在前15 ms一直保持减速趋势个例所占比例也大于地闪。云闪和地闪的起始阶段放电过程的发展方向有向上、向下和水平发展3种情形,可用于指示闪电始发位置的环境电场方向。