A comparison of gravity wave phase velocities calculated from the impedance equation and pressure data

Abstract

Marked oscillations in wind speed, wind direction and pressure with periods of between 5 and 15 minutes were recorded on a number of occasions and it is assumed that they were produced by atmospheric gravity waves. Gossard and Munk (1954) have shown that the phase velocity of a gravity wave can be calculated using wind and pressure fluctuation data from a single station by means of the “impedance equation”. It is also possible to measure the phase velocity of a gravity wave using an array of microbarovariographs. In this investigation, the phase velocity of gravity waves calculated using these two different methods are compared and the feasibility of deriving wave characteristics from single point measurements is discussed. The results show that, while in a few isolated cases the two velocities agree reasonably well, the application of the impedance equation to data from a single station is limited due to the superposition of local wind systems upon the wave-induced perturbations.     

Oblique wave diffraction by a flexible floating structure in the presence of a submerged flexible structure

Expansion formulae associated with the interaction of oblique surface gravity waves with a floating flexible plate in the presence of a submerged horizontal flexible structure are derived using Green’s integral theorem in water of finite and infinite water depths. The associated Green’s functions are derived using the fundamental solution associated with the reduced wave equation. The integral forms of the Green’s functions and the velocity potentials are advantageous over the eigenfunction expansion method in situation when the roots of the dispersion relation coalesce. As an application of the expansion formulae, diffraction of oblique waves by a finite floating elastic plate in the presence of a submerged horizontal flexible membrane is investigated in water of finite depth. The accuracy of the numerical computation is demonstrated by analysing the convergence of the complex amplitude of the reflected waves and the energy relation. Effect of the submerged membrane on the diffraction of surface waves is studied by analysing the reflection and transmission coefficients for various parametric values. Further, the derivation of long wave equation under shallow water approximation is derived in a direct manner in the appendix. The concept and methodology can be easily extended to deal with acoustic wave interaction with flexible structures and related problems of mathematical physics and engineering.     

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

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

一次暴雨过程中重力波参数演变特征的模拟结果

本文利用中尺度数值模式WRF对2003年7月4-5日淮河流域特大暴雨过程进行了数值模拟,并利用高时空分辨率模拟结果资料,提取了暴雨中心区大气重力波频率、周期、水平波长、垂直波长、水平相速和群速等特征参数,分析了暴雨过程中重力波参数随时间的演变特征.结果表明,对此次暴雨强降水过程影响较大的重力波主要是发展的中 α 尺度波和中 β 尺度波,暴雨后期随着重力波的频散,周期和水平波长有减小趋势,频率有增大趋势.非降水区的重力波参数特征和降水中心区有明显不同,大气中小振幅的中 α 尺度和中 β 尺度重力波是否发展和暴雨强降水的发生关系密切.     

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.     

用传输函数构建的大气重力波传播理论模式

本文根据考虑大气热传导和黏滞的重力波复色散关系,采用传输函数的概念,基于重力波的线性理论,构建了用于研究对流层内重力波激发源与电离层响应之间的传输函数数值模式.在相空间中讨论了传输函数振幅的分布特性,并以地面单位脉冲源为例,分析了从地面到300 km高空的响应,得到了物理量的时空分布特征.结果表明：（1）对内重力波的传播而言,大气相当于一个滤波器,只有波动周期在15～30 min,水平波长在200～450 km之间的重力波扰动才最容易到达300 km电离层高度;（2）电离层的响应主要在与地面的激发源之间相隔较远的水平距离上发生;（3）黏滞和热传导系数在低层对上传重力波的影响较小,随着高度的增加它们对重力波的影响越来越大;（4）在低层计算的波动频率与Row理论的计算结果比较一致,然而到了高层却相差较大.     

Application of wavelet transformation to determine wavelengths and phase velocities of gravity waves observed by lidar measurements

During the last two decades, important advances have been made in the investigation of gravity waves. However, more efforts are needed to study certain aspects of gravity waves. In the real atmosphere, gravity waves occur with different properties at different altitudes and, most often, simultaneously. In this case, when there is more than one dominant wave, the determination of gravity wave characteristics, such as the vertical wavelength and the phase velocity, is difficult. The interpretation of temperature perturbation plots versus the altitude and time as well as the application of the Fourier spectral analysis can produce errors.Exact knowledge of the wave characteristics is important both for determination of other characteristics, for example, the horizontal wave components, and for study of wave climatology. The wavelet analysis of vertical temperature profiles allows one to examine the wave's location in space. Up to now, gravity waves have been studied mainly by continuous wavelet transformation to determine dominant waves. We apply wavelet analysis to a time series of temperature profiles, observed by the ALOMAR ozone lidar at Andoya, Norway, and by the U. Bonn lidar system at ESRANGE, Sweden, both for determination of the dominant waves and for specifying the vertical wavelengths and the vertical component of the phase velocities. For this purpose, the wavelet amplitude spectra and the wavelet phase spectra are filtered and Hovmöller diagrams for dominant wavelengths are constructed. The advantage of this type of diagrams is that they give clear evidence for the localization of the dominant waves in space and time and for the development of their phase fronts.     

Frequency shift of acoustic gravity waves in a stratified,isothermal, turbulent atmosphere

We study the modification of the frequency of small amplitude acoustic gravity waves which propagate in an isothermal turbulent atmosphere that is stratified by a homogeneous gravitational field. Using a Green's function method, the dispersion relation for the frequency of the waves is formulated as an integral eigenvalue equation and it is solved by perturbation techniques. We draw the following main conclusions: (a) for an arbitrary turbulent correlation spectrum the dispersion relation has a root with a negative imaginary part in the unphysical Riemann sheet of the dispersion function, leading to wave attenuation, much in the same way as it happens for Landau damping; (b) the real part of this root differs from the frequency of an acoustic gravity wave propagating in a nonturbulent medium and, for all forms of the turbulent correlation spectrum, the absolute value of this difference increases if gravity increases; (c) for a Gaussian turbulent correlation spectrum, this difference is always positive; (d) conversely if this frequency difference is known from observations, the auto-correlation function of the temperature fluctuations can be calculated through a simple inversion formula.     

Characteristics of disturbances in the atmosphere and oceans

The characteristics of the disturbances in the atmosphere and oceans and in other stably stratified and rotating fluids are analyzed according to their phase and group velocities. It is shown that both stable stratification and rotation augment the velocity of the sound waves, and that the internal gravity waves and inertial waves are mutually exclusive when the Brunt-Väisälä frequency is different from the Coriolis parameter. It is also shown that both the barotropic and the internal Rossby waves are well separated from the gravity waves and that they can be represented accurately by the quasi-geostrophic potential vorticity equation, even close to the equator, except for the one member withn=0 which is coupled with an eastward propagating gravity wave.     

论不同模式大气中的静力适应过程

为了研究大气静力平衡适应过程的本质,利用波动理论和能量转换角度,分别对完全可压缩的等温大气模型、滞弹近似下的等温大气模型和层结中性大气模型进行研究比较.结果表明：大气静力平衡适应过程的本质是声波和混合声重力波对扰动能量的频散过程,滞弹近似模型和层结中性模型均不能完全描述此适应过程;在波动假设下,此三类大气模型中扰动物理量之间的偏振关系同波动的性质有关,气团的运动方程均为椭圆方程,声波和混合声重力波对气团运动的作用差异较显著.
大气静力平衡适应过程中扰动能量以有效势能、有效弹性势能、动能或波动能量的形式存在并相互转换;扰动有效势能与其他形式能量之间的转换与混合声重力波或者重力内波有关,扰动有效弹性势能与其他形式能量之间的转换与声波有关.在完全可压缩的等温大气模型中,扰动有效势能增加1个单位,其中69.9%来自扰动垂直动能,其余30.1%来自扰动有效弹性势能.     

利用瑞利激光雷达观测北京地区上平流层地形重力波活动

本文利用中国科学院空间科学与应用研究中心的瑞利激光雷达首次观测到了平流层地形重力波活动的现象,并结合美国国家环境预报中心(NCEP)的全球预报系统(GFS)的风场数据分析了该地形重力波的基本参数.与惯性重力波相比较,地形重力波的密度扰动没有下传的相位,在同一高度上,其扰动相位保持不变.北京空间科学与应用研究中心瑞利激光雷达自2012年开始观测实验以来,已经观测到多起地形重力波活动事件.本文以2013年11月11日的观测数据为例,研究北京上空的地形重力波活动,并结合GFS风场数据分析了北京上平流层地形重力波的波长、传播方向、传播速度等参量.通过分析得到在2013年11月11日北京上空存在一列传播方向为北偏西52.4°,水平波长为5.5km,平均垂直波长约为6.0km的地形重力波.     

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

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

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

Wind control of the propagation of acoustic gravity waves in the polar atmosphere

The relationship between the directions of polar acoustic gravity waves and a wind at 250–350 km altitudes has been studied based on an analysis of the Dynamics Explorer 2 satellite measurements. A method, which makes it possible to determine the direction of these waves relative to the satellite velocity vector based on one-point measurements of different neutral atmosphere parameters, is presented. It has been established that acoustic gravity waves observed over the polar caps systematically propagate upwind, which argues for their spatial wind filtering. In the polar regions, waves mainly propagate in two directions: toward magnetic noon and 15–16 MLT. Waves tend to move counterclockwise and clockwise over the northern and southern polar caps, respectively.     

台风“梅花”诱发平流层重力波的数值模拟与AIRS观测

为了分析台风这类强对流诱发平流层重力波的过程,本文利用中尺度数值模式WRF-ARW(V3.5)和卫星高光谱红外大气探测器AIRS数据对2011年第9号强热带气旋"梅花"的重力波特征进行了分析.首先,针对模式输出的垂直速度场资料的分析表明,台风在对流层各个方向上几乎都具有诱发重力波的能量,而在平流层内则呈现出只集中于台风中心以东的半圆弧状波动,且重力波到达平流层后其影响的水平范围可达1000km.此外,平流层波动与对流层雨带在形态、位置以及尺度上均具有一定的相似性.其次,对风场的分析结果表明,不同高度上波动形态的差异主要是由于重力波垂直上传的过程中受到了平流层向西传的背景风场以及风切变的调制作用,揭示了重力波逆着背景流垂直上传的特征.随后,基于FFT波谱分析的结果表明,"梅花"诱发的平流层重力波水平波长中心值达到了1000km,周期在15~25h,垂直波长主要在8~12km.最后,利用AIRS观测资料分析了平流层30~40km高度上的大气波动,得到了与数值模拟结果相一致的半圆弧状波动.对比结果也验证了WRF对台风诱发平流层重力波的波动形态、传播方向、不同时刻扰动强度的变化以及影响范围的模拟效果.此外,也揭示了多资料的结合对比有助于更加全面地了解台风诱发平流层重力波的波动特征.     

Ionospheric Response to the Acoustic Gravity Wave Singularity

An original model of atmospheric wave propagation from ground sources to the ionosphere in the atmosphere with a realistic high-altitude temperature profile is analyzed. Shaping of a narrow domain with elevated pressure in the resonance region where the horizontal phase wave velocity is equal to the sound velocity is examined theoretically within the framework of linearized Eq.s. Numerical simulations for the model profiles of atmospheric temperature and viscosity confirm analytical result for the special feature of wave fields. The formation of the narrow domain with plasma irregularities in the D and low E ionospheric layers caused by the acoustic gravity wave singularity is discussed.     

径向分层TI孔隙介质井孔中激发的模式波的数值研究

横向各向同性介质是地层中普遍存在的一种各向异性介质.本文对径向分层TI孔隙介质包围井孔中激发的斯通利波和弯曲波的传播特性进行了理论计算,发现模式波在低频时更多的是反应原状地层的信息,而随着频率的增加侵入带参数逐渐起控制作用;Biot理论描述的地层衰减比速度更容易受井壁附近地层参数的影响.利用灵敏度曲线定量研究了不同频率下地层各个参数对相速度和衰减系数的贡献大小,主要结果显示模式波的衰减受水平渗透率影响明显,而垂直渗透率的变化对模式波几乎无影响;斯通利波对水平向传播的横波速度比弯曲波的灵敏度高.从单极子和偶极子声源在井孔中激发的全波波形也可发现,声波测井仪器较宽的声源频带和合适的源距设置有利于对不同径向深度上的地层声学参数进行成像.     

Electromagnetic and plasma effects during a carrier rocket flight

The disturbance generation model for the total electron content of the ionosphere and formation of the narrowband spectrum of electromagnetic disturbance on the Earth during a rocket flight along the horizontal leg of the trajectory has been considered. It has been indicated that a change in the total electron content is caused by the propagation of an acoustic gravity wave pulse, generated during a rocket flight along the horizontal trajectory leg, in the ionosphere. This pulse forms horizontal inhomogeneities of ionospheric conductivity in the bottomside ionosphere. Electric currents, induced by the background electromagnetic field in these inhomogeneities, are emitters of discrete modes of coherent gyrotropic waves propagating horizontally in a conductive layer of a finite thickness in the bottomside ionosphere. The line spectrum of electromagnetic disturbances has been calculated. The calculation results agree with the observational data.     

Simulation of improved pure P-wave equation in transversely isotropic media with a horizontal symmetry axis

Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation.     

非对称圆弧片状声源在井旁地层中产生的声场

提出了一种可用于随钻三维反射声波成像测井的圆弧片状声源,推导了该声源声学性质在波数-频率域内的数学描述,并利用实轴积分的方法对该声源在井旁地层中产生的声场进行求解,考察了声源的线度、频率等因素对该声源向地层中辐射的声场的影响.研究结果表明,特定尺寸和特定频率下的圆弧片状声源产生的纵波场的水平和垂直指向性图中仅具有一个明显的主瓣,且主瓣三分贝角宽较窄,方位分辨率较高,指向性良好,适用于随钻三维反射成像测井;SV波场和SH波场的指向性图中存在两个或者多个角瓣,SV波场在垂直于井轴方向上辐射的能量为0,另外在井中发射和接收的反射SH波场和SV波场会相互干扰,所以难以利用SV波场和SH波场进行三维反射声波成像;声源的线度和频率等因素对其辐射声场影响较大,考虑到激励效率和方位分辨率等因素,对于本文描述的井孔模型,选择圆周角在75°和90°之间, 主频为12 kHz左右的声源是合适的.