Solar wind structure in the outer heliosphere

A solar wind parcel evolves as it moves outward, interacting with the solar wind plasma ahead of and behind it and with the interstellar neutrals. This structure varies over a solar cycle as the latitudinal speed profile and current sheet tilt change. We model the evolution of the solar wind with distance, using inner heliosphere data to predict plasma parameters at Voyager. The shocks which pass Voyager 2 often have different structure than expected; changes in the plasma and/or magnetic field do not always occur simultaneously. We use the recent latitudinal alignment of Ulysses and Voyager 2 to determine the solar wind slowdown due to interstellar neutrals at 80 AU and estimate the interstellar neutral density. We use Voyager data to predict the termination shock motion and location as a function of time.     

Possible physical properties of the flow in the vicinity of the heliopause

An interface between the fully ionized hydrogen plasma of the solar wind (SW) and the partially ionized hydrogen gas flow of the local interstellar medium (LISM) is formed as a region where there is a strong interaction between these two flows. The interface is bounded by the solar wind termination shock (TS) and the LISM bow shock (BS) and is separated on two regions by the heliopause (HP) separating the solar wind and charged component of the LISM (plasma component below). The BS is formed due to the deceleration of the supersonic LISM flow relative to the solar system. Regions of the interface between the TS and HP and between the HP and BS were in literature named as the inner and outer heliosheaths, respectively. An investigation of the structure and physical properties of the heliosheath is at present especially interested due to the fact that Voyager-1 and Voyager-2 have crossed the TS in December 2004 (Burlaga, L.F., Ness, N.F., Acuna, M.Y., et al. Crossing the termination shock into the the heliosheath. Magnetic fields. Science 309, 2027–2029, 2005; Fisk, L.A. Journey into the unknown beyond. Science 309, 2016–2017, 2005; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Voyager 1 in the foreshock, termination shock and heliosheath. Science 309, 2020–2024, 2005; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. Voyager 1 explores the termination shock region and the heliosheath beyond. Science 309, 2017–2020, 2005) and in September 2007 (Jokipii, J.R. A shock for Voyager 2. Nature 454, 38–39, 2008; Gurnett, D.A., Kurth, W.S. Intense plasma waves at and near the solar wind termination shock. Nature 454, 78–80, 2008. doi: 10.1038/nature07023; Wang, L., Lin, R.P., Larson, D.E., Luhmann, J.G. Domination of heliosheath pressure by shock-accelerated pickup ions from observations of neutral atoms. Nature 454, 81–83, 2008. doi: 10.1038/nature07068.14; Burlaga, L.F., Ness, N.F., Acuna, M.H., et al. Magnetic fields at the solar wind termination shock. Nature 454, 75–77, 2008. doi: 10.1038/nature07029; Richardson, J.D., Kasper, J.C., Wang, C., et al. Cool heliosheath plasma and deceleration of the upstream solar wind at the termination shock. Nature 454, 63–66, 2008. doi: 10.1038/nature07024; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. An asymmetric solar wind termination shock. Nature 454, 71–74, 2008. doi: 10.1038/nature07022; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Mediation of the solar wind termination shock by non-thermal ions. Nature 454, 67–70, 2008. doi: 10.1038/nature 07030), respectively, and entered to the inner heliosheath.     

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.     

Heliospheric asymmetries due to the action of the interstellar magnetic field

We discuss the asymmetry of the heliospheric discontinuities obtained from the analysis of 3D modeling of the solar wind (SW) interaction with local interstellar medium (LISM). The flow of charged particles is governed by the ideal MHD equations and the flow of neutral particles is described by the Boltzmann equation. The emphasis is made on the asymmetries of the termination shock (TS) and the heliopause under the combined action of the interstellar and interplanetary magnetic fields (ISMF and IMF) in the presence of neutral hydrogen atoms whose transport through the heliosphere is modeled kinetically, using a Monte Carlo approach. We show that the deflection of neutral hydrogen flow from its original direction in the unperturbed LISM is highly anisotropic and evaluate a possible angle between the hydrogen deflection plane measured in the SOHO SWAN experiment and the plane containing the ISMF and LISM velocity vectors for different ISMF strengths. It is shown that the ISMF of a strength greater than 4 μG can account for the 10 AU difference in the TS heliocentric difference observed during its crossing by the Voyager 1 and Voyager 2 spacecraft, which however results in a larger discrepancy between the calculated and observed velocity distributions. The effect of a strong ISMF on the distribution of plasma quantities in the inner heliosheath and on 2–3 kHz radio emission is discussed.     

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.     

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.     

On the dynamics of the heliosphere on intermediate and long time-scales

A time-dependent, three-dimensional model of the dynamics of the heliosphere as a result of solar activity and a time-varying local interstellar medium is presented. The model is based on a recent version of the well known ZEUS code and employs parallel processing. It includes the solar and interstellar plasma components as well as neutral atoms, and contains the heliospheric magnetic field in a kinematic fashion. We study the dynamics of the heliosphere due to solar activity on periods of months to years up to the so-called Schwabe (11-year) cycle as well as due to time variations of the local interstellar medium, all of which have drawn increasing attention during recent years, as the significance of their direct or indirect effect on the Earth and its environment is under lively debate.     

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.     

Spectral properties of backscattered solar Ly-α radiation in the heliosphere: A search for heliospheric boundary effects

We explore the imprints of the solar wind interaction with the local interstellar medium on spectral properties of backscattered solar Ly-α radiation. We employ a newly developed effective model for the interstellar H atom velocity distribution function in the heliosphere ( Katushkina and Izmodenov, 2011). The model takes into account both global effects of H atom perturbations at the heliospheric boundaries and local (i.e. within 10–20 AU from the Sun) effects of the solar ionization, charge exchange, solar gravitation and radiation pressure.     

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.     

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.     

Dynamic processes in the outer heliosphere: voyager observations and models

We review recent Voyager 2 observations in the vicinity of 70 AU. The character of the solar wind plasma data between 2002 and 2003 changed to a regime in which the speed, density and magnetic field magnitude are positively correlated. The average speed of the solar wind at Voyager 2 increased between early 2003 and mid-2004, which we attribute to a return of fast coronal hole flow. We use solar wind data at Earth as input to numerical models which include the effect of pickup ions to model the radial evolution of the solar wind. The model reproduces the basic features of the observations. As a specific example, we investigate the propagation of the Halloween (Oct.–Nov.), 2003 storms in the outer heliosphere. The model predictions are in reasonable agreement with Voyager 2 observations.     

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.     

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.     

Coupling of the interstellar and interplanetary magnetic fields at the heliospheric interface: The effect of neutral hydrogen atoms

We present numerical results showing the effect of neutral hydrogen atoms on the solar wind (SW) interaction with the local interstellar medium (LISM), where the interstellar magnetic field (ISMF) is coupled to the interplanetary magnetic field (IMF) at the surface of the heliopause. The IMF on the inner boundary surrounding the Sun is specified in the form of a Parker spiral and self-consistently develops in accordance with the SW motion inside the heliopause. The model of the SW–LISM interaction involves both plasma and neutral components which are treated as fluids. The configuration is chosen where the ISMF is orthogonal to the LISM velocity and tilted 60° to the ecliptic plane. This orientation of the magnetic field is a possible explanation of the 2–3 kHz emission data which is believed to originate ahead of the heliopause. It is shown that the topology of the heliospheric current sheet is substantially affected by the ISMF. The interaction pattern dependence on the neutral hydrogen density is analyzed.     

The Compton–Getting effect of energetic particles with an anisotropic pitch-angle distribution

This paper presents a simulation of anisotropy measurements by the low-energy charged particle (LECP) experiment on Voyager 1 for cases when the particle pitch-angle distribution function in the solar wind plasma reference frame is not isotropic. The model includes both the Compton–Getting anisotropy and perpendicular diffusion anisotropy that possibly exists in the upstream region of the termination shock. The results show that the Voyager 1 data cannot rule out either the model with zero solar wind speed or the one with a finite speed on qualitative basis. The determination of solar wind speed using the Compton–Getting effect is affected by the assumption of the magnetic field direction and perpendicular diffusion anisotropy. Because the pitch-angle distribution anisotropy is so large, a small uncertainty in the magnetic field direction can produce very different solar wind speeds ranging from ∼0 to >400 km/s. In fact, if the magnetic field is chosen to be in the Parker spiral direction, which is consistent with the magnetometer measurement on Voyager 1, the derived solar wind speed is still close to the supersonic value. Only the two lowest-energy channels of the LECP instrument may give a definitive answer to the solar wind speed. However, because these channels contain a very high level of cosmic ray background, an uncertainty of just a few percent in the background can entirely hamper the estimate of solar wind speed.     

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.     

Saturn kilometric radiation as a monitor for the solar wind?

Since the Voyager mission it is known that Saturn Kilometric Radiation (SKR) is strongly influenced by external forces, i.e., the solar wind and in particular the solar wind ram pressure. Recent studies using Cassini data essentially confirmed these findings for particular periods during the first Cassini orbit of Saturn. The data coverage of SKR by the Cassini/RPWS experiment for the period of six months prior to Saturn Orbit Insertion (July 1, 2004) is rather continuous, whereas there are gaps in the solar wind plasma data. The strong correlation of SKR with the solar wind may provide an indication on the variations of the solar wind plasma, specifically during the gap periods. These periods lacking solar wind data are substituted by Ulysses solar wind data which have been propagated over ∼4 AU, applying magnetohydrodynamic propagation models. Cross correlation studies showed that Ulysses solar wind data can be taken as a substitute for missing Cassini data. The use of SKR as monitor for solar wind variations is discussed. With the present set of observations the SKR proxy lacks significant reliability.     

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.