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

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

Secular variations of OH nightglow emission and of the OH intensity-weighted temperature induced by gravity–wave forcing in the MLT region

Secular variations of OH airglow (8,3) band and of the intensity-weighted temperature induced by a gravity wave packet at three latitudes in the northern hemisphere in the Mesosphere/Lower Thermosphere region are investigated using a 2-D, time-dependent, fully nonlinear OH chemistry–dynamics model. A net-cycled average survives due to wave transience and dissipation. The integrated OH (8,3) volume emission rates show a 20 percent increase induced by the wave packet. The corresponding intensity-weighted temperature shows an initial increase up to 2 percent then a gradual decrease for the remainder of the simulation time by the same wave packet. These secular variations could be mistaken as long-period or short-period waves in the airglow observations. Therefore, care must be taken when analyzing the data from observations.     

A comparison between CRISTA satellite data and Warner and McIntyre gravity wave parameterization scheme: Horizontal and vertical wavelength filtering of gravity wave momentum flux

An along-track analysis of temperature height profiles measured by the CRISTA satellite instrument was performed to obtain height profiles of absolute values of gravity wave momentum flux. Global distributions of momentum flux from CRISTA are compared to results of the Warner and McIntyre gravity wave parameterization scheme. Only a limited range of horizontal and vertical wavelengths is visible to the CRISTA instrument. Therefore the model results have to be filtered accordingly to make comparisons meaningful. We demonstrate the importance of such a wavelength filtering which significantly increases the correlation of observed and modeled momentum fluxes.     

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

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

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

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

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

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

Wave perturbations in the MLT at high northern latitudes in winter,observed by two SATI instruments

To investigate the Mesosphere and Lower Thermosphere (MLT) region, several ground-based instruments called SATI (Spectral Airglow Temperature Imager) were designed and built to measure airglow emission and temperature in the upper mesosphere. One SATI instrument was installed at Resolute Bay (74.7°N, 94.9°W) and has monitored the polar MLT region since November, 2001. In October 2007 another SATI instrument was installed at Eureka (80.0°N, 86.3°W) at the Polar Environment Atmospheric Research Laboratory (PEARL) as part of the Canadian Network for the Detection of Atmospheric Change (CANDAC) project. SATI is a spatial scanning Fabry–Perot spectrometer measuring column emission rates for several rotational lines of OH and O₂ airglow at 87 and 94 km height. The rotational temperatures are inferred from the ratios of these lines. The measurements are divided into 12 sectors with an annular field of view. The phase differences between the sectors yield information on the horizontal atmospheric wave direction and wavelength. Horizontal perturbations of 2–8 h period have correlatively been observed and investigated at both locations. Short-periodic oscillations identified as gravity waves with periods between 2 and 8 h propagate in southward and eastward directions, but in opposite directions in some cases. The wave propagation characteristics are often different at the two locations; the relationship with the lower mean wind is considered.     

日全食期间中间层大气臭氧光化过程试验

利用包括活O_x、HO_x、和NO_x等成份在内的一维时变光化模式,研究了日全食时上中层大气臭氧和O₂(1△g)的变化特性,并通过与观测资料的比较进行了光化模式的试验.结果表明,日全食时太阳辐照变化对1.27μm气晖影响的高度比曙暮期间明显偏低,利用1.27μm气晖光度计测量大气臭氧的高度范围约为52-67km.1997年3月9日漠河日全食时测量的大气臭氧含量变化幅度为50％.这与考虑O_x、HO_x、和NO_x等成份的光化模式理论模拟结果基本一致.它表明在52-67km高度范围不存在臭氧亏损现象.      

Comparison of calculation models for determination of the mesopause temperature using SATI images

The Spectral Airglow Temperature Imager is an instrument for ground-based spectroscopic measurements of the night-glow atmosphere emissions. This instrument was developed specially for gravity wave investigation. The measured airglow spectra are matched to synthetic spectra calculated in advance for determination of the temperature in the mesopause region where the radiation maximum of some О₂ emissions is situated. The synthetic spectra are transformed into a format which corresponds to the measured spectra in order to be matched. This transformation is based on the known values of the refractive index and the central wavelength of the interference filter used. A substantial part of the processing algorithms of the SATI images is connected with determination of these two filter parameters. The results of the original and newly-proposed algorithms for filter parameter calculation and their importance for the final results for temperature determination on the basis of the О₂ (864–868 nm) emission measurements are presented.     

Adelaide地区大气重力波的气辉观测

利用澳大利亚Adelaide(34.5°S,138.5°E)地区7年(1995-2001年)的OH和OI气辉观测数据,以及同地区中频雷达所测背景风场数据,通过最大熵谱分析方法,统计分析了该地区大尺度重力波的传播特性.结果表明,(1)在中层顶气辉观测区,利用OH气辉探测到48个重力波事件,利用OI气辉探测到29个重力波事件,这些重力波的周期和水平相速度分别集中在68 min和36 m/s左右,具有东南向的优势传播方向,平均仰角为8°～10°,接近水平方向传播;(2)通过OH气辉和OI气辉观测到了同一个重力波事件,并根据两个高度上观测到的扰动相位差判断该重力波为上行波.      

氧气A(0,0)波段气辉体发射率和临边辐射强度模拟与分析

临近空间大气参数如温度、密度、风场等对预报模型精度及航天器运行安全等有较大的影响，而气辉的辐射模拟是大气参数反演的重要过程.本文基于光化学模型计算了氧气A（0，0）波段气辉的体发射率和临边辐射强度.基于氧气A（0，0）波段气辉的光化学反应机制、大气动力学和光化学反应理论，建立产生O₂（b1Σ_g+）的光化学模型.计算气辉体发射率，基于临边探测几何路径进行气辉辐射强度模拟.体发射率计算结果与AURIC模型结果的辐射值及辐射高度均一致.基于计算和模拟结果，对氧气A波段气辉体发射率和辐射强度的影响因素进行了分析.      

A review on the upper atmospheric sodium observations from India: Insights

This review article briefly brings out the historical development of atmospheric sodium (Na) measurements over India and the importance of coordinated measurements with multiple techniques to address physical processes in the Earth's upper atmosphere. These measurements were initiated in the early 1970s by observing Na airglow emission intensities with broad band airglow photometer from Mt. Abu, a low-latitude hill station in India. Considerable amount of night-to-night variations in nocturnal emission intensities of the Na airglow were observed. Later, investigations regarding the dependence with the magnetic activities from the equatorial and low latitudes were carried out and double-humped structures in the nocturnal variation of intensities were reported. With the advent of Na lidar at Gadanki around 2005, the measurements of atmospheric neutral Na atom concentration became possible and more frequent occurrences of sporadic Na layers over the magnetic low latitude station compared to other latitudes were detected indicating the role played by electrodynamics. Later, a possible relationship between E-region field aligned plasma irregularities and the concentration of neutral Na atoms was investigated using coordinated measurements of VHF radar and Na lidar. Further, simultaneous measurements with Na lidar and a narrow band airglow photometer with narrow field of view brought out the importance of coordinated observation wherein the characterization of gravity waves could be carried out and also revealed the importance of collisional quenching due to ambient molecules in the Na airglow emission processes. In addition, combining the ground based measurements of Na lidar and meteor wind radar along with satellite measurements made possible to hypothesize the over-turning Kelvin–Helmholtz billow in the Na layer manifests “C-type” or inverted lambda shape structures in the height-time-concentration map of neutral Na atoms. This review paper presents a synoptic view mostly based on the previously reported observations of Na airglow emission, Na lidar and coordinated Na airglow and Na lidar observations from the Indian sector and highlights the importance of simultaneous measurements of mesospheric Na and its emissions along with satellite-borne measurements to address interesting geophysical processes in the Earth's upper atmosphere.     

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

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

135.6nm夜气辉峰值电子密度反演算法及误差

根据夜间135.6nm大气辉光光强与F₂层峰值电子密度N_mF₂平方成正比的物理机制,在前期夜间135.6nm气辉辐射激发模型研究的基础上建立了峰值电子密度的反演算法,把全球经纬度分成若干格点,每个格点的电离层及中性成分信息分别由IRI2000和MSISE90提供,将电离层及中性成分廓线输入夜气辉辐射激发模型,计算每个格点135.6nm气辉的辐射强度,然后将各个格点的135.6nm气辉辐射强度与电离层廓线输入的N_mF₂平方拟合得到气辉强度与N_mF₂的转换因子.利用此方法可获得不同地方时、季节和太阳活动周期的转换因子组成查算表,进而根据实际探测的135.6nm气辉辐射强度反演相应时空的N_mF₂.最后对该算法的反演误差进行了综合分析,为该算法适用的时空特性提供重要理论支撑.      

重力波波包在剪切风场中的非线性传播

采用二维全隐欧拉格式,对具有高斯分布的重力波波包在剪切风场中的非线性传播过程进行数据模拟和分析,数值分析结果表明,在没有出现临界层的情况下,尽管存在非线性效应,在整个传播过程中,波动的等相面向下运动,波包和波相关能量仍然能够自由地向上传输。波相关扰动速度随高度增加而增加,并且垂直波长随高度增加而减小。     

Thermospheric F-region travelling disturbances detected at low latitude by an OI 630 nm digital imager system

Gravity wave effects in the nocturnal thermospheric F-region domain are seldom detected in the intertropical region by optical (airglow) techniques, especially during geomagnetically quiet times, in part because the low inclination of the magnetic field, as opposed to the case of the mid-latitude region, does not favor significant vertical excursions of ionospheric plasma in response to meridional winds. Such difficulty in detecting gravity wave signatures in the F-region by means of optical techniques tends to increase in the absence of geomagnetic storms because of the lack of strong forcing mechanisms necessary to generate high intensity gravity waves. The purpose of this work is to show that during the quiet day of 9 August 1999, the Terminator may have been a source region of wave-like disturbances in the nocturnal F-region at the low-latitude station Cachoeira Paulista (22°41'S; 45°00W, dip 30°). A digital all-sky OI 630nm imager system located at that station has shown propagating wave-like spatial structures in the airglow intensity near the Terminator. This observation supports a previous study on the evidence of the presence of gravity waves during the post-sunset period at the same location by means of a scanning photometer system (1997, Sobral, J. Atmos. Terr. Phys. 59, 1611–1623). The absence of range-type spread-F as monitored by a local digisonde and the absence of radio wave scintillation as monitored by a local GPS receiver, excludes the hypothesis that the wave-like airglow structures are associated with the occurrence of the ionospheric plasma bubbles. Downwards motion of the iso-density real height contours at 22.0 ms⁻¹ and 33.1 ms⁻¹ are observed. The wave detection by the imager system is reported and discussed here.     

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

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

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.     

大气重力波引起的突发钠层的模拟研究

应用钠原子和中性大气分子的质量连续性方程来模拟突发纳层（SSL），垂直风场采用接近实际大气重力波的正弦行波模式，结果较好地反映了SSL的形成过程。SSL的形成时刻在5－15min之间并可持续到30min之后，形成高度大约在90－100km之间，峰宽为0.5-2km之间，这些都与实际观测SSL的特点相符，同时还进一步地研究了当重力波参数（主要指垂直波长和周期）、风速以及常态钠层半宽度发生变化时SSL的变化趋势。     

非线性效应对重力波特性的影响

通过分析二维情况下重力波波包在等温大气中非线性传播的数值模拟结果,对非线性效应对重力波特性的影响做了定量分析。分析结果表明;重力波在非线性传播过程中,一方面由于非线性作用,随着波振幅的增加,垂直波长减小,同时平流成份增加;另一方面,尽管存在非线性效应,线性重力波的色散关系与偏振关系仍然适用。