南半球环状模事件的准地转调整过程分析

本文利用NCEP/NCAR逐日再分析资料,分析了南半球环状模（SAM）事件生命过程中的准地转调整过程.由于SAM沿纬圈的水平尺度远远大于临界尺度罗斯贝变形半径,因此要求纬向风场在地转调整过程向气压场适应以达到地转平衡.研究结果表明,在纬向平均环流中,异常Ferrel环流强度的变化超前于SAM强度变化约2/16位相,异常Ferrel环流能够通过超前的整层大气质量的经向输运,改变中高纬度的质量分布装状况,导致中、高纬度地区间的位势高度梯度异常变化,而中、高纬度地区间的位势高度梯度异常发生变化就意味着SAM强度和位相发生变化;而当SAM强度和位相发生改变后,即中、高纬度地区南北方向上的位势高度梯度发生变化后,可破坏中纬度地区纬向风场与位势高度场之间的地转平衡,产生地转偏差;地转偏差产生后,又可驱动经向散度风场,造成Ferrel环流异常的变化,由此形成一个自我内部循环调整过程.     

数值模式中气压梯度力的算法试验

本文在差微差一致性二阶精度的坐标变换公式的基础上,推导出了气压梯度力计算的新公式。同时,根据静力扣除法原理,提出了初值静力扣除法,用以计算有地形数值模式中的气压梯度力。将这两种新的方法与二阶计算精度的回插法和局地等温法(Corby格式)相比后发现:新方法(文中称为经典修正法)与回插法精度相近,但可节省计算时间,比局地等温法精度高。静力扣除法则效果较差,但其原因可能在于地转风垂直插值方法有误差。本文的试验模式是上海台风研究所的业务数值预报模式,文中对模式及若干技术处理也作了简要介绍。     

一种初始方程中应用实测风进行风场初值化的方法

本文提出一种利用实测风场、高度场、温度场和地面三小时气压变化资料,考虑风场和质量场之间基本满足准地转平衡关系,大气具有准无辐散性,但又保留一部分与地面气压变化相应散度的地转偏差,将以实测的等压面风场插到σ坐标的风场进行订正的初值化方法,并用一个单时资料进行初步试验.文中还对风场初值化问题进行了讨论.     

A study of nonhydrostatic effects in idealized sea breeze systems

Nonhydrostatic effects in two-dimensional mesoscale sea breeze systems are investigated by numerical simulations. It is shown that nonhydrostatic effects are directly contributed by the vertical gradients of the vertical velocity variance as well as by the vertical accelerations. It is also shown that a K-type turbulence closure is not suitable in a nonhydrostatic primitive equation model, and a higher-order closure scheme should therefore be used. Results from hydrostatic and fully-nonhydrostatic models are compared for various surface and atmospheric background conditions, such as scale and strength of surface heating, geostrophic wind, stability, surface roughness contrast, Coriolis effect, etc. It is found that for strongly developed sea breeze cases, vertical gradients of vertical velocity variance contribute most to nonhydrostatic forcing in the lower layers, and that the resultant nonhydrostatic pressure gradient acts against the hydrostatic pressure gradient, so that nonhydrostatic simulations produce weaker systems than hydrostatic ones. For weak sea breeze systems, the difference between the two models tends to be small.     

Meso-beta scale numerical simulations of terrain drag-induced along-stream circulations Part I: Midtropospheric frontogenesis

Summary The problem of along-stream ageostrophic frontogenesis is studied by employing a numerical model at meso-alpha and meso-beta scales in simulations of the downstream circulations over the Front Range of the Rocky Mountains. Three-dimensional real data simulations at these two scales of motion are used to diagnose the transition from semigeostrophic cross-stream frontogenesis accompanying a propagating baroclinic upper-level jet streak to midtropospheric along-stream ageostrophic frontogenesis. This along-stream ageostrophic frontogenesis results from the perturbation of the jet streak by the Rocky Mountain range. The case study represents an example of internal wave dynamics which are forced by the drag of the Rocky Mountains on a strong jet streak in the presence of a low-level inversion.The simulation results indicate that, unlike semi-geostrophic frontogenesis, a front (which is alligned perpendicular to the axis of the jet stream) may form when significant adiabatic heating occurs within a stratified shear flow over horizontal length scales shorter than the Rossby radius of deformation. The mechanism responsible for the frontogenesis is the growth of the divergent along-stream wind velocity component which becomes coupled to the front's along-stream pressure gradient force. This nonlinear interaction produces hydrostatic mesoscale frontogenesis as follows: 1) vertical wind shear in the along-stream plane strengthens resulting in the increasingly nonuniform vertical variation of horizontal temperature advection as the ageostrophic wind component grows in magnitude downstream of the meso-scale terrain-induced adiabatic heating, 2) increasing along-stream differential vertical motions (i.e., along-stream thermally indirect circulation with warm air sinking to the west and cold air rising to the east) tilt the vertical gradient of isentropes into the horizontal as the vertical temperature gradient increases due to the previous process in proximity to horizontal gradients in the along-stream component of the ageostrophic wind, 3) as tilting motions act to increase the along-stream horizontal temperature gradient, the along-stream confluence acts to nonuniformly increase the along-stream frontal temperature gradient which increases the along-stream pressure gradient force resulting in further accelerations, ageostrophy, and frontal steepening as part of a scale contraction process.The evolution of the aforementioned processes results in the three-dimensional hydrostatic frontogenesis accompanying the overturning of isentropic surfaces. These adjustments act to turn air parcels to the right of the southwesterly geostrophic wind vector at successively lower atmospheric levels as the scale contraction continues. This simulated along-stream front is verified from diagnostic analysis of the profiler-derived temperature and wind fields.With 17 Figures     

The local spline vertical interpolation method of temperature and geopotential height fields and the time-dependent difference form of the hydrostatic equation

In numerical weather prediction (NWP), the accuracy of vertical interpolation of the initial data is a problem which is greatly concerned by people. In this paper, we specify vertical distributions of the temperature and the geopotential height fields and examine three interpolation methods, i.e. the Lagrangian polynomial inter-polation method (hereafter abbreviated to LP method), the linear interpolation method (LN method) and the local spline interpolation method (LS method) proposed by the author. The examination shows that when the vertical resolution of the initial data is high enough, for example, the number of the given data levels N is 10 or more, all the three methods get good accuracy of interpolation, especially, the LP and the LS methods have very little errors almost tending to zero, while the LN method has a little larger errors than the two formers and the errors at various levels have the same sign. When N is reduced to 5, the LP and the LS methods still have quite good accuracy and similar error distributions, while the LN method has less accuracy. If the geopo-tential height field needs to be adjusted in order to satisfy the hydrostatic equilibrium with the temperature field which is assumed fixed, then the LS method has minimum errors. The examination also indicates that the vertical resolution with at least 5 levels of initial data can keep the interpolation accuracy. Otherwise the accuracy will not be guaranteed no matter which method is used.It is also pointed out in this paper that the temperature and the geopotential height fields can be given inde-pendently in numerical prediction models in order to keep higher interpolation accuracy. However, the hydro-static equation should be finite differenced in other way which is somewhat different from the conventional one. In other words, the time dependent difference form of the equation should be used, so that the initial interpola-tion accuracy could have influence on the time integration.     

Accurate pressure gradient calculations in hydrostatic atmospheric models

Most models of atmospheric flow which use the primitive equations require a diagnostic equation to determine local total pressure. In hydrostatic models, this equation is the vertically integrated hydrostatic equation. A frequently used approximation to this integration is to hold the temperature constant within model layers yielding a linear proportionality between p or (Exner's function) and z. This procedure yields static pressures with errors on the order of 10^–3mb.If terrain following coordinates are used, terms arise in the horizontal momentum equations involving the gradient of total pressure along the coordinate surface, less a correction for the variation of the hydrostatic pressure along a sloped surface. Erroneous horizontal accelerations are common in these models which result from spurious pressure gradients that are due to inaccurate computation of the static pressure. This error may be amplified if the computation of the slope correction term of the horizontal pressure gradient is not consistent with the method of calculating the total pressure.We derive a methodology to be used in the vertical pressure integrations that is exact if the potential temperature lapse rate is constant between integration limits. The method is applied to both the integration of the hydrostatic equation and the computation of the slope correction term in the horizontal pressure gradient. The method employs a fixed vertical grid and a dynamic one defined by the significant levels in the vertical temperature distribution. With this methodology, the error in calculation of the horizontal pressure gradient acceleration is greatly reduced, especially in situations where the isothermal surfaces are not parallel to the vertical coordinate surfaces. The problem of aliasing and the treatment of significant temperature levels is described.     

Marginal sea surface temperature variation as a pre-cursor of heat waves over the Korean Peninsula

This study examines the role of the marginal sea surface temperature (SST) on heat waves over Korea. It is found that sea surface warming in the south sea of Korea/Japan (122-138°E, 24- 33°N) causes heat waves after about a week. Due to the frictional force, the positive geopotential height anomalies associated with the south sea warming induce divergent flows over the boundary layer. This divergent flow induces the southerly in Korea, which leads to a positive temperature advection. On the other hand, over the freeatmosphere, the geostrophic wind around high-pressure anomalies flows in a westerly direction over Korea during the south sea warming, which is not effective in temperature advection. Therefore, the positive temperature advection in Korea due to the south sea warming decreases with height. This reduces the vertical potential temperature gradient, which indicates a negative potential vorticity (PV) tendency over Korea. Therefore, the high-pressure anomaly over the south sea of Korea is propagated northward, which results in heat waves due to more incoming solar radiation.     

Pressure gradient force and calculation initial temperature in terrain-following coordinate models

Summary A method for temperature initialization in the terrainfollowing coordinate models from pressure gradient force interpolated from constant pressure surface is considered. A non-unique solution of initial temperature obtained from regularized interpolated pressure gradient force, is discussed. An approximative procedure for defining unique temperature is presented. The proposed method is tested using real data. The obtained results show advantages of the proposed technique compared with the initial temperature calculation based on interpolated geopotential via hydrostatic equation.

---

Zur Initialisierung des Temperaturfeldes mit Hilfe der Druckgradientkraft in einem der Topographie folgenden Koordinatensystem

Zusammenfassung Es wird eine neue Methode zur Initialisierung des Temperaturfeldes mit Hilfe der Druckgradientkraft in einem numerischen Modell, dessen Koordinaten der Topographie folgen, vorgestellt. Die Methode wurde anhand realen Datenmaterials getestet und es zeigen sich Vorteile gegenüber der konventionellen Initialisierung des Temperaturfeldes mit Hilfe der hydrostatischen Grundgleichung.

---

---

With 3 Figures     

地转动量近似下边界层中风场的调整

本文利用地转动量近似,并假设气压场为定常的圆形涡旋和初始风场不满足四力平衡(气压梯度力、科里奥利力、湍流粘性力和半地转惯性力)的条件下,求解了正压边界层中风场向气压场调整的初边值问题,得到了一些初步结论。本工作为利用四力平衡下的风速分布来诊断预报边界层风场提供了理论依据。     

多级扰动法——静力扣除及其在非静力模式设计中的可应用性

在尺度分析的基础上,利用大气垂直运动方程中各项量级的分布特点,提出了多级扰动法.通过它可以使该方程中量级最大的垂直气压梯度力项(VPGF)和重力项(G)的相应扰动项的量级随级次n的增加明显减小,而其他项的量级不变;而且,方程中最大的垂直截断误差的量级也随n的增加而减小,直到高级扰动项不再是扰动方程中仅有的最大项时为止....     

Geostrophic Wind, Gradient Wind, Thermal Wind And The Vertical Wind Profile – A Sample Analysis Within A Planetary Boundary Layer Over Arctic Sea-Ice

For the interpretation of many boundary-layer field experiments the geostrophic wind is needed as an external parameter. However, quite often and especially in remote areas this wind is not known at all or difficult to determine because there are not enough measurements of the surface pressure.Here it is shown how measurements carried out with the HELIPOD system, a helicopter-borne meteorological turbulence measuring system, may be used to evaluate the geostrophic wind. This is done by the analysis of the pressure field at different heights. An additional analysis of the temperature field in the same heights allows for the discussion of the quality of the derived geostrophic wind. An intercomparison with the vertical wind profile enables us to discuss the influence of the curvature of the isobars. From this, finally, also the gradient wind can be estimated.     

GRAPES业务全球集合控制预报初值的动力升尺度方法

以生成GRAPES全球集合预报业务系统的控制预报初值为目的,基于GRAPES全球模式,开展了控制预报初值生成方法研究,发展了高分辨率初值动力升尺度方法,并检验了不同方法的可行性。通过对比不同方法产生的初始场形态,证实了仅仅对高分辨率初始场进行二维水平插值存在不足,基于静力学方程对Exner气压变量进行三维插值至关重要。结果表明,动力升尺度方法利用静力平衡关系,对全场的温压场进行调整,使之协调平衡,可以改善二维水平插值方法导致的初始位势高度场和温度场的噪音问题,产生适用于GRAPES全球集合预报业务系统的控制预报初值。     

利用MM5模式及NCEP资料对沈阳一次污染天气的数值模拟

利用NCEP/NCAR再分析资料和中尺度天气模式MM5对2010年1月14—19日沈阳大气污染天气系统进行模拟分析。对此次天气过程的地面和高空气压场、地面至高空各高度层随时间变化的水平风场及垂直剖面风场、垂直方向温度廓线等气象要素进行分析和模拟,描述大气污染中天气系统的变化过程,分析造成大气污染的气象要素变化。结果表明：此次污染天气过程对应地面场为长白山高压、地形槽环流型;500 hPa高空天气形势为两槽一脊,地面风场主要受高压辐散气流控制;地面至高空不同高度的水平风场均有偏南风切变和偏西风切变,垂直剖面风场对应有下沉气流,地面至高空的温度廓线出现明显的逆温。这些气象条件共同造成了持续污染天气。而500 hPa位势高度场持续长时间两槽一脊的环流形势,是造成长时间污染天气的主要原因。     

A comparison of two methods for the assimilation of altimeter data into a shallow-water model

Two simple methods for assimilating ocean surface pressure data into a three-layer adiabatic primitive equation model are tested and compared using a twin-experiment approach. One of the methods is based on a recently presented direct insertion scheme using a quasi-geostrophic model. The method explicitly avoids making changes to the potential vorticity fields in the lower layers and the vertical current structure after assimilation is determined by this criterion. The modification required for a primitive equation framework are discussed and a comparison is made with the more familiar ‘nudging’ method of assimilation. We use root mean square errors to quantify the response of the model to the assimilation of a single complete surface pressure field and also to repeated or intermittent data available every 20 days for a period of 1 year. The two methods are almost equally effective in the short term, although the direct insertion scheme appears to be more effective in the longer term. It is suggested that this is because insertion maintains a geostrophic balance at all levels and thus avoids the generation of internal gravity waves which are needed when nudging is used to restore geostrophic currents in the deeper layers.     

Dynamics of the planetary boundary layer with fine structure in the large scale background field

Summary On the basis of Wu and Blumen's work (1982) on the geostrophic momentum approximation (GMA) in the planetary boundary layer (PBL) and Tan and Wu (1992, 1994) on the Ekman momentum approximation (EMA) in the PBL, some improvements about the eddy exchange coefficientK, the advective inertial force and the lower boundary condition of the PBL are developed in this paper: (1) apply theK which is a gradually varying function of height instead of a constant value in the Ekamn layer, and introduce a surface layer; (2) take the effect of the vertical advective inertial force into account; (3) the solution technique is extended from level terrain to orographically formed terrain. Under the condition of the equilibrium among four forces (the pressure — gradient force, Coriolis force, eddy viscous force and inertial force including horizontal and vertical advective inertial forces), we have obtained the analytical solutions of the distributions of the wind and the vertical velocity. The computation of an individual example shows that: (1) both the wind velocity near surface and the angle between which and the non-viscous wind are more consistent with usual observations than that of Wu and Blumen (1982); (2) comparing with the horizontal advective inertial force, the vertical advective inertial force can not be neglected, when the orography is considered, the effect of the latter is even more important than the former.With 3 Figures     

A numerical study of geostrophic adjustment and frontogenesis

ＡＮｕｍｅｒｉｃａｌＳｔｕｄｙｏｆＧｅｏｓｔｒｏｐｈｉｃＡｄｊｕｓｔｍｅｎｔａｎｄＦｒｏｎｔｏｇｅｎｅｓｉｓＭａｊｉｄＭ．ＦａｒａｈａｎｉａｎｄＷｕＲｏｎｇｓｈｅｎｇ（伍荣生）ＤｅｐａｒｔｍｅｎｔｏｆＡｔｍｏｓｐｈｅｒｉｃＳｃｉｅｎｃｅｓ,Ｎａｎｊｉ...     

An Online Model Correction Method Based on an Inverse Problem:PartⅡ——Systematic Model Error Correction

An online systematic error correction is presented and examined as a technique to improve the accuracy of real-time numerical weather prediction, based on the dataset of model errors(MEs) in past intervals. Given the analyses, the ME in each interval(6 h) between two analyses can be iteratively obtained by introducing an unknown tendency term into the prediction equation, shown in Part I of this two-paper series. In this part, after analyzing the 5-year(2001–2005) GRAPESGFS(Global Forecast System of the Global and Regional Assimilation and Prediction System) error patterns and evolution,a systematic model error correction is given based on the least-squares approach by firstly using the past MEs. To test the correction, we applied the approach in GRAPES-GFS for July 2009 and January 2010. The datasets associated with the initial condition and SST used in this study were based on NCEP(National Centers for Environmental Prediction) FNL(final) data.The results indicated that the Northern Hemispheric systematically underestimated equator-to-pole geopotential gradient and westerly wind of GRAPES-GFS were largely enhanced, and the biases of temperature and wind in the tropics were strongly reduced. Therefore, the correction results in a more skillful forecast with lower mean bias and root-mean-square error and higher anomaly correlation coefficient.     

Spectral nudging sensitivity experiments in a regional climate model

In this study, the scale selective bias correction (SSBC) method described by Kanamitsu et al. (2010) is further examined by considering the full wind nudging and the vertically weighted damping coefficient. The 2001 June?CJuly?CAugust RSM simulation over a relatively large domain covering much of the Asian continent, the northern part of Australia, and the Indian and western Pacific oceans was the main focus. The full wind nudging shows wind fields closer to the driving global analysis. However, it leads to significantly distorted fields (e.g., temperature and geopotential height) aloft, accompanying excessive precipitation over the western Pacific. The gradual reduction of vorticity nudging from the model top to the ground surface improves rainfall patterns without a discernible distortion of large-scale fields. Further evaluation of a 10-year-summer simulation over East Asia confirmed that this revised SSBC method improves the monsoonal rainfall against the method of Kanamitsu et al. It is therefore concluded that vorticity nudging alleviates largescale errors by maintaining the near geostrophic balance between mass and winds. The reduction of this nudging factor in the lower troposphere allows the ageostrophic component of wind to develop as in nature, which leads to the improvement of precipitation.