Aura/MLS与TIMED/SABER观测全球重力波特性

大气重力波是临近空间环境主要大气波动之一,对全球环流具有重要影响。卫星上搭载的临边探测器能够探测临近空间大气温度,可用于临近空间大气重力波研究。利用2012－2014年Aura的微波临边探测器（MLS）和TIMED的红外临边探测器（SABER）的探测数据,对20~50 km高度的大气重力波扰动分布特征开展了分析研究,两种观测重力波活动基本一致,重力波随季节、纬度及高度的变化显著。冬季半球高纬度重力波扰动较强,赤道和夏季半球近赤道地区上空也存在明显重力波活动区域,夏季半球高纬度重力波扰动最弱。重力波扰动强度随高度增加。TIMED/SABER重力波扰动强度数值比 Aura/MLS略强。      

Estimates of mesospheric gravity wave activity over convection from a global model

A new convective gravity wave source spectrum parameterization has been implemented in the Whole Atmosphere Community Climate Model version 2 (WACCM2). This parameterization specifies the momentum flux phase speed spectrum of gravity waves in the Tropics based on the properties of underlying convection; Hence, this parameterization provides realistic global estimates of gravity wave activity. In this paper, we show the estimated gravity wave phase speed spectra in the Tropics from a WACCM2 simulation, at the source level and at 85 km. Spatial distribution of gravity wave activity at 85 km is also presented. Subsequently, we discuss the factors that are primarily responsible for the estimated differences in gravity wave distribution across phase speeds with latitude and asymmetries in direction of gravity wave propagation in the mesosphere. We also examine which of the model assumptions can lead to uncertainties in our estimates of mesospheric gravity wave activity and we discuss how these assumptions provide challenges for comparison with observations of gravity waves in the mesosphere.     

The response of the high-latitude thermosphere to geomagnetic substorms

Enhanced ionization and electric fields associated with geomagnetic substorms lead to time-varying momentum (ion drag) and heat (Joule dissipation) sources in the high latitude thermosphere. These momentum and heat sources, in turn, cause changes in the structure of the wind, pressure and temperature fields, and, in addition, provide a mechanism for the generation of atmospheric gravity waves. We have developed a model which has been used to simulate, on a spatial scale of tens of kilometers, the response of the high latitude thermosphere to these momentum and heat sources, which constitute a major factor in the coupling between the magnetosphere, the ionosphere and the upper atmosphere. Simulations reported herein indicate that, because the coupling between ions and neutrals is highly dependent on the neutral-ion collision frequency, the response of the neutral atmosphere varies, not only as a function of total ionization, but also as a function of the details of the electron density profile.     

Seasonal variation of gravity wave sources from satellite observation

Recent high-resolution satellite observations of gravity waves in the middle atmosphere have shown correlations with the strength of the stratospheric jet stream, surface topography, and tropical convection. Seasonal variations of wave-induced stratospheric radiance variances are often the manifestations of modulations of these sources and refractive influences. In this paper, we focus on the seasonal climatology of gravity waves observed by the UARS MLS, while also showing some new results from GPS and AMSU instruments. Our analysis is aided by MWFM modeling of mountain waves at high latitudes and CMAP precipitation indices in the subtropics to provide a clearer picture of global gravity wave dynamics.     

对流层雷暴源所产生的柱形重力波

已有充分的证据表明, 大气对流层的雷暴迹象是大气中间层重力波活动的显著代表源.在雷暴迹象的上方, 通过火箭已观测到大气中间层出现的热效应, 也已通过雷 达探测到大气平流层出现的上行重力波, 从地面和卫星平台上观察到夜间气辉有序而 成环状的重力波波形. 所有这些实验结果都与位于观测点下方的雷暴活动有紧密联系. 此类雷暴通常主要集中在中国东部沿海以及地球其他沿海海湾地域. 关于此类雷暴对大 气中间层的影响尚未被充分研究和了解. 为能有效地探究其成因, 利用所开发的一个二维计算机数值模型模拟和研究大气对流层的雷暴源所引发的上行重力波, 进而揭示 此类重力波产生的基本物理机理, 及其在大气中间层的能量耗散. 通过模拟研究发现, 雷暴源可以大面积高强度地聚集和释放积雨云的能量, 当这种周期性的对流变化引发大气对流层的不稳定性后, 就会有圆柱体重力波的产生和传播.      

Numerical modeling of inhomogeneous orographic wave influence on planetary waves in the middle atmosphere

Parameterization of dynamical and thermal effects of stationary orographic gravity waves (OGWs) generated by the Earth’s surface topography is incorporated into a numerical model of general circulation of the middle and upper atmosphere. Responses of atmospheric general circulation and characteristics of planetary waves at altitudes from the troposphere up to the thermosphere to the effects of OGWs propagating from the earth surface are studied. Changes in atmospheric circulation and amplitudes of planetary waves due to variations of OGW generation and propagation in different seasons are considered. It is shown that during solstices the main OGW dynamical and heat effects occur in the middle atmosphere of winter hemispheres, where changes in planetary wave amplitudes due to OGWs may reach up to 50%. During equinoxes OGW effects are distributed more homogeneously between northern and southern hemispheres.     

Uncertainties in gravity wave parameters,momentum fluxes,and flux divergences estimated from multi-layer measurements of mesospheric nightglow layers

Measurements of dynamic parameters of atmospheric gravity waves, mainly the vertical wavelength, the momentum flux and the momentum flux divergence, are affected by large uncertainties crudely documented in the scientific literature. By using methods of error analysis, we have quantified these uncertainties for frequently observed temporal and spatial wave scales. The results show uncertainties of ～10%, ～35%, and ～65%, at least, in the vertical wavelength, momentum flux, and flux divergence, respectively. The large uncertainties in the momentum flux and flux divergence are dominated by uncertainties in the Brunt-Väisälä frequency and in spatial separation of the nightglow layers, respectively. The measured uncertainties in fundamental wave parameters such as the wave amplitude, intrinsic period, horizontal wavelength, and wave orientation are ～10% or less and estimated directly from our nightglow image data set. Other key environmental quantities such as the scale height and the Brunt-Väisälä frequency, frequently considered as constants in gravity wave parameter estimations schemes, are actually quite variable, presenting uncertainties of ～4% and ～9%, respectively, according to the several solar activity and seasonal atmosphere scenarios from the NRLMSISE-00 model simulated here.     

大气波动监测仪在岩石圈-大气层-电离层研究中的应用

针对北京大学空间物理与应用技术研究所研制的两台大气波动监测仪近三年的观测数据,对这一时段内所观测到的重力波和次声波周期尺度扰动的形态特征、谱结构特征及其在时间分布上的统计特征进行了分析. 给出了几例雷暴、地震等事件中观测到的大气扰动,并揭示了这些事件期间观测到的与地面大气扰动周期尺度相类似的电离层扰动. 结果表明,北京大学研制的大气波动监测仪可有效记录到地面的微弱大气扰动,观测数据可用于进行岩石圈-大气层-电离层之间的耦合研究.      

冕洞网络的大气模型

本文采用非径向磁流管位形的假设,计算了太阳冕洞网络部分的色球-日冕过渡区的能量平衡模型。所考虑的能量流包括辐射、传导、对流和机械波加热（如阿尔芬波）,计算结果表明在冕洞网络部分的过渡区中,电子温度T和密度N分别比宁静太阳中的值低60%和2倍,而其过渡区的厚度比宁静太阳中的大4倍。这种大气模型可满意地解释T≥105K范围的远紫外观测发射量度的分布。另外我们也发现在冕洞大气的过渡区中,阿尔芬波加热似乎不能忽略,尤其是在冕洞的过渡区底部,它的加热作用可能会超过热传导。在冕洞大气中,由于波动量的淀积而产生的对流能损耗也是重要的,在过渡区底部650km以上,对流能损耗逐渐超过辐射损耗。      

地形产生的山地波及其传播过程模拟研究

气流经过地形产生的山地波是大气重力波研究的重要类别之一.从大气运动的控制方程组出发,建立模拟地形产生的山地波及其传播过程的二维数值模式.利用水平背景风场、地形和垂直速度之间的关系,在模式中引入垂直速度扰动作为地形产生山地波的激发源.通过模拟该激发源引起的扰动即山地波在大气中的传播过程,再现了山地波的产生、传播及充分发展过程.通过分析水平波长、垂直波长、位温扰动、流线,在空间尺度上描述了山地波的产生、传播及充分发展的过程.在气流经过地形产生的山地波的传播过程中,其水平波长λ_x的范围为2.5~5km,垂直波长λ_z约为2.5km.这些结果与利用山地波线性理论计算的垂直波长一致,从而验证了本模式能够模拟地形产生的山地波及其传播过程,为深入了解山地波的产生过程及其在中高层大气中的传播机制和效应奠定了基础.      

卫星观测台风重力波数值模拟与直接对比验证

以ECMWF的T799资料作为WRF-ARW（V3.5）初始场,对卫星高光谱红外大气垂直探测器AIRS观测的2013年超强台风苏力激发平流层重力波过程进行数值模拟,并利用卫星观测对数值模拟结果进行了直接对比验证.数值模拟表明,该台风诱发的重力波在20~40km高度逆着东风背景流向东向上传播,在水平方向呈半圆弧状;大气的垂直扰动随着高度的增加而增强,在40km高度上达到0.5m·s-1.基于三维傅里叶变换的波谱分析表明,平流层重力波水平波长中心值在500km附近,周期为3~5h,垂直波长主要为10~26km.分析表明,在18~40km高度的净纬向动量通量为6.7×10-4~1.89×10-3Pa,背景流强迫计算值为-0.23~1.21m·s-1·d-1,且在18km和40km高度的数值较大.最后,基于辐射传输模式计算的直接对比表明,卫星观测与数值模拟同时揭示了激发的平流层波动可传至40km以上高度及距台风中心2000km以外的区域,且不同资料得到的波动形态、方位以及水平尺度具有较好的一致性.      

一种全天空气辉成像仪数据预处理方法和初步结果

大气重力波是大气中的基本波动形式之一,在中高层大气动力和热力学过程中起着十分重要的作用.全天空气辉成像仪是一种以大气气辉辐射为示踪物,能够有效对大气重力波成像的仪器.本文针对中国科学院国家空间科学中心空间天气学国家重点实验室中高层大气组自主研制的全天空气辉成像仪所观测的数据,提出了一种气辉图像预处理方法,进行平场校正、方位校正、星光去除和坐标映射等数据订正.利用该方法处理2015年5月17日21:00BLT至次日05:00BLT西宁台站（36.6°N,101.7°E）的OH气辉数据,发现一次重力波事件,分析并获得了该重力波的水平波长、观测水平相速度和传播方向（分别为17.72km,47m·s-1,339°）.研究结果表明该方法是可行的.      

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.     

损耗大气中惯性重力波的非线性相互作用方程及其初步讨论

根据弱相互作用理论,本文建立了损耗大气中极性重力波的非线性相互作用方程。这组方程在三个方面推广了前人的工作:考虑了波的空间传播;包含了粘滞产生的衰减;波谱可以是连续的。粘滞衰减率的大小与波的空间尺度以及传播方向有关。Coriolis力的引入使相互作用系数成为复量。根据这组方程,考察了惯性重力波的参量激发。结果表明:在共振条件满足时,主波存在一个阈值,阈值大小与次级波的损耗率成正比。当主波振幅大于这个阈值时,次级波将指数增长。在相互作用过程中,次级波的频率将发生变化,变化的大小与主波能量成正比。      

Evaluation of uncertainty in gravity wave potential energy calculations through GPS radio occultation measurements

The application of the Global Positioning System (GPS) radio occultation (RO) method to the atmosphere enables the determination of height profiles of temperature, among other variables. From these measurements, gravity wave activity is usually quantified by calculating the potential energy through the integration of the ratio of perturbation and background temperatures between two given altitudes in each profile. The uncertainty in the estimation of wave activity depends on the systematic biases and random errors of the measured temperature, but also on additional factors like the selected vertical integration layer and the separation method between background and perturbation temperatures. In this study, the contributions of different parameters and variables to the uncertainty in the calculation of gravity wave potential energy in the lower stratosphere are investigated and quantified. In particular, a Monte Carlo method is used to evaluate the uncertainty that results from different GPS RO temperature error distributions. In addition, our analysis shows that RO data above 30 km height becomes dubious for gravity waves potential energy calculations.     

AIRS观测资料研究全球平流层重力波特性

利用2012—2014年1月和7月AIRS（Atmospheric Infrared Sounder）第79通道的观测数据,分析了平流层重力波活动强弱的全球分布以及重力波发生频率的全球分布;分析了重力波活动随纬度和经度的变化特征,给出了重力波活动在全球范围内的热点区域及其活动强度;对比了白天与夜间的重力波活动强度及发生频率.研究表明重力波活动强度呈现出随纬度变化的特征,在低纬度地区（0°—30°）,冬季半球重力波活动强度低,夏季半球重力波活动高;在中高纬度地区,冬季半球重力波活动强度高而夏季半球重力波活动强度低.在1月,全球重力波活动有4个突出的热点区域,分别为50°N附近欧洲大陆与大西洋交接地带、北美洲与大西洋交接地,20°S附近南美洲与大西洋交接地区、非洲与印度洋交接地区.在7月,重力波活动突出的地方为巴塔哥尼亚至南极半岛地区,以及50°S和75°E附近的印度洋区域.重力波活动强度在夜间大于白天,但是夜间的强重力波活动区域小于白天.      

内重力波传播的3维传输函数模式研究

在考虑背景风场及大气耗散的条件下,建立了3维内重力波传输函数数值模式.分析了300 km高度3维传输函数在频率波数域的特性,并以近地面单位脉冲点源为激发源,得到了内重力波在3维空间中的时空分布.讨论了不同时空尺度地面方波源激发的内重力波在电离层高度的能量分布特征.结果表明,(1)对内重力波而言,背景大气相当于一个带通滤波器,只有波动周期和波长分别在15～30 min和200～400 km之间的重力波扰动最容易上传到300km高度;(2)在背景风场的作用下水平面上以同心圆扩散的波阵面以及垂直方向上成漏斗状的波阵面发生了变形,并且逆风方向比顺风方向更有利于声重力波由对流层向电离层高度传播;(3)300km高度对时间尺度和空间尺度分别在20～30 min和150～250 km之间的地面方波源响应的总能量最强.      

粘性耗散对重力波波包共振相互作用的影响

采用数值方法考察了粘性耗散对重力波波包共振相互作用的影响．结果表明,粘性耗散不仅衰减相互作用波的能量,而且强烈地影响波包之间局地能量的传递,使波包的形状发生进一步的改变．粘性耗散减缓相互作用过程的发展,阻止高波数的次级波与主波进行相互作用．在低热层高度上,耗散过程占支配地位,三波能量演变的特征发生了相当大的变化．大尺度主波通过与小尺度次级波的相互作用,更快地实现对大气的加热．      

武汉上空惯性重力波活动

利用武汉中心气象台(30.5°N,114.4°E)无线电探空仪在2001年1月到2003年12月之间的观测数据,研究了武汉上空1-9km和14-25km高度范围内惯性重力波参数的变化特征.数据分析结果表明,重力波活动存在明显的季节变化,冬季较频繁,夏季活动较弱,与急流强度的季节变化一致,这意味着重力波的激发与背景风密切相关.矢端曲线分析显示1-9km范围内能量向上传播和向下传播的波的比例大致相同,而14-25km范围内绝大部分波能量向上传播,这与最大急流强度的高度相对应;进一步分析表明,夏天14-25km范围内波几乎全向上传播,而冬天则有相当一部分波向下传播,这可能是反射引起的.重力波的本征频率集中在1-3倍的惯性频率之间,水平波长约数百公里;1-9 km范围内垂直波长集中在3-3.5 km,14-25 km高度内则集中在4.5-5.5 km.      

重力波波包在可压大气中的非线性传播

本文采用二维全隐欧拉(FICE)格式对具有高斯分布的重力波波包在等温、可压大气中的非线性传播过程进行数值模拟和分析.数值分析结果表明:尽管存在非线性效应,在整个传播过程中,波动的等相面向下运动,波包和波相关能量向上传输.波相关扰动速度随高度增加指数增长,并且波与平流会发生非线性相互作用,最后导致平均流场增强.这与线性重力波理论完全一致.重力波波包的传播路径与重力波线性射线理论预言非常接近,但平均水平群速度和平均垂直群速度均明显小于线性射线理论给出的结果,可见波动的非线性过程会改变波相关能量的传输速度.模拟结果首次定量地展示出非线性效应对重力波波包传播的影响,表明建立在线性理论框架中的重力波运动学定义的合理性.