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.     

Mean energy of response to galactic cosmic ray spectrum: IMP 8 and Climax neutron monitor

Solar modulations of galactic cosmic ray (GCR) intensity contain a wealth of information about their transport in the heliosphere. To extract this information from the data one studies the dependence of the observed modulations on the mean energy of response of detectors providing data for the analyses. There is a great deal of confusion about the detector energy response to GCR spectrum in the literature. We present a preliminary report on the computations of the mean energy of response for the Climax neutron monitor (CL/NM) and IMP 8 cosmic ray nuclear composition instrument to GCR protons for 1973–1998, covering the solar cycles 21 and 22. We find that for penetrating proton channel on IMP 8 the mean energy changes by a factor of over two whereas for the neutron monitor the change is only 21%. However, the corresponding change for the computed modulation function is a factor of about 3.5.     

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.     

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.     

The transport of galactic and jovian cosmic ray electrons in the heliosphere

An overview is given on what we know about the cosmic ray diffusion process from the modelling of low-energy (MeV) electron transport in the heliosphere. For energies below ∼300 MeV, these electrons give a direct indication of the average mean free paths because they do not experience large adiabatic energy changes and their modulation is largely unaffected by global gradient and curvature drifts. Apart from galactic cosmic ray electrons, the jovian magnetosphere at ∼5 AU in the ecliptic plane is also a relatively strong source of MeV electrons, with energies up to ∼30 MeV. Therefore, when modelling the transport of these particles in the inner heliosphere, a three-dimensional treatment is essential. By comparing these models to observations from the Ulysses, Pioneer and Voyager missions, important conclusions can be made on e.g., the relative contributions of the galactic and jovian electrons to the total electron intensity, the magnitude of the parallel and perpendicular transport coefficients, and the time dependant treatment thereof.     

Regional and temporal variability of solar activity and galactic cosmic ray effects on the lower atmosphere circulation

In this work we studied the spatial and temporal structure of long-term effects of solar activity (SA) and galactic cosmic ray (GCR) variations on the lower atmosphere circulation as well as possible reasons for the peculiarities of this structure. The study revealed a strong latitudinal and regional dependence of SA/GCR effects on pressure variations in the lower troposphere which seems to be determined by specific features of baric systems formed in different regions. The temporal structure of SA/GCR effects on the troposphere circulation at high and middle latitudes is characterized by a roughly 60-year periodicity which is apparently due to the epochs of the large-scale atmospheric circulation. It is suggested that a possible mechanism of long-term effects of solar activity and cosmic ray variations on the troposphere circulation involves changes in the evolution of the polar vortex in the stratosphere of high latitudes, as well as planetary frontal zones.     

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.     

Modelling of galactic Carbon in an asymmetrical heliosphere: Effects of asymmetrical modulation conditions

Observations of galactic cosmic rays (GCRs) from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that in addition to a possible global asymmetry in the north–south dimensions (meridional plane) of the heliosphere, it is also possible that different modulation (turbulence) conditions could exist between the two hemispheres of the heliosphere. We focus on illustrating the effects on GCR Carbon of asymmetrical modulation conditions combined with a heliosheath thickness that has a significant dependence on heliolatitude. To reflect different modulation conditions between the two heliospheric hemispheres in our numerical model, the enhancement of both polar and radial perpendicular diffusion off the ecliptic plane is assumed to differ from heliographic pole to pole. The computed radial GCR intensities at polar angles of 55° (approximating the Voyager 1 direction) and 125° (approximating the Voyager 2 direction) are compared at different energies and for both particle drift cycles. This is done in the context of illustrating how different values of the enhancement of both polar and radial perpendicular diffusion between the two hemispheres contribute to causing differences in radial intensities during solar minimum and moderate maximum conditions. We find that in the A > 0 cycle these differences between 55° and 125° change both quantitatively and qualitatively for the assumed asymmetrical modulation condition as reflected by polar diffusion, while in the A < 0 cycle, minute quantitative differences are obtained. However, when both polar and radial perpendicular diffusion have significant latitude dependences, major differences in radial intensities between the two polar angles are obtained in both polarity cycles. Furthermore, significant differences in radial intensity gradients obtained in the heliosheath at lower energies may suggest that the solar wind turbulence at and beyond the solar wind termination shock must have a larger latitudinal dependence.     

On the effects of dynamical turbulence on the perpendicular diffusion of low-energy cosmic ray electrons

Diffusion perpendicular to the heliospheric magnetic field plays an integral role in the transport of charged particles in the heliosphere. In this study the perpendicular diffusion coefficient of low-energy cosmic ray electrons is calculated, using an equation derived from the random ballistic decorrelation interpretation of nonlinear guiding centre theory. An observationally motivated 2D turbulence power spectrum is assumed and the effects of various turbulence inputs on the resulting perpendicular diffusion coefficient are investigated. The perpendicular diffusion coefficients are first determined at 1 AU, for both magnetostatic and dynamical turbulence conditions. These solutions are also evaluated for radial distances of 0.1 AU to 10 AU to further investigate the values of the perpendicular diffusion coefficients in the very inner heliosphere. The results of this study show that the dissipation range of the turbulence power spectrum provides a negligible contribution towards the perpendicular diffusion coefficient, and that solutions derived using only the energy containing range serve as good approximations for solutions derived assuming the full 2D turbulence power spectrum. Finally, it is shown that the effects of dynamical turbulence, as considered in the present study, do not affect the perpendicular diffusion coefficients derived from the scattering theory considered here.     

50–200 GV cosmic rays at various heliolatitudes

We consider regular motion of 50 – 200 GV particles in a large-scale interplanetary magnetic field model which contains a wavy neutral sheet responsible for the sector-structure. Numerical calculations based upon energy losses along various trajectories are carried out to obtain the predicted omni-directional density and anisotropy of cosmic rays at various solar latitudes. A marked difference is found between odd and even solar cycles. The post-1969 field configuration gives small radial and large latitudinal gradient: cosmic ray density increases toward the poles. The latitudinal gradient turns out smaller and of opposite sense for the pre-1969 epoch. Anisotropy changes dramatically as we move off the solar equator: corotation appears to be restricted to low latitudes.     

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 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.     

27-day variations of the galactic cosmic ray intensity and anisotropy in the minimum of the 23rd solar activity cycle

We study the 27-day variations of the solar wind velocity, galactic cosmic ray (GCR) intensity and anisotropy in the last minimum epoch of solar activity (2007–2009, A < 0). The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy (A27A) in the current minimum epoch of solar activity (2007–2009, A < 0) is lesser than in previous positive polarity period as it is expected from the drift theory. So, polarity dependence rule for the 27-day variation of the GCR anisotropy is fully kept. It is a universal principle for the amplitudes of the 27-day variation of the GCR anisotropy. At the same time, the average amplitude of the 27-day variation of the GCR intensity (A27I) remains at the same level as for previous minimum epoch 1995–1997 (A > 0) showing by the same token an violation of its polarity dependence rule established earlier. We assume that this phenomenon could be generally related with the well established 27-day variation of the solar wind velocity being in anti-correlation with the similar changes of the 27-day variation of the GCR intensity. Generally, a character of the heliolongitudinal asymmetry of spatial large-scale structure of the solar wind velocity (SWV) established in the recent minimum epoch, preferentially pronounces in the behavior of the 27-day variation of the GCR intensity than anisotropy. The formation of the 27-day variation of the GCR anisotropy preferentially takes place in a restricted disk like local vicinity in the helioequatorial region, whilst the 27-day variation of the GCR intensity is formed in the global three dimensional vicinity of the heliosphere.     

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.     

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.     

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.     

Low energy galactic cosmic rays in the heliosphere

The flux of galactic cosmic rays (GCR) extends over a wide range of energies (from 10⁸ to 10²⁰ eV); it has a strong dependence on particle energy. Given the large span of energies the detection techniques, transport mechanisms and other characteristics vary as energy increases. In the low energy region (<10¹² eV) the flux of GCR is modulated by the solar activity. Continuous registers are necessary to study intensity variations that must have their origin in the Sun. Detectors were designed and constructed for the purpose, they operate since the middle of the last century providing valuable information to study recurrent periodicities and their relationship to those of solar phenomena, but also to elucidate whose are the relevant transport mechanisms inside the heliosphere. A brief review of the advancement in the comprehension of these phenomena is presented.     

Galactic cosmic ray radiation levels in spacecraft on interplanetary missions.

Using the Langley Research Center galactic cosmic ray (GCR) transport computer code (HZETRN) and the computerized anatomical man (CAM) model, crew radiation levels inside manned spacecraft on interplanetary missions are estimated. These radiation-level estimates include particle fluxes, LET (linear energy transfer) spectra, absorbed dose, and dose equivalent within various organs of interest in GCR protection studies. Changes in these radiation levels resulting from the use of various different types of shield materials are presented.