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.     

Feasibility of developing an ionospheric E-region electron density storm model using TIMED/SABER measurements

We present a new technique for improving ionospheric models of nighttime E-region electron densities under geomagnetic storm conditions using TIMED/SABER measurements of broadband 4.3 μm limb radiance. The response of E-region electron densities to geomagnetic activity is characterized by SABER-derived NO⁺(v) 4.3 μm Volume Emission Rates (VER). A storm-time E-region electron density correction factor is defined as the ratio of storm-enhanced NO⁺(v) VER to a quiet-time climatological average NO⁺(v) VER, which will be fit to a geomagnetic activity index in a future work. The purpose of this paper is to demonstrate the feasibility of our technique in two ways. One, we compare storm-to-quiet ratios of SABER-derived NO⁺(v) VER with storm-to-quiet ratios of electron densities measured by Incoherent Scatter Radar. Two, we demonstrate that NO⁺(v) VER can be parameterized by widely available geomagnetic activity indices. The storm-time correction derived from NO⁺(v) VER is applicable at high-latitudes.     

Star visibility and tracking from the Space Station

Purpose of the present study is to provide algorithms for and examples of how to simulate star visibility and tracking by a Telescope attached to the main truss of the International Space Station (ISS).     

Evolution of coronal mass ejection/shock system in interplanetary space

The evolution of coronal mass ejection/shock system is investigated by numerically solving the usual set of two-dimensional single-fluid polytropic magnetohydrodynamic equations from 1 R_s to 1 AU in the meridian plane. The simulation result reveals that the coronal mass ejection/shock system formed near the sun evolves into the magnetic cloud/shock system near the earth’s orbit through the following three phases: the initial formation, the dominant latitudinal expansion and the similar expansion.     

Space weather and commercial airlines

Space weather phenomena can effect many areas of commercial airline operations including avionics, communications and GPS navigation systems. Of particular importance at present is the recently introduced EU legislation requiring the monitoring of aircrew radiation exposure, including any variations at aircraft altitudes due to solar activity. With the introduction of new ultra-long-haul “over-the-pole” routes, “more-electric” aircraft in the future, and the increasing use of satellites in the operation, the need for a better understanding of the space weather impacts on future airline operations becomes all the more compelling. This paper will present the various space weather effects, some provisional results of an ongoing 3-year study to monitor cosmic radiation in aircraft, and conclude by summarising some of the identified key operational issues, which must be addressed, with the help of the science community, if the airlines want to benefit from the availability of space weather services.     

CME-geometry and cosmic-ray anisotropy observed by a prototype muon detector network

We analyze the cosmic-ray anisotropy observed by a prototype network of muon detectors during geomagnetic storms associated with coronal mass ejections (CMEs). The network currently consists of multidirectional surface muon detectors at Nagoya (Japan) and Hobart (Australia), together with a prototype detector at São Martinho (Brazil) which has been in operation since March, 2001. In this report, we analyze the anisotropy recorded in both the muon detector and neutron monitor (the Spaceship Earth) networks and find significant enhancements of cosmic-ray anisotropy during geomagnetic storms. Following the analysis by Bieber and Evenson [Bieber, J.W., Evenson, P. CME geometry in relation to cosmic ray anisotropy. Geophys. Res. Lett. 25 (1998) 2955–2958] for the neutron monitor data at 10 GeV, we also derive cosmic-ray density gradients from muon data at higher-energy (50 GeV), possibly reflecting the larger-scale geometry of CMEs causing geomagnetic storms. We particularly find in some events the anisotropy enhancement clearly starting prior to the storm onset in both the muon and neutron data. This is the first result of the CME-geometry derived from simultaneous observations of the anisotropy with networks of multidirectional muon detectors and neutron monitors.     

Sporadic emission of ultra high energy gamma rays from crab pulsar

In a further study of sporadic emission from pulsars we find evidence for short lived intense emission from the Crab pulsar.     

Results from the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Observatory

The Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory covers the high energy gamma ray energy range, approximately 30 MeV to 30 GeV, with a sensitivity considerably greater than earlier high energy gamma-ray satellites. Thus far, 4 pulsars have been detected and their properties measured, including in 3 cases the energy spectrum as a function of phase. The details of the galactic plane are being mapped and a spectra of the center region has been obtained in good agreement with that expected from cosmic ray interactions. The Magellanic clouds have been examined with the Large Magellanic Cloud having been detected at a level consistent with it having a cosmic ray density compatible with quasi-stable equilibrium. Sixteen Active Galactic Nuclei (AGN's) have been seen thus far with a high degree of certainty including 12 quasars and 4 BL Lac objects, but no Seyferts. Time variation has been detected in some of these AGN's.     

Cometary Deuterium

Deuterium fractionations in cometary ices provide important clues to the origin and evolution of comets. Mass spectrometers aboard spaceprobe Giotto revealed the first accurate D/H ratios in the water of Comet 1P/Halley. Ground-based observations of HDO in Comets C/1996 B2 (Hyakutake) and C/1995 O1 (Hale-Bopp), the detection of DCN in Comet Hale-Bopp, and upper limits for several other D-bearing molecules complement our limited sample of D/H measurements. On the basis of this data set all Oort cloud comets seem to exhibit a similar ratio in H₂O, enriched by about a factor of two relative to terrestrial water and approximately one order of magnitude relative to the protosolar value. Oort cloud comets, and by inference also classical short-period comets derived from the Kuiper Belt cannot be the only source for the Earth's oceans. The cometary O/C ratio and dynamical reasons make it difficult to defend an early influx of icy planetesimals from the Jupiter zone to the early Earth. D/H measurements of OH groups in phyllosilicate rich meteorites suggest a mixture of cometary water and water adsorbed from the nebula by the rocky grains that formed the bulk of the Earth may be responsible for the terrestrial D/H. The D/H ratio in cometary HCN is 7 times higher than the value in cometary H₂O. Species-dependent D-fractionations occur at low temperatures and low gas densities via ion-molecule or grain-surface reactions and cannot be explained by a pure solar nebula chemistry. It is plausible that cometary volatiles preserved the interstellar D fractionation. The observed D abundances set a lower limit to the formation temperature of (30 ± 10) K. Similar numbers can be derived from the ortho-to-para ratio in cometary water, from the absence of neon in cometary ices and the presence of S₂. Noble gases on Earth and Mars, and the relative abundance of cometary hydrocarbons place the comet formation temperature near 50 K. So far all cometary D/H measurements refer to bulk compositions, and it is conceivable that significant departures from the mean value could occur at the grain-size level. Strong isotope effects as a result of coma chemistry can be excluded for molecules H₂O and HCN. A comparison of the cometary ratio with values found in the atmospheres of the outer planets is consistent with the long-held idea that the gas planets formed around icy cores with a high cometary D/H ratio and subsequently accumulated significant amounts of H₂ from the solar nebula with a low protosolar D/H. This revised version was published online in June 2006 with corrections to the Cover Date.     

Semi-empirical MHD Model of the Solar Wind and its Comparison with Ulysses

We have developed a 2D semi-empirical model (Sittler and Guhathakurta 1999) of the corona and the interplanetary medium using the time independent MHD equations and assuming azimuthal symmetry, utilizing the SOHO, Spartan and Ulysses observations. The model uses as inputs (1) an empirically derived global electron density distribution using LASCO, Mark III and Spartan white light observations and in situ observations of the Ulysses spacecraft, and (2) an empirical model of the coronal magnetic field topology using SOHO/LASCO and EIT observations. The model requires an estimate of solar wind velocity as a function of latitude at 1 AU and the radial component of the magnetic field at 1 AU, for which we use Ulysses plasma and magnetic field data results respectively. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, temperature T_eff, and heat flux Q_eff which are derived from the equations of conservation of mass, momentum and energy, respectively, in the rotating frame of the Sun. The term "effective" indicates possible wave contributions. The model can be used as a planning tool for such missions as Solar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind. This revised version was published online in June 2006 with corrections to the Cover Date.