Perspective on space radiation for space flights in 2020–2040

The Sun undergoes several well known periodicities in activity, such as the Schwabe 11 year cycle, the Gleissberg 80–90 year cycle, the Suess 200–210 year cycle and the Halstatt 2200–2300 year cycle. In addition, there is evidence that the 20th century levels of solar activity are unusually high. The years 2020–2040 are expected to coincide with increased activity in human space flight beyond low Earth orbit. The solar cycles and the present level of solar activity are reviewed and their activities during the years 2020–2040 are discussed with a perspective on space radiation and the future program of space flight. It is prudent to prepare for continuing levels of high solar activity as well as for the low levels of the current deep minimum, which has corresponded to high galactic cosmic ray flux.     

Comparison of the transport codes HZETRN,HETC and FLUKA for galactic cosmic rays

The HZETRN deterministic radiation code is one of several tools developed to analyze the effects of harmful galactic cosmic rays (GCR) and solar particle events on mission planning and shielding for astronauts and instrumentation. This paper is a comparison study involving the two Monte Carlo transport codes, HETC–HEDS and FLUKA and the deterministic transport code, HZETRN. Each code is used to transport an ion from the 1977 solar minimum GCR spectrum impinging upon a 20 g/cm² aluminum slab followed by a 30 g/cm² water slab. This research is part of a systematic effort of verification and validation to quantify the accuracy of HZETRN and determine areas where it can be improved. Comparisons of dose and dose equivalent values at various depths in the water slab are presented in this report. This is followed by a comparison of the proton and forward, backward and total neutron flux at various depths in the water slab. Comparisons of the secondary light ion ²H, ³H, ³He and ⁴He fluxes are also examined.     

An update about recent developments of the PHITS code

PHITS (Particle and Heavy-Ion Transport code System) is a general-purpose three-dimensional Monte Carlo code, developed and maintained by RIST, JAEA and KEK in Japan together with Sihver et al. at Chalmers in Sweden. PHITS can deal with the transports of all varieties of hadrons and heavy ions with energies up to around 100 GeV/nucleon, and in this paper the current status of PHITS is presented. We introduce a relativistically covariant version of JQMD, called R-JQMD, that features an improved ground state initialization algorithm, and we will present the introduction of electron and photon transport in PHITS using EGS5, which have increased the energy region for the photon and energy transport from up to around 3 GeV to up to several hundred GeV depending on the atomic number of the target. We show how the accuracy in dose and fluence calculations can be improved by using tabulated cross sections. Benchmarking of shielding and irradiation effects of high energy protons in different materials relevant for shielding of accelerator facilities is also presented. In particular, we show that PHITS can be used for estimating the dose received by aircrews and personnel in space. In recent years, many countries have issued regulations or recommendations to set annual dose limitations for aircrews. Since estimation of cosmic-ray spectra in the atmosphere is an essential issue for the evaluation of aviation doses, we have calculated these spectra using PHITS. The accuracy of the atmospheric propagation simulation of cosmic-ray performed by PHITS has been well verified by experimental cosmic-ray spectra taken under various conditions. Based on a comprehensive analysis of the simulation results, an analytical model called “PARMA” has been proposed for instantaneously estimating the atmospheric cosmic-ray spectra below the altitude of 20 km. We have also performed preliminary simulations of long-term dose distribution measurements at the ISS performed with the joint ESA-FSA experiment MATROSHKA-R (MTR-R) led by the Russian Federation Institute of Biomedical Problems (IMBP) and the ESA supported experiment MATROSHKA (MTR), led by the German Aerospace Center (DLR). For the purpose of examining the applicability of PHITS to the shielding design in space, the absorbed doses in a tissue equivalent water phantom inside an imaginary space vessel has been estimated for different shielding materials of different thicknesses. The results confirm previous results which indicate that PHITS is a suitable tool when performing shielding design studies of spacecrafts.     

Electrostatic space radiation shielding

For the success of NASA’s new vision for space exploration to Moon, Mars and beyond, exposures from the hazards of severe space radiation in deep space long duration missions is ‘a must solve’ problem. The payload penalty demands a very stringent requirement on the design of the spacecrafts for human deep space missions. The exploration beyond low Earth orbit (LEO) to enable routine access of space will require protection from the hazards of the accumulated exposures of space radiation, Galactic Cosmic Rays (GCR) and Solar Particle Events (SPE), and minimizing the production of secondary radiation is a great advantage. There is a need to look to new horizons for newer technologies. The present investigation revisits electrostatic active radiation shielding and explores the feasibility of using the electrostatic shielding in concert with the state-of-the-art materials shielding and protection technologies. The full space radiation environment has been used, for the first time, to explore the feasibility of electrostatic shielding. The goal is to repel enough positive charge ions so that they miss the spacecraft without attracting thermal electrons. Conclusions are drawn for the future directions of space radiation protection.     

Pion and electromagnetic contribution to dose: Comparisons of HZETRN to Monte Carlo results and ISS data

Recent work has indicated that pion production and the associated electromagnetic (EM) cascade may be an important contribution to the total astronaut exposure in space. Recent extensions to the deterministic space radiation transport code, HZETRN, allow the production and transport of pions, muons, electrons, positrons, and photons. In this paper, the extended code is compared to the Monte Carlo codes, Geant4, PHITS, and FLUKA, in slab geometries exposed to galactic cosmic ray (GCR) boundary conditions. While improvements in the HZETRN transport formalism for the new particles are needed, it is shown that reasonable agreement on dose is found at larger shielding thicknesses commonly found on the International Space Station (ISS). Finally, the extended code is compared to ISS data on a minute-by-minute basis over a seven day period in 2001. The impact of pion/EM production on exposure estimates and validation results is clearly shown. The Badhwar–O’Neill (BO) 2004 and 2010 models are used to generate the GCR boundary condition at each time-step allowing the impact of environmental model improvements on validation results to be quantified as well. It is found that the updated BO2010 model noticeably reduces overall exposure estimates from the BO2004 model, and the additional production mechanisms in HZETRN provide some compensation. It is shown that the overestimates provided by the BO2004 GCR model in previous validation studies led to deflated uncertainty estimates for environmental, physics, and transport models, and allowed an important physical interaction (π/EM) to be overlooked in model development. Despite the additional π/EM production mechanisms in HZETRN, a systematic under-prediction of total dose is observed in comparison to Monte Carlo results and measured data.     

Comparison of the transport codes HZETRN,HETC and FLUKA for a solar particle event

The protection of astronauts and instrumentation from galactic cosmic rays and solar particle events is one of the primary constraints associated with mission planning in low earth orbit or deep space. To help satisfy this constraint, several computational tools have been developed to analyze the effectiveness of various shielding materials and structures exposed to space radiation. These tools are now being carefully scrutinized through a systematic effort of verification, validation, and uncertainty quantification. In this benchmark study, the deterministic transport code HZETRN is compared to the Monte Carlo transport codes HETC-HEDS and FLUKA for a 30 g/cm² water target protected by a 20 g/cm² aluminum shield exposed to a parameterization of the February 1956 solar particle event. Neutron and proton fluences as well as dose and dose equivalent are compared at various depths in the water target. The regions of agreement and disagreement between the three codes are quantified and discussed, and recommendations for future work are given.     

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.     

空间辐射场蒙特卡罗模拟计算方法研究

对利用蒙特卡罗方法对由银河宇宙射线引起的空间辐射场各成分进行计算的方法进行了调研,对计算模型的建立以及计算过程中通常使用的方差减小技术进行分析,给出了美国的Roesler等人利用FLUKA程序以及加拿大Anid等人利用MCNPX程序计算得到的由银河宇宙射线引起的空间辐射场各量值及其与实验结果的比较,验证了计算方法与计算模型的可靠性。对任意航线空间辐射场剂量分布预评估方法进行分析,给出了由银河宇宙射线引起的空间辐射场的基本特征。     

A space radiation shielding model of the Martian radiation environment experiment (MARIE)

The 2001 Mars Odyssey spacecraft was launched towards Mars on April 7, 2001. Onboard the spacecraft is the Martian radiation environment experiment (MARIE), which is designed to measure the background radiation environment due to galactic cosmic rays (GCR) and solar protons in the 20–500 MeV/n energy range. We present an approach for developing a space radiation-shielding model of the spacecraft that includes the MARIE instrument in the current mapping phase orientation. A discussion is presented describing the development and methodology used to construct the shielding model. For a given GCR model environment, using the current MARIE shielding model and the high-energy particle transport codes, dose rate values are compared with MARIE measurements during the early mapping phase in Mars orbit. The results show good agreement between the model calculations and the MARIE measurements as presented for the March 2002 dataset.     

Symplectic integration of space debris motion considering several Earth’s shadowing models

In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones. Indeed, it has already been demonstrated that such objects could stay there for hundreds of years. Our model takes into account the Earth’s gravitational potential, luni-solar and planetary gravitational perturbations and direct solar radiation pressure. Based on the analysis of the energy conservation and on a comparison with a high order non-symplectic integrator, we show that our algorithm allows us to use large time steps and keep accurate results. We also propose an innovative method to model Earth’s shadow crossings by means of a smooth shadow function. In the particular framework of symplectic integration, such a function needs to be included analytically in the equations of motion in order to prevent numerical drifts of the energy. For the sake of completeness, both cylindrical shadows and penumbra transitions models are considered. We show that both models are not equivalent and that big discrepancies actually appear between associated orbits, especially for high area-to-mass ratios.     

Deterministic pion and muon transport in Earth’s atmosphere

An accurate understanding of the physical interactions and transport of space radiation is important for safe and efficient space operations. Secondary particles produced by primary particle interactions with intervening materials are an important contribution to radiation risk. Pions are copiously produced in the nuclear interactions typical of space radiations and can therefore be an important contribution to radiation exposure. Charged pions decay almost exclusively to muons. As a consequence, muons must also be considered in space radiation exposure studies. In this work, the NASA space radiation transport code HZETRN has been extended to include the transport of charged pions and muons. The relevant transport equation, solution method, and implemented cross sections are reviewed. Muon production in the Earth’s upper atmosphere is then investigated, and comparisons with recent balloon flight measurements of differential muon flux are presented. Muon production from the updated version of HZETRN is found to match the experimental data well.     

Comparison of organ dose and dose equivalent for human phantoms of CAM vs. MAX

For the evaluation of organ dose and dose equivalent of astronauts on space shuttle and the International Space Station (ISS) missions, the CAMERA models of CAM (Computerized Anatomical Male) and CAF (Computerized Anatomical Female) of human tissue shielding have been implemented and used in radiation transport model calculations at NASA. One of new human geometry models to meet the “reference person” of International Commission on Radiological Protection (ICRP) is based on detailed Voxel (volumetric and pixel) phantom models denoted for male and female as MAX (Male Adult voXel) and FAX (Female Adult voXel), respectively. We compared the CAM model predictions of organ doses to those of MAX model, since the MAX model represents the male adult body with much higher fidelity than the CAM model currently used at NASA. Directional body-shielding mass was evaluated for over 1500 target points of MAX for specified organs considered to be sensitive to the induction of stochastic effects. Radiation exposures to solar particle event (SPE), trapped protons, and galactic cosmic ray (GCR) were assessed at the specific sites in the MAX phantom by coupling space radiation transport models with the relevant body-shielding mass. The development of multiple-point body-shielding distributions at each organ made it possible to estimate the mean and variance of organ doses at the specific organ. For the estimate of doses to the blood forming organs (BFOs), data on active marrow distributions in adult were used to weight the bone marrow sites over the human body. The discrete number of target points of MAX organs resulted in a reduced organ dose and dose equivalent compared to the results of CAM organs especially for SPE, and should be further investigated. Differences of effective doses between the two approaches were found to be small (<5%) for GCR.     

Using high-energy proton fluence to improve risk prediction for consequences of solar particle events

The potential for exposure to large solar particle events (SPEs) with high energy levels is a major concern during interplanetary transfer and extra-vehicular activities (EVAs) on the lunar and Mars surface. Previously, we have used data from the last 5 solar cycles to estimate percentiles of dose to a typical blood-forming organ (BFO) for a hypothetical astronaut in a nominally shielded spacecraft during a 120-d lunar mission. As part of this process, we made use of complete energy spectra for 34 large historical SPEs to calculate what the BFO mGy-Eq dose would have been in the above lunar scenario for each SPE. From these calculated doses, we then developed a prediction model for BFO dose based solely on an assumed value of integrated fluence above 30 MeV (Φ₃₀) for an otherwise unspecified future SPE. In this study, we reasoned that since BFO dose is determined more by protons with higher energies than by those with lower energies, more accurate BFO dose prediction models could be developed using integrated fluence above 60 (Φ₆₀) and above 100 MeV (Φ₁₀₀) as predictors instead of Φ₃₀. However to calculate the unconditional probability of a BFO dose exceeding a pre-specified limit (“BFO dose risk”), one must also take into account the distribution of the predictor (Φ_30,Φ₆₀, or Φ₁₀₀), as estimated from historical SPEs. But Φ₆₀ and Φ₁₀₀ have more variability, and less available historical information on which to estimate their distributions over many SPE occurrences, than does Φ₃₀. Therefore, when estimating BFO dose risk there is a tradeoff between increased BFO dose prediction at a given energy threshold and decreased accuracy of models for describing the distribution of that threshold over future SPEs as the threshold increases. Even when taking the second of these two factors into account, we still arrived at the conclusion that overall prediction improves as the energy level threshold increases from 30 to 60 to 100 MeV. These results can be applied to the development of approaches to improve radiation protection of astronauts and the optimization of mission planning for future space missions.     

Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA’s radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ∼15 Tm and passive/structural shielding thicknesses greater than 40 g/cm² have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.     

Statistical validation of HZETRN as a function of vertical cutoff rigidity using ISS measurements

Measurements taken in Low Earth Orbit (LEO) onboard the International Space Station (ISS) and transit vehicles have been extensively used to validate radiation transport models. Primarily, such comparisons were done by integrating measured data over mission or trajectory segments so that individual comparisons to model results could be made. This approach has yielded considerable information but is limited in its ability to rigorously quantify and differentiate specific model errors or uncertainties. Further, as exploration moves beyond LEO and measured data become sparse, the uncertainty estimates derived from these validation cases will no longer be applicable. Recent improvements in the underlying numerical methods used in HZETRN have resulted in significant decreases in code run time. Therefore, the large number of comparisons required to express error as a function of a physical quantity, like cutoff rigidity, are now possible. Validation can be looked at in detail over any portion of a flight trajectory (e.g. minute by minute) such that a statistically significant number of comparisons can be made. This more rigorous approach to code validation will allow the errors caused by uncertainties in the geometry models, environmental models, and nuclear physics models to be differentiated and quantified. It will also give much better guidance for future model development. More importantly, it will allow a quantitative means of extrapolating uncertainties in LEO to free space. In this work, measured data taken onboard the ISS during solar maximum are compared to results obtained with the particle transport code HZETRN. Comparisons are made at a large number (∼77,000) of discrete time intervals, allowing error estimates to be given as a function of cutoff rigidity. It is shown that HZETRN systematically underestimates exposure quantities at high cutoff rigidity. The errors are likely associated with increased angular variation in the geomagnetic field near the equator, the lack of pion production in HZETRN, and errors in high energy nuclear physics models, and will be the focus of future work.     

How Galactic Cosmic Ray models affect the estimation of radiation exposure in space

The radiation environment in space is a major concern for human spaceflight because of the adverse effects of high levels of radiation on astronauts’ health. Therefore, it is essential to perform radiation risk assessments already during the concept studies of a manned mission. Galactic Cosmic Rays (GCR) have been identified to be one of the primary sources of radiation exposure in space.     

Human response to high-background radiation environments on Earth and in space

The main long-term objective of the space exploration program is the colonization of the planets of the Solar System. The high cosmic radiation equivalent dose rate represents an inescapable problem for the safe establishment of permanent human settlements on these planets. The unshielded equivalent dose rate on Mars ranges between 100 and 200 mSv/year, depending on the Solar cycle and altitude, and can reach values as high as 360 mSv/year on the Moon. The average annual effective dose on Earth is about 3 mSv, nearly 85% of which comes from natural background radiation, reduced to less than 1 mSv if man-made sources and the internal exposure to Rn daughters are excluded. However, some areas on Earth display anomalously high levels of background radiation, as is the case with thorium-rich monazite bearing sand deposits where values 200–400 times higher than the world average can be found. About 2% of the world’s population live above 3 km and receive a disproportionate 10% of the annual effective collective dose due to cosmic radiation, with a net contribution to effective dose by the neutron component which is 3–4 fold that at sea level. Thus far, epidemiological studies have failed to show any adverse health effects in the populations living in these terrestrial high-background radiation areas (HBRA), which provide an unique opportunity to study the health implications of an environment that, as closely as possibly achievable on Earth, resembles the chronic exposure of future space colonists to higher-than-normal levels of ionizing radiation. Chromosomal aberrations in the peripheral blood lymphocytes from the HBRA residents have been measured in several studies because chromosomal damage represents an early biomarker of cancer risk. Similar cytogenetic studies have been recently performed in a cohort of astronauts involved in single or repeated space flights over many years. The cytogenetic findings in populations exposed to high dose-rate background radiation on Earth or in space will be discussed.     

The effects of proton radiation on UHMWPE material properties for space flight and medical applications

Ultra High Molecular Weight Polyethylene (UHMWPE) is a polymer widely used as a radiation shielding material in space flight applications and as a bearing material in total joint replacements. As a long chain hydrocarbon based polymer, UHMWPE’s material properties are influenced by radiation exposure, and prior studies show that gamma irradiation is effective for both medical sterilization and increased wear resistance in total joint replacement applications. However, the effects of space flight radiation types and doses on UHMWPE material properties are poorly understood. In this study, three clinically relevant grades of UHMWPE (GUR 1020, GUR 1050, and GUR 1020 blended with Vitamin E) were proton irradiated and tested for differences in material properties. Each of the three types of UHMWPE was irradiated at nominal doses of 0 Gy (control), 5 Gy, 10 Gy, 20 Gy, and 35 Gy. Following irradiation, uniaxial tensile testing and thermal testing using Differential Scanning Calorimetry (DSC) and Dynamic Mechanical Analysis (DMA) were performed. Results show small but significant changes in several material properties between the control (0 Gy) and 35 Gy samples, indicating that proton irradiation could have a effect on the long term performance of UHMWPE in both medical and space flight applications.