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.     

Utilization of CAM,CAF, MAX,and FAX for space radiation analyses using HZETRN

To estimate astronaut health risk due to space radiation, one must have the ability to calculate various exposure-related quantities that are averaged over specific organs and tissue types. Such calculations require computational models of the ambient space radiation environment, particle transport, nuclear and atomic physics, and the human body. While significant efforts have been made to verify, validate, and quantify the uncertainties associated with many of these models and tools, relatively little work has focused on the uncertainties associated with the representation and utilization of the human phantoms. In this study, we first examine the anatomical properties of the Computerized Anatomical Man (CAM), Computerized Anatomical Female (CAF), Male Adult voXel (MAX), and Female Adult voXel (FAX) models by comparing the masses of various model tissues used to calculate effective dose to the reference values specified by the International Commission on Radiological Protection (ICRP). The MAX and FAX tissue masses are found to be in good agreement with the reference data, while major discrepancies are found between the CAM and CAF tissue masses and the reference data for almost all of the effective dose tissues. We next examine the distribution of target points used with the deterministic transport code HZETRN (High charge (Z) and Energy TRaNsport) to compute mass averaged exposure quantities. A numerical algorithm is presented and used to generate multiple point distributions of varying fidelity for many of the effective dose tissues identified in CAM, CAF, MAX, and FAX. The point distributions are used to compute mass averaged dose equivalent values under both a galactic cosmic ray (GCR) and solar particle event (SPE) environment impinging isotropically on three spherical aluminum shells with areal densities of 0.4 g/cm², 2.0 g/cm², and 10.0 g/cm². The dose equivalent values are examined to identify a recommended set of target points for each of the tissues and to further assess the differences between CAM, CAF, MAX, and FAX. It is concluded that the previously published CAM and CAF point distributions were significantly under-sampled and that the set of point distributions presented here should be adequate for future studies involving CAM, CAF, MAX, or FAX. It is also found that the errors associated with the mass and location of certain tissues in CAM and CAF have a significant impact on the mass averaged dose equivalent values, and it is concluded that MAX and FAX are more accurate than CAM and CAF for space radiation analyses.     

Simulations of absorbed dose on the phantom surface of MATROSHKA-R experiment at the ISS

The health risks associated with exposure to various components of space radiation are of great concern when planning manned long-term interplanetary missions, such as future missions to Mars. Since it is not possible to measure the radiation environment inside of human organs in deep space, simulations based on radiation transport/interaction codes coupled to phantoms of tissue equivalent materials are used. However, the calculated results depend on the models used in the codes, and it is therefore necessary to verify their validity by comparison with measured data. The goal of this paper is to compare absorbed doses obtained in the MATROSHKA-R experiment performed at the International Space Station (ISS) with simulations performed with the three-dimensional Monte Carlo Particle and Heavy-Ion Transport code System (PHITS). The absorbed dose was measured using passive detectors (packages of thermoluminescent and plastic nuclear track detectors) placed on the surface of the spherical tissue equivalent phantom MATROSHKA-R, which was exposed aboard the ISS in the Service Zvezda Module from December 2005 to September 2006. The data calculated by PHITS assuming an ISS shielding of 3 g/cm² and 5 g/cm² aluminum mass thickness were in good agreement with the measurements. Using a simplified geometrical model of the ISS, the influence of variations in altitude and wall mass thickness of the ISS on the calculated absorbed dose was estimated. The uncertainties of the calculated data are also discussed; the relative expanded uncertainty of absorbed dose in phantom was estimated to be 44% at a 95% confidence level.     

Estimation of Galactic Cosmic Ray exposure inside and outside the Earth’s magnetosphere during the recent solar minimum between solar cycles 23 and 24

The evidently low solar activity observed between solar cycles 23 and 24 during the years 2008–2010 led to a substantial increase in the Galactic Cosmic Ray (GCR) intensity in comparison with preceding solar minima. As the GCRs consist of highly-ionizing charged particles having the potential to cause biological damage, they are a subject of concern for manned missions to space. With the enhanced particle fluxes observed between 2008 and 2010, it is reasonable to assume that the radiation exposure from GCR must have also increased to unusually high levels. In this paper, the GCR exposure outside and inside the Earth’s magnetosphere is numerically calculated for time periods starting from 1970 to the end of 2011 in order to investigate the increase in dose levels during the years 2008–2010 in comparison with the last three solar minima. The dose rates were calculated in a water sphere, used as a surrogate for the human body, either unshielded or surrounded by aluminium shielding of 0.3, 10 or 40 g/cm².     

Comparisons of several transport models in their predictions in typical space radiation environments

We have used several transport codes to calculate dose and dose equivalent values as well as the particle spectra behind a slab or inside a spherical shell shielding in typical space radiation environments. Two deterministic codes, HZETRN and UPROP, and two Monte Carlo codes, FLUKA and Geant4, are included. A soft solar particle event, a hard solar particle event, and a solar minimum galactic cosmic rays environment are considered; and the shielding material is either aluminum or polyethylene. We find that the dose values and particle spectra from HZETRN are in general rather consistent with Geant4 except for neutrons. The dose equivalent values from HZETRN and Geant4 are not far from each other, but the HZETRN values behind shielding are often lower than the Geant4 values. Results from FLUKA and Geant4 are mostly consistent for considered cases. However, results from the legacy code UPROP are often quite different from the other transport codes, partly due to its non-consideration of neutrons. Comparisons for the spherical shell geometry exhibit the same qualitative features as for the slab geometry. In addition, results from both deterministic and Monte Carlo transport codes show that the dose equivalent inside the spherical shell decreases from the center to the inner surface and this decrease is large for solar particle events; consistent with an earlier study based on deterministic radiation transport results. This study demonstrates both the consistency and inconsistency among these transport models in their typical space radiation predictions; further studies will be required to pinpoint the exact physics modules in these models that cause the differences and thus may be improved.     

On the validity of the aluminum equivalent approximation in space radiation shielding applications

The origin of the aluminum equivalent shield approximation in space radiation analysis can be traced back to its roots in the early years of the NASA space programs (Mercury, Gemini and Apollo) wherein the primary radiobiological concern was the intense sources of ionizing radiation causing short term effects which was thought to jeopardize the safety of the crew and hence the mission. Herein, it is shown that the aluminum equivalent shield approximation, although reasonably well suited for that time period and to the application for which it was developed, is of questionable usefulness to the radiobiological concerns of routine space operations of the 21st century which will include long stays onboard the International Space Station (ISS) and perhaps the moon. This is especially true for a risk based protection system, as appears imminent for deep space exploration where the long-term effects of Galactic Cosmic Ray (GCR) exposure is of primary concern. The present analysis demonstrates that sufficiently large errors in the interior particle environment of a spacecraft result from the use of the aluminum equivalent approximation, and such approximations should be avoided in future astronaut risk estimates. In this study, the aluminum equivalent approximation is evaluated as a means for estimating the particle environment within a spacecraft structure induced by the GCR radiation field. For comparison, the two extremes of the GCR environment, the 1977 solar minimum and the 2001 solar maximum, are considered. These environments are coupled to the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport (HZETRN), which propagates the GCR spectra for elements with charges (Z) in the range 1 ? Z ? 28 (H–Ni) and secondary neutrons through selected target materials. The coupling of the GCR extremes to HZETRN allows for the examination of the induced environment within the interior of an idealized spacecraft as approximated by a spherical shell shield, and the effects of the aluminum equivalent approximation for a good polymeric shield material such as generic polyethylene (PE). The shield thickness is represented by a 25 g/cm² spherical shell. Although, one could imagine the progression to greater thickness, the current range will be sufficient to evaluate the qualitative usefulness of the aluminum equivalent approximation. Upon establishing the inaccuracies of the aluminum equivalent approximation through numerical simulations of the GCR radiation field attenuation for PE and aluminum equivalent PE spherical shells, we further present results for a limited set of commercially available, hydrogen rich, multifunctional polymeric constituents to assess the effect of the aluminum equivalent approximation on their radiation attenuation response as compared to the generic PE.     

Verification of shielding effect by the water-filled materials for space radiation in the International Space Station using passive dosimeters

The dose reduction effects for space radiation by installation of water shielding material (“protective curtain”) of a stack board consisting of the hygienic wipes and towels have been experimentally evaluated in the International Space Station by using passive dosimeters. The averaged water thickness of the protective curtain was 6.3 g/cm². The passive dosimeters consisted of a combination of thermoluminescent detectors (TLDs) and plastic nuclear track detectors (PNTDs). Totally 12 passive dosimeter packages were installed in the Russian Service Module during late 2010. Half of the packages were located at the protective curtain surface and the other half were at the crew cabin wall behind or aside the protective curtain. The mean absorbed dose and dose equivalent rates are measured to be 327 μGy/day and 821 μSv/day for the unprotected packages and 224 μGy/day and 575 μSv/day for the protected packages, respectively. The observed dose reduction rate with protective curtain was found to be 37 ± 7% in dose equivalent, which was consistent with the calculation in the spherical water phantom by PHITS. The contributions due to low and high LET particles were found to be comparable in observed dose reduction rate. The protective curtain would be effective shielding material for not only trapped particles (several 10 MeV) but also for low energy galactic cosmic rays (several 100 MeV/n). The properly utilized protective curtain will effectively reduce the radiation dose for crew living in space station and prolong long-term mission in the future.     

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.     

Method for the prediction of the effective dose equivalent to the crew of the International Space Station

This paper describes a methodology for assessing the pre-mission exposure of space crew aboard the International Space Station (ISS) in terms of an effective dose equivalent. In this approach, the PHITS Monte Carlo code was used to assess the particle transport of galactic cosmic radiation (GCR) and trapped radiation for solar maximum and minimum conditions through an aluminum shield thickness. From these predicted spectra, and using fluence-to-dose conversion factors, a scaling ratio of the effective dose equivalent rate to the ICRU ambient dose equivalent rate at a 10 mm depth was determined. Only contributions from secondary neutrons, protons, and alpha particles were considered in this analysis.     

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.     

Degradation of thermal control materials under a simulated radiative space environment

A spacecraft with a passive thermal control system utilizes various thermal control materials to maintain temperatures within safe operating limits. Materials used for spacecraft applications are exposed to harsh space environments such as ultraviolet (UV) and particle (electron, proton) irradiation and atomic oxygen (AO), undergo physical damage and thermal degradation, which must be considered for spacecraft thermal design optimization and cost effectiveness. This paper describes the effect of synergistic radiation on some of the important thermal control materials to verify the assumptions of beginning-of-life (BOL) and end-of-life (EOL) properties. Studies on the degradation in the optical properties (solar absorptance and infrared emittance) of some important thermal control materials exposed to simulated radiative geostationary space environment are discussed. The current studies are purely related to the influence of radiation on the degradation of the materials; other environmental aspects (e.g., thermal cycling) are not discussed. The thermal control materials investigated herein include different kind of second-surface mirrors, white anodizing, white paints, black paints, multilayer insulation materials, varnish coated aluminized polyimide, germanium coated polyimide, polyether ether ketone (PEEK) and poly tetra fluoro ethylene (PTFE). For this purpose, a test in the constant vacuum was performed reproducing a three year radiative space environment exposure, including ultraviolet and charged particle effects on North/South panels of a geostationary three-axis stabilized spacecraft. Reflectance spectra were measured in situ in the solar range (250–2500 nm) and the corresponding solar absorptance values were calculated. The test methodology and the degradations of the materials are discussed. The most important degradations among the low solar absorptance materials were found in the white paints whereas the rigid optical solar reflectors remained quite stable. Among the high solar absorptance elements, as such the change in the solar absorptance was very low, in particular the germanium coated polyimide was found highly stable.     

Interplanetary crew doses and dose equivalents: variations among different bone marrow and skin sites.

Previously, calculations of bone marrow dose from the large solar particle event (SPE) of July 2000 were carried out using the BRYNTRN space radiation transport code and the computerized anatomical man (CAM) model. Results indicated that the dose for a bone marrow site in the mid-thigh might be twice as large as the dose for a site in the pelvis. These large variations may be significant for space radiation protection purposes, which traditionally use an average of many (typically 33) sites throughout the body. Other organs that cover large portions of the body, such as the skin, may also exhibit similar variations with doses differing from site to site. The skin traditionally uses an average of 32 sites throughout the body. Variations also occur from site to site among the dose equivalents, which may be important in determining stochastic effects. In this work, the magnitudes of dose and dose equivalent variations from site to site are investigated. The BRYNTRN and HZETRN transport codes and the CAM model are used to estimate bone marrow and skin doses and dose equivalents as a function of position in the body for several large solar particle events and annual galactic cosmic ray spectra from throughout the space era. These position-specific results are compared with the average values usually used for radiation protection purposes. Various thicknesses of aluminum shielding, representative of nominal spacecraft, are used in the analyses.     

Cosmic ray H/H ratio measured from BESS in 2000 during solar maximum

The Balloon-borne Experiment with a Superconducting Spectrometer (BESS) was flown from Lynn Lake, Manitoba, Canada in August, 2000, during the maximum solar modulation period, with an average residual atmospheric overburden of 4.3 g/cm². Precise spectral measurements of cosmic ray hydrogen isotopes from 0.178 GeV/n to 1.334 GeV/n were made during the 28.7 h of flight. This paper presents the measured energy spectra and their ratio, ²H/¹H. The results are also compared with previous measurements and theoretical predictions.     

Human exposure to large solar particle events in space.

Whenever energetic solar protons produced by solar particle events traverse bulk matter, they undergo various nuclear and atomic collision processes which significantly alter the physical characteristics and biologically important properties of their transported radiation fields. These physical interactions and their effect on the resulting radiation field within matter are described within the context of a recently developed deterministic, coupled neutron-proton space radiation transport computer code (BRYNTRN). Using this computer code, estimates of human exposure in interplanetary space, behind nominal (2 g/cm2) and storm shelter (20 g/cm2) thicknesses of aluminum shielding, are made for the large solar proton event of August 1972. Included in these calculations are estimates of cumulative exposures to the skin, ocular lens, and bone marrow as a function of time during the event. Risk assessment in terms of absorbed dose and dose equivalent is discussed for these organs. Also presented are estimates of organ exposures for hypothetical, worst-case flare scenarios. The rate of dose equivalent accumulation places this situation in an interesting region of dose rate between the very low values of usual concern in terrestrial radiation environments and the high dose rate values prevalent in radiation therapy.     

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.     

Simulation of ASTROD I test mass charging due to solar energetic particles and interplanetary electrons

As ASTROD I travels through space, its test mass will accrue charge due to exposure of the spacecraft to high-energy particles. This test mass charge will result in Coulomb forces between the test mass and the surrounding electrodes. In earlier work, we have used the GEANT 4 toolkit to simulate charging of the ASTROD test mass due to cosmic-ray protons of energies between 0.1 and 1000 GeV at solar maximum and at solar minimum. Here we use GEANT 4 to simulate the charging process due to solar energetic particle events and interplanetary electrons. We then estimate the test mass acceleration noise due to these fluxes. The predicted charging rates range from 2247 e⁺/s to 47,055 e⁺/s, at peak intensity, for the four largest SEP events in September and October 1989. Although the noise due to charging exceeds the ASTROD I budget for the two larger events, it can be suppressed through continuous discharging. The acceleration noise during the two small events is well below the design target. The charging rate of the ASTROD I test mass due to interplanetary electrons in this simulation is about −11% of the cosmic-ray protons at solar minimum, and over −37% at solar maximum. In addition to the Monte Carlo uncertainty, an error of ±30% in the net charging rates should be added to account for uncertainties in the spectra, physics models and geometry implementations.     

Impact effects from small size meteoroids and space debris

When the impact risk from meteoroids and orbital debris is assessed the main concern is usually structural damage. With their high impact velocities of typically 10–20 km/s millimeter or centimeter sized objects can puncture pressure vessels and other walls or lead to destruction of complete subsystems or even whole spacecraft. Fortunately chances of collisions with such larger objects are small (at least at present). However, particles in the size range 1–100 μm are far more abundant than larger objects and every orbiting spacecraft will encounter them with certainty. Every solar cell (8 cm² area) of the Hubble Space Telescope encountered on average 12 impacts during its 8.25 years of space exposure. Most were from micron sized particles.     

Neuroimmune response and sleep studies after whole body irradiation with high-LET particles

In order to investigate the biological effects of galactic rays on astronaut cerebral functions after space flight, mice were exposed to different heavy ions (HZE) in whole-body conditions at doses comparable to the galactic flux: ¹²C, ¹⁶O and ²⁰Ne (95 MeV/u, at 42–76 mGy). Animals were also exposed to 42 mGy of ⁶⁰Co radiation for comparison with HZE. The neuroimmune response, evaluated by interleukin-1 (IL-1) measurement, showed that this cytokine was produced 3 h after irradiation by ¹⁶O or ⁶⁰Co. In contrast, neither ¹²C (56.7 mGy) nor ²⁰Ne (76 mGy) induced IL-1 production. However, immunohistochemical staining of ¹²C-irradiated mouse brain tissue showed 2 months later a marked inflammatory reaction in the hippocampus and a diffuse response in parenchyma. Sleep studies were realized before and after exposure to 42 mGy of ¹⁶O and 76 mGy of ²⁰Ne: only the ²⁰Ne radiation displayed a small effect. A slight decrease in paradoxical sleep, corresponding to a reduction in the number of episodes of paradoxical sleep, was manifested between 8 and 22 days after exposure. Exposure to ¹²C and ¹⁶O induced no changes either in cellularity of spleen or thymus, or in caspase 3 activity (as much as four months after irradiation). Taken together, these data indicate that the CNS could be sensitive to heavy ions and that responses to HZE impact depend on the nature of the particle, the dose threshold and the time delay to develop biological processes. Differences in responses to different HZE highlight the complex biological phenomena to which astronauts are submitted during space flight.     

A new approach to solar wind monitoring

The monitoring of solar wind parameters is a key problem of the space weather program. We are presenting a new solution of plasma parameter determination suitable for small and fast solar wind monitors. The first version will be launched during the SPECTR-R project into a highly elongated orbit with apogee ∼350,000 km. The method is based on simultaneous measurements of the total ion flux and ion integral energy spectrum by six identical Faraday cups. Three of them are dedicated to determination of the ion flow direction, the other three (equipped with control grids supplied by a retarding potential) are used for determination of the density, temperature, and speed of the plasma flow. The version under development is primarily designed for the measurements in the solar wind and tail magnetosheath, thus for velocities range from 270 to 750 km/s, temperatures from 1 to 30 eV, and densities up to 200 cm⁻³. However, the instrument design can be simply modified for measurements in other regions with a substantial portion of low-energy plasma as a subsolar magnetosheath, cusp or low-latitude boundary layer. Testing of the engineering model shows that the proposed method can provide reliable plasma parameters with a high time resolution (up to 8 Hz). The paper presents not only the method and its technical realization but it documents all advantages and peculiarities of the suggested approach.     

Acceleration of charged particles by fluctuating and steady electric fields in a X-type magnetic field

Release of stored magnetic energy via particle acceleration is a characteristic feature of astrophysical plasmas. Magnetic reconnection is one of the mechanisms for releasing energy from magnetized plasmas. Collisionless magnetic reconnection could provide both the energy release mechanism and the particle accelerator in space plasmas. Here we studied particle acceleration when fluctuating (in-time) electric fields are superposed on an static X-type magnetic field in collisionless hot solar plasma. This system is chosen to mimic the reconnective dissipation of a linear MHD disturbance. Our results are compared to particle acceleration from constant electric field superposed on an X-type magnetic field. The constant electric field configuration represents the effects of steady state magnetic reconnection. Time evolution of ion and electron distributions are obtained by numerically integrating particle trajectories. The frequencies of the electric field represent a turbulent range of waves. Depending on the frequency and amplitude of the electric field, electrons and ions are accelerated to different degrees and have energy distributions of bimodal form consisting of a lower energy part and a high energy tail. For frequencies (ω in dimensioless units) in the range 0.5 ? ω ? 1.0 a substantial fraction (20%–30%) of the proton distribution is accelerated to gamma-ray producing energies. For frequencies in the range 1 ? ω ? 100.0 the bulk of the electron distribution is accelerated to hard X-ray producing energies. The acceleration mechanism is important for solar flares and solar noise storms but it could be applicable to all collisionless astrophysical plasmas.