一个双波地形重力波拖曳参数化方案

当地形次尺度强迫的作用与显式的经典动力作用效应相当时,地形重力波拖曳力对于环流的维持,以及动量和热量通量输送的动力效应变得十分显著。这种地形次尺度拖曳作用项可通过参数化的方法,在动力方程中加入额外的小项而引入数值模式。目前成熟的地形重力波拖曳参数化方法,如第1代基于线性单波理论的参数化方案;以及侧重考虑了临界层作用等因素对拖曳力的额外贡献的第2代参数化方案,都无法有效表达风速垂直变化引起的波动应力随高度变化的特征。基于上述考虑,本文给出了一个双波参数化方案用于计算地形重力波拖曳中由线性自由传播重力波造成的波动应力的垂直分布。通过二阶WKB近似,它对由风速垂直变化引起的对波动应力的选择性临界层吸收过程和经典的临界层吸收过程做了显式处理;而在不发生临界层吸收现象的地区,则用两个单波同时在垂直方向上进行应力的传播,并利用波饱和标准进行应力耗散。进一步地在真实地形(以大别山地区为个例)条件下的测试结果表明,通过在不同理想风速廓线以及北半球冬季中纬度纬向平均风廓线下对波动应力垂直分布的计算,证明该方案确实能有效地给出应力随高度变化的特征。     

青藏高原地形重力波拖曳的初步分析及数值模拟研究

针对目前对青藏高原大地形激发的重力波拖曳相关问题还不十分清楚,在GRAPES_Meso模式中引入次网格地形重力波拖曳参数化方案,通过数值试验初步研究了青藏高原地区次网格地形重力波拖曳的一些相关参数,结果指出:(1)沿30°N地形重力波拖曳的垂直分布显示,阻塞拖曳主要存在于模式的低层(第1—5层),重力波拖曳主要存在于模式的第5—10层;从水平分布看,模式第3层以阻塞拖曳为主,主要位于青藏高原边缘地区,阻塞拖曳大值区沿喜马拉雅山脉走向和青藏高原东坡;模式第5层以重力波拖曳为主,主要位于青藏高原东部地区和云贵高原的北部边缘。(2)弗劳德数和气流绕流高度分析表明,在青藏高原喜马拉雅山脉一带和高原东部边缘地区,气流爬坡能力强,同时在这一地区绕流高度最高;弗劳德数越大的地区绕流高度距离地表越高。(3)采用次网格地形重力波拖曳参数化方案后,对于低层和高层地形重力波破碎的发生有更准确的描述,地形重力波是向上垂直传播的。(4)个例和批量试验检验结果表明,采用次网格地形重力波拖曳参数化方案对于风场和降水模拟有正效果,提高了模式预报的准确率。     

对流云街激发的重力波和波动阻力

大气边界层中的对流活动,可以在其上部稳定层中激发出重力波,并引起垂直动量输送,影响到对流层和平流层中的动量平衡过程.从二层模式中大气波动方程的线性解出发,得出了对流云街激发的重力波波阻解析表达式,并讨论了大气条件对波阻的影响.这些分析可有助于大气环流模式(GCM)中此类重力波波阻参数化表达式的建立和改进.     

Test of a convectively forced gravity wave drag parameterization in a general circulation model

The influence of gravity wave drag induced by cumulus convection (GWDC) on a simulated boreal summer climate was evaluated in a general circulation model. For this, the GWDC scheme developed by Chun and Baik was implemented into a version of the National Centers for Environmental Prediction (NCEP) global spectral model (GSM). Ensemble simulations with the two different convection schemes, the simplified Arakawa-Schubert (SAS) scheme and Community Climate Model (CCM) convection scheme, were conducted for the boreal summer of 1996. A cloud factor to modulate the stress intensity with respect to the cloud type was introduced in this study, in order to prevent unrealistic behaviors of the GWDC scheme in GSM. The effect of gravity wave drag on the zonal mean of wind and temperature fields was focused. On the whole, the effect of GWDC in this study is positive on the simulated seasonal climate. It is evident that biases in temperature in the polar region as well as in the zonal and meridional winds in the upper atmosphere are reduced. The percentage of reduction of the bias in zonal winds is about 10–20%. Such a response of the GWDC forcing widely appears not only in tropical regions but also in mid-latitude regions. These characteristics are prominent in the case of the SAS scheme, which is due to the various convective cloud types. The magnitude of GWDC forcing is generally small, but still positive, in the case of the CCM scheme, which is due to rather homogeneous cloud types. It is also found that the role of a particular GWDC forcing depends upon the inherent systematic biases of a particular model. It is concluded that incorporation of the GWDC parameterization in GCMs should be taken into account to improve the seasonal prediction.     

A NOTE ON THE APPLICATION OF LINEAR WAVE THEORY AT A CRITICAL LEVEL

The purpose of this note is to estimate the accuracy and practical limitations of applying linear theory at a critical level over a realistic range of atmospheric stabilities for an idealized surface terrain. These estimates are made by comparing the results of a linear model with a nonlinear numerical model at a critical level. Essentially similar results are obtained from each model for wave stress, wave breaking height and wave dissipation through the critical level. Because gravity waves can be either evanescent or internal depending on the relative sizes of the Scorer parameter and the wavenumber of the ground surface disturbance, the somewhat paradoxical result develops that wave breaking and non-linearity increase with increasing bulk Richardson number. It is recommended that steady linear wave theory be used in gravity wave drag parameterizations provided near real time profiles of background velocity and temperature are available.     

Diagnosis of the forcing of inertial-gravity waves in a severe convection system

The non-hydrostatic wave equation set in Cartesian coordinates is rearranged to gain insight into wave generation in a mesoscale severe convection system. The wave equation is characterized by a wave operator on the lhs, and forcing involving three terms—linear and nonlinear terms, and diabatic heating—on the rhs. The equation was applied to a case of severe convection that occurred in East China. The calculation with simulation data showed that the diabatic forcing and linear and nonlinear forcing presented large magnitude at different altitudes in the severe convection region. Further analysis revealed the diabatic forcing due to condensational latent heating had an important influence on the generation of gravity waves in the middle and lower levels. The linear forcing resulting from the Laplacian of potential-temperature linear forcing was dominant in the middle and upper levels. The nonlinear forcing was determined by the Laplacian of potential-temperature nonlinear forcing. Therefore, the forcing of gravity waves was closely associated with the thermodynamic processes in the severe convection case. The reason may be that, besides the vertical component of pressure gradient force, the vertical oscillation of atmospheric particles was dominated by the buoyancy for inertial gravity waves. The latent heating and potential-temperature linear and nonlinear forcing played an important role in the buoyancy tendency. Consequently, these thermodynamic elements influenced the evolution of inertial-gravity waves.     

不同阶曲线拟合扰动场对下平流层重力波气候特征影响研究

重力波参数气候特征是确定大气模式中重力波参数化方案的重要条件之一,高垂直分辨率探空资料扰动场是获取重力波参数气候特征的基础数据;目前,获取扰动场的方法较多,但基于不同方法计算的扰动场对重力波参数气候特征影响的研究较少。基于2014—2017年山西太原气象台高垂直分辨率探空资料,利用2—4阶曲线拟合方法获取下平流层（17—24 km高度）温度扰动场、纬向风扰动场和经向风扰动场,经统计发现2阶与3阶曲线拟合方法的扰动场相似程度较高;在此选取相似度较高的2阶、3阶曲线拟合方法的扰动场分别计算大气重力波参数,并对大气重力波参数间的气候差异特征进行研究。结果表明:（1）不同阶曲线拟合方法扰动场的变化振幅及随高度变化趋势存在差异,且扰动场间的相关较弱;（2）2阶、3阶曲线拟合方法扰动场得到的重力波参数大小、年内变化趋势及在不同区间范围内占有率均存在差异,且相关较弱;（3）1—12月,相对3阶曲线拟合方法的扰动场,基于2阶曲线拟合方法的扰动场得到的重力波群速、水平波长、垂直波长、周期、固有相速均较大,而重力波能量上传百分比在某些月份较大。因此,不同阶曲线拟合方法扰动场间存在差异,会导致计算得到的大气重力波参数气候特征存在差异,最终对研制大气模式中的大气重力波参数化方案产生影响。      

Gravity Wave Momentum Fluxes and Surface Pressure Drags due to the Pyrénées: Variation with Numerical Model Horizontal Resolution

Summary. ?A hydrostatic numerical model is used to simulate the lee wave event IOP3 (0000 GMT to 1200 GMT 15^th October 1990) from the PYREX mountain experiment. Results from integrations at different horizontal resolutions are used to investigate the effect on surface pressure drag and the vertical flux of horizontal momentum due to orographically forced gravity waves. In particular, results showing the dependence on resolution of the partitioning between resolved and parametrized wave drag and fluxes are presented. With the model horizontal gridlength changing from 50 km to 10 km the majority of wave momentum flux changes from being parametrized to becoming resolved. More significantly, there is a change in the profile of flux with height. At 50 km resolution the largest inferred mean flow decelerations are at lower stratospheric level due to the parametrization scheme. At 10 km resolution this is shifted, with less deceleration high up and more wave deceleration lower down in the troposphere. Numerical weather prediction models are now beginning to take account of such low level drag with beneficial results. Received March 2, 1999/Revised July 15, 1999     

Parametrization and influence of subgridscale orography in general circulation and numerical weather prediction models

Summary The problem of representing the drag due to subgridscale orography is examined. Results from model simulations are used to illustrate clear deficiencies in the global angular momentum budgets and possible ways of correcting for these deficiencies are considered. It is argued that a formulation for the stress due to subgridscale gravity waves is required, as was first recognized by Boer et al. (1984a, b), and the impact of a scheme based on Palmer et al. (1986) is presented. The scheme is improved by using directionally-dependent subgridscale orographic variances. Results from 90-day integrations using the ECMWF and UKMO models with similar resolutions are very much in accord both with and without a wave drag parametrization scheme; showing much improved wintertime circulations.The relationship between the wave drag and the model orography is examined with the use of idealized stress profiles whereby the wave drag is limited to either near the surface or in the stratosphere. A combination of parametrized wave drag with an envelope orography performs best at this stage of development.Results are presented from a substantial series of ten-day forecast experiments with the ECMWF operational model using mean and (1 ) envelope orographies, these show significant improvements in forecast skill.With 26 Figures     

地形重力波拖曳参数化对热带气旋强度和路径预报影响的研究

本文在GRAPES_TMM（Global/Regional Assimilation and Prediction System for Tropical Mesoscale Model）——中国南海台风模式版（面向南海和东南亚）中发展和引进了KA95（Kim and Arakawa,1995）地形重力波拖曳参数化方案（GWDO）,并对2012年主要的9个登陆台风进行了试验对比研究,考察了不同标准Richardson数（Ri_c）的GWDO试验对台风路径和强度预报的影响。结果表明,在引入地形重力波拖曳参数化过程后,模式对台风登陆时路径和强度的预报能力均要有提高,对台风预报时长越长,GWDO的影响也更为显著。对双台风“SAOLA”和“DAMREY”试验结果表明,GWDO对台风外围距台风中心150 km的对流层中下层风速减弱较为明显,减弱了GRAPES区域模式对台风强度预报偏强的现象,对台风强度长时间预报改善更为明显。不同标准Ri_c对重力波拖曳力的计算较为敏感,当Ri_c取1.0时,动能迅速的在低层被频散,能量无法有效地上传;Ri_c取0.25时,大部分的能量在中高层被频散。总的来说,Ri_c取0.75时对台风路径和强度预报改进更为显著,其结果可为业务预报提供指导意义。     

Generation mechanisms of convectively induced internal gravity waves in a three-dimensional framework

The generation mechanisms of convective gravity waves in the stratosphere are investigated in a three-dimensional framework by conducting numerical simulations of four ideal storms under different environmental conditions: one un-sheared and three constant low-level sheared basic-state winds with the depth of the shear layer of 6 km and the surface wind speeds (U_s) of 8, 18, and 28 m s^?1, using the Advanced Regional Prediction System (ARPS) model. The storms simulated under the un-sheared (U_s = 0 m s^?1), weakly sheared (U_s = 8 and 18ms^?1), and strongly sheared (U_s = 28ms^?1) basicstate winds are classified into single-cell, multicell, and supercell storms, respectively. For each storm, the wave perturbations in a control simulation, including nonlinearity and microphysical processes, are compared with those in quasi-linear dry simulations forced by diabatic forcing and nonlinear forcing that are obtained from the control simulation. The gravity waves generated by the two forcing terms in the quasi-linear dry simulations are out of phase with each other for all of the storms. The gravity waves in the control simulation are represented by a linear sum of the wave perturbations generated by the nonlinear forcing and diabatic forcing. This result is consistent with the results of previous studies in a two-dimensional framework. This implies that both forcing mechanisms are important for generating the convective gravity waves in the three-dimensional framework as well. The characteristics of the three-dimensional gravity waves in the stratosphere were determined by the spectral combination of the forcing terms and the wave-filtering and resonance factor that is determined from the basic-state wind and stability as well as the vertical structure of the forcing.     

青藏高原边坡地形重力波拖曳方案应用的个例研究

本文针对我国西南地区的一次持续性暴雨过程,利用WRF模式引入地形重力波拖曳方案进行120h的模拟,对比结果认为引入地形重力波拖曳的模拟效果总体看来略有优势,并详细分析了不同模拟时段的位势高度场,风场,海平面气压场与降水,其结果认为,引入与未引入地形重力波模拟的差异主要出现在48~72h之后,位势高度场方面,72h之后,在青藏高原主体上引入拖曳方案后的模拟结果相对未引入拖曳方案的结果是正偏差,高原下游地区为负偏差,引入后的模拟相对准确;风场模拟方面同样72h以内的模拟差距较小,72h之后高原上出现一处气旋性涡旋,其中心位置的模拟,引入拖曳方案后的结果相对准确,本次个例中出现的一例西南低涡,两实验都能够较好的模拟出700h Pa的气旋性涡旋;降水模拟的结果表明,在72h之后,对于四川地区降水中心位置的模拟,引入拖曳方案的模拟有所提高,但中心依然过强。这方面主要考虑降水的模拟牵涉模式内部更为复杂的处理过程和其他因素,地形重力波拖曳只是其中一个影响因素而不是核心因素。另外此方案对温度场及海平面气压场模拟的影响不明显。理论上讲引入地形重力波拖曳效应的模拟是物理过程相对更为完善的,在本次个例中模拟效果也略优,但在高原下坡复杂地形处长期应用的模拟效果仍需要更多的测试,将在后续工作中加以完善。      

A Dynamic Urban Canopy Parameterization for Mesoscale Models Based on Computational Fluid Dynamics Reynolds-Averaged Navier–Stokes Microscale Simulations

Urban canopy parameterizations (UCPs) are necessary in mesoscale modelling to take into account the effects of buildings on wind and turbulent structures. This study is focused on the dynamical part of UCPs. The main objective is twofold: first, computing important UCP input parameters (turbulent length scales and the sectional drag coefficient) by means of Reynolds-averaged Navier–Stokes (RANS) simulations of turbulent flow; and second, comparing UCP variables with spatially-averaged properties obtained from RANS simulations for the same configurations. The results show the importance of using a suitable parameterization of the drag force for different packing densities. An urban canopy parameterization that is a compromise between simplicity and accuracy is proposed. This scheme accounts for the variation of drag coefficients with packing densities, and has a parameterization of turbulent length scales. The technique adopted ensures that, at least for the simple configurations studied, the urban canopy parameterization gives values of spatially-averaged variables similar to those computed from a more complex simulation, such as RANS that resolves explicitly the flow around buildings.     

Forced Rossby wave packets in barotropic shear flows with critical layers

In this paper, we investigate the meridional propagation of a forced Rossby wave packet towards a critical layer in a zonal shear flow by solving the linearized barotropic vorticity equation. The forcing is applied north of the critical layer. Two approaches are employed for solving this problem. First, an analytic solution valid for large time is derived, using Fourier and Laplace transform techniques and asymptotic approximations. This solution exhibits the modification due to the wave packet of the solution obtained by Warn and Warn (1976) [Warn, T., Warn H., 1976. On the development of a Rossby wave critical level. J. Atmos. Sci., 33, 2021–2024.] in the monochromatic case. A numerical investigation is then carried out using a finite difference scheme and a time-dependent radiation condition. It is found that the forced wave packet is absorbed at the critical layer and the total momentum transferred to the mean flow as a result of the absorption is observed to be proportional to the length scale of the wave packet. We also consider the case of a north–south mean flow with a longitudinally propagating wave packet forced to the east or west of the critical layer. The monochromatic version of this problem has been used before (Geisler, J.E., Dickinson, R.E., 1975. Critical level absorption of barotropic Rossby waves in a north–south flow. J. Geophys. Res., 80, 3805–3811.) to examine the interaction of western boundary currents and oceanic Rossby waves.     

Generation of oceanic internal gravity waves by a cyclonic surface stress disturbance

Atmospheric cyclones with strong winds significantly impact ocean circulation, regional sea surface temperature, and deep water formation across the global oceans. Thus they are expected to play a key role in a variety of energy transport mechanisms. Even though wind-generated internal gravity waves are thought to contribute significantly to the energy balance of the deep ocean, their excitation mechanisms are only partly understood.The present study investigates the generation of internal gravity waves during a geostrophic adjustment process in a Boussinesq model with axisymmetric geometry. The atmospheric disturbance is set by an idealized pulse of cyclonic wind stress with a Rankine vortex structure. Strength, radius and duration of the forcing are varied. The effect upon wave generation of stratification with variable mixed-layer depth is also examined.Results indicate that internal gravity waves are generated after approximately one inertial period. The outward radial energy flux is dominated by waves having structure close to vertical mode-1 and with frequency close to the inertial frequency. Less energetic higher mode waves are observed to be generated close to the sea floor underneath the storm. The total radiated energy corresponds to approximately 0.02% of the wind input. Deeper mixed-layer conditions as well as weaker stratification reduce this fraction.The low energy transfer rates suggest that other processes that drive vertical motion like surface heat fluxes, turbulent motion, mixed region collapse and storm translation are essential for significant energy extraction by internal gravity waves to occur.     

对流边界层中过山气流的数值模拟

采用ARPS4.0非静力中尺度气象模式模拟了对流边界层中气流过山引起的地形波,讨论了地形及大气条件改变对其的影响.模拟表明,当大气边界层是对流边界层时,气流过山引起的地形强迫,仍能在上部稳定层结中造成足够的垂直扰动,产生向上传播的重力内波,重力内波引起的波动阻力仍不可忽略.     

IMPROVED FORECASTING OF COLD AIR OUTBREAKS OVER SOUTHERN CHINA THROUGH OROGRAPHIC GRAVITY WAVE DRAG PARAMETERIZATION

The parameterization of gravity wave drag induced by sub-grid scale orography(GWDO), which has been used in the regional model based on the Global/Regional Assimilation and Prediction System for Tropical Mesoscale Model(GRAPES_TMM), is extended to include the effect of mountain flow blocking drag(MBD). The extended scheme is evaluated against non-GWDO parameterization, including a cold air outbreak over southern China and a monthly verification in February 2012. The experiment results show that the GWDO and MBD parameterization both improves the forecasting of the cold air outbreaks over southern China, as well as alleviations of system bias of GRAPES_TMM.(1) The extended scheme alleviates the strong southerly wind and high surface temperature simulation during the cold air outbreak, especially over northern Guangxi and Guangdong(NGG) province, where local high surface temperature simulation reduces nearly 5 degree.(2) The MBD parameterization improves southerly wind simulations over NGG, as well as surface temperature forecasts improvement over Guangxi, Guizhou province and southern Yunnan-Guizhou plateau(YUP), and low level southerly wind simulation improvement over intertidal zone over south China.(3) The formation of MBD is mainly in the mountain area(Wuyi, Daba mountain, east of YUP) and coastal area. The MBD over plateau, which is mainly formed at the west of 105°E, is stronger and thicker than that over Nanling mountain.(4) The improvement of GWDO and MBD parameterization is stable in model physics. MBD parameterization demonstrates more overall improvements in the forecasts than GWDO, and the larger of the model forecast error is, the greater improvements of MBD contribute to. Overall, the extended GWDO scheme successfully improves the simulations of meteorological elements forecasting during cold air outbreaks.     

Momentum flux of stratospheric gravity waves generated by Typhoon Ewiniar(2006)

The momentum flux of stratospheric gravity waves generated by Typhoon Ewiniar (2006) is examined using a Weather Research and Forecasting (WRF) model. In the stratosphere, zonal momentum flux with a positive sign by eastward-propagating waves is significant during the northward moving of the typhoon, while both zonal and meridional momentum fluxes with positive signs are significant during the typhoon decaying stage in which the typhoon moves northeastward. The magnitude of the momentum flux is greater during the mature stage of the typhoon than the decaying stage, and the phase speeds of the dominant momentum flux are less than 30 m s^?1 with a peak at 10–16 m s^?1. Positive momentum flux decreases with height overall in the stratosphere for both zonal and meridional components. The resultant gravity-wave drag forcing plays a role to decelerate the easterly background wind in the stratosphere. This drag forcing is relatively large above z = 40 km and below z = 20 km, and lower stratospheric wave drag is expected to affect the typhoon dynamics by modifying the background wind shear and inducing the secondary circulation in the troposphere.     

中尺度大气波动的波谱和谱函数——数学模型和计算方法

作者得到了准二维Boussinesq方程组，并用其研究了中尺度大气波动的波谱和谱函数。在一定条件下对该方程组线性化并取标准模后，可将其初边值问题转化为矩阵的广义特征值问题来进行数值求解，这样就可知原问题波谱和谱函数的性质。当无基本流且取地转参数、层结参数为常数时，可求得其波谱和谱函数的解析解。此时该模式中仅包含有一对重力惯性内波模态，且各模态均是简谐波；模态越高，垂直波数越大则波动传播得越慢，所有的模态均为离散谱，并存在聚点。对此作者用数值解作了验算，结果表明，该数值求解方案合理可行，对不太高的模态其精度也令人满意。在无基本流然而考虑层结的垂直变化后，则一般无法求取解析解，为此进行了数值求解。这时该模式仍仅包含有一对重力惯性内波的离散谱模态，不过由于层结参数的变化，各模态结构与简谐波出现了偏差。