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.     

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.     

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.     

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.     

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.     

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.     

Heavy ion contributions to organ dose equivalent for the 1977 galactic cosmic ray spectrum

Estimates of organ dose equivalents for the skin, eye lens, blood forming organs, central nervous system, and heart of female astronauts from exposures to the 1977 solar minimum galactic cosmic radiation spectrum for various shielding geometries involving simple spheres and locations within the Space Transportation System (space shuttle) and the International Space Station (ISS) are made using the HZETRN 2010 space radiation transport code. The dose equivalent contributions are broken down by charge groups in order to better understand the sources of the exposures to these organs. For thin shields, contributions from ions heavier than alpha particles comprise at least half of the organ dose equivalent. For thick shields, such as the ISS locations, heavy ions contribute less than 30% and in some cases less than 10% of the organ dose equivalent. Secondary neutron production contributions in thick shields also tend to be as large, or larger, than the heavy ion contributions to the organ dose equivalents.     

A new dynamical atmospheric ionizing radiation (AIR) model for epidemiological studies.

A new Atmospheric Ionizing Radiation (AIR) model is currently being developed for use in radiation dose evaluation in epidemiological studies targeted to atmospheric flight personnel such as civilian airlines crewmembers. The model will allow computing values for biologically relevant parameters, e.g. dose equivalent and effective dose, for individual flights from 1945. Each flight is described by its actual three dimensional flight profile, i.e. geographic coordinates and altitudes varying with time. Solar modulated primary particles are filtered with a new analytical fully angular dependent geomagnetic cut off rigidity model, as a function of latitude, longitude, arrival direction, altitude and time. The particle transport results have been obtained with a technique based on the three-dimensional Monte Carlo transport code FLUKA, with a special procedure to deal with HZE particles. Particle fluxes are transformed into dose-related quantities and then integrated all along the flight path to obtain the overall flight dose. Preliminary validations of the particle transport technique using data from the AIR Project ER-2 flight campaign of measurements are encouraging. Future efforts will deal with modeling of the effects of the aircraft structure as well as inclusion of solar particle events.     

Galactic cosmic ray radiation levels in spacecraft on interplanetary missions.

Using the Langley Research Center galactic cosmic ray (GCR) transport computer code (HZETRN) and the computerized anatomical man (CAM) model, crew radiation levels inside manned spacecraft on interplanetary missions are estimated. These radiation-level estimates include particle fluxes, LET (linear energy transfer) spectra, absorbed dose, and dose equivalent within various organs of interest in GCR protection studies. Changes in these radiation levels resulting from the use of various different types of shield materials are presented.     

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.     

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.     

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.     

Galactic cosmic rays and cell-hit frequencies outside the magnetosphere.

An evaluation of the exposure of space travelers to galactic cosmic radiation outside the earth's magnetosphere is made by calculating fluences of high-energy primary and secondary particles with various charges traversing a sphere of area 100 microns2. Calculations relating to two shielding configurations are presented: the center of a spherical aluminum shell of thickness 1 g/cm2, and the center of a 4 g/cm2 thick aluminum spherical shell within which there is a 30 g/cm2 diameter spherical water phantom with the point of interest 5 g/cm2 from the surface. The area of 100 microns2 was chosen to simulate the nucleus of a cell in the body. The frequencies as a function of charge component in both shielding configurations reflects the odd-even disparity of the incident particle abundances. For a three-year mission, 33% of the cells in the more heavily shielded configuration would be hit by at least one particle with Z greater than 10. Six percent would be hit by at least two such particles. This emphasizes the importance of studying single high-Z particle effects both on cells which might be "at risk" for cancer induction and on critical neural cells or networks which might be vulnerable to inactivation by heavy charged particle tracks. Synergistic effects with the more numerous high-energy protons and helium ions cannot be ruled out. In terms of more conventional radiation risk assessment, the dose equivalent decreased by a factor of 2.85 from free space to that in the more heavily shielded configuration. Roughly half of this was due to the decrease in energy deposition (absorbed dose) and half to the decrease in biological effectiveness (quality factor).     

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.     

Analysis of radiation risk from alpha particle component of solar particle events.

The solar particle events (SPE) will contain a primary alpha particle component, representing a possible increase in the potential risk to astronauts during an SPE over the often studied proton component. We discuss the physical interactions of alpha particles important in describing the transport of these particles through spacecraft and body shielding. Models of light ion reactions are presented and their effects on energy and linear energy transfer (LET) spectra in shielding discussed. We present predictions of particle spectra, dose, and dose equivalent in organs of interest for SPE spectra typical of those occurring in recent solar cycles. The large events of solar cycle 19 are found to have substantial increase in biological risk from alpha particles, including a large increase in secondary neutron production from alpha particle breakup.     

Role of shielding in modulating the effects of solar particle events: Monte Carlo calculation of absorbed dose and DNA complex lesions in different organs.

Distributions of absorbed dose and DNA clustered damage yields in various organs and tissues following the October 1989 solar particle event (SPE) were calculated by coupling the FLUKA Monte Carlo transport code with two anthropomorphic phantoms (a mathematical model and a voxel model), with the main aim of quantifying the role of the shielding features in modulating organ doses. The phantoms, which were assumed to be in deep space, were inserted into a shielding box of variable thickness and material and were irradiated with the proton spectra of the October 1989 event. Average numbers of DNA lesions per cell in different organs were calculated by adopting a technique already tested in previous works, consisting of integrating into "condensed-history" Monte Carlo transport codes--such as FLUKA--yields of radiobiological damage, either calculated with "event-by-event" track structure simulations, or taken from experimental works available in the literature. More specifically, the yields of "Complex Lesions" (or "CL", defined and calculated as a clustered DNA damage in a previous work) per unit dose and DNA mass (CL Gy-1 Da-1) due to the various beam components, including those derived from nuclear interactions with the shielding and the human body, were integrated in FLUKA. This provided spatial distributions of CL/cell yields in different organs, as well as distributions of absorbed doses. The contributions of primary protons and secondary hadrons were calculated separately, and the simulations were repeated for values of Al shielding thickness ranging between 1 and 20 g/cm2. Slight differences were found between the two phantom types. Skin and eye lenses were found to receive larger doses with respect to internal organs; however, shielding was more effective for skin and lenses. Secondary particles arising from nuclear interactions were found to have a minor role, although their relative contribution was found to be larger for the Complex Lesions than for the absorbed dose, due to their higher LET and thus higher biological effectiveness.     

Radiation hazards on space missions outside the magnetosphere.

Future space missions outside the magnetosphere will subject astronauts to a hostile and unfamiliar radiation environment. An annual dose equivalent to the blood-forming organs (BFOs) of approximately 0.5 Sv is expected, mostly from heavy ions in the galactic cosmic radiation. On long-duration missions, an anomalously-large solar energetic particle event may occur. Such an event can expose astronauts to up to approximately 25 Gy (skin dose) and up to approximately 2 Sv (BFO dose) with no shielding. The anticipated radiation exposure may necessitate spacecraft design concessions and some restriction of mission activities. In this paper we discuss our model calculations of radiation doses in several exo-magnetospheric environments. Specific radiation shielding strategies are discussed. A new calculation of aluminum equivalents of potential spacecraft shielding materials demonstrates the importance of low-atomic-mass species for protection from galactic cosmic radiation.     

Estimates of HZE particle contributions to SPE radiation exposures on interplanetary missions.

Estimates of radiation doses resulting from possible HZE (high energy heavy ion) components of solar particle events (SPEs) are presented for crews of manned interplanetary missions. The calculations assume a model spectrum obtained by folding measured solar flare HZE particle abundances with the measured energy spectra of SPE alpha particles. These hypothetical spectra are then transported through aluminum spacecraft shielding. The results, presented as estimates of absorbed dose and dose equivalent, indicate that HZE components by themselves are not a major concern for crew protection but should be included in any overall risk assessment. The predictions are found to be sensitive to the assumed spectral hardness parameters.