地球同步轨道高能电子通量预报方法研究

统计分析了GOES卫星测量得到的E > 2MeV能道电子通量与地磁Ap指数及太阳风数据的关系, 构建了基于径向基函数RBF的神经网络模型框架, 对GOES-12卫星所处的地球同步轨道高能电子通量进行提前1天的预报, 其对2008-2010年数据预测的效果较好. 另外, 发现在GOES-12卫星观测的E >2MeV能道高能电子达到108 cm-2·d-1·sr-1以上时, FY-2D卫星的测量数据同时达到108 cm-2·d-1·sr-1以上的比例达到90%左右. 通过对FY-2D卫星E >2MeV能道电子通量与GOES卫星E>2MeV电子通量的相关性分析, 建立了FY-2D卫星高能电子预报模型, 预报结果与实测通量符合较好.      

基于同步轨道区磁场对相对论电子通量的预报

地球同步轨道区相对论电子通量的变化与该区域的磁场分量P(垂直轨道面指向北)有很好的相关性, 后者的变化相对于前者有1~2天的时间提前量, 这为相对论电子通量的预报提供了可能. 通过对二者相关性物理机制的研究, 结合GOES11 (135oW)卫星数据分析, 确定最佳时间提前量和最优相关系数, 并提出了一个定量的预报模型. 模型的输出参量为24 h之后的当地时子夜(23:31 LT---00:30 LT)、清晨(05:31 LT---06:30 LT)、正午(11:31 LT---12:30 LT)、傍晚(17:31 LT---18:30 LT) 4个特征时间段相对论电子通量1 h平均值, 预报的相对论电子通量有>0.6 MeV和>2 MeV两个谱段, 预报精度0.7左右. 这种预报模式对地球同步轨道卫星的自主安全运行具有较好的应用价值.      

基于经验模态分解的地球同步轨道高能电子通量预报

在磁暴恢复相期间,大量相对论(高能)电子从磁层的外辐射带渗透到地球同步轨道区.其中> 2 MeV的高能电子能够穿透卫星表面并聚积在材料内部,导致卫星无法正常运行或完全损坏.磁暴期间的高能电子通量变化的非平稳与非线性特征十分明显.通过实验发现,经验模态分解法能够极大地降低高能电子通量非平稳性问题造成的预报影响.以2008-2009年的数据作为训练集,2010-2013年数据作为测试集.结果表明:2010-2013年的预报率约为0.84;在太阳活动较为复杂的2013年,预报率达到0.81.引入经验模态分解后预报效率得到显著提高.     

地球同步轨道高能电子增强事件预报方法

分析了地球同步轨道高能电子通量增强事件的发生规律及其与太阳风和行星际磁场参数的关系，并在此基础上建立了基于人工神经网络的高能电子增强事件模式，经实测数据检验，预报模式可以对未来1天的高能电子通量进行预报，误差为8．2％，达到了较高水平．     

基于支持向量机方法的地球同步轨道相对论电子事件预报模型研究

提出了一种基于支持向量机方法(SVM)的地球同步轨道相对论电子事件预报模型. 模型以平均影响值(MIV)作为指标, 筛选出预报输入参量. 这些参量包括, 前一日的大于2MeV电子日积分通量、太阳风速度、太阳风密度、Dst指数和前二日的AE指数. 模型包含回归和分类两个部分, 可以分别对未来一天的电子日积分通量和相对论电子事件强度的级别做出预报. 对2008年样本进行测试, 在相对论电子通量的预报中, 预报值和实测值之间的线性相关系数为0.85, 预报效率为0.71; 对相对论电子事件级别预报的准确率为82%, 可以较准确区分开事件状态与非事件状态. 结果表明, SVM预报模型对相对论电子事件有较好的预报效果.      

强磁暴、能量粒子暴与热层大气密度涨落之间的相关关系

利用1997-2007年由GOES8, GOES11和GOES12星载高能粒子探测器在地球同步轨道高度上所探测到的高能质子和高能电子通量探测数据以及高度560km左右星载大气密度探测器所得的热层大气密度探测数据, 统计分析了强地磁扰动、高能粒子通量跃变和热层大气密度涨落之间的相关关系, 初步获得强地磁扰动期间, 地球同步轨道(外辐射带外环)均出现了增幅大于三个数量级的高能质子通量(尤其是E>1MeV)强增强现象, 随后热 层大气密度强烈上涨, 表明三者之间是正相关关系. 在时间上地球同步轨道高能质子通量强增强现象先于日均Ap值(地磁活动程度)上涨约一天左右, 而热层大气密度强涨落现象又明显滞后于强地磁扰动事件.      

中地球轨道高能电子辐射环境特性分析

采用中国中地球轨道卫星在太阳活动下降相到上升相的高能电子探测数据, 首次分析研究了该轨道高能电子环境的空间分布、通量强度、时序变化以及对地磁暴活动响应的特性. 结果表明, 中地球轨道高能电子的空间分布 范围稳定, 电子通量强度随能量升高而下降; 中地球轨道高能电子环境是 一个在不同时间尺度上剧烈变化的动态系统, 该系统可能间歇性地出现27天重 现性变化, 该系统变化受地磁暴事件调制, 但其对磁暴的响应呈现出非线性特征.      

GEO轨道相对论电子日积分通量预报统计建模

基于辐射带相对论电子哨声波局地加速理论,将地磁AE指数作为源电子通量和通量各向异性的指标,将地磁Dst指数作为损失机制的指标,利用滑动窗口线性滤波器方法,建立了一个地球静止轨道大于2MeV相对论电子预报模型.利用该模型开展了2000-2009年地球静止轨道相对论电子通量预报试验.研究发现,这10年总预报效率为0.818,2003年的预报效率（0.633）最低,2009年的预报效率（0.856）最高.模型预报效果与持续模型相比有很大提高,略低于利用太阳风参数作为输入的同类预报模型的预报效果.这说明即使在缺少太阳风参数的情况下,该模型利用地磁扰动参数也能取得较好的预报效果.当模型输入参数增加了太阳风速度时,即综合考虑了行星际扰动和磁层扰动对辐射带粒子加速过程的影响,模型逐年的预报效率进一步提升.其中,2005年的预报效率提升了9.5%,这10年的总预报效率增加到0.848,预报值与实测值之间的线性相关系数为0.918,均方根误差为0.422.      

太阳质子事件中短期预报模型研究

太阳质子事件通量的预测对航天器抗辐射加固设计和航天员出舱活动具有非常重要的意义.针对一年以下的航天任务,利用经验统计方法,确认太阳活跃年和太阳平静年期间,1——365天不同时间段内 > 10MeV,> 30MeV和 > 60MeV的太阳质子事件积分通量符合对数正态分布,且通量对数的标准偏差σ和期望值μ随任务期时间的变化满足对数函数形式.以此为基础,构建太阳质子通量的中短期预报模型.该模型能够针对太阳活跃年和太阳平静年,给出一定置信度下1——365天不同时间内 > 10MeV,> 30MeV和 > 60MeV的质子事件通量分布,从而为执行中短期航天任务提供太阳质子事件通量的预测,以规避不必要的风险.      

使用地磁脉动参数定量预报地球同步轨道相对论电子通量的建模研究

磁层超低频波(ULF波)对种子电子的加速机制是磁层相对论电子产生的一个重要机制, 而地磁脉动参数可以作为此机制的有效指标. 本文采用地磁脉动参数作为输入参数, 借鉴线性预测滤波器技术, 构建一个多参量非线性函数, 进而利用此函数以及卡尔曼滤波技术, 建立一个地球同步轨道相对论电子通量日积分值预报模式, 提供提前一天的预报值. 使用2004年数据对该模式进行训练, 预报结果的预报效率为0.73, 线性相关系数为0.85. 使用2005-2006年的数据对该模式进行测试, 预报值与实测值之间的线性相关系数为0.83, 预报效率为0.69, 相比Persistence模式具有较大提升, 与仿REFM模式的预报效率相当      

Prediction of geosynchronous electron fluxes using an artificial neural network driven by solar wind parameters

There are hundreds of satellites operating at the geosynchronous (GEO) orbit where relativistic electrons can cause severe damage. Thus, predicting relativistic electron fluxes is significant for spacecraft safety. In this study, using GOES satellite data during 2011–2020, we propose two neural network models with two hidden layers to predict geosynchronous relativistic electron fluxes at two energy channels (>0.8 MeV and > 2 MeV). The number of input neurons of the two channels (>0.8 MeV and > 2 MeV) are determined to be 36 and 44, respectively. The > 0.8 MeV model has 22 and 9 neurons in the hidden layers, while the > 2 MeV model has 25 and 15 neurons in the hidden layers. The input parameters include the north–south component of the interplanetary magnetic field, solar wind speed, solar wind dynamic pressure and solar wind proton density. Through the analysis of different time delays, we determine that the optimal time delays of two energy channels (>0.8 MeV and > 2 MeV) are 8 days and 10 days, respectively. The training set and validation set (Jan 2011-Dec 2018) are divided by the 10-fold cross-validation method, and the remaining data (Jan 2019-Feb 2020) is used to analyze the model performance as a test set. The prediction results of both energy channels show good agreement with satellite observations indicated by low RMSE (～0.3 cm^-2sr^-1s^?1), high PE (～0.8) and CC (～0.9). These results suggest that only using solar wind parameters is capable of obtaining reasonable predictions of geosynchronous relativistic electron fluxes.     

Topology of the high latitude magnetosphere during large magnetic storms and the main mechanisms of relativistic electron acceleration

One of the main endeavors of the “Space Weather” program is the prediction of the appearance of very large fluxes of relativistic electrons with energies larger than 1 MeV, because they represent a serious potential hazard for satellite missions. Large fluxes of relativistic electrons are formed in the outer radiation belt during the recovery phase of some storms. The formation of large fluxes is connected to a balance between the acceleration and loss processes. A two-step acceleration process is ordinarily analyzed. A “Seed” population with energies ∼hundreds of keV appeared during expansion phase of magnetospheric substorm. A “Seed” population is additionally accelerated obtaining relativistic energies by some other process. Several acceleration mechanisms have been proposed for the explanation of the electron acceleration, including radial diffusion and internal acceleration by wave-particle interactions. Nevertheless, none of them takes into account great changes of magnetospheric topology during a magnetic storm. Such changes are mainly connected with asymmetric and symmetric ring current development. We analyze the changes of magnetospheric topology during magnetic storms. We show that a change of the magnetospheric magnetic field can be the important factor determining the acceleration of relativistic electrons.     

Qualitative estimation of magnetic storm efficiency in producing relativistic electron flux in the Earth’s outer radiation belt using geomagnetic pulsations data

It is well known that during many but not all of the geomagnetic storms enhanced fluxes of high-energy electrons are observed in the outer radiation belt. Here we examine relativistic (>2 MeV) electron fluxes measured by GOES at the synchronous orbit and on-ground observations of two types of ULF pulsations during 30 magnetic storms occurred during 1996–2000. To characterize the effectiveness of the chosen magnetic storms in producing relativistic electron fluxes, following to (Reeves, G.D., McAdams, K.L., Friedel, R.H.W., O’Brien, T.R. Acceleration and loss of relativistic electrons during geomagnetic storms. Geophys. Res. Lett. 30, doi:10.1029/2002GL016513, 2003), we calculate a ratio of the maximum daily-averaged electron flux measured during the recovery phase, to the mean pre-storm electron flux. A storm is considered an effective one if its ratio exceeds 2. We compare behavior of Pi1 and Pc5 geomagnetic pulsations during effective and non-effective storms and find a tendency for a storm efficiency to be higher when the mid-latitude Pi1 pulsations are observed for a long time during the magnetic storm main phase. We note also that the prolonged powerful Pc5 pulsation activity during the recovery phase of a magnetic storm is the necessary condition for the storm effectiveness. To interpret the found dependences, we suggest that there are two prerequisites for generating relativistic electron populations during a storm: (1) the availability of seed electrons in the magnetosphere, and Pi1 emissions are indicators of the mid-energy electron interaction with the ionosphere and (2) acceleration of the seed electrons to MeV energies, and interaction of electrons with the MHD wave activity in the Pc5 range is one of the most probable mechanisms proposed in the literature for this purpose.     

A global empirical orthogonal function model of plasmaspheric electron content

Based on the method for establishing a global plasmaspheric model using observations from COSMIC and MetOp-A orbit determination GNSS receivers, Chen et al. (2017) obtained a global plasmaspheric total electron content product with a spatial resolution of 2.5° × 5° and a time resolution of 4 h. In this paper, we use those global plasmaspheric electron content product in 2008, 2010, 2011, 2012, 2014 for 1446 days to establish a global plasmaspheric empirical model based on empirical orthogonal function (EOF). The model can well characterize the spatiotemporal variation of plasmaspheric electron content (PEC) and the influence of solar radiation on it. Only the first four orders of EOF sequences can characterize the 98.43% features of the original PEC dataset. The principal component coefficient Pk is decomposed twice during modeling, and the combination of trigonometric function and linear function is used to model Pk to characterize the solar cycle, annual cycle, semi-annual cycle and quarter-cycle variation. We compare the PEC model values with the actual observation data, the results show that the empirical PEC model values are highly correlated with the actual observations. The correlation between the two is above 0.96, and the RMS maximum of the difference between the PEC model values and the observed values are 0.70 TECU, and the average of the difference between the PEC model values and the observed values are −0.18 TECU, respectively. In addition, we validate the reliability of the global plasmaspheric model established by two empirical orthogonal function decomposition method using actual observation data, according to the global distribution of the differences between the PEC model values and the observed values in low solar activity and high solar activity, it can be seen that under low solar activity and high solar activity conditions, the model has good adaptability.     

A numerical model approximating extreme energetic electron events involved in the physical and chemical processes of the middle atmosphere

Relativistic electrons (with energies >150 keV) which originate in the outer radiation belt and detected by the Russian ‘Meteor’ series of satellites have been correlated with the atmospheric total ozone data compiled by almost 90 stations located around the world within the latitude zone 40°–70°N. In more than 60% of the stations examined we have detected a clear decrease of the ozone 3–5 days after the electron flux excess. A numerical model has been applied to approximate this effect based on relativistic electron initiated nitric oxides creation in the upper mesosphere with subsequent atmospheric transport (both vertical and horizontal) towards the upper stratosphere. A first attempt of local and temporal prediction of ozone depletion because of energetic electrons impact in the middle atmosphere has been illustrated.