电离层暴时经验模型STORM在中国区域的适应性研究

利用中国区域内9个垂测站1976---1987年一个太阳活动周期的电离层暴时f0F2数据, 统计分析了电离层暴事件的等级, 以及不同等级的电离层暴随季节和地磁纬度的分布特征. 研究发现, 中小型电离层暴在春秋季发生的概率较大, 不同季节的发生次数与地磁纬度具有明显的关系. 利用STORM模型对电离层暴时f0F2和大型及特大型电离层暴时f0F2的预测值与月中值进行了比较. 结果表明, 除了冬季误差增大外, 发生电离层暴时STORM模型能够有效地改善月中值模型. 增加中国的暴时数据, 并提高对冬季的暴时参数f0F2的预测是改善STORM模型的重要因素. 建立合适的暴时指数来预测f0F2是未来研究的重点.      

正相电离层骚扰研究

本文对我国五个电离层站月报资料中的正相电离层骚扰(f₀F₂增加)作了分析研究,获得了正相骚扰时空变化的规律性,地磁暴与正相电离层骚扰的关系.对正相电离层骚扰的机理也作了讨论.      

基于电离层暴时f_0F_2经验模型Kalman滤波短期预报

利用时间累积地磁指数印ap(T),建立了强地磁扰动条件下电离层f0F2与月中值相对偏差经验模型.该经验模型只在春秋季节和夏季特强地磁扰动条件ap(T)>100,即时间累积地磁指数大于100时达到理想精度.尝试利用气象预报中常用的Kalman滤波方法对模型的系数进行实时修正,以提高预报精度,并对长春站1986-1995年近一个太阳周f0F2数据进行提前一小时预报试验.冬季预报均方根误差为0.76MHz,春秋季节为0.68MHz,夏季为0.61MHz.在特强地磁扰动条件下,预测误差在0.87～1.43MHz之间.该预报方法同时与包含暴时修正模型STORM的国际参考电离层IRI2001进行了比较,展示了Kalman滤波方法实时修正模型系数的能力和良好的应用前景.     

中低纬电离层不同季节地磁暴响应的事例分析

利用中国中低纬台站漠河（53.5°N，122.3°E）、北京（40.3°N，116.2°E）、武汉（30.5°N，114.2°E）和三亚（18.3°N，109.6°E）的电离层观测数据，对比分析了4个台站电离层参数在2015年不同季节4个地磁扰动事件期间的变化特征.结果表明，4个磁暴事件期间电离层的响应特征并不完全一致，有着明显的季节特征，春季、夏季和秋季电离层以负相扰动为主，冬季以正相扰动为主.分析发现，中性成分O/N₂的降低与电离层负相扰动有关，但三亚地区的负相扰动还与扰动发电机电场相关.正相扰动的机制在不同事件中并不相同，穿透电场可能是引起春季磁暴事件期间电离层短时正暴效应的原因，而冬季长时间的正暴效应则是扰动电场和中性风共同作用的结果.      

由环电流的宇宙线效应验证二维对称环电流模型

本文利用不同纬度超中子堆记录的宇宙线强度, 对不同磁暴类型下环电流的宇宙线效应进行了分析。加强的西向环电流使宇宙线截止刚度下降, 在某些特大的磁暴中, 中低纬字宙线截止刚度下降可达1GV以上, 从而使宇宙线强度增加。利用二维对称环电流模式, 从宇宙线的角度计算了环电流的磁暴效应。与实际观测值的对比显示, 除非特大的磁暴过程, 利用该模式所得结果与实际观测符合得较好。      

冬夏典型地磁-电离层暴期间太平洋区域F层行为

本文在前人已公布的暴时电离层动态资料的基础上, 进行了太平洋地区冬至与夏至期间电子总量和F层特征的分析研究, 取得了新的结果.认为在磁暴主相出现之时, 以电子总量为表征的电离层暴达到最明显的程度.以F层为代表的电离层的行为强烈地受控于经度(暴时)、纬度、地磁格局和季节(南北半球)的缔合关系;并且电子总浓度的消长关系需从电离增长率(磁壳层收缩)-电离消矢率(损失系数加大)-电离迁移率(磁共轭迁移和电离峰谷异动)作统一考虑.文中提出了较合理的电离层暴区划意见并讨论了暴情倒相问题.      

磁暴期间中低纬电离层暴形态和物理机制分析

本文利用东亚地区１２个低纬电离层台站的测高仪观测数据，对１９７８年８月２７日发生的一次曲型磁暴期间电离层峰值高度和密度的变化进行了分析。采用滑动平均区分开电离层中不同时间尺度的扰动，分析了影响中低纬度电离层暴的几种扰动形态特征，并对其物理机制进行了讨论。结果表明：伴随磁暴急始的磁层压缩，电离层中表现出峰值密度增加和峰值高度下降；磁暴主相期间热层大气暴环流及其所引起的中性大气成分变化控制着电离层的大     

基于深度学习递归神经网络的电离层总电子含量经验预报模型

利用行星际太阳风参数与太阳活动指数、地磁活动指数、电离层总电子含量格点化地图数据,首次基于一种能处理时间序列的深度学习递归神经网络（Recurrent Neural Network,RNN）,建立提前24h的单站电离层TEC预报模型.对北京站（40°N,115°E）的预测结果显示,RNN对扰动电离层的预测误差低于反向传播神经网络（Back Propagation Neural Network,BPNN）0.49～1.46TECU,将太阳风参数加入预报因子模型后对电离层正暴预测准确率的提升可达16.8%.RNN对2001和2015年31个强电离层暴预报的均方根误差比BPNN低0.2TECU,将太阳风参数加入RNN模型可使31个事件的平均预报误差降低0.36～0.47TECU.研究结果表明深度递归神经网络比BPNN更适用于电离层TEC的短期预报,且在预报因子中加入太阳风数据对电离层正暴的预报效果有明显改善.      

中纬电离层暴负相开始时间与磁暴主相开始时间的对应关系

本文假设在磁暴主相期间,由于极光椭圆带处的空气被加热上升,从而使高纬高空出现富含分子的气体。这些气体由于扩散及与中性风的相互作用会向中低纬移动,其所到之处电子消失系数增加,从而导致负相电离层暴的发生。计算给出了全球中纬电离层暴负相的开始时间与磁暴主相开始时间之间的关系,并讨论了负相电离层暴发生的"时间禁区"问题。结果与有关统计结果符合得很好。      

中纬电离层爆负相开始时间与磁暴主相开始时间的对应关系

本文假设以磁暴主相期间，由于极光椭圆带处的空气被加热上升，从而使高纬高空出现富含分子的气体，这些气体由于扩散及与中性风的相互作用会向低纬移动，其所到之处电子消失系数增加，从而导致负相电离层暴的发生，计算给出了全球中纬电离层暴负相的开始时间与磁暴主相开始时间之间的关系，并讨论了负相电离层暴发生的“时间禁区”问题，结果与有关统计结果符合得很好。     

宇宙线强度变化与磁扰K类型

本文把1966—1983年期间发生的679个地磁暴进行了分类,利用统计方法分析了各类磁暴发生前后宇宙线强度的变化特征.突发急始脉冲发生后,宇宙线强度没有出现显著的变化;缓始型暴发生后,宇宙线强度出现Forbush下降,但下降幅度很小;急始型暴发生后,宇宙线强度出现十分明显的Forbush下降.当把急始型暴按K指数大小和持续时间分为5种类型,发现它们伴随的宇宙线Forbush下降是不一样的,其下降幅度随磁暴的增强而加大,下降的速率随磁暴的增强而加快,扰动的持续时间随磁暴的减弱而增加.      

宇宙线变化特征在强地磁暴预报中的应用

利用宇宙线中子探测数据定性分析了地面宇宙线多台站之间的相互联系以及大磁暴与宇宙线之间的响应关系. 以Irkutsk和Oulu宇宙线台站为例, 运用小波去噪技术提高数据的稳定性. 结果表明, 相同世界时条件下, 两站宇宙线通量相关性在事件发生时较高; 而相同地方时条件下, 相关性则在平静期较高. 进一步采用相同地方时条件对不同宇宙线台站的通量在平静期和扰动期的相对变化进行分析, 选取2004年7月强地磁暴典型事例进行直观分析, 发现大地磁暴前Irkutsk和Oulu台站的宇宙线相对通量发生明显差异, 可以尝试作为强地磁暴宇宙线先兆特征. 通过对2001年3月至2005年5月的强磁暴和中强磁暴进行统计, 得到与强地磁暴相关的适当宇宙线相对差异阈值. 将得到的阈值对2005年9月至2011年12月所有强磁暴及中强磁暴进行验证, 总成功率达到87.5%, 误报率为35.7%, 结果较好.      

利用地面宇宙线强度变化预报地磁暴方法初步研究

分析了Nagoya宇宙线闪烁体望远镜探测数据的变化特点,定性地探讨了CME可能引起的地面宇宙线的变化特征,通过实例证实了地面宇宙线通量的异常波动是地磁暴发生的重要先兆特征,并且将8 h内宇宙线通量与该时间段内平均通量的偏差D8(t)参数应用到宇宙线数据分析中.通过数据分析与讨论,认为D8(t)参数达到一定阈值是地磁暴的重要先兆特征,但不是充分条件,虚假信号仍占多数;D8(t)参数与太阳质子事件探测结果相结合,对于大地磁暴的预报有较好的效果.      

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

Cosmic ray exposure factors for Shuttle altitudes derived from calculated cut-off rigidities.

The allowed cosmic radiation flux accessible to an earth-orbiting spacecraft is a complex function of the satellite position and the geomagnetic cutoff characteristics at each zenith and azimuth angle at each position. We have determined cosmic ray exposure factors for the galactic cosmic ray spectrum for typical shuttle altitudes and inclinations up to 50 degrees. We have utilized d world grid of trajectory-derived cutoff rigidity calculations at 400 km altitude to determine geomagnetic transmission functions that permit a simple and direct calculation of the allowed cosmic ray spectrum to a 400 km satellite orbit. If the interplanetary cosmic ray spectrum is multiplied by the orbit-averaged geomagnetic transmission function the result is the allowed cosmic ray spectrum at the spacecraft.     

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

Seasonal variations of storm-time TEC at European middle latitudes

Global Navigation Satellite System (GNSS) measurements of the Total Electron Content (TEC) from local (Dourbes, 50.1°N, 04.6°E) and European IGS (International GNSS Service) stations were used to obtain the TEC changes during the geomagnetic storms of the latest solar activity cycle. A common epoch analysis, with respect to geomagnetic storm intensity, season, and latitude, was performed on data representing nearly 300 storm events. In general, the storm-time behaviour of TEC shows clear positive and negative phases, relative to the non-storm (median) behaviour, with amplitudes that tend to increase during more intense storms. The most pronounced positive phase is observed during winter, while the strongest and yet shortest negative phase is detected during equinox. Average storm-time patterns in the TEC behaviour are deduced for potential use in ionosphere prediction services.     

Effect of Forbush decreases of cosmic ray flux on ozone at higher middle latitudes

The Forbush decreases of cosmic ray flux occur prevailingly together with geomagnetic storms, because these phenomena have a similar origin in solar/interplanetary processes. To study the effects of large Forbush decreases on total ozone at middle latitudes, we use the TOMS total ozone data along latitudinal circles 40°N and 50°N. The effects of Forbush decreases are found to occur or to be non-measurable under the same conditions as those of geomagnetic storms: certain effect occurs only at 50°N (not 40°N), in winter, under conditions of high solar activity and the east phase of the QBO. However, the effects of the analyzed Forbush decreases are weaker than the effects of strong geomagnetic storms.