Heliospheric terminal shock acceleration and modulation of the anomalous cosmic-ray component

It is shown that a non-drift, spherically-symmetric, terminal shock acceleration model gives reasonable fits to the observed anomalous oxygen spectra for 1976 and 1982 only when a compression ratio of 2 is used together with a moderate radial dependence for K_rr. Preliminary results with an axially-symmetric, terminal shock acceleration model, including drifts, show that despite a change in polarity the accelerated spectra, at 50 AU in the equatorial plane, differ negligibly. However, the calculated modulation of anomalous oxygen at 2 AU is significantly less for a qA < 0 epoch (1987) than when qA > 0 (1976).     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

Galactic cosmic ray composition and energy spectra.

Galactic cosmic ray nuclei represent a significant risk to long-duration spaceflight outside the magnetosphere. We review briefly existing measurements of the composition and energy spectra of heavy cosmic ray nuclei, pointing out which species and energy ranges are most critical to assessing cosmic ray risks for spaceflight. Key data sets are identified and a table of cosmic ray abundances is presented for elements from H to Ni (Z = 1 to 28). Because of the 22-year nature of the solar modulation cycle, data from the approaching 1998 solar minimum is especially important to reducing uncertainties in the cosmic ray radiation hazard. It is recommended that efforts to model this hazard take advantage of approaches that have been developed to model the astrophysical aspects of cosmic rays.     

Cosmic rays in the dynamic heliosheath

Cosmic ray modulation in the outer heliosphere is discussed from a modeling perspective. Emphasis is on the transport and acceleration of these particles at and beyond the solar wind termination shock in the inner heliosheath region and how this changes over a solar cycle. We will show that by using numerical models, and by comparing results to spacecraft observations, much can be learned about the dependence of cosmic ray modulation on solar cycle changes in the solar wind and heliospheric magnetic field. While the first determines the heliospheric geometry and shock structure, the latter results in a time-dependence of the transport coefficients. Depending on energy, both these effects contribute to cosmic ray intensities in the inner heliosheath changing over a solar cycle.     

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.     

Effects of the position of the solar wind termination shock and the heliopause on the heliospheric modulation of cosmic rays

The effects of changing the position of the solar wind termination shock and the position of the heliopause, and therefore the extent of the heliosheath, on the modulation of cosmic ray protons are illustrated. An improved numerical model with diffusive termination shock acceleration, a heliosheath and drifts is used. The modulation is computed in the equatorial plane and at 35 heliolatitude using recently derived diffusion coefficients applicable to a number of cosmic ray species during both magnetic polarity cycles of the Sun. It was found that qualitatively the modulation results for the different heliopause positions are similar although they differ quantitatively, e.g., clearly different radial gradients are predicted for the regions beyond the termination shock compared to inside the shock. The difference between the modulation for the two solar polarity cycles are less significant at a heliolatitude of 35° than in the equatorial plane. We found that moving the termination shock from 90 to 100 AU, with the heliopause fixed at 120 AU, caused only quantitative differences so that the exact position of the TS in the outer heliosphere seems not crucially important to global modulation. Moving the heliopause outwards, to represent the modulation in the tail region of the heliosphere, causes overall decreases in the cosmic ray intensities but not linearly as a function of energy, e.g., at 1 GeV the effect is insignificant. We conclude from this modelling that the modulation of protons in the heliospheric nose and tail regions are qualitatively similar although, clear quantitative and interesting differences occur.     

The dynamic heliosphere,solar activity,and cosmic rays

This brief review addresses the relation between solar activity, cosmic ray variations and the dynamics of the heliosphere. The global features of the heliosphere influence what happens inside its boundaries on a variety of time-scales. Galactic and anomalous cosmic rays are the messengers that convey vital information on global heliospheric changes in the manner that they respond to these changes. By observing cosmic rays over a large range of energies at Earth, and with various space detectors, a better understanding is gained about space weather and climate. The causes of the cosmic ray variability are reviewed, with emphasis on the 11-year and 22-year cycles, step modulation, charge-sign dependent modulation and particle drifts. Advances in this field are selectively discussed in the context of what still are some of the important uncertainties and outstanding issues.     

Galactic cosmic radiation model and its applications.

A model for the differential energy spectra of galactic cosmic radiation as a function of solar activity is described. It is based on the standard diffusion-convection theory of solar modulation. Estimates of the modulation potential based on fitting this theory to observed spectral measurements from 1954 to 1989 are correlated to the Climax neutron counting rates and to the sunspot numbers at earlier times taking into account the polarity of the interplanetary magnetic field at the time of observations. These regression lines then provide a method for predicting the modulation at later times. The results of this model are quantitatively compared to a similar Moscow State University (MSU) model. These model cosmic ray spectra are used to predict the linear energy transfer spectra, differential energy spectra of light (charge < or = 2) ions, and single event upset rates in memory devices. These calculations are compared to observations made aboard the Space Shuttle.     

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.     

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.     

Prediction of expected global climate change by forecasting of galactic cosmic ray intensity time variation in near future based on solar magnetic field data

A method of prediction of expected part of global climate change caused by cosmic ray (CR) by forecasting of galactic cosmic ray intensity time variation in near future based on solar activity data prediction and determined parameters of convection-diffusion and drift mechanisms is presented. This gave possibility to make prediction of expected part of global climate change, caused by long-term cosmic ray intensity variation. In this paper, we use the model of cosmic ray modulation in the Heliosphere, which considers a relation between long-term cosmic ray variations with parameters of the solar magnetic field. The later now can be predicted with good accuracy. By using this prediction, the expected cosmic ray variations in the near Earth space also can be estimated with a good accuracy. It is shown that there are two possibilities: (1) to predict cosmic ray intensity for 1–6 months by using a delay of long-term cosmic ray variations relatively to effects of the solar activity and (2) to predict cosmic ray intensity for the next solar cycle. For the second case, the prediction of the global solar magnetic field characteristics is crucial. For both cases, reliable long-term cosmic ray and solar activity data as well as solar magnetic field are necessary. For solar magnetic field, we used results of two magnetographs (from Stanford and Kitt Peak Observatories). The obtained forecasting of long-term cosmic ray intensity variation we use for estimation of the part of global climate change caused by cosmic ray intensity changing (influenced on global cloudiness covering).     

Time-dependent cosmic ray modulation

Time-dependent cosmic ray modulation is calculated over multiple solar cycles using our well established two-dimensional time-dependent modulation model. Results are compared to Voyager 1, Ulysses and IMP cosmic ray observations to establish compatibility. A time-dependence in the diffusion and drift coefficients, implicitly contained in recent expressions derived by , ,  and , is incorporated into the cosmic ray modulation model. This results in calculations which are compatible with spacecraft observations on a global scale over consecutive solar cycles. This approach compares well to the successful compound approach of Ferreira and Potgieter (2004). For both these approaches the magnetic field magnitude, variance of the field and current sheet tilt angle values observed at Earth are transported time-dependently into the outer heliosphere. However, when results are compared to observations for extreme solar maximum, the computed step-like modulation is not as pronounced as observed. This indicates that some additional merging of these structures into more pronounced modulation barriers along the way is needed.     

The role of radial perpendicular diffusion and latitude dependent acceleration along the solar wind termination shock

A numerical model, based on Parker’s transport equation, describing the modulation of anomalous cosmic rays and containing diffusive shock acceleration is applied. The role of radial perpendicular diffusion at the solar wind termination shock, and as the dominant diffusion coefficient in the outer heliosphere, is studied, in particular the role it plays in the effectiveness of the acceleration of anomalous protons and helium when its latitude dependence is changed. It is found that the latitudinal enhancement of radial perpendicular diffusion towards the heliospheric poles and along the termination shock has a prominent effect on the acceleration of these particles. It results in a ‘break’ in the energy spectrum for anomalous protons at ∼6.0 MeV, causing the spectral index to change from E^−1.38 to E^−2.23, but for anomalous helium at ∼3.0 MeV, changing the spectral index from E^−1.38 to E^−2.30. When approaching the simulated TS, the changes in the modulated spectra as they unfold to a ‘steady’ power law shape at energies below 50 MeV are much less prominent as a function of radial distances when radial perpendicular diffusion is increased with heliolatitude.     

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.     

Cosmic ray modulation in the outer heliosphere: Predictions for cosmic ray intensities up to the heliopause along Voyager 1 and 2 trajectories

Time dependent cosmic ray modulation in the outer heliosphere is calculated and results are compared to Voyager 1 and 2 observations using a two-dimensional time-dependent cosmic ray transport model. We predict possible future 133–242 MeV proton observations along the Voyager 1 and 2 spacecraft trajectories. Recent theoretical advances in cosmic ray transport parameters are introduced in order to provide a time-dependence for the assumed transport parameters used in the model. This leads to results that are in general compatible with the spacecraft observations in the inner and outer heliosphere over multiple solar cycles. However, for the outer heliosphere, we find that the Voyager 1 and 2 spacecraft observations cannot be fitted with an identical set of parameters along both trajectories. This indicates a possible asymmetric heliosphere or a symmetric heliosphere but with different diffusion parameters in the northern and southern hemispheres, respectively. Furthermore, results indicate that Voyager 2 observations are still under the influence of solar cycle related changes because of the large modulation volume between the heliopause and spacecraft location in contrast to Voyager 1 which shows a steady increase in cosmic ray intensities.     

Modulation of cosmic rays: Perpendicular diffusion and drifts in a heliosphere with a solar wind termination shock

In this study the roles of polar perpendicular diffusion and drifts are illustrated in a model containing a heliosheath and diffusive shock acceleration as applied to the solar wind termination shock. Of particular interest is the relation of polar perpendicular diffusion to particle drifts and how the effectiveness of the termination shock acceleration of galactic and anomalous protons is influenced by this relation. We found that drifts have a more prominent effect than the polar enhancement of perpendicular diffusion so that its omission from termination shock models would produce unrealistically large shock acceleration and consequently also larger modulation effects throughout the heliosphere. The computed spectra at a heliolatitude of 35° are almost similar for the two polarity epochs indicating that the two Voyager spacecraft might not observe differences between the two cycles in future.     

Rigidity spectrum of the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity

We study the temporal evolution of the power rigidity spectrum of the first (27 days) and the second (14 days) harmonics of the 27-day variation of the galactic cosmic ray intensity measured by neutron monitors in the period of 1965–2002. The rigidity spectrum of these variations can be approximated by a power law. We show the rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity have similar time profiles. These spectra are hard (γ ≈ 0.5 ± 0.1) and soft (γ ≈ 1.1 ± 0.2) during solar maximum and minimum activity, respectively. We ascribe this to the alternation of the sizes of the modulation regions responsible for the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity. Especially, the average radial sizes of the modulation regions of the 27-day variation and the heliolatitudinal extension of the heliolongitudinal asymmetry are smaller during solar minimum than during solar maximum. We show also, that the temporal changes of the power rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity are in a negative correlation with the changes of the rigidity spectrum of the corresponding 11-year variation.     

A transient MHD model applicable for the source of solar cosmic ray acceleration

A two-dimensional, time-dependent magnetohydrodynamic (MHD) model is used to describe the possible mechanisms for the source of solar cosmic ray acceleration following a solar flare. The hypothesis is based on the propagation of fast mode MHD shocks following a sudden release of energy. This model has already been used with some success for simulation of some major features of type II shocks and white light coronal transients. In this presentation, we have studied the effects of initial magnetic topology and strength on the formation of MHD shocks. We consider the plasma beta (thermal pressure/magnetic pressure) as a measure of the $\text{initial}$, relative strength of the field. During dynamic mass motion, the Alfvén Mach number is the more appropriate measure of the magnetic field's ability to control the outward motion. We suggest that this model (computed self-consistently) provides the shock wave and the disturbed mass motion behind it as likely sources for solar cosmic ray acceleration.     

Modeling ground and space based cosmic ray observations

After entering our local astrosphere (called the heliosphere), galactic cosmic rays, as charged particles, are affected by the Sun’s turbulent magnetic field. This causes their intensities to decrease towards the inner heliosphere, a process referred to as modulation. Over the years, cosmic ray modulation has been studied extensively at Earth, utilizing both ground and space based observations. Moreover, modelling cosmic ray modulation and comparing results with observations, insight can be gained into the transport of these particles, as well as offering explanations for observed features. We review some of the most prominent cosmic ray observations made near Earth, how these observations can be modelled and what main insights are gained from this modelling approach. Furthermore, a discussion on drifts, as one of the main modulation processes, are given as well as how drift effects manifest in near Earth observations. We conclude by discussing the contemporary challenges, fuelled by observations, which are presently being investigated. A main challenge is explaining observations made during the past unusual solar minimum.