The water-bearing numerical model and its operational forecasting experiments part ii: The operational forecasting experiments

The water-bearing numerical model is undergone all round examinations during the operational forecasting ex-periments from 1994 to 1996. A lot of difficult problems arising from the model’s water-bearing are successfully re-solved in these experiments through developing and using a series of technical measures. The operational forecasting running of the water-bearing numerical model is realized stably and reliably, and satisfactory forecasts are obtained.     

东北冷涡数值预报系统建立与业务应用

东北冷涡强降水一直是东北地区气象预报的难点和重点。从提高东北冷涡降水数值预报质量着手,选择模式中可能对东北冷涡强对流有影响的不同辐射方案、边界层方案和积云参数化方案等多种物理过程,对东北冷涡降水预报的性能进行批量试验;利用天气学检验方法,选择预报结果较好的RRTM长波辐射方案、Gayno Seaman边界层方案、Grell积云参数化方案和混合相湿过程显式方案等物理方案,组合建立了东北冷涡数值预报系统。经检验表明：系统对东北冷涡降水预报的水平较东北区域中尺度数值预报业务系统有所提高。同时利用天气学检验方法,从降水预报的强度、落区和移速等方面着手,评估东北冷涡数值预报产品对冷涡降水预报的性能和系统误差,得出降水强度易偏强、落区易偏西偏南和范围易偏小等定性结论,可直接应用于预报业务,取得了较好的应用效果。     

基于风暴数值模拟的冰雹临近预报方法研究

该文利用基于格点冰雹识别的研究成果,结合ARPS风暴数值模式的输出结果,提出基于风暴数值模式的冰雹临近预报方法,即用风暴数值预报的水物质场反演的反射率因子场作为冰雹的预报因子,并通过建立基于格点的强冰雹识别算法作为冰雹预报模型,从而对冰雹的落区及大小做出预报。与一般的冰雹预报模型相比,新的方法有以下特点：选取的冰雹预报因子物理意义更加明确,更加全面;建立的冰雹预报模型比较稳定;建立冰雹预报模型的过程相对简单。新的方法在一次强冰雹过程中得到了成功应用,在3h的临近预报中基本准确预报了强冰雹的落区位置。     

CAMS复杂云微物理方案与GRAPES模式耦合的数值试验

CAMS复杂云微物理方案是混合相双参数方案, 包括11个云物理变量和31个云物理过程, 能够同时预报水成物的比质量和数浓度。通过在GRAPES非静力中尺度模式中增加预报量并修改相关程序后, 实现了二者的耦合, 耦合后模式运行稳定。选取2005年8月15—17日我国华北地区一次暴雨过程, 利用耦合后的模式进行48 h模拟试验, 同时还选取了GRAPES模式中其他3个比较复杂的微物理方案进行模拟, 着重分析了降水和水成物分布的模拟结果。研究结果表明: CAMS方案能够模拟出与实测相接近的雨带分布特征, 并且对降水演变的模拟结果与其他方案比较一致, 对暴雨中心位置的模拟有待改进。CAMS方案模拟的水成物垂直分布与其他方案相比具有相似的总体特征, 各相态粒子的量级和分布合理, 不同方案的结果在量值上有所差别。个例分析结果显示出CAMS方案对降水和水成物的分布能够合理描述。今后应通过更多个例进行更为精细的模拟试验, 对新方案进行检验。     

The impact of microphysical schemes on hurricane intensity and track

During the past decade, both research and operational numerical weather prediction models [e.g. the Weather Research and Forecasting Model (WRF)] have started using more complex microphysical schemes originally developed for high-resolution cloud resolving models (CRMs) with 1–2 km or less horizontal resolutions. WRF is a next-generation meso-scale forecast model and assimilation system. It incorporates a modern software framework, advanced dynamics, numerics and data assimilation techniques, a multiple moveable nesting capability, and improved physical packages. WRF can be used for a wide range of applications, from idealized research to operational forecasting, with an emphasis on horizontal grid sizes in the range of 1–10 km. The current WRF includes several different microphysics options. At NASA Goddard, four different cloud microphysics options have been implemented into WRF. The performance of these schemes is compared to those of the other microphysics schemes available in WRF for an Atlantic hurricane case (Katrina). In addition, a brief review of previous modeling studies on the impact of microphysics schemes and processes on the intensity and track of hurricanes is presented and compared against the current Katrina study. In general, all of the studies show that microphysics schemes do not have a major impact on track forecasts but do have more of an effect on the simulated intensity. Also, nearly all of the previous studies found that simulated hurricanes had the strongest deepening or intensification when using only warm rain physics. This is because all of the simulated precipitating hydrometeors are large raindrops that quickly fall out near the eye-wall region, which would hydrostatically produce the lowest pressure. In addition, these studies suggested that intensities become unrealistically strong when evaporative cooling from cloud droplets and melting from ice particles are removed as this results in much weaker downdrafts in the simulated storms. However, there are many differences between the different modeling studies, which are identified and discussed.     

中尺度碘化银催化数值模式在人工影响天气业务中的应用试验

数值模式越来越多地应用于人工影响天气业务工作中。文中着眼北京地区人工影响天气业务对数值模式预报的需求,结合地基作业的特点,初步进行了碘化银（AgI）冷云催化数值模式在人工影响天气业务中的应用试验。试验采用同化了多种观测资料的北京睿图快速更新和多尺度分析预报系统-临近子系统（RMAPS-ST,Rapid-refresh Multi-scale Analysis and Prediction System-Short Term）的预报结果做为初始场,运用耦合了冷云催化模块的中尺度数值模式WRF（Weather Research and Forecasting Model）进行冷云催化数值试验,催化模块添加在Morrison双参数微物理方案中,考虑了水汽在人工冰核上的凝华核化、凝结冻结和接触冻结核化过程,模拟范围覆盖北京地区所有的地面作业站点以及华北人工影响天气飞机作业区域。业务试验包括运行对照试验以判定具有催化潜力的地基作业站点及催化信息、进行预催化模拟和实际作业后的催化模拟。个例应用显示,冷云催化数值模式的业务试验可以预判未来一定时段内人工影响天气相关物理量场的变化及作业效果,可为前期人工影响天气地面作业方案的制定提供参考,后期可进行实际作业信息输入进行催化模拟,对比验证与分析可以改进和优化流程算法,不断提升制定地面作业方案的科技水平。      

HLAFS显式云降水方案及其对暴雨和云的模拟(I)云降水显式方案

将三套显式云降水方案移植、耦合到国家气象中心业务有限区数值预报模式HLAFS中, 显式云降水方案包括双参数暖云方案、简化混合相云方案和双参数混合相云方案,明显地改进了原模式预报系统的湿物理过程,特别是对云降水过程的模拟能力.详细介绍了这三套云降水方案的物理过程,以及与动力框架的耦合.     

基于扰动模式的四维变分资料同化系统框架的设计完善和数值试验

为了建立一个应用于区域数值预报的四维变分资料同化（4DVar）系统,在近期开发的扰动预报模式GRAPES_PF基础上,开发完善增量四维变分同化系统框架。该框架中暂不包含物理过程（长短波辐射、边界层过程、对流参数化和云微物理等）。对比业务使用的GRAPES 3DVar系统,增加了温度控制变量。将无量纲Exner气压与流函数的线性风压平衡方程直接在地形追随垂直坐标面上求解,且通过广义共轭余差法（GCR）求解扰动亥姆霍兹（Helmholtz）伴随方程。利用人造“探空”资料对2015年10月台风“彩虹”进行了理想数值试验。试验结果表明,所开发的扰动四维变分同化框架得到了预期的结果,即同化更多资料并反复受到模式约束的四维变分同化系统能有效改善初值质量,进而改善区域数值预报。建立的区域四维变分同化框架合理可行,为进一步发展包含完整物理过程的区域四维变分同化系统奠定了研究基础。      

WRF模式不同微物理方案水凝物的预报能力检验与集成试验

文章利用TRMM(Tropical Rainfall Measuring)卫星TMI(TRMM Microwave Imager)探测反演的云水、雨水和冰水新版改进资料,定量化检验了WRF(Weather Research and Forecasting)模式6种微物理方案(Lin,WSM6,Thompson,Morrison 2-mom,CAM5.1,NSSL 2-mom)对江苏近海及周边海域上空云中水凝物的预报能力。19次个例统计检验结果表明,6种微物理方案基本都能预报出水凝物的量级及大致范围,对于云水的预报,除NSSL2-mom方案误差较大外,其余5种方案误差均较小;CAM5.1方案对雨水含量较大时的预报效果较好,但对于雨水含量较小时的预报误差较大;Lin方案对于冰水的预报效果较好,更接近于实况。在此基础上,选用等权集成(EMN)和消除偏差集成(BREM)两种方法开展集成预报试验,结果发现,两种方法均可以减小预报的误差,且消除偏差集成的方法比等权集成的改进效果更好。     

区域气候自忆模式中记忆系数的某些特征及汛期降水预报分析

根据制作中国汛期降水业务预报的需要, 调整了区域气候自忆预测模式的边界条件及计算格式并进行了数值试验; 讨论了模式中记忆系数的线性趋势及周期变化等特征; 对1988~1997年中国汛期降水分布预报试验和1998年实际预报结果进行了分析, 结果表明, 自忆模式对汛期降水年际预报具有一定的能力.     

东北区域空气质量数值预报模式应用研究进展

中国气象局沈阳大气环境研究在2014年基于CUACE和CMAQ建立了东北区域空气质量数值预报业务系统,本文介绍了空气质量数值模式的研究进展以及业务现状,为东北区域空气质量和雾霾预报提供了技术支撑。然而,随着预报精细化和更长预报时效的业务发展需求,存在预报准确率不高、计算资源短缺、科技创新能力不足等问题。本文基于存在的问题提出了东北区域空气质量数值预报未来发展建议与展望,包括加强大气污染源清单研究与技术规范制定、观测资料同化技术研究与业务应用、物理过程参数化方案的改进优化、发展数值预报产品订正技术、开发高分辨率7—10 d数值预报产品、加强人才引进和科技创新等。     

A bin-microphysics cloud model with high-order, positive-definite advection

An axisymmetric, anelastic model of a convective cloud is described. The model comprises prognostic equations for the azimuthal vorticity, the perturbation potential temperature, the perturbation water vapor mixing ratio, 44 categories of cloud condensation nuclei, and 100 categories of liquid-phase hydrometeors. Results from a control simulation show that the model is capable to reproduce realistically the life cycle of a convective cloud including the production of warm rain.A discussion of the role of advection in bin-microphysics models is presented and sensitivity tests were performed regarding the order of advection. The results show that, although the global characteristics of all simulated clouds were similar, significant differences occur with respect to their microstructure, particularly close to the cloud edges, when the order of the advective scheme changes. The conclusion is that intermediate-order advection schemes can indeed be used in cloud-resolving simulations, as far only as the gross characteristics of the cloud/cloud system are being investigated, but not poor, low-order schemes. On the other hand, the sensitivity with respect to the advection suggests that the evaluation of cloud phenomena that occur in fine-scales, such as entrainment and certain microphysical and radiational processes, must require the use of accurate, higher-order schemes.     

载水数值模式的实时预报试验

该文对载水数值模式进行了持续56天的实时预报试验，并对模式的常规预报和云场预报结果进行检验和分析。试验表明，模式运行稳定可靠；载水数值模式有能力合理而比较正确地预报大气中各种水物质的分布及随时间变化，预报云系、云团的宏、微观结构和云物理特征；与此同时，其形势预报及以降水为代表的要素预报能力，亦不低于一般的常规预报模式。     

Impact of horizontal diffusion, radiation and cloudiness parameterization schemes on fog forecasting in valleys

Fog and low stratus forecasting experiments have been carried out with the numerical weather prediction model ALADIN on a case of long lasting fog. The model has been used with different radiation, cloud diagnosing and horizontal diffusion schemes, different representation of orography, increased vertical resolution and with or without prognostic condensates and turbulent kinetic energy (TKE). Some of the numerical set-ups are able to reproduce the fog (low stratus) field as seen in the satellite images as well as the measured 2m temperature and relative humidity diurnal cycles. The results show that cloud diagnosing schemes and overlap assumptions play a more important role than a more sophisticated radiation scheme, or introduction of prognostic cloud water, ice, rain, snow or TKE. More realistic orography representation and a more physical horizontal diffusion scheme significantly improve the modelled low stratus and 2m temperature in the areas with variable orography.     

应用于双向套网格模式的一种变格距差分解法

本文提出了一种用于双向套网格模式的变格距差分计算方案。该方案在不同格距的网格区采用不同精度的差分格式，它自然地连接粗细网格，避免了一般套网格方案在粗细网格相重合点上进行的重复计算。用解析法和数值试验证明了:它与其它一些变格距差分格式相比，对短波的穿透能力有明显改进，虚假的反射也较小。应用该方案建立了正压原始方程双向套网格模式，并采用空间分解和时间分解计算方法。这不仅使二维问题转化为二个一维问题，而且二维套网格也可简化为一维均匀网格和一维套网格两部分，从而使计算和程序简化。用理想场为初值所做的一系列数值试验表明，该模式中的波可以自由进出粗细网格区，计算稳定。最后，还用该模式做了台风路径预报试验，给出了一些试验结果。     

LONG VALID TIME ENERGY PERFECT CONSERVATIVE FIDELITY SPECTRAL SCHEMES OF BAROTROPIC PRIMITIVE EQUATIONS*

In accordance with a new compensation principle of discrete computations,the traditional meteorological global (pseudo-) spectral schemes of barotropic primitive equation (s) are transformed into perfect energy conservative fidelity schemes,thus resolving the problems of both nonlinear computational instability and incomplete energy conservation,and raising the computational efficiency of the traditional schemes.As the numerical tests of the new schemes demonstrate,in solving the problem of energy conservation in operational computations,the new schemes can eliminate the (nonlinear) computational instability and,to some extent even the (nonlinear) computational diverging as found in the traditional schemes,Further contrasts between new and traditional schemes also indicate that,in discrete operational computations,the new scheme in the case of nondivergence is capable of prolonging the valid in-tegral time of the corresponding traditional scheme,and eliminating certain kind of systematical computational "climate drift",meanwhile increasing its computational accuracy and reducing its amount of computation.The working principle of this paper is also applicable to the problem concerning baroclinic primitive equations.     

云分析系统在台风莫拉菲数值模拟中的应用

借鉴美国CAPS开发的中尺度ARPS模式的资料分析系统ADAS,设计了基于GRAPES_Meso的云分析系统,实现了地面云观测资料、卫星云图、多普勒雷达反射率资料在GRAPRS模式中的综合融合应用。通过对登陆台风莫拉菲 (0906) 数值模拟检验云分析系统的性能。结果表明:云分析系统通过地面云观测资料、卫星红外云图、可见光云图、多普勒雷达反射率的同化,能够反演出合理的三维云覆盖状况;在三维云覆盖的基础上结合云底云顶高度,进而反演出云水、云冰、雨水、雪、霰等云微物理量,并显著改善模式初始湿度场;在台风登陆过程的模拟中,虽然对台风路径的预报有所偏差,但对比控制试验,在台风登陆地点、台风强度、路径的预报中云分析发挥了正作用;通过对云分析后初始场反演的雷达回波进行检验,验证了云分析的可行性与正确性,对整体的降水预报效果也有明显改善。     

The ECMWF medium-range prediction models development of the numerical formulations and the impact of increased resolution

Summary An outline is given of the numerical formulations that have been utilized for medium-range prediction at the European Centre for Medium-Range Weather Forecasts. The sigma-coordinate grid-point model adopted for the first phase of operational forecasting, and the hybrid-coordinate spectral model which replaced it, are described. Particular emphasis is placed on aspects of the time-stepping employed in the spectral model. It is shown how a semi-implicit treatment of the zonal advection of vorticity and specific humidity, and selective use of enhanced horizontal diffusion, enables time-steps to be used which are well over 50% longer than would otherwise be possible in a spectral model with a classical semi-implicit treatment of gravity-wave terms, with negligible computational cost and negligible impact on forecast accuracy, at T 106 resolution at least. Some discussion of the stability of the semi-implicit gravity-wave scheme is also included.The experimental evidence which led to the operational change from the grid-point to the spectral model is summarized, and the sensitivity of forecasts to the horizontal resolution of the spectral model is illustrated. Both the change in method and increases in resolution result in systematic improvements in the details of cyclone development, and they occasionally have major beneficial impact on the evolution of the forecast in the medium range. Conventional objective verification indicates a much smaller mean improvement of T 106 resolution over T 63 than of T 63 over T 42, although the advantage of T 106 is nevertheless clear-cut. This advantage is more than confirmed by synoptic assessment. Scope for further improvement is indicated by first results from global forecasts with a resolution of T 159. Indications of sensitivity to the magnitude of horizontal diffusion and the resolution of the model orography are also given.The spectral model was initially used operationally with a 16-level vertical resolution, and the uppermost level at which prognostic variables were held was 25 mb. The results which justified a subsequent change to 19 levels, with increased stratospheric resolution and a 10 mb top level, are presented. In this case the principal medium-range forecast improvement came about because better initial analyses could be produced when the 19-level model was used in data assimilation.With 24 Figures     

TEST OF A TROPICAL LIMITED AREA NUMERICAL PREDICTION MODEL INCLUDING EFFECT OF REAL TOPOGRAPHY

In this paper, a tropical limited numerical prediction p.e. model in sigma coordinate is developed. The predicted variables are deviations from a rest reference atmosphere. This transformation is of benefit in reducing truncation error and guarantees computational stability after the very steep real topography is introduced into the model. The numerical tests with this model show that the large-scale obstacle of the Tibetan Plateau is not negligible in forecasting tropical systems like the South China Sea typhoon.