云的移动轨迹成像模拟及云量变化分析

以全天空数字成像仪的等角投影成像原理为基础，将云型简化为正方体及圆柱体云体。模拟了相同云体在不同空间位置的移动轨迹情况，对其所占面积变化（云量）进行了计算，并对云在移动过程中云体侧面成像情况做了分析研究。结果表明，云量随空间位置变化情况与云的宽高比相关，当宽高比大于某一值时云量随天项角（云所处位置）的增大先增大而后减小，反之则随着天顶角的增大而减小。     

Radiation transfer in finite cylindrical clouds

Summary An analytical method is developed for determining radiation transport in finite cylindrical clouds. The cylindrical form is taken as the idealized shape of cumulus clouds. In the shortwave part of the spectrum illumination by the direct solar beam is considered. In the infrared the cloud's emission as well as the radiation emitted by an underlying ground is taken into account. Numerical computations are carried out to illustrate the applicability of the technique developed in the present study. For a cloud in an absorbing midlatitude summer atmosphere results are shown for the shortwave region including spectral integration. Cloud transmission and albedo are presented as a function of solar zenith angle and as a function of the cloud's dimensions. The results are in agreement with published results from Monte Carlo calculations as far as a comparison is possible. Vertical profiles of the solar heating rate in the cylindrical cloud indicate the effect of the cloud's finite dimensions. Similarity as well as differences are found from comparison with results obtained for a cuboidal cloud from a finite analytical model that was tested successfully against Monte Carlo calculations.With 5 Figures     

新疆乌鲁木齐地区冬季层积云研究——数值模拟研究

在分析实测结果的基础上,建立了一个边界层与云物理学的数值模式,模拟了云生成的过程。模式包括了边界层内的湍流混合、大气与地表的长波辐射和太阳短波辐射以及云层的辐射平衡等,同时也包括了水的汽、液、固三相物质在大气中的相互转换即云微物理过程。模拟的云及边界层的主要结构特征和演变趋势与实测情况大致吻合。模拟结果还提供了有关云结构及演变过程的一些详细图象。     

Cloud-radiation feedbacks in a general circulation model and their dependence on cloud modelling assumptions

The general circulation model (GCM) used in this study includes a prognostic cloud scheme and a rather detailed radiation scheme. In a preceding paper, we showed that this model was more sensitive to a global perturbation of the sea surface temperatures than most other models with similar physical parametrization. The experiments presented here show how this feature might depend on some of the cloud modelling assumptions. We have changed the temperature at which the water clouds are allowed to become ice clouds and analyzed separately the feedbacks associated with the variations of cloud cover and cloud radiative properties. We show that the feedback effect associated with cloud radiative properties is positive in one case and negative in the other. This can be explained by the elementary cloud radiative forcing and has implications concerning the use of the GCMs for climate sensitivity studies.     

三维强风暴动力-电耦合数值模拟研究--模式的初步检验

利用文献「１」所发展的三维强风暴动力－电耦合模式对ＣＣＯＰＥ（Ｃｏｏｐｅｒａｔｉｖｅ Ｃｏｎｖ－ｅｃｔｉｖｅＰｒｅｃｉｐｉｔａｔｉｏｎＥｘｐｅｒｉｍｅｎｔ,１９８１）计划中７月１９日的一次强雷暴过程进行了模拟检验,并对雷暴云宏观和微观两个方面实测和模拟值进行了对比。通过分析表明,该模式基本能反映出雷暴云在各发展阶段上动力、微物理和电学特征的时空演变规律,模拟结果与实际观测结果较为接近。     

MIXED-PHASE STRATIFORM CLOUD SYSTEM MODEL AND CASE MODELING ON TWO LOW-LEVEL MESOSCALE VORTICES

We introduced the two-parameter stratiform cloud model of Hu and Yan (1986) into the mesoscale model of Anthes et al.(1987),and reprogramed the latter,then constructed a three-dimensional stratiform cloud system model which includes three phases of water and detailed cloud physical processes.For the stability and accuracy of calculation in a larger time step,we accepted a set of hybrid-schemes for all and the time split scheme for some of the cloud physical processes,and proposed a parameterized method which calculates different types of phase change processes simultaneously,and designed the falling schemes of particles following the Lagrangian method.We used a dry model,a cumulus parameterization model,a two-phase explicit scheme model,and the model presented here to simulate two low-level mesoscale vortices,compared and analysed the simulating capability of these models.The results show that in simulation of the circulation structure of meso-vortex,the structure of cloud system,and surface precipitation,the model presented here is more reasonable and closer to the observations than other models.     

显式云物理方案的研究进展

回顾了近年来显式云物理方案的研究进展。显式云方案主要有体积水方法和详细微物理方法(分档法)。体积水方法有单参数和双参数两种谱描述方法,根据模式预报变量和物理过程的不同,可以分为暖云方案、简单冰相方案和复杂冰相方案。详细的微物理方法由于预报变量繁多、计算量巨大而一般多应用于研究工作。不同的模式,有不同的显式云方案,并不是粒子分类越复杂模拟效果就越好,需要根据研究的重点、计算资源的许可选择使用不同的物理方案。物理过程参数化需要建立在理论和实验研究的基础上,因此应加强这方面的理论和实验研究,使物理参数化具有更坚实的物理基础。     

云微物理参数化对华北降雪影响的数值模拟

对发生在华北地区的一次降雪过程进行了中尺度数值模拟。结果表明,高纬强冷空气南下和低纬倒槽的水汽输送是造成这次长时间降雪过程的主要原因。采用混合方案的中尺度数值模拟表明,这次降雪天气不是对流云造成的,而是稳定性的非对流云降雪。敏感性试验也表明,采用其他积云参数化方案对模拟的降雪量基本没有影响。控制试验模拟的24h降雪量与实际观测比较吻合。模拟结果表明,当采用Dudhia简单冰相方案时,会有过多的云冰、过冷却水及雪;当采用Reisner 1混合相方案时,会有过多的云冰和雪;修改的各个Reisner 2方案对此次降雪的预报改进不大,但各个Reisner 2方案的敏感性试验中云冰混合比、过冷却水混合比和雪混合比稍微有差异。     

Diagnosis and testing of low-level cloud parameterizations for the NCEP/GFS model using satellite and ground-based measurements

The objective of this study is to investigate the quality of clouds simulated by the National Centers for Environmental Prediction global forecast system (GFS) model and to examine the causes for some systematic errors seen in the simulations through use of satellite and ground-based measurements. In general, clouds simulated by the GFS model had similar spatial patterns and seasonal trends as those retrieved from passive and active satellite sensors, but large systematic biases exist for certain cloud regimes especially underestimation of low-level marine stratocumulus clouds in the eastern Pacific and Atlantic oceans. This led to the overestimation (underestimation) of outgoing longwave (shortwave) fluxes at the top-of-atmosphere. While temperature profiles from the GFS model were comparable to those obtained from different observational sources, the GFS model overestimated the relative humidity field in the upper and lower troposphere. The cloud condensed water mixing ratio, which is a key input variable in the current GFS cloud scheme, was largely underestimated due presumably to excessive removal of cloud condensate water through strong turbulent diffusion and/or an improper boundary layer scheme. To circumvent the problem associated with modeled cloud mixing ratios, we tested an alternative cloud parameterization scheme that requires inputs of atmospheric dynamic and thermodynamic variables. Much closer agreements were reached in cloud amounts, especially for marine stratocumulus clouds. We also evaluate the impact of cloud overlap on cloud fraction by applying a linear combination of maximum and random overlap assumptions with a de-correlation length determined from satellite products. Significantly better improvements were found for high-level clouds than for low-level clouds, due to differences in the dominant cloud geometry between these two distinct cloud types.     

IAP AGCM中短波辐射方案的改进研究I·引入Fu-Liou短波辐射方案

冰云和水云对短波辐射性质(消光系数、单次散射反照率及不对称因子)的影响很不相同,应分别计算。Fu-Liou短波辐射方案(以下称Fu-Liou code)就是对冰云和水云分别采用了不同的参数化方案,云的短波辐射性质直接由云的物理性质来确定。因此,Fu-Liou code在云的处理方面物理意义更清晰且很合理。作者将Fu-Liou code引入IAP AGCM-II中,称为Version 2。对当代气候场的模拟结果表明,Version 2的各个物理过程是协调匹配的,且其对气候场的模拟性能是好的,从而为进一步改进IAP AGCM的短波辐射方案提供了很好的模式基础。     

一个简化的混合相云降水显式方案

该文提出一个新的混合相云降水显式方案,它预报2个云物理量,即云水（冷区为过冷云水）和降水（冷区为冰雪,暖区为雨）,考虑了7种云物理过程。文中给出了详细的方程组,可以作为一个子程序供大、中尺度天气模式使用。该方案还与详细微物理显式方案和暖云方案作了实例模拟比较。     

A PDF-based hybrid prognostic cloud scheme for general circulation models

A new cloud parameterization based on prognostic equations for the subgrid-scale fluctuations in temperature and total water content is introduced for global climate models. The proposed scheme, called hybrid prognostic cloud (HPC) parameterization, employs simple probability density functions (PDFs) to the horizontal subgrid-scale inhomogeneity, allowing them to vary in shape in response to small-scale processes such as cumulus detrainment and turbulent mixing. Simple tests indicate that the HPC scheme is highly favorable as compared to a diagnostic scheme in terms of the cloud fraction and cloud water content under either uniform or non-uniform forcing. The relevance of the HPC scheme is investigated by implementing it in an atmospheric component model of the climate model MIROC with a coarse resolution of T42. A comparison of the short-term integrations between the T42 model and a global cloud resolving model (GCRM) reveals that the HPC scheme can reproduce, to a certain degree, the subgrid-scale variance and skewness of temperature and total water content simulated in the GCRM. It is also found that the HPC scheme significantly alters the climatological distributions in cloud cover, precipitation, and moisture, which are all improved from the model using a conventional diagnostic cloud scheme.     

A first order mixing–condensation scheme for nocturnal stratocumulus

A simple closure scheme for nocturnal stratocumulus is proposed. The scheme is formulated in conserved variables. Cloud fraction and cloud water amount are diagnosed assuming a top-hat distribution for total water. Conversion of cloud water into rain water is parameterized in terms of cloud water and the incoming rain flux. Turbulence transport in the cloud layer is accounted for by a first-order vertical diffusion scheme with a profile-type diffusivity. The length scale corresponds to the thickness of the cloud layer. The turbulent velocity scale is directly related to the long wave radiative flux divergence in the cloud. Entrainment at cloud top is implicitly treated by extending the in-cloud mixing profile slightly beyond cloud top. The excess height is derived from the buoyancy frequency at cloud top and a radiative–convective velocity scale. The scheme is capable of simulating realistic profiles of the conserved variables and cloud parameters for a case of nocturnal stratocumulus prepared on the basis of ASTEX data.     

用多普勒雷达反射率调整模式大气的云微物理变量

一种简单云分析方案, 用于由多普勒天气雷达反射率反演中尺度大气模式初值分析中的云微物理变量（云水混合比和雨水混合比）和空气湿度变量（比湿），使模式积分初始场反映出观测空间的云微物理特征以及哪些空间位置上的大气处于饱和状态。应用于2002年6月梅雨期安徽省马鞍山市一次降水过程的临近数值预报试验结果表明，模式预报的大气综合反射率与雷达观测的回波图像相近，由云微物理变量变化表示的模式云系演变与雷达观测的回波图像一致, 伴随模拟的中小尺度云系, 模式大气能很快调整出合理的中小尺度流场辐散、辐合结构；它们明显好于模式初始场不引入雷达反射率时的结果，即这种方法对改进临近数值天气预报准确率是有效的。     

Variational assimilation of satellite cloud water/ice path and microphysics scheme sensitivity to the assimilation of a rainfall case

Hydrometeor variables (cloud water and cloud ice mixing ratios) are added into the WRF three-dimensional variational assimilation system as additional control variables to directly analyze hydrometeors by assimilating cloud observations. In addition, the background error covariance matrix of hydrometeors is modeled through a control variable transform, and its characteristics discussed in detail. A suite of experiments using four microphysics schemes (LIN, SBU-YLIN, WDM6 and WSM6) are performed with and without assimilating satellite cloud liquid/ice water path. We find analysis of hydrometeors with cloud assimilation to be significantly improved, and the increment and distribution of hydrometeors are consistent with the characteristics of background error covariance. Diagnostic results suggest that the forecast with cloud assimilation represents a significant improvement, especially the ability to forecast precipitation in the first seven hours. It is also found that the largest improvement occurs in the experiment using the WDM6 scheme, since the assimilated cloud information can sustain for longer in this scheme. The least improvement, meanwhile, appears in the experiment using the SBU-YLIN scheme.     

Design and performance of a new cloud microphysics scheme developed for the ECHAM general circulation model

A new cloud microphysics scheme including a prognostic treatment of cloud ice (PCI) is developed to yield a more physically based representation of the components of the atmospheric moisture budget in the general circulation model ECHAM. The new approach considers cloud water and cloud ice as separate prognostic variables. The precipitation formation scheme for warm clouds distinguishes between maritime and continental clouds by considering the cloud droplet number concentration, in addition to the liquid water content. Based on several observational data sets, the cloud droplet number concentration is derived from the sulfate aerosol mass concentration as given from the sulfur cycle simulated by ECHAM. Results obtained with the new scheme are compared to satellite observations and in situ measurements of cloud physical and radiative properties. In general, the standard model ECHAM4 and also PCI capture the overall features, and the simulated results usually lie within the range of observed uncertainty. As compared to ECHAM4, only slight improvements are achieved with the new scheme. For example, the overestimated liquid water path and total cloud cover over convectively active regions are reduced in PCI. On the other hand, some shortcomings of the standard model such as underestimated shortwave cloud forcing over the extratropical oceans of the respective summer hemisphere are more pronounced in PCI.This paper was presented at the Third International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 4–9 Sept. 1995 under the auspices of the Max Planck Institute for Meteorology, Hamburg. Editor for these papers is L. Dümenil.     

The implementation of reverse Kessler warm rain scheme for radar reflectivity assimilation using a nudging approach in New Zealand

Reverse Kessler warm rain processes were implemented within the Weather Research and Forecasting Model (WRF) and coupled with a Newtonian relaxation, or nudging technique designed to improve quantitative precipitation forecasting (QPF) in New Zealand by making use of observed radar reflectivity and modest computing facilities. One of the reasons for developing such a scheme, rather than using 4D-Var for example, is that radar VAR scheme in general, and 4D-Var in particular, requires computational resources beyond the capability of most university groups and indeed some national forecasting centres of small countries like New Zealand. The new scheme adjusts the model water vapor mixing ratio profiles based on observed reflectivity at each time step within an assimilation time window. The whole scheme can be divided into following steps: (i) The radar reflectivity is firstly converted to rain water, and (ii) then the rain water is used to derive cloud water content according to the reverse Kessler scheme; (iii) The cloud water content associated water vapor mixing ratio is then calculated based on the saturation adjustment processes; (iv) Finally the adjusted water vapor is nudged into the model and the model background is updated. 13 rainfall cases which occurred in the summer of 2011/2012 in New Zealand were used to evaluate the new scheme, different forecast scores were calculated and showed that the new scheme was able to improve precipitation forecasts on average up to around 7 hours ahead depending on different verification thresholds.     

相变修正方案在GRAPES模式标量平流中的应用

如何更好地模拟水物质的空间分布和小尺度变化,对于数值天气预报效果的改进,特别是对于更好地模拟降水过程,具有重要的意义.计算机的飞速发展使数值模式的分辨率不断提高,云的显式计算成为可能,这样就要求水物质在平流的过程中必须要做到高精度、守恒、保形.水物质场是正定标量的场,具有空间和时间变化幅度大、存在强梯度甚至不连续的特点,水物质场的合理模拟一直是数值预报中的一个难题.GRAPES模式中的标量平流方案采用PRM分段有理函数方法,比较好地解决了该半拉格朗日模式中水物质平流的高精度、守恒、保形问题,但是当有凝结潜热发生时,由于半拉格朗日平流方案求解上游点时的插值,在云边缘区域会造成虚假的云水,进而导致不合理的相变过程.为了解决以上问题,本研究在GRAPES模式中PRM平流方案的案础上,加入了非线性半拉格朗日相变潜热的修正方案,旨在改进GRAPES模式对水物质平流问题的模拟,提高降水的预报效果.该研究通过理想试验,验证了非线性半拉格朗口相变修正方案可以有效地限制云边缘由于半拉格朗日平流方案插值产生的虚假柑变;然后将该方案加入GRAPES模式的PRM水物质平流方案中,通过实际个例模拟验证了加入非线性半拉格朗日方案以后,模式可以更好地模拟水物质的平流过程,且对云中热力场及水物质分布地模拟更加合理,同时预报出的雨带中心区与实况更加符合.     

WRF模式云参数化方案对一次深对流系统模拟的验证和改进

使用中尺度数值模式WRF中的双参数云微物理方案WDM6针对2008年台风“凤凰”登陆过程中造成的强降水进行数值模拟,通过卫星模拟器利用MTSAT-1R和TRMM卫星观测的红外云顶黑体亮温TBB、PR雷达反射率资料使用统计方法验证模拟结果。通过修改云水向雨水自动转化过程、冰晶核化过程、雪和霰的下落末速度、雪和霰的截距进行敏感性试验,减小模拟结果和卫星观测结果的差异。研究结果表明:WDM6方案模拟的台风“凤凰”登陆后的降水,强对流云系及对流柱状雷达回波基本符合实况,但模拟结果局部偏强。WDM6方案模拟产生了较多的浅对流云,低估了对流云系的出现频率。不同云类型模拟的雷达回波均偏强,对流云系雷达回波垂直分布接近观测。敏感性试验结果说明修改WDM6方案中云水向雨水自动转化率有效地改善了模拟效果。同时发现云滴初始数浓度影响云水向雨水自动转化率并最终影响云系结构和雷达反射率的模拟结果,过高的云滴初始数浓度会使模拟结果变差。     

晴空与有云大气辐射分布的数值模拟及其对全天空图像云识别的应用

全天空图像自动云识别的研究相对来说是一项较新的研究领域。当前国际上较为通用的全天空图像云识别方法主要依靠图像蓝红灰度值对比的阈值方法进行判断。然而非洁净大气中气溶胶的增多给云识别增加了难度,同时太阳高度角不同天空色度分布情况也不同。文中利用Libradtran辐射传输模式,计算了不同能见度不同太阳高度角情况下3个典型波长(450,550,650nm;蓝/绿/红)无云及有云大气的天空辐亮度分布情况,并进行了比较分析。结果表明,相同太阳高度角情况下,无云及有云大气中蓝红比值随能见度的下降呈单调下降趋势。在特定的云光学厚度和能见度情况下,天空色度彼此呈现出类似的分布状况。全天空图像阈值判断云识别自适应算法的建立需要与太阳高度角、地面能见度联系起来。当前尚无法建立一个判断阈值随太阳高度角以及能见度变化的函数关系式。较为可行的办法是建立典型能见度、典型太阳高度角情况下的辐射信息库,在具体云识别时,首先确定太阳高度角,而后根据天空辐射比情况确定天空能见度,并利用辐射信息标准库做云或非云判别。该工作为全天空云识别算法提供判别依据,同时建立云识别随能见度和太阳高度角变化的判别信息库。