对流层上传重力波的非线性演化

利用二维全隐欧拉格式对重力波在可压、非等温大气中的非线性传播过程进行了数值模拟和分析.分析结果表明，从对流层顶激发的重力波能稳定地经平流层传到中层顶，从而将能量和动量从一个区域带到另一个区域；在向上传播过程中，重力波经历了发展、位温翻转、对流直至最终破碎的演变；重力波的破碎是对流和小尺度波动的重要的源，对流不稳定和翻转是非线性现象的一个基本特征.计算还显示，扰动源的大小直接影响着重力波的非线性传播过程，当扰动源足够小时，重力波能稳定传播，而大振幅扰动可以加速重力波的破碎.     

重力波波包在向上传播过程中的破碎

采用二维全隐欧拉（FICE）格式对具有高斯分布的重力波波包在可压大气中传播时的饱和过程进行数值模拟和分析．数值计算结果表明,波振幅首先随高度增加而增加,但当波振幅接近于线性不稳定性给出的阈值时,不再增加,重力波波包达到饱和进而破碎．破碎出现的高度（86．50km）比线性理论预言的结果（84．59km）要高一些,并且一般都在波包的下游出现．波破碎过程能使波能量在空间重新分配,并对重力波能量有明显的耗散作用．并且波破碎会使波相关能量传输方向偏离线性射线理论的射线路径．     

重力波包在中层大气温度波导中传播的数值模拟

给出了重力波在中层大气温度波导中的反对称形式导制传播的线性理论模型，并采用二维非线性的数值模型对重力波波包在中层大气温度波导中的传播和演变过程进行了模拟研究.模拟的结果表明，下层大气激发的重力波能量进入波导区域后被俘获形成导制传播.重力波在波导内不停地来回反射，垂直方向的自由传播受到限制，能量在波导内沿着水平方向传输，模拟得到的波参数与理论值相近.重力波包在温度波导中传播时伴随着能量泄漏，且能量泄漏的速率随时间变缓，最终总有部分能量被限制在波导区域.重力波在水平方向上传播几百公里后，依然维持着良好的波结构，同时数值模拟也给出了重力波在波导区域内能量密度的时空分布.     

重力波非共振激发的数值研究

应用二阶算子分裂格式，采用全非线性模式，对可压大气中重力波波包非共振激发进行数值模拟和分析.研究表明：两个向上传播的重力波波包，通过非共振相互作用，也能激发出一个向下传播的重力波波包.生成波的能量几乎正比于主波的能量，也随次波的能量增大而增大.能量传递方向主要表现为从高频的主波向次波和生成波转移，这与共振相互作用不完全相同.非共振作用的特征时间由主次波的初始位置和波长频率决定，与主次波的初始振幅无关.生成波的波矢量并不严格满足弱非线性相互作用理论预言的波矢量匹配条件，而且生成波的波长和频率在传播的过程中会随时间发生变化，这表明了弱非线性相互作用理论的局限性.     

中层耗散大气中Rossby波包的非线性相互作用理论

采用弱非线性近似得出中层耗散大气连续谱Ｒｏｓｓｂｙ波包的非线性时空演化方程，讨论了Ｒｏｓｓｂｙ波包的三波相互作用问题．数值计算表明，耗散和非线性的共同效应决定了Ｒｏｓｓｂｙ波包的演变．当一个Ｒｏｓｓｂｙ波包通过大气传播时，它的振幅若超过某个阈值，空间尺度分别比它大和比它小的两个次级Ｒｏｓｓｂｙ波包的振幅会随时间增长．特别当这两个次级波包同时随时空变化时，仅当主波的振幅超过一个更大的阈值，且其群速度介于两次级波包的群速度之间时，两次级波包的振幅才会随时空同时增长，即出现绝对不稳定现象，耗散和３个波包的频率失配都会增大不稳定的阈值．     

损耗大气中重力波的非线性相互作用理论

基于中高层大气重力波动力学是由非线性过程和损耗过程共同决定的物理思想,本文采用弱非线性相互作用近似,推导出损耗大气中重力波的非线性相互作用方程.这组方程是研究固定相位和随机相位重力波相互作用问题的出发点.通过引入平均振幅,我们得到了损耗情况下离散重力波的三波相互作用方程,该方程描述了重力波波包非线性时空演变的规律.作为该方程的一个具体应用,我们考虑了由于波-波相互作用产生的不稳定性.当一大尺度大振幅的主重力波通过大气传播时,非线性相互作用可能导致两个次级波振幅随时间指数增长.由于分子损耗和频率失配,主波的振幅必须大于一个阈值,这种指数增长才可能出现.共振条件满足时,阈值变为最小.频率失配还会导致次级波本征频率发生改变,改变的大小是频率失配值的一半.     

重力波与潮汐之间的非线性相互作用

采用二维可压缩大气中重力波非线性传播的数值模式, 研究了重力波与潮汐之间的非线性相互作用. 结果表明, 重力波在潮汐背景中传播时, 先后在z = 75~85, z = 90~110和z = 15~130 km 3个高度上发生不稳定. 垂直波长首先由12 km变成27 km左右, 新生成的长波逐渐被压缩并再次生成20 km左右的长波. 长波和短波出现的高度分别对应于反向和同向背景风场区域的高度. 在重力波主要的破碎区域(90~110 km)以上, 仍有部分重力波继续上传. 重力波在上传过程中除了对背景风场加速之外, 还增大了潮汐的振幅, 特别是在重力波发生不稳定之后, 对潮汐振幅的放大作用更加明显.     

中层大气Rossby波与惯性重力波的非线性相互作用 （Ⅰ）相互作用方程

从包含Ｒｏｓｓｂｙ波和惯性重力波的大气运动方程组出发，采用弱非线性相互作用近似，推导出耗散大气中这两种尺度相差很大的波动之间的非线性相互作用方程．以此为基础，得到了描述窄角谱Ｒｏｓｓｂｙ波包和惯性重力波包的非线性时空演变规律的三波相互作用方程．数值分析表明，当一个Ｒｏｓｓｂｙ波包与两个惯性重力波包发生相互作用时，两个惯性重力波包之间进行快速的能量交换，同时与Ｒｏｓｓｂｙ波包之间进行缓慢的能量传输．从时间尺度上讲，惯性重力波可以看作Ｒｏｓｓｂｙ波包运动的背景噪声，因此上述非线性相互作用过程可以理解为大尺度Ｒｏｓｓｂｙ波包与背景噪声之间的能量交换过程．     

一次冰雹过程的惯性重力波观测及数值模拟

使用高灵敏度的电容式微压波传感器对1998年4月11日16时发生在贵州省普定县的一次降冰雹过程的重力波进行观测,利用WRF（Weather Research and Forecast）中尺度模式对这一过程进行数值模拟,使用Morlet小波方法对模拟结果进行分析,得出这一过程中惯性重力波的分布和变化规律,并分析急流、地形及切变线对惯性重力波的影响.观测发现：在降冰雹前,每隔1~4小时出现一次短周期重力波阵性增强的现象.数值模拟结果显示：在低空降冰雹前几个小时有强的短周期重力波出现,其中周期较长的出现早、存在时间长,周期较短的出现晚、存在时间短;强的低空急流和风速垂直切变触发对流或湍流的发生和加强,对流或湍流又激发了80~200 min的短周期重力波;短周期重力波更容易向垂直方向传播,而长周期重力波倾向于水平方向传播.长周期重力波在降冰雹后周期有明显变短现象,随高度越加明显.由地形形成的重力波在最高山峰上空振幅最大.     

重力波诱导OH气辉分布的变化

利用全隐欧拉格式和全球热层-电离层-中间层-电动力学环流耦合模式（TIME GCM） ,本文分析了非线性重力波在向上传播过程中对OH气辉廓线的影响,结果表明,当较小垂直 尺度重力波经过OH气辉发射层时,对其峰值分布有明显的影响,有时会出现双峰现象,说明 非线性重力波传播是引起气辉双峰分布的一个直接原因;重力波在中层顶区域的破碎 在OH气辉分布均有不同程度的表现,这不仅可以解释一些实际的观测,同时也对利用OH气辉 观测分析重力波传播特性有明显的指示意义.      

Full-correlation analysis of turbulent scattering layers in the mesosphere observed by the MU radar

We have applied a full-correlation analysis technique to the echo power fluctuations observed by the MU radar (35°N, 136°E), and analyzed the horizontal structure of the scattering pattern in the mesosphere as well as their horizontal motions. The velocity of the scattering pattern did not agree with the background wind velocity, but was associated with the horizontal propagating direction of a saturated inertia gravity wave identified in the wind field. The length of the long axis of the characteristic ellipse of the scattering pattern was approximately 50 km, and the direction was almost perpendicular to the propagating direction of the wave. The correlation time of the scattering pattern was approximately 700 s, which is much longer than the lifetime of the isolated turbulence itself. This implies that the observed scattering pattern is associated with a region where the saturated inertia gravity wave generates turbulence.     

Cause of winter gravity wave spectrum saturation

Gravity waves play a significant role in establishing the large-scale circulation and structure of the middle atmosphere. Through gravity wave saturation proc-esses, such motions are believed to cause turbulence, resulting in divergence of momentum flux and the diffusion of heat and constituents in the meso-sphere[1,2]. The mechanisms that contribute signifi-cantly to the gravity wave saturation are thought to be the dynamical and convective instabilities[3]. However, it is difficult to distin…     

空间尘埃等离子体中的重力波特性

建立尘埃等离子体中重力波的基本方程，推导尘埃等离子体中重力波的色散关系，分析地球极区中间层顶处尘埃等离子体层中的重力波特性，研究了重力波在电子密度垂直分层的尘埃等离子体中的反射. 结果表明尘埃等离子体改变了通常大气中的重力内波的色散关系，限制了小水平波数重力内波的传播，改变了波的能量特性，减小了重力波在不均匀大气中垂直向上传播时振幅的增长；在尘埃等离子体中传播时重力波可被电子密度垂直分层的结构反射而导致波能量的集中, 它产生的湍动所导致的空间电子密度的不均匀性分布是极区上空PMSEs的可能机制.     

The effect of horizontal resolution on gravity waves simulated by the GFDL “SKYHI” general circulation model

To examine the effects of horizontal resolution on internal gravity waves simulated by the 40-level GFDL SKYHI general circulation model, a comparison is made between the 3° and 1° resolution models during late December. The stratospheric and mesospheric zonal flows in the winter and summer extratropical regions of the 1° model are much weaker and more realistic than the corresponding zonal flows of the 3° model. The weaker flows are consistent with the stronger Eliassen-Palm flux divergence (EPFD).The increase in the magnitude of the EPFD in the winter and summer extratropical mesospheres is due mostly to the increase in the gravity wave vertical momentum flux convergence (VMFC). In the summer extratropical mesosphere, the increase in the resolvable horizontal wavenumbers accounts for most of the increase in the gravity wave VMFC. In the winter extratropical mesosphere, the increase of VMFC associated with large-scale eastward moving components also accounts for part of the increase in the gravity wave VMFC.The gravity waves in the summer and winter mesosphere of the 1° model are associated with a broader frequency-spectral distribution, resulting in a more sporadic time-distribution of their VMFC. This broadening is due not only to the increase in resolvable horizontal wavenumbers but also occurs in the large-scale components owing to wave-wave interactions. It was found that the phase velocity and frequency of resolvable small-scale gravity waves are severely underestimated by finite difference approximations.     

用数字测高仪漂移测量研究电离层声重波扰动

数字测高仪Digisonde中的漂移测量,常用来研究小尺度电离层扰动,如电离层小不均匀体的漂移。文中提出一种新的分析方法,利用漂移测量数据中多普勒频移和到达角参量的最大熵动态功率谱,估算声重波一类大尺度电离层扰动的水平传播速度和传播方向。作为实例,研究了Millstone Hill测高仪站的漂移观测资料,并对处理结果进行了初步分析。分析结果表明,从数字测高仪漂移测量数据中,可有效地提取声重波一类大尺度电离层扰动的传播参量,在电离层动力过程的研究中很有意义。     

Simultaneous observations of planetary waves from 30 to 220km

Planetary wave activity at quasi 16-, 10- and 5-day periods has been compared at various altitudes through the middle and upper atmosphere over Halley (76°S, 27°W), Antarctica, during the austral winters of 1997–1999. Observational data from the mesosphere, E-region ionosphere and F-region ionosphere have been combined with stratospheric data from the ECMWF assimilative operational analysis. Fourier and wavelet techniques have shown that the relationship between planetary wave activity at different altitudes is complex and during the winter eastward wind regime does not conform to a simple combination of vertical planetary wave propagation and critical filtering. Strong planetary wave activity in the stratosphere can coincide with a complete lack of wave activity at higher altitudes; conversely, there are also times when planetary wave activity in the mesosphere, E-region or F-region has no apparent link to activity in the stratosphere. The latitudinal activity pattern of stratospheric data tentatively suggests that when the stratospheric signatures are intense over a wide range of latitudes, propagation of planetary waves into the mesosphere is less likely than when the stratospheric activity is more latitudinally restricted. It is possible that, on at least one occasion, 16-day planetary wave activity in the mesosphere may have been ducted to high latitudes from the lower latitude stratosphere. The most consistent feature is that planetary wave activity in the mesosphere is almost always anti-correlated to planetary wave activity in the E-region even though the two are in close physical proximity. The oscillatory critical filtering of vertical gravity wave propagation by planetary waves and the re-generation of the planetary wave component at higher altitudes through subsequent critical filtering or breaking of the gravity waves may provide an explanation for some of these characteristics. Alternatively the nonlinear interaction between planetary waves and tides, indicated in the E-region data, may play a role.     

基于OH全天空气辉成像仪观测的中国低纬地区的重力波传播统计特征

利用位于海南富克(19.5°N,109.1°E)和广西桂平(23.4°N,110.1°E)两个台站两年多的OH全天空气辉成像仪观测数据,对中国低纬地区的重力波传播统计特征进行了研究.从富克和桂平的气辉成像观测中, 分别提取了65和86个重力波事件.研究结果表明,观测水平波长,观测周期和水平相速度分别集中分布在10~35 km, 4~14 min和20~90 m·s^-1范围.重力波传播方向,在夏季表现出很强的东北方向传播.然而,在冬季主要沿东南和西南方向传播. 同时,结合流星雷达风场观测和TIMED/SABER卫星的温度数据,也发现在中层-低热层中传播的大多数重力波表现为耗散传播.且低层-中层大气中背景风场的滤波作用和多普勒频移可能对纬向方向传播的重力波产生的各向异性起到重要的调制作用.然而,经向方向传播的重力波产生的各向异性可能同时被低层大气中波源的非均匀分布以及潮汐变化所影响.     

Gravity wave signatures at mesopause heights

The ALOMAR SOUSY Radar operated at 53.5 MHz has been used in a five-beam configuration to study dynamical processes at gravity wave periods in the summer polar mesosphere. A case study of a gravity wave with a period of about 9 min is presented and analysed in some detail. The three-dimensional wave number vector is determined from the phase information of the 9-min velocity oscillations obtained in all beam directions and all range gates. The horizontal wavelengths in the north and east direction are estimated to be about 60 and 50 km, respectively. The echo power variations, simultaneously observed at different beam pointing positions, are investigated using cross-correlation analysis. The results show that these wave-associated variations lead to horizontal and vertical deformation of the echoing layers. The horizontal wavelength in the north-south direction, derived from the echo power modification, is in good agreement with the one obtained from the velocity analysis.     

Gravity wave reflection: Case study based on rocket data

Since gravity waves significantly influence the atmosphere by transporting energy and momentum, it is important to study their wave spectrum and their energy dissipation rates. Besides that, knowledge about gravity wave sources and the propagation of the generated waves is essential. Originating in the lower atmosphere, gravity waves can move upwards; when the background wind field is equal to their phase speed a so-called critical layer is reached. Their breakdown and deposition of energy and momentum is possible. Another mechanism which can take place at critical layers is gravity wave reflection.In this paper, gravity waves which were observed by foil chaff measurements during the DYANA (DYnamics Adapted Network for the Atmosphere) campaign in 1990 in Biscarrosse (44°N, 1°W)—as reported by Wüst and Bittner [2006. Non-linear wave–wave interaction: case studies based on rocket data and first application to satellite data. Journal of Atmospheric and Solar-Terrestrial Physics 68, 959–976]—are investigated to look for gravity wave reflection processes. Following nonlinear theory, energy dissipation rates according to Weinstock [1980. Energy dissipation rates of turbulence in the stable free atmosphere. Journal of the Atmospheric Sciences 38, 880–883] are calculated from foil chaff cloud and falling sphere data and compared with the critical layer heights. Enhanced energy dissipation rates are found at those altitudes where the waves’ phase speed matches the zonal background wind speeds. Indication of gravity wave trapping is found between two altitudes of around 95 and 86 km.