Modulation of galactic cosmic rays in a north–south asymmetrical heliosphere

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock. Voyager 1 crossed this termination shock at ∼94 AU in 2004, while Voyager 2 crossed it in 2007 at a different heliolatitude, about 10 AU closer to the Sun. These different positions of the termination shock confirm the dynamic and cyclic nature of the shock’s position. Observations from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that apart from the dynamic nature caused by changing solar activity there also may exist a global asymmetry in the north–south (polar) dimensions of the heliosphere, in addition to the expected nose–tail asymmetry. This relates to the direction in which the heliosphere is moving in interstellar space and its orientation with respect to the interstellar magnetic field. In this paper we focus on illustrating the effects of this north–south asymmetry on modulation of galactic cosmic ray Carbon, between polar angles of 55° and 125°, using a numerical model which includes all four major modulation processes, the termination shock and the heliosheath. This asymmetry is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. When comparing the computed spectra between the two polar angles, we find that at energies E < ∼1.0 GeV the effects of the assumed asymmetry on the modulated spectra are insignificant up to 60 AU from the Sun but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > ∼1.0 GeV, these effects remain insignificant throughout the heliosphere even very close to the heliopause. Furthermore, we find that a higher local interstellar spectrum for Carbon enhances the effects of asymmetric modulation between the two polar angles at lower energies (E < ∼300 MeV). In conclusion, it is found that north–south asymmetrical effects on the modulation of cosmic ray Carbon depend strongly on the extent of the geometrical asymmetry of the heliosheath together with the assumed value of the local interstellar spectrum.     

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.     

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.     

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.     

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.     

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.     

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.     

Effects of various dissipation range onset models on the 26-day variations of low-energy galactic cosmic-ray electrons

The effect of various models presented by Leamon et al. (2000) for the dissipation range cutoff wavenumber on the 26-day variations of galactic cosmic-ray electrons in a Fisk-Parker hybrid field is investigated, by means of a three-dimensional steady-state numerical modulation code. Analytical expressions for the mean free paths parallel and perpendicular to the heliospheric magnetic field are adapted from the works of 31 and 28, respectively. Note that only solar minimum conditions are considered, and that only qualitative agreement with data is sought. Effective diffusion for galactic electrons pertaining to 26-day variations is found to be dominated by the ratio of the perpendicular to parallel mean free paths at low energies, and the relationship between changes in cosmic-ray intensities and the modulation parameter postulated by Zhang (1997) is found to no longer hold when this ratio drops below a critical value. Use of ion inertial scale dependent models for the dissipation range cutoff leads to possible second linearities in the relative amplitudes as functions of latitude gradient.     

Viewing radiation signatures of solar energetic particles in interplanetary space

A current serious limitation on the studies of solar energetic particle (SEP) events is that their properties in the inner heliosphere are studied only through in situ spacecraft observations. Our understanding of spatial distributions and temporal variations of SEP events has come through statistical studies of many such events over several solar cycles. In contrast, flare SEPs in the solar corona can be imaged through their radiative and collisional interactions with solar fields and particles. We suggest that the heliospheric SEPs may also interact with heliospheric particles and fields to produce signatures which can be remotely observed and imaged. A challenge with any such candidate signature is to separate it from that of flare SEPs. The optimum case for imaging high-energy (E > 100 MeV) heliospheric protons may be the emission of π⁰-decay γ-rays following proton collisions with solar wind (SW) ions. In the case of E > 1 MeV electrons, gyrosynchrotron radio emission may be the most readily detectible remote signal. In both cases we may already have observed one or two such events. Another radiative signature from nonthermal particles may be resonant transition radiation, which has likely already been observed from solar flare electrons. We discuss energetic neutrons as another possible remote signature, but we rule out γ-ray line and 0.511 MeV positron annihilation emission as observable signatures of heliospheric energetic ions. We are already acquiring global signatures of large inner-heliospheric SW density features and of heliosheath interactions between the SW and interstellar neutral ions. By finding an appropriate observable signature of remote heliospheric SEPs, we could supplement the in situ observations with global maps of energetic SEP events to provide a comprehensive view of SEP events.     

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.     

MHD modeling of the outer heliosphere: Achievements and challenges

An overview is presented of magnetic-field-related effects in the solar wind (SW) interaction with the local interstellar medium (LISM) and the different theoretical approaches used in their investigation. We discuss the possibility that the interstellar magnetic field (ISMF) introduces north–south and east–west asymmetries of the heliosphere, which might explain observational data obtained by the Voyager 1 and Voyager 2 spacecraft. The SW–LISM interaction parameters that are responsible for the deflection of the interstellar neutral hydrogen flow from the direction of propagation of neutral helium in the inner heliosheath are outlined. The possibility of a strong ISMF, which increases the heliospheric asymmetry and the H–He flow deflection, is discussed. The effect of the combination of a slow-fast solar wind during solar minimum over the Sun’s 11-year activity cycle is illustrated. The consequences of a tilt between the Sun’s magnetic and rotational axes are analyzed. Band-like areas of an increased magnetic field distribution in the outer heliosheath are sought in order to discover regions of possible 2–3 kHz radio emission.     

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.     

Where does the heliospheric modulation of galactic cosmic rays start?

The long outstanding question of where the heliospheric (solar) modulation of galactic cosmic rays actually begins, in terms of spatial position, as well as at what high kinetic energy, can now be answered. Both answers are possible by using the results of an advanced numerical model, together with appropriate observations. Voyager 1 has been exploring the outskirts of the heliosphere and is presently entering what can be called the very local interstellar medium. It has been generally expected, and accepted, that once the heliopause is crossed, the local interstellar spectrum (LIS) should be measured in situ by the Voyager spacecraft. However, we show that this may not be the case and that modulation effects on galactic cosmic rays can persist well beyond the heliopause. For example, proton observations at 100 MeV close to the heliopause can be lower by ∼$\sim$25% to 40% than the LIS, depending on solar modulation conditions. It is also illustrated quantitatively that significant solar modulation diminishes above ∼$\sim$50 GeV at Earth. It is found that cosmic ray observations above this energy contain less that 5%$5\%$ solar modulation effects and should therefore reflect the LIS for galactic cosmic rays. Input spectra, in other words the very LIS, for solar modulation models are now constrained by in situ observations and can therefore not any longer be treated arbitrarily. It is also possible for the first time to determine the lower limit of the very LIS from a few MeV/nuc to very high energies.     

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.     

The heliospheric modulation of 3–10 MeV electrons: Modeling of changes in the solar wind speed in relation to perpendicular polar diffusion

The discrepancy between cosmic ray model predictions representing solar minimum conditions in the heliosphere and the 3–10 MeV post-1998 electrons observations by the Kiel Electron Telescope (KET) onboard Ulysses suggests the need for consistent changes in model parameters with increasing solar activity. In order to reduce this discrepancy, an effort is made to model the KET observations realistically during periods of increased solar activity by applying an advanced three-dimensional, steady-state electron modulation model based on Parker’s transport equation including the Jovian electron source. Some elements of the diffusion tensor which were not previously emphasized are revisited. A new relation is also established between the latitudinal dependence of the solar wind speed and the perpendicular polar diffusion. Based on this relation, a transition of an average solar wind speed from solar minimum to solar maximum conditions, as observed on board the Ulysses spacecraft, is modeled on the concept of the time-evolution of large polar coronal holes. These changes are correlated to different scenarios of the enhancement of perpendicular polar diffusion. Effects of these scenarios are illustrated, as a series of steady-state solutions, on the computed 7 MeV Jovian and galactic electrons in comparison with 3–10 MeV electrons observed from the period 1998 to the end of 2003. It is shown that this approach improves compatibility with the KET observations but it also points to the need for a time-dependent electron modulation model to fully describe modulation during moderate to extreme solar maximum conditions.     

Detection of the high energy component of Jovian electrons in Low Earth Orbit with the PAMELA experiment

The PAMELA experiment is devoted to the study of cosmic rays in Low Earth Orbit with an apparatus optimized to perform a precise determination of the galactic antimatter component of c.r. It is constituted by a number of detectors built around a permanent magnet spectrometer. PAMELA was launched in space on June 15th 2006 on board the Russian Resurs-DK1 satellite for a mission duration of 3 years. The characteristics of the detectors, the long lifetime and the orbit of the satellite, will allow to address several aspects of cosmic-ray physics. In this work we discuss the observational capabilities of PAMELA to detect the electron component above 50 MeV. The magnetic spectrometer allows a detailed measurement of the energy spectrum of electrons of galactic and Jovian origin. Long term measurements and correlations with Earth–Jupiter 13 months synodic period will allow to separate these two contributions and to measure the primary electron Jovian component, dominant in the 50–70 MeV energy range. With this technique it will also be possible to study the contribution to the electron spectrum of Jovian e⁻ reaccelerated up to 2 GeV at the Solar Wind Termination Shock.     

High energy charged particle experiment – A payload proposal to KuaFu-B mission

High Energy Charged Particle Experiment (HECPE) is to measure the fluxes of MeV electrons and tens of MeV protons. The two satellites of KuaFu-B are in the same polar orbit with apogee 7.0R_E, perigee 1.8R_E. They can sweep large L values and pass through the inner and outer radiation belts. The high energy electrons and protons in the radiation belts are principal sources for failures of satellites and spacecrafts in the Earth orbits. The enhancements of the high energy electrons and protons, so-called energetic particle events, are important phenomena of the Space Weather. The energy ranges monitored by HECPE are 0.3–0.5 MeV, 0.5–1.0 MeV, 1.0–2.0 MeV, and E > 2.0 MeV for electrons, 5–10 MeV, 10–20 MeV, 20–40 MeV, and 40–80 MeV for protons.     

On the transport of cosmic rays in the distant heliosheath

This paper discusses the transport of energetic charged particles through a sectored magnetic field in distant regions of the inner heliosheath. As the plasma flow slows down on approach to the stagnation point on the heliopause, the distance between the folds of the current sheet decreases to the point where it becomes comparable to the cyclotron radius of a cosmic ray particle. Under these conditions a particle can effectively drift across the stack of magnetic sectors with a speed comparable with the particle’s velocity. For a random distribution of current sheet separation distances, a diffusive transport across the stack of sectors occurs instead. The proposed mechanism could have contributed to unusually high intensities of galactic cosmic rays measured by Voyager 1 in the heliosheath during 2009–2010.