Numerical simulation of the scavenging rates of ice crystals of various microphysical characteristics

Numerical models of trajectories of small aerosol spheres relative to oblate spheroids were used to determine ice crystal scavenging efficiencies. The models included the effects of aerodynamic flow about the ice particle, gravity, aerosol particle inertia and drag and electrostatic effects. Two electric configurations of the ice particle were investi-gated in detail. The first applied a net charge to the ice particle, of magnitude equal to the mean thunderstorm charge distribution, while the second applied a charge distribution, with no net charge, to the ice particle to model the electric multipole charge distribution. The results show that growing ice crystals with electric multipoles are better scavengers than single ice crystals with net thunderstorm charges, especially in the Greenfield gap (0.1 to 1.0 μm), and that larger single crystals are better scavengers than smaller single crystals. The results also show that the low density ice crystals are more effective scavengers with net charges than they are with charge distribution.     

Calculation of Ice Crystal Diffraction

The diffraction angular distributions are calculated for different ice crystals (e.g,solid and hollow hexagonal columns,solid plates and single hexagonal bullet) randomly oriented in space,based on the theory of Cai and Liou.Results of solid column and plate are compared with those from formulas of Wendling et al and Liou.The resull comparison shows that all three procedures discussed in this paper have good agreement.After considering the coin putation time and comparison results,it is suggested that the formula of Wendling et al.can be used to calculate the solid column diffraction,and the formula of Liou can be used to calculate the plate diffraction.The comparison also shows that the diffraction results of solid columu,hollow column and single bullet ice particles randomly oriented in space are very close.     

Calculation on the light scattering function of hexagonal ice crystals

Hexagonal ice crystal is a basic form of the particles in cirrus. The emphasis in this paper is on discussing the computational results of scattering phase matrices for infrared absorptive band by making use of the model developed by the authors. Comparisons have been made between infrared and visible light as well as between hexagonal columns and plates and equivalent spheriods and spheres. It is found that the effects of light absorption on phase function may cause the decreases of scattering energy in side direction and the directions near backscattering, and weaken halo and rainbow peaks which are induced by the refractions and internal reflections. There is a general agreement in the phase functions for ice crystals with different shapes. However, the scattering intensity in side direction for spheres and in backward and backside directions for spheroids is weaker as compared with hexagons. It is believed that these features would, doubtless, have influences on the detection of cirrus by laser and radiation transfer in the atmosphere.     

A general hydrodynamic theory for mixed-phase microphysics. Part I: drag and fall speed of hydrometeors

A general parameterization for solid and liquid hydrometeors is presented. hydrometeors basically are viewed as porous spheroids with the following variable parameters: diameter, axial ratio, mass, and porosity. Based on this parameterization a functional dependence on the Reynolds number of the drag of hydrometeors is derived, which is based on boundary layer theory. The basic form of this functional dependence is consistent with viscous theory and the inertial drag at low Reynolds numbers is predicted with good accuracy by matching the results from the boundary layer theory with Oseen's theory of creeping motion. Based on this solution a general semi-empirical expression for the Reynolds number and fall speed of particles is found. The results from the present theory are in remarkable agreement with experiments: The errors generally are < 5–10% for a wide variety of hydrometeors in the range of Reynolds numbers 0<N_Re<5×10⁵, including columnar and variously branched planar ice crystals, rimed and unrimed aggregates, lump, conical, and hexagonal graupel, hail, and rain drops. The present parameterisation aims far beyond the limits of the conventional methods since it is suitable for mixed-phase models of the microphysics of precipitation with continuously varying particle mass and shape characteristics and including processes such as depositional growth of ice crystals under varying environmental conditions, collisional growth of particles, and melting.     

不同环境温度下冰云粒子的单次散射参量

本文提出了直接Stokes参量变换的射线追踪方案,优化了作者过去建立的求解六角冰晶粒子散射的射线追踪模式.在冰晶有吸收(折射指数为复数)的情况下,新模式采用了有效折射指数处理粒子内的非单色折射波.其次,新模式还引入了有效单次散射参量,消除了来自几何透射线的前向δ散射效应.最后,根据不同环境温度下的冰云中冰晶粒子的粒径谱分布,本文计算了粒径谱积分平均的3D随机取向六角冰晶粒子的单次散射参量.文中给出的冰晶粒子单次散射参量的数值,对研究冰云辐射特性将有一定参考价值.     

SCATTERING PHASE MATRICES OF ICE CRYSTALS WITH HEXAGONAL PRISM AND TRIANGULAR PYRAMID FORM—A VECTOR RAY TRACING METHOD

The work presented previously by the authors(Cai and Liou,1982)has been extended in this paper.By making use of our improved model the calculations on scattering phase matrices of hexagonal prismice crystals(HPIC)have been conducted for monodisperse and polydisperse systems.Compared with themodel of Cai and Liou,the required computational quantity is decreased by about two orders of magni-tude and the errors of results are less for the new model.Meanwhile,the scattering phase matrices of triangularpyramid ice crystals(TPIC)are also computed in the paper,and the comparison between the scatterings of thetwo forms of ice crystals is performed.     

Single scattering from frozen hydrometeors at microwave frequencies

Scattering of electromagnetic waves from homogeneous or coated spheres can be computed in a mathematically exact way using the Mie theory. Therefore, for many approaches in remote sensing, frozen hydrometeors are parameterized as ice spheres. However, many frozen hydrometeors have non-spherical overall shapes and lack a spherically symmetric internal structure. They exist in a huge variety of shapes and exhibit different mixtures of ice, water and air. Therefore it is desirable to accurately compute scattering from non-spherical particles in order to clearly understand the effect the shape of a hydrometeor has on its scattering pattern.In this study, single scattering parameters like scattering cross section, absorption cross section, and asymmetry factor were calculated for frozen hydrometeors using the Discrete Dipole Approximation (DDA). The particles were modeled as hexagonal plates, columns, needles and dendrites by applying known dimensional relationships. The calculations were carried out over a wide range of centimeter and millimeter-wavelengths (1 GHz to 300 GHz), since millimeter-wave radiometers are highly sensitive to scattering by frozen hydrometeors in the atmosphere.The study results show that for size parameters < 1 (a ratio between wavelength and particle size) the scattering cross section of randomly orientated ice crystals is close to that of an equal volume ice sphere. Absorption cross section and asymmetry factor of non-spherical particles however are up to twice as high as that of equal volume ice spheres. Further the influence of the assumed model for the refractive index of ice at microwave wavelengths on the scattering parameters is investigated.     

冰晶产生的非催化触发机制研究

本文研究了用压缩空气爆破(薄膜爆破,激波管等)和火药爆炸非催化方法触发冰晶产生的作用机制。我们通过大量实验认为超声波和冲击波的作用不一定能够触发过冷水滴冻结,而超声速气流的绝热膨胀冷却作用则是触发冰晶产生的确切机制。然后设计了超声速气流发生装置并在过冷雾中进行触发冰晶产生的实验。结果表明,超声速气流(M=1.1)触发冰晶产生非常有效。阈温为-2℃。冰晶产生的临界压强为2.1atm,阈压随实验温度升高而增加,当总压超过阈压0.2—0.3atm时,冰晶大量产生。这些结果与超声速绝热膨胀理论相一致。     

六棱柱和三棱锥冰晶的散射相矩阵——一种矢量射线追踪方法

本文发展了作者过去的工作(Cai-Liou),利用改进后的数值模式对六角柱冰晶(匀谱和非匀谱)的散射相矩阵进行了计算,所需计算量比Cai-Liou的约小二个量级,且本文的数值模式的结果误差较小。同时,对三棱锥冰晶的散射相矩阵,本文也进行了计算,并与六角冰晶的结果进行了比较。     

Laboratory studies about the interaction of ammonia with ice crystals at temperatures between 0 and −20°C

In laboratory experiments the interactions of ammonia with ice crystals were studied within the temperature range between 0 and −20°C. In a first series of experiments dendritic ice crystals were grown from water vapor in presence of ammonia gas in various concentrations between 4 and 400 ppbv. In a second series of experiments pure ice crystals were exposed to a humidified ammonia–air mixture inside a horizontal flow tube. The influence of temperature, ammonia gas concentration (0.6, 1.5, and 10 ppmv), exposure time, and the presence of impurities such as sulfate on the ammonia uptake by the ice surface was investigated by determining the ammonium content in the melt water of the ice crystals by ion chromatography. During the growth of ice crystals significant amounts of ammonia (around 200 μg/l) were taken up even at small gas concentrations. In contrast, even at high gas concentrations the uptake of ammonia by non-growing ice crystals was lower by approximately one order of magnitude. The presence of sulfate on the ice surface affected an enhanced uptake of ammonia by a factor of 5–10. A model is presented which describes the uptake of ammonia by ice considering the chemical processes occurring in the ice surface layer and simultaneous diffusion of ammonia into bulk ice. Even the increased uptake of ammonia by growing ice is rather small compared to the uptake by water droplets; thus, the major process for scavenging of ammonia from the atmosphere via the ice phase might not be the direct uptake by ice crystals but the riming involving super-cooled droplets containing ammonia.     

Development of a continuous flow thermal gradient diffusion chamber for ice nucleation studies

A supercooled continuous flow, thermal gradient diffusion chamber has been developed to study the ice nucleating properties of natural or artificial aerosols. The chamber has concentric cylinder geometry with the cylinder axis alignment and airflow vertically downward. Sample airflow is 1 l min⁻¹ and occupies the central 10% of the annular lamina; it is separated from the ice-covered walls by filtered sheath air. The wall temperatures are independently controlled over the range from about −4°C to −25°C, so that the vapor concentration at the location of the sample lamina can be set to a well defined value between ice saturation and a few percent water supersaturation. There is a range of temperature and supersaturation values across the sample region; for lamina center conditions of −15°C and +1% with respect to water, the range is −14.6 to −15.4°C and +0.53 to +1.31%. Errors in temperature control produce variations estimated as ±0.1°C and ±0.23%. Typical sample residence time is about 10 s. Ice crystals which form on active nuclei are detected optically at the outlet end of the chamber. To enhance the size difference between ice crystals and cloud droplets, the downstream 25% of the warm ice wall is covered with a thermally insulating vapor barrier which reduces the vapor concentration to ice saturation at the cold wall temperature, so cloud droplets evaporate.A mathematical model was developed to describe the temperature and vapor fields and to calculate the growth, evaporation, and sedimentation of water and ice particles. At 1% water supersaturation, the model predicts that ice particles will grow to about 5 μm diameter, and cloud droplets will achieve about 1 μm before they reach the evaporation section of the chamber. A different model was developed to describe the steady state airflow profile and location of the sample lamina.Experimental tests of the chamber were performed to characterize the airflow, to assess the ability of the technique to detect silver iodide ice nucleating aerosols and to distinguish ice crystals from water droplets.     

一次积-层混合云系垂直结构和降水机制的飞机观测资料分析与数值模拟

积-层混合云是影响北京地区的重要降水云系,运用飞机探测资料结合中尺度数值模式WRF,对2014年9月23日发生在北京地区的一次积-层混合云系的垂直结构和降水机制进行了探测资料分析和数值模拟研究。通过分析云系的雷达回波演变,发现云中的对流泡没有出现爆发式增长,回波在垂直方向上增长不明显,此次过程属于积-层水平混合型云系降水。飞机探测资料分析显示,上、下午探测云系的液态水含量都不高（最大低于1 g/m3）;在云系不同高度,飞机探测到的冰晶形状主要有板状、针柱状、辐枝状和不规则状,由于云中过冷水含量相对较低,聚合冰晶的数量明显多于凇附冰晶,冰晶的聚合是云中粒子增长的主要过程。对模拟云系垂直微物理结构和降水粒子的源、汇项分析得到:高层,由凝华产生的冰晶和雪晶在过冷水含量较低的环境中不断聚并、长大并下落,云系中霰的含量很低,增大的冰晶和雪晶下落至0℃层附近融化是产生地面降水的主要机制。此外,融化层附近,雨滴捕获云滴不断长大并下降至地面也是地面降水的另一个重要来源。      

青海省秋季一次对流云人工增雨的数值模拟

利用中国气象科学研究院的三维对流云数值模式，模拟了青海省2002年秋季一次对流云过程，分析了青海省秋季对流云云体发展的动力学特点和微物理特征。数值模拟结果表明，青海省秋季的对流云降水几乎全部是由霰的融化形成的，而自然云中冰晶含量少、冰霰自动转化形成霰胚的过程非常微弱。但采用合适的方法催化以后，冰晶的核化、繁生量增加，通过冰霰自动转化过程形成大量霰胚，霰胚再通过其他冷云微物理过程迅速增长，融化成降水。催化后各种微物理机制都比催化前活跃，同时催化改变了云体的动力场分布，在动力过程和微物理过程的相互促进下，使增雨取得了很好的效果。     

Resolving the Effects of Aperture and Volume Restriction of the Flow by Semi-Porous Barriers Using Large-Eddy Simulations

The Regional Atmospheric Modelling System (RAMS)-based Forest Large-Eddy Simulation (RAFLES) model is used to simulate the effects of large rectangular prism-shaped semi-porous barriers of varying densities under neutrally buoyant conditions. RAFLES model resolves flows inside and above forested canopies and other semi-porous barriers, and it accounts for barrier-induced drag on the flow and surface flux exchange between the barrier and the air. Unlike most other models, RAFLES model also accounts for the barrier-induced volume and aperture restriction via a modified version of the cut-cell coordinate system. We explicitly tested the effects of the numerical representation of volume restriction, independent of the effects of the drag, by comparing drag-only simulations (where we prescribed neither volume nor aperture restrictions to the flow), restriction-only simulations (where we prescribed no drag), and control simulations where both drag and volume plus aperture restrictions were included. Previous modelling and empirical work have revealed the development of important areas of increased uplift upwind of forward-facing steps, and recirculation zones downwind of backward-facing steps. Our simulations show that representation of the effects of the volume and aperture restriction due to the presence of semi-porous barriers leads to differences in the strengths and locations of increased-updraft and recirculation zones, and the length and strength of impact and adjustment zones when compared to simulation solutions with a drag-only representation. These are mostly driven by differences to the momentum budget of the streamwise wind velocity by resolved turbulence and pressure gradient fields around the front and back edges of the barrier. We propose that volume plus aperture restriction is an important component of the flow system in semi-porous environments such as forests and cities and should be considered by large-eddy simulation (LES).     

区域冷锋降水微物理机制研究

20世纪利用一维层状云模式对2002年4月4~5日河南省冷锋降水过程进行了模拟。数值模拟结果显示,此次冷锋降水属于冷云降水过程,冷锋前后云中主要以冰相粒子为主,云中水质粒自上而下的空间分布依次为冰晶、雪、云水、霰、雨水。冷锋前后,各种水质粒有着不同的含量及数密度,但形成水质粒的主要微物理过程都表现为:冰晶数密度的增加主要依靠核化、繁生,大部分雪主要靠凝华、撞冻过冷云水和冰晶增长,霰的质量增加主要靠撞冻雪、过冷云水和雪自动转化而来,大部分的雨水是由霰融化而来,因而此次冷锋降水机制表现为“水汽—雪—霰—雨水”。     

对非球形冰晶94 GHz云雷达后向散射和衰减的研究

针对94 GHz毫米波云雷达的数据处理,基于离散偶极子近似法(DDA)计算了几种非球形冰晶的后向散射及衰减效率,探讨了不同冰云模型下冰云的雷达反射率因子(Z_e)和衰减系数(k)及冰水含量(IWC,记作W)的关系。结果表明:(1)形状对冰晶的散射及衰减效率与粒子大小有关。(2)在实际的冰云中,将六角形冰晶和椭圆冰晶看做同体积的球形粒子将低估其衰减和后向散射。将聚合物冰晶看做等体积球形,将高估其后向散射及衰减,子弹花冰晶的后向散射与同体积的球形相比有减小也有增大,但是衰减比同体积球形的小。(3)冰云模型对Z_e-k、Z_e-W关系具有较大影响,假设滴谱相同的条件下,得到了具体冰云模型下对应Z_e-k、Z_e-W关系的系数。这些探讨为我国W波段云雷达的数据处理提供了参考。     

对流云中冰相过程的三维数值模拟——Ⅱ:繁生过程作用

在可压缩三维冷云数值模式里对冰晶繁生过程进行了参数化,显式考虑了过冷滴冻结破裂和Hallett-Mossop两个繁生机制,并对冰晶浓度进行了订正。对暖底孤立积云的模拟实例表明繁生作用对对流系统宏观特征没有显著改变,地面总降水略有增加,最大雨强增大且时间提前。繁生过程的主要作用是改变了云的微物理结构特征,使云中冰化时间提早,冰化程度大大加强,引起云上部增暖并抬高云顶。在所考虑的两个繁生机制中,冻结破裂繁生机制是引起微物理结构发生变化的主要过程。还对各冰相微物理过程的相对作用进行了分析。     

Observations and parametrization of momentum transfer in heterogeneous terrain consisting of regularly spaced obstacles

Tethered balloon data collected over a heterogeneous, semi-rural area of southern England are presented. The terrain within which the observations were taken is dominated by regularly spaced bluff obstacles that form the boundaries to grass fields. The data are compared with previous surface data to assess how representative these fixed-point measurements are of area-averaged values. It is found that there is a vertical height above which the flow approximates to being homogeneous and observations made above it are representative of the area average. The effective roughness length for momentum is then compared with that derived from a model describing the drag partition around isolated obstacles that are either non-porous or porous to the incident flow. It is shown that the method chosen to partition the drag is of the right order to produce reasonable predictions of the effective roughness length for different roughness densities.     

1960年以来东亚季风区云-降水微物理的直接观测研究

云-降水的直接观测结果是云微物理参数化的重要依据。自1960年以来,处于东亚季风影响下的中国实施了大量对云-降水微物理参数的观测和研究,旨在加深对云-降水微物理过程的认识,从而改进数值模式中云微物理参数化方案和指导人工影响天气作业。云-降水微物理参数包括气溶胶、冰核、云滴、雨滴、冰晶、雪晶、冰雹等粒子浓度和谱分布,以及云滴、雨滴含水量等。中国已有云-降水微物理参数的成果可归纳为:(1)通常云-降水微物理粒子浓度变化较大,但总体变化有一定的范围;(2)采用Γ函数拟合云滴谱更接近实际谱,但不同拟合谱参数差异较大;(3)可用指数函数和Γ函数来拟合层状云降水雨滴谱,Γ函数拟合积云和层积混合云降水雨滴谱精度更高;(4)中国冰核浓度较高,冰核浓度随温度的降低近似成指数变化;(5)冰晶谱、雪晶谱、冰雹谱通常采用指数函数来描述;(6)通常使用荣格(Junge)和Γ函数来分段描述气溶胶粒子谱拟合误差更小。由于云-降水过程及其反馈作用描述不准确是数值模式预报结果不确定性的最大因素,中国正在不断地推进云降水的微物理观测研究,以期进一步加深对东亚季风区云-降水微物理特征的认识,从而为模式中微物理参数化方案的改进提供观测依据和科学指导。基于数值预报模式中云微物理过程参数化发展的需要,总结了中国1960年以来云-降水微物理直接观测的研究成果,可为东亚地区云-降水微物理研究及其模式参数化方案的改进提供观测依据。此外,针对云微物理参化发展的需求,结合过去已有的大量观测提出了几点建议,为今后云-降水物理综合性观测方案的设计提供参考。