A spectral fully compressible nonhydrostatic mesoscale model in hydrostatic sigma coordinates: Formulation and preliminary results

Summary The set of fully compressible nonhydrostatic equations governing a broad spectrum of atmospheric motion was transformed fromz coordinates to sigma coordinates under a hydrostatic base state. The hydrostatic base state may be either time-independent, such as a hydrostatic balance with-out motion or with motion such as a thermal wind balance, or time-dependent such as might be obtained from the result of integrating a hydrostatic model. The transformed set of equations can be used to predict and study all scales of at mospheric phenomena.The set of perturbation equations was also derived under the same condition. The computational sensitivity in computation of pressure gradient force in sigma coordinate can be improved by computing the pressure gradient interms of perturbations under a certain hydrostatic state.The hydrostatic regional spectral model developed by Juang and Kanamitsu (1991) was modified to be a nonhydrostatic spectral model based on the nonhydrostatic equations in sigma coordinates with time-dependent hydrostatic base states. A semi-implicit time integration scheme was used. Two experiments were performed to test this nonhydrostatic spectral model with acceptable results.With 7 Figures     

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.     

Evaluation of an alternative method for numerically modeling nonhydrostatic flows over irregular terrain

Modeling nonhydrostatic atmospheric flow requires the solution of the vertical equation of motion and a prognostic or diagnostic equation for pressure. If the nonhydrostatic components of the flow are relatively small, they can be approximated and incorporated into a purely hydrostatic model, which usually is conceptually simpler and computationally more efficient. A method to do this for a linear model of local thermally-induced circulations is further developed and adapted to a non-linear numerical model of the neutral atmospheric boundary layer. A hydrostatic model and the quasi-nonhydrostatic version were used to simulate neutral flow over simple terrain features. One set of observations taken over a simple change in roughness and another set taken over a change in both roughness and terrain were simulated by both models to assess the capabilities of the quasi-nonhydrostatic technique.It is found that (as expected) the pressure deviation from the hydrostatic state is negligible for the roughness change, but it is an important aspect of neutral flow over terrain. Thus, for flow encountering a simple roughness change, the hydrostatic approximation is good, even for small horizontal scales. However, the quasi-nonhydrostatic model qualitatively produces the features in the observations for flow over a terrain change that the hydrostatic model cannot produce.Journal Paper No. J-12737 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa. Project No. 2779.     

A nonhydrostatic model for simulation of airflow over mesoscale bell-shaped ridges

This paper describes a nonhydrostatic and incompressible mesoscale model formulation using a terrain-following coordinate system. A tensor transformation procedure is used to derive a diagnostic equation for the nonhydrostatic pressure field. The model features a simplified second-order turbulence closure scheme. The two-dimensional version of the nonhydrostatic model, as well as the corresponding hydrostatic model, are applied to simulate stably stratified airflow over mesoscale bell-shaped mountain ridges. The results show that the nonhydrostatic model is capable of simulating nonhydrostatic dynamics of mesoscale lee wave systems such as the trapped wave phenomenon.     

DEVELOPMENT OF MESOSCALE DISTURBANCE IN THREE-DIMENSIONAL INHOMOGENEOUS VORTEX

Starting from the anelastic equations describing deep convection in cylindrical coordinates, the WKBJ method is used to discuss the stability of asymmetric three-dimensional inhomogeneous vortex under the conditions of nonhydrostatic and non-equilibrium gradient wind. From the equation of wave action, the development of disturbance is qualitatively analyzed.     

强风暴的中尺度结构特征及强风暴与暴雨的大尺度结构差异

用滤波方法得到产生强风暴的中尺度天气系统。分析表明:强风暴一般发生在中尺度辐合线、辐合中心及中尺度气旋与反气旋之间,有较强的辐合上升运动,同时还讨论了强风暴与暴雨的大尺度结构差异。      

在大气运动中考虑重力加速度空间变化的问题

给出了在旋转坐标系中考虑重力加速度g的空间变化的大气控制方程组,并证明了大气总质量和总能量的守恒性,以及和取常值g时得到的结果相一致的动、位能和动、内能之间的转换关系.还讨论了在球坐标系中应用方程组时可能出现的困难,给出了在高度近似下在该系中考虑g的空间变化的方案,它可以用来建立完全弹性的非静力模式.     

贵州一次大暴雨过程的中尺度数值模拟与诊断分析

利用三重嵌套的非静力中尺度数值模式MM5V3.5,对2005年5月31日至6月1日贵州省发生的一次大暴雨天气过程进行数值模拟,并利用模拟结果对该过程进行诊断分析。结果表明:模式较好地模拟这次大暴雨过程,并对与暴雨过程相关的中尺度系统的发生发展做出了较成功的模拟,此次过程中,西南涡是造成大暴雨的主要影响系统。对中尺度系统的模拟表明:强降水与强上升运动区及正涡度区有很好的对应关系,低层辐合、高层辐散、西南低空急流、垂直运动增强等是此次暴雨维持和发展的重要机制之一。强降水与水汽辐合的大值密切相关,降水的强弱与辐合的强弱变化一致。     

对流参数在强对流天气潜势预测中的作用

文章结合强对流天气发生发展的物理机制，介绍了与浮力能和风切变有关的几个对流参数的物理概念和计算方法，并结合具体个例分析了这些参数在强对流天气潜势预测中的作用。     

地形追随垂直运动方程在南疆极端暴雨中的诊断分析

针对2021年6月15～17日发生在昆仑山脉北坡的南疆极端暴雨过程,本文综合考虑地形对暴雨发生、发展的作用后,利用地形追随坐标控制方程并采用Boussinesq近似推导建立了地形追随坐标的非静力平衡广义垂直运动方程。诊断结果表明,经向气压梯度力耦合经向散度项（项一）、垂直气压梯度力耦合纬向散度项（项二）和非绝热加热经向梯度项（项三）是激发暴雨垂直运动发展演变的三个主要强迫项。项一体现了偏北风逐渐增强,在昆仑山脉的阻挡下导致经向辐合增强,触发了垂直上升运动。经向气流辐合始终是对流活动最主要的强迫过程,其次为纬向气流辐合。在地形追随坐标形式下,经向和垂直气压梯度能够增强项一和项二。对流发展阶段,水汽辐合与非绝热加热过程增强了非绝热加热经向梯度,促进了垂直上升运动发展。在地形的影响下,对流层中高层西风过山气流波动特征明显。重力波活动导致的高层辐散进一步促进了山脉迎风坡对流活动。经向和纬向气流辐合、非绝热加热过程以及重力波活动等多个因素共同造成了此次南疆极端暴雨。     

PARAMETERIZATION OF MOBILIZATION AND TRANSPORT OF SAND-DUST DURING BLACK STORM AND MESOSCALE NUMERICAL EXPERIMENTS*

The typical black storms or sand-dust storms in the northwestern China are generated and developed through an interaction between the specific large scale circulation pattern and mesoseale systems.The passing by/over a huge sand-abundant desert of a strong cold front with intensive frontal zone at mid and lower levels is a necessary condition for the formation and development of a black storm or a severe sand-dust storm.In order to investigate the mechanism of the sand-dustmobilization,transport and sedimentation during the black or sand-dust storms,a parameterization scheme of sand-dust source-sink terms and an equation of transport for the sand-dust were proposed and incorporated into the MM4 mesoscale model.The modified MM4 model was applied to the "May 1993" black storm case and succeeded in reproducing the evolution of the weather systems associated with the black-storm,the sand-dust concentration at surface layer and its vertical distribution,and the sand-dust sedimentation and transport.Our results show that the numerical simulating method by using a mesoscale model,with inclusion of an equation of the sand-dust transport and a parameterization scheme of the sand-dust source-sink terms,is a promising approach to study the mechanism for sand-dust mobilization,transport and sedimentation during a sand-dust storm event.     

Atmosphere-ocean interaction in mid-latitude storms

Summary Surface fluxes of heat, latent heat, and momentum, and entrainment fluxes and vertical motion at the top of the boundary layer have been calculated for limited regions of several mid-latitude ocean storms. Results have been combined to describe distributions of boundary layer processes which are characteristic of such storms. Surface heat fluxes have important effects in the region west of cold or occluded fronts and are relatively unimportant within a band of about 200 km width east of fronts. Entrainment in pre-frontal regions is driven largely by vertical shear at the top of the boundary layer, while in post-frontal regions it is driven largely by surface heat flux. Boundary layers are well defined in regions more than roughly 200 km east or west of fronts; but closer to fronts boundary layers are not well defined due to the combined effects of entrainment, condensation, and vertical motion associated with the distribution of surface stress.With 12 Figures     

近海台风结构和能量平衡的对比研究

本文以8006和8309两个台风为例，分析了近海台风的结构，指出其重要特点是运动场和热力场的显著不对称和只有弱的暖心。本文还着重对它们在近海及登陆后的总动能和涡动动能收支及加热场情况，进行了对比诊断研究。指出近海台风的不对称结构中，隐含着消亡机制，这是由于不对称的暖区和上升区都偏离台风中心且不相一致，使台风区有效位能向动能转换的过程很弱甚至反向。另外，如台风登陆后继续维持热带系统特点，未能迅速将主要能源由潜热释放转变为斜压过程，也易迅速消亡。在加热场方面，登陆后迅速减弱的台风，最大加热强度较大，其所在层次也较高。     

广东大冰雹风暴单体的多普勒天气雷达特征

选取2004—2012年广东省12个大冰雹风暴单体为样本,利用多普勒天气雷达资料,计算了最大反射率因子及其高度等多个雷达参数,分析了三体散射、旁瓣回波和环境温度层上回波特征以及大冰雹与非冰雹风暴单体间的反射率因子垂直廓线差异。结果表明:大冰雹风暴单体发展均非常旺盛,最大反射因子多超过65 dBZ,对应高度几乎都达到5 km。除受周围大范围雷达回波影响外,大冰雹风暴单体均观测到了三体散射或旁瓣回波特征,并具有一定的预报提前量;在0℃和-20℃层高度上的最大反射率因子均超过54 dBZ。大冰雹风暴单体与非冰雹风暴单体相比,低层回波迅速增加,强核心区垂直伸展更深厚,回波垂直递减率更小。     

准平衡和非平衡流对台风百合（2001）内中尺度深厚湿对流的影响

利用台风百合(2001)的高分辨率模式资料, 应用PV－ω反演方法, 在分析台风准平衡和非平衡垂直环流基本特征的基础上, 诊断研究了热力和动力强迫对台风深厚湿对流结构的影响, 结果指出: 准平衡流能够描述台风中具有较长生命史过程的中尺度对流系统的环流结构, 中层大振幅垂直运动主要由准平衡ω方程中凝结潜热释放的热力强迫所决定; 台风低层由动力辐合强迫产生的弱对流, 对眼墙区深厚湿对流的形成起到触发作用; 而高层动力强迫产生的下沉运动则削弱了凝结潜热项的影响, 抑制了垂直运动向高层伸展。非平衡垂直环流结构上呈现的短波振荡和快时间尺度调整的时空分布特征表明, 它是与重力快波频散过程相联系的非定常环流。分析台风中深厚湿对流形成的物理模型可知, 准平衡动力强迫引起的低层弱对流达到一定振幅, 则会引起中层水汽相变引起的凝结潜热反馈, 使得准平衡和非平衡流叠加形成了眼墙区上升运动的大值中心, 加上非平衡垂直环流对动量和热量起到的补偿和调整作用, 构成了台风内中尺度深厚湿对流的组织化过程。     

The Advanced Regional Prediction System (ARPS) – A multi-scale nonhydrostatic atmospheric simulation and prediction model. Part I: Model dynamics and verification

Summary  A completely new nonhydrostatic model system known as the Advanced Regional Prediction System (ARPS) has been developed in recent years at the Center for Analysis and Prediction of Storms (CAPS) at the University of Oklahoma. The ARPS is designed from the beginning to serve as an effective tool for basic and applied research and as a system suitable for explicit prediction of convective storms as well as weather systems at other scales. The ARPS includes its own data ingest, quality control and objective analysis packages, a data assimilation system which includes single-Doppler velocity and thermodynamic retrieval algorithms, the forward prediction component, and a self-contained post-processing, diagnostic and verification package. The forward prediction component of the ARPS is a three-dimensional, nonhydrostatic compressible model formulated in generalized terrain-following coordinates. Minimum approximations are made to the original governing equations. The split-explicit scheme is used to integrate the sound-wave containing equations, which allows the horizontal domain-decomposition strategy to be efficiently implemented for distributed-memory massively parallel computers. The model performs equally well on conventional shared-memory scalar and vector processors. The model employs advanced numerical techniques, including monotonic advection schemes for scalar transport and variance-conserving fourth-order advection for other variables. The model also includes state-of-the-art physics parameterization schemes that are important for explicit prediction of convective storms as well as the prediction of flows at larger scales. Unique to this system are the consistent code styling maintained for the entire model system and thorough internal documentation. Modern software engineering practices are employed to ensure that the system is modular, extensible and easy to use. The system has been undergoing real-time prediction tests at the synoptic through storm scales in the past several years over the continental United States as well as in part of Asia, some of which included retrieved Doppler radar data and hydrometeor types in the initial condition. As the first of a two-part paper series, we describe herein the dynamic and numerical framework of the model, together with the subgrid-scale turbulence and the PBL parameterization. The model dynamic and numerical framework is then verified using idealized and realistic mountain flow cases and an idealized density current. Other physics parameterization schemes will be described in Part II, which is followed by verification against observational data of the coupled soil-vegetation model, surface layer fluxes and the PBL parameterization. Applications of the model to the simulation of an observed supercell storm and to the prediction of a real case are also found in Part II. In the latter case, a long-lasting squall line developed and propagated across the eastern part of the United States following a historical number of tornado outbreak in the state of Arkansas. Received April 14, 2000 Revised July 17, 2000     

A climatic study of severe storms over Bulgaria produced by Mediterranean cyclones in 1990−2001 period

Over south-eastern Europe, severe weather events are often associated with Mediterranean cyclones. This paper presents a climatic study of severe storms over Bulgaria produced by synoptic-scale Mediterranean cyclones, which are the main high-impact weather systems for the region during the winter season. The study is based on a synoptically oriented data set that contains systematic information about the pronounced Mediterranean cyclones including their life and trajectories over the Mediterranean area and the severe storms over Bulgaria produced by them. The definition of a severe storm is a storm in which the wind speed and precipitation exceed thresholds of 15 m/s and 30 mm/24 h, respectively. The observed severe storms were split into two groups by the number of districts where they have caused damages. During the last years a trend of decreasing numbers of initiating synoptic-scale Mediterranean cyclones has been observed. However, the number of those producing a high-impact weather phenomenon over Bulgaria has increased. In these high-impact cases, the observed paths of the cyclones are connected to the specific circulation conditions over the region. In the most severe cases, almost 80% of the cyclones move through the southernmost parts of Balkan Peninsula and for a large portion of them, this motion is associated with a blocking regime in the mid-level mass field.The development of such Mediterranean cyclone cases has been studied and results are presented here. The results illustrate the synoptic-scale mechanisms for intensification of a feeding flow of Mediterranean air towards the area affected by extremely severe weather.     

一次热带风暴造成辽宁暴雨的数值模拟研究

利用非静力平衡的中尺度模式MM5V3,对1997年8月20日Winnie台风减弱的热带风暴,进行了48h数值模拟。结果表明：MM5不仅比较好地模拟出高空形势场,还成功地模拟出了强降水的时间和强度。应用模拟出的物理量对这次过程进行了诊断分析,并对暴雨形成机制进行了讨论。     

2010年春季北京地区强沙尘暴过程的微气象学特征

利用北京大学校园地区PM10质量浓度观测资料、中国科学院大气物理研究所325m气象塔气象要素梯度和湍流观测资料,分析了北京地区2010年3月20～22日两次强沙尘暴过程微气象学要素和沙尘参量的时空演变以及湍流输送特征,为理解北京地区强沙尘暴天气沙尘输送规律和微气象学特征提供参考。结果表明:3月20～22日强沙尘暴过程前后不同高度温度先升后降,气压和相对湿度则相反。强沙尘暴来临时,高层风速先迅速增大,低层风速增加略有滞后,风切变明显加强,PM10浓度最大值和风速极大值出现时间较吻合。强沙尘暴过境时,不同高度向下的湍流动量输送、向上的湍流热量输送和湍流动能明显加强。与3月21日非沙尘暴日相比,强沙尘暴过程湍流动量通量增加,有利于沙尘粒子的水平和垂直输送过程;由于冷锋过境,水平热通量增大;垂直热通量因白天温度垂直梯度减小而减小,夜间因逆温层被破坏而增加;水平湍流动能对湍流动能占主要贡献,垂直湍流动能仅占水平湍流动能的10%～25%。