北京地区大气温度及重力波活动的季节变化

利用瑞利激光雷达观测数据,分析了北京地区35~70km高度范围内大气温度和重力波活动的季节变化.发现北京地区30~70km高度范围内的大气温度有明显的年周期变化:平流层顶最高温度出现在6,7月份,大约为270K;中间层70km高度最低温度也出现在6,7月份,大约为200K.以2014年10月14日晚数据为例,分析重力波势能密度,发现50km以下重力波势能存在耗散,而在50km以上重力波近乎无耗散地向上传播.通过对比35~50km高度范围内的平均势能密度,对北京地区重力波活动强弱的季节变化进行了研究.研究结果表明,北京上空重力波活动强度具有明显的年周期变化,冬季平均势能密度为18J·kg-1,夏季为8J·kg-1,且冬季重力波活动强度约为夏季的两倍.此外,还分析了春夏秋冬四个季节重力波势能密度随高度的变化.结果表明,不同季节和不同高度的重力波势能密度不同.      

2009年1月5日廊坊OH气辉观测图像的谱分析

2009年1月5日在中国河北省廊坊市韩村镇(39.4°N,116.6°E)进行了中国首次中间层顶大气重力波成像观测实验,记录了近4 h的OH气辉辐射数据,给出了依据该数据计算所得的该时间段内廊坊韩村附近水平方向约185 km范围内的中间层顶OH气辉扰动的谱分析结果.所得到的传播方向确定的二维功率谱和角度谱揭示了该时间段内重力波活动的波长和传播方向分布特点,表明波长下限较短,西向传播分量明显大于东向分量.     

重力波对中间层和低热层大气环流的影响

利用β通道准地转近似大气平均运动方程组,采用重力波线性饱和参量化方法,定性地研究大气重力 波对中间层和低热层大气环流的作用.模拟计算得到,大气重力波对平均东西风速可产生100m·s-1／d左右的作 用力和产生120 m2/s的湍流扩散,这些作用平衡了Coriolis扭力,导致大气的平均东西风速大大偏离辐射平衡 风场,中层顶附近的平均东西风速在冬季(夏季)甚至反转为东风(西风).平均东西风速计算结果与冬季和夏季 中频雷达东西风速观测值和大气模式剖面等大致一致.     

冷锋期间武汉MST雷达低平流层回波异常现象

武汉MST雷达是中国子午工程建设的两台中高层大气无线电探测雷达之一.该雷达探测频率在VHF频段，雷达回波在低平流层和对流层上部具有角谱特性，可为研究大气动态稳定性提供有效技术手段.本文利用武汉MST雷达2016年4月17日冷锋活动期间及平静天气的角谱实验数据，从雷达回波特性变化、风场空间分布、湍流生成机制以及内重力波影响四个方面分析并解释了MST雷达对流层顶上部区域出现持续异常强回波带的成因.分析结果表明，对流层冷锋的强对流作用诱发内重力波，内重力波向上传播至低平流层后受增强的剪切急流影响发生耗散甚至破碎，激励了长时间跨度的K-H不稳定性，进而导致水平反射层结构发生扰动生成湍流，使得雷达回波结构发生变化.      

重力波波包相互作用时空演变的数值研究

本文设计了一种研究二维空间中重力波波包共振相互作用演变的数值模式.利用此模式获得二个大振幅重力波波包通过碰撞而发生的完整的相互作用过程.数值结果表明,能量上行重力波可以通过波-波相互作用激发能量下行重力波;在相互作用的时空演变过程中,3个波包的能量收支情况和相互作用的强度具有局地性;由于考虑了空间传播效应,相互作用过程不再具有周期性,而是出现一个具有不同物理含义的相互作用特征时间,这个新的相互作用特征时间和相互作用最终的发展程度不仅由相互作用系数和初始波包振幅决定,还受到波包空间尺度和波包群速度的相对大小控制.      

红精灵发光光谱的数值模拟研究

红精灵是发生在雷暴层云顶的一类大气瞬态发光现象, 是能量由对流层耦合到中高层大气的直接证据. 其发光光谱研究是了解整个事件对中高层大气能量注入的重要手段, 有助于认识事件发生区域的大气电离度及事件过程的能量电子分布, 进一步为研究红精灵的产生机制提供重要信息, 同时为大气辐射背景资料研究提供重要依据. 本文利用Boltzmann方程求解了电场作用下弱电离气体中电子能量分布的时变函数, 以此为基础, 模拟计算了红精灵各典型发射带的光谱强度. 模拟计算结果表明, 约化电场E/N越强, 电子获得的能量就越多, 高能电子也就越多, 致使撞击中性大气产生的辐射光强就越强; 模拟显示红精灵光辐射谱分布从远紫外直至近红外.      

北京地区重力波活动及其波谱的季节分布特性研究

利用激光雷达对北京地区上空Na层进行持续观测,通过连续三年累积的夜间观测数据对北京地区重力波活动及其波谱进行研究.根据重力波的线性理论计算,得到北京地区上空的大气密度扰动规律、空间功率谱和时间频率谱.通过选择重力波波长在1~8 km,具有特定波长以及特定周期为60,45,25 min的重力波活动辅助研究重力波的季节变化规律,结果表明北京地区重力波大气密度扰动具有夏季大、冬季小的活动规律.结合波源与背景风场的季节性变化规律,分析得出北京上空重力波活动季节性变化的主要原因为青藏高原地形和对流因素与我国北方地区季节性背景风场共同作用的结果.     

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

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

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

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

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

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

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.     

A simple model for the interhemispheric coupling of the middle atmosphere circulation

The interhemispheric coupling of the middle atmosphere general circulation is characterized by a global anomaly pattern of the zonal-mean temperature. This pattern reflects an anomalous stratospheric and mesospheric residual circulation, in which a weaker (stronger) stratospheric winter circulation is linked to an upward (downward) shift of its upper mesospheric branch reaching from the summer to the winter pole. This phenomenon is robust in observational data and several middle atmosphere general circulation models. In the present study, the recently proposed mechanism of the interhemispheric coupling is unequivocally proven within the framework of a zonally symmetric model that excludes any additional effects due to resolved waves and non-zonally propagating gravity waves. Two simulations are conducted that differ in the strength of the polar vortex. A weaker polar vortex results in a downward shift of the winter mesospheric gravity wave drag. This leads to changes also in the summer upper mesosphere via a feedback solely between gravity wave breaking and the zonal-mean state. The accompanying temperature anomaly reproduces the pattern of the interhemispheric coupling.     

The identification of mesospheric frontal gravity-wave events at a mid-latitude site

Mesospheric frontal-type gravity waves are an uncommon type of wave disturbance that occurs in the mesospheric OH, Na, O₂, and O(¹S) nightglow. They are understood to be the result of gravity waves exhibiting various degrees of non-linear behavior. Despite their similar appearance in all-sky images, careful analysis reveals that there are at least two distinct types of frontal wave disturbances, each with completely different consequences in terms of vertical momentum transport and deposition. Therefore, a correct identification is important in order to characterize their propagation modes. In this report we present the frontal gravity wave activity that occurred during a twelve-month period at Millstone Hill (42.6°N, 172.5°W), a mid-latitude site, to illustrate their range of behaviors.     

First observation of upper mesospheric semi annual oscillations using ground based airglow measurements from Indian low latitudes

We summarize two years of Mesosphere Lower Thermosphere Photometer (MLTP) operation of mesospheric OH and O₂ emission monitoring. The deduced mesospheric OH and O₂ temperatures show large variability. Nightly temperature variations over Gadanki (13.5°N, 79.2°E) are dominated by the short period wave features, while tidal amplitudes are relatively small. Our measurements are the first to report a long period seasonal variation at two upper mesospheric altitudes simultaneously over the Indian sector. Our observations reveal the presence of a dominant semi-annual oscillation (∼6 months periodicity) together with a shorter period (∼2.5  months periodicity) oscillation in both OH and O₂ data.     

Analysis of gravity-wave induced instabilities and turbulence viscosity parameters from optical emissions

The effects of gravity waves on OH and 5577 0I emissions have been used to investigate (1) the possible production of “secondary” waves or “ripples” by a “primary” wave; (2) the possible application of such observed optical emissions for computing atmospheric parameters such as turbulence viscosity in the mesospheric regions.     

Simultaneous bi-station measurements of mesospheric waves from Indian low latitudes

Coordinated mesospheric night-airglow measurements have been carried out from two stations, Gadanki (13.5° N, 79.2° E) and Allahabad (25.5° N, 81.9° E), India during April 2009 to study the common gravity wave features. With two nights of coordinated measurements we find some of the wave periodicities to be similar at the two locations. Simultaneous OH and O₂ intensity measurements over Gadanki reveal these features to be upward propagating gravity waves while the coordinated OH intensity measurements of similar waves from Allahabad show the large spatial extent of these waves.     

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.      

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.     

武汉上空中层顶区域潮汐的MF雷达观测

利用武汉中频雷达观测数据进行分析,研究2至3月份武汉上空中层顶潮汐结构及其随高度和时间变化的特性.用Lomb-Scargle周期图方法计算的水平风场动态功率谱表明,武汉上空存在持续的周日潮汐,是中层顶区域风场结构的主要成分.周日潮的平均振幅随高度的增加呈先增后减的趋势.大多数情况下,潮汐谱峰对应的频率与定义值有一定的偏移.周日潮水平扰动速度矢量随时问和高度变化的轨迹表明,经向分量的相位比纬向分量的相位超前,潮汐能量向上传播,对应于向下的相位传播速度.计算得出的经向分量和纬向分量的垂直相速度分别为1.10和1.15 km/h.