尹金方(1984-),男,贵州盘县人,助理研究员,主要从事云—降水微物理和云微物理参数化方案研究. E-mail:yinjf@cams.cma.gov.cn
云微物理参数化方案在数值模式中起着重要的作用,是影响数值天气预报和气候预测准确性的最大因素。系统回顾了中尺度数值模式中云微物理参数化方案的研究进展,并统计分析了最近十余年云微物理参数化方案在中国范围内的敏感性试验研究成果。Lin方案和RutledgeHobbs方案奠定了中尺度模式中云微物理参数化方案的基础,其他方案都是直接或间接在这2个方案的基础上从多方面改进而形成的。这些改进主要体现在:①水凝物粒子分类数目;②冰核活化;③粒子谱分布描述函数;④粒子谱截距的取值;⑤粒子间相互转换阈值大小的设定。中国范围内云微物理参数化方案敏感性试验研究成果统计表明,使用WRF模式中Lin方案的模拟效果较好,MM5模式采用Goddard和Reisner方案效果较好。
: Cloud microphysical scheme plays a significant role in cloud and precipitation simulation, and it is however one of the main error originations in predictions. An overview of the progresses in the cloud microphysical parameterization in mesoscale models has been done, and a statistical analysis of the sensitivity experiment results on the microphysical parameterization schemes over China have been performed. The Lin and RutledgeHobbs schemes lay a solid foundation for the development of cloud microphysical parameterization, and all later developed schemes were directly or indirectly based on the Lin and RutledgeHobbs schemes with some improvements. These improvements mainly include:①number of hydrometeor classicication; ②ice nucleation; ③hydrometeor particle size distributions; ④given values of intercept for the hydrometeor particle size distributions, and ⑤the threshold values for autoconversion processes. According to the statistical results, the Lin scheme in the Weather Research and Forecasting (WRF) model performed well in simulations, and the Goddard and Reisner schemes in the PSU/NCAR Mesoscale Model (MM5) simulated well.
云在天气和气候的中起着重要的作用。云中水汽和各种水凝物之间的相互转换,以及它们的动力、热力效应直接影响天气过程(如,雷暴,冰雹,飑线等)的发生发展。其次,感热、潜热和动量输送等反作用于大尺度环境,直接影响大气温度、湿度场的垂直分布。因此,合理的描述云物理过程及云与大尺度环境之间的相互作用是准确模拟各种天气过程的关键之一。早在1970年代,许多学者在数值模式中引入了较复杂的云微物理过程[ 1, 2]。特别是Lin[ 3]方案和Rutledge-Hobbs[ 4, 5]方案(简称RH方案)提出之后,云微物理参数化方案得到了极大的发展,形成了多种各具特色的云微物理参数化方案[ 6~ 9]。
IPCC[ 10]明确地指出云物理过程及其反馈机制是导致数值模式预报不确定性的最大因素。不同的云微物理参数化方案“再现”不同类型天气过程的能力存在较大的差别,即使各方案均能够合理地“再现”出大气环流形势的演变、降水分布和强度,但对云中微物理过程的模拟能力却存在着较大的差别[ 11]。孙晶等[ 12]认为MM5模式中的Reisner[ 13]霰方案对华南暴雨具有较好的模拟能力,而饶纲伟等[ 14]则认为Goddard方案[ 15]对粤北冻雨天气的模拟优于Reisner方案。闫之辉等[ 16]对比分析指出使用WRF模式中Lin方案的总体预报结果较好,但不同的分辨率各方案表现出较大的差异[ 17]。总之,当前中尺度模式中微物理参数化方案各具优势和局限性,在选取微物理参数化方案时需要结合关注研究对象和研究目的。
不同云微物理参数化方案的描述方法存在较大的差异。不仅粒子分类数目、密度、分布谱、下落末速度等方面不同,而且在具体的云微物理过程的物理描述和处理方式上也存在较大的差异。Kessler[ 18]方案仅包含云水和雨水两种水凝物粒子,仅考虑云滴凝结、云雨自动转化、云滴同雨滴的碰并等基本过程;Dudhia[ 19]方案中引入了简单冰相过程;Tao 等[ 15]增加了饱和调整过程;McCumber等[ 20]则将雹或霰当作单独的一类冰相粒子;而Meyers[ 7]方案和Morrison[ 21]方案则采用Gamma函数描述粒子谱分布。除了各方案本身的差别之外,各方案关注的天气过程也不尽相同。Rutledge等[ 4]关注大尺度抬升和中尺度辐合造成的锋面降水;Morrison 等[ 8]强调云凝结核和冰核的重要性及其对长短波辐射的影响;Thompson等[ 22]的主要目标是改善对冻雨事件的预报,提供较好的航空危险天气预警。
许多研究意识到了云微物理参数方案的重要性。一些学者[ 23, 24]曾对中尺度模式中的云微物理参数化方案作了相应的介绍,但都未详细地分析不同参数化方案的联系和异同,也未对云微物理参数化方案在中国的使用效果做调研分析。尽管国外在这方面有较多的研究成果[ 25, 26],但是由于东亚季风区背景下的云微物理过程具有独特性[ 27],加之青藏高原大地形的影响,中尺度数值模式中的云微物理参数化方案是否适用于中国地区天气过程,有必要进一步的分析。
本文选取中尺度数值模式MM5,WRFκ,ARPSλ和GRAPESμ ( The PSU/NCAR Mesoscale Model.κ Weather Research and Forecasting Model.λ Advanced Regional Prediction System.μ Global/Regional Assimilation and Prediction System.) 中包含的云微物理参数化方案,详细地分析了云微物理参数化方案之间的联系和异同,以及它们主要适用的天气过程,并调研了最近十余年来中尺度模式云微物理参数化方案在中国范围内敏感性试验研究成果。以期充分认识不同云微物理参数化方案之间的异同,为参数化方案研究提供参考;同时为模拟研究中具体的天气过程(如,台风、华南暴雨、梅雨、东北冷涡等)选取合适的云微物理参数化方案提供依据。
根据云的相态将云微物理方案分为暖云方案和混合相云方案。暖云方案只包含云水和雨水两类水凝物,通常只考虑云滴凝结、云雨自动转化、云滴同雨滴的碰并、雨滴凝结蒸发和重力沉降等物理过程。混合相云方案包括云水、雨水、冰晶、雪晶、霰、雹等水凝物粒子,细致地考虑了冰相粒子的凝华、冻结、融化和相互碰并等过程。根据混合相云方案中冰相粒子数目的不同、冰相过程复杂程度等,混合相云方案可分为简单冰相方案和复杂混合相云方案两大类。
暖云是指由液相粒子构成的云,一般处于温度高于0℃的环境内。暖云方案假定冰相过程对云的热力和降水过程都不起重要作用,描述云水通过水汽达到饱和后的凝结过程产生,当云水含水量达到一定阈值,云水向雨水的转化,初生雨滴下落过程中不断碰并云滴增长,直至降落到地面的整个过程。暖云方案体现了最基本和最主要的云微物理过程,能够描述成云致雨的基本特征和主要的热力、动力反馈,有时推广应用于混合云降水过程。暖云方案的预报量一般包括云水和雨水比质量两种,通常考虑云滴凝结、云雨自动转化、雨滴碰并云滴、雨滴凝结蒸发和重力沉降等过程,复杂一些的还会考虑雨滴的自动破碎等。
自然云内温度通常在0℃以下,含有冰相粒子(冰晶、雪晶、霰和雹)。特别是在中高纬地区,云中包含大量冰相粒子,因此参数化方案中必须考虑云中的冰相过程。相对于暖云方案,简单冰相方案中增加了冰相粒子,其中云冰和雪晶是同一类水凝物,云水和雨水是另一类水凝物。当温度高于0℃时为云水和雨水,反之为云冰和雪晶。当云冰和雪晶下落通过0℃高度层时,云冰和雪晶由于融化而向云水和雨水转化;反之,当云水和雨水上升穿过0℃高度层后,云水和雨水通过凝华、冻结等过程向云冰和雪晶转化。
简单冰相方案抓住了自然云中的主要云微物理过程,并且节省计算资源,在气候模式和早期实时业务中尺度模式中应用比较广泛。简单冰相方案的主要缺点是没有考虑过冷水的存在区域,相变潜热导致冻结层附近可能出现温度反馈的不连续,并且融化过程只在冻结层的相邻层发生。其次,当液态水处于0 ℃以下环境时全部冻结,这与液态水在-40 ℃以下才全部自发冻结有较大的差别。再次,云冰和雪晶下落通过0℃层时立即全部融化,这种假定与实际不符,这正是导致模式计算过程中0 ℃层附近层出现“振荡”的主要原因。此外,简单冰相方案导致云滴浓度明显偏高,从而引起较大的辐射差异[ 28]。
Lin方案在云微物理参数化方案中占有重要的地位,对云微物理参数化方案的发展给予了极大的促进作用。Lin方案包含云水、雨水、冰晶、雪晶和雹五类水凝物粒子。Rutledge 等[ 5]在Lin方案的基础上发展了另一套参数化方案,将雹替换为霰粒子。为了节省计算资源,Dudhia[ 19]基于Lin方案和RH方案建立了一套简单混合相云参数化方案。Hong 等[ 29]对RH方案和Dudhia方案进行了改进,包括:①物理过程开始计算之前,格点上云水(冰)和雨(雪)的负值设为零;②引入了时步分离积分法;③改进了高层云量的算法。此后,Hong等[ 30]做了进一步改进,将雪晶谱截距设为温度的函数,而不是定值。其次,增加了冰晶重力沉降作用的新描述方法,改进了冰晶质量依赖于数浓度等冰晶过程,即WRF模式中的WSM3方案。Lim等[ 31]在WSM3方案的基础上增加了云滴和雨滴数浓度的预报方程,形成了WDM3双参方案。
随着云微物理理论和观测研究的不断深入,人们逐渐认识到自然云中的成冰机制及冰相过程远远复杂于简单冰相方案描述的过程。简单的冰相方案不能充分反映出冰相粒子动力和热力效应及其转化机制。鉴于简单冰相过程的不足,加之计算机计算能力的提高,冰相过程不断被细化,一些方案甚至包含高达110余种云微物理过程[ 32]。
2.3.1 冰相粒子种数
由于冰相粒子的形状和属性各异,与温度、湿度等气象因子的关系又不确定,冰相云过程的处理远比液相云困难。冰相粒子的下落末速度、碰并系数也有很大的不确定性,加之冰相粒子形状、大小各异,以及冰晶繁生等一系列不确定的因素,使得冷云中云微物理过程变得极其复杂。
由于不同冰相粒子下落末速度之间明显的差异给微物理参数化方案中带来了较大的困难,因此如何考虑冰相粒子的分类是云微物理参数化方案的关键科学问题之一[ 33]。一方面,为了增加参数化方案的灵活性和普适性,需要不断的细化粒子的种类。然而,不断细化粒子带来新的问题:①粒子的物理属性难以界定,②粒子之间相互转化的描述问题和随之而来的海量计算。因此,合理的选择粒子分类数目,是云微物理参数化方案中考虑的首要问题。
当前已有20余种冰相云微物理参数化方案。Ogura等[ 34]、Koenig等[ 35]和Orville等[ 36]等方案考虑一类或两类冰相粒子;Lin等[ 3]包含云冰、雪晶和雹三类冰相粒子;Rutledge等[ 4, 5]则将将雹粒子替换为霰粒子;其他一些方案[ 7, 37]同时引入了霰粒子和雹粒子。一些方案将所有的冰粒子统称为一类(如, Schultz[ 38]; Lin等[ 9])。相反,另一些方案则根据冰相粒子的形状、密度、属性等将冰相不断细化,如Straka等[ 33]将冰相粒子分为3类冰粒子,雪团,三类霰粒子,冻滴,两类冰雹,共计10种冰相粒子。McCumber等[ 20]研究指出,合理的模拟云中的冰相过程至少需要包含四类冰相水凝物粒子。当前参数化方案中通常将水凝物粒子分为云水、雨水、云冰、雪晶、霰和雹,其中冰晶、雪晶、霰和雹为冰相粒子。
2.3.2 冰核活化
“贝杰龙”过程在冰相云中对降水的形成有着重要的作用[ 39]。自然云通常是固、液相混合云,冰晶是形成其他冰相粒子的“胚胎”,在冰相过程中起着重要的作用。因此,合理参数化冰核活化过程尤为重要。当前冰核活化有不同的参数化方法,Fletcher[ 40]认为冰核浓度是温度的函数,Cooper[ 41]指出冰核浓度是饱和度的函数,胡志晋等[ 42]、Thompson等[ 22]、许焕斌等[ 43]等同时考虑饱和度和温度对冰核活化的影响。Yin等[ 44]统计分析指出中国北方地区的冰核浓度与其他地区的冰核特征有着明显的差别。因此,在选取方案时,需要根据实际观测结果调整参数化方案中参数的取值。
2.3.3 云微物理过程的发展
简单冰相方案的优点是极大地减少了计算量,在计算条件有限的条件下,能够利用简单冰相混合方案进行高分辨率的数值模拟。随着计算机计算能力的不断提高,无论是暖云方案,还是简单混合相云方案,都难以全面反映出真实大气中的云微物理过程。为了在数值模式中更加详细、合理地描述云微物理过程,需要在数值模式中引入复杂混合相云微物理参数化方案。
Lin方案为复杂混合相云方案的发展奠定了坚实的基础。该方案的主要优点是引入了雪晶,雪晶的形成主要通过大尺度冰晶与云(雨)的碰冻,凝华和淞附等方式。雪晶的增加导致了云水含量的减少,使得在云生命史早期雨水的形成减少。Tao等[ 15]在Lin方案的基础上增加了饱和调整方案,形成了Goddard方案;Chen等[ 45]基于Lin方案和RH方案,结合Goddard方案的饱和调整,构建出了一套新的复杂混合相云方案,即当前WRF模式中的Purden Lin方案。
当前WRF模式中的Goddard方案是经过几次调整后的方案。首先,McCumber等[ 20]增加了将雹或霰当作第三类冰物质的开关,霰粒子主要用于热带对流系统降水,而雹粒子适合于中高纬深对流天气系统降水,且指出这种分类依据主要是霰和雹的粒子分布谱存在明显的差异。霰和雹的这种差异不仅影响水物质数量和层状云砧区域的描述,还与云微物理-动力-辐射过程等密切相关。其次,增加了新的饱和调整方案,确保格点上晴空(多云)时不存在过饱和(未饱和)现象。第三,除融化、蒸发和升华之外的过程,其它过程都是基于格点热力学方程来计算,确保每一个过程都用同样的方式来处理。第四,每一种粒子降落到地面的总质量不能超过所有粒子质量总和,确保云微物理计算过程中水物质守恒。此外,Goddard方案提供仅包括冰和雪两类冰粒子的选项,通常只用于低分辨(格距大于5 km)大尺度天气系统模拟中,如冬季对流系统和锋面对流系统。Gilmore等[ 46]认为霰粒子谱分布的截距和斜率变化非常大,将截距固定为常数在很大程度上歪曲了物理过程,最近NSSL单参方案[ 47]采用相同的处理方法。
2.3.4 粒子谱的演变
采用不同函数描述水凝物粒子谱,不但粒子的分布形式不同,而且粒子谱参数的处理上也有很大的差异。早期绝大部分参数化方案假定粒子谱服从指数分布。其表达式为
N( D)= N0exp(- λD) , (1)
式中参数 N0和 λ分别为截距和斜率。通常固定 N0的值, λ的值由水凝物的比质量确定。谱型的演变仅由 λ来确定,称之为单参(One-moment)谱演变。1970年代,Srivastava[ 48]提出了双参数演变,但在实际参数化的过程中,往往假定 N0和 λ的比值为常数,这相当于比浓度不变,本质上为固定一个参数的前提下讨论粒子谱的演变[ 49]。
观测研究表明,Gamma(Γ)函数能够更好地描述云中水凝物粒子谱。Γ函数可表示为
N(D)=N0D μexp(-λD), (2)
式中 μ为谱形状参数。谱参数 N0和 λ可由比质量和比浓度方程联立求解得出,粒子谱可由 N0和 λ两个参数共同决定,称为双参(Two-moment)谱演变。Meyers等(Meyers方案)用Γ函数来描述水凝物粒子的分布被广泛的应用的云微物理参数化方案中[ 7]。Thompson等[ 22]以Reisner方案为基本框架,允许过冷水和冰晶同时存在低于0 ℃的环境中,并根据实际观测结果改进参数化方案。Thompson方案与其他方案的主要区别为:①假定雪晶谱分布依赖于冰水含量和温度;②采用指数函数和Γ函数来分段描述雪晶谱;③假定雪晶不是圆形的,且密度是最大维长度的函数。
Γ函数需要确定三个参数,但实际参数化过程中只能根据比质量和数浓度预报方程求解出 N0和 λ。因此, μ的取值常设取为定值。NSSL双参方案[ 47]中,提供了设置参数 μ的选项,模式使用者可以根据实际需要设置粒子谱参数 μ的取值。观测指出[ 27], μ并不是一个常数。Milbrandt等[ 37]为了改进双参谱演变中存在的不足,增加了雷达反射率的预报方程,实现了三参数演变。然而,雷达反射率本身不是独立的预报量,而是依赖于粒子比浓度和比质量的衍生变量,并不是真正意义上的三参演变。粒子谱的分布需要 N0, λ和 μ决定,预报这3个参数的参数化方案,称为三参(Three-moment)方案。然而,三参方案的缺点是求解谱参数时需要3个方程联立求解,而实际参数化过程中只有比质量和比浓度两个预报方程,还需要一个方程才能使方程组闭合。因此,在实际参数化过程中谱型参数往往根据观测结果设为常数,或诊断给出。在今后的参数化研究中,如何求解谱型参数将是参数化发展的重点之一。
![]() | 图1 中尺度模式中混合相云微物理参数化方案之间的相互联系Fig.1 Relationship among different bulk microphysical schemes in GRAPES, WRF, MM5, and ARPS models |
尽管不同的云微物理方案有所差别,但是这些方案之间存在着紧密的联系,有着共同的物理基础。图1显示的是常用中尺度模式GRAPES,WRF,MM5和ARPS中云微物理参数化方案之间的联系。可以看出,Lin方案和RH方案奠定了云微物理参数化方案的基础,其他方案是直接或间接在Lin方案的基础上,从不同的方向做了相应的改进,从而形成了各具特点的方案。
2.3.5 粒子谱截距的取值
单参方案中粒子谱截距取为定值,仅通过斜率的变化来反映谱变化。由于当前使用的微物理参数化方案绝大多数来源于国外,其参数的取值主要采用其它地区的观测结果,不同云微物理方案雪晶谱截距取值也不同(表1)。Yin等[ 27]对比东亚和其他地区粒子谱发现,东亚地区雪晶谱的截距比其他地区小,斜率也比其他地区略小;雨滴谱也类似。因此,在数值模拟研究中,可以根据研究对象的地理分布、天气过程特征等因素结合已有的观测结果对截距的值适当地调整。
![]() | 表1 不同云微物理方案雪晶谱取值 Table 1 Values of intercept for snow crystal size distribution |
2.3.6 粒子间相互转换阈值
Kessler型云水向雨水自动转化参数化过程中阈值的选取尤为重要,因为阈值的大小关系到自动转化的早迟。将阈值设定为常数的缺点是当空气中云水含量较低时,无云水向雨水的转化,从而低估了云的降水效率;相反,在云水含量较高时,大量云水瞬间转化为降水,而实际仅有部分含水量转换为降水,这样就高估了云的降水可能。因此,不同的方案中往往根据需要对阈值进行调整。表2中列举出了不同微物理参数化方案中云水向雨水自动转化和雪晶自动转换的阈值。可以看出,不同方案中阈值的大小存在一定的差异,这些阈值的差异很可能导致云中各微物理变量含量的明显不同。在数值模拟研究中,可以根据研究的地理范围,天气过程的不同对阈值的设定进行适当的调整,以便得到更为理想的模拟结果。
![]() | 表2 不同云微物理方案中自动转换过程中的阈值 Table 2 Threshold values for autoconversion processes in different schemes |
为了提高数值模式对东亚地区云和降水的预报能力,国内许多学者对云微物理参数化方案进行了系列研究。胡志晋等[ 52]改进了Kessler暖云方案,将云滴形成雨滴的转换过程分成酝酿和碰并增长两个阶段。许焕斌等[ 53]建立了包含雨滴和雹粒子的双参冰雹云模式。胡志晋等[ 54]在已有工作的基础上构建了对流云和层状云双参数方案。刘公波等[ 55]在前人的基础上提出了一个混合相层状云模式。刘奇俊等[ 56]吸收并改进了胡志晋[ 42]的层云和积雨云微物理方案,成功将双参数暖云方案、简化混合相云方案和双参数混合相云方案移植到国家气象中心业务有限区数值预报模式HLAFS中,提高了原模式预报系统的湿物理过程和对云-降水过程的预报能力。楼小凤[ 57]以胡志晋等[ 42, 54]早期的对流云模式和层状云模式为基础发展了一套混合相云双参数微物理方案,命名为CAMS复杂云微物理方案;孙晶等[ 58]将CAMS复杂云微物理方案耦合到GRAPES模式中;Gao等[ 59]进一步改进了CAMS复杂微物理参数化方案,耦合到WRF模式中。尹金方[ 60]基于东亚地区长时间云和降水微物理的观测结果,建立了一套较为复杂的混合相云双参数微物理方案。为了研究冰雹的形成机制,洪延超[ 61]将冻滴从霰粒子中分离出来,并建立相应的预报方程,形成了包含云滴、雨滴、冰晶、雪、霰、冻滴和冰雹7 种水凝物粒子的参数化方案。为了更好地描述雨滴谱,郭学良等[ 62]将体积水方案与分档法结合,对雨滴谱进行分档处理。
为了适应不同条件下对云降水的模拟研究,一些学者对中尺度模式中的云微物理方案作了相应的改进。孔凡铀等[ 63]综合多个云微物理参数化方案,构建出较为复杂的云微物理参数化方案。赵震等[ 64]在MM5中Reisner2方案的基础上加入了新的云微物理过程。刘卫国等[ 65]增加了冰晶数浓度、人工冰晶数浓度及比水量到ARPS模式云微物理方案中,使该方案的微物理预报量增加到9个。
在过去的30余年里,中国对云微物理体积水参数化方案的研究取得了一系列成果。这些参数化方案的改进有效地促进了对中国地区云和降水的数值模拟研究。当然,国内还有其它许多关于云微物理参数化研究的成果,由于本文只追踪了国内外主要气象相关期刊,难免会有遗漏。值得注意的是,国内研发和改进的这些参数化方案中,参数的取值几乎采用国外的观测结果,其在中国地区的实用性还需进一步验证。随着中国外场试验的广泛开展,有必要针对中国地区的观测结果对模式中的云微物理参数化方案进行修订。
云微物理参数化方案对云和降水的模拟结果有着重要的影响,就如何选取合适的物理参数化方案是有重要意义的科学问题。世界各地进行了大量云微物理参数化方案对降水的影响的研究[ 20, 66, 67],目的是为了能够选出一种理想的云微物理参数化方案。
表3列举出了中国区域内2000—2012年针对中尺度模式中云微物理参数化方案敏感性试验研究的结果。总体而言,WRF模式中 Lin方案的效果较好;MM5模式中,Goddard方案和Reisner方案效果较好。需要注意的是WRF和MM5模式中云微物理参数化方案的对比分析研究较多[ 68],而对GRAPES和ARPS模式中云微物理方案的对比分析研究较少,其原因是GRAPES模式发展时间较晚[ 69],ARPS模式早期版本中云微物理参数化方案使用较少。其次,这些研究以单参方案为主,对最近几年发展的双参云微物理参数化方案的研究较少,这也是今后工作重点的一部分。随着计算能力的不断提高和双参云微物理参数化的发展,双参云微物理参数化方案正被逐渐应用到数值模拟研究和业务预报系统中。Zhang等[ 70]研究认为,当网格距在20~50 km之间应该同时使用积云参数化和云微物理参数化方案。由于当前中尺度模式使用的背景场资料分辨率还不够高,还需要同时使用积云参数化和云微物理物理方案。因此,表3中同时列出了与云微物理参数化组合使用的对流参数化方案。
![]() | 表3 2000—2012年期间中国范围内云微物理参数化方案敏感性试验结果 Table 3 Sensitivity experiment results for the cloud microphysical parameterization schemes over China during the period from 2000-2012 |
本文总结了当前中国常用中尺度模式中的主要云微物理参数化方案的特征,归纳分析了不同方案之间的联系和差异。尽管不同的方案有所差别,但有着相同的物理基础,各方案之间有着一定的联系。Lin方案和Rutledge-Hobbs方案奠定了云微物理参数化方案的基础,其他方案都是直接或间接在这两个方案的基础上从不同方面改进而形成的。具体改进可以归纳为以下几方面。
(1)当前中尺度数值模式云微物理参数化方案中,水凝物粒子主要包括云水、雨水、冰晶、雪晶、霰和(或)雹。其中包含霰粒子的方案主要适用于热带地区降水(如,Reisner霰方案、WDM6方案等),而中高纬深对流云降水主要采用含雹粒子的方案(如Lin方案、Goddard方案、Thompson方案等)。
(2)云微物理参数化中通常可用指数函数和Gamma函数来描述粒子谱。许多观测资料统计表明,用Γ函数能够很好地描述云和降水粒子谱。因此,新兴的云微物理参数化方案均采用Gamma函数描述粒子谱。
(3)当前中尺度模式中包含单参和双参云微物理参数化方案。单参方案固定截距,通过斜率的变化来反映粒子谱演变,计算量小,且相对稳定;双参方案理论上更加合理地反映云和降水的物理过程,但当前双参方案粒子数浓度和比质量的预报方程之间缺乏约束,一定程度上导致了双参方案的不确定性。
(4)使用单参方案时,不同方案中粒子谱截距的取值不同。在数值模拟过程中,可以根据已有的观测结果对取值的设定进行适当的调整。此外,云水自动向雨水转换阈值、冰核活化参数的取值应根据观测结果进行调整。
不同云微物理参数化方案在处理具体物理过程中有所差别,在选择云微物理参数化方案时应当对方案有所了解。不同的方案研究的天气现象不同,在选择参数化方案时,应针对研究的具体天气过程,选择合适的参数化方案以得到更加合理的结果。中国范围内2000—2012年云微物理参数化方案敏感性试验研究结果表明,使用WRF模式中Lin方案的模拟效果较好,MM5模式采用Goddard和Reisner方案模拟效果较好。需要注意的是最近的一些研究也表明新Thompson双参方案对于中国区域内强降水落区和量级都表现出较好的预报能力。
随着云微物理观测和理论研究的深入[ 87, 88, 89, 90],进一步改进云微物理参数化方案研究势在必行。结合当前中尺度数值模式中的云微物理参数化方案特征及不足,提出以下建议,供今后云微物理参数化方案构建研究中方案设计作参考。①云中水凝物粒子分类数目是云微物理参数化方案构建中需要考虑的首要问题,需要针对研究的天气过程合理地对进行粒子分类。②冰晶是形成其他冰相粒子的“胚胎”,在冰相过程中起着重要的作用,需要根据大量实际观测结果合理参数化冰核活化过程。③Gamma函数能够较好地描述云中水凝物粒子谱,如何求解谱型参数将是参数化发展的重点之一。④水凝物粒子比质量和数浓度预报方程之间缺乏物理约束,常常导致两者之间的不匹配,今后需要根据大量观测建立比质量和数浓度两者之间约束关系,以改进双参方案的稳定性和不确定性[ 91]。⑤结合CloudSat卫星、双偏振雷达等观测资料,设定云微物理参数化方案中不同垂直高度层上自动转换阈值相同,改进自动转换过程。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[41] |
|
[42] |
|
[43] |
|
[44] |
|
[45] |
|
[46] |
|
[47] |
|
[48] |
|
[49] |
|
[50] |
|
[51] |
|
[52] |
|
[53] |
|
[54] |
|
[55] |
|
[56] |
|
[57] |
|
[58] |
|
[59] |
|
[60] |
|
[61] |
|
[62] |
|
[63] |
|
[64] |
|
[65] |
|
[66] |
|
[67] |
|
[68] |
|
[69] |
|
[70] |
|
[71] |
|
[72] |
|
[73] |
|
[75] |
|
[76] |
|
[77] |
|
[78] |
|
[79] |
|
[80] |
|
[81] |
|
[82] |
|
[83] |
|
[84] |
|
[85] |
|
[86] |
|
[87] |
|
[88] |
|
[89] |
|
[90] |
|
[91] |
|