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.     

航天员受银河宇宙线辐射的剂量计算

在近地空间(LEO)和深空探测中,航天员遭受的辐射风险主要来自于银河宇宙线(GCR)照射.银河宇宙线的辐射剂量是航天员辐射风险评价的基础.国际放射防护委员会(ICRP)于2013年提出了新的航天员空间辐射剂量估算方法,以更准确给出空间重离子辐射的剂量.基于此方法,开发了宇宙线粒子在物质中输运的蒙特卡罗程序,并在程序中实现用中国成年男性人体数字模型来仿真航天员.采用该程序计算了粒子(Z=1~92)各向同性照射航天员时器官的通量-器官剂量转换因数,并估算出航天员在近地轨道空间受银河宇宙线辐射的剂量.     

Modeling and experimental study of the 27-day variation of galactic cosmic-ray intensity for a solar wind velocity depending on heliolongitude

We develop a three-dimensional (3-D) model of the 27-day variation of galactic cosmic-ray (GCR) intensity with a spatial variation of the solar wind velocity. A consistent, divergence-free interplanetary magnetic field is derived by solving the corresponding Maxwell equations with a variable solar wind speed, which reproduces in situ observed experimental data for the time interval to be analyzed (24 August 2007–28 February 2008). We perform model calculations for the GCR intensity using the variable solar wind and the corresponding magnetic field. Results are compatible with experimental data; the correlation coefficient between our model predictions and observed 27-day GCR variation is 0.80 ± 0.05.     

Nuclear fragmentation database for GCR transport code development

A critical need for NASA is the ability to accurately model the transport of heavy ions in the Galactic Cosmic Rays (GCR) through matter, including spacecraft walls, equipment racks, etc. Nuclear interactions are of great importance in the GCR transport problem, as they can cause fragmentation of the incoming ion into lighter ions. Since the radiation dose delivered by a particle is proportional to the square of (charge/velocity), fragmentation reduces the dose delivered by incident ions. The other mechanism by which dose can be reduced is ionization energy loss, which can lead to some particles stopping in the shielding. This is the conventional notion of shielding, but it is not applicable to human spaceflight since the particles in the GCR tend to be too energetic to be stopped in the relatively thin shielding that is possible within payload mass constraints. Our group has measured a large number of fragmentation cross sections, intended to be used as input to, or for validation of, NASA’s radiation transport models. A database containing over 200 charge-changing cross sections and over 2000 fragment production cross sections has been compiled. In this report, we examine in detail the contrast between fragment measurements at large acceptance and small acceptance. We use output from the PHITS Monte Carlo code to test our assumptions using as an example ⁴⁰Ar data (and simulated data) at a beam energy of 650 MeV/nucleon. We also present preliminary analysis in which isotopic resolution was attained for beryllium fragments produced by beams of ¹⁰B and ¹¹B. Future work on the experimental data set will focus on extracting and interpreting production cross sections for light fragments.     

Galactic cosmic ray abundances and spectra behind defined shielding.

LET spectra have been measured for lunar missions and for several near Earth orbits ranging from 28 degrees to 83 degrees inclination. In some of the experiments the flux of GCR was determined separately from contributions caused by interactions in the detector material. Results of these experiments are compared to model calculations. The general agreement justifies the use of the model to calculate GCR fluxes. The magnitude of variations caused by solar modulation, geomagnetic shielding, and shielding by matter determined from calculated LET spectra is generally in agreement with experimental data. However, more detailed investigations show that there are some weak points in modeling solar modulation and shielding by material. These points are discussed in more detail.     

Flight safety implications of the extreme solar proton event of 23 February 1956

There is considerable speculation about the effects at aircraft altitudes resulting from extreme solar proton events. The ground level event (GLE) of 23 February 1956 (GLE 5), remains the largest solar proton event of the neutron monitor era in terms of its influence on count rates at monitors near sea level. During this GLE the count rate was increased by as much as 4760% (15-min average) at the Leeds monitor relative to the count rate from galactic cosmic radiation (GCR). Two modern models of the event cumulative solar proton spectrum for this event, a 6-parameter fit in energy and a 4-parameter Band fit in rigidity, are compared with 1-h of GCR at solar minimum. While effective doses calculated with CARI-7A for both models at low geomagnetic cutoff rigidities are indeed high when compared with GCR and can exceed recommended exposure limits, both GLE spectra exhibit a much stronger dependence on cutoff rigidity than GCR, and a larger fraction of the dose from neutrons. At locations with cutoff rigidities above 4.2 and 6.4?GV, respectively, the GLE effective doses are smaller than the GCR hourly dose. At locations with cutoff rigidities above about 4?GV, GCR was the dominant source of exposure in 10?h or less at all altitudes examined. This suggests that if a similar event occurs in the future, low- and mid-latitude flights at modern jet flight altitudes could be well-protected by Earth’s magnetic field.     

Sunspot numbers,interplanetary magnetic field,and cosmic ray intensity at earth: Nexus for the twentieth century

The pivotal role played by the interplanetary magnetic field (B) in modulating galactic cosmic ray (GCR) intensity in the heliosphere is described. We show that the inverse correlation observed by Forbush (1958) between GCRs and sunspot numbers (SSNs) is reflected in high correlation between SSNs and B (cc = 0.94). The SSN data are available since 1700 and the derived B data since 1835. The paleo-cosmic ray data are available for several millennia in the form of ¹⁰Be radionuclide sequestered in polar ice. The data of the ion chambers (ICs) at the Cheltenham–Fredericksburg–Yakutsk (CFY) sites are combined to create a data string for 1937–1988. In turn, these data are used to extend the measurements of the low energy GCR ions (>0.1 GeV) at balloon altitudes at high latitudes in Russia to 1937. These data are then correlated to B and the fit parameters are used to extend the low energy ion data to 1900, creating the instrumental era GCR time series for the twentieth century. The derived GCR time series is compared to ¹⁰Be measured at two sites in Greenland, namely Dye 3 and NGRIP for 1900–2000 to check the internal consistency of datasets for the long-term trend. We find that the annual mean rate (%) for 1965 at NGRIP is an outlier. We replace it with the mean of 1964 and 1965 rates and construct a new re-normalized time series at NGIP, improving the agreement with the derived instrumental era GCR time series for the twentieth century as well. This should encourage its use by heliophysics community for varied applications.     

A model for predicting the radiation exposure for mission planning aboard the international space station

The International Space Station Cosmic Radiation Exposure Model (ISSCREM) has been developed as a possible tool for use in radiation mission planning as based on operational data collected with a tissue equivalent proportional counter (TEPC) aboard the ISS since 2000. It is able to reproduce the observed trapped radiation and galactic cosmic radiation (GCR) contributions to the total dose equivalent to within ±20% and ±10%, respectively, as would be measured by the onboard TEPC at the Zvezda Service Module panel 327 (SM-327). Furthermore, when these contributions are combined, the total dose equivalent that would be measured at this location is estimated to within ±10%. The models incorporated into ISSCREM correlate the GCR dose equivalent rate to the cutoff rigidity magnetic shielding parameter and the trapped radiation dose equivalent rate to atmospheric density inside the South Atlantic Anomaly. The GCR dose equivalent rate is found to vary minimally with altitude and TEPC module location however, due to the statistics and data available, the trapped radiation model could only be developed for the TEPC located at SM-327. Evidence of the variation in trapped radiation dose with detector orientation and the East–West asymmetry were observed at this location.     

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.     

离子推力器羽流热效应仿真分析

离子推力器工作时向外喷出的羽流与航天器表面碰撞,会引起敏感材料热变形等热效应,严重时会导致航天任务失败。针对兰州空间技术物理研究所研制的LIPS-200型离子推力器羽流热效应进行了仿真分析。仿真中,使用粒子网格（PIC）方法处理等离子体运动,使用直接模拟蒙特卡罗（DSMC）方法处理粒子间碰撞,使用Maxwell模型处理粒子与壁面的能量交换,对电推进羽流热效应测量中的部分测点进行了数值模拟。结果表明,仿真结果与实验数据符合较好,离子推力器出口轴线上滞止热流仿真值与实验测量值误差小于17.0%。此外,热流计对流场的影响主要集中在热流计附近0.1 m范围内,对整体流场影响较小。      

Low fluence.

The question of the appropriate extrapolation to low dose has long been a subject of controversy. A linear no-threshold model is favored by regulatory bodies as the basis of RBE assignments and estimates of radiation hazards to the general population. This model is largely supported by extensive application of the linear-quadratic survival formula "fitted" statistically to a wide variety of experimental data obtained at doses typically exceeding 1 Gy, and then extrapolated to mGy for practical applications, and even to the prediction of hazards from single electrons. Such extrapolations are questionable at best, and may even prove hazardous for risk evaluations. Fluence and geometry rather than dose based data are proposed as a basis for a limiting "threshold" for a "low dose" extrapolation. The proposed threshold is one where the fluence of particles is one per square micron, where on average only 2/3 of the 1 micrometers2 pixels covering an irradiated area are traversed by one or more particles. The corresponding dose threshold is determined by the LET of the bombarding radiation. For relativistic electrons this dose is about 0.032 Gy.     

Neutron yields from interactions of GCR-like beams in stopping targets.

In order to accurately determine the radiation risk to astronauts from GCR, the nature of the secondary radiation field created by the fragmentation of GCR in shielding and tissue must be understood. Due to the their high penetrabilities, neutrons are an important component of the secondary radiation field, especially for astronauts protected by thick shielding on lunar or Martian bases. Neutron yields from 435A MeV and 272A MeV Nb stopping in Nb and Al targets are presented, along with some preliminary analysis of neutron yields from 155A MeV C stopping in Al. Energy spectra and angular distributions are shown for neutron energies above 20 MeV. The data provides some information about the dependence of the neutron yield on projectile energy and target mass. Comparisons of the data with BUU calculations are also shown.     

A model of galactic cosmic rays for use in calculating linear energy transfer spectra.

The galactic cosmic rays (GCR) contain fully stripped nuclei, from Hydrogen to beyond the Iron group, accelerated to high energies and are a major component of the background radiation encountered by satellites and interplanetary spacecraft. This paper presents a GCR model which is based upon our current understanding of the astrophysics of GCR transport through interstellar and interplanetary space. The model can be used to predict the energy spectra for all stable and long-lived radioactive species from H to Ni over an energy range from 50 to 50,000 MeV/nucleon as a function of a single parameter, the solar modulation level phi. The details of this model are summarized, phi is derived for the period 1974 to present, and results from this model during the 1990/1991 CRRES mission are presented.     

Solar modulation and nuclear fragmentation effects in galactic cosmic ray transport through shielding.

Crews of manned interplanetary missions may accumulate significant radiation exposures from the galactic cosmic ray (GCR) environment in space. Estimates of how these dose levels are affected by the assumed temporal and spatial variations in the composition of the GCR environment, and by the effects of the spacecraft and body self-shielding on the transported radiation fields are presented. In this work, the physical processes through which shielding alters the transported radiation fields are described. We then present estimates of the effects on model calculations of (1) nuclear fragmentation model uncertainties, (2) solar modulation, (3) variations between solar cycles, and (4) proposed changes to the quality factors which relate dose equivalent to absorbed dose.     

Fluence-related risk coefficients using the Harderian gland data as an example.

The risk of radiation-induced cancer to space travelers outside the earth's magnetosphere will be of concern on missions to the Moon and beyond to Mars. High energy galactic cosmic rays with high charge (HZE particles) will penetrate the spacecraft and the bodies of the astronauts, sometimes fragmenting into nuclear secondary species of lower charge but always ionizing densely, thus causing cellular damage which may lead to malignant transformation. To quantitate this risk, the concept of dose equivalent (in which a quality factor Q as a function of LET is assumed) may not be adequate, since different particles of the same LET may have different efficiencies for tumor induction. Also, RBE values on which quality factors are based depend on response to low-LET radiation at low doses, a very difficult region for which to obtain reliable experimental data. Thus, we introduce a new concept, a fluence-related risk coefficient (F), which is the risk of a cancer per unit particle fluence and which we call the risk cross section. The total risk is the sum of the risk from each particle type: sigma i integral Fi(Li) phi i(Li) dLi, where Li is the LET and phi i(Li) is the fluence-LET spectrum of the ith particle type. As an example, tumor prevalence data in mice are used to estimate the probability of mouse Harderian gland tumor induction per year on an extra-magnetospheric mission inside an idealized shielding configuration of a spherical aluminum shell 1 g/cm2 thick. The combined shielding code BRYNTRN/GCR is used to generate the LET spectra at the center of the sphere. Results indicate a yearly prevalence at solar minimum conditions of 0.06, with 60% of this arising from charge components with Z between 10 and 28, and two-thirds of the contribution arising from LET components between 10 and 200 keV/micrometers.     

Radiation issues for piloted Mars mission.

Man is now entering an era of colonizing the moon and exploration of Mars. The crewmembers of a piloted mission to Mars will be exposed to inner belt trapped protons, the outer trapped electrons, and the galactic cosmic radiation. In addition there is always the added risk of acute exposure to a solar particle event. Current radiation risk is estimated using the idea of absorbed dose and ICRP-26, LET-dependent quality factors. In a spacecraft with aluminum walls (2 g cm-2) at solar minimum the calculated dose equivalent is 0.73 Sv for a 406-day mission. Based on the current thinking this leads to an excess cancer mortality in a 35 year male of about 1%. About 75% of the dose equivalent is contributed by HZE particles and target fragments with average quality factors of 10.3 and 20, respectively. The entire concept of absorbed dose, quality factor, and dose equivalent as applied to such missions needs to be reexamined, in light of the fact that less than 50% of the nuclei in the body of the astronaut would have been traversed by a single GCR nuclei in the 406-day mission. Clearly, more biologically relevant information about the effects of heavy ions and target fragments is needed and fluence based risk estimation strategy developed for such long term stays in space.     

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.     

Recent results form measurements of the energy spectrum of cosmic-ray induced neutrons aboard an ER-2 airplane and on the ground.

Crews of future high-altitude commercial aircraft may be significantly exposed to atmospheric cosmic radiation from galactic cosmic rays (GCR). To help determine such exposures, the Atmospheric Ionizing Radiation Project, an international collaboration of 15 laboratories, made simultaneous radiation measurements with 14 instruments on a NASA ER-2 high-altitude aircraft. The primary instrument was a sensitive extended-energy multisphere neutron spectrometer, which was also used to make measurements on the ground. Its detector responses were calculated for neutrons and charged hadrons at energies up to 100 GeV using the radiation transport code MCNPX. We have now recalculated the detector responses including the effects of the airplane structure. We are also using new FLUKA calculations of GCR-induced hadron spectra in the atmosphere to correct for spectrometer counts produced by charged hadrons. Neutron spectra are unfolded from the corrected measured count rates using the MAXED code. Results for the measured cosmic-ray neutron spectrum (thermal to >10 GeV), total neutron fluence rate, and neutron dose equivalent and effective dose rates, and their dependence on altitude and geomagnetic cutoff generally agree well with results from recent calculations of GCR-induced neutron spectra.     

Interplanetary proton cumulated fluence model update

Solar particle events leading to important increase of particle fluxes at energies of order of magnitude ranging from MeV to GeV constitute an important hazard for space missions. They may lead to effects seen in microelectronics or damage to solar cells and constitute a potential hazard for manned missions. Cumulative damage is commonly expressed as a function of fluence which is defined as the integral of the flux over time. A priori deterministic estimates of the expected fluence cannot be made because over the time scale of a space mission, the fluence can be dominated by the contribution of a few rare and unpredictable high intensity events. Therefore, statistical approaches are required in order to estimate fluences likely to be encountered by a space mission in advance. This paper extends work done by Rosenqvist et al. [Rosenqvist, L., Hilgers, A., Evans, H., Daly, E., Hapgood, M., Stamper, R., Zwickl, R., Bourdarie, S., Boscher, D. Toolkit for updating interplanetary proton-cumulated fluence models. J. Spacecraft Rockets, 42(6), 1077–1090, 2005] to describe an updated predictive engineering model for the proton interplanetary fluence with energies >30 MeV. This model is derived from a complete list of solar proton fluences based on data from a number of calibrated sources covering almost three solar cycles.     

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.