Challenges to cosmic ray modeling: From beyond the solar wind termination shock

Voyager 1 crossed the solar wind termination shock on December 16, 2004 at a distance of 94 AU from the Sun, to become the first spacecraft to explore the termination shock region and to enter the heliosheath, the final heliospheric frontier. By the end of 2006, Voyager 1 will be at ∼101 AU, with Voyager 2 at ∼81 AU and still approaching the termination shock. Both spacecraft have been observing the modulation of galactic and anomalous cosmic rays since their launch in 1977. The recent observations close to or inside the heliosheath have provided several interesting ‘surprises’ with subsequent theoretical and modeling challenges. Examples are: what does the modulation of galactic cosmic rays amount to in this region?; how do the anomalous cosmic rays get accelerated and modulated?; why are there ‘breaks’ in the power-law slopes of the spectra of accelerated particles? Several numerical models have been applied to most of these topics over the years and comprehensive global predictions have been made the past decade, thought to be based on reasonable assumptions about the termination shock and the heliosheath. Examples of these predictions and assumptions are concisely discussed within the context of the main observed features of cosmic rays in the vicinity of the termination shock, ending with a discussion of some of the issues and challenges to cosmic ray modeling in particular.     

Cosmic ray anisotropies in the outer heliosphere

The observation of the directional distribution of energetic and cosmic ray particles has been done with the Voyager spacecraft over a long period. Since 2002, when the first flux enhancements of charged particles associated with the approach of Voyager 1 to the solar wind termination shock were observed, these anisotropy measurements have become of special interest. They play an important role to understand the magnetic field and shock structure and the basics of the modulation of cosmic ray and anomalous particles at and beyond the termination shock. They also serve as motivation to study the spatial behavior of galactic and anomalous cosmic ray anisotropies with numerical modulation models in order to illustrate how the radial anisotropy, at different energies, change from upstream to downstream of the termination shock. Observations made by Voyager 1 indicate that the termination shock is a complicated region than previously thought, hence the effects of the latitude dependence of the termination shock’s compression ratio and injection efficiency on the radial anisotropies of galactic and anomalous protons will be illustrated. We find that the magnitude and direction of the radial anisotropy strongly depends on the position in the heliosphere and the energy of particles. The effect of the TS on the radial anisotropy is to abruptly increase its value in the heliosheath especially in the A > 0 cycle for galactic protons and in both polarity cycles for anomalous protons. Furthermore, the global effect of the latitude dependence of the shock’s compression ratio is to increase the radial anisotropy for galactic protons throughout the heliosphere, while when combined with the latitude dependence of the injection efficiency this increase depends on modulation factors for anomalous protons and can even alter the direction of the radial anisotropy.     

Radial distribution of galactic cosmic rays at solar maximum

In this paper we analyze the spatial distribution of galactic cosmic rays during periods of maximum solar activity of the cycles 21, 22 and 23. We have used a two dimensional model to solve the cosmic ray transport equation. This model includes all relevant physical processes: diffusion, convection, drift and shock effects on cosmic ray propagation inside the heliosphere. We focused on the study of the radial distribution of galactic cosmic rays, and compare our results with the spacecraft observations for two energies (175 MeV H and 265 MeV/n He). Although the radial intensities of galactic cosmic rays can be explained qualitatively with all three local interstellar spectra (LISs) used in this work, we applied a reduced chi-squared analysis to investigate the best LIS that could fit the data.     

Calculated cosmic rays dynamics before the geomagnetic storms onset caused by the interplanetary shocks

Dynamics of the galactic cosmic ray intensity caused by their interactions with a shock front in the September 9, 1992 event has been determined. Corresponding variations of the cosmic ray intensity have been calculated for different stations of the world network of neutron monitors and muon telescopes of stations Nagoya and Sakashita. Comparison calculated results with observational data shows, in general, satisfactory consensus both on amplitude and in time. The developed method can be used for investigation dynamics of the solar wind disturbances precursors in the cosmic rays.     

On statistical relationship of solar, geomagnetic and human activities.

Data of galactic cosmic rays, solar and geomagnetic activities and solar wind parameters on the one side and car accident events (CAE) in Poland on the other have been analyzed in order to reveal the statistical relationships among them for the period of 1990-2001. Cross correlation and cross spectrum of the galactic cosmic ray intensity, the solar wind (SW) velocity, Kp index of geomagnetic activity and CAE in Poland have been carried out. It is shown that in some epochs of the above-mentioned period there is found a reliable relationship between CAE and solar and geomagnetic activities parameters in the range of the different periodicities, especially, 7 days. The periodicity of 7 days revealed in the data of the CAE has the maximum on Friday without any exception for the minimum and maximum epochs of solar activity. However, the periodicity of 7 days is reliably revealed in other parameters characterizing galactic cosmic rays, SW, solar and geomagnetic activities, especially for the minimum epoch of solar activity. The periodicity of 3.5 days found in the series of CAE data more or less can be completely ascribed to the social effects, while the periodicity of 7 days can be ascribed to the social effect or/to the processes on the Sun, in the interplanetary space and in the Earth's magnetosphere and atmosphere.     

Acceleration of interplanetary pick-up ions and anomalous cosmic rays

It may not be doubted anymore that anomalous cosmic rays (ACRs) are produced in the heliosphere from interplanetary pick-up ions through their acceleration at the solar wind termination shock. However, there is no general agreement in the community of heliospheric researchers concerning the mechanism of injection of the pick-up ions into the shock acceleration. We discuss here three possible ways for pick-up ions to be involved into the acceleration process at the termination shock: (1) preacceleration of pick-up ions in the whole region from the Sun up to the termination shock by solar wind turbulences and interplanetary shock waves, (2) local preacceleration of pick-up ions in a vicinity of the termination shock by shock surfing, and (3) formation of high-velocity tails in pick-up ion spectra consisting of secondary pick-up ions which are produced in the supersonic solar wind due to ionization of energetic neutral atoms entering from the inner heliosheath.     

The physics of cosmic-ray modulation

The theory of the modulation of galactic cosmic rays by the solar wind is reviewed. The basic transport equation is presented, interpreted and then applied to cosmic-ray transport in a model heliosphere immersed in a constant uniform bath of galactic cosmic rays. The results of numerical modelling are presented and the dominant physical effects analyzed. A variety of observational tests of the model which were reported over the last several years are summarized and shown, generally, to support a model in which particle drifts play an important role. Recent measurements which show that the latitudinal gradient of cosmic rays changed sign in the recent sunspot minimum (relative the last sunspot minimum) are shown to provide additional, strong, support for the model. A new picture of the interplanetary magnetic field is presented, which gives promise of improving considerably the agreement between the theory and observations in the few remaining problem areas.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.     

Cosmic ray modulation with a Fisk-type heliospheric magnetic field and a latitude-dependent solar wind speed

The effect of a latitude-dependent solar wind speed on a Fisk heliospheric magnetic field [Fisk, L. A. Motion of the footpoints of heliospheric magnetic field lines at the Sun: implications for recurrent energetic particle events at high heliographic latitudes. J. Geophys. Res. 101, 15547–15553, 1996] was first discussed by Schwadron and Schwadron and McComas [Schwadron, N.A. An explanation for strongly underwound magnetic field in co-rotating rarefaction regions and its relationship to footpoint motion on the the sun. Geophys. Res. Lett. 29, 1–8, 2002. and Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003]. Burger and Sello [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005] found a significant effect for a simplified 2D version of a latitude-dependent Fisk-type field while Miyake and Yanagita [Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, 445–448, 2007] found a smaller effect. The current report improves on a previous attempt Hitge and Burger [Hitge, M., Burger, R.A. The effect of a latitude-dependent solar wind speed on cosmic-ray modulation in a Fisk-type heliospheric magnetic field. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 449–450, 2007] where the global change in the solar wind speed and not the local speed gradient was emphasized. The sheared Fisk field of Schwadron and McComas [Schwadron, N.A., McComas, D.J. Heliospheric “FALTS”: Favored acceleration locations at the termination shock. Geophys. Res. Lett. 30, 41–1, 2003.) is similar to the current Schwadron–Parker hybrid field. Little difference is found between the effects of a Parker field and a Schwadron–Parker hybrid field on cosmic-ray modulation, in contrast to the results of Burger and Sello and Miyake and Yanagita [Burger, R.A., Sello, P.C. The effect on cosmic ray modulation of a Parker field modified by a latitudinal-dependent solar wind speed. Adv. Space Res. 35, 643–646, 2005 and Miyake, S., Yanagita, S. The effect of a modified Parker field on the modulation of the galactic cosmic rays. In: Proceedings of 30th International Cosmic Ray Conference. Merida, Mexico, vol. 1, pp. 445–448, 2007]. The two-dimensional approximation used by these authors is therefore inadequate to model the complexities of the actual three-dimensional field. We also show that a Fisk-type field with a latitude-dependent solar wind speed (Schwadron–Parker hybrid field) decreases both the relative amplitude of recurrent cosmic ray intensity variations and latitude gradients and yields similar constants of proportionality for these quantities as for the constant solar wind speed case.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.     

Modeling the acceleration and modulation of anomalous cosmic ray oxygen

After the solar wind termination shock crossings of the Voyager spacecraft, the acceleration of anomalous cosmic rays has become a very contentious subject. In this paper we examine several topics pertinent to anomalous cosmic ray oxygen acceleration and transport using a numerical cosmic ray modulation model. These include the effects of drifts on a purely Fermi I accelerated spectra, the effects of introducing higher charge states of oxygen into the modulation model, examining the viability of momentum diffusion as a re-acceleration process in the heliosheath and examining energy spectra, and intensity gradients, in the inner heliosphere during consecutive drift cycles.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

Self-similar solutions to stellar winds with cosmic rays

A two-fluid model is used to study the time evolution of stellar winds including the dynamical effect of cosmic rays. Neglecting the diffusion of cosmic rays, we seek self-similar solutions to spherically symmetric winds with a termination shock. The velocity upstream of the shock is taken to be zero. Physical solutions are those that can connect the shock to the star with the velocity approaching zero at the star. Two parameters govern the behaviour of the solutions, namely, the ratio of the upstream sound speed to the shock speed (in an inertial frame) and the gravitation potential of the central star. In some parameter regimes, no physical solution is possible.     

The composition of cosmic rays at high energies

Measurements of the composition of the cosmic rays at high energies, and of the energy spectra of the individual components provide the basis for the understanding of the sources, of the acceleration mechanism, and of the galactic containment of these particles. We briefly review the presently available information, and we describe a recent measurement on the Space Shuttle that we performed in order to substantially extend the range of energies in which the elemental composition is known. We present and discuss the results, and we also summarize and discuss recent data on the electron component of cosmic rays. The body of data now available contains several features that are difficult to explain within current models of galactic shock acceleration and “leaky box” containment. We emphasize the need for further measurements, and we briefly discuss possible opportunities for future work.     

The charge-sign dependent effect in the solar modulation of cosmic rays

The centennial anniversary of the discovery of cosmic rays was in 2012. Since this discovery considerable progress has been made on several aspects related to galactic cosmic rays in the heliosphere. It is known that they encounter a turbulent solar wind with an imbedded heliospheric magnetic field when entering the Sun’s domain. This leads to significant global and temporal changes in their intensity inside the heliosphere, a process known as the solar modulation of cosmic rays. The prediction of a charge-sign dependent effect in solar modulation in the late 1970s and the confirmatory observational discoveries can also be considered as a milestone. A short review is given of these predictions based on theoretical and numerical modelling work, the observational confirmation and related issues.     

Radiation environment due to galactic and solar cosmic rays during manned mission to Mars in the periods between maximum and minimum solar activity cycles.

A possibility of a manned mission to Mars without exceeding the current radiation standards is very doubtful during the periods of minimum solar activity since the dose equivalent due to galactic cosmic rays exceeds currently recommended standards even inside a radiation shelter with an equivalent of 30 g cm-2 aluminum. The radiation situation at the time of maximum solar activity is determined by the occurrence of major solar proton events which are exceedingly difficult to forecast. This paper discusses the radiation environment during a manned mission to Mars in the years between minimum and maximum solar activity when the galactic cosmic ray intensity is considerably reduced, but the solar flare activity has not yet maximized.     

Observation of enhanced sub-iron (Sc---Cr) to iron ration in low energy cosmic rays of 50–100 MeV/N in spacelab-3

Relative abundances of sub-iron (Sc-Cr) to iron nuclei in low energy (50–100 MeV/N) galactic cosmic rays have been determined from an analysis of about 100 events of heavy ions (Z = 10−28) recorded in a detector assembly flown in the Anuradha cosmic ray experiment in the Spacelab-3 on a six day mission in April–May 1985. The measured abundance ratio of (Sc-Cr)/Fe nuclei in 50–100 MeV/N energy range is 1.1 ± 0.3, and the present result of enhanced ratio of sub-iron to iron nuclei is in agreement with other experimental results in 200–800 MeV/N range. The over-abundance of iron secondaries at these low energies cannot be explained in the conventional models for propagation of cosmic rays. Available experimental data indicate a very different time history for the low energy iron-group, as compared to those of lighter nuclei in galactic cosmic rays.     

Cosmic-ray composition and its relation to shock acceleration by supernova remnants

An overview is given on the present status of the understanding of the origin of galactic cosmic rays. Recent measurements of charged cosmic rays and photons are reviewed. Their impact on the contemporary knowledge about the sources and acceleration mechanisms of cosmic rays and their propagation through the Galaxy is discussed. Possible reasons for the knee in the energy spectrum and scenarios for the end of the galactic cosmic-ray component are described.