大气中的内船舶波

通过波动方程的分析解求得稳定层结大气中孤立地形激发的大气内船舶波,证明这类波是由非拦截波和拦截波两部分构成,非拦截波中的辐散波和横波的形态,王要由地形形状和走向决定,它们都可向大气上层传输,而拦截背风波的波形主要由层结、风速及风速切变等大气条件决定.     

《地球物理学报》总目次

0001磁暴期间几种主要磁扰成分的演化特征………吴迎燕,徐文耀,陈耿雄,陈博,…刘晓灿0010由非相干散射雷达数据重建极光沉降粒子能谱…………………………………蔡红涛,…马淑英0018电离层暴时的foF2预报技术研究…………………………陈丹俊,吴健,…王先义0024重力波非共振激发的数值研究…………………………黄开明,张绍东,…易帆0034三维分层流动过双山脉地形激发的大气船舶重力波动力学理论和数值试验………………………李子良0043夜间稳定边界层中小尺度地形激发的形式阻力和波动阻力……………姜瑜君,桑建国,刘辉志,…刘树华0051平流…     

风垂直切变对中尺度地形对流降水影响的研究

针对长江中下游中尺度地形特点以及暴雨过程发生发展期间风垂直切变的主要观测特征,设计了一系列中尺度地形的三维理想数值试验,分析了干大气地形流和重力波特征,探讨了条件不稳定湿大气地形对流降水的模态分布,在此基础上研究了圆形、直线风垂直切变和切变厚度对中尺度地形对流降水强度和模态分布的影响.结果发现:在 Fr≈1的干大气条件下,气流遇到地形后分支、绕流和爬升现象同时存在,地形激发的重力波在水平和垂直方向上传播,其在迎风坡、背风坡、地形上游和下游的振幅不同,并组织出不同强度的垂直上升运动.在Fr > 1的条件不稳定湿大气下,地形对流降水主要存在三种模态,即迎风坡和背风坡准静止对流降水以及地形下游移动性对流降水,地形对流降水的形成与重力波在低层组织的上升运动密切相关.风垂直切变对地形对流降水的强度和模态分布有重要作用,其中圆形风垂直切变(风随高度旋转)不仅影响地形下游对流降水系统的移动方向,而且影响迎风坡和背风坡山脚处对流降水中心的分布和强度;直线风垂直切变(风随高度无旋转)主要影响地形对流降水的移动速度和强度.风随高度自下而上顺(逆)时针旋转,地形对流系统向下游传播时向右(左)偏移.风垂直切变主要通过影响地形重力波的结构和传播以及对流系统的形成、移动方向和速度,来影响地形对流降水的模态分布,其中对流层中低层的风垂直切变对地形对流降水强度和模态分布有重要影响.     

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

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

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

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

利用激光雷达测量重力波三维结构

激光雷达观测得到的密度、温度等数据被广泛应用于大气重力波研究.瑞利激光雷达可以获取激光路径上的大气密度、温度数据.对于大气中的三维波动而言,单条路径上的观测参量能提取得到的波动信息有限.本文首先以单色重力波为例,分析了利用激光雷达直接观测三维波动结构的可行性.激光雷达垂直观测即可得到重力波的垂直波长,当激光雷达以一定的天顶角斜向测量时,所得到的波长包含了重力波的垂直波长以及水平波长信息.因此,利用激光雷达同时以三个方向（垂直、向南（天顶角30°）以及向西（天顶角30°））测量,可以提取得到重力波的垂直波长和水平波长.本文利用中国科学院国家空间科学中心研制的车载532 nm瑞利激光雷达的经向系统和纬向系统同时以不同的指向角观测大气重力波,对利用激光雷达获取三维波动结构的方法进行了分析研究.本文给出了北京地区激光雷达观测重力波的诸多案例,分析了30~60 km高度范围内北京地区大气重力波的垂直及水平波长信息.并以2017年11月7日观测的准单色重力波为例,结合再分析资料的风场数据,分析了该重力波的水平波长,垂直波长及传播方向等信息.

     

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

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

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

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

线性缓变下垫面和耗散共同作用的Rossby孤立波包

本文从正压流体的准地转位涡方程出发,考虑到Rossby波随时间、空间演变的多尺度特性,采用多重尺度和摄动展开的方法推导了在线性缓变下垫面和耗散共同作用下的非线性Rossby波包演变满足推广的非线性Schdinger方程,并得到单个包络的孤立波解.通过分析包络的孤立波解,指出在耗散与下垫面共同作用下,大气中会出现双曲正割形状的Rossby孤立波包.同时得到,在切变基本流存在的情况下,地形强迫和耗散对Rossby包络孤立波的传播速度、波数和频率均有影响的结论.     

热带海盆对热力强迫的线性响应

通过对线性两层海洋模式进行正交模求解,得到了热带矩形海盆在热力强迫下的海洋动力场水平结构.在这个线性两层模式中,没有施加风应力,仅考虑了热力强迫下的Rayleigh摩擦和Newton冷却效应.在一种理想化的经向不均匀加热强迫下,动力场表现出类似于风生环流的特征:窄而强的西边界区,宽而弱的东边界区;具有双涡(double-gyre)结构.线性响应中斜压模态比正压模态大一个量级,在响应中占主要地位.     

On the atmospheric internal ship waves

The analytical solutions of the atmospheric internal ship waves induced by three-dimensional terrain are obtained by solving the atmospheric wave equation. The solutions show that the waves consist of the untrapped and trapped parts. The patterns of the diverging wave and transverse wave in the untrapped parts are mainly determined by the shape and orientation of the terrain. This kind of wave may transport the wave energy to the upper atmosphere. The patterns of trapped lee waves are decided by the atmospheric conditions such as stratification, mean wind speeds and wind shear.     

On terrain corrections in terrestrial heat flow

The theory of heat-flow variation over a Lees-type hill is well known for evaluating terrain effects on terrestrial heat flow. However, a hill cannot be converted into a valley by simply changing the sign of one of the terrain-defining parameters, nor can several Lees' hills be superposed to simulate a series of hills and valleys for correcting a terrain-induced disturbance of heat flow. The proper derivations for a Lees-type valley are presented and the superposition is compared with analytic solutions for two parallel semi-cylindrical valleys and two semi-spherical basins embedded in an otherwise planar ground. Generally, heat flux climaxes over central valleys or basins, and dips toward their margins where rapid change in topography occurs. Variation of heat flux induced by a three-dimensional terrain is relatively large, as compared to two-dimensional features, but its areal extent is relatively limited. The two-dimensional effects also extend relatively deeper. In applying two-dimensional algorithms to a three-dimensional terrain, the correction may be over- or underestimated depending on the distance from a borehole to prominent terrain features in the surrounding area.     

Effects of three-dimensionality on seismic response of Gaussian-shaped hills for simple incident pulses

In this paper the boundary element algorithm which uses the time-convoluted traction kernels is applied to a numerical parametric study on the seismic behavior of three-dimensional Gaussian-shaped hills subjected to vertically propagating incident waves. All calculations were executed in the time-domain and the medium was assumed to have a linear elastic constitutive behavior. Results are discussed in both time and frequency domain with respect to the dimensionless parameters. It was shown that wave length and site geometry, including shape and dimension ratios and, to some extent, wave type are the key independent parameters governing hill amplification behavior. Comparing two- and three-dimensional hills with similar shape ratios, two-dimensional hill had greater characteristic periods, where the three-dimensional hill had greater maximum amplification potential. Three-dimensionality has a strong effect on the seismic responses of the hill; however the rate of seismic response variation with the three-dimensionality factor depends on the shape ratio. It was shown that two-dimensional behavior was dominant in low height three-dimensional hills, however, as the shape ratio increased, three-dimensionality effects appeared and the seismic response of the hill tends toward the axisymmetric three-dimensional hill.     

弹性层状半空间中凸起地形对入射平面SH波的放大作用

对Wolf理论进行拓展,使之可解决凸起地形对波的散射问题,进而利用间接边界元法,求解了弹性层状半空间中凸起地形对入射平面SH波的放大作用问题。本文模型的显著特点之一是考虑了层状半空间的动力特性以及层状半空间和凸起地形的阻尼;特点之二是计算精度高。文中以基岩上单一土层中半圆凸起地形对入射平面SH波的放大作用为例进行了数值计算分析。研究表明,基岩上单一土层中凸起地形对入射平面SH波放大作用和均匀半空间中凸起地形有着本质的差别;土层动力特性不仅影响凸起地形地表位移的幅值,还会影响地表位移的频谱;阻尼会显著降低凸起地形对高频波的放大作用。     

WIND TUNNEL EXPERIMENTS OF AEOLIAN DUST DEPOSITION ALONG RANGES OF HILLS

Wind tunnel experiments of aeolian dust deposition on topographic scale models of ranges of hills were conducted. Different hill sizes and hill spacings were used, and comparisons with the deposition patterns over single, isolated hills were made. Dust profiles over ranges of hills differ from the profiles over identical, but isolated hills. On isolated hills the sedimentation maximum on the windward hillslope is always single and located on the concave part of the slope. In the case of ranges of hills, the maximum is either single or double, with the second peak on the convex part of the windward slope in the latter case. The local sedimentation maximum on the convex leeslope, which is rather unimportant on isolated hills, is much more developed in multiple-hill topography. Also, dust deposition on the leeslopes is significantly higher in multiple-hill topography than on isolated hills. Dust patterns on ranges of hills may be affected by the dust shadow created by the most upstream-located hill. If hills succeed each other quickly, they are located within the shadow zone and are protected from important dust deposition. The plume of high air dust concentration that is created by a hill largely determines the dust pattern on the next hill. As a result of the supply of dust from above by the descending plume, areas that are normally devoid of dust now experience significant dust deposition.     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.

     

Observations indicating parametric instabilities in internal gravity waves at thermospheric heights

Abstract

Theoretical studies predict a parametric instability of finite-amplitude internal gravity waves which hitherto has been observed only in laboratory experiments. The occurrence of this process in the atmosphere is of basic interest because finite-amplitude gravity waves, which are almost ubiquitous especially at upper atmospheric heights, would produce unstable flows even at large Richardson numbers. Maximum entropy power spectra of a strong internal gravity wave in the thermosphere, which was generated by a volcanic eruption and detected on records of the Doppler shift of high-frequency radio waves, in fact show good agreement with the spectra of synthetic Doppler records obtained from a calculated unstable gravity wave. The frequencies and wavenumbers observed in the gravity wave domain satisfy in particular the theoretically predicted resonance conditions. The observed Doppler records also show two significant lines in the acoustic domain which probably result from a nonlinear interaction with the basic gravity wave. It is suggested that acoustic double peaks, which are commonly observed in high-frequency Doppler spectra in the presence of nearby thunderstorms, represent parametric instabilities of internal gravity waves generated by penetrative cumulus convection.     

On the GPS-based ionospheric perturbation after the Tohoku earthquake of March 11, 2011

Based on the data from the GPS receiving networks in Japan and America which have a high time resolution (2 min), two-dimensional (2D) distributions of the variations in the ionospheric total electron content (TEC) are constructed both close to and far from of the epicenter of the submarine earthquake of March 11, 2011 in Japan. Above the epicenter, a diverging multi-period disturbance appears after the main shock due to the acoustic gravity waves. Far from the epicenter, the wave trains associated with the tsunamigenic atmospheric internal gravity waves are revealed. These atmospheric waves significantly advance the arrival of the tsunami signal initially on the Hawaiian islands and then on the western coast of North America. The presence of the tsunami precursor in the form of atmospheric gravity waves is supported by the numerical calculations and by the analysis of the dispersion relation for the waves in the atmosphere. The detected ionospheric responses close and far from the epicenter can be used in the early tsunami warning systems.