The effects of proton exposure on neurochemistry and behavior.

Future space missions will involve long-term travel beyond the magnetic field of the Earth, where astronauts will be exposed to radiation hazards such as those that arise from galactic cosmic rays. Galactic cosmic rays are composed of protons, alpha particles, and particles of high energy and charge (HZE particles). Research by our group has shown that exposure to HZE particles, primarily 600 MeV/n and 1 GeV/n 56Fe, can produce significant alterations in brain neurochemistry and behavior. However, given that protons can make up a significant portion of the radiation spectrum, it is important to study their effects on neural functioning and on related performance. Therefore, these studies examined the effects of exposure to proton irradiation on neurochemical and behavioral endpoints, including dopaminergic functioning, amphetamine-induced conditioned taste aversion learning, and spatial learning and memory as measured by the Morris water maze. Male Sprague-Dawley rats received a dose of 0, 1.5, 3.0 or 4.0 Gy of 250 MeV protons at Loma Linda University and were tested in the different behavioral tests at various times following exposure. Results showed that there was no effect of proton irradiation at any dose on any of the endpoints measured. Therefore, there is a contrast between the insignificant effects of high dose proton exposure and the dramatic effectiveness of low dose (<0.1 Gy) exposures to 56Fe particles on both neurochemical and behavioral endpoints.     

Evaluation of solar proton spectra using balloon cosmic ray observations and Monte Carlo simulation results

We have developed a method to evaluate the spectrum of solar energetic protons at the top of the Earth’s atmosphere from the measurements of our balloon cosmic ray experiment. By using the Monte Carlo PLANETOCOSMICS code based on Geant4 we compute the interaction of solar protons [10 MeV–10 GeV] with the Earth’s atmosphere. We obtain the angular and energy distributions of secondary particles (p, e⁻, e⁺, photons, muons) at different atmospheric levels as a function of primary proton spectra. By comparing the calculated depth dependence of the particle flux with the data obtained by our balloon experiment we can deduce the parameters of the solar proton spectrum that best fit the observations. In this paper we discuss our solar proton spectrum estimation method, and present results of its application to selected solar proton events from 2001 to 2005.     

The energy spectra of protons and helium measured with the ATIC experiment

The Advanced Thin Ionization Calorimeter (ATIC) balloon experiment is designed to investigate the composition and energy spectra of cosmic rays at the highest energies currently accessible by direct measurements, i.e., the region up to 100 TeV. The instrument consists of a silicon matrix for charge measurement, a graphite target (0.75 nuclear interaction length) to induce hadronic interactions, three layers of scintillator strip hodoscopes for triggering and trajectory reconstruction, and a Bismuth Germanate (BGO) crystal calorimeter (18 radiation lengths) to measure particle energies. ATIC has had two successful Long Duration Balloon (LDB) flights from McMurdo, Antarctica: one from 12/28/00 to 01/13/01 and the other from 12/29/02 to 01/18/03. We present the energy spectra of protons and helium extracted from the first flight, over the energy range from 100 GeV to 100 TeV, and compare them with the results from other experiments at both the lower and higher energies. ATIC-1 results do not indicate significant differences in spectral shape between protons and helium over the investigated energy range.     

On the production of highest energy solar protons at 20 January 2005

On January 20, 2005, 7:02–7:05 UT the Aragats Multidirectional Muon Monitor (AMMM) located at 3200 m a.s.l. registered enhancement of the high energy secondary muon flux (threshold ∼5 GeV). The enhancement, lasting for 3 min, has statistical significance of ∼4σ and is related to the X7.1 flare seen by the GOES, and very fast (>2500 km/s) CME seen by SOHO, and the Ground Level Enhancements (GLE) #69 detected by the world-wide network of neutron monitors and muon detectors. The energetic and temporal characteristics of the muon signal from the AMMM are compared with the characteristics of other monitors located at the Aragats Space-Environmental Center (ASEC) and with other neutron and muon detectors. Since secondary muons with energies >5 GeV are corresponding to solar proton primaries with energies 20–30 GeV we conclude that in the episode of the particle acceleration at 7:02–7:05 UT 20 January 2005 solar protons were accelerated up to energies in excess of 20 GeV.     

Countermeasures for space radiation induced adverse biologic effects

Radiation exposure in space is expected to increase the risk of cancer and other adverse biological effects in astronauts. The types of space radiation of particular concern for astronaut health are protons and heavy ions known as high atomic number and high energy (HZE) particles. Recent studies have indicated that carcinogenesis induced by protons and HZE particles may be modifiable. We have been evaluating the effects of proton and HZE particle radiation in cultured human cells and animals for nearly a decade. Our results indicate that exposure to proton and HZE particle radiation increases oxidative stress, cytotoxicity, cataract development and malignant transformation in in vivo and/or in vitro experimental systems. We have also shown that these adverse biological effects can be prevented, at least partially, by treatment with antioxidants and some dietary supplements that are readily available and have favorable safety profiles. Some of the antioxidants and dietary supplements are effective in preventing radiation induced malignant transformation in vitro even when applied several days after the radiation exposure. Our recent progress is reviewed and discussed in the context of the relevant literature.     

Comparison of genotoxic damage in monolayer cell cultures and three-dimensional tissue-like cell assemblies

Assessing the biological risks associated with exposure to the high-energy charged particles encountered in space is essential for the success of long-term space exploration. Although prokaryotic and eukaryotic cell models developed in our laboratory and others have advanced our understanding of many aspects of genotoxicity, in vitro models are needed to assess the risk to humans from space radiation insults. Such models must be representative of the cellular interactions present in tissues and capable of quantifying genotoxic damage. Toward this overall goal, the objectives of this study were to examine the effect of the localized microenvironment of cells, cultured as either 2-dimensional (2D) monolayers or 3-dimensional (3D) aggregates, on the rate and type of genotoxic damage resulting from exposure to Fe-charged particles, a significant portion of space radiation. We used rodent transgenic cell lines containing 50–70 copies of a LacI transgene to provide the enhanced sensitivity required to quantify mutational frequency and type in the 1100-bp LacI target as well as assessment of DNA damage to the entire 45-kbp construct. Cultured cells were exposed to high-energy Fe charged particles at Brookhaven National Laboratory’s Alternating Gradient Synchrotron facility for a total dose ranging from 0.1 to 2 Gy and allowed to recover for 0–7 days, after which mutational type and frequency were evaluated. The mutational frequency was found to be higher in 3D samples than in 2D samples at all radiation doses. Mutational frequency also was higher at 7 days after irradiation than immediately after exposure. DNA sequencing of the mutant targets revealed that deletional mutations contributed an increasingly high percentage (up to 27%) of all mutations in cells as the dose was increased from 0.5 to 2 Gy. Several mutants also showed large and complex deletions in multiple locations within the LacI target. However, no differences in mutational type were found between the 2D and the 3D samples. These 3D tissue-like model systems can reduce the uncertainty involved in extrapolating risk between in vitro cellular and in vivo models.     

On the quality of mutations in mammalian cells induced by high LET radiations

The deleterious effects of accelerated heavy ions as component of the space radiation environment on living cells are of increasing importance for long duration human space flight activities. The most important aspect of such densely ionizing particle radiation is attributed to the type and quality of biological damage induced by them. This issue is addressed by investigating cell inactivation and mutation induction at the Hprt locus (coding for hypoxanthine-guanine-phosphoribosyl-transferase) of cultured V79 Chinese hamster cells exposed to densely ionizing radiation (accelerated heavy ions with different LETs from oxygen to gold, specific energies ranging from 1.9 to 69.7 MeV/u, corresponding LET values range from 62 to 13,223 keV/μm) and to sparsely ionizing radiation (200 kV X-rays). 30 spontaneous, 40 X-ray induced and 196 heavy ion induced 6-thioguanine resistant Hprt mutant colonies were characterized by Southern technique using the restriction enzymes EcoRI, PstI and BglII and a full length Hprt cDNA probe isolated from the plasmid pHPT12. Restriction patterns of the spontaneous Hprt mutants were indistinguishable from the wild type pattern, as these mutants probably contain only small deletions or even point mutations in the Hprt locus. In contrast, the overall spectrum of heavy ion induced mutations revealed a majority of partial or total deletions of the Hprt gene. With constant particle fluence (3 × 10⁶ particles/cm²) the quality of heavy ion induced mutations in the Hprt locus depends on physical parameters of the beam (atomic number, specific energy, LET). This finding suggests a relationship between the type of DNA damage and track structure. The fraction of mutants with severe deletions in the Hprt locus after exposure to oxygen ions increases from 65% at 60 keV/μm up to a maximum (100%) at 300 keV/μm and declines with higher LET values to 75% at 750 keV/μm. With heavier ions (Ca- and Au-ions) and even higher LET-values this mutant fraction decreases to 58% at 13,200 keV/μm. Heavy ion induced DNA break points in the Hprt locus are not randomly distributed.     

Effect of random nature of cosmic ray sources – Supernova remnants – on cosmic ray intensity fluctuations,anisotropy, and electron energy spectrum

The random nature of sources (the supernova remnants) leads to the fluctuations of cosmic ray intensity in space and time. We calculate the expected fluctuations in a flat-halo diffusion model for particles with energies from 0.1 to 10³ TeV. The data on energy spectra and anisotropy of very high energy protons, nuclei and electrons, and the astronomical data on supernova remnants, the potential sources of cosmic rays, are used to constrain the value of the cosmic-ray diffusion coefficient and its dependence on energy.     

The Carrington event: Possible doses to crews in space from a comparable event

Recent ice core analyses suggest that the Carrington event of 1859 may have been the largest solar energetic particle event in the past several hundred years. Previous analyses of potential doses to humans and electronics from such an event suggested that a Carrington-like event, with a hard spectrum similar to that of the event of September 1989 could be catastrophic. Subsequent analyses of the ¹⁰Be concentration in the ice core data suggest that the spectral hardness of the Carrington event was softer and similar to the August 1972 event. In this work we review the earlier estimates of doses from a Carrington event, and present updated dose estimates for deep space crews and electronics using the Carrington event proton fluence ⩾30 MeV in combination with an event spectrum similar to that of the August 1972 event. Potential ramifications of these doses for humans and electronics on deep space missions are discussed.     

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.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.     

Modeling of effect of polarization on UV sky radiance during twilight

The spatial distribution of the vector of the Stokes parameters characterizing the radiance intensity and the radiance polarization describes the radiation field in the atmosphere. A simplified treatment of light as the scalar has only restricted application. A few studies compared previously results of the vector and scalar radiative transfer models and showed that scalar models are in error by up to 10% for many cases. Though several observational conditions were exploited, an effect of polarization on modeling of UV radiance has not been investigated yet for twilight. The paper presents a preliminary study of modeled UV radiance during twilight taking into account polarization. The intensity and the degree of linear polarization of the scattered UV radiance for two cases of the ground-based observations are discussed. In the first case, radiation incoming from the zenith for the solar zenith angles (SZA) from 90° to 98° is under investigation. Radiation in the solar principal plane for the beginning of twilight (SZA = 90.1°) was calculated in the second case. The study showed that the UV radiation field in the twilight atmosphere can be handled correctly only using the vector theory. The errors of scalar radiative transfer strongly depend on wavelength, line of an observation and solar position. The revealed distortion of the zenith radiance caused by using of the scalar approximation reaches maximum of 15% at 340 nm for the solar zenith angle (SZA) equal to 98°. The shorter wavelengths have the smaller errors, about 5% at 305 nm for SZA = 98°, due to the larger part of the single scattered radiance. The error of the scalar modeling may be as large as −17% for radiance incoming from the horizon for SZA = 90.1°. Scalar radiative transfer models underestimate the integral intensity in the principal plane up to 3–4% ± 0.5% at SZA = 90.1° for wavelengths from 320 to 340 nm. This should be taken into account in problems of radiative budget estimation and remote sensing of the atmosphere exploiting the twilight period.     

Simulation study of motion of charged particles trapped in Earth’s magnetosphere

This article aims to understand the motion of the charged particles trapped in the Earth’s inner magnetosphere. The emphasis is on identifying the numerical scheme, which is appropriate to characterize the trajectories of the charged particles of different energies that enter the Earth’s magnetosphere and get trap along the magnetic field lines. These particles perform three different periodic motions, namely: gyration, bounce, and azimuthal drift. However, often, the gyration of the particle is ignored, and only the guiding center of the particle is traced to reduce the computational time. It is because the simulation of all three motions (gyro, bounce, and drift) together needed a robust numerical scheme, which has less numerical dissipation. We have developed a three-dimensional test particle simulation model in which the relativistic equation of motion is solved numerically using the fourth and sixth-order Runge-Kutta methods. The stability of the simulation model is verified by checking the conservation of total kinetic energy and adiabatic invariants linked with each type of motion. We found that the sixth-order Runge-Kutta method is suitable to trace the complete trajectories of both proton and electron of a wide energy range, 5 keV to 250 MeV for L = 2 – 6. We have estimated the bounce and drift periods from the simulations, and they are found to be in good agreement with the theory. The study implies that a simulation model with sixth-order Runge-Kutta method can be applied to the time-vary, non-analytical form of magnetic configuration in future studies to understand the dynamics of charged particles trapped in Earth’s magnetosphere.     

Energetic particle transport in anisotropic magnetic turbulence

We study energetic particle transport in a magnetic field configuration which models the solar wind magnetic turbulence plus the background field. A power-law Fourier amplitude is used for the fully 3D turbulence model, and in order to model anisotropic turbulence, the constant amplitude surfaces in k space are ellipsoids. The turbulence correlation lengths parallel (perpendicular) to the background magnetic field l_∥ (l_⊥) are varied in a wide range, and proton energies from 1 MeV to 10 GeV are assumed. Considering propagation on a distance corresponding to 1 AU, it is found that transport parallel and perpendicular to the background field heavily depends on the turbulence anisotropy, that is on the ratio l_∥/l_⊥. The spatial distribution of energetic particle follows the shape of magnetic flux tube up to about 10 MeV, while for larger energies the structure of the magnetic flux tube is progressively washed out. The scatterplots of particle distribution show intermittent, non Gaussian structures for l_∥ ≪ l_⊥ (quasi slab turbulence), while a more diffusive, Gaussian structure is obtained for l_∥ ≫ l_⊥ (quasi 2D turbulence). The long time behavior of transport shows that anomalous (subdiffusive perpendicular and superdiffusive parallel) transport regimes are obtained for l_∥ ≪ l_⊥, while Gaussian diffusive transport is obtained for both l_∥ ≫ l_⊥ and the isotropic turbulence case.     

Proton versus electron heating in solar flares

Proton and electron heating of a flaring atmosphere is compared in a kinetic approach for the particles ejected from a non-neutral reconnecting current sheet (RCS) located above the top of reconnected flaring loops in a two-ribbon flare. Two kinds of high-energy particles are considered: particles accelerated by a super-Dreicer electric field and those ejected from the reconnection region as neutral outflows, or separatrix jets. The beam electrons are assumed to deposit their energy in Coulomb collisions and Ohmic heating of the ambient plasma particles by the electric field induced by the precipitating beams. The protons are assumed to deposit their energy in generation of kinetic Alfvén waves (KAWs), which, in turn, dissipate due to Cherenkov resonant scattering on the ambient plasma electrons. The beam electrons are found to provide a fast (within a few tenth of a second) heating of the atmosphere that is well spread in depth from the corona to the lower chromosphere. The protons are shown to precipitate to the lower atmosphere much slower (up to few seconds for beam and up to 10–20 s for slow jets). Slow jet protons provide heating of the two compact regions: the first located at the top of a flaring loop just below the RCS, and the second one appearing at the transition region (TR) and upper chromosphere; fast beam protons deposit their energy in the TR and chromosphere only.     

Single event upset investigations on the “Flying Laptop” satellite mission

The radiation effects in electronic parts are called single-event effects, which are deemed to be critical for space missions. This paper presents the Single Event Upsets that were observed in an onboard memory device of the Low Earth Orbit “Flying Laptop” satellite mission during its in-orbit operation. The Single Event Upsets were carefully mapped on the satellite orbital space itself and their root causes were investigated together with their rates of occurrence. Subsequently, the events were traced to show several root cause sources such as (i) trapped energetic protons leaking to low altitudes within the South Atlantic Anomaly, (ii) Solar Energetic Particles emitted by an impulsive event on 10 September 2017, and (iii) Galactic Cosmic Rays. A profound analysis was carried out on the observed flight data, and its corresponding results are actually in agreement with the standard energetic particle models. The presented results provide another important insight on the Single Event Upsets for future Low Earth Orbit satellite missions.     

Analysis of electromagnetic and electrostatic effects of particle impacts on spacecraft

Particle impacts on spacecraft can cause considerable damage, even leading to complete failure. A theory for the resulting vehicle potential changes and the electromagnetic radiation from impact-induced plasma has been published by Close et al. (2010). Here we compare this theory to impacts registered by the Radio and Plasma Wave Science instrumentation on the Cassini spacecraft. We study both low-velocity (16 km/s) large particles (2.6 μm radius) detected in Saturn’s rings and high-velocity (450 km/s) small particles (1 nm radius) in the solar wind. The agreement with the theory is quite good. We also apply these results to earth orbit and conclude that both Electrostatic Discharge and Electromagnetic Pulse radiation are likely and could lead to spacecraft failure.     

On some large solar energetic particle events: Direct impulsive component and broad associated coronal mass ejections

Using the proton intensity and X-ray flux data from the GOES, combined with the observations of the associated solar eruptions by the Large Angle and Spectrometric Coronagraph Experiment (LASCO) on board the Solar and Heliospheric Observatory (SOHO), 14 large SEP events occurring in the period 2000 January–2002 April have been studied. It is found that: (1) events with the SEPs increasing shortly after the maximum of their parent flares (<1 h; hereafter prompt events) have rapid and great (up to four orders of magnitude) SEP increments in high-energy channels (> ∼100 MeV); however, for events whose onset of the SEP injection lags the flare maximum for a long time (>3 h; hereafter delayed events), the high-energy SEPs show no obvious enhancements (within one order of magnitude); (2) peak intensity of the prompt events is distinctly larger than that of the delayed events; (3) CMEs associated with the poorly magnetically connected events (source region <W30°) in our survey are all halo CMEs. From these observational differences, we propose a special scenario of the production of the largest SEP events: both CMEs and flares are induced in the same coronal process; high-energy particles accelerated in the reconnection region can escape easily from the open field lines and/or be transported by fast CMEs into interplanetary space, indicating a direct impulsive component in large gradual SEP events. Meanwhile, the broad width of the associated CMEs implies that the CME width is more important in SEP events production than previously considered.     

Mice heterozygous for the ATM gene are more sensitive to both X-ray and heavy ion exposure than are wildtypes

Previous studies have shown that the eyes of ATM heterozygous mice exposed to low-LET radiation (X-rays) are significantly more susceptible to the development of cataracts than are those of wildtype mice. The findings, as well as others, run counter to the assumption underpinning current radiation safety guidelines, that individuals are all equally sensitive to the biological effects of radiation. A question, highly relevant to human space activities is whether or not, in similar fashion there may exist a genetic predisposition to high-LET radiation damage.Mice haplodeficient for the ATM gene and wildtypes were exposed to 325 mGy of 1 GeV/amu ⁵⁶Fe ions at the AGS facility of Brookhaven National Laboratory. The fluence was equivalent to 1 ion per lens epithelial cell nuclear area. Controls consisted of irradiated wildtype as well as unirradiated wildtype and heterozygous mice. Prevalence analyses for stage 0.5–3.0 cataracts indicated that not only cataract onset but also progression were accelerated in the mice haplo-deficient for the ATM gene.The data show that heterozygosity for the ATM gene predisposes the eye to the cataractogenic influence of heavy ions and suggest that ATM heterozygotes in the human population may also be radiosensitive. This may have to be considered in the selection of individuals who will be exposed to both HZE particles and low-LET radiation as they may be predisposed to increased late normal tissue damage.