地球等离子体层顶与地磁活动的关系研究

本文利用IMAGE卫星EUV相机观测的等离子体层图像,并采用最小L算法反演磁赤道面等离子体层顶位置.文中选取了2000-2002年间的3579幅等离子体层图像,并反演得到了时间间隔为1 h的等离子体层顶位置数据库,包含48899个等离子体层顶位置数据.利用该数据库统计研究了等离子体层顶位形随地磁活动的变化特性.统计发现等离子体层顶高度依赖地磁活动,与地磁指数Kp、Dst和AE均呈负相关,且等离子体层顶随地磁指数的变化趋势具有显著的MLT分布特性;亚暴活动对等离子体层顶演化的贡献在不同地磁活动期间有所不同,磁暴期间亚暴活动的贡献小,而地磁平静期亚暴的贡献大.本文研究工作为后续建立等离子体层顶模型和了解等离子体层顶的动态结构提供了重要研究基础.     

基于THEMIS卫星数据的等离子体层顶位置模型

利用2010~2011年的THEMIS卫星观测的卫星电位数据反演电子密度,并用来确定等离子层顶的位置,研究其与地磁指数的关系.利用观测数据得到了一个THEMIS卫星等离子体层顶位置数据库.利用该数据库建立了新的等离子体层顶位置的磁地方时经验模型,并比较了其与现有的四个等离子体层顶位置经验模型的优劣.发现新建立的等离子体层顶位置模型明显优于其余四个经验模型.在模型的三种输入参数中,Kp指数和Dst指数模型的预测效果相当,并明显优于AE指数的模型,这说明磁层的大尺度对流强度是影响等离子体层顶位置最主要的因素.     

EISCAT雷达观测到的极区电离层等离子体反常对流

综合分析EISCAT雷达与卫星当地测量数据,并利用磁层磁场模式对磁力线进行追踪,研究了发生在极光椭圆朝极盖边界附近电离层中,一例反常的背离太阳流动的强等离子体对流事件,及相关的太阳风-磁层-电离层耦合过程.结果表明,磁暴期间IMFBz指向南时观测到这一反常高速对流,及其相应的等离子体性态特征,很可能是向阳侧磁层顶磁重联过程在电离层中的印记.     

磁暴期间亚极光区极化流（SAPS）的DMSP观测与RAM模拟的比较研究

亚极光区极化流（Subauroral Polarization Streams, SAPS）为快速流动的西向等离子体流,位于昏侧-子夜前亚极光区,是磁层-电离层-热层耦合的重要过程之一.本文利用密西根大学的RAM （Ring current-Atmosphere Interaction Model）模型对一次典型磁暴期间发生的SAPS事件进行了模拟,并与DMSP卫星观测值进行了比较.结果表明：模拟结果能大致反映观测现象;模拟得到的SAPS峰值速度所在纬度随磁暴时间的变化与观测值有较大差别;SAPS速度观测值在约18∶00 UT和约20∶00 UT左右出现两个峰值,而模拟值只有一个峰值,出现在约18∶00 UT,主要原因是模型对亚暴过程的模拟存在不足.     

等离子体片离子与太阳风及地磁条件的关联研究

本文根据搭载于Cluster卫星的CIS/CODIF和RAPID仪器的观测数据,统计研究了等离子体片中的H⁺、O⁺离子在磁暴期间的时间变化特性,及其对太阳风条件的响应.观测结果表明：(1) 磁暴开始前,O⁺离子（0~40 keV）数密度保持在较低水平.随着磁暴的发展,O⁺数密度缓慢上升,其峰值出现在D_st极小值附近;H⁺离子(0~40 keV)数密度在磁暴开始之前的较短时间迅速增加并达到峰值,在磁暴开始之后迅速降低,并在整个主相和恢复相期间保持在相对较低水平.更高能量的离子则在磁暴开始后迅速增多,并在低能O⁺离子达到峰值之前达到峰值.因此我们推测磁暴初期从等离子体片注入环电流的主要是H⁺离子,主相后期O⁺离子可能扮演更为重要的角色.(2)在地磁活动时期,太阳风密度和动压强与等离子体片中的H⁺、O⁺数密度存在一定相关性.等离子体片中的H⁺离子对北向IMF B_z较为敏感,而IMF B_z南向条件下更有利于太阳风参数对等离子体片中O⁺数密度的影响.在地磁活动平静期,太阳风条件对等离子体片中的离子没有明显影响.     

亚暴期间等离子体对流速度及热层纬向风速度变化的统计学分析

本文利用2001年至2005年间CHAMP卫星及DMSP(F13,F15)卫星观测数据,对亚暴发生前后,亚暴初始位置所在磁地方时(Magnetic Local Time,MLT)东侧2h(+2h,MLT)至西侧4h(-4h,MLT)范围内等离子体对流速度(Vy)及热层纬向风速度(Uy)的变化进行了统计学分析.研究发现在亚暴发生后的1.5h内,所有MLT区间的Vy均明显增大,且峰值位置的地磁纬度向赤道侧移动,1.5h后,Vy减小,峰值的纬度向极区移动,表明亚暴的发生能显著增强等离子体对流速度;在亚暴发生位置的西侧0~2h内,Vy增幅最大,这表明亚暴对热层的影响主要在西侧,影响最大的区域是西侧0~2hMLT区间;Uy在亚暴发生后3h持续增大,其对亚暴的响应相较于Vy有1.5h的延迟.     

利用IMAGE卫星观测数据重建地球等离子体层的磁赤道面分布

美国利用IMAGE卫星的极紫外辐射（EUV）探测器对地球等离子体层进行了连续5年的遥感成像观测。由于IMAGE卫星数据是沿观测路径上的积分投影数据,并且存在地球“遮挡”、“阴影”、“数据缺失”等问题,无法直接利用传统的CT方法对等离子体层进行三维重建。本文利用地球磁场模型,基于地球等离子体层的物理性质,建立一个联系地球磁赤道面密度与投影数据的EUV成像模型,实现了从单个角度的EUV观测图像进行地球等离子层三维重构的方法。      

基于VAP卫星的等离子体层顶电子密度波动统计分析

基于范艾伦辐射带探测卫星的观测数据(2012年9月至2015年11月),收集了584个等离子体层顶密度波动事件,研究了这些事件分布随磁地方时、磁壳值以及地磁活动的变化关系,并使用快速傅里叶变换分析了全部事件平均功率谱.统计结果表明,等离子体层顶密度波动事件主要发生在磁地方时黄昏扇区,其分布与磁地方时和地磁活动具有明显的相关性.等离子体层顶密度波动在1～100mHz区间内具有接近-5/3的功率谱斜率,表明存在二维磁流体动力学湍流.本文统计结果将有助于进一步深入理解等离子体层顶密度波动在内磁层中能量传输的具体作用,并且将促进对等离子体层顶波动的激发、增强与传播过程的进一步研究.     

地磁扰动期间等离子体层顶结构的模拟研究

本文选取2001年6月8-10日的一个亚暴事件,模拟了在这期间等离子体层的结构演化过程.选取Weimer(2001模式)电场和Tsyganenko(1996模式)磁场作为背景电磁场,基于E×B的漂移运动计算磁赤道面内的带电粒子分布,模拟磁扰期间的等离子体层变化.模拟了等离子体层顶的结构和形状,结果有羽状、肩状和通道状结构,与同一时间点的EUV/IMAGE探测结果一致.     

GPS测量误差对大气掩星反演精度影响分析

在无线电GPS大气掩星观测中,GPS和LEO卫星的位置误差、速度误差、卫星钟的漂移、载波相位观测误差等会给其反演的大气参量带来一定的误差.为定量分析上述GPS测量误差对大气掩星反演精度的影响,本文采用模拟分析的方法,利用CHAMP实测轨道数据和大气与电离层经验模式,采用三维射线追踪方法模拟得到大气附加相位观测数据;然后根据各种GPS测量误差的不同性质,在模拟掩星观测大气附加相位数据中加入GPS测量误差并进行反演,将其反演结果与未加入误差时的反演结果进行比较,分析各测量误差所带来的掩星反演误差的特点;利用一天的151个分布全球的掩星事件的误差结果进行统计分析.结果表明:(1)±10 m的GPS径向位置误差引起的温度反演误差的平均误差在30 km的高度为(-+)0.8K,标准偏差为0.05K;±0.4 m的LEO径向位置误差引起的温度反演误差的平均误差在30 km的高度为(-+)1.2K,标准偏差为0.07K.(2)仅考虑卫星速度误差对掩星几何关系的影响,则±30 m/s的GPS速度误差引起的温度反演误差的平均误差最大为±0.02K,标准偏差为0.3K;±10 m/s的LEO速度误差引起的温度反演误差的平均误差最大为±0.4K,标准偏差为0.06K.(3)当卫星钟漂为1×10-9时,反演温度误差的平均误差在30 km的高度为263K,标准偏差为16K,当卫星钟漂为1×10-14时,反演温度误差的平均误差在30km的高度为0.019K,标准偏差为0.004K.(4)当L1和L2载波相位观测误差分别为1 mm和2 mm时,引起的温度反演的平均误差在30 km的高度为0.01K,标准偏差为0.4K.(5)利用较接近真实观测的三维大气和电离层背景下的模拟附加相位数据,同时加入国际最新指标的各种GPS测量误差进行反演,误差分析结果为:在30 km以下,折射率和密度的相对误差小于0.3%(1σ),压强的相对误差小于0.5%(1σ),温度误差小于1K(1σ).     

Interplanetary control of thermospheric densities during large magnetic storms

During the main phase of large magnetic storms significant energy can be deposited in the ionosphere but produce no commensurate magnetic perturbations on the ground. Consequently, models designed to predict and specify thermospheric energy budgets based on ground magnetic data are negatively impacted. To quantify these effects we compare thermospheric densities predicted by the MSIS model with those inferred from accelerometer measurements by the Gravity Recovery and Climate Experiment (GRACE) satellites during two magnetic storm periods in 2004. Although predictions and measurements are in substantial agreement during quiet times, the model significantly underpredicts densities during storms. Also, the model's maxima occur several hours after observed stormtime peaks. We show that polar cap potentials and magnetospheric electric fields derived from interplanetary parameters measured by the Advanced Composition Explorer satellite are roughly proportional to neutral densities observed by GRACE with lead times of ∼4 h. Finally, ion drift meter data from Defense Meteorological Satellite Program spacecraft suggest that unpredicted positive and negative spikes found in high latitude accelerometer data reflect encounters with strong head and tail thermospheric winds driven by anti-sunward convecting plasma.     

A new plasmapause location model based on THEMIS observations

A new empirical model of plasmapause location as functions of magnetic local time and geomagnetic indices has been developed based on the observations from THEMIS mission. We use the two-year data of electron density inferred from spacecraft potential to identify the plasmapause crossings and create a database of plasmapause locations. The database is further used to build up an empirical model of plasmapause related to magnetic local time based on the equation from O’Brien and Moldwin (2003). The new model is compared with previous plasmapause location models. It is found that our newly developed model is the best in predicting plasmapause locations among the existing models. The models based on Kp and Dst indices are better than the model based on AE index, suggesting that the plasmapause location is controlled by large scale convection of the magnetosphere.     

Indicators of an increase in the flux of relativistic electrons in geostationary orbit during geomagnetic storms

The relation of the fluxes of relativistic electrons in geostationary orbit during magnetic storms to the state of the magnetosphere and variations in the solar wind parameters is studied based on the GOES satellite data (1996–2000). It has been established that, in ～52–65% of all storms, the fluxes of electrons with energies higher than 0.6 and 2 MeV during the storm recovery phase are more than twice as high as the electron fluxes before a storm. It has been indicated that the probability of such cases is closely related to the prestorm level of fluxes and to a decrease in fluxes during the storm main phase. It has been found that the solar wind velocity on the day of the storm main phase and the geomagnetic activity indices at the beginning of the storm recovery phase are also among the best indicators of occurrence of storms with increased fluxes at the storm recovery phase.     

哨声模波对高能电子槽区和外辐射带的调节作用

本文利用磁层哨声模嘶声和合声波的幅度分布模型、近赤道背景电子(能量在eV量级)的数密度分布模型和IGRF10磁场模型建立了一个高能电子(能量大于50 keV)准线性扩散模型.模型的数值结果表明,在不同的地磁条件下,等离子体层顶位置的变化改变了磁层背景电子数密度的空间分布,从而改变了哨声模嘶声对高能电子有效的投掷角扩散（损失）区域,同时也改变了哨声模合声波对高能电子有效的动量扩散（加速）区域.哨声模嘶声对电子投掷角扩散区域的变化和RRES卫星探测到的高能电子的槽区变化是一致的,而合声波对电子的动量扩散区域的变化和卫星探测到外辐射带的变化相同.这种对应关系说明：在不同的地磁条件下,哨声模波对高能电子扩散区域的变化是造成高能电子槽区和外辐射带的空间位置和大小变化的一个重要因素.在一些强磁暴期间,由于嘶声对部分能量范围电子的投掷角扩散作用消失,这些电子的槽区也随之消失,从而使内外辐射带连接在一起.     

Geomagnetic storm intensity forecast caused by magnetic clouds of solar wind

Method of short-term forecast intensity of geomagnetic storms, expected by effect Solar wind magnetic clouds in the Earth’s magnetosphere is developed. The method is based calculation of the magnetic field clouds distribution, suitable to the Earth, the initial satellite measurements therein components of the interplanetary magnetic field in the solar ecliptic coordinate system. Conclusion about the magnetic storm intensity is expected on the basis of analysis of the dynamics of the reduced magnetic field Bz component clouds and established communication intensity of geomagnetic storms on Dst-index values and Bz component of the interplanetary magnetic field vector.     

Equatorial ionization anomaly and thermospheric meridional winds during two major storms over Brazilian low latitudes

The equatorial ionosphere responses over Brazil to two intense magnetic storms that occurred during 2001 are investigated. The equatorial ionization anomaly (EIA) and variations in the zonal electric field and meridional winds at different storms phases are studied using data collected by digisondes and GPS receivers. The difference between the F layer peak density (foF2) at an equatorial and a low latitude sites was used to quantify the EIA; while the difference between the true heights (h_F) at the equatorial and an off-equatorial site was used to calculate the magnetic meridional winds. The vertical drift was calculated as dh_F/dt. The results show prompt penetration electric fields causing unusual early morning development of the EIA, and disturbed dynamo electric field producing significant modification in the F region parameters. Variations to different degrees in the vertical drift, the thermospheric meridional winds and the EIA developments were observed depending on the storm phases.     

Studying ionospheric responses to magnetic storms depending on the local time of the storm main phase onset

A morphological analysis of vertical sounding data obtained in Irkutsk from 2003 to 2008 has been performed. The AE index was used to determine the geomagnetic activity level, and the storm main phase onset was registered based on the D _st index. The ionospheric response to a magnetic storm was estimated based on the relative deviation of the critical frequency and altitude of the ionospheric F2 region from the median values. Superstrong magnetic storms and storms without positive initial phases were not considered when the data were selected. We found that positive ionospheric disturbances, which were accompanied by an increase in the F2 region maximum altitude, predominated between the storm initial phase and main phases during all considered magnetic storms. Between these storm phases, negative disturbances were only registered at night. Predominance of positive ionospheric disturbances over negative ones can be related to the selection of storms for studies.     

磁层顶压缩事件的磁场分析

向日面磁层顶在平静太阳风条件下，处于10RE（RE为地球半径）左右．但在异常的太阳风条件下，即南向行星际磁场很强和（或）太阳风的动压很大时，会被压缩，甚至到达同步轨道附近．集中分析2001年4月11日的磁暴事件，研究当磁层顶发生强烈压缩以后。在地球空间和地面上产生的磁场影响．磁层顶位形选取Shue（1998）模型计算．当把计算结果与GOESl0卫星的观测数据对比时发现：磁层顶在强的太阳风条件下的确会被压缩到同步轨道以内．Shue（1998）模型的预测基本正确，通常的漏报可能是由于预报的位置误差所致．实际磁层顶电流片的位置和强度与我们假设的理想磁层顶间断面计算结果基本吻合．在分析大磁暴过程时，磁层顶压缩使磁层顶电流对于中低纬度地磁场扰动有突出的贡献，在2001年4月事件中，这个贡献可以大于50nT，占主相的1／6左右．这一贡献可以使Dst指数产生相应的误差．     

Ion drift in the earth’s inner plasmasphere during magnetospheric disturbances and proton temperature dynamics

Based on the thermal plasma measurements in the Earth’s inner plasmasphere on the INTER-BALL-2 and MAGION-5 satellites it has been indicated that the plasmaspheric ion temperature as a rule decreases during the main phase of magnetic storms; in this case the plasma density increases or remains at the level typical of undisturbed conditions. The physical mechanism by which the ion drift during a magnetic storm results in a temperature decrease is described. It is shown that the third adiabatic invariant also remains in processes with a characteristic time shorter than the period of charged particle drift around the Earth for cold equatorial plasma. The constructed model of the drift shell displacement from the Earth caused by a decrease in the magnetic field in the inner magnetosphere during the development of a magnetic storm satisfactorily describes the decrease in the proton temperature near the equatorial plane.     

Experimental evaluation of the effect of storm movement on peak discharge

The hypothesis that downstream moving storms with storm length less than watershed length(L_s/L< 1.0) magnify the peak discharges as indicated by kinematic-wave models in previous studies was evaluated in an analysis of the dimensionless peak discharge and dimensionless storm velocity.Previously unpublished experimental data collected for a V-shaped watershed in the Watershed Experimentation System(WES) at the University of Illinois at Urbana-Champaign,were used in comparison with the simulation results of a kinematic-wave model.It is found that downstream moving storms with L_s/L < 1.0 increase the peak discharges to a limited extent compared to stationary storms,and the kinematic-wave model overstates the increase in the peak flows resulting from downstream moving storms with L_s/L < 1.0.To evaluate the importance of the backwater effects in the experimental watershed,the accuracy of kinematic-wave and dynamic-wave models for the simulation of surface runoff resulting from upstream and downstream moving storms was evaluated utilizing the same experimental data.The kinematic-wave model simulates the upstream moving storms pretty well,i.e.Nash-Sutcliffe coefficient of model fit efficiency equal to 0.948 and 0.831 for storms lengths equal to and not equal to the watershed length,respectively.Whereas,the kinematic wave model substantially overestimates the peak discharge of downstream moving storms,and yields generally poorer fits than for upstream moving storm,i.e.NSE equal to 0.867 and 0.674 for storms with lengths equal to and not equal to the watershed length,respectively.The dynamic-wave model simulates the downstream moving storms pretty well,i.e.NSE equal to 0.843 and 0.879 for storms with lengths equal to and not equal to the watershed length,respectively,indicating backwater significantly affects runoff for even this simple experimental watershed.Considering that storm movement did not substantially magnify peak discharge,the assumption of stationary storms made in standard hydrologic design seems reasonable and adequate.