基于SAMPEX卫星观测的南大西洋异常区高能质子动态分布特征

利用SAMPEX卫星1992年7月至2004年6月19～27MeV高能质子数据对南大西洋异常区的分布特征进行研究, 发现南大西洋异常区高能质子分布随高度及F_10.7的变化十分显著. 在540±25km高度上, 地磁较为平静时期南大西洋异常区高能质子微分通量随着F_10.7的增大而减小, 同时在F_10.7≥115sfu时减小趋势较为平缓. 对中等及以上磁暴进行统计分析发现, 磁暴期间南大西洋异常区高能质子微分通量和SYM-H指数的绝对值存在明显的反相关关系, 且地磁暴对南大西洋异常区高能质子微分通量存在明显的持续影响效应. 磁暴发生期间高能质子微分通量明显减少. 磁暴恢复相及其之后高能质子微分通量呈现较为显著的恢复过程.      

FY-3A卫星与NOOA系列卫星高能带电粒子实测结果的比较

FY-3A卫星是运行于830 km高度的太阳同步轨道气象卫星,其搭载的空间环境监测器可观测3～300 MeV的高能质子和0.15～5.70 MeV的高能电子.FY-3A卫星在轨工作期间,太阳活动处于由谷年向峰年过渡期,空间环境非常平静,探测结果显示3～300 MeV的高能质子分布主要集中在南大西洋辐射带异常区,0.15～5.70 MeV的高能电子分布区域除南大西洋异常区外,还分布在南北两极高纬区域.FY-3A与NOAA卫星测量结果反映出带电粒子强度及分布区域随投掷角变化的空间各向异性特征.本文在充分考虑了带电粒子时间、空间分布差异以及比对探测器之间自身设计差异的前提下,经过归一化处理后,首次对两颗卫星同期探测结果进行相关性分析,验证了两颗卫星相同时空条件下高能带电粒子通量分布的一致性;说明FY-3A空间环境监测器不仅具备空间带电粒子辐射监测能力,且探测结果有效可靠,可用作辐射带环境数据源的组成部分,为发展新的模型,深入研究辐射带高能粒子的分布、起源和传输等提供新的观测依据.     

质子辐射带辐射中心区域模型

利用AP-8和CRRESPRO质子辐射带模式数据库, 比较了二者在磁赤道面上计算结果 的差异并给出其差异原因. 利用AP-8模式数据库数据, 建立起质子辐射带函数形式的辐射中心区域模式, 包括各能道全向微分通量峰值所对应的L值(L_c) 随能量E的变化模式以及各能道全向微分通量峰值J_max随能量E的变化模式. 利用RBSP A卫星REPT望远镜在磁赤道面上的高能质子观测数据, 分别与 AP-8模式、CRRESPRO模式及本文所得辐射中心区域模式计算结果进行比较, 发现在78.9, 102.6和208MeV三个能道上, RBSP A卫星观测所得各能道全向微分 通量明显偏大, 而L_c与AP-8或本文辐射中心模式所得结果基本一致; RBSP A卫星也观测到CRRESPRO Quiet模式所得的隐性第二质子辐射带结构.      

利用质子强度变化对太阳质子事件进行短期预报方法初步研究

采用GOSE-10卫星4～9 MeV(P2),9～15 MeV(P3),15～40 MeV(P4),40～80 MeV(P5)能段上的质子通量数据,结合质子能谱,对太阳质子事件发生前各能谱参数的变化特征进行分析,详细介绍利用能谱参数的变化特征及能量E>10 MeV的质子通量数据对太阳质子事件进行预报的新方法,并运用这种方法对2002-2006年期间太阳质子事件进行了预报.预报结果显示,预报提前量最多达到100 h以上,对质子事件的报准率达97.5%,预报方法具备一定的有效性和实用性.      

(70~100)MeV准单能中子参考辐射场设计

国内基于串列加速器建立了144keV~15MeV的单能中子参考辐射场，解决了20MeV以下能区中子探测器注量、剂量等参数的量值溯源，但尚未建立20MeV以上能区准单能中子参考辐射场。随着航天空间中子探测任务的需要以及地面高能质子加速器的发展，20MeV以上能区准单能中子参考辐射场需求越来越突出。发展20MeV以上能区高能中子参考辐射场对于满足国防需求、推动中子计量学科发展具有重要意义。     

(70~100)MeV准单能中子参考辐射场设计

国内基于串列加速器建立了144keV~15MeV的单能中子参考辐射场，解决了20MeV以下能区中子探测器注量、剂量等参数的量值溯源，但尚未建立20MeV以上能区准单能中子参考辐射场。随着航天空间中子探测任务的需要以及地面高能质子加速器的发展，20MeV以上能区准单能中子参考辐射场需求越来越突出。发展20MeV以上能区高能中子参考辐射场对于满足国防需求、推动中子计量学科发展具有重要意义。     

FY-3A卫星与NOAA系列卫星高能带电粒子实测结果的比较

FY-3A卫星是运行于830 km高度的太阳同步轨道气象卫星, 其搭载的空间环境监测器可观测3~300 MeV的高能质子和0.15~5.70 MeV的高能电子. FY-3A卫星在轨工作期间, 太阳活动处于由谷年向峰年过渡期, 空间环境非常平静, 探测结果显示3~300 MeV的高能质子分布主要集中在南大西洋辐射带异常区, 0.15~5.70 MeV的高能电子分布区域除南大西洋异常区外, 还分布在南北两极高纬区域. FY-3A与NOAA卫星测量结果反映出带电粒子强度及分布区域随投掷角变化的空间各向异性特征. 本文在充分考虑了带电粒子时间、空间分布差异以及比对探测器之间自身设计差异的前提下, 经过归一化处理后, 首次对两颗卫星同期探测结果进行相关性分析, 验证了两颗卫星相同时空条件下高能带电粒子通量分布的一致性; 说明FY-3A空间环境监测器不仅具备空间带电粒子辐射监测能力, 且探测结果有效可靠, 可用作辐射带环境数据源的组成部分, 为发展新的模型, 深入研究辐射带高能粒子的分布、起源和传输等提供新的观测依据.      

利用ACE卫星数据对太阳质子事件预警方法的研究

通过分析GOES和ACE卫星大于10MeV能量段的5min平均质子通量数据,发现两者有很好的相关性,最佳相关系数显示,GOES通量数据较ACE卫星数据有数十分钟至几小时的时间延迟.这为利用ACE数据进行质子事件预警提供了依据.本文提出一种利用ACE卫星大于10MeV能量段的通量数据对质子事件进行预警的方法.该方法在2001年8月7日至2012年2月20日这11年多时间的试验中,表现出较高的报准率(76.3%)和较低的虚报率(14.7%),预警提前时间在数十分钟至几小时不等.该方法的结果与美国国家海洋和大气管理局(NOAA)空间天气预测中心(SWPC)使用的预警模型对比,具有更高的报准率,较低的虚报率,在预警时间提前量方面两者相当.      

质子引发的单粒子翻转率估算的研究

根据质子单粒子翻转事件的物理本质,及空间高能质子环境特性,给出了一种基于重离子单粒子翻转实验数据基础之上的估算质子引发的单粒子翻转的方法.利用这种方法计算了几种器件的翻转情况,与地面模拟实验及飞行观测结果符合得非常好.      

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

统计分析了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卫星高能电子预报模型, 预报结果与实测通量符合较好.      

MONITORING OF ENERGETIC PARTICLE ENVIRONMENT INSIDE THE CHINA-BRAZIL EARTH RESOURCE SATELLITE

On 14 October 1999, the Chinese-Brazil earth resource satellite (CBERS-1) was launched in China. On board of the satellite there was an instrument designed at Peking University to detect the energetic particle radiation inside the satellite so the radiation fluxes of energetic particles in the cabin can be monitored continuously. Inside a satellite cabin, radiation environment consists of ether penetrated energetic particles or secondary radiation from satellite materials due to the interactions with primary cosmic rays.Purpose of the detectors are twofold, to monitor the particle radiation in the cabin and also to study the space radiation environment The data can be used to study the radiation environment and their effects on the electronics inside the satelhte cabin. On the other hand, the data are useful in study of geo-space energetic particle events such as solar proton events, particle precipitation and variations of the radiation belt since there should be some correlation between the radiation situation inside and outside the satellite.The instrument consists of two semi-conductor detectors for protons and electrons respectively. Each detector has two channels of energy ranges. They are 0.5-2MeV and ≥2MeV for electrons and 5-30MeV and 30-60MeV for protons. Counting rate for all channels are up to 104/(cm2@s)and power consumption is about 2.5 W. There are also the additional functions of CMOS TID (total integrated dose) effect and direct SEU monitoring. The data of CBMC was first sent back on Oct. 17 1999 and it's almost three years from then on. The detector has been working normally and the quality of data is good.The preliminary results of data analysis of CBMC not only reveal the effects of polar particle precipitation and radiation belt on radiation environment inside a satellite, but also show some important features of the geo-space energetic particle radiation.As one of the most important parameters of space weather, the energetic charged particles have great influences on space activities and ground tech nology. CBMC is perhaps the first long-term on-board special equipment to monitor the energetic particle radiation environment inside the satellite and the data it accnmulated are very useful in both satellite designing and space research.     

Neutron monitor asymptotic directions of viewing during the event of 13 December 2006

During the recent ground level enhancement of 13 December 2006, also known as GLE70, solar cosmic ray particles of energy bigger that ∼500 MeV/nucleon propagated inside the Earth’s magnetosphere and finally accessed low-altitude satellites and ground level neutron monitors. The magnitude and the characteristics of this event registered at different neutron monitor stations of the worldwide network can be interpreted adequately on the basis of an estimation of the solar particle trajectories in the near Earth interplanetary space. In this work, an extended representation of the Earth’s magnetic field was realized applying the Tsyganenko 1989 model. Using a numerical back-tracing technique the solar proton trajectories inside the magnetospheric field of the Earth were calculated for a variety of particles, initializing their travel at different locations, covering a wide range of energies. In this way, the asymptotic directions of viewing were calculated for a significant number of neutron monitor stations, providing crucial information on the Earth’s “magnetospheric optics” for primary solar cosmic rays, on the top of the atmosphere, during the big solar event of December 2006. The neutron monitor network has been treated, therefore, as a multidimensional tool that gives insights into the arrival directions of solar cosmic ray particles as well as their spatial and energy distributions during extreme solar events.     

Unexpected burst of solar activity recorded by neutron monitors during October–November 2003

During the extreme burst of solar activity in October–November 2003, a series of outstanding events distinguished by their magnitude and peculiarities were recorded by the ground based neutron monitor network. The biggest and most productive in 23rd solar cycle active region 486 generated the most significant series of solar flares among of which the flare X28/3B on November 4, 2003 was the mostly powerful over the history of X-ray solar observations. The fastest arrival of the interplanetary disturbance from the Sun after the flare event in August 1972 and the highest solar wind velocity and IMF intensity were observed during these events. In one-week period three ground level enhancements (GLEs) of solar cosmic rays were recorded by neutron monitor network (28, 29 October and 2 November 2003). Maximum proton energy in these events seems to be ranged from 5 to 10 GeV. Joint analysis of data from ground level stations (neutron monitors) and satellite measurements allows the estimation of the particle path length, the onset time of the injection on the Sun and some other proton flux characteristics.     

Real-time GLE alert in the ANMODAP Center for December 13, 2006

Within the last years, a real-time system to monitor high energy cosmic rays for space weather use has been operated at Athens cosmic ray station. Neutron monitors and satellite high resolution data in real time are used, making it possible to observe cosmic rays in dual energy range observations. In large solar energetic particle (SEP) events, ground level enhancement (GLE) can provide the earliest alert for the onset of the SEP event. This system watches for count rate increases recorded in real time by 23 neutron monitors, which triggers an alarm if a ground level enhancement (GLE) of cosmic ray intensity is detected.     

FEM code simulation of the magnetospheric proton fluxes

We report a study of the numeric solution to the diffusive transport equation for energetic protons magnetically trapped in the Earth's equatorial magnetosphere. The analysis takes into account the pertinent physical processes in this region, including deceleration of protons by Coulomb collisional interactions with free and bound electrons, the charge exchange process, cosmic ray albedo neutron decay source, and electric and magnetic radial diffusion. These results were obtained using the Finite Element Method with magnetic moment and geomagnetic L-shell as free variables. Steady state boundary conditions were imposed at L=1 as zero distribution function and at L=7 with proton distribution function extracted from ATS 6 satellite observations. The FEM-code yields unidirectional proton flux in the energy range of 0.1–1000 MeV at the equatorial top of the geomagnetic lines, and the results are found to be in satisfactorily agreement with the empirical NASA AP-8 model proton flux within the energy range of 0.5–100 MeV. Below 500 keV, the empirical AP-8 model proton fluxes are several orders of magnitude greater than those computed with the FEM-code at L<3. This discrepancy is difficult to explain by uncertainties of boundary spectrum parameters or transport coefficients.     

The role of drift and convection–diffusion mechanisms in small energy global cosmic ray modulation: Application to the hysteresis phenomenon of proton and alpha particle satellite data

The hysteresis effect for small energies of galactic cosmic rays is due to two effects. The first is the same as for neutron monitor energies – the delay of the interplanetary processes responsible for cosmic ray modulation with respect to the initiating solar processes, according to the effective velocity of solar wind and shock waves propagation. Then, the observed cosmic ray intensity is connected to the solar activity variations during many months before the time of cosmic ray measurement. The second is caused by the time delay of small energy cosmic ray diffusion from the boundary of modulation region to the Earth’s orbit. The model describing the connection between solar activity variation and cosmic ray convection–diffusion global modulation for neutron monitor energies is here developed by taking into account also the time-lag of the small energy particle diffusion in the Heliosphere. We use theoretical results on drifts and analytically approximate the dependences of drifts from tilt angle, and take into account the dependence from the sign of primary particles, and from the sign of polar magnetic field (A > 0 or A < 0). The obtained results are applied on proton and alpha-particle satellite data. We analyze satellite 5-min data of proton fluxes with energies >1 MeV, >2 MeV, >5 MeV, >10 MeV, >30 MeV, >50 MeV, >60 MeV, >100 MeV, and in intervals 10–30 MeV, 30–60 MeV, and 60–100 MeV during January 1986–December 1999. We exclude periods with great cosmic ray increases caused by particle acceleration in solar flare events. Then, we determine monthly averaged fluxes, as well as 5-month and 11-month smoothed data. We analyze also satellite 5-min data on alpha-particle fluxes in the energy intervals 60-160 MeV, 160–260 MeV and 330–500 MeV during January 1986–May 2000. We correct observation data for drifts and then compare with what is expected according to the convection–diffusion mechanism. We assume different dimensions of the modulation region (by the time propagation X₀ of solar wind from the Sun to the boundary of modulation region), for X₀ values from 1 to 60 average months, by one-month steps. For each value of X₀ we determine the correlation coefficient between variations of expected and observed cosmic ray intensities (the estimation of cosmic ray intensities values is given in Section 3 by Eq. (9), and the determination of correlation and regression coefficients in Section 3 by Eq. (8)). The dimension of modulation region is determined by the value of X_0 max, for which the correlation coefficient reaches the maximum value. Then the effective radial diffusion coefficient and residual modulation in small energy region can be estimated.