中纬电离层的低频等离子体不稳定性

本文研究了中纬电离层Rayleigh-Taylor不稳定性和漂移不稳定性的基本性质.如果不存在中性风和背景电场,中纬Rayleigh-Taylor不稳定性的增长率比赤道相应条件下的小.中性风和电场对中纬Rayleigh-Taylor不稳定性有重要影响.当中性风速度达到10m/s时,其作用就超过了重力.漂移不稳定性也受到重力、中性风和电场的影响,离子-中性粒子碰撞能降低漂移不稳定性的增长率.本文的分析表明,在赤道电离层容易产生大尺度的Rayleigh-Taylor不稳定性,在中纬电离层易产生小尺度的漂移不稳定性.     

带电粒子在中性线磁场中的运动

本文利用了小扰动方法、轨道法以及粒子的运动区域,比较系统地研究了带电粒子在中性线磁场中的运动。其结果是: 1.带电粒子的运动轨道可分为漂移轨道、波动轨道与8字形轨道三种形式。小扰动方法不能使用于中性线的区域内,在这个区域出现一个与小扰动漂移运动相反方向的运动。 2.在沿中性线方向的电场的作用下,中性线周围的部分粒子可以聚集到中性线附近。当粒子进入非小扰动区时,它们将被中性线磁场反射,并被电场加速。 3.计算出磁尾中性片的厚度为在这个区域内,大部分带正电的粒子的平均运动是沿晨昏方向,带负电的粒子的运动则相反。磁尾区出现一个等离子体中性片电流。     

赤道电离层R-T不稳定性发展的控制因素分析

本文从质量和电荷守恒方程出发,分析了控制电离层等离子体R-T不稳定性线性增长的各种因素,重点研究了热层风和背景电场的空间梯度对R-T不稳定性线性增长的影响.结果表明,热层风场和背景电场两者的空间梯度对R-T不稳定性线性增长有不可忽视的促进或抑制作用;对R-T不稳定性线性增长起促进作用还是起抑制作用,依赖热层风场和背景电场及其空间梯度的方向;对R-T不稳定性线性增长影响的显著程度主要依赖于热层风场和背景电场两者空间梯度的大小.数值计算结果表明,对典型的背景电离层条件,磁力线顶点高度为330 km时,对线性增长率的影响最高达到120%.     

行星际-赤道电场穿透效率的数值模拟研究

电场穿透作为强烈地磁活动期间太阳风-磁层-电离层电动力学耦合的重要形式,近年来一直是空间物理学最为活跃的课题之一.在过去的近40年里人们一直关注于观测证据的寻找以及物理过程的模拟,“电场穿透效率”这一重要定量化指标直到最近几年才得到重视.本文定义赤道电场增量与对应的行星际电场增量的比值为穿透效率,通过数值模拟的方法对穿透效率的地方时变化以及与中性风发电机的关系做了实验性讨论,在一定简化条件下结果表明：（1）不考虑跨极盖电势饱和的情况下,赤道电场增量与行星际电场增量呈线性关系,且中性风发电机并不影响电场穿透效率;（2）恒定重联线长度L的限定下,穿透效率具有明显地方时依赖性.例如对于L=2.6 R_E,在9LT至23LT之间,穿透效率维持在10％左右;0LT至7LT之间,穿透效率迅速从2％上升至30％后又迅速回落到原始水平,形成尖峰.这些结果基本符合观测特征.     

Global transport and localized layering of metallic ions in the upper atmosphere

A numerical model has been developed which is capable of simulating all phases of the life cycle of metallic ions, and results are described and interpreted herein for the typical case of Fe⁺ ions. This cycle begins with the initial deposition of metallics through meteor ablation and sputtering, followed by conversion of neutral Fe atoms to ions through photoionization and charge exchange with ambient ions. Global transport arising from daytime electric fields and poleward/downward diffusion along geomagnetic field lines, localized transport and layer formation through descending convergent nulls in the thermospheric wind field, and finally annihilation by chemical neutralization and compound formation are treated. The model thus sheds new light on the interdependencies of the physical and chemical processes affecting atmospheric metallics. Model output analysis confirms the dominant role of both global and local transport to the ions life cycle, showing that upward forcing from the equatorial electric field is critical to global movement, and that diurnal and semidiurnal tidal winds are responsible for the formation of dense ion layers in the 90–250 km height region. It is demonstrated that the assumed combination of sources, chemical sinks, and transport mechanisms actually produces F-region densities and E-region layer densities similar to those observed. The model also shows that zonal and meridional winds and electric fields each play distinct roles in local transport, whereas the ion distribution is relatively insensitive to reasonable variations in meteoric deposition and chemical reaction rates.     

On the possibility of ion-drag to induce dynamic instability in the lower thermosphere neutral gas

Strong wind shears may result in dynamic instability, often characterised by the Richardson number lying between zero and 0.25. The extent to which electric-field driven ion flow may induce such neutral wind shears is examined. Further, it is proposed that, in the ionosphere, it is possible for electric fields to drive ion winds such that the collisionally induced neutral air response may be comparable to viscous damping of neutral motion. We shall present an analogy to the Reynolds Number Re to quantify this effect. In the same way that Re may be used to evaluate the likelihood of a flow being turbulent, the analagous metric may also indicate where in the atmosphere plasma dynamics may be strong enough to destabilise the neutral dynamics.     

Longitudinal effects of ionospheric responses to substorms at middle and lower latitudes: a case study

An ionospheric model is used to simulate total electron content (TEC) disturbance events observed at middle and lower latitude sites near 75○W and 7○E longitudes. Within this longitudinal range, daytime TEC disturbances show patterns that are correlated with substrom activity seen in both auroral electrojet and ring current behavior. In modeling studies of the observed ionospheric effects, both electric field and neutral wind perturbations are examined as possible mechanisms. The morphological features of the required electric field perturbations near drawn and dusk are compared with those at other times to examine the local time characteristics of magnetospheric influence. Large-scale traveling atmospheric disturbances (TADs), an alternative candidate for the disturbance source, are also characterized and compared with known thermospheric behavior.     

Global aspects of solar wind–ionosphere interactions

Recent observations have quantified the auroral wind O⁺ outflow in response to magnetospheric inputs to the ionosphere, notably Poynting energy flux and precipitating electron density. For moderate to high activity periods, ionospheric O⁺ is observed to become a significant or dominant component of plasma pressure in the inner plasma sheet and ring current regions. Using a global circulation model of magnetospheric fields and its imposed ionospheric boundary conditions, we evaluate the global ionospheric plasma response to local magnetospheric conditions imposed by the simulation and evaluate magnetospheric circulation of solar wind H⁺, polar wind H⁺, and auroral wind O⁺. We launch and track the motions of millions of test particles in the global fields, launched at randomly distributed positions and times. Each particle is launched with a flux weighting and perpendicular and parallel energies randomly selected from defined thermal ranges appropriate to the launch point. One sequence is driven by a two-hour period of southward interplanetary magnetic field for average solar wind intensity. A second is driven by a 2-h period of enhanced solar wind dynamic pressure for average interplanetary field. We find that the simulated ionospheric O⁺ becomes a significant plasma pressure component in the inner plasma sheet and outer ring current region, particularly when the solar wind is intense or its magnetic field is southward directed. We infer that the reported empirical scalings of auroral wind O⁺ outflows are consistent with a substantial pressure contribution to the inner plasma sheet and plasma source surrounding the ring current. This result violates the common assumption that the ionospheric load is entirely confined to the F layer, and shows that the ionosphere is often an important dynamic element throughout the magnetosphere during moderate to large solar wind disturbances.     

热层风场的理论模拟与分析

通过求解中性大气Navier Stokes动量方程建立了一个时变的三维风场理论模式，利用目前新版的中性大气模式NRLMSISE 00及国际电离层参考模式IRI2000作为输入参数给出热层风场. 基于该模式，计算得到中等太阳活动年磁静日风场的变化形态及其受电场和离子曳力的影响. 同时，将Navier Stokes动量方程作不同形式的简化，并利用简化模式与本文的模式计算结果的对比，分析中性大气Navier Stokes动量方程中黏性项以及非线性项(U·Δ)U的作用. 结果表明，本文所建立伪三维风场模式给出的结果更为合理，而简化模式在某些地区尤其在低纬和赤道区不适用，黏性项及非线性项的作用不可忽略. 本文所建立的风场模式将对研究电离层动力学过程、电离层与热层的耦合过程以及空间天气学研究都有着重要意义.     

Morphology of slowly-varying geomagnetic external fields — a review

The morphology of slowly-varying geomagnetic external fields (such as Sq, L, and Dst) and the production mechanisms of electric current systems for those fields are briefly reviewed to provide background knowledge of the present state of research regarding the source fields for scientists concerned with the electromagnetic induction within the earth. It is concluded that both the Sq and Dst fields seem to have sources both in the magnetosphere and the ionosphere, while the classic idea of L-field production solely in the ionosphere by the wind dynamo is still acceptable.     

电场产生的赤道扩展F的时空演变

利用计算机模拟研究了非均匀电场产生的赤道扩展Ｆ的时空演变。非均匀电场能在赤道电离层Ｆ区底部触发Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性，导致等离子体泡形态。计算中非均匀电场的幅度取为０．２５－１．００ｍＶ／ｍ，所产生的等离子体泡在２０００ｓ以内就能穿过Ｆ峰达到５４０ｋｍ的高度。所得结果阐明了非均匀电场在赤道扩展Ｆ中的作用，指出了产生等离子体泡一种可能的扰动源。由于Ｅ区和Ｆ区电场是相互影响的，从而揭示了Ｅ区和Ｆ区扰动的相互联系。     

磁暴主相的經度不对称性

太阳带电粒子流通过磁边界层的中性点或阻塞点,由于极化电場的漂移运动或雷萊-泰勒（Rayleish-Taylor）不稳定运动,穿透磁层而进入地磁場被輻射带所捕获。这些被捕获的电子将向东漂移而形成西向电流。在它們围繞地球一周形成封閉的电流环以前,由于地球的自轉以及这个电流弧对地球各地的相对位置的不同,电流分布以及它在地面磁赤道附近所产生的磁場分布都将是不均匀的,对經度来耕是不对称的。本文假定电子繞地球的漂移周期为两天,計算了当磁暴开始时位于地方时00,06,12及18小时的台站記录的磁場随“暴时”T_st变化的曲线。理論結果与用佘山及洪伽幼地磁台站得到的統計結果是一致的。     

Evidence for nonlinear coupling of planetary waves and tides in the lower thermosphere over Bulgaria

Analyses of hourly values of zonal and meridional wind near 95 km observed by meteor radar at Yambol (42.5°N, 26.6°E) during January 1991–June 1992 indicate the presence of planetary waves with prevailing periods of 1.5–2.5, 4–6, 9–10 and 16–18 days. About 20% of the whole power of atmospheric motions is connected with these waves, so they play an important role in the dynamics of the mesosphere-lower thermosphere (MLT) region. By dynamic spectral analysis applied to the hourly neutral wind and to the calculated hourly values of tidal amplitudes it has been demonstrated that there is considerable modulation of tidal amplitudes by planetary waves in the neutral wind, as this process is better expressed in the semidiurnal tides. The nonlinear interaction between tides and planetary waves is studied by bispectral analysis. The results of these analyses indicate again that the nonlinear interactions between semidiurnal tides and planetary waves with periods 2–20 days are stronger than those of the diurnal tides and planetary waves. A peculiar feature of dynamics in the MLT region above Bulgaria is the presence of strong oscillations with periods of 20 and 30 h, which indicate significant nonlinear coupling between them.     

Analytical configurations of a force-free magnetic cylinder in the solar wind

The problem of a smooth field configuration, which should be an initial configuration in modeling (using the method of coarse particles) the problem of a stationary solar wind flow around a magnetic cloud in the case of a spatially two-dimensional statement (when a magnetic cloud is considered as a force-free magnetic cylinder with a finite radius) is considered. It has been indicated that such a statement is possible only when the magnetic field in the solar wind is parallel to the cylinder axis. The method for finding the magnetic field of a force-free cylinder with a finite radius, when some field component is specified and another component is determined based on this one (which makes it possible to construct fields with preassigned properties), has been proposed. The variant for constructing the initial field configuration in the transition region around a cylinder has been proposed. This variant makes it possible to gradually pass from homogeneous crossed fields in the solar wind to a force-free magnetic and zero electric fields within a cylinder, an electric field being potential and orthogonal to a magnetic field (in the reference system related to a magnetic cloud).     

Electric field effects on ionospheric and thermospheric parameters above the EISCAT station for summer conditions

Numerical calculations of the thermospheric and ionospheric parameters above EISCAT are presented for quiet geomagnetic conditions in summer. The Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) was used. The numerical results were obtained both with a self-consistent calculation of the electric fields of magnetospheric and dynamo-action origin and with the magnetospheric electric fields only. It was found that the dynamo-electric field has some effect on the ionospheric convection pattern during quiet geomagnetic conditions. It has a marked effect mainly on the zonal neutral wind component above EISCAT (±20m/s at 140 km altitude). We have studied the effects of various field-aligned current (FAC) distributions on thermosphere/ionosphere parameters and we show that a qualitative agreement can be obtained with region-I and -II FAC zones at 75° and 65° geomagnetic latitude, respectively. The maximum FAC intensities have been assumed at 03–21 MLT for both regions with peak values of 2.5 × 10^–7 Am^–2 (region I) and 1.25 × 10^–7 A m^–2 (region II). These results are in agreement with statistical potential distribution and FAC models constructed by use of EISCAT data. The lack of decreased electron density in the night-time sector as observed by the EISCAT radar was found to be due to the spatial distribution of ionospheric convection resulting from electric fields of magnetospheric origin.     

关于暴时电离层电流分布的南北半球不对称性

采用国际上广泛认可的高层大气和电离层经验模式提供的各种参数, 通过电离层电流连续方程, 计算出强磁暴条件下6月至日和12月至日内, 磁纬±72°和磁地方时00:00～24:00之间电离层电场、电流等的分布. 计算中考虑了地磁和地理坐标间的偏离; 除中性风场感生的发电机效应外, 还包含了磁层耦合(极盖区边界的晨昏电场和二区场向电流)的驱动外源. 结果表明, 6月至日时, 磁层扰动自极光区向中低纬的穿透情况在南、北半球内基本接近, 北半球内略强; 但12月至日时, 呈现明显的不对称性, 南半球的电流穿透远强于北半球, 而电场的穿透则是在北半球更强. 无论南北半球, 在中高纬地区, 午夜至黎明时段出现较强的东向电场分量, 其[WTHX]E×B[WTBZ]的向上漂移效应, 正是解释我们以往不少研究现象中所期盼的物理机制.     

Ground Motion Modeling for Seismic Hazard Analysis in the Near-source Regime: An Asperity Model

—A new, yet simple, method using the asperity model to estimate ground motion in the near-source regime for probabilistic seismic hazard analyses is proposed in this study. This near-source model differs from conventional empirical attenuation equations. It correlates peak ground motions with the local contributing source in terms of the static stress drop released non-uniformly on the causative fault plane rather than with the whole seismic source in terms of magnitude. Here the model is simplified such that ground motions at a rock or firm soil site near extended vertical strike-slip faults are dominated by direct shear waves. The proposed model is tested by comparing its predictions with strong ground motion observations from the 1979 Imperial Valley and the 1984 Morgan Hill earthquakes. The results have revealed that ground motions in the near-source region can be adequately predicted using the asperity model with appropriate calibration factors. The directivity effect of ground motion in the near-source region is negligible for high-frequency accelerations. The cut-off frequency (?_max?) at a site is an important parameter in the near-source region. Higher values of ?_max yield higher estimates of peak ground accelerations. For high-frequency structures, ?_max should be carefully estimated. In the near- source region both non-uniform and uniform source models can produce non-stationary high-frequency ground motions. Peak motions may not be caused by the nearest sections of the fault (even if the uniform source model is considered).     

中性风与等离子体行星波的同步振荡

本文是“赤道异常的行星波振荡”的续篇。提供新的事实以说明:由于中层大气的行星波振荡会调制到潮汐风上面,潮汐风上升到E层发电机区域,通过发电机效应产生电场及S_q电流体系的行星波振荡。电场沿磁力线传到F_2层,通过“喷泉效应”引起电离层赤道异常的行星波振荡。所有这些振荡来源相同,所以是同步的。     

Wavelet analysis on transient behaviour of tidal amplitude fluctuations observed by meteor radar in the lower thermosphere above Bulgaria

On the basis of bispectral analysis applied to the hourly data set of neutral wind measured by meteor radar in the MLT region above Bulgaria it was demonstrated that nonlinear processes are frequently and regularly acting in the mesopause region. They contribute significantly to the short-term tidal variability and are apparently responsible for the observed complicated behavior of the tidal characteristics. A Morlet wavelet transform is proposed as a technique for studying nonstationary signals. By simulated data it was revealed that the Morlet wavelet transform is especially convenient for analyzing signals with: (1) a wide range of dominant frequencies which are localized in different time intervals; (2) amplitude and frequency modulated spectral components, and (3) singular, wave-like events, observed in the neutral wind of the MLT region and connected mainly with large-scale disturbances propagated from below. By applying a Morlet wavelet transform to the hourly values of the amplitudes of diurnal and semidiurnal tides the basic oscillations with periods of planetary waves (1.5/20 days), as well as their development in time, are obtained. A cross-wavelet analysis is used to clarify the relation between the tidal and mean neutral wind variability. The results of bispectral analysis indicate which planetary waves participated in the nonlinear coupling with the atmospheric tides, while the results of cross-wavelet analysis outline their time intervals if these interactions are local.