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.     

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.     

Instrumentation for investigation of the depth-dose distribution by the Liulin-5 instrument of a human phantom on the Russian segment of ISS for estimation of the radiation risk during long term space flights.

Described is the Liulin-5 experiment and instrumentation, developed for investigation of the space radiation doses depth distribution in a human phantom on the Russian Segment of the International Space Station (ISS). Liulin-5 experiment is a part of the international project MATROSHKA-R on ISS. The experiment MATROSHKA-R is aimed to study the depth dose distribution at the sites of critical organs of the human body, using models of human body-anthropomorphic and spherical tissue-equivalent phantoms. The aim of Liulin-5 experiment is long term (4-5 years) investigation of the radiation environment dynamics inside the spherical tissue-equivalent phantom, mounted in different places of the Russian Segment of ISS. Energy deposition spectra, linear energy transfer spectra, flux and dose rates for protons and the biologically-relevant heavy ion components of the galactic cosmic radiation will be measured simultaneously with near real time resolution at different depths of the phantom by a telescope of silicon detectors. Data obtained together with data from other active and passive dosimeters will be used to estimate the radiation risk to the crewmembers, verify the models of radiation environment in low Earth orbit, validate body transport model and correlate organ level dose to skin dose. Presented are the test results of the prototype unit. The spherical phantom will be flown on the ISS in 2004 year and Liulin-5 experiment is planned for 2005 year.     

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.     

Assessment and requirements of nuclear reaction databases for GCR transport in the atmosphere and structures

The transport properties of galactic cosmic rays (GCR) in the atmosphere, material structures, and human body (self-shielding) are of interest in risk assessment for supersonic and subsonic aircraft and for space travel in low-Earth orbit and on interplanetary missions. Nuclear reactions, such as knockout and fragmentation, present large modifications of particle type and energies of the galactic cosmic rays in penetrating materials. We make an assessment of the current nuclear reaction models and improvements in these model for developing required transport code data bases. A new fragmentation data base (QMSFRG) based on microscopic models is compared to the NUCFRG2 model and implications for shield assessment made using the HZETRN radiation transport code. For deep penetration problems, the build-up of light particles, such as nucleons, light clusters and mesons from nuclear reactions in conjunction with the absorption of the heavy ions, leads to the dominance of the charge Z = 0, 1, and 2 hadrons in the exposures at large penetration depths. Light particles are produced through nuclear or cluster knockout and in evaporation events with characteristically distinct spectra which play unique roles in the build-up of secondary radiation's in shielding. We describe models of light particle production in nucleon and heavy ion induced reactions and make an assessment of the importance of light particle multiplicity and spectral parameters in these exposures.     

Investigation of dose and flux dynamics in the Liulin-5 dosimeter of the tissue-equivalent phantom onboard the Russian segment of the International Space Station.

Described is the Liulin-5 active dosimetric telescope designed for measurement of the space radiation dose depth-distribution in a human phantom on the Russian Segment of the International Space Station (ISS). The Liulin-5 experiment is a part of the international project MATROSHKA-R on ISS. The MATROSHKA-R project is aimed to study the depth-dose distribution at the sites of critical organs of the human body, using models of human body-anthropomorphic and spherical tissue-equivalent phantoms. The aim of Liulin-5 experiment is a long term (4-5 years) investigation of the radiation environment dynamics inside the spherical tissue-equivalent phantom, mounted in different compartments. Energy deposition spectra, linear energy transfer spectra, and flux and dose rates for charged particles will be measured simultaneously with near real time resolution at different depths of the phantom by means of three silicon detectors. Data obtained together with data from other active and passive dosimeters will be used to estimate the radiation risk to the crewmembers, which verify the models of radiation environment in low Earth orbit. Presented are the test results of the prototype unit. Liulin-5 will be flown on the ISS in the year 2003.     

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.     

Radiation measured during ISS-Expedition 13 with different dosimeters

Radiation in low Earth orbit (LEO) is mainly composed of galactic cosmic rays (GCR), solar energetic particles and particles in SAA (South Atlantic Anomaly). The biological impact of space radiation to astronauts depends strongly on the particles’ linear energy transfer (LET) and is dominated by high LET radiation. It is important to measure the LET spectrum for the space radiation field and to investigate the influence of radiation on astronauts. At present, the preferred active dosimeters sensitive to all LET are the tissue equivalent proportional counter (TEPC) and the silicon detectors in various configurations; the preferred passive dosimeters are CR-39 plastic nuclear track detectors (PNTDs) sensitive to high LET and thermoluminescence dosimeters (TLDs) as well as optically stimulated luminescence dosimeters (OSLDs) sensitive to low LET. The TEPC, CR-39 PNTDs, TLDs and OSLDs were used to investigate the radiation field for the ISS mission Expedition 13 (ISS-12S) in LEO. LET spectra and radiation quantities (fluence, absorbed dose, dose equivalent and quality factor) were measured for the space mission with different dosimeters. This paper introduces the role of high LET radiation in radiobiology, the operational principles for the different dosimeters, the LET spectrum method using CR-39 detectors, the method to combine the results measured with TLDs/OSLDs and CR-39 PNTDs, and presents the LET spectra and the radiation quantities measured and combined.     

Ground-based measurements of galactic cosmic ray fragmentation in shielding.

The mean free path for nuclear interactions of galactic cosmic-rays is comparable to shielding and tissue thicknesses present in human interplanetary exploration, resulting in a significant fraction of nuclear reaction products at depth. In order to characterize the radiation field, the energy spectrum, the angular distribution, and the multiplicity of each type of secondary particles must also be known as a function of depth. Reactions can take place anywhere in a thick absorber; therefore, it is necessary to know these quantities as a function of particle energy for all particles produced. HZE transport methods are used to predict the radiation field; they are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The most cost-effective method of validation is a comparison with ground-based experimental measurements. A research program to provide such validation measurements using neon, iron and other accelerated heavy ion beams will be discussed and illustrated using results from ongoing experiments and their comparison with current transport codes. The extent to which physical measurements yield radiobiological predictions will be discussed.     

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.     

An analytical model for the prediction of a micro-dosimeter response function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (δ ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCRs) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 s/ion. The ionizing radiation environment at LEO is represented by O’Neill’s GCR model (2004), covering charged particles in the 1 ? Z ? 28 range. O’Neill’s free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.     

Accommodations for earth-viewing payloads on the international space station

The design of the International Space Station (ISS) includes payload locations that are external to the pressurized environment. These external or attached payload accommodation locations will allow direct access to the space environment at the ISS orbit and direct viewing of the earth and space. NASA sponsored payloads will have access to several different types of standard external locations; the S3 Truss Sites, the Columbus External Payload Facility (EPF), and the Japanese Experiment Module Exposed Facility (JEM-EF). As the ISS Program develops, it may also be possible to locate external payloads at the P3 Truss Sites or at non-standard locations similar to the handrail-attached payloads that were flown during the MIR Program. Earth-viewing payloads may also be located within the pressurized volume of the US Lab in the Window Observational Research Facility (WORF). Payload accommodations at each of the locations will be described, as well as transport to and retrieval from the site.     

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.     

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.     

Analysis of the space radiation doses obtained simultaneously at two different locations outside the ISS

Space weather and related ionizing radiation has been recognized as one of the main health concerns for the International Space Station (ISS) crew. The estimation of the radiation effect on humans outside the ISS requires at first order accurate knowledge of their accumulated absorbed dose rates, which depend on the global space radiation distribution, solar cycle and local variations generated by the 3D mass distribution surrounding the ISS. The R3DE (Radiation Risks Radiometer-Dosimeter for the EXPOSE-E platform) on the European Technological Exposure Facility (EuTEF) worked successfully outside of the European Columbus module between February 2008 and September 2009. A very similar instrument named R3DR for the EXPOSE-R platform worked outside the Russian Zvezda module of the ISS between March 2009 and August 2010. Both are Liulin-type detectors, Bulgarian-built miniature spectrometer-dosimeters. The acquired approximately 5 million deposited energy spectra from which the flux and absorbed dose rate were calculated with 10 s resolution behind less than 0.41 g cm⁻² shielding. This paper analyses the spectra collected in 2009 by the R3DE/R instruments and the long-term variations in the different radiation environments of Galactic Cosmic Rays (GCR), inner radiation belt trapped protons in the region of the South Atlantic Anomaly (SAA) and relativistic electrons from the Outer Radiation Belt (ORB). The R3DE instrument, heavily shielded by the surrounding structures, measured smaller primary fluxes and dose rates from energetic protons from the SAA and relativistic electrons from the ORB but higher values from GCRs because of the contribution from secondary particles. The main conclusion from this investigation is that the dose rates from different radiation sources around the International Space Station (ISS) have a large special and temporal dynamic range. The collected data can be interpreted as possible doses obtained by the cosmonauts and astronauts during Extra Vehicular Activities (EVA) because the R3DE/R instruments shielding is very similar to the Russian and American space suits average shielding (,  and ). Fast, active measurements are required to assess accurately the dose accumulated by astronauts during EVA.     

The SOS-LUX-LAC-FLUORO-Toxicity-test on the International Space Station (ISS).

In the 21st century, an increasing number of astronauts will visit the International Space Station (ISS) for prolonged times. Therefore it is of utmost importance to provide necessary basic knowledge concerning risks to their health and their ability to work on the station and during extravehicular activities (EVA) in free space. It is the aim of one experiment of the German project TRIPLE-LUX (to be flown on the ISS) to provide an estimation of health risk resulting from exposure of the astronauts to the radiation in space inside the station as well as during extravehicular activities on one hand, and of exposure of astronauts to unavoidable or as yet unknown ISS-environmental genotoxic substances on the other. The project will (i) provide increased knowledge of the biological action of space radiation and enzymatic repair of DNA damage, (ii) uncover cellular mechanisms of synergistic interaction of microgravity and space radiation and (iii) examine the space craft milieu with highly specific biosensors. For these investigations, the bacterial biosensor SOS-LUX-LAC-FLUORO-Toxicity-test will be used, combining the SOS-LUX-Test invented at DLR Germany (Patent) with the commercially available LAC-FLUORO-Test. The SOS-LUX-Test comprises genetically modified bacteria transformed with the pBR322-derived plasmid pPLS-1. This plasmid carries the promoterless lux operon of Photobacterium leiognathi as a reporter element under control of the DNA-damage dependent SOS promoter of ColD as sensor element. This system reacts to radiation and other agents that induce DNA damages with a dose dependent measurable emission of bioluminescence of the transformed bacteria. The analogous LAC-FLUORO-Test has been developed for the detection of cellular responses to cytotoxins. It is based on the constitutive expression of green fluorescent protein (GFP) mediated by the bacterial protein expression vector pGFPuv (Clontech, Palo Alto, USA). In response to cytotoxic agents, this system reacts with a dose-dependent reduction of GFP-fluorescence. Currently, a fully automated miniaturized hardware system for the bacterial set up, which includes measurements of luminescence and fluorescence or absorption and the image analysis based evaluation is under development. During the first mission of the SOS-LUX-LAC-FLUORO-Toxicity-Test on the ISS, a standardized, DNA-damaging radiation source still to be determined will be used as a genotoxic inducer. A panel of recombinant Salmonella typhimurium strains carrying either the SOS-LUX plasmid or the fluorescence-mediating lac-GFPuv plasmid will be used to determine in parallel on one microplate the genotoxic and the cytotoxic action of the applied radiation in combination with microgravity. Either in addition to or in place of the fluorometric measurements of the cytotoxic agents, photometric measurements will simultaneously monitor cell growth, giving additional data on survival of the cells. The obtained data will be available on line during the TRIPLE-LUX mission time. Though it is the main goal during the TRIPLE-LUX mission to measure the radiation effect in microgravity, the SOS-LUX-LAC-FLUORO-Toxicity-test in principle is also applicable as a biomonitor for the detection and measurement of genotoxic substances in air or in the (recycled) water system on the ISS or on earth in general.     

Analysis of laser ranges and angular measurements data fusion for space debris orbit determination

In the framework of space debris, the orbit determination process is a fundamental step, both, for researchers and for satellite operators. The accurate knowledge of the orbit of space debris objects is needed to allow space debris characterization studies and to avoid unnecessary collision avoidance maneuvers.The accuracy of the results of an orbit determination process depends on several factors as the number, the accuracy, the kind of processed measurements, their distribution along the orbit, and the object-observer relative geometry. When the observation coverage of the target orbit is not homogeneous, the accuracy of the orbit determination can be improved processing different kind of observables. Recent studies showed that the satellite laser ranging technique can be successfully applied to space debris.In this paper, we will investigate the benefits of using laser ranges and angular measurements for the orbit determination process. We will analyze the influence of the number of used observations, of the covered arc of orbit, of each observable, and of the observation geometry on the estimated parameters. Finally, using data acquired on short observation arcs, we analyze the achievable accuracies for the orbital regimes with the highest space debris density, and to the consequences of the data fusion on catalog maintenance operations. The results shown are obtained using only real data (both angular and laser measurements) provided by sensors of the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald owned by the Astronomical Institute of the University of Bern (AIUB) and for some studies also using ranges provided from other stations of the International Laser Ranging Service (ILRS).