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

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

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

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

基于MST雷达数据的中纬对流层和低平流层大气行星波研究

通过分析中国河北香河站MST (Mesosphere-Stratosphere-Troposphere)雷达 2012-2014年的水平风场数据, 研究了北半球中纬地区对流层和低平流层 (Troposphere and Lower Stratosphere, TLS)区域大气行星波的特性. 谱分 析发现, 在这一区域准16天波和准10天波占据主导地位, 准16天波更为显著. 在 对流层区域, 行星波具有丰富的频谱成分, 活动具有间断性, 持续时间一般不 超过三个月, 并没有明显的季节性变化特征, 其中纬向分量的振幅大于经向分量. 在 平流层区域(高度17km以上), 行星波一般出现在冬季, 并且主要在纬向分量中. 通常平流层区域的振幅要小于对流层区域. 结合MERRA再分析资料分 析了强行星波传播特性, 结果表明: 2014年2-3月纬向分量中的准16天波垂 直向上传播, 垂直波长约为64km, 纬圈波数约为2, 纬向传播方向自西向东, 水平波长约为15324.7km, 对应的相速度为11.1m·s-1 (向东为正); 2014年5月纬向分量中的准10天波在10~18km高度范围内向下传播, 垂直波长约为50km, 纬圈波数约为1, 传播方向自西向东, 水平波长约为 30649.4km, 对应相速为35.5m·s-1.      

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

以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以外的区域,且不同资料得到的波动形态、方位以及水平尺度具有较好的一致性.      

星载全球大气波动成像仪研究

在卫星上对中间层顶区域的气辉辐射成像观测, 对于全球大气波动的监测和研究具有重要的意义. 利用TDICCD对O₂A (0-0)气辉进行成像观测, 计算了曝光积分时间、信噪比和能达到的最高空间分辨率, 分析了地球自转对空间分辨率的影响, 在此基础上提出一种星载全球大气波动成像仪方案, 该方案可以观测垂直波长大于10km的大气波动, 最高水平分辨率可以达到0.33km.      

冬季中层顶区域重力波谱特性

使用DYANA试验期间19枚Chaff火箭在81—105km高度范围测量的高分辨率水平速度数据，研究标量水平速度垂直波数谱在冬季的谱特性．分析结果显示冬季平均垂直波数谱有-3．02的谱斜率，这与各种饱和理论预测的-3谱斜率很一致，提供了饱和重力波谱在冬季的观测证据．与各种饱和理论比较表明，这个冬季谱的谱振幅与线性饱和理论预测的谱振幅完全一致，但是总是小于非线性饱和理论预期的、激光雷达测量的谱振幅．这个冬季谱与最近得到的夏季饱和重力波谱一道，共同表明我们的观测结果支持了线性饱和理论，而不管季节，地理位置和高度范围的不同．由这个冬季谱估计的盛行垂直波长为15km，它比夏季的12．8km盛行垂直波长稍大．     

极区中层顶区域重力波谱的季节变化

使用MAP/WINE和MAC/SINE两次试验中测量的25m高分辨率水平速度数据和1km低分辨率温度数据,研究极区中层顶区域重力波谱的季节变化.温度的直接测量使计算的谱振幅和Richardson数更接近真实大气.结果显示,极区中层顶区域水平速度垂直波数谱的斜率和振幅存在相当大的变率,这些大的观测变率用各种饱和模式及普适垂直波数谱不能解释.然而平均垂直波数谱显示了明显的季节变化,在夏季,平均谱具有饱和特性;在冬季,平均谱具有非饱和特性.这意味着饱和过程存在于夏季而不是冬季.因此,夏季比冬季应有更强的湍流.这个结果与湍流季节变化的观测大致一致.从Brunt-Vaisala频率N和水平风切变计算的Richardson数Ri剖面也显示出季节差异,Ri<1/4的动力不稳定区出现在夏季,而Ri>0.4的稳定区出现在冬季.这些不稳定区与夏季谱结合很好,而稳定区则与冬季谱结合很好.     

80-100km高度范围高分辨率大气水平速度扰动谱

使用1989/1990年冬季在DYANA试验中从5枚Chaff火箭测量获得的高分辨率水平速度数据,讨论80-100km高度范围内水平速度垂直波数谱与饱和重力波谱的一致性,并决定特征垂直尺度.      

中国中部低层大气行星波无线电探空仪的观测研究

通过分析武汉、宜昌和恩施气象局无线电探空仪2001-2003年的观测数据,研究了中国中部地区对流层和低平流层中行星波的特性.通过Lomb-Scargle(L-S)的周期图方法发现了周期为准16天和周期为准10天的谱分量占据着主导地位.观察发现,较大振幅的行星波振荡主要集中在5-15 km之间.准16天行星波沿纬圈向西传播,对应的纬圈波数大约为2,水平波长约为17 324.8 km,传播相速度约为-12.5 m·s-1(东向为正),通过计算准16天行星波在10 km以下相位随高度的改变可以得到其垂直波长大约为25-30 km,而在对波层顶附近其相位几乎没有发生改变,呈现出静态波特性.准10天行星波沿纬圈向东传播,对应的纬圈波数大约为4,水平波长约为8627.3 km,传播相速度约为10.0 m·s-1,垂直波长约为22-40 km.      

中纬度冬季低热层潮汐非线性相互作用的MF雷达观测

采用武汉(30°N,114°E)MF雷达在2001年冬季的风场观测数据研究中纬度低热层大气潮汐之间的二阶非线性相互作用.经向风场的Lomb-Scargle归一化振幅谱表明,周日、半日和8 h潮汐是中纬冬季中层顶区域占优势的大气扰动;此外6 h潮汐也清晰可见.双相干谱分析揭示大多数显著的双相干谱峰代表潮汐谐振分量之间的相位互相关或单个潮汐分量的自相关.对随时间变化的潮汐垂直波长的比较发现,实际观测的8h潮汐垂直波长与假定的由观测的24 h潮汐和12 h潮汐非线性相互作用产生的8 h潮汐的理论垂直波长具有明显的一致性.在94.0～98.0km高度范围,周日、半日和8h潮汐之间不仅存在明显的相位相关和垂直波数相关,且它们的振幅随时间变化也显示出振荡幅值相近、振荡相位同步或反相的相关性,表明它们之间已经发生了二阶非线性相互作用.但是在94.0 km以下,三个潮汐分量之间的各种相关性随高度的下降变得越来越弱,因此潮汐二阶相互作用更可能是一种局地和暂态的现象.      

Planetary and gravity waves in the mesosphere and lower thermosphere

Planetary and gravity waves contribute significantly to the variability of atmospheric parameters in the middle atmosphere. In the mesosphere and lower thermosphere the wave fluctuations are sufficiently large to often mask the prevailing or mean state of the atmosphere. This review summarizes current knowledge about the motion, temperature and density fields associated with both large and small scale waves and stresses improved understanding that has come from recent ground based, rocket and satellite investigations.     

Planetary and gravity waves in the mesosphere and lower thermosphere

Planetary and gravity waves contribute significantly to the variability of atmospheric parameters in the middle atmosphere. In the mesosphere and lower thermosphere the wave fluctuations are sufficiently large to often mask the prevailing or mean state of the atmosphere. This review summarizes current knowledge about the motion, temperature and density fields associated with both large and small scale waves and stresses improved understanding that has come from recent ground based, rocket and satellite investigations.     

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.     

Directional Filtering Due to Mesospheric Wind Shear on the Propagation of Acoustic-gravity Waves

Gravity waves with periods close to the Brunt-V?is?l? period of the upper troposphere are often observed at mesopause altitudes as short period, quasi-monochromatic waves. The assumption that these short period waves originate in the troposphere may be problematic because their upward propagation to the mesosphere and lower thermosphere region could be significantly impeded due to an extended region of strong evanescence above the stratopause. To reconcile this apparent paradox, an alternative explanation is proposed in this paper. The inclusion of mean winds and their vertical shears is sufficient to allow certain short period waves to remain internal above the stratopause and to propagate efficiently to higher altitudes. A time-dependent numerical model is used to demonstrate the feasibility of this and to determine the circumstances under which the mesospheric wind shears play a role in the removal and directional filtering of short period gravity waves. Finally this paper concludes that the combination of the height-dependent mean winds and the mean temperature structure probably explains the existence of short period, quasi-monochromatic structures observed in airglow images of mesopause region.      

一次暴雨激发平流层重力波的卫星观测与数值模拟

针对卫星Aqua/AIRS观测到的与2011年7月25日山东省乳山市特大暴雨相伴的一次平流层重力波过程,利用中尺度数值模式WRF进行暴雨诱发平流层重力波的数值模拟.对模式输出的垂直速度场和温度扰动场的分析表明,暴雨在平流层内的弧状波结构主要集中在降水云系东侧,水平影响范围大于1000km,且随着高度的增加,圆弧状结构趋于闭合,波动能量显著增强.此外,对垂直速度剖面结构分析表明,受高空东风和风切变的影响,重力波在上传过程中逆着背景风场向东传输,不同高度波动形态各异.基于快速傅里叶变换（FFT）的功率谱分析结果表明,此次暴雨激发的平流层重力波在35km高度的周期为7~20h,水平波长约为1000km,垂直波长为5~10km.通过分析动量通量的垂直输送,定量反映出重力波上传过程中的动力学变化特征.      

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.     

基于激光雷达实测和动力学仿真方法研究中间层顶垂直风扰动特性

利用中国廊坊台站钠荧光多普勒激光雷达82h垂直风和水平风观测数据,统计得到中间层顶区域中存在10m·-1量级的垂直风扰动和纬向风扰动,其中垂直风扰动远远超过平均风速为-0.015m·-1的背景垂直风速.根据三维准单色重力波的极化关系和色散关系,对高中低三种频率重力波产生的垂直风扰动进行仿真,结果显示在满足短周期、大纬向风扰动条件下,高频重力波能够产生最大10m·-1量级的垂直风扰动,中频重力波能够产生10m·-1以内的垂直风扰动,低频重力波能产生1m·-1以内的垂直风扰动.理论条件下准单色重力波能够产生10m·-1量级的垂直风扰动,钠激光雷达观测到的最大10m·-1量级的垂直风扰动真实存在.研究结果可对高层大气垂直风场探测、垂直风场模拟和重力波参数化提供依据.      

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.     

Gravity wave propagation studies using the Indian MST radar observations

The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.     

Gravity waves and wind-shear in the MLT at 23°S

We have used the technique suggested by Hocking [Hocking, W. A new approach to momentum flux determinations using SKiYMET meteor radars. Ann. Geophys. 23, 2005.] to derive short period wind variances in the 80–100 km region from meteor radar data. We find that these fluctuating winds, assumed to correspond to gravity waves and turbulence, are closely correlated with the vertical shear of the horizontal tidal winds. This close correlation suggests that in situ wind shear may be a major source of gravity waves and turbulence in the MLT. If this is the case, gravity waves generated in the troposphere and propagating up to the MLT region, generally assumed to constitute an important influence on the climatology of the region, may be a less important source of energy and momentum in the 80–100 km region than has been hitherto believed.