Sources of inner radiation zone energetic helium ions: Cross-field transport versus in-situ nuclear reactions

Radial transport theory for inner radiation zone MeV He ions has been extended by combining radial diffusive transport, losses due to Coulomb friction and charge exchange reaction with local generation of ³He and ⁴He ions due to nuclear reactions taking place on the inner edge of the inner radiation zone. From interactions between high energy trapped protons and upper atmospheric constituents we have included a nuclear reaction yield helium flux source that was numerically derived from a nuclear reaction model originally developed at the Institute of Nuclear Researches of Moscow, Russia and implemented in the computer system at the University of Campinas, Brazil. Magnetospheric transport computations have been made covering the L-shell range L=1.0 to 1.6 and the resulting MeV He ion flux distributions show a strong influence of the local nuclear source mechanism on the inner zone energetic He ion content.     

Dynamics of the low altitude secondary proton radiation belt

At the interface between the upper atmosphere and the radiation belt region, there exists a secondary radiation belt consisting mainly of energetic ions that have become neutralized in the ring current and the main radiation belt and then re-ionized by collisions in the inner exosphere. The time history of the proton fluxes in the 0.64 – 35 MeV energy range was traced in the equatorial region beneath the main radiation belts during the three year period from 21 February 1984 to 26 March 1987 using data obtained with the HEP experiment on board the Japanese OHZORA satellite. During most of this period a fairly small proton flux of −1.2 cm⁻² s⁻¹ sr⁻¹ was detected on geomagnetic field lines in the range 1.05 < L < 1.15. We report a few surprisingly deep and rapid flux decreases (flux reduction by typically two orders of magnitude). These flux decreases were also long in duration (lasting up to three months). We also registered abrupt flux increases where the magnitude of the proton flux enhancements could reach three orders of magnitude with an enhancement duration of 1–3 days. Possible reasons for these unexpected phenomena are discussed.     

Maps of hydrogen isotopes at low altitudes in the inner zone from sampex observations

The PET instrument aboard the SAMPEX satellite has provided us with long-term intra-calibrated observations of geomagnetically trapped protons and deuterons in the inner zone, suitable for use in constraining the low-altitude portions of radiation belt models being developed as successors to AP-8. These observations have been summarized elsewhere (Looper et al., 1996). Here we report a detection of geomagnetically-trapped tritium at energies from 14 to 35 MeV/nuc below L = 1.2, at about 1/8 the flux of deuterium previously reported at that location and at similar energy per nucleon. We also demonstrate the utility of the SAMPEX/PET observations for measuring the east-west anisotropy in the trapped particle flux at low altitudes, which is due to displacement of particle gyrocenters from the position of observation in a region of strong flux gradients. This anisotropy is implicitly ignored in omnidirectional radiation-flux models, but it can be important to mission planners considering how to distribute shielding over the surface of oriented spacecraft in low Earth orbit.     

A first step towards proton flux forecasting

We present a preliminary version of a potential tool for real time proton flux prediction which provides proton flux profiles and cumulative fluence profiles at 0.5 and 2 MeV of solar energetic particle events, from their onset up to the arrival of the interplanetary shock at the spacecraft position (located at 1 or 0.4 AU). Based on the proton transportation model by Lario et al. [Lario, D., Sanahuja, B., Heras, A.M. Energetic particle events: efficiency of interplanetary shocks as 50 keV E < 100 MeV proton accelerators. Astrophys. J. 509, 415–434, 1998] and the magnetohydrodynamic shock propagation model of Wu et al. [Wu, S.T., Dryer, M., Han, S.M. Non-planar MHD model for solar flare-generated disturbances in the Heliospheric equatorial plane. Sol. Phys. 84, 395–418, 1983], we have generated a database containing “synthetic” profiles of the proton fluxes and cumulative fluences of 384 solar energetic particle events. We are currently validating the applicability of this code for space weather forecasting by comparing the resulting “synthetic” flux profiles with those of several real events.     

Statistical properties of SEP event flux declines

The interplanetary space is not a passive medium, which merely constitutes a scene for the propagation of previously accelerated energetic particles, but influences the distribution of particles by changing their energies as well due to interactions with magnetic field inhomogeneities. Such processes manifest themselves in the energy spectra of solar energetic particle (SEP) events. In this paper the fluxes of protons with energies of 4–60 MeV are investigated on the basis of two data sets. Both sets are homogeneous, obtained by the CPME instrument aboard the IMP 8 satellite between 1974 and 2001. The first includes all SEP events where the integral fluxes of >4 MeV protons exceeded 2 particle/cm² s sr. The other set consists of fluxes recorded in differential energy windows between 0.5 and 48 MeV. Important characteristics of SEP events include the rates of decrease of particle flux, which, as well as peak flux time, is an integral feature of the interplanetary medium within a considerable region, surrounding the observation point. The time intervals selected cover the decay phases of SEP events following flares, CMEs and interplanetary shocks of different origin. Only those parts of declines were selected, that could reasonably be described by exponential dependence, irrespective of the gradual/impulsive character of the events. It is shown that the average values of characteristic decay time, τ, and energy spectral index, γ, are all changing with the solar activity phase. Distributions of τ and γ values are obtained in SEPs with and without shocks and during different phases of events: just after peak flux and late after maximum.     

Storm-time formation of a relativistic electron belt and some relevant phenomena in other magnetospheric plasma domains

An empirical formula relating the strength of a storm given by its |Dst|_max with the L-coordinate of the peak of storm-injected relativistic electrons is one of a few well-confirmed quantitative relations found in the magneto-spheric physics. We successively extended a dataset of the formula’s basic storms with several events of high Dst-amplitude up to the highest observed |Dst|_max = 600 nT. Possible applying of the formula to the predicting of the ring-current plasma-pressure distribution and the lowest westward electrojet position for a storm are discussed. We have also analyzed the 2000–2001 years’ data on relativistic electrons from our instruments installed on EXPRESS-A (geosynchronous orbit; E_e = 0.8–6 MeV), Molniya-3 (h = 500 × 40 000 km, i = 63°; E_e = 0.8–5.5 MeV) and GLONASS (h = 20 000 km, i = 64°; E_e  l MeV) along with other correlated measurements: GOES series (E_e > 2 MeV), geomagnetic indices (Dst, AE, AL) and interplanetary parameters (solar wind, IMF). The goal is to investigate which outer conditions are most responsible for the high/low output of the storm-injected relativistic electrons. For the geosynchronous orbit, two factors are found as the necessary condition of the highest electron output: high and long-lasting substorm activity on a storm recovery phase and high velocity of solar wind. On the contrary, extremely low substorm activity surely observed during whole the storm recovery phase constitutes a sufficient condition of the non-increased after-storm electron intensity. For the first time found cases of the increased after-storm electron intensity observed at the inner L-shells with no simultaneously seen increase in the geosynchronous distances are presented.     

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

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

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

利用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模式所得的隐性第二质子辐射带结构.      

Energy content of accelerated electrons and ions in intense solar flares

The energy content of nonthermal particles in solar flares is shared between accelerated electrons and ions. It isimportant for understanding the particle acceleration mechanism in solar flares. Yohkoh observed a few intense flares which produced both strong gamma-ray lines and electron bremsstrahlung continuum. We analyze energy spectra of X-class solar flares on October 27, 1991(X6.1), November 6, 1997 (X9.4), July 14, 2000 (X5.7) and November 24, 2000 (X2.3). The accelerated electron and proton spectra are derived from a spectral analysis of their high-energy photon emission and the energy contents in >1 MeV electrons and >10 MeV protons are estimated to be 6×l0²⁸ – 4×10³⁰ and 2×10²⁸ – 5×10²⁹ erg, respectively. We study the flare to flare variation in the energy content of >1 MeV electrons and >10 MeV protons for the four Yohkoh gamma-ray flares. Ratios of >1 MeV electron energy content to >10 MeV proton energy content are roughly within an order of magnitude.     

Spatial and temporal variations of energetic trapped helium ions measured with Japanese satellite OHZORA

The results of measurements of absolute flux values and long term temporal evolution of the spatial distribution of trapped He ions in the energy range 1.2 – 9.2 MeV/nucleon below L = 4 are reported. The observations were made with ion counter on board the Japanese OHZORA satellite during the period of January 1984 through March 1987.     

Measurements of high-energy neutron and proton fluxes on-board “Mir-Spectr” orbital complex

The measurements of high-energy neutron (with energies 30–300 MeV) and proton (with energies 1–200 MeV) fluxes are being conducted on-board “Mir-Spectr” orbital complex. Neutrons are detected by the undirected (FOV 4π sr) scintillator spectrometer, consisting of 4 identical CsI(Tl) detector units (the effective area for neutrons 30 cm²). The gamma-quanta, which can be also detected by this instrument, are separated from neutrons by the analysis of the scintillator output pulse shape. To exclude registration of charged particles an anticoincidence plastic scintillator shield is realized in each detector unit. The proton fluxes are measured by the telescope based on 3 semiconductor detectors with small geometry factor (1 cm²×sr). As the first result of the experiment the upper limit of the integral flux of local and albedo neutrons in the equatorial region (L<1.1) was estimated. The results of this measurements can be useful for the radiation security. Also, the neutrons of solar flares can be detected in this experiment.     

南大西洋异常区高能质子通量南北向差异

利用2012年8月NOAA N15,N16,N17,N18卫星中能质子和电子探测器（MEPED）>6.9MeV质子能档的数据,研究了800km高度处南大西洋异常区（SAA）质子分布的南北方向差异.数据显示,质子计数率沿经线随纬度的变化近似满足高斯分布,且向南运动的质子多于向北运动的质子.MEPED具有0°和90°两个探测方向,可对质子投掷角分布进行估算,从而消除在探测南向和北向质子时探测器方向角偏差所造成的影响.结果表明,在经度40°W、纬度13°S至23°S位置上,仍然存在南向质子比北向质子多20%~30%的现象.其原因可能与地磁异常使南向质子镜像点更低,更多质子散射进入大气层沉降有关.分析还发现,南大西洋异常区质子通量存在地方时变化.日侧比夜侧强约20%,这可能是由于地磁场日变化引起的.      

A model for stochastic acceleration of electrons during geomagnetic storms

The theory of resonant diffusion is extended to fully relativistic plasmas, and we examine resonant interactions between electrons and electromagnetic R mode (whistler) and L-mode (EMIC) waves. Resonant diffusion curves are constructed for plasma parameters representative of the Earth's storm time magnetosphere, both inside and outside the plasmapause. EMIC waves can resonate with electrons > 1 MeV, but the energies remain nearly constant along the diffusion curves. Storm-time EMIC waves can induce rapid pitch—angle scattering, but the waves are ineffective for stochastic acceleration of elections. Substantial energy change can occur along the diffusion curves for interactions between resonant electrons and whistler—mode waves, especially in regions of low plasma density. Specifically, whistlers can accelerate electrons from energies near 100 keV to above 1 MeV outside the plasmapause. A model is proposed comprising energy diffusion by whistler-mode chorus and pitch-angle scattering by EMIC waves to account for the gradual acceleration of electrons over the region 4 ≤ L ≤ 6 during the recovery phase of a geomagnetic storm.     

Solar wind dependence of the Pc1 wave activity

There are a host of factors influencing the excitation of Pc1 geomagnetic pulsations, which are ULF waves in the frequency range between 0.2 and 5 Hz. We have studied carefully the dependence of the pearl-type Pc1 activity at Sodankylä, Finland (L = 5.1) on the plasma density N in front of the magnetosphere, the bulk velocity V of the solar wind, and the intensity B of the IMF. The result is as follows: high values of N and reduced values of V are favorable to appearance of Pc1, whereas the dependence of Pc1 activity on B is practically absent. We also show that the probability of Pc1 occurrence decreases with the interplanetary electric field, and increases with solar wind impact pressure and with the plasma to magnetic pressure ratio “beta”.     

Creation of model of quasi-trapped proton fluxes below Earth’s radiation belt

The object of investigation is the phenomenon of proton (from tens keV to several MeV) flux enhancement in near-equatorial region (L < 1.15) at altitude up to ∼1300 km (the storm-time equatorial belt). These fluxes are quite small but the problem of their origin is more interesting than the possible damage they can produce. The well known sources of these protons are radiation belt and ring current. The mechanism of transport is the charge-exchange on neutral hydrogen of exosphere and the charge-exchange on oxygen of upper atmosphere. Therefore this belt is something like the ring current projection to low altitudes. Using the large set of satellites data we obtain the average energy spectrum, the approximation of spectrum using kappa-function, the flux dependence on L, B geomagnetic parameters. On the basis of more than 30 years of experimental observations we made the empiric model that extends model of proton fluxes below 100 keV in the region of small L-values (L < 1.15). The model was realized as the package of programs integrated into COSRAD system available via Internet. The model can be used for revision of estimation of dose that low-orbital space devices obtain.     

Some peculiarities of longitudinal distribution of proton fluxes at high latitudes

Dynamical features of proton fluxes at high and middle latitudes were studied based on data measured by Sun-synchronous low-altitude (800 km height) polar-orbiting three NOAA series satellites: POES 15, 16, and 17 during the geomagnetic storm on January, 21–22, 2005. Data from three satellites that passed the Northern hemisphere along different MLTs allow reconstructing the longitudinal distribution of the proton fluxes. Measurements of protons with energies of 30–80 keV and 80–240 keV (the ring current energy range) by 0- and 90-detectors were used to evaluate and compare the longitudinal asymmetry of proton flux distribution measured in the regions equatorward and poleward of the isotropic boundary. It was found that during all the phases of the geomagnetic storm distribution of the maximum flux of precipitating protons (0-detector data) is sufficiently asymmetric. The maximal flux position along MLT is moving from pre-midnight sector in quiet time to post-midnight one before and during SSC and moving back during recovery phase. The longitudinal distribution of precipitation maxima demonstrates the local increase in afternoon sector (approximately at 13:30 MLT) and decrease in the dusk one during SSC. These features are evident consequence of the magnetosphere compression. To identify the origin of the particles, the locations of maximum fluxes have been projected to the magnetosphere. It was determined that during geomagnetic storm main and recovery phases maximum fluxes were measured at latitudes poleward of the isotropic boundary. To evaluate the trapped particle flux asymmetry, the particle fluences (90-detector data) were calculated along the satellite orbit from L = 2 to the isotropic boundary. The total fluences of trapped particles calculated along the satellite orbit show regular asymmetry between dusk and dawn during main and recovery phases. The maximal intensity of proton fluxes of both investigated populations located poleward and equatorward of the isotropic boundary is achieved during SSC. The total flux measured during crossing the anisotropic region can be considered as a proxy for ring current injection rate.     

Solar cycle variation of the low-altitude trapped proton flux

Under NASA's Space Environment Effects (SEE) program, we are developing new models for the low-altitude (250–1000 km, L < 1.5) trapped radiation environment based on data from the TIROS/NOAA polar orbiting spacecraft. The unique features of this data base and model include the long time series (more than one complete solar cycle) obtained from the TIROS/NOAA data and the use of a coordinate system more applicable to the low-altitude environment. The data show a strong variation (as much as a factor of 10) over the solar cycle and a hysteresis effect between the rising and falling portions of the solar cycle. Both the solar cycle variation and the hysteresis are functions of L. In addition to the hysteresis effect, the flux during a given cycle appears to be a function of the previous cycle. Superimposed on the gradual variation over the solar cycle, transient effects, correlated with solar particle events (SPEs), can be clearly seen. Comparison with the AP8 models shows that the measured flux is a factor of 2–3 higher than the model. These data have important implications for the development and use of trapped radiation models, and will also contribute to our knowledge of the source and loss mechanisms at work in the inner zone.     

Energetic particle fluxes data base of “CORONAS-I” satellite observations

The data base DB SCR uses data obtained by the SCR instrument package on the CORONAS-I satellite. DB SCR contains information about fluxes of relativistic electrons (0.5–124 MeV), protons (1–300 MeV), nuclei (1–19 MeV/nucl) and γ-radiation (0.1–7 MeV) in the low altitude region (500 km). The time resolution of the data is 2.5 s. Magnetic field parameters (B, magnetic latitude and longitude), L-shell and local times (LT and MLT) are included in DB SCR. Since all parameters are equivalents it is possible to perform the multidimensional analysis for any set of DB SCR parameters. The additional DB SCR software packages may be used to develop different semi-empirical models.     

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

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

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.