2. 河北省人工影响天气办公室, 石家庄 050021
2. Weather Modi fication Of fice of Hebei Province, Shijiazhuang 050021
云在全球能量与水循环中具有非常重要的作用。云的反照率与云滴数浓度的关系十分密切,而云滴数浓度又决定于云凝结核 (CCN) 浓度和云内水汽状况[1]。
硫酸盐粒子是海洋上空CCN的主要来源之一。海洋中二甲基硫 (DMS) 的排放可以增加硫酸盐气溶胶的数浓度,进而改变云滴数浓度,从而增大海洋上云的反射率,而海洋上层云覆盖了海洋表面的25 %,它对于全球的辐射平衡有很大的影响,这就是所谓的DMS-云-气候假设[2~4]。最近的观测和研究发现,随着人类活动的影响,硫酸盐气溶胶有增加的趋势[5~6],全球已经形成北美、欧洲和东亚3个硫酸盐的高值中心[7]。Twomey曾做过云中液态水含量不随CCN数目而变化的假定[5]。人们普遍认为,硫酸盐气溶胶的增加,将导致云滴数目的增加和云滴半径的减小[1, 5]。
硫酸盐气溶胶对于云和气候的影响被人们广泛认识和研究,但是准确估计硫酸盐气溶胶对云特性的影响需要全面了解包括气溶胶和云的微物理过程。近期有研究表明海盐粒子在此过程中有不可忽略的作用[8]。本文研究的目的就是利用详细的气溶胶和云微物理分档模式,来定量评估海盐气溶胶和硫酸盐气溶胶对海洋上云微物理特征的影响。
国内在云的分档模式方面也做了一些研究。肖辉等[9~10]利用云分档模式研究了盐核谱和浓度、碰并过程、大气层结等对云滴谱形成的作用。许焕斌等[11]的研究指出一些重要的云-降水物理过程应当用分档方法来处理,这在学科发展上和实际应用中都是有必要的。赵春生等[12~13]建立的气溶胶分档模式最初是研究海洋边界层中大气气溶胶的发生发展,并应用到研究海洋DMS对于CCN的影响。本文利用该模式研究海洋上硫酸盐气溶胶和海盐气溶胶粒子对海洋云微物理特性的影响。
1 模式介绍本模式中的分档气溶胶模式最初由赵春生等[12~13]建立并用来研究遥远海洋大气边界层中气溶胶产生和发展的演变过程,后来改进的微物理模式可用来研究CCN的生长过程和云物理过程。模式考虑的微物理过程有凝结、蒸发过程,粒子的核化过程和粒子间的相互凝并过程。模式中粒子共分4种类型,分别为硫酸盐粒子、海盐粒子、混合粒子 (包括硫酸盐和海盐) 及云滴,各分为100档。其中前3种气溶胶粒子半径范围为0.01~4 μm,云滴粒子半径范围为0.02~1330 μm。模式中各微物理过程的处理如下:
1.1 CCN核化过程每个时步模式都由当时的水汽浓度和温度显式计算过饱和比,为了保证计算的稳定性,我们采用可变步长的方法。当过饱和比超过临界饱和比时,凝结核 (CN) 就被核化成相应大小的含水粒子,进行凝结生长过程。平衡过饱和比随半径变化的关系由Kohler方程计算得到[1] :
|
(1) |
式 (1) 中,er是液滴表面的饱和蒸汽压,esat,w是水面的饱和蒸汽压,M是水分子量,Ms是盐分的分子量,σw是液滴的表面张力,ρw是水密度,r是云滴半径,γ是溶解于水中盐的离子数,Φs是渗透系数,ms是盐粒子质量,ρs是液态盐溶液的密度。
平衡过饱和比最大值,即临界过饱和比为:
|
(2) |
其中,
粒子的凝结生长过程的处理方法采用赵春生等的计算方法[12~13]。云滴的生长速率决定于准稳态的质量和热量的扩散输送过程。经CCN活化的粒子通过凝结过程的生长速率采用Pruppacher和Klett的计算方法[1] :
|
(3) |
式 (3) 中,SS是每个时步由温度和水汽含量计算的过饱和比,ka是热扩散系数,L是水的凝结潜热,Dv水汽扩散系数。
1.3 碰并过程本模式考虑多组分粒子,碰并过程主要考虑两种:布朗运动和重力碰并 (由沉降速度不同引起)。布朗运动对小粒子作用明显而重力收集过程是大于几个微米的粒子间碰并的主要机制。不同碰并机制的过程在粒子浓度方程的碰并项中表现为不同的碰并核。本模式中碰并系数的确定和对于碰并过程的计算方法采用赵春生等的计算方案[13]。体积为v的粒子由于碰并带来的浓度在Δt时间内变化方程为:
|
(4) |
其中,C是粒子的数浓度,β是碰撞粒子的碰并核,f是碰并参数。
1.4 干沉降过程虽然重力沉降过程仅对于大粒子比较明显,但它对于一个完善的气溶胶云模式是非常重要的。本模式中不同大小粒子有着不同的下落速度,下落速度的计算公式采用Pruppacher等人的计算方法[1]。
2 海盐和硫酸盐气溶胶对云滴数目的影响 2.1 模拟设计O' Dowd等的观测结果[14]指出,硫酸盐的数浓度分布模态可以用平均半径0.08 μm,标准偏差1.4的对数正态分布来表示。后来O' Dowd等对海盐的观测[15]表明,海盐可以用3个对数正态分布来表示,分别为film模态、jet模态和spume模态,并且随着水平风速的变化,海盐粒子的浓度也会不同 (如表 1所示)。
|
|
表 1 海盐粒子浓度随水平风速的变化 |
根据以上观测结果,应用模式计算过程中,我们对硫酸盐浓度分别取100,300和1000 cm-3。由于海盐3个模态中第三个模态数量级相对前两个很小,因此海盐只考虑前两个模态。模式运行初始时刻,设置温度为280 K,大气压强为1000 hPa,相对湿度为100 %。
2.2 结果与分析图 1给出的是根据寇拉曲线计算的硫酸盐粒子活化比随最大过饱和比的变化曲线。由图 1可以看出:当最大过饱和比比较低 (0 %~0.3 %) 的时候,硫酸盐粒子随最大过饱和比的变化很显著。而当最大过饱和比比较高的时候,由于硫酸盐粒子绝大部分都已经被活化,所以基本不随最大过饱和比的变化而变化。
|
|
| 图 1. 由寇拉曲线计算出的不同过饱和比情况下硫酸盐粒子 (平均半径0.08 μm, 标准偏差1.4的对数正态分布) 的活化比 | |
图 2给出的是硫酸盐粒子浓度分别在100,500和1000 cm-3情况下,随着上升速度和海盐粒子浓度的不同,最大过饱和比的变化。由图 2可以看出:增加硫酸盐和海盐浓度都会导致最大过饱和比的降低,这是由于增加的粒子吸湿增长所造成的;增大上升速度将导致最大过饱和比的升高。
|
|
| 图 2. 硫酸盐粒子浓度分别为(a)100cm-3, (b)300cm-3, (c) 1000cm-3情况,最大过饱和比随上升速度和海盐粒子浓度 (以水平风速表示,见表 1)的变化 | |
图 3是在不同的海盐和硫酸盐情况下,活化粒子总数随上升速度的变化。从图中可以看出:当硫酸盐粒子浓度相对较低而且上升气流强,绝大部分硫酸盐粒子被活化,我们称为“水汽相对充足”的情况,此时增加海盐粒子的浓度将导致活化的粒子总数增加;而如果硫酸盐粒子浓度相对较高而且上升气流弱,相当数量的硫酸盐粒子不能被活化,我们称为“水汽相对不足”的情况,此时增加海盐粒子的浓度将导致活化的粒子总数减少。从物理意义上来讲,当硫酸盐粒子浓度低而上升气流强时,海盐粒子的活化数目显著增大,从而导致总的粒子活化数目的增加;而当硫酸盐粒子浓度高而上升气流弱时,海盐粒子相对很少,但是由于海盐粒子相对比硫酸盐粒子大,会先活化吸收相对不足的水汽从而导致最高过饱和比较低,也就使得硫酸盐的活化数目大大降低,从而使总的粒子活化数目降低。
|
|
| 图 3. 活化粒子总数与垂直风速和硫酸盐粒子浓度的关系 | |
3 海盐和硫酸盐气溶胶对云微物理特性的影响 3.1 模拟设计
将O'Dowd等的观测结果[14~15]作为初始条件来研究海盐和硫酸盐在云形成过程中的作用。在模式计算过程中,硫酸盐浓度取100,300,1000,5000 cm-3分别代表清洁大气、典型大气、污染大气和严重污染大气的情况。海洋上空最经常出现的层云和层积云,根据观测[16]高度约为1.2 km。因此,模式运行初始时刻,设置各物理参数 (云底) 为:温度为280 K,大气压强为880 hPa,相对湿度为100 %。据观测数据[16],云厚度约为0.5km左右,我们根据这个高度,再由云团上升速度的不同来决定模式模拟的时间。
3.2 模拟结果及分析(1) 首先,我们取云团上升速度为0.5 m·s-1运行模式,当云团上升到云顶即500 m时所得云内微物理特征和海盐、硫酸盐浓度的关系如表 2所示。其中最大过饱和比为云团上升过程中所达到的最大过饱和比,有效半径是云滴谱的三阶矩与二阶矩的比值,平均半径和标准偏差是将云滴在云顶时的尺度分布做一个对数正态分布的拟合所得到的正态分布的平均半径和标准偏差。
|
|
表 2 云团上升速度0.5 m·s-1时, 云顶各微物理特性随模拟初始条件的变化 |
由表 2可知:①云中液态水含量随初始海盐、硫酸盐粒子数目的变化非常小,最大值与最小值相差不到2 %,基本上可以视为与初始干粒子的浓度无关。这说明云中液态水含量主要取决于水汽条件。如果能确定云中的云滴数目和大气气象条件,即可以粗略地估计云滴的半径。这一特点对于研究云顶的辐射特性以及卫星应用有着重要的意义。②在0.5 m·s-1的上升速度情况下,水汽相对充足,因此云滴的数目主要取决于硫酸盐粒子和海盐粒子的数目之和。由于硫酸盐粒子数目明显多于海盐粒子数目,云滴数目跟硫酸盐数目关系最为明显,云滴数目跟初始硫酸盐粒子数目为同一量级,而且相差很小。海盐粒子的作用为简单地增加了云滴数目。③硫酸盐粒子的增加,即大气污染程度的增加,明显地增加了云滴数目,由于液态水含量不变,因此明显减小了云滴的半径,而海盐粒子对云滴半径的影响不大。
(2) 取云团上升速度为0.2 m·s-1运行模式,当云团上升到云顶即500 m时所得云内微物理特征和海盐、硫酸盐浓度的关系如表 3所示。
|
|
表 3 云团上升速度0.2 m·s-1, 云顶各微物理特性随模拟初始条件的变化 |
由表 3可知:①对于这种条件下,硫酸盐对云滴的影响与表 2相同,即增加硫酸盐粒子的结果是增加了云滴的数目,减小了云滴的半径。但是海盐粒子的影响比较复杂。对于清洁和典型大气,海盐的影响与表 2相同,海盐粒子的增加直接增加了云滴的数目。但是对于后两种情况,由于水汽相对不足,海盐粒子增加导致的过饱和比的减少使得硫酸盐粒子活化为云滴的数目减少,从而减少了云滴的数目。②与表 2相比,在减小了云团上升速度后,对于水汽相对充足的清洁大气和典型大气,云滴的数目和半径变化都不大。而对于后两种情况,云滴的数目和半径均有比较明显的减少。但是云内液态水含量仍然变化不大,最大值和最小值差别不到2 %。
4 结论综合所作的数值模拟结果,可以得出以下结论:
(1) 云中液态水含量基本上与海盐和硫酸盐粒子的初始数目,以及云团上升速度无关。由于液态水含量不变,随着云滴数目的增加,云滴的有效半径减小。
(2) 硫酸盐气溶胶粒子对云微物理特性影响很大,云滴数目与硫酸盐数目同量级。水汽相对充足情况下,数值上差别也不大;水汽相对不足情况下,由于水汽供应不足,云滴数目与硫酸盐数目有明显差异。海盐粒子的作用分两个方面:一方面海盐粒子半径较大,在云的形成初期可以活化为云滴;另一方面由于海盐的活化吸收水汽导致过饱和比的降低,又会减少由硫酸盐活化成的水滴数目。
(3) 硫酸盐粒子浓度特别强的情况下 (人类活动污染比较严重时),如果水汽相对不足,云滴数目会明显小于硫酸盐粒子浓度;海盐粒子的存在,加剧了水汽供应不足,从而可以在很大程度上进一步降低云滴数目。这一结果表明如果不考虑海盐气溶胶粒子,硫酸盐气溶胶粒子对于云微物理特性的影响可能被大大高估了。
| [1] | Pruppacher H R, Klett J D. Microphysics of clouds and precipitation. Kluwer academic Publishers, 1997. |
| [2] | Shaw G, Bio-controlled thermostasis involving the sulfur cycle. Climate Change, 1983, 5: 297–303. DOI:10.1007/BF02423524 |
| [3] | Nguyen B C, Bonsang B, Gaudry A, The role of the ocean in the global atmospheric sulfur cycle. J Geophys Res, 1983, 88: 10903–10914. DOI:10.1029/JC088iC15p10903 |
| [4] | Charlson R J, Lovelock J E, Andreae M O, et al. Oceanic phytoplankton, atmospheric sulphur, cloud albedo and climate. Nature, 1987, 326: 655–661. DOI:10.1038/326655a0 |
| [5] | Twomey S. Atmospheric Aerosols. Elsevier Scientific Publishing Company, 1973. |
| [6] | Charlson R J, Schwartz S E, Hales J M, et al. Climate forceing by anthropogenic aerosol. Science, 1992, 255: 423–430. DOI:10.1126/science.255.5043.423 |
| [7] | Langner J, Rodhe H, Crutzen P J, et al. Anthropogenic influence on the distribution of tropospheric sulfate aerosol. Nature, 1992, 359: 712–716. DOI:10.1038/359712a0 |
| [8] | Ghan S J, Guzman G, Razzak H A, Competition between sea salt and sulfate particles as cloud condensation nuclei. J Atmos Sci, 1998, 55: 3340–3347. DOI:10.1175/1520-0469(1998)055<3340:CBSSAS>2.0.CO;2 |
| [9] | 肖辉, 徐华英, 黄美元. 积云中云滴谱形成的数值模拟研究 (一). 大气科学, 1988, 12, (2): 121–130. |
| [10] | 肖辉, 徐华英, 黄美元. 积云中云滴谱形成的数值模拟研究 (二). 大气科学, 1988, 12, (3): 312–319. |
| [11] | 许焕斌, 段英. 云粒子谱演化研究中的一些问题. 气象学报, 1999, 57, (4): 450–460. |
| [12] | 赵春生.遥远海洋边界层中大气气溶胶的产生和演变过程的数值模拟:[博士学位论文].北京:北京大学, 1996. |
| [13] | 赵春生, 张代洲, 秦瑜. 海洋边界层中大气气溶胶产生和发展的数值模拟. 自然科学进展, 1998, 8, (4): 440–447. |
| [14] | O'Dowd C D, Jennings S G, Submicron particle, radon, and soot carbon characteristics over the northeast Atlantic. J Geophys Res, 1993, 98: 1123–1135. DOI:10.1029/92JD02387 |
| [15] | O'Dowd C D, Consterdine I E, Lowe J A, Marine aerosol, sea-salt, and the marine sulfur cycle:A short review. Atmos Environ, 1997, 31: 73–80. DOI:10.1016/S1352-2310(96)00106-9 |
| [16] | Hobbs P V. Aerosol-Cloud-Climate Interactions. Academic Press INC, 1993. |
2005, 16 (4): 417-425




