热源激发重力波特征以及波流作用的数值模拟研究

本文在二维等温可压大气中引入了一个随时间和空间变化的热源扰动,分别以静止风和中纬1月份月平均向东的纬向风急流为背景,对不同背景下热源激发的重力波的传播详细过程及其特性进行了数值模拟研究.热源激发出来的重力波在初始阶段有很宽的频谱范围,随后由于重力波的传播效应,水平波长和垂直波长分布范围随时间都有所减小.顺风传播的重力波...     

中国地区下平流层惯性重力波参数分布特征的资料分析

重力波参数特征的确定是全球大气环流模式中平流层重力波参数化的一个重要环节,这些参数的选取需要根据观测资料来确定.本文根据中国地区南北向8个站点(自北向南依次为嫩江、锡林浩特、北京、郑州、南阳、宜昌、怀化和南宁,纬度范围为22°~49°N)两年多的垂直高分辨率无线电探空资料,分析了下平流层(18~25km)惯性重力波活动的季节变化及其随纬度的分布特征.主要结果有:(1)与全球其他地区一样,中国地区下平流层重力波能量具有明显的季节变化和经向分布特征:冬季大夏季小,且随着纬度增加而减小;纬向和经向速度扰动大小一致,表明重力波能量在水平方向上是各向同性的.(2)重力波垂直波长随纬度没有明显变化,集中在1~3km,占所有样本的85%以上,平均值约为2.0km;水平波长80%以上集中在100~800km,平均值约为450km,有随纬度增加而降低的趋势(南北水平波长相差达40km左右),水平波长与垂直波长之比大约为200:1,表明下平流层重力波基本上沿水平方向传播,垂直方向的夹角极小.(3)固有频率和科氏参数之比有随纬度增加而减小的特征,集中在1~2,平均值约为1.5.(4)重力波能量主要向上传播,各站点重力波上传频率均在60%以上;水平传播方向有比较显著的方向性,主要是东西方向传播且与盛行风向有关,但是该各向异性随着纬度的增加有所减弱.     

新疆库尔勒地区秋季下平流层重力波特性

利用2011年秋季无线电探空数据,采用矢端曲线法首次分析了新疆库尔勒地区下平流层重力波特征参量,得到36组准单色重力波的结果.结果统计显示：库尔勒秋季下平流层重力波垂直波长、水平波长平均值分别为2.8 km和580 km,固有频率平均值为1.74f.垂直传播方向以上传为主,约占78%,其中下传重力波水平波长较短,固有频率较高.水平传播方向以西北和东南为主,各占1/3,其中上传（下传）重力波水平传播方向以西北（东南）居多,这与热带低纬站点和其他中纬站点观测结果不同.与其他站点比较,库尔勒地区ŵ/f最小,中高纬地区水平波长、垂直波长随纬度增加大致有减小的趋势,库尔勒地区偏离这一趋势,波长偏大.     

一次暴雨过程中重力波参数演变特征的模拟结果

本文利用中尺度数值模式WRF对2003年7月4-5日淮河流域特大暴雨过程进行了数值模拟,并利用高时空分辨率模拟结果资料,提取了暴雨中心区大气重力波频率、周期、水平波长、垂直波长、水平相速和群速等特征参数,分析了暴雨过程中重力波参数随时间的演变特征.结果表明,对此次暴雨强降水过程影响较大的重力波主要是发展的中 α 尺度波和中 β 尺度波,暴雨后期随着重力波的频散,周期和水平波长有减小趋势,频率有增大趋势.非降水区的重力波参数特征和降水中心区有明显不同,大气中小振幅的中 α 尺度和中 β 尺度重力波是否发展和暴雨强降水的发生关系密切.     

中层大气中重力波饱和机制的数值分析

通过对向上传播的重力波波包在中层大气中的非线性传播过程进行数值模拟,讨论中层大气中重力波的饱和机制.数值模拟结果表明,向上传播的重力波波包的扰动振幅在接近波包的本征水平相速度之前随高度单调增长,而当波振幅接近本征水平相速度时,在对流不稳定区域出现等位温面的翻转,同时波振幅的增长达到饱和（波振幅随高度不发生变化）.小尺度对流在等位温面的翻转和波饱和发生后产生,随后波包开始破碎,这些非线性过程的最终结果产生湍流.表明导致重力波饱和的关键因素是等位温面的翻转而不是诸如波破碎、湍流、波-流相互作用等其他的一些物理过程.     

利用瑞利激光雷达和无线电探空仪观测数据对武汉上空重力波特性的研究

本文利用武汉大学的瑞利激光雷达的瑞利散射回波数据(30～65km)来研究武汉地区上空(30.5°N,114.4°E)重力波的活动规律和统计特性.通过对2003年12月到2005年3月观测的200h数据反演的密度进行处理分析,得到了重力波的一些个例特征,并提取垂直波长为2 km以上的重力波进行统计分析.结果表明,最可几的垂直波长是3～4 km和17～20 km,重力波振幅的月平均值在冬季有较大值,夏季值较小.与武汉2004年无线电探空仪的密度扰动提取到的重力波做比较,发现与瑞利激光雷达得到的重力波振幅的月平均值有很强的相关性,也是冬季值比较大,夏季值比较小.通过无线电探空仪的风场数据,本文还得到了急流的季节变化规律、最大风剪切年变化规律,发现急流和最大风剪切与激光雷达的重力波统计结果有很强的相关性.     

基于轨道扰动引力谱分析的方法确定低低观测重力卫星反演重力场的空间分辨率

基于低低卫-卫跟踪重力卫星的轨道特性,从垂直和水平两个方向计算了重力卫星高空扰动引力,并根据其谱特性及星载加速度的测量噪声水平分析了重力卫星能反演重力场的阶数.利用EGM96重力场模型分别计算了400 km、450 km和500 km 轨道高度处重力卫星受到的扰动引力谱及扰动引力谱的平均量级,分析其垂直特性表明:在三个轨道高度处能分别能反演150、140和130阶的重力场模型.利用两颗同轨重力卫星相距220 km的特性,计算了400 km、450 km和500 km 轨道高度处纬度相差2°的两颗卫星纬向扰动引力差,即扰动引力水平分量,分析其谱特性,表明:重力卫星能反演至117阶的地球重力场模型.     

基于激光雷达手段的海南地区重力波与其波谱的季节分布特性研究

利用子午工程海南激光雷达对我国海南地区上空进行持续观测,通过3年的累积观测数据对我国低纬度地区重力波活动的季节分布特性进行研究,依据重力波线性理论对海南地区上空的大气密度扰动规律、空间功率谱及时间频率谱进行分析,并通过选择波长在1km至8km范围内具有特定波长以及具有波动周期为60 min至25min的特定频率的重力波辅助研究大气密度扰动的季节变化规律,总结得出海南地区重力波活动具有夏季大、春秋季小、而冬季依然频繁的季节性分布规律.结合海南地区特殊的地理位置与当地季节性气候特征分析得出海南地区上空重力波活动季节性变化的可能原因为青藏高原地形及我国南海地区存在的热带强对流与赤道潜流共同作用的结果.     

高纬地区低层大气行星波的无线电探空仪观测研究

本文通过分析美国阿拉斯加地区三个站点(Nome,64.50°N,165.43°W;McGrath,62.97°N,155.62°W;Fairbanks,64.82°N,147.87°W)无线电探空仪1998～2006年观测数据,研究了北半球高纬地区低层大气行星波特性.通过分析发现行星波主要存在于两个区域,一个在对流层顶附近,一个在冬季极夜急流附近,两个区域的行星波都具有明显的间断性,持续时间一般不超过2个月;三个分量中,温度扰动量的振幅最小,经向风扰动量的振幅最大.对流层顶附近的行星波没有明显的季节变化且谱成分较为复杂,5天波的振幅最小,10天波的振幅略强于16天波.极夜急流附近的行星波主要出现在冬季,波振幅比对流层顶附近小,主要为10天波和16天波,且16天波的振幅强于10天波.由折射指数可以看出,夏季在对流层上方有明显反射层,冬季则较弱甚至消失,很好地解释了平流层行星波主要在冬季出现的原因.对2003/2004年冬季三个站点行星波的细致分析发现对流层区域和极夜急流区域出现明显的准10天波和准16天波,准10天在垂直方向为驻波,温度分量垂直波长约为12km,经向风分量垂直波长大于26km,波自东向西传,纬向...     

地球南北半球的非对称性

依据新的计算分析和空间观测数据,进一步论述了地球南北半球的非对称性. 全球热散失量的计算得出,南半球高出北半球33髎;南半球地幔热散失量是北半球的2倍. 比较南北半球S波速度分布,得出南半球的上地幔为低速、高温,北半球的上地幔为高速、低温. 计算地幔各层的质心位置发现,地球的质心偏于北半球. 计算地球经、纬圈长度的年变化率表明,南半球在扩张,北半球在收缩. 用空间大地测量数据的检测结果证实,南半球处于扩张状态,北半球处于压缩状态. 对地球的非对称性作了初步的动力学解释.     

Seasonal–latitudinal variations of PMC particle size from SME measurements for the northern 1983 season and SNOE measurements for the northern 2000 and southern 2000/2001 seasons

The ultraviolet spectrometers (UVS) on the solar mesosphere explorer (SME) and student nitric oxide explorer (SNOE) measured scattered limb radiance at small and large scattering angles from polar mesospheric clouds (PMCs). The SME data are for the northern summer hemisphere (NH) in 1983 and 1984. The SNOE measurements are for the NH in 2000 and for the southern hemisphere (SH) in the 2000/2001 season. From this database, we deduce the modal particle size from the measured scattering angle asymmetry in radiance. This quantity is determined as a function of time within the PMC season, and latitude, assuming several scattering models depending on the adopted size distribution and particle shape. For assumed spherical particles with a Gaussian distribution of width of 14 nm, the results for SME show mode particle sizes that vary from about 35 to 60 nm throughout the season. The results for SNOE under the same assumption show that for high latitudes in the NH the particle size grows systematically from the seasonal onset, from about 25 nm to a maximum of about 45 nm at 30 days after solstice. Lower latitudes show a similar time dependence, but with smaller particle sizes. SH PMC particle sizes display a more complicated seasonal variability. Generally, variability in measured cloud height is anti-correlated with particle size for the seasons analyzed here. Particle sizes in the SH are generally smaller than those in the NH, consistent with the northern bias in PMC brightness, and with previous satellite studies. These results are interpreted in terms of our understanding of PMC microphysics and inter-hemispheric differences in temperature and dynamics. Our quantitative results for mode radius depend on the assumption of a constant distribution width. If the width varies with latitude or time, our calculated gradients of mode radius would be different.     

Inertial gravity wave parameters for the lower stratosphere from radiosonde data over China

Characterization of gravity wave(GW)parameters for the stratosphere is critical for global atmospheric circulation models.These parameters are mainly determined from measurements.Here,we investigate variation in inertial GW activity with season and latitude in the lower stratosphere(18-25 km)over China,using radiosonde data with a high vertical resolution over a 2-year period.Eight radiosonde stations were selected across China,with a latitudinal range of 22°-49°N.Analyses show that the GW energy in the lower stratosphere over China has obvious seasonal variation and a meridional distribution,similar to other regions of the globe.The GW energy is highest in winter,and lowest in summer;it decreases with increasing latitude.Velocity perturbations with longitude and latitude are almost the same,indicating that GW energy is horizontally isotropic.Typically,85%of the vertical wavelength distribution is concentrated between elevations of 1 and 3 km,with a mean value of 2 km;it is almost constant with latitude.Over 80%of all the horizontal wavelengths occur in the range 100-800 km,with a mean value of 450km;they show a weak decrease with increasing latitude,yielding a difference of about 40 km over the 22°-49°N range.The ratio of horizontal wavelength over vertical wavelength is about 200:1,which implies that inertial GWs in the lower stratosphere propagate along nearly horizontal planes.Ratios of their intrinsic frequency to the Coriolis parameter decrease with increasing latitude;most values are between 1 and 2,with a mean value of 1.5.Study of the propagation directions of GW energy shows that upward fractions account for over 60%at all stations.In contrast,the horizontal propagation direction is significantly anisotropic,and is mainly along prevailing wind directions;this anisotropy weakens with increasing latitude.     

Interhemispheric comparison of mesospheric ice layers from the LIMA model

The new LIMA/ice model is used to study interhemispheric temperature differences at the summer upper mesosphere and their impact on the morphology of ice particle related phenomena such as noctilucent clouds (NLC), polar mesosphere clouds (PMC), and polar mesosphere summer echoes (PMSE). LIMA/ice nicely reproduces the mean characteristics of observed ice layers, for example their variation with season, altitude, and latitude. The southern hemisphere (SH) is slightly warmer compared to the NH but the difference is less than 3 K at NLC/PMC/PMSE altitudes and poleward of 70N/S. This is consistent with in situ temperature measurements by falling spheres performed at 69N and 68S. Earth's eccentricity leads to a SH mesosphere being warmer compared to the NH by up to approximately 85 km and fairly independent of latitude. In general, NH/SH temperature differences in LIMA increase with decreasing latitude and reach at 50. The latitudinal variation of NH/SH temperature differences is presumably caused by dynamical forcing and explains why PMSE are basically absent at midlatitudes in the SH whereas they are still rather common at similar colatitudes in the NH. The occurrence frequency and brightness of NLC and PMC are larger in the NH but the differences decrease with increasing latitude. Summer conditions in the SH terminate earlier compared to NH, leading to an earlier weakening and end of the ice layer season. The NLC altitude in the SH is slightly higher by 0.6–1 km, whereas the NLC altitudes itself depend on season in both hemispheres. Compared to other models LIMA/ice shows smaller interhemispheric temperature differences but still generates the observed NH/SH differences in ice layer characteristics. This emphasizes the importance of temperature controlling the existence and morphology of ice particles. Interhemispheric differences in NLC/PMC/PMSE characteristics deduced from LIMA/ice basically agree with observations from lidars, satellites, and radars.     

Simulated impacts of atmospheric gravity waves on the initiation and optical emissions of sprite halos in the mesosphere

The correlation between atmospheric gravity waves(GWs) and Transient Luminous Events(TLEs) has been poorly studied using both synchronous observations and numerical simulations. To investigate the modulation effects of GWs on TLEs,a troposphere-mesosphere quasi-electrostatic field model is developed in three-dimensional Cartesian coordinates, and the effects of GW perturbations on the initiation and optical emissions of sprite halos are simulated using the model. Simulation results indicate that the atmospheric density at lower ionosphere altitudes becomes inhomogeneous due to GW perturbations, and sprite halos tend to initiate in the GW troughs due to the lower electric breakdown threshold. GW perturbations cause the deformation of sprite halos, strong luminous regions distribute mainly along the GW troughs while optical intensities along the GW peaks is relatively weak. Larger GW perturbations lead to more pronounced deformation of sprite halos, however, stronger lightning discharges in the troposphere result in less optical perturbations of sprite halos. The observed luminous intensities and optical morphology of sprite halos are also affected by the observing orientations and the lightning polarities.     

与东北冷涡相伴的高空急流诱发平流层重力波的数值模拟研究

使用中尺度数值模式WRF-ARW,针对2010年6月发生在中国东北地区一例伴随对流层高空西风急流（位于~9 km高度）演变过程出现的平流层重力波活动特征开展了数值模拟. 事件发生期间,对流层区域环流处在一个东北冷涡系统的控制之下. 模拟结果再现了该东北冷涡的发展和维持过程,以及与之相伴的高空急流的特征. 模拟结果揭示出在急流区域上空的平流层中存在显著重力波活动现象. 分析结果显示,重力波活动与急流存在紧密联系,在水平方向上,重力波呈显著的二维结构,出现在急流出口区上部并逆背景流向西传播. 功率谱分析结果表明盛行波动具有~700 km水平尺度、9~12 h时间尺度以及4~5 km垂直波长. 由于急流的存在,造成其与平流层中下部之间存在显著的水平风速垂直切变,与切变相伴的耗散使得上传的重力波动量通量数值随着高度升高而递减. 同时,在18~20 km高度间出现的西风-东风转换带极大地抑制了波动在垂直方向的传播,形成显著动量通量沉积效应. 估算结果表明,在11~20 km高度之间,这种效应的整体作用相当于对该层背景流施加强度为0.86 m·s^-1·day^-1的动力阻曳.     

Interpretation of SOFIE PMC measurements: Cloud identification and derivation of mass density,particle shape,and particle size

The Solar Occultation For Ice Experiment (SOFIE) was launched onboard the Aeronomy of Ice in the Mesosphere (AIM) spacecraft to measure polar mesospheric clouds (PMCs) and their environment. This work describes methods for identifying PMCs in SOFIE observations and determining mass density, particle shape, particle effective radius, and the parameters of a Gaussian size distribution. Results using SOFIE measurements from the northern summer of 2007 are compared with concurrent observations by the ALOMAR lidar in northern Norway. Ice particle properties determined from SOFIE are in good agreement with the lidar results, considering the differences in instrument characteristics.     

Phase functions of polar mesospheric cloud ice as observed by the CIPS instrument on the AIM satellite

The Cloud Imaging and Particle Size (CIPS) instrument on the Aeronomy of the Ice in the Mesosphere (AIM) spacecraft is a 4-camera nadir pointed imager with a bandpass centered at 265 nm and a field of view of 120°×80°. CIPS observes polar mesospheric clouds (PMCs) against the sunlit Rayleigh-scattered background. At individual polar locations approximately 5 km×5 km in area, CIPS observes the same volume of air seven times over a range of scattering angles from about 35° to 150°. These multi-angle observations allow the identification and extraction of the PMC scattered radiance from the Rayleigh-scattered background. We utilize the fact that the former has a highly asymmetric phase function about 90° scattering angle, while the latter has a phase function that is symmetric. The retrieved PMC phase function can then be interpreted to obtain PMC particle size distributions. We describe a technique for identification of PMCs in the CIPS observations through the separation of the Rayleigh and PMC radiances. PMC phase function results are shown for the first season of CIPS observations. Assuming the particles are oblate spheroids with an axial ratio of 2, and a Gaussian distribution of width 14 nm, we find the phase functions are consistent with mean radii between 50 and 60 nm. These results are similar to those discussed by Hervig et al. [2009. Interpretation of SOFIE PMC measurements: cloud identification and derivation of mass density, particle shape, and particle size. J. Atmos. Sol. Terr. Phys., in review.] in this issue from the Solar Occultation for Ice Experiment (SOFIE) which also flies on the AIM satellite.     

The global distribution of gravity wave energy in the lower stratosphere derived from GPS data and gravity wave modelling: Attempt and challenges

Five years of global temperatures retrieved from radio occultations measured by Champ (Challenging Minisatellite Payload) and SAC-C (Satelite de Aplicaciones Cientificas-C) are analyzed for gravity waves (GWs). In order to separate GWs from other atmospheric variations, a high-pass filter was applied on the vertical profile. Resulting temperature fluctuations correspond to vertical wavelengths between 400 m (instrumental resolution) and 10 km (limit of the high-pass filter). The temperature fluctuations can be converted into GW potential energy, but for comparison with parameterization schemes GW momentum flux is required. We therefore used representative values for the vertical and horizontal wavelength to infer GW momentum flux from the GPS measurements. The vertical wavelength value is determined by high-pass filtering, the horizontal wavelength is adopted from a latitude-dependent climatology. The obtained momentum flux distributions agree well, both in global distribution and in absolute values, with simulations using the Warner and McIntyre parameterization (WM) scheme. However, discrepancies are found in the annual cycle. Online simulations, implementing the WM scheme in the mechanistic COMMA-LIM (Cologne Model of the Middle Atmosphere—Leipzig Institute for Meteorology) general circulation model (GCM), do not converge, demonstrating that a good representation of GWs in a GCM requires both a realistic launch distribution and an adequate representation of GW breaking and momentum transfer.