气溶胶是指悬浮在气体中的固体和液体微粒与气体载体共同组成的多相体系，气溶胶可以通过吸收和散射太阳辐射直接影响地—气系统的辐射平衡，也可以作为云凝结核和冰核影响云微物理过程。目前，人们对气溶胶—云滴的相互作用了解较多，但对气溶胶—冰晶的相互作用了解较少，冰晶粒子由可溶性气溶胶液滴同质核化形成或在不可溶气溶胶粒子(冰核)帮助下异质核化形成，异质核化过程有多种机制，如凝华、碰触、浸润等，多种气溶胶可以充当冰核，如沙尘、黑碳、花粉等(Blyth and Latham，1993；王春明等，1997；Yin et al.,2000，2002；Murray et al.,2010；房文和郑国光，2011)。大量观测证明气溶胶浓度改变对冰晶含量及冰晶尺度均有影响，其中气溶胶对冰晶数浓度影响最为明显，气溶胶数浓度增加导致云中过冷水量增多，冰晶粒子增加，霰粒子浓度降低(van den Heever et al.,2006；Fan et al.,2007；Yang et al.,2011)；Khain(2009)使用二维云模式证明了海洋上空气溶胶的增多会增加雷暴云中过冷水滴含量、冰晶含量以及上升风速；在沙尘天气系统下，许多学者认为沙尘气溶胶作为有效冰核，在低于－8°C时通过接触冻结核化、低于－15°C时凝华核化，在－9～－30°C时依靠表面的化学物质通过浸润冻结核化(Pruppacher and Klett，1997; Khain et al.,2005)；近几年一些学者探讨了生物气溶胶(如细菌、真菌和藻类)作为冰核影响冰晶形成，尤其典型的冰核活性细菌能有效地促进其他冰核的异质核化过程(Morris et al.,2004；Gonçalves et al.,2012)。可见，气溶胶是参与冰晶形成的主要成分，它的分布可以改变云内的动力以及微物理过程，尤其是对冰晶的空间分布和谱变化影响最为突出。

目前云模式中与气溶胶数浓度和类型有关的冰晶核化方案尚少。初期的云模式中，冰晶数浓度由野外探测或云室实验得到的经验公式计算，冰核活化有不同的表达式，Fletcher(1962)认为冰核浓度是温度的函数，Cooper(1980)认为冰核浓度是饱和度的函数，许焕斌等(2004)和Thompson et al.(2008)则同时考虑空气饱和度和温度。但完全根据经验公式得到的冰晶谱与实际观测结果存在较大差异。Cooper(1980)、Saunders and Al-Naimi(1985)收集飞机观测和云室结果，他们的研究结果中并没有显示出冰晶数浓度和温度之间很好的相关性，此外，Ryan(2000)使用模型计算的冰晶数浓度在－17°C时为8×10⁵ m^-3，而在－47°C时减小到2～3×10⁴m^-3，他的结果显示冰晶数浓度并不是随着温度的降低而增加；Gultepe et al.(2000)在文章中总结了四个现场观测，尽管是在不同的时期和不同的地方进行的，但冰晶数浓度的分布范围与所对应的温区范围是相似的，在温区－46～－2°C时对应的冰晶数浓度量级为0.001～1×10⁶m^-3，但并未发现冰晶数浓度和温度之间很好的相关性。鉴于此，近几年许多学者用基于物理过程的参数化方案取代原有的经验公式，改进对冰云的模拟，如Barahona and Nenes(2009)给出的冰晶参数化方 案，该方案考虑了气溶胶的微物理化学特征以及同质核化和异质核化之间的竞争，根据不同粒度分布及化学组成的气溶胶计算了冰晶数浓度、粒度分布和最大过饱和度；Liu and Penner(2002)专门考虑了硫酸盐气溶胶的同质核化和烟尘气溶胶的异质核化并发现硫酸盐气溶胶液滴在205 K以下，呈现出类似于玻璃的特性，不再发生同质核化，并提出冰晶核化参数化方案可以很大程度上影响核化的发生概率，进而影响冰晶数浓度以及其他物理过程；Gultepe et al.(2000)提出风速、冰面过饱和度、气溶胶数浓度、温度均应被考虑到计算冰晶数浓度的参数化方案中，这和Heymsfield and Sabin(1989)对冰相卷云的研究结论一致。此外，风速是一个计算冰晶数浓度的重要参数，Blyth and Latham(1993)使用飞机观测夏季积云，他们发现不断增加的风速很可能导致低的冰晶数浓度，而在风速小于3 m s^-1 的地方冰晶数浓度很大。而冰晶对于冷云降水、云的辐射效应、水分循环甚至平流层的水汽含量、雷暴云的起放电等均有重要影响；冰晶谱分布能够改变冰相粒子的垂直结构，从而改变非感应起电率的时空分布，导致雷暴云内放电属性、放电位置以及闪电频次产生差异(Rawlins，1981；周志敏和郭学良，2009；谢屹然等，2005；Krämer，2009；Gonçalves et al.,2012)；冰晶粒子尺度和冰晶含量等均可影响卷云辐射特征，冰晶含量增加，尺度变小可以导致卷云在整个短波波段的反照率增大(Rangno and Hobbs，1994；毛节泰和李成才，2005；李娟和毛节泰，2006)；冰相物理过程增强，抑制雨水的生成，从而影响云的微物理过程(Pruppacher and Klett，1997；Morris et al.,2004；Khain et al.,2005；陈丽等，2007；王亚玲等，2012；Gonçalves et al.,2012)。因此，如何提高云模式中对冰晶核化过程的描述能力变得日益重要。

目前，国内现有云模式中冰晶核化参数方案大都采用依靠云室实验得出的经验公式，模拟所得的冰晶分布也与实际观测结果存在较大差异，除此之外，也难以探讨气溶胶浓度、化学组分等对冰晶数浓度、尺度、分布等的影响。基于此，本文在已有的二维对流云起、放电模式(谭涌波等，2006，2007，2012；Tan et al.,2006)中耦合了一个与气溶胶有关的冰晶核化方案，取代通过云室观测得到的经验公式；通过对比模拟试验，探讨新方案对冰晶比含水量空间分布的影响，以及新方案加入后云滴、雨滴、霰、雹比含水量空间分布发生的调整；分析新方案对五种水成物粒子比含水量、电荷结构的影响，为继后探讨气溶胶与雷暴云起放电过程、电荷结构之间的相互关系提供了有利条件。

本文利用马明和谭涌波等学者(马明，2004；谭涌波等，2006，2007，2012；Tan et al.,2006)在中国气象科学研究院发展的二维对流云数值模式(胡志晋和何观芳，1987；王谦和胡志晋，1990；于达维等，2001)的基础上，先后加入二维雷暴云起、放电参数化方案，建立的二维雷暴云起、放电数值模式，该模式中微物理过程根据积云中水粒子增长、下落速度以及物理特征的不同，将它们分成云滴、雨滴、冰晶(包括雪团)、霰(包括冻雨滴)、雹五种成分，模式包括21个预报量，与微物理过程相关的预报量有水汽、云滴、雨滴、冰晶、霰、雹的比含水量(Q_v、Q_c、Q_r、Q_i、Q_g、Q_h)和五种水成物粒子的比浓度(N_r、N_i、N_g、N_h)，单位分别为g kg^-1和kg^-1，其中云滴比浓度N_c取常数(4×10⁸ kg^-1)。云内的微物理过程考虑了云中凝结蒸发、碰并、自动转化、核化繁生、融化冻结等，包括积雨云中27种主要微物理过程，有冰、雨、云、霰、雹的凝结(凝华)和蒸发；冰、雨、霰、雹对云滴的碰并，雨和冰的碰并，霰、雹碰并雨，霰、雹碰并冰晶；冰晶核化、繁生；云雨转化，冰霰转化，霰雹转化，雨冻结成霰，霰、雹、冰融化成雨；冰晶碰并，雨滴碰并，雹的湿增长极限过程等。其中微物理过程参数化的具体方程组以及各种水成物粒子谱和下落末速度公式详见文献(胡志晋等，1987)，本文不再赘述。

在已有的微物理参数化方案的基础上，本文选用Liu and Penner(2005)发展的冰晶核化参数化方案(简称LP方案)取代原来模式中所采用的游来光和石安英(1964)通过云室观测得到的经验公式(简称YS方案)，并通过模拟对比试验，探讨两种方案在云微物理过程中的差异。因此除参数化方案外，模式中还需考虑气溶胶成分。在LP方案中气溶胶通过同质核化或作为冰核异质核化直接形成冰晶，此外，气溶胶也可充当云凝结核形成云滴，但气溶胶—云滴这部分工作已完成，详见文献(王宁宁等，2013a，2013b)，本文不再赘述。

参考Liu and Penner(2005)的做法，气溶胶暂作为背景场，只是为核化过程提供初始分布，而不作为模式中新的预报量。将气溶胶简单分为硫酸盐气溶胶、沙尘气溶胶和粉尘气溶胶，不考虑气溶胶的谱分布，仅考虑了气溶胶的数浓度。考虑到布朗运动和重力沉降作用，气溶胶随高度的增加按指数关系递减(Yin et al.,2000；王宁宁等，2013b)。其指数递减率为

其中，A₀表示高度z=0时的气溶胶数浓度，A_n表示不同高度对应的气溶胶数浓度。对气溶胶各个网格点初始化时，模式假设各层气溶胶数浓度是均匀分布的，所以只需给出地面气溶胶的初始浓度就可算出整个空间的气溶胶分布。

自然云通常是固、液相混合云，冰晶是形成其他大冰相粒子的“胚胎”，在冰相云微物理过程中有着重要作用(黄美元和徐华英，1999)，因此，合理确定冰晶数浓度就非常重要。LP方案考虑了与核化有关的诸多因素，下面简单介绍一下耦合进云模式中的新核化方案。LP方案结合了Meyers et al.(1992，1994)、Liu and Penner(2002)给出的接触冻结核化、浸润核化、沉积核化计算理论，Kärcher and Lohmann(2002a，2002b)通过实验室实验得到的硫酸盐同质核化方案，以及Kärcher and Lohmann(2003)在同质核化框架下得到的异质核化方案等。该方案中发生同质核化的为硫酸盐气溶胶，发生异质核化的为粉尘气溶胶和沙尘气溶胶。在一定温度和上升速度情况下，若有冰核的帮助，异质核化能发生在过饱和度(相对于冰)较低情况下，若冰核数目较少(小于10 L^-1)，生成的冰晶少，对水汽的消耗相对也小，过饱和度还可以继续上升，当过饱和度达到同质核化所要求时(至少0.4以上)，同质核化便可发生。通常过冷气溶胶溶液滴浓度很高，一旦有同质核化发生，便生成大量冰晶，过饱和度立即下降。因此，同质核化持续时间很短。若冰核的数浓度足够多，异质核化生成的冰晶可阻止上升到同质核化所需的过饱和度，同质核化不再发生，可见同质核化和异质核化之间还存在竞争关系(Liu and Penner，2005；Barahona and Nenes，2009)。在LP方案中，在给定的气溶胶数浓度和气块儿上升速度情况下，首先依据温度判断同质核化是否会发生，这个温度阈值(定义为T_c)是气溶胶数浓度和上升风速的函数，当气块儿温度高于T_c时，同质核化不可能发生，冰晶来自异质核化。当温度小于T_c－5时，同质核化发生，冰晶主要来自同质核化，仅由同质核化计算。当温度在T_c－5到T_c之间时，为过度区间，由同质核化和异质核化共同计算。下面简单介绍参数方案主要步骤。

(1)求温度阈值：

其中，N_dust+soot表示总的沙尘气溶胶和粉尘气溶胶的数浓度(m^-3)，w为上升速度(m s^-1)，T_c表示温度阈值，单位采用°C。

(2)若－64.0≤T≤T_c－5.0，可能发生同质核化：

H_Rw为湿度阈值，若气块湿度超过阈值H_Rw，同质核化才发生，方案中，将同质核化分成两类情况分别处理。一是冰晶快速增长，二是冰晶缓慢增长。依据温度和上升速度的关系，判断属于哪种情况。

(一)当T＞6.07×ln(w)－55.0时，冰晶快速增长，冰晶快增长公式为

(二)当T≤6.07×ln(w)－55.0时，冰晶增长缓慢，冰晶慢增长公式为

其中，a₁、a₂、b₁、b₂、b₃、c₁、c₂为给定系数或是与风速、温度有关的量，N_a表示硫酸盐气溶胶溶的数浓度(m^-3)，N_ia表示硫酸盐气溶胶同质核化生成的冰晶数浓度。

(3)若T_c≤T≤－5.0发生异质核化

(一)沙尘气溶胶发生沉积核化，计算公式为

N_id=1×10⁻³×exp{−0.639+0.1296[100×(H_Ri−1)]},      (7)

其中，H_Ri为冰面相对湿度，N_id表示沙尘气溶胶异质核化生成的冰晶数浓度。

(二)粉尘气溶胶发生浸润核化，计算公式为

其中，a₁₁、a₁₂、a₂₁、a₂₂、b₁₁、b₁₂、b₂₁、b₂₂为给定系数，N_s表示粉尘气溶胶的数浓度(m^-3)，N_is表示粉尘气溶胶异质核化生成的冰晶数浓度。

(4)T_c－5.0＜T＜T_c同质核化和异质核化同时发生

(一)T＜－40.0°C并且w＞1.0 m s^-1，同质核化占主导抑制异质核化的发生，只考虑同质核化。

(二)其他情况下两者同时发生，冰晶数浓度的计算公式为

其中N_i表示同质核化和异质核化同质发生时生成的冰晶数浓度。

(5)冰晶比浓度和比含水量的计算公式为

其中，N_ice为气溶胶核化生成的冰晶数浓度，ρ表示空气密度，m_i0=10^-10 g为单个冰晶的质量，H_i和Q_i分别表示最终计算得到的冰晶比浓度和比含水量。

图 1为2005年6月10日20:00(北京时，下同)长春站环境温、湿层结，0°C层位于685 hPa左右，地面温度为21°C左右，露点温度为14°C左右。低层大气相对湿度为40%～75%，在400～600 hPa间层结呈现弱的热力不稳定状态。周志敏等(2012)使用三维强风暴冰雹分档模式对此探空个例进行了模拟和分析证明为强降雹个例，且该次对流过程发展迅速且强烈，过程伴有雷鸣和闪电，表明云内电场能发展到击穿阂值，因此该探空曲线适合用来模拟一次雷暴过程。在初始时间给出了一个水平半径为5 km，垂直半径为1 km的湿热泡扰动，中心位于模拟区域中心左侧，高度为1 km的格点上，垂直扰动半径为40 km，水平扰动半径为6 km，温度扰动和湿度扰动分别取0.5 K和72%，从中心向外以余弦函数递减。模式的计算域为76 km×20 km，模拟时间为80 min，时间步长2 s。地面气溶胶的初始浓度参考2006年春季飞机在天津—唐山一带进行航测得到的观测资料(杨素英等，2010)，沙尘气溶胶和硫酸盐气溶胶的地面初始浓度取100 m^-3，粉尘气溶胶地面初始浓度取10 m^-3。

图 1

Fig. 1

图 1 环境温、湿层结：(a)露点温度和环境温度；(b)相对湿度 Fig. 1 Environment temperature and humidity stratification:(a)The dew point temperature and environmental temperature;(b)the relative humidity

在耦合新冰晶核化方案的二维起、放电云模式的基础上，对YS经验公式和LP方案进行对比验证。核化方案直接影响冰晶空间分布，其他微物理过程受冰晶核化的影响发生不同程度的变化，进而造成水成物粒子空间分布发生调整，因此有必要将模拟得到的冰晶空间分布与国内外一些模拟和观测结果进行对比验证。若LP方案得到的冰晶空间分布正确，本文认为由核化引起的其他云内水成物粒子(云滴、雨滴、霰、雹)空间分布变化也是合理的。

图 2给出了LP方案和YS方案Q_i在31、35、45、55、65、75 min的空间分布。从图 2可以明显看出：(1)在冰晶发生核化的初始阶段(雷暴云发展至31～40 min)，两种不同核化方案背景下，冰晶出现的高度差异明显，其中LP方案下的冰晶中心高度在6 km左右，而YS方案冰晶则出现在8 km左右，相比YS方案，LP方案中浓度大于0.1 g kg^-1的冰晶垂直空间分布区域更大；(2)40～65 min为核化发展最为旺盛阶段，LP方案和YS方案下模拟得到的冰晶空间分布情况存在明显区别，LP方案 中冰晶空间分布区域大于YS方案，除45 min在 10 km左右处的冰晶水平分布增大外，此阶段浓度大于0.5 g kg^-1的冰晶垂直高度区间同样增加明显，例如：在45 min时刻，耦合LP方案后浓度大于0.5 g kg^-1冰晶的垂直分布区间由YS方案中的6.0～8.8 km范围变为3.3～9.0 km，此外，在旺盛阶段LP方案中的冰晶最大浓度值大于YS方案；(3)在雷暴云消散阶段(65～75 min)，浓度大于0.5 g kg^-1的冰晶存在高度区间和最大浓度值依然大于YS方案。

图 2

Fig. 2

图 2 冰晶比含水量的空间分布：(a)冰晶核化参数化方案(LP方案)；(b)云室观测经验公式(YS方案) Fig. 2 The specific water content distribution of the ice crystal:(a)ice crystal nucleation empirical formula(LP scheme);(b)cloud chamber observation experience formula(YS scheme)

为了更加细致的了解两种不同冰晶核化方案对冰晶空间分布的影响，表 1给出了Q_i在不同时刻对应的空间分布高度和温度区间。从表 1可以看 出：(1)雷暴云发展至35 min时(冰晶核化发生初始阶段)，LP方案中浓度大于0.1 g kg^-1的冰晶分布在4.2～7.5 km高度区间以及－14.6～－40.7°C温度区间，YS方案则分布在7.2～8.4 km，－39.2～ －46.9°C，因此，LP方案和YS方案高度区间段和温区段分别为3.3 km，26.1°C和1.2 km，7.7°C，LP方案高度区间段是YS方案的2.8倍，温度区间段是YS方案的3.4倍；(2)40～65 min 为核化发展旺盛阶段，来自LP方案中的浓度大于0.5 g kg^-1的冰晶垂直空间分布为：3.3～9.0 km，－7.6～ －50.1°C，而YS方案中其分布在5.5～9.0 km，－24.2～－50.1°C，即LP方案和YS方案高度区间段和温区段分别为5.7 km，42.5°C和3.5 km，25.9°C，此时LP方案高度区间段和温度区间段是YS方案的1.6倍；(3)在雷暴云的消散阶段(65～80 min)，不同方案中浓度大于0.5 g kg^-1的冰晶空间分布同样存在差异，LP方案(高度区间：3.5～8.5 km，温度区间：－9.1～－49.1°C)和YS方案(高度区间：5.9 km～8.7 km，温度区间：－27.8～－48.5°C)高度区间段和温度区间段分别为5.0 km，40.0°C和2.8 km，20.7°C，LP方案高度区间段是YS方案的1.8倍，温度区间段是YS方案的1.9倍。从整体来看，LP方案得到的冰晶空间分布高度和温度区间分别在3.3～9.0 km，－50.1～－7.6°C范围内，YS方案背景下核化发生在5.5～9.0 km，－50.1～－24.2°C。

表 1(Table 1)

表 1 LP和YS方案值较大的Qi在不同时刻的分布高度和温区 Table 1 The height and temperature region of the larger value Qi of LP and YS schemes at different time 模拟时 间/min Qi对应的高度/km Qi对应的温区/°C LP方案 YS方案 LP方案 YS方案 35 4.2～7.5 7.2～8.4 －14.6～－40.7 －39.2～－46.9 40 5.0～8.5 7.5～9.0 －20.3～－49.1 －40.7～－50.1 45 3.3～9.0 6.0～8.8 －7.6～－50.1 －28.3～－49.1 50 3.5～8.6 5.7～8.5 －9.1～－48.2 －26.0～－49.1 55 3.7～8.5 5.5～8.0 －10.9～－49.1 －24.2～－44.0 60 4.0～8.5 5.6～8.7 －13.2～－49.1 －25.0～－48.5 65 3.6～8.2 5.7～7.6 －9.4～－45.2 －26.0～－43.9 70 3.5～8.5 5.9～7.7 －9.1～－49.1 －27.8～－41.8 75 3.6～7.7 6.0～7.8 －9.4～－41.8 －28.3～－42.8 80 3.5～7.9 6.0～8.7 －9.1～－43.2 －28.3～－48.5 注：表中的高度和温区仅考虑了值较大的Qi(模拟时间35 min为Qi大于0.1 g kg-1的分布，其余时刻均为Qi大于0.5 g kg-1的分布)

表 1 LP和YS方案值较大的Qi在不同时刻的分布高度和温区 Table 1 The height and temperature region of the larger value Qi of LP and YS schemes at different time

由此可见，在雷暴云整个发展过程中，与YS方案相比LP方案得到的冰晶高度以及温度分布区间要更大。YS方案产生的冰晶主要集中在低温区(小于－25°C)，在大于－25°C的温区几乎没有冰晶存在，而LP方案在高温区(大于－25°C)和低温区(小于－25°C)均有大量冰晶存在，且LP方案在整个雷暴云的发展过程中冰晶的分布高度、温度区间以及最大浓度值均大于YS方案。

模式中与Q_i有关的源项为核化过程、冰晶繁生以及新生冰晶和云滴的碰并，而核化产生的冰晶量远远大于其他两项，对Q_i产生影响的主要为核化过程(孔凡铀等，1991a，1991b)。模式中冰晶比浓度变化率和单个冰晶质量(模式中冰晶质量取定值)的乘积为比含水量变化率，所以核化产生的冰晶数浓度直接影响Q_i的空间分布。因此造成以上现象的主要原因是不同冰晶核化方案得到的冰晶谱不同进而造成Q_i空间分布差异很大(Fletcher，1962；Huffman，1973；Meyers et al.,1992；肖辉等，2012)。YS方案中基于云室观测得到的冰晶谱分布很难观测到每个温度下的冰晶数浓度，只是按照温度高低选取几个温度，记录冰晶的日变化，根据这几个高低不同温度所对应的冰晶数浓度的变化趋势建立与温度相关的函数关系，但由于气溶胶种类不同以及随高度变化的空间分布，温度、湿度、压强等随高度的变化，以及不同地区不同条件下冰核谱分布的差异等，所以传统的以温度为函数的YS经验公式存在不准确性(游来光和石安英，1964；杨磊等，2013)。原方案采用YS经验公式随温度的降低呈e指数递增，只是针对具体条件对某些参数进行了修订和限制如在云水含量较低或过饱和度为0的空间，核化不发生，在温度大于－40°C时，冰晶数浓度随着温度的降低而增加，低于－40°C时按照 －40°C计算，所以在高温区(大于－25°C)，核化产生的冰晶量极少，在－40°C左右为冰晶主要集中的温区，而LP方案中用气溶胶计算冰晶数浓度，生成的冰晶数浓度与温度、过饱和度和上升风速 等多个量相关，所以冰晶数浓度并不同经验公式随着温度的降低而增加，它与温度并无直接的相关性，而是与风速、相对与冰的过饱和度，和气溶胶数浓度、温度等多个物理量有关，考虑的更为全面(Heymsfield and Sabin，1989；Meyers et al.,1992；Gultepe et al.,2000)。

关于冰晶分布已有大量的模拟和探测资料，许多学者(Cooper，1980；Saunders and Al-Naimi，1985；Meyers et al.,1992；Heymsfield and Milosohevich，1993；Ryan，2000)收集飞机观测和云室结果，研究显示冰晶数浓度主要集中在10⁶ m^-3，它并不是随着温度的降低而呈e指数增加，并未发现冰晶数浓度和温度之间很好的相关性。由于YS方案冰晶数浓度随温度的降低呈e指数增加造成核化主要发生在低温区(小于－25°C)，显然YS方案高估了低温区冰晶的核化能力。此外，YS方案在高温区(大于－25°C)几乎无冰晶的存在，但研究发现气溶胶粒子在云水过饱和度更高的条件下，在大于－12°C时可自发发生核化产生高达10⁶ m^-3的冰晶浓度，沙尘气溶胶在低于－8°C时通过接触冻结核化、低于－15°C发生凝华核化(Rangno and Hobbs，1985；Levin et al.,2005；李磊，2007)，可见YS方案中基于云室观测得到的经验公式低估了高温区(大于－25°C)冰晶核化能力，此外许多探空观测也证实了这一点，如CCOPE实验、自动光学冰粒子计数器和二维云探针测量和观测以及一些资料统计结果证明冰晶在云中－15～－4°C，平均浓度在4×10⁵ m^-3左右，最高浓度可达10⁶ m^-3左右(Turner and Radke，1973；Knollenberg，1981)。

综上所述，LP方案中冰晶数浓度并不同YS方案随温度的降低呈e指数增加，而是与上升风速、相对与冰的过饱和度、气溶胶数浓度、气溶胶种类和温度等多个物理量有关，造成新核化方案和YS经验公式得到的Q_i空间分布存在明显差异，主要体现在LP方案和YS方案冰晶核化发生的主要高度和温度区间分别集中在3.3～9.0 km，－50～ －7.6°C左右和5.5～9.0 km，－50.1～－24.2°C左右，YS方案冰晶主要集中在低温区(小于－25°C)，在高温区(大于－25°C)几乎无冰晶的存在，而LP方案冰晶的空间分布高度和温度区间大于YS方案，在高温区和低温区均分布大量的冰晶，通过对比国内外的相关探测和模拟结果来看新的核化方案要比原经验公式更为合理。

雷暴云中起、放电特征，依赖于云中微物理过程的发展和各粒子的空间分布特征，其中在非感应和感应起电机制下粒子所得到的电荷量大小、主负电荷区高度、以及电荷载体等均与各水成物粒子的空间分布相关，由此可见微物理过程的发展和各水成物粒子的空间分布特征和雷暴云的起电存在密切而又复杂的关系(言穆弘等，1996；孙平安等，2000；谢屹然等，2005)，本文在二维对流云起、放电模式中耦合一个新冰晶核化方案，取代通过云室观测得到的经验公式后，致使核化过程发生变化，进而影响其他微物理过程的变化，水成物粒子空间分布和数值大小势必会改变？通过详细讨论两方案引起的不同微物理过程变化，所导致的各水成物粒子的空间分布的不同，能为继后分析雷暴云起电提供理论依据。

图 3给出了各水成物最大比含水量(Q_max)随模拟时间的变化；其中图 3a为以YS经验公式计算的冰晶数浓度各水成物最大比含水量随模拟时间的变化；图 3b为耦合了LP参数化方案后各水成物最大比含水量随模拟时间的变化。表 2给出了各水成物最大比含水量的峰值以及到达峰值的时间。

图 3

Fig. 3

图 3 LP方案和YS方案各水成物最大比含水量(Qmax)随模拟时间的变化：(a)LP方案；(b)YS方案。Qc、Qr、Qi、Qg、Qh分别代表云滴、雨滴、冰晶、霰、雹的比含水量 Fig. 3 The peak emergence time and value of maximum specific water content(Qmax)of hydrometeors of LP and YS schemes:(a)LP scheme;(b)YS scheme

表 2(Table 2)

表 2 LP方案和YS方案各水成物粒子最大比含水量的峰值以及到达的时间 Table 2 Peak and arrival time of the max mass of hydrometeors of LP and YS schemes 水成物 粒子 最大比含水量的峰值/g kg-1 到达峰值的时间/min LP方案 YS方案 LP方案 YS方案 云 6.19 6.26 31 30 雨 0.96 2.47 33 32 冰 3.48 2.18 48 50 霰 5.77 7.82 43 43 雹 0.0019 0.0637 44 36

表 2 LP方案和YS方案各水成物粒子最大比含水量的峰值以及到达的时间 Table 2 Peak and arrival time of the max mass of hydrometeors of LP and YS schemes

新冰晶核化方案对冰晶和霰的影响最大，进而对雨滴等含量也有很大影响。对云滴和雨滴的影响主要从17 min和30 min后开始，加入气溶胶后使云水消耗变快，雨滴在34 min后出现急剧下降56 min后又缓慢回升的过程，在42 min时雨滴出现0.0000012 g kg^-1的极小值。这主要是由于高温时(大于－25°C)LP方案中的冰晶数浓度高于YS方案中的，加入气溶胶后，在温度较高的过冷区中产生大量的冰晶，争食云中的水分，由于贝吉龙效应(Bergeron)，冰晶生长比过冷云滴和小雨滴快得多，使得大量过冷水提前冻结而转化到冰相粒子上，同时，减慢过冷云滴的生长及其向雨水的自动转化速率，所以云中高温冰晶数浓度高导致云中液态水提早消耗，阻止云滴的凝结增长和雨滴的长大(胡志晋等，1987)。此外，冰晶数浓度最大量级集中在35～60 min左右，因此在35～60 min时，冰晶数浓度是整个雷暴云中值最大的时段，此时有大量冰晶和雨滴发生碰撞，消耗掉雨滴，这是雨滴比含水量出现急剧下降的主要原因，随后的缓慢回升是由于高空中的霰和冰晶在0°C层下融化造成的(Pruppacher and Klett，1997；肖辉等，2004；刘校辰，2006；杨慧玲等，2011)。

与霰比含水量有关的微物理过程繁多，有冰晶碰并霰，冰晶碰并雨滴，雨滴碰并冰晶，霰碰并云滴，雨滴冻结成霰，冰晶转化为霰、云滴自动转化成霰等，对比YS方案，其中造成Q_g值变化较大的物理过程有：雨滴冻结成霰和冰晶和雨滴间的碰并。LP方案得到的Q_g值基本均小于YS方案，其主要原因是Q_r在这时段大幅度减小，致使雨滴冻结成霰和雨滴碰并冰晶两个物理过程要比YS方案减弱很多。

从图 3可看出：LP方案的Q_i最大值基本均大于YS方案中的Q_i。冰晶核化参数化方案可以很大程度上影响核化的发生概率，进而影响冰晶数浓度以及其他微物理过程，不同的核化方案得到冰晶谱不同，必然会造成Q_i在出现时间、峰值大小等方面差异，如Meyer et al.(1992)的参数化公式，在较暖的温区可给出较高的活化冰晶数浓度；而Fletcher(1962)和Huffman(1973)的参数化公式，在较低的温区则产生出较高的活化冰晶数浓度。

相比YS方案雹的比含水量有所减小，与雹的形成和增长有关的过程有：雹碰并霰，霰雹自动转化，雹的凝华增长，雨滴碰并雹，冰晶碰并雹等，大多数物理过程都有不同程度的减弱，其中减弱最为明显的是霰向雹的自动转化和雹的增长过程，前者减弱的主要原因是，气溶胶加入后使大霰粒子所占比例减小，从而造成霰向雹的转化率减小；后者减小的主要原因是，雹的增长过程与雹比浓度和雹的直径有关，气溶胶的加入使雹比浓度和直径均减小，所以造成Q_h的减小。

如图 4所示由LP方案得到的云滴和雨滴空间分布不同于YS方案，造成这种现象的原因同上述造成云滴和雨滴比含水量值不同的原因，新核化方案在温度较高的过冷区中产生大量的冰晶，争食云中的水分，由于贝吉龙效应，冰晶生长比过冷云滴和小雨滴快得多，使得大量的过冷水提前冻结而转化到各种冰相粒子上，同时，减慢过冷云滴的生长及其向雨水的自动转化速率，所以云中高温冰晶数浓度高导致云中液态水提早消耗，阻止云滴的凝结增长和雨滴的长大。此外，生成的冰晶和雨滴发生碰撞，消耗掉雨滴，这也是LP方案得到的雨滴在核化阶段小于由YS方案得到的雨滴的原因。随后雨滴的缓慢回升是由于高空中的霰和冰晶在0°C层下融化造成的。

图 4

Fig. 4

图 4 云滴、雨滴、霰比含水量的空间分布：(a、b、c)LP方案；(d、e、f)YS方案 Fig. 4 The specific water content distribution of the cloud drop，rain and graupel:(a，b，c)LP scheme;(d，e，f)YS scheme

除上述提到的Q_g值小于由YS方案计算得到的值外，从图 4可以看到另一个最大不同是Q_g在4 km以下的空间分布，LP方案得到的Q_g在40 min时全部分布在4 km以上，在45 min时Q_g全部分布在2.3 km以上，在50 min时Q_g全部分布在2.2 km以上，在55 min时Q_g全部分布在2.1 km以上，而YS方案得到的Q_g在40 min时全部分布在1.5 km以上，在45 min时Q_g全部分布在0.5 km以上，在50 min时Q_g全部分布在0.7 km以上，在55 min时Q_g全部分布在1.0 km以上。4 km以下Q_g的空间 分布主要受雨滴影响，从图 3可知由LP方案得到的雨滴在35～60 min为整个雷暴云峰谷阶段，此阶段LP方案的Q_r要小于YS方案4个数量级，因此大幅度的减弱了雨滴冻结成霰，雨滴碰并冰晶生成霰，雨滴碰并霰生成霰三个物理过程，这是导致Q_g在 4 km 以下的空间分布不同于YS方案的主要原因。

由于Q_h的值相比其他四种水成物粒子值极小且有的时段为0，所以图略。LP方案得到的雹比含水量空间分布范围和值均要小于YS方案，原因同上述，与雹有关的多个物理过程均有不同程度的减弱。

在起电过程中，冰晶—霰(软雹)的非感应碰撞过程是雷暴云中主要荷电机制，因此冰晶的物理特征、出现时间等影响电过程的强弱、主负电荷区的高度及分布区域、和电荷载体等(孙安平等2002，2004)；此外，冰晶—霰(软雹)的非感应碰撞过程是雷暴云中强电场产生的根本原因，因为感应起电需要在一定的雷暴环境电场作用下发生，因此冰晶粒子又间接影响感应起电过程，而气溶胶是参与冰晶形成的主要成分，它的分布特征不仅改变高层冰晶的分布，而且还通过核化过程影响云的微物理过程，而雷暴云的起电机制，依赖云中微物理过程的发展和各粒子的空间分布特征，尤其是冰相粒子的分布特征，由3.1节知，考虑气溶胶后的核化方案得到的冰晶空间分布和原经验公式得到的冰晶空间分布特征存在明显差别，霰、雨滴、云滴、雹四种粒子的空间分布也略有不同，雷暴云的起电机制与云微物理过程以及水成物粒子息息相关，因此考虑气溶胶后的核化方案必然对雷暴云内起电过程产生影响。

图 5给出了LP方案和YS方案雷暴云在39、45、55、70 min四个时刻的空间电荷结构分布，从图 5可以明显看出：(1)在雷暴云初期时刻，LP方案和YS方案模拟得到的电荷结构空间分布情况存在明显区别，LP方案中为四极性电荷结构，而YS方案中则为偶极性电荷结构；(2)在雷暴云发展旺盛与消散阶段(40～80 min)，两种不同核化方案背景下，电荷的空间分布差异明显，LP方案中电荷空间分布区域和电荷量均大于YS方案；(3)在整个雷暴云的发展过程中两方案模拟得到的主正电荷区和主负电荷区的中心高度不同。

图 5

Fig. 5

图 5 电荷结构的空间分布：(a)LP方案；(b)YS方案 Fig. 5 The charge structure distribution:(a)LP scheme;(b)YS scheme

由此可见，在雷暴云整个发展过程中，LP方案与YS方案在不同时刻得到的空间电荷在雷暴云初期结构不同，在发展旺盛与消散阶段LP方案中电荷空间分布区域和电荷量均大于YS方案，此外，在不同时刻主正电荷区和主负电荷区的中心高度不同。

利用耦合了新冰晶核化方案的二维对流云起、放电模式对长春地区降雹型强雷暴个例的微物理特征和起电过程进行了数值模拟，验证了新的核化方案较原计算冰晶数浓度的经验公式更为合理，分析了LP方案的加入造成雷暴云一些微物理过程与电荷结构不同于YS方案的原因，得到以下结论：

(1)冰晶数浓度在LP方案中与上升风速、相对与冰的过饱和度、气溶胶数浓度、气溶胶种类和温度等多个物理量有关，而在YS方案随温度的降低呈e指数增加，造成新核化方案和原经验公式得到的Q_i空间分布差异明显，主要体现在LP方案和YS方案冰晶核化发生的主要高度和温度区间分别集中在3.3～9.0 km，－50.1～－7.6°C左右和5.5～9.0 km，－50.1～－24.2°C左右，YS方案冰晶主要集中在低温区(小于－25°C)，在高温区(大于 －25°C)几乎无冰晶的存在，而LP方案在高温区和低温区均分布大量的冰晶，在整个雷暴云的发展过程中LP方案冰晶的分布高度、温度区间以及最大浓度值均大于YS方案。通过对比国内外的相关探测和模拟结果来看新的核化方案要比原经验公式更加合理。

(2)冰晶核化方案的改变使五种水成物粒子的空间分布发生了调整，相比YS方案的模拟结果，LP方案中高温冰晶数浓度高导致云中液态水提早消耗，阻止云滴的凝结增长和雨滴的长大；雨滴冻结成霰、雨滴碰并冰晶生成霰、雨滴碰并霰生成霰等多个物理过程发生变化，造成LP方案得到的Q_g值小于YS方案，在4 km以下的空间分布不同于YS方案的模拟结果；与雹的形成和增长有关的过程基本均有不同程度减弱，减弱最为明显的是霰向雹的自动转化和雹的增长两个微物理过程。

(3)雷暴云的起电机制与云微物理过程以及水成物粒子息息相关，尤其是冰相粒子的空间分布，LP方案与YS方案对云的微物理过程影响不同，得到的冰晶空间分布差异明显，造成在雷暴云整个发展过程中，LP方案与YS方案在不同时刻得到的空间电荷在雷暴云初期结构不同，在发展旺盛与消散阶段LP方案中电荷空间分布区域和电荷量均大于YS方案，此外，主正电荷区和主负电荷区的中心高度不同。

本文利用与气溶胶有关的冰晶核化方案取代原云室观测得到的经验公式后，通过对比模拟试验证明核化方案对云的微物理过程影响明显，特别是冰晶的空间分布，而雷暴云中的起电过程与冰相粒子的物理特征息息相关，因此本文利用耦合新核化方案的云模型初步探讨了气溶胶—冰晶对雷暴云内起电过程的影响，而气溶胶谱分布、化学成分等如何影响雷暴云起、放电过程将是下一步的工作重点。此外，气溶胶不仅能够充当云凝结核活化成为云滴，而且可以成为冰核核化生成冰晶粒子，而云滴和冰晶是雷暴云中主要荷电贡献粒子，因此气溶胶如何影响雷暴云起电过程、电荷结构等成为今后的主要工作方向。