Research plans in Europe for radiation health hazard assessment in exploratory space missions.

The European Space Agency (ESA) is currently expanding its efforts in identifying requirements and promoting research towards optimizing radiation protection of astronauts. Space agencies use common limits for tissue (deterministic) effects on the International Space Station. However, the agencies have in place different career radiation exposure limits (for stochastic effects) for astronauts in low-Earth orbit missions. Moreover, no specific limits for interplanetary missions are issued. Harmonization of risk models and dose limits for exploratory-class missions are now operational priorities, in view of the short-term plans for international exploratory-class human missions. The purpose of this paper is to report on the activity of the ESA Topical Team on space radiation research, whose task was to identify the most pertinent research requirements for improved space radiation protection and to develop a European space radiation risk model, to contribute to the efforts to reach international consensus on dose limits for deep space. The Topical Team recommended ESA to promote the development of a space radiation risk model based on European-specific expertise in: transport codes, radiobiological modelling, risk assessment, and uncertainty analysis. The model should provide cancer and non-cancer radiation risks for crews implementing exploratory missions. ESA should then support the International Commission on Radiological Protection to harmonize international models and dose limits in deep space, and guarantee continuous support in Europe for accelerator-based research configured to improve the models and develop risk mitigation strategies.

The origin of neutron biological effectiveness as a function of energy.

The understanding of the impact of radiation quality in early and late responses of biological targets to ionizing radiation exposure necessarily grounds on the results of mechanistic studies starting from physical interactions. This is particularly true when, already at the physical stage, the radiation field is mixed, as it is the case for neutron exposure. Neutron Relative Biological Effectiveness (RBE) is energy dependent, maximal for energies ~1 MeV, varying significantly among different experiments. The aim of this work is to shed light on neutron biological effectiveness as a function of field characteristics, with a comprehensive modeling approach: this brings together transport calculations of neutrons through matter (with the code PHITS) and the predictive power of the biophysical track structure code PARTRAC in terms of DNA damage evaluation. Two different energy dependent neutron RBE models are proposed: the first is phenomenological and based only on the characterization of linear energy transfer on a microscopic scale; the second is purely ab-initio and based on the induction of complex DNA damage. Results for the two models are compared and found in good qualitative agreement with current standards for radiation protection factors, which are agreed upon on the basis of RBE data.

Radiation environment at aviation altitudes and in space.

On the Earth, protection from cosmic radiation is provided by the magnetosphere and the atmosphere, but the radiation exposure increases with increasing altitude. Aircrew and especially space crew members are therefore exposed to an increased level of ionising radiation. Dosimetry onboard aircraft and spacecraft is however complicated by the presence of neutrons and high linear energy transfer particles. Film and thermoluminescent dosimeters, routinely used for ground-based personnel, do not reliably cover the range of particle types and energies found in cosmic radiation. Further, the radiation field onboard aircraft and spacecraft is not constant; its intensity and composition change mainly with altitude, geomagnetic position and solar activity (marginally also with the aircraft/spacecraft type, number of people aboard, amount of fuel etc.). The European Union Council directive 96/29/Euroatom of 1996 specifies that aircrews that could receive dose of >1 mSv y(-1) must be evaluated. The dose evaluation is routinely performed by computer programs, e.g. CARI-6, EPCARD, SIEVERT, PCAire, JISCARD and AVIDOS. Such calculations should however be carefully verified and validated. Measurements of the radiation field in aircraft are thus of a great importance. A promising option is the long-term deployment of active detectors, e.g. silicon spectrometer Liulin, TEPC Hawk and pixel detector Timepix. Outside the Earth's protective atmosphere and magnetosphere, the environment is much harsher than at aviation altitudes. In addition to the exposure to high energetic ionising cosmic radiation, there are microgravity, lack of atmosphere, psychological and psychosocial components etc. The milieu is therefore very unfriendly for any living organism. In case of solar flares, exposures of spacecraft crews may even be lethal. In this paper, long-term measurements of the radiation environment onboard Czech aircraft performed with the Liulin since 2001, as well as measurements and simulations of dose rates on and outside the International Space Station were presented. The measured and simulated results are discussed in the context of health impact.

Contribution of indirect effects to clustered damage in DNA irradiated with protons.

Protons are the dominant particles both in galactic cosmic rays and in solar particle events and, furthermore, proton irradiation becomes increasingly used in tumour treatment. It is believed that complex DNA damage is the determining factor for the consequent cellular response to radiation. DNA plasmid pBR322 was irradiated at U120-M cyclotron with 30 MeV protons and treated with two Escherichia coli base excision repair enzymes. The yields of SSBs and DSBs were analysed using agarose gel electrophoresis. DNA has been irradiated in the presence of hydroxyl radical scavenger (coumarin-3-carboxylic acid) in order to distinguish between direct and indirect damage of the biological target. Pure scavenger solution was used as a probe for measurement of induced OH· radical yields. Experimental OH· radical yield kinetics was compared with predictions computed by two theoretical models-RADAMOL and Geant4-DNA. Both approaches use Geant4-DNA for description of physical stages of radiation action, and then each of them applies a distinct model for description of the pre-chemical and chemical stage.

Dose distribution outside the target volume for 170-MeV proton beam.

Dose delivered outside the proton field during radiotherapy can potentially lead to secondary cancer development. Measurements with a 170-MeV proton beam were performed with passive detectors (track etched detectors and thermoluminescence dosemeters) in three different depths along the Bragg curve. The measurement showed an uneven decrease of the dose outside of the beam field with local enhancements. The major contribution to the delivered dose is due to high-energy protons with linear energy transfer (LET) up to 10 keV µm(-1). However, both measurement and preliminary Monte Carlo calculation also confirmed the presence of particles with higher LET.

Simulations of the MATROSHKA experiment at the international space station using PHITS.

Concerns about the biological effects of space radiation are increasing rapidly due to the perspective of long-duration manned missions, both in relation to the International Space Station (ISS) and to manned interplanetary missions to Moon and Mars in the future. As a preparation for these long-duration space missions, it is important to ensure an excellent capability to evaluate the impact of space radiation on human health, in order to secure the safety of the astronauts/cosmonauts and minimize their risks. It is therefore necessary to measure the radiation load on the personnel both inside and outside the space vehicles and certify that organ- and tissue-equivalent doses can be simulated as accurate as possible. In this paper, simulations are presented using the three-dimensional Monte Carlo Particle and Heavy-Ion Transport code System (PHITS) (Iwase et al. in J Nucl Sci Tech 39(11):1142-1151, 2002) of long-term dose measurements performed with the European Space Agency-supported MATROSHKA (MTR) experiment (Reitz and Berger in Radiat Prot Dosim 120:442-445, 2006). MATROSHKA is an anthropomorphic phantom containing over 6,000 radiation detectors, mimicking a human head and torso. The MTR experiment, led by the German Aerospace Center (DLR), was launched in January 2004 and has measured the absorbed doses from space radiation both inside and outside the ISS. Comparisons of simulations with measurements outside the ISS are presented. The results indicate that PHITS is a suitable tool for estimation of doses received from cosmic radiation and for study of the shielding of spacecraft against cosmic radiation.

Shielding of relativistic protons.

Protons are the most abundant element in the galactic cosmic radiation, and the energy spectrum peaks around 1 GeV. Shielding of relativistic protons is therefore a key problem in the radiation protection strategy of crewmembers involved in long-term missions in deep space. Hydrogen ions were accelerated up to 1 GeV at the NASA Space Radiation Laboratory, Brookhaven National Laboratory, New York. The proton beam was also shielded with thick (about 20 g/cm2) blocks of lucite (PMMA) or aluminium (Al). We found that the dose rate was increased 40-60% by the shielding and decreased as a function of the distance along the axis. Simulations using the General-Purpose Particle and Heavy-Ion Transport code System (PHITS) show that the dose increase is mostly caused by secondary protons emitted by the target. The modified radiation field after the shield has been characterized for its biological effectiveness by measuring chromosomal aberrations in human peripheral blood lymphocytes exposed just behind the shield block, or to the direct beam, in the dose range 0.5-3 Gy. Notwithstanding the increased dose per incident proton, the fraction of aberrant cells at the same dose in the sample position was not significantly modified by the shield. The PHITS code simulations show that, albeit secondary protons are slower than incident nuclei, the LET spectrum is still contained in the low-LET range (<10 keV/microm), which explains the approximately unitary value measured for the relative biological effectiveness.

Simulations of an accelerator-based shielding experiment using the particle and heavy-ion transport code system PHITS.

In order to estimate the biological effects of HZE particles, an accurate knowledge of the physics of interaction of HZE particles is necessary. Since the heavy ion transport problem is a complex one, there is a need for both experimental and theoretical studies to develop accurate transport models. RIST and JAERI (Japan), GSI (Germany) and Chalmers (Sweden) are therefore currently developing and bench marking the General-Purpose Particle and Heavy-Ion Transport code System (PHITS), which is based on the NMTC and MCNP for nucleon/meson and neutron transport respectively, and the JAM hadron cascade model. PHITS uses JAERI Quantum Molecular Dynamics (JQMD) and the Generalized Evaporation Model (GEM) for calculations of fission and evaporation processes, a model developed at NASA Langley for calculation of total reaction cross sections, and the SPAR model for stopping power calculations. The future development of PHITS includes better parameterization in the JQMD model used for the nucleus-nucleus reactions, and improvement of the models used for calculating total reaction cross sections, and addition of routines for calculating elastic scattering of heavy ions, and inclusion of radioactivity and burn up processes. As a part of an extensive bench marking of PHITS, we have compared energy spectra of secondary neutrons created by reactions of HZE particles with different targets, with thicknesses ranging from <1 to 200 cm. We have also compared simulated and measured spatial, fluence and depth-dose distributions from different high energy heavy ion reactions. In this paper, we report simulations of an accelerator-based shielding experiment, in which a beam of 1 GeV/n Fe-ions has passed through thin slabs of polyethylene, Al, and Pb at an acceptance angle up to 4 degrees.

Fragmentation of 1 GeV/nucleon iron ions in thick targets relevant for space exploration.

We have measured charged nuclear fragments produced by 1 GeV/nucleon 56Fe ions interacting with aluminium, polyethylene and lead. These materials are relevant for assessment of radiation risk for manned space flight. The data will be presented in a form suitable for comparison with models of nuclear fragmentation and transport, including linear energy transfer (LET) spectrum, fluence for iron and fragments, event-tack- and event-dose-averaged LET, total dose and iron contribution to dose.

Fragmentation of 200 and 244 MeV/u carbon beams in thick tissue-like absorbers.

Stacks consisting of thin CR-39 sheets sandwiched between thick lucite and water absorbers were perpendicularly bombarded by 12C ions at 200 and 244 MeV/u. Track radius distributions representing the charge composition of the fragmented beams were automatically measured by a particle track analysis system. After analysis of the nuclear charge distributions, the total charge removal cross-sections and elemental production cross-sections of fragments with atomic numbers from 5 to 3, were obtained down to the lower energies (approximately 50 and 100 MeV/u, respectively). It has been found that the measured total charge removal cross-section agrees with theoretical predictions within approximately 10% and very well with previous experiments in corresponding energy regions. Two model calculations for production of B fragment are in good agreement with our measured data while a third model overestimates it by approximately 12%. Theoretical cross-sections for Be and Li fragments differ strongly among the different models and from measured values.

Transport of light ions in matter.

A recent set of light ion experiments are analyzed using the Green's function method of solving the Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 code reasonably represents the fragmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model including shell structure and short range correlations appearing as tightly bound clusters in the light ion nucleus. The most recent NUCFRG2 code is augmented with a quasielastic alpha knockout model and semiempirical adjustments (up to 30 percent in charge removal) in the fragmentation process allowing reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross sections must await the full development of a coupled channel reaction model in which shell structure and clustering can be accurately evaluated.

Calculations of depth-dose distributions, cross sections and momentum loss.

The ability to know with precision the depth-dose, dose average LET, fluence and energy distributions is of great importance in many research fields, including therapeutic and diagnostic medicine when using heavy ion beams, as well as in space research. We have therefore developed a model and a computer code for calculating these distributions when using high energy proton or heavy ion beams. In this model, we use semi-empirical total reaction and partial cross section formulas developed by us and a new prescription to take into account the energy and momentum loss of the secondary nuclei. In this paper, we will also present an empirical equation for the total inelastic 4He-p cross section, as well as the partial cross sections for the production of 3He, 3H, 2H, p and n.

Nuclear fragmentation of high-energy heavy-ion beams in water.

As a part of the physical-technical program of the heavy-ion therapy project at GSI we have investigated the nuclear fragmentation of high-energy ion beams delivered by the heavy-ion synchrotron SIS, using water as a tissue-equivalent target. For a direct comparison of fragmentation properties, beams of 10B, 12C, 14N, and 16O were produced simultaneously as secondary beams from a primary 18O beam and separated in flight by magnetic beam analysis. The Z-distributions of beam fragments produced in the water target were measured via energy loss in a large ionisation chamber and a scintillator telescope. From these data we obtained both total and partial charge-changing cross sections. In addition we have performed Bragg measurements using two parallel-plate ionization chambers and a water target of variable length. The detailed shape of the measured Bragg curves and the measured cross sections are in good agreement with model calculations based on semi-empirical formulae.