High energy charged particle experiment – A payload proposal to KuaFu-B mission

High Energy Charged Particle Experiment (HECPE) is to measure the fluxes of MeV electrons and tens of MeV protons. The two satellites of KuaFu-B are in the same polar orbit with apogee 7.0R_E, perigee 1.8R_E. They can sweep large L values and pass through the inner and outer radiation belts. The high energy electrons and protons in the radiation belts are principal sources for failures of satellites and spacecrafts in the Earth orbits. The enhancements of the high energy electrons and protons, so-called energetic particle events, are important phenomena of the Space Weather. The energy ranges monitored by HECPE are 0.3–0.5 MeV, 0.5–1.0 MeV, 1.0–2.0 MeV, and E > 2.0 MeV for electrons, 5–10 MeV, 10–20 MeV, 20–40 MeV, and 40–80 MeV for protons.     

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.     

Long-term observations of the trapped high-energy proton population (L < 4) by the NOAA Polar Orbiting Environmental Satellites (POES)

The Space Environment Monitor (SEM) onboard the NOAA POES satellites has been measuring the near-Earth charged particle environment since 1978, providing an extensive database that can be used for studying the long-term behavior of this population of trapped particles. POES stands for Polar Orbiting Environmental Satellite. These satellites orbit at ∼840 km altitude and at an inclination of 98°. The SEM-1 instrument was flown on the POES satellites beginning in 1978 with TIROS-N and NOAA-6 in 1979 and continuing to NOAA-14 launched in 1995 with the exception of NOAA-9 and NOAA-11 (NOAA-13 failed shortly after launch). Its replacement, SEM-2, has flown on the POES NOAA-15, -16, and -17 satellites (from 1998). Here we present the results of a study on the statistical variations of the high-energy trapped proton environment. Among the detectors in SEM, the four SEM-2 omni-directional proton detectors for energies >16 MeV, >36 MeV, >70 MeV, and >140 MeV provide the data most relevant to this study.     

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.     

Near-Earth radiation model deficiencies as seen on CRRES.

The Space Radiation (SPACERAD) experiments on the Combined Release and Radiation Effects Satellite (CRRES) gathered 14 months of radiation particle data in an 18 degrees inclination orbit between 350 km and 36000 km from July 1990 to October 1991. When compared to the NASA radiation belt models AP8 and AE8, the data show the proton model (AP8) does not take into account a second belt formed after major solar flare/shock injection events, and the electron model (AE8) is misleading, at best, in calculating dose in near-Earth orbits. The second proton belt, although softer in energy than the main proton belt, can produce upsets in proton sensitive chips and would produce significant dose in satellites orbiting in it. The MeV electrons observed on CRRES show a significant particle population above 5 MeV (not in the AE8 model) which must be included in any meaningful dose predictions for satellites operating between L-shells of 1.7 and 3.0 RE.     

An analytical model for the prediction of a micro-dosimeter response function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (δ ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCRs) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 s/ion. The ionizing radiation environment at LEO is represented by O’Neill’s GCR model (2004), covering charged particles in the 1 ? Z ? 28 range. O’Neill’s free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.     

Energy spectrum of proton flux near geomagnetic equator at low altitudes

The results of proton energy (tens keV – several MeV) spectrum measurements near geomagnetic equator (L < 1.15) at low altitudes (<1000 km) are presented. We used data of experiments onboard ACTIVE, SAMPEX, NOAA TIROS-N satellites and SPRUT-VI (MIR station) and cover a time range of about 30 years (including previous measurements). It was found that the kappa-distribution function fits the experimental spectrum with the best correlation coefficient. A comparison of energy spectra of near-equatorial protons and ring-current protons was made. Using the estimation of the life time of near-equatorial protons we explain the difference in spectral indices of radiation belt and near-equatorial proton formation. We conclude that the ring current is the main source of the near-equatorial protons.     

On the spatial extent of the proton radiation belt from solar cell output variation of the Akebono satellite

Output current of silicon solar cells of Akebono satellite orbiting in the inner magnetosphere decreased from 13 A in 1989 to about 7 A in 2009, due to accumulated damage by energetic particles. A fair correlation between the monthly decrease rate of solar cell output current and trapped proton flux was found between 1989 and 1992, indicating that trapped energetic protons are responsible for the solar cell degradation. Simple variation from month to month was dominated by a fluctuation with one-year period and the correlation was not discernible after 1993. On the other hand, during orbiting the earth, another kind of decrease of output current emerged at low altitude above the sub-solar point, i.e., in the midst of sunlit condition. The decrease was larger in later years. The fact indicates more prominent temperature effect in the later years because of progress of the degradation. By removing data affected by the orbit condition above, variation from the same month in the previous year shows a clear oscillation due to orbit precession correlated with trapped energetic proton flux up to 1996. The amplitude of the annual variation oscillation tends to be larger than that expected from a degradation model based on energetic proton distribution of the NASA’s AP8 model. The larger oscillation amplitude suggests that the proton radiation belt was more sharply localized than given by the AP8 model throughout the early half of 1990s.     

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.     

Using cross-correlation techniques to assess the validity of AP8 at low altitude

Background effects, such as single event upsets, are used as proxy data sets to evaluate AP8 at low altitude. The approach is to use a 2-D (longitude and latitude) cross-correlation between the background data set and the current-epoch AP8 predicted fluxes. The technique can be used to determine the energy of the particle that is producing a particular effect. This cross-correlation technique shows that using AP8 with a present-epoch magnetic field model accurately predicts the present location of the South Atlantic Anomaly proton flux enhancement and that a dual-peaked intensity structure above 100 MeV in the model is an artifact of AP8 and its interpolation routines.     

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.     

Change in the rigidity dependence of the galactic cosmic ray modulation in 2008–2009

Observations of charged particle fluxes in the stratosphere of the polar regions represent the cosmic rays variations with energy above 100 MeV. At the end of 2009 these fluxes reached the highest level for the time of observations from mid 1957 and were by 17% higher than the previous extremum value of May 1965. In the mean time the ground-based neutron monitors showed the remarkably less count rate enhancement. These results argue for the significant change in the energy spectrum of incoming particles in 2008–2009 in the energy range of ∼100–1500 MeV/n.     

The hydroxyl radical as an indicator of SEP fluxes in the high-latitude terrestrial atmosphere

The low background values at nighttime of the mesospheric hydroxyl (OH) radical make it easier to single out the atmospheric response to the external solar forcing in Polar Regions. Because of the short lifetime of HO_x, it is possible to follow the trails of Solar Energetic Particle (SEP) events in the terrestrial atmosphere, as shown by Storini and Damiani (2008). The sensitivity of this indicator makes discernible not only extreme particle events with a flux peak of several thousand pfu [1 pfu = 1 particle/(cm² s sr)] at energies >10 MeV, but also those with lower flux up to about 300 pfu. Using data from the Microwave Limb Sounder (MLS) on board the EOS AURA satellite, we examined the correlation of OH abundance vs. solar proton flux for almost all the identified SEP events spanning from November 2004 to December 2006 (later on no more SEP events occurred during Solar Cycle no. 23). The channels at energies greater than 5 MeV and 10 MeV showed the best correlation values (r ∼ 0.90–0.95) at altitudes around 65–75 km whereas, as expected, the most energetic channels were most highly correlated at lower altitudes. Therefore, it is reasonably possible to estimate the solar proton flux from values of mesospheric OH (and viceversa) and it could be useful in studying periods with gaps in the records of solar particles.     

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.     

Solar proton events recorded in the stratosphere during cosmic ray balloon observations in 1957–2008

Long-term balloon observations have been performed by the Lebedev Physical Institute since 1957 up to the present time. The observations are taken several times a week at the polar and mid latitudes and allow us to study dynamics of galactic and solar cosmic ray as well as secondary particle fluxes in the atmosphere and in the near-Earth space. Solar energetic particles (120) – mostly protons – (SEP) events with >100 MeV proton intensity above 1 cm⁻² s⁻¹ s⁻¹ were recorded during 1958–2006. Before the advent of the SEP monitoring on spacecraft these results constituted the only homogeneous series of >100 MeV SEP events. The SEP intensities and energy spectra inferred from the Lebedev Physical Institute observations are consistent with the results taken in the adjacent energy intervals by the spacecraft and neutron monitors. Joint consideration of the SEP events series recorded by balloons and by neutron monitors during solar cycles 20–23 makes it possible to restore the probable number of events in solar cycle 19, which was not properly covered by observations. Some correlation was found between duration of SEP event production in a solar cycle and sunspot cycle characteristics.     

Relativistic electron fluxes and dose rate variations during April–May 2010 geomagnetic disturbances in the R3DR data on ISS

Space radiation has been monitored successfully using the Radiation Risks Radiometer-Dosimeter (R3D) installed at the ESA EXPOSE-R (R3DR) facility outside of the Russian Zvezda module of the International Space Station (ISS) between March 2009 and January 2011. R3DR is a Liulin type spectrometer–dosimeter with a single Si PIN detector 2 cm² of area and 0.3 mm thick. The R3DR instrument accumulated about 2 million measurements of the absorbed dose rate and flux of 10 s resolution. The total external and internal shielding before the detector of R3DR device is 0.41 g cm⁻². The calculated stopping energy of normally incident particles to the detector is 0.78 MeV for electrons and 15.8 MeV for protons. After the Coronal Mass Ejection (CME) at 09:54 UTC on 3 April 2010, a shock was observed at the ACE spacecraft at 0756 UTC on 5 April, which led to a sudden impulse on Earth at 08:26 UTC. Nevertheless, while the magnetic substorms on 5 and 6 of April were moderate; the second largest in history of GOES fluence of electrons with energy >2 MeV was measured. The R3DR data show a relatively small amount of relativistic electrons on 5 April. The maximum dose rate of 2323 μGy day⁻¹ was reached on 7 April; by 9 April, a dose of 6600 μGy was accumulated. By the end of the period on 7 May 2010 a total dose of 11,587 μGy was absorbed. Our data were compared with AE-8 MIN, CRESS and ESA-SEE1 models using SPENVIS and with similar observations on American, Japanese and Russian satellites.     

Highly Elliptical Orbits for Arctic observations: Assessment of ionizing radiation

The ionizing radiation environment was analyzed for a variety of potential Highly Elliptical Orbits (HEOs) with orbital periods ranging from 6 h to 24 h suitable to continuously monitor the Arctic region. Several models available from the ESA Space Environment Information System (SPENVIS) online tool were employed, including the new-generation AE9/AP9 model for trapped radiation. Results showed that the Total Ionizing Dose (TID) has a well-pronounced local minimum for the 14-h orbit, which is nearly identical to the overall minimum observed for the longest orbital period (24 h). The thickness of slab aluminum shielding required to keep the annual TID below 10, 5 and 3.33 krad (i.e. 150, 75 and 50 krad for 15 years of mission duration) for a 14-h orbit is 2.1, 2.7 and 3.1 mm respectively. The 16-h orbit requires an additional 0.5 mm of aluminum to achieve the same results, while the 24-h orbit requires less shielding in the order of 0.2–0.3 mm. Comparison between the AE8/AP8 and AE9/AP9 models was conducted for all selected orbits. Results demonstrated that differences ranged from −70% to +170% depending on orbit geometry.     

Photoelectron flux variations observed from the FAST satellite

This paper examines high resolution (ΔE/E = 0.15) photoelectron energy spectra from 10 eV to 1 keV, created by solar irradiances between 1.2 and 120 nm. The observations were made from the FAST satellite at ∼3000 km, equatorward of the auroral oval for the July–August, 2002 solar rotation. These data are compared with the solar irradiance observed by the Solar EUV Experiment (SEE) on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite and fluxes calculated using the Field Line Interhemispheric Plasma (FLIP) code. The 41 eV photoelectron flux, which corresponds to solar EUV fluxes near 20 nm, shows a clear solar rotation variation in very good agreement with the EUV flux measurements. This offers the possibility that the 41 eV photoelectron flux could be used as a check on measured solar EUV fluxes near 20 nm. Because of unexpected noise, the solar rotation signal is not evident in the integral photoelectron flux between 156 and 1000 eV corresponding to EUV wavelengths between 0.1 and 7 nm measured by the SEE instrument. Examination of daily averaged photoelectron fluxes at energies between 25 and 500 eV show significant changes in the photoelectron spectra in response X and M class flares. The intensity of photoelectrons produced in this energy region is primarily due to two very narrow EUV wavelength regions at 2.3 and 3 nm driving Auger photoionization in O at 500 eV and N₂ at ∼360 eV. Comparison of calculated and daily averaged electron fluxes shows that the HEUVAC model solar spectrum used in the FLIP code does not reproduce the observed variations in photoelectron intensity. In principle, the 21 discrete photoelectron energy channels could be used to improve the reliability of the solar EUV fluxes at 2.3 and 3 nm inferred from broad band observations. In practice, orbital biases in the way the data were accumulated and/or noise signals arising from natural and anthropogenic longitudinally restricted sources of ionization complicate the application of this technique.     

Simultaneous observations of Pc 1 micropulsation activity and stratospheric electrodynamic perturbations on 27 January 2003

The 2nd Polar Patrol Balloon campaign (2nd-PPB) was carried out at Syowa Station in Antarctica during 2002–2003. Identical stratospheric balloon payloads were launched as close together in time as allowed by weather conditions to constitute a cluster of balloons during their flights. A very pronounced negative ion conductivity enhancement was observed at 32 km in the stratosphere below the auroral zone on 27 January 2003 from 1500 to 2200 UT. During this event, the conductivity doubled for an interval of about 7 h. This perturbation was associated with an extensive Pc 1 or Pi 1 wave event that was observed by several Antarctic ground stations, balloon PPB 10, and the Polar spacecraft. No appreciable X-ray precipitation was observed in association with this event, which would point to >60 Mev proton precipitation as a possible magnetosphere–stratosphere coupling mechanism responsible for the conductivity enhancement. Such precipitation is consistent with the wave data. During the latter half of the event, E_z was briefly positive. There was a tropospheric Southern Ocean storm system underneath the balloon during this interval. If the event was associated with this storm system and not energetic proton precipitation, the observations imply an electrified Southern Ocean storm and major perturbations in stratospheric conductivity driven by a tropospheric disturbance. This event represents a poorly understood source for global circuit current. Precipitating energetic proton data from Akebono and NOAA POES spacecraft show significant >16 MeV precipitation was occurring at the location of PPB 8 but not PPB 10, suggesting that proton precipitation was, in fact, the responsible coupling mechanism.