探测大气温度和气溶胶的瑞利-拉曼-米氏散射激光雷达

研制了一台瑞利-拉曼-米氏散射激光雷达,实现了对流层和平流层大气温度和密度的探测。作为多参数大气探测系统,该激光雷达也实现了夜间至25km、白天至5km高度气溶胶的探测能力;其中激光雷达是探测平流层气溶胶最有效的手段之一。利用该激光雷达对目前合肥地区对流层温度、平流层逆温现象、对流层和平流层气溶胶做了探测和分析,并给出若干典型结果。分析表明,该激光雷达数据可靠,可用于大气温度、密度、气溶胶的常规观测和分析研究。     

用于大气温度廓线测量的瑞利-拉曼激光雷达

大气温度廓线及其时间演变特征资料在地球科学领域具有重要的应用,为获取高时空分辨的大气温度的垂直分布,建立了瑞利-拉曼温度测量激光雷达。介绍了瑞利-拉曼激光雷达进行温度测量的主要原理和研制的瑞利-拉曼激光雷达的主要参数;数据处理方面,通过背景噪声剔除和小波算法降噪提高系统的信噪比;使用研制的激光雷达对大气温度廓线进行观测,将观测结果与大气模式数据和卫星观测结果进行对比,均显示较好的一致性,证明了激光雷达温度测量结果的准确性,其温度测量数据可以用于气象学研究。     

基于纯转动拉曼谱线激光雷达的大气温度反演分析

由于气溶胶的影响,传统的瑞利散射法测量低空大气温度有一定的局限,为此开展了纯转动拉曼法测量低空大气温度。利用纯转动拉曼激光雷达在北京进行了2个月的大气温度观测,由观测数据反演了温度廓线。在基于N₂和O₂的纯转动拉曼谱线特征进行大气温度反演过程中,分析了平滑窗口、定标范围和定标常数对温度反演精度的影响。结果显示随着平滑窗口的增大,雷达和无线电探空仪测量的温度之间的平均绝对偏差先减小后增加,为有效去除信号中随机误差的影响,同时保留温度廓线的垂直结构,平滑窗口应选择600~1 200 m比较好。定标范围不同,雷达和无线电探空仪测量的温度之间的平均绝对偏差就不同,相对变化约为0.07 K。当定标常数a,b都增大或都减小时,雷达和无线电探空仪测量的温度之间的平均偏差增大,当一个增大另一个减小时,平均偏差相互抵消;a,b的变化不是等几率的,在符号上总是相反的;平均偏差对a的变化不敏感,对b的变化也不敏感,对a与b的整体变化敏感,约91.7%平均偏差落入-3~3 K之间。该研究分析结果对纯转动拉曼激光雷达数据反演中涉及的平滑窗口、定标范围的选择提供了理论依据,对激光雷达定标常数造成实际温度反演结果的误差提供了参考。     

北京一次小震级地震前异常声-重力波产生机理的研究

观测到北京一次小震级地震(ML=3.0)前的异常声-重力波,并同时记录到异常地表倾斜信号。二者相关系数比较高,在接近震中的地区相关系数达到最大值0.5,相关系数随着远离震中的距离增大而减小。经过研究表明,异常声-重力波的产生与震前的地表缓慢活动引发的山体缓慢晃动导致的气流风速波动变化有关。基于350 m高的气象铁塔观测的风速数据,对两层大气模型传播的声-重力波进行了数值计算,得到的相对山体非平稳气流引起的声-重力波与观测信号在特征周期、幅度和速度三个参数上近似。以地表倾斜为例的模型与观测数据的一致性表明该震前声-重力波与地震存在关联性。      

激光雷达探测北京城区夏季大气边界层

为了研究北京城区夏季大气边界层结构变化特征及大气边界层内气溶胶消光特性,2004年8月利用便携式米散射激光雷达对北京城区夏季大气边界层进行了系统观测。反演了观测站上空大气气溶胶的消光特性垂直分布以及大气边界层的高度。分析了气象条件和人类活动对大气边界层结构的影响。观测数据表明： 北京城区夏季大气边界层有明显的日变化特征,早晚比较低,日间有一个从低变高再变低的过程,中午前后达到最高。结合气象参数对测量数据进行的统计分析表明：北京城区夏季大气边界层高度相对稳定,多分布在1.8km以下,平均值为0.68km;大气边界层内存在浓度较高的气溶胶粒子,平均光学厚度（3km以内）在0.30左右。     

激光雷达测量大气湍流廓线

 介绍了一种利用成像激光光斑测量大气湍流廓线的激光雷达原理。通过对分层大气湍流的光束波面变化的测量,获取各分层大气湍流的相干长度,据此利用平面波近似算法反演湍流强度廓线。通过搭建的大气湍流廓线激光雷达实验系统获得了湍流廓线的实验数据,并且与系留飞艇搭载的温度脉动仪在同时段固定高度进行了对比实验,验证了近似递推算法的可行性。最后对实验中出现的误差进行了讨论。     

基于纯转动拉曼谱线激光雷达的大气温度反演分析（英文）

由于气溶胶的影响,传统的瑞利散射法测量低空大气温度有一定的局限,为此开展了纯转动拉曼法测量低空大气温度。利用纯转动拉曼激光雷达在北京进行了2个月的大气温度观测,由观测数据反演了温度廓线。在基于N_2和O_2的纯转动拉曼谱线特征进行大气温度反演过程中,分析了平滑窗口、定标范围和定标常数对温度反演精度的影响。结果显示随着平滑窗口的增大,雷达和无线电探空仪测量的温度之间的平均绝对偏差先减小后增加,为有效去除信号中随机误差的影响,同时保留温度廓线的垂直结构,平滑窗口应选择600~1 200m比较好。定标范围不同,雷达和无线电探空仪测量的温度之间的平均绝对偏差就不同,相对变化约为0.07K。当定标常数a,b都增大或都减小时,雷达和无线电探空仪测量的温度之间的平均偏差增大,当一个增大另一个减小时,平均偏差相互抵消;a,b的变化不是等几率的,在符号上总是相反的;平均偏差对a的变化不敏感,对b的变化也不敏感,对a与b的整体变化敏感,约91.7%平均偏差落入-3~3K之间。该研究分析结果对纯转动拉曼激光雷达数据反演中涉及的平滑窗口、定标范围的选择提供了理论依据,对激光雷达定标常数造成实际温度反演结果的误差提供了参考。     

全天时紫外高光谱瑞利测温激光雷达系统EI北大核心CSCD

紫外高光谱瑞利测温激光雷达是一种探测大气温度廓线的有效工具。目前,紫外高光谱瑞利测温激光雷达通常采用355nm波长的光,然而白天太阳背景光辐射会影响雷达系统的信噪比(SNR),进而制约温度探测的距离和精度。针对大气温度的全天时探测,提出了基于法布里-珀罗标准具的266nm紫外高光谱瑞利测温激光雷达系统。由于到达地面的太阳背景光辐射不包含266nm波长的光,只需考虑臭氧对266nm波长光吸收的影响,进而实现全天时大气温度的探测。基于脉冲能量、望远镜直径、望远镜接收视场角、臭氧浓度以及太阳背景光强度等主要影响参数,对266nm和355nm两个波长紫外高光谱瑞利测温激光雷达系统的谱宽、透过率、回波信号SNR以及温度偏差参数进行数值仿真和对比分析。结果表明,大气分子和气溶胶散射对266nm波长光的影响远大于对355nm波长光的影响。白天266nm紫外高光谱瑞利测温激光雷达系统的有效探测距离为4km左右,比355nm紫外高光谱瑞利测温激光雷达系统的有效探测距离远2.9km;夜间266nm紫外高光谱瑞利测温激光雷达系统有效探测距离为6km。探测距离小于5km时,白天266nm紫外高光谱瑞利测温激光雷达的探测温度偏差比355nm紫外高光谱瑞利测温激光雷达的探测温度偏差小10K。266nm紫外高光谱瑞利测温激光雷达可实现全天时大气温度的探测。     

基于Fizeau干涉仪的激光雷达测温方法研究

气温是描述大气状态的基本参数之一,温度的准确测量对天气预报、气候预测及其他气象参数的反演都至关重要。激光雷达作为一种遥感仪器,已经用于气象要素的探测中（风、温度、气溶胶的光学厚度等）。目前,测温激光雷达主要有拉曼激光雷达（振动和转动）、共振荧光激光雷达和Rayleigh散射激光雷达等,拉曼激光雷达需要大功率的激光器和复杂的背景滤波器;共振荧光激光雷达无法探测平流层内的温度;基于Rayleigh散射的测温激光雷达多应用于温度的相对测量,反演温度时需要建立响应函数和校准程序;基于固体腔扫描F-P干涉仪测量大气Rayleigh散射光谱来反演温度的方法,时间分辨率较低,并且该方法在测量过程中需要运动部件,所以不利于星载。在大气低层,分子的Rayleigh散射光谱会受到Brillouin散射的影响,两种散射信号叠加形成的Rayleigh-Brillouin散射光谱不再服从Gaussian分布,直接通过测量散射光谱的半高全宽来反演温度,会产生误差。基于回波能量的方法会受到气溶胶Mie散射信号的影响,所以在对流层中该方法并不适用。为了实现对流层内温度的高精度和高时间分辨率的测量,提出利用Fizeau干涉仪和PMT阵列对对流层内分子的Rayleigh-Brillouin散射光谱进行测量,并通过插值的方法来对回波信号中气溶胶Mie散射信号进行抑制,从而使Mie散射信号对温度反演的影响较小,最后将测量光谱和理论光谱进行全光谱匹配来实现温度的反演。除此之外,还对Fizeau干涉仪的自由光谱区、固体腔几何长度、腔体反射率、扫描间隔等参数进行了优化设计。为了验证本文提出方法的可行性,利用Matlab软件建立了一套仿真模型,通过模拟表明,在不考虑云、风和水汽含量的影响时,利用该方法测量对流层内的大气温度时,测量误差小于1 K。该测温方法可以对对流层内的大气温度廓线实现高精度、高时间分辨率的测量, 在测量过程中不需要使用运动部件,有较高的使用价值,并对同类高光谱激光雷达分光系统的研究具有借鉴意义, 为我国高光谱激光雷达陆基及星载应用提供了一套可行的技术方案和温度反演方法。     

北京一次小震级地震前异常声-重力波产生机理的研究

观测到北京一次小震级地震(ML=3.0)前的异常声-重力波,并同时记录到异常地表倾斜信号。二者相关系数比较高,在接近震中的地区相关系数达到最大值0.5,相关系数随着远离震中的距离增大而减小。经过研究表明,异常声-重力波的产生与震前的地表缓慢活动引发的山体缓慢晃动导致的气流风速波动变化有关。基于350 m高的气象铁塔观测的风速数据,对两层大气模型传播的声-重力波进行了数值计算,得到的相对山体非平稳气流引起的声-重力波与观测信号在特征周期、幅度和速度三个参数上近似。以地表倾斜为例的模型与观测数据的一致性表明该震前声-重力波与地震存在关联性。     

Simulation of Internal Gravity Wave Propagation Due to Sudden Stratospheric Warming

The results of modeling the effect of stratospheric internal gravity waves (IGWs) excited in the region of sudden stratospheric warming (SSP) on the state of the upper atmosphere are presented. The numerical experiment used a two-dimensional model of atmospheric wave propagation including the dissipative and nonlinear processes accompanying the wave propagation. As a perturbation source, perturbations of temperature and density in a localized region at stratospheric altitudes during the SSW were considered. The amplitude and frequency of the perturbation source were evaluated from the results of observations and IGW theory. The results of numerical calculations showed that the heat source localized in the stratosphere excites IGWs with periods of up to several hours, which reach the thermospheric altitudes within a few hours. The maximum relative perturbations created by these waves relative to the unperturbed conditions were noted at altitudes of 100–200 km and at distances of up to ~1000 km from the source center.     

Features of interannual variations of ozone in the middle stratosphere over Moscow according to observations at millimeter waves

The features of interannual variations of ozone in the middle stratosphere over Moscow in cold half-years since 1995 to 2015 are considered. These features are most pronounced in two separated groups of winters (six winters in each group) in December–January. It appeared that the seasonal variation of ozone averaged within each group is characterized by a higher ozone concentration \({C_{{O_3}}}\) in group I (in comparison with group II) in December and lower values in January. The differences in the seasonal variation of ozone between these groups at the level of 10 mbar in December and January exceed 2 ppm. A feature of the stratosphere circulation for half-year group I is the existence of a stable polar vortex and the absence of strong midwinter perturbations and sudden stratospheric warmings. On the contrary, for cold half-years of group II, the appearance of strong stratospheric warmings in January–February is in common. The results presented provide a quantitative estimate of the effect of these stratospheric warmings on ozone of the middle stratosphere overMoscow.     

Characteristics of atmospheric gravity waves observed using the MU (Middle and Upper atmosphere) radar and GPS (Global Positioning System) radio occultation

The wind velocity and temperature profiles observed in the middle atmosphere (altitude: 10–100 km) show perturbations resulting from superposition of various atmospheric waves, including atmospheric gravity waves. Atmospheric gravity waves are known to play an important role in determining the general circulation in the middle atmosphere by dynamical stresses caused by gravity wave breaking. In this paper, we summarize the characteristics of atmospheric gravity waves observed using the middle and upper atmosphere (MU) radar in Japan, as well as novel satellite data obtained from global positioning system radio occultation (GPS RO) measurements. In particular, we focus on the behavior of gravity waves in the mesosphere (50–90 km), where considerable gravity wave attenuation occurs. We also report on the global distribution of gravity wave activity in the stratosphere (10–50 km), highlighting various excitation mechanisms such as orographic effects, convection in the tropics, meteorological disturbances, the subtropical jet and the polar night jet.     

External regulation of convection in a weak electrically conducting non-Newtonian liquid with g-jitter

The effect of time-periodic temperature/gravity modulation at the onset of magneto-convection in weak electrically conducting second order liquids is investigated by making a linear stability analysis. The Venezian [J. Fluid Mech. 35 (1969) 243] approach is adopted in arriving at the critical Rayleigh and wave numbers for small amplitude temperature/gravity modulation. The temperature modulation is shown to give rise to subcritical motion and gravity modulation leads to delayed convection. Subcriticality or super-criticality depends on the sign of a correction Rayleigh number being negative or positive. An asymptotic analysis is also presented for small/large frequencies. A similar role is shown to be played by the Prandtl number. Comparison is made between the effects of temperature, gravity and combined (temperature+gravity) modulations.     

On the possible perturbation of stratospheric dynamics due to Pinatubo aerosols

Summary The injection of a large amount of SO₂ in the tropical lower stratosphere following the eruption of Mt. Pinatubo (15.1^oN, 120.4^oE; June 15, 1991) can produce a not negligible perturbation in the stratospheric dynamics. Sulfate aerosols formed by nucleation of H₂SO₄ vapour, which is in turn produced from the SO₂ plume, are responsible for a net heating of the order of 0.15K/day in the equatorial stratosphere, located where the thickest portion of the aerosol cloud is found. As shown by recent satellite data, this diabatic perturbation is responsible for a temperature increase of the order of (2÷3) K at 30m bar during the summer and fall of 1991, that is similar to what found for the 1982 eruption of El Chichón. In this work we use a stratospheric three-dimensional model to investigate the dynamical response of the middle atmosphere to this diabatic forcing. Few months after the eruption, change in the diabatic circulation larger than 10% is found in the lower stratosphere, along with a stronger planetary wave activity during the 1991/92 winter season. This supports some conclusions of a recent investigation by Dunkerton and Delisi on the effects of the 1982 eruption of El Chichón and is shown to have a substantial effect on trace species transport, comparable to that associated to the Quasi-Biennal Oscillation (QBO).     

A large coaxial lidar for elastic and inelastic scattering studies of the stratosphere

A very powerful, fixed zenith laser radar system, based upon a liquid dye laser, is described in detail. The laser is a linear flashlamp pumped dye system capable of an average power output of 6W at a prf of 3 s^–1 using rhodamine 6G dye and about 3W with the same prf using esculin monohydrate dye. The use of this laser as a transmitting source in the lidar system has permitted the detection of Stokes shifted vibration/rotation Raman back-scattered light from neutral nitrogen in the stratosphere and from these measurements, stratospheric temperature profiles have been derived. The lidar has also been used to study elastic backscattering from the stratosphere at two wavelengths and preliminary results of a new technique for studying stratospheric aerosol are presented.     

Effect of the stratospheric medium temperature on strong lines observed in the far infrared region

Linearization of the radiative transfer equation is applied to the stratosphere. Calculations show a strong change in the water-vapour line shape on saturation. To interpret this phenomenon, we suggest an equivalent temperature of the stratospheric medium. The calculations are applied to terrestrial limb sighting using a spectrometer with a 0.001 cm^-1 resolution onboard a satellite.     

Recent stratospheric warmings at mid-latitudes and planetary-scale wave activity

Summary Recent stratospheric warmings are herewith examined in association with the influence of long-period planetary waves. Investigations on height and time variations of geopotential and wind data obtained in the winters of 1980–87 at different tropo-stratospheric levels from messages of the radiosondes situated in the Central Europe take necessarily into account the propagation of quasi-stationary and free planetary waves which in proper conditions mutually interfere with each other. The definition and evolution of the different types of warming are connected with the structure of these waves and enable one to interpret more completely coupling phenomena among different regions of the neutal atmosphere during the most disturbed winters. Distinct spectral analysis methods permit to extract from geopotential data series prominent wave periods within 1÷3 weeks. Westerly winds in the stratosphere reach unusual amplitudes before a sudden stratospheric warming, and this amplification is seen to be in connection with transient normal modes which are able to interact even with the mean circulation of the mesosphere and lower thermosphere.     

次声波在平流层波导中传播的数值模拟

对次声波在大气中传播进行了建模。通过结合保色散关系空间差分格式和Runge-Kutta时间格式的数值方法,建立了次声波传播模型。应用该次声波模型,研究了在耗散的重力分层大气中次声波的平流层导行传播。数值模拟结果表明,当次声波波包在平流层高度上被反射时,反射区域存在焦散现象,在声波的声压下降的同时,声波的能量得到聚焦。通过数值模拟结果与射线计算结果的对比表明,大气中声波传播的轨迹的精确描述需要应用全波解。