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.     

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

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

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

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

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

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

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.     

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     

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.     

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

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

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

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

GRAPES模式标准初始化方案设计与实现

文章介绍了GRAPES模式标准初始化系统的方案设计及程序结构，对3种不同方案产生的初始场进行了分析。结果表明:水平风场、温度场、位温场对插值方法的精度并不是很敏感，而高度场以及气压Exner函数则对垂直插值方法的精度要求很高；在3种模式变量产生方案中，方案一的效果较差，但它计算速度快，在地形梯度不大的地方可以选择使用；方案二与方案三的效果则比较好，用户可根据自己的需求任意选择。GRAPES模式标准初始化系统在完成其主要功能的同时，通过灵活的程序设计，对标准初始化系统的区域、分辨率、垂直层次、水平插值方法、垂直插值方法等的任意选取，使得标准初始化系统可以很容易满足模式初值的各种需求，为模式开发研究实验提供了各种便利。     

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全球集合预报业务系统的控制预报初值。     

SL-AV atmospheric model version using σ-<Emphasis Type="Italic">p</Emphasis> hybrid vertical coordinates

The SL-AV atmospheric model version using hybrid vertical coordinates combies the advantages of sigma and isobaric coordinates. The formulation and discretization of model equations maintain the equivalency of the new model version to the basic sigma version in the special case, when hybrid coordinates coincide with sigma coordinates. The SL-AV model version with hybrid vertical coordinate is verified with medium-range weather forecasts. The decrease in the errors of predicted geopotential height and wind as compared to the sigma model version is demonstrated. The use of hybrid coordinates also leads to a certain increase in forecast skill scores for some meteorological parameters characterizing aviation significant weather.     

On the parameterization of the vertical velocity at the top of planetary boundary layer

In this paper, an equation of the vertical velocity at the top of PBL is derived by use of a PBL model which is based on an analytic and actual form of K. Results show that the vertical velocity is a function of geostrophic vorticity, geostrophic wind speed, Coriolis parameter and the roughness of the ground, thus improving Charney-Eliassen’s formula. The order of magnitude of the vertical velocity computed from our equation is in agreement with that from the latter, but more factors affecting the vertical velocity are included.     

An ice breeze mechanism for boundary-layer jets

The existence of a low-level (z=~1000 m) jet adjacent to a sea-ice boundary is investigated with a two-dimensional numerical model. A thermally-direct ice breeze circulation is induced by specifying an ice-sea surface temperature gradient, with the mean geostrophic wind parallel to the ice edge. Pressure changes associated with over-water mixed-layer development create an increase in geostrophic velocity that accounts for most of the increase in wind speed. A change in initial geostrophic wind direction has significant effects on location and intensity of the low-level jet; geostrophic winds parallel to the ice edge result in stronger jets than occur with cross-ice geostrophic winds. An inertial oscillation simulated by the model in 1-D makes a negligible contribution to the low-level jet.     

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.     

垂直变量配置对静力适应过程的影响

为了研究垂直变量配置对静力适应过程的影响,本文从描写静力适应过程的方程组出发,分别在将所有变量置于整层上的非跳点N网格;将垂直速度和温度放置在整层,水平速度、气压和密度等变量放置在半层的Charney Phillips跳点网格(CP网格);将水平速度、气压和温度放置在整层,将垂直速度和密度放置在半层的Lorenz跳点网格(L网格);将密度变量放置在整层的Charney Phillips跳点网格(CP_N网格);将密度放置在整层的Lorenz跳点网格(L_N网格)上进行离散,垂直格距分1 km、0.5 km、0.2 km和0.01 km,研究了在这5种网格上产生的频率和垂直群速的相对误差。结果表明:(1)L_N网格和CP网格是完全等效的两种网格。(2)不论垂直格距为多少,CP网格和L网格的误差都是最小,N网格次之,CP_N网格的误差最大。(3)随着垂直格距的减少,在这几种网格上产生的误差都在减小。对于CP网格、L网格和N网格,在水平长波和垂直短波处产生的误差较大。而CP_N网格对水平波长变化不敏感,垂直波长越短,误差越大。(4)当垂直格距为0.01 km时,这几种网格都对水平波长的变化不敏感了,仅对垂直波长敏感。(5)CP网格、L_N网格和L网格在描写静力适应过程和斜压地转适应过程都是误差最小的垂直变量配置方案,因此在非静力完全可压缩深层大气数值预报模式中应优先选择这3种方案。