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.     

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.     

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.     

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.     

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.     

Modeling of the outer heliosphere with the realistic solar cycle

Time-dependent kinetic-continuum model of the solar wind interaction with the two-component local interstellar cloud (LIC) has been developed recently [Izmodenov, V., Malama, Y.G., Ruderman, M.S. Solar cycle influence on the interaction of the solar wind with local interstellar cloud. Astron. Astrophys. 429, 1069–1080, 2005a.]. Here, we adopted this model to the realistic solar cycle, when the solar wind parameters at the Earth’s orbit are taken from space data. This paper focuses on the results related to the termination shock (TS) excursion with the solar cycle that may help to understand Voyager 1 data obtained at and after the crossing of the termination shock and to predict the time of the TS crossing by Voyager 2.     

Heliopause stability revisited: dispersion analysis and numerical simulations

The heliopause, a surface separating the tenuous hot heliosheath flow and the dense, strongly magnetized interstellar flow, is subject to instabilities of the Rayleigh–Taylor and Kelvin–Helmholtz types. The dynamic evolution of this discontinuity is of considerable importance for understanding the neutral atom and cosmic-ray filtration at the interface. Here, we investigate the stability of the upwind heliopause in the presence of charge exchange collisions using both an analytic (dispersion relation) approach and a numerical model that includes the interstellar magnetic field. The linear analysis yields a cubic dispersion relation that admits imaginary solutions for the full range of wavenumbers, implying that the stagnation point on the heliopause is unconditionally Rayleigh–Taylor unstable to small perturbations propagating parallel to the discontinuity surface. We confirm this result by following the nonlinear development of the instability with a time-dependent simulation using a four fluid MHD-neutral numerical code. For the typical solar wind and LISM conditions, we obtain cyclical evolution of the upwind heliopause with a period of the order of 100 years. We also identify two areas of space physics where the instability may have important implications.     

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.     

α–ω Effect and peculiarities of the 27-day variations of the galactic cosmic rays intensity,solar wind and solar activity parameters

We show that the amplitudes of the 27-day variations of galactic cosmic ray (GCR) intensity, solar wind and solar activity parameters have a periodicity with duration of three to four Carrington rotation periods (3–4 CRP). We assume that the general reason for this phenomenon may be related to similar cyclicity of topological structure of the solar magnetic field lines created owing to the asymmetry of turbulent solar dynamo and solar differential rotation transforming the Sun’s poloidal magnetic field to the toroidal (α–ω effect), and vice versa.     

Modelling of low-energy galactic electrons in the heliosheath

The modulation of cosmic ray electrons in the heliosphere plays an important role in improving our understanding and assessment of the processes applicable to low-energy galactic electrons. A full three-dimensional numerical model based on Parker’s transport equation is used to study the modulation of 10 MeV galactic electrons, in particular inside the heliosheath. The emphasis is placed on the role that perpendicular diffusion plays in causing the extraordinary large increase in the observed intensities of these electrons in the heliosheath. The modelling is compared with observations of 6–14 MeV electrons from the Voyager 1 mission. Results are shown for the radial intensity profiles of these electrons, as well as the modulation effects of varying the extent of the heliosheath by changing the location of the termination shock and the heliopause and the value of the local interstellar spectrum. We confirm that the heliosheath acts as a modulation ‘barrier’ for low-energy galactic electrons. The significance of this result depends on how wide the inner heliosheath is; on how high the very local interstellar spectrum is at these low energies (E < 100 MeV) and on how small perpendicular diffusion is inside the inner heliosheath.     

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.     

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.     

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.     

Dynamical characterization of the last prolonged solar minima

The planetary hypothesis of the solar cycle is an old idea in which the gravitational influence of the planets has a non-negligible effect on the causes of the solar magnetic cycle. The advance of this hypothesis is based on phenomenological correlations between dynamical parameters of the Sun’s movement around the barycentre of the Solar System and sunspots time series; and more especially, identifying relationships linking solar barycentric dynamics with prolonged minima (especially Grand Minima events). However, at present there is no clear physical mechanism relating these phenomena. The possible celestial influence on solar cycle modulation is of great importance not only in solar physics but also in Earth sciences, because prolonged solar minima have associated important climatic and telluric variations, in particular, during the Maunder and Dalton Minimum. In this work we looked for a possible causal link in relation with solar barycentric dynamics and prolonged minima events. We searched for particular changes in the Sun’s acceleration and concentrated on long-term variations of the solar cycle. We show how the orbital angular momentum of the Sun evolves and how the inclination of the solar barycentric orbit varies during the epochs of orbital retrogressions. In particular, at these moments, the radial component of the Sun’s acceleration (i.e., in the barycentre-Sun direction) had an exceptional magnitude. These radial impulses occurred at the very beginning of the Maunder Minimum, during the Dalton Minimum and also at the maximum of cycle 22 before the present extended minimum. We also found a strong correlation between the planetary torque and the observed sunspots international number around that maximum. We apply our results in a novel theory of Sun–planets interaction that it is sensitive to Sun barycentric dynamics and found a very important effect on the Sun’s capability of storing hypothetical reservoirs of potential energy that could be released by internal flows and might be related to the solar cycle. This process begins about 40 years before the solar angular momentum inversions, i.e., before Maunder Minimum, Dalton Minimum, and before the present extended minimum. Our conclusions suggest a dynamical characterization of peculiar prolonged solar minima. We discuss the possible implications of these results for the solar cycle including the present extended minimum.     

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.     

Saturn radio emission and the solar wind: Voyager-2 studies

Voyager 2 data from the Plasma Science experiment, the Magnetometer experiment and the Planetary Radio Astronomy experiment were used to analyze the relationship between parameters of the solar wind/interplanetary medium and the nonthermal Saturn radiation. Solar wind and interplanetary magnetic field properties were combined to form quantities known to be important in controlling terrestrial magnetospheric processes.The Voyager 2 data set used in this investigation consists of 237 days of Saturn preencounter measurements. However, due to the immersion of Saturn and the Voyager 2 spacecraft into the extended Jupiter magnetic tail, substantial periods of the time series were lacking solar wind data. To cope with this problem a superposed epoch method (CHREE analysis) was used. The results indicate the superiority of the quantities containing the solar wind density in stimulating the radio emission of Saturn — a result found earlier using Voyager 1 data — and the minor importance of quantities incorporating the interplanetary magnetic field.     

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.     

Correlations between neutral and ionized solar wind

We report results of a statistical study correlating ionized solar wind (ISW) fluxes observed by ACE during late 2000 and throughout 2001 with neutral solar wind (NSW) fluxes observed by IMAGE/LENA over the same period. The average correlation coefficient between the neutral and ionized solar wind is 0.66 with correlations greater than 0.80 occurring about 29% of the time. Correlations appear to be driven by high solar wind flux variability, similar to results obtained by in situ multi-spacecraft correlation studies. In this study, however, IMAGE remains inside the magnetosphere on over 95% of its orbits. As a function of day of year, or equivalently ecliptic longitude, the slope of the relationship between the neutral solar wind flux and the ionized solar wind flux shows an enhancement near the upstream direction, but the symmetry point appears shifted toward higher ecliptic longitudes than the interstellar neutral (ISN) flow direction by about 20°. The estimated peak interstellar neutral upstream density inside of 1 AU is about 7 × 10⁻³ cm⁻³.     

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.