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

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

非相干雷达探测低电离层的结果分析

位于波多黎各的Arecibo非相干雷达可以获得低电离层电子和离子密度, 利用此非相干雷达数据对中纬度低电离层的运动特征进行研究. 得到了电子密度随时间和高度的变化 情况, 结果显著呈现出周日变化特征, 并分析了电子密度随高度的变化规律. 进一步对数据进行频谱分析, 深入研究低电离层电子密度的周日变化效应. 得到电子密度的高度剖面, 发现从F层底部到E层有明显的等离子体沉降. 低电离层的层结构特征及电子密度变化表明, 在该区域还存在不同程度的等离子体扰动, 由此对低电离层的作用因素 进行分析, 认为大气潮汐或声重波可能对低电离层产生扰动, 即低电离层与大气存在一定程度的耦合作用.      

非线性大气重力波产生的行进电离层扰动

本文导出了大气重力波的非线性方程,研究了重力波的非线性性质。利用等离子体的连续性方程,计算了电离层对重力波的非线性响应,讨论了非线性电离层扰动(TID)的形态。结果表明,非线牲效应使TID表现为多种不规则形态。      

声重波激发的赤道 Spread—F 不稳定性

本文导出了声重波扰动与赤道电离层扩展F层相互作用的基本方程,并对数值模拟结果进行了讨论。结果表明当声重波扰动与Kavleigh-Taylor不稳定性共同作用于赤道F层时,将导致大尺度波状结构的Spread-F.其特征为以一定速度漂移的准周期不规则性,发展到后期伴有羽柱形状的bubble(电子密度低区)产生。这种现状的bubble可上升至F层顶部并发生倾斜。声重波波矢量的方向决定着羽柱的倾斜方向。这些结果与实验观测特征相符,说明声重波扰动不仅是大尺度Saread-F的触发源,而且对其整个形态发展过程起着重要的作用。      

强震前异常电场激发声重波对电离层的影响

震前地震孕育期地表异常增强的电场,通过大气电导率传输到电离层高度.该异常电场通过非稳态局部加热,可以在电离层高度激发声重波.基于该理论,利用一维时变中纬电离层物理模型,模拟了该扰动源对电离层电子密度的影响.结果表明,重力波引起的中性风速度扰动对电离层电子密度分布影响甚微,该机理无法解释震前电离层异常扰动现象.      

电离层F层中舒曼共振现象的检测

利用Ａｕｒｅｏｌ－３卫星在电离层Ｆ层高度上（６５０ｋｍ左右）对电离层电场扰动与电子密度扰动的观测进行分析处理,在高纬电离层Ｆ层高度上可以检测到舒曼共振现象。根据等离子体参数,可以认为,发生在这个高度上的舒曼共振现象与Ｌ模式的波传播有关,舒曼共振的基频是与电子密度的大尺度不规则性有关,谐波分量与正向密度梯度相关。     

中层大气行星波的跨赤道传播

使用一个全球原始方程半谱模式,模拟了行星波的跨赤道传播。结果指出:冬半球对流层顶附近准定常行星波的幅度起伏和相位变化,可以造成中层大气行星波的跨赤道传播,并且在夏半球形成瞬变行星波。所得结果为中层大气行星波的跨赤道传播、夏半球行星波的形成和南北半球之间的耦合提供了一种物理机制。      

中高层大气行星波与惯性重力内波的非线性相互作用

在弱非线性理论基础上,将三维大气中行星波和惯性重力波从原始非弹性近似方程中分离出来,讨论了典型的2天行星波与惯性重力内波的非线性相互作用过程.从共振曲面和参量不稳定增长率来看,行星波倾向于与空间尺度较大的惯性重力波发生相互作用.利用潮汐波的等价重力波假设,讨论了2天行星波与半日潮及9.6h惯性重力波的相互作用,三波相互作用时能量守恒.非线性相互作用使2天波和潮汐波的波幅受到长期调制.      

大气波动监测仪在岩石圈-大气层-电离层研究中的应用

针对北京大学空间物理与应用技术研究所研制的两台大气波动监测仪近三年的观测数据,对这一时段内所观测到的重力波和次声波周期尺度扰动的形态特征、谱结构特征及其在时间分布上的统计特征进行了分析. 给出了几例雷暴、地震等事件中观测到的大气扰动,并揭示了这些事件期间观测到的与地面大气扰动周期尺度相类似的电离层扰动. 结果表明,北京大学研制的大气波动监测仪可有效记录到地面的微弱大气扰动,观测数据可用于进行岩石圈-大气层-电离层之间的耦合研究.      

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

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

The ionospheric effect of atmospheric gravity waves excited prior to strong earthquake

We have computed perturbations in the nighttime mid-latitude F2 region ionosphere that could be produced by internal atmospheric gravity waves generated before strong earthquakes through ionospheric Joule heating due to the seismogenic electric field of short duration. There is a strong anisotropy of the atmospheric gravity wave effect with respect to the imminent earthquake epicentre, the electron density changes being maximum poleward and equatorward of the epicentre and being minimum eastward and westward of it. It should be noted that the duration of the electron density perturbation in the F2 region ionosphere is much longer than the duration of the primary precursor of an earthquake – the enhancement of the vertical electric field at the Earth’s surface, which initiates the atmospheric gravity wave generation. This fact is important from the practical point of view of predicting catastrophic earthquakes.     

Mathematical modeling of plasma drifts over equatorial low latitude regions

This paper presents a mathematical model to simulate ionospheric plasma drifts at equatorial low latitude regions by coupling of E- and F-regions. The governing non-linear differential equations (of elliptic and parabolic nature) are solved numerically through finite-difference schemes and obtained neutral winds and electric fields. The temperature and electron density profiles are generated utilizing MSIS-86 atmospheric model. The continuity equation is employed to obtain night-time E-region density profile using measured ionograms at Trivandrum (India). The computed vertical and zonal plasma drifts are comparable with measured Jacamarca plasma drifts with little variations during noon and evening times. The plasma drifts at Trivandrum (8.5° N, 76.5° E, dip 0.5° N) are compared with those of Jicamarca (12° S, 76.9° W, dip 2° N). Neutral wind simulations of present model agree well with those of horizontal wind model (HWM-93). The post-sunset enhancement and its reversal are also discussed.     

Theoretical simulation of electron density and temperature distribution at Indian equatorial and low latitude ionosphere

The electron density and temperature distribution of the equatorial and low latitude ionosphere in the Indian sector has been investigated by simultaneously solving the continuity, momentum and energy balance equations of ion and electron flux along geomagnetic field lines from the Northern to the Southern hemisphere. Model algorithm is presented and results are compared with the electron density and electron temperature measured in situ by Indian SROSS C2 satellite at an altitude of ∼500 km within 31°S–34°N and 75 ± 10°E that covers the Indian sector during a period of low solar activity. Equatorial Ionization Anomaly (EIA) observed in electron density, morning and afternoon enhancements, equatorial trough in electron temperature have been simulated by the model within reasonable limits of accuracy besides reproducing other normal diurnal features of density and temperature.     

论赤道扩展F的密度不规则体(Ⅱ)

本文对等离子体中的电子和离子分别采用流体和动力论描述, 导出适于描述夜间赤道F层现象的若干非线性模型。由此证明赤道扩展F小尺度区域可存在大量偶极涡、孤子和孤子链等类型的密度不规则体。      

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

Effect of the seed perturbation amplitude on the equatorial spread F initiation during solar minimum

The contribution of gravity wave (GW) to the initiation/development of spread F during a solar minimum year was investigated through the comparison of the observed precursory parameters and characteristics of the corresponding equatorial spread F (ESF) events. The ionospheric parameters were recorded at the magnetic equatorial station Sao Luis (2.3°S, 44°W, dip latitude 2°S) during March and October 2010. These data were used to estimate the influence of the relative gravity wave amplitude and the ambient ionospheric condition on the diurnal variation of the spread F initiation. The vertical velocity drift indicated a clear control and defines the threshold for the seasonal variability of the ESF occurrence. However, it was insufficient to solely determine or predict the day to day variation of ESF occurrence. Thus, few days with contrasting ambient ionospheric condition and magnitude of GW amplitude were analysed in order to investigate the role of the different precursory factors in the observed diurnal variation of the plasma irregularity development. The density scale length and gravity wave amplitude were shown to immensely contribute to the linear instability growth rate, especially during the days with a low post-sunset rise. Thus, the experimental observations have demonstrated the strong inter-dependence between the precursory factors and they have also highlighted the probable control of the ESF morphology.     

Signatures of strong geomagnetic storms in the equatorial latitude

Ionosonde data from two equatorial stations in the African sector have been used to study the signatures of four strong geomagnetic storms on the height – electron density profiles of the equatorial ionosphere with the objective of investigating the effects and extent of the effects on the three layers of the equatorial ionosphere. The results showed that strong geomagnetic storms produced effects of varying degrees on the three layers of the ionosphere. Effect of strong geomagnetic storms on the lower layers of the equatorial ionosphere can be significant when compared with effect at the F2-layer. Fluctuations in the height of ionization within the E-layer were as much as 0% to +20.7% compared to −12.5% to +8.3% for the F2-layer. The 2007 version of the International Reference Ionosphere, IRI-07 storm-time model reproduced responses at the E-layer but overestimated the observed storm profiles for the F1- and F2-layers.