海口上空中高层大气强烈的突发钠层观测

海南激光雷达站（20°N,110°E）2011年11月2日观测到强烈的突发钠层（SSL）.该SSL峰值密度达37087cm-3,半峰全宽仅为0.9km.对该SSL发生前后的钠层峰值密度和高度特性进行分析.统计激光雷达2010年5月至2013年12月共计377天观测到的SSL事件222次,其中仅有1次SSL峰值密度超过30000cm-3.对比距离最近的海南儋州（19.5°N,109.1°E）测高仪和VHF雷达观测到的突发E层（Es）事件,分析事件的相关性.Es最低高度与SSL峰值高度差均在5km以内,约有75%的Es与SSL时间差在±30min之内.实验中使用continuum的Nd:YAG激光器和泵浦染料激光器产生589nm激光,能量为45mJ,使用直径为1000mm的望远镜接收钠层的光子回波.      

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

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

重力波波包在真实大气中传播特性的数值研究

采用二维全隐欧拉（FICE）格式，对重力波波包在真实大气中的非线性传播和演变过程进行了数值模拟，模拟结果表明，在中层大气下部激发的向上传播的重力波波包在传播到中层顶之前，波相关能量沿着射线路径传播，非线性效应和背景温场对波能量传播路径的影响很小，当波包传播到低热层大气后，波包饱和波相关能量几乎完全沿着水平方向传播，垂直方向的能量传播受到抑制，这与在无耗大气下，WKB近似条件下的线性重力波理论的预言相差很大，深入的分析表明抑制重力波波包向上传输能量的关键因素是大气分子粘性的垂直非均匀性。非线性和背景温场的影响不足以完全抑制波能量的向上传播，此外，在波包的整个传播过程中，由于非线性，背景温度和背景耗散的共同作用，重力波波包的垂直波长随时间明显减小，这些结果说明大气的分子粘性特别是分子粘性的垂直不均匀性对重力波波包在中、高层大气的非线性传播过程起着重要作用。     

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

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

我国地区20—80km高空大气温度特征

利用Nimbus－7卫星1979－1981年的平流层和中间层大气温度探测（SAMS）数据，首次较全面地分析了我国上空20－80km高度大气温度的分布特征，并与COSPAR国际参考大气CIRA－1986进行了比较。结果指出，我国上空大气温度分布与COSPAR国际参考大气CIRA－1986的纬向平均温度分布存在一定差别，可供航空航天工作者参考。     

2012年7月海口与北京大气钠层观测

2012年7月海南激光雷达站（20°N,110°E）和北京激光雷达站（40°N,116°E）同时进行观测,对比低纬和中纬度中高层大气钠层分布及密度变化特性.两地使用相同的Nd:YAG激光器和泵浦染料激光器产生589nm激光,设定相同的激光参数,使用相同的直径1000mm的望远镜接收钠层的光子回波.两地均观测到突发钠层和双钠层事件.其中:海口观测到突发钠层事件9个,双钠层事件1个;北京观测到突发钠层事件1个,双钠层事件4个.分析海口和北京同一天夜间观测的大气钠层特性,对比两地同时观测到的钠层高度与密度的不同变化,发现当天夜间低纬度地区海口市上空钠层出现双钠层结构,持续时间123min,但是北京上空钠层没有发现该结构.      

三维重力波非线性传播数值模型的一种求解格式

采用交替方向隐格式（ADI）求解三维可压大气中重力波非线性传播的数值模型，经过几个波周期后，重力波波包括估然沿着线性理论射线方 稳定地传播，波包保持着很好的形状，内点和边界保持稳定；水平风与垂直运动的振幅随着高度的增加而增加。这说明ADI格式是模拟三维可压大气中重力波传播的一种可靠的格式，模拟结果表明，在非线性情况下，重力波仍可保持性传播的某些特征，非线性效应有现在对波包括传播速度的修正。     

极区夏季70—100km大气垂直速度扰动功率谱的观测

使用1987年夏季MAC/SINEcampaign期间在挪威北部的Andenes发射的5枚chaff火箭垂直速度数据研究了70—100km高度范围内大气垂直速度扰动垂直波数谱.发现垂直速度垂直波数谱在39×10-4m－1处有1个谱缺口.在大于这个谱缺口的垂直波数范围,垂直速度垂直波数谱与饱和谱模式有好的一致.在小于这个谱缺口的垂直波数范围,垂直速度垂直波数谱有谱斜率－1.87和谱振幅比饱和谱振幅低1个量级,认为在低波数,除重力波外的其它地球物理过程存在于垂直速度扰动中.      

武汉中层大气中频雷达及其初步探测结果

首先简要地讨论了武汉中频雷达观测原理和设备的组成，该雷达测量60-100km高度的大气风场和电子密度，风场采用分布天线测量技术和全相关分析方法得到，电子密度通过微分吸收和微分相位技术获得，初步观测结果表明：（1）武汉上空冬季60-100km高度的纬向风多为西风，风速为30-50m/s，经向风速为10-20m/s，垂直风速较小，一般在5m/s以内，（2）60-100km高度范围的大气风场和电子密度均有明显的日变化，风场在某些时段和高度区间有较强的风剪切出现。（3）80km以上高度大气的风场和电子密度存在较明显的扰动现象，它可能与大气波动过程有关。     

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

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

重力波波包在三维非等温大气和水平风场中非线性传播的数值研究

采用全隐欧拉格式（FICE）对重力波波包在三维非等温大气、均匀和剪切风场中的非线性传播进行了数值模拟，给出了重力波波包三维非线性传播的全过程，分析了重力波的传播特性及背景温度场、风场对重力波传播的影响。结果表明：波包扰动速度振幅的增长比在WKB条件下振幅的增长要慢；波包非线性传播的路径、能量传输速度不同于WKB近似下的结果，非线性效应导致了重力波的传播特性的改变；温度场的非均匀性会改变重力波传播的路径和速度；剪切风场使扰动速度振幅的增加变得缓慢，使垂直波长减小。     

The VEGA Venus balloon experiment

In June 1985, two instrumented balloons were placed in the atmosphere of Venus as part of the VEGA mission. Each balloon traveled about 30 percent of the way around the planet at a float altitude near 54 kilometers. In situsensors measured pressure, temperature, vertical wind velocity, cloud particle backscatter, ambient light level, and frequency of lightning. A ground-based network of 20 radio antennas tracked the balloons by very long baseline interferometry (VLBI) techniques to monitor the Venus winds. The history, organization, and principal characteristics of this international balloon experiment are described.     

用电离层特性参量提取等效风场信息

导出了利用中低纬电离层特性参量获取电离层F层峰区高度上等效风场（包含电场和风场信息在内）的基本方程，并尝试用该方法从电离层特性参量（峰高和临频）提取等效风场信息，利用武汉站DGS-256电离层数字测高仪数据及由美国Massachusetts Lowell大学最新版的剖面反演程序换算得到F层峰高，获得了武汉地区夏季至日点附近，冬季至日点附近，冬季地磁特别宁静的九天和冬季平均等效风场的初步特征，并利用Fejer经验电场模式计算冬季电场引起的垂直漂移，估计电场和风场对武汉地区的垂直等效风场的贡献大小，结果表明：等效风场呈现出白天与夜晚幅度和方向的差异。至日点附近冬季与夏季白天的幅度差异以及明显的凌晨凹陷现象；平均情况下，垂直等效风场幅度和方向的变化主要是由中性风引起，受电场的影响不大。     

General aspect of GPS data use for atmospheric science

The radio link between a GPS satellite and a GPS receiver is appropriate for bistatic radar sounding of the Earth's atmosphere, ionosphere, and ocean surface (latter in case of GPS reflection). Measurements of GPS phases and amplitudes are currently performed by spaceborne, airborne, mountain- and ground-based GPS receivers. In the present paper, an uniform approach based on geometrical optics and spherical symmetry of the atmosphere is applied to various GPS observation configurations. Atmospheric mapping function, influence and retrieval of ionospheric layers/disturbances, tropospheric water vapor, and possible measurement of vertical winds and wave velocities are investigated by use of simulation data of GPS phase path excess and bending angle.     

Initial estimates of Venus winds from ground-based radio tracking of the VEGA balloons

A global array of 20 radio observatories was used to measure the three-dimensional position and velocity of the two meteorological balloons that were injected into the equatorial region of the Venus atmosphere by the VEGA spacecraft. Initial analysis of only radial velocities indicates that each balloon was blown westward about 11,500 kilometers (8000 kilometers on the night side) by zonal winds with a mean speed of about 70 meters per second. Excursions of the data from a model of constant zonal velocity are generally less than 3 meters per second; however, a much larger variation is evident near the end of the flight of the second balloon. Consistent systematic trends in the residuals for both balloons indicate the possibility of a solar-fixed atmospheric feature.     

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

通过分析武汉、宜昌和恩施气象局无线电探空仪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.      

中高层大气对低层大气脉冲扰动的响应

从基本的大气运动方程出发,采用二维全隐欧拉格式(FICE)建立极坐标下二维声重波传播的数值模式.对小振幅高斯形态声重波传播过程的模拟结果表明,该数值模式能够再现声重波波包向上的稳定传播过程.应用此数值模式模拟了在没有背景风场和存在背景风场条件下中高层大气对10km高空处近似脉冲点源扰动的响应过程,获得了中高层大气对低层大气脉冲扰动响应的图像.结果表明,电离层高度对重力波扰动的响应不但出现在激发源激发之后,还与激发源具有较大距离,扰动的幅度随高度的增大而增大;有背景风时,顺风情况下重力波波动传播路径较逆风情况下平缓,且传播的水平距离比逆风情况下远,波动的振幅较逆风情况下弱,逆向背景风场能够加速重力波向上传播.同时,模拟还表明脉冲扰动正上方可以激发产生一种以6 min为主要周期,具有声学分支特性的快速波动.本文模拟的扰动结果与2004年12月26日苏门答腊-安达曼地震引起的电离层TEC扰动具有类似的图像.      

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

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

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.