On the off-ecliptic distributions of anomalous cosmic rays and their relation to the large-scale geometry of the heliospheric shock

The scenario explaining the origin of the anomalous component of cosmic rays (ACR) implies a close relation between these high energy particles and the solar wind termination shock representing their main acceleration region. Consequently, one should expect the ACR distributions in the heliosphere to reflect some information about the structure as well as the large-scale geometry of the shock. We study the influence of a non-spherically symmetric heliospheric shock on the off-ecliptic — i.e. high latitude — ACR distributions using a two-dimensional model including their anisotropic diffusion and drift in the heliospheric magnetic field as well as a solar wind flow dependent on the heliographic latitude. The model calculations are used to investigate the probability of a possible polar elongation of the heliospheric shock from observations of the distributions of the ACR at high latitudes during solar minimum conditions.     

Implications of a weak termination shock

Recent observations from the Voyager spacecraft have suggested that the spectrum of the anomalous cosmic ray component is relatively steep at the termination shock, which is believed to be responsible for accelerating these particles. This conclusion argues that the termination shock must be weak, which in turn requires that the upstream Mach number in the solar wind must be quite low, 2.4. It is pointed out that such conditions are unlikely to prevail at all locations along the shock front. However, it is possible for such conditions to exist at the interface between high speed streams at high heliographic latitudes and the region at low latitudes where high and low speed streams have interacted and come into equilibrium. This discussion suggests a preferred location for the injection of the anomalous component into the shock acceleration process.     

The Theory of Anomalous Cosmic Rays

Anomalous cosmic rays are a heliospheric phenomenon in which interstellar neutral atoms stream into the heliosphere, are ionized by either solar radiation or the solar wind, and are subsequently accelerated to very high energies, greater than 1 GeV. Current thinking has the bulk of the acceleration to very-high energies taking place, by the mechanism of diffusive shock acceleration, at the termination shock of the solar wind. Detailed two-dimensional numerical simulations and models based on this picture show broad agreement with a number of the observed properties of anomalous cosmic rays. Recent improvements to this picture include the observation of multiply charged cosmic rays and the suggestion that some "preacceleration" of the initially ionized particles occurs in the inner heliosphere.     

Stream Interaction Regions at High Heliographic Latitudes During Ulysses12/22/2004 6:25PM Second Polar Orbit

Ulysses observed well-defined stream interaction regions, SIRs, associated with solar wind stream structure up to a latitude of S65° and shocks to at least a latitude of S71° during the second polar orbit. These SIRs and shocks produced a substantial heliospheric processing of the solar wind. Only a subset of the SIRs recurred on successive solar rotations and only about half of the well-defined SIRs observed poleward of S9.8° were bounded by forward-reverse shock pairs. The majority of the SIRs had local magnetic topologies and azimuthal orientations similar to, but meridional tilts different from, those observed in the first polar orbit when most SIRs corotated with the Sun. The irregular meridional tilts presumably were a consequence of a complex coronal geometry and the temporally evolving nature of the solar wind flow at this time. A lack of reverse shocks poleward of S54° (with one exception) and a lack of well defined SIRs poleward of S65° is evidence that SIRs develop more slowly with distance at high latitudes. This revised version was published online in August 2006 with corrections to the Cover Date.     

The high latitude heliospheric magnetic field

As the Ulysses spacecraft approaches its first pass under the south pole of the sun, it is an appropriate time to review our current knowledge and predictions regarding the three dimensional behaviour of the heliospheric magnetic field, in particular at high heliographic latitudes. Optical techniques for measuring the photospheric magnetic field and observations of coronal brightness structures provide indications of the behaviour of the source of the heliospheric field in the corona. As the coronal fields are carried out into the heliosphere by the solar wind, from Parker's model we would expect that the spiral field observed in the equatorial plane should gradually unwind with latitude leading to open, approximately radial, field lines over the polar regions. Predictions of departures from, and models extending this simple picture are discussed. Both the Pioneer and Voyager spacecraft have spent brief periods in the regions above the maximum latitude of the heliospheric current sheet-relevant results from these missions are reviewed as well as results from the early stages of the out-of-ecliptic phase of the Ulysses mission. The configuration of the coronal magnetic field exhibits a strong dependence on the phase of the solar activity cycle. While the forthcoming Ulysses polar passes take place near to solar minimum, the different conditions which might be encountered on a second orbit of the sun at solar maximum are described.     

Interstellar and Inner Source Pickup Ions Observed with Swics on Ulysses

Many species of pickup ions, both of interstellar origin and from an inner, distributed source have been discovered using data from the Solar Wind Ion Composition Spectrometer (SWICS) on Ulysses. Velocity distribution functions of these ions were measured for the first time over heliocentric distances between 1.35 and 5.4 AU, both at high and low latitudes, and in the disturbed slow solar wind as well as the steady fast wind of the polar coronal holes. This has given us the first glance at plasma properties of suprathermal ions in various solar wind flows, and is enabling us to study the chemical and, in the case of He, the isotopic composition of the local interstellar cloud. Among the new findings are (a) the surprisingly weak pitch-angle scattering of low rigidity, suprathermal ions leading to strongly anisotropic velocity distributions in radial magnetic fields, (b) the efficient injection and consequent acceleration of pickup ions, especially He+ and H+, in the turbulent solar wind, and (c) the discovery of a new extended source releasing carbon, oxygen, nitrogen and possibly other atoms and molecules in the inner solar system. Pickup ion measurements are now used to study the characteristics of the local interstellar cloud (LIC) and, in particular, to determine accurately the abundance of atomic H, He, N, O, and Ne, the isotopes of He and Ne, as well as the ionization fractions of H and He in the LIC. Pickup ion observations allow us to infer the location of the termination shock and, in combination with measurements of anomalous cosmic rays, to investigate termination shock acceleration mechanisms. This revised version was published online in June 2006 with corrections to the Cover Date.     

Physics of the Solar Wind–Local Interstellar Medium Interaction: Role of Magnetic Fields

The interaction of the solar wind with the local interstellar medium is characterized by the self-consistent coupling of solar wind plasma, both upstream and downstream of the heliospheric termination shock, the interstellar plasma, and the neutral atom component of interstellar and solar wind origin. The complex coupling results in the creation of new plasma components (pickup ions), turbulence, and anomalous cosmic rays, and new populations of neutral atoms and their coupling can lead to energetic neutral atoms that can be detected at 1 AU. In this review, we discuss the interaction and coupling of global sized structures (the heliospheric boundary regions) and kinetic physics (the distributions that are responsible for the creation of energetic neutral atoms) based on models that have been developed by the University of Alabama in Huntsville group.     

Voyager observations of the magnetic field, interstellar pickup ions and solar wind in the distant heliosphere

Voyagers 1 and 2 are now observing the latitudinal structure of the heliospheric magnetic field in the distant heliosphere (the legion between - 30 AU and the termination shock). Voyager 2 is observing the influence of the interstellar medium on the solar wind. The pressure of the interstellar pickup protons, measured by their contribution to pressure balanced structures, is greater than or equal to the magnetic pressure and much greater than the thermal pressures of the solar wind protons and electrons in the distant heliosphere. The solar wind speed is observed to decrease and the proton temperature increase with increasing distance from the sun. This may result from the production of pickup ions by the charge exchange process with the interstellar neutrals. The introduction of the pickup ions into the dynamics of the magnetized solar wind plasma appears to be an important new process which must be considered in future theoretical studies of the termination shock and boundary with the local interstellar medium.     

The heliospheric magnetic field

The heliospheric magnetic field (HMF) is an important component of the heliospheric medium. It has been the subject of extensive studies for the past thirty five years. There is a very large observational data base, mostly from the vantage point of the ecliptic plane, but now also from the solar polar regions, from the Ulysses mission. This review aims to present its most important large scale characteristics. A key to understand the HMF is to understand the source functions of the solar wind and magnetic fields close to the sun. The development of new modelling techniques for determining the extent and geometry of the open magnetic field regions in the corona, the sources of the solar wind and the HMF has provided a new insight into the variability of the source functions. These are now reasonably well understood for the state of the corona near solar minimum. The HMF at low-to-medium heliolatitudes is dominated, near solar minimum, by the Corotating Interaction Regions (CIRs) which arise from the interaction of alternating slow and fast solar wind streams, and which, in turn, interact in the outer heliosphere to form the large scale Merged Interaction Regions. The radial component of the HMF is independent of heliolatitude; the average direction is well organised by the Parker geometry, but with wide distributions around the mean, due, at high latitudes, to the presence of large amplitude, Alfvénic fluctuations. The HMF at solar maximum is less well understood, due in part to the complexity of the solar source functions, and partly to the lack of three dimensional observations which Ulysses is planned to remedy at the next solar maximum. It is suggested that the in-ecliptic conditions in the HMF, largely determined by the dynamics of transients (Coronal Mass Ejections) may also be found at high latitudes, due to the wide latitude distribution of the CMEs.     

Heliospheric ENA fluxes: How Sensitive Are They to the Ionization State of LIC?

We explore the sensitivity of the fluxes of heliospheric energetic neutral atoms (ENA) at 1 AU to the ionization state of the local interstellar cloud (LIC). The solar wind plasma is compressed and heated in the termination shock transition. The shocked solar plasma is convected toward the heliospheric tail in the heliosheath, the region between the termination shock and the heliopause. The ENAs are produced in charge exchange of the plasma protons and background neutral gas and can be readily detected at 1 AU. The expected ENA fluxes depend on the shocked plasma density, temperature, and velocity in the heliosheath. The size and structure of the heliospheric interface region depend on the parameters of the interstellar medium. ENA fluxes would thus reveal the LIC parameters. We demonstrate the sensitivity of the heliospheric ENA fluxes to the ionization state of the LIC. The axi-symmetric model of the solar wind/LIC interaction includes the self-consistent treatment of the plasma-gas coupling and Monte Carlo simulations of the neutral gas distribution. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Heliospheric Magnetic Field and Its Extension to the Inner Heliosheath

The general structure of the heliospheric magnetic field is well known and has been extensively studied, mostly in the inner heliosphere, out to the orbit of Saturn. Beyond 10 AU, the Pioneer and now the Voyager spacecraft have provided a view of the outer heliosphere. Its structure is strongly affected by large-scale phenomena originating in the Sun’s activity, such as the pattern of fast and slow solar wind streams around solar minimum that lead to Corotating Interaction Regions, and the increased frequency and strength of Coronal Mass Ejections around solar maximum. The large current sheet that separates the dominant magnetic polarities in the heliospheric medium, the Heliospheric Current Sheet, provides a variable structure that evolves from a relatively simple geometry close to the solar equatorial plane to what is likely to be a highly complex and dynamic surface reaching to high heliolatitudes at high levels of solar activity. The magnetic field observed in a fluctuating, dynamical heliosheath differs considerably from that in a static heliosheath. In particular, the time between current-sheet crossings (sectors) is quite sensitive to the radial speed of the solar-wind termination shock. If an inwardly moving termination shock moves past an observer on a slowly moving spacecraft, the time between current-sheet crossings in the heliosheath becomes larger, and can become very large, for reasonably expected inward shock speeds. This effect may help to explain recent observations of the magnetic field from the Voyager 1 spacecraft, where, in the heliosheath, the magnetic field remained directed outward from the Sun for several months without a current-sheet crossing. The crossings finally resumed and now occur somewhat regularly. In addition, the magnetic fluctuations in the heliosheath are observed to be quite different from those in the supersonic upstream solar wind.     

The acceleration of pickup ions

The well-established association of pickup ions with anomalous cosmic rays shows that acceleration of pickup ions to energies above 1 GeV occurs. At present, diffusive shock acceleration of the pickup ions at the termination shock of the solar wind seems to be the best candidate for acceleration to the high energies of anomalous cosmic rays, accounting well for many of their observed properties. However, it is shown that acceleration of pickup ions from their initial energies by this process appears to be difficult at very strong, nearly perpendicular shocks such as the termination shock. This injection problem remains without a clear solution. A number of alternatives have been proposed for the initial acceleration of pickup ions to the point where diffusive acceleration at the termination shock can take over, but none of these processes has yet emerged as a clear favorite.     

Composition of Anomalous Cosmic Rays

The “classic” anomalous cosmic ray (ACR) component originates as interstellar neutral atoms that drift into the heliosphere, become ionized and picked up by the solar wind, and carried to the outer heliosphere where the pickup ions are accelerated to hundreds of MeV, presumably at the solar wind termination shock. These interstellar ACRs are predominantly singly charged, although higher charge states are present and become dominant above ～350 MeV. Their isotopic composition is like that of the solar system and unlike that of the source of galactic cosmic rays. A comparison of their energy spectra with the estimated flux of pickup ions flowing into the termination shock reveals a mass-dependent acceleration efficiency that favors heavier ions. There is also a heliospheric ACR component as evidenced by “minor” ACR ions, such as Na, Mg, S, and Si that appear to be singly-ionized ions from a source likely in the outer heliosphere.     

Particle acceleration at the Sun and in the heliosphere

Energetic particles are accelerated in rich profusion at sites throughout the heliosphere. They come from solar flares in the low corona, from shock waves driven outward by coronal mass ejections (CMEs), from planetary magnetospheres and bow shocks. They come from corotating interaction regions (CIRs) produced by high-speed streams in the solar wind, and from the heliospheric termination shock at the outer edge of the heliospheric cavity. We sample many populations near Earth, but can distinguish them readily by their element and isotope abundances, ionization states, energy spectra, angular distributions and time behavior. Remote spacecraft have probed the spatial distributions of the particles and examined new sources in situ. Most acceleration sources can be ‘seen’ only by direct observation of the particles; few photons are produced at these sites. Wave-particle interactions are an essential feature in acceleration sources and, for shock acceleration, new evidence of energetic-proton-generated waves has come from abundance variations and from local cross-field scattering. Element abundances often tell us the physics of the source plasma itself, prior to acceleration. By comparing different populations, we learn more about the sources, and about the physics of acceleration and transport, than we can possibly learn from one source alone. This revised version was published online in August 2006 with corrections to the Cover Date.     

Composition of anomalous cosmic rays and implications for the heliosphere

We use energy spectra of anomalous cosmic rays (ACRs) measured with the Cosmic Ray instrument on the Voyager 1 and 2 spacecraft during the period 1994/157-313 to determine several parameters of interest to heliospheric studies. We estimate that the strength of the solar wind termination shock is 2.42 (–0.08, +0.04). We determine the composition of ACRs by estimating their differential energy spectra at the shock and find the following abundance ratios: H/He = 5.6 (–0.5, ⁺0.6), C/He = 0.00048 ± 0.00011, N/He = 0.011 ± 0.001, O/He = 0.075 ± 0.006, and Ne/He = 0.0050 ± 0.0004. We correlate our observations with those of pickup ions to deduce that the long-term ionization rate of neutral nitrogen at 1 AU is 8.3 × 10^–7 s^–1 and that the charge-exchange cross section for neutral N and solar wind protons is 1.0 × 10^–15 cm² at 1.1 keV. We estimate that the neutral C/He ratio in the outer heliosphere is 1.8(–0.7, +0.9) × 10^–5. We also find that heavy ions are preferentially injected into the acceleration process at the termination shock.     

Latitudinal Variation of the Underlying Heliospheric Magnetic Field Direction: Comparison of the Ulysses First and Second Orbits

We discuss the underlying direction of the heliospheric magnetic field measured by Ulysses in the latitude range 6° S-65° S by examining distributions of the magnetic field azimuthal angle with respect to the simple Parker spiral model. During the first Ulysses traversal of this latitude range in 1992–1994, while solar activity was declining, the shape of the distributions obtained at high latitudes in the fast solar wind differed from that at lower latitudes. In the present data set, obtained during rising solar activity, both field polarities are present at all latitudes and the peaks of the distributions agree with the predicted spiral direction to first approximation. However, compared to the first orbit, a significantly greater percentage of the observed field vectors have large deviations from the spiral direction. This revised version was published online in August 2006 with corrections to the Cover Date.     

Energetic Ions Observed at Low to High Latitudes in the Southern Heliosphere During Declining and Rising Solar Activity

We present and compare observations of energetic protons during the two first transits of the Ulysses spacecraft from low to high latitudes in the southern heliosphere. Protons in the energy range 1.8–3.8 MeV from the COSPIN experiment are studied for global trends and in relation to some ambient structures in the solar wind (corotating interaction region, forward/reverse shock). The global trends show the large dependence on the heliospheric condition and solar activity, including indications of a larger ambient particle population during the rising phase of solar activity and more efficient solar wind particle accelerators during the declining phase. More enhancements in the proton flux intensity are time associated with forward shocks than reverse contrary to first pass. Recurrent structures are found even during the second transit. Some latitude dependent periodicities are observed that could relate to the differential solar rotation. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Heliospheric Termination Shock

The heliospheric termination shock is a vast, spheroidal shock wave marking the transition from the supersonic solar wind to the slower flow in the heliosheath, in response to the pressure of the interstellar medium. It is one of the most-important boundaries in the outer heliosphere. It affects energetic particles strongly and for this reason is a significant factor in the effects of the Sun on Galactic cosmic rays. This paper summarizes the general properties and overall large-scale structure and motions of the termination shock. Observations over the past several years, both in situ and remote, have dramatically revised our understanding of the shock. The consensus now is that the shock is quite blunt, is with the front, blunt side canted at an angle to the flow direction of the local interstellar plasma relative to the Sun, and is dynamical and turbulent. Much of this new understanding has come from remote observations of energetic charged particles interacting with the shock, radio waves and radiation backscattered from interstellar neutral atoms. The observations and the implications are discussed.     

Shock acceleration of energetic particles in the heliosphere

The theory of shock acceleration of energetic particles is briefly discussed and reviewed with an emphasis on clarifying the apparent distinction between the V × B and Fermi mechanisms. Attention is restricted to those situations in which the energetic particles do not themselves influence the given shock structure. In particular, application of the theory to the acceleration of energetic particles in corotating interaction regions (CIR) in the solar wind is presented. Here particles are accelerated at the forward and reverse shocks which bound the CIR by being compressed between the shock fronts and magnetic irregularities upstream from the shocks, or by being compressed between upstream irregularities and those downstream from the shocks. Particles also suffer adiabatic deceleration in the expanding solar wind, an effect not included in previous shock models for acceleration in CIRs. The model is able to account for the observed exponential spectra at Earth, the observed behavior of the spectra with radial distance, the observed radial gradients in the intensity, and the observed differences in the intensity and spectra at the forward and reverse shocks.Calculations and resulting energy spectra are also presented for shock acceleration of energetic particles in large solar flare events. Based on the simplifying assumption that the shock evolves as a spherically symmetric Sedov blast wave, the calculation yields the time-integrated spectrum of particles initially injected at the shock which eventually escape ahead of the shock into interplanetary space. The spectra are similar to those observed at Earth. Finally further applications are suggested.An invited paper presented at STIP Workshop on Shock Waves in the Solar Corona and Interplanetary Space, 15–19 June, 1980, Smolenice, Czechoslovakia.     

Coronal mass ejections at high heliographic latitudes: Ulysses

Nine coronal mass ejections (CMEs) have been detected in the solar wind by the Ulysses plasma experiment between 31° and 61° South. One of these events, which was also a magnetic cloud, was directly associated with an event observed by the soft X-ray telescope on Yohkoh in which large magnetic loops formed in the solar corona directly beneath Ulysses. This association suggests that the flux rope topology of the magnetic cloud resulted from reconnection between the legs of neighboring magnetic loops within the rising CME. The average CME speed (740 km s^–1) at these latitudes was comparable to that of the normal solar wind there and is much greater than average CME speeds observed either in the solar wind in the ecliptic plane or in the corona close to the Sun. We suggest that the same basic acceleration process applies to both slow CMEs and the normal solar wind at any latitude.