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.     

Nuclear model calculations and their role in space radiation research.

Proper assessments of spacecraft shielding requirements and concomitant estimates of risk to spacecraft crews from energetic space radiation requires accurate, quantitative methods of characterizing the compositional changes in these radiation fields as they pass through thick absorbers. These quantitative methods are also needed for characterizing accelerator beams used in space radiobiology studies. Because of the impracticality/impossibility of measuring these altered radiation fields inside critical internal body organs of biological test specimens and humans, computational methods rather than direct measurements must be used. Since composition changes in the fields arise from nuclear interaction processes (elastic, inelastic and breakup), knowledge of the appropriate cross sections and spectra must be available. Experiments alone cannot provide the necessary cross section and secondary particle (neutron and charged particle) spectral data because of the large number of nuclear species and wide range of energies involved in space radiation research. Hence, nuclear models are needed. In this paper current methods of predicting total and absorption cross sections and secondary particle (neutrons and ions) yields and spectra for space radiation protection analyses are reviewed. Model shortcomings are discussed and future needs presented.     

Ground-based simulations of galactic cosmic ray fragmentation and transport.

Since mean free paths for nuclear fragmentation are of the order of the ranges of primary Galactic Cosmic Ray (GCR) nuclei, determination of the radiation field produced by successive fragmentations of nuclei in material and tissue is essential to accurate assessment of GCR radiation risk to humans on long-duration missions outside the geomagnetosphere. We describe some recent measurements made at the Bevalac of heavy ion transport through materials, with representative results and examples of how they may be applied to aspects of the space radiation problem, including efforts to devise analytical tools for predicting biological effects and for designing spacecraft shielding.     

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.     

Fragmentation studies of relativistic iron ions using plastic nuclear track detectors.

We measured fluence and fragmentation of high-energy (1 or 5 A GeV) 56Fe ions accelerated at the Alternating Gradient Synchrotron or at the NASA Space Radiation Laboratory (Brookhaven National Laboratory, NY, USA) using solid-state CR-39 nuclear track detectors. Different targets (polyethylene, PMMA, C, Al, Pb) were used to produce a large spectrum of charged fragments. CR-39 plastics were exposed both in front and behind the shielding block (thickness ranging from 5 to 30 g/cm2) at a normal incidence and low fluence. The radiation dose deposited by surviving Fe ions and charged fragments was measured behind the shield using an ionization chamber. The distribution of the measured track size was exploited to distinguish the primary 56Fe ions tracks from the lighter fragments. Measurements of projectile's fluence in front of the shield were used to determine the dose per incident particle behind the block. Simultaneous measurements of primary 56Fe ion tracks in front and behind the shield were used to evaluate the fraction of surviving iron projectiles and the total charge-changing fragmentation cross-section. These physical measurements will be used to characterize the beam used in parallel biological experiments.     

Computational methods for the HZETRN code.

Asymptotic expansion has been used to simplify the transport of high charge and energy ions for broad beam applications in the laboratory and space. The solution of the lowest order asymptotic term is then related to a Green's function for energy loss and straggling coupled to nuclear attenuation providing the lowest order term in a rapidly converging Neumann series for which higher order collisions terms are related to the fragmentation events including energy dispersion and downshift. The first and second Neumann corrections were evaluated numerically as a standard for further analytic approximation. The first Neumann correction is accurately evaluated over the saddle point whose width is determined by the energy dispersion and located at the downshifted ion collision energy. Introduction of the first Neumann correction leads to significant simplification of the second correction term allowing application of the mean value theorem and a second saddle point approximation. The regular dependence of the second correction spectral dependence lends hope to simple approximation to higher corrections. At sufficiently high energy nuclear cross-section variations are small allowing non-perturbative methods to all orders and renormalization of the second corrections allow accurate evaluation of the full Neumann series.     

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.     

Removing Solar Radiation Based on the Empirical Mode Decomposition Method for Seismo-Ionospheric Anomaly Before the M9.0 Tohoku Earthquake

The ionospheric effect from solar activity can be seen as the background in the process of detecting the ionospheric precursor prior to strong earthquakes. The ionospheric variation induced by the forthcoming earthquake can be covered by the strong solar background during the period of high solar activity. The issue of how to remove the ionospheric effect from solar radiation is of outstanding significance. In this paper, a method of Empirical Mode Decomposition (EMD) is used to eliminate the solar background. As a case study, the global ionospheric map TEC before the M9.0 Tohoku earthquake on 11 March 2011 is analyzed. After the effect of solar radiation is removed using the EMD method, the precursor of the imminent earthquake is more obvious. The ionospheric anomaly had a local character and only appeared close to the earthquake epicenter while the useful signals were covered by the solar radiation background with traditional method, which implies that the EMD method is effective in eliminating solar radiation background.      

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.     

Solar particle event organ doses and dose equivalents for interplanetary crews: variations due to body size.

Proper assessments of spacecraft shielding requirements and concomitant estimates of risk to critical body organs of spacecraft crews from energetic space radiation require accurate, quantitative methods of characterizing the compositional changes in these radiation fields as they pass through the spacecraft and overlying tissue. When estimating astronaut radiation organ doses and dose equivalents it is customary to use the Computerized Anatomical Man (CAM) model of human geometry to account for body self-shielding. Usually, the distribution for the 50th percentile man (175 cm height; 70 kg mass) is used. Most male members of the U.S. astronaut corps are taller and nearly all have heights that deviate from the 175 cm mean. In this work, estimates of critical organ doses and dose equivalents for interplanetary crews exposed to an event similar to the October 1989 solar particle event are presented for male body sizes that vary from the 5th to the 95th percentiles. Overall the results suggest that calculations of organ dose and dose equivalent may vary by as much as approximately 15% as body size is varied from the 5th to the 95th percentile in the population used to derive the CAM model data.     

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.     

海洋二号卫星微波散射计自然扩展目标在轨定标

海洋二号卫星微波散射计(HY-2 Scatterometer, HY-2 SCAT)是一种旋转扫描笔形波束散射计, 能够对同一观测面元提供4次方位角和入射角的观测组合, 并通过地球物理模型反演海面风场. 为达到设计的风场反演精度, 要求其系统定标精度达到0.5dB. 利用不同区域自然扩展目标对HY-2 SCAT进行在轨外定标, 并与OSCAT散射计在同时期内的测量结果进行了对比. 定标结果显示可以消除因散射计天线指向偏差带来的方位向测量误差. 针对HY-2 SCAT方位向测量偏差进行了误差分析, 利用仿真方法以及海洋二号卫星雷达高度计同期测量数据的反演结果进行比对, 验证了误差来源.      

基于数据网格化方法的低轨辐射带建模技术

地球辐射带中的高能带电粒子是引起航天器材料和器件性能退化甚至失效的主要空间环境因素.因此,航天器设计中所采用的辐射带模型的准确程度对于航天器的生存能力和航天任务的完成质量至关重要.在利用我国自主辐射带高能粒子探测数据进行的辐射带建模中,探测数据的空间网格化是一项非常重要的工作.介绍了我国辐射带探测数据的情况,以及辐射带建模的方法和步骤;重点研究了不同插值方法在低地球轨道(LEO)空间辐射带建模数据网格化中的应用,并开展了误差分析.研究结果表明:在各种常用的插值方法中,反距离加权法、自然邻点法和最近邻点法适合工程化应用.其中,反距离加权法生成的数据网格对粒子通量的反演结果精度最高,该方法采用低阶距离时得到的反演结果更为合理.      

Space radiation dosimetry using bubble detectors.

Bubble detectors--a new development in radiation detection--has only recently been used for radiation measurements in space. One important characteristic of the bubble detector is that it operates on a phenomenon which bears considerable resemblance to biological response. Recent experimental results from irradiating bubble detectors with high-energy heavy ions point to the need to re-examine the methodology used for assessing space radiation and the relevance of conventional quantities such as dose equivalent for space dosimetry. It may be that biological hazard associated with the intensely ionizing events--associated with nuclear fragmentation but delivering relatively small dose equivalent--may be much more important than that associated with lightly ionizing events which comprise the bulk of the conventional radiation dose equivalent.     

Theoretical model of HZE particle fragmentation by hydrogen targets.

The fragmenting of high energy, heavy ions (HZE particles) by hydrogen targets is an important, physical process in several areas of space radiation research. In this work quantum mechanical optical model methods for estimating cross sections for HZE particle fragmentation by hydrogen targets are presented. The cross sections are calculated using a modified abrasion-ablation collision formalism adapted from a nucleus-nucleus collision model. Elemental and isotopic production cross sections are estimated and compared with report measurements for the breakup of neon, sulphur, and iron, nuclei at incident energies between 400 and 910 MeV/nucleon. Good agreement between theory and experiment is obtained.     

Inner radiation belt source of helium and heavy hydrogen isotopes

Nuclear interactions between inner zone protons and atoms in the upper atmosphere provide the main source of energetic H and He isotopes nuclei in the radiation belt. This paper reports on the specified calculations of these isotope intensities using various inner zone proton intensity models (AP-8 and SAMPEX/PET PSB97), the atmosphere drift-averaged composition and density model MSIS-90, and cross-sections of the interaction processes from the GNASH nuclear model code. To calculate drift-averaged densities and energy losses of secondaries, the particles were tracked in the geomagnetic field (modelled through IGRF-95) by integrating numerically the equation of the motion. The calculations take into account the kinematics of nuclear interactions along the whole trajectory of trapped proton. The comparison with new data obtained from the experiments on board RESURS-04 and MITA satellites and with data from SAMPEX and CRRES satellites taken during different phases of solar activity shows that the upper atmosphere is a sufficient source for inner zone helium and heavy hydrogen isotopes. The calculation results are energy spectra and angular distributions of light nuclear isotopes in the inner radiation belt that may be used to develop helium inner radiation belt model and to evaluate their contribution to SEU (single event upset) rates.     

A comparison of total reaction cross section models used in particle and heavy ion transport codes

Understanding the interactions and propagations of high energy protons and heavy ions are essential when trying to estimate the biological effects of Galactic Cosmic Rays (GCR) and Solar Particle Events (SPE) on personnel in space. To be able to calculate the shielding properties of different materials and radiation risks, particle and heavy ion transport codes are needed. In all particle and heavy ion transport codes, the probability function that a projectile particle will collide within a certain distance x in the matter depends on the total reaction cross sections, and the calculated partial fragmentation cross sections scale with the total reaction cross sections. It is therefore crucial that accurate total reaction cross section models are used in the transport calculations. In this paper, different models for calculating nucleon–nucleus and nucleus–nucleus total reaction cross sections are compared with each other and with measurements. The uncertainties in the calculations with the different models are discussed, as well as their overall performances with respect to the available experimental data. Finally, a new compilation of experimental data is briefly presented.